Potential lncRNA Biomarkers for HBV-Related Hepatocellular Carcinoma Diagnosis Revealed by Analysis on Coexpression Network

Nong, Shunqiang; Chen, Xiaohao; Wang, Zechen; Xu, Guidan; Wei, Wujun; Peng, Bin; Zhou, Lü; Wei, Liuzhi; Zhao, Jingjie; Wei, Qiuju; Deng, Yibin; Meng, Lingzhang

doi:https://doi.org/10.1155/2021/9972011

BioMed Research International

On this page

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

Special Issue

Interplay between Immune Cells and their Microenvironment Niches 2021

View this Special Issue

Research Article | Open Access

Volume 2021 | Article ID 9972011 | https://doi.org/10.1155/2021/9972011

Potential lncRNA Biomarkers for HBV-Related Hepatocellular Carcinoma Diagnosis Revealed by Analysis on Coexpression Network

Shunqiang Nong,^1,2,3,4Xiaohao Chen,^1,2,3,4Zechen Wang,¹Guidan Xu,^1,2,3,4Wujun Wei,^1,2,3,4Bin Peng,^1,2,3,4Lü Zhou,^1,2,3,4Liuzhi Wei,¹Jingjie Zhao,⁵Qiuju Wei,¹Yibin Deng,^2,3,4and Lingzhang Meng¹

Academic Editor: Hannes Stockinger

Received19 Mar 2021

Accepted31 Aug 2021

Published15 Oct 2021

Abstract

Background. Increasing evidence demonstrated that long noncoding RNA (lncRNA) could affect inflammatory tumor immune microenvironment by modulating gene expression and could be used as a biomarker for HBC-related hepatocellular carcinoma (HCC) but still needs further research. The aim of the present study was to determine an lncRNA signature for the diagnosis of HBV-related HCC. Methods. HBV-related HCC expression profiles (GSE55092, GSE19665, and GSE84402) were abstracted from the GEO (Gene Expression Omnibus) data resource, and R package limma and RobustRankAggreg were employed to identify common differentially expressed genes (DEGs). Using machine learning, optimal diagnostic lncRNA molecular markers for HBV-related HCC were identified. The expression of candidate lncRNAs was cross-validated in GSE121248, and an ROC (receiver operating characteristic) curve of lncRNA biomarkers was carried out. Additionally, a coexpression network and functional annotation was built, after which a PPI (protein-protein interaction) network along with module analysis were conducted with the Cytoscape open source software. Result. A total of 38 DElncRNAs and 543 DEmRNAs were identified with a fold change larger than 2.0 and a value < 0.05. By machine learning, AL356056.2, AL445524.1, TRIM52-AS1, AC093642.1, EHMT2-AS1, AC003991.1, AC008040.1, LINC00844, and LINC01018 were screened out as optional diagnostic lncRNA biosignatures for HBV-related HCC. The AUC (areas under the curve) of the SVM (support vector machine) model and random forest model were 0.957 and 0.904, respectively, and the specificity and sensitivity were 95.7 and 100% and 94.3 and 86.5%, respectively. The results of functional enrichment analysis showed that the integrated coexpressed DEmRNAs shared common cascades in the p53 signaling pathway, retinol metabolism, PI3K-Akt signaling cascade, and chemical carcinogenesis. The integrated DEmRNA PPI network complex was found to be comprised of 87 nodes, and two vital modules with a high degree were selected with the MCODE app. Conclusion. The present study identified nine potential diagnostic biomarkers for HBV-related HCC, all of which could potentially modulated gene expression related to inflammatory conditions in the tumor immune microenvironment. The functional annotation of the target DEmRNAs yielded novel evidence in evaluating the precise functions of lncRNA in HBV-related HCC.

1. Introduction

Hepatocellular carcinoma (HCC) is a primary tumor that originates from hepatocytes, which is one of the most frequent malignant tumors and ranks sixth in incidence and fourth in mortality globally [1]. Chronic hepatitis B virus (HBV) infection constitutes the primary risk factor in most high-risk HCC areas. Multiple studies have documented that HBV participates in the carcinogenesis and invasion as well as metastasis of liver cells, playing a crucial role in the onset and development of liver cancer [2–4]. Given the lack of distinct clinical manifestations, most diagnosed HBV-linked HCC patients are in the advanced stage and have a poor prognosis. Although some serum biomarkers for early diagnosis have been reported in the literature, the results were not very satisfactory [5–7]. Therefore, finding new molecular markers for the early diagnosis, treatment, and prognostic evaluation of HBV-related HCC is urgently required.

Long noncoding RNAs (lncRNAs) constitute a large group of transcribed RNA biomolecules that are longer than 200 nucleotides, but their protein-coding ability is limited since they lack functional open readings [8]. lncRNAs modulate gene expression of posttranslational, translational, posttranscriptional, and transcriptional as well as epigenetic levels in diverse ways [9, 10]. Moreover, lncRNAs have been documented to participate actively in the modulation of diverse aspects of tumor progress consisting of growth and metastasis, as well as recurrence [11, 12]. For example, HULC, a well-explored lncRNA that is associated with HBV infection along with HCC tumor growth, has been reported to be upregulated in HCC and is linked to its grades, metastasis, and drug resistance [11]. Upregulated-lncRNA ZEB1-AS1 is correlated with tumor growth along with metastasis in HCC via modulating the expression of epithelial-mesenchymal transition-mediated markers, and patients with high expressions of ZEB1-AS1 exhibit high metastatic relapse, as well as dismal survival [13]. lncRNA MALAT1 is responsible for proliferation, migration, and infiltration, as well as apoptosis through its interactions with miR-200a in hypoxic Hep3B cells. Evidence of lncRNA dysregulation, along with the occurrence and development of HCC, has been gradually increasing. However, studies on HBV-related HCC are relatively limited.

Herein, three microarray expression data of a large number of HBV-related HCC patients were abstracted from the NCBI-GEO data resource, and the DElncRNAs (differentially expressed lncRNAs), as well as DEmRNAs (differentially expressed mRNAs), between HBV-related HCC and nontumor tissues were screened. Feature selection along with classification models was conducted to identify optimal lncRNAs with diagnostic value for HBV-related HCC. GO (Gene Ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) analyses, as well as the PPI (protein-protein interaction) and hub gene identification, were also performed. Accordingly, this study provided novel biomarkers for diagnosis and identified promising treatment targets for the treating HBV-related HCC.

2. Methods and Materials

2.1. Microarray Data Analysis

The NCBI-GEO data resource (https://www.ncbi.nlm.nih.gov/geo/) is a free repository of microarray gene patterns and next-generation sequencing, from which HBV-related hepatocellular carcinomas and nontumor or neighboring tissue gene expression profiles of GSE55092, GSE19665, and GSE84402 were abstracted. The three microarray data were all based on GPL570 platforms (HG-U133_Plus_2) and the Affymetrix Human Genome U133 Plus 2.0 Array.GSE55092, GSE19665, and GSE84402 contained 91, 5, and 13 cases of nontumor tissues from hepatitis B patients, respectively, and 49, 5, and 13 cases of HBV-associated HCC tissues, respectively. In order to acquire the international standard gene name, the probe was reannotated to determine the expression levels of lncRNA, as well as mRNA potentially related to HBV-related HCC. Briefly, the chip sequence file was first downloaded in FASTA format from the Affymetrix web page (www. http://affymetrix.com). Second, we aligned the probe sets of the HG-U133_Plus_2.0 array to the human genome (GRCh38), as well as the lncRNA gene sequence from GENCODE (https://www.gencodegenes.org/) (release 30) using the SeqMap tool with no mismatches. Finally, 4,971 probes were matched with 3,750 lncRNAs, and 31,475 probes were matched with 15,435 mRNAs in the HG-U133 Plus 2.0 microarray. Subsequently, the gene expression data that was subjected to log2 transformation was normalized between the “Arrays” function in the R limma package (http://www.bioconductor.org/).

2.2. Screening for Differentially Expressed mRNAs and lncRNAs

The DElncRNAs along with the DEmRNAs between HBV-linked HCC and nontumor tissues were computed through the limma R package V3.5.2 in R software in each dataset. To avoid missing those that are rare but important in the regulation of pathways, we set the threshold of the expression of DElncRNAs, as well as DEmRNAs, to and . Afterwards, the RobustRankAggreg (RRA) R package was used to integrate and analyze the three gene lists that were subsequently sorted via the logFC value. The package RobustRankAggreg (RRA) is usually used for comparing sequenced gene lists based on reviewing each gene’s ranking position (in case that each gene is randomly arranged). We exported and saved the lists of remarkably upregulated, as well as downregulated, genes as Excel files. After that, the hierarchical clustering of DEmRNAs and DElncRNAs was carried out with the R "heatmap” package.

2.3. Identifying the Optimal Diagnostic lncRNA Molecular Markers for HBV-Related HCC

The random forest analysis was employed to perform the feature selection procedures which were conducted to determine the optimal diagnostic lncRNA biosignatures for HBV-related HCC. Additionally, the importance value of every DElncRNA was ranked as per the MDG (Mean Decrease Gini). The optimal number of features was then determined by subsequently adding one DElncRNA at a time in a top-down forward-wrapper manner. Through applying the SVM (support vector machine) at each increment, we explored the accuracy and uncovered the optimal diagnostic lncRNA molecular markers for HBV-related HCC. Selection of the optimal DElncRNAs in HBV-related HCC, as well as nontumor tissue, was then performed. The random forest model and the SVM model were established using the randomForest package and the e1071 package, respectively.

2.4. Cross-Validation

Furthermore, another microarray dataset was collected from the GEO database under accession number GSE121248 (platform: GPL570 [HG-U133_Plus_2] Affymetrix Human Genome U133 Plus 2.0 Array) to serve as the validation set. GSE121248 contained 37 cases of nontumor tissues from hepatitis B patients as well as 70 cases of HBV-linked HCC tissues. The diagnostic ability of each model and lncRNA molecular markers was explored by the ROC (receiver operating characteristic), AUC (area under curve), sensitivity, and specificity in the validation set.

2.5. DEmRNA-DElncRNA Coexpressed Analysis

The coexpression network searches for the core modulated genes in the network and locates the functions along with the signal transduction cascades in which lncRNA participates via the genes it modulates. Here, the coexpression network evidently exhibited the association of the lncRNAs with the genes and the overall cross-talk, facilitating the discovery of the modulatory correlation of the lncRNAs with the genes, as well as the evaluation of its potential functional impacts. The correlation between DElncRNAs and DEmRNAs was assessed by the pairwise Pearson correlation coefficient. Besides, the threshold for DElncRNA-DEmRNA coexpression pairs was adjusted to and . Afterwards, the HBV-related HCC-specific DElncRNA-DEmRNA coexpression network was generated using the Cytoscape platform (http://cytoscape.org).

2.6. Functional and Signal Pathway Enrichment Analysis

The DAVID version 6.8 data resource (http://david.ncifcrf.gov) is a web-accessible program for integrating biological data that consists of analysis tools, which provide an exhaustive set of functional annotation data of extensive lists of proteins or genes for researchers to comprehend the biological features. In order to reveal the biological roles along with the potential cascades of the genes coexpressed using optimal lncRNAs with a diagnostic value for HBV-related HCC, GO classification and KEGG pathway enrichment were conducted using DAVID, in which signified significant difference.

2.7. Protein-Protein Interaction Network and Module Analysis

The STRING online data resource (https://string-db.org/) was employed to develop the PPI network for screening genes, where the interaction score was set to >0.4. Cytoscape was employed to visualize the results, and the MCODE (Molecular Complex Detection) plug-in was employed to screen the modules or clusters of the PPI network. The MCODE parameters were set by default, and the KEGG as well as GO assessments of DEGs in modules were utilized in their investigation using DAVID, in which was set at the cut-off criterion.

3. Results

3.1. Differentially Expressed lncRNAs and mRNAs

Before data analysis, we performed principal component analysis (PCA) for quality control (Supplementary figure (available here)); the PCA plots showed that the three datasets explored in this study exhibited distinguished lncRNA profiles between normal and tumor biopsies. Three HBV-linked HCC gene expression patterns were abstracted from the GEO data resource. Afterward, we normalized the gene expression data, and the DEGs were determined with the limma package of R. Using and as cut-off criteria, 103, 333, and 158 differentially expressed lncRNAs (DElncRNAs) and 1182, 2147, and 1579 differentially expressed mRNAs (DEmRNAs) were extracted from the expression profile datasets GSE55092, GSE19665, and GSE84402, respectively, presented in heat map plots (Figure 1). On the basis of the integration along with the analysis of RRA, 38 DElncRNAs, including 25 upregulated and 13 downregulated lncRNAs, as well as 541 DEmRNAs, including 348 upregulated and 195 downregulated mRNAs, of invariably expressed genes were uncovered from the three profile cohorts.

(a)

(b)

(c)

3.2. Optimal Diagnostic lncRNA Biosignatures for HBV-Related HCC

To assess the optimal diagnostic lncRNA biosignatures for HBV-related HCC, classification models were carried out. The random forest assessment was employed to rank DElncRNAs from large to small, according to the decrease in Mean Decrease Gini as indicated in Figure 2(a). The 10-fold cross-validation indicated that the average accuracy rate of 9 DElncRNAs reached their highest point for the first time as illustrated in Figure 2(b). These 9 DElncRNAs consisting of five upregulated DElncRNAs (AL356056.2, AL445524.1, TRIM52-AS1, AC093642.1, and EHMT2-AS1) and four downregulated DElncRNAs (AC003991.1, AC008040.1, LINC00844, and LINC01018) were identified as the optimal lncRNA molecular markers for the HBV-related HCC diagnosis (Table 1). Therefore, we selected these 9 DElncRNAs as promising optimal lncRNA molecular markers for the diagnosis of HBV-linked HCC, and they were employed to develop the random forest, as well as the SVM models.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

(j)

(k)

Figure 2

Determination of the prospective lncRNA biosignatures for HBV-related HCC. (a) The importance value of every DElncRNA was ranked as per the Mean Decrease Gini through the random decision forest. (b) The variance rate of classification performance when increasing numbers of the predictive DElncRNAs. Box plot exhibiting the expression levels of (c) EHMT2-AS1, (d) LINC01018, (e) AL356056.2, (f) AL445524.1, (g) LINC00844, (h) AC093642.1, (i) TRIM52-AS1, (j) AC003991.1, and (k) AC008040.1 between HBV-related HCC and nontumor tissues. lncRNA: long non-coding RNA; HCC: hepatocellular carcinoma.

3.3. Cross-Validation

A box plot was used to display the levels of expression of the nine DElncRNAs between HBV-linked HCC and nontumor tissues in the validation set (Figures 3(c)–3(g)). The AUC of the random forest model was 0.904, and the specificity and sensitivity were 94.3 and 86.5%, respectively (Figure 3(j)). Moreover, the AUC of the SVM model was 0.957, while the specificity and sensitivity were 95.7 and 100%, respectively (Figure 3(k)). The AUC of the logistic regression model was 0.976, while the specificity and sensitivity were 89.2 and 95.7%, respectively (Figure 3(l)). The ROC of these nine DElncRNAs, AL356056.2, AL445524.1, TRIM52-AS1, AC093642.1, EHMT2-AS1, AC003991.1, AC008040.1, LINC00844, and LINC01018, were 0.750, 0.817, 0.905, 0.778, 0.927, 0.859, 0.722, 0.776, and 0.719, respectively (Figures 3(a)–3(i)), which indicated that these nine lncRNAs were related to HBV-linked HCC, playing a core role in estimating the onset of HBV-linked HCC.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

(j)

(k)

(l)

Figure 3

ROC assessment of nine HBV-related HCC-distinct lncRNA biosignatures. The ROC data of these nine diagnostic lncRNA biosignatures (AL356056.2, AL445524.1, TRIM52-AS1, AC093642.1, EHMT2-AS1, AC003991.1, AC008040.1, LINC00844, and LINC01018) and individual (a) AL356056.2, (b) AL445524.1, (c) TRIM52-AS1, (d) AC093642.1, (e) EHMT2-AS1, (f) AC003991.1, (g) AC008040.1, (h) LINC00844, and (i) LINC01018; their combination based on (j) random forest, (k) support vector machine model, and (l) logistic regression. AUC: area under the curve; HCC: hepatocellular carcinoma; ROC: receiver operating characteristic; lncRNA: long noncoding RNA.

3.4. DEmRNA-DElncRNA Coexpressed Analysis

In order to assess the possible functions of these dysregulated lncRNAs in HBV-related HCC, dysregulated DEmRNAs and DElncRNAs were selected to develop a coexpression network. Consequently, the optimal diagnostic lncRNAs were coexpressed with 126 DEmRNAs, accounting for 199 DEmRNA-DElncRNA coexpression pairs. The details pertaining to regulatory information is shown in Figure 4.

3.5. Functional and Signal Pathway Enrichment Analysis

The above-mentioned coexpression of 126 DEmRNAs with optimal diagnostic lncRNAs was employed to conduct the GO, as well as the KEGG enrichment analyses. GO analysis could classify genes into 3 categories: biological process (BP), cellular component (CC), and molecular function (MF). indicated statistical significance, in which the 3 parts of the GO results of the dysregulated genes are illustrated in Figure 5. The top 15 data generated from the GO enrichment of the dysregulated genes are indicated in Table 2.

KEGG analysis of DEGs demonstrated that they were majorly associated within the p53 signaling cascade, retinol metabolism, PI3K-Akt signaling cascade, chemical carcinogenesis, peroxisome, and caffeine metabolism. The network diagram was then drawn using the Cytoscape software, as shown in Figure 5. A total of 5 DEmRNAs, including STEAP3, CCNE1, CDK1, RRM2, and IGFBP3, were found to be enriched in the p53 signaling pathway. Meanwhile, four DEmRNAs, including CYP3A5, CYP2C18, ADH4, and HSD17B6, were observed to be remarkably enriched in the retinol metabolism pathway. FGFR2, CCNE1, ITGA9, ITGA6, GYS2, EFNA4, GHR, and SPP1 were enriched in the PI3K-Akt signaling pathway, while four DEmRNAs (CYP3A5, CYP2C18, ADH4, and NAT2) were significantly enriched in the chemical carcinogenesis pathway. Furthermore, four DEmRNAs (XDH, HAO2, ACSL4, and SLC27A2) and two DEmRNAs (XDH and NAT2) were enriched in the peroxisome pathway and caffeine metabolism pathway, respectively.

3.6. Protein-Protein Interaction (PPI) Network and Module Analyses

Subsequently, a PPI network was built in order to predict the interactions among proteins encoded by the 126 DEGs. As a result, 87 DEGs with a (medium confidence) were selected for the development of the PPI networks (Figure 6(a)). Afterwards, we clustered two remarkable gene modules by the MCODE APP in Cytoscape, in which module 1 constituted 18 genes/nodes and 127 edges (Figure 6(b)), whereas module 2 comprised 6 genes/nodes and 9 edges (Figure 6(c)). Accordingly, module 1 was primarily abundant in the p53 signaling cascade, cell cycle, and DNA replication, whereas genes in module 2 were primarily abundant in Salmonella infection, NF-kappa B signaling cascade, Toll-like receptor signaling cascade, and phagosome and tuberculosis cascade (Table 3).

(a)

(b)

(c)

4. Discussion

Recently, the diagnosis and treatment of hepatocellular carcinoma have been significantly improved. Although HCC patients benefit from a variety of treatment options, the survival time of HCC patients is still limited, regardless of the cancer stage at diagnosis. Hence, in regard to the significance of early diagnosis on survival, further investigations toward the determination of accurate, as well as specific, molecular markers of HCC are urgently needed. A growing number of researches have reported that lncRNA plays a crucial role in hepatocellular carcinoma; lncRNA also has a certain role in HBV-linked HCC [14]. Nevertheless, research on lncRNAs as predictive molecular markers as well as treatment targets remains limited.

Herein, three gene expression datasets were analyzed, in which 38 DElncRNAs along with 543 DEmRNAs were identified in HBV-linked HCC and neighboring samples. The feature selection along with the classification model analysis were employed to identify nine optimal diagnostic lncRNA molecular markers consisting of AL356056.2, AL445524.1, TRIM52-AS1, AC093642.1, EHMT2-AS1, AC003991.1, AC008040.1, LINC00844, and LINC01018 for further analysis. Accordingly, the possible role of these uncovered genes and their correlation with the occurrence and progress of HBV-related HCC were studied.

GO analyses uncovered the prospective biological roles of the abnormally expressed protein-coding genes, while the GO function enrichment results demonstrated that the dysregulated genes were associated with the monocarboxylic acid metabolic and steroid metabolic processes, cell division, mitotic cell cycle process (BP), basal plasma membrane, spindle pole, basolateral plasma membrane, basal part of cell, extracellular space, spindle (CC), cofactor binding, protein homodimerization activity, coenzyme binding, small molecule binding, and identical protein binding (MF). Moreover, the results showed that DEGs were involved in proliferation and the cell cycle of HCC cells induced by chronic hepatitis B.

The pathway analysis demonstrated that the integrated DGEs were particularly abundant in the p53 signaling cascade, retinol metabolism, PI3K-Akt signaling axis, chemical carcinogenesis, and peroxisome. Previous researches have verified that the p53 signaling cascade has a pivotal role in modulating HCC tumorigenesis, mediating HCC inhibition, apoptosis, and senescence, and DNA damage [15]. The PI3K/Akt signaling cascade is overactive in numerous cancers, and activation of this pathway contributes to cell proliferation, survival, and motility, as well as angiogenesis [16, 17]. Previous investigations have demonstrated that the PI3K/Akt signaling cascade is involved in HCC occurrence and metastasis [18]. Cancer is definitely linked to abnormal metabolic processes, and altered retinol metabolism is a pathway that participates in the process of liver fibrosis, where enzymes related to retinol metabolism are linked to liver cancer [19]. Increasing data have suggested that peroxisomes have a pivotal role in cancer. Enzymes involved in peroxisome lipid processing have been previously found to be elevated in a variety of cancers, including liver cancer [20]. Peroxisomes participate in the onset and progress of tumors through autophagic degradation, regulation of cellular immune response, and regulation of metabolism by crosstalk with the mitochondria [21]. Chemical carcinogenesis is a process in which cells undergo genetic and epigenetic variation based on oncogene and tumor suppressors [22, 23]. To this effect, studies have documented that chemical carcinogenesis participates in the onset and progress of liver cancer [24–26]. Therefore, exploring these cascades will aid and elucidate the mechanism of occurrence, proliferation, and invasion of HBV-linked HCC and help predict tumor progression.

The ROC data illustrated the significant diagnostic value of these nine diagnostic lncRNA biosignatures and their combinations for HBV-related HCC. Upregulated EHMT2-AS1 was found to be in the core of the coexpression network, with an AUC of 0.927 and a specificity along with a sensitivity of 90.0 and 91.9%, respectively. Genome-wide association studies also demonstrated that the variants of EHMT2-AS1 were linked to a risk of chronic hepatitis B [27, 28]. It has been reported that upregulated TRIM52-AS1 served an oncogenic role in driving HCC progression, which may be employed as a novel treatment target for HCC [29]. Herein, TRIM52-AS1 was observed to be elevated in HBV-linked HCC tissues in contrast to the neighboring tissues; however, the diagnostic significance of HBV-related HCC was shown to have an AUC of 0.905, as well as a specificity and sensitivity of 80.0 and 94.6%, respectively.

Research has shown TRIM52-AS1 upregulation in both HCC tissue and cell lines. TRIM52-AS1 silencing suppresses cell growth, migration, and infiltration, as well as epithelial-mesenchymal transition in vitro while repressing tumor growth in vivo. In a certain mechanism, TRIM52-AS1 sponged miR-514a-5p in promoting HCC progression by increasing MRPS18A expression. In contrast, another study demonstrated TRIM52-AS1 downregulation in renal cell carcinoma (RCC). Here, overexpression of TRIM52-AS1 inhibited cell migration as well as proliferation, triggered apoptosis of RCC cells in vitro, and interfered with its expression, leading to the opposite effect [30]. However, the underlying biological function of TRIM52-AS1 in cellular processes remains unclear. Further investigation may provide novel insight into the development along with the progression of HBV-linked HCC.

LINC01018 was found to be downregulated in this study. It has been reported that LINC01018 and FOXO1 decreased in HCC, in which LINC01018 sponged miR-182-5p targeted with FOXO1. Overexpression of LINC01018 sponging miR-182-5p repressed proliferation and enhanced HCC cell apoptosis via upregulation of FOXO1 [31]. Moreover, studies have shown that LINC01018 plays critical roles in HBV-induced liver cancer via ceRNA networks [32].

Additionally, investigations have demonstrated that LINC00844 affects drug metabolism and toxicity by regulating the expression of the drug-metabolizing enzyme [33]. Meanwhile, LINC00844 was reported to be downregulated in HCC, and its overexpression remarkably suppressed the proliferation, migration, and invasion of HCC cells [34].

Other than EHMT2-AS1, TRIM52-AS1, LINC01018, and LINC00844, the other five lncRNAs have never previously been documented, and their biological roles remain unclear. The coexpression assessment of lncRNA-mRNA constitute the most frequent approach for determining the prospective target genes of lncRNA and further studying its biological role in many diseases. Specific lncRNA biomarkers of HBV-linked HCC were employed to develop the coexpression network of DElncRNA-DEmRNA, and the role of DEmRNA coexpressed with the nine lncRNAs was annotated.

Accordingly, DEmRNAs and STEAP3 coexpressed with AC003991.1; CCNE1 coexpressed with EHMT2-AS1; CDK1 coexpressed with LINC01018, EHMT2-AS1, and TRIM52-AS1; RRM2 coexpressed with AC093642.1 and TRIM52-AS1; and IGFBP3 coexpressed with AL445524.1 were found to be enriched in the p53 signaling pathway. Cancer cells require considerable amounts of iron due to rapid proliferation; STEAP3 encodes protein functions as an iron transporter, which may regulate intracellular iron storage and help tumors grow in iron-deficient environments [33]. Additionally, STEAP3 and STEAP4 can form a heterodimer, which can affect metal homeostasis and apoptosis, as well as cell cycle regulation [35]. Cyclin E has been implicated in various carcinomas; specifically, CCNE1 (Cyclin E1) is upregulated by HOXC13 to promote proliferation and cell cycle progression [36]. CCNE1 overexpression is linked to the initiation of liver cancer [37] and dismal prognosis in patients suffering from HCC [38]. Moreover, Cyclin-dependent kinase 1 (CDK1) is a vital player in cell cycle modulation, where elevated CDK1 is linked to dismal overall survival in individuals with HCC [38]. CDK1/PDK1/β-Catenin cascades are activated in HCC, which may function as a prospective treatment target pathway [39]. RRM2 is a ribonucleotide reductase, a rate-limiting enzyme for cell replication, which catalyzes the reduction of ribonucleotides into deoxyribonucleotides during DNA synthesis. Research has shown that RRM2 modulates cancer cell growth, differentiation, and metastasis, as well as drug resistance, which is overexpressed in prostate cancer, breast cancer, and cervical cancer [40–43]. The IGFBP3 (insulin-like growth factor-binding protein-3) can competitively bind IGF1 (insulin-like growth factor-1) with the IGF1 receptor, thereby modulating cell growth and differentiation, as well as apoptosis. Furthermore, IGFBP3 is decreased in HCC, and low expressions of IGFBP3 are remarkably linked to the tumor size, portal vein invasion, and histological differentiation, as well as capsular invasion [44]. Serum levels of IGF1 along with IGFBP3 are also positively correlated with HCC incidence [45]. IGFBP3 silencing in various tumors is because of the methylation of its promoter, while the high IGFBP3 expression triggers apoptosis and inhibits cell survival, as well as growth [46, 47]. Therefore, the present study puts forward that AC003991.1, EHMT2-AS1, LINC01018, TRIM52-AS1, AC093642.1, and AL445524.1 have pivotal roles in the onset or progress of HBV-related HCC by modulating the p53 signaling cascade.

Cytochrome P450 enzymes participate in multiple types of cancer through a variety of mechanisms, including metabolic enzymes, catalysis of biological activity of chemical carcinogens, activating cancer treatment drugs, and as targets for cancer therapy [48, 49]. CYP3A5, cytochrome P450 family 3 subfamily A member 5, encodes a member of the cytochrome P450 superfamily of enzymes that participate in the metabolic processes of endogenous molecules, such as drugs and steroids, as well as exogenous carcinogens [50]. CYP3A4 along with CYP3A5 account for about 30% of hepatic cytochrome P450, and about half of the drugs oxidatively metabolized by P450 are CYP3A substrates [51]. In recent years, however, numerous studies have documented that CYP3A5 is aberrantly expressed in multiple tumors, and its aberrant expression is strongly related to the infiltration as well as metastasis of tumors [52, 53]. Previous studies have also shown that CYP3A5 suppresses HCC pathogenesis and metastasis through the modulation of the mTORC2/Akt signaling cascade, serving as a prognostic biosignature [54, 55]. Aflatoxin is one of the leading predisposing factors for HCC, and CYP3A5 polymorphism is linked to elevated mutagenic AFB1-exo-8,9-epoxide levels, which may control the individual risk of HCC [56]. In accordance with these findings, the data of the present study demonstrated that CYP3A5 was reduced in HCC. CYP3A5, which coexpressed with LINC01018, EHMT2-AS1, and AC008443.1, was observed to be enriched in both the retinol metabolism and chemical carcinogenesis signaling pathways. Therefore, it was hypothesized that LINC01018, EHMT2-AS1, and AC008443.1 may serve important roles in HBV-related HCC by modulating the expression of CYP3A5 or retinol metabolism and chemical carcinogenesis.

A KEGG pathway evaluation of genes was also conducted in the two most important modules according to the module analysis of the generated PPI network. The results demonstrated that the genes in module 1 were primarily abundant in the p53 signaling cascade, cell cycle, and DNA replication, whereas the genes in module 2 were primarily abundant in the Salmonella infection, NF-kappa B signaling axis, Toll-like receptor signaling cascade, and phagosome, as well as tuberculosis pathways. These data suggested that the onset along with the development of HBV-linked HCC may be correlated to these cascades, and the regulation of cell cycle-related biological processes or pathways may be an effective treatment for HBV-related HCC.

5. Conclusion

Here, a 9-lncRNA diagnosis model was constructed based on microarray datasets to highlight its diagnostic significance for HBV-linked HCC. In this regard, the functional annotation of DEmRNAs coexpression with nine lncRNA molecular markers of HBV-related HCC yielded valuable insights in assessing the pathogenesis of HBV-related HCC. Collectively, these data may provide an understanding of the regulatory mechanisms of lncRNAs in the onset, as well as the development of HBV-related HCC, and help identify potential diagnostic and treatment targets for HBV-linked HCC patients. Nevertheless, this work had some limitations. The sample size was small, and the samples were acquired from an independent patient cohort. Furthermore, this study only preliminarily predicts the role of lncRNAs in HBV-linked HCC, and further wet experiments are required in order to verify their mechanisms.

Data Availability

The datasets, as well as the code generated or analyzed in this work, can be accessed from the corresponding authors on reasonable request.

Conflicts of Interest

The authors have no conflicting interests to disclose.

Authors’ Contributions

Shunqiang Nong, Xiaohao Chen, and Zechen Wang contributed equally to this work.

Acknowledgments

This research was funded by grants from the Guangxi Natural Science Foundation (Nos. 2018JJA14004 and 2020GXNSFAA259050), Guangxi Clinic Medicine Research Center of Hepatobiliary Disease (No. AD17129025), Innovation Program of Guangxi Graduate Education (No. YCSW2019222), and Youjiang Medical University for Nationalities Research Project (No. yy2019bsky001).

Supplementary Materials

PCA analysis of quality control over three datasets used in this study. (Supplementary Materials)

References

F. Bray, J. Ferlay, I. Soerjomataram, R. L. Siegel, L. A. Torre, and A. Jemal, “Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries,” CA: a Cancer Journal for Clinicians, vol. 68, no. 6, pp. 394–424, 2018.
View at: Publisher Site | Google Scholar
Y. C. Liu, L. F. Lu, C. J. Li et al., “Hepatitis B virus X protein induces RHAMM-dependent motility in hepatocellular carcinoma cells via PI3K-Akt-Oct-1 signaling,” Molecular Cancer Research, vol. 18, no. 3, pp. 375–389, 2020.
View at: Publisher Site | Google Scholar
S. Zheng, H. Wu, F. Wang et al., “The oncoprotein HBXIP facilitates metastasis of hepatocellular carcinoma cells by activation of MMP15 expression,” Cancer Management and Research, vol. 11, pp. 4529–4540, 2019.
View at: Publisher Site | Google Scholar
Y. Zha, Q. Yao, J. S. Liu, Y. Y. Wang, and W. M. Sun, “Hepatitis B virus X protein promotes epithelial-mesenchymal transition and metastasis in hepatocellular carcinoma cell line HCCLM3 by targeting HMGA2,” Oncology Letters, vol. 16, pp. 5709–5714, 2018.
View at: Publisher Site | Google Scholar
O. M. Ali, H. A. el Amin, Y. L. Sharkawy, A. A. Mohamed Ali, E. Kholef, and W. Elsewify, “Golgi protein 73 versus alpha-fetoprotein as a new biomarker in early diagnosis of hepatocellular carcinoma,” International Journal of General Medicine, vol. 13, pp. 193–200, 2020.
View at: Publisher Site | Google Scholar
H. J. Cho, G. O. Baek, C. W. Seo et al., “Exosomal microRNA-4661-5p-based serum panel as a potential diagnostic biomarker for early-stage hepatocellular carcinoma,” Cancer Medicine, vol. 9, no. 15, pp. 5459–5472, 2020.
View at: Publisher Site | Google Scholar
X. L. Ma, W. G. Tang, M. J. Yang et al., “Serum STIP1, a novel indicator for microvascular invasion, Predicts Outcomes and Treatment Response in Hepatocellular Carcinoma,” Frontiers in Oncology, vol. 10, p. 511, 2020.
View at: Publisher Site | Google Scholar
M. K. Iyer, Y. S. Niknafs, R. Malik et al., “The landscape of long noncoding RNAs in the human transcriptome,” Nature Genetics, vol. 47, no. 3, pp. 199–208, 2015.
View at: Publisher Site | Google Scholar
X. Zhang, W. Wang, W. Zhu et al., “Mechanisms and functions of long non-coding RNAs at multiple regulatory levels,” International Journal of Molecular Sciences, vol. 20, no. 22, p. 5573, 2019.
View at: Publisher Site | Google Scholar
R. R. Pandey and C. Kanduri, “Transcriptional and posttranscriptional programming by long noncoding RNAs,” Progress in Molecular and Subcellular Biology, vol. 51, pp. 1–27, 2011.
View at: Publisher Site | Google Scholar
C. Klec, T. Gutschner, K. Panzitt, and M. Pichler, “Involvement of long non-coding RNA HULC (highly up-regulated in liver cancer) in pathogenesis and implications for therapeutic intervention,” Expert Opinion on Therapeutic Targets, vol. 23, no. 3, pp. 177–186, 2019.
View at: Publisher Site | Google Scholar
Z. Huang, J. K. Zhou, Y. Peng, W. He, and C. Huang, “The role of long noncoding RNAs in hepatocellular carcinoma,” Molecular Cancer, vol. 19, no. 1, p. 77, 2020.
View at: Publisher Site | Google Scholar
T. Lan, L. Chang, L. Wu, and Y. Yuan, “Downregulation of ZEB2-AS1 decreased tumor growth and metastasis in hepatocellular carcinoma,” Molecular Medicine Reports, vol. 14, no. 5, pp. 4606–4612, 2016.
View at: Publisher Site | Google Scholar
R. Chauhan and N. Lahiri, “Tissue- and serum-associated biomarkers of hepatocellular carcinoma,” Biomark Cancer, vol. 8s1, pp. BIC.S34413–BIC.S34455, 2016.
View at: Publisher Site | Google Scholar
T. Lin, P. F. Hou, S. Meng et al., “Emerging roles of p53 related lncRNAs in cancer progression: a systematic review,” International Journal of Biological Sciences, vol. 15, no. 6, pp. 1287–1298, 2019.
View at: Publisher Site | Google Scholar
P. Liu, H. Cheng, T. M. Roberts, and J. J. Zhao, “Targeting the phosphoinositide 3-kinase pathway in cancer,” Nature Reviews. Drug Discovery, vol. 8, no. 8, pp. 627–644, 2009.
View at: Publisher Site | Google Scholar
M. de Santis, V. Sala, M. Martini, G. B. Ferrero, and E. Hirsch, “PI3K Signaling in tissue hyper-proliferation: from overgrowth syndromes to kidney cysts,” Cancers, vol. 9, no. 12, p. 30, 2017.
View at: Publisher Site | Google Scholar
D. L. Ou, B. S. Lee, L. I. Lin et al., “Vertical blockade of the IGFR-PI3K/Akt/mTOR pathway for the treatment of hepatocellular carcinoma: the role of survivin,” Molecular Cancer, vol. 13, no. 1, p. 2, 2014.
View at: Publisher Site | Google Scholar
P. Pettinelli, B. M. Arendt, A. Teterina et al., “Altered hepatic genes related to retinol metabolism and plasma retinol in patients with non-alcoholic fatty liver disease,” PLoS One, vol. 13, no. 10, article e0205747, 2018.
View at: Publisher Site | Google Scholar
M. Cai, X. Sun, W. Wang et al., “Disruption of peroxisome function leads to metabolic stress, mTOR inhibition, and lethality in liver cancer cells,” Cancer Letters, vol. 421, pp. 82–93, 2018.
View at: Publisher Site | Google Scholar
M. S. Dahabieh, E. Di Pietro, M. Jangal et al., “Peroxisomes and cancer: the role of a metabolic specialist in a disease of aberrant metabolism,” Biochimica Et Biophysica Acta. Reviews on Cancer, vol. 1870, pp. 103–121, 2018.
View at: Publisher Site | Google Scholar
P. A. Oliveira, A. Colaço, R. Chaves, H. Guedes-Pinto, P. L. F. de-la-Cruz, and C. Lopes, “Chemical carcinogenesis,” Anais da Academia Brasileira de Ciências, vol. 79, no. 4, pp. 593–616, 2007.
View at: Publisher Site | Google Scholar
L. A. Loeb and C. C. Harris, “Advances in chemical carcinogenesis: a historical review and prospective,” Cancer Research, vol. 68, pp. 6863–6872, 2008.
View at: Publisher Site | Google Scholar
G. He, D. Dhar, H. Nakagawa et al., “Identification of liver cancer progenitors whose malignant progression depends on autocrine IL-6 signaling,” Cell, vol. 155, no. 2, pp. 384–396, 2013.
View at: Publisher Site | Google Scholar
J. Han, L. Ding, B. Yuan et al., “Hepatitis B virus X protein and the estrogen receptor variant lacking exon 5 inhibit estrogen receptor signaling in hepatoma cells,” Nucleic Acids Research, vol. 34, pp. 3095–3106, 2006.
View at: Publisher Site | Google Scholar
T. Uehara, G. R. Ainslie, K. Kutanzi et al., “Molecular mechanisms of fibrosis-associated promotion of liver carcinogenesis,” Toxicological Sciences, vol. 132, no. 1, pp. 53–63, 2013.
View at: Publisher Site | Google Scholar
Y. J. Kim, H. Y. Kim, J. H. Lee et al., “A genome-wide association study identified new variants associated with the risk of chronic hepatitis B,” Human Molecular Genetics, vol. 22, no. 20, pp. 4233–4238, 2013.
View at: Publisher Site | Google Scholar
D. K. Jiang, X. P. Ma, H. Yu et al., “Genetic variants in five novel loci including CFB and CD40 predispose to chronic hepatitis B,” Hepatology, vol. 62, no. 1, pp. 118–128, 2015.
View at: Publisher Site | Google Scholar
C. Zhou, Z. Chen, C. Peng, C. Chen, and H. Li, “Long noncoding RNA TRIM52-AS1 sponges miR-514a-5p to facilitate hepatocellular carcinoma progression through increasing MRPS18A,” Cancer Biotherapy & Radiopharmaceuticals, vol. 36, no. 2, pp. 211–219, 2021.
View at: Publisher Site | Google Scholar
Z. Liu, H. Y. Yan, S. Y. Xia, C. Zhang, and Y. C. Xiu, “Downregulation of long non-coding RNA TRIM52-AS1 functions as a tumor suppressor in renal cell carcinoma,” Molecular Medicine Reports, vol. 13, pp. 3206–3212, 2016.
View at: Publisher Site | Google Scholar
S. Wang, M. Xu, Z. Sun, X. Yu, Y. Deng, and H. Chang, “LINC01018 confers a novel tumor suppressor role in hepatocellular carcinoma through sponging microRNA-182-5p,” American Journal of Physiology. Gastrointestinal and Liver Physiology, vol. 317, no. 2, pp. G116–116G126, 2019.
View at: Publisher Site | Google Scholar
Q. Wen, D. Wang, Y. Yang et al., “Competing endogenous RNA screening based on long noncoding RNA-messenger RNA co-expression profile in hepatitis B virus-associated hepatocarcinogenesis,” Journal of Traditional Chinese Medicine, vol. 37, no. 4, pp. 510–521, 2017.
View at: Google Scholar
T. Isobe, E. Baba, S. Arita et al., “Human STEAP3 maintains tumor growth under hypoferric condition,” Experimental Cell Research, vol. 317, no. 18, pp. 2582–2591, 2011.
View at: Publisher Site | Google Scholar
W. Zhou, K. Huang, Q. Zhang et al., “LINC00844 promotes proliferation and migration of hepatocellular carcinoma by regulating NDRG1 expression,” PeerJ, vol. 8, article e8394, 2020.
View at: Publisher Site | Google Scholar
R. T. Scarl, C. M. Lawrence, H. M. Gordon, and C. S. Nunemaker, “STEAP4: its emerging role in metabolism and homeostasis of cellular iron and copper,” The Journal of Endocrinology, vol. 234, no. 3, pp. R123–123R134, 2017.
View at: Publisher Site | Google Scholar
Y. Yao, J. Luo, Q. Sun et al., “HOXC13 promotes proliferation of lung adenocarcinoma via modulation of CCND1 and CCNE1,” American Journal of Cancer Research, vol. 7, no. 9, pp. 1820–1834, 2017.
View at: Google Scholar
R. Sonntag, N. Giebeler, Y. A. Nevzorova et al., “Cyclin E1 and cyclin-dependent kinase 2 are critical for initiation, but not for progression of hepatocellular carcinoma,” Proceedings. National Academy of Sciences. United States of America, vol. 115, no. 37, pp. 9282–9287, 2018.
View at: Publisher Site | Google Scholar
X. Bisteau, M. J. Caldez, and P. Kaldis, “The complex relationship between liver cancer and the cell cycle: a story of multiple regulations,” Cancers, vol. 6, pp. 79–111, 2014.
View at: Publisher Site | Google Scholar
C. X. Wu, X. Q. Wang, S. H. Chok et al., “Blocking CDK1/PDK1/β-catenin signaling by CDK1 inhibitor RO3306 increased the efficacy of sorafenib treatment by targeting cancer stem cells in a preclinical model of hepatocellular carcinoma,” Theranostics, vol. 8, no. 14, pp. 3737–3750, 2018.
View at: Publisher Site | Google Scholar
B. L. Gan, L. J. Zhang, L. Gao et al., “Downregulation of miR-224-5p in prostate cancer and its relevant molecular mechanism via TCGA, GEO database and in silico analyses,” Oncology Reports, vol. 40, pp. 3171–3188, 2018.
View at: Publisher Site | Google Scholar
W. H. Liang, N. Li, Z. Q. Yuan, X. L. Qian, and Z. H. Wang, “DSCAM-AS1 promotes tumor growth of breast cancer by reducing miR-204-5p and up-regulating RRM2,” Molecular Carcinogenesis, vol. 58, pp. 461–473, 2019.
View at: Publisher Site | Google Scholar
N. Wang, T. Zhan, T. Ke et al., “Increased expression of RRM2 by human papillomavirus E7 oncoprotein promotes angiogenesis in cervical cancer,” British Journal of Cancer, vol. 110, pp. 1034–1044, 2014.
View at: Publisher Site | Google Scholar
H. Feng, Q. Wang, W. Xiao, B. Zhang, Y. Jin, and H. Lu, “lncRNA TTN-AS1 regulates miR-524-5p and RRM2 to promote breast cancer progression,” Oncotargets and Therapy, vol. 13, pp. 4799–4811, 2020.
View at: Publisher Site | Google Scholar
S. Aishima, Y. Basaki, Y. Oda et al., “High expression of insulin-like growth factor binding protein-3 is correlated with lower portal invasion and better prognosis in human hepatocellular carcinoma,” Cancer Science, vol. 97, no. 11, pp. 1182–1190, 2006.
View at: Publisher Site | Google Scholar
J. Wang, Y. C. Li, M. Deng et al., “Serum insulin-like growth factor-1 and its binding protein 3 as prognostic factors for the incidence, progression, and outcome of hepatocellular carcinoma: a systematic review and meta-analysis,” Oncotarget, vol. 8, pp. 81098–81108, 2017.
View at: Publisher Site | Google Scholar
A. A. Dar, S. Majid, M. Nosrati, D. de Semir, S. Federman, and M. Kashani-Sabet, “Functional modulation of IGF-binding protein-3 expression in melanoma,” The Journal of Investigative Dermatology, vol. 130, pp. 2071–2079, 2010.
View at: Publisher Site | Google Scholar
Y. Hashimoto, T. J. Zumwalt, and A. Goel, “DNA methylation patterns as noninvasive biomarkers and targets of epigenetic therapies in colorectal cancer,” Epigenomics, vol. 8, pp. 685–703, 2016.
View at: Publisher Site | Google Scholar
A. M. Alzahrani and P. Rajendran, “The multifarious link between cytochrome P450s and cancer,” Oxidative Medicine and Cellular Longevity, vol. 2020, Article ID 3028387, 18 pages, 2020.
View at: Publisher Site | Google Scholar
E. G. Hrycay and S. M. Bandiera, “Involvement of cytochrome P450 in reactive oxygen species formation and cancer,” Advances in Pharmacology, vol. 74, pp. 35–84, 2015.
View at: Publisher Site | Google Scholar
I. Fleming, “The cytochrome P450 pathway in angiogenesis and endothelial cell biology,” Cancer Metastasis Reviews, vol. 30, pp. 541–555, 2011.
View at: Publisher Site | Google Scholar
U. M. Zanger, M. Turpeinen, K. Klein, and M. Schwab, “Functional pharmacogenetics/genomics of human cytochromes P450 involved in drug biotransformation,” Analytical and Bioanalytical Chemistry, vol. 392, pp. 1093–1108, 2008.
View at: Publisher Site | Google Scholar
M. S. Islam, A. G. Mostofa, M. U. Ahmed, M. S. Bin Sayeed, M. R. Hassan, and A. Hasnat, “Association of CYP3A4, CYP3A5 polymorphisms with lung cancer risk in Bangladeshi population,” Tumour Biology, vol. 35, pp. 1671–1678, 2014.
View at: Publisher Site | Google Scholar
L. Wang, W. Chen, J. Zha et al., “miR-543 acts as a novel oncogene in oral squamous cell carcinoma by targeting CYP3A5,” Oncology Reports, vol. 42, pp. 973–990, 2019.
View at: Publisher Site | Google Scholar
F. Jiang, L. Chen, Y. C. Yang et al., “CYP3A5 functions as a tumor suppressor in hepatocellular carcinoma by regulating mTORC2/Akt signaling,” Cancer Research, vol. 75, no. 7, pp. 1470–1481, 2015.
View at: Publisher Site | Google Scholar
T. Yu, X. Wang, G. Zhu et al., “The prognostic value of differentially expressed CYP3A subfamily members for hepatocellular carcinoma,” Cancer Management and Research, vol. 10, pp. 1713–1726, 2018.
View at: Publisher Site | Google Scholar
L. Wojnowski, P. C. Turner, B. Pedersen et al., “Increased levels of aflatoxin-albumin adducts are associated with CYP3A5 polymorphisms in the Gambia, West Africa,” Pharmacogenetics, vol. 14, pp. 691–700, 2004.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2021 Shunqiang Nong 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

407

Downloads

1077

Citations