- Open Access
Up-regulated LINC01234 promotes non-small-cell lung cancer cell metastasis by activating VAV3 and repressing BTG2 expression
Journal of Hematology & Oncology volume 13, Article number: 7 (2020)
Long noncoding RNAs (lncRNAs) are known to regulate tumorigenesis and cancer progression, but their contributions to non-small-cell lung cancer (NSCLC) metastasis remain poorly understood. Our previous and other studies have revealed the involvement of upregulated LINC01234 in regulating gastric cancer and colon cancer cells proliferation, and we aimed to investigate whether LINC01234 overexpression also contribute to cancer cells metastasis in this study.
We collect the NSCLC tissues and adjacent non-tumor tissues and analyzed expression levels of LINC01234 by quantitative reverse-transcription PCR. LINC01234 were knocked down by using siRNAs or shRNAs, and overexpressed by transfection with overexpression vector; RNA levels of miRNA were downregulated or upregulated with inhibitors or mimics. Transwell assays were used to evaluate cell migration and invasive ability; in vivo metastasis experiments were performed to investigate the effect of LINC01234 on NSCLC cells metastasis. Luciferase reporter, RIP, and ChIP assays were used to determine the regulation of LINC01234 on its targets.
LINC01234 expression is increased in NSCLC tissues, and its upregulation is associated with metastasis and shorter survival in NSCLC. Downregulation of LINC01234 impairs cell migration and invasion in vitro, and inhibits cells metastasis in vivo by acting as a competing endogenous RNA for the miR-340-5p and miR-27b-3p. LINC01234 also interacts with the RNA-binding proteins LSD1 and EZH2, leading to histone modification and transcriptional repression of the anti-proliferative genes BTG2.
Taken together, our findings identify two oncogenic regulatory axes in NSCLC centering on LINC01234: one involving miR-340-5p/miR-27b-3p in the cytoplasm and the second involving EZH2, LSD1, and BTG2 in the nucleus. Our study indicates that these genes may be targeted to reduce or prevent NSCLC metastasis.
Lung cancer is the leading cause of cancer-related death in men world-wide . About 80% of all lung cancer diagnoses are non-small-cell lung cancers (NSCLCs), the major subtypes being lung adenocarcinoma and lung squamous cell carcinoma. Despite rapid advances in diagnostic techniques, molecular-targeted drugs, and immune checkpoint therapy, the 5-year overall survival (OS) time of NSCLC patients remains less than 15% due to lots of patients diagnosed with NSCLC had distant metastases which has been a consistent problem in tumor therapy. Several genetic alterations have been reported to be “drivers” of NSCLC progression, including mutations in epidermal growth factor receptor and mesenchymal–epidermal transition-related genes. Nevertheless, the mechanisms underlying NSCLC progression are not well understood [2,3,4,5]. Therefore, exploration of the molecular mechanisms involved in the progression of NSCLC is critical to improving treatment and the rates of patient survival.
Advanced sequencing techniques and bioinformatics methods have revealed that protein-coding genes only occupy a small proportion (~ 2%) of the whole human genome, while the rest encodes a large number of noncoding RNAs (ncRNAs) including small ncRNAs, pseudogenes, and long noncoding RNAs (lncRNAs) [6,7,8]. Accumulating evidences have uncovered that lncRNAs are pivotal regulators that affect tumorigenesis and cancer progression through affecting almost every aspect of cancer cells behaviors, including cell growth, apoptosis, autophagy, invasion, and metastasis. For example, Wang etc. reported that SATB2-AS1 suppressed colorectal carcinoma aggressiveness through recruiting p300 to the SATB2 promoter and upregulating its expression, which subsequently inhibited Snail transcription and epithelial-to-mesenchymal transition . Our previous studies revealed that overexpression of the lncRNA HOXA11-antisense (AS) facilitates gastric metastasis by interacting with HuR to regulate β-catenin mRNA stability in gastric cancer . As a result, lncRNAs have been highlighted as novel players in tumor metastasis. However, only a small part of these lncRNAs function and underlying mechanisms in cancer progression were clarified, while most of which remains largely unknown.
We speculate that there are still a great number of lncRNAs that are closely related to NSCLC metastasis that have not been explored. In our previous study, we identified a gastric cancer-associated lncRNA LINC01234, which promoted gastric tumorigenesis via sponging miR-204-5p to regulate CBFB expression . Interestingly, a recent study found that LINC01234 could also regulate cell invasion in esophageal cancer . In this study, we analyzed differentially expressed lncRNAs between NSCLC tissues with metastasis and those without metastasis and found that LINC01234 was significantly increased in the metastatic specimen. We further investigated the functional roles and characterized the molecular mechanisms of LINC01234 in NSCLC progression.
Materials and methods
NSCLC sample collection and cell lines
A total of 45 paired NSCLC specimens and adjacent noncancerous tissues were obtained from 45 patients with a histopathological diagnosis of NSCLC who underwent surgery at the Second Affiliated Hospital of Nanjing Medical University between 2010 and 2013. These patients did not receive chemotherapy or radiotherapy prior to surgery. Tissue samples were immediately snap-frozen in liquid nitrogen and stored at − 80 °C until required. This study was approved by the Research Ethics Committee of Nanjing Medical University and informed consent was obtained from all patients.
Four NSCLC adenocarcinoma cell lines (A549, SPC-A1, H1299, and PC9), a NSCLC squamous carcinoma cell line (H226), and a normal human bronchial epithelial cell line (16HBE) were purchased from the Institute of Biochemistry and Cell Biology of the Chinese Academy of Sciences (Shanghai, China). A549, H1299, and H226 cells were maintained in RPMI 1640 basic medium, and SPC-A1, PC9, and 16HBE cells were maintained in DMEM medium with 10% fetal bovine serum, 100 U/ml penicillin, and 100 μg/ml streptomycin (Invitrogen, Carlsbad, CA, USA). All cells were maintained in a humidified 5% CO2 atmosphere at 37 °C.
RNA extraction and quantitative reverse transcription-PCR assays
Total RNA was extracted from tissues or cultured cells using TRIzol reagent (Invitrogen) according to the manufacturer’s instructions. Aliquots of 1 μg RNA were reverse transcribed to cDNA in a final volume of 20 μl under standard conditions using a PrimeScript RT Reagent Kit (TaKara, Dalian, China). Real-time PCR analyses were performed with SYBR Premix Ex Taq (Takara) on an Applied Biosystems 7500 Real-Time PCR System. Expression of LINC01234 and other genes was normalized to that of glyceraldehyde-3-phosphate dehydrogenase (GAPDH) using the relative threshold cycle method, and then converted to fold changes. Specific primers are listed in Additional file 1: Table S1. Primers for miR-27b-3p, miR-340-5p, and U6 were purchased from GeneCopoeia (Rockville, MD, USA). Our quantitative reverse transcription-PCR (qRT-PCR) results were analyzed and expressed relative to threshold cycle (CT) values, and then converted to fold changes.
Plasmid construction and cell transfection
The LINC01234 sequence was synthesized according to the full-length cDNA of human LINC01234 and BTG2 sequences were synthesized according to its coding sequences. All sequences were cloned into the expression vector pCDNA3.1 (Invitrogen). Control and LINC01234-targeting short hairpin RNAs (shRNAs) were purchased from Invitrogen and inserted into the pLKO.1 vector. All final constructs were verified by sequencing. Plasmids were purified using DNA Midiprep Kits (Qiagen, Valencia, CA) and transfected into NSCLC cells using X-treme GENE HP DNA transfection reagent (Roche, Basel, Switzerland). Three LINC01234-targeting small interfering RNAs (siRNAs) obtained from Invitrogen and other gene-targeting siRNAs, miRNA mimics, and miRNA inhibitors (Genepharma, Shanghai, Chian) were transfected into NSCLC cells using Lipofectamine 2000 (Invitrogen) according to the manufacturer’s instructions. Nucleotide sequences for the siRNAs and shRNAs are listed in Additional file 1: Table S1. Cells were harvested 48 h after transfection and analyzed as indicated for the individual experiments.
Luciferase reporter assays
Online bioinformatics databases DIANA Tools (http://carolina.imis.athena-innovation.gr/diana_tools/web/index.php) and miRbase (http://www.mirbase.org/) were used to predict potential miRNA-binding sites in LINC01234. The putative binding sequences were synthesized, inserted into the pGL3-Basic luciferase reporter vector (Promega), and verified by sequencing. Vectors were transfected into human HEK293T cells for 48 h, and luciferase activity was then measured using a Dual Luciferase Kit (Promega), according to the manufacturer’s instructions. The data are presented as relative firefly luciferase activity normalized to the Renilla luciferase activity. All experiments were performed three times.
Cell migration and invasion assays
For cell migration and invasion assays, cells were collected at 48 h post-transfection, and 5 × 104 (for migration assay) or 1 × 105 (for invasion assay) cells in serum-free medium were placed into the upper chamber of an insert (8-μm pore size; Millipore, Billerica, MA, USA). Medium containing 10% fetal bovine serum was added to the lower chamber. After incubation for 24 h, cells in the upper chamber were removed with cotton swabs, and the cells on the lower membrane surface were fixed and stained with 0.5% crystal violet solution. Experiments were independently performed three times.
All protocols were approved by the Committee on the Ethics of Animal Experiments of the Nanjing Medical University and were carried out in strict accordance with the recommendations of the Guide for the Care and Use of Laboratory Animals of the National Institutes of Health. For metastasis assays, SPC-A1 and A549 cells stably transfected with control shRNA or sh-LINC01234 (3 × 106) were injected intravenously via the tail vein. Eight weeks post-injection, the mice were sacrificed and the lungs were removed and photographed. Tumors visible on the lung surface were counted, and the lungs were then stored in formalin.
Cytoplasmic and nuclear RNA were isolated and purified from NSCLC cells using a PARIS Kit (Life Technologies), according to the manufacturer’s instructions.
RNA immunoprecipitation (RIP) assays were performed using an EZ Magna RIP kit (Millipore) using the manufacturer’s protocol. A549 and SPC-A1 cells were lysed in complete lysis buffer, and the extracts were incubated with magnetic beads conjugated with the appropriate specific antibodies or control IgGs (Millipore) for 3–6 h at 4 °C. The beads were washed, incubated with proteinase K to remove proteins, and the purified RNA was eluted and analyzed for the presence of LINC01234 by qRT-PCR. Details of the antibodies and primers are given in Additional file 1: Table S1.
RNA pull-down assays
LINC01234 or control RNAs were transcribed in vitro from pcDNA3.1-LINC01234 using T7 RNA polymerase (Ambion Life) and purified using an RNeasy Mini Kit (Qiagen). One aliquot of transcribed LINC01234 RNA was biotinylated with a Biotin RNA Labeling Mix (Ambion Life). Positive control, negative control, nonbiotinylated, and biotinylated RNAs were incubated with A549 cell lysates. Streptavidin-conjugated magnetic beads were then added and the samples were incubated at room temperature. The beads were then washed, and the eluted proteins were examined by western blot analysis.
Chromatin immunoprecipitation assays
Chromatin immunoprecipitation (ChIP) assays were performed using a MagnaChIP Kit (Millipore) according to the manufacturer’s instructions, as described previously .
Western blot analysis
A549 and SPC-A1 cells were lysed with RIPA extraction reagent (Beyotime) supplemented with a protease inhibitor cocktail (Roche). Proteins in cell lysates were separated by 10% sodium dodecyl sulfate-polyacrylamide gel electrophoresis and transferred to 0.22 μm polyvinylidene fluoride membranes (Millipore). Membranes were probed with specific antibodies using standard methods. Specific protein bands were detected by incubation with ECL chromogenic substrate and quantified by densitometry (Quantity One software; Bio-Rad, Hercules, CA, USA). Antibodies against E-cadherin, N-cadherin, Vimentin, and GAPDH (1:1000) were purchased from Cell Signaling Technology; antibodies against VAV3, EZH2, LSD1, Ago2, and HuR were purchased from Millipore; antibody against BTG2 was purchased from Absin. GAPDH was probed as an internal control. Antibodies are listed in Additional file 1: Table S1.
Statistical analyses were performed using SPSS 20.0 (IBM, Armonk, NY, USA) and Prism software (GraphPad, La Jolla, CA, USA). LncRNA expression levels in primary solid tumors and normal solid tissue samples were compared using the Mann–Whitney U test. For the remaining assays, differences between groups were assessed by paired, two-tailed Student’s t test, Wilcoxon’s test, or χ2 test, as appropriate. Spearman’s correlation analysis was used to calculate the correlations between clinical factors and LINC01234, miR-27b-3p, miR-340-5p, BTG2, and VAV3 expression. All tests were two-sided, and a P value < 0.05 was considered statistically significant.
LINC01234 expression is upregulated in NSCLC and correlates with poor prognosis
We first analyzed lung adenocarcinoma and lung squamous cell carcinoma RNA sequencing datasets from TCGA and found LINC01234 was upregulated in NSCLC tissues compared with adjacent tissues (Fig. 1a). In addition, we found a significant correlation between LINC01234 expression and lung adenocarcinoma stage from TCGA dataset (Fig. 1b). Furthermore, we examined the expression level of LINC01234 in NSCLC tissues and cell lines. qRT-PCR analysis of 45 paired NSCLC and adjacent normal tissues indicated significant upregulation of LINC01234 (fold-change > 1; P < 0.01) in 78% (35/45) of cancerous tissues compared with normal tissues (Fig. 1c). Results also indicated that LINC01234 showed increased expression in NSCLC cell lines in comparison to normal bronchial epithelial cells (Fig. 1d).
Next, we explored the relationship between LINC01234 expression and clinicopathological features of patients with NSCLC. For this, the 45 NSCLC tumor tissues were assigned to two groups (high and low LINC01234 expression) based on the median value (Fig. 1e). Higher LINC01234 expression level was also significantly correlated with tumor size (Chi-square test, P = 0.042), TNM stage (Chi-square test, P = 0.015), and lymph node metastasis (Chi-square test, P = 0.011) in NSCLC patients, but not with other factors, including sex and age (Fig. 1f, g, Additional file 1: Table S2). Univariate Cox regression analyses identified histologic grade, lymphatic metastasis, TNM stage, and LINC01234 expression level as prognostic factors. Other clinicopathologic features such as sex and age were not statistically significant prognosis factors. Multivariate Cox regression analyses showed that LINC01234 expression level were independent prognostic factors for NSCLC patients (Additional file 1: Table S3). Kaplan–Meier survival analysis showed that the OS rates were lower for the high vs. low LINC01234 expression groups, which was supported by Kaplan–Meier Plotter analysis (www.kmplot.com)  (Fig. 1h).
LINC01234 modulates NSCLC cell migration and invasion in vitro and metastasis in vivo
To assess the mechanisms by which LINC01234 promotes the progression of NSCLC, we firstly performed transwell migration and invasion assays. We found that LINC01234 knockdown significantly decreased the migration and invasion ability of A549 and SPC-A1 cells compared with cells expressing si-NC, while overexpression of LINC01234 promoted A549 and SPC-A1 cells migration and invasion (Fig. 2a–c). To validate these results in vivo, we examined the metastatic potential of A549 and SPC-A1 cells stably transfected with empty vector or sh-LINC01234 after injection into nude mice. Consistent with the in vitro analyses, LINC01234 knockdown reduced the number of metastatic lung nodules compared with the control group. H&E staining of excised lung sections confirmed the lower frequency of metastases in the LINC01234-deleted tumors (Fig. 2d, e). However, we did not find metastases in the liver, kidney, intestine, spleen, and other organs (Additional file 1: Figure S1A). Interestingly, expression of the epithelial–mesenchymal transition (EMT) marker E-cadherin was increased, whereas N-cadherin and Vimentin were decreased in LINC01234 downregulated cells (Fig. 2f). In addition, E-cadherin protein levels were also found to be upregulated in A549 and SPC-A1 cells by immunofluorescence assay (Fig. 2g). These in vivo data therefore complement the results of the in vitro functional studies and establish LINC01234 as a regulator of NSCLC metastasis.
LINC01234 interacts with RNA-binding proteins to modulate target gene expression in NSCLC cells
To better understand the underlying mechanism of LINC01234 in NSCLC, we examined its distribution in NSCLC cells via subcellular fractionation. Interestingly, qRT-PCR analysis of RNA purified from nuclear and cytoplasmic cellular fractions revealed that LINC01234 RNA was localized both in nucleus and cytoplasm (Fig. 3a). We also performed ISH to detect the distribution of LINC01234 in NSCLC tissues. The results showed that LINC01234 was expressed both in nucleus and cytoplasm in NSCLC tissues (Additional file 1: Figure S1B). Multiple studies have indicated that lncRNAs could regulate gene expression at transcription level by interacting with RNA binding proteins, such as EZH2, SUZ12, and LSD1. lncRNAs could also affect gene expression at post-transcription level by binding with Ago2. RIP assays confirmed the presence of LINC01234 in EZH2, LSD1, and Ago2 immunoprecipitates from A549 and SPC-A1 cell lysates (Fig. 3b). We next performed RNA pull-down assays to identify LINC01234-associated proteins. We found that biotinylated LINC01234 RNA, but not labeled negative control or antisense RNAs, bound to EZH2, LSD1, and Ago2 (Fig. 3c). Collectively, these data demonstrate that LINC01234 binds to EZH2, LSD1, and Ago2 proteins in NSCLC cells.
LINC01234 functions as a ceRNA for miR-27b-3p and miR-340-5p in NSCLC cells
Emerging evidence suggests that some lncRNAs regulate target gene expression by functioning as ceRNAs, or sponges, for miRNAs, thereby preventing miRNA-mediated regulation of target mRNAs. Interestingly, one of the LINC01234-binding proteins identified above is Ago2, a component of the RNA-induced silencing complex involved in miRNA-mediated repression of mRNAs. Therefore, we investigated the possibility that LINC01234 may also function as a ceRNA. To test this hypothesis, we searched the online bioinformatics databases LncACTdb and miRcode for putative miRNA-binding sites in the LINC01234 sequence. Both databases identified nine putative sites (Fig. 3d), from which we analyzed five with known tumor suppressor functions (miR-27b-3p, miR-340-5p, miR-490-3p, miR-100-5p, and miR-153-5p). Luciferase reporter assays showed that co-transfection with miR-27b-3p and miR-340-5p significantly decreased the luciferase activity of HEK293T cells harboring a LINC01234-driven reporter construct compared with control miRNAs (Fig. 3d). Mutation of the putative miR-27b-3p and miR-340-5p binding sites in LINC01234 abolished their ability to suppress luciferase expression, confirming that the miRNAs specifically interact with LINC01234 (Fig. 3e, f). Consistent with this, Ago2 immunoprecipitates from A549 and SPC-A1 cell lysates were enriched in LINC01234, miR-27b-3p, and miR-340-5p compared with control IgG immunoprecipitates (Fig. 3g). Finally, miR-27b-3p and miR-340-5p levels in A549 and SPC-A1 cells were significantly increased and decreased by LINC01234 depletion and overexpression, respectively (Fig. 3h, i), whereas overexpression of miR-27b-3p or miR-340-5p had no effect on LINC01234 levels (Fig. 3j). These data indicate that LINC01234 acts as a sponge to directly interact with miR-27b-3p and miR-340-5p.
VAV3 is a target of miR-27b-3p and miR-340-5p and is indirectly regulated by LINC01234
To identify LINC01234-regulated target genes in NSCLC in an unbiased manner, we performed RNA-seq to compare the gene expression profiles of A549 cells expressing LINC01234 siRNA or si-NC (Fig. 4a). To identify target genes regulated by both miRNAs and LINC01234, we performed an integrative analysis of the miRDB database and the gene expression profiles of LINC01234-depleted and control A549 cells. Among the 266 genes significantly downregulated by LINC01234 silencing (log2fold-change > 1 and P < 0.05), 12 and 10 were target genes of miR-340-5p and miR-27b-3p, respectively. Interestingly, VAV3 was identified as a target gene for both miRNAs (Fig. 4b). We validated the LINC01234, miR-27b-3p, miR-340-5p, and VAV3 network using luciferase reporter assays. We cloned the wild-type 3′-untranslated region (3′UTR) of VAV3 (wt-VAV3) or one carrying mutations in the miR-27b-3p/miR-340-5p-binding sites (mut-VAV3) into the luciferase vector and co-transfected them with control miRNA or miR-27b-3p/miR-340-5p mimics into HEK293T cells. We found that miR-27b-3p and miR-340-5p significantly decreased luciferase activity driven by wt-VAV3 but not by mut-VAV3 (Fig. 4c), indicating that VAV3 mRNA levels are directly regulated by miR-27b-3p/miR-340-5p in NSCLC cells. Consistent with this, VAV3 mRNA levels in NSCLC cells were reduced or elevated by expression of miR-27b-3p/miR-340-5p mimics or inhibitors, respectively (Fig. 4d). Additionally, E-cadherin was increased, whereas Vimentin was decreased, in VAV3 knockdown cells (Fig. 4e). Knockdown of LINC01234 significantly reduced VAV3 mRNA and protein levels in A549 and SPC-A1 cells, which is consistent with the loss of miR-27b-3p and miR-340-5p sponging (Fig. 4f). In luciferase reporter assays, miR-27b-3p and miR-340-5p-mediated suppression of wt-VAV3-driven luciferase activity was partly reversed by co-expression of the wild-type LINC01234 sequence, but not by LINC01234 mutated in the miR-27b-3p/miR-340-5p binding sites (Fig. 4g). Finally, the suppression of VAV3 protein levels induced by LINC01234 silencing was effectively reversed by co-expression of miR-27b-3p/miR-340-5p inhibitors in NSCLC cells (Fig. 4h). Moreover, the expression levels of miR-27b-3p/miR-340-5p in NSCLC tissues were negatively related with VAV3 expression, while LINC01234 expression was positively related to VAV3 expression. Additionally, miR-27b-3p/miR-340-5p were negatively correlated with LINC01234 (Fig. 4i). Collectively, these data indicate that LINC01234 modulates the expression of VAV3 by acting as a ceRNA for miR-27b-3p and miR-340-5p in NSCLC cells.
miR-27b-3p, miR-340-5p, and VAV3 were involved in NSCLC progression
To determine whether miR-27b-3p and miR-340-5p have tumor suppressor activities, we transfected A549 and SPC-A1 cells with miRNA mimics or inhibitors and examined the effects on cell function (Fig. 5a). Expression of miR-27b-3p and miR-340-5p mimics significantly reduced cell invasion (Fig. 5b). Moreover, co-transfection with miR-27b-3p/miR-340-5p inhibitors partially reversed the inhibition of cell invasion induced by LINC01234 knockdown (Fig. 5c). Next, we silenced VAV3 expression in A549 and SPC-A1 cells by transfection with VAV3 siRNA, which was confirmed to be effective by qRT-PCR (Fig. 5d). Inhibition of VAV3 significantly inhibited NSCLC cell invasion (Fig. 5e). In addition, co-transfection with VAV3 siRNA partially reversed the promotion of A549 and SPC-A1 cell invasion induced by transfection with miR-27b-3p/miR-340-5p inhibitors (Fig. 5f). The increasing in VAV3 protein levels in miR-27b-3p/miR-340-5p inhibitor transfected A549 and SPC-A1 cells was also partly rescued by co-transfection with VAV3 siRNA (Fig. 5g). Collectively, these results suggest that miR-27b-3p, miR-340-5p, and VAV3 are involved in mediating the functional effects of LINC01234 in NSCLC cells.
BTG2 is critical downstream targets of LINC01234
According to the RNA-seq results, we selected eight genes with altered expression involved in cancer cell migration and invasion for further study (Fig. 4a). The identified gene expression changes were validated in A549 and SPC-A1 cells by qRT-PCR. BTG2 is of particular interest because of its remarkable expression fold change upon LINC01234 knockdown (Fig. 6a). Western blot analysis showed that knockdown of LINC01234 significantly increased BTG2 protein levels in NSCLC cells (Fig. 6b), which is consistent with the qRT-PCR results (Fig. 6a) and suggests that BTG2 may play a role in mediating the functional effects of LINC01234. Our findings indicate that, in addition to acting as a ceRNA in the cytoplasm, LINC01234 may also act via EZH2 and LSD1 to regulate transcription of target genes (such as BTG2) in the nucleus. To test this hypothesis, we examined the behavior of A549 and SPCA1 cells expressing EZH2 or LSD1 siRNA. Interestingly, siRNA-mediated knockdown of the LINC01234-binding proteins EZH2 and LSD1 induced BTG2 expression at both the mRNA and protein levels (Fig. 6c, d). Importantly, ChIP assays showed that EZH2 and LSD1 bound directly to the promoter regions of BTG2 and induced deposition of the trimethylated histone 4 lysine 37 (H3K27) and demethylated H3K4 marks (Fig. 6e), which are associated with transcriptional repression. These data suggest that LINC01234 acts, at least in part, by interacting with EZH2 and LSD1, which, in turn, repress transcription of BTG2.
The oncogenic role of LINC01234 is mediated by inhibition of BTG2 expression
Finally, we examined the potential role of BTG2 as a tumor suppressor in NSCLC. In support of this possibility, BTG2 protein was expressed at lower level in tumor tissues compared with normal tissues (Fig. 6f). Kaplan–Meier analysis of patient survival showed that BTG2 expression positively correlated with prognosis, which was supported by Kaplan–Meier Plotter analysis (www.kmplot.com)  (Fig. 6g). We analyzed the correlation between BTG2 and LINC01234 expression, and found that BTG2 expression in paired NSCLC tissues and normal lung tissues was negatively correlated with LINC01234 expression (Fig. 6h).
To better understand the biological function of BTG2, we performed gain-of-function assays. A549 and SPC-A1 cells were transfected with BTG2 overexpression vectors or siRNAs, and the mRNA and protein levels were confirmed to be significantly upregulated or downregulated by qRT-PCR and Western blotting (Fig. 7a, c). Importantly, overexpression of BTG2 led to impairment and invasion by A549 and SPC-A1 cells (Fig. 7b). Interestingly, expression of EMT marker E-cadherin was increased, whereas Vimentin was decreased, in BTG2 overexpressed cells (Fig. 7c). We next performed rescue experiments to determine whether LINC01234-mediated regulation of invasion is dependent on repression of BTG2. A549 cells were cotransfected with si-LINC01234 and si-BTG2, and SPC-A1 cells were cotransfected with LINC01234 and BTG2 vector. Indeed, co-transfection of cells with si-BTG2 partially rescued si-LINC01234-mediated impairment of invasion and EMT process. Conversely, BTG2 overexpression partially rescued the LINC01234-induced promotion of cell invasion and EMT progress (Fig. 7d, e). Taken together, these findings demonstrate that LINC01234 affects NSCLC cell development and progression through the epigenetic repression of BTG2, at least in part.
Multiple factors are involved in tumor progression, and numerous recent studies have implicated lncRNAs as critical regulators of these processes. In the present study, we identified a NSCLC metastasis-associated lncRNA LINC01234, which is highly expressed in metastatic NSCLC tissues and significantly associated with shorter survival time. In addition, modulation of LINC01234 expression revealed its oncogenic activity through promotion of cell migration, invasion, supporting a potential role for LINC01234 dysregulation in NSCLC progression. Indeed, complementary in vivo studies using a mouse model revealed that LINC01234 plays a key role in tumor metastasis.
A great number of recent studies have demonstrated that lncRNAs contribute to cancer progression through numerous mechanisms; for example, by recruiting histone modification enzymes (such as EZH2, SUZ12, and LSD1) that repress or activate gene transcription [15, 16] acting as competing endogenous RNAs (ceRNAs) or sponges to inhibit microRNA (miRNA) activity , interacting with RNA-binding proteins (e.g., STAU1, UPF1, and hnRNPL) to regulate mRNA stability [10, 17, 18] and encoding small active peptide . In the present study, we investigated the molecular mechanisms through which LINC01234 regulated the tumor progression-related behavior of NSCLC cells, and found that LINC01234 interacted with several RNA-binding proteins, including Ago2, EZH2, LSD1 and SUZ12. Accumulating evidence has revealed the existence of a widespread ceRNA interaction network in which lncRNAs compete with miRNAs for binding sites in the 3′-UTR of target mRNAs. For example, HOXA11-AS promotes gastric cancer cell growth by functioning as a ceRNA for miR-1297 , while HOXD-AS1 acts as a ceRNA for miR-130a-3p and facilitates liver cancer metastasis by regulating SOX4 . Here, we showed that LINC01234 is a ceRNA for miR-340-5p and miR-27b-3p and antagonizes their repression of VAV3 protein translation in NSCLC cells. MiR-340 and miR-27b have been reported to have tumor-suppressive functions in multiple cancers. For instance, Li et al. reported that miR-340 inhibits ovarian cancer metastasis via inactivation of NF-x03BA;B1 . Yan et al. found that miRNA-340 inhibits the invasion of esophageal cancer by targeting phosphoserine aminotransferase 1 . MiR-27b has also been shown to inhibit gastric cancer metastasis by targeting NR2F2 , and suppress NSCLC invasion by targeting SP1 . In addition to these findings, our results showed that overexpression of miRNA-340 and miR-27b repressed NSCLC cell invasion by targeting VAV3 expression. The VAV family of guanine nucleotide exchange factors participates in numerous important pathological processes, including oncogenesis and cell transformation. Recent studies showed that VAV3 expression is increased in breast, prostate, and colorectal cancer [25,26,27], and VAV3 promoted cell metastasis in gastric cancer . Consistent with this, we also found that VAV3 was upregulated in NSCLC, and its knockdown inhibited NSCLC cell invasion. These findings suggest that the LINC01234–miR-340-5p/miR-27b-3p–VAV3 axis plays important roles in the progression of NSCLC.
Our data showed that, in addition to functioning as a ceRNA in the cytoplasm, LINC01234 interacts with some well-known histone modification enzymes, such as EZH2, SUZ12, and LSD1, to repress target gene expression (BTG2) in the nucleus. EZH2 and SUZ12 are core subunits of the Polycomb repressive complex 2 (PRC2), which suppresses gene transcription by trimethylating H3K27. In human melanoma cells, loss of EZH2 partially interfered with invasion capacity . LSD1, one of the first discovered protein lysine demethylases, demethylates H3K4me2 to H3K4me1 or H3K4me0 . LSD1 has been found to contribute to invasion and metastasis of luminal breast cancer cells . We propose that LINC01234 acts as a scaffold and recruits EZH2 and LSD1 to the promoter regions of BTG2, thereby repressing its transcription in NSCLC cells. BTG2 is a newly identified tumor suppressor that belongs to the BTG/TOB family, and many studies have revealed that BTG2 is downregulated in various cancers, including breast cancer, osteosarcoma, and bladder cancer. BTG2 inhibited the hepatocellular carcinoma cell invasion and metastasis [32,33,34]. Our data show that BTG2 expression is reduced in NSCLC tissues compared with normal lung tissues, and is associated with shorter patient survival. In NSCLC, overexpression of BTG2 inhibited cell invasion, and rescue experiments confirmed that the oncogenic function of LINC01234 is partly dependent on repression of BTG2 transcription.
In summary, the present study identifies a NSCLC metastasis-associated lncRNA, LINC01234, which is upregulated in human NSCLC and associated with poorer prognosis. LINC01234 exerts its oncogenic function by promoting cell invasion and metastasis, which it achieves by acting as a miRNA sponge in the cytoplasm and a scaffold for histone modification enzymes in the nucleus. These findings advance our understanding of the lncRNA–miRNA–mRNA ceRNA network in NSCLC progression, and suggest that LINC01234 may have utility as a diagnostic marker and/or therapeutic target for NSCLC (Fig. 7f). However, whether LINC01234 displays similar functions, mechanisms of action, and targets in other cancers is unknown and should be investigated in the future.
Availability of data and materials
The dataset(s) supporting the findings of this study are included within the article.
Competing endogenous RNA
Long noncoding RNA
Non-small-cell lung cancer
Polycomb repressive complex 2
Torre LA, Bray F, Siegel RL, Ferlay J, Lortet-Tieulent J, Jemal A. Global cancer statistics, 2012. CA Cancer J Clin. 2015;65:87–108.
Koivunen JP, Mermel C, Zejnullahu K, Murphy C, Lifshits E, Holmes AJ, Choi HG, Kim J, Chiang D, Thomas R, et al. EML4-ALK fusion gene and efficacy of an ALK kinase inhibitor in lung cancer. Clin Cancer Res. 2008;14:4275–83.
Janku F, Stewart DJ, Kurzrock R. Targeted therapy in non-small-cell lung cancer—is it becoming a reality? Nat Rev Clin Oncol. 2010;7:401–14.
Roberts PJ, Stinchcombe TE. KRAS mutation: should we test for it, and does it matter? J Clin Oncol. 2013;31:1112–21.
Sadiq AA, Salgia R. MET as a possible target for non-small-cell lung cancer. J Clin Oncol. 2013;31:1089–96.
Borel C, Gagnebin M, Gehrig C, Kriventseva EV, Zdobnov EM, Antonarakis SE. Mapping of small RNAs in the human ENCODE regions. Am J Hum Genet. 2008;82:971–81.
Pei B, Sisu C, Frankish A, Howald C, Habegger L, Mu XJ, Harte R, Balasubramanian S, Tanzer A, Diekhans M, et al. The GENCODE pseudogene resource. Genome Biol. 2012;13:R51.
Harrow J, Frankish A, Gonzalez JM, Tapanari E, Diekhans M, Kokocinski F, Aken BL, Barrell D, Zadissa A, Searle S, et al. GENCODE: the reference human genome annotation for the ENCODE project. Genome Res. 2012;22:1760–74.
Wang YQ, Jiang DM, Hu SS, Zhao L, Wang L, Yang MH, Ai ML, Jiang HJ, Han Y, Ding YQ, Wang S. SATB2-AS1 suppresses colorectal carcinoma aggressiveness by inhibiting SATB2-dependent snail transcription and epithelial-mesenchymal transition. Cancer Res. 2019;79:3542-56.
Liu Z, Chen Z, Fan R, Jiang B, Chen X, Chen Q, Nie F, Lu K, Sun M. Over-expressed long noncoding RNA HOXA11-AS promotes cell cycle progression and metastasis in gastric cancer. Mol Cancer. 2017;16:82.
Chen X, Chen Z, Yu S, Nie F, Yan S, Ma P, Chen Q, Wei C, Fu H, Xu T, et al. Long noncoding RNA LINC01234 functions as a competing endogenous RNA to regulate CBFB expression by sponging miR-204-5p in gastric cancer. Clin Cancer Res. 2018;24:2002–14.
Ghaffar M, Khodahemmati S, Li J, Shahzad M, Wang M, Wang Y, Li C, Chen S, Zeng Y. Long non-coding RNA LINC01234 regulates proliferation, invasion and apoptosis in esophageal cancer cells. J Cancer. 2018;9:4242–9.
Wan L, Sun M, Liu GJ, Wei CC, Zhang EB, Kong R, Xu TP, Huang MD, Wang ZX. Long noncoding RNA PVT1 promotes non-small cell lung cancer cell proliferation through epigenetically regulating LATS2 expression. Mol Cancer Ther. 2016;15:1082–94.
Gyorffy B, Surowiak P, Budczies J, Lanczky A. Online survival analysis software to assess the prognostic value of biomarkers using transcriptomic data in non-small-cell lung cancer. PLoS One. 2013;8:e82241.
Chen Z, Chen X, Chen P, Yu S, Nie F, Lu B, Zhang T, Zhou Y, Chen Q, Wei C, et al. Long non-coding RNA SNHG20 promotes non-small cell lung cancer cell proliferation and migration by epigenetically silencing of P21 expression. Cell Death Dis. 2017;8:e3092.
Sun M, Nie F, Wang Y, Zhang Z, Hou J, He D, Xie M, Xu L, De W, Wang Z, Wang J. LncRNA HOXA11-AS promotes proliferation and invasion of gastric cancer by scaffolding the chromatin modification factors PRC2, LSD1, and DNMT1. Cancer Res. 2016;76:6299–310.
Kretz M, Siprashvili Z, Chu C, Webster DE, Zehnder A, Qu K, Lee CS, Flockhart RJ, Groff AF, Chow J, et al. Control of somatic tissue differentiation by the long non-coding RNA TINCR. Nature. 2013;493:231–5.
Niknafs YS, Han S, Ma T, Speers C, Zhang C, Wilder-Romans K, Iyer MK, Pitchiaya S, Malik R, Hosono Y, et al. The lncRNA landscape of breast cancer reveals a role for DSCAM-AS1 in breast cancer progression. Nat Commun. 2016;7:12791.
Corley M, Solem A, Phillips G, Lackey L, Ziehr B, Vincent HA, Mustoe AM, Ramos SBV, Weeks KM, Moorman NJ, Laederach A. An RNA structure-mediated, posttranscriptional model of human alpha-1-antitrypsin expression. Proc Natl Acad Sci U S A. 2017;114:E10244–53.
Wang H, Huo X, Yang XR, He J, Cheng L, Wang N, Deng X, Jin H, Wang N, Wang C, et al. STAT3-mediated upregulation of lncRNA HOXD-AS1 as a ceRNA facilitates liver cancer metastasis by regulating SOX4. Mol Cancer. 2017;16:136.
Li P, Sun Y, Liu Q. MicroRNA-340 induces apoptosis and inhibits metastasis of ovarian cancer cells by inactivation of NF-x03BA;B1. Cell Physiol Biochem. 2016;38:1915–27.
Yan S, Jiang H, Fang S, Yin F, Wang Z, Jia Y, Sun X, Wu S, Jiang T, Mao A. MicroRNA-340 inhibits esophageal cancer cell growth and invasion by targeting Phosphoserine aminotransferase 1. Cell Physiol Biochem. 2015;37:375–86.
Feng Q, Wu X, Li F, Ning B, Lu X, Zhang Y, Pan Y, Guan W. miR-27b inhibits gastric cancer metastasis by targeting NR2F2. Protein Cell. 2017;8:114–22.
Jiang J, Lv X, Fan L, Huang G, Zhan Y, Wang M, Lu H. MicroRNA-27b suppresses growth and invasion of NSCLC cells by targeting Sp1. Tumour Biol. 2014;35:10019–23.
Aguilar H, Urruticoechea A, Halonen P, Kiyotani K, Mushiroda T, Barril X, Serra-Musach J, Islam A, Caizzi L, Di Croce L, et al. VAV3 mediates resistance to breast cancer endocrine therapy. Breast Cancer Res. 2014;16:R53.
Lin KT, Gong J, Li CF, Jang TH, Chen WL, Chen HJ, Wang LH. Vav3-rac1 signaling regulates prostate cancer metastasis with elevated Vav3 expression correlating with prostate cancer progression and posttreatment recurrence. Cancer Res. 2012;72:3000–9.
Uen YH, Fang CL, Hseu YC, Shen PC, Yang HL, Wen KS, Hung ST, Wang LH, Lin KY. VAV3 oncogene expression in colorectal cancer: clinical aspects and functional characterization. Sci Rep. 2015;5:9360.
Xie M, Ma T, Xue J, Ma H, Sun M, Zhang Z, Liu M, Liu Y, Ju S, Wang Z, De W. The long intergenic non-protein coding RNA 707 promotes proliferation and metastasis of gastric cancer by interacting with mRNA stabilizing protein HuR. Cancer Lett. 2019;443:67–79.
Zingg D, Debbache J, Schaefer SM, Tuncer E, Frommel SC, Cheng P, Arenas-Ramirez N, Haeusel J, Zhang Y, Bonalli M, et al. The epigenetic modifier EZH2 controls melanoma growth and metastasis through silencing of distinct tumour suppressors. Nat Commun. 2015;6:6051.
Kim SA, Chatterjee N, Jennings MJ, Bartholomew B, Tan S. Extranucleosomal DNA enhances the activity of the LSD1/CoREST histone demethylase complex. Nucleic Acids Res. 2015;43:4868–80.
Hu X, Xiang D, Xie Y, Tao L, Zhang Y, Jin Y, Pinello L, Wan Y, Yuan GC, Li Z. LSD1 suppresses invasion, migration and metastasis of luminal breast cancer cells via activation of GATA3 and repression of TRIM37 expression. Oncogene. 2019;38:7017–34.
Takahashi F, Chiba N, Tajima K, Hayashida T, Shimada T, Takahashi M, Moriyama H, Brachtel E, Edelman EJ, Ramaswamy S, Maheswaran S. Breast tumor progression induced by loss of BTG2 expression is inhibited by targeted therapy with the ErbB/HER inhibitor lapatinib. Oncogene. 2011;30:3084–95.
Li YJ, Dong BK, Fan M, Jiang WX. BTG2 inhibits the proliferation and metastasis of osteosarcoma cells by suppressing the PI3K/AKT pathway. Int J Clin Exp Pathol. 2015;8:12410–8.
Wagener N, Bulkescher J, Macher-Goeppinger S, Karapanagiotou-Schenkel I, Hatiboglu G, Abdel-Rahim M, Abol-Enein H, Ghoneim MA, Bastian PJ, Muller SC, et al. Endogenous BTG2 expression stimulates migration of bladder cancer cells and correlates with poor clinical prognosis for bladder cancer patients. Br J Cancer. 2013;108:973–82.
This work was supported by grants from the National Natural Science Foundation of China (grant numbers 81871871 to Z.W., 81902333 to X.C., 81602013 to F. N., 81802275 to Q.C. and 81703056 to C.W.); Key Research and Development plan (Social development) of science and technology department of Jiangsu Province(No. BE2019760 to Z.W); the Medical Innovation Team Foundation of the Jiangsu Provincial Enhancement Health Project (No. CXTDA2017021 to Z.W.); the Key young medical talents of Jiangsu Province (No. QNRC2016662 to B.L.); the Science and technology development fund of Nanjing Medical University (No.NMUB2018035 to X.C.); and Postgraduate Research & Practice Innovation Program of Jiangsu Province (No. KYCX19_1170 to Z.C.)
Ethics approval and consent to participate
The procedures of this study were approved by the Institutional Review Board of Nanjing Medical University.
Consent for publication
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 1: Figure S1. (A) Metastases to other organs (Kidney, liver, spleen, and intestines). (B) Distribution of LINC01234 in NSCLC tumor tissues. Table S1. Primer, siRNA and shRNA sequences, antibodies. Table S2. Correlation between LINC01234 expression and clinicopathological characteristics of NSCLC patients (n = 45). Table S3. Univariate and multivariate analysis of clinicopathological factors for over-survival in NSCLC patients (n = 45).
About this article
Cite this article
Chen, Z., Chen, X., Lu, B. et al. Up-regulated LINC01234 promotes non-small-cell lung cancer cell metastasis by activating VAV3 and repressing BTG2 expression. J Hematol Oncol 13, 7 (2020). https://doi.org/10.1186/s13045-019-0842-2