Skip to main content
medRxiv
  • Home
  • About
  • Submit
  • ALERTS / RSS
Advanced Search

Neurofeedback-guided kinesthetic motor imagery training in Parkinson’s disease: Randomized trial

View ORCID ProfileSule Tinaz, Serageldin Kamel, Sai S. Aravala, Mohamed Elfil, Ahmed Bayoumi, Amar Patel, Dustin Scheinost, Rajita Sinha, Michelle Hampson
doi: https://doi.org/10.1101/2022.01.05.22268816
Sule Tinaz
aYale University School of Medicine, Department of Neurology, Division of Movement Disorders, 15 York St, LCI 710, New Haven, CT, 06510, USA
bYale University School of Medicine, Clinical Neuroscience Imaging Center, 789 Howard Ave, New Haven, CT, 06519, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Sule Tinaz
  • For correspondence: sule.tinaz@yale.edu
Serageldin Kamel
aYale University School of Medicine, Department of Neurology, Division of Movement Disorders, 15 York St, LCI 710, New Haven, CT, 06510, USA
bYale University School of Medicine, Clinical Neuroscience Imaging Center, 789 Howard Ave, New Haven, CT, 06519, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Sai S. Aravala
aYale University School of Medicine, Department of Neurology, Division of Movement Disorders, 15 York St, LCI 710, New Haven, CT, 06510, USA
bYale University School of Medicine, Clinical Neuroscience Imaging Center, 789 Howard Ave, New Haven, CT, 06519, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Mohamed Elfil
aYale University School of Medicine, Department of Neurology, Division of Movement Disorders, 15 York St, LCI 710, New Haven, CT, 06510, USA
bYale University School of Medicine, Clinical Neuroscience Imaging Center, 789 Howard Ave, New Haven, CT, 06519, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Ahmed Bayoumi
aYale University School of Medicine, Department of Neurology, Division of Movement Disorders, 15 York St, LCI 710, New Haven, CT, 06510, USA
bYale University School of Medicine, Clinical Neuroscience Imaging Center, 789 Howard Ave, New Haven, CT, 06519, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Amar Patel
aYale University School of Medicine, Department of Neurology, Division of Movement Disorders, 15 York St, LCI 710, New Haven, CT, 06510, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Dustin Scheinost
cYale School of Medicine, Department of Radiology and Biomedical Imaging, 300 Cedar St, New Haven, CT, 06519, USA
dYale University, Department of Biomedical Engineering, 17 Hillhouse Avenue, New Haven, CT, 06520, USA
eYale School of Medicine, Child Study Center, 230 South Frontage Road, New Haven, CT, 06519, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Rajita Sinha
fYale School of Medicine, Yale Stress Center, 2 Church St South, Suite 209, New Haven, CT, 06519, USA
gYale School of Medicine, Department of Psychiatry, 300 George St, New Haven, CT, 06511, USA
hYale School of Medicine, Department of Neuroscience, 333 Cedar St, SHM-L-200, New Haven, CT, 06510, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Michelle Hampson
cYale School of Medicine, Department of Radiology and Biomedical Imaging, 300 Cedar St, New Haven, CT, 06519, USA
dYale University, Department of Biomedical Engineering, 17 Hillhouse Avenue, New Haven, CT, 06520, USA
eYale School of Medicine, Child Study Center, 230 South Frontage Road, New Haven, CT, 06519, USA
gYale School of Medicine, Department of Psychiatry, 300 George St, New Haven, CT, 06511, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • Abstract
  • Full Text
  • Info/History
  • Metrics
  • Data/Code
  • Preview PDF
Loading

Abstract

Background Parkinson’s disease (PD) causes difficulty with maintaining the speed, size, and vigor of movements, especially when they are internally generated. We previously proposed that the insula is important in motivating intentional movement via its connections with the dorsomedial frontal cortex (dmFC). We demonstrated that subjects with PD can increase the right insula-dmFC functional connectivity using fMRI-based neurofeedback (NF) combined with kinesthetic motor imagery (MI). The current study is a randomized clinical trial testing whether NF-guided kinesthetic MI training can improve motor performance and increase task-based and resting-state right insula-dmFC functional connectivity in subjects with PD.

Methods We assigned nondemented subjects with mild PD (Hoehn & Yahr stage ≤ 3) to the experimental kinesthetic MI with NF (MI-NF, n=22) and active control visual imagery (VI, n=22) groups. Only the MI-NF group received NF-guided MI training (10-12 runs). The NF signal was based on the right insula-dmFC functional connectivity strength. All subjects also practiced their respective imagery tasks at home daily for 4 weeks. Post-training changes in 1) task-based and resting-state right insula-dmFC functional connectivity were the imaging outcomes, and 2) MDS-UPDRS motor exam and motor function scores were the clinical outcomes.

Results The MI-NF group did not show significant NF regulation and was not significantly different from the VI group in any of the imaging or clinical outcome measures. The MI-NF group reported subjective improvement in kinesthetic body awareness. There was significant and comparable improvement only in motor function scores in both groups. This improvement correlated with NF regulation of the right insula-dmFC functional connectivity only in the MI-NF group. Both groups showed specific training effects in whole-brain functional connectivity with distinct neural circuits supporting kinesthetic motor and visual imagery (exploratory outcome).

Conclusions The functional connectivity-based NF regulation was unsuccessful in our cohort with mild PD. However, kinesthetic MI practice by itself or in combination with other imagery techniques is a promising tool in motor rehabilitation in PD.

Highlights

  • Parkinson’s disease (PD) causes difficulty with sustained motor performance

  • Insula and dorsomedial frontal cortex (dmFC) are implicated in motivating movement

  • Regulation of insula-dmFC functional connectivity with neurofeedback (NF) failed

  • Motor imagery practice regardless of NF improved motor function and body awareness

  • Motor imagery is a promising strategy for motor rehabilitation in PD

1. Introduction

Parkinson’s disease (PD) is a neurodegenerative disorder that causes significant motor and nonmotor disability. Loss of dopaminergic neurons in the substantia nigra pars compacta is the pathological hallmark of the disease (Braak et al., 2003). The impaired ability to sustain a steady motor performance is a major cause of morbidity in patients with PD. This is characterized by a rapid progressive decrement in the speed, amplitude, or force of movements during continuous tasks (e.g., walking, writing) (Chee et al., 2009, Kang et al., 2011, Ling et al., 2012). The standard pharmacological (e.g., dopaminergic medications) and surgical therapies (e.g., deep brain stimulation) are ineffective in improving the decrement (Baraduc et al., 2013, Bologna et al., 2018, Espay et al., 2011). It has been shown that the decrement is most pronounced when patients with PD have to internally generate movement, and improves when they are provided external cues for movement (Demirci et al., 1997, Morris et al., 2008, Tinaz et al., 2016). This suggests that there is a problem with the internal motivation of movement in patients with PD (Mazzoni et al., 2007). Neuroimaging studies typically implicate the dysfunction of motor cortical-basal ganglia circuits as the neural underpinning of the difficulty with internally-generated movement in patients with PD (Berardelli et al., 2001). The dorsomedial frontal cortex (dmFC) regions, including the supplementary motor area (SMA), pre-SMA, and cingulate motor areas, are involved in intentional motor control including movement initiation and maintenance. Neuroimaging studies have shown reduced activation in these cortical regions and the basal ganglia during internally-generated sequential movements in patients with PD (Catalan et al., 1999, Jahanshahi et al., 1995, Nakamura et al., 2001, Sabatini et al., 2000, Samuel et al., 1997, Yu et al., 2007). However, dysfunction in limbic brain circuits pertaining to the internal drive behind intentional movement may also play a role in PD. In a previous paper, we proposed that the standard cortical-basal ganglia circuit model of motor dysfunction in PD needs to be expanded to include the insula (Tinaz et al., 2018). Briefly, we argued that the insula, by processing a wide range of sensory signals arising from the body and integrating them with the emotional and motivational context, provides the impetus for intentional movement via its connections with the dmFC regions. Furthermore, we demonstrated in a proof-of-concept experiment that the functional connectivity between the right insula and dmFC can be strengthened with neurofeedback (NF)-guided kinesthetic motor imagery using functional MRI (fMRI) in subjects with PD (Tinaz et al., 2018).

In this study, we aimed to further test the findings of the proof-of-concept experiment in a randomized trial to investigate whether kinesthetic motor imagery training combined with fMRI-based NF using the right insula-dmFC functional connectivity can improve motor performance in patients with PD. We also examined whether this connectivity-based training changes the task-based and intrinsic whole-brain functional connectivity as has been demonstrated in other studies (Megumi et al., 2015, Yamashita et al., 2017).

Motor imagery refers to the mental rehearsal of motor acts without overt body movements. Imagined movements share commonalities with real movements including similar neural substrates, autonomic responses, and duration that have led to the notion of “functional equivalence” (Guillot et al., 2014). Motor imagery has been used to enhance athletic performance (Guillot and Collet, 2008) and in stroke rehabilitation (Di Rienzo et al., 2014), but rarely in a systematic way in motor rehabilitation in PD even though it is considered a valid compensatory strategy, for example, for gait impairment in people with PD (Tosserams et al., 2020). Several studies using short-term motor imagery practice in PD cohorts have failed to demonstrate significant motor facilitation (Abraham et al., 2021, Caligiore et al., 2017).

However, in longer-term motor rehabilitation studies, it has been shown that patients with PD have the capacity to perform motor imagery of everyday actions to improve slowness (Tamir et al., 2007) and can increase motor imagery-related brain activation using fMRI-based NF (Subramanian et al., 2016, Tinaz et al., 2018).

Different types of motor imagery (e.g., kinesthetic vs visual) recruit different brain circuits (e.g., sensorimotor vs. visuospatial) (Guillot et al., 2009, Guillot et al., 2014, Lorey et al., 2009). Here, we specifically chose kinesthetic motor imagery training as the experimental intervention because this type of imagery preferentially recruits sensorimotor brain regions including the insula (Lorey et al., 2009, Marchesotti et al., 2016). The imagery training focused on gross motor tasks that belonged to the individual’s motor repertoire and activities of daily living (e.g., walking, gym exercises, shoveling snow). We used visual imagery training without NF as the active control condition. We hypothesized that kinesthetic motor imagery training with NF, but not visual imagery, would increase the task-based and resting-state right insula-dmFC functional connectivity strength and improve motor performance in patients with PD.

2. Methods

2.1. Overview

This study was designed as a phase 1 clinical trial (NCT03623386). Our design is in line with the recommendations in the “consensus on the reporting and experimental design of clinical and cognitive-behavioral neurofeedback studies (CRED-nf)” checklist (Ros et al., 2020). The sample size was determined based on moderate clinically meaningful improvement criteria as determined by the part III motor exam scores of the Movement Disorders Society-Unified Parkinson’s Disease Rating Scale (MDS-UPDRS) (Goetz et al., 2008). Two groups of subjects with PD were randomly assigned in parallel. The experimental group (MI-NF) received NF-guided kinesthetic motor imagery training and the active control group (VI) received visual imagery training without NF. The control group was not aware of the NF component in the experimental arm. Subjects in both groups were provided with explicit imagery strategies and practiced their respective imagery tasks daily for about four weeks. The primary clinical outcome assessor (A.P.) was masked to the intervention assignment.

Our choice of the control condition warrants further explanation: Our preliminary findings clearly suggested that frustration and mental disengagement during the task consistently generated negative NF in subjects with PD (Tinaz et al., 2018). Sham (i.e., noncontingent) NF is a commonly used control condition in NF paradigms. However, it can induce frustration and disengagement (Sorger et al., 2019, Stoeckel et al., 2014), hence, generate negative NF which might lead to overestimation of the effects of contingent NF. Therefore, we chose an active control condition instead of using sham NF to minimize this risk. We used visual imagery as the control task to account for the cognitive and attentional aspects of the experimental kinesthetic motor imagery task. This design also allowed us to match the duration of imagery practice between the groups throughout the study, and to facilitate motivation and active participation of the control group.

We used a bar graph representation of the functional connectivity strength between the right insula and dmFC as the NF modality. The measure of successful NF learning was the increase in this functional connectivity strength from the first (baseline) to the last (post-training) imagery scans that were performed without NF by all subjects. We used intermittent NF which has the following advantages over continuous NF: 1) avoiding potentially distracting confound of cognitive load associated with constant NF monitoring (Stoeckel et al., 2014), and 2) separating out brain activity related to NF evaluation and actual task performance (Johnson et al., 2012).

2.2. Subjects

We recruited subjects with PD defined according to the UK Brain Bank diagnostic criteria (Hughes et al., 1992) through the Yale Movement Disorders Clinic and via local PD support groups in Connecticut. All subjects participated in the study after giving written informed consent in accordance with the procedures approved by the Human Research Protection Office of the Yale School of Medicine. We conducted the study at the Yale Magnetic Resonance Research Center. Subjects underwent screening for MRI safety and were asked detailed questions about their clinical history to determine eligibility. Subjects with PD who were not fully independent, had a neurological or psychiatric disorder (other than PD and comorbid depression or anxiety), or a medical condition that might affect the central nervous system, history of alcohol or illicit drug abuse, head injury resulting in loss of consciousness, dementia (Montreal Cognitive Assessment (MoCA) score < 21), or contraindications for MRI were excluded (see Supplementary Material for details).

We enrolled the subjects continuously and randomized them to two groups using the “minimization” method (Zhao et al., 2015). After truly randomly allocating the first ten subjects, each subsequent subject was allocated to either group considering the factors age, gender, and disease stage. This method allowed us to minimize the imbalance between the groups in terms of demographic and clinical characteristics.

2.3. Sample size calculation

The sample size calculation was based on the clinically important difference (CID) in the primary motor outcome measure. Approximately a 5-point decrease in the mean MDS-UPDRS motor exam score after an intervention has been defined as a moderate CID (Schrag et al., 2006, Shulman et al., 2010). In this study, we expected a moderate CID. We used the mean and standard deviation of the MDS-UPDRS part III motor exam scores (24 ± 9) of a large group of subjects with PD (N = 347) as the population mean and standard deviation (Shulman et al., 2010). Assuming α = 0.05 and power = 0.80, 22 subjects were required for a moderate CID. Considering the possibility of a 15% drop-out rate, we planned to recruit 25 subjects for each group, a total of 50 subjects. We recruited 54 subjects (i.e., the exclusion rate was 18.5%). A total of 10 subjects were excluded: Two dropped out, one could not complete the visits due to Covid-related restrictions, one felt claustrophobic during the first MRI scan, one was found to have a cavernoma in the occipital cortex, and five (two in the VI group and three in the MI-NF group) were excluded due to homework noncompliance (i.e., < 50% completion rate of daily homework).

2.4. Clinical and psychometric evaluations

Subjects underwent all motor assessments in the medication “off” state (12-hr washout) because the “off” state reflects the disease severity more accurately. We assessed disease severity and stage using the MDS-UPDRS (Goetz et al., 2008) and the Hoehn and Yahr (H & Y) scale (Hoehn and Yahr, 1967). We videotaped the part III motor exam of the MDS-UPDRS for all subjects. A movement disorders neurologist (A.P.), who was masked to the group assignment and visit order of the subjects, viewed the videos and scored the exam features (except for rigidity, which was scored by S.T. during the actual exam).

The motor function tests included the two-minute endurance walking, timed up-and-go, five times sit-to-stand, 360-degree turning, and Physical Performance Test to measure the speed of gross movements. The Physical Performance Test measures the speed of everyday motor tasks (e.g., writing, simulated eating, putting on and removing a coat, climbing up the stairs, etc.) (Reuben et al., 1990).

To rule out dementia we administered the MoCA test (Nasreddine et al., 2005). In line with the MDS Task Force recommendations, we administered the Spielberger State-Trait Anxiety Inventory (STAI-S and -T) (Spielberger et al., 1983), Beck Depression Inventory-II (BDI-II) (Beck et al., 1996), and Starkstein Apathy Scale (Starkstein et al., 1992) to evaluate subjects for anxiety, depression, and apathy, respectively. In addition, we administered the Parkinson’s Fatigue Scale (PFS-16) (Brown et al., 2005) and PD Quality of Life Questionnaire (PDQ39) (Jenkinson et al., 1997). We used the 12-item Movement Imagery Questionnaire-3 (MIQ-3) (Williams et al., 2012) to assess subjects’ ability for visual and kinesthetic imagery of movements. For the VI group, we used the Scales for Outcomes of PD-cognition (Marinus et al., 2003) and additional surveys (proactive coping inventory (Greenglass & Schwarzer, 1998); self-efficacy (Chen et al., 2001); self-regulation (Schwarzer et al., 1999)) as “fillers” in visit 3.

2.5. MRI Scanning

We scanned the subjects in the morning or early afternoon after they took their dopaminergic medication. We collected the scans of 32 subjects in a 3.0 Tesla Siemens Trio TIM, and those of 12 subjects in a 3.0 Tesla Siemens Prisma scanner using a 32-channel head coil and identical sequences in both scanners. First, we collected high-resolution T1-weighted MPRAGE images (176 slices, voxel size: 1 mm3, FoV: 250 mm, matrix: 256 × 256, TR: 1900 ms, TE: 2.52 ms, TI: 900 ms, flip angle: 9 degrees) for an accurate localization of the functional MRI data in the beginning of each scan session. Then, we collected T1-weighted FLASH axial images (36 slices, voxel size: 0.9 × 0.9 × 4 mm, FoV: 224 mm, matrix: 256 × 256, TR: 300 ms, TE: 2.47 ms, flip angle: 60 degrees) as an intermediate scan to coregister MPRAGE and echo-planar functional images for the imagery scans. Finally, we obtained axial oblique T2*-weighted, echo-planar functional images (36 slices, voxel size: 3.5 × 3.5 × 4 mm, FoV: 224 mm, matrix: 64×64, TR: 2000 ms, TE: 25 ms, flip angle: 90 degrees). The number of acquisitions was 120 (4 min) for all imagery scans and 304 (10 min 8 s) for the resting-state scans. We always collected the resting-state scans first, during which subjects were instructed to keep their eyes closed, avoid any voluntary movement, and let their mind wander. We assessed wakefulness at the end of the scan by subjects’ report.

2.6. Study Protocol

The study protocol is summarized in Fig. 1. After screening, subjects were scheduled for an initial evaluation visit. Visit 2 was scheduled within five weeks after the initial visit. The subsequent visits were on average two weeks apart.

Fig. 1.
  • Download figure
  • Open in new tab
Fig. 1. Study Protocol

Visit 1: We performed the initial clinical and psychometric evaluations including the administration of the MDS-UPDRS, motor function tests, MoCA, BDI-II, STAI, Apathy Scale, PFS-16, PDQ39, and MIQ-3 in this visit. Upon completion of visit 1, we randomized the eligible subjects to either the VI or MI-NF group.

Visit 2: First, we reviewed the basic concepts and purpose of the respective imagery tasks with all subjects. We discussed the potential benefits of imagery practice in improving motor functioning and visuospatial attention with the MI-NF and VI groups, respectively, and emphasized the importance of consistent practice to achieve these benefits. We also provided the respective imagery strategies in written form.

MI-NF group: We determined each subject’s motor repertoire, identified their motor difficulties, and familiarized them with the kinesthetic motor imagery tasks. Subjects were primed for body awareness by engaging in a mindfulness body scan practice during which they listened to an audio recording guiding them to pay attention to sensations in different body parts. Subsequently, subjects practiced kinesthetic motor imagery of basic movements via audio-recorded scripts. Subjects first performed the movements (e.g., raise your knee, lift up your arm, tap your foot), then imagined performing the movements. Finally, subjects practiced kinesthetic motor imagery of complex whole-body movements of common activities (e.g., walking, balance exercises, calisthenics) via audio-recorded scripts and were instructed to focus on body sensations and positively-valenced emotions associated with imagined movements.

VI group: Subjects received guided visual imagery training via audio-recorded scripts. Subjects practiced visual imagery of static objects or scenery (e.g., tree, lake, sunset, mountain, special place) from a motionless first-person view focusing on sensory features such as colors, shapes, smells, sounds. There was no reference to body movements or sensations in the scripts and subjects were instructed to avoid movement in their imagery. Subjects were encouraged to evoke memories of familiar places or scenes during imagery and focus on positively-valenced emotions associated with the visual imagery content.

The MRI scans followed the practice session. All subjects completed the resting-state fMRI scan first. Subsequently, all subjects practiced their respective baseline imagery task in the scanner during which no NF was provided. Subjects in the MI-NF group then started the NF-guided kinesthetic motor imagery sessions in the scanner. They were free to choose the content of imagery, but told not to change it within a block. There were 4-5 NF training runs each lasting 4 min for the MI-NF group in Visit 2. The VI group did not receive any NF training sessions. All subjects were debriefed about their experience and imagery content after scanning. The MI-NF group was encouraged to use the “winning” strategies in NF sessions during homework practice.

Homework: Upon completion of scanning, we instructed subjects on how to perform and report their daily imagery homework using the Yale Qualtrics online survey platform. The MI-NF group was instructed to practice mindfulness body scan and kinesthetic motor imagery via audio-recorded scripts at home every day for two weeks until Visit 3. We structured the online homework survey based on the key components of motor imagery practice including the setting, type of activity, vividness, difficulty, duration, as well as the type and quality of body sensations and emotions evoked by the imagined movements (Collins and Carson, 2017). Subjects in the VI group were also instructed to continue practicing visual imagery exercises via audio-recorded scripts at home every day for two weeks until Visit 3. The online homework survey included questions about the setting, content, sensory modalities, vividness, difficulty, and duration of the visual imagery practice, as well as the emotions evoked by this practice. All subjects responded to the survey questions by selecting existing options for each respective imagery component and also entered free text. All subjects received daily reminders via emails to ensure adherence to daily imagery exercises.

Visit 3: First, we reviewed the homework of all subjects and discussed strategies for improvement and refinement of their imagery practice as needed. Then, the MI-NF group completed another set of 6-7 NF training runs in Visit 3. The VI group completed psychometric tests and surveys, which served as “fillers” to match the exposure time.

Both groups continued their respective imagery practices at home and completed the online homework surveys until Visit 4.

Visit 4: Subjects held the morning dose of their dopaminergic medications. The MDS-UPDRS part III motor exam, motor function tests, and MIQ-3 were repeated on all subjects. We discussed with all subjects in the MI-NF group whether they noticed any improvement in the specific motor difficulties that they reported in the first visit. After these assessments, subjects took their medication and completed a resting-state scan as well as an imagery scan during which they practiced kinesthetic motor and visual imagery, respectively, without NF.

2.7. Outcome measures

The primary imaging outcome measure was the change in right insula-dmFC functional connectivity strength 1) from the first to the last respective imagery scans without NF and 2) from the first to the last resting-state scans. We hypothesized that the MI-NF group would show a significant increase in both functional connectivity measures compared with the VI group.

The primary clinical outcome measure was the change in MDS-UPDRS part III motor exam scores from the first to the last visits. The secondary clinical outcome measure was the change in motor function test scores from the first to the last visits. We hypothesized that the MI-NF group would show a significant improvement in both clinical measures compared with the VI group.

We also examined the relationship between the imaging and clinical outcome measures in each group. We hypothesized that the increase in the right insula-dmFC functional connectivity would correlate with improvement in MDS-UDRS part III and motor function test scores only in the MI-NF group.

In addition, change in MIQ-3 scores and change in whole-brain resting-state and imagery task-based functional connectivity without NF were exploratory outcome measures. We hypothesized improvement in MIQ-3 scores in both groups and imagery task-specific changes in functional connectivity in both groups in the motor and visual networks, respectively.

2.8. Clinical and psychometric data analysis

We first assessed the normality of distribution of all clinical scores using the Shapiro-Wilk test. We compared the means and standard deviations of the normally distributed, and medians and median absolute deviations of the non-normally distributed scores with the population means or cutoff scores of the respective tests, when applicable, using one-sample t-tests (p < 0.05, two-tailed). For between-group comparisons, we used the two-sample t-test for normally distributed and nonparametric Mann-Whitney U test for non-normally distributed continuous variables (p < 0.05, two-tailed). We used the Chi-Square test for a between-group comparison of the categorical variables including handedness and disease onset side (p < 0.05, two-tailed). We used repeated measures ANOVA (within-subject: time with two levels, i.e., pre- and post-training; between-subject: group) to evaluate the between-group changes in the MDS-UPDRS part III exam scores, motor function test performances, and MIQ-3 scores as a result of training (p < 0.05, two-tailed). We used the SPSS 26 software for all statistical analyses.

2.9. Homework analysis

The number of entries and total time investment in homework, vividness, and difficulty of imagery were numerical values and averaged. Vividness and difficulty of imagery were rated on a Likert-like scale from 0 to 10, 0: not vivid/not difficult and 10: very vivid/very difficult. We compiled the selected and free-text responses in each imagery component under categories and calculated the percentages.

2.10. Imaging Data Analysis

2.10.1. Right insula – dmFC functional connectivity during motor imagery scans with NF

A subdued background image of a night sky instructed subjects to engage in kinesthetic motor imagery for 40-s blocks (total of five blocks). At the end of each kinesthetic motor imagery block, NF was plotted as a bar graph (blue: negative, red: positive feedback) and presented for 8 s. NF was defined as the functional connectivity between the right insula and dmFC. We created cubic anatomical masks (6 × 6 × 6 mm) centered at x = 44, y = 4, z = 8 in the right insula and at x = −4, y = 2, z = 62 in the dmFC (Fig. 2). These regions of interest were functionally determined using a silent heartbeat counting task as described before (Tinaz et al, 2018). These masks were created based on the standard Montreal Neurological Institute (MNI) brain template and then registered to each subject’s native functional space using a series of transformations. All transformations were estimated using the Bioimage Suite software and manually inspected for accuracy. Functional scans of each subject were motion-corrected in real-time using the algorithms described in Scheinost et al. (2013). Covariates of no interest were regressed out including linear and quadratic drifts; mean cerebral spinal fluid, white matter, and gray matter signals; and motion-related confounds (24-parameter motion model including six rigid-body motion parameters, six temporal derivatives, and these terms squared). The data were temporally smoothed with a zero-mean unit variance Gaussian filter (approximate cutoff frequency = 0.12 Hz). Then, the signal time course of the right insula and dmFC masks in a given subject were computed as the average time course across all voxels within each of these masks. Finally, the right insula-dmFC functional connectivity was obtained by correlating the time courses and Fisher z-transforming the r-values (Scheinost et al., 2020). A Matlab program plotted the z-values as a bar graph, which was presented to the subjects in the MI-NF group as NF.

Fig. 2.
  • Download figure
  • Open in new tab
Fig. 2. Neurofeedback brain-computer interface.

Subjects in the MI-NF group performed motor imagery in the scanner during task blocks and received NF in the form of red (positive) and blue (negative) bars. The height of the bars reflected the strength of the positive/negative right insula-dmFC functional connectivity (fc). During the first and last imagery sessions without NF, all subjects performed their respective imagery in the scanner during task blocks and received no NF (flat line). Subjects in the VI group never received NF.

2.10.2. Right insula – dmFC functional connectivity during motor and visual imagery scans without NF

We used the same task design and followed the same steps to compute the right insula-dmFC functional connectivity during imagery scans without NF as described in section 2.10.1. However, instead of bar graphs representing the z-values, a flat line (i.e., no NF) was presented. Subjects were explicitly told not to expect NF during these scans.

We calculated the percent change in the right insula-dmFC functional connectivity from the first to the last respective imagery scans without NF and performed a between-group comparison using the nonparametric Mann Whitney U test (p < 0.05, two-tailed).

2.10.3. Analysis of the resting-state fMRI scans

We used the Connectivity toolbox v17 for the resting-state data analysis (Whitfield-Gabrieli and Nieto-Castanon, 2012). Preprocessing steps included the removal of the first four scans to reach magnetization steady state, motion correction, outlier detection (frame-wise displacement above 0.9 mm or global signal changes above 5 standard deviations), coregistration of functional scans with the structural scan, normalization to the standard MNI brain template, and smoothing with a Gaussian kernel with a FWHM of 8 mm to account for inter-individual anatomical variability. De-noising steps included correction for physiological and other sources of noise by regressing out the principal components of the white matter and cerebrospinal fluid signal using the CompCor method (Chai et al., 2012), regression of motion artifacts and outliers before bandpass-filtering, and linear detrending. Global signal was not removed. Finally, we bandpass-filtered (0.008 < f < 0.09 Hz) the data to capture the fluctuations of the blood oxygenation level-dependent (BOLD) signal that typically occur within this frequency range at rest.

For the whole-brain functional connectivity analyses, we used the functionally defined nodes (n = 268) of the whole-brain Shen atlas (Shen et al., 2013). For each subject, we extracted the average BOLD signal time courses from these nodes and correlated them with each other using Pearson correlations. The r-values corresponded to the functional connectivity strength between node pairs. We Fisher z-transformed the r-values and obtained group-level functional connectivity maps for statistical analyses. We used a 2 × 2 mixed ANOVA (within-subject: time with two levels, i.e., pre- and post-training; between-subject: group) for a between-group comparison of the training effect. We used the false discovery rate (FDR) method for correction for multiple comparisons (p < 0.05, two-tailed) (Genovese et al., 2002).

For the right insula-dmFC functional connectivity analyses, we extracted and correlated the resting-state BOLD signal time courses from these regions as explained in the previous paragraph. We then calculated the percent change in the right insula-dmFC functional connectivity from the first to the last resting-state scan and performed a between-group comparison using the nonparametric Mann Whitney U test (p < 0.05, two-tailed).

2.10.4. Whole-brain connectivity analyses of imagery scans without NF

We applied the same preprocessing and denoising steps to the imagery scans as explained in section 2.10.3 with the following differences: The first three scans were removed and the BOLD signal was high-pass filtered [0.008 < f < Inf]. We then used the generalized psychophysiological interaction (gPPI) model in the Connectivity toolbox v17 (Whitfield-Gabrieli and Nieto-Castanon, 2012) to assess the imagery task-based whole-brain functional connectivity changes in each group. We convolved the imagery blocks and feedback periods separately with the canonical hemodynamic response function. We used the same functional connectivity calculations based on the imagery blocks and the 2 × 2 mixed ANOVA (within-subject: time with two levels, i.e., pre- and post-training; between-subject: group) as explained in section 2.10.3 for between-group comparison of the training effect (FDR-corrected p < 0.05, two-tailed).

Of note, to address the potential confounding effects of different scanners, we compared the voxel-wise global mean correlations (GCOR) for resting-state and imagery scans after denoising between subjects scanned in the Trio TIM and Prisma scanners.

2.10.5. Correlations between imaging and clinical outcomes

We calculated the average percent change across all motor function test scores and MDS-UPDRS part III scores for each group. We then performed Kendall’s tau-b nonparametric bivariate correlations with bootstrapping (x1000) between the average percent change in motor function and part III scores and percent change in the 1) imagery task-based and 2) resting-state right insula-dmFC functional connectivity separately for each group (p < 0.05, two-tailed).

3. Results

3.1. Clinical and psychometric data at baseline

There was no significant difference between the MI-NF and VI groups in any of the demographic, clinical, and psychometric variables listed in Table 1. Moreover, both groups scored significantly below the anxiety, apathy, minimal depression, and fatigue cutoff scores; significantly above the MoCA cutoff score for mild cognitive impairment, and had a significantly better quality of life compared with the PD population with the same disease stage of H & Y 2 (see Supplementary Material Section 2 and Tables S1 and S2). The details of the medication types and LEDD distribution are shown in Tables S3 and S4 in the Supplementary Material.

View this table:
  • View inline
  • View popup
Table 1. Clinical and psychometric data

3.2. Training effects on behavior

The repeated measures ANOVA revealed no significant main effect of group, no significant time × group interaction, but significant main effect of time indicating that both groups significantly improved their gross motor performance and motor imagery skills post-training (Physical Performance Test: F(1,42) = 8.967, p = 0.005, endurance walking: F(1,42) = 18.276, p = 0.000, gross motor combined: F(1,42) = 5.309, p = 0.026; MIQ-3: F(1,42) = 6.553, p = 0.014).

However, the repeated measures ANOVA on the MDS-UPDRS part III scores did not show any significant main effect of group (F(1,42) = 1.204, p = 0.279), time F(1,42) = 0.072, p = 0.790), or significant time × group interaction (F(1,42) = 1.085, p = 0.303) (Table 2). A power analysis based on the mean difference and standard deviation between the pre- and post-training MDS-UPDRS part III scores of the MI-NF group revealed an effect size of 0.169 (α = 0.05 and power = 0.118). In other words, for α = 0.05 and power = 0.80, 276 subjects would be needed to see a significant change in the part III scores.

View this table:
  • View inline
  • View popup
Table 2. Pre- and post-training behavioral data

3.3. Characteristics of the imagery homework

The overall engagement and time investment in imagery homework were comparable between the groups (two-sample t-test, p = 0.934) with less variability across subjects in the VI group. Imagery vividness and difficulty ratings were not significantly different between the groups (two-sample t-test, p = 0.406 and 0.899, respectively). The VI group reported significantly higher rates of positive emotions (e.g., relaxed, pleasant) compared with the MI-NF group (e.g., energized, pleasant) (Mann – Whitney U Test, p < 0.001) associated with the imagery practice (Table 3).

View this table:
  • View inline
  • View popup
  • Download powerpoint
Table 3. Imagery homework

Walking, aerobic and resistance exercises, and house chores were common themes in kinesthetic motor imagery practices. Special indoor places and landscapes with mountains, lakes, etc. were common themes in visual imagery practices. Sensations of movement, heartbeat, and breathing were the most commonly reported body sensations associated with kinesthetic motor imagery. Light, color, and shape were the most commonly reported sensory modalities associated with visual imagery (see Supplementary Material Tables S5 and S6 for details).

3.4. Debriefing at the last visit

All subjects in the MI-NF group reported enhanced body and kinesthetic awareness as well as more presence and mindfulness during everyday movements as a result of kinesthetic motor imagery practice. One subject described this practice as “body amnesia prevention.” Increased awareness of asymmetry between the most- and least-affected sides in movement speed and range of motion was reported most frequently. This awareness prompted the subjects to make conscious corrective efforts and led to improvement in limb symmetry and coordination in movements (e.g., improved arm swing and stride length while walking, improved hook and punch while boxing, more symmetrical strokes while swimming, more coordinated upper and lower limb movements during gym exercises). About one-third of the subjects also reported improved balance while walking and putting pants or socks on. One subject reported remarkable improvement in freezing episodes while walking and jogging and was pleased by his ability to complete a road race. On the other hand, there was no noticeable improvement in fine motor skills (e.g., writing, typing, shaving, brushing teeth, playing the guitar or piano) or tremor.

3.5. Imaging

3.5.1. Preprocessing

On average, subjects in both groups had submillimeter head motion in both resting-state and imagery scans without NF. Head motion did not differ significantly between the groups. The GCOR distribution was also not significantly different between subjects who were tested in different scanners suggesting that denoising successfully removed any potential scanner effects on the voxel-wise correlations (see Supplementary Material Tables S7 and S8).

3.5.2. Task-based and resting-state functional connectivity changes

The NF sessions with short breaks in-between to reduce fatigue were overall well tolerated. Negative NF was usually associated with distraction, fatigue, thinking of the movements rather than imagining them, weak body sensations, and frustration due to imagining limitations and negative body sensations during actual movements. Positive NF was usually associated with vivid imagery with strong and positive body sensations, sustained focus, feeling in a flow state, and positive emotions (e.g., sense of accomplishment). At times, subjects felt that the NF was noncontingent with their level of effort or the same imagery practice yielded positive NF in one block and negative NF in another.

The bar graphs in Fig. 3 display the Fisher z-transformed right insula-dmFC functional connectivity strength (mean ± standard error) during rest and imagery scans with and without NF and the percent change in the right insula-dmFC functional connectivity strength (mean ± standard error) from the first to the last rest and imagery scans without NF for both groups. There was significant between-subject variability in all functional connectivity measures in both groups. Between-group comparisons of the percent change in right insula-dmFC functional connectivity strength from the first to the last rest and imagery scans without NF did not reveal a significant difference (Mann-Whitney U test, p = 0.481 and 0.981, respectively).

Fig. 3.
  • Download figure
  • Open in new tab
Fig. 3. Right insula-dmFC functional connectivity.

Error bars show the standard error. FC: Functional connectivity, Im1 and Im2: first and last respective imagery scans without neurofeedback, NF1 and NF2: Neurofeedback scans on days 1 (visit 2) and 2 (visit 3), rs1 and rs2: first and last resting-state scans. MI-NF: Im1: 0.05 ± 0.08, NF1: −0.01 ± 0.05, NF2: 0.01 ± 0.06, Im2: 0.09 ± 0.06, real-time % change: 93.2 ± 84.1, rs1: 0.22 ± 0.05, rs2: 0.23 ± 0.07, resting-state % change: −2.9 ± 43.7. VI: Im1: −0.07 ± 1.00, Im2: 0.03 ± 0.07, real-time % change: 97.8 ± 106.6, rs1: 0.12 ± 0.05, rs2: 0.17 ± 0.06, resting-state % change: 74.2 ± 58.8.

3.5.3. Correlations between functional connectivity and motor performance changes

The Kendall’s tau-b correlation between the average percent change in motor performance scores and percent change in imagery task-based right insula-dmFC functional connectivity was significant in the MI-NF group (correlation coefficient = 0.307, p = 0.045, 95% CI [0.012, 0.553]) (Fig. S1) and not significant in the VI group (correlation coefficient = −0.108, p = 0.481, 95% CI [−0.392, 0.194]).

The Kendall’s tau-b correlation between the average percent change in motor performance scores and percent change in resting-state right insula-dmFC functional connectivity was not significant in the MI-NF (correlation coefficient = 0.273, p = 0.076, 95% CI [−0.26, 0.526]) or the VI (correlation coefficient = −0.195, p = 0.204, 95% CI [−0.464, 0.018]) group.

The Kendall’s tau-b correlation between the average percent change in MDS-UPDRS part III scores and percent change in imagery task-based right insula-dmFC functional connectivity was not significant in the MI-NF group (correlation coefficient = −0.035, p = 0.821, 95% CI [−0.328, 0.264]) or the VI group (correlation coefficient = −0.091, p = 0.554, 95% CI [−0.377, 0.211]).

The Kendall’s tau-b correlation between the average percent change in MDS-UPDRS part III scores and percent change in resting-state right insula-dmFC functional connectivity was not significant in the MI-NF (correlation coefficient = 0.017, p = 0.910, 95% CI [−0.280, 0.312]) or the VI (correlation coefficient = −0.091, p = 0.554, 95% CI [−0.377, 0.211]) group.

3.5.4. Whole-brain resting-state functional connectivity

The VI group did not show any significant pairwise whole-brain functional connectivity differences in the rest 2 > rest 1 contrast. The MI-NF group showed significantly reduced pairwise whole-brain functional connectivity in the rest 2 > rest1 contrast between several ventral and dorsal frontal nodes, and between the temporal pole and caudate nodes, and increased connectivity between a few cerebellar and fronto-temporal nodes. There was a significant group × time interaction (MI-NF > VI and rest 2 > rest 1 contrast) demonstrating decreased functional connectivity between only a few cerebellar nodes and increased functional connectivity between a cerebellar node and orbitofrontal nodes (see Supplementary Material Fig. S2 and Table S9).

3.5.5. Whole-brain imagery task-based functional connectivity

Both groups demonstrated significant pairwise functional connectivity changes as a result of training (imagery 2 > imagery 1 contrast) predominantly in the form of decreased functional connectivity. The main effects of group are summarized in the Supplementary Material Fig. S3 and S4 and Tables S10 and S11. Here, we focus on the interaction effects, namely the MI-NF > VI and imagery 2 > imagery 1 contrast (Fig. 4 and Table 4). As a result of kinesthetic motor imagery training, the MI-NF group showed stronger functional connectivity in cerebellar, occipito-temporal, premotor, and frontal nodes, as well as the insula and midbrain. As a result of visual imagery training, the VI group showed stronger functional connectivity in a different set of cerebellar, occipito-temporal, and frontal nodes, as well as posterior cingulate nodes and pons.

Fig. 4.
  • Download figure
  • Open in new tab
Fig. 4. Whole-brain imagery task-based functional connectivity changes.

Red: MI-NF > VI and blue: VI > MI-NF. Cb: Cerebellum, dPCC: Dorsal posterior cingulate cortex, FG: Fusiform gyrus, FP: Frontal pole, IFGo: Inferior frontal gyrus, orbital part; Ins: Insula, ITG: Inferior temporal gyrus, MTG: Middle temporal gyrus, OFC: Orbitofrontal cortex, PHG: Parahippocampal gyrus, PMC: Premotor cortex, V2: Secondary visual cortex, VA: Visual association cortex.

View this table:
  • View inline
  • View popup
Table 4.

4. Discussion

In summary, the primary imaging and clinical outcomes of our study did not reach statistical significance. However, we did see significant improvements in the secondary clinical outcome where both groups receiving motor imagery (MI-NF) and visual imagery (VI) improved their composite gross motor function scores comparably. Notably, there was a significant positive relationship between the change in real-time right insula-dmFC functional connectivity and improvement in gross motor function scores only in the MI-NF group. There was no significant relationship between other imaging and clinical outcomes. The exploratory outcome measures showed statistical significance. Both groups reported significant subjective improvement in their imagery skills as hypothesized. The whole-brain task-based functional connectivity patterns mapped onto distinct networks reflecting the specific type of imagery training (i.e., kinesthetic motor vs visual) in each group.

4.1. Imaging outcomes: Heterogeneity in NF regulation success

Contrary to our findings based on eight subjects in the pilot study (Tinaz et al., 2018), we did not find significant NF regulation in the MI-NF group. A high degree of variability in NF response in healthy controls and clinical populations is not uncommon. It is estimated that about 30% of participants fail to learn how to regulate their brain signal with NF (Sorger et al., 2019), and studies report mixed results regarding overall regulation success (about 30% of studies) and behavior change (about 50% of studies) (Thibault et al., 2018). There are numerous methodological and subjective factors contributing to this variability (Sorger et al., 2019). During post-scan debriefing, we identified nonspecific factors including mental fatigue and distraction, and subject-specific factors such as frustrated efforts during imagery as potential contributors to NF regulation failure. This is not surprising given that the insula is an essential hub in the appraisal of emotional valence (Craig 2009), and negative emotions during motor imagery would be expected to interfere with the regulation of insula activity and co-activity with the dmFC.

Another potential factor was the choice of mental strategy during NF training. We instructed the subjects to use an explicit kinesthetic motor imagery strategy emulating everyday complex motor activities with the aim of maximizing clinical benefit. However, using spontaneous mental strategies as opposed to following explicit instructions can be more efficient to achieve successful NF modulation (Kober et al., 2013). For example, NF training with explicit motor imagery instructions has been shown to be ineffective in up-regulation of the SMA activity in a small group of healthy controls (Sepulvada et al., 2016). Furthermore, in a graded NF paradigm combined with kinesthetic hand motor imagery, subjects could up-regulate the SMA activity to the specified level, but no interaction between target level and feedback condition was found suggesting a general motor imagery effect, but no specific feedback modulation in the SMA (Mehler et al., 2019).

Finally, the number of NF training runs (a total of 10-12 per subject) may not have been sufficient for all subjects to ensure successful regulation.

4.2. Clinical outcomes: Scale sensitivity and training-specific vs nonspecific effects

There was no clinically important difference in the MDS-UPDRS part III scores in the MI-NF group post-training. The total MDS-UPDRS part III score is the gold standard primary outcome measure in clinical trials in PD. Tremor and rigidity (muscle stiffness) are major components of the scale, which were not specific targets of the motor imagery training. Bradykinesia (slowness) is another major component, also a training target, but is rated using rather simple repetitive movements. Therefore, part III may not have been sensitive enough to the changes in complex movements. These changes were captured by the motor function tests that we used to measure the speed of complex movements and activities of daily living. Both groups improved their motor function scores significantly and comparably suggesting a nonspecific or placebo effect as a result of participating in an interventional study. This is particularly relevant for individuals with PD, who also demonstrate strong placebo effects in clinical trials with pharmacological and surgical interventions (Lidstone 2014). Our study design did not allow us to exclude the placebo effect, however, we think that the motor improvement in the MI-NF group may also reflect some specificity. The modestly significant positive correlation between the change in real-time right insula-dmFC functional connectivity and improvement in motor function scores was observed only in the MI-NF group suggesting training-specific effects. An alternative explanation for the motor function improvement in the VI group would be “spillover” effects. The visual imagery practice required sustained attention, multisensory integration, and visuospatial construction (Pearson 2019). Moreover, subjects in the VI group reported significantly higher rates of positive emotions (predominantly “relaxed”) associated with their daily practice compared with the MI-NF group. Taken together, the visual imagery practice may have indirectly enhanced motor functioning by promoting relaxed alertness and sustained visuospatial attention.

4.3. Exploratory outcomes: Training effects on behavior and brain connectivity

Though the MIQ-3 is a movement imagery questionnaire, the imagery task includes two visualization conditions from a first- and third-person perspective, and a kinesthetic condition. Therefore, it is not surprising that both groups reported significant improvement in the MIQ-3 post-training. The MI-NF group also reported enhanced kinesthetic awareness and mindfulness during everyday movements as a result of practice, which led to corrective efforts and improvement in movement size, speed, coordination, and balance. These were subjective accounts gathered at the post-training debriefing, therefore, could only be analyzed qualitatively. Many NF studies with neuropsychiatric populations use self-evaluation surveys as outcome measures of symptomatic improvement underlining the need to develop specific surveys to assess changes in kinesthetic awareness of complex movements quantitatively.

Notably, we did see specificity by the MI-NF versus the VI training (imagery training effects: imagery 2 > imagery 1 contrast) mapped onto separate whole-brain functional connectivity patterns during the respective imagery scans without NF. Of note, the whole-brain intrinsic functional connectivity changes (rest 2 > rest 1 contrast) were not robust and did not express a clear pattern suggesting that the training effects were specific to task performance. Post-training, the MI-NF group compared with the VI group exhibited stronger functional connectivity between the nodes involved in kinesthetic motor imagery including the insula and premotor cortex, whereas the VI group compared with the MI-NF group demonstrated stronger functional connectivity between the nodes involved in visual imagery including nodes in the ventral and dorsal visual streams. These results provide a manipulation check and validate that the training targeted the expected brain networks, and are also in line with previous reports of differential activation patterns during kinesthetic motor (Guillot et al., 2009; Hétu et al., 2013) and visual imagery (Pearson 2019) corroborating that subjects in both groups adhered to the specific imagery instructions in their practice.

Different parts of the cerebellum show intrinsic functional connectivity with different large-scale brain networks and are involved in a range of sensorimotor and cognitive tasks (Buckner et al., 2011). Different cerebellar nodes were also recruited in both imagery tasks demonstrating increased functional connectivity with the premotor and frontal nodes in the motor imagery task, and with the orbitofrontal (OFC) and dorsal posterior cingulate cortex (dPCC) in the visual imagery task. The OFC is considered a hub for integrating multimodal sensory information with hedonic value (Kringelbach 2005) and may have been involved in the retrieval of multimodal sensory information and integrating it with subjects’ emotional experience during visual imagery. The dPCC has rich anatomical connections with temporal, parietal, frontal regions and is thought to regulate attention according to cognitive demands (Leech and Sharp, 2014). During visual imagery, the dPCC may have been involved in monitoring the retrieval of episodic memories and visuospatial construction of scenes.

In addition to the insula and premotor nodes, increased functional connectivity in the motor imagery network also included the fusiform and parahippocampal nodes, which may have played a role in imagining one’s body movements (Olsson et al., 2008) in specific environmental contexts (Epstein et al., 2017), respectively. Finally, a brainstem node overlapping with the midbrain/red nucleus region also displayed increased functional connectivity with the premotor cortex during motor imagery. The midbrain/red nucleus is an important subcortical region in motor control especially of the limbs, motor coordination, and motor preparation through its connections with the motor cortices, basal ganglia, and cerebellum (Basile et al., 2021).

Taken together, our results demonstrate the neural bases of task-specific practice effects over the course of four weeks in the form of functional re-organization of distinct neural circuits supporting kinesthetic motor and visual imagery. However, we cannot rule out the potential large-scale neuromodulatory effects of NF training in the MI-NF group. At least indirectly, subjects’ experience during NF sessions and post-scan debriefing informed their practice strategies at home and may have helped them hone their kinesthetic motor imagery skills with practice.

Limitations

Involuntary muscle contraction as a result of motor imagery can potentially confound the BOLD signal (Thibault et al., 2018). We did not monitor muscle activity in the scanner using surface EMG electrodes. Though we did not observe overt movements while subjects were in the scanner and head motion levels in the scanner across both groups were equivalent, we could not rule out covert muscle contractions. We also did not track eye movements during imagery blocks, thus, could not assess whether subjects in both groups were matched in the amount of eye movement. The training and assessment of outcome measures were limited to a four-week period with no follow-up. Continued symptomatic benefits for weeks after NF have been demonstrated in clinical populations (Rance et al., 2018).

Conclusions

The functional connectivity-based NF regulation was unsuccessful in our cohort with mild PD. The heterogeneity in NF response in healthy controls and particularly in clinical populations is an ongoing challenge in the field, and concerted efforts are made to address the potential sources of heterogeneity (e.g., selecting the appropriate NF targets, dose, mental strategies, study designs, and meaningful outcome measures, and ways to predetermine potential success). Yet, kinesthetic motor imagery practice by itself is a promising intervention in PD. There is also evidence that multisensorial integrative imagery techniques (e.g., dynamic-cognitive imagery) (Abraham et al., 2018) and motor imagery combined with action observation (Caligiore et al., 2017) may improve motor functions in people with PD. Future studies are needed to develop specific motor imagery protocols and determine practice parameters (e.g., dose, duration, outcome measures), which then can be individualized and incorporated into clinical practice as compensatory rehabilitation strategies for motor impairment in people with PD.

CRediT authorship contribution statement

Sule Tinaz: Conceptualization, Methodology, Investigation, Formal Analysis, Visualization, Supervision, Writing – original draft, Funding Acquisition. Serageldin Kamel: Investigation, Formal Analysis, Data Curation, Visualization. Sai S. Aravala: Investigation, Formal Analysis.. Mohamed Elfil: Investigation, Formal Analysis. Ahmed Bayoumi: Investigation, Formal Analysis. Amar Patel: Investigation, Formal Analysis. Dustin Scheinost: Methodology, Software. Rajita Sinha: Conceptualization, Methodology, Supervision. Michelle Hampson: Conceptualization, Methodology, Software, Resources. All authors were involved in the review and editing of the manuscript.

Declaration of Competing Interest

The authors declare no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Data Availability

All data produced in the present study are available upon reasonable request to the authors

Supplementary Material

1. Eligibility

Inclusion Criteria:

  • Subjects with a diagnosis of idiopathic PD defined according to the UK Brain Bank diagnostic criteria

  • Age ≥ 40 years

  • Expected to be on a stable dopaminergic medication regimen throughout the study period

Exclusion Criteria:

  • Age < 40 years

  • Non-English speaking

  • Pregnancy

  • Breastfeeding

  • Excessive alcohol consumption (> 7 drinks per week for women, > 14 drinks per week for men) or illicit substance use

  • History of a neurological disorder such as a brain tumor, stroke, central nervous system infection, multiple sclerosis, movement disorder (other than PD), or seizures

  • History of schizophrenia, bipolar disorder, attention deficit disorder, or obsessive-compulsive disorder

  • History of head injury with loss of consciousness longer than a minute

  • Metallic surgical implants or traumatically implanted metallic foreign bodies

  • Inability to lie flat for about an hour in the MRI scanner

  • Discomfort being in small, enclosed spaces

  • Dementia (Montreal Cognitive Assessment score < 21)

  • Depression (Beck Depression Inventory-II score > 28 indicating more than moderate depression)

  • Hoehn & Yahr stage > 3

  • Focal neurological findings on exam that suggest cerebral pathology other than that associated with parkinsonism

  • Motor symptoms that could potentially introduce too much motion artifact in the imaging data (e.g., Movement Disorders Society-Unified Parkinson’s Disease Rating Scale resting tremor score > 2 in limbs, head/chin tremor, or more than mild dyskinesia by history or exam).

  • Poor homework compliance (< 50% completion rate).

2. Clinical and psychometric results

View this table:
  • View inline
  • View popup
Table S1. Normative mean and cutoff scores for each test

2.1. Comparisons with normative data

MI-NF group: The distribution of the MoCA and apathy scores was not normal. We used the median scores of MoCA and apathy scale, and mean scores of the rest for comparison with normative data:

MoCA: The median score (29) was significantly above the MCI cutoff (p = 0.000).

STAI-S: Males (mean: 30.33) were not significantly different from (p = 0.058) and females (mean: 28.10) showed significantly less state anxiety (p = 0.014) than the age-matched normal population.

STAI-T: Males (mean: 33.83) and females (mean: 32.70) were not significantly different from the age-matched normal population in trait anxiety (p = 0.989 and 0.717, respectively).

BDI-II: The mean score was significantly below the mild depression cutoff (p = 0.000). Apathy: The median score (5) was significantly below the apathy cutoff (p = 0.000).

PFS-16: The mean score was significantly below the fatigue cutoff (p = 0.002).

PDQ39-SI: The mean score was significantly below the cutoff for poor quality of life (p = 0.000).

VI group: The distribution of the MoCA, STAI-S, STAI-T, BDI-II, PFS-16, and PDQ39-SI was not normal. We used the median scores for comparison with normative data:

MoCA: The median score (29) was significantly above the MCI cutoff (p = 0.000).

STAI-S: Males (median: 31.50) were not significantly different (p = 0.192) and females (median: 26.00) showed significantly less state anxiety (p = 0.006) than the age-matched normal population.

STAI-T: Males (median: 34.00) and females (median: 30.50) were not significantly different from the age-matched normal population in trait anxiety (p = 0.953 and 0.389, respectively). BDI-II: The median score (6) was significantly below the mild depression cutoff (p = 0.000). Apathy: The mean score was significantly below the apathy cutoff (p = 0.000).

PFS-16: The median score (2.095) was significantly below the fatigue cutoff (p=0.000).

PDQ-39 SI: The median score (12.32) was significantly below the cutoff for poor quality of life (p=0.000).

2.2. Between-group comparisons

The two-sample t-test on normally distributed variables, Mann-Whitney U test on non-normally distributed variables, and Chi-Square test for handedness and disease onset side did not reveal any significant difference between the MI-NF and VI groups.

View this table:
  • View inline
  • View popup
  • Download powerpoint
Table S2. Between-group statistics p values
View this table:
  • View inline
  • View popup
  • Download powerpoint
Table S3. Medications
View this table:
  • View inline
  • View popup
Table S4. LEDD distribution
View this table:
  • View inline
  • View popup
  • Download powerpoint
Table S5. MI-NF group imagery homework
View this table:
  • View inline
  • View popup
Table S6. VI group imagery homework
View this table:
  • View inline
  • View popup
  • Download powerpoint
Table S7. Head Motion
View this table:
  • View inline
  • View popup
  • Download powerpoint
Table S8. Voxel-wise global mean correlations (GCOR) per scanner after denoising
Fig. S1.
  • Download figure
  • Open in new tab
Fig. S1. Relationship between the imaging and motor outcomes in the MI-NF group.

The graph displays the relationship between the average percent change in motor performance scores and percent change in real-time imagery task-based right insula-dmFC functional connectivity in the MI-NF group.

Fig. S2.
  • Download figure
  • Open in new tab
Fig. S2. Whole-brain resting-state functional connectivity: MI-NF > VI and rest2 > rest1

Red: increased, blue: decreased. Cb: Cerebellum, FG: Fusiform gyrus, IFGo: Inferior frontal gyrus, orbital part.

View this table:
  • View inline
  • View popup
  • Download powerpoint
Table S9. MI-NF > VI resting-state functional connectivity
Fig. S3.
  • Download figure
  • Open in new tab
Fig. S3. Whole-brain imagery task-based functional connectivity main effect of MI-NF (imagery 2 > imagery 1).

Red: increased, blue: decreased. Cb: Cerebellum, Cd: Caudate, dlPFC: Dorsolateral prefrontal cortex, dPCC: Dorsal posterior cingulate cortex, FEF: Frontal eye fields, FP: Frontal pole, HP: Hippocampus, IFGo: Inferior frontal gyrus, orbital part; ITG: Inferior temporal gyrus, MTG: Middle temporal gyrus, PHG: Parahippocampal gyrus, Put: Putamen, RSC: Retrosplenial cortex, SMA: Supplementary motor area, SMG: Supramarginal gyrus, STG: Superior temporal gyrus, TP: Temporal pole, vPCC: Ventral posterior cingulate cortex.

View this table:
  • View inline
  • View popup
Table S10. MI-NF imagery task-based functional connectivity
Fig. S4.
  • Download figure
  • Open in new tab
Fig. S4. Whole-brain imagery task-based functional connectivity main effect of VI (imagery 2 > imagery 1).

Red: Increased, blue: decreased. Amyg: Amygdala, Bs: Brainstem, Cb: Cerebellum, FEF: Frontal eye fields, FP: Frontal pole, IFGo: Inferior frontal gyrus, orbital part; Ins: Insula, MTG: Middle temporal gyrus, OFC: Orbitofrontal cortex, PHG: Parahippocampal gyrus, PMC: Premotor cortex, SMG: Supramarginal gyrus, TP: Temporal pole, vPCC: Ventral posterior cingulate cortex.

View this table:
  • View inline
  • View popup
Table S11. VI imagery task-based functional connectivity

Acknowledgment

This work was supported by the National Center for Advancing Translational Science, a component of the National Institutes of Health (grant numbers UL1TR001863, KL2TR001862) and by the National Institute of Neurological Disorders and Stroke (grant number K23NS099478). We thank Jessica Lemere and Sayan Basu for their help with data collection, and Jessica Lemere for her help with imagery audio recordings.

Footnotes

  • R. Sinha and M. Hampson are both senior authors.

References

  1. ↵
    Abraham A, Duncan RP, Earhart GM. The Role of Mental Imagery in Parkinson’s Disease Rehabilitation. Brain Sci. 2021 Feb 2;11(2):185. doi: 10.3390/brainsci11020185. PMID: 33540883; PMCID: PMC7913152.
    OpenUrlCrossRefPubMed
  2. ↵
    Abraham A, Hart A, Andrade I, Hackney ME. Dynamic Neuro-Cognitive Imagery Improves Mental Imagery Ability, Disease Severity, and Motor and Cognitive Functions in People with Parkinson’s Disease. Neural Plast. 2018 Mar 14;2018:6168507. doi: 10.1155/2018/6168507. PMID: 29725348; PMCID: PMC5872663.
    OpenUrlCrossRefPubMed
  3. ↵
    Baraduc P, Thobois S, Gan J, Broussolle E, Desmurget M. A common optimization principle for motor execution in healthy subjects and parkinsonian patients. J Neurosci. 2013 Jan 9;33(2):665–77. doi: 10.1523/JNEUROSCI.1482-12.2013. PMID: 23303945.
    OpenUrlAbstract/FREE Full Text
  4. ↵
    Basile GA, Quartu M, Bertino S, Serra MP, Boi M, Bramanti A, Anastasi GP, Milardi D, Cacciola A. Red nucleus structure and function: from anatomy to clinical neurosciences. Brain Struct Funct. 2021 Jan;226(1):69–91. doi: 10.1007/s00429-020-02171-x. Epub 2020 Nov 12. PMID: 33180142; PMCID: PMC7817566.
    OpenUrlCrossRefPubMed
  5. ↵
    Beck depression inventory-II. San Antonio, TX: The Psychological Corporation. Beck, A.T., Steer, R.A., Brown, G.K., 1996. Manual for the Beck depression inventory-II. San Antonio, Texas: Psychol. Corp.
  6. ↵
    Berardelli A, Rothwell JC, Thompson PD, Hallett M. Pathophysiology of bradykinesia in Parkinson’s disease. Brain. 2001 Nov;124(Pt 11):2131–46. Review. PMID: 11673316
    OpenUrlCrossRefPubMedWeb of Science
  7. ↵
    Bologna M, Guerra A, Paparella G, Giordo L, Alunni Fegatelli D, Vestri AR, et al. Neurophysiological correlates of bradykinesia in Parkinson’s disease. Brain. 2018 Jun 12. doi: 10.1093/brain/awy155. [Epub ahead of print] PMID: 29901693.
    OpenUrlCrossRefPubMed
  8. ↵
    Braak H, Del Tredici K, Rüb U, de Vos RA, Jansen Steur EN, Braak E. Staging of brain pathology related to sporadic Parkinson’s disease. Neurobiol Aging. 2003 Mar-Apr;24(2):197–211. doi: 10.1016/s0197-4580(02)00065-9. PMID: 12498954.
    OpenUrlCrossRefPubMedWeb of Science
  9. ↵
    Brown RG, Dittner A, Findley L, Wessely SC. The Parkinson fatigue scale. Parkinsonism Relat Disord. 2005 Jan;11(1):49–55. doi: 10.1016/j.parkreldis.2004.07.007. PMID: 15619463.
    OpenUrlCrossRefPubMedWeb of Science
  10. ↵
    Buckner RL, Krienen FM, Castellanos A, Diaz JC, Yeo BT. The organization of the human cerebellum estimated by intrinsic functional connectivity. J Neurophysiol. 2011 Nov;106(5):2322–45. doi: 10.1152/jn.00339.2011. Epub 2011 Jul 27. PMID: 21795627; PMCID: PMC3214121.
    OpenUrlCrossRefPubMedWeb of Science
  11. ↵
    Caligiore D, Mustile M, Spalletta G, Baldassarre G. Action observation and motor imagery for rehabilitation in Parkinson’s disease: A systematic review and an integrative hypothesis. Neurosci Biobehav Rev. 2017 Jan;72:210–222. doi: 10.1016/j.neubiorev.2016.11.005. Epub 2016 Nov 16. PMID: 27865800.
    OpenUrlCrossRefPubMed
  12. ↵
    Catalan MJ, Ishii K, Honda M, Samii A, Hallett M. A PET study of sequential finger movements of varying length in patients with Parkinson’s disease. Brain. 1999 Mar;122 (Pt 3):483–95. PMID: 10094257.
    OpenUrlCrossRefPubMedWeb of Science
  13. ↵
    Chai XJ, Castañón AN, Ongür D, Whitfield-Gabrieli S. Anticorrelations in resting state networks without global signal regression. Neuroimage. 2012 Jan 16;59(2):1420–8. doi: 10.1016/j.neuroimage.2011.08.048. Epub 2011 Aug 26. PMID: 21889994; PMCID: PMC3230748.
    OpenUrlCrossRefPubMedWeb of Science
  14. ↵
    Chee R, Murphy A, Danoudis M, Georgiou-Karistianis N, Iansek R. Gait freezing in Parkinson’s disease and the stride length sequence effect interaction. Brain. 2009 Aug;132(Pt 8):2151–60. doi: 10.1093/brain/awp053. Epub 2009 May 11. PMID: 19433440.
    OpenUrlCrossRefPubMedWeb of Science
  15. ↵
    Chen G, Gully SM, Eden D. Validation of a new general self-efficacy scale. Organizational Research Methods. 2001; 4(62). doi.org/10.1177/109442810141004.
    OpenUrlCrossRefWeb of Science
  16. ↵
    Collins D, Carson HJ. The future for PETTLEP: a modern perspective on an effective and established tool. Curr Opin Psychol. 2017 Aug;16:12–16. doi:10.1016/j.copsyc.2017.03.007. Epub 2017 Mar 23. Review. PMID: 28813334.
    OpenUrlCrossRefPubMed
  17. ↵
    Craig AD. How do you feel--now? The anterior insula and human awareness. Nat Rev Neurosci. 2009 Jan;10(1):59–70. doi: 10.1038/nrn2555. PMID: 19096369.
    OpenUrlCrossRefPubMedWeb of Science
  18. ↵
    Demirci M, Grill S, McShane L, Hallett M. A mismatch between kinesthetic and visual perception in Parkinson’s disease. Ann Neurol. 1997 Jun;41(6):781–8. PMID: 9189039.
    OpenUrlCrossRefPubMedWeb of Science
  19. ↵
    Di Rienzo F, Collet C, Hoyek N, Guillot A. Impact of neurologic deficits on motor imagery: a systematic review of clinical evaluations. Neuropsychol Rev. 2014 Jun;24(2):116–47. doi: 10.1007/s11065-014-9257-6. Epub 2014 Apr 4. Review. PMID: 24700315.
    OpenUrlCrossRefPubMed
  20. ↵
    Epstein RA, Patai EZ, Julian JB, Spiers HJ. The cognitive map in humans: spatial navigation and beyond. Nat Neurosci. 2017 Oct 26;20(11):1504–1513. doi: 10.1038/nn.4656. PMID: 29073650; PMCID: PMC6028313.
    OpenUrlCrossRefPubMed
  21. ↵
    Espay AJ, Giuffrida JP, Chen R, Payne M, Mazzella F, Dunn E, et al. Differential response of speed, amplitude, and rhythm to dopaminergic medications in Parkinson’s disease. Mov Disord. 2011 Dec;26(14):2504–8. doi: 10.1002/mds.23893. Epub 2011 Sep 23. PMID: 21953789; PMCID: PMC3318914.
    OpenUrlCrossRefPubMed
  22. ↵
    Genovese CR, Lazar NA, Nichols T. Thresholding of statistical maps in functional neuroimaging using the false discovery rate. Neuroimage. 2002 Apr;15(4):870–8. doi: 10.1006/nimg.2001.1037. PMID: 11906227.
    OpenUrlCrossRefPubMedWeb of Science
  23. ↵
    Goetz CG, Tilley BC, Shaftman SR, Stebbins GT, Fahn S, Martinez-Martin P, et al. Movement Disorder Society UPDRS Revision Task Force. Movement Disorder Society-sponsored revision of the Unified Parkinson’s Disease Rating Scale (MDS-UPDRS): scale presentation and clinimetric testing results. Mov Disord. 2008 Nov 15;23(15):2129–70. doi:10.1002/mds.22340. PMID: 19025984.
    OpenUrlCrossRefPubMedWeb of Science
  24. ↵
    1. R. Scwarzer
    Greenglass E & Schwarzer R (1998). The Proactive Coping Inventory (PCI). In R. Scwarzer (Ed), Advances in Health Psychology Research. Berlin: Free University of Berlin. Institut for Arbeits, Organizations-und Gesundheitpsychologie (ISBN 3-00-002776-9).
  25. ↵
    Guillot A, Collet C. Construction of the motor imagery integrative model in sport: a review and theoretical investigation of motor imagery use. International Review of Sport and Exercise Psychology. 2008 Mar 1;1:31–44.
    OpenUrl
  26. ↵
    Guillot A, Collet C, Nguyen VA, Malouin F, Richards C, Doyon J. Brain activity during visual versus kinesthetic imagery: an fMRI study. Hum Brain Mapp. 2009 Jul;30(7):2157–72. doi: 10.1002/hbm.20658. PMID: 18819106.
    OpenUrlCrossRefPubMedWeb of Science
  27. ↵
    Guillot A, Di Rienzo F, Collet C. The neurofunctional architecture of motor imagery. 2014 in Advanced Brain Neuroimaging Topics in Health and Disease - Methods and Applications, Intech pp. 433–456
  28. ↵
    Hétu S, Grégoire M, Saimpont A, Coll MP, Eugène F, Michon PE, Jackson PL. The neural network of motor imagery: an ALE meta-analysis. Neurosci Biobehav Rev. 2013 Jun;37(5):930–49. doi: 10.1016/j.neubiorev.2013.03.017. Epub 2013 Apr 10. PMID: 23583615.
    OpenUrlCrossRefPubMed
  29. Hoehn MM, Yahr MD. Parkinsonism: onset, progression, and mortality. 1967. Neurology. 1998 Feb;50(2):318–43. PMID: 9484345.
    OpenUrlCrossRefPubMed
  30. ↵
    Hughes AJ, Daniel SE, Kilford L, Lees AJ. Accuracy of clinical diagnosis of idiopathic Parkinson’s disease: a clinico-pathological study of 100 cases. J Neurol Neurosurg Psychiatry. 1992 Mar;55(3):181–4. PMID: 1564476; PMCID: PMC1014720.
    OpenUrlAbstract/FREE Full Text
  31. ↵
    Jahanshahi M, Jenkins IH, Brown RG, Marsden CD, Passingham RE, Brooks DJ. Self-initiated versus externally triggered movements. I. An investigation using measurement of regional cerebral blood flow with PET and movement-related potentials in normal and Parkinson’s disease subjects. Brain. 1995 Aug;118 (Pt4):913–33. PMID: 7655888.
    OpenUrlCrossRefPubMedWeb of Science
  32. ↵
    Jenkinson C, Fitzpatrick R, Peto V, Greenhall R, Hyman N. The Parkinson’s Disease Questionnaire (PDQ-39): development and validation of a Parkinson’s disease summary index score. Age Ageing. 1997 Sep;26(5):353–7. doi: 10.1093/ageing/26.5.353. PMID: 9351479.
    OpenUrlCrossRefPubMedWeb of Science
  33. ↵
    Johnson KA, Hartwell K, LeMatty T, Borckardt J, Morgan PS, Govindarajan K, et al. Intermittent “real-time” fMRI feedback is superior to continuous presentation for a motor imagery task: a pilot study. J Neuroimaging. 2012 Jan;22(1):58–66. doi: 10.1111/j.1552-6569.2010.00529.x. Epub 2010 Oct 26. PMID: 20977537; PMCID: PMC3070198.
    OpenUrlCrossRefPubMed
  34. ↵
    Kang SY, Wasaka T, Shamim EA, Auh S, Ueki Y, Dang N, et al. The sequence effect in de novo Parkinson’s disease. J Mov Disord. 2011 May;4(1):38–40. doi: 10.14802/jmd.11006. Epub 2011 Apr 30. PMID: 24868390; PMCID: PMC4027704.
    OpenUrlCrossRefPubMed
  35. ↵
    Kober SE, Witte M, Ninaus M, Neuper C, Wood G. Learning to modulate one’s own brain activity: the effect of spontaneous mental strategies. Front Hum Neurosci. 2013 Oct 18;7:695. doi: 10.3389/fnhum.2013.00695. PMID: 24151462; PMCID: PMC3798979.
    OpenUrlCrossRefPubMed
  36. ↵
    Kringelbach ML. The human orbitofrontal cortex: linking reward to hedonic experience. Nat Rev Neurosci. 2005 Sep;6(9):691–702. doi: 10.1038/nrn1747. PMID: 16136173.
    OpenUrlCrossRefPubMedWeb of Science
  37. ↵
    Leech R, Sharp DJ. The role of the posterior cingulate cortex in cognition and disease. Brain. 2014 Jan;137(Pt 1):12–32. doi: 10.1093/brain/awt162. Epub 2013 Jul 18. PMID: 23869106; PMCID: PMC3891440.
    OpenUrlCrossRefPubMedWeb of Science
  38. ↵
    1. Benedetti F.,
    2. Enck P.,
    3. Frisaldi E.,
    4. Schedlowski M
    Lidstone SC. (2014) Great Expectations: The Placebo Effect in Parkinson’s Disease. In: Benedetti F., Enck P., Frisaldi E., Schedlowski M. (eds) Placebo. Handbook of Experimental Pharmacology, vol 225. Springer, Berlin, Heidelberg. https://doi.org/10.1007/978-3-662-44519-8_8
  39. ↵
    Ling H, Massey LA, Lees AJ, Brown P, Day BL. Hypokinesia without decrement distinguishes progressive supranuclear palsy from Parkinson’s disease. Brain. 2012 Apr;135(Pt 4):1141–53. doi: 10.1093/brain/aws038. Epub 2012 Mar 6. PMID: 22396397; PMCID: PMC3326257.
    OpenUrlCrossRefPubMedWeb of Science
  40. ↵
    Lorey B, Bischoff M, Pilgramm S, Stark R, Munzert J, Zentgraf K. The embodied nature of motor imagery: the influence of posture and perspective. Exp Brain Res. 2009 Apr;194(2):233–43. doi: 10.1007/s00221-008-1693-1. Epub 2009 Jan 13. PMID: 19139856.
    OpenUrlCrossRefPubMedWeb of Science
  41. ↵
    Marchesotti S, Bassolino M, Serino A, Bleuler H, Blanke O. Quantifying the role of motor imagery in brain-machine interfaces. Sci Rep. 2016 Apr 7;6:24076. doi: 10.1038/srep24076. PMID: 27052520; PMCID: PMC4823701.
    OpenUrlCrossRefPubMed
  42. ↵
    Marinus J, Visser M, Verwey NA, Verhey FR, Middelkoop HA, Stiggelbout AM, van Hilten JJ. Assessment of cognition in Parkinson’s disease. Neurology. 2003 Nov 11;61(9):1222–8. PMID: 14610124.
    OpenUrlAbstract/FREE Full Text
  43. ↵
    Mazzoni P, Hristova A, Krakauer JW. Why don’t we move faster? Parkinson’s disease, movement vigor, and implicit motivation. J Neurosci. 2007 Jul 4;27(27):7105–16. PMID: 17611263.
    OpenUrlAbstract/FREE Full Text
  44. ↵
    Megumi F, Yamashita A, Kawato M, Imamizu H. Functional MRI neurofeedback training on connectivity between two regions induces long-lasting changes in intrinsic functional network. Front Hum Neurosci. 2015 Mar 30;9:160. doi: 10.3389/fnhum.2015.00160. PMID: 25870552; PMCID: PMC4377493.
    OpenUrlCrossRefPubMed
  45. ↵
    Mehler DMA, Williams AN, Krause F, Lührs M, Wise RG, Turner DL, Linden DEJ, Whittaker JR. The BOLD response in primary motor cortex and supplementary motor area during kinesthetic motor imagery based graded fMRI neurofeedback. Neuroimage. 2019 Jan 1;184:36–44. doi: 10.1016/j.neuroimage.2018.09.007. Epub 2018 Sep 8. PMID: 30205210; PMCID: PMC6264383.
    OpenUrlCrossRefPubMed
  46. ↵
    Morris ME, Iansek R, Galna B. Gait festination and freezing in Parkinson’s disease: pathogenesis and rehabilitation. Mov Disord. 2008;23 Suppl 2:S451–60. doi: 10.1002/mds.21974. Review. Erratum in: Mov Disord. 2008 Aug 15;23(11):1639-40. PMID: 18668618.
    OpenUrlCrossRefPubMedWeb of Science
  47. ↵
    Nakamura T, Ghilardi MF, Mentis M, Dhawan V, Fukuda M, Hacking A, et al. Functional networks in motor sequence learning: abnormal topographies in Parkinson’s disease. Hum Brain Mapp. 2001 Jan;12(1):42–60. PMID: 11198104.
    OpenUrlCrossRefPubMedWeb of Science
  48. ↵
    Nasreddine ZS, Phillips NA, Bédirian V, Charbonneau S, Whitehead V, Collin I, et al. The Montreal Cognitive Assessment, MoCA: a brief screening tool for mild cognitive impairment. J Am Geriatr Soc. 2005 Apr;53(4):695–9. PMID: 15817019.
    OpenUrlCrossRefPubMedWeb of Science
  49. ↵
    Olsson CJ, Jonsson B, Nyberg L. Learning by doing and learning by thinking: an FMRI study of combining motor and mental training. Front Hum Neurosci. 2008 Jun 18;2:5. doi: 10.3389/neuro.09.005.2008. PMID: 18958206; PMCID: PMC2525972.
    OpenUrlCrossRefPubMed
  50. ↵
    Pearson J. The human imagination: the cognitive neuroscience of visual mental imagery. Nat Rev Neurosci. 2019 Oct;20(10):624–634. doi: 10.1038/s41583-019-0202-9. PMID: 31384033.
    OpenUrlCrossRefPubMed
  51. ↵
    Rance M, Walsh C, Sukhodolsky DG, Pittman B, Qiu M, Kichuk SA, Wasylink S, Koller WN, Bloch M, Gruner P, Scheinost D, Pittenger C, Hampson M. Time course of clinical change following neurofeedback. Neuroimage. 2018 Nov 1;181:807–813. doi: 10.1016/j.neuroimage.2018.05.001. Epub 2018 May 2. PMID: 29729393; PMCID: PMC6454268.
    OpenUrlCrossRefPubMed
  52. ↵
    Reuben DB, Siu AL. An objective measure of physical function of elderly outpatients. The Physical Performance Test. J Am Geriatr Soc. 1990 Oct;38(10):1105–12. PMID: 2229864.
    OpenUrlCrossRefPubMedWeb of Science
  53. ↵
    Ros T, Enriquez-Geppert S, Zotev V, Young KD, Wood G, Whitfield-Gabrieli S, et al. Consensus on the reporting and experimental design of clinical and cognitive-behavioural neurofeedback studies (CRED-nf checklist). Brain. 2020 Jun 1;143(6):1674–1685. doi: 10.1093/brain/awaa009. PMID: 32176800; PMCID: PMC7296848.
    OpenUrlCrossRefPubMed
  54. ↵
    Sabatini U, Boulanouar K, Fabre N, Martin F, Carel C, Colonnese C, et al. Cortical motor reorganization in akinetic patients with Parkinson’s disease: a functional MRI study. Brain. 2000 Feb;123 (Pt 2):394–403. PMID: 10648446.
    OpenUrlCrossRefPubMedWeb of Science
  55. ↵
    Samuel M, Ceballos-Baumann AO, Blin J, Uema T, Boecker H, Passingham RE, et al. Evidence for lateral premotor and parietal overactivity in Parkinson’s disease during sequential and bimanual movements. A PET study. Brain. 1997 Jun;120 (Pt 6):963–76. PMID: 9217681.
    OpenUrlCrossRefPubMedWeb of Science
  56. ↵
    Scheinost D, Hampson M, Qiu M, Bhawnani J, Constable RT, Papademetris X. A graphics processing unit accelerated motion correction algorithm and modular system for real-time fMRI. Neuroinformatics. 2013 Jul;11(3):291–300. doi:10.1007/s12021-013-9176-3. PMID: 23319241; PMCID: PMC3644525.
    OpenUrlCrossRefPubMed
  57. ↵
    Scheinost D, Hsu TW, Avery EW, Hampson M, Constable RT, Chun MM, Rosenberg MD. Connectome-based neurofeedback: A pilot study to improve sustained attention. Neuroimage. 2020 May 15;212:116684. doi: 10.1016/j.neuroimage.2020.116684. Epub 2020 Feb 27. PMID: 32114151; PMCID: PMC7165055.
    OpenUrlCrossRefPubMed
  58. ↵
    Schrag A, Sampaio C, Counsell N, Poewe W. Minimal clinically important change on the unified Parkinson’s disease rating scale. Mov Disord. 2006 Aug;21(8):1200–7. PMID: 16673410.
    OpenUrlCrossRefPubMedWeb of Science
  59. ↵
    Schwarzer R, Diehl M, Schmitz GS (1999). Self-Regulation Scale. Online source: http://www.fu-berlin.de/gesund/skalen.
  60. ↵
    Sepulveda P, Sitaram R, Rana M, Montalba C, Tejos C, Ruiz S. How feedback, motor imagery, and reward influence brain self-regulation using real-time fMRI. Hum Brain Mapp. 2016 Sep;37(9):3153–71. doi: 10.1002/hbm.23228. Epub 2016 Jun 6. PMID: 27272616; PMCID: PMC6867497.
    OpenUrlCrossRefPubMed
  61. ↵
    Shen X, Tokoglu F, Papademetris X, Constable RT. Groupwise whole-brain parcellation from resting-state fMRI data for network node identification. Neuroimage. 2013 Nov 15;82:403–15. doi: 10.1016/j.neuroimage.2013.05.081. Epub 2013 Jun 4. PMID: 23747961; PMCID: PMC3759540.
    OpenUrlCrossRefPubMedWeb of Science
  62. ↵
    Shulman LM, Gruber-Baldini AL, Anderson KE, Fishman PS, Reich SG, Weiner WJ. The clinically important difference on the unified Parkinson’s disease rating scale. Arch Neurol. 2010 Jan;67(1):64–70. doi: 10.1001/archneurol.2009.295. PMID: 20065131.
    OpenUrlCrossRefPubMedWeb of Science
  63. ↵
    Sorger B, Scharnowski F, Linden DEJ, Hampson M, Young KD. Control freaks: Towards optimal selection of control conditions for fMRI neurofeedback studies. Neuroimage. 2019 Feb 1;186:256–265. doi: 10.1016/j.neuroimage.2018.11.004. Epub 2018 Nov 10. PMID: 30423429; PMCID: PMC6338498.
    OpenUrlCrossRefPubMed
  64. ↵
    Spielberger CD, Gorsuch RL, Lushene R, Vagg PR, Jacobs GA, 1983. Manual for the State-Trait Anxiety Inventory. Palo Alto, California: Consulting Psychologists Press.
  65. ↵
    Starkstein SE, Mayberg HS, Preziosi TJ, Andrezejewski P, Leiguarda R, Robinson RG. Reliability, validity, and clinical correlates of apathy in Parkinson’s disease. J Neuropsychiatry Clin Neurosci. 1992 Spring;4(2):134–9. doi: 10.1176/jnp.4.2.134. PMID: 1627973.
    OpenUrlCrossRefPubMedWeb of Science
  66. ↵
    Stoeckel LE, Garrison KA, Ghosh S, Wighton P, Hanlon CA, Gilman JM, et al. Optimizing real time fMRI neurofeedback for therapeutic discovery and development. Neuroimage Clin. 2014 Jul 10;5:245–55. doi:10.1016/j.nicl.2014.07.002. eCollection 2014. Review. PMID: 25161891; PMCID: PMC4141981.
    OpenUrlCrossRefPubMed
  67. ↵
    Subramanian L, Morris MB, Brosnan M, Turner DL, Morris HR, Linden DE. Functional Magnetic Resonance Imaging Neurofeedback-guided Motor Imagery Training and Motor Training for Parkinson’s Disease: Randomized Trial. Front Behav Neurosci. 2016 Jun 8;10:111. doi: 10.3389/fnbeh.2016.00111. eCollection 2016. PMID: 27375451; PMCID: PMC4896907.
    OpenUrlCrossRefPubMed
  68. ↵
    Tamir R, Dickstein R, Huberman M. Integration of motor imagery and physical practice in group treatment applied to subjects with Parkinson’s disease. Neurorehabil Neural Repair. 2007 Jan-Feb;21(1):68–75. PMID: 17172556.
    OpenUrlCrossRefPubMedWeb of Science
  69. ↵
    Thibault RT, MacPherson A, Lifshitz M, Roth RR, Raz A. Neurofeedback with fMRI: A critical systematic review. Neuroimage. 2018 May 15;172:786–807. doi: 10.1016/j.neuroimage.2017.12.071. Epub 2017 Dec 27. PMID: 29288868.
    OpenUrlCrossRefPubMed
  70. ↵
    Tinaz S, Para K, Vives-Rodriguez A, Martinez-Kaigi V, Nalamada K, Sezgin M, et al. Insula as the Interface Between Body Awareness and Movement: A Neurofeedback-Guided Kinesthetic Motor Imagery Study in Parkinson’s Disease. Front Hum Neurosci. 2018 Dec 7;12:496. doi: 10.3389/fnhum.2018.00496. PMID: 30581383; PMCID: PMC6292989.
    OpenUrlCrossRefPubMed
  71. ↵
    Tinaz S, Pillai AS, Hallett M. Sequence Effect in Parkinson’s Disease Is Related to Motor Energetic Cost. Front Neurol. 2016 May 24;7:83. doi: 10.3389/fneur.2016.00083. eCollection 2016. PMID: 27252678; PMCID: PMC4877367.
    OpenUrlCrossRefPubMed
  72. ↵
    Tosserams A, Nijkrake MJ, Sturkenboom IHWM, Bloem BR, Nonnekes J. Perceptions of Compensation Strategies for Gait Impairments in Parkinson’s Disease: A Survey Among 320 Healthcare Professionals. J Parkinsons Dis. 2020;10(4):1775–1778. doi: 10.3233/JPD-202176. PMID: 32925103; PMCID: PMC7683052.
    OpenUrlCrossRefPubMed
  73. ↵
    Whitfield-Gabrieli S, Nieto-Castanon A. Conn: a functional connectivity toolbox for correlated and anticorrelated brain networks. Brain Connect. 2012;2(3):125–41. doi: 10.1089/brain.2012.0073. Epub 2012 Jul 19. PMID: 22642651.
    OpenUrlCrossRefPubMed
  74. ↵
    Williams SE, Cumming J, Ntoumanis N, Nordin-Bates SM, Ramsey R. Further validation and development of the movement imagery questionnaire. J Sport Exerc Psychol. 2012 Oct;34(5):621–46. PMID: 23027231.
    OpenUrlPubMed
  75. ↵
    Yamashita A, Hayasaka S, Kawato M, Imamizu H. Connectivity Neurofeedback Training Can Differentially Change Functional Connectivity and Cognitive Performance. Cereb Cortex. 2017 Oct 1;27(10):4960–4970. doi: 10.1093/cercor/bhx177. PMID: 28922830.
    OpenUrlCrossRefPubMed
  76. ↵
    Yu H, Sternad D, Corcos DM, Vaillancourt DE. Role of hyperactive cerebellum and motor cortex in Parkinson’s disease. Neuroimage. 2007 Mar;35(1):222–33. Epub 2007 Jan 12. PMID: 17223579; PMCID: PMC1853309.
    OpenUrlCrossRefPubMedWeb of Science
  77. ↵
    Zhao W, Hill MD, Palesch Y. Minimal sufficient balance-a new strategy to balance baseline covariates and preserve randomness of treatment allocation. Stat Methods Med Res. 2015 Dec;24(6):989–1002. doi: 10.1177/0962280212436447. Epub 2012 Jan 26. PMID: 22287602; PMCID: PMC3474894.
    OpenUrlCrossRefPubMed
Back to top
PreviousNext
Posted January 06, 2022.
Download PDF
Data/Code
Email

Thank you for your interest in spreading the word about medRxiv.

NOTE: Your email address is requested solely to identify you as the sender of this article.

Enter multiple addresses on separate lines or separate them with commas.
Neurofeedback-guided kinesthetic motor imagery training in Parkinson’s disease: Randomized trial
(Your Name) has forwarded a page to you from medRxiv
(Your Name) thought you would like to see this page from the medRxiv website.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Share
Neurofeedback-guided kinesthetic motor imagery training in Parkinson’s disease: Randomized trial
Sule Tinaz, Serageldin Kamel, Sai S. Aravala, Mohamed Elfil, Ahmed Bayoumi, Amar Patel, Dustin Scheinost, Rajita Sinha, Michelle Hampson
medRxiv 2022.01.05.22268816; doi: https://doi.org/10.1101/2022.01.05.22268816
Reddit logo Twitter logo Facebook logo LinkedIn logo Mendeley logo
Citation Tools
Neurofeedback-guided kinesthetic motor imagery training in Parkinson’s disease: Randomized trial
Sule Tinaz, Serageldin Kamel, Sai S. Aravala, Mohamed Elfil, Ahmed Bayoumi, Amar Patel, Dustin Scheinost, Rajita Sinha, Michelle Hampson
medRxiv 2022.01.05.22268816; doi: https://doi.org/10.1101/2022.01.05.22268816

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Subject Area

  • Neurology
Subject Areas
All Articles
  • Addiction Medicine (228)
  • Allergy and Immunology (504)
  • Anesthesia (110)
  • Cardiovascular Medicine (1238)
  • Dentistry and Oral Medicine (206)
  • Dermatology (147)
  • Emergency Medicine (282)
  • Endocrinology (including Diabetes Mellitus and Metabolic Disease) (531)
  • Epidemiology (10021)
  • Forensic Medicine (5)
  • Gastroenterology (499)
  • Genetic and Genomic Medicine (2453)
  • Geriatric Medicine (238)
  • Health Economics (479)
  • Health Informatics (1643)
  • Health Policy (752)
  • Health Systems and Quality Improvement (636)
  • Hematology (248)
  • HIV/AIDS (533)
  • Infectious Diseases (except HIV/AIDS) (11864)
  • Intensive Care and Critical Care Medicine (626)
  • Medical Education (252)
  • Medical Ethics (75)
  • Nephrology (268)
  • Neurology (2280)
  • Nursing (139)
  • Nutrition (352)
  • Obstetrics and Gynecology (454)
  • Occupational and Environmental Health (536)
  • Oncology (1245)
  • Ophthalmology (377)
  • Orthopedics (134)
  • Otolaryngology (226)
  • Pain Medicine (157)
  • Palliative Medicine (50)
  • Pathology (324)
  • Pediatrics (730)
  • Pharmacology and Therapeutics (313)
  • Primary Care Research (282)
  • Psychiatry and Clinical Psychology (2280)
  • Public and Global Health (4834)
  • Radiology and Imaging (837)
  • Rehabilitation Medicine and Physical Therapy (491)
  • Respiratory Medicine (651)
  • Rheumatology (285)
  • Sexual and Reproductive Health (238)
  • Sports Medicine (227)
  • Surgery (267)
  • Toxicology (44)
  • Transplantation (125)
  • Urology (99)