- 1Division of Biokinesiology and Physical Therapy, University of Southern California, Los Angeles, CA, United States
- 2Department of Biomedical Engineering, University of Southern California, Los Angeles, CA, United States
Parkinson's disease (PD) is typically diagnosed and evaluated on the basis of overt motor dysfunction, however, subtle changes in the frequency spectrum of neural drive to muscles have been reported as well. During dynamic actions, coactive muscles of healthy adults often share a common source of 6–15 Hz (alpha-band) neural drive, creating synchronous alpha-band activity in their EMG signals. Individuals with PD commonly exhibit kinetic action tremor at similar frequencies, but the potential relationship between the intermuscular alpha-band neural drive seen in healthy adults and the action tremor associated with PD is not well-understood. A close relationship is most tenable during voluntary dynamic tasks where alpha-band neural drive is strongest in healthy adults, and where neural circuits affected by PD are most engaged. In this study, we characterized the frequency spectrum of EMG synchronization (intermuscular coherence) in 16 participants with PD and 15 age-matched controls during two dynamic motor tasks: (1) rotation of a dial between the thumb and index finger, and (2) dynamic scaling of isometric precision pinch force. These tasks produce different profiles of coherence between the first dorsal interosseous and abductor pollicis brevis muscles. We sought to determine if alpha-band intermuscular coherence would be amplified in participants with PD relative to controls, if such differences would be task-specific, and if they would correlate with symptom severity. We found that relative to controls, the PD group displayed amplified, but similarly task-dependent, coherence in the alpha-band. The magnitude of coherence during the rotation task correlated with overall symptom severity as per the UPDRS rating scale. Finally, we explored the potential for our coherence measures, with no additional information, to discriminate individuals with PD from controls. The area under the Receiver Operating Characteristic curve (AUC) indicated a clear separation between groups (AUC = 0.96), even though participants with PD were on their typical medication and displayed only mild-moderate symptoms. We conclude that a task-dependent, intermuscular neural drive within the alpha-band is amplified in PD. Its quantification via intermuscular coherence analysis may provide a useful tool for detecting the presence of PD, or assessing its progression.
Introduction
The evaluation of Parkinson's disease (PD) is currently dominated by subjective clinical ratings of symptom severity, such as the Unified Parkinson's Disease Rating Scale (UPDRS). The coarse nature of this examination, along with the diversity of possible symptoms, has driven a search for more direct, quantitative measures of neural dysfunction which can objectively assess the early presence and progression of the disease.
The frequency spectrum of neural activity within the motor system is altered in PD. Neural oscillations are ubiquitous in the healthy motor system, but PD is characterized by a particularly large variety of abnormal oscillations. Tremor at rest (rest tremor) typically has a frequency of ~3–6 Hz, while “action tremor” extends from ~6–15 Hz and occurs during voluntary static (postural action tremor) or dynamic (kinetic action tremor) muscle activation (1–6). Further, PD is associated with increased corticomuscular drive in the range of 15–30 Hz during static contractions (7–11), as well as reduced neural drive to muscles in the 30–50 Hz range in unmedicated patients (12).
During motor behavior, coactive muscles often share a portion of their neural drive which synchronizes their activities at different frequencies. This entrainment of muscle activity by a common source of oscillatory neural drive can be quantified in the frequency domain by calculating the coherence between their EMG signals (13–20). Intermuscular coherence associated with postural action tremor (21) or rest tremor (22) in PD suggests that the distribution of tremulous neural drive across muscles could be an important feature of the disease, and one which cannot be assessed reliably within the context of typical clinical evaluations. Clinical evaluations, and indeed most scientific studies, have focused on the visible/overt forms of tremor, but it has been known for decades that in PD, dynamic voluntary activity consistently evokes kinetic action tremor in the alpha-band (~6–15 Hz) which simultaneously affects multiple muscles and is observable in EMG even when no visible tremor is apparent (1).
In fact, this kinetic action tremor is the strongest and most consistently-evoked form of non-overt tremor in PD (1, 4, 23–25), and can be found in most individuals with PD. While kinetic action tremor is often described as an amplification of the ~6–15 Hz physiological tremor seen in healthy adults due to its similarity in frequency and the fact that this frequency doesn't change with loading (4, 26), its underlying neurophysiology in PD is not well-understood, especially since ~6–15 Hz neural drive to muscles can come from a variety of different sources (27–36). Also, compared with other manifestations of PD, kinetic action tremor has received relatively little attention.
Most previous studies of action tremor in PD have focused on forces, motions, or individual muscles rather than the coherence of tremor-generating drive across muscles. If the relevant descending drive is fundamentally intermuscular, then action tremor should not only depend on action, but also on the dynamic coordination among muscles required by a given task, as this is known to influence the strength of ~6–15 Hz intermuscular neural drive (16, 19, 32, 37, 38). While elevated intermuscular coherence between anatomically-synergistic muscles has been found during static voluntary tasks in PD (7, 11), these findings may not extend to dynamic actions, or to functionally-different muscles whose coordination can change depending on task.
The aim of this study was therefore to characterize intermuscular coherence during voluntary dynamic tasks in participants with PD, compared to age-matched controls. Specifically, we tested two tasks: (1) rotation of a dial between the thumb and index finger, and (2) dynamic scaling of isometric precision pinch force. These tasks were chosen because they evoke different levels of coherence (especially between 6 and 15 Hz) between the first dorsal interosseous (FDI) and abductor pollicis brevis (APB) muscles (19), and because dynamic multifinger manipulation with the fingertips evokes strong functional coupling among the fingers (39).
Our primary hypotheses were that (1) intermuscular coherence would be larger in PD, especially at ~10 Hz during both dynamic tasks, (2) the task-dependent modulation of intermuscular coherence seen in controls would be preserved in PD, and (3) the amplification of alpha-band coherence would correlate with clinical severity, since dynamic tasks should preferentially utilize neural pathways known to deliver alpha-band drive to muscles, and which are known to be disrupted in PD, such as the cerebello-thalamo-cortical circuit (30, 40–42). Given that intermuscular coherence analysis has been suggested as a potential avenue for biomarker development (11), and that we lack simple, cost-effective methods for detecting the presence of PD and the severity of neural damage, a secondary aim was to determine the ability of coherence measures to discriminate between PD and control groups, as this would justify future efforts to develop clinically-applicable metrics using EMG.
Methods
Participants
We recruited 16 individuals with mild-moderate severity Parkinson's disease (age: 62.3 +/– 8.6, 8 male) and 15 control participants (age: 60.5 +/– 10.3, 9 male). The details of the patient population are shown in Table 1. All participants with PD were on their normal medication, and in the ON state at the time of testing. All were diagnosed with idiopathic Parkinson's disease, and all but one were on dopaminergic medication. All participants understood the task and scored >23 on the mini mental status exam. Our sample size is intended to be sufficient for detection of large differences and correlations, and is in line with similar recent literature [e.g., Flood et al., (11)] where such effects were found.
All participants gave written informed consent prior to participation and all procedures were approved by the University of Southern California Institutional Review Board.
Experimental Setup
We asked participants to pinch or rotate a custom-made dial (diameter: 3 cm) between the thumb and index finger, as described in Laine and Valero-Cuevas (19) (see Figure 1). Both tasks generate a physiological tremor in the muscles of the thumb and index finger, and coherence between their EMG signals varies across tasks even without related changes in the shape of their power spectra. Therefore, these tasks alter the extent to which alpha-band neural drive is shared among muscles rather than simply altering its amplitude. Briefly, the dial held a potentiometer to track rotation angle and a miniature load cell (ELB4-10, Measurement Specialties, Hampton, VA) under the index finger to measure pinch force. Surface EMG sensors (Biometrics, Newport, UK) were placed over the first dorsal interosseous (FDI) and abductor pollicis brevis (APB) muscles of each hand. All signals were acquired at 1,000 Hz using a Biometrics LTD DataLINK system (Biometrics, Newport, UK). Visual feedback of rotation or force was provided using custom software designed in MATLAB (The Math Works, Natick, MA). We instructed participants to prioritize the production of a smooth force or rotation effort, guided/paced by the target sinusoid. This instruction was intended to reduce the likelihood that participants with involuntary tremors would focus on counteracting them rather than executing the prescribed slow voluntary action. Each participant completed four, 3-min tasks with each hand, (2 trials for each of 2 tasks, described below). Practice trials were given prior to recordings, and breaks between each task were given to prevent fatigue. The order of tasks and hands was randomized for each participant, and subjects did not find these simple tasks fatiguing.
Figure 1. Behavioral Task: Participants grasped a small dial between the thumb and index finger and either (1) slowly scaled their pinch force using visual feedback to track a 0.1 Hz sinusoid spanning 1–3 N, or (2) rotated the dial back and forth over 45 degrees at 0.1 Hz using visual feedback of dial rotation. A miniature load cell at the index finger pad measured pinch force, and surface EMG was recorded over the first dorsal interosseous (FDI) and abductor pollicis brevis (APB) muscles. The EMG signals were used to calculate intermuscular coherence. The bottom panels show coherence profiles for each task from three example individuals in the control group (top) and in those with Parkinson's disease (PD, bottom). The displayed identification numbers correspond to the participant descriptions in Table 2.
Task 1: Dynamic modulation of isometric pinch force. With each hand separately, participants pinched the dial and slowly varied their pinch force between 1 and 3 N by tracking a sinusoidal target displayed on screen. The vertical height of the cursor was controlled by pinch force and while the horizontal position moved left-to-right across the screen automatically as a function of time, taking 30 s before looping back to the left. The sinusoidal target had a peak-to-peak period of 10 s, such that the frequency of force modulation was 0.1 Hz. We chose this frequency for ease of tracking and used it for all participants to avoid potential effects of movement speed on our EMG measurements. Practice trials were given to familiarize each participant with the task prior to recordings.
Task 2: Dial rotation. In this task, subjects rotated the dial back and forth +/– 22.5 degrees with each hand separately. Visual feedback was provided as before, but with the vertical position of the cursor controlled by rotation angle. To ensure that pinch force remained similar across tasks, the cursor color changed to indicate if pinch force exited the 1–3 N range during the task. All participants were able to maintain pinch force within this range and made few errors after initial practice.
EMG signals were high-pass filtered using a zero-phase 4th order Butterworth filter with a cutoff at 250 Hz, then rectified and normalized to unit variance as in our previous study (19). This follows the general recommendations for accentuating the timing and density of motor unit action potentials within the surface EMG signal (43–47). It should be noted that the extreme high-pass filtering may be a precaution more than a universal necessity, and that surface EMG signals can be expected to contain some degree of noise or amplitude cancellation which can distort coherence measures when assessing low frequency neural drive (e.g., <5 Hz) or using high contraction levels (48), neither of which are a concern in the present study. The two trials for each task for each hand were concatenated, yielding a total of 6 min of data per hand. The first few (~5) seconds of each trial were trimmed manually to ensure that stable tracking had been obtained prior to analysis.
Coherence Analysis
Coherence between EMG signals describes the frequency content of their synchronized activity. Coherence between the EMG signals of the FDI and APB muscles was calculated using the “mscohere” function in MATLAB, specifying segment sizes of 2 s, tapered with a Hann window, and overlapped by 0.5 s. Prior to statistical comparisons, the raw coherence values (C) were first converted to Fisher's Z values using the formula Fz = atanh [sqrt(C)]. Then, for better comparison with previous work, and to provide a more standard index of statistical strength, we converted the Fz values to standard Z-scores using the formula Z = Fz / √ (1/ 2L)—bias. In this formula, L is the degrees of freedom derived from the number of segments used in the coherence calculation (49–51), and the bias was calculated as the mean uncorrected Z-score between 100 and 500 Hz, since this frequency range contains no physiological coherence (19, 52). Coherence profiles for three individuals from each group are shown in the bottom panels of Figure 1.
Statistical Testing
Statistical evaluation of coherence is often simplified by binning the frequency spectrum into a few common bands of interest. The exact boundaries for each band can vary across studies, and such boundaries may change when evaluating pathology. Therefore, to address our main hypotheses, we allowed relevant frequency bands to be defined from the data itself, using a non-parametric version of statistical parametric mapping (SPM), as described previously (53, 54). This is a random-permutation test that assigns p-values to regions of interest (corrected for multiple comparisons) within a “map” of cross-group differences calculated over space, time, or (in our case) frequencies (55). To map group differences in coherence across frequencies, we used an effect size measure, Cohen's D, as our initial statistic. This was then smoothed over frequencies using a gaussian window spanning 4 Hz. We defined regions of interest, or “clusters,” as any group of consecutive frequencies (min width = 3 Hz) exceeding a threshold. The threshold can be set arbitrarily but we automated this by using the 95% confidence interval for our Cohen's D values, calculated with respect to their mean and standard deviation from 100 to 500 Hz, where no true group level differences in coherence should exist. The above-threshold area of each cluster was tested for statistical significance based on a 10,000 iteration random permutation test, as described previously (55).
Using the above procedure, we identified differences in intermuscular coherence between controls and participants with PD for (1) the rotation task, (2) the scaling task, and (3) the difference between rotation and scaling (rotation-scaling, per hand). The latter addresses whether any effects of PD on intermuscular coherence are task-specific. For each, we also created box and whisker plots to visualize how the average coherence within each statistically significant frequency band varied across individuals. To better assess variability across individuals rather than hands, data from both hands were averaged per individual. A Cohen's D effect size as well as a p-value (derived from a random-permutation test) were also calculated for this binned/averaged data. Since it is possible that PD may have stronger effects on one hand, we confirmed the appropriateness of averaging across hands by calculated an absolute laterality index for each individual (the absolute difference in coherence between hands divided by their sum), and comparing across groups, again using a random-permutation test.
Correlation Analysis
The intermuscular coherence values above were tested for correlation with the total UPDRS score, as well as with its subsection II (activities of daily living) and subsection III (motor evaluation). To be conservative, this was conducted using a non-parametric Spearman's rank correlation. The correlation coefficients and associated p-values were obtained using the “corr” function in MATLAB. This analysis tested for a non-zero correlation between coherence and symptom measures. With so few participants, an exact magnitude of correlation cannot be determined with high precision. We therefore calculated a 95% confidence interval around each correlation measure using a 10,000 iteration bootstrap procedure.
Discriminability
We evaluated the extent to which the data above could classify a participant as being a member of the PD or control group. To do this, we constructed Receiver Operating Characteristic curves (ROC curves) for each task [for a brief overview, see Eng (56)]. Each point on an ROC curve describes the fraction of patients who could be correctly identified (true positive rate, y-axis) using a particular threshold value for discrimination, while at the same time misclassifying some proportion of the controls (false positive rate, x-axis). Defining these proportions for every possible threshold produces the ROC curve. The area under the curve (AUC) is 1 for perfect discrimination and 0.5 for random chance. Our SPM analysis identified that the greatest difference of coherence between PD and control groups occurred in the alpha-band for both tasks (5.8–16.6 Hz for the rotation task and 4.4–15.6 for the scaling task). Using this information, we extracted a single measurement from each coherence spectrum, herein referred to as the alpha-ratio, by calculating the fraction of total intermuscular coherence (between 2 and 80 Hz) falling within the alpha-band (defined per task, as described above). This ratio-based normalization strategy has been used previously for reporting tremor measures (23, 57), and in our case, it reduces inter-subject variability from sources that could influence all frequencies at once (e.g., noise or cross-talk) while emphasizing the overall shape of the coherence spectrum. ROC curves were first constructed for each task. Then, for a final combined analysis, we created an ROC curve after averaging all 4 alpha-ratios obtained for each participant (2 hands × 2 tasks).
Task Performance
All participants could execute the task with both hands. Because systematic differences in the overall frequency/magnitude of error corrections could influence alpha-band drive to individual muscles (58, 59) we used a 10,000 iteration permutation test to determine if fluctuations in force/rotation about the target sinusoid differed significantly between groups. Here, performance was quantified as the standard deviation of pinch force or rotation angle after filtering out the 0.1 Hz modulation associated with the voluntary task of tracking the target sinusoid. Task performance was also tested for correlation with alpha-band drive using Spearman's rank correlation.
Results
For both tasks, alpha-band coherence between the FDI and APB muscles differed between the PD and control group. For the rotation task, coherence differed within the frequency range of 5.8 to 16.6 Hz, with p < 0.001 (see Figures 2A,G). The scaling task showed significant difference between groups from 4.4 to 15.6 Hz, with p < 0.001 (Figures 2B,H). The difference between the two tasks (rotation—scaling) showed a range of interest between 7.8 and 12.2 Hz (Figures 2C,I), but was not statistically significant (p = 0.58). Intersubject variability was generally high (Figures 2D–F).
Figure 2. Parkinson's disease (PD) affects intermuscular coherence in the alpha-band (6–15 Hz). The three columns represent analyses of group differences in intermuscular coherence (Control, black, 15 participants; Parkinson's, red, 16 participants) during the rotation task (left), the force scaling task (middle), as well as the difference between the two tasks (rotation minus scaling, right). All plots include data from both hands of each participant. (A–C) shows the grand average FDI:APB coherence in each task. (D–F) shows the interquartile range (IQR) of coherence z-scores across subjects at each individual frequency sample. (G–I) shows the statistical difference between groups (expressed as Cohen's D), calculated at every frequency. The horizontal dashed lines represent a statistical threshold for identifying frequency bands of interest for further statistical testing (see Methods). The gray bars show the bands of interest identified for each condition, with an * indicating that the band as a whole differed significantly between groups (p < 0.05, corrected for multiple comparisons).
Figures 3A–C show the consistency of alpha-band differences in intermuscular coherence between controls and participants with PD, for each task separately, as well as and their difference. For these plots, the coherence values are averaged across both hands per individual, since we found that there were no group differences (PD vs. CT) in the laterality of coherence between right and left hands (p = 0.55, p = 0.36, and p = 0.57 for rotation, scaling, and their difference, respectively). For consistency, we defined the alpha-band per task according to the precise frequency ranges identified in our SPM analysis. Using a generic 6–15 Hz frequency range for all tasks produced nearly identical results. Generally, a substantial degree of overlap between groups was apparent due to high inter-subject variability. Nonetheless, this binned/averaged alpha-band coherence measure showed significant differences between groups (p = 0.04, p = 0.02) for rotation and scaling tasks, respectively. Within each group, coherence in the alpha-band was task-dependent, differing significantly between scaling and rotation tasks (p = 0.008 and p = 0.005 for control and PD groups, respectively), but again, the magnitude of this difference did not differ between groups (p = 0.95). Coherence differences between groups had Cohen's D effect sizes of 0.75, 0.91, and 0.58, for the rotation task, the scaling task, and their difference, respectively.
Figure 3. Individual alpha-band coherence profiles differ between groups and correlate with UPDRS. (A–C) Box plots summarizing the mean coherence values for each participant calculated over the frequency range of significant difference (both alpha-band, see Figure 2). Both rotation (A) and scaling (B), but not the difference between tasks (C), showed significantly (p < 0.05) larger values in the Parkinson's Disease group relative to controls. (D–F) Total UPDRS scores vs. mean coherence (as plotted above). The correlation between UPDRS scores and alpha-band coherence was strongest for the rotation task and statistically significant (D, rho = 0.058, p = 0.018, see Table 2). + symbols indicate outliers. ° symbols show the individual data from each participant, as summarized by the box plots. * symbols show individual data from each participant with PD.
Within the PD group, it is possible that coherence within the alpha-band covaries with symptom severity. We therefore tested the hypothesis of a non-zero correlation between coherence (as plotted in Figures 3A–C) and the total UPDRS score for each individual (Figures 3D–F and Table 2). Table 2 shows the Spearman correlation coefficients between coherence in the alpha-band and the total UPDRS score, as well as the two main subsections of the UPDRS score: the motor symptoms score (Part III), and the self-evaluation of daily activities (part II). The 95% confidence interval around each correlation value is shown as well. Moderate correlations were found for coherence within the rotation task but not the force scaling task or their difference. Overall, the correlation between intermuscular coherence and symptom severity was not restricted to (or strongest for) the motor-only portion of the scale (UPDRS III), suggesting a more general relationship with disease state.
Table 2. Spearman's correlation (rho) between mean coherence within the alpha-band, and patient symptom severities, as per the UPDRS, along with 95% confidence interval (CI) boundaries.
To determine if the accuracy with which controls and participants with PD performed the visuo-motor tracking related to coherence metrics, we compared the standard deviation of rotation angle or pinch force (after removing the 0.1 Hz target frequency), separately across groups. Task performance did not differ across groups (p = 0.67 and p = 0.86, respectively). This implies that any changes in alpha-band neural drive did not substantially contribute to tracking error, and indeed, neither group showed a significant correlation between tracking performance and alpha-band coherence in either task.
Finally, we sought to determine if a single coherence index calculated for each individual could be used to discriminate between groups. For this analysis, we calculated the proportion of total coherence (2 to 80 Hz) falling within the alpha-band for each task/hand (the “alpha-ratio”). This is a simple way to reduce inter-subject sources of variability and focus on the contribution of alpha-band coherence to the overall shape of the coherence spectra. We confirmed that the alpha-ratio showed no group differences in laterality between hands (p = 0.41 and p = 0.8 for rotation and scaling tasks, respectively), and therefore averaged across hands to obtain a single ratio for each individual, per task. Figures 4A,B shows a clear separation between PD and control groups for the rotation (p < 0.001, Cohen's D = 1.96) and scaling task (p = 0.001, Cohen's D = 1.22). Since the purpose of this test was to determine the potential for coherence measures to discriminate between groups, we did not include an evaluation of the change in coherence between tasks, as there were no group differences in this measure. Instead, we averaged the alpha-ratio values obtained from the scaling and rotation tasks together (Figure 4C) to obtain a single, combined alpha-ratio per individual, which also differed significantly between groups (p < 0.001 Cohen's D = 1.84). We constructed Receiver Operating Characteristic Curves for the alpha-ratios derived for each task and their combination (Figures 4D–F). To quantify the overall discriminability, we calculated the area under the curve (AUC), which yielded values of 0.9, 0.84, and 0.96 for rotation, scaling, and their combination, respectively. These high values indicate that excellent separation between groups was possible in our study population, and that nearly all patients could be correctly classified, with practically no false positives (misclassified controls).
Figure 4. Alpha-ratio (i.e., proportion of total intermuscular coherence within the alpha-band) can discriminate participants with PD from controls. Top row: Box plots showing the average alpha-ratio (across hands) for each individual. Results are shown for rotation (A), scaling (B), individually, as well as their combination (C). Both rotation (A) and scaling (B) showed reasonable separation between groups, but the combination of both (average) per individual resulted in the least overlap (C). (D–F) Receiver Operating Characteristic (ROC) curves for the metrics shown above. Area Under the Curve (AUC) values indicate strong-to-excellent discriminability in all cases.
Discussion
In this study, we describe a robust and consistent PD-related amplification of alpha-band (~6–15 Hz) intermuscular coherence between finger muscles evoked during simple dynamic actions. Slow rotation of a dial between the thumb and index finger produced the strongest coherence, but the same effect was also observed during isometric modulation of precision pinch force. Clinical ratings of symptom severity (as per UPDRS) correlated significantly with the increased magnitude of intermuscular coherence during the rotation task only. The global alpha-ratio (proportion of total coherence within the alpha-band) provided an index which allowed excellent discrimination between controls and participants with PD within our study sample. This provides valuable information for future development of simple, practical measures of neural dysfunction in PD.
Kinetic Action Tremor in PD
The disruptions in neural drive that we have characterized in this study appear to be an intermuscular component of Parkinsonian kinetic action tremor (1, 4, 23, 60). In 1963, Lance et al. noted anecdotally that in PD a 5–15 Hz tremor was (i) always present at the beginning of a muscle contraction, (ii) was sometimes sustained during the static portion of a contraction, and (iii) was usually visible within EMG traces as a synchronous “grouping” of action potentials within and across contracting muscles—even when the tremor itself was not detectable by eye. Action tremor is distinct from rest tremor [or its re-emergence during steady contraction (2)] in that is higher in frequency (6–15 vs. 3–5 Hz), enhanced during dynamic action rather than reduced, and it is not attenuated by dopamine replacement in PD (4).
There is currently no standard clinical procedure for quantifying 6–15 Hz action tremor, perhaps because it is often invisible without special equipment and can be difficult to distinguish from “normal” physiological tremor (discussed further below). Most previous literature quantifies the magnitude of tremulous activity at a particular muscle or joint [e.g., (4, 60–63)]. However, if action tremor stems from an inherently common neural drive to multiple muscles, then it makes sense to use intermuscular coherence to quantify it and avoid reliance on limb/task/person-specific forces and motions.
Although the idea that tremor-generating neural drive in PD is inherently distributed to multiple muscles has received some support (21, 22), there has been less effort to identify the particular muscles/tasks in which an intermuscular action tremor is most reliably evoked and detected. So far, an amplified intermuscular coherence within the alpha-band has been reported only for static tasks such as isometric knee extension (11) or across two sides of a puckered upper lip (7). However, the alpha-band drive underlying action tremor in PD is best observed during dynamic actions (1, 4), where additional neural control circuitry may be required (64). Also, previous findings have described amplified coherence across anatomically synergistic muscles which naturally share a substantial amount of their neural drive (17). This makes it difficult to disambiguate changes in overall signal strength from an effect of PD on the process of “binding” of muscles into functional synergies through neural synchronization (65). Our results seem to suggest that a natural intermuscular alpha-band drive is amplified in PD. Further, the modification of this drive across different tasks is preserved, indicating that the amplification is not especially task-specific. This may indicate that the neural circuit which distributes alpha-band drive across muscles is still functional, even if overdriven.
Neural Origins of Parkinsonian Kinetic Action Tremor
Action tremor has been described as an amplification of naturally present physiological tremor rather than something that exclusively emerges in PD (4, 60); a notion that is well-aligned with our current data. Because the origins of physiological tremor are likely multifactorial (27–36), the effects of PD on each potential source of alpha-band neural drive to muscles must be considered. There are at least four potential sources for alpha-band neural drive in healthy individuals, and they are not mutually exclusive: (1) the stretch reflex system, (2) the motor cortex, (3) the brainstem (specifically reticulospinal output), and (4) the cerebellum (via the cerebello-thalamo-cortical circuit). We favor the latter, according to the following rationale.
If alpha-band drive is measured within the output of a single muscle, it may simply reflect cycles of excitation around the monosynaptic Ia afferent reflex loop of that muscle (27, 33, 34, 58, 66). However, this mechanism would not easily explain coherent ~10 Hz activity across different muscles, which can be observed even between muscles in different hands in healthy adults (32, 38). Mechanical coupling among muscles may serve to link afferent feedback without need for monosynaptic connectivity, and a potential role for spinal interneuron pools in coordinating afferent activity across multiple muscles cannot be ruled out, but to our knowledge, there is no evidence that either mechanism could explain the generation of intermuscular coherence between the FDI and APB within the context of our task. Also, an origin in the peripheral reflex system would imply that stretch reflex or H-reflex amplitudes should be elevated in PD, which is not the case (67, 68). PD can, in some cases, disrupt task-dependent stretch reflex modulation (68), suggesting upstream mismanagement rather than pure peripheral dysfunction.
Studies of neural drive to individual muscles have indicated potential cortical involvement in alpha-band Parkinsonian tremor (12–18 Hz) during steady contraction of the wrist extensors (69, 70), although a subsequent delay analysis suggested that the cortical activity might not generate this drive, but instead, it may receive an efference copy from a subcortical source (69). Also, steady muscle contractions do not usually produce alpha-band corticomuscular coherence (EEG-EMG), even when intermuscular coherence (EMG-EMG) in this range is simultaneously observed (38, 71). For these reasons, it seems unlikely that our findings reflect an amplification of corticospinal output.
The reticulospinal pathway may also produce an alpha-band drive to muscles. Acoustic startle stimulates the reticulospinal pathway (72) and produces a brief wave of 12–16 Hz (73) coherent activity in bilateral pairs of muscles. However, involuntary reactions to acoustic startle may not directly predict reticulospinal contributions to intermuscular coherence during voluntary tasks. The potential association between acoustic startle responses and intermuscular coherence during voluntary tasks has never been fully explored, and current evidence is indirect. For example, spasticity in chronic stroke is thought to relate to overdriven reticulospinal drive to muscles (74, 75), and both alpha-band intermuscular coherence (76, 77) and responses to acoustic startle are increased in spastic muscles of stroke survivors (78, 79). However, if PD amplifies reticulospinal output, startle responses should be increased in magnitude or consistency relative to controls, but they are not (80). The preparation of brainstem output by upstream circuits may, however, be abnormal in PD since acoustic startle does not produce the usual acceleration of reaction times (the “StartReact” paradigm) in people with PD-associated freezing of gate (81). Therefore, the amplified alpha-band drive in PD may involve reticulospinal pathways, but probably not through an overall change in the excitability of the pontomedullary reticular formation.
Finally, the frequency spectrum of neural drive to muscles may be substantially influenced by the cerebellum via its actions within the cerebello-thalamo-cortical (CTC) circuit. The cerebellum oscillates in synchrony with the motor cortex at frequencies between 10 and 40 Hz (82), and slow movements of the finger produce coherent alpha-band oscillations between muscle activity and the cerebellum, thalamus, and motor cortex (30). In addition to a possible influence on the cortex and corticospinal output to muscles, the deep cerebellar nuclei may also transmit a 10 Hz oscillation to muscles through modulation of the brainstem and reticulospinal tract output (82–84). In theory, output from the deep cerebellar nuclei to the cortex and brainstem could be amplified (disinhibited) in PD since Purkinje cell damage is associated with the disease (42, 85, 86). Regardless of the precise pathway to muscles, a wide variety of studies have implicated general CTC dysfunction in PD (40–42, 87–89). Further, dysfunction of the CTC is consistently associated with the generation of tremor in both PD (90–92) and the more common, but often-associated (93) condition known as essential tremor (94–99). In fact, medication in PD may reduce rest tremor by reconnecting cerebellar communication with the cortex through the ventrolateral thalamus (92), and a similar mechanism may explain the levodopa-reversible loss of 40 Hz corticomuscular drive in PD (12), since inactivation of the cerebellar cortex disrupts gamma-band cortical activity (100). In this view, basal ganglia dysfunction, as in PD, almost inevitably implies cerebellar network dysfunction. In fact, these structures have been described recently as nodes within a larger, tightly-interconnected network (101).
Taken together, it is clear that damage to this cerebellar circuit is not only plausible within our study population, but would very likely produce changes in alpha-band drive to muscles. Given the strong associations between cerebellar activity, cortical activity, and alpha-band spindle discharge during slow movements (102), it is likely that cerebellar-circuit dysfunction could both directly and indirectly influence the frequency spectrum of neuromuscular oscillations in PD. In addition to shaping neural oscillations, cerebellar circuit dysfunction can be expected to influence a variety of non-motor functions as well (86, 91, 103–106) and this might explain why correlations between coherence and UPDRS scores were not restricted to the direct motor evaluation in our study. That said, a specific or exclusive role for the CTC in the present study not possible to determine, and will require further research.
Intermuscular Coherence as a Potential Clinically-Practical Biomarker in PD
We found that the proportion of total coherence within the ~6–15 Hz range was sufficient to statistically separate groups in either task. When all ratio values obtained from each individual were averaged together, it allowed for excellent discrimination between our two groups (AUC of 0.96). If this level of discriminability were to generalize to the larger population, it would meet or exceed the performance of the best MRI-based diagnostic biomarkers to date (107, 108) and yet not require expensive neuroimaging, invasive collection of cerebrospinal fluid, extensive patient history, etc. To be clear though, much larger studies would be needed to evaluate the readiness of any putative biomarker for clinical translation. Our intent in exploring discriminability between groups was not to establish an application-ready feature set or threshold for optimal real-world classification, but rather to provide evidence that EMG-based measures have clinical potential, since they have rarely been explored in this capacity. In fact, the only other EMG-derived measure which has achieved such strong discrimination between PD and controls was calculated from the cross-trial distribution of EMG burst durations in the biceps during fast elbow flexions in temporarily-off-medication patients (57). The neural mechanisms that contributed to that effect were not identified, and the measurement equipment and procedure itself was somewhat specialized and not well-suited for widespread clinical implementation. Our simpler method, with further refinements, seems more practical for common application and to complement other non-invasive biomarkers, such as those derived from olfaction (109, 110), or finger movements (111, 112).
Limitations
An important consideration in this study is that we recorded from participants who were on their normal medication, as this was a sample-of-convenience. To our knowledge, there is no evidence that medication could explain the increased alpha-band neural drive we observed, especially since tremor in this frequency range is not enhanced or attenuated by medication (4, 113). Of course, if it were attenuated, then presumably unmedicated participants would show even stronger differences from controls, strengthening our classification ability. Medication may have normalized 30–50 Hz intermuscular coherence in our study, as this frequency band is dopamine-dependent (12), but again, if our analysis were applied to unmedicated participants we would expect better discrimination between PD and controls, not worse. In PD, dopamine depletion in the basal ganglia is already extreme by the time of symptom onset (e.g., ~80% loss in the putamen) making it difficult to detect or study the early stages of the disease (114). Thus, a method capable of detecting both low dopamine and amplified physiological tremor could contribute to such efforts, especially if combined with PD-sensitive metrics of finger movement speed/timing (111, 112). This would, in turn, enable the early initiation of neuroprotective strategies once they become available. Of course, any claims concerning translational/diagnostic capabilities at this point is speculative. Our sample size does not allow us to accurately characterize the strength and variability of effects across the total population of people with PD. Effect size estimates such as Cohen's D may be inflated in smaller studies, while skewed distributions or outliers may have the opposite effect. That said, our study is appropriately powered to detect the strong effects that are most likely to be of scientific and translational relevance, as well as to justify/enable their investigation in larger cohorts. Additionally, determining the optimal battery of task parameters will require future work. For example, it may be that features of visual feedback, the speed of the required dynamic action, the particular muscles measured, etc. may all be important variables. Also, this study focused on kinetic action tremor, however, the postural form of action tremor which emerges during static contractions may yield complementary information. A final limitation regarding intermuscular coherence analysis is that it measures the relative proportion of shared vs. total EMG signal variance, and thus changes to both intramuscular (muscle-specific) and intermuscular (shared) neural drive can alter coherence. Separation of shared and muscle-specific neural drive can be accomplished using methods based on motor unit analysis [e.g., Laine et al., (17)], and may represent a valuable direction for future investigation. At present though, it is clear that tasks designed to evoke kinetic action tremor from functionally-coordinated muscles provide a clear window into the nervous system, which holds value for the assessment of PD and its progression.
Data Availability Statement
The datasets generated for this study are available on request to the corresponding author.
Ethics Statement
The studies involving human participants were reviewed and approved by University of Southern California Internal Review Board. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
CL and FV-C contributed to the study design, interpretation of data, editing, and approval of the final manuscript. CL contributed to collection of experimental data, data analysis, and drafting of the manuscript.
Funding
This project was supported by NIH: NIAMS R01 AR-050520 and R01 AR-052345 to FV-C, and by: Department of Defense CDMRP Grant MR150091 and Award W911NF1820264 from the DARPA-L2M program to FV-C. The content of this endeavor is solely the responsibility of the authors and does not represent the official views of the National Institutes of Health, or the Department of Defense.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
The authors thank Catherine Parsekian, Lena Sulzenbacher, Shaohui Qian, Jamie Flores, and Kian Jalaleddini for assistance with subject recruitment and data collection, and Dr. Jun Yong Shin for creating the 3D-printed dial used in the study.
References
1. Lance JW, Schwab RS, Peterson EA. Action tremor and the cogwheel phenomenon in Parkinson's disease. Brain. (1963) 86:95–110. doi: 10.1093/brain/86.1.95
2. Jankovic J, Schwartz KS, Ondo W. Re-emergent tremor of Parkinson's disease. J Neurol Neurosurg Psychiatry. (1999) 67:646–50. doi: 10.1136/jnnp.67.5.646
3. Timmermann L, Gross J, Dirks M, Volkmann J, Freund H-J, Schnitzler A. The cerebral oscillatory network of parkinsonian resting tremor. Brain J Neurol. (2003) 126:199–212. doi: 10.1093/brain/awg022
4. Raethjen J, Pohle S, Govindan RB, Morsnowski A, Wenzelburger R, Deuschl G. Parkinsonian action tremor: interference with object manipulation and lacking levodopa response. Exp Neurol. (2005) 194:151–60. doi: 10.1016/j.expneurol.2005.02.008
5. Jankovic J. Parkinson's disease: clinical features and diagnosis. J Neurol Neurosurg Psychiatry. (2008) 79:368–76. doi: 10.1136/jnnp.2007.131045
6. Rana AQ, Saleh M. Relationship between resting and action tremors in Parkinson's disease. J Neurosci Rural Pract. (2016) 7:232–7. doi: 10.4103/0976-3147.176192
7. Caviness JN, Liss JM, Adler C, Evidente V. Analysis of high-frequency electroencephalographic-electromyographic coherence elicited by speech and oral non-speech tasks in Parkinson's disease. J Speech Lang Hear Res. (2006) 49:424. doi: 10.1044/1092-4388(2006/033)
8. Hammond C, Bergman H, Brown P. Pathological synchronization in Parkinson's disease: networks, models and treatments. Trends Neurosci. (2007) 30:357–64. doi: 10.1016/j.tins.2007.05.004
9. Little S, Pogosyan A, Kuhn AA, Brown P. Beta band stability over time correlates with Parkinsonian rigidity and bradykinesia. Exp Neurol. (2012) 236:383–8. doi: 10.1016/j.expneurol.2012.04.024
10. Little S, Brown P. The functional role of beta oscillations in Parkinson's disease. Parkinsonism Relat. Disord. (2014) 20(Suppl.1):S44–8. doi: 10.1016/S1353-8020(13)70013-0
11. Flood MW, Jensen BR, Malling A-S, Lowery MM. Increased EMG intermuscular coherence and reduced signal complexity in Parkinson's disease. Clin Neurophysiol. (2019) 130:259–69. doi: 10.1016/j.clinph.2018.10.023
12. McAuley JH. Levodopa reversible loss of the Piper frequency oscillation component in Parkinson's disease. J Neurol Neurosurg Psychiatry. (2001) 70:471–5. doi: 10.1136/jnnp.70.4.471
13. Kilner JM, Baker SN, Salenius S, Jousmäki V, Hari R, Lemon RN. Task-dependent modulation of 15–30 Hz coherence between rectified EMGs from human hand and forearm muscles. J Physiol. (1999) 516:559. doi: 10.1111/j.1469-7793.1999.0559v.x
14. Boonstra TW, van Wijk B, Praamstra P, Daffertshofer A. Corticomuscular and bilateral EMG coherence reflect distinct aspects of neural synchronization. Neurosci Lett. (2009) 463:17–21. doi: 10.1016/j.neulet.2009.07.043
15. Boonstra TW. The potential of corticomuscular and intermuscular coherence for research on human motor control. Front Hum Neurosci. (2013) 7:855. doi: 10.3389/fnhum.2013.00855
16. DeMarchis C, Severini G, Castronovo AM, Schmid M, Conforto S. Intermuscular coherence contributions in synergistic muscles during pedaling. Exp Brain Res. (2015) 233:1907–19. doi: 10.1007/s00221-015-4262-4
17. Laine CM, Martinez-Valdes E, Falla D, Mayer F, Farina D. Motor neuron pools of synergistic thigh muscles share most of their synaptic input. J Neurosci. (2015) 35:12207–16. doi: 10.1523/JNEUROSCI.0240-15.2015
18. Farina D, Negro F, Muceli S, Enoka RM. Principles of motor unit physiology evolve with advances in technology. Physiology. (2016) 31:83–94. doi: 10.1152/physiol.00040.2015
19. Laine CM, Valero-Cuevas FJ. Intermuscular coherence reflects functional coordination. J Neurophysiol. (2017) 118:1775–83. doi: 10.1152/jn.00204.2017
20. Boonstra TW, Danna-Dos-Santos A, Xie HB, Roerdink M, Stins JF, Breakspear M. Muscle networks: connectivity analysis of EMG activity during postural control. Sci Rep. (2015) 5:17830. doi: 10.1038/srep17830
21. van der Stouwe AMM, Conway BA, Elting JW, Tijssen MAJ, Maurits NM. Usefulness of intermuscular coherence and cumulant analysis in the diagnosis of postural tremor. Clin Neurophysiol. (2015) 126:1564–9. doi: 10.1016/j.clinph.2014.10.157
22. He X, Hao M-Z, Wei M, Xiao Q, Lan N. Contribution of inter-muscular synchronization in the modulation of tremor intensity in Parkinson's disease. J NeuroEngineering Rehabil. (2015) 12:108. doi: 10.1186/s12984-015-0101-x
23. Beuter A, Edwards R. Kinetic tremor during tracking movements in patients with Parkinson's disease. Parkinsonism Relat Disord. (2002) 8:361–8. doi: 10.1016/S1353-8020(01)00051-7
24. Kraus PH, Lemke MR, Reichmann H. Kinetic tremor in Parkinson's disease—an underrated symptom. J Neural Transm. (2006) 113:845–53. doi: 10.1007/s00702-005-0354-9
25. Gironell A, Pascual-Sedano B, Aracil I, Marín-Lahoz J, Pagonabarraga J, Kulisevsky J. Tremor types in Parkinson disease: a descriptive study using a new classification. Park Dis. (2018) 2018:5. doi: 10.1155/2018/4327597
26. Forssberg H, Ingvarsson PE, Iwasaki N, Johansson RS, Gordon AM. Action tremor during object manipulation in Parkinson's disease. Mov Disord. (2000) 15:244–54. doi: 10.1002/1531-8257(200003)15:2<244::AID-MDS1007>3.0.CO;2-H
27. Lippold OCJ. Oscillation in the stretch reflex arc and the origin of the rhythmical, 8–12 c/s component of physiological tremor. J Physiol. (1970) 206:359–82. doi: 10.1113/jphysiol.1970.sp009018
28. Elble RJ, Randall JE. Motor-unit activity responsible for 8- to 12-Hz component of human physiological finger tremor. J Neurophysiol. (1976) 39:370–83. doi: 10.1152/jn.1976.39.2.370
29. McAuley JH, Marsden CD. Physiological and pathological tremors and rhythmic central motor control. Brain. (2000) 123:1545–67. doi: 10.1093/brain/123.8.1545
30. Gross J, Timmermann L, Kujala J, Dirks M, Schmitz F, Salmelin R, et al. The neural basis of intermittent motor control in humans. Proc Natl Acad Sci USA. (2002) 99:2299–302. doi: 10.1073/pnas.032682099
31. Elble RJ. Characteristics of physiologic tremor in young and elderly adults. Clin Neurophysiol. (2003) 114:624–35. doi: 10.1016/S1388-2457(03)00006-3
32. Evans CMB, Baker SN. Task-dependent intermanual coupling of 8-Hz discontinuities during slow finger movements. Eur J Neurosci. (2003) 18:453–6. doi: 10.1046/j.1460-9568.2003.02751.x
33. Christakos CN, Papadimitriou NA, Erimaki S. Parallel neuronal mechanisms underlying physiological force tremor in steady muscle contractions of humans. J Neurophysiol. (2006) 95:53–66. doi: 10.1152/jn.00051.2005
34. Erimaki S, Christakos CN. Coherent motor unit rhythms in the 6–10 Hz range during time-varying voluntary muscle contractions: neural mechanism and relation to rhythmical motor control. J Neurophysiol. (2008) 99:473–83. doi: 10.1152/jn.00341.2007
35. Lakie M, Vernooij CA, Osborne TM, Reynolds RF. The resonant component of human physiological hand tremor is altered by slow voluntary movements. J Physiol. (2012) 590:2471–83. doi: 10.1113/jphysiol.2011.226449
36. Vernooij CA, Reynolds RF, Lakie M. A dominant role for mechanical resonance in physiological finger tremor revealed by selective minimization of voluntary drive and movement. J Neurophysiol. (2013) 109:2317–26. doi: 10.1152/jn.00926.2012
37. Nazarpour K, Barnard A, Jackson A. Flexible cortical control of task-specific muscle synergies. J Neurosci. (2012) 32:12349–60. doi: 10.1523/JNEUROSCI.5481-11.2012
38. de Vries IEJ, Daffertshofer A, Stegeman DF, Boonstra TW. Functional connectivity in the neuromuscular system underlying bimanual coordination. J Neurophysiol. (2016) 116:2576–85. doi: 10.1152/jn.00460.2016
39. Racz K, Brown D, Valero-Cuevas FJ. An involuntary stereotypical grasp tendency pervades voluntary dynamic multifinger manipulation. J Neurophysiol. (2012) 108:2896–911. doi: 10.1152/jn.00297.2012
40. Sen S, Kawaguchi A, Truong Y, Lewis MM, Huang X. Dynamic changes in cerebello-thalamo-cortical motor circuitry during progression of Parkinson's disease. Neuroscience. (2010) 166:712–9. doi: 10.1016/j.neuroscience.2009.12.036
41. Wu T, Hallett M. The cerebellum in Parkinson's disease. Brain. (2013) 136:696–709. doi: 10.1093/brain/aws360
42. Muthuraman M, Raethjen J, Koirala N, Anwar AR, Mideksa KG, Elble R, et al. Cerebello-cortical network fingerprints differ between essential, Parkinson's and mimicked tremors. Brain. (2018) 141:1770–81. doi: 10.1093/brain/awy098
43. Potvin JR, Brown SHM. Less is more: high pass filtering, to remove up to 99% of the surface EMG signal power, improves EMG-based biceps brachii muscle force estimates. J Electromyogr Kinesiol Off J Int Soc Electrophysiol Kinesiol. (2004) 14:389–99. doi: 10.1016/j.jelekin.2003.10.005
44. Riley ZA, Terry ME, Mendez-Villanueva A, Litsey JC, Enoka RM. Motor unit recruitment and bursts of activity in the surface electromyogram during a sustained contraction. Muscle Nerve. (2008) 37:745–53. doi: 10.1002/mus.20978
45. Brown SHM, Brookham RL, Dickerson CR. High-pass filtering surface EMG in an attempt to better represent the signals detected at the intramuscular level. Muscle Nerve. (2009). doi: 10.1002/mus.21470
46. Staudenmann D, Roeleveld K, Stegeman DF, van Dieën JH. Methodological aspects of SEMG recordings for force estimation—A tutorial and review. J Electromyogr Kinesiol. (2010) 20:375–87. doi: 10.1016/j.jelekin.2009.08.005
47. Boonstra TW, Breakspear M. Neural mechanisms of intermuscular coherence: implications for the rectification of surface electromyography. J Neurophysiol. (2012) 107:796–807. doi: 10.1152/jn.00066.2011
48. Dideriksen JL, Farina D. Amplitude cancellation influences the association between frequency components in the neural drive to muscle and the rectified EMG signal. PLOS Comput Biol. (2019) 15:e1006985. doi: 10.1371/journal.pcbi.1006985
49. Rosenberg JR, Amjad AM, Breeze P, Brillinger DR, Halliday DM. The Fourier approach to the identification of functional coupling between neuronal spike trains. Prog Biophys Mol Biol. (1989) 53:1–31. doi: 10.1016/0079-6107(89)90004-7
50. Terry K, Griffin L. How computational technique and spike train properties affect coherence detection. J Neurosci Methods. (2008) 168:212–23. doi: 10.1016/j.jneumeth.2007.09.014
51. Kattla S, Lowery MM. Fatigue related changes in electromyographic coherence between synergistic hand muscles. Exp Brain Res Exp Hirnforsch Expérimentation Cérébrale. (2010) 202:89–99. doi: 10.1007/s00221-009-2110-0
52. Baker SN, Pinches EM, Lemon RN. Synchronization in monkey motor cortex during a precision grip task. II effect of oscillatory activity on corticospinal output. J Neurophysiol. (2003) 89:1941–53. doi: 10.1152/jn.00832.2002
53. Holmes AP, Blair RC, Watson JDG, Ford I. Non-parametric analysis of statistic images from functional mapping experiments. J Cereb Blood Flow Metab. (1996) 16:7–22. doi: 10.1097/00004647-199601000-00002
54. Nichols TE, Holmes AP. Non-parametric permutation tests for functional neuroimaging: a primer with examples. Hum Brain Mapp. (2002) 15:1–25. doi: 10.1002/hbm.1058
55. Maris E, Oostenveld R. Non-parametric statistical testing of EEG- and MEG-data. J Neurosci Methods. (2007) 164:177–90. doi: 10.1016/j.jneumeth.2007.03.024
56. Eng J. Receiver operating characteristic analysis. Acad Radiol. (2005) 12:909–16. doi: 10.1016/j.acra.2005.04.005
57. Robichaud JA, Pfann KD, Leurgans S, Vaillancourt DE, Comella CL, Corcos DM. Variability of EMG patterns: a potential neurophysiological marker of Parkinson's disease? Clin Neurophysiol. (2009) 120:390–7. doi: 10.1016/j.clinph.2008.10.015
58. Laine CM, Yavuz SU, Farina D. Task-related changes in sensorimotor integration influence the common synaptic input to motor neurones. Acta Physiol. (2014) 211:229–39. doi: 10.1111/apha.12255
59. Laine CM, Nagamori A, Valero-Cuevas FJ. The dynamics of voluntary force production in afferented muscle influence involuntary tremor. Front Comput Neurosci. (2016) 10:86. doi: 10.3389/fncom.2016.00086
60. Christakos CN, Erimaki S, Anagnostou E, Anastasopoulos D. Tremor-related motor unit firing in Parkinson's disease: implications for tremor genesis. J Physiol. (2009) 587:4811–27. doi: 10.1113/jphysiol.2009.173989
61. Brown P, Corcos DM, Rothwell JC. Does parkinsonian action tremor contribute to muscle weakness in Parkinson's disease? Brain. (1997) 120:401–8. doi: 10.1093/brain/120.3.401
62. Teräväinen H, Calne DB. Action tremor in Parkinson's disease. J Neurol Neurosurg Psychiatry. (1980) 43:257–63. doi: 10.1136/jnnp.43.3.257
63. Brown P, Corcos DM, Rothwell JC. Action tremor and weakness in Parkinson's disease: a study of the elbow extensors. Mov Disord. (1998) 13:56–60. doi: 10.1002/mds.870130114
64. Oya T, Takei T, Seki K. Distinct sensorimotor feedback loops for dynamic and static control of primate precision grip. bioRxiv. (2019). doi: 10.1101/640201
65. Farmer SF. Rhythmicity, synchronization and binding in human and primate motor systems. J Physiol. (1998) 509:3. doi: 10.1111/j.1469-7793.1998.003bo.x
66. Nagamori A, Laine CM, Valero-Cuevas FJ. Cardinal features of involuntary force variability can arise from the closed-loop control of viscoelastic afferented muscles. PLOS Comput Biol. (2018) 14:e1005884. doi: 10.1371/journal.pcbi.1005884
67. Cody FWJ, Macdermott N, Matthews PBC, Richardson HC. Observations on the genesis of the stretch reflex in Parkinson's disease. Brain. (1986) 109:229–49. doi: 10.1093/brain/109.2.229
68. Hayashi R, Tokuda T, Tako K, Yanagisawa N. Impaired modulation of tonic muscle activities and H-reflexes in the soleus muscle during standing in patients with Parkinson's disease. J Neurol Sci. (1997) 153:61–7. doi: 10.1016/S0022-510X(97)00175-5
69. Salenius S, Avikainen S, Kaakkola S, Hari R, Brown P. Defective cortical drive to muscle in Parkinson's disease and its improvement with levodopa. Brain. (2002) 125:491–500. doi: 10.1093/brain/awf042
70. Caviness JN, Shill HA, Sabbagh MN, Evidente VGH, Hernandez JL, Adler CH. Corticomuscular coherence is increased in the small postural tremor of Parkinson's disease. Mov Disord. (2006) 21:492–9. doi: 10.1002/mds.20743
71. Reyes A, Laine CM, Kutch JJ, Valero-Cuevas FJ. Beta band corticomuscular drive reflects muscle coordination strategies. Front Comput Neurosci. (2017) 11:17. doi: 10.3389/fncom.2017.00017
72. Rothwell JC. The startle reflex, voluntary movement, and the reticulospinal tract. Suppl Clin Neurophysiol. (2006) 58:223–31. doi: 10.1016/S1567-424X(09)70071-6
73. Grosse P, Brown P. Acoustic startle evokes bilaterally synchronous oscillatory EMG activity in the healthy human. J Neurophysiol. (2003) 90:1654–61. doi: 10.1152/jn.00125.2003
74. Dewald JPA, Pope PS, Given JD, Buchanan TS, Rymer WZ. Abnormal muscle coactivation patterns during isometric torque generation at the elbow and shoulder in hemiparetic subjects. Brain. (1995) 118:495–510. doi: 10.1093/brain/118.2.495
75. Krakauer JW, Carmichael ST. Broken Movement: The Neurobiology of Motor Recovery After Stroke. Cambridge, MA: The MIT Press (2017). doi: 10.7551/mitpress/9310.001.0001
76. Lan Y, Yao J, Dewald JPA. Reducing the impact of shoulder abduction loading on the classification of hand opening and grasping in individuals with poststroke flexion synergy. Front Bioeng Biotechnol. (2017) 5:39. doi: 10.3389/fbioe.2017.00039
77. Chen Y-T, Li S, Magat E, Zhou P, Li S. Motor overflow and spasticity in chronic stroke share a common pathophysiological process: analysis of within-limb and between-limb EMG-EMG coherence. Front Neurol. (2018) 9:795. doi: 10.3389/fneur.2018.00795
78. Li S, Chang S-H, Francisco GE, Verduzco-Gutierrez M. Acoustic startle reflex in patients with chronic stroke at different stages of motor recovery: a pilot study. Top Stroke Rehabil. (2014) 21:358–70. doi: 10.1310/tsr2104-358
79. Li S, Bhadane M, Gao F, Zhou P. The reticulospinal pathway does not increase its contribution to the strength of contralesional muscles in stroke survivors as compared to ipsilesional side or healthy controls. Front Neurol. (2017) 8:627. doi: 10.3389/fneur.2017.00627
80. Kofler M, Müller J, Wenning GK, Reggiani L, Hollosi P, Bösch S, et al. The auditory startle reaction in parkinsonian disorders. Mov Disord. (2001) 16:62–71. doi: 10.1002/1531-8257(200101)16:1<62::AID-MDS1002>3.0.CO;2-V
81. Nonnekes J, Geurts ACH, Oude Nijhuis LB, van Geel K, Snijders AH, Bloem BR, et al. Reduced StartReact effect and freezing of gait in Parkinson's disease: two of a kind? J Neurol. (2014) 261:943–50. doi: 10.1007/s00415-014-7304-0
82. Soteropoulos DS. Cortico-cerebellar coherence during a precision grip task in the monkey. J Neurophysiol. (2005) 95:1194–206. doi: 10.1152/jn.00935.2005
83. Soteropoulos DS, Baker SN. Bilateral representation in the deep cerebellar nuclei. J Physiol. (2008) 586:1117–36. doi: 10.1113/jphysiol.2007.144220
84. Williams ER, Soteropoulos DS, Baker SN. Spinal interneuron circuits reduce approximately 10-Hz movement discontinuities by phase cancellation. Proc Natl Acad Sci USA. (2010) 107:11098–103. doi: 10.1073/pnas.0913373107
85. Kishore A, Meunier S, Popa T. Cerebellar influence on motor cortex plasticity: behavioral implications for Parkinson's disease. Front Neurol. (2014) 5:68. doi: 10.3389/fneur.2014.00068
86. O'Callaghan C, Hornberger M, Balsters JH, Halliday GM, Lewis SJG, Shine JM. Cerebellar atrophy in Parkinson's disease and its implication for network connectivity. Brain. (2016) 139:845–55. doi: 10.1093/brain/awv399
87. Lewis MM, Slagle CG, Smith DB, Truong Y, Bai P, McKeown M, et al. Task specific influences of Parkinson's disease on the striato-thalamo-cortical and cerebello-thalamo-cortical motor circuitries. Neuroscience. (2007) 147:224–35. doi: 10.1016/j.neuroscience.2007.04.006
88. Zhang J, Wei L, Hu X, Xie B, Zhang Y, Wu G-R, et al. Akinetic-rigid and tremor-dominant Parkinson's disease patients show different patterns of intrinsic brain activity. Parkinsonism Relat Disord. (2015) 21:23–30. doi: 10.1016/j.parkreldis.2014.10.017
89. Piccinin CC, Campos LS, Guimarães RP, Piovesana LG, Santos MCA, dos Azevedo PC, et al. Differential pattern of cerebellar atrophy in tremor-predominant and akinetic/rigidity-predominant Parkinson's disease. Cerebellum. (2017) 16:623–8. doi: 10.1007/s12311-016-0834-5
90. Dirkx MF, Ouden H, den Aarts E, Timmer M, Bloem BR, Toni I, et al. The cerebral network of Parkinson's tremor: an effective connectivity fMRI study. J Neurosci. (2016) 36:5362–72. doi: 10.1523/JNEUROSCI.3634-15.2016
91. Caligiore D, Pezzulo G, Baldassarre G, Bostan AC, Strick PL, Doya K, et al. Consensus paper: towards a systems-level view of cerebellar function: the interplay between cerebellum, basal ganglia, and cortex. Cerebellum. (2017) 16:203–29. doi: 10.1007/s12311-016-0763-3
92. Dirkx MF, Ouden D, M HE, Aarts E, Timmer MHM, Bloem BR, et al. Dopamine controls Parkinson's tremor by inhibiting the cerebellar thalamus. Brain. (2017) 140:721–34. doi: 10.1093/brain/aww331
93. Tarakad A, Jankovic J. Essential tremor and Parkinson's disease: exploring the relationship. Tremor Hyperkinetic Mov. (2019) 8:589. doi: 10.7916/D8MD0GVR
94. Marsden JF, Ashby P, Limousin-Dowsey P, Rothwell JC, Brown P. Coherence between cerebellar thalamus, cortex and muscle in manCerebellar thalamus interactions. Brain. (2000) 123:1459–70. doi: 10.1093/brain/123.7.1459
95. Minks E, Mareček R, Pavlík T, Ovesná P, Bareš M. Is the cerebellum a potential target for stimulation in Parkinson's disease? results of 1-Hz rTMS on upper limb motor tasks. Cerebellum. (2011) 10:804–11. doi: 10.1007/s12311-011-0290-1
96. Bareš M., Husárová I, Lungu OV. Essential tremor, the cerebellum, and motor timing: towards integrating them into one complex entity. Tremor Hyperkinetic Mov. (2012) 2.
97. Buijink AWG, Broersma M, van der Stouwe AMM, van Wingen GA, Groot PFC, Speelman JD, et al. Rhythmic finger tapping reveals cerebellar dysfunction in essential tremor. Parkinsonism Relat Disord. (2015) 21:383–8. doi: 10.1016/j.parkreldis.2015.02.003
98. Louis ED. Essential tremor: a common disorder of purkinje neurons? Neurosci Rev J Bringing Neurobiol Neurol Psychiatry. (2016) 22:108–18. doi: 10.1177/1073858415590351
99. Louis ED. Essential tremor then and now: how views of the most common tremor diathesis have changed over time. Parkinsonism Relat Disord. (2018) 46:S70–4. doi: 10.1016/j.parkreldis.2017.07.010
100. Popa D, Spolidoro M, Proville RD, Guyon N, Belliveau L, Léna C. Functional role of the cerebellum in gamma-band synchronization of the sensory and motor cortices. J Neurosci. (2013) 33:6552–6. doi: 10.1523/JNEUROSCI.5521-12.2013
101. Bostan AC, Strick PL. The basal ganglia and the cerebellum: nodes in an integrated network. Nat Rev Neurosci. (2018) 19:338–50. doi: 10.1038/s41583-018-0002-7
102. Schieber MH, Thach WT. Trained slow tracking. II. Bidirectional discharge patterns of cerebellar nuclear, motor cortex, and spindle afferent neurons. J Neurophysiol. (1985) 54:1228–70. doi: 10.1152/jn.1985.54.5.1228
103. Middleton FA, Strick PL. Basal ganglia and cerebellar loops: motor and cognitive circuits. Brain Res Rev. (2000) 31:236–50. doi: 10.1016/S0165-0173(99)00040-5
104. Tedesco AM, Chiricozzi FR, Clausi S, Lupo M, Molinari M, Leggio MG. The cerebellar cognitive profile. Brain. (2011) 134:3672–86. doi: 10.1093/brain/awr266
105. Gao L, Zhang J, Hou Y, Hallett M, Chan P, Wu T. The cerebellum in dual-task performance in Parkinson's disease. Sci Rep. (2017) 7:45662. doi: 10.1038/srep45662
106. Schmahmann JD. The cerebellum and cognition. Neurosci Lett. (2019) 688:62–75. doi: 10.1016/j.neulet.2018.07.005
107. Castellanos G, Fernández-Seara MA, Lorenzo-Betancor O, Ortega-Cubero S, Puigvert M, Uranga J, et al. Automated neuromelanin imaging as a diagnostic biomarker for Parkinson's disease. Mov Disord Off J Mov Disord Soc. (2015) 30:945–52. doi: 10.1002/mds.26201
108. Le W, Dong J, Li S, Korczyn AD. Can biomarkers help the early diagnosis of Parkinson's disease? Neurosci Bull. (2017) 33:535–42. doi: 10.1007/s12264-017-0174-6
109. Mesholam RI, Moberg PJ, Mahr RN, Doty RL. Olfaction in neurodegenerative disease: a meta-analysis of olfactory functioning in Alzheimer's and Parkinson's diseases. Arch Neurol. (1998) 55:84–90. doi: 10.1001/archneur.55.1.84
110. Haugen J, Müller MLTM, Kotagal V, Albin RL, Koeppe RA, Scott PJH, et al. Prevalence of impaired odor identification in Parkinson disease with imaging evidence of nigrostriatal denervation. J Neural Transm. (2016) 123:421–4. doi: 10.1007/s00702-016-1524-7
111. Arroyo-Gallego T, Ledesma-Carbayo MJ, Sánchez-Ferro Á, Butterworth I, Mendoza CS, Matarazzo M, et al. Detection of motor impairment in Parkinson's disease via mobile touchscreen typing. IEEE Trans Biomed Eng. (2017) 64:1994–2002. doi: 10.1109/TBME.2017.2664802
112. Shirani A, Newton BD, Okuda DT. Finger tapping impairments are highly sensitive for evaluating upper motor neuron lesions. BMC Neurol. (2017) 17:55. doi: 10.1186/s12883-017-0829-y
113. Kulisevsky J, Avila A, Barbanoj M, Antonijoan R, Torres J, Arcelus R. Levodopa does not aggravate postural tremor in Parkinson's disease. Clin Neuropharmacol. (1995) 18:435–42. doi: 10.1097/00002826-199510000-00006
Keywords: coherence, EMG, biomarker, manual tasks, alpha-band, kinetic tremor, action tremor
Citation: Laine CM and Valero-Cuevas FJ (2020) Parkinson's Disease Exhibits Amplified Intermuscular Coherence During Dynamic Voluntary Action. Front. Neurol. 11:204. doi: 10.3389/fneur.2020.00204
Received: 05 September 2019; Accepted: 09 March 2020;
Published: 03 April 2020.
Edited by:
Francesco Negro, University of Brescia, ItalyReviewed by:
Tjeerd W. Boonstra, Maastricht University, NetherlandsJakob Dideriksen, Aalborg University, Denmark
Juan Álvaro Gallego, Northwestern University, United States
Copyright © 2020 Laine and Valero-Cuevas. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Francisco J. Valero-Cuevas, dmFsZXJvJiN4MDAwNDA7dXNjLmVkdQ==