If you don't remember your password, you can reset it by entering your email address and clicking the Reset Password button. You will then receive an email that contains a secure link for resetting your password
If the address matches a valid account an email will be sent to __email__ with instructions for resetting your password
Heart & Brain Group, Brainclinics Foundation, Nijmegen, The NetherlandsDepartment of Cognitive Neuroscience, Faculty of Psychology and Neuroscience, Maastricht University, Maastricht, The NetherlandsNeurowave, Amsterdam, The Netherlands
Heart & Brain Group, Brainclinics Foundation, Nijmegen, The NetherlandsDepartment of Cognitive Neuroscience, Faculty of Psychology and Neuroscience, Maastricht University, Maastricht, The Netherlands
Heart & Brain Group, Brainclinics Foundation, Nijmegen, The NetherlandsDepartment of Cognitive Neuroscience, Faculty of Psychology and Neuroscience, Maastricht University, Maastricht, The NetherlandsSynaeda Psycho Medisch Centrum, Leeuwarden, The Netherlands
Department of Psychiatry, University of Toronto, Toronto, Ontario, CanadaTemerty Centre for Therapeutic Brain Intervention, Centre for Addiction and Mental Health, Toronto, ON, Canada
Heart & Brain Group, Brainclinics Foundation, Nijmegen, The NetherlandsDepartment of Cognitive Neuroscience, Faculty of Psychology and Neuroscience, Maastricht University, Maastricht, The Netherlands
Neuro-cardiac-guided transcranial magnetic stimulation (NCG-TMS) employs repetitive TMS (rTMS)-induced heart rate deceleration to confirm activation of the frontal-vagal pathway. Here we test a novel NCG-TMS method that utilizes Heart-Brain-Coupling (HBC) in order to quantify rTMS-induced entrainment of the inter-beat-interval as a function of TMS cycle-time. Since prior NCG-TMS studies indicated no association between motor and frontal excitability threshold, we also introduce the approach of using HBC to establish individualized frontal excitability thresholds for optimally dosing frontal TMS.
Methods
In study 1A and 1B we validated intermittent theta-burst stimulation (iTBS)-induced HBC (2s iTBS-on; 8s off: HBC=0.1Hz) in fifteen (1A) and twenty-two (1B) patients with major depressive disorder from two double-blind placebo-controlled studies. In study 2, HBC was measured in ten healthy subjects during the 10Hz "dash”-protocol (5s 10Hz-on; 11s off: HBC=0.0625Hz) applied with fifteen increasing intensities to four evidence-based TMS locations.
Results
Using blinded electrocardiogram-based HBC-analysis, we successfully identified sham from real-iTBS sessions (accuracy Study 1A=83%, Study 1B=89.5%) and found a significantly stronger HBC at 0.1Hz in active compared to sham iTBS (d=1.37) (Study 1A). In study 2, clear dose-dependent entrainment (p=.002) was observed at 0.0625Hz in a site-specific manner.
Conclusions
We demonstrated rTMS-induced HBC as a function of TMS cycle-time for two commonly used clinical protocols (iTBS and 10Hz-dash). These preliminary results supported individual site-specificity and dose-response effects, indicating that this is a potentially valuable method for clinical rTMS site-stratification and frontal-thresholding. Further research should control for TMS side-effects, such as pain of stimulation, to confirm these findings.
Repetitive transcranial magnetic stimulation (rTMS) is a non-invasive neuromodulation technique, which is increasingly used as intervention for the treatment of refractory major depressive disorder (MDD) (
). rTMS for MDD is often targeted to the dorsolateral prefrontal cortex (DLPFC), and clinical response is thought to be mediated by network connectivity between the DLPFC and the subgenual anterior cingulate cortex (sgACC)(
). There is great potential to optimize rTMS parameters and move towards an individualized approach of rTMS therapy, including optimizing target engagement and parameters such as rTMS pattern, intensity, and frequency of stimulation (
It is well known that MDD is associated with an increased risk of cardiovascular disease, and heart rate (HR) is often dysregulated in MDD, quantified by a higher HR and lower HR variability (HRV) (
). Relatedly, we recently proposed the frontal-vagal network theory for MDD, stating that major hubs such as the DLPFC, sgACC and vagus nerve (VN) share overlap with the networks involved in autonomic control as well as in MDD(5). The neuro-anatomical framework of the heart-brain connection in MDD is provided in (
), called neuro-cardiac guided TMS (NCG-TMS). In NCG-TMS, HR-deceleration in response to 10Hz-rTMS or intermittent theta burst stimulation (iTBS) is used to confirm activation of the frontal-vagal network in a site-specific manner, with HR-deceleration at sites often used in TMS (e.g. F3, FC3) and HR-acceleration in control sites such as the motor (C3, C4) or parietal (Pz) cortex (
), and therefore holds promise as a possible biomarker of response and target-engagement method allowing for determination of the best prefrontal rTMS target in MDD treatment(
). As an intermediate step, we recently proposed that NCG-TMS could be employed as an enhanced rTMS stratification technique to select between one of two evidence-based DLPFC-sites: the Beam-cluster (Beam F3 and F4) and 5cm-rule (5CM) method (
), but it also resulted in more side effects, such as light-headedness, emotional reactions, and painfulness. For this reason, we set-out to improve the current 10Hz NCG-TMS method. Recently, the FDA approved the 10Hz “dash”-protocol(
), where the inter-train interval (ITI) was shortened to 11s, allowing for more rapid delivery of the stimulation trains during the session.
Given that effects of TMS-induced HR-deceleration are nearly immediate, since the vagus nerve is myelinated (whereas the sympathetic inputs to the heart are unmyelinated; for review see (
)), we hypothesized that specific TMS-protocol parameters could induce entrainment of the cardiac rhythm as a function of TMS cycle-time (Figure 1), as described in (
) for rTMS and EEG alpha oscillations. If entrained, the 2s (iTBS) or 5s (10Hz) stimulation would decelerate HR during and immediately after each train, followed by a recovery of the HR during the ITI. Over the session with repeated trains, an entrainment of the inter-beat-interval at a specific frequency could be observed. The heart-brain directionality of TMS would thus be top-down from the brain on the HR(13). This hypothesis is visualized in Figure 1 for the 10Hz dash-protocol (top) and iTBS (bottom). This novel ‘Heart-Brain-Coupling’ (HBC) marker is proposed, since HBC is independent of ‘baseline’ HR state and other influences on HR, such as respiratory sinus arrhythmia (RSA), making it a more reliable and stable biomarker. Furthermore, RSA is generally faster (0.15–0.4 Hz) than these TMS cycle-times (
Figure 1rTMS-induced Heart-Brain-Coupling. This figure visualizes rTMS-induced entrainment of the inter-beat-interval as a function of rTMS cycle-time. The cycle-time of the 10Hz dash-protocol comprises one rTMS train every 16 seconds (5s stimulation on, 11s stimulation off). This results in a specific entrainment frequency of 1/16 = 0.0625Hz. For iTBS, the cycle-time is 10 seconds (2s on and 8s off), resulting in an entrainment frequency of 1/10 = 0.1Hz. The neuro-anatomical framework of the heart-brain connection illustrated on the left is provided in (
We thus set-out to test the above predictions for iTBS and the 10Hz dash-protocol in three independent studies. Finally, as prior NCG-TMS studies indicate that the frontal excitability threshold (FT) was not associated with the motor threshold (MT, which is used to ‘dose’ prefrontal rTMS treatment) (
), we also tested if this method could be utilized to establish an FT. Two previous studies investigated the dose-response relationship of stimulation intensity of DLPFC-iTBS and clinical response. In healthy individuals, the largest neurophysiological changes were found at 75%MT, compared to 50% and 100%MT (
Subthreshold stimulation intensity is associated with greater clinical efficacy of intermittent theta-burst stimulation priming for Major Depressive Disorder.
). These studies imply that higher TMS intensities could lead to attenuated effects and research is needed to define the individual sweet spot of stimulation intensity.
METHODS AND MATERIALS
Study 1: iTBS induced Heart-Brain-Coupling
To test the specific entrainment of iTBS on the cardiac rhythm, we used data from two previous iTBS studies: Baseline ECG data collected at a Canadian health science center (author FVR) from a double-blind placebo controlled iTBS trial (CARTBIND; ClinicalTrials.gov Identifier: NCT02729792), reported in more detail in Blumberger et al. (
A randomized sham controlled comparison of once vs twice-daily intermittent theta burst stimulation in depression: A canadian rTMS treatment and biomarker network in depression (CARTBIND) study.
) (Study 1A), and ECG data from a double-blind randomized controlled trial performed at the Department of Psychiatry at Stanford University (SNT: ClinicalTrials.gov Identifier: NCT03068715), reported in more detail in Cole et al.(
) (Study 1B). In summary, the methods comprised the following:
Participants
Fifteen (Study 1A) and twenty-two (Study 1B) MDD patients were included, and all participants provided written informed consent.
iTBS device and protocol
iTBS-treatment localization was MRI-based, and rTMS was applied with a MagPro X100 system (Magventure, Denmark) equipped with a B70 figure-of-eight coil (Study 1A) and a double-sided Cool-B65 A/P coil (Study 1B).
Study 1A: All subjects received two sessions of both sham and active iTBS. The active condition comprised stimulation at the DLPFC location with triplet 50Hz bursts, repeated at 5Hz; 2s on and 8s off; 600 pulses per session. A sham (internally shielded) coil without active electrical stimulation was positioned over the vertex during the sham condition. Subjects were randomized to one of two treatment arms (arm A: sham-active, 54 min. pause, active-sham, or arm B: sham-sham, 54 min. pause, active-active). ECG was acquired simultaneous with iTBS stimulation using a Biopac MP 150 system (Biopac Systems Inc., Goleta, CA, USA) comprising modular hardware and “AcqKnowledge” software.
Study 1B: Subjects received either active or sham “Stanford Neuromodulation Therapy” (SNT): A high-dose iTBS protocol over the DLPFC, which consisted of ten sessions (18,000 pulses) per day, on 5 consecutive days. Both sham and active iTBS were targeted at the left DLPFC. Participants and study staff were blinded to treatment assignments. All iTBS sessions utilized the same stimulation coil, with no indication of active or sham orientation. During sham, participants wore noise-cancelling earphones connected to a sham noise generator to simulate the stimulation noise pattern. Additionally, lidocaine was applied to the stimulation site to reduce sensation. Spontaneous side effects were recorded daily. ECG was acquired simultaneous with iTBS using NCG-engage (neuroCare), only the first recording for each participant was used for analysis.
Coil positioning and MT
Study 1A
The target location was specified by reverse co-registration from a stereotaxic coordinate on the standard Montreal neurological Institute (MNI-152) template brain, onto each individual anatomical MRI. MNI coordinates for left DLPFC were [x-38 yþ44 zþ26], drawn from a study identifying this site as optimal based on clinical outcomes and resting-state functional connectivity (
During active and sham stimulation, the coil was positioned over the DLPFC using MRI-guided neuronavigation. The stimulation intensity at 90% MT was adjusted for depth of the identified MRI-target. For safety, stimulation intensity never exceeded 120% MT.
) using a frequency range of 0.02 to 0.18 Hz in steps of 5e-4Hz, with three cycles yielding results in high time resolution (Figure 2B) and ten cycles for high frequency resolution (Figure 2C). Before TFR-analysis, data was padded with the first block of rTMS and rest period of data at the beginning, and the last block of data at the end to allow for analysis at these low frequencies. To compute HBC, the mean power (μV2) at 0.1Hz was computed for each block and averaged.
Figure 2HBC report of active and sham iTBS condition. This overview visualizes the effect of one active and one sham iTBS session on HR (panel A) over 3-min time (x-axis), the 0.1Hz high time resolution (panel B) and the 0.1Hz high frequency resolution (panel C). The vertical lines represent the start of a stimulation train. Panel A visualizes that during stimulation, the HR consistently decreases, and normalizes during the ITI, in line with the hypothesized effects in Figure 1. This rTMS-induced cardiac rhythm of 0.1Hz (green line) is also clearly visible as increased power (red) in the high time and high frequency plots. HBC=Heart-brain-coupling, iTBS=intermittent theta-burst stimulation, HR=heart rate.
For three subjects, insufficient ECG data was available, resulting in n=12 included in the analysis. Researcher HvD - blinded to group assignment and stimulation details - processed the ECG data based on the above hypothesis of 0.1Hz entrainment in RR-signal and made predictions if ECG data belonged to sham or real stimulation within each subject. These predictions were examined after data-lock by an unblinded researcher (MA). Next, a repeated measures ANOVA was conducted to test the effect of active and sham iTBS on the HBC marker (0.1Hz).
Study 1B
Three participants were excluded for low quality ECG, resulting in n=19. The ECG was analysed by blinded researcher HvD, who used the same method as in Study 1A to classify their active/sham status.
Study 2: DASH (10Hz) NCG-TMS
Participants
Ten healthy participants (seven female, aged 25-46 [mean 34.0] years) were included in this study after providing written informed consent and meeting safety criteria for rTMS. Exclusion criteria for both studies were (
Figure 3Visualization of the intensity-sweep in study 2. This figure shows the intensity-sweep that was applied in this pilot study 2, using fifteen increasing intensities with 2%MSO-steps. Step 15 corresponds to 120%MT. This rTMS-protocol was applied on four sites: the Beam and 5CM location on both hemispheres. %MSO=% of machine outpout, MT=motor threshold. rTMS=repetitive transcranial magnetic stimulation, 5CM location=the location 5cm anterior to the scalp position for optimal activation of the first dorsal interosseus muscle in a para-sagittal line.
The rTMS-coil was positioned over the Beam and 5CM location at a 45° angle relative to the parasagittal plane (the coil handle pointing posteriorly), which seems to be the optimal angle to stimulate frontal areas (
). If the distance between these two locations was smaller than one cm, the Beam and 5CM location on that hemisphere were taken as one location (n=1). The 5CM location is defined as the location 5cm anterior to the scalp position for optimal activation of the first dorsal interosseus muscle in a para-sagittal line. Beam locations were defined using the Beam-F3 algorithm and software (
). MT intensity was defined as the lowest stimulation intensity that in four trials, induced at least two visible twitches in the contralateral hand.
Analysis
Data were analyzed using the custom-built analysis package described above, now including the determination of the RR intervals. For this, ECG was bandpass-filtered with a bidirectional 4th order Butterworth filter between 5 and 49Hz, after which signal.find_peaks (Scipy (
)) was used to detect R-peaks. RR intervals were then corrected for ectopic beats and subsequent analysis was performed as described above, in this case for 0.0625Hz. HBC was quantified across the intensity-sweep and four sites and visualized in a custom HBC-marker report (Figure 4 and supplement). Primary outcome measure was the HBC-marker for the 10Hz dash-protocol (0.0625Hz). Using the HBC reports, the ‘best’ location was determined for every subject, based on the location with the highest average oscillatory power (μV2) at 0.0625Hz (i.e. over all intensities) (supplement). Data for all target locations and for the best location was subsequently analyzed in SPSS with a repeated measures ANOVA, using within subject factors Intensity (0-15) and Location (BF3, BF4 and 5CM left and right).
Figure 4HBC report example. This figure shows HBC-power of 0.0625Hz over time during the intensity-sweep for one location in two different individuals (Panel A1 and A2). Step 0-15 indicates the intensity-sweep of increasing stimulator output, where step 0 is no stimulation. HBC is seen later in time for individual A1 compared to individual A2, suggesting individual differences in low and high threshold response. Panel B shows the summary of HBC-power of 0.0625Hz during the intensity-sweep for all stimulation locations for individual A1. Here, BF4 would be the selected best location given the highest average oscillatory power at 0.0625 Hz (“pow”) at that location. HBC=Heart-brain-coupling, BF4=Beam F4 location.
The blinded classifications on sham and active iTBS conditions from the ECG of each subject were highly accurate with 10/12 predicted correctly by the blinded researcher (accuracy=83%), suggesting successful unblinding of stimulation only by inspection of the ECG-based HBC-marker.
Two examples of HBC-output are seen in Figure 2 for a representative example of sham and real iTBS in the same subject. This overview visualizes the effect of one 3-min iTBS session on mean-HR (panel A) and 0.1Hz HBC, visualized using a high time resolution (panel B) and high frequency resolution (panel C). Vertical lines represent the start of each stimulation train. During active iTBS, HR decelerates during every stimulation train (2s) and normalizes afterwards (8s), until the next stimulation train, resulting in an rTMS-entrained HR oscillation at 0.1Hz (horizontal green line) also visible as increased power (red) in the high time and high frequency plots. This is not seen during sham stimulation.
Repeated measures ANOVA yielded a main effect of Stimulation (F(1,10)=16.318; p=.002, d=1.37) and a main effect of Time (F(1,10)=11.767, p=.006; d=0.54), but no Stimulation X Time interaction (p=.318). Median and mean self-rated pain scores during active stimulation, on a scale from 1 [no pain] to 10 [intolerable pain], were higher (4.0 and 4.3, IQR=2.7-6.0, SD=2.0) than during sham stimulation (1.0 and 1.1, IQR=1.0-1.1, SD=0.1)
Study 1B
The generated HBC values predicted active/sham conditions with 89.5% accuracy. These results show that iTBS stimulation can be clearly distinguished from sham iTBS, based on the ECG-based HBC-marker only, also when sham and active stimulation are targeted at the exact same location. Of measured spontaneous side effects, there was only a higher incidence of headache in the active iTBS group compared to sham (Fisher’s exact test, p<0.06). Participants did not guess their treatment allocation beyond chance in both groups.
Study 2
The repeated measures ANOVA on data of all locations showed a significant effect of Intensity (F(15,135)=3.540, p<.001) on the HBC-marker, and no effect of Location (p=0.997) or Location x Intensity interaction (p=0.955). The repeated measures ANOVA conducted on data for the best location per individual, demonstrated a significant effect of Intensity (F(15,135)=3.718, p<.001). The HBC-peak for the best location, defined as the stimulation intensity at which the HBC reached the highest oscillatory power at 0.0625 Hz, reaches its maximum at intensity 8, with a large effect size (d = 2.40) (Figure 5). The data in Figure 5 is derived from the individual HBC reports (supplement), which were normalized to a 0-1 scale. The Beam F3 and F4 were the individual best location for two subjects each, the 5CM left and right for three subjects each. We investigated whether order of stimulation had an impact on the HBC effect at the different rTMS sites. No order-effect was found, implying that HBC effects were independent of which site was stimulated first.
Figure 5Dose-response effect of increasing stimulation intensities on HBC. Dose-response effect of increasing intensities of machine output on the HBC-marker in all locations (
). The HBC-effect for the best location has a higher maximum than the separate locations. The data in Figure 5 is derived from the individual HBC reports (supplement A), which were normalized to a 0-1 scale. HBC=Heart-brain-coupling.
The results indicate a clear dose-response effect of rTMS-induced HBC, with no difference between sites at the group level. This is as expected, as we are investigating a method to individualize rTMS targeting, based on the hypothesis that large inter-individual differences exist. When the best site was chosen for each individual participant, a strong dose-response effect was observed. The dose-response relationship follows an inverted U-curve. More importantly, the data visualized in Figure 5 suggest large inter-individual differences in FT, also confirmed by individual trajectories as visualized in Figure 4 and 6. The peak effect of stimulation on the HBC-marker arises at different stimulation strengths for individual subjects (Figure 6), with the greatest difference between subject 10 (intensity 7) and 1 (intensity 15). These HBC-peaks correspond with a relative %MT of 81% for subject 10 and 120% for subject 1. Importantly, for most of the individual subjects, the greatest effect arises sooner than step 15, which is currently the standard stimulation strength in clinical practice applied to frontal rTMS.
Figure 6Individual trajectories of the HBC-marker over time. Dose-response effect of increasing intensities of machine output on the HBC-marker in all subjects. HBC=Heart-brain-coupling.
This study is the first to test and validate a novel NCG-TMS method utilizing Heart-Brain-Coupling, that quantifies rTMS-induced entrainment of the heart inter-beat-interval as a function of TMS cycle-time (Figure 1). The proposed model makes clear predictions about two commonly used clinical rTMS protocols used for the treatment of depression (
Effectiveness of theta burst versus high-frequency repetitive transcranial magnetic stimulation in patients with depression (THREE-D): a randomised non-inferiority trial.
): iTBS with a cycle-time of 10s and 10Hz-dash with a cycle-time of 16s. First, we validated this model for iTBS, in which the cycle-time resulted in a 0.1Hz entrainment, where we could successfully differentiate between active and sham iTBS with a large effect size and unblind two datasets with high accuracy. The effect of iTBS on HBC did not attenuate over time in study 1A (opposed to the effects of HR-deceleration in the older NCG-TMS 1.0 (
)), which could be caused by neuroplastic effects of the first iTBS exposure. We further validated this finding for a second rTMS protocol (10Hz-dash), where stimulation resulted in specific 0.0625Hz HR-entrainment.
In addition to target engagement results described above, study 2 also demonstrated clear dose-response effects on the HBC-marker with large inter-individual variability in FT. Prior studies already demonstrated that NCG-TMS strength was not correlated to %MT but to %MSO, suggesting %MT is a poor proxy for the FT (
). In our limited sample, the greatest difference in HBC-peak was found between two subjects with a relative %MT of 81%MT and 120%MT (Figure 6). These different TMS-HBC distributions could reflect a combination of TMS intensity and possibly another factor, such as pain, target engagement, or individual differences in susceptibility to TMS. In larger samples it is likely that larger differences will be found, necessitating studies to investigate this method further as an FT-technique. The non-linear and non-sigmoidal dose-response effects suggest possible ‘over-stimulation’, which is in line with a recent case of TMS-induced syncope (
Subthreshold stimulation intensity is associated with greater clinical efficacy of intermittent theta-burst stimulation priming for Major Depressive Disorder.
). Compared to those studies, we used more and smaller incremental steps in the current study and investigated HBC on the individual level, enabling the determination of an individual fine-grained dose-response curve.
Some potential limitations should be acknowledged. In all three studies, small sample sizes were included, decreasing the power of the studies. Study 1A and 1B are hard to compare, as they differ in the number of TMS pulses (1A: 1200 pulses in 2 sessions; 1B: 90000 pulses in 50 sessions) and the design (1A: within-subject; 1B: between-subject design). Also, in study 1B and study 2 subjective unpleasantness of stimulation was not measured, study 1A lacked a good control condition (vertex vs. DLPFC) and active condition was more painful than sham, and study 1B was not a within-subject design (participants could not compare sham and active condition directly). Therefore, distinctive scalp sensations might have contributed to differences between active and sham condition. The sensory effects of TMS, such as the pain and discomfort associated with stimulation, and its influence on heart rate and HBC should be controlled for in future TMS studies, by measuring pain and discomfort during stimulation and use that as a covariate, or by conducting stimulation inside an MRI to demonstrate correlation between BOLD changes and HR changes. This is important, as it is known that frontal TMS is especially uncomfortable (
). This was not sufficiently controlled for in the current study, necessitating future research to investigate the potential relevance of these results. Finally, it should be acknowledged that much is still unknown about the FT, hence it is unclear whether the individualized FT has any associations or implications for clinical outcomes of rTMS treatment.
Despite of the limitations, there are at least three good reasons to believe that the TMS protocol has direct effects on the frontal-vagal system and that non-specific effects of rTMS, such as somatosensory stimulation, are unlikely to have affected the results: (
) in study 2, the HBC-effect demonstrated clear site-specificity. For non-specific effects, a similar effect on the HR would have been expected for all four stimulated locations. In addition, (
), rTMS was applied in anesthetized rats and they found HR-decelerations after active stimulation but not sham, emphasizing active signaling, rather than just the effect of pain.
Our findings suggest that the use of rTMS-induced HBC might be a valuable method for rTMS target engagement of the frontal-vagal pathway. Given this novel approach deviates substantially from the earlier published NCG-TMS studies, we propose to refer to this method as NCG-TMS 2.0 or Heart-Brain-Coupling. In previous studies, it was shown that the effects of rTMS on HR are similar between healthy subjects and MDD patients (
), which demonstrates the possible clinical implications of this novel method for target stratification in MDD patients. In line with the suggestions by (
), this method could be specifically used to stratify MDD patients between one of two evidence-based stimulation clusters: The Beam and 5CM cluster. Large effectiveness studies have demonstrated comparable response (47-58%) and remission (29-37%) rates for these clusters (
Effectiveness of theta burst versus high-frequency repetitive transcranial magnetic stimulation in patients with depression (THREE-D): a randomised non-inferiority trial.
Transcranial magnetic stimulation (TMS) for major depression: A multisite, naturalistic, observational study of acute treatment outcomes in clinical practice.
), thus NCG TMS 2.0 could be employed to effectively address inter-individual differences between MDD patients and stratify patients to their optimal rTMS target. This could enhance response and remission rates on the individual level, while adhering to validated rTMS-protocols, as both targets are widely implemented in clinical practice.
Future studies should focus on replicating these findings in larger samples and investigate whether NCG TMS 2.0 should be used for rTMS target stratification and probing the FT in MDD patients to improve treatment outcome. Further research should focus on optimization of other rTMS parameters as well, such as frequency of stimulation. Also, the potential of NCG-TMS 2.0 could be investigated further for TMS in other (cognitive) applications. The current study will hopefully aid to individualize rTMS treatment and might have important implications for the field of stratified psychiatry. NCG TMS may become for the prefrontal cortex, what the thumb-twitch is for the motor system.
Subthreshold stimulation intensity is associated with greater clinical efficacy of intermittent theta-burst stimulation priming for Major Depressive Disorder.
A randomized sham controlled comparison of once vs twice-daily intermittent theta burst stimulation in depression: A canadian rTMS treatment and biomarker network in depression (CARTBIND) study.
Effectiveness of theta burst versus high-frequency repetitive transcranial magnetic stimulation in patients with depression (THREE-D): a randomised non-inferiority trial.
Transcranial magnetic stimulation (TMS) for major depression: A multisite, naturalistic, observational study of acute treatment outcomes in clinical practice.
FVR receives research support from CIHR, Brain Canada, Michael Smith Foundation for Health Research, Vancouver Coastal Health Research Institute, and Weston Brain Institute for investigator-initiated research. Philanthropic support from Seedlings Foundation. In-kind equipment support for an investigator-initiated trial from MagVenture. He has received an honorarium for participation in an advisory board for Allergan.
DMB receives research support from the Canadian Institutes of Health Research (CIHR), National Institutes of Health – US (NIH), Brain Canada Foundation and the Temerty Family through the CAMH Foundation and the Campbell Family Research Institute. He received research support and in-kind equipment support for an investigator-initiated study from Brainsway Ltd. and he was the site principal investigator for three sponsor-initiated studies for Brainsway Ltd. He received in-kind equipment support from Magventure for investigator-initiated studies. He received medication supplies for an investigator-initiated trial from Indivior. He has participated in an advisory board for Janssen and for Welcony Inc.
JD reports research grants from CIHR, the National Institute for Mental Health, Brain Canada, the Canadian Biomarker Integration Network in Depression, the Ontario Brain Institute, the Klarman Family Foundation, the Arrell Family Foundation, and the Edgestone Foundation; reports travel stipends from Lundbeck and ANT Neuro, reports in-kind equipment support for investigator-initiated trials from MagVenture, and is an advisor for BrainCheck, NeuroStim TMS and Salience Neuro Health.
NW is a named inventor on Stanford-owned intellectual property relating to accelerated TMS pulse pattern sequences and neuroimaging-based TMS targeting; he has served on scientific advisory boards for Otsuka, NeuraWell, Nooma, and Halo Neuroscience; and he has equity/stock options in Magnus Medical, NeuraWell, and Nooma
ATS is Chief Scientific Advisor of PlatoScience and Alphasys, CEO of Neurowear Medical B.V., received equipment support from MagVenture, Deymed, and MagStim Company, and is Scientific Director of the International Clinical TMS Certification Course (www.tmscourse.eu).
MA is unpaid chairman of the Brainclinics Foundation, holds equity/stock in neurocare, serves as consultant to neurocare, and is named inventor on neurocare owned patent and intellectual property related to neuro-cardiac-guided TMS, but receives no royalties; Research Institute Brainclinics received research and consultancy support from neurocare, Brainify, UrgoTech and equipment support from neuroconn and Deymed.
All other authors report no biomedical financial interests or potential conflicts of interest.