Long non-coding RNA PSMA3-AS1 enhances cell proliferation, migration and invasion by regulating miR-302a-3p/RAB22A in glioma

Abstract Glioma is the most prevalent solid tumor in the central nervous system (CNS). Recently, it has been indicated that long non-coding RNAs (lncRNAs) substantially adjust the development of a variety of human cancers. In the present study, it was found and verified via microarray analysis that lncRNA PSMA3-AS1 exhibited a high expression in glioma tissues and cell lines. Then CCK-8, 5-Ethynyl-2′-deoxyuridine (EdU) staining, plate clone assay, Transwell assay, Western blotting and nude mouse model were adopted to verify PSMA3-AS1’s effects on glioma. Knockdown of PSMA3-AS1 inhibited the migration, proliferation and invasion of glioma cells in vivo and in vitro. Besides, PSMA3-AS1 bound to miR-302a-3p directly reduced the expression of miR-302a-3p, thus functioning as an endogenous sponge confirmed by luciferase reporter assay and bioinformatics analysis. PSMA3-AS1 knockdown remarkably enhanced the role of miR-302a-3p overexpression in cell behaviors in glioma. Moreover, these assays also confirmed that RAB22A was a target of miR-302a-3p. In this research, therefore, the PSMA3-AS1/miR-302a-3p/RAB22A pathway regulatory axis may be revealed in the pathogenesis of glioma, and PSMA3-AS1 can be used as an underlying target for the treatment and prognosis of glioma.


Introduction
As the most common solid tumor, Glioma is the most common solid tumor in the central nervous system (CNS) [1,2].Glioma in the CNS has been classified into the low and high grade by the World Health Organization (WHO), based on which low-grade glioma is further sub-divided into Grade I and II, while high-grade glioma into Grade III and IV [3].Generally, cases of Grade IV glioma (glioblastomas; GBM) die within two years, and cases of Grade III anaplastic glioma and Grade II glioma have 2-5-year and 2-year prognosis, respectively [4].Hence, investigating the potential molecular mechanisms of glioma development and progression should be the top priority for enhancing the diagnosis and treatment of the tumor.
Long non-coding RNAs (lncRNAs) are a type of RNAs with length over 200 nucleotides that have no coding ability [5].Besides, multiple biological processes, including cell apoptosis, proliferation and metastasis in cancers are modulated by lncRNAs [6][7][8].There are many lncRNAs related to glioma that have been reported.For example, Xia L et al. confirmed that FER1L4 adjusts the cycle and proliferation of glioma cells [9], Ni W. et al proposed that FoxD2-AS1 adjusts the PI3K/AKT signaling pathway as well as the miR-185-5P/HMGA2 axis, so as to accelerate the progression of glioma [10], and Yang C et al. found that lncRNA HERC2P2 is a tumor suppressor in glioma [11].However, there are still many gaps in our current understanding of lncRNA function.PSMA3-AS1 (ENSG00000257621) is a lncRNA located in Chromosome 14q23.1, and no reports examine the role of PSMA3-AS1 in glioma.
In recent years, it has been widely investigated that lncRNAs interact with miRNAs, Downloaded from http://portlandpress.com/bioscirep/article-pdf/doi/10.1042/BSR20191571/892629/bsr-2019-1571.pdf by guest on 18 September 2020 becoming a mechanism associated with the tumor biology mediated by lncRNAs.MiRNAs are relatively conservative with 18-22 nt in length, which are different from lncRNAs [12].MiRNAs boost the degradation of RNAs in mammals and bind to the 3'-untranslated region (3'UTR) of target mRNAs, thus modulating negative genes [13].Different biological processes contained and are controlled by a variety of miRNAs.By reference to the definition of competing endogenous RNAs (ceRNAs), RNAs compete for shared miRNAs to interact with each other, which implies another method of post-transcriptional regulation [14].Then the miR-302a-3p was selected based on the StarBase database through bioinformatic prediction..Moreover, miR-302a-3p dampens the initiation and development of cancers, such as pancreatic cancer and hepatocellular carcinoma by interaction with various downstream targets [15,16].However, the potential mechanism in which miR-302a-3p functions is still rarely reported.
As such, it was speculated that PSMA3-AS1 contributed to malignant glioma behaviors by adjusting miR-302a-3p/RAB22A.To validate this, the expression pattern of PSMA3-AS1 in glioma tissues and cell lines was firstly studied, followed by investigation of its biological effects and clinical value in glioma.Ultimately, the PSMA3-AS1/miR-302a-3p/RAB22A pathway was validated to be related to the growth, invasion and migration of glioma, so it was considered as an underlying target in the diagnosis and treatment of glioma.

Ethical compliance
The research was approved by the Ethics Committee of Heze NO.3 People's Hospital (Shandong, China).All Population-related study was carried out in accordance with the World Medical Association Declaration of Helsinki, and that all subjects written informed consent.

Clinical tissue samples
Twenty pairs of glioma tissue and adjacent normal tissue samples were extracted from cases who underwent routine surgery with glioma confirmed by pathology in Heze NO.3 People's Hospital from January 2017 to December 2018.All the cases underwent no radiotherapy or chemotherapy prior to surgery, and they signed the informed consent.Immediately after extraction, all tissue samples were frozen in liquid nitrogen for later assays.Detailed information is presented in Table Table1.This research was approved by the Ethics Committee of Heze NO.3 People's Hospital.

Cell lines and cell culture
A172, LN229, U87 and U251, glioma cell lines, and NHA, a human normal astrocyte cell line, were provided by the Cell Bank of CAS (Shanghai, China), purchased from Thermo Fisher Scientific and maintained in DMEM (Gigliomao, Carlsbad, CA, USA) supplemented with 10% FBS (Gigliomao) under the condition of 5% CO 2 at 37°C.

Subcellular fractionation
The nuclear and cytosolic fractions of LN229 and U251 cells were separated using the PARIS Kit (Invitrogen, USA) according to the manufacturer's instructions.RNA was extracted from both fractions.qRT-PCR was then performed using GAPDH as the cytosolic control, and U6 as the nuclear control.

CCK-8 assay
Transfected glioma cells were cultured for 48h and collected in the logarithmic phase.The cells were generated and seeded in96-well plates at a density of 5000 cells per well and cells were cultured at 37°C for 24h.Then 10 μL CCK-8 solution (Solarbio, Beijing, China) was added into each well for 2 h and the absorbance (450 nm) was measured by means of a microplate reader (BioTek Instruments, USA).Finally, the cell growth curve was plotted in accordance with the absorbance at each time point.

5-Ethynyl-2'-deoxyuridine (EdU) staining
Cell-Light EdU DNA Cell Proliferation Kit (RiboBio Co., Ltd, Guangzhou, China) was applied to assess cell proliferation based on the manufacturer's scheme.Specifically, cells were inoculated in 96-well plates at a density of 5,000 cells/well after transfection, and each well were in triplicate.Following culture for 12 h, 4% formaldehyde was added to fix the cells and 100 μl of 1× Apollo® reaction buffer for 30 min of cell incubation.After that, Hoechst 33342 was utilized to stain the cell nuclei for 30 min.The images of cells stained with EdU and Hoechst were observed under an optical microscope (Nikon, Tokyo, Japan), and EdU-positive cells were counted using ImageJ software.

Colony formation assays
In the first place, cell cultures were trypsinized to prepare single-cell suspensions, followed by inoculation into 6-well plates at 150 cells/well and incubation at 37°C for 14 days.Following formation, methanol was added to fix visible colonies and 0.5% crystal violet was used to stain them, followed by cell count in each colony.Colonies formed with over 50 cells were calculated and recorded for statistical analysis.

Cell migration and invasion assays
The 24-well Transwell chamber with polycarbonate membrane (pore size: 8 μm) (EMD Millipore, Billerica, MA, USA) was utilized to evaluate the invasion and migration abilities of glioma cells in vitro.In invasion assay, the cells were added with trypsin after transfection and with medium for suspension, and 20 μL Matrigel (BD Biosciences, San Jose, CA, USA) was used to pre-coat the upper surface of the filter.In migration assay, the cells received the same treatment as the invasion assay except that the upper chamber was added with 2×10 4 cells in medium free of serum, while the lower chamber was added with medium containing 10% FBS.Subsequently, the cells were incubated for 24 h, after which those migrated to the lower chamber were subjected to fixation by 4% paraformaldehyde and dying with 0.05% crystal violet.At least 3 fields of view were selected, in which the stained cells were photographed and calculated.

Luciferase activity assay
Some DNA sequences of PSMA3-AS1 and RAB22A 3'-UTRs containing the binding sites for wild-type (WT) or mutant (MUT) miR-302a-3p were subjected to PCR amplification and cloned into pmirGLO luciferase vectors (Promega, Madison, USA) to produce RAB220A WT, RAB22A MUT, PSMA3-AS1 WT and PSMA3-AS1 MUT reporter plasmids, respectively, so as to elucidate whether PSMA3-AS1 and RAB22A 3'UTRs were directly targeted by miR-302a-3p.During the luciferase assay, the constructed reporter plasmids and miR-302a-3p or miR-NC were applied to independently co-transfect 293T cells by means of Lipofectamine 2000 (Invitrogen, CA, USA), and the Dual-Luciferase® Reporter Assay kit (Promega) was utilized to examine the luciferase activity at 48 h after transfection based on the manufacturer's scheme.

Tumor formation experiment in nude mice
Nude mice at the age of 6 weeks were bought from Cancer Institute of the Chinese Academy of Medical Science.Both sides of the flank area of mice (n=3) in each group were injected with about 1×10 6 U251 that stably transfected by sh-NC, sh-PSMA3-AS1 or sh-PSMA3-AS1+miR-302a-3p inhibitor lentivirus in 0.1 mL PBS.The tumor volume was calculated using a caliper from the 7 th day by reference to the formula V = (L×W 2 )/2 (L: length, W: width of the tumor) every three days.Based on the average tumor volume, the growth curve was drawn in each group.31 days later, mice were intraperitoneally injected with 3% pentobarbital sodium and were killed by excessive anesthesia with a dose of 90 ml/kg, and the tumors were removed for follow-up study.Animal experiments took place in SPF Animal Laboratory at Taishan Medical University.Subsequently, qRT-PCR was carried out to detect the expression of RAB22A or miR-302a-3p.All assays obtained the approval from the Animal Care and Use committee of the Heze NO.3 People's Hospital and Taishan Medical University with the approval number of 201710012-1.And all in vivo experiments were performed in Taishan Medical University.SPSS 22.0 (IBM, SPSS, Chicago, IL, USA) was adopted for statistical processing.The paired-samples t-test was conducted to assess the differences in the expression between glioma tissues and paired adjacent normal tissues, and features of glioma tissues and the expression of PSMA3-AS1 were studied via the χ 2 test.In three independent experiments, data were expressed as mean ± SD.Besides, the Student's t-test and one-way ANOVA were adopted for assessment of intergroup or intragroup differences, respectively.p<0.05 represented a different that was statistically significant.

Characteristic and expression of PSMA3-AS1 in glioma
Through analyzing the GSE103227 microarray, PSMA3-AS1 was found to be highly expressed in glioma.Differentially expressed lncRNAs were identified according to the criteria of p<0.05 and fold change >2.The differently expressed lncRNAs including PSMA3-AS1 with most significance was identified in glioma.To ascertain these results, PSMA3-AS1 expression was assessed in 20 pairs of glioma tissues and adjacent normal tissues via qRT-PCR.In comparison with paired adjacent tissues, PSMA3-AS1 exhibited a notably raised expression in glioma tissues (Figure 1A).The PSMA3-AS1 was elevated in patients with advanced stages of gliomas (III + IV vs I + II, P < 0.05).Detailed information is presented in Table 1.Identically, NHA cells had a lower PSMA3-AS1 expression than human glioma cells (Figure 1B).Notably, LN229 and U251 had high expressions compared with the other glioma cell lines, so the two cell lines were chosen for the loss-of-function assays.We performed the nuclear and cytoplasmic separation experiment and found that PSMA3-AS1 mainly located in the cytoplasm (Figure 1C).This result revealed that PSMA3-AS1 will mainly take part in the post-transcriptional regulation of genes related to glioma.

PSMA3-AS1 down-regulation impeded the migration, proliferation and invasion of glioma cells
Since PSMA3-AS1 in glioma cell lines was raised, whose biological function was further unfolded in vitro.LN229 and U251 were transfected with si-PSMA3-AS1-1, si-PSMA3-AS1-2 or si-PSMA3-AS13, respectively.It was found in qRT-PCR that PSMA3-AS1 expressions in LN229 and U251 were reduced after transfection, among which the sharpest decrease in PSMA3-AS1 expression was caused by si-PSMA3-AS1-2 transfection (Figure 2A).According to CCK-8 assay, the growth of LN229 and U251 was remarkably blocked by si-PSMA3-AS1-2 transfection instead of si-NC transfection (Figure 2B).Subsequently, EdU-positive cells were confirmed to be down-regulated in LN229 and U251 by EdU assay (Figure 2C).The results of the plate clone experiment were consistent with those of the previous experiments, and the colony formation of glioma cells after si-PSMA3-AS1-2 transfection was significantly inhibited (Figure 2D).Functions of PSMA3-AS1 in cell migration and invasion were then examined via Transwell assay.It was found that knockdown of PSMA3-AS1 decreased the number of migrated cells in LN229 and U251 (Figure 3A).Similarly, PSMA3-AS1 down-regulation suppressed the invasion of LN229 and U251 (Figure 3B).As epithelial-mesenchymal transition (EMT) exerts a pivotal effect on glioma cell invasion and migration, the role of PSMA3-AS1 in epithelial characteristics was further investigated.The expressions of ZO-1 (an epithelial marker), Vimentin (a mesenchymal marker) and Snail1 and Twist (EMT-related transcription factors) were measured.As shown in Figure 3C, si-PSMA3-AS1-2 promoted ZO-1 at the mRNA and protein levels, while inhibiting Vimentin, Snail1 and Twist.It can be seen that reduction of PSMA3-AS1 evidently suppresses the migration, proliferation and invasion of glioma cells.

MiR-302a-3p adjusted the proliferation, migration and invasion of glioma cells by targeting PSMA3-AS1's 3'UTR
There has been more and more evidence that lncRNAs are able to adjust miRNAs' effects by sponging them.StarBase revealed that miR-302a-3p and 3'-UTR of PSMA3-AS1 were complementary as they shared binding sites (Figure 4A), which was then validated via the dual luciferase reporter assay.According to the results, it was not the luciferase activity of PSMA3-AS1 MUT but that of PSMA3-AS1 WT that was markedly weakened by miR-302a-3p (Figure 4B).The miR-302a-3p mRNA level was decreased in glioma tissues and was negatively correlated with PSMA3-AS1(Figure 4C and D).At the same time, the content of miR-302a-3p in glioma cell lines was measured by qRT-PCR, which was lower than that of NHA (Figure 4E).Furthermore, miR-302a-3p in LN229 and U251 exhibited a higher expression in PSMA3-AS1 knockdown group than si-NC group (Figure 4F).Subsequently, the roles of PSMA3-AS1/miR-302a-3p in the migration, proliferation and invasion of glioma cells were continued to be discussed.It could be discovered that miR-302a-3p overexpression prominently blocked cell proliferation (Figure 4G) and decreased the cells suffering from metastasis and invasion (Figure 4H), and PSMA3-AS1 knockdown in U251 cells could enhance these effects.

PSMA3-AS1 partially negatively adjusted miR-302a-3p in vivo to cause tumor formation
Findings in vitro were then ascertained by the establishment of a nude mouse model of xenograft.Substantial differences were observed in tumor formation in the sh-NC, sh-PSMA3-AS1 or sh-PSMA3-AS1+miR-302a-3p inhibitor groups (Figure 6A).Additionally, it was demonstrated in the time-dependent analysis that sh-PSMA3-AS1 group had smaller tumor volumes than sh-NC group, and sh-PSMA3-AS1+miR-302a-3p inhibitor could rescue the effect of sh-PSMA3-AS1 as a tumor suppressor.(Figure 6B).Identically, tumor weight in all groups had the same results as the tumor volume (Figure 6C).MiR-302a-3p (Figure 6D) and RAB22A (Figure 6E-F) expressions in xenografts were monitored via qRT-PCR and western blot, the results of which illustrated that PSMA3-AS1 knockdown evidently lowered the expression level of RAB22A but elevated that of miR-302a-3p, which was rescued by simultaneous miR-302a-3p inhibition, in consistent with findings in vitro.

Discussion
In recent years, an enormous body of research has revealed that lncRNAs are not inhibited in malignant tumors and become targets for diagnosing and treating cancers, including glioma [17][18][19].However, there is almost no research related to PSMA3-AS1.Interestingly, it was found in this study that PSMA3-AS1 might be associated with glioma by microarray analysis, which was then further verified and studied in population and cells.Moreover, PSMA3-AS1 was discovered to be remarkably raised in glioma samples and cells, and the reduced PSMA3-AS1 greatly inhibited glioma cell proliferation, migration and invasion, which was verified by assays in vitro and in vivo.These results suggest that PSMA3-AS1 might be a new idea for glioma treatment.
There is a close correlation between lncRNA's biological function and miRNAs and their interactions substantially adjust the progression of numerous cancers [20][21][22].In this research, PSMA3-AS1 was ascertained to directly bind to miR-302a-3p and reduce the expression of the latter in glioma cells.Identical to the results of PSMA3-AS1 knockdown, the migration, proliferation and invasion of glioma cells were notably impeded by miR-302a-3p overexpression, but the function of PSMA3-AS1 knockdown in the glioma cell behaviors was remarkably stimulated by it.Actually, miR-302a-3p suppresses multiple tumors, such as hepatoblastoma [16] and gastric cancer [23], so it was hypothesized that PSMA3-AS1, a host gene of snoRNAs, targets miR-302a-3p to accelerate glioma progression.
It has been found that RAB22A, a member of RAS oncogene family, is markedly raised in numerous human cancers [24].In recent years, epigenetical mechanism has been discovered to adjust RAB22A expression in human cancers, and multiple miRNAs directly target RAB22A to inhibit tumors [25][26][27].Moreover, Yin Y et al. demonstrated that miR-204-5p decreases RAB22A so as to strengthen the chemotherapeutic sensitivity [28].Similarly, it was found that RAB22A was upregulated in glioma.Moreover, RAB22A was regarded as miR-302a-3p's target in glioma examined by bioinformatics analysis.In vivo assays were carried out, which further ascertained that PSMA3-AS1 knockdown was capable of down-regulating RAB22A expression by raising miR-302a-3p, thus slowing down glioma progression.

Figure legends Figure 1 :
Figure legendsFigure 1: Relative PSMA3-AS1 expressions in glioma tissues and cell lines.(A) The relative expressions of PSMA3-AS1 in glioma tissues and adjacent normal tissues (n=20) detected via qRT-PCR.(B) Relative PSMA3-AS1 expressions in one normal cell line and glioma cell lines.(C)Nuclear and cytoplasmic separation experiment of U251 and LN299.Three independent experiments are conducted.Error bars represent mean ± SD of at least three experiments.* p<0.05; ** p<0.01; *** p<0.001.

Figure 2 :
Figure 2: Regulation of PSMA3-AS1 on the growth and proliferation of LN229 and U251.(A) Assessment of the expressions of PSMA3-AS1 in LN229 and U251

Figure 3 :
Figure 3: Regulation of PSMA3-AS1 on the migration and invasion of LN229 and U251.(A-B) The percentage of cell suffering from migration and invasion number in LN229 and U251 after si-PSMA3-AS1-2 transfection is detected.(C) The protein and mRNA levels of ZO-1, Vimentin, Snail1 and Twist in LN229 and U251 upon down-expression of PSMA3-AS1.Three independent experiments are conducted.Error bars represent mean ± SD of at least triplicate experiments.* p<0.05; ** p<0.01; *** p<0.001.