LncRNA RPSAP52 Upregulates TGF-β1 to Increase Cancer Cell Stemness and Predict Postoperative Survival in Glioblastoma.

INTRODUCTION: Ribosomal protein SA pseudogene 52 (RPSAP52) has been characterized as an oncogenic lncRNA in pituitary tumors. Analysis of TCGA dataset revealed the upregulation of RPSAP52 in glioblastoma (GBM). We, therefore, investigated the roles of RPSAP52 in GBM.
METHODS: A total of 50 GBM patients (33 males and 20 females; 54-75 years; mean age: 61.8±5.8 years) were selected from the 89 cases of GBM patients. Under the guidance of MRI, brain biopsy was performed to collect GBM tissues from each patient for the diagnosis of GBM. U-373 MG cells were employed and had transient transfections. qRNA, Western blot, and a series of experiments were performed to characterize their associations.
RESULTS: The results showed that RPSAP52 was upregulated in GBM patients, and its high expression levels predicted poor survival. In GBM tissues, expression levels of RPSAP52 were significantly and positively correlated with that of TGF-β1. In GBM tissues, RPSAP52 positively regulated the expression of TGF-β1. Cell stemness assay showed that, compared to C and NC groups, overexpression of RPSAP52 and TGF-β1 led to increased, while silencing of RPSAP52 led to decreased CD133+ cells. Overexpression of TGF-β1 attenuated the effects of RPSAP52 siRNA silencing.
CONCLUSION: Therefore, RPSAP52 upregulates TGF-β1 to increase cancer cell stemness and predict postoperative survival in GBM.

Introduction

Glioblastoma (GBM) is a rare type of malignant brain tumor that only affects 3.6 out of 100,000 people worldwide.1 However, GBM has been considered as a major cause of cancer deaths due to its highly aggressive nature and the lack of effective therapeutic approaches.1,2 With active treatment, such as surgical resection, radiation therapy and chemotherapy, the median survival time of GBM patients after the initial diagnosis is 15–18 months, and only less than 5% GBM patients can live long than 5 years.3 Targeted therapies are emerging therapeutic approaches for GBM,4,5 However, the development of novel targeted therapies is limited because of the unclear molecular pathogenesis.6

TGF-β signaling plays critical roles in nearly all aspects of cancer development and progression.7 Inhibition of the TGF-β pathway is considered as a promising therapeutic target for the treatment of different types of cancer including GBM.8,9 It has been well established that the TGF-β signaling may interact with long (>200 nt) non-coding RNAs (lncRNAs) to participate in cancer biology.10,11 Ribosomal protein SA pseudogene 52 (RPSAP52) is a recently characterized oncogenic lncRNA in pituitary tumors,12 while its functions in other types of cancer are unknown. Through the analysis of TCGA dataset, we found that RPSAP52 was only upregulated in GBM tissues but not in non-tumor tissues. This observation triggered our interest to investigate the roles of RPSAP52 in GBM. In this study, we found that RPSAP52 may interact with the TGF-β signaling to regulate GBM cell stemness.

Materials and Methods

GBM Patients

This study passed the ethics review board of the Sixth Medical Center, Chinese PLA General hospital, and was carried out according to the Declaration of Helsinki. A total of 50 GBM patients (33 males and 20 females; 54 to 75 years; mean age: 61.8±5.8 years) were selected from the 89 cases of GBM patients admitted to the aforementioned hospital between May 2012 and July 2014. Inclusion criteria: 1) patients confirmed by histopathological exam; 2) newly diagnosed GBM cases. Exclusion criteria: 1) other clinical disorders were diagnosed; 2) therapies were initiated; 3) recurrent GBM. All patients were informed of the details of this project and signed informed consent.

Biopsy

Under the guidance of MRI, brain biopsy was performed to collect GBM tissues from each patient for the diagnosis of GBM. Adjacent (within 3 cm around the tumors) non-tumor tissues were also collected from each patient. All tissue samples were confirmed by performing a histopathological exam.

Treatment and Follow-Up

Based on patients’ conditions, all patients received surgical resection combined with chemotherapy, radiation therapy or targeted drug therapy. From the day of admission, all patients were followed up for 5 years to record patients’ survival conditions. Only the patients who completed the follow-up or died of GBM were included in the survival analysis.

U-373 MG Cells and Transient Transfections

Human GBM cell line U-373 MG (Sigma-Aldrich, USA) was used as the cell model of GBM. Cells were cultivated in a 5% CO2 incubator at 37 ºC with 95% humidity. Cell culture medium was composed of 10% FBS and 90% DMEM. U-373 MG cells were harvested at 75–85% confluence before transfections. RPSAP52 and TGF-β1 expression vectors were constructed using pcDNA3.1 vector (GeneCopoeia, Guangzhou, China) as the backbone. Negative control (NC) siRNA and RPSAP52 siRNA were designed and synthesized by GenePharma (Shanghai, China). U-373 MG cells were counted and 3×106 cells were transfected with 10 nM vector (empty vector as NC group) or 50 nM siRNA (NC siRNA as NC group) through lipofectamine 2000 (GenePharma)-mediated transient transfection. Following experiments were performed using cells harvested at 24 h post-transfection. Untransfected cells were used as Control (C) cells.

RNA Extractions and qPCR

U-373 MG cells were collected at 24 h post-transfection and cell number was counted. Cell pellets containing 3×104 cells were resuspended in 1 mL Ribozol reagent (Sigma-Aldrich, USA) to extract total RNAs. Biopsies were ground in liquid nitrogen, followed by the addition of Ribozol reagent (Sigma-Aldrich, USA) to extract total RNAs. All operations were performed following the instructions from Sigma-Aldrich. All RNA samples were digested with DNase I (Sigma-Aldrich) at 37 ºC for 90 min to remove genomic DNAs. With ploy (T) as a primer, reverse transcriptions (RTs) were performed using SSRT III system (Thermo Fisher Scientific). With cDNA as a template, qPCR mixtures were prepared using Power SYBR® Green PCR Master Mix (Applied Biosystems) to measure the expression levels of RPSAP52 and TGF-β1 mRNA. GAPDH was used as the endogenous control. All qPCR reactions were performed in a triplicate manner and 2−ΔΔCT method was used for data normalization.

Western Blot

Cell pellets containing 3×104 cells were resuspended in 1 mL RIPA solution (GenePharma) to extract total proteins. Protein samples were quantified using BCA assay (GenePharma). Protein samples were first denatured in a 95°C incubator for 10 min, followed by electrophoresis using 10% SDS-PAGE gel. Proteins were transferred to PVDF membranes, followed by blocking in PBS containing 5% fat-free milk at 25°C for 2 h. After that, TGF-β1 (1:1500, ab92486, Abcam) and GAPDH (1:1500, ab37168, Abcam) rabbit polyclonal primary antibodies were used to incubate with the membranes at 4°C for 18 h. Membranes were further incubated with HRP (IgG) secondary antibody (goat anti-rabbit, 1:2000; ab6721; Abcam) at 25°C for 2 h. Signal production was performed using ECL Substrate (Thermo Fisher Scientific). Expression levels of TGF-β1 were normalized to endogenous control GAPDH using Image J v1.48 software.

Cell Stemness Assay

U-373 MG cells were collected at 24 h post-transfection and were counted after trypsinization. After that, 3×105 cells were incubated with immunoglobulin (Ig) G1-PE (cat. no. 130-093-193; Miltenyi Biotec, Bergisch Gladbach, Germany) or phycoerythrin (PE)-conjugated CD133 (cat. no. 566593, Biosciences, USA) antibody at 4°C for 30 min. Flow cytometry was performed to separate CD133+ cells. Flow cytometry data were processed using the Cell Quest software v5.1 (BD Biosciences, San Jose, CA, USA).

Statistical Analysis

All experiments were performed in 3 biological replicates. Mean values of 3 replicates were calculated and used in the following data analysis. Correlations were explored using linear regression. Survival analysis was performed by dividing the 50 GBM patients into high and low RPSAP52 level groups (n = 25), following by using K-M plotter to plot survival curves and log-rank test to compare survival curves. Paired t test was used to explore differences between non-tumor and GBM tissues. ANOVA (one-way) combined with Tukey’s test were used to explore differences among multiple transfection groups. P < 0.05 was statistically significant.

Results

Upregulation of RPSAP52 Predicted Poor Survival of GBM Patients

The differential expression of RPSAP52 in GBM was first explored using TCGA dataset. It was observed that RPSAP52 can only be detected in GBM tissues but not in non-tumor tissues (0.1 vs 0, ). qPCR was performed to measure the expression levels of RPSAP52 in both GBM and non-tumor tissues from the 50 GBM patients included in this study. Different from TCGA dataset, we detected the expression of RPSAP52 in both non-tumor tissues and GBM tissues. However, expression levels of RPSAP52 were significantly higher in GBM tissues in comparison to non-tumor tissues (Figure 1A, p < 0.0001). Survival curves were plotted and compared using the aforementioned methods. Compared to patients in low RPSAP52 level group, patients in high RPSAP52 level group experienced a significantly lower overall survival rate (Figure 1B).

Figure 1

Upregulation of RPSAP52 predicted poor survival of GBM patients. qPCR was performed to measure the expression levels of RPSAP52 in both GBM and non-tumor tissues from the 50 GBM patients included in this study. Data were compared between two types of cells by paired t test (A). Survival analysis was performed by dividing the 50 GBM patients into high and low RPSAP52 level groups (n = 25), following by using K-M plotter to plot survival curves and log-rank test to compare survival curves (B). PCRs were repeated 3 times and mean values were presented, ***p < 0.0001.

RPSAP52 Significantly and Positively Correlated with TGF-β1 in GBM Tissues

Expression levels of TGF-β1 in both GBM and non-tumor tissues from the 50 GBM patients included were also measured by qPCR. Data were compared between two types of tissues using paired t test. Compared to non-tumor tissues, expression levels of TGF-β1 were significantly higher in GBM tissues (Figure 2A, p < 0.0001). The correlations between TGF-β1 and RPSAP52 were analyzed by linear regression. It showed that expression levels of TGF-β1 were significantly and positively correlated with the expression levels of RPSAP52 across GBM tissues (Figure 2B) but not non-tumor tissues (Figure 2C).

Figure 2

RPSAP52 was significantly and positively correlated with TGF-β1 in GBM tissues. Expression levels of TGF-β1 in both GBM and non-tumor tissues from the 50 GBM patients included in this study were also measured by qPCR. Data were compared between two types of tissues using paired t test (A). Correlations between TGF-β1 mRNA and RPSAP52 across GBM tissues (B) and non-tumor tissues (C) were analyzed by linear regression. PCRs were repeated 3 times and mean values were presented, ***p < 0.0001.

RPSAP52 Positively Regulated TGF-β1 in U-373 MG Cells

U-373 MG cells were transfected with RPSAP52, TGF-β1 expression vectors, or RPSAP52 siRNA to further investigate the interactions between RPSAP52 and TGF-β1. Compared to C and NC groups, expression levels of RPSAP52 and TGF-β1 were significantly altered at 24 h post-transfection (Figure 3A, p < 0.05), indicating successful transfections. Compared to C and NC groups, overexpression of RPSAP52 led to upregulated TGF-β1 at both mRNA (Figure 3B) and protein (Figure 3C) levels (p < 0.05). In addition, silencing of RPSAP52 led to downregulated TGF-β1 at both mRNA (Figure 3B) and protein (Figure 3C) levels (p < 0.05). In contrast, overexpression of TGF-β1 did not affect the expression of RPSAP52 (Figure 3D).

Figure 3

RPSAP52 positively regulated TGF-β1 in U-373 MG cells. U-373 MG cells were transfected with RPSAP52 and TGF-β1 expression vectors as well as RPSAP52 siRNA to further investigate the interactions between RPSAP52 and TGF-β1 mRNA. Overexpression of RPSAP52 and TGF-β1 as well as silencing of RPSAP52 were confirmed by qPCR at 24 h post-transfection (A). qPCR and Western blot were performed to evaluate the effects of overexpressing and silencing of RPSAP52 on the expression of TGF-β1 at mRNA (B) and protein (C) levels, respectively. qPCR was performed to evaluate the effects of overexpressing TGF-β1 on the expression of RPSAP52 (D). Experiments were repeated 3 times and data were expressed as mean values, *p < 0.05.

RPSAP52 Increased U-373 MG Cell Stemness Through TGF-β1

Cell stemness assay was performed to assess the effects of overexpressing RPSAP52 or TGF-β1, or silencing of RPSAP52 on the stemness of U-373 MG cells. Cell stemness assay showed that, compared to C and NC (NC siRNA or empty pcDNA3.1 vector transfection) groups, overexpression of RPSAP52 and TGF-β1 led to increased, while silencing of RPSAP52 led to decreased percentage of CD133+ cells. In addition, overexpression of TGF-β1 attenuated the effects of RPSAP52 siRNA silencing (Figure 4, p < 0.05).

Figure 4

RPSAP52 increased U-373 MG cell stemness through TGF-β1. Cell stemness assay was performed to evaluate the effects of overexpressing RPSAP52 and TGF-β1 as well as silencing of RPSAP52 on the stemness of U-373 MG cells. Experiments were repeated 3 times and data were expressed as mean values, *p < 0.05.

Discussion

This study mainly investigated the expression pattern and functionality of RPSAP52 in GBM. We showed that RPSAP52 was overexpressed in GBM and predicted poor survival of GBM patients. In addition, RPSAP52 may positively regulate TGF-β1 to increase the stemness of GBM cells.

The functions of RPSAP52 have only been investigated in pituitary tumors.12 In pituitary tumors, RPSAP52 is significantly overexpressed and can regulate high mobility group A (HMGA) protein by serving as a miRNA sponge to promote cancer cell proliferation.12 We analyzed TCGA dataset and observed the expression of RPSAP52 in GBM tissues but not in adjacent non-tumor tissues. However, the present study detected the expression of RPSAP52 in non-tumor tissues as well from GBM patients included in this study. This is possibly due to the higher sensitivity of qPCR compared to RNA-seq. Through the analysis of TCGA data, we observed the differential expression of RPSAP52 in most types of cancer, such as colon cancer (0.26 vs 0.07) and thyroid cancer (1.57 vs 0.06). The differential expression may indicate its involvement in these types of malignancies. Future studies are needed to investigate the roles of RPSAP52 in other types of cancer.

TGF-β signaling plays different roles at different stages of cancer development.7,13 For instance, activation of TGF-β inhibits cancer cell proliferation at early stages to play tumor-suppressive roles, while it inhibits cancer cell invasion and migration at advanced cancer stages to promote cancer progression.7,13 Recent studies have also shown that TGF-β is critical for the stemness of cancer cells. For example, the activation of TGF-β has been proved to be critical for the maintenance of glioma-initiating cells.14 In effect, TGF-β is a therapeutic target to decrease cancer cell stemness in high-grade gliomas.15 Consistently, our study also showed the enhancing effects of TGF-β1 in regulating GBM cell stemness. Interestingly, our study showed that RPSAP52 can positively regulate the expression of TGF-β1 to participate in the regulation of the stemness of GBM cells. However, the underlying mechanisms are still unclear. Our preliminary bioinformatics analysis showed that RPSAP52 may sponge multiple miRNAs, such as miR-663, which can target TGF-β1.16 Therefore, RPSAP52 may serve as a miRNA sponge to positively regulate TGF-β1.

In conclusion, RPSAP52 is upregulated in GBM and predicts poor survival. In addition, RPSAP52 can positively regulate TGF-β1 to increase GBM cell stemness.