RASAL2, a RAS GTPase-activating protein, inhibits stemness and epithelial–mesenchymal transition via MAPK/SOX2 pathway in bladder cancer

Muscle-invasive or metastatic bladder cancer (BCa) is associated with a very poor prognosis, and the underlying mechanism remains poorly understood. In this study, we demonstrate RASAL2, a RAS GTPase-activating protein (RAS GAP), acts as a tumor suppressor in BCa. First, RASAL2 was downregulated in BCa specimens and inversely correlated with pathological grades and clinical stages. Furthermore, we observed that RASAL2 could inhibit BCa stemness and epithelial–mesenchymal transition (EMT) based on our gain-of-function and loss-of-function experiments. Mechanistically, we found that mitogen-activated protein kinase/SOX2 signaling had a critical role for maintaining the stemness and mesenchymal properties of RASAL2-deficient BCa cells because inhibition of ERK activity or knockdown of SOX2 could reverse these phenotypes. Also, RASAL2 could inhibit BCa tumorigenesis and distant metastasis in vivo. Moreover, there was an inverse correlation between RASAL2 expression and the stemness/EMT status in subcutaneous xenograft and human BCa specimens. Taken together, our data indicate that RASAL2 is a tumor suppressor in BCa, and modulates cancer stemness and EMT for BCa recurrence and metastasis.

Bladder cancer (BCa) is one of the most common urological cancers worldwide, and most of BCa are transitional cell carcinoma. 1 More than 70% cases are non-muscle-invasive BCa (NMIBC), which can be treated with transurethral resection. However, 450% cases will relapse and 10-30% cases will progress into muscle-invasive BCa (MIBC), and finally result in distant metastasis and poor prognosis of patients. 2 So it is urgent to understand the molecular mechanisms of BCa recurrence and metastasis, and explore novel therapeutic targets to improve patient survival.
Nowadays, cancer stemness has been identified in many cancer types. Cancer stem cells (CSCs) share some properties with somatic stem cells, including self-renewal and multipotent differentiation potential. 3,4 Previous studies have revealed that cancer stemness is responsible for tumorigenicity, therapeutic resistance, relapse and metastasis in BCa. 5,6 Also, accumulating evidences suggest that epithelialmesenchymal transition (EMT) has an important role in the enrichment of cells with CSC properties, which are believed to be the origin of cancer progression. 7,8 Importantly, recent whole-genome studies have discovered the intrinsic basal and luminal MIBC subtypes associated with different chemotherapy sensitivity or patient prognosis, in which a comprehensive molecular alteration involving CSC and EMT has been well characterized. [9][10][11] All these findings provide new insights about invasion and metastasis of BCa.
RASAL2 is a member of RAS GTPase-activating proteins (RAS-GAPs), which can catalyze GTP into GDP and inactivate RAS. 12,13 Several preliminary studies have reported the aberrant expression of RASAL2 in different cancer types, including breast cancer, ovarian cancer, lung cancer and nasopharyngeal carcinoma. [14][15][16][17][18] However, its tumor-suppressive or oncogenic roles in cancer development remain controversial, 16,18 especially, its expression and function in BCa are completely unknown. In this study, we will determine the dysregulation of RASAL2 expression in BCa specimens, and explore its unique roles in modulating stemness and EMT phenotypes of BCa.

Results
RASAL2 is downregulated in human BCa tissues and inversely correlated with tumor grades and clinical stages. To investigate the expression pattern of RASAL2 in BCa tissues, we performed immunohistochemical staining in our cohort including 132 cases of BCa and 10 cases of normal bladder epithelium specimens ( Table 1). As shown in Figures 1a, b and e, we found that the expression of RASAL2 protein was not only lower in BCa tissues than normal epithelia, but also negatively correlated with tumor grades (*Po0.05, **Po0.01). Moreover, the expression of RASAL2 in MIBC tissues was significantly lower than that in NMIBC tissues (Figures 1a, c and e, *Po0.05). Furthermore, we found that RASAL2 expression in high-grade NMIBC tissues are lower than that in low-grade NMIBC tissues (Figures 1d and e, *Po0.05). Also, we analyzed RASAL2 expression in another cohort of BCa tissues profiled by microarray from NCBI GEO. Consistently, we further observed that RASAL2 mRNA was downregulated in human BCa tissues, and there was a negative correlation between RASAL2 and BCa grades and stages (Supplementary Figure 1A, *Po0.05, ***Po0.001), indicating that RASAL2 loss may be associated with BCa progression.
RASAL2 modulates the stemness properties of BCa cells. Cancer stemness is one of the hotspot mechanisms leading to BCa development. 5 Therefore, we applied both gain-of-function and loss-of-function strategies to further  explore the roles of RASAL2 in BCa stemness. As screening in Figure 2a by quantitative RT-PCR and western blotting assays, RASAL2 is highly expressed in SVHUC-1, an immortalized normal bladder epithelial cell line, whereas a relatively lower expression was detected in all BCa cell lines. Fortunately, we successfully established the stable 5637 sublines with endogenous RASAL2 knockdown and T24 sublines with ectopic RASAL2 overexpression ( Figure 2b). Indeed, we noticed that knockdown of RASAL2 expression in 5637 cell lines enhanced the ability of tumorsphere and colony formation, meanwhile overexpression of RASAL2 in T24 cells could decrease this ability (Figures 2c and d Figure 2). These data suggest that RASAL2 could suppress BCa migration and invasion through targeting EMT. RASAL2 regulates stemness and EMT via MAPK/SOX2 signaling pathway in BCa. Next, it is the key to dissect the mechanisms of RASAL2 in regulating BCa stemness and EMT. As RASAL2 is a family member of RAS-GAP, we examined the activities of mitogen-activated protein kinase (MAPK) signaling. Indeed, we found that phosphorylated-ERK (p-ERK) levels increased in 5637 sublines with RASAL2 knockdown but decreased in T24 sublines with RASAL2 overexpression (Figure 4a). So we treated RASAL2-deficient 5637 cells with a specific MEK1/2 inhibitor U0126, and observed that these cells treated with U0126 exhibited a decrease of tumorsphere and colony formation ( RASAL2 modulates BCa cell tumorigenicity and distant metastasis in vivo. To further verify the tumor-suppressive role of RASAL2 in BCa in vivo, we established the subcutaneous xenograft using 5637 sublines. We observed that RASAL2 downregulation resulted in an increased tumor weight and size compared with control (Figures 5a and b). Furthermore, we also next compared the expression of E-cadherin and CD44 in the xenograft tissues by immunohistochemistry (IHC) staining. In consistent with our in vitro results, 5637/shRASAL2 tumors with lower RASAL2 expression presented a higher CD44 staining but a decreased E-cadherin staining compared with control ( Figure 5c, **Po0.01). Also, we further established the tail-vein injection metastasis model using T24-L sublines to explore the role of RASAL2 on BCa metastasis in vivo. As described in our previous studies, 20 tail-vein injection of T24-L sublines could form lung metastasis, whereas overexpression of RASAL2 in T24-L sublines markedly abolished the incidence of lung metastasis in vivo (Figure 5d, Po0.0001). All these results suggested that RASAL2 could modulate tumorigenicity and distant metastasis along with modulation of stemness and EMT in vivo.
RASAL2 is correlated with the status of stemness and EMT in BCa specimens. We also utilized our clinical samples to strengthen our finding in cell lines and xenografts. Indeed, we found that higher CD44 and vimentin expression but lower E-cadherin expression was detected in MIBC than   Figure 6a, RASAL2 is negatively correlated with CD44 and vimentin expression (Pearson's correlation coefficient is − 0.61 between RASAL2 and CD44, − 0.53 between RASAL2 and vimentin), which was also supported by other microarray data (GSE3167) from NCBI GEO (Figure 6b). These data supported RASAL2 as a critical regulator in stemness and EMT of BCa.

Discussion
Patients with MIBC or metastatic BCa have a very poor prognosis, and now fewer therapeutic strategies are available to prolong survival. 21 Therefore, it is crucial and urgent to dissect the molecular mechanisms leading to BCa invasion and metastasis. Herein, we provide evidence that RASAL2 acts as a tumor suppressor in BCa, and modulates the phenotypes of cancer stemness and EMT through MAPK/ SOX2 pathway. RASAL2 belongs to the RAS GAP family members, in which neurofibromin (NF1) and DAB2IP have been most extensively studied. In most instances, NF1 and DAB2IP behave as the classical tumor-suppressor genes. 12 For example, our and other previous studies have demonstrated that DAB2IP has distinct cellular functions such as concurrently modulating different oncogenic pathways associated with cell proliferation, survival, apoptosis and metastasis. [22][23][24] However, the expression and functions of RASAL2 in different cancer types remain largely unknown. Especially, the controversial effects of RASAL2 in breast cancer development have been reported, 16,18 indicating its role as a double-edged sword in certain condition.
Initially, RASAL2 was identified as a potential RAS GAP tumor suppressor in a functional cell-based screen. 25 McLaughlin et al. 18 further demonstrated that RASAL2 appeared to be more frequently silenced by epigenetic mechanisms in breast cancer, and RASAL2 suppression promoted breast tumor development as a consequence of activating K-and H-RAS based on the studies with human xenografts and genetically engineered mouse models. Moreover, in a luminal model of breast cancer, RASAL2 mutations promoted metastasis and correlated with recurrence and poor overall survival of patients with luminal B cancers. Instead, Feng et al. 16 showed that RASAL2 was overexpressed in a subset of triple-negative or estrogen receptor-negative (ERnegative) breast tumors, and activated small GTPase RAC1 signaling to drive mesenchymal invasion and metastasis, indicating its context dependency in breast cancer. In this study, we found that RASAL2 was consistently downregulated either in NMIBC or MIBC tissues, and inversely correlated with pathological grades and clinical stages. Also, by performing gain-of-function and loss-of-function studies in vitro and in vivo, we demonstrated that RASAL2 could inhibit BCa stemness and EMT, which were critical for BCa invasion and metastasis. Moreover, to be consistent with the previous studies reported in luminal breast cancer and ovarian cancer, 17,18 we also observed that RASAL2 could modulate BCa cells stemness and EMT via MAPK pathway, in which the transcription factor SOX2 acted as an important bridge. Indeed, amount of studies have highlighted the central roles of SOX2 in maintaining stemness or cellular plasticity in cancer cells, [26][27][28][29] in particular, it has been reported that SOX2 could be regulated by MAPK signaling for maintaining side population or CSCs in human NMIBC. 30 Taken together, our study revealed a novel mechanism of BCa invasion and metastasis, in which RASAL2 loss could facilitate BCa cell migration, invasion, stemness and tumorigenesis. Furthermore, our findings provide evidence for the critical roles of downstream MAPK/SOX2 signaling in BCa development. Therefore, RASAL2 could be a potential prognostic marker and drug target for BCa diagnosis and treatment.
Plasmid or siRNA transfection or lentiviral infection. RASAL2 shRNA or cDNA plasmids (GenePharma, Shanghai, China) were used to stably silence or overexpress the expression of RASAL2. RASAL2 shRNA backbone was pGPH1. The sequence of RASAL2 shRNAs was as follows: shRASAL2-1: 5ʹ-CAC CGCATGCATCTGTCATGCTTGATTCAAGAGATCAAGCATGACAGATGCATGCTTT TTTG-3ʹ, shRASAL2-2: 5ʹ-CACCGCCAAAGGCCTCTATAGATTCTTCAAGAGAGA ATCTATAGAGGCCTTTGGCTTTTTTG-3ʹ. siRNAs (RiboBio,Guangzhou, China) were used to transiently silence the expression of RASAL2 and SOX2. The sequence of siRNAs for RASAL2 and SOX2 was as follows: siRASAL2-1: 5ʹ-TTT GCTCGTACAACCAGCA-3ʹ, siRASAL2-2: 5ʹ-GGATCGTTGTGGAGAGCAT-3ʹ; Colony formation assay. Cells were seeded into a six-well plate (1000 cells per well), and incubated with fresh medium in a humidified atmosphere at 37°C with 5% CO 2 for 14 days. The plates were washed by phosphate-buffered saline (PBS), fixed in 4% formalin, stained in crystal violet solution for 15 min and then washed with PBS to remove excess dye. The number of colony was counted for each sample.
Migration and invasion assay. Migration and invasion were tested by Boyden chamber assay, obtained from Millipore (Schaffhausen, Switzerland). For migration, the cells were harvested and seeds into the upper chamber (8 μm pore polycarbonate membrane filters) with 3000 cells per well in 0.3 ml serum-free RPMI-1640 for 5637 cells or DMEM for T24, 253J and 253J-BV cells, and 1 ml RPMI-1640 or DMEM medium containing 10% FBS was added to the lower chamber. After 12-h incubation, the upper surface of the chambers were wiped with a Q-tip and fixed with 4% formalin for 15 min, then stained with crystal violet for 15 min followed by washing three times with PBS. The cell number was counted in five random fields in the × 100 magnification. For the invasion assay, the upper chamber was coated with Matrigel (Sigma, St. Louis, MO, USA) and incubated in 37°C with 5% CO 2 , 5 h later, 80 000 cells were seeded into the upper chamber in 0.3 ml serum-free RPMI-1640 for 5637 cells or DMEM for T24, 253J and 253J-BV cells. The plates were incubated for 36 h and the rest of the protocol conducted in a similar manner as the migration assay.
Western blotting. Cells lysates were prepared with RIPA (50 mM Tris, PH 8.0), 150 mM NaCl, 0.1% SDS, 1% NP40 and 0.5% sodium deoxycholate) containing proteinase inhibitors (1% inhibitors cocktail and 1 mM PMSF) (Sigma). Proteins were separated by 12% SDS-PAGE, and transferred onto nitrocellulose membranes. After blocking with 5% non-fat dry milk in TBST, membrane was incubated with primary antibody at 4°C overnight. Then washed by TBST, the membranes were incubated with fluorescent secondary antibodies (Licor, Rockford, IL, USA) at room temperature for 1 h. After washing by TBST, the membranes were scanned by Odyssey Detection system (Licor). Loading differences were normalized using GAPDH.
Tumor xenograft model. Male athymic nude mice were used according to the protocols approved by the ethical committee of Xi'an Jiaotong University. In all, 100 μl serum-free RPMI-1640 medium containing Matrigel (1 : 1, v/v) with 1 × 10 6 5637 sublines (scramble or shRASAL2-2) were injected subcutaneously into both flanks. The mice were killed and xenografts were harvested at day 30, and tumors were weighed and measured for the tumor diameter, and then were stained by immunohistochemistry.
Tail-vein injection metastasis model were generated as described in previous studies. 20,31 Female athymic BALB/c nu/nu mice at age of 4-6 weeks were used according to the protocols approved by the ethical committee of Xi'an Jiaotong university. In all, 5 × 10 6 T24-L cells were injected i.v. via the tail vein. Then, bioluminescence imaging (BLI) was performed to monitor lung metastases with injection of 450 mg/kg D-luciferin substrate (Biosynth, Naperville, IL, USA) in PBS into anesthetized mice. 20 Clinical specimens and immunohistochemistry. All the primary BCa tissues (n = 132) were obtained from the Department of Urology, The First Affiliated Hospital of Xi'an Jiaotong University, Xi'an, China. All samples were used after written consent was obtained from patients. The IHC was performed using the EnVision system (Dako, Carpinteria, CA, USA), Tissue section were de-paraffinized, rehydrated and subjected to 5-min pressure cooker antigen retrieval methods, 10-min of endogenous enzyme block, incubated with the primary antibody overnight at 4°C, 30min Dako Cytomation EnVision-HRP reagent incubation for primary antibodies. Then, the signals were detected by diaminobenzidine (DAB) followed by hematoxylin counterstaining. The result was evaluated according to the intensity of the staining (0, 1 +, 2+ and 3+) and the percentage of positive cells, which were separated by 0 (0%), 1 (1-25%), 2 (26-50%), 3 (51-75%) and 4 (76-100%). Finally, the result is considered by the staining score and the percentage of staining level negative (0 score), weak (1-4 score), moderate (5-8 score) and strong (9-12 score).
Statistical analysis. All the statistical analyses were performed by GraphPad Prism version 5.0 software (GraphPad Software, La Jolla, CA, USA). All data were reported as mean ± S.E.M., and the differences between two groups were compared by the two-tailed Student's t-test. For analyzing gene expression profiling and correlation, we used the data set GSE3167 32 and GSE31684, 33,34 which were downloaded from NCBI GEO database. *Po0.05 was considered statistically significant.

Conflict of Interest
The authors declare no conflict of interest.