Keywords
transcranial direct current stimulation (tDCS), transcranial magnetic stimulation (TMS), non-invasive brain stimulation (NIBS), major depressive disorder (MDD), Head models, computational modelling ,magnetic resonance imaging (MRI)
transcranial direct current stimulation (tDCS), transcranial magnetic stimulation (TMS), non-invasive brain stimulation (NIBS), major depressive disorder (MDD), Head models, computational modelling ,magnetic resonance imaging (MRI)
Non-invasive brain stimulation (NIBS) techniques such as transcranial direct current stimulation (tDCS) and transcranial magnetic stimulation (TMS) have been used to investigate the relationship between activity in different cortical regions and cognitive processes. A key advantage of NIBS is that it allows direct manipulation of neural excitability. Therefore when used carefully, it allows causal interpretation of how specific brain regions might be involved in mental phenomena such as perception, working memory, attention, decision-making or emotional regulation. In addition, NIBS has been used as a clinical tool to obtain symptom reduction in several neurological and psychiatric conditions1–4. Importantly, the same stimulation protocol can result in different neural effects across individuals because the distribution of stimulation-induced electric fields (E-fields) in the brain is strongly contingent on anatomical variability5,6. This can manifest in strong variability in the effects of NIBS on cognitive performance in healthy individuals (e.g., working memory7,8 and on treatment outcomes in patients (e.g., depression9–13).
Given the heterogeneity in the efficacy of NIBS protocols in modulating behavior and clinical symptoms, there has been a move towards computational modelling of the spatial distribution of E-fields in the brain to understand how stimulation parameters such as electrode placement, rotation or type affects current flow in the neural tissue. After an initial phase during which simplistic spherical head models were used14, the focus has shifted to very detailed, realistic head models created from individual structural magnetic resonance (MR) images using freely available tools (e.g. SimNIBS15 or ROAST16). Creating individual head models remains challenging, however, due to (1) the requirement of high-quality structural MR images for each study participant and (2) the need to manually improve the automatic segmentation of the MR images into the different tissue types (i.e., bone, cerebrospinal fluid, white and grey matter, air, etc). As a consequence, individual head models are rarely used in practice. Instead, researchers usually rely on E-fields induced in a reference individual which is generalized to all participants. This approach, of course, neglects the importance of individual brain anatomy which can have strong influences on E-field distribution5.
To circumvent this issue, the New York (NY) Head model17 was created from the ICBM152 (v6 template18,19 and v2009b template20), which is a non-linear average of 152 individual MRIs and is extended to the full head and neck region by fusing it with an average of an additional 26 brains. Therefore, this head model represents an unbiased population average and should be a better representation of the individual participants than any randomly picked reference individual. However, calculation of the electric field on a single head model does not allow to quantify variation across individuals that can be substantial both in cortical regions directly below the electrodes and those that are farther away from the stimulation sites. Also, given that the NY Head was created using healthy individuals, possible systematic differences between patient groups may not be noticed.
To improve this situation, we present 38 head models created from MR images of 19 healthy participants and 19 patients with major depressive disorder (MDD). These head models were manually checked and edited to improve their accuracy in giving a faithful representation of cortical regions. Our head models can aid future research in at least three ways: First, calculating electric fields across a large sample allows to get a sense of the variability of the induced E-field due to anatomical differences. Second, our head models also enable comparison of the neural effects on NIBS protocols in healthy vs. depressed brains. Given the different protocols used for left dorsolateral prefrontal cortex (lDLPFC) targeting with NIBS21,22, computational modelling using these head models will enable the development of protocols for more selective lDLPFC stimulation in this disorder. Third, EEG source localization also relies on head models to calculate the spatial location of possible current sources in the brain. Usually, boundary element models (BEM) are used which have the limitation of less anatomical detail. Therefore, finite element models (FEM) derived from high resolution MR images have become more widespread because they are able to incorporate more tissue types, increasing the precision of EEG source localization. Our head models can be used for source localization using open-source software23. Thus, our head models can help researchers to optimize NIBS protocols and EEG source localization methods, and to test them on a larger sample (including both healthy and patient data).
High-resolution head models were created from T1-weighted anatomical images that were collected in a separate functional magnetic resonance imaging (fMRI) study24. The data was obtained from the OpenfMRI database (accession number:ds000171). Structural scans of 19 healthy adult participants with no history of depression or other psychiatric disorders (11 females; mean ± SD age: 28.79 ± 10.86; range: 18–59) and 19 individuals diagnosed with MDD and experiencing a depressive episode at the time of the scanning (11 females; mean ± SD age: 33.52 ± 13.35; range: 18–56) were used. Data of one control participant (’sub-control20’) was excluded due to technical problems with head model creation. Patients did not suffer from current or past manic episodes, current comorbid anxiety disorders or current alcohol dependence or abuse. At the time of data collection, all patients were unmedicated, but 6 received antidepressant pharmacotherapy in the past. For full details regarding demographic data, we refer to the supporting information of the original paper24.
As the very first step, we inspected scans of all participants, and manually removed signals corresponding to the MRI marker placed on the forehead of each subject using FreeSurfer 5.3.025. Automated tissue segmentation was performed in SPM1226 for skin, skull, eyeballs and CSF, and in FreeSurfer for gray and white matter. Subsequently, segmented images of each participant were visually inspected and manually corrected with FreeSurfer (done by investigator G.Cs., verified by O.P.). Manual corrections were primarily restricted to the skull-CSF boundary, but in some cases also involved the skin-skull interface. In addition, during manual corrections we verified that the segmentation of the cortical gray matter corresponded to the anatomical scans except for medial temporal lobe structures (i.e., the parahippocampal gyrus and the hippocampus proper). Moreover, the resulting masks were not corrected for inconsistencies relating to subcortical nuclei and thus, these head models are not suitable for estimating stimulation-related E-fields in structures such as the thalamus, basal ganglia, amygdala, or brainstem nuclei. Additionally, because the original dataset was de-faced and did not include the neck/shoulder region, our head models do not include these regions. This limits their usability regarding the simulation of tDCS montages with extracephalic return electrodes. Nevertheless, they can be used for all tES protocols using scalp electrodes and most TMS protocols.
Head models were created with an extended, pre-release version of SimNIBS 2.115, a freely available software package for simulating the effects of NIBS techniques. The final head mesh of each participant consisted of a total number of approximately 3,200,000 tetrahedral elements, assigned to six tissue types (Figure 1). For comparability with other datasets, we also report measurements for head size: The distance between the nasion and inion (mean ± SD: 19.2 cm ±1.01 cm; range: 16.7 cm – 21.3 cm) and the distance between the right and left pre-auricular points (mean ± SD: 14.8 cm ±0.65 cm; range: 13.6 cm – 16.3 cm). The total volume of the brain was 1.22 dm3 ±0.11 dm3 (range: 1.02 dm3 – 1.49 dm3). Individual measurements are found at our data repository27
Tissue conductivities were set as follows: 0.465 S/m (skin), 0.01 S/m (skull), 0.5 S/m (eyeballs), 1.654 S/m (cerebrospinal fluid), 0.275 S/m (gray matter), 0.126 S/m (white matter). Although our head models do not account for white matter anisotropy, this property has been shown to primarily influence current density in deeper structures, leaving superficial gray matter relatively unaffected28. The accuracy of tissue segmentation and the good correspondence between anatomical scans and the resulting head models for 8 individuals are shown in Figure 2.
Except for two manual steps (removal of MRI markers from the forehead, manual correction of tissue segmentations), the process of head model creation was automated using a pre-release version of SimNIBS 2.1 that employed FreeSurfer 5.3.0 for brain segmentation (as described in 29) and SPM12 for segmentation of the remaining tissues (similarly to 30). This new SimNIBS pipeline provides more accurate tissue segmentation relative to other protocols. The scripts used for automated simulation of tDCS effects for all head models (as shown in Figure 3) are available for download at our data repository27.
Please note the large degree of variability in E-field magnitude (both protocols) and in the lateralization of effects (bipolar protocol).
For validating the reliability of our head models, we compared the effects of two tDCS protocols targeting the lDLPFC (one conventional bipolar montage and one multi-electrode 4×1 setup) against the effects observed in the NY Head17. The mesh for the NY Head (abaqus format; https://www.parralab.org/nyhead/) has been reformatted to be compatible with SimNIBS and is also available for download in our data repository27. For each head model, tDCS electrodes of appropriate size (bipolar montage: 5 × 5 cm, circular connectors (diameter: 0.5 cm) at the middle of the electrode pads; 4×1 montage: diameter of 1.2 cm) and thickness (1 mm for all electrodes + 2.5 mm sponge pocket/gel layer for the bipolar and 4×1 montages, respectively) were placed at scalp locations corresponding to electrode positions of the International 10/10 system (bipolar montage: anode - F3, cathode - F4; 4×1 montage: anode: F3, cathodes: C3, FT7, Fp1, Fz). Stimulation intensity for the anode was set to 2 mA, with equal distribution of return currents for the 4 cathodes (0.5 mA for each) in the 4×1 protocol. Results of the simulations were visualized using Gmsh31 (Figure 3).
In our previous study6, we reported stronger stimulation-induced E-fields in the lDLPFC for the bipolar montage used by Brunoni and colleagues (2013)9 than for the 4×1 protocol, albeit the bipolar montage was also associated with reduced focality (i.e., more intensive stimulation of other cortical areas). Therefore, we extracted three measures for both tDCS protocols for our 38 head models: stimulation strength (the norm of the electric field vector, ’normfield’) in the lDLPFC was assessed by extracting individual mean (calculated across all nodes in this region) and maximum (peak) E-field values, whereas spatial focality of the stimulation was analyzed by calculating the focality-index (FI), with the target region as reference. FI was quantified as the proportion of highest-intensity nodes (nodes within the upper 1% percentile of all E-field values) in the lDLPFC relative to the whole cortex. Mean and peak E-field values corresponding to both tDCS montages were extracted by reconstructing the two-dimensional cortical surface (more precisely, the middle of the cortical sheet) of each individual along with the corresponding E-field cortical map in FreeSurfer, and applying automated atlas-based parcellation of the frontal lobe to delineate the lDLPFC region in each brain32.
As a result, we show that (1) both protocols induce strong E-fields in the DLPFC (with symmetrical effects for the bipolar montage and unilateral E-field distribution for the 4×1 protocol), (2) effects are similar for our head models and for the NY Head, (3) all E-field measures (peak and mean strength, FI) show great degree of variability, and (4) montage-specific effects are consistent with previous results reported in the literature regarding both the spatial distribution and the magnitude of E-fields33–37 (Figure 3, Figure 4).
Both peak (left panel) and mean (middle panel) E-field magnitudes in the left dorsolateral prefrontal cortex (lDLPFC) are stronger for the bipolar montage, whereas the 4×1 protocol yields more focal stimulation of the target region (right panel). Group means (red: healthy individuals, green: major depressive disorder patients) do not differ substantially, but large degree of inter-individual variability can be observed in both groups. The triangles show data for the New York Head. Horizontal lines within boxes represent median values, whereas lower and upper box hinges correspond to the first and third quartiles (25th and 75th percentiles). Lengths of upper/lower whiskers extend to the largest/smallest values that do not exceed 1.5* the inter-quartile range; dots represent individual data.
For detailed discussion of the results presented in Figure 4 we refer the reader to our accompanying paper6.
Open Science Framework (OSF): Dataset 1. Head models of healthy and depressed adults for simulating the effects of non-invasive brain stimulation, http://doi.org/10.17605/OSF.IO/EXBD527
License: CC0 1.0 Universal
At our data repository27, the following data are available for download for all subjects (healthy adults: ’sub-control’, patients: ’sub-mdd’): T1-weighted anatomical scans registered to FreeSurfer conform space (’nii.gz’ files), the corresponding segmentation masks for the 6 tissue types (’nii.gz’ files), the final head models (’msh’ files), the files containing electrode coordinates (International 10/10 system) for all participants (’txt’ files can be used for performing new script-based simulations, whereas the ’geo’ files can be used to plot electrode positions directly onto the mesh files in Gmsh31), and files organized in 2 folders (’fs_*.tar.gz’ and ’m2m_*.tar.gz’) that enable creating simulation outputs in average FreeSurfer space (’fsaverage’). We also included a README file with a detailed description of the data and scripts.
Our head models are compatible with SimNIBS 2.0 (http://simnibs.de/) for simulating the effects of tDCS and TMS protocols. This software package has an easy-to-use graphical user interface (GUI) for setting all stimulation parameters (scalp location, intensity, etc.) for both NIBS techniques. By using our custom-written script27, it is also possible to run tDCS simulations for any given montage for all participants at once. Although the current version of SimNIBS does not allow exporting the resulting E-field maps into FreeSurfer, this option will be available in an upcoming version. This will enable registering data to an average surface (’fsaverage’) and creating group averaged data, as we have shown previously6. In addition, using the upcoming new version of SimNIBS, researchers will have the opportunity to extract E-field components that are either radial (normal) or tangential relative to the cortical surface, and have been associated with different cellular effects38. At our data repository27 we also provide the manually corrected segmentations for the different tissue types for those who would like to create high-quality meshes of their own using open-source software such as iso2mesh39. Finally, our meshes can be used for improving the anatomic precision of EEG source localization, using open-source tools23.
This work was supported by the Northern Norway Regional Health Authority (grant no. PFP1237-15) for G.Cs. and M.M. and by the Lundbeck foundation (grant no. R118-A11308) and the Novonordisk foundation (grant no. NNF14OC0011413) for O.P. and A.T. The publication charges for this article have been funded by a grant from the publication fund of UiT The Arctic University of Norway.
The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Views | Downloads | |
---|---|---|
F1000Research | - | - |
PubMed Central
Data from PMC are received and updated monthly.
|
- | - |
Is the rationale for creating the dataset(s) clearly described?
Yes
Are the protocols appropriate and is the work technically sound?
Yes
Are sufficient details of methods and materials provided to allow replication by others?
Partly
Are the datasets clearly presented in a useable and accessible format?
Yes
References
1. Huang Y, Liu A, Lafon B, Friedman D, et al.: Measurements and models of electric fields in the in vivo human brain during transcranial electric stimulation. eLife. 2017; 6. Publisher Full TextCompeting Interests: No competing interests were disclosed.
Is the rationale for creating the dataset(s) clearly described?
Yes
Are the protocols appropriate and is the work technically sound?
Yes
Are sufficient details of methods and materials provided to allow replication by others?
Yes
Are the datasets clearly presented in a useable and accessible format?
Yes
Competing Interests: No competing interests were disclosed.
Alongside their report, reviewers assign a status to the article:
Invited Reviewers | ||
---|---|---|
1 | 2 | |
Version 2 (revision) 15 Nov 18 |
read | |
Version 1 06 Jun 18 |
read | read |
Provide sufficient details of any financial or non-financial competing interests to enable users to assess whether your comments might lead a reasonable person to question your impartiality. Consider the following examples, but note that this is not an exhaustive list:
Sign up for content alerts and receive a weekly or monthly email with all newly published articles
Already registered? Sign in
The email address should be the one you originally registered with F1000.
You registered with F1000 via Google, so we cannot reset your password.
To sign in, please click here.
If you still need help with your Google account password, please click here.
You registered with F1000 via Facebook, so we cannot reset your password.
To sign in, please click here.
If you still need help with your Facebook account password, please click here.
If your email address is registered with us, we will email you instructions to reset your password.
If you think you should have received this email but it has not arrived, please check your spam filters and/or contact for further assistance.
Comments on this article Comments (0)