Abstract
Reprogramming somatic cells into induced pluripotent stem cells (iPSCs) has provided huge insight into the pathways, mechanisms and transcription factors that control differentiation. Here we use high-throughput RT–PCR technology to take a snapshot of splicing changes in the full spectrum of high- and low-expressed genes during induction of fibroblasts, from several donors, into iPSCs and their subsequent redifferentiation. We uncover a programme of concerted alternative splicing changes involved in late mesoderm differentiation and controlled by key splicing regulators MBNL1 and RBFOX2. These critical splicing adjustments arise early in vertebrate evolution and remain fixed in at least 10 genes (including PLOD2, CLSTN1, ATP2A1, PALM, ITGA6, KIF13A, FMNL3, PPIP5K1, MARK2 and FNIP1), implying that vertebrates require alternative splicing to fully implement the instructions of transcriptional control networks.
Similar content being viewed by others
Introduction
Intense efforts are dedicated to the study of stem cells because of the potential medical benefits and the insights into differentiation that they offer1,2,3,4. Pluripotent stem cells can be induced into ectoderm, endoderm and mesoderm, recapitulating the early events in embryogenesis5,6. Early differentiation involves the induction of cell line-specific transcription factors, whereas the late stages of differentiation are less well characterized4,7,8,9. Alternative splicing is a means by which the genome can control the expression of different protein isoforms that can drive fundamental cellular changes and this is mostly achieved through the differential expression of RNA-binding proteins10,11,12. Ever since the discovery that alternative splicing controls Drosophila sex differentiation over 20 years ago, alternative splicing has been suggested to be a potential mechanism whereby different protein isoforms could be delivered to implement mammalian differentiation13,14,15,16.
Several studies have also hinted at an involvement of alternative splicing in pluripotency and pluripotent stem cell differentiation along various pathways17,18. For example, different isoforms of the forkhead transcription factor Foxp1 have differential effects on the induction of key pluripotency genes such as Oct4 and Nanog19. Similarly, alternative splice forms of DNMT3B are specific to stem cells, implying that layered and integrated regulation of gene expression occurs at the levels of transcription and splicing20. Current analyses of genome-wide changes in alternative splicing rely mostly on two techniques: microarrays and next-generation sequencing. However, these techniques are heavily biased towards the most highly expressed genes and cannot quantify alternative splicing of medium- and low-expressed genes that can be expressed at 10,000 times lower levels than housekeeping genes21. We recently developed a unique high-throughput RT–PCR platform to probe the entire repertoire of over 3,000 alternative splice events annotated in the RefSeq database without bias against medium- and low-expressed genes22, and further used this technique to probe the involvement of 81 different RNA-binding proteins in epithelial-to-mesenchymal transition by comparing the splicing profiles of epithelial and mesenchymal tissues with those of the cells lines treated with siRNAs against the 81 proteins23.
Here we have used high-throughput RT–PCR to implicate the RNA-binding protein Muscleblind (MBNL1) in differentiation of stem cells, and we show that essential induction of MBNL1 late during mesoderm differentiation controls a splicing programme that induces and represses specific exons at the RNA level.
Results
Alternative splicing programs in fibroblasts and stem cells
To analyse the programme of alternative splicing changes occurring in human pluripotent stem cells during stem cell induction, maintenance and differentiation, we designed a large-scale screening strategy. Using high-throughput PCR, we previously screened the entire NCBI Reference Sequence (RefSeq) database for the presence or absence of both of two predicted isoforms at thousands of loci. From these, detected, alternatively spliced events we selected 303 high-quality PCRs from across the alternative exons, using only purity and yield of the two PCR products in 10 different cell lines as parameters We used these high-quality PCRs to study alternative splicing in pluripotent stem cells in more detail because they occur in an essentially random and functionally diverse cross-section of genes22. Alternative splicing was studied in human skin fibroblasts and their induced pluripotent stem cells (iPSCs) and then also in fibroblasts redifferentiated from those iPSCs24. We observed a reversible programme of alternative splicing (Fig. 1a correlation R=−0.88). Equally remarkably, the inclusion levels of 15 of the 303 alternative exons shifted >50 percentage points on average between pluripotent stem cells and their parental and daughter fibroblasts (and >40% in both transitions Supplementary Data 1). Such extreme shifts in alternative splicing are rarely described and are defined as ‘switch-like’25,26,27. Thus, using a randomly selected set of alternative splicing events we have uncovered an unprecedented programme of reversible regulated switch-like alternative splicing in stem cells, whereby 5% of a functionally random snapshot of all alternative splicing events can be induced to shift reversibly from predominantly one isoform to nearly exclusively the other, between the fibroblastic and pluripotent states. Furthermore, there is a very good correlation (R=0.824, Supplementary Fig. S2) between percent-spliced-in values obtained with RNA-Seq data28 and our PCR-based methodology for those events in which RNA-seq data had enough read depth to produce confident PSI estimates (104 out of 303 mapped events; Supplementary Data 1). To establish if these alternative splicing changes were a general feature of stem cells we performed manual end-point PCRs with the 15 switch-like splicing events on three different pluripotent stem cells and their resulting differentiated and parental fibroblasts. Samples included the well-characterized H9 embryonic stem cells and their redifferentiated fibroblasts. In all cases we observed a near complete and fully reversible shift from one isoform to the other between stem cells and fibroblasts (Fig. 1b).
Modulation of alternative splicing programs
To create a tool to understand the systems biology of splicing, we recently knocked down 81 potential splicing factors in various cell lines and performed RT–PCRs across 47 alternatively spliced regions in different genes23. Here we monitored these 47 splicing events using RNA extracted from iPSCs and their parental fibroblasts24 with the aim of identifying the splicing factors involved in pluripotency and reprogramming (Supplementary Data 2). The splicing profile of these 47 PCRs showed an equivalent near-perfect anti-correlation of R=−0.91 in reversible stem cell induction and redifferentiation, demonstrating that this smaller set of 47 alternative splice events can be used to investigate the splicing control mechanisms involved in stemness and maintenance of pluripotency (Fig. 1c). The differences in percent-spliced-in values, in these 47 alternative splice events, induced in pluripotency (iPSC values-original fibroblasts) were then compared with the differences induced when we knocked down 81 splicing potential factors (knockdown-control, Supplementary Data 3 and 4). Strikingly, the profile of MBNL1 knockdown correlated most strongly (R=0.6) with the splicing profile of the induction of pluripotency (Fig. 2a, lines 1 and 2).
MBNL1 is a key gene involved in the myotonic dystrophies, DM1 and DM2, which were the first defined RNA diseases. One normal function of Muscleblind is to modulate splicing during muscle and heart development29,30,31,32,33,34. The second most correlated knockdown was RBFOX2, which is a well-characterized splicing factor that we and others have implicated in epithelial-to-mesenchymal transition and invasion22,23,34,35. To confirm the association of these two proteins with stem cell differentiation, we knocked down MBNL1 and RBFOX2 in fibroblasts, alone or in combination, and performed automated RT–PCR for the 303 original alternative splicing events. We observed a significant correlation of splicing changes between individual MBNL1 and RBFOX2 knockdowns and the changes occurring upon induction of the original fibroblasts into the pluripotent state (R=0.4). Remarkably, the correlation was even higher (R=0.6) between splicing changes associated with pluripotent stem cell derivation and those occurring in fibroblasts and a double-knockdown of MBNL1 and RBFOX2 (Fig. 2b). Manual PCR for the 15 switch-like splice events further confirmed these results, with five of the exons (PLOD2, OSBPL3, CTTN, CLSTN and SYNE2) strongly enhanced by RBFOX2 (that is, inclusion of these exons was strongly inhibited by knockdown of RBFOX2), whereas the remaining 10 exons were either enhanced or inhibited by MBNL1 (Fig. 3a,b).
Control of alternative splicing programs by MBNL1 and RBFOX2
One of the major questions in developmental biology is to understand the events related to the specification of the different cell lineages in the embryo. As pluripotent stem cell differentiation is considered to occur in two phases (before and after commitment)4, we sought to establish whether our splicing programme was associated with specification to committed stem cells or to the subsequent differentiation step. We studied alternative splicing of the 15 switch-like targets in H9 human embryonic stem cells (hESCs) that differentiated into early ectoderm, endoderm and mesoderm precursors and surprisingly found only very minor changes in alternative splicing of the 5 RBFOX2-controlled exons (first 5 lines of Supplementary Fig. S1a) and no change at all for the 10 MBNL1-controlled exons. We confirmed that the pluripotent stem cells had indeed committed to the three germ layers using two marker genes for each type (Supplementary Fig. S1b, primers listed in Supplementary Table S1). To confirm that the MBNL- and RBFOX-controlled splicing programme was a late event in differentiation we took samples of stem cells at 1, 2 and 3 weeks into the 5-week differentiation protocol. qPCR analysis demonstrated an early upregulation of RBFOX2 and a later upregulation of MBNL1 (Fig. 3c). Consistent with this, manual end-point PCR demonstrated a gradual change in the five RBFOX2 targets, starting in the first week of differentiation, whereas the 10 MBNL1 targets changed later and more sharply, mostly between 2–3 weeks of differentiation (Fig. 3d). Changes in splicing were essentially complete after 3 weeks of differentiation consistent with those changes occurring during the later stages of mesoderm differentiation.
Next we checked whether MBNL1 and RBFOX2 were indeed downregulated in pluripotent stem cells and we observed almost complete abrogation of MBNL1 expression and significant downregulation of RBFOX2 at both the RNA and protein levels in all stem cell types, irrespective of their provenance (Fig. 4; primers listed in Supplementary Table S1). Downregulation of MBNL and RBFOX genes was also confirmed by retrospective analysis of our previously published microarray data for all three MBNL and RBFOX orthologues24. Both MBNL1 and MBNL2 mRNAs were downregulated in pluripotent stem cells with MBNL1 being the more abundant in fibroblasts, whereas RBFOX2 was essentially the only RBFOX gene present in these cell types (Supplementary Fig. S3a). We also confirmed upregulation of the fibroblast-specific exons of the RBFOX2-enhanced exon of PLOD2 and the MBNL1-enhanced exon of PALM by raising specific antibodies against the alternatively spliced peptides (Supplementary Fig. S3b). We conclude that committed stem cells undergo a massive reversible reprogramming of splicing from inclusion or exclusion of exons to the opposite splice isoform during differentiation (>5% of exons assayed) and that the Muscleblind and Fox proteins (mostly MBNL1 and RBFOX2) control the vast majority of this programme.
MBNL1 and RBFOX2 alternative splicing programs are conserved
To confirm the specificity of alternative splicing changes for late mesoderm differentiation we compared splicing changes for the 15 genes in mouse and human stem cells that we fully differentiated into late ectoderm, endoderm and mesoderm5. In this case we noted that 12 of the alternative splicing events (all excluding OSBPL3, CTTN and CD47) shifted in human and mouse mesoderm and very few of these events changed in ectoderm or endoderm (Fig. 5). Control markers for ectoderm, endoderm and mesoderm differentiation were verified by qPCR (Supplementary Fig. S4). We conclude that these 12 alternative splice events undergo a conserved programme that qualitatively regulates mesodermal gene expression. These results led us to suspect that co-regulation of MBNL1 and RBFOX2 was necessary for the pluripotent stem cell splicing signature and for their differentiation. To test this we inhibited their re-expression in pluripotent stem cells after 2 weeks of redifferentiation. Although we did not observe any effect on the expression of specific markers of pluripotency OCT4, NANOG and SOX2 (Supplementary Fig. S5a), we confirmed that splicing changes occurring late during differentiation can be reversed, either by individual knockdown of RBFOX2 and MBNL1 or by double-knockdown of both splicing factors (Fig. 6). The knockdowns in these samples were consistent with an early induction of RBFOX2 and a late induction of MBNL1.
Examination of the potential physiological roles of the 12 switch-like exons in mammalian mesoderm revealed that nearly all of them encode parts of proteins involved in membrane dynamics, cell adhesion, migration and polarity (Fig. 7a). To support the notion that the observed alternative isoforms have important roles in controlling cellular morphology, we examined their conservation further across evolution. All the genes probably fulfil generic cell functions as they were already present in early metazoans or coelomates. Equally strikingly we found that 8 of the 12 alternatively spliced exons arose concurrently with gene duplications, mostly between jawless and jawed vertebrates (Fig. 7a). For the majority of genes, we found the presence of alternative exons associated with alternative splicing after duplications, whereas no exon homologue could be detected before duplications, suggesting that exon gain and alternative splicing were innovative features selected together to control development in jawed vertebrates. Similar findings on exonic splicing elements and duplications were drawn from genome-wide studies36. Consistent with the function of these conserved splicing events in membrane dynamics, cell adhesion, migration and polarity, we observed a halving of cell size when either MBNL1 or RBFOX2 re-expression was inhibited during stem cell redifferentiation, countering the normal growth in cell size during differentiation (Fig. 7b). Cells had less stress fibres with more actin bundles (Supplementary Fig. S5b). Cells also had reduced level of N-CADHERIN and interestingly increased level of E-CADHERIN, an epithelial marker (Fig. 7c). Consistent with E-CADHERIN re-expression, the RBFOX2 and MBNL1 knockdowns trigger a reduction of the transcription factor TWIST, which represses the expression of E-CADHERIN37 (Supplementary Fig. S5c). We conclude therefore that cooperation between MBNL1 and RBFOX2 is essential for committed stem cell differentiation.
Discussion
Gene expression programs are known to change during germ cell differentiation. In this study we implicated a programme of genetic changes, which occur at the level of alternative splicing of genes that is essential for differentiation into one of the primary gene layers. The genes concerned are involved in important cellular functions including membrane dynamics, cell adhesion, migration and polarity (Fig. 7a). Of the stem cell switch-like targets we identified, several are also involved in diseases, which the isoforms described here likely impact. MBNL1 itself is responsible for muscular dystrophy, SYNE2 and ATP2A1 are involved in the muscle diseases Emery–Dreyfuss muscular dystrophy (OMIM 612999) and Brody myopathy (OMIM 601003) respectively, ITGA6 is implicated in epidermolysis bullosa junctionalis (OMIM 226730) and PLOD2 is involved in Bruck syndrome (OMIM 609220). Most of the target genes described here are novel targets of the MBNL1 splicing factor. Exceptions include ATP2A1 exon 22 splicing, which was previously shown to be controlled by MBNL1 (refs 38, 39), as was the auto-regulation of MBNL1 exon 5 (ref. 40).
The high-throughput PCR method we used identifies many robust splicing shifts starting from all the alternative splice events in the highly curated NCBI Reference Sequence database41. As the PCR method used (35 cycles) amplifies all expressed genes, even low- and medium-expressed ones, in the precise regions of alternative splicing, it creates lists that are complementary to those generated by RNA-Seq and our findings that stem cell splicing mirrors knockdown of MBNL1 have also been observed by this method28.
Interestingly, although genes differentially expressed between iPSCs and fibroblasts are enriched in the same physiological pathways42, genes regulated at the expression level did not overlap with any of the switch-like target exons identified here. Thus, we propose that regulated splicing cooperates with up- and downregulation of genes43 to cement the stem cell-induction during reprogramming and differentiation. Both MBNL1 (refs 29, 44, 45) and RBFOX2 (refs 22, 46, 47) have been shown to act to enhance exons when they bind downstream and inhibit exons when they bind upstream. Analysis of previously published CLIP-seq data29 showed that 4 out of 6 (67%) MBNL1 targets identified here with expression in mouse C2C12 cells (MARK2, FMNL3, CD47 and MBNL1 itself) have evidence of MBNL1-binding clusters, either in or directly upstream of the alternative exon, consistent with a role for MBNL1 downregulation in enhancing these exons in stem cells (Supplementary Data 1). It is therefore likely that binding sites for these proteins have been selected soon after duplications, mostly between jawless and jawed vertebrates to control exons important for cell differentiation in vertebrates. In summary, our study reinforces the idea that splicing factors determine the specific nature of a tissue and that the alternative splicing process has been selected during evolution to control the later stages of differentiation following cell cycle control and commitment.
Methods
Cell culture and in vitro differentiation assays
Human neonatal HFFSM3, line 1 and line 2 (ref. 24) fibroblasts (kind gift from J.M. Lemaitre, IGF Montpellier, France) were maintained in Dulbecco’s Modified Eagle Medium (DMEM, Life Technologies) containing 10% heat-inactivated fetal bovine serum (FBS, PAA), 2 mM L-Glutamine, 1% penicillin and streptomycin (all from Life Technologies). H9 hESCs (WiCell Research) and iPSCs24 (kind gift from J.M. Lemaitre, IGF Montpellier, France) were maintained on matrigel (BD Biosciences) with chemically defined mTeSR (Stemcell Technologies) medium. E14 mouse Embryonic Stem Cells (kind gift from T. Bouschet, IGF Montpellier, France) were maintained on feeder-free culture on gelatin (Sigma-Aldrich) with chemically defined ESGRO complete PLUS clonal grade medium (Millipore).
In vitro differentiation
For the three germ layer differentiation inductions, H9 hESCs were used and differentiated using the protocol described below and in a previous study5.
Ectoderm model: H9 hESCs were seeded at 300,000 cells per well in six-well plates pre-coated with matrigel (BD Biosciences) with chemically defined mTeSR (Stemcell Technologies) medium for 2 days. From days 3 to 9, H9 cells were cultured in the presence of N2 media (Life Technologies) and specific inhibitor of SMAD pathway, SB431542 (TGFbeta receptor inhibitor) and Noggin (Life Technologies). The medium was changed daily from days 3 to 5 and every 2 days from day 5 to day 9.
Endoderm model: H9 hESCs were seeded at 100,000 cells per well in six-well plates pre-coated with matrigel (BD Biosciences) with chemically defined mTeSR (Stemcell Technologies) medium for 1 day. Differentiation was carried out in endoderm priming medium (DMEM 5.5 mM glucose, 10% KO-SR, 2 mM glutamine, 0.1 mM non essential amino-acid, 0,1 mM β-mercaptoethanol) supplemented with Activin A (100 ng ml−1, Life Technologies), LY296002 (30 μM, Sigma-Aldrish) and FGF2 (10 ng ml−1, Peprotech) from days 2 to 4, and Wnt3a (25 ng ml−1, Life Technologies) from days 3 to 4. Then, cells were treated in hepatocyte-priming medium (DMEM 5.5 mM glucose, 2% KO-SR, 2 mM glutamine, 0.1 mM non essential amino-acid, 0,1 mM β-mercaptoethanol), supplemented with BMP4 (50 ng ml−1, Life Technologies), FGF4 (25 ng ml−1, Life Technologies) and 1% DMSO from days 5 to 9. Medium was changed daily in both step of differentiation.
Mesoderm model: H9 hESCs were seeded at 100,000 cells cm−2 in six-well plates pre-coated with matrigel (BD Biosciences) with chemically defined mTeSR (Stemcell Technologies) medium for 1 day. Cells were cultured in DMEM or F12 medium supplemented with 10% FBS (Hyclone Laboratories), supplemented with VEGF (20 ng ml−1, Life Technologies). The medium was changed every 2 days.
For early differentiation assay, differentiated cells were stopped at 10 days. For late differentiation assay, differentiated cells were re-seeded at 10 days and cultured in the same conditions for 10 more days.
For fibroblast differentiation, we used mesoderm model during 3 weeks during which cells were re-seeded once a week. After 3 weeks, VEGF was removed and derivated fibroblasts were cultured as describe above.
The same protocols were used for mouse differentiation assays.
Knockdown of RBFOX2 and MBNL1
HFFSM3 fibroblasts and fibroblast-differentiated hESC H9 cells were transfected with siRNA ON-TARGET plus SMART pools (Dharmacon, RBFOX2: L-020616-01-0005, MBNL1 L-014136-00-005), and scrambled siRNA (Dharmacon, D-001810-10-05). HFFSM3 fibroblasts were transfected using INTERFERin (Polyplus Transfection). Fibroblast-differentiated H9 cells were transfected twice every 2 days using DharmaFECT (Thermo Scientific).
RNA preparation and analysis
RNA was isolated from cells using TRI reagent (Sigma-Aldrich). cDNA was generated by reverse transcription with Maxima Reverse Transcriptase (Thermo Scientific). For qPCR assays, gene expression analysis was performed using 5 ng of cDNA and SYBR Green I Master (Roche). mRNA levels were normalized using the mean expression level of GAPDH, RPLO, RPL13A and PBGD. End-point RT–PCR was performed as in ref. 9. Briefly, reactions were carried out in the Eppendorf Mastercycler PCR Cycler. A first cycle of 15 min at 95 °C was followed by 35 cycles of 30 s at 94 °C, 30 s at 55 °C, and 1 min at 72 °C. The reaction was ended with the extension step of 10 min at 72 °C. Visualization and analysis of amplified products were done using the LabChip HT DNA assay on an automated microfluidic station (Caliper, Hopkinton, MA, USA). For the 15 RBFOX2/MBNL1 selected targets, the amplified products were loaded on 1.5% agarose gel. All primers are detailed in Supplementary Table S1 and Supplementary Data 1.
Immunofluorescence analysis
Cells plated on cover slips were fixed for 10 min in 3.7% formalin (in PBS) followed by a 2-min permeabilization with 0.1% Triton X-100 (in PBS) and incubation in PBS containing 0.1% bovine serum albumin. F-actin was revealed using rhodamin-conjugated phalloidin (Sigma-Aldrich) and nuclei were stained using Hoechst 33342 (Sigma-Aldrich). Cells were washed in PBS and water, mounted with DAKO mounting medium and observed under the fluorescence microscope.
Cell imaging was performed with a Leica DM6000 (Leica, Wetzlae, Germany) with PL APO grade oil × 20 objective. Images were captured with a Coolsnap HQ2 camera (Roper Scientific Inc.) driven by Metamorph (Molecular Devices) and processed using FiJi software and cell surface area was determined using the measure module of FiJi.
Immunoblot analysis
For whole-cell extracts, cells were lysed in RIPA buffer or urea buffer depending on the antibody used (RIPA extract for MBNL1, RBFOX2, α-TUBULIN; urea extract for PLOD2, PALM, actin). For RIPA extracts, cells were lysed in RIPA buffer (50 mM Tris–HCl pH 8, 1% (w/v) NP-40, 1% (w/v) SDS, 0.5% (w/v) Na-deoxycholate (DOC), 2 mM EDTA) supplemented with ‘protease inhibitor cocktail’ (Sigma), incubated for 30 min at 4 °C, centrifuged for 30 min at 14,000 g and the protein amounts were quantifed using Bradford method. For urea extracts, cells were lysed in urea buffer (63 mM Tris–HCl pH7.5, 2% (w/v) SDS, 5% (v/v) 2-mercaptoethanol, 8 M Urea), supplemented with ‘protease inhibitor cocktail’. Cell lysates were next sonicated and centifuged. Supernatants were resolved by SDS–PAGE and transferred to PVDF membranes (Amersham Biosciences). Membranes were blocked in PBS containing 5% (w/v) milk powder and 0.05% (v/v) Tween-20 (PBST) for 1 h and incubated with primary antibodies at 4 °C overnight. After washing with PBST, membranes were incubated with either anti-mouse or anti-rabbit secondary antibodies conjugated to horseradish peroxidase (HRP) for 1 h at room temperature. Signals were detected with ‘chemiluminescent’ (Thermo Scientific), anti-RBFOX2 (1:100, ab51361, Abcam), anti-MBNL1 (gift from Wolfson Centre for Inherited Neuromuscular Disease MB1a-Mouse Anti-MBNL1 (ref. 48)), anti-β-actin (1:20,000, AC-74 clone, Sigma), anti-α tubulin (1:15,000, DM1A clone Sigma), anti-mouse IgG-HRP (1:10,000, 0867428, MP Biomedicals) or anti-rabbit IgG-HRP (1:10,000, 0861202, MP Biomedicals). Polyclonal antibodies to human PLOD2 and PALM alternative peptides were generated by CovalAb, Lyon-France by injecting rabbits with synthetic peptides spanning amino acids 446-TLQREKDSPTPETFQMLSPPK-466 and 175-TVEKDKVTGETRVLSST-191, respectively. Antibodies were affinity-purified using a Sepharose bead column conjugated to the relevant peptide.
Additional information
How to cite this article: Venables, J. P. et al. MBNL1 and RBFOX2 cooperate to establish a splicing programme involved in pluripotent stem cell differentiation. Nat. Commun. 4:2480 doi: 10.1038/ncomms3480 (2013).
References
Park, I.-H. et al. Disease-specific induced pluripotent stem cells. Cell 134, 877–886 (2008).
Robinton, D. A. & Daley, G. Q. The promise of induced pluripotent stem cells in research and therapy. Nature 481, 295–305 (2012).
Wu, S. M. & Hochedlinger, K. Harnessing the potential of induced pluripotent stem cells for regenerative medicine. Nat. Cell Biol. 13, 497–505 (2011).
Evseenko, D. et al. Mapping the first stages of mesoderm commitment during differentiation of human embryonic stem cells. Proc. Natl Acad. Sci. USA 107, 13742–13747 (2010).
Ramirez, J.-M. et al. Brief report: benchmarking human pluripotent stem cell markers during differentiation into the three germ layers unveils a striking heterogeneity: all markers are not equal. Stem Cells 29, 1469–1474 (2011).
Kroon, E. et al. Pancreatic endoderm derived from human embryonic stem cells generates glucose-responsive insulin-secreting cells in vivo. Nat. Biotechnol. 26, 443–452 (2008).
Neph, S. et al. Circuitry and dynamics of human transcription factor regulatory networks. Cell 150, 1274–1286 (2012).
Loh, Y.-H. et al. The Oct4 and Nanog transcription network regulates pluripotency in mouse embryonic stem cells. Nat. Genet. 38, 431–440 (2006).
Boyer, L. A. et al. Core transcriptional regulatory circuitry in human embryonic stem cells. Cell 122, 947–956 (2005).
Maniatis, T. & Tasic, B. Alternative pre-mRNA splicing and proteome expansion in metazoans. Nature 418, 236–243 (2002).
Ankö, M.-L. & Neugebauer, K. M. RNA-protein interactions in vivo: global gets specific. Trends Biochem. Sci. 37, 255–262 (2012).
Barash, Y. et al. Deciphering the splicing code. Nature 465, 53–59 (2010).
Venables, J. P., Tazi, J. & Juge, F. Regulated functional alternative splicing in Drosophila. Nucleic Acids Res. 40, 1–10 (2012).
Black, D. L. Mechanisms of alternative pre-messenger RNA splicing. Annu. Rev. Biochem. 72, 291–336 (2003).
Ule, J. et al. An RNA map predicting Nova-dependent splicing regulation. Nature 444, 580–586 (2006).
Xu, X. et al. ASF/SF2-regulated CaMKIIdelta alternative splicing temporally reprograms excitation-contraction coupling in cardiac muscle. Cell 120, 59–72 (2005).
Salomonis, N. et al. Alternative splicing regulates mouse embryonic stem cell pluripotency and differentiation. Proc. Natl Acad. Sci. USA 107, 10514–10519 (2010).
Nelles, D. A. & Yeo, G. W. Alternative splicing in stem cell self-renewal and diferentiation. Adv. Exp. Med. Biol. 695, 92–104 (2010).
Gabut, M. et al. An alternative splicing switch regulates embryonic stem cell pluripotency and reprogramming. Cell 147, 132–146 (2011).
Gopalakrishna-Pillai, S. & Iverson, L. E. A DNMT3B alternatively spliced exon and encoded peptide are novel biomarkers of human pluripotent stem cells. PLoS One 6, e20663 (2011).
Schwanhäusser, B. et al. Global quantification of mammalian gene expression control. Nature 473, 337–342 (2011).
Venables, J. P. et al. Cancer-associated regulation of alternative splicing. Nat. Struct. Mol. Biol. 16, 670–676 (2009).
Venables, J. P. et al. RBFOX2 is an important regulator of mesenchymal tissue-specific splicing in both normal and cancer tissues. Mol. Cell Biol. 33, 396–405 (2013).
Lapasset, L. et al. Rejuvenating senescent and centenarian human cells by reprogramming through the pluripotent state. Genes Dev. 25, 2248–2253 (2011).
Xing, Y. & Lee, C. Alternative splicing and RNA selection pressure—evolutionary consequences for eukaryotic genomes. Nat. Rev. Genet. 7, 499–509 (2006).
Xing, Y. & Lee, C. J. Protein modularity of alternatively spliced exons is associated with tissue-specific regulation of alternative splicing. PLoS Genet. 1, e34 (2005).
Wang, E. T. et al. Alternative isoform regulation in human tissue transcriptomes. Nature 456, 470–476 (2008).
Han, H. et al. MBNL proteins repress ES-cell-specific alternative splicing and reprogramming. Nature 498, 241–245 (2013).
Wang, E. T. et al. Transcriptome-wide regulation of pre-mRNA splicing and mRNA localization by muscleblind proteins. Cell 150, 710–724 (2012).
Rinaldi, F. et al. Aberrant splicing and expression of the non muscle myosin heavy-chain gene MYH14 in DM1 muscle tissues. Neurobiol. Dis. 45, 264–271 (2012).
Lin, X. et al. Failure of MBNL1-dependent post-natal splicing transitions in myotonic dystrophy. Hum. Mol. Genet. 15, 2087–2097 (2006).
Kalsotra, A. et al. A postnatal switch of CELF and MBNL proteins reprograms alternative splicing in the developing heart. Proc. Natl Acad. Sci. USA 105, 20333–20338 (2008).
Fugier, C. et al. Misregulated alternative splicing of BIN1 is associated with T tubule alterations and muscle weakness in myotonic dystrophy. Nat. Med. 17, 720–725 (2011).
Venables, J. P., Vignal, E., Baghdiguian, S., Fort, P. & Tazi, J. Tissue-specific alternative splicing of Tak1 is conserved in deuterostomes. Mol. Biol. Evol. 29, 261–269 (2012).
Shapiro, I. M. et al. An EMT-driven alternative splicing program occurs in human breast cancer and modulates cellular phenotype. PLoS Genet. 7, e1002218 (2011).
Zhang, Z. et al. Divergence of exonic splicing elements after gene duplication and the impact on gene structures. Genome Biol. 10, R120 (2009).
Yang, J. et al. Twist, a master regulator of morphogenesis, plays an essential role in tumor metastasis. Cell 117, 927–939 (2004).
Kimura, T. et al. Altered mRNA splicing of the skeletal muscle ryanodine receptor and sarcoplasmic/endoplasmic reticulum Ca2+-ATPase in myotonic dystrophy type 1. Hum. Mol. Genet. 14, 2189–2200 (2005).
Hino, S.-I. et al. Molecular mechanisms responsible for aberrant splicing of SERCA1 in myotonic dystrophy type 1. Hum. Mol. Genet. 16, 2834–2843 (2007).
Gates, D. P., Coonrod, L. A. & Berglund, J. A. Autoregulated splicing of muscleblind-like 1 (MBNL1) Pre-mRNA. J. Biol. Chem. 286, 34224–34233 (2011).
NCBI Resource Coordinators. Database resources of the National Center for Biotechnology Information. Nucleic Acids Res. 41, D8–D20 (2013).
Boué, S., Paramonov, I., Barrero, M. J. & Izpisúa Belmonte, J. C. Analysis of human and mouse reprogramming of somatic cells to induced pluripotent stem cells. What is in the plate? PLoS One 5, e12664 (2010).
Pan, Q. et al. Revealing global regulatory features of mammalian alternative splicing using a quantitative microarray platform. Mol. Cell 16, 929–941 (2004).
Goers, E. S., Purcell, J., Voelker, R. B., Gates, D. P. & Berglund, J. A. MBNL1 binds GC motifs embedded in pyrimidines to regulate alternative splicing. Nucleic Acids Res. 38, 2467–2484 (2010).
Charizanis, K. et al. Muscleblind-like 2-mediated alternative splicing in the developing brain and dysregulation in myotonic dystrophy. Neuron 75, 437–450 (2012).
Yeo, G. W. et al. An RNA code for the FOX2 splicing regulator revealed by mapping RNA-protein interactions in stem cells. Nat. Struct. Mol. Biol. 16, 130–137 (2009).
Zhang, C. et al. Defining the regulatory network of the tissue-specific splicing factors Fox-1 and Fox-2. Genes Dev. 22, 2550–2563 (2008).
Holt, I. et al. Defective mRNA in myotonic dystrophy accumulates at the periphery of nuclear splicing speckles. Genes Cells 12, 1035–1048 (2007).
Acknowledgements
This work was supported by the collaborative laboratory Splicos therapeutics, and Institut National du Cancer (Cancéropôle Ile-de-France). L.L. was supported by a fellowship from Fondation pour la Recherche Médicale (FRM) grant (Equipe FRM 2011 –no. DEQ20111223745, DCR20091217183), l’Association pour la Recherche contre le Cancer ARC (A09/4/5062 and SFI20111203957) and la Ligue régionale contre le Cancer (LNCC PYRENEES-ORIENTALES 2010, 2012). This work was also funded by grants from the Canadian Cancer Society and the Canadian Institute for Health Research. S. Abou Elela is the Canada Research Chair in RNA Biology and Cancer Genomics. B. Chabot is the Canada Research Chair in Functional Genomics. This research using hESC has been authorized by the French Biomedical Agency (ABM) for JM Lemaitre Laboratory, (NOR: SJSB0721833S). We thank Mathieu Durand and Elvy Lapointe for technical assistance. We thank Wolfson Centre for inherited Neuromuscular Disease for providing anti-MBNL1 antibody. We thank T. Bouschet for the gift of E14 mouse ES line. We are grateful to S. Alaoui from Covalab for design of immunogenic peptides and production of specific antibodies against PLOD2 and PALM splicing isoforms. We thank B.J. Blencowe for comments on the manuscript. We acknowledge Montpellier RIO Imaging (MRI) for the imaging acquisition and analysis.
Author information
Authors and Affiliations
Contributions
L.L. perfomed pluripotent stem cell experiments, differentiation assays and siRNA experiments on pluripotent stem cells. J.P.V. performed RT-PCR and G.G. performed siRNA experiments on fibroblasts. L.L. and G.G. performed immunofluorescence experiments, qPCR and analysis. E.V. contributed to the analysis of cell area and statistical analysis. J.P.V., P.T., R.K., P.P. and B.C. designed high-throughput RT–PCR experiments. J.P.V., R.K., P.T., P.P., B.C. and S.A.E. contributed to the analysis of the data. M.I. analysed RNA-seq and CLIP-seq data. P.F. analysed splicing evolution and gene ontology. J.P.V., L.L., G.G., P.R., J.-M.L. and J.T. were involved in study design and manuscript preparation, performed experiments and analysed data. All authors discussed the results and commented on the manuscript.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing financial interests.
Supplementary information
Supplementary Figures and Tables
Supplementary Figures S1-S7 and Supplementary Tables S1-S2 (PDF 2315 kb)
Supplementary Data 1
Splice data from 303 exons in stem cells and fibroblasts. (XLSX 115 kb)
Supplementary Data 2
Data from 47 alternative splice events in stem cells and fibroblasts. (XLSX 18 kb)
Supplementary Data 3
Original data for knockdown of the 81 potential splicing factors using a selected siRNA (XLSX 154 kb)
Supplementary Data 4
Median splicing shifts for each RBP. (XLSX 57 kb)
Rights and permissions
About this article
Cite this article
Venables, J., Lapasset, L., Gadea, G. et al. MBNL1 and RBFOX2 cooperate to establish a splicing programme involved in pluripotent stem cell differentiation. Nat Commun 4, 2480 (2013). https://doi.org/10.1038/ncomms3480
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/ncomms3480
This article is cited by
-
Effect of RNA-binding proteins on osteogenic differentiation of bone marrow mesenchymal stem cells
Molecular and Cellular Biochemistry (2024)
-
RBFOX2 deregulation promotes pancreatic cancer progression and metastasis through alternative splicing
Nature Communications (2023)
-
Myotonic dystrophy RNA toxicity alters morphology, adhesion and migration of mouse and human astrocytes
Nature Communications (2022)
-
Dynamics of alternative splicing during somatic cell reprogramming reveals functions for RNA-binding proteins CPSF3, hnRNP UL1, and TIA1
Genome Biology (2021)
-
Splicing reprogramming of TRAIL/DISC-components sensitizes lung cancer cells to TRAIL-mediated apoptosis
Cell Death & Disease (2021)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.