Abstract
Free full text
Motor Disorder in Huntington’s Disease Begins as a Dysfunction in Error Feedback Control
A steady progression of motor dysfunction takes place in Huntington’s disease1 (HD). The origin of this disturbance with relation to the motor control process is not understood. Here we studied reaching movements in asymptomatic HD gene-carriers (AGCs) and subjects with manifest HD. We found that movement jerkiness, which characterizes the smoothness and efficiency of motion, was a sensitive indicator of presymptomatic HD progression. A large fraction of AGCs displayed elevated jerk even when more than 7 years remained until predicted disease onset. Movement termination was disturbed much more than initiation yet was highly variable from trial to trial. Analysis of this variability revealed that the sensitivity of end-movement jerk to subtle, self-generated early-movement errors was greater in HD subjects than in controls. Additionally, we found that HD corrective responses to externally-generated force pulses were dramatically disturbed, indicating that HD subjects display aberrant responses to both external and self-generated errors. Because feedback corrections are driven by error and are delayed such that they predominantly effect movement termination, these findings suggest that a dysfunction in error correction characterizes the motor control deficit in early HD. This dysfunction may be observed years before clinical disease onset and it worsens with disease progression.
Huntington’s disease is an autosomal dominant inherited neurologic disorder caused by a glutamate repeat expansion in the IT15 gene2. Disease symptoms appear in the fourth or fifth decade of life1, and then progress steadily during the succeeding 10-20 years. Significant, progressive atrophy of the basal ganglia is evident in HD3, and may even be detectable before clinical onset4. It’s hallmark motor sign is chorea: the occurrence of rapid, irregular, and arrhythmic complex involuntary movements. Along with chorea, substantial impairment of voluntary movement also occurs1, and it may be of greater functional importance in the lives of patients5. Two control processes take part in the execution of arm movements. Feedforward control is the generation of motor commands based a priori on the desired action and an internal model of the system’s response6,7, while feedback control entails the ‘mid-flight’ correction of these commands based on errors detected during their execution8. To study these processes in HD, we used a high performance manipulandum9 to record visually guided reaching arm movements in patients with HD and presymptomatic individuals with the HD mutation up to 18 years from predicted disease onset.
Figure 1 displays hand paths of the two most regular and two least regular movements in each direction, for several subjects after 200 practice movements. The range of movement quality is generally much larger in HD patients than controls. Some movements made by symptomatic HD subjects appear normal while others are markedly irregular. Specifically, many HD movements have large changes in direction, smooth or abrupt, when approaching the target, while almost all movements fail to stop efficiently and smoothly. One way to characterize these irregularities is to quantify the smoothness of each movement. Smoothness can be defined as the lack of abrupt change, so a trajectory that minimizes abrupt changes in a variable will maximize its smoothness. Human reaching movements are of near maximal acceleration smoothness, quantified by minimal cumulative squared acceleration change10, also referred to as minimum total squared jerk.
Figure 2a shows that all HD patients and several AGCs made movements with above normal jerk. Notably, all HD subjects reduced their movement jerk with practice, although not nearly to control levels. The smoothness of movements increased over time even in subjects who were markedly more jerky than normal. In order to compare the initiation and completion of movement, we split each movement at the peak in the speed profile. This is roughly the point at which movement toward the target switches from acceleration to deceleration. Comparison of the total squared jerk before and after the peak in the speed profile (figure 2b) reveals that while fewer than half of the HD patients have high jerk during both parts of the movement, all have above normal post-peak jerk. Only 2 of 16 AGCs show above normal pre-peak jerk, but a majority have high post-peak jerk. To assess disease progression in AGCs, we used an estimate of disease onset age based on each subject’s parental onset age and glutamate repeat length4. The amount of post-peak jerk correlated significantly (r=-0.62) with estimated time to disease onset for AGCs while the pre-peak jerk did not (r=0.02). Post-peak jerk was above normal in 4 of 5 close to onset (<7 years) subjects and in 3 of 9 far from onset subjects. As groups, both HD patients and AGCs had significantly higher than normal post-peak jerk (p<10-6, p<0.00012) (see figure 2e). Moreover, AGCs who were close to predicted disease onset (< 7 years) had significantly higher jerk in the third set than far from onset (> 7 years) subjects (p<0.0062) who in turn had significantly higher jerk than controls (p<0.013).
Figure 2c shows that initial aiming is not dramatically disturbed in HD. Average directional aiming bias (see methods) is normal in 9 of 11 HD patients. Aiming variability is in the control range for a majority (6 of 11) of patients, although as a group, subjects with manifest HD display increased aiming variability. All AGCs have normal aiming biases of 3-7 degrees and none has above normal aiming variability.
These results suggest that HD movements often begin normally, but become jerky and irregular at some point during their course. Comparison of the average time course of raw squared jerk profiles between groups for movements in different speed ranges (figure 3), reveals a strikingly consistent pattern. At each speed range, the HD jerk profile closely matches the control profile during the beginning of movement, but begins to separate from it 200-300 msec after onset. Note that the end movement squared jerk is 10-30 times greater in HDs than controls.
Why might HD movements begin to become irregular 200-300ms into their course and not before? Corrective actions based on visual11,12 and proprioceptive13 information acquired during reaching movements begin to take place at about the time at which HD movements become irregular, so one possibility is that the system that generates these corrective actions is disturbed. Alternatively, there may be a disturbance in decelerating movement or another end-movement process independent of error correction.
To distinguish between these possiblities, we characterized the performance of the error feedback control system in HD. If the system that generates corrective actions during movement (the feedback controller) is affected in HD, then a strong relationship should exist between the magnitude of error early in the movement and size of disturbance later. Error, the difference between desired system state and actual state, is the primary input to a feedback controller. If the error during a movement is small, then large corrective actions are not necessary and dysfunction in feedback control might have only minimal effects. However, if errors are large, substantial corrective actions are required, and dysfunctional feedback control may cause significant problems with movement.
While we can measure actual movements quite accurately, we have no reliable way of estimating the desired motion state within a single movement. This makes direct estimate of error infeasible. If errors are symmetrically distributed, then the average movement profile should approximate the average motion plan. However, subjects may plan to move a bit faster on one trial then slower on the next, and these variations in the plan may be as large as the typical errors between plan and action. Although the motor plan might change across trials, a few of its properties are invariant. Velocity profiles of reaching movements are symmetric and unimodal: the peak in the velocity profile occurs very close to the movement’s midway point, regardless of amplitude, direction, or speed14. Therefore, the distance between the midway point and the position where the peak speed occurs can be used as an indicator of error early in the movement.
The relationship between error early in movement and jerk late in movement is shown in figure 4. The probability distributions of peak speed positions (PSPs) do not differ dramatically between groups (figure 4a,c,e,g,i). This implies that amount of error early in the movement is comparable between groups. However, the reaction to this error is very different between groups. Figure 4b,d,f,h,j show that the jerk that occurs after the hand reaches the PSP varies considerably with PSP in all subject groups. When the PSP is far from the midway point, indicating large error early in the movement, post-peak jerk is high. However, the sensitivity of post-peak jerk to PSP is much greater in HD patients and AGCs than controls. Figure 4h shows that the difference in mean post-peak jerk between controls and AGCs is much greater at large PSP displacements than at small displacements. Subjects with manifest HD (p=0.006) and AGCs (p=0.04) display significantly higher sensitivity in their response of end-movement jerk to early movement error than do controls, as measured by slope in this relationship. The sensitivity of post-peak jerk to PSP displacement indicates the degree to which end-movement smoothness depends on early-movement error. The increase in this sensitivity in HD subjects with and without clinically manifest symptoms suggests that an error dependent control process is disturbed early in the disease course and further deteriorates with disease progression.
Previous reports have shown that cortical sensorimotor pathways are affected in patients with manifest HD. While short-loop reflexes that are mediated by spinal mechanisms are normal, long-loop reflexes, which involve the transfer of proprioceptive information through cortical pathways, are reduced or absent in HD15,16. Cortical responses to peripheral nerve stimulation, as measured by somatosensory evoked potentials (SEPs), are also reduced in HD15,17,18, suggesting that diminished cortical sensory input is responsible for the long-loop reflex reduction16. Additionally, a strong correlation exists between SEP deterioration and striatal glucose metabolism early in the course of HD17, hinting at a possible involvement of the striatum in the pathology of the SEPs and long-loop reflexes. Since cortical sensorimotor pathways play a large role in mediating error correction during voluntary movement, disruption of these pathways could lead to dysfunctional feedback control during movement.
Our analysis of the response to internal, self-generated errors suggested that error correction in general might be disturbed in HD. To test this hypothesis we externally imposed errors upon movements of these subjects via an occasional force pulse (70ms) shortly after movement initiation. These pulses were given randomly on a minority of trials, and could be in 8 different directions and of 3 different magnitudes: 6, 12, and 18 Newtons peak force. In addition to subjects with HD, we studied the performance of individuals with cerebellar deficits on this task. Hand paths of two symptomatic HD, one AGC, one control, and three cerebellar subjects for movements perturbed with 12 and 18 Newton force pulses are shown in figure 5a. These perturbations are substantial: trajectories can be altered by up to 5 cm (movement distance is 10cm). Control subjects correct the perturbations smoothly and efficiently, but extremely large successive overshoots are seen during certain HD movements. Not all HD movement corrections appear irregular, but some are so dramatically disturbed that they bear little qualitative similarity to any movements that we recorded from controls.
Several movements made by cerebellar subjects appear to have irregular and inefficient corrections, but close inspection of many of these trajectories reveals a striking similarity between the pattern of overshoots in these movements and the overshoots present in unperturbed movements in the same direction. This suggests that the irregularity that exists in these perturbed movements has a stronger relationship to the direction of movement than to the presence of perturbation. Note the movements toward 3:00, 4:30, and 7:30 for subjects 5 and 7.
HD subjects appear to have dysfunctional reactions to movement errors caused by external perturbations. During the post perturbation movement segment (the corrective period), jerk and path length are relatively increased over unperturbed movements significantly more for HD subjects than controls (see figure 5b) for both 12 and 18 N perturbations. The error-correction performance of AGCs generally falls between that of controls and subjects with manifest HD. Cerebellar subjects, like symptomatic HDs, had worse unperturbed movement performance than controls, but the decrement in their performance when perturbations were given was generally more like that of controls than HD subjects. This suggests that subjects with HD generally have greater deficits in error feedback control than do cerebellar patients.
The existence of performance deficits on tests of behavioral function before the clinical onset of HD is controversial. To date, assessments of motor, cognitive, or psychiatric function in HD have revealed only subtle deficits in presymptomatic subject groups19,20,21 if any22,23. When changes have been detected, they were not sufficiently specific to permit discrimination between mutation positive and mutation negative individuals, or even reliable identification of people with early stage manifest HD. In contrast, low glucose metabolism in the basal ganglia has been reported in about two thirds of asymptomatic at risk individuals tested17,24,25 and basal ganglia volumes may be reduced years before clinical onset4, suggesting that brain pathology precedes manifestation of the behavioral dysfunction previously studied. However, end-movement jerk appears to be a sensitive indicator of HD progression, suggesting the existence of a direct behavioral correlate to the early brain pathology of HD.
Our results, which suggest relatively spared feedforward control but dysfunctional feedback control in HD, may help explain the pattern of motor learning deficits reported. Learning on rotary pursuit tasks, which involve long continuous movements largely under closed-loop feedback control, is impaired in HD26,27. In contrast, learning to mirror trace polygons, which involves multiple short, discrete movements that rely more heavily on open-loop feedforward control, is intact in HD.
Real time error feedback control presents a formidable challenge to the CNS because of the large sensorimotor delays that occur28. One way to improve the performance of such a time-delayed system is to include a component which can effectively predict away the delay and allow efficient error correction. This component could combine delayed sensory information with a knowledge of system’s physical dynamics and recent motor output to predict future motion state. This sort of predictor has been referred to as a forward model of system dynamics8,29. Neurons in the basal ganglia have been shown to predict reward30, hinting that predictive capacity may be a general feature of some basal ganglia structures. The main output of the basal ganglia modulates the action of the thalamus, which relays sensory information to the cortex. This information stream is likely to participate in error feedback control, thus the error correction process may be modulated by the basal ganglia, although the mechanism through which this occurs, including whether predictive capacity in basal ganglia structures is involved, is not yet understood.
Methods
Subjects
11 patients positive for the IT-15 mutation and symptomatic with Huntington’s Disease, 16 mutation-positive presymptomatic subjects, 3 mutation-negative subjects who had a parent with HD, and 12 other age-matched controls participated in the first experiment. 5 symptomatic and 9 presymptomatic HD subjects, 8 age-matched controls, and 6 subjects with cerebellar lesions participated in the second experiment. All subjects used their dominant hand, and all but one presymptomatic in the first experiment were right handed. The direct gene test for IT-15 mutation was conducted at the Johns Hopkins Huntington’s Disease Project. The number of CAG (glutamate) trinucleotide repeats was determined, and subjects with >37 repeats were called mutation-positive. Subjects with <34 were called mutation-negative. All cerebellar subjects that we studied had been diagnosed clinically with cerebellar dysfunction, and all had lesions localized to the cerebellum on MRI. 4 patients had generalized cerebellar atrophy, and 2 had suffered strokes of the right posterior inferior cerebellar artery (PICA). One of these patients also had a left PICA and a right superior cerebellar artery stroke.
Task
Subjects made quick reaching movements to targets spaced 10cm away while grasping a lightweight two-joint manipulandum. The 1cm square targets and a small cursor indicating the subject’s hand position were displayed on a computer monitor in front of the subject6. We define the end of the movement as the end of the first time interval after movement onset when hand velocity always remained below a threshold of 0.03m/s for 200ms. During the first experiment, the training took place in two sessions and within sessions it was subdivided into 100-movement sets. In the first session, subjects completed three sets of training. After a 3 to 4 hour break, they began the second session, which consisted of a single set. The robot arm remained passive for all movements in both sets. During the second experiment the robot produced a 70 ms bell shaped force pulse on a minority of randomly pre-selected trials (with a probability of 1 in 4). The force pulse could be in any one of 8 directions and of magnitude 6, 12, or 18N. One force pulse in each direction and of each magnitude was given for each direction of movement. During both experiments we used a sling suspended from the ceiling to support the subject’s upper arm in the horizontal plane. This helped regularize the subjects’ arm position and minimize the effort required to support the arm against gravity.
Analysis
We characterized aiming direction by defining angular aiming error as difference between target direction and the direction of travel up to peak speed point. Note that aiming error for a given movement can be positive or negative. Aiming bias for each subject was defined as the average of magnitude of the average aiming error in each direction. Similarly, aiming variability for each subject was defined as the average of the standard deviation of aiming error in each direction
Jerk is defined as the rate of change of acceleration with respect to time. In order to minimize the effect of discretization noise on the differentiation of the velocity signal, jerk was estimated by applying a Savitsky-Golay filter. The minimum total squared jerk (TSJ) required for a point-to-point movement is proportional to the fifth power of the movement speed and inversely proportional to the third power of movement excursion. Because of the strength of the relationship between these variables and TSJ, it is critical to appropriately normalize movement speed and excursion when comparing the amount of jerk between two movements. To account for the effect of speed and excursion on a movement’s typical jerk, we normalized the TSJ measured for each movement by the average TSJ for movements made after practice by a large separate group of controls (n=35), at the appropriate speed and excursion. Similarly, we normalized the TSJ in the pre-peak and post-peak movement segments by the average values of these quantities for control movements of the same peak speed and excursion, and we refer to these quantities as pre-peak jerk and post-peak jerk.
Since the trajectories of perturbed movements were often quite irregular, and the peak speed was often strongly influenced by the presence and direction of the perturbing force pulse, normalization of jerk by peak speed was no longer appropriate for these movements. Instead, we normalized the jerk measured after perturbation offset for a given movement by the minimum possible jerk required to make the movement state transition characteristic of that movement segment. We refer to this quantity as post-perturbation jerk. Similarly, we define post-perturbation path length as the path length between the point of perturbation offset and the end point of movement divided by the straight-line distance between these two points. We compare the values of these quantities for movements during which perturbations were given to unperturbed movements. For unperturbed movements, we define the time point of “perturbation offset” as that time when perturbation offset would have occurred had a perturbation been given.
Suppliment to Methods (internet)
Position and velocity of the joints of the robot arm were recorded at 100 Hz from absolute and relative joint position encoders with resolution of 5.5*10-3 degrees and 8.0*10-4 degrees, respectively. This produced estimates of hand position and velocity in Cartesian coordinates with accuracy greater than 0.1mm and 1.3mm/sec, respectively.
Jerk was estimated by applying a fourth order Savitsky-Golay filter on a 250ms window of velocity data. This filter is equivalent to taking the second derivative at the window’s center of the continuous least squares best fit fourth order polynomial. This fourth order polynomial fit is a low-pass filter with cutoff frequency of 6.83Hz. Power spectra of mean subtracted velocity profiles of very fast 10cm reaching movements show that 99.9% of the power is below 6Hz.
We assessed motion state transition efficiency using cumulative squared jerk in order help to characterize the efficiency of recovery after perturbation offset. To accomplish this, we compared the amount of jerk that occurred between two different motion states within the same movement, with the jerk that would occur for a maximally smooth transition between those two states in the elapsed time. The minimum jerk trajectory between two motion states (state = [position, velocity, acceleration]) is given by a fifth order polynomial in time:
Where position, velocity, and acceleration are represented by x(t),v(t), and a(t) respectively, t is time and tf is the final time. Ck’s are parameters that depend on the boundary motion states and on the time between them, tf. They can be found by solving the following linear system of equations.
Once the coefficients are determined the cumulative squared jerk can be computed by simply integrating the squared jerk profile.
Acknowledgements
This work was greatly enriched by interactions with K. Thoroughman and N. Bhushan. The work was supported by grants from the Whitaker Foundation (to R.S.) and the National Institutes of Health (to R.S. and J.B.), and a pre-doctoral fellowship from the National Institute of General Medical Sciences to M.S. R.S. and M.S. conceived of the experiments, M.S. collected the data, and designed and carried out the data analysis guided by interaction with R.S.,J.B.,N.B., and K.T., M.S. and R.S. wrote the manuscript. J.B. and the Hopkins HD Project arranged patient visits and clinical assessment of HD patients.
References
Full text links
Read article at publisher's site: https://doi.org/10.1038/35000576
Read article for free, from open access legal sources, via Unpaywall: http://www.seas.harvard.edu/motorlab/Reprints/nature00.pdf
Citations & impact
Impact metrics
Citations of article over time
Smart citations by scite.ai
Explore citation contexts and check if this article has been
supported or disputed.
https://scite.ai/reports/10.1038/35000576
Article citations
The gut microbiota modulate locomotion via vagus-dependent glucagon-like peptide-1 signaling.
NPJ Biofilms Microbiomes, 10(1):2, 16 Jan 2024
Cited by: 4 articles | PMID: 38228675 | PMCID: PMC10791613
Rapamycin reduces neuronal mutant huntingtin aggregation and ameliorates locomotor performance in Drosophila.
Front Aging Neurosci, 15:1223911, 26 Sep 2023
Cited by: 5 articles | PMID: 37823007 | PMCID: PMC10562706
Comparison of Various Smoothness Metrics for Upper Limb Movements in Middle-Aged Healthy Subjects.
Sensors (Basel), 23(3):1158, 19 Jan 2023
Cited by: 2 articles | PMID: 36772197 | PMCID: PMC9919347
A single bout of moderate-intensity aerobic exercise improves motor learning in premanifest and early Huntington's disease.
Front Psychol, 14:1089333, 08 Mar 2023
Cited by: 3 articles | PMID: 36968757 | PMCID: PMC10032374
The Impact of Dysregulated microRNA Biogenesis Machinery and microRNA Sorting on Neurodegenerative Diseases.
Int J Mol Sci, 24(4):3443, 08 Feb 2023
Cited by: 5 articles | PMID: 36834853 | PMCID: PMC9959302
Review Free full text in Europe PMC
Go to all (190) article citations
Similar Articles
To arrive at the top five similar articles we use a word-weighted algorithm to compare words from the Title and Abstract of each citation.
Intact ability to learn internal models of arm dynamics in Huntington's disease but not cerebellar degeneration.
J Neurophysiol, 93(5):2809-2821, 29 Dec 2004
Cited by: 295 articles | PMID: 15625094
Quantitative Electroencephalographic Analysis Provides an Early-Stage Indicator of Disease Onset and Progression in the zQ175 Knock-In Mouse Model of Huntington's Disease.
Sleep, 39(2):379-391, 01 Feb 2016
Cited by: 27 articles | PMID: 26446107 | PMCID: PMC4712407
Effects of chronic stress on the onset and progression of Huntington's disease in transgenic mice.
Neurobiol Dis, 71:81-94, 01 Aug 2014
Cited by: 19 articles | PMID: 25088714
Huntington's disease - clinical signs, symptoms, presymptomatic diagnosis, and diagnosis.
Handb Clin Neurol, 100:3-13, 01 Jan 2011
Cited by: 17 articles | PMID: 21496568
Review
Funding
Funders who supported this work.
NIGMS NIH HHS (1)
Grant ID: T32 GM007057
NINDS NIH HHS (2)
Grant ID: R01 NS037422
Grant ID: R01 NS037422-02