Transcriptional profiling of differentially vulnerable motor neurons at pre-symptomatic stage in the Smn 2b/- mouse model of spinal muscular atrophy
Acta Neuropathologica Communications volume 3, Article number: 55 (2015)
The term motor neuron disease encompasses a spectrum of disorders in which motor neurons are the lost. Importantly, while some motor neurons are lost early in disease and others remain intact at disease end-stage. This creates a valuable experimental paradigm to investigate the factors that regulate motor neuron vulnerability. Spinal muscular atrophy is a childhood motor neuron disease caused by mutations or deletions in the SMN1 gene. Here, we have performed transcriptional analysis on differentially vulnerable motor neurons from an intermediate mouse model of Spinal muscular atrophy at a presymptomatic time point.
We have characterised two differentially vulnerable populations, differing in the level neuromuscular junction loss. Transcriptional analysis on motor neuron cell bodies revealed that reduced Smn levels correlate with a reduction of transcripts associated with the ribosome, rRNA binding, ubiquitination and oxidative phosphorylation. Furthermore, P53 pathway activation precedes neuromuscular junction loss, suggesting that denervation may be a consequence, rather than a cause of motor neuron death in Spinal muscular atrophy. Finally, increased vulnerability correlates with a decrease in the positive regulation of DNA repair.
This study identifies pathways related to the function of Smn and associated with differential motor unit vulnerability, thus presenting a number of exciting targets for future therapeutic development.
Motor neuron diseases (MNDs) are a heterogeneous group of neurodegenerative disorders that are caused by a diverse array of factors, including both genetic and sporadic. The clinical severity can also vary widely, but MNDs are frequently very severe, causing fatality within months to years of diagnosis. Despite a range of causes and severities, motor neuron diseases are united by the common vulnerability of motor neurons. The reason why these cells are selectively vulnerable to the genetic or environmental insult is unknown. Importantly, however, not all motor neurons are equally affected . There are some pools of motor neuron in which there are high levels of pathology throughout the motor unit, and high levels of motor neuron loss [6, 22, 27, 16]. In other pools in the same individual, there will be minimal evidence of motor unit pathology, even at late stages of disease. The reasons for this selective vulnerability are currently unknown, however this phenomenon creates an exciting and valuable opportunity to investigate the factors governing motor neuron vulnerability and pathology. By comparing different motor neuron groups we can investigate the molecular mechanisms underlying the disease, the molecular mechanisms underlying motor neuron pathology and the mechanisms that regulate motor neuron vulnerability.
Spinal muscular atrophy (SMA) is a childhood motor neuron disease caused by mutations and deletions within the survival motor neuron 1 gene (SMN1). There is a second partially functional copy of the SMN gene, termed SMN2. Due to a point mutation, SMN2 predominantly produces a form of SMN that lacks exon 7, which is rapidly degraded. This gene therefore produces only low levels of the full length SMN transcript. The copy number of SMN2 can vary, making it an important phenotypic modifier, with disease severity closely correlating to the number of copies of the SMN2 gene. For this reason, SMA has a range of clinical severities, which are classified into types 0 to 4 based on age of onset and motor milestones achieved. Despite a degree of controversy in the field as to the cell types affected in SMA, it is clear that motor neurons are particularly vulnerable to reduced SMN levels. Indeed, SMA is typified by the loss of lower motor neurons and the atrophy of associated skeletal musculature.
In SMA patients there is evidence of pathology at the neuromuscular junction (NMJ) [37, 53]. Indeed in models of SMA, structural pathology in the motor neuron is first observed at the NMJ . NMJs are lost early in the disease, with denervation evident in some muscles prior to symptom onset . Additional defects at the NMJ include neurofilament accumulation, poor terminal arborisation, delayed post-synaptic maturation, defects in calcium handing and a disruption in synaptic vesicle release [10, 25, 28, 29, 34, 40, 44]. Importantly, not all NMJs are equally affected. This has been well characterised in mouse models of SMA, where there appears to be high levels of NMJ loss in some muscles and very low levels of NMJ loss in others [39, 40, 33]. This selective vulnerability has even been observed within single muscles, where there are areas of high levels of NMJ loss, and other areas showing no apparent denervation . The reasons for this differential vulnerability are unclear. One study has investigated the various aspects of differentially vulnerable motor units, including motor unit size, muscle fibre type, NMJ size, branching patterns and Terminal Schwann cell number but found no correlation with differential vulnerability . There is some evidence that differential vulnerability might correlate with a differential sprouting competence, presenting the possibility that relative plasticity can be an important modifier. Furthermore, recent work has shown that there are defects in synaptic remodelling at the NMJ in a mouse model of SMA . However, this has not been shown to be a definitive modulator of vulnerability. The reasons for the variability in NMJ vulnerability are therefore currently unclear.
The reasons why motor neurons are selectively vulnerable to a reduction in Smn are currently undefined and the cellular functions for Smn have long been debated. The best described function for Smn is in pre-mRNA splicing [8, 57]. Growing evidence suggests that a reduction in Smn does not result in widespread splicing defects however it has been suggested that reduced Smn levels may cause splicing defects in a small subset of RNAs [36, 3]. Other roles for Smn have also been proposed, including transport of mRNAs  and a potential role as a translational regulator . Ultimately, the precise cellular function of Smn that, when compromised in SMA, causes motor neuron pathology is unknown.
SMA research has been greatly facilitated by a number of mouse models . These models are traditionally concentrated at the severe end of the spectrum, with a life expectancy of less than 14 days of age [31, 38, 23]. Despite their limitations, they have been instrumental in understanding the pathophysiology of SMA and in therapeutic development. The development of the Smn 2B/- mouse model of SMA has been an important addition to existing mouse models [5, 20]. This mouse was created by introducing a mutation in a splice enhancer site in the murine Smn gene, resulting in production of 10-15 % of normal Smn levels . This mouse model displays all the hallmark pathologies associated with SMA but, with a phenotypic onset of around 10 days and a life span of around 28 days, has a slightly milder phenotype than most other models . Importantly, unlike other mouse models, the Smn 2B/- mouse has a prolonged pre-symptomatic time period, thus allowing analysis of the events preceding motor neuron loss in SMA.
In this study, we have used the Smn 2B/- mouse model of SMA to investigate the transcriptional changes that occur in differentially vulnerable motor neurons. We have identified and defined muscles with differential levels of NMJ pathology and identified the time points associated with the onset of degeneration. At the pre-degenerative time points, we have used retrograde tracers to identify the motor neuron cell bodies which correspond to these differentially vulnerable NMJs and isolated them by laser capture micro-dissection. Vulnerable and less vulnerable motor neuron populations were isolated from both SMA and wild-type control mice. We then performed RNAseq on these differentially vulnerable motor neuron pools. Comparative bioinformatics analysis and functional clustering was performed to identify the transcriptional changes that correlate with reduced Smn levels, increased vulnerability and increased pathology. We demonstrate that reduced Smn levels are correlative with a decrease in transcripts implicated in ribosome, rRNA binding, ubiquitin and oxidative phosphorylation. We demonstrate a selective up regulation of cell death pathways in selectively vulnerable motor neurons, and demonstrate that an increase in these transcripts is observed prior to NMJ loss. Finally, we demonstrate that there is a decrease in markers of DNA repair in selectively vulnerable motor neurons. Overall this work details a four way comparative analysis of differentially vulnerable motor neurons using high resolution transcriptional profiling.
Materials and methods
The Smn 2B/- mice  were established in our laboratory and maintained in the University of Ottawa vivarium on a C57BL/6 x CD1 hybrid background. Mice were sacrificed by exposure to rising CO2 levels or by cervical dislocation. All procedures were performed in accordance with institutional guidelines (Animal Care and Veterinary Services, University of Ottawa). Smn −/− ;SMN2 mice  were maintained in the animal facilities at the University of Edinburgh and were sacrificed by overdose of anaesthetic. Tissues were provided by generous agreement with Prof Tom Gillingwater. All procedures were carried out in accordance with the procedures approved and licenced by the Home Office, United Kingdom.
Motor neuron labelling
Rhodamine conjugated dextran (RhDextran; 3000 MW; Molecular Probes) was administered under general anaesthesia . For immunostaining, a fixable analogue of the RhDextran was used. Mothers were pre-dosed with buprenorphine prior to surgery. Mice were anesthetised by inhalation of isofluorane (2 % in 1:1 N2O/O2). For abdominal muscle labelling, a small incision between the last rib and xyphoid process was made and 5 μl of 5 % RhDextran was injected into the space between the external oblique muscle and the transversus abdominis muscle. For cranial muscle labelling, a small incision was made on the back of the neck, and 5 μl of RhDextran was injected into the space between the levator auris longus and auricularis superior muscles. Mice were allowed to recover from anaesthetic before being returned to standard cages. For muscle analysis, mice we sacrificed 24 h later and muscles were dissected, fixed in paraformaldehyde (PFA) and mounted on slides. For spinal cords and brainstems, mice were sacrificed 48 h later and spinal cords and brainstems were removed and snap frozen in liquid nitrogen in 50 % Tissue Tek O.C.T, 15 % sucrose in PBS.
Laser capture microdissection of RhDextran labelled motor neurons
Spinal cords or brainstems were frozen and embedded as described above. Tissues were then sectioned using a cryostat at a thickness of 12 μm and mounted onto uncoated, uncharged glass slides and immediately frozen at −20 °C before being stored at −80 °C. Fluorescent motor neurons were then identified and imaged on an inverted epifluorescent microscope. During this process, slides were kept frozen using a freezing aerosol. Images were assembled and montaged using Adobe Photoshop and printed onto overhead transparencies. In the interim, slides were nissl stained and dehydrated as follows: 75 % EtOH, 30 s; ddH2O, 30 s; 1 % Toludine blue, 60 s; ddH2O, 30 s; 75 % EtOH, 30 s; 95 % EtOH, 30 s; 100 % EtOH x 3, 60 s; Xylene, 5 min; Air dry, 5 min. Steps were taken throughout to minimise RNAse exposure and DEPC water was used throughout.
LCM was performed using an Arcturus XT laser capture microscope from Applied Biosystems. Labelled motor neurons were identified by realigning the printed images of the fluorescent motor neurons with the nissl stained motor neurons on the LCM computer screen. This step was necessary due to the water soluble nature of the RhDextran, which is lost during LCM processing. Laser settings were optimised for motor neuron cell body size. Captured cells were snap frozen in qiazol (Qiagen) and stored at −80 °C until extraction of RNA.
RNA was extracted using a Qiagen microRNA micro RNeasy kit as per manufacturer’s instructions. RNA was amplified using the ovation RNAseq system version 2 from Nugen according to manufacturer instructions and template DNA library construction was performed with the Encore Rapid Library System (Nugen). Four cycles of PCR were performed using reagents from the Ovation Ultralow kit (Nugen). After that, 36 cycles of single-end sequencing was performed with the Genome Analyzer IIx (Illumina). Reads were mapped to the mouse mm9 assembly using tophat (v1.4.0) using the transcripts from Ensembl release 67 to guide mapping. Quality control was performed using RNAseQC and FastQC. Relative transcript levels were compared using CuffDiff software v1.3 using the UCSC transcript model. Significance was considered with a Q-value of less that 0.05. The statistic package R was used to identify transcriptional changes which were common in more than one data set i.e. identifying transcriptional changes occurring in SMAv and SMAr motor neurons compared to their respective wild-types. Fold changes were converted to log2. Therefore no change, which is equivalent to a one fold change is equal to log21, which equals 0. Therefore all numbers greater than 0 imply an up-regulation. All numbers less that 0 imply a down-regulation. As the biological consequences for a specific magnitude of change are not known, we have not applied any fold change restrictions to the data.
Functional clustering was performed using the functional annotation clustering tool available on the Database for Annotation, Visualization and Integrated Discovery (DAVID) online software from the National Institute of Allergy and Infectious Diseases . Using this tool, enrichment analysis was performed using over 40 annotation categories. An algorithm was then applied which clustered the functional annotations based upon the degree of co-associated genes. Only those functional annotations that had a significant non-adjusted P value, were clustered. Those clusters with an enrichment score of >1.3 are considered significant. The enrichment score is a Log10 representation of the non-adjusted P value. Therefore an enrichment score of 1.3 is reflective of a P value of 0.05. This software also provides Kyoto Encyclopaedia of Genes and Genomes (KEGG) pathway annotation to identify specific cellular pathways that contain an enrichment for differentially expressed genes.
For q-RT-PCR on laser captured motor neurons, motor neuron cell bodies were isolated and RNA was extracted and amplified as above and cDNA was diluted to approximately 100 ng/μl. Transcripts were selected for validation based on having a high relative expression (based upon the read counts from the RNAseq results, preferably a read value of >50) balanced with a large change in expression level (ideally >1.5 fold). For whole spinal cord, RNA was extracted using a micro RNeasy kit (Qiagen) and 1 ug of RNA was used to perform reverse transcription using the RT2 First Strand kit (Qiagen). SYBR gene based Q-RT-PCR was performed using pre-optimised primers purchased from BioRad. Amplification was performed using KAPA SYBR fast universal PCR mastermix as per manufacturer instructions on a BioRad CFX connect real-time PCR detection system. Relative gene expression was calculated using the 2-ΔΔcT formula . β-actin and Y-Whaz were used as reference genes.
For NMJ labelling, muscles were immediately dissected from recently sacrificed mice and fixed in 4 % PFA (Electron Microscopy Science) in PBS for 15 min. Post-synaptic AChRs were labelled with α-bungarotoxin (BTX) for 30 min. Muscles were permeabilised in 2 % Triton X-100 in PBS for 30 min, then blocked in 4 % bovine serum albumin (BSA)/1 % Triton X-100 in PBS for 30 min before incubation overnight in primary antibodies [Neurofilament (NF; 2H3) - Developmental Studies Hybridoma Bank; synaptic vesicle protein 2 (SV2) - Developmental Studies Hybridoma Bank; S100 – Dako; all 1:250] and visualised with Cy3-conjugated secondary antibodies [Cy3 goat anti-mouse; 1:250, Jackson]. Muscles were then whole-mounted in Dako Fluorescent mounting media. Images were taken with a Zeiss LSM-510 confocal microscope.
For spinal cord sections, spinal cords were removed from recently sacrificed mice and fixed overnight in 4 % PFA. Tissues were equilibrated in 30 % sucrose, embedded in 50 % Tissue Tek O.C.T, 15 % sucrose in PBS and sectioned at 12 μm on a cryostat. Sections were washed in PBS, permeabilised in 0.1 % Triton X-100 for 10 min and blocked in 4 % BSA for 30 min before exposure to primary antibodies (rabbit anti-pH2AX, Cell Signalling Technology; mouse anti-neurofilament heavy chain [SMI-32] BioLegend) overnight. Secondary antibodies (AlexaFluor 555 Goat anti-Rabbit, AlexaFluor 488 Rabbit anti-Mouse, both Life Technologies) were applied for a period of 2–4 h at a dilution of 1:200. Sections were counterstained in DAPI (Life Technologies) and NeuroTrace® 500/525 (Life Technologies) as per manufacturer instructions and mounted in mowoil® (Sigma Aldrich). Sections were imaged on a Zeiss LSM-510 confocal microscope. All laser settings were kept constant between images.
Quantification and statistics
The percentage of fully occupied endplates was determined by classifying each endplate in a given field of view either fully occupied (pre-synaptic terminal (SV2 and NF) completely overlies endplate (BTX)), partially occupied (pre-synaptic terminal only partially covers endplate (BTX)), or vacant (no pre-synaptic label overlies endplate). At least 4 fields of view were analysed per muscle totalling >100 endplates per muscle.
For quantification of the number of labelled motor neurons following RhDextran injection, longitudinal sections of the anterior horn of the thoracic spinal cord or pons and medulla of the brainstem were cut and mounted sequentially on glass slides. The number of fluorescent motor neurons that were visible was counted per section. Quantification was done visually using a Zeiss Axiovert 200 M microscope.
pH2AX quantification was done using Image-J software. The average background intensity was subtracted from each image, leaving only the pixels with an intensity above background levels. The number of pH2AX foci per neuron was then quantified. All quantification was done with the experimenter blind to the experimental group.
All data was assembled and analysed using Microsoft Excel and GraphPad Prism. All figures were assembled using Adobe Photoshop.
Results and discussion
Differential motor unit vulnerability in the Smn2B/- mouse model of SMA
Here we have used the Smn 2B/- mouse to investigate the transcriptional differences between differentially vulnerable motor units in SMA. To do this, we first identified two differentially vulnerable and experimentally accessible pools of motor neurons. We aimed to identify one pool of motor neurons in which we saw a high degree of NMJ loss, and another pool of motor neurons in which there was minimal or no evidence of NMJ loss at end stage of disease. For this study, it was important to us to use muscles that fulfilled the following criteria. Firstly, we aimed to use thin and flat muscles, which are easy to perform immunofluorescence on, thus allowing comprehensive analysis of the whole muscle. This was important so we could accurately quantify the level of morphological pathology present, and to ensure uniformity of pathology throughout the muscle. Secondly, muscles had to be experimentally accessible. This was important so tracer injections could be performed to specifically target the muscle group of interest with minimal trauma or risk of injection side effects, which may occur in the diaphragm or intercostal for example. Lastly, we felt it was preferable to use local groups of muscles, rather than single muscles. Clearly there is significant molecular variability between muscles. By using muscle groups, we minimised the contribution of inter-muscular variability. After investigating NMJ pathology in a range of muscles situated throughout the body, we defined our vulnerable group as those motor neurons innervating a group of abdominal muscles, specifically transversus abdominis (TVA), rectus abdominis (RA) and external oblique (EO). Collectively 20.5 % of NMJs were denervated by P21 in these muscles, and all these muscles display significant NMJ defects (Fig. 1). We observed defects such as neurofilament accumulation and poor terminal arborisation. Endplates appeared immature and there was a high number of vacant denervated endplates. Importantly, pathology was uniform throughout each muscle (data not shown). For our less vulnerable group of motor neuron, we chose to use motor neurons innervating a group of cranial muscles, encompassing levator auris longus (LAL), auricularis superior (AS) and adductor auris longus (AAL). Collectively, 97.9 % of NMJs were fully occupied at P21 in these muscles (Fig. 1b). Pre-synaptic terminals appeared more elaborate compared to those in abdominal muscles, and endplates displayed more complexity consistent with a more mature phenotype. Again, pathology was generally consistent throughout each muscle. One exception to this was a slight increase in NMJ pathology in the rostral band of the LAL muscle compared to the caudal band, where we observed 95.6 % of denervated endplates compared to 98.3 % in the LALc. Intriguingly, this is actually the opposite to what is observed in more severe mouse models of SMA . Importantly, however, this denervation was still comparatively mild compared to abdominal muscles.
It is important to note that our goal was to investigate the transcriptional differences between differentially vulnerable motor neurons prior to the onset of degeneration. We therefore aimed to identify a time point just prior to the onset of NMJ loss in vulnerable muscles. We performed a time course analysis of NMJ loss in the TVA, RA, IO and EO muscles. This revealed that the latest time point at which there was no evidence of NMJ loss was P10 (Fig. 2). Analysis of NMJs at P10 revealed no denervation in any of the abdominal muscles analysed (Fig. 2a). We have therefore defined this as a pre-degenerative time point and aimed to use this time point for all subsequent transcriptional analysis.
Motor unit tracing
To identify the motor neuron cell bodies that correspond to the differentially vulnerable pools of NMJs, we used rhodamine-conjugated dextran as a retrograde tracer. For abdominal muscle, RhDextran was injected into the space between the TVA and RA/EO muscles. For cranial muscles, RhDextran was injected into the space between the LAL and AS/AAL muscles. For initial experiments, mice were sacrificed 24 h later and muscles were analysed for the presence of RhDextran. Analysis of abdominal muscles revealed strong staining in the superior parts of the EO, RA and TVA (data not shown). Importantly, the dye remained relatively local to the site of injection, and did not leach to surrounding intrinsic back or appendicular muscles. Equivalent analysis was performed in the cranial muscles. Strong staining was observed in the LAL, AS and AAL muscles and again dye was observed to remain local to the site of injection, with minimal leakage to surrounding muscles or to the contralateral side. This analysis confirmed that we were labeling the targeted muscles specifically and robustly.
For subsequent experiments, mice were sacrificed 48 h following RhDextran injection and spinal cords and brainstems were analysed for the presence of RhDextran staining. This revealed a discrete cluster of motor neurons in the thoracic spinal cord, which correspond to the motor neurons innervating the abdominal muscles, and a cluster of labeled motor neurons in the facial nucleus of the brainstem, corresponding to those motor neurons innervating the cranial muscles (Fig. 3). Quantification of the number of labeled motor neurons between Smn 2B/- mice and wild-type control animals revealed no difference between groups (Fig. 3d). This suggests that labeling efficiency is not compromised in Smn 2B/- mice. This was important to consider given previous work reporting defects in axon transport in SMA mouse models .
Transcriptional analysis of differentially vulnerable motor neurons
Following identification of suitable differentially vulnerable motor neuron populations, and validation of the RhDextran labeling methodology, laser capture micro-dissection (LCM) was used to isolate the cell bodies corresponding to these neurons. The motor neuron cell bodies were isolated from Smn 2B/- and wild-type mice, with approximately 250 motor neurons from 2–3 mice per biological replicate. This resulted in four experimental groups. Vulnerable (abdominal) motor neurons from Smn 2B/- and wild-type mice (SMAv and WTv respectively) and less vulnerable (cranial) motor neurons from Smn 2B/- and wild-type mice (SMAr and WTr), respectively. Transcriptional analysis was performed using RNAseq and expression profiles were compared between each experimental group (Additional file 1). For each comparison, around 17,500 genes were identified and compared. From this number, there were around 1200 to 1500 transcriptional changes between each of our four comparative groups (Fig. 4). Encouragingly, preliminary analysis of the data revealed a number of changes that have been previously implicated in pathology in SMA or in motor neuron vulnerability. For example, we observed a significant 4.8 Log2 fold up regulation of cyclin dependant kinase inhibitor 1A (CDKN1a) between SMAv and WTv motor neurons, and a 2.8 Log2 fold up regulation of CDKN1a between SMAr and WTr motor neurons. Interestingly, a significant up-regulation of CDKN1a has been reported in numerous other models of SMA, including in other mouse models of SMA  and embryonic stem cell derived motor neurons from SMA mouse models . We revealed a significant 1.3 Log2 fold down-regulation of chondrolectin (Chodl), which has previously been described in other mouse models of SMA [3, 60]. There was a significant 1.03 fold Log2 up-regulation of the Fused in Sarcoma transcript (Fus), which has been implicated in motor neuron pathology in amyotrophic lateral sclerosis (ALS), and the protein has been shown to interact with Smn [18, 59]. We also observed a significant 3.5 Log2 fold down regulation of insulin like growth factor 1 (IGF1) which has previously been reported in serum and liver of a SMA mouse model . Interestingly, we did not observe any change in Fus, IGF1 or Chodl in SMAr motor neurons compared to WTr. As these genes have also been implicated in motor neuron pathology in amyotrophic lateral sclerosis [9, 30, 52, 56], this perhaps indicates that these genes are involved in motor neuron pathology, which occurs as a secondary effect to the loss of Smn protein. These genes are therefore potential common downstream regulators of neuronal pathology in a variety of different motor neuron diseases.
What transcriptional changes occur in motor neurons when Smn is absent: transcriptional differences in vulnerable motor neurons between SMA and WT mice
We first aimed to use this data to investigate the cellular processes disrupted when Smn levels are decreased. We reasoned that Smn levels are reduced in both SMAv and SMAr motor neurons compared to their respective wild-type counterparts. We therefore looked for the transcriptional changes that occurred in both of these populations. In SMAr motor neurons, there were 1557 transcripts, which were statistically up or down regulated compared to WTr motor neurons. In SMAv motor neurons, there were 1277 transcripts, which were statistically up or down regulated compared to WTv motor neurons. Of these changes, 273 were found in both comparisons and the majority (222) showed a common directional change, with 121 being up-regulated and 101 being down regulated (Fig. 4, Fig. 5a, Additional file 2).
Functional clustering of these transcripts was performed using DAVID bioinformatics resources version 6.7. This revealed a number of significantly altered functional clusters (Table 1). As we were looking for transcripts implicated in the function of Smn, we decided to focus on transcripts that are down-regulated in SMAv and SMAr motor neurons. The most significant clusters which were down-regulated comprised of transcripts involved in the ribosome and cytosolic ribosome. This was closely followed by transcripts implicated in rRNA binding. Indeed there was considerable overlap between these three functional clusters, with the majority of rRNA binding proteins also implicated in the ribosome (Fig. 6). In order to confirm this down regulation of transcripts involved in rRNA binding and ribosome, we performed qPCR validation on new biological replicates of cDNA samples prepared from motor neurons isolated as described above, and compared transcript levels between SMAv and WTv samples. This confirmed a down-regulation of RPL7, RPLP1, RPS27a and MRPL20 (Fig. 6b). The only transcript that did not validate was LSM5.
The data presented above therefore indicated that there is a down regulation of factors involved in the ribosome, particularly in those involved in rRNA binding. Smn has long been established as an RNA binding protein, with work generally focused on its role as an mRNA binding protein, and its role on pre-mRNA splicing and mRNA transport. Less is known about its potential role in rRNA metabolism. It is easy to speculate that reduced Smn levels may disrupt the production or assembly for rRNA and lead to downstream defects in ribosomal function. Furthermore, Smn has recently been identified as a negative regulator of translation . Smn has also been shown to localise to the nucleolus and interact with non-ribosomal nucleolar proteins, suggesting that it has a role in ribonucleoprotein complex assembly [32, 54]. Interestingly, recent work has identified that mutations in genes involved in ribosome subunit assembly and maturation lead to motor neuron specific defects in Spinal Muscular Atrophy with Respiratory Distress (SMARD) [7, 12, 19]. This work therefore creates an intriguing parallel between SMA and SMARD and warrants further work into ribosome assembly and function in Smn depleted motor neurons.
The functional clustering also revealed a down regulation of transcripts involved in ubiquitin metabolism (Table 1, Fig. 7). The down regulation of ubiquitin (Ubb) was confirmed by qPCR (Fig. 7b). A disruption in ubiquitin homeostasis has recently been reported in SMA mouse models. Specifically, the authors observed a disruption in proteins implicated in ubiquitination in proteomic analysis on synapses from P1 SMA mice and a marked reduction in ubiquitin activating enzyme 1 (UBA1) levels in spinal cord and muscle from SMA mice . They suggest this down-regulation in UBA1 leads to an increase in beta-catenin levels, and inhibition of beta-catenin signalling rescued motor neuron specific defects in zebrafish and mouse models of SMA. This work therefore confirms the disruption of this pathway and promotes further work to dissect the mechanism by which a reduction in Smn levels lead to a disruption in ubiquitin homeostasis.
We also observed enrichment for transcripts involved in oxidative phosphorylation (Table 1, Fig. 8). Encouragingly, oxidative phosphorylation was also the top hit on a proteomic screen on synapses from a P1 SMA mouse model . Of particular note, we reveal a down regulation of haemoglobin beta 1 (Hbb-b1) which was confirmed by qPCR between SMAv and WTv motor neurons (Fig. 8b). Interestingly, we also noted a strong increase in Hbb-b1 levels in WTr motor neurons compared to WTv. This suggests that Hbb levels may be higher in less vulnerable motor neurons. This makes haemoglobin an exciting candidate as a regulator of motor neuron vulnerability. Indeed, recent work has demonstrated that haemoglobin subunits are expressed in neurons in both rats and humans [41, 43] and colocalised with the mitochondria inner membrane . They have been observed to be up-regulated in response to ischemic-reperfusion injury in rats and in an in vitro ischemic injury model in primary cortical neuron cultures . Furthermore, haemoglobin chains have been shown to be down-regulated in a rat model of chronic stress , in microglia, astrocytes and brain mitochondria from aged mice [42, 47]. The role of haemoglobin in neurons is unclear, although it been suggested to act as an oxygen carrier or to be protective against oxidative stress [2, 46]. We may speculate about the relationship between Smn and haemoglobin, however it is perhaps more likely that the down-regulation of Hbb transcripts that we observed is a downstream consequence of the cellular defects caused by a reduction in Smn levels. However, as outlined above, there is strong evidence that haemoglobin can influence pathways implicated in neuronal vulnerability, and this makes it an exciting candidate for further investigation.
Up-regulation of transcripts involved in programmed cell death in selectively vulnerable motor neurons that precedes NMJ pathology
Functional clustering of the transcriptional changes occurring between SMAv and WTv motor neurons revealed a number of clusters, most of which were common with the functional clustering of changes between SMAr and WTr motor neurons (Table 2). One functional cluster that was only observed in SMAv vs. WTv motor neurons, was an up-regulation of factors involved in cell death (Fig. 9). This functional cluster was not significantly altered in equivalent analysis on SMAr vs. WTr motor neurons. There were 13 transcripts that are associated with cell death pathways, which were up-regulated in SMAv motor neurons compared to SMAr (Fig. 9a). Of these, only 2 were statistically up-regulated in SMAr motor neuron compared to WTr (Fig. 9a). Of particular note, 5 of these transcripts were up-regulated in SMAv motor neurons compared to SMAr (Fig. 9a). KEGG pathway analysis revealed that a number of these up-regulated transcripts pertain to the P53 signalling pathway. Specifically, the RNAseq results indicated a 4.8 Log2 fold and 5.9 Log2 fold up-regulation of PMAIP and CDKN1a respectively. This up-regulation was confirmed by qPCR on cDNA from laser captured motor neurons. An up-regulation of Fas, PMAIP and CDKN1a was also seen in cDNA deriving from whole spinal cord from P10 Smn 2B/- mice (Fig. 9c). A similar up-regulation was observed in spinal cord from pre-symptomatic (P2) Smn −/− ;SMN2 mice, a more severe mouse model of SMA (Fig. 9c).
The up-regulation of factors that are strongly implicated in apoptotic pathways may at first glance appear unsurprising. However, we must remember that this up-regulation is occurring in two different mouse models of SMA at a pre-degenerative time point, prior to any NMJ loss. This suggests that cell death pathways are activated at the cell body prior to pathology at the NMJ. A central debate in the SMA research field has been whether SMA is due to a loss of a central housekeeping function for Smn, or whether it is due to a loss of a specific axonal or synaptic role for Smn [8, 15]. The observation that NMJs are lost so early in the disease has often been used as evidence for the latter. The data presented here suggests that cell death pathways are activated before degenerative events can be observed at the axon or NMJ. Clearly, this work does not eliminate a role for Smn in the axon or synapse. It also remains possible that defects in this synaptic or axonal role leads to cell death activation at the cell body, which is followed by the withdrawal of synaptic and axonal compartments. However, it is intriguing to draw parallels to congenic myasthenia syndromes, which are caused by specific synaptic defects . In these conditions there is profound synaptic dysfunction and denervation, however this does not lead to activation of cell death pathways and motor neuron cell body loss. Further work is clearly required to determine the time course and specific location of motor unit pathology in SMA and investigate whether synaptic defects in SMA are a cause or consequence of cell death pathway activation at the cell body.
What transcriptional changes correlate with motor neuron vulnerability: Transcriptional differences selectively occurring in vulnerable motor neurons
In this study we were particularly keen to identify regulators of motor neuron vulnerability. For this analysis, we reasoned that any changes that occurred only in SMAv motor neurons compared to either WTv motor neurons, or to SMAr motor neurons, could be specifically implicated in NMJ and motor unit pathology. We therefore looked for transcriptional changes that occurred in SMAv motor neurons compared to SMAr and compared to WTv. By comparing SMAv and SMAr motor neurons, we identified 1299 transcripts that were up or down-regulated (Fig. 4). As highlighted above, when comparing SMAv and WTv motor neurons, we identified 1277 transcripts that were up or down-regulated (Fig. 4). Of the changes that were identified in these two comparisons, 292 were common, and notably 94 % occurred with the same directional regulation (Fig. 5b). Indeed we found 140 and 134 transcriptional changes which were up or down regulated respectively in SMAv motor neurons compared to both WTv and SMAr (Additional file 3). Of particular note, both IGF1 and IGF2 were significantly down-regulated in SMAv motor neurons. IGF1 has previously been reported to be down-regulated in SMA models, and increasing IGF1 levels have been shown to be phenotypically beneficial [4, 50, 51]. IGF2 levels have been associated with differential neuronal vulnerability in an ALS model, down-regulated in an experimental model of stress, and suggested to have neuroprotective qualities against excitotoxicity [1, 22]. RNAseq results also implied that PMAIP and CDKN1a were selectively up-regulated in SMAv motor neurons. For PMAIP, this result was confirmed by qPCR (Fig. 9).
Functional clustering of these changes highlighted a number of cellular pathways, which are detailed in Table 3. Due to their enrichment selectively in vulnerable motor neurons, these clusters are all potentially of interest and worthy of further investigation. We were particularly interested to note the decrease in factors involved in the positive regulation of DNA repair. The transcripts within this functional cluster are listed in Table 4. This presents the possibility that the response to DNA damage is decreased in selectively vulnerable motor neurons. In order to investigate this further, we looked at phospho-histone H2AX (pH2AX) levels in differentially vulnerable motor neurons (Fig. 10). pH2AX is one of the first proteins to be recruited to sites of DNA damage, and recruits other proteins involved in DNA repair to form a DNA repair complex. Immunohistochemical staining of differentially vulnerable motor neurons with antibodies against pH2AX revealed a marked increase in foci number in WTr motor neurons compared to WTv. This suggests that there is a basal increase in the level of DNA repair occurring in less vulnerable motor neurons. The number of pH2AX foci was reduced in both SMAv and SMAr motor neurons compared to WTr, suggesting a further decrease in the levels of DNA repair in Smn 2B/- motor neurons. Interestingly, there remained a trend for increased DNA repair in SMAr motor neurons compared to SMAv. This data indicated that the number of DNA repair complexes correlates with relative vulnerability.
The data presented above provide preliminary evidence that an increase in DNA repair can be neuroprotective. It is unclear whether this is because there is a greater requirement for DNA repair in less vulnerable motor neurons, or whether the requirement is similar in both cell populations, but DNA repair is just more efficient in less vulnerable cells. The DNA damage repair system is a crucial system to maintain genomic integrity, which is especially relevant for terminally differentiated and long-lived cells such as neurons. It is easy to speculate how an increase in the activity in the basic cellular repair mechanisms could be neuroprotective. Indeed, DNA damage can be a key mechanism by which apoptosis is induced. The increase in DNA repair in less vulnerable motor neurons may reduce the likelihood of apoptotic activation. A number of hypotheses state that stressed cells cope better with additional stress. It is possible that this basal higher level of cellular stress primes selective motor neuron populations to cope better with the additional stresses when Smn levels are reduced. Clearly this idea requires further validation in this context, however this has important implications both for development of neuroprotective strategies, and for understanding the mechanisms by which sub-populations of neurons are rendered more vulnerable to a variety of insults.
In summary, the data presented above represent detailed transcriptional analysis on differentially vulnerable motor neurons from an SMA mouse model at a pre-symptomatic time point. They highlight a number of pathways which are disrupted upon Smn depletion, including a reduction of transcripts involved with the ribosome, rRNA binding, ubiquitination and oxidative phosphorylation. Subsequent work is required to ascertain the mechanisms by which reduced Smn levels impact upon these pathways. We have also revealed an early and selective up-regulation of cell death pathways, and it will now be key to understand how cell death pathway activation relates to the time course of pathology within the motor unit. Finally, we show that, among a number of pathways of potential interest, an increase in DNA damage repair complexes correlate with a reduced vulnerability. Future efforts dissecting the mechanism of this protection are now warranted.
Andrus BM, Blizinsky K, Vedell PT, Dennis K, Shukla PK, Schaffer DJ, Radulovic J, Churchill GA, Redei EE (2012) Gene expression patterns in the hippocampus and amygdala of endogenous depression and chronic stress models. Mol Psychiatry 17:49–61. doi:10.1038/mp.2010.119
Ascenzi P, Gustincich S, Marino M (2014) Mammalian nerve globins in search of functions. IUBMB Life 66:268–276. doi:10.1002/iub.1267
Baumer D, Lee S, Nicholson G, Davies JL, Parkinson NJ, Murray LM, Gillingwater TH, Ansorge O, Davies KE, Talbot K (2009) Alternative splicing events are a late feature of pathology in a mouse model of spinal muscular atrophy. PLoS Genet 5, e1000773. doi:10.1371/journal.pgen.1000773
Bosch-Marce M, Wee CD, Martinez TL, Lipkes CE, Choe DW, Kong L, Van Meerbeke JP, Musaro A, Sumner CJ (2011) Increased IGF-1 in muscle modulates the phenotype of severe SMA mice. Hum Mol Genet 20:1844–1853. doi:10.1093/hmg/ddr067
Bowerman M. MLM, Beauvais A., Pinheiro B., Kothary R. (2011) A critical Smn threshold in mice dictates onset of an intermediate Spinal Muscular Atrophy phenotype associated with a distinct neuromuscular junction pathology. Neuromuscul Disord. In Press
Brockington A, Ning K, Heath PR, Wood E, Kirby J, Fusi N, Lawrence N, Wharton SB, Ince PG, Shaw PJ (2013) Unravelling the enigma of selective vulnerability in neurodegeneration: motor neurons resistant to degeneration in ALS show distinct gene expression characteristics and decreased susceptibility to excitotoxicity. Acta Neuropathol 125:95–109. doi:10.1007/s00401-012-1058-5
Butterfield RJ, Stevenson TJ, Xing L, Newcomb TM, Nelson B, Zeng W, Li X, Lu HM, Lu H, Farwell Gonzalez KD, Wei JP, Chao EC, Prior TW, Snyder PJ, Bonkowsky JL, Swoboda KJ (2014) Congenital lethal motor neuron disease with a novel defect in ribosome biogenesis. Neurology 82:1322–1330. doi:10.1212/WNL.0000000000000305
Coady TH, Lorson CL (2011) SMN in spinal muscular atrophy and snRNP biogenesis. Wiley interdisciplinary reviews RNA 2:546–564. doi:10.1002/wrna.76
Corti S, Locatelli F, Papadimitriou D, Del Bo R, Nizzardo M, Nardini M, Donadoni C, Salani S, Fortunato F, Strazzer S, Bresolin N, Comi GP (2007) Neural stem cells LewisX+ CXCR4+ modify disease progression in an amyotrophic lateral sclerosis model. Brain j neurol 130:1289–1305. doi:10.1093/brain/awm043
Dachs E, Hereu M, Piedrafita L, Casanovas A, Caldero J, Esquerda JE (2011) Defective neuromuscular junction organization and postnatal myogenesis in mice with severe spinal muscular atrophy. J Neuropathol Exp Neurol 70:444–461. doi:10.1097/NEN.0b013e31821cbd8b
Dale JM, Shen H, Barry DM, Garcia VB, Rose FF Jr, Lorson CL, Garcia ML (2011) The spinal muscular atrophy mouse model, SMADelta7, displays altered axonal transport without global neurofilament alterations. Acta Neuropathol 122:331–341. doi:10.1007/s00401-011-0848-5
de Planell-Saguer M, Schroeder DG, Rodicio MC, Cox GA, Mourelatos Z (2009) Biochemical and genetic evidence for a role of IGHMBP2 in the translational machinery. Hum Mol Genet 18:2115–2126. doi:10.1093/hmg/ddp134
Dennis G Jr, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, Lempicki RA (2003) DAVID: Database for Annotation, Visualization, and Integrated Discovery. Genome Biol 4:P3
Engel AG, Shen XM, Selcen D, Sine SM (2015) Congenital myasthenic syndromes: pathogenesis, diagnosis, and treatment. Lancet Neurol 14:420–434. doi:10.1016/S1474-4422(14)70201-7
Fallini C, Bassell GJ, Rossoll W (2012) Spinal muscular atrophy: the role of SMN in axonal mRNA regulation. Brain Res 1462:81–92. doi:10.1016/j.brainres.2012.01.044
Filezac de L'Etang A, Maharjan N, Cordeiro Brana M, Ruegsegger C, Rehmann R, Goswami A, Roos A, Troost D, Schneider BL, Weis J, Saxena S (2015) Marinesco-Sjogren syndrome protein SIL1 regulates motor neuron subtype-selective ER stress in ALS. Nat Neurosci 18:227–238. doi:10.1038/nn.3903
Fritzsch B (1993) Fast axonal diffusion of 3000 molecular weight dextran amines. J Neurosci Methods 50:95–103
Groen EJ, Fumoto K, Blokhuis AM, Engelen-Lee J, Zhou Y, van den Heuvel DM, Koppers M, van Diggelen F, van Heest J, Demmers JA, Kirby J, Shaw PJ, Aronica E, Spliet WG, Veldink JH, van den Berg LH, Pasterkamp RJ (2013) ALS-associated mutations in FUS disrupt the axonal distribution and function of SMN. Hum Mol Genet 22:3690–3704. doi:10.1093/hmg/ddt222
Guenther UP, Handoko L, Laggerbauer B, Jablonka S, Chari A, Alzheimer M, Ohmer J, Plottner O, Gehring N, Sickmann A, von Au K, Schuelke M, Fischer U (2009) IGHMBP2 is a ribosome-associated helicase inactive in the neuromuscular disorder distal SMA type 1 (DSMA1). Hum Mol Genet 18:1288–1300. doi:10.1093/hmg/ddp028
Hammond SM, Gogliotti RG, Rao V, Beauvais A, Kothary R, DiDonato CJ (2010) Mouse survival motor neuron alleles that mimic SMN2 splicing and are inducible rescue embryonic lethality early in development but not late. PLoS One 5, e15887. doi:10.1371/journal.pone.0015887
He Y, Hua Y, Liu W, Hu H, Keep RF, Xi G (2009) Effects of cerebral ischemia on neuronal hemoglobin. official journal of the International Society of Cerebral Blood F. 29:596–605. doi:10.1038/jcbfm.2008.145
Hedlund E, Karlsson M, Osborn T, Ludwig W, Isacson O (2010) Global gene expression profiling of somatic motor neuron populations with different vulnerability identify molecules and pathways of degeneration and protection. Brain j neurol 133:2313–2330. doi:10.1093/brain/awq167
Hsieh-Li HM, Chang JG, Jong YJ, Wu MH, Wang NM, Tsai CH, Li H (2000) A mouse model for spinal muscular atrophy. Nat Genet 24:66–70. doi:10.1038/71709
Hua Y, Sahashi K, Rigo F, Hung G, Horev G, Bennett CF, Krainer AR (2011) Peripheral SMN restoration is essential for long-term rescue of a severe spinal muscular atrophy mouse model. Nature 478:123–126. doi:10.1038/nature10485
Jablonka S, Beck M, Lechner BD, Mayer C, Sendtner M (2007) Defective Ca2+ channel clustering in axon terminals disturbs excitability in motoneurons in spinal muscular atrophy. J Cell Biol 179:139–149. doi:10.1083/jcb.200703187
Kanning KC, Kaplan A, Henderson CE (2010) Motor neuron diversity in development and disease. Annu Rev Neurosci 33:409–440. doi:10.1146/annurev.neuro.051508.135722
Kaplan A, Spiller KJ, Towne C, Kanning KC, Choe GT, Geber A, Akay T, Aebischer P, Henderson CE (2014) Neuronal matrix metalloproteinase-9 is a determinant of selective neurodegeneration. Neuron 81:333–348. doi:10.1016/j.neuron.2013.12.009
Kariya S, Park GH, Maeno-Hikichi Y, Leykekhman O, Lutz C, Arkovitz MS, Landmesser LT, Monani UR (2008) Reduced SMN protein impairs maturation of the neuromuscular junctions in mouse models of spinal muscular atrophy. Hum Mol Genet 17:2552–2569. doi:10.1093/hmg/ddn156
Kong L, Wang X, Choe DW, Polley M, Burnett BG, Bosch-Marce M, Griffin JW, Rich MM, Sumner CJ (2009) Impaired synaptic vesicle release and immaturity of neuromuscular junctions in spinal muscular atrophy mice. J Neurosci 29:842–851. doi:10.1523/JNEUROSCI.4434-08.2009
Kwiatkowski TJ Jr, Bosco DA, Leclerc AL, Tamrazian E, Vanderburg CR, Russ C, Davis A, Gilchrist J, Kasarskis EJ, Munsat T, Valdmanis P, Rouleau GA, Hosler BA, Cortelli P, de Jong PJ, Yoshinaga Y, Haines JL, Pericak-Vance MA, Yan J, Ticozzi N, Siddique T, McKenna-Yasek D, Sapp PC, Horvitz HR, Landers JE, Brown RH Jr (2009) Mutations in the FUS/TLS gene on chromosome 16 cause familial amyotrophic lateral sclerosis. Science 323:1205–1208. doi:10.1126/science.1166066
Le TT, Pham LT, Butchbach ME, Zhang HL, Monani UR, Coovert DD, Gavrilina TO, Xing L, Bassell GJ, Burghes AH (2005) SMNDelta7, the major product of the centromeric survival motor neuron (SMN2) gene, extends survival in mice with spinal muscular atrophy and associates with full-length SMN. Hum Mol Genet 14:845–857. doi:10.1093/hmg/ddi078
Lefebvre S, Burlet P, Viollet L, Bertrandy S, Huber C, Belser C, Munnich A (2002) A novel association of the SMN protein with two major non-ribosomal nucleolar proteins and its implication in spinal muscular atrophy. Hum Mol Genet 11:1017–1027
Ling KK, Gibbs RM, Feng Z, Ko CP (2012) Severe neuromuscular denervation of clinically relevant muscles in a mouse model of spinal muscular atrophy. Hum Mol Genet 21:185–195. doi:10.1093/hmg/ddr453
Ling KK, Lin MY, Zingg B, Feng Z, Ko CP (2010) Synaptic defects in the spinal and neuromuscular circuitry in a mouse model of spinal muscular atrophy. PLoS One 5, e15457. doi:10.1371/journal.pone.0015457
Livak KJ, Schmittgen TD (2001) Analysis of relative gene expression data using real-time quantitative PCR and the 2(−Delta Delta C(T)) Method. Methods 25:402–408. doi:10.1006/meth.2001.1262
Lotti F, Imlach WL, Saieva L, Beck ES, le Hao T, Li DK, Jiao W, Mentis GZ, Beattie CE, McCabe BD, Pellizzoni L (2012) An SMN-dependent U12 splicing event essential for motor circuit function. Cell 151:440–454. doi:10.1016/j.cell.2012.09.012
Martinez-Hernandez R, Bernal S, Also-Rallo E, Alias L, Barcelo MJ, Hereu M, Esquerda JE, Tizzano EF (2013) Synaptic defects in type I spinal muscular atrophy in human development. J Pathol 229:49–61. doi:10.1002/path.4080
Monani UR, Sendtner M, Coovert DD, Parsons DW, Andreassi C, Le TT, Jablonka S, Schrank B, Rossoll W, Prior TW, Morris GE, Burghes AH (2000) The human centromeric survival motor neuron gene (SMN2) rescues embryonic lethality in Smn(−/−) mice and results in a mouse with spinal muscular atrophy. Hum Mol Genet 9:333–339.
Murray LM, Beauvais A, Bhanot K, Kothary R (2013) Defects in neuromuscular junction remodelling in the Smn(2B/-) mouse model of spinal muscular atrophy. Neurobiol Dis 49:57–67. doi:10.1016/j.nbd.2012.08.019
Murray LM, Comley LH, Thomson D, Parkinson N, Talbot K, Gillingwater TH (2008) Selective vulnerability of motor neurons and dissociation of pre- and post-synaptic pathology at the neuromuscular junction in mouse models of spinal muscular atrophy. Hum Mol Genet 17:949–962. doi:10.1093/hmg/ddm367
Ohyagi Y, Yamada T, Goto I (1994) Hemoglobin as a novel protein developmentally regulated in neurons. Brain Res 635:323–327
Orre M, Kamphuis W, Osborn LM, Melief J, Kooijman L, Huitinga I, Klooster J, Bossers K, Hol EM (2014) Acute isolation and transcriptome characterization of cortical astrocytes and microglia from young and aged mice. Neurobiol Aging 35:1–14. doi:10.1016/j.neurobiolaging.2013.07.008
Richter F, Meurers BH, Zhu C, Medvedeva VP, Chesselet MF (2009) Neurons express hemoglobin alpha- and beta-chains in rat and human brains. J Comp Neurol 515:538–547. doi:10.1002/cne.22062
Ruiz R, Casanas JJ, Torres-Benito L, Cano R, Tabares L (2010) Altered intracellular Ca2+ homeostasis in nerve terminals of severe spinal muscular atrophy mice. J Neurosci 30:849–857. doi:10.1523/JNEUROSCI.4496-09.2010
Sanchez G, Dury AY, Murray LM, Biondi O, Tadesse H, El Fatimy R, Kothary R, Charbonnier F, Khandjian EW, Cote J (2013) A novel function for the survival motoneuron protein as a translational regulator. Hum Mol Genet 22:668–684. doi:10.1093/hmg/dds474
Schelshorn DW, Schneider A, Kuschinsky W, Weber D, Kruger C, Dittgen T, Burgers HF, Sabouri F, Gassler N, Bach A, Maurer MH (2009) Expression of hemoglobin in rodent neurons. J cereb blood flow metab : official journal of the International Society of Cerebral Blood Flow and Metabolism 29:585–595. doi:10.1038/jcbfm.2008.152
Shephard F, Greville-Heygate O, Marsh O, Anderson S, Chakrabarti L (2014) A mitochondrial location for haemoglobins--dynamic distribution in ageing and Parkinson's disease. Mitochondrion 14:64–72. doi:10.1016/j.mito.2013.12.001
Sleigh JN, Gillingwater TH, Talbot K (2011) The contribution of mouse models to understanding the pathogenesis of spinal muscular atrophy. Dis Model Mech 4:457–467. doi:10.1242/dmm.007245
Thomson SR, Nahon JE, Mutsaers CA, Thomson D, Hamilton G, Parson SH, Gillingwater TH (2012) Morphological characteristics of motor neurons do not determine their relative susceptibility to degeneration in a mouse model of severe spinal muscular atrophy. PLoS One 7, e52605. doi:10.1371/journal.pone.0052605
Tsai LK, Chen CL, Ting CH, Lin-Chao S, Hwu WL, Dodge JC, Passini MA, Cheng SH (2014) Systemic administration of a recombinant AAV1 vector encoding IGF-1 improves disease manifestations in SMA mice. Mol ther: the journal of the American Society of Gene Therapy 22:1450–1459. doi:10.1038/mt.2014.84
Tsai LK, Chen YC, Cheng WC, Ting CH, Dodge JC, Hwu WL, Cheng SH, Passini MA (2012) IGF-1 delivery to CNS attenuates motor neuron cell death but does not improve motor function in type III SMA mice. Neurobiol Dis 45:272–279. doi:10.1016/j.nbd.2011.06.021
Vance C, Rogelj B, Hortobagyi T, De Vos KJ, Nishimura AL, Sreedharan J, Hu X, Smith B, Ruddy D, Wright P, Ganesalingam J, Williams KL, Tripathi V, Al-Saraj S, Al-Chalabi A, Leigh PN, Blair IP, Nicholson G, de Belleroche J, Gallo JM, Miller CC, Shaw CE (2009) Mutations in FUS, an RNA processing protein, cause familial amyotrophic lateral sclerosis type 6. Science 323:1208–1211. doi:10.1126/science.1165942
Wadman RI, Vrancken AF, van den Berg LH, van der Pol WL (2012) Dysfunction of the neuromuscular junction in spinal muscular atrophy types 2 and 3. Neurology 79:2050–2055. doi:10.1212/WNL.0b013e3182749eca
Wehner KA, Ayala L, Kim Y, Young PJ, Hosler BA, Lorson CL, Baserga SJ, Francis JW (2002) Survival motor neuron protein in the nucleolus of mammalian neurons. Brain Res 945:160–173
Wishart TM, Mutsaers CA, Riessland M, Reimer MM, Hunter G, Hannam ML, Eaton SL, Fuller HR, Roche SL, Somers E, Morse R, Young PJ, Lamont DJ, Hammerschmidt M, Joshi A, Hohenstein P, Morris GE, Parson SH, Skehel PA, Becker T, Robinson IM, Becker CG, Wirth B, Gillingwater TH (2014) Dysregulation of ubiquitin homeostasis and beta-catenin signaling promote spinal muscular atrophy. J Clin Invest 124:1821–1834. doi:10.1172/JCI71318
Wootz H, Fitzsimons-Kantamneni E, Larhammar M, Rotterman TM, Enjin A, Patra K, Andre E, Van Zundert B, Kullander K, Alvarez FJ (2013) Alterations in the motor neuron-renshaw cell circuit in the Sod1(G93A) mouse model. J Comp Neurol 521:1449–1469. doi:10.1002/cne.23266
Workman E, Kolb SJ, Battle DJ (2012) Spliceosomal small nuclear ribonucleoprotein biogenesis defects and motor neuron selectivity in spinal muscular atrophy. Brain Res 1462:93–99. doi:10.1016/j.brainres.2012.02.051
Wu CY, Whye D, Glazewski L, Choe L, Kerr D, Lee KH, Mason RW, Wang W (2011) Proteomic assessment of a cell model of spinal muscular atrophy. BMC Neurosci 12:25. doi:10.1186/1471-2202-12-25
Yamazaki T, Chen S, Yu Y, Yan B, Haertlein TC, Carrasco MA, Tapia JC, Zhai B, Das R, Lalancette-Hebert M, Sharma A, Chandran S, Sullivan G, Nishimura AL, Shaw CE, Gygi SP, Shneider NA, Maniatis T, Reed R (2012) FUS-SMN protein interactions link the motor neuron diseases ALS and SMA. Cell rep 2:799–806. doi:10.1016/j.celrep.2012.08.025
Zhang Z, Pinto AM, Wan L, Wang W, Berg MG, Oliva I, Singh LN, Dengler C, Wei Z, Dreyfuss G (2013) Dysregulation of synaptogenesis genes antecedes motor neuron pathology in spinal muscular atrophy. Proc Natl Acad Sci U S A 110:19348–19353. doi:10.1073/pnas.1319280110
We would like to thank Caroline Vergette, Christopher Porter, Ted Perkins and Sophia Rahimi for assistance with experiments and data analysis. We would like to thank Dr. Erin Basset for constructive suggestions and all members of the Kothary laboratory for helpful discussions.
This work was supported by grants from Cure SMA Canada (Formerly Families of SMA Canada, KT1112 to LM and RK) and Cure SMA (Formerly Families of SMA, grant number MU1415 to LM). LM is the recipient of a Muscular Dystrophy Association Development Grant (grant number 294433) and an Emerging Investigator award from Fight SMA and the Gwendolyn Strong Foundation. RK is supported by grants from the Canadian Institutes of Health Research (grant number MOP–130279) and the Muscular Dystrophy Association (grant number 294568).
All applicable national and institutional guidelines for the care and use of animals were followed.
Compliance with ethical standards
Research involving animals: All animal procedures have been carried out in accordance with the guidelines set out by the Animal Care and Veterinary Committee at the University of Ottawa and the UK home office, where appropriate.
The authors declare that they have no conflict of interest. The funders (identified above) had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
LM conceived of the study, and directed it, performed the NMJ analysis, gather material for LCM and performed qPCR validation. AB and SG developed techniques for LCM, participated in the collection of motor neuron and performed qPCR validation. NC performed qPCR validation and performed the statistical analysis. RK conceived of the study, participated in its design and coordination and helped to draft the manuscript. All authors read and approved the final manuscript.
Table shows a summary of the sequencing data obtained from the RNAseq analysis, as analysed using Tophat software. Table show the number of reads (reads), the percentage which mapped to a unique location (% unique), the percentage which mapped to a distinct location (% distinct), the number of mapped location (mapped locations), the number of mapped reads (mapped reads) and the percentage of 10.1186/s40478-015-0231-1 reads which were mapped (% mapped) for each of the 8 samples. (DOCX 53 kb)
Table show a list of genes which are differentially expressed between SMAr and SMAv motor neurons when compared to their respective wild-types. (DOCX 113 kb)
Table show a list of genes which are differentially expressed between SMAv motor neurons when compared to SMAr or Wtv motor neurons. (DOCX 113 kb)
About this article
Cite this article
Murray, L.M., Beauvais, A., Gibeault, S. et al. Transcriptional profiling of differentially vulnerable motor neurons at pre-symptomatic stage in the Smn 2b/- mouse model of spinal muscular atrophy. acta neuropathol commun 3, 55 (2015). https://doi.org/10.1186/s40478-015-0231-1