Shared and specific signatures of locomotor ataxia in mutant mice
Abstract
Several spontaneous mouse mutants with deficits in motor coordination and associated cerebellar neuropathology have been described. Intriguingly, both visible gait alterations and neuroanatomical abnormalities throughout the brain differ across mutants. We previously used the LocoMouse system to quantify specific deficits in locomotor coordination in mildly ataxic Purkinje cell degeneration mice (pcd; Machado et al., 2015). Here, we analyze the locomotor behavior of severely ataxic reeler mutants and compare and contrast it with that of pcd. Despite clearly visible gait differences, direct comparison of locomotor kinematics and linear discriminant analysis reveal a surprisingly similar pattern of impairments in multijoint, interlimb, and whole-body coordination in the two mutants. These findings capture both shared and specific signatures of gait ataxia and provide a quantitative foundation for mapping specific locomotor impairments onto distinct neuropathologies in mice.
Introduction
Visibly ataxic mouse mutants exhibit varying patterns of neuropathology throughout the brain (Cendelin, 2014; Fortier et al., 1987; Goldowitz et al., 1997; Lalonde and Strazielle, 2007; Lalonde and Strazielle, 2019; Mullen et al., 1976; Walter et al., 2006). Although their motor coordination deficits are generally attributed to abnormal cell patterning within the cerebellum (Arshavsky et al., 1983; Orlovsky et al., 1999), these lines have distinct patterning defects within the cerebellum, varying degrees of extracerebellar involvement, and differences in age of onset (Cendelin, 2014; Lalonde and Strazielle, 2019). The nature of the motor deficits exhibited by these mice also varies, and can often be distinguished by trained observers (Berman, 2018; Brooks and Dunnett, 2009; Hoogland et al., 2015; Lalonde and Strazielle, 2019; Schiffmann et al., 1999; Stroobants et al., 2013; Vinueza Veloz et al., 2015). However, analysis of motor coordination is often limited to low dimensional descriptions of limited specificity that fail to distinguish between related behavioral phenotypes (Brooks and Dunnett, 2009; Lalonde and Strazielle, 2019). Analysis of locomotor kinematics can provide higher dimensional readouts of locomotor behavior (Cendelín et al., 2010; Gabriel et al., 2009; Zörner et al., 2010), but can still suffer from a lack of specificity due to an abundance of highly correlated measures that ultimately reflect non-specific features such as changes in walking speed or body size (Batka et al., 2014; Cendelín et al., 2010; Machado et al., 2015). A quantitative understanding of the specific nature of gait ataxia in mutants with well-described abnormalities in circuit architecture could provide important clues into neural mechanisms of motor coordination (Anderson and Perona, 2014; Bastian et al., 1996; Berman, 2018; Brown and de Bivort, 2018; Darmohray et al., 2019; Datta et al., 2019; Kiehn, 2016; Morton and Bastian, 2007; Powell et al., 2015; Sarnaik and Raman, 2018; Udo et al., 1980).
We previously used the LocoMouse system (Machado et al., 2015) to analyze the locomotor coordination of mildly ataxic Purkinje cell degeneration (pcd) mice, in which neural degeneration, particularly early in postnatal development, is largely restricted to cerebellar Purkinje cells, effectively disconnecting the output of the cerebellar cortex (Chen et al., 1996; Fernandez-Gonzalez et al., 2002; Le Marec and Lalonde, 1997). We found that locomotor deficits in pcd were restricted to specific aspects of multijoint, interlimb, and whole-body coordination, while the forward trajectories of individual paws were spared (Machado et al., 2015). We further found that the tail movements of pcd mice reflected the passive consequences of limb movement (Machado et al., 2015). However, it remained unclear to what extent these features represented fundamental features of cerebellar ataxia, or were specific to pcd mice.
Reeler mice are a classic ataxic mutant (Cendelin, 2014; Curran and D'Arcangelo, 1998; D'Arcangelo et al., 1999; D'Arcangelo et al., 1995; Falconer, 1951) with an autosomal recessive mutation in the reelin gene, which is important for neural cell migration (Beckers et al., 1994; Hack et al., 2002). Its loss causes several defects, in particular aberrant localization of neurons and failure of neuronal layer formation. Several brain regions are affected, including cerebellum (Hamburgh, 1963; Terashima et al., 1983), hippocampus (Stanfield et al., 1979), neocortex (Mikoshiba et al., 1980), inferior olive (Blatt and Eisenman, 1988) and substantia nigra (Kang et al., 2010). Neuropathology in these mice is particularly striking within the cerebellum, where severe irregularities in cellular localization are also associated with corresponding aberrant synaptic connectivity between cell types, abnormal foliation, and hypoplasia. Although their locomotor kinematics and whole-body coordination have not been reported, homozygous reeler mutants have been described as having a severely ataxic, ‘reeling’ gait, with difficulties in maintaining their hindquarters upright (Cendelin, 2014; Lalonde et al., 2004; Lalonde and Strazielle, 2019). Like most ataxic mutants, reelers also exhibit poor performance in rotarod, stationary beam and water maze tests (Lalonde et al., 2004).
Thus, pcd and reeler mice share grossly abnormal cerebellar circuitry, but exhibit marked differences in synaptic connectivity within the cerebellum and across the brain. We wondered whether these similarities and differences on the anatomical level might be associated with similarly shared and distinct features of motor behavior. Here we analyze the locomotor behavior of reeler mutants and compare it quantitatively to that of the more mildly ataxic pcd mice (Machado et al., 2015). Detailed comparison of locomotor kinematics and linear discriminant analysis reveals both shared and distinct features of gait ataxia in these two mutants. This approach provides a quantitative foundation for mapping specific locomotor impairments onto distinct neuropathologies.
Results
Reeler mice have impaired hindlimb control and exhibit increased variability of movement
Reeler mice exhibited visible and severe gait ataxia when walking on the LocoMouse setup (Video 1). Like pcd mice (Machado et al., 2015), reelers were smaller and walked more slowly than control littermates (Materials and methods; Figure 1D). However, the locomotor phenotypes of reeler and pcd mice were clearly distinguishable by eye, with reeler mice appearing much more severely ataxic than the mildly ataxic pcd mice (Video 1; Lalonde and Strazielle, 2007; Lalonde and Strazielle, 2019; Machado et al., 2015).
We analyzed the locomotor phenotype of reeler mice using the quantitative framework for locomotor coordination that we established previously (Figure 1A–C; Machado et al., 2015). First, the equations we previously generated with mixed-effects linear models (Machado et al., 2015) to predict stride parameters based on walking speed and body size were able to accurately predict the individual stride parameters (including stride length, cadence, stance duration) of reeler mice (Figure 1. E-G), as we had previously shown for pcd (Machado et al., 2015). Moreover, the continuous forward trajectories of reeler front paws were similar to those of size-matched controls, across walking speeds (Figure 1H, top). Notably, the forward trajectories of reeler hind paws exhibited lower forward velocities compared to size and speed-matched controls (Figure 1H, bottom), and they had accompanying increases in swing duration (from size-matched controls = 98.78 ± 6.48 ms to reeler = 114.72 ± 8.21 ms). Comparison of 3D trajectories revealed clear differences between reeler and control mice in both the side-to-side and the vertical paw movements of both front and hind paws (Figure 1I–K). Finally, and perhaps surprisingly given their severe ataxia, reeler mice did not exhibit an increased width of base of support (Figure 1K).
Direct comparison of average paw kinematics in the two mutants reveals remarkable similarities (Figure 2). Both reeler and pcd exhibited altered 3D trajectories of all paws when compared to those of speed and size-matched controls (Figure 2E–G). In particular, the off-axis (side-to-side and vertical) movements of all paws showed nearly identical alterations in the two mutants (Figure 2F,G). In addition, forward hind paw trajectories were profoundly affected in reelers, while more subtle effects were seen on the forward movement of the front paws (Figure 2B–E).
In contrast to the broad similarities in averaged paw trajectories, there were clear differences in the variability of paw kinematics between reeler and pcd (Figure 2H,I). Despite their ataxia, neither the front nor hind limb trajectories of pcd mice exhibited increased variability (Figure 2H; Machado et al., 2015). Paw movement was generally more variable in reeler, including the forward motion of the hind paws (Figure 2H-bottom) and the vertical movements of both front and hind paws (Figure 2I-top and bottom).
Impaired interlimb coordination and increased front paw support in Reeler
Mice typically walk in a symmetrical trot pattern across a wide range of walking speeds (Figure 3A; Machado et al., 2015). The normal pattern of interlimb coordination was markedly disrupted in reeler, due to specific and consistent changes in the phase relationship between front and hind limbs (Figure 3B–C). Remarkably, the alterations in front-hind limb stance phasing in reelers were identical on average to those of pcd mice; in both mutants, hind paw touch downs were delayed relative to their diagonal partners (Figure 3B,C; Machado et al., 2015). In contrast, relative left-right stance phasing of both the front and hind limbs was intact in both reelers and pcd (Figure 3D; Machado et al., 2015). Consistent with the increased variability of hind limb movements, the front-hind limb phasing was more variable on average in reeler, but not pcd (size-matched controls = 0.25 ± 0.05%, pcd = 0.29 ± 0.09% F1,82 = 1.02, p=0.32, reeler = 0.48 ± 0.11% F1,82 = 5.79, **p=9.5×10−3), and left-right phasing was not more variable in either mutant (size-matched controls = 0.25 ± 0.06%, pcd = 0.19 ± 0.07% F1,54 = 0.23, p=0.63, reeler = 0.22 ± 0.07%, F1,74 = 0.11, p=0.75).
We also observed changes in support patterns (ie, the configuration of paws that are in stance at any given time, Gorska et al., 1998) in reeler mice. At most natural walking speeds, wildtype mice typically have a single diagonal pair of paws on the ground at any given time (Machado et al., 2015). Reeler mice exhibited an increase in 3-paw support patterns (Figure 3E). They also spent more time in unstable support configurations such as non-diagonal 2-paw support (Figure 3F) and 2-front paw supports (Figure 3G). This increased instability was also observed in pcd (Machado et al., 2015) and is consistent with impaired interlimb coordination rather than a simple switch to a different gait pattern (Bellardita and Kiehn, 2015).
Although the alterations in interlimb coordination were similar in reeler and pcd, there were some notable differences between the two mutants. In particular, reeler mice spent more time with both front paws on the ground than both controls and pcd, either as sole supports (Figure 3G), or as part of a 3 or four paw support configuration (Figure 3H). Interestingly, reelers also spent less time with both hind paws on the ground than pcd or control mice (Figure 3I). The ratio of front to hind paw double support was higher in reeler (front/hind = 0.43) than in pcd (front/hind = 0.14) and control mice (front/hind = 0.10). This decreased hind paw double support is likely to be a consequence of the impaired hindlimb control described in Figure 1 and Figure 2.
Impaired whole-body coordination in reeler reflects passive consequences of front paw motion
Like pcd (Machado et al., 2015), reeler mice exhibited large side-to-side oscillations of the tail and nose with respect to the stride cycle (Figure 4A–F, green) when compared to controls (black). Also like pcd, both tail (Figure 4A,C) and nose (Figure 4D,F) movements became increasingly phase-lagged relative to the locomotor cycle at faster speeds.
We previously showed that the tail and nose movements of pcd mice could be successfully modelled as a passive consequence of hind limb movement (Machado et al., 2015). A model that converted hindlimb oscillation into nose and tail trajectories using solely the geometric relationships between body parts moving with fixed time delays accurately predicted the side-to-side tail and nose movements of pcd mice across walking speeds.
At first glance, the similarities in averaged tail trajectories suggest that the reeler tail and nose movements, like pcd, might also reflect passive consequences of limb movement during locomotion. However, the specific phase relationships of both the tail and nose with respect to the locomotor cycle were dramatically different in reeler compared to pcd (Figure 4B,E). We wondered whether these phase differences could be accounted for by the shift in support patterns towards the front limbs in reelers (Figure 3F–H). To our surprise, simply aligning the reeler tail oscillations to the front limbs, rather than hind limbs (Figure 4C), immediately revealed that the phase relationship between the reeler tail and the front limbs was nearly identical to that of the pcd tail to the hind limbs (Figure 4C). In other words, the tail movements of reeler mice have the same quantitative relationship to front limb movement that the pcd tail had to the hind limbs.
Nose trajectories in reeler were similarly pronounced, but also phase-shifted, compared to those of pcd (Figure 4D,E). Side-to-side nose movements in reeler were almost perfectly out of phase with the base of the tail (Figure 4F), suggesting that both the tail and nose movements of reeler mice oscillate along a single, straight body axis with each stride (Figure 4G). One possible explanation for the differences in relative tail and nose phasing between pcd and reeler (Figure 4E) is that the shift of supports to front paws leads to the loss of a front-limb steering component in reeler (Figure 4G). Consistent with this idea, reelers, but not pcd, exhibited larger side-to-side movements of the body center while walking (Figure 4H, Video 1).
To test the idea that the nose and tail movements of reeler mice might reflect the passive consequence of front, rather than hind, limb movement, we built an analytical model that computed predicted lateral trajectories of the tail and nose directly from the forward-backward oscillations of the front limbs (Figure 4I; Materials and methods). The model and its parameters were analogous to the geometrical model described in Machado et al., 2015, but with a shift to the front limbs and a single body axis for reeler (see Materials and methods). The model output accurately reproduced the phases as well as the actual trajectories of the tail and nose of reeler mice, across a range of walking speeds (Figure 4J and Figure 4—figure supplement 1A–C). Notably, hindlimb-driven models, although optimal for pcd (Figure 4—figure supplement 1H,I), performed less well at predicting reeler tail and nose movements, even allowing for the possibility of substantially longer time delays (Figure 4—figure supplement 1D-F). Similarly, unlike reeler, a single body-axis model was not appropriate for the pcd data.
Thus, in reeler, like pcd, the tail and nose appear to move as a passive consequence of forward limb motion. The altered phasing in the two mutants suggests that the shift in the support patterns towards the front paws in reeler (Figure 3, and presumably resulting from impaired hindlimb control shown in Figure 1) causes the tail and body axis to oscillate as a passive consequence of the front, rather than hind limbs.
Linear discriminant analysis reveals shared and specific features of gait ataxia
The results so far have described a comprehensive set of locomotor features and highlighted similarities and differences between reeler and pcd mice. Finally, we sought an unbiased way to summarize and conceptualize these findings (Brown and de Bivort, 2018; Berman, 2018; Datta et al., 2019). To do this we turned to linear discriminant analysis (LDA; Fisher, 1936; James et al., 2013), as a means to quantitatively distinguish the three groups of mice (wildtype, pcd and reeler; Figure 5).
Forty-five variables representing the various features of paw, nose, and tail movements that we measured during locomotion were extracted for this analysis (Figure 5—figure supplement 1). Because this relatively high-dimensional dataset contained many highly correlated variables (Figure 5—figure supplement 2A), we first applied principal component analysis (PCA) to account for inter-variable correlations and avoid overfitting (see Materials and methods, Figure 5—figure supplement 2B–F).
LDA successfully captured meaningful differences in an abstract feature subspace to separate the three genotypes based on their locomotor phenotypes across walking speeds (Figure 5). Unexpectedly, this analysis revealed that two distinct, orthogonal axes effectively separate controls from mutants (LD1), and the two mutants from each other (LD2) (Figure 5). The locomotor features contributing most strongly to these two linear discriminants thus correspond to the shared (LD1) and distinct (LD2) features of ataxia in reeler and pcd mice.
Inspection of the contributions of each gait parameter to the two LDs reveals that LD1, which separates controls from ataxic mutants, is highly influenced by variables representing 3D paw trajectories and interlimb and whole-body coordination (Figure 5, x-axis bar graph). The features contributing most strongly to LD2, which captured the differences in locomotor phenotype between the two mutants, included measurements of hindpaw movement (and subsequent changes in support patterns), variability, and relative phase of tail/nose movements (Figure 5, y-axis bar graph). Most variables relating to the averaged forward motion of individual limbs, which are largely spared in both mutants, do not strongly influence either LD.
The LDA provides a quantitative summary that captures the essential elements of the similarities and differences in locomotor behavior across groups of mice. In so doing, it shows that a high-dimensional set of gait parameters contains a subspace of mixed features in which multiple neuropathologies are represented behaviorally. The results, together with our previous description of the locomotor phenotype of pcd mice (Machado et al., 2015), reveal multijoint, interlimb, and whole-body coordination as shared, and possibly fundamental, features of locomotor ataxia.
Discussion
When assessing motor coordination deficits in mice there is often a tradeoff between specificity and interpretability. For example, approaches such as rotarod testing provide measurable and intuitive low dimensional outputs, but lack specificity. Others, such as the CatWalk system (Gabriel et al., 2009) can provide many detailed measurements of locomotor behavior, but their meaning is not always readily apparent. With LocoMouse we have tried to provide both a comprehensive, quantitative description of locomotor behavior as well as a conceptual framework within which to interpret that high-dimensional data.
Here we analyzed the locomotor behavior of severely ataxic reeler mice and quantitatively compared it with that of Purkinje cell degeneration mutants (Machado et al., 2015). Detailed comparison of locomotor kinematics and linear discriminant analysis revealed both shared and distinct features of gait ataxia in the two mouse lines. Although the generality of the conclusions that can be reached from the analysis of two mutants is of course limited, the approach described here provides a comprehensive and quantitative way to map complex patterns of locomotor features onto partially overlapping neuropathologies. Extending it to additional models and manipulations could enable association of specific movement features with increasingly precise alterations in underlying neural circuits.
Our first finding lies in capturing specific quantitative differences in gait ataxia between pcd and reeler mice. Rather than simple differences in the level of severity in the same set of affected features, the visible gait differences in the two mutants (captured by the second linear discriminant) appear to stem specifically from an increase in movement variability and additional hindlimb involvement (Cendelin, 2014) in reelers. As a likely consequence of those deficits, analysis of support patterns and nose and tail movements suggests that reelers, unlike controls and pcd, use their front paws as their main supports while walking. This difference in support patterns fully accounts for differences in tail movements observed in the two mutants (Figure 4). It could also explain the inability of reelers to walk in a straight line (Video 1, Figure 4H), likely because the front limbs are unable to provide both support and steering control to keep the body moving forward. Thus, while all of these features in combination contribute to the ability of human observers to visually distinguish the overall differences in walking patterns in the two mice, the quantitative LocoMouse analysis is able to distill them down to reveal fundamental underlying alterations in variability and hindlimb control.
Capturing the essential differences in locomotor control between ataxic phenotypes provides a necessary starting point for understanding the contributions of individual circuit elements to this complex, whole-body behavior. In pcd, the main anatomical phenotype is cerebellar, with striking degeneration of Purkinje cells that effectively removes cerebellar cortical input to the cerebellar nuclei. In contrast, reeler mice exhibit aberrant cell localization throughout the brain (Terashima et al., 1983; Stanfield et al., 1979; Mikoshiba et al., 1980; Blatt and Eisenman, 1988; Kang et al., 2010), which could also contribute to their gait abnormalities. Moreover, in reeler, abnormal cell migration during development results in aberrant circuit wiring that alters, rather than removes, Purkinje cell activity (Dupont et al., 1983; Curran and D'Arcangelo, 1998; D'Arcangelo et al., 1995; Lalonde and Strazielle, 2019).
Locomotor differences between the two mutants (captured by LD2 in Figure 5) thus likely represent a combination of specific differences in circuitry within the cerebellum itself plus extracerebellar effects. The differences in variability may be of cerebellar origin (Walter et al., 2006; Medina and Lisberger, 2007). Movement variability is often considered to be a hallmark of cerebellar ataxia, and yet surprisingly, pcd mice were less variable than controls on several movement measures (Machado et al., 2015). We speculate that the opposite effects on variability in the two mutants may arise as a consequence of altered Purkinje cell output in reeler, compared to a lack of Purkinje cell output in pcd (Medina and Lisberger, 2007). Meanwhile, differences in hindlimb involvement, and the resulting reliance on front paw supports (Figure 3G,H), may be extracerebellar, perhaps resulting from cell positioning defects in reeler spinal cord (Yip et al., 2000; Phelps et al., 2002). Future experiments analyzing locomotor phenotypes in additional cerebellar and extracerebellar models will be crucial for establishing mechanistic links between specific features of locomotor behavior and precise alterations in underlying neural circuits.
Despite the more severe behavioral and anatomical phenotype in reeler, we found the overall pattern of affected locomotor features in pcd and reeler to be surprisingly similar. The first linear discriminant pulled out a set of shared impairments in multi-joint, interlimb, and whole-body coordination that may reflect core quantitative features of mouse cerebellar ataxia. We propose that these features comprise the gestalt impression of clearly cerebellar gait phenotypes that trained observers can readily identify, despite substantial variation in specific manifestations across models.
In particular, like pcd (Machado et al., 2015), the tail and nose movements of reelers were also successfully modeled as passive consequences of limb movement, with a shift of supports to the front, rather than hind, paws. We had previously interpreted the pattern of coordination deficits in pcd mice, and particularly the passive tail oscillation, as consistent with the lack of an internal forward model that predicts and compensates for the consequences of movement (Bastian et al., 1996; Ebner and Pasalar, 2008; Ito, 2008; Kennedy et al., 2014; Wolpert et al., 1998). Specifically for locomotor control, a forward model can predict how the movement of one part of the body (e.g., a limb or a paw) will affect the movement of another part (e.g., another limb, or the tail), and inject a compensatory control signal to eliminate those consequences.
The idea that the cerebellum could provide a forward model for motor control is often invoked when considering the role of the cerebellum in motor learning (Ito, 2008; Shadmehr and Krakauer, 2008; Wolpert et al., 1998). Recently, we showed that both pcd and reeler mice were unable to learn to restore gait symmetry on a locomotor learning task that requires predictive control (Darmohray et al., 2019; Morton and Bastian, 2006; Reisman et al., 2005). That finding, together with the passive tail oscillations of both mutants, reinforces the idea that these animals lack the ability, provided by an intact cerebellum, to predict and compensate unintended movements during locomotion. The consistent lack of compensatory predictive mechanisms across mutants and behavioral paradigms suggests that they may represent core features of cerebellar contributions to coordinated locomotion. Extending the current findings and applying a similar approach across a broader range of circuit manipulations could become a key element in understanding how various neural components work together to control complex, whole-body behaviors.
Materials and methods
Animals
All procedures were reviewed and performed in accordance with the Champalimaud Centre for the Unknown Ethics Committee guidelines and approved by the Portuguese Direcção Geral de Veterinária (Ref. No. 0421/000/000/2015).
Heterozygous reeler (Relnrl) mice on a C57BL7/6 background were obtained from Jackson labs (#000235 B6C3Fe a/a-Relnrl/J). Data were collected from homozygous reeler mice (n = 9392 strides; N = 7 mice; two females; five males; 8–18 g; 35–52 days old; average weight = 13.57 ± 3.3 g) and their littermates (n = 9524 strides; N = 12 mice; eight females, four males; 12–25 g; 36–52 days old; average weight = 18.67 ± 3.5 g). Data from Purkinje cell degeneration mice (#0537 B6.BRAgtpbp1pcd/J) and additional wildtype C57BL7/6 (Jackson #000664) mice used for size-matching was previously collected and described in Machado et al., 2015. Mice were housed on a reversed light cycle 12 hr light/12 hr dark cycle, in standard cages with no more than five animals per cage. They had access to water and food ad libitum.
Experimental procedures
LocoMouse overground walking setup
Request a detailed protocolThe LocoMouse system for overground locomotion was composed by a glass corridor, 66.5 cm long and 4.5 cm wide with a mirror placed at 45 deg under the corridor. A single high-speed camera (AVT Bonito, 1440 × 250 pixels, 400 frames per second) recorded both bottom and side views of walking mice. Infrared sensors were used to automatically detect when mice entered and exited the corridor and trigger data collection, which was performed in LABVIEW 2012 (Machado et al., 2015).
Data collection
Request a detailed protocolAnimals were handled and acclimated to the overground setup for several sessions before data collection. Individual trials consisted of single crossings of the corridor, in which mice walked freely between two dark boxes. No food or water restriction or reward was used. Ten to twenty-five corridor trials were collected in each session for five consecutive days. An average of n = 1342 ± 637 strides were collected per reeler mouse (n = 348 ± 116 strides per paw) and n = 635 ± 341 strides per littermate mouse (n = 162 ± 86 strides per animal per paw) were collected.
To track the paws, nose and tail of locomoting mice we used the previously described noninvasive, markerless LocoMouse tracking system (Machado et al., 2015; https://github.com/careylab/LocoMouse), with additional, subsequent updates to the tail tracking, in order to handle the more erratic tail movements of the reeler mice, which often left the field of view of the videos. The new tail tracking algorithm was implemented using Matlab and the Signal Processing, Image Processing and Statistics and Machine Learning toolboxes. Tail tracking started with the side view using binary thresholding followed by a skeletonization operation for finding candidate tail segments based on length and position. These points were then projected onto the bottom view starting from the distal tail segment. The bottom view image was convoluted with a hamming window with a kernel width representative of a mouse tail. Subsequent points were identified iteratively towards the proximal tail, stopping at the base of the tail. In the case of tail segmentation after skeletonization the additional step of looking for tail points towards the distal tail was taken. The tail was then divided into 15 tail points (referred to as segments) with constant Euclidian distance (in 3D) between them, similarly to Machado et al., 2015. Matlab code for the updated tail tracker is available at https://github.com/careylab/LocoMouse_Dev.
Data analysis and statistics
Matlab 2012b and 2015a were used to process and analyze the data. Paw, nose and tail tracks (x,y,z) were obtained from the LocoMouse tracker (Machado et al., 2015). All tracks were divided in strides cycles. Stride cycle was defined as the period from stance onset to subsequent stance onset. For each stride, average walking speed was calculated. All data were sorted into speed bins (0.05 m/s bin width). Individual limb movements and interlimb coordination were calculated as follows:
Individual limb parameters
Request a detailed protocolWalking speed: x displacement of the body center during that stride divided by the stride duration.
Stride duration: time between two consecutive stance onsets.
Cadence: inverse of stride duration.
Swing velocity: x displacement of single limb during swing phase divided by swing duration.
Stride length: x displacement from touchdown to touchdown of single limb.
Stance duration: time in milliseconds that foot is on the ground during stride.
Duty factor: stance duration divided by stride duration.
Trajectories: (x,y,z) trajectories were aligned to swing onset and resampled to 100 equidistant points using linear interpolation. Interpolated trajectories were then binned by speed and the average trajectory was computed for each individual animal and smoothed with a Savitzky-Golay first-order filter with a 3-point window size.
Instantaneous swing velocity: the derivative of swing trajectory.
Model predictions: equations that were previously generated with mixed-effects models (Machado et al., 2015) to predict basic stride parameters.
Interlimb and whole-body coordination parameters
Request a detailed protocolBase of support: width between the two front and two hind paws during stance phase.
Body y displacement: y displacement of the body center during that stride.
Stance phase: relative timing of limb touchdowns to stride cycle of reference paw (FR). Calculated as: stance time - stance timereference paw/stride duration.
Supports: Support types were categorized by how many and which paws were on the ground, expressed as a percentage of the total stride duration for each stride. Paw support categories include 3-paw, 2-paw diagonal, 2-paw other/non diagonal (homolateral and homologous), and 2-paw front (only)supports.
Double support for each limb is defined as the percentage of the stride cycle between the touch down of a reference paw to lift-off of the contralateral paw. Because at higher speeds (running), the opposing limb lifts off before the reference paw touches down, we included negative double support by looking backwards in time, up to 25% of the stride cycle duration. Positive values of double support indicate that contralateral lift-off occurred after reference paw touch down, and negative values indicate that contralateral lift-off occurred before reference paw touch down. Note that front paw double support percentages include 2-paw front (only) support patterns as well as 3 and 4-paw support patterns in which both front paws were on the ground.
Tail and nose phases: For each speed bin we correlate the stridewise tail and nose trajectories with the trajectory given by the difference between the forward position of the right paw and the forward position of the left paw (also normalized to the stride). We do this both for front limbs (for the analysis of reeler mice) and hind limbs (for the analysis of previous pcd data). The phase is then calculated by the delay in which this correlation is maximized.
Tail and nose peak-to-peak amplitude: the change between peak (highest amplitude value) and trough (lowest amplitude value) in y or z during a stride duration.
Variability: All variability analyses were based on coefficients of variation (CV).
Geometric model of the tail and nose
Request a detailed protocolThe analytical model of the nose and tail was a simpler version of our previously described geometric model (Machado et al., 2015). The current model transforms the forward movements of the front limbs into predicted lateral oscillations of tail and nose. The model is described by the equation:
where yst and ynt are the lateral positions of tail segment s and nose, respectively; Gs and Gn are gains, obtained from fitting the data, that transform the limb oscillation amplitude to the amplitude of tail segment s and nose movements, respectively; xrt and xlt are the positions of front right and front left limbs at time t obtained from average trajectories of limb movements during strides at different speeds. Ds and Dn are the delays of tail segment s and nose, where the delay D1 is the delay of the base of the tail obtained by fitting the data.
As in Machado et al., 2015, delays between subsequent tail segments decreased according to the equation:
where s is the segment number, starting at the base of the tail. The delay of the nose was the same as the delay of the base of the tail (e.g. Dn = D1).
Principal component and linear discriminant analyses
Request a detailed protocolThe dataset consisted of a matrix of 109 observations of 45 features. Each observation was data from one mouse locomoting at a certain speed (binned) and features are z-scored gait parameters. LDA assumes independence within the feature space, which we knew to be violated due to the high speed-dependence of many gait features (Machado et al., 2015). Therefore, PCA was applied to address inter-variable correlation and avoid overfitting in LDA. PCA was performed by eigenvalue decomposition of the data covariance matrix. The first 10 PCs explained 88% of the variance and the data projected onto these 10 PCs were used as input to the LDA. The end contributions of the initial gait parameters to the two LD axes were obtained by multiplying the PCA mapping by the LDA mapping. LDA output is presented for each speed bin to verify that the pattern of differences across groups was captured, across all walking speeds.
Statistical analyses
Request a detailed protocolAll statistics can be found in Table S1. Statistical analyses were done in Matlab with the Statistics toolbox. An independent samples t-test was used to test for differences in walking speed distributions (Figure 1D and Figure 2A). For all other gait parameters, analysis was performed on animal averages binned by speed using mixed-effects models (Bates et al., 2013). Fixed-effects terms included speed and genotype; animals were included as random terms. We report F statistics from mixed ANOVAs with Satterthwaite degrees of freedom correction. Differences were considered significant at *p < 0.05, **p < 0.01, and ***p < 0.001; asterisks report main effects of genotype.
Data availability
Source data have been provided for all Figures.
References
-
The cerebellum and control of rhythmical movementsTrends in Neurosciences 6:417–422.https://doi.org/10.1016/0166-2236(83)90191-1
-
Cerebellar ataxia: abnormal control of interaction torques across multiple jointsJournal of Neurophysiology 76:492–509.https://doi.org/10.1152/jn.1996.76.1.492
-
Softwarelme Linear mixed-effects models using Eigen and S4, version 1.0-5R Package Version.
-
The need for speed in rodent locomotion analysesThe Anatomical Record 297:1839–1864.https://doi.org/10.1002/ar.22955
-
Topographic and zonal organization of the olivocerebellar projection in the Reeler mutant mouseThe Journal of Comparative Neurology 267:603–615.https://doi.org/10.1002/cne.902670412
-
Tests to assess motor phenotype in mice: a user's guideNature Reviews Neuroscience 10:519–529.https://doi.org/10.1038/nrn2652
-
Ataxic gait analysis in a mouse model of the olivocerebellar degenerationBehavioural Brain Research 210:8–15.https://doi.org/10.1016/j.bbr.2010.01.035
-
From mice to men: lessons from mutant ataxic miceCerebellum & Ataxias 1:4.https://doi.org/10.1186/2053-8871-1-4
-
Impaired classical eyeblink conditioning in cerebellar-lesioned and purkinje cell degeneration (pcd) mutant miceThe Journal of Neuroscience 16:2829–2838.https://doi.org/10.1523/JNEUROSCI.16-08-02829.1996
-
Role of reelin in the control of brain developmentBrain Research. Brain Research Reviews 26:285–294.https://doi.org/10.1016/s0165-0173(97)00035-0
-
Cerebellum predicts the future motor stateThe Cerebellum 7:583–588.https://doi.org/10.1007/s12311-008-0059-3
-
The use of multiple measurements in taxonomic problemsAnnals of Eugenics 7:179–188.https://doi.org/10.1111/j.1469-1809.1936.tb02137.x
-
Locomotor deficits in the mutant mouse, lurcherExperimental Brain Research 66:271–286.https://doi.org/10.1007/BF00243304
-
The CatWalk method: assessment of mechanical allodynia in experimental chronic painBehavioural Brain Research 198:477–480.https://doi.org/10.1016/j.bbr.2008.12.018
-
Cerebellar disorganization characteristic of Reeler in scrambler mutant mice despite presence of reelinThe Journal of Neuroscience 17:8767–8777.https://doi.org/10.1523/JNEUROSCI.17-22-08767.1997
-
Overground locomotion in intact rats: interlimb coordination, support patterns and support phases durationActa Neurobiologiae Experimentalis 59:131–144.
-
Control of mental activities by internal models in the cerebellumNature Reviews Neuroscience 9:304–313.https://doi.org/10.1038/nrn2332
-
Decoding the organization of spinal circuits that control locomotionNature Reviews Neuroscience 17:224–238.https://doi.org/10.1038/nrn.2016.9
-
Cerebellar contributions to locomotor adaptations during splitbelt treadmill walkingJournal of Neuroscience 26:9107–9116.https://doi.org/10.1523/JNEUROSCI.2622-06.2006
-
Mechanisms of cerebellar gait ataxiaThe Cerebellum 6:79–86.https://doi.org/10.1080/14734220601187741
-
Evidence for a cell-specific action of reelin in the spinal cordDevelopmental Biology 244:180–198.https://doi.org/10.1006/dbio.2002.0580
-
Interlimb coordination during locomotion: what can be adapted and stored?Journal of Neurophysiology 94:2403–2415.https://doi.org/10.1152/jn.00089.2005
-
A computational neuroanatomy for motor controlExperimental Brain Research 185:359–381.https://doi.org/10.1007/s00221-008-1280-5
-
The organization of certain afferents to the Hippocampus and dentate gyrus in normal and reeler miceThe Journal of Comparative Neurology 185:461–483.https://doi.org/10.1002/cne.901850304
-
Increased gait variability in mice with small cerebellar cortex lesions and normal rotarod performanceBehavioural Brain Research 241:32–37.https://doi.org/10.1016/j.bbr.2012.11.034
-
Distribution and morphology of corticospinal tract neurons in reeler mouse cortex by the retrograde HRP methodThe Journal of Comparative Neurology 218:314–326.https://doi.org/10.1002/cne.902180307
-
Cerebellar control of gait and interlimb coordinationBrain Structure and Function 220:3513–3536.https://doi.org/10.1007/s00429-014-0870-1
-
Internal models in the cerebellumTrends in Cognitive Sciences 2:338–347.https://doi.org/10.1016/S1364-6613(98)01221-2
Article and author information
Author details
Funding
H2020 European Research Council (Starting Grant 640093)
- Megan R Carey
Howard Hughes Medical Institute (International Early Career Scientist 55007413)
- Megan R Carey
Fundação para a Ciência e a Tecnologia (FCT-PTDC/MED-NEU/30890/2017)
- Megan R Carey
Fundação para a Ciência e a Tecnologia (SFRH/BD/51210/2010)
- Ana S Machado
Fundação para a Ciência e a Tecnologia (SFRH/BPD/119404/2016)
- Hugo G Marques
Fundação para a Ciência e a Tecnologia (PD/BD/141643/2018)
- Diogo F Duarte
Fundação para a Ciência e a Tecnologia (SFRH/BD/86265/2012)
- Dana M Darmohray
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank Tracy Pritchett for maintenance of mouse lines and help with histology and imaging and João Fayad for optimization and maintenance of the LocoMouse tracker. Olivia Carmo assisted with data acquisition and Alex Azinheira made the video. We are grateful to the members of Carey Lab and the Champalimaud Neuroscience Program for helpful discussion.
This work was supported by a Howard Hughes Medical Institute International Early Career Scientist Grant #55007413 (to MRC), European Research Council Starting Grant #640093 (to MRC), FCT grant FCT-PTDC/MED-NEU/30890/2017 (to MRC), Congento LISBOA-01–0145-FEDER-022170, and fellowships from the Portuguese Fundação para a Ciência e a Tecnologia: SFRH/BD/51210/2010 (to ASM), SFRH/BD/86265/2012 (to DMD), PD/BD/141643/2018 (to DFD), and SFRH/BPD/119404/2016 (to HGM).
Ethics
Animal experimentation: All procedures were reviewed and performed in accordance with the Champalimaud Centre for the Unknown Ethics Committee guidelines and approved by the Portuguese Direcção Geral de Veterinária (Ref. No. 0421/000/000/2015).
Copyright
© 2020, Machado 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,500
- views
-
- 249
- 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
Time estimation is an essential prerequisite underlying various cognitive functions. Previous studies identified ‘sequential firing’ and ‘activity ramps’ as the primary neuron activity patterns in the medial frontal cortex (mPFC) that could convey information regarding time. However, the relationship between these patterns and the timing behavior has not been fully understood. In this study, we utilized in vivo calcium imaging of mPFC in rats performing a timing task. We observed cells that showed selective activation at trial start, end, or during the timing interval. By aligning long-term time-lapse datasets, we discovered that sequential patterns of time coding were stable over weeks, while cells coding for trial start or end showed constant dynamism. Furthermore, with a novel behavior design that allowed the animal to determine individual trial interval, we were able to demonstrate that real-time adjustment in the sequence procession speed closely tracked the trial-to-trial interval variations. And errors in the rats’ timing behavior can be primarily attributed to the premature ending of the time sequence. Together, our data suggest that sequential activity maybe a stable neural substrate that represents time under physiological conditions. Furthermore, our results imply the existence of a unique cell type in the mPFC that participates in the time-related sequences. Future characterization of this cell type could provide important insights in the neural mechanism of timing and related cognitive functions.
-
- Neuroscience
Granule cells of the cerebellum make up to 175,000 excitatory synapses on a single Purkinje cell, encoding the wide variety of information from the mossy fibre inputs into the cerebellar cortex. The granule cell axon is made of an ascending portion and a long parallel fibre extending at right angles, an architecture suggesting that synapses formed by the two segments of the axon could encode different information. There are controversial indications that ascending axon (AA) and parallel fibre (PF) synapse properties and modalities of plasticity are different. We tested the hypothesis that AA and PF synapses encode different information, and that the association of these distinct inputs to Purkinje cells might be relevant to the circuit and trigger plasticity, similar to the coincident activation of PF and climbing fibre inputs. Here, by recording synaptic currents in Purkinje cells from either proximal or distal granule cells (mostly AA and PF synapses, respectively), we describe a new form of associative plasticity between these two distinct granule cell inputs. We show for the first time that synchronous AA and PF repetitive train stimulation, with inhibition intact, triggers long-term potentiation (LTP) at AA synapses specifically. Furthermore, the timing of the presentation of the two inputs controls the outcome of plasticity and induction requires NMDAR and mGluR1 activation. The long length of the PFs allows us to preferentially activate the two inputs independently, and despite a lack of morphological reconstruction of the connections, these observations reinforce the suggestion that AA and PF synapses have different coding capabilities and plasticity that is associative, enabling effective association of information transmitted via granule cells.