- Open Access
Long non-coding RNA HOTAIR, a c-Myc activated driver of malignancy, negatively regulates miRNA-130a in gallbladder cancer
Molecular Cancervolume 13, Article number: 156 (2014)
Protein coding genes account for only about 2% of the human genome, whereas the vast majority of transcripts are non-coding RNAs including long non-coding RNAs. A growing volume of literature has proposed that lncRNAs are important players in cancer. HOTAIR was previously shown to be an oncogene and negative prognostic factor in a variety of cancers. However, the factors that contribute to its upregulation and the interaction between HOTAIR and miRNAs are largely unknown.
A computational screen of HOTAIR promoter was conducted to search for transcription-factor-binding sites. HOTAIR promoter activities were examined by luciferase reporter assay. The function of the c-Myc binding site in the HOTAIR promoter region was tested by a promoter assay with nucleotide substitutions in the putative E-box. The association of c-Myc with the HOTAIR promoter in vivo was confirmed by chromatin immunoprecipitation assay and Electrophoretic mobility shift assay. A search for miRNAs with complementary base paring with HOTAIR was performed utilizing online software program. Gain and loss of function approaches were employed to investigate the expression changes of HOTAIR or miRNA-130a. The expression levels of HOTAIR, c-Myc and miRNA-130a were examined in 65 matched pairs of gallbladder cancer tissues. The effects of HOTAIR and miRNA-130a on gallbladder cancer cell invasion and proliferation was tested using in vitro cell invasion and flow cytometric assays.
We demonstrate that HOTAIR is a direct target of c-Myc through interaction with putative c-Myc target response element (RE) in the upstream region of HOTAIR in gallbladder cancer cells. A positive correlation between c-Myc and HOTAIR mRNA levels was observed in gallbladder cancer tissues. We predicted that HOTAIR harbors a miRNA-130a binding site. Our data showed that this binding site is vital for the regulation of miRNA-130a by HOTAIR. Moreover, a negative correlation between HOTAIR and miRNA-130a was observed in gallbladder cancer tissues. Finally, we demonstrate that the oncogenic activity of HOTAIR is in part through its negative regulation of miRNA-130a.
Together, these results suggest that HOTAIR is a c-Myc-activated driver of malignancy, which acts in part through repression of miRNA-130a.
Gallbladder cancer (GBC) is the common biliary tract cancer and the fifth most common gastrointestinal malignancy . The outcome of patients with more advanced disease is dismal with 5-year survival rates ranging from 20% to 40% . Chemo-resistance and progression are the most remarkable characteristics of GBC . Although previous studies identified accumulated genomic damage promote the progression of GBC [4–6], the pathophysiological mechanism contributing to GBC are still largely unknown. Therefore, it is of paramount importance to understand the roles of novel molecules involved in this process.
It is well known that protein-coding genes account for only about 2% of the human genome, whereas the vast majority of transcripts are non-coding RNAs , and both microRNAs (21–24 nt) (miRNAs) and long ncRNAs (>200 nt) (lncRNAs) are now emerging as mammalian transcription key regulators [8, 9]. A number of studies over the past five years have identified a large number of miRNAs differentially expressed in GBC that are correlated with malignancy [10–12]. A growing volume of literature has proposed that lncRNAs are important players in cancer [13–15]. Although it is well known that miRNAs can target a number of protein-coding genes, little is known whether miRNAs/lncRNAs can also target lncRNAs/miRNAs. Recently, a competitive RNA (ceRNA) hypothesis has been proposed and several studies have suggested the interaction between lncRNA and miRNA in cancer [16–19], adding another piece of puzzle to the miRNA and lncRNA regulatory networks.
HOTAIR, is a 2158-bp lnRNA that is located at the antisense strand of the HOXC gene locus in chromo-some 12, flanked by HOXC11 and HOXC12. HOTAIR is negative prognostic factor for a variety of carcinomas, and HOTAIR expression levels are correlated with tumor metastases while loss of HOTAIR has been linked with decrease in cancer invasiveness [15, 20, 21]. The activity of HOTAIR is due, in part, to its induction of genome-wide targeting of the polycomb repressive complex 2 (PRC2), leading to an altered methylation of histone H3 lysine 27 (H3K27) and genes expression. In this study, we also find that HOTAIR is upregulated in GBC compared to adjacent normal tissues. However, although HOTAIR has been shown to play a vital role in cancer, the factors that contribute to its upregulation and the interaction between HOTAIR and miRNAs are largely unknown.
In the present study, we demonstrate that HOTAIR is a direct target of c-Myc through interaction with putative c-Myc target response element (RE) in the upstream region of HOTAIR. Our data also reveal that HOTAIR negatively regulates miRNA-130a and the oncogenic activity of HOTAIR is, at least in part, through the negative regulation of miRNA-130a, which may function as a part of the “competitive endogenous RNAs (ceRNA)” network .
Material and method
Four human GBC cell lines (GBC-SD, SGC-996,NOZ and EH-GB2) were used in this study. GBC-SD and SGC-996 were purchased from Cell Bank of the Chinese Academy of Science (Shanghai, China). NOZ was purchased from the Health Science Research Resources Bank (Osaka, Japan). EH-GB2 was a generous gift from Eastern Hepatobiliary Surgical Hospital and Institute, The Second Military University, Shanghai . The cell lines were cultured in Dulbecco’s modified Eagle’s medium (DMEM, Gibco BRL), containg 10% fetal calf serum (FBS, HyClone) as well as 100 U/ml penicillin and 100 μg/ml streptomycin (Invitrogen). Cells were maintained in a humidified incubator at 37°C in the presence of 5% CO2. All cell lines have been passaged for fewer than 6 months.
Patients and samples
Sixty-five GBC tissue samples and neighboring noncancerous gallbladder tissue samples (collected postoperatively from August 2007 to September 2010) used in this study were obtained from Xinhua Hospital, Shanghai, China. Upon removal of the surgical specimen, surgical pathology faculty performed a gross analysis of the specimen and selected gallbladder tissues that appeared to be cancerous and gallbladder tissues that appeared to be normal for analysis. Each sample was snap-frozen in liquid nitrogen and stored at -80°C prior to RNA isolation and qRT-PCR analysis. All patients recruited to this study did not receive any pre-operative treatments. GBC patients were staged according to the tumor node metastasis (TNM) staging system (the 7th edition) of the American Joint Committee on Cancer (AJCC) staging system. The data do not contain any information that could identify the patients. All patients provided written informed consent. Complete clinicopathological follow-up data of the GBC patients from which the specimens were collected were available. This study was approved by the Human Ethics Committee of Xinhua Hospital at Shanghai Jiaotong University (Shanghai, China).
RNA preparation, reverse transcription and qPCR
Total RNA from tissues and cells was extracted using Trizol reagent (Invitrogen, CA). RNA was reversed transcribed into cDNAs using the Primer-Script™ one step RT-PCR kit (TaKaRa, Dalian, China). The cDNA template was amplified by real-time RT-PCR using the SYBR® Premix Dimmer Eraser kit (TaKaRa, Dalian, China). GAPDH was used as an internal control, and mRNA values were normalized to GAPDH. Real-time RT-PCR reactions were performed by the ABI7500 system (Applied Biosystems, CA). The relative expression fold change of mRNAs was calculated by the 2-ΔΔCt method. The primer sequences were listed in the Additional file 1.
Expression plasmids for HOTAIR, c-Myc or corresponding mutants of HOTAIR by mutating the has-miRNA-130a seed region binding site (seed sequence binding fragment 5’-GACTTTGCACT -3’ changed to 5’-TTGTAACGTGA-3’) were created using PCR amplification with human genomic DNA as templates. The primer sequences used are listed in Additional file 1. The PCR product was verified and subcloned into the mammalian expression vector pcDNA3.1 (Invitrogen). Plasmids were transfected into cells with lipofectamine 2000 (Invitrogen, CA).
The promoter region of HOTAIR was PCR-amplified by PrimerStar polymerase (TaKaRa) with the primers 5’-ACTGGTACCTAAGCGGAGAGAGTCCC-3’ (forward) and 5’-ACTAAGCTTGAGTCAGAGTTCCCCAC-3’ (reverse) and was subcloned into the pGL3 basic firefly luciferase reporter. The pGL3 construct containing the HOTAIR promoter with a point mutation in the E-box element was PCR-amplified by PrimerStar polymerase (TaKaRa) with the primers 5’-ACTGGTACCTAAGCGGAGAGAGTCCCACACAGG-3’ (forward) and 5’-ACTAAGCTTGAGTCAGAGTTCCCCAC-3’ (reverse) and was subcloned into the pGL3 basic firefly luciferase reporter.
HOTAIR siRNA, c-Myc siRNA and Allstars Negative Control siRNA were purchased from Qiagen, Hilden, Germany. Target sequences are listed in the Additional file 1. Hsa-miRNA-130a mimic/negative control mimic and hsa-miRNA-130a inhibitor/negative control inhibitor were purchased from Genechem, shanghai, China.
Cells were grown on six-well plates to confluency and transfected using Lipofectamine 2000 (Invitrogen) according to the manufacturer’s instructions. Forty-eight hours after transfection, cells were harvested for qRT-PCR or western blot analyses.
Luciferase reporter assay
For the promoter activity of HOTAIR, the promoter/luciferase reporter construct and pcDNA3.1-c-Myc or c-Myc siRNA were cotransfected into cultured cells by Lipofectamine-mediated gene transfer. Each sample was cotransfected with the pRL-TK plasmid, which expressed renilla luciferase to monitor transfection efficiency (Promega, Madison, WI, USA). The relative luciferase activity was normalized with renilla luciferase activity.
For the dual luciferase activity, LncRNA-HOTAIR (lncRNA-HOTAIR-wt) or its mutant devoid of specific miRNA binding sites (lncRNA-HOTAIR-mu) was cloned into 3’UTR of the Renilla luciferase gene in the vector pRL-TK (Promega, Madison, WI, USA). Each plasmid was transfected into cells, together with specific miRNAs mimics or with a negative control mimic (RiboBio, Guangzhou, China). Firefly luciferase gene in the vector pGL3-control (Promega) was used as a control for transfection efficiency. Luciferase assays were performed using the dual-luciferase reporter assay system kit (Promega) according to the manufacturer’s instructions. Luciferase expression was analyzed by Modulus single-tube multimode reader (Promega). The relative luciferase expression equals the expression of Renilla luciferase (pRL-TK) divided by the expression of firefly luciferase.
Chromatin immunoprecipitation (ChIP) assay
Chromatin immunoprecipitation (ChIP) assays were performed according to the EZ ChIP Chromatin Immunoprecipitation Kit (Millipore, Bedford, MA, USA). Briefly, crosslinked chromatin was sonicated into 200- to 1,000-bp fragments. Anti-c-Myc and anti-Max antibodies (Cell Signal Technology, USA) were used to precipitate DNA–protein complexes. Normal mouse immunoglobulin G (IgG) was used as a negative control. ChIP-derived DNA was quantified using qRT-PCR with SYBR-Green incorporation (Applied Biosystems, Foster City, CA, USA). The primers are listed in Additional file 1.
Electrophoretic mobility shift assay (EMSA)
The following double-stranded oligonucleotides were used (wild type and mutant binding sites are underlined): E-box, 5’-CGAGCGCAGTGGCGCATGGCTGTAATCCCA-3’; E-box mutant, 5’-CGAGCGCAGTGGCATGGGGCTGTAATCCCA-3’. Oligonucleotide labeling was performed using the Biotin 3’ End Labeling Kit (Pierce, USA). EMSA was performed using a light shift chemiluminescent EMSA kit (Pierce). Nuclear proteins from gallbladder cancer cells were prepared using NE-PER nuclear and cytoplasmic extraction reagents (Pierce). Nuclear extract proteins (4 μg) were incubated in with binding buffer (2.5% glycerol, 5 mM MgCl2, 0.05% NP-40, 1 μg poly(dI-dC), 10 mM Tris, 50 mM KCl, 1 mM DTT, pH 7.5). Samples were electrophoresed on a 6% polyacrylamide gel in 0.5× TBE buffer (45 mM Tris, 45 mM boric acid, 1 mM EDTA, pH 8.3). For competition assays, samples were preincubated with a 200-fold excess (4 pmol) of the unlabeled wild type competitors for 20 min. For the supershift reaction, 1 μg of each anti-c-Myc antibody was preincubated with the nuclear extracts in the absence of poly (dI · dC) for 1 hour at 4°C. Subsequently, poly (dI · dC) was added and incubated for 5 min, followed by the incubation of the biotin labeled probes (20 fmol) for 20 min. Samples were separated by electrophoresis on a 6% non-denaturing acrylamide gel in 0.5 × TBE, transferred to positively charged nylon membranes, and visualized by streptavidin-horseradish peroxidase followed by chemiluminescent detection.
Northern blot analysis
Total RNA (15 μg) from samples were separated on 15% denaturing polyacrylamide gels, transferred onto GeneScreen Plus membranes (PerkinElmer), and hybridized using UltraHyb-Oligo buffer (Ambion). Oligonucleotides complementary to mature miR-130a (5’-AGCAAAAATGTGCTAGTGCCAAA-3’) were end-labeled with T4 Kinase (Invitrogen) and used as probes. Following hybridization at 42°C overnight, the membranes washed twice in 0.1 × SSPE and 0.1% SDS at 42°C for 15 min each. Membranes were then exposed to a storage phosphor screen (GE Healthcare Bio-Sciences) for 8 h and imaged using a Typhoon 9410 Variable Mode Imager (GE Healthcare Bio-Sciences). Northern blots hybridized with a 5S ribosomal RNA (rRNA) cDNA were used as controls.
RNA pull-down assay
To determine whether HOTAIR is associated with the RNA-induced silencing complex (RISC) complex, we performed RNA pull-down assay using synthesized biotin-labeled HOTAIR as a probe and then detected Ago 2 from the pellet by western or detected miRNA-130a by quantitative RT-PCR (qRT-PCR).
RNA pull-down were performed as described previously . Briefly, the DNA fragment covering has-miRNA-130a seed region binding site of HOTAIR was PCR-amplified using a T7 containing primer and then cloned into pCR8 (Invitrogen). In addition, lncRNA loc285194  was also cloned and used in RNA Pull-Down Assay as a positive control. The resultant plasmid DNA was linearized with restriction enzyme Not I. Biotin-labeled RNAs were in vitro transcribed with the Biotin RNA Labeling Mix (Roche Diagnostics, Indianapolis, IN) and T7 RNA polymerase (Roche), treated with RNase-free DNase I (Roche), and purified with the RNeasy Mini Kit (Qiagen, Inc.,Valencia, CA). Cell nuclear extract (2 ug) was mixed with biotinylated RNA (100 pmol). Washed Streptavidin agarose beads (100 ml) were added to each binding reaction and further incubated at room temperature for 1 h. Beads were washed briefly three times and boiled in SDS buffer, and the retrieved protein was detected by standard western blot technique.
The Ago2 antibodies used for RIP are purchased from Abcam (Abcam, Cambridge, MA). The coprecipitated RNAs were detected by reverse transcription PCR. Total RNAs and controls were also assayed to demonstrate that the detected signals were from RNAs specifically binding to Ago2.
Western blot analysis to assess c-Myc, Ago2 and GADPH expression was carried out as described previously . GADPH primary antibodies were purchased from Sigma (MO, USA).
Cell invasion assay
For the invasion assays, 48 h after transfection, 5 × 104 cells in serum-free media were placed into the upper chamber of an insert (8.0 μm, Millipore, MA) coated with Matrigel (Sigma, USA). The chambers were then incubated for 24 h in culture medium with 10% FBS in the bottom chambers before examination. The cells on the upper surface were scraped and washed away, whereas the invaded cells on the lower surface were fixed and stained with 0.05% crystal violet for 2 h. Finally, invaded cells were counted under a microscope and the relative number was calculated. Experiments were independently repeated in triplicate.
Flow cytometric analysis
Cells were seeded at a density of 1 × 106 cells/well in six-well plates. After 24 h, cells were washed with PBS and fixed in ice-cold 70% ethanol for 1 h and then treated with 100 uL of 50 mg/L propidium iodide for 30 min at 4°C in the dark. The cell-cycle profiles were assayed using the Elite ESP flow cytometer at 488 nm, and data were analyzed with the CELL Quest software (BD Biosciences,San Jose, CA, USA).
All statistical analyses were performed using SPSS 17.0 (SPSS, Chicago, USA). The gene expression level of HOTAIR in tumors was compared with adjacent normal tissues utilizing paired samples t-test. The expression differences between high/low grades, high/low stages, cell lines, the expression changes after transfection, luciferase activity, cell cycle and cell migration assays were analyzed using independent samples t-test. All data were presented as mean ± standard error. A two-sided p value of less than 0.05 was considered to be statistically significant.
HOTAIR is upregulated in gallbladder cancer tissues
The expression level of HOTAIR was examined using real-time PCR in 65 pairs of gallbladder cancer tissues and matched adjacent normal tissues. Detailed clinical features are presented in Table 1. As illustrated in Figure 1A, the HOTAIR transcripts were expressed at higher levels in the tumor tissues compared with adjacent normal tissues (p < 0.0001, Figure 1A), indicating that HOTAIR was frequently up-regulated in GBC. Next, we examined the expression level of HOTAIR with clinical characteristics in 65 gallbladder cancer tissues (Figure 1B, C). HOTAIR was more highly expressed in tumors extending beyond the gallbladder (T3 + T4) compared with tumors only detected in the gallbladder (T1 + T2) (Figure 1B) and more highly expressed in tumors spread to regional lymph nodes (N1) compared with tumors localized only in the gallbladder (N0) (Figure 1C). We then examined the expression level of HOTAIR in a series of gallbladder cancer cell lines (Figure 1D). We selected gallbladder cancer cell line GBC-SD as our experimental model as GBC-SD harbors a moderate expression level of HOTAIR, which makes it easier for manipulation.
c-Myc directly binds to the promoter regions of HOTAIR and upregulates its expression
To determine how transcription of HOTAIR was regulated, we then performed a search for possible transcription factor-binding sites using online software programs.
MatInspector (http://www.genomatix.de/online_help/help_matinspector/matinspector_help.html) and TFSEARCH (http://www.cbrc.jp/research/db/TFSEARCH.html) in the promoter region of HOTAIR and found that one putative E-box element at ~1053 upstream of HOTAIR which could be recognized by c-Myc (Additional file 2: Figure S1). c-Myc has been reported to be an oncoprotein and is deregulated in gallbladder cancer [5, 24].
Firstly, we would like to evaluate whether c-Myc regulates the expression of HOTAIR. We manipulated the expression of c-Myc by transfecting pcDNA3.1-c-Myc (Figure 2A, B) or siRNA into gallbladder cancer cells GBC-SD (Figure 2D, E). At 48 h after treatment, we measured the expression level of HOTAIR. Our results indicated that ectopic expression of c-Myc enhanced the expression of HOTAIR (Figure 2C) while inhibition of c-Myc decreased the expression of HOTAIR (Figure 2F).
To further explore the mechanism of c-Myc-induced HOTAIR upregulation in GBC cells, we cloned the promoter of HOTAIR into the pGL3 basic firefly luciferase reporter (Promega) and cotransfected the construct with pcDNA3.1-c-Myc or c-Myc specific siRNA into GBC-SD cells. Overexpression of c-Myc dramatically increased the luciferase activity of construct while inhibition of c-Myc decreased the luciferase activity (Figure 3A, B). We then mutated the E-box element, as depicted in Figure 3C, and cotransfected cells with pcDNA3.1-c-Myc or c-Myc siRNA. Our results indicated that mutation of the E-box element abolished the effects of c-Myc on the promoter activity of HOTAIR (Figure 3D, E). The data suggest that c-Myc may activate the HOTAIR promoter. c-Myc binds to DNA sequence elements called E-boxes to allow for target gene transcription with its obligate heterodimerization partner Max. In a ChIP assay, c-Myc and Max immunoprecipitates were highly enriched in the DNA fragments compared with negative control IgG immunoprecipitates (Figure 4A). A random region which does not contain an E-box site, showed no significant enrichment (Figure 4B). To confirm the physical interaction of c-Myc with the putative E-box element in HOTAIR promoter, we performed the electrophoretic mobility shift assay using biotin-labeled, synthetic double-stranded oligonucleotides corresponding to the E-box. As shown in Figure 4C, a binding complex was formed between c-Myc and the labeled wild-type HOTAIR oligonucleotides. However, the binding was not observed with mutated probe. Furthermore, the bands were supershifted by specific c-Myc antibodies. Together, these results suggest that c-Myc interacts with the c-Myc responsive element in the HOTAIR promoter to induce its transcription.
Reciprocal repression of miRNA-130a and HOTAIR
The function of lncRNAs in human diseases may have to do with their ability to regulate gene expression. Recently, increasing evidence suggests that non-coding RANs may participate in ceRNA regulatory network . For example, there is a negative correlation between loc285194 and miRNA-211 . We performed a search for miRNAs with complementary base paring with HOTAIR utilizing online software program starbase v2.0 (http://starbase.sysu.edu.cn/mirLncRNA.php) . We found that 20 miRNAs formed complementary base pairing with HOTAIR (Table 2). To determine whether any of them is truly regulated by HOTAIR as predicted, we profiled the expression of the 20 miRNAs in GBC-SD transfected with si-NC or si-HOTAIR. The initial profiling identified two miRNAs (miRNA-326, miRNA-130a) that were with a fold-change greater than 2 compared to the control (Table 3). We then focused on miRNA-130a (Figure 5A), which is of the greatest fold-change in response to HOTAIR knockdown.HOTAIR siRNAs significantly reduced the endogenous HOTAIR (Figure 5B); at the same time, HOTAIR siRNA increased the miRNA-130a level (Figure 5C, D). In contrast, ectopic expression of HOTAIR increased the expression level of HOTAIR (Figure 6A) while dramatically suppressed miRNA-130a (Figure 6B, C). To determine whether this suppression is through the potential interaction at the putative miRNA-130a-binding site, we generated a HOTAIR mutant (Figure 6D). This mutant HOTAIR clone revealed no significant suppression of miRNA-130a compared with wide-type of HOTAIR (Figure 6B, C). To determine whether miRNA-130a is able to negatively regulate HOTAIR, we also transfected miRNA-130a mimic into GBC-SD cells. As shown in Figure 7A, miRNA-130a mimic reduced the HOTAIR level by approximately 64%. Furthermore, miRNA-130a inhibitor increased the expression of HOTAIR (Figure 7B). To further confirm that the miR-130a target site is functional, luciferase reporter constructs were generated (Figure 7C). WT HOTAIR or HOTAIR mutant devoid of the miR-130a binding site was cloned downstream of Renilla luciferase gene and transfected into 293 T cells together with specific miR-130a mimics or the negative control mimic. The data revealed that luciferase expression was obviously reduced in cells transfected with HOTAIR and miR-130a mimics compared with that in cells transfected with HOTAIR and negative control. However, luciferase expression in cells transfected with HOTAIR mutant and the miR-130a mimics was comparable to that of control cells (Figure 7D). These data demonstrates that the binding sites are vital for the reciprocal repression of HOTAIR and miRNA-130a.
To explore the mechanism of the negative regulation of miRNA-130a by HOTAIR, we examined the effect of knockdown of HOTAIR on the expression level of mature miRNA-130a, pri-miRNA-130a and pre-miRNA-130a. As demonstrated in Figure 8A, while HOTAIR siRNA induced a significant upregulation of mature miRNA-130a, it had no effect on pri-miRNA-130a or pre-miRNA-130a, implying this negative regulation might be through a post-transcriptional mechanism. It is well known that miRNAs exert their gene silencing functions through a ribonucleoprotein complex called the RNA induced silencing complex (RISC) . Potential microRNA targets can be isolated from this complex after Ago2 co-immunprecipitation [16, 27] as Ago2 is a vital component of RISC complex necessary for siRNA or miRNA-mediated gene silencing. It does not exclude the possibility that miRNA-130a might be in separate RISC complexes even if we detected both miRNA-130a and HOTAIR in the Ago2 pellet with the Ago2 co-immunprecipitation assay. To determine whether miRNA-130a and HOTAIR are in the same RISC complex, we performed RNA pull-down experiments using HOTAIR probe and then examined Ago2 and miRNA-130a simultaneously as described previously [16, 17]. A biotin-labeled HOTAIR RNA probe was synthesized and mixed with cellular extract . We next performed in vitro RNA pull-down to validate the association between HOTAIR and Ago2. As a result, we detected Ago2 (Figure 8B). What’s more, we detected miRNA-130a in the same pellet (Figure 8C). Thus, these results indicate that both HOTAIR and miRNA-130a are probably in the same Ago2 complex.
HOTAIR expression is positively correlated with c-Myc and negatively correlated with miRNA-130a in gallbladder cancer tissues
As c-Myc upregulates the expression of HOTAIR, we investigated whether a correlation exists between c-Myc and HOTAIR expression levels in gallbladder cancer tissues. We examined the expression level of c-Myc transcript in gallbladder cancer tissues from Figure 1. The c-Myc mRNA levels in cancer tissues were significantly higher than those in adjacent normal tissues (p < 0.0001, Figure 9A). A significant positive correlation was observed between HOTAIR and c-Myc mRNAs (r = 0.7063, p < 0.0001, Figure 9B), supporting the role of c-Myc in the expression of HOTAIR. We also determined the expression levels of miRNA-130a in gallbladder cancer tissues from Figure 1. The miRNA-130a mRNA was markedly downregulated in gallbladder cancer tissues compared to adjacent normal tissues (p < 0.0001, Figure 9C). HOTAIR expression was negatively correlated with miRNA-130a in gallbladder cancer tissues (r = -0.6398, p < 0.0001, Figure 9D), providing evidence to the reciprocal negative regulation of HOTAIR and miRNA-130a.
HOTAIR’s oncogenic activity is in part through negative regulation of miRNA-130a
To investigate the biological roles of HOTAIR and miRNA-130a in gallbladder cancer, we employed gain-of-function and loss-of-function approaches. We demonstrated that knockdown of HOTAIR inhibited the invasion of gallbladder cancer cells while miRNA-130a inhibitor reversed the decrease in invasiveness (Figure 10A, B). Flow cytometric analysis indicated knockdown of HOTAIR suppressed cancer cells proliferation (S-phase fraction) in vitro while miRNA-130a inhibitor rescued the proliferation (Figure 10C, D). These results may imply that the oncogenic activity of HOTAIR is partly through negative regulation of miRNA-130a.
As a new class of non-coding RNAs, long non-coding RNAs were found to be deregulated in a variety of diseases, especially cancer [13–15]. Understanding the precise molecular mechanism by which lncRNAs function is vital for exploring new potential strategies for early diagnosis and therapy. In this study, we present evidence that HOTAIR is a direct target of c-Myc and exhibits oncogenic activity partly through negative regulation of miRNA-130a. This study provides experimental evidence to the existence of the ceRNA regulatory network  where HOTAIR and miRNA-130a negatively regulate each other.
HOTAIR, a long non-coding RNA initially identified in breast cancer, was shown to be upregulated in a variety of carcinomas [15, 28, 29]. A large number of studies have focused on its biological role and association with clinical prognosis, yet the precise factors regulating its expression remains largely unknown expect that HOTAIR is transcriptionally regulated by estradiol in breast cancer , which is quite tumor specific. In this study, we predicted a putative binding site of c-Myc in the promoter region of HOTAIR. c-Myc is a well-known transcriptional factor, with its obligate heterodimerization partner Max, binds to DNA sequence elements called E-boxes to allow for target gene transcription . Recent studies have suggested that c-Myc can regulate numerous protein-coding and non-coding genes, expecially miRNAs .
In this study, we demonstrated that c-Myc induced HOTAIR expression through direct interaction with the E-box in the HOTAIR promoter region. Ectopic expression of c-Myc increased HOTAIR expression and its promoter activity, while knockdown of c-Myc reduced HOTAIR expression and its promoter activity. Nucleotide mutant in the E-box element in the promoter region abrogated c-Myc-dependent promoter activation. The association of c-Myc with the HOTAIR promoter was confirmed by chromatin immunoprecipitation assays. What’s more, a positive correlation between c-Myc and HOTAIR mRNAs was observed in gallbladder cancer tissues, providing additional evidence to c-Myc’s regulation of HOTAIR. At present, two important lncRNAs have been shown to be c-Myc transcription targets [31, 33], involved in the c-Myc mediated cellular process. Our study would strengthen the notion that, lncRNAs are also an important part of c-Myc regulatory network.
Previous studies suggested that HOTAIR interacts with PRC2 (polycomb repressive complex 2) and histone demethylases LSD1/CoREST/REST complexes through it 5’- and 3’-end, respectively. HOTAIR acts as a bridge coordinating the targeting of PRC2 and LSD1 complexes to chromatin for coupled histone H3K27 methylation and H3 lysine-4 (H3K4) demethylation processes, which leads to silencing of target genes [15, 34]. A growing volume of recent work has established that lncRNAs can also regulate other non-coding RNAs, in particular miRNAs, and miRNAs may have an effect on the regulation of lncRNAs [18, 19]. In support of this notion, we demonstrate that HOTAIR-mediated oncogenic activity is at least partly through suppression of miRNA-130a. Knockdown of HOTAIR induced the upregulation of miRNA-130a. Ectopic expression of HOTAIR reduced the miRNA-130a level and the miRNA-130a-binding site is vital for the HOTAIR-mediated repression. On the other hand, miRNA-130a inhibitor upregulated HOTAIR level while miRNA-130a mimic repressed HOTAIR level. HOTAIR and miRNA-130a may form a reciprocal repression feedback loop. In addition, a negative correlation was observed between HOTAIR and miRNA-130a in gallbladder cancer tissues, providing supporting evidence to such a feedback loop. In addition, we explored the mechanism of such a feedback loop. We found that HOTAIR and miRNA-130a bind to the same RISC complex. As miRNAs are known to mediate post-transcriptional control of gene expression by binding to the 3’-untranslated regions of protein coding genes, we suppose that the way that miR-130a promoted the downregulation of HOTAIR is somewhat similar to the miRNA-mediated silencing of protein-coding genes. HOTAIR is well-known for its induction of genome-wide targeting of the polycomb repressive complex 2 (PRC2), leading to an altered methylation of histone H3 lysine 27 (H3K27) and genes expression. Thus, we may hypothesize that HOTAIR may lead to the downregulation of miR-130a via increasing the methylation status of the promoter of miRNA-130a as Vrba et al.  demonstrated that the downregulation of miR-130a is linked to increased promoter methylation.
In this regard, HOTAIR may function as the endogenous sponge, similar to what has been reported for lncRNA loc285194 , GAS5  and HULC . For example, loc285184 is identified as a tumor-suppressive lncRNA in colon cancer. Loc285194 can downregulate miRNA-211 through its interaction with miRNA-211. Though much of the focus in ncRNA research is directed towards understanding the regulation of protein-coding genes mediated by them, it may seem that ncRNAs could form a well-orchestrated regulatory interaction network [18, 22]. Our study further demonstrated that the reciprocial repression of HOTAIR and miRNA-130a is likely through the pathway involving RISC complex.
miRNA-130a has been found to be downregulated in a variety of carcinomas and exhibits tumor-suppressive activity [37–39] while HOTAIR was demonstrated be an oncogene and upregulated in carcinomas [28–30]. We then then studied that biological function of HOTAIR and miRNA-130a in gallbladder cancer cells. HOTAIR was shown to promote invasiveness and proliferation of cancer cells [28, 29]. Our data showed that while knockdown of HOTAIR inhibited the invasiveness and proliferation of gallbladder cancer cells, miRNA-130a inhibitor reversed the effects that knockdown of HOTAIR exerted. These results provide additional evidence to the reciprocal repression loop in a functional aspect.
Our study suggests another layer of regulation that involves the ncRNAs. A better understanding of the ncRNA interaction regulatory network would deepen our understanding of the pathophysiological mechanism of various aspects of tumorigenesis including tumor growth, invasion, metastasis and chemo-resistance in gallbladder cancer.
Ming-zhe Ma, Chun-xiao Li and Yan Zhang are co-first authors
Zhu AX, Hong TS, Hezel AF, Kooby DA: Current management of gallbladder carcinoma. Oncologist. 2010, 15: 168-181. 10.1634/theoncologist.2009-0302
Gold DG, Miller RC, Haddock MG, Gunderson LL, Quevedo F, Donohue JH, Bhatia S, Nagorney DM: Adjuvant therapy for gallbladder carcinoma: the Mayo Clinic Experience. Int J Radiat Oncol Biol Phys. 2009, 75: 150-155. 10.1016/j.ijrobp.2008.10.052
Caldow Pilgrim CH, Groeschl RT, Quebbeman EJ, Gamblin TC: Recent advances in systemic therapies and radiotherapy for gallbladder cancer. Surg Oncol. 2013, 22: 61-67. 10.1016/j.suronc.2012.12.001
Jiao Y, Pawlik TM, Anders RA, Selaru FM, Streppel MM, Lucas DJ, Niknafs N, Guthrie VB, Maitra A, Argani P, Offerhaus GJ, Roa JC, Roberts LR, Gores GJ, Popescu I, Alexandrescu ST, Dima S, Fassan M, Simbolo M, Mafficini A, Capelli P, Lawlor RT, Ruzzenente A, Guglielmi A, Tortora G, de Braud F, Scarpa A, Jarnagin W, Klimstra D, Karchin R: Exome sequencing identifies frequent inactivating mutations in BAP1, ARID1A and PBRM1 in intrahepatic cholangiocarcinomas. Nat Genet. 2013, 45: 1470-1473. 10.1038/ng.2813
Ooi A, Suzuki S, Nakazawa K, Itakura J, Imoto I, Nakamura H, Dobashi Y: Gene amplification of Myc and its coamplification with ERBB2 and EGFR in gallbladder adenocarcinoma. Anticancer Res. 2009, 29: 19-26.
Nagahashi M, Ajioka Y, Lang I, Szentirmay Z, Kasler M, Nakadaira H, Yokoyama N, Watanabe G, Nishikura K, Wakai T, Shirai Y, Hatakeyama K, Yamamoto M: Genetic changes of p53, K-ras, and microsatellite instability in gallbladder carcinoma in high-incidence areas of Japan and Hungary. World J Gastroenterol. 2008, 14: 70-75. 10.3748/wjg.14.70
Djebali S, Davis CA, Merkel A, Dobin A, Lassmann T, Mortazavi A, Tanzer A, Lagarde J, Lin W, Schlesinger F, Xue C, Marinov GK, Khatun J, Williams BA, Zaleski C, Rozowsky J, Röder M, Kokocinski F, Abdelhamid RF, Alioto T, Antoshechkin I, Baer MT, Bar NS, Batut P, Bell K, Bell I, Chakrabortty S, Chen X, Chrast J, Curado J: Landscape of transcription in human cells. Nature. 2012, 489: 101-108. 10.1038/nature11233
Kugel JF, Goodrich JA: Non-coding RNAs: key regulators of mammalian transcription. Trends Biochem Sci. 2012, 37: 144-151. 10.1016/j.tibs.2011.12.003
Fachel AA, Tahira AC, Vilella-Arias SA, Maracaja-Coutinho V, Gimba ER, Vignal GM, Campos FS, Reis EM, Verjovski-Almeida S: Expression analysis and in silico characterization of intronic long noncoding RNAs in renal cell carcinoma: emerging functional associations. Mol Cancer. 2013, 12: 140- 10.1186/1476-4598-12-140
Peng HH, Zhang YD, Gong LS, Liu WD, Zhang Y: Increased expression of microRNA-335 predicts a favorable prognosis in primary gallbladder carcinoma. Onco Targets Ther. 2013, 6: 1625-1630.
Qiu Y, Luo X, Kan T, Zhang Y, Yu W, Wei Y, Shen N, Yi B, Jiang X: TGF-β upregulates miR-182 expression to promote gallbladder cancer metastasis by targeting CADM1. Mol Biosyst. 2014, 10: 679-685. 10.1039/c3mb70479c
Kono H, Nakamura M, Ohtsuka T, Nagayoshi Y, Mori Y, Takahata S, Aishima S, Tanaka M: High expression of microRNA-155 is associated with the aggressive malignant behavior of gallbladder carcinoma. Oncol Rep. 2013, 30: 17-24.
Yang F, Huo XS, Yuan SX, Zhang L, Zhou WP, Wang F, Sun SH: Repression of the long noncoding RNA-LET by histone deacetylase 3 contributes to hypoxia-mediated metastasis. Mol Cell. 2013, 49: 1083-1096. 10.1016/j.molcel.2013.01.010
Huang JF, Guo YJ, Zhao CX, Yuan SX, Wang Y, Tang GN, Zhou WP, Sun SH: Hepatitis B virus X protein (HBx)-related long noncoding RNA (lncRNA) down-regulated expression by HBx (Dreh) inhibits hepatocellular carcinoma metastasis by targeting the intermediate filament protein vimentin. Hepatology. 2013, 57: 1882-1892. 10.1002/hep.26195
Gupta RA, Shah N, Wang KC, Kim J, Horlings HM, Wong DJ, Tsai MC, Hung T, Argani P, Rinn JL, Wang Y, Brzoska P, Kong B, Li R, West RB, van de Vijver MJ, Sukumar S, Chang HY: Long non-coding RNA HOTAIR reprograms chromatin state to promote cancer metastasis. Nature. 2010, 464: 1071-1076. 10.1038/nature08975
Liu Q, Huang J, Zhou N, Zhang Z, Zhang A, Lu Z, Wu F, Mo YY: LncRNA loc285194 is a p53-regulated tumor suppressor. Nuclei Acids Res. 2013, 41: 4976-4987. 10.1093/nar/gkt182.
Zhang Z, Zhu Z, Watabe K, Zhang X, Bai C, Xu M, Wu F, Mo YY: Negative regulation of lncRNA GAS5 by miR-21. Cell Death Differ. 2013, 20: 1558-1568. 10.1038/cdd.2013.110
Jalali S, Bhartiya D, Lalwani MK, Sivasubbu S, Scaria V: Systematic transcriptome wide analysis of lncRNA-miRNA interactions. PLoS One. 2013, 8: e53823- 10.1371/journal.pone.0053823
Juan L, Wang G, Radovich M, Schneider BP, Clare SE, Wang Y, Liu Y: Potential roles of microRNAs in regulating long intergenic noncoding RNAs. BMC Med Genomics. 2013, 6 (Suppl 1): S7-10.1186/1755-8794-6-S1-S7.
Kim K, Jutooru I, Chadalapaka G, Johnson G, Frank J, Burghardt R, Kim S, Safe S: HOTAIR is a negative prognostic factor and exhibits pro-oncogenic activity in pancreatic cancer. Oncogene. 2013, 32: 1616-1625. 10.1038/onc.2012.193
Li X, Wu Z, Mei Q, Li X, Guo M, Fu X, Han W: Long non-coding RNA HOTAIR, a driver of malignancy, predicts negative prognosis and exhibits oncogenic activity in oesophageal squamous cell carcinoma. Br J Cancer. 2013, 109: 2266-2278. 10.1038/bjc.2013.548
Karreth FA, Pandolfi PP: ceRNA cross-talk in cancer: when ce-bling rivalries go awry. Cancer Discov. 2013, 3: 1113-1121. 10.1158/2159-8290.CD-13-0202
Wang JH, Li LF, Yu Y, Li B, Jin HJ, Shen DH, Li J, Jiang XQ, Qian QJ: Establishment and characterization of a cell line, EH-GB2, derived from hepatic metastasis of gallbladder cancer. Oncol Rep. 2012, 27: 775-782.
Yukawa M, Fujimori T, Hirayama D, Idei Y, Ajiki T, Kawai K, Sugiura R, Maeda S, Nagasako K: Expression of oncogene products and growth factors in early gallbladder cancer, advanced gallbladder cancer, and chronic cholecystitis. Hum Pathol. 1993, 24: 37-40. 10.1016/0046-8177(93)90060-T
Li JH, Liu S, Zhou H, Qu LH, Yang JH: starBase v2.0: decoding miRNA-ceRNA, miRNA-ncRNA and protein-RNA interaction networks from large-scale CLIP-Seq data. Nucleic Acids Res. 2014, 42 (Database issue): D92-D97.
Gregory RI, Chendrimada TP, Cooch N, Shiekhattar R: Human RISC couples microRNA biogenesis and posttranscriptional gene silencing. Cell. 2005, 123: 631-640. 10.1016/j.cell.2005.10.022
Karginov FV, Conaco C, Xuan Z, Schmidt BH, Parker JS, Mandel G, Hannon GJ: A biochemical approach to identifying microRNA targets. Proc Natl Acad Sci U S A. 2007, 104: 19291-19296. 10.1073/pnas.0709971104
Liu XH, Liu ZL, Sun M, Liu J, Wang ZX, De W: The long non-coding RNA HOTAIR indicates a poor prognosis and promotes metastasis in non-small cell lung cancer. BMC Cancer. 2013, 13: 464- 10.1186/1471-2407-13-464
Chen FJ, Sun M, Li SQ, Wu QQ, Ji L, Liu ZL, Zhou GZ, Cao G, Jin L, Xie HW, Wang CM, Lv J, De W, Wu M, Cao XF: Upregulation of the long non-coding RNA HOTAIR promotes esophageal squamous cell carcinoma metastasis and poor prognosis. Mol Carcinog. 2013, 52: 908-915. 10.1002/mc.21944
Bhan A, Hussain I, Ansari KI, Kasiri S, Bashyal A, Mandal SS: Antisense transcript long noncoding RNA (lncRNA) HOTAIR is transcriptionally induced by estradiol. J Mol Biol. 2013, 425: 3707-3722. 10.1016/j.jmb.2013.01.022
Barsyte-Lovejoy D, Lau SK, Boutros PC, Khosravi F, Jurisica I, Andrulis IL, Tsao MS, Penn LZ: The c-Myc oncogene directly induces the H19 noncoding RNA by allele-specific binding to potentiate tumorigenesis. Cancer Res. 2006, 66: 5330-5337. 10.1158/0008-5472.CAN-06-0037
Tao J, Zhao X, Tao J: c-MYC-miRNA circuitry: A central regulator of aggressive B-cell malignancies. Cell Cycle. 2014, 13: 191-198. 10.4161/cc.27646
Yang F, Xue X, Bi J, Zheng L, Zhi K, Gu Y, Fang G: Long noncoding RNA CCAT1, which could be activated by c-Myc, promotes the progression of gastric carcinoma. J Cancer Res Clin Oncol. 2013, 139: 437-445. 10.1007/s00432-012-1324-x
Tsai MC, Manor O, Wan Y, Mosammaparast N, Wang JK, Lan F, Shi Y, Segal E, Chang HY: Long noncoding RNA as modular scaffold of histone modification complexes. Science. 2010, 329: 689-693. 10.1126/science.1192002
Vrba L, Muñoz-Rodríguez JL, Stampfer MR, Futscher BW: miRNA gene promoters are frequent targets of aberrant DNA methylation in human breast cancer. PLoS One. 2013, 8: e54398- 10.1371/journal.pone.0054398
Wang J, Liu X, Wu H, Ni P, Gu Z, Qiao Y, Chen N, Sun F, Fan Q: CREB up-regulates long non-coding RNA, HULC expression through interaction with microRNA-372 in liver cancer. Nucleic Acids Res. 2010, 38: 5366-5383. 10.1093/nar/gkq285
Boll K, Reiche K, Kasack K, Mörbt N, Kretzschmar AK, Tomm JM, Verhaegh G, Schalken J, von Bergen M, Horn F, Hackermüller J: MiR-130a, miR-203 and miR-205 jointly repress key oncogenic pathways and are downregulated in prostate carcinoma. Oncogene. 2013, 32: 277-285. 10.1038/onc.2012.55
Zhang X, Huang L, Zhao Y, Tan W: Downregulation of miR-130a contributes to cisplatin resistance in ovarian cancer cells by targeting X-linked inhibitor of apoptosis (XIAP) directly. Acta Biochim Biophys Sin (Shanghai). 2013, 45: 995-1001. 10.1093/abbs/gmt113
Acunzo M, Visone R, Romano G, Veronese A, Lovat F, Palmieri D, Bottoni A, Garofalo M, Gasparini P, Condorelli G, Chiariello M, Croce CM: miR-130a targets MET and induces TRAIL-sensitivity in NSCLC by downregulating miR-221 and 222. Oncogene. 2012, 31: 634-642.
The authors thank the Eastern Hepatobiliary Surgical Hospital and Institute, The Second Military University, Shanghai for their generous help. This work was supported by National Natural Science Foundation of China (grant no. 81272747) and Doctorial innovation fund of Shanghai Jiaotong University School of Medicine. The funding sources had no role in the study design, in the collection, analysis and interpretation of data; in the writing of the manuscript and in the decision to submit the manuscript for publication.
The authors declare that they have no competing interests.
MZM conceived of the study and participated in its design and coordinated and helped to draft the manuscript. MZM, CXL, YZ and MZW performed the experiments. MDZ and YYQ participated in the design of the study and performed the statistical analysis. MZM, ZWQ and WG wrote the paper. All authors read and approved the final manuscript.
Ming-zhe Ma, Chun-xiao Li, Yan Zhang contributed equally to this work.
Electronic supplementary material
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.