- 1Department of Neurology, University Hospital Würzburg, Würzburg, Germany
- 2Academic Unit of Neurology, Trinity College Dublin, The University of Dublin, Dublin, Ireland
- 3Department of Neurology, University Medicine of Greifswald, Greifswald, Germany
- 4Paris Brain Institute, Paris, France
- 5Functional Imaging Unit, Center for Diagnostic Radiology, University Medicine Greifswald, Greifswald, Germany
- 6Department of Neurology, University Medical Center of the Johannes Gutenberg University Mainz, Mainz, Germany
The execution of voluntary movements is primarily governed by the cerebral hemisphere contralateral to the moving limb. Previous research indicates that the ipsilateral motor network, comprising the primary motor cortex (M1), supplementary motor area (SMA), and premotor cortex (PM), plays a crucial role in the planning and execution of limb movements. However, the precise functions of this network and its interplay in different task contexts have yet to be fully understood. Twenty healthy right-handed participants (10 females, mean age 26.1 ± 4.6 years) underwent functional MRI scans while performing biceps brachii representations such as bilateral, unilateral flexion, and bilateral flexion-extension. Ipsilateral motor evoked potentials (iMEPs) were obtained from the identical set of participants in a prior study using transcranial magnetic stimulation (TMS) targeting M1 while employing the same motor tasks. The voxel time series was extracted based on the region of interest (M1, SMA, ventral PM and dorsal PM). Directed functinal connectivity was derived from the extracted time series using time-resolved partial directed coherence. We found increased connectivity from left-PMv to both sides M1, as well as right-PMv to both sides SMA, in unilateral flexion compared to bilateral flexion. Connectivity from left M1 to left-PMv, and left-SMA to right-PMd, also increased in both unilateral flexion and bilateral flexion-extension compared to bilateral flexion. However, connectivity between PMv and right-M1 to left-PMd decreased during bilateral flexion-extension compared to unilateral flexion. Additionally, during bilateral flexion-extension, the connectivity from right-M1 to right-SMA had a negative relationship with the area ratio of iMEP in the dominant side. Our results provide corroborating evidence for prior research suggesting that the ipsilateral motor network is implicated in the voluntary movements and underscores its involvement in cognitive processes such as movement planning and coordination. Moreover, ipsilateral connectivity from M1 to SMA on the dominant side can modulate the degree of ipsilateral M1 activation during bilateral antagonistic contraction.
1 Introduction
A long-standing hypothesis that the voluntary limb movement is dominated by the contralateral hemisphere of the moving limb. It finds support in study demonstrating that electrical stimulation of the motor cortex leads to limb movement on the side opposite to the site of stimulation (Fleming, 1938). Furthermore, unilateral damage to the motor cortex often results in motor deficits on the contralateral side of the body (Duncan et al., 2018). Conversely, during unimanual limb movements, several studies have documented activation of motor areas in the same side or ipsilateral hemisphere, including the primary motor (M1) and premotor areas (PM) as opposed to the observed contralateral neural activity (Bundy and Leuthardt, 2019).
Functional imaging investigations have demonstrated the activation of the motor cortex in the same hemisphere as the moving limb. In contrast to the activation observed in the contralateral M1 during unilateral limb movements, the ipsilateral M1 is also active but with a smaller area of activation (Grafton et al., 1992; Kim et al., 1993). The strength of the ipsilateral motor activity has been found to associate with task complexity and the demands for accuracy (Rao et al., 1993; Buetefisch et al., 2014). In the context of unilateral isometric contraction from the dominant hand, there was documented evidence of a discernible bilateral network (Mayhew et al., 2017). Despite the evidence demonstrating the involvement of ipsilateral motor cortex during limb movement, only a limited number of studies have dived deeper into the brain network within the motor cortex. Most commonly, dynamic causal modelling (DCM), which enables the estimation of brain connectivity within a hypothesis-driven framework, is frequently employed to examine the relationships between regions of interest (ROIs). Grefkes and colleagues observed reduced connectivity towards ipsilateral motor areas during unimanual movements and increased connectivity towards the contralateral M1, with highlighted brain region in supplementary area (SMA) (Grefkes et al., 2008). In another study (Volz et al., 2015), the inhibitory effects on the ipsilateral M1 were found strongly influenced by PM and contralateral M1 homologue. However, connectivity analysis based on data-driven approach, which offers valuable alternative by leveraging observed data to infer the connectivity pattern and dynamics of the brain network, is sparsely applied to neuroimaging modalities to investigate the human motor network.
Furthermore, the recruitment of ipsilateral motor pathways can be examined neurophysiologically with transcranial magnetic stimulation (TMS) (Armand and Kuypers, 1980; Wassermann et al., 1994; Ziemann et al., 1999). In contrast to contralateral motor evoked potential (cMEP), ipsilateral motor evoked potential (iMEP) exhibits distinctive characteristics, including a delayed onset, higher activation threshold, and reduced amplitude. These features suggest a weaker and potentially indirect pathway to spinal alpha motor neurons. Notably, iMEPs are more readily elicited in proximal muscles compared to distal muscles (Wassermann et al., 1994; Bawa et al., 2004) and are typically observed when the target muscle is pre-activated (Chen et al., 2003; Bawa et al., 2004). These characteristics underscore the potential significance of iMEPs in the context of bimanual or postural motor interactions. The excitability of iMEPs depending on the task context was previously suggested by Tazoe and colleagues (Tazoe and Perez, 2014). In their study, they identified the largest iMEPs in bilateral contraction of heteronymous muscles, which involved the contraction of the non-dominant arm’s biceps brachii (elbow flexion) concurrently with the contraction of the dominant arm’s triceps brachii (elbow extension). Additionally, iMEP obtained during bilateral elbow flexion from the non-dominant side were smaller compared to those observed during unilateral elbow flexion (Perez et al., 2013; Seusing et al., 2023). However, in a similar paradigm to our recent work (Seusing et al., 2023), we found that the size of the iMEP was not influenced by the different upper limb movements. Hence, the functional significance of iMEP and its precise underlying neural mechanism continue to be unclear and subject to ongoing debate.
Understanding the role of ipsilateral motor activity is crucial for improving our knowledge of how the brain facilitates voluntary movements in healthy individuals and for understanding the potential role played by the contralesional hemisphere during the recovery of ipsilateral motor function after brain injuries such as stroke or traumatic brain injury (Bundy and Leuthardt, 2019). By integrating iMEP assessments and neuroimaging techniques, we can uncover how the excitability of ipsilateral motor activity, as assessed by iMEP, relates to brain connectivity patterns detected through fMRI. This investigation enables a deeper understanding of the neural mechanisms underlying motor control, interhemispheric interactions, and the role of the ipsilateral motor network. It sheds light on the intricate relationship between neural excitability and functional brain activity.
The present study aimed to build upon our prior research and investigate the neural network interactions within the motor cortex during various task contexts, while also examining their association with iMEP as depicted in Figure 1. We chose three conditions that are comparable to previous studies (Tazoe and Perez, 2014; Seusing et al., 2023) to access whether unilateral, bilateral homologous and bilateral heterologous/antagonistic contexts would reveal differential modulation based on the ROIs (M1, PMv, PMd, SMA) defined by Human motor area template (HMAT). We estimated directed functional connectivity to examine the information flow between these regions using time-resolved partial directed coherence (TPDC) based on fMRI data (Anwar et al., 2013). Finally, we employed a correlation analysis between the directed connectivity and iMEP obtained from TMS (Seusing et al., 2023) to explore the relationship between iMEP and network connectivity. We hypothesis that the connectivity patterns within the motor cortex would reflect directed neural modulations under different task contexts. Furthermore, we anticipate that the involvement of ipsilateral motor areas in network modulation would be more pronounced during complex motor tasks and demonstrate association with iMEP specifically within the complex task contexts.
FIGURE 1. Experiment and analysis pipeline. 1. Task based MRI experiment (top: bilateral flexion; middle: unilateral flexion; bottom: bilateral flexion-extension) 2. MRI pre-processing pipeline, head motion parameters, WM and CSF were used in the covariate regression 3. Time series extracted from ROls. 4 Directed functional connectivity estimation using TPDC based on extracted data. 5 Statistical analysis including ANOVA and correlation analysis. Abbreviation: WM white matter, CSF Cerebrospinal fluid, TS time series, HMAT human motor area template, ROI region of interest, M1 primary motor cortex, SMA supplementary motor area, PMd dorsal premotor, PMv ventral premotor. TPDC time-resolved partial directed coherence.
2 Materials and methods
2.1 Participant
Twenty healthy volunteers (26.1 ± 4.6 years old, 10 female) enrolled in the present study. All participants were right-handed according to Edinburgh Handedness inventory (mean score 94.3) (Oldfield 1971). None of the participants had any neurological disorder or took any medication. All participants gave written informed consent to the experimental procedures, which were approved by the local ethics committee at the University Medicine of Greifswald (BB139/18).
2.2 Experiment procedure
The experiments were implemented as blocked and event-related designs during the MRI scanning. The blocked design was 12 min, consisting of interleaving 12 task blocks and 11 rest blocks with a 30 s duration (Figure 1). The participants performed three different motor tasks: 1) bilateral flexion (homologous contraction): bilateral contraction of both biceps brachii; 2) unilateral flexion: contraction of left biceps brachii and relaxation of the right arm; 3) bilateral flexion-extension (antagonistic contraction): contraction of the left biceps brachii and extension of the right arm. Participants were instructed to perform an isometric contraction or extension at 50% of their maximal voluntary movement and keep their heads straight during the experiment. Muscle force during the MRI was monitored visually and with a dynamometer. Real-time visual feedback was provided to the participants on a screen display to allow them to observe their muscle activity. In addition, a dynamometer was used to measure and quantify the actual force exerted by the muscles during the fMRI sessions. This dual approach provided a comprehensive means of monitoring and controlling the muscle force throughout the study.
2.3 iMEP acquisition
The present study included participants who had participated in our prior investigation, which involved a TMS employing identical motor tasks as described in the perimental procedure section, with the aim of obtaining iMEP. During the TMS experiment, participants were seated comfortably in a chair and connected to an electromyography (EMG) system, while also undergoing frameless neuronavigation for registration purposes. The EMG activity was recorded from the left and right biceps brachii muscles using surface electrodes arranged in a tendon-belly-montage configuration, utilizing 10 mm Ag/AgCl electrodes. The recorded EMG signals were subsequently amplified (CED 1902; Cambridge Electronic Design, United Kingdom), band-pass filtered within the range of 20–1,000 Hz, and sampled at a rate of 2 kHz (CED 1401). The data were stored for offline analysis using Signal software (V6.0, CED). The figure-eight coil was applied tangentially to the scalp, oriented at a 45° angle to induce current flow in a posterior to anterior direction. TMS was delivered starting from the anatomical hand knob landmark. Throughout the experiment, participants were instructed to maintain a straight head position at all times and were encouraged to perform an isometric contraction at 50% of their maximum voluntary contraction level in the biceps brachii muscle on the side ipsilateral to the stimulation. The muscle force exerted was visually monitored using a dynamometer, and participants were granted sufficient rest periods within the tasks as needed. For a more comprehensive description of the TMS experiment and setup, we direct the reader to consult the original publication (Seusing et al., 2023).
2.4 MRI acquisition
Whole-brain imaging data from all subjects were collected on a 3T scanner (Verio, Siemens, Erlangen, Germany) with a 32-channel head coil. T1-weighted images were acquired using three-dimensional magnetization-prepared rapid gradient-echo imaging with following parameters: repetition time (TR) = 1.69 m; echo time (TE) = 2.52 m; field of view (FoV) = 256 × 256 mm; number of slices = 176, voxel size = 1 × 1 × 1 mm3; fMRI data consisted of 360 volumes using single-shot echo-planar imaging (EPI) sequence with TR = 2000 m; TE = 30 m ; FA = 90°; FoV = 256 mm ;voxel size = 3 × 3 × 3 mm3.
2.5 MRI pre-processing
fMRI data pre-processing steps were performed following Statistical Parametric Mapping (SPM12) guidelines (Esteban et al., 2019; Diaz et al., 2021). For each subject, the first five volumes were removed to account for T1 relaxation effects. Realignment was performed to remove head motion using rigid body translation and rotation followed by slice timing correction and co-registration with T1 image. Scans were then spatially normalized to Montreal neurological institute (MNI) space using the deformation matrices obtained during MRI pre-processing using the CAT12 toolbox (Structural Brain Mapping Group, Jena University Hospital, Jena, Germany) (Gaser and Dahnke, 2016). The scans were smoothed by convolving with Gaussian kernel of fixed width (6 mm full-width half maximum kernel) to suppress noise. In the end, the nuisance covariate regression was performed including six motion correction parameters and averaged white matter (WM) and cerebrospinal fluid signals (CSF). The fMRI time series data underwent visual inspection by experienced neuroscientists to ensure its quality for subsequent analysis. It worth to mention that we did not employ a data-driven artifacts removal method, such as ICA-AROMA (Pruim et al., 2015) in our preprocessing pipeline. This decision was based on our experimental design, which incorporated short task performance runs, thereby reducing the impact of motion artifacts. Consequently, this approach can avoid the aggressive removal of meaningful signal relevant to our task from the independent components. Thus, we considered the inclusion of motion parameters regression, including WM and CSF, as adequate for this study. This regression approach has been widely utilized in fMRI connectivity studies in previous research (Blasi et al., 2020; Fleischer et al., 2020; Diaz et al., 2021; Ding et al., 2022) and was found to better capture temporal degrees of freedom compared to other advanced artifacts removal methods (Blasi et al., 2020). Finally, the time series extraction based on HMAT (Archer et al., 2018) atlas was carried out. Based on previous findings in cortical physiology of limb movements (Bundy and Leuthardt, 2019), our analysis involves both hemispheres and focuses on cortical areas, namely, M1, PMv, PMd, SMA.
2.6 Time solved partial directed coherence
In the present study, we used time-resolved partial directed coherence (TPDC) method to estimate directed functional connectivity estimation. This method enables us to focus on the temporal dynamics of a signal and analyze directional influence at any specific frequency band. It relies on dual-extended Kalman filtering (DEKF) (Wan and Nelson, 2001), which is one of the most widely used estimation algorithms for nonlinear systems, to calculate time-dependent multivariate autoregressive (MVAR) coefficients at each time point. Briefly, one extended Kalman filter (EKF) estimates the states and feeds this information to the other. The second EKF estimates the model parameters and shares this information with the first. By using two Kalman filters working in parallel, we can estimate both the states and model parameters of the system at each time instant. For detailed DEKF algorithm, we refer reader to the original publication (Wan and Nelson, 2001). Subsequently, Fourier transformation was applied to the MVAR coefficients and computed partial directed coherence (PDC) for each time point. PDC from time-series xj to xi at each time point can be calculated by:
where Aij refers to the Fourier transformed MVAR coefficients matrix and N refers to the number of the connections. In the fMRI time series, we extracted the frequency band of interests from 0.009 to 0.08 Hz and averaged across each time point to obtain robust connectivity values between different motor areas. The choice of this frequency range is based on several factors. Firstly, it is known that neuronal activity in the brain exhibits low-frequency oscillations (<0.1 Hz) and believed to reflect functional networks that are active in that region during task performance (Biswal et al., 1995; Lowe et al., 1998; Cordes et al., 2000; Raichle et al., 2001; Greicius et al., 2004; Roy et al., 2009; Biswal et al., 2010). Secondly, the BOLD signal is relatively slow, with changes occurring over several seconds, which is believed to be related to underlying neuronal activity and functional connectivity in the brain (Fox and Raichle, 2007). The low-frequency range (0.009–0.08 Hz) has been used in numerous fMRI studies to investigate functional connectivity in the human brain (Biswal et al., 1995; Fox et al., 2005; Fox and Raichle, 2007; Ding et al., 2022). Since the precise distribution of the MVAR coefficients is not known, we used a surrogate approach to check the significance of the results. In short, we randomly shuffled the order of these task blocks to create a new time series. The TPDC value is calculated based on a randomly shuffled time series for 1,000 times and the 95th percentile of the connectivity value was taken as the significance threshold. This process is performed separately for each subject.
2.7 Statistical analysis
Statistical analysis of the connectivity data was performed in RStudio (R version 4.1.2). Each connectivity was compared between conditions using non-parametric repeated measure ANOVA. The level of statistical significance against the null-hypothesis was set to p < 0.05 (two-tailed). Wilcoxon test was performed as post hoc method with Bonferroni correction.
2.8 Correlation analysis
In order to examine the correlation between the task-related directed functional connectivity and the iMEP parameters, Spearman’s correlation (Cohen et al., 2009) was performed on the directed functional connectivity and the iMEP measurements (iMEP duration, amplitude, and area) obtained from our recent publication. The level of statistical significance against the null hypothesis was set to p < 0.01. It is worth mentioning that the iMEP area ratio was calculated using the following formula: [area of rectified electromyography (EMG) in iMEP duration/(mean prestimulus EMG*iMEP duration) *100] (Tazoe and Perez, 2014), thus iMEP area ratio expressing the relative size of iMEP compared to the prestimulus EMG. For more detailed iMEP data processing, we refer the reader to the original publication (Seusing et al., 2023).
3 Results
3.1 Ipsilateral motor evoked potential parameters
As shown in Table 1, unilateral flexion shows the highest mean iMEP duration, iMEP area and area ratio. The highest iMEP amplitude was observed during bilateral flexion-extension task.
TABLE 1. iMEP parameters under different motor tasks. Abbreviation: IMEP (ipsilateral motor evoked potential), ms millisecond, μV μV, ‰ permille, % percentage.
3.2 Difference of neural modulation in unilateral flexion, bilateral flexion-extension, in comparison with bilateral flexion
Within selected motor areas, non-parametric repeated ANOVA revealed six directed functional connectivity that were statistically significant between three different motor tasks as shown in Figure 2. In the unilateral flexion task, where the left arm performed elbow flexion and the right arm was at rest, we observed a significant increase in information flow from the left PMv to both sides M1 (left M1: p = 0.0014; right M1: p = 0.0043, Bonferroni corrected) compared to bilateral contraction. Likewise, the information flow from right PMv to both left SMA (p = 0.0408, Bonferroni corrected) and right SMA (p = 0.3624, Bonferroni corrected) were found increased as well. Moreover, ipsilateral connectivity from left M1 to left PMv (p = 0.0362, Bonferroni corrected) and inter-hemispheric connectivity from left SMA to right PMd (p = 0.0110, Bonferroni corrected) were also increased compared to bilateral flexion. Furthermore, these two ECs were also found to be significantly higher in bilateral flexion-extension task compared with bilateral flexion: left M1 to left PMv (p = 0.0095, Bonferroni corrected) and left SMA to right PMd (p = 0.0249, Bonferroni corrected) as shown in Figure 2B. The detailed results are summarized in the Supplementary Table S1.
FIGURE 2. Statistical significance summary of directed functional connectivity between different motor tasks. (A) Significant directed functional connectivity between different motor tasks revealed by Wilcoxon test (Bonferroni corrected). In this illustration, the arrows stroked with white colour indicates the distinct directed functional connectivity between unilateral and bilateral flexion, while a golden stroked arrow indicates both bilateral flexion-extension and unilateral flexion as distinct from bilateral flexion. (B) Boxplot of significant directed functional connectivity summary. Abbreviations: Ml: primary motor cortex; PMv: ventral premotor area; PMd: dorsal premotor area; SMA: supplementary motor area.; white line (within the box): median level. (*): p < 0.05; (**): p < 0.01.
3.3 Difference between unilateral flexion and bilateral flexion-extension
To further reveal the difference in connectivity between unilateral flexion and bilateral flexion-extension, post hoc analysis with Wilcoxon test revealed a decreased level of inter-hemispheric connectivity from left PMv to right PMv (p = 0.0283, Bonferroni corrected) and right M1 to left PMd (p = 0.0304, Bonferroni corrected) during bilateral flexion-extension compared with unilateral flexion as shown in Figure 3.
FIGURE 3. Different neural modulation between unilateral flexion and bilateral flexion-extension. (A) Significant directed functional connectivity between unilateral flexion and bilateral flexion-extension task revealed by Wilcoxon test (Bonferroni corrected). (B) Boxplot of significant directed functional connectivity summary. Abbreviations: Ml: primary motor cortex; PMv: ventral premotor area; PMd: dorsal premotor area; SMA: supplementary motor area.; white line (within the box): median level. (*): p < 0.05; (**): p < 0.01.
3.4 Correlation between directed functional connectivity and iMEP parameters
Figure 4 shows a negative correlation between directed functional connectivity in right M1 to right SMA and iMEP area ratio (p = 0.0094, Bonferroni correction) during bilateral flexion-extension task. No statistical significance was observed during bilateral and unilateral flexion task. The rest correlation results are summarized in the (Supplementary Figure S1).
FIGURE 4. Spearman correlation analysis between directed functional connectivity and iMEP parameters. (A): Scatter plot shows the significant correlation between iMEP area ratio and the connectivity (right M1 to the right SMA). (B): Boxplot shows the iMEP area ratio during different motor tasks.
4 Discussion
Our study contributes to expanding our understanding of the modulatory mechanism involved in voluntary movements. By demonstrating the involvement of ipsilateral motor network and highlighting functional connectivity related MEP. We provide novel insights into the active and regulatory role of this network across different motor tasks. Consequently, our findings enhance our understanding of the intricate dynamics that govern voluntary movements.
The neural modulations within the motor cortex during different motor tasks are not yet fully understood. Precisely quantifying directed functional connectivity, which is based on interactions between relevant brain regions within the motor cortex, can provide a deeper understanding of the function of the ipsilateral hemisphere in voluntary movements. Previous research suggested that the iMEP may be task-specific and spatially distinct from the contralateral hotspot (Seusing et al., 2023). However, our study identified distinct connectivity patterns for bilateral, unilateral, and bilateral flexion-extension movements involving both hemispheres. Additionally, a specific directed functional connectivity pattern (right M1 to right SMA) derived from fMRI time series was found to have a negative correlation with the MEP area ratio, offering further insights into the complex motor networks involved.
The premotor cortex has long been studied in primates, yet the functions of the premotor cortex are diverse and not fully understood (Wise, 1985). In primates, PMv projects to the hand and arm fields of M1, which contains the largest and most detailed representation of hand movements among all cortical motor areas. Additionally, its sensory representation has been shown to contribute to the planning and execution of movements (Kandel et al., 2012). Between bilateral and unilateral flexion, we found that not only the contralateral but also the ipsilateral M1 connections desynchronized from PMv as the strength of the information flow decreased. This highlights the engagement of ipsilateral M1 in the representation of biceps brachii contraction. Previous studies reported that the involvement of ipsilateral M1 in hand movement is often associated with more complex tasks (Chen et al., 1997a; Chen et al., 1997b). However, recent studies in both animals and humans have confirmed the conventional view that ipsilateral projections from M1 to the upper limb are directed to truncal and proximal muscles, which are mainly located in the upper arm. The observation of connectivity from the left PMv to both sides of M1 indicates that M1 functions bilaterally while still being driven by the PMv, even during unilateral flexion task. Furthermore, the lower median level of connectivity towards right M1 (left PMv to right M1) compared with left M1 (left PMv to left M1) might be the indicative of laterality (i.e., as in our study in right-handed individuals).
Interestingly, our results highlighted another bi-hemispheric relationship from right PMV to both sides of SMA as the level of connectivity was significantly decreased during the bilateral flexion task compared with unilateral flexion. In SMA, neurons often relate to complex actions such as postural adjustments and movements that require switching between actions, plans and strategies (Kandel et al., 2012). The projection from PMv to SMA was recently reported during grasping movements. This study proposed a role for integration for force variation at the high-level processing stage (Haller et al., 2009). Due to its direct projections to the hand and motor neurons (Dum and Strick, 1996; Maier et al., 2002), SMA can influence the motor programming without being transmitted to M1 and PMv (Bencivenga et al., 2021). In this study, we observed a decreased connectivity (right PMv to bilateral SMA) during bilateral flexion compared to unilateral flexion, which indicates less modulation recruited from right PMv, suggesting an independent functioning during biceps brachii contraction. Given the role of SMA mainly involved in sensory processing, cognitive and fine motor tasks, the right PMv seem to drive the continuous adjustment and coordination of the voluntary movements.
In addition, we identified two connectivity (left M1 to left PMv; left SMA to right PMd as shown in Figure 2) that distinguish the right arm at flexion, rest and extension, while left arm is in flexion all the time. The projections between PMv and M1 in motor representation was documented in previous studies (Dum and Strick, 2005; Sel et al., 2021). M1 and PM including both dorsal and ventral areas, were previously thought to form a densely interconnected network of cortical areas involved in the generation and control of hand movements. In a Go and no-go motor task, both projections from PMv to M1 and from M1 to PMv were observed, indicating a neural modulation in the reciprocal direction. Likewise, our results demonstrate a higher information flow from M1 to PMv during unilateral flexion and bilateral flexion-extension tasks compared with bilateral flexion, suggesting not only the proxy role of PMv driving M1 as discussed above, but might also play a role in voluntary movements that are dependent on motor context as well.
Furthermore, another inter-hemispheric connection, from left SMA to right PMd, varied in response to changes in the task context. This shows the importance and sensitivity of the PMd in directed functional connectivity analysis, as it was influenced by SMA. In early studies, single-pulse TMS interfered with ipsilateral PMd in patients with acquired lesions, showing an association with the degree of motor impairment and involvement of complex performance (Johansen-Berg et al., 2002; Lotze et al., 2006; Groppa et al., 2012). This association might indicate a compensatory role of ipsilateral PMd in patients with a movement disorder (Stoeckel and Binkofski, 2010). These findings underscore not only the involvement of ipsilateral PMd in the complex task but also a secondary role in the movement when there is dysfunction in the motor network. In addition, one study in finger movement suggests a relationship between SMA and right PMd to relate to bimanual coordination. A recent study (Verstraelen et al., 2021) assessed the planning and executing of bimanual movements via rTMS indicating left PMd increases movement speed and improves movement accuracy whereas right PMd induced deterioration of movement stability. Their results suggest right PMd fulfils a role in continuous adjustment processes of movement. Thus, the involvement of inter-hemispheric connectivity between SMA and PMd is likely to be maintaining an ongoing movement which requires continuous adjustment.
A number of previous studies have documented a reciprocal and interhemispheric relationship between the left and the right PMv (Rushworth et al., 2003; Fujiyama et al., 2016; Verstraelen et al., 2021). As compared to unilateral and bilateral flexion, the bilateral flexion-extension task reduced unidirectional connectivity (left PMv to right PMv), suggesting a specific lateralized modulation during the extension context. It is noteworthy that this connectivity did not differ between bilateral and unilateral flexion. Thus, it is likely to be a distinct feature between homologous and antagonistic contraction.
As mentioned in the above section, PMd was shown to possess the potential to compensate for motor dysfunction and involvement in motor planning. However, the function of PMd is lateralized as left PMd is dominant in motor planning regardless of which hand is moved (Schluter et al., 1998; Schluter et al., 2001; Rushworth et al., 2003; Fujiyama et al., 2013; Fujiyama et al., 2016), whereas right PMd is modulated by task complexity (Groppa et al., 2012). Given the association between M1 and PMd and bimanual performance (Fujiyama et al., 2016; Babaeeghazvini et al., 2019), our finding in the reduction of connectivity from right M1 to left PMd during bilateral flexion-extension suggests an extensive adjustment in the ongoing movement compared with bilateral flexion.
The iMEP parameters derived from TMS experiments exhibited a trend that was similar to the directed functional connectivity findings. Specifically, unilateral flexion demonstrated a greatest iMEP area ratio than the other tasks. Moreover, a noteworthy finding was that the ipsilateral directed functional connectivity from right M1 to right SMA was negatively associated with the iMEP area ratio during the bilateral flexion-extension task, which exhibited an iMEP area ratio higher than that in bilateral flexion but lower than that in unilateral flexion. This negative correlation implies that the modulation from M1 to SMA affects the level of ipsilateral M1 activation during the bilateral antagonistic contraction task, and that the iMEP area ratio may rise exponentially at certain levels of directed functional connectivity. It is noteworthy that no significant correlation was observed between directed functional connectivity and MEP in the context of bilateral and unilateral flexion tasks. This lack of association may be explained by the possibility that the bilateral flexion-extension task necessitates greater interhemispheric communication and coordination, resulting in heightened dependence on the ipsilateral hemisphere and increased modulation from M1 to SMA. Moreover, the unique neural mechanisms engaged in the execution of the bilateral flexion-extension task may contribute to the observed correlation. Nonetheless, additional investigation is warranted to gain a more comprehensive understanding of the underlying mechanisms driving this particular relationship.
This study has several limitations that should be considered when interpreting the results. The directed connectivity methods based on Granger causality require high-temporal resolution data. However, the study relies on low-temporal resolution fMRI data, a significant limitation of the presented results. In addition, fMRI relies on the assumption that the BOLD signal is a reliable proxy for neuronal activity. There are several factors can influence the BOLD signal, such as vascular changes, motion artifacts, and individual variations in neurovascular coupling. These factors can introduce noise and confounds into the fMRI data, potentially affecting the interpretation of the results. In this study, we focus on the low-frequency range commonly used in fMRI connectivity studies, which avoids more susceptible noise and physiological artifacts at high-frequency range. Nevertheless, there is growing evidence demonstrating meaningful network patterns at high-frequency range via high-speed fMRI (Lin et al., 2015; Trapp et al., 2018; Cabral et al., 2023). Therefore, exploring brain network organization in the high-frequency range during upper limb movement might offer valuable insights into the underlying neuronal mechanism, providing complementary information to enhance our understanding. As exploring high-frequency range connectivity is beyond the scope of the present study, further research is needed to fully investigate and elucidate the role of high-frequency components in the context of upper limb motor control.
In summary, the findings of our study shed light on the intricate neural mechanisms underlying upper-limb movements and the role of the ipsilateral motor region in voluntary movements. The observed negative correlation between ipsilateral connectivity from M1 to SMA and iMEP area ratio during the bilateral flexion-extension task suggests a potential modulatory effect of this neural pathway on the degree of ipsilateral M1 activation. These results contribute to a deeper understanding of the functional organization of the motor cortex and provide a foundation for future research on the neural mechanisms underlying voluntary movements.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.
Ethics statement
The studies involving human participants were reviewed and approved by local ethics committee at the University Medicine of Greifswald. The patients/participants provided their written informed consent to participate in this study.
Author contributions
HD: data analysis, interpretation, manuscript writing. NS, SS, ML and MG: experimental design and data acquisition, manuscript review and critique. A.R.A.: data interpretation. BN: review and critique. SG and MM: study design, data interpretation, review and critique. All authors contributed to the article and approved the submitted version.
Funding
The article contribution received funding support from multiple sources. It was supported by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under Project-ID 424778381-TRR 295, and MU 4354/1-1 by MM. In addition, the Fondazione Grigioni per il Morbo di Parkinson to MM, also provided financial support for this research.
Acknowledgments
We extend our gratitude to all participants who participated in this study. We also express our thanks to Miss Anna Grzybowska for her invaluable assistance in graphic design for the manuscript. This research was partially carried out using the Mogon supercomputer and with the support of the advisory services provided by Johannes Gutenberg University Mainz (hpc.uni-mainz.de), which is a member of the AHRP and the Gauss Alliance e.V. We are thankful for the computing time made available on Mogon.
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2023.1199338/full#supplementary-material
References
Anwar, A. R., Muthalib, M., Perrey, S., Galka, A., Granert, O., Wolff, S., et al. (2013). “Comparison of causality analysis on simultaneously measured fMRI and NIRS signals during motor tasks,” in 2013 35th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC), Osaka, Japan, 03-07 July 2013. doi:10.1109/EMBC.2013.6610079
Archer, D. B., Vaillancourt, D. E., and Coombes, S. A. (2018). A template and probabilistic atlas of the human sensorimotor tracts using diffusion MRI. Cereb. Cortex 28 (5), 1685–1699. doi:10.1093/cercor/bhx066
Armand, J., and Kuypers, H. G. J. M. (1980). Cells of origin of crossed and uncrossed corticospinal fibers in the cat: A quantitative horseradish peroxidase study. Exp. Brain Res. 40 (1), 23–34. doi:10.1007/BF00236659
Babaeeghazvini, P., Rueda-Delgado, L. M., Adab, H. Z., Gooijers, J., Swinnen, S., and Daffertshofer, A. (2019). A combined diffusion-weighted and electroencephalography study on age-related differences in connectivity in the motor network during bimanual performance. Hum. Brain Mapp. 40 (6), 1799–1813. doi:10.1002/hbm.24491
Bawa, P., Hamm, J. D., Dhillon, P., and Gross, P. A. (2004). Bilateral responses of upper limb muscles to transcranial magnetic stimulation in human subjects. Exp. Brain Res. 158 (3), 385–390. doi:10.1007/s00221-004-2031-x
Bencivenga, F., Sulpizio, V., Tullo, M. G., and Galati, G. (2021). Assessing the effective connectivity of premotor areas during real vs imagined grasping: A DCM-PEB approach. Neuroimage 230, 117806. doi:10.1016/j.neuroimage.2021.117806
Biswal, B., Yetkin, F. Z., Haughton, V. M., and Hyde, J. S. (1995). Functional connectivity in the motor cortex of resting human brain using echo-planar mri. Magn. Reson Med. 34 (4), 537–541. doi:10.1002/mrm.1910340409
Biswal, B. B., Eldreth, D. A., Motes, M. A., and Rypma, B. (2010). Task-dependent individual differences in prefrontal connectivity. Cereb. Cortex N. Y. N. Y. 20 (9), 2188–2197. doi:10.1093/cercor/bhp284
Blasi, B. D., Caciagli, L., Storti, S. F., Galovic, M., Koepp, M., Menegaz, G., et al. (2020). Noise removal in resting-state and task fMRI: Functional connectivity and activation maps. J. Neural Eng. 17 (4), 046040. doi:10.1088/1741-2552/aba5cc
Buetefisch, C. M., Revill, K. P., Shuster, L., Hines, B., and Parsons, M. (2014). Motor demand-dependent activation of ipsilateral motor cortex. J. Neurophysiol. 112 (4), 999–1009. doi:10.1152/jn.00110.2014
Bundy, D. T., and Leuthardt, E. C. (2019). The cortical physiology of ipsilateral limb movements. Trends Neurosci. 42 (11), 825–839. doi:10.1016/j.tins.2019.08.008
Cabral, J., Fernandes, F. F., and Shemesh, N. (2023). Intrinsic macroscale oscillatory modes driving long range functional connectivity in female rat brains detected by ultrafast fMRI. Nat. Commun. 14 (1), 375. doi:10.1038/s41467-023-36025-x
Chen, R., Cohen, L. G., and Hallett, M. (1997a). Role of the ipsilateral motor cortex in voluntary movement. Can. J. Neurol. Sci. J. Can. Des. Sci. Neurol. 24 (4), 284–291. doi:10.1017/s0317167100032947
Chen, R., Gerloff, C., Hallett, M., and Cohen, L. G. (1997b). Involvement of the ipsilateral motor cortex in finger movements of different complexities. Ann. Neurol. 41 (2), 247–254. doi:10.1002/ana.410410216
Chen, R., Yung, D., and Li, J. Y. (2003). Organization of ipsilateral excitatory and inhibitory pathways in the human motor cortex. J. Neurophysiol. 89 (3), 1256–1264. doi:10.1152/jn.00950.2002
Cohen, I., Huang, Y., Chen, J., and Benesty, J. (2009). Noise reduction in speech processing. Springer Science and Business Media, 1–4.
Cordes, D., Haughton, V. M., Arfanakis, K., Wendt, G. J., Turski, P. A., Moritz, C. H., et al. (2000). Mapping functionally related regions of brain with functional connectivity MR imaging. Ajnr Am. J. Neuroradiol. 21 (9), 1636–1644.
Diaz, C. V. T., González-Escamilla, G., Ciolac, D., García, M. N., Rivas, P. P., Sola, R. G., et al. (2021). Network substrates of centromedian nucleus deep brain stimulation in generalized pharmacoresistant epilepsy. Neurotherapeutics 18, 1665–1677. doi:10.1007/s13311-021-01057-y
Ding, H., Droby, A., Anwar, A. R., Bange, M., Hausdorff, J. M., Nasseroleslami, B., et al. (2022). Treadmill training in Parkinson’s disease is underpinned by the interregional connectivity in cortical-subcortical network. Npj Park Dis. 8 (1), 153. doi:10.1038/s41531-022-00427-3
Dum, R. P., and Strick, P. L. (2005). Frontal lobe inputs to the digit representations of the motor areas on the lateral surface of the hemisphere. J. Neurosci. 25 (6), 1375–1386. doi:10.1523/JNEUROSCI.3902-04.2005
Dum, R. P., and Strick, P. L. (1996). Spinal cord terminations of the medial wall motor areas in macaque monkeys. J. Neurosci. 16 (20), 6513–6525. doi:10.1523/JNEUROSCI.16-20-06513.1996
Duncan, P. W., Goldstein, L. B., Matchar, D., Divine, G. W., and Feussner, J. (2018). Measurement of motor recovery after stroke. Outcome assessment and sample size requirements. Stroke 23 (8), 1084–1089. doi:10.1161/01.str.23.8.1084
Esteban, O., Markiewicz, C. J., Blair, R. W., Moodie, C. A., Isik, A. I., Erramuzpe, A., et al. (2019). fMRIPrep: a robust preprocessing pipeline for functional MRI. Nat. methods 16 (1), 111–116. doi:10.1038/s41592-018-0235-4
Fleischer, V., Muthuraman, M., Anwar, A. R., Gonzalez-Escamilla, G., Radetz, A., Gracien, R. M., et al. (2020). Continuous reorganization of cortical information flow in multiple sclerosis: A longitudinal fMRI effective connectivity study. Sci. Rep-uk 10 (1), 806. doi:10.1038/s41598-020-57895-x
Fleming, G. W. T. H. (1938). Somatic motor and sensory representation in the cerebral cortex of man as studied by electrical stimulation (Brain, vol Lx, p 389, dec, 1937) penfield W, and boldrey E. J. Ment. Sci. 84 (352), 868. doi:10.1192/bjp.84.352.868-a
Fox, M. D., and Raichle, M. E. (2007). Spontaneous fluctuations in brain activity observed with functional magnetic resonance imaging. Nat. Rev. Neurosci. 8 (9), 700–711. doi:10.1038/nrn2201
Fox, M. D., Snyder, A. Z., Vincent, J. L., Corbetta, M., Essen, D. C. V., and Raichle, M. E. (2005). The human brain is intrinsically organized into dynamic, anticorrelated functional networks. Proc. Natl. Acad. Sci. 102 (27), 9673–9678. doi:10.1073/pnas.0504136102
Fujiyama, H., Hinder, M. R., and Summers, J. J. (2013). Functional role of left PMd and left M1 during preparation and execution of left hand movements in older adults. J. Neurophysiol. 110 (5), 1062–1069. doi:10.1152/jn.00075.2013
Fujiyama, H., Soom, J. V., Rens, G., Cuypers, K., Heise, K. F., Levin, O., et al. (2016). Performing two different actions simultaneously: The critical role of interhemispheric interactions during the preparation of bimanual movement. Cortex 77, 141–154. doi:10.1016/j.cortex.2016.02.007
Gaser, C., and Dahnke, R. (2016). Cat - a computational anatomy toolbox for the analysis of structural MRI data.
Grafton, S., Mazziotta, J., Presty, S., Friston, K., Frackowiak, R., and Phelps, M. (1992). Functional anatomy of human procedural learning determined with regional cerebral blood flow and PET. J. Neurosci. 12 (7), 2542–2548. doi:10.1523/JNEUROSCI.12-07-02542.1992
Grefkes, C., Eickhoff, S. B., Nowak, D. A., Dafotakis, M., and Fink, G. R. (2008). Dynamic intra- and interhemispheric interactions during unilateral and bilateral hand movements assessed with fMRI and DCM. Neuroimage 41 (4), 1382–1394. doi:10.1016/j.neuroimage.2008.03.048
Greicius, M. D., Srivastava, G., Reiss, A. L., and Menon, V. (2004). Default-mode network activity distinguishes Alzheimer’s disease from healthy aging: Evidence from functional MRI. Proc. Natl. Acad. Sci. 101 (13), 4637–4642. doi:10.1073/pnas.0308627101
Groppa, S., Schlaak, B. H., Münchau, A., Werner-Petroll, N., Dünnweber, J., Bäumer, T., et al. (2012). The human dorsal premotor cortex facilitates the excitability of ipsilateral primary motor cortex via a short latency cortico-cortical route. Hum. Brain Mapp. 33 (2), 419–430. doi:10.1002/hbm.21221
Haller, S., Chapuis, D., Gassert, R., Burdet, E., and Klarhöfer, M. (2009). Supplementary motor area and anterior intraparietal area integrate fine-graded timing and force control during precision grip. Eur. J. Neurosci. 30 (12), 2401–2406. doi:10.1111/j.1460-9568.2009.07003.x
Johansen-Berg, H., Rushworth, M. F. S., Bogdanovic, M. D., Kischka, U., Wimalaratna, S., and Matthews, P. M. (2002). The role of ipsilateral premotor cortex in hand movement after stroke. Proc. Natl. Acad. Sci. 99 (22), 14518–14523. doi:10.1073/pnas.222536799
Kandel, R. E., Schwartz, H. J., and Jessell, M. T. (2012). Principles of neural science. Fifth edition. New York, N.Y.: McGraw-Hill Education LLC.
Kim, S. G., Ashe, J., Georgopoulos, A. P., Merkle, H., Ellermann, J. M., Menon, R. S., et al. (1993). Functional imaging of human motor cortex at high magnetic field. J. Neurophysiol. 69 (1), 297–302. doi:10.1152/jn.1993.69.1.297
Lin, F. H., Chu, Y. H., Hsu, Y. C., Lin, J. F. L., Tsai, K. W. K., Tsai, S. Y., et al. (2015). Significant feed-forward connectivity revealed by high frequency components of BOLD fMRI signals. Neuroimage 121, 69–77. doi:10.1016/j.neuroimage.2015.07.036
Lotze, M., Markert, J., Sauseng, P., Hoppe, J., Plewnia, C., and Gerloff, C. (2006). The role of multiple contralesional motor areas for complex hand movements after internal capsular lesion. J. Neurosci. Off. J. Soc. Neurosci. 26 (22), 6096–6102. doi:10.1523/JNEUROSCI.4564-05.2006
Lowe, M. J., Mock, B. J., and Sorenson, J. A. (1998). Functional connectivity in single and multislice echoplanar imaging using resting-state fluctuations. Neuroimage 7 (2), 119–132. doi:10.1006/nimg.1997.0315
Maier, M. A., Armand, J., Kirkwood, P. A., Yang, H. W., Davis, J. N., and Lemon, R. N. (2002). Differences in the corticospinal projection from primary motor cortex and supplementary motor area to macaque upper limb motoneurons: An anatomical and electrophysiological study. Cereb. Cortex 12 (3), 281–296. doi:10.1093/cercor/12.3.281
Mayhew, S. D., Porcaro, C., Tecchio, F., and Bagshaw, A. P. (2017). fMRI characterisation of widespread brain networks relevant for behavioural variability in fine hand motor control with and without visual feedback. NeuroImage 148, 330–342. doi:10.1016/j.neuroimage.2017.01.017
Perez, M. A., Butler, J. E., and Taylor, J. L. (2013). Modulation of transcallosal inhibition by bilateral activation of agonist and antagonist proximal arm muscles. J. Neurophysiol. 111 (2), 405–414. doi:10.1152/jn.00322.2013
Pruim, R. H. R., Mennes, M., Rooij, D., Llera, A., Buitelaar, J. K., and Beckmann, C. F. (2015). ICA-AROMA: A robust ICA-based strategy for removing motion artifacts from fMRI data. Neuroimage 112, 267–277. doi:10.1016/j.neuroimage.2015.02.064
Raichle, M. E., MacLeod, A. M., Snyder, A. Z., Powers, W. J., Gusnard, D. A., and Shulman, G. L. (2001). A default mode of brain function. Proc. Natl. Acad. Sci. 98 (2), 676–682. doi:10.1073/pnas.98.2.676
Rao, S. M., Binder, J. R., Bandettini, P. A., Hammeke, T. A., Yetkin, F. Z., Jesmanowicz, A., et al. (1993). Functional magnetic resonance imaging of complex human movements. Neurology 43 (11), 2311–2318. doi:10.1212/wnl.43.11.2311
Roy, A. K., Shehzad, Z., Margulies, D. S., Kelly, A. M. C., Uddin, L. Q., Gotimer, K., et al. (2009). Functional connectivity of the human amygdala using resting state fMRI. Neuroimage 45 (2), 614–626. doi:10.1016/j.neuroimage.2008.11.030
Rushworth, M. F. S., Johansen-Berg, H., Göbel, S. M., and Devlin, J. T. (2003). The left parietal and premotor cortices: Motor attention and selection. Neuroimage 20, S89–S100. doi:10.1016/j.neuroimage.2003.09.011
Schluter, N. D., Krams, M., Rushworth, M. F. S., and Passingham, R. E. (2001). Cerebral dominance for action in the human brain: The selection of actions. Neuropsychologia 39 (2), 105–113. doi:10.1016/s0028-3932(00)00105-6
Schluter, N. D., Rushworth, M. F., Passingham, R. E., and Mills, K. R. (1998). Temporary interference in human lateral premotor cortex suggests dominance for the selection of movements. A study using transcranial magnetic stimulation. Brain 121 (5), 785–799. doi:10.1093/brain/121.5.785
Sel, A., Verhagen, L., Angerer, K., David, R., Klein-Flügge, M. C., and Rushworth, M. F. S. (2021). Increasing and decreasing interregional brain coupling increases and decreases oscillatory activity in the human brain. P Natl. Acad. Sci. U. S. A. 118 (37), e2100652118. doi:10.1073/pnas.2100652118
Seusing, N., Strauss, S., Fleischmann, R., Nafz, C., Groppa, S., Muthuraman, M., et al. (2023). The excitability of ipsilateral motor evoked potentials is not task-specific and spatially distinct from the contralateral motor hotspot. Exp. Brain Res. (Under review).
Stoeckel, M. C., and Binkofski, F. (2010). The role of ipsilateral primary motor cortex in movement control and recovery from brain damage. Exp. Neurol. 221 (1), 13–17. doi:10.1016/j.expneurol.2009.10.021
Tazoe, T., and Perez, M. A. (2014). Selective activation of ipsilateral motor pathways in intact humans. J. Neurosci. 34 (42), 13924–13934. doi:10.1523/JNEUROSCI.1648-14.2014
Trapp, C., Vakamudi, K., and Posse, S. (2018). On the detection of high frequency correlations in resting state fMRI. Neuroimage 164, 202–213. doi:10.1016/j.neuroimage.2017.01.059
Verstraelen, S., Dun, K., Depestele, S., Hoornweder, S. V., Jamil, A., Ghasemian-Shirvan, E., et al. (2021). Dissociating the causal role of left and right dorsal premotor cortices in planning and executing bimanual movements – a neuro-navigated rTMS study. Brain Stimul. 14 (2), 423–434. doi:10.1016/j.brs.2021.02.006
Volz, L. J., Eickhoff, S. B., Pool, E. M., Fink, G. R., and Grefkes, C. (2015). Differential modulation of motor network connectivity during movements of the upper and lower limbs. Neuroimage 119, 44–53. doi:10.1016/j.neuroimage.2015.05.101
Wan, E. A., and Nelson, A. T. (2001). “Dual extended Kalman filter methods,” in Kalman filtering and neural networks. Editor S. Haykin (John Wiley and Sons), 123.
Wassermann, E. M., Pascual-Leone, A., and Hallett, M. (1994). Cortical motor representation of the ipsilateral hand and arm. Exp. Brain Res. 100 (1), 121–132. doi:10.1007/BF00227284
Wise, S. P. (1985). The primate premotor cortex: Past, present, and preparatory. Annu. Rev. Neurosci. 8 (1), 1–19. doi:10.1146/annurev.ne.08.030185.000245
Keywords: ipsilateral motor network, upper limb, ipsilateral motor evoked potential, voluntary movement, bilateral motor network
Citation: Ding H, Seusing N, Nasseroleslami B, Anwar AR, Strauss S, Lotze M, Grothe M, Groppa S and Muthuraman M (2023) The role of ipsilateral motor network in upper limb movement. Front. Physiol. 14:1199338. doi: 10.3389/fphys.2023.1199338
Received: 04 April 2023; Accepted: 23 June 2023;
Published: 03 July 2023.
Edited by:
Gisela E. Hagberg, University of Tübingen, GermanyReviewed by:
Vinod Kumar, Max Planck Institute for Biological Cybernetics, GermanyFranca Tecchio, National Research Council (CNR), Italy
Copyright © 2023 Ding, Seusing, Nasseroleslami, Anwar, Strauss, Lotze, Grothe, Groppa and Muthuraman. 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: Sergiu Groppa, c2Vncm9wcGFAdW5pLW1haW56LmRl; Muthuraman Muthuraman, TXV0aHVyYW1hbl9NQHVrdy5kZQ==
†These authors have contributed equally to this work and share last authorship