# Gene expression changes in spinal motoneurons of the SOD1(G93A) transgenic model for ALS after treatment with G-CSF

Alexandre Henriques, Stefan Kastner, Eva Chatzikonstantinou, Claudia

Pitzer, Christian Plaas, Friederike Kirsch, Oliver Wafzig, Carola Krüger, Robert Spoelgen, Jose-Luis Gonzalez de Aguilar, et al.

► To cite this version:

Alexandre Henriques, Stefan Kastner, Eva Chatzikonstantinou, Claudia Pitzer, Christian Plaas, et al.. Gene expression changes in spinal motoneurons of the SOD1(G93A) transgenic model for ALS after treatment with G-CSF. Frontiers in Cellular Neuroscience, 2015, 8 (Art. n°464), 10.3389/fncel.2014.00464. hal-05255543

HAL Id: hal-05255543

https://univoak.hal.science/hal-05255543v1

Submitted on 15 Sep 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

# Gene expression changes in spinal motoneurons of the SOD1<sup>G93A</sup> transgenic model for ALS after treatment with G-CSF

Alexandre Henriques<fnref n="1" />,<fnref n="2" />†, Stefan Kastner<fnref n="3" />†, Eva Chatzikonstantinou<fnref n="3" />, Claudia Pitzer<fnref n="3" />, Christian Plaas<fnref n="3" />, Friederike Kirsch<fnref n="3" />, Oliver Wafzig<fnref n="3" />, Carola Krüger<fnref n="3" />, Robert Spoelgen<fnref n="3" />, Jose-Luis Gonzalez De Aguilar<fnref n="1" />,<fnref n="2" />, Norbert Gretz<fnref n="4" /> and Armin Schneider<fnref n="3" />*

**Edited by:**

Shawn Hayley, Carleton University,
Canada

**Reviewed by:**

Rafael Linden, Federal University of
Rio de Janeiro, Brazil
Florence Evelyne PERRIN,
University of Montpellier, France

*Correspondence:

Armin Schneider, Sygnis Bioscience,
Im Neuenheimer Feld 515, 69120
Heidelberg, Germany
e-mail: schneider@sygnis.de

† These authors have contributed
equally to this work.

**Background:** Amyotrophic lateral sclerosis (ALS) is an incurable fatal motoneuron disease with a lifetime risk of approximately 1:400. It is characterized by progressive weakness, muscle wasting, and death ensuing 3–5 years after diagnosis. Granulocyte-colony stimulating factor (G-CSF) is a drug candidate for ALS, with evidence for efficacy from animal studies and interesting data from pilot clinical trials. To gain insight into the disease mechanisms and mode of action of G-CSF, we performed gene expression profiling on isolated lumbar motoneurons from SOD1<sup>G93A</sup> mice, the most frequently studied animal model for ALS, with and without G-CSF treatment.

**Results:** Motoneurons from SOD1<sup>G93A</sup> mice present a distinct gene expression profile in comparison to controls already at an early disease stage (11 weeks of age), when treatment was initiated. The degree of deregulation increases at a time where motor symptoms are obvious (15 weeks of age). Upon G-CSF treatment, transcriptomic deregulations of SOD1<sup>G93A</sup> motoneurons were notably restored. Discriminant analysis revealed that SOD1<sup>G93A</sup> mice treated with G-CSF has a transcriptom close to presymptomatic SOD1<sup>G93A</sup> mice or wild type mice. Some interesting genes modulated by G-CSF treatment relate to neuromuscular function such as CCR4-NOT or Prss12.

**Conclusions:** Our data suggest that G-CSF is able to re-adjust gene expression in symptomatic SOD1<sup>G93A</sup> motoneurons. This provides further arguments for G-CSF as a promising drug candidate for ALS.

**Keywords:** motoneuron, ALS, G-CSF, neurodegeneration, mouse model, gene expression, laser microdissection

## INTRODUCTION

Amyotrophic lateral sclerosis (ALS) is a devastating neurodegenerative disease that results in progressive loss of motoneurons, motor weakness, and death within 1–5 years after disease onset. Although more than 130 years have passed since the disease was first described by Charcot its etiology is still enigmatic, precluding a truly rational drug development. Motoneurons in the spinal cord and the motor unit are the common endpoint in all hypotheses concerning pathophysiology (Gordon, 2013; Turner et al., 2013). The transgenic SOD-1<sup>G93A</sup> mice are a model for ALS and recapitulate most of the ALS pathology, including loss of motor units. The first evidence of stress on the motor axis in this model occurs early before appearance of motor deficit (Frey et al., 2000).

G-CSF was primarily described as a hematopoietic factor able to stimulate the proliferation and differentiation of myeloid precursors. It is used for the treatment of chemotherapy-induced neutropenia and for bone marrow transplants (Kobbe et al., 2009;

Yang et al., 2012). G-CSF has more recently been uncovered as a neurotrophic factor able to rescue motoneurons from apoptosis, protect motor functions and enhance overall survival of SOD1<sup>G93A</sup> mice through the preservation of motor units (Pitzer et al., 2008; Henriques et al., 2010a,b, 2011).

Pilot clinical studies have investigated the potential of G-CSF as a drug candidate in ALS. Administrations were mainly subcutaneous with limited duration (unique or repeated injection for up to 5 days) and enrollment ranged from 10 to 39 patients. Adverse effects were minor and the drug was well tolerated in all studies. Effects of G-CSF on disease progression are unclear since all studies suffer from small cohort sizes. One study concluded that G-CSF treatment did not bring benefit for patients (Nefussy et al., 2010), while others either reported slower disease progression upon treatment (Martinez et al., 2009; Zhang et al., 2009), reduced damage of white matter tracts (Duning et al., 2011) or lower level of pro-inflammatory markers in the cerebrospinal fluid (Chio et al., 2011).

<footnote n="1">¹ INSERM, U1118, Mécanismes Centraux et Péripheriques de la Neurodégénérescence, Strasbourg, France</footnote>

<footnote n="2">² UMRS1118, Fédération de Médecine Translationnelle de Strasbourg, Université de Strasbourg, France</footnote>

<footnote n="3">³ Sygnis Bioscience GmbH & Co KG, Heidelberg, Germany</footnote>

<footnote n="4">⁴ Medical Research Center, Medical Faculty Mannheim, University of Heidelberg, Mannheim, Germany</footnote>
Effects of G-CSF treatment in SOD-1<sup>G93A</sup> mice have been examined at the clinical, electrophysiological, and histological level. Here we sought to determine if and what changes on the gene expression level in alpha motoneurons would be related to G-CSF treatment. Therefore, our objective was to (1) define changes on the mRNA level that occur in spinal motoneurons in an animal model of ALS, the SOD1<sup>G93A</sup> mice and (2) characterize the effects of G-CSF on detected changes in the transcriptome of those motoneurons. We combined laser microdissection (LMD), RNA amplification, and gene expression profiling on DNA arrays to study these cell-type specific expression changes.

# MATERIALS AND METHODS

## ALS MODEL

We used mice transgenic for the SOD1<sup>G93A</sup> mutation (Gurney et al., 1994) on a C57/Bl6 background (B6.Cg-Tg(SOD1-G93A)1Gur/J strain; Jackson laboratory, Bar Harbor, Maine, USA), harboring the high copy number of the mutant allele human SOD1. The hemizygous line was maintained by mating transgenic males with C57BL/6 wild-type females. Transgenic littermate females were used in all experiments. They were either treated with vehicle (250 mM sorbitol, 0.004% Tween-80 and 10 mM sodium acetate buffer) or with G-CSF (filgrastim dissolved in vehicle, 30 µg/kg body weight/day) with an osmotic, subcutaneously implanted, paravertebrally located minipump (Alzet Minipump, model 2004; ALZET Osmotic Pumps, Cupertino, CA, USA) (Pitzer et al., 2008). Wild type littermates served as control group (group 3). Each mouse was housed under standard conditions in cages together with 4-6 littermates and fed *ad libitum*, in a room that also housed male mice, thus resulting in likely estrous cycle synchronization. All animal experiments were approved by the Regierungspräsidium Karlsruhe (Karlsruhe, Germany).

## TISSUE ISOLATION

All kinds of surgical instruments (forceps, scissors) were cleaned with RNase Zap spray or 100 % ethanol before use. Directions for handling RNA or tissue for RNA isolation (e.g., gloves, etc.) were followed. For anesthetizing purpose, a mixture of ketamine (120 mg/kg boday mass) and xylasine (Rompun, 16 mg/kg body mass was injected intraperitoneally. After complete absence of peripheral nociceptive reflexes the animals were perfused transscardially via the left ventricle with ice-cold HBSS. The pump was run at 10 ml/min for 3 min. Care was taken not to injure any blood vessel during the preparation. The perfused mouse was fixed onto a metal platform cooled by wet ice to avoid RNA degradation. After removal of the skin, cranial bones and vertebrae were carefully removed dorsally and spinal cords were dissected from the remaining opened cranium and vertebral canal. The thoracic position T9 was marked with a blue dye spot. Tissue was frozen on dry ice for about 5–10 min and stored at –80°C.

## LASER MICRODISSECTION

Spinal cords were cryosectioned into 8 µm coronal sections, starting at the lumbar position L2. Thirty sections per mice were collected. Sections taken for microdissection were 100 µm apart.

The tissue on the POL-membrane of a frame slide was overlaid with thionin staining solution so that the whole sample was completely covered. After 20 s ultraPURE™ water was poured onto the sample to wash away excess thionin. The frame slide with the tissue sample on the POL-membrane was mounted onto the microscope. Samples were collected into the cap of a 0.65 ml LMD tube filled with 35 µl of buffer RLT (Qiagen RNeasy® Micro Kit). Using the program “Leica Laser Microdissection V 5.0,” cells were marked and then automatically cut by a moving laser beam (power 40–50, speed 4, specimen/balance 11, offset 45). A total of 300 motoneurons per mouse spinal cord have been sampled by laser microdissection (Supplementary Figure 1), from a total of 35 mice (week 11 wt: *n* = 6; SOD *n* = 8); week 15 [wt *n* = 7; SOD (vehicle): *n* = 7; SOD (G-CSF treated): *n* = 7)], and subjected amplified RNA to array hybridization. Cutting was done at 40x magnification. Experimenters involved in laser microdissection were blind to the identity of the animal.

## RNA AMPLIFICATION

Each RNA sample was linearly amplified by two consecutive rounds of T7 polymerase-driven amplification. The precipitated total RNA of each LMD sample was reverse transcribed with a first-strand synthesis mix comprising 10 U Superscript III (Invitrogen, Karlsruhe, Germany), 1 U Superasin (Ambion, Huntingdon, UK) and 0.5 pmol of T7-T18V primer by incubating for 2 h at 50°C after second-strand synthesis with 6 U DNA polymerase I, 2 h at 22°C with 0.2 U RNaseH, and 15 min at 16°C with 4 U T4 DNA polymerase I (Invitrogen). The T7 promoter tagged cDNA template was purified by standard phenol extraction and ethanol precipitation. This cDNA was used for first round amplified antisense RNA production applying a MegaScript kit (Ambion) and purified with an RNeasy Mini Kit (Qiagen) according to the manufacturer's instructions. Second round amplification was performed accordingly but using T7-tagged random primer for first strand cDNA synthesis and T18V primer for second-strand synthesis, yielding amplified sense RNA. UV measurements at 280 and 260 nm and 2100 Bioanalyzer runs using the RNA-6000 Nano-LabChip (Agilent, Böblingen, Germany) were used to ensure the quality and quantity of the sense RNA produced. Only samples showing RNA integrity numbers (RIN) above 8 were used. The yield varied from 40 to 80 µg.

## ARRAY ANALYSES

The arrays were performed with the GeneChip® mouse genome 430 2.0 array from Affymetrix (Affymetrix, Santa Clara, CA, USA). Expression console software (Affymetrix, Santa Clara, CA, USA) was used to normalize arrays with the gcrma algorithm (Harr and Schlotterer, 2006) and to perform quality control. The database was filtered with intensity and standard deviation. We removed genes when one of the following conditions was not met: less than 4 arrays had intensity higher than the centile 20% of all genes or the standard deviation was lower than centile 30% of all genes.

Discriminant analysis has been performed with JMP 11.0.0 (SAS Institute, Cary, North Carolina) with genes presenting a minimum standard deviation greater than 1.2 in all individuals.
Class prediction has been computed using a non-parametric
k-Nearest Neighbors algorithm.

Differential gene expression among experimental groups has
been assessed by ANOVA followed by a post-hoc Tukey's Honest,
with genes presenting a minimum fold-change of 1.5-fold. We
considered a *p*-value lower than 0.05 with a FDR lower than 10%
as significant difference. Analyses were performed using BRB-
ArrayTools developed by Dr. Richard Simon and BRB-ArrayTools
Development Team (http://linus.nci.nih.gov).

WEB-based Gene SeT AnaLysis Toolkit was used to detect
enrichment of biological functions (WebGestalt, http://bioinfo.
vanderbilt.edu/webgestalt/) (Zhang et al., 2005; Wang et al.,
2013).

REAL-TIME QUANTITATIVE PCR

One µg total RNA used for the microarray was used to synthe-
size cDNA using superscript II reverse transcriptase (Invitrogen)
and oligo-dT primers, accordingly to manufacturer's recom-
mendation. Quantitative reverse transcription-polymerase chain
reaction was performed using the Lightcycler system with SYBR-
Green reagent (Roche). Cycling conditions were as follows:
10 min at 95°C; 5 s at 95°C, 10 s at the annealing temperature of
64°C, 30 s at 72°C for 50 cycles and 10 min at 95°C. Relative quan-
tification were calculated after normalization to the reference gene
cyclophilin B.

RESULTS

DESIGN OF THE EXPERIMENT

A first group of SOD1<sup>G93A</sup> and WT mice was included in the study at week 11 of age when SOD1<sup>G93A</sup> mice present no signs of motor dysfunction but subtle signs of denervation detectable by electromyography (Pitzer et al., 2008). The second cohort of mice was treated with G-CSF or vehicle from week 11 to week 15. At the time of study completion, SOD1<sup>G93A</sup> mice presented clear motor impairment and motoneuron degeneration is documented (Pitzer et al., 2008; Henriques et al., 2011).

During laser microdissection, there is a risk of collecting cells
surrounding motoneurons and therefore analyzing non-neuronal
transcripts (Supplementary Figure 1A, Perrin et al., 2005;
Bandyopadhyay et al., 2013). To assess the proportion of con-
taminating RNA, we have studied the expression values of highly
specific markers for astrocytes, oligodendrocytes, microglia,
endothelial cells, and neurons (Cahoy et al., 2008). All exper-
imental groups show high expression of the specific neuronal
markers, while markers for glial and endothelial cells are present
at very low levels or not detected (Supplementary Figure 1B).
These data suggest that contamination with RNA from other cells
than motoneurons is minimal and irrelevant for the following
analyses.

The design of the present study should be informative on genes
altered in motoneurons of SOD1<sup>G93A</sup> mice from the clinically
non-symptomatic to an early symptomatic stage, and give insight
into genes influenced by G-CSF treatment (Figure 1A). Results
were analyzed by unsupervised and supervised approaches
(Figure 1B), and biological functions for relevant genes analyzed
based on gene ontology.

<figure><img src="image_3.png" /></figure>

GENE REGULATION IN SOD1<sup>G93A</sup> MOTONEURONS

At 11 weeks of age, when SOD1<sup>G93A</sup> mice do not show motor
symptoms despite early electromyographic changes we detected
985 transcripts with a significantly altered expression, of which
725 genes were upregulated and 260 dowregulated (Table 1).
At 15 weeks of age, when motor symptoms are obvious and
muscular denervation signs are prominent by electromyogra-
phy, the number of significantly deregulated transcripts was two
times higher, with 790 genes upregulated and 1245 dowregulated
(Table 2). Functional annotation using gene ontology revealed
that metabolic process, response to stimulus/stress and nucleic
acid binding were the main biological and molecular functions
altered in SOD-1 motoneurons (Supplementary Figure 2).

At the presymptomatic disease stage (week 11), genes which
were most downregulated were the collectin sub-family member
11 (colec11) and vacuolar protein sorting 4b (vpsb4) (Table 1).
Expression of colec11 has been reported in neurons and is
involved in innate immunity (Motomura et al., 2008; Ohtani
et al., 2012). vpsb4 is known to interact with the pro-apoptotic
pathway (Cui et al., 2014) and the autophagosome in neurode-
generative conditions (Han et al., 2012). The most strongly upreg-
ulated genes were related to response to stress and injury, such
as activating transcription factor 3 (atf3) and small proline-rich
<table><caption>Table 1 | List of the main transcripts whose expression are altered in SOD1<sup>G93A</sup> mice at 11 week of age.</caption><thead><tr><th>Probe set ID</th><th>Symbol</th><th>Name</th><th>FC</th><th>p-value</th></tr></thead><tbody><tr><td>1458345_s_at</td><td>Colec11</td><td>Collectin sub-family member 11</td><td>0.28</td><td>0.005</td></tr><tr><td>1457067_at</td><td>Vps4b</td><td>Vacuolar protein sorting 4b (yeast)</td><td>0.32</td><td>0.016</td></tr><tr><td>1444869_at</td><td>–</td><td>–</td><td>0.34</td><td>0.001</td></tr><tr><td>1427351_s_at</td><td>Ighm</td><td>Immunoglobulin heavy constant mu</td><td>0.34</td><td>0.000</td></tr><tr><td>1419332_at</td><td>Egfl6</td><td>EGF-like-domain, multiple 6</td><td>0.38</td><td>0.000</td></tr><tr><td>1423608_at</td><td>Itm2a</td><td>Integral membrane protein 2A</td><td>0.38</td><td>0.000</td></tr><tr><td>1438763_at</td><td>Dnahc2</td><td>Dynein, axonemal, heavy chain 2</td><td>0.39</td><td>0.000</td></tr><tr><td>1435815_at</td><td>Ldoc1</td><td>Leucine zipper, down-regulated in cancer 1</td><td>0.40</td><td>0.000</td></tr><tr><td>1423390_at</td><td>Siah1a</td><td>Seven in absentia 1A</td><td>0.41</td><td>0.000</td></tr><tr><td>1418916_a_at</td><td>Spp2</td><td>Secreted phosphoprotein 2</td><td>0.43</td><td>0.000</td></tr><tr><td>1441114_at</td><td>9330156P08Rik</td><td>RIKEN cDNA 9330156P08 gene</td><td>0.44</td><td>0.004</td></tr><tr><td>1420386_at</td><td>Seh1l</td><td>SEH1-like</td><td>0.45</td><td>0.005</td></tr><tr><td>1419435_at</td><td>Aox1</td><td>Aldehyde oxidase 1</td><td>0.46</td><td>0.001</td></tr><tr><td>1445767_at</td><td>Ptprd</td><td>Protein tyrosine phosphatase, receptor type, D</td><td>0.46</td><td>0.001</td></tr><tr><td>1452276_at</td><td>Smarcad1</td><td>SWI/SNF-related, matrix-associated actin-dependent regulator of chromatin, subfamy a, containing DEAD/H box1</td><td>0.46</td><td>0.013</td></tr><tr><td>1437796_at</td><td>–</td><td>–</td><td>0.46</td><td>0.000</td></tr><tr><td>1427196_at</td><td>Wnk4</td><td>WNK lysine deficient protein kinase 4</td><td>0.46</td><td>0.016</td></tr><tr><td>1451440_at</td><td>Chodl</td><td>Chondrolectin</td><td>0.46</td><td>0.000</td></tr><tr><td>1439011_at</td><td>Gm20245</td><td>Predicted gene, 20245</td><td>0.46</td><td>0.044</td></tr><tr><td>1438143_s_at</td><td>Atxn2</td><td>Ataxin 2</td><td>0.47</td><td>0.000</td></tr><tr><td>1449662_at</td><td>Reep1</td><td>Receptor accessory protein 1</td><td>4.0</td><td>0.0000</td></tr><tr><td>1434442_at</td><td>Stbd1</td><td>Starch binding domain 1</td><td>4.0</td><td>0.0001</td></tr><tr><td>1426851_a_at</td><td>Nov</td><td>Nephroblastoma overexpressed gene</td><td>4.1</td><td>0.0003</td></tr><tr><td>1426808_at</td><td>Lgals3</td><td>Lectin, galactose binding, soluble 3</td><td>4.2</td><td>0.0000</td></tr><tr><td>1440142_s_at</td><td>Gfap</td><td>Glial fibrillary acidic protein</td><td>4.4</td><td>0.0000</td></tr><tr><td>1419665_a_at</td><td>Nupr1</td><td>Nuclear protein 1</td><td>4.5</td><td>0.0000</td></tr><tr><td>1423427_at</td><td>Adcyap1</td><td>Adenylate cyclase activating polypeptide 1</td><td>4.9</td><td>0.0002</td></tr><tr><td>1449519_at</td><td>Gadd45a</td><td>Growth arrest and DNA-damage-inducible 45 alpha</td><td>5.1</td><td>0.0000</td></tr><tr><td>1417022_at</td><td>Slc7a3</td><td>Solute carrier family 7, member3</td><td>5.1</td><td>0.0000</td></tr><tr><td>1451285_at</td><td>Fus</td><td>Fusion, derived from t(12;16) malignant liposarcoma (human)</td><td>5.6</td><td>0.0060</td></tr><tr><td>1417868_a_at</td><td>Ctsz</td><td>cathepsin Z</td><td>5.7</td><td>0.0000</td></tr><tr><td>1422916_at</td><td>Fgf21</td><td>Fibroblast growth factor 21</td><td>6.1</td><td>0.0000</td></tr><tr><td>1437621_x_at</td><td>Phgdh</td><td>3-phosphoglycerate dehydrogenase</td><td>8.3</td><td>0.0000</td></tr><tr><td>1437232_at</td><td>Bpifc</td><td>BPI fold containing family C</td><td>8.4</td><td>0.0000</td></tr><tr><td>1454714_x_at</td><td>Phgdh</td><td>3-phosphoglycerate dehydrogenase</td><td>8.7</td><td>0.0000</td></tr><tr><td>1434129_s_at</td><td>Lhfpl2</td><td>Lipoma HMGIC fusion partner-like 2</td><td>10.6</td><td>0.0000</td></tr><tr><td>1449153_at</td><td>Mmp12</td><td>Matrix metallopeptidase 12</td><td>12.3</td><td>0.0002</td></tr><tr><td>1449133_at</td><td>Sprr1a</td><td>small proline-rich protein 1A</td><td>13.7</td><td>0.0000</td></tr><tr><td>1449363_at</td><td>Atf3</td><td>activating transcription factor 3</td><td>14.4</td><td>0.0000</td></tr></tbody></table>

repeat protein 1A (sprr1A). Upregulation of atf3 has already been noted in motoneurons after axotomy (Tsujino et al., 2000) and in response to deficiency in neurotrophic factors (Hyatt Sachs et al., 2007). Sprr1A was the second most stimulated gene and, similar to atf3, its expression is enhanced in motoneurons upon axotomy in order to promote axonal growth (Bonilla et al., 2002; Li and Strittmatter, 2003; Starkey et al., 2009). Interestingly, the fus (fused in sarcoma) gene was strongly upregulated in the presymptomatic disease stage. Fus protein is a RNA-binding protein and has important functions in RNA metabolism (e.g. splicing, formation of stress granule, transport and local translation). Mutations in fus are found in familial forms of ALS, and

inclusions of the protein are present in patients with frontotemporal dementia, a syndrome that is connected to ALS (Ling et al., 2013).

At the symptomatic disease stage, the most overexpressed genes were mainly associated with motoneuron stress response and support of axonal regeneration, similar to the presymptomatic disease stage (Table 2). Beside sprr1A and atf3, the growth associated protein 43 (gap43) shows increased expression (9-fold) which is involved in axonal regeneration after motoneuron injury (Chen et al., 2010; Grasselli et al., 2011). Interestingly, the transcription factor responsible for sprr1A expression in motoneurons, sox11, is strongly upregulated at
<table><caption>Table 2 | List of the main transcripts whose expression are altered in SOD1<sup>G93A</sup> mice at 15 week of age.</caption><thead><tr><th>Probe set ID</th><th>Symbol</th><th>Name</th><th>FC</th><th>p-value</th></tr></thead><tbody><tr><td>1457587_at</td><td>Kcnq5</td><td>Potassium voltage-gated channel, subfamily Q, member5</td><td>0.08</td><td>0</td></tr><tr><td>1423608_at</td><td>Itm2a</td><td>Integral membrane protein 2A</td><td>0.1</td><td>0</td></tr><tr><td>1427329_a_at</td><td>Ighm</td><td>Immunoglobulin heavy constant mu</td><td>0.12</td><td>0</td></tr><tr><td>1427351_s_at</td><td>Ighm</td><td>Immunoglobulin heavy constant mu</td><td>0.13</td><td>0</td></tr><tr><td>1451440_at</td><td>Chodl</td><td>Chondrolectin</td><td>0.13</td><td>0</td></tr><tr><td>1451047_at</td><td>Itm2a</td><td>Integral membrane protein 2A</td><td>0.14</td><td>0</td></tr><tr><td>1429759_at</td><td>Rps6ka6</td><td>Ribosomal protein S6 kinase polypeptide 6</td><td>0.15</td><td>0</td></tr><tr><td>1445268_at</td><td>-</td><td>-</td><td>0.15</td><td>0</td></tr><tr><td>1455291_s_at</td><td>Znrf2</td><td>Zinc and ring finger 2</td><td>0.17</td><td>0.001</td></tr><tr><td>1455238_at</td><td>Mum11l</td><td>Melanoma associated antigen (mutated) 1-like 1</td><td>0.17</td><td>0</td></tr><tr><td>1449155_at</td><td>Polr3g</td><td>Polymerase (RNA) III (DNA directed) polypeptide G</td><td>0.17</td><td>0</td></tr><tr><td>1433898_at</td><td>-</td><td>-</td><td>0.17</td><td>0</td></tr><tr><td>1441801_at</td><td>Kctd4</td><td>Potassium channel tetramerisation domain containing 4</td><td>0.18</td><td>0</td></tr><tr><td>1452366_at</td><td>Csgalnact1</td><td>Chondroitin sulfate n-acetylgalactosaminyltransferase 1</td><td>0.18</td><td>0</td></tr><tr><td>1418469_at</td><td>Nrip1</td><td>Nuclear receptor interacting protein 1</td><td>0.18</td><td>0.005</td></tr><tr><td>1434102_at</td><td>Nfib</td><td>Nuclear factor i/b</td><td>0.18</td><td>0</td></tr><tr><td>1433827_at</td><td>Atp8a1</td><td>ATPase, aminophospholipid transporter (APLT), class I, type 8A, member 1</td><td>0.19</td><td>0</td></tr><tr><td>1438214_at</td><td>Trps1</td><td>Trichorhinophalangeal syndrome I (human)</td><td>0.19</td><td>0.003</td></tr><tr><td>1422880_at</td><td>Sypl</td><td>Synaptophysin-like protein</td><td>0.19</td><td>0</td></tr><tr><td>1460227_at</td><td>Timp1</td><td>Tissue inhibitor of metalloproteinase 1</td><td>4.8</td><td>0.0001</td></tr><tr><td>1434376_at</td><td>Cd44</td><td>CD44 antigen</td><td>5.2</td><td>0.0013</td></tr><tr><td>1434129_s_at</td><td>Lhfpl2</td><td>Lipoma HMGIC fusion partner-like 2</td><td>5.6</td><td>0</td></tr><tr><td>1423427_at</td><td>Adcyap1</td><td>Adenylate cyclase activating polypeptide 1</td><td>5.6</td><td>0.0002</td></tr><tr><td>1426852_x_at</td><td>Nov</td><td>Nephroblastoma overexpressed gene</td><td>5.9</td><td>0</td></tr><tr><td>1448471_a_at</td><td>Ctla2a</td><td>Cytotoxic T lymphocyte-associated protein 2 alpha</td><td>5.9</td><td>0.0065</td></tr><tr><td>1429051_s_at</td><td>Sox11</td><td>SRY-box containing gene 11</td><td>6</td><td>0.0003</td></tr><tr><td>1454714_x_at</td><td>Phgdh</td><td>3-phosphoglycerate dehydrogenase</td><td>6.4</td><td>0</td></tr><tr><td>1437726_x_at</td><td>C1qb</td><td>Complement component 1, q subcomponent, beta polypeptide</td><td>6.6</td><td>0</td></tr><tr><td>1436996_x_at</td><td>Lyz1</td><td>Lysozyme 1</td><td>6.7</td><td>0.0054</td></tr><tr><td>1450792_at</td><td>Tyrobp</td><td>TYRO protein tyrosine kinase binding protein</td><td>7.1</td><td>0.0002</td></tr><tr><td>1437621_x_at</td><td>Phgdh</td><td>3-phosphoglycerate dehydrogenase</td><td>7.8</td><td>0</td></tr><tr><td>1423760_at</td><td>Cd44</td><td>CD44 antigen</td><td>7.9</td><td>0.0003</td></tr><tr><td>1434046_at</td><td>AA467197</td><td>Expressed sequence AA467197</td><td>8.3</td><td>0.0001</td></tr><tr><td>1426808_at</td><td>Lgals3</td><td>Lectin, galactose binding, soluble 3</td><td>8.9</td><td>0</td></tr><tr><td>1423537_at</td><td>Gap43</td><td>Growth associated protein 43</td><td>8.9</td><td>0</td></tr><tr><td>1449363_at</td><td>Atf3</td><td>Activating transcription factor 3</td><td>10</td><td>0</td></tr><tr><td>1440142_s_at</td><td>Gfap</td><td>Glial fibrillary acidic protein</td><td>13</td><td>0</td></tr><tr><td>1449133_at</td><td>Sprr1a</td><td>Small proline-rich protein 1A</td><td>20.6</td><td>0</td></tr></tbody></table>

this disease stage. Potassium voltage-gated channel, subfamily Q, member 5 (kcnq5) was the most downregulated gene (downregulated 14-fold) at the symptomatic disease stage. Its function in motoneurons is not described, but should be related to the control of resting and release properties of synapses (Huang and Trussell, 2011). The expression of chondrolectin (chodl) was also strongly repressed in motoneurons of SOD1<sup>G93A</sup> mice. Interestingly, this protein promotes motor axon growth and growth cone interactions, and its downregulation is associated with motoneuron degeneration in an animal model of spinal muscular atrophy (Zhong et al., 2012; Sleigh et al., 2014). As shown in Supplementary Table 1, 83 transcripts were commonly deregulated at both asymptomatic

and symptomatic disease stages. We noted deregulation at the two disease stages for several growth factors (FGF21, VGF, TGFa), and genes related to DNA/RNA metabolism (wars, eef1a1, eif4e).

## G-CSF PARTIALLY RESTORES THE MOTONEURON TRANSCRIPTIONAL SIGNATURE IN SOD1<sup>G93A</sup> MICE

In order to better understand the effects of G-CSF in the SOD1<sup>G93A</sup> mice, we have monitored transcriptional changes in motoneurons at a symptomatic disease stage upon a 4 weeks treatment (from week 11 to week 15). We performed a discriminant analysis with transcripts showing a standard deviation greater than 1.2 across all samples (734 transcripts).
As shown in a canonical plot, control groups are close to each other and SOD1<sup>G93A</sup> mice are distant from non-transgenic animals (Figure 2A). The SOD1 mice treated with G-CSF (green dots in Figure 2) tended to have a transcriptomic signature different from the sympatomatic SOD1 mice (red dots in Figure 2), that was closer to non-symptomatic mice. This observation was further confirmed by class prediction using a non-parametric k-Nearest Neighbors algorithm, used to predict the class of unknown samples. The algorithm was built with the signatures of all samples with the exclusion of mice treated with G-CSF. There were 4 defined classes: wild type 11 weeks, SOD-1<sup>G93A</sup> 11 weeks (presymptomatic), wild type 15 weeks and SOD-1<sup>G93A</sup> mice 15 weeks (symptomatic). SOD-1<sup>G93A</sup> mice treated with G-CSF were considered as unknown and the algorithm sorted them into the defined classes based on their transcriptome signature. The mean percent of correct classification reached 96% of efficacy for the defined classes, with a sensitivity above 0.875 and a specificity of 0.955.

The group treated with G-CSF was composed by seven SOD-1<sup>G93A</sup> mice. Three out of seven were sorted as symptomatic SOD-1<sup>G93A</sup>, one as asymptomatic SOD-1<sup>G93A</sup> 11 weeks, and three as wild type 15 weeks (Figure 2B). These data suggest that G-CSF can modulate the gene expression in motoneurons of 15 weeks old SOD-1<sup>G93A</sup> mice and can prevent disease-related changes to approach the transcriptomic signature of healthy mice or asymptomatic SOD1<sup>G93A</sup> animals.

GENES REGULATED BY G-CSF IN MOTONEURONS OF SOD1<sup>G93A</sup> MICE

Next, we focused on the identity of differentially expressed genes in motoneurons after G-CSF treatment. G-CSF had an effect on the regulation of 463 transcripts in total, and was able to adjust the expression of 187 transcripts that were either

up- or downregulated in SOD1<sup>G93A</sup> mice compared to controls (Tables 3, 4, Figure 3).

Two of the most downregulated transcripts in SOD1<sup>G93A</sup> mice, kcnq5 and larp4, had their expression enhanced by a factor 3 upon G-CSF treatment as compared to placebo group. Larp4 has no known specific functions in neurons but is known to regulates mRNA stability (Yang et al., 2011). G-CSF prevented the upregulation of NPY, a well-known neuropeptide involved in energetic metabolism and food intake but also expressed by neurons in response to stress such as spinal cord inflammation (Taylor et al., 2014). Upregulation of cklf, a cytokine with lymphocyte chemoattractive properties, was also inhibited by treatment (Han et al., 2001). Interestingly, Cpt1a and elovl4, two transcripts directly associated with lipid metabolism, were downregulated in SOD1<sup>G93A</sup> motoneurons, and upregulated by G-CSF. Cpt1a controls the entrance of fatty acids into mitochondria and is therefore the rate limiting enzyme for fatty acid oxidation, and Elovl4 elongates long fatty acids without distinction on their saturation levels, suggesting that G-CSF could enhance the processing of fatty acids in spinal cord. Additionally, G-CSF modulated the expression of 251 previously unaltered transcripts. The effect of G-CSF was prominent on gene ontology biological functions deregulated in SOD1<sup>G93A</sup> mice, including metabolic process, and response to stimulus (Supplementary Figure 2).

In order to assess the reliability of the analysis, we selected several genes for validation by quantitative real time PCR (qPCR). Those genes were selected because of their altered expression in SOD1<sup>G93A</sup> mice at 15 weeks of age and by a significant effect of G-CSF on their regulation. An additional criterion was their potential role in the neurodegenerative process of ALS based on the literature. We chose to study the regulation of CNOT10, CTSB, Prss12 and BMPR1a. CNOT10 is a member of the CCR4-NOT complex that acts as a chaperon platform (Collart and Panasenko, 2012). The CCR4-NOT complex controls and promotes neuromuscular junction formation in drosophila (Pradhan et al., 2012). Prss12, also known as motopsin, is a protease specific to neurons and its expression is downregulated upon stress, such as axotomy, in motoneurons (Numajiri et al., 2009). Its regulation is therefore seen as a marker of denervation. We chose ctsb, as meta-analysis of transcriptomic studies in ALS have shown its expression systematically altered, both in human patients and SOD1<sup>G93A</sup> mice (Saris et al., 2013). Ctsb codes for cathepsin B, a protein associated with neuroinflammation, motoneuron degeneration, and tissue destruction (Kikuchi et al., 2003). BMPR1a is a receptor serine/threonine kinase and regulates apoptosis and cell differentiation (Mishina et al., 2004). In the central nervous system, BMP signaling is required during formation of large-transmitting synapses (Xiao et al., 2013), build-up of neuromuscular junctions (Higashi-Kovtun et al., 2010) and is stimulated upon spinal cord and sciatic nerve injury (Tsujii et al., 2009; Sahni et al., 2010).

Relative expressions for these transcripts were highly similar to the relative abundance derived from microarray analyses (Figure 4). In particular, upregulation of CNOT10 present in SOD1<sup>G93A</sup> mice is completely reversed upon G-CSF treatment. The regulation of bmpr1a showed a trend to the wild type value but this was not significant in the qPCR analysis.

<figure><img src="image_4.png" /><figcaption>FIGURE 2 | Canonical plot shows dispersion of treated animals off the non-treated groups. Multivariate analyses have been conducted with genes having a standard deviation greater than 1.2 (734 transcripts). (A) Canonical plot shows that symptomatic SOD-1 mice (red dots) are distant from other groups, including asymptomatic SOD1 mice (blue dots) and when considered as unknown, SOD-1 mice treated with G-CSF (green dots) show a strong dispersion and tend to be close to wild types at 11 weeks (gray dots) or 15 weeks of age (black dots). (B) After G-CSF treatment, discriminant analysis classified only three samples of SOD1 motoneurons as symptomatic motor neurons. One sample was considered has presymptomatic and three as wild type motor neurons.</figcaption></figure>
<table><caption>Table 3 | List of the main transcripts whose expression is downregulated by G-CSF treatment in SOD1<sup>G93A</sup> mice at 15 week of age.</caption><thead><tr><th>Probe set ID</th><th>Gene symbol</th><th>p-value</th><th>WT</th><th>SOD1</th><th>SOD1 G-CSF</th><th>FC WTvsSOD1</th><th>FC WT vs. SOD1-GCSF</th></tr></thead><tbody><tr><td>1421979_at</td><td>Phex</td><td>0.000</td><td>1.994</td><td>4.091</td><td>1.983</td><td>4.279</td><td>0.992</td></tr><tr><td>1421316_at</td><td>Lce1g</td><td>0.001</td><td>2.505</td><td>4.563</td><td>2.415</td><td>4.167</td><td>0.940</td></tr><tr><td>1419127_at</td><td>Npy</td><td>0.000</td><td>4.025</td><td>6.074</td><td>4.075</td><td>4.138</td><td>1.035</td></tr><tr><td>1456084_x_at</td><td>Fmod</td><td>0.000</td><td>2.536</td><td>4.400</td><td>2.752</td><td>3.640</td><td>1.161</td></tr><tr><td>1437720_at</td><td>Eif2d</td><td>0.000</td><td>4.914</td><td>6.458</td><td>5.661</td><td>2.916</td><td>1.679</td></tr><tr><td>1438856_x_at</td><td>Serpinb5</td><td>0.018</td><td>2.880</td><td>4.356</td><td>2.842</td><td>2.781</td><td>0.974</td></tr><tr><td>1441836_x_at</td><td>Gtsf1</td><td>0.003</td><td>2.524</td><td>3.958</td><td>2.718</td><td>2.703</td><td>1.144</td></tr><tr><td>1444780_at</td><td>Nav2</td><td>0.001</td><td>3.251</td><td>4.684</td><td>3.440</td><td>2.700</td><td>1.140</td></tr><tr><td>1418778_at</td><td>Ccdc109b</td><td>0.001</td><td>3.279</td><td>4.642</td><td>3.071</td><td>2.571</td><td>0.865</td></tr><tr><td>1437775_at</td><td>Dlst</td><td>0.003</td><td>4.229</td><td>5.570</td><td>4.142</td><td>2.532</td><td>0.941</td></tr><tr><td>1448048_at</td><td>Nmrk1</td><td>0.006</td><td>2.889</td><td>4.219</td><td>2.927</td><td>2.515</td><td>1.027</td></tr><tr><td>1458056_at</td><td>Srek1</td><td>0.003</td><td>2.954</td><td>4.280</td><td>3.260</td><td>2.507</td><td>1.236</td></tr><tr><td>1454557_at</td><td>6720454L07Rik</td><td>0.000</td><td>3.338</td><td>4.651</td><td>3.503</td><td>2.485</td><td>1.122</td></tr><tr><td>1417491_at</td><td>Ctsb</td><td>0.000</td><td>6.668</td><td>7.909</td><td>7.128</td><td>2.364</td><td>1.376</td></tr><tr><td>1436242_a_at</td><td>Cklf</td><td>0.019</td><td>6.877</td><td>8.107</td><td>6.819</td><td>2.346</td><td>0.961</td></tr><tr><td>1458335_x_at</td><td>Urm1</td><td>0.000</td><td>6.757</td><td>7.911</td><td>7.240</td><td>2.225</td><td>1.398</td></tr><tr><td>1428514_at</td><td>Cpne3</td><td>0.000</td><td>2.936</td><td>4.085</td><td>3.124</td><td>2.218</td><td>1.139</td></tr><tr><td>1456816_at</td><td>Zmiz1</td><td>0.001</td><td>3.413</td><td>4.530</td><td>3.732</td><td>2.170</td><td>1.248</td></tr><tr><td>1433733_a_at</td><td>Cry1</td><td>0.001</td><td>2.270</td><td>3.378</td><td>2.549</td><td>2.155</td><td>1.213</td></tr><tr><td>1418467_at</td><td>Smarcd3</td><td>0.000</td><td>4.533</td><td>5.617</td><td>5.038</td><td>2.121</td><td>1.419</td></tr><tr><td>1429688_at</td><td>Arntl2</td><td>0.014</td><td>2.893</td><td>3.971</td><td>2.846</td><td>2.111</td><td>0.968</td></tr><tr><td>1460578_at</td><td>Fgd5</td><td>0.000</td><td>3.593</td><td>4.669</td><td>3.192</td><td>2.108</td><td>0.757</td></tr><tr><td>1457147_at</td><td>Etl4</td><td>0.001</td><td>3.500</td><td>4.570</td><td>3.118</td><td>2.099</td><td>0.767</td></tr><tr><td>1437822_at</td><td>Yme1l1</td><td>0.004</td><td>3.426</td><td>4.485</td><td>3.395</td><td>2.084</td><td>0.979</td></tr><tr><td>1458669_at</td><td>Agpat5</td><td>0.003</td><td>3.253</td><td>4.304</td><td>3.132</td><td>2.073</td><td>0.919</td></tr><tr><td>1443883_at</td><td>Sys1</td><td>0.000</td><td>4.305</td><td>5.316</td><td>4.112</td><td>2.015</td><td>0.875</td></tr><tr><td>1433882_at</td><td>Cnot10</td><td>0.000</td><td>5.577</td><td>6.566</td><td>5.363</td><td>1.985</td><td>0.862</td></tr><tr><td>1437718_x_at</td><td>Fmod</td><td>0.000</td><td>2.890</td><td>3.862</td><td>2.837</td><td>1.961</td><td>0.964</td></tr><tr><td>1440398_at</td><td></td><td>0.001</td><td>2.641</td><td>3.608</td><td>2.866</td><td>1.954</td><td>1.169</td></tr><tr><td>1445861_at</td><td>Usp25</td><td>0.014</td><td>3.428</td><td>4.393</td><td>3.399</td><td>1.952</td><td>0.980</td></tr><tr><td>1443242_at</td><td>D5Ertd121e</td><td>0.005</td><td>2.186</td><td>3.150</td><td>2.205</td><td>1.950</td><td>1.013</td></tr></tbody></table>

Taken together, these data demonstrate that G-CSF promotes transcriptional changes in motoneurons in symptomatic SOD1<sup>G93A</sup> mice and can correct the alterations of many disease related gene expression changes.

## DISCUSSION

ALS is a fatal condition characterized by the constant death of motoneurons. Although the triggering factors for ALS remain unknown, the step that induces symptoms is the disruption of the motor units, with the dismantlement of neuromuscular junctions and motoneuron stress and death. The study of global changes in the transcriptome of specific cells in the motor axis is a promising approach to address the complexity of ALS (Henriques and Gonzalez De Aguilar, 2011).

## TRANSCRIPTOMIC CHANGES OF MOTONEURONS FROM SOD1<sup>G93A</sup> MICE

Here, we studied the transcriptome of motoneurons of SOD1<sup>G93A</sup> mice, at two disease stages and after a pharmacological treatment with G-CSF, a trophic factor with neuroprotective properties

(Pitzer et al., 2008; Henriques et al., 2010a,b, 2011). Our aim was to characterize the transcriptomic response of motor neurons to a neuroprotective treatment based on G-CSF.

In 2005, a first study was published describing the transcriptome of laser microdissected motoneurons of SOD1<sup>G93A</sup> mice at three different disease stages (Perrin et al., 2005). The conclusion of the authors was that few changes occurred at presymptomatic (28 deregulated transcripts) and early symptomatic (148 deregulated transcripts) stages, and that symptomatic stages were characterized by few expression changes in cell-death-associated genes. They found deregulation of related to energetic metabolism, cell growth and/or maintenance (Perrin et al., 2005). Many of these genes also showed an altered expression in our study. In 2007, another group published data on the transcriptome of dissected motoneurons from SOD1<sup>G93A</sup> mice, and reported dysregulations on genes related to transcription and nuclear proteins, protein synthesis, metabolism and cell cycle/growth (Ferraiuolo et al., 2007). Both studies identified lower deregulation events than our study. This difference can most likely be explained by higher statistical power in our study
<table><caption>Table 4 | List of the main transcripts whose expression is upregulated by G-CSF treatment in SOD1<sup>G93A</sup> mice at 15 week of age.</caption><thead><tr><th>Probe set ID</th><th>Gene symbol</th><th>p-value</th><th>WT</th><th>SOD1</th><th>SOD1 G-CSF</th><th>FC WTvsSOD1</th><th>FC WT vs. SOD1-GCSF</th></tr></thead><tbody><tr><td>1457587_at</td><td>Kcnq5</td><td>0.000</td><td>7.713</td><td>4.026</td><td>5.814</td><td>0.078</td><td>0.268</td></tr><tr><td>1455102_at</td><td>Larp4</td><td>0.000</td><td>6.999</td><td>4.642</td><td>6.371</td><td>0.195</td><td>0.647</td></tr><tr><td>1455381_at</td><td>Marf1</td><td>0.002</td><td>7.906</td><td>5.585</td><td>7.319</td><td>0.200</td><td>0.666</td></tr><tr><td>1449861_at</td><td>Nek4</td><td>0.001</td><td>5.754</td><td>3.478</td><td>5.055</td><td>0.207</td><td>0.616</td></tr><tr><td>1460440_at</td><td>Lphn3</td><td>0.000</td><td>6.848</td><td>4.670</td><td>6.104</td><td>0.221</td><td>0.597</td></tr><tr><td>1419555_at</td><td>Elf5</td><td>0.000</td><td>5.646</td><td>3.526</td><td>5.185</td><td>0.230</td><td>0.727</td></tr><tr><td>1455960_at</td><td>Megf9</td><td>0.000</td><td>6.959</td><td>4.882</td><td>6.446</td><td>0.237</td><td>0.701</td></tr><tr><td>1429475_at</td><td>Ubash3b</td><td>0.001</td><td>6.351</td><td>4.279</td><td>5.680</td><td>0.238</td><td>0.628</td></tr><tr><td>1420388_at</td><td>Prss12</td><td>0.000</td><td>7.504</td><td>5.467</td><td>6.385</td><td>0.244</td><td>0.461</td></tr><tr><td>1449348_at</td><td>Mpp6</td><td>0.000</td><td>8.104</td><td>6.114</td><td>7.248</td><td>0.252</td><td>0.552</td></tr><tr><td>1438435_at</td><td>Acer3</td><td>0.001</td><td>4.492</td><td>2.529</td><td>4.626</td><td>0.257</td><td>1.098</td></tr><tr><td>1429107_at</td><td>Ubr3</td><td>0.001</td><td>6.415</td><td>4.532</td><td>5.996</td><td>0.271</td><td>0.748</td></tr><tr><td>1451217_a_at</td><td>Immp1l</td><td>0.000</td><td>6.378</td><td>4.530</td><td>6.677</td><td>0.278</td><td>1.230</td></tr><tr><td>1419357_at</td><td>Isy1</td><td>0.000</td><td>5.599</td><td>3.778</td><td>5.114</td><td>0.283</td><td>0.715</td></tr><tr><td>1423597_at</td><td>Atp8a1</td><td>0.000</td><td>8.998</td><td>7.187</td><td>8.151</td><td>0.285</td><td>0.556</td></tr><tr><td>1418973_at</td><td>Blzf1</td><td>0.000</td><td>6.188</td><td>4.426</td><td>5.637</td><td>0.295</td><td>0.683</td></tr><tr><td>1454937_at</td><td>B630005N14Rik</td><td>0.000</td><td>6.420</td><td>4.661</td><td>5.956</td><td>0.295</td><td>0.725</td></tr><tr><td>1435224_at</td><td>Crebbp</td><td>0.015</td><td>7.833</td><td>6.076</td><td>7.621</td><td>0.296</td><td>0.863</td></tr><tr><td>1455728_at</td><td>Pten</td><td>0.001</td><td>5.930</td><td>4.204</td><td>5.737</td><td>0.302</td><td>0.875</td></tr><tr><td>1416468_at</td><td>Aldh1a1</td><td>0.000</td><td>8.034</td><td>6.314</td><td>7.206</td><td>0.304</td><td>0.563</td></tr><tr><td>1425493_at</td><td>Bmpr1a</td><td>0.000</td><td>4.559</td><td>2.907</td><td>4.253</td><td>0.318</td><td>0.809</td></tr><tr><td>1417564_at</td><td>Med7</td><td>0.000</td><td>7.151</td><td>5.551</td><td>6.891</td><td>0.330</td><td>0.835</td></tr><tr><td>1423298_at</td><td>Add3</td><td>0.000</td><td>7.243</td><td>5.681</td><td>6.749</td><td>0.339</td><td>0.710</td></tr><tr><td>1421904_at</td><td>Tgs1</td><td>0.001</td><td>8.757</td><td>7.283</td><td>8.236</td><td>0.360</td><td>0.697</td></tr><tr><td>1416655_at</td><td>C1galt1c1</td><td>0.001</td><td>9.662</td><td>8.190</td><td>9.617</td><td>0.360</td><td>0.969</td></tr><tr><td>1460409_at</td><td>Cpt1a</td><td>0.000</td><td>8.140</td><td>6.686</td><td>7.440</td><td>0.365</td><td>0.616</td></tr><tr><td>1451308_at</td><td>Elovl4</td><td>0.000</td><td>6.986</td><td>5.539</td><td>6.436</td><td>0.367</td><td>0.683</td></tr><tr><td>1423490_at</td><td>Fbxo3</td><td>0.000</td><td>8.932</td><td>7.495</td><td>8.538</td><td>0.369</td><td>0.761</td></tr><tr><td>1455807_at</td><td>Tspyl5</td><td>0.001</td><td>7.369</td><td>5.965</td><td>6.909</td><td>0.378</td><td>0.727</td></tr></tbody></table>

as both studies used fewer animals (n = 3/group) than our study (at least n = 6 per group). Higher degree of contaminating RNA from motoneuron-surrounding cells could also account for the differences between studies.

Based on our results, the transcriptome of isolated motoneurons of SOD1<sup>G93A</sup> mice is characterized by the deregulation of genes, starting at presymptomatic disease stages, related to response to stress, to regulation of transcription and to metabolic processes.

Response to stress is stimulated early on, at presymptomatic disease stage. Trophic signaling seems also to be stimulated, with the upregulation of the trophic factors such as fgf21 and vgf at presymptomatic and symptomatic disease stages. FGF21 belongs to the fibroblast growth factor family, known to counteract motoneuron death after axotomy-induced apoptosis (Cuevas et al., 1995), but also shown to regulate activation of astrocytes in animal models of ALS (Cassina et al., 2005). Additionally, motoneurons of SOD1<sup>G93A</sup> show an upregulation of transcription factors related to stress response, that are also closely connected to trophic factors, such as ATF3 (Averill et al., 2004).

During laser-microdissection of motoneurons, small amount of astrocyte RNA may have been captured. Although this

contamination is minimal in our study, genes that are expressed at low abundance in motoneurons but at high levels in glial cells could be detected as differentially regulated, provided that the difference between two groups is very strong. An indication for that could be the detection of GFAP among the deregulated transcripts in our study (upregulated in SOD1 mice). Presence of transcripts of glial origin in studies working on laser-captured motoneurons has been described previously (Perrin et al., 2005; Ferraiuolo et al., 2007). Therefore, it is possible that the alterations in the transcriptome linked to energetic/lipid metabolism arise from astrocytes surrounding motoneurons. G-CSF was indeed interestingly able to stimulate the expression of cpt1a in SOD1<sup>G93A</sup> mice, which is surprising since fatty acids are not the main energetic source for neurons. It is therefore possible that upregulation of cpt1a takes place in astrocytes and not in motoneurons, presumably to sustain higher energetic needs for neurons. Indeed, during energetic stress, astrocytes are able to convert lipid to ketone bodies that will be used as energetic fuel by surrounding neurons (Guzman and Blazquez, 2004). Interestingly, recent studies have shown increased ketone bodies in the cerebrospinal fluid of patients and fatty acid and ketone body supplementations have
<figure><img src="image_5.png" /></figure>

**FIGURE 3** | G-CSF treatment modulates transgene-related gene expression changes. (A) Chart showing the number of differentially regulated genes after G-CSF treatment in the motoneurons of SOD1 mice. (B–C) Given is the relative expression of the 30 most significantly modified transcripts re-adjusted by G-CSF treatment and previously found downregulated (B) or upregulated (C) in symptomatic SOD1 mice.

neuroprotective effects on motoneurons in animal models of ALS (Blasco et al., 2010; Schmitt et al., 2014). Energetic metabolism represents one promising field of research in ALS. Patients present with hypermetabolism, higher energetic needs and a higher level of circulating lipids is associated with better survival (Schmitt et al., 2014). Two recent clinical reports strongly suggest benefit after high lipid/caloric diets in ALS patients with a strong effect on survival (Dorst et al., 2013; Wills et al., 2014).

We found deregulation of several genes associated to RNA
metabolism, such as Wars, Eef1a1 and Eif4e. Growing evidence
suggests a role for defect RNA processing in motoneuron death.
ALS patients present pathological inclusion of TDP-43 and FUS,
two proteins that closely interact with RNA metabolism (Ling
et al., 2013). Mutations in TARDBP (coding for TDP-43) and
FUS are described in ALS patients with a family history but also
in patients with fronto-temporal dementia (FTD) with a high
risk (around 30%) to develop ALS. Expanded repeats in the gene
C9orf72 are largely present in sporadic and familial ALS and
FTD cases. Interestingly, whose expanded repeats are believed to
induce RNA toxicity through a novel negative gain-of-function.
Our results therefore provide interesting new clues when related
to most recent pathophysiological hypotheses in ALS research.

<figure><img src="image_6.png" /></figure>

**FIGURE 4** | Validation of the expression of selected transcripts by qPCR. (A–D) Expression levels from microarray analyses of Bmpr1a (A), CNOT10 (B), CTSB (C) and PRSS12 (D) at the symptomatic disease stage. (E–H) Expression levels derived from quantitative PCR of Bmpr1a (E), CNOT10 (F), CTSB (G) and PRSS12 (H) at the symptomatic disease stage. Data are presented as mean ± standard variation of the mean. \*p < 0.05 (ANOVA, followed by multicomparison *t*-test).

G-CSF ACTION ON MOTONEURONS

Our transcriptomic study describes for the first time the tran-
scriptomic changes in motoneurons after a neuroprotective/ -
regenerative treatment. We chose to perform our study with
SOD1G93A mice at 15 week of age with an administration of
G-CSF for 4 weeks. At this age, motor symptom are clearly
established in untreated animals (Pitzer et al., 2008; Henriques
et al., 2011) and G-CSF has a significant beneficial effect on
motor functions and preserves motoneuron counts (Pitzer et al.,
2008). In our study, the altered transcriptome of SOD1G93A alpha
motoneurons has been modified by G-CSF and the expressions
of 187 genes were partially restored. Genes involved in the maintenance of the motor-units, such as Prss12 and CNOT10, were modulated upon treatment, but effect of G-CSF was also observed on other different molecular functions such as neuroinflammation, RNA processing, and energetic metabolism. Completely unsupervised analyses identified more than half of treated mice as asymptomatic based on their transcriptomic profile, confirming the positive effect of the treatment on disease progression.

Modification of the transcript was most likely due to a direct effect of G-CSF on the motor neurons. Motoneurons express the receptor of G-CSF in basal condition and can upregulate it after an acute stress in neonatal age (Pitzer et al., 2008; Henriques et al., 2010a). We have also demonstrated that G-CSF can rescue motoneurons in a pure apoptotic condition, induced by sciatic nerve axotomy (Henriques et al., 2010a) and promotes reinnervation after sciatic nerve crush (Henriques et al., 2011). However, we cannot formally conclude whether all or most of the observed changes are due to direct actions of G-CSF on motoneurons or reflect improved health status of the alpha motoneurons by indirect treatment effects. Indeed, G-CSF was first described as a hematopoietic growth factor involved in the proliferation and differentiation of neutrophils. Despite an abundant literature on the involvement of neuroinflammation in ALS (Philips and Robberecht, 2011), there is no evidence that neutrophils have a role in the initiation or progression of ALS, neither in animal model nor in ALS patients. Additionally, a recent report suggests that G-CSF might be involved in response to muscle damage, which is also a characteristic of ALS (Kruger et al., 2014).

## CONCLUSION

G-CSF was able to re-adjust the expression of genes important for motoneuron survival and function but also genes from other physiological pathways. The global changes occurring after treatment strengthens the notion of multimodal activity of G-CSF for preservation of motor units in ALS. Our results complement the existing data on characterization of G-CSF's action in mouse ALS models (Pitzer et al., 2008; Henriques et al., 2010a, 2011; Naumenko et al., 2011; Pollari et al., 2011) on the gene expression level. This increases confidence in G-CSF as a promising drug candidate for ALS.

## AUTHOR CONTRIBUTIONS

Armin Schneider designed the study. Stefan Kastner, Eva Chatzikonstantinou, Friederike Kirsch, Carola Krüger, Christian Plaas performed animal experiments, histological work, laser microdissection, RNA preparation and amplification. Norbert Gretz performed array hybridizations. Alexandre Henriques performed data analysis. Armin Schneider, Alexandre Henriques, Stefan Kastner, Claudia Pitzer, Carola Krüger, Robert Spoelgen, Jose-Luis Gonzalez De Aguilar and Armin Schneider were involved in data evaluation and discussions. Alexandre Henriques and Armin Schneider wrote the manuscript. All authors read and approved the final manuscript.

## ACKNOWLEDGMENTS

The authors would like to acknowledge the excellent technical assistance of Ulrike Bolz, Frank Herzog, and Gisela Eisenhardt.

This work was supported by funds from European Community's Health Seventh Framework Program under grant agreement no. 259867 (FP7/2007-2013). Alexandre Henriques is a research fellow receiving funds from FP7/2007-2013. Jose-Luis Gonzalez De Aguilar is recipient of a "Chaire d'Excellence INSERM/Université de Strasbourg."

## SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at: http://www.frontiersin.org/journal/10.3389/fncel.2014.00464/abstract

**Supplementary Figure 1 | Laser microdissection and neuronal specificity.**
(A) Representative pictures of spinal cord section before and after laser capture. Sections are stained with thionin. (B) Given is the relative expression of genes that are accepted as specific markers of astrocytes (n = 6 markers), neurons (n = 7 markers), oligodendrocytes (n = 7), microglia (n = 4) and endothelial cells (n = 4).

**Supplementary Figure 2 | Biological and molecular functions of genes altered in SOD1G93A mice.** Given are the main biological (A) and molecular (B) functions of genes altered in SOD1G93A motoneurons at presymptomatic and symptomatic disease stages, and re-adjusted after G-CSF treatment.

**Availability of Supporting Data |** Microarray data files have been uploaded to the Gene Expression Omnibus (GEO) database repository, under the accession number GSE60856 (www.ncbi.nlm.nih.gov/geo/).

## REFERENCES

Averill, S., Michael, G. J., Shortland, P. J., Leavesley, R. C., King, V. R., and Bradbury, E. J. (2004). NGF and GDNF ameliorate the increase in ATF3 expression which occurs in dorsal root ganglion cells in response to peripheral nerve injury. *Eur. J. Neurosci.* 19, 1437-1445. doi: 10.1111/j.1460-9568.2004.03241.x

Bandyopadhyay, U., Cotney, J., Nagy, M., Oh, S., Leng, J., Mahajan, M., et al. (2013). RNA-Seq profiling of spinal cord motor neurons from a presymptomatic SOD1 ALS mouse. *PLoS ONE* 8:e53575. doi: 10.1371/journal.pone.0053575

Blasco, H., Corcia, P., Moreau, C., Veau, S., Fournier, C., Vourc'h, P., et al. (2010). 1H-NMR-based metabolomic profiling of CSF in early amyotrophic lateral sclerosis. *PLoS ONE* 5:e13223. doi: 10.1371/journal.pone.0013223

Bonilla, I. E., Tanabe, K., and Strittmatter, S. M. (2002). Small proline-rich repeat protein 1A is expressed by axotomized neurons and promotes axonal outgrowth. *J. Neurosci.* 22, 1303-1315.

Cahoy, J. D., Emery, B., Kaushal, A., Foo, L. C., Zamanian, J. L., Christopherson, K. S., et al. (2008). A transcriptome database for astrocytes, neurons, and oligodendrocytes: a new resource for understanding brain development and function. *J. Neurosci.* 28, 264-278. doi: 10.1523/JNEUROSCI.4178-07.2008

Cassina, P., Pehar, M., Vargas, M. R., Castellanos, R., Barbeito, A. G., Estevez, A. G., et al. (2005). Astrocyte activation by fibroblast growth factor-1 and motor neuron apoptosis: implications for amyotrophic lateral sclerosis. *J. Neurochem.* 93, 38-46. doi: 10.1111/j.1471-4159.2004.02984.x

Chen, L. J., Ren, Y. H., Liu, L., Zhang, X. Q., Zhao, Y., Wu, W. T., et al. (2010). Upregulated expression of GAP-43 mRNA and protein in anterior horn motoneurons of the spinal cord after brachial plexus injury. *Arch. Med. Res.* 41, 513-518. doi: 10.1016/j.arcmed.2010.10.005

Chio, A., Mora, G., La Bella, V., Caponnetto, C., Mancardi, G., Sabatelli, M., et al. (2011). Repeated courses of granulocyte colony-stimulating factor in amyotrophic lateral sclerosis: clinical and biological results from a prospective multicenter study. *Muscle Nerve* 43, 189-195. doi: 10.1002/mus.21851

Collart, M. A., and Panasenko, O. O. (2012). The Ccr4-not complex. *Gene* 492, 42-53. doi: 10.1016/j.gene.2011.09.033
Cuevas, P., Carceller, F., and Gimenez-Gallego, G. (1995). Acidic fibroblast growth factor prevents post-axotomy neuronal death of the newborn rat facial nerve. *Neurosci. Lett.* 197, 183-186. doi: 10.1016/0304-3940(95)11926-N

Cui, G., Wang, Y., Yu, S., Yang, L., Li, B., Wang, W., et al. (2014). The expression changes of vacuolar protein sorting 4B (VPS4B) following middle cerebral artery occlusion (MCAO) in adult rats brain hippocampus. *Cell. Mol. Neurobiol.* 34, 83-94. doi: 10.1007/s10571-013-9989-5

Dorst, J., Cypionka, J., and Ludolph, A. C. (2013). High-caloric food supplements in the treatment of amyotrophic lateral sclerosis: a prospective interventional study. *Amyotroph. Lateral Scler. Frontotemporal Degener.* 14, 533-536. doi: 10.3109/21678421.2013.823999

Duning, T., Schiffbauer, H., Warnecke, T., Mohammadi, S., Floel, A., Kolpatzik, K., et al. (2011). G-CSF prevents the progression of structural disintegration of white matter tracts in amyotrophic lateral sclerosis: a pilot trial. *PLoS ONE* 6:e17770. doi: 10.1371/journal.pone.0017770

Ferraiuolo, L., Heath, P. R., Holden, H., Kasher, P., Kirby, J., and Shaw, P. J. (2007). Microarray analysis of the cellular pathways involved in the adaptation to and progression of motor neuron injury in the SOD1G93A mouse model of familial ALS. *J. Neurosci.* 27, 9201-9219. doi: 10.1523/JNEUROSCI.1470-07.2007

Frey, D., Schneider, C., Xu, L., Borg, J., Spooren, W., and Caroni, P. (2000). Early and selective loss of neuromuscular synapse subtypes with low sprouting competence in motoneuron diseases. *J. Neurosci.* 20, 2534-2542.

Gordon, P. H. (2013). Amyotrophic lateral sclerosis: an update for 2013 clinical features, pathophysiology, management and therapeutic trials. *Aging Disease* 4, 295-310. doi: 10.14336/AD.2013.0400295

Grasselli, G., Mandolesi, G., Strata, P., and Cesare, P. (2011). Impaired sprouting and axonal atrophy in cerebellar climbing fibres following in vivo silencing of the growth-associated protein GAP-43. *PLoS ONE* 6:e20791. doi: 10.1371/journal.pone.0020791

Gurney, M. E., Pu, H., Chiu, A. Y., Dal Canto, M. C., Polchow, C. Y., Alexander, D. D., et al. (1994). Motor neuron degeneration in mice that express a human Cu,Zn superoxide dismutase mutation. *Science* 264, 1772-1775. doi: 10.1126/science.8209258

Guzman, M., and Blazquez, C. (2004). Ketone body synthesis in the brain: possible neuroprotective effects. *Prostaglandins Leukot. Essent. Fatty Acids* 70, 287-292. doi: 10.1016/j.plefa.2003.05.001

Han, J. H., Ryu, H. H., Jun, M. H., Jang, D. J., and Lee, J. A. (2012). The functional analysis of the CHMP2B missense mutation associated with neurodegenerative diseases in the endo-lysosomal pathway. *Biochem. Biophys. Res. Commun.* 421, 544-549. doi: 10.1016/j.bbrc.2012.04.041

Han, W., Lou, Y., Tang, J., Zhang, Y., Chen, Y., Li, Y., et al. (2001). Molecular cloning and characterization of chemokine-like factor 1 (CKLF1), a novel human cytokine with unique structure and potential chemotactic activity. *Biochem. J.* 357(Pt 1), 127-135. doi: 10.1042/0264-6021:35:70127

Harr, B., and Schlotterer, C. (2006). Comparison of algorithms for the analysis of Affymetrix microarray data as evaluated by co-expression of genes in known operons. *Nucleic Acids Res.* 34:e8. doi: 10.1093/nar/gnj010

Henriques, A., and Gonzalez De Aguilar, J. L. (2011). Can transcriptomics cut the gordian knot of amyotrophic lateral sclerosis? *Curr. Genomics* 12, 506-515. doi: 10.2174/138920211797904043

Henriques, A., Pitzer, C., Dittgen, T., Klugmann, M., Dupuis, L., and Schneider, A. (2011). CNS-targeted viral delivery of G-CSF in an animal model for ALS: improved efficacy and preservation of the neuromuscular unit. *Mol. Therapy* 19, 284-292. doi: 10.1038/mt.2010.271

Henriques, A., Pitzer, C., Dupuis, L., and Schneider, A. (2010a). G-CSF protects motoneurons against axotomy-induced apoptotic death in neonatal mice. *BMC Neurosci.* 11:25. doi: 10.1186/1471-222010-11-25

Henriques, A., Pitzer, C., and Schneider, A. (2010b). Neurotrophic growth factors for the treatment of amyotrophic lateral sclerosis: where do we stand? *Front. Neurosci.* 4:32. doi: 10.3389/fnins.2010.00032

Higashi-Kovtun, M. E., Mosca, T. J., Dickman, D. K., Meinertzhagen, I. A., and Schwarz, T. L. (2010). Importin-beta11 regulates synaptic

phosphorylated mothers against decapentaplegic, and thereby influences synaptic development and function at the Drosophila neuromuscular junction. *J. Neurosci.* 30, 5253-5268. doi: 10.1523/JNEUROSCI.3739-09.2010

Huang, H., and Trussell, L. O. (2011). KCNQ5 channels control resting properties and release probability of a synapse. *Nat. Neurosci.* 14, 840-847. doi: 10.1038/nn.2830

Hyatt Sachs, H., Schreiber, R. C., Shoemaker, S. E., Sabe, A., Reed, E., and Zigmond, R. E. (2007). Activating transcription factor 3 induction in sympathetic neurons after axotomy: response to decreased neurotrophin availability. *Neuroscience* 150, 887-897. doi: 10.1016/j.neuroscience.2007.10.008

Kikuchi, H., Yamada, T., Furuya, H., Doh-ura, K., Ohyagi, Y., Iwaki, T., et al. (2003). Involvement of cathepsin B in the motor neuron degeneration of amyotrophic lateral sclerosis. *Acta Neuropathol.* 105, 462-468. doi: 10.1007/s00401-002-0667-9

Kobbe, G., Bruns, I., Fenk, R., Czibere, A., and Haas, R. (2009). Pegfilgrastim for PBSC mobilization and autologous haematopoietic SCT. *Bone Marrow Transplant.* 43, 669-677. doi: 10.1038/bmt.2009.59

Kruger, K., Pilat, C., Schild, M., Lindner, N., Frech, T., Muders, K., et al. (2014). Progenitor cell mobilization after exercise is related to systemic levels of G-CSF and muscle damage. *Scand. J. Med. Sci. Sports* doi: 10.1111/sms.12320

Li, S., and Strittmatter, S. M. (2003). Delayed systemic Nogo-66 receptor antagonist promotes recovery from spinal cord injury. *J. Neurosci.* 23, 4219-4227.

Ling, S. C., Polymenidou, M., and Cleveland, D. W. (2013). Converging mechanisms in ALS and FTD: disrupted RNA and protein homeostasis. *Neuron* 79, 416-438. doi: 10.1016/j.neuron.2013.07.033

Martinez, H. R., Gonzalez-Garza, M. T., Moreno-Cuevas, J. E., Caro, E., Gutierrez-Jimenez, E., and Segura, J. J. (2009). Stem-cell transplantation into the frontal motor cortex in amyotrophic lateral sclerosis patients. *Cytotherapy* 11, 26-34. doi: 10.1080/14653240802644651

Mishina, Y., Starbuck, M. W., Gentile, M. A., Fukuda, T., Kasparcova, V., Seedor, J. G., et al. (2004). Bone morphogenetic protein type IA receptor signaling regulates postnatal osteoblast function and bone remodeling. *J. Biol. Chem.* 279, 27560-27566. doi: 10.1074/jbc.M404222200

Motomura, W., Yoshizaki, T., Ohtani, K., Okumura, T., Fukuda, M., Fukuzawa, J., et al. (2008). Immunolocalization of a novel collectin CL-K1 in murine tissues. *J. Histochem. Cytochem.* 56, 243-252. doi: 10.1369/jhc.7A7312.2007

Naumenko, N., Pollari, E., Kurronen, A., Giniatullina, R., Shakirzyanova, A., Magga, J., et al. (2011). Gender-Specific mechanism of synaptic impairment and its prevention by GCSF in a mouse model of ALS. *Front. Cellular Neurosci.* 5:26. doi: 10.3389/fncel.2011.00026

Nefussy, B., Artamonov, I., Deutsch, V., Naparstek, E., Nagler, A., and Drory, V. E. (2010). Recombinant human granulocyte-colony stimulating factor administration for treating amyotrophic lateral sclerosis: a pilot study. *Amyotroph. Lateral Scler.* 11, 187-193. doi: 10.3109/17482960902933809

Numajiri, T., Mitsui, S., Hisa, Y., Ishida, T., Nishino, K., and Yamaguchi, N. (2009). The expression of a motoneuron-specific serine protease, motopsin (PRSS12), after facial nerve axotomy in mice. *J. Plast. Reconstr. Aesthet. Surg.* 59, 393-397. doi: 10.1016/j.jbps.2005.04.018

Ohtani, K., Suzuki, Y., and Wakamiya, N. (2012). Biological functions of the novel collectins CL-L1, CL-K1, and CL-P1. *J. Biomed. Biotechnol.* 2012:493945. doi: 10.1155/2012/493945

Perrin, F. E., Boisset, G., Docquier, M., Schaad, O., Descombes, P., and Kato, A. C. (2005). No widespread induction of cell death genes occurs in pure motoneurons in an amyotrophic lateral sclerosis mouse model. *Hum. Mol. Genet.* 14, 3309-3320. doi: 10.1093/hmg/ddi357

Philips, T., and Robberecht, W. (2011). Neuroinflammation in amyotrophic lateral sclerosis: role of glial activation in motor neuron disease. *Lancet Neurol.* 10, 253-263. doi: 10.1016/S1474-4422(11)70015-1

Pitzer, C., Kruger, C., Plaas, C., Kirsch, F., Dittgen, T., Muller, R., et al. (2008). Granulocyte-colony stimulating factor improves outcome in a mouse model of amyotrophic lateral sclerosis. *Brain* 131(Pt 12), 3335-3347. doi: 10.1093/brain/awn243

Pollari, E., Savchenko, E., Jaronen, M., Kanninen, K., Malm, T., Wojciechowski, S., et al. (2011). Granulocyte colony stimulating factor attenuates inflammation in a mouse model of amyotrophic lateral sclerosis. *J. Neuroinflammation* 8:74. doi: 10.1186/1742-2094-8-74
Pradhan, S. J., Nesler, K. R., Rosen, S. F., Kato, Y., Nakamura, A., Ramaswami, M., et al. (2012). The conserved P body component HPat/Pat1 negatively regulates synaptic terminal growth at the larval Drosophila neuromuscular junction. *J. Cell Sci.* 125(Pt 24), 6105–6116. doi: 10.1242/jcs.113043

Sahni, V., Mukhopadhyay, A., Tysseling, V., Hebert, A., Birch, D., McGuire, T. L., et al. (2010). BMPR1a and BMPR1b signaling exert opposing effects on gliosis after spinal cord injury. *J. Neurosci.* 30, 1839–1855. doi: 10.1523/JNEUROSCI.4459-09.2010

Saris, C. G., Groen, E. J., Koekkoek, J. A., Veldink, J. H., and van den Berg, L. H. (2013). Meta-analysis of gene expression profiling in amyotrophic lateral sclerosis: a comparison between transgenic mouse models and human patients. *Amyotroph. Lateral Scler. Frontotemporal Degener.* 14, 177–189. doi: 10.3109/21678421.2012.729842

Schmitt, F., Hussain, G., Dupuis, L., Loeffler, J. P., and Henriques, A. (2014). A plural role for lipids in motor neuron diseases: energy, signaling and structure. *Front. Cell. Neurosci.* 8:25. doi: 10.3389/fncel.2014.00025

Sleigh, J. N., Barreiro-Iglesias, A., Oliver, P. L., Biba, A., Becker, T., Davies, K. E., et al. (2014). Chondrolectin affects cell survival and neuronal outgrowth in in vitro and in vivo models of spinal muscular atrophy. *Hum. Mol. Genet.*, 23, 855–869. doi: 10.1093/hmg/ddt477

Starkey, M. L., Davies, M., Yip, P. K., Carter, L. M., Wong, D. J., McMahon, S. B., et al. (2009). Expression of the regeneration-associated protein SPRR1A in primary sensory neurons and spinal cord of the adult mouse following peripheral and central injury. *J. Comp. Neurol.* 513, 51–68. doi: 10.1002/cne.21944

Taylor, B. K., Fu, W., Kuphal, K. E., Stiller, C. O., Winter, M. K., Chen, W., et al. (2014). Inflammation enhances Y1 receptor signaling, neuropeptide Y-mediated inhibition of hyperalgesia, and substance P release from primary afferent neurons. *Neuroscience* 256, 178–194. doi: 10.1016/j.neuroscience.2013.10.054

Tsujii, M., Akeda, K., Iino, T., and Uchida, A. (2009). Are BMPs involved in normal nerve and following transection?: a pilot study. *Clin. Orthop. Relat. Res.* 467, 3183–3189. doi: 10.1007/s11999-009-1009-1

Tsujino, H., Kondo, E., Fukuoka, T., Dai, Y., Tokunaga, A., Miki, K., et al. (2000). Activating transcription factor 3 (ATF3) induction by axotomy in sensory and motoneurons: a novel neuronal marker of nerve injury. *Mol. Cell. Neurosci.* 15, 170–182. doi: 10.1006/mcne.1999.0814

Turner, M. R., Bowser, R., Bruijn, L., Dupuis, L., Ludolph, A., McGrath, M., et al. (2013). Mechanisms, models and biomarkers in amyotrophic lateral sclerosis. *Amyotroph. Lateral Scler. Frontotemporal Degener.* 14, 19–32. doi: 10.3109/21678421.2013.778554

Wang, J., Duncan, D., Shi, Z., and Zhang, B. (2013). WEB-based GEne SeT AnaLysis Toolkit (WebGestalt): update 2013. *Nucleic Acids Res.* 41, W77–W83. doi: 10.1093/nar/gkt439

Wills, A. M., Hubbard, J., Macklin, E. A., Glass, J., Tandan, R., Simpson, E. P., et al. (2014). Hypercaloric enteral nutrition in patients with amyotrophic lateral sclerosis: a randomised, double-blind, placebo-controlled phase 2 trial. *Lancet. 383*, 2065–2072. doi: 10.1016/S0140-6736(14)60222-1

Xiao, L., Michalski, N., Kronander, E., Gjoni, E., Genoud, C., Knott, G., et al. (2013). BMP signaling specifies the development of a large and fast CNS synapse. *Nat. Neurosci.* 16, 856–864. doi: 10.1038/nn.3414

Yang, B. B., Savin, M. A., and Green, M. (2012). Prevention of chemotherapy-induced neutropenia with pegfilgrastim: pharmacokinetics and patient outcomes. *Chemotherapy* 58, 387–398. doi: 10.1159/000345626

Yang, R., Gaidamakov, S. A., Xie, J., Lee, J., Martino, L., Kozlov, G., et al. (2011). La-related protein 4 binds poly(A), interacts with the poly(A)-binding protein MLLE domain via a variant PAM2w motif, and can promote mRNA stability. *Mol. Cell. Biol.* 31, 542–556. doi: 10.1128/MCB.01162-10

Zhang, B., Kirov, S., and Snoddy, J. (2005). WebGestalt: an integrated system for exploring gene sets in various biological contexts. *Nucleic Acids Res.* 33, W741–W748. doi: 10.1093/nar/gki475

Zhang, Y., Wang, L., Fu, Y., Song, H., Zhao, H., Deng, M., et al. (2009). Preliminary investigation of effect of granulocyte colony stimulating factor on amyotrophic lateral sclerosis. *Amyotroph. Lateral Scler.* 10, 430–431. doi: 10.3109/17482960802588059

Zhong, Z., Ohnmacht, J., Reimer, M. M., Bach, I., Becker, T., and Becker, C. G. (2012). Chondrolectin mediates growth cone interactions of motor axons with an intermediate target. *J. Neurosci.* 32, 4426–4439. doi: 10.1523/JNEUROSCI.5179-11.2012

**Conflict of Interest Statement:** Some authors are employees of Sygnis Bioscience. Armin Schneider is inventor on patent applications claiming the use of G-CSF for the treatment of neurodegenerative conditions. The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Received: 21 October 2014; accepted: 20 December 2014; published online: 20 January 2015.

Citation: Henriques A, Kastner S, Chatzikonstantinou E, Pitzer C, Plaas C, Kirsch F, Wafzig O, Krüger C, Spoelgen R, Gonzalez De Aguilar J-L, Gretz N and Schneider A (2015) Gene expression changes in spinal motoneurons of the SOD1G93A transgenic model for ALS after treatment with G-CSF. *Front. Cell. Neurosci.* 8:464. doi: 10.3389/fncel.2014.00464

This article was submitted to the journal *Frontiers in Cellular Neuroscience*.

Copyright © 2015 Henriques, Kastner, Chatzikonstantinou, Pitzer, Plaas, Kirsch, Wafzig, Krüger, Spoelgen, Gonzalez De Aguilar, Gretz and Schneider. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.