STAT3-mediated upregulation of lncRNA HOXD-AS1 as a ceRNA facilitates liver cancer metastasis by regulating SOX4
Molecular Cancervolume 16, Article number: 136 (2017)
Several of the thousands of human long noncoding RNAs (lncRNAs) have been functionally characterized, yet their potential involvement in hepatocellular carcinoma (HCC) remains poorly understood.
LncRNA-HOXD-AS1 was identified by microarray and validated by real-time PCR. The clinicopathological significance of HOXD-AS1 was analyzed by Kaplan-Meier method. Chromatin immunoprecipitation was conducted to examine the mechanism of HOXD-AS1 upregulation. The role of HOXD-AS1 in HCC cells was assessed both in vitro and in vivo. ceRNA function of HOXD-AS1 was evaluated by RNA immunoprecipitation and biotin-coupled miRNA pull down assays.
In this study, we found that HOXD-AS1 was significantly upregulated in HCC tissues. Clinical investigation demonstrated high expression level of HOXD-AS1 was associated with poor prognosis and high tumor node metastasis stage of HCC patients, and was an independent risk factor for survival. Moreover, our results revealed that STAT3 could specifically interact with the promoter of HOXD-AS1 and activate HOXD-AS1 transcription. Knockdown of HOXD-AS1 significantly inhibited migration and invasion of HCC cells in vitro and distant lung metastasis in vivo. Additionally, HOXD-AS1 was enriched in the cytoplasm, and shared miRNA response elements with SOX4. Overexpression of HOXD-AS1 competitively bound to miR-130a-3p that prevented SOX4 from miRNA-mediated degradation, thus activated the expression of EZH2 and MMP2 and facilitated HCC metastasis.
In summary, HOXD-AS1 is a prognostic marker for HCC patients and it may play a pro-metastatic role in hepatocarcinogenesis.
Hepatocellular carcinoma (HCC) is one of the most frequently diagnosed malignancies and the third most common cause of cancer-related death worldwide. Despite the use of innovative therapeutic strategies for HCC, survival rate is still poor for HCC patients. It is mainly due to the high rates of recurrence and metastasis after surgical resection [1, 2]. Although great progress has been achieved in the identification of protein-coding genes or small noncoding RNAs, i.e., microRNAs, in development and progression of HCC [3, 4], little is known about the roles of long noncoding RNAs (lncRNAs) in hepatocarcinogenesis.
LncRNAs are a class of noncoding RNA transcripts that have recently attracted great research interest. So far, a large range of functions have been attributed to lncRNAs, such as control of muscle differentiation , reprogramming of induced pluripotent stem cells , and modulation of cell apoptosis and invasion . Indeed, accumulating evidence indicates that alteration and dysfunction of lncRNAs have been shown to result in aberrant genes expression that promotes tumor formation, progression and metastasis of many cancer types [8,9,10,11,12,13]. Several recent reports have described that lncRNAs are important cis- or trans-regulators of genes activities through a variety of mechanisms. In the nucleus, lncRNAs can function as scaffolds to bring proteins to form ribonucleoprotein complexes and as guides to recruit chromatin-modifying complexes to target genes [14,15,16]. In the cytoplasm, lncRNAs also can function as competing endogenous RNAs (ceRNAs) by competitively binding miRNAs, thereby modulating the derepression of miRNAs targets [17, 18]. To date, only few studies have documented lncRNAs in HCC [19,20,21,22] and the underlying mechanisms remain largely unknown.
SOX4, a member of the SOX (SRY-related HMG-box) gene family that containing the highly conserved HMG-domain responsible for specific DNA binding was identified as a common transcription factor that is involved in many cancer progression [23,24,25]. However, its direct target genes that mediate cancer progression are not well defined.
In the current study, using gene expression profiling analysis we discovered a lncRNA HOXD-AS1 was upregulated in HCC and significantly correlated with poor prognosis of HCC patients. Meanwhile, we demonstrated that HOXD-AS1 promoted HCC metastasis by acting as a ceRNA to upregulate the expression of SOX4 and activated the expression of EZH2 and MMP2, two direct target genes of SOX4.
Human HCC cell lines SNU449, Huh28, SMMC-7721, Huh7, Hep3B, and human embryonic kidney cell line HEK-293 T were cultured in Dulbecco’s Modified Eagle’s Medium (DMEM), supplemented with 10% fetal bovine serum, 100U/ml penicillin, and 100 μg/ml streptomycin. All of cell lines were cultured in a humidified chamber with 5% CO2 at 37 °C. The cells were tested regularly for mycoplasma (R&D Systems’ new MycoProbe Mycoplasma Detection Kit).
Patients and specimens
Primary human HCC cancerous tissues (C) and corresponding adjacent noncancerous liver tissues (N) were acquired from the surgical specimen archives of Zhongshan Hospital, Shanghai, China. HCC was on the basis of CT, ultrasound or MRI characteristics and biochemistry (AFP serology and liver function enzymes), and was confirmed by histopathology, according to the American Association for the Study of Liver Diseases guidelines . The inclusion criteria were that the samples contained matched tumors (percentage of tumor cells >70%) and corresponding normal liver tissue (>5 cm laterally from the edge of the cancerous region), the patient had a single primary lesion, and no neoadjuvant therapy. The patients who had the history of other solid tumors, had the secondary liver cancer from other primary regions, or did not have the follow-up information were excluded from this study.
Tumor differentiation was graded using the Edmondson grading system. Clinical staging was performed according to the tumor-node-metastasis (TNM) staging system. The follow-up procedures and postoperative treatments were based on a uniform guideline and have been described previously . Overall survival (OS) was calculated from the date of surgery to the date of death. Data were censored at the last follow-up visit or at the time of a patient’s death without relapse. Ethical approval was obtained from the Zhongshan Hospital Research Ethics Committee, and written informed consent was obtained from each patient.
Total RNAs were isolated from the paired tissue samples of 14 HCC patients and purified using TRIzol reagent (Invitrogen, Carlsbad, CA) and RNeasy mini kit (Qiagen Inc, Valencia, CA) according to the manufacturer's protocol. Clinicopathological characteristics of the 14 HCC patients are provided in Additional file 1: Table S1. Following RNA isolation and cDNA synthesis, biotin-labeled cRNA was prepared by Enzo® BioArrayTM HighYieldTM RNA transcript label kit (Enzo life sciences, INC), then hybridized to the Affymetrix Gene Chip Human Genome U133 Plus 2.0 Array (Santa Clara, CA). Differentially expressed lncRNAs were analyzed as described previously .
5’and 3’ rapid amplification of cDNA ends (RACE) analyses
5’RACE and 3’RACE analyses were performed with 5 μg of total RNA. The SMARTer™ RACE cDNA kit (Clontech) was used according to the manufacturer’s instructions. The gene specific primers used for PCR are presented in Additional file 2: Table S2.
Nuclear fractionation was performed with a PARIS™ Kit (Ambion, Austin, TX). For nuclear fractionation, 1 × 107 cells (Huh7 or SMMC-7721) were collected and resuspended in the cell fraction buffer and incubated on ice for 10 min. After centrifugation, supernatant and nuclear pellet were preserved for RNA extraction using a cell disruption buffer according to the manufacturer’s instructions.
Fluorescence in situ hybridization analysis
Probe for HOXD-AS1 was presented in Additional file 2: Table S2. SMMC-7721 and Huh7 cells were used for RNA-FISH analysis. Cell suspension was diluted to 100 cells/μL and seeded in the autoclaved glass slides. Slides were treated with 0.2 mol/L HCl for 20 min at room temperature, and washed in 2 × SSC buffer for 5 min. Slides were then incubated in 1%NaSCN for 30 min at 80 °C, and washed in 2 × SSC buffer for 5 min. Slides were incubated in 4% pepsase (2500 ~ 3000U/mg) for 10 min at 37 °C, and washed in 2 × SSC buffer for 5 min, and then fixed in 4% paraformaldehyde for 10 min at room temperature, washed in 2 × SSC buffer for 5 min. Dried the slides and prehybridized with prehybridization buffer for 2 h at 50 °C. Hybridization using HOXD-AS1 probe was performed overnight at 37 °C, and slides were rinsed in 2 × SSC with 0.3%NP-40 (pH7.0 ~ 7.5) for 30 min at 72 °C, washed in 2 × SSC with 0.3%NP-40 for 20 min at room temperature, and counterstained with 4’-6’diamidino-2-phenylindole (DAPI) for 5 min. The images were acquired using a confocal microscope (Leica).
RNA immunoprecipitation (RIP) experiments were performed with a Magna RIP™ RNA-Binding Protein Immunoprecipitation Kit (Millipore, Billerica, MA) according to the manufacturer’s instructions. AGO2 antibody was used for RIP (Cell Signaling Technology, Beverly, MA). Co-precipitated RNAs were detected by real-time PCR. Total RNAs (input controls) and IgG controls were assayed simultaneously to demonstrate that the detected signals were the result of RNAs specifically binding to AGO2. The gene-specific primers used for detecting HOXD-AS1 and SOX4 were presented in Additional file 2: Table S2.
A chromatin immunoprecipitation (ChIP) assay was performed using the EZ-ChIPTM kit (Millipore, Billerica, MA) according to the manufacturer’s instructions. The following antibodies were utilized to immunoprecipitate crosslinked protein-DNA complexes: rabbit anti-STAT3 (10253-2-AP, Protein tech), rabbit anti-SOX4 (ARP38234, AVIVA) and normal rabbit IgG (12–370, Millipore). The immunoprecipitated DNA was purified for quantitative PCR analyses with specific primers. The primers were listed in Additional file 2: Table S2.
Pull down of biotin-coupled miRNA
Biotin was attached to the 3’-end of miR-130a-3p or negative control mimics. 1 × 106 HCC cells were transfected with 100 pmol Bi-miR-130a-3p, or negative control using Lipofectamine 2000 according to the manufacturer’s protocol. Forty-eight hours later, cells were pelleted at 1000 rpm. After washing twice with PBS, cell pellets were resuspended in 0.7 ml lysis buffer (5 mM MgCl2, 100 mM KCl, 20 mM Tris (pH7.5), 0.3% NP-40, 50U of RNase OUT (Invitrogen, USA)), complete protease inhibitor cocktail (Roche Applied Science, IN), and incubated on ice for 10 min. The cell lysate was isolated by centrifugation at 10,000 g for 10 min. miRNA biotin pull down experiments were performed according to previous reports . The level of HOXD-AS1 or SOX4 in the pull down of Biotin-miR-130a-3p or negative control was quantified by real-time PCR.
Small interfering (si) RNA duplexes were designed and synthesized by Genepharma (Shanghai, China). MiR-130a-3p mimics, inhibitors and corresponding negative control (NC) were synthesized by Ribobio (Guangzhou, China). The sequences used are shown in Additional file 2: Table S2. Cells were transfected using Lipofectamine 2000.
Total RNA was extracted from cells using Trizol reagent (Invitrogen) according to the manufacturer’s instruction. A total of 1 μg of RNA was subjected to reverse transcription using HiScript II Q RT SuperMix for qPCR (+gDNA wiper) (Vazyme, Nanjing, China) and real-time PCR was carried out using AceQ qPCR SYBR Green Master Mix (Vazyme, Nanjing, China). The PCR reaction conditions were as follows: 95 °C for 15 s followed by 40 cycles of 95 °C for 5 s and 60 °C for 30s. The expression levels were normalized against those of the internal reference gene β-actin, the relative expression levels were determined by the following equation: 2-ΔΔCt (ΔCt = ΔCttarget - ΔCtβ-actin). A list of primers used for real-time PCR experiments were in Additional file 2: Table S2.
Western blot analysis
Proteins from cell lysates were prepared in 1 × sodium dodecyl sulfate buffer, separated by sodium dodecyl sulfate–polyacrylamide gel electrophoresis (SDS-PAGE) and transferred to a nitrocellulose membrane (Bio-Rad, Hercules, CA, USA). The membranes were blocked with 5% non-fat milk and incubated with the appropriate antibody. Antigen-antibody complex was detected with enhanced chemiluminescence reagents (Pierce, Rockford, IL). Antibodies used in this study were shown in Additional file 3: Table S3.
In vitro assays for migration and invasion
For transwell migration assay, 5 × 104 cells were plated in the top chamber of each insert (BD Biosciences, NJ, USA) with a non-coated membrane. For invasion assay, 8 × 104 cells were placed in the upper chamber of each Matrigel-coated insert (BD Biosciences, Billerica, MA, USA). After 18 h of incubation at 37 °C, cells that migrated or invaded were fixed and stained in dye solution containing 0.1% crystal violet and 20% methanol. The number of cells that had migrated or invaded was counted and imaged using an IX71 inverted microscope (Olympus Corp, Tokyo, Japan).
In vivo assay for metastasis
For in vivo metastasis assays, 3 × 106 Huh7 cells infected with shRNA-HOXD-AS1 or empty vector were suspended in 300 μL of serum-free DMEM per male BALB/c-nu/nu mice, and injected into nude mice through the tail vein (12 mice per group). After six weeks, mice were sacrificed, and their lungs were dissected, fixed with phosphate-buffered neutral formalin, and prepared for standard histological examination. Mice were handled and housed according to protocols approved by the Shanghai Medical Experimental Animal Care Commission.
Luciferase reporter assay
HEK-293 T cells were seeded in 96-well plates at a density of 5000 cells per well. After 24 h, the cells were transiently transfected with a mixture of 5 ng of pRL-CMV Renilla luciferase reporter, 50 ng of the firefly luciferase reporter, and 5pmol small RNA (siRNA or miRNA mimics). After 48 h, luciferase activity was measured using the dual-luciferase reporter assay system (Promega, Madison, WI, USA).
Results are presented as the means ± SD. Data were analyzed using Student’s t-test (two-tailed, with p < 0.05 considered significant) unless otherwise specified (paired t-test, χ2 test or Spearman correlation). Cumulative survival was evaluated using the Kaplan-Meier method (log-rank test). All the statistical analyses were performed using the SPSS 16.0 software (SPSS, Inc., Chicago, IL).
HOXD-AS1 is upregulated and positively associated with poor prognosis in HCC patients
To uncover the potential involvement of lncRNAs in HCC, we performed global gene expression analysis for total RNAs from the paired samples of cancerous tissues (C) and corresponding adjacent noncancerous liver tissues (N) of 14 HCC patients with Affymetrix GeneChip Human Genome U133 Plus 2.0 Array (Santa Clara, CA) (GEO Submission: GSE84402) . In total, we identified 26 lncRNAs that were significantly differentially expressed (fold change >22 or <2-2, p < 0.05) (Fig. 1a, Additional file 4: Table S4). In all 26 dysregulated lncRNAs, AW340112 (termed HOXD-AS1) was most upregulated in HCC tissues (fold change >8, Additional file 4: Table S4), it was expressed at a significantly high level in most of HCC tissues according to probe signal calculation (Fig. 1b, Additional file 5: Figure S1). We next assessed the expression of HOXD-AS1 in 89 paired samples of HCC tissues and corresponding adjacent noncancerous liver tissues. The results showed that HOXD-AS1 was significantly increased in HCC tissues, the upregulated expression of HOXD-AS1 was observed in 73% (65/89) of HCC (p < 0.0001) (Fig. 1c).
To further investigated the clinicopathological and prognostic significance of HOXD-AS1 levels in HCC patients, 120 HCC tissues (all of which had detailed clinical and follow-up data available) were examined. Our results showed that upregulation of HOXD-AS1 was significant correlated with TNM stage in HCC patients (Table 1). Furthermore, Kaplan-Meier analysis revealed that high level of HOXD-AS1 was associated with poor overall survival (OS) in HCC patients (Fig. 1d; p = 0.0179). Multivariate regression revealed that the expression level of HOXD-AS1 was an independent and significant factor for OS (HR: 0.552; 95% CI: 0.321–0.950; p = 0.032) (Table 2).
HOXD-AS1 is transcribed in antisense orientation of the protein-coding genes HOXD1. We obtained the full-length sequence of HOXD-AS1 by RACE analyses, it is composed of 3956 nucleotides and consists of four exons and three introns with 3’ polyadenylate tail (Additional file 6: Figure S2). Furthermore, subcellular fractionation analysis showed that HOXD-AS1 was predominantly abundant in the cytoplasm of HCC cells (Fig. 1e). In situ hybridization analysis also showed that HOXD-AS1 was located mainly in the cytoplasm (Fig. 1f).
HOXD-AS1 is regulated by the transcription factor STAT3
To explore the mechanism of HOXD-AS1 upregulation in HCC, we applied online bioinformatical software programs JASPAR (http://jaspar.genereg.net/cgi-bin/jaspar_db.pl) to analyze promoter region of HOXD-AS1 and 13 potential sites of STAT3 binding were found. Next, we examined the mRNA expression level of STAT3 in 68 HCC specimens. We found that expression level of STAT3 was significantly upregulated, and was positively correlated with HOXD-AS1 in HCC tissues (Fig. 2a-b). To further investigate the correlation between HOXD-AS1 and STAT3, we detected the transcriptional level of HOXD-AS1 and the protein level of phosphorylated STAT3 in the same group of HCC tissues. We found that active STAT3 was also positively correlated with HOXD-AS1 in 30 HCC tissues (Additional file 7: Figure S3A-B). Then, we measured mRNA levels of HOXD-AS1 in the case of STAT3 knockdown. The results revealed that STAT3 knockdown led to a significant decrease of HOXD-AS1 in Hep3B and Huh7 cells (Fig. 2d).
To confirm HOXD-AS1 was a transcriptional target of STAT3, serial truncations of HOXD-AS1 promoter were cloned into pGL3-basic vector, and luciferase activity was measured after transfection of these constructs into HEK-293 T cells. The highest activities were associated with -2000 nt ~ -1159 nt and -1158 nt ~ -746 nt (Fig. 2e), indicating that these two fragments contained regulatory elements and they were critical for the transcription of HOXD-AS1. Next, we cotransfected luciferase reporter and siRNA against STAT3 into HEK-293 T cells. The results showed that STAT3 knockdown significantly reduced luciferase activity of -1158 nt ~ -746 nt fragment, however, luciferase activity of reporter containing -2000 nt ~ -1159 nt was not affected (Fig. 2f). These results indicate that the region between -1158 bp to -746 bp on the HOXD-AS1 promoter is responsible for STAT3-mediated activation of HOXD-AS1. Sequence analysis of -1158 nt ~ -746 nt fragment uncovered three putative STAT3 binding sites, located at -828 nt ~ -818 nt (site 1), -938 nt ~ -928 nt (site 2), and -1121 nt ~ -1111 nt (site 3). To determine which STAT3 binding sites were responsive to STAT3-mediated transcriptional activation of HOXD-AS1. ChIP assay was performed. Our results showed that STAT3 could directly bind to the site 2 on the HOXD-AS1 promoter in HCC cell lines (Fig. 2g). These findings indicate that the upregulation of HOXD-AS1 is mediated by STAT3 in HCC.
HOXD-AS1 promotes invasion and metastasis of HCC cells in vitro and in vivo
To investigate the function of HOXD-AS1 on cell biological behavior, we first determined the expression levels of HOXD-AS1 in various HCC cell lines. The results showed that HOXD-AS1 expression was relatively high in SMMC-7721 and Huh7 cells, whereas it was relatively low in Huh28 and SNU449 cells (Additional file 8: Figure S4A). We therefore constructed a lentivirus vector harboring shRNA-HOXD-AS1 and established two stable knockdown cell lines in SMMC-7721 and Huh7 (Additional file 8: Figure S4B). Meanwhile, we stably overexpressed HOXD-AS1 in SNU449 and Huh28 cell lines (Additional file 8: Figure S4C). Transwell assays showed that migratory and invasive abilities of SMMC-7721 and Huh7 cells were significantly reduced when HOXD-AS1 was decreased (Fig. 3a-b). In contrast, exogenous expression of HOXD-AS1 dramatically promoted migration and invasion in Huh28 and SNU449 cells (Fig. 3c-d). However, HOXD-AS1 did not affect HCC cell proliferation and apoptosis (data not shown). To explore the effect of HOXD-AS1 on tumor metastasis in vivo, we implanted either control or Huh7 cells with stable HOXD-AS1 knockdown into nude mice via the tail vein. The results showed that HOXD-AS1 knockdown significantly decreased the number of metastatic lung nodules compared with the control cells (n = 12, p = 0.035) (Fig. 3e).
HOXD-AS1 is involved in migration of HCC cells through regulating SOX4
To identify the mechanism of HOXD-AS1 in HCC cell metastasis, we analyzed the co-expression patterns of HOXD-AS1 and core transcript factors (TFs) in the 14 HCC tissues used in global gene expression analysis. Five TFs were highly correlated with the level of HOXD-AS1, which are SOX4, E2F7, E2F8, NFYA, and TCF19 (Fig. 4a, Additional file 9: Figure S5A). We therefore hypothesized that HOXD-AS1 might regulate these TFs in HCC cells. To test this hypothesis, we first investigated whether HOXD-AS1 knockdown affected expression levels of these five TFs. Our results showed that the expression level of SOX4 was downregulated after HOXD-AS1 reduction in Huh7 cells (Fig. 4b), however, expression levels of E2F7, E2F8, NFYA, and TCF19 were not affected significantly (Additional file 9: Figure S5B). We then examined levels of five TFs in SNU449 cells with HOXD-AS1 overexpression. As a result, we found that only SOX4 significantly upregulated when HOXD-AS1 was ectopically overexpressed in SNU449 cells (Fig. 4c, Additional file 9: Figure S5C). To further determine functional association between HOXD-AS1 and SOX4, cell migration after overexpression or knockdown of SOX4 was examined. We found that the suppression of migratory ability by HOXD-AS1 knockdown was reversed when simultaneously overexpressing SOX4 in Huh7 cells (Fig. 4d). In addition, knockdown of SOX4 abrogated the effects of HOXD-AS1 overexpression on SNU449 cell migration (Fig. 4e). These results indicate that SOX4 is necessary for HOXD-AS1 to regulate HCC metastasis.
HOXD-AS1 regulates SOX4 expression through a miRNA-dependent mechanism
Because HOXD-AS1 transcripts were abundant in the cytoplasm and the expression levels of HOXD-AS1 and SOX4 changed synchronously in HCC cells, we hypothesized that HOXD-AS1 may have the role as a ceRNA to regulate the expression of SOX4. To test this hypothesis, we first performed RNA immunoprecipitation (RIP) with an antibody against Argonaute (AGO2) in SMMC-7721 and Huh7 cells. We observed an enrichment of HOXD-AS1 and SOX4 with AGO2 antibody (Fig. 5a). These results indicated that HOXD-AS1 and SOX4 were recruited to AGO2 related RNA-induced silencing complexes (RISC) and might functionally interact with miRNAs. Additionally, we also evaluated the expression of HOXD-AS1 and SOX4 in Dicer-deficient Huh7 cells. We found that Dicer knockdown partially rescued the reduction of SOX4 by HOXD-AS1 knockdown (Fig. 5b). These results suggest that miRNAs may play essential roles in HOXD-AS1 mediated regulation of SOX4.
To investigate the potential miRNAs, we used bioinformatic tool miRanda to search for miRNAs that target the transcripts of HOXD-AS1 and SOX4. According to the prediction, miR-130a-3p, miR130b-3p, and miR454-3p have the same putative binding sites mapped to HOXD-AS1 and SOX4. We next analyzed expression status of three miRNAs in TCGA HCC cohort. The results showed that only miR-130a-3p was greatly reduced (Additional file 10: Figure S6A). In addition, previous study reported that miR-130a-3p could suppress cell migration and invasion in HCC cells . These results suggest that miR-130a-3p may be a critical regulatory miRNA for HOXD-AS1 and SOX4.
To confirm HOXD-AS1 and SOX4 were regulated by miR-130a-3p, we constructed luciferase reporters containing wild type (WT) and mutated (MUT) putative binding sites of HOXD-AS1 or SOX4 transcripts, respectively (Additional file 10: Figure S6B-C). We found that reporter with wild type binding sites of HOXD-AS1 or SOX4 showed a markedly lower luciferase activity in HEK-293 T cells with overexpression of miR-130a-3p, however, we did not observe variations of luciferase activity in reporter with the mutated binding sites (Fig. 5c). To further confirm miRNAs-SOX4 interaction, protein level of SOX4 was evaluated in HCC cells transfected with miR-130a-3p. We found that miR-130a-3p suppressed endogenous SOX4 in HCC cells (Fig. 5d). Additionally, to validate the direct binding ability of miR-130a-3p on HOXD-AS1 and SOX4, we performed biotin-coupled miRNA pull down to capture HOXD-AS1 and SOX4 using streptavidin-coated beads from cells transfected with 3’end- biotinylated miR-130a-3p (Bi-has-miR-130a-3p). Control samples were transfected with a biotinylated C.elegans miRNA (Bi-cel-miR67). After 48 h transfection, we found that HOXD-AS1 was enriched to more than fourfold in Bi-miR-130a-3p pull down of HCC cells as compared to Bi-cel-miR67 transfection, and also we observed more than fivefold enrichment of SOX4 in the miR-130a-3p captured fraction compared with the negative control (Fig. 5e). Furthermore, we detected competitive binding activities of HOXD-AS1 and SOX4 to miR-130a-3p in SNU449 cells with HOXD-AS1 overexpression and Huh7 cells with HOXD-AS1 knockdown. Results showed that overexpression of HOXD-AS1 in SNU449 cells led to the decreased enrichment of SOX4 transcripts on miR-130a-3p. However, HOXD-AS1 knockdown in Huh7 cells caused a significant increase in the recruitment of SOX4 to miR-130a-3p (Additional file 11: Figure S7A-B).
To further explore the ceRNA role of HOXD-AS1, the expression levels of SOX4 after overexpression or knockdown of miR-130a-3p were examined. As shown in Fig. 5f, overexpression of miR-130a-3p could reduce the levels of SOX4, which rescued by the overexpression of HOXD-AS1. In addition, miR-130a-3p inhibitors nullified the inhibitory effect of HOXD-AS1 on the expression of SOX4 (Fig. 5f). These results demonstrate that HOXD-AS1 regulates SOX4 partly in miR-130a-3p dependent manner by sharing miRNAs response elements with SOX4.
SOX4 is required for expression of MMP2 and EZH2 mediated by HOXD-AS1
Because HOXD-AS1 promoted metastasis of HCC cells, it is important to explore whether some genes involved in this progression. We assessed the expression levels of 15 metastasis-related genes (Additional file 3: Table S3) in Huh7 cells with HOXD-AS1 knockdown. Our results showed that inhibition of HOXD-AS1 suppressed MMP13, MAPK1, HDAC1, EZH2, and MMP2 at both mRNA and protein levels compared with negative control (Fig. 6a). Conversely, overexpression of HOXD-AS1 increased mRNA and protein levels of these five genes in SNU449 cells (Fig. 6b). Next, we detected functional role of SOX4 on HOXD-AS1 regulation of these five genes. We found that knockdown of SOX4 strongly decreased expression levels of EZH2 and MMP2 that were enhanced by HOXD-AS1 overexpression (Fig. 6c), while the expression of MMP13, MAPK1, and HDAC1 were not affected (data not shown). Given that SOX4 is a transcription factor, we decided to investigate whether EZH2 and MMP2 is the target genes of SOX4. Chromatin immunoprecipitation (ChIP) with antibody against SOX4 or IgG was performed with Huh28 cell extracts of stable HOXD-AS1 overexpression or control. As expected, promoter regions of MMP2 and EZH2 were recruited to SOX4 in the control cells. Notably, the recruitment of the promoters of MMP2 and EZH2 to SOX4 was markedly increased in Huh28 cells with HOXD-AS1 overexpression (Fig. 6d). These results suggest that EZH2 and MMP2 are the direct target genes of SOX4, and HOXD-AS1 can mediate the expression of EZH2 and MMP2 through SOX4.
In the present study, we identified a number of lncRNAs that are aberrantly expressed in human HCC. Among them, HOXD-AS1 was most upregulated in HCC. High level of HOXD-AS1 expression was associated with significantly reduced overall survival, it may represent an independent prognostic biomarker in patients with HCC. LncRNAs with differential expression in cancer were correlated with good or bad prognosis, making them promising prognostic biomarkers . For example, GAS5 was downregulated in several cancers [33,34,35], and low level of GAS5 indicated a poor prognosis in HCC . Upregulation of HOTAIR was shown to be a marker of poor prognosis in a number of cancers, including HCC [36,37,38,39,40,41]. Moreover, increased HOTAIR was a prognostic biomarker of tumor recurrence following liver transplantation . HOTTIP, HOXA13, and HEIH also act as independent prognostic factors associated with recurrence in HCC [43, 44], suggesting that lncRNAs have considerable prognostic potential in HCC.
Here, we identified the mechanism responsible for HOXD-AS1 upregulation in HCC cells. We first found that STAT3 could specific interact with the HOXD-AS1 promoter through the binding site located -938 nt ~ -928 nt and illustrated the mechanism by which STAT3 upregulated the expression of HOXD-AS1. The present results noticed that the region between -2000 nt to -746 nt on the HOXD-AS1 promoter contained regulatory elements for the transcription of HOXD-AS1. However, only the region between -1158 nt ~ -746 nt had the STAT3 binding site and contributed to mediate transcription of HOXD-AS1. This result indicates that there may be additional transcriptional factors in the region between -2000 to -1159 that regulate the transcription of HOXD-AS1.
It has been described that lncRNAs play an important role in hepatocarcinogenesis [45, 46]. In this study, we have shown that HOXD-AS1 can promote migration and invasion of HCC cells in vitro and contribute to distant lung metastasis in vivo. To our knowledge, this is the first study to report that HOXD-AS1 regulates cellular metastasis in HCC. In bladder cancer cells, knockdown of HOXD-AS1 could suppress cell proliferation and migration, and increased the rate of apoptotic cell . Additionally, HOXD-AS1 was involved in angiogenesis and inflammation and controlled cell differentiation in neuroblastoma . These findings further support the potential pro-oncogenic role of HOXD-AS1 in cancers. However, molecular mechanisms of HOXD-AS1 in pro-oncogenesis have not been clarified in human cancers, including HCC.
Based on global genome microarray data, we integrated gene co-expression patterns and identified positive correlation between HOXD-AS1 and SOX4 in HCC tissues. Synchronous change in expression levels of HOXD-AS1 and SOX4 was also detected in HCC cell lines. Previous studies reported that lncRNAs could crosstalk with protein-coding genes in a miRNA-dependent manner in HCC progression [13, 49, 50]. In this study, we revealed that HOXD-AS1 could function as a ceRNA that sponge miRNA130a-3p to protect SOX4 against degradation. Importantly, Liao et al. observed SOX4 potentiates metastasis in HCC . The identification of a model of miRNA/lncRNA interaction may promote the understanding of the underlying mechanism of HCC metastasis.
We have also shown that 5 metastasis-related genes MMP13, MAPK1, HDAC1, EZH2, and MMP2 are required for HOXD-AS1 in the role of HCC metastasis. Among these five genes, we further found that EZH2 and MMP2 were the direct target genes of SOX4. Our results demonstrate that HOXD-AS1 is important for induction of SOX4, which result in enhanced expression of EZH2 and MMP2 and, in turn, lead to HCC metastasis. Additionally, the potential mechanism of MMP13, MAPK1, and HDAC1 regulated by HOXD-AS1 will be further analyze.
In conclusion, HOXD-AS1 is regulated by the transcriptional factor STAT3, it significantly upregulated in HCC and can be used as a prognosis biomarker for HCC patients. HOXD-AS1 functions as a ceRNA that competitively binds to miR-130a-3p, then upregulates SOX4 and promotes HCC cell metastasis. These findings provide a new mechanism for understanding HCC metastasis and HOXD-AS1 may be a potential candidate in the prevention and treatment of HCC.
In this study, we reported that a number of lncRNAs are differentially expressed in HCC tissues. One of these lncRNAs, HOXD-AS1 was markedly upregulated in HCC. The high expression level of HOXD-AS1 was associated with poor prognosis and high tumor node metastasis stage in HCC patients and it may represent an independent prognostic biomarker in patients with HCC. Moreover, our results suggested that the transcription factor STAT3 could combine to the promoter of HOXD-AS1 and activate the transcription of HOXD-AS1. We further found that HOXD-AS1 facilitated HCC metastasis in vitro and in vivo. Furthermore, we had evidenced that HOXD-AS1 shared miRNA response elements with SOX4. Overexpression of HOXD-AS1 competitively bound to miR-130a-3p that prevented SOX4 from miRNA-mediated degradation, thus activated the expression of EZH2 and MMP2 and facilitated HCC metastasis. In conclusion: HOXD-AS1 is a prognostic marker for HCC patients and it may play a pro-metastatic role in hepatocarcinogenesis.
Competing endogenous RNAs
Dulbecco’s modified eagle’s medium
E2F transcription factor 7
E2F transcription factor8
Enhancer of zest homolog 2
Fluorescence in situ hybridization
Gene expression omnibus
Histone deacetylase 1
Long noncoding RNA
Mitogen-activated protein kinase 1
Matrix metallopeptidase 13
Matrix metallopeptidase 2
Nuclear transcription factor Y subunit alpha
Rapid amplification of cDNA ends
RNA-induced silencing complexes
Short hairpin RNA
Signal transducer and activator of transcription 3
Transcription factor 19
Forner A, Llovet JM, Bruix J. Hepatocellular carcinoma. Lancet. 2012;379:1245–55.
Portolani N, Coniglio A, Ghidoni S, Giovanelli M, Benetti A, Tiberio GA, Giulini SM. Early and late recurrence after liver resection for hepatocellular carcinoma: prognostic and therapeutic implications. Ann Surg. 2006;243:229–35.
Ding J, Huang S, Wu S, Zhao Y, Liang L, Yan M, Ge C, Yao J, Chen T, Wan D, et al. Gain of miR-151 on chromosome 8q24.3 facilitates tumour cell migration and spreading through downregulating RhoGDIA. Nat Cell Biol. 2010;12:390–9.
Yao J, Liang L, Huang S, Ding J, Tan N, Zhao Y, Yan M, Ge C, Zhang Z, Chen T, et al. MicroRNA-30d promotes tumor invasion and metastasis by targeting Galphai2 in hepatocellular carcinoma. Hepatology. 2010;51:846–56.
Cesana M, Cacchiarelli D, Legnini I, Santini T, Sthandier O, Chinappi M, Tramontano A, Bozzoni I. A long noncoding RNA controls muscle differentiation by functioning as a competing endogenous RNA. Cell. 2011;147:358–69.
Loewer S, Cabili MN, Guttman M, Loh YH, Thomas K, Park IH, Garber M, Curran M, Onder T, Agarwal S, et al. Large intergenic non-coding RNA-RoR modulates reprogramming of human induced pluripotent stem cells. Nat Genet. 2010;42:1113–7.
Khaitan D, Dinger ME, Mazar J, Crawford J, Smith MA, Mattick JS, Perera RJ. The melanoma-upregulated long noncoding RNA SPRY4-IT1 modulates apoptosis and invasion. Cancer Res. 2011;71:3852–62.
Gupta RA, Shah N, Wang KC, Kim J, Horlings HM, Wong DJ, Tsai MC, Hung T, Argani P, Rinn JL, et al. Long non-coding RNA HOTAIR reprograms chromatin state to promote cancer metastasis. Nature. 2010;464:1071–6.
Prensner JR, Iyer MK, Sahu A, Asangani IA, Cao Q, Patel L, Vergara IA, Davicioni E, Erho N, Ghadessi M, et al. The long noncoding RNA SChLAP1 promotes aggressive prostate cancer and antagonizes the SWI/SNF complex. Nat Genet. 2013;45:1392–8.
Trimarchi T, Bilal E, Ntziachristos P, Fabbri G, Dalla-Favera R, Tsirigos A, Aifantis I. Genome-wide mapping and characterization of Notch-regulated long noncoding RNAs in acute leukemia. Cell. 2014;158:593–606.
Xing Z, Lin A, Li C, Liang K, Wang S, Liu Y, Park PK, Qin L, Wei Y, Hawke DH, et al. lncRNA directs cooperative epigenetic regulation downstream of chemokine signals. Cell. 2014;159:1110–25.
Yang L, Lin C, Jin C, Yang JC, Tanasa B, Li W, Merkurjev D, Ohgi KA, Meng D, Zhang J, et al. lncRNA-dependent mechanisms of androgen-receptor-regulated gene activation programs. Nature. 2013;500:598–602.
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:666–81.
Engreitz JM, Pandya-Jones A, McDonel P, Shishkin A, Sirokman K, Surka C, Kadri S, Xing J, Goren A, Lander ES, et al. The Xist lncRNA exploits three-dimensional genome architecture to spread across the X chromosome. Science. 2013;341:1237973.
Huarte M, Guttman M, Feldser D, Garber M, Koziol MJ, Kenzelmann-Broz D, Khalil AM, Zuk O, Amit I, Rabani M, et al. A large intergenic noncoding RNA induced by p53 mediates global gene repression in the p53 response. Cell. 2010;142:409–19.
Simon MD, Pinter SF, Fang R, Sarma K, Rutenberg-Schoenberg M, Bowman SK, Kesner BA, Maier VK, Kingston RE, Lee JT. High-resolution Xist binding maps reveal two-step spreading during X-chromosome inactivation. Nature. 2013;504:465–9.
Poliseno L, Salmena L, Zhang J, Carver B, Haveman WJ, Pandolfi PP. A coding-independent function of gene and pseudogene mRNAs regulates tumour biology. Nature. 2010;465:1033–8.
Salmena L, Poliseno L, Tay Y, Kats L, Pandolfi PP. A ceRNA hypothesis: the Rosetta stone of a hidden RNA language? Cell. 2011;146:353–8.
Chen ZZ, Huang L, Wu YH, Zhai WJ, Zhu PP, Gao YF. LncSox4 promotes the self-renewal of liver tumour-initiating cells through Stat3-mediated Sox4 expression. Nat Commun. 2016;7:12598.
Wang X, Sun W, Shen W, Xia M, Chen C, Xiang D, Ning B, Cui X, Li H, Li X, et al. Long non-coding RNA DILC regulates liver cancer stem cells via IL-6/STAT3 axis. J Hepatol. 2016;64:1283–94.
Zhang H, Xing Z, Mani SK, Bancel B, Durantel D, Zoulim F, Tran EJ, Merle P, Andrisani O. RNA helicase DEAD box protein 5 regulates Polycomb repressive complex 2/Hox transcript antisense intergenic RNA function in hepatitis B virus infection and hepatocarcinogenesis. Hepatology. 2016;64:1033–48.
Zhou CC, Yang F, Yuan SX, Ma JZ, Liu F, Yuan JH, Bi FR, Lin KY, Yin JH, Cao GW, et al. Systemic genome screening identifies the outcome associated focal loss of long noncoding RNA PRAL in hepatocellular carcinoma. Hepatology. 2016;63:850–63.
Castillo SD, Matheu A, Mariani N, Carretero J, Lopez-Rios F, Lovell-Badge R, Sanchez-Cespedes M. Novel transcriptional targets of the SRY-HMG box transcription factor SOX4 link its expression to the development of small cell lung cancer. Cancer Res. 2012;72:176–86.
Vervoort SJ, van Boxtel R, Coffer PJ. The role of SRY-related HMG box transcription factor 4 (SOX4) in tumorigenesis and metastasis: friend or foe? Oncogene. 2013;32:3397–409.
Bilir B, Osunkoya AO, Wiles WG, Sannigrahi S, Lefebvre V, Metzger D, Spyropoulos DD, Martin WD, Moreno CS. SOX4 is essential for prostate tumorigenesis initiated by PTEN ablation. Cancer Res. 2016;76:1112–21.
Bruix J, Sherman M, Practice Guidelines Committee AAftSoLD. Management of hepatocellular carcinoma. Hepatology. 2005;42:1208–36.
Gao Q, Wang XY, Qiu SJ, Yamato I, Sho M, Nakajima Y, Zhou J, Li BZ, Shi YH, Xiao YS, et al. Overexpression of PD-L1 significantly associates with tumor aggressiveness and postoperative recurrence in human hepatocellular carcinoma. Clin Cancer Res. 2009;15:971–9.
Kawaguchi A, Yajima N, Komohara Y, Aoki H, Tsuchiya N, Homma J, Sano M, Natsumeda M, Uzuka T, Saitoh A, et al. Identification and validation of a gene expression signature that predicts outcome in malignant glioma patients. Int J Oncol. 2012;40:721–30.
Lal A, Thomas MP, Altschuler G, Navarro F, O’Day E, Li XL, Concepcion C, Han YC, Thiery J, Rajani DK, et al. Capture of microRNA-bound mRNAs identifies the tumor suppressor miR-34a as a regulator of growth factor signaling. PLoS Genet. 2011;7:e1002363.
Jin H, Wang C, Jin G, Ruan H, Gu D, Wei L, Wang H, Wang N, Arunachalam E, Zhang Y, et al. Regulator of Calcineurin 1 Gene Isoform 4, Downregulated in hepatocellular carcinoma, prevents proliferation, migration, and invasive activity of cancer cells and growth of orthotopic tumors by inhibiting nuclear translocation of NFAT1. Gastroenterology. 2017. doi:10.1053/j.gastro.2017.05.045.
Liu Y, Li Y, Wang R, Qin S, Liu J, Su F, Yang Y, Zhao F, Wang Z, Wu Q. MiR-130a-3p regulates cell migration and invasion via inhibition of Smad4 in gemcitabine resistant hepatoma cells. J Exp Clin Cancer Res. 2016;35:19.
Yarmishyn AA, Kurochkin IV. Long noncoding RNAs: a potential novel class of cancer biomarkers. Front Genet. 2015;6:145.
Lu X, Fang Y, Wang Z, Xie J, Zhan Q, Deng X, Chen H, Jin J, Peng C, Li H, Shen B. Downregulation of gas5 increases pancreatic cancer cell proliferation by regulating CDK6. Cell Tissue Res. 2013;354:891–6.
Mourtada-Maarabouni M, Pickard MR, Hedge VL, Farzaneh F, Williams GT. GAS5, a non-protein-coding RNA, controls apoptosis and is downregulated in breast cancer. Oncogene. 2009;28:195–208.
Sun M, Jin FY, Xia R, Kong R, Li JH, Xu TP, Liu YW, Zhang EB, Liu XH, De W. Decreased expression of long noncoding RNA GAS5 indicates a poor prognosis and promotes cell proliferation in gastric cancer. BMC Cancer. 2014;14:319.
Huang L, Liao LM, Liu AW, Wu JB, Cheng XL, Lin JX, Zheng M. Overexpression of long noncoding RNA HOTAIR predicts a poor prognosis in patients with cervical cancer. Arch Gynecol Obstet. 2014;290:717–23.
Ishibashi M, Kogo R, Shibata K, Sawada G, Takahashi Y, Kurashige J, Akiyoshi S, Sasaki S, Iwaya T, Sudo T, et al. Clinical significance of the expression of long non-coding RNA HOTAIR in primary hepatocellular carcinoma. Oncol Rep. 2013;29:946–50.
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–25.
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.
Nakagawa T, Endo H, Yokoyama M, Abe J, Tamai K, Tanaka N, Sato I, Takahashi S, Kondo T, Satoh K. Large noncoding RNA HOTAIR enhances aggressive biological behavior and is associated with short disease-free survival in human non-small cell lung cancer. Biochem Biophys Res Commun. 2013;436:319–24.
Wu ZH, Wang XL, Tang HM, Jiang T, Chen J, Lu S, Qiu GQ, Peng ZH, Yan DW. Long non-coding RNA HOTAIR is a powerful predictor of metastasis and poor prognosis and is associated with epithelial-mesenchymal transition in colon cancer. Oncol Rep. 2014;32:395–402.
Yang Z, Zhou L, Wu LM, Lai MC, Xie HY, Zhang F, Zheng SS. Overexpression of long non-coding RNA HOTAIR predicts tumor recurrence in hepatocellular carcinoma patients following liver transplantation. Ann Surg Oncol. 2011;18:1243–50.
Quagliata L, Matter MS, Piscuoglio S, Arabi L, Ruiz C, Procino A, Kovac M, Moretti F, Makowska Z, Boldanova T, et al. Long noncoding RNA HOTTIP/HOXA13 expression is associated with disease progression and predicts outcome in hepatocellular carcinoma patients. Hepatology. 2014;59:911–23.
Yuan SX, Yang F, Yang Y, Tao QF, Zhang J, Huang G, Yang Y, Wang RY, Yang S, Huo XS, et al. Long noncoding RNA associated with microvascular invasion in hepatocellular carcinoma promotes angiogenesis and serves as a predictor for hepatocellular carcinoma patients’ poor recurrence-free survival after hepatectomy. Hepatology. 2012;56:2231–41.
Xie H, Ma H, Zhou D. Plasma HULC as a promising novel biomarker for the detection of hepatocellular carcinoma. Biomed Res Int. 2013;2013:136106.
Braconi C, Kogure T, Valeri N, Huang N, Nuovo G, Costinean S, Negrini M, Miotto E, Croce CM, Patel T. microRNA-29 can regulate expression of the long non-coding RNA gene MEG3 in hepatocellular cancer. Oncogene. 2011;30:4750–6.
Li J, Zhuang C, Liu Y, Chen M, Chen Y, Chen Z, He A, Lin J, Zhan Y, Liu L, et al. Synthetic tetracycline-controllable shRNA targeting long non-coding RNA HOXD-AS1 inhibits the progression of bladder cancer. J Exp Clin Cancer Res. 2016;35:99.
Yarmishyn AA, Batagov AO, Tan JZ, Sundaram GM, Sampath P, Kuznetsov VA, Kurochkin IV. HOXD-AS1 is a novel lncRNA encoded in HOXD cluster and a marker of neuroblastoma progression revealed via integrative analysis of noncoding transcriptome. BMC Genomics. 2014;15 Suppl 9:S7.
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–83.
Wang L, Guo ZY, Zhang R, Xin B, Chen R, Zhao J, Wang T, Wen WH, Jia LT, Yao LB, Yang AG. Pseudogene OCT4-pg4 functions as a natural micro RNA sponge to regulate OCT4 expression by competing for miR-145 in hepatocellular carcinoma. Carcinogenesis. 2013;34:1773–81.
Liao YL, Sun YM, Chau GY, Chau YP, Lai TC, Wang JL, Horng JT, Hsiao M, Tsou AP. Identification of SOX4 target genes using phylogenetic footprinting-based prediction from expression microarrays suggests that overexpression of SOX4 potentiates metastasis in hepatocellular carcinoma. Oncogene. 2008;27:5578–89.
This work was supported by grants from National Key Basic Research Program of China (973 Program: 2015CB553905), National Natural Science Foundation of China (81201626, 81301759, 81402278, 81421001, 81472676 and 81672839), National High Technology Research and Development Program (863 Program) of China (2015AA020401), State Key Program of National Natural Science of China (81530077), Project of Special Research Fund for Health (201402003), Shanghai Natural Science Foundation of China (16ZR1434700), and Shanghai Jiao Tong University School of Medicine (YG2015QN34).
Availability of data and materials
All data generated or analysed during this study are included in this published article and its Additional files.
HW, XH designed the study, analyzed and interpreted the data, and wrote the paper; HW, XH, and JH performed most of the experiments with crucial help from LC, NW, XD, HJ, NW, CW, and FZ; JF and MY made critical advices to the conception and design of the study; XY and JF provided HCC patients’ samples and/or clinical data. WQ designed the study, supervised the project, and revised the manuscript. All authors vouch for the respective data and analysis, approved the final version and agreed to submit the manuscript.
The authors declare that they have no competing interests.
Consent for publication
We have obtained patients’ consent for publication.
Ethics approval and consent to participate
The research protocol was reviewed and approved by the Human Research Ethics Committee of Zhongshan Hospital, and written informed consent was obtained from each patient included in the study.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Clinicopathological characteristics of 14 HCC patients in global gene expression analysis. (DOCX 18 kb)
Oligonucleotide sequences in this study. (DOCX 19 kb)
Antibodies used in this study. (DOCX 15 kb)
The 26 differentially expressed lncRNAs in microarray analysis. (DOCX 18 kb)
The fold change of HOXD-AS1 expression in 14 HCC tissues. (A) probe: 239182_a. (B) probe: 242042_s_at. (TIF 998 kb)
RACE analysis of the full length of HOXD-AS1. (A) A garose gel electrophoresis of PCR products from 5’-RACE and 3’-RACE analysis of HOXD-AS1. (B) The nucleotide sequence of the full-length HOXD-AS1. (TIF 5631 kb)
Analysis for HOXD-AS1 expression and phosphorylated STAT3 in HCC tissues. (A) The transcriptional expression level of HOXD-AS1 and phosphorylated STAT3 in 30 HCC tissues. (B) Correlation analysis between p-STAT3 (x) and HOXD-AS1 (y) in 30 HCC tissues (R = 0.53, p = 0.0023). (TIF 2920 kb)
The expression level of HOXD-AS1 in HCC cells. (A) Real-time PCR analysis of HOXD-AS1 expression in 8 different HCC cell lines and normal liver cell L02. (B) Expression of HOXD-AS1 was quantified by Real-time PCR after knockdown of HOXD-AS1 in SMMC-7721 and Huh7 cells. (C) Expression levels of HOXD-AS1 in Huh28 and SNU449 cells that had been stably transfected with lentivirus encoding HOXD-AS1. (TIF 1314 kb)
Expression correlation between HOXD-AS1 and transcription factors. (A) Expression correlation between HOXD-AS1 (x) and E2F7, E2F8, NFYA, and TCF19 (y) in cancerous tissues of the 14 HCC patients used in global gene expression analysis. Data are depicted as probe signal value of global gene expression microarray. (B) Real-time analysis for changes of E2F7, E2F8, NFYA, and TCF19 in Huh7 cells with HOXD-AS1 knockdown. (C) Real-time PCR analysis for changes of E2F7, E2F8, NFYA, and TCF19 in SNU449 cells with HOXD-AS1 overexpression. (TIF 1665 kb)
Putative binding sites of HOXD-AS1 and SOX4 with miR-130a-3p. (A) Expression of potential miRNAs in HCC in TCGA cohorts. (B) Comparison summary of miR-130a-3p target sites in HOXD-AS1 and SOX4. The red nucleotides (target sites) were deleted in the mutant constructs. (C) pGL3 luciferase reporter constructs containing wild type and mutated putative binding sites of HOXD-AS1 or SOX4 transcripts were shown. (TIF 1598 kb)
Competitive binding activities of HOXD-AS1 and SOX4 to miR-130a-3p. (A) Binding activities of HOXD-AS1 and SOX4 to miR-130a-3p in SNU449 cells with HOXD-AS1 overexpression. (B) Binding activities of HOXD-AS1 and SOX4 to miR-130a-3p in Huh7 cells with HOXD-AS1knockdown. (TIF 973 kb)