Regional electric field induced by electroconvulsive therapy in a realistic finite element head model: Influence of white matter anisotropic conductivity
Highlights
► Anatomically-realistic head models allow ROI analysis of ECT electric field. ► Bilateral ECT produces the strongest electric field overall and in hippocampus. ► Neglecting white matter anisotropy results in electric field errors up to 39%. ► Computational models could help optimize ECT efficacy and tolerability.
Introduction
Electroconvulsive therapy (ECT) is a therapeutic intervention in which electric current is applied through scalp electrodes to induce a generalized seizure in anesthetized patients (Abrams, 2002, APA, 2001). Although ECT plays a vital role in the treatment of medication-resistant psychiatric disorders, such as major depression, the use of ECT has been limited by its cognitive side effects (particularly amnesia (Sackeim et al., 2007, Squire, 1986)), by cardiac complications (Nuttall et al., 2004), by the need for general anesthesia, as well as by the high rate of relapse (Sackeim et al., 1990). Despite the introduction of various improvements of ECT technique, there is still limited knowledge of how to optimally select electrode placement (Kellner et al., 2010b) or stimulus current parameters (Peterchev et al., 2010) for maximal efficacy and tolerability. Indeed, the therapeutic action and adverse side effects of ECT are highly dependent upon electrode placement and stimulus current parameters, but a complete mechanistic explanation for these relationships is still lacking. For instance, right unilateral (RUL) ECT leads to fewer cognitive side effects than bilateral frontotemporal (BL) ECT (Sackeim et al., 2000), but it is not known whether this is by virtue of lower electric field (E-field) strength in hippocampus and other regions crucial for memory. Furthermore, alternative ECT electrode configurations such as bifrontal (BF) (Abrams, 2002) and focal electrically administered seizure therapy (FEAST) (Sackeim, 2004) have been proposed with the goal of preferentially targeting frontal brain regions to reduce memory impairment, but the frontal E-field strength relative to the rest of the brain and relative to other electrode placements has not been quantified.
To understand the underlying biophysical mechanisms of ECT, a few early studies undertook measurements of the E-field generated by ECT in human cadavers (Lorimer et al., 1949, Smitt and Wegener, 1944) and in an electrolytic tank containing a human half-skull (Rush and Driscoll, 1968). However, the electrolytic tank measurements did not account for the geometry and conductivity properties of the scalp and the brain. The intracerebral cadaver measurements were carried out after an uncontrolled interval of time following death, potentially resulting in altered conductivity profile of the head tissues, and the tissues were damaged in the process of inserting the recording probes, potentially altering the paths of current flow generated by the scalp electrodes. Furthermore, neither of the studies produced a high-resolution map of the E-field or the current density distributions in the brain.
In order to provide more detailed field maps, a number of computational studies have simulated the distribution of the E-field or the current density field (which equals the product of E-field and conductivity) in the brain using a volume conductor model of the head. The representation of the head in computational ECT models ranges in detail from concentric spheres (Deng et al., 2009, Deng et al., 2011, Weaver et al., 1976) to low-resolution realistically-shaped representations (Bai et al., 2011, Sekino and Ueno, 2002, Sekino and Ueno, 2004) to high-resolution anatomically-accurate models (Nadeem et al., 2003, Szmurło et al., 2006). Furthermore, a substantial number of E-field/current density modeling studies have been published in the context of other transcranial electric stimulation paradigms, again ranging from simplified to realistic head representations (Datta et al., 2008, Datta et al., 2009, Grandori and Rossini, 1988, Holdefer et al., 2006, Im et al., 2008, Lee et al., 2009b, Miranda et al., 2006, Miranda et al., 2007, Nathan et al., 1993, Oostendorp et al., 2008, Parazzini et al., 2011, Rush and Driscoll, 1968, Sadleir et al., 2010, Salvador et al., 2010, Saypol et al., 1991, Stecker, 2005, Suh et al., 2009, Suh et al., 2010, Suihko, 2002, Wagner et al., 2007). However, these studies have various limitations. The spherical and simplified geometry models do not fully account for tissue inhomogeneity and anisotropy, and the complex geometries of head tissues, including orifices in the skull such as the auditory canals and the orbits. The published anatomically-accurate ECT models (Nadeem et al., 2003, Szmurło et al., 2006) consider only isotropic tissue conductivity, explore only a limited set of electrode configurations (BL and RUL), and do not perform region of interest (ROI) analysis of the field distribution in the brain. The computational models of non-ECT transcranial electric stimulation offer some insights into the biophysics of the problem, but do not provide data specific to ECT electrode configurations and stimulus current parameters.
For realistic models of the E-field generated by ECT, the inclusion of anisotropic conductivity of the white matter (WM) may be of particular importance since the E-field induced by ECT is typically widespread and reaches deep brain regions (Deng et al., 2011, Nadeem et al., 2003), and since depression itself is associated with regionally specific abnormalities of the WM fractional anisotropy (Korgaonkar et al., 2011, Wu et al., 2011). Our and other groups have previously incorporated tissue anisotropic conductivity in models of electroencephalography and magnetoencephalography (Gullmar et al., 2006, Gullmar et al., 2010, Hallez et al., 2008, Hallez et al., 2009, Haueisen et al., 2002, Kim et al., 2003, Lee et al., 2008, Lee et al., 2009a, Marin et al., 1998, Rullmann et al., 2009, Wolters et al., 2006), transcranial direct current stimulation (tDCS) (Lee et al., 2009b, Oostendorp et al., 2008, Suh et al., 2009, Suh et al., 2010), deep brain stimulation (Butson et al., 2007), transcranial magnetic stimulation (TMS) (De Lucia et al., 2007, Thielscher et al., 2011), and electrical impedance tomography (Abascal et al., 2008). These studies demonstrate that anisotropic conductivity of the brain tissue can have a non-negligible effect on the electromagnetic field solutions. However, computational models of ECT have not incorporated tissue conductivity anisotropy to date.
No direct and non-invasive in vivo measurement of brain conductivity anisotropy is available, but the similarity between the transportation processes of electrical charge carriers and water molecules enables estimation of the effective electrical conductivity tensors from the water self-diffusion tensors which can be non-invasively acquired with diffusion tensor magnetic resonance imaging (DT-MRI) (Basser et al., 1994b). Several methods have been proposed to derive the WM anisotropic conductivity from the measured diffusion tensors. In the effective medium approach (Tuch et al., 1999, Tuch et al., 2001), the WM anisotropic conductivity tensors were directly calculated by a linear scaling of the diffusion tensors using an empirically determined scaling factor (Haueisen et al., 2002, Tuch et al., 1999, Tuch et al., 2001). However, Rullmann et al., 2009, Gullmar et al., 2010 have pointed out that using this linear scaling approach may lead to extremely large anisotropic ratios in the resulting conductivity tensors. An alternative is the volume constraint approach where the WM anisotropic conductivity tensors are computed with a fixed anisotropic ratio in each WM voxel, under the assumption that the shape of the WM diffusion tensors is prolate (cigar-shaped), rotationally symmetric ellipsoid (Shimony et al., 1999, Wolters et al., 2006). With this method, the fixed anisotropic conductivity ratio of the WM tissue can be obtained from direct measurements, e.g., 10:1 for parallel:normal orientation relative to the nerve fibers (Nicholson, 1965). Another anisotropy modeling technique is based on the linear conductivity-to-diffusivity relationship in combination with a constraint on the magnitude of the electrical conductivity tensor (Hallez et al., 2008, Hallez et al., 2009). A “volume fraction algorithm” considering the partial volume effects of the cerebrospinal fluid (CSF) and the intravoxel fiber crossing structure has also been suggested (Wang et al., 2008), but no further studies using this approach have been reported.
In summary, existing studies of the E-field or current density resulting from ECT have investigated few electrode configurations in realistic-geometry head models, have not incorporated tissue conductivity anisotropy, and have not carried out analysis of the E-field strength in specific brain ROIs. Addressing these limitations, in the present study we develop an anatomically-accurate finite element (FE) model of the human head incorporating tissue heterogeneity and WM anisotropic conductivity, based on individual structural MRI and DT-MRI scans. We use the head model to simulate the E-field generated in the brain by the BL, BF, RUL, and FEAST ECT electrode configurations. We quantify the differences in E-field strength among the various ECT electrode configurations in brain ROIs that have putative role in the therapeutic action and/or adverse side effects of ECT. This analysis enables us, for example, to explore whether forms of ECT associated with fewer cognitive side effects induce lower E-field strengths in hippocampus, and to evaluate the degree to which frontal electrode configurations (BF and FEAST) achieve focal frontal stimulation. We also investigate how the WM conductivity anisotropy affects the E-field distribution in the brain. This study demonstrates the utility of anatomically-realistic computational models to provide clinically salient analysis and recommendations for the optimization of ECT. Preliminary results from this study were previously presented in part in conference proceedings (Lee et al., 2010, Lee et al., 2011).
Section snippets
Materials and methods
The steps of the E-field modeling and analysis are diagrammed in Fig. 1 and described below.
3D finite element head model
The human head model used for the E-field simulation is displayed in Fig. 2. Fig. 2(a) shows a cut-away 3D rendering of the head model. The cropped section illustrates the five segmented tissue types using the color convention defined in Fig. 1. Fig. 2(b) shows a transaxial conductivity map with the principal orientation of the electrical conductivity tensors (corresponding to the orientation of the WM fibers) projected as black bars onto the WM regions. For clarity, a portion of the
Discussion and conclusions
Even though the E-field spatial distribution is a key aspect of dosage in ECT, it has not been accurately characterized. The spatial distribution of the E-field strength is a determinant of which brain regions are directly activated by the electric stimulation delivered by various ECT electrode configurations. This work represents the first quantitative study investigating the regional differences in E-field strength resulting from variations in the ECT electrode configuration in an
Acknowledgment
This work was supported by the National Institutes of Health under grant R01MH091083, the National Science Foundation through TeraGrid resources provided by National Center for Supercomputing Applications under grant TG-MCB100050, and the Korea Science and Engineering Foundation (KOSEF) grant funded by the Korean government (MEST) (No. 2009-0075462). We would like to thank Drs. Richard Weiner and Andrew Krystal from the Department of Psychiatry and Behavioral Sciences at Duke University for
References (128)
- et al.
Use of anisotropic modelling in electrical impedance tomography: description of method and preliminary assessment of utility in imaging brain function in the adult human head
Neuroimage
(2008) - et al.
Estimation of the effective self-diffusion tensor from the NMR spin-echo
J. Magn. Reson.
(1994) - et al.
MR diffusion tensor spectroscopy and imaging
Biophys. J.
(1994) - et al.
Selective frontal, parietal, and temporal networks in generalized seizures
Neuroimage
(2003) - et al.
Provably good sampling and meshing of surfaces
Graph. Model.
(2005) - et al.
Patient-specific analysis of the volume of tissue activated during deep brain stimulation
Neuroimage
(2007) - et al.
A structurally detailed finite element human head model for simulation of transcranial magnetic stimulation
J. Neurosci. Methods
(2009) - et al.
Diffusion tensor MRI-based estimation of the influence of brain tissue anisotropy on the effects of transcranial magnetic stimulation
Neuroimage
(2007) - et al.
A full subtraction approach for finite element method based source analysis using constrained Delaunay tetrahedralisation
Neuroimage
(2009) - et al.
Hypothalamic abnormalities in schizophrenia: sex effects and genetic vulnerability
Biol. Psychiatry
(2007)
Thalamic relay functions and their role in corticocortical communication: generalizations from the visual system
Neuron
Influence of anisotropic electrical conductivity in white matter tissue on the EEG/MEG forward and inverse solution. a high-resolution whole head simulation study
Neuroimage
A tractography analysis of two deep brain stimulation white matter targets for depression
Biol. Psychiatry
The influence of brain tissue anisotropy on human EEG and MEG
Neuroimage
Predicted current densities in the brain during transcranial electrical stimulation
Clin. Neurophysiol.
The effects of ECT stimulus dose and electrode placement on the ictal electroencephalogram: an intraindividual crossover study
Biol. Psychiatry
Influence of white matter anisotropic conductivity on EEG source localization: comparison to fMRI in human primary visual cortex
Clin. Neurophysiol.
Path of current distribution in brain during electro-convulsive therapy; preliminary report
Electroencephalogr. Clin. Neurophysiol.
MRI-based topographic parcellation of human cerebral white matter and nuclei II. Rationale and applications with systematics of cerebral connectivity
Neuroimage
Deep brain stimulation for treatment-resistant depression
Neuron
Anterior cingulate cortex: an MRI-based parcellation method
Neuroimage
Focal network involvement in generalized seizures: new insights from electroconvulsive therapy
Epilepsy Behav.
Modeling the current distribution during transcranial direct current stimulation
Clin. Neurophysiol.
Determination of current-density distributions generated by electrical-stimulation of the human cerebral-cortex
Electroencephalogr. Clin. Neurophysiol.
Specific impedance of cerebral white matter
Exp. Neurol.
Manual MRI parcellation of the frontal lobe
Psychiatry Res.
EEG source analysis of epileptiform activity using a 1 mm anisotropic hexahedra finite element head model
Neuroimage
Convulsant and anticonvulsant properties of electroconvulsive therapy: towards a focal form of brain stimulation
Clin. Neurosci. Res.
Stimulus intensity, seizure threshold, and seizure duration: impact on the efficacy and safety of electroconvulsive therapy
Psychiatr. Clin. North Am.
Effects of pulse width and electrode placement on the efficacy and cognitive effects of electroconvulsive therapy
Brain Stimul.
Transcranial direct current stimulation (tDCS) in a realistic head model
Neuroimage
Electroconvulsive Therapy
Diencephalic stimulation and the effects of ECT in endogenous depression
Br. J. Psychiatry
Conductivities of three-layer live human skull
Brain Topogr.
The Practice of Electroconvulsive Therapy: Recommendations for Treatment, and Privileging: A Task Force Report of the American Psychiatric Association
Effect of conductivity uncertainties and modeling errors on EEG source localization using a 2D model
IEEE Trans. Biomed. Eng.
A computational model of direct brain excitation induced by electroconvulsive therapy: comparison among three conventional electrode placements
Brain Stimul.
Experimental validation of the influence of white matter anisotropy on the intracranial EEG forward solution
J. Comput. Neurosci.
Guidelines for precise and accurate computational models of tDCS
Brain Stimul.
Analysis of the quasi-static approximation for calculating potentials generated by neural stimulation
J. Neural Eng.
Changes in brain volume in major depression
Depress. Mind Brain
Transcranial current stimulation focality using disc and ring electrode configurations: FEM analysis
J. Neural Eng.
Gyri-precise head model of transcranial direct current stimulation: improved spatial focality using a ring electrode versus conventional rectangular pad
Brain Stimul.
Effect of anatomical variability on neural stimulation strength and focality in electroconvulsive therapy (ECT) and magnetic seizure therapy (MST)
Conf. Proc. IEEE Eng. Med. Biol. Soc.
Improving the focality of electroconvulsive therapy: the roles of current amplitude, and electrode size and spacing
J. ECT
Electric field strength and focality in electroconvulsive therapy and magnetic seizure therapy: a finite element simulation study
J. Neural Eng.
A multiple electrode scheme for optimal non-invasive electrical stimulation
Segmentation of skull and scalp in 3-D human MRI using mathematical morphology
Hum. Brain Mapp.
The effects of repetitive transcranial magnetic stimulation in the treatment of depression
Expert Rev. Med. Devices
Development of volume conductor and source models to localize epileptic foci
J. Clin. Neurophysiol.
Cited by (94)
Electroconvulsive therapy effects on anhedonia and reward circuitry anatomy: A dimensional structural neuroimaging approach
2022, Journal of Affective DisordersCitation Excerpt :While our clinical understanding of ECT is significant, it contrasts with how little we know about its neurobiological mechanisms of action. Given how diffuse the induced electric fields are (Lee et al., 2012) and the fact that ECT leads to a generalized seizure (i.e., the entire brain seizes), ECT has been traditionally considered a diffuse and non-specific neuromodulation treatment (Abrams, 1991; Ottosson, 1960). However, previous neuroimaging ECT research has shown ECT is more focal than previously considered (Cano et al., 2017, 2019).
Electric field strength induced by electroconvulsive therapy is associated with clinical outcome
2021, NeuroImage: ClinicalThe ictal EEG in ECT: A systematic review of the relationships between ictal features, ECT technique, seizure threshold and outcomes
2020, Brain StimulationCitation Excerpt :This view also seems consistent with the differences in both the magnitude and distribution of stimulation when comparing bilateral and unilateral electrode placements, or brief and ultra-brief pulse width sessions, respectively. Such differences have been demonstrated most clearly in recent computational modelling studies - the fact that ictal features appear to “saturate” at a certain intensity above threshold (and that this ceiling point will most likely vary between forms of ECT) appears to reflect the resulting field of stimulation [65]. While there are limitations when attempting to distinguish treatment at different intensities between patients, two studies examined suggested monitoring the ictal EEG within a patient may be particularly relevant for detecting rises in seizure threshold during a course of ECT, with clear implications for dosing decisions [20,53].