Next Article in Journal
STIM Protein-NMDA2 Receptor Interaction Decreases NMDA-Dependent Calcium Levels in Cortical Neurons
Next Article in Special Issue
Poison-Exon Inclusion in DHX9 Reduces Its Expression and Sensitizes Ewing Sarcoma Cells to Chemotherapeutic Treatment
Previous Article in Journal
Increased O-GlcNAcylation of c-Myc Promotes Pre-B Cell Proliferation
Previous Article in Special Issue
An Intricate Connection between Alternative Splicing and Phenotypic Plasticity in Development and Cancer
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Cells
Volume 9
Issue 1
10.3390/cells9010159
cells-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Splicing Players Are Differently Expressed in Sporadic Amyotrophic Lateral Sclerosis Molecular Clusters and Brain Regions

by
Valentina La Cognata
1,
Giulia Gentile
1,
Eleonora Aronica
2 and
Sebastiano Cavallaro
1,*
1
Institute for Biomedical Research and Innovation, National Research Council, 95126 Catania, Italy
2
Department of (Neuro)Pathology, Amsterdam UMC, University of Amsterdam, Amsterdam Neuroscience, Meibergdreef 9, 1105 Amsterdam, The Netherlands
*
Author to whom correspondence should be addressed.
Cells 2020, 9(1), 159; https://doi.org/10.3390/cells9010159
Submission received: 26 November 2019 / Revised: 23 December 2019 / Accepted: 4 January 2020 / Published: 8 January 2020
(This article belongs to the Special Issue Alternative Splicing: Recent Insights into Mechanisms and Functional Roles)
Browse Figures
Versions Notes

Abstract

:
Splicing is a tightly orchestrated process by which the brain produces protein diversity over time and space. While this process specializes and diversifies neurons, its deregulation may be responsible for their selective degeneration. In amyotrophic lateral sclerosis (ALS), splicing defects have been investigated at the singular gene level without considering the higher-order level, involving the entire splicing machinery. In this study, we analyzed the complete spectrum (396) of genes encoding splicing factors in the motor cortex (41) and spinal cord (40) samples from control and sporadic ALS (SALS) patients. A substantial number of genes (184) displayed significant expression changes in tissue types or disease states, were implicated in distinct splicing complexes and showed different topological hierarchical roles based on protein–protein interactions. The deregulation of one of these splicing factors has a central topological role, i.e., the transcription factor YBX1, which might also have an impact on stress granule formation, a pathological marker associated with ALS.

Graphical Abstract

1. Introduction

Amyotrophic lateral sclerosis (ALS) is a progressive and fatal disease of motor function, characterized by the selective degeneration of upper and lower motor neurons. Alternative splicing plays an important role in neuronal differentiation, and the dysregulation of RNA processing and splicing is a key contributor to ALS pathogenesis [1,2]. Indeed, not only do a vast number of alternative splicing defects involve ALS associated genes, but several of their encoded products participate in RNA processing and splicing (FUS, TDP-43, C9ORF72, ATXN2, TAF15, hnRNPA1, hnRNPA2B1, MATR3, EWSR1, SETX) [1,2].
RNA splicing is one of the major drivers of eukaryotic RNA metabolism, allowing the removal of introns and the joining of exons from precursor mRNA in a constitutive or alternative manner. Nearly all human multi-exon genes undergo splicing, but this process is particularly widespread in the brain, where splicing programs allow for cellular protein diversity [3,4,5]. Although RNA splicing represents an extremely economical means for increasing protein diversity, in the terms of both time and space, finely tuning the genomic information to meet the unique needs of each neuron, it adds an additional layer of complexity and vulnerability [6]. In fact, while this process specializes and diversifies neurons, its deregulation may be responsible for their selective degeneration.
To date, splicing defects in ALS have been investigated at the singular gene level without considering the higher-order level involving the splicing machinery [2]. Indeed, RNA splicing is mediated by the spliceosome, a dynamic macromolecular machinery comprised of more than 100 individual components [7]. Unlike many enzymes, the spliceosome does not have a preformed catalytic center. Instead, it is constructed de novo on the pre-mRNA in a stepwise manner for each round of splicing and undergoes multiple rearrangements to configure the catalytic core [8]. The major spliceosome machine is composed of five small nuclear ribonucleoprotein complexes, (snRNPs), U1, U2, U4, U5, and U6, each containing a cognate U-rich noncoding small nuclear RNA (snRNA). A small subset of unusual introns containing non-consensus termini are processed by a much less abundant spliceosomal machine (minor spliceosome) which contains the U5 snRNP and analogs of the remaining snRNPs (U11, U12, U4atac, and U6atac) [8,9,10,11,12]. Following transcription, snRNAs are exported through the nuclear pore complex (NPC) into the cytoplasm where snRNPs maturate inside a large assemblysome (the survival motor neuron (SMN) complex); core snRNPs are then imported into the nucleus, assembled in Cajal bodies and are ready to be recruited by the nascent pre-mRNA for splicing [8,13]. The basepairing of snRNPs with pre-mRNA relies on a series of conserved regulatory sequence elements (the splice donor site, the splice acceptor site, the branch site, and the polypyrimidine tract) [14]. Additional layers of regulation are introduced by some auxiliary motifs (exon and intron splicing enhancers or silencers) that may stimulate exon inclusion or inhibit splicing with the mediation of RNA-binding proteins (RBPs), mainly of which belong to the families of serine/arginine-rich splicing factor (SRSF) proteins and heterogeneous nuclear ribonucleoproteins (hnRNPs), acting as splicing enhancers or repressors, respectively [15,16]. For a complete overview of splicing machinery, intermediated complexes and auxiliary regulatory factors, readers are invited to refer to some complete reviews [8,17,18,19].
To investigate the role of the spliceosome in ALS, here we analyzed the entire spectrum (396) of genes encoding splicing factors in the motor cortex and spinal cord samples from control and sporadic ALS (SALS) patients. These genes encode both the core spliceosomal components involved in constitutive splicing, as well as hnRNPs (usually splice site inhibitors), SRSFs (usually splice site activators) and additional splicing factor regulators. A substantial number of genes displayed significant expression changes in tissue types or disease states and showed different topological hierarchical roles based on protein–protein network interactions. One of these, the transcription factor YBX1, might have an impact on stress granule formation, pathological markers associated with ALS.

2. Materials and Methods

2.1. Data Source and Gene Selection

For this study, we referred to two transcriptome data sets, both available at ArrayExpress (http://www.ebi.ac.uk/arrayexpress/) with the accession numbers E-MTAB-2325 (https://www.ebi.ac.uk/arrayexpress/experiments/E-MTAB-2325/) (previously described as reported in [20]) and E-MTAB-8635 (https://www.ebi.ac.uk/arrayexpress/experiments/E-MTAB-8635/) (new submission). These data sets consisted of the mRNA expression profiles of 41,059 genes from 41 motor cortex samples and 40 spinal cord samples of SALS and control subjects (Table 1) hybridized onto 4 × 44K Whole Human Genome Oligo expression microarrays (Agilent Technologies, Turin, Italy). A detailed description of the study design, subject characteristics and experimental procedures was reported in the original publication [20]. In order to identify the common or differential genomic changes between the spinal cords and motor cortexes belonging to the same cohort of patients, the sample preparation, RNA extraction and microarray processing of the spinal cord tissues obtained from patients of the same cohort were performed as previously described [20]. The transcriptomic data obtained from the spinal cord were combined with data from the motor cortex to create a complete dataset of 81 whole transcriptomic observations, as listed in Table 1. Once obtained, the raw intensity values were thresholded to 1, log2-transformed, normalized, and baselined to the median of all samples by using GeneSpring GX (Agilent Technologies).

2.2. Gene List Filtering and Differential Expression Analysis

The whole-genome expression data were filtered to include probes targeting genes involved in the splicing core machinery or the auxiliary regulatory factors. In particular, an a priori splicing gene list was generated using query keywords in the Gene Ontology (GO) database (http://www.geneontology.org/). Overall, 396 genes were selected from the following GOs: Spliceosomal complex (GO:0005681), Regulation of mRNA splicing, via spliceosome (GO:0048024), mRNA splicing, via spliceosome (GO: 0000398). Moreover, a number of additional genes (n = 15) were filtered from SPLICE AID-F (http://srv00.recas.ba.infn.it/SpliceAidF/) [21] and included in the gene panel. A full list of the a priori splicing selected genes is reported in Supplementary Data S1.
The statistical analysis of the filtered microarray data was performed with the GeneSpring GX v13.1 software package (Agilent Technologies), applying a two-way analysis of variance (ANOVA) followed by Benjamini–Hochberg false discovery rate (FDR) correction to identify the differentially expressed genes (DEGs) and minimize false-positive cases. A two-way ANOVA test was used, since it offers the possibility to identify genes (i) significantly deregulated across the first variable (disease state: Control or SALS subtypes), (ii) significantly deregulated across the second variable (Tissue types: motor cortex or spinal cord), and (iii) showing a significant interaction between the two variables. The list of probes that were statistically significant, together with their expression levels in the motor cortex and spinal cord of the control and SALS patients, is reported in Supplementary Data S2.
Probes with an adjusted p-value < 0.05 and an FC value greater than ±2 were used as candidate deregulated genes for further exploration in 7 pairwise correlations relevant to the study’s aim (Cortex SALS1 vs. Cortex CTRL; Cortex SALS2 vs. Cortex CTRL; Spinal Cord SALS1 vs. Spinal Cord CTRL; Spinal Cord SALS2 vs. Spinal Cord CTRL; Cortex SALS1 vs. Spinal Cord SALS1; Cortex SALS2 vs. Spinal Cord SALS2; Cortex CTRL vs. Spinal Cord CTRL). For further analysis, the probes were converted into the corresponding genes. When two or more probes targeted the same gene, the mean expression level was considered to address the up- or down-regulation.
To visualize the results, Venn diagrams were generated by using the Interactivenn web tool (http://www.interactivenn.net/) [22].

2.3. Functional Characterization of Differentially Expressed Splicing Factors

To understand which splicing complexes were implicated in ALS, the differential expressed genes were enriched by their cellular component gene ontology (GO) and analyzed statistically by applying Fisher’s Exact test with the use of an overrepresentation test in the web-based, open-access PANTHER database, using the microarray gene list as a reference (see Supplementary Data S3). The p-value were adjusted by using Bonferroni correction for multiple testing and p < 0.05 was set as the significance threshold.

2.4. Protein–Protein Interaction Network Analysis

To evaluate the potential protein–protein interaction (PPI) in the motor cortex, we generated two human tissue-specific networks using NetworkAnalyst v.3.0 [23] and the data belonging to the human tissue-specific db DifferentialNet [24], using as input data two DEG lists in SALS1 vs. control cortex and SALS2 vs. control cortex. The former contained 79 DEGs, while the latter was composed of 176 DEGs. The two gene lists were used as the input data in NetworkAnalyst v.3.0, together with their logFC values, in order to indicate the direction of the expression change, and duplicates were filtered using the log FC mean. Using the NetworkAnalyst tool, the input data were matched with their corresponding human brain cortex-specific proteins using the differential protein–protein interactions human tissue-specific database, i.e., DifferentialNet db, filtered by a percentile threshold of 15%, and were then used to generate tissue-specific first-order PPI networks.
The same procedure was applied to the spinal-cord SALS1 and SALS2 DEG lists, made up, respectively, by 8 and 5 DEGs. Their data were used to create spinal-cord specific first order PPI networks filtered by a percentile threshold of 15%.
The data referring to the descriptive statistics of the tissue-specific networks for both SALS clusters were obtained using the PHStat2 add-in for excel and are shown in the Supplementary Data S5.
To calculate the cut-off for hub-bottlenecks and nonhub-bottlenecks, we referred to the method described in [25], identifying the nodes with a degree greater than or equal to the sum of the mean and twice the standard deviation (mean + 2*S.D. of the degree distribution) as hub-bottlenecks. While nonhub-bottleneck screening was performed by sorting the node betweenness distribution in descending order and selecting the nodes with values lying in the top 50% of the betweenness distribution. The nodes ranging from those with a degree less than the degree cut-off for hub-bottlenecks, directly connected to at least two hub-bottlenecks, and with a betweenness value, were selected as nonhub-bottlenecks.
Venn diagrams for both the hub-bottleneck and nonhub-bottleneck categories in the tissue-specific networks for both ALS clusters were obtained using the Interactivenn web tool [22] (http://www.interactivenn.net/).
A literature-based network was obtained using the Genomatix Pathways System (GePS, www.genomatix.de) using as the input list the tissue-specific DEGs involved in cytoplasmic stress granule formation and belonging to GO:0010494 (CASC3, FMR1, DDX1, PABPC1, YBX1, HNRNPK, MBNL1, CIRBP, TARDBP, and RBPMS) in Homo sapiens. The gene expression log FC values were added using the Metadata import option.

3. Results and Discussion

3.1. Identification of Differentially Expressed Genes in Tissue Types and Disease States

To identify the complete repertoire of the splicing components differentially regulated in SALS, we investigated the expression profiles of 396 splicing factors in the motor cortex (n = 41) and spinal cord (n = 40) samples from control and SALS patients. In two previous studies, we were able to successfully discriminate controls from SALS patients and distinguish these latter in two greatly divergent SALS subtypes (SALS1 and SALS2) by using unsupervised hierarchical clustering [20,26]. Each of these subtypes were characterized by specific transcriptomic alterations associated with divergent deregulated genes and pathways, allowing a molecular taxonomy of patients. The existence of distinct SALS molecular clusters has also been recently confirmed in an RNA-seq study [27]. Since the union of SALS subgroups into a single disease entity would conceal the great existing variability in ALS, in the present study, patients and tissues (both cortex and spinal cord) were analyzed taking into account the previously demonstrated separation. Therefore, we performed the differential gene expression analysis by comparing the controls to each of the two SALS subtypes previously characterized. Figure 1 shows globally the analysis methods used in the present study to characterize the differentially expressed genes encoding splicing factors and investigate their functional significance.
A two-way analysis of variance (ANOVA) with multiple test correction (Benjamini–Hochberg’s FDR test) was applied to measure the gene expression changes influenced by each of the two variables, tissue type (motor cortex and spinal cord) and disease state (control, SALS1, and SALS2), independently or simultaneously. Through this analysis, 259 genes out of 396 were identified as having a significant differential expression (p < 0.05) with respect to the tissue type, disease state or both (Supplementary Data S2). The Venn diagram in Figure 2 displays the number of statistically significant genes influenced independently and/or simultaneously by the two variables. A total of 31 genes were differentially expressed across the tissue types, 231 genes were significantly deregulated across the disease states, whereas 213 genes showed an interaction effect and changed in the two tissue types along with the disease states (Figure 2). Twenty-two genes were deregulated in the disease state, tissue type and also in their interaction effect (Figure 2).
Figure 3 shows the supervised cluster of gene entities (n = 302 probes corresponding to 213 genes) deregulated in the two tissue types, along with the disease state (interaction effect, Figure 2). These transcriptomic profiles confirm the great divergence between SALS1 and SALS2 motor cortex samples, as previously reported in [20], and a less pronounced difference in spinal cord samples.
Next, we focused on genes (n = 184) out of 213 showing major expression changes (fold change > |2|) in different pairwise correlations that were relevant to the study’s aim (Figure 1 and Supplementary Data S3).
The results reveal that major changes according to the disease state in the different tissue types occurred, particularly in the SALS2 motor cortex vs. the control cortex (DEGs = 139, the majority of which were down-regulated), followed by the SALS1 cortex vs. the control (DEGs = 63, the majority were up-regulated), while minor expression change events occurred in the spinal cord of both SALS1 and SALS2 compared to the controls (DEGs = 8 and 5, respectively) (Figure 4 and Supplementary Data S3). YBX1 was the only gene commonly deregulated in the comparison of tissue types; TRA2A was down-regulated in both SALS1 and SALS2 spinal cords vs. the controls and up-regulated in the SALS2 cortex compared to the controls. SRRM2 and RNF113B were differentially expressed in both SALS1 and SALS2 cortexes, and SALS1 spinal cords (Figure 4a,b).
The major changes according to tissue type in the different disease states are shown in Figure 4c,d and Supplementary Data S3. A small number of genes (n = 14) were differentially expressed in the motor cortex vs. spinal cord of control patients, whereas a significant deregulation was observed in SALS1 (DEGs = 102, mostly up-regulated in the cortex) and SALS2 (DEGs = 137, mostly down-regulated in the cortex). Four genes were deregulated in all tissue type comparisons (A2BP1, HNRNPA1, POLR2D, RNF113B).
The differential expression of genes encoding the subset of splicing factors in motor cortex of SALS1 and SALS2 patients was compared with their corresponding whole transcriptome analysis previously published [20]. This comparison showed a complete overlap in terms of genes and their regulation for splicing-related DEGs and, at the same time, highlighted the contribution of the gene set focused analysis in putting in light the deregulation of other splicing-related genes not previously assessed (Supplementary Data S3). Moreover, a further comparison between DEGs and their copy number variants (CNVs), previously screened using a neurological exon-custom aCGH [28], was also performed, highlighting the SMN1 recurrent loss (penetrance > 10%) in SALS2 patients (chr5:70220768-70249769) with a concordance between gene expression deregulation and CNV (Supplementary Data S3).

3.2. Deregulated Splicing Complexes

DEGs encode factors that are involved in different splicing complexes. functional enrichment analysis of their GOs allowed us to group these factors into different splicing complexes and identify which of these were significantly deregulated in specific disease states and tissue types (Table 2 and Supplementary Data S3). Functional enrichment analysis was performed starting from DEGs in each of the seven pairwise comparisons. To reduce redundancy and simplify their comprehension, the most significant overrepresented GOs are summarized in Table 2 and are listed in Supplementary Data S3bis. In the following paragraphs, we will discuss the expression changes of these splicing factors by focusing on the main functional groups (splicing complexes), highlighting common features and specific differences.

3.2.1. Spliceosomal Complex Core

U1 snRNP

Among the five canonical ribonucleoprotein complexes, U1 snRNP is the most studied in ALS. Genes composing U1 snRNP were globally downregulated in SALS2 motor cortexes compared to the control, and with a divergent expression trend in SALS1 and SALS2 tissue type comparisons. U1 is constituted of 18 components and is well known to interact with the ALS-associated FET (FUS, EWSR1, TAF15) proteins, with MATR3, as well as with the SMN complex, required for intranuclear gem formation (Gemini of coiled bodies) and snRNP biogenesis [29,30,31,32]. Mutations in FUS determined a poor interaction with U1, causing a mis-localization to the cytoplasm and a global loss of splicing activity [33,34]. In addition, the knockdown of any of the U1 snRNP-specific proteins resulted in a dramatic loss of SMN-containing gems and in motor axon truncations in zebrafish models [34].

U2 snRNP

Several DEGs in the motor cortexes of both disease states and in SALS tissue type comparisons were overrepresented in U2-coordinated spliceosomal complexes (Table 2). Enriched GOs included the U2-type pre-spliceosome (formed by association of U1 snRNP to the 5′ splice site and U2 snRNP to the branch point sequence), the precatalytic spliceosome (determined by the recruitment of the preassembled U4/U6.U5 tri-snRNP), the catalytic step 1 (or GT-AG activated spliceosome) and the catalytic step 2 spliceosome (an intermediate complex in which the first catalytic cleavage has already occurred).
The disruption of U2 snRNP localization, activity and function was described in ALS patients carrying the expanded C9ORF72. This latter encodes aberrant dipeptide repeat (DPR) proteins, which prevent spliceosome assembly through interaction with U2 particles [35], blocking them in the cytoplasm far away from splicing sites and producing mis-spliced cassette exons in C9ORF72 patient brains [35].
A number of U2 subunit DEGs (SF3B5, SNRNP200, SF3B1, MFAP1, PABPC1) interact with, or are genetic modulators of, ALS pathogenic genes, including FUS, TDP-43, and VCP [30,36,37,38,39], while others (SYF2) have been associated with neuronal apoptosis [40]. Noteworthy, a consistent dysregulation of the SRRM2 splicing factor was previously described in substantia nigra, amygdala and the peripheral blood of parkinsonian patients, drawing attention to the role of this gene in neurodegeneration [41].

U4, U5, U6 snRNPs and tri-snRNPs Complex

U4/U6.U5 tri-snRNP is a 1.5 MDa pre-assembled spliceosomal complex comprising U5 snRNA, extensively base-paired U4/U6 snRNAs and >30 proteins. This tri-snRNP combines with a pre-mRNA substrate bound to U1 and U2 snRNPs and transforms into a catalytically active spliceosome, following extensive compositional and conformational changes triggered by the unwinding of the U4/U6 snRNAs [42].
U6 snRNA expression levels are regulated by TDP-43 and its amount significantly decreased after TDP-43 depletion, leading to neuronal cell death [43]. Two U6 components (LSM7 and LSM4) were involved in FUS toxicity suppression and exhibited alterations in ALS fibroblasts [44,45], while SNRPN and SNRPC expression changed following FUS depletion and toxic C9ORF72-derived DPR exposure [35,46].
Other interesting splicing factors include PRPF6 and PRPF8, which were previously identified as hub genes in a network/pathway-based analysis conducted to investigate differential pathogenic mechanisms in sporadic ALS [47]. Interestingly, PRPF6 is located inside an ALS critical genetic region, carrying SNPs significantly associated with ALS pathogenesis in Japanese people [48].

3.2.2. U12-Dependent Minor Splicing Pathway

Despite the majority of overrepresented GOs being related to the U2-dependent splicing pathway, some DEGs are involved in the minor splicing pathway (U12-dependent). The minor counterpart of the splicing machinery regulates the splicing of introns with atypical AT–AC terminal dinucleotides, and has received increasing attention during the last few years as a novel pathomechanistic player in both neurodevelopmental and neurodegenerative diseases [9].
It has been recently suggested as a possible pathomechanism by which mutated FUS inhibits the correct minor splicing pathway. Mutated FUS seems to forms cytoplasmic aggregates, which are able to trap U11 and U12 snRNAs inside [49]. Moreover, a TDP-43 loss-of-function was shown to determine a disturbance of U11/U12 spliceosomes in both tissues affected by ALS and cellular models [50].
YBX1 (or YB-1) deserves a particular mention, as it is the only gene commonly deregulated in both tissue types and disease states. This gene, a highly conserved transcription factor known to be a potential drug target in cancer therapy [51,52], is likely to have important local implications for motor neuron activity. Its mRNA was recently reported as being abundantly enriched in distal axons of motor neurons compared to the soma compartment, where it may support the local regulation of distal axonal processes [53]. It has also been described as an essential inducer-splicing of the human muscle-specific receptor tyrosine kinase (MuSK) gene, a postsynaptic transmembrane molecule that mediates the clustering of acetylcholine receptors [54].

3.2.3. Exon–Exon Junction Complex (EJC)

The exon junction complex (EJC) is a protein complex deposited onto mRNAs following splicing, and is able to both stably bind mRNAs and function as an anchor for further processing steps. Indeed, EJC acts as a central node of post-transcriptional gene expression networks and regulates mRNAs’ cellular localization, translation and non-sense mediated decay (NMD) [55].
The dysregulation of EJC was found in both SALS1 and SALS2 motor cortexes versus controls, and in SALS tissue type comparisons, but with opposite expression trends (Table 2). Interestingly, one of the central EJC effectors (UPF1) is a key-player in ALS and has been proposed as a novel therapeutic target for treatment. Its expression preserves forelimb function in a rat model of TDP-43-induced motor paralysis [56] and reduces toxicity in primary mammalian neurons expressing mutant FUS or TDP-43 [57,58], as well as efficiently blocking neurotoxicity caused by C9ORF72 DPR [59,60]. The depletion of the EJC member SRSF1 inhibits the nuclear exporting of pathological C9ORF72 transcripts, the production of DPR and alleviates neurotoxicity in several human neuronal models [61,62]. Other EJC constituents (UPF3B and EIF4A3) were upregulated in ALS cell lines and were shown to interact with FUS and TDP-43 [58,63].

3.2.4. CRD-Mediated mRNA Stability Complex

The CRD-mediated mRNA stability complex is a protein complex that binds to and promotes the stabilization of mRNA molecules containing the coding region instability determinant (CRD). Among the complex members, HNRNPU was previously identified as one of five top-ranked RNA binding proteins, showing significant alterations in the motor neurons, cerebellum and spinal cords of ALS patients compared to the control, and co-localizes with cytoplasmic TDP-43 aggregates [64,65]. Moreover, we observed in SALS2 motor cortexes vs. control a down-regulation of SYNCRIP, whose overexpression was previously found to ameliorate spinal muscular atrophy pathological phenotypes [66].

3.2.5. RNA Polymerase II, Core Complex

In SALS2 motor cortexes and both SALS tissue type comparisons, we observed the deregulation of genes involved in the RNA polymerase II core complex. Thanks to the association of U1 snRNP with the RNAP II machinery, this transcription potently enhances splicing, and vice versa, by a mechanism of reciprocal coupling [67,68].
RNAP2 is the largest subunit of RNA polymerase II, and its phosphorylation state at the C-terminal domain is orchestrated by FUS. When FUS is missing, RNAP2 accumulates at the transcription start site and causes a shift in mRNA transcripts expression toward early polyadenylation sites [69]. Moreover, defects in FUS and TDP-43 recruitment influence RNAPII termination and determine R-loop (DNA:RNA hybrids) accumulation, leading to elevated DNA damage at transcription terminators [70]. Also, the recently described ALS-associated gene MATR3 interacts with POLR2A, POLR2B, and POLR2C at the early steps of transcription [71].

3.2.6. Prp19 Complex

The conserved Prp19 complex (Prp19C)—also known as NineTeen Complex (NTC)—acts in several processes for cellular homeostasis maintenance, but its best-characterized function is in splicing. NTC/Prp19C is a non-snRNP splicing complex crucial that plays multiple roles: it associates with the assembling spliceosome during or after the dissociation of the U4 snRNP, stabilizes the U5/U6 snRNP in the activated spliceosomal complex (B-act) and remains linked with the spliceosome during the second step of splicing [72]. The DEGs involved in the Prp19 complex were specifically deregulated in SALS2 cortexes vs. control cortexes and in SALS2 tissue type comparison (Table 2).
A Prp19 complex member (HSPA8) encodes a constitutive heat shock protein which is able to associate with TDP-43 [36,73]. The level of this protein was found reduced at the neuromuscular junctions in both mutated TDP-43 mice and human models, as well as in motor neurons carrying the C9ORF72 repeat expansion [74].

3.2.7. SMN–Sm Protein Complex

The SMN–Sm protein complex is a stable multiprotein complex present in all eukaryotic cells that surveys the correct identity of the target RNAs and mediates the biogenesis of snRNPs. It is composed of SMN (the survival motor neuron) protein, at least six additional proteins named Gemins2–7, which associate with both Sm proteins and snRNAs [75,76].
SMN1, the core of the complex, has received great attention in the study of ALS, especially because it is located in a genomic region that frequently undergoes copy number abnormalities [77]. A number of studies have demonstrated that high or low SMN1 copy numbers increased the risk of developing ALS, although data in the literature are conflicting (all reviewed in [77]). Decreased levels of SMN protein cause a reduction of the chaperone-like activity following SOD1-mediated toxicity, and thus contribute to free radical injury and oxidative stress [78]. On the contrary, increased SMN protein levels protect motor neurons from stress- or ALS mutation-induced death, raising important clinical implications for SMN therapeutics in motor neuronal diseases [79].

3.2.8. Supraspliceosomal Complex

The supraspliceosome is a huge (21 MDa) and highly dynamic RNP machine composed of four active spliceosomes connected to each other by the pre-mRNA [80,81]. This multi-task regulator harbors all five spliceosomal U snRNPs at all splicing stages, and additional pre-mRNA processing components, such as the 5′-end and 3′-end processing factors, and the RNA editing enzymes ADAR1 and ADAR2 [80]. ADARs, in particular, catalyzes the hydrolytic deamination of adenosine to inosine in double-stranded RNA (dsRNA), and this process is referred to as A-to-I RNA editing. ADAR isoform down-regulation was previously observed in neurons carrying C9ORF72 mutations and in TDP43-deficient cell lines; moreover these proteins are incorporated into RNA-rich nuclear foci and are targeted directly by TDP43 [82].

3.3. Tissue-Specific Networks of Spliceosome-Related Gene Products

To investigate the functional relevance of splicing-related proteins and their direct interactors in both tissues, we generated first-order PPI networks using DEGs which were known to be specific for those tissues, and expanding them to their direct interactions, which were also tissue-specific. After that, we investigated their hierarchical role through the identification of the topological significance of their nodes, as described in the following.
Network topological properties have been described in the literature in terms of functional significance, making a difference between a) highly connected nodes (hubs), and b) nodes which are highly passed through as the shortest paths (bottlenecks). While the former are characterized by a high degree of connection, the latter are identified by high betweenness values, and, in turn, can be divided into nonhub-bottlenecks and hub-bottlenecks, depending on their degree of connection. Since in particular cases of PPI networks with indirect edges, such as signaling transduction and permanent interactions (i.e., protein complexes), betweenness tend to be more essential than degree [83], we analyzed first-order tissue-specific networks of spliceosome-related gene products, looking for bottlenecks categories with high and low degrees of connection. In fact, both categories have their importance, not only in highlighting network centrality (hub-bottlenecks), which could disrupt the network if removed, but also in finding peripheral proteins mediating different pathways, connecting different complexes and/or being involved in cross-talking (nonhub-bottlenecks), both having key roles.
In the next paragraph, we will illustrate the PPI networks of splicing factors deregulated in SALS subtypes, either in the cortex or in the spinal cord.

3.3.1. Cortex-Specific PPI Networks

The PPI networks of splicing factors deregulated in SALS1 and SALS2 cortexes and their direct (first-order) interaction partners are illustrated in Figure 5 (network properties and topological significance analysis are described in Supplementary Data S5). Using 55 deregulated tissue-specific genes as seeds for SALS1, we obtained 19 hub- and 338 nonhub-bottlenecks, while 35 hub- and 441 nonhub-bottlenecks were identified for SALS2, starting with 122 seeds (see Supplementary Data S5 for details about hub- and nonhub-bottlenecks). The profound divergence of SALS1 and SALS2 subtypes is evident in Figure 5. Among hub-bottlenecks (depicted in Figure 5 as spheres higher in size), 10 are common to the two SALS subtypes and most of them (TARDBP, HNRNPA1, HNRNPC, HNRNPD, HNRNPK, SNWI, YBX1, HSPA8, DHX9 and PCBP1) shows an opposite expression trend in the cortex networks, whereas only SSRM2 is overexpressed in both.

3.3.2. Spinal Cord-Specific PPI Networks

The PPI networks of splicing factors that are deregulated in SALS1 and SALS2 spinal cords and their direct (first-order) interaction partners are illustrated in Figure 6 (network properties and topological significance analysis are described in Supplementary Data S5). According to DEGs data, spinal cord-specific first-order PPI networks (Figure 6) show a less characterized distinction between SALS1 and SALS2 when compared to those obtained in the cortex (Figure 5). Using six deregulated tissue-specific genes as seeds for SALS1, we obtained five hub- and 92 nonhub-bottlenecks, while only one hub-bottleneck was identified for SALS2 starting with three seeds (see Supplementary Data S5 for details about hub- and nonhub-bottlenecks). As shown in Figure 6, only one hub-bottleneck is common to both disease subtypes, YBX1, with a down-regulation of its transcript in both, whereas SSRM2 appears to be a SALS1-specific hub-bottleneck with an up-regulation of its transcript.

3.3.3. The Most Relevant Splicing-Related Genes

Among the 22 most relevant splicing-related genes (deregulated independently and simultaneously by the two variables, disease state, and tissue type), ten had a relevant topological role (hub- or non-hub-bottleneck) in at least one of the four disease state pairwise comparisons, either in the motor cortex or the spinal cord (see Supplementary Data S5): YBX1, HNRNPU, HNRNPC, RBM3, PABPC1, DNAJC8, HNRNPA3, HNRNPUL1, SNRPN, and YTHDC1. Interestingly, three of these (HNRNPUL1, SNRPN, and YTHDC1) have been recently reported as differentially expressed in peripheral blood sample transcription profiling in a huge heterogeneous ALS cohort (not only sporadic cases), but their up or down regulation was not specified [84]. Moreover, recent evidence has shown that YBX1, HNRNPC, PABPC1, HNRNPA3, and YTHDC1 RNAs bound to TDP-43 protein in human SH-SY5Y neuroblastoma cells using eCLIP-seq [27].
YBX1 deserves a particular mention, as not only is it differentially expressed in all four disease state comparisons, but is also known to participate in several steps of the spliceosomal machinery, playing a central topological role as a hub-bottleneck. Recent evidences highlighted the role of YBX1 in regulating protein recruitment for the formation of stress granules [85], which represents a pathological marker associated with multiple neurodegenerative disorders [86]. The ALS-associated proteins TDP-43 and FUS are known to aggregate in stress granules and can be released by YBX1, which represents a key source of protein aggregation regulation for assembling/disassembling stress granules in motor neurons [87]. Based on this evidence, we generated a regulatory network of DEGs implicated in cytoplasmic stress granule formation (GO: 0010494: CASC3, FMR1, DDX1, PABPC1, YBX1, HNRNPK, MBNL1, CIRBP, and TARDBP) to investigate the hierarchical flow coordinating stress granule formation. As shown in Figure 7, the hub-bottleneck YBX1 is the top-gene regulating all the tissue-specific DEGs involved in stress granule formation. The majority of these have a topological significance and are deregulated in at least one condition. Interestingly, the transcription rate of five of them (HNRNPK, FMR1, CIRBP, MBNL1, and PABPC1) may be regulated by YBX1, since they own in their promoters the corresponding predicted transcription factor binding site. Given the role of YBX1 in a number of spliceosomal complexes, its deregulation at different genomic levels may not only alter splicing but also influence gene expression levels of proteins recruited in stress granules, causing aberrations in the assembly/disassembly metabolism.

4. Conclusions

A number of scientific studies support the role of alternative splicing dysregulation in ALS. Fluctuations in the concentration of core spliceosomal proteins, as well as additional regulatory factors, may alter gene expression and underly the selective neurodegeneration in ALS.
To date, splicing defects in ALS have been investigated at the singular gene level, without considering the higher-order level involving the splicing machinery. In this study, we analyzed the entire spectrum of genes encoding splicing factors in the motor cortex and spinal cord samples from control and SALS patients, previously characterized as two distinct molecular subgroups (SALS1 and SALS2) [20,26]. The data presented show that a substantial number of splicing-related genes displayed significant expression changes according to tissue type and disease state. Overexpression was a trend for the SALS1 group, while down-regulation was a trend for SALS2, although the DEGs were more pronounced for the motor cortex regions than the spinal cord. A number of DEGs in both SALS1 and SALS2 motor cortexes were involved in spliceosomal complexes or step reaction regulation (i.e., U2-type precatalytic and catalytic step reactions, the U4/U6.5 tri-snRNP complex, the U12-type minor spliceosomal pathway, the CRD-mediated mRNA stability complex and the exon–exon junction complex) while other overrepresented GOs were specific to SALS2 (i.e., supraspliceosomal complex, SMN–Sm protein complex, RNA polymerase II core complex). Data integration between gene expression profiling and topological significance (hub-bottlenecks and nonhub bottlenecks) in network analysis aided us to stratify genes based on their deregulation impact, identifying those that are likely to be more relevant in the balance between physiological and pathological status. Deregulation of one of these bottlenecks, including YBX1, might impact stress granule formation, a pathological marker associated with ALS.
Despite it not being possible to completely distinguish the causative factors from secondary degenerative changes ongoing in the diseased brain, it is interesting to note that most of the deregulated splicing factors are known to interact with the ALS-causative mendelian genes, thus suggesting that perturbations in a small subset of splicing components could trigger pathological mechanisms and influence motoneuronal degeneration in sporadic cases.
The existence of transcriptional differences separating SALS patients in two (or more) distinct sub-clusters might explain why several compounds showing promising results in preclinical studies failed to translate into successful clinical trials [27,88]. Indeed, the lack of therapeutic progress may mainly due to an insufficient understanding of the complexity and heterogeneity underlying ALS [89,90]. The investigation of splicing deregulation in both the motor cortex and spinal cord samples into molecularly separated disease groups revealed significant alterations that would be hidden considering ALS as a single disease entity. Moreover, elucidating the impact of splicing machinery changes on global RNAs deregulation may provide novel mechanistic insights underlying the selective degeneration of motor neurons.
Future functional and clinical investigation will certainly be necessary to assess the potential role of candidates’ splicing factors in affecting the origin or modifying the progression of the disease. Notably, the recent technological advancements have made possible the development and use of RNA-targeting molecules as potential therapeutic strategies to increase the expression or mediate the reduction of a given RNA [1,91,92]. This kind of approach may be used in the future to specifically target upstream regulators (e.g., splicing factors or RNA-binding proteins), avoiding large RNAs downstream deregulation that could create more extended damage to cellular physiology.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4409/9/1/159/s1. Data S1: List of splicing-related genes. Data S2: List of significant differentially expressed probes, as resulted from Two Way-ANOVA test and post-hoc correction. Data S3:. List of genes included in the Agilent 4x44K array and used as reference for GO enrichment (sheet 1); up and down deregulated genes in each of the seven pair comparisons (sheet 2); Comparisons between splicing related DEGs, whole transcriptome DEGs and CNVs (sheets 3 and 4). Data S3bis: Over-represented GOs cellular component enrichment of DEGs listed in S3 for each of the seven pairwise comparisons. Data S4: Brain cortex-specific PPI network nodes for both SALS1 and SALS2 clusters (sheet 1 and sheet 2) and spinal cord-specific PPI network nodes (sheet 3 and sheet 4). Data S5: Statistics of degree distribution underlying topological analysis for both brain cortex-specific and spinal cord-specific networks for SALS1 and SALS2 disease clusters.

Author Contributions

V.L.C. was responsible for data curation, formal analysis, data interpretation, and wrote the manuscript, G.G. was responsible for PPI and TF network analysis, and participated in data interpretation. E.A. conducted tissue sample preparation, S.C. designed the concept of the study and supervised it. G.G. and S.C. critically revised the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by the grant (CIP 2014.IT.05.SFOP.014/3/10.4/9.2.10/0008) from the European Social Fund operational program for the Sicily region (Italy) “Development and application of biosensoristic technologies in genomics” PO FESR Sicily 2014-2020. EA is supported by the Stichting ALS Nederland, “The Dutch ALS Tissue Bank”.

Acknowledgments

Authors gratefully acknowledge Cristina Calì, Alfia Corsino, Maria Patrizia D’Angelo and Francesco Marino for administrative and technical support. We acknowledge the team who helped in the collection of ALS tissue samples (D. Troost, M. de Visser, A.J. van der Kooi and J. Raaphorst) and Lucia Mesarosova and Jasper Anink for technical support. This work is dedicated to the memory of our collegue Barbara Maino.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Butti, Z.; Patten, S.A. RNA Dysregulation in Amyotrophic Lateral Sclerosis. Front. Genet. 2018, 9, 712. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Perrone, B.; La Cognata, V.; Sprovieri, T.; Ungaro, C.; Conforti, F.L.; Andò, S.; Cavallaro, S. Alternative Splicing of ALS Genes: Misregulation and Potential Therapies. Cell. Mol. Neurobiol. 2019, 40, 1–14. [Google Scholar] [CrossRef]
  3. Su, C.-H.; Tarn, W.-Y. Alternative Splicing in Neurogenesis and Brain Development. Front. Mol. Biosci. 2018, 5, 12. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Raj, B.; Blencowe, B.J. Alternative Splicing in the Mammalian Nervous System: Recent Insights into Mechanisms and Functional Roles. Neuron 2015, 87, 14–27. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Weyn-Vanhentenryck, S.M.; Feng, H.; Ustianenko, D.; Duffié, R.; Yan, Q.; Jacko, M.; Martinez, J.C.; Goodwin, M.; Zhang, X.; Hengst, U.; et al. Precise temporal regulation of alternative splicing during neural development. Nat. Commun. 2018, 9, 2189. [Google Scholar] [CrossRef] [Green Version]
  6. Porter, R.S.; Jaamour, F.; Iwase, S. Neuron-specific alternative splicing of transcriptional machineries: Implications for neurodevelopmental disorders. Mol. Cell. Neurosci. 2018, 87, 35–45. [Google Scholar] [CrossRef]
  7. Hoskins, A.A.; Moore, M.J. The spliceosome: A flexible, reversible macromolecular machine. Trends Biochem. Sci. 2012, 37, 179–188. [Google Scholar] [CrossRef] [Green Version]
  8. Dvinge, H. Regulation of alternative mRNA splicing: Old players and new perspectives. FEBS Lett. 2018, 592, 2987–3006. [Google Scholar] [CrossRef]
  9. Jutzi, D.; Akinyi, M.V.; Mechtersheimer, J.; Frilander, M.J.; Ruepp, M.-D. The emerging role of minor intron splicing in neurological disorders. Cell Stress 2018, 2, 40–54. [Google Scholar] [CrossRef]
  10. Verma, B.; Akinyi, M.V.; Norppa, A.J.; Frilander, M.J. Minor spliceosome and disease. Semin. Cell Dev. Boil. 2018, 79, 103–112. [Google Scholar] [CrossRef]
  11. La Cognata, V.; Iemmolo, R.; D’Agata, V.; Scuderi, S.; Drago, F.; Zappia, M.; Cavallaro, S. Increasing the Coding Potential of Genomes Through Alternative Splicing: The Case of PARK2 Gene. Curr. Genom. 2014, 15, 203–216. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Turunen, J.J.; Niemelä, E.H.; Verma, B.; Frilander, M.J. The significant other: Splicing by the minor spliceosome. Wiley Interdiscip. Rev. RNA 2013, 4, 61–76. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Patel, S.B.; Bellini, M. The assembly of a spliceosomal small nuclear ribonucleoprotein particle. Nucleic Acids Res. 2008, 36, 6482–6493. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. La Cognata, V.; D’Agata, V.; Cavalcanti, F.; Cavallaro, S. Splicing: Is there an alternative contribution to Parkinson’s disease? Neurogenetics 2015, 16, 245–263. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Dvinge, H.; Kim, E.; Abdel-Wahab, O.; Bradley, R.K. RNA splicing factors as oncoproteins and tumour suppressors. Nat. Rev. Cancer 2016, 16, 413–430. [Google Scholar] [CrossRef]
  16. Zhou, Z.; Fu, X.-D. Regulation of splicing by SR proteins and SR protein-specific kinases. Chromosoma 2013, 122, 191–207. [Google Scholar] [CrossRef]
  17. Lee, Y.; Rio, D.C. Mechanisms and Regulation of Alternative Pre-mRNA Splicing. Annu. Rev. Biochem. 2015, 84, 291–323. [Google Scholar] [CrossRef] [Green Version]
  18. Ule, J.; Blencowe, B.J. Alternative Splicing Regulatory Networks: Functions, Mechanisms, and Evolution. Mol. Cell 2019, 76, 329–345. [Google Scholar] [CrossRef]
  19. Wang, Y.; Liu, J.; Huang, B.O.; Xu, Y.M.; Li, J.; Huang, L.F.; Lin, J.; Zhang, J.; Min, Q.H.; Yang, W.M.; et al. Mechanism of alternative splicing and its regulation. Biomed. Rep. 2015, 3, 152–158. [Google Scholar] [CrossRef] [Green Version]
  20. Aronica, E.; Baas, F.; Iyer, A.; Asbroek, A.L.T.; Morello, G.; Cavallaro, S. Molecular classification of amyotrophic lateral sclerosis by unsupervised clustering of gene expression in motor cortex. Neurobiol. Dis. 2015, 74, 359–376. [Google Scholar] [CrossRef] [Green Version]
  21. Giulietti, M.; Piva, F.; D’Antonio, M.; D’Onorio De Meo, P.; Paoletti, D.; Castrignano, T.; D’Erchia, A.M.; Picardi, E.; Zambelli, F.; Principato, G.; et al. SpliceAid-F: A database of human splicing factors and their RNA-binding sites. Nucleic. Acids Res. 2013, 41, D125–D131. [Google Scholar] [CrossRef] [PubMed]
  22. Heberle, H.; Meirelles, G.V.; Da Silva, F.R.; Telles, G.P.; Minghim, R. InteractiVenn: A web-based tool for the analysis of sets through Venn diagrams. BMC Bioinform. 2015, 16, 169. [Google Scholar] [CrossRef]
  23. Zhou, G.; Soufan, O.; Ewald, J.; Hancock, R.E.W.; Basu, N.; Xia, J. NetworkAnalyst 3.0: A visual analytics platform for comprehensive gene expression profiling and meta-analysis. Nucleic Acids Res. 2019, 47, W234–W241. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Basha, O.; Shpringer, R.; Argov, C.M.; Yeger-Lotem, E. The DifferentialNet database of differential protein-protein interactions in human tissues. Nucleic Acids Res. 2018, 46, D522–D526. [Google Scholar] [CrossRef] [PubMed]
  25. Rakshit, H.; Rathi, N.; Roy, D. Construction and Analysis of the Protein-Protein Interaction Networks Based on Gene Expression Profiles of Parkinson’s Disease. PLoS ONE 2014, 9, e103047. [Google Scholar] [CrossRef] [PubMed]
  26. Morello, G.; Spampinato, A.G.; Cavallaro, S. Molecular Taxonomy of Sporadic Amyotrophic Lateral Sclerosis Using Disease-Associated Genes. Front. Neurol. 2017, 8, 3. [Google Scholar] [CrossRef] [Green Version]
  27. Tam, O.H.; Rozhkov, N.V.; Shaw, R.; Kim, D.; Hubbard, I.; Fennessey, S.; Propp, N.; Fagegaltier, D.; Ostrow, L.; Phatnani, H.; et al. Postmortem Cortex Samples Identify Distinct Molecular Subtypes of ALS: Retrotransposon Activation, Oxidative Stress, and Activated Glia. Cell Rep. 2019, 29, 1164–1177.e5. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Morello, G.; Guarnaccia, M.; Spampinato, A.G.; Salomone, S.; D’Agata, V.; Conforti, F.L.; Aronica, E.; Cavallaro, S. Integrative multi-omic analysis identifies new drivers and pathways in molecularly distinct subtypes of ALS. Sci. Rep. 2019, 9, 9968. [Google Scholar] [CrossRef]
  29. Yamazaki, T.; Chen, S.; Yu, Y.; Yan, B.; Haertlein, T.C.; Carrasco, M.A.; Tapia, J.C.; Zhai, B.; Das, R.; Lalancette-Hebert, M.; et al. FUS-SMN protein interactions link the motor neuron diseases ALS and SMA. Cell Rep. 2012, 2, 799–806. [Google Scholar] [CrossRef] [Green Version]
  30. Groen, E.J.; Fumoto, K.; Blokhuis, A.M.; Engelen-Lee, J.; Zhou, Y.; Heuvel, D.M.V.D.; Koppers, M.; Van Diggelen, F.; Van Heest, J.; Demmers, J.A.; et al. ALS-associated mutations in FUS disrupt the axonal distribution and function of SMN. Hum. Mol. Genet. 2013, 22, 3690–3704. [Google Scholar] [CrossRef] [Green Version]
  31. Therrien, M.; Rouleau, G.A.; Dion, P.A.; Parker, J.A. FET proteins regulate lifespan and neuronal integrity. Sci. Rep. 2016, 6, 25159. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Chi, B.; O’Connell, J.D.; Yamazaki, T.; Gangopadhyay, J.; Gygi, S.P.; Reed, R. Interactome analyses revealed that the U1 snRNP machinery overlaps extensively with the RNAP II machinery and contains multiple ALS/SMA-causative proteins. Sci. Rep. 2018, 8, 8755. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Sun, S.; Ling, S.-C.; Qiu, J.; Albuquerque, C.P.; Zhou, Y.; Tokunaga, S.; Li, H.; Qiu, H.; Bui, A.; Yeo, G.W.; et al. ALS-causative mutations in FUS/TLS confer gain- and loss-of-function by altered association with SMN and U1-snRNP. Nat. Commun. 2015, 6, 6171. [Google Scholar] [CrossRef] [PubMed]
  34. Yu, Y.; Chi, B.; Xia, W.; Gangopadhyay, J.; Yamazaki, T.; Winkelbauer-Hurt, M.E.; Yin, S.; Eliasse, Y.; Adams, E.; Shaw, C.E.; et al. U1 snRNP is mislocalized in ALS patient fibroblasts bearing NLS mutations in FUS and is required for motor neuron outgrowth in zebrafish. Nucleic Acids Res. 2015, 43, 3208–3218. [Google Scholar] [CrossRef]
  35. Yin, S.; Lopez-Gonzalez, R.; Kunz, R.C.; Gangopadhyay, J.; Borufka, C.; Gygi, S.P.; Gao, F.-B.; Reed, R. Evidence that C9ORF72 Dipeptide Repeat Proteins Associate with U2 snRNP to Cause Mis-splicing in ALS/FTD Patients. Cell Rep. 2017, 19, 2244–2256. [Google Scholar] [CrossRef] [Green Version]
  36. Freibaum, B.D.; Chitta, R.K.; High, A.A.; Taylor, J.P. Global Analysis of TDP-43 Interacting Proteins Reveals Strong Association with RNA Splicing and Translation Machinery. J. Proteome Res. 2010, 9, 1104–1120. [Google Scholar] [CrossRef] [Green Version]
  37. Pons, M.; Miguel, L.; Miel, C.; Avequin, T.; Juge, F.; Frebourg, T.; Campion, D.; Lecourtois, M. Splicing factors act as genetic modulators of TDP-43 production in a new autoregulatory TDP-43 Drosophila model. Hum. Mol. Genet. 2017, 26, 3396–3408. [Google Scholar] [CrossRef]
  38. Rode, S.; Ohm, H.; Zipfel, J.; Rumpf, S. The spliceosome-associated protein Mfap1 binds to VCP in Drosophila. PLoS ONE 2017, 12, e0183733. [Google Scholar] [CrossRef] [Green Version]
  39. Kim, H.J.; Raphael, A.R.; LaDow, E.S.; McGurk, L.; Weber, R.A.; Trojanowski, J.Q.; Lee, V.M.; Finkbeiner, S.; Gitler, A.D.; Bonini, N.M. Therapeutic modulation of eIF2alpha phosphorylation rescues TDP-43 toxicity in amyotrophic lateral sclerosis disease models. Nat. Genet. 2014, 46, 152–160. [Google Scholar] [CrossRef] [Green Version]
  40. Xu, W.; Cao, M.; Zheng, H.; Tan, X.; Li, L.; Cui, G.; Xu, J.; Cao, J.; Ke, K.; Wu, Q. Upregulation of SYF2 is associated with neuronal apoptosis caused by reactive astrogliosis to neuroinflammation. J. Neurosci. Res. 2014, 92, 318–328. [Google Scholar] [CrossRef]
  41. Shehadeh, L.A.; Yu, K.; Wang, L.; Guevara, A.; Singer, C.; Vance, J.; Papapetropoulos, S. SRRM2, a Potential Blood Biomarker Revealing High Alternative Splicing in Parkinson’s Disease. PLoS ONE 2010, 5, e9104. [Google Scholar] [CrossRef] [PubMed]
  42. Nguyen, T.H.D.; Galej, W.P.; Bai, X.-C.; Savva, C.G.; Newman, A.J.; Scheres, S.H.W.; Nagai, K. The architecture of the spliceosomal U4/U6.U5 tri-snRNP. Nature 2015, 523, 47–52. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Yahara, M.; Kitamura, A.; Kinjo, M. U6 snRNA expression prevents toxicity in TDP-43-knockdown cells. PLoS ONE 2017, 12, e0187813. [Google Scholar] [CrossRef] [PubMed]
  44. Sun, Z.; Diaz, Z.; Fang, X.; Hart, M.P.; Chesi, A.; Shorter, J.; Gitler, A.D. Molecular determinants and genetic modifiers of aggregation and toxicity for the ALS disease protein FUS/TLS. PLoS Boil. 2011, 9, e1000614. [Google Scholar] [CrossRef] [PubMed]
  45. Raman, R.; Allen, S.P.; Goodall, E.F.; Kramer, S.; Ponger, L.-L.; Heath, P.R.; Milo, M.; Hollinger, H.C.; Walsh, T.; Highley, J.R.; et al. Gene expression signatures in motor neurone disease fibroblasts reveal dysregulation of metabolism, hypoxia-response and RNA processing functions. Neuropathol. Appl. Neurobiol. 2015, 41, 201–226. [Google Scholar] [CrossRef]
  46. Lagier-Tourenne, C.; Polymenidou, M.; Hutt, K.R.; Vu, A.Q.; Baughn, M.; Huelga, S.C.; Clutario, K.M.; Ling, S.-C.; Liang, T.Y.; Mazur, C.; et al. Divergent roles of ALS-linked proteins FUS/TLS and TDP-43 intersect in processing long pre-mRNAs. Nat. Neurosci. 2012, 15, 1488–1497. [Google Scholar] [CrossRef]
  47. Han, L.-H.; Fan, X.-Y.; Guo, H.; Wei, W.; Chen, M.-M.; Yan, S.-F. Network strategy to investigate differential pathways in sporadic amyotrophic lateral sclerosis. J. Cancer Res. Ther. 2018, 14, 1057. [Google Scholar] [CrossRef]
  48. Iida, A.; Takahashi, A.; Kubo, M.; Saito, S.; Hosono, N.; Ohnishi, Y.; Kiyotani, K.; Mushiroda, T.; Nakajima, M.; Ozaki, K.; et al. A functional variant in ZNF512B is associated with susceptibility to amyotrophic lateral sclerosis in Japanese. Hum. Mol. Genet. 2011, 20, 3684–3692. [Google Scholar] [CrossRef] [Green Version]
  49. Reber, S.; Stettler, J.; Filosa, G.; Colombo, M.; Jutzi, D.; Lenzken, S.C.; Schweingruber, C.; Bruggmann, R.; Bachi, A.; Barabino, S.M.; et al. Minor intron splicing is regulated by FUS and affected by ALS-associated FUS mutants. EMBO J. 2016, 35, 1504–1521. [Google Scholar] [CrossRef]
  50. Ishihara, T.; Ariizumi, Y.; Shiga, A.; Kato, T.; Tan, C.-F.; Sato, T.; Miki, Y.; Yokoo, M.; Fujino, T.; Koyama, A.; et al. Decreased number of Gemini of coiled bodies and U12 snRNA level in amyotrophic lateral sclerosis. Hum. Mol. Genet. 2013, 22, 4136–4147. [Google Scholar] [CrossRef] [Green Version]
  51. Shibata, T.; Tokunaga, E.; Hattori, S.; Watari, K.; Murakami, Y.; Yamashita, N.; Oki, E.; Itou, J.; Toi, M.; Maehara, Y.; et al. Y-box binding protein YBX1 and its correlated genes as biomarkers for poor outcomes in patients with breast cancer. Oncotarget 2018, 9, 37216–37228. [Google Scholar] [CrossRef] [PubMed]
  52. Kuwano, M.; Shibata, T.; Watari, K.; Ono, M. Oncogenic Y-box binding protein-1 as an effective therapeutic target in drug-resistant cancer. Cancer Sci. 2019, 110, 1536–1543. [Google Scholar] [CrossRef] [PubMed]
  53. Nijssen, J.; Benitez, J.A.; Hoogstraaten, R.; Kee, N.; Hedlund, E. Axon-seq decodes the motor axon transcriptome and its modulation in response to ALS. Stem Cell Rep. 2018, 11, 321596. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Nasrin, F.; Rahman, M.A.; Masuda, A.; Ohe, K.; Takeda, J.-I.; Ohno, K. HnRNP C, YB-1 and hnRNP L coordinately enhance skipping of human MUSK exon 10 to generate a Wnt-insensitive MuSK isoform. Sci. Rep. 2014, 4, 6841. [Google Scholar] [CrossRef] [Green Version]
  55. Le Hir, H.; Saulière, J.; Wang, Z. The exon junction complex as a node of post-transcriptional networks. Nat. Rev. Mol. Cell Biol. 2016, 17, 41–54. [Google Scholar] [CrossRef]
  56. Jackson, K.L.; Dayton, R.D.; Orchard, E.A.; Ju, S.; Ringe, D.; Petsko, G.A.; Maquat, L.E.; Klein, R.L. Preservation of forelimb function by UPF1 gene therapy in a rat model of TDP-43-induced motor paralysis. Gene. Ther. 2015, 22, 20–28. [Google Scholar] [CrossRef] [Green Version]
  57. Barmada, S.J.; Ju, S.; Arjun, A.; Batarse, A.; Archbold, H.C.; Peisach, D.; Li, X.; Zhang, Y.; Tank, E.M.H.; Qiu, H.; et al. Amelioration of toxicity in neuronal models of amyotrophic lateral sclerosis by hUPF1. Proc. Natl. Acad. Sci. USA 2015, 112, 7821–7826. [Google Scholar] [CrossRef] [Green Version]
  58. Kamelgarn, M.; Chen, J.; Kuang, L.; Jin, H.; Kasarskis, E.J.; Zhu, H. ALS mutations of FUS suppress protein translation and disrupt the regulation of nonsense-mediated decay. Proc. Natl. Acad. Sci. USA 2018, 115, E11904–E11913. [Google Scholar] [CrossRef] [Green Version]
  59. Xu, W.; Bao, P.; Jiang, X.; Wang, H.; Qin, M.; Wang, R.; Wang, T.; Yang, Y.; Lorenzini, I.; Liao, L.; et al. Reactivation of nonsense-mediated mRNA decay protects against C9orf72 dipeptide-repeat neurotoxicity. Brain 2019, 142, 1349–1364. [Google Scholar] [CrossRef] [Green Version]
  60. Sun, Y.; Eshov, A.; Guo, J.U. C9orf72 Dipeptide Repeats Inhibit UPF1-Mediated RNA Decay Independent of Stress Granule Formation. BioRxiv 2019, 623769. [Google Scholar] [CrossRef] [Green Version]
  61. Hautbergue, G.M.; Castelli, L.M.; Ferraiuolo, L.; Sanchez-Martinez, A.; Cooper-Knock, J.; Higginbottom, A.; Lin, Y.-H.; Bauer, C.S.; Dodd, J.E.; Myszczynska, M.A.; et al. SRSF1-dependent nuclear export inhibition of C9ORF72 repeat transcripts prevents neurodegeneration and associated motor deficits. Nat. Commun. 2017, 8, 16063. [Google Scholar] [CrossRef]
  62. Balendra, R.; Isaacs, A.M. C9orf72-mediated ALS and FTD: Multiple pathways to disease. Nat. Rev. Neurol. 2018, 14, 544–558. [Google Scholar] [CrossRef] [PubMed]
  63. Mihevc, S.P.; Baralle, M.; Buratti, E.; Rogelj, B. TDP-43 aggregation mirrors TDP-43 knockdown, affecting the expression levels of a common set of proteins. Sci. Rep. 2016, 6, 33996. [Google Scholar] [CrossRef] [PubMed]
  64. Bakkar, N.; Kovalik, T.; Lorenzini, I.; Spangler, S.; Lacoste, A.; Sponaugle, K.; Ferrante, P.; Argentinis, E.; Sattler, R.; Bowser, R. Artificial intelligence in neurodegenerative disease research: Use of IBM Watson to identify additional RNA-binding proteins altered in amyotrophic lateral sclerosis. Acta Neuropathol. 2018, 135, 227–247. [Google Scholar] [CrossRef] [Green Version]
  65. Suzuki, H.; Shibagaki, Y.; Hattori, S.; Matsuoka, M. Nuclear TDP-43 causes neuronal toxicity by escaping from the inhibitory regulation by hnRNPs. Hum. Mol. Genet. 2015, 24, 1513–1527. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Rizzo, F.; Nizzardo, M.; Vashisht, S.; Molteni, E.; Melzi, V.; Taiana, M.; Salani, S.; Santonicola, P.; Di Schiavi, E.; Bucchia, M.; et al. Key role of SMN/SYNCRIP and RNA-Motif 7 in spinal muscular atrophy: RNA-Seq and motif analysis of human motor neurons. Brain 2019, 142, 276–294. [Google Scholar] [CrossRef] [PubMed]
  67. Das, R.; Dufu, K.; Romney, B.; Feldt, M.; Elenko, M.; Reed, R. Functional coupling of RNAP II transcription to spliceosome assembly. Genome Res. 2006, 20, 1100–1109. [Google Scholar] [CrossRef] [Green Version]
  68. Yu, Y.; Reed, R. FUS functions in coupling transcription to splicing by mediating an interaction between RNAP II and U1 snRNP. Proc. Natl. Acad. Sci. USA 2015, 112, 8608–8613. [Google Scholar] [CrossRef] [Green Version]
  69. Schwartz, J.C.; Ebmeier, C.C.; Podell, E.R.; Heimiller, J.; Taatjes, D.J.; Cech, T.R. FUS binds the CTD of RNA polymerase II and regulates its phosphorylation at Ser2. Genes Dev. 2012, 26, 2690–2695. [Google Scholar] [CrossRef] [Green Version]
  70. Zhao, R.Y.; Ni, Z.; Pu, S.; Zhong, G.; Schmitges, F.W.; Braunschweig, U.; Blencowe, B.J.; Greenblatt, J.F.; Zhao, Y. Regulation of transcription termination by FUS and TDP-43. BioRXiv 2019, 788778. [Google Scholar] [CrossRef]
  71. Boehringer, A.; Garcia-Mansfield, K.; Singh, G.; Bakkar, N.; Pirrotte, P.; Bowser, R. ALS Associated Mutations in Matrin 3 Alter Protein-Protein Interactions and Impede mRNA Nuclear Export. Sci. Rep. 2017, 7, 14529. [Google Scholar] [CrossRef] [PubMed]
  72. Chanarat, S.; Sträßer, K. Splicing and beyond: The many faces of the Prp19 complex. Biochim. Biophys. Acta (BBA) 2013, 1833, 2126–2134. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  73. Kalmar, B.; Greensmith, L. Cellular Chaperones As Therapeutic Targets in ALS to Restore Protein Homeostasis and Improve Cellular Function. Front. Mol. Neuro Sci. 2017, 10, 251. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  74. Coyne, A.N.; Lorenzini, I.; Chou, C.-C.; Torvund, M.; Rogers, R.S.; Starr, A.; Zaepfel, B.L.; Levy, J.; Johannesmeyer, J.; Schwartz, J.C.; et al. Post-transcriptional Inhibition of Hsc70-4/HSPA8 Expression Leads to Synaptic Vesicle Cycling Defects in Multiple Models of ALS. Cell Rep. 2017, 21, 110–125. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  75. Battle, D.; Kasim, M.; Yong, J.; Lotti, F.; Lau, C.-K.; Mouaikel, J.; Zhang, Z.; Han, K.; Wan, L.; Dreyfuss, G. The SMN Complex: An Assembly Machine for RNPs. Cold Spring Harb. Symp. Quant. Boil. 2006, 71, 313–320. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  76. Gubitz, A. The SMN complex. Exp. Cell Res. 2004, 296, 51–56. [Google Scholar] [CrossRef]
  77. Morello, G.; Guarnaccia, M.; Spampinato, A.G.; La Cognata, V.; D’Agata, V.; Cavallaro, S. Copy Number Variations in Amyotrophic Lateral Sclerosis: Piecing the Mosaic Tiles Together through a Systems Biology Approach. Mol. Neurobiol. 2018, 55, 1299–1322. [Google Scholar] [CrossRef] [Green Version]
  78. Zou, T.; Ilangovan, R.; Yu, F.; Xu, Z.; Zhou, J. SMN Protects cells against mutant SOD1 toxicity by increasing chaperone activity. Biochem. Biophys. Res. Commun. 2007, 364, 850–855. [Google Scholar] [CrossRef] [Green Version]
  79. Rodríguez-Muela, N.; Litterman, N.K.; Norabuena, E.M.; Mull, J.L.; Galazo, M.J.; Sun, C.; Ng, S.-Y.; Makhortova, N.R.; White, A.; Lynes, M.M.; et al. Single-Cell Analysis of SMN Reveals Its Broader Role in Neuromuscular Disease. Cell Rep. 2017, 18, 1484–1498. [Google Scholar] [CrossRef]
  80. Shefer, K.; Sperling, J.; Sperling, R. The Supraspliceosome—A Multi-Task Machine for Regulated Pre-mRNA Processing in the Cell Nucleus. Comput. Struct. Biotechnol. J. 2014, 11, 113–122. [Google Scholar] [CrossRef] [Green Version]
  81. Azubel, M.; Habib, N.; Sperling, R.; Sperling, J. Native Spliceosomes Assemble with Pre-mRNA to Form Supraspliceosomes. J. Mol. Boil. 2006, 356, 955–966. [Google Scholar] [CrossRef] [PubMed]
  82. Barmada, S.J. Linking RNA Dysfunction and Neurodegeneration in Amyotrophic Lateral Sclerosis. Neurotherapeutics 2015, 12, 340–351. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  83. Yu, H.; Kim, P.M.; Sprecher, E.; Trifonov, V.; Gerstein, M. The Importance of Bottlenecks in Protein Networks: Correlation with Gene Essentiality and Expression Dynamics. PLoS Comput. Boil. 2007, 3, e59. [Google Scholar] [CrossRef] [PubMed]
  84. Van Rheenen, W.; Diekstra, F.P.; Harschnitz, O.; Westeneng, H.-J.; Van Eijk, K.R.; Saris, C.G.J.; Groen, E.J.N.; Van Es, M.A.; Blauw, H.M.; Van Vught, P.W.J.; et al. Whole blood transcriptome analysis in amyotrophic lateral sclerosis: A biomarker study. PLoS ONE 2018, 13, e0198874. [Google Scholar] [CrossRef] [PubMed]
  85. Guarino, A.M.; Di Mauro, G.; Ruggiero, G.; Geyer, N.; Delicato, A.; Foulkes, N.S.; Vallone, D.; Calabrò, V. YB-1 recruitment to stress granules in zebrafish cells reveals a differential adaptive response to stress. Sci. Rep. 2019, 9, 9059. [Google Scholar] [CrossRef] [PubMed]
  86. Berchtold, D.; Battich, N.; Pelkmans, L. A Systems-Level Study Reveals Regulators of Membrane-less Organelles in Human Cells. Mol. Cell 2018, 72, 1035–1049.e5. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  87. Abrakhi, S.; Kretov, D.A.; Desforges, B.; Dobra, I.; Bouhss, A.; Pastré, D.; Hamon, L. Nanoscale Analysis Reveals the Maturation of Neurodegeneration-Associated Protein Aggregates: Grown in mRNA Granules then Released by Stress Granule Proteins. ACS Nano 2017, 11, 7189–7200. [Google Scholar] [CrossRef]
  88. Morello, G.; Spampinato, A.G.; Conforti, F.L.; D’Agata, V.; Cavallaro, S. Selection and Prioritization of Candidate Drug Targets for Amyotrophic Lateral Sclerosis through a Meta-Analysis Approach. J. Mol. Neurosci. 2017, 61, 563–580. [Google Scholar] [CrossRef] [Green Version]
  89. Morello, G.; Cavallaro, S. Transcriptional analysis reveals distinct subtypes in amyotrophic lateral sclerosis: Implications for personalized therapy. Future Med. Chem. 2015, 7, 1335–1359. [Google Scholar] [CrossRef] [Green Version]
  90. Morello, G.; Conforti, F.L.; Parenti, R.; D’Agata, V.; Cavallaro, S. Selection of Potential Pharmacological Targets in ALS Based on Whole-Genome Expression Profiling. Curr. Med. Chem. 2015, 22, 2004–2021. [Google Scholar] [CrossRef]
  91. Yin, W.; Rogge, M. TargetingRNA: A Transformative Therapeutic Strategy. Clin. Transl. Sci. 2019, 12, 98–112. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  92. Nussbacher, J.K.; Tabet, R.; Yeo, G.W.; Lagier-Tourenne, C. Disruption of RNA Metabolism in Neurological Diseases and Emerging Therapeutic Interventions. Neuron 2019, 102, 294–320. [Google Scholar] [CrossRef] [PubMed]
Cells 09 00159 g001 550
Figure 1. The flow chart illustrates the analysis methods used in the present study to characterize the differentially expressed genes encoding splicing factors (left green box) and investigate their functional significance (right orange box).
Figure 1. The flow chart illustrates the analysis methods used in the present study to characterize the differentially expressed genes encoding splicing factors (left green box) and investigate their functional significance (right orange box).
Cells 09 00159 g001
Cells 09 00159 g002 550
Figure 2. Venn diagram of differentially expressed splicing factor genes in tissue districts (motor cortex and spinal cord) of control and SALS subtype patients. Detailed information for the lists of gene probes is provided in Supplementary Data S2.
Figure 2. Venn diagram of differentially expressed splicing factor genes in tissue districts (motor cortex and spinal cord) of control and SALS subtype patients. Detailed information for the lists of gene probes is provided in Supplementary Data S2.
Cells 09 00159 g002
Cells 09 00159 g003 550
Figure 3. Supervised hierarchical clustering was used to visualize the gene expression changes in differentially expressed genes (302 entities) in the motor cortexes and spinal cords of control and ALS (SALS1 and SALS2 subtypes) patients. As shown in the color bar, red indicates up-regulation, and blue down-regulation.
Figure 3. Supervised hierarchical clustering was used to visualize the gene expression changes in differentially expressed genes (302 entities) in the motor cortexes and spinal cords of control and ALS (SALS1 and SALS2 subtypes) patients. As shown in the color bar, red indicates up-regulation, and blue down-regulation.
Cells 09 00159 g003
Cells 09 00159 g004 550
Figure 4. Venn diagrams of the total number of deregulated genes in SALS subtypes versus controls (a) and in motor cortex versus spinal cord (c). Tables show the numbers of up and down-regulated genes in SALS subtypes versus controls (b) and in motor cortex versus spinal cord (d), respectively.
Figure 4. Venn diagrams of the total number of deregulated genes in SALS subtypes versus controls (a) and in motor cortex versus spinal cord (c). Tables show the numbers of up and down-regulated genes in SALS subtypes versus controls (b) and in motor cortex versus spinal cord (d), respectively.
Cells 09 00159 g004
Cells 09 00159 g005 550
Figure 5. PPI networks in (a) SALS1 and (b) SALS2 cortexes. To better visualize the generated networks we used the Force Atlas layout algorithm. Sphere size is proportional to the degree of connection, whereas the color (red for up-regulated and green for down-regulated in SALS subtypes) represents the expression logFC value. Details about nodes properties are listed in Supplementary Data S4.
Figure 5. PPI networks in (a) SALS1 and (b) SALS2 cortexes. To better visualize the generated networks we used the Force Atlas layout algorithm. Sphere size is proportional to the degree of connection, whereas the color (red for up-regulated and green for down-regulated in SALS subtypes) represents the expression logFC value. Details about nodes properties are listed in Supplementary Data S4.
Cells 09 00159 g005
Cells 09 00159 g006 550
Figure 6. Spinal cord-specific protein–protein interaction (PPI) networks are shown for both (a) SALS1 and (b) SALS2. To better visualize generated networks, we used the Auto layout algorithm. Sphere size is proportional to the degree of connection, whereas the color (red for up-regulated and green for down-regulated in SALS subtypes) represents the expression logFC value. Details about node properties are listed in Supplementary Data S4.
Figure 6. Spinal cord-specific protein–protein interaction (PPI) networks are shown for both (a) SALS1 and (b) SALS2. To better visualize generated networks, we used the Auto layout algorithm. Sphere size is proportional to the degree of connection, whereas the color (red for up-regulated and green for down-regulated in SALS subtypes) represents the expression logFC value. Details about node properties are listed in Supplementary Data S4.
Cells 09 00159 g006
Cells 09 00159 g007 550
Figure 7. Literature-based network of tissue-specific deregulated genes involved in stress granule formations generated by the Genomatix Pathways System (GePS). The hierarchical layout has been used to highlight the main direction or information flow of the network. The legend below reports information about elements in figure.
Figure 7. Literature-based network of tissue-specific deregulated genes involved in stress granule formations generated by the Genomatix Pathways System (GePS). The hierarchical layout has been used to highlight the main direction or information flow of the network. The legend below reports information about elements in figure.
Cells 09 00159 g007
Table
Table 1. Number of samples processed for sporadic amyotrophic lateral sclerosis (SALS) patients’ molecular clusters and regions on the central nervous system (CNS).
Table 1. Number of samples processed for sporadic amyotrophic lateral sclerosis (SALS) patients’ molecular clusters and regions on the central nervous system (CNS).
Molecular ClustersMotor CortexSpinal Cord
SALS11817
SALS21313
CTRL1010
TOT4140
Table
Table 2. Enriched gene ontologies (GOs) of differentially expressed genes (DEGs) in both disease state and tissue types pairwise comparisons.
Table 2. Enriched gene ontologies (GOs) of differentially expressed genes (DEGs) in both disease state and tissue types pairwise comparisons.
Enriched GO Cellular ComponentMajor Disease State Changes in Tissue TypesMajor Tissue Type Changes in Disease States
SALS1 Motor Cortex vs ControlSALS2 Motor Cortex vs ControlSALS1 Spinal Cord vs ControlSALS2 Spinal Cord vs ControlCtrlSALS1SALS2
U1 snRNP Cells 09 00159 i0011.55 × 109 Cells 09 00159 i0023.70 × 107 Cells 09 00159 i0015.18 × 108
A2AF Cells 09 00159 i0021.24 × 102 Cells 09 00159 i0032.87× 102 Cells 09 00159 i0032.82× 102
U2 snRNP Cells 09 00159 i0011.59 × 109 Cells 09 00159 i0028.44 × 109 Cells 09 00159 i0011.38 × 109
U2-type prespliceosome Cells 09 00159 i0018.35 × 1010 Cells 09 00159 i0028.33 × 106 Cells 09 00159 i0017.23 × 1010
U2-type precatalytic spliceosome Cells 09 00159 i0022.84 × 106 Cells 09 00159 i0018.57 × 1028 Cells 09 00159 i0022.80 × 1015 Cells 09 00159 i0012.89 × 1026
U2-type catalytic step 1 spliceosome Cells 09 00159 i0024.21× 102 Cells 09 00159 i0011.46 × 104 Cells 09 00159 i0022.65 × 103 Cells 09 00159 i0014.01 × 106
U2-type catalytic step 2 spliceosome Cells 09 00159 i0022.29 × 107 Cells 09 00159 i0016.89 × 1013 Cells 09 00159 i0026.47 × 1011 Cells 09 00159 i0012.09 × 1011
U4 snRNP Cells 09 00159 i0013.36 × 103 Cells 09 00159 i0013.22 × 103
U5 snRNP Cells 09 00159 i0026.80 × 107 Cells 09 00159 i0011.54 × 108 Cells 09 00159 i0025.47 × 106 Cells 09 00159 i0015.55 × 107
U6 snRNP Cells 09 00159 i0011.38 × 106 Cells 09 00159 i0011.26 × 106
U4/U6.U5 tri-snRNP complex Cells 09 00159 i0023.39 × 104 Cells 09 00159 i0011.61 × 1017 Cells 09 00159 i0028.01 × 1011 Cells 09 00159 i0015.64 × 1016
U12-type spliceosomal complex Cells 09 00159 i0023.75× 106 Cells 09 00159 i0016.35 × 109 Cells 09 00159 i0034.30 × 102 Cells 09 00159 i0021.09 × 106 Cells 09 00159 i0011.81 × 107
Exon-exon junction complex Cells 09 00159 i0029.41 × 105 Cells 09 00159 i0012.07 × 109 Cells 09 00159 i0024.60 × 107 Cells 09 00159 i0015.09 × 1011
CRD-mediated mRNA stability complex Cells 09 00159 i0021.90 × 104 Cells 09 00159 i0011.02 × 103 Cells 09 00159 i0022.70 × 102 Cells 09 00159 i0012.35 × 105
RNA polymerase II, core complex Cells 09 00159 i0012.84 × 1010 Cells 09 00159 i0021.65 × 104 Cells 09 00159 i0012.45 × 1010
Prp19 complex Cells 09 00159 i0012.44 × 104 Cells 09 00159 i0032.38 × 104
SMN-Sm protein complex Cells 09 00159 i0011.90 × 105 Cells 09 00159 i0025.93 × 103 Cells 09 00159 i0011.64 × 105
SMN complex Cells 09 00159 i0013.31 × 103 Cells 09 00159 i0013.18 × 103
pICln-Sm protein complex Cells 09 00159 i0012.68 × 105 Cells 09 00159 i0012.39 × 105
Supraspliceosomal complex Cells 09 00159 i0012.84 × 104 Cells 09 00159 i0021.64 × 104 Cells 09 00159 i0011.43 × 102
Spliceosomal complex Cells 09 00159 i0011.11 × 104
Catalytic step 2 spliceosome Cells 09 00159 i0016.15 × 103 Cells 09 00159 i0019.80 × 103
Post-mRNA release spliceosomal complex Cells 09 00159 i0016.29 × 103 Cells 09 00159 i0023.74 × 103 Cells 09 00159 i0012.34 × 104
For each enriched GO, the Table shows GO cellular component name and the enrichment p-value. Red or green arrow indicates if the GO is up or down-regulated.

Share and Cite

MDPI and ACS Style

La Cognata, V.; Gentile, G.; Aronica, E.; Cavallaro, S. Splicing Players Are Differently Expressed in Sporadic Amyotrophic Lateral Sclerosis Molecular Clusters and Brain Regions. Cells 2020, 9, 159. https://doi.org/10.3390/cells9010159

AMA Style

La Cognata V, Gentile G, Aronica E, Cavallaro S. Splicing Players Are Differently Expressed in Sporadic Amyotrophic Lateral Sclerosis Molecular Clusters and Brain Regions. Cells. 2020; 9(1):159. https://doi.org/10.3390/cells9010159

Chicago/Turabian Style

La Cognata, Valentina, Giulia Gentile, Eleonora Aronica, and Sebastiano Cavallaro. 2020. "Splicing Players Are Differently Expressed in Sporadic Amyotrophic Lateral Sclerosis Molecular Clusters and Brain Regions" Cells 9, no. 1: 159. https://doi.org/10.3390/cells9010159

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop