Next Article in Journal
A Potential Nutraceutical Candidate Lactucin Inhibits Adipogenesis through Downregulation of JAK2/STAT3 Signaling Pathway-Mediated Mitotic Clonal Expansion
Next Article in Special Issue
Alternative Splicing: Recent Insights into Mechanisms and Functional Roles
Previous Article in Journal
Mouse Embryonic Stem Cell-Derived Ureteric Bud Progenitors Induce Nephrogenesis
Previous Article in Special Issue
Splicing Players Are Differently Expressed in Sporadic Amyotrophic Lateral Sclerosis Molecular Clusters and Brain Regions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Poison-Exon Inclusion in DHX9 Reduces Its Expression and Sensitizes Ewing Sarcoma Cells to Chemotherapeutic Treatment

by
Ramona Palombo
1,
Veronica Verdile
1,2 and
Maria Paola Paronetto
1,2,*
1
Laboratory of Cellular and Molecular Neurobiology, IRCCS Fondazione Santa Lucia, 00143 Rome, Italy
2
Department of Movement, Human and Health Sciences, Università degli Studi di Roma “Foro Italico”, Piazza Lauro de Bosis, 15, 00135 Rome, Italy
*
Author to whom correspondence should be addressed.
Cells 2020, 9(2), 328; https://doi.org/10.3390/cells9020328
Submission received: 27 December 2019 / Revised: 27 January 2020 / Accepted: 29 January 2020 / Published: 31 January 2020

Abstract

:
Alternative splicing is a combinatorial mechanism by which exons are joined to produce multiple mRNA variants, thus expanding the coding potential and plasticity of eukaryotic genomes. Defects in alternative splicing regulation are associated with several human diseases, including cancer. Ewing sarcoma is an aggressive tumor of bone and soft tissue, mainly affecting adolescents and young adults. DHX9 is a key player in Ewing sarcoma malignancy, and its expression correlates with worse prognosis in patients. In this study, by screening a library of siRNAs, we have identified splicing factors that regulate the alternative inclusion of a poison exon in DHX9 mRNA, leading to its downregulation. In particular, we found that hnRNPM and SRSF3 bind in vivo to this poison exon and suppress its inclusion. Notably, DHX9 expression correlates with that of SRSF3 and hnRNPM in Ewing sarcoma patients. Furthermore, downregulation of SRSF3 or hnRNPM inhibited DHX9 expression and Ewing sarcoma cell proliferation, while sensitizing cells to chemotherapeutic treatment. Hence, our study suggests that inhibition of hnRNPM and SRSF3 expression or activity could be exploited as a therapeutic tool to enhance the efficacy of chemotherapy in Ewing sarcoma.

1. Introduction

Alternative splicing is a sophisticated mechanism by which exons are joined in different combinations to generate multiple mRNA variants, thus fine-tuning gene expression programs [1,2]. Regulation of alternative splicing enacts cell growth, differentiation, and survival [3] and is deeply involved in the development and progression of human pathological conditions [4].
A surveillance mechanism of regulation that limits the expression of aberrant transcripts in the cells is nonsense-mediated mRNA decay (NMD) [5], which is triggered by the presence of a premature stop codon (PTC) in non-last exons. Such “poison” cassette exons are normally skipped, while their inclusion in the mature mRNA targets it to NMD [6,7,8]. Indeed, the presence of a PTC located more than 50 nucleotides upstream of the final exon–exon junction is thought to signal the premature or aberrant nature of a transcript [6,7,8]. In addition to poison cassette exons, other classes of alternative splicing events, such as alternative 5′- and 3′-splice sites or intron retention, also contribute to the alternative splicing-NMD post-transcriptional regulation [9], thus limiting the accumulation of unproductive splicing variants. Noteworthy, the expression of many RNA binding proteins (RBPs) involved in splicing regulation, such as the serine/arginine-rich (SR) proteins, is fine-tuned by NMD [6,8].
Changes in the RNA polymerase II (RNAPII) elongation rate contribute to the regulation of NMD-linked alternative splicing, thus coordinating the cellular requirements of splicing factors and RBPs in response to internal and external cues. RNAPII slowing down and occupancy over the intronic regions flanking the regulated exons favor the recognition of suboptimal splice sites, such as those in poison exons [10]. We previously identified a poison exon (exon 6A) in the DHX9 gene, whose inclusion targets the transcript to NMD [11]. Inclusion of exon 6A is normally repressed, thus insuring high expression levels of DHX9. However, reduction in the RNAPII elongation rate within the DHX9 transcription unit favors exon 6A inclusion and targets the transcript to NMD [11]. Both UV light irradiation and etoposide treatment induced this event by slowing down the RNAPII [11], with the consequent decrease in DHX9 expression, thus leading to higher sensitivity of Ewing sarcoma cells to genotoxic stress [11,12]. Nevertheless, the mechanism by which exon 6A inclusion is normally repressed in Ewing sarcoma cells is currently unknown.
DHX9 is a member of the DExH subgroup of RNA helicases, which play important roles in several aspects of RNA metabolism [12]. DHX9 is involved in the regulation of gene expression by acting as a scaffold for the interaction of breast cancer 1 (BRCA1) [13] and cyclic adenosine monophosphate (AMP) response element-binding protein-binding protein (CBP) [14] with the RNAPII holoenzyme, thus modulating their activity and regulating transcription. Moreover, DHX9 is involved in the maintenance of genomic stability [15,16,17]. In Ewing sarcoma, DHX9 forms a complex with the EWS-FLI1 oncoprotein and modulates EWS-FLI1-dependent transcription [18]. In particular, the functional interaction between DHX9 and EWS-FLI1 enhances the engagement of the transcriptional machinery at responsive promoters, induces local changes in chromatin structure, and unwinds the DNA. DHX9 also interacts with the RBP Sam68 and with the promoter-associated noncoding RNA pncCCND1b to form an RNA-protein complex inhibiting CCND1 transcription in Ewing sarcoma cells [19].
The EWS-FLI1/DHX9 complex represents a good therapeutic target for Ewing sarcoma [11,18,20,21,22,23]. Thus, understanding the regulation of the DHX9 poison-exon 6A inclusion might pave the way for novel splicing-directed strategies to inhibit DHX9 gene expression and EWS-FLI1 oncogenic activity. Herein, we screened a library of siRNAs targeting RBPs to identify factors that regulate DHX9 alternative splicing. We identified hnRNPM and SRSF3 as key factors required to suppress exon 6A inclusion and maintain high DHX9 expression in Ewing sarcoma cells. Importantly, downregulation of SRSF3 or hnRNPM sensitized Ewing sarcoma cells to doxorubicin, a genotoxic agent used in Ewing sarcoma chemotherapy. Therefore, our study suggests that inhibition of hnRNPM or SRSF3 expression could be exploited as a therapeutic tool in Ewing sarcoma.

2. Materials and Methods

2.1. Cell Cultures and Drug Treatment

Ewing sarcoma cell lines TC-71 (RRID: CVCL_2213 and SK-N-MC RRID: CVCL_0530) were purchased from DSMZ (Braunschweig, Germany). LAP-35 (RRID: CVCL_A096) was a generous gift from Drs. Katia Scotlandi and Cristina Manara. The absence of mycoplasma contamination was verified every two months by PCR analysis. Cells were maintained in culture in Iscove’s modified Dulbecco’s medium (IMDM) (GIBCO—Thermo Fisher Scientific, Waltham, USA, Massachusetts), supplemented with 10% fetal bovine serum, and penicillin and streptomycin (GIBCO) and maintained at 37 °C in humidified 5% CO2 atmosphere. For doxorubicin treatment, Ewing sarcoma cells were treated for the indicated time with either DMSO or the indicated concentrations of doxorubicin (ranging from 0.1 nM to 150 nM).

2.2. Transfections

Lipofectamine RNAiMax reagent (Thermo Fisher Scientific, Waltham, MA, USA) was used for siRNA transfections. Briefly, 20,000 TC-71 cells were subjected to double pulse of reverse-transfection by using 2 μL of Lipofectamine RNAiMAX, and cells were collected or re-plated for further experiments 24 h after the last pulse of transfection. siRNAs and primers oligonucleotides were purchased from Sigma–Aldrich (Milan, Italy). Sequences are listed in Supplementary Tables S1 and S2, respectively.

2.3. SDS–PAGE and Western Blot Analyses

For protein extract preparation, cells were washed twice with ice-cold phosphate-buffered saline (PBS), resuspended in RIPA lysis buffer (150 mM NaCl, 50 mM Tris-HCl pH 7.5, 2 mM EDTA, 0.1 % in sodium dodecyl sulfate (SDS), 0.5% sodium deoxycolate,1mM dithiothreitol, 0.5 mM Na-orthovanadate, 1%, 10 mM β-glycerolphosphate, 10 mM sodium fluoride, 1% NP-40 and Protease-Inhibitor Cocktail (Sigma–Aldrich)) and kept on ice for 10 min. Soluble protein extracts were separated by centrifugation at 12,000 rpm for 10 min and diluted in Laemlli sample buffer. The obtained cell lysates were resolved on SDS–polyacrylamide gels (SDS-PAGE) and transferred on PVDF membrane Hybond TM-P (Amersham Bioscience, Buckinghamshire, UK). Membranes were saturated with 5% BSA at room temperature and incubated with the following primary antibodies at 4 °C overnight: mouse GAPDH (SC-32233), rabbit DHX9 (SC-66997), mouse SRSF1 (SC-33652), mouse HNRNPM (SC-20002), rabbit HNRNPK (SC-25373), mouse FUS (SC-47711), and β-actin (SC-47778) from Santa Cruz Biotechnology Inc. (Dallas, TX, USA), and mouse SRSF3 (Abnova, Taipei, Taiwan, (H00006428-MO8)). Secondary anti-mouse or anti-rabbit IgGs conjugated to horseradish peroxidase (Amersham Bioscience) were incubated with the membranes for 1 h at room temperature at a 1:10,000 dilution. Immunostained bands were detected by a chemiluminescent method (Thermo Fisher Scientific).

2.4. Real-time Quantitative PCR Analyses (RT-qPCR)

RNA was isolated and DNase digested using either Trizol (Thermo Fisher Scientific) or RNeasy kit (Qiagen, Hilden, Germany). Total RNA (1 µg) was reverse transcribed by using M-MLV Reverse Transcriptase (Promega, Madison, WI, USA) following the manufacturer’s instructions. RT reaction was used as a template together with the different primers listed in Supplementary Table S2. Primers were designed using Primer 3 Plus (http://www.bioinformatics.nl/cgi-bin/primer3plus/primer3plus.cgi) and Primer-Blast using the reference and the alternative RefSeq accession numbers. For RT-qPCR, the primers were designed such that their annealing temperature was 60 °C, generating single-amplification products in the range of 60- to 120-base-pairs (bp) long. PCR amplification was carried out with 1 μL of the 1:10 diluted reverse transcription sample with 10 μL of 2X SYBR Green Master Mix (Roche, Basel, Germany) and 4 pmol of specific gene primer pairs in a 20 μL total volume in 96-well microtiter plates. PCR reactions were run in triplicates on a LightCycler 480 system (Roche). Each experiment was performed at least in triplicate; data are represented as the mean ± standard deviation (SD). For all experiments, no-RT controls have been performed.

2.5. CLIP Assays

Cross-linked and immunoprecipitation (CLIP) assays were performed as previously described [24]. In brief, TC-71 cells were irradiated once with 400 mJ/cm2 in a Stratalinker 2400 at 254 nm. Cell suspension was centrifuged at 4000 rpm for 3 min, and the pellet was incubated for 10 min on ice in lysis buffer (100 mM NaCl, 10 mM MgCl2, 30 mM Hepes pH 7.6, 2 mM EDTA pH 8, 10% Glycerol, 0.5% Tryton-X100, RNase inhibitor, cocktail protease inhibitor (Sigma–Aldrich), 1 mM dithiothreitol, 0.5 mM Na-ortovanadate, 1%, 10 mM β-glycerolphosphate, and RNase inhibitor (Promega). Samples were briefly sonicated and incubated with 10 μL of 1/1000 RNase I (Ambion, AM2295, Thermo Fisher Scientific) dilution and 2 μL Turbo DNase (Ambion, AM2238, Thermo Fisher Scientific) for 3 min at 37 °C shaking at 1100 rpm, and then centrifuged at 13,000 rpm for 10 min at 4 °C. One point five milligrams of the extract was immunoprecipitated overnight using anti-hnRNPM, anti-SRSF1, anti-SRSF3, anti-hnRNPK, anti-FUS antibodies or purified rabbit or mouse IgGs (negative control) in the presence of protein A/G magnetic Dynabeads (Life Technologies–Thermo Fisher Scientific). Immunoprecipitates were incubated overnight at 4 °C under constant rotation. After stringent washes with high salt buffer (300 mM NaCl, 10 mM MgCl2, 30 mM Hepes pH 7.6, 2 mM EDTA pH 8, 10% Glycerol, 0.5% Tryton-X100) beads were equilibrated with PK buffer (100 mM Tris-HCl, pH 7.4, 50 mM NaCl, 10 mM EDTA). An aliquot (10%) was kept as a control of immunoprecipitation, while the rest was treated with 50 µg Proteinase K and incubated for 20 min at 37 °C shaking at 1100 rpm. Seven molars of urea was added to the PK buffer, and incubation was performed for a further 20 min at 37 °C and 1100 rpm. The solution was collected, and phenol/CHCl3 (Ambion, 9722—Thermo Fisher Scientific) was added. After incubation for 5 min at 30 °C shaken at 1100 rpm, phases were separated by spinning for 5 min at 13,000 rpm at room temperature. The aqueous layer was transferred into a new tube and precipitated by the addition of 0.5 μL glycoblue (Ambion, 9510—Thermo Fisher Scientific), 3 M sodium acetate pH 5.5, and 100% ethanol. After mixing, the solution containing retained RNA was precipitated overnight at −20 °C. Purified RNA was used for qPCR analysis. Primers used are listed in Supplementary Table S2.

2.6. MTS Proliferation Assay

Cell proliferation was determined using the Cell Titer A96 3-(4,5-dimethylthiazol-2-yl)-5-(3-carboxymethoxyphenyl)-2-(4-sulfophenyl)-2H-tetrazolium, inner salt (MTS) method according to the manufacturer’s instructions (Promega) by plating 5,000 cells/well in 96-well culture plates.

2.7. Colony Formation Assay

Cells were plated in 35 mm plates at a density of 3,000 cells/plate. After one day, cells were transfected with either scrambled or siRNAs oligonucleotides. Twenty-four hours after transfection, cells were treated with the indicated reagents and incubated at 37 °C in a humidified atmosphere containing 5% CO2 for 12 days, replacing medium every two days. At the end of the incubation period, cells were washed with PBS, fixed in methanol for 10 min at RT, and stained 30 min at room temperature with 0.05% Crystal Violet in distilled water on a rotating shaker. After staining, cells were washed twice with tap water and air-dried overnight. Clones were counted, and the percentage of control was calculated.

2.8. Patient Datasets

The patient dataset used for the analyses in Figure 1 and Figure 2 derives from public available Affymetrix Human Genome U133 Plus 2.0 microarray data of a total of 88 Ewing sarcoma samples (gse17679; Ethical Review Board of Helsinki University Central Hospital no. 329 HUS/E0/05 and 328 HUS/W0/05) [25]. Data are available online at the following link [26]. Detailed patient information is available at the following link: https://www.ncbi.nlm.nih.gov/geo/geo2r/?acc=GSE17679.

2.9. Bioinformatic Analysis

Event-free and overall survival analyses were performed on a dataset composed of 64 Ewing sarcoma, 4 Askin, and 20 PNET (Peripheral Primitive Neuro-Ectodermal Tumors) patients (dataset: GSE17679 [25]) using R2: Genomics Analysis and Visualization Platform website (http://r2.amc.nl).

2.10. Statistical Analysis

Statistical analysis for biological assays was performed by using Graphpad Prism software to calculate EC50, 2-tailed unpaired t-test or one-way or two-way ANOVA, as appropriate. All data are presented as mean ± S.D.

3. Results

3.1. DHX9 Expression Correlates with Worse Prognosis in Ewing Sarcoma Patients

Given the functional relevance of DHX9 in genomic instability and in cancer [12], we asked whether its expression represents a prognostic factor in Ewing sarcoma. Analysis of event-free survival and overall survival in a dataset comprising 64 Ewing sarcoma, 4 Askin, and 20 perypheral primitive neuroectodermal tumors [25], which all belong to the family of Ewing tumors, showed that high DHX9 expression correlates with worse prognosis of the patients (Figure 1A,B). No significant differences in DHX9 expression were observed between Askin, Ewing, and PNET tumors (Figure 1C). However, we observed a significantly higher expression of DHX9 in metastatic tumors versus primary tumors (p = 6.2 × 103). These results support a role of DHX9 in Ewing sarcoma malignancy and highlight its potential contribution to the metastatization process (Figure 1D).

3.2. A siRNA Library Identifies Regulators of DHX9 Alternative Splicing

DHX9 pre-mRNA can undergo alternative splicing to produce either a transcript translated into the full-length protein (NM001357) or a noncoding transcript (NR033302), containing the poison exon 6A, which is targeted to NMD (Figure 3A; [11]). To identify endogenous regulators of this alternative exon, we screened a siRNA library for splicing factors belonging to the serine–arginine (SR) rich and the heterogeneous nuclear ribonucleoprotein (hnRNP) families, which represent the main regulators of alternative splicing [27,28]. Oligonucleotides to knockdown either SR proteins or hnRNPs were transfected in TC-71 Ewing sarcoma cells, and RNA was extracted 48 h later. RT-qPCR analysis was performed to verify the downregulation of the RBP transcripts (Figure 3B,C) and to assess the level of exon 6A inclusion in cells transfected with the siRNA library. The level of exon 6A inclusion was normalized for both a constitutive exon (DHX9 exon 4; Figure 3D,E) and the exon junction exon 6-exon 7 (Supplementary Figure S1A,B). Increased inclusion of exon 6A was observed upon knockdown of SRSF1, SRSF3, SRSF10, HNRNPM, and HNRNPP (FUS), with respect to cells transfected with a control siRNA. On the other hand, knockdown of the HNRNPK transcript promoted further skipping of the alternative exon. Downregulation of the most significant regulators (SRSF1, SRSF3, hnRNPM, hnRNPK, and FUS) was then confirmed at the protein level by Western blot analysis (Supplementary Figure S1C). Given the high similarity between FUS and EWS, we also ascertained that FUS downregulation did not affect neither EWS nor EWS-FLI1 expression (Supplementary Figure S1D). Collectively, these results suggest that several specific SR proteins and hnRNPs contribute to the regulation of DHX9 alternative splicing in TC-71 cells.

3.3. CLIP Assay Unveils Direct Binding of Specific SR Proteins and HnRNPs to DHX9 Pre-mRNA

SR proteins and hnRNPs generally regulate splice site selection by direct binding to the pre-mRNA [2,27]. In silico analysis was performed by querying the Splice Aid database (http://193.206.120.249/splicing_tissue.html) to identify potential binding sites for the main regulators of exon 6A splicing (SRSF1, SRSF3, SRSF10, hnRNPK, hnRNPM, and FUS) in the region of the pre-mRNA encompassing the regulated exon. Remarkably, we found that DHX9 exon 6A contains consensus motifs for all of them, except SRSF10 (Figure 4A). To verify the binding of these factors in live cells, we performed cross-linked and immunoprecipitation (CLIP) experiments of SRSF1, SRSF3, hnRNPM, hnRNPK, and FUS from UV-cross-linked TC-71 cell extracts. RT-qPCR analysis of the RNA associated with these splicing factors revealed that all of them were able to bind in vivo to the DHX9 pre-mRNA within the exon 6A region, although with different affinities (Figure 4B). On the other hand, as predicted from the in silico analysis (Figure 4A), only hnRNPM was also able to bind in the upstream intron 6 region (Figure 4B). These experiments indicate that SRSF1, SRSF3, hnRNPM, hnRNPK, and FUS are endogenous regulators of DHX9 exon 6A alternative splicing in TC-71 cells.

3.4. Expression of DHX9, HNRNPM, and SRSF3 is Positively Correlated in Ewing Sarcoma Patients

To test whether the expression of splicing factors involved in exon 6A regulation was correlated with DHX9 mRNA levels in Ewing sarcoma patients, we performed Pearson correlation analyses using a dataset of 64 Ewing sarcoma patients, 4 Askin, and 20 PNET tumors. A significant correlation was observed with all the splicing regulators of exon 6A, with highest R values observed with SRSF3 (R = 0.759, p = 3.46 × 10−23 and HNRNPM (R = 0.755, p = 7.59 × 10−23) transcripts. Remarkably, RBPs that did not regulate DHX9 splicing in our screening, such as hnRNPL and SRSF5 (Figure 4), did not correlate with DHX9 expression (Supplementary Figure S2). Since high DHX9 expression correlates with poor prognosis (Figure 1), these results suggest that regulation of DHX9 expression by SRSF3 and hnRNPM could represent a prognostic factor in Ewing sarcoma pathogenesis.

3.5. Depletion of SRSF3 and HnRNPM Affects the Expression of EWS-FLI1 Target Genes

Since the inclusion of exon 6A targets the DHX9 transcript to NMD, we performed Western blot analysis to evaluate whether SRSF3 and hnRNPM knockdown affected DHX9 expression. A significant decrease (35% for siHNRNPM and 50% for siSRSF3) in the DHX9 protein level was observed upon the silencing of these two splicing factors in TC-71 cells (Figure 5A,B). Similar results were also obtained in SK-N-MC and LAP-35 Ewing sarcoma cells (Supplementary Figure S3), even though HNRNPM and SRSF3 knockdown in LAP-35 weakly affected DHX9 alternative splicing.
As previously reported [11,18,20,21,22,23], disruption of the DHX9-EWS-FLI1 interaction strongly impacts on the EWS-FLI1-driven transcriptional program in Ewing sarcoma cells. Thus, we asked whether knockdown of SRSF3 and hnRNPM had a similar effect on EWS-FLI1 target genes. RT-qPCR analysis showed that expression of c-MYC, EZH2, and ID2, was downregulated upon silencing of the two RBPs, with stronger effects elicited by depletion of SRSF3. Moreover, expression of NR0B1 and CCND1 was affected only by SRSF3 depletion (Figure 5C).
To evaluate the effect of these splicing factors on Ewing sarcoma cell proliferation, we performed colony assays with cells transfected with control siRNA or with si-HNRNPM and si-SRSF3 oligonucleotides. Depletion of hnRNPM, and even more of SRSF3, caused a significant reduction in the colony formation potential of Ewing Sarcoma cells (Figure 5D). These findings demonstrate that hnRNPM and SRSF3 promote DHX9 expression and Ewing sarcoma cells proliferation and viability, suggesting that targeting their expression or activity may have beneficial effects.

3.6. Depletion of SRSF3 and hnRNPM Increases Doxorubicin Sensitivity of Ewing Sarcoma Cells

DHX9 helicase has been associated to cell protection from genome instability [11,12,16,17,29,30]. To test whether SRSF3 and HNRNPM depletion had an effect on Ewing sarcoma sensitivity to chemotherapeutic agents, we treated TC-71 cells with doxorubicin. This drug is commonly used in the clinic as a genotoxic agent and is included in the chemotherapeutic regimen for Ewing sarcoma patients [31]. Doxorubicin induces DNA double-strand breaks and can cause a nonreversible checkpoint arrest or trigger cell death, thus curbing the rapid proliferation of cancer cells. MTS assays showed a statistically significant decrease in the viability of TC-71 cells silenced for HNRNPM (red) and SRSF3 (blue) versus control (black) (Figure 6A). In particular, the EC50 for doxorubicin activity was 87.82 nM in control transfected TC-71 cells, whereas it was reduced to 48.88 nM and 12.12 nM in siHNRNPM and siSRSF3 cells, respectively (Figure 6A). Similar results were obtained in SK-N-MC (Supplementary Figure S4A), whereas LAP-35 Ewing sarcoma cells displayed only a slight decrease in the EC50, in line with the more modest effect on DHX9 exon 6A inclusion (Supplementary Figure S4B). To validate these results by a different assay, we performed clonogenic assays in TC-71 cells exposed to an increasing amount of doxorubicin (from 0.1 nM to 10 nM). The knockdown of hnRNPM and SRSF3 significantly reduced the ability of TC-71 cells to form colonies in the presence of the drug (Figure 6B). In particular, at 10 nM concentration, we observed a reduction in the percentage of colonies from 54% for control, to 37.5% and 32% for siHNRMPM and siSRSF3, respectively.
Collectively these results unveil a novel role for SRSF3 and hnRNPM in the regulation of DHX9 alternative splicing, which also impacts on Ewing sarcoma cell sensitivity to chemotherapeutic treatments.

4. Discussion

Ewing sarcoma treatment relies on a multidisciplinary approach that combines multi-drug chemotherapy with surgery and local radiation therapy. Such treatments effectively reduce the risk of recurrence and increase the overall survival of patients to 65%. However, pediatric patients have to face severe long-term toxicities due to these invasive therapies. Thus, the development of tailored therapies based on new valuable prognostic markers and therapeutic targets is urgently needed. In this regard, the splicing signature of human cancers is emerging as a powerful tool to distinguish tumor subtypes and stratify patients [32,33]. Notably, although alternative splicing dysregulation has also been reported in Ewing sarcoma [24,34,35,36,37,38], limited information is available regarding the RBPs responsible for this process and their possible association with prognosis. An important regulator of RNA metabolism with strong implications in Ewing sarcoma malignancy is the DNA/RNA helicase DHX9, which interacts with the oncogene EWS-FLI1 and promotes its transcriptional activity. Herein, by querying Ewing sarcoma datasets of patients, we found that the expression of DHX9 is positively associated with disease progression and worse prognosis. Notably, DHX9 interacts with several splicing factors and is itself regulated by alternative splicing [10,11]. Indeed, the inclusion of a poison-exon in the DHX9 transcript targets it to NMD and lowers DHX9 expression level [11]. Nevertheless, the mechanism that suppresses the inclusion of this poison exon and supports high DHX9 expression in Ewing sarcoma cells was unknown. In this study, by employing multiple approaches, we have now identified several splicing factors (SRSF1, SRSF3, and SRSF10, hnRNPK, hnRNPM, and FUS) as regulators of DHX9 splicing in Ewing sarcoma cells. Moreover, our work indicates that the silencing of hnRNPM and SRSF3, the strongest repressors of exon 6A, significantly reduces DHX9 expression and Ewing sarcoma cell survival, thus suggesting the functional relevance of their effects on DHX9 splicing regulation.
In our study, we focused on members of the SR proteins and hnRNPs families that were previously shown to play oncogenic functions [27,39,40,41]. Bioinformatic analysis of exon 6A and flanking introns identified the presence of putative binding sites for all the SR and hnRNP proteins that displayed a significant regulation of DHX9 exon 6A splicing. CLIP experiments confirmed the binding of SRSF1, SRSF3, hnRNPK, hnRNPM, and FUS to the exon 6A of DHX9 (Figure 3). We also found that the expression of hnRNPM and SRSF3, which showed the strongest repression of exon 6A inclusion, is highly associated with the expression of DHX9 in Ewing sarcoma patients. Remarkably, silencing of SRSF3 and HNRNPM, impacts on DHX9 expression as well as on Ewing sarcoma cell viability and proliferation, while it reduces survival of Ewing sarcoma cells to doxorubicin treatment.
Several studies have previously implicated hnRNPM and SRSF3 in human cancers. hnRNPM is a well-known splicing factor that binds to GU-rich cis-elements in target RNAs [42,43,44,45]. In breast cancer, hnRNPM contributes to metastasis by activating alternative splicing changes that promote epithelial–mesenchymal transition (EMT) [46]. In Ewing sarcoma, we previously documented that the inhibition of the mTOR/AKT/PI3K pathway sets in motion an hnRNPM-dependent alternative splicing program that limits the therapeutic efficacy of these pharmacologic inhibitors [24]. On this basis, we indicated hnRNPM as a new potential therapeutic target to counteract Ewing sarcoma malignancy [24]. Notably, our transcriptomic analysis upon inhibition of the mTOR/AKT/PI3K pathway documented the upregulation of both DHX9 and hnRNPM, and downregulation of DHX9 exon 6A [24], suggesting a direct regulation of DHX9 alternative splicing by hnRNPM. Herein, we demonstrate that hnRNPM directly binds to intron 5 and exon 6A in the DHX9 pre-mRNA, resulting in repression of exon 6A inclusion and stabilization of DHX9 expression. Thus, our results suggest that hnRNPM contributes to Ewing sarcoma malignancy by promoting the high expression of a key regulator of the EWS-FLI1 oncogene.
SRSF3 is the smallest member of the SR family of proteins and plays important roles in the regulation of RNA metabolism, including splicing [47], RNA export [48] and polyadenylation [49], protein translation [50,51], pri-miRNA processing [52], and genome stability [53]. SRSF3 is frequently upregulated in cancer and displays pro-oncogenic functions [54,55,56], whereas tumors displaying reduced SRSF3 expression exhibited slower growth and higher apoptosis [56]. Accordingly, SRSF3 depleted cells were less tumorigenic in nude mice [56]. In particular, SRSF3 was shown to regulate the alternative splicing of genes related to cell proliferation and cell cycle progression in osteosarcoma [55]. Our study now suggests an oncogenic role of SRSF3 also in Ewing sarcoma cells. SRSF3 binds and regulates the DHX9 exon 6A similarly to hnRNPM, thus also promoting DHX9 expression. Notably, the effect of SRSF3 knockdown on DHX9 protein expression and Ewing sarcoma cell viability was stronger than that elicited by hnRNPM, although displaying a milder effect on the inclusion of the alternative poison exon. This apparent discrepancy could be explained by a direct role played by SRSF3 in the NMD process. Indeed, it was previously shown that SRSF3 depletion in pluripotent cells leads to the downregulation of hundreds of mRNAs that are also regulated at the level of NMD, suggesting a link between SRSF3 and mRNA surveillance pathway [57,58]. However, the exact role of SRSF3 in NMD has not been elucidated yet.
Our study documents that SRSF3 and hnRNPM downregulation reduces Ewing sarcoma cell proliferation and increases Ewing sarcoma sensitivity to doxorubicin treatment (Figure 6). Ewing sarcoma patients are generally treated with a potent cocktail of five drugs, including doxorubicin. Thus, targeting SRSF3 and hnRNPM expression, could have a potential therapeutic value by sensitizing cells to drug treatment, and could, therefore, be exploited in a combination therapy. SRSF3 depletion displayed a stronger effect than HNRNPM depletion on Ewing sarcoma sensitivity to chemotherapy treatment, in line with its stronger impact on DHX9 expression. Notably, the chromosome region bearing the SRSF3 locus on chromosome 6p is commonly amplified in cancer [59], thus causing aberrant SRSF3 overexpression. Moreover, the downregulation of SRSF3 by antisense oligonucleotides (ASO) was shown to sensitize oral squamous cell carcinoma and breast cancer cells to paclitaxel treatment [60]. These findings suggest that the depletion of SRSF3 expression may represent a valuable therapeutic tool also for Ewing sarcoma. Noteworthy, although we focused on the DHX9 poison-exon as a splicing target of SRSF3 and hnRNPM with direct relevance for Ewing sarcoma, our findings do not rule out possible effects of these RBPs on other splicing variants of the same or other genes in this disease.

5. Conclusions

In conclusion, our study shows that DHX9 expression in Ewing sarcoma correlates with worse outcome in patients. SRSF3 and hnRNPM regulate poison-exon inclusion in DHX9 and sensitize Ewing sarcoma cells to chemotherapy. Thus, we suggest that modulation of SRSF3 and hnRNPM expression and their splicing signature could represent a novel therapeutic opportunity for combined and less aggressive therapy in Ewing sarcoma.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4409/9/2/328/s1, Supplementary Figure S1 RT-qPCR analysis to monitor the inclusion of DHX9 exon 6A normalized to the skipped isoform (junction exon 6-exon 7), upon knockdown of the indicated SR (A) and hnRNP (B) proteins, Supplementary Figure S2 Pearson correlation analysis on Ewing sarcoma patients, between the expression of DHX9 and either hnRNPL (A) or SRSF5 (B), Supplementary Figure S3. HNRNPM and SRSF3 affect DHX9 exon 6A inclusion in SK-N-MC and LAP-35 Ewing sarcoma cells, Supplementary Figure S4. Depletion of SRSF3 and hnRNPM increases doxorubicin sensitivity of SK-N-MC and LAP-35 Ewing sarcoma cells, Supplementary Table S1 List of oligonucleotides used in the siRNA library, Supplementary Table S2 List of primers.

Author Contributions

Conceptualization, M.P.P.; methodology, R.P., V.V.; software, R.P., V.V., and M.P.P.; validation, R.P. and V.V.; formal analysis, R.P., V.V., and M.P.P.; investigation, R.P., V.V., and M.P.P.; data curation, R.P., V.V., and M.P.P.; writing—original draft preparation, M.P.P.; writing—review and editing, R.P., V.V., and M.P.P.; supervision, M.P.P.; project administration, M.P.P.; funding acquisition, M.P.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Associazione Italiana Ricerca sul Cancro (AIRC), grant number IG17278 and the Ministry of Health “Ricerca Finalizzata”, grant number RF-2016-02363460 to M.P.P. R.P. was supported by a scholarship from Fondazione Umberto Veronesi.

Acknowledgments

The authors wish to thank Gabriele Ferrante and Gloria Guizzo for technical assistance and Prof. Claudio Sette for helpful discussion.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Graveley, B.R. Alternative splicing: Increasing diversity in the proteomic world. Trends Genet. 2001, 17, 100–107. [Google Scholar] [CrossRef]
  2. Black, D.L. Mechanisms of alternative pre-messenger RNA splicing. Annu. Rev. Biochem 2003, 72, 291–336. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Paronetto, M.P.; Passacantilli, I.; Sette, C. Alternative splicing and cell survival: From tissue homeostasis to disease. Cell Death Differ. 2016, 23, 1919–1929. [Google Scholar] [CrossRef] [PubMed]
  4. Cartegni, L.; Krainer, A.R. Disruption of an SF2/ASF-dependent exonic splicing enhancer in SMN2 causes spinal muscular atrophy in the absence of SMN1. Nat. Genet. 2002, 30, 377–384. [Google Scholar] [CrossRef] [PubMed]
  5. Baker, K.E.; Parker, R. Nonsense-mediated mRNA decay: Terminating erroneous gene expression. Curr. Opin. Cell Biol. 2004, 16, 293–299. [Google Scholar] [CrossRef] [PubMed]
  6. 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] [Green Version]
  7. Ruiz-Echevarría, M.J.; González, C.I.; Peltz, S.W. Identifying the right stop: Determining how the surveillance complex recognizes and degrades an aberrant mRNA. EMBO J. 1998, 17, 575–589. [Google Scholar] [CrossRef] [Green Version]
  8. Valacca, C.; Bonomi, S.; Buratti, E.; Pedrotti, S.; Baralle, F.E.; Sette, C.; Ghigna, C.; Biamonti, G. Sam68 regulates EMT through alternative splicing-activated nonsense-mediated mRNA decay of the SF2/ASF proto-oncogene. J. Cell Biol. 2010, 191, 87–99. [Google Scholar] [CrossRef]
  9. Jacob, A.G.; Smith, C.W.J. Intron retention as a component of regulated gene expression programs. Hum. Genet. 2017, 136, 1043–1057. [Google Scholar] [CrossRef] [Green Version]
  10. Ip, J.Y.; Schmidt, D.; Pan, Q.; Ramani, A.K.; Fraser, A.G.; Odom, D.T.; Blencowe, B.J. Global impact of RNA polymerase II elongation inhibition on alternative splicing regulation. Genome Res. 2011, 21, 390–401. [Google Scholar] [CrossRef] [Green Version]
  11. Fidaleo, M.; Svetoni, F.; Volpe, E.; Miñana, B.; Caporossi, D.; Paronetto, M.P. Genotoxic stress inhibits Ewing sarcoma cell growth by modulating alternative pre-mRNA processing of the RNA helicase DHX9. Oncotarget 2015, 6, 31740–31757. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Fidaleo, M.; De Paola, E.; Paronetto, M.P. The RNA helicase A in malignant transformation. Oncotarget 2016, 7, 28711–28723. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Anderson, S.F.; Schlegel, B.P.; Nakajima, T.; Wolpin, E.S.; Parvin, J.D. BRCA1 protein is linked to the RNA polymerase II holoenzyme complex via RNA helicase A. Nat. Genet. 1998, 19, 254–256. [Google Scholar] [CrossRef]
  14. Aratani, S.; Fujii, R.; Oishi, T.; Fujita, H.; Amano, T.; Ohshima, T.; Hagiwara, M.; Fukamizu, A.; Nakajima, T. Dual roles of RNA helicase A in CREB-dependent transcription. Mol. Cell Biol. 2001, 21, 4460–4469. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Jain, A.; Bacolla, A.; Chakraborty, P.; Grosse, F.; Vasquez, K.M. Human DHX9 helicase unwinds triple-helical DNA structures. Biochemistry 2010, 49, 6992–6999. [Google Scholar] [CrossRef] [Green Version]
  16. Jain, A.; Bacolla, A.; Del Mundo, I.M.; Zhao, J.; Wang, G.; Vasquez, K.M. DHX9 helicase is involved in preventing genomic instability induced by alternatively structured DNA in human cells. Nucleic Acids Res. 2013, 41, 10345–10357. [Google Scholar] [CrossRef] [Green Version]
  17. Chakraborty, P.; Grosse, F. WRN helicase unwinds Okazaki fragment-like hybrids in a reaction stimulated by the human DHX9 helicase. Nucleic Acids Res. 2010, 38, 4722–4730. [Google Scholar] [CrossRef] [Green Version]
  18. Toretsky, J.A.; Erkizan, V.; Levenson, A.; Abaan, O.D.; Parvin, J.D.; Cripe, T.P.; Rice, A.M.; Lee, S.B.; Uren, A. Oncoprotein EWS-FLI1 activity is enhanced by RNA helicase A. Cancer Res. 2006, 66, 5574–5581. [Google Scholar] [CrossRef] [Green Version]
  19. Palombo, R.; Frisone, P.; Fidaleo, M.; Mercatelli, N.; Sette, C.; Paronetto, M.P. The Promoter-Associated Noncoding RNA. Cancer Res. 2019, 79, 3570–3582. [Google Scholar] [CrossRef] [Green Version]
  20. Mercatelli, N.; Fortini, D.; Palombo, R.; Paronetto, M.P. Small molecule inhibition of Ewing sarcoma cell growth via targeting the long non coding RNA HULC. Cancer Lett. 2019. [Google Scholar] [CrossRef]
  21. Erkizan, H.V.; Kong, Y.; Merchant, M.; Schlottmann, S.; Barber-Rotenberg, J.S.; Yuan, L.; Abaan, O.D.; Chou, T.H.; Dakshanamurthy, S.; Brown, M.L.; et al. A small molecule blocking oncogenic protein EWS-FLI1 interaction with RNA helicase A inhibits growth of Ewing’s sarcoma. Nat. Med. 2009, 15, 750–756. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Hong, S.H.; Youbi, S.E.; Hong, S.P.; Kallakury, B.; Monroe, P.; Erkizan, H.V.; Barber-Rotenberg, J.S.; Houghton, P.; Üren, A.; Toretsky, J.A. Pharmacokinetic modeling optimizes inhibition of the ‘undruggable’ EWS-FLI1 transcription factor in Ewing Sarcoma. Oncotarget 2014, 5, 338–350. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Zöllner, S.K.; Selvanathan, S.P.; Graham, G.T.; Commins, R.M.T.; Hong, S.H.; Moseley, E.; Parks, S.; Haladyna, J.N.; Erkizan, H.V.; Dirksen, U.; et al. Inhibition of the oncogenic fusion protein EWS-FLI1 causes G. Sci. Signal. 2017, 10. [Google Scholar] [CrossRef] [Green Version]
  24. Passacantilli, I.; Frisone, P.; De Paola, E.; Fidaleo, M.; Paronetto, M.P. hnRNPM guides an alternative splicing program in response to inhibition of the PI3K/AKT/mTOR pathway in Ewing sarcoma cells. Nucleic Acids Res. 2017, 45, 12270–12284. [Google Scholar] [CrossRef]
  25. Savola, S.; Klami, A.; Tripathi, A.; Niini, T.; Serra, M.; Picci, P.; Kaski, S.; Zambelli, D.; Scotlandi, K.; Knuutila, S. Combined use of expression and CGH arrays pinpoints novel candidate genes in Ewing sarcoma family of tumors. BMC Cancer 2009, 9, 17. [Google Scholar] [CrossRef] [Green Version]
  26. Available online: http://cbio.mskcc.org/Public/sarcoma_array_data/ (accessed on 1 December 2011).
  27. Smith, C.W.; Valcárcel, J. Alternative pre-mRNA splicing: The logic of combinatorial control. Trends Biochem. Sci. 2000, 25, 381–388. [Google Scholar] [CrossRef]
  28. Cáceres, J.F.; Kornblihtt, A.R. Alternative splicing: Multiple control mechanisms and involvement in human disease. Trends Genet. 2002, 18, 186–193. [Google Scholar] [CrossRef]
  29. Chakraborty, P.; Huang, J.T.J.; Hiom, K. DHX9 helicase promotes R-loop formation in cells with impaired RNA splicing. Nat. Commun. 2018, 9, 4346. [Google Scholar] [CrossRef]
  30. Aktaş, T.; Avşar Ilık, İ.; Maticzka, D.; Bhardwaj, V.; Pessoa Rodrigues, C.; Mittler, G.; Manke, T.; Backofen, R.; Akhtar, A. DHX9 suppresses RNA processing defects originating from the Alu invasion of the human genome. Nature 2017, 544, 115–119. [Google Scholar] [CrossRef]
  31. Balamuth, N.J.; Womer, R.B. Ewing’s sarcoma. Lancet Oncol. 2010, 11, 184–192. [Google Scholar] [CrossRef]
  32. Trincado, J.L.; Sebestyén, E.; Pagés, A.; Eyras, E. The prognostic potential of alternative transcript isoforms across human tumors. Genome Med. 2016, 8, 85. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Wu, H.Y.; Peng, Z.G.; He, R.Q.; Luo, B.; Ma, J.; Hu, X.H.; Dang, Y.W.; Chen, G.; Pan, S.L. Prognostic index of aberrant mRNA splicing profiling acts as a predictive indicator for hepatocellular carcinoma based on TCGA SpliceSeq data. Int. J. Oncol. 2019, 55, 425–438. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Selvanathan, S.P.; Graham, G.T.; Grego, A.R.; Baker, T.M.; Hogg, J.R.; Simpson, M.; Batish, M.; Crompton, B.; Stegmaier, K.; Tomazou, E.M.; et al. EWS-FLI1 modulated alternative splicing of ARID1A reveals novel oncogenic function through the BAF complex. Nucleic Acids Res. 2019, 47, 9619–9636. [Google Scholar] [CrossRef] [PubMed]
  35. Selvanathan, S.P.; Graham, G.T.; Erkizan, H.V.; Dirksen, U.; Natarajan, T.G.; Dakic, A.; Yu, S.; Liu, X.; Paulsen, M.T.; Ljungman, M.E.; et al. Oncogenic fusion protein EWS-FLI1 is a network hub that regulates alternative splicing. Proc. Natl. Acad. Sci. USA 2015, 112, E1307–E1316. [Google Scholar] [CrossRef] [Green Version]
  36. Sand, L.G.; Szuhai, K.; Hogendoorn, P.C. Sequencing Overview of Ewing Sarcoma: A Journey across Genomic, Epigenomic and Transcriptomic Landscapes. Int. J. Mol. Sci. 2015, 16, 16176–16215. [Google Scholar] [CrossRef] [Green Version]
  37. Paronetto, M.P.; Bernardis, I.; Volpe, E.; Bechara, E.; Sebestyén, E.; Eyras, E.; Valcárcel, J. Regulation of FAS exon definition and apoptosis by the Ewing sarcoma protein. Cell Rep. 2014, 7, 1211–1226. [Google Scholar] [CrossRef] [Green Version]
  38. Sanchez, G.; Bittencourt, D.; Laud, K.; Barbier, J.; Delattre, O.; Auboeuf, D.; Dutertre, M. Alteration of cyclin D1 transcript elongation by a mutated transcription factor up-regulates the oncogenic D1b splice isoform in cancer. Proc. Natl. Acad. Sci. USA 2008, 105, 6004–6009. [Google Scholar] [CrossRef] [Green Version]
  39. 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]
  40. Obeng, E.A.; Stewart, C.; Abdel-Wahab, O. Altered RNA Processing in Cancer Pathogenesis and Therapy. Cancer Discov. 2019, 9, 1493–1510. [Google Scholar] [CrossRef] [Green Version]
  41. Anczuków, O.; Krainer, A.R. Splicing-factor alterations in cancers. RNA 2016, 22, 1285–1301. [Google Scholar] [CrossRef] [Green Version]
  42. Swanson, M.S.; Dreyfuss, G. RNA binding specificity of hnRNP proteins: A subset bind to the 3′ end of introns. EMBO J. 1988, 7, 3519–3529. [Google Scholar] [CrossRef] [PubMed]
  43. Datar, K.V.; Dreyfuss, G.; Swanson, M.S. The human hnRNP M proteins: Identification of a methionine/arginine-rich repeat motif in ribonucleoproteins. Nucleic Acids Res. 1993, 21, 439–446. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Hovhannisyan, R.H.; Carstens, R.P. Heterogeneous ribonucleoprotein m is a splicing regulatory protein that can enhance or silence splicing of alternatively spliced exons. J. Biol. Chem. 2007, 282, 36265–36274. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Huelga, S.C.; Vu, A.Q.; Arnold, J.D.; Liang, T.Y.; Liu, P.P.; Yan, B.Y.; Donohue, J.P.; Shiue, L.; Hoon, S.; Brenner, S.; et al. Integrative genome-wide analysis reveals cooperative regulation of alternative splicing by hnRNP proteins. Cell Rep. 2012, 1, 167–178. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Xu, Y.; Gao, X.D.; Lee, J.H.; Huang, H.; Tan, H.; Ahn, J.; Reinke, L.M.; Peter, M.E.; Feng, Y.; Gius, D.; et al. Cell type-restricted activity of hnRNPM promotes breast cancer metastasis via regulating alternative splicing. Genes Dev. 2014, 28, 1191–1203. [Google Scholar] [CrossRef] [Green Version]
  47. Zahler, A.M.; Lane, W.S.; Stolk, J.A.; Roth, M.B. SR proteins: A conserved family of pre-mRNA splicing factors. Genes Dev. 1992, 6, 837–847. [Google Scholar] [CrossRef] [Green Version]
  48. Huang, Y.; Steitz, J.A. Splicing factors SRp20 and 9G8 promote the nucleocytoplasmic export of mRNA. Mol. Cell 2001, 7, 899–905. [Google Scholar] [CrossRef]
  49. Lou, H.; Neugebauer, K.M.; Gagel, R.F.; Berget, S.M. Regulation of alternative polyadenylation by U1 snRNPs and SRp20. Mol. Cell Biol. 1998, 18, 4977–4985. [Google Scholar] [CrossRef] [Green Version]
  50. Kim, J.; Park, R.Y.; Chen, J.K.; Jeong, S.; Ohn, T. Splicing factor SRSF3 represses the translation of programmed cell death 4 mRNA by associating with the 5′-UTR region. Cell Death Differ. 2014, 21, 481–490. [Google Scholar] [CrossRef]
  51. Bedard, K.M.; Daijogo, S.; Semler, B.L. A nucleo-cytoplasmic SR protein functions in viral IRES-mediated translation initiation. EMBO J. 2007, 26, 459–467. [Google Scholar] [CrossRef] [Green Version]
  52. Auyeung, V.C.; Ulitsky, I.; McGeary, S.E.; Bartel, D.P. Beyond secondary structure: Primary-sequence determinants license pri-miRNA hairpins for processing. Cell 2013, 152, 844–858. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Loomis, R.J.; Naoe, Y.; Parker, J.B.; Savic, V.; Bozovsky, M.R.; Macfarlan, T.; Manley, J.L.; Chakravarti, D. Chromatin binding of SRp20 and ASF/SF2 and dissociation from mitotic chromosomes is modulated by histone H3 serine 10 phosphorylation. Mol. Cell 2009, 33, 450–461. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. He, X.; Arslan, A.D.; Pool, M.D.; Ho, T.T.; Darcy, K.M.; Coon, J.S.; Beck, W.T. Knockdown of splicing factor SRp20 causes apoptosis in ovarian cancer cells and its expression is associated with malignancy of epithelial ovarian cancer. Oncogene 2011, 30, 356–365. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Ajiro, M.; Jia, R.; Yang, Y.; Zhu, J.; Zheng, Z.M. A genome landscape of SRSF3-regulated splicing events and gene expression in human osteosarcoma U2OS cells. Nucleic Acids Res. 2016, 44, 1854–1870. [Google Scholar] [CrossRef]
  56. Jia, R.; Li, C.; McCoy, J.P.; Deng, C.X.; Zheng, Z.M. SRp20 is a proto-oncogene critical for cell proliferation and tumor induction and maintenance. Int. J. Biol. Sci. 2010, 6, 806–826. [Google Scholar] [CrossRef] [PubMed]
  57. Hurt, J.A.; Robertson, A.D.; Burge, C.B. Global analyses of UPF1 binding and function reveal expanded scope of nonsense-mediated mRNA decay. Genome Res. 2013, 23, 1636–1650. [Google Scholar] [CrossRef] [Green Version]
  58. Ratnadiwakara, M.; Archer, S.K.; Dent, C.I.; Ruiz De Los Mozos, I.; Beilharz, T.H.; Knaupp, A.S.; Nefzger, C.M.; Polo, J.M.; Anko, M.L. SRSF3 promotes pluripotency through. Elife 2018, 7. [Google Scholar] [CrossRef]
  59. Santos, G.C.; Zielenska, M.; Prasad, M.; Squire, J.A. Chromosome 6p amplification and cancer progression. J. Clin. Pathol. 2007, 60, 1–7. [Google Scholar] [CrossRef] [Green Version]
  60. Sun, Y.; Yan, L.; Guo, J.; Shao, J.; Jia, R. Downregulation of SRSF3 by antisense oligonucleotides sensitizes oral squamous cell carcinoma and breast cancer cells to paclitaxel treatment. Cancer Chemother Pharm. 2019, 84, 1133–1143. [Google Scholar] [CrossRef]
Figure 1. DHX9 expression is a prognostic factor for Ewing sarcoma malignancy. Kaplan–Meier plots of overall (A) and event-free survival (B) probabilities generated by gene expression profiling ([25]; GSE17679) in R2 Genomics (https://hgserver1.amc.nl/cgi-bin/r2/main.cgi). The dataset analyzed is composed of 64 Ewing sarcoma, 4 Askin, and 20 PNET (Peripheral Primitive Neuro-Ectodermal Tumors) patients. P-values and Bonferroni post-hoc corrections are indicated on the bottom. (C) The box plot shows DHX9 expression in Ewing sarcoma (64 patients), Askin tumor (4 patients), and PNET (20 patients). Statistical analysis was performed by one-way ANOVA. (D) DHX9 expression was monitored by stratifying the previous dataset in primary and metastatic tumors ([25] GSE17679). Statistical analysis was performed by Student’s t-test (p-value is indicated on the top of the graph).
Figure 1. DHX9 expression is a prognostic factor for Ewing sarcoma malignancy. Kaplan–Meier plots of overall (A) and event-free survival (B) probabilities generated by gene expression profiling ([25]; GSE17679) in R2 Genomics (https://hgserver1.amc.nl/cgi-bin/r2/main.cgi). The dataset analyzed is composed of 64 Ewing sarcoma, 4 Askin, and 20 PNET (Peripheral Primitive Neuro-Ectodermal Tumors) patients. P-values and Bonferroni post-hoc corrections are indicated on the bottom. (C) The box plot shows DHX9 expression in Ewing sarcoma (64 patients), Askin tumor (4 patients), and PNET (20 patients). Statistical analysis was performed by one-way ANOVA. (D) DHX9 expression was monitored by stratifying the previous dataset in primary and metastatic tumors ([25] GSE17679). Statistical analysis was performed by Student’s t-test (p-value is indicated on the top of the graph).
Cells 09 00328 g001
Figure 2. SRSF3 and hnRNPM expression in ES patients exhibit a significant correlation with DHX9 expression. Pearson correlation analysis on Ewing sarcoma patients, performed between the expression of DHX9 and SRSF1 (A), SRSF3 (B), SRSF10 (C), hnRNPK (D), hnRNPM (E), or FUS (F). Values are expressed as base 2-logarithm. In each panel, the correlation value (R) and the relative p-value are indicated.
Figure 2. SRSF3 and hnRNPM expression in ES patients exhibit a significant correlation with DHX9 expression. Pearson correlation analysis on Ewing sarcoma patients, performed between the expression of DHX9 and SRSF1 (A), SRSF3 (B), SRSF10 (C), hnRNPK (D), hnRNPM (E), or FUS (F). Values are expressed as base 2-logarithm. In each panel, the correlation value (R) and the relative p-value are indicated.
Cells 09 00328 g002
Figure 3. Regulation of DHX9 alternative splicing. (A) Schematic representation of the alternative splicing of exon 6A in DHX9 pre-mRNA. Exclusion of the poison exon 6A in DHX9 mRNA leads to the main transcript (NM001357) encoding the full-length DHX9 protein (upper part). Exon 6A inclusion leads to the alternative noncoding transcript NR033302, containing a premature stop codon (PTC) and targeted to the NMD machinery. (BE) Histograms represents RT-qPCR analysis of a siRNA library to downregulate the expression of SR proteins (B) (one-way ANOVA p-value < 0.0001; Bonferroni correction for siSRSF6 p-value < 0.05; siSRSF1 and siSRSF5 p-value < 0.005; siSRSF3, siSRSF4, siSRSF7, siSRSF9 p-value < 0.001, siSRSF2 and siSRSF10 p-value > 0.05) or hnRNPs (C) (ANOVA p-value < 0.0001; Bonferroni correction for sihnRNPF and sihnRNPH1 p-value < 0.05; siHNRNPAB, siHNRNPAO, siHNRNPC, siHNRNPK, siHNRNPI, siHNRNPM, siHNRNPP, and siHNRNPR p-value < 0.001, siHNRNPU p-value > 0.05). DHX9 Exon 6A (Ex6a) inclusion was monitored by RT-qPCR and normalized to the constitutive exon 4 (D,E). Reported values represent the average (± S.D.) of at least three independent experiments. Statistical analysis was performed by one-way ANOVA, with p-value < 0.0001 (D,E) and with Bonferroni post-hoc test. (p-value: *** < 0.001, ** < 0.01, *< 0.05).
Figure 3. Regulation of DHX9 alternative splicing. (A) Schematic representation of the alternative splicing of exon 6A in DHX9 pre-mRNA. Exclusion of the poison exon 6A in DHX9 mRNA leads to the main transcript (NM001357) encoding the full-length DHX9 protein (upper part). Exon 6A inclusion leads to the alternative noncoding transcript NR033302, containing a premature stop codon (PTC) and targeted to the NMD machinery. (BE) Histograms represents RT-qPCR analysis of a siRNA library to downregulate the expression of SR proteins (B) (one-way ANOVA p-value < 0.0001; Bonferroni correction for siSRSF6 p-value < 0.05; siSRSF1 and siSRSF5 p-value < 0.005; siSRSF3, siSRSF4, siSRSF7, siSRSF9 p-value < 0.001, siSRSF2 and siSRSF10 p-value > 0.05) or hnRNPs (C) (ANOVA p-value < 0.0001; Bonferroni correction for sihnRNPF and sihnRNPH1 p-value < 0.05; siHNRNPAB, siHNRNPAO, siHNRNPC, siHNRNPK, siHNRNPI, siHNRNPM, siHNRNPP, and siHNRNPR p-value < 0.001, siHNRNPU p-value > 0.05). DHX9 Exon 6A (Ex6a) inclusion was monitored by RT-qPCR and normalized to the constitutive exon 4 (D,E). Reported values represent the average (± S.D.) of at least three independent experiments. Statistical analysis was performed by one-way ANOVA, with p-value < 0.0001 (D,E) and with Bonferroni post-hoc test. (p-value: *** < 0.001, ** < 0.01, *< 0.05).
Cells 09 00328 g003
Figure 4. SRSF1, SRSF3, hnRNPK, hnRNPM, and FUS bind DHX9 exon 6A. (A) In silico analysis was performed using the Splice Aid database to identify putative consensus motifs for RBPs in the exon 6A sequence. In (A), schematic representation of DHX9 alternative exon 6A, with the position of SRSF1 (blue), SRSF3 (red), hnRNPK (green), hnRNPM (grey), and FUS (yellow) binding sites. Arrows indicate the primers’ positions used for the amplicons along the sequence. (B) TC-71 cells were UV-crosslinked, and protein-RNA extracts were immunoprecipitated with either control IgGs or hnRNPM, hnRNPK, FUS, SRSF1, and SRSF3 antibodies. Experiments were performed at least three times. Histograms show RT-qPCR analysis of the cross-linked and immunoprecipitation (CLIP) assay to analyze the binding of the indicated RNA binding proteins (RBPs) to intron 6 and exon 6A of DHX9 pre-mRNA. Results are expressed as a percentage of input (± S.D.). Statistical analysis was performed by one-way ANOVA with Bonferroni post-hoc test (p-value: *** < 0.001).
Figure 4. SRSF1, SRSF3, hnRNPK, hnRNPM, and FUS bind DHX9 exon 6A. (A) In silico analysis was performed using the Splice Aid database to identify putative consensus motifs for RBPs in the exon 6A sequence. In (A), schematic representation of DHX9 alternative exon 6A, with the position of SRSF1 (blue), SRSF3 (red), hnRNPK (green), hnRNPM (grey), and FUS (yellow) binding sites. Arrows indicate the primers’ positions used for the amplicons along the sequence. (B) TC-71 cells were UV-crosslinked, and protein-RNA extracts were immunoprecipitated with either control IgGs or hnRNPM, hnRNPK, FUS, SRSF1, and SRSF3 antibodies. Experiments were performed at least three times. Histograms show RT-qPCR analysis of the cross-linked and immunoprecipitation (CLIP) assay to analyze the binding of the indicated RNA binding proteins (RBPs) to intron 6 and exon 6A of DHX9 pre-mRNA. Results are expressed as a percentage of input (± S.D.). Statistical analysis was performed by one-way ANOVA with Bonferroni post-hoc test (p-value: *** < 0.001).
Cells 09 00328 g004
Figure 5. DHX9 exon 6A inclusion affects Ewing sarcoma cells’ viability and proliferation. Western blot analysis of TC-71 cell extracts knocked down for either hnRNPM (A) or SRSF3 (B). DHX9 protein levels were monitored and normalized to β-actin. Histogram represents the densitometric analysis of DHX9 expression from three independent experiments (average ± S.D.). Statistical analysis was performed by Student t-test (p-value *** < 0.001). (C) RT-qPCR analysis to detect the expression of EWS-FLI1 target genes (CCND1, c-MYC, EZH2, ID2, NR0B1) upon silencing of HNRNPM (upper part) and SRSF3 (lower part). (D) Representative clonogenic assay of cells transfected with either scrambled or siRNA oligonucleotides targeting HNRNPM and SRSF3 transcripts. On the bottom, the bar graph shows the percentage of colonies in each condition versus scrambled. Student t-test was used for statistical analysis (p-value: *** < 0.001, *<0.05; ##<0.01 siHNRNPM vs siSRSF3).
Figure 5. DHX9 exon 6A inclusion affects Ewing sarcoma cells’ viability and proliferation. Western blot analysis of TC-71 cell extracts knocked down for either hnRNPM (A) or SRSF3 (B). DHX9 protein levels were monitored and normalized to β-actin. Histogram represents the densitometric analysis of DHX9 expression from three independent experiments (average ± S.D.). Statistical analysis was performed by Student t-test (p-value *** < 0.001). (C) RT-qPCR analysis to detect the expression of EWS-FLI1 target genes (CCND1, c-MYC, EZH2, ID2, NR0B1) upon silencing of HNRNPM (upper part) and SRSF3 (lower part). (D) Representative clonogenic assay of cells transfected with either scrambled or siRNA oligonucleotides targeting HNRNPM and SRSF3 transcripts. On the bottom, the bar graph shows the percentage of colonies in each condition versus scrambled. Student t-test was used for statistical analysis (p-value: *** < 0.001, *<0.05; ##<0.01 siHNRNPM vs siSRSF3).
Cells 09 00328 g005
Figure 6. Depletion of SRSF3 and hnRNPM increases doxorubicin sensitivity of Ewing sarcoma cells. (A) Dose-response curve of scrambled (black), siHNRNPM (red), and siSRSF3 (blue) TC-71 cells after treatment with increasing concentration of doxorubicin (from 0 to 150 nM). Cells were collected at 72 h after treatment. EC50 values are reported on the bottom. (B) Colony assay of scrambled, siSRSF3, and siHNRNPM TC-71 cells treated with increasing concentration of doxorubicin (from 0.1 nM to 10 nM in DMSO). Histogram shows the percentage of colonies determined 12 days after treatment. Statistical analysis was performed by two-way ANOVA. Asterisks indicate significance with Bonferroni post-hoc test (p-value: ***< 0.001, ** < 0.01, *< 0.05). (C) Graphical representation of the hypothetical regulatory mechanism driving DHX9 alternative splicing. On the left, hnRNPM and SRSF3 bind DHX9 pre-mRNA to induce the skipping of exon 6A. The translated full-length DHX9 protein can interact with EWS-FLI1 to promote transcription of target genes involved in cell proliferation and transformation. On the right, in the absence of either SRSF3 or hnRNPM, DHX9 pre-mRNA is processed to the exon 6a-included noncoding transcript targeted to NMD. The consequent reduction in DHX9 protein levels impacts Ewing sarcoma cells proliferation and sensitizes Ewing sarcoma cells to doxorubicin treatment.
Figure 6. Depletion of SRSF3 and hnRNPM increases doxorubicin sensitivity of Ewing sarcoma cells. (A) Dose-response curve of scrambled (black), siHNRNPM (red), and siSRSF3 (blue) TC-71 cells after treatment with increasing concentration of doxorubicin (from 0 to 150 nM). Cells were collected at 72 h after treatment. EC50 values are reported on the bottom. (B) Colony assay of scrambled, siSRSF3, and siHNRNPM TC-71 cells treated with increasing concentration of doxorubicin (from 0.1 nM to 10 nM in DMSO). Histogram shows the percentage of colonies determined 12 days after treatment. Statistical analysis was performed by two-way ANOVA. Asterisks indicate significance with Bonferroni post-hoc test (p-value: ***< 0.001, ** < 0.01, *< 0.05). (C) Graphical representation of the hypothetical regulatory mechanism driving DHX9 alternative splicing. On the left, hnRNPM and SRSF3 bind DHX9 pre-mRNA to induce the skipping of exon 6A. The translated full-length DHX9 protein can interact with EWS-FLI1 to promote transcription of target genes involved in cell proliferation and transformation. On the right, in the absence of either SRSF3 or hnRNPM, DHX9 pre-mRNA is processed to the exon 6a-included noncoding transcript targeted to NMD. The consequent reduction in DHX9 protein levels impacts Ewing sarcoma cells proliferation and sensitizes Ewing sarcoma cells to doxorubicin treatment.
Cells 09 00328 g006

Share and Cite

MDPI and ACS Style

Palombo, R.; Verdile, V.; Paronetto, M.P. Poison-Exon Inclusion in DHX9 Reduces Its Expression and Sensitizes Ewing Sarcoma Cells to Chemotherapeutic Treatment. Cells 2020, 9, 328. https://doi.org/10.3390/cells9020328

AMA Style

Palombo R, Verdile V, Paronetto MP. Poison-Exon Inclusion in DHX9 Reduces Its Expression and Sensitizes Ewing Sarcoma Cells to Chemotherapeutic Treatment. Cells. 2020; 9(2):328. https://doi.org/10.3390/cells9020328

Chicago/Turabian Style

Palombo, Ramona, Veronica Verdile, and Maria Paola Paronetto. 2020. "Poison-Exon Inclusion in DHX9 Reduces Its Expression and Sensitizes Ewing Sarcoma Cells to Chemotherapeutic Treatment" Cells 9, no. 2: 328. https://doi.org/10.3390/cells9020328

APA Style

Palombo, R., Verdile, V., & Paronetto, M. P. (2020). Poison-Exon Inclusion in DHX9 Reduces Its Expression and Sensitizes Ewing Sarcoma Cells to Chemotherapeutic Treatment. Cells, 9(2), 328. https://doi.org/10.3390/cells9020328

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