Back to Journals » OncoTargets and Therapy » Volume 12

Identification Of Novel lncRNAs For Detection Of HBV-Associated Hepatocellular Carcinoma

Authors Li Y , Chen X, Huang H, Li G, Liao L, Yuan T, Deng S

Received 9 September 2019

Accepted for publication 2 November 2019

Published 27 November 2019 Volume 2019:12 Pages 10199—10211

DOI https://doi.org/10.2147/OTT.S230377

Checked for plagiarism Yes

Review by Single anonymous peer review

Peer reviewer comments 3

Editor who approved publication: Dr Carlos E Vigil



Yuwei Li,1,* Xia Chen,1,* Hengliu Huang,1,* Guangyao Li,2 Ling Liao,1 Tao Yuan,2 Shaoli Deng1

1Department of Clinical Laboratory, Daping Hospital, Army Medical University (Third Military Medical University), Chongqing 400042, People’s Republic of China; 2Department of Hepatobiliary Surgery, Army Medical University (Third Military Medical University), Chongqing 400042, People’s Republic of China

*These authors contributed equally to this work

Correspondence: Shaoli Deng
Department of Clinical Laboratory, Daping Hospital, Army Medical University (Third Military Medical University), No. 10, Changjiang Zhilu, DaPing, Yuzhong District, Chongqing 400042, People’s Republic of China
Tel +862368757602
Email [email protected]
Tao Yuan
Department of Hepatobiliary Surgery, Daping Hospital, Army Medical University (Third Military Medical University), No. 10, Changjiang Zhilu, DaPing, Yuzhong District, Chongqing 400042, People’s Republic of China
Tel +862368767966
Email [email protected]

Purpose: This study was conducted to investigate the differentially expressed profiles of long non-coding RNAs (lncRNAs) in HBV-associated HCC (HBV-HCC), which may serve as potential diagnostic biomarkers and therapeutic targets.
Methods: To examine the differentially expressed profiles of lncRNAs and mRNAs using microarray analysis, we collected 15 specimens: five HBV-associated HCC tissues, five paired adjacent peritumoral liver tissues (APLT), and five distant peritumoral liver tissues (DPLT). Then, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were performed to predict the biological roles and potential signaling pathways of these dysregulated mRNAs. In addition, lncRNA-mRNA co-expression network and signal transduction pathway network (Signal-net) were employed to further explore the potential target genes and roles of dysregulated lncRNAs in HBV-HCC pathogenesis. Finally, quantitative real-time PCR (qRT-PCR) was used to confirm the expression of six selected dysregulated lncRNAs.
Results: A total number of 719 lncRNAs and 3438 mRNAs were significantly more dysregulated in HBV-HCC tissues than in APLT and DPLT (fold change > 2, P < 0.05, FDR < 0.05). Additionally, 337 GO terms and 53 KEGG pathways were established to be significantly enriched. These dysregulated mRNAs were mainly enriched in metabolism-related biological processes. Additionally, lncRNA-mRNA coexpression network analysis showed that NONHSAT053785 is at the core of the network. Furthermore, the Signal-net analysis showed that CYP3A4 was gene with the highest degree. Finally, the data of five of the six selected differentially expressed lncRNAs were in agreement with the microarray data obtained by qRT-PCR verification.
Conclusion: Our study revealed the differentially expressed profiles of lncRNAs and mRNAs for HBV-HCC, and five novel dysregulated lncRNAs were identified in HBV-HCC tissues. The aforementioned dysregulated lncRNAs may represent potential diagnostic biomarkers and therapeutic targets of HBV-HCC, which needs to be validated in future studies.

Keywords: long non-coding RNA, hepatocellular carcinoma, microarray, biomarker

Introduction

Hepatocellular carcinoma (HCC), which accounts for almost 90% of primary liver cancer cases, is one of the most prevalent malignancies and leading causes of cancer-related death worldwide.14 Globally, an estimated 841,080 (4.7%) new HCC cases and 781,631 (8.2%) cancer-related deaths occur annually.3 Chronic hepatitis B virus (HBV) infection is the major risk factor in most high-risk HCC areas.5 Over 90% of the cases of HCC globally have a viral etiology, and their majority has been caused by chronic hepatitis B infection.6 Although HBV is a known hepatocarcinogen, the molecular mechanisms that determine HBV-associated HCC (HBV-HCC) development remain elusive.7 The lack of effective diagnostic and prognostic biomarkers leads limits the available options for HBV-HCC treatment. Thus, the further exploration of the underlying molecular pathogenesis of HBV-HCC and the subsequent development of new effective early diagnosis biomarkers and therapeutic targets for HBV-HCC are urgently needed.

The development of high-throughput RNA sequencing technology has enabled the discovery of thousands of non-coding RNA (ncRNA) genes. Generally, ncRNAs shorter than 200 nucleotides are usually described as small ncRNAs, whereas ncRNAs longer than 200 nucleotides are referred to as long non-coding RNAs (lncRNAs). The secrets of long noncoding RNAs are gradually elucidated.8 LncRNAs were recently shown to have abundant functions in multiple biologic processes. The overexpression, deficiency, or mutation of lncRNA genes has been implicated in numerous human diseases.9 Studies on cancer have revealed the prominent role of lncRNAs in various cancer forms by their promotion of growth and metastatic potential.10 It is now clear that lncRNAs are also critically involved in HCC pathogenesis. Some lncRNAs, such as the HULC, MALAT1 and H19, are overexpressed and participate in the tumorigenesis of HCC.1113 In 2007, by the use of an HCC-specific cDNA microarray platform, HULC was identified as the first lncRNA with highly specific upregulation in HCC.11 In another examination, Lin et al found that MALAT1 was closely related to HCC neoplasia and progression.13 In addition, more recent research indicates that lncRNA H19 can serve as a ceRNA to sponge miR-326 and modulate TWIST1 levels in HCC pathogenesis and may serve as a novel HCC target.12 The evidence of a correlation between the changes in the expression levels of lncRNAs and HCC development is increasing. However, currently, the research on lncRNAs in HBV-HCC is quite limited. Therefore, it is highly necessary and extremely urgent to investigate the association between lncRNAs and HBV-HCC, which may play an important role in the tumorigenesis and pathogenesis, and may potentially be used as a diagnostic marker or therapeutic target.

In past decades, the mechanism of HCC development has been investigated, and different clinical therapies to overcome this disease have been examined. In particular, the application of high-throughput techniques in HCC research has generated a vast amount of data, which has provided important resources for further exploration of the molecular mechanisms and identification of HCC-related molecules. Doubtless, these data will help to improve the diagnostic, prognostic, and therapeutic classification. Therefore, in the present study, we applied microarray technology to analyze the expression profiles of lncRNAs and mRNAs in hepatocellular carcinoma with HBV infection. Additionally, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were performed to predict the biological roles and potential signaling pathways of these differentially expressed mRNAs. Moreover, the lncRNA-mRNA coexpression network and signal transduction pathway network (Signal-net) were conducted to further explore the potential roles of differentially expressed lncRNAs in HBV-HCC pathogenesis and analysis, followed by quantitative real-time PCR (qRT-PCR) verification.

Materials And Methods

Patient Samples

From July 2016 to December 2016, samples were randomly collected for transcriptome microarray profile analysis from the Daping Hospital of the Army Medical University (Third Military Medical University), including five cases HBV-associated HCC tissues, paired adjacent peritumoral liver tissues (APLT, less than 2 cm away from the cancerous nest) and distant peritumoral liver tissues (DPLT, at least 5 cm away from the cancerous nest). Then, other 20 groups of specimens were collected from March 2016 to July 2017 for expressional validation (Table 1). Serological tests pathologically confirmed that all patients had HCC and were infected with HBV. The clinical specimens were frozen in liquid nitrogen immediately after collection. Written informed consent was obtained from the patients, and the study was approved by the Medical Ethics Committee of Daping Hospital.

Table 1 Information Of Clinical Specimens For Microarray And qRT-PCR

RNA Isolation

Total RNA was extracted from each liver tissue sample using TRIzol (Invitrogen, Carlsbad, CA, USA) per the manufacturer’s protocol. RNA quantity and quality were measured using a NanoDrop ND-1000 spectrophotometer (Thermo, Wilmington, DE, USA) and RNA integrity was evaluated by standard denaturing agarose gel electrophoresis.

Microarray Profiling

Microarray analysis was performed using GeneChip® Human Transcriptome Array 2.0 (Affymetrix, Santa Clara, CA, USA), which is designed for the global profiling of 245,349 protein coding transcripts and 40,914 non-protein coding transcripts from public transcriptome databases, including RefSeq, Ensembl, UCSC (known genes and lincRNA transcripts), as well as landmark publications and was carried out by Beijing Cnkingbio Biotechnology Co, Ltd (Beijing, China). Affymetrix® Transcriptome Analysis Console (TAC) Software was employed to analyze the results from Human Transcriptome Array 2.0.

GO Enrichment Analysis

The Gene Ontology (GO) is a structured and controlled vocabulary of terms. The terms are subdivided in three non-overlapping ontologies, Molecular Function (MF), Biological Process (BP) and Cellular Component (CC), and now is used widely for annotating genes and gene products.14 To detect the overrepresentation of Gene-Ontology annotations among the large list of genes from high-throughput experiments such as microarray that is effective in targeting downstream events and molecules. Parent-Child-Intersection method was applied for enrichment analysis and Benjamini-Hochberg was utilized for multiple tests correction, and curated association (ie, all evidence codes except IEA, ND, and NR) was employed for enrichment analysis.15

KEGG Pathway Analysis

Kyoto Encyclopedia of Genes and Genomes (KEGG) is a knowledge base utilized for systematic analysis of gene functions, linking genomic information with higher order functional information, which is now widely used for pathway-related analysis.16 Hypergeometric distribution was used to calculate the pathway enrichment, and FDR (false discovery rate) was implemented to adjust the P-values in the multiple comparisons.

Co-Expression Analysis

The co-expression network between lncRNA and genes was constructed by combining the expression profile data of differentially expressed lncRNA components with the expression profile data of differentially expressed genes. Coexpression network can search for the core regulated genes in the network and locate the functions and signal transduction pathways in which lncRNA may participate through genes it regulates. The co-expression network clearly reflects the relationship between lncRNA and genes as well as the overall interaction, enabling the discovery of the regulatory relationship between lncRNA and genes and the exploration of its possible functional effects.

Global Signal Transduction Network

After parsing the whole KEGG database, all study genes involved pathways were extracted, and the study pathway network was generated with the help of the pathway topology in the KEGG database. The specific gene network of one pathway was generated based on the pathway topology analysis, and the study gene network was generated after mapping to the generated reference KEGG gene network (ie a collection of gene networks from each pathway).

Quantitative Real-Time PCR (qRT-PCR)

The expression levels of selected lncRNAs were validated by qRT-PCR. The specimens of 20 groups were used for further validation. cDNA was synthesized by EvoScript Universal cDNA Master RT reagent Kit (Roche, Germany) under the manufacturer’s document. The lncRNA expression level was quantified using TaqMan® Fast Advanced Master Mix (ABI, USA) on the CFX96TM Real-Time system (BIO-RAD, USA). All experiments were performed in triplicate (glyceraldehyde 3-phosphate dehydrogenase (GAPDH) as the internal control). Finally, the 2−∆∆Ct method was used to calculate the relative expression. The primers, whose sequences are listed in Table 2, were obtained from Sangon Biotech (Shanghai) Co, Ltd (Shanghai, China).

Table 2 Primer Sequences Of The qRT-PCR Used In This Study

Statistical Analysis

Statistical analyses were performed using GraphPad Prism 5.0 (San Diego, CA, USA). The qRT-PCR data were statistically analyzed by Student’s t-test with SPSS 21.0; differences were considered statistically significant at P < 0.05.

Results

lncRNA And mRNA Expression Profiles In HBV-HCC

Fifteen specimens (five HCC, five APLT, and five DPLT) were subjected to microarray analysis of lncRNA-mRNA. We found that 1035 lncRNAs and 4313 mRNAs were differentially expressed between HCC and APLT, whereas 1126 lncRNAs and 4634 mRNAs were differentially expressed between HCC and DPLT (fold change > 2, P < 0.05, FDR < 0.05; Table 3). The Venn diagram (Figure 1GH) showed that 719 of the differentially expressed lncRNAs and 3438 of the differentially expressed mRNAs were identified in both sets. They were considered to be candidate tumor-derived RNAs and subjected to subsequent analysis. The differential lncRNA expression was jointly evaluated by heatmap, volcano plot, and scatter plot, and categorized in the HCC, APLT, and DPLT (Figure 1AF).

Table 3 Number Of Significant Dysregulated Genes On Microarray In Groups

Figure 1 Differentially expressed lncRNAs (fold change >2, P < 0.05) in the microarray profiles. (A) Heatmap of five paired HCC and APLT. (B) Heatmap of five paired HCC and DPLT. (C) Volcano plot of five paired HCC and APLT. (D) Volcano plot of five paired HCC and DPLT. (E) Scatter plot of five paired HCC and APLT. (F) Scatter plot of five paired HCC and DPLT. (G) Venn diagram indicating that 719 lncRNAs are consistently in both microarray analyses. (H) Venn diagram indicating that 3438 mRNAs are consistently in both microarray analyses.

GO Enrichment And KEGG Pathway Analysis

To predict the functions of the differentially expressed genes, we conducted functional enrichment analysis of the 3438 dysregulated mRNAs, and 337 GO terms were identified as significantly enriched. As can be seen in Figure 2, the GO analysis revealed that the most enriched GO terms were small-molecule metabolic process (ontology: biological process, GO: 0044281), extracellular exosome (ontology: cellular component, GO: 0070062) and poly(A) RNA binding (ontology: molecular function, GO:0044822). GO network analysis was performed to generate the interaction network among these enriched GO terms, of which the vitamin metabolic process (degree = 3) had the highest degrees (Figure 3). Furthermore, the KEGG pathway analysis indicated 53 KEGG pathways were significantly enriched and the dysregulated mRNAs were mainly involved in metabolism-related biological processes, such as the regulation of systemic lupus erythematosus, complement and coagulation cascades, valine, leucine, and isoleucine degradation, alcoholism, and carbon metabolism (Figure 4). We performed pathway network analysis to generate an interaction network covering significantly changed pathways (Figure 5). The alanine, aspartate and glutamate metabolism (degree = 19) showed the highest degrees, suggesting that it may plays a core role.

Figure 2 GO terms for the dysregulated mRNAs (P < 0.01). (A) GO terms for the dysregulated mRNAs in biological process. (B) GO terms for the dysregulated mRNAs in cellular component. (C) GO terms for the dysregulated mRNAs in molecular function.

Figure 3 GO network. Significantly changed GO functions were connected in a GO network to show the relationship between these functions.

Figure 4 KEGG pathways analysis (P < 0.05). Major pathways of the dysregulated mRNAs.

Figure 5 Pathway network. Significantly changed pathways were connected in a pathway network to show the relationship between these pathways.

lncRNA-mRNA Coexpression Network Analysis

To investigate the potential roles of these dysregulated lncRNAs in HCC, we selected dysregulated mRNAs with both significantly changed GOs and pathways through GO and KEGG analysis, and constructed a coexpression network of these mRNAs with dysregulated lncRNAs. As can be observed in Figure 6, the whole coexpression network profile consisted of 336 network nodes and 416 connections among 127 differentially expressed lncRNAs and 209 differentially expressed mRNAs. There were 350 positive and 66 negative interactions within the network. Moreover, our data showed that one lncRNA may correlate with 1–28 mRNAs. Among them, NONHSAT053785, RP1-81D8.6, UGT2B25P, HIST1H2APS4, and RP11-371E8.2 were at the core of the co-expression network and were related to at least 20 mRNAs. The top five mRNAs related to NONHSAT053785, RP1-81D8.6, UGT2B25P, HIST1H2APS4, and RP11-371E8.2 ranked by their correlation value, are listed in Table 4.

Table 4 Top 5 mRNAs Ranked By Correlation After Coexpression Network Analysis

Figure 6 lncRNA-mRNA coexpression network. Square nodes represent dysregulated lncRNAs, round nodes represent dysregulated mRNAs. The dashed lines between lncRNAs and mRNAs indicate a negative correlation, while the solid lines indicate a positive correlation.

Global Signal Transduction Network

The signal transduction network analysis was used to define the core key regulatory mRNA genes which had a strong capacity to modulate adjacent mRNA genes. Based on the significantly regulated GOs and pathways, we performed Signal-net analysis to screen the important candidate genes involved in HCC associated with hepatitis B virus infection. As illustrated in Figure 7, 995 genes were obtained in the Signal-net. The resulting network indicated that genes with high degrees were CYP3A4, PRKX, CYP1A1, CYP2C9, CYP1A2, CYP2A6, MAPK1, MAPK3, PIK3R1, and ITGA2 (Table 5).

Table 5 Top Ten Genes Ranked By Degree After Analysis Of Signal-Net

Figure 7 Global signal transduction network. The area of the circle represents the degree. The larger area indicates a higher degree. Interaction between the genes is denoted as follows: a: activation, b(a): binding/association, c: compound; dep: dephosphorylation, dis: dissociation, ex: expression, ind(e): indirect effect, inh: inhibition, mis(i): missing interaction, pho: phosphorylation, rep: repression, s(c): state change, and ubi: ubiquitination.

The Validation Of Novel lncRNAs That May Become Potential Biomarkers

Based on the high normalized intensity in the raw data, the results of bioinformatics analysis, high fold change (fold change > 5), significant P-value (P < 0.01) and sorted by adjustment of the P-value, we selected six lncRNAs for qRT-PCR verification to confirm the reliability of the microarray data. Sorted by the above parameters, the six lncRNAs were all ranked in the front position. Among these six dysregulated lncRNAs, three (ENST00000514608.1, NONHSAT138156, and ENST00000560295.1) were upregulated, whereas the other three (NONHSAT024276, NONHSAT124357.2, and NONHSAT053785) were downregulated. Then, we analyzed their expression levels by qRT-PCR in another 60 specimens (20 HCC, 20 APLT, and 20 DPLT samples). The result of qRT-PCR indicated that five lncRNAs were consistent with the microarray analysis and showed the same trends of up- or downregulation for each lncRNA with statistical significance. ENST00000514608.1, NONHSAT138156 were significantly upregulated and NONHSAT024276, NONHSAT124357, NONHSAT053785 were significantly downregulated in HCC cases (all P < 0.05 and detection rate = 100%). No significant differences in ENST00000560295.1 were found among the three groups. The △CT of each lncRNA is presented in Figure 8.

Figure 8 The △CT of the Screening lncRNAs in HCC, APLT, and DPLT by qRT-PCR. (A) △CT of lncRNA ENST00000514608.1. (B) △CT of lncRNA NONHSAT138156. (C) △CT of lncRNA ENST00000560295.1. (D) △CT of lncRNA NONHSAT024276. (E) △CT of lncRNA NONHSAT124357.2. (F) △CT of lncRNA NONHSAT053785. HCC group, n = 20; APLT group, n = 20; DPLT, n = 20. *P < 0.05, **P < 0.005,***P < 0.001,****P < 0.0001.

Discussion

In recent years, numerous studies have shown that HBV is involved in the carcinogenesis, invasion, and metastasis of liver cells, playing a key role in the occurrence and development of liver cancer.1719 LncRNA was reported to be involved in the development of various hepatocellular carcinomas.20,21 LncRNA was also found to participate in HBV-HCC pathogenesis; HULC promoted hepatoma cell proliferation via downregulating p18, thereby exacerbating the malignancy of HCC.22 More studies on lncRNAs in HBV-related HCC are necessary to understand the process of HCC pathogenesis. Currently, high-throughput gene chip technology enabled the exploration of the occurrence and development of diseases at the whole-genome or transcriptome level, and has been widely used in some analysis of disease gene expression profile, gene cloning, and the search for disease-specific molecular markers.2325 Thus, in our study, we utilized microarray and bioinformatics to seek crucial lncRNAs and reveal their potential functions in association with HBV-HCC. We identified six crucial lncRNAs and used qRT-PCR to validate their expression levels.

In our present research, using microarray analysis, we found aberrant expression in 719 lncRNAs and 3438 mRNAs. Using bioinformatics analyses, including GO, KEGG pathway, lncRNA-mRNA network, and signal-net analyses, we studied the potential role of these identified genes and their relationship with the development of HCC. GO analysis identified the significant biological functions of the aberrantly expressed protein-coding genes. As a result, 337 GO terms were significantly enriched, which were related to the biological processes, cellular component, and molecular functions. The most significant GO terms were the small-molecule metabolic process (GO:0044281), extracellular exosome (GO: 0070062), and poly(A) RNA binding (GO:0044822), respectively, in these three aspects, which were all important pathophysiological and morphological changes in HCC. HCC patients with high PKM2 expression suffered from more severe metabolic damage, transcriptional regulation imbalance, and poor prognosis.26 The pathway analysis revealed 53 enriched pathways, including the systemic lupus erythematosus signaling pathway, complement and coagulation cascades signaling pathway, and valine, leucine, and isoleucine degradation signaling pathway. To further investigate the complicated regulatory network between lncRNAs and mRNAs, we conducted the lncRNA-mRNA coexpression network analysis based on the aberrantly expressed lncRNAs and mRNAs. A total number of 127 lncRNAs and 209 mRNAs were included in this network, involving 336 network nodes and 416 connections. These numerous connections, which were lncRNA-mRNA coexpressed pairs (positive or negative correlation), suggested that the aberrantly expressed lncRNAs interacted with aberrantly expressed mRNAs. NONHSAT053785 was at the core of the co-expression network and was associated with 28 mRNAs, which had not been reported up to that time. In the future, we intend to study its role in HCC cell proliferation and recurrence. To study the regulatory relationships among these mRNAs, we used the intersecting genes between GO analysis and pathway analysis to perform the signal-net analysis. Our results showed that genes with high degrees were CYP3A4, PRKX, CYP1A1, CYP2C9, CYP1A2, CYP2A6, MAPK1, MAPK3, PIK3R1, and ITGA2. We selected lncRNAs in the genomic interzone, and finally determined six lncRNAs based on the calibrated P-value, which were validated by qRT-PCR analysis. The result of qRT-PCR indicated that the data of five lncRNAs were consistent with those of the microarray analysis. ENST00000514608.1 and NONHSAT138156 were significantly upregulated, whereas NONHSAT024276, NONHSAT124357, and NONHSAT053785 were significantly downregulated in HCC cases. On the other hand, ENST00000560295.1 (also called CASC7) showed no significant difference among the three groups. However,lncRNA CASC7 was reported to inhibit the proliferation and migration of colon cancer cells via microRNA-21 suppression.27 Apart from CASC7, all other (five) lncRNAs were novel and not yet reported. LncRNA SNHG20 expression was more upregulated in tumor tissues of HBV-positive HCC patients than in tumor tissues of HBV-negative HCC patients.28 Nonetheless, more experiments are required to establish whether the differences in the screened lncRNAs are indeed associated with HBV. The current study on long non-coding RNA provided novel insights into the expression profile in liver cancer and predicted the potential roles of long non-coding RNAs in hepatocellular carcinogenesis.

Despite the interesting profiling of lncRNA expression in HBV-HCC, several limitations in our study should be acknowledged. First of all, the small numbers of the samples in the microarray assay limited its reliability and biases might have occurred, affecting the results of the further analyses. Moreover, our study predicted mainly the functions of differentially expressed lncRNAs, but we were unable to exactly reveal the regulatory roles of these lncRNAs. In the future, we aim to further investigate that subject, focusing on lncRNA molecular mechanisms in HCC. At the same time, we will also focus on the differential expression of HBV-HCC and non HBV-HCC.

In summary, in this study, we comprehensively analyzed the profiles of key lncRNA in the progression of HBV-HCC. We discovered new lncRNAs that may play an important role in the progression of HBV-HCC, such as NONHSAT053785. Nevertheless, further studies on the role and mechanism of the above key lncRNA in the occurrence and development of HBV-HCC are expected to provide new molecular and potential new targets for the diagnosis, treatment and prognosis of HBV-HCC.

Conclusion

Our study identified five novel dysregulated lncRNAs in HCC with HBV infection. Among them, NONHSAT053785 was at the core of the lncRNA-mRNA co-expression network, which may play an important role in the tumorigenesis of HBV-HCC. Therefore, our study provides an important foundation for further exploration of the molecular mechanisms of HCC and identification of the related molecules.

Abbreviations

lncRNA, long non-coding RNA; HCC, hepatocellular carcinoma; HBV, hepatitis B virus; HBV-HCC, HBV-associated HCC; HULC, highly upregulated in liver carcinoma; MALAT1, metastasis associated lung adenocarcinoma transcript 1; ncRNA, non-coding RNA; APLT, adjacent peritumoral liver tissues; DPLT, distant peritumoral liver tissues; GO, Gene Ontology; KEGG, Kyoto Encyclopedia of Genes and Genomes; MF, molecular function; BP, biological process; CC, cellular component; Signal-net, signal transduction pathway network; FC, fold change; FDR, false discovery rate; qRT-PCR, quantitative real-time PCR.

Ethical Approval And Informed Consent

All procedures performed in studies involving human participants were in accordance with the ethical standards of the institutional and/or national research committee and with the 1964 Helsinki Declaration and its later amendments or comparable ethical standards. Informed consent was obtained from all individual participants, and the study was approved by the Medical Ethics Committee of Daping Hospital.

Availability Of Data And Material

The datasets analyzed during the current study are available from the corresponding author on reasonable request.

Acknowledgments

This work was supported by the major project of Army Medical University (2018 XLC 2027) and excellent talents fund of Daping Hospital. Yuwei Li, Xia Chen and Hengliu Huang are the co-first authors of this article.

Author Contributions

All authors contributed to data analysis, drafting, and revising the article, gave final approval of the version to be published, and agree to be accountable for all aspects of the work.

Disclosure

The authors report no conflicts of interest in this work.

References

1. Forner A, Reig M, Bruix J. Hepatocellular carcinoma. Lancet. 2018;391(10127):1301–1314. doi:10.1016/S0140-6736(18)30010-2

2. Llovet JM, Zucman-Rossi J, Pikarsky E, et al. Hepatocellular carcinoma. Nat Rev Dis Primers. 2016;2:16018.

3. Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;68(6):394–424. doi:10.3322/caac.v68.6

4. Xu RH, Wei W, Krawczyk M, et al. Circulating tumour DNA methylation markers for diagnosis and prognosis of hepatocellular carcinoma. Nat Mater. 2017;16(11):1155–1161. doi:10.1038/nmat4997

5. Chimed T, Sandagdorj T, Znaor A, et al. Cancer incidence and cancer control in Mongolia: results from the National Cancer Registry 2008–12. Int J Cancer. 2017;140(2):302–309. doi:10.1002/ijc.v140.2

6. MacLachlan JH, Locarnini S, Cowie BC. Estimating the global prevalence of hepatitis B. Lancet. 2015;386(10003):1515–1517. doi:10.1016/S0140-6736(15)61116-3

7. Seehawer M, Heinzmann F, D’Artista L, et al. Necroptosis microenvironment directs lineage commitment in liver cancer. Nature. 2018;562(7725):69–75. doi:10.1038/s41586-018-0519-y

8. Yang X, Xie X, Xiao YF, et al. The emergence of long non-coding RNAs in the tumorigenesis of hepatocellular carcinoma. Cancer Lett. 2015;360(2):119–124. doi:10.1016/j.canlet.2015.02.035

9. Quinn JJ, Chang HY. Unique features of long non-coding RNA biogenesis and function. Nat Rev Genet. 2016;17(1):47–62. doi:10.1038/nrg.2015.10

10. Lorenzen JM, Thum T. Long noncoding RNAs in kidney and cardiovascular diseases. Nat Rev Nephrol. 2016;12(6):360–373. doi:10.1038/nrneph.2016.51

11. Panzitt K, Tschernatsch MM, Guelly C, et al. Characterization of HULC, a novel gene with striking up-regulation in hepatocellular carcinoma, as noncoding RNA. Gastroenterology. 2007;132(1):330–342.

12. Wei LQ, Li L, Lu C, Liu J, Chen Y, Wu H. Involvement of H19/miR-326 axis in hepatocellular carcinoma development through modulating TWIST1. J Cell Physiol. 2019;234(4):5153–5162. doi:10.1002/jcp.27319

13. Lin R, Maeda S, Liu C, Karin M, Edgington TS. A large noncoding RNA is a marker for murine hepatocellular carcinomas and a spectrum of human carcinomas. Oncogene. 2007;26(6):851–858. doi:10.1038/sj.onc.1209846

14. Du Plessis L, Skunca N, Dessimoz C. The what, where, how and why of gene ontology–a primer for bioinformaticians. Brief Bioinform. 2011;12(6):723–735. doi:10.1093/bib/bbr002

15. Bauer S, Gagneur J, Robinson PN. GOing Bayesian: model-based gene set analysis of genome-scale data. Nucleic Acids Res. 2010;38(11):3523–3532. doi:10.1093/nar/gkq045

16. Kanehisa M, Goto S, Sato Y, Furumichi M, Tanabe M. KEGG for integration and interpretation of large-scale molecular data sets. Nucleic Acids Res. 2012;40(Database issue):D109–D114. doi:10.1093/nar/gkr988

17. Tang L, Covert E, Wilson E, Kottilil S. Chronic hepatitis B infection: a review. JAMA. 2018;319(17):1802–1813. doi:10.1001/jama.2018.3795

18. Kanda T, Goto T, Hirotsu Y, Moriyama M, Omata M. Molecular mechanisms driving progression of liver cirrhosis towards hepatocellular carcinoma in chronic hepatitis B and C infections: a review. Int J Mol Sci. 2019;20(6). doi:10.3390/ijms20061358

19. Feld JJ, Krassenburg L. What comes first: treatment of viral hepatitis or liver cancer? Dig Dis Sci. 2019;64(4):1041–1049. doi:10.1007/s10620-019-05518-5

20. Abbastabar M, Sarfi M, Golestani A, Khalili E. lncRNA involvement in hepatocellular carcinoma metastasis and prognosis. Excli J. 2018;17:900–913. doi:10.17179/excli2018-1541

21. Heo MJ, Yun J, Kim SG. Role of non-coding RNAs in liver disease progression to hepatocellular carcinoma. Arch Pharm Res. 2019;42(1):48–62. doi:10.1007/s12272-018-01104-x

22. Du Y, Kong G, You X, et al. Elevation of highly up-regulated in liver cancer (HULC) by hepatitis B virus X protein promotes hepatoma cell proliferation via down-regulating p18. J Biol Chem. 2012;287(31):26302–26311. doi:10.1074/jbc.M112.342113

23. Woo HG, Park ES, Thorgeirsson SS, Kim YJ. Exploring genomic profiles of hepatocellular carcinoma. Mol Carcinog. 2011;50(4):235–243. doi:10.1002/mc.v50.4

24. Marquardt JU, Galle PR, Teufel A. Molecular diagnosis and therapy of hepatocellular carcinoma (HCC): an emerging field for advanced technologies. J Hepatol. 2012;56(1):267–275. doi:10.1016/j.jhep.2011.07.007

25. Hung JH, Weng Z. Analysis of microarray and RNA-seq expression profiling data. Cold Spring Harb Protoc. 2017;2017(3). doi:10.1101/pdb.top093104

26. Lv WW, Liu D, Liu XC, et al. Effects of PKM2 on global metabolic changes and prognosis in hepatocellular carcinoma: from gene expression to drug discovery. BMC Cancer. 2018;18(1):1150. doi:10.1186/s12885-018-5023-0

27. Zhang Z, Fu C, Xu Q, Wei X. Long non-coding RNA CASC7 inhibits the proliferation and migration of colon cancer cells via inhibiting microRNA-21. Biomed Pharmacother. 2017;95:1644–1653. doi:10.1016/j.biopha.2017.09.052

28. Tu W, Yang Y, Song Y, Zhu W. Hepatitis B virus x protein accelerated the proliferation of hepatocellular carcinoma cell through lncRNA SNHG20/PTEN pathway. J Biochem. 2019;165(5):423–431.

Creative Commons License © 2019 The Author(s). This work is published and licensed by Dove Medical Press Limited. The full terms of this license are available at https://www.dovepress.com/terms.php and incorporate the Creative Commons Attribution - Non Commercial (unported, v3.0) License. By accessing the work you hereby accept the Terms. Non-commercial uses of the work are permitted without any further permission from Dove Medical Press Limited, provided the work is properly attributed. For permission for commercial use of this work, please see paragraphs 4.2 and 5 of our Terms.