Functional neuroimaging of the baboon during concurrent image-guided transcranial magnetic stimulation. Neuroimage. 2011;57:1393-1401. [PMID: 21664276 DOI: 10.1016/j.neuroimage.2011.05.065]

What can neuroimaging of neuromodulation reveal about the basis of circuit therapies for psychiatry?

Neuromodulation is increasingly becoming a therapeutic option for treatment resistant psychiatric disorders. These non-invasive and invasive therapies are still being refined but are clinically effective and, in some cases, provide sustained symptom reduction. Neuromodulation relies on changing activity within a specific brain region or circuit, but the precise mechanisms of action of these therapies, is unclear. Here we review work in both humans and animals that has provided insight into how therapies such as deep brain and transcranial magnetic stimulation alter neural activity across the brain. We focus on studies that have combined neuromodulation with neuroimaging such as PET and MRI as these measures provide detailed information about the distributed networks that are modulated and thus insight into both the mechanisms of action of neuromodulation but also potentially the basis of psychiatric disorders. Further we highlight work in nonhuman primates that has revealed how neuromodulation changes neural activity at different scales from single neuron activity to functional connectivity, providing key insight into how neuromodulation influences the brain. Ultimately, these studies highlight the value of combining neuromodulation with neuroimaging to reveal the mechanisms through which these treatments influence the brain, knowledge vital for refining targeted neuromodulation therapies for psychiatric disorders.

Transcranial magnetic stimulation in non-human primates: A systematic review

Transcranial magnetic stimulation (TMS) is widely employed as a tool to investigate and treat brain diseases. However, little is known about the direct effects of TMS on the brain. Non-human primates (NHPs) are a valuable translational model to investigate how TMS affects brain circuits given their neurophysiological similarity with humans and their capacity to perform complex tasks that approach human behavior. This systematic review aimed to identify studies using TMS in NHPs as well as to assess their methodological quality through a modified reference checklist. The results show high heterogeneity and superficiality in the studies regarding the report of the TMS parameters, which have not improved over the years. This checklist can be used for future TMS studies with NHPs to ensure transparency and critical appraisal. The use of the checklist would improve methodological soundness and interpretation of the studies, facilitating the translation of the findings to humans. The review also discusses how advancements in the field can elucidate the effects of TMS in the brain.

Mapping Brain Motor Functions Using Transcranial Magnetic Stimulation with a Volume Conductor Model and Electrophysiological Experiments

Transcranial magnetic stimulation (TMS) activates brain cells in a noninvasive manner and can be used for mapping brain motor functions. However, the complexity of the brain anatomy prevents the determination of the exact location of the stimulated sites, resulting in the limitation of the spatial resolution of multiple targets. The aim of this study is to map two neighboring muscles in cortical motor areas accurately and quickly. Multiple stimuli were applied to the subject using a TMS stimulator to measure the motor-evoked potentials (MEPs) in the corresponding muscles. For each stimulation condition (coil location and angle), the induced electric field (EF) in the brain was computed using a volume conductor model for an individualized head model of the subject constructed from magnetic resonance images. A post-processing method was implemented to determine a TMS hotspot using EF corresponding to multiple stimuli, considering the amplitude of the measured MEPs. The dependence of the computationally estimated hotspot distribution on two target muscles was evaluated (n = 11). The center of gravity of the first dorsal interosseous cortical representation was lateral to the abductor digiti minimi by a minimum of 2 mm. The localizations were consistent with the putative sites obtained from previous EF-based studies and fMRI studies. The simultaneous cortical mapping of two finger muscles was achieved with only several stimuli, which is one or two orders of magnitude smaller than that in previous studies. Our proposal would be useful in the preoperative mapping of motor or speech areas to plan brain surgery interventions.

A Baboon Brain Atlas for Magnetic Resonance Imaging and Positron Emission Tomography Image Analysis

The olive baboon (Papio anubis) is phylogenetically proximal to humans. Investigation into the baboon brain has shed light on the function and organization of the human brain, as well as on the mechanistic insights of neurological disorders such as Alzheimer's and Parkinson's. Non-invasive brain imaging, including positron emission tomography (PET) and magnetic resonance imaging (MRI), are the primary outcome measures frequently used in baboon studies. PET functional imaging has long been used to study cerebral metabolic processes, though it lacks clear and reliable anatomical information. In contrast, MRI provides a clear definition of soft tissue with high resolution and contrast to distinguish brain pathology and anatomy, but lacks specific markers of neuroreceptors and/or neurometabolites. There is a need to create a brain atlas that combines the anatomical and functional/neurochemical data independently available from MRI and PET. For this purpose, a three-dimensional atlas of the olive baboon brain was developed to enable multimodal imaging analysis. The atlas was created on a population-representative template encompassing 89 baboon brains. The atlas defines 24 brain regions, including the thalamus, cerebral cortex, putamen, corpus callosum, and insula. The atlas was evaluated with four MRI images and 20 PET images employing the radiotracers for [11C]benzamide, [11C]metergoline, [18F]FAHA, and [11C]rolipram, with and without structural aids like [18F]flurodeoxyglycose images. The atlas-based analysis pipeline includes automated segmentation, registration, quantification of region volume, the volume of distribution, and standardized uptake value. Results showed that, in comparison to PET analysis utilizing the "gold standard" manual quantification by neuroscientists, the performance of the atlas-based analysis was at >80 and >70% agreement for MRI and PET, respectively. The atlas can serve as a foundation for further refinement, and incorporation into a high-throughput workflow of baboon PET and MRI data. The new atlas is freely available on the Figshare online repository (https://doi.org/10.6084/m9.figshare.16663339), and the template images are available from neuroImaging tools & resources collaboratory (NITRC) (https://www.nitrc.org/projects/haiko89/).

The baboon in epilepsy research: Revelations and challenges

The baboon offers a natural model for genetic generalized epilepsy with photosensitivity. In this review, we will summarize some of the more important clinical, neuroimaging, and elctrophysiological findings form recent work performed at the Southwest National Primate Research Center (SNPRC, Texas Biomedical Research Institute, San Antonio, Texas), which houses the world's largest captive baboon pedigree. Due to the phylogenetic proximity of the baboon to humans, many of the findings are readily translatable, but there may be some important differences, such as the mutlifocality of the ictal and interictal epileptic discharges (IEDs) on intracranial electroencephalography (EEG) and greater parieto-occipital connectivity of baboon brain networks compared to juvenile myoclonic epilepsy in humans. Furthermore, there is still limited knowledge of the natural history of the epilepsy, which could be transformative for research into epileptogenesis in genetic generalized epilepsy (GGE) and sudden unexpected death in epilepsy (SUDEP).

Altered Topological Organization in the Sensorimotor Network After Application of Different Frequency rTMS

The application of repetitive transcranial magnetic stimulation (rTMS) over the primary motor cortex (M1) could influence the intrinsic brain activity in the sensorimotor network (SMN). However, how rTMS modulates the topological organization of the SMN remains unclear. In this study, we employed resting-state fMRI to investigate the topological alterations in the functional SMN after application of different frequency rTMS over the left M1. To accomplish this, we collected MRI data from 45 healthy participants who were randomly divided into three groups based on rTMS frequency (HF, high-frequency 3 Hz; LF, low-frequency 1 Hz; and SHAM). Individual large-scale functional SMN was constructed by correlating the mean time series among 29 regions of interest (ROI) in the SMN and was fed into graph-based network analyses at multiple levels of global organization and nodal centrality. Our results showed that compared with the network metrics before rTMS stimulation, the left paracentral lobule (PCL) exhibited reduced nodal degree and betweenness centrality in the LF group after rTMS, while the right supplementary motor area (SMA) exhibited reduced nodal betweenness centrality in the HF group after rTMS. Moreover, rTMS-related alterations in nodal metrics might have been attributable to the changes in connectivity patterns and local activity of the affected nodes. These findings reflected the potential of using rTMS over M1 as an effective intervention to promote motor function rehabilitation.

Mapping Structure-Function Relationships in the Brain

Mapping the structural and functional connectivity of the brain is a major focus of systems neuroscience research and will help to identify causally important changes in neural circuitry responsible for behavioral dysfunction. Several methods for examining brain activity in humans have been extended to rodent and monkey models in which molecular and genetic manipulations exist for linking to human disease. In this review, which is part of a special issue focused on bridging brain connectivity information across species and spatiotemporal scales, we address mapping brain activity and neural connectivity in rodents using optogenetics in conjunction with either functional magnetic resonance imaging or optical intrinsic signal imaging. We chose to focus on these techniques because they are capable of reporting spontaneous or evoked hemodynamic activity most closely linked to human neuroimaging studies. We discuss the capabilities and limitations of blood-based imaging methods, usage of optogenetic techniques to map neural systems in rodent models, and other powerful mapping techniques for examining neural connectivity over different spatial and temporal scales. We also discuss implementing strategies for mapping brain connectivity in humans with both basic and clinical applications, and conclude with how cross-species mapping studies can be utilized to influence preclinical imaging studies and clinical practices alike.

A high-resolution computational localization method for transcranial magnetic stimulation mapping

BACKGROUND

Transcranial magnetic stimulation (TMS) is used for the mapping of brain motor functions. The complexity of the brain deters determining the exact localization of the stimulation site using simplified methods (e.g., the region below the center of the TMS coil) or conventional computational approaches.

OBJECTIVE

This study aimed to present a high-precision localization method for a specific motor area by synthesizing computed non-uniform current distributions in the brain for multiple sessions of TMS.

METHODS

Peritumoral mapping by TMS was conducted on patients who had intra-axial brain neoplasms located within or close to the motor speech area. The electric field induced by TMS was computed using realistic head models constructed from magnetic resonance images of patients. A post-processing method was implemented to determine a TMS hotspot by combining the computed electric fields for the coil orientations and positions that delivered high motor-evoked potentials during peritumoral mapping. The method was compared to the stimulation site localized via intraoperative direct brain stimulation and navigated TMS.

RESULTS

Four main results were obtained: 1) the dependence of the computed hotspot area on the number of peritumoral measurements was evaluated; 2) the estimated localization of the hand motor area in eight non-affected hemispheres was in good agreement with the position of a so-called "hand-knob"; 3) the estimated hotspot areas were not sensitive to variations in tissue conductivity; and 4) the hand motor areas estimated by this proposal and direct electric stimulation (DES) were in good agreement in the ipsilateral hemisphere of four glioma patients.

CONCLUSION(S)

The TMS localization method was validated by well-known positions of the "hand-knob" in brains for the non-affected hemisphere, and by a hotspot localized via DES during awake craniotomy for the tumor-containing hemisphere.

Anaesthesia, not number of sessions, influences the magnitude and duration of an aHF-rTMS in dogs

BACKGROUND

Currently, the rat has been a useful animal model in brain stimulation research. Nevertheless, extrapolating results from rodent repetitive Transcranial Magnetic Stimulation (rTMS) research to humans contains several hurdles. This suggests the desperate need for a large animal model in translational rTMS research. The dog would be a valid choice, not only due to the fact that humans and dogs share a neurophysiological background, but a similar neuropathological background as well.

HYPOTHESIS

In order to evaluate the feasibility of the canine rTMS animal model, this study aimed to evaluate the neurophysiological response in dogs on a, clinically used, accelerated high frequency (aHF) rTMS protocol. This aHF-rTMS (20 Hz) protocol was performed under anaesthesia or sedation and either 20 sessions or 5 sessions were given to each dog.

METHODS

21 healthy dogs were randomly subjected to one of the four aHF-rTMS protocols (1 sham and 3 active protocols). For each dog, the perfusion indices (PI), of a [99mTc]HMPAO scan at 4 time points, for the left frontal cortex (stimulation target) were calculated for each protocol.

RESULTS

Concerning sham stimulation, the average PI remained at the baseline level. The main result was the presence of a direct transitory increase in rCBF at the stimulation site, both under anaesthesia and sedation. Nevertheless the measured increase in rCBF was higher but shorter duration under sedation. The magnitude of this increase was not influenced by number of sessions. No changes in rCBF were found in remote brain regions.

CONCLUSION

This study shows that, despite the influence of anaesthesia and sedation, comparable and clinically relevant effects on the rCBF can be obtained in dogs. Since less methodological hurdles have to be overcome and comparable results can be obtained, it would be acceptable to put the dog forward as an alternative translational rTMS animal model.

Repetitive Transcranial Magnetic Stimulation Educes Frequency-Specific Causal Relationships in the Motor Network

BACKGROUND

Repetitive transcranial magnetic stimulation (rTMS) has the potential to treat brain disorders by modulating the activity of disease-specific brain networks, yet the rTMS frequencies used are delivered in a binary fashion - excitatory (>1 Hz) and inhibitory (≤1 Hz).

OBJECTIVE

To assess the effective connectivity of the motor network at different rTMS stimulation rates during positron-emission tomography (PET) and confirm that not all excitatory rTMS frequencies act on the motor network in the same manner.

METHODS

We delivered image-guided, supra-threshold rTMS at 3 Hz, 5 Hz, 10 Hz, 15 Hz and rest (in separate randomized sessions) to the primary motor cortex (M1) of the lightly anesthetized baboon during PET imaging. Each rTMS/PET session was analyzed using normalized cerebral blood flow (CBF) measurements. Path analysis - using structural equation modeling (SEM) - was employed to determine the effective connectivity of the motor network at all rTMS frequencies. Once determined, the final model of the motor network was used to assess any differences in effective connectivity at each rTMS frequency.

RESULTS

The exploratory SEM produced a very well fitting final network model (χ(2) = 18.04, df = 21, RMSEA = 0.000, p = 0.647, TLI = 1.12) using seven nodes of the motor network. 5 Hz rTMS produced the strongest path coefficients in four of the seven connections, suggesting that this frequency is the optimal rTMS frequency for stimulation the motor network (as a whole); however, the premotor cerebellum connection was optimally stimulated at 10 Hz rTMS and the supplementary motor area caudate connection was optimally driven at 15 Hz rTMS.

CONCLUSION(S)

We have demonstrated that 1) 5 Hz rTMS revealed the strongest path coefficients (i.e. causal influence) on the nodes of the motor network, 2) stimulation at "excitatory" rTMS frequencies did not produce increased CBF in all nodes of the motor network, 3) specific rTMS frequencies may be used to target specific none-to-node interactions in the stimulated brain network, and 4) more research needs to be performed to determine the optimum frequency for each brain circuit and/or node.

Resting-state functional connectivity in the baboon model of genetic generalized epilepsy

OBJECTIVE

The baboon provides a natural model of genetic generalized epilepsy (GGE). This study compares the intrinsic connectivity networks of epileptic and healthy control baboons using resting-state functional magnetic resonance imaging (rs-fMRI) and data-driven functional connectivity mapping.

METHODS

Twenty baboons, matched for gender, age, and weight, were classified into two groups (10 epileptic [EPI], 10 control [CTL]) on the basis of scalp electroencephalography (EEG) findings. Each animal underwent one MRI session that acquired one 5-min resting state fMRI scan and one anatomic MRI scan-used for registration and spatial normalization. Using independent component analysis, we identified 14 unique components/networks, which were then used to characterize each group's functional connectivity maps of each brain network.

RESULTS

The epileptic group demonstrated network-specific differences in functional connectivity when compared to the control animals. The sensitivity and specificity of the two groups' functional connectivity maps differed significantly in the visual, motor, amygdala, insular, and default mode networks. Significant increases were found in the occipital gyri of the epileptic group's functional connectivity map for the default mode, cingulate, intraparietal, motor, visual, amygdala, and thalamic regions.

SIGNIFICANCE

This is the first study using resting-state fMRI to demonstrate intrinsic functional connectivity differences between epileptic and control nonhuman primates. These results are consistent with seed-based GGE studies in humans; however, our use of a data-driven approach expands the scope of functional connectivity mapping to include brain regions/networks comprising the whole brain.

Computational and experimental analysis of TMS-induced electric field vectors critical to neuronal activation

OBJECTIVE

Transcranial magnetic stimulation (TMS) represents a powerful technique to noninvasively modulate cortical neurophysiology in the brain. However, the relationship between the magnetic fields created by TMS coils and neuronal activation in the cortex is still not well-understood, making predictable cortical activation by TMS difficult to achieve. Our goal in this study was to investigate the relationship between induced electric fields and cortical activation measured by blood flow response. Particularly, we sought to discover the E-field characteristics that lead to cortical activation.

APPROACH

Subject-specific finite element models (FEMs) of the head and brain were constructed for each of six subjects using magnetic resonance image scans. Positron emission tomography (PET) measured each subject's cortical response to image-guided robotically-positioned TMS to the primary motor cortex. FEM models that employed the given coil position, orientation, and stimulus intensity in experimental applications of TMS were used to calculate the electric field (E-field) vectors within a region of interest for each subject. TMS-induced E-fields were analyzed to better understand what vector components led to regional cerebral blood flow (CBF) responses recorded by PET.

MAIN RESULTS

This study found that decomposing the E-field into orthogonal vector components based on the cortical surface geometry (and hence, cortical neuron directions) led to significant differences between the regions of cortex that were active and nonactive. Specifically, active regions had significantly higher E-field components in the normal inward direction (i.e., parallel to pyramidal neurons in the dendrite-to-axon orientation) and in the tangential direction (i.e., parallel to interneurons) at high gradient. In contrast, nonactive regions had higher E-field vectors in the outward normal direction suggesting inhibitory responses.

SIGNIFICANCE

These results provide critical new understanding of the factors by which TMS induces cortical activation necessary for predictive and repeatable use of this noninvasive stimulation modality.

Electric Field Model of Transcranial Electric Stimulation in Nonhuman Primates: Correspondence to Individual Motor Threshold

OBJECTIVE

To develop a pipeline for realistic head models of nonhuman primates (NHPs) for simulations of noninvasive brain stimulation, and use these models together with empirical threshold measurements to demonstrate that the models capture individual anatomical variability.

METHODS

Based on structural MRI data, we created models of the electric field (E-field) induced by right unilateral (RUL) electroconvulsive therapy (ECT) in four rhesus macaques. Individual motor threshold (MT) was measured with transcranial electric stimulation (TES) administered through the RUL electrodes in the same subjects.

RESULTS

The interindividual anatomical differences resulted in 57% variation in median E-field strength in the brain at fixed stimulus current amplitude. Individualization of the stimulus current by MT reduced the E-field variation in the target motor area by 27%. There was significant correlation between the measured MT and the ratio of simulated electrode current and E-field strength (r(2) = 0.95, p = 0.026). Exploratory analysis revealed significant correlations of this ratio with anatomical parameters including of the superior electrode-to-cortex distance, vertex-to-cortex distance, and brain volume (r(2) > 0.96, p < 0.02). The neural activation threshold was estimated to be 0.45 ±0.07 V/cm for 0.2-ms stimulus pulse width.

CONCLUSION

These results suggest that our individual-specific NHP E-field models appropriately capture individual anatomical variability relevant to the dosing of TES/ECT. These findings are exploratory due to the small number of subjects.

SIGNIFICANCE

This study can contribute insight in NHP studies of ECT and other brain stimulation interventions, help link the results to clinical studies, and ultimately lead to more rational brain stimulation dosing paradigms.

Value of evoked potential in study of functional bowel disorders

In recent years, the development of neurogastroenterology and the application of neurophysiologic examinations have gradually revealed the association of gastrointestinal activity with cortical activity through the efferent and afferent pathways. The state of nerve conduction pathway between the brain and gastrointestinal tract is closely related with specific functions of the anus and rectum. Scholars have put forward the "brain gut axis" and "brain gut interaction" theories to explain the bidirectional interaction between the gastrointestinal tract and central nervous system. Evoked potentials have an important role in the diagnosis and research of electrophysiological changes in various parts of the neural system, which provides practical information for the study of the brain-gut pathway, promotes the diagnosis and understanding of diseases related to the brain-gut axis abnormalities, and provides the basis for developing new treatment methods. In this paper, we summarize the roles of various evoked potential techniques in the study of functional bowel disorders.

Concurrent TMS to the primary motor cortex augments slow motor learning

Transcranial magnetic stimulation (TMS) has shown promise as a treatment tool, with one FDA approved use. While TMS alone is able to up- (or down-) regulate a targeted neural system, we argue that TMS applied as an adjuvant is more effective for repetitive physical, behavioral and cognitive therapies, that is, therapies which are designed to alter the network properties of neural systems through Hebbian learning. We tested this hypothesis in the context of a slow motor learning paradigm. Healthy right-handed individuals were assigned to receive 5 Hz TMS (TMS group) or sham TMS (sham group) to the right primary motor cortex (M1) as they performed daily motor practice of a digit sequence task with their non-dominant hand for 4 weeks. Resting cerebral blood flow (CBF) was measured by H2(15)O PET at baseline and after 4 weeks of practice. Sequence performance was measured daily as the number of correct sequences performed, and modeled using a hyperbolic function. Sequence performance increased significantly at 4 weeks relative to baseline in both groups. The TMS group had a significant additional improvement in performance, specifically, in the rate of skill acquisition. In both groups, an improvement in sequence timing and transfer of skills to non-trained motor domains was also found. Compared to the sham group, the TMS group demonstrated increases in resting CBF specifically in regions known to mediate skill learning namely, the M1, cingulate cortex, putamen, hippocampus, and cerebellum. These results indicate that TMS applied concomitantly augments behavioral effects of motor practice, with corresponding neural plasticity in motor sequence learning network. These findings are the first demonstration of the behavioral and neural enhancing effects of TMS on slow motor practice and have direct application in neurorehabilitation where TMS could be applied in conjunction with physical therapy.

Repetitive transcranial magnetic stimulation elicits rate-dependent brain network responses in non-human primates

BACKGROUND

Transcranial magnetic stimulation (TMS) has the potential to treat brain disorders by tonically modulating firing patterns in disease-specific neural circuits. The selection of treatment parameters for clinical repetitive transcranial magnetic stimulation (rTMS) trials has not been rule based, likely contributing to the variability of observed outcomes.

OBJECTIVE

To utilize our newly developed baboon (Papio hamadryas anubis) model of rTMS during position-emission tomography (PET) to quantify the brain's rate-response functions in the motor system during rTMS.

METHODS

We delivered image-guided, suprathreshold rTMS at 3 Hz, 5 Hz, 10 Hz, 15 Hz and rest (in separate randomized sessions) to the primary motor cortex (M1) of the lightly anesthetized baboon during PET imaging; we also administered a (reversible) paralytic to eliminate any somatosensory feedback due to rTMS-induced muscle contractions. Each rTMS/PET session was analyzed using normalized cerebral blood flow (CBF) measurements; statistical parametric images and the resulting areas of significance underwent post-hoc analysis to determine any rate-specific rTMS effects throughout the motor network.

RESULTS

The motor system's rate-response curves were unimodal and system wide--with all nodes in the network showing highly similar rate response functions--and an optimal network stimulation frequency of 5 Hz.

CONCLUSION(S)

These findings suggest that non-invasive brain stimulation may be more efficiently delivered at (system-specific) optimal frequencies throughout the targeted network and that functional imaging in non-human primates is a promising strategy for identifying the optimal treatment parameters for TMS clinical trials in specific brain regions and/or networks.

Repetitive magnetic stimulation induces functional and structural plasticity of excitatory postsynapses in mouse organotypic hippocampal slice cultures

Repetitive transcranial magnetic stimulation (rTMS) is a noninvasive brain stimulation technique that can alter cortical excitability in human subjects for hours beyond the stimulation period. It thus has potential as a therapeutic tool in neuropsychiatric disorders associated with alterations in cortical excitability. However, rTMS-induced neural plasticity remains insufficiently understood at the cellular level. To learn more about the effects of repetitive magnetic stimulation (rMS), we established an in vitro model of rMS using mouse organotypic entorhino-hippocampal slice cultures. We assessed the outcome of a high-frequency (10 Hz) rMS protocol on functional and structural properties of excitatory synapses in mature hippocampal CA1 pyramidal neurons. Whole-cell patch-clamp recordings, immunohistochemistry, and time-lapse imaging techniques revealed that rMS induces a long-lasting increase in glutamatergic synaptic strength, which is accompanied by structural remodeling of dendritic spines. The effects of rMS on spine size were predominantly seen in small spines, suggesting differential effects of rMS on subpopulations of spines. Furthermore, our data indicate that rMS-induced postsynaptic changes depend on the NMDA receptor-mediated accumulation of GluA1-containing AMPA receptors. These results provide first experimental evidence that rMS induces coordinated functional and structural plasticity of excitatory postsynapses, which is consistent with a long-term potentiation of synaptic transmission.

Where does transcranial magnetic stimulation (TMS) stimulate? Modelling of induced field maps for some common cortical and cerebellar targets

Computational models have been be used to estimate the electric and magnetic fields induced by transcranial magnetic stimulation (TMS) and can provide valuable insights into the location and spatial distribution of TMS stimulation. However, there has been little translation of these findings into practical TMS research. This study uses the International 10-20 EEG electrode placement system to position a standard figure-of-eight TMS coil over 13 commonly adopted targets. Using a finite element method and an anatomically detailed and realistic head model, this study provides the first pictorial and numerical atlas of TMS-induced electric fields for a range of coil positions. The results highlight the importance of subject-specific gyral folding patterns and of local thickness of subarachnoid cerebrospinal fluid (CSF). Our modelling shows that high electric fields occur primarily on the peaks of those gyri which have only a thin layer of CSF above them. These findings have important implications for inter-individual generalizability of the TMS-induced electric field. We propose that, in order to determine with accuracy the site of stimulation for an individual subject, it is necessary to solve the electric field distribution using subject-specific anatomy obtained from a high-resolution imaging modality such as MRI.