Long noncoding RNA HEGBC promotes tumorigenesis and metastasis of gallbladder cancer via forming a positive feedback loop with IL-11/STAT3 signaling pathway

Background Gallbladder cancer (GBC) is a highly malignant cancer with poor prognosis. Several long noncoding RNAs (lncRNAs) have been reported to be involved in the tumorigenesis and progression of GBC. However, the expressions, clinical significances, and roles of most other lncRNAs in GBC are still unknown. Methods The differentially expressed lncRNAs in GBC were screened through re-analyzing the public available microarray datasets. The expression of lncRNA high expressed in gallbladder cancer (lncRNA-HEGBC) in GBC was measured by qRT-PCR. The correlations between HEGBC with clinicopathological characteristics and prognosis were analyzed by Pearson chi-square test and log-rank test. A series of in vitro and in vivo, gain-of and loss-of function assays were performed to investigate the roles of HEGBC in GBC cell proliferation, apoptosis, migration, tumor growth and metastasis. The interactions between HEGBC and IL-11/STAT3 signaling were explored using chromatin isolation by RNA purification (ChIRP), chromatin immunoprecipitation (ChIP), enzyme linked immunosorbent assay (ELISA), qRT-PCR, western blot, and luciferase reporter assays. Results We identified a novel lncRNA HEGBC, which is upregulated in GBC and positively associated with advanced TNM stages and poor prognosis of GBC patients. Overexpression of HEGBC increased GBC cell viability, inhibited GBC cell apoptosis, promoted GBC cell migration, and promoted GBC tumor growth and metastasis in vivo. Conversely, depletion of HEGBC decreased GBC cell viability, promoted GBC cell apoptosis, inhibited GBC cell migration, and inhibited GBC tumor growth and metastasis in vivo. Mechanistic investigations showed that HEGBC bound to the promoter of IL-11, increased IL-11 transcription, induced IL-11 autocrine, and activated IL-11/STAT3 signaling pathway. Furthermore, STAT3 also bound to the promoter of HEGBC and activated HEGBC expression. Thus, HEGBC/IL-11/STAT3 formed a positive regulatory loop in GBC. Depletion of IL-11 attenuated the oncogenic roles of HEGBC in GBC. Conclusions Our findings identified a novel lncRNA HEGBC, which is upregulated and indicts poor prognosis of GBC. HEGBC exerts oncogenic roles in GBC via forming a positive regulatory loop with IL-11/STAT3 signaling. Our data suggested that HEGBC could be a potential prognostic biomarker and therapeutic target for GBC. Electronic supplementary material The online version of this article (10.1186/s13046-018-0847-7) contains supplementary material, which is available to authorized users.


Background
Gallbladder cancer (GBC) is the most common biliary tract malignant cancer, which is highly lethal and has an extremely poor prognosis [1,2]. Due to its non-specific symptoms and highly invasive property, most GBC patients are diagnosed at advanced stages [3]. Hence, most GBC patients at advanced stages are not candidates for surgical resection [4]. Unfortunately, until now the only curative treatment for GBC is still surgery [5]. Therefore, the mean survival time for GBC ranges from 13.2 months to 19 months [6]. The poor outcome of GBC and the lack of efficient therapies request better understanding of molecular mechanisms underlying GBC tumorigenesis and metastasis, and developing more efficient targeted therapies for GBC.
Although many aberrantly expressed and mutated molecular events have been identified in GBC, most of these focused on protein-coding genes [7]. Recently, with the great progressions of genome and transcriptome sequencing, many non-protein-coding genes have been identified, which accounts for about 70% of the genome, with only 2% of the genome encoding proteins [8]. Most of these non-protein-coding genes transcribe long non-coding RNAs (lncRNAs) [9]. lncRNA is a class of RNA with limited protein coding potential and more than 200 nucleotides in length [10]. Accumulating evidences revealed that lncRNAs are frequently involved in many pathophysiological processes, including cancers [11][12][13][14][15]. Many lncRNAs are dysregulated in cancers [16][17][18][19]. Furthermore, many lncRNAs are revealed to control cell proliferation, cell cycle, cell apoptosis, senescence, cell migration, cell invasion, drug resistance of cancer cells and so on [20][21][22][23]. For example, lncRNA-PAGBC is reported to be up-regulated in GBC and promote GBC tumorigenesis via competitively binding miR-133b and miR-511 [24]. LncRNA-CCAT1 is reported to be up-regulated in GBC and promote GBC development via negative regulating miR-218-5p [25]. LncRNA HOXA-AS2 is also up-regulated in GBC and promotes GBC proliferation and epithelial-mesenchymal transition [26]. LncRNA GCASPC is reported to be down-regulated in GBC and inhibit pyruvate carboxylase-dependent cell proliferation of GBC cells [27]. LncRNA H19 is reported to be up-regulated and have oncogenic roles in GBC via modulating miR-342-3p and FOXM1 [28]. In our previous study, we investigated the expression and roles of lncRNA SPRY4-IT1 in GBC and found that lncRNA SPRY4-IT1 is upregulated in GBC and promotes GBC cell proliferation, migration, and invasion [29]. Although several lncRNAs have been reported to be involved in GBC, the expression and roles of most of lncRNAs in GBC are still unclear.
In this study, we searched differentially expressed lncRNAs in GBC via analyzing public available microarray data and identified a novelly differently expressed lncRNA in GBC. We further analyzed the expression, clinical significances, roles, and mechanisms of action of this lncRNA in GBC.

Clinical specimens
A total of 102 pairs of GBC tissues and adjacent non-tumor tissues were obtained from GBC patients with written informed consent who underwent surgery at Eastern Hepatobiliary Surgery Hospital (Shanghai, China). All these GBC patients did not receive any pre-operative treatments. The tissue specimens were confirmed by histopathological diagnosis. All resected specimens were immediately snap-frozen in liquid nitrogen and stored at − 80°C until RNA extraction. The Review Board of Eastern Hepatobiliary Surgery Hospital reviewed and approved this study.

Cell lines and treatments
The human non-tumorigenic biliary epithelial cell line H69 and GBC cell lines SGC-996, NOZ, GBC-SD, and EH-GB2 were obtained from the Institute of Biochemistry and Cell Biology of the Chinese Academy of Sciences (Shanghai, China) or maintained in our hospital [30]. The cells were maintained in Dulbecco's Modified Eagle's Medium (Gibco BRL, Grand Island, NY, USA) supplemented with 10% fetal bovine serum (Gibco) in a humidified incubator containing 5% CO 2 at 37°C. Where indicated, the GBC cells were treated with 5 ng/mL doxorubicin (Selleck, Houston, TX, USA) for 24 h, 5 μM p-STAT3 inhibitor SC144 (Selleck) for 72 h, or 20 ng/mL IL-11 (Gibco) for 72 h.
For construction of HEGBC stably overexpressed GBC cells, pcDNA3.1-HEGBC or pcDNA3.1 was transfected into SGC-996 and NOZ cells, and then the cells were selected with 800 μg/mL neomycin for 4 weeks. For construction of HEGBC stably depleted GBC cells, shHEGBC-1, shHEGBC-2, or shControl was transfected into GBC-SD and EH-GB2 cells, and then the cells were selected with 1 μg/mL puromycin for 4 weeks.

Cell proliferation, apoptosis, and migration assays
Cell proliferation was detected with Glo cell viability assay and Ethynyl deoxyuridine (EdU) incorporation assay. For Glo cell viability assay, 2000 indicated GBC cells were plated into 96-well plates and cultured for indicated time. At each indicated time, the luminescence values were detected using the CellTiter-Glo® Luminescent Cell Viability Assay (Promega, Madison, WI, USA) in accordance with the manufacturer's instruction. EdU incorporation assay was performed with the EdU kit (Roche) in accordance with the manufacturer's instruction. Results were acquired using the Zeiss fluorescence photomicroscope (Carl Zeiss, Oberkochen, Germany) and quantified via counting at least five random fields. Cell apoptosis was detected with terminal deoxynucleotidyl transferase (TdT)-mediated dUTP nick end labelling (TUNEL) assay. After being treated with 5 ng/mL doxorubicin for 24 h, indicated GBC cells were used to perform TUNEL assay with the In-Situ Cell Death Detection Kit (Roche) in accordance with the manufacturer's instruction. Results were acquired using the Zeiss fluorescence photomicroscope (Carl Zeiss) and quantified via counting at least five random fields. Cell migration was detected with transwell assay. For transwell assay, 40,000 indicated GBC cells re-suspended in serumfree media with 1 μg/ml Mitomycin C to inhibit cell proliferation were placed into the upper chamber of transwell insert (8-μm pore size; Millipore, Bedford, MA, USA). Medium containing 10% FBS was added to the lower chamber. After incubation for 48 h, the GBC cells remaining on the upper membrane were removed with cotton wool, and whereas the cells migrating through the membrane were fixed with methanol, stained with 0.1% crystal violet, imaged using the Zeiss fluorescence photomicroscope (Carl Zeiss), and quantified via counting at least five random fields.

In vivo tumorigenesis and metastasis assays
To establish in vivo tumorigenesis model, 2 × 10 6 indicted GBC cells in 50 μL phosphate buffered saline mixed with 50% matrigel (Invitrogen) were subcutaneously injected into the flanks of 6 weeks old male athymic nude mice (SLRC Laboratory Animal Center, Shanghai, China). Subcutaneous tumor volumes were detected every 7 days with caliper for 28 days and calculated as a × b 2 × 0.5 (a, longest diameter; b, shortest diameter). At the 28th day after injection, the mice were sacrificed and subcutaneous tumors were resected and weighted. The resected subcutaneous tumors were fixed in formalin, paraffin embedded, deparaffinized, rehydrated, and antigen retrieved. For immunohistochemical staining of Ki67, the sections were incubated with primary antibody for Ki67 (Cell Signaling Technology, Boston, MA, USA), followed by being incubated with a horseradish peroxidase-conjugated secondary antibody (Cell Signaling Technology). Finally, the slides were visualized with 3, 3-diaminobenzidine. For detection of cell apoptosis of subcutaneous tumors, the sections were used to perform TUNEL assay with the In-Situ Cell Death Detection Kit (Roche) in accordance with the manufacturer's instruction. Liver metastasis model was established with intra-splenic injection of 2× 10 6 indicated GBC cells. Mice were allowed to grow for 6 weeks, and then the mice were sacrificed, and the livers were resected. The liver metastatic foci number was counted via HE staining. The Review Board of Eastern Hepatobiliary Surgery Hospital reviewed and approved the use of animals.

Isolation of cytoplasmic and nuclear RNA
Cytoplasmic and nuclear RNA were isolated and purified using the Cytoplasmic & Nuclear RNA Purification Kit (Norgen, Belmont, CA, USA) in accordance with the manufacturer's instruction. The isolated RNA was detected by qRT-PCR as described above.

RNA immunoprecipitation (RIP) assay
RIP assay was performed with the Magna RIP RNA-Binding Protein Immunoprecipitation Kit (Millipore) and p-STAT3 antibody (5 μg per reaction; Cell Signaling Technology), STAT3 antibody (5 μg per reaction; Cell Signaling Technology), RPLP0 antibody (5 μg per reaction; Abcam, Hong Kong, China), or negative control IgG in accordance with the manufacturer's instruction. RIP-derived RNA was quantified using qRT-PCR to detect enrichment of lncRNAs.

Western blot analysis
Total proteins were extracted from indicated GBC cells using RIPA buffer (Beyotime, Shanghai, China) and separated by 10% sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE), followed by being transferred to NC membrane. After being blocked with 5% bovine serum albumin, the membranes were incubated with primary antibodies against p-STAT3 (Cell Signaling Technology), STAT3 (Cell Signaling Technology), or β-actin (Sigma-Aldrich, Saint Louis, MO, USA). After being washed, the membranes were incubated with IRDye 800CW goat anti-rabbit IgG or IRDye 700CW goat anti-mouse IgG (Li-Cor, Lincoln, NE, USA), and detected using Odyssey infrared scanner (Li-Cor).

Luciferase reporter assays
The promoter of HEGBC containing the predicted p-STAT3 binding sites was PCR amplified using Thermo Scientific Phusion Flash High-Fidelity PCR Master Mix (Thermo-Fisher Scientific) and subcloned into the Kpn I and Xho I sites of the pGL3-basic vector (Promega), termed as pGL3-HEGBC-pro. The sequences of the primers were as follows: 5'-GGGGTACCCTATTGCTGCACTCACACA CCC-3′ (forward) and 5'-CCGCTCGAGCGCCAGAGC CCAAGCTATC-3′ (reverse). The empty vector pGL3basic was used as negative control. The p-STAT3 binding sites mutated HEGBC promoter was synthesized by Gen-Script (Nanjing, China) and subcloned into the Kpn I and Xho I sites of the pGL3-basic vector, termed as pGL3-HEGBC-pro-mut. The constructed luciferase reporter plasmids were cotransfected with the pRL-TK plasmid expressing renilla luciferase into NOZ cells. 12 h later after transfection, the NOZ cells were treated with 5 μM SC144 or 20 ng/mL IL-11 for 72 h. Then the luciferase activity was measured using Dual-Luciferase® Reporter Assay System (Promega) in accordance with the manufacturer's instruction.

Statistical analysis
Statistical analyses were carried out using the GraphPad Prism Software. For comparisons, Wilcoxon signed-rank test, Mann-Whitney test, Pearson chi-square test, Log-rank test, Student's t test, or Pearson's correlation analysis was performed as indicated. P < 0.05 was regarded as statistically significant.

Results
A newly identified lncRNA HEGBC was increased in GBC and associated with poor survival of GBC patients To search the lncRNAs involved in GBC, we analyzed the expression of lncRNAs in GBC using public available dataset from the National Center for Biotechnology Information (NCBI) Gene Expression Omnibus (GEO) (https:// www.ncbi.nlm.nih.gov/geo/). We noted a microarray data comparing the expression of lncRNAs and mRNAs in nine GBC tissues and nine corresponding non-tumor tissues with the GEO accession number GSE76633. We re-analyzed the microarray results and the 60 most differentially expressed lncRNAs and mRNAs are shown in Additional file 1: Table S1. Among the differentially expressed lncRNAs, several lncRNAs have been reported. A new lncRNA ENST00000414772 (LOC401585), which was 5.69-fold higher in GBC tissues than corresponding nontumor tissues, caught our attention and was named as lncRNA high expressed in gallbladder cancer (lncRNA--HEGBC). The HEGBC gene has one transcript in the NCBI database (https://www.ncbi.nlm.nih.gov/; NCBI Reference Sequence: NR_125365.1). The full-length of HEGBC was confirmed using rapid amplification of cDNA ends (RACE) assays (Additional file 2: Figure S1). The Open Reading Frame Finder from NCBI (https://www.ncbi.nlm.nih.gov/ orffinder/) failed to predict a protein of more than 19 amino acids. The Coding Potential Calculator (http:// cpc.cbi.pku.edu.cn/programs/run_cpc.jsp) also predic ted HEGBC as a noncoding RNA with coding potential score − 1.26. Furthermore, RIP assays showed no interaction between HEGBC and ribosomal protein RPLP0, a component of the 60S subunit of ribosome (Additional file 2: Figure S2). These data suggested that HEGBC had no protein-coding capability.
To confirm the expression pattern of HEGBC in GBC, qRT-PCR was carried out on 102 pairs of GBC tissues and adjacent non-tumor tissues. The results showed that HEGBC was significantly increased in GBC tissues compared with adjacent non-tumor tissues (Fig. 1a). HEGBC was more highly expressed in GBC extending beyond the gallbladder (T3 + T4) than that in GBC only detected in the gallbladder (T1 + T2) (Fig. 1b). HEGBC was also more highly expressed in GBC with lymph nodes metastasis (N1 + N2) than that in GBC without lymph nodes metastasis (N0) (Fig. 1c). Furthermore, HEGBC was more highly expressed in GBC with advanced clinical stages (III-IV) than that in GBC with early clinical stages (I-II) (Fig. 1d). Analyzing of the correlation between the expression of HEGBC and clinicopathologic characteristics of these 102 GBC patients indicated that high HEGBC expression is positively correlated with extending extent (T), lymph node metastasis (N), and TNM stages ( Table 1). Kaplan-Meier survival analysis showed that GBC patients with high HEGBC expression had worse survival than those with low HEGBC expression (Fig. 1e). Furthermore, in GBC with clinical stages I-II or III-IV, high HEGBC expression also indicated poor survival (Fig. 1f, g).
The expression of HEGBC was measured in human non-tumorigenic biliary epithelial cell line H69 and GBC cell lines SGC-996, NOZ, GBC-SD, and EH-GB2 by qRT-PCR. The results showed that HEGBC was also markedly increased in GBC cell lines compared with that in biliary epithelial cell line (Fig. 1h). Collectively, these results showed that HEGBC was upregulated in GBC and associated with poor survival of GBC patients.

Ectopic expression of HEGBC promoted the proliferation and migration of GBC cells
To investigate the biological roles of HEGBC in GBC, we constructed HEGBC stably overexpressed SGC-996 and NOZ cells via transfecting HEGBC overexpression plasmid pcDNA3.1-HEGBC (Fig. 2a, b). Glo cell viability assays indicated that ectopic expression of HEGBC significantly promoted cell proliferation of SGC-996 and NOZ cells (Fig. 2c, d). EdU incorporation assays also indicated that ectopic expression of HEGBC markedly promoted cell proliferation of SGC-996 and NOZ cells (Fig. 2e). TUNEL assays indicated that ectopic expression of HEGBC markedly inhibited cell apoptosis of SGC-996 and NOZ cells (Fig. 2f ). Transwell assays indicated that ectopic expression of HEGBC significantly promoted cell migration of SGC-996 and NOZ cells (Fig. 2g). Collectively, these results showed that ectopic expression of HEGBC promoted the proliferation and migration of GBC cells in vitro.

Depletion of HEGBC inhibited the proliferation and migration of GBC cells
To completely elucidate the biological roles of HEGBC on GBC cell proliferation and migration, HEGBC stably depleted GBC-SD and EH-GB2 cells were constructed via transfecting two dependent HEGBC specific shRNAs (Fig. 3a, b). Glo cell viability assays indicated that depletion of HEGBC significantly inhibited cell proliferation of GBC-SD and EH-GB2 cells (Fig. 3c, d). EdU incorporation assays also indicated that depletion of HEGBC significantly inhibited cell proliferation of GBC-SD and EH-GB2 cells (Fig. 3e). TUNEL assays indicated that depletion of HEGBC markedly promoted cell apoptosis of GBC-SD and EH-GB2 cells (Fig. 3f). Transwell assays indicated that depletion of HEGBC significantly inhibited cell migration of GBC-SD and EH-GB2 cells (Fig. 3g). Collectively, these results showed that depletion of HEGBC inhibited the proliferation and migration of GBC cells in vitro.

Ectopic expression of HEGBC promoted tumorigenesis and metastasis of GBC in vivo
Next, we investigated the roles of HEGBC on GBC tumorigenesis and metastasis in vivo. HEGBC stably overexpressed and control NOZ cells were subcutaneously injected into nude mice. The results indicated that ectopic expression of HEGBC significantly increased tumor growth rate (Fig. 4a). Also, HEGBC stably overexpressed NOZ cells formed larger subcutaneous tumors than that formed by control NOZ cells (Fig. 4b, c). Proliferation marker Ki67 immunohistochemical staining of the subcutaneous tumors indicated that the tumors formed by HEGBC overexpressed NOZ cells had higher proportion of Ki67-positive cells than that formed by control NOZ cells (Fig. 4d). TUNEL staining of the subcutaneous tumors indicated that the tumors formed by HEGBC overexpressed NOZ cells had lower proportion of apoptotic cells than that formed by control NOZ cells (Fig. 4e). To investigate the roles of HEGBC on GBC metastasis, HEGBC stably overexpressed and control NOZ cells were injected through the spleen to establish liver metastasis model in nude mice. The results indicated that ectopic expression of HEGBC increased liver metastatic foci formed by NOZ cells (Fig. 4f ). Collectively, these results showed that ectopic expression of HEGBC promoted GBC tumorigenesis and metastasis in vivo.

Depletion of HEGBC inhibited tumorigenesis and metastasis of GBC in vivo
To completely elucidate the biological roles of HEGBC on GBC tumorigenesis and metastasis in vivo, HEGBC stably depleted EH-GB2 cells were subcutaneously injected into nude mice. The results indicated that depletion of HEGBC significantly decreased tumor growth rate (Fig. 5a). Also, HEGBC stably overexpressed EH-GB2 cells formed larger subcutaneous tumors than that formed by control EH-GB2 cells (Fig. 5b, c). Ki67 immunohistochemical staining of the subcutaneous tumors indicated that the tumors formed by HEGBC depleted EH-GB2 cells had lower proportion of Ki67-positive cells than that formed by control EH-GB2 cells (Fig. 5d). TUNEL staining of the subcutaneous tumors indicated that the tumors formed by HEGBC depleted EH-GB2 cells had higher proportion of apoptotic cells than that formed by control EH-GB2 cells (Fig. 5e). HEGBC stably depleted and control EH-GB2 cells were also injected through the spleen to establish liver metastasis model in nude mice. The results indicated that depletion of HEGBC significantly decreased liver metastatic foci formed by EH-GB2 cells (Fig. 5f ). Collectively, these results showed that depletion of HEGBC inhibited GBC tumorigenesis and metastasis in vivo.

HEGBC bound to IL-11 promoter and activated IL-11/ STAT3 signaling pathway
To investigate the mechanisms mediating the roles of HEGBC in GBC, we first detected the subcellular distribution of HEGBC in GBC cells. Cytoplasmic and nuclear RNA isolation indicted that HEGBC was mainly localized in the nucleus (Fig. 6a). Several nuclear lncRNAs have been reported to directly bind the promoter of target genes and regulate the transcription of target genes [31]. The activation of STAT3 signaling has been reported to be involved in GBC initiation and progression [32,33]. Therefore, we investigated whether HEGBC regulated the expression of STAT3 signaling inducer via binding their promoters. Intriguingly, using Basic Local Alignment Search Tool (BLAST) from NCBI, we found a HEGBC binding site on the promoter of IL11 which is well-known to activate STAT3 signaling (Fig. 6b). Moreover, IL-11 was 4.88-fold higher in GBC tissues than corresponding non-tumor tissues from GSE76633 (Additional file 1: Table S1). To investigate whether HEGBC bound to the promoter of IL11, ChIRP Pearson chi-square tests were used to analyze the correlation between the expression levels of HEGBC and clinical features assays were carried out using antisense probe sets against HEGBC or LacZ (negative control). The results showed that the promoter of IL11, but not ACTB was markedly enriched in the DNA pulled down by HEGBC antisense probe compared with LacZ antisense probe and no probe (Fig. 6c), suggesting that HEGBC bound to the promoter of IL11. To investigate the effects of HEGBC on IL-11 expression, qRT-PCR was carried out on HEGBC stably overexpressed and control NOZ cells.
The results showed that ectopic expression of HEGBC up-regulated IL-11 expression (Fig. 6d). Conversely, depletion of HEGBC significantly decreased IL-11 expression in EH-GB2 cells (Fig. 6e). In addition, IL-11 concentrations in the cell supernatants of HEGBC stably The expression of HEGBC in HEGBC stably overexpressed and control NOZ cells was detected using qRT-PCR. c Cell proliferation of HEGBC stably overexpressed and control SGC-996 cells was detected using Glo cell viability assay. d Cell proliferation of HEGBC stably overexpressed and control NOZ cells was detected using Glo cell viability assay. e Cell proliferation of HEGBC stably overexpressed and control SGC-996 and NOZ cells was detected using EdU incorporation assay. The red color represents EdU-positive cells. Scale bars, 100 μm. f HEGBC stably overexpressed and control SGC-996 and NOZ cells were treated with 5 ng/mL doxorubicin for 24 h, and then cell apoptosis was detected using TUNEL assay. g Cell migration of HEGBC stably overexpressed and control SGC-996 and NOZ cells was detected using transwell assay. Scale bars, 100 μm. Results are shown as mean ± s.d. of 3 independent experiments. *P < 0.05, **P < 0.01, ***P < 0.001 by Student's t test overexpressed NOZ cells and HEGBC stably depleted EH-GB2 cells were measured by ELISA assays. The results showed that ectopic expression of HEGBC up-regulated IL-11 concentration, and while depletion of HEGBC decreased IL-11 concentration (Fig. 6f, g). Next, we investigated the effects of HEGBC on STAT3 signaling activation. Western blot assays indicated that ectopic expression of HEGBC up-regulated the phosphorylation level of STAT3, and while depletion of HEGBC decreased the phosphorylation level of STAT3 (Fig. 6h, i). Moreover, ectopic expression of HEGBC up-regulated the expression of the target genes of STAT3, and while depletion of HEGBC decreased the expression of the target genes of STAT3 (Fig. 6j, k). Collectively, these results showed that HEGBC bound to the promoter of IL11, up-regulated the expression of IL-11, induced the secretion of IL-11, and activated STAT3 signaling. The expression of HEGBC in HEGBC stably depleted and control EH-GB2 cells was detected using qRT-PCR. c Cell proliferation of HEGBC stably depleted and control GBC-SD cells was detected using Glo cell viability assay. d Cell proliferation of HEGBC stably depleted and control EH-GB2 cells was detected using Glo cell viability assay. e Cell proliferation of HEGBC stably depleted and control GBC-SD and EH-GB2 cells was detected using EdU incorporation assay. The red color represents EdU-positive cells. Scale bars, 100 μm. f HEGBC stably depleted and control GBC-SD and EH-GB2 cells were treated with 5 ng/mL doxorubicin for 24 h, and then cell apoptosis was detected using TUNEL assay. g Cell migration of HEGBC stably depleted and control GBC-SD and EH-GB2 cells was detected using transwell assay. Scale bars, 100 μm. Results are shown as mean ± s.d. of 3 independent experiments. *P < 0.05, **P < 0.01, ***P < 0.001 by Student's t test

P-STAT3 bound to HEGBC promoter and activated HEGBC expression
Once activated by cytokines and growth factors, such as IL-11, the phosphorylated STAT3 (p-STAT3) translocates to the cell nucleus, binds to the promoters of target genes, and modulates the expression of target genes. To investigate whether HEGBC is the target gene of p-STAT3, we searched the binding sites of p-STAT3 on the promoter of HEGBC using JASPAR [34] and identified four potential binding sites (Fig. 7a). ChIP assays were performed on NOZ cells using p-STAT3 specific antibody or control IgG. The results showed that p-STAT3 specifically bound to the four predicted sites on the promoter of HEGBC, but not a distal non-binding site (Fig. 7b). Several lncRNAs have been reported to physically bind STAT3, such as TSLNC8 [35]. To investigate whether HEGBC direct binds STAT3, RIP assays were performed using p-STAT3 or STAT3 specific antibody. The results showed that neither p-STAT3 nor STAT3 physically bound HEGBC (Fig. 7c). To investigate whether p-STAT3 regulates promoter activity of HEGBC, the promoter of HEGBC was subcloned into the luciferase reporter pGL3. After transfection of the constructed luciferase reporter, NOZ cells were treated with 5 μM p-STAT3 inhibitor SC144 for 72 h. Dual luciferase reporter assays showed that inhibiting p-STAT3 decreased the promoter activity of HEGBC, but not p-STAT3 binding sites mutated HEGBC promoter (Fig. 7d). Conversely, activation of p-STAT3 using IL-11 increased the promoter activity of HEGBC, but not p-STAT3 binding sites mutated HEGBC promoter (Fig. 7e). Treatment of NOZ cells with SC144 Fig. 4 Ectopic expression of HEGBC promoted GBC growth and metastasis in vivo. a HEGBC stably overexpressed and control NOZ cells were subcutaneously injected into nude mice. Tumor volumes were detected every 7 days. b The subcutaneous tumor weights were detected at the 28th day after injection. c Representative images of the subcutaneous tumors formed by HEGBC stably overexpressed and control NOZ cells. Scale bars, 1 cm. d Immunohistochemical staining of Ki67 in subcutaneous tumors formed by HEGBC stably overexpressed and control NOZ cells. Scale bars, 100 μm. e TUNEL staining of subcutaneous tumors formed by HEGBC stably overexpressed and control NOZ cells. Scale bars, 100 μm. f HEGBC stably overexpressed and control NOZ cells were injected through the spleen to establish liver metastasis. The number of liver metastatic foci was detected at the 42th day after intra-splenic injection using HE staining. Scale bars, 400 μm. Results are shown as mean ± s.d. of 5 mice in each group. *P < 0.05, **P < 0.01 by Mann-Whitney test also up-regulated the expression of HEGBC (Fig. 7f ), and while treatment of NOZ cells with IL-11 down-regulated the expression of HEGBC (Fig. 7g). The expression of IL-11 in the 102 GBC tissues in Fig. 1a was measured by qRT-PCR. Correlation analysis showed that the expression of HEGBC was significantly positively correlated with the expression of IL-11 in GBC tissues (Fig. 7h), supporting the positive modulation between HEGBC and IL-11. Collectively, these results showed that IL-11/STAT3 signaling activated the expression of HEGBC, and HEGBC and IL-11/STAT3 formed double positive regulatory loop.

Depletion of IL-11 attenuated the oncogenic roles of HEGBC in GBC cells
To investigate whether the activation of IL-11/STAT3 signaling mediates the oncogenic roles of HEGBC in GBC, we knocked-down IL-11 expression in HEGBC stably overexpressed NOZ cells via transfecting IL-11 specific shRNA (Fig. 8a). Glo cell viability assays indicated that depletion of IL-11 attenuated the proliferation promoting roles of HEGBC on NOZ cells (Fig. 8b). EdU incorporation assays also indicated that depletion of IL-11 attenuated the proliferation promoting roles of HEGBC on NOZ cells (Fig. 8c). TUNEL assays indicated that depletion of IL-11 attenuated the apoptosis inhibitory roles of HEGBC on NOZ cells (Fig. 8d). Transwell assays indicated that depletion of IL-11 attenuated the migration promoting roles of HEGBC on NOZ cells (Fig. 8e). Collectively, these results showed that depletion of IL-11 attenuated the oncogenic roles of HEGBC.

Discussion
Human transcriptome sequencing identified 58,648 lncRNAs, of which 79% were unannotated [36]. Among The subcutaneous tumor weights were detected at the 28th day after injection. c Representative images of the subcutaneous tumors formed by HEGBC stably depleted and control EH-GB2 cells. Scale bars, 1 cm. d Immunohistochemical staining of Ki67 in subcutaneous tumors formed by HEGBC stably depleted and control EH-GB2 cells. Scale bars, 100 μm. e TUNEL staining of subcutaneous tumors formed by HEGBC stably depleted and control EH-GB2 cells. Scale bars, 100 μm. f HEGBC stably depleted and control EH-GB2 cells were injected through the spleen to establish liver metastasis. The number of liver metastatic foci was detected at the 42th day after intra-splenic injection using HE staining. Scale bars, 400 μm. Results are shown as mean ± s.d. of 6 mice in each group. *P < 0.05, **P < 0.01 by Mann-Whitney test Fig. 6 HEGBC bound to the promoter of IL11 and activated IL-11/STAT3 signaling pathway. a The distribution of HEGBC in the cytoplasmic and nuclear fractions of NOZ cells was detected using cytoplasmic and nuclear RNA isolation followed by qRT-PCR. U6 and β-actin were used as nuclear and cytoplasmic controls, respectively. b Schematic outline of the predicted binding site for HEGBC on the promoter of IL11. c ChIRP assay in NOZ cells was performed with antisense probe sets against HEGBC or LacZ (negative control). The bound DNA was detected using qRT-PCR with specific primers against the promoters of IL11 or ACTB. d The expression of IL-11 in HEGBC stably overexpressed and control NOZ cells was detected using qRT-PCR. e The expression of IL-11 in HEGBC stably depleted and control EH-GB2 cells was detected using qRT-PCR. f The concentration of IL11 in the culture medium from HEGBC stably overexpressed and control NOZ cells was detected using ELISA. g The concentration of IL11 in the culture medium from HEGBC stably depleted and control EH-GB2 cells was detected using ELISA. h STAT3 phosphorylation level in HEGBC stably overexpressed and control NOZ cells was detected using western blot. i STAT3 phosphorylation level in HEGBC stably depleted and control EH-GB2 cells was detected using western blot. j The expression of the target genes of STAT3 in HEGBC stably overexpressed and control NOZ cells was detected using qRT-PCR. k The expression of the target genes of STAT3 in HEGBC stably depleted and control EH-GB2 cells was detected using qRT-PCR. Results are shown as mean ± s.d. of 3 independent experiments. **P < 0.01, ***P < 0.001 by Student's t test Fig. 7 p-STAT3 bound to the promoter of HEGBC and activated HEGBC expression. a Schematic outline of the predicted binding site for p-STAT3 on the promoter of HEGBC. b ChIP assay in NOZ cells was performed using p-STAT3 specific antibody or negative control IgG. The bound DNA was detected using qRT-PCR with specific primers against the predicted binding sites on the promoter of HEGBC or a distal non-binding site (negative control, NC). c RIP assay in NOZ cells was performed using p-STAT3 specific antibody, STAT3 specific antibody, or negative control IgG. The bound RNA was detected using qRT-PCR with specific primers against HEGBC or TSLNC8 (positive control). d After co-transfection of renilla luciferase reporter pRL-TK with firefly luciferase reporters containing HEGBC promoter, p-STAT3 binding sites mutated HEGBC promoter, or nothing into NOZ cells, the NOZ cells were treated with 5 μM p-STAT3 inhibitor SC144 for 72 h. Luciferase activities in these NOZ cells were detected. Data are presented as the relative ratio of firefly luciferase activity to renilla luciferase activity. e After co-transfection of renilla luciferase reporter pRL-TK with firefly luciferase reporters containing HEGBC promoter, p-STAT3 binding sites mutated HEGBC promoter, or nothing into NOZ cells, the NOZ cells were treated with 20 ng/mL IL-11 for 72 h. Luciferase activities in these NOZ cells were detected. Data are presented as the relative ratio of firefly luciferase activity to renilla luciferase activity. f The expression of IL-11 in NOZ cells treated with 5 μM SC144 for 72 h was detected using qRT-PCR. g The expression of IL-11 in NOZ cells treated with 20 ng/mL IL-11 for 72 h was detected using qRT-PCR. For b-g, results are shown as mean ± s.d. of 3 independent experiments. ns, not significant, **P < 0.01, ***P < 0.001 by Student's t test. h The correlation between HEGBC and IL-11 expression level in 102 GBC tissues. r = 0.7996, P < 0.0001 by Pearson's correlation analysis these lncRNAs, only a few were investigated in cancers. As to GBC, only several lncRNAs were investigated, including SPRY4-IT1 [29], CCAT1 [25], MALAT1 [37], H19 [28], GCASPC [27], UCA1 [38], HOXA-AS2 [26], PAGBC [24], AFAP1-AS1 [39], LET [40], LOC344887 [41], HOTAIR [42], and ROR [43]. Due to the huge number of lncRNAs Fig. 8 Depletion of IL-11 attenuated the oncogenic roles of HEGBC in GBC cells. a The expression of IL-11 in HEGBC stably overexpressed and control NOZ cells transfected with IL-11 specific shRNA was detected using qRT-PCR. b Cell proliferation of HEGBC stably overexpressed and control NOZ cells transfected with IL-11 specific shRNA was detected using Glo cell viability assay. c Cell proliferation of HEGBC stably overexpressed and control NOZ cells transfected with IL-11 specific shRNA was detected using EdU incorporation assay. The red color represents EdU-positive cells. Scale bars, 100 μm. d HEGBC stably overexpressed and control NOZ cells transfected with IL-11 specific shRNA were treated with 5 ng/ml doxorubicin for 24 h, and then cell apoptosis was detected using TUNEL assay. e Cell migration of HEGBC stably overexpressed and control NOZ cells transfected with IL-11 specific shRNA was detected using transwell assay. Scale bars, 100 μm. Results are shown as mean ± s.d. of 3 independent experiments. **P < 0.01, ***P < 0.001 by Student's t test in human cells, we could not rule out other lncRNAs which may also play critical roles in GBC. Therefore, we re-analyzed the public available microarray data about the differentially expressed lncRNAs in GBC. Intriguingly, we identified a novel lncRNA HEGBC, which is significantly up-regulated in GBC tissues and cell lines, compared with adjacent non-tumor tissues and non-tumorigenic biliary epithelial cell line, respectively. Increased expression of HEGBC is positively correlated with GBC extension, lymph node metastasis, and TNM stages. Furthermore, increased expression of HEGBC predicts poor outcome of GBC patients. Thus, this study identified a novel GBC-associated lncRNA HEGBC which indicates poor prognosis of GBC patients. The expression pattern and clinical significances of HEGBC in other cancers need further investigation. Subsequently, in vitro functional assays showed that enhanced expression of HEGBC increased cell viabilities, promoted cell proliferation and migration, and inhibited cell apoptosis of GBC cells. Conversely, depletion of HEGBC decreased cell viabilities, inhibited cell proliferation and migration, and induced cell apoptosis of GBC cells. In vivo xenograft assays showed that enhanced expression of HEGBC promoted GBC tumor growth and liver metastasis, and while depletion of HEGBC inhibited GBC tumor growth and liver metastasis. Therefore, our data suggested that HEGBC functions as an oncogene in GBC, and implied that HEGBC may be a potential therapeutic target for GBC. Whether the oncogenic role of HEGBC is GBC specific or tumor general needs further study.
The excessive activation of STAT3 signaling pathway has been found in many cancers, including GBC [44,45]. In response to stimulation of cytokines or growth factors, STAT3 is phosphorylated by receptor associated kinases, translocates to cell nucleus, and activates the transcription of many target genes that regulating cell viability, apoptosis, and so on. The activated STAT3 has oncogenic roles in many cancers, including GBC [46,47]. In this study, we found that HEGBC activated STAT3 signaling pathway. Further mechanistic investigations revealed that the nucleus-localized HEGBC directly bound the promoter of IL11, which stimulates the activation of STAT3. Via binding the promoter of IL11, HEGBC activated the transcription of IL11, and further upregulated the expression and secretion of IL-11. The auto-secreted IL-11 further stimulated the activation of STAT3. Intriguingly, except the activation of IL-11/STAT3 signaling by HEGBC, we also found that STAT3 directly bound the promoter of HEGBC and activated the transcription of HEGBC. Thus, HEGBC/ IL-11/STAT3 form a positive regulatory loop in GBC. The double positive regulatory roles amplified the aberrant expressions and roles of the participators. Loss-of-function assays showed that depletion of IL-11 attenuated the oncogenic roles of HEGBC in GBC. In addition, the expression of HEGBC is positively associated with that of IL-11 in GBC tissues. These data supported the positive interaction between HEGBC and IL-11 and suggested that HEGBC exerts its oncogenic roles at least partially via activating IL-11/STAT3 signaling. Although the oncogenic roles of STAT3 are well known, the roles of IL-11 in cancer are relatively unknown. Previous studies have revealed the critical roles of IL-11 in tumorigenesis of liver, colon, and gastric cancers [48,49]. In this study, we further verified that IL-11 also has critical roles in GBC. Therefore, our data suggest that HEGBC/IL-11/STAT3 regulatory loop may be potential therapeutic targets for GBC.

Conclusions
In summary, our findings identified a novel lncRNA HEGBC, which is upregulated in GBC and associated with poor prognosis of GBC patients. HEGBC promotes GBC cell proliferation and migration, inhibited GBC cell apoptosis, and promoted tumor growth and metastasis of GBC via forming a positive regulatory loop with IL-11/STAT3. Our data suggested that HEGBC may be a potential prognostic biomarker and therapeutic target for GBC.

Additional files
Additional file 1: Table S1. Top 60 differentially expressed genes in GBC. (DOCX 23 kb) Additional file 2: Figure S1. The full-length sequence of HEGBC. Representative image of PCR products from 5'-RACE and 3'-RACE assays are shown. Figure S2. RIP assay in NOZ cells was performed using RPLP0 specific antibody or negative control IgG. The bound RNA was detected using qRT-PCR with specific primers against HEGBC. Results are shown as mean ± s.d. of 3 independent experiments. ns, not significant by Student's t test. (DOCX 15 kb) Abbreviations ChIP: chromatin immunoprecipitation; ChIRP: chromatin isolation by RNA purification; EdU: ethynyl deoxyuridine; ELISA: enzyme linked immunosorbent assay; GBC: gallbladder cancer; GEO: gene expression omnibus; HEGBC: high expressed in gallbladder cancer; lncRNA: long noncoding RNA; NCBI: national center for biotechnology information; qRT-PCR: quantitative real-time polymerase chain reaction; RACE: rapid amplification of cDNA ends; RIP: RNA immunoprecipitation; TUNEL: terminal deoxynucleotidyl transferase-mediated dUTP nick end labelling

Availability of data and materials
The data used and analyzed during this study are available from the corresponding author on request.
Authors' contributions FF, KX, and LY designed the experiments; LY, QG, and XW performed the experiments; FF, KX, LY, QG, and XW analyzed the data. FF and LY wrote the manuscript. All authors read and approved the final manuscript.
Ethics approval and consent to participate Use of patient tissue samples and nude mice were approved by the Review Board of Eastern Hepatobiliary Surgery Hospital (Shanghai, China).

Consent for publication
The authors agree for publication.