Next Article in Journal
The Impact of Aging on Cardio and Cerebrovascular Diseases
Next Article in Special Issue
Allopregnanolone Alters the Gene Expression Profile of Human Glioblastoma Cells
Previous Article in Journal
The Role of Trio, a Rho Guanine Nucleotide Exchange Factor, in Glomerular Podocytes
Previous Article in Special Issue
When Immune Cells Turn Bad—Tumor-Associated Microglia/Macrophages in Glioma
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

CircSMARCA5 Inhibits Migration of Glioblastoma Multiforme Cells by Regulating a Molecular Axis Involving Splicing Factors SRSF1/SRSF3/PTB

by
Davide Barbagallo
1,
Angela Caponnetto
1,
Matilde Cirnigliaro
1,
Duilia Brex
1,
Cristina Barbagallo
1,
Floriana D’Angeli
1,
Antonio Morrone
2,
Rosario Caltabiano
2,
Giuseppe Maria Barbagallo
2,3,
Marco Ragusa
1,3,†,
Cinzia Di Pietro
1,†,
Thomas Birkballe Hansen
4,5,*,† and
Michele Purrello
1,3,*,†
1
Department of Biomedical and Biotechnological Sciences—Section of Biology and Genetics, University of Catania, 95123 Catania, Italy
2
Department of Medical, Surgical Sciences and Advanced Technologies and Biotechnological Sciences G.F. Ingrassia, University of Catania, 95123 Catania, Italy
3
Multidisciplinary Research Center on Brain Tumors Diagnosis and Therapy, University of Catania, 95123 Catania, Italy
4
Department of Molecular Biology and Genetics (MBG), Aarhus University, 8000 Aarhus C, Denmark
5
Interdisciplinary Nanoscience Center (iNANO), Aarhus University, 8000 Aarhus C, Denmark
*
Authors to whom correspondence should be addressed.
These are Senior Authors.
Int. J. Mol. Sci. 2018, 19(2), 480; https://doi.org/10.3390/ijms19020480
Submission received: 11 January 2018 / Revised: 1 February 2018 / Accepted: 3 February 2018 / Published: 6 February 2018
(This article belongs to the Special Issue Glioma Cell Invasion)

Abstract

:
Circular RNAs (circRNAs) have recently emerged as a new class of RNAs, highly enriched in the brain and very stable within cells, exosomes and body fluids. To analyze their involvement in glioblastoma multiforme (GBM) pathogenesis, we assayed the expression of twelve circRNAs, physiologically enriched in several regions of the brain, through real-time PCR in a cohort of fifty-six GBM patient biopsies and seven normal brain parenchymas. We focused on hsa_circ_0001445 (circSMARCA5): it was significantly downregulated in GBM biopsies as compared to normal brain tissues (p-value < 0.00001, student’s t-test), contrary to its linear isoform counterpart that did not show any differential expression (p-value = 0.694, student’s t-test). Analysis of a public dataset revealed a negative correlation between the expression of circSMARCA5 and glioma’s histological grade, suggesting its potential negative role in the progression to malignancy. Overexpressing circSMARCA5 in U87MG cells significantly decreased their migration, but not their proliferation rate. In silico scanning of circSMARCA5 sequence revealed an enrichment in binding motifs for several RNA binding proteins (RBPs), specifically involved in splicing. Among them, serine and arginine rich splicing factor 1 (SRSF1), a splicing factor known to be a positive controller of cell migration and known to be overexpressed in GBM, was predicted to bind circSMARCA5 by three different prediction tools. Direct interaction between circSMARCA5 and SRSF1 is supported by enhanced UV crosslinking and immunoprecipitation (eCLIP) data for SRSF1 in K562 cells from Encyclopedia of DNA Elements (ENCODE). Consistently, U87MG overexpressing circSMARCA5 showed an increased expression of serine and arginine rich splicing factor 3 (SRSF3) RNA isoform containing exon 4, normally skipped in a SRSF1-dependent manner, resulting in a non-productive non-sense mediated decay (NMD) substrate. Interestingly, SRSF3 is known to interplay with two other splicing factors, polypyrimidine tract binding protein 1 (PTBP1) and polypyrimidine tract binding protein 2 (PTBP2), that positively regulate glioma cells migration. Collectively, our data show circSMARCA5 as a promising druggable tumor suppressor in GBM and suggest that it may exert its function by tethering the RBP SRSF1.

Graphical Abstract

1. Introduction

Circular RNAs (circRNAs) have recently been established as a comprehensive class of RNAs generated by non-linear back-splicing [1,2,3,4,5,6]. CircRNAs are characterized by covalently joined 5′- and 3′-ends: this peculiarity renders them intrinsically resistant to degradation mediated by exonucleases and, consequently, more stable than linear isoform counterparts, both inside and outside the cells [7]. CircRNAs show species-, cell type- and developmental stage-specific expression patterns [8,9,10]: in animals, they are highly enriched in neural tissues [11,12] and furthermore detectable in several biological fluids, thereby useful as non-invasive robust biomarkers [13,14,15,16,17]. Notwithstanding dysregulation of circRNA expression has been highlighted in several types of cancer, most deregulated circRNAs remain functionally undisclosed till today [18,19,20]. Most studies describe circRNAs as microRNA (miRNA) sponges, however, for some of them, a role as modulators of RNA binding proteins (RBPs) has also been suggested [21,22,23,24,25,26,27,28]. To pinpoint functional involvement of circRNAs in glioblastoma multiforme (GBM), we chose twelve circRNAs abundantly expressed in human brain and whose expression is known to be modulated during biological processes linked to oncogenesis (e.g., differentiation, epithelial-to-mesenchymal transition, etc.). In this study, we focused on circSMARCA5, the most downregulated circRNA (among the twelve analyzed) in GBM samples as compared to healthy brain parenchyma, and characterized it as a novel tumor-suppressor, regulating the migration of GBM cells. Our data indicate that circSMARCA5 may exert its function by modulating the RNA binding protein serine and arginine rich splicing factor 1 (SRSF1), a known oncoprotein involved in positive regulation of cell migration [29].

2. Results

2.1. CircRNA Expression Profile

Expression profile of twelve selected circRNAs in a training set of five fresh-frozen GBM samples matched with adjacent non-tumor tissues revealed a significant downregulation of hsa_circ_0001649/circSHPRH (median fold-change = −2.18, p-value = 0.015) and hsa_circ_0001445/circSMARCA5 (median fold-change = −2.42, p-value = 0.017) (Figure 1A). Here, we focus on circSMARCA5, because of its higher abundance in several regions of the brain as well in GBM biopsies with respect to circSHPRH, as reported by circBase and by Song et al. [30]. Downregulation of circSMARCA5 was confirmed in a test set cohort of fifty-six formalin-fixed paraffin embedded (FFPE) GBM biopsies and seven controls (median fold-change = −5.83, p-value < 0.00001). Interestingly, the expression of the linear isoform counterpart did not vary, suggesting a specific dysregulation of the circular isoform (Figure 1B). CircSMARCA5 down expression was also observed in five GBM cell lines with respect to healthy astrocytes, suggesting a contribution of this cell type to circSMARCA5 dysregulation in GBM (Figure S1) (see Section 3).
2.2. circSMARCA5 Expression Is Linked to Progression through Malignancy
Data reported in the previous paragraph were reproducible in an independent dataset made of twenty GBM, three grade III glioma, four grade II glioma, thirteen normal brain cortex and six normal cerebellum [30]: specifically, this dataset showed that a decreased expression of circSMARCA5 matched with an increase in glioma grade malignancy (Figure 2).

2.3. CircSMARCA5 Cloning and Expression in GBM Cells

pcDNA3_circSMARCA5 transfected in U87MG cells produced the circular isoform of SMARCA5, as revealed by Northern blotting (Figure S2A). The introns flanking the two exons of SMARCA5 pre-mRNA involved in circularization contain a 132 nucleotides long conserved inverted repeat (cloned within pcDNA3_circSMARCA5 vector) with a 77% of complementarity that is most probably involved in circSMARCA5 biogenesis (Figure S2B). Based on our analyses, a very small (30 nucleotides long), probably not functional, inverted repeat appeared in genomes of lagomorpha and rodentia for the first time during evolution, whereas a longer, most probably functional, inverted repeat emerged in primate genomes (Table S1). Consistently, circSMARCA5 has not been annotated in mouse according to CircBase, suggesting that circSMARCA5 is primate specific and that biogenesis depends on the long inverted repeat.

2.4. CircSMARCA5 Inhibits GBM Cells Migration without Altering Cell Viability

Overexpression of pcDNA3_circSMARCA5 in U87MG revealed significant decrease of migration rate as compared to negative control (U87MG transfected with pcDNA3) (average fold change = −2.41, p-value = 0.015) (Figure 3A). Viability of U87MG appeared not to be affected under the same experimental conditions (Figure 3B).

2.5. CircSMARCA5 Is Predicted to Function as Modulator of Several RBPs

To investigate the mechanism through which circSMARCA5 may exert its function, we scanned its sequence for motifs predicted to be bound by miRNAs or RBPs. Probably due to the short length (only 269 nucleotides), no significant enrichment of miRNA binding sites was observed. On the contrary, circSMARCA5 appeared to be enriched in RBP binding sites. Among the RBPs, SRSF1 was predicted to bind circSMARCA5 at multiple sites, by three different tools (Table S2). This prediction was also supported by enhanced UV crosslinking and immunoprecipitation (eCLIP) data on K562, available in Encyclopedia of DNA Elements (ENCODE) (Figure 4).

2.6. Splicing of Serine and Arginine Rich Splicing Factor 3 (SRSF3) Is Regulated by CircSMARCA5

To investigate the possible functional role of circSMARCA5 as regulator of splicing, we transfected U87MG cells with pcDNA3_circSMARCA5 vector and analyzed the expression of two isoforms of (mRNA) SRSF3, with or without exon 4 skipped. It is known that the inclusion of exon 4 within SRSF3 pre-mRNA is controlled by SRSF1 during splicing [31]. We observed an increase in the expression specifically of the (mRNA) SRSF3 isoform including exon 4 in U87MG overexpressing circSMARCA5 (Figure 5A). These data are consistent with a significant upregulation of (mRNA) SRSF3 with skipped exon 4 (SRSF3 No Ex4) in the same cohort of GBM biopsies exhibiting a significant downregulation of circSMARCA5 (Figure 5B). Unexpectedly, (mRNA) SRSF3 including exon 4 (SRSF3 Ex4) was also upregulated in GBM biopsies (Figure 5B) but, consistently with data obtained in U87MG cells overexpressing circSMARCA5, the SRSF3 Ex4/SRSF3 No Ex4 ratio positively correlated with the expression of circSMARCA5 in GBM biopsies (Figure 5C) (see Section 3).

3. Discussion

CircRNAs are on the forefront of basic and applied research but, notwithstanding the accumulating evidence of their potential interest as critical regulators of gene expression, only a few studies describe the molecular mechanisms through which they work. In this paper, we focused on hsa_circ_0001445/circSMARCA5, an exonic 269-nucleotide-long circRNA, highly enriched in human brain (http://www.circbase.org/). We found that its expression decreases in GBM cells with respect to normal brain parenchyma, independently from the expression level of the linear isoform. Moreover, astrocytes (which is the most representative cell type in GBM biopsies) most likely contribute to altered expression of circSMARCA5 in GBM (Figure S1), although it is not possible to exclude that also other cell types—e.g., Tumor Associated Macrophages (TAMs)—may impact on circSMARCA5 dysregulation. This result matches with a previously described increase of circular to linear transcript ratios, observed during differentiation of cardiomyocytes [32] and, together with the inverse relationship between circSMARCA5 expression and glioma malignancy grade, this collectively suggests a functional involvement in GBM. These data, together with Receiver Operating Characteristic (ROC) curves data (Figure S3), also suggest circSMARCA5 as a candidate GBM biomarker, although the use of the term biomarker is referred to a scanty cohort of patients at present and should be re-evaluated in a wider case-control study. To characterize the effects of circSMARCA5 expression, we first cloned its sequence, together with a segment of the introns upstream and downstream the exons that circularize, into an expression vector. Flanking intronic sequences contained an inverted repeat of 132 nucleotides showing 77% complementarity: these sequences are most probably responsible for the biogenesis of circSMARCA5 and are conserved among primates (Figure S4), suggesting a primate-specific expression of this circRNA. According to functional analysis, circSMARCA5 negatively regulates GBM cell migration. Based on further in vitro and in silico characterization, we also speculate that circSMARCA5 may exert its function by disrupting splicing of GBM cells, specifically by modulating the SRSF1 splicing factor. This hypothesis is supported by eCLIP data from ENCODE (Figure 4) as well as by the circSMARCA5-mediated changes in (mRNA) SRSF3 splicing pattern (Figure 5A) and by data on GBM biopsies (Figure 5B,C). Based on our data of circSMARCA5 modulation (Figure 5A) and according to Jumaa et al. and Jihua et al. [31,33], we hypothesize that circSMARCA5 downregulation in GBM biopsies, associated with upregulation of SRSF1 (Figure S5), may promote the skipping of exon 4 in SRSF3 pre-mRNA. At the same time, increased levels of functional SRSF3 protein may determine the increased expression of (mRNA) SRSF3 containing exon 4 (which results in a non-productive non-sense mediated decay (NMD) substrate [34,35]), in a self-regulatory manner (Figure 5B). Despite the increased expression of both (mRNA) SRSF3 isoforms in GBM biopsies compared to normal brain parenchyma, we observe a significant positive correlation between SRSF3 Ex4/SRSF3 No Ex4 ratio and circSMARCA5 expression (Figure 5C). These data support our hypotheses that circSMARCA5 may indirectly regulate the expression of (mRNA) SRSF3 isoforms, by tethering SRSF1, and, in turn, functional SRSF3 protein (synthesized from SRSF3 No Ex4) may positively regulate the expression of SRSF3 Ex4, which is normally expressed at very low levels within cells (Figure S6). Often the networks involved in gene regulation within cells are intricate, so we do not exclude that other nodes (in addition to circSMARCA5) may contribute in regulating the expression of (mRNA) SRSF3 isoforms within GBM cells. Interestingly, SRSF1 is known to regulate a plethora of biological functions, further to splicing, within cells [36] and is located in the middle of a self-regulatory network involving several splicing factors such as PTBP1 and SRSF3 that act as known oncoproteins in different types of cancer [37]. Most specifically, splicing pattern of (mRNA) SRSF3 is defined by SRSF1 and SRSF3 itself regulates the expression of PTBP1 whose overexpression is known to positively regulate GBM cells migration [38]. Conclusively, circSMARCA5 is significantly downregulated in GBM and it may exert its function by modulating the activity of SRSF1 with subsequent effects on SRSF3 and PTBP1 splicing and expression. An in-depth study of the interaction between circSMARCA5 and SRSF1 and its downstream network within GBM cells and to strengthen the possible use of circSMARCA5 as GBM biomarker by extending the case-control study, will represent the natural continuation of this project.

4. Materials and Methods

4.1. GBM Specimens and Cell Lines

GBM biopsies from five fresh-frozen (training set) and fifty-six FFPE (test set) samples were diagnosed by at least three experienced pathologists, according to the 2007 WHO criteria [39]. Each sample contained more than 90% of cancer tissue. Normal brain parenchyma was obtained from area surrounding the tumor that appeared negative to 5-aminolevulinic acid (5-ALA) fluorescence during surgery, from seven patients. Written informed consent was received from participants prior to inclusion in the study. Absence of infiltration of cancer cells within normal brain parenchyma was also confirmed by microscopy analysis performed by at least three experienced pathologists. Additionally, the commercially available FirstChoice® Human Brain Reference RNA (Ambion, Austin, TX, USA) and Human Astrocyte Total RNA (ScienCell Research Laboratories®, San Diego, CA, USA) were used as control healthy tissues. Age, sex and clinical features of the analyzed tumor samples and healthy controls are reported in Table 1. Human GBM cell lines A172, CAS1, DBTRG, SNB19 and U87MG were cultured as previously described [19].

4.2. Immunohistochemistry

Immunohistochemical analysis, incubation with biotinylated anti-rabbit secondary antibody and immunoreaction visualization were performed as previously described [40]. Rabbit polyclonal antibody against GFAP (DAKO, Glostrup, Denmark) was used at a dilution 1:7000. GFAP-positive (GFAP+) cells percentage was evaluated in the highest immunoreactivity fields. It was determined by dividing the number of positive staining cells by 1000 cells. Cells were considered to be positive if there was any cytoplasmic staining present.

4.3. Constructs and Transfections

A 1450 bp PCR amplified fragment of the human genome, comprising exons 15 and 16 of SMARCA5 (the two exons that circularize generating circSMARCA5), the intron between them and 440 bp upstream exon 15 and 546 bp downstream exon 16, was inserted in pcDNA3 vector (Figure S7) after digestion with BamHI and ApaI. U87MG cells were transfected with Lipofectamine 2000 (Thermo Fisher Scientific, Waltham, MA, USA), according to supplier’s protocol. Transfection efficiencies of more than 80% were achieved.

4.4. RNA Extraction

Within 30 min from surgical resection, biopsies were collected in sterile tubes on ice, washed in cold sterile PBS to eliminate any blood residue and stored in RNAlater® Stabilization Solution (Thermo Fisher Scientific) for 16 h at 4 °C. Biopsies were then either fresh-frozen or formalin-fixed and paraffin embedded and stored until their processing. RNA from fresh-frozen and FFPE biopsies was extracted by using TRIzol® and RecoverAll® kits (Thermo Fisher Scientific), respectively, according to manufacturer’s instructions and quantified both through spectrophotometer and Qubit™ fluorometer (Thermo Fisher Scientific).

4.5. CircRNA Candidates Selection and Primer Design

Manual curation of scientific papers, allowed us to select candidate circRNAs, based on their association with at least one of the following biological processes: (i) regulation of neuronal differentiation; (ii) epithelial-to-mesenchymal (EMT) transition; and (iii) cell proliferation (specifically referred to GBM cells or as general process). All selected circRNA transcripts had to be expressed (and preferentially enriched) in human brain, according to RNA Seq data stored in circBase (http://www.circbase.org/). Candidate circRNAs assayed in this study are reported in Table 2. Both convergent (detecting linear isoforms) and divergent (detecting circular RNAs) primers were designed through the tool NCBI primer blast. Primer pairs were tested in silico and all those primer pairs that recognized more than one circular isoforms produced by the same host gene were discarded. Primer pairs used in this study are reported in Table S3.

4.6. Northern Blotting

Ten micrograms of whole cell RNA were loaded on a 1.2% denaturing agarose gel. Transfer to the membrane, labelling of probes, hybridization and detection were conducted as previously described [21].

4.7. qRT-PCR

Total RNA was reverse-transcribed in cDNA and amplified using gene-specific primers either in two steps, by MLV-RT (Thermo Fisher Scientific) using random hexamer primers and LightCycler® 480 SYBR Green I Master (Roche Molecular Systems, Inc., Pleasanton, CA, USA), or in one step, using RNA-to-Ct™ 1-Step Kit (Thermo Fisher Scientific). Quantitative PCR was performed on either 7900HT Fast Real-Time PCR System (Thermo Fisher Scientific) or LightCycler 480 (Roche Molecular Systems, Inc.), according to supplier’s protocol. TBP and GAPDH were used as endogenous control genes. The relative amount of gene expression and the fold change for each transcript were calculated according to 2−ΔΔCt method [41].

4.8. Cell Migration Assay

Cell migration was assessed using Oris™ Cell Migration Assay (Platypus Technologies, Madison, WI, USA), as previously described [42]. Briefly, thirty-five thousand cells per well were seeded in the migration ninety-six wells plate with the detection mask attached at the bottom and the stoppers placed. Twenty-four hours after seeding, cells were transfected and grown for other 24 h in a serum free medium. After this time, the stoppers were removed and migration rate was revealed in the detection zone after 24 h by staining the cells with Hoechst® 33342 at a final concentration of 5 µg/mL. Migrated cells were quantified through ImageJ software v. 1.51.

4.9. Cell Viability Assay

Cell viability was measured using MTT (3(4,5-dimethylthiazol-2-yl)2,5-diphenyl-tetrazoliumbromide), as previously described [43]. Briefly, twenty-four thousand cells per well were seeded in a ninety-six wells plate, grown for 24 h and then transfected for other 24 h. Viability was assessed reading absorbance at 580 nm 24 h, 48 h, 72 h and 96 h after transfection.

4.10. In Silico Analysis

RBP binding site predictions within circSMARCA5 sequence were done by providing the circSMARCA5 FASTA sequence as input in RBPmap [44], ATTRACT (https://attract.cnic.es/index) and catRapid (http://service.tartaglialab.com/page/catrapid_group) online tools, respectively. Inverted repeats were retrieved using EMBOSS einverted tool (http://emboss.bioinformatics.nl/cgi-bin/emboss/einverted). UCSC (https://genome.ucsc.edu/) and ClustalW2 (https://www.ebi.ac.uk/Tools/msa/clustalw2/) were used to assess conservation among sequences.

4.11. eCLIP Analysis

Bam-files (accessions ENCFF913MGS, ENCFF374XWQ, and ENCFF833PXB) were downloaded from http://www.encodeproject.org and indexed with samtools. Read densities were visualized using the IGV browser and analyzed using the pysam module in python.

4.12 TCGA Data

TCGA GBM RNA-seq data were retrieved from Cancer RNA-Seq Nexus database (http://syslab4.nchu.edu.tw/).

4.13 Nomenclature of Genes and Proteins

Genes, mRNAs and proteins are named according to the guidelines set by Human Gene Nomenclature [45].

4.14. Statistical Analyses

Unless otherwise specified, statistical significance was assessed by two-tailed student’s t-test.

Supplementary Materials

Supplementary materials can be found at https://www.mdpi.com/1422-0067/19/2/480/s1.

Acknowledgments

This work was supported by a grant from Lega Italiana per la Lotta contro i Tumori (LILT), section of Catania. We wish to thank Massimo Libra (LILT) and Colleagues for his interest in the project. We also acknowledge Fondazione Umberto Veronesi (FUV) for granting Davide Barbagallo a post-doctoral fellowship funding during his stay at Aarhus University. We thank Francesco Certo and Giuseppe Broggi for helping in the management of clinical data; and Lotte Stagsted, Karoline Ebbesen and Katrine Mandrup for their kind inputs. Finally, we thank the Scientific Bureau of the University of Catania for language support.

Author Contributions

Michele Purrello, Thomas Birkballe Hansen, Cinzia Di Pietro, Marco Ragusa and Davide Barbagallo designed and conceived the experiments. Giuseppe Maria Barbagallo and Rosario Caltabiano obtained, characterized and curated clinical data from patients’ biopsies. Davide Barbagallo, Angela Caponnetto and Antonio Morrone performed the experiments. Matilde Cirnigliaro, Duilia Brex, Cristina Barbagallo and Floriana D’Angeli contributed to the acquisition, analysis, and interpretation of data. Michele Purrello and Davide Barbagallo wrote the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Memczak, S.; Jens, M.; Elefsinioti, A.; Torti, F.; Krueger, J.; Rybak, A.; Maier, L.; Mackowiak, S.D.; Gregersen, L.H.; Munschauer, M.; et al. Circular RNAs are a large class of animal RNAs with regulatory potency. Nature 2013, 495, 333–338. [Google Scholar] [CrossRef] [PubMed]
  2. Chen, I.; Chen, C.Y.; Chuang, T.J. Biogenesis, identification, and function of exonic circular RNAs. Wiley Interdiscip. Rev. RNA 2015, 6, 563–579. [Google Scholar] [CrossRef] [PubMed]
  3. Zhang, Y.; Zhang, X.O.; Chen, T.; Xiang, J.F.; Yin, Q.F.; Xing, Y.H.; Zhu, S.; Yang, L.; Chen, L.L. Circular intronic long noncoding RNAs. Mol. Cell 2013, 51, 792–806. [Google Scholar] [CrossRef] [PubMed]
  4. Li, Z.; Huang, C.; Bao, C.; Chen, L.; Lin, M.; Wang, X.; Zhong, G.; Yu, B.; Hu, W.; Dai, L.; et al. Exon-intron circular RNAs regulate transcription in the nucleus. Nat. Struct. Mol. Biol. 2015, 22, 256–264. [Google Scholar] [CrossRef] [PubMed]
  5. Vicens, Q.; Westhof, E. Biogenesis of Circular RNAs. Cell 2014, 159, 13–14. [Google Scholar] [CrossRef] [PubMed]
  6. Conn, S.J.; Pillman, K.A.; Toubia, J.; Conn, V.M.; Salmanidis, M.; Phillips, C.A.; Roslan, S.; Schreiber, A.W.; Gregory, P.A.; Goodall, G.J. The RNA binding protein quaking regulates formation of circRNAs. Cell 2015, 160, 1125–1134. [Google Scholar] [CrossRef] [PubMed]
  7. Suzuki, H.; Tsukahara, T. A view of pre-mRNA splicing from RNase R resistant RNAs. Int. J. Mol. Sci. 2014, 15, 9331–9342. [Google Scholar] [CrossRef] [PubMed]
  8. Salzman, J.; Chen, R.E.; Olsen, M.N.; Wang, P.L.; Brown, P.O. Cell-type specific features of circular RNA expression. PLoS Genet. 2013, 9, e1003777. [Google Scholar] [CrossRef]
  9. Jeck, W.R.; Sorrentino, J.A.; Wang, K.; Slevin, M.K.; Burd, C.E.; Liu, J.; Marzluff, W.F.; Sharpless, N.E. Circular RNAs are abundant, conserved, and associated with ALU repeats. RNA 2013, 19, 141–157. [Google Scholar] [CrossRef] [PubMed]
  10. Lasda, E.; Parker, R. Circular RNAs: Diversity of form and function. RNA 2014, 20, 1829–1842. [Google Scholar] [CrossRef] [PubMed]
  11. Rybak-Wolf, A.; Stottmeister, C.; Glazar, P.; Jens, M.; Pino, N.; Giusti, S.; Hanan, M.; Behm, M.; Bartok, O.; Ashwal-Fluss, R.; et al. Circular RNAs in the Mammalian Brain Are Highly Abundant, Conserved, and Dynamically Expressed. Mol. Cell 2015, 58, 870–885. [Google Scholar] [CrossRef] [PubMed]
  12. Gruner, H.; Cortes-Lopez, M.; Cooper, D.A.; Bauer, M.; Miura, P. CircRNA accumulation in the aging mouse brain. Sci. Rep. 2016, 6, 38907. [Google Scholar] [CrossRef] [PubMed]
  13. Bahn, J.H.; Zhang, Q.; Li, F.; Chan, T.M.; Lin, X.; Kim, Y.; Wong, D.T.; Xiao, X. The landscape of microRNA, Piwi-interacting RNA, and circular RNA in human saliva. Clin. Chem. 2015, 61, 221–230. [Google Scholar] [CrossRef] [PubMed]
  14. Li, Y.; Zheng, Q.; Bao, C.; Li, S.; Guo, W.; Zhao, J.; Chen, D.; Gu, J.; He, X.; Huang, S. Circular RNA is enriched and stable in exosomes: A promising biomarker for cancer diagnosis. Cell Res. 2015, 25, 981–984. [Google Scholar] [CrossRef] [PubMed]
  15. Memczak, S.; Papavasileiou, P.; Peters, O.; Rajewsky, N. Identification and Characterization of Circular RNAs as a New Class of Putative Biomarkers in Human Blood. PLoS ONE 2015, 10, e0141214. [Google Scholar] [CrossRef] [PubMed]
  16. Di Pietro, V.; Ragusa, M.; Davies, D.; Su, Z.; Hazeldine, J.; Lazzarino, G.; Hill, L.J.; Crombie, N.; Foster, M.; Purrello, M.; et al. MicroRNAs as Novel Biomarkers for the Diagnosis and Prognosis of Mild and Severe Traumatic Brain Injury. J. Neurotrauma 2017, 34, 1948–1956. [Google Scholar] [CrossRef] [PubMed]
  17. Romano, G.L.; Platania, C.B.M.; Drago, F.; Salomone, S.; Ragusa, M.; Barbagallo, C.; Di Pietro, C.; Purrello, M.; Reibaldi, M.; Avitabile, T.; et al. Retinal and Circulating miRNAs in Age-Related Macular Degeneration: An In vivo Animal and Human Study. Front. Pharmacol. 2017, 8, 168. [Google Scholar] [CrossRef] [PubMed]
  18. Kristensen, L.S.; Hansen, T.B.; Veno, M.T.; Kjems, J. Circular RNAs in cancer: Opportunities and challenges in the field. Oncogene 2018, 37, 555–565. [Google Scholar] [CrossRef] [PubMed]
  19. Barbagallo, D.; Condorelli, A.; Ragusa, M.; Salito, L.; Sammito, M.; Banelli, B.; Caltabiano, R.; Barbagallo, G.; Zappala, A.; Battaglia, R.; et al. Dysregulated miR-671-5p/CDR1-AS/CDR1/VSNL1 axis is involved in glioblastoma multiforme. Oncotarget 2016, 7, 4746–4759. [Google Scholar] [CrossRef] [PubMed]
  20. Ragusa, M.; Barbagallo, C.; Statello, L.; Condorelli, A.G.; Battaglia, R.; Tamburello, L.; Barbagallo, D.; Di Pietro, C.; Purrello, M. Non-coding landscapes of colorectal cancer. World J. Gastroenterol. 2015, 21, 11709–11739. [Google Scholar] [CrossRef] [PubMed]
  21. Hansen, T.B.; Jensen, T.I.; Clausen, B.H.; Bramsen, J.B.; Finsen, B.; Damgaard, C.K.; Kjems, J. Natural RNA circles function as efficient microRNA sponges. Nature 2013, 495, 384–388. [Google Scholar] [CrossRef] [PubMed]
  22. Zheng, Q.; Bao, C.; Guo, W.; Li, S.; Chen, J.; Chen, B.; Luo, Y.; Lyu, D.; Li, Y.; Shi, G.; et al. Circular RNA profiling reveals an abundant circHIPK3 that regulates cell growth by sponging multiple miRNAs. Nat. Commun. 2016, 7, 11215. [Google Scholar] [CrossRef] [PubMed]
  23. Ashwal-Fluss, R.; Meyer, M.; Pamudurti, N.R.; Ivanov, A.; Bartok, O.; Hanan, M.; Evantal, N.; Memczak, S.; Rajewsky, N.; Kadener, S. circRNA biogenesis competes with pre-mRNA splicing. Mol. Cell 2014, 56, 55–66. [Google Scholar] [CrossRef] [PubMed]
  24. Abdelmohsen, K.; Panda, A.C.; Munk, R.; Grammatikakis, I.; Dudekula, D.B.; De, S.; Kim, J.; Noh, J.H.; Kim, K.M.; Martindale, J.L.; et al. Identification of HuR target circular RNAs uncovers suppression of PABPN1 translation by CircPABPN1. RNA Biol. 2017, 14, 361–369. [Google Scholar] [CrossRef] [PubMed]
  25. Du, W.W.; Yang, W.; Chen, Y.; Wu, Z.K.; Foster, F.S.; Yang, Z.; Li, X.; Yang, B.B. Foxo3 circular RNA promotes cardiac senescence by modulating multiple factors associated with stress and senescence responses. Eur. Heart J. 2017, 38, 1402–1412. [Google Scholar] [CrossRef] [PubMed]
  26. Du, W.W.; Yang, W.; Liu, E.; Yang, Z.; Dhaliwal, P.; Yang, B.B. Foxo3 circular RNA retards cell cycle progression via forming ternary complexes with p21 and CDK2. Nucleic Acids Res. 2016, 44, 2846–2858. [Google Scholar] [CrossRef] [PubMed]
  27. Holdt, L.M.; Stahringer, A.; Sass, K.; Pichler, G.; Kulak, N.A.; Wilfert, W.; Kohlmaier, A.; Herbst, A.; Northoff, B.H.; Nicolaou, A.; et al. Circular non-coding RNA ANRIL modulates ribosomal RNA maturation and atherosclerosis in humans. Nat. Commun. 2016, 7, 12429. [Google Scholar] [CrossRef] [PubMed]
  28. Yang, Z.G.; Awan, F.M.; Du, W.W.; Zeng, Y.; Lyu, J.; Gupta, S.; Yang, W.; Yang, B.B. The Circular RNA Interacts with STAT3, Increasing Its Nuclear Translocation and Wound Repair by Modulating Dnmt3a and miR-17 Function. Mol. Ther. 2017, 25, 2062–2074. [Google Scholar] [CrossRef] [PubMed]
  29. Ghigna, C.; Giordano, S.; Shen, H.; Benvenuto, F.; Castiglioni, F.; Comoglio, P.M.; Green, M.R.; Riva, S.; Biamonti, G. Cell motility is controlled by SF2/ASF through alternative splicing of the Ron protooncogene. Mol. Cell 2005, 20, 881–890. [Google Scholar] [CrossRef] [PubMed]
  30. Song, X.; Zhang, N.; Han, P.; Moon, B.S.; Lai, R.K.; Wang, K.; Lu, W. Circular RNA profile in gliomas revealed by identification tool UROBORUS. Nucleic Acids Res. 2016, 44, e87. [Google Scholar] [CrossRef] [PubMed]
  31. Jumaa, H.; Nielsen, P.J. The splicing factor SRp20 modifies splicing of its own mRNA and ASF/SF2 antagonizes this regulation. EMBO J. 1997, 16, 5077–5085. [Google Scholar] [CrossRef] [PubMed]
  32. Siede, D.; Rapti, K.; Gorska, A.A.; Katus, H.A.; Altmuller, J.; Boeckel, J.N.; Meder, B.; Maack, C.; Volkers, M.; Muller, O.J.; et al. Identification of circular RNAs with host gene-independent expression in human model systems for cardiac differentiation and disease. J. Mol. Cell. Cardiol. 2017, 109, 48–56. [Google Scholar] [CrossRef] [PubMed]
  33. Guo, J.; Jia, J.; Jia, R. PTBP1 and PTBP2 impaired autoregulation of SRSF3 in cancer cells. Sci. Rep. 2015, 5, 14548. [Google Scholar] [CrossRef] [PubMed]
  34. Lareau, L.F.; Inada, M.; Green, R.E.; Wengrod, J.C.; Brenner, S.E. Unproductive splicing of SR genes associated with highly conserved and ultraconserved DNA elements. Nature 2007, 446, 926–929. [Google Scholar] [CrossRef] [PubMed]
  35. Ni, J.Z.; Grate, L.; Donohue, J.P.; Preston, C.; Nobida, N.; O’Brien, G.; Shiue, L.; Clark, T.A.; Blume, J.E.; Ares, M., Jr. Ultraconserved elements are associated with homeostatic control of splicing regulators by alternative splicing and nonsense-mediated decay. Genes Dev. 2007, 21, 708–718. [Google Scholar] [CrossRef] [PubMed]
  36. Das, S.; Krainer, A.R. Emerging functions of SRSF1, splicing factor and oncoprotein, in RNA metabolism and cancer. Mol. Cancer Res. 2014, 12, 1195–1204. [Google Scholar] [CrossRef] [PubMed]
  37. Dvinge, H.; Kim, E.; Abdel-Wahab, O.; Bradley, R.K. RNA splicing factors as oncoproteins and tumour suppressors. Nat. Rev. Cancer 2016, 16, 413–430. [Google Scholar] [CrossRef] [PubMed]
  38. Cheung, H.C.; Hai, T.; Zhu, W.; Baggerly, K.A.; Tsavachidis, S.; Krahe, R.; Cote, G.J. Splicing factors PTBP1 and PTBP2 promote proliferation and migration of glioma cell lines. Brain 2009, 132 Pt 8, 2277–2288. [Google Scholar] [CrossRef] [PubMed]
  39. Louis, D.N.; Ohgaki, H.; Wiestler, O.D.; Cavenee, W.K.; Burger, P.C.; Jouvet, A.; Scheithauer, B.W.; Kleihues, P. The 2007 WHO classification of tumours of the central nervous system. Acta Neuropathol. 2007, 114, 97–109. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Caltabiano, R.; Puzzo, L.; Barresi, V.; Ieni, A.; Loreto, C.; Musumeci, G.; Castrogiovanni, P.; Ragusa, M.; Foti, P.; Russo, A.; et al. ADAM 10 expression in primary uveal melanoma as prognostic factor for risk of metastasis. Pathol. Res. Pract. 2016, 212, 980–987. [Google Scholar] [CrossRef] [PubMed]
  41. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef] [PubMed]
  42. Ragusa, M.; Statello, L.; Maugeri, M.; Majorana, A.; Barbagallo, D.; Salito, L.; Sammito, M.; Santonocito, M.; Angelica, R.; Cavallaro, A.; et al. Specific alterations of the microRNA transcriptome and global network structure in colorectal cancer after treatment with MAPK/ERK inhibitors. J. Mol. Med. (Berl.) 2012, 90, 1421–1438. [Google Scholar] [CrossRef] [PubMed]
  43. Barbagallo, D.; Piro, S.; Condorelli, A.G.; Mascali, L.G.; Urbano, F.; Parrinello, N.; Monello, A.; Statello, L.; Ragusa, M.; Rabuazzo, A.M.; et al. miR-296-3p, miR-298-5p and their downstream networks are causally involved in the higher resistance of mammalian pancreatic alpha cells to cytokine-induced apoptosis as compared to beta cells. BMC Genom. 2013, 14, 62. [Google Scholar] [CrossRef] [PubMed]
  44. Paz, I.; Kosti, I.; Ares, M., Jr.; Cline, M.; Mandel-Gutfreund, Y. RBPmap: A web server for mapping binding sites of RNA-binding proteins. Nucl. Acids Res. 2014, 42, W361–W367. [Google Scholar] [CrossRef] [PubMed]
  45. Wain, H.M.; Bruford, E.A.; Lovering, R.C.; Lush, M.J.; Wright, M.W.; Povey, S. Guidelines for human gene nomenclature. Genomics 2002, 79, 464–470. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Expression of candidate circRNAs in GBM and control (NORM) samples. (A) Heat map representing the expression profile of twelve candidate circRNAs in GBM and control samples. Expression values (reported as DCt) higher and lower than those of the control are shown in yellow and blue, respectively. circSHPRH and circSMARCA5 are in bold, due to their significant downregulation (p-value circSHPRH = 0.015, p-value circSMARCA5 = 0.017, n = 5, Student’s t-test). (B) Box-plots representing the expression of circSMARCA5 and linearSMARCA5 in GBM and control samples. Expression values are reported as −DCt relative to (mRNA) TBP (p-values are shown in the figure, nGBM = 56, nNORM = 7, Student’s t-test).
Figure 1. Expression of candidate circRNAs in GBM and control (NORM) samples. (A) Heat map representing the expression profile of twelve candidate circRNAs in GBM and control samples. Expression values (reported as DCt) higher and lower than those of the control are shown in yellow and blue, respectively. circSHPRH and circSMARCA5 are in bold, due to their significant downregulation (p-value circSHPRH = 0.015, p-value circSMARCA5 = 0.017, n = 5, Student’s t-test). (B) Box-plots representing the expression of circSMARCA5 and linearSMARCA5 in GBM and control samples. Expression values are reported as −DCt relative to (mRNA) TBP (p-values are shown in the figure, nGBM = 56, nNORM = 7, Student’s t-test).
Ijms 19 00480 g001
Figure 2. Correlation between circSMARCA5 expression and glioma grade malignancy. Box-plot representing circSMARCA5 expression in three grade tumor groups (GBM, Grade III—GIII- and Grade II—GII-gliomas) and two different control groups (Cortex and Cerebellum). Expression values are reported as reads per million mapped reads (RPM) (* p-value < 0.05, ANOVA test).
Figure 2. Correlation between circSMARCA5 expression and glioma grade malignancy. Box-plot representing circSMARCA5 expression in three grade tumor groups (GBM, Grade III—GIII- and Grade II—GII-gliomas) and two different control groups (Cortex and Cerebellum). Expression values are reported as reads per million mapped reads (RPM) (* p-value < 0.05, ANOVA test).
Ijms 19 00480 g002
Figure 3. Effects of circSMARCA5 overexpression on migration and viability of U87MG cells. (A) Migration of untransfected (UC), pcDNA3 and pcDNA3_circSMARCA5 (circSMARCA5) transfected U87MG cells. Cells were stained with Hoechst® 33342 at a final concentration of 5 µg/mL (upper panel). Quantitative data are reported as the number of pixel within the detection area (lower panel, * p-value < 0.05, n = 3, Student’s t-test). (B) Viability of untransfected (UC), pcDNA3 and pcDNA3_circSMARCA5 (circSMARCA5) transfected U87MG cells. Data are reported as baseline-corrected absorbance at 580 nm (n = 6).
Figure 3. Effects of circSMARCA5 overexpression on migration and viability of U87MG cells. (A) Migration of untransfected (UC), pcDNA3 and pcDNA3_circSMARCA5 (circSMARCA5) transfected U87MG cells. Cells were stained with Hoechst® 33342 at a final concentration of 5 µg/mL (upper panel). Quantitative data are reported as the number of pixel within the detection area (lower panel, * p-value < 0.05, n = 3, Student’s t-test). (B) Viability of untransfected (UC), pcDNA3 and pcDNA3_circSMARCA5 (circSMARCA5) transfected U87MG cells. Data are reported as baseline-corrected absorbance at 580 nm (n = 6).
Ijms 19 00480 g003
Figure 4. SRFS1 eCLIP data from ENCODE on K562 cells. (A) eCLIP read density within the SMARCA5 gene locus (chr4:144461450-144466050, hg19), containing the circSMARCA5 and the immediate flanking exons, visualized by the Integrative Genomics Viewer (IGV). (B) Max read density obtained on circSMARCA5 and in the immediate flanking exons (as shown in a) from SRSF1 eCLIP replicates 1 and 2.
Figure 4. SRFS1 eCLIP data from ENCODE on K562 cells. (A) eCLIP read density within the SMARCA5 gene locus (chr4:144461450-144466050, hg19), containing the circSMARCA5 and the immediate flanking exons, visualized by the Integrative Genomics Viewer (IGV). (B) Max read density obtained on circSMARCA5 and in the immediate flanking exons (as shown in a) from SRSF1 eCLIP replicates 1 and 2.
Ijms 19 00480 g004
Figure 5. circSMARCA5 stimulates inclusion of exon 4 within SRSF3 mRNA. Graphs representing the expression of the mRNA of the splicing factor SRFS3 in: (A) U87MG overexpressing circSMARCA5; and (B,C) GBM biopsies. (A) U87MG overexpressing circSMARCA5 show a significant increased expression of SRSF3 containing exon 4 (SRSF3 Ex4) (** p-value < 0.01, n = 3, Student’s t-test). Data are reported as fold-change (FC) versus untransfected cells (UC). (B) Both SRSF3 No Ex4 and SRSF3 Ex4 are significantly overexpressed in GBM biopsies (GBM) with respect to normal brain parenchyma (NORM) (* p-value < 0.05; ** p-value < 0.01, nGBM = 56, nNORM = 7, Student’s t-test). Expression values are reported as −DCt relative to (mRNA) TBP.(C) Correlation between SRSF3 Ex4/SRSF3 No Ex4 ratio and circSMARCA5 (r-value = 0.36, p-value = 0.004, Spearman’s correlation test). Data are represented as DCt relative to (mRNA) TBP in a scatter plot.
Figure 5. circSMARCA5 stimulates inclusion of exon 4 within SRSF3 mRNA. Graphs representing the expression of the mRNA of the splicing factor SRFS3 in: (A) U87MG overexpressing circSMARCA5; and (B,C) GBM biopsies. (A) U87MG overexpressing circSMARCA5 show a significant increased expression of SRSF3 containing exon 4 (SRSF3 Ex4) (** p-value < 0.01, n = 3, Student’s t-test). Data are reported as fold-change (FC) versus untransfected cells (UC). (B) Both SRSF3 No Ex4 and SRSF3 Ex4 are significantly overexpressed in GBM biopsies (GBM) with respect to normal brain parenchyma (NORM) (* p-value < 0.05; ** p-value < 0.01, nGBM = 56, nNORM = 7, Student’s t-test). Expression values are reported as −DCt relative to (mRNA) TBP.(C) Correlation between SRSF3 Ex4/SRSF3 No Ex4 ratio and circSMARCA5 (r-value = 0.36, p-value = 0.004, Spearman’s correlation test). Data are represented as DCt relative to (mRNA) TBP in a scatter plot.
Ijms 19 00480 g005aIjms 19 00480 g005b
Table 1. Clinical data of GBM and control samples.
Table 1. Clinical data of GBM and control samples.
Type of SamplesN° of SamplesMean Age (Years ± StdDev)SexMean Overall Survival (Months)Preoperative Tumor Volume (mm3)
MF
Training set (Fresh-frozen biopsies)1060.5 ± 12.15519.1 ± 7.239.9 ± 11.8
Test set (FFPE biopsies)5662 ± 12.7332317 ± 14.228.2 ± 9.2
Normal Brain Parenchyma760.4 ± 11.125
FirstChoice® Human Brain Reference Total RNA1 (commercially available)68.3 ± 151310
Table 2. Candidate circRNAs.
Table 2. Candidate circRNAs.
#Candidate circRNA (circBase ID)Parental Gene SymbolKnown Modulation or Function of circRNA Expression in Specific Cell ContextSource (PMID)
1hsa_circ_0028270ATXN2Upregulated during EMT. This circRNA is also highly expressed in several SNC districts (see PMID: 25921068)25768908
2hsa_circ_0008702GNB1Downregulated during EMT. This circRNA is also highly expressed in several SNC districts (see PMID: 25921068)25768908
3hsa_circ_0000284HIPK3Involved in cell growth. Highly expressed in normal Brain27050392
4hsa_circ_0008002POLE2Upregulated during EMT. This circRNA is also highly expressed in several SNC districts (see PMID: 25921068)25768908
5hsa_circ_0132250RIMS1Downregulated in GBM vs. Normal Brain (another dataset)26873924
6hsa_circ_0099634RmstLncRNA Rmst (host gene of the same name circ_Rmst) regulates neuronal differentiation in mouse25921068
7hsa_circ_0054598RTN4Upregulated during neuronal differentiation both in humans and mice25921068
8hsa_circ_0001649SHPRHUpregulated during EMT. This circRNA is also highly expressed in several SNC districts (see PMID: 25921068)25768908
9hsa_circ_0003694SMAD2Upregulated during EMT. This circRNA is also highly expressed in several SNC districts (see PMID: 25921068)25768908
10hsa_circ_0001445SMARCA5Upregulated during EMT. This circRNA is also highly expressed in several SNC districts (see PMID: 25921068)25768908
11hsa_circ_0073237VCANUpregulated in GBM vs. Normal Brain (another dataset)26873924
12hsa_circ_0004383ZNF292Upregulated in HUVEC under hypoxia. Its silencing reduces endothelial cell proliferation and suppresses tube formation by inhibiting glioma cell proliferation and cell cycle progression in human glioma U87MG and U251 cells26377962; 27613831

Share and Cite

MDPI and ACS Style

Barbagallo, D.; Caponnetto, A.; Cirnigliaro, M.; Brex, D.; Barbagallo, C.; D’Angeli, F.; Morrone, A.; Caltabiano, R.; Barbagallo, G.M.; Ragusa, M.; et al. CircSMARCA5 Inhibits Migration of Glioblastoma Multiforme Cells by Regulating a Molecular Axis Involving Splicing Factors SRSF1/SRSF3/PTB. Int. J. Mol. Sci. 2018, 19, 480. https://doi.org/10.3390/ijms19020480

AMA Style

Barbagallo D, Caponnetto A, Cirnigliaro M, Brex D, Barbagallo C, D’Angeli F, Morrone A, Caltabiano R, Barbagallo GM, Ragusa M, et al. CircSMARCA5 Inhibits Migration of Glioblastoma Multiforme Cells by Regulating a Molecular Axis Involving Splicing Factors SRSF1/SRSF3/PTB. International Journal of Molecular Sciences. 2018; 19(2):480. https://doi.org/10.3390/ijms19020480

Chicago/Turabian Style

Barbagallo, Davide, Angela Caponnetto, Matilde Cirnigliaro, Duilia Brex, Cristina Barbagallo, Floriana D’Angeli, Antonio Morrone, Rosario Caltabiano, Giuseppe Maria Barbagallo, Marco Ragusa, and et al. 2018. "CircSMARCA5 Inhibits Migration of Glioblastoma Multiforme Cells by Regulating a Molecular Axis Involving Splicing Factors SRSF1/SRSF3/PTB" International Journal of Molecular Sciences 19, no. 2: 480. https://doi.org/10.3390/ijms19020480

APA Style

Barbagallo, D., Caponnetto, A., Cirnigliaro, M., Brex, D., Barbagallo, C., D’Angeli, F., Morrone, A., Caltabiano, R., Barbagallo, G. M., Ragusa, M., Di Pietro, C., Hansen, T. B., & Purrello, M. (2018). CircSMARCA5 Inhibits Migration of Glioblastoma Multiforme Cells by Regulating a Molecular Axis Involving Splicing Factors SRSF1/SRSF3/PTB. International Journal of Molecular Sciences, 19(2), 480. https://doi.org/10.3390/ijms19020480

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop