2010 Bilsland R2 Supporting Information · 10. Motulsky H, Christopoulos A (2004) Fitting models to...

19
1 Supporting Information Supplementary Materials and Methods Materials. Reagents were from Sigma-Aldrich unless otherwise specified. AlexaFluor555- maleimide, Texas Red-dextran (3,000 Da) and AlexaFluor488 α-bungarotoxin were from Invitrogen. The ChAT antibody was from Millipore and anti-myelin basic protein (MBP) from Roche. Animals. Transgenic mice carrying a human wild-type (WT) SOD1 (B6SJL- Tg[SOD1]2Gur/J) or mutant SOD1 G93A (TgN[SOD1-G93A]1Gur) gene were obtained from Jackson Laboratories. Colonies were maintained by breeding male heterozygous carriers with female (C57BL/6 x SJL) F1 hybrids. Mice were identified by genotyping for the human SOD1 transgene using DNA extracted from tail snips (1). Female thy1-mitoCFP-S (MitoMouse) (2) were bred with male SOD1 G93A to generate double heterozygous mice. Animals were housed in a controlled temperature and humidity environment and maintained on a 12 h light/dark cycle with access to food and water provided ad libitum. All experiments were carried out under license from the UK Home Office in accordance with the Animals (Scientific Procedures) Act 1986 and following approval from the Cancer Research UK and Institute of Neurology Ethical Review Committees. Neuronal cultures. Primary MN cultures were prepared from E12.5 WT mouse embryos as previously described (3). Cultures from 5 d in vitro were incubated with 40 nM H C 555 for 30 min at 37 o C, washed twice in Dulbecco’s modified Eagle’s medium, without phenol red, supplemented with 30 mM HEPES-NaOH, pH 7.4, and imaged using time-lapse confocal microscopy (4, 5).

Transcript of 2010 Bilsland R2 Supporting Information · 10. Motulsky H, Christopoulos A (2004) Fitting models to...

Page 1: 2010 Bilsland R2 Supporting Information · 10. Motulsky H, Christopoulos A (2004) Fitting models to biological data using linear and nonlinear regression. A practical guide to curve

1

Supporting Information

Supplementary Materials and Methods

Materials. Reagents were from Sigma-Aldrich unless otherwise specified. AlexaFluor555-

maleimide, Texas Red-dextran (3,000 Da) and AlexaFluor488 α-bungarotoxin were from

Invitrogen. The ChAT antibody was from Millipore and anti-myelin basic protein (MBP)

from Roche.

Animals. Transgenic mice carrying a human wild-type (WT) SOD1 (B6SJL-

Tg[SOD1]2Gur/J) or mutant SOD1G93A (TgN[SOD1-G93A]1Gur) gene were obtained from

Jackson Laboratories. Colonies were maintained by breeding male heterozygous carriers with

female (C57BL/6 x SJL) F1 hybrids. Mice were identified by genotyping for the human

SOD1 transgene using DNA extracted from tail snips (1). Female thy1-mitoCFP-S

(MitoMouse) (2) were bred with male SOD1G93A to generate double heterozygous mice.

Animals were housed in a controlled temperature and humidity environment and maintained

on a 12 h light/dark cycle with access to food and water provided ad libitum. All experiments

were carried out under license from the UK Home Office in accordance with the Animals

(Scientific Procedures) Act 1986 and following approval from the Cancer Research UK and

Institute of Neurology Ethical Review Committees.

Neuronal cultures. Primary MN cultures were prepared from E12.5 WT mouse embryos as

previously described (3). Cultures from 5 d in vitro were incubated with 40 nM HC555 for 30

min at 37oC, washed twice in Dulbecco’s modified Eagle’s medium, without phenol red,

supplemented with 30 mM HEPES-NaOH, pH 7.4, and imaged using time-lapse confocal

microscopy (4, 5).

Page 2: 2010 Bilsland R2 Supporting Information · 10. Motulsky H, Christopoulos A (2004) Fitting models to biological data using linear and nonlinear regression. A practical guide to curve

2

DRG neurons were isolated from adult SOD1G93A mice at 4 different disease stages: pre-

symptomatic, early symptomatic, symptomatic and late symptomatic and plated onto poly-L-

lysine and laminin-coated MatTek dishes (MatTek Inc) and cultured as previously described

(6). After 2 d in vitro, DRG neurons were incubated with 40 nM HC555 and 50 nM

Mitotracker Green for 30 min at 37˚C. Cells were washed twice in Dulbecco’s modified

Eagle’s medium (without phenol red, supplemented with 30 mM HEPES-NaOH, pH 7.4) and

imaged using time-lapse confocal microscopy (4, 5).

Histochemistry and immunofluorescence. To assess the distribution of HC555 at the NMJ

and its uptake into sciatic nerve axons, this probe was injected into the TA and GC muscles.

In selected cases, a sciatic nerve ligation at the mid-thigh level was performed before the

intramuscular injection. 6 h later (or 3, 5 and 8 h in the absence of nerve ligation), the animals

were transcardially perfused and the muscles and sciatic nerves removed, post-fixed and

cryoprotected prior to serial sectioning (muscles, 30 µm; sciatic nerve, 10 µm). Specimens

were allowed to dry and incubated in AlexaFluor488-α bungarotoxin (1:200 in distilled

water) for 1 h. Alternatively, confocal images of transverse sections of sciatic nerve were

taken and the number of axons containing HC555 that counterstained with the ChAT antibody

(1:250) were calculated as a percentage of the total number of axons. Myelin sheaths were

immunostained for MBP (1:500).

Silver cholinesterase staining. This technique, which combines the use of potassium

ferricyanide for the staining of cholinesterase at NMJs and silver impregnation (modified

from 7), allowed the detection of motor axons and nerve terminals and was used to monitor

the percentage of denervation. Muscles were removed and immediately fixed in a slightly

stretched position using buffered-formaldehyde (4% formaldehyde, 1% CaCl2, 0.5% MgCl2

6H2O, 0.1% CdCl2 2.5H2O) in veronal-acetate buffer (pH 6.45) for 6 h. Subsequently,

Page 3: 2010 Bilsland R2 Supporting Information · 10. Motulsky H, Christopoulos A (2004) Fitting models to biological data using linear and nonlinear regression. A practical guide to curve

3

muscles were incubated in 10% sucrose overnight at 4˚C. Longitudinal sections (40 µm

thick) were then cut on a cryomicrotome and held in distilled water on ice for 10 min,

followed by incubation in an acetylthiocholine iodide solution (1.6 mM acetylthiocholine

iodide in 0.06 M maleate buffer (pH 6.0), containing 5 mM sodium citrate, 3 mM copper

sulphate, 0.5 mM potassium ferricyanide and 0.4 M sucrose) for 20 min on ice. Sections were

then rinsed in distilled water for 30 s and immersed in potassium ferricyanide (7.3 mM in

distilled water) for 10 min at room temperature to allow an iodide precipitate to develop.

Sections were washed three times in distilled water for 5 min, immersed in absolute ethanol

for 1 h then washed in distilled water twice for 15 min. Sections were then placed in a silver

solution (0.1% CaCO3, 0.05% CuSO4 5H2O, 10% AgNO3 in distilled water) for 35 min at

37˚C to visualise the axons. After a rinse in distilled water, sections were immersed in a

reducer solution (1% hydroquinone, 10% Na2SO3 in distilled water), rinsed again and

mounted onto slides coated with gelatin. Over 500 end plates were examined and quantified

for each muscle.

MN survival. MN survival was monitored on serial sections (20 µm) of lumbar spinal cord

(L2-L6) were cut on a cryostat and stained with cresyl fast violet, a Nissl stain (8). The

number of surviving MNs within the sciatic motor pool of each ventral horn was assessed by

counting the number of Nissl-stained MNs in every third spinal cord section (n=60). This

method avoided counting the same cell twice in consecutive sections. Morphological criteria

were also employed, so that only large polygonal neurons in which the nucleolus was clearly

visible at high magnification were included in the counts. This method of counting has

previously been used to assess MN survival (9). To assess the extent of MN degeneration,

mean MN survival in the sciatic motor pool of SOD1G93A mice was compared to the number

surviving in WT littermates and SOD1WT mice.

Page 4: 2010 Bilsland R2 Supporting Information · 10. Motulsky H, Christopoulos A (2004) Fitting models to biological data using linear and nonlinear regression. A practical guide to curve

4

Uptake by lumbar MNs and DRG neurons. Following each in vivo imaging experiment,

mice were terminally anaesthetized with sodium pentobarbitone and perfused transcardially

with 4% paraformaldehyde (PFA; TAAB Ltd) in 77 mM Na2HPO4, 23 mM NaH2PO4. The

TA and GC muscles, the imaged sciatic nerve and the L2-L6 section of the spinal cord were

removed, post-fixed for 4 h and then cryoprotected in 30% sucrose in PBS at 4˚C. Serial

sections (20 µm) of the L2-L6 region of the spinal column were cut and unstained sections

were used to examine HC555 localization in lumbar MNs and DRGs using confocal

microscopy. Sections were excited using the 543 nm laser line and emitted fluorescence

collected using the 560-615 nm bandpass filter.

TrkB antibody studies. The sciatic nerve was ligated at the mid-thigh level of both WT and

SOD1G93A mice at 75 d of age. Subsequently 8 µl of Alexafluor555-conjugated TrkB

antibody (Millipore) together with BDNF was injected into the TA and GC muscles.

Approximately 6 h later the sciatic nerve was gently squeezed and the extruded axoplasm

snap frozen. TrkB antibody was detected by western blotting using a horseradish peroxidase-

conjugated swine anti-rabbit antibody (Dako) and quantified by ImageJ. Data were

normalised for GAPDH levels using a horseradish peroxidase-conjugated anti-GAPDH rabbit

polyclonal antibody (Cell Signaling).

Statistical analysis. The parameters in Fig. 1-5, S2-S5 and S7 {AREA, SD, mean speed},

their standard errors and confidence intervals were estimated using non-linear regression (10,

11) for each sample separately, regressing the relative frequencies against the lower bound of

the speed bin, fitting against the model: Freq=(AREA/(SD*(2*pi)^0.5))*exp(-0.5*((speed-

mean speed)/SD)^2). The parameter 'mean speed', along with its standard error, were used in

two-tailed z-tests to infer the equivalence or otherwise the location of the peaks amongst the

pairwise comparisons. The p-values were then corrected using the conservative Bonferroni

Page 5: 2010 Bilsland R2 Supporting Information · 10. Motulsky H, Christopoulos A (2004) Fitting models to biological data using linear and nonlinear regression. A practical guide to curve

5

correction method. An omnibus F-test of the 'mean speed' across all the groups was passed at

an alpha of 0.05 (p<0.0001 F=173.6 with 6:329 degrees of freedom – the p-value could not

be determined to more precision than this inequality by our software). R2 statistics of the non-

linear regression were all found to be >0.9, and residuals were visually checked, and

determined to be free of any systematic effects that would invalidate our approach.

The estimated values of the mean speed parameter, for samples labelled WT or early

symptomatic SOD1G93A, were then used to build a linear discriminant classifier (12) to find

the boundary between the two classes, and then predict the classes of the unknown samples

(pre-symptomatic data, both 25 d and 36 d SOD1G93A mice) by seeing on which side of the

mean speed boundary they lie. The difference between WT and early symptomatic SOD1G93A

mean speeds were tested using two slightly different approaches. In the 'Separate' strategy,

each of the peaks was summarised (as before) using non-linear regression to extract the

modal speeds, and the n=5 (WT), n=6 (early symptomatic SOD1G93A) replicates in each

group were subjected to a α=0.05 two-tailed Student's t-test, which passed with p=1.23e-15.

Graphs didn't provide any visual evidence that the requirement for normality of group means

was violated. In the 'Together' approach, we estimated two curves through the dataset, one

through the WT replicates (n=5), and one through the early symptomatic SOD1G93A replicates

(n=6) with non-linear regression, and an F-test at α=0.5 indicated a significant difference in

the modal speeds of the two groups (p=6.53e-75, F=471.6 on 1:494 degrees of freedom).

Thus both approaches found a significant difference between WT and early symptomatic

SOD1G93A mean speeds at α=0.05. Calculations were performed using GraphPad Prism

version 4.03 for Windows (GraphPad Software).

Statistical analysis of MN counts, TrkB transport, muscle denervation and mitochondrial

pausing were performed using an unpaired two-tailed Student’s t-test with equal variance

Page 6: 2010 Bilsland R2 Supporting Information · 10. Motulsky H, Christopoulos A (2004) Fitting models to biological data using linear and nonlinear regression. A practical guide to curve

6

using Kaleidagraph or SigmaStat (2.03; SyStat Software) and the normality of the data was

assessed by Sigma Stat. Significance for all statistical analysis was set at p<0.05.

References

1. Gurney ME, et al. (1994) Motor neuron degeneration in mice that express a human

Cu,Zn superoxide dismutase mutation. Science 264:1772-1775.

2. Misgeld T, Kerschensteiner M, Bareyre FM, Burgess RW, Lichtman JW (2007)

Imaging axonal transport of mitochondria in vivo. Nat Methods 4:559-561.

3. Bilsland LG, Nirmalananthan N, Yip J, Greensmith L, Duchen MR (2008) Expression

of mutant SOD1 in astrocytes induces functional deficits in motoneuron

mitochondria. J Neurochem 107:1271-1283.

4. Lalli G, Gschmeissner S, Schiavo G (2003) Myosin Va and microtubule-based motors

are required for fast axonal retrograde transport of tetanus toxin in motor neurons. J

Cell Sci 116:4639-4650.

5. Kieran D, et al. (2005) A mutation in dynein rescues axonal transport defects and

extends the life span of ALS mice. J Cell Biol 169:561-567.

6. Deinhardt K, et al. (2006) Rab5 and Rab7 control endocytic sorting along the axonal

retrograde transport pathway. Neuron 52:293-305.

7. Namba T, Nakamura T, Grob D (1967) Staining for nerve fiber and cholinesterase

activity in fresh frozen sections. Am J Clin Pathol 47:74-77.

8. Culling CFA (1963) Handbook of Histopathological Techniques. 2nd Ed. pp 362-363.

9. Bilsland LG, et al. (2006) Increasing cannabinoid levels by pharmacological and

genetic manipulation delay disease progression in SOD1 mice. Faseb J 20:1003-

1005.

Page 7: 2010 Bilsland R2 Supporting Information · 10. Motulsky H, Christopoulos A (2004) Fitting models to biological data using linear and nonlinear regression. A practical guide to curve

7

10. Motulsky H, Christopoulos A (2004) Fitting models to biological data using linear

and nonlinear regression. A practical guide to curve fitting Oxford University Press,

New York.

11. Bates DM, Watts DG (1980) Non linear regression analysis and its applications.

Wiley, New York.

12. Fisher RA (1936) The use of multiple measurements in the taxonomic problems. Ann

Eugenics 7:179-188.

13. Chiu AY, et al. (1995) Age-dependent penetrance of disease in a transgenic mouse

model of familial amyotrophic lateral sclerosis. Mol Cell Neurosci 6:349-362.

14. Mourelatos Z, Gonatas NK, Stieber A, Gurney ME, Dal Canto MC (1996) The Golgi

apparatus of spinal cord motor neurons in transgenic mice expressing mutant Cu,Zn

superoxide dismutase becomes fragmented in early, preclinical stages of the disease.

Proc Natl Acad Sci U S A 93:5472-5477.

15. Pun S, Santos AF, Saxena S, Xu L, Caroni P (2006) Selective vulnerability and

pruning of phasic motoneuron axons in motoneuron disease alleviated by CNTF. Nat

Neurosci 9:408-419.

16. Hegedus J, Putman CT, Gordon T (2007) Time course of preferential motor unit loss

in the SOD1(G93A) mouse model of amyotrophic lateral sclerosis. Neurobiol Dis

28:154-164.

17. Tu PH, et al. (1997) Selective degeneration fo Purkinje cells with Lewy body-like

inclusions in aged NFHLACZ transgenic mice. J Neurosci 17:1064-1074.

18. Johnston JA, Dalton MJ, Gurney ME, Kopito RR (2000) Formation of high molecular

weight complexes of mutant Cu, Zn-superoxide dismutase in a mouse model for

familial amyotrophic lateral sclerosis. Proc Natl Acad Sci U S A 97:12571-12576.

Page 8: 2010 Bilsland R2 Supporting Information · 10. Motulsky H, Christopoulos A (2004) Fitting models to biological data using linear and nonlinear regression. A practical guide to curve

8

19. Fischer LR, et al. (2004) Amyotrophic lateral sclerosis is a distal axonopathy:

evidence in mice and man. Exp Neurol 185:232-240.

20. Frey D, et al. (2000) Early and selective loss of neuromuscular synapse subtypes with

low sprouting competence in motoneuron diseases. J Neurosci 20:2534-2542.

Supplementary Figure Legends

Figure S1: Internalization of HC555 by motor axons. Following intramuscular injection,

HC555 accumulated at the NMJ within 3 h (A; in red), and then declined (B, 5 h; C, 8 h; scale

bar, 5 µm) due to uptake by axons (arrowhead) in the sciatic nerve (D; dashed box has been

magnified; scale bar, 10 µm). The motor end plate has been counterstained with

AlexaFluor488 α-bungarotoxin (A-D; in green).

Figure S2: The rate of retrograde transport of HC555 was not influenced by age. HC555 was

injected intramuscular into female WT mice of varying ages. Following a recovery period,

the mice were anaesthetized, the sciatic nerve exposed and retrograde transport imaged.

Figure S3: Kinetic analysis of HC555 retrograde transport. To control for the effects of

SOD1G93A overexpression, WT littermates and SOD1WT mice were used as controls (A-C).

No significant difference in MN survival was found between WT and SOD1WT mice (A;

p>0.05; n=3). Quantitative analysis of retrograde transport (B) however, indicated that

transport of HC555 was significantly faster in SOD1WT mice (200 carriers; 30 axons; n=3)

compared to WT mice (258 carriers; 44 axons; n=5). This increase in speed in SOD1WT mice

was due to a higher percentage of fast carriers (C). Displacement of HC555 positive carriers

tracked during the course of a representative experiment from WT (D) and SOD1G93A mice at

pre-symptomatic (E) and early symptomatic (F) stages of disease. The start of tracking for

Page 9: 2010 Bilsland R2 Supporting Information · 10. Motulsky H, Christopoulos A (2004) Fitting models to biological data using linear and nonlinear regression. A practical guide to curve

9

each carrier, which is displayed by a single dotted line, was set to time = 0. The slope of the

line is proportional to the speed of transport. (G) HC555 accumulated at the NMJ (labeled

with AlexaFluor488 α-bungarotoxin in merged image; asterisks) and was taken up into

innervating axons (arrowhead; see also Fig. S1) of SOD1G93A mice at the late symptomatic

stage of disease. In the sciatic nerve, HC555 was found in axons co-labeled for ChAT,

indicating that they were motor axons (H; arrowheads). In addition, a punctate HC555 signal

was detected in MN cell bodies located in the spinal cord of late symptomatic mice (I). Scale

bars, 20 µm. Error bars represent S.E.M..

Figure S4: Retrograde transport was unaffected in DRG neurons from SOD1G93A mice.

HC555 was internalized by sensory neurons projecting in the sciatic nerve (A). Following

intramuscular injection, sections (20 µm) of the spinal column revealed accumulations of

HC555 in the DRG somas (arrowheads; scale bar, 50 µm). (B) DRGs were removed from

adult SOD1G93A mice at varying stages of disease and maintained in culture. Kinetic analysis

of HC555 transport revealed that there was no significant difference in the rates of transport in

SOD1G93A DRG neurons (56–235 carriers; 10-29 axons; n ≥ 2 independent cultures) at any

disease stage examined, compared to WT.

Figure S5: The onset of retrograde transport deficits in pre-symptomatic SOD1G93A mice. (A)

Analysis of the rate of retrograde transport in pre-symptomatic SOD1G93A mice at 36 d

following intramuscular injection, identified two distinct disease-stage sub-populations; one

(grey circles; 140 carriers; 30 axons) characterized by a peak speed of retrograde HC555

carriers statistically indistinguishable from that of WT animals (black; 258 carriers; 44 axons)

and another pool (cyan; 204 carriers; 22 axons) in which the peak speed did not differ

significantly from early symptomatic SOD1G93A mice (red; 193 carriers; 39 axons). This

suggests that this pre-symptomatic stage represents the approximate age at which retrograde

Page 10: 2010 Bilsland R2 Supporting Information · 10. Motulsky H, Christopoulos A (2004) Fitting models to biological data using linear and nonlinear regression. A practical guide to curve

10

transport deficits became evident in SOD1G93A mice. Indeed, analysis of retrograde transport

at an earlier stage of disease (25 d; purple), revealed peak speeds that did not differ

significantly from WT. Error bars represent S.E. of the mean speed. (B) Speed distribution

curves of HC555 transport for WT and SOD1G93A mice at 25 d (164 carriers; 24 axons). Error

bars represent S.E.M..

Figure S6: Muscle denervation in SOD1G93A mice at 36 d of age. Fixed muscle sections from

36 d old WT and SOD1G93A mice underwent silver cholinesterase staining to label the motor

end plates and the innervating axons. Muscle sections from WT mice; scale bar, 100 µm (A)

and 25 µm (B). The number of denervated motor end plates was quantified and displayed in

(C). There was no significant difference in the innervation of NMJ between WT and

SOD1G93A mice at a pre-symptomatic stage. 7.0±0.9% (n=6 muscles) and 7.5±1.1% (n=6

muscles) of NMJs were denervated in WT EDL and soleus muscles, respectively, whereas in

SOD1G93A mice the relative values were 9.0±0.7% (n=6 muscles; EDL) and 9.4±1.4 (n=6

muscles; soleus). Error bars represent S.E.M.; n ≥ 500 end plates from at least 4 mice.

Figure S7: Axonal transport of mitochondria was unaffected in DRG neurons isolated from

SOD1G93A mice. DRG neurons were removed from both WT and SOD1G93A mice at 36 d and

maintained in culture. Cultures were loaded with Mitotracker Green (50 nM) for 30 min at

37˚C. Axonal transport of Mitotracker was monitored by time-lapse confocal microscopy.

Quantitative kinetic analysis revealed that there was no significant difference in the rates of

transport of mitochondria in either the anterograde (a) or retrograde (b) direction in DRG

neurons from SOD1G93A mice (in green) compared to WT (in black). Error bars represent

S.E.M..

Page 11: 2010 Bilsland R2 Supporting Information · 10. Motulsky H, Christopoulos A (2004) Fitting models to biological data using linear and nonlinear regression. A practical guide to curve

11

Table legend

Table S1: Summary of the key pathological signs in SOD1G93A mice during disease

progression. This table attempts to correlate our observed changes in axonal retrograde

transport and MN survival (in bold) with the documented selective vulnerabilities of different

MN pools and with overt signs of disease seen in SOD1G93A mice. FF, fast, fatiguable; FR,

fast, resistant.

Movies

Movie 1: Retrograde transport of HC555 in a single axon in the intact sciatic nerve of WT

mice. HC555 was injected intramuscularly into WT mice. After a recovery period, the mouse

was terminally anaesthetized and the sciatic nerve exposed. Retrograde transport of HC555

was observed in the direction of the spinal cord (to the left) within individual axons of the

intact sciatic nerve. The video consists of 100 frames and is played at 5 frames/s. The image

size is 71.6 µm x 12.7 µm.

Movie 2: Retrograde transport of HC555 in DRG axons. This video shows a DRG axon that

has been incubated with 40 nM HC555 for 30 min at 37˚C and imaged by time-lapse

microscopy. This video consists of 100 frames and is played at 5 frames/s. The image size is

70.0 µm x 12.1 µm.

Page 12: 2010 Bilsland R2 Supporting Information · 10. Motulsky H, Christopoulos A (2004) Fitting models to biological data using linear and nonlinear regression. A practical guide to curve

Table S1. Key pathological signs in SOD1G93A mice during disease progression

Disease stage

Disease signs

MN pathology

Axonal transport

Pre-symptomatic (36±0.95 d)

No overt disease signs.

Ultrastructural pathology including abnormal mitochondria, Golgi fragmentation and vacuolation (13, 14).

No MN loss.

Synaptic vesicle stalling (35-38 d).

Loss of synaptic vesicle markers in FF MNs (38-46 d)(15).

Deficit in retrograde transport of HC555 (this work).

Deficit in slow anterograde and fast axonal transport (5, 15).

Early symptomatic

(74±1.4 d) Loss of tetanic muscle force in fast muscles (16).

Loss of approximately 23% of lumbar MNs.

FF innervated muscle fibres have become denervated. Onset of denervation of FR innervated muscle fibres (15)

Protein aggregation and neurofilament inclusions in MNs (17, 18)

Greater deficit in retrograde transport of HC555 as a consequence of vulnerable fast MN degeneration (this work).

Symptomatic (95±1.4 d)

Fine hindlimb tremor, muscle spasticity and hyperreflexia (13).

Further loss of muscle force/grip strength.

Decline in constant rotarod performance (19)

Loss of approximately 39% of lumbar MNs.

Innervation by fast MNs predominantly lost. Compensatory sprouting of slow MNs to reinnervate previously fast NMJs (15).

Deficit in retrograde transport of HC555, but less severe than at early-symptomatic stage due to death of fast MNs and apparent compensation by slow, relatively unaffected MNs, and DRGs (this work).

Late symptomatic

(113±3.5 d)

Weight loss.

Hindlimb paralysis and muscle atrophy (13).

Death at approximately 130 d (1).

Loss of approximately 56% of lumbar MNs.

Denervation of soleus muscle innervated by slow MNs (20).

Deficit in retrograde transport of HC555 as at the symptomatic disease stage (95 d; this work).

Page 13: 2010 Bilsland R2 Supporting Information · 10. Motulsky H, Christopoulos A (2004) Fitting models to biological data using linear and nonlinear regression. A practical guide to curve

α-bungarotoxin HC555 MergeA

B

C

HC555

Merge

3 h

5 h

8 h

Figure S1 - Bilsland et al.

Dα-bungarotoxin

Page 14: 2010 Bilsland R2 Supporting Information · 10. Motulsky H, Christopoulos A (2004) Fitting models to biological data using linear and nonlinear regression. A practical guide to curve

Rel

ativ

e fre

quen

cy

Speed (μm/s)-0.4 0 0.4 0.8 1.2 1.6 2.0 2.4 2.8 3.2

38 d WT81 d WT101 d WT125 d WT

0

0.02

0.04

0.06

0.08

0.1

0.12

0.14

0.16

Figure S2 Bilsland et al.

Page 15: 2010 Bilsland R2 Supporting Information · 10. Motulsky H, Christopoulos A (2004) Fitting models to biological data using linear and nonlinear regression. A practical guide to curve
Page 16: 2010 Bilsland R2 Supporting Information · 10. Motulsky H, Christopoulos A (2004) Fitting models to biological data using linear and nonlinear regression. A practical guide to curve

Speed (μm/s)

Rel

ativ

e fre

quen

cy

-0.4 0 0.4 0.8 1.2 1.6 2.0 2.40

0.05

0.10

0.15

0.20

0.25

0.30

WT

Pre-symptomatic SOD1G93A

Early symptomatic SOD1G93A

Symptomatic SOD1G93A

Late symptomatic SOD1G93A

PhaseHC555

DRGs

A

B

Figure S4 - Bilsland et al.

Page 17: 2010 Bilsland R2 Supporting Information · 10. Motulsky H, Christopoulos A (2004) Fitting models to biological data using linear and nonlinear regression. A practical guide to curve

Figure S5 - Bilsland et al.

WT

(25 d)

Pre-sym

ptomati

c

Early s

ympto

matic

(36 d)

(74 d)

Pre-sym

ptomati

c

1.2

0.8

1.0

1.4

1.6

Spee

d (µ

m/s

)

0

0.02

0.04

0.06

0.08

0.10

0.12

0.14

WTSOD1G93A

-0.4 0 0.4 0.8 1.2 1.6 2.0 2.4 2.8 3.2

Rel

ativ

e fre

quen

cy

Speed (μm/s)

0.16

A

B

Page 18: 2010 Bilsland R2 Supporting Information · 10. Motulsky H, Christopoulos A (2004) Fitting models to biological data using linear and nonlinear regression. A practical guide to curve

EDL SOL

WT

2

4

6

8

10

12

0

% o

f den

erva

tion

CA

B

Figure S6 - Bilsland et al

SOD1G93A

Page 19: 2010 Bilsland R2 Supporting Information · 10. Motulsky H, Christopoulos A (2004) Fitting models to biological data using linear and nonlinear regression. A practical guide to curve

Figure S7 - Bilsland et al

A

B

0 0.4 0.8 1.2 1.6 2.0 2.4Speed (μm/s)

0

0.1

0.2

0.3

0.4

Rel

ativ

e fre

quen

cy

WTSOD1G93A

- 0.4

Rel

ativ

e fre

quen

cy

0

0.1

0.2

0.3

0.4

WTSOD1G93A

0 0.4 0.8 1.2 1.6 2.0 2.4- 0.4

Retrograde transport

Anterograde transport

Mitochondrial speed

Speed (μm/s)