- Research article
Transcranial focused ultrasound stimulation of motor cortical areas in freely-moving awake rats
BMC Neurosciencevolume 19, Article number: 57 (2018)
Low-intensity transcranial focused ultrasound (tFUS) has emerged as a new non-invasive modality of brain stimulation with the potential for high spatial selectivity and penetration depth. Anesthesia is typically applied in animal-based tFUS brain stimulation models; however, the type and depth of anesthesia are known to introduce variability in responsiveness to the stimulation. Therefore, the ability to conduct sonication experiments on awake small animals, such as rats, is warranted to avoid confounding effects of anesthesia.
We developed a miniature tFUS headgear, operating at 600 kHz, which can be attached to the skull of Sprague–Dawley rats through an implanted pedestal, allowing the ultrasound to be transcranially delivered to motor cortical areas of unanesthetized freely-moving rats. Video recordings were obtained to monitor physical responses from the rat during acoustic brain stimulation. The stimulation elicited body movements from various areas, such as the tail, limbs, and whiskers. Movement of the head, including chewing behavior, was also observed. When compared to the light ketamine/xylazine and isoflurane anesthetic conditions, the response rate increased while the latency to stimulation decreased in the awake condition. The individual variability in response rates was smaller during the awake condition compared to the anesthetic conditions. Our analysis of latency distribution of responses also suggested possible presence of acoustic startle responses mixed with stimulation-related physical movement. Post-tFUS monitoring of animal behaviors and histological analysis performed on the brain did not reveal any abnormalities after the repeated tFUS sessions.
The wearable miniature tFUS configuration allowed for the stimulation of motor cortical areas in rats and elicited sonication-related movements under both awake and anesthetized conditions. The awake condition yielded diverse physical responses compared to those reported in existing literatures. The ability to conduct an experiment in freely-moving awake animals can be gainfully used to investigate the effects of acoustic neuromodulation free from the confounding effects of anesthesia, thus, may serve as a translational platform to large animals and humans.
Over the past few decades, various brain stimulation techniques have significantly contributed to enhancing our current understanding of neural/neuronal function and offered non-pharmacological options for the treatment of neurological and neuropsychiatry diseases [1,2,3]. Approaches, such as deep brain stimulation (DBS) or epidural cortical stimulation (EpCS) , allow for stimulating brain regions with excellent spatial specificity, but require invasive surgical procedures. Transcranial direct current stimulation (tDCS) and transcranial magnetic stimulation (TMS) provide non-invasive alternatives to the surgical procedures, but may not reach deep brain areas with a centimeter-scale area for stimulation, limiting spatial specificity [1, 2]. Optogenetic techniques are capable of modulating cellular level activity of the brain ; however, the necessary genetic modification of neurons to gain light-sensitivity and limited transcranial penetration of stimulatory light may obstruct its translational application in humans.
Focused ultrasound (FUS) technique allows for the non-invasive, focal delivery of mechanical pressure waves to regional biological tissues [5,6,7], measuring a few millimeters in diameter and length. The advances in FUS techniques have further enabled the transcranial delivery of acoustic energy to specific regions of the brain [8,9,10]. This transcranial FUS (tFUS) technique has been utilized for non-invasive functional neurosurgery by thermally ablating localized deep brain structures, whereby the ultrasound waves are delivered at high acoustic intensities [11, 12]. tFUS has also been applied to temporarily open the blood-brain barrier (BBB) when combined with intravascular administration of microbubbles (detailed review can be found in ). In addition to these therapeutic potentials, tFUS, given in a train of pulses at a low-intensity (under the threshold for heat generation), has been shown to reversibly modulate regional brain excitability [14,15,16,17]. Taking advantage of the exquisite ability to transcranially reach deep brain areas [18, 19] as well as cortical areas [20,21,22,23,24,25] with high spatial selectivity, low-intensity tFUS has rapidly gained momentum as a new mode of non-invasive brain stimulation [26, 27].
FUS has shown to modulate excitability in motor/visual cortical areas in rabbits , stimulated various motor cortices in mice [16, 28,29,30,31,32], suppressed epileptic seizure electroencephalographic (EEG) activities , and altered the extracellular neurotransmitter level [34, 35] and anesthesia time in rats . Investigations have also been conducted to study the effect of varying acoustic parameters  and spatial profile of neuromodulation using a rat model [38, 39]. Additionally, tFUS has stimulated the motor and visual cortices in sheep and elicited corresponding electrophysiological responses . The majority of these studies, conducted on anesthetized animals, showed a degree of variability in response to the stimulation, depending on the types and depths of anesthesia [24, 28, 31, 37, 40]. To examine the behavioral responses to FUS, without the confounding effects from anesthesia, experimentations in an awake setting are desired, and several recent studies on non-human primates and human subjects started to demonstrate the feasibility of tFUS in brain stimulation without the use of anesthesia [18, 20,21,22,23, 25, 41, 42].
We were motivated to develop a technique that will allow tFUS to be applied among unanesthetized, freely-moving small animals. Typically, a FUS transducer, much larger in size than the animal’s head, is maneuvered with optional image/visual-guidance for its stereotactic application during anesthesia [17, 24, 28, 30, 31, 37, 43]. To enable the experimentation in freely-moving small animals, one critical technical element is to make the transducer wearable. Accordingly, we developed a miniaturized, light-weight FUS transducer that can be worn (and detachable) by Sprague–Dawley rats (anesthetized) and demonstrated that FUS can be delivered to their primary somatosensory areas, with possibility for inducing long-term neuromodulatory effects . A 3D-printed applicator that is designed to adjust the position of the transducer was attached to a pedestal, which was implanted onto the rat skull. The design enabled the individual adjustment of location/depth/orientation of the sonication focus. Recently, Li et al.  developed a dual-channel miniature FUS system that can stimulate two separate regions of the mice brain, and observed stimulation-mediated behaviors and extracellular neural action potentials. In their study, the transducers were surgically-fixed to the skull, which granted the use of the system among freely-moving mice. In the present study, we applied our wearable tFUS platform to stimulate motor cortical areas of freely-moving awake rats, and examined sonication-related behavioral responses from three different experimental conditions—(1) freely-moving awake status, (2) ketamine/xylazine anesthesia, and (3) isoflurane anesthesia. The response rates and latencies to the sonication were compared. After the completion of the sonication sessions, histological analysis was conducted on the rat brains to assess the presence of any undesirable tissue damage.
All animal experiments were conducted under the approval of the local Institutional Animal Care and Use Committee.
Preparation of the miniature FUS transducer/headgear
A small (16 mm in diameter, 12 mm in height) and light (~ 6 g in weight) FUS transducer was built in-house (Fig. 1a) . A disc-shape zirconate titanate (PbZrxTi(1−x)O3; PZT) ceramic (American Piezo Ceramics, Mackeyville, PA) was used and fitted (air-backed) inside of a custom-built plastic housing. The plastic housing and back-lid of the transducer was designed (using CAD software; Solidworks Corp., Concord, MA) and printed by three-dimensional (3D) printing (Form2; FormLabs Inc., Somerville, MA). The back-lid of the transducer contained a ball-shape structure to fit the socket of an applicator (also 3D-printed), and held the transducer at a desired location/orientation (Fig. 1a). Both the transducer and applicator constituted the miniature tFUS headgear, and were attached to a pedestal (also 3D-printed), which was implanted on the skull of Sprague–Dawley rat (Charles River Laboratories, Wilmington, MA; see following section). Two set-screws were used to fasten the FUS headgear to the pedestal, ensuring a reproducible placement and orientation via lock-and-key mechanism. To accommodate the differences in individual neuroanatomy and cranial structures, applicators were customized with different ‘Arm’ and ‘Drop’ lengths (Fig. 1a).
Surgical implantation of a pedestal on the rat skull
To apply the miniature tFUS headgear in a wearable form, a pedestal was surgically implanted on the anterior region of the rat’s skull. During the surgery, we measured the relative coordinates between the mounted pedestal and major skull anatomies (i.e., aural meatus, bregma, and lambda) to provide coordinates for the later FUS targeting. Two small screws were inserted (via burr holes) to the skull around the pedestal’s base to provide support along with a medical-grade adhesive (Loctite 18690; Henkel, Rocky Hill, CT). The skin around the pedestal (while exposing the top portion) was sutured back (using Vicryl 5-0 polyglactin 910 suture; Ethicon Inc., Somerville, NJ). After undergoing these surgical procedures, the rats were housed for at least 2 weeks to recover from the surgery prior to the tFUS sessions. The pedestal remained in place and provided long-term mechanical stability over 8 months.
Actuation and characterization of the miniature FUS transducer
A fundamental frequency (FF) of 600 kHz was used to actuate the miniature transducer, and the acoustic intensity profile of the FUS transducer was characterized along the sonication direction as well as on the transversal plane at the focus (Fig. 1b). The detailed methods for the characterization process are described elsewhere . The input signal was a sinusoidal wave generated by a function generator (33210A; Agilent, Santa Clara, CA) and amplified by a class-A linear amplifier (240 L; Electronics and Innovations Ltd., Rochester, NY) with an impedance-matching circuit. At the focus, the miniature transducer was capable of generating over 20 W/cm2 spatial-peak pulse-average intensity (Isppa). The acoustic focus was formed ~ 10 mm away from the exit plane of the transducer. The size of the focus, measured at full-width at half-maximum (FWHM) of acoustic intensity profile, was 11.5 mm in length and 3.0 mm in diameter. When it was measured at full-width at 90%-maximum (FW90%M), previously reported as the spatial dimension of the FUS-mediated neuromodulatory area [38, 39], the focal area was 3.5 mm in length and 1.0 mm in diameter.
Acoustic coupling using PVA gel
A cone-shaped, polyvinyl alcohol (PVA) hydrogel (7–9% weight per volume; two freeze–thaw cycles, U228-08; Avantor, Center Valley, PA) was manufactured in-house for acoustic coupling between the transducer and scalp (Fig. 1c, right) (the detailed method can be found elsewhere ). The hydrogel showed negligible pressure attenuation on the order of 1%. A plastic cone [28, 32] or a bag [37, 39] containing degassed water has been typically used to couple the acoustic path, but could not be used for freely-moving awake animals due to the possibility of water escaping out of the coupling path/container depending on the rat’s dynamic behaviors (such as head-shaking and grooming).
Animal preparation for tFUS sessions
For the tFUS sessions using anesthesia, the Sprague–Dawley rats (all male, n = 7) were anesthetized with either ketamine/xylazine (80:10 mg/kg; intraperitoneal; i.p.) or isoflurane (initial induction with 3–4% followed by 0.5% for the maintenance, at oxygen flow rate of 2 L per min; inhalation). An attempt was made to decrease the maintenance isoflurane concentration under 0.1%, as used by previous investigations in mice [28, 29], but rats emerged from the anesthesia prematurely, and therefore, not used in the present study. The fur on the head was shaved prior to each sonication to prevent any potential blocking of the sonication. The rats were then placed on a custom-built plastic platform in a prone posture with their limbs and tail freely hanging. After positioning the headgear and the accompanying PVA hydrogel, a generic ultrasound gel (Aquasonic; Parker Laboratories, Fairfield, NJ) was applied at each interface. Subsequently, we used the transducer geometry to estimate the virtual focal spot of sonication in space, and aligned the acoustic focus to the motor areas of the tail, limbs, or whiskers (Fig. 1e) while referencing the functional atlas of the rat motor cortex [47, 48]. Once an adequate level of anesthetic plane was detected, such as irregular breathing, the sonication session was conducted. We allowed for slight adjustment in the orientation of the transducer (Fig. 1a) for eliciting motor responses. Also, tFUS was intentionally delivered to off-target locations (lateral or caudal to the target, few millimeters away and including unilateral auditory areas) to examine the spatial specificity in stimulation. After each sonication session, the FUS headgear was removed, and the rats were returned to the housing facility for a minimum of 48 h before the next session (Fig. 1c, left).
To conduct the tFUS experiment in an awake state, we applied the same experimental procedures with the following steps. To shave the fur and apply the tFUS headgear (with the coupling hydrogel), the animals were lightly anesthetized using isoflurane (induction with 3–4%) for ~ 5 min. Then, the rats were moved to an empty cage and allowed to recover until they fully regained their pre-anesthetic behaviors (we determined that ~ 20 min was sufficient across the animals). No additional anesthesia was given to detach the FUS headgear from the pedestal.
Experimental setup compatible with anesthetized/awake rats and data acquisition settings
We established experimental setups that accommodated both anesthetized and awake rats. The schematics of the implemented wearable tFUS headgear, with the transducer actuation systems, are shown in Fig. 1d. A swivel connector (slip ring with flange-736; Adafruit, New York, NY) was located above the middle of the cage/platform, granting unrestricted motion and access to a power source for actuating the transducer during the awake tFUS sessions. A data acquisition system (PowerLab 8/30 and LabChart 7; ADInstruments, Colorado Springs, CO) was used to acquire time-series data of sonication events (onset timing and duration), being synchronized with a video recording (29.97 frames per second; FPS, by QTH44; Q-See; Anaheim, CA) to analyze the location and onset timing of the movement elicited by the sonication. Additionally, a light-emitting diode (LED), turned on in-sync with each sonication event, was placed within the field-of-view of the video recording as a visual indicator of the sonication timing (shown in Fig. 2a–c, upper panels).
Sonication parameters for repeated tFUS sessions with anesthetized/awake rats
We conducted repeated tFUS sessions using a pulsed sonication scheme across all conditions. Based on our previous studies , we used the acoustic parameters (Fig. 1f) as follows: pulse repetition frequency (PRF) of 500 Hz, tone burst duration (TBD) of 1 ms (i.e., a duty cycle of 50%), and sonication duration of 300 ms, with a 5–10 s inter-stimulation interval (ISI), with varying acoustic output (see below). The sonication was administered to the motor areas in the left or right (side randomized) hemisphere of the rat brain. At the initial phase of this study, we gave stimulatory tFUS to each rat brain, starting from an acoustic intensity of 2.1 W/cm2 Isppa, increasing in increments of ~ 1 W/cm2, until the stimulatory response (i.e., movements from the tail, limbs or whiskers) was observed from the ketamine/xylazine as well as awake sessions. We determined that 14.9 W/cm2 Isppa (for ketamine/xylazine anesthesia) and 8.8 W/cm2 Isppa (for awake condition) were most suitable to elicit motor responses (regardless of their type) across all animals. These intensities were used in subsequent measurement of response rates. Acoustic intensity values at the target were estimated based on applying 17% of acoustic pressure attenuation through the rat skull .
Response rates comparison across the repeated different anesthetic/awake conditions
We examined the response rates to the sonication from the same group of animals (n = 7, named as ‘R1’ to ‘R7’) through three repeated tFUS sessions, under each experimental condition. The sequence of these experimental sessions was randomized and balanced. Each tFUS session consisted of a total of 10 sonication events, targeting the tail, limb, or whisker motor areas in the brain. The individual animal’s mean response rates were compared using one-way analysis of variance (ANOVA) within each condition. The grand mean response rates were compared by repeated measures ANOVA and paired t-test across the conditions, with two-sample F-test for the equality of group variances.
Analysis of the FUS-mediated movement location and onset latency
The location of FUS-mediated movement and the onset latency, across all the sonication parameters, were analyzed with high-resolution videos frame-by-frame using video analysis software (Quintic Player v29; Quintic Consultancy Ltd., Sutton Coldfield, UK) by three investigators. The onset of the tFUS was identified from the frame that showed the LED light turned on. A period greater than 500 ms before and after the tFUS onset (i.e., ≥ 15 frames) was examined for each sonication event. Only frames that showed distinctive movements were used to identify the type of movement and to measure the response latency with respect to the tFUS onset. Spontaneous movements from the body (for example, breathing-related movements) or a pattern of whisker movements were excluded to isolate stimulation-specific responses for the analysis.
Examination of potential thermal effect
Potential thermal effect from the sonication was estimated using a formula of ΔT = 2αIt/(ρb∙Cp); where α = the absorption coefficient (0.014 cm−1 at ~ 600 kHz) , I = the intensity of ultrasound in the focal region, t = the ultrasound pulse duration, ρb = the density of brain tissue, and Cp = the specific heat of the brain tissue, where ρb∙Cp is 3.796 J∙cm−3∙ °C [50, 51]. Using the equation, 0.016 °C was the estimated thermal increase, but considering a long ISI (≥ 5 s) (Fig. 1f) and subsequent heat dissipation, in conjunction with the small size of the acoustic focus, this temperature increase was considered to be negligible. An acoustic intensity level that corresponds to the mechanical index (MI) of 1.9, maximum allowed for diagnostic ultrasound device according to the food and drug administration (FDA)-guideline , was 46.5 W/cm2 Isppa at 600 kHz.
Post-sonication behavior monitoring and histological assessment
The biological effects of the repeated sonication sessions were examined across the experimental conditions (awake, ketamine/xylazine, and isoflurane). During the resting and survival periods after the sonication sessions, we regularly monitored the behavior and body condition of the animal for detecting any signs that indicated undesired neurological sequelae, including pain or distress. To examine the potential tissue damage, the animals were sacrificed at short-term (sacrificed within 0.7 ± 1.2 days; n = 3 rats) and long-term (41.5 ± 0.6 days; n = 4 rats) after the end of the last sonication session using the systemic cardiac perfusion of 10% formaldehyde (i.e., the method used to euthanize the animals) under ketamine/xylazine anesthesia, and the fixed brains were harvested. The brains were sectioned along the motor cortical areas, and the presence of hemorrhage, edema, ischemia, gliosis, inflammations were examined through histological analysis. Hematoxylin & eosin (H&E) staining was used to detect cell necrosis or local recruitment of inflammatory cells. Vanadium acid fuchsin (VAF)-toluidine blue staining was used to visualize ischemic neurons. Immunohistochemistry (IHC) of glial fibrillary acidic protein (GFAP) and caspase-3 staining were performed to examine glia infiltration or signs of neurodegeneration and to detect any apoptotic activity at and around the sonicated area, respectively. Two rats belonging to the short-term assessment underwent tail vein injection of the trypan blue dye, within 1 h after the end of the last sonication session to examine the presence of BBB disruption .
Types of elicited responses from anesthetized/awake rats
The average weight of the same group of rats (n = 7, ‘R1–R7’) was 412.7 ± 33.8 g, 395.3 ± 55.0 g, and 388.3 ± 39.6 g (mean ± SD) in the awake, ketamine/xylazine, and isoflurane conditions, respectively (no significant differences, paired t-test, two-tailed, all p > 0.01). Table 1 shows the types of responses elicited by sonication from the wearable tFUS headgear across the conditions. The range of acoustic intensities used for the experiment was 2.3–14.9 W/cm2 Isppa for the awake sessions, 7.5–14.9 W/cm2 Isppa for the ketamine/xylazine sessions, and 9.0–14.9 W/cm2 Isppa for the isoflurane sessions.
The responses were observed above a certain threshold of acoustic intensities, i.e., 3.4 ± 1.8 W/cm2 Isppa (mean ± SD, n = 7) for the awake condition, 10.2 ± 2.4 W/cm2 Isppa (n = 7) for the ketamine/xylazine condition, and 12.4 ± 2.8 W/cm2 Isppa (n = 6) for the isoflurane condition. The acoustic threshold levels from the awake condition were significantly lower than those observed from both anesthetic conditions (t-test, one-tailed, both p < 0.001) while there was no statistical difference between the ketamine/xylazine and isoflurane conditions (t-test, one-tailed, p > 0.05). Also, when tFUS was delivered to off-target locations (including auditory areas) or given under the effective Isppa, no responses were detected.
The elicited movements were seen from either of the tail/limbs/whiskers across all experimental conditons. These movements were similar with previous rodent studies involving ketamine/xylazine anesthesia [16, 31, 37]. We also observed twitches of the head/neck/ears and chewing behaviors in the awake and ketamine/xylazine conditions (listed as ‘other responses’ in Table 1), individually or accompanying the movements from the tail/limbs/whiskers. Under isoflurane anesthesia, the head/neck/ears movements and chewing behaviors were not seen. In terms of their qualitative evaluation, the range of the elicited movement was generally smaller in the case of the awake condition, than those observed from the anesthetic conditions (e.g., video-frame analysis from the tail response; Fig. 2a–c; Additional files 1, 2, 3). The head/neck/ears movements and chewing behaviors in the awake condition can be found in Additional files 4, 5 and 6.
Response rates across the different conditions
The response rate was calculated from each sonication session per each rat (‘R1’–’R7’), and averaged across three sessions. Each animal’s mean response rates (and its standard errors) are shown in Fig. 3 across the three different conditions of (1) awake (Fig. 3a), (2) ketamine/xylazine (Fig. 3b), and (3) isoflurane sessions (Fig. 3c). In the isoflurane condition, one animal (‘R2’) did not show any responses to the sonication.
The mean response rate in the awake sessions ranged 56.7%–86.7% while anesthetic conditions showed much wider ranges (i.e., 36.7%–96.7% in the ketamine/xylazine sessions and 0–96.7% in the isoflurane sessions). To evaluate the inter-animal variability in mean response rate, a one-way ANOVA was performed across the animals, and showed that the mean responses were not significantly different for the awake sessions (p = 0.25). On the other hand, during the anesthetic sessions, the ratio of FUS stimulation events resulted in motor response were significantly different among the animals (one-way ANOVA, p < 0.001 for both ketamine/xylazine and isoflurane conditions). Therefore, the data implicate that response rates were relatively even across the animals during the awake condition compared to those during the anesthetic conditions.
The overall response rate representing each condition was calculated by taking a grand mean of the response rates pooled from all rats (Fig. 3d), and revealed that both awake and ketamine/xylazine conditions showed significantly higher response rates than the isoflurane condition (repeated measures ANOVA, p < 0.05; augmented by paired t-test, one-tailed, p ≤ 0.01 for both awake and ketamine/xylazine sessions compared to the isoflurane sessions). Comparisons of the grand mean response rates between the awake and ketamine/xylazine sessions did not show statistical differences (paired t-test, one-tailed, p = 0.25). Meanwhile, the variability of the grand mean response rate (i.e., variances or dispersions) from the awake condition was significantly decreased compared to those from both anesthetic conditions (two-sample F-test, one-tailed, both p < 0.05), while there was no significant difference between the ketamine/xylazine and isoflurane sessions (p = 0.43).
Onset latency of the elicited movements across the different conditions
The number of events describing the successful tFUS stimulation (resulting in the movement of the tail/limbs/whiskers) and the onset latency were assessed for each condition using a histogram (Fig. 4a–c). Regardless of the experimental conditions, most (> 93%) of these responses were observed within a time frame of ~ 400 ms after the sonication onset. An average latency in motor responses was 139.1 ± 111.1 ms in the awake condition (n = 510), 212.8 ± 127.2 ms under ketamine/xylazine anesthesia (n = 821), and 282.9 ± 103.2 ms under isoflurane anesthesia (n = 293), while these latency values were significantly different to each other across the conditions (one-way ANOVA, p < 0.001; post hoc Tukey test, all p < 0.001). It is notable that the average latency of responses from the tail/limbs/whiskers in the awake condition was shorter than those under the anesthetic conditions.
In the awake and ketamine/xylazine conditions, we observed movements from the head/neck/ears as well as chewing behaviors (Table 1), and the same type of histogram showing its latency distributions was separately constructed (Fig. 4d and e; note that none were detected during the isoflurane sessions). The average latency of 111.9 ± 116.0 ms in the awake condition (n = 592) was also significantly shorter than the latency observed under ketamine/xylazine anesthesia (287.5 ± 178.0 ms; n = 181; t-test, one-tailed, p < 0.001).
To examine the presence of movement that is believed to be associated with acoustic startle responses (ASR) having short latencies (on the order of 10 ms [53,54,55]), we calculated the ratio of responses that occurred within 33 ms after the sonication onset (the limit of the video time frame based on 29.97 FPS), with respect to the total number of observed responses. For the tail/limbs/whiskers movements, the ratio was 32.0% in the awake condition, 14.6% under ketamine/xylazine anesthesia, and 3.8% under isoflurane anesthesia. For the head/neck/ears movements and chewing behaviors, the ratio was 55.1% in the awake condition, and 9.4% under ketamine/xylazine anesthesia. These data demonstrate that a greater portion of the responses occurred at a short latency range (< 33 ms) during the awake sessions.
Post-sonication behavioral monitoring and histological analysis
All animals showed normal behavior and health status after the sonication experiments. The histological analysis (H&E, VAF-toluidine blue, GFAP, and caspase-3 staining) performed on the sonicated brain tissues at a short-term (0.7 ± 1.2 days, n = 3 rats) or long-term (41.5 ± 0.6 days, n = 4 rats) after the last FUS session showed no apparent signs of damage (Fig. 5 shows example slides from rat ‘R6’). The two rats that underwent the tail-vein trypan blue perfusion procedure did not show any signs of BBB disruption.
A miniature FUS transducer was developed in a wearable configuration and transcranially stimulated the motor cortical areas in rats. The transducer unit was attached to an implanted pedestal for each experimental session and detached prior to returning the rats to the animal housing. The location of the acoustic focus was adjusted by the transducer applicator, having different sizes (via 3D-printing) to fit the individual cranial anatomy of the rats. The setup enabled the tFUS experiments to be conducted repeatedly in both awake and anesthetized conditions (either i.p. injection of ketamine/xylazine or isoflurane inhalation). Subsequently, it allowed for systematic condition-specific comparisons of neuromodulatory outcomes, in terms of their physical representations, and response rates/variability with onset latencies. To our knowledge, this is the first study to demonstrate the efficacy of tFUS brain stimulation in awake rats, while having comparisons with two different anesthetic conditions.
Types of elicited responses
The tFUS sonication elicited various physical motor responses across the study. Regardless of the experimental conditons, the elicited movements were seen from either of the tail/limbs/whiskers, demonstrating similiarities with previous rodent studies involving anesthesia [16, 28, 30, 31, 37]. In addition to these FUS-mediated movements, we also observed twitches from the head/neck/ears and chewing behaviors (which are new types of tFUS stimulation-related movement) in the awake and ketamine/xylazine conditions (listed as ‘other responses’ in Table 1). We conjecture that these new-found responses may be associated with the stimulation of corresponding motor areas due to the spatial proximity or overlap with intended motor regions for the whisker and forelimb [47, 48]. For example, imperfections in applying the sonication (e.g., mechanical slippage during application or due to the growth of cranium) can result in slight misalignments of the sonication target. Acoustic reverberation inside a small cavity of the rat skull [40, 56] with the potential to create multiple sonication foci may be another possible cause. It is also plausible that the twitches from the head/neck/ears and chewing behaviors were not seen in the previous studies due to the weight of transducer/coupling devices (water bags or plastic standoffs were used along with much bigger/heavier transducers), which became detectable in the present study using a light-weight wearable tFUS apparatus.
Under isoflurane anesthesia, a previous mice study  did report neck twitching behaviors, however, head/neck/ears movements and chewing behaviors were not seen in the present study. Although the definite causes for this discrepancy is difficult to ascertain, we conjecture that the given anesthetic setting (i.e., 0.5% isoflurane) did not allow sufficient motor neuron recruitment for overt movement. Provision of adequate anesthetic planes, e.g., accomodation of much lower isoflurane concentration using sophisticated anesthetic devices supported by body temperature control [28, 29], will allow for the further exploration of physical responses to tFUS stimulation.
Acoustic intensity to elicit the responses
We found that thresholds existed, in terms of acoustic intensity, in eliciting motor responses. This is congruent with previous studies involving rodents [24, 28, 37] as well as in large animals  and in humans [21, 22]. The threshold acoustic intensity that started to elicit motor responses among the awake rats was much lower than those from anesthetic conditions. This finding is well-aligned with the notion that anesthesia generally suppresses neuronal excitability or dissociate the neural signal connectivity , which may elevate the threshold for excitation. The use of a lower acoustic intensity (in the awake condition), which will reduce overall dosimetry for the sonication, would be particularly advantageous for long and repeated FUS stimulation sessions.
Qualitative examination of the range of the elicited movements
In terms of the qualitative evaluation of the range of the elicited movement, a tail movement, for example, was smaller in the case of the awake condition than those observed from the anesthetic conditions. We speculate that the observation may be attributed to the presence of residual muscle tension during awake state or the animal’s crawling postures that imposed weight to each of the limbs, which may hinder overt motor responses. Further study using measurements of strength of electromyography (EMG) or motor evoked potentials (MEP) is warranted to ascertain the electrophysiological information from FUS-mediated motor responses, especially in freely-moving awake animals.
Response rates and their variability across the different conditions
We found that there were degrees of variability in the response rates among the animals and across the experimental conditions. Existence of such variabilities in the responsiveness were congruent with previous FUS-mediated studies reporting that the types/depths of anesthesia as well as individual differences can alter response rates [24, 28, 31, 37, 40]. Further analysis of inter-animal variability on response rates, measured from the movement data for the tail/limbs/whiskers, showed that the animals during the awake sessions manifested a more consistent level of responses compared to those during the anesthetic conditions. As to the causes for this reduced variability of responses in awake condition, individual-specific responsiveness/susceptibility to the anesthetic agents  as well as the method of its delivery (e.g., i.p. injection of ketamine/xylazine) might have played an important role. Regarding the grand mean response rate, although there were no statistical differences between the awake and ketamine/xylazine sessions, a significant difference did exist for the awake and isoflurane settings. Taken together, the awake condition offers the advantages of higher and more consistent/reproducible response rates compared to the anesthetic conditions.
Onset latency of the elicited movements
Regarding the movement onset latency, most of the elicited responses, either from the tail/limbs/whiskers or from the head/neck/ears and chewing behaviors, were distributed within ~ 400 ms after the onset of the sonication event. An average latency in motor responses (from the tail/limbs/whiskers) was 139.1 ± 111.1 ms for the awake condition, 212.8 ± 127.2 ms for ketamine/xylazine, and 282.9 ± 103.2 ms for isoflurane. We note that the average onset latencies in awake rats were shorter compared to the ones from the anesthetic conditions, which may implicate that the use of anesthesia delays the onset timing of these elicited movements.
In the analysis of onset latency, intriguingly, a greater portion of responses were elicited within ~ 33 ms in the awake condition (over 30% for the tail/limbs/whiskers and over 50% for the head area) compared to below 15% in the anesthetic conditions. These responses having short latencies may be associated with the acoustic startle responses (ASR), known to be occurring within ~ 10 ms after the onset of the acoustic stimuli in rats [53,54,55]. Recently, Sato and colleagues reported a mice study that both ultrasound and audible sound showed similar brain activation patterns and motor response (consistent with a startle reflex) that were reduced by the chemical deafening of the animals , indicating that ultrasound may have an indirect link to acoustic-related (startle) effects and the elicitation of short latency responses. In this perspective, it is not surprising that awake animals, supposedly more susceptible to any external stimuli, showed a higher ratio of responses having short latencies than the anesthetic conditions. Wattiez and colleagues recently reported that cell-level acoustic neuromodulation occurs with an onset latency ≥ ~ 30 ms , lending further support to the idea that responses to the sonication below this latency could be related to startle effects. In the present study, most of the stimulation-related movements were observed at much longer latency, which cannot be explained solely by the ASR. In addition, the stimulation of the auditory areas did not produce any stimulation-related movement. Taken together, our data suggests that one should be aware of the presence of ASR-like phenomena, and exert caution when interpreting the physical responses to the acoustic stimulation.
In reviewing the execution of experimental settings, only the behavioral data was analyzed using video recording due to the lack of measurement of electrophysiological signals, such as EMG. As briefly discussed above, the small range of the elicited movements from awake animals made their detection difficult, which might have possibly contributed to the reduced response rates. These limitations warrant the integration of EMG measurement in future studies using freely-moving awake animals to ascertain the elicitation of the FUS-mediated motor responses. For enabling the EMG measurement from freely-moving awake animals, subdermal wires need to be implanted to the desired body/muscle parts (such as limbs or tail base) , whereby these wires are connected to a multi-channel electrode head pedestal that is compatible with our wearable tFUS headgear. Additional experimental modifications, such as the use of a high-speed camera, could also help to examine the response latencies with a higher time resolution.
We also note that the focal area, 3.5 mm in length and 1.0 mm in diameter measured at FW90%M of its intensity profile, can stimulate the brain regions outside the intended target (the motor cortex), reaching deeper brain structure. Since the present study did not have sufficient spatial resolutions in stimulating discrete rodent functional brain anatomy, the detailed effects of the stimulation on the response rate or the latencies could not be ascertained. We contemplate that use of large animal models (such as ovine, and corresponding larger neuroanatomy) will increase the relative spatial specificity of stimulation compared to that acquired from the rodent model, improving the assessment of region-specific effects of acoustic neuromodulation.
Safety and non-thermal mechanism
In terms of the safety profile, all the animal behaviors were normal, with no brain damage or hemorrhage, after the repeated sonication sessions during a long-term period of ~ 5–8 months. In our previous rat study examining sonication parameters , H&E histology on a rat’s brain exposed to 22.4 W/cm2 Isppa (corresponding to a spatial-peak temporal-average intensity of 11.2 W/cm2 Ispta with peak rarefactional pressure of 0.81 MPa, MI of 1.38) showed hemosiderin indicating potential earlier bleeding, while such signs were not observed in the present study with 14.9 W/cm2 Isppa (7.5 W/cm2 Ispta, 0.67 MPa, MI of 0.86). We conjectured that the use of longer ISIs (≥ 5 s vs. previously 2 s) and lower MI, with a miniature tFUS transducer having a smaller acoustic focus, compared to those used in the previous studies, possibly prevented the occurrence of sonication related brain hemorrhage. Also, the estimated potential thermal increase of 0.016 °C (see Methods), which is believed to be negligible considering heat dissipation during the ISI (≥ 5 s) and the small size of acoustic focus, supports that the biophysical mechanism behind the tFUS stimulation of neural cells could be linked with non-thermal mechanical factors . The present work utilized the sonication parameters that are compliant with safety guidelines for the diagnostic ultrasound equipment (with an exception of the maximum MI of 0.23 for ophthalmological applications). However, we note that there is neither clear consensus nor the data on the sonication parameters (such as the acoustic intensity and the MI) for safe brain tissue stimulation. Further studies are, therefore, urgently needed to establish the safety guidelines for the acoustic neuromodulation.
We demonstrated the application of FUS brain stimulation in a freely-moving rat model, utilizing a wearable tFUS headgear. The awake rats showed an increased response rate with reduced variability and shorter latency to FUS, in comparison with the neurostimulatory outcomes under the anesthetic conditions. Our analysis of latency distribution of responses suggests the possible involvement of ASR-like phenomena mixed with the stimulation-related physical movement. The use of small animal models, without confounding factors from anesthesia (including its unclear mechanism of action ), would be beneficial not only to gain further knowledge for reducing the variability (thus, may increase the reproducibility) in responsiveness to FUS but to gain more informative data regarding the potential presence of ASR. The ability to conduct FUS-mediated brain stimulation in awake small animals provides unprecedented opportunities for investigations that are not possible with anesthesia, such as sociobehavioral studies (e.g., self-administered brain stimulation ), or for the studies dealing with disease models that are influenced by anesthesia (e.g., epilepsy ).
deep brain stimulation
epidural cortical stimulation
transcranial direct current stimulation
transcranial magnetic stimulation
transcranial focused ultrasound
full-width at half-maximum
full-width at 90%-maximum
frames per second
pulse repetition frequency
tone burst duration
analysis of variance
food and drug administration
hematoxylin & eosin
vanadium acid fuchsin
glial fibrillary acidic protein
acoustic startle responses
motor evoked potentials
Fregni F, Pascual-Leone A. Technology insight: noninvasive brain stimulation in neurology-perspectives on the therapeutic potential of rTMS and tDCS. Nat Clin Pract Neurol. 2007;3(7):383–93.
George MS, Aston-Jones G. Noninvasive techniques for probing neurocircuitry and treating illness: vagus nerve stimulation (VNS), transcranial magnetic stimulation (TMS) and transcranial direct current stimulation (tDCS). Neuropsychopharmacology. 2010;35(1):301–16.
Hoy KE, Fitzgerald PB. Brain stimulation in psychiatry and its effects on cognition. Nat Rev Neurol. 2010;6(5):267–75.
Miesenböck G. The optogenetic catechism. Science. 2009;326(5951):395–9.
Fry WJ, Barnard JW, Fry FJ, Brennan JF. Ultrasonically produced localized selective lesions in the central nervous system. Am J Phys Med. 1955;34(3):413–23.
Jolesz FA, Hynynen K, McDannold N, Tempany C. MR imaging-controlled focused ultrasound ablation: a noninvasive image-guided surgery. Magn Reson Imaging Clin N Am. 2005;13(3):545–60.
Lynn JG, Zwemer RL, Chick AJ, Miller AE. A new method for the generation and use of focused ultrasound in experimental biology. J Gen Physiol. 1942;26(2):179–93.
Hynynen K, Clement GT, McDannold N, Vykhodtseva N, King R, White PJ, Vitek S, Jolesz FA. 500-element ultrasound phased array system for noninvasive focal surgery of the brain: a preliminary rabbit study with ex vivo human skulls. Magn Reson Med. 2004;52(1):100–7.
Pinton G, Aubry J-F, Bossy E, Muller M, Pernot M, Tanter M. Attenuation, scattering, and absorption of ultrasound in the skull bone. Med Phys. 2012;39(1):299–307.
White PJ, Clement GT, Hynynen K. Longitudinal and shear mode ultrasound propagation in human skull bone. Ultrasound Med Biol. 2006;32(7):1085–96.
Elias WJ, Huss D, Voss T, Loomba J, Khaled M, Zadicario E, Frysinger RC, Sperling SA, Wylie S, Monteith SJ, et al. A pilot study of focused ultrasound thalamotomy for essential tremor. N Engl J Med. 2013;369(7):640–8.
Martin E, Jeanmonod D, Morel A, Zadicario E, Werner B. High-intensity focused ultrasound for noninvasive functional neurosurgery. Ann Neurol. 2009;66(6):858–61.
Cammalleri A, Croce P, Lee W, Yoon K, Yoo S-S. Therapeutic potentials of localized blood-brain barrier disruption by non-invasive transcranial focused ultrasound: A technical review. J Clin Neurophysiol. 2018; (in press).
Bachtold MR, Rinaldi PC, Jones JP, Reines F, Price LR. Focused ultrasound modifications of neural circuit activity in a mammalian brain. Ultrasound Med Biol. 1998;24(4):557–65.
Rinaldi PC, Jones JP, Reines F, Price LR. Modification by focused ultrasound pulses of electrically evoked responses from an in vitro hippocampal preparation. Brain Res. 1991;558(1):36–42.
Tufail Y, Matyushov A, Baldwin N, Tauchmann ML, Georges J, Yoshihiro A, Tillery SIH, Tyler WJ. Transcranial pulsed ultrasound stimulates intact brain circuits. Neuron. 2010;66(5):681–94.
Yoo S-S, Bystritsky A, Lee J-H, Zhang Y, Fischer K, Min B-K, McDannold NJ, Pascual-Leone A, Jolesz FA. Focused ultrasound modulates region-specific brain activity. Neuroimage. 2011;56(3):1267–75.
Legon W, Ai L, Bansal P, Mueller JK. Neuromodulation with single-element transcranial focused ultrasound in human thalamus. Hum Brain Mapp. 2018;39(5):1995–2006.
Monti MM, Schnakers C, Korb AS, Bystritsky A, Vespa PM. Non-invasive ultrasonic thalamic stimulation in disorders of consciousness after severe brain injury: a first-in-man report. Brain Stimul. 2016;9(6):940–1.
Lee W, Chung YA, Jung Y, Song I-U, Yoo S-S. Simultaneous acoustic stimulation of human primary and secondary somatosensory cortices using transcranial focused ultrasound. BMC Neurosci. 2016;17(1):68.
Lee W, Kim H, Jung Y, Song I-U, Chung YA, Yoo S-S. Image-guided transcranial focused ultrasound stimulates human primary somatosensory cortex. Sci Rep. 2015;5:8743.
Lee W, Kim H-C, Jung Y, Chung YA, Song I-U, Lee J-H, Yoo S-S. Transcranial focused ultrasound stimulation of human primary visual cortex. Sci Rep. 2016;6:34026.
Lee W, Kim S, Kim B, Lee C, Chung YA, Kim L, Yoo S-S. Non-invasive transmission of sensorimotor information in humans using an EEG/focused ultrasound brain-to-brain interface. PLoS ONE. 2017;12(6):e0178476.
Lee W, Lee SD, Park MY, Foley L, Purcell-Estabrook E, Kim H, Fischer K, Maeng L-S, Yoo S-S. Image-guided focused ultrasound-mediated regional brain stimulation in sheep. Ultrasound Med Biol. 2016;42(2):459–70.
Legon W, Sato TF, Opitz A, Mueller J, Barbour A, Williams A, Tyler WJ. Transcranial focused ultrasound modulates the activity of primary somatosensory cortex in humans. Nat Neurosci. 2014;17(2):322–9.
Bystritsky A, Korb AS, Douglas PK, Cohen MS, Melega WP, Mulgaonkar AP, DeSalles A, Min B-K, Yoo S-S. A review of low-intensity focused ultrasound pulsation. Brain Stimul. 2011;4(3):125–36.
Yoo S-S, Lee W, Jolesz FA. Chapter 23. FUS-mediated image-guided neuromodulation of the brain. In: Chen Y, Kateb B, editors. Neurophotonics and brain mapping. Boca Raton: CRC Press; 2017. p. 443–55.
King RL, Brown JR, Newsome WT, Pauly KB. Effective parameters for ultrasound-induced in vivo neurostimulation. Ultrasound Med Biol. 2013;39(2):312–31.
King RL, Brown JR, Pauly KB. Localization of ultrasound-induced in vivo neurostimulation in the mouse model. Ultrasound Med Biol. 2014;40(7):1512–22.
Li G-F, Zhao H-X, Zhou H, Yan F, Wang J-Y, Xu C-X, Wang C-Z, Niu L-L, Meng L, Wu S. Improved anatomical specificity of non-invasive neuro-stimulation by high frequency (5 MHz) ultrasound. Sci Rep. 2016;6:24738.
Mehić E, Xu JM, Caler CJ, Coulson NK, Moritz CT, Mourad PD. Increased anatomical specificity of neuromodulation via modulated focused ultrasound. PLoS ONE. 2014;9(2):e86939.
Ye PP, Brown JR, Pauly KB. Frequency dependence of ultrasound neurostimulation in the mouse brain. Ultrasound Med Biol. 2016;42(7):1512–30.
Min B-K, Bystritsky A, Jung K-I, Fischer K, Zhang Y, Maeng L-S, Park SI, Chung Y-A, Jolesz F, Yoo S-S. Focused ultrasound-mediated suppression of chemically-induced acute epileptic EEG activity. BMC Neurosci. 2011;12:23.
Min B-K, Yang PS, Bohlke M, Park S, Vago DR, Maher TJ, Yoo S-S. Focused ultrasound modulates the level of cortical neurotransmitters: potential as a new functional brain mapping technique. Int J Imaging Syst Technol. 2011;21(2):232–40.
Yang PS, Kim H, Lee W, Bohlke M, Park S, Maher TJ, Yoo S-S. Transcranial focused ultrasound to the thalamus is associated with reduced extracellular GABA levels in rats. Neuropsychobiology. 2012;65(3):153–60.
Yoo S-S, Kim H, Min B-K, Franck E, Park S. Transcranial focused ultrasound to the thalamus alters anesthesia time in rats. NeuroReport. 2011;22(15):783–7.
Kim H, Chiu A, Lee SD, Fischer K, Yoo S-S. Focused ultrasound-mediated non-invasive brain stimulation: examination of sonication parameters. Brain Stimul. 2014;7(5):748–56.
Kim H, Lee SD, Chiu A, Yoo S-S, Park S. Estimation of the spatial profile of neuromodulation and the temporal latency in motor responses induced by focused ultrasound brain stimulation. NeuroReport. 2014;25(7):475–9.
Kim H, Park M-A, Wang S, Chiu A, Fischer K, Yoo S-S. PET∕CT imaging evidence of FUS-mediated (18)F-FDG uptake changes in rat brain. Med Phys. 2013;40(3):033501.
Younan Y, Deffieux T, Larrat B, Fink M, Tanter M, Aubry J-F. Influence of the pressure field distribution in transcranial ultrasonic neurostimulation. Med Phys. 2013;40(8):082902.
Deffieux T, Younan Y, Wattiez N, Tanter M, Pouget P, Aubry J-F. Low-intensity focused ultrasound modulates monkey visuomotor behavior. Curr Biol. 2013;23(23):2430–3.
Wattiez N, Constans C, Deffieux T, Daye PM, Tanter M, Aubry J-F, Pouget P. Transcranial ultrasonic stimulation modulates single-neuron discharge in macaques performing an antisaccade task. Brain Stimul. 2017;10(6):1024–31.
Kim H, Chiu A, Park S, Yoo S-S. Image-guided navigation of single-element focused ultrasound transducer. Int J Imaging Syst Technol. 2012;22(3):177–84.
Yoo S-S, Yoon K, Croce P, Cammalleri A, Margolin RW, Lee W. Focused ultrasound brain stimulation to anesthetized rats induces long-term changes in somatosensory evoked potentials. Int J Imaging Syst Technol. 2018;28(2):106–12.
Li G, Qiu W, Zhang Z, Jiang Q, Su M, Cai R, Li Y, Cai F, Deng Z, Xu D et al. Noninvasive ultrasonic neuromodulation in freely moving mice. IEEE Trans Biomed Eng. 2018; (in press).
Lee W, Lee SD, Park MY, Yang J, Yoo S-S. Evaluation of polyvinyl alcohol cryogel as an acoustic coupling medium for low-intensity transcranial focused ultrasound. Int J Imaging Syst Technol. 2014;24(4):332–8.
Fonoff ET, Pereira JF, Camargo LV, Dale CS, Pagano RL, Ballester G, Teixeira MJ. Functional mapping of the motor cortex of the rat using transdural electrical stimulation. Behav Brain Res. 2009;202(1):138–41.
Tennant KA, Adkins DL, Donlan NA, Asay AL, Thomas N, Kleim JA, Jones TA. The organization of the forelimb representation of the C57BL/6 mouse motor cortex as defined by intracortical microstimulation and cytoarchitecture. Cereb Cortex. 2011;21(4):865–76.
Goss SA, Frizzell LA, Dunn F. Ultrasonic absorption and attenuation in mammalian tissues. Ultrasound Med Biol. 1979;5(2):181–6.
Elwassif MM, Kong Q, Vazquez M, Bikson M. Bio-heat transfer model of deep brain stimulation-induced temperature changes. J Neural Eng. 2006;3(4):306–15.
O’Brien WD. Ultrasound—biophysics mechanisms. Prog Biophys Mol Biol. 2007;93(1–3):212–55.
Duck FA. Medical and non-medical protection standards for ultrasound and infrasound. Prog Biophys Mol Biol. 2007;93(1–3):176–91.
Gómez-Nieto R, Horta-Júnior JDAC, Castellano O, Millian-Morell L, Rubio ME, López DE. Origin and function of short-latency inputs to the neural substrates underlying the acoustic startle reflex. Front Neurosci. 2014;8:216.
Pilz PK, Schnitzler HU. Habituation and sensitization of the acoustic startle response in rats: amplitude, threshold, and latency measures. Neurobiol Learn Mem. 1996;66(1):67–79.
Yeomans JS, Li L, Scott BW, Frankland PW. Tactile, acoustic and vestibular systems sum to elicit the startle reflex. Neurosci Biobehav Rev. 2002;26(1):1–11.
Yoon K, Lee W, Croce P, Cammalleri A, Yoo S-S. Multi-resolution simulation of focused ultrasound propagation through ovine skull from a single-element transducer. Phys Med Biol. 2018;63(10):105001.
Århem P, Klement G, Nilsson J. Mechanisms of anesthesia: towards integrating network, cellular, and molecular level modeling. Neuropsychopharmacology. 2003;28(S1):S40–7.
Sato T, Shapiro MG, Tsao DY. Ultrasonic neuromodulation causes widespread cortical activation via an indirect auditory mechanism. Neuron. 2018;98(5):1031–41.
Towne C, Montgomery KL, Iyer SM, Deisseroth K, Delp SL. Optogenetic control of targeted peripheral axons in freely moving animals. PLoS ONE. 2013;8(8):e72691.
Kubanek J, Shukla P, Das A, Baccus SA, Goodman MB. Ultrasound elicits behavioral responses through mechanical effects on neurons and ion channels in a simple nervous system. J Neurosci. 2018;38(12):3081–91.
Olds J. Pleasure centers in the brain. Sci Am. 1956;195(4):105–17.
WL, PC, RWM, AC, KY, SSY participated in study design, data acquisition and analysis. WL, PC, AC, RWM performed animal surgical procedures. WL, PC, RWM, SSY prepared equipment settings. All participated in manuscript writing. All authors read and approved the final manuscript.
We thank Dr. Yongzhi Zhang for the helpful technical advice on rat surgical procedures. The initial help by Mr. Michael Y. Park is also acknowledged. We thank the Center for Comparative Medicine (CCM) of the Brigham and Women’s Hospital (BWH) for the trainings on handling anesthetized and freely-moving awake animals.
The authors declare that they have no competing interests.
Availability of data and materials
All data generated or analyzed during this study are included in this published article and its supplementary information files. The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Consent for publication
Ethics approval and consent to participate
All animal experiments were conducted under the approval of the local Institutional Animal Care and Use Committee (IACUC) at the Brigham and Women’s Hospital (Protocol number: 2016N000495).
This study was supported by the Focused Ultrasound Surgery Foundation (FUS 461, to W. Lee and S.S. Yoo) and the National Institutes of Health (NIH, R01 MH111763, to S.S. Yoo).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.