Retraction Note: Long noncoding RNA neuroblastoma-associated transcript 1 gene inhibits malignant cellular phenotypes of bladder cancer through miR-21/SOCS6 axis

Bladder cancer (BC) is one of the most common tumors in the urinary system. Noncoding RNAs are considered to take part in cellular phenotypes and are emerging as diagnostic and prognostic biomarkers of BC. The aim of this study is to investigate the clinical significance of neuroblastoma- associated transcript 1 (NBAT1) gene and its effects on malignant cellular phenotypes in BC. NBAT1 gene was low-expressed in BC tissues and cell lines and its low-expression was related with high pathological grade and metastasis of BC. Upregulation of NBAT1 gene depressed cell viability and invasiveness of KK47 and T24 cells and arrested KK47 and T24 cells at G1 stage. In addition, NBAT1 could target silence the expression of miR-21-5p in RNA-induced silencing complex-dependent manner. KK47 and T24 cells with miR-21-5p knockdown showed reduced cell viability, G1-stage arrest, and depressed invasiveness. MiR-21-5p mediates the regulatory effects of NBAT1 on malignant cellular phenotypes of BC cells. Moreover, SOCS6 gene was a target gene of miR-21-5p, and miR-21-5p modulated malignant cellular phenotypes of KK47 and T24 cells through targeted silencing of SOCS6. In conclusion, low-expression of NBAT1 is associated with the progress and metastasis of BC, and NBAT1 inhibits malignant cellular phenotypes through miR-21-5p/SOCS6 axis in BC. Our findings help to elucidate the tumorigenesis of BC, and future study will provide a novel therapeutic target for BC.


Introduction
Bladder cancer (BC) is a malignant tumor originating from bladder mucosa. In the urinary system, BC is the most common malignant tumor in China and the second most common tumor worldwide after prostatic cancer 1,2 . BC mainly included urothelial carcinoma, squamous cell carcinoma, and adenocarcinoma, as well as other rare types such as small cell carcinoma, carcinoid, malignant melanoma, and so on. Bladder urothelial carcinoma is the main type of BC, accounting for 95% of all BC.
BC can be divided into non-muscle-invasive bladder cancer (NMIBC) and muscle-invasive bladder cancer (MIBC) according to whether it invaded the muscular layer of the bladder wall. NMIBC, including Ta, T1, and Tis-stage BC, is also known as superficial BC. MIBC has intruded into the muscular layer of the bladder wall (T2-T4 stage) and is more likely to have lymphatic or distant metastasis. Despite the comprehensive treatment based on surgery, the recurrence rate of MIBC is high and the prognosis is poor 3,4 . Therefore, it is important and necessary to elucidate the underlying mechanism of BC growth and metastasis as well as find new therapeutic targets.
Noncoding RNAs (ncRNAs) consist of long noncoding RNAs (lncRNAs) and short noncoding RNA which include microRNAs, piwi-interacting RNAs, and short interfering RNAs. NcRNAs have become the focus of life science, especially oncology research in recent years. NcRNAs have been confirmed to participate in various complex diseases of human, especially malignant cancers [5][6][7] . It is well known that ncRNAs play important roles in tumorigenesis through modulating multiple important cellular biological phenotypes, such as cell proliferation, invasiveness, chemoresistance, and so on [8][9][10] . Recent studies found that ncRNAs were biomarkers for diagnosis and prognosis of some malignant cancers and might be new therapeutic targets in the future 11,12 .
Neuroblastoma-associated transcript 1 (NBAT1) gene is a newly identified functional lncRNA gene located at chromosome 6p22.3 and is identified and named by Pandey GK in the risk research of neuroblastoma in 2014 13 . Heretofore, the research on NBAT1 and tumorigenesis is rare. Recent studies had found that NBAT1 gene downregulated and acted as a tumor suppressor gene in osteosarcoma and breast cancer 14,15 . However, the expression level and roles in BC remain unclear.
MiR-21-5p originates from 5′ end of pre-miR-21 which is mapped at chromosome 17q23.1. MiR-21-5p was confirmed to be highly expressed and plays its oncogene roles in a variety of tumors, including BC [16][17][18] . For example, miR-21-5p advanced migration and invasion of cervical carcinoma cells through targeting von Hippel-Lindau tumor suppressor (VHL) gene 19 . But the effects of miR-21-5p on malignant cellular phenotypes of BC are not very clear. Wu Y reported that formononetin could inhibit the invasiveness of BC cells and decrease the expression of miR-21 20 , but the correlation of miR-21-5p expression and the growth and metastasis of BC was not certain.
Suppressor of cytokine signaling 6 (SOCS6) gene is located at chromosome 18q22.2 and encodes a protein containing 535 amino acids. SOCS6 protein belongs to a suppressor of cytokine signaling family and is a cytokine-inducible negative regulator of cytokine signaling. SOCS6 gene has been proved to be a tumor suppressor gene in many malignant tumors, including prostate cancer, non-small-cell lung cancer, cervical cancer, and so on [21][22][23] . SOCS6 can control cell signal transduction by inducing ubiquitination degradation of signal protein 24 . There are no reports of SOCS6 and BC. Nevertheless, the function of SOCS6 in metastasis of BC is still unknown.
Together with the previous study that ncRNA NBAT1 could negatively modulate growth and metastasis of osteosarcoma cells through suppression of miR-21 14 , and SOCS6 is targeted by significantly downregulated miR-21-5p in secondary progressive multiple sclerosis 25 , this study will explore the clinical significance of NBAT1 expression, the role of NBAT1 in the regulation of cellular phenotypes, and its molecular mechanism in BC.

Clinical specimens
In total, 76 BC patients in this study were diagnosed and treated in the Department of Urology of Shengjing Hospital between July 2015 and May 2017. BC tissue and the corresponding normal bladder tissue (NBT) specimens were obtained through radical cystectomy and transurethral resection of bladder tumor. All the tissue specimens were confirmed by pathological diagnosis, and the clinical data were complete. All the patients were not treated with radiotherapy or chemotherapy before operation. This study got the approval of the Ethics Committee of Shengjing Hospital, and each participant provided an informed consent.

Cell lines and culture
Human normal bladder epithelial cell line (HCV29), BC cell line KK47 (low-grade NMIBC), and T24 (high-grade MIBC) were obtained from the Cell Resource Center of Chinese Academy of Medical Sciences (Beijing, China). Those cells were cultured in Dulbecco's modified Eagle's medium (DMEM) containing 10% fetal bovine serum (EXcell, Shanghai, China) in a 95% air/5% CO 2 incubator at 37°C.

Cells transfection
The expression plasmid of NBAT1 (pUC-NBAT1) and its negative control (pUC-NC) were synthesized by Cyagen Biosciences Inc. (Santa Clara, USA). The plasmids were transfected into BC cells via HiPerFect reagent (QIAGEN, Hilden, Nordrhein-Westfalen, Germany) in a six-well culture plate according to instructions. The stable transfected cells were selected using Gibco TM Geneticin (Thermo Fisher Scientific, Waltham, MA, USA).

Cell proliferation assay
The cell proliferation was examined using Enhanced Cell Counting Kit-8 (Beyotime, Beijing, China) according to instructions. Briefly, 2000 cells in 100 μl of medium were added into one pore of 96-well plates, addition of 10 μl enhanced CCK-8 solution to each pore, and incubated for 1 h. The absorbance value was detected with a Microplate Reader at 450 nm.

Flow cytometry assay
Cell Cycle Analysis Kit (Beyotime, Beijing, China) was used to examine the cell cycle (Beyotime, Beijing, China) in accordance with instructions. Cells were fixed for 2 h by 70% ethanol, centrifugated, and cleaned up. In total, 0.5 mL of propidium iodide (PI) staining solution was added and resuspended the cells, and then incubated for 30 min. The cell cycle was determined by the FACScan flow cytometry with Diva 8.0 software (Becton Dickinson, Franklin Lakes, NJ, USA).
Cell apoptosis was examined using Annexin V-FITC/PI apoptosis detection kit (Jiancheng, Nanjing, Jiangsu, China) according to instructions. In total, 2 × 10 5 cells were resuspended in 500 µL of binding buffer, added 5 µL of Annexin V-FITC, and 5 µl of PI in order, and incubated at 25°C for 10 min. The apoptosis rate was detected and analyzed by FACScan flow cytometry with Diva 8.0 software. The apoptosis rate was the percent of cells with FITC-Annexin V positive/PI negative in the right lower quadrant.

Cell invasion assay
Cell invasiveness in vitro was assayed using Invasion Chamber with Matrigel Matrix (Corning, Corning, NY, USA). The lower chamber was added to 600 μL of DMEM medium containing 10% FB; the upper chamber was added to 100 μL of serum-free DMEM medium with 5 × 10 4 cells. Then the chamber was cultured at 37°C for 48 h. Cells located at the upper side of the membrane were removed. The membrane was fixed with methanol and stained with Giemsa. Cells invaded to the lower side of the membrane were counted under a microscope in five randomly chosen fields and the average number was calculated.

Western blotting
Protein of cells was extracted using Protein Extraction Kit (Beyotime, Beijing, China), and quantified using Bradford Protein Assay Kit (Beyotime, Beijing, China). Totally, 30 μg of protein were processed, separated by sodium dodecyl sulfatepolyacrylamide gel electrophoresis (SDS-PAGE), transferred to a polyvinylidene fluoride (PVDF) membrane, blocked with nonfat milk, hybridized with SOCS6 and cyclin D1 antibodies (ab197335 and ab134175, Abcam, Cambridge, MA, USA), and incubated with a secondary antibody. The PVDF membrane was treated with BeyoECL Plus reagent (Beyotime, China). Then, the bands were detected and analyzed by Image J software (NIH, Bethesda, MD, USA).

RNA pull-down assay
The biotinylated probes for NBAT1 and miR-21-5p (bio-NBAT1-W and bio-miR-21-5p-W, containing wildtype binding site) as well as their negative control (bio-NBAT1-M and bio-miR-21-5p-M, containing mutant binding site) were synthetized by GenePharma Co. (Shanghai, China). Probes were dissolved in the buffer and incubated with Dynabeads M-280 Streptavidin (Thermo Fisher Scientific, Waltham, MA, USA) for 10 min at 25°C to form probe-coated beads. Those probe-coated beads were incubated with the lysates from KK47 and T24 cells, and eluted with the washing buffer. The pull-down RNAs were detected with qRT-PCR.

RNA immunoprecipitation (RIP) assay
RIP assay was operated using Magna RIP™ RNA-Binding Protein Immunoprecipitation Kit (Millipore Sigma, Burlington, MA, USA) according to previous literature 26 . Briefly, cells were lysed using RIP buffer and incubated with magnetic beads conjugated with human anti-Ago2 antibody or negative control normal IgG. Then the immunoprecipitated RNA was isolated and detected by qRT-PCR to demonstrate the presence of the binding targets.

Luciferase reporter assay
The luciferase reporter plasmids (pmiR-SOCS6-W, containing wild-type binding site; pmiR-SOCS6-M, containing mutant binding site) were synthesized by the Genescript Co. (Piscataway, NJ, USA). HEK-293T cells were co-transfected with the luciferase reporter plasmids and microRNAs, respectively. The Luciferase Reporter Kit (Beyotime, Beijing, China) was applied to detect the luciferase activity 48 h later in accordance with instructions.

Statistical analysis
The SPSS 22.0 (IBM, USA) was used for statistical analysis. All data were expressed as mean ± standard deviation and analyzed with one-way ANOVA and Student's t test. If the P-value is less than 0.05, the difference was considered statistically significant.

NBAT1 gene was low-expressed in BC and correlated with progress and metastasis of BC
The expression of NBAT1 gene in BC specimens was much higher than that in NBT specimens (Fig. 1a). Compared with HCV29 cells, NBAT1 gene was highly expressed in KK47 and T24 cells (Fig. 1b).
Clinical parameters analysis showed that the decreased expression of NBAT1 was correlated with high pathological grade and smoking history, while it was independent of other parameters, including age and gender of BC patients (Table 1). Furthermore, the low-expression of NBAT1 was associated with the muscle invasion and lymph node metastasis of BC patients (Table 1), and the NBAT1 expression in T24 cells was much lower than that in KK47 cells (Fig.1b). These results revealed that NBAT1 gene was involved in the progress and metastasis of BC.

Upregulation of NBAT1 inhibited malignant cellular phenotypes of BC cells
Further, a series of gain-of-function assays was applied to examine the effect of NBAT1 on malignant cellular phenotypes of BC cells. First, KK47 and T24 cells were transfected with pUC-NBAT1 to upregulate the expression of NBAT1 (Fig.1c). And, the enhanced CCK8, flow cytometry, and cell invasion assays showed that upregulation of NBAT1 depressed cell viability of KK47 and T24 cells (Fig.1d), arrested KK47 and T24 cells at G1 stage (Fig.1e), and inhibited invasiveness of KK47 and T24 cells (Fig.1f). Meanwhile, western blotting assay showed that the expression of G1-stage checkpoint protein cyclin D1 was significantly increased when KK47 and T24 cells occurred in G1-stage block (Fig.1g). But, the cell apoptosis did not change significantly in the KK47 and T24 cells with NBAT1 enhancement (Fig.1h).

miR-21-5p was high-expressed in BC and was necessary in NBT1-induced malignant cellular phenotypes of BC cells
The expression of miR-21-5p in BC was higher than that in NBT specimens (Fig. 3a). The expression of miR-21-5p in KK47 and T24 cells was highly expressed compared with HCV29 cells (Fig. 3b). Clinical parameters analysis displayed that the high-expression of miR-21-5p was correlated with high pathological grade, muscle invasion, and lymph node metastasis of BC patients (Table 1). In addition, the miR-21-5p expression in T24 cells was much higher than that in KK47 cells (Fig.1b). These results suggested that miR-21-5p also takes part in the progress and metastasis of BC. AntagomiR-21-5p was transfected into T24 cells to depress the expression of miR-21-5p (Fig. 3c). And, the cellular phenotype detection confirmed that KK47 and T24 cells with miR-21-5p knockdown showed reduced cell viability (Fig. 3d), G1-stage arrest (Fig. 3e), and depressed invasiveness (Fig. 3f). Also, antagomiR-21-5p could decrease the expression level of G1-stage checkpoint protein cyclin D1 (Fig. 3g).
To clarify the effects of miR-21-5p in the regulation process of NBAT1 on malignant cellular phenotypes, agomiR-21-5p was transfected into T24 cells with stably high-expressed NBAT1 (Fig. 4a). Compared with the group of pUC-NBAT1 and agomiR-NC, pUC-NBAT1 and agomiR-21-5p group exhibited the stronger proliferation and invasiveness, as well as the weaker G1-stage arrest, and miR-21-5p enhancement rescued the inhibitory effect of NBAT1 upregulation on cellular phenotypes (Fig. 4b-e). In view of the above, miR-21-5p mediated the regulatory effects of NBAT1 on malignant cellular phenotypes in BC cells.
Furthermore, western blotting assay showed that NBAT1 enhancement upregulated the expression of SOCS6 protein in KK47 and T24 cells, and agomiR-21-5p could rescue this upregulation (Fig. 5d), which proved that NBAT1 positively regulated the expression of SOCS6 through interacting with miR-21-5p.

MiR-21-5p modulated the malignant cellular phenotypes of BC cells by silencing SOCS6 expression
To determine whether the miR-21-5p-mediated malignant cellular phenotypes of BC cells were regulated by SOCS6, T24 cells were co-transfected with antagomiR-21-5p and pS-SOCS6 or pS-NC. Figure 7a demonstrated co-transfection with antagomiR-21-5p and pS-SOCS6 led to an obvious increase of SOCS6 expressions in T24 cells. The cellular phenotype detection found that, compared with antagomiR-21-5p and pS-NC group, the antagomiR-21-5p and pS-SOCS6 group displayed stronger proliferation and invasiveness as well as the weaker G1-stage arrest (Fig. 7b-e).
These findings indicated that miR-21-5p modulated the malignant cellular phenotypes of BC cells by silencing SOCS6 expression.

Discussion
Accumulating research found that lncRNAs were involved in modulation of malignant cellular phenotypes in almost all malignant tumors and could be predictive biomarkers for metastasis and survival in various cancers. For instance, gastric cancer-associated transcript 3 (GACAT3) gene advanced cell proliferation and The relative luciferase reporter assay of HEK 293 T cells co-transfected with pmiR-SOCS6-W or pmiR-SOCS6-M and agomiR-21-5p or agomiR-NC. **P < 0.01 vs. pmiR-SOCS6-W + agomiR-NC group. c The expression of SOCS6 protein in KK47 and T24 cells. *P < 0.05 vs. agomiR-NC group, **P < 0.01 vs. agomiR-NC group; # P < 0.05 vs. antagomiR-NC group, ## P < 0.01 vs. antagomiR-NC group. d The expression of SOCS6 protein in KK47 and T24 cells. **P < 0.01 vs. pUC-NC + agomiR-NC group; ## P < 0.01 vs. pUC-NBAT1 + agomiR-NC group invasion of colorectal cancer cells; 27 LINC00152 was a predictive biomarker of metastasis and survival in various cancers; 28 our previous study found that growth arrestspecific 5 (GAS5) gene inhibited the malignant proliferation and doxorubicin resistance of BC, and was an independent prognostic biomarker for BC 11 .
Till now, there is no research on the expression level and functional role of NBAT1 in BC. In our study, NBAT1 gene was low-expressed in BC tissues and cell lines, and its low-expression was positively related with high pathological grade, lymphatic, and distant metastasis of BC, which suggested that NBAT1 gene was involved in the progress and metastasis of BC. Yang C reported that NBAT1 inhibited the growth and metastasis of osteosarcoma cells 14 . Hu P found that NBAT1 depressed the migration and invasion of breast cancer cells 15 . Nevertheless, it is still unknown whether NBAT1 takes part in the regulation of malignant cellular phenotypes in BC.
To verify the roles of NBAT1 on malignant cellular phenotypes, the expression of NBAT1 was upregulated in BC cells to carry out a series of gain-of-function assays. Upregulation of NBAT1 inhibited cell proliferation and invasiveness of BC cells and arrested them at G1 stage, which showed that NBAT1 enhancement restricted malignant cellular phenotypes of BC cells. Nevertheless, the underlying mechanism is unknown.  The expression of SOCS6 protein in KK47 and T24 cells. *P < 0.05 vs. antagomiR-NC + pS-NC group; ## P < 0.01 vs. agomiR-21-5p + pS-NC group. b The cell viability of KK47 and T24 cells. **P < 0.01 vs. antagomiR-NC + pS-NC group; ## P < 0.01 vs. agomiR-21-5p + pS-NC group. c The cell cycle of KK47 and T24 cells. *P < 0.05 vs. antagomiR-NC + pS-NC group; # P < 0.01 vs. agomiR-21-5p + pS-NC group. d The number of transmembrane KK47 and T24 cells. Scale bars represent 20 μm. **P < 0.01 vs. antagomiR-NC + pS-NC group; ## P < 0.01 vs. agomiR-21-5p + pS-NC group Little attention has been paid to the miR-21-5passociated malignant cellular phenotypes of BC, and its roles of miR-21-5p are still unclear. Therefore, the correlation between the miR-21-5p expression and BC was examined, and the results found that its high-expression was related with progress and metastasis of BC patients. Furthermore, BC cells with miR-21-5p knockdown showed reduced cell viability and invasiveness as well as G1-stage arrest. Our previous research elucidated that miR-21-5p promoted invasion of renal cancer via TCF21-KISS1 pathway 29 .
Some lncRNAs could combine with microRNAs and form complex regulatory networks, and then modulate the expression and function of microRNAs 30,31 . In addition, noncoding RNAs ordinarily form ribonucleoprotein (RNP) complexes with their partner proteins to exert their functions and miRNAs assemble with argonaute (Ago) family proteins into the effector complex called RISC that mediates the target gene silencing 32 . For instance, X-inactive-specific transcript (XIST) could combine and silence miR-152 to exert its tumor-suppressive functions in glioblastoma 33 . Together with online bioinformatics databases and the published reference 14 , there may be the above lncRNA-microRNA regulatory model between NBAT1 and miR-21-5p.
Then, a series of experiments were carried out to test this hypothesis. First, there was a negative correlation between the expression of NBAT1 and miR-21-5p in BC tissues, and NBAT1 enhancement inhibited the expression of miR-21-5p in BC cells. Second, NBAT1 could specifically combine with miR-21-5p and downregulated the miR-21-5p expression. Third, miR-21-5p enhancement rescued the inhibitory effect of NBAT1 upregulation on cellular phenotypes. On these grounds, miR-21-5p mediated the regulatory effects of NBAT1 on malignant cellular phenotypes in BC cells.
It is well known that microRNAs exert their biological functions by targeted regulation of the expression of its target genes 29,34 ; miR-21-5p might also regulate malignant cellular phenotypes through silencing its target gene. Subsequently, this conjecture was confirmed based on the following: (1) SOCS6 was confirmed to be a target gene of miR-21-5p; (2) upregulation of SOCS6 inhibited malignant cellular phenotypes of BC cells; (3) silencing of SOCS6 separately reversed the regulatory roles of miR-21-5p knockdown on malignant cellular phenotypes in BC cells. These findings verified that miR-21-5p modulated the malignant cellular phenotypes of BC cells by silencing SOCS6 expression. Similarly, Li ZB reported that miR-21 and miR-183 could simultaneously silence the expression of SOCS6, and then regulate cell viability and invasiveness of hepatocellular carcinoma cells 35 . Moreover, we found that NBAT1 positively regulated the expression of SOCS6 through interacting with miR-21-5p. Accordingly, NBAT1 inhibited malignant cellular phenotypes of BC through miR-21/SOCS6 axis.
Recent research found that miR-19 could activate JAK2/ STAT3 signaling pathway via silencing SOCS6 and in osteosarcoma and promote osteosarcoma growth in vitro and in vivo, including the reduction of G1-S arrest and an increase of the S phase 36 . JAK2/STAT3 signaling pathway is one of the critical signaling pathways involved in the development and homeostasis in mammals and recent has been reported to involves in the pathogenesis of various cancers including BC 37,38 . Therefore, JAK2/STAT3 signaling pathway might also be involved in suppressing the modulation of NBAT1/miR-21/SOCS6 axis on malignant cellular phenotypes in BC.
In conclusion, low-expression of NBAT1 is associated with the progress and metastasis of BC, and NBAT1 inhibits malignant cellular phenotypes through miR-21-5p/SOCS6 axis in BC. Our findings help to elucidate the tumorigenesis of BC, and future study will provide a novel therapeutic target for BC.