ORIGINAL RESEARCH article
Sec. Medical Physics and Imaging
Volume 11 - 2020 | https://doi.org/10.3389/fphys.2020.01042
Neuronavigated Repetitive Transcranial Ultrasound Stimulation Induces Long-Lasting and Reversible Effects on Oculomotor Performance in Non-human Primates
- 1Institute of Brain and Spinal Cord, UMRS 975 Inserm, CNRS 7225, UMPC, Paris, France
- 2Rogue Research Inc., Montreal, QC, Canada
- 3Physics for Medicine Paris, Inserm, ESPCI Paris, CNRS, PSL Research University, Univ Paris Diderot, Sorbonne Paris Cite, Paris, France
- 4Université de Paris, Institute of Psychiatry and Neuroscience of Paris (IPNP), Inserm U1266, Team Pathophysiology of Psychiatric Disorders, Paris, France
- 5GHU Paris Psychiatrie et Neurosciences, Site Sainte-Anne, Service Hospitalo-Universitaire, Paris, France
- 6Institut de la Vision CNRS, Inserm, Sorbonne Université, Paris, France
Since the late 2010s, Transcranial Ultrasound Stimulation (TUS) has been used experimentally to carryout safe, non-invasive stimulation of the brain with better spatial resolution than Transcranial Magnetic Stimulation (TMS). This innovative stimulation method has emerged as a novel and valuable device for studying brain function in humans and animals. In particular, single pulses of TUS directed to oculomotor regions have been shown to modulate visuomotor behavior of non-human primates during 100 ms ultrasound pulses. In the present study, a sustained effect was induced by applying 20-s trains of neuronavigated repetitive Transcranial Ultrasound Stimulation (rTUS) to oculomotor regions of the frontal cortex in three non-human primates performing an antisaccade task. With the help of MRI imaging and a frame-less stereotactic neuronavigation system (SNS), we were able to demonstrate that neuronavigated TUS (outside of the MRI scanner) is an efficient tool to carry out neuromodulation procedures in non-human primates. We found that, following neuronavigated rTUS, saccades were significantly modified, resulting in shorter latencies compared to no-rTUS trials. This behavioral modulation was maintained for up to 20 min. Oculomotor behavior returned to baseline after 18–31 min and could not be significantly distinguished from the no-rTUS condition. This study is the first to show that neuronavigated rTUS can have a persistent effect on monkey behavior with a quantified return-time to baseline. The specificity of the effects could not be explained by auditory confounds.
Transcranial Ultrasound Stimulation (TUS) is a safe, focused and non-invasive method of brain stimulation that has emerged as a novel and valuable tool for studying brain function in humans and animals. The use of TUS gathered momentum in the first decade of the new millennium, using low frequency, low intensity ultrasound waves to stimulate rodent primary motor cortex that generated motor responses (Yoo et al., 2011; King et al., 2013; Younan et al., 2013; Li et al., 2016; Ye et al., 2016) without damaging brain tissue (Yoo et al., 2011). In the same vein, TUS has been shown to immediately alter electromyographic and electroencephalographic measurements in sheep (Lee et al., 2016) and humans (Legon et al., 2014; Lee et al., 2016), and suppress somatosensory evoked potentials in swine (Dallapiazza et al., 2018). Recently, our group has shown in awake non-human primates that single pulses of TUS can modulate visuomotor behavior (Deffieux et al., 2013) with brief modulating effects (∼100 ms) from TUS in single neuronal responses (Wattiez et al., 2017). Whether such modulations due to TUS can be temporally extended and controlled remains to be demonstrated.
Accurate transducer positioning is crucial when sonicating target sites for desired stimulation effects. Presently, TUS studies of higher cognitive function demonstrate high inter-subject and between-group variability since anatomical variability and function may not correlate well in many brain areas (Brett et al., 2002). Clearly, localization of a subject’s individual anatomical brain region as well as precise placement and angle of the transducer is indispensable for successful targeting. However, skull differences as well as individual variability of the cerebral sulci have shown variations of up to 20 mm in the different axes, with some electrode positions having larger variability than others (Herwig et al., 2003; Okamoto et al., 2004). A different method is sometimes used to position the TMS coil by employing a functional-guided approach. If the TMS effect on the performance of a certain behavioral task is known, this task can serve as a “functional” probe to position the coil with subsequent targets or tasks (Göbel et al., 2001). Such “hunting” procedures can, however, be time-consuming due to the fact that different locations need to be tested by trial-and-error. Furthermore, the interpretation of the results with respect to brain anatomy is limited by the fact that the coil position is determined functionally and not anatomically.
To better account for inter-individual anatomical differences, image-guided frameless stereotaxic neuronavigation systems (SNS) have been used with Transcranial Magnetic Stimulation (TMS) in both humans and non-human primates (Schoenbaum et al., 2003; Sparing et al., 2007; Gerits et al., 2011). SNS use the subject’s individual MRI for navigation via a subject-image co-registration procedure based on facial/cranial landmarks. Although there can be technical limitations to this procedure as the quality of the MRI and limitations of the position sensor, the targeting error is within several millimeters (Neggers et al., 2004; Schönfeldt-Lecuona et al., 2005). In non-human primates, the accuracy of transducer placement can be improved to within 1 mm when using rigid fiducial markers during the co-registration of the subject to their scan (Frey et al., 2004; Figure 1A). SNS have been used in several ways: first, the gross anatomy of the cerebral cortex itself may serve as a reference system. However, reliable anatomical landmarks exist for very few target sites and individual differences of gyral folding and cortical layering have to be taken into consideration. The functional identification of a brain area with SNS requires that the same individual participates in both the fMRI and TUS study. Positioning of the coil over the local maxima of the fMRI activation cluster no longer assumes a fixed relationship between anatomical landmarks and task-related functional activations. Alternatively, functional neuroimaging data obtained from the literature can be used for navigation (“probabilistic approach”) (Paus et al., 1997). Although there are undoubtably slight anatomical variabilities amongst inter-individual regions in primate brain there is still predictable functional anatomical regions within primates (Amiez et al., 2019); this method takes advantage of and that has a relatively high consistency in the location of task-related “activations” across individuals.
Figure 1. Experimental set-up. (A) List of navigation tools to carry out a monkey registration: Brainsight computer (a), position sensor (b), calibration block (c), tool tracker (d), navigation pointer (e), subject tracker (f), and transducer tracker (g). (B) Subject registration. Top: registration can be based on facial landmarks which are defined as the bridge of the nose, the tip of nose and the notch above the tragus of the ears. Bottom: In the present study, we used fiducial markers which are fixed to the headpost of the subject resulting in higher co-registration accuracy between the subject and the MRI scan. These fiducial markers create hyperintensities on the MR images that are easily distinguishable. (C) Antisaccade paradigm. The animal starts the trial by fixating a central spot (green square) for 500–1000 ms. A peripheral cue (a red square at 16 degrees of eccentricity) is then presented either at the left or right of the visual field on the horizontal meridian for 1000 ms. The animal is rewarded by making an eye movement to the opposite location of the peripheral cue (antisaccade). The brain area stimulated by rTUS is indicated by the red circle (here, left FEF). Each animal practices a block of 100 trials before rTUS (pre-rTUS block) and then six post-rTUS trials. Each block contains 100 trials (50 right and 50 left) and lasts on average 5 min. The rTUS train lasts 20 s at a frequency of 10 Hz with stimulation pulses of 320 kHz. (D) Experimental set-up. The ultrasonic transducer is manually guided toward the Region of Interest (ROI) by a neuro-navigation system (Brainsight Vet). The cigar-shaped focal stimulation (bottom panel) targeted the FEF and SEF location [indicated by red dots: FEF (#) ans SEF (∧)]. We also performed sham stimulation sessions by placing the transducer over the primary motor or visual cortex [V1 (*) and M1 (**) indicated by blue dots]. The insets represent the MRI images used to locate ROI (red squares) and control areas (blue squares). The intersection of horizontal and vertical green lines in each quadrant indicates the point at which the focused ultrasound stimulations were directed.
Many studies optimized TMS setups in order to induce a modulatory effect which outlasts the stimulation, either for stand-alone treatments (Rossini et al., 2015), or combined with neuroimaging (Bestmann et al., 2008; Siebner et al., 2009).
Producing inhibitory effects that can outlast the period of stimulation for minutes or hours is crucial for rehabilitative purposes (Rossi and Rossini, 2004). It is typically achieved with repetitive stimulation (Wassermann and Lisanby, 2001; Hoffman and Cavus, 2002), each repeated stimulation increasing the magnetic energy delivered to the target. Increasing the acoustic energy is constrained by thermal considerations, as acoustic energy can be transferred to tissue (McDannold et al., 2010) and bone (Pinton et al., 2011), especially when considering that numerical models showed that the minimum activation energy was six orders of magnitude higher to elicit an action potential with ultrasound as compared to direct current (Krasovitski et al., 2011; Plaksin et al., 2014). Using repetitive pulses of ultrasound decreases the thermal rise (Constans et al., 2018). For these reasons, repetitive pulses have been proposed here and tested.
In this study, we present a detailed account of the effects of repetitive Transcranial Ultrasound Stimulation (rTUS) on the oculomotor regions of awake macaque monkeys while the animals were performing an antisaccade task. In primates, the frontal eye fields (FEF) and the supplementary eye fields (SEF) are two sophisticated cortical brain areas that play important roles in the control of visual attention and eye movements (Schall, 2015). Electrophysiological approaches exploring these areas have provided causal evidence about the role of the FEF, with direct bearing on oculomotor functions (Bruce and Goldberg, 1985), while a more indirect role has been proposed for the SEF (Schall, 2004, 2015). In addition to seeking new visual targets, an important part of saccadic behavior is to suppress eye movements that would be made to novel but behaviorally irrelevant stimuli. To investigate such control of voluntary vs. reflexive saccades, a test paradigm called the antisaccade task has been developed (Condy et al., 2005). In this task the subject is required to suppress a saccade toward a stimulus that appears in their periphery but instead must generate a voluntary saccade toward the opposite visual field (antisaccade).
In humans and macaques, functional imaging studies have shown that both the FEF and SEF are activated bilaterally during antisaccades (Schlag-Rey et al., 1997; Brown et al., 2006, 2007). Furthermore, patients with FEF lesions display a normal percentage of errors on the antisaccade task, but their correct antisaccades exhibit an increased latency (Pierrot-Deseilligny et al., 1991). Stimulation over FEF with single pulses of TUS or TMS have reported saccade latencies to be significantly delayed (Grosbras and Paus, 2002; Deffieux et al., 2013). These results have been interpreted as a modulation of saccadic inhibition in the contralateral visual field due to disruption of neural processing in the stimulated FEF area. Using repetitive TMS (rTMS), previous studies have described facilitating effects (Crevits et al., 2005) when compared to single pulse parameters (Olk et al., 2005), speculating that repetitive stimulation may more strongly interfere with the processing of inhibition, particularly active in antisaccades compared to single pulses (Everling and Munoz, 2000). Our current results establish rTUS as an innovative neurostimulation methodology by showing (1) that rTUS induced a reduction of saccade latencies only on ipsilateral antisaccades; (2) that this effect was reversible with a return to baseline after 18–31 min; and (3) that the latency decrements were dependent on neurostimulation locations.
Materials and Methods
The Stereotaxic Navigation System (SNS)
The Brainsight Vet stereotaxic navigation system (SNS, Figure 1A) was created by Rogue Research Inc. (Qc, Canada), and is composed of a standard computer with the neuronavigation software, an optical position sensor (Vicra-Northern Digital Inc., Ont, Canada), an articulated arm to hold the transducer, a subject and transducer tracker, a navigation pointer, as well as a calibration block for registering the acoustic focus of the transducer. Once the subject’s MRI data are transferred to the navigation computer, 3D reconstructions of the brain are created, atlas spaces can be introduced (e.g., MNI or Tailairach), and regions of interest are selected along with trajectories to brain targets. The primary reason for using the Brainsight navigation system was to aid in sonicating intended brain regions, and although TUS is seen as a non-invasive procedure, reducing the number of unnecessary stimulation sessions. An important advantage in using the MRI guided navigation system is that targets and trajectories to targets can be planned well in advance of the actual experimental procedures; a refinement to existing sonication protocols. During a navigated TUS session, the researcher is able to refer to the saved sonication locations in real-time on the MRI reconstructed images of the subject, simplifying data collection.
There are generally two methods to register the subject to their high resolution T1-weighted MRI brain scan prior to a neuronavigated TUS procedure. The first is based on skin landmarks while the second method uses a fiducial marker system that is rigidly attached to the animal’s skull. The first procedure, mostly used with human subjects, identifies anatomical landmarks on the 3D reconstructed skin within the software and homologous points on the subject’s skin prior to a sonication session. Recommended landmarks are defined as the bridge of the nose, the tip of nose, and the notch above the tragus of the ears (Figure 1B). Unfortunately, the accuracy of this type of registration leads to an error of approximately 2–5 mm with one of the difficulties being identifying the exact landmark on both the skin and the MRI scan (Ruohonen and Karhu, 2010). This method can also be used for non-human primate subjects. A second procedure is to fix fiducial markers on the subject before the MR scan; this method is typically more accurate as the landmarks are rigid. The present study used rigid markers for registration and is the most common method used with neuronavigation with non-human primate researchers. Prior to MRI, the animal is anesthetized and fitted with an MRI compatible head-fixation implant to ensure head immobilization during sonication. Attached to this headpost is an adaptor that is fitted with MRI compatible disks that hold adhesive radio-opaque fiducial markers (Figure 1B). These fiducial markers create hyperintensities on the MR images (Figure 1B) and are also easily identified on the animal. The neuronavigation software uses the positions of these fiducial markers to compute the coordinates of the neuronavigated TUS session and are necessary for subsequent co-registrations of the subject to the MR scan. The fiducial markers are only needed for the MRI and during TUS registration in the testing suite. This method has an error of approximately 1 mm based on phantom MRI tests (Frey et al., 2004). In order to calibrate the focal point of the transducer, a tool tracker is rigidly fixed to the TUS device. A jig is placed on the Brainsight calibration block to hold the transducer steady above the calibration pin in front of the optical position sensor (Figure 1A). Once calibrated, the TUS transducer can be moved over the head of the subject corresponding to desired pre-planned stimulation sites inside of the brain. The corresponding locations are visualized and updated in real-time on the computer monitor as the transducer is navigated into the correct orientation and trajectory of the desired stimulation sites.
Data were collected from three captive-born male macaques (one Macaca fascicularis; monkey G, 8 kg; two Macaca mulatta; monkey L and monkey S, 13, and 8.5 kg, respectively). Monkeys were paired-housed and handled in strict accordance with the recommendations of the Weatherall Report on good animal practice. Monkey housing conditions, surgical procedures and experimental protocols were all carried out in strict accordance with the National Institutes of Health guidelines (1996) and the recommendations of the EEC (86/609) and the French National Committee (87/848). Our experiments were authorized by the Animal Health and Veterinary Medication Division of the Department of Public Veterinary Health, Nutrition and Food Safety of the French Ministry of Health (last renewal N°: DTPP 2010-424). Our routine laboratory procedures included an environmental enrichment program where monkeys had access to toys, mirrors and swings. Monkeys also had visual, auditory and olfactory contact with other animals and, when appropriate, could touch and groom each other.
Under isoflurane and aseptic conditions, we surgically implanted a titanium headpost on each animal. Any possible pain associated with surgeries was pharmacologically ameliorated by means of a daily injection of ketofen (0.03 ml/kg) or buprenorphine (0.067 ml/kg). An institutional veterinarian regularly monitored the environment and wellbeing of the monkeys. Before participating in this study, animals had been periodically chaired, head-posted and trained to perform a series of tasks for a period of 3–6 months, until they became regular and proficient performers.
Experimental Apparatus for Non-human Primate Antisaccade Task Sessions
The animals were seated in a primate chair (Crist Instrument, Hagerstown, MD, United States) with their head fixed in a darkened booth. Animals were positioned 52 cm in front of a 60 Hz monitor. Eye movements were recorded with an infrared eye tracker (Eyelink 1k, SR-Research, Ont, Canada), and eye position was digitized and sampled at 1000 Hz and stored for off-line analysis. Visual paradigms and data acquisition were controlled by a computer running a real-time data acquisition system (Rexeno software; for further details see Krasovitski et al., 2011). Saccades were detected using custom Matlab scripts that first searched for significantly elevated velocity (>30°s–1). Saccade initiation and termination were then defined as the beginning and end of the monotonic change in eye position lasting 12 ms before and after the high-velocity gaze shift.
Prior to the first experimental session, animals were specifically trained in an anti-saccadic paradigm (Figure 1C), in which they were required to initially fix on a central brown target. Between 500 and 1000 ms after fixation onset, simultaneous to the disappearance of the central brown stimulus, and without a time gap, a red square appeared for 1000 ms at a location 16° to the right or left of it. Monkeys were trained not to look at this peripheral target but instead initiate a saccade in the opposite direction as soon as possible (i.e., an antisaccade). Monkeys were rewarded if the saccade fell within a 12°× 12° window centered at the mirror location of the visual target. Failure to trigger a saccade within 1000 ms after target onset canceled the trial and was considered an error trial. To analyze the performance in the antisaccade task, we divided the behavioral data into unambiguous categories. An eye movement centered in the mirror location of the 12°× 12° window visual target was defined as a correct trial, any other eye movement was defined as an error. The antisaccade task was chosen since prior human and monkey TMS experiments have revealed prosaccade paradigms to be much less sensitive to single pulse TMS interference than antisaccades (Olk et al., 2005). We computed saccade latency on rewarded antisaccades trials only and all non-rewarded trials were excluded from these latency comparisons.
Positioning the rTUS Transducer
We used a single element focused ultrasound transducer (H115, Sonic Concept, Bothell, WA, United States; central frequency 250 kHz, diameter 64 mm, FD# 1) in the current experiment. Even though the manufacturer specifies a 250 kHz central frequency, the transducer can be operated at 250, 320, 850, and 1380 kHz (Constans et al., 2017). We decided to use a 320 kHz frequency, as previously validated for modulating non-human primates behavior (Deffieux et al., 2013). A coupling cone (C103, Sonic Concepts, Bothell, WA, United States) filled with degassed water was placed between the transducer and the animal’s head (Figure 1D). The transducer was fixed on a mechanical arm with four rotational axes (Viewmaster LCD, Osmond Ergonomics, Wimborne, United Kingdom) to enable flexible positioning of the transducer over the head. A thin layer of echographic gel (Aquasonic 100, Parker Laboratories Inc., Fairfield, NJ, United States) was applied to the skin and the membrane of the coupling cone to ensure acoustic coupling (Figure 1D). Based on simulations done in water at 320 kHz, the −6 dB ovoidal focal spot was cigar shaped with a size of 4.7 × 4.7 × 33 mm (Figure 1D).
Acoustic and Thermal Modeling
The acoustic wave propagation of our focused ultrasound protocol was simulated using a k-space pseudospectral method-based solver, k-Wave (Cox et al., 2007) to obtain estimates for the pressure amplitude, peak intensity, spatial distribution, and thermal impact at steady state. 3D maps of the skull were extracted from a monkey CT scan (monkey L (Constans et al., 2017), 0.14 mm slice resolution, 0.33 mm slice distance). Water acoustic values were ρwater = 1000 kg.m−3 and cwater = 1500 m.s−1 and soft tissues were assumed to be homogeneous, with the following acoustic values: ρtissue = 1030 kg.m−3 and ctissue = 1560 m.s−1.
In the bone, a linear relationship between the Hounsfield Units (HU) from the CT scan and the sound speed, as well as the density, was used. The power law model for attenuation was where the porosity Φ is defined by in the skull (Aubry et al., 2003). The attenuation coefficients for the acoustic propagation αmin and αmax depended on the frequency: αmin = αmin0fb with αmin0 = 0.2dB.cm−1.MHz−b and αmax = αmax0fb with αmax0 = 8dB.cm−1.MHz−b (Aubry et al., 2003). We set the parameters to ρmax = 2200kg.m−3, cmax = 3100m.s−1 (Younan et al., 2013), β = 0.5 (Aubry et al., 2003), b = 1.1 (Constans et al., 2018). The attenuation coefficients in bone accounted for both absorption and scattering (Pinton et al., 2011). In a consistent manner, in soft tissues, attenuation coefficient for the acoustic propagation αtissue depended on the frequency: αtissue = αtissue0fb with αtissue0 = 0.6dB.cm−1.MHz−b.
The propagation simulation was performed at 320 kHz with a 150 μs-long pulse signal (enough to reach a steady state). The transducer was modeled as a spherical section (63.2 mm radius of curvature and 64 mm active diameter). The acoustic focus of the transducer was positioned on the target (Left FEF). The simulated pulses were spatially apodised (r = 0.35) on the spherical section. Ultrasound propagates first through water before entering the skull cavity with the geometrical focal point located below the surface, inside the brain. Simulations were performed in free water, and the maximum amplitude obtained was used to rescale the results in skull. The thermal modeling is based on the bio-heat equation (Legon et al., 2014):
where T, ρ, C, κ, and q are the temperature, density, specific heat, thermal conductivity and rate of heat production, respectively. Heat production is defined as , αabs being the absorption coefficient and P the peak negative pressure. κ is set to 0.528 W.m–1.K–1 in soft tissue and 0.4 W.m–1.K–1 in the skull; C is set to 3600 J.kg–1.K–1 in soft tissue and 1300 J.kg–1.K–1 in the skull. In the tissue, the absorption coefficient was set to (Goss et al., 1979). In the skull the longitudinal absorption coefficient is proportional to the density with (Pinton et al., 2011). The last term corresponds to the perfusion process: w,ρb,Cb, and Ta correspond to the blood perfusion rate, blood density, blood specific heat and blood ambient temperature, respectively. These parameters are assumed homogeneous over the brain, although a more detailed description of the brain cooling processes can be found in the literature (Wang et al., 2016). The perfusion parameters are based on previous reports (Pulkkinen et al., 2011): w = 0.008 s–1; rb = 1030 kg.m–3; Cb = 3620 J.kg–1.K–1 and Ta = 37°C.
The bioheat equation is solved by using a 3D finite-difference scheme in MATLAB (Mathworks, Natick, United States) with Dirichlet boundary conditions. Initial temperature conditions were 37°C in the brain, skull and tissue, and 24°C in the water coupling cone. Simulations were run over 1 min pre-sonication, followed by 20 s of sonication, and 5 min post-sonication, closely following the experimental procedure.
To quantify the pressure amplitude, peak intensities, spatial distribution, and potential temperature changes in the monkey brain associated with the TUS protocol targeting the FEF used in this study we simulated the acoustic wave propagation and its thermal effect in a whole head finite element model based on a high-resolution monkey CT scan. As estimated by these numerical simulations, the maximum spatial-peak pulse-averaged intensity (Isppa) at the acoustic focus point was 21.2 W/cm2 for the left FEF target (spatial peak temporal average intensities, Ispta: 6.4 W/cm2). Given that the skull is more acoustically absorbing than soft tissue, the highest thermal increase is located in the skull itself, estimated by the simulation to be 1.53°C. For an approximate 0.5 mm thickness of the dura the maximum temperature in the brain below the dura was 37.4°C. The maximal thermal increase at the geometrical focus of the sonic transducer was less than 0.1°C.
Repetitive Transcranial Ultrasound Stimulation Parameters
The ultrasound frequency was set to 320 kHz. The pulse duration was 30 ms with a rise and fall time set to 1 ms to avoid abrupt changes in pressure that generate brief white noise. The pulse repetition frequency was set to 10 Hz and the total sonication time was 20 s. The signal was generated by a TiePie generator (Handyscope HS5). A 75-watt amplifier (75A250A, Amplifier Research, Souderton, PA) was used to deliver the required power to the transducer and the input voltage of the transducer was monitored using a voltage probe (P6139A, Tektronix, Melrose, MA) connected to a TiePie oscilloscope. The amplifier gain was set to deliver an output voltage Vout = 173 V peak-to-peak, corresponding to a pressure amplitude of 0.76 MPa in water (calibrated with an in-house heterodyne interferometer; Constans et al., 2017). In order to estimate the skull attenuation of the ultrasound beam, a clean and degassed primate (Macaca mulatta) skull specimen was placed in front of the transducer in a degassed water tank and the pressure at the focus was estimated using a custom-built heterodyne interferometer (Royer and Casula, 1994) which uses a laser beam to detect the vibration of the ultrasound wave on a Mylar membrane. The amplitude of the vibration is then converted to pressure with a 10–4 Å/(Hz)1/2 sensitivity and a 20 MHz bandwidth (Royer and Dieulesaint, 1986). The transmission of the pressure through the degassed primate skull was assessed at seven different points arbitrarily chosen on the skull. The transmission ratio was found to be 58 ± 8% (standard deviation). The in-situ pressure delivered to the monkey brain transcranially was subsequently estimated at 0.44 ± 0.06 MPa. The equivalent mechanical index (MI) value was 1.3 with an intensity spatial peak pulse average (ISPPA) of 19 W/cm2 in free water. These values were attenuated to MI = 0.8 ± 0.1 and ISPPA = 6.5 ± 1.8 W/cm2, respectively inside the primate skull. Taking into account the pulse duration and PRF (respectively, 30 ms and 10 Hz, corresponding to a 30% duty cycle) during the neuro-stimulation sequence, the intensity spatial peak time average (ISPTA) was estimated to be 5.7 W/cm2 in free water and 1.9 W/cm2 inside the primate skull.
We targeted the FEF, SEF, cortical motor, and visual cortex with the navigated transducer using Brainsight Vet neuronavigation (Rogue Research Inc., QC, Canada). For each animal, the target areas were identified based on the anatomical MRI data. The FEF was identified as the fundus of the arcuate sulcus in front of the spur, along the anterior bank of the arcuate sulcus (Schall et al., 2009) and SEF neurostimulation site was anatomically centered along the dorso-medial part of the frontal convexity cortex (Parthasarathy et al., 1992). In monkey L, the target ROI was placed in the right FEF, whereas the control was placed in the right hemisphere of the visual cortex (V1). In monkey S, the target ROI was set in the left FEF and the control region in the left hemisphere of the visual cortex (V1). In monkey G, the target ROI was set in the SEF of the left hemisphere, while the control region was set closer to the ROI in motor cortex (left M1) due to the large mass of muscle on the animal’s skull. SEF, motor and visual cortex areas were defined on the MRI visually and according to stereotaxic coordinates (Figure 1D).
In this study, animals performed a total of 40 sessions (10 stimulated sessions “rTUS” and 10 non-stimulated sessions “no-rTUS” in the regions of interest “ROIs” and in the control regions). Each experimental session contained one practice block, followed by a total of 7 blocks of 100 antisaccades (50 on each side randomly distributed). Navigated repetitive Transcranial Focused Ultrasound Stimulation (rTUS) was only delivered after the first block of trials (1 of 7) was completed. In this study, the first block is referred to as pre-rTUS and the six following blocks as post-rTUS (Figure 1C). In order to keep conditions across all experimental sessions as similar as possible, the transducer was placed on the target regions (ROI and control) in both rTUS and non-rTUS sessions.
Saccade latencies were compared across different factors: (1) the side of the antisaccade (i.e., contralateral or ipsilateral to the stimulation location), (2) rTUS stimulation or no stimulation and (3) the stimulated area (ROI vs. Control), by first normalizing the saccade latencies with the transformation log 1/RT where RT defined the reaction time or saccade latency. We excluded from theses analyses all the saccades in which the latencies were inferior to 100 ms. We performed a 3-way ANOVA on the normalized RTs using the three factors defined previously (Side, rTUS or the area) for each animal. For the timeline comparison of the mean latencies for no-rTUS and rTUS trials, we used a non-parametrical test with a correction for multiple comparisons (Benjamini and Hochberg, 1995).
Global Behavioral Performance During Antisaccade Task
All three monkeys performed the task across a total of 20 behavioral sessions for each targeted area (ROI or control) with an average error rate between 3 and 14% (Figure 2A). Monkey L and monkey G presented significantly higher error rates on the ipsilateral side of the FEF/SEF stimulation. Nevertheless, these increases in error were very low (monkey G from 2.5% to 3.5% for SEF ipsilateral side, and monkey L from 4.6 to 7.2% for FEF ipsilateral side), suggesting that the animals were performing the task uniformly during both rTUS and no-rTUS sessions. In addition, it suggests that rTUS did not disturb the animals’ general performance when either the frontal or the primary cortices were the targets of the neurostimulation. Monkey S, however, behaved differently and presented a higher error rate in all conditions (more than 9% for contra/ipsi side and control/tested areas). Our interpretation is that monkey S was using a different strategy with a higher guessing rate than the other animals for the antisaccade target location. In the context of antisaccades, a guessing strategy generates high error rates and often a bimodal latency distribution.
Figure 2. Behavioral performances during the antisaccade task. (A) Percentage of error rates for all three animals (monkey G; left quadrant, monkey L; middle quadrant and monkey S; right quadrant). Each bar represents the mean error rates for the post-rTUS blocks for the different conditions (violet/orange for contralateral vs. ipsilateral antisaccades, respectively, and dark/light and violet/green for rTUS vs. no-rTUS, respectively). Error bars indicate the standard error of the mean (SEM). We compared rTUS vs. no-rTUS mean error rates for each condition using a Wilcoxon rank sum test (*p < 0.05). Note that the performance for all conditions (ROI, Control, Contra, and Ipsi) were similar (p > 0.05) for all three animals except for monkeys G and L who made significantly more errors when rTUS was directed toward FEF/SEF regions for ipsilateral antisaccades. (B,C) Antisaccade latencies (ms) for all conditions [Contra/Ipsi, rTUS/no-rTUS for each animal (Monkeys G, L, and S), when the transducer was directed toward the oculomotor (B) or the control (C) areas]. The latencies were pooled over the pre-stimulation block (practice) and the next six post-rTUS blocks for each condition. The percentages indicate the percentage of anticipated saccades (<100 ms, dashed lines) over the total number of trials for each condition. For each histogram, the downside arrow indicates the median latency (ms) with the cutoff of anticipations.
To confirm our theory, we analyzed the latency distributions of all monkeys (Figures 2B,C). The antisaccade latency for each individual trial was calculated as the time between the stimulus presentation and the onset of the eye saccade velocity (30°s–1). We compared saccade latencies between trials depending on different conditions: (1) whether it was a rTUS or no-rTUS session; (2) whether it was an ipsilateral or a contralateral antisaccade; and (3) whether the transducer was directed toward the oculomotor (FEF or SEF) or the control areas (V1 or M1). Figure 2B represents all the histograms for the different conditions (rTUS/no-rTUS: dark orange/violet vs. light orange/green; contra/ipsi: orange vs. violet/green; and ROI/Control: Figure 2B vs. Figure 2C) for all three animals.
Monkeys G and L both showed a unimodal distribution of the saccade latencies for ROI and control areas with very few anticipated saccades (i.e., saccade with latencies < 100 ms). Monkey G almost never anticipated the saccade (only 2 out of 26,821 trials) and monkey L had a relatively low number of anticipated saccades (<4%). In comparison, the distribution of saccade latencies for monkey S differed significantly presenting larger proportions of anticipated saccades (14–21%), suggesting that monkey S used a different strategy than the other two animals. The distribution of saccade latencies for monkey S confirmed that this animal was guessing and anticipating the appearance of the target.
We computed an ANOVA on normalized saccade latencies for each animal with the following main factors: the stimulation type (rTUS/no-rTUS), the side of the antisaccade (contra/ipsi) and the area (ROI/Control). Based on the ANOVA, we first reported the results of the main factors. We noticed that monkey G and monkey S showed an important lateralized bias for the contralateral side with saccades with significant shorter latencies (ANOVA, p < 0.001). For monkey L, a similar bias was present but only for the ipsilateral side, with saccades having significant shorter latencies (ANOVA, p < 0.001). Concerning the effect of rTUS stimulation, rTUS induced shorter latencies in all three animals (p < 0.001 for monkey G and monkey S, and p = 0.0064 for monkey L). Monkey S showed similar latencies for saccades when the transducer was directed toward FEF or the visual cortex (ANOVA, p = 0.192) although the saccades were made significantly earlier when rTUS was focused on the frontal areas (FEF/SEF) compared to the control sites (visual or motor cortex, ANOVA, p < 0.001). Concerning the interactions, for monkey G and monkey S, rTUS did shorten the latencies for contralateral and ipsilateral sides for both ROI and control areas (ANOVA, interaction Bonferroni, p < 0.05) although monkey L had no significant interaction for rTUS on saccade latencies.
Sustained Effect of rTUS on Ipsilateral Anti-saccades
We showed that rTUS has a main effect of inducing shorter saccades. We now wanted to investigate how long the effect of rTUS lasted. Each animal performed a complete block of trials before rTUS stimulation (block #0). We then collected the saccade latencies during six additional blocks of trials (block #1 to block #6). The timeline of the mean latencies for all three monkeys for ipsilateral and contralateral antisaccades are displayed in Figure 3. We compared for each block of trials the no-rTUS and r-TUS saccades latencies with a non-parametric Wilcoxon test using multiple corrections (Benjamini and Hochberg, 1995). For monkey G, a significant effect of rTUS was observed in the first four blocks of trials for ipsilateral antisaccades only when the transducer was directed toward frontal regions (Figure 3, monkey G – SEF). For monkey L, rTUS showed no effect in the first block of trials (block #1) but was significant in the next three blocks (block #2 to block #4) with shorter latencies saccades on the ipsilateral side. On average, four blocks of trials correspond to a duration of between 21 and 30 min (mean ± SEM; 21 min ± 1 for monkey G and 30 min ± 0.4 for monkey L). After this time period, the saccade latencies were statistically indistinguishable from the latencies recorded in the sham condition (no-rTUS). For contralateral trials, only a few sporadic blocks of trials showed significantly shorter latencies (block #4 for monkey G and none for monkey L) for the rTUS condition. No significant behavioral change was seen when the transducer was directed toward the control regions (Figure 3, control areas; M1 and V1). For both monkey G and monkey L, as seen in Figure 2, we did not notice any difference in any of the six post-rTUS blocks for both contralateral and ipsilateral trials (Wilcoxon rank sum test, p > 0.05). These results confirm that rTUS during an antisaccade task only affect the behavioral performance of the animals when targeted in their frontal cortex and not their visual/motor cortex.
Figure 3. Time course of rTUS effect on mean latencies. Mean latencies per block when rTUS was directed toward the oculomotor areas (SEF/FEF; two left columns) or the control areas (M1/V1; two right columns) in three animals (monkey G; top row, monkey L; middle row and monkey S; bottom row). Each bar represents the mean latency averaged over all successful antisaccades in each block of trials for each side (contra/ipsi), error bars indicate SEM. Open bars indicate non-significant differences (Wilcoxon rank sum test, p > 0.05); filled bars indicate significant differences (* for p < 0.05 and *** for p < 0.001, corrected for multiple comparisons). Gray backgrounds indicate when error rates were superior to 15% although light backgrounds indicate when error rates were inferior to 15%. Black arrows schematize when rTUS was performed between blocks #0 and #1.
We then wanted to better understand the behavioral performance of monkey S as we noticed that, contrary to the two others monkeys, monkey S performed significantly more anticipated saccades, suggesting the use of a different strategy. We focused our analysis on saccade latencies that were superior to 175 ms. We choose this threshold based on the histograms of the saccade latencies for monkey S (Figures 2B,C) as they were more similar to the ones for the other animals (even though approximately a third of the saccades were excluded from this analysis, 30% for the FEF neuronavigated rTUS sessions, and 38–39% for V1 neuronavigated rTUS sessions). We next wanted to check whether the effect of neuronavigated rTUS on early blocks after rTUS was present in these antisaccades. Similarly to the two other animals, the analyses of the time courses per block (Figure 3 bottom panel) showed a significant decrease of mean latencies for ipsilateral antisaccades when the transducer was directed toward the FEF locations for the first blocks of trial post-rTUS (block #1 to block #3; 18 min ± 0.5). No variations on ipsilateral saccades on early blocks were noticed with the visual cortex stimulations. However, monkey S showed significant differences for late blocks (block #4 to block #6) when the transducer was directed toward V1 triggering slower saccades for the no-rTUS sessions. Indeed, the mean latencies for these three last blocks of trials increased significantly compared to block #0. Moreover, monkey S showed poor performance during these blocks of trials (error rate > 15%). We also noticed slower latencies for late blocks for contralateral saccades for V1 stimulation sessions (blocks #2, #4, and #5).
This study demonstrates the feasibility of using focused ultrasound to modulate visual behavior for a sustained period of several minutes in the awake non-human primate brain. Here we show that continuous low frequency (320 kHz) and low-pressure amplitude (estimated at 0.4 MPa in situ corresponding mechanical index of 0.7) pulses of ultrasound lasting 20 s can be used in order to achieve these long-lasting effects. We have demonstrated that rTUS set with such parameters, applied to oculomotor regions of the frontal cortex, produces a subtle significant shortening of ipsilateral saccade latencies and a subtle increase of error rates, as compared to non-stimulated trials. Together, these results show that neuronavigated rTUS applied to oculomotor regions can have a sustained effect on saccades when they occur on the ipsilateral side of stimulation.
Single vs. Repetitive Stimulation
Compared to the modulation of response times reported with repetitive Transcranial Magnetic Stimulation (rTMS), our results appear to be less affected by baseline variability of response times as no normalization was required (Thut and Pascual-Leone, 2009). In addition, no discomfort (muscle twitching) was observed during ultrasonic neurostimulation, in contrast to that observed during rTMS sessions. In previous studies, while stimulating FEF, saccadic latencies after single pulses of TUS or TMS have been reported to be significantly delayed (Grosbras and Paus, 2002; Deffieux et al., 2013). As described in this study, the use of repetitive stimulation facilitates antisaccade triggering. These results are in accordance with findings observed following rTMS, demonstrating facilitation effects (Crevits et al., 2005) while compared to single pulse TMS (Olk et al., 2005). We propose that repetitive stimulation may significantly interfere, compared to single pulse, with the processing of inhibition, which is particularly active in antisaccade tasks (Everling and Munoz, 2000). Although the exact mechanisms are still unknown, ultrasonic neurostimulation is believed to be mechanical rather than thermal, as illustrated by experimental evidence on rodents, where lower frequencies yield smaller motor thresholds (Ye et al., 2016). Also, we speculate that our results following rTUS over primary visual cortex may interfere with the processing of visual activity in the superior colliculi and therefore may provoke an increased rate of anticipated saccades (here, observed with monkey S who performed a high rate of saccades with short latencies).
Altogether, our results support the use of this approach to study brain function and non-invasive neurostimulation for exploratory and therapeutic purposes with unprecedented spatial resolution. The spatial resolution is proportional to the ultrasonic wavelength (1.5 mm at 1 MHz) and decreases with the ultrasonic frequency. Increasing the frequency thus improves the precision of the stimulation but, unfortunately, it also induces a higher thermal rise due to increased tissue absorption (Constans et al., 2018). Furthermore, aberrations induced by the skull increase with frequency (Maimbourg et al., 2018b; Leung et al., 2019). Therefore, targeting brain regions non-invasively is easier with low frequency ultrasound (Yin and Hynynen, 2005). For frequencies higher than 300 kHz, CT-based (Marquet et al., 2009; Pulkkinen et al., 2014) or MR-based (Wintermark et al., 2014; Miller et al., 2015) aberration correction needs to be performed. This can be achieved with multielement arrays of transducers (Jeanmonod et al., 2012; Chauvet et al., 2015) or with single element transducers combined with an acoustic lens (Maimbourg et al., 2018a, 2019). With these latest developments, rTUS appears to be very promising and offers a novel and competitive neurostimulation technique. When compared to rTMS, rTUS offers a higher spatial resolution (Pinton et al., 2012), a larger targeting envelope (Krishna et al., 2018) and the absence of noise or mechanical vibration during stimulation (Mohammadjavadi et al., 2019). This technique shows promise in the non-invasive exploration of cognition with non-human primate models. Used with ultrasound contrast agents (Magnin et al., 2015; Pouliopoulos et al., 2019), neuronavigated rTUS can be used to deliver allosteric modulators locally as well as to transiently modulate brain activity (McDannold et al., 2015; Constans et al., 2019).
Application on Humans
Although the exact mechanisms of focused ultrasound neurostimulation are still unknown, the ability to modulate precisely defined deep brain structures non-invasively opens exciting possibilities in the treatment of psychiatric and neurological disorders in humans, such as treatment-resistant depression, anxiety disorders, Parkinson’s disease, essential tremor, and disorders of consciousness. The promise of TUS also lies in its unique ability to explore safely brain circuitry and mechanisms underlying complex cognitive functions, as our team has previously demonstrated in non-human primates (Fouragnan et al., 2019; Khalighinejad et al., 2020).
In current rTMS clinical trials, neuronavigation can be seen as a necessary add-on to take advantage of inter-subject anatomical variability and allow more precise targeting (Fitzgerald et al., 2009; Dollfus et al., 2017). Even though several cortical-rTUS human studies have used non-navigated procedures (Mueller et al., 2014; Sanguinetti et al., 2020), neuronavigation should be seen not only as a valuable tool in deep-rTUS, but also a requirement. Indeed, given its millimeter accuracy and its capability of hitting deep targets (without stimulating unintended brain structures dorsal to stimulation), neuronavigated rTUS may be the ideal non-invasive deep brain stimulation tool. These two innovations have already been coupled in several proof-of-concept human studies in healthy volunteers (Legon et al., 2018a, b) and patients (Brinker et al., 2020). Interesting is the fact that clinicians are now using fMRI and neuronavigated rTMS to modulate different target areas in the brain of patients depending on the neural circuits of these individuals (Siddiqi et al., 2020). The use of fMRI, along with navigated rTUS in humans, may enhance even more our ability to provide personalized and supportive-based neuromodulation therapy. Our team took advantage of this technique and used fMRI data along with navigated rTUS in non-human primates (Khalighinejad et al., 2020). Like any new technique, the safety of the ultrasound parameters needs to be further investigated before translating these novel stimulation protocols to the clinic. In particular, ultrasound induces a higher thermal rise at higher frequencies due to increased tissue absorption (Duck, 1990). Numerical models estimating the thermal rise in tissues has been developed (Constans et al., 2018) and experimentally validated (Ozenne et al., 2019) to help design novel TUS systems. The safety of ultrasonic neurostimulation has additionally been assessed by histology in rodents (Yoo et al., 2011; Yang et al., 2012; Mehić et al., 2014; Li et al., 2016), pigs (Dallapiazza et al., 2018), sheep (Gaur et al., 2020) and monkeys (Blackmore et al., 2019; Verhagen et al., 2019; Gaur et al., 2020).
This study demonstrates that long lasting (20 s) repetitive transcranial ultrasound stimulation modulates visual behavior of awake non-human primates for a sustained period of several minutes. rTUS was specifically targeted to the Frontal Eye Field and the Supplementary Eye field with a frameless neuronavigation system. Neuronavigated rTUS holds promise in exploring brain circuitry by offering novel treatment modalities based on sustained neuromodulation.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
The animal study was reviewed and approved by the recommendations of the EEC (86/609) and the French National Committee (87/848). Our experiments were authorized by the Animal Health and Veterinary Medication Division of the Department of Public Veterinary Health, Nutrition and Food Safety of the French Ministry of Health (last renewal N°: DTPP 2010-424).
J-FA and PP designed the study. HA, JC, and CC performed the animal experiments. J-FA, CC, and DA performed the acoustic calibrations. SF, PP, and HA set up the neuronavigator. FA, HA, and PP analyzed the data. FA, PP, SF, DA, and J-FA wrote the manuscript with input from all authors. All authors contributed to the article and approved the submitted version.
This work was supported by the Bettencourt Schueller Foundation, Paris, France (UltraBrain project), the “Agence Nationale de la Recherche” under the program “Future Investments” with the reference ANR-10-EQPX-15, IHU FOReSIGHT (ANR-18-IAHU-0001) and the European Research Council SYNERGY Grant scheme (HELMHOLTZ, ERC Grant Agreement #610110).
Conflict of Interest
SF was a full time employee and partner of Rogue Research Inc.
The remaining 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.
Amiez, C., Sallet, J., Hopkins, W. D., Meguerditchian, A., Hadj-Bouziane, F., Hamed, S., et al. (2019). Sulcal organization in the medial frontal cortex provides insights into primate brain evolution. Nat. Commun. 10:3437. doi: 10.1038/s41467-019-11347-x
Aubry, J. F., Tanter, M., Pernot, M., Thomas, J. L., and Fink, M. (2003). Experimental demonstration of noninvasive transskull adaptive focusing based on prior computed tomography scans. J. Acoustical Soc. Am. 113, 84–93. doi: 10.1121/1.1529663
Benjamini, Y., and Hochberg, Y. (1995). Controlling the false discovery rate: a practical and powerful approach to multiple testing. R. Stat. Soc. 57, 289–300. doi: 10.1111/j.2517-6161.1995.tb02031.x
Bestmann, S., Ruff, C. C., Blankenburg, F., Weiskopf, N., Driver, J., and Rothwell, J. C. (2008). Mapping causal interregional influences with concurrent TMS–fMRI. Exp. Brain Res. 191, 383–402. doi: 10.1007/s00221-008-1601-8
Blackmore, J., Shrivastava, S., Sallet, J., Butler, C. R., and Cleveland, R. O. (2019). Ultrasound neuromodulation: a review of results, mechanisms and safety. Ultrasound Med. Biol. 45, 1509–1536. doi: 10.1016/j.ultrasmedbio.2018.12.015
Brett, M., Johnsrude, I. S., and Owen, A. M. (2002). The problem of functional localization in the human brain. Nat. Rev. Neurosci. 3, 243–249. doi: 10.1038/nrn756
Brinker, S. T., Preiswerk, F., White, P. J., Mariano, T. Y., McDannold, N. J., and Bubrick, E. J. (2020). Focused ultrasound platform for investigating therapeutic neuromodulation across the human hippocampus. Ultrasound Med. Biol. 46, 1270–1274. doi: 10.1016/j.ultrasmedbio.2020.01.007
Brown, M. R. G., Goltz, H. C., Vilis, T., Ford, K. A., and Everling, S. (2006). Inhibition and generation of saccades: rapid event-related fMRI of prosaccades, antisaccades, and nogo trials. Neuroimage 33, 644–659. doi: 10.1016/j.neuroimage.2006.07.002
Brown, M. R. G., Vilis, T., and Everling, S. (2007). Frontoparietal activation with preparation for antisaccades. J. Neurophysiol. 98, 1751–1762. doi: 10.1152/jn.00460.2007
Bruce, C. J., and Goldberg, M. E. (1985). Primate frontal eye fields. I. Single neurons discharging before saccades. J. Neurophysiol. 53, 603–635. doi: 10.1152/jn.19184.108.40.2063
Chauvet, D., Marsac, L., Pernot, M., Boch, A.-L., Guillevin, R., Salameh, N., et al. (2015). Targeting accuracy of transcranial magnetic resonance–guided high-intensity focused ultrasound brain therapy: a fresh cadaver model. J. Neurosurg. 118, 1046–1052. doi: 10.3171/2013.1.JNS12559
Condy, C., Wattiez, N., Rivaud-Péchoux, S., and Gaymard, B. (2005). Ketamine-induced distractibility: an oculomotor study in monkeys. Biol. Psychiatry 57, 366–372. doi: 10.1016/j.biopsych.2004.10.036
Constans, C., Ahnine, H., Santin, M., Lehericy, S., Tanter, M., Pouget, P., et al. (2019). Non-invasive ultrasonic modulation of visual evoked response by GABA delivery through the blood brain barrier. J. Controlled Release 318, 223–231. doi: 10.1016/j.jconrel.2019.12.006
Constans, C., Deffieux, T., Pouget, P., Tanter, M., and Aubry, J.-F. (2017). A 200–1380-kHz quadrifrequency focused ultrasound transducer for neurostimulation in rodents and primates: transcranial In Vitro calibration and numerical study of the influence of skull cavity. IEEE Trans. Ultrason Ferroelect. Freq. Contr. 64, 717–724. doi: 10.1109/TUFFC.2017.2651648
Constans, C., Mateo, P., Tanter, M., and Aubry, J.-F. (2018). Potential impact of thermal effects during ultrasonic neurostimulation: retrospective numerical estimation of temperature elevation in seven rodent setups. Phys. Med. Biol. 63, 25003–25018. doi: 10.1088/1361-6560/aaa15c
Cox, B. T., Kara, S., Arridge, S. R., and Beard, P. C. (2007). k-space propagation models for acoustically heterogeneous media: application to biomedical photoacoustics. J. Acoustical Soc. Am. 121:3453. doi: 10.1121/1.2717409
Crevits, L., Van den Abbeele, D., Audenaert, K., Goethals, M., and Dierick, M. (2005). Effect of repetitive transcranial magnetic stimulation on saccades in depression: a pilot study. Psychiatry Res. 135, 113–119. doi: 10.1016/j.psychres.2003.10.008
Dallapiazza, R. F., Timbie, K. F., Holmberg, S., Gatesman, J., Lopes, M. B., Price, R. J., et al. (2018). Noninvasive neuromodulation and thalamic mapping with low-intensity focused ultrasound. J. Neurosurg. 128, 875–884. doi: 10.3171/2016.11.JNS16976
Deffieux, T., Younan, Y., Wattiez, N., Tanter, M., Pouget, P., and Aubry, J.-F. (2013). Low-intensity focused ultrasound modulates monkey visuomotor behavior. Curr. Biol. 23, 2430–2433. doi: 10.1016/j.cub.2013.10.029
Dollfus, S., Jaafari, N., Guillin, O., Trojak, B., Plaze, M., Saba, G., et al. (2017). High-frequency neuronavigated rTMS in auditory verbal hallucinations: a pilot double-blind controlled study in patients with schizophrenia. Schizophrenia Bull. 44, 505–514. doi: 10.1093/schbul/sbx127
Duck, F. A. (ed.) (1990). “Thermal properties of tissue,” in Physical Properties of Tissues, Amsterdam: Elsevier, 9–42. doi: 10.1016/B978-0-12-222800-1.50006-1
Everling, S., and Munoz, D. P. (2000). Neuronal correlates for preparatory set associated with pro-saccades and anti-saccades in the primate frontal eye field. J. Neurosci. 20, 387–400. doi: 10.1523/jneurosci.20-01-00387.2000
Fitzgerald, P. B., Hoy, K., McQueen, S., Maller, J. J., Herring, S., Segrave, R., et al. (2009). A randomized trial of rTMS targeted with MRI based neuro-navigation in treatment-resistant depression. Neuropsychopharmacol 34, 1255–1262. doi: 10.1038/npp.2008.233
Fouragnan, E. F., Chau, B. K. H., Folloni, D., Kolling, N., Verhagen, L., Klein-Flügge, M., et al. (2019). The macaque anterior cingulate cortex translates counterfactual choice value into actual behavioral change. Nat. Neurosci. 22, 797–808. doi: 10.1038/s41593-019-0375-6
Frey, S., Comeau, R., Hynes, B., Mackey, S., and Petrides, M. (2004). Frameless stereotaxy in the nonhuman primate. Neuroimage 23, 1226–1234. doi: 10.1016/j.neuroimage.2004.07.001
Gaur, P., Casey, K. M., Kubanek, J., Li, N., Mohammadjavadi, M., Saenz, Y., et al. (2020). Histologic safety of transcranial focused ultrasound neuromodulation and magnetic resonance acoustic radiation force imaging in rhesus macaques and sheep. Brain Stimul. 13, 804–814. doi: 10.1016/j.brs.2020.02.017
Gerits, A., Ruff, C. C., Guipponi, O., Wenderoth, N., Driver, J., and Vanduffel, W. (2011). Transcranial magnetic stimulation of macaque frontal eye fields decreases saccadic reaction time. Exp. Brain Res. 212, 143–152. doi: 10.1007/s00221-011-2710-3
Göbel, S., Walsh, V., and Rushworth, M. F. S. (2001). The mental number line and the human angular gyrus. Neuroimage 14, 1278–1289. doi: 10.1006/nimg.2001.0927
Goss, S. A., Frizzell, L. A., and Dunn, F. (1979). Ultrasonic absorption and attenuation in mammalian tissues. Ultrasound Med. Biol. 5, 181–186. doi: 10.1016/0301-5629(79)90086-3
Grosbras, M.-H., and Paus, T. (2002). Transcranial magnetic stimulation of the human frontal eye field: effects on visual perception and attention. J. Cogn. Neurosci. 14, 1109–1120. doi: 10.1162/089892902320474553
Herwig, U., Satrapi, P., and Schönfeldt-Lecuona, C. (2003). Using the international 10-20 EEG system for positioning of transcranial magnetic stimulation. Brain Topogr. 16, 95–99. doi: 10.1023/b:brat.0000006333.93597.9d
Hoffman, R. E., and Cavus, I. (2002). Slow transcranial magnetic stimulation, long-term depotentiation, and brain hyperexcitability disorders. AJP 159, 1093–1102. doi: 10.1176/appi.ajp.159.7.1093
Jeanmonod, D., Werner, B., Morel, A., Michels, L., Zadicario, E., Schiff, G., et al. (2012). Transcranial magnetic resonance imaging-guided focused ultrasound: noninvasive central lateral thalamotomy for chronic neuropathic pain. Neurosurg. Focus 32:E1. doi: 10.3171/2011.10.FOCUS11248
Khalighinejad, N., Bongioanni, A., Verhagen, L., Folloni, D., Attali, D., Aubry, J.-F., et al. (2020). A basal forebrain-cingulate circuit in macaques decides it is time to act. Neuron 105:370–384.e8. doi: 10.1016/j.neuron.2019.10.030
King, R. L., Brown, J. R., Newsome, W. T., and Pauly, K. B. (2013). Effective parameters for ultrasound-induced In Vivo neurostimulation. Ultrasound Med. Biol. 39, 312–331. doi: 10.1016/j.ultrasmedbio.2012.09.009
Krasovitski, B., Frenkel, V., Shoham, S., and Kimmel, E. (2011). Intramembrane cavitation as a unifying mechanism for ultrasound-induced bioeffects. Proc. Natl. Acad. Sci. U.S.A. 108, 3258–3263. doi: 10.1073/pnas.1015771108
Krishna, V., Sammartino, F., and Rezai, A. (2018). A review of the current therapies, challenges, and future directions of transcranial focused ultrasound technology. JAMA Neurol. 75, 246–254. doi: 10.1001/jamaneurol.2017.3129
Lee, W., Lee, S. D., Park, M. Y., Foley, L., Purcell-Estabrook, E., Kim, H., et al. (2016). Image-guided focused ultrasound-mediated regional brain stimulation in sheep. Ultrasound Med. Biol. 42, 459–470. doi: 10.1016/j.ultrasmedbio.2015.10.001
Legon, W., Ai, L., Bansal, P., and Mueller, J. K. (2018a). Neuromodulation with single-element transcranial focused ultrasound in human thalamus. Hum. Brain Mapp. 39, 1995–2006. doi: 10.1002/hbm.23981
Legon, W., Bansal, P., Tyshynsky, R., Ai, L., and Mueller, J. K. (2018b). Transcranial focused ultrasound neuromodulation of the human primary motor cortex. Sci. Rep. 8:10007. doi: 10.1038/s41598-018-28320-1
Legon, W., Sato, T. F., Opitz, A., Mueller, J., Barbour, A., Williams, A., et al. (2014). Transcranial focused ultrasound modulates the activity of primary somatosensory cortex in humans. Nat. Publishing Group 17, 322–329. doi: 10.1038/nn.3620
Leung, S. A., Webb, T. D., Bitton, R. R., Ghanouni, P., and Butts Pauly, K. (2019). A rapid beam simulation framework for transcranial focused ultrasound. Sci. Rep. 9:7965. doi: 10.1038/s41598-019-43775-6
Li, G.-F., Zhao, H.-X., Zhou, H., Yan, F., Wang, J.-Y., Xu, C.-X., et al. (2016). Improved anatomical specificity of non-invasive neuro-stimulation by high frequency (5MHz) ultrasound. Sci. Rep. 6:24738. doi: 10.1038/srep24738
Magnin, R., Rabusseau, F., Salabartan, F., Mériaux, S., Aubry, J.-F., Le Bihan, D., et al. (2015). Magnetic resonance-guided motorized transcranial ultrasound system for blood- brain barrier permeabilization along arbitrary trajectories in rodents. J. Ther. Ultrasound 3:22. doi: 10.1186/s40349-015-0044-5
Maimbourg, G., Houdouin, A., Deffieux, T., Tanter, M., and Aubry, J.-F. (2018a). 3D-printed adaptive acoustic lens as a disruptive technology for transcranial ultrasound therapy using single-element transducers. Phys. Med. Biol. 63:025026. doi: 10.1088/1361-6560/aaa037
Maimbourg, G., Houdouin, A., Deffieux, T., Tanter, M., and Aubry, J.-F. (2019). Steering capabilities of an acoustic lens for transcranial therapy: numerical and experimental studies. IEEE Trans. Biomed. Eng. 67, 27–37. doi: 10.1109/TBME.2019.2907556
Maimbourg, G., Houdouin, A., Santin, M., Lehericy, S., Tanter, M., and Aubry, J.-F. (2018b). Inside/outside the brain binary cavitation localization based on the lowpass filter effect of the skull on the harmonic content: a proof of concept study. Phys. Med. Biol. 63, 135012–135018. doi: 10.1088/1361-6560/aaca21
Marquet, F., Pernot, M., Aubry, J. F., Montaldo, G., Marsac, L., Tanter, M., et al. (2009). Non-invasive transcranial ultrasound therapy based on a 3D CT scan: protocol validation and in vitroresults. Phys. Med. Biol. 54, 2597–2613. doi: 10.1088/0031-9155/54/9/001
McDannold, N., Clement, G. T., Black, P., Jolesz, F., and Hynynen, K. (2010). Transcranial magnetic resonance imaging– guided focused ultrasound surgery of brain tumors. Neurosurgery 66, 323–332. doi: 10.1227/01.NEU.0000360379.95800.2F
McDannold, N., Zhang, Y., Power, C., Arvanitis, C. D., Vykhodtseva, N., and Livingstone, M. (2015). Targeted, noninvasive blockade of cortical neuronal activity. Sci. Rep. 5:16253. doi: 10.1038/srep16253
Mehić, E., Xu, J. M., Caler, C. J., Coulson, N. K., Moritz, C. T., and Mourad, P. D. (2014). Increased anatomical specificity of neuromodulation via modulated focused ultrasound. PLoS One 9:e86939. doi: 10.1371/journal.pone.0086939
Miller, G. W., Eames, M., Snell, J., and Aubry, J.-F. (2015). Ultrashort echo-time MRI versus CT for skull aberration correction in MR-guided transcranial focused ultrasound: in vitrocomparison on human calvaria. Med. Phys. 42, 2223–2233. doi: 10.1118/1.4916656
Mohammadjavadi, M., Ye, P. P., Xia, A., Brown, J., Popelka, G., and Pauly, K. B. (2019). Elimination of peripheral auditory pathway activation does not affect motor responses from ultrasound neuromodulation. Brain Stimul. 12, 901–910. doi: 10.1016/j.brs.2019.03.005
Mueller, J., Legon, W., Opitz, A., Sato, T. F., and Tyler, W. J. (2014). Transcranial focused ultrasound modulates intrinsic and evoked EEG dynamics. Brain Stimul. 7, 900–908. doi: 10.1016/j.brs.2014.08.008
Neggers, S. F. W., Langerak, T. R., Schutter, D. J. L. G., Mandl, R. C. W., Ramsey, N. F., Lemmens, P. J. J., et al. (2004). A stereotactic method for image-guided transcranial magnetic stimulation validated with fMRI and motor-evoked potentials. Neuroimage 21, 1805–1817. doi: 10.1016/j.neuroimage.2003.12.006
Okamoto, M., Dan, H., Sakamoto, K., Takeo, K., Shimizu, K., Kohno, S., et al. (2004). Three-dimensional probabilistic anatomical cranio-cerebral correlation via the international 10–20 system oriented for transcranial functional brain mapping. Neuroimage 21, 99–111. doi: 10.1016/j.neuroimage.2003.08.026
Olk, B., Chang, E., Kingstone, A., and Ro, T. (2005). Modulation of antisaccades by transcranial magnetic stimulation of the human frontal eye field. Cereb. Cortex 16, 76–82. doi: 10.1093/cercor/bhi085
Ozenne, V., Constans, C., Bour, P., Santin, M. D., Valabrègue, R., Ahnine, H., et al. (2019). MRI monitoring of temperature and displacement for transcranial focus ultrasound applications. Neuroimage 204:116236. doi: 10.1016/j.neuroimage.2019.116236
Parthasarathy, H. B., Schall, J. D., and Graybiel, A. M. (1992). Distributed but convergent ordering of corticostriatal projections: analysis of the frontal eye field and the supplementary eye field in the macaque monkey. J. Neurosci. 12, 4468–4488. doi: 10.1523/JNEUROSCI.12-11-04468.1992
Paus, T., Jech, R., Thompson, C. J., Comeau, R., Peters, T., and Evans, A. C. (1997). Transcranial magnetic stimulation during positron emission tomography: a new method for studying connectivity of the human cerebral cortex. J. Neurosci. 17, 3178–3184. doi: 10.1523/JNEUROSCI.17-09-03178.1997
Pierrot-Deseilligny, C., Rivaud, S., Gaymard, B., and Agid, Y. (1991). Cortical control of reflexive visually-guided saccades. Brain 114(Pt 3), 1473–1485. doi: 10.1093/brain/114.3.1473
Pinton, G., Aubry, J.-F., Bossy, E., Muller, M., Pernot, M., and Tanter, M. (2011). Attenuation, scattering, and absorption of ultrasound in the skull bone. Med. Phys. 39, 299–307. doi: 10.1118/1.3668316
Pinton, G., Aubry, J.-F., Fink, M., and Tanter, M. (2012). Numerical prediction of frequency dependent 3D maps of mechanical index thresholds in ultrasonic brain therapy. Med. Phys. 39, 455–467. doi: 10.1118/1.3670376
Plaksin, M., Shoham, S., and Kimmel, E. (2014). Intramembrane cavitation as a predictive bio-piezoelectric mechanism for ultrasonic brain stimulation. Phys. Rev. X 4:011004. doi: 10.1103/PhysRevX.4.011004
Pouliopoulos, A. N., Wu, S.-Y., Burgess, M. T., Karakatsani, M. E., Kamimura, H. A. S., and Konofagou, E. E. (2019). A clinical system for non-invasive blood–brain barrier opening using a neuronavigation-guided single-element focused ultrasound transducer. Ultrasound Med. Biol. 46, 73–89. doi: 10.1016/j.ultrasmedbio.2019.09.010
Pulkkinen, A., Huang, Y., Song, J., and Hynynen, K. (2011). Simulations and measurements of transcranial low-frequency ultrasound therapy: skull-base heating and effective area of treatment. Phys. Med. Biol. 56, 4661–4683. doi: 10.1088/0031-9155/56/15/003
Pulkkinen, A., Werner, B., Martin, E., and Hynynen, K. (2014). Numerical simulations of clinical focused ultrasound functional neurosurgery. Phys. Med. Biol. 59, 1679–1700. doi: 10.1088/0031-9155/59/7/1679
Rossi, S., and Rossini, P. M. (2004). TMS in cognitive plasticity and the potential for rehabilitation. Trends Cogn. Sci. 8, 273–279. doi: 10.1016/j.tics.2004.04.012
Rossini, P. M., Burke, D., Chen, R., Cohen, L. G., Daskalakis, Z., Di Iorio, R., et al. (2015). Non-invasive electrical and magnetic stimulation of the brain, spinal cord, roots and peripheral nerves: basic principles and procedures for routine clinical and research application. An updated report from an I.F.C.N. Committee. Clin. Neurophysiol. 126, 1071–1107. doi: 10.1016/j.clinph.2015.02.001
Royer, D., and Casula, O. (1994). “Quantitative imaging of transient acoustic fields by optical heterodyne interferometry,” in Proceedings of IEEE Ultrasonics Symposium 2, Cannes, 1153–1162.
Royer, D., and Dieulesaint, E. (1986). Optical probing of the mechanical impulse response of a transducer. Appl. Phys. Lett. 49, 1056–1058. doi: 10.1063/1.97473
Ruohonen, J., and Karhu, J. (2010). Navigated transcranial magnetic stimulation. Neurophysiol. Clin. Clin. Neurophysiol. 40, 7–17. doi: 10.1016/j.neucli.2010.01.006
Sanguinetti, J. L., Hameroff, S., Smith, E. E., Sato, T., Daft, C. M. W., Tyler, W. J., et al. (2020). Transcranial focused ultrasound to the right prefrontal cortex improves mood and alters functional connectivity in humans. Front. Hum. Neurosci. 14:52. doi: 10.3389/fnhum.2020.00052
Schall, J. D. (2004). On the role of frontal eye field in guiding attention and saccades. Vis. Res. 44, 1453–1467. doi: 10.1016/j.visres.2003.10.025
Schall, J. D. (2015). Visuomotor functions in the frontal lobe. Annu. Rev. Vis. Sci. 1, 469–498. doi: 10.1146/annurev-vision-082114-035317
Schall, J. D., Morel, A., and Kaas, J. H. (2009). Topography of supplementary eye field afferents to frontal eye field in macaque: implications for mapping between saccade coordinate systems. Vis. Neurosci. 10, 385–393. doi: 10.1017/S0952523800003771
Schlag-Rey, M., Amador, N., Sanchez, H., and Schlag, J. (1997). Antisaccade performance predicted by neuronal activity in the supplementary eye field. Nature 390, 398–401. doi: 10.1038/37114
Schoenbaum, G., Setlow, B., Saddoris, M. P., and Gallagher, M. (2003). Encoding predicted outcome and acquired value in orbitofrontal cortex during cue sampling depends upon input from basolateral amygdala. Neuron 39, 855–867. doi: 10.1016/S0896-6273(03)00474-4
Schönfeldt-Lecuona, C., Thielscher, A., Freudenmann, R. W., Kron, M., Spitzer, M., and Herwig, U. (2005). Accuracy of stereotaxic positioning of transcranial magnetic stimulation. Brain Topogr. 17, 253–259. doi: 10.1007/s10548-005-6033-1
Siddiqi, S. H., Taylor, S. F., Cooke, D., Pascual-Leone, A., George, M. S., and Fox, M. D. (2020). Distinct symptom-specific treatment targets for circuit-based neuromodulation. Am. J. Psychiatry 177, 435–446. doi: 10.1176/appi.ajp.2019.19090915
Siebner, H. R., Bergmann, T. O., Bestmann, S., Massimini, M., Johansen-Berg, H., Mochizuki, H., et al. (2009). Consensus paper: combining transcranial stimulation with neuroimaging. Brain Stimul. 2, 58–80. doi: 10.1016/j.brs.2008.11.002
Sparing, R., Buelte, D., Meister, I. G., Paus, T., and Fink, G. R. (2007). Transcranial magnetic stimulation and the challenge of coil placement: a comparison of conventional and stereotaxic neuronavigational strategies. Hum. Brain Mapp. 29, 82–96. doi: 10.1002/hbm.20360
Thut, G., and Pascual-Leone, A. (2009). A review of combined TMS-EEG studies to characterize lasting effects of repetitive TMS and assess their usefulness in cognitive and clinical neuroscience. Brain Topogr. 22, 219–232. doi: 10.1007/s10548-009-0115-4
Verhagen, L., Gallea, C., Folloni, D., Constans, C., Jensen, D. E., Ahnine, H., et al. (2019). Offline impact of transcranial focused ultrasound on cortical activation in primates. eLife 8:875. doi: 10.7554/eLife.40541
Wang, H., Kim, M., Normoyle, K. P., and Llano, D. (2016). Thermal regulation of the brain—an anatomical and physiological review for clinical neuroscientists. Front. Neurosci. 9:528. doi: 10.3389/fnins.2015.00528
Wassermann, E. M., and Lisanby, S. H. (2001). Therapeutic application of repetitive transcranial magnetic stimulation: a review. Clin. Neurophysiol. 112, 1367–1377. doi: 10.1016/s1388-2457(01)00585-5
Wattiez, N., Constans, C., Deffieux, T., Daye, P. M., Tanter, M., Aubry, J.-F., et al. (2017). Transcranial ultrasonic stimulation modulates single-neuron discharge in macaques performing an antisaccade task. Brain Stimul. 10, 1024–1031. doi: 10.1016/j.brs.2017.07.007
Wintermark, M., Tustison, N. J., Elias, W. J., Patrie, J. T., Xin, W., Demartini, N., et al. (2014). T1-weighted MRI as a substitute to CT for refocusing planning in MR-guided focused ultrasound. Phys. Med. Biol. 59, 3599–3614. doi: 10.1088/0031-9155/59/13/3599
Yang, P. S., Kim, H., Lee, W., Bohlke, M., Park, S., Maher, T. J., et al. (2012). Transcranial focused ultrasound to the thalamus is associated with reduced extracellular GABA levels in rats. Neuropsychobiology 65, 153–160. doi: 10.1159/000336001
Ye, P. P., Brown, J. R., and Pauly, K. B. (2016). Frequency dependence of ultrasound neurostimulation in the mouse brain. Ultrasound Med. Biol. 42, 1512–1530. doi: 10.1016/j.ultrasmedbio.2016.02.012
Yin, X., and Hynynen, K. (2005). A numerical study of transcranial focused ultrasound beam propagation at low frequency. Phys. Med. Biol. 50, 1821–1836. doi: 10.1088/0031-9155/50/8/013
Yoo, S.-S., Bystritsky, A., Lee, J.-H., Zhang, Y., Fischer, K., Min, B.-K., et al. (2011). Focused ultrasound modulates region-specific brain activity. Neuroimage 56, 1267–1275. doi: 10.1016/j.neuroimage.2011.02.058
Younan, Y., Deffieux, T., Larrat, B., Fink, M., Tanter, M., and Aubry, J.-F. (2013). Influence of the pressure field distribution in transcranial ultrasonic neurostimulation. Med. Phys. 40, 82902–82912. doi: 10.1118/1.4812423
Keywords: focused ultrasound, neurostimulation, neuronavigation, brain, transcranial, non-human primates
Citation: Pouget P, Frey S, Ahnine H, Attali D, Claron J, Constans C, Aubry J-F and Arcizet F (2020) Neuronavigated Repetitive Transcranial Ultrasound Stimulation Induces Long-Lasting and Reversible Effects on Oculomotor Performance in Non-human Primates. Front. Physiol. 11:1042. doi: 10.3389/fphys.2020.01042
Received: 15 April 2020; Accepted: 29 July 2020;
Published: 19 August 2020.
Edited by:Hermes Arytto Salles Kamimura, Columbia University, United States
Reviewed by:Suliann Ben Hamed, UMR 5229 Institut des Sciences Cognitives Marc Jeannerod, France
Hyungmin Kim, Korea Institute of Science and Technology (KIST), South Korea
Copyright © 2020 Pouget, Frey, Ahnine, Attali, Claron, Constans, Aubry and Arcizet. 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: Fabrice Arcizet, email@example.com; firstname.lastname@example.org
†These authors have contributed equally to this work