Long Intergenic Noncoding RNA 00261 Acts as a Tumor Suppressor in Non-Small Cell Lung Cancer via Regulating miR-105/FHL1 Axis.

As the third general cancer among human malignant tumors in 2015 1 , the incidence rate of NSCLC remains high both in male and in female worldwide. In the last decades, although molecular targeted therapy is available for NSCLC patients, only a small part of patients' benefits from those discovered driver mutations 2 . Therefore, further exploring the mechanisms of these genomic changes in NSCLC is urgently required. LINC00261 has been shown to exert as a molecular biomarker in several cancers. For example, LINC00261 may be a useful biomarker of the prognosis for HCC patients, based on bioinformatics analysis 3 . Down regulation of LINC00261 may contribute to the carcinogenesis and development of laryngeal carcinoma 4 . However, it remains unknown the roles of LINC00261 in NSCLC. Recently, plentiful studies have revealed the interaction between lncRNAs and miRNAs in several kinds of tumors which probably was associated with many oncogenes that had been known. LncRNA NEAT1 enhances the progression of colorectal cancer via competitively binding miR-34a with SIRT1 5 . SP1-induced lncRNA-ZFAS1 is implicated in the progression of colorectal cancer via miR-150/VEGFA 6 . LncRNA DGCR5 inhibits hepatocellular carcinoma development via regulating the miR-346/KLF14 axis 7 . Shen et al . revealed that LncRNA MEG3 acts as a competing endogenous RNA to modulate HOXA11 via sponging miR-181a in myeloma 8 . LINC00888 enhances tumorigenicity of melanoma through miR-126/CRK 9 . However, the effect of this interaction on NSCLC progression is still unclear. In this study, LINC00261 was found decreased in NSCLC samples. Besides, it inhibited invasion and proliferation in vitro . What's more, we further discovered the interaction between LINC00261 and miR-105 as well as the possible mechanism and confirmed the inhibition effect of LINC00261 on tumorgenesis of NSCLC in vivo .

60 NSCLC patients were enrolled for human tissues who received surgery at Affiliated Hospital of Jiangnan University. Before operation, written informed consent was achieved. The exclusion criteria were prior chemotherapy or radiotherapy. Tumor tissues and adjacent tissues were collected during the surgery and then immediately stored at -80°C. All tissues were analyzed by an experienced pathologist and the clinical data of those patients were recorded and analyzed (as shown in Table 1 ). The overall survival was calculated from the day of primary surgery to death or last follow-up. This study conforms to requirements of the Ethics Committee of Affiliated Hospital of Jiangnan University. A549, SPCA1, H1975, H1299 and PC-9, and 16HBE (human normal bronchial epithelial cell) and 293T embryonic kidney cell line (American Type Culture Collection) were used in this study. Culture medium was consisted of penicillin, RPMI-1640 medium (Thermo Fisher Scientific, USA) and 10% fetal bovine serum (FBS; Invitrogen Life Technologies). Besides, cells were cultured at 37 ˚C in a humidified incubator containing 5% CO 2 . After synthesized, lentiviral small hairpin RNA (shRNA) targeting LINC00261 was cloned into the pLenti-EF1a-EGFP-F2A-Puro vector (BiosettiaInc. San Diego, CA, USA). 293T cells were adopted for packaging LINC00261 lentivirus (LINC00261), the viruses and the empty vector (control). A549 cells were transfected by miRNA mimics and inhibitor provided by Genepharma Co., Ltd. (Shanghai, China). siRNA Universal Negative Control (Sigma-Aldrich) was used as controls. Total RNA was extracted from macrodissected tissue samples or cells by using TRIzol reagent (Invitrogen, USA), and then was reverse-transcribed to cDNAs via reverse Transcription Kit (Takara Biotechnology Co., Ltd., Dalian, China). The β-actin and U6 served as the internal controls for LINC00261 and mi-RNA. Available primers used for SYBR Green real-time PCR were as follows: LINC00261, forwards 5'- ACATTTGGTAGCCCGTGGAG -3' and reverse 5'- TCTTCCCCGGAGAACTAGCA-3'; miR-105, forward 5′- GCCCTCG AGATACCATATCTATCCCCTTTTTCA-3′, and reverse 5′- GCCGAATTC CAACCATGAAGATACGAATTGATG-3′;β-actin, forward 5'-TTGTTACAGGAAGTCCCTTGCC-3' and reverse 5'-ATGCTATCACCTCCCCTGTGTG-3'. U6, forward 5'-CTCGCTTCGGCAGCACATATACT-3' and reverse 5'-ACGCTTCACGAATTTGCGTGTC-3'; Thermal cycle was as follows: 30 sec at 95˚C, 5 sec at 95˚C for 40 cycles, 35 sec at 60˚C. Reagent RIPA (Beyotime) was utilized to extract protein from cells. BCA protein assay kit (Takara) was chosen for quantifying protein concentrations. The target proteins were separated by SDS-PAGE. Then they were incubated with antibodies after replaced to the polyvinylidene fluoride (PVDF) membrane. Rabbit anti-β-actin and rabbit anti-FHL1, as well as goat anti-rabbit secondary antibody were purchased from Cell Signaling Technology (CST, USA). Chemiluminescent film was applied for assessment of protein expression with Image J software. In our study, pGL3 vector (Promega) was used as a backbone for the construction of of 3ʹ-UTR of FHL1 or LINC00261, wild-type (WT) 3ʹ-UTR. Quick-change site-directed mutagenesis kit (Stratagene, Cedar Creek, USA) was adopted for site-directed mutagenesis of miR-105 binding site in FHL1 or LINC00261 3ʹ-UTR, mutant (MUT) 3ʹ-UTR. WT-3ʹ-UTR or MUT-3ʹ-UTR and miR-ctrl or miR-105 mimics was used for cell transfection. 48h later, the dual luciferase reporter assay system (Promega) was utilized to perform the luciferase assays. RIP assay was performed utilizing Magna RIP RNA-Binding Protein Immunoprecipitation Kit (Millipore). Co-precipitated RNAs were detected via RT-qPCR. Proliferating cells were determined by using the 5-ethynyl-2'-deoxyuridine (EdU) abeling/detection kit (Ribobio, Guangzhou, China) according to the manufacturer's protocol. After 48 h of cultivation, the treated cells were incubated with 50 μM EdU labeling medium for 2 h at 37°C under 5% CO 2 and fixed with 4% paraformaldehyde (pH 7.4) for 30 min, and then treated with 0.5% Triton X-100 for 20 min at room temperature. After washing with PBS, staining with anti-EdU working solution was performed at room temperature for 30 min. Cells were incubated with 100 μL Hoechst 33342 (5 μg/mL) at room temperature for 30 min, followed by observation under a fluorescent microscope. The percentage of EdU-positive cells was calculated from five random fields in three wells. After cultured with FBS for 14 days in a six-well plate, all the cells were fixed with methanol and stained with 0.1% crystal violet. Meanwhile, numbers of colonies were counted for comparison. A549 and SPCA1 cell lines were collected for wound healing assay. Firstly, cells were transferred into 6-well plates, and then cultured in DMEM medium overnight. After scratched with a plastic tip, cells were cultured in serum-free DMEM. Wound closure was viewed at 24 h and 48 h later. Each assay was independently repeated in triplicate. 5 ×10 4 cells in 200µL serum-free RPMI-1640 were transformed to top chamber of an insert (8μm pore size; Millipore) which was coated with 50 µg Matrigel (BD Biosciences). And the bottom chamber was added RPMI-1640 and FBS. 48h later, after wiped by cotton swab, the top surface of chambers was immersed for 10 min with precooling methanol. Following were stain in crystal violet for 30 min. Three fields were used to count the data for invasion membrane. The research was approved by the Animal Ethics Committee of Shandong University Animal Center. A549 cells (6×10 5 /mL) transfected with LINC00261 or control were implanted into both axillae of NOD/SCID mice (4-5 weeks old) subcutaneously. Tumor diameters were detected every 5 days. Tumor volume was calculated as the formula (Volume was calculated as length × width 2 × 1/2). Mice were sacrificed and the grafts were removed after 4 weeks. Afterwards, the tumor tissues were got and prepared for Western Blot and qRT-PCR assay to detect the protein level of FHL1 and the RNA level of miR-105 in mice from different groups. SPSS 17.0 (SPSS, Chicago, IL, USA) was utilized to statistical analyze. Results were presented as mean ± standard error of the mean. Chi-square test, ANOVA test and Kaplan-Meier method were selected for clinical data analysis when appropriate. The data in figures were analyzed using a paired Student's t test. p values < 0.05 were considered significant.

Firstly, RT-qPCR was conducted for detecting LINC00261 expression in 60 patients' tumor tissues and 5 NSCLC cells. As the results, LINC00261 was significantly downregulated in tumor tissue samples compared to the adjacent tissues (Fig. 1 a). Meanwhile, LINC00261 level was lower in NSCLC cells than that in normal human bronchial epithelial cell (Fig. 1 b). Analysis of clinicopathological features in those patients demonstrated that downregulated LINC00261 obviously correlated to tumor size, lymph node metastasis and tumor stage (Table 1 ). Moreover, the patients were divided into two groups according to the median expression of LINC00261, and then the survival curves were estimated using the Kaplan-Meier analysis. The results showed that patients had a better disease-free survival with higher LINC00261 level (Fig. 1 c). According to LINC00261 expression in cancer cells, we chose A549 and SPCA1 cancer cells for overexpression of LINC00261. The LINC00261 lentiviruses (LINC00261) and the empty vector (control) were synthetized and transduced into these two cells. Then the LINC00261 expression was determined by qRT-PCR (Fig. 2 a). Furthermore, the results of Edu assay showed that cell proliferation of NSCLC cells was inhibited after LINC00261 overexpressed (Fig. 2 b, c). Meanwhile, cancer formation assay showed that numbers of colonies were less in LINC00261 group than in control group (Fig. 2 d). Then, we performed wound healing assay and transwell assay, and found that overexpressed LINC00261 suppressed NSCLC cell migration and invasion (Fig. 2 e, f). MicroRNAs participate in biological process through their downstream targets. Jin et al . had reported that Mcl-1 was a target gene of miR-105, which contributed into NSCLC progression 11 . In this work, by using StarBase Predicted data, we predicted that FHL1 was another novel downstream target of miR-105 (Fig. 4 a). The luciferase assay revealed that co-transfection of FHL1-WT and miR-105 reduced luciferase activity, which could be partly restored by LINC00261-WT (Fig. 4 b). RT-PCR results further demonstrated that expression of FHL1 in NSCLC cells which transfected with miR-105 mimics was significantly downregulated (Fig. 4 c). Besides, it was found that FHL1 could be inhibited at protein level by miR-105 mimics through western blot assay (Fig. 4 d, e). To explore the interaction between LINC00261 and FHL1, FHL1 expression was detected in both cells and tissues. As the result, FHL1 expression of LINC00261 cells was much higher when compared to that of the control cells (Fig. 4 f). The linear correlation analysis revealed that the FHL1 expression positively correlated to LINC00261 expression in NSCLC tissues (Fig. 4 g). Subcutaneous tumor mice model was established to assess cell growth. A549 cells transfected with LINC00261 or control were injected into mice subcutaneously to induce subcutaneous tumor mice model. Tumor volume was significantly smaller in LINC00261 groups compared with control groups (Fig. 5 a-c). Then tumor tissues were collected to verify the expression of LINC00261 firstly. As shown in Fig. 5 d, the expression of LINC00261 were overexpressed successfully in mice model of LINC00261 group. Further PCR and western blot assay showed the down-regulated of miR-105 and the higher expression of FHL1 both in RNA and protein levels in LINC00261 group (Fig. 5 e-g), which revealed the effect of LINC00261 on miR-105/FHL1 axis.

Evidence has proved that lncRNAs participate in tumorigenesis and development. In the present study, LINC00261 was found downregulated in tissue samples and NSCLC cell lines, and the expression of LINC00261 was significant correlated with tumor size, tumor stage, lymph node metastasis and patients' prognosis. The results were consistent with previous studies, in which decreased expression of LINC00261 is a prognostic marker for NSCLC patients 12 . Zhang et al . found that LINC00261 is a novel prognostic marker for pancreatic cancer 13 . Furthermore, overexpression of LINC00261 expression could inhibit cell proliferation and invasion in NSCLC cells. Moreover, LINC00261 was also confirmed could suppress NSCLC cell growth in vivo . These data indicated that LINC00261 act as a tumor suppressor in NSCLC. Recent studies reveal that LINC00261 participates in tumor progression in many cancers. In esophageal cancer, LINC00261 can induce chemosensitization to 5-fluorouracil via mediating methylation-dependent repression of DPYD 14 . LINC00261 controls endometrial carcinoma progression via regulating miRNA/FOXO1 15 . Moreover, LINC00261 decreases cisplatin resistance and promotes the anti-cancer effect of cisplatin in colon cancer 16 . LINC00261 represses cell growth and migration in endometriosis 17 . Another study revealed LINC00261 inhibits cell proliferation and invasion and promotes cell apoptosis in choriocarcinoma 18 . Here, we demonstrated that LINC00261 was downregulated in NSCLC tissue samples and cells, and significantly correlated with tumor size, tumor stage, lymph node metastasis and patients' prognosis. In vitro, upregulated LINC00261 inhibited cell growth and invaded ability. Latest studies revealed the functions of lncRNAs in NSCLC progression by binding to miRNAs. For example, LncRNA H19 acts as a ceRNA and participates in the development of NSCLC via controlling miR-107 19 . LINC00978 enhances NSCLC cell proliferation and invasion via repressing miR-6754 20 . It has been revealed that knockdown of lncRNA PVT1 promotes radiosensitivity in NSCLC via sponging miR-195 21 . LINC00673 influences NSCLC proliferation, invasion and epithelial mesenchymal transition via sponging miR-150 22 . Recent research found LINC00472 suppresses proliferation and promotes apoptosis of lung adenocarcinoma cells by modulating miR-24-3p/DEDD 23 . Our study revealed that miR-105 could directly bind to LINC00261 through a luciferase assay. And miR-105 was significantly enriched by LINC00261 RIP assay. In addition, miR-105 expression could be downregulated with overexpression of LINC00261. In NSCLC tissues, the expression of miR-105 was downregulated by LINC00261 expression. All the results above suggest that LINC00261 might inhibit tumorigenesis of NSCLC via sponging miR-105. High miR-105 expression was associated with aggressive phenotype of colorectal cancer, and the enhanced expression of miR-105 was required for TNF-α-induced epithelial-mesenchymal transition (EMT) 24 . Cancer-secreted miR-105 destroyed vascular endothelial barriers to promote metastasis 25 . MiR-105 promoted EMT of NSCLC via upregulating Mcl-1 26 . However, miR-105 suppressed cell proliferation and inhibited PI3K/AKT signaling in hepatocellular carcinoma 27 . miR-105 inhibited prostate tumor growth by suppressing CDK6 levels 28 . Reduced miR-105-1 was related to poor survival of patients with NSCLC 29 . Collectively, this inconsistent finding may be due to the tumor type or contamination of normal tissue in the tumor samples or lack of samples or others. Thus, additional studies are required to clarify the role of miR-105 in NSCLC. MicroRNAs participate in biological process through their downstream targets. Jin et al. had reported that Mcl-1 was a target gene of miR-105, which contributed into NSCLC progression 11 . In this work, we identified that FHL1 was another novel downstream target of miR-105. The four and a half LIM domains (FHL) family participates in cell proliferation, apoptosis and differentiation in skeletal and cardiac muscle growth, and the member FHL1, located on chromosome Xq26 30 , 31 . Notably, FHL1 as a tumor suppressor gene is reduced in the malignant cancers, including gastric 32 , 33 , breast, kidney and prostate cancers 34 . Niu et al. confirmed that FHL1 was downregulation and has an inhibitory role in the progression of lung cancer 35 . miRNAs can regulate biological behaviors by targeting one or more genes, which indicates the complicated processes of tumors are owing to not only one molecular but also two or more factors, even a network. Thus, further understanding the regulatory mechanism of multiple molecules in tumors progression is needed.

Above data identified that in NSCLC tissues, LINC00261 was remarkably downregulated and was better related to disease-free survival of NSCLC patients. Besides, LINC00261 could inhibit cell proliferation and invasion through miR-105/FHL1 axis both confirmed in vitro and in vivo . These findings suggested that LINC00261 might contribute to therapy for NSCLC as a candidate target.