Skip to main content

Quantified assessment of deep brain stimulation on Parkinson’s patients with task fNIRS measurements and functional connectivity analysis: a pilot study



Deep brain stimulation (DBS) has proved effective for Parkinson’s disease (PD), but the identification of stimulation parameters relies on doctors’ subjective judgment on patient behavior.


Five PD patients performed 10-meter walking tasks under different brain stimulation frequencies. During walking tests, a wearable functional near-infrared spectroscopy (fNIRS) system was used to measure the concentration change of oxygenated hemoglobin (HbO2) in prefrontal cortex, parietal lobe and occipital lobe. Brain functional connectivity and global efficiency were calculated to quantify the brain activities.


We discovered that both the global and regional brain efficiency of all patients varied with stimulation parameters, and the DBS pattern enabling the highest brain efficiency was optimal for each patient, in accordance with the clinical assessments and DBS treatment decision made by the doctors.


Task fNIRS assessments and brain functional connectivity analysis promise a quantified and objective solution for patient-specific optimization of DBS treatment.

Trial registration

Name: Accurate treatment under the multidisciplinary cooperative diagnosis and treatment model of Parkinson’s disease. Registration number is ChiCTR1900022715. Date of registration is April 23, 2019.


Parkinson’s disease (PD) is a neurodegenerative disease caused by the progressive loss of nigrostriatal dopaminergic neurons in substantia nigra pars compacta [1]. The loss of dopaminergic neurons induces severe motor symptoms such as tremor, rigidity, bradykinesia and dyskinesia, as well as non-motor symptoms such as constipation, fatigue, anxiety, cognitive dysfunction, and dementia [2, 3]. Deep brain stimulation (DBS) has proved an effective therapy for symptom improvement after PD, especially in the late stage when medication is less effective [4, 5].

With DBS operation, electrodes are implanted into specific target brain locations, such as the subthalamic nucleus or globus pallidus internal, as shown in Fig. 1. The electrodes work as the neurostimulator and send electrical stimulation for the treatment of movement disorders [6, 7]. DBS programming, the identification of stimulation parameters of the implanted neurostimulator for symptom management, is crucial for successful and optimal treatment. However, the functioning neural mechanism of DBS remains unclear, and DBS programming in current clinical practice is typically conducted by doctors according to their observation of patient behavior, strongly relying on the doctors’ skill, experience, and subjective judgment [810]. An objective approach that uses personalized neurophysiological measurements to optimize DBS programming is therefore highly demanded.

Fig. 1
figure 1

Deep brain stimulation: neurologists implant a neurostimulator to send electrical stimulation through the implanted electrodes to specific targets. After neurosurgery, DBS programming was conducted to optimize the stimulation parameters for the patient

As a neurodegenerative disease, PD damages the central nervous system and its function, leading to movement disorders. While DBS improves motor symptoms, changes with respect to brain function definitely happen. Brain functional connectivity (BFC), which refers to the statistical correlation between physiological signals from PD-related brain regions, might characterize DBS-induced functional variation and provide a quantified and objective measure for DBS treatment. Different brain regions communicate and coordinate during function fulfillment and task execution, and brain functional connectivity represents this effort [1113]. Algorithms have been developed for BFC analysis of neurological disorders [14, 15]. With the DBS treatment, the BFC strength indicates the “cost” of the brain while the patient trying to complete a specific task. The lower the cost, the higher communication efficiency among the brain regions, and the better the DBS parameters might be. In this paper, we verify this idea by clinical tests of PD patients.

To investigate the brain activation of PD patients, previous studies have applied SPECT [16, 17], PET[18], fMRI[19, 20], and EEG [21, 22], mainly by comparison between patients and the healthy controls for diagnosis, but none for DBS programming. Moreover, none of fMRI, SPECT, or PET could allow “task-state” measurements while the patient trying to complete a motor task, and the preparation process for EEG measurement is tedious and challenging for PD patients after DBS surgery, restricting the possibilities for these imaging modalities to be applied in DBS programming. Instead, functional near-infrared spectroscopy (fNIRS) is an optical functional neuroimaging technique that uses near-infrared light (700–900 nm spectral interval) to perform continuous and non-invasive monitoring of blood hemoglobin changes related to brain functions [23]. The fNIRS-based functional connectivity has been successfully used in clinical applications. In [24], Qitao et al. assessed resting-state functional connectivity in cerebral infarction patients with fNIRS. In [25], Didem et al. used functional connectivity features to perform clinical binary classification in patients with fibromyalgia. A further advantage for our study is that the fNIRS equipment can be made portable, enabling “task-state” measurements during clinical evaluation such as walking.

With a wearable fNIRS system, we measured five PD patients in the clinical DBS programming process, recording their brain activities during standard walking tests under different brain stimulation frequencies that were specified by the doctors. The fNIRS-based brain functional connectivity analysis was conducted on three brain regions: the prefrontal cortex, parietal lobe, and occipital lobe. The coordination and communication among these brain regions are essentially evolved during the walking test in clinical DBS programming and their functional connectivity can characterize DBS-induced improvement of neural function.



The patients of consideration were clinically diagnosed idiopathic Parkinson’s disease without other plus syndromes. Exclusion criteria for patient recruitment include the following: (1) unnecessary to re-adjust the stimulation parameters; (2) unable to stand or walk for 90 s at a time; (3) any factors affecting their gait performance, such as idiopathic scoliosis and leg injury; (4) any mental diseases, such as neuropsychiatric comorbidity, schizophrenia, and personality disorders; and (5) age > 70 years. Five PD patients that received DBS surgery and qualified for the study were recruited. These PD patients had bilateral symptoms and were treated with bilateral stimulation. Each patient was fully informed of the experimental purpose and procedures and provided written consent prior to the measurement. The clinical characteristics and initial stimulation parameters in the DBS surgery of these patients are shown in Table 1.

Table 1 Clinical characteristics of the five PD patients for DBS programming

Functional near-infrared spectroscopy

A wearable, wireless, continuous-wave fNIRS system (Nirsmart, Danyang Huichuang Medical Equipment Co, Ltd, China) [26] was used to monitor the concentration change of oxygenated hemoglobin (HbO2). The wavelength of the near-infrared light was 760 nm, and the sampling rate was 10 Hz. Six regions, i.e., left and right prefrontal cortex (L/R-PFC), parietal lobe (L/R-PL), and occipital lobe (L/R-OL) were chosen as the areas of interest for recording, and 34 fNIRS electrodes including 16 sources and 18 detectors were placed to the selected region, as shown in Fig. 2. The prefrontal cortex is implicated in cognitive control and information processing for complex behavior. The parietal lobe plays an important role in motor function, working memory, and the integration of multiple sensory information. The occipital lobe is responsible for visual processing, working memory, and modulation of different sensory stimulation.

Fig. 2
figure 2

Deployment of the fNIRS sources (16 orange rectangles) and detectors (18 blue dots). Ci indicates the ith channel. The six brain regions of interest, i.e., L/R-PFC, L/R-PL, and L/R-OL, are separated by green frames

Experimental design

The crucial brain stimulation parameters in DBS programming consist of the location of electrode contact, voltage amplitude, and frequency. In our experiments, the DBS frequency was varied while the location of electrode contact and voltage amplitude were invariant. The reasons of only varying DBS frequency included the following: (1) the locations of electrode contact had been optimized by doctors with MRI scans and 3D reconstruction technique. The voltage amplitudes were fixed by the doctors with clinical diagnosis. (2) DBS frequency was related to the improvement of gait and balance [2730]. After varying the DBS frequency, the patients performed the walking test, as shown in Fig. 3. It has to be noted that DBS patients are normally not strong enough to take too many walking tests. For the five PD patients of this study, we limited the number of walking tests to be 4 at most. The process of the walking test was as follows: (1) the doctor performed the frequency adjustment, (2) the patient sat on a chair for 5 min to ensure that the new DBS paradigm actually took effect, (3) the patient stood up from the chair and stood quite for 30 s, (4) the patient performed the 10-m walking task, and (5) the patient stood quite for 30 s again. Instructions of “Standing,” “Walking,” “Stop,” and “Finish” were given by the doctors during each test. PD patients were tested under medicine off condition.

Fig. 3
figure 3

The experimental procedure. Left: each DBS patient performed n test sessions, n is the number of frequency adjustment. Each session contains frequency adjustment and walking test. The walking test involves 30 s of quiet standing at the beginning and end, intermediated with 30 s of consecutive walking. Instructions about the walking test ("Standing for 30 s," "Walking," "Stop," and "Finish") were given by the doctor. Right: the patients wore a portable fNIRS measurement equipment

Data processing


Firstly, the collected fNIRS measurement data were processed with a 0.01–0.2 bandpass filter to remove the instrumental and physiological noises (e.g., heartbeats, respirations and Mayer waves) [3133]. Then, the HbO2 of each channel was calculated with the filtered data according to the modified Beer-Lambert law [34]. Further, motion artifacts were removed based on moving standard deviation and spline interpolation [35].

Functional connectivity and global efficiency

Firstly, the Pearson’s correlation coefficient \(\phantom {\dot {i}\!}P_{C_{x}C_{y}}\) between two channels was calculated as follows:

$$ \begin{aligned} P_{C_{x}C_{y}}&=\frac{cov(C_{x},C_{y})}{\sigma_{C_{x}}\sigma_{C_{y}}}\\ &=\frac{\Sigma_{i=1}^{m}(C_{x,i}-\bar{C_{x}})(C_{y,i}-\bar{C_{y}})}{\sqrt{\Sigma_{i=1}^{m}(C_{x,i}-\bar{C_{x}})^{2}}\sqrt{\Sigma_{i=1}^{m}(C_{y,i}-\bar{C_{y}})^{2}}} \end{aligned} $$

where Cx and Cy are the measurements (HbO2) of the xth and yth channels, cov(Cx,Cy) is the covariance between Cx and \(\phantom {\dot {i}\!}C_{y}, \sigma _{C_{x}}\) and \(\phantom {\dot {i}\!}\sigma _{C_{y}}\) are the standard deviation of Cx and Cy, m is the length of the measurement, Cx,i and Cy,i are the ith measurement (HbO2) of the xth and yth channels, and \(\bar {C_{x}}=\frac {1}{m}\sum _{i=1}^{m}C_{x,i}\) and \(\bar {C_{y}}=\frac {1}{m}\sum _{i=1}^{m}C_{y,i}\) are the average values of Cx and Cy.

Then, Fisher’s z-transformation was applied to decrease the skewness of \(\phantom {\dot {i}\!}P_{C_{x}C_{y}}\) and normalize its distribution:

$$ F_{C_{x}C_{y}} = artanh(P_{C_{x}C_{y}})=\frac{1}{2}ln\left(\frac{1+P_{C_{x}C_{y}}}{1-P_{C_{x}C_{y}}}\right) $$

where \(\phantom {\dot {i}\!}F_{C_{x}C_{y}}\) is the connection strength between channel Cx and Cy, and artanh(·) is the inverse hyperbolic tangent function. The \(\phantom {\dot {i}\!}F_{C_{x}C_{y}}\) values were used to construct the connectivity matrix M, which is defined as follows:

$$ M = \left[\begin{array}{cccc} F_{C_{N}C_{1}} & F_{C_{N}C_{2}}& \cdots & F_{C_{N}C_{N}}\\ \vdots & \vdots & \ddots & \vdots \\ F_{C_{2}C_{1}} & F_{C_{2}C_{2}} & \cdots & F_{C_{2}C_{N}}\\ F_{C_{1}C_{1}} & F_{C_{1}C_{2}} & \cdots & F_{C_{1}C_{N}} \end{array}\right] $$

where N is the channel number of global brain regions. It should be noted that higher connection strength in the connectivity matrix M corresponds to lower brain communication efficiency [3638].

Further, the global efficiency (GE) describes the overall communication efficiency:

$$ GE = \frac{1}{\frac{1}{N} \sum_{x=1}^{N} \frac{\sum_{y=1,y\neq x}^{N}F_{C_{x}C_{y}}}{N-1}} $$

where N is the channel number of global brain regions. Cx and Cy indicate the measurement (HbO2) of the xth and yth channels. \(\phantom {\dot {i}\!}F_{C_{x}C_{y}}\) indicates the connection strength between channel Cx and Cy. Higher GE scores represent lower communication strength of global regions.

Local strength

Local strength represents the communication strength between different brain regions. The local strength of PFC, PL, and OL are respectively defined as LSPFC,LSPL, and LSOL, i.e.,

$$ \begin{aligned} {LS}_{PFC} &= \frac{1}{2}\left(\frac{1}{N_{LPFC}\times N} \sum_{x=1}^{N_{LPFC}} \sum_{y=1,y\neq x}^{N} F_{C_{x}C_{y}} \right.\\ &+\left. \frac{1}{N_{RPFC}\times N} \sum_{x=1}^{N_{RPFC}} \sum_{y=1,y\neq x}^{N} F_{C_{x}C_{y}}\right) \end{aligned} $$
$$ \begin{aligned} {LS}_{PL} &= \frac{1}{2}\left(\frac{1}{N_{LPL}\times N} \sum_{x=1}^{N_{LPL}} \sum_{y=1,y\neq x}^{N} F_{C_{x}C_{y}} \right.\\ &+\left. \frac{1}{N_{RPL}\times N} \sum_{x=1}^{N_{RPL}} \sum_{y=1,y\neq x}^{N}F_{C_{x}C_{y}}\right) \end{aligned} $$
$$ \begin{aligned} {LS}_{OL} &= \frac{1}{2}\left(\frac{1}{N_{LOL}\times N} \sum_{x=1}^{N_{LOL}} \sum_{y=1,y\neq x}^{N} F_{C_{x}C_{y}}\right. \\ &+\left. \frac{1}{N_{ROL}\times N} \sum_{x=1}^{N_{ROL}} \sum_{y=1,y\neq x}^{N}F_{C_{x}C_{y}}\right) \end{aligned} $$

where N is the channel number of global brain regions. NLPFC,NRPFC,NLPL,NRPL,NLOL, and NROL indicate the channel number of LPFC, RPFC, LPL, RPL, LOL, and ROL, respectively. Moreover, the averaged local strength (Laver) is defined as follows:

$$ L_{aver} = \frac{N_{PFC}}{N}\times L_{PFC}^{\prime} + \frac{N_{PL}}{N}\times L_{PL}^{\prime} + \frac{N_{OL}}{N}\times L_{OL}^{\prime} $$

where NPFC,NPL and NOL indicate the channel number in PFC, PL, and OL. \(\phantom {\dot {i}\!}L_{PFC}^{\prime }, L_{PL}^{\prime }\), and \(\phantom {\dot {i}\!}L_{OL}^{\prime }\) represent the normalized LPFC,LPL, and LOL, respectively. \(\phantom {\dot {i}\!}L_{PFC}^{\prime }, L_{PL}^{\prime }\), and \(\phantom {\dot {i}\!}L_{OL}^{\prime }\) are defined as follows:

$$ \small L_{PFC}^{'} = \frac{L_{PFC}-L_{PFC}^{min}}{L_{PFC}^{m}ax-L_{PFC}^{min}} $$
$$ \small L_{PL}^{'} = \frac{L_{PL}-L_{PL}^{min}}{L_{PL}^{m}ax-L_{PL}^{min}} $$
$$ \small L_{OL}^{'} = \frac{L_{OL}-L_{OL}^{min}}{L_{OL}^{m}ax-L_{OL}^{min}} $$

where \(L_{PFC}^{min}, L_{PL}^{min}\) and \(L_{OL}^{min}\) indicate the minimum values of LPFC,LPL and LOL of each PD patient. \(L_{PFC}^{max}, L_{PL}^{max}\) and \(L_{OL}^{max}\) indicate the maximum values of LPFC,LPL, and LOL of each PD patient.


Brain functional connectivity and global efficiency

The measurements of fNIRS during the 10-m walking of each patient were recorded, and the brain connectivity matrix (CM) as well as the global efficiency (GE) were calculated with respect to each tested frequency.

All results of the 5 patients are shown in Fig. 4. Each of the square color frame presented the CM item at the specific frequency labeled below (the corresponding GE value was also listed following the frequency value). The x and y axes of the square frame were the channel number of the fNIRS, and the color at (x,y) indicated the value of the CM, which represented the connectivity strength between Channel-x and Channel-y. The higher the connectivity strength (the warmer the color), the lower the communication efficiency. Moreover, the global efficiency (GE) is defined to describe the overall communication efficiency. And the frequency corresponding to the coldest color distribution, i.e., the largest GE value, is promisingly optimal for each patient, respectively.

Fig. 4
figure 4

The brain functional connectivity matrices calculated from the fNIRS measurements of the five PD patients (P1, P2, P3, P4, P5) at different DBS frequencies. The color bar corresponding to the values of CM was presented at the bottom. The frequency values in red are promisingly optimal for each patient

In addition to the global connectivity, local connection strengths of PFC, PL, and OL, i.e., LSPFC,LSPL, and LSOL, as well as the averaged local strength LSavg were also calculated, as shown in Table 2. It can be seen that for each patient, the frequency at which the smallest average local connectivity was obtained was exactly the one where the highest GE was achieved. The results on local connectivity and global efficiency were consistent. Therefore, the frequencies of 125Hz, 130Hz, 100Hz, 160Hz, and 130Hz are promisingly optimal for patients P1, P2, P3, P4, and P5, respectively.

Table 2 The local functional connectivity analysis

Comparison with clinical assessments

The proposed quantified assessment based on task fNIRS measurements and brain functional analysis was compared with clinical assessments including the on-site DBS programming decisions and post independent MDS-UPDRS ratings on the recorded videos of walking tests. The DBS programming decision, MDS-UPDRS ratings, and analysis on brain functional efficiency were conducted independently from each other, and the relative results were not revealed to the personnel performing other analysis until the entire study was accomplished. The results are shown in Table 3.

Table 3 Comparison of brain efficiency with clinical assessments including the on-site DBS programming decisions and post independent MDS-UPDRS ratings based on the recorded videos

For all the five patients, the doctors’ on-site decision on DBS frequency were consistent with the optimal frequency by the brain functional connectivity analysis, i.e., the one that corresponds to both the highest global efficiency and lowest local connection strength. The doctors made the decisions according to the following clinical process [39, 40]: (1) assessing the rigidity and tremor of PD patients; (2) measuring the gait performance with MDS-UPDRS ratings; (3) communicating with patients about their comprehensive feeling such as comfortability, etc; and (4) choosing the best DBS frequency considering the rigidity, tremor, gait, and patients’ feeling. In the experiments, the doctors’ on-site decision on DBS frequency was consistent with frequencies that patients felt the most comfortable. One-month follow-up after DBS programming reported that all patients were satisfied with the DBS treatment, and no adverse effect or feelings was reported or noticed by the patients, their families, or doctors. The dopaminergic medication was unchanged during the one-month follow-up period.

MDS-UPDRS, the Unified Parkinson’s Disease Rating Scale by the Movement Disorder Society, has been universally used in clinical assessment on motor and non-motor aspects of Parkinson’s disease [41]. The videos of the five patients while performing the walking tests were sent to two qualified specialists, S1 and S2, who are experienced on gait assessment of PD patients and independent from this study. The MDS-UPDRS scores include 5 subratings, i.e., 0, 1, 2, 3, and 4, indicating normal, slight, mild, moderate, and severe symptoms, respectively. For all the five patients, both specialists gave the same MDS-UPDRS scores of 1, indicating slight movement disorder, for the patients’ walking performance under the applied DBS treatments that were also consistent with the best brain efficiency. This verified the efficacy of the DBS treatments and brain functional analysis. Nevertheless, the MDS-UPDRS scores were unable to further discriminate less significant differences in motor performance. For all the tested frequencies of patients P2, P4, and P5, both specialists gave the same MDS-UPDRS scores. For the performance differences that could be distinguished by the MDS-UPDRS scores, as patients P1 and P3, the results were all consistent with the brain functional analysis.


Starting from the fact that Parkinson’s disease is a neurological disorder and implicates motor performance, we proposed in this paper an objective and quantitative assessment method for DBS programming with task fNIRS measurements and brain functional connectivity analysis. To the best of our knowledge, this is the first fNIRS-based study on assessment and optimization of DBS therapy via recording of brain activation while performing motor tasks and analysis on global as well as regional brain efficiency. The methods were developed for post-operative DBS programming, but also has the potential for intra-operative assessment on correct targeting.

Brain functional analysis has been a long-time research focus of Parkinson’s disease. In [18], Rascol et al. measured the regional cerebral blood flow changes with PET during the execution of a finger-to-thumb opposition motor task in the cerebellar hemisphere of parkinsonian patients. Compared with healthy controls, Parkinson’s patients had increased brain activation in ipsilateral cerebellar hemisphere. In [42], Sabatini et al. analyzed the cortical change of PD patients in a complex sequential motor task with fMRI. Compared with normal controls, PD patients had a significant bilateral increase of fMRI signals in the primary sensorimotor cortex, lateral premotor cortex, inferior parietal cortex, caudal part of SMA, and anterior cingulate cortex. In [19], Zhang et al. analyzed the functional connectivity of ventral intermediate nucleus of thalamus (Vim) in tremor-dominant (TD) and akinetic-/rigid-dominant (ARD) PD patients with fMRI. In TD patients, the Vim nucleus had an increase of brain connectivity with dentate nucleus, primary motor cortex (M1), SMA, globus pallidus, premotor cortex, and parietal cortex compared with normal controls. In ARD patients, the Vim nucleus only exhibited increased connectivity with globus pallidus and limbic lobe compared with normal controls. In [20], Hou et al. evaluated the functional connectivity of default mode network (DMN) with resting-state fMRI data and found significantly increased connectivities of anterior DMN and prefrontal regions. The common finding of these studies is that PD patients are characterized with hypoactivation of SMA and hyperactivation of cortical motor regions (e.g., primary motor cortex, premotor cortex, parietal cortex) compared with normal controls. DBS can relatively normalize the hypoactivation of SMA and hyperactivation of other cortical regions and optimize the network profile toward healthy controls [43, 44]. Our findings are consistent with these studies. Effective DBS parameters could induce strong normalization and decrease the “extra” brain activation of PD patients in order to optimize the network profile toward healthy controls.

Gait performance is a primary concern for the patients and doctors. Therefore, we took the clinical 10-m walking test as the motor task for assessment, and accordingly, the DBS frequency as the varying parameter since it is directly associated with gait [28]. Other DBS parameters, the location of electrode contact, voltage amplitude, and impulse duration, as well as other motor and non-motor functions, are also important and can be addressed in future studies. The type and number of motor tasks are constrained by the patients’ physical condition.

Besides fNIRS, EEG can also measure brain activation for functional analysis of Parkinson’s disease, but typically for non-motor functions such as cognition [45] and emotion [46]. Technical challenges for EEG measurement during motor tasks include motion artifact and noise removal, source localization, fast set-up, etc. Although fMRI prohibits entry of DBS patients due to the electromagnetic fields, its excellent localization accuracy can facilitate brain functional analysis of PD patients without DBS for non-motor and motor tasks, as it has done to investigate rehabilitation induced brain reorganization after stroke with upper and lower extremity movements assisted by special mechatronic systems [47, 48].

A major limitation of this study is the small number of patients. Nevertheless, our experiment design and analysis method were based on the current understanding of the neurological mechanism of Parkinson’s disease, and the results were in line with this knowledge and also consistent across all the five patients. Moreover, our fNIRS system could calculate the results in 90.47 s, which is effective to assist the doctors in selecting the stimulus parameters. We hope this work can encourage more study, and more clinical evidence will promisingly enable quantified and individualized optimization of deep brain stimulation therapy for each Parkinson’s patient.

In the future, we will try to expand the number of patients and conduct a more detail investigation of DBS treatment optimization with different contacts, voltage, and pulse width. And we will try to explore more objective evaluation indicators and find appropriate gait analysis equipment for our experiments. Moreover, we will consider the change of total electrical energy delivered (TEED) [49, 50]. TEED is a comprehensive parameter considering frequency, voltage, pulse width, and impedance, which is directly related to the power consumption and battery drainage rate of implanted impulse generator.


This was a pilot study on quantified assessment of DBS programming. For the first time, we recorded the brain signals of PD patients in clinical DBS programming process with a wearable fNIRS system and analyzed the collected signals for brain functional connectivity. Experimental results showed that fNIRS assessments and brain functional connectivity analysis promised an objective solution for patient-specific optimization of DBS treatment.

Availability of data and materials

Data that support the findings and software codes developed for the data analysis in this paper will be made available upon reasonable request to the corresponding authors.



Deep brain stimulation


Parkinson’s disease


Functional near-infrared spectroscopy


Prefrontal cortex


Parietal lobe


Occipital lobe

GE :

Global efficiency

CM :

Connectivity matrix


  1. Lang AE, Lozano AM. Parkinson’s disease. N Engl J Med. 1998; 339(16):1130–43.

    Article  CAS  PubMed  Google Scholar 

  2. Jankovic J. Parkinson’s disease: clinical features and diagnosis. J Neurol Neurosurg Psychiatry. 2008; 79(4):368–76.

    Article  CAS  PubMed  Google Scholar 

  3. Zis P, Erro R, Walton CC, Sauerbier A, Chaudhuri KR. The range and nature of non-motor symptoms in drug-naive Parkinson’s disease patients: a state-of-the-art systematic review. NPJ Parkinson’s Dis. 2015; 1(1):1–8.

    Google Scholar 

  4. Bouthour W, Mégevand P, Donoghue J, Lüscher C, Birbaumer N, Krack P. Biomarkers for closed-loop deep brain stimulation in parkinson disease and beyond. Nat Rev Neurol. 2019; 15(6):343–52.

    Article  PubMed  Google Scholar 

  5. Lozano AM, Lipsman N, Bergman H, Brown P, Chabardes S, Chang JW, Matthews K, McIntyre CC, Schlaepfer TE, Schulder M, et al. Deep brain stimulation: current challenges and future directions. Nat Rev Neurol. 2019; 15(3):148–60.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Benabid AL. Deep brain stimulation for Parkinson’s disease. Curr Opin Neurobiol. 2003; 13(6):696–706.

    Article  CAS  PubMed  Google Scholar 

  7. Perlmutter JS, Mink JW. Deep brain stimulation. Ann Rev Neurosci. 2006; 29:229–57.

    Article  CAS  PubMed  Google Scholar 

  8. Groiss S, Wojtecki L, Südmeyer M, Schnitzler A. Deep brain stimulation in Parkinson’s disease. Ther Adv Neurol Disord. 2009; 2(6):379–91.

    Article  Google Scholar 

  9. Chen S, Gao G, Feng T, Zhang J. Chinese expert consensus on programming deep brain stimulation for patients with Parkinson’s disease. Transl Neurodegener. 2018; 7(1):11.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Koeglsperger T, Palleis C, Hell F, Mehrkens JH, Boetzel K. Deep brain stimulation programming for movement disorders: current concepts and evidence-based strategies. Front Neurol. 2019; 10:410.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Caspi A, Moffitt TE. Gene–environment interactions in psychiatry: joining forces with neuroscience. Nat Rev Neurosci. 2006; 7(7):583–90.

    Article  CAS  PubMed  Google Scholar 

  12. Van Den Heuvel MP, Pol HEH. Exploring the brain network: a review on resting-state fMRI functional connectivity. Eur Neuropsychopharmacol. 2010; 20(8):519–34.

    Article  CAS  PubMed  Google Scholar 

  13. Friston KJ. Functional and effective connectivity: a review. Brain Connectivity. 2011; 1(1):13–36.

    Article  PubMed  Google Scholar 

  14. Lynall M-E, Bassett DS, Kerwin R, McKenna PJ, Kitzbichler M, Muller U, Bullmore E. Functional connectivity and brain networks in schizophrenia. J Neurosci. 2010; 30(28):9477–87.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Niethammer M, Tang CC, Vo A, Nguyen N, Spetsieris P, Dhawan V, Ma Y, Small M, Feigin A, During MJ, et al. Gene therapy reduces parkinson’s disease symptoms by reorganizing functional brain connectivity. Sci Transl Med. 2018; 10(469):0713.

    Article  CAS  Google Scholar 

  16. Benamer HT, Patterson J, Wyper DJ, Hadley DM, Macphee G, Grosset DG. Correlation of Parkinson’s disease severity and duration with 123I-FP-CIT SPECT striatal uptake. Mov Disord. 2000; 15(4):692–8.

    Article  CAS  PubMed  Google Scholar 

  17. Sommer U, Hummel T, Cormann K, Mueller A, Frasnelli J, Kropp J, Reichmann H. Detection of presymptomatic Parkinson’s disease: combining smell tests, transcranial sonography, and SPECT. Mov Disord. 2004; 19(10):1196–202.

    Article  PubMed  Google Scholar 

  18. Rascol O, Sabatini U, Fabre N, Brefel C, Loubinoux I, Celsis P, Senard J, Montastruc J, Chollet F. The ipsilateral cerebellar hemisphere is overactive during hand movements in akinetic parkinsonian patients. Brain. 1997; 120(1):103–10.

    Article  PubMed  Google Scholar 

  19. Zhang J-R, Feng T, Hou Y-N, Chan P, Wu T. Functional connectivity of vim nucleus in tremor-and akinetic-/rigid-dominant Parkinson’s disease. CNS Neurosci Ther. 2016; 22(5):378–86.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Hou Y, Luo C, Yang J, Ou R, Liu W, Song W, Gong Q, Shang H. Default-mode network connectivity in cognitively unimpaired drug-naive patients with rigidity-dominant Parkinson’s disease. J Neurol. 2017; 264(1):152–60.

    Article  CAS  PubMed  Google Scholar 

  21. Han C-X, Wang J, Yi G-S, Che Y-Q. Investigation of EEG abnormalities in the early stage of Parkinson’s disease. Cogn Neurodyn. 2013; 7(4):351–9.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Markser A, Maier F, Lewis C, Dembek T, Pedrosa D, Eggers C, Timmermann L, Kalbe E, Fink G, Burghaus L. Deep brain stimulation and cognitive decline in Parkinson’s disease: the predictive value of electroencephalography. J Neurol. 2015; 262(10):2275–84.

    Article  CAS  PubMed  Google Scholar 

  23. Scholkmann F, Kleiser S, Metz AJ, Zimmermann R, Pavia JM, Wolf U, Wolf M. A review on continuous wave functional near-infrared spectroscopy and imaging instrumentation and methodology. Neuroimage. 2014; 85:6–27.

    Article  PubMed  Google Scholar 

  24. Tan Q, Zhang M, Wang Y, Zhang M, Wang Y, Xin Q, Wang B, Li Z. Frequency-specific functional connectivity revealed by wavelet-based coherence analysis in elderly subjects with cerebral infarction using NIRS method. Med Phys. 2015; 42(9):5391–403.

    Article  PubMed  Google Scholar 

  25. Gokcay D, Eken A, Baltaci S. Binary classification using neural and clinical features: an application in fibromyalgia with likelihood-based decision level fusion. IEEE J Biomed Health Informatics. 2019; 23(4):1490–8.

    Article  Google Scholar 

  26. Wang Z, Zhou Y, Chen L, Gu B, Yi W, Liu S, Xu M, Qi H, He F, Ming D. BCI monitor enhances electroencephalographic and cerebral hemodynamic activations during motor training. IEEE Trans Neural Syst Rehabil Eng. 2019; 27(4):780–7.

    Article  PubMed  Google Scholar 

  27. Jia F, Guo Y, Wan S, Chen H, Hao H, Zhang J, Li L. Variable frequency stimulation of subthalamic nucleus for freezing of gait in Parkinson’s disease. Parkinsonism Relat Disord. 2015; 21(12):1471–2.

    Article  PubMed  Google Scholar 

  28. Xie T, Vigil J, MacCracken E, Gasparaitis A, Young J, Kang W, Bernard J, Warnke P, Kang UJ. Low-frequency stimulation of STN-DBS reduces aspiration and freezing of gait in patients with PD. Neurology. 2015; 84(4):415–20.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Zibetti M, Moro E, Krishna V, Sammartino F, Picillo M, Munhoz RP, Lozano AM, Fasano A. Low-frequency subthalamic stimulation in Parkinson’s disease: long-term outcome and predictors. Brain Stimulation. 2016; 9(5):774–9.

    Article  PubMed  Google Scholar 

  30. Jia F, Hu W, Zhang J, Shukla AW, Almeida L, Meng F. -g., Okun MS, Li L. Variable frequency stimulation of subthalamic nucleus in Parkinson’s disease: rationale and hypothesis. Parkinsonism Relat Disord. 2017; 39:27–30.

    Article  PubMed  Google Scholar 

  31. Hu X, Zhuang C, Wang F, Liu Y-J, Im C-H, Zhang D. fNIRS evidence for recognizably different positive emotions. Front Hum Neurosci. 2019; 13.

  32. Klein F, Kranczioch C. Signal processing in fNIRS: a case for the removal of systemic activity for single trial data. Front Hum Neurosci. 2019; 13:331.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Huo C, Li X, Jing J, Ma Y, Li W, Wang Y, Liu W, Fan Y, Yue S, Wang Y, et al. Median nerve electrical stimulation–induced changes in effective connectivity in patients with stroke as assessed with functional near-infrared spectroscopy. Neurorehabil Neural Repair. 2019; 33(12):1008–17.

    Article  PubMed  Google Scholar 

  34. Cope M, Delpy DT. System for long-term measurement of cerebral blood and tissue oxygenation on newborn infants by near infra-red transillumination. Med Biol Eng Comput. 1988; 26(3):289–294.

    Article  CAS  PubMed  Google Scholar 

  35. Scholkmann F, Spichtig S, Muehlemann T, Wolf M. How to detect and reduce movement artifacts in near-infrared imaging using moving standard deviation and spline interpolation. Physiol Meas. 2010; 31(5):649.

    Article  CAS  PubMed  Google Scholar 

  36. Latora V, Marchiori M. Efficient behavior of small-world networks. Phys Rev Lett. 2001; 87(19):198701.

    Article  CAS  PubMed  Google Scholar 

  37. Rubinov M, Sporns O. Complex network measures of brain connectivity: uses and interpretations. Neuroimage. 2010; 52(3):1059–69.

    Article  PubMed  Google Scholar 

  38. Bullmore E, Sporns O. Complex brain networks: graph theoretical analysis of structural and functional systems. Nat Rev Neurosci. 2009; 10(3):186–98.

    Article  CAS  PubMed  Google Scholar 

  39. Bronstein JM, Tagliati M, Alterman RL, Lozano AM, Volkmann J, Stefani A, Horak FB, Okun MS, Foote KD, Krack P, et al. Deep brain stimulation for Parkinson disease: an expert consensus and review of key issues. Arch Neurol. 2011; 68(2):165.

    Article  PubMed  Google Scholar 

  40. Picillo M, Lozano AM, Kou N, Munhoz RP, Fasano A. Programming deep brain stimulation for Parkinson’s disease: the Toronto Western Hospital algorithms. Brain Stimulation. 2016; 9(3):425–37.

    Article  PubMed  Google Scholar 

  41. Goetz CG, Tilley BC, Shaftman SR, Stebbins GT, Fahn S, Martinez-Martin P, Poewe W, Sampaio C, Stern MB, Dodel R, et al. Movement Disorder Society-sponsored revision of the Unified Parkinson’s Disease Rating Scale (MDS-UPDRS): scale presentation and clinimetric testing results. Mov Disord. 2008; 23(15):2129–70.

    Article  PubMed  Google Scholar 

  42. Sabatini U, Boulanouar K, Fabre N, Martin F, Carel C, Colonnese C, Bozzao L, Berry I, Montastruc J, Chollet F, et al. Cortical motor reorganization in akinetic patients with Parkinson’s disease: a functional MRI study. Brain. 2000; 123(2):394–403.

    Article  PubMed  Google Scholar 

  43. Grafton S, Turner R, Desmurget M, Bakay R, Delong M, Vitek J, Crutcher M. Normalizing motor-related brain activity: subthalamic nucleus stimulation in Parkinson disease. Neurology. 2006; 66(8):1192–9.

    Article  CAS  PubMed  Google Scholar 

  44. Horn A, Wenzel G, Irmen F, Huebl J, Li N, Neumann W-J, Krause P, Bohner G, Scheel M, Kühn AA. Deep brain stimulation induced normalization of the human functional connectome in Parkinson’s disease. Brain. 2019; 142(10):3129–43.

    Article  PubMed  Google Scholar 

  45. Hassan M, Chaton L, Benquet P, Delval A, Leroy C, Plomhause L, Moonen AJ, Duits AA, Leentjens AF, van Kranen-Mastenbroek V, et al. Functional connectivity disruptions correlate with cognitive phenotypes in Parkinson’s disease. NeuroImage: Clin. 2017; 14:591–601.

    Article  CAS  Google Scholar 

  46. Yuvaraj R, Murugappan M, Acharya UR, Adeli H, Ibrahim NM, Mesquita E. Brain functional connectivity patterns for emotional state classification in Parkinson’s disease patients without dementia. Behav Brain Res. 2016; 298:248–60.

    Article  CAS  PubMed  Google Scholar 

  47. Yu N, Estévez N, Hepp-Reymond M-C, Kollias SS, Riener R. fMRI assessment of upper extremity related brain activation with an MRI-compatible manipulandum. Int J Comput Assisted Radiol Surg. 2011; 6(3):447–55.

    Article  Google Scholar 

  48. Jaeger L, Marchal-Crespo L, Wolf P, Riener R, Kollias S, Michels L. Test-retest reliability of fMRI experiments during robot-assisted active and passive stepping. J Neuroeng Rehabil. 2015; 12(1):102.

    Article  PubMed  PubMed Central  Google Scholar 

  49. Koss AM, Alterman RL, Tagliati M, Shils JL. Calculating total electrical energy delivered by deep brain stimulation systems. Ann Neurol. 2005; 58(1):168.

    Article  PubMed  Google Scholar 

  50. Hui D, Murgai AA, Gilmore G, Mohideen SI, Parrent AG, Jog MS. Assessing the effect of current steering on the total electrical energy delivered and ambulation in Parkinson’s disease. Sci Rep. 2020; 10(1):1–11.

    Article  CAS  Google Scholar 

Download references


The authors would like to thank Ms. Yuanyuan Cheng and Ms. Yue Wang from Tianjin Huanhu Hospital for their help in assessment of the task performance videos. The authors would express their sincere gratitude to all the patients that participated in this study and made the discovery possible.


This work was supported by the National Natural Science Foundation of China (U1913208, 61873135, 61720106012) and the fundamental research funds for the central universities.

Author information

Authors and Affiliations



Authors’ contributions

N. Yu, J. Lu, and S. Liang designed the experiments and interpreted the results. N. Yu, J. Lu, Z. Shu, H. Li, Y. Yu, and S. Liang performed the experiments with input from J. Wu and J. Han. N. Yu, J. Lu, and J. Han analyzed the data and wrote the manuscript. All authors provided feedback on the manuscript. The author(s) read and approved the final manuscript.

Authors’ information

Not applicable.

Corresponding authors

Correspondence to Jialing Wu or Jianda Han.

Ethics declarations

Ethics approval and consent to participate

The study was approved by the Ethical Committee of Tianjin Huanhu Hospital (2019-35) and has been registered in Chinese Clinical Trial Registry (ChiCTR1900022715). Each patient was fully informed of the experimental purpose and procedures and provided written consent prior to the measurement.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Yu, N., Liang, S., Lu, J. et al. Quantified assessment of deep brain stimulation on Parkinson’s patients with task fNIRS measurements and functional connectivity analysis: a pilot study. Chin Neurosurg Jl 7, 34 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: