Deregulated expression of cytoskeleton related genes in the spinal cord and sciatic nerve of presymptomatic SOD1G93A Amyotrophic Lateral Sclerosis mouse model

Early molecular events related to cytoskeleton are poorly described in Amyotrophic Lateral Sclerosis (ALS), especially in the Schwann cell (SC), which offers strong trophic support to motor neurons. Database for Annotation, Visualization and Integrated Discovery (DAVID) tool identified cytoskeleton-related genes by employing the Cellular Component Ontology (CCO) in a large gene profiling of lumbar spinal cord and sciatic nerve of presymptomatic SOD1G93A mice. One and five CCO terms related to cytoskeleton were described from the spinal cord deregulated genes of 40 days (actin cytoskeleton) and 80 days (microtubule cytoskeleton, cytoskeleton part, actin cytoskeleton, neurofilament cytoskeleton, and cytoskeleton) old transgene mice, respectively. Also, four terms were depicted from the deregulated genes of sciatic nerve of 60 days old transgenes (actin cytoskeleton, cytoskeleton part, microtubule cytoskeleton and cytoskeleton). Kif1b was the unique deregulated gene in more than one studied region or presymptomatic age. The expression of Kif1b [quantitative polymerase chain reaction (qPCR)] elevated in the lumbar spinal cord (40 days old) and decreased in the sciatic nerve (60 days old) of presymptomatic ALS mice, results that were in line to microarray findings. Upregulation (24.8 fold) of Kif1b was seen in laser microdissected enriched immunolabeled motor neurons from the spinal cord of 40 days old presymptomatic SOD1G93A mice. Furthermore, Kif1b was dowregulated in the sciatic nerve Schwann cells of presymptomatic ALS mice (60 days old) that were enriched by means of cell microdissection (6.35 fold), cell sorting (3.53 fold), and primary culture (2.70 fold) technologies. The gene regulation of cytoskeleton molecules is an important occurrence in motor neurons and Schwann cells in presymptomatic stages of ALS and may be relevant in the dying back mechanisms of neuronal death. Furthermore, a differential regulation of Kif1b in the spinal cord and sciatic nerve cells emerged as key event in ALS.


INTRODUCTION
Amyotrophic Lateral Sclerosis (ALS) is a progressive, rapid and fatal neurodegenerative disease that affects motor neurons of the spinal cord, brainstem, and cerebral cortex (Tripathi and Al-Chalabi, 2008). The mortality is often due to a respiratory failure (Shaw et al., 2001).
ALS pathogenesis is still unknown. Nevertheless, the mechanisms underlying neurodegeneration in ALS seem multifactorial and take place in neurons and non-neuronal cells (Boillee et al., 2006a,b;Yamanaka et al., 2008;Wang et al., 2013). Recent analyses have showed the involvement of cytoskeleton, leading a disruption of intracellular function, and intercellular communication, with relevance to the triggering of motor neuron death (Guipponi et al., 2010). In fact, those events are especially important to motor neurons, highly polarized cells that establish contact with their target and surrounding Schwann cells through long axons.
The steady bidirectional flux of molecules and organelles in the motor neuron axons is necessary for cell survival and maintenance (Liu et al., 2013;Vinsant et al., 2013). In this context, cytoskeleton impairments might account for the described ALS mechanisms as regarding axonal/mitochondrial alteration, signaling endosome dysfunction, protein aggregation and apoptosis (Boillee et al., 2006a,b;Ferraiuolo et al., 2011;Kiernan et al., 2011;Usuki et al., 2012).
The presence of Schwann cell-expressing the distress biomarker ATF-3 in spinal nerves (Malaspina et al., 2010) before symptom onset suggests the contribution of those cells to ALS pathogenesis (Keller et al., 2009). Remarkably, axonal retraction and motor neuron disconnection from neuromuscular joints are ALS early events (Fischer et al., 2004;Parkhouse et al., 2008) that seem to be induced by Schwann cell mechanisms (Vinsant et al., 2013). For instance, distal Schwann cells produce semaphorin 3, a chemorepellent molecule for terminal axons (De Winter et al., 2006). Furthermore, the expression of the glial intermediate filament protein GFAP in Schwann cells of the peripheral nerve implies a dynamic alteration of cytoskeleton and turnover of myelin sheath (Hanyu et al., 1982). Moreover, an accumulation of iNOS immunoreactivity at the paranodal regions of Schwann cell myelin sheaths of peripheral nerves of presymptomatic ALS mice gives additional evidence for the impaired paracrine mechanisms between motor neuron and Schwann cell (Chen et al., 2010). Thus, it should be considered that the early peripheral events related to cytoskeleton of motor neurons and Schwann cells may contribute to neuronal dying back via disruption of peripheral neurotrophic stimuli (Keller et al., 2009;Dadon-Nachum et al., 2011;Gould and Kendall, 2011;Gould and Oppenheim, 2011;Liu et al., 2013).
As a short lasting disease, the challenge on ALS investigation is the employment of an adequate experimental model to evaluate presymptomatic mechanisms triggering motor neuron death. With this regard, it is known that transgene mice expressing human mutant copper/zinc superoxide dismutase 1 (SOD1 G93A ) develop clinical and pathological features similar to those seen in human ALS and are considered an excellent model to study the pathogenic mechanisms of the disease (Gama Sosa et al., 2012). The model is particularly useful to evaluate the events related to motor neuron degeneration prior neurological symptoms (Alves et al., 2011).
Large-scale microarray-based gene expression has been trying to identify new molecular cues potentially involved in the ALS pathogenesis both in animal models and postmortem tissue (Olsen et al., 2001;Hensley et al., 2002;Yoshihara et al., 2002;Dangond et al., 2004;Perrin et al., 2005;Ferraiuolo et al., 2007Ferraiuolo et al., , 2009Fukada et al., 2007;Lobsiger et al., 2007;Vargas et al., 2008;Kudo et al., 2010;Boutahar et al., 2011;Cooper-Knock et al., 2012). However, there is a lack of investigation on the analysis of cytoskeleton-related gene profiling. The evaluation of deregulated genes in specific enriched cells obtained by in vitro purification, single cell laser microdissection or cell sorting might contribute to refine the alterations of gene expression-related to cytoskeleton molecules on specific cells of peripheral motor neuron unit.
By means of a high-density oligonucleotide microarray-linked to specific tools capable to identify cellular components, the aim of this work was to identify the regulation of cytoskeleton-related genes in the presymptomatic stage in the spinal cord and sciatic nerve of the SOD1 G93A mouse model. The work has also evaluated the modulation of Kif1b in the enriched spinal cord motor neurons and sciatic nerve Schwann cells.

ANIMAL AND TISSUE SAMPLE
Transgene SOD1 G93A mice (The Jackson Laboratory, Bar Harbor, ME, USA) were crossbred and the colony was maintained in a specific pathogen-free environment of the animal facility of University of São Paulo Medical School (São Paulo, Brazil) as described previously (Gurney, 1994;Scorisa et al., 2010;Alves et al., 2011). Animals were kept under controlled temperature and humidity conditions with a standardized light-dark cycle (lights on at 7.00 a.m. and off at 7.00 p.m.) and free access to food pellets and tap water. Mice were genotyped by PCR amplification of DNA extracted from their tails in order to identify the SOD1 mutation (Gurney, 1994;Scorisa et al., 2010;Alves et al., 2011). The study was conducted under protocols approved by the Animal Care and Use of Ethic Committee at University of São Paulo and in accordance to the Guide for the Care and Use of Laboratory Animals adopted by the National Institutes of Health. Forty, 60, and 80 days old presymptomatic specific pathogenfree male SOD1 G93A mice and their age-paired wild-type controls (20-25 g body weight) were used in the experiments. No motor neuron death was seen in those animal ages (Alves et al., 2011) so that they were chosen for the present presymptomatic analyses. Animals were killed by decapitation. Lumbar spinal cords (40 and 80 days old mice) and sciatic nerves (60 days old mice) were removed, frozen, and stored at −80 • C until use. Four-five animals per group were used in the microarray experiments. The quantitative polymerase chain reaction (qPCR) analyses of lumbar spinal cords (40 days old mice), and sciatic nerves (60 days old mice) as well as of enriched cells samples (60 days old mice) employed four mice of each transgene and wild-type groups.

RNA ISOLATION AND MICROARRAY EXPERIMENTS
The procedures of microarray experiments and statistical analysis of the mouse spinal cords were described in our previous publication which has employed a Whole Mouse Genome Oligo 4 × 44 K microarray platform (Agilent Technologies, USA) (De Oliveira et al., 2013). Regarding the sciatic nerve samples, total RNA was isolated using the Miniprep kit (Zymo, USA). The procedure was performed according to the manufacturer's instructions. The quantity and integrity of RNA were determined by spectrophotometer (Nanodrop, Thermo Scientific, USA) and microfluidics-based electrophoresis (Agilent 2100 Bioanalyzer, Agilent Technologies, USA), respectively. RNA samples with OD 260/280 of approximately 2.0 and RIN >7.0 were used for microarray experiments and qPCR. A pool of RNAs from neonatal organs (heart, kidney, liver) was used as reference sample. A representative electropherogram from Bioanalyzer evaluation of RNA integrity of the sciatic nerve samples is shown in the supplementary material ( Figure S1).
In the case of sciatic nerve analysis, RNAs of samples (25 ng) and reference (100 ng) were reverse transcribed by the Low-input RNA Linear Amplification kit and then transcribed to Cy3labelled (samples) or Cy5-labelled (reference) according to the manufacturer's instructions (Agilent Technologies, USA) and to previous descriptions (De Oliveira et al., 2013Oliveira et al., , 2014. A total of 300 ng of Cy3-labelled cRNA was hybridized together with the same amount of Cy5-labelled reference to Whole Mouse Genome Oligo 8 × 60 K. After an overnight hybridization at 65 • C, the slides were washed and treated with a Stabilizing and Drying Solution (Agilent Technologies, USA) and scanned (Agilent Microarray Scanner). All steps were performed according to the manufacturer's instructions (Agilent Technologies, USA).
The raw data from hybridizations and experimental conditions are available on the Gene Expression Omnibus website under accession numbers GSE50642 (spinal cord analysis, according to De Oliveira et al., 2013) and GSE56926 (sciatic nerve analysis).

MICROARRAY ANALYSIS
cord in our previous study (De Oliveira et al., 2013), sciatic nerve microarray raw data (.txt files) were transferred to R v. 3.0.1 software (Team RDC, 2012) and analyzed with the Bioconductor (Gentleman et al., 2004) package limma (Smyth, 2005). Finally, the probes were tested for differential expression using a linear model followed by Bayes moderated t-test (Smyth, 2005) for the comparisons of interest. P-values < 0.05 were accepted as differentially expressed genes.

COMPLEMENTARY DNA MICROARRAY DATA ANALYSIS
The Database for Annotation, Visualization and Integrated Discovery (DAVID) v6.7b functional tool (http://david.abcc. ncifcrf.gov/) (Huang Da et al., 2008) was used to identify genes related to cytoskeleton through the Gene Ontology (GO) annotation database. DAVID analysis focused in the Cellular Component Ontology (CCO). The analysis was conducted on the lists containing the up-regulated and down-regulated genes for each experimental group. High stringency (EASE score set to 0.05) parameters were selected to improve confidence on the terms to be pointed as enriched. Cellular component terms related to cytoskeleton gene lists were then organized. The BioVenn tool (http://www.cmbi.ru.nl/cdd/biovenn/) (Hulsen et al., 2008) was used to identify common and exclusively expressed genes between groups.

LASER CAPTURE MICRODISSECTION OF MOTOR NEURONS FROM SPINAL CORD
Immunolabeled motor neurons of lumbar mouse spinal cord (SOD1 G93A and wild-type groups) were microdissected as described previously (De Oliveira et al., 2009. Spinal cord sections were rinsed for 3 min in phosphate buffered saline (PBS) containing 3% Triton X-100 (Sigma, USA) and then incubated overnight with a polyclonal goat anti-choline acetyltransferase (ChAT, 1:100; Abcam, USA) diluted in 0.3% Triton X-100 containing 1% bovine serum albumin (BSA; Sigma, USA), 1 mM dithiothreitol (DTT; Invitrogen, CAN), and 0.1 U/μl RNAse inhibitor (Invitrogen, CAN). Sections were then washed in PBS (3 × 15 s) and then incubated for 1 h in the dark and at room temperature with an Alexa 594-conjugated donkey anti-goat antibody (Invitrogen, USA) diluted (1:100) in the solution described above. Sections were rinsed carefully three times with PBS for 15 s and immediately submitted to single cell laser microdissection procedures. The ChAT immunofluorescence profiles for specific identification of motor neurons in the microdissection procedure are illustrated in Figures 1A,B.
About 100 motor neurons were isolated from each lumbar spinal cord using P.A.L.M. Microlaser Technologies (Zeiss). RNA was extracted from the microdissected motor neurons using the PicoPure RNA isolation kit (Arcturus, USA). Linear amplification of RNA was performed following Eberwine's procedure (Van Gelder et al., 1990) using the RiboampHSplus kit (Arcturus, USA) according to the manufacturer's protocol. The quantity and quality of the amplified RNA were analyzed as described above. Laser microdissected motor neuron samples were submitted to PCRs for verification of sample enrichment and the results are shown in the supplementary material ( Figure S2).

FIGURE 1 | Photomicrographs illustrating motor neuron (A,B) and
Schwann cell (C,D) laser microdissection process. The quick ChAT and S100β immunofluorescence procedures allow recognizing the motor neuron (A,B) and Schwann cell (C,D) profiles (arrows). Profiles were then selected for microdissection (A,C). After laser firing and microdissection, selected cell profiles (arrows) can no longer be visualized in the tissue (B,D). Scale bars: 20 μm.

LASER CAPTURE MICRODISSECTION OF SCHWANN CELLS FROM SCIATIC NERVE
Sciatic nerve of 60 days old mice (SOD1 G93A and wild-type groups) were rapidly removed and frozen in ice cold isopentane at −45 • C and stored at −80 • C until use. The labeling procedure was performed according to adaptation of a previous description (De Oliveira et al., 2009Oliveira et al., , 2013. Taw-mounted mouse sciatic nerve sections (5 μm) were rapidly defrosted for 30 s and fixed with ice-cold acetone, for 3 min. Sections were then rinsed (3 min) in PBS containing 3% Triton X-100 and incubated with a polyclonal rabbit anti-S100β antibody (1:200; Dako, USA) diluted in PBS containing 0.3% Triton X-100 (Sigma, USA) and 1% BSA (Sigma, USA) for 5 min. Sections were washed in PBS (3 × 15 s) and then incubated (5 min) in the dark and at room temperature with a texas red-conjugated goat anti-rabbit antibody (Jackson ImmunoResearch Laboratories, USA) diluted (1:50) in the solution described above. Sections were then rinsed carefully three times with PBS for 15 s and immediately submitted to single cell laser microdissection procedures. Schwann cell S100β immunofluorescence profiles identified in the microdissection procedure are illustrated in Figures 1C,D.
About 200 Schwann cells were isolated from each sciatic nerve using the P.A.L.M. Microlaser Technologies. The RNA was extracted from the cells and amplified as described above.
The quantity and quality of amplified RNA were analyzed as described above. Laser microdissected Schwann cell profiles were also submitted to PCRs for verification of sample enrichment. The detailed protocol and results are shown in the supplementary material ( Figure S3).

FLOW CYTOMETRY SORTING SCHWANN CELLS
Schwann cells were isolated by means of flow cytometry sorting from the sciatic nerve explants of 60 days old SOD1 G93A mice and their age-paired wild-type controls. Briefly, animals were deeply anaesthetized with sodium pentobarbital 3% (100 mg/kg, ip) and their sciatic nerves were dissected under aseptic conditions. Nerves were then placed in 60 mm dishes containing Leibovitz-15 medium (Gibco, USA), divested of their epineurial sheaths and chopped into 1 mm pieces. The fragments were then transferred to new 60 mm dishes containing D-10 culture medium [composed by DMEM (Gibco, USA) supplemented with 10% fetal bovine serum and 1% penicillin/streptomycin (Sigma, USA)] and were maintained there in 5% CO 2 at 37 • C for 5 days. The sciatic nerve fragments were then transferred to 30 mm dishes containing 2.5 ml Hanks' Balanced Salt solution (Sigma, USA), 0.05% trypsin (Gibco, USA), and 1 mg collagenase (Worthington, USA). The fragments were kept in that solution for 2 h in 5% CO 2 at 37 • C. Tissue fragments were washed with D-10 and dissociated by trituration through a 200 μl-pipette and a 19-gauge sterile needle. The suspension was centrifuged at 1500 rpm for 5 min at 4 • C and the cells were resuspended in D-10 medium. This step was repeated and cells were passed through a 70 μm-cell strainer (BD Bioscience, USA).
The cells were centrifuged at 1,500 rpm for 5 min at 4 • C and the pellets were resuspended in PBS containing 10% fetal bovine serum and 0.01% sodium azide, (Sigma, USA). Sciatic nerve-derived cell suspension was incubated with a fluorescein isothiocyanate (FITC)-conjugated mouse p75NGF receptor antibody (Abcam, USA) diluted in the buffer solution (1/200) for 1 h at room temperature as mentioned above. The p75NGF receptor labeling was employed in the cell sorting experiments because it is a well-characterized surface marker for Schwann cells (Niapour et al., 2010). The samples were then centrifuged (300 × g for 5 min at 4 • C). The pellets were washed two times and resuspended in the PBS described above (500 μl). Cells were then analyzed for type and specificity as well as separated on a FACSAria III Cell Sorter (BD Biosciences, USA). A maximum of 10 6 cells were resuspended in 500 μl of buffer. Flow cytometry dot plot Schwann cell profiles are shown in Figure 2. Details of flow cytometry procedures for cell specificity are described in the supplementary material ( Figure S4).
RNA of enriched cells was extracted using Trizol (Life Technologies, USA) according to the manufacturer's protocol. The quantity (NanoDrop 1000 Spectrophotometer) and quality (Agilent 2100 bioanalyser, RNA 6000 Pico LabChip) of RNA were analyzed as described above. Also, the Schwann cell samples were submitted to PCRs in order to access contamination from other cell types. Protocol and results regarding specificity of separated Schwann cell samples are presented in the supplementary material (Figures S3).

PRIMARY SCHWANN CELL CULTURE
Highly purified Schwann cell cultures were obtained from sciatic nerve explants taken from 60 days old SOD1 G93A and wild-type mice as described above. Nerve pieces were transferred weekly to new 60 mm dishes filled with 1 ml of D-10 for 5 weeks. Dishes were replaced every other day with a fresh medium (Oudega et al., 1997). After that period, explants were replated onto 35 mm dishes containing a solution of 1.25 U/ml dispase (Boehringer Mannheim, Germany), 0.05% collagenase (Worthington, USA), and D-10, and were kept under overnight incubation in 5% CO2 at 37 • C. Following, explants were washed in D-10 and dissociated. The resulting cells were treated with a Thy1.2 antibody (BD Bioscience, USA) and a rabbit serum complement (Calbiochem, USA) for 30 min at room temperature for fibroblast elimination. The protocol for cell enrichment was described elsewhere (Brockes et al., 1979;White et al., 1983;Dong et al., 1999) and was modified according to our experience. The obtained Schwann cells were then seeded onto laminin (Sigma, USA) coated 100 mm dishes for expansion. Twenty-four hours later, the culture medium was replaced by a D-10 medium supplemented with 2 mM forskolin (Sigma, USA) and 20 mg/ml pituitary extract (Gibco, USA). Cells were allowed to expand in that medium until confluence has reached. The medium was changed every other day in the expansion period. The cells of the third passage were used for experiments. Samples of the primary Schwann cell cultures were fixed and immunostained with S100β antibody and nuclei were stained with diamidino-2-phenylindole (DAPI) for cell type verification, as showed in Figure 3.
Schwann cell total RNA was extracted using Trizol (Life Technologies) according to the manufacturer's protocol. The quantity (NanoDrop 1000 Spectrophotometer) and quality (Agilent 2100 bioanalyser, RNA 6000 Pico LabChip) of RNA were analyzed as described above. Cultured Schwann cell RNA samples were submitted to PCRs in order to access fibroblast contamination; the protocol and results are shown in the supplementary material ( Figure S3).

QUANTITATIVE PCR
Microarray analyses identified the differentially expressed Kif1b in the spinal cord (40 days old mice) and sciatic nerve (60 days old mice) and it was the only gene that its product has been described in the context of ALS (Conforti et al., 2003;Pantelidou et al., 2007). The gene was then selected for verification by qPCR in the whole spinal cord (40 days old mice) and sciatic nerve (60 days) as well as in the enriched motor neurons and Schwann cells. The qPCR verification was performed on independent samples.
Spinal cord cDNA was synthesized from DNAse-treated 1 μg total RNA by employing a TaqMan reverse transcription kit (Applied Biosystems, USA). Sciatic nerve cDNA was synthesized from 100 ng of total RNA by using a Maxima First Strand cDNA Synthesis Kit (Thermo Scientific, USA) according to manufacturers.
The Kif1b expression was also evaluated in the enriched microdissected Schwann cells and motor neurons as well as in the Schwann cells enriched by means of cell sorting and primary culture procedures. cDNA samples of microdissected cells and of cultured Schwann cells were synthesized from 1 μg of amplified RNA as described previously (De Oliveira et al., 2013). The cDNA of flow cytometry sorting Schwann cells was synthesized from 100 ng of total RNA.
The cycling for SYBR reactions was composed by an initial denaturation at 95 • C for 10 min. Templates were amplified by 40 cycles of 95 • C for 15 s and of 60 • C for 30 s. A dissociation curve was then generated to ensure amplification of a single product and absence of primer dimers. A standard curve was generated for each primer pair in order to determine the efficiency of the PCR reaction over a range of template concentrations from 0.032 to 20 ng/μl, using cDNA synthesized from reference mouse RNA. The efficiency for each set of primers was 100 ± 5%. Gene expressions, which were normalized by Gapdh, could be determined using the Ct mathematical model (ABI PRISM 7700 Sequence Detection System protocol; Applied Biosystems). Gapdh was chosen as a housekeeping gene to normalize the qPCR values because the microarray analysis showed no alteration in the gene expression across samples.

STATISTICAL ANALYSIS
The statistical method employed in the microarray analysis is described above (Hulsen et al., 2008). Furthermore, one-tailed unpaired t-test was used to determine the statistical significance of differences in gene expression [Graphpad Prism 5 (San Diego, CA)] in the qPCR analyses.

MICROARRAY ANALYSIS
The DAVID analysis of differentially expressed genes of 40 days old SOD1 G93A mice pointed 34 enriched GO terms under high stringency conditions. The CCO indicated only one GO term and immunofluorescence (B,C) microscopy. Cultured Schwann cells were evidenced by means of S100 immunofluorescence (greenish color), and the cell nuclei (bluish color) were stained by DAPI (B,C). The different morphology of Schwann cells (arrowhead) and fibroblasts (arrow) is observed (A). S100 positive immature Schwann cells (arrow) and DAPI positive nuclei of cells lacking cytoplasmic S100 labeling (arrowhead) are seen in a 24 h plating non-purified culture (B). Vast majority of S100 immunolabeled Schwann cells possess a homogeneous morphology in a 7-days purified culture after fibroblast elimination [arrow, (C)]. Scale bars: 50 μm. related to cytoskeleton, the actin cytoskeleton, with 14 genes (six down and eight upregulated), which are shown in Table 1. DAVID also pointed 63 enriched terms from the differentially expressed genes of spinal cord of 80 days old SOD1 G93A mice ( Table 2). The CCO indicated five GO terms related to cytoskeleton (Table 2) in the spinal cord of 80 days old SOD1 G93A mice, specifically the microtubule cytoskeleton (35 genes), cytoskeleton part (53 genes), actin cytoskeleton (16 genes), neurofilament cytoskeleton (three genes), and cytoskeleton (76 genes). Those genes are overlapped with the 76 deregulated genes of cytoskeleton category (25 down and 51 upregulated).
The DAVID analysis also pointed 55 enriched terms from the deregulated genes of sciatic nerve of 60 days old SOD1 G93A mice. Furthermore, the CCO indicated four GO terms related with cytoskeleton, specifically the actin cytoskeleton (43 genes), cytoskeleton part (101 genes), microtubule cytoskeleton (64 genes), and cytoskeleton (146 genes). The 146 genes of the cytoskeleton GO term (74 down and 72 upregulated) are overlapped with all other GO terms ( Table 3).
From the above lists of differentially expressed genes related to cytoskeleton pointed by CCO,nine genes (Dync1li2,Kif1b,Kif1c,Lcp1,Lima1,Myh11,Myo1c,Srcin,Twf1) were deregulated in both spinal cord and sciatic nerve of SOD1 G93A mice at the presymptomatic ages 40 and 60 days, respectively, and 10 genes (Abl1,Calm3,Clic5,Dynll2,Krt7,Myl4,Myo1d,Ncs1,Pdlim2,Strbp) were deregulated in the above regions of 80 and 60 days mice, respectively, as shown by the Venn diagram (Figure 4). Furthermore, only the Maea was seen differentially expressed in the spinal cord of 40 and 80 days old presymptomatic SOD1 G93A mice and no gene appeared to repeat in the three lists of the studied regions of the presymptomatic mice (Figure 4). It should be pointed out that only the Kif1b from the above described differentially expressed genes related to cytoskeleton has been already mentioned in the context of ALS and has been detailed in the present cellular and molecular analyses.

MOTOR NEURON AND SCHWANN CELL ENRICHMENT
In order to analyze the modulation of gene expression in specific cell populations possibly involved in the pathogenic mechanisms in ALS, spinal cord motor neurons (lumbar regions) were obtained by means of single cell laser microdissection and the sciatic nerve Schwann cells were achieved by means of laser microdissection, flow cytometry cell sorting and cell culture. The levels of cDNA specific cell type marker that demonstrated the enrichment for each cell type obtained by respective technique are shown in the Supplementary Material (S2, S3).

Kif1b REGULATION EVIDENCED BY qPCR
qPCR analyses of Kif1b expression showed an upregulation (1.21 fold) of the gene in spinal cord of presymptomatic 40 days old SOD1 G93A mice ( Figure 5A) and a downregulation (1.57 fold) in the sciatic nerve of presymptomatic 60 days old transgene mice ( Figure 5B). These regulations were coincident and supported the microarray findings. Additionally, qPCR analyses also demonstrated the differentially expression of Kif1b in enriched cell assays using two cycle amplified RNA. Upregulation of Kif1b (24.8 fold) was seen in laser microdissected motor neurons from 40 days old SOD1 G93A mice ( Figure 5C); this regulation was in the same direction to that found in the whole spinal cord preparation of 40 days transgene mice by means of microarray and qPCR analyses. Remarkably, Kif1b was downregulated in the enriched sciatic nerve Schwann cells (60 days old SOD1 G93A mice) by means of single cell laser microdissection (6.35 fold), cell sorting (3.53 fold), and cell culture (2.70 fold) (Figure 5D), regulations that were in the same direction to that found in the whole sciatic nerve preparation of 60 days transgene mice by means of microarray and qPCR analyses.

CYTOSKELETON ENRICHMENT ANALYSIS OF GENE PROFILING IN PRESYMPTOMATIC ALS
Gene regulation of cytoskeleton-related molecules employing microarray technology has been described in the spinal cord and/or microdissected survival neurons from post mortem material of ALS patients (Jiang et al., 2005;Offen et al., 2009;Cox et al., 2010;Tanaka et al., 2012), thus reflecting the cytoskeleton responses to injury instead its role on neurodegenerative triggering. Sporadic and unsystematic results on deregulated genes-related to cytoskeleton in ALS animal models have been demonstrated in several stages of the disease (Perrin et al., 2005;Tanaka et al., 2006;Ferraiuolo et al., 2007;Guipponi et al., 2010;Kudo et al., 2010). The present study extends previous descriptions by detailing gene profiling in expanded categories of cytoskeleton-related genes in the presymptomatic ages of the ALS mouse model. The DAVID analysis was applied by using the full lists of deregulated genes of the spinal cord of 40 and 80 days old presymptomatic ALS mice we have published recently (De Oliveira et al., 2013). Our study has also demonstrated for the first time the gene profiling of cytoskeleton category in the peripheral nerve (sciatic nerve) of the ALS model in a presymptomatic period of the disease. It should be highlighted the large number of deregulated genes of the cytoskeleton-related category in the spinal cord of presymptomatic ALS mouse of the present study compared to previous descriptions that employed different methodologies (Ferraiuolo et al., 2007;Guipponi et al., 2010;Kudo et al., 2010). In fact, from 13 deregulated genes grouped in categories of cytoskeleton and transport described in the presymptomatic spinal cord of a late onset ALS animal model (Guipponi et al., 2010), the down regulation of kinesin light chain 2 (Klc2) gene was in agreement to our analysis. Also, from 11 deregulated genes of the cytoskeleton and motor activity categories described in microdissected toluidine blue-labeled neurons from the spinal cord of asymptomatic SOD1 G93A mice (Ferraiuolo et al., 2007), the upregulated Kif1b and Gsn were found, respectively, in spinal cord of 40 and 80 days old mice of our analysis. Furthermore, similar study that employed a categorization of enriched pathways in asymptomatic ALS mice has not identified cytoskeleton category but pointed six differentially expressed cytoskeleton genes (Kudo et al., 2010).
It is interesting the robust regulation of gene expression in the spinal cord of the ALS mouse model in the presymptomatic phases of the disease, thus underlining the early events on modulation of cytoskeleton elements before the death of spinal cord motor neurons. The genes of cellular components related to cytoskeleton pointed by DAVID enriched analysis showed 14 and 76 deregulated genes in the spinal cord of 40 and 80 days old presymptomatic animals respectively, and 146 deregulated genes in the sciatic nerve of 60 days presymptomatic ALS mice.
In fact, early events regarding impairment of axonal transport, cytoplasm aggregation and neurite/axonal abnormalities have been described before the onset of ALS symptoms (Warita et al., 1999;Williamson and Cleveland, 1999;Magrane and Manfredi, 2009;Rothstein, 2009) and might contribute to the triggering of motor neuron death. Additional novelty of the present analysis was the demonstration of deregulated transcripts-related to cytoskeleton by DAVID categorization in the sciatic nerve of the presymptomatic ALS mice. Gene profiling in the sciatic nerve accounts for axonal (minority) and Schwann cell (majority) transcripts (Baraban et al., 2013;Malmqvist et al., 2014), thus adding an important contribution to molecular analysis on dying back hypothesis of ALS pathology (Dadon-Nachum et al., 2011;De Oliveira et al., 2014).
Despite the evidence of protein synthesis in growing or regenerating axons in vitro and of mRNA axonal transport and local translation in developing zebrafish (Baraban et al., 2013), it is still not clear whether axons of adult motor neurons contain ribosomes and other elements that are necessary for protein translation (Jablonka et al., 2004). Altogether, this manuscript highlights the importance to evaluate the cytoskeleton changes in the Schwann cells in the still poorly described peripheral pathology in ALS (Xiao et al., 2006) and their contribution to impair paracrine trophic actions to motor neurons (unpublished results from our laboratory, presented as an abstract form in the 2013 Society for Neuroscience Meting, San Diego, USA). In order to address this issue, we have developed and presented here different methods to enrich specific cell types for molecular analysis.

DEREGULATED GENES OF CYTOSKELETON MOLECULES ALREADY POINTED IN ALS
The description of 14 differentially expressed genes (six down and eight upregulated) in the spinal cord of 40 days old presymptomatic ALS mice already indicates a very early presymptomatic event related to cytoskeleton with a possible implication to physiopathological mechanisms of the disease onset. From those genes, only Kif3a and Kif1b or related molecules have been studied in the context of kinesin dysfunction or impaired anterograde transport of cargos, like neurofilament, in ALS (Dupuis et al., 2000;Conforti et al., 2003;Pantelidou et al., 2007). The Kif3a (downregulation) and Kif1b (upregulation) deregulated genes seen in the spinal cord of 40 days old SOD1 G93A mice were in line with descriptions of reduction KIF3Aβ in motor cortices of ALS human and animal model (Pantelidou et al., 2007) and of KIF3-associated proteins in ALS rodents (Dupuis et al., 2000), thus underlining the presence of an early and complex mechanism involved in the impairment of the fast anterograde axonal transport machinery in ALS prior motor neuron degeneration. From the 76 deregulated genes (25 down and 51 upregulated) in the spinal cord of 80 days old presymptomatic ALS mice, only eight genes (Actg1, Adora, Akt1, App,Dctn1,Kif1a,Sirt2,and Stmn1) or related molecules have been studied in the context of ALS.
The upregulation of Kif1a is in accordance to previous description (Dupuis et al., 2000) and might represent a regulatory mechanism in order to compensate the impaired anterograde transport in neurons. KIF1 is divided into KIF1A, responsible for transport of synaptic vesicle precursors (Okada et al., 1995), and KIF1B, described above, a monomeric motor responsible for the anterograde transport of mitochondria (Nangaku et al., 1994).
Actg1 deregulation described in the spinal cord of 80 days old ALS mice was already mentioned in a previous ALS publication (Baciu et al., 2012) and the related actin product of the gene might impair dendritic spine plasticity with potential implication to motor neuron toxicity in ALS (Sunico et al., 2011). Such a mechanism may also involve an impairment of purinergic receptor-mediated actin cytoskeleton remodeling (Goldman et al., 2013), which is in line with the upregulation of Adora1 transcription codifying adenosine A1 receptor, described in this and previous studies (Gundlach et al., 1990).
Dctn1 expression was upregulated in presymptomatic spinal cord (80 days old) of ALS mouse model, thus denoting dynactin impairment as a mechanism in the presymptomatic phase of the disease. Mutation in Dctn1 gene has been associated to motor neuron degeneration in ALS (Hafezparast et al., 2003) and downregulation of the gene was described in residual motor neurons of postmortem material of ALS patients (Jiang et al., 2007).   There is a lack of information on dynactin regulation before clinical onset of ALS despite the fact that dynein-dynactin complex, the only retrograde transport motor, contributes to formation of SOD1 inclusions in the disease (Strom et al., 2008). Furthermore, deregulation of Schwann cell genes related to neurothropindependent mechanisms in association to an impaired axonal transport may enhance motor neuron vulnerability in ALS (Koh et al., 2005;Niewiadomska et al., 2011). The upregulation of Sirt2 expression in the spinal cord of 80 days old mice might trigger toxicicity to motor neurons in the late stage of the presymptomatic age by increasing deacetylation of alfa tubulin (Korner et al., 2013;Taes et al., 2013). Moreover, the upregulation of Stmn1 in the ALS spinal cord before neuronal death is in agreement to previous description on stathmin protein accumulation in spinal cord motor neurons leading to Gogi apparatus fragmentation and collapse of microtubule network (Strey et al., 2004). Furthermore, the formation of perikaryal/axonal intermediate filament inclusions, neurofilament abnormalities and genetic defects in microtubule-based transport that may facilitate the elevation of the toxic amyloid beta precursor in ALS (Spadoni et al., 2009;Bryson et al., 2012) might correlate the downregulation of App seen in the spinal cord of 80 days presymptomatic phase to a neuroprotective regulation before the neuronal death onset.
Finally, the downregulation of Vapa in the spinal cord of 80 days old ALS mice is also an original and interesting finding of the present analysis. VAMP/synaptobrevin-associated proteins A and B (VAPA and VAPB) are both enriched on endoplasmic reticulum and Golgi membranes and are capable to interact with cytoskeleton elements in order to maintain the organelle morphology (Nishimura et al., 1999). Despite a lack of information on VAPA function in the central nervous system, VAPB mutation is associated to a familial form of ALS (Nishimura et al., 2004).
It should be pointed out the regulation Maea, the unique gene that was overexpressed in the spinal cord of both 40 and 80 days old ALS mice, indicating its long lasting involvement in the presymptomatic events in the ALS spinal cord. The role of Maea in ALS is unknown but it could participate in the immunomodulatory signaling of non-neuronal cells-induced toxicity in ALS (Levine et al., 1999;McGeer and McGeer, 2002;Pasinelli and Brown, 2006).
From the 146 deregulated genes (74 down and 72 upregulated) in the sciatic nerve of 60 days presymptomatic old SOD1 G93A old mice, only 10 genes or related molecules have been studied in context of ALS.
The downregulation of the Aif1, Ccnb1, and Mapt in the presymptomatic ALS mice is likely to participate in the early events in the ALS peripheral nerve pathology. Aif1 encodes the allograft inflammatory factor-1 (AIF-1) and AIF-1 positive microglia/macrophages are among the earliest cells to respond to nerve injury (Schwab et al., 2001). It is likely that AIF-1 may act as an initiator of the early microglial/macrophage-induced immunomodulation leading a motor axon retraction and neuromuscular junction disruption before neuronal degeneration (Dibaj et al., 2011). Furthermore, Ccnb1 deregulation-induced cytoskeleton disorganization (Husseman et al., 2000) and also altered neuronal cytoskeleton protein Tau encoded by Maptinduced microtubule stabilization and assembly deregulation  (Aronov et al., 2002) are possible mechanisms related to early axonal retraction taking place in presymptomatic phases of the disease (Aronov et al., 2002). Actn3, Als2, Kif5a, Kif5c, Nos3, and Tmod3 were found upregulated in the sciatic nerve of 60 days old presymptomatic ALS mice and their related molecules have been mentioned in the context of ALS mechanisms. ACTN3, one of the four human alpha-actinin isoforms, has been associated to ALS progression in human muscle (Pradat et al., 2012;Bernardini et al., 2013). Moreover, the upregulated Tmod3, which codifies the tropomodulin (TMOD), an actin-capping protein for the slow-growing end of filamentous actin (Ito et al., 1995), may represent a need for the dynamic polymerization of actin cytoskeleton, probably in the Schwann cells of ALS nerve. Interestingly, the mutation of ubiquitously expressed TMOD3 protein is responsible for type 5 familial ALS (Cox and Zoghbi, 2000). Furthermore, despite ALS2 deficiency accounts for ALS2 familial form (Hadano et al., 2001;Yang et al., 2001), the upregulation of Als2 seen in the 60 days old presymptomatic ALS sciatic nerve could reflect a transient Schwann cell neuroprotective paracrine response (Hadano et al., 2010). In fact, ALS2/alsin, a guanine nucleotide exchange factor for GTPase Rab5, is involved in endosome fusion/trafficking, neurite outgrowth and corticospinal axon integrity (Deng et al., 2007;Lai et al., 2009), probably by interfering with the accumulation of immature vesicles and misfolded proteins (Lai et al., 2009).
We have also found a deregulation of the Nos2 and Nos3 expression in sciatic nerve of 60 days old ALS mice, without changes in presymptomatic spinal cord. The synthesis of inducible NOS in the spinal cord and peripheral nerve of ALS model in the presymptomatic phase of the disease has been mentioned (Almer et al., 1999;Chen et al., 2010).
Important finding of the present study was also the regulation of several genes of kinesin molecules in the sciatic nerve of 60 days old presymptomatic ALS mice. The differential Kif5 regulation in the sciatic nerve observed in our microarray analysis may in fact represent the gene regulation in the Schwann cells where the KIF5 participate in the myelin integrity (Bolis et al., 2009). In this context, KIF5B might be of substantial interest because it also expresses in non-neuronal cells and its regulation/activity has not been explored in ALS mechanisms. Further evidence for KIF5 mechanisms in non-neuronal cells were obtained from the absence of Kif5 expression in the spinal cord of presymptomatic SOD1 G93A mice by gene profiling study (this work) and qPCR analysis (Kuzma-Kozakiewicz et al., 2013).

TOPOGRAPHIC AND CELLULAR MODULATION OF Kif1b OF KINESIN FAMILY
The downregulation of Kif1b in the sciatic nerve showed by microarray and qPCR was in the opposite direction to the upregulation of the gene seen in the spinal cord of presymptomatic ALS mice. Importantly, only the Kif1b was deregulated in the two evaluated regions and also was reported previously in the context of ALS (Ferraiuolo et al., 2007;Kuzma-Kozakiewicz et al., 2013). Furthermore, it is the first time Kif1b expression or its protein has been described in Schwann cells. The development of technology to enrich Schwann cells allowed the present analysis.
Kif1b deregulation in presymptomatic ALS seems to be an important event, specially in the Schwann cells because KIF1B is required for adequate myelination process by oligodendrocytes (Lyons et al., 2009). Presently, myelin pathology is not clear in ALS peripheral nerves (Heads et al., 1991). Nonetheless, it is still undefined whether peripheral myelin morphological alteration in ALS is a consequence of axonal degeneration (Perrie et al., 1993).
There is a lack of information on myelin alterations in peripheral nerves of presymptomatic ALS mice and it is uncertain whether impairments in the Schwann cell function could contribute to ALS axonal pathology and dying back events. That is actually an important issue in the pathogenic mechanism of the disease because myelin cell function overtakes action potential conduction along peripheral axons (Monk and Talbot, 2009). In fact, the control of axoplasmic Ca 2+ and posttranslational modifications of local trafficking proteins are part of trophic support signaling provided by myelinating cells. We might speculate that an impaired crosstalk between Schwann cell and motor axon in presymptomatic stages of the disease could trigger axonal retraction and Wallerian degeneration (Lyons et al., 2009;Kiryu-Seo and Kiyama, 2011;Gentil and Cooper, 2012).
The above discussion indicates the importance to evaluate the events in specific cell types notably in the context of peripheral nerve-induced neuropathology in ALS, because a differential regulation can occur specifically in Schwann cells and also altered axonal transport might modify the local traffic of RNAs (Ticozzi et al., 2010).
Single cell laser microdissection has been employed by our group and other researchers to evaluate gene expression on enriched cell types (Ferraiuolo et al., 2007;De Oliveira et al., 2009;Guipponi et al., 2010;Kudo et al., 2010;Tanaka et al., 2012). We have developed and showed here for the first time the method to immunolabel mouse Schwann cells and motor neurons and also the procedures to enrich cell samples for molecular analyses by means of single cell laser microdissection. We have also developed and demonstrated in this work the methodology to enrich mouse Schwann cells by means of primary cell culture and cell sorting. We still do not know precisely the limitations of the enrichment techniques employed here. Nevertheless, it should be emphasized the agreement of the results obtained by the different methods.
It is of substantial importance that the methodology allowed the observation of a differential regulation of Kif1B in the enriched ChAT immunolabeled motor neurons (upregulation) and Schwann cells (downregulation). The results were coincident to those obtained in qPCR and microarray analyses of whole spinal cord and sciatic nerve of presymptomatic ALS mice.
A differential gene regulation in specific cell types in a neuroglial unit, the motor neuron-Schwann cell unit in this case, highlights the complexity of cellular and molecular mechanisms of ALS, remarkably before clinical onset. Kif1b upregulation in immunolabeled motor neurons was in line to a previous work that employed toluidine blue-enriched putative motor neurons of presymptomatic ALS mice (Ferraiuolo et al., 2007). This finding indicates an involvement of a motor protein of kinesin family in the axonal trafficking before the death of motor neurons and the appearance of neurological symptoms. The elevation of the protein in ALS neurons might be a substrate for an increased kinesin-1 phosphorylation and a diminution of kinesin-1 function with a subsequent defect of fast axonal transport (Morfini et al., 2013).
The Kif1b downregulation in the sciatic nerve and also in the enriched Schwann cells of presymptomatic ALS mice is a major original contribution of the present work. The recent description on the role of KIF1B for the adequate function of central myelinating cells (Lyons et al., 2009;Gentil and Cooper, 2012) opens up the possibly for the existence of KIF1B mechanisms in the paracrine trophic actions of Schwann cells to peripheral motor neurons. Deregulated KIF1B in Schwann cells highlights the possibility of no autonomous cell toxicity of Schwann cells to motor neurons in ALS, mechanisms that should be investigated in details in future works. In fact, the non-autonomous cell toxicity of central glia to motor neurons has been described (Boillee et al., 2006a) and the related molecular pathways are under investigation.
In conclusion, the present work demonstrated cytoskeleton gene regulation as an important occurrence in motor neurons and Schwann cells in the presymptomatic stages of ALS and may be of importance in the dying back mechanisms of neuronal death in the neurodegenerative disease. The differential regulation of Kif1b in the spinal cord (upregulation) and sciatic nerve (downregulation) was coincident to that found in the enriched motor neurons and Schwann cells and emerged as an important event in the pathogenic mechanism of ALS.