Skip to main content

ORIGINAL RESEARCH article

Front. Mol. Neurosci., 04 June 2020
Sec. Neuroplasticity and Development
Volume 13 - 2020 | https://doi.org/10.3389/fnmol.2020.00070

Brain Metabolite Changes After Anodal Transcranial Direct Current Stimulation in Autism Spectrum Disorder

  • 1Department of Pediatrics, Faculty of Medicine, Khon Kaen University, Khon Kaen, Thailand
  • 2Department of Psychiatry, Faculty of Medicine, Khon Kaen University, Khon Kaen, Thailand
  • 3Department of Radiology, Faculty of Medicine, Khon Kaen University, Khon Kaen, Thailand
  • 4Department of Physiology, Faculty of Medicine, Khon Kaen University, Khon Kaen, Thailand
  • 5Department of Anesthesiology, Faculty of Medicine, Khon Kaen University, Khon Kaen, Thailand
  • 6Integrated Epilepsy Research Group, Faculty of Medicine, Khon Kaen University, Khon Kaen, Thailand
  • 7Faculty of Physical Therapy, Mahidol University, Salaya, Thailand

Objectives: Previous research has provided evidence that transcranial direct current stimulation (tDCS) can reduce severity of autism spectrum disorder (ASD); however, the exact mechanism of this effect is still unknown. Magnetic resonance spectroscopy has demonstrated low levels of brain metabolites in the anterior cingulate cortex (ACC), amygdala, and left dorsolateral prefrontal cortex (DLPFC) in individuals with ASD. The aim of this study was to investigate the effects of anodal tDCS on social functioning of individuals with ASD, as measured by the social subscale of the Autism Treatment Evaluation Checklist (ATEC), through correlations between pretreatment and posttreatment concentrations of brain metabolites in the areas of interest (DLPFC, ACC, amygdala, and locus coeruleus) and scores on the ATEC social subscale.

Methods: Ten participants with ASD were administered 1 mA anodal tDCS to the left DLPFC for 20 min over five consecutive days. Measures of the ATEC social subscale and the concentrations of brain metabolites were performed before and immediately after the treatment.

Results: The results showed a significant decrease between pretreatment and immediately posttreatment in the ATEC social subscale scores, significant increases in N-acetylaspartate (NAA)/creatine (Cr) and myoinositol (mI)/Cr concentrations, and a decrease in choline (Cho)/Cr concentrations in the left DLPFC and locus coeruleus after tDCS treatment. Significant associations between decreased ATEC social subscale scores and changed concentrations in NAA/Cr, Cho/Cr, and mI/Cr in the locus coeruleus were positive.

Conclusion: Findings suggest that beneficial effects of tDCS in ASD may be due to changes in neuronal and glia cell activity and synaptogenesis in the brain network of individuals with ASD. Further studies with larger sample sizes and control groups are warranted.

Introduction

The abnormalities of brain functions in individuals with autism spectrum disorder (ASD) are not fully understood (Trottier et al., 1999). There is some evidence of a specific local increase in the thickness of the fusiform gyrus associated with face processing impairment in individuals with ASD (Schultz, 2005; Dziobek et al., 2010), as well as evidence of abnormal synaptic maturation, resulting in neurobiological connectivity defects (Levy et al., 2009). Several studies have proposed that decreased cortical plasticity may play an important role in ASD pathogenesis (Markram et al., 2007; Oberman et al., 2012, 2016). At present, there is no effective treatment for ASD, although behavioral therapy is a standard recommendation (Myers et al., 2007); however, the outcomes of this therapy are unsatisfactory. In severe cases with attention-deficit disorder, pharmacological therapies may cause adverse effects such as drowsiness, dry mouth, nausea, agitation, behavioral activation, and sleep problems (Oswald and Sonenklar, 2007).

Transcranial direct current stimulation (tDCS) is one of the non-invasive brain stimulation techniques that was shown to be effective in reducing severity of ASD (Schneider and Hopp, 2011; Amatachaya et al., 2014, Amatachaya et al., 2015; D’Urso et al., 2015; Gómez et al., 2017; English et al., 2018; Esse Wilson et al., 2018). In a previous tDCS study with individuals with ASD, low-voltage stimulation (1–2 mA) was applied via anodal electrode to the left dorsolateral prefrontal cortex (DLPFC) for 20–30 min, and a positive effect was found immediately after the treatment (Schneider and Hopp, 2011). Likewise, our previous study showed significant improvement on the Autism Treatment Evaluation Checklist (ATEC): total scores decreased from 67.25 pretreatment to 58 posttreatment (p < 0.001). We also found a significant decrease in the mean ATEC social subscale scores in the tDCS group (14.45) compared to the sham group (17.7) (p = 0.015) (Amatachaya et al., 2014, Amatachaya et al., 2015). A study in which cathodal 1-mA tDCS over the left DLPFC was applied, a significant decrease in symptoms of ASD lasting 1 month was found (Gómez et al., 2017). D’Urso et al. (2015) used 10 sessions of cathodal tDCS over the left DLPFC and found 26.7% reduction of the total Aberrant Behavior Checklist score in ASD. Taken together, this evidence suggests that tDCS could be effective for individuals with ASD. However, there were no clinical studies that investigate the mechanism of action of tDCS in reducing severity of ASD symptom.

Magnetic resonance spectroscopy (MRS) is a non-invasive technique for measuring biochemical changes in different organs including the brain. It is able to quantify steady-state metabolic levels of neurotransmitters such as N-acetylaspartate (NAA), which is a marker of axonal integrity; the small molecules such as choline (Cho) in the cell membrane composition, which is a marker for membrane turnover and reflects neuronal connection (Cudalbu et al., 2012); glutamine combined glutamate (Glx), an important excitatory neurotransmitter; myoinositol (mI), an osmolyte and astrocyte marker, which reflects synaptogenesis (Walecki et al., 2003; Kubas et al., 2012); and creatine (Cr), which is often used as an internal standard that normalizes other metabolites (Cudalbu et al., 2012). Previous studies that used MRS provided evidence of abnormal brain metabolites in individuals with ASD (Levitt et al., 2003; DeVito et al., 2007; Endo et al., 2007; Kleinhans et al., 2007; Bernardi et al., 2011; Bejjani et al., 2012; Mori et al., 2013).

The DLPFC is the brain area for which positive findings after anodal tDCS were reported. For example, Fujii et al. (2010) found that NAA/Cr was significantly decreased, whereas Cho/Cr was significantly increased in the left DLPFC in pediatric participants with ASD compared with control group participants (Fujii et al., 2010). Horder et al. (2014) reported lower concentrations of NAA and no significant differences in concentrations of Cho and Glx in the left DLPFC of adult participants with ASD compared with control group participants (Horder et al., 2014). They suggested that NAA/Cr deficits in the brain region responsible for executive functions may be associated with social and communication difficulties in ASD (Fujii et al., 2010; Shepherd and Freiwald, 2018). There were also reports showing low levels of Cho in the left inferior anterior cingulate cortex (ACC) (Levitt et al., 2003), low levels of Glx in the right ACC (Bernardi et al., 2011), and decreased NAA, Cho, and mI in the ACC(Goji et al., 2017). Other MRS studies of ASD reported low levels of NAA in the left amygdala (Mori et al., 2013) and low levels of NAA and Glx in the left frontal cortex (Kleinhans et al., 2007). Recently, there was a hypothesis that explored the possibility of explaining the symptoms found in ASD in terms of inefficient neuromodulation using the functioning of the locus coeruleus and norepinephrine as exemplars (London, 2018).

These evidences showed that there were many brain areas involved in determining the ASD core symptoms such as DLPFC, ACC, amygdala, and locus coeruleus. Magnetic resonance spectroscopy is a safe method that provides in-depth information about the abnormal brain metabolites in such brain areas in ASD. To the best of our knowledge, there is no study about the relationship between these metabolite changes after tDCS, so every type of brain metabolites might provide a preliminary investigation of biological changed following tDCS. We hypothesize that tDCS can modulate the brain function in ASD, and then the brain metabolite concentrations would change after tDCS treatment and the symptoms of ASD would decrease. Therefore, the primary objective of this study was to find out concentrations of brain metabolite changes in the areas of interest after anodal tDCS. The secondary objective was to study the associations between brain metabolite changes and the ATEC social subscale score changes.

Materials and Methods

Study Design

This study was a pilot study composed of three phases: (1) baseline evaluation consisted of the baseline characteristic assessments, vital signs, ATEC, Childhood Autism Rating Scale (CARS), and baseline MRS; (2) treatment period consisted of 1 mA anodal tDCS for 20 min over five consecutive days and the post–tDCS treatment MRS, which was immediately performed after the fifth tDCS session; (3) the 2 weeks’ follow-up performed at days 7 and 14 after tDCS. The ATEC, CARS, adverse events of tDCS, and MRS were assessed at the first and second weeks after tDCS. Participants were asked to continue their routine medication regimen throughout the duration of the study. The schematic representation of timeline of the study is presented in Figure 1.

FIGURE 1
www.frontiersin.org

Figure 1. Schematic representation of experimental timeline.

Participant Recruitment and Informed Consent

Study participants were recruited by advertisement at the Pediatric Outpatient Children’s Psychiatric Clinic; Child Neurology Clinic; Child Development Clinic of Srinagarind Hospital, Faculty of Medicine, Khon Kaen University; and Khon Kaen Special Education Center Region 9, Thailand. The study procedures were described to caregivers who participating in the study by clinic physicians. Autism spectrum disorder diagnosis was confirmed by a child psychiatrist following a clinical review of the Diagnostic and Statistical Manual of Mental Disorders, Fifth Edition criteria (American Psychiatric Association, 2013).

The inclusion criteria were the same as in our previous study: (a) male participants with ASD; (b) age between 5 and 8 years; and (c) severe symptoms of ASD (CARS score ≥37). The exclusion criteria were the following: (a) having a pacemaker or metallic device; (b) severe neurological disorders such as brain tumor or intracranial infection; (c) drug abuse; (d) epilepsy; (e) skull defect; (f) use of herbal remedies or other alternative therapies; and (g) uncooperative parents and caregivers.

A routine brain magnetic resonance imaging (MRI) was performed on all participants to confirm the absence of organic disease, and MRS measurements were performed immediately afterward. Anesthesia was induced with 2 mg/kg intravenous propofol. After a laryngeal mask airway was inserted, participants received 2% sevoflurane and a 40:60 air/oxygen mixture for anesthesia maintenance at a total gas flow rate of 2 L/min. All participants were managed by the anesthetist team throughout the sedation period, and the procedures were performed following the guidelines for monitoring and management of pediatric patient sedation published by the American Academy of Pediatrics (Coté and Wilson, 2016).

The study was conducted in accordance with the Declaration of Helsinki and was approved by the Ethics Committee of Khon Kaen University (identifier no. HE 561188). The written informed consents were obtained from all participants and caregivers before participation.

Measures

Two main outcomes were assessed in this study, the primary outcome was the ATEC social subscale, and the secondary outcome was the levels of brain metabolites in the brain areas of interest such as DLPFC, amygdala, ACC, and locus coeruleus. We also assessed vital signs, physical examination, and neurological examination to evaluate the possibility that tDCS might have adverse effects.

Autism Treatment Evaluation Checklist

The ATEC is a questionnaire used to evaluate the effectiveness of treatments for individuals with ASD, which is completed by caregivers. It consists of four subscales: (1) speech/language/communication subscale (14 items; ceiling score 28); (2) social subscale (20 items; ceiling score 40); (3) sensory and cognitive awareness subscale (18 items; ceiling score 36); and (4) health/physical/behavior subscale (25 items; ceiling score 75). The total score ranges from 0 to 179; a higher score indicates worsening, whereas a lower score indicates improvement (Geier et al., 2012). The participants were assessed by the ATEC at baseline, immediately postsession on day 5, week 1, and week 2 after tDCS treatment.

Since there was research showing that NAA/Cr levels in the left DLPFC were related to social function (Fujii et al., 2010; Shepherd and Freiwald, 2018), we aimed to study associations between the change in the ATEC social subscale and brain metabolite posttreatment.

Childhood Autism Rating Scale

The CARS is a well-established measure of ASD severity (Schopler et al., 1980; Chlebowski et al., 2010). The participants were evaluated using the CARS conducted by three investigators (N.P., C.S., and P.A.) who observed the participants and interviewed the parents and were unaware as to the treatment status of the participants. The details of the CARS evaluation were provided in as our previous study (Amatachaya et al., 2014, Amatachaya et al., 2015). The CARS is a 15-item behavioral rating scale developed to identify ASD and quantitatively describe its severity. The 15 items are as follows: relating to people, imitative behavior, emotional response, body use, object use, adaptation to change, visual response, listening response, other sensory responses, fear or anxiety, verbal communication, non-verbal communication, activity level, level and consistency of intellectual response, and general impressions (Rellini et al., 2004). Assessment was performed at baseline and immediately postsession on day 5, week 1, and week 2 after tDCS treatment.

Brain Metabolite Levels

Our areas of interest were both DLPFC, both amygdalae, both ACC, and locus coeruleus. We used MRS to assess baseline and posttreatment levels of NAA, Glx, Cho, mI, and Cr in the regions of interest. Single-voxel proton magnetic resonance spectra were obtained from the areas of interest using a Philips Achieva 3.0 T (Philips Healthcare, Best, Netherlands) running Release 2.6.3.3 MR Workspace software. The procedures were performed under propofol-induced sedation to ensure that participants would lie still for approximately 60 min while the data were collected. Single-voxel proton magnetic resonance spectra were acquired and quantified with LC Model to determine metabolite concentration ratios. Magnetic resonance spectroscopy voxels (2 cm × 2 cm × 2 cm) were positioned on the coronal, sagittal, and axial images of the DLPFC, amygdala, ACC, and locus coeruleus. The spectra used a point resolved spectroscopy (PRESS) sequentially with a repetition time of 2 s, short echo time of 35 ms, spectral width of 2,000 Hz, 1,024 time points, and partial water suppression. Shimming was performed using manufacturer-supplied shimming procedures. The analysis of the metabolite concentrations was performed using LC-Model (Stephen Provencher Inc., Oakville, ON, Canada). Levels of NAA, Glx, mI, Cr, and Cho were analyzed by fitting a linear combination of a basis set of metabolite model spectra to the data. The analyzing spectrum was set between 3.8 and 0.2 ppm with no eddy-current correction and water scaling. The metabolite concentrations were expressed as nm and ratios relative to Cr peak. The metabolite concentrations and metabolite-to-Cr ratios were determined in NAA, Cho, mI, and Glx spectra for each subject.

Adverse Events and Safety

Although the electric current induced by tDCS is weak (1 mA well below the pain threshold), it was applied continuously for five sessions of 20 min each (over the course of 5 days). The participants were asked to report any adverse events as well as other signs and symptoms immediately after each stimulation session. The participants were also closely observed by physicians during the study sessions. They were also followed up for adverse effects 7 and 14 days after the treatment.

Transcranial Direct Current Stimulation

Transcranial direct current stimulation was applied via 0.9% NaCl-soaked pair of surface sponge electrodes (35 cm2) and delivered through battery-driven power supply. The constant direct current stimulator had a maximum output of 10 mA (Model 1224-B; Soterix Medical, New York, NY, United States). The anode electrode was placed over the left DLPFC, which was located using the 10–20 international system of electrode placement. A 1-mA current was applied for 20 min once a day for five consecutive days during the treatment period. The cathode, or the reference electrode, was placed on the right shoulder contralateral to the anode. The current was applied gradually by increasing it until the necessary current level was reached and then decreasing it after the stimulation.

Data Analysis

Analyses were performed using the Stata software, version 10.0 (StataCorp, College Station, TX, United States). Data are presented as mean and standard deviation. Because of small sample size, the normality test was considered, and we found that most data were non-normal distribution. Therefore, we used the non-parametric Wilcoxon signed-ranks test to evaluate differences between pretreatment and posttreatment. Pre and post social subscale scores of ATEC and CARS were calculated in percent change. We examined the associations between the differences in brain metabolite levels and the social subscale of the ATEC by computing Pearson correlation coefficients between pretreatment to posttreatment changes in the measures of these domains. The effect sizes were determined by the mean difference between pretreatment and posttreatment and then dividing the result by the pooled standard deviation (Cohen d). Data are presented as means and SD. p < 0.05 was considered significant.

Results

Ten children with ASD were screened for possible participation between November 2014 and February 2016, and all met the study inclusion criteria. Six right-handed and four left-handed participants completed the entire protocol without any adverse events. The mean age of the participants was 6.60 years (SD = 0.84 years); the mean age at diagnosis was 29.1 months (SD = 6.03 months). The demographic data are presented in Table 1.

TABLE 1
www.frontiersin.org

Table 1. Demographic data of participants (n = 10).

Autism Treatment Evaluation Checklist Scores

The changes on the ATEC subscale scores from baseline to immediately posttreatment revealed a decrease of social subscale scores from 12.5 to 10.7 (14.4% change, p = 0.023); health subscale scores from 12.0 to 10.5 (12.5% change, p = 0.010), language subscale scores from 8.1 to 7.4 (8.6% change, p = 0.102), and sensory and cognitive function subscales from 23.3 to 21.3 (8.6% change, p = 0.016), as well as in the total ATEC score from 55.9 to 45.9 (17.9% change, p = 0.011). Data of ATEC social subscale scores at immediately posttreatment, week 1, and week 2 after treatments are presented in Table 2.

TABLE 2
www.frontiersin.org

Table 2. Change of ASD symptoms at immediately posttreatment, week 1, and week 2 after treatment (n = 10).

Childhood Autism Rating Scale Scores

At day 5 of anodal tDCS, the average absolute difference of CARS score was 2.7 (SD = 0.75, SEM = 0.24). The average percent change of CARS score was 6.14% (3.5–8.75%) A significant decrease in symptoms of ASD from baseline to immediately posttreatment was found (p = 0.005). With respect to the effect size (Cohen d), the pretreatment to immediately posttreatment change in CARS score was 0.84 (clinical significance). Data of CARS scores at immediately posttreatment, week 1, and week 2 after treatment are presented in Table 2.

Changes in Brain Metabolites

Data of metabolites concentration at the baseline, immediately posttreatment in the left DLPFC, right DLPFC, left ACC, right ACC, left amygdala, right amygdala, and locus coeruleus are shown in Table 3. To study the changes in brain metabolites in the seven areas of interest, in the left DLPFC, Wilcoxon signed-ranks tests revealed a significant increase in NAA/Cr (p = 0.016) and mI/Cr (p = 0.012) and significant decrease in Cho/Cr (p = 0.037) (Figure 3). The effect sizes associated with the changes in NAA/Cr, Cho/Cr, and mI/Cr were 0.86, 0.82, and 0.82, respectively (clinically significant). The effect size associated with the change in Glx/Cr was 0.10 (negative). In the locus coeruleus, Wilcoxon signed-ranks tests revealed significant increase in NAA/Cr (p = 0.005) and mI/Cr (p = 0.025), decrease in Cho/Cr (p = 0.003), but no significant change in Glx/Cr (p = 0.120) at the locus coeruleus, from baseline to posttreatment (Figure 4). The effect sizes associated with the changes in NAA/Cr, Cho/Cr, mI/Cr, and Glx/Cr were 1.6, 0.69, 0.81, and 0.59, respectively.

TABLE 3
www.frontiersin.org

Table 3. Metabolites concentration in the areas of interest at baseline and immediately posttreatment (n = 10).

FIGURE 2
www.frontiersin.org

Figure 2. The voxels were placed at (A) left DLPFC, (B) right DLPFC, (C) right ACC, (D) left ACC, (E) right amygdala, (F) left amygdala, and (G) locus coeruleus. The yellow box indicates the location of a single voxel in the areas of interest. Abbreviations: A, anterior; F, feet; H, head; L, left; P, posterior; R, right.

FIGURE 3
www.frontiersin.org

Figure 3. Metabolite changes in the left dorsolateral prefrontal cortex. Data are presented as means of NAA/Cr, Cho/Cr, mI/Cr, and Glx/Cr concentrations, compared between baseline (before tDCS) and after a 5-day tDCS treatment (day 5). Vertical lines represent SD. Represents a significant difference; p < 0.05.

FIGURE 4
www.frontiersin.org

Figure 4. Metabolite changes in the locus coeruleus. Data are presented as means of NAA/Cr, Cho/Cr, mI/Cr, and Glx/Cr concentrations, compared between baseline (before tDCS) and after a 5-day tDCS treatment (day 5). Vertical lines represent SD. Represents significant difference; ∗∗p < 0.01, p < 0.05.

The Associations Between the ATEC Social Subscale and Metabolite Changes

To study the associations between the social subscale of the ATEC and metabolite changes in the seven areas of interest, we computed a series of Pearson correlation coefficients. Our data showed significant associations between increased NAA/Cr and decreased ATEC social subscale (r = −0.684; p = 0.015), decreased Cho/Cr and decreased ATEC social subscale (r = 0.780; p = 0.004), and increased mI/Cr and decreased ATEC social subscale (r = −0.872; p = 0.001) in the locus coeruleus. The associations between the ATEC social subscale and metabolite changes in other areas of interest are shown in Table 4.

TABLE 4
www.frontiersin.org

Table 4. Association between metabolite changes in the areas of interest and decreased ATEC social subscale scores (n = 10).

Discussion

Our study confirmed that tDCS had 17.9% reduction of ASD core symptoms assessed by ATEC total scores. The preliminary investigation of biological change following tDCS showed positive effect in the left DLPFC and the locus coeruleus. The significant associations between decreased ATEC social subscale and NAA/Cr, Cho/Cr, and mI/Cr concentration changes in the locus coeruleus were positive. N-acetylaspartate/Cr was tentatively increased in the right DLPFC. However, there were no significant changes in other brain metabolites either in the right DLPFC or in the left and right ACC and amygdala after the tDCS treatment.

Our study also showed clinical significant decrease in CARS scores from 44.0 to 41.30 (6.14% change), corresponding to the improvement of social–emotional understanding, emotional response, visual response, and verbal communication. The outcomes were similar to our previous study, which also revealed a significant decrease in CARS score from pretreatment to posttreatment. The effect of tDCS was maintained for 7 days in participants with the active tDCS condition relative to those in the sham tDCS condition (Amatachaya et al., 2014). Moreover, there was a report that showed a higher rate of chance corrected agreement (sensitivity of 0.89) between the CARS and clinical judgment for ASD (Chlebowski et al., 2010). The mechanisms of improved clinical outcomes of ASD after tDCS are not fully understood (Schneider and Hopp, 2011; Amatachaya et al., 2014, Amatachaya et al., 2015; D’Urso et al., 2015; Gómez et al., 2017). In this study, we surveyed the seven brain areas that have evidence related to ASD. We found changes in brain function of two parts, for example, the left DLPFC and the locus coeruleus detected by MRS, in which the brain metabolite changes were in the same pattern.

N-acetylaspartate is an amino acid derivative found exclusively in the central nervous system (CNS) neurons. It is considered as a marker of global neuronal health and attenuation (Verma et al., 2016). Lower levels of NAA are found in a variety of CNS disorders that cause neuronal cell destruction (Gonen et al., 2000; Rigotti et al., 2007) including ASD (Fujii et al., 2010). In general, NAA level quantified with 1H-MRS is thought to reflect neural density and indicate neuronal integrity and metabolism (Ford and Crewther, 2016). The decrease in parietal axon density and marked decrease in white matter NAA in children with ASD are associated with deficits in social functioning and memory (Levitt et al., 2003). Similarly, reduced NAA/Cr in the ACC is associated with poor social functioning (Fujii et al., 2010).

The results of this study support those of the previous study that NAA/Cr was lower in children with ASD (Fujii et al., 2010). The exact mechanism of increase in NAA after tDCS is not yet known. However, there is a possibility that the anodal tDCS can generally augment the neuronal excitability, and the current induces a sustainable response in the form of a long-term potentiation such as plasticity (Auvichayapat and Auvichayapat, 2011; Pelletier and Cicchetti, 2015).

Measured cortical Cho levels by 1H-MRS indicate membrane phosphatidylcholine breakdown. Choline levels in children are generally reduced, particularly in cortical gray matter, temporal regions, and the left thalamus, suggesting a decrease in membrane phospholipid turnover. In children with ASD, Cho/Cr level is reported to be increased. It represents abnormal regional increased membrane phospholipid turnover (Vasconcelos et al., 2008). Our result revealed a decrease in Cho/Cr after anodal tDCS treatment. This finding may suggest that tDCS could decrease regional membrane phospholipid breakdown (Pelletier and Cicchetti, 2015).

Myoinositol is thought to be a marker of astrocytes and plays an important role in glial cell proliferation. It also maintains cell metabolism and signaling as an intracellular postreceptor second messenger system. This second messenger system is linked to several receptors including glutamate receptors in the CNS. At present, there is no report of mI/Cr level in children with ASD (Ford and Crewther, 2016). However, our study showed an increase in mI/Cr as well as a positive association between the increased mI/Cr and decreased ATEC social subscale score after tDCS over the left DLPFC. It may be postulated that tDCS could reduce severity of ASD by increased synaptogenesis at the stimulation site (presynaptic neuron) or increased neuronal signaling to other brain areas (postsynaptic neuron) (Pelletier and Cicchetti, 2015).

There is still a lack of evidence about Glx concentration in individuals with ASD. Glutamate is involved in neurotransmitter regulation and detoxification. Reduced Glx in the ACC may relate to executive function deficits such as decision-making, impulse control, empathy, and emotion. However, our study showed no significant change in Glx/Cr after tDCS. This is compatible with a study that showed no effect of anodal tDCS on Glx/Cr levels in the left posterior superior temporal gyrus (Dwyer et al., 2018). This negative finding could be postulated that N-methyl-D-aspartate receptor and enzyme synthesis, which modulate the availability of neuronal, glial, and synaptic glutamate and glutamine, go undetected in 1H-MRS, thus inherently affecting the level of detected Glx (Ford and Crewther, 2016), or the Glx concentration did not change after anodal stimulation (Stagg et al., 2009).

Changes in brain function in the area of interest that were far away from the stimulation site should be explained by two concepts. First is the putative molecular mechanism of tDCS. Direct current electrical field (DCEF) can guide neural cell migration, growth cone direction, differentiation, and metabolism as a phenomenon called electrotaxis (Pelletier and Cicchetti, 2015). It is suspicious why a change has been found at the locus coeruleus, the most far away from the stimulation site. This can be explained by the greater total number of neurite branches found at the area near the cathode (McCaig et al., 2005). In this study, our cathode electrode was placed at the right shoulder, in which the locus coeruleus was the most nearby. Other areas of interest were far away from the cathode than the locus coeruleus, so the putative molecular mechanism of tDCS might occur difficultly and found the negative result for the brain metabolite changes. Second, there was more connectivity of the left DLPFC to locus coeruleus than the left DLPFC to other areas of interest.

Much of the recent research on ASD has focused not only on the brain structures but also on neuronal connectivity, and both underconnectivity and overconnectivity were reported based on connection lengths (London, 2018). The apparent abnormal neuronal connectivity in individuals with autistic is not fully understood. In general, the core symptoms of ASD are impaired social communication and restricted and repetitive behaviors. Neuromodulatory changes in ASD are suspected of involving many brain regions such as the frontal, temporal, and temporoparietal lobe, insular cortex, ACC, locus coeruleus, and limbic system (London, 2018; Shepherd and Freiwald, 2018). Input from a small number of neuromodulatory cells can abruptly interrupt the activity of neural networks and reorganize the elements into new functional networks; a single neuron can participate in several networks, and a single anatomical network can mediate multiple functions (London, 2018).

After stimulation at the left DLPFC, the results of our experiment showed a significant change in brain metabolites in the locus coeruleus.

The locus coeruleus is a brain area that receives inputs from a number of other brain regions including the prefrontal cortex; this connection is constant and excitatory and increases in strength with raised activity level of the individual (Ford and Crewther, 2016). The change in brain metabolites in the locus coeruleus in our study could be explained by the mechanism found in the DLPFC. The changes in the locus coeruleus may be a result of the changes in presynaptic neurons in the DLPFC. However, the associations between decreased scores on the ATEC social subscale and metabolite changes in the locus coeruleus were more positive than those in the left DLPFC. However, at present, there is little evidence of the locus coeruleus’ role in ASD. Nevertheless, there is evidence that atypical noradrenergic activity in the locus coeruleus is the underlying mechanism of aberrant attentional function in children with ASD (Bast et al., 2018). Thus, it is possible that the function of the locus coeruleus may be important in determining severity of ASD.

Taken together, changes in brain metabolites in the locus coeruleus after anodal tDCS can imply that the stimulation over the left DLPFC not only causes neuromodulation in the stimulated area but also modulates the locus coeruleus, the long-distance connectivity, which is impaired in individuals with ASD. The increased NAA/Cr and mI/Cr and decreased Cho/Cr in the locus coeruleus can reflect the mechanism of neuromodulation by increased neuronal integration, connectivity (Cudalbu et al., 2012), and synaptogenesis (Walecki et al., 2003; Kubas et al., 2012).

Limitations

The present study has some limitations that should be considered when interpreting the results. First, as a pilot study, it had a small sample size that has been unable to detect significant effects in all areas of interest. Therefore, the findings presented here should be viewed as tentative and in need for replication by studies with larger samples. Second, it lacks a control or sham group. Further study using a control group would also be useful.

Summary and Conclusion

To the best of our knowledge, this is the first study that addresses the beneficial effects of anodal tDCS on amelioration of ASD core symptoms. The study provides a preliminary investigation of biological changes following tDCS and the statistical changes of brain metabolites found in the left DLPFC and locus coeruleus, while there were no brain metabolite changes for the right DLPFC, left amygdala, right amygdala, left ACC, and right ACC. Increased NAA/Cr represents increase in neuronal function, increased mI/Cr represents increased glial cell proliferation and synaptogenesis, and decreased Cho/Cr represents the decreased cell membrane breakdown. The mechanism by which anodal tDCS improves the symptoms of ASD may related to the changes of these brain metabolites. However, a study with a larger sample size and a control group is warranted.

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 the Ethics Committee of Khon Kaen University (Identifier number: HE 561188). Written informed consent to participate in this study was provided by the participants’ legal guardian/next of kin.

Author Contributions

PA, NA, and ST contributed conception and design of the study. BA, WK, and WP organized the database. OT and PA performed the statistical analysis. PA and KT wrote the first draft of the manuscript. PA, NP, and CS assessed autistic severity and abnormal behaviors. WP assessed brain MRI and MRS. WB was an MRS technician. AS was an anesthesiologist. PA, KT, and NA wrote sections of the manuscript. All authors contributed to manuscript revision, read and approved the submitted version.

Funding

This work was supported by a grant of the National Research Council of Thailand (NRCT), a grant of the Research and Service Institute for Autism, and a grant of Integrated Epilepsy Research Group.

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.

Acknowledgments

The authors thank Dr. Anuwat Amatachaya, Assist Prof Dr. Niran Ngernyam for participant recruitment.

References

Amatachaya, A., Auvichayapat, N., Patjanasoontorn, N., Suphakunpinyo, C., Ngernyam, N., Aree-Uea, B., et al. (2014). Effect of anodal transcranial direct current stimulation on autism: a randomized double-blind crossover trial. Behav. Neurol. 2014:173073. doi: 10.1155/2014/173073

PubMed Abstract | CrossRef Full Text | Google Scholar

Amatachaya, A., Jensen, M. P., Patjanasoontorn, N., Auvichayapat, N., Suphakunpinyo, C., Janjarasjitt, S., et al. (2015). The short-term effects of transcranial direct current stimulation on electroencephalography in children with autism: a randomized crossover controlled trial. Behav. Neurol. 2015:928631. doi: 10.1155/2015/928631

PubMed Abstract | CrossRef Full Text | Google Scholar

American Psychiatric Association (2013). DSM-5 Task Force. Diagnostic and Statistical Manual of Mental Disorders: DSM-5. Available online at: https://books.google.co.th/books?hl=en&lr=&id=-JivBAAAQBAJ&oi=fnd&pg=PT18&dq=American+Psychiatric+Association:+Diagnostic+and+Statistical+Manual+of+Mental+Disorders+2013+5th+ed+Arlington+American+Psychiatric+Association+&ots=ceTN14KEs9&sig=-OOTSotIoKSNaOQBquaRMaOYPyo&redir_esc=y#v=onepage&q=American Psychiatric Association%3A Diagnostic and Statistical Manual of Mental Disorders 2013 5th ed Arlington American Psychiatric Association&f=false (accessed July 16, 2019).

Google Scholar

Auvichayapat, P., and Auvichayapat, N. (2011). Basic knowledge of transcranial direct current stimulation. J. Med. Assoc. Thail. 94, 518–527.

Google Scholar

Bast, N., Poustka, L., and Freitag, C. M. (2018). The locus coeruleus-norepinephrine system as pacemaker of attention - a developmental mechanism of derailed attentional function in autism spectrum disorder. Eur. J. Neurosci. 47, 115–125. doi: 10.1111/ejn.13795

PubMed Abstract | CrossRef Full Text | Google Scholar

Bejjani, A., O’Neill, J., Kim, J. A., Frew, A. J., Yee, V. W., Ly, R., et al. (2012). Elevated glutamatergic compounds in pregenual anterior cingulate in pediatric autism spectrum disorder demonstrated by 1H MRS and 1H MRSI. PLoS One 7:e38786. doi: 10.1371/journal.pone.0038786

PubMed Abstract | CrossRef Full Text | Google Scholar

Bernardi, S., Anagnostou, E., Shen, J., Kolevzon, A., Buxbaum, J. D., Hollander, E., et al. (2011). In vivo 1H-magnetic resonance spectroscopy study of the attentional networks in autism. Brain Res. 1380, 198–205. doi: 10.1016/j.brainres.2010.12.057

PubMed Abstract | CrossRef Full Text | Google Scholar

Chlebowski, C., Green, J. A., Barton, M. L., and Fein, D. (2010). Using the childhood autism rating scale to diagnose autism spectrum disorders. J. Autism Dev. Disord. 40, 787–799. doi: 10.1007/s10803-009-0926-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Coté, C. J., and Wilson, S. (2016). Guidelines for monitoring and management of pediatric patients before, during, and after sedation for diagnostic and therapeutic procedures: update 2016. Pediatrics 138:e20161212. doi: 10.1542/peds.2016-1212

PubMed Abstract | CrossRef Full Text | Google Scholar

Cudalbu, C., Lanz, B., Duarte, J. M., Morgenthaler, F. D., Pilloud, Y., Mlynárik, V., et al. (2012). Cerebral glutamine metabolism under hyperammonemia determined in vivo by localized 1 H and 15 N NMR spectroscopy. J. Cereb. Blood Flow Metab. 32, 696–708. doi: 10.1038/jcbfm.2011.173

PubMed Abstract | CrossRef Full Text | Google Scholar

DeVito, T. J., Drost, D. J., Neufeld, R. W. J., Rajakumar, N., Pavlosky, W., Williamson, P., et al. (2007). Evidence for cortical dysfunction in autism: a proton magnetic resonance spectroscopic imaging study. Biol. Psychiatry 61, 465–473. doi: 10.1016/j.biopsych.2006.07.022

PubMed Abstract | CrossRef Full Text | Google Scholar

D’Urso, G., Bruzzese, D., Ferrucci, R., Priori, A., Pascotto, A., Galderisi, S., et al. (2015). Transcranial direct current stimulation for hyperactivity and noncompliance in autistic disorder. World J. Biol. Psychiatry 16, 361–366. doi: 10.3109/15622975.2015.1014411

PubMed Abstract | CrossRef Full Text | Google Scholar

Dwyer, G. E., Craven, A. R., Hirnstein, M., Kompus, K., Assmus, J., Ersland, L., et al. (2018). No effects of anodal tDCS on local GABA and Glx levels in the left posterior superior temporal gyrus. Front. Neurol. 9:1145. doi: 10.3389/fneur.2018.01145

PubMed Abstract | CrossRef Full Text | Google Scholar

Dziobek, I., Bahnemann, M., Convit, A., and Heekeren, H. R. (2010). The role of the fusiform-amygdala system in the pathophysiology of autism. Arch. Gen. Psychiatry 67, 397–405. doi: 10.1001/archgenpsychiatry.2010.31

PubMed Abstract | CrossRef Full Text | Google Scholar

Endo, T., Shioiri, T., Kitamura, H., Kimura, T., Endo, S., Masuzawa, N., et al. (2007). Altered chemical metabolites in the amygdala-hippocampus region contribute to autistic symptoms of autism spectrum disorders. Biol. Psychiatry 62, 1030–1037. doi: 10.1016/j.biopsych.2007.05.015

PubMed Abstract | CrossRef Full Text | Google Scholar

English, M. C. W., Kitching, E. S., Maybery, M. T., and Visser, T. A. W. (2018). Modulating attentional biases of adults with autistic traits using transcranial direct current stimulation: a pilot study. Autism Res. 11, 385–390. doi: 10.1002/aur.1895

PubMed Abstract | CrossRef Full Text | Google Scholar

Esse Wilson, J., Trumbo, M. C., Wilson, J. K., and Tesche, C. D. (2018). Transcranial direct current stimulation (tDCS) over right temporoparietal junction (rTPJ) for social cognition and social skills in adults with autism spectrum disorder (ASD). J. Neural Transm. 125, 1857–1866. doi: 10.1007/s00702-018-1938-1935

PubMed Abstract | CrossRef Full Text | Google Scholar

Ford, T. C., and Crewther, D. P. (2016). A comprehensive review of the 1H-MRS metabolite spectrum in autism spectrum disorder. Front. Mol. Neurosci. 9:14. doi: 10.3389/fnmol.2016.00014

PubMed Abstract | CrossRef Full Text | Google Scholar

Fujii, E., Mori, K., Miyazaki, M., Hashimoto, T., Harada, M., and Kagami, S. (2010). Function of the frontal lobe in autistic individuals: a proton magnetic resonance spectroscopic study. J. Med. Invest. 57, 35–44. doi: 10.2152/jmi.57.35

PubMed Abstract | CrossRef Full Text | Google Scholar

Geier, D. A., Kern, J. K., and Geier, M. R. (2012). A prospective cross-sectional cohort assessment of health, physical, and behavioral problems in autism spectrum disorders. Maedica 7, 193–200.

PubMed Abstract | Google Scholar

Goji, A., Ito, H., Mori, K., Harada, M., Hisaoka, S., Toda, Y., et al. (2017). Assessment of anterior cingulate cortex (ACC) and left cerebellar metabolism in asperger’s syndrome with proton magnetic resonance spectroscopy (MRS). PLoS One 12:e0169288. doi: 10.1371/journal.pone.0169288

PubMed Abstract | CrossRef Full Text | Google Scholar

Gómez, L., Vidal, B., Maragoto, C., Morales, L. M., Berrillo, S., Vera Cuesta, H., et al. (2017). Non-invasive brain stimulation for children with autism spectrum disorders: a short-term outcome study. Behav. Sci. 7:63. doi: 10.3390/bs7030063

PubMed Abstract | CrossRef Full Text | Google Scholar

Gonen, O., Catalaa, I., Babb, J. S., Ge, Y., Mannon, L. J., Kolson, D. L., et al. (2000). Total brain N-acetylaspartate: a new measure of disease load in MS. Neurology 54, 15–19. doi: 10.1212/WNL.54.1.15

PubMed Abstract | CrossRef Full Text | Google Scholar

Horder, J., Lavender, T., Mendez, M. A., O’Gorman, R., Daly, E., Craig, M. C., et al. (2014). Reduced subcortical glutamate/glutamine in adults with autism spectrum disorders: a [(1)H]MRS study. Transl. Psychiatry 4:e364. doi: 10.1038/tp.2014.7

PubMed Abstract | CrossRef Full Text | Google Scholar

Kleinhans, N. M., Schweinsburg, B. C., Cohen, D. N., Müller, R.-A., and Courchesne, E. (2007). N-acetyl aspartate in autism spectrum disorders: regional effects and relationship to fMRI activation. Brain Res. 1162, 85–97. doi: 10.1016/j.brainres.2007.04.081

PubMed Abstract | CrossRef Full Text | Google Scholar

Kubas, B., Kułak, W., Sobaniec, W., Tarasow, E., Łebkowska, U., and Walecki, J. (2012). Metabolite alterations in autistic children: a 1H MR spectroscopy study. Adv. Med. Sci. 57, 152–156. doi: 10.2478/v10039-012-0014-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Levitt, J. G., O’Neill, J., Blanton, R. E., Smalley, S., Fadale, D., McCracken, J. T., et al. (2003). Proton magnetic resonance spectroscopic imaging of the brain in childhood autism. Biol. Psychiatry 54, 1355–1366. doi: 10.1016/S0006-3223(03)00688-687

PubMed Abstract | CrossRef Full Text | Google Scholar

Levy, S. E., Mandell, D. S., and Schultz, R. T. (2009). Autism. Lancet 374, 1627–1638. doi: 10.1016/S0140-6736(09)61376-3

PubMed Abstract | CrossRef Full Text | Google Scholar

London, E. B. (2018). Neuromodulation and a reconceptualization of autism spectrum disorders: using the locus coeruleus functioning as an exemplar. Front. Neurol. 9:1120. doi: 10.3389/fneur.2018.01120

PubMed Abstract | CrossRef Full Text | Google Scholar

Markram, H., Rinaldi, T., and Markram, K. (2007). The intense world syndrome–an alternative hypothesis for autism. Front. Neurosci. 1, 77–96. doi: 10.3389/neuro.01.1.1.006.2007

PubMed Abstract | CrossRef Full Text | Google Scholar

McCaig, C. D., Rajnicek, A. M., Song, B., and Zhao, M. (2005). Controlling cell behavior electrically: current views and future potential. Physiol. Rev. 85, 943–978. doi: 10.1152/physrev.00020.2004

PubMed Abstract | CrossRef Full Text | Google Scholar

Mori, K., Toda, Y., Ito, H., Mori, T., Goji, A., Fujii, E., et al. (2013). A proton magnetic resonance spectroscopic study in autism spectrum disorders: amygdala and orbito-frontal cortex. Brain Dev. 35, 139–145. doi: 10.1016/j.braindev.2012.09.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Myers, S. M., and Johnson, C. P., American Academy of Pediatrics Council on Children With Disabilities (2007). Management of children with autism spectrum disorders. Pediatrics 120, 1162–1182. doi: 10.1542/peds.2007-2362

PubMed Abstract | CrossRef Full Text | Google Scholar

Oberman, L., Eldaief, M., Fecteau, S., Ifert-Miller, F., Tormos, J. M., and Pascual-Leone, A. (2012). Abnormal modulation of corticospinal excitability in adults with Asperger’s syndrome. Eur. J. Neurosci. 36, 2782–2788. doi: 10.1111/j.1460-9568.2012.08172.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Oberman, L. M., Ifert-Miller, F., Najib, U., Bashir, S., Heydrich, J. G., Picker, J., et al. (2016). Abnormal mechanisms of plasticity and metaplasticity in autism spectrum disorders and fragile X syndrome. J. Child Adolesc. Psychopharmacol. 26, 617–624. doi: 10.1089/cap.2015.0166

PubMed Abstract | CrossRef Full Text | Google Scholar

Oswald, D. P., and Sonenklar, N. A. (2007). Medication use among children with autism spectrum disorders. J. Child Adolesc. Psychopharmacol. 17, 348–355. doi: 10.1089/cap.2006.17303

PubMed Abstract | CrossRef Full Text | Google Scholar

Pelletier, S. J., and Cicchetti, F. (2015). Cellular and molecular mechanisms of action of transcranial direct current stimulation: evidence from in vitro and in vivo models. Int. J. Neuropsychopharmacol. 18:pyu047. doi: 10.1093/ijnp/pyu047

PubMed Abstract | CrossRef Full Text | Google Scholar

Rellini, E., Tortolani, D., Trillo, S., Carbone, S., and Montecchi, F. (2004). Childhood autism rating scale (CARS) and autism behavior checklist (ABC) correspondence and conflicts with DSM-IV criteria in diagnosis of autism. J. Autism Dev. Disord. 34, 703–708. doi: 10.1007/s10803-004-5290-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Rigotti, D. J., Inglese, M., and Gonen, O. (2007). Whole-Brain N-Acetylaspartate as a surrogate marker of neuronal damage in diffuse neurologic disorders. Am. J. Neuroradiol. 28, 1843–1849. doi: 10.3174/ajnr.A0774

PubMed Abstract | CrossRef Full Text | Google Scholar

Schneider, H. D., and Hopp, J. P. (2011). The use of the bilingual aphasia test for assessment and transcranial direct current stimulation to modulate language acquisition in minimally verbal children with autism. Clin. Linguist. Phon. 25, 640–654. doi: 10.3109/02699206.2011.570852

PubMed Abstract | CrossRef Full Text | Google Scholar

Schopler, E., Reichler, R. J., DeVellis, R. F., and Daly, K. (1980). Toward objective classification of childhood autism: childhood autism rating scale (CARS). J. Autism Dev. Disord. 10, 91–103. doi: 10.1007/BF02408436

PubMed Abstract | CrossRef Full Text | Google Scholar

Schultz, R. T. (2005). Developmental deficits in social perception in autism: the role of the amygdala and fusiform face area. Int. J. Dev. Neurosci. 23, 125–141. doi: 10.1016/J.IJDEVNEU.2004.12.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Shepherd, S. V., and Freiwald, W. A. (2018). Functional networks for social communication in the macaque monkey. Neuron 99, 413–420.e3. doi: 10.1016/j.neuron.2018.06.027

PubMed Abstract | CrossRef Full Text | Google Scholar

Stagg, C. J., Best, J. G., Stephenson, M. C., O’Shea, J., Wylezinska, M., Kineses, Z. T., et al. (2009). Polarity-sensitive modulation of cortical neurotransmitters by transcranial stimulation. J. Neurosci. 29, 5202–5206. doi: 10.1523/JNEUROSCI.4432-08.2009

PubMed Abstract | CrossRef Full Text | Google Scholar

Trottier, G., Srivastava, L., and Walker, C. D. (1999). Etiology of infantile autism: a review of recent advances in genetic and neurobiological research. J. Psychiatry Neurosci. 24, 103–115.

PubMed Abstract | Google Scholar

Vasconcelos, M. M., Brito, A. R., Domingues, R. C., da Cruz, L. C. H., Gasparetto, E. L., Werner, J., et al. (2008). Proton magnetic resonance spectroscopy in school-aged autistic children. J. Neuroimaging 18, 288–295. doi: 10.1111/j.1552-6569.2007.00200.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Verma, A., Kumar, I., Verma, N., Aggarwal, P., and Ojha, R. (2016). Magnetic resonance spectroscopy - Revisiting the biochemical and molecular milieu of brain tumors. BBA Clin. 5, 170–178. doi: 10.1016/j.bbacli.2016.04.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Walecki, J., Tarasów, E., Kubas, B., Czemicki, Z., Lewko, J., Podgórski, J., et al. (2003). Hydrogen-1 MR spectroscopy of the peritumoral zone in patients with cerebral glioma: assessment of the value of the method. Acad. Radiol. 10, 145–153.

PubMed Abstract | Google Scholar

Keywords: autism spectrum disorder, transcranial direct current stimulation, brain metabolites, magnetic resonance spectroscopy, locus coeruleus

Citation: Auvichayapat N, Patjanasoontorn N, Phuttharak W, Suphakunpinyo C, Keeratitanont K, Tunkamnerdthai O, Aneksan B, Klomjai W, Boonphongsathian W, Sinkueakunkit A, Punjaruk W, Tiamkao S and Auvichayapat P (2020) Brain Metabolite Changes After Anodal Transcranial Direct Current Stimulation in Autism Spectrum Disorder. Front. Mol. Neurosci. 13:70. doi: 10.3389/fnmol.2020.00070

Received: 30 December 2019; Accepted: 08 April 2020;
Published: 04 June 2020.

Edited by:

Claudio Grassi, Catholic University of the Sacred Heart, Italy

Reviewed by:

Franca Tecchio, National Research Council (Cnr), Italy
Maria Vittoria Podda, Catholic University of the Sacred Heart, Italy

Copyright © 2020 Auvichayapat, Patjanasoontorn, Phuttharak, Suphakunpinyo, Keeratitanont, Tunkamnerdthai, Aneksan, Klomjai, Boonphongsathian, Sinkueakunkit, Punjaruk, Tiamkao and Auvichayapat. 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: Paradee Auvichayapat, aparad@kku.ac.th

Members of Non-invasive Brain Stimulation Research Group of Thailand

Download