Individualized rTMS Treatment for Depression using an fMRI-Based Targeting Method

Summary

The present protocol describes the application of repetitive transcranial magnetic stimulation (rTMS), where a subregion of the dorsolateral prefrontal cortex (DLPFC) with the strongest functional anticorrelation with the subgenual anterior cingulate cortex (sgACC) was located as the stimulation target under the assistance of a fMRI-based neuronavigation system.

Abstract

To achieve greater clinical efficacy, a revolution in treatment for major depressive disorder (MDD) is highly anticipated. Repetitive transcranial magnetic stimulation (rTMS) is a non-invasive and safe neuromodulation technique that immediately changes brain activity. Despite its wide application in the treatment for MDD, the treatment response remains different among individuals, which may be attributable to the inaccurate positioning of the stimulation target. Our study aims to examine whether the functional magnetic resonance imaging (fMRI)-assisted positioning improves the efficacy of rTMS in treating depression. We intend to identify and stimulate the subregion of dorsolateral prefrontal cortex (DLPFC) in MDD with strongest anti-correlation with the subgenual anterior cingulate cortex (sgACC), and to conduct a comparative investigation of this novel method and the traditional 5-cm rule. To achieve more precise stimulation, both methods were applied under the guidance of neuronavigation system. We expected that the TMS treatment with individualized positioning based on resting state functional connectivity may show better clinical efficacy than the 5-cm method.

Introduction

Major depressive disorder (MDD) is characterized by significant and persistent depression, and in more severe cases, patients can encounter hallucinations and/or delusions^1,2. Compared with the general population, the risk of suicide among MDD patients is approximately 20 times higher³. While medication is currently the most used treatment for MDD, 30% - 50% of the patients lack adequate response to antidepressants⁴. For the responders, the symptom improvement tends to appear after a relatively long latent period and is accompanied by side effects. Psychotherapy, although effective for some patients, is costly and time-consuming. A safer and more effective treatment for MDD is therefore urgently required.

Repetitive transcranial magnetic stimulation (rTMS)is a non-invasive and safe technique and has been approved for the treatment of various mental disorders^5,6,7. Although its therapeutic mechanism remains unclear, rTMS was speculated to work by regulating the activity of the stimulated brain regions and the neural plasticity^8,9,10, thus normalizing specific functional networks^10,11,12. rTMS also causes network effect, which evokes changes in remote brain areas through connection pathways, leading to an amplified therapeutic effect¹³. Although rTMS changes brain activity immediately and robustly, its response rate in the treatment of MDD is only about 18%¹⁴. The main reason may be the inaccurate location of stimulation targets¹⁵.

The subgenual anterior cingulate cortex (sgACC) is mainly responsible for emotional processing and plays a role in regulating the response to stressful events, emotional response to internal and external stimuli, and emotional expression^16,17,18. This subregion of ACC shares substantial structural and functional connectivity with the cerebral cortex and the limbic system^19,20. Interestingly, studies have shown that the post-stimulation activity of this area is closely related to the clinical efficacy of TMS. For instance, the blood flow of sgACC decreased after a course of TMS targeted on the right dorsolateral prefrontal cortex (DLPFC), which was associated with the alleviation of depressive symptoms²¹. Vink et al.⁸ found that stimulation targeted on DLPFC was propagated to sgACC, and suggested that sgACC activity can be a biomarker of the treatment response of TMS. According to previous researches, Fox and colleagues²² proposed that targeting on a subregion of DLPFC that shows strongest functional anti-connectivity with sgACC (MNI coordinate: 6, 16, -10) enhances the antidepressant effect. Here, we demonstrate a study protocol aimed to examine this hypothesis.

Protocol

Inform all participants about the study and ask them to sign the informed consent form prior to the start of the study. The present protocol was approved by the Research Ethics Committee of the Affiliated Brain Hospital of Guangzhou Medical University.

NOTE: In this double-blind study, patients with depression were randomly divided into two groups. In the experimental group, stimulation targets are located by the DLPFC-sgACC-based individualized location method (Please see 3.3 for detailed description). The targets of the control group are obtained by the average 5-cm method (i.e. (-41, 16, 54))²².

1. Participants' selection

1. Recruit patients with a diagnosis of MDD as confirmed by an expert psychiatrist.
   NOTE: Confirm the diagnosis with the standardized MINI-International Neuropsychiatric Interview (M.I.N.I.)²³ . The total score of Montgomery-Asberg Depression Rating Scale (MADRS)²⁴ should be no less than 22.
2. Exclude patients who meet the exclusion criteria: (1) serious physical diseases such as malignant tumor, acute heart failure, multiple organ failure, or severe neurological conditions including but not limited to epilepsy, stroke, encephalitis, brain trauma; (2) comorbidity of other mental illness, or a history of substance use disorder; (3) having metallic implants, especially in brain or heart; (4) women during pregnancy or lactation; (5) had suicidal behaviors or attempted suicide in the past six months; and (6) a diagnosis of bipolar depression or psychotic depression.
   ​NOTE: Recruit at least 36 subjects for each group to ensure statistical power. A balanced demographic profile between the two groups is recommended.

2. Preparation of Magnetic Resonance Imaging (MRI) and TMS

1. Obtain fMRI images by a 3T MRI scanner before performing TMS.
   1. Reconfirm that the patient has no contraindications before an MRI scanning. Instruct the patient to try to lie still and think of nothing during the scan.
   2. Conduct a resting-state fMRI (rs-fMRI) scan using the FFE-EPI sequence with the following parameters: TR/TE = 2000/30 ms, FA = 90°, field-of-view = 220 x 220 x 256 mm³, matrix = 64 x 64, voxel size = 3.44 x 3.44 x 4 mm³, gap = 0.6 mm, number of signal averages = 1, volumes = 240, number of slices = 33.
   3. Conduct a structural MRI scan using the sagittal T1 weighted 3D turbo field echo (T1W 3D TFE) sequence with the following parameters: field of view = 256 x 256 mm², TR/TE = 8.2/3.8 ms, view matrix= 256 x 256, slice thickness= 1 mm.
2. Set TMS parameters.
   NOTE: The protocol of TMS in our study is the intermittent theta-burst stimulation (iTBS). A daily treatment session includes 60 cycles of 10 bursts of 3 pulses at 50 Hz delivered at 100% RMT in 2-s trains, with an interval of 8 s. The whole treatment consists of 10 sessions carried out on weekdays of two consecutive weeks.

3. Treatment (Figure 1)

1. Conduct MRI scans and clinical assessments of symptoms and cognitive performance one day before the treatment.
2. Assign the patient randomly to one of the two groups, after scanning.
3. For the experimental group, identify the subregion of DLPFC that shows strongest functional anti-connectivity with sgACC. For the control group, simply locate the target in the standard-space using the average 5-cm method, then convert it to the individual space coordinates.
   1. rs-fMRI data preprocessing
      1. Preprocess the rs-fMRI data using an MRI analysis software: (a) remove the first 10 volumes; (b) conduct  the slice timing correction; (c) correct the head motion; (d) co-register EPI images to T1 images; (e) perform segmentation; (f) perform normalization using T1 images; (g) smooth the normalized images with a 6-mm Gaussian kernel of full-width half maximum (FWHM); (h) band-pass filter (0.01 - 0.08 Hz); and (i) perform nuisance regression (head motion effects, linear trends, white matter, cerebrospinal fluid, and global mean time course).
   2. Functional connectivity (FC) of the sgACC
      1. Select the sgACC (MNI coordinate: 6, 16, −10; Fox et al.) as the region of interest (ROI)²⁵ with a 10-mm radius.
      2. Remove the white matter and cerebrospinal fluid in the ROI based on the Harvard-Oxford cortical atlas (http://www.cma.mgh.harvard.edu/), using a gray matter probability threshold of 0.25.
      3. Extract the average time course of the ROI.
      4. To generate FC map, compute Pearson's correlation coefficients between the ROI (sgACC) and DLPFC in a voxel-wise manner. Normalize each correlation coefficient using the Fisher's r-to-z Transformation.
         NOTE: The DLPFC mask is a combination of 20mm radius spheres centered along the left hemisphere at BA9 (x=-36, y=39, z=43), BA46 (x=-44, y=40, z=29), the 5-cm approach site (x=-41, y=16, z=54), and the F3 Beam group-average stimulation site (x=-39, y=26, z=49)²⁶.
      5. According to the FC map, identify the peak coordinate in DLPFC that has the largest Pearson's anti-correlation coefficient with sgACC. This is the subregion of DLPFC with the strongest negative FC with sgACC, which will be later targeted in the TMS treatment for the experimental group.
4. Determine the resting motor threshold (RMT) for each subject and record the hotspot.
   1. Instruct the patient to sit back and relax, then put two recording electrodes on the thenar of the right hand and a reference electrode on the bony part of the wrist.
   2. Stimulate the motor hotspot with 10 consecutive stimulations with different intensities; in the meanwhile, record the times of thenar muscle contraction.
   3. Identify the minimum TMS intensity at which a motor evoked potential (MEP) ≥ 50 µV is recorded at least 5 times. Define it as the patient's RMT.
5. Assess the severity of depression using clinical scales as described in Clinical Data Collection.
6. Perform TMS treatment twice a day for 10 days.
   NOTE: For a subject who did not receive treatments as planned, perform additional stimulations after the end of treatment course as needed. However, any subject who misses the treatment for more than four consecutive days should be excluded.
   1. Create a new patient entry.
      1. Select the option Create New Patient. Input the patient's ID number or name in the textbox.
   2. Overlay the structural MRI images onto the navigation system.
      1. Select Import patient MRI, and then import the structural image of the patient and select the image type.
   3. Create individual head model and define the stimulation target.
      1. Press the button Specify MRI Fiducials.
      2. Place the crosshair on these spots in the MRI image: (1) Fiducial Markers: nasion, both left and right tragi; (2) AC-PC Markers for Talairach: anterior commissure, posterior commissure, inter-hemispheric point; (3) Talairach Markers: anterior point, posterior point, superior point, inferior point, left point, and right point.
         NOTE: The "Talairach Markers" mark the borders of the brain.
      3. Press Create Head Model. Select Manual Brain Segmentation and adjust the threshold of the scalp, the lower brain, and the upper brain.
      4. Click Define Target to proceed.
      5. Select Target Marker page. Click … to input the coordinate of the treatment target as identified in step 3.3, and then press Go to. Press Add Marker to name the point.
         NOTE: The coordinates of the control group will be (-41, 16, 54).
   4. Coil calibration
      1. Click Proceed to Neuronavigation. In the textbox, select the right type of tools to be used in the treatment. Ensure all the tools of reference are in the view of the infrared camera.
      2. Press Validate Coil. Put the tip of pointer on the marked coil point. Press Validate (or the green button on the remote control) when the indicator of each tool turns green.
   5. Select patient and target.
      1. Select patient's name or ID on the Select Patient page. Click Select Targets on the next page.
      2. Choose Read target markers to browse the file of targets. Import the file and select the target as in step 3.6.3.5.
   6. Define the coordinate system.
      1. Click Define Coordinate System. Put a headband with a reference tool on the patient. Make sure the coil tracker and the reference tool are in the vision of the navigation system.
      2. Place pointer's tip on the nasion and both tragi in turn. Press the green button on the remote control each time when the indicator of a marker turns green.
   7. Head shape generation
      1. Continuously move the pointer's tip on the top of the head. Press the button on the remote control (or Fit) to continue.
         NOTE: One can press Pause to stop the process, and resume by pressing the Start button again, once the pointer has been placed correctly.
   8. Neuronavigation and stimulation
      1. Press Neuronavigation. On the Active Coil page, set the stimulation intensity to 100% RMT. Choose Stimulate at targets to see the target on the head model online.
      2. When the coil matches the target crosshair, that is, when the indicator text turns green, stimulate.
   9. Prepare and start the treatment from step 3.6.4 directly if the patient's entry had been created before.
7. Conduct follow-up assessments on Day 1, Day 28, and Day 56 after the whole treatment course.

4. Clinic data collection (Figure 1b)

1. Perform clinical assessments using MADRS²⁴, Hamilton Depression Rating Scale (HAMD)²⁷, Beck Depression Inventory-II (BDI-II)²⁸, Hamilton Anxiety Scale (HAMA)²⁹, Clinical Global Impression (CGI)³⁰ and MATRICS Consensus Cognitive Battery (MCCB)^31,32.
   NOTE: MINI and MADRS are used for the screening. All the above scales are applied for pre- and post-treatment clinical assessment.

Results

ROI-wise FC analysis should show that sgACC is significantly anti-correlated with DLPFC, in which the strongest negative correlation is the stimulus target to be chosen. Significant anti-correlation between the sgACC-DLPFC functional connectivity and the treatment response should be found in the correlation analysis³³.

The current protocol is based on an innovative TMS targeting method that no previous studies have applied. Here we present results from an fMRI-guided TM...

Discussion

The sgACC is responsible for emotional processing and plays an important role in stress regulation^16,17,18. A study suggests that targeting on a subregion of DLPFC that shows strongest functional anti-connectivity with sgACC (6, 16, -10) may enhance the antidepressant effect²⁵. Therefore, precisely locating this target is the critical step of this protocol. Before the stimulation, the borders of the brain...

Materials

List of materials used in this article

Name	Company	Catalog Number	Comments
3T Philips Achieva MRI scanner	Philips
Harvard/Oxford cortical template	http://www.cma.mgh.harva rd.edu/
MATLAB	MathWorks
SPM12	http://www.fil.ion.ucl.ac.uk/spm
The Visor2 system	ANT Neuro		The Visor2 software, the optical tracking system, tracking tools and calibration board are part of the visor2 system.
TMS device	Magstim, Carmarthenshire, UK