Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Rapid and robust patterns of spontaneous locomotor deficits in mouse models of Huntington’s disease

  • Taneli Heikkinen ,

    Contributed equally to this work with: Taneli Heikkinen, Timo Bragge

    Roles Investigation, Methodology, Supervision, Writing – original draft, Writing – review & editing

    Affiliation Charles River Discovery Services, Kuopio, Finland

  • Timo Bragge ,

    Contributed equally to this work with: Taneli Heikkinen, Timo Bragge

    Roles Conceptualization, Data curation, Formal analysis, Investigation, Software, Writing – original draft, Writing – review & editing

    Affiliation Charles River Discovery Services, Kuopio, Finland

  • Niina Bhattarai,

    Roles Investigation

    Affiliation School of Pharmacy, Faculty of Health Sciences, University of Eastern Finland, Kuopio, Finland

  • Teija Parkkari,

    Roles Project administration

    Affiliation Charles River Discovery Services, Kuopio, Finland

  • Jukka Puoliväli,

    Roles Resources

    Affiliation Charles River Discovery Services, Kuopio, Finland

  • Outi Kontkanen,

    Roles Project administration

    Affiliation Charles River Discovery Services, Kuopio, Finland

  • Patrick Sweeney,

    Roles Resources

    Affiliation Naason Science Inc., Chungcheongbuk-do, South Korea

  • Larry C. Park,

    Roles Conceptualization, Funding acquisition, Supervision, Writing – original draft

    Affiliations Naason Science Inc., Chungcheongbuk-do, South Korea, CHDI Management/CHDI Foundation, Los Angeles, California, United States of America

  • Ignacio Munoz-Sanjuan

    Roles Funding acquisition, Resources, Supervision, Writing – original draft, Writing – review & editing

    ignacio.munoz@chdifoundation.org

    Affiliation CHDI Management/CHDI Foundation, Los Angeles, California, United States of America

Abstract

Huntington's disease (HD) is an inherited neurodegenerative disorder characterized by severe disruption of cognitive and motor functions, including changes in posture and gait. A number of HD mouse models have been engineered that display behavioral and neuropathological features of the disease, but gait alterations in these models are poorly characterized. Sensitive high-throughput tests of fine motor function and gait in mice might be informative in evaluating disease-modifying interventions. Here, we describe a hypothesis-free workflow that determines progressively changing locomotor patterns across 79 parameters in the R6/2 and Q175 mouse models of HD. R6/2 mice (120 CAG repeats) showed motor disturbances as early as at 4 weeks of age. Similar disturbances were observed in homozygous and heterozygous Q175 KI mice at 3 and 6 months of age, respectively. Interestingly, only the R6/2 mice developed forelimb ataxia. The principal components of the behavioral phenotypes produced two phenotypic scores of progressive postural instability based on kinematic parameters and trajectory waveform data, which were shared by both HD models. This approach adds to the available HD mouse model research toolbox and has a potential to facilitate the development of therapeutics for HD and other debilitating movement disorders with high unmet medical need.

Introduction

Mammalian locomotion is composed of a sequence of musculoskeletal positions to exert force [1,2] and, when compromised, organisms can adapt and uncouple these components to prioritize efficiency [35]. These adaptations for the economy of movement can be measured as changes in gait. In patients, the degeneration of critical brain regions for movement or posture control, e.g., basal ganglia, leads to progressive alterations in gait [610]. Neuroprotection for such disorders is likely to be most successful when targeted to the early phases of neural dysfunction before significant neuronal loss, placing an emphasis on the development of disease models that demonstrate early and progressive phenotypes in movement and gait. Such phenotypes have been observed in genetic, toxin-induced and surgical models of neurological disorders [1117]. However, measuring quantitative changes in the complex adaptive movement of aging diseased mice is challenging.

Huntington's disease (HD) is a progressive neurodegenerative disorder caused by an expansion of the CAG trinucleotide repeat region in the huntingtin (HTT) gene on chromosome 4 [18]. Although HD is characterized by multiple movement, psychiatric and cognitive disturbances, its diagnosis is ultimately confirmed on the basis of the positive genetic test and onset of motor symptoms as measured by the Unified HD Rating Scale (UHDRS) total motor score (TMS) [19,20]. UHDRS-TMS comprises several criteria, including the presence of bradykinesia and altered gait pattern [21]. Bradykinesia—expressed as decreased velocity and stride length—can be detected in HD gene-expansion carriers (HDGECs) already in the presymptomatic phase [22,23]. As disease progresses, these manifestations worsen in HDGECs, and their cadence, joint angle range and stride length decrease further [24,25]. Significantly increased intra-individual variability of gait parameters such as stride length, stance time and swing time are typical features of HD [2527].

Multiple genetically engineered mouse models expressing mutated HTT have been developed that show aspects of HD-like pathology and neurological symptoms. These mouse lines can be divided broadly into three types based on the mutant HTT (mHTT) construct: transgenic models expressing either N-terminal fragments of the mHTT gene or full-length mHTT, and knock-in (KI) models, in which CAG repeats are inserted into the endogenous mouse Htt gene (for review, see [28]). Neurological symptoms have been most extensively studied in HD mouse models using tests such as rotarod, open field, grip strength and footprint test for the evaluation of gross cadence changes [29]. Alterations in several gait parameters in these mouse models could be correlated to the observations in patients. For example, shorter stride length was noted in 24-month-old HdhQ111 KI heterozygous and homozygous mice and 12-month-old homozygous Hdh(CAG140) mice [30,31]. The mean stride and base lengths did not differ in 12-month-old heterozygous Hdh(CAG)150 KI mice and wild-type counterparts but the intra-individual variability of these indices was considerably higher in the mutants [32]. Being potentially useful for phenotypic and pharmacological screening, the footprint test nonetheless does not characterize the full spectrum of gait deficits of HD mouse models. Moreover, in some widely used HD lines, gait changes revealed by the footprint test are minor and transient (e.g., as in BACHD mice [29]) and/or appear very late in life (13.5 months in YAC128 line [33] and 24 months in the HdhQ111 KI line [30]).

To provide a more robust method of assessing locomotor behavior relevant to HD symptoms, we have developed a novel kinematic analysis approach based on our previous work [13] and tested it in two popular HD models, R6/2 and Q175. R6/2 mice express exon 1 of mHTT carrying a 141–157 CAG repeat, which triggers rapid progression of the pathology: the mutant mice start to lose body weight and develop motor and cognitive deficits at 6–8 weeks of age and die early [34]. Q175 mice have the human exon 1 sequence with approximately 190 CAG repeats knocked into the endogenous Htt gene; they show age- and genotype-dependent progression of motor, cognitive, and electrophysiological deficits, accompanied by decreased brain volumes and striatal metabolic changes, which develop more slowly than in the R6/2 mice [35,36]. Motor changes in Q175 mice typically present as hypoactivity in homozygotes from 2 months and in heterozygotes from 4 months of age [36]. We selected R6/2 and Q175 mice to test our new approach for the following reasons. R6/2 mice are by far the most popular HD mouse model with a clear gait phenotype revealed by different tests [13,29,37,38], therefore it could be a suitable benchmark. In contrast to the situation with R6/2 mice, motor symptoms in Q175 KI mice progress over several months, particularly in heterozygous animals, and their gait has not been analyzed in detail previously. Judging from at best moderate gait phenotypes in other KI HD mouse models [2932], Q175 mice would provide a good sensitivity test of our method. In our high precision kinematic analysis, we measured 79 parameters of gait, body posture and fine motor movements during early, middle and advanced phases of the disease in R6/2 and Q175 KI mice of both sexes. Five principal components of mouse movements were identified that captured main features of mouse gait. In addition, two generalized phenotypic scores based on kinematic parameters and trajectory waveform patterns, respectively, have been derived. We expect that these indices can be more broadly applied to mouse models of other movement disorders to determine treatment effects and potentially be utilized in neurotoxicology applications.

Materials and methods

Animals

Experiments were conducted at Charles River Discovery Services, Kuopio, Finland according to the National Institute of Health (NIH) Guidelines for the Care and Use of Laboratory Animals and approved by the State Provincial Office of Southern Finland. R6/2 and respective wild-type (WT) littermates (120 CAGs) were obtained from the Jackson Laboratory (JAX Stock No: 006494, The Jackson Laboratory, Bar Harbor, ME, USA) [34]. Heterozygous Q175 mice (strain CHDI-81003003 (neo+), JAX Stock No: 027410) were intercrossed to produce homozygous (HOM), heterozygous (HET) Q175 mice and wild-type (WT) littermates (Charles River Laboratories, Sulzfeld, Germany). PCR confirmed the genotype of each mouse analyzed (Charles River Finland and Laragen Inc., Culver City, CA, USA). Mice were housed in groups of up to 5 per cage, in a temperature (22 ± 1°C) and humidity (30–70%) controlled environment with a normal light-dark cycle (7:00–20:00). The group sizes at different ages are described in Table 1. Both females and males were included, and all cohorts comprised a comparable (±1) number of animals of each genotype per respective age group (Table 1). There were no significant differences between sexes, and the data are shown for groups with animals of both sexes combined. All mice were housed in cages (dimensions: length 35 cm × width 19 cm × height 13 cm) with clean wood shavings covering the ground, changed weekly to provide the animals with dry bedding. A red mouse igloo was placed in each cage to provide environmental enrichment and shelter. Food (Purina Lab Diet 5001) and water were available ad libitum (Table 1).

thumbnail
Table 1. Number, age and sex of R6/2 and Q175 KI mice in experimental cohorts.

https://doi.org/10.1371/journal.pone.0243052.t001

Fine motor skill and gait analysis

Fine motor skills and gait of the mice were captured by a high-speed camera (300 fps) under a brightly illuminated plexiglas corridor (153 × 5 × 10 cm, Motorater, TSE-systems GmbH, Bad Homburg, Germany). A few days before the test sessions, under light isoflurane anesthesia, the fur of the limbs was shaved. On the day of testing, the joints and tail of each mouse were highlighted with a non-toxic paint. The mice were analyzed while walking along the corridor and five or six complete strides were analyzed for each mouse. Only strides with continuous ambulatory movement were analyzed. Mirrors enabled the performance of the mouse to be detected simultaneously from left, right and ventral aspects. The movement of 24 points on each mouse were analyzed from the videos (Simi Reality Motion Systems, Unterschleissheim, Germany), and the trajectories of markers (see Table 2 for details) were analyzed by customized scripts.

thumbnail
Table 2. Anatomical landmarks marked and tracked in each mouse.

https://doi.org/10.1371/journal.pone.0243052.t002

Because the progression of behavioral changes in the transgenic mice was strain-specific, the time of each measurement was selected based on the results of standard behavioral tests of motor function in previous studies [3436]. Specifically, we selected the initial ages of 4 weeks for R6/2 and Q175 HOM mice, and 3 months for Q175 HET, because those were the age points at which each model displays either mild or no motor deficits. We selected the older ages of 10 weeks for R6/2 mice, 3 and 6 months for Q175 HOM, and 6 and 10 months for Q175 HET animals because for each respective strain, there is documented temporal progression from mild to moderate and further, to severe stages of motor symptoms. Motor deficits are clearly apparent at 10 weeks in R6/2 mice, and at 6 months, they can be classified as advanced in Q175 HOM mice and as moderate in Q175 HET mice [36]. The data were combined from two or three experimental groups for each genotype, each experiment consisting of measurements at one or two age points from Q175 HET, Q175 HOM or R6/2 mice and the corresponding age-matched wild-type (WT) littermates.

Measuring individual components of a complex spontaneous behavior

The components of spontaneous mouse motor activity were carefully defined. A full locomotor cycle, a stride, was defined as the period between two consecutive left hindlimb floor contacts. Parameter determination procedure was performed as follows: first, the initial ground contacts and the onsets of the swing phases for each limb were detected and strides were determined. Then, spatio-temporal parameters (e.g., stance and swing time, stride distance and mean speed), as well as parameters measuring inter-limb coordination, body posture and joint angles, and properties of limb trajectories during swing phase were determined. Altogether, 79 different parameters were established (Table 3).

Strides were excluded from analysis if the animal stopped for any reason, or if the stride was unusually slow (slower than 50% of the median speed). The experimental setup measured both left, right and ventral aspects of each mouse simultaneously. All data are presented as an average of both left and right-side movements for each mouse strain.

Data analysis

Principal component analysis (PCA) can measure the interdependence of kinematic parameters of the same movement and reduce the dimensionality of multivariate data set [3943]. First, 26 principal components of z-score normalized kinematic gait parameter data were obtained, corresponding to 90% of total variance in the 63 kinematic parameters. After Varimax rotation [44], five components (PC# 1, 2, 3, 5, and 11) were chosen based on the explained variance and the ability to differentiate the three genotypes from respective WT mice. The corresponding PC scores were obtained first for each stride of each individual, and finally, the mean PC scores were calculated over the strides.

To gain a comprehensive perspective on the spatio-temporal features of mouse gait, we established a method for marker trajectory data analysis and reconstruction to capture any correlated associations between different marker trajectories. The method is based on the PCA of an empirical data correlation matrix of marker trajectory data (also known as Karhunen-Loève transform) [41] and is independent of the subjectively pre-selected set of kinematic parameters. Similar approaches utilizing kinematic marker trajectory (waveform) data have been presented earlier [39,45,46]. In our approach, two-dimensional marker data of 22 selected markers of each stride (bilaterally: iliac crest, hip, knee, ankle, heel, hindpaw, shoulder, elbow, wrist and forepaw as well as single points of the chin and tail tip, middle, and base) were first referenced to the left iliac crest marker. The marker data were temporally normalized to 50 data points by interpolation, corresponding to 0–100% of the stride duration, adding up to 2,000 points per data vector. Empirical data correlation matrix (XTX/n) was formed from the data matrix X, consisting of data vectors from each stride of each mouse, and 50 PCs were computed. Finally, the PC scores were obtained and further incorporated into the phenotypic score.

Fine motor phenotype scores

The phenotypic score was determined as follows: first, the average PC score of the WT group was subtracted from all PC scores. Second, the PC scores for each individual mouse were orthogonally projected onto a unit direction vector of the multidimensional line starting from the average of WT (zero after subtraction) and pointing towards R6/2 or Q175 groups within the PC space. Finally, the Phenotypic score was obtained based on the distances from WT. By definition, the average control group individual scored zero. The average disease model phenotype has a positive score. The magnitude of that value corresponds to the magnitude of phenotype-specific deviations in the overall gait pattern.

This approach used 63 separate parameters and trajectory waveform patterns to produce two separate overall phenotypic scores for each mouse. The unrotated eigenvectors were used to maximize the information available. The dimension of the PC score space was 26 and 50 for the parameter and marker trajectory data based on projected distances, respectively. Because there was also notable age progression in many parameters, the multidimensional lines connecting WT and TG group averages were determined separately for each age.

Statistical analysis

Each of the 79 parameters was analyzed separately for every mouse strain at each age by using standard software (IBM SPSS Statistics 19 and StatsDirect). The comparisons of mutant genotype groups to the corresponding WT were performed by the unpaired t-test or the Mann-Whitney U-test. All values are presented as the mean ± standard error of the mean (SEM). Differences were considered significant when p < 0.05.

Results

Individual parameters of gait and fine motor features

Each mouse was assessed across 79 movement-related parameters per time point (Table 3). The overall timing and speed of motion were altered in R6/2 mice from 4 weeks of age. The stride time was longer, but the stride distance was shorter at 10 weeks, reflecting a longer stance time in both the hindlimbs and forelimbs compared to that in WT mice (p < 0.05, Table 4). In the forelimbs, the peak swing speed was slower at 4 weeks of age (p < 0.05, Table 4). In Q175 mice, the mean speed as well as the mean and peak swing speeds of the forepaws were slower from 6 months in the HET mice and from 3 months in the HOM mice (p < 0.05, Table 4). Furthermore, an ataxic movement at the peak of the swing phases was observed in the forepaws of R6/2 mice at 10 weeks of age. Representative example videos of R6/2 and corresponding WT mice at 10 weeks of age as well as of Q175 WT, HET and HOM mice at 6 months of age are included in Supporting Information (S1 Appendix).

thumbnail
Table 4. Summary of fine motor and gait deficits of R6/2 and Q175 mice during walking.

https://doi.org/10.1371/journal.pone.0243052.t004

The cadence and interlimb coordination in R6/2 mice were significantly affected at 10 weeks of age. R6/2 mice took wider steps using all limbs, and the hindlimb step width was more variable (p ≤ 0.01; Table 4; Fig 1A and 1B). The diagonal cadence, the dominant form of cadence in normal healthy rodents, was substituted by the homolateral cadence (p < 0.05; Table 4; Fig 1E and 1F). Q175 mice also had deficits in the cadence and interlimb coordination but to a lesser extent than R6/2 mice. Specifically, hindlimb steps were wider in HET mice from 6 months and in HOM mice at 3 months (Fig 1C and 1D; p < 0.05). Forelimb steps were wider only in HOM mice at 6 months of age (p < 0.05, Table 4), with smaller diagonal cadence in HET mice at 10 months and in HOM mice at 6 months (p < 0.05; Fig 1G and 1H; Table 4).

thumbnail
Fig 1. Gait phenotypes of R6/2 and Q175 mice during spontaneous motor activity.

Histograms of the measurements from R6/2 (dark blue), Q175 heterozygous (HET, green), Q175 homozygous mice (HOM, purple) and wild-type (WT) littermates (turquoise) of step width (A–D), diagonal inter limb coordination (E–H) and height of tail tip (I–L) over time. Data are presented as the mean ± SEM. Statistical significance of the differences is indicated as follows: *p < 0.05, **p < 0.01, ***p < 0.001 (R6/2 / Q175 KI HET / HOM vs. WT, unpaired t-test). Information about the numbers of each group at specific ages is given in Table 1.

https://doi.org/10.1371/journal.pone.0243052.g001

Pelvic height was significantly lower in both R6/2 and Q175 KI mice. In the R6/2 mice, this was expressed, in particular, as decreased iliac crest height at mid-stance and mean hip height as well as increased protraction from 4 weeks of age (p < 0.05; Table 4). In the Q175 mice, iliac crest height and hip height were both decreased in HET mice from 6 months and in HOM mice from 3 months (p < 0.05; Table 4). Hindlimb protraction was increased from 6 months in both HET and HOM mice (p < 0.05; Table 4).

Tail tip position and movement of the tail were among the most robust motor changes in both R6/2 and Q175 mice. Tail tip height was lower at 10 weeks and the minimum height already from 4 weeks in R6/2 mice compared to the values in WT mice (p < 0.05; Fig 1I and 1J; Table 4). In addition, the tail tip spent longer on the ground, and the vertical and three-dimensional movement of the tail tip was increased during walking at 10 weeks (p < 0.05, Table 4). Tail tip heights were lower in HET mice at 6 months, and in HOM mice from 3 months (p < 0.05; Fig 1K and 1L). Similarly, the tail tips moved more in two and three dimensions at the same ages in both HET and HOM Q175 mice (p < 0.05; Table 4). The tail tip spent less time over the hip in both HET and HOM Q175 mice at 6 months (p < 0.05, Table 4).

In R6/2 mice, the range of motion in the knee was broader from 4 weeks to 10 weeks of age (p < 0.05; Table 4) and more variable in the hip (p < 0.05), knee (p < 0.05 at 10 weeks) and hind ankle (p < 0.05 at 4 weeks; Table 4). In the Q175 mice, the range of motion was broader in the ankle at 10 months in HET (p < 0.05) and at 6 months in HOM mice (p < 0.05). Further, it was more variable in the hip (p < 0.05) in HET mice at 10 months and in HOM mice at 3 months (p < 0.05), and for the ankle from 6 months in both HET and HOM mice (p ≤ 0.05; Table 4).

The trajectory shapes of front and hind paws during the swing phase of each stride, were significantly compromised in both R6/2 and Q175 mice. Lower trajectory profile of forepaws was observed in R6/2 mice at 4 weeks, in HET Q175 mice at 10 months and in HOM Q175 mice at 3 months of age (p < 0.05; Table 4). The major characteristics of the mouse lines—movement speed, interlimb coordination, pelvic height and paw trajectory shape—deteriorated in the groups of older mice (Table 4).

Kinematic parameter PCA

From the 63 interrelated movement parameters measured for each mouse, five principal components represented CAG repeat-associated pathological features: PC#1—overall slowness; PC#2—stride distance, swing speed; PC#3—pelvic height; PC#4—tail tip height; PC#5—swing trajectory profile (Figs 2 and 3). All transgenic mice moved slower than WT littermates (R6/2 from 4 weeks, Q175 HET from 6 months, and Q175 HOM from 3 months; p < 0.05 for PC#1 in all groups). R6/2 mice at 10 weeks, Q175 HET mice at 6 months, and Q175 HOM mice from 3 months old had slower swing speeds and shorter strides than WT mice (p < 0.05 for PC#2 in all groups). Pelvic height was the component that most consistently differentiated HD mice from their WT littermates. Both R6/2 and older Q175 mice walked with their body closer to the ground than WT littermates (p < 0.05 for PC#3 in all groups), although for the youngest Q175 age groups (3-month-old HET mice and 1-month-old HOM mice), no difference was observed. In addition to hip, iliac crest and tail base height parameters, the body posture PC#3 also represented inversely correlated knee and ankle ranges of movements and pro/retractions, meaning that a reduced PC#3 indicated both lower body posture and increased vertical movement, or worse body control.

thumbnail
Fig 2. Fine motor and gait deficits using PCA.

Principal component (PC) scores PC#1–5 of the five selected Varimax-rotated PCs are illustrated. The corresponding PCs (eigenvectors) are shown in Fig 3. Percentage in each panel describes the proportion of variation in the whole data set that each PC comprises. Data are shown separately for R6/2 (A), Q175 KI HET (B), and Q175 KI HOM mice (C). Data are presented as the mean ± SEM. Statistical significance of the differences is indicated as follows: *p < 0.05, **p < 0.01, ***p < 0.001 (R6/2 / Q175 KI HET / HOM vs. WT, unpaired t-test). Information about the numbers of each group at specific ages is given in Table 1.

https://doi.org/10.1371/journal.pone.0243052.g002

thumbnail
Fig 3. Principal components of gait phenotypes in R6/2 and Q175 mice.

The five Varimax-rotated principal components are presented as a heat map, illustrating the interdependence of kinematic parameters within each component (red = positive correlation, blue = negative correlation, black = no correlation).

https://doi.org/10.1371/journal.pone.0243052.g003

PC#4 represented the height of the tail tip, which was significantly lower in R6/2 at 10 weeks (PCA; p < 0.05), in Q175 HET at 6 months (p < 0.05), and most distinctly in Q175 HOM from 3 months (p < 0.05). Overall, the tail tips of HD mice were closer to the ground than those of WT animals.

PC#5 represented the limb swing trajectory profile. HD mutants demonstrated lower hind limb trajectories and greater protraction. The maximal toe clearance was not different between HD and WT mice, but the paw trajectory shape was different in the mutants, namely, the hindpaw was closer to the ground during the swing phase. This principal component showed deficits in R6/2 mice at all ages (p < 0.05), in Q175 HET mice at 10 months (p < 0.05) and in Q175 HOM mice from 3 months (p < 0.05).

Finally, two Fine Motor Phenotype Scores were calculated by combining the kinematic parameters (Fig 4A–4C) and the interpolated marker trajectory data (Fig 4D–4F). R6/2 mice demonstrated altered gait parameter and marker trajectory data-based scores at both 4 and 10 weeks of age (t-test; p < 0.05) (Fig 4A and 4D). Similarly, Q175 HET and HOM mice had altered gait and marker trajectory scores from as early as 3 months and 1 month of age, respectively (t-test; p < 0.05 for all comparisons) (Fig 4B, 4C, 4E and 4F). Importantly, no single parameter was as sensitive as the Fine Motor Phenotype Scores for detecting early phenotypes in the movement of the mice: whereas no single parameter allowed for reliable differentiation of all age groups in the three cohorts (Table 4), both phenotypic scores were significantly different for mutants of all ages compared to the values in respective WT counterparts (Fig 4). Furthermore, power analyses demonstrated that the Fine Motor Phenotype Score based on the marker trajectory data required fewer mice to detect a prospective 50% treatment effect in young mice that the score based derived from the kinematic parameters (Table 5).

thumbnail
Fig 4. Gait and trajectory data phenotypic scores in R6/2 and Q175 mice worsen over time.

Graphs of the scores from each R6/2 (dark blue), Q175 heterozygous (green), Q175 homozygous mice (purple) and wild-type (WT) littermates (turquoise) after PCA for either 63 gait parameters (A–C) or 24 parameters of marker trajectory data (D–F) over time. Data are presented as the mean ± SEM. Statistical significance of the differences is indicated as follows: *p < 0.05, ***p < 0.001 (R6/2 / Q175 KI HET / HOM vs. WT, unpaired t-test). Information about the numbers of each group at specific ages is given in Table 1.

https://doi.org/10.1371/journal.pone.0243052.g004

thumbnail
Table 5. Power analysis for the fine motor phenotype scores computed from kinematic parameter and marker trajectory data.

https://doi.org/10.1371/journal.pone.0243052.t005

Discussion

Gait disturbances are important manifestations of HD, therefore gait features are typically assessed in standardized HD motor scores, such as UHDRS-TMS, which are frequently used in HD clinical trials [47]. Although various motor assessments are a conventional element of mouse HD model phenotyping, specific analysis of gait is seldom carried out. Here, to provide an integrative tool for drug discovery, we measured 79 locomotor parameters in each individual mouse to comprehensively characterize motor phenotypes in two genetically engineered HD mouse models during disease progression. From these individual constituents of mouse locomotion, we identified five central themes (principal components) and derived two functional phenotypic scores based on kinematic and marker trajectory indices that robustly quantified locomotor deficits. We showed that our Fine Motor Phenotype Scores more sensitively detected changes in locomotion in mutants than any single kinematic parameter and therefore, these scores may be applied for characterization of other HD mouse lines, as well as mouse models of other diseases associated with locomotor impairments.

Motor dysfunction in HD mouse models has been commonly detected by the open field and rotarod tests as part of preclinical evaluations of numerous potential HD therapeutics [48,49]. In addition, the utility of specific gait measurements in HD mouse models is becoming increasingly acknowledged, because gait alterations in HD patients may have different sensitivity to the available treatments compared to other motor manifestations, such as chorea [50,51]. Treatment-induced improvements of gait abnormalities in HD mouse models have been described in several published reports. We have recently observed a positive effect of intracerebroventricular infusions of an antisense oligonucleotide against the expanded CAG repeat in HTT mRNA on stride distance, nose height and peak swing speed of R6/2 mice [13]. Intrastriatal injections of a calmodulin fragment and dantrolene administration with food significantly increased the reduced stride length in R6/2 mice [52] and YAC128 mice [33], respectively. Intracerebroventricular infusion of ganglioside GM1 shortened stride duration, corrected stance‐to‐stride ratio and restored interlimb coupling in Q140 mice [53]. Notably, transplantation of embryonic stem cell-derived neural progenitors into the N171-82Q mouse model of HD rescued changes in stride length, print length and print area as revealed by the CatWalk gait assay, whereas conventional grip strength and rotarod tests failed to demonstrate positive effect of the treatment [54].

Progressive decreases in movement speed and lower stride distance observed by us in R6/2 and Q175 mice (Table 4) are in agreement with previous assessments of locomotor functions in these [13,37,38,5458] and other [2931,33] mouse models of HD. Here, we replicated our previous findings of lower stride speed, stride distance and peak swing speed in 10-week-old R6/2 mice [13] and extended some of them to the younger age of 4 weeks. Shorter stride distance was also noted in 10-week-old R6/2 mice in the CatWalk [54]. In contrast, experiments in the Digigait automated treadmill video capture system failed to reveal differences in stride distance, stride time or swing time in 17-week-old R6/2 mice and their WT counterparts although expected shorter stride distance was revealed in R6/2 mice in that study by using conventional footprint test [38]. Moreover, the stance time was shorter in the mutants in that study, whereas in our experiments, it was longer. These discrepancies are likely explained by the forced nature of treadmill testing and differential adaptation to it by WT mice and mutants at the advanced stage of the pathology [38].

In addition to comparing individual parameters, we also calculated two Fine Motor Phenotype Scores based, respectively, on the kinematic indices (Fig 4A–4C) and interpolated marker trajectory data (Fig 4D–4F). Similar approaches have measured patterns of gait changes in patients: scores such as the Gillette Gait Index and Gait Deviation Index are based on the principal components of pathological gait or kinematic waveform data, respectively, without predetermined parameters, and they can be used to assess treatment efficacy [3,39,59]. Recently, a computational approach that utilized support vector machines was applied to a large-scale analysis of a series of HD KI mouse line, including HET Q175 animals, to reveal behavioral features that would allow clear differentiation of mice with CAG repeats of different lengths [60]. Multiple spatio-temporal gait features of the KI lines were determined in the NeuroCube apparatus, where animals walked for 5 min, and the change in base width was among top 10 parameters that allowed to reliably differentiate six KI mouse lines with the number of CAG repeats ranging from 20 to 175 [60]. As detailed descriptions of the individual parameters and their changes were not explicitly reported in that study [60], it is impossible to compare adequately the sensitivity of the NeuroCube and our Motorater-based approaches to detect gait changes in Q175 mice. A potential advantage of the NeuroCube-based method is that it integrates more parameters over longer period of spontaneous walking. On the other hand, as it is not clear whether and how different anatomical parts were marked in that study, it is likely that our approach allowed better estimation of the movement trajectories.

In our present experiments, the strongest principal component of movement in both HD models was the slowing of the time between strides, consistent with similar findings in rodent models of other neurodegenerative diseases [6165], and people with either Huntington’s or Parkinson’s diseases [66]. Our approach allows us to measure the interactions between slower movements and other gait parameters, especially those related to stride tempo. Slower movement correlated with lowered hind limb paws and shorter steps, consistent with the data in the Hdh(CAG)150 KI HD mouse model [32]. Walking speed is linked to limb coordination [67], and the limbs of healthy mice are remarkably coordinated during walking [68]. Our data were consistent with evidence from mouse models of HD [54] as well as neurological disorders and traumatic brain injury [6975] in revealing decreased limb coordination in the mutants. In contrast, parameters related to the body posture–for example, the range of joint movement–were uncoupled from walking speed in our experiments, consistent with observations of R6/2 mice in a forced walking test [76] or during spontaneous movement [37]. In addition, the emergence of forelimb ataxia in the R6/2 mice is in accord with case reports of ataxia symptoms in HD patients [77,78], although mouse models of spinocerebellar ataxia exhibit different kinematic patterns [79,80]. Another important characteristic in both R6/2 and Q175 mice was the lowering of the rear body. The lower base of the tail was related to several other body posture parameters, including lower iliac crest, hip and hind limb paw clearance, and changes in the range of joint movements. Decreased height of the iliac crest was also reported in male (but not female) Q140 mice [54]. These data suggest postural instability in the mutant mice, consistent with clinical examinations of HD patients, for whom falls are a common problem [81,82]. Further studies are required to determine whether these movement phenotypes are primarily due to altered brain circuits, motor neuron dysfunction or whole-body changes in metabolism. The body weight of both R6/2 and Q175 mice decreased over time, but we did not include this factor in our analyses: in Q175 mice, this weight loss began after the observed motor deficits, indicating that it was unlikely a direct cause of the observed locomotor changes [36,37].

The earliest changes in the whole brain, cortical or striatal volumes are detected at 2, 3 and 4 months of age in R6/2, Q175 HOM and Q175 HET mice, respectively [36,83]. That fact that our phenotypic scores allowed distinguishing significant phenotypes in the same models at 1, 1 and 3 month, respectively (Fig 4), indicates that the approach suggested by us can detect changes between mutant and WT mice even before neuronal loss starts. Power analysis revealed that the phenotypic score based on marker trajectories was more sensitive than that based on gait parameters in differentiating the transgenic mice from their WT littermates, indicating that the trajectory waveform patterns are more useful, albeit harder to interpret, than a set of parameters derived from the same raw data. Moreover, the trajectory data-based method can be considered more objective, because the parameters used were not subjectively selected. However, the kinematic parameter-based approach may be easier to interpret, as the individual parameters can be traced and emphasized in the final score. Another limitation in the trajectory data-based phenotype score is the loss of absolute speed information, because the method is based on interpolated trajectory data.

In summary, our hypothesis-free analysis of mouse movements both improved the statistical power of motor function assessment for therapeutic studies and detected early gait changes more sensitively than the conventional motor tests. The behavioral test described here uses data obtained during spontaneous locomotion and requires neither extensive training nor pharmacological induction. We show that robust phenotypes in mice can be observed early, within 1 month of age for R6/2 mice, shortening in vivo studies and potentially accelerating drug discovery. The data capture and analysis are semi-automated, maximizing study objectivity while minimizing chances for data handling mistakes. In addition, this behavioral test takes less than 20 seconds per mouse, enabling multiple cohorts to be tested by the same handler on the same day on the same apparatus, minimizing potential confounding effects [84]. Data analysis and generation of phenotypic scores is rapid, and interpretation of phenotypic changes is easy. Overall, our method improves the translational toolbox to evaluate HD therapeutics and is applicable to other movement disorders.

Acknowledgments

The authors thank Oxana Lavrova and Brenda Lager from CHDI for their help and support for this project, and Sanna Gustafsson from Charles River Discovery Services, Kuopio, Finland for her excellent technical assistance. The authors thank also Dr. Maksym Kopanitsa for his invaluable help.

References

  1. 1. Barclay OR. Some aspects of the mechanics of mammalian locomotion. J Exp Biol. 1953;30: 116–120.
  2. 2. Muybridge E. The attitudes of animals in motion. Palo Alto, CA; 1881.
  3. 3. Schwartz MH, Rozumalski A. The Gait Deviation Index: a new comprehensive index of gait pathology. Gait Posture. 2008;28: 351–357. pmid:18565753
  4. 4. Sparrow WA. The efficiency of skilled performance. J Mot Behav. 1983;15: 237–261. pmid:15151872
  5. 5. Croft JL, Schroeder RT, Bertram JEA. The landscape of movement control in locomotion: cost, strategy, and solution. Front Psychol. 2019;10: 716. pmid:31024381
  6. 6. Li KZH, Bherer L, Mirelman A, Maidan I, Hausdorff JM. Cognitive involvement in balance, gait and dual-tasking in aging: a focused review from a neuroscience of aging perspective. Front Neurol. 2018;9: 913. pmid:30425679
  7. 7. Mirelman A, Bonato P, Camicioli R, Ellis TD, Giladi N, Hamilton JL, et al. Gait impairments in Parkinson’s disease. Lancet Neurol. 2019;18: 697–708. pmid:30975519
  8. 8. Vuong K, Canning CG, Menant JC, Loy CT. Gait, balance, and falls in Huntington disease. Handb Clin Neurol. 2018;159: 251–260. pmid:30482318
  9. 9. Grande G, Triolo F, Nuara A, Welmer A-K, Fratiglioni L, Vetrano DL. Measuring gait speed to better identify prodromal dementia. Exp Gerontol. 2019;124: 110625. pmid:31173841
  10. 10. Takakusaki K. Functional neuroanatomy for posture and gait control. J Mov Disord. 2017;10: 1–17. pmid:28122432
  11. 11. Broom L, Ellison BA, Worley A, Wagenaar L, Sörberg E, Ashton C, et al. A translational approach to capture gait signatures of neurological disorders in mice and humans. Sci Rep. 2017;7: 3225. pmid:28607434
  12. 12. Heikkinen T, Bragge T, Turkia T, Nurmi A, Hodgson R. Psychoactive compounds induce distinctive fine motor deficits and gait disturbance in mice. Program No 311.27 / M10. San Diego, CA: Society for Neuroscience [2016]; Online.
  13. 13. Datson NA, González-Barriga A, Kourkouta E, Weij R, van de Giessen J, Mulders S, et al. The expanded CAG repeat in the huntingtin gene as target for therapeutic RNA modulation throughout the HD mouse brain. PLoS ONE. 2017;12: e0171127. pmid:28182673
  14. 14. Rostosky CM, Milosevic I. Gait analysis of age-dependent motor impairments in mice with neurodegeneration. J Vis Exp. 2018. pmid:29985360
  15. 15. Fiander MD, Chedrawe MA, Lamport A-C, Akay T, Robertson GS. Sagittal plane kinematic gait analysis in C57BL/6 mice subjected to MOG35-55 induced experimental autoimmune encephalomyelitis. J Vis Exp. 2017. pmid:29155752
  16. 16. Balkaya M, Kröber JM, Rex A, Endres M. Assessing post-stroke behavior in mouse models of focal ischemia. J Cereb Blood Flow Metab. 2013;33: 330–338. pmid:23232947
  17. 17. Zörner B, Filli L, Starkey ML, Gonzenbach R, Kasper H, Röthlisberger M, et al. Profiling locomotor recovery: comprehensive quantification of impairments after CNS damage in rodents. Nat Methods. 2010;7: 701–708. pmid:20836253
  18. 18. The Huntington’s Disease Collaborative Research Group. A novel gene containing a trinucleotide repeat that is expanded and unstable on Huntington’s disease chromosomes. Cell. 1993;72: 971–983. pmid:8458085
  19. 19. McColgan P, Tabrizi SJ. Huntington’s disease: a clinical review. Eur J Neurol. 2018;25: 24–34. pmid:28817209
  20. 20. Mestre TA, Forjaz MJ, Mahlknecht P, Cardoso F, Ferreira JJ, Reilmann R, et al. Rating Scales for motor symptoms and signs in Huntington’s disease: critique and recommendations. Mov Disord Clin Pract. 2018;5: 111–117. pmid:30363393
  21. 21. Huntington Study Group. Unified Huntington’s Disease Rating Scale: reliability and consistency. Mov Disord. 1996;11: 136–142. pmid:8684382
  22. 22. Biglan KM, Ross CA, Langbehn DR, Aylward EH, Stout JC, Queller S, et al. Motor abnormalities in premanifest persons with Huntington’s disease: the PREDICT-HD study. Mov Disord. 2009;24: 1763–1772. pmid:19562761
  23. 23. Rao AK, Mazzoni P, Wasserman P, Marder K. Longitudinal change in gait and motor function in pre-manifest Huntington’s disease. PLoS Curr. 2011;3: RRN1268. pmid:22008726
  24. 24. Delval A, Krystkowiak P, Delliaux M, Dujardin K, Blatt J-L, Destée A, et al. Role of attentional resources on gait performance in Huntington’s disease. Mov Disord. 2008;23: 684–689. pmid:18175353
  25. 25. Delval A, Krystkowiak P, Blatt J-L, Labyt E, Dujardin K, Destée A, et al. Role of hypokinesia and bradykinesia in gait disturbances in Huntington’s disease: a biomechanical study. J Neurol. 2006;253: 73–80. pmid:16096818
  26. 26. Pyo SJ, Kim H, Kim IS, Park Y-M, Kim M-J, Lee HM, et al. quantitative gait analysis in patients with Huntington’s disease. J Mov Disord. 2017;10: 140–144. pmid:28851209
  27. 27. Gaßner H, Jensen D, Marxreiter F, Kletsch A, Bohlen S, Schubert R, et al. Gait variability as digital biomarker of disease severity in Huntington’s disease. J Neurol. 2020;267: 1594–1601. pmid:32048014
  28. 28. Ramaswamy S, McBride JL, Kordower JH. Animal models of Huntington’s disease. ILAR J. 2007;48: 356–373. pmid:17712222
  29. 29. Menalled L, El-Khodor BF, Patry M, Suárez-Fariñas M, Orenstein SJ, Zahasky B, et al. Systematic behavioral evaluation of Huntington’s disease transgenic and knock-in mouse models. Neurobiol Dis. 2009;35: 319–336. pmid:19464370
  30. 30. Wheeler VC, Gutekunst C-A, Vrbanac V, Lebel L-A, Schilling G, Hersch S, et al. Early phenotypes that presage late-onset neurodegenerative disease allow testing of modifiers in Hdh CAG knock-in mice. Hum Mol Genet. 2002;11: 633–640. pmid:11912178
  31. 31. Menalled LB, Sison JD, Dragatsis I, Zeitlin S, Chesselet M-F. Time course of early motor and neuropathological anomalies in a knock-in mouse model of Huntington’s disease with 140 CAG repeats. J Comp Neurol. 2003;465: 11–26. pmid:12926013
  32. 32. Lin CH, Tallaksen-Greene S, Chien WM, Cearley JA, Jackson WS, Crouse AB, et al. Neurological abnormalities in a knock-in mouse model of Huntington’s disease. Hum Mol Genet. 2001;10: 137–144. pmid:11152661
  33. 33. Chen X, Wu J, Lvovskaya S, Herndon E, Supnet C, Bezprozvanny I. Dantrolene is neuroprotective in Huntington’s disease transgenic mouse model. Mol Neurodegener. 2011;6: 81. pmid:22118545
  34. 34. Mangiarini L, Sathasivam K, Seller M, Cozens B, Harper A, Hetherington C, et al. Exon 1 of the HD gene with an expanded CAG repeat is sufficient to cause a progressive neurological phenotype in transgenic mice. Cell. 1996;87: 493–506. pmid:8898202
  35. 35. Menalled LB, Kudwa AE, Miller S, Fitzpatrick J, Watson-Johnson J, Keating N, et al. Comprehensive behavioral and molecular characterization of a new knock-in mouse model of Huntington’s disease: zQ175. PLoS One. 2012;7: e49838. pmid:23284626
  36. 36. Heikkinen T, Lehtimäki K, Vartiainen N, Puoliväli J, Hendricks SJ, Glaser JR, et al. Characterization of neurophysiological and behavioral changes, MRI brain volumetry and 1H MRS in zQ175 knock-in mouse model of Huntington’s disease. PLoS One. 2012;7: e50717. pmid:23284644
  37. 37. Carter RJ, Lione LA, Humby T, Mangiarini L, Mahal A, Bates GP, et al. Characterization of progressive motor deficits in mice transgenic for the human Huntington’s disease mutation. J Neurosci. 1999;19: 3248–3257. pmid:10191337
  38. 38. Pallier PN, Drew CJG, Morton AJ. The detection and measurement of locomotor deficits in a transgenic mouse model of Huntington’s disease are task- and protocol-dependent: influence of non-motor factors on locomotor function. Brain Res Bull. 2009;78: 347–355. pmid:19010400
  39. 39. Federolf PA, Boyer KA, Andriacchi TP. Application of principal component analysis in clinical gait research: identification of systematic differences between healthy and medial knee-osteoarthritic gait. J Biomech. 2013;46: 2173–2178. pmid:23910389
  40. 40. Schutte LM, Narayanan U, Stout JL, Selber P, Gage JR, Schwartz MH. An index for quantifying deviations from normal gait. Gait Posture. 2000;11: 25–31. pmid:10664482
  41. 41. Jolliffe IT. Principal component analysis. 2nd ed. New York: Springer; 2002.
  42. 42. Abdi H, Williams LJ. Principal component analysis. WIREs Computational Statistics. 2010;2: 433–459.
  43. 43. Richman MB. Rotation of principal components. Journal of Climatology. 1986;6: 293–335.
  44. 44. Kaiser HF. The varimax criterion for analytic rotation in factor analysis. Psychometrika. 1958;23: 187–200.
  45. 45. Milovanović I, Popović DB. Principal component analysis of gait kinematics data in acute and chronic stroke patients. Comput Math Methods Med. 2012;2012: 649743. pmid:22400054
  46. 46. Gløersen Ø, Myklebust H, Hallén J, Federolf P. Technique analysis in elite athletes using principal component analysis. J Sports Sci. 2018;36: 229–237. pmid:28287028
  47. 47. Carlozzi NE, Miciura A, Migliore N, Dayalu P. Understanding the outcomes measures used in Huntington disease pharmacological trials: a systematic review. J Huntingtons Dis. 2014;3: 233–252. pmid:25300328
  48. 48. Abada Y-SK, Ellenbroek BA. Of rodents and men: understanding the emergence of motor and cognitive symptoms in Huntington disease. Behav Pharmacol. 2016;27: 403–414. pmid:26886208
  49. 49. William Yang X, Gray M. Mouse models for validating preclinical candidates for Huntington’s disease. In: Lo DC, Hughes RE, editors. Neurobiology of Huntington’s Disease: Applications to Drug Discovery. Boca Raton (FL): CRC Press/Taylor & Francis; 2011. Available: http://www.ncbi.nlm.nih.gov/books/NBK56001/
  50. 50. Kegelmeyer DA, Kloos AD, Fritz NE, Fiumedora MM, White SE, Kostyk SK. Impact of tetrabenazine on gait and functional mobility in individuals with Huntington’s disease. J Neurol Sci. 2014;347: 219–223. pmid:25456459
  51. 51. Ferrara JM, Mostile G, Hunter C, Adam OR, Jankovic J. Effect of tetrabenazine on motor function in patients with huntington disease. Neurol Ther. 2012;1: 5. pmid:26000211
  52. 52. Dai Y, Dudek NL, Li Q, Fowler SC, Muma NA. Striatal expression of a calmodulin fragment improved motor function, weight loss, and neuropathology in the R6/2 mouse model of Huntington’s disease. J Neurosci. 2009;29: 11550–11559. pmid:19759302
  53. 53. Alpaugh M, Galleguillos D, Forero J, Morales LC, Lackey SW, Kar P, et al. Disease-modifying effects of ganglioside GM1 in Huntington’s disease models. EMBO Mol Med. 2017;9: 1537–1557. pmid:28993428
  54. 54. Zimmermann T, Remmers F, Lutz B, Leschik J. ESC-derived BDNF-overexpressing neural progenitors differentially promote recovery in Huntington’s disease models by enhanced striatal differentiation. Stem Cell Reports. 2016;7: 693–706. pmid:27693427
  55. 55. Paldino E, Balducci C, La Vitola P, Artioli L, D’Angelo V, Giampà C, et al. Neuroprotective effects of doxycycline in the R6/2 mouse model of Huntington’s disease. Mol Neurobiol. 2020;57: 1889–1903. pmid:31879858
  56. 56. Deng Y, Wang H, Joni M, Sekhri R, Reiner A. Progression of basal ganglia pathology in heterozygous Q175 knock-in Huntington’s disease mice. Journal of Comparative Neurology. n/a. pmid:32869871
  57. 57. Abd-Elrahman KS, Hamilton A, Hutchinson SR, Liu F, Russell RC, Ferguson SSG. mGluR5 antagonism increases autophagy and prevents disease progression in the zQ175 mouse model of Huntington’s disease. Sci Signal. 2017;10. pmid:29259100
  58. 58. Reiner A, Lafferty DC, Wang HB, Del Mar N, Deng YP. The group 2 metabotropic glutamate receptor agonist LY379268 rescues neuronal, neurochemical and motor abnormalities in R6/2 Huntington’s disease mice. Neurobiol Dis. 2012;47: 75–91. pmid:22472187
  59. 59. Deluzio KJ, Astephen JL. Biomechanical features of gait waveform data associated with knee osteoarthritis: an application of principal component analysis. Gait Posture. 2007;25: 86–93. pmid:16567093
  60. 60. Alexandrov V, Brunner D, Menalled LB, Kudwa A, Watson-Johnson J, Mazzella M, et al. Large-scale phenome analysis defines a behavioral signature for Huntington’s disease genotype in mice. Nat Biotechnol. 2016;34: 838–844. pmid:27376585
  61. 61. Hsieh T-H, Chen J-JJ, Chen L-H, Chiang P-T, Lee H-Y. Time-course gait analysis of hemiparkinsonian rats following 6-hydroxydopamine lesion. Behav Brain Res. 2011;222: 1–9. pmid:21435355
  62. 62. Lee H-Y, Hsieh T-H, Liang J-I, Yeh M-L, Chen J-JJ. Quantitative video-based gait pattern analysis for hemiparkinsonian rats. Med Biol Eng Comput. 2012;50: 937–946. pmid:22707230
  63. 63. Wang XH, Lu G, Hu X, Tsang KS, Kwong WH, Wu FX, et al. Quantitative assessment of gait and neurochemical correlation in a classical murine model of Parkinson’s disease. BMC Neurosci. 2012;13: 142. pmid:23151254
  64. 64. Mancuso R, Oliván S, Osta R, Navarro X. Evolution of gait abnormalities in SOD1(G93A) transgenic mice. Brain Res. 2011;1406: 65–73. pmid:21733494
  65. 65. Vandeputte C, Taymans J-M, Casteels C, Coun F, Ni Y, Van Laere K, et al. Automated quantitative gait analysis in animal models of movement disorders. BMC Neurosci. 2010;11: 92. pmid:20691122
  66. 66. Hausdorff JM, Cudkowicz ME, Firtion R, Wei JY, Goldberger AL. Gait variability and basal ganglia disorders: stride-to-stride variations of gait cycle timing in Parkinson’s disease and Huntington’s disease. Mov Disord. 1998;13: 428–437. pmid:9613733
  67. 67. Bellardita C, Kiehn O. Phenotypic characterization of speed-associated gait changes in mice reveals modular organization of locomotor networks. Curr Biol. 2015;25: 1426–1436. pmid:25959968
  68. 68. Leblond H, L’Esperance M, Orsal D, Rossignol S. Treadmill locomotion in the intact and spinal mouse. J Neurosci. 2003;23: 11411–11419. pmid:14673005
  69. 69. Lazarus A, Stelmach GE. Interlimb coordination in Parkinson’s disease. Mov Disord. 1992;7: 159–170. pmid:1584238
  70. 70. Nanhoe-Mahabier W, Snijders AH, Delval A, Weerdesteyn V, Duysens J, Overeem S, et al. Walking patterns in Parkinson’s disease with and without freezing of gait. Neuroscience. 2011;182: 217–224. pmid:21382449
  71. 71. Tanahashi T, Yamamoto T, Endo T, Fujimura H, Yokoe M, Mochizuki H, et al. Noisy interlimb coordination can be a main cause of freezing of gait in patients with little to no parkinsonism. PLoS ONE. 2013;8: e84423. pmid:24391953
  72. 72. Garry MI, van Steenis RE, Summers JJ. Interlimb coordination following stroke. Hum Mov Sci. 2005;24: 849–864. pmid:16343673
  73. 73. Thibaudier Y, Hurteau M-F, Dambreville C, Chraibi A, Goetz L, Frigon A. Interlimb coordination during tied-belt and transverse split-belt locomotion before and after an incomplete spinal cord injury. J Neurotrauma. 2017;34: 1751–1765. pmid:27219842
  74. 74. Neumann M, Wang Y, Kim S, Hong SM, Jeng L, Bilgen M, et al. Assessing gait impairment following experimental traumatic brain injury in mice. J Neurosci Methods. 2009;176: 34–44. pmid:18805438
  75. 75. Hetze S, Römer C, Teufelhart C, Meisel A, Engel O. Gait analysis as a method for assessing neurological outcome in a mouse model of stroke. J Neurosci Methods. 2012;206: 7–14. pmid:22343052
  76. 76. Kale A, McCue S, Amende I, Hampton T. Hyperkinetic and rigid gait in R6/2 mice, a model of Huntington’s disease. B-C2079. Montreal, Canada; 2008. Available: https://www.can-acn.org/meeting2008/abstracts/352.htm
  77. 77. Rodríguez-Quiroga SA, Gonzalez-Morón D, Garretto N, Kauffman MA. Huntington’s disease masquerading as spinocerebellar ataxia. BMJ Case Rep. 2013;2013. pmid:23853009
  78. 78. Dong Y, Sun Y-M, Liu Z-J, Ni W, Shi S-S, Wu Z-Y. Chinese patients with Huntington’s disease initially presenting with spinocerebellar ataxia. Clin Genet. 2013;83: 380–383. pmid:22775398
  79. 79. Takeuchi E, Sato Y, Miura E, Yamaura H, Yuzaki M, Yanagihara D. Characteristics of gait ataxia in δ2 glutamate receptor mutant mice, ho15J. PLoS ONE. 2012;7: e47553. pmid:23077633
  80. 80. Shuvaev AN, Hosoi N, Sato Y, Yanagihara D, Hirai H. Progressive impairment of cerebellar mGluR signalling and its therapeutic potential for cerebellar ataxia in spinocerebellar ataxia type 1 model mice. J Physiol (Lond). 2017;595: 141–164. pmid:27440721
  81. 81. Grimbergen YAM, Knol MJ, Bloem BR, Kremer BPH, Roos RAC, Munneke M. Falls and gait disturbances in Huntington’s disease. Mov Disord. 2008;23: 970–976. pmid:18381643
  82. 82. Waln O, Jankovic J. An update on tardive dyskinesia: from phenomenology to treatment. Tremor Other Hyperkinet Mov (N Y). 2013;3. pmid:23858394
  83. 83. Stack EC, Kubilus JK, Smith K, Cormier K, Del Signore SJ, Guelin E, et al. Chronology of behavioral symptoms and neuropathological sequela in R6/2 Huntington’s disease transgenic mice. J Comp Neurol. 2005;490: 354–370. pmid:16127709
  84. 84. Gulinello M, Mitchell HA, Chang Q, Timothy O’Brien W, Zhou Z, Abel T, et al. Rigor and reproducibility in rodent behavioral research. Neurobiol Learn Mem. 2019;165: 106780. pmid:29307548