Skip to main content

LINC01287/miR-298/STAT3 feedback loop regulates growth and the epithelial-to-mesenchymal transition phenotype in hepatocellular carcinoma cells

Abstract

Background

The long non-coding RNAs (lncRNAs) have participated in the promotion of hepatocellular carcinoma (HCC) initiation and progression. Nevertheless, the biological role and underlying mechanism of LINC01287 in HCC has never been reported.

Methods

The TGCA database was used to explore the abnormal expression of lncRNAs in HCC. Real-time PCR and in situ hybridization assays were used to examine the expression of LINC01287 in HCC tissues. The clinicopathological characteristics of HCC patients in relation to LINC01287 expression were then analyzed. Infection of cells with the si-LINC01287 lentiviral vector was performed to down-regulate LINC01287 expression in HCC cells. MTT and colony formation assays were performed to examine cell growth ability, and FACS analysis was performed to examine the cell cycle and apoptosis. A Boyden assay was used to examine HCC cell invasion ability, and RNA immunoprecipitation tested the interaction between LINC01287 and miR-298. A luciferase reporter assay was used to examine whether STAT3 was a direct target of miR-298, and chromatin immunoprecipitation (ChIP) was used to examine the potential binding of c-jun to the miR-298 promoter.

Results

We revealed that the expression of LINC01287 was increased in HCC cell lines, as well as tissues. Knockdown of LINC01287 decreased HCC cell growth and invasion both in vitro and in vivo. LINC01287 can negatively regulate miR-298 expression by acting as a ceRNA. miR-298 directly targeted STAT3 and inhibited its expression. LINC01287 exerted its function via the miR-298/STAT3 axis in HCC. Interestingly, STAT3 elevated LINC01287 expression via c-jun, which bound to the LINC01287 promoter. A feedback loop was also discovered between LINC01287 and the miR-298/STAT3 axis.

Conclusions

Our data indicated that LINC01287 played an oncogenic role in HCC growth and metastasis and that this lncRNA might serve as a novel molecular target for the treatment of HCC.

Background

Ranking the fifth most common cancer, Hepatocellular carcinoma (HCC) is the third leading cause of cancer-related death worldwide [1]. Lots of HCC patients are at the late stages of the disease when they are diagnosed. In addition, there is a high frequency of tumor recurrence in HCC patients although they have surgical resection and this contributes to the poor prognosis of HCC [2]. Therefore, it is urgent in need to develop novel strategies for the diagnosis and treatment of HCC.

Approximately 2% of the human genome accounts for protein coding genes. However, the majority of transcripts consists of non-coding RNAs (ncRNAs) [3] that can be grouped into the following classes depending on their transcript size: long non-coding RNAs (lncRNAs) and small ncRNAs [4]. LncRNAs are transcripts with a length greater than 200 nucleotides (nt). The abnormal expression of lncRNAs often contributes to tumor initiation, growth and metastasis [5]. Many studies have demonstrated that lncRNAs are dysregulated in many cancers, including HCC [6]. For instance, lncRNA CRNDE is upregulated in HCC and significantly associated with poor clinical outcomes, and knockdown of its expression impairs cell proliferation and invasion [7]. In addition, lncRNA TUG1 contributes to HCC cells proliferation, migration and tumorigenesis via interacting with miR-144 [8]. Moreover, NEAT1 upregulates TGF-β1 to induce HCC progression by sponging hsa-miR-139-5p [9]. Although large numbers of lncRNAs have been annotated, the role and molecular regulatory mechanisms of lncRNAs in HCC still require further clarification.

MicroRNAs (miRNAs), a group of small and non-coding RNAs, regulate down-stream targets expression via modulation of post-transcriptional [10].Mounting evidence has demonstrated that miRNAs regulate diverse biological process, including cancer cell proliferation, apoptosis and invasion [11]. It is worth to note that lncRNAs can function as ceRNAs that compete for miRNA binding, and therefore, they derepress the expression of miRNA-targeted mRNAs [12]. The lncRNA-miRNA-mRNA regulatory network has been implicated to regulate tumorigenesis [13, 14].

In the current study, we explore the role and underlying mechanism of LINC01287 in HCC. We also explore the interaction between LINC01287 and miR-298.

Methods

Cell culture and collection of HCC patient samples

HCC cell lines (HepG-2, Huh7, Bel7402 and Hep3B) and the normal liver epithelial cell line LO2 were purchased from the Institute of Biochemistry and Cell Biology of the Chinese Academy of Sciences (Shanghai, China). All cell lines were maintained at 37 °C in a humidified 5% CO2 atmosphere in RPMI-1640 medium supplemented with 10% fetal bovine serum.

HCC tissues and matched normal tissue samples (98 pairs, formalin-fixed and paraffin-embedded) were obtained from patients at Gaozhou People’s Hospital. These patients were diagnosed with hepatocellular carcinoma in our hospital from Sep. 2011 to Nov. 2016. Written informed consent was obtained from all patients, and the project was approved by the Ethical and Scientific Committees of Gaozhou People’s Hospital. All experimental protocols were approved by the Clinical Research Ethics Committees of Gaozhou People’s Hospital.

Cell transfection, lentivirus production and transduction

The cell transfections were performed using lipofectamine 2000 reagent (Invitrogen) according to the manufacturer’s instructions. siRNA against STAT3 (Ruibio, Guangzhou, China) and a non-targeting siRNA control (Ruibio, Guangzhou, China) were used to knock down gene expression. Briefly, the oligonucleotides and plasmids were transfected using Lipofectamine 2000 (Invitrogen) according the manufacturer’s protocol. 48 h later, G418 was added into the medium to select of stable clones. The pcDNA3.1-STAT3 and pcDNA3.1 vectors were purchased from Santa Cruz Biotechnology (USA). miR-298, miR-ctrl, anti-miR-ctrl and anti-miR-298 were purchased from Genechem (Shanghai, China).

The shRNA sequence that targeted LINC01287 was AAGCATTGTAGACCTGGCTGCTGAA. The sequences were cloned into the pGFP-C-shLenti vector according to the manufacturer’s instructions (Origene). Then, the viruses were packaged in 293 T cells according to a standard protocol. HCC cells were infected with virus particles plus 6 μg/ml Polybrene.

Quantitative real-time PCR (q-RT-PCR)

Total RNA was extracted from HCC cells using TRIzol reagent (Invitrogen) according to the manufacturer’s instructions. The RNA was reverse transcribed to cDNA, which was followed by real-time PCR analyses. The primers used in the study are listed in Table 1.

Table 1 Primers used in the study

To determine the sub-cellular distribution of LINC01287, nuclear and cytoplasmic fractions of cells were separated using the PARIS Kit (Life Technologies) according to the manufacturer’s instructions. Then, the RNA was extracted from both fractions. Subsequently, RT-PCR was performed to examine the expression ratios of specific RNA molecules between the nuclear and cytoplasmic fractions.

MTT, colony formation and cell cycle assays

For the MTT assay, HCC cells were seeded in 96-well plates. After 24 h, 5 mg/ml MTT was added to each well. The HCC cells were then incubated for another 6 h. Subsequently, DMSO was added to each well, and the absorbance was measured at 490 nm by spectrophotometry.

For the colony formation assay, HCC cells were seeded in 6-well culture plates. Two weeks later, the cells were fixed in paraformaldehyde, stained with hematoxylin solution and counted under a microscope.

For the cell cycle assay, cells were harvested from the culture dishes and washed three times in cold PBS. The cells were then fixed in 70% ice-cold ethanol at 4 °C overnight. Finally, the cells were incubated in propidium iodide (supplemented with RNase A). A FACScaliber flow cytometry system (BD Biosciences) was used to examine the DNA content of labeled cells.

Boyden assay

For the Boyden assay, HCC cells in serum-free DMEM were seeded in the upper chamber, which was inserted into a 24-well plate. RPMI-1640 supplemented with 10% FBS was added to each well. After 24 h, the HCC cells that remained on the upper surface of the membrane were removed, while the cells that migrated to the lower membrane were fixed in paraformaldehyde, stained with crystal violet and counted.

Western blot and immunofluorescence (IF) assays

In the Western blot assay, the protein samples were transferred onto a PVDF membrane, which was followed by an incubation overnight at 4 °C in a 1:500 dilution of primary antibodies. Subsequently, the membrane was incubated with HRP-conjugated rabbit or mouse secondary antibodies for 1 h at room temperature and were then developed using a chemiluminescence reagent.

For the IF assay, cells were plated on culture slides. After 24 h, when the cells had adhered to the slides, the cells were rinsed in phosphate-buffered saline (PBS) three times and were then fixed in ice-cold methanol-acetone for 10 min. Subsequently, the cells were blocked for 10 min in 5% BSA in PBS, followed by incubation with the primary antibodies in PBS for 1.5 h at room temperature. After three washes in PBS, the slides were incubated with the secondary antibodies for 40 min. After three additional washes, the slides were stained with 4-, 6-diamidino-2-phenylindole (DAPI) for 10 min and were examined using an Olympus confocal imaging system.

In situ hybridization assays

The in situ detection of lncRNA was performed as previously described [15].

Chromatin immunoprecipitation (ChIP) assay

We performed the ChIP assay using a ChIP assay kit (Millipore, catalog: 17–371) according to the manufacturer’s instructions. Briefly, the cells were fixed with 1% formaldehyde to covalently crosslink the proteins to DNA, after which the chromatin was harvested from the cells. Subsequently, the crosslinked DNA was sheared by sonication (sheared to 200–1000 base pairs in length) and was subjected to an immunoselection process. Then, PCR was performed to measure the enrichment of DNA fragments of the putative c-JUN-binding sites within the LINC01287 promoter.

Luciferase reporter assay

We cloned the full-length STAT3 cDNA (lacking the 3-UTR) into the eukaryotic expression vector pcDNA3.1 (Invitrogen). Subsequently, the 3’-UTR of STAT3 was amplified and cloned downstream of the firefly luciferase gene in the pGL3 vector (Promega); this vector was termed the wild type (WT) STAT3–3’-UTR. Using the GeneTailor™ Site-Directed Mutagenesis System (Invitrogen), we established site-directed mutagenesis of the miR-298 binding sites in the STAT3 3’-UTR. This vector was termed the mutant type (MUT) STAT3–3’-UTR. Subsequently, we co-transfected the HCC cells with the wt or mut STAT3–3’-UTR vector and the miR-298 mimic or inhibitor. Finally, we performed a luciferase assay using a dual Luciferase reporter assay system (Promega) 36 h after transfection. To perform LINC01287 promoter luciferase assays, HCC cells were seeded into 24-well plates and were co-transfected with plasmids that contained the LINC01287 promoter, the pRL-TK-Renilla plasmid (Promega, USA) and pcDNA.3-c-jun.

In vivo tumor growth and invasion assay

All procedures involving animals were approved by the Institutional Committee on Animal Care of Gaozhou People’s Hospital. For tumor growth study, five mice were included into each group. Sh-ctrl, sh-LINC01287 or sh-LINC01287/anti-miR-298 cells were injected subcutaneously into both flanks of nude mice. Four weeks after implantation, the xenografts were removed from the mice and weighed. The tumor volume was calculated according to the following formula: 4π/3 × (width/2)2 × (length/2). The invasion assay was performed as previously described [16]. After 3 weeks, nude mice were evaluated for lung colonization capacity.

Statistical analysis

SPSS 13.0 and Graph Pad Prism 5.0 software were used for the statistical analysis. The values are shown as the mean ± the standard error of the mean (S.E.M). Analyses of different groups were performed using one-way ANOVA or two-tailed Student’s t-test. P < 0.05 was considered statistically significant.

Results

LINC01287 expression was elevated in HCC tissues and cell lines

With the use of the online software program circlncRNAnet [17], we identified several lncRNAs that were significantly elevated in HCC tissues. Among these lncRNAs, LINC01287 was the most significantly different (Additional file 1: Figure S1A). Thus, we focused on LINC01287 for further analyses.

We further analyzed the expression level of LINC01287 in a larger cohort of 117 HCC patients. We revealed that the expression level of LINC01287 was up-regulated in HCC tissues compared with adjacent tissues (Fig. 1a). Interestingly, the expression level of LINC01287 was up-regulated in advanced-stage HCC (Fig. 1b) in comparison with early-stage HCC. In parallel, the expression level of LINC01287 was up-regulated in HCC cell lines compared with the normal human liver cell line LO2 (Fig. 1c).

Fig. 1
figure 1

The LINC01287 expression level was up-regulated in HCC cell lines and tissues. a LINC01287 expression level was significantly up-regulated in primary HCC tissues. b The expression level of LINC01287 was up-regulated in patients with advanced stage HCC. c LINC01287 expression was higher in HCC cells compared with a normal liver cell line. d and e High expression of LINC01287 was associated with a shorter overall survival and a shorter disease-free survival time of HCC patients (blue and green curves represent low and high expression of LINC01287, respectively). * represents P value < 0.05

We then questioned whether LINC01287 is indeed a non-coding RNA. An online bioinformatics analysis revealed that LINC01287 has no coding capability (http://cpc.cbi.pku.edu.cn/programs/run_cpc.jsp). The in vitro translation assay confirmed this finding (Additional file 1: Figure S1B). Subcellular fractionation and real-time PCR analysis were then performed to show that LINC01287 was primarily located in the cytoplasm (Additional file 1: Figure S1C).

High expression of LINC01287 was associated with unfavorable prognosis

The association between LINC01287 expression and the clinicopathological features of the patients is shown in Table 2. The median value of LINC01287 expression was selected as the cut-off value for the separation of HCC patients with a low level of expression (49/98, 50%) from patients with a high level of expression of LINC01287 (49/98, 50%). Although LINC01287 expression was not associated with parameters such as age (P = 0.311), gender (P = 0.285) or HBV infection (P = 0.223), high LINC01287 expression was significantly correlated with tumor size (P = 0.001), lymph node metastasis (P = 0.026) and late clinical stage (P = 0.026). Interestingly, a Kaplan-Meier analysis found that a high expression level of LINC01287 was correlated with a shorter overall survival time and disease-free survival time in patients with HCC (Fig. 1d and e, P < 0.05).

Table 2 Associations between lncRNA LINC01287 expression and patients’ clinicopathological features

Inhibition of LINC01287 decreased HCC cell proliferation and invasion

Since Huh7 and Bel7402 cells exhibited the highest expression of LINC01287, we selected them for further studies. To investigate the biological function of LINC01287 in HCC cells, we established Huh7 and Bel7402 cells in which LINC01287 was stably knocked down (sh-LINC01287) (Fig. 2a). It was revealed that LINC01287 inhibition significantly decreased cell proliferation, as determined by MTT assay (Fig. 2b). By colony formation assay, we found that oncogenic survival was significantly inhibited in sh-LINC01287 cells compared with sh-ctrl cells (Fig. 2c, Additional file 2: Figure S2A). Flow cytometry analysis showed that sh-LINC01287 cells presented a significantly higher frequency of cells in G1 phase and a lower frequency of cells in S phase (Fig. 2d, Additional file 2: Figure S2B). Interestingly, the Western blot assay revealed that the expression levels of G1/S phase checkpoint proteins such as cyclin D1, CDK4 and CDK6 were significantly decreased in sh-LINC01287 cells (Fig. 2e).

Fig. 2
figure 2

LINC01287 inhibition decreased HCC cell proliferation and invasion. a LINC01287 expression in Huh-7 and Bel7402 cells transduced with the control shRNA vector (sh-ctrl) or the LINC01287 shRNA vector (sh-LINC01287). b The MTT assay revealed that LINC01287 down-regulation significantly decreased cell proliferation. c The colony formation assay demonstrated that oncogenic survival was significantly decreased in sh-LINC01287 cells compared with sh-ctrl cells. d A higher frequency of sh-LINC01287 cells were in G1 phase and a lower frequency of these cells were in S phase. e LINC01287 down-regulation affected the expression of G1/S phase checkpoint proteins. f LINC01287 down-regulation decreased the invasion ability of HCC cells, as revealed by the Boyden assay. g The expression levels of MMP-2, MMP-3 and MMP-9 proteins were lower in sh-LINC01287 cells. h The expression level of E-cadherin was increased, while the expression levels of N-cadherin and Vimentin were decreased in sh-LINC01287 cells

Subsequently, we asked whether LINC01287 affected HCC cell invasiveness. The Boyden assay found that inhibition of LINC01287 decreased HCC cell invasiveness (Fig. 2f Additional file 2: Figure S2C). Consistently, according to the Western blot results, the expression levels of protein markers of invasion, including MMP-2, MMP-3 and MMP-9, were markedly decreased in sh-LINC01287 cells (Fig. 2g). Interestingly, the expression of epithelial to mesenchymal transition-related markers was also altered when LINC01287 was inhibited; E-cadherin expression was elevated, while N-cadherin and Vimentin expression was downregulated in sh-LINC01287 cells (Fig. 2h).

Overall, these data suggested that LINC01287 promotes HCC cell proliferation and invasion in vitro.

LncRNA LINC01287 functioned as a ceRNA to regulate miR-298 expression

LncRNA can exert its function by acting as a molecular sponge of miRNA [18]. We thus hypothesized that LINC01287 interacts with miRNAs in HCC. Using the online software miRDB (http://mirdb.org/miRDB/index.html), we searched for miRNAs with base pairs that are complementary to the LINC01287 sequence. We identified several microRNAs that may be regulated by LINC01287 (Additional file 2: Figure S2K). We first tested the difference in the expression of these microRNAs between the sh-ctrl and sh-LINC01287 groups of Huh7 cells. We found that only miR-298, miR-4308 and miR-23c were significantly altered between the sh-ctrl and sh-LINC01287 groups (Additional file 2: Figure S2L). In the sh-LINC01287 group, the level of miR-298 expression was increased approximately ten-fold higher compared with the sh-ctrl group, and consequently, we chose miR-298 for further study. The binding sites of miR-298 on LINC01287 are indicated in Fig. 3a. MiRNAs exert their gene silencing functions through a ribonucleoprotein complex, which is termed the RNA-induced silencing complex (RISC), the core component of which is Ago2 [19]. A radioimmunoprecipitation (RIP) assay was performed to examine whether LINC01287 and miR-298 are in the same RISC complex. It was revealed that LINC01287 and miR-298 were enriched in Ago2 immunoprecipitates compared with control IgG immunoprecipitates (Fig. 3b, P < 0.05). We then confirmed that when LINC01287 was down-regulated, the expression level of miR-298 was elevated in Huh7 and Bel7402 cells (Fig. 3c, P < 0.05). These data suggested that LINC01287 may negatively regulate miR-298. Subsequently, we examined whether miR-298 inhibited LINC01287 expression. The predicted miR-298 binding site of LINC01287 (LINC01287-WT) and a mutated miR-298 binding site of LINC01287 (LINC01287-MUT) were cloned into a reporter plasmid. We then performed luciferase assay to reveale whether miR-298 bound at LINC01287 and regulated LINC01287 expression. Co-transfection of miR-298 and LINC01287 decreased the luciferase activity, while co-transfection of miR-ctrl and LINC01287-WT and co-transfection of miR-298 and LINC01287-MUT did not change the luciferase activity (Fig. 3d, P < 0.05). In addition, we examined whether miR-298 regulated LINC01287 expression. It was revealed that, compared with the miR-ctrl, miR-298 inhibited LINC01287 expression. However, anti-miR-298 increased LINC01287 expression compared with anti-miR-ctrl (Fig. 3e, P < 0.05). Taken together, our data suggested a reciprocal interaction between LINC01287 and miR-298.

Fig. 3
figure 3

LncRNA LINC01287 acted as a ceRNA to regulate miR-298 expression. a The binding sites of miR-298 on LINC01287. b The RIP assay revealed that LINC01287 and miR-298 were enriched in the same Ago2 immunoprecipitates. c MiR-298 expression was increased in sh-LINC01287 cells compared with sh-ctrl cells. d Co-transfection of miR-298 and LINC01287-Wt strongly decreased the luciferase activity, while co-transfection of miR-ctrl and LINC01287-Wt did not change the luciferase activity. Co-transfection of miR-298 and LINC01287-Mut also did not change the luciferase activity. e MiR-298 decreased LINC01287 expression, while anti-miR-298 increased LINC01287 expression

LINC01287 regulated STAT3 expression through miR-298

With the use of the software Targetscan, we identified a list of genes that may be down-stream targets of miR-298. STAT3 is a key master regulator that controls cancer proliferation and invasion, and thus we selected this protein for further study. The 3’-UTR region of STAT3 mRNA, including the predicted miR-298 recognition site (wild type) or the mutated sequence (mutant type), were subcloned into luciferase reporter plasmids (Fig. 4a). It was observed that miR-298 decreased luciferase activity in the wild type vector instead of that in the mutant type (Fig. 4b). Moreover, miR-298 decreased the expression of both STAT3 mRNA and protein (Fig. 4c and d). Thus, our data revealed that STAT3 is a potential target of miR-298.

Fig. 4
figure 4

STAT3 was a down-stream target of miR-298. a The binding sites of miR-298 on STAT3. b The luciferase assay showed that cells transfected with miR-298 had less luciferase activity than those transfected with miR-ctrl. c miR-298 repressed STAT3 mRNA expression in HCC cells. d miR-298 repressed STAT3 protein expression in HCC cells. e Anti-miR-298 treatment led to the restoration of STAT3 expression in sh-LINC01287 cells. f The MTT assay revealed that sh-LINC01287 cells grew more slowly than sh-ctrl cells, while the overexpression of STAT3 rescued this effect. g The colony formation assay showed that sh-LINC01287 cells formed smaller and fewer colonies than the sh-ctrl cells, which was counteracted by the overexpression of STAT3. h LINC01287 down-regulation affected the cell cycle distribution, which was counteracted by the overexpression of STAT3. i LINC01287 down-regulation inhibited the invasion ability of HCC cells, which was rescued by the overexpression of STAT3

Subsequently, we wanted to know whether LINC01287 exerts its function through the miR-298/STAT3 axis. We observed that LINC01287 inhibition decreased STAT3 expression. However, when miR-298 was inhibited, the effect of LINC01287 down-regulation on STAT3 expression was abolished (Fig. 4e). The overexpression of STAT3 was able to counteract the effect of LINC01287 down-regulation on cell growth, cell cycle distribution and invasion (Fig. 4f-i, Additional file 2: Figure S2D-F). We also knocked down miR-298 expression in sh-LINC01287 cells. We revealed that knockdown of miR-298 could counteract LINC01287 down-regulation’s effect on HCC cells (Additional file 2: Figure S2G-J).

Taken together, these data suggested that LINC01287 could exert its function through the miR-298/STAT3 axis.

STAT3 elevated LINC01287 expression via c-Jun

We next explored the existence of a feedback loop between STAT3 and LINC01287. We used the UCSC and PROMO bioinformatics software programs to analyze an 800-bp region upstream of the transcriptional start site (TSS) of LINC01287. Four c-jun-binding motifs at − 79 to − 85, − 251 to − 257, − 364 to − 370 and − 674 to − 680 were identified inside the putative promoter region upstream of the LINC01287 TSS. We named these transcription factor-binding sites (TFBSs) A, B, C and D (Fig. 5a). A previous study revealed that STAT3 expression elevated c-jun expression [20], and we assumed that STAT3 may regulate LINC01287 via c-jun. We revealed that c-jun expression was down-regulated in si-STAT3 cells. However, forced overexpression of STAT3 elevated c-jun expression (Fig. 5b). These data suggested that STAT3 positively regulates c-jun expression.

Fig. 5
figure 5

STAT3 elevated LINC01287 expression via c-jun. a The putative transcription factor-binding sites of c-jun in the LINC01287 promoter region. b Western blot revealed that STAT3 positively regulated c-jun expression. c Knockdown of c-jun decreased the expression of miR-298 in both Huh7 and Bel7402 cells. d An up-regulation of wild-type miR-298 promoter luciferase activity was observed upon up-regulation of c-jun in Huh cells. e When c-jun binding sites on LINC01287 were mutated, the effect of c-jun on LINC01287 in Huh7 cells was abolished. f In si-STAT3 cells, the forced overexpression of c-jun increased LINC01287 expression. g In HCC cells in which STAT3 was overexpressed, the knockdown of c-jun decreased LINC01287 expression. h LINC01287 expression was positively correlated with STAT3 expression in HCC

To explore whether c-jun increased LINC01287 expression, c-jun was down-regulated by specific siRNAs. The knockdown of c-jun decreased the expression of LINC01287 in both Huh7 and Bel7402 cells (Fig. 5c). To further confirm that c-jun binds to the LINC01287 promoter, chromatin immunoprecipitation and q-PCR assays were performed. An up-regulation of wild-type LINC01287 promoter luciferase activity was observed upon up-regulation of c-jun in Huh7 cells (Fig. 5d). Interestingly, when c-jun binding sites on LINC01287 were mutated, we observed that the effect of c-jun on LINC01287 was abolished (Fig. 5e). Considering all these findings, these data suggested that c-jun binds to the LINC01287 promoter and increases its expression.

In si-STAT3 cells, forced overexpression of c-jun increased the LINC01287 expression level (Fig. 5f). However, when STAT3 was overexpressed, the knockdown of c-jun decreased LINC01287 expression (Fig. 5g). Taken together, these findings suggested that STAT3 regulates LINC01287 expression via c-jun.

We further revealed that LINC01287 expression was positively correlated with STAT3 expression (Fig. 5h, Spearman’s correlation, R = 0.5399).

Knockdown of lncRNA LINC01287 inhibited tumor growth in vivo

We then asked whether LINC01287 down-regulation inhibited tumor growth in vivo. When LINC01287 was inhibited, HCC cell-derived tumors grew more slowly. However, when miR-298 was inhibited in sh-LINC01287 cells, the LINC01287 down-regulation effect was partly abolished (Fig. 6a). The mean weight of the xenograft tumors was lower in the sh-LINC01287 group than in the sh-ctrl group and this effect was partly abolished by anti-miR-298 treatment (Fig. 6b). Tumor sections were stained for Ki-67 expression to quantitatively assess the proliferation index in the xenograft tumors; the proliferation index was lower in the sh-LINC01287 group and this effect was partly abolished by anti-miR-298 treatment (Fig. 6c). In addition, the knockdown of LINC01287 decreased the ability of these tumors to metastasize to the lungs in vivo and this effect was partly abolished by anti-miR-298 treatment (Fig. 6d). Taken together, these data supported the growth-promoting effect of LINC01287 and suggested that LINC01287 may exerted its function via miR-298 in vivo.

Fig. 6
figure 6

LINC01287 inhibition decreased tumor growth and invasion in vivo. a Compared with sh-ctrl cell-derived xenograft tumors, sh- LINC01287 cell-derived xenograft tumors grew more slowly, while anti-miR-298 treatment partly abolished this effect. b The mean weight of sh-LINC01287 cell-derived xenograft tumors was significantly lower than that of sh-ctrl cell-derived xenograft tumors, while anti-miR-298 treatment partly abolished this effect. c The knockdown of LINC01287 significantly decreased the percentage of Ki-67-positive cells in the tumors compared with the negative control group, while anti-miR-298 treatment partly abolished this effect. d The knockdown of LINC01287 decreased lung metastasis in vivo, while anti-miR-298 treatment partly abolished this effect

Discussion

LncRNAs regulate multiple biological processes, including cancer cell growth and invasion [21]. Recent studies have suggested that lncRNAs can serve as effective therapeutic targets for cancer treatment. Recently, a list of lncRNAs have been reported to be implicated in regulating HCC proliferation, migration and invasion [22,23,24]. These literatures suggested that lncRNAs may be a useful strategy for HCC treatment.

In the present study, we used the TGCA database to explore lncRNAs that may be involved in HCC progression. LINC01287 was identified as an oncogene whose expression is significantly increased in HCC tissues. Furthermore, we revealed that LINC01287 was up-regulated in HCC cell lines and tissues. The biological role of LINC01287 has not been investigated before. We then performed functional study to explore the role of LINC01287 in HCC cells. The knockdown of LINC01287 inhibited cell proliferation and colony formation in vitro. The in vivo study revealed that LINC01287 down-regulation decreased tumor growth. In terms of its mechanism, LINC01287 down-regulation contributed to cell cycle arrest in G1 stage. The G1/S phase checkpoint proteins (e.g., c-myc, cyclin D1 and CDK4) were altered when LINC01287 was inhibited. These data suggested that LINC01287 may lead to cell growth via changes in cell cycle progression. We also revealed that LINC01287 down-regulation inhibited GC cell invasion in vitro and decreased lung metastasis in vivo. EMT plays a vital role in promoting cancer cell invasion [25]. We thus asked whether LINC01287 is involved in the EMT phenotype. It was revealed that LINC01287 down-regulation increased the expression of the epithelial marker E-cadherin, while it decreased the expression of the mesenchymal markers N-cadherin and Vimentin. These data suggested that LINC01287 may promote an EMT phenotype and thus lead to HCC cell invasion.

A previous study revealed that lncRNAs may act as endogenous molecular sponges that compete for miRNAs and negatively regulate miRNA expression [26, 27]. The interaction between lncRNAs and microRNAs have well been documented in cancer. For instance, lncRNA-PAGBC competitively binds to the tumour suppressive microRNAs miR-133b and miR-511 [28]. LncRNA FAL1 promotes cell proliferation and migration by acting as a ceRNA of miR-1236 in HCC cells [29]. These literatures prompted us to ask whether there was interaction between LINC01287 and miRNAs. Using online software, we identified several miRNAs that may interact with LINC01287. Among these miRNAs, we selected miR-298 for further study, since the expression level of miR-298 was significantly increased in sh-LINC01287 cells. Previous documents revealed that miR-298 was implicated in regulating cancer progression. MiR-298 was frequently down-regulated in cancer tissues and acted as a tumor suppressor by inhibiting cell proliferation and migration [30,31,32]. The role of miR-298 has seldom been reported in HCC, but we revealed that miR-298 may be a tumor suppressor in HCC. We also confirmed the regulatory relationship between LINC01287 and miR-298 based on the following data: 1) LINC01287 down-regulation increased miR-298 expression; 2) the luciferase activity assay confirmed the direct binding ability of the predicted miR-298 binding site on LINC01287; 3) The RIP assays found that LINC01287 and miR-298 were in the same RISC.

Emerging evidence has demonstrated the role of STAT3 in cancer progression, invasion and metastasis [33]. STAT3 was identified as a down-stream target of miR-298 in our study, and our data revealed that LINC01287 inhibition decreased STAT3 expression. However, when miR-298 was inhibited, the effect of LINC01287 down-regulation on STAT3 expression was abolished. The overexpression of STAT3 can counteract the effect of LINC01287 on HCC cells. These data suggested that LINC01287 exerts its function through the miR-298/STAT3 axis.

We further sought to determine the presence of a positive feedback loop between STAT3 and LINC01287. The recruitment of specific transcription factors often leads to abnormal lncRNA expression [34]. The transcription factor c-jun is a key regulator of cell growth [35] and metastasis [36] in cancer. We revealed four putative binding sites of c-jun in the region upstream of the LINC01287 locus. A subsequent experiment demonstrated that c-jun could positively regulate LINC01287 expression by directly binding to its promoter. Our study showed that STAT3 increased c-jun expression and therefore regulated LINC01287 expression. We further confirmed a positive correlation between STAT3 and LINC01287 in HCC tissues. Taken together, our findings revealed a feedback loop within the LINC01287/miR-298/STAT3 axis.

Overall, our data provided the first evidence that the LINC01287/miR-298/STAT3 axis controls cell growth and invasiveness of HCC cells. Therapeutics that target LINC01287 may therefore improve the treatment of HCC.

Conclusions

In conclusion, our study demonstrated that LINC01287 was upregulated in HCC and may function as a ceRNA to increase STAT3 expression by sponging miR-298, which consequently contributes to HCC growth and metastasis. Our findings indicated that LINC01287 may be a potential therapeutic target for HCC treatment.

References

  1. Marquardt JU, Thorgeirsson SS. SnapShot: Hepatocellular carcinoma. Cancer Cell. 2014;25(4):550–1.

    Article  PubMed  CAS  Google Scholar 

  2. Yuan JH, Yang F, Wang F, Ma JZ, Guo YJ, Tao QF, Liu F, Pan W, Wang TT, Zhou CC, et al. A long noncoding RNA activated by TGF-beta promotes the invasion-metastasis cascade in hepatocellular carcinoma. Cancer Cell. 2014;25(5):666–81.

    Article  PubMed  CAS  Google Scholar 

  3. Ma MZ, Chu BF, Zhang Y, Weng MZ, Qin YY, Gong W, Quan ZW. Long non-coding RNA CCAT1 promotes gallbladder cancer development via negative modulation of miRNA-218-5p. Cell Death Dis. 2015;6:e1583.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  4. Jie Sun ZG, Zhang B, Deng X, Ma L, Zhengtao G, Yang H. MiR-137 inhibits proliferation and angiogenesis of human glioblastoma cells by targeting EZH2. J Neuro-Oncol. 2014;122(3):481–9.

    Article  CAS  Google Scholar 

  5. Cheetham SW, Gruhl F, Mattick JS, Dinger ME. Long noncoding RNAs and the genetics of cancer. Br J Cancer. 2013;108(12):2419–25.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  6. Lu S, Zhou J, Sun Y, Li N, Miao M, Jiao B, Chen H. The noncoding RNA HOXD-AS1 is a critical regulator of the metastasis and apoptosis phenotype in human hepatocellular carcinoma. Mol Cancer. 2017;16(1):125.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Huang JL, Cao SW, Ou QS, Yang B, Zheng SH, Tang J, Chen J, Hu YW, Zheng L, Wang Q. The long non-coding RNA PTTG3P promotes cell growth and metastasis via up-regulating PTTG1 and activating PI3K/AKT signaling in hepatocellular carcinoma. Mol Cancer. 2018;17(1):93.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Lv J, Kong Y, Gao Z, Liu Y, Zhu P, Yu Z. LncRNA TUG1 interacting with miR-144 contributes to proliferation, migration and tumorigenesis through activating the JAK2/STAT3 pathway in hepatocellular carcinoma. Int J Biochem Cell Biol. 2018;101:19–28.

    Article  PubMed  CAS  Google Scholar 

  9. Tu J, Zhao Z, Xu M, Lu X, Chang L, Ji J. NEAT1 upregulates TGF-beta1 to induce hepatocellular carcinoma progression by sponging hsa-mir-139-5p. J Cell Physiol. 2018.

  10. Deng X, Ma L, Wu M, Zhang G, Jin C, Guo Y, Liu R. miR-124 radiosensitizes human glioma cells by targeting CDK4. J Neuro-Oncol. 2013;114(3):263–74.

    Article  CAS  Google Scholar 

  11. Pei ZJ, Zhang ZG, Hu AX, Yang F, Gai Y. miR-122-5p inhibits tumor cell proliferation and induces apoptosis by targeting MYC in gastric cancer cells. Die Pharmazie-Int J Pharm Sci. 2017;72(6):344–7.

    CAS  Google Scholar 

  12. Dasgupta P, Kulkarni P, Majid S, Shahryari V, Hashimoto Y, Bhat NS, et al. MicroRNA-203 inhibits long noncoding RNA HOTAIR and regulates tumorigenesis through epithelial-to-mesenchymal transition pathway in renal cell carcinoma. Mol Cancer Ther. 2018;17:1061–9.

    Article  PubMed  CAS  Google Scholar 

  13. Shang M, Wang X, Zhang Y, Gao Z, Wang T, Liu R. LincRNA-ROR promotes metastasis and invasion of esophageal squamous cell carcinoma by regulating miR-145/FSCN1. OncoTargets Ther. 2018;11:639–49.

    Article  Google Scholar 

  14. Xie CR, Wang F, Zhang S, Wang FQ, Zheng S, Li Z, Lv J, Qi HQ, Fang QL, Wang XM, et al. Long noncoding RNA HCAL facilitates the growth and metastasis of hepatocellular carcinoma by acting as a ceRNA of LAPTM4B. Mol Ther Nucleic Acids. 2017;9:440–51.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  15. Liu SJ, Nowakowski TJ, Pollen AA, Lui JH, Horlbeck MA, Attenello FJ, He D, Weissman JS, Kriegstein AR, Diaz AA, et al. Single-cell analysis of long non-coding RNAs in the developing human neocortex. Genome Biol. 2016;17:67.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  16. Zhao X, He L, Li T, Lu Y, Miao Y, Liang S, Guo H, Bai M, Xie H, Luo G, et al. SRF expedites metastasis and modulates the epithelial to mesenchymal transition by regulating miR-199a-5p expression in human gastric cancer. Cell Death Differ. 2014;21(12):1900–13.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  17. Wu SM, Liu H, Huang PJ, Chang IY, Lee CC, Yang CY, Tsai WS, Tan BC. circlncRNAnet: an integrated web-based resource for mapping functional networks of long or circular forms of noncoding RNAs. GigaScience. 2018;7(1):1–10.

    CAS  Google Scholar 

  18. Cheng Z, Li Z, Ma K, Li X, Tian N, Duan J, Xiao X, Wang Y. Long non-coding RNA XIST promotes glioma Tumorigenicity and angiogenesis by acting as a molecular sponge of miR-429. J Cancer. 2017;8(19):4106–16.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Grigoriev A, Bonini NM. Age-dependent patterns of microRNA RISC loading. Aging. 2014;6(9):705–6.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  20. Yang E, Lerner L, Besser D, Darnell JE Jr. Independent and cooperative activation of chromosomal c-fos promoter by STAT3. J Biol Chem. 2003;278(18):15794–9.

    Article  PubMed  CAS  Google Scholar 

  21. Chen DL, Ju HQ, Lu YX, Chen LZ, Zeng ZL, Zhang DS, Luo HY, Wang F, Qiu MZ, Wang DS, et al. Long non-coding RNA XIST regulates gastric cancer progression by acting as a molecular sponge of miR-101 to modulate EZH2 expression. J Exp Clin Cancer Res. 2016;35(1):142.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  22. Xiao JN, Yan TH, Yu RM, Gao Y, Zeng WL, Lu SW, Que HX, Liu ZP, Jiang JH. Long non-coding RNA UCA1 regulates the expression of Snail2 by miR-203 to promote hepatocellular carcinoma progression. J Cancer Res Clin Oncol. 2017;143(6):981–90.

    Article  PubMed  CAS  Google Scholar 

  23. Li T, Xie J, Shen C, Cheng D, Shi Y, Wu Z, Deng X, Chen H, Shen B, Peng C, et al. Upregulation of long noncoding RNA ZEB1-AS1 promotes tumor metastasis and predicts poor prognosis in hepatocellular carcinoma. Oncogene. 2016;35(12):1575–84.

    Article  PubMed  CAS  Google Scholar 

  24. Zhang JY, Weng MZ, Song FB, Xu YG, Liu Q, Wu JY, Qin J, Jin T, Xu JM. Long noncoding RNA AFAP1-AS1 indicates a poor prognosis of hepatocellular carcinoma and promotes cell proliferation and invasion via upregulation of the RhoA/Rac2 signaling. Int J Oncol. 2016;48(4):1590–8.

    Article  PubMed  CAS  Google Scholar 

  25. Weng J, Xiao J, Mi Y, Fang X, Sun Y, Li S, Qin Z, Li X, Liu T, Zhao S, et al. PCDHGA9 acts as a tumor suppressor to induce tumor cell apoptosis and autophagy and inhibit the EMT process in human gastric cancer. Cell Death Dis. 2018;9(2):27.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  26. Xiao Y, Jiao C, Lin Y, Chen M, Zhang J, Wang J, Zhang Z. lncRNA UCA1 contributes to Imatinib resistance by acting as a ceRNA against miR-16 in chronic myeloid leukemia cells. DNA Cell Biol. 2017;36(1):18–25.

    Article  PubMed  CAS  Google Scholar 

  27. Liu Q, Guo X, Que S, Yang X, Fan H, Liu M, Li X, Tang H. LncRNA RSU1P2 contributes to tumorigenesis by acting as a ceRNA against let-7a in cervical cancer cells. Oncotarget. 2017;8(27):43768–81.

    PubMed  Google Scholar 

  28. Wu XS, Wang F, Li HF, Hu YP, Jiang L, Zhang F, Li ML, Wang XA, Jin YP, Zhang YJ, et al. LncRNA-PAGBC acts as a microRNA sponge and promotes gallbladder tumorigenesis. EMBO Rep. 2017;18(10):1837–53.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  29. Li B, Mao R, Liu C, Zhang W, Tang Y, Guo Z. LncRNA FAL1 promotes cell proliferation and migration by acting as a CeRNA of miR-1236 in hepatocellular carcinoma cells. Life Sci. 2018;197:122–9.

    Article  PubMed  CAS  Google Scholar 

  30. Zhou F, Chen J, Wang H. MicroRNA-298 inhibits malignant phenotypes of epithelial ovarian cancer by regulating the expression of EZH2. Oncol Lett. 2016;12(5):3926–32.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  31. Zhao H, Zhao D, Tan G, Liu Y, Zhuang L, Liu T. Bufalin promotes apoptosis of gastric cancer by down-regulation of miR-298 targeting bax. Int J Clin Exp Med. 2015;8(3):3420–8.

    PubMed  PubMed Central  CAS  Google Scholar 

  32. Bao L, Hazari S, Mehra S, Kaushal D, Moroz K, Dash S. Increased expression of P-glycoprotein and doxorubicin chemoresistance of metastatic breast cancer is regulated by miR-298. Am J Pathol. 2012;180(6):2490–503.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  33. Xue J, Liao L, Yin F, Kuang H, Zhou X, Wang Y. LncRNA AB073614 induces epithelial-mesenchymal transition of colorectal cancer cells via regulating the JAK/STAT3 pathway. Cancer Biomark. 2018;21:849–58.

    Article  PubMed  CAS  Google Scholar 

  34. Cao L, Zhang P, Li J, Wu M. LAST, a c-Myc-inducible long noncoding RNA, cooperates with CNBP to promote CCND1 mRNA stability in human cells. eLife. 2017;6:e30433.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Chakraborty A, Diefenbacher ME, Mylona A, Kassel O, Behrens A. The E3 ubiquitin ligase Trim7 mediates c-Jun/AP-1 activation by Ras signalling. Nat Commun. 2015;6:6782.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  36. Lamb RF, Hennigan RF, Turnbull K, Katsanakis KD, MacKenzie ED, Birnie GD, Ozanne BW. AP-1-mediated invasion requires increased expression of the hyaluronan receptor CD44. Mol Cell Biol. 1997;17(2):963–76.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

Download references

Acknowledgements

We thank Dr. Xiaoming Lv, of Sun Yat-sen University, who provided the pcDNA.3-STAT3 and pcDNA.3-c-jun plasmids.

Funding

None declared.

Availability of data and materials

Data sharing not applicable to this article as no datasets were generated or analyzed during the current study.

Author information

Authors and Affiliations

Authors

Contributions

Study concept and design: SC and YM; acquisition of data: LH, ZL and ZW; analysis and interpretation of data: QC, SP, LL and DL; drafting of the manuscript: JH and FX. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Siyao Che.

Ethics declarations

Ethics approval and consent to participate

The study on HCC cancer samples was approved and supervised by the Research Ethics Committee of Gaozhou people’s Hospital. Written Informed Consents were obtained from all patients. The animal experiments were performed in strict accordance with the guidelines of the Research Animal Care Committee Gaozhou people’s Hospital.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Figure S1. (A) The TGCA database revealed that LINC01287 was up-regulated in HCC tissues. (B) The in vitro translation assay revealed that LINC01287 did not have coding ability. (C) LINC01287 was primarily expressed in the cytoplasm, as determined by RT-PCR. (TIF 4712 kb)

Additional file 2:

Figure S2. (A) The colony formation assay demonstrated that oncogenic survival was significantly decreased in sh-LINC01287 cells compared with sh-ctrl cells. (B) sh-LINC01287 cells were significantly more likely to be in G1 phase and were less likely to be in S phase. (C) LINC01287 down-regulation decreased the invasion ability of HCC cells, as revealed by the Boyden assay. (D) The colony formation assay showed that sh-LINC01287 cells formed smaller and fewer colonies than the sh-ctrl cells, which was counteracted by the overexpression of STAT3. (E) LINC01287 down-regulation affected the cell cycle distribution, which was counteracted by the overexpression of STAT3. (F) LINC01287 down-regulation inhibited the invasion ability of HCC cells, which was rescued by the overexpression of STAT3. (G) The MTT assay revealed that LINC01287 down-regulation significantly decreased cell proliferation and the effect was counteracted by anti-miR-298 treatment. (H) The colony formation assay showed that sh-LINC01287 cells formed smaller and fewer colonies than the sh-ctrl cells, which was counteracted by anti-miR-298 treatment. (I) LINC01287 down-regulation affected the cell cycle distribution, which was counteracted by anti-miR-298 treatment. (J) LINC01287 down-regulation inhibited the invasion ability of HCC cells, which was rescued by anti-miR-298 treatment. (K) The putative microRNAs that may be regulated by LINC01287. (L) The expression levels of miR-298, miR-4308 and miR-23c were increased in sh-LINC01287 cells. (TIF 18965 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Mo, Y., He, L., Lai, Z. et al. LINC01287/miR-298/STAT3 feedback loop regulates growth and the epithelial-to-mesenchymal transition phenotype in hepatocellular carcinoma cells. J Exp Clin Cancer Res 37, 149 (2018). https://doi.org/10.1186/s13046-018-0831-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13046-018-0831-2

Keywords