Next Article in Journal
Alterations of Thalamic Nuclei Volumes and the Intrinsic Thalamic Structural Network in Patients with Multiple Sclerosis-Related Fatigue
Previous Article in Journal
Corticomotor Plasticity Underlying Priming Effects of Motor Imagery on Force Performance
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Increased Inertia Triggers Linear Responses in Motor Cortices during Large-Extent Movements—A fNIRS Study

1
Department of Rehabilitation Medicine, Ruijin Hospital, Shanghai Jiao Tong University School of Medicine, Shanghai 200025, China
2
School of Medicine, Shanghai Jiao Tong University, Shanghai 200025, China
3
State Key Laboratory of Cognitive Neuroscience and Learning, Beijing Normal University, Beijing 100091, China
*
Author to whom correspondence should be addressed.
Brain Sci. 2022, 12(11), 1539; https://doi.org/10.3390/brainsci12111539
Submission received: 13 October 2022 / Revised: 9 November 2022 / Accepted: 10 November 2022 / Published: 13 November 2022
(This article belongs to the Section Neurorehabilitation)

Abstract

:
Activities of daily living consist of accurate, coordinated movements, which require the upper limbs to constantly interact with environmental loads. The magnitude of the load was shown to affect kinematic outcomes in healthy subjects. Moreover, the increase in load facilitates the recovery of motor function in patients with neurological disorders. Although Brodmann Areas 4 and 6 were found to be active during loaded movements, it remains unclear whether stronger activation can be triggered simply by increasing the load magnitude. If such a linear relationship exists, it may provide a basis for the closed-loop adjustment of treatment plans in neurorehabilitation. Fourteen healthy participants were instructed to lift their hands to their armpits. The movements were grouped in blocks of 25 s. Each block was assigned a magnitude of inertial loads, either 0 pounds (bare hand), 3 pounds, or 15 pounds. Hemodynamic fNIRS signals were recorded throughout the experiment. Both channel-wise and ROI-wise analyses found significant activations against all three magnitudes of inertia. The generalized linear model revealed significant increases in the beta coefficient of 0.001673/pound in BA4 and 0.001338/pound in BA6. The linear trend was stronger in BA6 (conditional r2 = 0.9218) than in BA4 (conditional r2 = 0.8323).

1. Introduction

Accurate, coordinated movements in the upper extremities are essential for activities of daily living (ADL). During ADL tasks such as eating, dressing, and self-grooming [1], the upper limbs interact constantly with loads from food, clothes, or tools. The magnitude of the load determines both the organization of muscle activation [2] and the resulting performance [3] in healthy individuals; moreover, when the magnitude is increased progressively, it facilitates the recovery of motor function in patients with neurological diseases, e.g., stroke [4] and Parkinson’s disease [5,6]. Therefore, gauging the progress of rehabilitation requires an understanding of how the brain reacts to progressively heavier loads, which is a common protocol in clinical rehabilitation [7]. On the other hand, if stronger activation in Brodmann Area 4 (BA4) and Brodmann Area 6 (BA6) can be triggered by adding loads, this potentially provides a viable path toward closed-loop, individualized therapy.
Mounting evidence has delineated the role of BA4 and BA6 for upper-limb movements against loads [8,9]. BA4 mainly contains the primary motor cortex (M1), which assumes the function of motor execution; BA6 mainly contains the supplementary motor cortex (SMA) and the premotor cortex (PMC), which assume the functions of motor planning and motor learning [10]. Experiments with fMRI have identified that the activation volume of contralateral BA4 and BA6 increases with increasing force during hand griping [11], the isometric contraction of the index finger [12], and hand squeezing [13], whereas another study claimed that only contralateral BA4 increased monotonically with force during isometric hand griping [14]. However, most fMRI studies have focused on movements of minimal joint excursions, and the load magnitude has been shown to be relatively small, partly due to the intolerance of head movements during imaging. Among other modalities of imaging, electroencephalogram (EEG) studies have found that the motor-activity-related cortical potential is linearly correlated with the isometric force output of hand grip [15] and elbow-flexion contractions [16]. Positron emission tomography (PET) and magnetoencephalography (MEG) studies reported that the activation of contralateral BA4 increased with the index-finger flexion [17] and extension force [18]. The movements in these studies were also restricted to small extents.
Functional near-infrared spectroscopy (fNIRS) is an emerging technology for brain imaging featuring high ecological validity [19], which makes it suitable for the detection of cortical responses during large-extent movements that may drag the head or trunk [20]. With regards to upper-limb movements, a recent study showed stronger activation in SMA, PMC, and M1 during voluntary movement compared to passive movement [8]. In lower-limb movements, fNIRS detected an exercise-induced decrease of deoxygenated hemoglobin in PMC [21]. Few studies have examined how mechanical load affects fNIRS responses during large-extent movements in the upper limbs.
In this study, we focused on whether increases in load inertia (i.e., the mass of an inertial load) would lead to commensurate responses in motor-related cortices within BA4 and BA6 in healthy subjects. The range, tempo, and mechanical load of the movement were chosen such that the movement was large-extent and meaningful for ADL. Our central hypothesis was that the activation of contralateral BA4 and BA6 would correlate positively with load inertia. We focused on inertial load because it is a ubiquitous type arising from interactions with daily objects, e.g., dumbbells, heavy doors, etc. [22]. This study provides an essential step toward the closed-loop adjustment of treatment plans for the rehabilitation of motor functionalities.

2. Materials and Methods

2.1. Participants

Fourteen healthy subjects participated in the study (10 males, 4 females; age = 26.5 ± 3.2; 13 right-handed, 1 mixed-handed, according to the Edinburgh handedness inventory [23]). None of the participants reported any history of neurological or musculoskeletal disorders. Before the experiment, each participant gave written consent as approved by the Ethics Committee of Ruijin Hospital, School of Medicine, Shanghai Jiao Tong University.

2.2. Task Design

Each subject sat on an armless chair while lifting the right hand from the hanging position to the armpit; the forearm and upper arm were required to move within the sagittal plane (Figure 1A). Subjects lifted and descended their hands repeatedly at 0.5 Hz following a metronome. Before the experiment, all subjects practiced the tasks for 3–5 min. Videos are available in Supplementary Materials for reference regarding the movement task.
The lifting movements were performed under three different load conditions: 0-pound (bare hand), 3-pound (1.36 kg), and 15-pound (6.80 kg). The 3-pound condition was chosen because it was the entry-level weight for strength training [24]; the 15-pound condition was chosen in reference to previous studies [25], and also because this was the maximum weight that allowed pilot subjects to finish the protocol without excessive fatigue.
Movements were grouped in blocks during fNIRS data collection. Each block contained 12–14 movements against a specific load magnitude. Each block lasted 25 s, and the subjects were required to rest for 40 s between adjacent blocks. The completion of all blocks took about 13 min. The sequence of the blocks was randomized. The procedure for each block was as follows:
(1)
The experimenter reminded the subject of the correct dumbbell for the next block;
(2)
The subject picked up the assigned dumbbell (if any);
(3)
The computer played a starting sound;
(4)
The subject performed lifting movements following the metronome tempo;
(5)
The computer played an ending sound;
(6)
The subject put down the dumbbell and rested.
(7)
Brain-hemodynamic data were acquired.
Brain-hemodynamic signals were captured using a multichannel fNIRS device (model ETG-4100, HITACHI Inc., Tokyo, Japan) at a 10-hertz sampling rate. The device evaluated the absorption of near-infrared light at two wavelengths (695 nm and 830 nm) and computed the corresponding hemoglobin and deoxyhemoglobin density following the modified Beer–Lambert Law [26]. ROIs included contralateral BA4 (mainly M1) and BA6 (encompassing SMA and PMC). BA4 was chosen because of its role in motor execution, and BA6 were chosen because of its functionality in motor planning [27]. Both BA4 and BA6 were known to be active during motor-related tasks in previous studies [28,29].
Optodes were placed on the head using a high-density 3 × 10 probe with 16 emitters and 14 detectors, which formed 44 channels covering the BA4 and BA6 of both hemispheres. The configuration of the probe is shown in Figure 1B. During the probe placement, the international 10–20 system [30] was adopted to ensure coverage of the ROIs. Specifically, we first placed the holders over subjects’ heads with channel 23 located at the CZ point, after which we set the front brim of the probe parallel to the coronal plane of the skull. After the scan, the locations of channels were measured using a 3D magnetic-space digitizer (Polhemus Patriot, Polhemus Inc., Colchester, VT, USA) to estimate the anatomical brain region of each channel. The MNI coordinates of each electrode were calculated using NFRI functions [31]. Examples of anatomic labeling and the corresponding probabilistic registrations are listed in Table 1. The flow chart of the experiment is shown in Figure 1C.

2.3. Brain-Hemodynamic-Signal Processing

The NIRS-KIT toolbox (version 2.0) [32] was used for brain hemodynamic signal processing. Oxy-Hb signals were used to quantify cortical activity in this study because they were sensitive to regional-cerebral-blood-flow fluctuations [33]. The pre-processing of oxy-Hb consisted of 3 steps: (1) the raw oxy-Hb signals were detrended using the linear-detrending method; (2) motion artifacts were removed using temporal derivative distribution repair (TDDR); (3) the oxy-Hb signals were band-pass filtered with cut-off frequencies at 0.01 Hz and 0.08 Hz (third-order Butterworth) to remove physiological artifacts (such as heartbeats, breath, and Mayer wave) and high-frequency noise.
The general linear model (GLM) approach [34] was used to characterize the hemodynamic response of oxy-Hb in all 44 channels under 3 conditions. Beta coefficients, of which the sign and magnitude indicate the direction (positive/negative) and intensity of change in oxy-Hb, were generated for each subject, channel, and condition. Beta coefficients are referred to interchangeably with cortical activation in this manuscript unless otherwise specified.

2.4. Statistical Analysis

Beta coefficients were first compared against 0 to detect cortical activation using one-sample two-tailed t-tests with FDR correction. Effects of weight on beta coefficients were analyzed using mixed-effect linear models, as follows:
b e t a   ~   w e i g h t + ( 1 | s u b j e c t )
where weight represents the mass of the dumbbell and the term ( 1 | s u b j e c t ) accounts for subject-specific intercepts due to repeated measures. The significance level in the analyses was set at 0.05. All statistical analyses were performed using R (version 4.2.1). The models were fitted using the ‘lmer’ package [35], and the significances were calculated using the ‘lmerTest’ package [36]. The conditional r-squares were calculated using the ‘MuMIn’ package [37]; conditional r-square is the proportion of total variance explained through both fixed and random effects.

3. Results

We first examined the MNI coordinates and spatial registration of each channel. According to the registration records, the channels most likely to be situated in the contralateral BA4 and BA6 are marked in Figure 2. The 3D locations of all the channels are also depicted in Figure 2. The spatial registrations of all the channels for a representative subject are listed in Table 1.

3.1. Activation Analysis (Channel-Wise)

In the 0-pound condition (Figure 3A), significant activations were detected in five out of seven BA4 channels (7, 8, 16, 17, 18, all FDR-adjusted p < 0.05) as well as two out of seven BA6 channels (26, 35, all FDR-adjusted p < 0.05). In the 3-pound condition (Figure 3B), we discovered that all the BA4-channels (all FDR-adjusted p < 0.05) and six out of seven BA6 channels (24, 25, 26, 34, 35, 36, all FDR-adjusted p < 0.05) were significantly activated. In the 15-pound condition (Figure 3C), all the BA4 channels (all FDR-adjusted p < 0.05) were also significantly activated, and six out of seven BA6 channels (24, 25, 26, 34, 35, and 36, all FDR-adjusted p < 0.05) were activated. The overview of the activation in the ROIs map, interpolated from the averaged beta coefficients, is shown in Figure 4.

3.2. Activation Analysis (ROI-Wise)

To improve the spatial consistency across channels, we employed a group analysis based on ROI [38]. The three channels with the highest probability of registration were chosen for the weight adjustment of signals based on ROI. Therefore, the weight-adjusted hemodynamic signal for each ROI was:
o x y H b R O I = i = 1 3 P i o x y i i = 1 3 P i
where o x y H b R O I represented the oxy-Hb signals of each ROI, and P i represented the likelihood of channel i out of three, ranked by the likelihood in the probabilistic registration.
In the 0-pound condition (Figure 5), both ROIs were significantly activated (BA4: t13 = 2.821, p < 0.05 BA6: t13 = 2.4619, p < 0.05). In the 3-pound condition, both BA4 and BA6 were significantly activated (BA4: t13 = 4.3567, p < 0.001; BA6: t13 = 3.7747, p < 0.01). In the 15-pound condition, the activation was also significant in both BA4 and BA6 (BA4: t13 = 4.6549, p < 0.001; BA6: t13 = 4.3725, p < 0.001).

3.3. Correlation Analysis (Channel-Wise)

We performed a traditional channel-wise analysis to test the relationship between the beta coefficients and the magnitude of load inertia. All the BA6 channels exhibited statistically significant positive correlations between load magnitudes and beta coefficients (p < 0.05); six out of the seven BA4 channels exhibited significant positive correlations between loads and beta coefficients (p < 0.05). Figure 6 shows the linear relationships from the three most significant channels. The channel-wise details of the slope, T-value, p-value, and location are listed in Table 2.

3.4. Correlation Analysis (ROI-Wise)

We ran a linear model in each ROI to examine the impact of load magnitudes on the cortical activation. As shown in Figure 7A, we discovered significantly positive correlations between the load magnitudes and the beta. The beta coefficient increased by 0.001673/pounds (conditional r2 = 0.8323, p < 0.01) in BA4. The regression lines of each subject were displayed to show individual trends. As shown in Figure 7B, the activation level of BA6 also significantly correlated with the magnitude of the load by 0.001338/pounds (conditional r2 = 0.9218, p < 0.0001). BA6 show a stronger linear trend than BA4, as suggested by the mixed-effect linear model.

4. Discussion

In this study, we investigated the correlation between cortical activation and the magnitude of inertial load using fNIRS. The subjects were instructed to lift their right hands to their armpits against three different loads: 0 pounds (bare hand), 3 pounds, and 15 pounds. Contralateral BA4 and BA6 activations were detected in all three conditions. Linear correlations between the beta coefficients and the magnitude of the load were identified in both channel-wise and ROI-wise analyses.
There have been debates on whether the correlation between BA4/BA6 activities and load magnitude is linear. One fMRI study claimed that only the BOLD signal in BA4, but not in BA6, increased linearly with the force output during hand gripping [14]; by contrast, another study reported that the BA6 activity showed monotonic relationships with increased force output [13]. Our results showed that the linear trends were significant in both BA4 and BA6 between the cortical activity and the load inertia. In addition, we added to existing knowledge with the observation that the linear trend seemed stronger in BA6. The difference between our results and previous reports might have been due to the larger extent of the movements in our study, i.e., larger range, higher velocity, and higher loads.
A stronger linear trend was found in BA6 (conditional r2 = 0.9218) than in BA4 (conditional r2 = 0.8323) with increased inertia. This might imply that the main effect of increasing the load inertia was to engage motor planning or learning, the suggested role of SMA and PMC within BA6 [39,40,41], rather than to engage only movement execution, which is arguably the primary role of BA4 [42]. The importance of BA6 in clinical rehabilitation has been confirmed in multiple studies [43,44]. As a result, it is interesting to ask whether a patient shows sufficient responses in the BA6 against each extra kg of load, which is a potential benchmark for clinical rehabilitation.
Our results are compatible with several studies showing that BA4 and BA6 activate more strongly against increased load during upper-limb movements. One fNIRS study reported that resistance exercise induced higher cortical activation than non-resistance exercise [9], but the load perceived by the subject was relatively small and the MNI coordinates of the activation regions were not reported. Another fNIRS study claimed that oxy-Hb exhibited an increasing trend with increasing hand-grip force [45], but the movements were restricted to a small range of motion. In general, this study provided clear data to suggest that both BA4 and BA6 linearly correlated with the magnitude of the load, sampled within the 0–15-pounds range.
In terms of clinical significance, this study might: (1) help physicians to understand how the brain reacts to progressively heavier loads, which are commonly used in clinical rehabilitation [46]; (2) provide a viable path toward the closed-loop and individualized adjustment of treatment plans; (3) provide a potential benchmark for rehabilitation assessment. Our long-term goal requires the current protocol to be generalized to patients with neurological disorders, such as stroke or traumatic brain injury. The limitations of this study include the following: (1) the patients may not be capable of lifting the arm holding the dumbbell, particularly for multiple repetitions; (2) the movement trajectories and joint angles were not motion-captured, which might be essential for the explanation of differences between healthy subjects and patients.

5. Conclusions

We found that during large-extent upper-limb movements, the increased magnitude of the inertial load correlated linearly with the activation in the motor-related cortices BA4 and BA6. The linear trend was stronger in BA6 than in BA4. This raises the question of whether this linear relationship is also present in patients with neurological disorders.

Author Contributions

Conceptualization, C.Z. and C.M.N.; formal analysis, Z.C.; funding acquisition, Q.X. and C.M.N.; investigation, Z.C., Y.Q. and J.Y.; methodology, Z.C., X.S., C.Z., and C.M.N.; project administration, Q.X. and C.M.N.; resources, Q.X.; supervision, C.Z. and C.M.N.; visualization, Z.C.; writing—original draft, Z.C.; writing—review and editing, C.M.N. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China, grant number: 81971722; Shanghai Municipal Health Commission Project, grant number: 2019SY004; Shanghai Science and Technology Commission Project, grant number: 22490711000; Shanghai Municipal Key Clinical Specialty, grant number: shslczdzk02701.

Institutional Review Board Statement

This study was approved by the Ethics Committee of Ruijin Hospital, School of Medicine, Shanghai Jiao Tong University (no. 209 of 2020).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

The authors thank Qiang Xu for his help with the multi-device data acquisition.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Edemekong, P.F.; Bomgaars, D.L.; Sukumaran, S.; Schoo, C. Activities of Daily Living; StatPearls Publishing: Treasure Island, FL, USA, 2022. [Google Scholar]
  2. Gottlieb, G.L. On the voluntary movement of compliant (inertial-viscoelastic) loads by parcellated control mechanisms. J. Neurophysiol. 1996, 76, 3207–3229. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Shapiro, M.B.; Niu, C.M.; Poon, C.; David, F.J.; Corcos, D.M. Proprioceptive feedback during point-to-point arm movements is tuned to the expected dynamics of the task. Exp. Brain Res. 2009, 195, 575–591. [Google Scholar] [CrossRef] [PubMed]
  4. Flansbjer, U.; Miller, M.; Downham, D.; Lexell, J. Progressive resistance training after stroke: Effects on muscle strength, muscle tone, gait performance and perceived participation. Acta. Derm. Venereol. 2008, 40, 42–48. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Lima, L.O.; Scianni, A.; Rodrigues-De-Paula, F. Progressive resistance exercise improves strength and physical performance in people with mild to moderate Parkinson’s disease: A systematic review. J. Physiother. 2013, 59, 7–13. [Google Scholar] [CrossRef] [Green Version]
  6. Corcos, D.M.; Robichaud, J.A.; David, F.; Leurgans, S.E.; Vaillancourt, D.E.; Poon, C.; Rafferty, M.R.; Kohrt, W.M.; Comella, C.L. A two-year randomized controlled trial of progressive resistance exercise for Parkinson’s disease. Mov. Disord. 2013, 28, 1230–1240. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Wist, S.; Clivaz, J.; Sattelmayer, K.M. Muscle strengthening for hemiparesis after stroke: A meta-analysis. Ann. Phys. Rehabil. Med. 2016, 59, 114–124. [Google Scholar] [CrossRef] [PubMed]
  8. Xia, W.; Dai, R.; Xu, X.; Huai, B.; Bai, Z.; Zhang, J.; Jin, M.; Niu, W. Cortical mapping of active and passive upper limb training in stroke patients and healthy people: A functional near-infrared spectroscopy study. Brain Res. 2022, 1788. [Google Scholar] [CrossRef]
  9. Shi, P.; Li, A.; Yu, H. Response of the Cerebral Cortex to Resistance and Non-resistance Exercise Under Different Trajectories: A Functional Near-Infrared Spectroscopy Study. Front. Neurosci. 2021, 15. [Google Scholar] [CrossRef]
  10. Strotzer, M. One Century of Brain Mapping Using Brodmann Areas*. Clin. Neuroradiol. 2009, 19, 179–186. [Google Scholar] [CrossRef]
  11. Dai, T.H.; Liu, J.Z.; Sahgal, V.; Brown, R.W.; Yue, G.H. Relationship between muscle output and functional MRI-measured brain activation. Exp. Brain Res. 2001, 140, 290–300. [Google Scholar] [CrossRef]
  12. van Duinen, H.; Renken, R.; Maurits, N.M.; Zijdewind, I. Relation between muscle and brain activity during isometric contractions of the first dorsal interosseus muscle. Hum. Brain Mapp. 2007, 29, 281–299. [Google Scholar] [CrossRef] [PubMed]
  13. Cramer, S.C.; Weisskoff, R.M.; Schaechter, J.D.; Nelles, G.; Foley, M.; Finklestein, S.P.; Rosen, B.R. Motor cortex activation is related to force of squeezing. Hum. Brain Mapp. 2002, 16, 197–205. [Google Scholar] [CrossRef] [PubMed]
  14. Keisker, B.; Hepp-Reymond, M.-C.; Blickenstorfer, A.; Meyer, M.; Kollias, S.S. Differential force scaling of fine-graded power grip force in the sensorimotor network. Hum. Brain Mapp. 2009, 30, 2453–2465. [Google Scholar] [CrossRef] [PubMed]
  15. Liu, J.Z.; Yang, Q.; Yao, B.; Brown, R.; Yue, G.H. Linear correlation between fractal dimension of EEG signal and handgrip force. Biol. Cybern. 2005, 93, 131–140. [Google Scholar] [CrossRef] [PubMed]
  16. Siemionow, V.; Yue, G.H.; Ranganathan, V.K.; Liu, J.Z.; Sahgal, V. Relationship between motor activity-related cortical po-tential and voluntary muscle activation. Exp. Brain Res. 2000, 133, 303–311. [Google Scholar] [CrossRef]
  17. Dettmers, C.; Fink, G.R.; Lemon, R.N.; Stephan, K.M.; Passingham, R.E.; Silbersweig, D.; Holmes, A.; Ridding, M.C.; Brooks, D.; Frackowiak, R. Relation between cerebral activity and force in the motor areas of the human brain. J. Neurophysiol. 1995, 74, 802–815. [Google Scholar] [CrossRef]
  18. Toyoshima, T.; Yazawa, S.; Murahara, T.; Ishiguro, M.; Shinozaki, J.; Ichihara-Takeda, S.; Shiraishi, H.; Matsuhashi, M.; Shimohama, S.; Nagamine, T. Load effect on background rhythms during motor execution: A magnetoencephalographic study. Neurosci. Res. 2016, 112, 26–36. [Google Scholar] [CrossRef]
  19. Boas, D.A.; Elwell, C.E.; Ferrari, M.; Taga, G. Twenty years of functional near-infrared spectroscopy: Introduction for the special issue. NeuroImage 2013, 85, 1–5. [Google Scholar] [CrossRef]
  20. Dans, P.; Foglia, S.; Nelson, A. Data Processing in Functional Near-Infrared Spectroscopy (fNIRS) Motor Control Research. Brain Sci. 2021, 11, 606. [Google Scholar] [CrossRef]
  21. Seidel, O.; Carius, D.; Roediger, J.; Rumpf, S.; Ragert, P. Changes in neurovascular coupling during cycling exercise measured by multi-distance fNIRS: A comparison between endurance athletes and physically active controls. Exp. Brain Res. 2019, 237, 2957–2972. [Google Scholar] [CrossRef]
  22. Baechle, T.R.; Earle, R.W. Essentials of Strength Training and Conditioning; Human Kinetics: Champaign, IL, USA, 2008. [Google Scholar]
  23. Oldfield, R.C. The assessment and analysis of handedness: The Edinburgh inventory. Neuropsychologia 1971, 9, 97–113. [Google Scholar] [CrossRef]
  24. Epstein, L.H.; Wing, R.R. Aerobic exercise and weight. Addict. Behav. 1980, 5, 371–388. [Google Scholar] [CrossRef]
  25. Petrofsky, J.S.; Phillips, C.A. The Use of Functional Electrical Stimulation for Rehabilitation of Spinal Cord Injured Patients. Central Nerv. Syst. Trauma 1984, 1, 57–74. [Google Scholar] [CrossRef] [PubMed]
  26. Delpy, D.T.; Cope, M.; Van Der Zee, P.; Arridge, S.; Wray, S.; Wyatt, J. Estimation of optical pathlength through tissue from direct time of flight measurement. Phys. Med. Biol. 1988, 33, 1433–1442. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Cheney, P.D. Role of Cerebral Cortex in Voluntary Movements: A Review. Phys. Ther. 1985, 65, 624–635. [Google Scholar] [CrossRef]
  28. Holper, L.; Biallas, M.; Wolf, M. Task complexity relates to activation of cortical motor areas during uni- and bimanual performance: A functional NIRS study. NeuroImage 2009, 46, 1105–1113. [Google Scholar] [CrossRef] [Green Version]
  29. Bajaj, S.; Drake, D.; Butler, A.J.; Dhamala, M. Oscillatory motor network activity during rest and movement: An fNIRS study. Front. Syst. Neurosci. 2014, 8, 13. [Google Scholar] [CrossRef] [Green Version]
  30. Homan, R.W. The 10–20 Electrode System and Cerebral Location. Am. J. EEG. Technol. 1988, 28, 269–279. [Google Scholar] [CrossRef]
  31. Singh, A.K.; Okamoto, M.; Dan, H.; Jurcak, V.; Dan, I. Spatial registration of multichannel multi-subject fNIRS data to MNI space without MRI. NeuroImage 2005, 27, 842–851. [Google Scholar] [CrossRef]
  32. Hou, X.; Zhang, Z.; Zhao, C.; Duan, L.; Gong, Y.; Li, Z.; Zhu, C. NIRS-KIT: A MATLAB toolbox for both resting-state and task fNIRS data analysis. Neurophotonics 2021, 8, 010802. [Google Scholar] [CrossRef]
  33. Jalalvandi, M.; Alam, N.R.; Sharini, H.; Hashemi, H.; Nadimi, M. Brain Cortical Activation during Imagining of the Wrist Movement Using Functional Near Infrared Spectroscopy (fNIRS). J. Biomed. Phys. Eng. 2021, 11, 583–594. [Google Scholar] [CrossRef] [PubMed]
  34. Friston, K.J.; Holmes, A.P.; Worsley, K.J.; Poline, J.-P.; Frith, C.D.; Frackowiak, R.S.J. Statistical Parametric Maps in Functional Imaging: A General Linear Approach. Hum. Brain Mapp. 1994, 2, 189–210. [Google Scholar] [CrossRef]
  35. Bates, D. Fitting Linear Mixed Models in R. R News 2005, 5, 27–30. [Google Scholar]
  36. Kuznetsova, A.; Brockhoff, P.B.; Christensen, R.H.B. lmerTest Package: Tests in linear mixed effects models. J. Stat. Softw. 2017, 82, 1–26. [Google Scholar] [CrossRef] [Green Version]
  37. Barton, K.; Barton, M.K. Package ‘Mumin’. Vienna, Austira. 2015, Volume 1, p. 439. Available online: https://www.r-project.org/ (accessed on 12 October 2022).
  38. Okamoto, M.; Tsuzuki, D.; Clowney, L.; Dan, H.; Singh, A.K.; Dan, I. Structural atlas-based spatial registration for functional near-infrared spectroscopy enabling inter-study data integration. Clin. Neurophysiol. 2009, 120, 1320–1328. [Google Scholar] [CrossRef]
  39. Abe, M.; Hanakawa, T. Functional coupling underlying motor and cognitive functions of the dorsal premotor cortex. Behav. Brain Res. 2009, 198, 13–23. [Google Scholar] [CrossRef]
  40. Rizzolatti, G.; Fogassi, L.; Gallese, V. Motor and cognitive functions of the ventral premotor cortex. Curr. Opin. Neurobiol. 2002, 12, 149–154. [Google Scholar] [CrossRef]
  41. Tanji, J. The supplementary motor area in the cerebral cortex. Neurosci. Res. 1994, 19, 251–268. [Google Scholar] [CrossRef]
  42. Wang, L.; Ma, L.; Yang, J.; Wu, J. Human Somatosensory Processing and Artificial Somatosensation. Cyborg Bionic Syst. 2021, 2021, 9843259. [Google Scholar] [CrossRef]
  43. Fujimoto, H.; Mihara, M.; Hattori, N.; Hatakenaka, M.; Kawano, T.; Yagura, H.; Miyai, I.; Mochizuki, H. Cortical changes underlying balance recovery in patients with hemiplegic stroke. NeuroImage 2014, 85, 547–554. [Google Scholar] [CrossRef]
  44. Mihara, M.; Miyai, I.; Hattori, N.; Hatakenaka, M.; Yagura, H.; Kawano, T.; Kubota, K. Cortical control of postural balance in patients with hemiplegic stroke. NeuroReport 2012, 23, 314–319. [Google Scholar] [CrossRef] [PubMed]
  45. Wriessnegger, S.C.; Kirchmeyr, D.; Bauernfeind, G.; Müller-Putz, G.R. Force related hemodynamic responses during execution and imagery of a hand grip task: A functional near infrared spectroscopy study. Brain Cogn. 2017, 117, 108–116. [Google Scholar] [CrossRef] [PubMed]
  46. Flansbjer, U.; Lexell, J.; Brogårdh, C. Long-term benefits of progressive resistance training in chronic stroke: A 4-year follow-up. J. Rehabil. Med. 2012, 44, 218–221. [Google Scholar] [CrossRef] [PubMed]
Figure 1. (A) Actual scenes of the experiment. The subject is seated on an armless chair, wearing loose clothes. (B) The high-density fNIRS probe. (C) The flow chart of experiments. Each subject passed through signal-quality screening (about 5 min), pre-scanning training (about 3 min), and task-state measurement (about 15 min).
Figure 1. (A) Actual scenes of the experiment. The subject is seated on an armless chair, wearing loose clothes. (B) The high-density fNIRS probe. (C) The flow chart of experiments. Each subject passed through signal-quality screening (about 5 min), pre-scanning training (about 3 min), and task-state measurement (about 15 min).
Brainsci 12 01539 g001
Figure 2. Three-dimensional location of each channel. We used NFRI fNIRS tools incorporated in NIRS-SPM to calculate the spatial registration of each channel from real space to MNI (Montreal Neurological Institute templates) coordinates and plotted them on the standard brain. We considered channels 6, 7, 8, 15, 16, 17, and 18 (red box) as the region of contralateral BA4 and channels 24, 25, 26, 33, 34, 35, and 36 (blue box) as the region of contralateral BA6 based on the probabilistic registration.
Figure 2. Three-dimensional location of each channel. We used NFRI fNIRS tools incorporated in NIRS-SPM to calculate the spatial registration of each channel from real space to MNI (Montreal Neurological Institute templates) coordinates and plotted them on the standard brain. We considered channels 6, 7, 8, 15, 16, 17, and 18 (red box) as the region of contralateral BA4 and channels 24, 25, 26, 33, 34, 35, and 36 (blue box) as the region of contralateral BA6 based on the probabilistic registration.
Brainsci 12 01539 g002
Figure 3. Channel-wise activation. Box-plot shows the median and percentile beta value of channels located in ROIs, and jitters showing the beta value of each subject. (A) The 0-pound condition. Five out of seven channels were activated in BA4 and two out of seven channels were activated in BA6. (B) The 3—pound condition. All channels in BA4 were activated, and six of seven channels were activated in BA6. (C) The 15—pound condition. The activated channels were the same as in the 3—pound condition but with a higher beta value. * FDR-adjusted p < 0.05 from one-sample t-test.
Figure 3. Channel-wise activation. Box-plot shows the median and percentile beta value of channels located in ROIs, and jitters showing the beta value of each subject. (A) The 0-pound condition. Five out of seven channels were activated in BA4 and two out of seven channels were activated in BA6. (B) The 3—pound condition. All channels in BA4 were activated, and six of seven channels were activated in BA6. (C) The 15—pound condition. The activated channels were the same as in the 3—pound condition but with a higher beta value. * FDR-adjusted p < 0.05 from one-sample t-test.
Brainsci 12 01539 g003
Figure 4. Averaged beta on a standardized brain model (n = 14). (A) The 0—pound condition, in which a typical contralateral activation pattern was observed. (B) The 3—pound condition, in which BA4 and BA6 show higher activation levels. (C) The 15—pound condition, in which BA4 and BA6 exhibited the highest activation among the three conditions.
Figure 4. Averaged beta on a standardized brain model (n = 14). (A) The 0—pound condition, in which a typical contralateral activation pattern was observed. (B) The 3—pound condition, in which BA4 and BA6 show higher activation levels. (C) The 15—pound condition, in which BA4 and BA6 exhibited the highest activation among the three conditions.
Brainsci 12 01539 g004
Figure 5. ROI-wise (Region of interest, ROI) activation. (A) BA4 was activated in all three conditions. (B) BA6 was activated in all three conditions. * p < 0.05, indicating beta-values were significantly greater than 0.
Figure 5. ROI-wise (Region of interest, ROI) activation. (A) BA4 was activated in all three conditions. (B) BA6 was activated in all three conditions. * p < 0.05, indicating beta-values were significantly greater than 0.
Brainsci 12 01539 g005
Figure 6. Representative examples of channel-wise correlation. (A) Channel 18, located in BA4. (B) Channel 25, located in BA6. (C) Channel 36, located in BA6. Representative channels were chosen based on the p-values of the mixed-effect model. We plotted the sample-regression line to display the overall tendency.
Figure 6. Representative examples of channel-wise correlation. (A) Channel 18, located in BA4. (B) Channel 25, located in BA6. (C) Channel 36, located in BA6. Representative channels were chosen based on the p-values of the mixed-effect model. We plotted the sample-regression line to display the overall tendency.
Brainsci 12 01539 g006
Figure 7. ROI-wise correlation between beta coefficients and weight. We plotted the sample-regression line to display the overall tendency. (A) Pooled and individual data in BA4. (B) Pooled and individual data in BA6. Effects of load on the cortical activation level (measured by the mean beta-value of ROIs) were fitted using mixed-effect linear models. The corresponding p-value and conditional r-square from mixed-effect linear models were placed above the regression line for each ROI. Significantly positive correlations were found in both ROIs.
Figure 7. ROI-wise correlation between beta coefficients and weight. We plotted the sample-regression line to display the overall tendency. (A) Pooled and individual data in BA4. (B) Pooled and individual data in BA6. Effects of load on the cortical activation level (measured by the mean beta-value of ROIs) were fitted using mixed-effect linear models. The corresponding p-value and conditional r-square from mixed-effect linear models were placed above the regression line for each ROI. Significantly positive correlations were found in both ROIs.
Brainsci 12 01539 g007
Table 1. Location of channels for subject 09.
Table 1. Location of channels for subject 09.
ChannelAnatomical Label in BAMNI CoordinateProbability
154—Primary Motor Cortex(−13.67, −31.33, 79.33)0.868
164—Primary Motor Cortex(−23.67, −31.33, 71)0.72
246—Pre-Motor and Supplementary Motor Cortex(−14.67, −14.33, 78)0.818
256—Pre-Motor and Supplementary Motor Cortex(−27.33, −16.33, 74.33)0.853
336—Pre-Motor and Supplementary Motor Cortex(−14.67, 0.33, 74)1
346—Pre-Motor and Supplementary Motor Cortex(−26.33, 0.33, 70.33)1
356—Pre-Motor and Supplementary Motor Cortex(−38.33, −1.67, 64.67)0.982
366—Pre-Motor and Supplementary Motor Cortex(−50.33, −2.67, 55.67)0.926
BA, Brodmann Area Atlas; MNI, Montreal Neurological Institute templates.
Table 2. Statistical details of channel-wise correlation.
Table 2. Statistical details of channel-wise correlation.
ChannelT-Value (df = 27)Slopep-ValueLocation
63.4640.0021290.0018BA4
73.6900.0021250.0010BA4
82.5760.0017970.0158BA4
152.9470.0009880.0065BA4
162.7670.0009950.0101BA4
183.7780.0017240.0008BA4
244.1810.0015050.0003BA6
254.8620.0017840.0000BA6
264.0420.0012240.0004BA6
333.7470.0017390.0009BA6
342.8560.0010070.0082BA6
354.5340.0012600.0001BA6
364.4260.0015850.0001BA6
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chen, Z.; Song, X.; Qiao, Y.; Yan, J.; Zhu, C.; Xie, Q.; Niu, C.M. Increased Inertia Triggers Linear Responses in Motor Cortices during Large-Extent Movements—A fNIRS Study. Brain Sci. 2022, 12, 1539. https://doi.org/10.3390/brainsci12111539

AMA Style

Chen Z, Song X, Qiao Y, Yan J, Zhu C, Xie Q, Niu CM. Increased Inertia Triggers Linear Responses in Motor Cortices during Large-Extent Movements—A fNIRS Study. Brain Sciences. 2022; 12(11):1539. https://doi.org/10.3390/brainsci12111539

Chicago/Turabian Style

Chen, Zhi, Xiaohui Song, Yongjun Qiao, Jin Yan, Chaozhe Zhu, Qing Xie, and Chuanxin M. Niu. 2022. "Increased Inertia Triggers Linear Responses in Motor Cortices during Large-Extent Movements—A fNIRS Study" Brain Sciences 12, no. 11: 1539. https://doi.org/10.3390/brainsci12111539

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop