Identification of a Novel Eight-lncRNA Prognostic Signature for HBV-HCC and Analysis of Their Functions Based on Coexpression and ceRNA Networks

Zhao, Xiaonan; Bai, Zhenzi; Li, Chenghua; Sheng, Chuanlun; Li, Hongyan

doi:https://doi.org/10.1155/2020/8765461

BioMed Research International

On this page

Abstract Introduction Materials and Methods Results Discussion Conclusion Data Availability Conflicts of Interest Authors’ Contributions References Copyright Related Articles

Research Article | Open Access

Volume 2020 | Article ID 8765461 | https://doi.org/10.1155/2020/8765461

Identification of a Novel Eight-lncRNA Prognostic Signature for HBV-HCC and Analysis of Their Functions Based on Coexpression and ceRNA Networks

Xiaonan Zhao,¹Zhenzi Bai,¹Chenghua Li,¹Chuanlun Sheng,¹and Hongyan Li¹

Academic Editor: Gerald J. Wyckoff

Received01 Nov 2019

Revised09 Feb 2020

Accepted12 Feb 2020

Published15 Apr 2020

Abstract

Studies have demonstrated the prognosis potential of long noncoding RNAs (lncRNAs) for hepatocellular carcinoma (HCC), but specific lncRNAs for hepatitis B virus- (HBV-) related HCC have rarely been reported. This study was aimed at identifying a lncRNA prognostic signature for HBV-HCC and exploring their underlying functions. The sequencing dataset was collected from The Cancer Genome Atlas database as the training set, while the microarray dataset was obtained from the European Bioinformatics Institute database (E-TABM-36) as the validation set. Univariate and multivariate Cox regression analyses identified that eight lncRNAs (TSPEAR-AS1, LINC00511, LINC01136, MKLN1-AS, LINC00506, KRTAP5-AS1, ZNF252P-AS1, and THUMPD3-AS1) were significantly associated with overall survival (OS). These eight lncRNAs were used to construct a risk score model. The Kaplan-Meier survival curve results showed that this risk score can significantly differentiate the OS between the high-risk group and the low-risk group. Receiver operating characteristic curve analysis demonstrated that this risk score exhibited good prediction effectiveness (area under the curve for the training set; for the validation set). Furthermore, this lncRNA risk score was identified as an independent prognostic factor in the multivariate analysis after adjusting other clinical characteristics. The crucial coexpression (LINC00511-CABYR, THUMPD3-AS1-TRIP13, LINC01136-SFN, LINC00506-ANLN, and KRTAP5-AS1/TSPEAR-AS1/MKLN1-AS/ZNF252P-AS1-MC1R) or competing endogenous RNA (THUMPD3-AS1-hsa-miR-450a-TRIP13) interaction axes were identified to reveal the possible functions of lncRNAs. These genes were enriched into cell cycle-related biological processes or pathways. In conclusion, our study identified a novel eight-lncRNA prognosis signature for HBV-HCC patients and these lncRNAs may be potential therapeutic targets.

1. Introduction

Liver cancer, 90% of which is hepatocellular carcinoma (HCC), is one of the most common malignancies and the leading cause of cancer-related deaths worldwide [1, 2]. Despite great advances made in therapeutic approaches, relapse, progression, and metastasis rates remain high, which leads to poor prognostic outcomes [3]. Hepatitis B virus (HBV) infection is recognized as the most important risk factor associated with the biological aggressiveness of HCC and patients’ survival [4]. Therefore, it may be an important issue to screen prognostic biomarkers and therapeutic targets for patients with HBV-HCC in order to optimize treatment, prevent progression, and improve the prognosis.

Recently, more and more evidence has highlighted the important roles of long noncoding RNAs (lncRNAs) in HBV-HCC via directly influencing the expression of their neighboring genes (coexpression hypothesis) or acting as microRNA (miRNAs) sponges to negatively regulate the expression of miRNA target genes (competing endogenous RNA (ceRNA) hypothesis) [5]. For example, Chen et al. found that LINC01152 was significantly upregulated in HBV-positive HCC tissues and cells treated by HBV X protein (HBx). Overexpression of LINC01152 could increase HCC cell proliferation and promote tumor formation in nude mice by directly binding to the promoter region of interleukin- (IL-) 23 to increase its transcription [6]. Hu et al. identified that lncRNA WEE2-AS1 was overexpressed in HCC tissues and was positively correlated to HBV infection. Compared with the low expression and HBV-negative groups, HBV-HCC patients with high expression of WEE2-AS1 exhibited the worst prognosis. WEE2-AS1 accelerated the proliferation, migration, invasion, and cell cycle progression of HCC cells by upregulating the downstream Fermitin family member 3 (FERMT3) [7]. Lin et al. observed that the expression of miR143HG was markedly decreased in HCC tissues and its expression was associated with the presence of HBV infection. Highly expressed miR143HG can serve as an independent prognostic factor to predict a good prognosis. In vitro analysis revealed that miR143HG may exert its tumor suppressor roles by sponging miR-155 and then promoting the expression of miR-155 target gene adenomatous polyposis coli (APC) which can inhibit the Wnt/β-catenin signaling pathway [8]. Fan et al. demonstrated that the expression of lncRNA n335586 was significantly increased in HBV-positive HCC tissues and cells. lncRNA n335586 may contribute to the migration, invasion, and metastasis of HCC cells by facilitating the expression of its host gene creatine kinase, mitochondrial 1A (CKMT1A), through competitively binding to miR-924 [9]. These findings suggest that lncRNAs may be potential therapeutic targets and prognostic biomarkers for HBV-HCC patients. However, the lncRNAs specifically associated with HBV-HCC remain rarely reported and most of these studies focused on the prognostic roles of single lncRNA, and the attempt to identify a lncRNA prognostic signature for HBV-HCC was limited [10].

The present study is aimed at further identifying a lncRNA prognostic signature for HBV-HCC patients through using the sequencing data in The Cancer Genome Atlas (TCGA) database and a microarray dataset deposited in the European Bioinformatics Institute (EMBL-EBI) which had the larger sample size compared with the dataset used by Liu et al. [10] (44 vs. 19). Moreover, the potential functions of lncRNAs were predicted by constructing a ceRNA regulatory network, in addition to a coexpression network which was used in the study of Liu et al. [10]. Hereby, our study may provide a novel lncRNA prognostic signature and identify potential therapeutic targets for HBV-HCC patients.

2. Materials and Methods

2.1. Data Source

The mRNA- and miRNA-sequencing data (level 3, Illumina HiSeq 2000 RNA Sequencing platform) of HCC were obtained from TCGA database portal (https://portal.gdc.cancer.gov/) on March 25, 2019, which included 373 HCC samples and 50 normal liver controls. A total of 95 HCC samples were used in this analysis because they (1) sequenced for both mRNAs and miRNAs, (2) belonged to HBV alone-positive samples, and (3) had complete clinical data (including survival). This dataset was used as the training set.

Furthermore, a microarray dataset was also collected from the EMBL-EBI database (https://www.ebi.ac.uk/arrayexpress/) under accession number E-TABM-36 (platform: GPL 96[HG-U133A] Affymetrix GeneChip Human Genome HG-U133A) [11] to serve as the validation set. This dataset included 65 samples, 44 of which were HBV-positive and had survival information.

2.2. Differential Expression Analysis

The lncRNAs and mRNAs in RNA-sequencing data were annotated based on the HUGO Gene Nomenclature Committee (HGNC; http://www.genenames.org/) [12]. The Linear Models for Microarray Data (LIMMA) method (version 3.34.7; https://bioconductor.org/packages/release/bioc/html/limma.html) [13] in the Bioconductor R package (version 3.4.1; http://www.R-project.org/) was applied for the identification of differentially expressed genes (DEGs), lncRNAs (DELs), and miRNAs (DEMs) which were defined as the false discovery rate and . Bidirectional hierarchical clustering was performed using the pheatmap R package (version: 1.0.8; https://cran.r-project.org/web/packages/pheatmap) based on centered Pearson’s correlation to generate heat maps of differentially expressed RNAs.

2.3. Identification of a Prognostic lncRNA Signature Using the Training Set

The survival package (version 2.41-1; http://bioconductor.org/packages/survivalr/) was used for all statistical analyses. First, univariate Cox regression analysis was performed to mine DEGs, DELs, and DEMs that were associated with overall survival (OS) in the training set (log-rank value < 0.05). Second, multivariate Cox regression analysis was conducted to further select independent prognostic lncRNAs, which were used to construct a lncRNA prognostic signature model according to the following formula: where Exp_lncRNA was the expression level of DELs and was the regression coefficient for the DELs in multivariate Cox hazard model analysis.

Third, using the median risk score as the threshold, the HBV-HCC patients were divided into the low-risk and high-risk groups. The Kaplan-Meier (K-M) survival curves were used to calculate the differences in OS between the high- and low-risk groups. The receiver operating characteristic (ROC) curve and area under the curve (AUC) were also applied to assess the predictive performance of the prognostic risk score.

2.4. Validation of the Prognostic Power of the Risk Score Using the Validation Dataset

The risk score identified in the training dataset was further calculated for the validation dataset. Patients were also divided into the high-risk and low-risk groups according to the risk score. The K-M survival and ROC curves were analyzed.

2.5. The Prognostic Independence of the Risk Score

Univariate and multivariate Cox regression analyses were conducted to assess whether the prognostic performance of signature risk score was independent of other clinical characteristics (including age, gender, neoplasm histologic grade, pathologic TNM stage, vascular invasion, and recurrence) for patients with HBV-HCC.

2.6. Construction of lncRNA Coexpression and ceRNA Networks

The DIANA-LncBase (version 2.0; http://carolina.imis.athena-innovation.gr/diana_tools/web/index.php?r=lncbasev2/index-predicted) [14] database was used to predict the interactions between signature DELs and prognostic DEMs. The starBase database (version 2.0; http://starbase.sysu.edu.cn/) [15] was used to predict the interactions between prognostic DEMs and prognostic DEGs. The DEL-DEM and DEM-DEG interactors which had opposite expression trend were integrated to construct the ceRNA network, which was visualized using the Cytoscape software (version 3.6.1; http://www.cytoscape.org/) [16].

The coexpression network was constructed based on the correlation between signature DELs and DEGs which was calculated using the tcor.test function (https://stat.ethz.ch/R-manual/R-devel/library/stats/html/cor.test.html) in R to generate the Pearson correlation coefficients (PCC). Only the coexpression pairs with were selected to draw the network using Cytoscape (version 3.4; http://www.cytoscape.org/) [16].

2.7. Function and Pathway Enrichment

Gene Ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed using clusterProfiler (version 3.6.0; http://bioconductor.org/packages/release/bioc/html/clusterProfiler.html) to reveal the functions of DEGs in the lncRNA-related coexpression and ceRNA networks [17]. was set as the cut-off value.

3. Results

3.1. Differential Expression Analysis

After removing the expression value of zero, a total of 13,454 mRNAs, 1,238 lncRNAs, and 1,037 miRNAs were annotated. Using the LIMMA method, 1,898 differentially expressed RNAs were identified between the 95 HBV-HCC tissues and 50 controls, including 1,214 DEGs (318 downregulated; 896 upregulated), 584 DELs (37 downregulated; 547 upregulated), and 100 DEMs (42 downregulated; 58 upregulated) (Figure 1(a)). The heat map showed that the HBV-HCC and normal tissues were well distinguished by these differentially expressed RNAs (Figure 1(b)).

(a)

(b)

3.2. Identification of a lncRNA Signature with Prognostic Value

A total of 260 differentially expressed RNAs were identified to be significantly associated with OS by univariate Cox regression analysis, including 189 DEGs, 39 DELs, and 32 DEMs. To determine the optimal prognostic lncRNAs, multivariate Cox proportional hazards regression was further adopted to evaluate the prognostic independence of the above-identified 39 prognostic lncRNAs. As a result, 8 lncRNAs (TSPEAR-AS1, LINC00511, LINC01136, MKLN1-AS, LINC00506, KRTAP5-AS1, ZNF252P-AS1, and THUMPD3-AS1) were found to be still significant in the multivariate analysis (Table 1), indicating that these lncRNAs were independent prognostic biomarkers. As shown in Table 1, LINC00511, LINC01136, LINC00506, and THUMPD3-AS1 may be prognostic risk factors for OS because of negative coefficient and (that is, patients with high expression of lncRNAs had shorter survival) while TSPEAR-AS1, MKLN1-AS, KRTAP5-AS1, and ZNF252P-AS1 may be prognostic protective factors due to negative coefficient and (that is, patients with high expression of lncRNAs had longer survival).

The risk score model of this eight-lncRNA signature was developed for each patient according to the following formula: . The patients were classified into the high-risk and low-risk groups based on the median risk score. The K-M survival curve results showed that this eight-lncRNA signature can significantly differentiate the OS between the predicted two risk groups, with worse OS in HBV-HCC patients with high-risk scores (training dataset: , , Figure 2(a); validation dataset: , , Figure 2(c)). ROC analysis was subsequently performed to evaluate the prediction accuracy of this risk score of the eight-lncRNA signature. The results showed that the AUC of the ROC curve was, respectively, 0.990 and 0.903 for the training (Figure 2(b)) and validation datasets (Figure 2(d)), suggesting a good prediction performance.

(a)

(b)

(c)

(d)

Figure 2

Prognosis prediction of the risk score model for overall survival in patients with HBV-HCC. (a) Kaplan-Meier plots of overall survival for patients with a high- or low-risk score of the training dataset. (b) ROC curves of the risk score based on 8-lncRNA signature using the training dataset. (c) Kaplan-Meier plots of overall survival for patients with a high- or low-risk score of the validation dataset. (d) ROC curves of the risk score based on an 8-lncRNA signature using the validation dataset. ROC: receiver operating characteristic; AUC: area under the ROC curve; HR: hazard ratio; TCGA: The Cancer Genome Atlas.

3.3. Independence of the Eight-lncRNA Signature for Survival Prediction

Univariate and multivariate Cox regression analyses were also performed to further assess whether the prognostic value of the eight-lncRNA signature was independent of other clinical variables. As expected, the risk score status was identified as a significant prognostic factor in both of the univariate and multivariate Cox regression analyses (Table 2), indicating the independence of the eight-lncRNA signature for prognostic prediction.

3.4. Functional Analysis of lncRNAs

A total of 7 downregulated prognosis-related DEMs were predicted to interact with 4 upregulated signature DELs (LINC00511, LINC00506, MKLN1-AS, and THUMPD3-AS1) in the DIANA-LncBase database, which constituted 9 interaction pairs, while 30 upregulated prognostic DEGs were predicted to interact with 5 downregulated prognostic DEMs (hsa-miR-154, hsa-miR-187, hsa-miR-450a, hsa-miR-299, and hsa-miR-133b) by the starBase database, which constitute5d 39 interaction pairs. These interaction pairs were integrated to construct a lncRNA ceRNA network (Figure 3), in which THUMPD3-AS1-related ceRNA axes, such as THUMPD3-AS1-hsa-miR-450a-TRIP13 (thyroid hormone receptor interactor 13), may be potentially credible because the expression and prognosis trend was similar between this lncRNA and mRNA, but opposite between this lncRNA and miRNA (Table 1).

Furthermore, 7 upregulated signature DELs (LINC00511, LINC01136, MKLN1-AS, LINC00506, KRTAP5-AS1, ZNF252P-AS1, and THUMPD3-AS1) were predicted to coexpress with 75 prognostic DEGs according to the threshold of , which constituted 117 interaction pairs, and were used to construct the coexpression network (Figure 4). According to the principle of consistent expression and prognosis trend, the following coexpression pairs were suggested to be important, including LINC00511-CABYR (calcium binding tyrosine phosphorylation regulated) (), LINC01136-SFN (stratifin) (), THUMPD3-AS1-TRIP13 (), and KRTAP5-AS1-MC1R (melanocortin 1 receptor) () (Table 1; Figure 4). Furthermore, LINC00506-ANLN (anillin actin binding protein) (), TSPEAR-AS1 ()/MKLN1-AS ()/ZNF252P-AS1 ()-MC1R (melanocortin 1 receptor) coexpression relationships also complied with the principle, although their PCC were less than 0.4.

Function enrichment analysis screened 35 GO terms and 4 KEGG pathways (Table 3; Figure 5) for the DEGs of the ceRNA and coexpression networks. All of the DEGs included in the above crucial ceRNA axes or coexpression pairs were enriched, including the GO:0000279~M phase (ANLN, MC1R, and TRIP13), GO:0022402~cell cycle process (ANLN, MC1R, and TRIP13), GO:0006974~response to DNA damage stimulus (SFN, TRIP13, and ANLN), GO:0005856~cytoskeleton (CABYR), hsa04110:cell cycle (SFN), and hsa04115:p53 signaling pathway (SFN).

(a)

(b)

4. Discussion

In the present study, an eight-lncRNA signature (TSPEAR-AS1, LINC00511, LINC01136, MKLN1-AS, LINC00506, KRTAP5-AS1, ZNF252P-AS1, and THUMPD3-AS1) was identified as an independent prognostic biomarker for HBV-positive HCC patients. ROC curve analysis showed that the risk score established by this eight-lncRNA signature had high prediction accuracy for OS, with AUC of 0.990 for the training set and 0.903 for the validation set, respectively. Further mechanism analyses uncovered upregulated THUMPD3-AS1 may function as a ceRNA to sponge hsa-miR-450a and then regulate TRIP13. In addition, THUMPD3-AS1 may also directly regulate TRIP13 by coexpression with it. The other lncRNAs may exert protective or risk roles for the development of HBV-positive HCC mainly by coexpression with their downstream genes, such as LINC00511-CABYR, LINC01136-SFN, LINC00506-ANLN, and KRTAP5-AS1/TSPEAR-AS1/MKLN1-AS/ZNF252P-AS1-MC1R. These genes were involved in cell cycle-related GO terms or KEGG pathways.

Although TCGA data with essentially the same samples were also used as the training set to screen the lncRNA signature in the study of Liu et al. [10], there were no common lncRNAs between our identified lncRNA signature and that of Liu et al. [10] (DGCR9, GBA3, HCG4, NAT8B, NBR2, PART1, RFPL1S, SLC22A18AS, and TCL6). This may be attributed to the analysis method difference (i.e., only lncRNAs in partial module were used for univariate Cox regression in Liu et al. [10], but all DELs in our study). Furthermore, only TSPEAR-AS1 (upregulated, protective) [18] and LINC00511 (upregulated, risk) [19], with the prognosis trend in line with ours, were reported as prognostic factors for the whole HCC, while KRTAP5-AS1 (upregulated, protective; similar to ours) [20] and THUMPD3-AS1 (upregulated, protective; contrast to ours) [21] were identified to be associated with OS in papillary thyroid carcinoma and lung adenocarcinoma, respectively. No studies investigated the prognostic values of LINC01136, MKLN1-AS, LINC00506, and ZNF252P-AS1. These findings indicated that our lncRNA signature may be a first-time identified biomarker of the prognosis for HBV-HCC.

Besides, the predictive accuracy of our eight-lncRNA signature was comparable with that of 9-lncRNA signature for HBV-HCC identified by Liu et al. (0.990 vs. 0.953) [10] and obviously higher than that of the published studies focusing on the whole HCC, such as Shi et al. (8-lncRNA signature, , 0.708, and 0.75 for training, testing, and entire set) [22], Yan et al. (7-lncRNA signature, ) [23], Wang et al. (4-lncRNA signature, ) [24], Ma et al. (4-lncRNA signature, ) [25], Sui et al. (4-lncRNA signature, ) [26], and Zhang et al. (10-lncRNA signature, ) [27]. These results suggested that our eight-lncRNA signature may be a more effective, robust biomarker for HCC prognosis, particularly HBV-positive ones.

Several in vitro studies had demonstrated the functions of LINC00511 in cancer, but most of them focused on the ceRNA mechanism. The study of Zhang et al. suggested that highly expressed LINC0051 directly interacted with miR-29c to suppress its expression and then promoted the upregulation of CDK6, a direct target of miR-29c, leading to an increase in the cell viability of breast cancer cells [28]. Lu et al. concluded that LINC00511 promoted the proliferation, sphere-formation ability, and tumor growth in breast cancer cells by functioning as a ceRNA for miR-185-3p to positively recover E2F1 protein and promote the transcription of Nanog gene [29]. Furthermore, LINC00511-hsa-miR-29b-3p-VEGFA (vascular endothelial growth factor A) [30], LINC00511-miR-765-LAMC2 (laminin subunit gamma 2) [31], and LINC00511-miR-124-3p-CCND2 (cyclin D2) [32] ceRNA axes were also reported for pancreatic ductal adenocarcinoma, tongue squamous cell carcinoma, and glioma, respectively. But the coexpression mechanism of LINC00511 was rarely reported [33]. In this study, we speculated LINC00511 may coexpress with CABYR. Although no experiments were performed to confirm their relationships, previous studies on the roles of CABYR may indirectly explain our conclusion. For example, Li et al. found that the mRNA and protein levels of CABYR-c were significantly higher in HCC tissues than those in the adjacent noncancerous tissues. Silencing of CABYR-c by antisense oligonucleotides significantly inhibited the cell growth of HepG2 cells and arresting the cell cycle [34]. Similarly, knockdown of CABYR-a/b was also reported to increase the cell apoptosis and enhance the chemosensitivity [35]. In line with these studies, our study also identified that LINC00511 and CABYR were upregulated and predicted a poor prognosis.

THUMPD3-AS1 was not included in the ceRNA network of Li et al. [21], and thus, its potential function remained unknown. In this study, we predicted that upregulated THUMPD3-AS1 may directly coexpress with TRIP13 to promote its upregulated expression or regulate it by sponging hsa-miR-450a. Although no research was conducted to verify their interaction relationships, previous studies on the expression and roles of these miRNA and mRNA may indirectly prove our speculation. Ju et al. detected that the expression of TRIP13 was upregulated in HCC tissues and cell lines. Patients with higher expression levels of TRIP13 had significantly shorter survival periods [36]. Downregulation of TRIP13 inhibited HCC cell proliferation, migration, and invasion, promoted apoptosis and cell cycle arrest at S-phase in vitro [36], and suppressed the formation of tumor in vivo [37]. The study of Weng et al. revealed that miR-450a was significantly downregulated in HCC tissues and cells. Ectopic expression of miR-450a in HepG2 cells caused an inhibition of cell proliferation, the mechanism of which was related to downregulation of DNA methyltransferase 3a [38].

LINC01136 was a first-time identified oncogene in our study, and its role was unknown. We predicted it may exert tumor-promoting roles by coexpressing with SFN which had been demonstrated to be oncogenic previously as follows: SFN protein was previously identified to be upregulated in HCC via proteomics analyses [39]. SFN guided cell migration and invasion of breast cancer cells by stabilizing a complex of soluble actin [40]. Suppression of SFN expression by siRNA significantly reduced proliferation activity and the S-phase subpopulation of human lung adenocarcinoma cells and blocked tumor development in mice [41]. By analysis of ten datasets, Hu et al. found that high expression of SFN was significantly associated with worse OS in patients with ovarian cancer [42]. In agreement with these studies, we also found that SFN was upregulated and predicted a poor prognosis in HBV-HCC.

The function of LINC00506 was also not reported previously. We predicted it coexpressed with ANLN to promote carcinogenesis. ANLN was detected to be highly expressed in HCC tissues and cells infected with HBV. Higher ANLN expression was associated with a worse prognosis. Moreover, inhibition of ANLN resulted in growth restraint, colony formation reduction, and apoptosis activation [43] as well as survival time increase [44]. In accordance with these studies, we also found ANLN was highly expressed in HBV-HCC tissues and related to a poor prognosis.

The functions of TSPEAR-AS1 and KRTAP5-AS1 had been predicted based on the ceRNA network previously. For example, TSPEAR-AS1 may sponge miR-424 to regulate the expression of G protein subunit alpha L (GNAL) [18] and lncRNA KRTAP5-AS1 could act as a ceRNA to affect the functions of Claudin-4 [45]. However, their roles remain unclear. In this study, our study suggested that TSPEAR-AS1 and KRTAP5-AS1 as well as MKLN1-AS/ZNF252P-AS1 may exert protective roles for HBV-HCC by coexpressing with MC1R. The studies on MC1R mainly focused on melanoma, and the results indicated that patients had significantly longer OS in melanoma tumors with lower expression of MC1R [46], which was opposite with our expectation. Thus, further studies on MC1R in other cancers especially HBV-HCC are necessary.

There were some limitations in this study. First, the sample size was small. Therefore, large HBV-HCC cohorts should be collected to confirm the prognosis values of our lncRNA signature. Second, this study was only preliminary to predict the functions of lncRNAs in HBV-HCC and further wet experiments should be performed to verify their ceRNA (luciferase assay, overexpression or silencing) or coexpression (coimmunoprecipitation) mechanisms. Third, previous studies have demonstrated that the gene profiles induced by HBV, HCV, and other etiologies were different [47–49]. Thus, theoretically, a distinct signature or results using the same signature may be obtained when a similar analysis was performed [50], which should be confirmed in subsequent studies. Fourth, due to the limited description of treatment regimens in TCGA data, we could not investigate the influence of therapies (hepatectomy, radiofrequency ablation, transarterial embolization, and others) on the prognosis of HBV-HCC patients.

5. Conclusion

In the present study, we constructed an 8-lncRNA risk score model for survival prediction of HBV-HCC based on RNA-Seq and microarray datasets, with the prediction accuracy more than 90%. These signature lncRNAs may be involved in HBV-HCC by coexpressing with mRNAs (LINC00511-CABYR, THUMPD3-AS1-TRIP13, LINC01136-SFN, LINC00506-ANLN, and KRTAP5-AS1/TSPEAR-AS1/MKLN1-AS/ZNF252P-AS1-MC1R) or functioning as a ceRNA (THUMPD3-AS1-hsa-miR-450a-TRIP13) to regulate the cell cycle-related biological processes or pathways.

Data Availability

The origin RNA-seq data used in our study were all downloaded from the TCGA (https://portal.gdc.cancer.gov/) and EMBL-EBI (https://www.ebi.ac.uk/arrayexpress/).

Conflicts of Interest

The authors declare that they have no competing interests.

Authors’ Contributions

XNZ and HYL conceived the design of the original study; XNZ and ZZB performed the acquisition of data; XNZ, CHL, and CLS conducted the statistical analysis. ZZB and HYL were involved with the interpretation of the data; XNZ drafted the manuscript; HYL participated in critical revisions of the manuscript; all authors read and approved the final manuscript.

References

W. Chen, R. Zheng, P. D. Baade et al., “Cancer statistics in China, 2015,” CA: a Cancer Journal for Clinicians, vol. 66, no. 2, pp. 115–132, 2016.
View at: Publisher Site | Google Scholar
R. L. Siegel, K. D. Miller, and A. Jemal, “Cancer statistics, 2018,” CA: a Cancer Journal for Clinicians, vol. 68, no. 1, pp. 7–30, 2018.
View at: Publisher Site | Google Scholar
J. Li, L. Huang, J. Yan, M. Qiu, and Y. Yan, “Liver resection for hepatocellular carcinoma: personal experiences in a series of 1330 consecutive cases in China,” ANZ Journal of Surgery, vol. 88, no. 10, pp. E713–E717, 2018.
View at: Publisher Site | Google Scholar
F. Yang, L. Ma, Y. Yang et al., “Contribution of hepatitis B virus infection to the aggressiveness of primary liver cancer: a clinical epidemiological study in eastern China,” Frontiers in Oncology, vol. 9, 2019.
View at: Publisher Site | Google Scholar
L. Qiu, T. Wang, X. Xu, Y. Wu, Q. Tang, and K. Chen, “Long non-coding RNAs in hepatitis B virus-related hepatocellular carcinoma: regulation, functions, and underlying mechanisms,” International Journal of Molecular Sciences, vol. 18, no. 12, p. 2505, 2017.
View at: Publisher Site | Google Scholar
T. Chen, J. Pei, J. Wang et al., “HBx-related long non-coding RNA 01152 promotes cell proliferation and survival by IL-23 in hepatocellular carcinoma,” Biomedicine & Pharmacotherapy, vol. 115, p. 108877, 2019.
View at: Publisher Site | Google Scholar
Z. Hu, P. Huang, Y. Yan, Z. Zhou, J. Wang, and G. Wu, “Hepatitis B virus X protein related lncRNA WEE2-AS1 promotes hepatocellular carcinoma proliferation and invasion,” Biochemical and Biophysical Research Communications, vol. 508, no. 1, pp. 79–86, 2019.
View at: Publisher Site | Google Scholar
X. Lin, H. Xiaoqin, C. Jiayu, F. Li, L. Yue, and X. Ximing, “Long non-coding RNA miR143HG predicts good prognosis and inhibits tumor multiplication and metastasis by suppressing mitogen-activated protein kinase and Wnt signaling pathways in hepatocellular carcinoma,” Hepatology Research, vol. 49, no. 8, pp. 902–918, 2019.
View at: Publisher Site | Google Scholar
H. Fan, P. Lv, T. Mu et al., “LncRNA n335586/miR-924/CKMT1A axis contributes to cell migration and invasion in hepatocellular carcinoma cells,” Cancer Letters, vol. 429, pp. 89–99, 2018.
View at: Publisher Site | Google Scholar
H. Liu, P. Zhao, X. Jin et al., “A 9-lncRNA risk score system for predicting the prognosis of patients with hepatitis B virus-positive hepatocellular carcinoma,” Molecular Medicine Reports, vol. 20, no. 1, pp. 573–583, 2019.
View at: Publisher Site | Google Scholar
S. M. Kim, S. H. Leem, I. S. Chu et al., “Sixty-five gene-based risk score classifier predicts overall survival in hepatocellular carcinoma,” Hepatology, vol. 55, no. 5, pp. 1443–1452, 2012.
View at: Publisher Site | Google Scholar
S. Povey, R. Lovering, E. Bruford, M. Wright, M. Lush, and H. Wain, “The HUGO Gene Nomenclature Committee (HGNC),” Human Genetics, vol. 109, no. 6, pp. 678–680, 2001.
View at: Publisher Site | Google Scholar
M. E. Ritchie, B. Phipson, D. Wu et al., “limma powers differential expression analyses for RNA-sequencing and microarray studies,” Nucleic Acids Research, vol. 43, no. 7, p. e47, 2015.
View at: Publisher Site | Google Scholar
M. D. Paraskevopoulou, G. Georgakilas, N. Kostoulas et al., “DIANA-LncBase: experimentally verified and computationally predicted microRNA targets on long non-coding RNAs,” Nucleic Acids Research, vol. 41, pp. D239–D245, 2013.
View at: Publisher Site | Google Scholar
J. H. Yang, J. H. Li, P. Shao, H. Zhou, Y. Q. Chen, L. H. Qu et al., Nucleic Acids Research, vol. 39, pp. D202–D209, 2011.
View at: Publisher Site
M. Kohl, S. Wiese, and B. Warscheid, “Cytoscape: software for visualization and analysis of biological networks,” Methods in Molecular Biology, vol. 696, pp. 291–303, 2011.
View at: Publisher Site | Google Scholar
M. V. Kuleshov, M. R. Jones, A. D. Rouillard et al., “Enrichr: a comprehensive gene set enrichment analysis web server 2016 update,” Nucleic Acids Research, vol. 44, no. W1, pp. W90–W97, 2016.
View at: Publisher Site | Google Scholar
J. Ye, J. Zhang, Y. Lv et al., “Integrated analysis of a competing endogenous RNA network reveals key long noncoding RNAs as potential prognostic biomarkers for hepatocellular carcinoma,” Journal of Cellular Biochemistry, vol. 120, no. 8, pp. 13810–13825, 2019.
View at: Publisher Site | Google Scholar
R. P. Wang, J. Jiang, T. Jiang, Y. Wang, and L. X. Chen, “Increased long noncoding RNA LINC00511 is correlated with poor prognosis and contributes to cell proliferation and metastasis by modulating miR-424 in hepatocellular carcinoma,” European Review for Medical and Pharmacological Sciences, vol. 23, no. 8, pp. 3291–3301, 2019.
View at: Publisher Site | Google Scholar
X. You, S. Yang, J. Sui et al., “Molecular characterization of papillary thyroid carcinoma: a potential three-lncRNA prognostic signature,” Cancer Management and Research, vol. 10, pp. 4297–4310, 2018.
View at: Publisher Site | Google Scholar
D. S. Li, J. L. Ainiwaer, I. Sheyhiding, Z. Zhang, and L. W. Zhang, “Identification of key long non-coding RNAs as competing endogenous RNAs for miRNA-mRNA in lung adenocarcinoma,” European Review for Medical and Pharmacological Sciences, vol. 20, no. 11, pp. 2285–2295, 2016.
View at: Google Scholar
Y. M. Shi, Y. Y. Li, J. Y. Lin, L. Zheng, Y. M. Zhu, and J. Huang, “The discovery of a novel eight-mRNA-lncRNA signature predicting survival of hepatocellular carcinoma patients,” Journal of Cellular Biochemistry, vol. 120, no. 5, pp. 7539–7550, 2018.
View at: Publisher Site | Google Scholar
J. Yan, C. Zhou, K. Guo, Q. Li, and Z. Wang, “A novel seven-lncRNA signature for prognosis prediction in hepatocellular carcinoma,” Journal of Cellular Biochemistry, vol. 120, no. 1, pp. 213–223, 2019.
View at: Publisher Site | Google Scholar
Z. Wang, Q. Wu, S. Feng, Y. Zhao, and C. Tao, “Identification of four prognostic LncRNAs for survival prediction of patients with hepatocellular carcinoma,” PeerJ, vol. 5, p. e3575, 2017.
View at: Publisher Site | Google Scholar
Y. Ma, T. Luo, D. Dong, X. Wu, and Y. Wang, “Characterization of long non-coding RNAs to reveal potential prognostic biomarkers in hepatocellular carcinoma,” Gene, vol. 663, pp. 148–156, 2018.
View at: Publisher Site | Google Scholar
J. Sui, Y. Miao, J. Han et al., “Systematic analyses of a novel lncRNA-associated signature as the prognostic biomarker for hepatocellular carcinoma,” Cancer Medicine, vol. 7, no. 7, pp. 3240–3256, 2018.
View at: Publisher Site | Google Scholar
Z. Zhang, Y. Ouyang, Y. Huang et al., “Comprehensive bioinformatics analysis reveals potential lncRNA biomarkers for overall survival in patients with hepatocellular carcinoma: an on-line individual risk calculator based on TCGA cohort,” Cancer Cell International, vol. 19, no. 1, 2019.
View at: Publisher Site | Google Scholar
H. Zhang, B. Zhao, X. Wang, F. Zhang, and W. Yu, “LINC00511 knockdown enhances paclitaxel cytotoxicity in breast cancer via regulating miR-29c/CDK6 axis,” Life Sciences, vol. 228, pp. 135–144, 2019.
View at: Publisher Site | Google Scholar
G. Lu, Y. Li, Y. Ma et al., “Long noncoding RNA LINC00511 contributes to breast cancer tumourigenesis and stemness by inducing the miR-185-3p/E2F1/Nanog axis,” Journal of Experimental & Clinical Cancer Research, vol. 37, no. 1, p. 289, 2018.
View at: Publisher Site | Google Scholar
X. Zhao, Y. Liu, Z. Li et al., “Linc00511 acts as a competing endogenous RNA to regulate VEGFA expression through sponging hsa-miR-29b-3p in pancreatic ductal adenocarcinoma,” Journal of Cellular and Molecular Medicine, vol. 22, no. 1, pp. 655–667, 2018.
View at: Publisher Site | Google Scholar
J. Ding, C. Yang, and S. Yang, “LINC00511 interacts with miR-765 and modulates tongue squamous cell carcinoma progression by targeting LAMC2,” Journal of Oral Pathology & Medicine, vol. 47, no. 5, pp. 468–476, 2018.
View at: Publisher Site | Google Scholar
C. Li, H. Liu, J. Yang et al., “Long noncoding RNA LINC00511 induced by SP1 accelerates the glioma progression through targeting miR-124-3p/CCND2 axis,” Journal of Cellular and Molecular Medicine, vol. 23, no. 6, pp. 4386–4394, 2019.
View at: Publisher Site | Google Scholar
C. C. Sun, S. J. Li, G. Li, R. X. Hua, X. H. Zhou, and D. J. Li, “Long intergenic noncoding RNA 00511 acts as an oncogene in non-small-cell lung cancer by binding to EZH2 and suppressing p57,” Molecular Therapy - Nucleic Acids, vol. 5, no. 11, 2016.
View at: Publisher Site | Google Scholar
H. Li, L. Fang, X. Xiao, and L. Shen, “The expression and effects the CABYR-c transcript of CABYR gene in hepatocellular carcinoma,” Bulletin du Cancer, vol. 99, no. 3, pp. E26–E33, 2012.
View at: Publisher Site | Google Scholar
Z. Qian, M. Li, R. Wang et al., “Knockdown of CABYR-a/b increases chemosensitivity of human non-small cell lung cancer cells through inactivation of Akt,” Molecular Cancer Research, vol. 12, no. 3, pp. 335–347, 2014.
View at: Publisher Site | Google Scholar
L. Ju, X. Li, J. Shao, R. Lu, Y. Wang, and Z. Bian, “Upregulation of thyroid hormone receptor interactor 13 is associated with human hepatocellular carcinoma,” Oncology Reports, vol. 40, no. 6, pp. 3794–3802, 2018.
View at: Publisher Site | Google Scholar
J. Yao, X. Zhang, J. Li et al., “Silencing TRIP13 inhibits cell growth and metastasis of hepatocellular carcinoma by activating of TGF-β1/smad3,” Cancer Cell International, vol. 18, no. 1, 2018.
View at: Publisher Site | Google Scholar
Z. Weng, D. Wang, W. Zhao et al., “microRNA-450a targets DNA methyltransferase 3a in hepatocellular carcinoma,” Experimental and Therapeutic Medicine, vol. 2, no. 5, pp. 951–955, 2011.
View at: Publisher Site | Google Scholar
H. Reis, C. Pütter, D. A. Megger et al., “A structured proteomic approach identifies 14-3-3Sigma as a novel and reliable protein biomarker in panel based differential diagnostics of liver tumors,” Biochimica et Biophysica Acta, vol. 1854, no. 6, pp. 641–650, 2015.
View at: Publisher Site | Google Scholar
A. Boudreau, K. Tanner, D. Wang, F. C. Geyer, J. S. Reis-Filho, and M. J. Bissell, “14-3-3σ stabilizes a complex of soluble actin and intermediate filament to enable breast tumor invasion,” Proceedings of the National Academy of Sciences of the United States of America, vol. 110, no. 41, pp. E3937–E3944, 2013.
View at: Publisher Site | Google Scholar
A. Shiba-Ishii, Y. Kim, T. Shiozawa et al., “Stratifin accelerates progression of lung adenocarcinoma at an early stage,” Molecular Cancer, vol. 14, no. 1, 2015.
View at: Publisher Site | Google Scholar
Y. Hu, Q. Zeng, C. Li, and Y. Xie, “Expression profile and prognostic value of SFN in human ovarian cancer,” Bioscience Reports, vol. 39, no. 5, 2019.
View at: Publisher Site | Google Scholar
Y. F. Lian, Y. L. Huang, J. L. Wang et al., “Anillin is required for tumor growth and regulated by miR-15a/miR-16-1 in HBV-related hepatocellular carcinoma,” Aging, vol. 10, no. 8, pp. 1884–1901, 2018.
View at: Publisher Site | Google Scholar
S. Zhang, L. H. Nguyen, K. Zhou et al., “Knockdown of anillin actin binding protein blocks cytokinesis in hepatocytes and reduces liver tumor development in mice without affecting regeneration,” Gastroenterology, vol. 154, no. 5, pp. 1421–1434, 2018.
View at: Publisher Site | Google Scholar
Y. X. Song, J. X. Sun, J. H. Zhao et al., “Non-coding RNAs participate in the regulatory network of CLDN4 via ceRNA mediated miRNA evasion,” Nature Communications, vol. 8, no. 1, p. 289, 2017.
View at: Publisher Site | Google Scholar
T. Budden and N. A. Bowden, “MC1R CpG island regulates MC1R expression and is methylated in a subset of melanoma tumours,” Pigment Cell & Melanoma Research, vol. 32, no. 2, pp. 320–325, 2019.
View at: Publisher Site | Google Scholar
Q. Feng, J. E. Stern, S. E. Hawes, H. Lu, M. Jiang, and N. B. Kiviat, “DNA methylation changes in normal liver tissues and hepatocellular carcinoma with different viral infection,” Experimental and Molecular Pathology, vol. 88, no. 2, pp. 287–292, 2010.
View at: Publisher Site | Google Scholar
M. Bellodi-Privato, M. S. Kubrusly, J. T. Stefano et al., “Differential gene expression profiles of hepatocellular carcinomas associated or not with viral infection,” Brazilian Journal of Medical and Biological Research, vol. 42, no. 12, pp. 119–1127, 2009.
View at: Publisher Site | Google Scholar
G. Diaz, R. E. Engle, A. Tice et al., “Molecular signature and mechanisms of hepatitis D virus-associated hepatocellular carcinoma,” Molecular Cancer Research, vol. 16, no. 9, pp. 1406–1419, 2018.
View at: Publisher Site | Google Scholar
N. Iizuka, M. Oka, H. Yamada-Okabe et al., “Molecular signature in three types of hepatocellular carcinoma with different viral origin by oligonucleotide microarray,” International Journal of Oncology, vol. 24, no. 3, pp. 565–574, 2004.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2020 Xiaonan Zhao 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.

PDF Download Citation

Download other formats

Order printed copies

Views

641

Downloads

1285

Citations