Stage-dependent remodeling of projections to motor cortex in ALS mouse model revealed by a new variant retrograde-AAV9
Abstract
Amyotrophic lateral sclerosis (ALS) is characterized by the progressive degeneration of motoneurons in the primary motor cortex (pMO) and in spinal cord. However, the pathogenic process involves multiple subnetworks in the brain and functional MRI studies demonstrate an increase in functional connectivity in areas connected to pMO despite the ongoing neurodegeneration. The extent and the structural basis of the motor subnetwork remodeling in experimentally tractable models remain unclear. We have developed a new retrograde AAV9 to quantitatively map the projections to pMO in the SOD1(G93A) ALS mouse model. We show an increase in the number of neurons projecting from somatosensory cortex to pMO at presymptomatic stages, followed by an increase in projections from thalamus, auditory cortex and contralateral MO (inputs from 20 other structures remains unchanged) as disease advances. The stage- and structure-dependent remodeling of projection to pMO in ALS may provide insights into the hyperconnectivity observed in ALS patients.
https://doi.org/10.7554/eLife.36892.001Introduction
Amyotrophic lateral sclerosis (ALS) is classically described as a disease of upper and lower motor neurons (Ravits et al., 2007). However, postmortem studies have demonstrated that pathological hallmarks of ALS (such as TDP-43 inclusions) appear in motor cortex as well as in multiple cortical and subcortical structures during disease progression, depicting a propagation pattern (Brettschneider et al., 2013; Braak et al., 2013). Thus, ALS has been growingly conceptualized as a disease of the motor subnetwork or, in more advanced stages, a multi-network disease affecting motor, premotor and sensory areas. To which extent these networks are damaged by the neurodegeneration and remodel themselves during the disease progression, it is a subject of active investigation.
Diffusion-tensor imaging (DTI) studies have revealed that the degeneration of axons is detected in more and more white-matter tracts affecting different structures (corticospinal, corticorubral, corticostriatal and perforant tracts) at different stages of the propagation scheme (Neuroimaging Society in ALS (NiSALS) DTI Study Group et al., 2016; Kassubek et al., 2014; Kassubek et al., 2018) as well as in interhemispheric connections (Filippini et al., 2010) and in spinal cord (Cohen-Adad et al., 2013). However, analysis of resting-state networks by functional magnetic resonance imaging (fMRI) has revealed increased functional connectivity has been described in symptomatic ALS patients in the sensorimotor (Menke et al., 2016; Agosta et al., 2018) somatosensory (Agosta et al., 2011) and cortico-striatal networks (Fekete et al., 2013). Notably, an increase in functional connectivity in the motor subnetwork was already detectable in presymptomatic and early-stage ALS patients (Schulthess et al., 2016).
Although these data indicate that a significant architectural remodeling of the motor subnetwork takes place in ALS patients in the face of ongoing neurodegeneration, the extent of such remodeling and the structural changes underlying this effect in the motor network are unclear. Of particular interest is the extent of remodeling of the input network to pMO: in fact, the type and amount of projection to pMO determines not only the output from MO, which is strongly related to the execution of movements; in addition, the synaptic inputs are strong regulators of neuronal activity pattern and may influence activity-dependent pathogenic factors (Roselli and Caroni, 2015; Bading, 2017). To date, connectivity maps have been obtained at mesoscale (Oh et al., 2014; Mitra, 2014) or at cellular resolution (Mao et al., 2011; Wertz et al., 2015) only for wild-type models.
In order to address the remodeling of the subnetwork in an experimentally tractable condition, we have selected to quantitatively map the projections to the primary motor cortex network in the SOD1(G93A) ALS mouse model using a newly developed AAV variant endowed with retrograde tracing abilities. We have demonstrated that selected cortical areas display an increased projection to pMO already at presymptomatic stage and the increase in projection expands during disease progression to involve additional cortical and subcortical structures. We show that the pattern observed in the mouse model bears significant similarities with fMRI functional connectivity data gathered from human ALS patients. Thus, we have identified one structural component of the early and selective remodeling of the large-scale architecture of the motor subnetwork, which may contribute to explain the functional connectivity changes observed in patients.
Results
AAV9-SLR variant is a new, efficient viral tool for retrograde connectivity tracing
In order to identify a variant of AAV9 suitable for retrograde tracing, we considered three variants previously generated (AAV9-SLR, AAV9-NSS and AAV9-RGD) as part of an ongoing screening of the properties of newly generated AAV9 variants. These variants displayed a high infection efficiency in vitro in multiple cell lines (human coronary artery endothelial cells (HCAECs), human coronary artery smooth muscle cells (HCASMC), HEK293T, HeLa, 911, HepG2; Varadi et al., 2012), suggesting a potential for infecting neurons and neuronal processes with high efficiency; in particular, we evaluated them for highly desired applications such as high local infectivity and retrograde infectivity. We first tested infection rate of these AAV9 variants in vivo (Figure 1A). Five independent groups of tdTomato-ROSA26 reporter mice (N = 3 for each group, age P20) were injected in the dorsal striatum (DS) with 500 nl of viral suspension (9*1013 genomes/ml) of AAV9-SLR, AAV9-NSS, AAV9-RGD and, for comparison, WT-AAV9 and WT-AAV2 (Figure 1B). The total number of infected neurons in the injection area (local infectivity) and in regions projecting to DS was assessed at 15 days post-injection (DPI). One-way ANOVA revealed a significant difference (F(4,10)=198.8, p<0.0001) in the number of infected neurons per area unit (infection rate) in the injection ROI. Compared to WT-AAV9, infection rate was significantly higher for AAV9-SLR (post-hoc p=0.0001) but not for AAV-NSS, (p=0.0910), and was significantly lower for both AAV-RGD (p=0.0001) and WT-AAV2 (p=0.0001; detailed statistics are reported in Supplementary file 1a). In agreement with the limited toxicity of AAVs, no morphological sign of distress (axonal beading, dendrites fragmentation, fluorescent fragments of cell bodies) were identified in neurons infected with the AAV variants or with the WT serotypes.
-
Figure 1—source data 1
- https://doi.org/10.7554/eLife.36892.003
When screened for their ability to retrogradely infect neurons projecting to the injection site (in DS), the AAV variants differed significantly from WT-AAV9 in terms of retrogradely infected neurons located in substantia nigra (SNc; one-way ANOVA F(4,10)=135.1, p>0.0001; detailed statistics are reported in Supplementary File 1b) and motor cortex (MO; F(4,10)=46.5, p<0.0001; Supplementary File 1c). Injection of AAV9-SLR and AAV9-NSS resulted in a larger number of tdTomato+/tyrosine-hydroxylase-positive (TH+) neurons in SNc than WT-AAV9 (AAV9-SLR: p=0.0001; AAV9-NSS: p=0.0001; Figure 1C), whereas injection of WT-AAV2 and AAV-RGD resulted in no retrograde infection to SNc. Increased retrograde infection from DS to MO (Figure 1D) was detected for AAV9-SLR (p=0.0001) but not for AAV9-NSS (p=0.0909), whereas WT-AAV2 (p=0.1153) and AAV9-RGD (p=0.0018) had a significantly lower retrograde infection rate. Finally, we evaluated the retrograde infection rate from visual cortex (V1) to lateral geniculate nucleus (LGN, Figure 1E). Injection of 500 nl of suspension of each of the five AAV variants in V1 of reporter tdTomato/ROSA26 (N = 3) resulted in a significantly different number of tdTomato + neurons in LGN (one-way ANOVA, F(4,10)= 17.8p=0.0002; Figure 1F). AAV9-SLR and AAV9-NSS infected a significantly larger number (up to 10-fold) of LGN neurons than WT-AAV9 (AAV9-SLR p=0.0005; AAV9-NSS: p=0.0446; detailed statistics are reported in Supplementary file 1d), whereas, neither WT-AAV2 (p=0.7941) nor AAV9-RGD (p=0.7271) resulted in any infection in the LGN. Taken together, these data identified AAV9-SLR as a new AAV variant endowed with robust, reproducible and broad retrograde infection capability.
Mapping input to primary motor cortex by AAV9-SLR
We then exploited the AAV9-SLR for mapping the forebrain and subcortical neurons projecting to primary (pMO, Figure 2A) or secondary (sMO) motor cortex in WT mice. Retrogradely labeled neurons in the cortical and subcortical structures of the forebrain were quantified, anatomically annotated and normalized for the volume of the injection site (Figure 2B,C). Overall, we identified 28 distinct anatomical structures providing direct input to the injection site in pMO, including thalamic nuclei (TH; Figure 2D), followed by homolateral somatosensory cortex (SS; Figure 2D), contralateral secondary motor cortex (cMOs) and contralateral primary motor cortex (cMOp; Figure 2D). Values for raw neuronal counts (before normalization) are displayed in Figure 2—figure supplement 1A, detailed statistic provided in Supplementary file 1e. Detailed statistic after normalization for injection site volume are in Supplementary file 1f. When projecting neurons are expressed in terms of percentage of the total population of projecting neurons (Figure 2—figure supplement 1B), TH represented 32 ± 7% of the total input to pMO, SS 22 ± 7%, cMOs 8 ± 3% and cMOp 6 ± 2% (full details are provided for all structures in Supplementary file 1g). Additional projections were provided by auditory cortex (AUD 6 ± 1%; Figure 2D) and contralateral somatosensory cortex (cSS 5 ± 5%) and, among subcortical structures, by caudoputamen (CP 6 ± 2%), claustrum (CLA 3 ± 1%) and basolateral amygdala (BLA 0.4 ± 0.3%). We consistently identified a small contingent of neurons projecting to pMO located in hypothalamus (HY 0.6 ± 0.2% of the total pool of neurons projecting to pMO in the forebrain); out of these hypothalamic neurons, not previously reported, the majority was located in the lateral zone of hypothalamus (zona incerta 23 ± 10% of the total of hypothalamic neurons projecting to pMO; lateral preoptic area 23 ± 19% of the total of hypothalamic neurons projecting to pMO; lateral hypothalamic area 14 ± 18% of the total of hypothalamic neurons projecting to pMO). In addition to previously reported input structures to pMO (Hooks et al., 2013; Oh et al., 2014; Mao et al., 2011), we identified additional inputs from paraventricular hypothalamic area (PAL 0.8 ± 0.2) and gustatory area (GU 0.7 ± 0.3%).
-
Figure 2—source data 1
- https://doi.org/10.7554/eLife.36892.009
The input to the neighboring sMO was qualitatively similar to pMO, including all of the structures projecting to pMO (all in all, 20 distinct regions were found to project to sMO compared to 28 for pMO). Values for raw neuronal counts (before normalization) are displayed in Figure 2—figure supplement 1C, detailed statistic provided in Supplementary file 1e). The trend remained unchanged after normalization for injection site volume (Figure 2E–G; detailed statistic Supplementary file 1f). When projecting neurons are expressed in terms of percentage of the total population of projecting neurons (Supplementary file 1, Figure 2D), SS emerged as contributing a significantly larger share of input to sMO than to pMO (33 ± 3%, p<0.0001 vs pMO), followed by structures whose input to pMO and sMO was comparable (detailed values reported in Supplementary file 1g): TH (28 ± 1%, p=0.5965), AUD (10 ± 2%, p=0.1914), cMOs (7 ± 3%, p>0.9999), cMOp (4 ± 1%, p=0.9815), orbital area (ORB: 2 ± 0.4%, p>0.9999) and BLA (2 ± 0.7%, p>0.9999).
Since the thalamus appeared to provide a significant fraction of the projections to pMO and sMO, we annotated the labelled neurons to thalamic nuclei (raw and normalized neuronal counts are reported in Supplementary file 1h and 1i, respectively). In total, 19 thalamic nuclei were found to project to either pMO or sMO motor cortex (Figure 2—figure supplement 2A–F). When considered as a fraction of the total number of thalamic neurons projecting to pMO or sMO, the ventral anterior-lateral complex and the ventral medial nucleus appeared to be the most relevant source of projection to MO, and their contribution to pMO or sMO was comparable (VAL to pMO 29 ± 6%, to sMO 34 ± 6%, two-way ANOVA p=0.6496) and ventral medial nucleus (VM to pMO 24 ± 7%, to sMO 25 ± 12%, p>0.9999; detailed values reported in Supplementary file 1j).
Increased input to pMO in early presymptomatic mSOD
We then exploited the retrograde AAV9-SLR to quantitatively explore the projections to pMO in mSOD mice and WT littermates during disease progression. First, we focused on early-presymptomatic stages: AAV9-SLR was injected in pMO at P20 and animals were sacrificed at P35 (Figure 3A). Retrogradely labeled neurons in the cortical and subcortical structures of the forebrain were quantified, anatomically annotated and the raw counts normalized according to the volume of the injection site (Figure 3B,C). Strikingly, we found an overall increase in the number of neurons projecting to pMO in mSOD compared to WT both in raw counts (Supplementary file 2a and Figure 3—figure supplement 1A) and after normalization (Two-way ANOVA, brain regions F(26,189)=53.0 p<0.0001, genotype F(1,189)=16.3 p<0.0001). Post-hoc analysis (Figure 3B; detailed numerical values in Supplementary file 2b) revealed a statistically significant increase in the mean number of neurons projecting to pMO from SS (effect size: 1.6; p<0.0001) and from cMOs (effect size: 1.6; p=0.0189). Notably, among the 28 anatomical regions considered, all others displayed a comparable number of neurons projecting to pMO in WT and mSOD (p>0.05 WT vs mSOD for all the brain regions). For all regions considered, the median value of neurons projecting to pMO was close to the mean value, and the skewness value never exceeded 1.7 (with 16/28 regions in WT and 28/28 regions in mSOD having values ranging between −1.0 and +1.0), discounting the possibility of strongly biased or non-normal distributions.
-
Figure 3—source data 1
- https://doi.org/10.7554/eLife.36892.013
We further investigated the pattern of neurons projecting to pMO from structures with increased projection (SS, cMOs) and unaltered projection (cMOp, AUD) annotating them by cortical layers (detailed values before and after normalization for injection site volume and statistical analysis reported in Supplementary file 2c and d). Interestingly, despite the overall increase in the number of neurons projecting to pMO, a significant loss of projecting neurons was detected in SS layer V of mSOD mice (effect size: −1.5; p=0.0324, Supplementary file 2a) mirrored by a reciprocal trend toward increase projection from layer II/III (effect size: 1.2; p=0.2734) but not in layer VI (p=0.9898). (Figure 3D; values expressed in terms of percentage of the total population of projecting neurons, detailed statistic Supplementary file 2e). On the other hand, the distribution of projecting neurons across cortical layers was comparable in WT and mSOD for cMOs, cMOp and AUD (data before normalization displayed in Figure 3—figure supplement 1B).
To rule out any infectivity bias of the virus in WT vs mSOD mice, we verified that the selective increase in projection from SS to pMO could be detected by an independent method. We injected pMO with fluorescently labeled cholera toxin B (CTb) and assessed the number of CTb-labeled neurons in SS and cMOp (showing increased or unchanged projection to pMO, respectively, in the viral tracing experiment); neuronal counts were normalized for injection volume (Figure 4A–C; counts before normalization are depicted in Figure 4—figure supplement 1A, and numerical values in Supplementary file 2f). Confirming the viral tracing data, CTb labeled a larger number of neurons in the SS in the mSOD compared to WT (both N = 3; one-way ANOVA F(3,8)=17.0, p=0.0008), whereas no difference was detected in cMOp between the two genotypes (SS effect size: 1.4; p=0.0156, cMOp effect size: −0.4; p=0.9971; Figure 4A–C and Supplementary file 2g).
-
Figure 4—source data 1
- https://doi.org/10.7554/eLife.36892.016
We further sought to rule out any differential infectivity of the virus in mSOD model by quantifying the neurons projecting to primary visual cortex (V1) from LGN and MO in P20 mice (Figure 4D–F). Injection of AAV9-SLR in V1 resulted in a comparable number of labelled neurons (either before either after normalization for the injection site volume: Supplementary file 2h and 2i) in WT (N = 3) and mSOD (N = 3) both in LGN (p=0.9263) and in MO (p=0.6360). Counts before normalization are depicted in Figure 4—figure supplement 1B.
In contrast to our observation in the pMO, input to sMO was comparable in WT (N = 3) and mSOD (N = 4) mice (either before either after normalization; for normalized two- way ANOVA, genotype F(1,100)=0.7572 p=0.3863. Figure 5A–C and Supplementary file 2j and 2k for non-normalized and normalized counts, respectively). After normalization, a trend toward increased projection to sMO in mSOD was detected only in cMOs (effect size: 1.3; p=0.0708). In conclusion, retrograde tracing of the input to pMO revealed an unexpected increase in the number of neurons projecting from SS and cMOs to pMO as early as P20. This reveals that remodeling of cortical connectivity is already initiated in presymptomatic mSOD mice.
-
Figure 5—source data 1
- https://doi.org/10.7554/eLife.36892.018
Cortical and subcortical structures display an increased input to pMO in mSOD during disease progression
We then investigated how the projections to pMO changed during disease progression. We compared the number and anatomical locations of neurons projecting to pMO in mSOD and WT mice at the early-symptomatic age, P60 (WT N = 3, mSOD N = 4) and symptomatic age, P90 (WT N = 5, mSOD N = 5; Figure 6A, B); investigation of later timepoints was prevented by the observation, in a preliminary study, of high intraoperative mortality (3 out of 3 mSOD mice, possibly due to respiratory failure) in mice injected at P115. The number of neurons projecting to pMO from the different structures was expressed as fold-change of the age-match WT (detailed values reported in Supplementary file 3a). Distribution of the projections to pMO was stable over time in WT mice (two-way ANOVA; brain regions F(19,209)=0.001 p>0.99). However, when compared to their age-matched WT littermates at P60, mSOD mice displayed an increase in neurons projecting to pMO (two-way ANOVA; brain regions F(19,100)=2.3 p=0.0034) from SS (effect size: 1.6; p=0.0019; unchanged compared to P20 mSOD, p>0.9999) and cMOs (effect size: 1.6; p=0.0454; unchanged compared to P20 mSOD, p>0.9999). In addition, now an increased number of neurons projecting to pMO was found in cMOp (effect size: 1.4; p=0.0077 in mSOD vs WT littermates,+2.4 ± 0.9 vs P20 mSOD, p=0.0110). Besides these three, no other structure, among the 28 analyzed, showed a significant difference between WT and mSOD (Figure 5A–B and Supplementary file 3a). Investigation of the symptomatic stage (P90) revealed a further, significant modification of the abnormal pattern of projections to pMO (two-way ANOVA; brain regions F(19, 160)=7.6, p<0.0001), affecting both cortical and subcortical regions. Projections to pMO from SS was still increased (effect size:1.6; p=0.0095; same increase registered for P20 and P60 mSOD, p>0.9999 for both), cMOp (effect size: 1.7; p<0.0001, also increased when compared to the P60 timepoint:+1.6 ± 0.5 vs P60 mSOD, p=0.0124) and cMOs (Effect size: 1.4; p<0.0001, comparable to the increase registered in mSOD at P20 (p=0.7658) and P60 (p=0.3917)). Moreover, increased number of projecting neurons was detected from TH (effect size: 1.7; p=0.0014;+1.8 ± 0.4 vs P20 mSOD, p=0.0064), AUD (effect size: 1.6; p<0.0001;+1.4 ± 0.5 vs P20 mSOD, p<0.0001), and CP (effect size: 1.7; p<0.0001;+2.0 ± 0.5 vs P20 mSOD, p=0.0086). As for the P20 timepoint, median values were close to the mean value of relative change in the number of projecting neurons and skewness ranged between −1.0 and +1.0, indicating the absence of biased distributions of the data.
-
Figure 6—source data 1
- https://doi.org/10.7554/eLife.36892.022
To exclude an effect of a widespread neuronal loss in the pMO of mSOD on the number of axons projecting to pMO, we performed a NeuN staining in the pMO of P90 mSOD mice (n = 3) revealing no difference in the total number of neurons in our injection area (p=0.4987, Figure 6—figure supplement 1A).
Taken together, these findings revealed a stage- and structure-dependent remodeling of input to pMO, characterized by the progressive increase in the number of neurons projecting to the pMO, with cortical areas affected at earlier stages and subcortical areas only at later stages.
Pyramidal layer V neurons across multiple areas of the motor subnetwork display simultaneous loss of spines on basal dendrites
Having demonstrated an increase in projection to pMO from specific cortical and subcortical structures with a progression-dependent pattern, we sought to determine if this phenotype could be correlated with pathogenic processes ongoing in the projecting neurons themselves. In order to verify if the neurons projecting to pMO already displayed signs of ALS-related pathogenic process, we elected to use dendritic spine density as readout. In fact, dendritic spine loss has been reported to be one of the earliest morphological signs of pathological involvement of several cortical and subcortical neuronal subpopulations in ALS (Ozdinler et al., 2011; Fogarty et al., 2016; Fogarty et al., 2017). Pyramidal neurons located in layers II/III and V were acquired from: SS (whose projection to pMO was increased in early stage of the disease), cMOp (increased projection starting from intermediate stage of the disease) and AUD (increased projection only at the later stages). We also included a sample of layer V and layer II/III neurons from ipsilateral pMO (known to be affected by the disease already at early stages; Fogarty et al., 2017).
At early pre-symptomatic stage (when increased projections from SS, but not from AUD, are already detectable), overall spine density was comparable in neurons in SS, cMO and AUD projecting to pMO from WT and mSOD mice (one-way ANOVA, F(13,108)=0.8, p=0.6317; Figure 7A and Supplementary file 3b). Interestingly, layer-V neurons projecting to the pMO in SS, AUD and cMO displayed a decrease in spine density from P60 (Figure 7B), (F(13,119)=7.8, p<0.0001); at the same time- point, a decrease in dendritic spine density in neurons from layer V of pMO itself became evident (p=0.0071; in agreement with what previously reported; Ozdinler et al., 2011; detailed numerical values are reported in Supplementary file 3c). Unlike layer V, we found no statistically significant change in the spine density in layer II/III at P60 in any of the cortical areas considered. At P90, a significant decrease in spine density affected neurons projecting to pMO from all structures under scrutiny either in layer II/III either in layer V (effect size between −1.6 and −1.8; F(13,128)=47.8, p<0.0001, Figure 7C and Supplementary file 3d).
-
Figure 7—source data 1
- https://doi.org/10.7554/eLife.36892.024
In conclusion, the loss of spines did not temporally correlate with the increased-projections phenotype, suggesting that the latter may not be interpreted as a consequence of pathology in projecting neurons (within the limits of the spine-loss readout). In addition, we show that spine loss affects not only vulnerable neurons in pMO, but rather neurons across the whole motor subnetwork.
misfSOD expression is not related to network remodeling nor to spine pathology
Next, we investigated if the build-up of misfolded SOD1 (misfSOD1) in neurons projecting to pMO may be related to the increased-projection phenotype. MisfSOD1 has been reported to be a hallmark of mSOD-ALS pathobiochemistry (Bosco et al., 2010) and has been identified as the earliest disease marker in vulnerable spinal MN (Saxena et al., 2013). To this aim, mice were injected with fluorescently labeled CTb in pMO at P20 and sacrificed at P35, in order to label the neuronal subpopulations projecting to pMO; misfSOD1 was revealed by the B8H10 monoclonal antibody against one conformational epitope of misfolded SOD (Pickles et al., 2013). In pMO, where upper MN are mainly (but not exclusively) localized, already at this age 41 ± 7% of neurons in pMO layer V were positive for misfSOD +compared to only 1 ± 0.5% of layer II/III neurons. MisfSOD1 burden (assessed by the immunofluorescence intenty) in misfSOD1 +neurons was comparable in pMO layer II/III (1798±355) and layer V (2189 ± 393; p>0.05). We then identified the neurons in SS and AUD projecting to pMO by CTb retrograde labelling (the non-viral approach was used to avoid possible artifacts due to the over-expression of ZSGreen) and assessed the mSOD burden by immunostaining (Figure 8A–C). Both SS and AUD included a misfSOD +subset of neurons in layer II/III (20 ± 5%) and layer V (54 ± 6%). However, when we considered the CTb + and CTb- subpopulations of neurons, we found that, in layers II/III and V the population of neurons projecting to pMO (CTb+) displayed a significantly lower burden of misfSOD1 than the overall neuronal population of the CTb- neurons (effect size between −0.9 and −1.4; one-way ANOVA, F(11,1242) = 35.3, p<0.0001; Figure 8A–C and detailed numerical values in Supplementary file 3e) in both cortical layers of SS (p<0.0001) as well as in both cortical layers of AUD (p<0.0001). In order to verify that this finding was not due to a confounding effect of CTb itself, we traced a different subpopulation of SS neurons, namely those involved in inter-hemispheric projections; we injected fluorescently labeled CTb in contralateral SS of mSOD1 mice at P30 and we assessed the burden of misf SOD in CTb + and CTb- neurons (Figure 8D; E). In contrast to what observed in the subpopulation connected to pMO, SS subpopulation projecting to contralateral SS displayed a range of misfSOD levels comparable to non-projecting, nearby neurons (p=0.2696) including, in both cases, neurons with high burden and neurons with low burden of mSOD (Figure 8D–E and Supplementary file 3f). Thus, taken together, these findings suggest that although within regions part of the motor subnetwork there is a mixed population of neurons with high and low burden of misfSOD, the neurons directly projecting to pMO constitute a subpopulation with comparatively low misfSOD burden.
-
Figure 8—source data 1
- https://doi.org/10.7554/eLife.36892.026
Next, we asked if the spine loss phenotype observed in neurons projecting to pMO at later timepoints could be correlated with misfSOD build-up. To this aim, P60 mSOD/ZsGreen mice were injected with rAA9-SLR in pMO (N = 3) and the burden of misfSOD was assessed in SS ZsGreen + neurons (layer II/III and layer V) by immunostaining with the B8H10 monoclonal antibody (Figure 8F; G). Based on misfSOD intensity thresholds (the 10th and 90th percentile values of the single-neurons distribution shown in Figure 8B), we discriminated two population of ZsGreen + neurons (high and low misfSOD expression, either in layer II/III either in layer V) for which the spine density on basal dendrites was counted. Interestingly, despite the extreme difference in misfSOD burden, spine density was comparable in two groups and was within the range of WT mice for layer II/III neurons and reduced (compared to WT) for layer V neurons (in agreement with the values in Figure 7B; statistic in Supplementary file 3g), suggesting that spine loss may be a non-cell-autonomous phenomenon not directly related to the burden of mistSOD.
Pattern of functional connectivity alterations in ALS patients displays similarity to the changes in projections in ALS mouse model
Having identified a subgroup of cortical areas and subcortical structures whose projection to pMO was increased in mSOD mice, we set out to verify if a similar set of areas display increased connectivity to pMO in ALs patients. To this aim, we analyzed the connectivity pattern in a large cohort of 71 ALS patients compared to 28 healthy controls by using ‘resting-state’ fMRI. The connectivity patterns were studied using a connectome-based approach in order to investigate region-to-region functional connectivity in a manner allowing comparison of results between the patient and mSOD mouse studies. Each of the a priori defined seed regions allowed for the identification of well-known functional brain networks in healthy humans (Smith et al., 2009) (Figure 9A), which indicated robust seed locations for the functional connectome analysis. Overall functional connectivity analysis demonstrated no statistically significant differences between groups (t = −0.51, p=0.613), which rules out a potential bias in the overall level of connectivity (van den Heuvel et al., 2017). However, network-based statistics of functional connectome data revealed significantly increased functional connectivity (hyper-connectivity) between regions frontal to the pMO, namely the dorsolateral prefrontal association cortex and the pMO itself, both homolaterally and contralaterally (t ≥ 2.81, p≤0.0065, corrected; Figure 9B). Functional connectivity between left and right primary somatosensory cortex was decreased (t = 2.93, p=0.0048) and no effect was detected between thalamic nuclei and pMO or in visual cortex. Thus, the fMRI data from ALS patients confirmed the occurrence of a large-scale remodeling in the motor subnetwork (as observed in the murine dataset and confirming previously published data; Schulthess et al., 2016) and showed a predominant increase in connectivity between pMO and both homolateral and contralateral frontal-lobe areas, showing similarity with the increase projection from contralateral rostral structures, such as sMO, observed in the mouse model.
Discussion
In the present work, we have demonstrated, using viral-tracing approaches, that the motor subnetwork in an ALS mouse model undergoes distinct remodeling during disease progression with stage and structure-dependent effects: projection from SS is increased already in presymptomatic mice and remains stable, whereas projections from AUD and cMOs and cMOp either appear or increase further in the symptomatic stage. We show that this pattern displays similarities with the network remodeling in the form of increase functional connectivity changes observed by fMRI in ALS patients. Thus, we provide cellular-resolution data that contribute to analysis of large-scale architecture of motor circuits in ALS mouse model, which may provide entry points to understand changes occurring in human patients.
In order to map the input to the MO in mSOD and WT mice, a retrograde viral-tracing strategy was considered. We have introduced a new variant of AAV9 which displays significantly higher capacity for retrograde infection of neurons in vivo. Although retrograde tracing can be achieved with non-viral approaches, these methods suffer from limited sensitivity, stability and do not provide genetic access to the target neurons (Tervo et al., 2016). Retrograde tracing has been previously achieved in vivo using rabies-derived vectors (Wickersham et al., 2007), herpes simplex virus (Ugolini et al., 1987) or canine adenovirus (Soudais et al., 2001). These vectors, however, have been plagued by the intrinsic toxicity of the expressed viral proteins (as in the case of rabies - Ginger et al., 2013 - although new rabies variants may have overcome this issue, see Ciabatti et al., 2017) or by the tendency to activate inflammation (Muruve, 2004). Recently, the development of retrograde AAV vectors have been undertaken to allow retrograde tracing with genetic access and low toxicity; in particular, a retrograde AAV2 variant has been reported (rAAV2; Tervo et al., 2016). Our AAV9 variant extends and complements the existing toolset of AAV vectors endowed with retrograde infectivity. Although rAAV2 belongs to the AAV2 serotype and may strongly bind heparan sulphate, our AAV9-derived retrograde AAV variant has the advantage of reduced binding to heparin sulfate (Shen et al., 2011) and therefore it is less affected by changes in brain extracellular matrix (e.g. due to neuroinflammation associated with neurodegeneration); AAV9 is only infrequently targeted by pre-existing immunity to AAVs (Zhang et al., 2011; Saraiva et al., 2016) and may provide the entry point to enhanced vectors for circuit-specific gene therapy in humans since AAV9 has been already successfully tested for this application (Boutin et al., 2010).
We first validated the tool by mapping the input connectivity to pMO and sMO at cellular resolution; notably, the discreet pattern of regions containing retrogradely labelled neurons separated by others containing or no labeled cell, in a reproducible pattern and matching the list of brain structures projecting to pMO (e.g. Hooks et al., 2013) and the Allen Brain Connectivity atlas) suggest that the infection of distant neurons is not an artifact of passive viral diffusion or infection via blood supply (which would have resulted in a continuous or erratic distribution of labelled neurons). We further validated the reliability of our approach in the context of the mSOD1 mouse, verifying that viral-tracing results on SS and pMO could be reproduced by non-viral tracing (CTb) and that no abnormalities could be found in the projections to visual cortex, an area spared the pathogenic process. Nevertheless, AAV vectors relie on retrograde transport to reach the cell body (where expression of the genome takes place; Castle et al., 2014b) in process dependent on dynein/dynactin (Castle et al., 2014a) and one may speculate that the transport slow-down and dysfunction reported in ALS (Bilsland et al., 2010; Marinkovic et al., 2012) may bias the quantitative evaluation of the tracing experiments. However, the slowing of the transport would result in the decrease or delayed arrival of the AAV to the neuronal cell body, actually decreasing the chances of Cre expression and possibly resulting in decreased number of projecting neurons; therefore, this mechanism is unlikely to underlie the observed increase in the number of projecting neurons.
We exploited this new retrograde AAV9 to investigate in ALS models the qualitative and quantitative changes in the structure of the motor subnetwork, in particular of the projections to the pMO. A major finding of the present work is the increase in the numbers of neurons projecting from several cortical and subcortical regions to pMO; the increase is selective for some cortical areas only and the number of regions involved increases over time. The increase in the number of neurons projecting to a fixed, small volume of pMO is most likely to be interpreted as due to the expansion of the axonal arborizations of these neurons in pMO as disease progresses and it is in agreement with the increased frequency of EPSPs observed early in disease course (Fogarty et al., 2015). What may be the driving force of this expansion? One may hypothesize that a pathogenic process ongoing in projecting neurons may drive the expansion; however, neither the burden of misfolded SOD1 nor the signs of structural disturbances (i.e. the decrease in spine density) were detectable in the areas displaying increased projection to pMO. On the other hand, it has been demonstrated in several models (e.g. in the neuromuscular unit), that silencing of either the pre- or post-synaptic site is sufficient to induce axonal sprouting (Caroni et al., 1994; Caroni, 1997). Since reduced firing of vulnerable spinal MN has been shown to precede degeneration (Martínez-Silva et al., 2018; Roselli and Caroni, 2015), it is tempting to speculate that the expansion of the projections may be driven by abnormal activity of vulnerable post-synaptic neurons in the pMO (including vulnerable corticospinal MN; Ozdinler et al., 2011). Indeed, in mSOD mice excitability changes appear to affect a larger set of neurons in pMO than previously thought (Kim et al., 2017).
When the progressive increase in projections to pMO detected in the tracing experiments are compared with the changes in connectivity reported here and by others (Agosta et al., 2018; Menke et al., 2016; Schulthess et al., 2016), several similarities emerge: (i) both in human and mouse data, we observe a progressive remodeling of the motor subnetwork, (ii) the remodeling of the motor subnetwork is detectable before the appearance of overt clinical or behavioral abnormalities, (iii) the remodeling of the motor subnetwork involves the primary motor cortex and more rostral, homolateral or contralateral structures (in the frontal lobe, prefrontal cortex or secondary motor cortex in human and mouse, respectively) and projecting to pMO, and (iv) both in human and mouse data the primary somatosensory cortex appears involved in the motor subnetwork remodeling (either in terms of increased projection to pMO in the mouse, or in terms of connectivity between left and right SS, SS, sensory-motor or SS to other cortical areas; Douaud et al., 2011; Agosta et al., 2011). Nevertheless, thalamocortical connectivity does not appear to be altered in the functional connectivity studies.
Based on these similarities, might the increase in projection to pMO be a leading mechanism underlying the functional hyperconnectivity in ALS patients? Although the relationship between structural and functional connectivity is not precise, severals studies have demonstrated that, at large-scale level, structural connectivity is a strong predictor of functional connectivity (Honey et al., 2009; Honey et al., 2007; Damoiseaux and Greicius, 2009; Hermundstad et al., 2013) and maps of functional and structural connectivity display high degrees of similarity (van den Heuvel et al., 2009). Moreover, studies in non-human primates comparing functional connectivity with anatomical tracing have revealed that anatomical connectivity strongly (although not exclusively) shape the functional connectivity (Shen et al., 2012; Shen et al., 2015). On the other hand, deducing structural connectivity from functional connectivity data can be fraught with inaccuracies, since indirect connections may contribute to it (Adachi et al., 2013). Thus, on the assumption that the observations made on the remodelling of the motor subnetwork in mouse models could be extrapolated to humans, the expansion of projections to pMO may contribute a structural substrate of the functional hyperconnectivity observed in fMRI studies. Notably, alternative models have been proposed to explain the hyperconnectivity phenotype. Changes in functional connectivity could also result from an alteration in the local excitation/inhibition balance due to the dysfunction of inhibitory interneurons (Douaud et al., 2011): in fact, pathological studies have found that calbindin-positive GABAergic interneurons undergo degeneration in ALS (Maekawa et al., 2004), lending support to the hypothesis that loss of inhibitory tone contributes to clinical dysfunction in ALS (Turner and Kiernan, 2012). The two models are actually not mutually exclusive: the pattern of either increased either decreased connectivity in resting-state (Jelsone-Swain et al., 2010; Mohammadi et al., 2009; Zhou et al., 2013; Agosta et al., 2013; Douaud et al., 2011) and task-based fMRI (Heimrath et al., 2014; Konrad et al., 2002; Lulé et al., 2009) may be the result of ongoing degeneration and loss of axons (Dadon-Nachum et al., 2011; Sarica et al., 2014) accompanied by potentially compensative phenomenon (Bernier et al., 2017). The observed increase in projection from cortical and subcortical areas might be an early adaptive response to ongoing cell dysfunction up to the point where a critical cell loss is reached and a disconnection syndrome begins to emerge (Gorges et al., 2017).
In addition to the study of large-scale patterns of projections to vulnerable cortical areas, the use of retrograde viral vectors has allowed the access to morphological features of a distinct population of neurons selected by connectivity rules (i.e. all neurons projecting to pMO). Our analysis has revealed a stage-dependent loss of spines in all neurons projecting to pMO (first in layer V, then in layer II/III), which appears independent of misfSOD burden and involves neuronal populations in areas whose vulnerability was previously unrecognized, such as AUD neurons. Spine loss has been detected as a morphological counterpart of the early stages of the degenerative process ongoing in the layer V of the motor cortex of mSOD mice (Fogarty et al., 2015). The longitudinal study of multiple brain areas in Golgi-staining preparations (Fogarty et al., 2016) has revealed the early and progressive spine loss affecting the basal dendrites of pyramidal neurons from MO, medial pre-frontal (MPFC) cortex; notably, SS neurons appeared to be affected only very late in disease progression (P120; Fogarty et al., 2016). By identifying a neuronal subpopulation of SS based on connectivity rules (i.e. only the neurons projecting to pMO), we report an earlier loss of spines than previously though (P50 vs P120). This finding not only underscores the importance of assessing disease manifestations in homogeneous (and characterized by connectivity rules) populations but also suggests that spine loss may be a phenomenon occurring in a larger set of neurons, possibly across the whole motor subnetwork.
Our study displays some limitations. Firstly, in order to maintain tractability while still addressing the core question of the large-scale remodeling of the structure of projections in the motor subnetwork, we decided to limit the analysis to the forebrain, avoiding the analysis of brainstem nuclei and cerebellum. Although these regions provide relevant input to the motor subsystem, they are not usually identified in fMRI studies of ALS and therefore their contribution to the definition of the hyperconnectivity phenotype may be more limited. Second, although it is desirable to have large number of biological replicates, the need to maintain the genotype blinding and to manually identify and anatomically annotate each neuron (more than 150,000 for the present work) did not allow for a n > 3–6, thus underpowering the study for the comparison of structures with fewer number of projecting neurons. Although the present study is in line with the numerosity currently manageable in manual reconstructions (Hooks et al., 2013; Yamawaki et al., 2016; DeNardo et al., 2015; Tervo et al., 2016; Kohl et al., 2018), the advancement of machine learning and automated neuronal annotation (Fürth et al., 2018) may allow a more comprehensive quantitative analysis of input network to pMO.
Third, we did not consider ipsilateral sMO and pMO for cell counting in our analysis. Due to the proximity of pMO and sMO, we could not have conclusively excluded the possibility that infected neurons in the ipsilateral sMO were derived from retrograde infection or virus diffusion, introducing ambiguity to the connectivity map. Further refinement of the retrograde-tracing approach is ongoing to provide a whole brain, full resolution map of the connectivity of the motor subsystem in ALS.
Taken together, our circuit-tracing data in mice and imaging data in humans converge to show that a significant and progressive remodeling of cortical circuits takes places in ALS and that this phenotype may be amenable to investigation in animal models. The mechanistic drivers of this remodeling are currently unclear, although it is possible to speculate that homeostatic mechanisms may be activated by alterations in excitability of pMO neurons, including corticospinal motor neurons (Kim et al., 2017). In fact, pathological decrease or increase in the firing properties of vulnerable neuronal subpopulations appears to be a shared phenotype of multiple neurodegenerative conditions (Roselli and Caroni, 2015) and may be one of the driving forces of adaptive cortical remodeling in neurodegeneration (LANDSCAPE Consortium et al., 2015; Rosskopf et al., 2017; Rosskopf et al., 2018). Thus, the combination of the analysis of large-scale projection patterns (enabled by retrograde viral tracing) with chronic manipulation of excitability and firing (such as chemiogenetics; Saxena et al., 2013) may prove fundamental to link functional and structural cortical phenotypes.
Materials and methods
Animals
All animal experiments in this study were performed in agreement with the guidelines for the welfare of experimental animals issued by the Federal Government of Germany and by the local ethics committee (Ulm University) and authorized under licence no. 1312. The following strains of transgenic mice were obtained from Jackson Laboratories: B6SJL73 Tg(SOD1*G93A)1Gur/J (high-copy, henceforth mSOD), B6.Cg-Gt(ROSA)26Sortm6/(CAG74 ZsGreen)Hze/J (henceforth ZsGreen-ROSA26) and B6.Cg-Gt(ROSA)26Sortm6/(CAG75 TdTomato)Hze/J (henceforth tdTomato-ROSA26). For the generation of double-transgenic lines, heterozygous male mice from the mSOD line were crossed with females homozygous for the reporter cassette (either ZsGreen or tdTomato lines). Resulting progeny was 50% WT/ROSA26ZsGreen+ and 50% mSOD/ROSA26ZsGreen+. Under the repression of the STOP cassette in the ROSA26 locus, the fluorescent protein (either ZsGreen either tdTomato) is not expressed in basal conditions in the reporter lines; however, upon expression of the Cre recombinase (in our experimental design, via an AAV vector) leads to the excision of the STOP cassette (since it is flanked by the lox recombination sites) and results in the expression of the reporter (Gerfen et al., 2013).
Male mice were used for experimental purposes at post-natal day (P)20, P60 and P90 (sacrificed at P35, P75 and P105). Mice were housed at 3–4 animals per cage, with unlimited access to food and water, a light/dark cycle of 14/10 hr and humidity between 40% and 60%. All mice expressing mSOD were routinely tested for motor impairment and euthanized in case of overt motor disability. All male mice from a given litter were employed and the experimenter was kept blind to their genotype: injection, perfusion, brain section, imaging and neuronal counts and anatomical annotation; the genotype was revealed only when each mouse was assigned to a group for the statistical analysis.
All efforts were made to comply to the ‘3R’ guidelines. In order to maintain the data analysis tractable and in absence of previous data on quantitative analysis of input connectivity in ALS mice, we estimated the sample size based on the mean number and standard deviation of thalamic neurons at P90, anticipating a genotype effect of at least 30–40%, with α = 0.05% and 90% power, resulting in group size between 3 and 5.
Generation and production of AAVs
Request a detailed protocolAAV-vectors were generated using a three plasmid standard protocol (Jungmann et al., 2017) with a pscAAV-CMV-Cre vector genome plasmid, derived from pscAAV-CMVEnh/MLC0.26-Cre (Werfel et al., 2014), pDGdeltaVP providing helper sequences, and a distinct rep-cap-helper plasmid with either AAV2 or AAV9 wild-type capsid genomes (resulting in WT-AAV2-CMV-Cre or WT-AAV9-CMV-Cre, respectively), or AAV9 variants containing distinct heptapeptide motifs reported previously (Varadi et al., 2012), resulting in AAV9-RGDLRVS-CMV-Cre, AAV9-SLRSPPS-CMV-Cre and AAV9-NSSRFTP-CMV-Cre. Vectors were purified and concentrated as described before (Jungmann et al., 2017).
Intracerebral injections
Request a detailed protocolIntracerebral injection of AAV vectors was performed as previously reported (Karunakaran et al., 2016). Briefly, mice were administered buprenorphine (0.01 mg/Kg; Reckitt Beckiser Healthcare, Berkshire, UK) and meloxicam (1.0 mg/Kg; Böhringer Ingelheim, Biberach an der Riß, Germany) 30 min before being put into a stereotaxic frame (Bilaney Consultants GmbH D-40211 Düsseldorf, Germany) under continuous isoflurane anesthesia (4% in O2 at 800 ml/min). Skin scalp was incised to expose the underlying bone and a hand-held micro drill was used to drill the burr hole in the opportune location. Stereotactic coordinates were chosen for each target area based on The Mouse Brain in Stereotaxic Coordinates (as in Oh et al., 2014) and corresponded to primary motor cortex (x =+ 1.5, y =+ 1.0, z = −0.5), secondary motor cortex (x =+ 0.7, y =+ 1.0, z=-0–5), striatum (x =+ 1.0, y =+ 1.5, z = −3.5), visual cortex (x = −3.0, y =+ 2.0, z = −0.6) and primary somatosensory cortex(x =+ −1.0, y = −1.0, z = 0.5). 500 nl of viral suspension, mixed with an equal volume of 1% Fast Green solution (100 mg of Fast Green in 10 mL of H2O), was injected using a pulled-glass capillary connected to a Picospritzer microfluidic device. The Fast-Green visible dye, used to visually monitor the injection, is quickly washed away in PBS and does not interfere with the imaging procedures. Injection was performed with 10 msec pulses every 30 s over a span of 5 min. The capillary was kept in place for 10 more minutes before being withdrawn to prevent backflow of the virus. Scalp skin was then stitched with Prolene 7.0 surgical thread and the animals were transferred to single cages with facilitated access to water and food for recovery. Animals were monitored for eventual neurological impairment for the following 72 hr and were administered additional doses of buprenorphine as needed.
Each of the following AAV vectors was independently injected at the titre of 5–9*1013/ml: wild-type (WT) AAV9-CMV-Cre (WT-AAV9), WT-AAV2-CMV-Cre (WT-AAV2), AAV9-RGDLRVS-CMV-Cre (AAV9-RGD), AAV9-SLRSPPS-CMV-Cre (AAV9-SLR) and AAV9-NSSRFTP-CMV-Cre (AAV9-NSS).
For cholera toxin injections, Alexa-488-conjugated or Alexa-647-conjugated cholera toxin B (CTb-488, Invitrogen) was diluted according to the manufacturer’s instructions to 1 µg/µl by gentle mixing (no vortexing). 1.0 µl of CTb was injected using the same procedure as reported above.
Immunohistochemistry
Request a detailed protocolMice were terminally anesthetized with cloralium hydrate and transcardially perfused with 50 ml ice-cold PBS followed by 2.5–3 ml/g of 4% paraformaldehyde (PFA) in phosphate buffered saline (PBS, pH 7.4). Brain samples were quickly dissected and post-fixed in 4% PFA in PBS at 4°C for 18 hr, washed in PBS and cryo-protected in 30% sucrose in PBS for 36 hr. Samples were then snap-frozen in OCT and sectioned at −15°C in a cryotome (Leica CM1950) to 70µm-thick sections serially collected. Free-floating brain sections were incubated for 2 hr at 24°C in blocking buffer (5%-donkey serum, 3% BSA, 0.3% Triton X-100), followed by incubation with primary antibodies, such as: rabbit anti-RFP (1:1000, Rockland Immunochemicals, Limerick, PA) for 24 hr at 4°C; mouse anti-Tyrosine Hydroxylase (1:2000, Sigma) for 24 hr at 4°C; mouse anti NeuN (1:100, Millipore) for 24 hr at 4°C; mouse monoclonal anti-misfolded SOD1 (B8H10, 1:1000, MediMabs) for 48 hr at 4°C. Brain sections were thereafter washed 3 × 45 min in PBS, 0.1% Triton, incubated for 2 hr at 24°C with the opportune secondary antibody combination in blocking buffer (donkey anti-rabbit Alexa-568; donkey anti-mouse Alexa-488, 1:500; Life Technologies), together with DAPI (0.1 µg/ml, Life Technologies), and washed again 3 × 45 min in PBS, 0.1% Triton. Brain sections were mounted using ProLong Gold Antifade (ThermoFisher) mounting medium.
Image acquisition and analysis
Request a detailed protocolIn order to analyze the local infection rate of the virus in the dorsal striatum (DS), we identified the core of injection according to the brightness level. Four regions of interest (ROI, 300 × 300 µm2 each) were drawn at its cardinal points and 500 µm3 of tissue was analyzed in total. For each ROI, reporter +neurons and DAPI +nuclei were counted. Local infectivity for each AAV variant was calculated by averaging ROI values for every mouse (N = 3 for each viral variant) and the values are expressed as a percentage of AAV9 variants over WT- AAV9. For analysis of the retrograde infectivity in motor cortex (MO), substantia nigra (SNc), thalamus (TH) and lateral geniculate nucleus (LGN), brain slices corresponding to our target brain regions were cut (N = 3 for each viral variant). Reporter +soma were counted for every region and the values are expressed in terms of percentage over WT-AAV9.
For analysis of input connectivity to primary motor cortex (pMO, considering all time points analyzed: WT N = 11, SOD N = 12) and visual cortex (VIS, WT N = 3, SOD N = 3) via AAV9-SLR variant, all sections corresponding to each brain were scanned using the BZ-9000 fluorescence microscope (KEYENCE) equipped with a 4x objective, with exposure time set at 200 msec. All images were acquired in a 12bit format. Images were processed using the ImageJ software suite: images of individual brain sections were loaded, background was subtracted offline via the appropriate tool to unambiguously identify reporter +cells. Labeled cells were positively identified as neurons (in contrast to non-neuronal cells or to dense clusters of axons) by the detection of a round soma surrounded by basal dendrites and, in cortical neurons, by the presence of an apical dendritic shaft. The soma of fluorescently labeled neurons was manually identified and annotated for anatomical brain regions according to Allen Brain Atlas coordinates. Total numbers of neurons were normalized (Oh et al., 2014) by the volume of the injection site in pMO (AAV9- SLR: 1.3 ± 0.6 mm3, N = 23) and VIS (AAV9-SLR: 1.1 ± 0.5 mm3, N = 6). Cortical layers were discriminated according to positional criteria and using the tissue autofluorescence (which highlights nuclei in negative, since they are not strongly autofluorescent) and the region of interest (ROI, rectangular in shape, approximately 200 µm high) was placed approximately 100 µm (upper border) ventral to the cortical (pial) surface for layer II/III, whereas for layer V and VI the lower border of the ROI was placed 300 µm dorsal of, or in contact with, respectively, the boundary between gray and white matter. The correctness of this positions was checked in a subset of NeuN and DAPI-stained sections.
For analysis projecting neurons we injected a small volume (500 nl out of a suspension of a titre of 5–9*1013/ml of virus suspension) of AAV9-SLR virus suspension in the pMO corresponding to the whisker area (according to Zingg et al., 2014) or in the neighboring secondary motor cortex of tdTomato-ROSA26 or ZsGreen-ROSA26 reporter mice (age P20). After sacrificing the animals 15 days post-injection, brains were serially sectioned to 70 µm sections, and each section was digitized for manual annotation. An average of 5000 + neurons were manually identified (N = 6). In order to make the analysis tractable, cerebellum and brainstem nuclei were not included in the present study. Moreover, although ipsilateral sMO is a well-known input to pMO, its analysis was not performed because of technical issues: its proximity to the injection site made it difficult to distinguish between retrograde infection and spreading of the virus.
For analysis of input connectivity to pMO via cholera toxin (CTb Alexa Fluor 488 conjugate, Thermofisher, WT N = 3, SOD N = 3), brain sections were imaged with a Zeiss 710, 20x objective (NA 0.8), 0.9 optical zoom at 12 bit depth. CTb + soma were manually counted in contralateral pMO (cMOp) and somatosensory cortex (SS) following anatomical coordinates according to the Allen Brain Atlas and neuronal counting was normalized to the injection site volume (1.2 ± 07 mm3, N = 6).
For analysis of neuronal cells number in pMO, 70 µm brain slices were cut for WT (N = 2) and SOD (N = 2) mice at the age of P90. Slices were stained for NeuN (Millipore). Images were acquired with AF6000 fluorescence microscope (Leica) equipped with a 10x objective, with exposure time set at 200 msec. All images were acquired in a 12-bit format. Single optical section images of pMO were analyzed for NeuN +soma via Imaris 7.6.5 software (Bitplane AG, ‘Spot’ tool, followed by manual correction).
For spine density analysis, fluorescent images were acquired with a Zeiss 710 confocal microscope, with a 63x oil-immersion objective (NA 1.40) and 0.9 optical zoom at 12bit depth. Cortical neurons expressing the reporter protein, therefore part of the motor network, were acquired in their full dimension, that is dendrites did not exit the plane of the coronal brain slice before reaching the dendritic terminus. Criteria for selection of neurons were anatomical and morphological: selected neurons were located in the desired cortical layer, showing the distinctive pyramidal shape and displaying at least three uncorrupted basal dendrites. For each neuron, three to six, 30–100 μm long, intermediate dendritic segments (10 µm distant from the tip) were acquired for spine counts. After image collection, each segment was analyzed for spine density using the Imaris 7.6.5 software (Bitplane AG, ‘filament tracer’ tool followed by manual correction). Only protrusions with a clear connection of the head of the spine to the dendrite shaft were counted as spines; small processes were classified as a spine only if they were <3 µm long and <0.8 µm in cross- sectional diameter (Harris, 1999). Spine density was then expressed as the number of spines per 10 μm dendrite length (Saba et al., 2016; Jara et al., 2012; Vinsant et al., 2013). Number of mice N = 3 for both WT and SOD for each time point.
To analyze misfSOD1 expression distribution within the motor network, we performed pMO injection of mSOD animals at P20 (sacrificed at P35, N = 3) with CTb. Brain slices were stained with the mouse monoclonal anti misfolded SOD1 B8H10 antibody (MediMabs, Montreal, Canada). Imaging was performed with a Zeiss 710, 20x objective (NA 0.8), 0.9 optical zoom at 12bit depth and images were analyzed using ImageJ. For misfSOD1 +neurons assessment, an ROI of 165 mm2 was drawn for each cortical layer, a threshold was set and only neurons above the threshold were considered positive. For counting CTb + neurons in layer II/III, we selected a total range of 200 µm, 100 µm away from the Pia and for layer V, we considered a total range of 200 µm, 400 µm away from the grey matter. The perimeter of misfSOD +nuclei was manually drawn for both CTb + and CTb- in order to calculate and compare their mean grey value.
Statistical analysis - mouse experiments
Request a detailed protocolStatistical analysis was performed with the GraphPad Prism six software suite. For comparing normalized neuronal counts between viral variants, ordinary one-way ANOVA with Dunnett multiple comparison test was used. The same test was applied for analysis of input to V1 using AAV-SLR and to pMO using CTb. For pMO and sMO input connectivity via AAV9-SLR, including comparison between the different cortical layers, either raw and normalized neuronal counts was compared by two-way ANOVA with Sidak correction for multiple comparison. Data regarding NeuN +neuronal count in primary motor cortex and misfSOD expression in CTb-labeled neurons have been analyzed via unpaired t-Test with Mann-Whitney correction. For the characterization of the variability of the dataset and for the description of the distribution of the data points, coefficient of variation (standard deviation/mean) and Skewness coefficient have been calculated (and are provided in tabulated format along with the mean, median and standard deviation).
Skewness value provides a descriptor of the slope of data distribution. In a perfect normal distribution curve, data distribute symmetrically on both sides of the peak, and the value of the skewness is zero. In a non-gaussian distribution, skewness value is positive when data piles up on the peak’s left side and negative when data piles up on the peak’s right side (and the tail points left). Skewness value between −2 and +2 are accepted as indication of non-biased data distribution (George and Mallery, 2010).
The magnitude of the difference between WT and mSOD was described using the parameter known as the effect size ((mean mSOD-mean WT)/SD of either group; calculated according to Sullivan and Feinn, 2012). Unlike significance tests, effect size is independent of sample size (Ferguson, 2009).
Data are displayed as mean ± SD and values for median, skewness and effect size are provided in tables as supplementary file. Statistical significance was set at p<0.05 before multiple-comparisons correction: statistical significance is indicated by p<0.05 and the magnitude of the effect is highlighted by the effect size. Percentage changes are reported in relation to relative WT.
Functional connectivity data validation in humans
Ethical approval of human study
Request a detailed protocolAll subjects included in the human study provided written informed consent according to institutional guidelines; the consent includes the declaration of the understanding of the study design, the agreement to the participation to the study, to the publication of the results, and to the data protection and anonymization procedures (under the chaptes Einwilligungserklaerung’, ‘Probandeninformation’, ‘Darstellung der Experimente’, ‘Datenschutzerklärung’). The study was approved by the Ethics Committee of Ulm University, Ulm, Germany (reference #19/12) and was performed in accordance with the ethical standards laid down in the 1964 Declaration of Helsinki and its later amendments.
Data acquisition in humans
Request a detailed protocol‘Resting-state’ functional magnetic resonance imaging (rs-fMRI) data were acquired from 71 symptomatic patients with ALS, diagnosed according to the recently revised El Escorial diagnostic criteria (for The WFN Research Group On ALS/MND et al., 2015), and 28 matched healthy controls (see Table 1 for clinical and demographic characteristics). Human whole-brain based echo planar images (EPI) were acquired for all subjects at a 3T Siemens Allegra MR scanner (Siemens Medical, Erlangen, Germany) using a blood oxygen level dependent (BOLD) sensitized rs-fMRI sequence (30 transversal slices, gap 1 mm, 3 × 3×4 mm3 voxels, 64 × 64×30 matrix, FOV 192 × 192×149 mm3, TE 30 ms, TR2000ms, flip angle 90°, 300 volumes).
Preprocessing of functional data
Request a detailed protocolThe human functional data analysis followed a standardized procedure (Gorges et al., 2014) and was performed using the Tensor Imaging and Fiber Tracking (TIFT) software package (Müller et al., 2007). Preprocessing included (1) spatial upsampling to an 1 × 1×1 mm3 isogrid (matrix, 256 × 256×256), (2) motion correction using a rigid brain transformation (six degrees of freedom), (3) stereotaxic Montreal Neurological Institute (MNI) normalization (Brett et al., 2002) using a landmark-based deformation approach (Müller and Kassubek, 2013), (4) temporal demeaning and linear detrending, (5) temporal bandpass filter (0.01 < f < 0.08 Hz), and (6) spatial smoothing using a 7 mm 3-dimensional full-width at half maximum (FWHM) Gaussian blur filter (LANDSCAPE Consortium et al., 2016; LANDSCAPE Consortium et al., 2015).
Functional data analysis
Request a detailed protocolTen spherical a-priori-defined bilateral seed regions (r = 1 mm) were chosen based on the anatomical structures in order to compute functional connectivity maps of the motor system and the visual association system (reference) as follows: (1) primary motor system (seed voxels; MNI coordinates (x, y, z):±15, 30, 69), primary somatosensory cortex (±15,–30, 69), dorsolateral prefrontal association cortex (±22, 7, 55), Thalamus (±12,–29, 0), and visual association cortex (±51,–70, −13).
Functional connectome reconstruction
Request a detailed protocolWe computed a functional connectome of the motor and visual association system using the defined seed regions by computing pairwise region-to-region functional connectivity and by computing Pearson’s product moment correlation coefficient (r-value) pairwise between voxel time series. Correlation coefficient was Fisher's r- to z-transformed to normally distributed z(r) score for further statistical analysis.
Interference statistics
Request a detailed protocolOverall functional connectivity was computed and used as a regressor (van den Heuvel et al., 2017) prior to subjecting the data to a two-sided parametric Student's t-test in order to test for pairwise differences between ALS patients and controls using the network-based statistics approach (Zalesky et al., 2010).
Data availability
All the murine data generated or analysed during this study are included in the manuscript and supporting files. The raw images are deposited on the Dataverse database (https://doi.org/10.7910/DVN/5VNSXE)
-
Replication Data for: Stage-dependent remodeling of large-scale architecture of projections to primary motor cortex in ALS mouse model revealed by a new variant retrograde-AAV9Publicly available at the Harvard dataverse (DOI: 10.7910/DVN/5VNSXE).
References
-
Divergent brain network connectivity in amyotrophic lateral sclerosisNeurobiology of Aging 34:419–427.https://doi.org/10.1016/j.neurobiolaging.2012.04.015
-
Therapeutic targeting of the pathological triad of extrasynaptic NMDA receptor signaling in neurodegenerationsThe Journal of Experimental Medicine 7:jem.20161673.https://doi.org/10.1084/jem.20161673
-
Wild-type and mutant SOD1 share an aberrant conformation and a common pathogenic pathway in ALSNature Neuroscience 13:1396–1403.https://doi.org/10.1038/nn.2660
-
Amyotrophic lateral sclerosis--a model of corticofugal axonal spreadNature Reviews Neurology 9:708–714.https://doi.org/10.1038/nrneurol.2013.221
-
The problem of functional localization in the human brainNature Reviews Neuroscience 3:243–249.https://doi.org/10.1038/nrn756
-
Stages of pTDP-43 pathology in amyotrophic lateral sclerosisAnnals of Neurology 74:20–38.https://doi.org/10.1002/ana.23937
-
Involvement of spinal sensory pathway in ALS and specificity of cord atrophy to lower motor neuron degenerationAmyotrophic Lateral Sclerosis and Frontotemporal Degeneration 14:30–38.https://doi.org/10.3109/17482968.2012.701308
-
The "dying-back" phenomenon of motor neurons in ALSJournal of Molecular Neuroscience 43:470–477.https://doi.org/10.1007/s12031-010-9467-1
-
Connectivity of mouse somatosensory and prefrontal cortex examined with trans-synaptic tracingNature Neuroscience 18:1687–1697.https://doi.org/10.1038/nn.4131
-
An effect size primer: a guide for clinicians and researchersProfessional Psychology: Research and Practice 40:532–538.https://doi.org/10.1037/a0015808
-
An interactive framework for whole-brain maps at cellular resolutionNature Neuroscience 21:139–149.https://doi.org/10.1038/s41593-017-0027-7
-
BookSPSS for Windows Step by Step: A Simple Guide and Reference, 17.0 UpdateNeedham Heights, MA: Allyn & Bacon, Inc.
-
Revealing the secrets of neuronal circuits with recombinant Rabies virus technologyFrontiers in Neural Circuits 7:2.https://doi.org/10.3389/fncir.2013.00002
-
Additional resources and the default mode network: evidence of increased connectivity and decreased white matter integrity in amyotrophic lateral sclerosisAmyotrophic Lateral Sclerosis and Frontotemporal Degeneration 15:537–545.https://doi.org/10.3109/21678421.2014.911914
-
Network structure of cerebral cortex shapes functional connectivity on multiple time scalesProceedings of the National Academy of Sciences 104:10240–10245.https://doi.org/10.1073/pnas.0701519104
-
Organization of cortical and thalamic input to pyramidal neurons in mouse motor cortexJournal of Neuroscience 33:748–760.https://doi.org/10.1523/JNEUROSCI.4338-12.2013
-
Protocol for efficient generation and characterization of Adeno-Associated viral vectorsHuman Gene Therapy Methods 28:235–246.https://doi.org/10.1089/hgtb.2017.192
-
Imaging the pathoanatomy of amyotrophic lateral sclerosis in vivo: targeting a propagation-based biological markerJournal of Neurology, Neurosurgery & Psychiatry 89:374–381.https://doi.org/10.1136/jnnp-2017-316365
-
MRI-based functional neuroimaging in ALS: an updateAmyotrophic Lateral Sclerosis 10:258–268.https://doi.org/10.3109/17482960802353504
-
Increased functional connectivity common to symptomatic amyotrophic lateral sclerosis and those at genetic riskJournal of Neurology, Neurosurgery & Psychiatry 87:580–588.https://doi.org/10.1136/jnnp-2015-311945
-
Changes of resting state brain networks in amyotrophic lateral sclerosisExperimental Neurology 217:147–153.https://doi.org/10.1016/j.expneurol.2009.01.025
-
The innate immune response to adenovirus vectorsHuman Gene Therapy 15:1157–1166.https://doi.org/10.1089/hum.2004.15.1157
-
Diffusion tensor magnetic resonance imaging in the analysis of neurodegenerative diseasesJournal of Visualized Experiments 28:.https://doi.org/10.3791/50427
-
A large-scale multicentre cerebral diffusion tensor imaging study in amyotrophic lateral sclerosisJournal of Neurology, Neurosurgery & Psychiatry 87:570–579.https://doi.org/10.1136/jnnp-2015-311952
-
Mitochondrial damage revealed by immunoselection for ALS-linked misfolded SOD1Human Molecular Genetics 22:3947–3959.https://doi.org/10.1093/hmg/ddt249
-
Hyperconnective and hypoconnective cortical and subcortical functional networks in multiple system atrophyParkinsonism & Related Disorders 49:75–80.https://doi.org/10.1016/j.parkreldis.2018.01.012
-
Gene therapy for the CNS using AAVs: the impact of systemic delivery by AAV9Journal of Controlled Release 241:94–109.https://doi.org/10.1016/j.jconrel.2016.09.011
-
Information processing architecture of functionally defined clusters in the macaque cortexJournal of Neuroscience 32:17465–17476.https://doi.org/10.1523/JNEUROSCI.2709-12.2012
-
Network structure shapes spontaneous functional connectivity dynamicsJournal of Neuroscience 35:5579–5588.https://doi.org/10.1523/JNEUROSCI.4903-14.2015
-
Terminal N-linked galactose is the primary receptor for adeno-associated virus 9The Journal of Biological Chemistry 286:13532–13540.https://doi.org/10.1074/jbc.M110.210922
-
Using effect Size-or why the P value is not enoughJournal of Graduate Medical Education 4:279–282.https://doi.org/10.4300/JGME-D-12-00156.1
-
Does interneuronal dysfunction contribute to neurodegeneration in amyotrophic lateral sclerosis?Amyotrophic Lateral Sclerosis 13:245–250.https://doi.org/10.3109/17482968.2011.636050
-
A revision of the el escorial criteria - 2015Amyotrophic Lateral Sclerosis and Frontotemporal Degeneration 16:291–292.https://doi.org/10.3109/21678421.2015.1049183
-
A corticocortical circuit directly links retrosplenial cortex to M2 in the mouseThe Journal of Neuroscience 36:9365–9374.https://doi.org/10.1523/JNEUROSCI.1099-16.2016
Article and author information
Author details
Funding
Deutsche Forschungsgemeinschaft (GRK1789)
- David Bayer
Deutsche Forschungsgemeinschaft (SFB1149)
- Francesco Roselli
- Tobias M Boeckers
Ulm University Medical Faculty (Baustein program)
- Francesco Roselli
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
This work has been supported by the Baustein program of Ulm University Medical Faculty, by the Deutsche Forschungs Gemeinschaft (DFG) as part of the Graduate School in Cellular and Molecular Mechanisms in Aging and by the Synapsis Foundation. FR is also supported by the DFG as part of the Collaborative Research Center 1149 ‘Danger Response, Disturbance Factors and Regenerative Potential after Acute Trauma’. This study was also supported by the German Research Foundation (Deutsche Forschungsgemeinschaft, DFG Grant Number LU 336/15–1), and the German Network for Motor Neuron Diseases (BMBF 01GM1103A). The authors thank the Ulm University Center for Translational Imaging MoMAN for its support. BC and DB are members of the International Graduate School in Molecular Medicine at Ulm University; DB is part of the Graduate School in Cellular and Molecular Mechanisms in Aging at Ulm University. The authors wish to acknowledge Mrs. Sonja Fuchs for MRI data acquisition of the human cohort; Ms. Silvia Cursano and Ms. Najwa Ouali Alami for the support and collaboration; Dr. Kim Moorwood and Mrs. Akila Chandrasekar for proofreading.
Ethics
Human subjects: All subjects included in the human study provided written informed consent according to institutional guidelines; the consent includes the declaration of the understanding of the study design, the agreement to the participation to the study, to the publication of the results, and to the data protection and anonymization procedures (under the chaptes Einwilligungserklaerung", "Probandeninformation", "Darstellung der Experimente", "Datenschutzerklärung"). The study was approved by the Ethics Committee of Ulm University, Ulm, Germany (reference #19/12) and was performed in accordance with the ethical standards laid down in the 1964 Declaration of Helsinki and its later amendments.
Animal experimentation: All animal experiments have been approved by the Regierungpraesidium Tubingen under the licence no. 1312. All animals were handled according to the federal regulations on animal experimentations and under the supervision of the local veterinary office. Every effort was made to adhere to the 3R guidelines and to minimize suffering.
Copyright
© 2018, Commisso et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 2,650
- views
-
- 440
- downloads
-
- 22
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Neuroscience
Neuromodulatory inputs to the hippocampus play pivotal roles in modulating synaptic plasticity, shaping neuronal activity, and influencing learning and memory. Recently, it has been shown that the main sources of catecholamines to the hippocampus, ventral tegmental area (VTA) and locus coeruleus (LC), may have overlapping release of neurotransmitters and effects on the hippocampus. Therefore, to dissect the impacts of both VTA and LC circuits on hippocampal function, a thorough examination of how these pathways might differentially operate during behavior and learning is necessary. We therefore utilized two-photon microscopy to functionally image the activity of VTA and LC axons within the CA1 region of the dorsal hippocampus in head-fixed male mice navigating linear paths within virtual reality (VR) environments. We found that within familiar environments some VTA axons and the vast majority of LC axons showed a correlation with the animals’ running speed. However, as mice approached previously learned rewarded locations, a large majority of VTA axons exhibited a gradual ramping-up of activity, peaking at the reward location. In contrast, LC axons displayed a pre-movement signal predictive of the animal’s transition from immobility to movement. Interestingly, a marked divergence emerged following a switch from the familiar to novel VR environments. Many LC axons showed large increases in activity that remained elevated for over a minute, while the previously observed VTA axon ramping-to-reward dynamics disappeared during the same period. In conclusion, these findings highlight distinct roles of VTA and LC catecholaminergic inputs in the dorsal CA1 hippocampal region. These inputs encode unique information, with reward information in VTA inputs and novelty and kinematic information in LC inputs, likely contributing to differential modulation of hippocampal activity during behavior and learning.
-
- Computational and Systems Biology
- Neuroscience
Animal behaviour alternates between stochastic exploration and goal-directed actions, which are generated by the underlying neural dynamics. Previously, we demonstrated that the compositional Restricted Boltzmann Machine (cRBM) can decompose whole-brain activity of larval zebrafish data at the neural level into a small number (∼100-200) of assemblies that can account for the stochasticity of the neural activity (van der Plas et al., eLife, 2023). Here, we advance this representation by extending to a combined stochastic-dynamical representation to account for both aspects using the recurrent temporal RBM (RTRBM) and transfer-learning based on the cRBM estimate. We demonstrate that the functional advantage of the RTRBM is captured in the temporal weights on the hidden units, representing neural assemblies, for both simulated and experimental data. Our results show that the temporal expansion outperforms the stochastic-only cRBM in terms of generalization error and achieves a more accurate representation of the moments in time. Lastly, we demonstrate that we can identify the original time-scale of assembly dynamics by estimating multiple RTRBMs at different temporal resolutions. Together, we propose that RTRBMs are a valuable tool for capturing the combined stochastic and time-predictive dynamics of large-scale data sets.