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

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 [3][4][5]. 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 [6][7][8][9][10]. 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 a1111111111 a1111111111 a1111111111 a1111111111 a1111111111 their gait has not been analyzed in detail previously. Judging from at best moderate gait phenotypes in other KI HD mouse models [29][30][31][32], 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.

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

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.
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 [34 -36]. 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 [39][40][41][42][43]. 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 (X T X/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 Step Width (hind, fore) The distance between forepaws or hindpaws when both are touching the ground during stance, perpendicular to midline Step Width Deviation (hind, fore) The deviation of step width between strides

Posture
Toe Clearance (hind, fore) Maximum distance of the paw from the ground during swing 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.

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).
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).
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-

PLOS ONE
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 linesmovement 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;  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 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.
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).

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 . 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

PLOS ONE
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.  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

PLOS ONE
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,54-58] and other [29-31,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 [61][62][63][64][65], 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 [69][70][71][72][73][74][75] in revealing decreased limb coordination in the mutants. In contrast, parameters related to the body posture-for example, the range of joint movementwere 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.