Skip to main content

ORIGINAL RESEARCH article

Front. Hum. Neurosci., 19 March 2019
Sec. Brain Imaging and Stimulation
This article is part of the Research Topic Techniques Advances and Clinical Applications in Fused EEG-fNIRS View all 9 articles

Cortical Hemodynamic Response and Connectivity Modulated by Sub-threshold High-Frequency Repetitive Transcranial Magnetic Stimulation

\r\nRihui LiRihui Li1Thomas PotterThomas Potter1Jun WangJun Wang2Zhixi ShiZhixi Shi2Chushan WangChushan Wang2Lingling YangLingling Yang3Rosa ChanRosa Chan4Yingchun Zhang*Yingchun Zhang1*
  • 1Department of Biomedical Engineering, University of Houston, Houston, TX, United States
  • 2Guangdong Provincial Work Injury Rehabilitation Hospital, Guangzhou, China
  • 3Department of Computer Science, Sun Yat-sen University, Guangzhou, China
  • 4Department of Electrical Engineering, City University of Hong Kong, Hong Kong, China

Repetitive transcranial magnetic stimulation (rTMS) at sub-threshold intensity is a viable clinical strategy to enhance the sensory and motor functions of extremities by increasing or decreasing motor cortical excitability. Despite this, it remains unclear how sub-threshold rTMS modulates brain cortical excitability and connectivity. In this study, we applied functional near-infrared spectroscopy (fNIRS) to investigate the alterations in hemodynamic responses and cortical connectivity patterns that are induced by high-frequency rTMS at a sub-threshold intensity. Forty high-frequency (10 Hz) trains of rTMS at 90% resting motor threshold (RMT) were delivered through a TMS coil placed over 1–2 cm lateral from the vertex. fNIRS signals were acquired from the frontal and bilateral motor areas in healthy volunteers (n = 20) during rTMS administration and at rest. A significant reduction in oxygenated hemoglobin (HbO) concentration was observed in most defined regions of interest (ROIs) during the stimulation period (p < 0.05). Decreased functional connectivity within prefrontal areas as well as between symmetrical ROI-pairs was also observed in most participants during the stimulation (p < 0.05). Results suggest that fNIRS imaging is able to provide a reliable measure of regional cortical brain activation that advances our understanding of the manner in which sub-threshold rTMS affects cortical excitability and brain connectivity.

Introduction

Transcranial magnetic stimulation (TMS) is a therapeutic technique that applies a rapid alternating current to generate magnetic fields, which penetrate scalp and skull to reach the brain and induce secondary electrical fields that activate the cortical neurons (Wassermann and Lisanby, 2001). When TMS is provided periodically at a specific frequency, it is referred to as repetitive TMS (rTMS; Kobayashi and Pascual-Leone, 2003). Specifically, rTMS inhibits cortical excitability at lower pulse frequencies (less than the 1 Hz) and enhances the cortical excitability at higher pulse frequencies (5 Hz or higher; Ridding and Ziemann, 2010; Valero-Cabré et al., 2017). Since its introduction in 1985 (Barker et al., 1985), rTMS has demonstrated the ability to induce long-term changes in cortical activity (Fitzgerald et al., 2006b), leading to its adoption as a potential treatment for a variety of neurological and mental disorders, including depression and motor impairment (Takeuchi et al., 2005; Fitzgerald et al., 2006a; O’Reardon et al., 2007; Schwerin et al., 2011). The underlying mechanism of action for rTMS, however, remains poorly understood. In order to maximize the therapeutic potential of rTMS, it is necessary to investigate the effects that stimulation has on both brain activity and cortical connectivity.

To date, a number of neuroimaging modalities have been used to investigate the functional changes associated with TMS. The most dominant technological investigations have combined TMS with functional MRI (fMRI; Bohning et al., 1999; Bestmann et al., 2008; Wang et al., 2017), electroencephalography (EEG; Bonato et al., 2006; Pigorini et al., 2011; Chung et al., 2015), and magnetoencephalography (MEG; Shibasaki, 2008), or functional near infrared spectroscopy (fNIRS; Mochizuki et al., 2006; Koenraadt et al., 2011; Cao et al., 2013). Each modality has its strengths and weaknesses. FMRI, EEG, and MEG all measure primary or secondary electromagnetic signals that are associated with various brain activities. The magnetic and electrical fields produced by TMS may interfere with these measurements and introduce noise and error to assessments of TMS response. FNIRS, on the other hand, is an emerging noninvasive optical imaging technique that offers an alternative way to investigate TMS-evoked brain activity (Boas et al., 2004). It estimates activity-linked changes in tissue hemoglobin (Hb) concentration based on the different absorption coefficients for oxygenated hemoglobin (HbO) and deoxygenated hemoglobin (HbR; Bunce et al., 2006; Ernst et al., 2012). fNIRS devices with electromagnetism-compatible designs do not interfere with the electromagnetic fields produced by TMS. Consequently, fNIRS holds several advantages over other imaging modalities: a high signal-to-noise ratio, non-interference with electromagnetic fields, portability, and reduced sensitivity to motion artifacts (Li et al., 2017). These features make fNIRS a promising tool to investigate the cortical activation and connectivity associated with rTMS (Kozel et al., 2009; Cao et al., 2013; Groiss et al., 2013).

Previous studies have demonstrated the feasibility of using fNIRS to characterize the effects of rTMS on the cortical excitability and connectivity. Mochizuki et al. examined the HbO, HbR, and total hemoglobin (HBT) change when 1 Hz rTMS was delivered to the motor cortex at 100%, 120%, and 140% of the active motor threshold (AMT, defined as the lowest TMS intensity necessary to induce reliable motor evoke potentials in the target muscle during slightly active contraction; Mochizuki et al., 2006). Results suggested that hemodynamic responses to lower stimulation intensities showed no significant change, while hemodynamic responses at 140% of AMT revealed a significant reduction in both HbR and HbT. Another study focusing on the effects of 1 Hz rTMS on non-motor areas (prefrontal cortex) demonstrated that there was no significant change induced by TMS at 90% and 110% of the resting motor threshold (RMT, defined as the lowest TMS intensity necessary to induce reliable motor evoke potentials in the target muscle at resting), while a significant drop in HbO was induced by 130% RMT stimulation (Thomson et al., 2011b). The impact of rTMS on the cortical connectivity has also been investigated in previous studies. Kozel et al. (2009) measured the fNIRS signals from ipsilateral and contralateral brain regions, during which subjects underwent stimulation (1 Hz, 120% RMT) of the motor cortex and, subsequently, the prefrontal cortex. They reported that stimulation to either the motor and prefrontal cortices resulted in a significant reduction in HbO concentration in both the ipsilateral and contralateral cortices, and the ipsilateral and contralateral changes showed high temporal consistency.

While the majority of previous studies reported the observable regional cerebral blood flow (rCBF) at specific brain regions during rTMS intervention at supra-threshold intensity, the broad impact of sub-threshold rTMS stimulation on multiple brain regions has not been well documented. This may be due to the mild effect of sub-threshold rTMS on brain activation as well as the variation in analytical approaches (stimulation parameters and data analysis methods; Noguchi et al., 2003). In addition, it remains unknown whether sub-threshold rTMS has a significant effect on the connectivity between multiple brain regions. This presents a clear need to explore the full extent of hemodynamic responses to TMS at sub-threshold intensities.

In this study, we aimed to investigate the effects of high-frequency (10 Hz) sub-threshold rTMS on brain cortical excitability and connectivity across multiple brain regions. The global hemodynamic trend of the brain was monitored and change in regional connectivity was evaluated by statistically testing the Pearson correlation coefficients before and during the rTMS intervention. We hypothesized that alterations in brain excitability and connectivity can be induced in multiple cortical regions by high-frequency rTMS at sub-threshold intensity. We also expected that the quantified index of hemodynamic responses proposed in this study could be adjusted to optimize the setting during rTMS treatment and translated to clinical applications.

Materials and Methods

Participants

Twenty healthy, right-handed volunteers (16 male, four female, age of 25.65 ± 1.69) were recruited from the local community to take part in this study. All participants were free of current neurological or psychiatric disorders. None of the participants had previous experience with the study paradigm employed in this study. Each subject was fully informed about the purpose and risks of the research, and provided written, informed consent prior to the start of the experiment. All protocols and procedures were approved by the local ethics committee (Guangdong Provincial Work Injury Rehabilitation Center, China) and performed in accordance with the Declaration of Helsinki.

Concurrent fNIRS/rTMS Configuration

A concurrent fNIRS-rTMS measurement setup was employed in this study. A near infrared spectroscopy system (NIRx Medical Technologies, LLC, Glen Head, NY, USA) consisting of 16 emitters and 15 detectors was utilized to measure hemodynamic responses during rTMS administration. An fNIRS channel was defined as a pair comprised of an emitter and a detector, and all emitters and detectors were mounted by plastic grommets and holders in an elastic cap according to the international 10–20 EEG placement system, as shown in Figure 1A. The inter-optode distance was about 3 cm and a total of 46 fNIRS channels were equidistantly distributed throughout the frontal and bilateral motor areas and were steadily held in place with a plastic strap (Figure 1A). The wavelengths used for the detection of the concentration changes of HbO and HbR were 760 and 850 nm, and all fNIRS signals were acquired at a sampling rate of 3.91 Hz.

FIGURE 1
www.frontiersin.org

Figure 1. Experiment setup. (A) The configuration of concurrent functional near-infrared spectroscopy (fNIRS)-repetitive transcranial magnetic stimulation (rTMS) measurement. (B) The protocol of the rTMS. The whole stimulation period lasted about 20 min.

A MagStim Super Rapid magnetic stimulator (MagStim Company, Whitland, Wales, UK) equipped with a commercially available 70-mm diameter figure-eight coil was used to perform rTMS in this study. As shown in Figures 1A, 2, the TMS coil was steadily held and fixed by a metal arm over 1–2 cm lateral from the vertex (Cz channel) of the participant as they lay in a supine position, aligning with the leg representation of the motor cortex (M1; Groppa et al., 2012). For both motor threshold determination and subsequent rTMS, magnetic stimulation was delivered by steadily placing the TMS coil over the fNIRS probes.

FIGURE 2
www.frontiersin.org

Figure 2. The overview of the whole experimental scheme, including the TMS machine (A), TMS coil (B,D) and fNIRS system (C). The subject included in the figure provided written, informed consent for the publication of this figure.

RTMS Protocol

The RMT of each subject was measured and determined prior to the experiment. The RMT was defined as the lowest stimulus intensity that generated motor evoked potentials (MEPs) in the right tibialis anterior muscle (TAM) that were larger than 50 microvolts in 5 of 10 consecutive stimulations (Groppa et al., 2012). The intensity of rTMS in this study was then set at 90% of the RMT intensity. As shown in Figure 1B, the formal experimental session began with a 5 min baseline measurement, during which the rTMS coil was set up and placed over 1–2 cm lateral from the vertex without stimulation delivery. After the baseline measurement, sub-threshold rTMS was performed. Forty trains of TMS were administered over 20 min at 90% RMT. Each train lasted for 30 s and consisted of a 2 s stimulation at 10 Hz frequency and a 28 s break period. FNIRS data were simultaneously collected throughout the whole session, including the baseline and stimulation period. The rTMS protocol used in this study was a clinically standard protocol approved by the collaborated hospital. A complete paradigm of the experiment setup is illustrated in Figure 2.

Data Preprocessing

In this study, the global temporal change in the hemodynamic properties induced by rTMS was investigated. The global trend was analyzed based on the time course of measurements over the whole TMS experiment (20 min). Acquired fNIRS data were processed and analyzed for each individual subject using MATLAB (2014a, MathWorks, Natick, MA, USA). Considering the global trend was an extremely slow fluctuation (<0.01 Hz), raw fNIRS data were first filtered with a 4th order Butterworth low-pass filter with a high cutoff frequency of 0.2 Hz. In addition, spline interpolation was employed to remove any motion artifact contamination from the fNIRS signal (Scholkmann et al., 2010). The baseline period was defined as the 3 min period immediately before the onset of rTMS administration. Time courses of the changes in HbO, HbR, and HbT concentration throughout the whole stimulation session were calculated relative to the baseline period using the Modified Beer-Lambert Law (Scholkmann et al., 2014). The filtered data were further detrended to remove potential linear drift artifacts from the time series signals.

Brain Regions Selection

The effects of rTMS on eight regions of interest (ROIs) of the brain were investigated and assessed in this study. These regions included the left prefrontal area (LPF), middle prefrontal area (MPF), right prefrontal (RPF), supplementary motor area (SMA), left premotor cortex (LPMC), right premotor cortex (RPMC), left primary motor cortex (LM1), and right primary motor cortex (RM1). The details of these ROIs and the corresponding fNIRS channels are summarized in Figure 3. A fixed number of fNIRS channels were assigned to each ROI consistently across all subjects.

FIGURE 3
www.frontiersin.org

Figure 3. Summary of the segmented regions of interest (ROIs) and assigned fNIRS channels. LPF, Left prefrontal; RPF, Right prefrontal; MPF, Middle prefrontal; SMA, Supplementary motor area; LPMC, Left premotor cortex; RPMC, Right premotor cortex; LM1, Left primary motor cortex; RM1, Right primary motor cortex.

Statistical Analysis of Hemodynamic Response

Statistical analysis focused primarily on the HbO signals; while both HbO and HbR were recorded, HbO is considered to be a more robust and sensitive fNIRS parameter than HbR to stimulus-associated changes, and shows stronger correlation with fMRI BOLD response (Plichta et al., 2006; Cui et al., 2011). The signals from the fNIRS channels in each ROI were pooled and averaged to characterize the regional hemodynamic responses. A sliding window approach (60 s window with a 30 s overlap) was then used to smooth the pooled HbO time series of all ROIs for global trend analysis. This method revealed the common temporal pattern of the cortical hemodynamic response by minimizing transient temporal variations. Finally, the regional hemodynamic response induced by rTMS was quantified by calculating the mean area under curve (MAUC) over the entire 20-min pooled HbO time series. The hemodynamic response of the resting condition was also computed as the MAUC within the 3-min baseline period. The difference of MAUC between the two conditions at each region was subsequently tested by paired t-test with Bonferroni multiple comparison correction (pcorrected = 0.0063).

Functional Connectivity Analysis

In addition to exploring the temporal changes in HbO, the changes in regional connectivity induced by rTMS were evaluated via correlation analysis. Connectivity was calculated during the 3-min resting period and 20-min stimulation session, aligning with the lengthy analytical periods used in previous studies (Mesquita et al., 2010; Rosenbaum et al., 2017). The Pearson correlation coefficient between each ROI pair was computed from the HbO time courses during both the resting period and stimulation session, respectively, yielding 28 correlation values for each condition and each subject. The differences in regional connectivity between the resting and stimulation conditions were statistically evaluated using a series of two-tailed t-tests with Bonferroni multiple comparison correction to reflect the effects of rTMS on regional brain connectivity (pcorrected = 0.0018). Prior to statistical testing, all Pearson correlation coefficients were transformed using the Fisher z transformation (Fisher, 1915).

Results

Global Effect of the rTMS

Figure 4 shows the global temporal pattern and 95% confidence interval of the HbO concentration in each ROI during the resting period and 20-min rTMS session. As shown in Figure 4, during the rTMS session, the HbO concentrations generally exhibited a noticeable reduction followed by a recovery phase in most ROIs, indicating the effects of rTMS on the hemodynamic response. Statistical testing of the global HbO change (MAUC) revealed a highly significant difference between resting period and the entire 20-min stimulation period in all ROIs (p < pcorrected), with the single exception of SMA (p = 0.045), as shown in Figure 5.

FIGURE 4
www.frontiersin.org

Figure 4. The group-averaged temporal patterns of oxygenated hemoglobin (HbO) concentration changes (black lines) during rest period and the entire 20-min rTMS session in all ROIs. The upper and lower dotted lines denote 95% confidence intervals, while the vertical dash lines denote the start and end of the stimulation.

FIGURE 5
www.frontiersin.org

Figure 5. A summary of the mean area under curves (MAUCs) for global temporal HbO concentration changes before and during rTMS. The asterisk “*” indicates a significant difference between the rest period and stimulation session (p < 0.05).

Functional Connectivity Analysis

Prior to stimulation, the average correlation coefficient of all ROI pairs was 0.674 ± 0.116. During the whole stimulation, however, the average correlation coefficient of all ROI pairs was reduced to 0.615 ± 0.096. This indicated that most ROIs exhibited a decrease in inter-regional connectivity during rTMS administration. The alteration of functional connectivity induced by rTMS was then assessed for each ROI-pair. Figure 6 shows the connections that were significantly suppressed by the rTMS administration. Briefly, the connectivity between the RPF and MPF (p = 0.0011) as well as between the RPF and left frontal (LPF; p < 0.001) was significantly suppressed by the rTMS, as indicated by the white arrows in Figure 6. In addition to the connections within the prefrontal area, suppression of connectivity was also observed between LPF, RPF and SMA (pLPF-SMA = 0.0017, pRPF-SMA = 0.001), respectively.

FIGURE 6
www.frontiersin.org

Figure 6. Statistical test results of functional connectivity analysis. The circles denote the defined ROIs, while white arrows denote connections that were significantly decreased during the rTMS session relative to rest period.

Discussion

In the current study, the time-sensitive changes in global brain hemodynamics induced by sub-threshold high frequency rTMS were investigated at multiple cortical areas using fNIRS. A significant global reduction of HbO concentration was observed at all selected ROIs except SMA. In addition, the functional connectivity changes induced by rTMS were also examined for each ROI-pair and the pattern of these changes was characterized. The reported results demonstrate the feasibility and capability of using fNIRS to noninvasively assess changes in cerebral hemodynamic responses and inter-regional cortical connectivity during rTMS administration.

The hemodynamic response induced by sub-threshold rTMS in single trial is generally very mild and can be masked when averaging epochs, presenting a challenge to fNIRS-based measurements (Thomson et al., 2011b). In this study, the global changes in HbO concentration were used to investigate the effect of sub-threshold rTMS on the hemodynamic responses and connectivity of the brain. Significant HbO changes were observed in most of the specified brain regions during the rTMS session. The observed regional HbO concentrations typically decreased after the onset of the stimulation and gradually recovered back to the baseline over time, indicating an accumulative effect of the sub-threshold rTMS. This finding is in line with previous fMRI studies, where multiple areas surrounding the stimulated area were deactivated (Jung et al., 2016). A potential explanation for the observed reduction of HbO might be the byproduct effect of rTMS, in which the sensory responses induced by an rTMS pulse may capture the participant’s attention and thereby interrupt normal conscious activity. Another possible explanation for the global reduction in HbO is the concept of “vascular steal”; if the sub-threshold rTMS produced an increase in oxygen consumption at the directly stimulated area that exceeded the local supply, oxygenated blood may have been shunted to the region, resulting in the net drop in the HbO observed at other brain regions (Hu and Huang, 2015). Thomson et al. also reported a similar phenomenon when supra-threshold (120% RMT) 1 Hz rTMS was delivered to PFC, though they only observed a mild early reduction instead of a significant reduction in HbO when a sub-threshold (80% RMT) 1 Hz rTMS was applied (Thomson et al., 2012). This discrepancy may be due to the differences in the stimulus intensity and target location for the studies. A more comprehensive investigation may be necessary to address this in the future.

Interestingly, it is noteworthy that the reduction of the HbO induced by rTMS did not last for the entire rTMS session, as shown in Figure 4. The gradual recovery of the HbO at the end of the rTMS session suggests that the stimulation became less impactful on cortical brain regions as it continued. This may be due to the brain regions gradually adapting to the impact of rTMS or the restoration of regular vascular oxygenation through constant cerebral perfusion. As described in previous studies (Pecuch et al., 2000; Vernieri et al., 2009), vascular innervation represents one of the main regulating factors of cerebral circulation. The activation of the sympathetic nervous system under certain stimulation maintains a constant CBF through vasoconstriction, which might prevent the potential of dilatation of cerebral vessels induced by rTMS in this study. We might hypothesize that long-term rTMS administration might not be beneficial in clinical therapy. However, previous rTMS studies have demonstrated that the effect of stimulation on cortical excitability could linger for a few minutes after stimulation (Romero et al., 2002; Aydin-Abidin et al., 2006). As the direct relationship between hemodynamic response and cortical excitability during TMS administration is absent, whether fNIRS is able to evaluate the duration of rTMS should be further validated.

As previous studies have reported, the effects of rTMS extend beyond the “hot spot” area to influence the functional connectivity among distant brain regions (Watanabe et al., 2014; Valero-Cabré et al., 2017). For instance, Watanabe et al. reported that the functional connectivity between multiple regions assessed by fMRI was significantly reduced by the high frequency rTMS (Watanabe et al., 2014). Our results confirmed this finding, as the application of rTMS significantly decreased the functional connectivity between multiple ROI-pairs in most subjects, demonstrating the ability of rTMS to regulate global brain networks. In addition, the significant disruptions observed within the prefrontal areas (LPF, MPF and RPF) suggested the distant prefrontal cortex may be more sensitive to stimulation compared to motor areas, which is partly in line with the finding of a previous study that observed no alteration among motor network during the vertex stimulation (Jung et al., 2016). Furthermore, the suppressed connections between RPF, LPF and SMA regions imply that rTMS seemed to modify connectivity between symmetric cortical regions in our stimulation protocol. A previous fNIRS-based functional connectivity study reported a symmetrical correlation between the two hemispheres in healthy controls, which represents a regular brain connectivity pattern during resting state (Mesquita et al., 2010). In our study, the alterations in connectivity between symmetrical ROI pairs induced by rTMS may be used to evaluate the efficiency of the rTMS treatment in a clinical setting. For example, potential clinical applications include using rTMS and fNIRS to assess the alterations in patterns of brain activity and interhemispheric connectivity that accompany specific brain disorders (e.g., mental disorder or stroke). In this context, the degree of the characterized activity and connectivity changes could be used to determine the optimal settings of the TMS treatment for specific patients. Similarly, fNIRS and rTMS may be applied to better understand the underlying mechanism of disease progression and monitor the progress of recovering patients.

While the present study provides important details regarding the effects of TMS application on regional blood flow and connectivity, one limitation is that a sham group was not recruited in the experiment, making it difficult to validate if the alterations of hemodynamic response and the connectivity were specifically associated with the rTMS. However, because the hemodynamic response from the sham group has been repeatedly confirmed in previous studies (Thomson et al., 2011a,b; Park et al., 2017), we believe the observed abnormal hemodynamic response could very likely be caused by the rTMS. Considering that the limited sample size in this article may have been insufficient to draw a solid conclusion for its clinical value, we will include a sham group and more subjects in the future to validate the current findings and the feasibility of translating the proposed approach to clinical applications. Additionally, due to a limited number of optodes, in this study our fNIRS device was not able to provide full head coverage to measure comprehensive hemodynamic response signals. It is therefore suggested that future research should utilize an improved setup with more fNIRS optodes. It should also be noted that fNIRS is only able to measure the hemodynamic response on the cortical surface (1–3 cm) due to the inherent limited penetration depth; it would thereby be difficult to investigate any potential hemodynamic changes in subcortical or deep brain areas.

Conclusion

Global alterations in the hemodynamic response and connectivity between multiple brain regions induced by sub-threshold 10 Hz rTMS were investigated in this study. The HbO concentration in most regions was found to decrease after the onset of the stimulation and gradually recover to the baseline as stimulation continued. Results also showed that rTMS modified the functional connectivity within the prefrontal area and between symmetrical ROI-pairs. Results demonstrated the feasibility of using the combination of rTMS and fNIRS to evaluate both cortical excitability and connectivity, which may pave the way for optimized, patient-specific rTMS treatment programs.

Author Contributions

RL, JW, ZS and YZ: study design. RL, CW and ZS: data acquisition. RL, TP, and LY: analysis and interpretation. RL, TP, RC, and YZ: manuscript drafting. RL, TP, JW, ZS, CW, LY, RC, and YZ: final approval.

Funding

This work was supported in part by the University of Houston and the Guangdong Provincial Work Injury Rehabilitation Center.

Conflict of Interest Statement

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.

References

Aydin-Abidin, S., Moliadze, V., Eysel, U. T., and Funke, K. (2006). Effects of repetitive TMS on visually evoked potentials and EEG in the anaesthetized cat: dependence on stimulus frequency and train duration. J. Physiol. 574, 443–455. doi: 10.1113/jphysiol.2006.108464

PubMed Abstract | CrossRef Full Text | Google Scholar

Barker, A. T., Jalinous, R., and Freeston, I. L. (1985). Non-invasive magnetic stimulation of human motor cortex. Lancet 1, 1106–1107. doi: 10.1016/s0140-6736(85)92413-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Bestmann, S., Ruff, C. C., Blankenburg, F., Weiskopf, N., Driver, J., and Rothwell, J. C. (2008). Mapping causal interregional influences with concurrent TMS-fMRI. Exp. Brain Res. 191, 383–402. doi: 10.1007/s00221-008-1601-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Boas, D. A., Dale, A. M., and Franceschini, M. A. (2004). Diffuse optical imaging of brain activation: approaches to optimizing image sensitivity, resolution, and accuracy. Neuroimage 23, S275–S288. doi: 10.1016/j.neuroimage.2004.07.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Bohning, D. E., Shastri, A., McConnell, K. A., Nahas, Z., Lorberbaum, J. P., Roberts, D. R., et al. (1999). A combined TMS/fMRI study of intensity-dependent TMS over motor cortex. Biol. Psychiatry 45, 385–394. doi: 10.1016/s0006-3223(98)00368-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Bonato, C., Miniussi, C., and Rossini, P. M. (2006). Transcranial magnetic stimulation and cortical evoked potentials: a TMS/EEG co-registration study. Clin. Neurophysiol. 117, 1699–1707. doi: 10.1016/j.clinph.2006.05.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Bunce, S. C., Izzetoglu, M., Izzetoglu, K., Onaral, B., and Pourrezaei, K. (2006). Functional near-infrared spectroscopy—an emerging neuroimaging modality. IEEE Eng. Med. Biol. 25, 54–62. doi: 10.1109/MEMB.2006.1657788

CrossRef Full Text | Google Scholar

Cao, T. T., Thomson, R. H., Bailey, N. W., Rogasch, N. C., Segrave, R. A., Maller, J. J., et al. (2013). A near infra-red study of blood oxygenation changes resulting from high and low frequency repetitive transcranial magnetic stimulation. Brain Stimul. 6, 922–924. doi: 10.1016/j.brs.2013.04.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Chung, S. W., Rogasch, N. C., Hoy, K. E., and Fitzgerald, P. B. (2015). Measuring brain stimulation induced changes in cortical properties using TMS-EEG. Brain Stimul. 8, 1010–1020. doi: 10.1016/j.brs.2015.07.029

PubMed Abstract | CrossRef Full Text | Google Scholar

Cui, X., Bray, S., Bryant, D. M., Glover, G. H., and Reiss, A. L. (2011). A quantitative comparison of NIRS and fMRI across multiple cognitive tasks. Neuroimage 54, 2808–2821. doi: 10.1016/j.neuroimage.2010.10.069

PubMed Abstract | CrossRef Full Text | Google Scholar

Ernst, L. H., Schneider, S., Ehlis, A. C., and Fallgatter, A. J. (2012). Functional near infrared spectroscopy in psychiatry: a critical review. J. Near Infrared Spec. 20, 93–105. doi: 10.1255/jnirs.970

CrossRef Full Text | Google Scholar

Fisher, R. A. (1915). Frequency distribution of the values of the correlation coefficient in samples from an indefinitely large population. Biometrika 10, 507–521. doi: 10.2307/2331838

CrossRef Full Text | Google Scholar

Fitzgerald, P. B., Benitez, J., de Castella, A., Daskalakis, Z. J., Brown, T. L., and Kulkarni, J. (2006a). A randomized, controlled trial of sequential bilateral repetitive transcranial magnetic stimulation for treatment-resistant depression. Am. J. Psychiatry 163, 88–94. doi: 10.1176/appi.ajp.163.1.88

PubMed Abstract | CrossRef Full Text | Google Scholar

Fitzgerald, P. B., Fountain, S., and Daskalakis, Z. J. (2006b). A comprehensive review of the effects of rTMS on motor cortical excitability and inhibition. Clin. Neurophysiol. 117, 2584–2596. doi: 10.1016/j.clinph.2006.06.712

PubMed Abstract | CrossRef Full Text | Google Scholar

Groiss, S. J., Mochizuki, H., Furubayashi, T., Kobayashi, S., Nakatani-Enomoto, S., Nakamura, K., et al. (2013). Quadri-pulse stimulation induces stimulation frequency dependent cortical hemoglobin concentration changes within the ipsilateral motor cortical network. Brain Stimul. 6, 40–48. doi: 10.1016/j.brs.2011.12.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Groppa, S., Oliviero, A., Eisen, A., Quartarone, A., Cohen, L. G., Mall, V., et al. (2012). A practical guide to diagnostic transcranial magnetic stimulation: report of an IFCN committee. Clin. Neurophysiol. 123, 858–882. doi: 10.1016/j.clinph.2012.01.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Hu, D. W., and Huang, L. M. (2015). Negative hemodynamic response in the cortex: evidence opposing neuronal deactivation revealed via optical imaging and electrophysiological recording. J. Neurophysiol. 114, 2152–2161. doi: 10.1152/jn.00246.2015

PubMed Abstract | CrossRef Full Text | Google Scholar

Jung, J., Bungert, A., Bowtell, R., and Jackson, S. R. (2016). Vertex stimulation as a control site for transcranial magnetic stimulation: a concurrent TMS/fMRI study. Brain Stimul. 9, 58–64. doi: 10.1016/j.brs.2015.09.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Kobayashi, M., and Pascual-Leone, A. (2003). Transcranial magnetic stimulation in neurology. Lancet Neurol. 2, 145–156. doi: 10.1016/S1474-4422(03)00321-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Koenraadt, K. L. M., Munneke, M. A. M., Duysens, J., and Keijsers, N. L. W. (2011). TMS: a navigator for NIRS of the primary motor cortex? J. Neurosci. Methods 201, 142–148. doi: 10.1016/j.jneumeth.2011.07.024

PubMed Abstract | CrossRef Full Text | Google Scholar

Kozel, F. A., Tian, F. H., Dhamne, S., Croarkin, P. E., McClintock, S. M., Elliott, A., et al. (2009). Using simultaneous repetitive Transcranial Magnetic Stimulation/functional Near Infrared Spectroscopy (rTMS/fNIRS) to measure brain activation and connectivity. Neuroimage 47, 1177–1184. doi: 10.1016/j.neuroimage.2009.05.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, R. H., Potter, T., Huang, W. T., and Zhang, Y. C. (2017). Enhancing performance of a hybrid EEG-fNIRS system using channel selection and early temporal features. Front. Hum. Neurosci. 11:462. doi: 10.3389/fnhum.2017.00462

PubMed Abstract | CrossRef Full Text | Google Scholar

Mesquita, R. C., Franceschini, M. A., and Boas, D. A. (2010). Resting state functional connectivity of the whole head with near-infrared spectroscopy. Biomed. Opt. Express 1, 324–336. doi: 10.1364/boe.1.000324

PubMed Abstract | CrossRef Full Text | Google Scholar

Mochizuki, H., Ugawa, Y., Terao, Y., and Sakai, K. L. (2006). Cortical hemoglobin-concentration changes under the coil induced by single-pulse TMS in humans: a simultaneous recording with near-infrared spectroscopy. Exp. Brain Res. 169, 302–310. doi: 10.1007/s00221-005-0149-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Noguchi, Y., Watanabe, E., and Sakai, K. L. (2003). An event-related optical topography study of cortical activation induced by single-pulse transcranial magnetic stimulation. Neuroimage 19, 156–162. doi: 10.1016/s1053-8119(03)00054-5

PubMed Abstract | CrossRef Full Text | Google Scholar

O’Reardon, J. P., Solvason, H. B., Janicak, P. G., Sampson, S., Isenberg, K. E., Nahas, Z., et al. (2007). Efficacy and safety of transcranial magnetic stimulation in the acute treatment of major depression: a multisite randomized controlled trial. Biol. Psychiatry 62, 1208–1216. doi: 10.1016/j.biopsych.2007.01.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Park, E., Kang, M. J., Lee, A., Chang, W. H., Shin, Y. I., and Kim, Y. H. (2017). Real-time measurement of cerebral blood flow during and after repetitive transcranial magnetic stimulation: a near-infrared spectroscopy study. Neurosci. Lett. 653, 78–83. doi: 10.1016/j.neulet.2017.05.039

PubMed Abstract | CrossRef Full Text | Google Scholar

Pecuch, P. W., Evers, S., Folkerts, H. W., Michael, N., and Arolt, V. (2000). The cerebral hemodynamics of repetitive transcranial magnetic stimulation. Eur. Arch. Psychiatry Clin. Neurosci. 250, 320–324. doi: 10.1007/s004060070007

PubMed Abstract | CrossRef Full Text | Google Scholar

Pigorini, A., Casali, A. G., Casarotto, S., Ferrarelli, F., Baselli, G., Mariotti, M., et al. (2011). Time-frequency spectral analysis of TMS-evoked EEG oscillations by means of Hilbert-Huang transform. J. Neurosci. Methods 198, 236–245. doi: 10.1016/j.jneumeth.2011.04.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Plichta, M. M., Herrmann, M. J., Baehne, C. G., Ehlis, A. C., Richter, M. M., Pauli, P., et al. (2006). Event-related functional near-infrared spectroscopy (fNIRS): are the measurements reliable? Neuroimage 31, 116–124. doi: 10.1016/j.neuroimage.2005.12.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Ridding, M. C., and Ziemann, U. (2010). Determinants of the induction of cortical plasticity by non-invasive brain stimulation in healthy subjects. J. Physiol. 588, 2291–2304. doi: 10.1113/jphysiol.2010.190314

PubMed Abstract | CrossRef Full Text | Google Scholar

Romero, J. R., Anschel, D., Sparing, R., Gangitano, M., and Pascual-Leone, A. (2002). Sub-threshold low frequency repetitive transcranial magnetic stimulation selectively decreases facilitation in the motor cortex. Clin. Neurophysiol. 113, 101–107. doi: 10.1016/s1388-2457(01)00693-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Rosenbaum, D., Haipt, A., Fuhr, K., Haeussinger, F. B., Metzger, F. G., Nuerk, H. C., et al. (2017). Aberrant functional connectivity in depression as an index of state and trait rumination. Sci. Rep. 7:2174. doi: 10.1038/s41598-017-02277-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Scholkmann, F., Kleiser, S., Metz, A. J., Zimmermann, R., Mata Pavia, J., Wolf, U., et al. (2014). A review on continuous wave functional near-infrared spectroscopy and imaging instrumentation and methodology. Neuroimage 85, 6–27. doi: 10.1016/j.neuroimage.2013.05.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Scholkmann, F., Spichtig, S., Muehlemann, T., and Wolf, M. (2010). How to detect and reduce movement artifacts in near-infrared imaging using moving standard deviation and spline interpolation. Physiol. Meas. 31, 649–662. doi: 10.1088/0967-3334/31/5/004

PubMed Abstract | CrossRef Full Text | Google Scholar

Schwerin, S. C., Yao, J., and Dewald, J. P. A. (2011). Using paired pulse TMS to facilitate contralateral and ipsilateral MEPs in upper extremity muscles of chronic hemiparetic stroke patients. J. Neurosci. Methods 195, 151–160. doi: 10.1016/j.jneumeth.2010.11.021

PubMed Abstract | CrossRef Full Text | Google Scholar

Shibasaki, H. (2008). Human brain mapping: hemodynamic response and electrophysiology. Clin. Neurophysiol. 119, 731–743. doi: 10.1016/j.clinph.2007.10.026

PubMed Abstract | CrossRef Full Text | Google Scholar

Takeuchi, N., Chuma, T., Matsuo, Y., Watanabe, I., and Ikoma, K. (2005). Repetitive transcranial magnetic stimulation of contralesional primary motor cortex improves hand function after stroke. Stroke 36, 2681–2686. doi: 10.1161/01.str.0000189658.51972.34

PubMed Abstract | CrossRef Full Text | Google Scholar

Thomson, R. H., Daskalakis, Z. J., and Fitzgerald, P. B. (2011a). A near infra-red spectroscopy study of the effects of pre-frontal single and paired pulse transcranial magnetic stimulation. Clin. Neurophysiol. 122, 378–382. doi: 10.1016/j.clinph.2010.08.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Thomson, R. H., Maller, J. J., Daskalakis, Z. J., and Fitzgerald, P. B. (2011b). Blood oxygenation changes resulting from supra-threshold transcranial magnetic stimulation. Brain Stimul. 4, 165–168. doi: 10.1016/j.brs.2010.10.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Thomson, R. H., Rogasch, N. C., Maller, J. J., Daskalakis, Z. J., and Fitzgerald, P. B. (2012). Intensity dependent repetitive transcranial magnetic stimulation modulation of blood oxygenation. J. Affect. Disord. 136, 1243–1246. doi: 10.1016/j.jad.2011.08.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Valero-Cabré, A., Amengual, J. L., Stengel, C., Pascual-Leone, A., and Coubard, O. A. (2017). Transcranial magnetic stimulation in basic and clinical neuroscience: a comprehensive review of fundamental principles and novel insights. Neurosci. Biobehav. Rev. 83, 381–404. doi: 10.1016/j.neubiorev.2017.10.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Vernieri, F., Maggio, P., Tibuzzi, F., Filippi, M. M., Pasqualetti, P., Melgari, J. M., et al. (2009). High frequency repetitive transcranial magnetic stimulation decreases cerebral vasomotor reactivity. Clin. Neurophysiol. 120, 1188–1194. doi: 10.1016/j.clinph.2009.03.021

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, W. T., Xu, B., and Butman, J. A. (2017). Improved SNR for combined TMS-fMRI: a support device for commercially available body array coil. J. Neurosci. Methods 289, 1–7. doi: 10.1016/j.jneumeth.2017.06.020

PubMed Abstract | CrossRef Full Text | Google Scholar

Wassermann, E. M., and Lisanby, S. H. (2001). Therapeutic application of repetitive transcranial magnetic stimulation: a review. Clin. Neurophysiol. 112, 1367–1377. doi: 10.1016/s1388-2457(01)00585-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Watanabe, T., Hanajima, R., Shirota, Y., Ohminami, S., Tsutsumi, R., Terao, Y., et al. (2014). Bidirectional effects on interhemispheric resting-state functional connectivity induced by excitatory and inhibitory repetitive transcranial magnetic stimulation. Hum. Brain Mapp. 35, 1896–1905. doi: 10.1002/hbm.22300

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: repetitive transcranial magnetic stimulation, near infrared spectroscopy, hemodynamic response, cortical excitability, brain connectivity

Citation: Li R, Potter T, Wang J, Shi Z, Wang C, Yang L, Chan R and Zhang Y (2019) Cortical Hemodynamic Response and Connectivity Modulated by Sub-threshold High-Frequency Repetitive Transcranial Magnetic Stimulation. Front. Hum. Neurosci. 13:90. doi: 10.3389/fnhum.2019.00090

Received: 24 November 2018; Accepted: 25 February 2019;
Published: 19 March 2019.

Edited by:

Zhen Yuan, University of Macau, China

Reviewed by:

Haleh Aghajani, Baylor Scott & White Research Institute (BSWRI), United States
Reza Kazemi, Atieh Clinical Neuroscience Center, Iran
Banghe Zhu, University of Texas Health Science Center at Houston, United States

Copyright © 2019 Li, Potter, Wang, Shi, Wang, Yang, Chan and Zhang. 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: Yingchun Zhang, yzhang94@uh.edu

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.