- 1Institute of Biomedical and Neural Engineering, Reykjavik University, Reykjavik, Iceland
- 2Department of Engineering, University of Campania L. Vanvitelli, Aversa, Italy
- 3Department of Electrical, Electronic and Information Engineering, University of Bologna, Cesena, Italy
- 4Department of Pharmacy and Biotechnology, University of Bologna, Bologna, Italy
- 5Department of Civil Engineering and Computer Science Engineering, Tor Vergata University of Rome, Rome, Italy
- 6Engineering for Health Research Centre, Aston University, Birmingham, United Kingdom
- 7Laboratory of Psychophysiology and Cognitive Neuroscience, Department of Systems Medicine, Tor Vergata University of Rome, Rome, Italy
- 8IRCCS Fondazione Santa Lucia, Rome, Italy
- 9Department of Anatomy, University of Iceland, Reykjavik, Iceland
- 10Department of Science, Landspitali University Hospital, Reykjavik, Iceland
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., 2017; 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., 2007). 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., 2007). Middle-stage balance and coordination are affected: a moderate-to-severe disability that affects daily life (Goetz et al., 2007). Later-stage subjects have difficulty standing and walking even with aids. Patients in this stage have severe disability (Goetz et al., 2007). Figure 1 reports the general classification of PD stages according to its main symptoms.
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. Quantitative and qualitative aspects in Parkinson's Disease diagnosis (Maffoni et al., 2017).
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.
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., 2021).
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, 2015). Single positron emission tomography (SPECT) is also used as method to confirm a Parkinson's diagnosis by highlighting cerebral blood flow and dopamine transporters in the brain (Lu and Yuan, 2015). In patients with Parkinson's disease, a distinct intensity pattern can be noted in the brain region that is deeply affected by degeneration, namely the basal ganglia that controls movement. Magnetic resonance imaging (MRI) is also used to diagnose PD in early onset subjects. MRIs can show small changes and damage in the brain tissue that can indicate PD. Often, these markers are present even before symptoms of PD begin. Transcranial sonography (TCS) has been established as a valuable supplementary tool in the diagnosis of PD. Alterations in the area of the hypoechogenic mesencepahlic brainstem can be visualized in about 90% of PD patients, which is measured planimetrically to determine the magnitude of the change (Berg et al., 2008). Increased iron levels contribute to this sonographic abnormality and indicate that iron can be responsible for the change in the echo signal. Iron accumulation can be a very early indicator in the pathogenesis of PD. Electroencephalography (EEG) can detect damage in the central nervous system and alterations in neurophysiological activity associated with PD. In recent studies, quantitative analysis of EEG data identified significant differences in PD patients versus healthy subjects. In particular, the anterior cingulate and temporal lobe are areas with an established pathology in PD. Changes in cortico-cortical and cortico-thalamic coupling were observed as excessive EEG beta coherence in PD patients (Waninger et al., 2020).
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., 2010; Gandal et al., 2012; Kheiri et al., 2012). The combination of new analysis methods and EEG signal processing can contribute to the detection of early-stage PD. EEG reveals more important information underlying brain dysfunctions, which would be lost if analysis were restricted to traditional methods. Nowadays, many novel methods are suggested for EEG signal processing.
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., 2013). Most EEG electrodes showed an increase in θ-band relative power for PD patients, while several other electrodes decreased, such as in the frontal and occipital cortex (Fp1, Fp2, F7, F3, Fz, Oz). Moreover, an increase in δ-band relative powers were reported, and a decrease in α-band and β-band relative powers for PD patients compared with healthy patients. Other studies present higher spectral power in the low frequency domain of EEG, compared with controls. Also in these cases, subjects were in the resting awake condition with the eyes closed (Moazami-Goudarzi et al., 2008; Serizawa et al., 2008).
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., 2020). PD usually interferes in this regulatory system, as can be clearly demonstrated by most motor symptoms, but to date, no study has yet been conducted on the analysis of postural kinematics in movement disorders. Our aim is to (i) investigate the postural strategy adopted in PD individuals and in healthy subjects; (ii) describe adaptation and how the brain adapts to the induced movement of a platform and visual stimuli using virtual reality (VR). However, postural control and adaptation have been extensively studied in healthy and blind subjects. In a recent study (Barollo et al., 2020), postural kinematics from HD-EEG have been measured during a postural perturbation applied to calf muscles. The main changes in cortical activity were found in Absolute Spectral Power (ASP) over four frequency bands. For postural adaptation, increases in the θ band in the frontal-central region for closed-eyes trials, and in the θ and β bands in the parietal region for open-eyes trials were reported. In habituation of the stance, no significant variations in ASP were observed during closed-eyes trials, whereas an increase in the θ, α, and β bands were observed with open eyes (Curtis et al., 2001). Furthermore, open-eyed trials generally yielded a greater number of significant differences across all bands during both adaptation and habituation, suggesting that cortical activity during postural perturbation may be regulated with visual feedback. This clearly shows a correspondence in cortical activity and postural kinematics during postural perturbation, and could also be developed for pathological postural control.
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., 2015). Moreover, the increase of the ASP in the central region is demonstrated during high-demand postural correction, such as balance maintenance without allowing corrective foot placement (Barollo et al., 2020). In accordance, the increase in θ activity in the frontal-central regions implies the processing of postural stability during balance control. Thus, ASP differences in the θ band signify the planning of corrective steps and the analysis of the consequences of the subject falling. Instead, the significant differences in the α band reflect an inhibition of error detection within the cingulate cortex due to habituation. Other studies have been carried out with the blind (including both congenital blindness and acquired blindness). Congenitally blind subjects had poorer postural control (anterior-posterior and medio-lateral body swing) compared to sighted subjects. They use a more efficient mechanism for maintaining balance control through joint stiffness. These findings demonstrate that motor coordination, localization, or perception of body segments and movements in visually impaired individuals can be compensated by enhancing the proprioceptive and vestibular systems. Blindness leads to impaired postural balance and imbalance in static and dynamic tasks (Parreira et al., 2017). An overview of EEG studies in postural control is provided in Figure 3.
Figure 3. ASP bands analysis in different studies (Neufeld et al., 1988; Sinanović et al., 2005; Bosboom et al., 2006; Han et al., 2013; Chu et al., 2021).
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., 2022). The advantage of using BioVRSea is that we are measuring quantitative signals associated with postural control in a challenging environment. The experiment is a prototype and the purpose of the current research with BioVRSea is to gather as much data as possible with many simultaneous measurements in order to extrapolate the most relevant features which could then be used in a clinical setting, with a lower-profile machine that could be accessed easily by those with mobility problems. Current diagnosis of PD relies primarily on the presence of motor symptoms in the patient (such as MDS diagnosis criteria (Postuma et al., 2015) and generally lacks any quantitative measurement such as we perform in the BioVRSea experiment.
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, 1993).
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, 2016). Multi-scale entropy measurements include features such as complexity index (CI), which indicate the complexity of the CoP signal as calculated using multi-scale entropy methods. We extract several parameters from the stabilogram for evaluating the postural control response of the subject during the experiment. The list of features extracted from the CoP is outlined in Table 3. Statistical analysis normality was checked through the Shapiro-Wilk test and visual inspection of the variables, and comparisons were made between the healthy and PD groups for all features in the PRE and POST phases using the t-test with Welch's correction and the Mann-Whitney U-test with significance level p < 0.05. Effect sizes were calculated using Cliff's delta. Sway profiles were also outlined using 95% confidence ellipses in PRE and POST, as seen in Figure 5.
Table 3. 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.
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.
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. 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 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 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., 2022).
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., 2005). Different muscles are involved at different times in two cohorts under the same experiment and this may be pathology dependent. The most significant activity was found in the tibialis anterior on the right leg, with some significant activity in the right soleus and the left tibialis anterior. However, only features obtained from the right tibialis anterior (TA) statistically differentiated the two groups. This may be related to the fact that tibialis anterior is the primary dorsiflexor of the foot, and is critical in gait to lift the foot during the swing phase.
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, 2021).
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., 2018). All of those linked to a reduced mobility confidence may account for the reduced sway.
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., 2017; O'Keeffe et al., 2019). It is unclear why the complexity index is higher in the pathological group for our study.
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, 2007; Maki and McIlroy, 2007; Papegaaij et al., 2014). Understanding the network of cortical structures involved in a disease such as Parkinson's and how sensory information are processed can be an important step in diagnostics.
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., 2022). Figure 6 shows the topological plots of absolute PSD for theta and alpha frequency bands. Each row represents the PSD difference of task (PRE and POST). The PSD differences between Parkinson and healthy subjects are compared for significant electrodes with Benjamin-Hochberd FDR procedure.
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, 2010). Although the alpha band has been shown to be strongly correlated within postural task conditions, it currently remains less known whether the theta band shows an association with increased postural task difficulty (Kahya et al., 2022). One of the most recent publications on neurophysiology in healthy subjects showed that the theta band has important electrodes located mainly in the parietal scalp, associated with a slight decrease in PSD (Aubonnet et al., 2022). The parietal lobe is activated to plan and process the orientation of the body and sensory information, demonstrating the remarkable role of the theta band in postural strategies.
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., 2016; Kahya et al., 2022). As the PD group has lower overall alpha power (as seen in Figure 7), this may indicate higher demand to cope with the balance task compared to the healthy group. In another publication on the neurophysiology of healthy subjects using the BioVRSea experiment (Aubonnet et al., 2022), the alpha band is important for balance control across the whole scalp. Our study is consistent with these results.
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., 2022; Rasmussen et al., 2023). We speculate that BioVRSea protocol may be added to the routine assessment of individuals with a high frailty index to identify early postural control deficit and address early intervention to prevent falls. Furthermore, individuals with other neurodegenerative disorders (i.e., multiple sclerosis - Comber et al., 2018; Molhemi et al., 2021 - or dementia - Nyman et al., 2019; Chepisheva, 2023), or neurological conditions (i.e., stroke - Chen et al., 2016 - or brain injury - Perez et al., 2020) may find advantage in the use BioVRSea protocol to explore the postural control alterations.
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.
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
Agliardi, C., Meloni, M., Rosa Guerini, F., Zanzottera, M., Bolognesi, E., Baglio, F., 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. doi: 10.1016/j.nbd.2020.105185
Aubonnet, R., Shoykhet, A., Jacob, D., Lorenzo, G. D., Petersen, H., and Gargiulo, P. (2022). Postural control paradigm (biovrsea): towards a neurophysiological signature. Physiol. Meas. 43, 115002. doi: 10.1088/1361-6579/ac9c43
Baars, B. J., and Gage, N. M. (2010). Cognition, Brain, and Consciousness: Introduction To Cognitive Neuroscience. Cambridge MA: Academic Press. doi: 10.1016/B978-0-12-375070-9.00008-5
Barollo, F., Friḥriksdóttir, R., Edmunds, K. J., Karlsson, G. H., Svansson, H. ,°, 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. doi: 10.1109/TNSRE.2020.2988585
Berg, D., Godau, J., and Walter, U. (2008). Transcranial sonography in movement disorders. Lancet Neurol. 7, 1044–1055. doi: 10.1016/S1474-4422(08)70239-4
Bosboom, J. L. W., Stoffers, D., Stam, C. J., Dijk, B. W. V., Verbunt, J., Berendse, H. W., et al. (2006). Resting state oscillatory brain dynamics in parkinson's disease: an meg study. Clini. Neurophysiol. 117, 2521–2531. doi: 10.1016/j.clinph.2006.06.720
Busa, M. A., and Emmerik, R. E. A. (2016). Multiscale entropy: a tool for understanding the complexity of postural control. J. Sport Health Sci. 5, 44–51. doi: 10.1016/j.jshs.2016.01.018
Chen, L., Lo, W. L. A., Mao, Y. R., Ding, M. H., Lin, Q., Li, H., 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. doi: 10.1155/2016/7309272
Chepisheva, M. K. (2023). Spatial orientation, postural control and the vestibular system in healthy elderly and alzheimer's dementia. PeerJ. 11, e15040. doi: 10.7717/peerj.15040
Chu, C., Zhang, Z., Wang, J., Liu, S., Wang, F., Sun, Y., 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. doi: 10.1088/1741-2552/ac40a0
Cliff, N. (1993). Dominance statistics: ordinal analyses to answer ordinal questions. Psychol. Bull. 114, 494. doi: 10.1037/0033-2909.114.3.494
Comber, L., Sosnoff, J. J., Galvin, R., and Coote, S. (2018). Postural control deficits in people with multiple sclerosis: a systematic review and meta-analysis. Gait Posture. 61, 445–452. doi: 10.1016/j.gaitpost.2018.02.018
Curtis, D. R., and Andersen, P. (2001). Sir john carew eccles, ac 27 january 1903-2 may 1997 Biogr. Mem. Fellows R. Soc. 47, 161–187. doi: 10.1098/rsbm.2001.0010
Efstathiou, M. A., Giannaki, C. D., Roupa, Z., Hadjisavvas, S., and Stefanakis, M. (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. doi: 10.1515/sjpain-2021-0205
Gandal, M. J., Edgar, J. C., Klook, K., and Siegel, S. J. (2012). Gamma synchrony: towards a translational biomarker for the treatment-resistant symptoms of schizophrenia. Neuropharmacology. 62, 1504–1518 doi: 10.1016/j.neuropharm.2011.02.007
Gandolfi, M., Valè, N., Filippetti, M., Dimitrova, E. K., Geroin, C., Picelli, A., et al. (2018). “Postural control in individuals with parkinson's disease,” in Different Areas of Physiotherapy. London: IntechOpen. doi: 10.5772/intechopen.81098
Goetz, C. G., Fahn, S., Martinez-Martin, P., Poewe, W., Sampaio, C., Stebbins, G. 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. doi: 10.1002/mds.21198
Habtemariam, D., Zhou, J., Iloputaife, I., Lipsitz, L., and Manor, B. (2017). The complexity of standing postural sway associates with future falls in community-dwelling older adults: the mobilize boston study. Sci. Rep. 7, 4. doi: 10.1038/s41598-017-03422-4
Hampel, H., Frank, R., Broich, K., Teipel, S. J., Katz, R. G., Hardy, J., et al. (2010). Biomarkers for alzheimer's disease: academic, industry and regulatory perspectives. Nat. Rev. Drug Discovery. 9, 560–574. doi: 10.1038/nrd3115
Han, C.-X., Wang, J., Yi, G.-S., and Che, Y.-Q. (2013). Investigation of eeg abnormalities in the early stage of parkinson's disease. Cogn. Neurodyn. 7, 351–359 doi: 10.1007/s11571-013-9247-z
Hsieh, Y.-L., and Abbod, M. F. (2021). Gait analyses of parkinson's disease patients using multiscale entropy. Electronics. 10, 2604. doi: 10.3390/electronics10212604
Hülsdünker, T., Mierau, A., Neeb, C., Kleinöder, H., and Strüder, H. K. (2015). Cortical processes associated with continuous balance control as revealed by eeg spectral power. Neurosci. Lett. 592, 1–5. doi: 10.1016/j.neulet.2015.02.049
Hülsdünker, T., Mierau, A., and Strüder, K. H. (2016). Higher balance task demands are associated with an increase in individual alpha peak frequency. Front. Hum. Neurosci. 9, 695. doi: 10.3389/fnhum.2015.00695
Jacob, D., Kristensen, I. S. U., Aubonnet, R., Recenti, M., Donisi, L., Ricciardi, C., et al. (2022). Towards defining biomarkers to evaluate concussions using virtual reality and a moving platform (biovrsea). Sci. Rep. 12, 8996. doi: 10.1038/s41598-022-12822-0
Jacobs, J. V., and Horak, F. (2007). Cortical control of postural responses. J. Neural Transm. 114, 1339–1348. doi: 10.1007/s00702-007-0657-0
Kahya, M., Liao, K., Gustafson, K. M., Akinwuntan, A. E., Manor, B., and Devos, H. (2022). Cortical correlates of increased postural task difficulty in young adults: A combined pupillometry and eeg study. Sensors. 22, 5594. doi: 10.3390/s22155594
Kheiri, F., Bragin, A., Almajano, J., Winden, E., et al. (2012). Non-linear classification of heart rate parameters as a biomarker for epileptogenesis. Epilepsy Res. 100, 59–66. doi: 10.1016/j.eplepsyres.2012.01.008
Loram, D. I., Maganaris, C. N., and Lakie, M. (2005). Human postural sway results from frequent, ballistic bias impulses by soleus and gastrocnemius. J. Physiol. 564, 295–311. doi: 10.1113/jphysiol.2004.076307
Lu, F. M., and Yuan, Z. (2015). Pet/spect molecular imaging in clinical neuroscience: recent advances in the investigation of cns diseases. Quant. Imaging Med. Surg. 5, 433. doi: 10.3978/j.issn.2223-4292.2015.03.16
Maffoni, M., Giardini, A., Pierobon, A., Ferrazzoli, D., and Frazzitta, G. (2017). Stigma experienced by parkinson's disease patients: a descriptive review of qualitative studies. Parkinson's Dis. 2017, 7203259. doi: 10.1155/2017/7203259
Maki, B. E., and McIlroy, W. E. (2007). Cognitive demands and cortical control of human balance-recovery reactions. J. Neural Transm. 114, 1279–1296. doi: 10.1007/s00702-007-0764-y
Malisova, D. V., Petrenko, T. I., Kondratenko, A. V., and Bazanova, O. M. (2017). “Alpha eeg and emg indices of optimal music-performing movement and postural control.” in XXIII (Moscow), 967–968.
María Maitín, A., José García-Tejedor, A., and Pablo Romero Mu noz, J. (2020). Machine learning approaches for detecting parkinson's disease from eeg analysis: a systematic review. Applied Sciences 10, 8662. doi: 10.3390/app10238662
Martignon, C., Pedrinolla, A., Ruzzante, F., Giuriato, G., Giuseppe Laginestra, F., Bouça-Machado, R., 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. doi: 10.1007/s40520-020-01612-1
Moazami-Goudarzi, M., Sarnthein, J., Michels, L., Moukhtieva, R., and Jeanmonod, D. (2008). Enhanced frontal low and high frequency power and synchronization in the resting eeg of parkinsonian patients. Neuroimage. 41, 985–997. doi: 10.1016/j.neuroimage.2008.03.032
Mochizuki, G., Boe, S., Marlin, A., and McIlRoy, W. E. (2010). Perturbation-evoked cortical activity reflects both the context and consequence of postural instability. Neuroscience. 170, 599–609. doi: 10.1016/j.neuroscience.2010.07.008
Molhemi, F., Monjezi, S., Mehravar, M., Shaterzadeh-Yazdi, M.-J., Salehi, R., Hesam, S., 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. doi: 10.1016/j.apmr.2020.09.395
Neufeld, M. Y., Inzelberg, R., and Korczyn, A. D. (1988). Eeg in demented and non-demented parkinsonian patients. Acta Neurol. Scand. 78, 1–5. doi: 10.1111/j.1600-0404.1988.tb03609.x
Nyman, S. R., Ingram, W., Sanders, J., Thomas, P. W., Thomas, S., Vassallo, M., 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. doi: 10.2147/CIA.S228931
O'Keeffe, C., Taboada, L., Feerick, N., Gallaghar, L., Lynch, T., and Reily, R. (2019). Complexity based measures of postural stability provide novel evidence of functional decline in fragile x premutation carriers. J. Neuro. Engineering Rehabil. 16, 6. doi: 10.1186/s12984-019-0560-6
Papadatou-Pastou, M., Ntolka, E., Schmitz, J., Martin, M., Munafò, M. R., Ocklenburg, S., et al. (2020). Human handedness: a meta-analysis. Psychol. Bull. 146, 481. doi: 10.1037/bul0000229
Papegaaij, S., Taube, W., Baudry, S., Otten, E., and Hortobágyi, T. (2014). Aging causes a reorganization of cortical and spinal control of posture. Front. Aging Neurosci. 6, 28. doi: 10.3389/fnagi.2014.00028
Parreira, R. B., Grecco, L. A. C., and Oliveira, C. S. (2017). Postural control in blind individuals: A systematic review. Gait &Posture. 57, 161–167. doi: 10.1016/j.gaitpost.2017.06.008
Pedroni, A., Bahreini, A., and Langer, N. (2019). Automagic: standardized preprocessing of big eeg data. Neuroimage. 200, 460–473. doi: 10.1016/j.neuroimage.2019.06.046
Percio, C. D., Brancucci, A., Bergami, F., Marzano, N., Fiore, A., Ciolo, E. 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. doi: 10.1016/j.neuroimage.2007.02.054
Perez, O. H., Green, R. E., and Mochizuki, G. (2020). Spectral analysis of centre of pressure identifies altered balance control in individuals with moderate-severe traumatic brain injury. Disabil. Rehabil. 42, 519–527. doi: 10.1080/09638288.2018.1501101
Pfurtscheller, G. (2003). Induced oscillations in the alpha band: functional meaning. Epilepsia. 44:2–8. doi: 10.1111/j.0013-9580.2003.12001.x
Politis, M. (2014). Neuroimaging in parkinson disease: from research setting to clinical practice. Nat. Rev. Neurol. 10, 708–722. doi: 10.1038/nrneurol.2014.205
Porciuncula, F., Wasserman, P., Marder, K. S., and Rao, A. K. (2020). Quantifying postural control in premanifest and manifest huntington disease using wearable sensors. Neurorehabil. Neural Repair. 34, 771–783. doi: 10.1177/1545968320939560
Postuma, B. R., Berg, D., Stern, M., Poewe, W., Olanow, C. W., Oertel, W., et al. (2015). Mds clinical diagnostic criteria for Parkinson's disease. Mov. Dis. 30, 1591–1601. doi: 10.1002/mds.26424
Rasmussen, N. H.-H., Dal, J., Jensen, M. H., Kvist, A. V., Bergh, J. v. d., et al. (2023). Impaired postural control in diabetes a predictor of falls? Arch. Osteoporos. 18, 1–16. doi: 10.1007/s11657-022-01188-5
Recenti, M., Ricciardi, C., Aubonnet, R., Picone, I., Jacob, D., Svansson, H., 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. doi: 10.3389/fbioe.2021.635661
Serizawa, K., Kamei, S., Morita, A., Hara, M., Mizutani, T., Yoshihashi, H., et al. (2008). Comparison of quantitative eegs between parkinson disease and age-adjusted normal controls. J. Clini. Neurophysiol. 25, 361–366. doi: 10.1097/WNP.0b013e31818f50de
Sinanović, O., Kapidzić, A., Kovacević, L., Hudić, J., and Smajlović, D. (2005). Eeg frequency and cognitive dysfunction in patients with Parkinson's disease. Med. Arh. 59, 286–287.
Singh, A., Cole, R. C., Espinoza, A. I., Brown, D., Cavanagh, J. F., and Narayanan, N. S. (2020). Frontal theta and beta oscillations during lower-limb movement in parkinson's disease. Clini. Neurophysiol. 131, 694–702. doi: 10.1016/j.clinph.2019.12.399
Sipp, A. R., Gwin, J. T., Makeig, S., and Ferris, D. P. (2013). Loss of balance during balance beam walking elicits a multifocal theta band electrocortical response. J. Neurophysiol. 110, 2050–2060. doi: 10.1152/jn.00744.2012
Slobounov, S. M., Teel, E., and Newell, K. M. (2013). Modulation of cortical activity in response to visually induced postural perturbation: combined vr and eeg study. Neurosci. Lett. 547, 6–9. doi: 10.1016/j.neulet.2013.05.001
Taguchi, C. K., Menezes, P. L., Melo, A. C. S., Santana, L. S., Conceiç ao, W. R., Souza, G. F., et al. (2022). “Frailty syndrome and risks for falling in the elderly community,” in CoDAS. SciELO Brasil. doi: 10.1590/2317-1782/20212021025en
Tolosa, E., Wenning, G., and Poewe, W. (2006). The diagnosis of Parkinson's disease. Lancet Neurol. 5, 75–86. doi: 10.1016/S1474-4422(05)70285-4
Tornero-Quinones, I., Sáez-Padilla, J., Díaz, A. E., Robles, M. T. A., and Robles, Ṡ. (2020). Functional ability, frailty and risk of falls in the elderly: relations with autonomy in daily living. Int. J. Environ. Res. Public Health 17, 1006. doi: 10.3390/ijerph17031006
Vargha, A., and Delaney, D. 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. doi: 10.3102/10769986025002101
Vinik, A. I., Camacho, P., Reddy, S., Valencia, W. M., Trence, D., Matsumoto, A. M., et al. (2017). Aging, diabetes, and falls. Endocrine Prac. 23, 1120–1142. doi: 10.4158/EP171794.RA
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.
Edited by:
Federico Giove, Centro Fermi – Museo storico della fisica e Centro studi e ricerche Enrico Fermi, ItalyReviewed by:
Davide Balos Cappon, Harvard Medical School, United StatesGrigorios Nasios, University of Ioannina, Greece
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