- Open Access
Molecular interplay between linc01134 and YY1 dictates hepatocellular carcinoma progression
Journal of Experimental & Clinical Cancer Research volume 39, Article number: 61 (2020)
Revealing the mechanical role of long non-coding RNAs (lncRNAs) in tumorigenesis can contribute to novel therapeutic target for cancers. The regulatory role of linc01134 in hepatocellular carcinoma (HCC) has not been studied yet.
Materials and methods
qRT-PCR and western blot were conducted to measure relevant RNA and protein expressions. CCK-8, colony formation, EdU, flow cytometry, wound-healing, transwell assays and xenograft experiments were performed to determine the role of linc01134 in HCC. ChIP and luciferase reporter assays were performed to analyze the effects of Yin Yang-1 (YY1) on linc01134 transcription activity. Relevant mechanical experiments were performed to verify interaction between relative genes.
YY1 enhanced linc01134 transcription by interacting with linc01134 promoter. Knockdown of linc01134 inhibited proliferation, migration and epithelial-mesenchymal transition (EMT), yet promoting apoptosis in HCC cells. Mechanically, linc01134 acted as miR-324-5p sponge and interacted with insulin-like growth factor 2 mRNA binding protein 1 (IGF2BP1) to increase the stability of YY1 mRNA expression. Up-regulated YY1 continuously stimulated linc01134 expression by enhancing linc01134 promoter activity, forming a positive feedback loop.
Linc01134/miR-324-5p/IGF2BP1/YY1 feedback loop mediates HCC progression, which possibly provide prognosis and treatment target of HCC.
With approximately 850, 000 new diagnosis cases each year, liver cancer is the second leading cause of cancer-related death globally . Hepatocellular carcinoma (HCC) is the most prevalent subtype of liver cancer, accounting for about 90% amongst liver cancer cases. Although advance in therapeutic methods has been made in recent years, the mortality of HCC still ranks the third among all cancers. Worse still, this rate is still on the rise worldwide . Frequent recurrence and metastasis are major reasons for the high mortality of HCC. From molecular level, the unclear pathogenesis behind HCC is partially responsible for the clinical treatment stagnancy. Hence, in-depth study of the molecular mechanisms behind HCC tumorigenesis is of vital significance.
In human genome, the majority of transcripts are non-coding RNAs, while only approximately 1.2% of these transcripts represent protein-coding genes . The numerous non-protein-coding transcripts have been discovered to play important regulatory role in diseases, including cancers. Among these, long non-coding RNAs (lncRNAs), which are over 200 nucleotides in length, have emerged as key regulators in biological processes for their diverse and complex work mechanisms . LncRNAs have been reported to affect an array of cellular functions, such as cellular growth, apoptosis, migration and invasion, as well as EMT progression [5,6,7,8]. Also, increasing lncRNAs have been detected to be aberrantly expressed in cancers. Besides, the expression of lncRNAs could be impacted by transcription factors or pioneer factor. Transcription factor or pioneer factor induced lncRNA promoter region activation or silence to further alter the cellular function elicited by lncRNAs. For instance, SP1 transcriptionally activated lncRNA ZFAS1 to accelerate colorectal cancer progression via targeting downstream miR-150-5p/VEGFA axis . Yin Yang-1 (YY1) is a multifunctional transcription factor that can promote or suppress the promoter activity of various genes. Additionally, YY1 has been reported to be involved in the EMT progression and biological functions in colorectal cancer . YY1 was also erratically expressed in HCC and aggravate its progression. For instance, YY1 enrichment contributed to EZH2 recruitment for H3K27me3-regulated microRNAs silence, consequently activating NF-κB signaling in HCC carcinogenesis .
Recently, growing number of lncRNAs have been reported to function as miRNAs sponges and further mediate the target of miRNAs, constituting a competing endogenous RNA (ceRNA) network. This post-transcriptional molecular mechanism of lncRNAs has gained increasing attention [12, 13]. Cytoplasmic lncRNA can restore the expression of target gene via sponging miRNA, thus mediating the cellular function of target gene. Of note, lncRNAs are also elucidated to modulate gene expression via interacting with RNA binding proteins (RBPs). For instance, LINC01093 curbed HCC progression by interacting with IGF2BP1 and interfering the binding between IGF2BP1 and GLI1, thus impacting expression of GLI1 downstream molecules associated with HCC progression . Linc01134 is a novel identified long non-coding RNA without the potential to encode proteins. However, its molecular role in cancer has not been studied yet. Whether linc01134 regulates the progression of HCC via mediating the expression of target genes remain unclear. In this study, we investigated underlying mechanism of linc01134 in HCC tumorigenesis and progression.
Material and methods
This research was conducted under the approval and supervision of the Ethics Committee of Shandong Provincial Hospital Affiliated to Shandong University 46 recruited patients without any preoperative chemotherapy or radiation therapy had signed the written informed consent before study. After surgical resection, both HCC tissues and para-tumor tissues were instantly stored at − 80 °C in liquid nitrogen till RNA extraction.
Cell culture and treatment
Normal liver cell (QSG-7701) and HCC cells (SNU-449, HepG2, HLF, Huh-7) were purchased from the American Type Culture Collection (ATCC; Manassas, VA, USA) and cultured at 37 °C in a humidified atmosphere containing 5% CO2. DMEM (Gibco, Grand Island, NY, USA) adding antibiotics (Gibco) and 10% FBS (Gibco) were utilized for culturing cells. The RNA polymerase II inhibitor, α-amanitin (100 nM, Sigma-Aldrich, St. Louis, MO, USA) was added into culture medium for treating cells.
Specific shRNAs against linc01134 (sh-linc01134#1/2), YY1 (sh-YY1) or IGF2BP1 (sh-IGF2BP1) and the non-specific shRNAs as negative control (sh-NC), along with the pcDNA3.1/YY1, pcDNA3.1/linc01134 or pcDNA3.1/IGF2BP1 and the empty vectors as control, were constructed by Genechem (Shanghai, China). MiR-324-5p mimics and NC mimics were from GenePharma (Shanghai, China). SNU-449 or HepG2 cells were transfected individually with these plasmids via Lipofectamine 3000 (Invitrogen, Carlsbad, CA, USA). The transfected cells were collected and the stably transfected cells were selected out after 48 h of 2 μg/ml puromycin treatment.
Quantitative real-time polymerase chain reaction (qRT-PCR)
SNU-449 or HepG2 cells were re-suspended in 1 mL of TRIzol (Invitrogen), total cellular RNA was isolated, and reverse transcribed into cDNA. Gene expression was detected by qRT-PCR. GAPDH or U6 acted as endogenous control.
Cell counting Kit-8 (CCK-8) assay
1 × 103 transfected SNU-449 or HepG2 cells were added to 96-well plates before CCK-8 (Beyotime, Jiangsu, China) was added. Absorbance at 450 nm was read via microplate reader (bio-rad, Hercules, CA, USA)
Colony formation assay
Transfected SNU-449 or HepG2 cells were planted in 24-well plates (500 cells/well). Colonies were fixed for 15 min by 100% methanol (Sigma-Aldrich) following 14-day culturing. Cells were stained for 20 min utilizing 0.5% crystal violet (Sigma-Aldrich). Colonies with > 50 cells were counted manually.
EdU staining was undertaken with a Click iT™ EdU cell proliferation assay kit (Invitrogen). Cell nuclei were double-dyed for 10 min using DAPI (Sigma-Aldrich) and observed via florescence microscope (Olympus, Tokyo, Japan).
Annexin V-FITC/PI Apoptosis Detection Kit (Sigma-Aldrich) was applied to stain transfected SNU-449 or HepG2 cells in the dark. The apoptotic rate was assayed with flow cytometer (BD Biosciences, Franklin Lakes, NJ, USA) and CellQuest software (BD Biosciences).
Animal study was undertaken with the approval of the Animal Research Ethics Committee of Shandong Provincial Hospital Affiliated to Shandong University BALB/c nude mice (6 weeks old) were bought from Shi Laike Company (Shanghai, China). SNU-449 cells transfected with sh-linc01134 or sh-NC were injected subcutaneously to animals. Tumors volume was assessed every 4 days. Following 4 weeks, mice were killed. Tumors were harvested and weighed.
Xenograft tumor tissue samples were immunostained for Ki-67 with the specific antibody from Santa Cruz Biotechnology (Dallas, TX, USA) as per the established protocol.
Transfected cells in serum-free medium were planted into the Transwell upper chambers (Corning, Cambridge, MA, USA). 10% FBS-medium was added to the bottom chambers. In invasion assay, upper chambers were pre-coated with Matrigel (BD Biosciences). 48 h later, migratory cells were fixed, stained and counted under microscope (Olympus).
Transfected SNU-449 or HepG2 cells were cultured in 6-well plates to 90% confluence. Wound was made utilizing pipette tip. The widths of the wounds were examined at 0 and 24 h.
Transfected SNU-449 or HepG2 cells were incubated in sequence with the specific primary antibodies to E-cadherin (Abcam) and N-cadherin (Abcam), as well as secondary antibodies. Hoechst 33342 (Thermo Fisher Scientific, Waltham, MA, USA) was used for counterstaining nuclei.
Total cellular protein samples were separated with 10% SDS PAGE and transferred electrophoretically onto PVDF membranes (Millipore, Billerica, MA, USA). Primary antibodies against E-cadherin (ab40772), N-cadherin (ab76057), Vimentin (ab8978, Abcam), IGF2BP1 (ab82968, Abcam), YY1 (ab109228, Abcam), CDK4 (ab199728, Abcam), cyclin D1 (ab16663, Abcam), CDK2 (ab32147, Abcam) and GAPDH (ab9485, Abcam), along with HRP-conjugated secondary antibody were obtained and employed. GAPDH served as internal control.
Haematoxylin and eosin (HE) staining
Xenograft tumor tissue samples were fixed and embedded by paraffin. Sample sections with 4 μm of thickness were utilized for HE staining (Sigma-Aldrich).
Chromatin Immunoprecipitation (ChIP)
Chromatin was crosslinked and sonicated to fragments (500 bp). Immunoprecipitation was implemented by magnetic beads coated with anti-YY1 or anti-IgG (Millipore). Precipitated DNA fragments were studied with qRT-PCR following total RNA extraction.
Luciferase reporter assay
The pGL3-linc01134 promoter vector was constructed via cloning the wild-type or mutant interacting sequences of YY1 in linc01134 promoter into pGL3-Basis vector (Promega, Madison, WI, USA), co-transfected into with sh-NC and sh-YY1 or pcDNA 3.1 and pcDNA3.1/YY1. The wild-type or mutant sequences of miR-324-5p in linc01134 or IGF2BP1 were sub-cloned into pmirGLO dual-luciferase vector (Promega), co-transfected with miR-324-5p mimics or NC mimics into HEK-293 T, SNU-449 or HepG2 cells. Fourty-eight hour later, dual luciferase reporter assay system (Promega) was applied for examination of luciferase activities.
Cytoplasmic and Nuclear RNA Purification Kit was acquired commercially from Norgen (Thorold, ON, Canada). Expressions of GAPDH, U6 and linc01134 in nuclear and cytoplasm fractions of SNU-449 or HepG2 cells were examined via qRT-PCR.
Fluorescence in situ hybridization (FISH)
Cy3-labeled Linc01134 probe was synthesized by RiboBio. Fluorescent in Situ Hybridization Kit (RiboBio) was used to assess fluorescence signal. Cell nuclei were counterstained using Hoechst.
Cell lysates were individually incubated with biotinylated RNAs, followed by adding magnetic beads. The RNA enrichment was evaluated by qRT-PCR or western blotting.
RNA Immunoprecipitation (RIP)
RIP was undertaken in SNU-449 or HepG2 cells using Magna RIP™ RNA-Binding Protein Immunoprecipitation Kit (Millipore) and antibody against normal mouse IgG as control. Finally, qRT-PCR analysis determined the RNA enrichment in the immunoprecipitated complex.
Cell cycle analysis
A single-cell suspension was prepared for SNU449 and HepG2 cells in ice-cold phosphate-buffered saline (PBS). Next, the cells were fixed with 70% ethanol for 12 h at 4 °C. Cell clumps were removed by passing-through a cell strainer (BD BioSciences, CA, USA). Four-hundred microliters of propidium iodide (0.1 mg/mL) were incubated for 10 min and detected by flow cytometry and DNA content data was detected using BD Accuri Cytometry (BD BioSciences, Franklin Lakes, NJ, USA).
Data of at leastthree independent replications were analyzed by SPSS Vision 19.0 (SPSS, Chicago, IL, USA) and showed as means ± SD. Pearson’s correlation coefficient was adopted for evaluating the gene expression correlations. Differences in groups were evaluated by Student’s t-test or ANOVA, with p-value less than 0.05 as statistical significance.
Knockdown of linc01134 can dampen HCC cell proliferation in vitro and suppress tumor growth in vivo
To begin with, we analyzed the expression status of linc01134 in liver hepatocellular carcinoma (LIHC) tissues and para-tumor tissues samples from starBase. It can be noticed that linc01134 was remarkably overexpressed in LIHC tissues compared with normal tissue (Fig. 1a). Several databases manifested the significantly low expression of linc01134 in normal liver tissues (Fig. S1A, S1B, S1C). Besides, linc01134 was presented by LncBook (https://bigd.big.ac.cn/lncbook/index) to have no protein-coding potential (Fig. S1D). Kaplan-Meier survival analysis found that high expression of linc01134 was closely associated with shorter overall survival, while low expression of linc01134 had the opposite outcome (Fig. 1b). We then analyzed the expression of linc01134 in 46 paired collected HCC tissues and normal tissue. qRT-PCR result manifested that linc01134 was overtly up-regulated in HCC tissues (Fig. 1c). Likewise, qRT-PCR data also demonstrated an aberrant up-regulation of linc01134 in HCC cell lines (SNU-449, HepG2, HLF and Huh-7) compared with normal liver epithelial cell line (QSG-7701) (Fig. 1d). Based on the above data, we assumed that linc01134 might play a role in HCC. To verify this assumption, we knocked down linc01134 in SNU-449 and HepG2 cells to study its influence on cellular behavior. By transfecting sh-linc01134#1/2, we got low expression of linc01134 in SNU-449 and HepG2 (Fig. 1e). According to CCK-8 result, knockdown of linc01134 remarkably reduced cell viability compared with control group (Fig. 1f). Meanwhile, silencing linc01134 could suppress cell proliferation by counting colony cells and calculating EdU positive cell rate (Fig. 1g, h). Besides, PI-FACS analysis showed that linc01134 knockdown increased the percentage of cells at the G0/G1 phase, while decreased that in S and G2/M phases (Fig. S2A). We measured the mRNA and protein levels of CDK4, cyclin D1, CDK2, which were associated with cell cycle. Data presented that all these three genes were down-regulated by linc01134 depletion, suggesting that sh-linc01134 could induce cell cycle arrest (Fig. S2B). We performed flow cytometry to examine the effect of sh-linc01134 on cell apoptosis. An increase in the fraction of apoptotic cells was observed in cells transfected with sh-linc01134 (Fig. 1i). We used xenografts to study the effects of linc01134 on tumor growth in vivo. The tumor growth speed slowed down in sh-linc01134 group in comparison with control group (Fig. 1j). We found that the tumors in sh-linc01134 group were obviously smaller and lighter compared with those in control group (Fig. 1k, l). Consistently, immunohistochemistry (IHC) result revealed that Ki-67 positive cells were conspicuously decreased in sh-linc01134 group than control group (Fig. 1m). In short, linc01134 was up-regulated in HCC tissues and cell lines. Linc01134 down-regulation hampered HCC cell proliferation in vitro and repressed tumor growth in vivo.
Knockdown of linc01134 can curb HCC cell migration, invasion and EMT process in vitro and impede metastasis in vivo
Since migration and metastasis are crucial for HCC deterioration, we evaluated the effects of sh-linc01134 on HCC cell migration and invasion abilities. Transwell assay was conducted to evaluate cell invasion. The result found that knockdown of linc01134 greatly decreased the invaded cells (Fig. 2a). Transwell migration and wound-healing assays demonstrated that silencing linc01134 inhibit cell migration ability to a large extent (Fig. 2b, c).
EMT is a crucial phase that facilitates tumor cell migration and metastasis in HCC. In IF assay, down-regulating linc01134 evidently augmented the expression level of E-cadherin, whereas reduced the level of N-cadherin (Fig. 2d). Based on this phenomenon, we supposed that sh-linc01134 could exert inhibitory effect on EMT process. To further confirm this notion, western blot was performed to measure EMT-related proteins. We identified an increase of E-cadherin, yet a reduction of N-cadherin and Vimentin protein levels (Fig. 2e). Experiments in vivo were conducted to further validate the effects of linc01134 on migration and metastasis. HE staining result demonstrated the decreased metastasis nodes in sh-linc01134 group than control group (Fig. 2f). Altogether, linc01134 knockdown can curb HCC cell migration, invasion and EMT process in vitro and impede metastasis in vivo.
YY1 stimulates linc01134 transcription
Transcription factor has been found to induce or repress the transcription activity of lncRNA . Given the high expression of linc01134 in HCC, we presumed that transcription factor might play a role. We found potential transcription factor candidate that might interact with linc01134 promoter region using UCSC and Promo database, and noticed that YY1 was the only common transcription factor predicted by UCSC, Promo and JASPAR (Fig. 3a). We then detected the expression of YY1 in HCC tissues and found that YY1 was significantly overexpressed in HCC tissues (Fig. S3A). Correlation analysis revealed a positive relevance between YY1 and linc01134 in HCC tissues (Fig. S3B). To investigate the impacts of YY1 on linc01134 expression, firstly, we overexpressed the expression of YY1 (Fig. 3b). Strikingly, linc01134 expression was increased significantly due to YY1 overexpression (Fig. 3c). Similarly, low YY1 expression was obtained by transfecting sh-YY1 (Fig. 3d) and data showed a significant decline in linc01134 expression after knockdown of YY1 (Fig. 3e). This observation indicated that YY1 indeed regulated linc01134 expression. Secondly, ChIP assay results disclosed that a striking enrichment of linc01134 promoter in anti-YY1 group than anti-IgG control group (Fig. 3f). Moreover, knockdown of YY1 suppressed the promoter luciferase activity, while YY1 overexpression enhanced that of linc01134 promoter, evidenced by luciferase reporter assay (Fig. 3g). This phenomenon indicated that YY1 played a positively regulatory role in linc01134 transcription. Subsequently, we aimed to find out the precise binding site for YY1 at linc01134 promoter region since that two possible binding sites (− 773 to − 762, − 1127 to-1116) were predicted by JASPAR (Fig. 3h). The ChIP data manifested that YY1 could bind to linc01134 P3 promoter region (Fig. 3i). Next, we determined to verify whether − 773 to − 762 sequence existing in P3 that bind to YY1. We constructed Wt-P3 and Mut-P3 for the purpose of luciferase reporter assay. After co-transfecting pcDNA3.1/YY1 with Wt-P3 or Mut-P3, we observed that the promoter activity was greatly enhanced in Wt-P3 group, while no significant change in Mut-P3 group (Fig. 3j). Collectively, we found that YY1 stimulated the transcriptional activity of linc01134 by interacting with linc01134 promoter at − 773 to − 762 site upstream TSS.
Linc01134 functions as ceRNA to up-regulate IGF2BP1 via sponging miR-324-5p
To explore how linc01134 induce cell phenotype changes in HCC, we examined the underlying function mechanism. The ceRNA hypothesis has revealed an important interaction mechanism between RNAs. We testified whether linc01134 functions as a competing endogenous gene in HCC. We searched from starBase (http://starbase.sysu.edu.cn/) and DIANA (http://carolina.imis.athena-innovation.gr/diana_tools/web/index.php?r=site%2Ftools), and detected two shared miRNAs (miR-324-5p and miR-6512-3p) of linc01134 from two databases (Fig. 4a). By knocking down the expression of linc01134, miR-324-5p expression increased sharply but miR-6512-3p expression exhibited no significant change (Fig. 4b). More importantly, a significant down-regulation of miR-324-5p was discovered in HCC tissues, while no distinct difference between the expression of miR-6512-3p in HCC tissues and that of normal one (Fig. S3C). A negative correlation between miR-324-5p and linc01134 in HCC tissues was revealed by Pearson correlation analysis (Fig. S3D), indicating the sponge role of linc01134 possible. Therefore, we selected miR-324-5p as a candidate miRNA to study. We identified a putative binding site of miR-324-5p with linc01134 by browsing starBase (Fig. 4c). Next, subcellular fractionation and FISH assays were conducted to determine the subcellular location of linc01134. The findings showed that linc01134 was principally situated in cytoplasm (Fig. 4d). Luciferase reporter verified the combination between linc01134 and miR-324-5p. It indicated that miR-324-5p mimics significantly dampened the luciferase activity of linc01134-WT, but not linc01134-MUT (Fig. 4e). RNA pull down assay further proved this physical interaction. We observed enriched miR-324-5p was pulled down by biotinylated linc01134-WT compared with linc01134-MUT group (Fig. 4f). Further, we screened the target gene of miR-324-5p from DIANA and found that only five target genes prediction score exceeds 0.8. We transfected miR-324-5p mimics in cells and found that IGF2BP1 expression lessened most dramatically compared with balance mRNA candidates (Fig. 4g). qRT-PCR revealed that IGF2BP1 exhibited significantly high expression in HCC tissues, which is the same status with that of linc01134 (Fig. 4h). In addition, Pearson correlation analysis manifested positive association between IGF2BP1 and linc01134 (Fig. 4i), while a negative one between IGF2BP1 and miR-324-5p (Fig. S3E). Besides, we also discovered that IGF2BP1 expression in HCC tissues was proportional to the level of YY1, the upstream regulator of linc01134, as well (Fig. S3F). The potential binding site of miR-324-5p with IGF2BP1 was found by starBase (Fig. 4j). Luciferase reporter assay manifested that miR-324-5p up-regulation could impair the luciferase activity of IGF2BP1-WT (Fig. 4k). RIP assay indicated the co-existence of linc01134, miR-324-5p and IGF2BP1 in RNA-induced silencing complexes (RISCs) (Fig. 4l). We then examined the overexpression efficiency of linc01134 (Fig. 4m). Interestingly, RIP-based qRT-PCR found that the enrichment of IGF2BP1 in anti-Ago2 group was significantly decreased when linc01134 was overexpressed (Fig. 4n), indicating that linc01134 competitively bound with miR-324-5p against IGF2BP1 in the same RISCs. Furthermore, IGF2BP1 mRNA expression reduced sharply after miR-324-5p was overexpressed, but the falling tendency was reversed by linc01134 up-regulation, which was consistent with changes of IGF2BP1 protein expression by western blot (Fig. 4o). We further found that both IGF2BP1 and linc01134 were enriched in biotinylated miR-324-5p-WT group (Fig. S4A, S4B). Notably, overexpressing linc01134 markedly diminished the abundance of IGF2BP1 in biotinylated miR-324-5p-WT group (Fig. S4C). More importantly, this phenomenon seemed to be resulted by the enhanced enrichment of linc01134 precipitated in anti-Ago2 group when overexpressing linc01134, while Ago2 was the core component of miR-324-5p-guided RISCs (Fig. S4D). Above findings uncovered that cytoplasmic linc01134 functioned as ceRNA to release IGF2BP1 suppressed by miR-324-5p.
IGF2BP1 enhances the stability of stabilizes of YY1 mRNA in HCC
YY1 has been reported to be overexpressed in multiple cancers, including HCC . Functionally, YY1 could silence tumor-suppressive miRNAs in HCC . Data from starBase found that the expression of IGF2BP1 was positively related to that of YY1 in HCC tissues (Fig. 5a), so was data in 46 HCC tissues obtained clinically. IGF2BP1, a kind of RNA-binding protein (RBP), has been revealed to prevent the degradation of c-myc RNA . To study whether IGF2BP1 affected YY1 expression via preventing its mRNA degradation process, we firstly overexpressed and knocked down IGF2BP1 to observe the change of YY1 mRNA expression. Data displayed that YY1 expression was elevated when IGF2BP1 was up-regulated while was decreased with the down-regulation of IGF2BP1 (Fig. 5b-c). These results verified the positive regulator role of IGF2BP1 on YY1. Thereafter, we examined that whether IGF2BP1 could also modulate YY1 mRNA stability. We knocked down IGF2BP1 in SNU-449 and HepG2 cells, then used qRT-PCR to evaluate the expression changes of YY1 and GAPDH over 24 h period after 48-h α-amanitin (50 μM) treatment. Results demonstrated that knockdown of linc01134 shortened the half-life of YY1 mRNA compared with control group (Fig. 5d). To sum up, IGF2BP1 can enhance the stability of YY1 mRNA and up-regulate its mRNA expression consequently.
Linc01134 binds to IGF2BP1 protein and strengthen the physical interaction between YY1 mRNA and IGF2BP1 protein
Lately, it has been reported that lncRNAs could affect mRNA stability to regulate gene expression through physical interaction with IGF2BP1 protein . Considering that IGF2BP1 is an RBP, we preliminarily speculated that linc01134 may also bind to IGF2BP1 and impact the association between IGF2BP1 and YY1 mRNA to exert its function in HCC. To verify this hypothesis, we conducted RIP with anti-IGF2BP1 initially. Evident enrichment of linc01134 and YY1 mRNA was observed in anti-IGF2BP1 group compared with anti-IgG group (Fig. 5e). In addition, pull down results revealed that IGF2BP1 protein expression significantly was enriched in linc00134 or YY1 group (sense) compared with negative control group (antisense) (Fig. 5f). Above results manifested the physical interaction between IGF2BP1 and linc01134, as well as between IGF2BP1 and YY1 mRNA.
To investigate whether linc00134 impacted the interaction between IGF2BP1 and YY1 mRNA by binding to IGF2BP1, we silenced linc01134 and conducted RIP assay with anti-IGF2BP1. YY1 mRNA co-precipitated by anti-IGF2BP1 was decreased in linc00134 knockdown group compared with control (Fig. 5g), suggesting that linc01134 enhanced the physical interaction between IGF2BP1 protein and YY1 mRNA.
To further testify whether linc01134 regulated YY1 mRNA stability and mRNA expression via recruiting IGF2BP1, we first performed qRT-PCR after blocking new RNA synthesis with α-amanitin (50 μM) among sh-NC, sh-linc01134 and sh-linc01134 + pcDNA3.1/IGF2BP1 groups. The half-life of YY1 mRNA was significantly shortened by sh-linc01134, but was extended to some extent after co-transfection of pcDNA3.1/IGF2BP1 compared with control group (Fig. 5h). Next, we knocked down IGF2BP1 in pcDNA3.1 or pcDNA3.1/linc01134 transfected cells and found that YY1 mRNA expression was barely changed (Fig. 5i). In a summary, the increase of YY1 mRNA stability regulated by linc01134 was dependent on recruitment of IGF2BP1. Linc01134 interacted with IGF2BP1 protein and strengthen the physical interaction between YY1 mRNA and IGF2BP1 protein.
YY1 rescues the oncogenic function of linc01134 in HCC
A series of rescue experiments were conducted to investigate whether YY1 affect the modulating function of linc01134 in HCC. Results of qRT-PCR and western blot showed that linc01134 knockdown could down-regulate the expression of YY1 mRNA and protein levels, and then overexpression of YY1 reversed this repressing influence (Fig. 6a). CCK-8 result indicated that up-regulation of YY1 could offset the inhibitory effects of sh-linc01134 on cell viability (Fig. 6b). Additionally, colony formation and EdU assays further proved that YY1 enrichment markedly counteracted the obstructive effects of sh-linc01134 on cell proliferation (Fig. 6c, d). Also, cell cycle arrest induced by sh-linc01134 was remedied by overexpressing YY1 (Fig. S5A, S5B). Flow cytometry results manifested that up-regulation of YY1 greatly neutralized the promoting impacts of sh-linc01134 on cell apoptosis (Fig. 6e). On the contrary, YY1 overexpression could greatly countervail the suppressing effects of linc01134 knockdown on cell migration and invasion abilities (Fig. 6f, g, h). YY1 up-regulation restored inhibitory impacts imposed by linc01134 down-regulation on EMT progression, as evidenced by the changes of E-cadherin and N-cadherin levels (Fig. 6i). Besides, western blot results further validated that overexpression of YY1 rescued the impeditive function of linc01134 knockdown on EMT progression (Fig. 6j). These results elucidated that linc01134 affected the proliferative and migrating phenotype of HCC cells via engaging in the miR-324-5p/IGF2BP1/YY1 axis.
In a summary, YY1 protein can bind to the promoter region of linc01134 and activate the transcription of linc01134 in nucleus. Linc01134 finalizes transcription and transfers into cytoplasm. Cytoplasmic linc01134 competitively binds to miR-324-5p and up-regulates the expression of IGF2BP1. Further, linc01134 associates with IGF2BP1 protein and strengthens the interaction between IGF2BP1 protein and YY1 mRNA, consequently increasing the stability and expression of YY1 mRNA. YY1 mRNA up-regulation induces its protein enrichment. Up-regulated YY1 continuously stimulated linc01134 expression by binding to linc01134 promoter, forming a positive feedback loop in HCC cells (Fig. 7). Linc01134/miR-324-5p/IGF2BP1/YY1 feedback loop promotes cell proliferation, migration and EMT process in HCC. Therefore, targeting linc01134/miR-324-5p/IGF2BP1/YY1 might contribute to providing a new treatment strategy for HCC.
In present study, we examined the expression of the novel long non-coding linc01134 in HCC tissues and para-tumor tissues as well as tumor cells lines and normal cell line. The results manifested that linc01134 was up-regulated in HCC tissues and cell lines, which was attributed to the transcription activity stimulation by YY1. On the one hand, knockdown of linc01134 sharply dampened HCC cell proliferation, yet promoting apoptosis. On the other hand, linc01134 silence deterred cell migration, invasion and EMT progression in vitro. Experiments in vivo further proved that silencing linc01134 could slow down implanted tumor growth and metastasis. Therefore, our findings revealed that linc01134 aggravated the malignant cell phenotype in HCC.
YY1 can function as transcription factor and plays an important regulatory role in various biological processes by binding to DNA and numerous proteins, consequently impacting the pathologic processes of cancers [19,20,21,22]. YY1 has been reported to induce lncRNA MCM3AP-AS1 overexpression in lung cancer, and overexpressed MCM3AP-AS1 promoted angiogenesis via targeting KPNA4 through shared miR-340-5p . Moreover, YY1-activated LINC00673 has been revealed to promote proliferative activity of breast cancer cells via regulating MARK4 through sponging miR-515-5p . In current study, we investigated whether YY1 induced linc01134 transcription and affected its molecular regulatory mechanism.
IGF2BP1 belongs to RBPs that can bind to transcripts and regulate RNA stability. Predominantly expressed in the embryo, IGF2BP1 is an essential mediator through the interaction with RNAs in various cancers [25,26,27]. IGF2BP1 exerts synthesized effects, including cell proliferation, adhesion [28, 29]. In this study, we found that IFG2BP1 was the target of miR-324-5p and its expression was negatively correlated with miR-324-5p expression in HCC tissues. Additionally, linc01134 competitively bound to miR-324-5p against IFG2BP1. The expression of IFG2BP1 mRNA and protein levels were significantly down-regulated after overexpressing miR-324-5p but increased after up-regulating linc01134. Cytoplasmic linc01134 acted as miR-324-5p sponge, consequently elevating the expression of IFG2BP1 mRNA and protein levels. Previous study validated that IGF2BP1 was a crucial protein in modulating EMT progression. IGF2BP1 promoted mesenchymal-like cell properties in tumor-derived cells, and was experimentally proved to serve as a pro-mesenchymal factor . Therefore, it was reasonable to deduce that the linc01134 knockdown repressed EMT progression via diminishing IGF2BP1 expression released from the miR-324-5p induced RISCs. And these evidences suggested that IGF2BP1 was responsible for the acceleration function of linc01134 on EMT progression in HCC.
Apart from the classical ceRNA gene modulation network, increasing lncRNAs have been validated to cooperate with RBPs to modulate mRNA protein stability. For instance, lncRNA THOR was found to up-regulate c-myc expression via strengthening its combination with TGF2BP1 . Interestingly, we found that linc01134 bond to IGF2BP1 and further enhanced physical interaction between IGF2BP1 and YY1 mRNA, finally enhancing the effects of IGF2BP1 on promoting YY1 mRNA stability. Linc01134 depended on IGF2BP1 in stabilizing YY1 mRNA stability. The enhanced YY1 mRNA stability led to the up-regulation of YY1 protein, thus continuously stimulating the transcription of linc00134, forming a positive feedback loop. Besides, the facilitating role of YY1 in the tumorigenesis and EMT of HCC has also been indicated previously [32, 33], which may also explains the contribution of linc01134 to HCC progression. However, due to time limitation, we only explored the YY1 mRNA stability and expression profile changes after the interaction between linc01134 and IGF2BP1. IGF2BP1 can bind to a wide array of other RBPs and forms a distinct ribonucleoprotein complex . Whether linc01134 also impacted other biological functions of IGF2BP1 via physical interaction with IGF2BP1 requires further investigation.
Taken together, our findings indicate that linc00134 is a meaningful prognostic factor for HCC. Mechanically, YY1-induced linc01134 transcription can promote YY1 mRNA stability and expression via releasing miR-324-5p-targeted IGF2BP1 mRNA and recruiting IGF2BP1, which forms a positive feedback loop in HCC cells. YY1-induced linc01134 up-regulation can promote HCC cell proliferation, migration and invasion, while suppress apoptosis, suggesting that linc01134 might be an essential molecular marker for prognosis and a target for HCC treatment.
Availability of data and materials
Analysis of variance
American type culture collection
Cell counting kit-8
Competing endogenous RNA
Dulbecco’s modified eagle medium
Enhancer of zeste homolog 2
Fetal bovine serum
Fluorescence in situ hybridization
Glioma-associated oncogene homolog 1
Histone3 lysine27 tri-methylation
- HE Staining:
Haematoxylin and eosin
Insulin-like growth factor 2 mRNA binding protein 1
Liver hepatocellular carcinoma
Long non-coding RNAs
Microtubule affinity regulating kinase 4
Minichromosome maintenance complex component 3 associated protein antisense RNA 1
Nuclear factor kappa B
Quantitative real-time polymerase chain reaction
RNA binding proteins
- SDS PAGE:
Sodium dodecyl sulfate polyacrylamide gel electrophoresis
Transcription termination site
Vascular endothelial growth factor
Zinc finger antisense 1
Hepatocellular carcinoma. Nat Rev Dis Primers. 2016;2:16019. https://doi.org/10.1038/nrdp.2016.19.
Bertuccio P, Turati F, Carioli G, Rodriguez T, La Vecchia C, Malvezzi M, et al. Global trends and predictions in hepatocellular carcinoma mortality. J Hepatol. 2017;67(2):302–9.
ENCODE Project Consortium. An integrated encyclopedia of DNA elements in the human genome. Nature. 2012;489(7414):57–74.
Wilusz JE, Sunwoo H, Spector DL. Long noncoding RNAs: functional surprises from the RNA world. Genes Dev. 2009;23(13):1494–504.
Qi C, Xiaofeng C, Dongen L, Liang Y, Liping X, Yue H, et al. Long non-coding RNA MACC1-AS1 promoted pancreatic carcinoma progression through activation of PAX8/NOTCH1 signaling pathway. J Exp Clin Cancer Res. 2019;38(1):344.
Wang N, Zhang C, Wang W, Liu J, Yu Y, Li Y, et al. Long noncoding RNA DANCR regulates proliferation and migration by epigenetically silencing FBP1 in tumorigenesis of cholangiocarcinoma. Cell Death Dis. 2019;10(8):585.
Shen SN, Li K, Liu Y, Yang CL, He CY, Wang HR. Down-regulation of long non-coding RNA PVT1 inhibits esophageal carcinoma cell migration and invasion and promotes cell apoptosis via microRNA-145-mediated inhibition of FSCN1. Mol Oncol. 2019;13(12):2554–73.
Jia Y, Duan Y, Liu T, Wang X, Lv W, Wang M, et al. LncRNA TTN-AS1 promotes migration, invasion, and epithelial mesenchymal transition of lung adenocarcinoma via sponging miR-142-5p to regulate CDK5. Cell Death Dis. 2019;10(8):573.
Chen X, Zeng K, Xu M, Hu X, Liu X, Xu T, et al. SP1-induced lncRNA-ZFAS1 contributes to colorectal cancer progression via the miR-150-5p/VEGFA axis. Cell Death Dis. 2018;9(10):982.
Ye Y, Gu B, Wang Y, Shen S, Huang W. YY1-induced Upregulation of long noncoding RNA ARAP1-AS1 promotes cell migration and invasion in colorectal Cancer through the Wnt/beta-catenin signaling pathway. Cancer Biother Radiopharm. 2019;(8):519–528.
Tsang DP, Wu WK, Kang W, Lee YY, Wu F, Yu Z, et al. Yin Yang 1-mediated epigenetic silencing of tumour-suppressive microRNAs activates nuclear factor-kappaB in hepatocellular carcinoma. J Pathol. 2016;238(5):651–64.
Salmena L, Poliseno L, Tay Y, Kats L, Pandolfi PP. A ceRNA hypothesis: the Rosetta stone of a hidden RNA language? Cell. 2011;146(3):353–8.
Cesana M, Cacchiarelli D, Legnini I, Santini T, Sthandier O, Chinappi M, et al. A long noncoding RNA controls muscle differentiation by functioning as a competing endogenous RNA. Cell. 2011;147(2):358–69.
He J, Zuo Q, Hu B, Jin H, Wang C, Cheng Z, et al. A novel, liver-specific long noncoding RNA LINC01093 suppresses HCC progression by interaction with IGF2BP1 to facilitate decay of GLI1 mRNA. Cancer Lett. 2019;450:98–109.
Chen L, Gong X, Huang M. YY1-activated long noncoding RNA SNHG5 promotes Glioblastoma cell proliferation through p38/MAPK signaling pathway. Cancer Biother Radiopharm. 2019;34(9):589–96.
Kim JS, Son SH, Kim MY, Choi D, Jang IS, Paik SS, et al. Diagnostic and prognostic relevance of CP2c and YY1 expression in hepatocellular carcinoma. Oncotarget. 2017;8(15):24389–400.
Weidensdorfer D, Stohr N, Baude A, Lederer M, Kohn M, Schierhorn A, et al. Control of c-myc mRNA stability by IGF2BP1-associated cytoplasmic RNPs. RNA (New York, NY). 2009;15(\):104–115.
Yang F, Xue X, Zheng L, Bi J, Zhou Y, Zhi K, et al. Long non-coding RNA GHET1 promotes gastric carcinoma cell proliferation by increasing c-Myc mRNA stability. FEBS J. 2014;281(3):802–13.
Wu S, Wang H, Li Y, Xie Y, Huang C, Zhao H, et al. Transcription factor YY1 promotes cell proliferation by directly activating the pentose phosphate pathway. Cancer Res. 2018;78(16):4549–62.
Li M, Liu Y, Wei Y, Wu C, Meng H, Niu W, et al. Zinc-finger protein YY1 suppresses tumor growth of human nasopharyngeal carcinoma by inactivating c-Myc-mediated microRNA-141 transcription. J Biol Chem. 2019;294(15):6172–87.
Pentland I, Campos-Leon K, Cotic M, Davies KJ, Wood CD, Groves IJ, et al. Disruption of CTCF-YY1-dependent looping of the human papillomavirus genome activates differentiation-induced viral oncogene transcription. PLoS Biol. 2018;16(10):e2005752.
Wang Y, Wu S, Huang C, Li Y, Zhao H, Kasim V. Yin Yang 1 promotes the Warburg effect and tumorigenesis via glucose transporter GLUT3. Cancer Sci. 2018;109(8):2423–34.
Li X, Yu M, Yang C. YY1-mediated overexpression of long noncoding RNA MCM3AP-AS1 accelerates angiogenesis and progression in lung cancer by targeting miR-340-5p/KPNA4 axis. J Cell Biochem. 2019;(3):2258–226.
Qiao K, Ning S, Wan L, Wu H, Wang Q, Zhang X, et al. LINC00673 is activated by YY1 and promotes the proliferation of breast cancer cells via the miR-515-5p/MARK4/hippo signaling pathway. J Exp Clin Cancer Res. 2019;38(1):418.
Yang H, Fu G, Liu F, Hu C, Lin J, Tan Z, et al. LncRNA THOR promotes tongue squamous cell carcinomas by stabilizing IGF2BP1 downstream targets. Biochimie. 2019;165:9–18.
Wang L, Aireti A, Aihaiti A, Li K. Expression of microRNA-150 and its target gene IGF2BP1 in human osteosarcoma and their clinical implications. Pathol Oncol Res. 2019;25(2):527–33.
Xu Y, Zheng Y, Liu H, Li T. Modulation of IGF2BP1 by long non-coding RNA HCG11 suppresses apoptosis of hepatocellular carcinoma cells via MAPK signaling transduction. Int J Oncol. 2017;51(3):791–800.
Manieri NA, Drylewicz MR, Miyoshi H, Stappenbeck TS. Igf2bp1 is required for full induction of Ptgs2 mRNA in colonic mesenchymal stem cells in mice. Gastroenterology. 2012;143(1):110–21 e10.
Stohr N, Kohn M, Lederer M, Glass M, Reinke C, Singer RH, et al. IGF2BP1 promotes cell migration by regulating MK5 and PTEN signaling. Genes Dev. 2012;26(2):176–89.
Zirkel A, Lederer M, Stohr N, Pazaitis N, Huttelmaier S. IGF2BP1 promotes mesenchymal cell properties and migration of tumor-derived cells by enhancing the expression of LEF1 and SNAI2 (SLUG). Nucleic Acids Res. 2013;41(13):6618–36.
Shang Y. LncRNA THOR acts as a retinoblastoma promoter through enhancing the combination of c-myc mRNA and IGF2BP1 protein. Biomed Pharmacother. 2018;106:1243–9.
Han J, Meng J, Chen S, Wang X, Yin S, Zhang Q, et al. YY1 complex promotes quaking expression via super-enhancer binding during EMT of hepatocellular carcinoma. Cancer Res. 2019;79(7):1451–64.
Li Y, Kasim V, Yan X, Li L, Meliala ITS, Huang C, et al. Yin Yang 1 facilitates hepatocellular carcinoma cell lipid metabolism and tumor progression by inhibiting PGC-1β-induced fatty acid oxidation. Theranostics. 2019;9(25):7599–615.
The help of all laboratory members in this study is appreciated.
This study was funded by Shandong Provincial Science Foundation (Grant No. ZR2013HQ026).
Ethics approval and consent to participate
This research was conducted under the approval and supervision of the Ethics Committee of Shandong Provincial Hospital Affiliated to Shandong University. Animal study was undertaken with the approval of the Animal Research Ethics Committee of Shandong Provincial Hospital Affiliated to Shandong University.
Consent for publication
All authors agreed upon the publication of this manuscript.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 3: Supplementary Figure 3. A. Relative expression of YY1 in HCC and para-tumor tissues was measured by qRT-PCR. B. Pearson correlation analysis between the expression of YY1 and linc01134 in HCC tissues. C. Relative expression of miR-324-5p and miR-6512-3p in HCC and para-tumor tissues was detected by qRT-PCR. D-F. Pearson correlation analysis was used to study the expression correlation between miR-324-5p and linc01134, between IGF2BP1 and miR-324-5p as well as between IGF2BP1 and YY1. **P < 0.01. n.s. meant no significance.
Additional file 4: Supplementary Figure 4. A-B. RNA pull down assays measured the enrichment of IGF2BP1/linc01134 in Bio-miR-324-5p-WT/MUT group. C. RNA pull down assays the enrichment of IGF2BP1 in Bio-miR-324-5p-WT/MUT group when overexpressing linc01134. D. RIP assay detected the enrichment of linc01134 in anti-Ago2 group. *P < 0.05, **P < 0.01.
Additional file 5: Supplementary Figure 5. A. The cell cycle was detected after co-transfecting pcDNA3.1/YY1 into sh-linc01134 transfected HCC cells. B. The mRNA and protein changes of CDK4, cyclin D1 and CDK2 in differently transfected groups were detected by qRT-PCR and western blot assays. *P < 0.05, **P < 0.01.
About this article
Cite this article
Rong, Z., Wang, Z., Wang, X. et al. Molecular interplay between linc01134 and YY1 dictates hepatocellular carcinoma progression. J Exp Clin Cancer Res 39, 61 (2020). https://doi.org/10.1186/s13046-020-01551-9
- Hepatocellular carcinoma