- Research article
- Open Access
Computational extraction of a neural molecular network through alternative splicing
BMC Research Notes volume 7, Article number: 934 (2014)
Generally, the results of high throughput analyses contain information about gene expressions, and about exon expressions. Approximately 90% of primary protein-coding transcripts undergo alternative splicing in mammals. However, changes induced by alternative exons have not been properly analyzed for their impact on important molecular networks or their biological events. Even when alternative exons are identified, they are usually subjected to bioinformatics analysis in the same way as the gene ignoring the possibility of functionality change because of the alteration of domain caused by alternative exon. Here, we reveal an effective computational approach to explore an important molecular network based on potential changes of functionality induced by alternative exons obtained from our comprehensive analysis of neuronal cell differentiation.
From our previously identified 262 differentially alternatively spliced exons during neuronal cell differentiations, we extracted 241 sets that changed the amino acid sequences between the alternatively spliced sequences. Conserved domain searches indicated that annotated domain(s) were changed in 128 sets. We obtained 49 genes whose terms overlapped between domain description and gene annotation. Thus, these 49 genes have alternatively differentially spliced in exons that affect their main functions. We performed pathway analysis using these 49 genes and identified the EGFR (epidermal growth factor receptor) and mTOR (mammalian target of rapamycin) signaling pathway as being involved frequently. Recent studies reported that the mTOR pathway is associated with neuronal cell differentiation, vindicating that our approach extracted an important molecular network successfully.
Effective informatics approaches for exons should be more complex than those for genes, because changes in alternative exons affect protein functions via alterations of amino acid sequences and functional domains. Our method extracted alterations of functional domains and identified key alternative splicing events. We identified the EGFR and mTOR signaling pathway as the most affected pathway. The mTOR pathway is important for neuronal differentiation, suggesting that this in silico extraction of alternative splicing networks is useful. This preliminary analysis indicated that automated analysis of the effects of alternative splicing would provide a rich source of biologically relevant information.
Approximately 23,000 human protein coding genes have been identified; however, this is a much smaller number than the expected over 200,000 human proteins . Alternative splicing changes the use of exons, producing multiple transcripts from a single gene, and enhances proteomic diversity to support complexity in higher eukaryotes . Indeed, it was reported that approximately 90% of human genes undergo alternative splicing . There are five basic models of alternative splicing: exon skipping, mutually exclusive type, 5’ splice site selection, 3’ splice site selection, and intron retention. Additionally, alternative promoters and alternative polyadenylation sites can produce alternative isoforms of transcripts. These exons are sometimes continuously located on the genome, resulting in complex alternative splicing .
According to physiological or environmental changes, some alternative splicing occurs in a spatiotemporal manner, regulated by alternative splicing regulators. Several regulators, such as Nova1, RBFox1 and nPTB, and their targets, have been identified in neural tissues and cells [5–7]. Many isoforms produced by alternative splicing have distinctly different functions, and play important biological roles [8, 9]. Therefore, it is important to determine critical alternative splicing networks based on biological phenomena.
Over the last decade, many high throughput analyses have been performed, and their information has accumulated in databases. Both gene level expression and exon level expression are available from analyses using large-scale sequencing technologies. Furthermore, recently developed standard microarrays may provide information on both gene and exon expressions . Certain differentially expressed genes significantly affect biological phenomena or are useful as molecular markers; therefore, information on gene level expressions is generally well analyzed. Moreover, critical gene networks or target genes can be identified from gene expression information. However, this is not the case for exons.
Complicated species of alternative isoforms with small numbers of nucleotide changes can be predicted compared with genes. Thus, the quantitative credibility of an exon’s information is generally much less than that of genes. Moreover, the analysis of exons is more complex than that of genes. Although some differentially alternatively spliced exons have been investigated and extracted, the annotated genes of these exons may be used to examine the critical networks or targets [11–13]. However, significant changes resulting from alternative exons may occur in protein domain(s). Advanced molecular dynamics techniques were also applied to investigate alternative isoforms of particular individual isoforms [14–18]. These methods are not currently applicable to the genome-wide investigation of alternative isoforms at a time. Of course, several specialized programs support informatics analyses of alternative exons [19–22]. However, these do not support analyses of the protein domains encoded on alternative exons. In the case of individual small-scale analysis, researchers generally check alterations of amino acid sequences and of functional protein domains according to usages of alternative exons. AltAnalyze and DomainGraph could provide protein domain information from alternative exons although this type of analysis of exons has been rarely applied comprehensively .
Previously, we analyzed comprehensively the differentially alternatively spliced exons during neuronal differentiation of P19 mouse embryonic carcinoma cells . Validation by reverse-transcription polymerase chain reaction (RT-PCR) suggested that 87% of the obtained 262 exons were differentially alternatively spliced in neuronal cells compared with undifferentiated cells. Moreover, many of the genes that were annotated by 262 exons were suggested to be involved in neural events. Thus, the 262 exons were plausible as neural splicing exons, and we considered that these exons could be a good example group to investigate alterative splicing networks.
In this article, we searched for a network that involved these alternative splicing events using functional domain information. Although individual studies have focused on the functionalities of alternative spliced domains, all comprehensively obtained exons were used as objects in the domain functionality search. Ultimately, we identified the EGFR (epidermal growth factor receptor) and mTOR (mammalian target of rapamycin) signaling pathway as the molecular network most associated with neural alternative splicing. The EGFR/mTOR pathway is an important cellular signaling pathway that controls cell growth and proliferation [16, 24]. This pathway has been intensively studied, and many of the proteins involved have been identified. Recently, studies have found that the mTOR pathway plays important roles in the maintenance of neural stem cells and the differentiation of neuronal cells [reviewed in . However, the relationship between the EGFR/mTOR pathway and alternative splicing remains unknown. Our investigation suggested that gene expression and exon expression of the transcripts involved in this pathway were dramatically changed during neural differentiation of P19 cells. Indeed, our trial showed that searching for molecular network according to the functionalities of alternatively spliced domains was an effective strategy. Thus, it will help to precisely analyze mega data that are available from comprehensive analyses.
Data set of the exon array during neuronal differentiation of P19 cells
The accession number of the exon Array data is GSE23710, which we analyzed and reported previously . Information on the DAS exons, such as probeset sequences, is available in the supplemental table of that report , as is the information on the differentially expressed genes.
Collection of annotated sequence data from the probeset sequences
All probeset sequences were analyzed in the UCSC Blast-like alignment (blat) search tool (http://genome.ucsc.edu/cgi-bin/hgBlat) , and we determined alternative exons or region sequences compared with the annotated Refseq sequence, mRNA sequence or EST sequences. Similarly, we manually selected the most typical and representative sequences in which the alternative exon (or region) was either joined to or excised from a transcript. All transcript sequences were translated into amino acid sequences using the ExPASy translate tool (Swiss Institute of Bioinformatics). Sequences of the determined DAS exons were also translated into amino acid sequences. Whole amino acid sequences, excluding or including alternative exons or regions, were compared and validated in the UCSC blat search.
Conserved domain search for alternatively spliced isoforms
The obtained amino acid sequences were analyzed in the NCBI conserved domain search , and alternatively spliced domains were determined. Additionally, text descriptions of these domains were retrieved.
Gene ontology analysis, pathway analysis and text mining
The GO analyses were performed for 128 DAS genes with altered domain(s) . The GO term(s) were compared with the text descriptions of the domain(s). Overlaps between GO terms and domain descriptions were found in 49 out of 128 genes. These 49 genes were subjected to pathway analysis and statistically-significant (Fisher’s Exact Test p ≤ 0.01) biological process terms were obtained using PathwayStudio® (Ariadne Genomics Inc., Rockville, MD, USA) [28, 29]. Although the schematic representation of pathway analysis was based on this analysis, their relationships were validated by KEGG pathway analysis (http://www.genome.jp/kegg/pathway.html) and by previous articles for Arap1, Ep400 and Arhgef12 [30–32].
Cell culture and RNA purification
P19 cells were maintained in α minimum essential medium (α-MEM; Sigma–Aldrich, St. Louis, MO, USA) supplemented with 10% fetal bovine serum (FBS; Sigma–Aldrich) . To induce neuronal cell differentiation, P19 cells (1 × 105 cells/mL) were treated with 1 μM all trans-retinoic acid at 4°C in a 10 cm petri dish (Falcon) with α-MEM containing 10% FBS, as described previously [33, 34]. Total RNAs were collected from undifferentiated cells (Day 0), cells during neuronal induction (Day 1–4), differentiation (Day 5–9), and from cells in the early glial stage (Day 10–12) using the Trizol reagent (Life Technologies, Grand Island, NY, USA).
The total RNAs of P19 cells were prepared as described above. The cDNAs were produced from 2 μg of total RNAs using SuperScript III (Life technologies) and 0.5 μg of oligo dT primer in a 20 μL reaction mixture. The cDNAs are used as templates for PCR. GoTaq polymerase (Promega, Fitchburg, WI, USA) with specific primers performed the PCR reactions. The DNA primer sequences, number of cycles and annealing temperatures for the candidates are described in Additional file 1: Table S1. Primers and conditions for β-Actin and GluR1 (controls) were described in a previous report . The PCR products were analyzed on a 6% polyacrylamide gel. The gels were stained with SYBR Green I (Takara Bio Inc., Otsu, Japan), and a LAS-3000 (GE Healthcare, Fairfield, CT, USA) was used to analyze the images. The sequences of the PCR products were confirmed in a 3100 DNA sequencer (Life technologies). MultiGauge v 3.0 software (GE Healthcare) was used to perform the densitometry. Each experiment was performed at least three times to confirm reproducibility.
As for Ethics, this research did not involve any human subject, human material, or human data, and was not performed on any animals. This research involved Recombinant DNA Experiments and was approved by Life science committee of Japan Advanced Institute of Science and Technology.
Results and discussion
Extraction of important genes regulated by alternative splicing
Previously, we obtained 262 candidate exons (236 genes) that were differentially alternatively spliced during neuronal differentiation of P19 cells . Validation by RT-PCR indicated that the expressions of 87% of them actually changed; therefore, we used the 262 exons as good example for subsequent analysis (Figure 1). This was performed in a mouse Exon Array; therefore, intact information of these exons was provided as probeset sequences (probeset seqs). These probeset seqs were analyzed in a genome viewer; human blat analysis. According to annotated Refseqs, mRNAs or ESTs, the 262 alternative exons were analyzed and the two entire amino acid sequences that the alternative exon candidate joined or excised in its transcripts were determined (Figure 1, Additional file 2: Table S2). Occasionally, several Refseqs were annotated at the genome locus. In these cases, we selected the most plausible and representative transcript that involved in the alternative exon. Although 81 probeset seqs were located in the untranslated regions (UTRs), 52 probeset seqs indicated exons that contained protein-coding sequences (Figure 2A). Moreover, some alternative terminations by non-coding UTR exons could cause changes in the amino acid sequence. Altogether, 241 entire transcript sequences producing changes of amino acid sequence by including and excluding the alternative exon were identified (Figures 1 and 2A, Additional file 2: Table S2).To clarify the effects on functional domains of the alternative exons, we performed a conserved domain search at NCBI. Among 241 variable coding regions, 128 might cause changes to the annotated domains (Figures 1 and 2A). The others were suggested to affect N-terminal regions, C-terminal regions and linker region (inter-domains) (Figure 2A). In the case of alternative promoters or polyadenylation, their upstream or downstream exons could be excluded from its transcript, respectively. Therefore, some alternative exons are responsible for longer amino acid sequences and affect large parts of the amino acid sequence. These case are shown in the upper right area of Figure 2B and these drastic changes of amino acid sequences mostly caused functional domain changes (blue dots). In addition, many alternative splicings changed < 100 amino acids and < 20% of the full-length protein, as shown in lower left area of Figure 2B. Some of them also produced alterations in functional protein domains (Figure 2B).For the 128 regions that changing their domains, we extracted the descriptions of their domains from the NCBI conserved domain search. We also performed Gene Ontology (GO) analysis of the corresponding 128 genes. We compared the domain descriptions with the GO terms, and found 49 genes (54 probeset seqs, including five probeset seqs that indicated the same exons and/or domains) that had words in common, not including conjunctions and prepositions (Figures 1 and 2A). Thus, the critically important molecular functions of 49 genes could be affected or regulated by their differentially alternatively spliced exons.
Pathway analysis based on the domain change
To find molecular networks or pathways that involved the differentially alternatively spliced exons, we subjected the 49 genes to pathway analysis. As a result, we identified the EGFR/mTOR signaling pathway as the most significantly linked network by these alternative exons (Figures 1 and 3). There were seven alternative splicing events in this pathway. For example, the alternative exon of Egfr (epidermal growth factor receptor) suggested that the extracellular protein isoform was expressed via alternative polyadenylation located in an adjacent upstream exon that encoded the transmembrane domain . Although the expected extracellular isoform of Egfr is interesting, the distinct functions of this isoform are unclear. The alternative exon of mTor (mammalian target of rapamycin) is part of the N-terminal predicted short protein isoform predictively, which has been reported to be regulated by an alternative splicing regulator, SAM68, during adipogenesis. It was also suggested to be involved in nonsense mediated mRNA decay . Meanwhile, Tsc2 (tuberous sclerosis 2) was in the original 262 candidate exons but was eliminated from the 49 exons (Figure 3)  because the region affected by the alternative exon did not encode a functional domain.
Generally, differentially expressed genes affect biological phenomena. Thus, we also checked and extracted the differentially expressed genes associated with the EGFR/mTOR pathway, and assessed with their increased or decreased expression (Figure 3) . Additionally, it was reported that protein expression of p53 could increase during the neuronal differentiation of P19 cells . Including gene expressions and exons expression, various transcripts were observed to have altered in this pathway. However, it is speculated that upregulation of Egfr gene expression increases the potency of its signal, while increased expression of its alternative exon interferes with the signal via its imperfect extracellular isoform.
Validations of gene and exon expressions in EGFR/mTOR pathway
Speculations based on gene expressions or exon expressions may contradict each other. We suspected that some comprehensive deposited data include partial or uncertain information; we confirmed experimentally the gene and exon expressions during neuronal differentiation of P19 cells. First, we examined the expression of alternatively spliced exons of Egfr using semi-quantitative RT-PCR. The proportions of the transcript that encodes extracellular isoform and full-length products changed from about 1.5:8.5 to 4:6 during neuronal differentiation of P19 cells (Figure 4). This was suggested by the splicing index value (SI) in our previous study . The loss of domain type, such as an extracellular isoform, increased; therefore, we marked EGFR the inner red color on Figure 3. There is no experimental evidence of extracellular isoform’s functionality compared with the full-length product; therefore, it is difficult to determine the functional changes induced in these isoforms. In addition, the amounts of extracellular and full-length type transcripts increased by approximately 6.8-fold and 1.7-fold, respectively (Figure 4). This change also correlated with the change that was suggested by fold change (FC) in neuronal gene expression in a previous comprehensive analysis . We indicated increased gene-expression of EGFR with the outer blue color on Figure 3.
In the case of another key factor in this network, the SI value suggested that the N-terminal short isoform-encoding transcript of mTor increased during neuronal differentiation, and the FC value suggested that expression of the gene did not increase . Actually, the relative amounts of N-terminal isoform products did increase and the amount of total transcripts from the mTor gene did not remarkably change during differentiation (Figure 4). We indicated the increased loss of domain type, such as N-terminal short isoform of mTor, with the inner red color in Figure 3. Although the function of the mTor N-terminal isoform is unclear during differentiation, its transcript is degraded during adipogenesis . Perhaps a truncated isoform that loses a functional protein domain loses its function. As mentioned above, mTOR signaling promotes cell proliferation [24, 38], and the N-terminal isoform lacks certain domains present in the full-length protein. Cell differentiation and cell proliferation are generally contrary phenomena. Therefore, this N-terminal alternative spliced isoform of mTor may be important for the neuronal differentiation of P19 cells.
Although extracellular isoform of Egfr may negatively affect this pathway, similar to the mTor N-terminal isoform, the full-length type Egfr transcripts increased. Therefore, we tested the expressions of genes and their exons that exist between EGFR and mTOR in the pathway and that were suggested to change by their SI and FC values. Typically, the Akt1 (RAC-alpha serine/threonine-protein kinase) transcript increased 7-fold during neuronal differentiation (Figure 4). The change of Akt1 alternative isoforms was not suggested by its SI value. The decrease of gene expression of Tsc1 was suggested by its FC and was observed to change slightly during differentiation (Figure 4). Although the suggested change of Tsc2 was validated, the effect of its alternative splicing is unclear. Another suggested downstream gene, Shc1 (SHC-transforming protein 1), was validated and decreased during differentiation (Figure 4). Thus, gene expressions of Egfr and Akt1 increased, but the expressions of other genes in this pathway did not. Besides mTor or Tscs, various target proteins of Akt1 have been identified [reviewed in . We speculated that other target proteins may be dramatically upregulated by enhanced Akt1 potency, according to their gene expressions. In the case of mTor, its alternative N-terminal product may slightly negatively modulate the enhanced input signals.
Meanwhile, some upstream genes of the EGFR pathway, such as Hbegf (Heparin-binding EGF-like growth factor) and Arhgef12 (Rho guanine nucleotide exchange factor 12) increased (Figures 3 and 4). The alternative isoform of Arap1 (ArfGAP with RhoGAP domain, ankyrin repeat and PH domain 1) was not detected; however, the expression of its full-length product decreased during differentiation (Figure 4). In the case of Ep400 (E1A binding protein p400), its expression increased and the alternatively spliced product may positively affect EGFR pathway according to its changed domain. Additionally, changes in indirectly related genes, such as Kitlg (KIT Ligand) and Cdkn2a (cyclin-dependent kinase inhibitor 2A) were also validated their gene- and exon- expressions during differentiation (Figures 3 and 4). However, the functional change caused by the domain modification by alternative splicing of Kitlg was unclear. It was difficult to determine the effects on EGFR signaling because upstream genes showed positive and negative changes.
Time-course analysis of gene and exon expressions of Egfr
Autonomous regulation represents one possible explanation of the contradictory expressions of the Egfr gene and exons. For example, increased gene expression indirectly affects its alternative splicing to repress gene function. To test this hypothesis, we prepared a series of cDNAs from the undifferentiated stage to the glial cell stage of P19 cells. We then examined the most changed time points of Egfr gene or exon expression. The relative amount of the extracellular Egfr isoform increased in the day 1 sample during neural induction by retinoic acid (RA) under aggregation culture (Figure 5). Similar fast responses by alternative splicing have been observed in other genes during differentiation [40, 41]. The previously reported genes were also regulated by RA and/or cell aggregation. At this time point, the total amount of Egfr transcript showed little change (Figure 5). The proportion of the extracellular isoform then decreased, before rising again from day 4 to day 9: from neural induction to the neuronal cell stage of P19 cells (Figure 5). On the other hand, increased Egfr gene expression was observed from day 3 to day 8 (Figure 5). In other words, the elevation of Egfr gene expression was started 1 day before the proportional increase of its extracellular isoform and stopped 1 day before the isoform began to decrease again. Similarly, we tested mTor using a series of cDNAs. Its N-terminal isoform was elevated during the neural induction and differentiation stages, which was completely different from Egfr. Although the expression of the Egfr gene and its alternative isoforms expressions fitted with the autonomous regulation hypothesis and may rebalance their functions, we lack proof of the hypothesis. This should be addressed in our future work.
Changes in alternative exons affect protein functions via alterations of amino acid sequences and functional domains. Our method extracted alterations of functional domains and identified key alternative splicing events. We found that the EGFR/mTOR signaling pathway was the most affected pathway. The importance of mTOR in neuronal differentiation has been reported, suggesting that this in silico extraction of alternative splicing network is a useful strategy. This strategy for the analysis of alternative splicing should be automated. The experimental validations of exon and gene expression by RT-PCR suggested that the increase in the products of genes such as Egfr and Akt1 might increase the signal through the pathway during neuronal differentiation of P19 cells, meanwhile the alternative splicing events in mTor might control this pathway repressively. We speculate that our method will contribute to future studies of new molecular networks of alternative splicing regulation.
SA and HTTP are graduate students, and MT, TT and HS are faculties in Japan Advanced Institute of Science and Technology. The corresponding author is HS, who also belong to Center for Nano Materials and Technology, Japan Advanced Institute of Science and Technology. MO was a summer-time internship student in Japan Advanced Institute of Science and Technology and belongs to Department of Chemicals and Engineering, Miyakonojo National College of Technology. KK is a president of World Fusion Co., Ltd.
- P19 cells:
P19 mouse embryonic carcinoma cells
Differentially alternatively spliced
Epidermal growth factor receptor
Mammalian target of rapamycin
RAC-alpha serine/threonine-protein kinase
SHC-transforming protein 1
Heparin-binding EGF-like growth factor
Rho guanine nucleotide exchange factor 12
Cyclin-dependent kinase inhibitor 2A
ArfGAP with RhoGAP domain ankyrin repeat and PH domain 1
E1A binding protein p400.
International Human Genome Sequencing Consortium: Finishing the euchromatic sequence of the human genome. Nature. 2004, 431: 931-945. 10.1038/nature03001.
Graveley BR: Alternative splicing: increasing diversity in the proteomic world. Trends Genet. 2001, 17: 100-107. 10.1016/S0168-9525(00)02176-4.
Wang ET, Sandberg R, Luo S, Khrebtukova I, Zhang L, Mayr C, Kingsmore SF, Schroth GP, Burge CB: Alternative isoform regulation in human tissue transcriptomes. Nature. 2008, 456: 470-476. 10.1038/nature07509.
Park JW, Graveley BR: Complex alternative splicing: alternative splicing in the postgenomic Era. Advances in Experimental Medicine and Biology, Volume 623. Edited by: Blencowe BJ, Graveley BR. 2007, New York: LANDES Bioscience, 50-63.
Coutinho-Mansfield GC, Xue Y, Zhang Y, Fu XD: PTB/nPTB switch: a post-transcriptional mechanism for programming neuronal differentiation. Genes Dev. 2007, 21: 1573-1577. 10.1101/gad.1575607.
Ule J, Stefani G, Mele A, Ruggiu M, Wang X, Taneri B, Gaasterland T, Blencowe BJ, Darnell RB: An RNA map predicting Nova-dependent splicing regulation. Nature. 2006, 444: 580-586. 10.1038/nature05304.
Zhang C, Zhang Z, Castle J, Sun S, Johnson J, Krainer AR, Zhang MQ: Defining the regulatory network of the tissue-specific splicing factors Fox-1 and Fox-2. Genes Dev. 2008, 22: 2550-2563. 10.1101/gad.1703108.
Black DL: Mechanisms of alternative pre-messenger RNA splicing. Annu Rev Biochem. 2003, 72: 291-336. 10.1146/annurev.biochem.72.121801.161720.
Matlin AJ, Clark F, Smith CW: Understanding alternative splicing: towards a cellular code. Nat Rev Mol Cell Biol. 2005, 6: 386-398. 10.1038/nrm1645.
Hallegger M, Llorian M, Smith CW: Alternative splicing: global insights. FEBS J. 2010, 277: 856-866. 10.1111/j.1742-4658.2009.07521.x.
Weyn-Vanhentenryck SM, Mele A, Yan Q, Sun S, Farny N, Zhang Z, Xue C, Herre M, Silver PA, Zhang MQ, Krainer AR, Darnell RB, Zhang C: HITS-CLIP and integrative modeling define the rbfox splicing-regulatory network linked to brain development and autism. Cell Rep. 2014, 6: 1139-1152. 10.1016/j.celrep.2014.02.005.
Suzuki H, Osaki K, Sano K, Alam AH, Nakamura Y, Ishigaki Y, Kawahara K, Tsukahara T: Comprehensive analysis of alternative splicing and functionality in neuronal differentiation of P19 cells. PLoS One. 2011, 6: e16880-10.1371/journal.pone.0016880.http://www.plosone.org/article/info%3Adoi%2F10.1371%2Fjournal.pone.0016880,
Voineagu I, Wang X, Johnston P, Lowe JK, Tian Y, Horvath S, Mill J, Cantor RM, Blencowe BJ, Geschwind DH: Transcriptomic analysis of autistic brain reveals convergent molecular pathology. Nature. 2011, 474: 380-384. 10.1038/nature10110.
Reddy PP, Raghuram V, Hradsky J, Spilker C, Chakraborty A, Sharma Y, Mikhaylova M, Kreutz MR: Molecular dynamics of the neuronal EF-hand Ca2+-sensor Caldendrin. PLoS One. 2014, 9: e103186-10.1371/journal.pone.0103186.
Soni S, Tyagi C, Grover A, Goswami SK: Molecular modeling and molecular dynamics simulations based structural analysis of the SG2NA protein variants. BMC Res Notes. 2014, 7: 446-10.1186/1756-0500-7-446.
Barbany M, Morata J, Meyer T, Lois S, Orozco M, de la Cruz X: Characterization of the impact of alternative splicing on protein dynamics: the cases of glutathione S-transferase and ectodysplasin-a isoforms. Proteins. 2012, 80: 2235-2249. 10.1002/prot.24112.
Rajendran V, Purohit R, Sethumadhavan R: In silico investigation of molecular mechanism of laminopathy caused by a point mutation (R482W) in lamin A/C protein. Amino Acids. 2012, 43: 603-615. 10.1007/s00726-011-1108-7.
Suzuki H, Takeuchi M, Sugiyama A, Alam AH, Vu LT, Sekiyama Y, Dam HC, Ohki SY, Tsukahara T: Alternative splicing produces structural and functional changes in CUGBP2. BMC Biochem. 2012, 13: 6-10.1186/1471-2091-13-6.
Chang TY, Li YY, Jen CH, Yang TP, Lin CH, Hsu MT, Wang HW: easyExon—a Java-based GUI tool for processing and visualization of Affymetrix exon array data. BMC Bioinforma. 2008, 9: 432-10.1186/1471-2105-9-432.
Xing Y, Stoilov P, Kapur K, Han A, Jiang H, Shen S, Black DL, Wong WH: MADS: a new and improved method for analysis of differential alternative splicing by exon-tiling microarrays. RNA. 2008, 14: 1470-1479. 10.1261/rna.1070208.
Yates T, Okoniewski MJ, Miller CJ: X:Map: annotation and visualization of genome structure for affymetrix exon array analysis. Nucleic Acids Res. 2008, 36: D780-D786.
Purdom E, Simpson KM, Robinson MD, Conboy JG, Lapuk AV, Speed TP: FIRMA: a method for detection of alternative splicing from exon array data. Bioinformatics. 2008, 24: 1707-1714. 10.1093/bioinformatics/btn284.
Emig D, Salomonis N, Baumbach J, Lengauer T, Conklin BR, Albrecht M: AltAnalyze and DomainGraph: analyzing and visualizing exon expression data. Nucleic Acids Res. 2010, Web Server issue: W755-W762.
Wullschleger S, Loewith R, Hall MN: TOR signaling in growth and metabolism. Cell. 2006, 124: 471-484. 10.1016/j.cell.2006.01.016.
Magri L, Galli R: mTOR signaling in neural stem cells: from basic biology to disease. Cell Mol Life Sci. 2013, 70: 2887-2898. 10.1007/s00018-012-1196-x.
Kuhn RM, Karolchik D, Zweig AS, Wang T, Smith KE, Rosenbloom KR, Rhead B, Raney BJ, Pohl A, Pheasant M, Meyer L, Hsu F, Hinrichs AS, Harte RA, Giardine B, Fujita P, Diekhans M, Dreszer T, Clawson H, Barber GP, Haussler D, Kent WJ: The UCSC genome browser database: update 2009. Nucleic Acids Res. 2009, 37 (Database issue): D755-D761.
Marchler-Bauer A, Lu S, Anderson JB, Chitsaz F, Derbyshire MK, DeWeese-Scott C, Fong JH, Geer LY, Geer RC, Gonzales NR, Gwadz M, Hurwitz DI, Jackson JD, Ke Z, Lanczycki CJ, Lu F, Marchler GH, Mullokandov M, Omelchenko MV, Robertson CL, Song JS, Thanki N, Yamashita RA, Zhang D, Zhang N, Zheng C, Bryant SH: CDD: a conserved domain database for the functional annotation of proteins. Nucleic Acids Res. 2011, 39 (Database issue): D225-D229.
Nikitin A, Egorov S, Daraselia N, Mazo I: Pathway studio–the analysis and navigation of molecular networks. Bioinformatics. 2003, 19: 2155-2157. 10.1093/bioinformatics/btg290.
Novichkova S, Egorov S, Daraselia N: MedScan, a natural language processing engine for MEDLINE abstracts. Bioinformatics. 2003, 19: 1699-1706. 10.1093/bioinformatics/btg207.
Yoon HY, Kales SC, Luo R, Lipkowitz S, Randazzo PA: ARAP1 association with CIN85 affects epidermal growth factor receptor endocytic trafficking. Biol Cell. 2011, 103: 171-184. 10.1042/BC20100154.
Flinterman MB, Mymryk JS, Klanrit P, Yousef AF, Lowe SW, Caldas C, Gäken J, Farzaneh F, Tavassoli M: p400 function is required for the adenovirus E1A-mediated suppression of EGFR and tumour cell killing. Oncogene. 2007, 26: 6863-6874. 10.1038/sj.onc.1210497.
Bourguignon LY, Gilad E, Brightman A, Diedrich F, Singleton P: Hyaluronan-CD44 interaction with leukemia-associated RhoGEF and epidermal growth factor receptor promotes Rho/Ras co-activation, phospholipase C epsilon-Ca2+ signaling, and cytoskeleton modification in head and neck squamous cell carcinoma cells. J Biol Chem. 2006, 281: 14026-14040. 10.1074/jbc.M507734200.
Jones-Villeneuve EM, McBurney MW, Rogers KA, Kalnins VI: Retinoic acid induces embryonal carcinoma cells to differentiate into neurons and glial cells. J Cell Biol. 1982, 94: 253-262. 10.1083/jcb.94.2.253.
Hakim NH, Kounish T, Alam AHMK, Tsukahara T, Suzuki H: Alternative splicing of Mef2c promoted by Fox-1 during neural differentiation in P19 cells. Genes Cells. 2010, 15: 255-267. 10.1111/j.1365-2443.2009.01378.x.
Sharma SV, Bell DW, Settleman J, Haber DA: Epidermal growth factor receptor mutations in lung cancer. Nat Rev Cancer. 2007, 7: 169-181. 10.1038/nrc2088.
Huot MÉ, Vogel G, Zabarauskas A, Ngo CT, Coulombe-Huntington J, Majewski J, Richard S: The Sam68 STAR RNA-binding protein regulates mTOR alternative splicing during adipogenesis. Mol Cell. 2012, 46: 187-199. 10.1016/j.molcel.2012.02.007.
Ostrakhovitch EA, Semenikhin OA: p53-mediated regulation of neuronal differentiation via regulation of dual oxidase maturation factor 1. Neurosci Lett. 2011, 494: 80-85. 10.1016/j.neulet.2011.02.061.
Jozwiak J, Jozwiak S, Grzela T, Lazarczyk M: Positive and negative regulation of TSC2 activity and its effects on downstream effectors of the mTOR pathway. Neuromol Med. 2005, 7: 287-296. 10.1385/NMM:7:4:287.
Emamian ES: AKT/GSK3 signaling pathway and schizophrenia. Front Mol Neurosci. 2012, 5: 33-
Alam AH, Suzuki H, Tsukahara T: Retinoic acid treatment and cell aggregation independently regulate alternative splicing in P19 cells during neural differentiation. Cell Biol Int. 2010, 34: 631-643. 10.1042/CBI20090332.
Alam AH, Suzuki H, Tsukahara T: Expression analysis of Fgf8a &Fgf8b in early stage of P19 cells during neural differentiation. Cell Biol Int. 2009, 33: 1032-1037. 10.1016/j.cellbi.2009.06.015.
We thank Mr. K. Osaki for his helpful suggestions. We acknowledge Edanz Group Ltd. for their English editing service. HS was supported in part by a Grant-in-Aid for Young Scientists (B) from the Japan Society for the Promotion of Science (JSPS; 24700974), and in part by an Intramural Research Grant (25–5) for Neurological and Psychiatric Disorders of National Center of Neurology and Psychiatry. TT was supported in part by a Grant-in-Aid for Scientific Research (B) from the JSPS.
Although KK is affiliated with a company, his contribution in this article was purely academic. We declare that we have no competing interests.
HS and AS carried designed the study, performed certain informatics analyses of alternative spliced products, helped to perform the RT-PCR experiments and drafted the manuscript. HTTP and MO carried out the RT-PCR experiments. KK evaluated the GO analysis and pathway analysis. MT and TT helped to draft the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 2: Table S2: Results of the analyses in 262 DAS exons. Using filter conditions of the last column, 49 gene with 54 probesets are listed. (XLSX 255 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Alam, S., Phan, H.T.T., Okazaki, M. et al. Computational extraction of a neural molecular network through alternative splicing. BMC Res Notes 7, 934 (2014). https://doi.org/10.1186/1756-0500-7-934
- Comprehensive analysis
- Neuronal differentiation
- Alternative splicing