1H MR spectroscopy of the motor cortex immediately following transcranial direct current stimulation at 7 Tesla

Abstract

Transcranial direct current stimulation (tDCS) is a form of non-invasive brain stimulation that may modulate cortical excitability, metabolite concentration, and human behaviour. The supplementary motor area (SMA) has been largely ignored as a potential target for tDCS neurorehabilitation but is an important region in motor compensation after brain injury with strong efferent connections to the primary motor cortex (M1). The objective of this work was to measure tissue metabolite changes in the human motor cortex immediately following tDCS. We hypothesized that bihemispheric tDCS would change levels of metabolites involved in neuromodulation including N-acetylaspartate (NAA), glutamate (Glu), and creatine (tCr). In this single-blind, randomized, cross-over study, fifteen healthy adults aged 21–60 participated in two 7T MRI sessions, to identify changes in metabolite concentrations by magnetic resonance spectroscopy. Immediately after 20 minutes of tDCS, there were no significant changes in metabolite levels or metabolite ratios comparing tDCS to sham. However there was a trend toward increased NAA/tCr concentration (p = 0.08) in M1 under the stimulating cathode. There was a strong, positive correlation between the change in the absolute concentration of NAA and the change in the absolute concentration of tCr (p<0.001) suggesting an effect of tDCS. Both NAA and creatine are important markers of neurometabolism. Our findings provide novel insight into the modulation of neural metabolites in the motor cortex immediately following application of bihemispheric tDCS.

About the Authors:

Kayla Ryan

Roles Conceptualization, Data curation, Formal analysis, Methodology, Writing – original draft

Affiliations Department of Medical Biophysics, The University of Western Ontario, London, Ontario, Canada, Centre for Functional and Metabolic Mapping, Robarts Research Institute, The University of Western Ontario, London, Ontario, Canada

ORCID logo http://orcid.org/0000-0001-9572-7410

Krzysztof Wawrzyn

Roles Data curation, Formal analysis, Resources

Affiliation: Department of Physics and Astronomy, University of Western Ontario, London, ON, Canada

Joseph S. Gati

Roles Conceptualization, Methodology, Resources, Supervision

Affiliation: Centre for Functional and Metabolic Mapping, Robarts Research Institute, The University of Western Ontario, London, Ontario, Canada

Blaine A. Chronik

Roles Resources, Software, Writing – review & editing

Affiliation: Department of Physics and Astronomy, University of Western Ontario, London, ON, Canada

Dickson Wong

Roles Formal analysis, Writing – review & editing

Affiliations Department of Medical Biophysics, The University of Western Ontario, London, Ontario, Canada, Centre for Functional and Metabolic Mapping, Robarts Research Institute, The University of Western Ontario, London, Ontario, Canada

ORCID logo http://orcid.org/0000-0002-4584-7251

Neil Duggal

Roles Conceptualization, Supervision, Writing – review & editing

Affiliations Department of Medical Biophysics, The University of Western Ontario, London, Ontario, Canada, Centre for Functional and Metabolic Mapping, Robarts Research Institute, The University of Western Ontario, London, Ontario, Canada, Department of Clinical Neurological Sciences, University Hospital, London Health Sciences Centre, London, Ontario, Canada

Robert Bartha

Roles Conceptualization, Project administration, Resources, Supervision, Writing – review & editing

* E-mail: [email protected]

Affiliations Department of Medical Biophysics, The University of Western Ontario, London, Ontario, Canada, Centre for Functional and Metabolic Mapping, Robarts Research Institute, The University of Western Ontario, London, Ontario, Canada

ORCID logo http://orcid.org/0000-0002-7118-9096

Introduction

Transcranial direct current stimulation (tDCS) is a form of non-invasive brain stimulation that has shown promise in modulating cortical excitability and behaviour in humans [1–4]. However, many facets of its use remain controversial. For example, the optimization of stimulation parameters (current level, duration, electrode montage, etc.), characterization of individual response variability, and the physiological mechanism of action are still under active investigation.

Several potential mechanisms of action have been proposed based on pharmacological, behavioural, and imaging studies [1–3, 5, 6]. A review by Stagg and colleagues summarized the evidence linking the effects of tDCS to a physiological mechanism. They concluded that tDCS modulates cortical excitability by altering membrane potential during the stimulation period, while the after effects are driven by synaptic modification and plasticity[6]. At current levels of ~1 mA, tDCS is thought to induce alterations of the membrane potential, with anodal tDCS making it more likely, and cathodal tDCS making it less likely, for an action potential to fire [7, 8]. However, the level of current may also impact the neuronal response to stimulation. At 2 mA, cathodal stimulation has been shown to have an excitatory influence on membrane potential. This has been thought to occur due to an increased release of calcium (Ca2+) at the higher current [9]. Furthermore, higher current penetrates deeper into cortical tissues, potentially causing dendritic depolarization at a sufficient level to modulate excitation of adjacent neuronal structures [9]. Magnetic resonance spectroscopy (MRS) provides a means to investigate the effects of tDCS on cellular metabolism and synaptic transmission as it can be used to non-invasively quantify cerebral metabolites in vivo, including glutamate (Glu) and gamma-aminobutyric acid (GABA). Previous MRS studies have shown changes in excitatory and inhibitory neurotransmitter levels, minutes after tDCS, with current levels ranging from 1–2 mA [6, 7, 10–14]. Other studies have suggested that creatine may have an important role in bioenergetics and neuromodulation [15–17]. For example, Rae et al. found an increase in adenosine-triphosphate (ATP) synthesis, with a decrease in the concentration of phosphocreatine in the left temporo-frontal region following anodal tDCS to the left dorsolateral prefrontal cortex [17]. In another study, 2mA of anodal tDCS to the right parietal cortex caused an increase in both Glx and total N-aceytl-aspartate (NAA + NAAG) relative to sham, measured from the parietal cortex [10], while a study by Stagg and colleagues found that 1 mA of cathodal stimulation to left M1 decreased Glx under the electrode[7]. Other studies have found no effect. For example, Kim et. al. found no changes after 1.5 mA of cathodal tDCS to left M1 in any metabolite measured under the stimulating electrode [18]. Similarly, using 1mA of current in an M1-M1 bihemispheric montage, Tremblay et. al. found no significant changes in any metabolite in left M1 [19]. These conflicting results are difficult to interpret, and leads to uncertainty with regards to the implementation of an optimum stimulation paradigm.

The application of tDCS to improve motor performance and recovery in neurological disorders requires optimization of stimulation parameters. Bihemispheric tDCS can enhance both behaviour and physiological responses in healthy and neurologically injured individuals [20–22]. The supplementary motor area (SMA) has proven to be an important area of the brain during the execution of bimanual hand movements [23], and plays a compensatory role during the recovery of both stroke and spinal cord injury [24–26]. With its strongest efferent projections to M1 and the corticospinal tract, SMA is a unique target for tDCS [27]. Support for this notion comes from a recent study that showed enhanced motor performance by targeting the left SMA with 0.4 mA of anodal tDCS for 90 min over three days [28]. By targeting both SMA and M1 with 2mA of tDCS, it may be possible to induce additive effects on M1 excitability via interhemispheric connections, which are thought to be more focal than those associated with M1-supraortibal stimulation [19, 20].

The purpose of the current study was to demonstrate the feasibility of concurrent tDCS and 7T MRI, and to determine whether targeting both SMA and M1 using a bihemispheric tDCS montage would produce immediate changes in metabolite concentrations in M1 measured using ultra high-field (7T) MRS. To our knowledge, this is the first study to examine the metabolic changes after bihemispheric tDCS, delivered in the MR environment, at an ultra-high magnetic field strength. Based on previous studies, we hypothesized that 2 mA bihemispheric tDCS would enhance synaptic and metabolic activity [10, 11, 15, 17]. As such, metabolites involved in neurometabolism such as NAA, glutamate, and creatine would be altered by stimulation.

Methods

Participants and study design

15 healthy adults aged 21–60 years (mean ± standard deviation: 28 ± 10, 9 female), with no reported history of mental or neurological illness, participated in two sessions on a 7 Tesla (Siemens, Erlangen) head–only MRI scanner. All participants had 1H MRS in this single blind, sham controlled, cross-over design. Participants were randomized to receive tDCS stimulation or sham stimulation on their initial visit, and the contrary on their second visit, at least 7 days apart. Informed written consent was obtained for all procedures according to the Declaration of Helinski (World Medical Association, 2008) and the study was approved by the Western University Health Sciences Research Ethics Board.

Transcranial direct current stimulation

Using an MR-compatible DC-STIMULATOR (NeuroConn, Germany), 2 mA of current was applied to bihemispheric motor areas in the MRI scanner, for a total of 20 minutes. Electrodes were 3x3 cm2, providing a total current density of 0.22 mA/cm2 and a total charge with respect to time of 0.27 C/cm2. For use inside the scanner, electrodes were fit with 5 kOhm resistors placed next to the electrode to minimize the possibility of eddy currents induced in the leads during MRS acquisition. Electrodes were positioned on each participant outside the magnet using the electroencephalography (EEG) 10–10 system, which has been shown to be a reliable localization tool [29]. The cathode was placed on the left primary motor cortex (C3), anode on the right supplementary motor area (FC2). For stimulation, current was ramped up over 10 seconds to reach 2 mA and held constant for 20 minutes, followed by a 10 s ramp down period. During sham stimulation, current was ramped up over 10 s and then immediately turned off. This method has been shown to be a reliable blinding method, as participants are unable to distinguish the difference between sham and true tDCS using this paradigm. We used this measurement as a baseline comparison [3, 30].

Temperature monitoring

To ensure the safety of the participants during tDCS in the MRI, temperature was monitored on all subjects throughout the duration of the scan (approximately an hour and 15 minutes). Specifically, four T1C 1.7 mm diameter fibre optic temperature sensors (Neoptix, Quebec, Canada) were located under both electrode pads and the nearest cable chokes. Temperature was monitored in real time with a calibrated Reflex signal conditioner (Neoptix, Quebec, Canada) and a custom data collection program written in LabVIEW 2010 (National Instruments).

Magnetic resonance image acquisition and analysis

A 7 Tesla Siemens (Erlangen, Germany), head-only MRI (Magnetom) was used to acquire spectroscopy and imaging data. Data were acquired using an 8 channel transmit and 32 channel receive coil array. T1-weighted MP2RAGE anatomical images (TE/TR = 2.83/6000 ms and 750 μm isotropic resolution) were acquired and used for voxel positioning. These images were also used to estimate white-matter (WM), gray-matter (GM) and cerebrospinal fluid (CSF) fractions for partial volume correction when determining metabolite concentration. The MRS acquisition began immediately following the completion of the stimulation to capture alterations in metabolite concentration due to tDCS. Water suppressed (64 averages) and unsuppressed (8 averages) 1H MR spectra were acquired from a single voxel (1.6×2.0×1.8 cm3) located in the left primary motor cortex (under the cathode) (Fig 1) using the semi Localization by Adiabatic Selective Refocusing (semi-LASER) pulse sequence [31]: TE/TR = 60/7500 ms, voxel size = 1.6×2.0×1.8 cm3,total MRS acquisition time was approximately 10 minutes. A localized B0 and B1 shim were applied prior to data acquisition. The B0 shim was optimized using a two-echo gradient recalled echo (GRE) shimming technique [32] and the B1 field was optimized such that the phases of the transmit channels added constructively within the MRS voxel. Spectra were lineshape corrected using combined QUALITY deconvolution and eddy current correction (QUECC) with 400 QUALITY points [33]. Simulated prior knowledge metabolite lineshapes were fitted to post-processed spectra using the fitMAN software developed in-house (Fig 2) [34]. Metabolite concentrations were examined as ratios normalized to creatine and also as absolute concentrations using unsuppressed water as an internal reference standard as previously described [35]. Measurement of tissue partial volume with the voxel was made using the MP2RAGE images in FMRIB Software Library (FSL) [36] to obtain the fraction of WM, GM and CSF within the voxel. In addition, relaxation rates of the metabolites were incorporated into the quantification to correct for T1 and T2 relaxation induced signal loss [37–40].

[Figure omitted. See PDF.]

Fig 1. Voxel positioning.

Typical MP2RAGE anatomical images used for voxel placement were brain extracted using FSL. The voxel shown in green (2x2x2 cm3) was placed over the left primary motor cortex (under the cathode).

https://doi.org/10.1371/journal.pone.0198053.g001

[Figure omitted. See PDF.]

Fig 2. Spectrum of left motor cortex.

7 Tesla semi-LASER 1H MRS (TE = 60 ms) of the left primary motor cortex. The spectrum (grey) is overlaid on the fitted result (black) with the residual shown above (black). Select metabolite peaks are identified.

https://doi.org/10.1371/journal.pone.0198053.g002

Metabolites measured with a group coefficient of variation of less than <30% in the sham condition were included in statistical analyses. To identify differences in metabolite levels and metabolite ratios between sham and tDCS conditions, repeated measures MANOVA was performed in SPSS (IBM SPSS Statistics Version 25). The main factors were the type of stimulation (two levels: sham and tDCS) and the metabolite (six levels: N-acetyl aspartate (NAA), myo-inositol (mI), creatine (Cr), choline (Cho), glutamate (Glu), glutathione (GSH) or metabolite ratio (five levels: NAA/Cr, mI/Cr, Cho/Cr, Glu/Cr, GSH/Cr) measured. In addition, differences between sham and tDCS conditions were compared separately for each metabolite and metabolite ratio using paired t-tests. Correlation analysis was completed using the Pearson Product Correlation Coefficient and the coefficient of determination (R2).

Results

Temperature monitoring

The tDCS was safely and successfully applied in the 7T MRI environment in all subjects. The average temperature change in all four probes was 4.3 ± 0.2°C throughout the duration of the experiment. This temperature increase was largely due to warming of the bore and from the participant’s natural body heating. Once equilibrium was established, small fluctuations on the order of 1°C were observed during periods when radio frequency (RF) fields were turned on.

Metabolite ratio changes

Spectral quality measures including signal to noise ration and linewidth are summarized in Table 1 for all participants. There were no age or gender related effects. When examining the metabolite ratios the repeated measures MANOVA indicated a trend for the effect of stimulation (F(1,14) = 3.52, p = 0.08). In addition, there was a significant main effect of metabolite (F(12,3) = 343.35, p<0.001). There was no main interaction effect (F(12,3) = 1.25, p = 0.33) (Table 2). Post-hoc t-tests (uncorrected for multiple comparisons) were performed to confirm alterations in metabolite ratios. Our results showed a trend toward a 4% increase in the NAA/tCr ratio between sham and tDCS conditions (p = 0.08, Cohen’s d = 0.52) and no significant changes in any other metabolite ratios (Fig 3, Table 2).

[Figure omitted. See PDF.]

Fig 3. Metabolic ratios–sham vs stim.

tDCS increases NAA/tCr ratio when measured immediately after a 20 minute stimulation period (* p = 0.08). No difference was observed in Glu/Cr. Error bars indicate SEM.

https://doi.org/10.1371/journal.pone.0198053.g003

[Figure omitted. See PDF.]

Table 1. Spectral quality and quantification: Characterization of spectral quality and voxel tissue composition.

Signal to noise ratio represents the intensity of the NAACH3 peak divided by the standard deviation of the baseline noise after Fourier transformation of the initial 0.3 seconds. The linewidth represents the full width at half maximum (FWHM) of the unsuppressed water signal. The water area represents the area of the unsuppressed water spectrum. The voxel tissue partial volume is provided for gray matter (GM), white matter (WM) and cerebral spinal fluid (CSF). Data is presented as mean ± standard error of the mean. Repeated measured t-tests were conducted on all spectral parameters; no significant changes between sham and stimulation were observed.

https://doi.org/10.1371/journal.pone.0198053.t001

[Figure omitted. See PDF.]

Table 2. Analysis of metabolic ratios: Metabolite ratios relative to total creatine.

The p-values were calculated in post-hoc analysis of individual metabolite ratios using paired t-tests.

https://doi.org/10.1371/journal.pone.0198053.t002

Metabolite concentration changes

Using repeated measures MANOVA, there was no main effect of bihemispheric M1-SMA tDCS on the absolute concentration of any metabolite (F(1,14) = 1.55, p = 0.23). There was also no significant interaction effect of metabolite and condition (F(4,11) = 1.42, p = 0.29). Table 3 displays the absolute metabolite concentrations. Post-hoc comparisons using paired t-tests (uncorrected for multiple comparisons) showed a trend toward decreased tCr (p = 0.07, Cohen’s d = 0.42) and decreased mI (p = 0.08, Cohen’s d = 0.48), with no significant changes in any metabolites.

[Figure omitted. See PDF.]

Table 3. Analysis of absolute metabolic concentration: Absolute concentration of metabolites measured by MRS.

The p-values were measured by post-hoc analysis of individual metabolites using paired t-tests.

https://doi.org/10.1371/journal.pone.0198053.t003

Correlation between NAA and tCr

We observed a strong, positive correlation between the change in the absolute concentration of NAA and the change in the absolute concentration of tCr (stimulation–sham, R2 = 0.64, p < 0.001, Fig 4).

[Figure omitted. See PDF.]

Fig 4. NAA and tCr coupling.

The association between the Δ in the absolute concentration of NAA and tCr (stimulation − sham). We observed a strong, positive correlation, indicating NAA and tCr both change in the same direction after stimulation (R2 = 0.64, p <0.001).

https://doi.org/10.1371/journal.pone.0198053.g004

Discussion

Bihemispheric tDCS was safely and successfully performed in the 7T MRI environment with minimal heating effects. However, the bihemispheric tDCS of M1-SMA produced no significant metabolite level changes in the left primary motor cortex (M1) immediately after 20 minutes of stimulation measured by 7T MR spectroscopy. Post-hoc analysis did show a trend toward increased NAA/tCr and decreased levels of tCr. In addition, we observed a strong association between the change in absolute concentration of NAA and the change in absolute concentration of tCr that may indicate a coupling between these metabolites following tDCS.

The trends toward lower NAA/tCr and tCr observed in the current study may be due to changes in brain activity induced by tDCS. Brain activity has been shown to decrease both NAA and tCr levels. Specifically, NAA is associated with metabolic and mitochondrial activity [16, 41]. Following visual stimulation, Baslow and colleagues found that the concentration of NAA decreased by approximately 13% in the visual cortex [42]. Similarly, Castellano and colleagues observed a 20% decrease in NAA after visual stimulation [43]. This decrease in NAA was attributed to a lower rate of NAA synthesis compared to hydrolysis during periods of cortical activation, suggesting that the brain used NAA faster than it could be synthesized [42, 43]. NAA is the precursor for the synthesis of N-acetylaspartylglutamate (NAAG), a modulator of glutamate and GABA neurotransmitter release. When neural activity is increased, there is an increased release of NAAG from the synapse [43, 44]. It has been suggested that the reduction in NAA upon neural activation is due to increased demand for NAAG. In support of this hypothesis, both Landim et al. and Castellano et al. observed a decrease in NAA concentration with a subsequent increase in NAAG upon stimulation [43, 44].

Creatine (Cr) may also be altered as a consequence of neuronal stimulation due to its role in energy metabolism through its conversion to phosphocreatine (PCr) [15, 45]. In the central nervous system, Cr and PCr are involved in maintaining the high energy levels necessary for the maintenance of membrane potentials, ion gradients, calcium homeostasis, and intracellular signalling [46]. Cr has also been observed as a potential modulator of neurotransmission [15, 45]. The Cr peak measured by MR spectroscopy represents intracellular contributions from both Cr and PCr (tCr). Rango et. al. have also discovered a transient decrease in PCr after short bursts of visual stimulation, concluding that functional activation reduces PCr [47]. Furthermore, Cr is released from the neuron in an action potential dependent manner. An increase in the resting membrane potential, induced by tDCS, may result in release of Cr from the neuron to act as a co-transmitter. Early studies on rodents indicate Cr may modulate postsynaptic neurotransmitters such as GABA, inhibiting its action [48, 49]. Release of Cr from intracellular stores would decrease the measureable concentration of Cr.

A trend towards a decrease in absolute concentration of myo-inositol was also observed in the current study. To our knowledge only one other study examined a change in myoinositol after tDCS. Contrary to our results, Rango et al. observed an increase in the concentration of myo-inositol 30 minutes after 1.5 mA of anodal tDCS to the right motor cortex [12]. The use of different current levels, electrode montage, and area measured, could all contribute to this discrepancy.

The observed association between ΔNAA and ΔtCr indicates that individuals that had a decrease in NAA following tDCS relative to sham also had a decrease in tCr levels (Fig 4). One possible explanation for the correlation between NAA and tCr is related to the neuronal response to excitatory stimulation. Specifically, in response to excitatory potentials, glutamate, NAAG, and creatine are released from the neuron in a Ca2+ dependent manner [15, 44]. A change in visible tCr may occur for two reasons. First, when creatine is released, it acts on GABAa receptors and consequently MRS visibility may be reduced. Second, in areas of high energy demand, such as neuronal activation, PCr is used to convert adenosine diphosphate (ADP) to adenosine triphosphate (ATP), reducing intracellular stores of PCr and increasing Cr [15, 17, 45, 49]. A change in NAA may also occur for several reasons. First, NAA may be used to synthesize NAAG for neurotransmission. Second, recent studies have linked NAA to mitochondrial energy metabolism and indirectly to neuronal activation [50]. Third, in a study by Yan et al., NAA was found to be involved in excitation of the neuron by acting on glutamate receptors (mGluR) to induce an inward current [51]. Although the current study was not designed to elucidate the biological basis for the observed association between NAA and tCr, previous studies have shown that both metabolites do have roles in neuronal excitation.

The after effects of tDCS are thought to be dependent on alterations of the membrane potential and changes in glutamate and GABA signalling, relating to synaptic plasticity [3, 6]. As such, we expected to observe changes in glutamate and GABA following tDCS. However, the literature presents conflicting results [7, 10–12, 17, 19, 52, 53]. In a study observing metabolite concentration both during and after tDCS, Bachtiar et. al. observed a significant decrease of GABA concentration in left M1 after anodal tDCS to the same area compared to sham, but no differences between sham and anodal tDCS during the stimulation period [52]. This indicates that the alteration of neurotransmitters that occurs due to tDCS is predominantly evident after the stimulation period. It is possible that our measurement of metabolite concentration occurred outside the optimal window of neurotransmitter modulation, and instead we observed upstream events. Further studies are required to identify the mechanism by which GABA and Glu are altered to enhance or depress synaptic activity and the optimal time to observe the peak change in these neurotransmitters.

The current study is the first to measure metabolite concentrations of the motor cortex using a bihemispheric montage in an ultra-high-field MRI (7T). Currently, there are only two studies that have examined the metabolism of the motor cortex following tDCS at 7T, both using the conventional M1-supraorbital (unihemispheric) montage, and both stimulating outside of the scanner [7, 18]. Both studies examined the effects of 1 mA of cathodal stimulation over left M1 for 15 [18] and 10 [7] minutes with differing results. Stagg et. al. found a decrease in Glu/Cr after cathodal stimulation, while Kim and colleagues reported no significant change in Glu concentration following cathodal stimulation. Kim et. al. did report a significant reduction in GABA following anodal stimulation, and no changes in other key metabolites [18]. The current study applied 2 mA of current for 20 minutes. It has recently been shown that cathodal stimulation, which is believed to be inhibitory, reverses its polarity at 2 mA and becomes excitatory [9]. The higher current used in our study compared to previous studies may explain the differing results. Increasing the current to 2 mA delays the time of peak metabolic change from immediately after stimulation, to up to 90–120 minutes after stimulation [9]. This delay may explain why we did not observe changes in Glu and GABA in the current study.

Only one other study has examined bihemispheric (M1-M1) tDCS on motor cortex metabolism [19]. Using this montage, and 1mA of current for 20 minutes, Tremblay and colleagues reported no significant modulation in any metabolite concentration at 3T [19], consistent with the current study. They concluded that the complex relationship between excitatory and inhibitory mechanisms within and between the primary motor cortices resulted in high inter individual variability and response to tDCS stimulation. The ultra-high field MRS used in the current study provided greater signal to noise ratio and spectral dispersion compared to Tremblay et al.,[19] increasing metabolite measurement precision. This greater measurement precision may be responsible for the observation of trends towards increased NAA/tCr and decreased tCr in the current study.

Although there have been few studies observing the metabolic and functional changes following M1-M1 tDCS, the current study is the first to incorporate the SMA as a potential target for bihemispheric tDCS. The SMA has been studied as a potential tDCS target in behavioural studies of posture and visuomotor learning [28, 54]. Its anatomical positioning and strong connections to M1 make it a well-suited target for motor network modulation. In addition, the SMA has strong efferent connections to the corticospinal tract, making it an ideal candidate as a target for neurorehabilitation [27]. Various studies have shown the importance of the SMA and associated non primary motor areas after brain or spinal cord injury [24–26]. Neural recruitment is an important aspect of recovery. SMA should be considered to enhance synaptic connections of bilateral M1, subcortical structures, and further downstream to the corticospinal tract.

One important limitation of the current study was the omission of a within session baseline measurement. The MRS data acquired for this study was part of a longer imaging protocol that incorporated anatomical and resting-state fMRI measurements (to be published elsewhere). Therefore, time constraints prevented the inclusion of a baseline spectroscopy measurement. However, the cross-over design of this study including a separate spectroscopy measurement during sham stimulation on a separate day has previously been shown to be an acceptable approach [30]. However, the inclusion of a baseline measurement in the future would likely reduce inter-subject variability. A recent study at 7T has provided estimates of the reliability of metabolite measurements taken on separate days (31). In addition, a comparable protocol, using a sham control without baseline measurements observed no changes in metabolite concentration after 1 mA of bilateral dorsolateral prefrontal cortex (DLPFC), measured from the left DLPFC [11]. Another limitation of the current study is that metabolite measurements were not made in the SMA, again due to time constraints. Future studies would also benefit from examining metabolite changes in this brain region following stimulation. Finally, the current study was not designed to examine GABA levels. Although, previous studies have shown changes in GABA concentration, we were not able to measure GABA with sufficient reproducibility using the spectroscopy method applied in the current study. Future studies using GABA editing methods could help elucidate the modulation of GABA by tDCS.

In conclusion, bihemispheric transcranial direct current stimulation with anode over SMA and cathode over M1 was safely applied during 7T MRI for 20 minutes at 2 mA. Immediately following stimulation there were no changes in metabolite levels measured by 1H MR spectroscopy of the left primary motor cortex in this sham controlled cross-over study. However, when comparing stimulation to sham conditions, there was a significant positive association between the change in N-acetyl aspartate and the change in creatine in the same region.

Citation: Ryan K, Wawrzyn K, Gati JS, Chronik BA, Wong D, Duggal N, et al. (2018) 1H MR spectroscopy of the motor cortex immediately following transcranial direct current stimulation at 7 Tesla. PLoS ONE 13(8): e0198053. https://doi.org/10.1371/journal.pone.0198053

References

1. Nitsche MA, Paulus W. Excitability changes induced in the human motor cortex by weak transcranial direct current stimulation. J Physiol. 2000;527 Pt 3:633–9. Epub 2000/09/16. pmid:10990547; PubMed Central PMCID: PMCPMC2270099.

2. Nitsche MA, Paulus W. Sustained excitability elevations induced by transcranial DC motor cortex stimulation in humans. Neurology. 2001;57(10):1899–901. Epub 2001/11/28. pmid:11723286.

3. Nitsche MA, Paulus W. Transcranial direct current stimulation—update 2011. Restorative Neurology and Neuroscience. 2011;29(6):463–92. pmid:22085959

4. Stagg CJ, O’Shea J, Kincses ZT, Woolrich M, Matthews PM, Johansen-Berg H. Modulation of movement-associated cortical activation by transcranial direct current stimulation. European Journal of Neuroscience. 2009;30(7):1412–23. pmid:19788568

5. Yoon KJ, Oh BM, Kim DY. Functional improvement and neuroplastic effects of anodal transcranial direct current stimulation (tDCS) delivered 1 day vs. 1 week after cerebral ischemia in rats. Brain Research. 2012;1452:61–72. pmid:22444278

6. Stagg CJ, Nitsche MA. Physiological basis of transcranial direct current stimulation. Neuroscientist. 2011;17(1):37–53. Epub 2011/02/24. pmid:21343407.

7. Stagg CJ, Best JG, Stephenson MC, O'Shea J, Wylezinska M, Kincses ZT, et al. Polarity- sensitive modulation of cortical neurotransmitters by transcranial stimulation. J Neurosci. 2009;29(16):5202–6. Epub 2009/04/24. pmid:19386916.

8. Stagg CJ, Bachtiar V, Johansen-Berg H. The role of GABA in human motor learning. Curr Biol. 2011;21(6):480–4. Epub 2011/03/08. pmid:21376596; PubMed Central PMCID: PMCPMC3063350.

9. Batsikadze G, Moliadze V, Paulus W, Kuo MF, Nitsche MA. Partially non-linear stimulation intensity-dependent effects of direct current stimulation on motor cortex excitability in humans. J Physiol. 2013;591(7):1987–2000. Epub 2013/01/23. pmid:23339180; PubMed Central PMCID: PMCPMC3624864.

10. Clark VP, Coffman BA, Trumbo MC, Gasparovic C. Transcranial direct current stimulation (tDCS) produces localized and specific alterations in neurochemistry: a (1)H magnetic resonance spectroscopy study. Neurosci Lett. 2011;500(1):67–71. Epub 2011/06/21. pmid:21683766.

11. Hone-Blanchet A, Edden RA, Fecteau S. Online Effects of Transcranial Direct Current Stimulation in Real Time on Human Prefrontal and Striatal Metabolites. Biol Psychiatry. 2016;80(6):432–8. Epub 2016/01/18. pmid:26774968.

12. Rango M, Cogiamanian F, Marceglia S, Barberis B, Arighi A, Biondetti P, et al. Myoinositol content in the human brain is modified by transcranial direct current stimulation in a matter of minutes: a 1H-MRS study. Magn Reson Med. 2008;60(4):782–9. Epub 2008/09/26. pmid:18816828.

13. Stagg CJ, Jayaram G, Pastor D, Kincses ZT, Matthews PM, Johansen-Berg H. Polarity and timing-dependent effects of transcranial direct current stimulation in explicit motor learning. Neuropsychologia. 2011;49(5):800–4. Epub 2011/02/22. pmid:21335013; PubMed Central PMCID: PMCPMC3083512.

14. Stagg Charlotte J, Bachtiar V, Johansen-Berg H. The Role of GABA in Human Motor Learning. Elsevier; 2011. p. 480–4.

15. Beard E, Braissant O. Synthesis and transport of creatine in the CNS: importance for cerebral functions. J Neurochem. 2010;115(2):297–313. Epub 2010/08/28. pmid:20796169.

16. Rae CD. A Guide to the Metabolic Pathways and Function of Metabolites Observed in Human Brain 1H Magnetic Resonance Spectra. Neurochem Res. 2013;39(1):1–36 pmid:24258018

17. Rae CD, Lee VH, Ordidge RJ, Alonzo A, Loo C. Anodal transcranial direct current stimulation increases brain intracellular pH and modulates bioenergetics. Int J Neuropsychopharmacol. 2013;16(8):1695–706. Epub 2013/03/12. pmid:23473040.

18. Kim S, Stephenson MC, Morris PG, Jackson SR. tDCS-induced alterations in GABA concentration within primary motor cortex predict motor learning and motor memory: a 7 T magnetic resonance spectroscopy study.Neuroimage.2014; 99:237–243 pmid:24904994

19. Tremblay S, Lafleur LP, Proulx S, Beaule V, Latulipe-Loiselle A, Doyon J, et al. The effects of bi-hemispheric M1-M1 transcranial direct current stimulation on primary motor cortex neurophysiology and metabolite concentration. Restor Neurol Neurosci. 2016;34(4):587–602. Epub 2016/05/28. pmid:27232951.

20. Waters S, Wiestler T, Diedrichsen J. Cooperation Not Competition: Bihemispheric tDCS and fMRI Show Role for Ipsilateral Hemisphere in Motor Learning. J Neurosci. 2017;37(31):7500–12. Epub 2017/07/05. pmid:28674174; PubMed Central PMCID: PMCPMC5546115.

21. Waters-Metenier S, Husain M, Wiestler T, Diedrichsen J. Bihemispheric transcranial direct current stimulation enhances effector-independent representations of motor synergy and sequence learning. Journal of Neuroscience. 2014;34(3):1037–50. pmid:24431461

22. Floel A. tDCS-enhanced motor and cognitive function in neurological diseases. Neuroimage. 2014;85 Pt 3:934–47. Epub 2013/06/04. pmid:23727025.

23. Grefkes C, Eickhoff SB, Nowak DA, Dafotakis M, Fink GR. Dynamic intra- and interhemispheric interactions during unilateral and bilateral hand movements assessed with fMRI and DCM. Neuroimage. 2008;41(4):1382–94. Epub 2008/05/20. pmid:18486490.

24. Carey LM, Abbott DF, Egan GF, Bernhardt J, Donnan GA. Motor impairment and recovery in the upper limb after stroke: behavioral and neuroanatomical correlates. Stroke. 2005;36(3):625–9. Epub 2005/01/29. pmid:15677574.

25. Wong WW, Chan ST, Tang KW, Meng F, Tong KY. Neural correlates of motor impairment during motor imagery and motor execution in sub-cortical stroke. Brain Inj. 2013;27(6):651–63. Epub 2013/03/22. pmid:23514275.

26. Ryan K, Goncalves S, Bartha R, Duggal N. Motor Network Recovery in Patients with Chronic Spinal Cord Compression: a Longitudinal Study Following Decompression Surgery J Neurosurg: Spine. 2018; 28(4):379–388. Epub 2018 Jan 19. pmid:29350595

27. Luppino G, Matelli M, Camarda R, Rizzolatti G. Corticocortical connections of area F3 (SMA-proper) and area F6 (pre-SMA) in the Macaque monkey. Journal of Comparative Neurology. 1993;338(1):114–40. pmid:7507940

28. Hupfeld KE, Ketcham CJ, Schneider HD. Transcranial direct current stimulation (tDCS) to the supplementary motor area (SMA) influences performance on motor tasks. Exp Brain Res. 2017;235(3):851–9. Epub 2016/12/03. pmid:27909747.

29. Jurcak V, Tsuzuki D, Dan I. 10/20, 10/10, and 10/5 systems revisited: their validity as relative head-surface-based positioning systems. Neuroimage. 2007;34(4):1600–11. Epub 2007/01/09. pmid:17207640.

30. Gandiga PC, Hummel FC, Cohen LG. Transcranial DC stimulation (tDCS): a tool for double-blind sham-controlled clinical studies in brain stimulation. Clin Neurophysiol. 2006;117(4):845–50. Epub 2006/01/24. pmid:16427357.

31. Oz G, Tkac I. Short-echo, single-shot, full-intensity proton magnetic resonance spectroscopy for neurochemical profiling at 4 T: validation in the cerebellum and brainstem. Magn Reson Med. 2011;65(4):901–10. Epub 2011/03/18. pmid:21413056; PubMed Central PMCID: PMCPMC3827699.

32. Schar M, Kozerke S, Boesiger P. Navigator gating and volume tracking for double- triggered cardiac proton spectroscopy at 3 Tesla. Magnetic Resonance in Medicine2004. p. 1091–5. pmid:15170826

33. Bartha R, Drost DJ, Menon RS, Williamson PC. Spectroscopic lineshape correction by QUECC: combined QUALITY deconvolution and eddy current correction. Magn Reson Med. 2000;44(4):641–5. Epub 2000/10/12. pmid:11025521.

34. Bartha R, Drost DJ, Williamson PC. Factors affecting the quantification of short echo in- vivo 1H MR spectra: prior knowledge, peak elimination, and filtering. NMR Biomed. 1999;12(4):205–16. Epub 1999/07/28. pmid:10421912.

35. Rupsingh R, Borrie M, Smith M, Wells JL, Bartha R. Reduced hippocampal glutamate in Alzheimer disease. Neurobiol Aging. 2011;32(5):802–10. Epub 2009/06/09. pmid:19501936.

36. Jenkinson M, Beckmann CF, Behrens TE, Woolrich MW, Smith SM. FSL. Neuroimage. 2012;62(2):782–90. Epub 2011/10/08. pmid:21979382.

37. Andreychenko A, Klomp DW, de Graaf RA, Luijten PR, Boer VO. In vivo GABA T2 determination with J-refocused echo time extension at 7 T. NMR Biomed. 2013;26(11):1596–601. Epub 2013/07/31. pmid:23893556.

38. Marjanska M, Auerbach EJ, Valabregue R, Van de Moortele PF, Adriany G, Garwood M. Localized 1H NMR spectroscopy in different regions of human brain in vivo at 7 T: T2 relaxation times and concentrations of cerebral metabolites. NMR Biomed. 2012;25(2):332–9. Epub 2011/07/29. pmid:21796710; PubMed Central PMCID: PMCPMC3357544.

39. Xin L, Schaller B, Mlynarik V, Lu H, Gruetter R. Proton T1 relaxation times of metabolites in human occipital white and gray matter at 7 T. Magn Reson Med. 2013;69(4):931–6. Epub 2012/06/01. pmid:22648904.

40. Kreis R, Slotboom J, Hofmann L, Boesch C. Integrated data acquisition and processing to determine metabolite contents, relaxation times, and macromolecule baseline in single examinations of individual subjects. Magn Reson Med. 2005;54(4):761–8. Epub 2005/09/15. pmid:16161114.

41. Patel T, Blyth JC, Griffiths G, Kelly D, Talcott JB. Moderate relationships between NAA and cognitive ability in healthy adults: implications for cognitive spectroscopy. Front Hum Neurosci. 2014;8:39. Epub 2014/03/05. pmid:24592224; PubMed Central PMCID: PMCPMC3924143.

42. Baslow MH, Hrabe J, Guilfoyle DN. Dynamic relationship between neurostimulation and N-acetylaspartate metabolism in the human visual cortex: evidence that NAA functions as a molecular water pump during visual stimulation. J Mol Neurosci. 2007;32(3):235–45. Epub 2007/09/18. pmid:17873369.

43. Castellano G, Dias C, Foerster B, Li L, Covolan R. NAA and NAAG variation in neuronal activation during visual stimulation. Braz J Med Biol Res. 452012. p. 1031–6.

44. Landim RCG. Investigation of NAA and NAAG dynamics underlying visual stimulation using MEGA-PRESS in a functional MRS experiment. 2016;34(3):239–45. pmid:26656908; PubMed Central PMCID: PMCPMC5510168.

45. Royes LF, Fighera MR, Furian AF, Oliveira MS, Fiorenza NG, Ferreira J, et al. Neuromodulatory effect of creatine on extracellular action potentials in rat hippocampus: role of NMDA receptors. Neurochem Int. 2008;53(1–2):33–7. Epub 2008/06/11. pmid:18541342.

46. Wyss M, Kaddurah-Daouk R. Creatine and creatinine metabolism. Physiol Rev. 2000;80(3):1107–213. Epub 2000/07/14. pmid:10893433.

47. Rango M, Castelli A, Scarlato G. Energetics of 3.5 s neural activation in humans: a 31P MR spectroscopy study. Magn Reson Med. 1997;38(6):878–83. Epub 1997/12/24. pmid:9402187.

48. Deyn PPD, Macdonald RL. Guanidino compounds that are increased in cerebrospinal fluid and brain of uremic patients inhibit GABA and glycine responses on mouse neurons in cell culture. Annals of Neurology. 1990;28(5):627–33. pmid:2135950

49. Almeida LS, Salomons GS, Hogenboom F, Jakobs C, Schoffelmeer ANM. Exocytotic release of creatine in rat brain. Synapse. 2006;60(2):118–23. pmid:16715490

50. Moffett JR, Ross B, Arun P, Madhavarao CN, Namboodiri AM. N-Acetylaspartate in the CNS: from neurodiagnostics to neurobiology. Prog Neurobiol. 2007;81(2):89–131. Epub 2007/02/06. pmid:17275978; PubMed Central PMCID: PMCPMC1919520.

51. Yan HD, Ishihara K, Serikawa T, Sasa M. Activation by N-acetyl-L-aspartate of acutely dissociated hippocampal neurons in rats via metabotropic glutamate receptors. Epilepsia. 2003;44(9):1153–9. Epub 2003/08/16. pmid:12919386.

52. Bachtiar V, Near J, Johansen-Berg H, Stagg CJ. Modulation of GABA and resting state functional connectivity by transcranial direct current stimulation. Elife. 2015;4:e08789. Epub 2015/09/19. pmid:26381352; PubMed Central PMCID: PMCPMC4654253.

53. Horvath JC, Forte JD, Carter O. Evidence that transcranial direct current stimulation (tDCS) generates little-to-no reliable neurophysiologic effect beyond MEP amplitude modulation in healthy human subjects: A systematic review. Neuropsychologia. 2015;66:213–36. Epub 2014/12/03. pmid:25448853.

54. Vollmann H, Conde V, Sewerin S, Taubert M, Sehm B, Witte OW, et al. Anodal transcranial direct current stimulation (tDCS) over supplementary motor area (SMA) but not pre-SMA promotes short-term visuomotor learning. Brain Stimul. 2013;6(2):101–7. Epub 2012/06/05. pmid:22659022.