Abstract

The role of long noncoding RNAs- (lncRNAs-) associated competing endogenous RNA (ceRNA) in the field of hepatocellular carcinoma (HCC) biology is well established, but the involvement of lncRNAs competing interactions in the progression of liver cirrhosis to HCC is still unclear. We aimed to explore the differential expression profiles of lncRNAs, microRNAs (miRNA), and messenger RNAs (mRNAs) to construct a functional ceRNA network in cirrhotic HCC. The lncRNA, miRNA, and mRNA expression datasets were obtained from Gene Expression Omnibus and The Cancer Genome Atlas. Based on miRanda and TargetScan, the HCC-specific ceRNA network was constructed to illustrate the coexpression regulatory relationship of lncRNAs, miRNAs, and mRNAs. The potential prognostic indicators in the network were confirmed by survival analysis and validated by qRT-PCR. A total of 74 lncRNAs, 36 intersection miRNAs, and 949 mRNAs were differentially expressed in cirrhotic HCC samples compared with cirrhosis samples. We constructed a ceRNA network, including 47 lncRNAs, 35 miRNAs, and 168 mRNAs. Survival analysis demonstrated that 2 lncRNAs (EGOT and SERHL), 4 miRNAs, and 40 mRNAs were significantly associated with the overall survival of HCC patients. Two novel regulatory pathways, EGOT-miR-32-5p-XYLT2 axis and SERHL-miR-1269a/miR-193b-3p-BCL2L1/SYK/ARNT/CHST3/LPCAT1 axis, were built up and contribute to the underlying mechanism of HCC pathogenesis. The higher-expressed SERHL was associated with a higher risk of all-cause death. The expressions of SERHL-miR-1269a-BCL2L1 were significantly different using qRT-PCR in vitro studies. lncRNAs EGOT and SERHL might serve as effective prognostic biomarkers and potential therapeutic targets in cirrhotic HCC treatment.

1. Introduction

Liver cancer is the sixth most common incident carcinoma and the fourth most common cause of malignant tumor mortality [1], being estimated to be 18.1 million new patients and 9.6 million cancer-related deaths throughout the world in 2018 [2]. Accounting for over 85% of primary liver cancers, hepatocellular carcinoma (HCC) is not evenly distributed among geographic regions, being accompanied by the increasing overall disease burden of HCC [3]. According to the etiology of cirrhosis, the risk factors for HCC mainly included hepatitis B virus infection, hepatitis C virus infection, alcohol, primary biliary cholangitis, autoimmune hepatitis, primary sclerosing cholangitis, and cryptogenic cirrhosis [4]. Therefore, liver cirrhosis plays a vital role in the pathogenesis of HCC. Based on the cirrhosis affecting progression to HCC, specific molecular biomarkers and potential regulatory mechanisms are essential and meaningful to the early diagnosis, treatment strategies, and the evaluation of prognosis.

Long noncoding RNAs (lncRNAs) are endogenous noncoding RNA molecules more than 200 nucleotides in length and can be subclassified into exonic, intronic, overlapping, and intergenic lncRNAs in nuclear or cytosolic fractions [5]. Given their regulation of cell proliferation, apoptosis, autophagy, epithelial-mesenchymal transition, and angiogenesis, cancer-specific lncRNAs are involved in the initiation, aggravation, and metastasis of HCC [6] and can serve as detectable diagnostic markers and effective therapeutic targets for HCC [7]. In the competing endogenous RNA (ceRNA) hypothesis, lncRNAs could act as competitive binding sites for the target microRNA (miRNA) response elements and then regulate the expression of corresponding messenger RNAs (mRNAs) [8]. The HCC-related deregulated ceRNA network was built to reveal the candidate prognostic cytoplasmic lncRNAs by integrated analysis in 371 HCC tissues and 50 healthy tissues from The Cancer Genome Atlas (TCGA) database [9]. However, the comprehensive analysis of the lncRNA-miRNA-mRNA ceRNA regulatory network in the progression of liver cirrhosis to HCC remains poorly characterized and be worth being further researched.

Given this context, we compared the differential expression of lncRNAs, miRNAs, and mRNAs between cirrhosis tissues and cirrhotic HCC tissues in the National Center for Biotechnology Information Gene Expression Omnibus (GEO) database and constructed the lncRNA-miRNA-mRNA ceRNA coexpression network. Finally, two lncRNAs were significantly correlated with the overall survival of cirrhotic HCC by the univariate Cox regression analysis, and their pair subnetworks were related to HCC-related signaling pathways. The previous correlational studies obtained the ceRNA expression in HCC samples and adjacent normal liver samples and explored the etiopathogenesis of normal liver to HCC. Hence, an integrated analysis of the ceRNA network to elaborate the biological processes and pathways on cancer-specific lncRNAs in cirrhotic HCC is still lacking. Our present study enrolled cirrhosis with and without HCC samples to focus on the pathologic process of cirrhosis to HCC, illuminating the underlying HCC pathogenesis of ceRNA coexpression network in patients with cirrhosis.

2. Materials and Methods

2.1. Microarray Datasets and Preprocessing

The lncRNAs, miRNAs, and mRNAs expression datasets were retrieved and downloaded from the GEO database by searching liver cirrhosis with and without HCC. GSE17967 (63 patients, 16 cirrhotic tissues with HCC, and 47 cirrhotic tissues without HCC, http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE17967) was included in our study and be used to analyze the differentially expressed lncRNAs and mRNAs [10]. Annotation pipeline of lncRNA transcripts was developed and obtained with Affymetrix Human Genome U133A 2.0 microarray. GSE21362 (73 patients, http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE21362) and GSE63046 (15 patients, http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE63046) were enrolled to profile and overlap the differentially expressed intersection miRNAs [11, 12]. The characteristics of each dataset were shown in Table S1. The expressions of specific lncRNAs, miRNAs, and mRNAs in cirrhotic HCC patients were downloaded from the National Institutes of Health TCGA Data Portal. TCGA-LIHC (377 patients, https://portal.gdc.cancer.gov/projects/TCGA-LIHC) was obtained for survival analysis. Differentially expressed lncRNAs, miRNAs, and mRNAs were extracted with restricted criteria, which were set at and . Figure S1 depicts the flowchart for bioinformatics analysis, pointing out a cluster of lncRNA, miRNA, and mRNA.

2.2. Construction of the ceRNA Network

The regulatory network of lncRNAs, miRNAs, and mRNAs was constructed to demonstrate the regulatory relation. The ceRNA hypothesis is that lncRNAs, acting as the specific miRNA sponges, can bind to target miRNAs and then regulate corresponding mRNAs expression. We performed the lncRNAs-mRNAs network to construct the coexpression regulatory relation by estimating across the array (Fig. S2). The miRanda (http://www.microrna.org/) and TargetScan (http://www.targetscan.org/) were performed to predict the target miRNA, which play a negative role in regulating the expression of lncRNA and mRNA. Next, we overlapped the intersection datasets both in miRanda and TargetScan. Finally, we combined the lncRNA-miRNA and miRNA-mRNA pairs to build the ceRNA network, which was visualized by using Cytoscape v2.8.2 [13].

2.3. Survival Analysis

To investigate the prognostic significance of differentially expressed lncRNAs, miRNAs, and mRNAs, we performed the survival analysis by collecting the ceRNA expression and clinical information in 377 cirrhotic HCC samples from TCGA-LIHC database. We used the survival data to evaluate the association between overall survival in cirrhotic HCC and candidate corresponding RNA expression, which split into high-expression and low-expression groups using the median. Kaplan-Meier survival analysis with log-rank test was applied and compared to estimate the significant difference by a cutoff of . The statistically significant lncRNAs in the survival analysis were used to reconstruct a new ceRNA network by Cytoscape v2.8.2. Baseline characteristics for lncRNAs EGOT and SERHL were depicted in Table S2. To clarify the associations with overall survival and lncRNAs EGOT and SERHL, we analyzed the clinical characteristics of enrolled patients, including age, sex, race, pathologic stage, Child-Pugh score, neoplasm histologic grade, Ishak fibrosis score, vascular invasion, residual tumor grade, radiation therapy, neoadjuvant therapy, embolization therapy, family history of cancer, risk factors, and biomarkers (AFP and albumin) (Table S3).

2.4. Cell Culture and Quantitative Real-Time Polymerase Chain Reaction (qRT-PCR)

We obtained the Huh7 human hepatoma cell line and LX-2 human hepatic fibrosis cell line from the Shanghai Institute of Biochemistry and Cell Biology (Shanghai, China). Huh7 cells were cultured with the Dulbecco’s Modified Eagle’s Medium (Thermo Fisher Scientific, Waltham, MA) containing 10% fetal bovine serum (FBS), penicillin (100 U/mL), L-glutamine (2 mM), and streptomycin (100 μg/mL) at 37°C in a humidified incubator (5% CO2). LX-2 cells were cultured in Medium 199 (Sigma-Aldrich, St. Louis, Missouri, USA) containing 10% FBS in a 5% CO2 atmosphere at 37°C. All the cell lines were seeded in microplates and maintained for 14 days. Total RNA was collected by a miRNeasy Mini Kit (Qiagen, Hilden, Germany) under the manufacturer’s instructions. Then, the quantity, integrity, and purity were assessed using a NanoDrop 2000 (Thermo Fisher Scientific, MA, USA). 50 ng of RNA was reverse-transcribed to cDNA by an iScript cDNA Synthesis Kit (170-8891, Bio-RAD). Quantitative PCR was performed using the 7900HT Rapid Real-Time Biosystems using SYBR-GREEN (170-8882AP, BIO-RAD) technology. Primers are listed in Table S4.

2.5. Statistical Analyses

Categorical variables were expressed as frequencies and proportions and compared using chi-square, Pearson’s chi-squared, or Fisher exact test. Continuous variables were expressed as means with standard deviation (SD) or medians with interquartile ranges (IQRs) and compared using Student’s -tests or Wilcoxon rank-sum test. The absolute values of and were considered to indicate significant difference. The univariate and multivariate Cox proportional hazard regression models were conducted for discovering the connections between ceRNAs and overall survival. Backward stepwise selection was applied by using Akaike’s information criterion to identify variables for multivariable Cox proportional hazard regression analysis. Then, we constructed a nomogram based on the multivariable Cox regression analyses to visualize the probability of 5-year survival in HCC patients. Variance inflation factor (VIF) assessed the multicollinearity of independent variables in the Cox regression model ( suggests multicollinearity). Hazard ratios (HRs) were presented with their 95% confidence intervals (CIs). All analyses were implemented using the SPSS software (version 21.0, IBM, Chicago, Illinois, USA) and the R software (version 3.6.2, R Development Core Team, Vienna, Austria).

3. Results

3.1. Differential Expression of lncRNAs, miRNAs, and mRNAs in Cirrhotic HCC

With the tumorigenesis of HCC in cirrhotic patients, a total of 74 lncRNAs from GSE17967, 36 intersection miRNAs from GSE21362 and GSE63046, and 949 mRNAs from GSE17967 were identified with the restricted criteria of and by bioinformatics analysis (Figure 1). The 74 HCC-specific differentially expressed lncRNAs were composed of 22 upregulated and 52 downregulated lncRNAs in GSE17967 (Figure 1(a); Table 1; Table S5). For 949 differentially expressed mRNAs in the GSE17967 dataset, there were 291 upregulated and 658 downregulated HCC-specific mRNAs (Figure 1(b); Table S6). We initially identified 121 miRNAs from GSE63046 (Figure 1(c)) and 106 miRNAs from the GSE21362 data set (Figure 1(d)), respectively. The intersections of the two datasets of differentially expressed miRNAs consisted of 36 (18.8%) miRNAs (5 upregulated and 31 downregulated, Figure 1(e); Table S7). These data provided a new clue to reveal the essential genes involved in the sequential progression from liver cirrhosis to HCC.

3.2. Construction of ceRNA Network

The lncRNA-mRNA coexpression network was built by differentially expressed lncRNAs and mRNAs to investigate the critical regulatory relations. A total of 16 lncRNAs (6 upregulated lncRNAs; 10 downregulated lncRNAs) and 219 mRNAs (96 upregulated mRNAs; 123 downregulated mRNAs) were enrolled in this network for further analysis. Next, based on miRanda and TargetScan analysis, the intersection differentially expressed miRNAs were used to predict the miRNAs negative targeted mRNAs, identifying 35 miRNAs negatively correlated to 360 mRNAs. Subsequently, we confirmed 36 miRNAs and the corresponding negative regulation to 69 lncRNAs. The lncRNA-miRNA-mRNA ceRNA coexpression network was built up based on the targeted interactions. We identified that 47 lncRNAs had regulatory effects on 35 intersection miRNAs (Table 2). Next, 35 miRNAs were interacting with 168 specific mRNAs (Table 3). Finally, the numbers of lncRNAs, miRNAs, and mRNAs were decreased to 47, 35, and 168, respectively (Figure 2).

3.3. Prognosis Prediction for Differentially Expressed lncRNAs, miRNAs, and mRNAs

To identify the HCC-specific ceRNAs with prognostic characteristics, we enrolled 377 cirrhotic HCC patients from the TCGA-LIHC database, and the univariate Cox regression model was implemented based on the expression of these lncRNAs, miRNAs, and mRNAs, which illustrated the relationship between key ceRNAs and overall survival of HCC patients. Our results demonstrated that 2 lncRNAs, 4 miRNAs, and 40 mRNAs were significantly correlated to the overall survival in patients with cirrhotic HCC (Fig. S3). The expressions of top two lncRNAs (EGOT and SERHL), top two miRNAs (hsa-miR-139-5p and hsa-miR-139-3p), and top two mRNAs (LPCAT1 and CASP2) were depicted in Figure 3 (). And Kaplan-Meier survival analyses revealed that these two lncRNAs, EGOT and SERHL, had negative correlation with overall survival in the progression of HCC. According to pathway analysis, we further demonstrated the expression pattern of EGOT and SERHL lncRNAs as well as their pair subnetworks (EGOT-miR-32-5p-XYLT2 axis and SERHL-miR-1269a/miR-193b-3p-BCL2L1/SYK/ARNT/CHST3/LPCAT1 axis) (Fig. S4).

3.4. The Associations between Key lncRNAs and Clinical Characteristics

Multivariate Cox analysis was implemented to investigate the predictors of a clinical prognosis of cirrhotic HCC (Table S8). Due to less than 10 points of all VIF in the multivariate Cox model, there was no evidence of multicollinearity in independent variables. Compared with stage I, a later pathologic stage of stages III and IV was associated with a lower overall survival rate (HR 1.59, 95% CI 1.05-2.43, and and HR 4.80, 95% CI 1.49-15.45, and , respectively), while there was no significant different in the pathologic stage of stage II (HR 0.93, 95% CI 0.58-1.50, ). In patients with cirrhotic HCC, a higher alpha-fetoprotein (AFP) (HR 1.73, 95% CI 1.18-2.55, ) and a higher albumin (HR 2.45, 95% CI 1.61-3.72, ) were related to a lower overall survival rate. In Cox hazard modeling for survival free from total mortality, the unadjusted HR for high-expressed EGOT group was 1.48 (95% CI 1.05-2.09, ), which decreased to 1.35 (95% CI 0.95-1.92, ) in a multivariable model that included pathologic stage, AFP, and albumin (Table S3). The unadjusted HR for high-expressed SERHL group was 1.97 (95% CI 1.39-2.81, ), which decreased to 1.74 (95% CI 1.22-2.49, ) in the adjusted multivariable model. A nomogram was then constructed by using AFP, serum albumin, pathologic stage, and lncRNA EGOT (Fig. S5A) or SERHL (Fig. S5B) to visualize the prediction model of 5-year survival probability.

3.5. qRT-PCR Verification

lncRNA EGOT-miR-32-5p-XYLT2 axis and lncRNA SERHL-miR-1269a/miR-193b-3p-BCL2L1/SYK/ARNT/CHST3/LPCAT1 axis were selected to establish the validity and reliability of our results. The expressions of two key lncRNAs and their pair subnetworks were evaluated in Huh7 cells compared to LX-2 cells. As shown in Figure 4(a), lncRNA EGOT was significantly higher in Huh7 cells than in LX-2 cells (), suggesting that lncRNA EGOT was upregulated in human hepatoma cells compared with human hepatic fibrosis cells in vitro studies. However, the expressions of miR-32-5p and XYLT2 were not significantly different in Huh7 cells and LX-2 cells ( and , respectively). Similarly, lncRNA SERHL and BCL2L1 were significantly higher expressed ( and , respectively), and miR-1269a was significantly lower expressed () in Huh7 cells than in LX-2 cells (Figure 4(b)). However, no significant differences were found in the other subnetwork’s expression, including miR-193b-3p, ARNT, CHST3, LPCAT1, and SYK (Figure 4(b)).

4. Discussion and Conclusions

Our study demonstrated the differential expression profiling of lncRNAs, miRNAs, and mRNAs in cirrhotic HCC. Next, based on miRanda and TargetScan, we explored the coexpression regulatory relationships between lncRNA, miRNA, and mRNA. Then, we built a functional ceRNA network to provide a global view of HCC-associated differentially expressed RNAs with interregulated 47 lncRNAs, 35 miRNAs, and 168 mRNAs. We also uncovered that two differentially expressed lncRNAs (EGOT and SERHL), and their pair subnetworks could affect the prognosis of cirrhotic HCC and predict the overall survival of patients with HCC. This finding provides a novel insight into the initiation and development of cirrhotic HCC.

Emerging evidences indicate that the aberrant expression of lncRNAs, developing into a highly active research hotspot, plays significant roles in tumorigenesis and acts as novel molecular biomarkers for the diagnosis and prognosis of HCC. More recently, a growing body of research shows that the construction of the ceRNA network reveals its diagnostic and prognostic value in HCC by comprehensive analysis [1417]. For example, lncRNA SSTR5-AS1, acting as a ceRNA, sponge mir-15b-5p to regulate carbonic anhydrase 2 (CA2) function in the progression of hepatitis B virus-related hepatocellular carcinoma [18].

In the present study, we highlighted lncRNAs EGOT and SERHL might be two critical lncRNAs that act as diagnostic and prognostic biomarkers in the progression of liver cirrhosis to HCC. Previous work has demonstrated that lncRNA EGOT is involved in the progress of breast cancer [19], gastric cancer [20], head and neck squamous cell carcinomas [21], renal cell carcinoma [22], and hepatitis C virus-induced liver tumors [23]. Moreover, higher expression of lncRNA SERHL is significantly correlated with patient’s shorter overall survival in HCC, not involving in further research on the functional implications of SERHL [24]. Our results suggest that lncRNA EGOT-miR-32-5p-XYLT2 axis and lncRNA SERHL-miR-1269a/miR-193b-3p-BCL2L1/SYK/ARNT/CHST3/LPCAT1 axis could contribute to the regulatory mechanism underlying the pathogenetic process of HCC. lncRNA EGOT synergistically increased the expression of hsa-miR-32-5p and then decreased the expression of XYLT2, relating to significant regulation of key HCC-related signaling pathways such as glycosaminoglycan biosynthesis and metabolic pathways (Fig. S4). lncRNA SERHL synergistically upregulated the expression of hsa-miR-1269a and hsa-miR-193b-3p and then downregulated the expression of BCL2L1, SYK, ARNT, CHST3, and LPCAT1, which are associated with the significant regulation of key HCC-related signaling pathways such as Ras signaling pathway, PI3K-Akt signaling pathway, NF-kappa B signaling pathway, and JAK-STAT signaling pathway (Fig. S4). Consistent with previous research, lncRNAs EGOT and SERHL are negatively associated with the overall survival of patients with cirrhotic HCC. These two lncRNAs may be exploited as potential effective biomarkers in the diagnosis and prognosis of HCC.

In conclusion, we identified the differential expression of cancer-specific lncRNAs, miRNAs, and mRNAs in the progression of liver cirrhosis to HCC and constructed the ceRNA network to reveal their genetic interactions. The higher-expressed SERHL was associated with a higher risk of all-cause death. The expression of SERHL-miR-1269a-BCL2L1 was significantly related to cirrhotic HCC in vitro studies. Importantly, as candidate biomarkers for HCC prognosis prediction, two potential unexplored lncRNAs EGOT and SERHL showed significantly negative associations with overall survival, and their pair subnetworks could provide a novel clue to future study into the underlying regulatory mechanism of cirrhotic HCC.

Abbreviations

HCC:Hepatocellular carcinoma
lncRNAs:Long noncoding RNAs
ceRNA:Competing endogenous RNA
miRNA:MicroRNA
mRNA:Messenger RNA
TCGA:The Cancer Genome Atlas
GEO:Gene Expression Omnibus
FC:Fold change
FBS:Fetal bovine serum
SD:Standard deviation
IQR:Interquartile range
VIF:Variance inflation factor
HR:Hazard ratios
CI:Confidence interval.

Data Availability

The datasets of lncRNAs, miRNAs, and mRNAs have been contributed by GSE17967, GSE21362, GSE63046, and TCGA-LIHC. The data of the current study are available from the corresponding author on reasonable request.

Conflicts of Interest

There are no potential conflicts of interest.

Authors’ Contributions

Yuli Zhang performed the experiments, analyzed the data, and wrote the manuscript. Dinggui Chen, Miaomiao Yang, Xianfeng Qian, and Chunmei Long performed the data processing and statistical analysis. Zhongwei Zheng designed the study, supervised the experiments, and wrote the manuscript.

Acknowledgments

The datasets of lncRNAs, miRNAs, and mRNAs have been contributed by GSE17967, GSE21362, GSE63046, and TCGA-LIHC. We thank all investigators for sharing these published data. This work was supported by the Changzhou Guidance Program of Applied Basic Research (CJ20209026 to Y.L.Z.).

Supplementary Materials

The file of Supplementary Materials. (Supplementary Materials)