Effects of sonication parameters on transcranial focused ultrasound brain stimulation in an ovine model
Kyungho Yoon, Wonhye Lee, Ji Eun Lee, Linda Xu, Phillip Croce, Lori Foley, Seung-Schik Yoo
PLOS ONE 2019, 14, e0224311 · 10.1371/journal.pone.0224311
Abstract
Low-intensity focused ultrasound (FUS) has significant potential as a non-invasive brain stimulation modality and novel technique for functional brain mapping, particularly with its advantage of greater spatial selectivity and depth penetration compared to existing non-invasive brain stimulation techniques. As previous studies, primarily carried out in small animals, have demonstrated that sonication parameters affect the stimulation efficiency, further investigation in large animals is necessary to translate this technique into clinical practice. In the present study, we examined the effects of sonication parameters on the transient modification of excitability of cortical and thalamic areas in an ovine model. Guided by anatomical and functional neuroimaging data specific to each animal, 250 kHz FUS was transcranially applied to the primary sensorimotor area associated with the right hind limb and its thalamic projection in sheep (n = 10) across multiple sessions using various combinations of sonication parameters. The degree of effect from FUS was assessed through electrophysiological responses, through analysis of electromyogram and electroencephalographic somatosensory evoked potentials for evaluation of excitatory and suppressive effects, respectively. We found that the modulatory effects were transient and reversible, with specific sonication parameters outperforming others in modulating regional brain activity. Magnetic resonance imaging and histological analysis conducted at different time points after the final sonication session, as well as behavioral observations, showed that repeated exposure to FUS did not damage the underlying brain tissue. Our results suggest that FUS-mediated, non-invasive, region-specific bimodal neuromodulation can be safely achieved in an ovine model, indicating its potential for translation into human studies.
Abstract via europepmc.
Exposures
Exposure 1: M1 excitatory sonication (parameter sets EP1-EP30)
Target: primary motor cortex — “M1 (primary motor area, right hind limb representation)”
Device: Ultran · Ultran Group · GPS200-400128 ✓
| Waveform | pulsed, continuous | |
|---|---|---|
| Fundamental frequency (kHz) | 250 | ✓✓✓ |
| Pulse duration (ms) | 0.5, 1, 2, 3swept | ✓✓✓ |
| Pulse repetition frequency (Hz) | 600, 300, 150, 100, 1,000, 500, 250, 167, 1,400, 700, 350, 233swept | ✓✓✓ |
| Duty cycle (%) | 30, 50, 70, 100swept | ✓✓✓ |
| Sonication duration (s) | 0.06, 0.1, 0.14, 0.2swept | ✓✓✓ |
| Free-field pressure (kPa) | not reported | |
|---|---|---|
| Free-field Isppa (W/cm²) | not reported | |
| Free-field Ispta (W/cm²) | not reported | |
| In-situ estimate | deratingsingle value | |
| In-situ pressure (kPa) | not reported | |
| In-situ Isppa (W/cm²) | 15.8, 18.2swept | ✓✓✓ |
| In-situ Ispta (W/cm²) | 4.7, 5.5, 7.9, 9.1, 11.1, 12.7, 15.8, 18.2swept | ✓✓✓ |
In more detail, for the excitatory sonication in pulsed mode, we used a sonication duration of 200 ms. We combined four tone-burst durations (0.5, 1, 2, and 3 ms), three DCs (30, 50, and 70%, and corresponding PRFs), and two levels of in situ Isppa (15.8 and 18.2 W/cm2) for a total of 24 parameter sets (named EP1-EP24, see Table 1). Six additional sets (EP25-EP30), operating under continuous mode sonication (100% DC), were also administered with three sonication durations (60, 100, and 140 ms). Sonication in stimulation trials was given every 5 s (ISI = 5 s). The maximum response rate (12.4 +/- 9.2%) was observed with EP17 (0.5 ms tone-burst duration, 70% DC, 15.8 W/cm2 Isppa) stimulating M1.
Exposure 2: S1 suppressive sonication (parameter sets SP1-SP8)
Target: primary somatosensory cortex — “S1 (primary somatosensory area, right hind limb representation)”
Device: Ultran · Ultran Group · GPS200-400128 ✓
| Waveform | pulsed | |
|---|---|---|
| Fundamental frequency (kHz) | 250 | ✓✓✓ |
| Pulse duration (ms) | 0.5, 1swept | ✓✓✓ |
| Pulse repetition frequency (Hz) | 100, 60, 50, 30swept | ✓✓✓ |
| Duty cycle (%) | 3, 5swept | ✓✓✓ |
| Sonication duration (s) | 120 | ✓✓✓ |
| Free-field pressure (kPa) | not reported | |
|---|---|---|
| Free-field Isppa (W/cm²) | not reported | |
| Free-field Ispta (W/cm²) | not reported | |
| In-situ estimate | deratingsingle value | |
| In-situ pressure (kPa) | not reported | |
| In-situ Isppa (W/cm²) | 5.4, 11.6swept | ✓✓✓ |
| In-situ Ispta (W/cm²) | 0.16, 0.27, 0.35, 0.58swept | ✓✓✓ |
To evaluate the suppressive FUS parameters, eight parameter sets of varying tone-burst durations (0.5 and 1.0 ms), DCs (3 and 5%), and levels of in situ Isppa (5.4 and 11.6 W/cm2) were examined (named SP1-SP8, see Table 2). Each suppressive sonication set was given for a duration of 2 min once per session (therefore, ISI did not apply). The use of SP1 and SP3 (tone-burst duration of 0.5 ms and Isppa of 5.4 W/cm2, acquired at 3 and 5% DCs) temporarily suppressed the SEP (18-35% reduction).
Exposure 3: Thalamus (VL/VPL): excitatory sonication (EP1-EP30) and suppressive sonication (SP1-SP8)
Target: ventral lateral nucleus, ventral posterolateral nucleus — “thalamus (ventrolateral nucleus VL / ventral posterolateral nucleus VPL, treated as a single sonication target)”
Device: Ultran · Ultran Group · GPS200-400128 ✓
| Waveform | pulsed, continuous | |
|---|---|---|
| Fundamental frequency (kHz) | 250 | ✓✓✓ |
| Pulse duration (ms) | 0.5, 1, 2, 3, 1swept | ✓✓✓ |
| Pulse repetition frequency (Hz) | 600, 300, 150, 100, 1,000, 500, 250, 167, 1,400, 700, 350, 233, 60, 50, 30swept | ✓✓✓ |
| Duty cycle (%) | 30, 50, 70, 100, 3, 5swept | ✓✓✓ |
| Sonication duration (s) | 0.06, 0.1, 0.14, 0.2, 120swept | ✓✓✓ |
| Free-field pressure (kPa) | not reported | |
|---|---|---|
| Free-field Isppa (W/cm²) | not reported | |
| Free-field Ispta (W/cm²) | not reported | |
| In-situ estimate | deratingsingle value | |
| In-situ pressure (kPa) | not reported | |
| In-situ Isppa (W/cm²) | 15.8, 18.2, 5.4, 11.6swept | ✓✓✓ |
| In-situ Ispta (W/cm²) | 0.16, 0.27, 0.35, 0.58, 4.7, 5.5, 7.9, 9.1, 11.1, 12.7, 15.8, 18.2swept | ✓✓✓ |
The proximity of the VL and VPL in ovine neuroanatomy was beyond the spatial accuracy of sonication; therefore, a single spatial coordinate in the thalamus was assigned as a sonication target for both the excitatory parameter sets (EP1-EP30, same as used for M1) and the suppressive parameter sets (SP1-SP8, same as used for S1). Thalamic stimulation with excitatory parameters was highest with EP18, EP19, and EP20 (maximum response rate 6.9%); suppressive sonication of the thalamus with SP1 and SP3 produced a 31.4-34.3% reduction in SEP magnitude.
Flags from extraction
exposures[2]— Per the 'one exposure per target x frequency' rule, the thalamus exposure combines two very different protocols that share only the target (thalamus) and frequency (250 kHz): a short (200 ms/60-140 ms), high-duty-cycle (30-100%) excitatory protocol (EP1-EP30) and a long (2 min), low-duty-cycle (3-5%) suppressive protocol (SP1-SP8). The resulting list-valued timing/intensity fields span both protocols and should not be read as paired combinations; see protocol_description and exp_m1/exp_s1 for the disaggregated values.exposures[2].target— The paper treats a single spatial coordinate in the thalamus as the sonication target because the proximity of VL (motor) and VPL (sensory) nuclei was below the spatial accuracy of sonication/registration; both target ids (ventral_lateral_nucleus, ventral_posterolateral_nucleus) are listed since the intended targets for excitatory vs suppressive sessions differed (VL vs VPL) even though a single physical coordinate was used.*.free_field— The paper reports only 'in situ' (derated) Isppa/Ispta/MI values in Tables 1-2; the pre-derating (free-field, degassed-water) intensities are not given as numbers in the text, only the attenuation percentages (16.3% coupling gel, 69.4% skull) used to derive the in-situ values, so free_field fields are left not_reported.n_sessions_per_subject— 66 total sonication sessions were conducted across 10 sheep with 'up to four' sessions per animal and variable numbers of animals per parameter set (see Tables 1-2); not reduced to a single number per subject.sham_type— For the suppressive (SEP) experiments, an active control condition was used: 'we also delivered sonication to the S1 and thalamic area in the right hemisphere' (functionally non-target, contralateral-pathway hemisphere) using SP3/SP7 parameters, classified as active_control_site. No sham/control condition is described for the excitatory (EMG) experiments.auditory_control— Discussion notes that audible buzzing sound was sometimes produced by the coupling hydrogel during excitatory (but not suppressive) sonication, raising a possible auditory confound, but no explicit auditory masking or control procedure is reported.exposures[2] (thalamus)— Source quotes for this combined exposure are keyed primarily to the excitatory Table 1 (EP1-EP30) HTML row content; the suppressive Table 2 (SP1-SP8) content, quoted separately under exposures[1], additionally contributes the duty_cycle_pct=[3,5], pulse_repetition_frequency_hz=[30,50,60,100] and sonication_duration_s=120 s values in this exposure's lists but is not re-keyed here because only one source quote is stored per field path.
Notes: Sample sizes and number of parameter sets per animal vary considerably by target and parameter set (see Tables 1-2 'Number of animals' columns, ranging 8-10); a single n_subjects=10 (total sheep enrolled) is reported at the top level. Frequency (250 kHz) and transducer are constant across all three exposures. Isppa/Ispta values in Tables 1 and 2 are already the paper's derated 'in situ' estimates (intensity transmission level of 25.6% applied to free-field/water-calibrated values), hence in_situ.method='derating'.