Brain connectivity during simulated balance in older adults with and without Parkinson's disease.

Individuals with Parkinson's disease often experience postural instability, a debilitating and largely treatment-resistant symptom. A better understanding of the neural substrates contributing to postural instability could lead to more effective treatments. Constraints of current functional neuroimaging techniques, such as the horizontal orientation of most MRI scanners (forcing participants to lie supine), complicates investigating cortical and subcortical activation patterns and connectivity networks involved in healthy and parkinsonian balance control. In this cross-sectional study, we utilized a newly-validated MRI-compatible balance simulator (based on an inverted pendulum) that enabled participants to perform balance-relevant tasks while supine in the scanner. We utilized functional MRI to explore effective connectivity underlying static and dynamic balance control in healthy older adults (n = 17) and individuals with Parkinson's disease while on medication (n = 17). Participants performed four tasks within the scanner with eyes closed: resting, proprioceptive tracking of passive ankle movement, static balancing of the simulator, and dynamic responses to random perturbations of the simulator. All analyses were done in the participant's native space without spatial transformation to a common template. Effective connectivity between 57 regions of interest was computed using a Bayesian Network learning approach with false discovery rate set to 5%. The first 12 principal components of the connection weights, binomial logistic regression, and cross-validation were used to create 4 separate models: contrasting static balancing vs {rest, proprioception} and dynamic balancing vs {rest, proprioception} for both controls and individuals with Parkinson's disease. In order to directly compare relevant connections between controls and individuals with Parkinson's disease, we used connections relevant for predicting a task in either controls or individuals with Parkinson's disease in logistic regression with Least Absolute Shrinkage and Selection Operator regularization. During dynamic balancing, we observed decreased connectivity between different motor areas and increased connectivity from the brainstem to several cortical and subcortical areas in controls, while individuals with Parkinson's disease showed increased connectivity associated with motor and parietal areas, and decreased connectivity from brainstem to other subcortical areas. No significant models were found for static balancing in either group. Our results support the notion that dynamic balance control in individuals with Parkinson's disease relies more on cortical motor areas compared to healthy older adults, who show a preference of subcortical control during dynamic balancing.

Introduction

Postural instability is a debilitating symptom of Parkinson’s disease that is often treatment-resistant and a precursor to falls (Boonstraet al., 2008, Grimbergenet al., 2009, Crouseet al., 2016, Fasanoet al., 2017). Effective treatment of postural instability in Parkinson’s disease is significantly impeded by a relatively poor understanding of both the neural networks involved in healthy balance and the pathophysiology underlying balance deficits in individuals with Parkinson’s disease (Grimbergenet al., 2009, Kimet al., 2013, Crouseet al., 2016). Balance control in healthy individuals likely involves an integrated network of cortical and subcortical structures (Takakusaki, 2017), however, the specific networks involved remain unresolved. Moreover, the extent to which changes in the activation of, and connection between, cortical and subcortical structures contribute to balance deficits in Parkinson’s disease has not been established.

Investigating neural activation patterns and connectivity networks contributing to healthy and abnormal balance control is hindered by the constraints of current functional neuroimaging options. Although portable neuroimaging techniques, such as EEG and functional near-infrared spectroscopy (fNIRS), allow for functional neuroimaging in participants standing upright (albeit frequently with unacceptable levels of movement artifact), their ability to record from subcortical structures is limited. Both cortical and subcortical activity can be measured using PET and functional MRI (fMRI) scanners, however, these scanners are almost exclusively horizontal, requiring the participant to be lying supine. Wearable PET scanners have recently been developed (Baueret al., 2016, Melroyet al., 2017), but are limited by poor temporal resolution and the weight of the wearable PET device, potentially interfering with balance control.

Motor imagery of static and dynamic balance tasks offers a less-than-ideal approach to probing balance networks but has been successfully employed in healthy participants (Malouinet al., 2003, Jahnet al., 2004, Jahnet al., 2008, Zwergalet al., 2012, Ferrayeet al., 2014, Taubeet al., 2015, Bhattet al., 2018, Mouthonet al., 2018, Gilatet al., 2019) and individuals with Parkinson’s disease (Petersonet al., 2014, Gilatet al., 2019). Unsurprisingly, some brain areas are more strongly, or selectively, activated during actual motor execution of a task compared to motor imagery of the same task, and vice versa (Guillotet al., 2012, O'Shea and Moran, 2017). The ability to perform motor imagery varies greatly among individuals (Saimpont et al., 2015) and declines with age for complex movements (Saimpontet al., 2013, Kalicinskiet al., 2015), making it less suitable for investigating balance networks in healthy older participants and individuals with Parkinson’s disease.

A few studies have attempted performance of balance-related tasks while recording fMRI data in participants lying supine (Karimet al., 2014, de Lima-Pardiniet al., 2017). Karim et al. (2014) constructed an MRI-compatible force platform and had participants use visual feedback to generate anterior-posterior ankle torque using feed-forward volitional control. Likewise, De Lima-Pardini et al. (2017) developed an MRI compatible force measurement system to investigate anticipatory postural adjustments during single leg raises to simulate step initiation. Although these prior studies had participants perform balance-related tasks while lying down, neither truly simulated free-standing balance.

Recently, we developed and validated a novel MRI-compatible balance simulator able to detect postural instability in individuals with Parkinson’s disease (Pasman et al., 2019). The balance simulator was shown to be relatively easy to control using muscles around the ankle joint for both healthy participants and individuals with Parkinson’s disease after only a few minutes of practice outside the scanner. During validation, lab studies verified that qualitatively similar balance behaviour was seen in participants when maintaining balance of their body during upright standing and when maintaining balance using the balance simulator while supine with the eyes open and eyes closed. The balance simulator was also sensitive enough to detect deficits in both static and dynamic reactive balance control in individuals with Parkinson’s disease compared to healthy age-matched participants, independent of vision (Pasman et al., 2019).

Here we use this newly-validated balance simulator during fMRI scanning to investigate effective connectivity during static and dynamic balance control tasks in individuals with Parkinson’s disease and healthy older adults. As abnormal brain connectivity patterns have previously been found in individuals with Parkinson’s disease at rest and during a variety of motor and non-motor tasks (Filippi et al., 2018), we hypothesized differences in brain connectivity networks between individuals with Parkinson’s disease and healthy elderly during performance of simulated free-standing static and dynamic balance tasks.

Materials and methods

Participants

Twenty individuals with Parkinson’s disease and 19 age-matched elderly controls, recruited from the community and through the Pacific Parkinson’s Research Centre in Vancouver, participated in this study. Exclusion criteria for individuals with Parkinson’s disease were any of the following medical issues (self-reported during initial telephone screening): presence of atypical parkinsonism; any prior neurosurgical procedures such as deep brain stimulation; excessive levodopa-induced dyskinesias that impaired their balance; botulinum toxin injections in lower leg muscles within the last 3 months; documented proprioceptive loss (e.g., abnormal vibratory sense, altered joint position sense, etc.); dementia precluding informed consent; history of other neurological disease (e.g., stroke, seizures); and medical issues (other than Parkinson’s disease) that influenced their balance. Controls were excluded if during initial phone screening they self-reported any medical conditions that influenced their balance. Participants in both groups were also excluded if they exceeded the height (max height: 182 cm) and weight (max weight participant and balance simulator combined: 136 kg) restrictions of the MRI scanner, and if they had any contraindications precluding them from undergoing MRI scanning. All participants were fluent in English. All participants provided written informed consent prior to testing and followed experimental procedures that were approved by the UBC Clinical Research Ethics Board, the Vancouver Coastal Health Research Institute, and the UBC MRI Research Centre. Data collection took place at the UBC MRI Research Centre from April 2016 until October 2016. Five participants were excluded due to: inability to perform the simulated balancing tasks during the familiarization session in the laboratory (n = 2), inability to complete the protocol in the MRI scanner (n = 2), and finding of an incidental abnormality on their anatomical MRI scan (n = 1). Therefore, 17 individuals with Parkinson’s disease and 17 controls were included in the final data set (Table 1).

Table 1

Baseline participant characteristics.

	Parkinson	Control	Statistics
General information
Sample size	17	17
Age (years)	67.6 (1.1)	68.1 (1.3)	p = 0.810
Number of women (%)	8 (47%)	10 (59%)	p = 0.732
Height (cm)	168.7 (1.5)	166.2 (2.2)	p = 0.350
Weight (kg)	71.2 (3.2)	65.8 (3.7)	p = 0.282
Fallers (≤6 months)	4 (24%)	1 (6%)	p = 0.335

Data are displayed as mean (SE) or number of persons (percentage between parentheses).

Individuals with Parkinson’s disease were examined approximately one hour after intake of their regular antiparkinson medication to coincide with their subjectively best clinical ‘on’ condition and specifically assess dopa-unresponsive balance effects. All participants completed a brief balance-oriented medical history survey, including the occurrence of prior falls within the past 6 months. The Hoehn and Yahr (H&Y) scale (Goetz et al., 2004) and Unified Parkinson’s Disease Rating Scale motor examination (UPDRS-ME) (Goetz et al., 2008) were administered to individuals with Parkinson’s disease (Table 2). The clinical assessments were done on a separate day prior to MRI scanning (range 0–31 days, mean ± SD 9 ± 8 days) when participants were familiarized, in the laboratory, with the balance simulator and the tasks they would be performing in the MRI scanner.

Table 2

Clinical characteristics individuals with Parkinson’s disease.

PDON	Age (years)	Disease duration (years)	UPDRS-ME score	Hoehn & Yahr score	Levodopa equivalent dose	Antiparkinson medication
1	68	2	8	1	225	Levodopa/carbidopa
2	72	2	31	2	300	Levodopa/carbidopa
3	74	4	25	2	325	Levodopa/carbidopa, Rasagiline
4	67	3	44	2	850	Levodopa/carbidopa
5	73	6	36	2	900	Levodopa/carbidopa, Rasagiline, Pramipexole
6	60	2	7	1	300	Levodopa/carbidopa
7	61	9	18	2	750	Levodopa/carbidopa
8	64	14	46	3	1150	Levodopa/carbidopa, Pramipexole
9	75	4	31	2	775	Levodopa/carbidopa, Rasagiline
10	70	7	31	3	1450	Levodopa/carbidopa, Rasagiline, Amantadine
11	68	10	20	2	900	Levodopa/carbidopa
12	59	7	19	1	450	Levodopa/carbidopa
13	68	6	20	2	948	Levodopa/carbidopa, Entacapone
14	69	1	34	2	150	Levodopa/carbidopa
15	65	12	16	2	1015	Levodopa/carbidopa, Rasagiline, Ropinirol
16	69	4	28	2	200	Levodopa/carbidopa
17	67	6	39	2	600	Levodopa/carbidopa, Rasagiline
Range	59–75	1–14	7–46	1–3	150–1450
Mean (SE)	67.6 (1.1)	5.8 (0.9)	26.6 (2.8)	–	664.0 (92.5)

PDON, Individuals with Parkinson’s disease; Maximum Unified Parkinson’s Disease Rating Scale motor examination (UPDRS-ME) score is 108; Maximum Hoehn & Yahr (H&Y) score is 5.

Apparatus

During the simulated stance trials participants were asked to control a customized balance simulator (“simulator”) (Pasman et al., 2019). Briefly, in the simulator the participant lay supine with their feet placed against a footplate that controlled a free‐standing inverted pendulum in the anterior-posterior direction about an axis aligned with the ankle joints (Fig. 1A and 1B). The mass on the simulator was adjusted such that the estimated load stiffness was about 60% of the ankle load stiffness seen during upright quiet stance (range 38–81%, mean ± SD 60 ± 8.7%). Participants controlled the movement of the inverted pendulum by contracting ankle plantar-flexors and dorsi-flexors, similar to the control of upright stance (Billot et al., 2010). To mimic the gravitational force that pulls the body towards the ground during upright stance, participants had tightly-fastened adjustable straps wrapped around the waist and shoulders which were then attached to the base of the simulator and tightened as much as possible without discomfort, to increase the force between the participant’s feet and the simulator foot-plate. The straps helped ensure the participant’s heels would always remain in contact with the simulator footplate when contracting their plantar-flexor muscles and therefore movement of the simulator could generate relevant sensory input at the foot and ankle (e.g., inputs from joint receptors, Golgi tendon organs, muscle spindles, and foot sole cutaneous receptors) during the simulated stance trials. Mechanical stops were used to limit the simulator to a range of ±17° to ensure participant safety. In addition, the experimenter stood beside the simulator at all times throughout the testing to assist if needed.

Fig. 1

Experimental set-up of the balance simulator in the MRI environment (A and B), and functional MRI scanning experimental protocol (C) with the order of the static and dynamic balancing tasks counterbalanced across participants and the proprioceptive task always presented in between the balancing tasks.

Experimental tasks

Participants performed 4 different tasks: resting, proprioceptive, static balancing (SB), and dynamic balancing (DB). Since Parkinson’s disease effects on simulated SB performance were previously found to be independent of vision (Pasman et al., 2019), all tasks were performed with the eyes closed to prevent the activation of cortical regions associated with processing of visual information. The resting and proprioceptive tasks were used as control tasks. During the resting state portion, participants were instructed to remain awake and lay still. During the proprioceptive task, the simulator was moved by the experimenter and participants were instructed to continuously track the passive movement of the ankle joints by moving their left index finger. Using the proprioceptive task as one of the control tasks ensured brain connectivity seen during the balancing tasks was due to balance related activation, not solely proprioception. In the SB task, participants were instructed to keep the simulator as still as possible. In the DB task, participants were instructed to keep the simulator balanced while responding to transient, random, anterior-posterior perturbations applied to the simulator by the experimenter using a hand-held bar.

During the familiarization session in the laboratory, all participants were given the opportunity to practice the proprioceptive, SB, and DB tasks until they felt comfortable. Subsequently, they performed the tasks for 2-minute trials, corresponding to the experimental protocol used during the MRI scanning session (Fig. 1C). Additionally, the height and weight of the participants were measured to ensure the height and weight restrictions of the MRI scanning table would not be exceeded. Finally, participants who were unable to perform the SB and/or DB balancing tasks during the familiarization session were excluded from MRI scanning participation.

During the MRI scanning session, the resting task was performed first (Fig. 1C). The order of the SB and DB tasks was counterbalance across participants, but the proprioceptive task was always presented in between the balancing tasks to minimize fatigue. In addition, we allowed 2 min of rest between each task.

MRI scanning

MRI data were collected on a 3 Tesla scanner (Philips Achieva 3.0 T R3.2; Philips Medical Systems, The Netherlands) equipped with an 8-channel head-coil. Head motion was minimized by a strap placed around the participant's head within the head coil, and foam wedged between the participant's head and the head coil. Participants wore ear plugs and earmuffs to minimize scanner noise. An MRI safe sandbag was placed on top of the participant's pelvic area to further stabilize the hips during the balancing tasks.

During fMRI scanning, single-shot gradient-echo echo‐planar T2*‐weighted images with blood oxygenation level-dependent contrast were taken of the whole brain, including the cerebellum and brainstem. Scanning parameters were: repetition time 2000 ms, echo time 30 ms, flip angle 90°, field of view 240 × 240 × 143 mm, matrix size 80 × 80, pixel size 3.0 × 3.0 mm. Thirty‐six axial slices of 3.97 mm thickness were collected in each volume, with a gap thickness of 1 mm. A high resolution, 3‐dimensional T1‐weighted image of the whole brain was acquired, using a multi-shot Turbo Field Echo sequence (repetition/echo time 7.7/3.5 ms; field of view 256 × 200 × 170 mm; matrix 256 × 200; voxel size 1.0 mm3; 170 axial slices; no fat suppression; scanning time ~ 6.5 min), to facilitate anatomical localization of activation for each participant. The duration of each functional run was 8 min for the resting task and 2 min for the other three tasks.

Functional MRI preprocessing

The first three volumes from each functional run were discarded for all trials. The fMRI data collected were preprocessed using the AFNI software package. Preprocessing steps performed on the whole brain included despiking, slice time correction, and 3-D isotropic reslicing (3 mm in each dimension). Motion correction using rigid body alignment was performed to correct for any major head motion during scanning. As the brainstem can move independently from the rest of the brain, a separate motion correction of the brainstem was performed. The amount of head motion was assessed using mean volume-to-volume framewise displacement (Power et al., 2012). Fifty-seven regions of interest (ROIs) were defined using FreeSurfer software (Harvard, MA, USA), and HMAT atlas for motor areas (Mayka et al., 2006), from the T1-weighted scans. Each participant's structural scan was then co-registered to the fMRI scan using rigid registration.

All analyses were done in the individual’s native space rather than transforming all data to a common template, to prevent excessive error (Nieto-Castanonet al., 2003, Ozcanet al., 2005), particularly in small subcortical regions such as the basal ganglia, as we have specifically shown in individuals with Parkinson’s disease (Chen et al., 2009, Nget al., 2009). Nuisance regression was then used to remove several sources of variance such as head motion parameters, their temporal derivatives and their squares, white-matter signal, and the cerebrospinal fluid signal. The fMRI signal was then detrended and iteratively smoothed until it reached 6 full width at half maximum of smoothness. Finally, the fMRI signal was high-pass filtered using a 0.01 Hz filter cut-off. The preprocessed time courses of the voxels within each ROI were averaged as the overall activity of the ROI.

Effective brain connectivity analysis

Effective connectivity, which captures the causal and dynamic influence brain regions exert over one another and reveals the strength and directionality of information flow between brain regions (Appel-Cresswellet al., 2010, Friston, 2011), was computed between the 57 FreeSurfer-derived ROIs using the PCfdr (Peter Spirtes and Clark Glymour, false discovery rate) algorithm (Li et al., 2008a, Li and Wang, 2009). This algorithm is specifically designed to assess connectivity when a relatively large number of ROIs and relatively few time points are available. It explores the dependence/independence between pairs of ROIs conditional on an exhaustive search of all other ROI combinations. The FDR threshold was set to 5%. We pooled the Parkinson’s disease and control groups together and computed the significant connections amongst ROIs for all tasks. The strength and direction of significant connections between ROIs were determined using a dynamic Bayesian network (DBN) group analysis approach (Liet al., 2006, Liet al., 2007, Li et al., 2008b).

Four separate statistical models were created utilizing binomial logistic regression: contrasting SB vs {rest, proprioception} and DB vs {rest, proprioception} for both individuals with Parkinson’s disease and control groups. Principal component analysis (PCA) was employed to achieve dimensionality reduction for the significant connections detected by the PCfdr/DBN method. PCA was performed separately for controls and individuals with Parkinson’s disease. The first 12 principal components (PCs), explaining at least 75% of the variance in both groups, were selected. Overall model fits were determined by comparing the full model to an intercept only model using the log-likelihood test.

Only the DB vs {rest, proprioception} models for controls and individuals with Parkinson’s disease were statistically significant (see RESULTS), and thus these models were investigated further. The β regression coefficients and their Wald χ2 statistic were examined to determine which PCs were making a significant contribution to the prediction of the task. Connections with the largest effect on each significant PC were identified by converting the PCA loadings into Z-scores and selecting all connections with an absolute Z-score ≥ 1.96. To evaluate the ability of the logistic regression models to correctly predict the task category of observed cases, contingency tables and leave-one-out cross validation were used, followed by two-tailed Fisher’s Exact Tests of independence to examine the relationship between performed and predicted tasks.

In order to directly compare relevant connections between controls and individuals with Parkinson’s disease, we first took the union of the connections significant in the DB task in either individuals with Parkinson’s disease or controls. We then performed logistic regression with Least Absolute Shrinkage and Selection Operator (LASSO) regularization to identify connections predicting group membership. An α < 0.05 was used for all statistical comparisons.

Participant baseline characteristics and head motion parameters

For both participant baseline characteristics and head motion parameters, assumptions of normality were validated using Shapiro-Wilk’s test and inspection of histograms and quantile–quantile plots. Baseline characteristics between controls and individuals with Parkinson’s disease were compared using independent t-tests or Fisher’s Exact Tests where appropriate. Mean volume-to-volume framewise displacements, used to assess head motion, were log-transformed due to non-normality and subsequently analyzed using a 2 × 4 mixed design ANOVA with group (Parkinson’s disease, controls) and task (rest, proprioception, SB, DB) as independent variables. In the case where Mauchly’s tests of sphericity (p < 0.05) was significant, the Greenhouse-Geisser ε statistic was used. An α < 0.05 was used for all statistical comparisons. In case of significant main and interaction effects, post hoc comparisons were performed after adjusting for multiple comparisons using a Bonferroni correction.

Data availability

De-identified data are available from the corresponding author upon reasonable request.

Results

Head motion

Minimal head motion occurred during the MRI scans; for both groups mean volume-to-volume framewise displacement was <0.3 mm in all tasks (Table 3). This is comparable to head motion reported previously during resting-state fMRI in older adults (Li et al., 2017) and fMRI during a foot-pedalling task in individuals with Parkinson’s disease (Matar et al., 2019). No significant main effects of task or group, or interaction effects, were found for mean framewise displacement.

Table 3

Mean volume-to-volume framewise displacement.

Task	Parkinson	Control
Rest	0.16 (0.01)	0.19 (0.02)
Proprioception	0.26 (0.05)	0.17 (0.02)
Static balancing	0.16 (0.02)	0.17 (0.02)
Dynamic balancing	0.27 (0.03)	0.23 (0.02)

Data are displayed as mean (SE)

Effective connectivity

The PCfdr/DBN method detected 164 significant connections between ROIs, which represent ~5.1% of all possible (57 × 56 = 3,192) directional connections. PCA of these 164 significant connections showed that the first 12 PCs contained 76% and 90% of the total variability for controls and individuals with Parkinson’s disease respectively.

Static balancing task

In both controls and individuals with Parkinson’s disease, the full model with 12 predictors was not significantly better at predicting the tasks compared to the constant only model (controls: χ2[12, N = 51] = 18.632, p = 0.098; individuals with Parkinson’s disease: χ2[12, N = 51] = 10.706, p = 0.554).

Dynamic balancing task

In controls, the log-likelihood test indicated that the full model with 12 predictors was significantly better at predicting the tasks compared to the constant only model (χ2[12, N = 51] = 43.286, p = 2e-5). The full model was able to correctly classify 97% of the resting/proprioceptive tasks and 94% of the DB tasks, for an overall success rate of 96%. Fisher’s Exact Test indicated there was a statistically significant relationship between the performed and predicted tasks (p < 0.0001, odds ratio [95% asymptotic confidence interval] = 528.000 [30.991, 8995.500]). During leave-one-out cross validation, the full model was able to correctly classify 88% of the resting/proprioceptive tasks and 47% of the DB tasks, for an overall success rate of 75%. There was a statistically significant relationship between the performed and predicted tasks (p = 0.012, 6.667 [1.623, 27.377]). Examination of the β coefficients and their Wald χ2 statistic showed 4 PCs making a significant contribution to the prediction of the DB task (Table 4). In general, the connections with the largest effect on the significantly contributing PCs connect the following areas (number of corresponding connections in brackets): brainstem to either cerebellum (2), limbic (1), thalamus (1), or other brainstem structures (1); limbic to other limbic (2) or basal ganglia structures (2); motor to other motor areas (4); frontal to other frontal (2) or basal ganglia structures (1); temporal to frontal (4) or limbic areas (3); and parietal to frontal (1) or temporal areas (1) (Table 5, Fig. 2A and 2B).

Table 4

Multiple binomial logistic regression models.

Predictor	β	SE β	Wald’s χ2	df	p	OR
Elderly controls - DB task (H-DB)
Constant	−3.023	1.259	5.765	1	0.016	0.049
H-DB PC-1	0.236	0.277	0.723	1	0.395	1.266
H-DB PC-2	0.457	0.531	0.740	1	0.390	1.579
H-DB PC-3	0.814	0.565	2.078	1	0.149	2.258
H-DB PC-4	−0.078	0.622	0.016	1	0.901	0.925
H-DB PC-5	−0.452	0.709	0.405	1	0.524	0.637
H-DB PC-6	1.954	1.061	3.390	1	0.066	7.058
H-DB PC-7	2.600	1.077	5.829	1	0.016	13.462
H-DB PC-8	2.046	1.383	2.189	1	0.139	7.735
H-DB PC-9	−6.384	2.719	5.512	1	0.019	0.002
H-DB PC-10	0.047	0.952	0.002	1	0.961	1.048
H-DB PC-11	−6.428	2.646	5.902	1	0.015	0.002
H-DB PC-12	−4.827	2.288	4.450	1	0.035	0.008
Individuals with Parkinson's disease - DB task (P-DB)
Constant	−3.363	1.047	10.327	1	0.001	0.035
P-DB PC-1	−1.300	0.549	5.610	1	0.018	0.272
P-DB PC-2	−1.545	0.504	9.383	1	0.002	0.213
P-DB PC-3	1.061	0.446	5.658	1	0.017	2.890
P-DB PC-4	−0.236	0.293	0.650	1	0.420	0.790
P-DB PC-5	0.185	0.494	0.141	1	0.708	1.203
P-DB PC-6	1.514	0.860	3.097	1	0.078	4.545
P-DB PC-7	0.991	1.090	0.826	1	0.364	2.693
P-DB PC-8	8.977	2.834	10.031	1	0.002	7917.000
P-DB PC-9	0.802	1.528	0.275	1	0.600	2.230
P-DB PC-10	0.670	1.019	0.431	1	0.511	1.953
P-DB PC-11	−2.996	1.524	3.866	1	0.049	0.050
P-DB PC-12	−2.750	1.448	3.606	1	0.058	0.064

β = beta coefficient, SE β = standard error of beta coefficient, Wald’s χ2 = Wald’s chi-square statistic, df = degrees of freedom, p = significance Wald’s χ2, OR = odds ratio (equal to exp(β)). Bolded predictors made a significant contribution to the prediction of the task.

Table 5

Connections with largest effect on the significant principal components predicting dynamic balance task in controls.

PC	Connection from	Connection to	From area	To area	Loading	Z-score	Strength RL/PL	Strength DL	Δ Strength	Direction
PC-7
	Pons	L cerebellum	BS	CER	0.455	5.810	0.376	0.453	−0.077	↑
	R superior frontal ctx	R dorsolateral prefrontal ctx	FRON	FRON	0.191	2.448	0.286	0.543	−0.257	↑
	Midbrain	L insular ctx	BS	LIMB	−0.206	−2.630	0.607	0.993	−0.386	↑
	Midbrain	L thalamus	BS	THAL	−0.246	−3.141	0.531	0.780	−0.249	↑
	L superior temporal ctx	L amygdala	TEMP	LIMB	−0.257	−3.282	0.392	0.398	−0.006	↑
	L middle temporal ctx	L amygdala	TEMP	LIMB	0.397	5.078	0.346	0.335	0.010	↓
	Midbrain	L cerebellum	BS	CER	−0.181	−2.310	0.358	0.265	0.093	↓
	R pre-SMA	L pre-SMA	MOT	MOT	−0.184	−2.347	0.584	0.469	0.115	↓
PC-9
	Pons	L cerebellum	BS	CER	0.400	5.170	0.376	0.453	−0.077	↑
	R anterior cingulate ctx	R caudate	LIMB	BG	−0.165	−2.059	0.201	0.489	−0.288	↑
	R inferior temporal ctx	R lateral orbitofrontal ctx	TEMP	FRON	−0.189	−2.363	0.178	0.297	−0.119	↑
	R middle temporal ctx	R lateral orbitofrontal ctx	TEMP	FRON	0.352	4.549	0.393	0.375	0.019	↓
	R dorsal premotor area	R SMA	MOT	MOT	0.168	2.204	0.393	0.295	0.098	↓
	L inferior parietal ctx	L middle temporal ctx	PAR	TEMP	0.166	2.171	0.310	0.264	0.046	↓
	L inferior parietal ctx	L dorsolateral prefrontal ctx	PAR	FRON	0.163	2.130	0.247	0.222	0.025	↓
	Midbrain	L cerebellum	BS	CER	−0.238	−2.997	0.358	0.265	0.093	↓
PC-11
	R superior temporal ctx	L amygdala	TEMP	LIMB	0.300	3.794	0.317	0.384	−0.067	↑
	R lateral orbitofrontal ctx	R medial orbitofrontal ctx	FRON	FRON	0.242	3.054	0.530	0.574	−0.045	↑
	L middle temporal ctx	L lateral orbitofrontal ctx	TEMP	FRON	−0.163	−2.121	0.201	0.310	−0.109	↑
	R amygdala	L amygdala	LIMB	LIMB	−0.163	−2.124	0.457	0.478	−0.021	↑
	L anterior cingulate ctx	L caudate	LIMB	BG	−0.208	−2.693	0.168	0.304	−0.135	↑
	R anterior cingulate ctx	R caudate	LIMB	BG	−0.221	−2.862	0.201	0.489	−0.288	↑
	L SMA	L pre-SMA	MOT	MOT	−0.226	−2.921	0.209	0.391	−0.182	↑
	L dorsal premotor area	L pre-SMA	MOT	MOT	0.305	3.855	0.445	0.441	0.005	↓
	Pons	Midbrain	BS	BS	0.211	2.655	0.302	0.206	0.096	↓
	L middle temporal ctx	L amygdala	TEMP	LIMB	−0.240	−3.109	0.346	0.335	0.010	↓
	Midbrain	L cerebellum	BS	CER	−0.261	−3.379	0.358	0.265	0.093	↓
PC-12
	R lateral orbitofrontal ctx	R medial orbitofrontal ctx	FRON	FRON	0.247	3.073	0.530	0.574	−0.045	↑
	R superior temporal ctx	L amygdala	TEMP	LIMB	0.212	2.625	0.317	0.384	−0.067	↑
	R insular ctx	L insular ctx	LIMB	LIMB	0.172	2.113	0.546	0.572	−0.027	↑
	R middle temporal ctx	R medial orbitofrontal ctx	TEMP	FRON	−0.147	−1.981	0.200	0.307	−0.107	↑
	Pons	L cerebellum	BS	CER	−0.160	−2.149	0.376	0.453	−0.077	↑
	R inferior temporal ctx	R lateral orbitofrontal ctx	TEMP	FRON	−0.193	−2.571	0.178	0.297	−0.119	↑
	R dorsolateral prefrontal ctx	R caudate	FRON	BG	−0.196	−2.614	0.353	0.357	−0.004	↑
	L superior temporal ctx	L amygdala	TEMP	LIMB	−0.245	−3.233	0.392	0.398	−0.006	↑
	Midbrain	L cerebellum	BS	CER	0.333	4.176	0.358	0.265	0.093	↓
	R middle temporal ctx	R lateral orbitofrontal ctx	TEMP	FRON	0.268	3.347	0.393	0.375	0.019	↓
	R dorsal premotor area	R SMA	MOT	MOT	0.214	2.650	0.393	0.295	0.098	↓
	L inferior parietal ctx	L middle temporal ctx	PAR	TEMP	−0.151	−2.036	0.310	0.264	0.046	↓

PC = principal component, R = right, L = left, ctx = cortex, SMA = supplementary motor area, BS = brainstem, CER = cerebellum, TEMP = temporal, LIMB = limbic, FRON = frontal, MOT = motor, THAL = thalamus, PAR = parietal

Fig. 2

Effective connectivity network found in healthy older adults during the dynamic balancing task.

In individuals with Parkinson’s disease, the log-likelihood test showed that the full model was significantly better at predicting the outcomes compared to the constant only model (χ2[12, N = 51] = 46.982, p = 5e-6). The full model correctly classified 97% of the resting/proprioceptive tasks and 88% of the DB tasks, for an overall success rate of 94%. Fisher’s Exact Test indicated there was a statistically significant relationship between the performed and predicted tasks (p < 1e-4, 247.500 [20.792, 2946.100]). During leave-one-out cross validation, the full model was able to correctly classify 88% of the resting/proprioceptive tasks and 71% of the DB tasks, for an overall success rate of 82%. There was a statistically significant relationship between the performed and predicted tasks (p = 4e-5, 18.000 [4.116, 78.711]). Examination of the β coefficients and their Wald χ2 statistic showed 5 PCs making a significant contribution to the prediction of the DB task (Table 4). In general, the connections with the largest effect on the significantly contributing PCs connect the following areas (number of corresponding connections in brackets): brainstem to either cerebellum (2), limbic (2), thalamus (1), or basal ganglia structures (2); limbic to other limbic (2) or basal ganglia structures (1); motor to other motor (2) or frontal areas (1); frontal to other frontal (2) or basal ganglia structures (1); temporal to either frontal (2), limbic (2), or other temporal areas (1); and parietal to other parietal areas (1) (Table 6, Fig. 3A and 3B).

Table 6

Connections with largest effect on the significant principal components predicting dynamic balance task in individuals with Parkinson’s disease.

PC	Connection from	Connection to	From area	Toarea	Loading	Z-score	Strength RL/PL	Strength DL	Δ Strength	Direction
PC-1
	Midbrain	L amygdala	BS	LIMB	0.987	12.566	2.600	1.133	1.467	↓
PC-2
	Midbrain	R pallidum	BS	BG	−0.173	−2.269	1.147	2.194	−1.047	↑
	Midbrain	L nucleus accumbens	BS	BG	0.947	12.051	1.716	0.960	0.757	↓
PC-3
	Midbrain	R pallidum	BS	BG	0.923	11.764	1.147	2.194	−1.047	↑
	Pons	R insular ctx	BS	LIMB	0.222	2.725	0.969	1.206	−0.237	↑
	Midbrain	L nucleus accumbens	BS	BG	0.181	2.189	1.716	0.960	0.757	↓
PC-8
	L superior temporal ctx	L amygdala	TEMP	LIMB	0.299	3.722	0.371	0.486	−0.115	↑
	R lateral orbitofrontal ctx	R medial orbitofrontal ctx	FRON	FRON	0.231	2.851	0.507	0.578	−0.071	↑
	R middle temporal ctx	R lateral orbitofrontal ctx	TEMP	FRON	0.176	2.134	0.273	0.498	−0.225	↑
	R anterior cingulate ctx	L anterior cingulate ctx	LIMB	LIMB	0.172	2.088	0.578	0.790	−0.212	↑
	R superior frontal ctx	R dorsolateral prefrontal ctx	FRON	FRON	0.165	1.996	0.379	0.476	−0.097	↑
	L middle temporal ctx	L amygdala	TEMP	LIMB	−0.259	−3.464	0.314	0.400	−0.086	↑
	R middle temporal ctx	R medial orbitofrontal ctx	TEMP	FRON	−0.259	−3.464	0.308	0.346	−0.038	↑
	Midbrain	L cerebellum	BS	CER	0.197	2.404	0.465	0.426	0.039	↓
	Midbrain	L thalamus	BS	THAL	0.190	2.323	0.786	0.677	0.109	↓
	R dorsal premotor area	R dorsolateral prefrontal ctx	MOT	FRON	−0.170	−2.319	0.373	0.310	0.064	↓
	R insular ctx	R amygdala	LIMB	LIMB	−0.179	−2.432	0.688	0.584	0.105	↓
PC-11
	Pons	L cerebellum	BS	CER	0.287	3.716	0.399	0.722	−0.323	↑
	L dorsal premotor area	L pre-SMA	MOT	MOT	0.175	2.282	0.449	0.788	−0.339	↑
	R superior parietal ctx	R precuneus	PAR	PAR	0.161	2.096	0.216	0.250	−0.033	↑
	R middle temporal ctx	L middle temporal ctx	TEMP	TEMP	−0.170	−2.123	0.340	0.384	−0.043	↑
	R anterior cingulate ctx	R caudate	LIMB	BG	−0.198	−2.489	0.183	0.410	−0.226	↑
	L SMA	L pre-SMA	MOT	MOT	−0.213	−2.676	0.239	0.304	−0.065	↑
	Midbrain	L cerebellum	BS	CER	0.434	5.593	0.465	0.426	0.039	↓
	R dorsolateral prefrontal ctx	R caudate	FRON	BG	0.213	2.771	0.311	0.245	0.066	↓
	R insular ctx	R amygdala	LIMB	LIMB	0.176	2.287	0.688	0.584	0.105	↓
	Midbrain	L thalamus	BS	THAL	0.157	2.050	0.786	0.677	0.109	↓

PC = principal component, R = right, L = left, ctx = cortex, SMA = supplementary motor area, BS = brainstem, BG = basal ganglia, CER = cerebellum, TEMP = temporal, LIMB = limbic, FRON = frontal, MOT = motor, THAL = thalamus, PAR = parietal

Fig. 3

Effective connectivity network found in individuals with Parkinson’s disease during the dynamic balancing task.

To directly compare between groups, 34 connections, making a significant contribution to predicting the DB task in controls and/or individuals with Parkinson’s disease, were entered in the logistic LASSO regression. Five connections survived, indicating that they are important to predicting group membership (control or individual with Parkinson’s disease) during the DB task (Table 7). The connections predicting group membership involve the following areas: midbrain to left nucleus accumbens, left anterior cingulate cortex to left caudate, right insular cortex to left insular cortex, left dorsal premotor area to left pre-supplementary motor area, and right middle temporal cortex to right lateral orbitofrontal cortex.

Table 7

Multiple binomial logistic LASSO regression model predicting group membership in the dynamic balance task.

	Predictor	Connection from	Connection to	From	To	β
DB	Constant					−0483
	Conn-1	L dorsal premotor area	L pre-SMA	MOT	MOT	1.246
	Conn-2	R middle temporal ctx	R lateral orbitofrontal ctx	TEMP	FRON	0.618
	Conn-3	Midbrain	L nucleus accumbens	BS	BG	−0.067
	Conn-4	L anterior cingulate ctx	L caudate	LIMB	BG	−0.629
	Conn-5	R insular ctx	L insular ctx	LIMB	LIMB	−0.640

DB = dynamic balance, R = right, L = left, ctx = cortex, SMA = supplementary motor area, BS = brainstem, BG = basal ganglia, TEMP = temporal, LIMB = limbic, MOT = motor, β = beta coefficient.

Discussion

We used a novel MRI-compatible balance simulator (Pasman et al., 2019) to investigate effective brain connectivity patterns related to static and dynamic balance-related tasks in individuals with Parkinson’s disease and healthy older adults. The simulator provides a unique opportunity to study brain activation in participants during balance-like movements in the MRI scanner, as it requires the static and dynamic reactive control of an unstable inverted pendulum system through the activation of anti-gravity muscles around the ankle joint. This ensures that both motor and proprioceptive sensory pathways are engaged during the task in a way that is similar to upright standing, unlike prior studies that relied on motor imagery of balance. We showed that a network of cortical and subcortical neural structures, including frontal, parietal, and temporal cortices as well as basal ganglia, cerebellum, and brainstem, was active during dynamic balancing in healthy older adults. In individuals with Parkinson’s disease, a similar network of cortical and subcortical neural structures was active during dynamic balancing, however, the strength and direction of connections between certain areas were different.

Effective connectivity network for dynamic balancing in healthy older adults

We observed decreased connectivity between different motor areas and increased connectivity from the brainstem to several cortical and subcortical areas in healthy older adults during dynamic balancing compared to the rest/proprioception conditions. Taken together, this indicates a preference of subcortical over motor cortical control of dynamic balancing in healthy older adults. Our results are generally consistent with prior results from motor imagery studies that identified several subcortical (e.g., mesencephalic locomotor region, thalamus, cerebellum, and basal ganglia) and cortical (e.g., frontal, temporal, parietal, and cingulate cortices) areas as part of a distributed dynamic balance network (Ferrayeet al., 2014, Taubeet al., 2015, Bhattet al., 2018, Mouthonet al., 2018). Additionally, our results provide novel insight into the directionality of the connections between these areas and highlighting the potential context- and task-specificity of supplementary motor area (SMA) activity during dynamic balance tasks.

Connectivity involving subcortical and brainstem regions increases in healthy older adults during dynamic balancing

As shown in Fig. 2B, connectivity increased from the dorsolateral prefrontal and anterior cingulate cortices to the basal ganglia, more specifically the caudate nucleus, in healthy elderly during dynamic balancing. This finding is consistent with fMRI evidence of increased activation in the prefrontal and anterior cingulate cortices, as well as in the basal ganglia (putamen and globus pallidus), during motor imagery of dynamic balance tasks (Ferrayeet al., 2014, Taubeet al., 2015, Bhattet al., 2018, Mouthonet al., 2018). Increased activity in prefrontal and anterior cingulate cortices during dynamic balance responses has also been implicated with EEG and fNIRS (Miharaet al., 2008, Sippet al., 2013, Bogostet al., 2016, Solis-Escalanteet al., 2019).

The involvement of the prefrontal cortex during dynamic balancing that we observed is consistent with its role in error detection (Bogost et al., 2016) and monitoring postural stability (Solis-Escalante et al., 2019) through sensory integration and allocation of attentional resources (Miharaet al., 2008, Teoet al., 2018). The anterior cingulate also plays an important role in monitoring motor error which is critical for detecting loss of balance during challenging posture and gait tasks (Sippet al., 2013, Marlinet al., 2014, Bhattet al., 2018, Goelet al., 2019). The basal ganglia have long been linked to balance control, particularly in terms of sensorimotor integration, gain control of balance correcting responses, and proper selection and execution of context-specific balance correcting responses (Visser and Bloem, 2005, Jacobs and Horak, 2007, Ferrayeet al., 2014, Taubeet al., 2015). Therefore, the observed increase in connectivity from the prefrontal and anterior cingulate cortices to the basal ganglia during our simulated dynamic balance task adds further evidence for the critical role of these areas in maintaining balance control.

We also found increased connectivity from pons to cerebellum, and midbrain to thalamus and insular cortex (Fig. 2B). These results are consistent with prior findings from motor imagery of dynamic balance tasks that have demonstrated increased activation of the pons and mesencephalic locomotor region of the midbrain (Ferraye et al., 2014), and increased activation of the cerebellum, thalamus, and insula (Ferrayeet al., 2014, Taubeet al., 2015, Bhattet al., 2018). Increased connectivity between brainstem structures, the thalamus, and cerebellum supports growing evidence for the importance of the corticopontocerebellar-thalamocortical loop in integrating feedforward commands with sensory information to shape and modulate subsequent motor commands (MacKinnon, 2018) and allows for postural responses to be modified by prior experience and changes in central set (Jacobs and Horak, 2007, Bhattet al., 2018).

Several areas in the pons are thought to be involved in balance control, including the pontomedullary reticular formation, locus coeruleus, and vestibular nuclei (Takakusaki, 2017, MacKinnon, 2018). The vestibular nuclei contribute to the modulation of dynamic balance responses through descending vestibular-spinal pathways, ascending projections to the thalamus and insula, and bidirectional connections to the vermis and flocculonodular nuclei of the cerebellum (Takakusaki, 2017, MacKinnon, 2018). The increased connectivity observed during dynamic balancing from pons to cerebellum, supports extensive evidence for the vestibulo-cerebellar networks crucial role in maintaining postural control (Carpenteret al., 2001, Morton and Bastian, 2004, Schnieppet al., 2017), allowing for efferent and afferent proprioceptive and vestibular inputs to be compared in order to rapidly identify and correct movement error associated with a balance disturbance (Peterka, 2018). Nuclei within the midbrain have also been linked to postural control, including the mesencephalic locomotor region (composed of the cuneiform nucleus and pedunculopontine nucleus) and caudal raphe nuclei (Takakusaki, 2017, MacKinnon, 2018). The pedunculopontine nucleus has connections with the basal ganglia and limbic areas, thalamus, cerebellum, brainstem, spinal cord, and cerebral cortex (Alam et al., 2011). The observed connection from midbrain to thalamus makes sense given the role of the posterolateral thalamus in maintaining upright posture (Masdeu and Gorelick, 1988), and verticality perception (Jacobs and Horak, 2007, Barraet al., 2010). The parieto-insular vestibular cortex, which includes part of the posterior insula and parietal operculum, as well as part of the temporo-perisylvian cortex in the superior temporal gyrus, has strong connections with other vestibular-related cortical areas, and receives converging sensory inputs (Indovinaet al., 2015, Takakusaki, 2017). This region is thought to be essential for sensory integration during postural tasks (Jacobs & Horak, 2007), contributing to perception of visual and perceived gravitational vertical (Jacobs and Horak, 2007, Takakusaki, 2017) and spatial orientation and self-motion perception (Taubeet al., 2015, Bhattet al., 2018).

Healthy older adults also showed increased connectivity from temporal cortices to both subcortical and other cortical areas (Fig. 2A and 2B). Connectivity increased from the superior temporal cortex to the amygdala, as well as middle temporal cortex and inferior temporal cortex to the lateral and medial orbitofrontal cortices. The superior temporal gyrus was found to be activated in previous studies during motor imagery of dynamic balancing (Taube et al., 2015), and its temporo-perisylvian cortex is part of the parieto-insular vestibular cortex (Indovina et al., 2015).

Connectivity involving motor and parietal regions decreases in healthy older adults during dynamic balancing

Our results indicate that healthy older adults, when performing simulated dynamic balance, had decreased connectivity from the dorsal premotor area to SMA and pre-SMA, as well as the right to left pre-SMA (Fig. 2A). Our observation of reduced connectivity between and within motor cortical areas is in contrast to prior evidence of increased premotor and SMA activity recorded with MRI during motor imagery of dynamic balance tasks (Ferrayeet al., 2014, Taubeet al., 2015, Bhattet al., 2018, Mouthonet al., 2018). This discrepancy between the current study and prior motor imagery studies could be due to differences in the level of automaticity associated with the dynamic balance tasks. The dynamic balance task used in the current study resembles an externally triggered movement, that required rapid ankle responses to transient, unpredictable perturbations in random directions, that have components considered to be more automatic and reflexive (Carpenter et al., 2001). Reduced activation of pre-SMA, premotor areas, and superior and inferior parietal areas has been reported in healthy older adults as tasks become more automatic (Wu and Hallett, 2005a). In contrast to the dynamic balance-simulator task used in the current study, Ferraye et al. (2014) used a voluntary sway task which requires feedforward coordination of postural adjustments to the voluntary sway movements. In addition, internally and externally triggered movements are associated with different brain activation patterns, with increased activation in SMA during internally driven compared to externally triggered movements (Filyushkina et al., 2019). The SMA and premotor areas are most critical in preparing planned or complex motor programs (Takakusaki, 2017, MacKinnon, 2018) and thus more likely to contribute to balance tasks that require more feedforward, less automatic balance control (Ferrayeet al., 2014, Mierauet al., 2015, Takakuraet al., 2015), and/or later phases of the compensatory balance reactions such as stepping or reaching responses (Miharaet al., 2008, Marlinet al., 2014).

Dynamic balancing relies more on cortical than subcortical control in individuals with Parkinson’s disease

When individuals with Parkinson’s disease were asked to perform a dynamic balancing task, compared to the rest/proprioception tasks, they showed an effective connectivity pattern largely similar to healthy older adults. However, there were some notable differences. Individuals with Parkinson’s disease had increased connectivity associated with motor control and parietal areas, and decreased connectivity from brainstem to other subcortical areas, suggesting that dynamic balance control in individuals with Parkinson’s disease relies relatively more on cortical motor areas. The increased cortical involvement we observed in individuals with Parkinson’s disease during dynamic balancing compared to healthy older adults is consistent with the ‘posture second strategy’ observed in individuals with Parkinson’s disease, whereby engagement in a mental task (e.g., mental arithmetic) results in disproportionate decreases in postural control (Bloem et al., 2006).

Cortical connectivity was generally increased in individuals with Parkinson’s disease during dynamic balancing

In contrast to healthy older adults who showed decreased connectivity associated with motor control areas and the inferior parietal cortex (Fig. 2A), individuals with Parkinson’s disease had increased connectivity from the SMA and premotor areas to the pre-SMA, and from the superior parietal cortex to precuneus (Fig. 3A). Increased activation in the pre-SMA and SMA was previously found in individuals with Parkinson’s disease ‘on’ medication, compared to healthy older adults, during active, but not passive, ankle movements (Katschnig et al., 2011). Our findings are also in line with previous observations of reduced activation in motor control and parietal areas in healthy older adults (Wu and Hallett, 2005a), but increased activation in the cerebellum, motor, parietal, and prefrontal areas in individuals with Parkinson’s disease ‘off’ medication (Wu and Hallett, 2005b), during automatic movements. Therefore, it seems as if individuals with Parkinson’s disease ‘off’ medication compensate for their inefficient brain activity with increased premotor-parietal region activation during execution of automatic movements (Wu and Hallett, 2005b). The qualitative differences in connectivity associated with motor control areas were supported with direct statistical comparison between the groups. In particular, individuals with Parkinson’s disease were best discriminated from the healthy older adults by an abnormal increase in connectivity from the dorsal premotor area to the pre-SMA.

While cortical connectivity increased in motor and parietal areas, connectivity from the dorsal prefrontal and anterior cingulate cortices to the caudate nucleus were reduced in individuals with Parkinson’s disease compared to healthy older adults. Reduced activity in the superior parietal and anterior cingulate cortices was found previously during motor imagery of gait in individuals with Parkinson’s disease ‘off’ medication compared to healthy older adults (Snijders et al., 2011). Direct statistical comparison between the groups confirmed that a lack of normal increase in connectivity strength from the anterior cingulate cortex to the caudate nucleus was important in discriminating individuals with Parkinson’s disease from healthy older adults. In addition, an abnormal increase in connectivity from the middle temporal cortex to lateral orbitofrontal cortex, and lack of normal increase between the bilateral insular cortices were also important in best discriminating individuals with Parkinson’s disease from health older adults during dynamic balancing.

Subcortical connectivity was altered in individuals with Parkinson’s disease during dynamic balancing.

Reduced connectivity was seen in individuals with Parkinson’s disease from the midbrain to thalamus, cerebellum, nucleus accumbens, and amygdala (Fig. 3B). These findings corroborate previous work which found that the tendency to shift execution of automatic movements to subcortical areas is less clear in individuals with Parkinson’s disease ‘off’ medication compared to healthy older adults (Wu et al., 2010). Some of the subcortical connections identified in predicting the dynamic balancing task in individuals with Parkinson’s disease, such as the connections from the midbrain to thalamus and cerebellum, were common to both groups. However, unique connections in individuals with Parkinson’s disease included reduced connectivity from midbrain to nucleus accumbens and amygdala, and increased connectivity from midbrain to the pallidum and pons to insular cortex. Direct statistical comparison between the groups confirmed an abnormal decrease in connectivity from the midbrain to nucleus accumbens was important in discriminating individuals with Parkinson’s disease from health older adults.

Effective connectivity network for static balancing

In both healthy older adults and individuals with Parkinson’s disease we were unable to determine an effective connectivity network associated with the static balance task, compared to rest and proprioception tasks. Our results suggest that the brain connectivity network associated with the simulated static balancing task might not have been distinct enough from that associated with the proprioception task, involves a more distributed network, or may be controlled using predominantly passive stiffness and spinal reflex mechanisms (Winter et al., 1998). Our current protocol prioritized proprioceptive feedback, since the task was performed with the eyes closed, and no head rotational or translational inputs would be detected by vestibular receptors due to participants lying supine with the head fixed. However, this does not suggest that the vestibular system is inactive during our simulated balance task, since otoliths would still be activated by gravity, and the vestibular nuclei also integrate non-labyrinthine inputs from cortical, subcortical, and spinal cord pathways, including proprioceptive inputs from the limbs (McCall et al., 2017). Although congruent sensory information is normally available from multiple sensory systems during upright stance, the simulated static balance task used in the current study is still relevant, since balance can also be maintained by the central nervous system even in the absence of reliable vestibular and visual inputs. This is due to redundancy and the ability of the central nervous system to re-weigh sensory information (Peterka, 2018). Furthermore, similar balance deficits with Parkinson’s disease have been observed during real and simulated static and dynamic balance tasks with and without vision (Pasman et al., 2019). Future studies could introduce more balance-relevant sensory inputs by incorporating real-time visual feedback (Pasman et al., 2019) or electrical vestibular stimulation to determine whether additional sources of sensory information provide more distinct connectivity models for static balance control.

Limitations and future directions

For both the static and dynamic balancing tasks, the connectivity of cortical regions involved in sensory integration, such as parietal areas, could have been affected by the limited availability of balance-relevant sensory information, other than proprioceptive feedback. Future studies should investigate how the addition of visual information (i.e., real-time visual feedback from a virtual reality scene driven by the output from a potentiometer attached to the axis of the simulator) and/or vestibular information (i.e., using electrical vestibular information) changes the connectivity of these particular regions during the simulated balance tasks.

While all participants were familiarized with the balance simulator in a session separate from the MRI scanning session, differences between individuals with Parkinson's disease and healthy older adults regarding motor learning and retention may have contributed to the difference in cortical and subcortical control between the groups during the novel simulated balance task (Marinelli et al., 2017).

The conscious nature of the proprioceptive (i.e., finger tracking) task used as the control task may have resulted in a downplay of the role of cortical cognitive centers during the simulated balance tasks.

Using the simulator in its current form, it is not possible to assess medial–lateral balance control. However, while different muscle/joints are involved in controlling anterior-posterior and medial–lateral postural sway, there is no evidence to suggest they are controlled using different cortical or subcortical structures.

The risk of head movements was increased due to the clinical population investigated and the tasks with lower leg motion used. During data collection, head motion was minimized by immobilizing the participant's head using straps and foam, as well as stabilizing participant's hips using weights (see METHODS). Head motion, as assessed by mean volume-to-volume framewise displacement, was comparable to head motion reported previously (Liet al., 2017, Mataret al., 2019).

The individuals with Parkinson’s disease included in this study consisted mostly of individuals with moderate disease (Table 2). Future work is needed to investigate whether the neural substrates associated with postural instability in individuals with Parkinson’s disease change with disease severity. In addition, included individuals with Parkinson’s disease also had relatively good cognition, which was required to comply with the task requirements. However, cognitive decline contributes to balance deficits in individuals with Parkinson’s disease, and excluding cognitively impaired individuals – as we and many others have done (Domingos et al., 2015) – may create an incorrect perspective of the severest balance deficits in Parkinson’s disease.

Future work is needed to correlate the neural substrates of balance control with the balance behaviour exhibited by participants when maintaining balance using the simulator or during upright standing. This provides the opportunity to investigate whether the neural substrates associated with postural instability in individuals with Parkinson’s disease change with the severity of their balance deficits as measured using either clinical balance or posturographic measures.

All individuals with Parkinson’s disease were tested during their subjectively best clinical ‘on’ condition as postural instability is often unaltered by dopaminergic medication (Grimbergen et al., 2009) and testing individuals with Parkinson's disease ‘on’ medication avoids any potential confounds of fatigue, anxiety, and cumbersome bradykinesia/rigidity that may accompany the ‘off’ phase. Future work is needed to investigate the effect of dopaminergic medication on the brain connectivity patterns associated with balance control by testing individuals with Parkinson’s disease during both the ‘on’ and ‘off’ condition, as previous studies have shown levodopa reduces the utilization of motor reserve for compensation in individuals with Parkinson’s disease (Palmer et al., 2009)

Conclusion

In conclusion, the results of this study indicate that a network of cortical and subcortical neural structures, including frontal, parietal, and temporal cortices as well as basal ganglia, thalamus, cerebellum, and brainstem, was active during dynamic balancing in healthy older adults. In individuals with Parkinson’s disease a similar network of cortical and subcortical neural structures was found. However, dynamic balancing was more reliant on motor cortical control in individuals with Parkinson’s disease compared to healthy older adults. The increased understanding of the neural substrates contributing to postural instability in Parkinson’s disease provided by the novel MRI-compatible balance simulator could lead to new targets for improved pharmacological and neurosurgical interventions.

CRediT authorship contribution statement

Elizabeth P. Pasman: Conceptualization, Project administration, Formal analysis, Investigation, Methodology, Visualization, Writing - original draft. Martin J. McKeown: Conceptualization, Project administration, Methodology, Resources, Writing - review & editing. Saurabh Garg: Software, Data curation, Formal analysis, Writing - review & editing. Taylor W. Cleworth: Investigation, Formal analysis, Writing - review & editing. Bastiaan R. Bloem: Conceptualization, Writing - review & editing. J. Timothy Inglis: Conceptualization, Writing - review & editing. Mark G. Carpenter: Conceptualization, Resources, Methodology, Investigation, Writing - review & editing, Funding acquisition, Supervision.

Declaration of Competing Interest

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