Abstract
For human protection, the internal electric field is used as a dosimetric quantity for electromagnetic fields lower than 5–10 MHz. According to international standards, in this frequency range, electrostimulation is the main adverse effect against which protection is needed. One of the topics to be investigated is the quantification of the internal electric field threshold levels of perception and pain. Pain has been reported as a side effect during transcranial magnetic stimulation (TMS), especially during stimulation of the Broca’s (speech) area of the brain. In this study, we designed an experiment to conduct a dosimetry analysis to quantify the internal electric field corresponding to perception and pain thresholds when targeting the Broca’s and M1 areas from magnetic stimulator exposure. Dosimetry analysis was conducted using a multi-scale analysis in an individualized head model to investigate electrostimulation in an axonal model. The main finding is that the stimulation on the primary motor cortex has higher perception and pain thresholds when compared to Broca’s area. Also, TMS-induced electric field applied to Broca’s area exhibited dependence on the coil orientation at lower electric field threshold which was found to be related to the location and thickness of pain fibers. The derived dosimetry quantities provide a scientific rationale for the development of human protection guidelines and the estimation of possible side effects of magnetic stimulation in clinical applications.
Introduction
There has been concern about human safety under exposure to electromagnetic fields. To protect humans from electromagnetic exposure, safety guidelines/standards have been developed by international standardization bodies (; ; ). In the guidelines/standards for exposure to up to 100 kHz (5–10 MHz for brief pulse exposures), electrostimulation is the primary effect against which protection is needed. The IEEE ICES has published a research agenda for low-frequency exposure (). One of the items for the peripheral nervous system is target tissue for electrostimulation.
One difficulty to conduct assessment of protection limits is that exposure levels from conventional appliances are well below the limit prescribed in the guidelines/standards and may not cause any stimulation at all. However, electrostimulation can occur during medical treatment; one example is transcranial magnetic stimulation (TMS), which is often used for the diagnosis of brain functions, neuro-rehabilitation, therapy for depression, and so on (; ; ; ; ; ; ).
During magnetic stimulation, pain has been reported as a side effect, especially during stimulation of the Broca’s (speech) area of the brain. Pain thresholds for both Broca’s area and M1 are different and significantly lower than the motor threshold as measured via motor evoked potentials (MEPs; ), which suggests that experiments might evoke pain sensation at the site of stimulation. This side effect generally becomes relevant when the intended target is in deep regions (; ; ). Safety and recommendations for TMS have been also published ().
In this study, we designed an experiment to explore the perception and pain thresholds when targeting the Broca’s and M1 areas for exposure to magnetic stimulators. The corresponding internal electric fields in the skin and muscles were then computed. In addition, multi-scale modeling, i.e., the combination of electromagnetics and nerve activation modeling, has been used to clarify the stimulation site on the tissue.
Materials and Methods
Participants
Twelve participants were recruited for each of the BA (21.2 ± 0.8 years, seven male and five female) and M1 (23.2 ± 3.1 years, nine male and three female). Four participants took part in both conditions. None of the participants had any contraindications to TMS, took any medication on a regular basis, or had a history of psychiatric or neurological diseases through questionnaires. Written informed consent was obtained from all participants before their participation. The study was approved by the ethical committee of the Hamamatsu University School of Medicine and was conducted in accordance with the Declaration of Helsinki.
Magnetic Resonance Imaging
The head models of the 12 participants were constructed from T1- and T2-weighted images with acquisition parameters as follows: T1 MPRAGE sequence with TR/TE/FA/FOV/voxel size/slice number = 7.172 ms/2.12 ms/15°/256 mm/1.0 mm × 1.0 mm × 1.0 mm/196, and T2 with TR/TE/FOV/voxel size/slice number = 2502 ms/76.404 ms/256 mm/1.0 mm × 1.0 mm × 1.0 mm/196.
Anatomical Head Model
The volume conductor of the head model was obtained by estimating the electrical conductivity values of the brain and non-brain tissues by applying deep learning based on MR images, as shown in Figure 1. The head tissues were segmented using in-house software equipped with the FreeSurfer brain imaging software package, as described previously (), and the FreeSurfer image analysis software () was used to reconstruct the surfaces of the gray and white matter. Non-brain tissues were segmented from T1- and T2-weighted MRI using a semi-automatic procedure of region-growing and thresholding techniques.
FIGURE 1
Measurement Protocol
The procedure to detect perception and pain thresholds was based on our previous study (). In brief, all subjects were awake, sat comfortably in a reclining chair in a quiet environment, and were requested to relax. Single-pulse TMS was applied with a Magstim 2002 magnetic stimulator (Magstim Co., United Kingdom) connected to a double alpha BI coil (60 mm, Magstim Co., Ltd., United Kingdom). Perception/pain thresholds derived at two stimulation sites corresponding to stimulation of left Brodmann area 44 for Broca’area and the center of the hand knob () area on the Brodmann area 4 for M1. The cortical areas were anatomically identified by the construction of a 3D cortical surface model of the individual participants using a frameless navigation system (Brainsight, Rogue Research Inc, Canada) based on the individual T1 image. A total of seven coil rotations were measured with approximately 5 min of rest between coil-orientation conditions for each target region. Figure 2 shows the definition of coil orientations over both target regions.
FIGURE 2
The subjects were instructed to report the presence or absence of scalp perception/pain after each stimulation. For scalp perception, the participants were instructed to report the presence or absence of a sensation of pressure or force. For the perception of pain, the participants were instructed to report the presence or absence of scalp-pain after each stimulation, regardless of the magnitude or type of pain. We used an adaptive staircase method () to determine the thresholds. The intensity decreased when perception/pain was reported and increased when perception/pain was not reported. A perception/pain threshold was defined as the minimum intensity that induced perception/pain in at least five of the ten trials. The thresholds are given as the percentage of the maximum stimulation output of the device (%MSO). This procedure was repeated for each coil orientation and the target area (M1 and Broca’s).
Each participant completed the experiment on two separate days. The perceptual/pain thresholds for the seven coil-orientation conditions were measured in a randomized order on the first day. The order of the coil-orientation conditions on the second day was set in reverse to the first day. For each coil-orientation condition, perceptual/pain thresholds were averaged from the two days to reduce the potential influence of sensory adaptation or fatigue. The total number of stimulations varied between 250 and 450 stimuli depending on the participants.
Electromagnetic Computation
We assumed that the electric displacement current is negligible when compared to the conduction current (magneto-quasi-static approximation) and that the induced current does not perturb the external magnetic field (; ; ). First, the electric scalar potential induced in the brain ϕ was determined using the following equation:
where σ is the electric conductivity, as a scalar piecewise constant conductivity, using the same values in Aonuma et al. (2018). The time derivative of the magnetic vector potential was determined from the model using the thin-wire approximation and the Biot–Savart law. The coil consists of two wings with eleven equally spaced concentric current loops mimicking the coil in the experiment. The inner and outermost loops of the coil were 1.75 and 3.8 cm in diameter, respectively. Second, the induced electric field was calculated as follows:
The computed induced EF and scalar electric potential corresponded to temporal peak values at the coil operating frequency of 3 kHz (). Eq. (1) was solved numerically by the finite-element method with first-order cubical elements (0.488-mm voxel of the volume conductor model), using an in-house software, as described in , which has been applied in different TMS studies (; ; ).
We used an anatomical volume conductor model for each participant to compute the individualized induced electric fields. The coil position and orientation and stimulation intensity in the simulation were the same as those in the experiments. Coil orientation and position were recorded using a navigation system.
Nerve Model
The effects of the computed extracellular electric scalar field on nerve model axons placed over the Broca’s area were investigated. The smallest stimulator output intensity was obtained to elicit an action potential in a conductance-based nerve model using the following equation:
where cm is the membrane capacitance, Vm,n is the membrane potential at position n along the axon, and the variable R denotes the intra-axonal resistance between the centers of two adjacent compartments. The list of values of the nerve parameters can be found in . The spatial structure of the myelinated neuron consists of internodes (segments ensheathed by myelin) concatenated with nodes of Ranvier (ionic channels). The ionic current is passive in the internodes and nonlinear in the nodes, as described by the Chiu–Ritchie–Rogart–Stagg–Sweeney model ().
The right-hand side of Eq. 3 corresponds to Δ2ϕ = ϕ(n−1)−2ϕ(n) + ϕ(n + 1), which describes the driving term of the activation. To incorporate the total induced electric field E of Eq. 2 in Eq. 3, a local linear integral of the electric field along the trajectory of the cable model is used to obtain the extracellular driving potential ϕ :
The terminal condition used for the proximal axon node at the seed point was clamped at 0 V.
Results
Experimental Pain and Perception Thresholds
In Figure 3, the measured average pain and perceptual thresholds are shown for different coil orientations. The threshold is given by the percentage of the maximum stimulation output of the device (MSO%). One-way analysis of variance (ANOVA) with repeated measures revealed a significant main effect of coil angle for the Broca’s area (perception, F6, 66 = 3.58; p < 0.01; pain, F6, 66 = 4.38; p < 0.001), but not for the M1 region (perception, F6, 66 = 1.21; p = 0.31; pain, F6, 66 = 1.40; p = 0.26). We found that pain and perceptual thresholds were significantly lower in the 60° than the 180° conditions (Bonferonni post-hoc tests; p < 0.05). No significant differences were not observed between the other conditions. These results indicate a coil angle dependency only for TMS targeting the Broca’s area. The pain/perception thresholds were not the same for the target areas. In addition, the pain threshold (%MSO) was approximately twice that of perception.
FIGURE 3
Dosimetry of Pain and Perception Thresholds
Computational electromagnetic dosimetry was conducted to determine the induced electric field using the same stimulation conditions in the experiment (MSO, coil position/orientation, and individual head model). Figure 4 illustrates the induced electric field in the skin for different coil positions over the Broca’s and M1 areas.
FIGURE 4
To investigate the consistency of the coil dependency of the measured thresholds, the internal electric field in the skin and muscle was computed using the same experimental conditions (measured MSO% threshold, coil orientation/position, and individualized head models) for each subject. As shown in Figure 5, the dependency on the coil orientation of the induced electric showed similar trends in the Broca’s and M1 areas to the measured thresholds (Figure 3).
FIGURE 5
The induced electric field strength corresponding to the pain threshold varies between 85 and 104 V/m in skin tissue and 64 and 79 V/m in the muscle for stimulation over the Broca’s area. In addition, the perception threshold varies between 46 and 53 V/m in the skin and 35 and 42 V/m in the muscle for stimulation over the Broca’s area. In the case of M1, the induced electric field thresholds were higher than those of the Broca’s area. The variation of perception thresholds between Broca’s area and M1 indicates that different perception/pain fibers are involved.
Nerve Stimulation in Broca’s Area
We also derived the stimulation thresholds directly from multiscale modeling by incorporating a nerve model in the muscle over the Broca’s area (Supplementary Table 1). We then investigated the most effective orientation of the fiber below the Broca’s area and analyzed the coefficient of determination between induced electric field thresholds based on a nerve stimulation model and computational dosimetry based on the experimental measurements. We identified one effective orientation of the nerve model (45°) in which the coefficient of determination reaches the maximum value, as shown in Figure 6. For this orientation, the axon threshold lies between 31 and 79 V/m, which agrees with the dosimetry analysis of perception thresholds in the muscle (Figure 5). This indicates that pain thresholds may be derived from not only the skin tissue but also deeper tissues in Broca’s area where muscle tissue exists.
FIGURE 6
Discussion
There is a large body of evidence that motor and speech representations in the primary motor cortex and Broca’s area have an optimal coil rotation angle (i.e., the direction of the induced current), which depends on individual neuroanatomy (; ; ; Sollmann et al., 2015; ). Pain thresholds for both Broca’s area and M1 are significant differently (). Moreover, the pain threshold shows a significant difference between coil orientations in the Broca’s area in contrast to M1 (). This study conducted a dosimetry analysis of the induced electric field in non-brain tissues related to perception and pain thresholds for the first time to quantify the internal values and clarified differences observed between Broca’s area and M1. The computation of the induced electric field is based on the same conditions of TMS experiments to detect perception/pain thresholds when targeting the Broca’s and M1 areas, including the head model considering subject anatomy. In addition, multiscale analysis was performed to investigate the activation thresholds of neural fibers in non-brain tissue.
We confirmed that the induced electric field thresholds for pain and perception depend on the coil orientation for Broca’s area stimulation, but not for M1 stimulation. Moreover, Broca’s area stimulation has smaller induced electric field thresholds than M1 stimulation. Both observations suggest that different perception/pain fibers are involved during stimulation of the Broca’s and M1 areas. First, the non-dependence of coil orientation during M1 stimulation could be related to the stimulation of small fibers (e.g., Aδ- and C-fibers) that are approximately normal to the skin surface. Second, smaller thresholds when targeting the Broca’s area could result because thicker fibers underneath the skin are stimulated with lower thresholds following the inverse relationship between fiber thickness and stimulation threshold in deeper tissues (e.g., muscle fibers). For the first argument, the induced electric fields on the skin during M1 stimulation present similar values to Aδ-fiber stimulation thresholds for perception in the range of 65 to 130 V/m, as found in . For the second argument, we characterized peripheral stimulation based on multiscale computation by modeling fibers on the muscle tissue that suggest fiber orientation similar to muscle found under the Broca’s area (e.g., temporoparietal muscle fibers). These observations indicate that pain thresholds may be derived from not only the skin tissue, but also deeper tissues, such as the muscle tissue when the Broca’s area is targeted. This was not the case for M1 due to the lack of muscle tissue, and because the perception/threshold is mostly driven by small fibers on the skin. The participants identified the pain location on the scalp as near below the center of the coil during Broca’areas stimulation, which discards pain originated from other areas, such as the eyes (K. Tani, personal communication, 2021).
Finally, the exposure level of the internal electric field in non-brain tissues was quantified for the first time in perception and pain conditions. The minimum thresholds were 40 V/m in the Broca’s area and 80 V/m in the M1 area. These values are useful when considering the perception and pain levels in standardization and medical applications. Studies aiming to achieve considerable induced electric field values in deep brain regions, such as deep TMS, need to consider the high induced fields in non-brain tissues. It has been shown that the electric field level in deep brain regions corresponds to a maximum of 50% of cortical values and 25% of scalp values (). This large difference indicates the importance of considering the side effects of stimulation and the importance of dosimetry in non-brain tissues. Factors affecting the perception/pain thresholds is the TMS coil design, stimulation waveform, pain assessment for evaluation of multiple levels, and interindividual differences. For the coil design, the quantification of internal electric field thresholds reduces its effects. As future work, we will expand on dosimetry for brain tissues.
Statements
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
The studies involving human participants were reviewed and approved by Hamamatsu University School of Medicine. The patients/participants provided their written informed consent to participate in this study.
Author contributions
AH and ST conceived and designed the study. ST and KT conducted the experiments. JG-T and KH conducted the simulation experiments. JG-T, KT, and KH processed the data. All authors analyzed the data, wrote the manuscript, and read and approved the manuscript.
Funding
This work was supported by the Ministry of Internal Affairs and Communications (grant number JPMI10001).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnins.2021.644951/full#supplementary-material
References
1
AonumaS.Gomez-TamesJ.LaaksoI.HirataA.TakakuraT.TamuraM.et al (2018). A high-resolution computational localization method for transcranial magnetic stimulation mapping.Neuroimage17285–93. 10.1016/j.neuroimage.2018.01.039
2
BarchanskiA.De GersemH.GjonajE.WeilandT. (2005). Impact of the displacement current on low-frequency electromagnetic fields computed using high-resolution anatomy models.Phys. Med. Biol.50N243–N249. 10.1088/0031-9155/50/19/N02
3
BarkerA. T.JalinousR.FreestonI. L. (1985). Non-invasive magnetic stimulation of human motor cortex.Lancet (London, Engl.)11106–1107. 10.1016/s0140-6736(85)92413-4
4
BashirS.PerezJ. M.HorvathJ. C.Pascual-LeoneA. (2013). Differentiation of motor cortical representation of hand muscles by navigated mapping of optimal TMS current directions in healthy subjects.J. Clin. Neurophysiol.30390–395. 10.1097/WNP.0b013e31829dda6b
5
DengZ. D.LisanbyS. H.PeterchevA. V. (2014). Coil design considerations for deep transcranial magnetic stimulation.Clin. Neurophysiol.1251202–1212. 10.1016/J.CLINPH.2013.11.038
6
FischlB. (2012). FreeSurfer.Neuroimage62774–781. 10.1016/j.neuroimage.2012.01.021
7
Gomez-TamesJ.HamasakaA.HirataA.LaaksoI.LuM.UenoS. (2020a). Group-level analysis of induced electric field in deep brain regions by different TMS coils.Phys. Med. Biol.65:025007. 10.1088/1361-6560/ab5e4a
8
Gomez-TamesJ.LaaksoI.MurakamiT.UgawaY.HirataA. (2020b). TMS activation site estimation using multiscale realistic head models.Iopscience.Iop.Org.17:036004. 10.1088/1741-2552/ab8ccf
9
Gomez-TamesJ.TarnaudT.MiwaK.HirataA.Van de SteeneT.MartensL.et al (2019). Brain cortical stimulation thresholds to different magnetic field sources exposures at intermediate frequencies.IEEE Trans. Electromagnetic Compat.611944–1952. 10.1109/TEMC.2019.2943138
10
GuggisbergA. G.DubachP.HessC. W.WüthrichC.MathisJ. (2001). Motor evoked potentials from masseter muscle induced by transcranial magnetic stimulation of the pyramidal tract: the importance of coil orientation.Clin. Neurophysiol. Off. J. Int. Federat. Clin. Neurophysiol.1122312–2319. 10.1016/s1388-2457(01)00677-0
11
HirataA.ItoF.LaaksoI. (2013). Confirmation of quasi-static approximation in SAR evaluation for a wireless power transfer system.Phys. Med. Biol.58N241–N249. 10.1088/0031-9155/58/17/N241
12
ICNIRP (1998). Guidelines for limiting exposure to time-varying electric, magnetic, and electromagnetic fields (up to 300 GHz). International Commission on Non-Ionizing Radiation Protection.Health Phys.74494–522.
13
IEEE Std C95.1-2019 (2019). IEEE Standard for Safety Levels with Respect to Human Exposure to Electric, Magnetic, and Electromagnetic Fields, 0 Hz to 300 GHz.Piscataway, NJ: IEEE.
14
International Commission on Non-Ionizing Radiation Protection (2020). Guidelines for limiting exposure to electromagnetic fields (100 KHz to 300 GHz).Health Phys.118483–524. 10.1097/hp.0000000000001210
15
LaaksoI.HirataA. (2012). Fast multigrid-based computation of the induced electric field for transcranial magnetic stimulation.Phys. Med. Biol.577753–7765. 10.1088/0031-9155/57/23/7753
16
LaaksoI.MurakamiT.HirataA.UgawaY. (2018). Where and what TMS activates: experiments and modeling.Brain Stimul.11166–174. 10.1016/J.BRS.2017.09.011
17
LaaksoI.TanakaS.KoyamaS.De SantisV.HirataA. (2015). Inter-subject variability in electric fields of motor cortical TDCS.Brain Stimul.8906–913. 10.1016/j.brs.2015.05.002
18
LernerA. J.WassermannE. M.TamirD. I. (2019). Seizures from transcranial magnetic stimulation 2012–2016: results of a survey of active laboratories and clinics.Clin. Neurophysiol.1301409–1416. 10.1016/j.clinph.2019.03.016
19
NieminenJ. O.KoponenL. M.IlmoniemiR. J. (2015). Experimental characterization of the electric field distribution induced by TMS devices.Brain Stimul.8582–589. 10.1016/j.brs.2015.01.004
20
Pascual-LeoneA.GatesJ. R.DhunaA. (1991). Induction of speech arrest and counting errors with rapid-rate transcranial magnetic stimulation.Neurology41697–702. 10.1212/WNL.41.5.697
21
Pascual-LeoneA.CammarotaA.WassermannE. M.Brasil-NetoJ. P.CohenL. G.HallettM. (1993). Modulation of motor cortical outputs to the reading hand of braille readers.Ann. Neurol.3433–37. 10.1002/ana.410340108
22
PlonseyR.HeppnerD. B. (1967). Considerations of quasi-stationarity in electrophysiological systems.Bull. Math. Biophys.29657–664. 10.1007/bf02476917
23
RaffinE.PellegrinoG.Di LazzaroV.ThielscherA.SiebnerH. R. (2015). Bringing transcranial mapping into shape: sulcus-aligned mapping captures motor somatotopy in human primary motor hand area.Neuroimage120164–175. 10.1016/j.neuroimage.2015.07.024
24
ReillyJ. P.HirataA. (2016). Low-frequency electrical dosimetry: research agenda of the IEEE International Committee on Electromagnetic Safety.Phys. Med. Biol.61R138–R149. 10.1088/0031-9155/61/12/R138
25
RossiS.AntalA.BestmannS.BiksonM.BrewerC.BrockmöllerJ.et al (2020). Safety and recommendations for TMS use in healthy subjects and patient populations, with updates on training, ethical and regulatory issues: expert guidelines.Clin. Neurophysiol.132269–306. 10.1016/j.clinph.2020.10.003
26
RossiniP. M.BurkeD.ChenR.CohenL. G.DaskalakisZ.Di IorioR.et al (2015). Non-invasive electrical and magnetic stimulation of the brain, spinal cord, roots and peripheral nerves: basic principles and procedures for routine clinical and research application. An updated report from an I.F.C.N. committee.Clin. Neurophysiol.1261071–1107. 10.1016/j.clinph.2015.02.001
27
SollmannN.IlleS.ObermuellerT.NegwerC.RingelF.MeyerB.et al (2015). The impact of repetitive navigated transcranial magnetic stimulation coil positioning and stimulation parameters on human language function.Eur. J. Med. Res.20:47. 10.1186/s40001-015-0138-0
28
StephaniC.PaulusW.SommerM. (2016). The effect of current flow direction on motor hot spot allocation by transcranial magnetic stimulation.Physiol. Rep.4:e12666. 10.14814/phy2.12666
29
SweeneyJ. D.MortimerJ. T.DurandD. (1987). “Modeling of mammalian myelinated nerve for functional neuromuscular electrostimulation,”in Proceedings of the IEEE 97th Annual Conference of the Engineering in Medicine Biology Society,Boston91577–1578.
30
TanakaS.Gomez-TamesJ.WasakaT.InuiK.UenoS.HirataA. (2020). Electrical characterisation of Aδ-Fibres based on human in vivo electrostimulation threshold.Front. Neurosci.14:1305. 10.3389/FNINS.2020.588056
31
TanakaS.SandriniM.CohenL. G. (2011). Modulation of motor learning and memory formation by non-invasive cortical stimulation of the primary motor cortex.Neuropsychol. Rehab.21650–675. 10.1080/09602011.2011.605589
32
TaniK.HirataA.TanakaS. (2020). Quantitative assessment of pain threshold induced by a single-pulse transcranial magnetic stimulation.Front. Neurosci.14:559. 10.3389/fnins.2020.00559
33
TeraoY.UgawaY. (2002). Basic mechanisms of TMS.J. Clin. Neurophysiol. Off. Pub. Am. Electroencephalogr. Soc.19322–343. 10.1097/00004691-200208000-00006
34
YousryT. A.SchmidU. D.AlkadhiH.SchmidtD.PeraudA.BuettnerA.et al (1997). Localization of the motor hand area to a knob on the precentral gyrus. A new landmark.Brain J. Neurol.120(Pt 1)141–157. 10.1093/brain/120.1.141
35
ZiemannU. (2011). Transcranial magnetic stimulation at the interface with other techniques: a powerful tool for studying the human cortex.Neuroscientist17368–381. 10.1177/1073858410390225
Summary
Keywords
perception threshold, pain threshold, dosimetry, nerve model, TMS, standardization, side-effects
Citation
Gomez-Tames J, Tani K, Hayashi K, Tanaka S, Ueno S and Hirata A (2021) Dosimetry Analysis in Non-brain Tissues During TMS Exposure of Broca’s and M1 Areas. Front. Neurosci. 15:644951. doi: 10.3389/fnins.2021.644951
Received
22 December 2020
Accepted
02 February 2021
Published
19 February 2021
Volume
15 - 2021
Edited by
Ravi L. Hadimani, Virginia Commonwealth University, United States
Reviewed by
Shalini Narayana, University of Tennessee Health Science Center (UTHSC), United States; Janakinadh Yanamadala, MathWorks (United States), United States
Updates
Copyright
© 2021 Gomez-Tames, Tani, Hayashi, Tanaka, Ueno and Hirata.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jose Gomez-Tames, jgomez@nitech.ac.jp
This article was submitted to Neural Technology, a section of the journal Frontiers in Neuroscience
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.