A Novel lncRNA IHS Promotes Tumor Proliferation and Metastasis in HCC by Regulating the ERK- and AKT/GSK-3β-Signaling Pathways.

Long noncoding RNAs (lncRNAs) are involved in a variety of biological processes such as tumor proliferation and metastasis. A close relationship between hepatitis B virus X protein (HBx) and SMYD3 in promoting the proliferation and metastasis of hepatocellular carcinoma (HCC) was recently reported. However, the exact oncogenic mechanism of HBx-SMYD3 remains unknown. In this study, by performing lncRNA microarray analysis, we identified a novel lncRNA that was regulated by both HBx and SMYD3, and we named it lncIHS (lncRNA intersection between HBx microarray and SMYD3 microarray). lncIHS was overexpressed in HCC and decreased the survival rate of HCC patients. Knockdown of lncIHS inhibited HCC cell migration, invasion, and proliferation, and vice versa. Further study showed that lncIHS positively regulated the expression of epithelial mesenchymal transition (EMT)-related markers c-Myc and Cyclin D1, as well as the activation of the ERK- and AKT-signaling pathways. lncIHS exerted its oncogenic effect through ERK and AKT signaling. Moreover, results from transcriptome-sequencing analysis and mass spectrometry showed that lncIHS regulated multiple genes that were the upstream molecules of the ERK- and AKT-signaling pathways. Therefore, our findings suggest a regulatory network of ERK and AKT signaling through lncIHS, which is downstream of HBx-SMYD3, and they indicate that lncIHS may be a potential target for treating HCC.

Introduction

Hepatocellular carcinoma (HCC) is the third major cause of cancer mortality, and the global morbidity of the disease is increasing year by year. Among the therapeutic methods of HCC, surgical resection is the main treatment option, but about 80% of patients are preliminarily diagnosed with advanced stage and lose the operative chance.2, 3 Besides, the rate of postoperative recurrence and metastasis of HCC is high. Therefore, understanding the occurrence and development mechanism of HCC may provide new strategies for the diagnosis and treatment of this life-threatening disease.

Hepatitis B virus X protein (HBx) is a multifunctional regulator during HCC occurrence and development. Our previous studies have confirmed that HBx plays a key role in the molecular pathogenesis of hepatitis B virus (HBV)-HCC. On the other hand, our data and recent studies have indicated that SMYD3, an important methyltransferase in malignant tumor, is greatly associated with liver carcinogenesis and poor prognosis of HCC patients.7, 8, 9 Interestingly, HBx has been shown to not only upregulate SMYD3 expression in HCC but also interact with SMYD3 through a specific domain.10, 11 However, the exact mechanism of how HBx-SMYD3 promotes HCC occurrence and development remains unknown.

Long noncoding RNA (lncRNA) is a kind of RNA molecule that is over 200 nt, without a protein-encoding function. Studies have shown that lncRNAs can affect multiple biological processes of malignant tumor, such as cell proliferation and metastasis, and play an important role in the occurrence and development of HCC. Interestingly, HBx was recently reported to upregulate lncRNAs by recruiting EZH2, another important histone methyltransferase in malignant tumor. We hypothesize that HBx-SMYD3 may also promote HCC progression via lncRNAs.

In this study, by conducting microarray using HBx- and SMYD3-overexpressing HCC cells, we gained a novel lncRNA (LOC285740) and named it lncIHS (lncRNA intersection between HBx microarray and SMYD3 microarray). The role of lncIHS in HCC was then studied, which may provide a basis for clinical diagnosis and treatment of HCC.

Results

A Novel lncRNA Was Selected by Microarray

To explore lncRNAs regulated by both HBx and SMYD3, we first established HBx- and SMYD3-overexpressing HCC cells. We chose BEL-7402, Huh7, and SK-HEP-1 to establish SMYD3-overexpressing cell lines and MHCC-97H to establish the SMYD3 stable knockdown cell line (Figures S1A and S1B). SK-HEP-1 and HepG2 without HBV genome were chosen to establish HBx-overexpressing cell lines (Figure S1B). We confirmed the regulation of SMYD3 by HBx (Figure S1B).

Next, we used HBx- and SMYD3-overexpressing HCC cells to screen out lncRNAs that are regulated by both HBx and SMYD3 in HCC, and we found that five lncRNAs were upregulated by HBx-SMYD3, whereas three lncRNAs were downregulated by HBx-SMYD3 (Figure 1A). Considering the major function of SMYD3 is to trimethylate H3K4 and subsequently activate target gene expression, we chose lncRNAs upregulated by HBx-SMYD3 for the next investigation (Figure 1A). We first detected the expression of the five lncRNAs in HCC with tumor thrombus, HCC with non-tumor thrombus, and normal liver. We found that, compared with the other four lncRNAs, the expression of LOC285740 among the three groups was the most significant (Figure 1B), implying that LOC285740 was associated with HCC metastasis, which is the major oncogenic function of SMYD3 and HBx.8, 16 We therefore chose lncRNA LOC285740 for further investigation. Results from HBx-overexpressing, SMYD3-overexpressing, and SMYD3-knockdown HCC cells further confirmed that LOC285740 was positively regulated by HBx and SMYD3 (Figures 1C and S1C). Furthermore, knockdown of SMYD3 significantly reduced HBx-dependent expression of LOC285740 (Figure 1D). Because LOC285740 was initially screened out by both HBx microarray and SMYD3 microarray, we named LOC285740 as lncIHS (intersection between HBx microarray and SMYD3 microarray).

Figure 1

A Novel lncRNA Is Selected by Microarray

(A) lncRNAs were screened out by microarray using HBx- and SMYD3-overexpresseing HCC cell lines. (B) The expression of lncIHS in HCC with tumor thrombus was significantly more than others. (C) lncIHS was positively regulated by HBx and SMYD3. (D) HBx promoted lncIHS expression through SMYD3. Data are presented as mean ± SD for three independent experiments. *p < 0.05.

lncIHS Is Overexpressed in HCC Tissues and Correlates with Patient Clinicopathological Features and Unfavorable Survival

To explore the significance of lncIHS in HCC, we investigated its expression in HCC tissues. We observed that lncIHS expression was significantly increased in HCC tissues in comparison with para-tumor tissues (Figure 2A). High lncIHS expression in HCC was significantly associated with HBV infection, large tumor size, portal vein invasion, and high TNM stage (Table 1). Furthermore, HCC patients with high lncIHS expression had shorter overall and disease-free survival (Figures 2B and 2C). In addition, lncIHS expression was higher in HCC cell lines than normal liver cell line (Figure 2D). These observations indicated that lncIHS may act as an oncogene in HCC. The complete sequence of lncIHS was then identified by rapid amplification of cDNA ends (RACE) assay, and it is shown in Figure S2A. We used three different prediction tools, CPC2, coding potential assessment tool (CPAT), and PORTRAIT, to calculate the coding potential of lncIHS. lncIHS was classified as a noncoding sequence by the three prediction tools (coding probability, CPC2, 14%; CPAT, 17%; and PORTRAIT, 5.78%). In vitro translation analysis further confirmed the noncoding potential of lncIHS (Figure S2B).

Figure 2

lncIHS Is Overexpressed in HCC and Correlates with Unfavorable Survival

(A) Compared with para-tumor tissues, the expression of lncRNA was significantly increased in HCC tissues. (B and C) Kaplan-Meier analysis of overall survival (B) and disease-free survival (C) based on lncIHS expression levels in 105 cases of patients with HCC. The median level of lncIHS is used as the cutoff. (D) The expression of lncRNA was detected in a normal liver cell line and different HCC cell lines. (E) The levels of lncIHS mRNA in purified nuclear or cytoplasmic RNAs were detected using real-time PCR. U6 and β-actin served as nuclear and cytoplasmic controls, respectively. (F) lncIHS was located in the cytoplasm and nucleus, but mainly in the nucleus, as determined by FISH assays. Data are presented as mean ± SD for three independent experiments. *p < 0.05, **p < 0.01.

Table 1

Correlation between Expression of lncIHS and Clinicopathological Parameters in HCC Patients

Characteristics		lncIHS
	N = 105	Low Expression	High Expression	χ2	p Value
Age (Years)				0.06	0.801
<59	50	17	33
≥59	55	20	35
Gender				1.53	0.216
Male	51	21	30
Female	54	16	38
AFP				2.6	0.105
≤25 μg/mL	54	23	31
>25 μg/mL	51	14	37
HBV Infection				9.11	0.002*
(+)	55	12	43
(−)	50	25	25
Tumor Encapsulation				3.46	0.063
Present	64	27	37
Absent	41	10	31
Tumor Size (cm)				5.51	0.018*
≤5	49	23	26
>5	56	14	42
Tumor Stage				9.34	0.002*
I, II stage	47	24	23
III, IV stage	58	13	45
Portal Vein Thrombus				8.95	0.003*
Present	52	11	41
Absent	53	26	27

*p < 0.05.

lncIHS Is Located in the Nucleus of HCC Cells

To further investigate the oncogenic function of lncIHS in HCC, we first detected the location of lncIHS in HCC cell lines. RNA was extracted respectively from the cytoplasm and nucleus of SMMC-7721 for detecting the expression level of lncIHS. By comparing with U6 and β-actin, lncIHS was determined to be located mainly in the nucleus (Figure 2E). Our fluorescence in situ hybridization (FISH) assays further showed that lncIHS was located in the cytoplasm and nucleus, but mainly in the nucleus (Figure 2F). These data suggested that lncIHS may exert its function mainly in the nucleus.

lncIHS Promotes Cell Migration, Invasion, and Proliferation In Vitro and In Vivo

The exact function of lncIHS in HCC cells is unknown. We therefore performed gain- and loss-of-function study. For loss-of-function study, we chose SMMC-7721 and MHCC-97H to establish lncIHS stable knockdown cell lines (Figure 3A; Figure S3A). Because the above results showed significantly high expression of lncIHS in HCC with tumor thrombus, we first focused on the association between lncIHS and the migration or invasive ability of HCC cells. Migration and invasion assays showed that knockdown of lncIHS significantly suppressed cell migration and invasion (Figure 3B; Figure S3B). We also measured the effect of lncIHS on the proliferation of HCC cells. Results from the MTT assay and colony formation assay showed that lncIHS knockdown significantly decreased the proliferation of SMMC-7721 and MHCC-97H (Figures 3C and 3D; Figures S3C and S3D). Furthermore, knockdown of lncIHS hindered cell cycle transition from G1 to S stage (Figure 3E; Figure S3E). We confirmed these results using transient transfection in both cell lines (Figures S4 and S5).

Figure 3

Knockdown of lncIHS Inhibits SMMC-7721 Proliferation, Migration, and Invasion In Vitro

(A) lncIHS knockdown was confirmed by real-time PCR in SMMC-7721. (B) Knockdown of lncIHS suppressed the migration and invasion abilities of SMMC-7721, as determined by transwell assay. (C and D) In vitro proliferative ability was decreased in SMMC-7721 with lncIHS knockdown, as determined by MTT assay (C) and colony formation assay (D). (E) Knockdown of lncIHS inhibited cell cycle transition from G1 to S in HCC cells. Data are presented as mean ± SD for three independent experiments. *p < 0.05.

To further confirm the above in vitro results, in vivo experiments were performed. lncIHS knockdown decreased xenografted tumor growth, tumor weight, and tumor size in vivo (Figures 4A–4C). In addition, results from Ki67 staining of xenograft tumor tissues indicated that tumors derived from lncIHS-silenced cells showed significantly reduced proliferation (Figures 4D and 4E). In vivo effects of lncIHS on HCC metastasis were also detected. The incidence of lung metastasis in the lncIHS-knockdown group was significantly less than that in the control group (Figure 4F). These results indicated that lncIHS is required for HCC cell migration, invasion, and proliferation.

Figure 4

Knockdown of lncIHS Inhibits Tumor Growth and Lung Metastasis of HCC Cells in Nude Mice

(A) SMMC-7721 cells (5 × 106) with stable knockdown of lncIHS were inoculated into nude mice, and the effect of lncIHS on HCC tumor growth was examined after 6 weeks (n = 5). A photograph of the tumors is presented. (B and C) Effect of lncIHS on HCC growth was described by tumor weight (B) and tumor size (C) in the two groups. (D) Representative images of immunohistochemistry (IHC) staining of Ki67 showed that lncIHS inhibition decreased tumor proliferation in xenografted tumors. (E) The expression of lncRNA in the tumor xenografts of the sh-lncIHS group was significantly lower than in the sh-NC group. (F) Representative H&E staining of lung metastases between SMMC-7721 sh-lncIHS cells and control cells. Data are presented as mean ± SD for three independent experiments. *p < 0.05, **p < 0.01.

Next, we questioned whether overexpression of lncIHS could strengthen the ability of migration, invasion, and proliferation of HCC cells. BEL-7402 and PLC/PRF/5 cells were transiently transfected with a vector containing the complete sequence of lncIHS (Figure 5A). Consistent with the results from lncIHS-knockdown cell lines, lncIHS overexpression significantly enhanced the migration, invasion, and proliferation capacities of BEL-7402 and PLC/PRF/5 cells (Figures 5B–5E). These results were also confirmed in stable lncIHS-overexpressing cells (Figures S6A–S6C). Taken together, these results suggested that lncIHS had an important role in promoting HCC cell migration, invasion, and proliferation.

Figure 5

Overexpression of lncIHS Promotes HCC Cell Migration, Invasion, and Proliferation

(A) Overexpression of lncIHS was confirmed by PCR. (B and C) Overexpression of lncIHS promoted the migration and invasion abilities of HCC cells, as determined by wound-healing (B) and transwell assays (C). (D) In vitro proliferative ability of HCC cells was significantly promoted by lncIHS by MTT assay. (E) lncIHS overexpression stimulated cell cycle transition from G1 to S in HCC cells. Data are presented as mean ± SD for three independent experiments. *p < 0.05, **p < 0.01.

lncIHS Regulates the Expressions of the EMT and Cell Cycle-Associated Genes and the Activation of ERK- and AKT/GSK-3β-Signaling Pathways

As our above results showed that lncIHS significantly promoted HCC cell growth, migration, and invasion, we sought to determine whether lncIHS regulated the expression of the epithelial mesenchymal transition (EMT) and cell cycle-associated genes in HCC cells. Knockdown of lncIHS decreased the expression of c-Myc, Cyclin D1, as well as mesenchymal markers N-cadherin and Vimentin, and it increased the expression of epithelial marker E-cadherin in SMMC-7721 and MHCC-97H (Figure 6A; Figure S7A). In contrast, overexpression of lncIHS increased the expression of c-Myc, Cyclin D1, N-cadherin, and Vimentin, and it deceased E-cadherin expression (Figure 6B; Figure S7B).

Figure 6

lncIHS Regulates the Expression of the EMT and Cell Cycle-Associated Genes and the Activation of the ERK- and AKT/GSK-3β-Signaling Pathways

(A and B) Knockdown of lncIHS reduced the expression of c-Myc, Cyclin D1, N-cadherin, and Vimentin, and it increased the level of E-cadherin in SMMC-7721 and MHCC-97H (A), and vice versa in BEL-7402 and PLC/PRF/5 (B). (C and D) Knockdown of lncIHS decreased the expression of p-MEK1/2, p-ERK1/2, p-AKT, and p-GSK-3β in SMMC-7721 and MHCC-97H (C), and vice versa in BEL-7402 and PLC/PRF/5 (D).

AKT- and ERK-signaling pathways have been well known as the regulators of cell cycles and the EMT.17, 18, 19 We detected the effect of lncIHS on PI3K/AKT- and ERK-signaling pathways. Upon knockdown of lncIHS, the expressions of phosphorylated MEK1/2, ERK1/2, AKT, and GSK-3β were downregulated, whereas total MEK1/2, ERK1/2, AKT, and GSK-3β was kept unchanged (Figure 6C; Figure S7C). On the contrary, lncIHS overexpression enhanced the phosphorylation of MEK1/2, ERK1/2, AKT, and GSK-3β (Figure 6D; Figure S7D). These results indicated that lncIHS may promote HCC cell migration, invasion, and proliferation through the activation of ERK and AKT signaling and their downstream molecules, such as c-Myc, Cyclin D1, and EMT-associated genes.

lncIHS Promotes HCC Cell Migration, Invasion, and Proliferation through ERK- and AKT-Signaling Pathways

PI3K inhibitor Ly294002 and MEK1/2 inhibitor U0126 were applied to study the role of ERK and AKT signaling in lncIHS-mediated migration, invasion, and proliferation of HCC. The effects of Ly294002 and U0126 on the phosphorylation of AKT, GSK-3β, MEK1/2, and ERK1/2 were validated (Figure 7A; Figure S8A). Inactivation of ERK and AKT signaling in lncIHS-overexpressing cells not only reversed lncIHS-dependent decrease in E-cadherin expression and lncIHS-dependent increases in N-cadherin, Vimentin, c-Myc, and Cyclin D1 expression but also reduced lncIHS-enhanced cell migration, invasion, and proliferation (Figures 7B–7D; Figure S8B). To further verify the effects of ERK and AKT signaling on mediating the oncogenic function of lncIHS, we performed tumor xenograft experiments using stable lncIHS-overexpressing cell lines with or without LY294002 and U0126 treatment. ERK- and AKT-signaling inhibition significantly suppressed the tumor growth induced by lncIHS overexpression (Figures S6D–S6F). These data suggested that lncIHS promoted HCC cell migration, invasion, and proliferation through the MEK/ERK- and AKT/GSK-3β-signaling pathways.

Figure 7

lncIHS Promotes HCC Cell Migration, Invasion, and Proliferation through the ERK- and AKT-Signaling Pathways

(A) The effects of Ly294002 and U0126 on the phosphorylation of AKT, GSK-3β, MEK1/2, and ERK1/2 were validated. (B) Treatments with Ly294002 and U0126 in HCC cells reversed the lncIHS-dependent decrease of E-cadherin expression and the lncIHS-dependent increase of N-cadherin, Vimentin, c-Myc, and Cyclin D1 expressions. (C and D) Inactivation of ERK and AKT signaling significantly reduced lncIHS-enhanced cell migration, invasion, and proliferation, as determined by transwell assay (C) and MTT assay (D). Data are presented as mean ± SD for three independent experiments. *p < 0.05, **p < 0.01.

Next-Generation Sequencing Reveals a lncIHS-Dependent Regulatory Network of ERK- and AKT-Signaling Pathways

Our above results showed that lncIHS positively regulated ERK- and AKT-signaling pathways. However, our above findings also indicated that lncIHS was mainly located in the nucleus of HCC cells, which suggested that lncIHS did not directly regulate ERK and AKT signaling. We hypothesized that lncIHS could regulate genes upstream of the ERK and AKT pathways on the mRNA level. We therefore performed next-generation sequencing (NGS) using SMMC-7721 sh-lncIHS and sh-NC cells to screen out genes that are regulated by lncIHS on the mRNA level. We constructed Gene Ontology (GO) trees and pathway-act-network based on the differentially expressed genes between the two groups.

Consistent with our above findings, results of the GO trees showed there were many genes regulated by lncIHS enriched in mitogen-activated protein kinase (MAPK) signaling, cell proliferation, and cell migration (Figure S9). Pathway-act-network revealed that the MAPK-signaling pathway and PI3K/AKT-signaling pathway were core nodes in the lncIHS-associated pathway network (Figure S10). Genes previously reported to be the upstream molecules of the ERK or AKT pathway were then selected and validated by real-time PCR in both lncIHS-knockdown and -overexpression HCC cells (Figures 8A and 8B).

Figure 8

NGS Reveals a lncIHS-Dependent Regulatory Network of the ERK- and AKT-Signaling Pathways

(A) NGS screened out several genes that were upregulated or downregulated by lncIHS. (B) The lncIHS-related regulation of genes that may be the upstream molecules of the AKT- and ERK-signaling pathways was validated by real-time PCR. (C and D) Knockdown of lncIHS decreased MAP3K8 expression and increased DUSP5 and DUSP10 expressions in SMMC-7721 and MHCC-97H (C), and vice versa in BEL-7402 and PLC/PRF/5 (D). Data are presented as mean ± SD for three independent experiments. *p < 0.05.

We paid special attention to mitogen-activated protein kinase kinase kinase 8 (MAP3K8) and dual specificity phosphatase 5 and 10 (DUSP5 and DUSP10), which are kinase or phosphatase and have been well known as the upstream regulators of ERK- or AKT-signaling pathways.20, 21, 22, 23 The positive regulation of MAP3K8 and negative regulation of DUSP5/DUSP10 by lncIHS were further validated by western blot (Figures 8C and 8D; Figures S11A and S11B). We also observed a positive correlation between lncIHS and MAP3K8, as well as a negative correlation between lncIHS and DUSP5 or DUSP10 in HCC tissues (Figure S12). Furthermore, we employed lncRNA-WEE2-AS1, which we recently reported could promote HCC invasion through activating the PI3K/AKT/GSK-3β-signaling pathway, to see whether MAP3K8, DUSP5, and DUSP10 were regulated by lncRNA-WEE2-AS1. Overexpression of lncRNA-WEE2-AS1 did not affect the expressions of MAP3K8, DUSP5, and DUSP10 (Figure S13), suggesting the functional specificity of lncIHS. Taken together, our data suggested a complex regulation mechanism of the ERK- and AKT-signaling pathways by lncIHS in HCC.

lncIHS Regulated MAP3K8, DUSP5, and DUSP10 through Interacting with YBX1

To further identify the mechanisms underlying lncIHS-induced ERK and AKT pathway activation, we first examined whether lncIHS affected the mRNA expression levels of its cis genes. lncIHS is located on human chromosome 6q24.2, in antisense of the intron of a protein-coding gene PHACTR2 (Figure S14A). We found the expression of PHACTR2 was not regulated by lncIHS (Figures S14B and S14C). Besides, the expressions of other neighboring genes, such as FUCA2, LTV1, ADAT2, and PEX2, were also not regulated by lncIHS (Figures S14A–S14C). These results implied that lncIHS does not function in cis but in trans to activate the ERK and AKT pathways.

We then performed RNA pull-down assay using a biotin-labeled lncIHS RNA probe to identify potential binding proteins of lncIHS in HCC. The proteins pulled down by lncIHS were then analyzed by mass spectrometry. Among the proteins identified by mass spectrometry, we paid special attention to YBX1 (Figure 9A). YBX1 is a nuclear and cytoplasmic protein that can bind DNA and RNA by recognizing special nucleic acid sequences.25, 26 YBX1 has been demonstrated to serve as a transcriptional activator as well as mRNA stabilizer in cancer.27, 28 We further confirmed the interaction of lncIHS with YBX1 by RNA pull-down and RNA immunoprecipitation (RIP) assays (Figures 9B and 9C).

Figure 9

lncIHS Interacts with YBX1 in HCC Cells

(A) YBX1 was a potential interactive candidate of lncIHS, as determined by RNA pull-down assay, and the enriched products were eluted and separated by SDS-PAGE and Coomassie blue staining. (B) Western blot of the proteins from antisense lncIHS and lncIHS pull-down assays. (C) The interactions of lncIHS with YBX1 were verified by an RIP assay. (D–F) The mRNA expressions of MAP3K8 (D), DUSP5 (E), and DUSP10 (F) were detected in BEL-7402-lncIHS and its control cells with or without YBX1 knockdown. (G–I) RNA stability assays were performed in BEL-7402-lncIHS and its control cells with or without YBX1 knockdown, and the degradation rates of the MAP3K8 (G), DUSP5 (H), and DUSP10 (I) mRNA were measured at the indicated time points. (J and K) YBX1 potential binding sites on the MAP3K8 gene promoter (J) and the DUSP5/DUSP10 3′ UTR (K) are shown. (L) YBX1 was mainly located in the nucleus in lncIHS-overexpressing cells. Data are presented as mean ± SD for three independent experiments. *p < 0.05, **p < 0.01.

Next, we questioned whether YBX1 mediates the lncIHS-dependent regulation of MAP3K8, DUSP5, and DUSP10. Results from real-time PCR showed that knockdown of YBX1 decreased the expressions of MAP3K8, DUSP5, and DUSP10 in HCC cells without lncIHS overexpression (Figures 9D–9F). On the other hand, in cells with lncIHS overexpression, YBX1 knockdown more significantly decreased MAP3K8 expression, but not DUSP5 and DUSP10 (Figures 9D–9F). We then employed actinomycin D to inhibit gene transcription. Consistent with the real-time PCR results, DUSP5 and DUSP10 mRNA half-lives were significantly decreased in cells either with lncIHS overexpression or with YBX1 konckdown (Figures 9H and 9I). Interestingly, the MAP3K8 mRNA half-life was influenced by neither lncIHS overexpression nor YBX1 knockdown (Figure 9G).

Considering the nuclear localization of lncIHS and the molecular function of YBX1, we hypothesized that lncIHS could bind and sequester YBX1 in the nucleus, facilitating the recruitment of YBX1 on the promoter of the MAP3K8 gene, leading to transcriptional activation of MAP3K8. Meanwhile, the cytoplasmic expression of YBX1 was downregulated, leading to the instability of DUSP5/DUSP10 mRNA. Consistently, we found many YBX1 potential binding sites on the MAP3K8 gene promoter region and DUSP5/DUSP10 3′ UTR (Figures 9J and 9K) and that overexpression of lncIHS resulted in YBX1 nuclear accumulation (Figure 9L). Taken together, our results suggested that lncIHS may upregulate MAP3K8 and downregulate DUSP5 and DUSP10 through its binding protein YBX1 in HCC.

Discussion

HCC is one of the most common malignant tumors worldwide. In China, nearly 80% of HCC cases are attributable to HBV. HBx, one of the central viral proteins produced by HBV, is now well established to be involved in hepatocarcinogenesis due to pleiotropic effects on several signaling pathways regulating cell death, proliferation, differentiation, invasion, and metastasis; oxidative stress; and DNA repair.5, 30 Recently, many lncRNAs (e.g., Dreh and HOTAIR) have been reported to be associated with HBx in HCC, suggesting lncRNAs play an important role in the oncogenic process of HBx. SMYD3 is a methyltransferase, which mainly leads to trimethylation of the lysine 4 of histone H3 (H3K4), makes the spatial structure of chromosomes in an open state, and finally promotes target gene transcription.8, 33 SMYD3 is highly expressed in a variety of malignant tumors, such as HCC, leukemia, bladder cancer, and glioma. The overexpression of SMYD3 in malignant cells is associated with the proliferation, anti-apoptosis, autophagy, migration, and invasion abilities of tumor cells.38, 39

Given the potent roles of HBx and SMYD3 in HCC, it is important to enrich the oncogenic mechanism of both proteins. Interestingly, studies have shown that the expression of SMYD3 is regulated by HBx in HCC cells. Moreover, recent studies also indicate that HBx interacts with SMYD3 through a specific domain and the interaction between HBx and SMYD3 activates downstream signaling through AP-1. These studies suggest that HBx-SMYD3 is an important pathway in promoting tumorigenesis and progression of HCC. However, the downstream mechanism of HBx-SMYD3 in HCC needs to be further elucidated.

Many studies have shown that lncRNAs are closely associated with the occurrence, development, invasion, and metastasis of human malignant tumors.12, 13 However, there are few studies focusing on lncRNAs regulated by both HBx and SMYD3 in HCC. By conducting lncRNA microarray screening on HCC cell models overexpressing HBx or SMYD3 and subsequent validation in HCC cell lines, we found a novel lncRNA that was upregulated by HBx-induced overexpression of SMYD3, named it lncIHS, and studied its role in HCC.

It is well known that PI3K/AKT and MEK/ERK signaling are the major signaling pathways in the tumorigenesis and tumor progression of HCC. Numerous researchers have demonstrated that the activation of the PI3K/AKT- and ERK-signaling pathways is an important cause of EMT and cell cycle transition of G1/S in cancer. In this study, we showed that lncIHS could activate PI3K/AKT and ERK signaling. Moreover, by employing the inhibitors Ly294002 and U0126, we found that the activation of AKT and ERK was required for the lncIHS-dependent enhancement of HCC cell proliferation, migration, and invasion, as well as the EMT and cell cycle-related protein expressions, indicating that lncIHS is an important oncogenic noncoding RNA in HCC.

Regulation of the MEK/ERK and PI3K/AKT pathways is mediated by a series of kinases, phosphatases, and other factors. The synergistic alteration of these kinases or phosphatases will finally result in the pathway activation or inactivation. Based on this background knowledge, we applied NGS to screen out genes that are regulated by lncIHS on the mRNA level and, meanwhile, are the upstream regulators of the AKT- and ERK-signaling pathways. We found many genes that have been reported to be associated with the AKT- and ERK-signaling pathways in the NGS results. Among them, MAP3K8, DUSP5, and DUSP10 have been well studied as the upstream molecules of the ERK- and AKT-signaling pathways. MAP3K8 is a member of the serine or threonine protein kinase family. MAP3K8 controls several signaling pathways, including MEK/ERK and PI3K/AKT, and it functions as an oncogene in many human cancers. In contrast, DUSP5 and DUSP10, both of which are members of the DUSP family, mainly contribute to the dephosphorylation of ERK.23, 43 Both DUSP5 and DUSP10 serve as tumor suppressors in human cancers.23, 44 Hence, our results indicated that lncIHS positively regulates the AKT- and ERK-signaling pathways through multiple mechanisms.

lncRNAs function in a wide range of biological processes, and they regulate gene expression in cis or in trans via various mechanisms. Actually, many recent studies have reported a series of lncRNAs that function in trans, but not in cis.45, 46 Consistently, our data showed the funtions of lncIHS may be independent of its neighboring genes but dependent upon its binding protein, YBX1. YBX1 is not only a DNA-binding protein that recognizes the Y-box sequence in the nucleus but also a major component of cytoplasmic messenger ribonucleoprotein particles (mRNPs), with roles in pre-mRNA splicing and mRNA stability. Thus, different localization of YBX1 in cells exhibits different functions.

In this study, we found that YBX1 may be a functional partner of lncIHS according to the following findings: (1) lncIHS and YBX1 could bind to each other; (2) knockdown of YBX1 under the condition of low lncIHS expression only slightly influenced MAP3K8 mRNA expression, but it significantly downregulated DUSP5 and DUSP10 mRNA stability and expression; in contrast, under the condition of high lncIHS expression, knockdown of YBX1 significantly decreased MAP3K8 expression, but it had no effect on DUSP5 and DUSP10 mRNA stability and expression; (3) YBX1 did not regulate MAP3K8 mRNA stability; (4) many YBX1 potentail binding sites were found on the MAP3K8 gene promoter and DUSP5/DUSP10 3′ UTRs; and (5) overexpression of lncIHS accumulated YBX1 in the nucleus. All these clues lead to the hypothesis that lncIHS could bind and sequester YBX1 in the nucleus, thus leading to transcriptional activation of MAP3K8 and the instability of DUSP5/DUSP10 mRNA. These reasons could partly explain how lncIHS regulates ERK and AKT/GSK-3β signaling in HCC. However, further study is needed to illustrate a more detail mechanism.

In conclusion, our study reported a novel oncogenic lncRNA, lncIHS, which is positively regulated by HBx-induced overexpression of SMYD3 in HCC. The expression of lncIHS is upregulated in HCC tissues and is markedly associated with a poor prognosis in HCC patients. lncIHS promotes HCC cell proliferation and metastasis by activating the AKT/GSK-3β- and ERK-signaling pathways. More importantly, lncIHS, together with its binding protein YBX1, regulates ERK and AKT/GSK3-β signaling through multiple mechanisms, forming a regulatory network. Thus, lncIHS may be a candidate biomarker for HCC prognosis and a potential therapeutic target, especially in HBV-associated HCC with positive SMYD3 expression.

Materials and Methods

Cell Culture and Reagents

Human hepatocyte LO2 and HCC cell lines SMMC-7721, MHCC-97H, BEL-7402, PLC/PRF/5, SK-HEP-1, and HepG2 were purchased from the Chinese Academy of Sciences Cell Bank Type Culture Collection. The cells were cultured with DMEM and RPMI-1640 (Gibco, Carlsbad, CA) together with 10% fetal bovine serum (Gibco) at 37°C with a humidity of 5% CO2. PI3K inhibitor (LY294002) and ERK1/2 inhibitor (U0126) were purchased from Selleck Chemicals.

Patient Sample

105 fresh HCC tissues and para-tumor tissues were obtained from HCC patients at the time of 2010–2012 after surgery from Sun Yat-Sen Memorial Hospital. This study was conducted in accordance with the Declaration of Helsinki. The study protocol was approved by the research ethics committee of Sun Yat-Sen Memorial Hospital, and written informed consent was obtained from each patient. The diagnosis of HCC was based on pathology and dynamic contrast-enhanced imaging (computed tomography scan or MRI).

cDNA Preparation and Real-Time PCR

Total RNA was extracted from HCC cell lines, 105 fresh HCC tissues, and the para-tumor tissues with Trizol (Takara, China), according to the manufacturer’s instructions. Reverse transcription was performed using Prime Script RTase (Takara, China), according to the manufacturer’s protocol. The expression levels of mRNA were determined with real-time PCR using SYBR green (Takara, China), according to the manufacturer’s instructions. The primers used for real-time PCR are given in Table S1.

The Sequence of lncIHS

5′-RACE and 3′-RACE were performed using SMART RACE cDNA amplification kit (Clontech Laboratories, USA), according to the manufacturer’s instructions.

Establishment of Stable Knockdown and Overexpression Cell Lines

The HepG2 and SK-HEP-1 cell lines were used to construct HBx-overexpressing cells. The BEL-7402 cell line was used to construct SMYD3-overexpressing cells. The MHCC-97H cell line was used to construct SMYD3-knockdown cells. HEK293T packaging cells were transfected with the appropriate retroviral construct. Culture supernatants were collected at 48 h after transfection, and we infected the target cells. Cells were then selected with puromycin for 2 weeks. The target for SMYD3 small hairpin RNA (shRNA) was 5′-AGCCTGATTGAAGATTTGATTCT-3′.

lncIHS shRNA vectors were purchased from Obio Technology (Shanghai, China). The lentivirus of lncIHS shRNA was produced from HEK293T packaging cells and infected HCC cells (SMMC-7721 and MHCC-97H), as mentioned above. The targets for lncIHS shRNA were sh-IHS1, 5′-CCATCCTGTGGCTATGATT-3′; and sh-IHS2, 5′-GCACGAAACATCTCTCCAA-3′.

The full lengths of lncIHS sequences were sub-cloned into pEX-2 plasmid, and the recombined plasmid was referred to as pEX-IHS. pEX-IHS and pEX-2 were, respectively, transfected into BEL-7402 cells using the Lipofectamine 3000 reagent, according to the manufacturer’s instructions. Stably transfected cells were obtained after selection in culture medium containing 0.5 mg/mL G418 (Invitrogen, USA) for approximately 6 weeks.

Nuclear and Cytoplasmic RNA Isolation

HCC cells were prepared in a 150-mm plate, centrifuged, and the supernatant was discarded. RNA subcelluar isolation kit was purchased from Active Motif (Shanghai, 25501). Following the assay protocol, cytoplasmic RNA and nuclear RNA were extracted separately. The expression levels of RNA were determined by real-time PCR using SYBR green.

FISH Assays

4% paraformaldehyde was added to an SMCC-7721 cell sample; then the fixed cell was placed on the glass slide, immersed in ethanol solution, and dried. Further, 10 μL hybridization buffer and 1 μL probes were added to the sample on the sheet glass at the 46°C incubator. After 1.5 h of hybridization, the sheet glass was quickly put into the 48°C washing buffer for 30 min. Then the sheet glass was rinsed with ultrapure water and dried. The sample was covered with the proper amount of anti-fading reagent and the cover glass, and then it was placed under the confocal laser-scanning microscope for observation.

MTT Assays

2,000 HCC cells were added to each of the 5 wells of the 96-well plates. Then the cells were routinely cultured for 1–5 days according to the experimental time. At each time point, MTT solution was added to each well and incubated for 4 h at 37°C. DMSO was then used to dissolve the purple crystal substances. The absorbance at 570 nm was measured with a microplate reader, and the growth curve was drawn and analyzed.

Colony Formation Assay

1,000 HCC cells were placed in six-well plates. The cells were cultured with growth medium and 10% fetal bovine serum (FBS) for 2 weeks. The cells were fixed and stained with paraformaldehyde and crystal violet, respectively.

Cell Cycle Analysis

HCC cells were gathered, washed with PBS twice, fixed in 75% ice-cold alcohol, and then incubated with PBS containing propidium iodide (PI, 10 μg/mL) and RNase A (0.5 mg/mL). The cells were analyzed with a flow cytometer.

Wound-Healing Assay

HCC cells were cultured in six-well plates with 5% CO2 at 37°C. When the cells were full, three parallel scratch wounds were made across each well using a P-100 pipette tip. The cells were continuously incubated with fresh growth medium without FBS. Wound closure was observed at 0 and 48 h under a microscope.

Transwell Migration and Invasion Assays

For the migration assay, 1 × 105 cells in medium without FBS were placed on a polycarbonate membrane insert in a transwell apparatus (Corning, USA). After 24-h incubation, the bottom cells were removed with a cotton swab. The cells on the lower side of the chamber were analyzed under a microscope. The process for the invasion assay was the same as the migration assay, except that the transwell membranes were pre-treated with 24 mg/mL Matrigel (R&D Systems, USA).

In Vivo Experiments

The cells were placed in 10-cm-diameter plates, collected, and re-suspended with PBS. A total of 5 × 106 cells was subcutaneously injected into the right shoulder of BALB/c nude mice. Tumor growth was monitored with tumor weight and tumor volume, which was calculated using the formula V = W2 × L × 0.5. Meanwhile, 1 × 106 cells were injected into the caudal vein of nude mice. After 6 weeks, the mice were euthanized, tumors were excised and weighed, and the lung tissues were collected for H&E staining. To further explore the potential role of inhibitors in HCC progress, a dose of 10 mg/kg body weight U0126 or 75 mg/kg body weight LY294002 was injected intraperitoneally into mice 1 week after inoculation, three times per week for 4 weeks. The animal experiments of this study were approved by the animal ethics committee of Sun Yat-sen University.

Western Blot

All cells were lysed by cell lysis buffer (p0013 RIPA lysis, Beyotime). The BCA method was used to detect protein concentration. The proteins were separated by the SDS-PAGE and transferred to polyvinylidene fluoride (PVDF) membranes (Millipore, USA). Then, the PVDF membranes were incubated at 4°C overnight with the following primary antibodies: AKT, p-AKT, GSK-3β, p-GSK-3β, ERK, p-ERK, E-cadherin, Vimentin, N-cadherin (Cell Signaling Technology, 1:1,000), MAP3K8, DUSP5, DUSP10 (ABclone, 1:500), c-Myc, Cyclin D1 (Santa Cruz Biotechnology, 1:500), and GAPDH (Proteintech, 1:500). Then PVDF membranes were incubated with the appropriate secondary antibody (Santa Cruz Biotechnology, 1:10,000) for 1 h at room temperature. Exposure was detected by the chemiluminescence solution.

RNA Stability Assay

To analyze RNA stability, cells were treated with actinomycin D (1 μg/mL). Cells were collected at different time points, and the RNA was extracted using Trizol reagent and detected by real-time PCR.

RNA Pull-Down Assay

Cell nuclear lysates were incubated with biotinylated full-length lncIHS and streptavidin beads for RNA pull-down incubation. Then beads were collected by centrifugation. RNA-associated proteins were eluted and resolved by SDS-PAGE, followed by silver staining.

RIP

RIP assays were performed using an RNA-binding protein immunoprecipitation kit (Millipore), according to the manufacturer’s instructions. RIP products were analyzed by real-time PCR using total RNA as input controls.

Statistical Analysis

All data analysis was performed with SPSS software (version 18.0). The experimental data are represented by mean ± SD. Kaplan-Meier analysis and log rank test were used to determine cumulative survival. Student’s t test and Fisher’s exact test were used to analyze quantitative data and categorical data, respectively. p < 0.05 was considered statistically significant.

Author Contributions

C.H., Z.Z., and Z.X. conceived and designed the study. Z.C., W.Y., and Q.Z. carried out the experiments. H.J. and W.Y. collected and analyzed the clinical data. Z.C. and Z.Z. wrote the manuscript. D.H., J.W., and J.Z. revised the manuscript. All authors read and approved the final manuscript.

Conflicts of Interest

The authors declare no competing interests.