Dose-response in modulating brain function with transcranial direct current stimulation: From local to network levels.

Understanding the dose-response relationship is crucial in studying the effects of brain stimulation techniques, such as transcranial direct current stimulation (tDCS). The dose-response relationship refers to the relationship between the received stimulation dose and the resulting response, which can be described as a function of the dose at various levels, including single/multiple neurons, clusters, regions, or networks. Here, we are focused on the received stimulation dose obtained from computational head models and brain responses which are quantified by functional magnetic resonance imaging (fMRI) data. In this randomized, triple-blind, sham-controlled clinical trial, we recruited sixty participants with methamphetamine use disorders (MUDs) as a sample clinical population who were randomly assigned to receive either sham or active tDCS. Structural and functional MRI data, including high-resolution T1 and T2-weighted MRI, resting-state functional MRI, and a methamphetamine cue-reactivity task fMRI, were acquired before and after tDCS. Individual head models were generated using the T1 and T2-weighted MRI data to simulate electric fields. In a linear approach, we investigated the associations between electric fields (received dose) and changes in brain function (response) at four different levels: voxel level, regional level (using atlas-based parcellation), cluster level (identifying active clusters), and network level (task-based functional connectivity). At the voxel level, regional level, and cluster level, no FDR-corrected significant correlation was observed between changes in functional activity and electric fields. However, at the network level, a significant positive correlation was found between frontoparietal connectivity and the electric field at the frontopolar stimulation site (r = 0.42, p corrected = 0.02; medium effect size). Our proposed pipeline offers a methodological framework for analyzing tDCS effects by exploring dose-response relationships at different levels, enabling a direct link between electric field variability and the neural response to tDCS. The results indicate that network-based analysis provides valuable insights into the dependency of tDCS neuromodulatory effects on the individual's regional current dose. Integration of dose-response relationships can inform dose optimization, customization, or the extraction of predictive/treatment-response biomarkers in future brain stimulation studies.

Are we really targeting and stimulating DLPFC by placing transcranial electrical stimulation (tES) electrodes over F3/F4?

In many clinical trials involving transcranial electrical stimulation (tES), target electrodes are typically placed over DLPFC with the assumption that this will primarily stimulate the underlying brain region. However, our study aimed to evaluate the electric fields (EF) that are actually delivered and identify prefrontal regions that may be inadvertently targeted in DLPFC tES. Head models were generated from the Human Connectome Project database's T1 + T2-weighted MRIs of 80 healthy adults. Two common DLPFC montages were simulated; symmetric-F4/F3, and asymmetric-F4/Fp1. Averaged EF was extracted from (1) the center of the target electrode (F4), and (2) the top 1% of voxels showing the strongest EF in individualized EF maps. Interindividual variabilities were quantified with the standard deviation of EF peak location/value. Similar steps were repeated with 66 participants with methamphetamine use disorder (MUDs) as an independent clinical population. In healthy adults, the group-level location of EF peaks was situated in the medial-frontopolar, and the individualized EF peaks were positioned in a cube with a volume of 29 cm3 /46 cm3 (symmetric/asymmetric montages). EFs in the frontopolar area were significantly higher than EF "under" the target electrode in both symmetric (peak: 0.41 ± 0.06, F4:0.22 ± 0.04) and asymmetric (peak: 0.38 ± 0.04, F4:0.2 ± 0.04) montages (Heges'g > 0.7). Similar results with slight between-group differences were found in MUDs. We highlighted that in common DLPFC tES montages, in addition to interindividual/intergroup variability, the frontopolar received the highest EFs rather than DLPFC as the main target. We specifically recommended considering the potential involvement of the frontopolar area as a mechanism underlying the effectiveness of DLPFC tES protocols.

Transcranial direct current stimulation to modulate fMRI drug cue reactivity in methamphetamine users: A randomized clinical trial.

Transcranial direct current stimulation (tDCS) has been studied as a therapeutic option to alter maladaptive brain functions associated with chronic substance use. We present a randomized, triple-blind, sham-controlled, clinical trial to determine the neural substrates of tDCS effects on drug craving. Sixty participants with methamphetamine use disorder were assigned to two groups: active tDCS (5 x 7 cm2 , 2 mA, 20 min, anode/cathode over the F4/Fp1) and sham stimulation. Neuroimaging data of a methamphetamine cue reactivity task were collected immediately before and after stimulation. There was a significant reduction in self-reported craving after stimulation without any significant effect of time-by-group interaction. Our whole-brain analysis demonstrated that there was a global decrease in brain reactivity to cues following sham but not active tDCS. There were significant time-by-group interactions in five main clusters in middle and inferior frontal gyri, anterior insula, inferior parietal lobule, and precuneus with higher activations after active stimulation. There was a significant effect of stimulation type in the relationship between electrical current at the individual level and changes in task-modulated activation. Brain regions with the highest electric current in the prefrontal cortex showed a significant time-by-group interaction in task-modulated connectivity in the frontoparietal network. In this trial, there was no significant effect of the one session of active-F4/Fp1 tDCS on drug craving self-report compared to sham stimulation. However, activation and connectivity differences induced by active compared to sham stimulation suggested some potential mechanisms of tDCS to modulate neural response to drug cues.

Converging Evidence for Frontopolar Cortex as a Target for Neuromodulation in Addiction Treatment.

Noninvasive brain stimulation technologies such as transcranial electrical and magnetic stimulation (tES and TMS) are emerging neuromodulation therapies that are being used to target the neural substrates of substance use disorders. By the end of 2022, 205 trials of tES or TMS in the treatment of substance use disorders had been published, with heterogeneous results, and there is still no consensus on the optimal target brain region. Recent work may help clarify where and how to apply stimulation, owing to expanding databases of neuroimaging studies, new systematic reviews, and improved methods for causal brain mapping. Whereas most previous clinical trials targeted the dorsolateral prefrontal cortex, accumulating data highlight the frontopolar cortex as a promising therapeutic target for transcranial brain stimulation in substance use disorders. This approach is supported by converging multimodal evidence, including lesion-based maps, functional MRI-based maps, tES studies, TMS studies, and dose-response relationships. This review highlights the importance of targeting the frontopolar area and tailoring the treatment according to interindividual variations in brain state and trait and electric field distribution patterns. This converging evidence supports the potential for treatment optimization through context, target, dose, and timing dimensions to improve clinical outcomes of transcranial brain stimulation in people with substance use disorders in future clinical trials.

Cortical Morphology in Cannabis Use Disorder: Implications for Transcranial Direct Current Stimulation Treatment.

Introduction: Transcranial direct current stimulation (tDCS) has been studied as an adjunctive treatment option for substance use disorders (SUDs). Alterations in brain structure following SUD may change tDCS-induced electric field (EF) and subsequent responses; however, group-level differences between healthy controls (HC) and participants with SUDs in terms of EF and its association with cortical architecture have not yet been modeled quantitatively. This study provides a methodology for group-level analysis of computational head models to investigate the influence of cortical morphology metrics on EFs. Methods: Whole-brain surface-based morphology was conducted, and cortical thickness, volume, and surface area were compared between participants with cannabis use disorders (CUD) (n=20) and age-matched HC (n=22). Meanwhile, EFs were simulated for bilateral tDCS over the dorsolateral prefrontal cortex. The effects of structural alterations on EF distribution were investigated based on individualized computational head models. Results: Regarding EF, no significant difference was found within the prefrontal cortex; however, EFs were significantly different in left-postcentral and right-superior temporal gyrus (P<0.05) with higher levels of variance in CUD compared to HC [F(39, 43)=5.31, P<0.0001, C=0.95]. Significant differences were observed in cortical area (caudal anterior cingulate and rostral middle frontal), thickness (lateral orbitofrontal), and volume (paracentral and fusiform) between the two groups. Conclusion: Brain morphology and tDCS-induced EFs may be changed following CUD; however, differences between CUD and HCs in EFs do not always overlap with brain areas that show structural alterations. To sufficiently modulate stimulation targets, whether individuals with CUD need different stimulation doses based on tDCS target location should be checked.

Modeling of cerebellar transcranial electrical stimulation effects on hand tremor in Parkinson's disease.

Introduction: Parkinson's disease (PD) is a neurodegenerative disorder with different motor and neurocognitive symptoms. Tremor is a well-known symptom of this disease. Increasing evidence suggested that the cerebellum may substantially contribute to tremors as a clinical symptom of PD. However, the theoretical foundations behind these observations are not yet fully understood. Methods: In this study, a computational model is proposed to consider the role of the cerebellum and to show the effectiveness of cerebellar transcranial alternating current stimulation (tACS) on the rest tremor in participants with PD. The proposed model consists of the cortex, cerebellum, spinal circuit-muscular system (SC-MS), and basal ganglia blocks as the most critical parts of the brain, which are involved in generating rest tremors. The cortex, cerebellum, and SC-MS blocks were modeled using Van der Pol oscillators that interacted through synchronization procedures. Basal ganglia are considered as a regulator of the coupling weights defined between oscillators. In order to evaluate the global behavior of the model, we applied tACS on the cerebellum of fifteen PD patients for 15 min at each patient's peak frequency of their rest tremors. A tri-axial accelerometer recorded rest tremors before, during, and after the tACS. Results and Discussion: The simulation of the model provides a suggestion for the possible role of the cerebellum on rest tremors and how cerebellar tACS can affect these tremors. Results of human experiments also showed that the online and offline effects of cerebellar tACS could lead to the reduction of rest tremors significantly by about %76 and %68, respectively. Our findings suggest that the cerebellar tACS could serve as a reliable, therapeutic technique to suppress the PD tremor.

Optimizing Individual Targeting of Fronto-Amygdala Network with Transcranial Magnetic Stimulation (TMS): Biophysical, Physiological and Behavioral Variations in People with Methamphetamine Use Disorder.

Background: Previous studies in people with substance use disorders (SUDs) have implicated both the frontopolar cortex and amygdala in drug cue reactivity and craving, and amygdala-frontopolar coupling is considered a marker of early relapse risk. Accumulating data highlight that the frontopolar cortex can be considered a promising therapeutic target for transcranial magnetic stimulation (TMS) in SUDs. However, one-size-fits-all approaches to TMS targets resulted in substantial variation in both physiological and behavioral outcomes. Individualized TMS approaches to target cortico-subcortical circuits like amygdala-frontopolar have not yet been investigated in SUDs. Objective: Here, we (1) defined individualized TMS target location based on functional connectivity of the amygdala-frontopolar circuit while people were exposed to drug-related cues, (2) optimized coil orientation based on maximizing electric field (EF) perpendicular to the individualized target, and (3) harmonized EF strength in targeted brain regions across a population. Method: MRI data including structural, resting-state, and task-based fMRI data were collected from 60 participants with methamphetamine use disorders (MUDs). Craving scores based on a visual analog scale were collected immediately before and after the MRI session. We analyzed inter-subject variability in the location of TMS targets based on the maximum task-based connectivity between the left medial amygdala (with the highest functional activity among subcortical areas during drug cue exposure) and frontopolar cortex using psychophysiological interaction (PPI) analysis. Computational head models were generated for all participants and EF simulations were calculated for fixed vs. optimized coil location (Fp1/Fp2 vs. individualized maximal PPI location), orientation (AF7/AF8 vs. orientation optimization algorithm), and stimulation intensity (constant vs. adjusted intensity across the population). Results: Left medial amygdala with the highest (mean ± SD: 0.31±0.29) functional activity during drug cue exposure was selected as the subcortical seed region. Amygdala-to-whole brain PPI analysis showed a significant cluster in the prefrontal cortex (cluster size: 2462 voxels, cluster peak in MNI space: [25 39 35]) that confirms cortico-subcortical connections. The location of the voxel with the most positive amygdala-frontopolar PPI connectivity in each participant was considered as the individualized TMS target (mean ± SD of the MNI coordinates: [12.6 64.23 -0.8] ± [13.64 3.50 11.01]). Individual amygdala-frontopolar PPI connectivity in each participant showed a significant correlation with VAS scores after cue exposure (R=0.27, p=0.03). Averaged EF strength in a sphere with r = 5mm around the individualized target location was significantly higher in the optimized (mean ± SD: 0.99 ± 0.21) compared to the fixed approach (Fp1: 0.56 ± 0.22, Fp2: 0.78 ± 0.25) with large effect sizes (Fp1: p = 1.1e-13, Hedges'g = 1.5, Fp2: p = 1.7e-5, Hedges'g = 1.26). Adjustment factor to have identical 1 V/m EF strength in a 5mm sphere around the individualized targets ranged from 0.72 to 2.3 (mean ± SD: 1.07 ± 0.29). Conclusion: Our results show that optimizing coil orientation and stimulation intensity based on individualized TMS targets led to stronger electric fields in the targeted brain regions compared to a one-size-fits-all approach. These findings provide valuable insights for refining TMS therapy for SUDs by optimizing the modulation of cortico-subcortical circuits.

The cerebellum is causally involved in episodic memory under aging.

Episodic memory decline is a major signature of both normal and pathological aging. Many neural regions have been implicated in the processes subserving both episodic memory and typical aging decline. Here, we demonstrate that the cerebellum is causally involved episodic memory under aging. We show that a 12-day neurostimulation program delivered to the right cerebellum led to improvements in episodic memory performance under healthy aging that long outlast the stimulation period - healthy elderly individuals show episodic memory improvement both immediately after the intervention program and in a 4-month follow-up. These results demonstrate the causal relevance of the cerebellum in processes associated with long-term episodic memory, potentially highlighting its role in regulating and maintaining cognitive processing. Moreover, they point to the importance of non-pharmacological interventions that prevent or diminish cognitive decline in healthy aging.

Closing the loop between brain and electrical stimulation: towards precision neuromodulation treatments.

One of the most critical challenges in using noninvasive brain stimulation (NIBS) techniques for the treatment of psychiatric and neurologic disorders is inter- and intra-individual variability in response to NIBS. Response variations in previous findings suggest that the one-size-fits-all approach does not seem the most appropriate option for enhancing stimulation outcomes. While there is a growing body of evidence for the feasibility and effectiveness of individualized NIBS approaches, the optimal way to achieve this is yet to be determined. Transcranial electrical stimulation (tES) is one of the NIBS techniques showing promising results in modulating treatment outcomes in several psychiatric and neurologic disorders, but it faces the same challenge for individual optimization. With new computational and methodological advances, tES can be integrated with real-time functional magnetic resonance imaging (rtfMRI) to establish closed-loop tES-fMRI for individually optimized neuromodulation. Closed-loop tES-fMRI systems aim to optimize stimulation parameters based on minimizing differences between the model of the current brain state and the desired value to maximize the expected clinical outcome. The methodological space to optimize closed-loop tES fMRI for clinical applications includes (1) stimulation vs. data acquisition timing, (2) fMRI context (task-based or resting-state), (3) inherent brain oscillations, (4) dose-response function, (5) brain target trait and state and (6) optimization algorithm. Closed-loop tES-fMRI technology has several advantages over non-individualized or open-loop systems to reshape the future of neuromodulation with objective optimization in a clinically relevant context such as drug cue reactivity for substance use disorder considering both inter and intra-individual variations. Using multi-level brain and behavior measures as input and desired outcomes to individualize stimulation parameters provides a framework for designing personalized tES protocols in precision psychiatry.

DLPFC stimulation alters large-scale brain networks connectivity during a drug cue reactivity task: A tDCS-fMRI study.

Transcranial direct current stimulation (tDCS) is a promising intervention for reducing craving/consumption in individuals with substance use disorders. However, its exact mechanism of action has not yet been well explored. We aimed to examine the network-based effects of tDCS while people with methamphetamine use disorders (MUDs) were exposed to drug cues. In a randomized, double-blind sham-controlled trial with a crossover design, 15 participants with MUDs were recruited to receive 20 min of active/sham tDCS with an anode/cathode over F4/F3. MRI data, including structural and task-based functional MRI during a standard drug cue-reactivity task, were collected immediately before and after stimulation sessions. Craving scores were also recorded before and after MRI scans. Individualized head models were generated to determine brain regions with strong electric fields (EFs). Using atlas-based parcellation of head models, averaged EFs were extracted from the main nodes of three large-scale networks that showed abnormalities in MUDs; executive control (ECN), default mode (DMN), and ventral attention (VAN) networks. Main nodes with high EF intensity were used as seed regions for task-based functional connectivity (FC) [using generalized psychophysiological interaction (gPPI)] and activity [using a general linear model (GLM)] calculations. Subjective craving showed a significant reduction in immediate craving after active (-15.42 ± 5.42) compared to sham (-1 ± 2.63). In seed-to-whole brain results, the PFC node in ECN showed an enhanced PPI connectivity with precuneus and visual cortex; the cluster center in MNI (6, -84, -12); the PFC node in DMN showed a decreased PPI connectivity with contralateral parietal cortex;(-48, -60, 46). ROI-to-ROI results showed increased PPI connectivity within/between ECN-VAN while connectivity between ECN-DMN decreased. In line with connectivity, functional activity in the right PFC node in DMN decreased after tDCS while activity in PFC nodes of ECN/VAN increased. EF calculations in PFC nodes revealed that EF in DMN was outward, while the direction of EFs was inward in ECN/VAN. This study provides new insight into neural circuitry underlying MUDs that can be modulated by tDCS at the network level and specifically suggests that bilateral tDCS increases cortical excitability in ECN and VAN, while it has opposite effects on DMN that may be related to the direction of EFs.

How structural and functional MRI can inform dual-site tACS parameters: A case study in a clinical population and its pragmatic implications.

BACKGROUND: Abnormalities in frontoparietal network (FPN) were observed in many neuropsychiatric diseases including substance use disorders. A growing number of studies are using dual-site-tACS with frontoparietal synchronization to engage this network. However, a computational pathway to inform and optimize parameter space for frontoparietal synchronization is still lacking. In this case study, in a group of participants with methamphetamine use disorders, we proposed a computational pathway to extract optimal electrode montage while accounting for stimulation intensity using structural and functional MRI. METHODS: Sixty methamphetamine users completed an fMRI drug cue-reactivity task. Four main steps were taken to define electrode montage and adjust stimulation intensity using 4x1 high-definition (HD) electrodes for a dual-site-tACS; (1) Frontal seed was defined based on the maximum electric fields (EF) predicted by simulation of HD montage over DLPFC (F3/F4 in EEG 10-10), (2) frontal seed-to-whole brain context-dependent correlation was calculated to determine connected regions to frontal seeds, (3) center of connected cluster in parietal cortex was selected as a location for placing the second set of HD electrodes to shape the informed montage, (4) individualized head models were used to determine optimal stimulation intensity considering underlying brain structure. The informed montage was compared to montages with large electrodes and classic frontoparietal HD montages (F3-P3/F4-P4) in terms of tACS-induced EF and ROI-to-ROI task-based/resting-state connectivity. RESULTS: Compared to the large electrodes, HD frontoparietal montages allow for a finer control of the spatial peak fields in the main nodes of the FPN at the cost of lower maximum EF (large-pad/HD: max EF[V/m] = 0.37/0.11, number of cortical sub-regions that EF exceeds 50% of the max = 77/13). For defining stimulation targets based on EF patterns, using group-level head models compared to a single standard head model results in comparable but significantly different seed locations (6.43 mm Euclidean distance between the locations of the frontal maximum EF in standard-space). As expected, significant task-based/resting-state connections were only found between frontal-parietal locations in the informed montage. Cue-induced craving score was correlated with frontoparietal connectivity only in the informed montage (r = -0.24). Stimulation intensity in the informed montage, and not in the classic HD montage, needs 40% reduction in the parietal site to reduce the disparity in EF between stimulation sites. CONCLUSION: This study provides some empirical insights to montage and dose selection in dual-site-tACS using individual brain structures and functions and proposes a computational pathway to use head models and functional MRI to define (1) optimum electrode montage for targeting FPN in a context of interest (drug-cue-reactivity) and (2) proper transcranial stimulation intensity.

Cue-induced craving and negative emotion disrupt response inhibition in methamphetamine use disorder: Behavioral and fMRI results from a mixed Go/No-Go task.

BACKGROUND: Drug-related cue-reactivity, dysfunctional negative emotion processing, and response-disinhibition constitute three core aspects of methamphetamine use disorder (MUD). These phenomena have been studied independently, but the neuroscientific literature on their interaction in addictive disorders remains scant. METHODS: 62 individuals with MUD were scanned when responding to the geometric Go or No-Go cues superimposed over blank, neutral, negative-emotional and drug-related background images. Neural correlates of drug and negative-emotional cue-reactivity, response-inhibition and their interactions were estimated, and methamphetamine cue-reactivity was compared between individuals with MUD and 23 healthy controls. Relationships between behavioral characteristics and observed activations were investigated. RESULTS: Individuals with MUD had longer reaction times and more errors in drug and negative-emotional compared to blank blocks, and more omission errors in drug compared to neutral blocks. They showed higher drug cue-reactivity than controls across prefrontal, fusiform, and visual regions (Z > 3.1, p-corrected<0.05). Response-inhibition was associated with precuneal, inferior parietal, anterior cingulate, temporal, and inferior frontal activations (Z > 3.1, p-corrected<0.05). Response-inhibition in drug cue blocks coincided with higher activations in the visual cortex and lower activations in the paracentral lobule and superior and inferior frontal gyri, while inhibition during negative-emotional blocks led to higher superior parietal, fusiform, and lateral occipital activations (Z > 3.1, p-corrected<0.05). CONCLUSION: Drug cue-reactivity may impair response inhibition partly through activating dis-inhibitory regions, while temporal and parietal activations associated with response-inhibition in negative blocks suggest compensatory activity. Results suggest that drug and negative-emotional cue-reactivity influence response-inhibition, and the study of these interactions may aid mechanistic understanding of methamphetamine use disorder.

Online closed-loop real-time tES-fMRI for brain modulation: A technical report.

Recent studies suggest that transcranial electrical stimulation (tES) can be performed during functional magnetic resonance imaging (fMRI). The novel approach of using concurrent tES-fMRI to modulate and measure targeted brain activity/connectivity may provide unique insights into the causal interactions between the brain neural responses and psychiatric/neurologic signs and symptoms, and importantly, guide the development of new treatments. However, tES stimulation parameters to optimally influence the underlying brain activity may vary with respect to phase difference, frequency, intensity, and electrode's montage among individuals. Here, we propose a protocol for closed-loop tES-fMRI to optimize the frequency and phase difference of alternating current stimulation (tACS) for two nodes (frontal and parietal regions) in individual participants. We carefully considered the challenges in an online optimization of tES parameters with concurrent fMRI, specifically in its safety, artifact in fMRI image quality, online evaluation of the tES effect, and parameter optimization method, and we designed the protocol to run an effective study to enhance frontoparietal connectivity and working memory performance with the optimized tACS using closed-loop tES-fMRI. We provide technical details of the protocol, including electrode types, electrolytes, electrode montages, concurrent tES-fMRI hardware, online fMRI processing pipelines, and the optimization algorithm. We confirmed the implementation of this protocol worked successfully with a pilot experiment.

Group and individual level variations between symmetric and asymmetric DLPFC montages for tDCS over large scale brain network nodes.

Two challenges to optimizing transcranial direct current stimulation (tDCS) are selecting between, often similar, electrode montages and accounting for inter-individual differences in response. These two factors are related by how tDCS montage determines current flow through the brain considered across or within individuals. MRI-based computational head models (CHMs) predict how brain anatomy determines electric field (EF) patterns for a given tDCS montage. Because conventional tDCS produces diffuse brain current flow, stimulation outcomes may be understood as modulation of global networks. Therefore, we developed a network-led, rather than region-led, approach. We specifically considered two common "frontal" tDCS montages that nominally target the dorsolateral prefrontal cortex; asymmetric "unilateral" (anode/cathode: F4/Fp1) and symmetric "bilateral" (F4/F3) electrode montages. CHMs of 66 participants were constructed. We showed that cathode location significantly affects EFs in the limbic network. Furthermore, using a finer parcellation of large-scale networks, we found significant differences in some of the main nodes within a network, even if there is no difference at the network level. This study generally demonstrates a methodology for considering the components of large-scale networks in CHMs instead of targeting a single region and specifically provides insight into how symmetric vs asymmetric frontal tDCS may differentially modulate networks across a population.

fMRI and transcranial electrical stimulation (tES): A systematic review of parameter space and outcomes.

The combination of non-invasive brain stimulation interventions with human brain mapping methods have supported research beyond correlational associations between brain activity and behavior. Functional MRI (fMRI) partnered with transcranial electrical stimulation (tES) methods, i.e., transcranial direct current (tDCS), transcranial alternating current (tACS), and transcranial random noise (tRNS) stimulation, explore the neuromodulatory effects of tES in the targeted brain regions and their interconnected networks and provide opportunities for individualized interventions. Advances in the field of tES-fMRI can be hampered by the methodological variability between studies that confounds comparability/replicability. In order to explore variability in the tES-fMRI methodological parameter space (MPS), we conducted a systematic review of 222 tES-fMRI experiments (181 tDCS, 39 tACS and 2 tRNS) published before February 1, 2019, and suggested a framework to systematically report main elements of MPS across studies. Publications dedicated to tRNS-fMRI were not considered in this systematic review. We have organized main findings in terms of fMRI modulation by tES. tES modulates activation and connectivity beyond the stimulated areas particularly with prefrontal stimulation. There were no two studies with the same MPS to replicate findings. We discuss how to harmonize the MPS to promote replication in future studies.