Abstract
Background: Deep brain stimulation targeting the subcallosal cingulate gyrus (SCG DBS) improves the symptoms of treatment-resistant depression in some patients, but not in others. We hypothesized that there are pre-existing structural brain differences between responders and nonresponders to SCG DBS, detectable using structural MRI.
Methods: We studied preoperative, T1-weighted MRI scans of 27 patients treated with SCG DBS from 2003 to 2011. Responders (n = 15) were patients with a > 50% improvement in Hamilton Rating Scale for Depression score following 12 months of SCG DBS. Preoperative subcallosal cingulate gyrus grey matter volume was obtained using manual segmentation by a trained observer blinded to patient identity. Volumes of hippocampus, thalamus, amygdala, whole-brain cortical grey matter and white matter volume were obtained using automated techniques.
Results: Preoperative subcallosal cingulate gyrus, thalamic and amygdalar volumes were significantly larger in patients who went on to respond to SCG-DBS. Hippocampal volume did not differ between groups. Cortical grey matter volume was significantly smaller in responders, and cortical grey matter:white matter ratio distinguished between responders and nonresponders with high sensitivity and specificity.
Limitations: Normalization by intracranial volume nullified some between-group differences in volumetric measures.
Conclusion: There are structural brain differences between patients with treatment-resistant depression who respond to SCG DBS and those who do not. Specifically, the structural integrity of the subcallosal cingulate gyrus target region and its connected subcortical areas, and variations in cortical volume across the entire brain, appear to be important determinants of response. Structural MRI shows promise as a biomarker in deep brain stimulation for depression, and may play a role in refining patient selection for future trials.
Introduction
Deep brain stimulation (DBS) is a well-established therapy for patients with movement disorders such as Parkinson disease, essential tremor and dystonia.1 Recent attempts have been made to apply DBS to patients with refractory psychiatric disorders, including obsessive–compulsive disorder2 and treatment-resistant depression (TRD).3,4 In particular, DBS for TRD has generated considerable enthusiasm because of an anticipated increase in depression-related disease burden,5 and the immense socioeconomic costs of depression worldwide.6 Multiple small open-label trials of DBS in patients with TRD have been published to date, using various stimulation targets including the ventral capsule/ventral striatum,7 nucleus accumbens,8 medial forebrain bundle,9 inferior thalamic peduncle,10 lateral habenula11 and subcallosal cingulate gyrus with its subjacent white matter (referred to as the SCG from here on).12–15 These open-label studies largely report improvement in depression severity scores with long-term stimulation irrespective of target, although at present there are insufficient data to conclusively demonstrate the superiority of any one target over another.
To date, the therapeutic effectiveness of DBS for TRD has not yet been conclusively established in randomized clinical trials using sham stimulation and blinded outcome assessment. In fact, a recent trial of ventral capsule/ventral striatum DBS failed to show a significant difference in TRD response rate between active stimulation and sham groups at the end of a 16-week controlled phase.16 Similarly, a randomized, double-blind, sham-controlled trial of SCG DBS failed to demonstrate superior antidepressant efficacy in the active stimulation group at 6 months, and, following an interim futility analysis, the trial was stopped prematurely by its industry sponsor.17 However, clinically meaningful antidepressant responses were observed at long-term follow-up in patients who continued to receive open-label, active stimulation beyond the initial 6-month randomized phase (48% response, 25% remission at 2 years). As a result, SCG DBS for TRD is still considered to have promise, although it is clear that we lack a complete understanding of the factors that may influence antidepressant response.17 Indeed, the inherent heterogeneity of TRD makes it unlikely that all TRD patients will respond to DBS, and broadly offering DBS to all patients who meet criteria for a diagnosis of TRD may result in a substantial rate of treatment failure. Because DBS is an invasive and costly therapy with the potential for serious adverse effects, it would be useful to have reliable, objectively measured biomarkers of expected antidepressant response that could be measured before surgical intervention.
A small number of studies have examined potential factors influencing response to SCG DBS, including DBS electrode targeting,18,19 stimulation parameters,20 variations in specific electroencephalographic features21,22 and pre-DBS neurocognitive function.23 We hypothesized that variations in preoperative neuroanatomy might similarly moderate differences in response to SCG DBS; this hypothesis follows from previous findings suggesting that pretreatment hippocampal volume may predict medication treatment response in patients presenting with major depressive disorder.24 Anatomically, the SCG grey matter is composed of cortex belonging to Brodmann area 25, as well as the caudal portions of areas 24 and 32,25 and is cytoarchitectonically distinct from the more dorsal anterior cingulate region.26 Projections to and from the SCG have largely been inferred from anatomic studies in non-human primates, and a small number of diffusion tensor imaging tractography studies in humans (see Hamani and colleagues25 for a comprehensive review of SCG projections). In brief, the SCG has prominent connections with prefrontal, orbitofrontal and anterior cingulate cortices. Subcortically, strong afferent and efferent connections have consistently been found between the SCG and the hippocampus and amygdala, consistent with the key roles of these structures in emotional processing as components of the limbic system. Additionally, there are connections between the SCG and multiple thalamic and hypothalamic nuclei.
In this study, we applied quantitative morphometric analysis techniques to preoperatively acquired structural MRI scans from 27 SCG DBS patients (15 responders, 12 nonresponders) treated for TRD. Our prediction was that structural brain differences between responders and nonresponders would be detectable on MRI scans, and would reflect the relative integrity of the SCG and closely connected brain regions. Our approach involved 3 different analyses: measurement of SCG grey matter volume as a surrogate measure of stimulation target integrity; volumetric measurements of 3 key subcortical structures—the hippocampus, amygdala and thalamus—which are connected to the SCG region and which have been implicated in depression; and measurement of whole brain cortical grey and white matter volume.
Methods
Patient selection, DBS parameters and follow-up
This retrospective cohort study was approved by the research ethics boards of the University Health Network and the Centre for Addiction and Mental Health, Toronto, Canada. We searched the Toronto Western Hospital DBS program database to identify all patients who had undergone implantation of a bilateral DBS system targeted to the SCG between January 1, 2003, and December 31, 2011. In total, we identified 43 patients, all with a diagnosis of TRD. Of these patients, 8 had missing or poor-quality preoperative standard volumetric stereotactic T1-weighted MRI scans (see MRI acquisition and preprocessing, below). We then also excluded patients who were treated as part of a multicentre randomized controlled trial (n = 8), because trial patients routinely underwent a period of sham stimulation. This left a total of 27 patients for analysis. Waiver of informed consent for the present study was granted by the institutional research ethics board, because it involved only the retrospective use of de-identified imaging data; however, all patients had originally provided written informed consent to proceed with bilateral SCG DBS after having been made aware of the potential risks and benefits of the procedure.
All analyzed patients had been selected for SCG DBS based on standard inclusion and exclusion criteria at our centre.27 Of note, all patients had undergone at least 1 course of psychotherapy before SCG DBS implantation, but not all patients had undergone electroconvulsive therapy. Technical aspects of the bilateral SCG DBS implantation procedure have been detailed in previous publications.14,27 Similarly, guidelines for stimulation parameter selection have also been published previously.27 Responders were defined as patients who, after 12 months of SCG DBS, had achieved a 50% or greater reduction in Hamilton Rating Scale for Depression-17 (HRSD-17) scores compared with baseline. For details related to patient selection, DBS programming and psychiatric follow-up, see Appendix 1, Supplemental Materials, available at jpn.ca/180207-a1.
MRI acquisition and preprocessing
The MRI scans analyzed in this study were standard preoperative stereotactic scans obtained on the morning of DBS implantation surgery. All MRI scans were acquired on the same 1.5 T GE Signa EXCITE scanner (GE Healthcare). A T1-weighted, 3-dimensional spoiled gradient recalled volumetric sequence was obtained with the following parameters: repetition time 11.9 ms, echo time 5.0 ms, inversion time 300 ms, flip angle 20°, field of view 25.9 cm × 25.9 cm, matrix 256 × 256, reformatted into 1.4 mm axial slices (with 0.7 mm overlap to prevent staircase artifact). Effective voxel size was 0.5 × 0.5 × 1.4 mm = 0.35 mm3. All scans underwent N3 intensity nonuniformity correction28 and realignment to the anterior–posterior commissure plane using Medical Image Processing and Visualization (MIPAV) version 5.2.1 (National Institutes of Health). Images were resampled into isotropic 1 mm3 voxels using trilinear interpolation, as per our standard laboratory protocol.29
Volumetric analysis of SCG grey matter
To assess whether the anatomic attributes of the stimulation target affect response to SCG DBS, we measured SCG grey matter volume on preoperative structural MRI scans. A trained observer blinded to patient identity and clinical outcome measured left and right SCG volumes by manual segmentation. All manual segmentations were performed using ITK-SNAP version 2.4.0 (www.itksnap.org). The SCG grey matter was segmented according to the protocol of Drevets and colleagues (Fig. 1).30 Segmentations were performed on coronal sections. The anterior boundary of the SCG was the anterior-most coronal slice containing the corpus callosum, and the posterior boundary was the anterior-most slice where the internal capsule no longer divided the striatum. All grey matter belonging to the first full gyrus below the corpus callosum was included in the SCG volume. Occasionally, the SCG appeared as 2 contiguous small double gyri instead of a single larger gyrus; double gyri were always included together in the total SCG volume. Segmented volumes were inspected on sagittal slices to ensure a consistent inferior boundary along the antero-posterior extent of the SCG (Fig. 1B).
To assess inter-rater reliability, a second blinded observer independently segmented the SCG in 14 randomly selected patients. Reliability was quantified using the intraclass correlation coefficient. We report an absolute agreement standard for intraclass correlation coefficient, which for each structure represents the variance of the measurement, divided by the sum of the variance of the measurement and the variance over patients.31
Volumetric analysis of hippocampus, amygdala and thalamus
To test the hypothesis that the integrity of brain structures connected to the SCG may predict eventual response to SCG DBS, we measured the volumes of the hippocampus, amygdala and thalamus, which have strong afferent and efferent connections to the SCG.25 Volumetry of these structures was performed by automated segmentation using FSL-FIRST, part of the FSL toolkit32 (http://fsl.fmrib.ox.ac.uk/fsl/fslwiki/). FIRST is a model-based segmentation tool that uses shape- and appearance-based models constructed from manually segmented images.33 The manual labels are parameterized as surface meshes and modelled as a point distribution model. Deformable surfaces are used to automatically parameterize the volumetric labels using constraints to preserve vertex correspondence across the training data. Based on learned models, FIRST searches through linear combinations of shape modes of variation for the most probable shape instance given the observed intensities in a T1-weighted image. FIRST was implemented using FSL version 4.1.9. FIRST segmentations were performed using the run_first_all script according to the FIRST user guide (http://fsl.fmrib.ox.ac.uk/fsl/fslwiki/FIRST/UserGuide).
All segmentations generated by FIRST were visually inspected for gross errors by a single observer blinded to patient identity and clinical outcome.
Whole brain grey and white matter analysis
To assess whether global brain structure differences could influence response to SCG DBS in depression, we compared whole brain grey and white matter volumes between responders and nonresponders. Brain tissue volume, normalized for patient head size, was estimated with SIENAX,34 part of FSL.32 SIENAX starts by extracting brain and skull images from the single whole-head input data.35 The brain image is then affine-registered to MNI152 space36,37 (using the skull image to determine the registration scaling); this is primarily to obtain the volumetric scaling factor, to be used as a normalization for head size. Next, tissue-type segmentation with partial volume estimation is carried out38 to calculate the total volume of brain tissue (including separate estimates of volumes of grey matter, white matter, peripheral grey matter and cerebrospinal fluid). We ran SIENAX with the “-r” option to calculate cortical grey matter volumes (i.e., excluding subcortical and cerebellar grey matter volumes). We used cortical grey matter:white matter ratio as a composite global structural metric in comparisons between responders and nonresponders.
Normalization by intracranial volume
A common approach in quantitative volumetric imaging studies of the brain is to use intracranial volume (ICV) as a normalization factor when comparing the volumes of various brain structures and regions between groups.39 Theoretically, this ensures that head size does not confound between-group comparisons of structural volumes. A surrogate measure of ICV was determined for each patient using the scaling factor (VSCALING) generated by SIENAX, which represents the amount by which each patient’s skull volume must be multiplied to be transformed into a standard image space (MNI 152). A smaller scaling factor implies a larger ICV. All raw structural volumes were multiplied by the scaling factor and then subjected to further statistical analysis.
Statistical analysis
Given our small sample size, we compared the volumes of the SCG, hippocampus, amygdala and thalamus SCG DBS responders and nonresponders using the nonparametric Mann–Whitney test. Similarly, we compared cortical grey matter volume, whole brain white matter volume, and cortical grey matter:white matter ratio between responders and nonresponders using the Mann–Whitney test. All analyses were performed for both raw volumes and volumes normalized by SIENAX ICV scaling factor. To assess the discriminatory power of cortical grey matter:white matter ratio for response to SCG DBS, we constructed receiver operating characteristic curves, calculating the areas under the curves and the corresponding 95% confidence intervals using the method of DeLong and colleagues.40 We performed correlation analysis between structural variables of interest and HRSD-17 scores using Spearman correlation. In all statistical tests, we used p < 0.05 as the threshold for statistical significance. We performed Bonferroni correction to account for multiple comparisons in the volumetric analysis of the hippocampus, amygdala and thalamus (i.e., cut-off for significant p value 0.05/3 = 0.017). Given the exploratory nature of the study, the small sample size and the use of nonparametric statistics, we did not control for confounding variables in any analyses.
Results
Clinical and demographic variables
A summary comparison of clinical and demographic variables between responder and nonresponder groups is shown in Table 1. We found no differences for age at SCG DBS implantation; sex; HRSD-17 score; number of previous major depressive episodes; years since onset of major depressive disorder; duration of current depressive episode; number of antidepressant medication trials before DBS; or the proportion of patients previously treated with electroconvulsive therapy, treated with DBS during a first major depressive episode, previously treated with antipsychotics, or with a family history of major depressive disorder. Predictably, responders demonstrated a significantly larger median percentage improvement in HRSD-17 score than nonresponders after 1 year of stimulation (U = 7; p < 0.0001).
Subcallosal cingulate gyrus volume
Inter-rater reliability for SCG volume was excellent, with an intraclass correlation coefficient of 0.91. Volumetric results for SCG are summarized in Fig. 2. The median volume of the left SCG was significantly larger in responders to SCG DBS than in nonresponders (Fig. 2A; U = 44.5; p = 0.028). Median volume of the average SCG was also significantly larger in responders (Fig. 2C; U = 42.0; p = 0.021). Overall, there was no difference in left versus right median SCG volume in either the responder or the nonresponder group.
Average SCG volume was not significantly correlated with percent change in HRSD-17 scores at 12 months.
Hippocampus, amygdala and thalamus volume
Fig. 3 summarizes the results of automated volumetry for the hippocampus, amygdala and thalamus across responders and nonresponders. Visual inspection by an observer blinded to patient identity and clinical outcome identified grossly erroneous hippocampal segmentations in 3 of 27 patients (2 responders, 1 nonresponder), amygdalar segmentations in 3 of 27 patients (2 responders, 1 nonresponder), and thalamic segmentations in 2 of 27 patients (1 responder, 1 nonresponder). These patients were excluded from subsequent between-group statistical comparisons.
Left, right or average hippocampal volumes were not significantly different between responders and nonresponders. In contrast, average amygdala volume was significantly larger in responders compared with nonresponders (U = 28.00; p = 0.013) after correcting for multiple comparisons; on their own, left and right amygdala volumes showed strong trends toward being larger in responders, but these findings did not survive Bonferroni correction (p = 0.05 and p = 0.037, respectively). Left, right, and average thalamic volume were all significantly larger in responders after correction for multiple comparisons (U = 30.00, p = 0.011 for all).
We found a significant correlation between pre-DBS average thalamic volume and percent change in 12-month HRSD-17 scores (Spearman’s r = 0.54, p = 0.0057). We found no significant correlations between hippocampal or amygdala volume and percent change in HRSD-17 scores.
Whole brain grey and white matter analysis
Total cortical grey matter volume was significantly lower in responders than in nonresponders (U = 37.00; p = 0.010), but there was no difference in total white matter volume between both groups (Fig. 4A, B). Cortical grey matter:white matter ratio was also significantly lower in responders (U = 27.00; p = 0.002, Fig. 4C). Overall, the discriminatory power of cortical grey matter:white matter ratio for response to SCG DBS was high: the area under the receiver operating characteristic curve for cortical grey matter:white matter ratio was 0.85 (Fig. 4D). Cortical grey matter volume and cortical grey matter:white matter ratio were both well correlated with percent improvement in HRSD-17 scores after 12 months of DBS (Spearman’s r = 0.39, p = 0.043 and Spearman’s r = 0.41, p = 0.036, respectively).
Effect of normalization by intracranial volume
Normalization by ICV nullified any significant differences in SCG, hippocampus, amygdala and thalamus volume between responders and nonresponders. To further explore this finding, we compared median SIENAX VSCALING factor between responders and nonresponders and found that it was significantly smaller in responders (U = 23.00; p = 0.0019; Appendix 1, Figure S1). Put differently, responders had larger ICV than nonresponders on average. When we employed a different method for measuring ICV (Appendix 1, Supplemental Methods), we again found that ICV was significantly larger in responders than in nonresponders (U = 35.00; p = 0.0078; Appendix 1, Figure S2), with excellent correlation between ICV measurements by both methods (Appendix 1, Figure S3).
Discussion
In patients with TRD, we identified differences in preoperative brain structure between responders and nonresponders to SCG DBS therapy. Specifically, we found evidence of larger SCG, amygdalar and thalamic volume in responders, but no difference in hippocampal volume between groups. Preoperative thalamic volume was particularly well correlated with improvement in depression scores. At the brain-wide level, we observed a greater volume of cortical grey matter in nonresponders, and the cortical grey matter:white matter ratio performed well in discriminating between responders and nonresponders. Taken together, these novel findings suggest that neuroanatomical variability may influence eventual response to SCG DBS, as may clinical factors or DBS electrode position.
Cortical and subcortical volumetric predictors of response to SCG DBS
Several volumetric MRI studies have examined the SCG region — and in particular its grey matter volume — in unipolar and bipolar depression.30,41–50 Although the definition of the anatomic boundaries of the SCG varies across studies and effect sizes are variable, existing data suggest that SCG volume is decreased in unipolar depression, with more severe atrophy in the left SCG than in the right. Such SCG volume loss is observed in depressed patients as early as the initial episode of major depression and is more severe in patients with a family history of depression, suggesting that SCG atrophy may reflect susceptibility to — rather than a secondary consequence of — ongoing depression.50 Accordingly, our finding of reduced SCG volume in DBS nonresponders, particularly on the left, may suggest that nonresponders are neuroanatomically susceptible to TRD. Another possibility is that smaller SCG volume in nonresponders may reflect poorer integrity of the target region for DBS, leading to poorer clinical response. However, it has been argued that the putative SCG DBS target may actually be the SCG white matter, with preliminary evidence that optimal antidepressant effects require stimulation of the confluence of 3 white matter bundles connecting the SCG with the medial frontal cortex (via forceps minor and uncinate fasciculus), the rostral and dorsal cingulate cortex (via the cingulum bundle) and subcortical nuclei, respectively.19 Unfortunately, we did not acquire diffusion tensor imaging scans in our preoperative imaging protocol, so we were unable to assess the preoperative integrity of these key white matter bundles. The extent to which grey and white matter alterations in the SCG region are correlated, and how this relationship might influence DBS response, remains to be determined.
An association between depression and hippocampal atrophy has consistently been reported across several studies, including 5 separate meta-analyses of MRI volumetric studies in depression.51–55 The exact nature of the association remains uncertain. There is some evidence that a smaller hippocampus may predispose to depression,56 further supported by the association between smaller hippocampal volume and both a family history of depression and early-life adversity.57 Hippocampal atrophy is also correlated with longer overall duration of depressive illness,52,55 while hippocampal volume appears more significantly reduced in the depressed rather than the remitted state.53,56,58 Finally, smaller hippocampal volume appears to predict poorer response to antidepressant treatment.59 We found no differences in pre-DBS hippocampal volume between responders and nonresponders to SCG DBS in our population of TRD patients. It is possible that the similar proportion of patients with a family history of depression, similar duration of illness and similar duration of antidepressant therapy (see Table 1) may have attenuated any subtle inter-group pre-DBS variability in hippocampal volume. Furthermore, the hippocampus is a highly neuroplastic structure whose volume can be influenced by a multitude of factors, including aerobic exercise, glucocorticoid levels (elevated in depression), electroconvulsive therapy and mood-altering medications (notably lithium).60 These factors may have affected hippocampal volume to varying degrees across our patient population, possibly confounding any pre-DBS differences between responders and nonresponders.
We found significantly smaller thalamic and amygdalar volumes in nonresponders to SCG DBS. A recent meta-analysis of volumetric neuroimaging studies reported that the thalamus volume is consistently reduced in major depressive disorder.53 Thalamic atrophy is present even in early, treatment-naïve depression, suggesting that it is potentially a trait marker of susceptibility to depression.61 Further supporting this are recent data demonstrating that thalamic volume loss may be a key diagnostic feature (at the individual patient level) of pediatric unipolar depression, regardless of disease duration or extent of previous psychotropic medication use.62 There is also growing evidence that thalamic volume may be modulated by genetic factors known to be important in depression, including the Val66Met polymorphism in the brain-derived neurotrophic factor gene, as well as the serotonin transporter gene-linked polymorphism.63 Analogous to our SCG volume findings, thalamic atrophy in nonresponders may indicate a more neuroanatomically severe form of depression resistant to beneficial effects of SCG DBS. Amygdalar atrophy in nonresponders may reflect the same phenomenon, although the existing data on amygdalar volume in depression are variable and inconclusive,53,56,64–67 despite clear evidence for the key role of the amygdala in mood disorders.
Whole brain structural predictors of response to SCG DBS
Arguably, the most intriguing finding from our structural neuroimaging analysis was significantly reduced cortical grey matter volume and grey matter:white matter ratio in responders, with cortical grey matter:white matter ratio in particular revealing itself as a robust means of segregating responders from nonresponders. At first glance, these results appear counterintuitive, since many studies demonstrate significant thinning across various frontal cortical regions in patients with chronic depression, including the orbitofrontal cortex,68–70 dorsal anterior cingulate cortex,71 rostral middle frontal gyrus58,69 and dorsolateral prefrontal cortex.72 It also appears that some cortical regions — primarily in the frontal lobe but also encompassing the inferior temporal gyrus — may increase in thickness with successful medical therapy that leads to remission.58 However, recent work has shown increased thickness in the orbitofrontal cortex in people at increased familial risk of depression,73 and other studies have demonstrated increased orbitofrontal74 and cingulate/temporopolar cortical thickness in patients with first-episode, medication-naïve depression compared with controls.75 These findings support the notion that greater cortical thickness in certain brain regions may be a trait marker in depression, and therefore in line with the possibility, as suggested above, that nonresponders exhibit a neuroanatomically more severe form of depression. This neuroanatomical profile may also reflect a relative failure of cortical pruning during adolescent brain development in nonresponders, which has been linked to the development of depression.76 To our knowledge, there have been no previously published attempts to use the SIENAX tool in characterizing outcome or response to treatment. Most cortical thickness studies in depression have employed a vertex-based approach in which regionally specific clusters of difference in cortical thickness are reported between groups. However, this approach is insensitive to smaller variations in cortical thickness when averaged over larger regions or the entire brain.58 Our approach generates a simple scalar measure (i.e., cortical grey matter:white matter ratio) which can serve as a rapidly measured, automatically generated biomarker at the individual patient level.
The potential influence of ICV on our results is a key issue that merits discussion. Briefly, we found that median ICV — measured using 2 different methods — was significantly different between responders and nonresponders to SCG DBS. As a result, when normalizing by ICV — to correct for the potential confounding effect of head size and premorbid brain size differences between responders and nonresponders — significant between-group differences in volume of the SCG, amygdala and thalamus were abolished. Following a thorough literature search, we could not find any direct evidence of differences in average ICV between depressed patients and controls, or between responders and nonresponders to antidepressant therapy. However, there are some links between head size and inherited genetic polymorphisms (e.g., locus 17q2177,78) implicated in the development of depression. In our particular cohort, the ICV may have differed between response groups because of an uneven distribution of males and females in each group (i.e., the male:female ratio was higher in responders), although, strictly speaking, this did not meet the threshold for a statistically significant inter-group difference.
Accounting for ICV in volumetric MRI studies is a controversial topic. An important concern is that the specific method used to compute ICV can significantly influence between-group volumetric comparisons, with different ICV methods leading to divergent conclusions.79 Another concern is that the relationship between ICV and regional brain volumes is not necessarily linear, so that simple division of a structure’s volume by ICV — as performed in most studies — may be inappropriate.80 This limitation could be addressed by using ICV as a covariate in a general linear model; however, given our small sample size we chose to use nonparametric statistical methods that did not easily lend themselves to adjustment for covariates.81 For these reasons, we believe that our raw volumetric results should not necessarily be dismissed out of hand because of the effect of ICV normalization, although they will require validation in a larger study. It is also worth mentioning that cortical grey matter and associated cortical grey matter:white matter ratio results are already head size–corrected by the SIENAX algorithm, and are therefore unaffected by between-group differences in ICV.
Limitations
In addition to the influence of ICV on our volumetric results, there are several other reasons to interpret our results with caution. Most obviously, the sample size is small, which limits the robustness of any reported differences between responder and nonresponder groups. Both the small sample size and the obligatory use of nonparametric statistical methods precluded any meaningful attempt to control for confounding variables during inter-group comparisons. Second, MRI data were collected over a long period (8 yr), during which there were multiple software upgrades to the MRI scanner that could result in minor changes to image contrast, in turn influencing volumetric analysis (although the same scanner was used throughout the study). As well, all MRI data were resampled to generate isotropic voxels, which may have resulted in interpolation errors affecting image quality.82
Conclusion
Our findings suggest that TRD patients who do not respond to SCG DBS may exhibit a more neuroanatomically severe form of depression that may not necessarily be reflected in demographic or clinical features. Larger studies are required to determine if structural MRI can be used to improve patient selection in future trials of DBS in TRD.
Footnotes
Funding: T. Sankar’s work on this manuscript was supported by a Canadian Institutes of Health Research (CIHR) fellowship award. A. Lozano was supported by a Canada Research Chair in Neuroscience and the R.R. Tasker Chair in Functional Neurosurgery.
Competing interests: M. Chakravarty is a member of the JPN editorial board; he was not involved in the decision-making on this manuscript. P. Giacobbe and S. Kennedy have received honoraria from St. Jude Medical, Inc. H. Mayberg has received consulting and intellectual property fees from St. Jude Medical, Inc. C. Hamani is a consultant for St. Jude Medical, Inc. A. Lozano is a consultant to Medtronic, Inc., St. Jude Medical, Inc., and Boston Scientific, Inc.; serves on the scientific advisory board of Ceregene, Codman, Neurophage, Aleva and Alcyone Life Sciences; is co-founder of Functional Neuromodulation Inc.; and holds intellectual property in the field of deep brain stimulation. T. Sankar, N. Jawa, S. Li and S. Rizvi report no biomedical financial interests or other conflicts of interest.
Contributors: T. Sankar, C. Hamani and A. Lozano designed the study. T. Sankar, N. Jawa, S. Li, P. Giacobbe, S. Kennedy and S. Rizvi acquired the data, which T. Sankar, M. Chakravarty, P. Giacobbe, S. Kennedy and H. Mayberg analyzed. T. Sankar, M. Chakravarty, S. Kennedy and S. Rizvi wrote the article, which all authors reviewed. All authors approved the final version to be published and can certify that no other individuals not listed as authors have made substantial contributions to the paper.
- Received October 31, 2018.
- Revision received March 20, 2019.
- Accepted April 8, 2019.