Loss of BICD2 in muscle drives motor neuron loss in a developmental form of spinal muscular atrophy.

Autosomal dominant missense mutations in BICD2 cause Spinal Muscular Atrophy Lower Extremity Predominant 2 (SMALED2), a developmental disease of motor neurons. BICD2 is a key component of the cytoplasmic dynein/dynactin motor complex, which in axons drives the microtubule-dependent retrograde transport of intracellular cargo towards the cell soma. Patients with pathological mutations in BICD2 develop malformations of cortical and cerebellar development similar to Bicd2 knockout (-/-) mice. In this study we sought to re-examine the motor neuron phenotype of conditional Bicd2-/- mice. Bicd2-/- mice show a significant reduction in the number of large calibre motor neurons of the L4 ventral root compared to wild type mice. Muscle-specific knockout of Bicd2 results in a similar reduction in L4 ventral axons comparable to global Bicd2-/- mice. Rab6, a small GTPase required for the sorting of exocytic vesicles from the Trans Golgi Network to the plasma membrane is a major binding partner of BICD2. We therefore examined the secretory pathway in SMALED2 patient fibroblasts and demonstrated that BICD2 is required for physiological flow of constitutive secretory cargoes from the Trans Golgi Network to the plasma membrane using a VSV-G reporter assay. Together, these data indicate that BICD2 loss from muscles is a major driver of non-cell autonomous pathology in the motor nervous system, which has important implications for future therapeutic approaches in SMALED2.

Introduction

Spinal Muscular Atrophy Lower Extremity Predominant (SMALED) is a disease of lower motor neurons, principally affecting the lower limbs. Affected individuals often present at birth with contractures of the lower limbs. Autosomal dominant missense mutations in DYNC1H1 and BICD2 are the only two known genetic causes of SMALED and determine two indistinguishable forms of this disease [1-5].

DYNC1H1 encodes the cytoplasmic dynein heavy chain, a key component and force-generating subunit of the dynein/dynactin retrograde transport complex. This complex is responsible for the transport of intracellular cargoes towards the minus-end of microtubules located in the cell soma [6]. In contrast, anterograde transport towards the positive end of microtubules is driven by the kinesin family of motor proteins [7].

Motor neurons have extremely long axons, which make them preferentially susceptible to deficits in axonal transport [8]. It has long been thought that mutations in BICD2 and DYNC1H1 impair axonal transport leading to motor neuron degeneration. This is supported by work in Drosophila demonstrating a reduction of the in vitro run length in flies expressing disease-causing mutant DYNC1H1 compared to wild type [9]. This conclusion has been brought into question by the contrasting findings that several disease-causing mutations in BICD2 increase its binding affinity to the dynein/dynactin complex resulting in an increase in run length [10].

BICD2 is a cargo adaptor protein, comprising a N-terminal region that mediates binding to the dimerization domain of DYNC1H1 and a C-terminal cargo binding domain [11]. Missense mutations throughout BICD2 have been shown to cause SMALED2 [12], yet the disease mechanism is unclear. Mutations in the N-terminal domain increase the affinity to the dynein complex [10], whereas the p.Glu774Gly mutation in the C-terminal domain has no effect on DYNC1H1 binding, but disrupts the interaction with its cargo Rab6, causing a loss of function phenotype [5]. As BICD2 forms dimers, such mutations may also impair the cargo binding ability of wild type/mutant BICD2 heterodimers, leading to near complete loss of function as opposed to haploinsufficiency. A loss of function effect in SMALED2 was thought to be unlikely as the Bicd2 mouse was originally reported to lack a motor phenotype [13]. However, recent studies of SMALED2 patients with severe mutations in BICD2 have revealed additional cortical and cerebellar phenotypes similar to the Bicd2 mouse. These findings suggest that pathological BICD2 mutations may indeed induce a loss of function [14]. In this study, we sought to re-examine the motor neuron phenotype in the Bicd2 mouse to conclusively assess the molecular basis of the pathomechanism of SMALED2.

Materials and methods

Animals

Mice homozygous for an allele in which the loxP sequence is inserted into intron 1 and the 3’UTR of the endogenous Bicd2 gene [13] were crossed with heterozygous deleter-Cre mice (C57BL/6NTac-Gt(ROSA)26Sor) [15] to generate Bicd2 heterozygous knockout mice expressing a Cre-recombinase transgene. These mice were subsequently crossed with C57Bl/6 J mice to generate Bicd2 mice with no Cre recombinase transgene. Genotypes of mice were determined by PCR of ear clip DNA (P1 = AATGGAGAAGATCTCATCTTGGCAGG, P2 = GTGTAGCACTTCAGGAACATCCATGC, P3 = TGTCAGCAAACTCCATCTCTAGCCTC, P261 = CGGCGGCATCAGAGCAGCCGATTG).

To generate muscle-specific Bicd2 knockout mice, homozygous Bicd2 mice (C57Bl/6 J) were crossed with knockin MyoD-Cre mice (FVB.Cg-Myod1/J, Jackson laboratory) [16] and back-crossed to the original Bicd2 progeny for six generations to produce a congenic strain. To generate motor neuron-specific Bicd2 knockout mice, homozygous Bicd2 mice (C57Bl/6 J) were crossed with ChAT-IRES-Cre knock-in mice (B6;129S6-Chat/J, Jackson laboratory) [17] and then back-crossed to the original Bicd2 progeny to generate a homozygous Bicd2 background.

Histopathology

Mice were terminally anaesthetised with intraperitoneal injection of pentobarbitone followed by thoracotomy and transcardial perfusion with 10 ml 0.9% NaCl followed by 20 ml 4% paraformaldehyde (PFA, Fisher Scientific) in phosphate buffered saline (PBS). L4 DRG and L4 dorsal and ventral roots were dissected from PFA-fixed mice and post-fixed in PFA/glutaraldehyde buffer at 4 °C (20 ml 0.9% saline, 10 ml 10% PFA, 10 ml 10% glutaraldehyde (Sigma-Aldrich), 20% dextran (20,000 MW, Sigma-Aldrich) made up to 100 ml with 0.1 M PIPES-NaOH, pH 7.6, for 24 h before fixation in 1% osmium tetroxide (Agar Scientific) and processing into araldite CY212 epoxy resin (Agar Scientific) through graded alcohols and propylene oxide solutions using a standard protocol. Semi-thin sections (0.8 μm) were cut on an Ultracut E ultramicrotome (Leica), stained with 1% toluidine blue containing 1% borax (BDH), and examined with a Leica light microscope using an oil immersion 100x lens. Images were generated by image stitching using HUGIN-panorama photo stitching. Axon diameter was calculated using an in-house semi-automated thresholding programme using Definiens image analysis software.

Immunohistochemistry (For a summary of antibodies see Table 1)

Muscle fibre subtype and area quantification

P21 mice were euthanized by cervical dislocation and skeletal muscles (gastrocnemius) were dissected and snap frozen in liquid nitrogen-cooled isopentane and stored at − 80 °C for cryosectioning. Immunofluorescent staining was carried out on 10 μm frozen sections. Sections were blocked in PBS 0.2% Triton X-100, 5% goat serum for 1 h. Primary and secondary antibodies were diluted in PBS 0.2% Triton X-100, 2% goat serum. Primary antibodies: Type 1 BA-D5 (mouse IgG2b) DSHB (1 in 100), Type IIA SC-71 (mouse IgG1) DSHB (1 in 100), Type IIB BF-F3 (mouse IgM) DSHB (1 in 100) and laminin (rabbit) L9393 Sigma-Aldrich (1 in 100) were incubated for 1 h at 37 °C. Sections were then washed three times with PBS followed by goat anti-mouse IgG2b conjugated to AlexaFluor488, goat anti-mouse IgG1 conjugated to AlexaFluor647, goat anti-mouse IgM conjugated to AlexaFluor568 and goat anti-rabbit conjugated to AlexaFluor410 (all from Life Technology, 1:500). Secondary antibodies were incubated for 30 min at 37 °C and washed three times with PBS. Cover slips were mounted on stained sections using Dako fluorescent mounting media and dried for 48 h. Z-stack stitched images were acquired using a Zeiss LSM 710 confocal microscope using previously described parameters [18]. Quantification of the fibre number, area and type were performed using the muscle-J automated freeware and Image J as previously described [18].

NMJ analyses

Muscles were processed for immunohistochemistry and NMJ phenotypes scored as previously described [19]. 100 NMJs were scored per mouse for poly-innervation and occupancy counts, whilst 16–20 NMJs were used for area measurements.

Muscle spindle assessment

Soleus muscles were dissected and processed, and spindles analysed as described previously with minor modifications [20]. Briefly, soleus muscles were dissected then fixed overnight in 4% PFA in PBS prior to overnight equilibration in 20% (w/v) sucrose in PBS before freezing in Tissue-Tek O.C.T. (Sakura Finetek). 20 μm transverse serial sections throughout the entire muscle were cut onto three parallel slides for immunohistochemistry and analyses using an OTF Cryostat (Bright Instruments).

DRG dissection and staining

Mice were perfused with 4% PFA in PBS before L4 DRG were dissected as previously described [21]. DRG were sectioned and stained as previously published [20]. Briefly, DRG were post-fixed overnight in 4% PFA in PBS, before embedding in Tissue-Tek O.C.T., and sectioning at 10 μm across four parallel polysine-coated slides (VWR, 631–0107) with an OTF Cryostat. DRG sections were permeabilised for 30 min in PBS containing 0.3% Triton X-100 and blocked for 30 min in 10% bovine serum albumin (BSA) and 0.3% Triton X-100 in PBS, before probing overnight at 4 °C with primary antibodies (1:500 mouse anti-NF200 [N0142; Sigma-Aldrich] and 1:500 rabbit anti-peripherin (AB1530; Merck Millipore) in blocking solution. Sections were then washed with PBS for 30 min, probed for 2 h with secondary antibodies (1:1000) in PBS, before washing in PBS, flooding with fluorescence mounting medium (S3023; Dako) and covering with a 22 × 50 mm cover glass (VWR).

VSVG assay

Frozen vials of human skin fibroblasts were thawed and grown in T75 flasks for seven days in DMEM media supplemented with 20% foetal calf serum and L-glutamine (all from Gibco). Cells were trypsinised and plated onto 13 mm untreated glass coverslips (50,000 per coverslip). After 24 h, fibroblasts were transfected with a GFP-VSV-G ts045 plasmid using Lipofectamine 3000 (Invitrogen) as per manufacturer’s instructions and grown at 37 °C (VSV = vesicular stomatitis virus). After 8 h, cells were incubated at 40 °C for 14 h before being transferred to an incubator at 32 °C. Cells were fixed at 0, 60, 120, 180, 240 and 300 min after incubation at 32 °C in PFA 4% in PBS. Non-permeabilised cells were then incubated with PBS 1% BSA for 15 min prior to overnight incubation with mouse anti-VSV-G [8G5F11; Kerafast] antibody (1:2000) at 4 °C. Cells were washed three times with PBS before incubation with a donkey anti-mouse AlexaFluor568 secondary antibody (Life Technology, 1:500) for 1 h at room temperature. Cells were then washed, stained with DAPI (1:2000) and coverslips mounted and fixed with DAKO fluorescent mounting media. Ten fibroblasts per condition were imaged with a Zeiss LSM 510 confocal microscope and the ratio of the fluorescence between plasma membrane VSV-G (594 nm excitation) to the total VSV-G (GFP signal, 488 nm excitation) was quantified as previously described [22].

Brefeldin a assay

Human fibroblasts were cultured as previously described above, but were plated on 20 mm coverslips in 6 well plates at a density of 100,000 per well. After 24 h in culture, cells were transfected with a plasmid encoding human galactosyltransferase tagged with GFP (GT–GFP) (kind gift from Masayuki Murata, Tokyo). 24 h post-transfection, cells were imaged using a Zeiss LSM 710 confocal microscope. Briefly, cells were imaged at 37 °C, 5% CO2 with a 40x objective. The Golgi apparatus was identified by GFP staining. Brefeldin A ((BFA), 10 μg/ml final concentration; Sigma-Aldrich) was re-suspended in DMEM, 10% FCS, 20 mM HEPES-NaOH (final concentration), pH 7.3, and added to the coverslip at time 0. Time-lapse imaging at 5 s intervals was carried out and the time recorded to the beginning and end of the ‘Golgi blush’, defined as the stereotypical changes in Golgi morphology leading to its resorption into the ER.

Statistical analysis

Data were assumed to be normally distributed unless evidence to the contrary was provided by the D’Agostino and Pearson omnibus normality test. Data were statistically analysed using an unpaired t test, Mann-Whitney U-test, or one-way ANOVA with post-hoc Dunnett’s test against a control mean. GraphPad Prism 6 software was used for all statistical analyses. Means ± standard error of the mean (SEM) were plotted for all graphs.

Results and discussion

Bicd2 mice display a significant motor neuron loss

Due to the similarities between severe SMALED2 patients and Bicd2 mice [13, 14] and the observation that missense mutations in the cargo binding domain of BICD2 impair Rab6 binding with no effect on dynein binding [5], we predicted a loss of function pathomechanism in SMALED2. Following this hypothesis, loss of BICD2 function should result in motor neuron loss. We therefore re-examined the motor neuron phenotype of Bicd2 mice by examining the L4 dorsal and ventral nerve roots in the spinal cord of wild type and knockout mice at postnatal day 21 (p21). We were unable to examine these mice at later time points as knockout mice die by four weeks of age from obstructive hydrocephalus [13]. The L4 nerve root was chosen as it is the largest lumbar nerve root. At p21, there is a significant reduction in the total number of motor axons in Bicd2−/− mice compared to wild type and this is restricted to a subpopulation of motor axons with a diameter of 2.5–4 μm (Fig. 1a-c; Bicd2 mean = 926 ± 22.7 (SEM, n = 7), Bicd2 = 806 ± 25.0 (n = 7), unpaired t-test, p = 0.0041). In contrast, there was no reduction in the number of L4 dorsal root sensory axons between wild type (1778 ± 35.8; n = 6) and Bicd2 (1792 ± 60.4; n = 6) mice (unpaired t-test, p = 0.85) (Fig. 1d-f). The lack of a sensory phenotype was corroborated by further analyses of the percentage of medium-to-large, NF200+ neurons, and small peripherin+ neurons, in the L4 dorsal root ganglia (DRG), which did not reveal any difference between wild type and Bicd2 mice (Supplementary Figure 1).

Fig. 1

Loss of motor axons in Bicd2 mice at 21 days of age. a shows the number of motor axons in the L4 ventral root of Bicd2 (wild type) and Bicd2 (knockout) mice (**unpaired 2-sided t-test, p = 0.0041). b shows a histogram of the total number of L4 motor axons classified using 0.2 μm bins (blue = Bicd2, orange = Bicd2), *p < 0.001; multiple t-tests corrected for multiple comparisons using the Holm-Sidak method. c Representative images from 0.8 μm cross-sections of 1% toluene blue stained L4 ventral roots of Bicd2 and Bicd2 mice. d shows the number of sensory axons in the L4 dorsal root of Bicd2 and Bicd2 mice. e shows a histogram of the total number of L4 sensory axons classified as in b (blue = Bicd2, orange = Bicd2). f Representative images from 0.8 μm cross sections of 1% toluene blue stained L4 dorsal roots of Bicd2 and Bicd2 mice. Error bars = standard error of the mean (SEM). Scale bars = 50 μm

Loss of motor neurons in Bicd2 mice is caused by a muscle non-cell autonomous process

The cerebellar hypoplasia phenotype observed in Bicd2 mice is due to a non-cell autonomous process arising in Bergmann glia [13]. Therefore, we hypothesised that the motor neuron loss we observed in the same knockout mice may also be due to a non-cell autonomous mechanism, but in this instance, arising from muscle tissue, which has previously been shown to have important trophic support functions for motor neurons [23]. According to the neurotrophin hypothesis, during early development, an excess of motor neurons reach their target muscle and compete for muscle-secreted survival factors (e.g. neurotrophins) [24]. A stochastic process of programmed cell death then follows, as those motor neurons that receive insufficient neurotrophins undergo apoptosis. To test this hypothesis, we generated mice in which Bicd2 is selectively knocked out only in skeletal muscle. Homozygous Bicd2 mice, in which the endogenous Bicd2 allele is flanked by intronic loxP sequences [13], were crossed with mice expressing cre-recombinase driven by the endogenous Myod promoter to generate mice lacking both Bicd2 alleles in muscle tissue alone [13, 16]. These mice had a normal life span and no obvious gait abnormality. Analyses of these mice at p21 shows a significant reduction in motor axons in the L4 ventral nerve root compared to wild type mice (Fig. 2a and c; Bicd2 = 928 ± 24 (n = 7), Myod-Cre = 817 ± 39 (n = 7), one-way ANOVA p = 0.019, Dunnett’s t-test, p = 0.034). Furthermore, the subpopulation of motor axons affected is the same as the Bicd2 mouse (2.5–4.4 μm diameter). A small reduction in the number of L4 motor axons between wild type (928 ± 24) and homozygous Bicd2loxP/loxP (888 ± 23, n = 7; Dunnett’s t-test, p = 0.705) mice was detected, but it did not reach statistical significance (Fig. 2a). To confirm that the loss of motor axons in Bicd2 mice was solely due to loss of Bicd2 in muscle, we generated mice lacking Bicd2 in motor neurons. Homozygous Bicd2 mice were crossed with mice expressing cre-recombinase driven by the endogenous ChAT promoter to generate mice lacking both Bicd2 alleles selectively in motor neurons. Analysis of the number of L4 motor axons in these genetically modified mice at p21 revealed no difference compared to Bicd2+/+ mice (Fig. 2a and b; Bicd2 = 928 ± 24; n = 7; ChAT-Cre = 924 ± 23, n = 3; one-way ANOVA p = 0.019, Dunnett’s t-test, p = 0.939).

Fig. 2

Motor axon loss in muscle specific Bicd2 knockout mice. a shows the total number of axons at 21 days of age in Bicd2 mice (blue), Bicd2 mice in which both Bicd2 alleles are flanked by loxP sites (pink), mice with motor neuron specific Bicd2 knockout (ChAT-Cre, red), mice with muscle-specific Bicd2 knockout (Myod-Cre, green), and Bicd2 mice (orange). A significant reduction in the total number of L4 ventral axons was found in Myod-Cre and Bicd2−/− mice compared to wild type (one-way ANOVA p = 0.0187, Dunnett’s 2-sided t-test with Bicd2 as control, **p = 0.019, *p = 0.034). b shows a histogram of the total number of L4 ventral axons at 21 days of age classified using 0.2 μm bins (pink = Bicd2, red = ChAT-Cre). c shows a histogram of the total number of L4 motor axons at 21 days of age classified using 0.2 μm bins (pink = Bicd2, green = Myod-Cre); multiple t-tests corrected for multiple comparisons using the Holm-Sidak method, ***p = 0.0009. d shows a histogram of the total number of L4 motor axons at 42 days of age classified using 0.4 μm bins (pink = Bicd2 (n = 4), green = Myod-Cre (n = 6); multiple t-tests corrected for multiple comparisons using the Holm-Sidak method,****p = 0.0002). Error bars = SEM

A similar analysis performed on the ‘legs at odd angles mouse’ (loa) strain, which bear a missense mutation in Dync1h1 (a binding partner of Bicd2), showed at a later time point of six weeks a reduction in small diameter motor axons, which were assumed to be gamma motor neurons [25, 26]. We therefore repeated this analysis in Myod-Cre mice at six weeks of age. A comparison of the axon diameter cumulative distribution curves shows a loss of motor axons in mice with muscle-specific knockout of Bicd2, revealing a decrease in a subset of motor axons of smaller diameter in the putative range of gamma motor neurons (Fig. 2d). Myod-Cre mice showed no obvious behavioural phenotype which would be in concordance with only a 10% loss of motor neurons.

Active denervation is absent in Bicd2 mice

Next, we sought to determine whether the loss of motor neurons was due to active degeneration of motor axons. Morphological analysis of semi-thin ventral root sections revealed no acute axon degeneration profiles (Fig. 1c). We therefore analysed the neuromuscular junctions (NMJs) from two different muscles (lumbricals and flexor digitorum brevis (FDB) of the hind paw). Wholemount preparations of these two thin muscles were chosen due to reliable quantification of NMJ degeneration without the need for sectioning [19]. Interestingly, all muscles showed a full NMJ innervation pattern, suggesting that at the time of examination, there was no significant denervation (Fig. 3a and c). Furthermore, there was no evidence to indicate that the post-natal developmental process of synapse elimination at the NMJ was affected in Bicd2 mice (Fig. 3b). In contrast, the NMJ area of the hindfoot lumbrical and FDB muscles was significantly smaller in Bicd2 compared to wild type mice (Fig. 3d). The significance of this is not clear but may be a consequence of the significant reduction in overall size of Bicd2 mice and their muscle mass compared to wild type animals (Supplementary Figure 2A and C).

Fig. 3

Normal NMJ analysis of Bicd2 mice at 21 days of age. a shows representative images of the NMJs of the FDB (flexor digitorum brevis) and feet lumbrical muscles stained with anti-SV2/2H3 antibodies (green) to visualise motor neurons and fluorescent alpha-bungarotoxin (red) to identify post-synaptic acetylcholine receptors on the muscle fibre surface. Scale bars = 50 μm. b shows the percentage of fully and partially innervated NMJs in Bicd2 (n = 4) and Bicd2 (n = 4) mice. c shows the percentage of poly-innervated (measure of immaturity) NMJs between Bicd2 (n = 4) and Bicd2 (n = 4) mice. d shows the area of the NMJ (area occupied by each single AchR cluster) in the FDB and lumbrical muscles in Bicd2 (mean 205 and 230 μm2, respectively; n = 4) and Bicd2 (mean 157 and 161 μm2, respectively; n = 4) mice, (multiple t-tests corrected for multiple comparisons using the Holm-Sidak method, *p = 0.05, **p = 0.004). Error bars = SEM. The normal NMJ analysis suggests that there is no active denervation in Bicd2−/− mice at 21 days of age

Muscles in Bicd2 mice display loss of muscle spindles and gamma motor neurons, but no gross morphological changes

Patients with SMALED2 show muscle biopsy abnormalities indicative of both neurogenic and myopathic processes [27]. We therefore examined the gastrocnemius muscle of wild type and Bicd2 mice, but did not identify any consistent differences on H&E-stained muscle or on electron microscopy (Supplementary Figure 3A-D) other than a non-significant reduction in muscle fibre diameter in Bicd2−/− compared to Bicd2+/+mice and a significant reduction in total muscle (and body) mass (Supplementary Figure 2A,C&D).

Motor neuron subtypes may be subdivided according to the fibre type they innervate [24]. For example, ‘fast’ alpha motor neurons innervate types IIa, IIb and IIx muscle fibres, whereas slow motor neurons innervate type 1 muscle fibres. We therefore examined the proportion of muscle fibres of the L4 innervated gastrocnemius muscle in wild type and Bicd2 mice, but we found no difference (Supplementary Figure 2F), indicating that Bicd2 ablation does not induce a significant muscle fibre switch.

As shown in Fig. 2c and d, at both three and six weeks of age, the population of motor axons lost in the L4 ventral roots of mice in which Bicd2 has been selectively knocked out in muscle falls within the 2.5–4 μm range also seen in Bicd2−/− mice. The lack of denervation of NMJs or switch in muscle fibre types, suggests that alpha motor neurons are not affected. At six weeks of age, based on axon diameter, the loss of axons falls within the (Fig. 2d) presumed gamma motor neuron population. This is the same population of motor neurons that is lost in mice with missense mutations in Dync1h1 [26], which model SMALED in humans. The selective loss of gamma motor neurons in these two models suggests a deficit in muscle spindle derived GDNF secretion as GDNF knockout mice show isolated gamma motor neuron loss [28]. Gamma motor neurons represent approximately 30% of motor neurons innervating muscle, but, unlike large calibre alpha motor neurons, they do not form NMJs with skeletal muscle fibres and instead contact muscle spindles. At three weeks, the distinction of motor neuron subtype based on axon diameter is unreliable [24]. To confirm whether the loss of L4 ventral axons at 3 weeks also correlated with a reduction in gamma motor neurons, we examined the total number of muscle spindles in the soleus muscle. The soleus muscle was chosen as it is a small muscle in which the total number of muscle spindles can be accurately quantified on serial transverse sections as previously described [20]. By staining for the SV2 and 2H3 antigens, muscle spindles can be identified, although gamma motor neuron efferents and 1a sensory afferents are indistinguishable. A comparison of the total number of muscle spindles revealed a 20% reduction in Bicd2 mice (9.4 ± 0.68) compared to wild type (12.5 ± 0.29, unpaired t-test p = 0.0065) (Fig. 4b). This correlates well with the reduction in smaller diameter motor axons at six weeks in Myod-Cre mice (Fig. 2d) and confirms that the motor axons lost at three weeks are also gamma motor neurons. The innervation pattern of the remaining muscle spindles was unperturbed, suggesting a developmental, as opposed to a degenerative process. Of note, there was no reduction in large diameter (presumed sensory 1a afferents) in the L4 dorsal nerve root of Bicd2−/− (Fig. 1d-e) and Myod-Cre mice (data not shown). It is possible that SV2 staining is specific to the gamma motor neuron endplate and does not stain for sensory 1a afferents that do not form synapses at the muscle spindle.

Fig. 4

Loss of muscle spindles in Bicd2 mice at 21 days of age. (a) shows example images of a cross section through a muscle spindle in Bicd2 and Bicd2 mice stained for nuclei (DAPI, blue), the neuronal marker SV2/2H3 (green) and laminin (muscle membrane, red). Scale bars = 10 μm. The concentric SV2/SH3 staining around laminin positive muscle fibres indicates a muscle spindle (b) shows the total number of muscle spindles in the soleus muscle of Bicd2 and Bicd2 mice, **p = 0.0065 (unpaired t-test, n = 4–5). (c) shows full innervation patterns in the muscle spindles of Bicd2 and Bicd2 mice. Error bars = SEM. The loss of muscle spindles correlates with the loss of presumed gamma motor neurons in Bicd2 compared to Bicd2 mice

Patient fibroblasts show evidence of impaired secretion

Rab6, a small GTPase interacting with BICD2, is an established regulator of the secretory pathway [29, 30] and controls the flow of secreted proteins that are transported from the Golgi to the plasma membrane via a microtubule-dependent process. As a consequence, loss of Rab6 results in a global reduction of protein secretion [22].

We therefore hypothesised that SMALED2-causing mutations in BICD2 result in an impairment of Rab6 function through its mislocalisation, impairing the targeting of secretory vesicles to the plasma membrane. This would explain why mutations in the N-terminal domain of BICD2, which enhance retrograde processive motility leading to accumulation of BICD2 at the centromere [5, 10], and C-terminal mutations that impair the recruitment of Rab6, but not dynein binding, to BICD2 cause an identical phenotype [31].

To verify this hypothesis, we performed an established VSV-G secretion assay [32] using fibroblasts from a patient with a BICD2I189F SMALED2 mutation and from an age-matched control [31]. These cells were transfected with a plasmid encoding for a temperature-sensitive mutant of vesicular stomatitis virus glycoprotein (ts0–45 VSV-G) tagged with GFP and incubated at 40 °C for 14 h. This high temperature treatment results in misfolding of ts0–45 VSV-G and its retention in the ER. Cooling of the cells to 32 °C allows the refolding of VSV-G and its targeting to the plasma membrane, which could then be quantified. This assay revealed a significant reduction in the rate of secretion of VSVG over time (Fig. 5b and c). Similar to previous experiments performed in Rab6 knockout cell lines [29, 32], secretion was significantly delayed, but not abolished. To confirm that this difference was not simply due to clonal differences, the accumulation of VSV-G at the plasma membrane at 240 min (the time point with the largest difference) was repeated for two additional, unrelated, age-matched controls and SMALED2 patients bearing the S107L and R501P mutations in BICD2. These additional experiments showed a significant reduction in VSV-G appearance on the plasma membrane at 240 min (Fig. 5c; control mean = 0.64 ± 0.017, n = 3; SMALED2 = 0.52 ± 0.043, n = 3; unpaired t-test p = 0.05).

Fig. 5

SMALED2 patient fibroblasts show delayed VSV-G secretion compared to controls. a is an example of a human fibroblast transfected with a plasmid encoding for the thermo-sensitive-GFP vesicular stomatitis virus glycoprotein (ts0–45 VSV-G) at 32 °C. The time prior to fixation is indicated on the top; staining with an anti-VSV-G [8G5F11] against a surface epitope of VSV-G in non-permeabilised cells (top row of panels); GFP staining is shown in the bottom row. Scale bars = 20 μm. At 0 mins, all ts0–45 VSVG is retained within the ER with no surface staining. At 30mins all ts0–45 VSV-G has been trafficked to the Golgi. By 240 mins all ts045-VSV-g has been trafficked to the plasma membrane and is evident in both the GFP (bottom) and anti-VSV-G surface epitope antibody (top)_ panels. b Kinetics of VSV-G secretion in fibroblasts isolated from a patient with SMALED2 (I189F mutation, orange) and an age-matched control (blue) are quantified as the ratio of total surface VSV-G staining to total GFP. The x axis shows the time in minutes at 32 °C prior to fixation (n = 10 cells per condition; **p = 0.008, *p = 0.009; multiple unpaired t-tests corrected for multiple comparisons using the Holm-Sidak method). c shows the average ratio of surface VSV-G to total GFP at 240 min for three independent healthy control and three unrelated SMALED2 (S107L, I189F and R501P) fibroblast cell lines (# p = 0.052, unpaired t-test). Error bars = SEM. The impaired secretion in SMALED2 patient fibroblasts suggests that a similar impairment of secretion may be evident in the muscle of Bicd2 mice and SMALED2 patients

Since BICD2 and Rab6 have also been involved in COPI-independent transport from the Golgi to the ER [33, 34], we investigated the retrograde flow from the Golgi to the ER using galactosyltransferase–GFP in Brefeldin A-treated control and SMALED2 fibroblasts, but did not detect any significant difference between patient and control cells (Supplementary Figure 4A&B).

In light of these results, we propose a non-cell autonomous mechanism of motor neuron loss in SMALED2, which is also likely to be applicable to SMALED caused by mutations in DYNC1H1 (Summarised in Fig. 6). We have shown that loss of BICD2 leads to death of motor neurons and that this is a non-cell autonomous process driven by the loss of BICD2 in muscle. Altogether, our results demonstrate that the BICD2 mouse is a good model of SMALED2 as Bicd2 mice and patients with severe SMALED2 mutations show similar clinical features including hydrocephalus, cerebellar hypoplasia and cortical migration defects [13, 14]. A non-cell autonomous pathomechanism is also supported by a previous study demonstrating that the cerebellar hypoplasia in Bicd2 is driven by non-cell autonomous deficits arising in Bergmann glia [13].

Fig. 6

A summary of the proposed non-cell autonomous mechanism in SMALED2. In normal muscle, BICD2 ensures the correct localisation of Rab6 at the trans Golgi surface, required for the efficient trafficking of secretory vesicles (containing neurotrophins) to the plasma membrane. In SMALED2 muscle, loss of BICD2 function impairs the localisation of Rab6 to the trans Golgi surface resulting in impaired trafficking of secretory vesicles to the plasma membrane and a reduction in neurotrophin release required for motor neuron survival

Quantitative analyses of secretion in patient fibroblasts compared to controls shows only a modest secretion deficit, as found in Rab6 knockout models [22]. How such a delay in secretion causes a significant loss of motor function can be explained in light of the neurotrophin hypothesis in which a surplus of motor neurons compete for muscle-secreted neurotrophic factors for survival [24]. In this model, a reduction in the rate of secretion would yield a decrease in neurotrophin availability, driving an excess of motor neuron death. This hypothesis predicts that no progressive motor neuron degeneration should occur in SMALED2 and muscle-deficient BICD2 mice, as demonstrated in this work and elsewhere [31, 35].

Our model may also explain why the loa, crawling and sprawling mouse models, which bear point mutations in Dync1h1 develop gamma motor neuron and sensory 1a afferent loss [25, 26, 36]. Both gamma motor neurons and 1a sensory afferents innervate muscle spindles and rely on muscle-derived neurotrophic factor secretion during early development. Therefore, an impairment of localised neurotrophin secretion from the muscle of dync1h1 mutant muscles would lead to a loss of gamma motor and sensory 1a neurons.

Our model does not align with previous results generated using a Drosophila model of SMALED2 [37], which led to the proposal that this pathology is determined by a cell autonomous effect of mutant BICD2 in motor neurons [37]. However, Drosophila have only a bicaudal D protein (BicD), unlike mice and humans, which express four bicaudal D proteins (BICD1, BICD2, BICDR1 and BICDR2) with different functions and interacting partners [38, 39]. As such, the model of Martinez-Carrera is expected to be unable to differentiate between the effects of BICD1 and 2, which is crucial, as we have previously shown that BICD1 has an essential cell autonomous role in the sorting of signalling endosomes in motor neurons [40, 41]. Furthermore, Drosophila display important differences to mammals in the neurotrophin pathway at the NMJ [42].

Our model provides strong evidence in support of a non-cell autonomous mechanism of motor neuron loss in SMALED2 which involves impaired secretion of muscle-derived neurotrophins during development. However, a number of important questions still need to be addressed in future studies. The nature of the secreted muscle-derived neurotrophic factor(s) is still unclear as well as whether the phenotypes of mice with loss-of-function mutations in Bicd2 are similar to that of Bicd2 mice. However, the most important point to clarify in future studies is likely to be the precise time point at which motor neuron loss occurs, which is important when considering future therapeutic interventions in SMALED2 patients.

Additional file 1: Supplementary Figures.

Table 1

Primary and secondary antibodies

Epitope	Antibody	Concentration	Species
Type 1 muscle fibres	BA-D5 (DSHB)	1 in 100	Mouse IgG2b
Type IIA muscle fibres	SC-71 (DSHB)	1 in 100	Mouse IgG1
Type IIB muscle fibres	BF-F3 (DSHB)	1 in 100	Mouse IgM
Laminin	L9393 (Sigma-Aldrich)	1 in 100	Rabbit
Mouse IgG2b	AlexaFluor 488 (Life Sciences)	1 in 500	Goat
Mouse IgG1	AlexaFluor 647 (Life Sciences)	1 in 500	Goat
Mouse IgM	AlexaFluor 568 (Life Sciences)	1 in 500	Goat
Rabbit IgG	AlexaFluor 410 (Life Sciences)	1 in 500	Goat
Mouse IgG	AlexaFluor 488 (Life Sciences)	1 in 250	Goat
2H3	DSHB (Supernatant)	1 in 50	Mouse IgG1
SV2	DSHB (Supernatant)	1 in 100	Mouse IgG1
	Tetramethylrhodamine-bungaratoxin (−BTX, Cambridge Bioscience BT00012)	1.5 μg/ml
NF200	N0142, SIGMA	1 in 500	Mouse IgG
Peripherin	Merck Millipore, AB1530	1 in 500	Rabbit IgG
VSV-G	8G5F11, Kerafast	1 in 2000	Mouse IgG

DSHB Developmental Studies Hybridoma Bank