Figures
Abstract
Alternative splicing patterns are known to vary between tissues but these patterns have been found to be predominantly peculiar to one species or another, implying only a limited function in fundamental neural biology. Here we used high-throughput RT-PCR to monitor the expression pattern of all the annotated simple alternative splicing events (ASEs) in the Reference Sequence Database, in different mouse tissues and identified 93 brain-specific events that shift from one isoform to another (switch-like) between brain and other tissues. Consistent with an important function, regulation of a core set of 9 conserved switch-like ASEs is highly conserved, as they have the same pattern of tissue-specific splicing in all vertebrates tested: human, mouse and zebrafish. Several of these ASEs are embedded within genes that encode proteins associated with the neuronal microtubule network, and show a dramatic and concerted shift within a short time window of human neural stem cell differentiation. Similarly these exons are dynamically regulated in zebrafish development. These data demonstrate that although alternative splicing patterns often vary between species, there is nonetheless a core set of vertebrate brain-specific ASEs that are conserved between species and associated with neural differentiation.
Citation: Madgwick A, Fort P, Hanson PS, Thibault P, Gaudreau M-C, Lutfalla G, et al. (2015) Neural Differentiation Modulates the Vertebrate Brain Specific Splicing Program. PLoS ONE 10(5): e0125998. https://doi.org/10.1371/journal.pone.0125998
Academic Editor: Emanuele Buratti, International Centre for Genetic Engineering and Biotechnology, ITALY
Received: February 13, 2015; Accepted: March 27, 2015; Published: May 19, 2015
Copyright: © 2015 Madgwick et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited
Data Availability: Data are available within the manuscript and supporting information and on our web site (http://rnomics.med.usherbrooke.ca:3000/pcrreactiongroup/list/336).
Funding: The research leading to these results was primarily funded by the People Programme (Marie Curie Actions) of the European Union's Seventh Framework Programme (FP7/2007-2013) under REA grant agreement awarded to JPV n° [PCIG13-GA-2013- 618114]. DJE was funded by Wellcome trust grants WT080368MA and WT089225/Z/09/Z. TM was supported by CIHR grant (MOP- 86516). M-C G was supported by a fellowship from the Canadian Institutes of Health Research (CIHR). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
How both complex and simple body plans (e.g. human and nematode) can be encoded by of the order of only 20,000 genes is a major conundrum in genome biology. One part of the solution to this is alternative splicing, a mechanism through which different RNAs and protein products can be made from a single gene by differential incorporation of regions of pre-mRNA. Alternative splicing is a key player in gene expression of complex organisms like mammals and in complex tissues like the brain as it increases transcriptome and thereby proteome diversity. For two decades, alternative splicing has been studied on a gene-by-gene basis. In the last fifteen years a significant amount of transcriptome-wide data from microarray and RNA-Seq experiments has increased the number of known alternative splice events (ASEs). Recent developments have brought to light that at least 80% of genes produce alternative mRNAs. Which of these ASEs are important for development is now a major question in the field. One route to answer this question is through evolutionary comparisons, on the assumption that natural selection has maintained the most important regulated ASEs.
High-throughput RT-PCR is a valuable method to identify splicing changes in many genes at once [1–3] and is an alternative to current methods including RNA-Seq [2,4]. We have recently used a cross-evolutionary approach coupled to high-throughput RT-PCR to identify the most conserved tissue-specific alternative exons, regulated across a variety of different organisms. Starting with a list of 40 human alternative exons we found one that was alternatively spliced in muscle across all deuterostomes, that is, including vertebrates, chordates, echinoderms and hemichordates [5]. The identification of highly conserved patterns of alternative splicing in a specific tissue within animals is important, since it reveals exceptions to other recent transcriptome-wide studies that show, using RNA-Seq, that alternative splicing is for the most part specific to individual species and not to tissues [6,7]. These results suggested that alternative splicing has only a minor role if any in formulating the vertebrate body plan.
Extreme changes, from near total exon inclusion to exon skipping, known as switch-like splicing, have been observed in human tissues, but not in all vertebrate orthologues [8]. Here we use high-throughput RT-PCR in a cross-evolutionary approach to identify a new panel of ASEs that are highly tissue-specific in vertebrates, and have identical splicing patterns in the brain of mammals and zebrafish. Using a stem cell model we find that these exons are a conserved component of alternative splicing changes that shift abruptly during neural stem cell differentiation, and might therefore be important for conserved features of vertebrate central nervous system development.
Results
RT-PCR mining of tissue-specific splicing reveals 93 switch-like brain-specific exons in mouse
In order to first profile the extent of tissue-specific alternative splicing pattern variations in humans, we studied a panel of 47 [1] alternative splicing events (ASEs) across 6 human tissue cDNAs (Clontech, Mountain View, CA) by RT-PCR. These tissues were liver, kidney, lung, muscle, heart and brain (S1 Fig). As expected from previous high-throughput studies [7,9] these splicing events could be used to define three groups: brain with the most distinctive profile, muscle and heart with another profile and liver, kidney and lung as the third group. The largest difference in splicing profiles was between brain and liver/kidney. Therefore to select for important ASEs across the entire transcriptome we then compared splicing profiles of mouse brain, liver and kidney with a panel of 1329 (functionally random) ASEs composed of all the simple ASEs in the mouse RefSeq database [10]. After quality control, 809 PCRs produced robust data for all three samples. 70% of these (560/809) were simple cassette exons. The other events fell into three types: alternative 5’ sites (100), alternative 3’ sites (102) and more complex splice events (47). 93 ASEs showed a 50% shift in exon inclusion between brain and both other tissues (Fig 1; all splicing data is available numerically in S1 Table and all original electopherograms are at http://rnomics.med.usherbrooke.ca:3000/pcrreactiongroup/list/336). 84% of these brain-specific ASEs (78/93) were simple cassette exons, therefore conserved alternative splices are highly enriched in alternative exons compared to other types of alternative splicing (p = 0.00005 Fisher’s exact test).
PCR primers were designed across all the simple ASEs in the mouse RefSeq collection [10] and RT-PCR was performed on three mouse tissue cDNA;. 809 PCRs gave good data; the psi values are shown for mouse brain, kidney and liver. The ASEs have been clustered according to the shift between brain and the two other tissues, indicated by the difference in psi values. 93 ASE psi values shifted more than 50% between brain and both of the other tissues (see the top and bottom of the heat map). Numerical data are given in S1 Table.
Nine alternative splice events are regulated in a switch-like brain-specific manner from man to zebrafish
To test the conservation of the splicing patterns of these alternative events between mouse and humans, primers were designed to the 89 orthologous human ASEs, we could find, of the 93 brain-regulated mouse ASEs. Splicing profiles were then interrogated by RT-PCRs in six human tissues, brain, liver, kidney, lung, muscle and heart. From these, 62 ASEs produced robust data in all tissues (Fig 2). Combining the percent-spliced-in data for these 62 ASEs from mouse and human (from Figs 1 and 2 respectively) clustered human and mouse tissues together, demonstrating that splicing of this group of ASEs is tissue-specific rather than species-specific (S2 Fig). Furthermore, of all the tissue-specific ASEs that showed similar splicing patterns in mouse and humans, 15 ASEs had a shift of more than 50% in splicing between brain and all of the other tissues. The conservation of the switch-like brain-specific splicing of these 15 ASEs was then further investigated in zebrafish brain, kidney and liver, and 9 of the ASEs were confirmed to be brain-specific in zebrafish: EML4, KIF2A, MAP2, APLP2, NEO1, CLIP2, SMG7, PHF21A and KIAA0513 (Fig 3). We mapped the location of these 9 ASEs onto their respective protein structures, and found exons encoding both known structural domains and inter-domain regions (illustrated diagrammatically in Fig 4a).
Of the 93 mouse candidate ASEs for conserved brain-specificity, 89 had orthologous regions in human. Of these, 62 ASEs gave good data by PCR on 6 human tissue cDNAs: brain, heart, muscle, kidney, liver and lung. The psi values are shown in the heat map; the genes are clustered on both axes according to their psi values. X-axis clustering shows the brain, heart/muscle and liver/kidney/lung have 3 distinct splicing profiles. Y axis clustering groups the genes based on their psi values’ patterns across tissues; the brain-specific ASEs cluster at the bottom of the heat map but the clustering also distinguishes some ASEs in the centre that splice similarly in brain/heart/muscle as distinct from liver/kidney/lung.
RT-PCR was performed on 15 genes across human, mouse and zebrafish. The 9 genes shown have conserved switch-like splicing in all three vertebrate species. Brain-specific splice forms are indicated with a red arrow. The alternative kidney/liver-specific forms are indicated with a double-headed blue arrow. The different, expected and found, PCR sizes for the long and short form of each gene in each species are given in S1 Table.
a. Primary structures of proteins encoded by the 9 human genes with vertebrate conserved brain-specific splicing. Shown are the annotated human proteins with the regions of the nine splice events indicated by red boxes. TM: transmembrane region; β-APP C-ter: C-terminus of the β Amyloid Precursor Protein (pf10515); CAP-Gly: Cytoskeletal-Associated protein (pf00225); HELP: Hydrophibic EMAP-Like Protein (pf03451); WD40: β-transducin repeat (pf00400); SiP: Signal peptide; Ig: Immunoglobulin-like domain (pf00047); FN3: Fibronectin type III domain (pf00041); AT-hook: DNA-binding for A/T-rich regions (pf02178); PHD finger: Plant HomeoDomain (Cys)4-His-(Cys)3 (pf00628). Gene names are labelled with a 1 if exclusion of the alternatively splice regions directly affects structural domains. Note all ASEs are multiples of 3 nucleotides, thus all the alternative splicing events confer in frame peptide omission or insertion. b. Brain-specific alternative splicing is conserved in vertebrates, and possibly beyond, in microtubule-associated genes. Metazoan genomes in Ensembl were searched for paralogs and orthologs of each target gene and for the presence (yellow) or absence (gray) of potential ASEs (Accession numbers are shown in S1 Table). Yellow indicates that alternatively spliced mRNAs were detected in EST databases. Orange indicates that there were too few ESTs to conclude. Green indicates the absence of both genomic and EST data. Duplications are indicated by thick lines along with the names of the duplicated genes. An indication of the function of each gene is given on the right.
Next, we investigated the extent to which the vertebrate brain-specific splicing process might be conserved, as alternative spliced regions, further in evolution than bony vertebrates. The analysis of EST databases by peptide homology searching revealed alternative splicing in genes for NEO1, MAP2 and KIF2A in jawless vertebrates and the latter was also conserved in tunicates, the sister taxon of vertebrates (Fig 4b).
The conserved splices are dynamically regulated in diverse physiological contexts from man to fish
To understand how alternative splicing of these conserved ASEs might be modulated within an intact organism, we monitored them during the first 48 hours of embryogenic zebrafish development (Fig 5a). Each showed a dynamic pattern over this time frame. If these brain-specific splicing events are important for brain, they would also be expected to change during the course of neural stem cell differentiation. To test this hypothesis, we turned to an in vitro differentiation system (N1997 cells) that converts stem cells to a mixture of neurones, astrocytes and oligodendrocytes [11,12]. We tested all 89 mammalian brain-specific ASEs and observed a sharp change at 6–10 days of differentiation with 13 ASEs shifting over 50% (i.e. in a switch-like manner) in one direction, and all at or around this time point (Fig 5b).
a. The 9 conserved brain-specific ASEs shifts during zebrafish embryogenesis. Time-course of splicing of the nine vertebrate brain-specific ASEs during the first two days of zebrafish embryogenesis, spanning 2 to 48 hours post fertilisation (hpf). The red arrows indicate the commencement of expression of the brain-specific AS transcript. b. Splicing of 13 genes pivots from predominantly one AS transcript isoform to the other upon neural differentiation of stem cells between days 6–10. 93 ASEs were assayed by PCR and their psi values were evaluated during neural stem cell differentiation. 56 ASEs gave good data for all time points (see S1 Table). Chart showing the psi values of 13 ASEs whose psi values shift more than 50% during the time-course.
Conserved alternative splicing targets the cytoskeleton
The gene encoding Echinoderm Microtubule-Associated Protein-Like 4 (EML4), which affects the assembly dynamics of microtubules [13], was identified as undergoing among the largest shifts in both human stem cell differentiation and during zebrafish embryogenesis. Therefore we investigated the effect of splicing of EML4 at the protein level by Western blotting and confirmed a shift from the exon-included to the exon-skipped form of EML4 protein during stem cell differentiation (S3 Fig).
As over half of the evolutionarily conserved splicing changes identified in this study are involved in microtubule and or cytoskeleton dynamics we tested a panel of microtubule-related genes for their tissue-specific splicing patterns. The strongest shifts were displayed by EML1 (a paralog of EML4) and KIF13A, which is a gene encoding a kinesin microtubule-based motor, like KIF2A (S4 Fig). Another kinesin, KIF1b had already been shown to be alternatively spliced in brain [14,15] and it is regulated along with KIF2A and KIF5A by PTBP1 which is down-regulated in brain [16–18]. We thus analysed the splicing of 65 exons in 29 kinesin genes (S1 Table and http://rnomics.med.usherbrooke.ca:3000/data/related/2149). Just 17 ASEs gave good data in all six tissue tested. This analysis confirmed brain-specific splicing of KIF2A and KIF13B.
Discussion
In this study, we identify a panel of alternative splice events (ASEs) that switch dramatically in their inclusion levels between brain and other tissues of three distantly related vertebrate species: human, mouse and zebrafish. The brain has long been known to have a highly complex transcriptome due to alternative splicing [19] but the conservation of regulation of these ASEs has not been thoroughly investigated. Early indications, from microarray technology, indicated that even in our closest cousins the chimpanzees, there was relatively little conservation of actual observed alternative splicing patterns, even though the genomic regions encoding ASEs were near identical [20]. Many studies, therefore, have documented brain-specific ASEs in one species or another but ours is the first study to verify a group of such switch-like [8] ASEs across vertebrate taxa. Previous studies comparing ASEs across taxa found relatively few brain-specific ASEs [6,21,22], however two ASEs, in the NEO1 and APLP2 genes, have previously been confirmed to be have conserved switch-like brain-specific alternative splicing in vertebrates [21]. Our study confirmed these two switch-like vertebrate brain-specific splicing events but also adds seven others (Fig 3). These brain-specific splices are likely of importance for neural biology as four of them: NEO1, KIF2A, APLP2, and MAP2 were shown to be alternatively spliced during neuronal differentiation of P19 mouse embryonal carcinoma cells [23] and the latter three ASEs were also shown to be targets of PTBP1, which is a key splicing factor that is down-regulated during neurogenesis [16,17].
Several of the affected genes have functions in the cytoskeleton and they are particularly associated with microtubules. Microtubule-associated proteins act as important regulators in the staged development of neurite, axon and dendrite formation, and the appearance of alternative splicing of these genes following neural stem cell differentiation indicates that these are key events in these processes [24]. These results also support the hypothesis that the main function of alternative splicing is to enact a program of changes in the cytoskeleton between different cell types [25–27]. The microtubule-associated genes with conserved brain-specific splicing are EML4 (echinoderm microtubule associated protein like 4), MAP2 (microtubule-associated protein 2), KIF2A (microtubule-depolymerizing kinesin) [28] and CLIP2 which associates with the ends of growing microtubules [29]. Furthermore, APLP2 is involved in actin dynamics [30] and KIAA0513 interacts with the cytoskeleton [31]. Conservation analysis indicated that alternative exons of two of the microtubule-associated genes (MAP2, KIF2A) are even more ancient than vertebrate radiation (Fig 4b) and therefore regulation of their splicing may be an important ancient feature of sophisticated neural circuitry in deuterostomes. These alternatively spliced genes also have important functions in the brain. Although mice lacking MAP2 were normal, combined knockdown of MAP2 and MAP1B led to dendritic outgrowth problems [32]. KIF2A knockdown mice show aberrant axonal branching [33]. CLIP2 knockout mice have hippocampal dysfunction and motor coordination defects [34].
Five of the nine conserved ASEs (EML4, CLIP2 SMG7, PHF21A, KIAA0513) have not previously been reported to be alternatively spliced in any species’ brain. The ASE detected in the gene for the microtubule-associated protein EML4 was the most strongly altered event in all our validation assays. EML4 splicing shifts nearly completely during zebrafish embryogenesis (Fig 5a) as well as in stem cell differentiation (Fig 5b) where we confirmed the change at the protein level (S3 Fig). EML4 is currently under intense investigation in non small cell lung cancer where it is frequently found fused to the ALK gene [35]. In a further screen of microtubule-associated splice events, we also found brain-specific splicing of its homologue EML1 and of, another kinesin microtubule motor protein gene, KIF13A (S4 Fig). We previously found KIF13A exon inclusion during stem cell differentiation to primary mesoderm under the control of the MBNL1 and RBFOX2 splicing factors [2]. As all our observed ultra-conserved splicing changes are in-frame events (Fig 4a), this is consistent with the idea that the function of alternative splicing is to regulate protein interaction networks by the removal and addition of protein interaction domains [36] which might be especially important during the complex cytoskeletal changes occurring during neurogenesis [24]. To conclude, our study reveals the existence of at least 9 vertebrate-conserved switch-like brain-specific ASEs, which very likely have a direct bearing on the establishment of the neuronal microtubule network.
Methods
Ethics statement
Zebrafish embryogenesis was specifically approved by the Direction Sanitaire et Vétérinaire de l’Hérault and Comité d’Ethique pour l’Expérimentation Animale under reference CEEA-LR-13007. Zebrafish were handled according to European Union guidelines for the handling of laboratory animals (http://ec.europa.eu/environment/chemicals/lab_animals/home_en.htm). Adult mice were sacrificed by carbon dioxide inhalation as per standard operating procedures approved by the IRCM ACC and the CCAC. All efforts were made to minimize the number of animals used and to reduce their suffering. The protocols for the mice used in this study were reviewed and specifically approved by the IRCM Animal Care Committee (ACC, protocol number: 2009–12). Mice were sacrificed and organs were prepared according to institutional rules put in place by the IRCM Animal Care Committee (ACC Quebec Comité de protection des animaux), which follows the regulations and requirements of the Canadian Council on Animal Care (CCAC, www.ccac.ca)
Animals
Fish maintenance, staging and husbandry were as described previously [37] and performed under standard conditions [38]. The F1 golden zebrafish mutant [39] originally purchased from Antinea (http://groupe-antinea.fr/) was used. Embryos were obtained from pairs of adult fish by natural spawning and raised at 28.5°C. Embryos and larvae were staged as described previously [40]. Mice were housed under specific pathogen-free conditions and certified animal technicians regularly observed the mice. All mice were from the C57BL/6 genetic background and were offspring from animals originally obtained from The Jackson Laboratory (Bar Harbor, Maine, USA).
High-throughput splicing experiments.
were performed as in [3]. For analysis of splicing data, a quality control was used whereby PCRs were only considered with purity >75%, that is with >75% of the total molarity of the found peaks, found at the expected molecular weights. Percent spliced in (psi) values were calculated as the molarity of the long form divided by the molarity of the long and short combined. Manual PCRs and electrophoresis were performed as in [2]. All data is visible on our website http://rnomics.med.usherbrooke.ca:3000/pcrreactiongroup/list/336. Each section contains clickable graphics which leads to the primary electrophoresis data. Links are also available to download all the calculated splicing data in Excel format.
Human neuronal cells.
were generated from the human neuronal precursor stem cell (hNPSC) line N1997 [11,12]. The hNPSCs were seeded onto geltrex-coated cell culture plates in proliferation media (Advanced Dulbecco’s Modified Eagle’s Medium F12—Life Technologies) supplemented with 1% N-2 solution (Life Technologies), 2% B27 solution (Life Technologies), 1% ITS (10 μg ml-1 insulin, 5.5 μg ml-1 transferrin, 5 ng ml-1 sodium selenite solution), 1% PSG (100 units ml-1 penicillin, 100 μg ml-1 streptomycin, 2 mM glutamine), 2.5 μg ml-1 fungizone (Amphotericin B, Life Technologies), 1% MEM non-essential amino acids (100 X, Life Technologies), 2% sodium bicarbonate (7.5% solution), 1.5% D-glucose (45% in H2O), 20 mM L-ascorbic acid, 5 μg ml-1 heparin, 20 ng ml-1 epidermal growth factor (EGF, Sigma Aldrich), 20 ng ml-1 basic fibroblast growth factor (FGF2, R and D Systems) and 10 ng ml-1 leukaemia inhibitory factor (LIF, Sigma Aldrich). The hNPSC were cultured in proliferation media for 4 days prior to the start of differentiation into mature neural cells. To induce terminal differentiation, the proliferation media was withdrawn and replaced with differentiation media, which was prepared as the proliferation media but without the addition of EGF, FGF2 and LIF and with the addition of 10% fetal bovine serum (Sigma Aldrich). The hNPSCs were differentiated for 14 days with differentiation media being refreshed every other day. At days 0, 2, 6, 10 and 14 of differentiation, the cells were harvested for PCR by the addition of 500 μl of Trizol reagent (Sigma Aldrich) per well, after which the cells were scraped from the surface and collected. The protein samples for western blotting were collected by washing the cells with Dulbecco’s phosphate buffered saline prior to lysis with native lysis buffer.
Western blotting.
was performed by the addition of 10 μg of protein from the stem cell lysates, determined by a modified Bradford assay, to the lanes of 4–12% Bis-Tris Midi gels (Life Technologies). The proteins were separated by electrophoresis, then transferred to a nitrocellulose membrane using the iBlot gel transfer system (Life Technologies). The nitrocellulose membrane was blocked with tris-buffered saline with 0.05% TWEEN 20 and 5% w/v non-fat dried skimmed milk, prior to probing with primary antibodies EML4 (HPA036688 Sigma Aldrich) and GAPDH, labelled with HRP (sc-25778 HRP Santa-Cruz), all at 1:4000 dilutions. The anti-EML4 antibody was detected with sheep anti-rabbit antibody, labelled with HRP (ab6795 Abcam). The substrate ECL 2 (Fisher Scientific) was added, and the membranes were exposed to FujiFilm Super RX films (Fisher Scientific) prior to development in Kodak GBX developer and fixer solutions.
Conservation analysis.
was performed as in [2]. Orthologous genes and exons corresponding to alternatively spliced variants were searched at http://ensembl.org (vertebrates, tunicates, D. melanogaster and C. elegans) and at http://www.metazome.net/ (cnidarians). Splicing variants were identified in EST databases (http://blast.ncbi.nlm.nih.gov/). Accession numbers are in S1 Table.
Supporting Information
S1 Fig. High-throughput RT-PCR shows brain has the most distinct splicing profile amongst six human tissues
. PCR was performed across regions of alternative splicing in 47 diverse genes on 6 human tissue cDNA libraries. The percent-spliced-in (psi) values were calculated and plotted as a heat map. The psi values are indicated by a colour going from yellow to blue depending on the extent to which the alternative exon is spliced in for that tissue. The six tissues were clustered, based on their psi values, which shows brain has the most distinct splicing profile in human.
https://doi.org/10.1371/journal.pone.0125998.s001
(TIF)
S2 Fig. A group of mammalian tissue-specific ASEs.
Data from human ASEs in Fig 2 was combined with data for the orthologous mouse ASEs from Fig 1. The psi values are shown in the heat map; the genes were ordered according to the shift between human brain and the closest of the other five human tissues. Note that brain clusters separately from the other tissues, irrespective of species. Note also, the 4 top-most and 11 bottom-most ASEs psi values shift more than 50% between human brain and the other five tissues.
https://doi.org/10.1371/journal.pone.0125998.s002
(TIF)
S3 Fig. Splicing shifts to the short form of the EML4 protein during stem cell differentiation.
A. Representative manual PCRs for the most significant shifts in either direction. Times are shown in days. B. Western blot showing EML4 expression during stem cell differentiation. Biological triplicate samples were probed with anti-EML4 and anti-GAPDH antibody as a loading control. Note the shift in splicing predicted by the PCR experiments is verified at the protein level. Note also, a shorter unidentified protein is visible in all samples except the fully matured 14 day samples.
https://doi.org/10.1371/journal.pone.0125998.s003
(TIF)
S4 Fig. Brain-specific splicing of microtubule-associated genes.
94 exons in 74 microtubule-associated genes were assayed in 6 human tissues. Heat-map showing percent-spliced-in values for 39 alternative exons that gave good data for all six tissues. The genes were clustered, depending on their psi values patterns across the 6 tissues. Note, brain showed the most significant tissue-specific splicing.
https://doi.org/10.1371/journal.pone.0125998.s004
(TIF)
S1 Table.
Tab1. Splicing data for cross-evolutionary comparison. All data superior to the 75% purity cut off (see methods) in all tissues is shown. The nine ultra-conserved regulated human ASEs are shown in red. Column a. Human ASE names (gene names followed by a letter if more than one ASE occurs in this gene). b. Mouse gene name. c. Mouse ASE name (gene names followed by the ASE size if more than one ASE occurs in this gene). d. ASE type. e. Mouse in-house primer names separated by a double space. f,g. Mouse forward and reverse primer sequences. h. Mouse upstream intron sequence (150ntds). i. ASE sequence. j. Mouse downstream intron sequence 150 ntds. k-m. Mouse small and large PCR product and ASE sizes (l-k) respectively. n-p. Mouse brain, kidney and liver psi values respectively. q. Shift between brain psi and the nearest of liver or kidney. r. Human gene name again. s. Human in-house primer names separated by a double space. t,u. Human forward and reverse primer sequences. v-x. human small and large PCR product and ASE sizes (w-v) respectively. y-ad. Human brain, kidney, liver, lung, muscle and heart psi values respectively. ae-ai. Human neural stem cell psi values at 0, 2, 6, 10 and 14 days of differentiation. aj,ak. Zebrafish forward and reverse primer sequences. al,am. Zebrafish small and large PCR product sizes. Tab2. Splicing data for microtubule-related genes (rows 2–44) and Kif genes (rows 45–58). a. Gene name b,c. Short and long PCR product sizes. d. Alternative splice event type. e. In-house primer names. f,g. Forward and reverse primer sequences. h-m. Percent-spliced-in values in the indicated tissues. Tab3. Accession numbers for evolutionary analysis of alternative splice events.
https://doi.org/10.1371/journal.pone.0125998.s005
(XLS)
Acknowledgments
We are indebted to Marie-Claude Lavallée, Mélanie St-Germain and Jade Dussureault for animal care. Thanks also to Elvy Lapointe, Mathieu Durand and Roscoe Klinck of the Laboratoire de Génomique Fonctionnelle de l'Université de Sherbrooke high-throughput RT-PCR platform.
Author Contributions
Conceived and designed the experiments: GL TM CMM JPV. Performed the experiments: AM PF PSH PT M-CG GL BC. Analyzed the data: AM PF PT JPV. Wrote the paper: AM PF PSH GL TM SAE BC DJE CMM JPV.
References
- 1. Venables JP, Brosseau JP, Gadea G, Klinck R, Prinos P, Beaulieu JF, et al. (2013) RBFOX2 is an important regulator of mesenchymal tissue-specific splicing in both normal and cancer tissues. Mol Cell Biol 33: 396–405. pmid:23149937
- 2. Venables JP, Lapasset L, Gadea G, Fort P, Klinck R, Irimia M, et al. (2013) MBNL1 and RBFOX2 cooperate to establish a splicing programme involved in pluripotent stem cell differentiation. Nat Commun 4: 2480. pmid:24048253
- 3. Klinck R, Bramard A, Inkel L, Dufresne-Martin G, Gervais-Bird J, Madden R, et al. (2008) Multiple alternative splicing markers for ovarian cancer. Cancer Res 68: 657–663. pmid:18245464
- 4. Han H, Irimia M, Ross PJ, Sung HK, Alipanahi B, David L, et al. (2013) MBNL proteins repress ES-cell-specific alternative splicing and reprogramming. Nature 498: 241–245. pmid:23739326
- 5. Venables JP, Vignal E, Baghdiguian S, Fort P, Tazi J (2012) Tissue-specific alternative splicing of Tak1 is conserved in deuterostomes. Mol Biol Evol 29: 261–269. pmid:21873631
- 6. Merkin J, Russell C, Chen P, Burge CB (2012) Evolutionary dynamics of gene and isoform regulation in Mammalian tissues. Science 338: 1593–1599. pmid:23258891
- 7. Barbosa-Morais NL, Irimia M, Pan Q, Xiong HY, Gueroussov S, Lee LJ, et al. (2012) The evolutionary landscape of alternative splicing in vertebrate species. Science 338: 1587–1593. pmid:23258890
- 8. Wang ET, Sandberg R, Luo S, Khrebtukova I, Zhang L, Mayr C, et al. (2008) Alternative isoform regulation in human tissue transcriptomes. Nature 456: 470–476. pmid:18978772
- 9. Sugnet CW, Srinivasan K, Clark TA, O'Brien G, Cline MS, Wang H, et al. (2006) Unusual intron conservation near tissue-regulated exons found by splicing microarrays. Plos Computational Biology 2: 22–35.
- 10. Ehrmann I, Dalgliesh C, Liu Y, Danilenko M, Crosier M, Overman L, et al. (2013) The tissue-specific RNA binding protein T-STAR controls regional splicing patterns of neurexin pre-mRNAs in the brain. PLoS Genet 9: e1003474. pmid:23637638
- 11. Kurzawa-Akanbi M, Hanson PS, Blain PG, Lett DJ, McKeith IG, Chinnery PF, et al. (2012) Glucocerebrosidase mutations alter the endoplasmic reticulum and lysosomes in Lewy body disease. J Neurochem 123: 298–309. pmid:22803570
- 12. Nisar R, Hanson PS, He L, Taylor RW, Blain PG, Morris CM (2015) Diquat Causes Caspase-Independent Cell Death in SH-SY5Y Cells by Production of ROS Independently of Mitochondria. Arch Toxicol in press.
- 13. Pollmann M, Parwaresch R, Adam-Klages S, Kruse ML, Buck F, Heidebrecht HJ (2006) Human EML4, a novel member of the EMAP family, is essential for microtubule formation. Exp Cell Res 312: 3241–3251. pmid:16890222
- 14. Matsushita M, Yamamoto R, Mitsui K, Kanazawa H (2009) Altered motor activity of alternative splice variants of the mammalian kinesin-3 protein KIF1B. Traffic 10: 1647–1654. pmid:19744141
- 15. Dillman AA, Hauser DN, Gibbs JR, Nalls MA, McCoy MK, Rudenko IN, et al. (2013) mRNA expression, splicing and editing in the embryonic and adult mouse cerebral cortex. Nat Neurosci 16: 499–506. pmid:23416452
- 16. Xue Y, Zhou Y, Wu T, Zhu T, Ji X, Kwon YS, et al. (2009) Genome-wide analysis of PTB-RNA interactions reveals a strategy used by the general splicing repressor to modulate exon inclusion or skipping. Mol Cell 36: 996–1006. pmid:20064465
- 17. Boutz PL, Stoilov P, Li Q, Lin CH, Chawla G, Ostrow K, et al. (2007) A post-transcriptional regulatory switch in polypyrimidine tract-binding proteins reprograms alternative splicing in developing neurons. Genes Dev 21: 1636–1652. pmid:17606642
- 18. Yap K, Lim ZQ, Khandelia P, Friedman B, Makeyev EV (2012) Coordinated regulation of neuronal mRNA steady-state levels through developmentally controlled intron retention. Genes Dev 26: 1209–1223. pmid:22661231
- 19. Stamm S, Zhang MQ, Marr TG, Helfman DM (1994) A sequence compilation and comparison of exons that are alternatively spliced in neurons. Nucleic Acids Res 22: 1515–1526. pmid:8202349
- 20. Calarco JA, Xing Y, Caceres M, Calarco JP, Xiao X, Pan Q, et al. (2007) Global analysis of alternative splicing differences between humans and chimpanzees. Genes Dev 21: 2963–2975. pmid:17978102
- 21. Jelen N, Ule J, Zivin M, Darnell RB (2007) Evolution of Nova-dependent splicing regulation in the brain. PLoS Genet 3: 1838–1847. pmid:17937501
- 22. Barbosa-Morais NL, Carmo-Fonseca M, Aparicio S (2006) Systematic genome-wide annotation of spliceosomal proteins reveals differential gene family expansion. Genome Res 16: 66–77. pmid:16344558
- 23. Suzuki H, Osaki K, Sano K, Alam AH, Nakamura Y, Ishigaki Y, et al. (2011) Comprehensive analysis of alternative splicing and functionality in neuronal differentiation of P19 cells. PLoS One 6: e16880. pmid:21365003
- 24. Poulain FE, Sobel A (2010) The microtubule network and neuronal morphogenesis: Dynamic and coordinated orchestration through multiple players. Mol Cell Neurosci 43: 15–32. pmid:19660553
- 25. Das D, Clark TA, Schweitzer A, Yamamoto M, Marr H, Arribere J, et al. (2007) A correlation with exon expression approach to identify cis-regulatory elements for tissue-specific alternative splicing. Nucleic Acids Res 35: 4845–4857. pmid:17626050
- 26. Venables JP, Klinck R, Koh C, Gervais-Bird J, Bramard A, Inkel L, et al. (2009) Cancer-associated regulation of alternative splicing. Nat Struct Mol Biol 16: 670–676. pmid:19448617
- 27. Zhang C, Zhang Z, Castle J, Sun S, Johnson J, Krainer AR, et al. (2008) Defining the regulatory network of the tissue-specific splicing factors Fox-1 and Fox-2. Genes Dev 22: 2550–2563. pmid:18794351
- 28. Uehara R, Tsukada Y, Kamasaki T, Poser I, Yoda K, Gerlich DW, et al. (2013) Aurora B and Kif2A control microtubule length for assembly of a functional central spindle during anaphase. J Cell Biol 202: 623–636. pmid:23960144
- 29. Akhmanova A, Hoogenraad CC, Drabek K, Stepanova T, Dortland B, Verkerk T, et al. (2001) Clasps are CLIP-115 and -170 associating proteins involved in the regional regulation of microtubule dynamics in motile fibroblasts. Cell 104: 923–935. pmid:11290329
- 30. Pandey P, Rachagani S, Das S, Seshacharyulu P, Sheinin Y, Naslavsky N, et al. (2014) Amyloid precursor-like protein 2 (APLP2) affects the actin cytoskeleton and increases pancreatic cancer growth and metastasis. Oncotarget.
- 31. Lauriat TL, Dracheva S, Kremerskothen J, Duning K, Haroutunian V, Buxbaum JD, et al. (2006) Characterization of KIAA0513, a novel signaling molecule that interacts with modulators of neuroplasticity, apoptosis, and the cytoskeleton. Brain Res 1121: 1–11. pmid:17010949
- 32. Teng J, Takei Y, Harada A, Nakata T, Chen J, Hirokawa N (2001) Synergistic effects of MAP2 and MAP1B knockout in neuronal migration, dendritic outgrowth, and microtubule organization. J Cell Biol 155: 65–76. pmid:11581286
- 33. Homma N, Takei Y, Tanaka Y, Nakata T, Terada S, Kikkawa M, et al. (2003) Kinesin superfamily protein 2A (KIF2A) functions in suppression of collateral branch extension. Cell 114: 229–239. pmid:12887924
- 34. Hoogenraad CC, Koekkoek B, Akhmanova A, Krugers H, Dortland B, Miedema M, et al. (2002) Targeted mutation of Cyln2 in the Williams syndrome critical region links CLIP-115 haploinsufficiency to neurodevelopmental abnormalities in mice. Nat Genet 32: 116–127. pmid:12195424
- 35. Pillai RN, Ramalingam SS (2012) The biology and clinical features of non-small cell lung cancers with EML4-ALK translocation. Curr Oncol Rep 14: 105–110. pmid:22311682
- 36. Ellis JD, Barrios-Rodiles M, Colak R, Irimia M, Kim T, Calarco JA, et al. (2012) Tissue-specific alternative splicing remodels protein-protein interaction networks. Mol Cell 46: 884–892. pmid:22749401
- 37. Alibaud L, Rombouts Y, Trivelli X, Burguiere A, Cirillo SL, Cirillo JD, et al. (2011) A Mycobacterium marinum TesA mutant defective for major cell wall-associated lipids is highly attenuated in Dictyostelium discoideum and zebrafish embryos. Mol Microbiol 80: 919–934. pmid:21375593
- 38.
Westerfield M (2007) The Zebrafish Book: A Guide for the Laboratory Use of Zebrafish (Danio rerio): Eugene: University of Oregon Press.
- 39. Lamason RL, Mohideen MA, Mest JR, Wong AC, Norton HL, Aros MC, et al. (2005) SLC24A5, a putative cation exchanger, affects pigmentation in zebrafish and humans. Science 310: 1782–1786. pmid:16357253
- 40. Kimmel CB, Ballard WW, Kimmel SR, Ullmann B, Schilling TF (1995) Stages of embryonic development of the zebrafish. Dev Dyn 203: 253–310. pmid:8589427