PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of nihpaAbout Author manuscriptsSubmit a manuscriptHHS Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
 
Neuroimage. Author manuscript; available in PMC 2013 February 1.
Published in final edited form as:
PMCID: PMC3495594
NIHMSID: NIHMS332538

Regional electric field induced by electroconvulsive therapy in a realistic finite element head model: Influence of white matter anisotropic conductivity

Abstract

We present the first computational study investigating the electric field (E-field) strength generated by various electroconvulsive therapy (ECT) electrode configurations in specific brain regions of interest (ROIs) that have putative roles in the therapeutic action and/or adverse side effects of ECT. This study also characterizes the impact of the white matter (WM) conductivity anisotropy on the E-field distribution. A finite element head model incorporating tissue heterogeneity and WM anisotropic conductivity was constructed based on structural magnetic resonance imaging (MRI) and diffusion tensor MRI data. We computed the spatial E-field distributions generated by three standard ECT electrode placements including bilateral (BL), bifrontal (BF), and right unilateral (RUL) and an investigational electrode configuration for focal electrically administered seizure therapy (FEAST). The key results are that (1) the median E-field strength over the whole brain is 3.9, 1.5, 2.3, and 2.6 V/cm for the BL, BF, RUL, and FEAST electrode configurations, respectively, which coupled with the broad spread of the BL E-field suggests a biophysical basis for observations of superior efficacy of BL ECT compared to BF and RUL ECT; (2) in the hippocampi, BL ECT produces a median E-field of 4.8 V/cm that is 1.5–2.8 times stronger than that for the other electrode configurations, consistent with the more pronounced amnestic effects of BL ECT; and (3) neglecting the WM conductivity anisotropy results in E-field strength error up to 18% overall and up to 39% in specific ROIs, motivating the inclusion of the WM conductivity anisotropy in accurate head models. This computational study demonstrates how the realistic finite element head model incorporating tissue conductivity anisotropy provides quantitative insight into the biophysics of ECT, which may shed light on the differential clinical outcomes seen with various forms of ECT, and may guide the development of novel stimulation paradigms with improved risk/benefit ratio.

Keywords: Electroconvulsive therapy, Finite element method, Computational model, White matter anisotropy, Electric field, Magnetic resonance imaging, Brain stimulation

1. 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, 2011; Weaver et al., 1976) to low-resolution realistically-shaped representations (Bai et al., 2011; Sekino and Ueno, 2002, 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., 2009; Datta et al., 2008; Grandori and Rossini, 1988; Holdefer et al., 2006; Im et al., 2008; Lee et al., 2009b; Miranda et al., 2007; Miranda et al., 2006; 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., 2010; 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., 2009; Hallez et al., 2008; 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, 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, 2001). However, Rullmann et al. (2009) and 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., 2009; Hallez et al., 2008). 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).

2. Materials and methods

The steps of the E-field modeling and analysis are diagrammed in Fig. 1 and described below.

Fig. 1
Schematic illustration of the methods for generating a realistic finite element (FE) head model incorporating white matter (WM) anisotropic conductivity for E-field simulation and region of interest (ROI) analysis of ECT. T1-weighted MRI and diffusion-weighted ...

2.1. MRI and DT-MRI data acquisition and preprocessing

To construct an individual volume conductor model of the head, T1-weighted MRI and diffusion weighted MRI data sets of a healthy subject (28 year old male) were acquired with a 3T MRI scanner (Magnum 3.0, Medinus Inc., Republic of Korea). The MRI images captured the head above the level of the auditory canal, including the complete brain (see Fig. 2(a) and Fig. 3(top row)). The T1-weighted structural MRI data were obtained using a standard anatomical MR imaging sequence (TR = 35 ms; TE = 7 ms; 180 slices; 1×1×1 mm3 voxel; FOV = 256 mm). The diffusion weighted MRI data were acquired using a single-shot spin-echo echo-planar imaging sequence (TR = 8280 ms; TE = 70 ms; 70 contiguous slices; 1.75×1.75×2 mm3 voxel; 2 averages). The diffusion sensitizing gradients with a b-value of 600 s/mm2 were applied in 45 non-collinear directions. The diffusion weighted images were corrected to remove eddy current and subject motion artifacts. We then computed the diffusion tensor for each voxel of the DT-MRI data set based upon a mono-exponential relationship between the signal attenuation and the diffusion tensor matrix (Basser et al., 1994a). We also extracted the region corresponding to the subject’s head from the background noise and artifacts in the T1-weighted MRI data set by applying a morphological processing technique that included opening and closing of the head binary masks.

Fig. 2
(a) Partial volume rendering of the human head model. The cropped section shows the five segmented tissue compartments. (b) A transaxial conductivity map with the principal orientations (the largest eigenvectors) of the WM conductivity tensors projected ...
Fig. 3
Cut-away 3D rendering of the head model (top row) and the E-field magnitude spatial distribution in the anisotropic head model for BL, BF, RUL, and FEAST electrode configurations (second to bottom rows, respectively) with 800 mA current. Columns from ...

2.2. Tissue segmentation

For FE volume conductor modeling, the T1-weighted MRI was first coregistered to the DT-MRI using SPM5 (http://www.fil.ion.ucl.ac.uk/spm/) in order to avoid the transformation of the diffusion tensors (Abascal et al., 2008; De Lucia et al., 2007; Hallez et al., 2009; Hallez et al., 2008). In this process, a voxel similarity-based registration technique was applied to acquire coregistered T1-weighted MRI by maximizing mutual information that measures the degree of mutual dependence between the image intensities of corresponding voxels in both images (Wells et al., 1996). To generate the FE meshes from the coregistered structural MRI images, the MRI images were segmented into five different sub-regions including scalp, skull, CSF, gray matter, and WM. BrainSuite2 (Shattuck and Leahy, 2002) was used to extract brain tissue compartments (gray matter and WM) as well as CSF regions including the ventricles. We then segmented the skull and scalp regions using a skull extraction algorithm (Dogdas et al., 2005) based on a combination of thresholding and mathematical morphological operations including opening and closing.

2.3. ECT electrode configurations

We added realistically-shaped ECT electrodes to the 3D head model by intersecting cylindrical or rectangular solid geometries with the head rendering. We modeled three conventional electrode placements (BL, BF, and RUL) (Abrams, 2002; APA, 2001) and an investigational configuration (FEAST) (Sackeim, 2004; Spellman et al., 2009) diagrammed in Fig. 1. Standard round electrodes (5 cm diameter) were modeled for the BL, BF, and RUL electrode configurations. For BL ECT, the two electrodes were centered bilaterally at the frontotemporal positions located 2.5 cm above the midpoint of the line connecting the external canthus and tragus (APA, 2001). For BF ECT, the electrodes were placed bilaterally 5 cm above the outer angle of the orbit on a line parallel to the sagittal plane (Abrams, 2002). For RUL ECT, one electrode was centered 2.5 cm to the right of vertex and the second electrode was placed in the homologous right frontotemporal position (APA, 2001). For FEAST, a wide rectangular electrode (2.5 cm × 6.3 cm) was placed over the right motor strip and a small circular electrode (2 cm diameter) was placed over the right eyebrow. Since the conductivity of the steel electrodes is more than five orders of magnitude higher than that of the scalp (see Sec. 2.5), the electrode surface is effectively equipotential for any practical electrode thickness. Therefore, we did not accurately model the electrode thickness but, for simplicity, kept the outer surface of the electrodes flat, resulting in the electrode thickness between the flat outer surface and the scalp curvature varying from 5 mm to 12 mm. The electrode wires were modeled as current sources respectively sourcing and sinking current in the centers of the outer electrode surface, reflecting the wiring of conventional ECT electrode paddles, even though the high conductivity and the resulting equipotential surface of the electrodes make inconsequential the exact point of injection of current from the wires to the electrodes. Finally, we did not model the electrolyte gel applied to the ECT electrode surface, since the gel’s purpose is to stabilize the impedance of the electrode–scalp interface, which is already assumed in the computational model, and since the gel is largely displaced when the electrodes are pressed against the head.

2.4. Finite element meshing

For the FE mesh generation, we utilized the Computer Geometry Algorithm Library (CGAL) (http://www.cgal.org). The mesh generator is based on the labeled voxel-volume meshing technique (Boissonnat and Oudot, 2005) and allows generation of FE tetrahedral meshes which contain one sub-mesh for each sub-domain and surface meshes that approximate the boundaries of the domain and sub-domain. The triangulation algorithm provides a discretized approximation of tissue compartments and their surface boundaries according to the restricted Delaunay triangulation paradigm (Boissonnat and Oudot, 2005; Pons et al., 2007), resulting in 3D meshes of each tissue domain and conformal surface meshes for all tissue boundaries and subdividing surfaces. The resulting FE mesh of the human head and the ECT electrodes consists of approximately 1.6 million tetrahedral elements.

2.5. Electrical conductivity assignment

To evaluate the effect of WM anisotropy on the E-field, we simulated the E-field in a model with WM anisotropic conductivity as well as in a model with fully isotropic conductivity. In the isotropic head model, the electrical conductivities were assigned to 9.8×105 S/m for the steel electrodes, 0.33 S/m for scalp, 0.0132 S/m for skull, 1.79 S/m for CSF, 0.33 S/m for gray matter, and 0.14 S/m (nominal) for WM (Awada et al., 1998; Kim et al., 2003; Wolters et al., 2006). In the anisotropic model, all tissues except WM were assumed to be isotropic, since the WM has the most significant anisotropic microstructure (Geddes and Baker, 1967; Nicholson, 1965). Another tissue that is often treated as anisotropic is the skull (Fuchs et al., 2007; Marin et al., 1998; Rush and Driscoll, 1968; Suh et al., 2009; Suh et al., 2010). However, the effective skull anisotropy originates from the macroscopic skull structure consisting of a soft (spongiform) bone layer enclosed by two hard (compact) bone layers (Akhtari et al., 2002; Fuchs et al., 2007). Therefore, the most accurate approach is to model the skull as three discrete isotropic layers (Holdefer et al., 2006; Rampersad et al., 2011; Sadleir and Argibay, 2007). In the present study, however, we modeled the skull as a single isotropic layer, since the low MRI signal from bone makes accurate skull layer segmentation difficult (Chen and Mogul, 2009), and since the mostly radial flow of current in the skull justifies the compound three-layer conductivity to be approximated by an effective radial conductivity value (Rampersad et al., 2011).

To estimate the WM anisotropic conductivity tensors, we used the assumption that the conductivity tensors share eigenvectors with the measured diffusion tensors (Basser et al., 1994b). We deployed the volume constraint approach to estimate the WM anisotropic conductivity with a fixed anisotropy ratio in each WM voxel (Wolters et al. 2006). The WM anisotropic conductivity tensor σWM was modeled to be prolate

equation M1
(1)

where S denotes the orthogonal matrix of unit length eigenvectors of the measured DTs at the barycenter of the WM tetrahedral elements. Parameters σlong and σtrans are the conductivity eigenvalues longitudinal (parallel) and transverse (perpendicular) to the WM fiber direction, respectively, with σlong ≥σtrans. We computed σlong and σtrans from the WM isotropic conductivity value of σiso = 0.14 S/m and the anisotropic factor k (Gullmar et al., 2010; Wolters et al., 2006)

equation M2
(2)

equation M3
(3)

In the present study, we used a fixed anisotropy ratio of 10:1 (k = 10) (Nicholson, 1965), yielding σlong = 0.65 S/m and σtrans = 0.065 S/m.

To investigate how well isotropic models approximate the anisotropic model, we also simulated models with isotropic WM conductivity ranging from the transverse (low) conductivity estimate, σiso = 0.065 S/m, to the longitudinal (high) conductivity estimate, σiso = 0.65 S/m, including the nominal volume-constraint value of σiso = 0.14 S/m.

2.6. E-field computation

The E-field distribution was computed with the FE method software ANSYS (ANSYS Inc., PA, USA). For all electrode configurations the current was set to 800 mA, corresponding to the conventional setting used with the MECTA Spectrum 5000Q ECT device (MECTA Corp., OR, USA). Since the current waveform frequencies in ECT are relatively low (< 10 kHz), the E-field was obtained by solving the quasi-static Laplace equation with no internal sources (Bossetti et al., 2008; Plonsey and Heppner, 1967; Rattay, 1988)

equation M4
(4)

where Φ and σ are the electrical potential and the tissue conductivity tensor, respectively. The Neumann boundary conditions apply on the surface of the model

equation M5
(5)

where [n with circumflex] is the unit vector normal to the outer surface of the model, except at the outer centers of the electrodes where current is injected. For each of the electrode configurations, the system of linear equations of the FE method was solved using the preconditioned conjugate gradient solver with a relative tolerance of 10−8. Finally, the E-field distribution was computed by taking the gradient of the scalar potential Φ.

2.7. E-field analysis

The E-field was sampled in specific brain ROIs thought to be associated with therapeutic action of ECT, including frontal pole, orbitofrontal cortex (OFC), dorsolateral prefrontal cortex (DLPFC), thalamus, hypothalamus, and subcallosal cingulate cortex (SCC) (Abrams, 2002; Abrams and Taylor, 1976; Krystal and Weiner, 1994; Krystal et al., 1993; Mayberg, 2009; Mayberg et al., 2005; Sackeim, 2004), or with side effects of ECT, including hippocampus and insula (Kellner et al., 2010a; Lisanby et al., 2000; Sackeim et al., 2008; Squire et al., 1981). We also examined the E-field in primary motor cortex (specifically, the first dorsal interosseous (FDI) motor area) and brainstem, which are thought to be relevant to seizure initiation and motor manifestations of the seizure, respectively (Blumenfeld et al., 2003; Sackeim and Mukherjee, 1986). These ten anatomically defined ROIs were manually segmented from coronal MRI sections, based on human brain atlases and definitions in the literature (Goldstein et al., 2007; John et al., 2007; Makris et al., 1999; McCormick et al., 2006; Nakamura et al., 2008; Ranta et al., 2009; Weiss et al., 2005; Wible et al., 1997), and were verified with the BrainParser software (Tu et al., 2008). We then computed descriptive statistics of the E-field strength (median, 25th and 75th percentiles, and minimum and maximum) in each of these ROIs as well as in the whole left and right hemispheres for the BL, BF, RUL, and FEAST electrode configurations.

We assessed the effect of the WM anisotropic conductivity on the induced E-field qualitatively by plotting maps of the E-field magnitude and current density vectors, and quantitatively by calculating the relative difference in E-field magnitude between the two solutions in various ROIs. The difference between the isotropic and anisotropic solutions was quantified by the statistical measure of relative error, defined as

equation M6
(6)

where n is the number of samples in the respective ROI, and Eiso and Eaniso denote the E-field magnitude in the isotropic and anisotropic models, respectively (Drechsler et al., 2009; Kim et al., 2003; Klepfer et al., 1997; Lee et al., 2006; Lee et al., 2009b).

3. Results

3.1. 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 conductivity map framed in white in Fig. 2(b) is magnified in Fig. 2(c).

3.2. Comparison of ECT electrode configurations

Fig. 3 shows a set of results for the spatial E-field magnitude distribution in the head model incorporating a fixed WM anisotropy ratio of 10:1. A qualitative comparison of the spatial E-field distributions in Fig. 3 indicates that the different ECT electrode configurations result in distinct E-field distributions in the brain. As expected, the strongest E-field tends to occur in the brain volume under and between the two electrodes, although the detailed distribution is complex and depends heavily on the head anatomy. The symmetric ECT electrode configurations (BL and BF) have comparable E-field strength in both hemispheres, whereas the lateralized configurations (RUL and FEAST) generate higher E-field magnitude in the right hemisphere. Similarly, configurations with anterior electrodes (BF and FEAST) induce stronger E-fields in the anterior portions of the brain than the other configurations with more posterior electrodes (BL and RUL).

Fig. 4 shows descriptive statistics of the E-field strength (median, 25th and 75th percentiles, and minimum and maximum) in various ROIs in the left and right hemispheres of the anisotropic head model for each ECT electrode configuration. The median E-field strength in the whole brain is 3.9, 1.5, 2.3, and 2.6 V/cm, and the right-to-left hemisphere median E-field ratio is 1.0, 1.0, 2.1, and 1.9 for the BL, BF, RUL, and FEAST configurations, respectively. Besides generating the strongest E-field overall, BL ECT produces the highest median E-field in the hippocampi (4.8 V/cm), which is 1.5–2.8 times stronger than the hippocampal E-field of the other electrode configurations. In comparison to BL, RUL ECT has weaker median E-fields in all regions except in the right FDI motor area, where it is 1.4 times stronger. The ratio of the frontal (frontal pole, OFC, and DLPFC) to temporal (hippocampus) median E-field is 0.9, 3.5, 0.9, and 1.7 for the BL, BF, RUL, and FEAST configurations, respectively. Due to its frontal electrode placement, BF ECT has the highest frontal-to-temporal field ratio (3.5) among the configurations, inducing the strongest E-field in the frontal pole as well as substantial field strength in the OFC and DLPFC. Compared to the BF and RUL electrode placements, FEAST produces stronger median E-field stimulation in thalamus, hypothalamus, SCC, and insula. Moreover, compared to all of the other configurations, FEAST induces 1.2–2.0, 1.4–2.4, and 1.6–2.2 times stronger E-fields in the right OFC, DLPFC, and SCC, respectively.

Fig. 4
Descriptive statistics of the regional E-field magnitude generated by the four ECT electrode configurations in the left and right hemispheres of the anisotropic head model. The E-field strength (y-axis) is shown on a logarithmic scale. The boxes indicate ...

3.3. Effect of white matter anisotropic conductivity

To demonstrate the effect of WM anisotropy on the E-field and current density solutions, in Fig. 5 we plotted maps of the E-field magnitude and the current density vector field in a coronal slice (same as in Fig. 3(middle column)) for BL and RUL ECT in the nominal isotropic (WM σiso = 0.14 S/m) and anisotropic head models. In this slice, the BL and RUL configurations generate, respectively, predominantly mediolateral and inferosuperior current flow. In both the isotropic and anisotropic models, the high conductivity of CSF-filled structures tends to channel current flow, as illustrated by the larger arrows in the lateral ventricles in Figs. 5(c) and (d) and in the longitudinal fissure in Figs. 5(e) and (f). Despite the higher current density in CSF-filled structures, the E-field magnitude there is relatively low (blue color) due to the high conductivity. Compared to the isotropic models (Figs. 5(c) and (e)), in the anisotropic models (Figs. 5(d) and (f)) additional channeling effects emerge as a result of the WM fiber orientation. For example, when the current flow is partially aligned with the orientation of the WM fibers, the anisotropic model enhances the current alignment with the fibers, as demonstrated by the curving of the current flow along the WM fibers in the corpus callosum and in the regions inferior to the lateral ventricles in BL ECT (Fig. 5(d) versus Fig. 5(c)). Furthermore, long stretches of fibers aligned with the current flow can result in denser current along the fibers, as demonstrated by the concentrated current flow along the inferosuperior-oriented internal capsule fibers in both hemispheres in RUL ECT (Fig. 5(f) versus Fig. 5(e)). On the other hand, current flowing in direction transverse to the WM fiber orientation generates a relatively high E-field magnitude due to the lower electrical conductivity across the WM fibers, as exemplified by the regions well demarcated in red color in the BL E-field map (Fig. 5(d) versus Fig. 5(c)) corresponding to lateromedial current flow across inferosuperior and anteroposterior oriented pyramidal tract fibers.

Fig. 5
Comparison of the isotropic and anisotropic model simulations. (a) A coronal slice conductivity map using the same display conventions as in Fig. 2. (b) Enlarged view of the region framed in white in (a). (c–f) View corresponding to (b) of the ...

Quantitatively, the E-field magnitude relative error over the whole brain of the nominal isotropic model (WM σiso = 0.14 S/m) compared to the anisotropic model is 18%, 7%, 7%, and 6% for BL, BF, RUL, and FEAST ECT, respectively. Fig. 6 breaks these data down for distinct ROIs in the left and right hemispheres. The largest ROI error (39%) is observed in the left FDI motor area for RUL ECT, and the lowest error (1%) occurs in the left frontal pole for BL ECT. In prefrontal and frontal areas including frontal pole, OFC, and DLPFC in Figs. 6(b)–(d), the maximal error of 19% occurs with RUL ECT in the left OFC. The SCC in Fig. 6(g) exhibits relatively larger differences compared to other ROIs, ranging from 14% (FEAST) to 39% (RUL). In hippocampus and insula, the largest errors are 25% and 34%, respectively.

Fig. 6
Relative error of the E-field magnitude in the isotropic versus the anisotropic head model for the various brain ROIs and ECT electrode configurations. Relative error is defined in equation (6).

Fig. 7 shows the relative error in E-field magnitude over the whole brain between the isotropic and anisotropic models as a function of the WM conductivity in the isotropic model. The largest errors occur with the highest isotropic conductivity of 0.65 S/m for all configurations. The lowest errors, indicating the “best” isotropic model match to the anisotropic model, are at σiso = 0.11 S/m for BL and BF ECT and at σiso = 0.18 S/m for RUL and FEAST ECT, which are close to but distinct from the volume-constraint value of σiso = 0.14 S/m.

Fig. 7
E-field magnitude relative error in the whole brain between the anisotropic head model (WM σlong = 0.65 S/m and σtrans = 0.065 S/m) and the isotropic model with WM conductivity ranging from σiso = 0.065 S/m to 0.65 S/m. The four ...

4. 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 anatomically realistic head model.

4.1. Implications for ECT technique

It has been suggested that the flow of electric current and a resultant robust seizure expression in the prefrontal cortex is requisite to the antidepressant effects of ECT, as demonstrated, for example, by the superior efficacy of BL ECT compared to low dose RUL ECT (Abrams and Taylor, 1976; Ottosson, 1960; Sackeim, 2004). This view has motivated electrode configurations designed to target the prefrontal cortex, such as BF and FEAST. The results from our model confirm that BL ECT generates stronger E-field in anterior portions of the brain compared to RUL ECT. Our results further indicate that BF and FEAST generally produce higher E-field strengths in prefrontal structures compared to BL and RUL ECT (see Figs. 3 and and4).4). In particular, compared to the other electrode configurations, BF produces the highest E-field magnitude in the frontal poles. Furthermore, FEAST produces the strongest E-field in the right OFC, DLPFC, and SCC—likely a consequence of current flow through the orbits which represent low impedance paths into the brain. Significantly, DLPFC and SCC are targets for other forms of brain stimulation with reported antidepressant efficacy. The right DLPFC has been targeted with low-frequency (inhibitory) repetitive TMS for the treatment of depression (Fitzgerald and Daskalakis, 2011). There is evidence that focal stimulation of the SCC via deep brain stimulation reduces elevated SCC activity and normalizes aberrant network activity in depression, with resultant antidepressant effect (Mayberg et al., 2005; McNeely et al., 2008). The frontal pole and OFC have strong anatomical and functional connectivity with the SCC and the rest of the dysfunctional brain network associated with depression (Gutman et al., 2009; Johansen-Berg et al., 2008; Mayberg, 2009; Mayberg et al., 2005) and therefore are potential targets for stimulation as well.

Another prevailing hypothesis on the mechanism of action of ECT focuses on diencephalic stimulation. This hypothesis states that for optimal antidepressant efficacy, the seizure must be sufficiently generalized to involve diencephalon centers, particularly the thalamus and hypothalamus (Abrams, 2002; Abrams and Taylor, 1976). Cortical–thalamocortical interactions are crucial for the initiation, propagation, and behavioral manifestations of generalized seizures (Guillery and Sherman, 2002; McNally and Blumenfeld, 2004). Thalamic processes that inhibit cortical function have been hypothesized to play a role in the antidepressant effect of ECT (Krystal and Weiner, 1994; Krystal et al., 1993; McNally and Blumenfeld, 2004; Sackeim et al., 1991). Our results show that BL ECT produces stronger E-field in the thalamus and hypothalamus compared to BF and RUL ECT, consistent with the superior antidepressant efficacy of BL ECT as compared with low dose RUL ECT, but does not explain how adequately-dosed RUL ECT matches the efficacy of BL ECT. Notably, FEAST produces E-field strengths in the diencephalon centers that are higher than in BF and RUL, and approach those in BL. However, whether this confers high antidepressant efficacy to FEAST is unknown since FEAST clinical data has not been reported.

It has been suggested that retrograde amnesia may be related to seizure activity in the medial temporal lobe (Lisanby et al., 2000; Luber et al., 2000; Squire, 1986) and consequently it was proposed that techniques inducing seizures with reduced current spread to this area may cause less memory side effects (Sackeim, 2004). The medial temporal structure specifically associated with amnesia is the hippocampus, which has been shown to be uniquely affected by induced seizures, and to have a low seizure threshold (Lisanby et al., 2000). Of all configurations in our model, BL ECT produces the strongest E-field in the hippocampi, consistent with its greater acute cognitive impairment (Kellner et al., 2010a; Sackeim et al., 2008) and short- and long-term retrograde amnesia (Lisanby et al., 2000; Squire et al., 1981) compared to BF and RUL ECT. Our results may be particularly useful for predicting the clinical effects of ECT paradigms for which clinical data is not yet available, such as FEAST. For example, the E-field strength generated by FEAST is lower than RUL in the right hippocampus, but is higher than RUL in the left hippocampus, and is higher compared to BF in both hippocampi. Therefore, one might predict that FEAST may have memory side effects lower than in BL but higher than in BF (assuming identical current levels in all configurations); this is yet to be tested clinically.

The E-field model could also inform the mechanisms behind cardiac effects of ECT. Accumulating evidence indicates that cardiovascular regulation receives significant input from cortical structures, especially the insula (Oppenheimer, 2006). Direct electrical stimulation of the left caudal anterior insula leads to increased occurrence of bradycardia and depressor responses as a result of parasympathetic activation, whereas right anterior insular stimulation leads to tachycardia and diastolic blood pressure elevation as a result of sympathetic activation (Oppenheimer, 2006; Oppenheimer et al., 1992). Data from patients with insular lesions are consistent with left insular cardioinhibitory representation and right insular cardioexcitatory representation (Oppenheimer, 2006). During the stimulus delivery phase of ECT, a higher likelihood of heart rate reduction and longer duration of asystole were observed with RUL compared to BL and BF ECT (Nagler, 2010). Since RUL ECT produces nearly three times higher E-field strength in the right insula compared to the left insula (Fig. 4(i)), we hypothesize that the right insula is preferentially inhibited during stimulation compared to the left insula, resulting in sympathetic withdrawal. In contrast, BL and BF ECT induce similar E-field strengths in left and right insula, thus the sympathovagal balance is less affected, consistent with the lower rates of asystole observed clinically (Loo et al., 2009; Nagler, 2010). Other brain regions, such as the hypothalamus and brainstem, are also known to be involved in cardiac regulation (Welch and Drop, 1989); therefore, the E-field characteristics in these regions (Figs. 4(f) and (k)) may also contribute to cardiac effects.

Our results indicate that the E-field strengths induced by RUL are exclusively weaker than in BL ECT except in the right motor strip where the E-field magnitude induced by RUL ECT is about 38% stronger compared to BL. This observation is consistent with the lower seizure threshold of RUL compared to BL ECT, since the motor strip is thought be the likely site for seizure initiation (Sackeim and Mukherjee, 1986). We also found that BL and FEAST generate higher E-field strengths compared to BF and RUL in the brainstem which mediates the motor manifestations of generalized tonic-clonic seizures (Blumenfeld et al., 2003).

The above observations illustrate how the realistic ECT E-field model can contribute to a biophysical explanation of reported clinical differences among conventional electrode placements (BL, BF, and RUL), as well as to the evaluation and optimization of investigational configurations (e.g., FEAST). For example, the realistic head model can be coupled with an optimization algorithm to select scalp electrode locations and current strengths to target specific brain structures (Dmochowski et al., 2011; Park et al., 2011). Ultimately, this work may guide the development of novel stimulation paradigms with improved risk/benefit ratio.

4.2. Influence of tissue conductivity anisotropy

The second objective of this study was to evaluate the influence of WM conductivity anisotropy on the E-field simulation. Our results, summarized in Fig. 5 and Fig. 6, indicate that neglecting WM conductivity anisotropy and using instead the volume-constraint WM isotropic conductivity value of 0.14 S/m leads to relative errors in the E-field magnitude up to 18% for the whole brain and up to 39% within the considered ROIs. The maximum relative error was found in the left FDI motor area for RUL ECT. In addition, our results indicate that the E-field differences between the isotropic and anisotropic models depend upon the specific ECT electrode configuration used. For example, ignoring WM anisotropy produces errors in the FDI motor area E-field up to only 8% for the BF configuration, but up to 39% for the RUL configuration. Furthermore, the WM isotropic conductivity value that gives the best approximation of the anisotropic case, in terms of the lowest overall relative error, varies with electrode configuration, as shown in Fig. 7. Specifically, the optimal WM isotropic conductivity value for the lateralized configurations (0.18 S/m for RUL and FEAST) is 64% higher than that for the bilaterally symmetric configurations (0.11 S/m for BL and BF). A possible explanation for this difference is that in the lateralized electrode configurations the electrical current is more aligned, on average, with the orientation of the WM fibers than in the bilaterally symmetric configurations. The results in Fig. 7 also indicate that the volume-constraint WM isotropic conductivity (0.14 S/m) yields relative errors within 0.7% of the optima achievable with the isotropic models and is, therefore, a reasonable choice for isotropic models. Generally, the E-field magnitude errors and their sensitivity to the electrode configuration motivate the inclusion of the WM conductivity anisotropy in computational E-field models for increased accuracy.

Another important observation is that the error in the E-field strength between the isotropic and anisotropic models usually increases for brain regions that are farther away from the ECT electrodes. For example, the E-field errors in deep brain structures such as thalamus, hypothalamus, insula, and SCC (10%–39%) are higher than the errors in more superficial areas such as frontal pole and DLPFC (1%–12%) which lie closer to the scalp electrodes. Furthermore, although the overall relative error in the whole brain is higher for BL than for BF, the relative error is larger for BF than for BL in ROIs such as hypothalamus, hippocampus, insula, and brainstem that are farther away from the BF electrodes. Similarly, in the lateralized electrode configurations (RUL and FEAST), the E-field errors for highly lateralized ROIs (whole hemisphere, DLPFC, insula, and FDI motor area) are significantly larger on the left side than on the right side (where the electrodes are placed). The E-field error increase away from the electrodes could be explained in terms of the longer paths that the electrical current has to traverse from the electrodes to distant brain regions, which results in a larger cumulative error of all the differences in the conductivity tensors along the current path. Thus, incorporation of WM conductivity anisotropy in ECT E-field models is crucial for analysis of the E-field characteristics in brain regions that we try to avoid stimulating by placing the electrodes away from them. Often these are brain regions thought to be associated with adverse side effects of ECT, and thus the degree to which they are stimulated is of particular relevance to studies that evaluate ECT techniques aimed at improving safety. This observation further supports the inclusion of WM anisotropic conductivity in ECT models.

The comparison of the isotropic and anisotropic model simulations in Fig. 5 provides insight into how the orientation of WM fibers affects the current density and E-field distributions. The current flow generally follows the path of least electrical impedance. Consequently, compared to the isotropic model, current flow along the WM fibers in the anisotropic model is denser and follows the fiber orientation. In some cases the channeling of current along WM fiber tracts increases the local E-field strength, but in other cases this effect is offset by the lower impedance along the fibers, leading to reduced E-field strength. The current may steer away from segments of increased impedance resulting from the passage of WM fibers perpendicular to the current flow. However, when current comes across wide stretches of fibers perpendicular to its flow, the current cannot steer away and is forced to cross this high impedance barrier, resulting in high E-field magnitude. Thus, increased E-field strength in the anisotropic model compared to the isotropic model can result from either concentration of current flow along WM fibers that is not offset by the low impedance along the fibers, or from increased impedance for current flow perpendicular to WM fibers.

In support of our findings, other bioelectric head modeling studies have reported comparable effect sizes associated with the inclusion of WM conductivity anisotropy. We previously studied the impact of WM anisotropy on the current density distribution generated by tDCS and found relative errors of 53% and 19% in current density magnitude in the WM and gray matter, respectively (Lee et al., 2009b). Similarly, Sadleir et al. (2010) concluded that the inclusion of the WM anisotropic conductivity in a tDCS model would result in differences up to 39% in the median current density magnitude. De Lucia et al. (2007) found that WM anisotropy contributes a difference of 10% in the peak TMS-induced E-field, which is consistent with the findings of Thielscher et al. (2011). (It should be noted that the TMS E-field is substantially more superficial than the E-field in ECT and tDCS; therefore, the lower error associated with WM anisotropic conductivity in TMS is expected). Finally, Wolters et al. (2006) reported that the inclusion of WM anisotropic conductivity results in up to 30% difference in the magnitude of electroencephalographic scalp potentials generated by dipole sources within the brain. Thus, the differences between the isotropic and anisotropic model results in these studies are comparable to the relative errors up to 39% found in our study.

It has been reported that depression itself is associated with changes in the volume of specific brain structures (Bremner, 2005; Tsopelas et al., 2011) and with regionally specific abnormalities of the WM fractional anisotropy (Korgaonkar et al., 2010; Wu et al., 2011). Patients with major depressive disorder showed a fractional anisotropy reduction up to 14% in WM, with an 8% decrease in the superior longitudinal fasciculus associated with the DLPFC (Korgaonkar et al., 2010). Neglecting WM anisotropy in our model, which is associated with fractional anisotropy reduction of 100%, results in up to 39% difference in E-field strength. Extrapolating from these data, a pathological 14% decrease in fractional anisotropy could result in approximately 6% difference in the E-field strength. Even though this effect is relatively small, tissue conductivity anisotropy and other pathological brain structure changes should be considered as a potential source of E-field variability, as their compounded effect may be significant.

4.3. Model validity and limitations

Validation of brain stimulation E-field simulations remains a challenging problem due to the unavailability of methods for high-resolution in vivo E-field measurements. Nevertheless, our results in the brain (E-field strength median of 1.5–3.9 V/cm and range of 0.1–300 V/cm) are in good agreement with published modeling and experimental measurements. In the references to other studies below, the electrode current was scaled to 800 mA to allow comparison with our data. In their anatomically-realistic computational model of BL ECT, Nadeem et al. (2003) did not report statistics on the E-field distribution, but the figures suggest a brain E-field magnitude ranging from approximately 0.1 V/cm to as high as 700 V/cm. Rescaling the current from Sadleir et al.’s (2010) and Parazzini et al.’s (2011) anatomically-realistic tDCS models and converting current density to E-field strength, the average of their reported median values is 2.5 V/cm and 2.7 V/cm, respectively. Studies using simplified model geometries or low-resolution experimental measurement techniques generally reported lower E-field strengths (Deng et al., 2011; Lorimer et al., 1949; Rush and Driscoll, 1968; Smitt and Wegener, 1944). For instance, a recent simulation study by our group of BL, BF, RUL, and FEAST ECT in a spherical head model obtained median and maximum brain E-field strength ranges of 0.6–1.2 V/cm and 2.1–2.5 V/cm, respectively (Deng et al., 2011). Rush and Driscoll’s (1968) measurements of a frontal-occipital electrode configuration in a half-skull immersed in an electrolytic tank yielded maximum E-field in the range of 1.5–2.5 V/cm. Intracerebral measurements in cadaver heads with stimulation current applied through bifrontotemporal electrodes produced maximum E-field strength estimates of 0.7–1.8 V/cm (Lorimer et al., 1949; Smitt and Wegener, 1944). A likely factor contributing to the lower E-field strengths reported in these studies is the effective averaging out of the spatial distribution of the E-field resulting from the simplified head models and/or from the low-resolution spatial sampling of the electric potentials in the cadaver measurements. Indeed, a study comparing head models for TMS found a ~ 51% increase of the maximum E-field strength in the anatomically-realistic model compared to simplified spherical models (Thielscher et al., 2011).

A limitation of the present study that may impact the E-field strength and distribution is the truncation of the head model. The T1-weighted MRI data was only acquired for the portion of the head above the level of the auditory canal. The truncated head model eliminates shunting of the ECT stimulus current in the lower portion of the head, resulting in increased E-field magnitude in the brain. Further, the head model truncation has differential effects for various ECT electrode placements. For example, in our ECT analysis using a spherical head model (data not shown), we observed that truncation of the head model results in a 52% and 27% increase in the median E-field and a 41% and 14% increase in the peak E-field for BL and RUL ECT, respectively. Therefore, future ECT head models should be based on structural MRI and DT-MRI data set of the whole head including the skull base and a portion of the neck underneath.

Uncertainty about the thickness, structure, and conductivity of the various tissues in the model can contribute to inaccuracies of the simulated E-field and to discrepancies among various models and experimental measurements (Bikson and Datta, 2011), but it is difficult to assess the extent of these uncertainties since, as discussed in the beginning of this section, there are no adequate empirical data to compare simulations to. For example, we observed overall comparable E-field strength in the left and right hemispheres for the symmetric ECT electrode configurations (BL and BF). However, certain regions such as hypothalamus, SCC, and FDI motor area produced asymmetric results of the median E-field magnitude (see Fig. 4). This may be due to intrinsic anatomical asymmetry and/or errors in the tissue segmentation and the ROI boundary definition between the two sides of the head, but it is difficult to determine how much each of these factors is contributing.

There are two distinct sources of uncertainty that confound the conclusions of modeling studies: first, naturally occurring anatomical variability in the population and, second, errors in the tissue segmentation and tissue conductivity within the modeled individual(s). We have previously investigated the effect on the induced E-field in a spherical head model of varying the thickness and conductivity of various tissue layers within ranges reported in the literature (Deng et al., 2009). That study reported that, for example, individual male scalp and skull thickness variation, which can be as high as 58% and 34% of the average, resulted in up to 76% and 20% changes of the peak brain E-field, respectively (Deng et al., 2009). The present study does not account for anatomical variability in the population since it is based on imaging data of a single individual; this limitation has to be considered when applying these results to interpret clinical ECT data obtained from various individuals. Nevertheless, the methods presented in this paper are general and can be applied to any individual with appropriate MRI data.

The second source of uncertainty is modeling errors of the tissue structure and electrical properties within an individual model. For example, the average thickness of the scalp, skull, and CSF underlying the skull are about 5.5 mm, 7 mm, and 3 mm, respectively (Deng et al., 2009), and the T1-weighted structural MRI data used to create our model has spatial resolution of 1 mm. Therefore, the MRI resolution contributes potential error of approximately 18%, 14%, and 33% to the thickness of the scalp, skull, and CSF, corresponding to estimated errors in the peak brain E-field up to 24%, 11%, and 14%, respectively, based on our perturbation data in the spherical male head model (Deng et al., 2009). Further, even though a single skull layer is a reasonable choice for spherical models, the nonuniformity of the compact and spongiform layers in a real skull motivates the segmentation of these layers in a high-accuracy model (Rampersad et al., 2011), which may require coregistering the MRI scan with a computed tomography scan that provides substantially stronger signal from bone (Chen and Mogul, 2009). There is also a wide spread of tissue conductivity values reported in the literature, especially for the skull (Deng et al., 2009), and it is not known how much of this variability is attributable to individual variation and how much—to measurement error. For example, in this study we used skull conductivity value of 0.0132 S/m whereas other studies have used lower values, e.g., 0.0083 S/m (Deng et al., 2011). Compared to a simulation of our realistic head model with the lower skull conductivity of 0.0083 S/m (data not shown), the nominal model with skull conductivity of 0.0132 S/m increases the median E-field magnitude in the brain by 13%, 10%, 16%, and 16% for BL, BF, RUL, and FEAST ECT, respectively. These estimates of E-field variability due to potential tissue thickness and conductivity errors are commensurate with the E-field relative errors up to 18% overall and up to 39% in specific ROIs that result from neglecting WM anisotropy in the present study. Therefore, it could be argued that accounting for WM anisotropy is as important to E-field strength estimation as accurate tissue segmentation and conductivity assignment. Beyond that, anisotropic models may give more accurate E-field directionality information, as suggested by Fig. 5, which may be particularly relevant if the E-field data were coupled with neural models which are direction sensitive (Radman et al., 2009; Rushton, 1927).

In our anisotropic volume conductor modeling, we adopted the volume constraint algorithm to estimate the WM anisotropic conductivity tensors with the assumption of a fixed anisotropy ratio of 10:1 in each WM voxel (Wolters et al., 2006). However, this approach may overestimate the actual ratio of the WM anisotropic conductivity tensors (Hallez et al., 2008). In reality, the ratio of longitudinal to transverse WM conductivity varies. For example, the fractional anisotropy map shown in Fig. 1 indicates that strong anisotropy is present in the pyramidal tracts and corpus callosum. On the other hand, cortical brain regions include lower degree of WM anisotropy, which is associated with lower anisotropic conductivity ratios. Recently, Bangera et al. (2010) conducted an experimental validation of anisotropic head models by measuring intracranial electric potentials generated by stimulation with an implanted dipole source in the human brain (Bangera et al., 2010). Two different anisotropic models using, respectively, the effective medium approach (Tuch et al., 1999, 2001) and the volume constraint approach with fixed 10:1 anisotropy ratio (Wolters et al., 2006) were compared revealing that the former results in a better fit to the experimental data than the latter. In the present study, we did not examine alternative approaches for estimating the anisotropic conductivity tensors such as the effective medium approach (Tuch et al., 1999, 2001) and its constrained version (Hallez et al., 2009; Hallez et al., 2008). Our model allows a relatively uncomplicated incorporation of various anisotropy estimation approaches which could be investigated in future studies.

It should be also acknowledged that the present simulations address only the E-field distribution and not the direct neural activation and the resultant seizure, the topography of which is also considered to be a major contributor to clinical outcome. At present there are no computational models that can realistically simulate the induction and propagation of seizures throughout the whole brain. Because we cannot realistically simulate the neural response to ECT, we did not explore the effect of various parameters of the ECT stimulus current such as the current amplitude, pulse width, and frequency. Nevertheless, since the E-field strength is directly proportional to the stimulus current amplitude, our data can be straightforwardly scaled to other current intensities. Furthermore, systematic data on the clinical effects of various current amplitudes is presently lacking as ECT is done exclusively with fixed current amplitudes of 800 mA or 900 mA (Peterchev et al., 2010). We have previously used the assumption of a single neural activation threshold throughout the brain to explore the effect of the stimulus current amplitude on the focality of direct neural activation in a spherical head model, and have suggested that reduction and individualization of the current may be productive strategies for better targeting of ECT (Deng et al., 2010; Deng et al., 2011; Peterchev et al., 2010). However, more empirical data linking the E-field characteristics to seizure induction have to be accumulated to support the incorporation of neural dynamics in the realistic head model, which would, in turn, enable explanation of the effect of the stimulus current parameters.

Acknowledgement

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 their helpful comments and suggestions.

References

  • Abascal JFPJ, Arridge SR, Atkinson D, Horesh R, Fabrizi L, De Lucia M, Horesh L, Bayford RH, Holder DS. 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;43:258–268. [PubMed]
  • Abrams R. Electroconvulsive therapy. 4th ed. New York: Oxford Unversity Press; 2002.
  • Abrams R, Taylor MA. Diencephalic stimulation and the effects of ECT in endogenous depression. Br. J. Psychiatry. 1976;129:482–485. [PubMed]
  • Akhtari M, Bryant HC, Marnelak AN, Flynn ER, Heller L, Shih JJ, Mandelkern M, Matlachov A, Ranken DM, Best ED, DiMauro MA, Lee RR, Sutherling WW. Conductivities of three-layer live human skull. Brain Topography. 2002;14:151–167. [PubMed]
  • APA. Washington: American Psychiatric Association; 2001. The Practice of Electroconvulsive Therapy: Recommendations for Treatment, and Privileging: A Task Force Report of the American Psychiatric Association.
  • Awada KA, Jackson DR, Baumann SB, Williams JT, Wilton DR, Fink PW, Prasky BR. Effect of conductivity uncertainties and modeling errors on EEG source localization using a 2D model. IEEE Trans. Biomed. Eng. 1998;45:1135–1145. [PubMed]
  • Bai S, Loo C, Abed AA, Dokos S. A computational model of direct brain excitation induced by electroconvulsive therapy: comparison among three conventional electrode placements. Brain Stimulation. 2011 in press. [PubMed]
  • Bangera NB, Schomer DL, Dehghani N, Ulbert I, Cash S, Papavasiliou S, Eisenberg SR, Dale AM, Halgren E. Experimental validation of the influence of white matter anisotropy on the intracranial EEG forward solution. J. Comput. Neurosci. 2010;29:371–387. [PMC free article] [PubMed]
  • Basser PJ, Mattiello J, Lebihan D. Estimation of the effective self-diffusion tensor from the NMR spin-echo. J. Magn. Reson. 1994a;103:247–254. [PubMed]
  • Basser PJ, Mattiello J, Lebihan D. MR diffusion tensor spectroscopy and imaging. Biophy. J. 1994b;66:259–267. [PubMed]
  • Bikson M, Datta A. Guidelines for precise and accurate computational models of tDCS. Brain Stimulation. 2011 in press. [PubMed]
  • Blumenfeld H, Westerveld M, Ostroff RB, Vanderhill SD, Freeman J, Necochea A, Uranga P, Tanhehco T, Smith A, Seibyl JP, Stokking R, Studholme C, Spencer SS, Zubal IG. Selective frontal, parietal, and temporal networks in generalized seizures. Neuroimage. 2003;19:1556–1566. [PubMed]
  • Boissonnat JD, Oudot S. Provably good sampling and meshing of surfaces. Graphical Models. 2005;67:405–451.
  • Bossetti CA, Birdno MJ, Grill WM. Analysis of the quasi-static approximation for calculating potentials generated by neural stimulation. J. Neural Eng. 2008;5:44–53. [PubMed]
  • Bremner JD. Changes in brain volume in major depression. Depression: Mind and Brain. 2005;2:38–46.
  • Butson CR, Cooper SE, Henderson JM, McIntyre CC. Patient-specific analysis of the volume of tissue activated during deep brain stimulation. Neuroimage. 2007;34:661–670. [PMC free article] [PubMed]
  • Chen M, Mogul DJ. A structurally detailed finite element human head model for simulation of transcranial magnetic stimulation. J. Neurosci. Methods. 2009;179:111–120. [PubMed]
  • Datta A, Bansal V, Diaz J, Patel J, Reato D, Bikson M. Gyri-precise head model of transcranial direct current stimulation: improved spatial focality using a ring electrode versus conventional rectangular pad. Brain Stimulation. 2009;2:201–207. [PMC free article] [PubMed]
  • Datta A, Elwassif M, Battaglia F, Bikson M. Transcranial current stimulation focality using disc and ring electrode configurations: FEM analysis. J. Neural Eng. 2008;5:163–174. [PubMed]
  • De Lucia M, Parker GJM, Embleton K, Newton JM, Walsh V. Diffusion tensor MRI-based estimation of the influence of brain tissue anisotropy on the effects of transcranial magnetic stimulation. Neuroimage. 2007;36:1159–1170. [PubMed]
  • Deng ZD, Lisanby SH, Peterchev AV. 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. 2009;1:682–688. [PubMed]
  • Deng ZD, Lisanby SH, Peterchev AV. Improving the focality of electroconvulsive therapy: the roles of current amplitude, and electrode size and spacing. J. ECT. 2010;26(2):151.
  • Deng ZD, Lisanby SH, Peterchev AV. Electric field strength and focality in electroconvulsive therapy and magnetic seizure therapy: a finite element simulation study. J. Neural Eng. 2011;8 016007. [PubMed]
  • Dmochowski JP, Bikson M, Datta A, Su Y, Parra LC. A multiple electrode scheme for optimal non-invasive electrical stimulation; 5th Int. IEEE/EMBS Conf. Neural Eng; 2011. pp. 29–35.
  • Dogdas B, Shattuck DW, Leahy RM. Segmentation of skull and scalp in 3-D human MRI using mathematical morphology. Human Brain Mapping. 2005;26:273–285. [PubMed]
  • Drechsler F, Wolters CH, Dierkes T, Si H, Grasedyck L. A full subtraction approach for finite element method based source analysis using constrained Delaunay tetrahedralisation. Neuroimage. 2009;46:1055–1065. [PubMed]
  • Fitzgerald PB, Daskalakis ZJ. The effects of repetitive transcranial magnetic stimulation in the treatment of depression. Expert Rev. Med. Devices. 2011;8:85–95. [PubMed]
  • Fuchs M, Wagner M, Kastner J. Development of volume conductor and source models to localize epileptic foci. J. Clin. Neurophys. 2007;24:101–119. [PubMed]
  • Geddes L, Baker L. The specific resistance of biological material. A compendium of data for the biomedical engineer and physiologist. Med.Biol.Eng. 1967;5:271–293. [PubMed]
  • Goldstein JM, Seidman LJ, Makris N, Ahern T, O'Brien LM, Caviness VS, Kennedy DN, Faraone SV, Tsuang MT. Hypothalamic abnormalities in schizophrenia: sex effects and genetic vulnerability. Biol. Psychiatry. 2007;61:935–945. [PubMed]
  • Grandori F, Rossini P. Electrical stimulation of the motor cortex: theoretical considerations. Ann. Biomed. Eng. 1988;16:639–652. [PubMed]
  • Guillery RW, Sherman SM. Thalamic relay functions and their role in corticocortical communication: generalizations from the visual system. Neuron. 2002;33:163–175. [PubMed]
  • Gullmar D, Haueisen J, Eiselt M, Giessler F, Flemming L, Anwander A, Knoesche TR, Wolters CH, Dumpelmann M, Tuch DS, Reichenbach JR. Influence of anisotropic conductivity on EEG source reconstruction: investigations in a rabbit model. IEEE Trans. Biomed. Eng. 2006;53:1841–1850. [PubMed]
  • Gullmar D, Haueisen J, Reichenbach JR. 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. 2010;51:145–163. [PubMed]
  • Gutman DA, Holtzheimer PE, Behrens TEJ, Johansen-Berg H, Mayberg HS. A tractography analysis of two deep brain stimulation white matter targets for depression. Biol. Psychiatry. 2009;65:276–282. [PubMed]
  • Hallez H, Staelens S, Lemahieu I. Dipole estimation errors due to not incorporating anisotropic conductivities in realistic head models for EEG source analysis. Phys. Med. Biol. 2009;54:6079–6093. [PubMed]
  • Hallez H, Vanrumste B, Van Hese P, Delputte S, Lemahieu I. Dipole estimation errors due to differences in modeling anisotropic conductivities in realistic head models for EEG source analysis. Phys. Med. Biol. 2008;53:1877–1894. [PubMed]
  • Haueisen J, Tuch DS, Ramon C, Schimpf PH, Wedeen VJ, George JS, Belliveau JW. The influence of brain tissue anisotropy on human EEG and MEG. Neuroimage. 2002;15:159–166. [PubMed]
  • Holdefer RN, Sadleir R, Russell MJ. Predicted current densities in the brain during transcranial electrical stimulation. Clin. Neurophys. 2006;117:1388–1397. [PMC free article] [PubMed]
  • Im CH, Jung HH, Choi JD, Lee SY, Jung KY. Determination of optimal electrode positions for transcranial direct current stimulation (tDCS) Phys. Med. Biol. 2008;53:N219–N225. [PubMed]
  • Johansen-Berg H, Gutman DA, Behrens TEJ, Matthews PM, Rushworth MFS, Katz E, Lozano AM, Mayberg HS. Anatomical connectivity of the subgenual cingulate region targeted with deep brain stimulation for treatment-resistant depression. Cerebral Cortex. 2008;18:1374–1383. [PubMed]
  • John JP, Yashavantha BS, Gado M, Veena R, Jain S, Ravishankar S, Csernansky JG. A proposal for MRI-based parcellation of the frontal pole. Brain Struct. Funct. 2007;212:245–253. [PubMed]
  • Kellner CH, Knapp R, Husain MM, Rasmussen K, Sampson S, Cullum M, McClintock SM, Tobias KG, Martino C, Mueller M, Bailine SH, Fink M, Petrides G. Bifrontal, bitemporal and right unilateral electrode placement in ECT: randomised trial. Br. J. Psychiatry. 2010a;196:226–234. [PMC free article] [PubMed]
  • Kellner CH, Tobias KG, Wiegand J. Electrode placement in electroconvulsive therapy (ECT): a review of the literature. J. ECT. 2010b;26:175–180. [PubMed]
  • Kim S, Kim TS, Zhou YX, Singh M. Influence of conductivity tensors on the scalp electrical potential: study with 2-D finite element models. IEEE Trans. Nucl. Sci. 2003;50:133–139.
  • Klepfer RN, Johnson CR, Macleod RS. The effects of inhomogeneities and anisotropies on electrocardiographic fields: A 3-D finite-element study. IEEE Trans. Biomed. Eng. 1997;44:706–719. [PubMed]
  • Korgaonkar MS, Grieve SM, Koslow SH, Gabrieli JD, Gordon E, Williams LM. Loss of white matter integrity in major depressive disorder: evidence using tract-based spatial statistical analysis of diffusion tensor imaging. Human Brain Mapping. 2010 [PubMed]
  • Krystal AD, Weiner RD. ECT seizure therapeutic adequacy. Convuls. Ther. 1994;10:153–164. [PubMed]
  • Krystal AD, Weiner RD, McCall WV, Shelp FE, Arias R, Smith P. The effects of ECT stimulus dose and electrode placement on the ictal electroencephalogram: an intraindividual crossover study. Biol. Psychiatry. 1993;34:759–767. [PubMed]
  • Lee WH, Deng ZD, Kim TS, Laine AF, Lisanby SH, Peterchev AV. Regional electric field induced by electroconvulsive therapy: A finite element simulation study. Conf. Proc. IEEE Eng. Med. Biol. Soc. 2010;1:2045–2048. [PubMed]
  • Lee WH, Deng ZD, Laine AF, Lisanby SH, Peterchev AV. Influence of white matter conductivity anisotropy on electric field induced by electroconvulsive therapy. Conf. Proc. IEEE Eng. Med. Biol. Soc. 2011;1:5473–5476. [PubMed]
  • Lee WH, Kim TS, Cho MH, Ahn YB, Lee SY. Methods and evaluations of MRI content-adaptive finite element mesh generation for bioelectromagnetic problems. Phys. Med. Biol. 2006;51:6173–6186. [PubMed]
  • Lee WH, Kim TS, Kim AT, Lee SY. 3-D diffusion tensor MRI anisotropy content-adaptive finite element head model generation for bioelectromagnetic imaging; Conf. Proc. IEEE Eng. Med. Biol. Soc; 2008. pp. 4003–4006. [PubMed]
  • Lee WH, Liu ZM, Mueller BA, Lim K, He B. Influence of white matter anisotropic conductivity on EEG source localization: comparison to fMRI in human primary visual cortex. Clin. Neurophysiol. 2009a;120:2071–2081. [PMC free article] [PubMed]
  • Lee WH, Seo HS, Kim SH, Cho MH, Lee SY, Kim TS. Influence of white matter anisotropy on the effects of transcranial direct current stimulation: a finite element study; Int. Conf. Biomed. Eng; 2009b. pp. 460–464.
  • Lisanby SH, Maddox JH, Prudic J, Devanand DP, Sackeim HA. The effects of electroconvulsive therapy on memory of autobiographical and public events. Archives of General Psychiatry. 2000;57:581–590. [PubMed]
  • Loo C, Stewart P, MacPherson R, Hadzi-Pavlovic D. Lower incidence of asystole and bradycardia with ultrabrief pulsewidth stimulation and bifrontal electrode placement. J. ECT. 2009;25:149–150.
  • Lorimer FM, Segal MM, Stein SN. Path of current distribution in brain during electro-convulsive therapy; preliminary report. Electroencephalogr. Clin. Neurophysiol. 1949;1:343–348. [PubMed]
  • Luber B, Nobler MS, Moeller JR, Katzman GP, Prudic J, Devanand DP, Dichter GS, Sackeim HA. Quantitative EEG during seizures induced by electroconvulsive therapy: relations to treatment modality and clinical features. II. Topographic analyses. J. ECT. 2000;16:229–243. [PubMed]
  • Makris N, Meyer JW, Bates JF, Yeterian EH, Kennedy DN, Caviness VS. MRI-based topographic parcellation of human cerebral white matter and nuclei II. Rationale and applications with systematics of cerebral connectivity. Neuroimage. 1999;9:18–45. [PubMed]
  • Marin G, Guerin C, Baillet S, Garnero L, Meunier G. Influence of skull anisotropy for the forward and inverse problem in EEG: simulation studies using FEM on realistic head models. Human Brain Mapping. 1998;6:250–269. [PubMed]
  • Mayberg HS. Targeted electrode-based modulation of neural circuits for depression. J. Clin. Invest. 2009;119:717–725. [PMC free article] [PubMed]
  • Mayberg HS, Lozano AM, Voon V, McNeely HE, Seminowicz D, Hamani C, Schwalb JM, Kennedy SH. Deep brain stimulation for treatment-resistant depression. Neuron. 2005;45:651–660. [PubMed]
  • McCormick LM, Ziebell S, Nopoulos P, Cassell M, Andreasen NC, Brumm M. Anterior cingulate cortex: an MRI-based parcellation method. Neuroimage. 2006;32:1167–1175. [PubMed]
  • McNally KA, Blumenfeld H. Focal network involvement in generalized seizures: new insights from electroconvulsive therapy. Epilepsy Behav. 2004;5:3–12. [PubMed]
  • McNeely HE, Mayberg HS, Lozano AM, Kennedy SH. Neuropsychological impact of Cg25 deep brain stimulation for treatment-resistant depression: preliminary results over 12 months. J. Nerv. Ment. Dis. 2008;196:405–410. [PubMed]
  • Miranda PC, Correia L, Salvador R, Basser PJ. Tissue heterogeneity as a mechanism for localized neural stimulation by applied electric fields. Phys. Med. Biol. 2007;52:5603–5617. [PubMed]
  • Miranda PC, Lomarev M, Hallett M. Modeling the current distribution during transcranial direct current stimulation. Clin. Neurophysiol. 2006;117:1623–1629. [PubMed]
  • Nadeem M, Thorlin T, Gandhi OP, Persson M. Computation of electric and magnetic stimulation in human head using the 3-D impedance method. IEEE Trans. Biomed. Eng. 2003;50:900–907. [PubMed]
  • Nagler J. Absence of asystole during bifrontal stimulation in electroconvulsive therapy. J. ECT. 2010;26:100–103. [PubMed]
  • Nakamura M, Nestor PG, Levitt JJ, Cohen AS, Kawashima T, Shenton ME, McCarley RW. Orbitofrontal volume deficit in schizophrenia and thought disorder. Brain. 2008;131:180–195. [PMC free article] [PubMed]
  • Nathan SS, Sinha SR, Gordon B, Lesser RP, Thakor NV. Determination of current-density distributions generated by electrical-stimulation of the human cerebral-cortex. Electroencephalogra. Clin. Neurophysiol. 1993;86:183–192. [PubMed]
  • Nicholson PW. Specific impedance of cerebral white matter. Exp. Neurology. 1965;13:386–401. [PubMed]
  • Nuttall GA, Bowersox MR, Douglass SB, McDonald J, Rasmussen LJ, Decker PA, Oliver WC, Rasmussen KG. Morbidity and mortality in the use of electroconvulsive therapy. J. ECT. 2004;20:237–241. [PubMed]
  • Oostendorp TF, Hengeveld YA, Wolters CH, Stinstra J, van Elswijk G, Stegeman DF. Modeling transcranial DC stimulation; Conf. Proc. IEEE Eng. Med. Biol. Soc; 2008. pp. 4226–4229. [PubMed]
  • Oppenheimer S. Cerebrogenic cardiac arrhythmias: cortical lateralization and clinical significance. Clin. Auton. Res. 2006;16:6–11. [PMC free article] [PubMed]
  • Oppenheimer SM, Gelb A, Girvin JP, Hachinski VC. Cardiovascular effects of human insular cortex stimulation. Neurology. 1992;42:1727–1732. [PubMed]
  • Ottosson JO. Experimental studies of memory impairment after electroconvulsive therapy. the role of the electrical stimulation and of the seizure studied by variation of stimulus intensity and modification by lidocaine of seizure discharge. Acta. Psychiatr. Scand. Suppl. 1960;35:103–131. [PubMed]
  • Parazzini M, Fiocchi S, Rossi E, Paglialonga A, Ravazzani P. Transcranial direct current stimulation: estimation of the electric field and of the current density in an anatomical human head model. IEEE Trans. Biomed. Eng. 2011;58:1773–1780. [PubMed]
  • Park JH, Hong SB, Kim DW, Suh M, Im CH. A novel array-type transcranial direct current stimulation (tDCS) system for accurate focusing on targeted brain areas. IEEE Trans. Magn. 2011;47:882–885.
  • Peterchev AV, Rosa MA, Deng ZD, Prudic J, Lisanby SH. Electroconvulsive therapy stimulus parameters: rethinking dosage. J. ECT. 2010;26:159–174. [PMC free article] [PubMed]
  • Plonsey R, Heppner DB. Considerations of quasi-stationarity in electrophysiological systems. Bull. Math. Biophys. 1967;29:657–664. [PubMed]
  • Pons JP, Segonne E, Boissonnat JD, Rineau L, Yvinec M, Keriven R. High-quality consistent meshing of multi-label datasets. Inf. Process. Med. Imaging. 2007;20:198–210. [PubMed]
  • Radman T, Ramos RL, Brumberg JC, Bikson M. Role of cotical cell type and morphology in subthreshold and suprathreshold uniform electric field sitmulation in vitro. Brain Stimulation. 2009;2:215–228. [PMC free article] [PubMed]
  • Rampersad S, Stegeman D, Oostendorp T. On handling the layered structure of the skull in transcranial direct current stimulation models; Conf. Proc. IEEE Eng. Med. Biol. Soc; 2011. pp. 1989–1992. [PubMed]
  • Ranta ME, Crocetti D, Clauss JA, Kraut MA, Mostofsky SH, Kaufmann WE. Manual MRI parcellation of the frontal lobe. Psychiatry Res. 2009;172:147–154. [PMC free article] [PubMed]
  • Rattay F. Modeling the excitation of fibers under surface electrodes. IEEE Trans. Biomed. Eng. 1988;35:199–202. [PubMed]
  • Rullmann M, Anwander A, Dannhauer M, Warfield SK, Duffy FH, Wolters CH. EEG source analysis of epileptiform activity using a 1 mm anisotropic hexahedra finite element head model. Neuroimage. 2009;44:399–410. [PMC free article] [PubMed]
  • Rush S, Driscoll DA. Current distribution in brain from surface electrodes. Anesth. Analg. 1968;47:717–723. [PubMed]
  • Rushton WAH. The effect upon the threshold for nervous excitation of the length of nerve exposed, and the angle between current and nerve. J. Physiol. 1927;9:357–377. [PubMed]
  • Sackeim HA. Convulsant and anticonvulsant properties of electroconvulsive therapy: towards a focal form of brain stimulation. Clin. Neurosci. Res. 2004;4:39–57.
  • Sackeim HA, Devanand DP, Prudic J. Stimulus intensity, seizure threshold, and seizure duration: impact on the efficacy and safety of electroconvulsive therapy. Psychiatr. Clin. North Am. 1991;14:803–843. [PubMed]
  • Sackeim HA, Mukherjee S. Neurophysiological variability in the effects of the ECT Stimulus. Convuls. Ther. 1986;2:267–276. [PubMed]
  • Sackeim HA, Prudic J, Devanand DP, Decina P, Kerr B, Malitz S. The impact of medication resistance and continuation pharmacotherapy on relapse following response to electroconvulsive therapy in major depression. J. Clin. Psychopharmacol. 1990;10:96–104. [PubMed]
  • Sackeim HA, Prudic J, Devanand DP, Nobler MS, Lisanby SH, Peyser S, Fitzsimons L, Moody BJ, Clark J. A prospective, randomized, double-blind comparison of bilateral and right unilateral electroconvulsive therapy at different stimulus intensities. General Psychiatry. 2000;57:425–434. [PubMed]
  • Sackeim HA, Prudic J, Fuller R, Keilp J, Lavori PW, Olfson M. The cognitive effects of electroconvulsive therapy in community settings. Neuropsychopharmacology. 2007;32:244–254. [PubMed]
  • Sackeim HA, Prudic J, Nobler MS, Fitzsimons L, Lisanby SH, Payne N, Berman RM, Brakerneier EL, Perera T, Devanand DP. Effects of pulse width and electrode placement on the efficacy and cognitive effects of electroconvulsive therapy. Brain Stimulation. 2008;1:71–83. [PMC free article] [PubMed]
  • Sadleir RJ, Argibay A. Modeling skull electrical properties. Ann. Biomed. Eng. 2007;35:1699–1712. [PMC free article] [PubMed]
  • Sadleir RJ, Vannorsdall TD, Schretlen DJ, Gordon B. Transcranial direct current stimulation (tDCS) in a realistic head model. Neuroimage. 2010;51:1310–1318. [PubMed]
  • Salvador R, Mekonnen A, Ruffini G, Miranda PC. Modeling the electric field induced in a high resolution realistic head model during transcranial current stimulation. Conf. Proc. IEEE Eng. Med. Biol. Soc. 2010;1:2073–2076. [PubMed]
  • Saypol JM, Roth BJ, Cohen LG, Hallett M. A theoretical comparison of electric and magnetic stimulation of the brain. Ann. Biomed. Eng. 1991;19:317–328. [PubMed]
  • Sekino M, Ueno S. Comparison of current distributions in electroconvulsive therapy and transcranial magnetic stimulation. J. Appl. Phys. 2002;91:8730–8732.
  • Sekino M, Ueno S. FEM-based determination of optimum current distribution in transcranial magnetic stimulation as an alternative to electroconvulsive therapy. IEEE Trans. Magn. 2004;40:2167–2169.
  • Shattuck DW, Leahy RM. BrainSuite: an automated cortical surface identification tool. Med. Image Anal. 2002;6:129–142. [PubMed]
  • Shimony JS, McKinstry RC, Akbudak E, Aronovitz JA, Snyder AZ, Lori NF, Cull TS, Conturo TE. Quantitative diffusion-tenser anisotropy brain MR imaging: normative human data and anatomic analysis. Radiology. 1999;212:770–784. [PubMed]
  • Smitt JW, Wegener CF. ON ELECTRIC CONVULSIVE THERAPY with particular regard to a parietal application of electrodes, controlled by intracerebral voltage measurements. Acta. Psychiatrica. Et. Neurologica. 1944;19:529–549.
  • Spellman T, Peterchev AV, Lisanby SH. Focal electrically administered seizure therapy: a novel form of ECT illustrates the roles of current directionality, polarity, and electrode configuration in seizure induction. Neuropsychopharmacology. 2009;34:2002–2010. [PubMed]
  • Squire LR. Memory functions as affected by electroconvulsive-therapy. Ann. New York Acad. Sci. 1986;462:307–314. [PubMed]
  • Squire LR, Slater PC, Miller PL. Retrograde-amnesia and bilateral electroconvulsive-therapy-long-term follow-up. General Psychiatry. 1981;38:89–95. [PubMed]
  • Stecker MM. Transcranial electric stimulation of motor pathways: a theoretical analysis. Comput. Biol. Med. 2005;35:133–155. [PubMed]
  • Suh HS, Kim SH, Lee WH, Kim TS. Realistic simulation of transcranial direct current stimulation via 3-D high-resolution finite element analysis: effect of tissue anisotropy; Conf. Proc. IEEE Eng. Med. Biol. Soc; 2009. pp. 638–641. [PubMed]
  • Suh HS, Lee WH, Cho YS, Kim JH, Kim TS. Reduced spatial focality of electrical field in tDCS with ring electrodes due to tissue anisotropy. Conf. Proc. IEEE Eng. Med. Biol. Soc. 2010;1:2053–2056. [PubMed]
  • Suihko V. Modelling the response of scalp sensory receptors to transcranial electrical stimulation. Med. Biol. Eng. Comput. 2002;40:395–401. [PubMed]
  • Szmurło R, Sawicki B, Starzynski J, Wincenciak S. A comparison of two models of electrodes for ECT simulations. IEEE Trans. Magn. 2006;42:1395–1398.
  • Thielscher A, Opitz A, Windhoff M. Impact of the gyral geometry on the electric field induced by transcranial magnetic stimulation. Neuroimage. 2011;54:234–243. [PubMed]
  • Tsopelas C, Stewart R, Savva GM, Brayne C, Ince P, Thomas A, Matthews FE. Neuropathological correlates of late-life depression in older people. Br. J. Psychiatry. 2011;198:109–114. [PubMed]
  • Tu Z, Narr KL, Dollar P, Dinov I, Thompson PM, Toga AW. Brain anatomical structure segmentation by hybrid discriminative/generative models. IEEE Trans. Med. Imaging. 2008;27:495–508. [PMC free article] [PubMed]
  • Tuch DS, Wedeen VJ, Dale AM, George JS, Belliveau JW. Conductivity mapping of biological tissue using diffusion MRI. Ann. N. Y. Acad. Sci. 1999;888:314–316. [PubMed]
  • Tuch DS, Wedeen VJ, Dale AM, George JS, Belliveau JW. Conductivity tensor mapping of the human brain using diffusion tensor MRI. Proc. Natl. Acad. Sci. USA. 2001;98:11697–11701. [PubMed]
  • Wagner T, Fregni F, Fecteau S, Grodzinsky A, Zahn M, Pascual-Leone A. Transcranial direct current stimulation: a computer-based human model study. Neuroimage. 2007;35:1113–1124. [PubMed]
  • Wang K, Zhu S, Mueller BA, Lim KO, Liu ZM, He B. A new method to derive white matter conductivity from diffusion tensor MRI. IEEE Trans. Biomed. Eng. 2008;55:2481–2486. [PMC free article] [PubMed]
  • Weaver L, Williams R, Rush S. Current density in bilateral and unilateral ECT. Biol. Psychiatry. 1976;11:303–312. [PubMed]
  • Weiss AP, DeWitt I, Goff D, Ditman T, Heckers S. Anterior and posterior hippocampal volumes in schizophrenia. Schizophrenia Res. 2005;73:103–112. [PubMed]
  • Welch CA, Drop LJ. Cardiovascular effects of ECT. Convulsive Ther. 1989;5:35–43. [PubMed]
  • Wells WM, 3rd, Viola P, Atsumi H, Nakajima S, Kikinis R. Multi-modal volume registration by maximization of mutual information. Med. Image Anal. 1996;1:35–51. [PubMed]
  • Wible CG, Shenton ME, Fischer IA, Allard JE, Kikinis R, Jolesz FA, Iosifescu DV, McCarley RW. Parcellation of the human prefrontal cortex using MRI. Psychiatry Res. 1997;76:29–40. [PubMed]
  • Wolters CH, Anwander A, Tricoche X, Weinstein D, Koch MA, MacLeod RS. Influence of tissue conductivity anisotropy on EEG/MEG field and return current computation in a realistic head model: a simulation and visualization study using high-resolution finite element modeling. Neuroimage. 2006;30:813–826. [PubMed]
  • Wu F, Tang Y, Xu K, Kong L, Sun W, Wang F, Kong D, Li Y, Liu Y. White matter abnormalities in medication-naive subjects with a single short-duration episode of major depressive disorder. Psychiatry Res. 2011;191:80–83. [PMC free article] [PubMed]