Abstract
Introduction:
There is accumulating evidence that many pathological conditions affecting human balance are consequence of postural control (PC) failure or overstimulation such as in motion sickness. Our research shows the potential of using the response to a complex postural control task to assess patients with early-stage Parkinson's Disease (PD).
Methods:
We developed a unique measurement model, where the PC task is triggered by a moving platform in a virtual reality environment while simultaneously recording EEG, EMG and CoP signals. This novel paradigm of assessment is called BioVRSea. We studied the interplay between biosignals and their differences in healthy subjects and with early-stage PD.
Results:
Despite the limited number of subjects (29 healthy and nine PD) the results of our work show significant differences in several biosignals features, demonstrating that the combined output of posturography, muscle activation and cortical response is capable of distinguishing healthy from pathological.
Discussion:
The differences measured following the end of the platform movement are remarkable, as the induced sway is different between the two groups and triggers statistically relevant cortical activities in α and θ bands. This is a first important step to develop a multi-metric signature able to quantify PC and distinguish healthy from pathological response.
1. Introduction
Parkinson's Disease (PD) is a progressive disorder of the nervous system characterized by muscle tremors, muscle rigidity, decreased mobility (bradykinesia), stooped posture, slow voluntary movements, and a mask-like facial expression. It may take time to diagnose because some of its symptoms are associated with the natural process of aging (Tolosa et al., 2006). Globally, disability and death in PD are increasing faster than any other neurological disorder. The World Health Organization (WHO) reports that the prevalence of PD has doubled in the past 25 years and world estimates count over 8.5 million individuals with PD in 2019. In people with early-onset PD, the initial symptoms can arise between the ages of 21 and 40 years, while the first symptoms in juvenile-onset disease occur before the age of 20 years. Nowadays, a standard criterion in the evaluation of PD is still one of the main goals for clinicians. Finding the right category for the progression of the disease is necessary to prescribe the best treatment. Specific signs, symptoms, or test results can help in the classification of the disease. Over the years, accuracy has been improved by new diagnostic protocols that consider qualitative and quantitative aspects (Maffoni et al., ; María et al., 2020). Defining early-stage Parkinson's subjects when the symptoms are silent or weak remains a challenge.
PD stages are identified based on clinical observations: according to the Hoehn-Yahr staging system, stages are based primarily on motor symptoms (Goetz et al., ). Pre-clinical is characterized by the absence of signs or symptoms - genetic testing and counseling are available to identify risk factors. Prodromal corresponds to a stage of neurodegenerative changes. Symptoms are unspecific but the identification of early changes allows to intervene with initial therapies. Early-stage symptoms include mild tremors and some walking difficulty. It can affect only one side of the body and produce a decrease of facial expressions. These symptoms do not interfere with daily life much and are not always obvious (Goetz et al., ). Middle-stage balance and coordination are affected: a moderate-to-severe disability that affects daily life (Goetz et al., ). Later-stage subjects have difficulty standing and walking even with aids. Patients in this stage have severe disability (Goetz et al., ). Figure 1 reports the general classification of PD stages according to its main symptoms.
Figure 1
A preliminary analysis can be carried out by clinicians taking into account qualitative and quantitative aspects of PD, as shown in Figure 2.
Figure 2
Gait mobility and gait impairment is also seen to evolve with the progression of the disease are included as parameters of investigation. Issues with gait initiation, freezing of gait, reduced balance, and difficulties in postural control (Martignon et al., 2021) are some of common symptoms. Table 1 summarizes the main tests used in mobility and gait analyses.
Table 1
| Assessment | Test | Parameter investigated |
|---|---|---|
| Balance | Timed Up and Go Test | Functional mobility |
| Tinetti Balance and Gait Test | Static and dynamic balance | |
| Retropulsion Test | Postural stability | |
| One-leg Stance | Static balance | |
| Åstrand-Rhyming protocol during Graded Exercise Test | Dynamic balance | |
| Balance Evaluation System Test | Balance systems | |
| Mini Balance Evaluation Systems Test | Dynamical balance | |
| Berg Balance Scale | Static and dynamic balance | |
| Endurance | 2-min Walk Test | Functional capacity, walking ability |
| 6-min Walk Test | Physical capacity and gait | |
| Two-minute Step Test | Aerobic capacity | |
| Modified Bruce Protocol during Graded exercise test | Cardiac functional capacity | |
| Åstrand-Rhyming protocol during Graded exercise test | Maximal functional capacity | |
| Borg Ratio Scale | Physical capacity | |
| Resistance | Handgrip Strength Test | Upper limb strength |
| Isokinetic Strength Test | Upper/lower limb strength | |
| Manual Muscle Test | Individual or grouped muscle strength | |
| Arm Curl Test | Upper limb strength | |
| Chair Stand Test | Lower limb strength | |
| Five Time Sit to Stand Test | Lower limb strength | |
| One Repetition Maximum Test | Maximum lower or upper limb isotonic | |
| Flexibility | Goniometer | General joint flexibility |
| Inclinometer | Angles of slope measurement | |
| Leighton flexometer | Joint flexibility | |
| Sit and Reach Test | Lower back and hamstring muscles tightness | |
| Back Scratch Test | Shoulder range of motion |
Gold standard tests for assessing gait ability.
Among the innovative techniques, blood tests show potential to be used for the detection of early-stage PD. Researchers identified a molecular profile that defines the disease but it is still under investigation and not yet available clinically (Agliardi et al.,
In recent years, non-invasive brain imaging techniques have become more accurate for the detection of differences in brain morphology and functional activities in Parkinson's subjects (Politis, 2014). Brain positron emission tomography (PET) can estimate the disease progression and can be used to confirm the clinical diagnosis of PD. With specific radioactive drugs (18F-DOPA and 18-FDG) absorbed into the bloodstream, PET can provide very precise brain region and activation in PD subjects (Lu and Yuan,
Map structure and functions of the brain are obtained measuring the signals produced by neural activity. Each region can have a particular influence according to the disease and the activation of an area can be considered important in the understanding of the progression of the disease. Although cortical EEG coherence can serve as a reliable measure of disease severity, the use of EEG to study PD has not been fully investigated. Neurophysiological signals provide instantaneous information and can aid in improving the accuracy of the diagnosis.
EEG signals have different specific frequency bands. Features in sub-bands are particularly important to characterize different brain states. The standard frequency bands of interest are δ-band (0–4 Hz), θ-band (4–8 Hz), α-band (8–13 Hz), and β-band (13–30 Hz). Moreover, the quantification of EEG rhythms could provide an important biomarker for different neuropsychiatric and neurological disorders, such as schizophrenia, Alzheimer's disease, epilepsy, and Parkinson's disease (Hampel et al.,
A recent study analyzed the EEG signals from 15 early-stage PD patients and 15 age-matched healthy controls during eyes-closed resting state (Han et al.,
Postural control (PC) and adaptation are part of a complex system to maintain or restore balance from any position or during motor activity. The central nervous system is fundamental in PC strategies and electroencephalography can underline the different cortical brain activities under different postural perturbations (Mochizuki et al., 2010; Barollo et al.,
Other studies show similar results in healthy subjects, suggesting cortex activity as the main change in the frontal-central and frontal-parietal cortical regions during balance perturbation, specifically within α and θ frequencies (Sipp et al., 2013; Hülsdünker et al.,
Figure 3

ASP bands analysis in different studies (Neufeld et al., 1988; Sinanović et al., 2005; Bosboom et al.,
Our novel BioVRSea setup introduces a unique multi-biometric system that combines virtual reality and a moving platform to evaluate the postural control response. The system is designed to imitate the sensation of being at sea on a small boat, a situation which involves different balancing strategies. During the experiment, there are six phases (see Table 2, Figure 4) in which different biosignals are measured such as electromyography (EMG), center of pressure (CoP), and electroencephalography (EEG). Some of our recent studies emphasize the importance of BioVRSea setup allowing cohort differentiation and pathology assessment (Recenti et al., 2021; Aubonnet et al.,
Table 2
| Time (s) | Segment | VR Scene | Position of Hands | Platform |
|---|---|---|---|---|
| 0-120 | Baseline | Mountains | By side | Stationary |
| 120-160 | PRE | Sea | By side | Stationary |
| 160-200 | 25% | Sea | On bars | Moving |
| 200–240 | 50% | Sea | On bars | Moving |
| 240-280 | 75% | Sea | On bars | Moving |
| 280-320 | POST | Sea | By side | Stationary |
BioVRSea experimental paradigm.
Figure 4

BioVRSea experimental setup.
In this paper, candidates with early-stage Parkinson's Disease undergo a postural control task in the BioVRSea environment. We focus on the differences between the two cohorts in the PRE and POST phases. Our focus on PRE and POST will allow the study of early stage PD while they are adapting to a motor stimulus.
2. Materials and methods
2.1. Participants
Nine early-stage (recently diagnosed) PD participants (6 male, 3 female, between 56 and 76 years of age) and 29 healthy subjects (17 male, 12 female, between 50 and 73 years old) took part in the BioVRSea experiment. Few of them showed physical evidence of early stage PD such as weak tremor or onset of postural instability. All were taking the drug Levodopa as part of their treatment.
2.2. BioVRSea experiment
A 64-channel wet EEG was used record brain response to VR and motion stimulation. Wireless EMG sensors were placed on the tibialis anterior (TA), gastrocnemius lateral (GL), and soleus (S) muscles of both legs. A heart rate sensor strapped around the chest. For the experiment, the participant were asked to stand onto the force plates embedded in the platform. Finally, the participants donned the VR goggles. The experimental protocol was then explained to the participant. Participants stood quietly on the platform with their hands by their side observing a mountain view for the first 2 minutes of the experiment (Baseline). Then, the scene in the VR goggles changed, beginning the sea simulation but no platform movement. The participants remained standing quietly with their hands by their side for the first 35 seconds of the sea simulation (PRE). After the PRE phase the platform began synchronized movement with the sea scene in the VR goggles, increasing from 25% to 75% of maximal wave amplitude. For a total of 120 seconds the participants held the bars of the platform while continuing to observe the sea simulation. Finally, the platform stopped moving while the sea simulation is still showing and the participant was asked to remove their hands from the bars and attempt to stand quietly with their hands by their side for the final 40 seconds of the experiment. This is called the POST phase of the experiment; it is performed identically to the PRE phase but after the participant has performed movement in the central part of the procedure. A table of the VR experiment protocol is shown below in Table 2, shows a schematic of the experimental setup. Each participant took part in a single trial according to the experimental protocol.The subject undergoes different stimuli: visual stimulus (PRE-phase), motor stimulus (movement phase), and balance control (POST-phase).Measured data was post-processed in Matlab and analysis was made in the PRE and POST phases of the experiment. Each analysis pipeline for a particular measurement is explained below.
2.2.1. Heart rate
Heart rate was measured using a chest heart sensor (Polar Electro, Kempele, Finland, sampling frequency 1 Hz). The average and standard deviation for the HR for each section of the experiment was calculated.
2.2.2. EMG analysis
Muscle electrical activities from the lower limbs were acquired using six wireless EMG sensors (sampling frequency of 1600 Hz) placed on the tibialis anterior (TA), gastrocnemius lateral (GL), and soleus (S) muscles of each leg (Kiso ehf, Reykjavik, Iceland). EMG data processing was performed using Matlab 2021b. EMG data were filtered using a 4th-order Butterworth filter. Seven features were computed in the frequency domain and thirty-six features in the time domain for each muscle and each phase of the experiment. These features are listed in Table 5.
2.2.2.1. Statistical analysis
The Shapiro-Wilk test along with visual inspection of the distribution of each variable were used to test the normality of the data. Statistical comparisons between the healthy and PD groups in both the PRE and POST phases were carried out using the t-test with Welch's correction for the normally distributed variables and the Mann-Whitney U-test for the non-normally distributed variables, with a significance value of p < 0.05. Effect sizes were calculated through the non-parametric Cliff's delta using the R package “effsize” (Torchiano, 2020). Cliff's delta ranges from +1 if all observations in the first group are larger than all observations in the second group, to -1 if all observations in the first group are smaller than all observations in the second group (Cliff,
2.2.3. CoP analysis
CoP measurements were made using 4 sensors located under each foot platform. The sensors give information about the center of mass in the Antero-Posterior and Medio-Lateral axis (Virtualis, Clapiers, France, sampling frequency 90Hz). The processing of the CoP data was performed using Matlab 2021b. During the experiment, the force platform records the movement of the Centre of Pressure (CoP), a projection of the center of mass of the subject on the plane of the machine, also called stabilogram. The CoP data was filtered with a Savitsky-Golay filter with window size 7. Included in the CoP analysis were a number of multi-scale entropy measurements, which have been shown to have great importance in the analysis of CoP data in discriminating between pathological subjects (Busa and Emmerik,
Table 3
| Feature | PRE p-value | POST p-value | Cliff delta |
|---|---|---|---|
| Consecutive movement samples on the support plane (TOTEX) | - | - | 0.379 |
| Consecutive movement samples on ML plane(TOTEX-ML) | - | - | 0.402 |
| Consecutive movement samples on AP plane(TOTEX-AP) | - | - | 0.371 |
| Square root distance between a point and the plane origin (RD) | - | 0.0173 | 0.494 |
| Mean distance in medio-lateral direction (MDIST-ML) | - | 0.0339 | 0.494 |
| Mean distance in antero-posterior direction (MDIST-AP) | - | 0.0173 | 0.494 |
| Mean velocity on support plane (MVELO) | - | - | 0.379 |
| Mean velocity on ML plane (MVELO-ML) | - | - | 0.402 |
| Mean velocity on AP plane (MVELO-AP) | - | - | 0.371 |
| Root mean square distance respect to origin (RDIST) | - | 0.0115 | 0.540 |
| Root mean square distance in medio-lateral direction (RDIST-ML) | - | 0.0256 | 0.517 |
| Root mean square distance in antero-posterior direction (RDIST-AP) | - | 0.0115 | 0.533 |
| Medio-lateral sample entropy (ML-SampEn) | - | 0.0002 | -0.709 |
| Antero-posterior Sample Entropy (AP-SampEn) | - | - | -0.333 |
| Medio-lateral complexity index (ML-CI) | - | 0.0009 | -0.793 |
| Antero-posterior Complexity Index (AP-CI) | - | - | -0.325 |
| Ellipse area | - | 0.0083 | 0.571 |
| Ellipse angle | - | - | 0.057 |
| Ellipse main axis length | - | 0.0067 | 0.571 |
| Ellipse minor axis length | - | 0.0127 | 0.571 |
| Standard deviation in antero-posterior direction (SD AP) | - | 0.0115 | 0.532 |
| Standard deviation in medio-lateral direction (SD ML) | - | 0.0257 | 0.517 |
| SD magnitude | - | 0.0074 | 0.563 |
| SD direction | 0.0128 | - | -0.256 |
| Magnitude entropy | - | - | -0.510 |
| Direction entropy | - | - | -0.249 |
| Multivariate complexity index (multivariate CI) | - | 0.0128 | -0.639 |
| Antero magnitude | - | - | 0.379 |
| Antero angle | - | - | 0.019 |
| Postero magnitude | - | 0.0141 | 0.548 |
| Postero angle | - | - | 0.065 |
| Left magnitude maximum | - | 0.0406 | 0.417 |
| Left angle | - | - | -0.065 |
| Right magnitude maximum | - | 0.0282 | 0.448 |
| Right angle | - | - | -0.494 |
CoP features calculated - bold shows the features that were significantly different in the POST phase of the experiment when comparing the PD and healthy groups.
Figure 5

Ellipse areas comparison between Parkinson and healthy subjects in PRE-POST phases.
2.2.4. EEG analysis
The CA-204-64 wet electrode cap, EegoTM mylab with sampling frequency of 4096 Hz, measured the brain electrical activity in 64 electrodes. Raw EEG signals were processed using Matlab 2022b, Brainstorm, EEGlab 2022.1 and Automagic toolboxes (Pedroni et al., 2019). The signals were divided in segments for each phase of the experiment, then the signals were down-sampled from 4,096 to 1,024 Hz. During pre-processing, different settings were applied to the EEG, such as ICA MARA artefact removal and high pass and low pass filters respectively set at 1 Hz and 45 Hz. The data were interpolated finding the locations of bad channels. The EEG data set can be displayed as electrode channel plots, allowing a quick overview of data quality. Then, the absolute power spectral density (PSD) was calculated and compared between PD and Healthy groups in each phase of the experiment for each of the delta, theta alpha and delta bands. A Mann-Whitney U-test with significance level (α = 0.05) was used to determine significance. False detection rate (FDR) correction was applied to each electrode.
3. Results
The following results are reported for the analysis of the PRE and POST phases of the experiment with the aim of distinguishing between the PD and healthy groups based on their biosignal responses. Our experiment was able to identify changes in many of the analyzed domains.
The protocol is a visual-motor simulation of being on a boat and part of the subjects experienced the feeling of seasickness. Just under half of PD subjects experienced actual discomfort with various symptoms (reported on questionnaires), and a smaller percentage of them reported a self-assessment of motion sickness in daily life.
3.1. Heart rate
The results of the heart rate analysis are shown in Table 4. An increase in beats per minute (bpm) was measured in PD subjects, although not statistically significant. No statistically significant differences were found between groups.
Table 4
| Beats per minute (BPM) | Parkinson's | Controls |
|---|---|---|
| PRE-phase | 88.76 ± 18.09 | 82.48 ± 14.63 |
| POST-phase | 94.19 ± 20.97 | 82.54 ± 14.35 |
Average and standard deviation for Heart Rates inside cohorts in PRE and POST phases.
3.2. EMG
The right TA muscle showed a number of statistically significant features in the POST phase, with a p value (p < 0.05) and the corresponding effect sizes for each variable shown in bold in Table 5. The right side could be considered the dominant leg in the prevalence of the group. Significant changes were found also in the left soleus (MN - mean, p = 0.003, cliff delta = -0.64) in the POST phase and the soleus right which had one significant feature (MD -median, p = 0.007, cliff delta = -0.586) in the PRE phase.
Table 5
| Feature | PRE p-value | POST p-value | cliff delta POST |
|---|---|---|---|
| Total Power (PT) | - | - | 0.425 |
| Maximum Power (Pmax) | - | 0.0394 | 0.464 |
| Maximum frequency (Fmax) | - | - | 0.141 |
| Median Frequency (FMD) | - | - | 0.153 |
| Mean frequency(FMN) | - | - | 0.191 |
| Frequency Kurtosis (Fkurt) | - | - | 0.073 |
| Frequency skewness (Fskew) | - | - | 0.080 |
| Average amplitude change (AAC) | - | 0.0256 | 0.502 |
| Average energy (AE) | - | - | 0.425 |
| Absolute value of the summation of the exponential root(ASM) | - | 0.0195 | 0.524 |
| Absolute value of the summation of the square root (ASS) | - | 0.0234 | 0.0509 |
| Coefficient of variation (CV) | - | - | -0.172 |
| Difference absolute mean value (DAMV) | - | 0.0256 | 0.502 |
| Difference absolute standard deviation value (DASDV) | - | 0.0362 | 0.471 |
| Difference variance value (DVARV) | - | 0.0362 | 0.471 |
| Enhanced mean absolute value (EMAV) | - | 0.0234 | 0.510 |
| Enhanced wavelength (EWL) | - | 0.0162 | 0.0.540 |
| New zero crossing (FZC) | - | - | 0.333 |
| Kurtosis (KURT) | - | 0.0162 | -0.540 |
| Integrated EMG (IEMG) | - | 0.0214 | 0.517 |
| Interquartile range (IQR) | - | 0.0394 | 0.463 |
| Log CV(LCV) | - | - | |
| Log Detector (LD) | - | 0.0428 | -0.02 |
| Log DAMV (LDAMV) | - | 0.0256 | 0.502 |
| Log DASDV (LDASDV) | - | 0.0362 | 0.417 |
| Log Teager Kaiser energy operator (LTKEO) | - | 0.0394 | 0.464 |
| Mean absolute deviation (MAD) | - | 0.0256 | 0.502 |
| Mean absolute value (MAV) | - | 0.0256 | 0.502 |
| Maximum fractal length (MFL) | - | 0.0362 | 0.471 |
| Mean (MN) | - | - | 0.210 |
| Median (MD) | - | - | -0.241 |
| Modified mean absolute value (MMAV) | - | 0.0234 | 0.510 |
| Modified mean absolute value 2 (MMAV2) | - | 0.0428 | 0.455 |
| Mean value of the square root (MSR) | - | 0.0195 | 0.524 |
| Root mean square (RMS) | - | - | |
| Standard deviation (SD) | - | - | 0.425 |
| Skewness (SKEW) | - | - | 0.333 |
| Single square integral (SSI) | - | - | 0.425 |
| Absolute value of temporal moment (TM) | - | - | 0.433 |
| Variance (VAR) | - | - | 0.425 |
| Variance of EMG (VARE) | - | - | 0.425 |
| Variance order (VO) | - | - | 0.425 |
| Waveform length (WL) | - | 0.0256 | 0.502 |
EMG features—bold shows features that were significantly different between PD and healthy groups in the right tibialis anterior muscle in the POST phase.
3.3. CoP
3.3.1. Sway profile
Figure 5 highlights the CoP evolution between PD and healthy subjects on two of its main characteristics area and axis length of the sway ellipse. Sway is greater in healthy than PD participants.
The only significant feature (p < 0.05) for CoP in the PRE phase is the Direction Entropy (Nats), while the statistically significant ones in the POST are listed below Figure 5 with p-values and cliff delta values listed.
3.4. EEG
The topological plots were computed for all frequency bands during the phases of the acquisition. Each of them displays the difference of power spectral density between PD and healthy cohorts, only for the statistically significant electrodes (p ≤ 0.05, represented by a green point in the figure). Theta and alpha bands presented several significant electrodes in different locations of the brain. In the theta band, significant electrodes are found mostly in the temporal lobe (T7, T8, C6, FT7), one in the frontal lobe (AF3) and one in the occipital (PO6). In the alpha band, significant electrodes are found mostly in the temporal lobe (FC5, T7, T8, FT8), one in the parietal lobe (P4) and one in the occipital (PO6). The p-values of each electrode of the theta and alpha bands are shown in Tables 6, 7, comparing the two cohorts in the PRE and POST phases. They highlight the differences in brain activity in the two phases and underline the significant difference in the POST phase of the experiment between the PD and Healthy groups.
Table 6
| Electrodes | PRE phase | POST phase |
|---|---|---|
| M1 | 0.124 | 0.0326 |
| T7 | 0.0740 | 0.0169 |
| T8 | 0.139 | 0.0215 |
| AF3 | 0.0947 | 0.0383 |
| C6 | 0.0784 | 0.0247 |
| PO6 | 0.0983 | 0.0409 |
| FT7 | 0.112 | 0.0476 |
Electrodes for theta band with corresponding p-values comparing the two cohorts in POST and PRE phases-bold shows they were significantly different in the POST phase of the experiment.
Table 7
| Electrodes | PRE phase | POST phase |
|---|---|---|
| FC5 | 0.0597 | 0.0247 |
| T7 | 0.122 | 0.0187 |
| T8 | 0.0703 | 0.0247 |
| P4 | 0.0596 | 0.0391 |
| PO6 | 0.276 | 0.0247 |
| FT8 | 0.0674 | 0.0165 |
Electrodes for alpha band with corresponding p-values comparing the two cohorts in POST and PRE phases- bold shows they were significantly different in the POST phase of the experiment.
4. Discussion
In previous work, we compared the identical PRE and POST phases in order to characterize different cohorts while they performed unassisted standing under conditions of sensory conflict after exposure to the complex postural control task of the moving platform and VR scene. EEG delta and theta power spectrum analysis and EMG activity in the soleus muscle proved to be strong discriminators between groups (Jacob et al.,
In this paper, we use a similar multi-factorial approach to characterize PD and healthy participants on the basis of their postural control response during the BioVRSea experiment. In particular, we compared the balance response after a visual stimulus only–PRE phase: VR visual sea motion simulation; to the balance response obtained after a complete immersive sensory experience–POST phase: VR visual and correspondent motion stimulation.
This study is part of an extensive work in which a larger population is monitored. We have collected data from 324 volunteers (females 183, males 141, general age 33 ± 14). The overall population is between the ages of 18 and 29. Recruiting a larger group of older adults may open a new way to highlight age-related postural control strategies using the same protocol.
4.1. Heart rate
Although heart rate per minute appears to increase in subjects with Parkinson's from PRE to POST phase, there is no significant result to underscore the difference in the two cohorts.
4.2. Muscle activation
The lower leg muscles are involved in postural and balance control strategies (Loram et al.,
We can hypothesize that in the group of subjects considered in this study (Parkinson's and healthy) has overall right-handed prevalence. This is supported in one of the latest reported studies conducted on the human hand, from which it was stated that the precise prevalence is the right hand in the world population. The prevalence of left-handedness is between 10–20% (Papadatou-Pastou et al., 2020). Thus, we can explain the significant results obtained for the right side through the dominance of the legs.
4.3. Center of pressure
A common way to evaluate PD is based on gait analysis with accelerometer and force sensors inside the shoe. Gait analysis has revealed higher frequency values for PD compared to healthy controls. However, until now Parkinson was not assessed by measuring the force in a standing position (Hsieh and Abbod,
Our results highlight a reduced sway in PD subjects during the task ant that may be related to a multitude of factors. Some diseases, including PD, interfere in the ability to maintain balance. PD patients have less coordination of agonist and antagonist muscles, making it challenging to maintaining stability. They also frequently suffer from limb and axial rigidity that may reduce mobility (Gandolfi et al.,
We found significant differences in the POST phase between the PD and healthy cohorts. This is significant as the POST phase is the stage after a motor stimulus and it can be a good index of pathology progression. The most discriminating feature was found to be the complexity index in the medio-lateral direction as seen in Figure 6. The effect size for this variable is also large (-0.793 - the negative value indicating that most of the higher values for this variable were in the PD group) which is classified as a large effect as per (Vargha and Delaney, 2000) which ranks a delta value greater than 0.42 as a large effect. The mediolateral CI in the PD group is higher than in the healthy cohort, which is contrary to a number of studies which show that the complexity of postural dynamics tends to decrease in disease and aging (Habtemariam et al.,
Figure 6

Mediolateral Complexity Index of healthy (blue) and Parkinson (orange) groups.
4.4. Neural response
The strong involvement of the cerebral cortex in postural control responses to perturbation is well-known, but the correlations with pathologies affecting mobility are still poorly understood (Jacobs and Horak,
Each band can be associated with a particular neural function and in our research, differences in alpha and theta response prove to be statistically significant when comparing between the two groups. The result for alpha band shows that the activity in the healthy subjects is greater than in the Parkinson subjects in both the PRE and POST phases. On the other hand, the theta band shows different behavior depending on the brain area and phase. Theta and alpha bands are involved in the regulation of the posture, in particular when a visual feedback is altered (Kahya et al.,
4.4.1. Theta waves
The theta rhythm is one of the slowest oscillations in the normal waking state, just above the delta rhythm that dominates slow wave sleep. Theta waves are involved in attention and memory processes, especially in memory retrieval episodes (Baars and Gage,
In our results, parietal activity of the theta band from PRE to POST phase shows a decrease of PSD in healthy subjects, confirming what was found in the work just mentioned. Instead for the frontal lobe, the theta band shows an increased PSD in the healthy group in both tasks. Theta brain rhythms are associated with cognitive and motor functions, and patients with PD would have irregular theta rhythms during lower-limb activations (Singh et al., 2020).
4.4.2. Alpha waves
Different studies have shown that the performance of a generic balance task results in simultaneous changes in the amplitude of alpha oscillations (Slobounov et al., 2013; Malisova et al., 2017). Alpha activation is associated with cognitive events and has been found to increase during intentional tasks such as mental arithmetic and working memory. Planning actions and their execution also generate alpha (Pfurtscheller, 2003). Therefore, alpha waves play a functional role in human cognition and that it does not represent only an ‘idling rhythm, as many scientists believed until recently. Maintaining balance is an active process and requires constant awareness of any external stimuli. The alpha band has shown interesting results in postural control studies where a decrease in alpha power was associated with an increased task difficulty during upright stance in young adults (Percio et al., 2007; Hülsdünker et al.,
Figure 7

Absolute PSD bands analysis in Parkinson's vs. healthy cohort for theta and alpha frequency band.
Both the bands (theta and alpha) can be considered parameters to discriminate PD subjects during a complex postural control task and confirm the activation of the frontal, parietal and visual lobe in healthy subjects, underlining the difficulty the PD group experiences when making postural adjustments.
5. Conclusions and limitations
Our paper confirms that a comprehensive multi-factorial approach (our unique BioVRSea paradigm) is useful in discriminating early PD subjects. In addition, the ability to look at multiple parameters at once introduces the ability to further analyse the correlation and timing of our set of specific features changes. In fact, there are no previous multi-metric experiments, such as BioVRSea, and there are no experiments with larger patient cohorts, therefore a larger study is foreseen to provide more definitive conclusion. Several neural and motor strategies difference have been highlighted, and are in line with known literature.
It is also worth providing some elements for consideration, specifically in the PD diagnosis. Firstly, the 'early stage classification' has no absolute value, it is a subjective evaluation and may suffers from a degree of variability. This influences the time of the diagnosis, which, in turn, may also present early stage cohorts with non-homogeneous symptoms and impairment. This BioVRsea paradigm seems very promising for comprehensive quantitative assessments and may pave the way for highlighting the most relevant features in the specific motor diseases analyzed. This may favor the introduction of less complex quantitative analyses, specifically for clinical operators, to replicate only the most poignant aspects of our paradigm. Our results suggest that a simpler experimental design including concurrent EEG, bilateral lower limbs EMG and CoP analyses with a balance challenge (not necessarily VR led) can discriminate early stage PD and has the potential to stratify further stages of the pathology.
Postural control alteration is one of the major risk factors to facilitate the occurrence of falls in the elderly. The frailty index, which is the expression of the health status of older individuals, is related to falls (Tornero-Quinones et al., 2020; Taguchi et al., 2022). Higher levels of frailty indicate the presence of chronic diseases (Vinik et al., 2017), such as diabetes, chronic pain, and polypharmacy. Taken separately, those clinical conditions have been established to change the postural control and increase falls (Efstathiou et al.,
BioVRSea is an innovative model to analyze postural control response that can be used for different clinical applications. It may also prove useful to classify other related diseases and conditions of movement disorders, such as Progressive Supranuclear Palsy (PSP) or Huntington's disease (Porciuncula et al., 2020) and many other disorders that involve changes in postural control such as loss of balance and slow movement, where a similar performance can be expected. Having new collaborations with different disease categories will add more value and amplify the research.
Statements
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
The studies involving human participants were reviewed and approved by Icelandic National Bioethics Committee (Application no: 17-183-S1). The patients/participants provided their written informed consent to participate in this study.
Author contributions
The idea for the paper was conceived by PG, HP, and VM. FP, LG, DJ, and SP processed the raw data from EEG, EMG, CoP, and HR and made the figures and tables. CG performed statistical analysis on the processed data. FP, LG, DJ, CG, PG, GD, and AF contributed to the text for the manuscript. All authors contributed to the article and approved the submitted version.
Acknowledgments
The authors would like to thank all those who participated in the experiment, particularly those of the Parkinson's group. We also thank RANNIS and their support for the project under grant number 239612-051, postural control signature: Advanced assessment and diagnostics. We also thank the Landspitali Innovation Fund.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1
AgliardiC.MeloniM.Rosa GueriniF.ZanzotteraM.BolognesiE.BaglioF.et al. (2021). Oligomeric α-syn and snare complex proteins in peripheral extracellular vesicles of neural origin are biomarkers for Parkinson's disease. Neurobiol. Dis. 148, 105185. 10.1016/j.nbd.2020.105185
2
AubonnetR.ShoykhetA.JacobD.LorenzoG. D.PetersenH.GargiuloP. (2022). Postural control paradigm (biovrsea): towards a neurophysiological signature. Physiol. Meas. 43, 115002. 10.1088/1361-6579/ac9c43
3
BaarsB. J.GageN. M. (2010). Cognition, Brain, and Consciousness: Introduction To Cognitive Neuroscience. Cambridge MA: Academic Press. 10.1016/B978-0-12-375070-9.00008-5
4
BarolloF.FriḥriksdóttirR.EdmundsK. J.KarlssonG. H.SvanssonH. ,°et al. (2020). Postural control adaptation and habituation during vibratory proprioceptive stimulation: an hd-eeg investigation of cortical recruitment and kinematics. IEEE Trans. Neural Syst. Rehabilitat. Eng. 28, 1381–1388. 10.1109/TNSRE.2020.2988585
5
BergD.GodauJ.WalterU. (2008). Transcranial sonography in movement disorders. Lancet Neurol. 7, 1044–1055. 10.1016/S1474-4422(08)70239-4
6
BosboomJ. L. W.StoffersD.StamC. J.DijkB. W. V.VerbuntJ.BerendseH. W.et al. (2006). Resting state oscillatory brain dynamics in parkinson's disease: an meg study. Clini. Neurophysiol. 117, 2521–2531. 10.1016/j.clinph.2006.06.720
7
BusaM. A.EmmerikR. E. A. (2016). Multiscale entropy: a tool for understanding the complexity of postural control. J. Sport Health Sci. 5, 44–51. 10.1016/j.jshs.2016.01.018
8
ChenL.LoW. L. A.MaoY. R.DingM. H.LinQ.LiH.et al. (2016). Effect of virtual reality on postural and balance control in patients with stroke: a systematic literature review. Biomed Res. Int. 2016, 7309272. 10.1155/2016/7309272
9
ChepishevaM. K. (2023). Spatial orientation, postural control and the vestibular system in healthy elderly and alzheimer's dementia. PeerJ. 11, e15040. 10.7717/peerj.15040
10
ChuC.ZhangZ.WangJ.LiuS.WangF.SunY.et al. (2021). Deep learning reveals personalized spatial spectral abnormalities of high delta and low alpha bands in eeg of patients with early parkinson's disease. J. Neural Eng. 18, 066036. 10.1088/1741-2552/ac40a0
11
CliffN. (1993). Dominance statistics: ordinal analyses to answer ordinal questions. Psychol. Bull. 114, 494. 10.1037/0033-2909.114.3.494
12
ComberL.SosnoffJ. J.GalvinR.CooteS. (2018). Postural control deficits in people with multiple sclerosis: a systematic review and meta-analysis. Gait Posture. 61, 445–452. 10.1016/j.gaitpost.2018.02.018
13
CurtisD. R.AndersenP. (2001). Sir john carew eccles, ac 27 january 1903-2 may 1997 Biogr. Mem. Fellows R. Soc. 47, 161–187. 10.1098/rsbm.2001.0010
14
EfstathiouM. A.GiannakiC. D.RoupaZ.HadjisavvasS.StefanakisM. (2022). Evidence of distorted proprioception and postural control in studies of experimentally induced pain: a critical review of the literature. Scand. J. Pain. 22, 445–456. 10.1515/sjpain-2021-0205
15
GandalM. J.EdgarJ. C.KlookK.SiegelS. J. (2012). Gamma synchrony: towards a translational biomarker for the treatment-resistant symptoms of schizophrenia. Neuropharmacology. 62, 1504–151810.1016/j.neuropharm.2011.02.007
16
GandolfiM.ValèN.FilippettiM.DimitrovaE. K.GeroinC.PicelliA.et al. (2018). “Postural control in individuals with parkinson's disease,” in Different Areas of Physiotherapy. London: IntechOpen. 10.5772/intechopen.81098
17
GoetzC. G.FahnS.Martinez-MartinP.PoeweW.SampaioC.StebbinsG. T.et al. (2007). Movement disorder society-sponsored revision of the unified Parkinson's disease rating scale (mds-updrs): process, format, and clinimetric testing plan. Mov. Dis. 22, 41–47. 10.1002/mds.21198
18
HabtemariamD.ZhouJ.IloputaifeI.LipsitzL.ManorB. (2017). The complexity of standing postural sway associates with future falls in community-dwelling older adults: the mobilize boston study. Sci. Rep. 7, 4. 10.1038/s41598-017-03422-4
19
HampelH.FrankR.BroichK.TeipelS. J.KatzR. G.HardyJ.et al. (2010). Biomarkers for alzheimer's disease: academic, industry and regulatory perspectives. Nat. Rev. Drug Discovery. 9, 560–574. 10.1038/nrd3115
20
HanC.-X.WangJ.YiG.-S.CheY.-Q. (2013). Investigation of eeg abnormalities in the early stage of parkinson's disease. Cogn. Neurodyn. 7, 351–35910.1007/s11571-013-9247-z
21
HsiehY.-L.AbbodM. F. (2021). Gait analyses of parkinson's disease patients using multiscale entropy. Electronics. 10, 2604. 10.3390/electronics10212604
22
HülsdünkerT.MierauA.NeebC.KleinöderH.StrüderH. K. (2015). Cortical processes associated with continuous balance control as revealed by eeg spectral power. Neurosci. Lett. 592, 1–5. 10.1016/j.neulet.2015.02.049
23
HülsdünkerT.MierauA.StrüderK. H. (2016). Higher balance task demands are associated with an increase in individual alpha peak frequency. Front. Hum. Neurosci. 9, 695. 10.3389/fnhum.2015.00695
24
JacobD.KristensenI. S. U.AubonnetR.RecentiM.DonisiL.RicciardiC.et al. (2022). Towards defining biomarkers to evaluate concussions using virtual reality and a moving platform (biovrsea). Sci. Rep. 12, 8996. 10.1038/s41598-022-12822-0
25
JacobsJ. V.HorakF. (2007). Cortical control of postural responses. J. Neural Transm. 114, 1339–1348. 10.1007/s00702-007-0657-0
26
KahyaM.LiaoK.GustafsonK. M.AkinwuntanA. E.ManorB.DevosH. (2022). Cortical correlates of increased postural task difficulty in young adults: A combined pupillometry and eeg study. Sensors. 22, 5594. 10.3390/s22155594
27
KheiriF.BraginA.AlmajanoJ.WindenE.et al. (2012). Non-linear classification of heart rate parameters as a biomarker for epileptogenesis. Epilepsy Res. 100, 59–66. 10.1016/j.eplepsyres.2012.01.008
28
LoramD. I.MaganarisC. N.LakieM. (2005). Human postural sway results from frequent, ballistic bias impulses by soleus and gastrocnemius. J. Physiol. 564, 295–311. 10.1113/jphysiol.2004.076307
29
LuF. M.YuanZ. (2015). Pet/spect molecular imaging in clinical neuroscience: recent advances in the investigation of cns diseases. Quant. Imaging Med. Surg. 5, 433. 10.3978/j.issn.2223-4292.2015.03.16
30
MaffoniM.GiardiniA.PierobonA.FerrazzoliD.FrazzittaG. (2017). Stigma experienced by parkinson's disease patients: a descriptive review of qualitative studies. Parkinson's Dis. 2017, 7203259. 10.1155/2017/7203259
31
MakiB. E.McIlroyW. E. (2007). Cognitive demands and cortical control of human balance-recovery reactions. J. Neural Transm. 114, 1279–1296. 10.1007/s00702-007-0764-y
32
MalisovaD. V.PetrenkoT. I.KondratenkoA. V.BazanovaO. M. (2017). “Alpha eeg and emg indices of optimal music-performing movement and postural control.” in XXIII (Moscow), 967–968.
33
María MaitínA.José García-TejedorA.Pablo Romero Mu nozJ. (2020). Machine learning approaches for detecting parkinson's disease from eeg analysis: a systematic review. Applied Sciences10, 8662. 10.3390/app10238662
34
MartignonC.PedrinollaA.RuzzanteF.GiuriatoG.Giuseppe LaginestraF.Bouça-MachadoR.et al. (2021). Guidelines on exercise testing and prescription for patients at different stages of parkinson's disease. Aging Clin. Exp. Res. 33, 221–246. 10.1007/s40520-020-01612-1
35
Moazami-GoudarziM.SarntheinJ.MichelsL.MoukhtievaR.JeanmonodD. (2008). Enhanced frontal low and high frequency power and synchronization in the resting eeg of parkinsonian patients. Neuroimage. 41, 985–997. 10.1016/j.neuroimage.2008.03.032
36
MochizukiG.BoeS.MarlinA.McIlRoyW. E. (2010). Perturbation-evoked cortical activity reflects both the context and consequence of postural instability. Neuroscience. 170, 599–609. 10.1016/j.neuroscience.2010.07.008
37
MolhemiF.MonjeziS.MehravarM.Shaterzadeh-YazdiM.-J.SalehiR.HesamS.et al. (2021). Effects of virtual reality vs conventional balance training on balance and falls in people with multiple sclerosis: a randomized controlled trial. Arch. Phys. Med. Rehabil. 102, 290–299. 10.1016/j.apmr.2020.09.395
38
NeufeldM. Y.InzelbergR.KorczynA. D. (1988). Eeg in demented and non-demented parkinsonian patients. Acta Neurol. Scand. 78, 1–5. 10.1111/j.1600-0404.1988.tb03609.x
39
NymanS. R.IngramW.SandersJ.ThomasP. W.ThomasS.VassalloM.et al. (2019). Randomised controlled trial of the effect of tai chi on postural balance of people with dementia. Clin. Interv. Aging. 14, 2017–2029. 10.2147/CIA.S228931
40
O'KeeffeC.TaboadaL.FeerickN.GallagharL.LynchT.ReilyR. (2019). Complexity based measures of postural stability provide novel evidence of functional decline in fragile x premutation carriers. J. Neuro. Engineering Rehabil. 16, 6. 10.1186/s12984-019-0560-6
41
Papadatou-PastouM.NtolkaE.SchmitzJ.MartinM.MunafòM. R.OcklenburgS.et al. (2020). Human handedness: a meta-analysis. Psychol. Bull. 146, 481. 10.1037/bul0000229
42
PapegaaijS.TaubeW.BaudryS.OttenE.HortobágyiT. (2014). Aging causes a reorganization of cortical and spinal control of posture. Front. Aging Neurosci. 6, 28. 10.3389/fnagi.2014.00028
43
ParreiraR. B.GreccoL. A. C.OliveiraC. S. (2017). Postural control in blind individuals: A systematic review. Gait &Posture. 57, 161–167. 10.1016/j.gaitpost.2017.06.008
44
PedroniA.BahreiniA.LangerN. (2019). Automagic: standardized preprocessing of big eeg data. Neuroimage. 200, 460–473. 10.1016/j.neuroimage.2019.06.046
45
PercioC. D.BrancucciA.BergamiF.MarzanoN.FioreA.CioloE. D.et al. (2007). Cortical alpha rhythms are correlated with body sway during quiet open-eyes standing in athletes: A high-resolution eeg study. Neuroimage. 36, 822–829. 10.1016/j.neuroimage.2007.02.054
46
PerezO. H.GreenR. E.MochizukiG. (2020). Spectral analysis of centre of pressure identifies altered balance control in individuals with moderate-severe traumatic brain injury. Disabil. Rehabil. 42, 519–527. 10.1080/09638288.2018.1501101
47
PfurtschellerG. (2003). Induced oscillations in the alpha band: functional meaning. Epilepsia. 44:2–8. 10.1111/j.0013-9580.2003.12001.x
48
PolitisM. (2014). Neuroimaging in parkinson disease: from research setting to clinical practice. Nat. Rev. Neurol. 10, 708–722. 10.1038/nrneurol.2014.205
49
PorciunculaF.WassermanP.MarderK. S.RaoA. K. (2020). Quantifying postural control in premanifest and manifest huntington disease using wearable sensors. Neurorehabil. Neural Repair. 34, 771–783. 10.1177/1545968320939560
50
PostumaB. R.BergD.SternM.PoeweW.OlanowC. W.OertelW.et al. (2015). Mds clinical diagnostic criteria for Parkinson's disease. Mov. Dis. 30, 1591–1601. 10.1002/mds.26424
51
RasmussenN. H.-H.DalJ.JensenM. H.KvistA. V.BerghJ. v. d.et al. (2023). Impaired postural control in diabetes a predictor of falls?Arch. Osteoporos. 18, 1–16. 10.1007/s11657-022-01188-5
52
RecentiM.RicciardiC.AubonnetR.PiconeI.JacobD.SvanssonH.et al. (2021). Toward predicting motion sickness using virtual reality and a moving platform assessing brain, muscles, and heart signals. Front. Bioeng. Biotechnol. 9, 635661. 10.3389/fbioe.2021.635661
53
SerizawaK.KameiS.MoritaA.HaraM.MizutaniT.YoshihashiH.et al. (2008). Comparison of quantitative eegs between parkinson disease and age-adjusted normal controls. J. Clini. Neurophysiol. 25, 361–366. 10.1097/WNP.0b013e31818f50de
54
SinanovićO.KapidzićA.KovacevićL.HudićJ.SmajlovićD. (2005). Eeg frequency and cognitive dysfunction in patients with Parkinson's disease. Med. Arh. 59, 286–287.
55
SinghA.ColeR. C.EspinozaA. I.BrownD.CavanaghJ. F.NarayananN. S. (2020). Frontal theta and beta oscillations during lower-limb movement in parkinson's disease. Clini. Neurophysiol. 131, 694–702. 10.1016/j.clinph.2019.12.399
56
SippA. R.GwinJ. T.MakeigS.FerrisD. P. (2013). Loss of balance during balance beam walking elicits a multifocal theta band electrocortical response. J. Neurophysiol. 110, 2050–2060. 10.1152/jn.00744.2012
57
SlobounovS. M.TeelE.NewellK. M. (2013). Modulation of cortical activity in response to visually induced postural perturbation: combined vr and eeg study. Neurosci. Lett. 547, 6–9. 10.1016/j.neulet.2013.05.001
58
TaguchiC. K.MenezesP. L.MeloA. C. S.SantanaL. S.Conceiç aoW. R.SouzaG. F.et al. (2022). “Frailty syndrome and risks for falling in the elderly community,” in CoDAS. SciELO Brasil. 10.1590/2317-1782/20212021025en
59
TolosaE.WenningG.PoeweW. (2006). The diagnosis of Parkinson's disease. Lancet Neurol. 5, 75–86. 10.1016/S1474-4422(05)70285-4
60
TorchianoM. (2020). effsize: Efficient Effect Size Computation.
61
Tornero-QuinonesI.Sáez-PadillaJ.DíazA. E.RoblesM. T. A.RoblesṠ. (2020). Functional ability, frailty and risk of falls in the elderly: relations with autonomy in daily living. Int. J. Environ. Res. Public Health17, 1006. 10.3390/ijerph17031006
62
VarghaA.DelaneyD. H. (2000). A critique and improvement of the cl common language effect size statistics of mcgraw and wong. J. Educ. Behav. Stat. 25, 101–132. 10.3102/10769986025002101
63
VinikA. I.CamachoP.ReddyS.ValenciaW. M.TrenceD.MatsumotoA. M.et al. (2017). Aging, diabetes, and falls. Endocrine Prac. 23, 1120–1142. 10.4158/EP171794.RA
64
WaningerS.BerkaC.Stevanovic KaricM.KorszenS.MozleyP. D.HenchcliffeC.et al. (2020). Neurophysiological biomarkers of parkinson's disease. J. Parkinson's Dis. 10, 471–480. 10.3233/JPD-191844
Summary
Keywords
postural control, early-stage Parkinson's disease, quantitative neurophysiology, BioVRSea, balance control
Citation
Jacob D, Guerrini L, Pescaglia F, Pierucci S, Gelormini C, Minutolo V, Fratini A, Di Lorenzo G, Petersen H and Gargiulo P (2023) Adaptation strategies and neurophysiological response in early-stage Parkinson's disease: BioVRSea approach. Front. Hum. Neurosci. 17:1197142. doi: 10.3389/fnhum.2023.1197142
Received
30 March 2023
Accepted
28 June 2023
Published
17 July 2023
Volume
17 - 2023
Edited by
Federico Giove, Centro Fermi – Museo storico della fisica e Centro studi e ricerche Enrico Fermi, Italy
Reviewed by
Davide Balos Cappon, Harvard Medical School, United States; Grigorios Nasios, University of Ioannina, Greece
Updates

Check for updates
Copyright
© 2023 Jacob, Guerrini, Pescaglia, Pierucci, Gelormini, Minutolo, Fratini, Di Lorenzo, Petersen and Gargiulo.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Paolo Gargiulo paolo@ru.is
†These authors have contributed equally to this work and share first authorship
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.