- Research article
- Open Access
- Open Peer Review
Role of the splicing factor SRSF4 in cisplatin-induced modifications of pre-mRNA splicing and apoptosis
BMC Cancervolume 15, Article number: 227 (2015)
Modification of splicing by chemotherapeutic drugs has usually been evaluated on a limited number of pre-mRNAs selected for their recognized or potential importance in cell proliferation or apoptosis. However, the pathways linking splicing alterations to the efficiency of cancer therapy remain unclear.
Next-generation sequencing was used to analyse the transcriptome of breast carcinoma cells treated by cisplatin. Pharmacological inhibitors, RNA interference, cells deficient in specific signalling pathways, RT-PCR and FACS analysis were used to investigate how the anti-cancer drug cisplatin affected alternative splicing and the cell death pathway.
We identified 717 splicing events affected by cisplatin, including 245 events involving cassette exons. Gene ontology analysis indicates that cell cycle, mRNA processing and pre-mRNA splicing were the main pathways affected. Importantly, the cisplatin–induced splicing alterations required class I PI3Ks P110β but not components such as ATM, ATR and p53 that are involved in the DNA damage response. The siRNA-mediated depletion of the splicing regulator SRSF4, but not SRSF6, expression abrogated many of the splicing alterations as well as cell death induced by cisplatin.
Many of the splicing alterations induced by cisplatin are caused by SRSF4 and they contribute to apoptosis in a process requires class I PI3K.
Chemotherapy with platinum-based compounds is used extensively for the treatment of a wide range of solid tumours, including breast cancers resistant to first line therapy, ovarian, non-small cell lung, testis, endometrial, head and neck and colorectal cancers. Cisplatin (cis-diamine platinum (II) dichloride), the founding member of this class of agents, covalently binds to DNA and induces the formation of bulky DNA adducts consisting of intra-strand cross-links preferentially formed between adjacent guanine residues and, to a lower extent, inter-strand DNA lesions [1,2]. Cell toxicity is linked to these adducts that interfere with DNA replication and transcription. Intra-strand cross-links are mainly processed by removal of platinum adducts via the nucleotide excision repair, and inter-strand cross-links are removed via nucleotide excision repair, translesion polymerase and homologous recombination. Cisplatin activates various signalling pathways that include the DNA damage response (DDR) and the PI3K-Akt pathways . The DDR pathway detects and corrects DNA defects. However, when alterations are too numerous or too severe, cells are committed to death and eliminated. The DDR machinery relies on the activity of three enzymes that belong to the phosphatidyl inositol-3 kinases (PI3K) family: DNA-protein kinase (DNA-PK), Ataxia and Telangectasia Mutated (ATM) and Ataxia Telangiectasia and Rad 3-Related (ATR) . These kinases trigger specific and overlapping cascades of signalling events that result in cell cycle arrest, DNA repair or cell death .
Alternative splicing (AS) occurs in more than 90% of multi-exons primary transcripts [4,5]. Proteins produced through AS can have markedly different and sometimes opposite functions, as exemplified by a number of factors involved in apoptosis or cell survival . In other instances, AS controls the level of proteins by producing transcripts carrying premature termination codons that are degraded by non-sense mediated RNA decay (NMD) . Splicing decisions result from an interplay between highly degenerated cis-acting sequences and a large number of trans-acting factors that include the arginine- and serine-rich proteins (SR-proteins) and the heterogenous nuclear ribonucleoproteins (hnRNPs) families . The participation of these factors in splicing control is often regulated by post-translational modifications such as phosphorylation and acetylation which affect their localisation and their interaction with other proteins .
Aberrant AS occurs in cancer and a growing number of studies have reported a functional link between splicing anomalies and the evolution of the disease [9-12]. Several groups, including ours, have shown that chemotherapeutic drugs can affect the AS of a large number of transcripts [13-16]. However, the impact of these changes on the cancer cell is still poorly understood. Here, we analyse the transcriptome of cisplatin-treated cancer cells, and use AS changes to identify pathways that link cisplatin with the cellular response.
Cell culture, authentication, reagent and survival assay
MCF7, MDA-MB-231, HT1080, BT549, RD, HDF1 and HDF2, MG-63, MSU and AT5BIVA (deficient in ATM, Coriell Cell Repository, Camden, NJ, USA) cells were cultured in Dulbecco’s Modified Eagle’s Medium (DMEM, Lonza, Verviers, Belgium) supplemented with non-essential amino-acids (NEAA) (1%), penicillin and streptomycine (1%), gentamycin (0.1%), fungizone (0.1%) and 10% FCS (Lonza). Ishikawa cells (human endometrial adenocarcinoma cell line) were cultured in RPMI 1640-glutamax (Lonza) supplemented with NEAA (1%), sodium pyruvate (1%), penicillin and streptomycine (1%), fungizone (0.1%) and 10% FCS, GM09607 cells (deficient in ATM, Coriell Cell Repository) in EMEM (Lonza) supplemented with 10% FCS and 1% NEAA, and MO59J cells (glioblastoma cell line, deficient in the catalytic subunit of DNA-PK) in DMEM/F12 supplemented as DMEM.
The study conforms to the principles outlined in the Declaration of Helsinki and was approved by the ethic committee of Liège University Hospital (B707201110973).
MCF7 and Ishikawa cells were authenticated by DSMZ (Braunschweig, Germany). Although no authentication of the other cell lines was made, the deficiency in ATM of GM09607 and AT5BIVA was ascertained by western blotting, and that of p53 in MG-63 was confirmed by RT-PCR.
Cisplatin (cis-diamine platinum (II) dichloride), wortmannin, caffein, and triciribine were from Sigma-Aldrich (St-Louis, MO, USA), oxaliplatin from Santa Cruz Biotechnology (Santa Cruz, CA, USA), ATM kinase inhibitor from Calbiochem EMD biosciences (La Jolla, CA, USA), NU7026 from Merck Millipore (Darmstadt, Germany), TGX221, IC87114 and MK2206 from Selleckchem (Munich, Germany) and PX866 from LC Laboratories (Woburn, MA, USA).
Cell survival and apoptosis/necrosis were measured, respectively, by trypan blue exclusion in blind tests and by FACS analysis as described in .
RNA isolation, RT-PCR and RT-qPCR
RNAs were purified from cultured cells using the High Pure RNA isolation kit (Roche, Mannheim, Germany) and quantified by spectrometry. Gene expression was measured by RT-qPCR. Details according to the Minimum Information for Quantitative RT-PCR Experiment (MIQE) guidelines , are given in Additional files 1 and 2. For analysis of exon inclusion/exclusion, primers were chosen on exons surrounding the sequences potentially alternatively spliced. Primers, protocols and amplification products sizes are detailed in Additional file 2. Splice variants were discriminated by electrophoresis as described .
RNA sequencing analysis
RNA libraries and sequencing were performed on total RNA samples at the GIGA Genomics facility, University of Liège, Belgium. The quality of RNA was checked with BioAnalyser 2100 (Agilent technologies, CA, USA) that indicated a RQI score >8. The libraries were prepared with Truseq® mRNA Sample Prep kit (Illumina, CA, USA) from 1 microgram of total RNA following manufacturer’s instructions. mRNAs were isolated by poly-A selection and fragmented (8 minutes at 94°C). Fragmented mRNAs (around 170 nucleotide-long in average) were used for reverse-transcription in the presence of Superscript II (Invitrogen, Oregon, USA) and random primers. After second strand synthesis, end-repair, A-tailing and purification, the double strand cDNA fragments were ligated to Truseq® adapters containing the index sequences. Fifteen cycles of PCR in the presence of dedicated PCR primers and PCR master mix were applied to generate the final libraries. Libraries were sequenced in pair-end sequencing runs on the Illumina GAIIx in multiplexed 2 × 76 base protocols. The raw data was generated through CASAVA 1.6 suite (Illumina, CA, USA). TopHat (http://ccb.jhu.edu/software/tophat/index.shtml) software was used to align RNA-Seq reads to the reference genome (hg19, UCSA) and discover transcript splice sites. Cufflinks (http://cole-trapnelllab.github.io/cufflinks/) used the resulting alignment files to quantify the gene expression levels, identify up- and down-regulated transcripts and find the alternative splice junctions.
SpliceSeq(1.2) (http://bioinformatics.mdanderson.org/main/SpliceSeq:Overview) was used for a focused AS analysis. Using alignment database and Bowtie, SpliceSeq aligns reads from RNA-Seq data to a reference collection of splice variants [19,20].
Lists of genes modulated in term of expression and splicing were imported in the ToppGene Suite for analysis .
Antibodies and Western blotting
Antibodies directed against Akt, phospho-Akt (ser473) and β-actin were purchased from Cell Signalling (Beverly, MA, USA). Cells were lysed in Laemmli buffer containing 50 mM DTT. Lysates were briefly sonicated, incubated at 65°C for 15 min and analyzed by SDS-PAGE. Proteins were electroblotted and detected as described in . Probing of β-actin was performed as a control of protein loading.
SMARTpool siGENOME (Dharmacon by Thermo Fisher Scientific, Lafayette, CO, USA), consisting of four siRNA duplexes, were used to target SRSF4 and SRSF6 mRNA. siRNA targeting ATR were from Ambion (Life technologies). The 5′-UUGCAUACAGGACUCGUUATT-3′ and 5′-UAACGAGUCCUGUAUGCAATT-3′ oligoribonucleotides were used as control siRNA (siSCR) that does not target any known human transcript . Cells were transfected by siRNAs as previously described .
The means and standard deviation were calculated from three or four independent experiments. The significance of differences was determined using t-test or ratio paired t-test of Student.
Cisplatin alters alternative splicing
In vitro treatment of cells by cisplatin induces alterations of splicing in various transcripts [15,24,25]. Following treatment of MCF7 and Ishikawa cell lines with cisplatin (Figure 1A-B), the RT-PCR analysis of MDM2, a negative regulator of p53, showed a reduction in the full length product and the appearance of smaller splicing variants. The smallest variant had the expected size of MDM2-ALT1 splice variant. This splicing shift was maximal at 50 – 100 μM of cisplatin (Figure 1A-B) and after 24 – 48 hours (Figure 1E). A similar dose-dependent shift was obtained with VEGF where cisplatin decreased the expression of VEGF-165 and concomitantly increased the production of the VEGF-111 splice variant (Figure 1C-D). Similar alterations of splicing of MDM2 were also observed in MDA-MB-231 (breast adenosarcoma), BT549 (breast carcinoma), HT1080 (fibrosarcoma), RD (rhabdomyosarcoma), MG-63 (osteosarcoma), MSU (fibrobastic cell line) and HDF1 and 2 (primary dermal fibroblasts) cells treated with 50 μM cisplatin for 24 hours (Figure 1F). The cisplatin analog oxaliplatin induced similar effects on MDM2 splicing, suggesting that this splicing alteration is generalized to platinum-based agents (Figure 1G).
Poly A+ RNA from MCF7 cells untreated or treated with 50 μM cisplatin for 24 hours was isolated and prepared for next-generation sequencing analysis. No significant cell death over untreated samples was noted in these conditions (as measured by trypan blue exclusion). The average number of reads approached or exceeded 20 millions in both samples. Alignment of transcripts to the genome indicated that 16733 and 16969 genes were expressed in the control and cisplatin-treated samples, respectively. Sequencing data are available in the ArrayExpress database (www.ebi.ac.uk/arrayexpress) under accession number E-MTAB-2663. The global gene expression in the two conditions was highly correlated, with a Pearson correlation coefficient = 0.835 (p = 0.000000).
Effect of cisplatin on gene expression
For differential gene expression, the following filters were applied: absolute fold change >2 and q-value < 0.05. Five hundred fifty-three genes were regulated (111 up and 442 down, Additional file 3). The top 20 up-regulated and down-regulated genes are listed in Table 1 with SERPINB5 (126×) and GPHN (159×) being first in each category, respectively. The expression of a panel of genes commonly used as calibrators was not significantly affected (GAPDH: 1.20; ACTB: 0.79; ACTG: 0.93; PPIA (cyclophylin A): 0.94; PPIB (cyclophylin B): 1.07). RT-qPCR was performed to confirm the expression level of 9 up-, down- or non-regulated genes on the samples used for RNA-Seq. Selected genes were either conventional calibrators (GAPDH, ACTB, β2M) or encoded proto-oncogenes (MYB, SERPINB5, JAK2), anti-oncogenes (BRCA1, RB1) and a factor regulating apoptosis (FAS). RT-qPCR analysis correlated with RNA-seq data with a Pearson coefficient of 0.98 (p = 0.000004), validating the RNA-seq data. Moreover, these changes in gene expression noted by RT-qPCR were confirmed in three independent experiments using MCF7 cells, and in two independent experiments in Ishikawa cells, indicating that the changes are reproducible and not restricted to MCF7 cells.
Gene ontology analysis was performed using the ToppFun Suite software. Significantly (p < 0.05) affected biological processes were identified (Table 2). Many genes regulated by cisplatin belong to two main groups: cell cycle and proliferation. Surprisingly, neurogenesis also appeared as a regulated category.
ToppFun identified genes matching annotations for transcription factors PITX2 (38 genes), E2F (18 genes, Table 3) and FOXF2 (19 genes). For genes matching with PITX2 and FOX2, no significant difference in the proportion of up- and down-regulated genes (versus total numbers of up- and down-regulated genes, respectively) was observed. In contrast, the 18 genes matching with E2F were all down-regulated. As the expression of E2F transcription factors themselves was not significantly changed, this suggests that cisplatin may affect their activity.
The list of genes regulated by cisplatin was compared to lists of oncogenes (http://www.uniprot.org/uniprot/?query=keyword:KW-0656) and tumor suppressors (http://www.uniprot.org/uniprot/?query=keyword:KW-0043). Data showed that cisplatin reduced the level of some tumor suppressor genes and of oncogenes while inducing others (Table 4). Strikingly, the expression of the two AP-1 members FOS and JUN was found to be increased.
Effect of cisplatin on post-transcriptional events
Potential modifications of splicing by cisplatin were investigated from the RNA-seq data. The SpliceSeq software identified 717 AS events occurring in 619 primary transcripts (Additional file 4). Only 5 genes (UGDH, SLC38A1, RETSAT, PDE8A, NASP) (0.44%) were affected simultaneously at transcriptional and post-transcriptional levels. Changes in splicing were grouped based on the type of events being affected: 79 changes involved cassette exon inclusion events, 166 were cassette exon exclusion events (of which 49% were not annotated as alternative exons in NCBI), 243 changes affected alternative 5′ or 3′ splice site selection events, 144 involved alternative promoters, 83 indicated alternative terminations and 2 were splicing changes attributed to mutually exclusive exon. Significantly affected biological processes identified by ToppFun Suite software on affected genes were “RNA splicing and processing” and “cell cycle” (Table 2).
For validation purpose, 16 splicing events identified by SpliceSeq as affected by cisplatin were evaluated by RT-PCR. These events were chosen such as to cover the range of AltSplice RPKM (reads per kilobase per million reads) values from 0 to 40. RT-PCR confirmed the alternative splicing of the ten exon skipping events and four of the six exon inclusion events in MCF7 and Ishikawa cells, indicating a good concordance with RefSeq data (drawings of nine splicing events in MCF-7 cells are shown in Figure 2).
PI3K pathway, but not DNA damage response and p53, is involved in the alteration of splicing by cisplatin
As cisplatin induces DNA damage, the contribution of three main actors of the DDR pathway (ATM, ATR and DNA-PK) in cisplatin-induced exon inclusion/exclusion was investigated. Cisplatin affected the splicing in ATM-deficient AT5BIVA cells (illustrated for HNRNPDL exon 6 exclusion and exon 8 inclusion in Figure 3A) and GM9607 cells (not illustrated). Similarly to its effect in MCF7 cells (Figure 2), identical data were also found in MO59J cells that lack the catalytic subunit of DNA-PK (Figure 3B), and in MCF7 cells transfected with a siRNA targeting ATR (not shown). Moreover, specific inhibition of ATM (using ATM kinase inhibitor) or of DNA-PK (using NU7026) did not reverse the splicing induced by cisplatin in MCF7 cells (Figure 3C-D). Combined inhibition experiments were performed to evaluate whether the three DDR members might functionally compensate each other. Inhibition of ATM and ATR (using caffeine) and DNA-PK (using NU7026) failed to reverse the splicing induced by cisplatin (Figure 3E-F). The same results were obtained in MO59J cells (deficient in DNA-PK activity) treated with caffeine (not shown). Together, these data strongly suggest that the DDR does not participate in the splicing change of HNRNPDLe6 and AMZ2 induced by cisplatin. The DNA damage-activated protein p53 is similarly not involved since cisplatin induces an alteration of MDM2 splicing in the p53-deficient cell line MG-63 (Figure 1F).
As cisplatin is also known to activate PI3K in several cell types , the implication of the PI3Ks in the AS changes induced by cisplatin was evaluated. MCF7 cells were pre-treated with wortmannin three hours prior to adding cisplatin. At the concentrations used, wortmannin inhibits class I and III PI3Ks as well as PI3KC2b, but not ATM, ATR, DNA-PK. Dose-dependent inhibition of the effect of cisplatin on the splicing events was observed (Figure 3G-J). A similar impact was observed with the wortmannin derivative PX866 (Figure 3K-L). To gain further insights into the identity of PI3Ks involved, inhibitors specifically targeting the class I PI3Ks (TGX211, IC87114) were used (Figure 4A-D). As observed with wortmannin, these inhibitors significantly reversed the cisplatin-induced splicing changes, while no effect was observed in the absence of cisplatin. These results suggest that the class I PI3Ks are involved in the cisplatin-mediated response. RNA-Seq data indicate that p110α and p110β, but not p110γ and p110δ, are expressed in MCF7 cells. P110α, but not p110β, is activated by insulin. At the concentration used, TGX221 and IC87114 did not reduce the phosphorylation of Akt induced by insulin (not shown). Cisplatin treatment did not induce Akt phosphorylation on Ser473 under our experimental conditions (not shown). Moreover, the Akt inhibitors triciribine and MK2206, while efficiently reducing the insulin-induced phosphorylation of Akt, did not affect the cisplatin-induced changes in splicing in AMZ2 and HNRNPDL-E6 (Figure 4E-H). Finally, insulin did not induce a change in AS that was similar to cisplatin in conditions that increase Akt phosphorylation (not shown). Together, these observations strongly suggest that the splicing alterations elicited by cisplatin require p110β, but are independent of Akt.
Cisplatin-induced alteration of splicing involves SRSF4
Using a siRNA screen targeting 57 splicing factors, we identified SRSF4 as a regulator of hnRNPDL exon 6 splicing in MCF7 cells in basal growth conditions. We evaluated the role of SRSF4 in mediating the effect of cisplatin on AS by using a siRNA targeting SRSF4 and a control siRNA targeting the splicing factor SRSF6. Reduced levels of SRSF4 mRNA (81 ± 7% reduction, n = 3) and SRSF6 mRNA (80 ± 10%, n = 3) were confirmed by end-point RT-PCR (Figure 5A-B). siSRSF4 alone or when combined to siSRSF6 (siSRSF4/6) partly abrogated the splicing changes induced by cisplatin in the events tested (Figure 5C-F). The control siSCR and siSRSF6 alone had no effect. A similar reduction of the cisplatin-induced HNRNPDL exon 6 exclusion and exon 8 inclusion after knock-down of SRSF4 was observed in the breast cancer cell line BT549 (not illustrated).
Cisplatin-induced cell death involves SRSF4
To address the potential involvement of SRSF4-dependent splicing events induced by cisplatin in cell death or cell survival, MCF7 cells transfected with the control siRNA siSCR or with siSRSF4 were treated or not with cisplatin for 48 hours. Apoptosis and necrosis were measured by FACS (Figure 6A-D). Although the downregulation of SRSF4 had no effect on growth when cisplatin was absent, it strongly reduced cell death observed in the presence of cisplatin (15 ± 4% for siSRSF4 versus 6 ± 2% for siSCR; p = 0.02), which represents a 62 ± 9% reduction of cisplatin-induced cell death by SRSF4 repression (Figure 6E). These data were further confirmed by using the trypan blue exclusion assay (53 ± 3% reduction, p ≤ 0.01, n = 3, Figure 6E).
The development of chemotherapeutic agents has enabled tremendous progress in cancer therapy. However, the success of these treatments is offset by the development of drug resistance and by toxic side-effects on healthy cells and tissues. The development of this resistance is encouraged by several processes, including decreased access and increased efflux of the drug from the tumor, altered expression of oncogenes, reduced apoptosis and increased DNA repair . In order to evaluate the role of AS in the efficiency of cisplatin, we performed a transcriptome analysis of breast cancer cell line because platinum-based chemotherapy is used as second and third-line of treatment against resistant metastatic breast cancer [28,29]. Moreover, MCF7 cells are well-characterized notably in terms of their response to chemotherapeutic drugs. Our results indicate that cisplatin affects the expression level (absolute fold change >2) of more than 500 genes and provokes changes in at least 700 splicing events, thereby extending previous observations that chemotherapeutic agents affect AS [6,13,16]. This splicing reprogramming also occurs in other transformed cell lines including the breast cancer cell lines MDA-MB-231 and BT549, the endometrial adenocarcinoma cell line Ishikawa and in primary fibroblasts.
Many of the genes whose expression is altered by cisplatin have functions in cell cycle. Cisplatin-induced changes also affect the expression of tumor suppressor genes, oncogenes and genes involved in determining cell fate (Table 4). Strikingly, the list lacks genes encoding splicing factors, suggesting that the impact on splicing control principally stems from post-transcriptional and/or post-translational events affecting their expression, localization and activity. In contrast, cisplatin affected the AS of many splicing factors. Accordingly, our gene ontology analysis suggests that splicing function may be one of the pathways most affected by cisplatin.
We observed that other chemotherapeutic drugs, namely camptothecin and doxorubicin, induce the same changes in AS as those elicited by cisplatin (unpublished work). As these drugs all induce DNA damage, it is tempting to speculate that activation of the DDR pathway may be involved in promoting these splicing alterations. In contrast to this prediction, the genetic depletion and/or the specific inhibition of p53, ATM, ATR and DNA-PK failed to suppress AS re-programming upon cisplatin treatment. These data contrast with those of Shkreta et al.  who observed that the shift in Bcl-x splicing induced by oxaliplatin or cisplatin in HEK-293 cells was abrogated by inhibiting ATM, ATR or p53. However, no significant change in Bcl-x splicing by cisplatin was recorded here by deep sequencing or RT-PCR in MCF7 cells (not illustrated), consistent with the very small shift previously observed in MCF7 cells . These discrepancies may be related to the different cell lines used, which may display different thresholds to elicit the DNA damage response.
Previous reports indicate that the PI3K/Akt axis can affect the AS of many primary transcripts at least in part by activation of SRPK and the phosphorylation of SR proteins [31-33]. We investigated the role of this pathway in the AS changes induced by cisplatin by using a panel of inhibitors. Our results indicate that cisplatin alters AS in a process that requires the PI3K subunit p110β. The link between p110β and the splicing events altered by cisplatin remains unclear but is independent of Akt. An intriguing possibility is that cisplatin affects the nuclear activity of p110β, which in turn may directly affect the activity of splicing factors. A role for p110β is not totally unexpected since there is mounting evidence indicating that nuclear lipids can regulate nuclear functions including splicing [34,35]. While phosphoinositides associate with nuclear membranes, they also co-localize in nuclear speckles  and interact with various proteins or ribonucleoprotein complexes including the spliceosome components U2 snRNP, U4/U6 snRNP and SF3A1.
We observed that knocking down SRSF4, but not SRSF6, abrogated the cisplatin-induced changes in splicing. CLIP analysis followed by high-throughput sequencing identified GA rich pentamers with G/AAAA/GA sequence as a consensus motif for the binding of SRSF4 to RNA . Moreover, SRSF4 preferentially binds to exons, with a peak of binding ~50 nucleotides upstream of the 5′ splice site. Sequences matching with these sequences are observed in the exons that were skipped in response to cisplatin. However, that they represent binding sites for SRSF4 remains to be tested.
Although SRSF4 may also have an indirect function, for example by regulating the splicing of other splicing factors, we believe that this scenario is unlikely to explain the rapid changes in the steady state levels of splice variants imposed by cisplatin. Nevertheless, portions of the RS-rich regions of SRSF3 and SRSF7 are truncated due to exon skipping (SRSF7) or alternative termination (SRSF3) in response to cisplatin treatment, thereby possibly affecting the phosphorylation of these proteins and their association with other splicing partners.
A link between altered splicing and the efficacy of cancer treatment is suggested by several findings. In lymphocytes of patients with chronic lymphocytic leukemia, mutations in the gene encoding the splicing factor SF3B1 are more frequent after treatment, suggesting a chemotherapy-driven clonal selection for cells being affected in splicing [38,39]. The efficacy of chemotherapetic agents may act at least in part through reprogrammation of AS Consistent with this view, treatments of human 293 cell line with a panel of chemotherapeutic agents induced splicing shifts that encouraged the production of pro-apoptotic variants of Bcl-x, caspase-9 and survivin . Moreover, altering the ratio of splice variants of caspase-9 reduced the resistance of non-small lung cancer cells to various chemotherapeutical agents . On the contrary, splicing switches toward anti-apoptotic versions, as in the conversion from FAS to sFAS, have also been observed ( and personal observation). As high sFAS levels correlate with poor survival in patients with T-cell leukemia and gynecological malignancies [41,42], sFAS may contribute to the acquisition of drug resistance and a chemotherapy designed to revert splicing to FAS may increase treatment efficiency .
GO terms related to apoptosis were not highlighted by hierarchization analysis of the transcripts alternatively spliced upon cisplatin. We compared a list of transcripts related to apoptosis (GSEA [44,45]) with the list of transcripts with splicing affected by cisplatin treatment. Twenty-six actors involved in the regulation of apoptosis were common to both lists, as for example BAX, caspase-6, caspase-8 (pro-apoptotic) and MADD, API5 (anti-apoptotic). These examples illustrate that cisplatin-induced alterations of splicing may have both anti- and pro-apoptotic effects, and the net effect cannot be estimated on a theoretical basis.
Here, we observed that knocking down SRSF4 reduced the impact of cisplatin on cell death, suggesting an overall therapeutic benefit associated with the expression of SRSF4. Thus, while the pharmacological alterations of splicing induced by chemotherapic agents may fuel therapeutic efficiency, preventing these alterations by inhibiting SRSF4-regulated splicing may help cells to resist the cisplatin treatment. This situation is likely to be more complex given the large number of splicing regulators, their combinatorial mode of regulation and the diversity of their targets. A growing list of pharmacological agents that can modulate splicing is now emerging, with some demonstrating anti-tumor activity [46-48]. Pladienolide, spliceostatin and herboxidiene modulate the function of the spliceosome by binding to the SF3B core component protein [49,50]. A link between splicing alterations and inhibition of cancer cell proliferation was established , supporting the concept of using splicing to improve anti-cancer therapy. Another example is provided by the anti-hypertensive agent amiloride that also affects the level and/or the phosphorylation of splicing factors, alters the splicing of cancer genes in various tumor cell lines and sensitizes chronic myelogenous leukemia cells to imatinib . Similarly, dietary agents possessing anticancer activities as curcumin, resveratrol and epigallocatechin-gallate, have been shown to affect splicing, at least in part through modulation of splicing factors levels [52-55].
We showed that the reprogramming of splicing induced by cisplatin makes a large contribution to its anti-cancer property, and that its action requires class I PI3K p110β and the splicing factor SRSF4. In this context, our data have two major implications. They suggest that pharmacologically modulating AS can potentially affect the success of chemotherapy. Moreover, they raise the interesting possibility that molecules or conditions (as drugs used for non-tumoral diseases, food components and redox status) that modify AS may influence the response to anti-cancer treatments.
Wang D, Lippard SJ. Cellular processing of platinum anticancer drugs. Nat Rev Drug Discov. 2005;4(4):307–20. doi:10.1038/nrd1691.
Kartalou M, Essigmann JM. Mechanisms of resistance to cisplatin. Mutat Res. 2001;478(1–2):23–43.
Yang J, Yu Y, Hamrick HE, Duerksen-Hughes PJ. ATM, ATR and DNA-PK: initiators of the cellular genotoxic stress responses. Carcinogenesis. 2003;24(10):1571–80. doi:10.1093/carcin/bgg137.
Blencowe BJ. Alternative splicing: new insights from global analyses. Cell. 2006;126(1):37–47. doi:10.1016/j.cell.2006.06.023.
Burge C, Karlin S. Prediction of complete gene structures in human genomic DNA. J Mol Biol. 1997;268(1):78–94. doi:10.1006/jmbi.1997.0951.
Shkreta L, Froehlich U, Paquet ER, Toutant J, Elela SA, Chabot B. Anticancer drugs affect the alternative splicing of Bcl-x and other human apoptotic genes. Mol Cancer Ther. 2008;7(6):1398–409. doi:10.1158/1535-7163.mct-08-0192.
McGlincy NJ, Smith CW. Alternative splicing resulting in nonsense-mediated mRNA decay: what is the meaning of nonsense? Trends Biochem Sci. 2008;33(8):385–93. doi:10.1016/j.tibs.2008.06.001.
Zhou Z, Fu XD. Regulation of splicing by SR proteins and SR protein-specific kinases. Chromosoma. 2013;122(3):191–207. doi:10.1007/s00412-013-0407-z.
Tang JY, Lee JC, Hou MF, Wang CL, Chen CC, Huang HW, et al. Alternative splicing for diseases, cancers, drugs, and databases. Sci World J. 2013;2013:703568. doi:10.1155/2013/703568.
Venables JP. Unbalanced alternative splicing and its significance in cancer. Bioessays. 2006;28(4):378–86. doi:10.1002/bies.20390.
Shkreta L, Bell B, Revil T, Venables JP, Prinos P, Elela SA, et al. Cancer-associated perturbations in alternative Pre-messenger RNA splicing. Cancer Treat Res. 2013;158:41–94. doi:10.1007/978-3-642-31659-3_3.
Venables JP, Klinck R, Koh C, Gervais-Bird J, Bramard A, Inkel L, et al. Cancer-associated regulation of alternative splicing. Nat Struct Mol Biol. 2009;16(6):670–6. doi:10.1038/nsmb.1608.
Solier S, Barb J, Zeeberg BR, Varma S, Ryan MC, Kohn KW, et al. Genome-wide analysis of novel splice variants induced by topoisomerase I poisoning shows preferential occurrence in genes encoding splicing factors. Cancer Res. 2010;70(20):8055–65. doi:10.1158/0008-5472.can-10-2491.
Mineur P, Colige AC, Deroanne CF, Dubail J, Kesteloot F, Habraken Y, et al. Newly identified biologically active and proteolysis-resistant VEGF-A isoform VEGF111 is induced by genotoxic agents. J Cell Biol. 2007;179(6):1261–73. doi:10.1083/jcb.200703052.
Chandler DS, Singh RK, Caldwell LC, Bitler JL, Lozano G. Genotoxic stress induces coordinately regulated alternative splicing of the p53 modulators MDM2 and MDM4. Cancer Res. 2006;66(19):9502–8. doi:10.1158/0008-5472.can-05-4271.
Dutertre M, Sanchez G, Barbier J, Corcos L, Auboeuf D. The emerging role of pre-messenger RNA splicing in stress responses: sending alternative messages and silent messengers. RNA Biol. 2011;8(5):740–7. doi:10.4161/rna.8.5.16016.
Neutelings T, Lambert CA, Nusgens BV, Colige AC. Effects of mild cold shock (25 degrees C) followed by warming up at 37 degrees C on the cellular stress response. PLoS One. 2013;8(7):e69687. doi:10.1371/journal.pone.0069687.
Bustin SA, Benes V, Garson JA, Hellemans J, Huggett J, Kubista M, et al. The MIQE guidelines: minimum information for publication of quantitative real-time PCR experiments. Clin Chem. 2009;55(4):611–22. doi:10.1373/clinchem.2008.112797.
Langmead B, Trapnell C, Pop M, Salzberg SL. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009;10(3):R25. doi:10.1186/gb-2009-10-3-r25.
Ryan MC, Cleland J, Kim R, Wong WC, Weinstein JN. SpliceSeq: a resource for analysis and visualization of RNA-Seq data on alternative splicing and its functional impacts. Bioinformatics (Oxford, England). 2012;28(18):2385–7. doi:10.1093/bioinformatics/bts452.
Chen J, Bardes EE, Aronow BJ, Jegga AG. ToppGene Suite for gene list enrichment analysis and candidate gene prioritization. Nucleic Acids Res. 2009;37(Web Server issue):W305–11. doi:10.1093/nar/gkp427.
Ho TT, Merajver SD, Lapiere CM, Nusgens BV, Deroanne CF. RhoA-GDP regulates RhoB protein stability. Potential involvement of RhoGDIalpha. J Biol Chem. 2008;283(31):21588–98. doi:10.1074/jbc.M710033200.
Deroanne C, Vouret-Craviari V, Wang B, Pouyssegur J. EphrinA1 inactivates integrin-mediated vascular smooth muscle cell spreading via the Rac/PAK pathway. J Cell Sci. 2003;116(Pt 7):1367–76.
Avery-Kiejda KA, Zhang XD, Adams LJ, Scott RJ, Vojtesek B, Lane DP, et al. Small molecular weight variants of p53 are expressed in human melanoma cells and are induced by the DNA-damaging agent cisplatin. Clin Cancer Res. 2008;14(6):1659–68. doi:10.1158/1078-0432.ccr-07-1422.
Edmond V, Moysan E, Khochbin S, Matthias P, Brambilla C, Brambilla E, et al. Acetylation and phosphorylation of SRSF2 control cell fate decision in response to cisplatin. EMBO J. 2011;30(3):510–23. doi:10.1038/emboj.2010.333.
Muscella A, Urso L, Calabriso N, Ciccarese A, Migoni D, Fanizzi FP, et al. Differential response of normal, dedifferentiated and transformed thyroid cell lines to cisplatin treatment. Biochem Pharmacol. 2005;71(1–2):50–60. doi:10.1016/j.bcp.2005.10.022.
Lu HP, Chao CC. Cancer cells acquire resistance to anticancer drugs: an update. Biomed J. 2012;35(6):464–72. doi:10.4103/2319-4170.104411.
Boulikas T, Vougiouka M. Recent clinical trials using cisplatin, carboplatin and their combination chemotherapy drugs (review). Oncol Rep. 2004;11(3):559–95.
Cobleigh MA. Other options in the treatment of advanced breast cancer. Semin Oncol. 2011;38 Suppl 2:S11–6. doi:10.1053/j.seminoncol.2011.04.005.
Shkreta L, Michelle L, Toutant J, Tremblay ML, Chabot B. The DNA damage response pathway regulates the alternative splicing of the apoptotic mediator Bcl-x. J Biol Chem. 2011;286(1):331–40. doi:10.1074/jbc.M110.162644.
Blaustein M, Pelisch F, Srebrow A. Signals, pathways and splicing regulation. Int J Biochem Cell Biol. 2007;39(11):2031–48. doi:10.1016/j.biocel.2007.04.004.
Eisenreich A, Malz R, Pepke W, Ayral Y, Poller W, Schultheiss HP, et al. Role of the phosphatidylinositol 3-kinase/protein kinase B pathway in regulating alternative splicing of tissue factor mRNA in human endothelial cells. Circ J. 2009;73(9):1746–52.
Zhou Z, Qiu J, Liu W, Zhou Y, Plocinik RM, Li H, et al. The Akt-SRPK-SR axis constitutes a major pathway in transducing EGF signaling to regulate alternative splicing in the nucleus. Mol Cell. 2012;47(3):422–33. doi:10.1016/j.molcel.2012.05.014.
Okada M, Ye K. Nuclear phosphoinositide signaling regulates messenger RNA export. RNA Biol. 2009;6(1):12–6.
Martelli AM, Ognibene A, Buontempo F, Fini M, Bressanin D, Goto K, et al. Nuclear phosphoinositides and their roles in cell biology and disease. Crit Rev Biochem Mol Biol. 2011;46(5):436–57. doi:10.3109/10409238.2011.609530.
Chen R, Kang VH, Chen J, Shope JC, Torabinejad J, DeWald DB, et al. A monoclonal antibody to visualize PtdIns(3,4,5)P(3) in cells. J Histochem Cytochem. 2002;50(5):697–708.
Anko ML, Muller-McNicoll M, Brandl H, Curk T, Gorup C, Henry I, et al. The RNA-binding landscapes of two SR proteins reveal unique functions and binding to diverse RNA classes. Genome Biol. 2012;13(3):R17. doi:10.1186/gb-2012-13-3-r17.
Landau DA, Carter SL, Stojanov P, McKenna A, Stevenson K, Lawrence MS, et al. Evolution and impact of subclonal mutations in chronic lymphocytic leukemia. Cell. 2013;152(4):714–26. doi:10.1016/j.cell.2013.01.019.
Rozovski U, Keating M, Estrov Z. The significance of spliceosome mutations in chronic lymphocytic leukemia. Leuk Lymphoma. 2013;54(7):1364–6. doi:10.3109/10428194.2012.742528.
Shultz JC, Goehe RW, Murudkar CS, Wijesinghe DS, Mayton EK, Massiello A, et al. SRSF1 regulates the alternative splicing of caspase 9 via a novel intronic splicing enhancer affecting the chemotherapeutic sensitivity of non-small cell lung cancer cells. Mol Cancer Res. 2011;9(7):889–900. doi:10.1158/1541-7786.mcr-11-0061.
Kamihira S, Yamada Y, Tomonaga M, Sugahara K, Tsuruda K. Discrepant expression of membrane and soluble isoforms of Fas (CD95/APO-1) in adult T-cell leukaemia: soluble Fas isoform is an independent risk factor for prognosis. Br J Haematol. 1999;107(4):851–60.
Konno R, Takano T, Sato S, Yajima A. Serum soluble fas level as a prognostic factor in patients with gynecological malignancies. Clin Cancer Res. 2000;6(9):3576–80.
Eblen ST. Regulation of chemoresistance via alternative messenger RNA splicing. Biochem Pharmacol. 2012;83(8):1063–72. doi:10.1016/j.bcp.2011.12.041.
Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 2005;102(43):15545–50. doi:10.1073/pnas.0506580102.
Mootha VK, Lindgren CM, Eriksson KF, Subramanian A, Sihag S, Lehar J, et al. PGC-1alpha-responsive genes involved in oxidative phosphorylation are coordinately downregulated in human diabetes. Nat Genet. 2003;34(3):267–73. doi:10.1038/ng1180.
Bonnal S, Vigevani L, Valcarcel J. The spliceosome as a target of novel antitumour drugs. Nat Rev Drug Discov. 2012;11(11):847–59. doi:10.1038/nrd3823.
Webb TR, Joyner AS, Potter PM. The development and application of small molecule modulators of SF3b as therapeutic agents for cancer. Drug Discov Today. 2013;18(1–2):43–9. doi:10.1016/j.drudis.2012.07.013.
Butler MS. Remediating cancer via splicing modulation. J Med Chem. 2013;56(17):6573–5. doi:10.1021/jm401289z.
Gao Y, Vogt A, Forsyth CJ, Koide K. Comparison of splicing factor 3b inhibitors in human cells. Chembiochem. 2013;14(1):49–52. doi:10.1002/cbic.201200558.
Yokoi A, Kotake Y, Takahashi K, Kadowaki T, Matsumoto Y, Minoshima Y, et al. Biological validation that SF3b is a target of the antitumor macrolide pladienolide. FEBS J. 2011;278(24):4870–80. doi:10.1111/j.1742-4658.2011.08387.x.
Chang WH, Liu TC, Yang WK, Lee CC, Lin YH, Chen TY, et al. Amiloride modulates alternative splicing in leukemic cells and resensitizes Bcr-AblT315I mutant cells to imatinib. Cancer Res. 2011;71(2):383–92. doi:10.1158/0008-5472.can-10-1037.
Fang HY, Chen SB, Guo DJ, Pan SY, Yu ZL. Proteomic identification of differentially expressed proteins in curcumin-treated MCF-7 cells. Phytomedicine. 2011;18(8–9):697–703. doi:10.1016/j.phymed.2010.11.012.
Markus MA, Marques FZ, Morris BJ. Resveratrol, by modulating RNA processing factor levels, can influence the alternative splicing of pre-mRNAs. PLoS One. 2011;6(12):e28926. doi:10.1371/journal.pone.0028926.
Kim MH. Protein phosphatase 1 activation and alternative splicing of Bcl-X and Mcl-1 by EGCG + ibuprofen. J Cell Biochem. 2008;104(4):1491–9. doi:10.1002/jcb.21725.
Jiang M, Huang O, Zhang X, Xie Z, Shen A, Liu H, et al. Curcumin induces cell death and restores tamoxifen sensitivity in the antiestrogen-resistant breast cancer cell lines MCF-7/LCC2 and MCF-7/LCC9. Molecules (Basel, Switzerland). 2013;18(1):701–20. doi:10.3390/molecules18010701.
We thank Dr. Michael C. Ryan (In Silico Solutions, Fairfax, VA) for the help with the SpliceSeq® software. The assistance of Raafat Stefan (Cell Imaging and Flow cytometry platform, GIGA, University of Liège, Belgium) with FACS analysis is acknowledged. We thank Pr. Betty Nusgens for her careful review of the manuscript. B.C. is the Canada Research Chair in Functional Genomics. This work was supported by a grant from the National Fund for Scientific Research, Belgium (F.N.R.S-Télévie, # 7.4634.10), from Belspo, from Research Concerted Action (# ARC 10/15-02) and the Fonds Léon Frédericq of the University of Liège, Belgium, and from the Canadian Institutes of Health Research (MOP136948 to B.C.).
The authors declare that they have no competing interests.
MG performed most experiments and data analysis, drafted and edited the manuscript. YD participated in the treatment of raw data of RNA sequencing and sequence alignment. RK and BC performed the siRNA screen. AD performed part of the experiments. YH and JP helped to design the experiments and to draft the manuscript. BH carried out the RNA sequencing analysis. AC and CL conceived and coordinated the studies, designed the experiments, and drafted the manuscript. All authors read, edited and approved the final manuscript.
Alain Colige and Charles Lambert contributed equally to this work.
Detailed real-time RT-qPCR procedure, according to the MIQE guidelines.
Sequences of the primers used for RT-qPCR and RT-PCR analyses.
Differential gene expression of controls and cisplatin-treated MCF7 cells.
Post-transcriptional events, from RNA-Seq in control and cisplatin-treated MCF7 cells.