Skip to main content

Main menu

  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Current Issue
    • Issue Archive
    • Collections
    • Podcast
  • ALERTS
  • FOR AUTHORS
    • Information for Authors
    • Fees
    • Journal Clubs
    • eLetters
    • Submit
  • EDITORIAL BOARD
  • ABOUT
    • Overview
    • Advertise
    • For the Media
    • Rights and Permissions
    • Privacy Policy
    • Feedback
  • SUBSCRIBE

User menu

  • Log in
  • My Cart

Search

  • Advanced search
Journal of Neuroscience
  • Log in
  • My Cart
Journal of Neuroscience

Advanced Search

Submit a Manuscript
  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Current Issue
    • Issue Archive
    • Collections
    • Podcast
  • ALERTS
  • FOR AUTHORS
    • Information for Authors
    • Fees
    • Journal Clubs
    • eLetters
    • Submit
  • EDITORIAL BOARD
  • ABOUT
    • Overview
    • Advertise
    • For the Media
    • Rights and Permissions
    • Privacy Policy
    • Feedback
  • SUBSCRIBE
PreviousNext
Featured ArticleResearch Articles, Cellular/Molecular

Ultrasound Elicits Behavioral Responses through Mechanical Effects on Neurons and Ion Channels in a Simple Nervous System

Jan Kubanek, Poojan Shukla, Alakananda Das, Stephen A. Baccus and Miriam B. Goodman
Journal of Neuroscience 21 March 2018, 38 (12) 3081-3091; DOI: https://doi.org/10.1523/JNEUROSCI.1458-17.2018
Jan Kubanek
1Department of Molecular and Cellular Physiology and
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Jan Kubanek
Poojan Shukla
1Department of Molecular and Cellular Physiology and
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Alakananda Das
1Department of Molecular and Cellular Physiology and
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Alakananda Das
Stephen A. Baccus
2Department of Neurobiology, Stanford, California 94305
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Stephen A. Baccus
Miriam B. Goodman
1Department of Molecular and Cellular Physiology and
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Miriam B. Goodman
  • Article
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF
Loading

Abstract

Focused ultrasound has been shown to stimulate excitable cells, but the biophysical mechanisms behind this phenomenon remain poorly understood. To provide additional insight, we devised a behavioral-genetic assay applied to the well-characterized nervous system of Caenorhabditis elegans nematodes. We found that pulsed ultrasound elicits robust reversal behavior in wild-type animals in a pressure-, duration-, and pulse protocol-dependent manner. Responses were preserved in mutants unable to sense thermal fluctuations and absent in mutants lacking neurons required for mechanosensation. Additionally, we found that the worm's response to ultrasound pulses rests on the expression of MEC-4, a DEG/ENaC/ASIC ion channel required for touch sensation. Consistent with prior studies of MEC-4-dependent currents in vivo, the worm's response was optimal for pulses repeated 300–1000 times per second. Based on these findings, we conclude that mechanical, rather than thermal, stimulation accounts for behavioral responses. Further, we propose that acoustic radiation force governs the response to ultrasound in a manner that depends on the touch receptor neurons and MEC-4-dependent ion channels. Our findings illuminate a complete pathway of ultrasound action, from the forces generated by propagating ultrasound to an activation of a specific ion channel. The findings further highlight the importance of optimizing ultrasound pulsing protocols when stimulating neurons via ion channels with mechanosensitive properties.

SIGNIFICANCE STATEMENT How ultrasound influences neurons and other excitable cells has remained a mystery for decades. Although it is widely understood that ultrasound can heat tissues and induce mechanical strain, whether or not neuronal activation depends on heat, mechanical force, or both physical factors is not known. We harnessed Caenorhabditis elegans nematodes and their extraordinary sensitivity to thermal and mechanical stimuli to address this question. Whereas thermosensory mutants respond to ultrasound similar to wild-type animals, mechanosensory mutants were insensitive to ultrasound stimulation. Additionally, stimulus parameters that accentuate mechanical effects were more effective than those producing more heat. These findings highlight a mechanical nature of the effect of ultrasound on neurons and suggest specific ways to optimize stimulation protocols in specific tissues.

  • Caenorhabditis elegans
  • mechanosensation
  • mechanosensitive ion channels
  • neurostimulation
  • thermosensation
  • ultrasound

Introduction

Low-intensity, focused ultrasound affects the function of excitable cells in the central (Fry et al., 1958; Meyers et al., 1959; Foster and Wiederhold, 1978; Gavrilov et al., 1996; Tufail et al., 2011; Yoo et al., 2011; Deffieux et al., 2013; King et al., 2013; Menz et al., 2013) and peripheral (Mihran et al., 1990; Tsui et al., 2005; Colucci et al., 2009) nervous systems, and the heart (Harvey, 1929; Buiochi et al., 2012). Because it propagates deep into tissue while retaining spatial focus, it has garnered considerable attention for its potential as a noninvasive tool for stimulation of the brain and the heart (Tufail et al., 2011; Lee et al., 2016). Despite these investigations, how ultrasound stimulates excitable cells has been a mystery since the discovery of its stimulatory effects >8 decades ago (Harvey, 1929).

Physical mechanisms for ultrasound-dependent tissue excitation have been divided into thermal and nonthermal effects (Dalecki, 2004; Naor et al., 2016; Sassaroli and Vykhodtseva, 2016; Ye et al., 2016). The latter category encompasses mechanical effects, such as membrane oscillation, cavitation, or radiation force. Ultrasound stimulation generally occurs under conditions that heat tissues by <1°C (Tufail et al., 2010; Yoo et al., 2011; Menz et al., 2013). Nonetheless, even small temperature changes could activate certain classes of ion channels, such as temperature-sensitive TRP cation channels (Diaz-Franulic et al., 2016) and TREK potassium channels (Schneider et al., 2014). Moreover, the rate of temperature change may also contribute to the stimulatory effects (Rabbitt et al., 2016).

The mechanical effects of ultrasound could be converted into ionic currents and changes in electrical excitability by increasing mechanical strain in a manner that directly activates ion channels (Tyler, 2011) or by changes in membrane thickness and capacitance (Plaksin et al., 2014). For instance, increases in membrane tension are thought to catalyze conformational change and activate ion channels in which the open state has a larger cross-sectional area than the closed state (Sukharev and Corey, 2004). Thus, if ultrasound stimulation were to change membrane tension, then it could alter the activity of mechanosensitive channels, including voltage-gated sodium and potassium channels that exhibit membrane tension-dependent gating (Beyder et al., 2010; Brohawn et al., 2014). Consistent with this idea, ultrasound stimulation has been shown to increase currents carried by two-pore domain (K2P) potassium channels and voltage-gated sodium channels expressed in Xenopus oocytes (Kubanek et al., 2016).

This body of work has raised the question whether or not ultrasound affects nervous systems via thermal or nonthermal effects (Iversen et al., 2017) and, if nonthermal effects are dominant, how is ultrasound transduced into membrane strain (tension) and channel activation. To contribute additional insight into these questions, we developed a behavioral-genetic assay using Caenorhabditis elegans roundworms. The nervous system of this animal has been characterized in its entirety, the animal performs simple behaviors that are easily monitored using video tracking and machine vision tools (Husson et al., 2005), and the model is amenable to experimental studies involving tens to hundreds of individual animals. Two additional aspects of the worm's biology are exploited in this study. First, C. elegans has extraordinarily sensitive thermosensory and mechanosensory neurons able to detect thermal fluctuations of ≤0.05°C (Ramot et al., 2008) and applied forces of 50 nN (O'Hagan et al., 2005), respectively. Second, mutants exist in which thermosensory and mechanosensory neurons are disabled independently.

We found that pulsed ultrasound stimulation evokes avoidance responses whose probability increases with acoustic pressure and stimulus duration and shows optima for specific ultrasound pulsing protocols. Mutants lacking neurons required for thermosensation responded similar to wild-type animals, whereas those lacking neurons or ion channels required for mechosensation failed to mount avoidance responses to ultrasound stimulation. The mechanical nature of the effect led us to a detailed characterization of the involved neurons and ion channels as well as a characterization of the optimal stimulation parameters.

Materials and Methods

Animals and strains.

The following strains were used in this study: N2 (RRID:WB-STRAIN:N2(ancestral));CB1338 mec-3(e1338) IV (RRID:WB-STRAIN:CB1338); CB1611 mec-4(e1611) X (RRID:WB-STRAIN:CB1611); TU253 mec-4(u253) X (RRID:WB-STRAIN:TU253); IK597 gcy-23(nj37)gcy-8(oy44)gcy-18(nj38) IV (RRID:WB-STRAIN:IK597); VC1141 trp-4(ok1605) I (RRID:WB-STRAIN:VC1141); VC818 trp-4(gk341) I (RRID:WB-STRAIN:VC818); TQ296 trp-4(sy695) I (RRID:WB-STRAIN:TQ296); GN716 trp-4(ok1605) I, outcrossed four times from VC1141. All mutants were derived from the N2 (Bristol) background, which serves as the wild-type strain in this study. Except for GN716, which was prepared explicitly for this study, and TQ296, which was a gift of X.Z.S. Xu (University of Michigan), strains were obtained either from a frozen repository maintained in our laboratory or from the Caenorhabditis Genetics Center.

The three trp-4 alleles we studied all encode deletion or null alleles of the trp-4 gene, which encodes the key pore-forming subunit of a mechanosensory ion channel expressed in the CEP texture-sensing neurons (Kang et al., 2010). TRP-4 is orthologous to the mechanically-gated NOMPC channel from Drosophila.

The ok1605 allele encodes an in-frame, 1 kb deletion that removes exons 12–14 of the trp-4 gene and is predicted to result in the loss of ankyrin repeats 16–21 from the TRP-4 protein. The gk341 allele contains a small deletion encompassing exon 2 and is predicted to cause a frameshift in the transcript and the introduction of an early stop codon. The sy695 allele contains an unmapped 3 kb deletion in the 3′ region of the gene. This deletion is thought to disrupt the transmembrane pore-forming domain of TRP-4 (Li et al., 2006). The GN716 trp-4(ok1605) strain was derived by outcrossing VC1141 trp-4(ok1605) with wild-type (N2) in four rounds. We used PCR to verify that all trp-4 mutant strains harbored the expected genetic deletions in the trp-4 locus.

Imaging and transducer control.

For each assay, we transferred a single adult animal from a growth plate to a 4-mm-thick NGM agar slab that was free of bacteria. Because the agar slab consists of a 2% agar solution in saline, it is a good approximation of the acoustic properties of many biological tissues (Altman et al., 1974). To retain animals within the camera's field of view, we created a boundary consisting of a filter paper ring saturated by a copper sulfate (500 mm) solution, as described previously (Ramot et al., 2008).

To deliver ultrasound stimuli, we used a commercially available piezoelectric ultrasonic transducer (A327S-SU-CF1.00IN-PTF, Olympus, 1-inch line-focused) positioned 1 inch (2.54 cm) below the top of the agar slab and oriented perpendicular to the surface of the agar slab (see Fig. 1A). We filled the space between the transducer surface and the bottom of the agar slab with degassed water, contained within a plastic cone mounted on the transducer. The water was degassed by boiling for 30 min and stored in air-tight tubes. The slab did not appear to attenuate ultrasound pressure, according to measurements with a hydrophone (data not shown).

We used oblique illumination via a circular array (20 cm in diameter) of red LEDs to provide the optical contrast between animals and the surface of the agar slab needed to track animal movement using the Parallel Worm Tracker (Ramot et al., 2008). The image was magnified 3× (zoom lens, Navitar). The camera's chip FOV was 5.6 × 4.2 mm. Image contrast was optimal when the plane of the LED array was ∼1 cm above the top of the agar slab. We also used a blue LED, controlled by an Arduino Uno board and mounted 5 cm above the agar slab, to deliver an optical signal synchronized to the stimulus onset.

To generate signals driving the ultrasound transducer, we used a function generator (HP 8116A, Hewlett-Packard) controlled by an Arduino Uno board and an amplifier with a gain of 50 dB (ENI-240L, ENI). The acoustic pressures we generated were measured in free field using a calibrated hydrophone (HGL-0200, Onda) combined with a preamplifier (AG-2020, Onda). The hydrophone measurements were performed at the peak spatial pressure. The hydrophone manufacturer's calibration values around the frequency of 10 MHz were steady and showed only minimal level of noise.

Behavioral recordings.

For each trial, a freely moving animal was monitored via a digital video camera (SME-B050-U, Mightex) operated in a live-video mode until it approached the ultrasound focus head first. We started recording videos ∼5 s before the predicted approach to the ultrasound focus and continued recording for ∼10 s following the delivery of each stimulus. Each individual animal was tested in 10 trials with an intertrial interval of at least 20 s. All animals were hermaphrodites, assayed blind to genotype and as adults. We applied stimuli at the following pressures (in MPa) by controlling the output of the function generator (voltage in mV, shown in parentheses 0 (also known as sham treatment, the amplifier was operated but not connected to the transducer), 0.2 (60), 0.4 (120), 0.6 (180), 0.8 (240), and 1.0 (300). We were not able to deliver stimuli more intense than 1.0 MPa because the transducer was damaged by sustained operation at this level. The protocol for determining the effects of stimulus duration, duty cycle, and pulse repetition frequency (PRF) was analogous, except that we varied the levels of the respective quantities.

Each animal's movement was recorded at 20 frames per second at a resolution of 576 × 592 pixels. We recorded ∼350 frames per video. The resolution and frame rate were chosen to be high enough to provide reliable movement characterization while maintaining acceptable size of the stored videos.

Quantification of response frequency and baseline response frequency.

To detect bona fide ultrasound-evoked reversals, we measured the velocity vector over the interval from 250 ms to 1 s following the ultrasound onset and that during a 1 s period immediately preceding the ultrasound onset. Next, we computed the vector difference and evaluated the magnitude of that difference. We asked whether this metric differed from the null distribution constructed over all baseline recordings (same time windows, just shifted 1 s back in time so that there could be no effect of ultrasound) available for a given animal. If the vector difference was unlikely (p < 0.01) to have been drawn from the null distribution, we classified the response as a reversal. We computed the proportion of significant responses over the 10 stimulus repetitions for each animal and refer to this metric as the response frequency.

The computation of the baseline response frequency (dashed lines in the plots) was analogous to the computation of the response frequency with the exception that the metrics were taken in time windows before the ultrasound could have any impact (i.e., before the ultrasound was turned on). In particular, the velocity difference was computed by comparing a 1 s time window immediately preceding the ultrasound to a 1 s time window preceding the ultrasound onset by 1 s. The baseline distribution used the same time windows, just shifted back in time by 1 s. The baseline response frequency was indistinguishable across the tested animal strains (F(4,95) = 0.28, p = 0.90, one-way ANOVA).

Simulation of the relationship between reversal frequency and duty cycle.

The simulation shown in Figure 5B was derived as follows. First, we assumed that the signal relevant to modulation of reversal behavior is the envelope modulating the carrier frequency of the ultrasound transducer. Next, for each duty cycle value we tested, we converted this signal into the frequency domain using the function fft in MATLAB (The MathWorks). Finally, this signal was convolved with the filter computed by Eastwood et al. (2015), and the effective (root mean squared or rms) value of the resulting signal was taken as the model's output. Thus, this simulation rests on the delivered stimulus waveform and the previously proposed mechanical filter. The only adjustment was a linear scaling factor used to plot the results on common graph. We note that the filter reported by Eastwood et al. (2015) was defined over the range from 1 Hz to 3 kHz, which we extrapolated to 10 kHz.

Temperature measurements.

We used a dual sensing fiberoptic hydrophone (Precision Acoustics) to measure temperature (Morris et al., 2009). This device uses optical interferometry to record acoustically and thermally induced changes in thickness of its sensing membrane by ultrasound. When performing the measurements, we immersed the tip of the optic fiber into the agar at the location of maximal pressure. The device, which recorded the temperature at a 200 Hz sampling rate, rapidly registered the changes in temperature to ultrasound onset. As expected, temperatures were highest at the end of the ultrasound burst; consequently, we report the difference in temperature between the start and end of the ultrasound burst.

Simulation of acoustic radiation force.

We used the k-Wave simulation tool (Treeby and Cox, 2010) to estimate the acoustic radiation force that acted on the animals in our setup. The simulation used the same geometry and media as our setup, including water, agar (4-mm-thick), worm on the agar (0.05 × 0.05 × 1.0 mm), and air. The speed of sound (m s−1), density (kg m−3), and acoustic attenuation coefficient (dB cm−1 MHz−1) for these media were, respectively, set to the following: (1540, 1000, 0.0022), (1548, 1000, 0.40), (1562, 1081, 1.2), and (343, 1.2, 1.6) (Barber et al., 1970; Kremkau et al., 1981). The simulation grid was computed in steps of 1 ns in time and 10 μm in space. We set the stimulus level such that the simulated pressure field had an amplitude of 1 MPa at focus. We obtained the steady-state time-average radiation force density field from the steady-state time-average intensity field provided by the simulation. To obtain the net radiation force, we integrated the force density field over a volume that approximated the animal's head (0.05 × 0.05 × 0.2 mm). This resulted in a net force magnitude of 873 nN. Integrating the force density over a larger volume had only a small influence on the results because the force field was strongest near the animal's head (see Fig. 1A). In Figure 1B, the force density is integrated within cubes of 0.05 × 0.05 × 0.05 mm, using acoustic parameters of the worm.

Experimental design and statistical analysis.

The C. elegans hermaphrodites used in this study were age-synchronized by hypochlorite treatment (Stiernagle, 2006) and cultivated at 15°C or 20°C. There were no detectable effects of cultivation temperature, ambient temperature, or humidity on the ultrasound-evoked responses, which were collected over a period of months by two members of the research team. All behavioral recordings were performed blind to genotype, and we determined whether or not a given trial included a reversal event as described above. The response rate was computed as the fraction of 10 trials that evoked a reversal, and results were pooled from n = 20 C. elegans subjects.

We used ANOVAs to assess the effect of stimulus pressure, duration, PRF, and duty cycle. In this linear model (see Figs. 3, 5), the dependent variable is the response frequency as described above. The independent variable (the one factor) is pressure (see Fig. 3), or PRF or duty cycle (see Fig. 5). We report the p values as well as the F statistic and the degrees of freedom of the tests.

Results

As a first step toward determining how animals and excitable tissues detect and respond to ultrasound stimulation, we placed single adult wild-type (N2) nematodes on sterile agar slabs and tracked their movement using a digital video camera and the Parallel Worm Tracker (Ramot et al., 2008), adapted for tracking single animals (Fig. 1A). We subjected each animal to pulsed ultrasound (10 MHz frequency, 200 ms duration, 1 kHz pulse repetition frequency (PRF) at 50% duty cycle; Fig. 1B) when it approached the ultrasound focus and found that this stimulus elicits similar reversal behaviors over the course of 10 trials (Fig. 2). The effect was due to ultrasound stimulation per se because sham stimuli (0 MPa) did not increase reversals above their basal or unstimulated rate (Fig. 2A,B).

Figure 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 1.

A system for delivering pulsed ultrasound to C. elegans nematodes. A, Schematic side view of the setup, showing a single wild-type adult hermaphrodite crawling on the surface of agar slab, tracked by a digital video camera, and maintained within the field of view by a copper sulfate boundary. A piezoelectric ultrasound transducer (10 MHz carrier frequency, line-focused) is coupled directly to the bottom of the agar slab by a column of degassed water. B, Simulation of the distribution of radiation force expected from the line-focused 10 MHz ultrasound transducer (see Materials and Methods). Animals were stimulated only during forward movement, as they entered the zone corresponding to the highest expected pressure. C, Schematic of a typical stimulus consisting of ultrasound pulses delivered at PRF of 1 kHz for a total duration of 200 ms at 50% duty cycle. In this study, we systematically varied the applied pressure, pulse duration, PRF, and duty cycle.

Figure 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 2.

Ultrasound elicits reversal behavior in a pressure- and stimulus time-dependent manner in wild-type C. elegans. A, B, Raster plots showing the response of 20 animals (10 trials/animals) to a 200 ms sham stimulus (0 MPa pressure, A) and a bona fide stimulus (1.0 MPa, B). Heading angle is encoded in color such that headings similar to the average angle in the 1 s window immediately preceding stimulus onset are blue and reversals are encoded in yellow. Rows correspond to single trials and blocks are 10 trials delivered to each animal; traces were smoothed with a zero-lag rectangular sliding 150 ms window. Top, Silhouettes depict representative responses to sham (A) and 1.0 MPa stimuli (B). C, Reversal frequency increases with applied pressure. Points indicate mean ± SEM (n = 20) for animals stimulated at each of the six pressure values for a total of 10 trials. Solid line indicates a Boltzmann fit to the data with an P1/2 of 0.71 MPa, a slope factor of 0.15, and a maximum probability of 83%. Dotted line indicates the unstimulated reversal rate. Stimulus parameters: 1 kHz, 50% duty cycle, 200 ms pulse duration, variable pressure. D, Reversal probability increases with stimulus duration. Points indicate mean ± SEM (n = 20). Smooth line indicates an exponential fit to the data with a time constant of 90 ms. Stimulus parameters: 1 kHz, 50% duty cycle, variable pulse duration, 1.0 MPa pressure. C, D, Dotted line indicates baseline rate of responding (see Materials and Methods). Smooth line indicates an exponential fit to the data with a time constant of 90 ms.

Behavioral responses were robust across trials for a given individual and among all animals tested (Fig. 2B). We determined whether each animal's response to ultrasound stimulation was statistically different from spontaneous changes in direction (see Materials and Methods). For each animal, we quantified the proportion of significant responses over the 10 stimulus trials, and refer to this metric as response frequency.

The response frequency increased with ultrasound pressure (Fig. 2C) and was indistinguishable from the spontaneous reversal rate for the sham stimulus (Fig. 2A, dotted line; p = 0.52, t test, n = 20). The reversal rate increased above the baseline at a pressure of 0.6 MPa (p < 10−6). At 1 MPa (Fig. 2B), there was a significant response in 77.5% of trials, on average. The relationship between response frequency and applied pressure was well-described by a sigmoid function as follows: Embedded Image where F is the response frequency, P is the pressure, and base is the spontaneous reversal frequency.

For wild-type animals stimulated by a 200 ms pulse (1 kHz PRF, 50% duty cycle), the best fit parameters were Fmax = 83%; P1/2 = 0.71 MPa, slope = 0.15, base = 5%. Thus, the half-activation pressure equals 0.71 MPa. A one-way ANOVA also detected a significant modulation of the response frequency by pressure (F(5,114) = 103.4, p < 10−39), reinforcing the idea that the probability of ultrasound-induced reversal depends on stimulus pressure.

We also tested the effect of varying the total duration of the ultrasound stimulus (Fig. 2D), holding all other parameters (i.e., pressure, duty cycle, PRF) constant. In agreement with a previous study (Ibsen et al., 2015), responses were weak or absent when the stimulus was brief. There was a significant modulation of the response frequency by stimulus duration (one-way ANOVA, F(3,76) = 30.8, p < 10−12). Stimuli of 100 ms in duration or longer produced substantial effects. The response frequency did not increase substantially beyond stimulus duration of 200 ms (response frequency at 200 ms vs 400 ms: p = 0.24, paired t test, n = 20). Therefore, we used a stimulus duration of 200 ms for subsequent experiments.

In principle, ultrasound-evoked behaviors could depend on thermosensation, mechanosensation, or both. We used a genetic approach to distinguish among these possibilities, leveraging mutants deficient in thermosensation or mechanosensation. To test for thermal effects, we compared ultrasound-evoked behaviors in wild-type and gcy-23(nj37)gcy-8(oy44)gcy-18(nj38) that lack a trio of receptor guanylate cyclases expressed exclusively in the AFD thermoreceptor neurons and are defective in thermotaxis (Garrity et al., 2010; Glauser and Goodman, 2016). Although these mutants have an intact AFD thermoreceptor neuron, they lack the ability to sense tiny (<0.05°C) thermal fluctuations in temperature (Ramot et al., 2008; Wasserman et al., 2011). As shown in Figure 3A, the response of these thermosensory-defective mutants was indistinguishable from that of wild-type animals. The mutants retained modulation by stimulus pressure, as assessed by one-way ANOVA (F(5,114) = 80.7, p < 10−35). Furthermore, as expected from the plot, a two-way ANOVA with factors animal strain and pressure failed to detect a significant difference between the strains (F(1,228) = 0.02, p = 0.89) as well as the strain × pressure interaction (F(5,228) = 1.40, p = 0.23). Thus, the ability to sense tiny (<0.05°C) thermal fluctuations is not required for ultrasound-induced reversal behaviors.

Figure 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 3.

Loss of mechanosensation, but not thermosensation, disrupts ultrasound-evoked reversals. A, B, Pressure-response curves of wild-type N2 animals (blue) compared with a thermosensation-defective mutant (orange) and three mechanosensation-defective mutants (black): mec-3(e1338), mec-4(e1611), and mec-4(u253). Points indicate mean ± SEM (n = 20 animals tested in 10 trials/animals). Smooth lines indicate fit to the data according to a sigmoidal function. The data and fit for wild-type are the same as in Figure 2C. Fitting parameters for gcy-23(nj37)gcy-8(oy44)gcy-18(nj38) are (Fmax, P1/2, slope, base): 80%, 0.76 MPa, 0.10, 9%. Dotted lines indicate the average baseline response rate (see Materials and Methods) for each case. There was no significant effect of genotype on baseline reversal rates. Stimulus parameters: 1 kHz pulse frequency, 50% duty cycle, 200 ms duration, variable pressure. C, Response rate as a function of genotype. Error bars indicate mean (SEM; n = 20) reversal rate. Annotations below the graph indicate the nature of the known sensory deficit associated with each genotype. Animals were tested as young adult hermaphrodites and blind to genotype.

Having established that thermosensation is dispensable for ultrasound-evoked reversals, we compared responses in wild-type animals and mutants defective in mechanosensation. Specifically, we sought to quantify ultrasound-evoked responses in mutants in which selected mechanoreceptor neurons fail to properly differentiate during development, degenerate, or lack essential, pore-forming subunits of known sensory mechano-electrical transduction channels: MEC-4 (O'Hagan et al., 2005), TRP-4 (Kang et al., 2010). The goal of these experiments was to identify the neurons most likely to serve as the first responders to ultrasound stimulation and to determine whether or not such sensitivity relied upon known mechano-electrical transduction channels.

The mec-3(e1338) mutants fail to generate three sets of neurons known to participate in gentle and harsh touch sensation and are insensitive to both gentle and harsh touch (Way and Chalfie, 1989). The six touch receptor neurons (TRNs) are required for sensing gentle touch, and the two pairs of multidendritic PVD and FLP neurons act as polymodal sensors of mechanical and nociceptive stimuli (Schafer, 2015). We found that mec-3 mutants are insensitive to ultrasound stimulation (Fig. 3B): the mec-3 mutants showed no significant modulation of the response frequency by pressure (F(5,114) = 1.18, p = 0.32, one-way ANOVA). A two-way ANOVA detected both a highly significant difference between the strains (F(1,228) = 246.1, p < 10−37) and a highly significant strain × pressure interaction (F(5,228) = 56.8, p < 10−37). This effect was specific to ultrasound-evoked reversals; mutants moved at an average speed that was indistinguishable from wild-type animals (speed measured during 1 s period preceding ultrasound onset; wild-type: 0.21 mm/s; mec-3: 0.17 mm/s; p = 0.11, n = 20, t test). These average speeds are within the range of values reported previously for wild-type animals (Ramot et al., 2008). This result shows that the mec-3-dependent mechanoreceptor neurons are required for ultrasound-evoked reversals and suggests that ultrasound can exert forces on neural tissue sufficient to activate these neurons.

We narrowed the search to a subset of the mec-3-dependent mechanoreceptor neurons by testing ultrasound-evoked behavior in mec-4(e1611) mutants in which the six TRNs degenerate and the PVD and FLP mechanoreceptor neurons are intact (Driscoll and Chalfie, 1991). As found in mec-3 mutants, ultrasound failed to evoked reversals in mec-4(e1611) (Fig. 3C), and there was no significant modulation of the response frequency by the ultrasound pressure amplitude in these animals (Fig. 3C; F(5,114) = 1.47, p = 0.20). Moreover, a two-way ANOVA detected a significant difference between the mutant and wild-type strains and a highly significant strain × pressure interaction (both p < 10−36). Thus, the TRN neurons, which can detect forces as small as 50 nN (O'Hagan et al., 2005), are required for behavioral responses to ultrasound stimulation in C. elegans.

Next, we investigated proteins expressed in the TRN neurons that might mediate the effect. Of particular interest, the TRNs express mec-4, which encodes a non–voltage-gated sodium channel of the DEG/ENaC/ASIC family required for touch-evoked reversals. MEC-4 is expressed exclusively in the TRNs and is an essential pore-forming subunit of the mechanosensitive ion channel activated by mechanical loads applied directly to the animal's skin (O'Hagan et al., 2005). Like mec-3 and mec-4(e1611) mutants, mec-4(u253) null mutants are insensitive to ultrasound stimulation (Fig. 3D). These animals have intact TRNs but lack the MEC-4 protein required for mechanotransduction and showed no significant modulation of the response frequency by the ultrasound pressure (F(5,114) = 0.37, p = 0.87). A two-way ANOVA revealed a significant effect of genotype× pressure interaction (both p < 10−35, two-way ANOVA). Although the response of Mec mutants appeared as if it might be modulated by pressure (Fig. 3B–D), this apparent modulation was not significant (p > 0.09, one-way ANOVA). Collectively, these results establish that behavioral responses to focused ultrasound depend on the TRN neurons and the MEC-4 protein.

Thus far, we have shown that focused ultrasound evokes reversal behaviors in freely moving C. elegans nematodes in a pressure- and stimulus duration-dependent manner (Fig. 2) and that such responses depend on the animal's ability to detect mechanical, but not thermal, stimuli (Fig. 3). Together, these findings imply that ultrasound exerts its effects on excitable tissues via a nonthermal, mechanical mechanism. Although additional work is needed to determine how ultrasound produces these effects, a leading possibility is the generation of mechanical strain in neurons expressing mechanosensitive ion channels, such as MEC-4.

MEC-4 is not the only protein thought to form a mechanosensitive ion channel in C. elegans nematodes. The TRP-4 protein is expressed in the CEP mechanoreceptor neurons and is an ortholog of the Drosophila NOMPC channel (Li et al., 2006) known to form mechanosensitive ion channels (Yan et al., 2013). A previous study showed that C. elegans responds to ultrasound-induced cavitation of microbubbles and proposed that these responses were TRP-4-dependent (Ibsen et al., 2015). To determine whether TRP-4 also contributed to ultrasound-evoked behaviors elicited in the absence of microbubbles, we analyzed the same trp-4 strain used by Ibsen et al. (2015): VC1141 trp-4(ok1605). In agreement with the prior report (Ibsen et al., 2015), we observed a modest deficit in ultrasound-evoked behavior (Fig. 4A). A two-way ANOVA detected both a main effect of strain (F(1,228) = 17.8, p < 0.0001) and a significant strain × pressure interaction (F(5,228) = 4.8, p = 0.0003). The defect in these mutants was not specific for ultrasound-evoked behaviors, however: trp-4 mutants had a lower average speed than wild-type mutants under baseline conditions (0.17 vs 0.21 mm/s, p = 0.0086, t test, n = 20).

Figure 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 4.

Strains carrying deletions in the trp-4 NOMPC channel gene differ in their response to ultrasound stimulation. A, Pressure-response curves of wild-type (blue), VC1141 trp-4(ok1605) (magenta) mutants used in a previous study (Ibsen et al., 2015). Smooth curve fit to the trp-4(ok1605) data yielded Fmax = 65%; P1/2 = 0.83 MPa, slope = 0.13, base = 8%. B, Pressure response curves of three other trp-4 deletion mutants were indistinguishable from wild-type. VC818 trp-4(gk341), TQ296 trp-4(sy695), and GN716 trp-4(ok1605) mutants, which was derived from VC1141 by outcrossing four times with wild-type (N2) animals. C, Response rate in four trp-4 mutant strains. Error bars indicate mean ± SEM reversal rate evoked by ultrasound stimulation with the following parameters: 1 kHz, 50% duty cycle, 200 ms pulse duration, 1.0 MPa. Dotted line indicates the average baseline response rate (see Materials and Methods). The number of animals analyzed across 10 trials is indicated in parentheses. We used PCR to verify that all strains harbored the expected deletions in the trp-4 locus (Materials and Methods). Wild-type data are from Figure 2C.

Because the ok1605 allele encodes a partial in-frame deletion and because we also observed that these mutants grew slowly compared with wild-type animals, we tested two additional deletions in the trp-4 gene: gk341 and sy695. All three alleles (ok1605, gk341, and sy695) encode deletions in the trp-4 gene, which we verified by PCR analysis of genomic DNA (see Materials and Methods). Despite the expectation that that the three trp-4 alleles would have the same ultrasound phenotype, we found that gk341 and sy695 mutants responded to ultrasound just like wild-type animals (Fig. 4B; two-way ANOVAs, main effects and interactions, p > 0.29).

These findings suggested that the deficit in the VC1141 trp-4(ok1605) animals might be due to a mutation present in the genetic background. We tested this idea by outcrossing the trp-4(ok1605) animals against wild-type (N2) animals four times while tracking the trp-4 mutation via PCR. Animals from this new strain, GN716 trp-4(ok1605), had ultrasound-evoked behaviors that were indistinguishable from wild-type (Fig. 4B; two-way ANOVA, main effect and interaction, p > 0.23). These results are summarized for the pressure of 1 MPa in Figure 4C and suggest that the defect we and others (Ibsen et al., 2015) have observed in VC1141 trp-4(ok1605) animals is due to mutation/s in the genetic background of this strain.

The finding that mechanosensation is an essential component of the behavioral effects of ultrasound suggests that there might be an optimal frequency of the ultrasound delivery that matches the mechanical properties of the tissue. We investigated this possibility by varying PRFs in the range from 30 Hz to 10 kHz while keeping pulse duration, pressure, and duty cycle constant. We found that ultrasound indeed evoked reversals in a PRF-dependent manner (Fig. 5A). Response increased with PRF, reached a maximal value in the range of 300–1000 Hz, and decreased at higher frequencies. The shape of the curve is reminiscent of the prediction (Fig. 5A, green) of a model linking indentation to mechanical strain and MEC-4-dependent channel activation (Eastwood et al., 2015). We note that, because stimuli were delivered at 50% duty cycle at all the tested frequencies, the same amount of energy was delivered at all PRFs. If the behavioral responses were the result of tissue heating, little or no modulation by the PRF would be expected. Yet, the plot shows, and an ANOVA confirms, a strong modulation of the response by the PRF (F(5,114) = 10.8, p < 10−7).

Figure 5.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 5.

Ultrasound efficacy depends on pulse repetition frequency and duty cycle. A, The response (black circles, mean ± SEM) of wild-type animals to a train of ultrasound pulses plotted as a function of pulse repetition frequency. The duty cycle was held constant at 50% in all cases, ensuring that all stimuli deliver the same amount of energy. The smooth curve (green) shows a simulation of the sensitivity of TRN currents to sinusoidal mechanical indentations (Eastwood et al., 2015). B, The response (black circles, mean ± SEM) of wild-type animals as a function of duty cycle. Because the pulse repetition rate was 1 kHz in all cases, a duty cycle of 5, 10, 25, 50, 75, 100% corresponds to a pulse width of 50 μs, 100 μs, 250 μs, 500 μs, 750 μs, and 1 ms (continuous wave, no off epochs), respectively. The green curve shows the level of TRN activation expected from the frequency filtering shown in panel A (Materials and Methods for modeling details). Temperature increase (orange circles, mean ± SEM, n = 5) as a function of duty cycle. The smooth curve (orange) is a quadratic fit to the data and is included for visual clarity. A, B, The carrier frequency was 10 MHz, the stimulus amplitude was 1 MPa, and the stimulus duration 200 ms. The dotted line shows the unstimulated reversal frequency.

We further hypothesized that discrete pulses may be more potent in eliciting mechanical effects because discrete pulses deliver multiple discrete mechanical events into the tissue. To test this idea, we varied the duty cycle while holding stimulus duration, PRF, and pressure values constant at 200 ms, 1 kHz, and 1 MPa, respectively. Figure 5B shows the relationship between response frequency and duty cycle. It reveals that a duty cycle of 50% was more than threefold more potent than a 100% duty cycle (77.5% compared with 24.0%, p < 10−12, t test), even though a continuous stimulation delivers twice as much energy into the tissue as the pulsed protocol of 50% duty. In line with this finding, pulsed ultrasound stimulation has been found more effective than continuous stimulation in eliciting motor responses in rats (Kim et al., 2014). That study also found the value of 50% duty to be optimal. It is important to note that the 24.0% response rate for the continuous stimulus (100% duty) is above baseline (p < 0.001, t test, n = 20). Thus, although a pulsed stimulus is more effective than a continuous stimulus, pulsing the ultrasound is not necessary to elicit a detectable response. The figure further shows that the width of the individual mechanical events associated with the ultrasound can be quite brief, just 50 μs (5% of duty), and still trigger appreciable behavioral responses (response rate of 34.0%, significantly different from baseline at p < 0.0001, t test, n = 20). This is even though the energy delivered into the tissue is only 1/10th of that delivered at 50% duty.

Whereas heating increases linearly with duty cycle (Fig. 5B, orange plot), behavioral response frequency had a nonlinear dependence on duty cycle. The heating effect is expected from the fact that the energy delivered in the tissue increases with duty cycle. To ask whether or not the dependence on duty cycle could be explained by the frequency dependence of TRN activation, we simulated the frequency distribution expected as a function of duty cycle and combined with this the model from Eastwood et al. (2015). This model matched the experimental results (Fig. 5B, green line). Collectively, the effect of duty cycle reinforces the idea that behavioral responses to ultrasound are mediated by mechanical effects and not by heating.

Discussion

We sought to illuminate the biophysical mechanisms that underlie ultrasound stimulation of excitable cells. To do so, we used C. elegans as a model, harnessing its well-characterized and compact nervous system and comprehensive library of animals with specific genetic interventions. This animal has an extraordinary ability to detect tiny thermal fluctuations and mechanical stimuli (O'Hagan et al., 2005; Ramot et al., 2008). We found that ultrasound elicits robust reversal behaviors and that the response probability depends on stimulus intensity, duration, and specific pulsing protocols. Sensitivity to ultrasound and its modulation by pressure are preserved in mutants deficient in thermosensation and eliminated in mutants defective in mechanosensation. These findings are in agreement with a report (Zhou et al., 2017), showing that tax-4 mutants, which are defective in thermosensing, retain the ability to respond to brief pulses of high-frequency ultrasound. Consistent with nonthermal, mechanical activation of sensory neurons linked to reversal behaviors, the response probability exhibited optima in both PRF and duty cycle.

Ultrasound-evoked reversal responses required expression of MEC-4, a key subunit of a touch-activated mechanosensitive ion channel. This finding implies that ultrasound can activate neurons by acting on mechanosensitive ion channels. Notably, because the MEC-4 ion channel complex is insensitive to changes in membrane voltage or capacitance, our findings are inconsistent with the hypothesis that ultrasound acts by inducing changes in membrane capacitance (Krasovitski et al., 2011; Plaksin et al., 2014).

We tested whether another mechanosensitive channel might contribute to ultrasound-evoked behaviors by analyzing strains carrying deletions in the trp-4 gene. Prior work showed that TRP-4 is a pore-forming subunit of mechanosensitive channel expressed by texture-sensing neurons in the worm's head (Kang et al., 2010) and suggested that this protein could sensitize neurons to ultrasound stimulation. Yet, we did not detect any effect of the loss of TRP-4 channels on ultrasound-evoked responses in three independent strains carrying validated deletions in the trp-4 gene. We did detect a decrease in ultrasound sensitivity in a fourth strain (VC1141) that was used in a previous study (Ibsen et al., 2015). However, this phenotype is not due to loss of trp-4 function because four rounds of outcrossing eliminated it. Rather, the partial loss of ultrasound sensitivity is likely to be due to an unidentified mutation(s) in the VC1141 strain. Additional investigations will be needed to identify the affected gene(s), an effort that could reveal additional genetic factors regulating sensitivity to ultrasound stimulation.

This study provides evidence that ultrasound can stimulate neurons through its mechanical mode of action. Within the mechanical domain, there can be several specific candidate mechanisms at play. First, ultrasound may elicit cavitation, a phenomenon characterized by formation and collapse of gaseous bodies in liquid media or soft tissues. However, for frequencies >1 MHz, cavitation requires pressures >5 MPa and the cavitation threshold for 10 MHz is even higher (Nightingale et al., 2015). Thus, both the 10 MHz transducer and low pressures we used make cavitation unlikely. Second, the incident tissue, such as a cell membrane, experiences oscillations with period equal to the ultrasound carrier frequency. The pressures used for neuromodulation can cause appreciable particle displacement (on the order of 0.01–0.1 μm) (Gavrilov et al., 1976). Nonetheless, the displacement is distributed in sinusoidal fashion along the wavelength (∼100 μm at 10 MHz) of the propagating wave. This creates a very small displacement gradient (e.g., 0.1 μm per 100 μm). It is questionable whether such a small gradient can cause significant enough deformation of a pore segment of an ion channel with regard to the channel dimensions. Moreover, the primary pressure oscillations, which occur at a specific carrier frequency, cannot explain the frequency dependence of the responses (Fig. 5A). The third and most probable form of mechanical energy underlying the effects in this study is the acoustic radiation force (Trahey et al., 2004; Sarvazyan et al., 2010; Iversen et al., 2017). Acoustic radiation forces result from differences in acoustic intensities at individual points in space. The differences can be caused by ultrasound absorption, scattering, reflection, or other phenomena, and lead to net forces on the tissue (Duck et al., 1998). Acoustic radiation force exerts a steady pressure on a target throughout the time of ultrasound application. This steady pressure may stretch a cell membrane to an extent that affects conformation states of ion channels or other active molecules tied to the membrane. A simulation of the propagating ultrasound field, for the pressure of 1 MPa, revealed that a net acoustic radiation force of 873 nN can act on the animal's head (see Materials and Methods). This exceeds the animal's sensitivity threshold to mechanical forces, which is 50–100 nN (Petzold et al., 2013). Thus, the acoustic radiation force expected for a 1 MPa stimulus is sufficient to engage the animal's mechanosensation. It is worth noting that the radiation force acts during the On epochs of the ultrasound (Fig. 1C, black rectangles), and not during the Off epochs when the ultrasound amplitude is zero. This way, pulsed ultrasound delivers force pulses at a specific PRF, and there can therefore be a modulation by the PRF (Fig. 5A).

Our findings implicate mechanical force as a major physical effect of ultrasound on neurons and their ion channels, delineating a complete pathway from mechanical force to activation of excitable cells. In this light, it is tempting to reiterate a potential unifying mechanism linking ultrasound to activation of excitable cells (Tyler, 2011). Suppose that ultrasound exerts similar mechanical effects in complex nervous tissues as it does in C. elegans. In addition, suppose that ultrasound deforms tissue and generates mechanical strain in neurons sufficient to activate mechanosensitive ion channels, as in C. elegans. Ion channels likely to subserve this function in mammals include the intrinsically mechanosensitive K2P family of potassium channels (Brohawn et al., 2014) and piezo channels (Syeda et al., 2016) known to be expressed in the brain. In support of this idea, the activity of K2P channels, including TREK-1, TREK- 2, and TRAAK, is potentiated by ultrasound stimulation in heterologous cells (Kubanek et al., 2016). Sensitivity to the mechanical effects of ultrasound might not be limited to primarily mechanosensitive channels. For instance, voltage-gated sodium channels have been implicated in activation of neurons by ultrasound (Tyler, 2011, 2012; Kubanek et al., 2016) and are known to be sensitive to membrane tension (Beyder et al., 2010).

Exactly how ultrasound stimulation is translated into local mechanical strain will depend on the material properties of the excitable tissues under study. In the case of C. elegans, we found that the PRF and duty cycle dependence of ultrasound-evoked behaviors agreed with a model that links tissue indentation associated with a mechanical stimulus to the activation of MEC-4-dependent channels (Eastwood et al., 2015). The correspondence (Fig. 5A) suggests that ultrasound-induced radiation force elicits a profile of mechanical strain similar to that produced by indentation with a physical probe. The model (Eastwood et al., 2015) also provides insight into why stimuli of 50% duty are the most potent (Fig. 5B): stimuli delivered at lower or higher duty cycle values contain a majority of their energy in high-frequency harmonics (relative to the 1 kHz PRF), and these high-frequency harmonics are filtered out by the tissue (Fig. 5A). However, the relationship between bursts of ultrasound pulses and the mechanical strain we infer that it generates is not currently known, and future measurements of ultrasound-induced strain will be needed to fill this knowledge gap for worms and other tissues. Consistent with the proposal that tissues differ in their mechanical filtering properties, ultrasound-evoked behaviors have an optimal PRF near 500 Hz in C. elegans and show strong effects at ∼3000 Hz in mice (King et al., 2013). These results indicate that an improved understanding of mechanical filtering by soft tissues will be needed to further the long-term goal of applying ultrasound as a noninvasive modality to stimulate excitable cells.

This study reveals that pulsatile forces associated with ultrasound are potent enough to activate mechanically sensitive ion channels in living animals. Given that many ion channels expressed in neurons and glia are mechanically sensitive (Ostrow and Sachs, 2005; Tyler, 2012), this study illuminates one way that ultrasound could influence activity in the brain. Another way would be to sensitize specific neurons to ultrasound through ectopic expression of a channel known to be mechanosensitive, a strategy referred to as sonogenetics (Ibsen et al., 2015). In both scenarios, our work underscores the importance of tuning stimulus parameters to maximize acoustic radiation force and may help to further the development of ultrasound as a noninvasive and spatially precise tool to study the nervous system and potentially to ameliorate neurological disorders.

Footnotes

  • This work was supported by National Institutes of Health Grant K99NS100986 to J.K., Grants R01NS047715 and R01NS092099 to M.B.G., Stanford Neuroscience Institute and Mathers Foundation to M.B.G. and S.A.B., and a Stanford Medicine Dean's Fellowship to J.K. We thank our colleagues at Stanford, Zhiwen Liao, Michael Menz, Patrick Ye, and Martin Prieto for laboratory and technical assistance and input on the experimental design; and Merritt Maduke, Pierre Khuri-Yakub, Kimberly Butts-Pauly, and Julian Brown for helpful comments; Zachary Pincus (Washington University, St. Louis) for providing animals for initial pilot experiments; and Alessandro Sanzeni and Massimo Vergassola for additional model calculations.

  • The authors declare no competing financial interests.

  • Correspondence should be addressed to Dr. Miriam B. Goodman, B-111 Beckman Center, 279 Campus Drive, Stanford, CA 94305. mbgoodman{at}stanford.edu

References

  1. ↵
    1. Altman PL,
    2. Dittmer DS,
    3. Zwemer RL
    (1974) Biology data book. Bethesda, MD: Federation of American Societies for Experimental Biology.
  2. ↵
    1. Barber TW,
    2. Brockway JA,
    3. Higgins LS
    (1970) The density of tissues in and about the head. Acta Neurol Scand 46:85–92. doi:10.1111/j.1600-0404.1970.tb05606.x pmid:4983875
    OpenUrlCrossRefPubMed
  3. ↵
    1. Beyder A,
    2. Rae JL,
    3. Bernard C,
    4. Strege PR,
    5. Sachs F,
    6. Farrugia G
    (2010) Mechanosensitivity of nav1. 5, a voltage-sensitive sodium channel. J Physiol 588:4969–4985. doi:10.1113/jphysiol.2010.199034 pmid:21041530
    OpenUrlCrossRefPubMed
  4. ↵
    1. Brohawn SG,
    2. Su Z,
    3. MacKinnon R
    (2014) Mechanosensitivity is mediated directly by the lipid membrane in TRAAK and TREK1 K+ channels. Proc Natl Acad Sci U S A 111:3614–3619. doi:10.1073/pnas.1320768111 pmid:24550493
    OpenUrlAbstract/FREE Full Text
  5. ↵
    1. Buiochi EB,
    2. Miller RJ,
    3. Hartman E,
    4. Buiochi F,
    5. Bassani RA,
    6. Costa ET,
    7. O'Brien WD
    (2012) Transthoracic cardiac ultrasonic stimulation induces a negative chronotropic effect. IEEE Trans Ultrason Ferroelectr Freq Control 59:2655–2661. doi:10.1109/TUFFC.2012.2506 pmid:23221214
    OpenUrlCrossRefPubMed
  6. ↵
    1. Colucci V,
    2. Strichartz G,
    3. Jolesz F,
    4. Vykhodtseva N,
    5. Hynynen K
    (2009) Focused ultrasound effects on nerve action potential in vitro. Ultrasound Med Biol 35:1737–1747. doi:10.1016/j.ultrasmedbio.2009.05.002 pmid:19647923
    OpenUrlCrossRefPubMed
  7. ↵
    1. Dalecki D
    (2004) Mechanical bioeffects of ultrasound. Annu Rev Biomed Eng 6:229–248. doi:10.1146/annurev.bioeng.6.040803.140126 pmid:15255769
    OpenUrlCrossRefPubMed
  8. ↵
    1. Deffieux T,
    2. Younan Y,
    3. Wattiez N,
    4. Tanter M,
    5. Pouget P,
    6. Aubry JF
    (2013) Low-intensity focused ultrasound modulates monkey visuomotor behavior. Curr Biol 23:2430–2433. doi:10.1016/j.cub.2013.10.029 pmid:24239121
    OpenUrlCrossRefPubMed
  9. ↵
    1. Diaz-Franulic I,
    2. Poblete H,
    3. Min-o-Galaz G,
    4. González C,
    5. Latorre R
    (2016) Allosterism and structure in thermally activated transient receptor potential channels. Annu Rev Biophysics 45:371–398. doi:10.1146/annurev-biophys-062215-011034 pmid:27297398
    OpenUrlCrossRefPubMed
  10. ↵
    1. Driscoll M,
    2. Chalfie M
    (1991) The mec-4 gene is a member of a family of Caenorhabditis elegans genes that can mutate to induce neuronal degeneration. Nature 349:588–593. doi:10.1038/349588a0 pmid:1672038
    OpenUrlCrossRefPubMed
  11. ↵
    1. Duck FA,
    2. Baker AC,
    3. Starritt HC
    (1998) Ultrasound in medicine. Boca Raton, FL: CRC.
  12. ↵
    1. Eastwood AL,
    2. Sanzeni A,
    3. Petzold BC,
    4. Park SJ,
    5. Vergassola M,
    6. Pruitt BL,
    7. Goodman MB
    (2015) Tissue mechanics govern the rapidly adapting and symmetrical response to touch. Proc Natl Acad Sci U S A 112:E6955–E6963. doi:10.1073/pnas.1514138112 pmid:26627717
    OpenUrlAbstract/FREE Full Text
  13. ↵
    1. Foster KR,
    2. Wiederhold ML
    (1978) Auditory responses in cats produced by pulsed ultrasound. J Acoust Soc Am 63:1199–1205. doi:10.1121/1.381829 pmid:649878
    OpenUrlCrossRefPubMed
  14. ↵
    1. Fry F,
    2. Ades HW,
    3. Fry WJ
    (1958) Production of reversible changes in the central nervous system by ultrasound. Science 127:83–84. doi:10.1126/science.127.3289.83 pmid:13495483
    OpenUrlFREE Full Text
  15. ↵
    1. Garrity PA,
    2. Goodman MB,
    3. Samuel AD,
    4. Sengupta P
    (2010) Running hot and cold: behavioral strategies, neural circuits, and the molecular machinery for thermotaxis in C. elegans and Drosophila. Genes Dev 24:2365–2382. doi:10.1101/gad.1953710 pmid:21041406
    OpenUrlAbstract/FREE Full Text
  16. ↵
    1. Gavrilov LR,
    2. Gersuni GV,
    3. Ilyinsky OB,
    4. Sirotyuk MG,
    5. Tsirulnikov EM,
    6. Shchekanov EE
    (1976) The effect of focused ultrasound on the skin and deep nerve structures of man and animal. Prog Brain Res 43:279–292. doi:10.1016/S0079-6123(08)64360-5 pmid:1257484
    OpenUrlCrossRefPubMed
  17. ↵
    1. Gavrilov LR,
    2. Tsirulnikov EM,
    3. Davies IA
    (1996) Application of focused ultrasound for the stimulation of neural structures. Ultrasound Med Biol 22:179–192. doi:10.1016/0301-5629(96)83782-3 pmid:8735528
    OpenUrlCrossRefPubMed
  18. ↵
    1. Glauser DA,
    2. Goodman MB
    (2016) Molecules empowering animals to sense and respond to temperature in changing environments. Curr Opin Neurobiol 41:92–98. doi:10.1016/j.conb.2016.09.006 pmid:27657982
    OpenUrlCrossRefPubMed
  19. ↵
    1. Harvey EN
    (1929) The effect of high frequency sound waves on heart muscle and other irritable tissues. Am J Physiol 91:284–290. doi:10.1152/ajplegacy.1929.91.1.284
    OpenUrlCrossRef
  20. ↵
    1. Husson SJ,
    2. Costa WS,
    3. Schmitt C,
    4. Gottschalk A
    (2005) Keeping track of worm trackers. WormBook 22:1–17. doi:10.1895/wormbook.1.156.1 pmid:23436808
    OpenUrlCrossRefPubMed
  21. ↵
    1. Ibsen S,
    2. Tong A,
    3. Schutt C,
    4. Esener S,
    5. Chalasani SH
    (2015) Sonogenetics is a non-invasive approach to activating neurons in Caenorhabditis elegans. Nat Commun 6:8264. doi:10.1038/ncomms9264 pmid:26372413
    OpenUrlCrossRefPubMed
  22. ↵
    1. Iversen MM,
    2. Christensen DA,
    3. Parker DL,
    4. Holman HA,
    5. Chen J,
    6. Frerck MJ,
    7. Rabbitt RD
    (2017) Low-intensity ultrasound activates vestibular otolith organs through acoustic radiation force. J Acoust Soc Am 141:4209–4219. doi:10.1121/1.4984287 pmid:28618821
    OpenUrlCrossRefPubMed
  23. ↵
    1. Kang L,
    2. Gao J,
    3. Schafer WR,
    4. Xie Z,
    5. Xu XZ
    (2010) C. elegans TRP family protein TRP-4 is a pore-forming subunit of a native mechanotransduction channel. Neuron 67:381–391. doi:10.1016/j.neuron.2010.06.032 pmid:20696377
    OpenUrlCrossRefPubMed
  24. ↵
    1. Kim H,
    2. Chiu A,
    3. Lee SD,
    4. Fischer K,
    5. Yoo SS
    (2014) Focused ultrasound-mediated non-invasive brain stimulation: examination of sonication parameters. Brain Stimul 7:748–756. doi:10.1016/j.brs.2014.06.011 pmid:25088462
    OpenUrlCrossRefPubMed
  25. ↵
    1. King RL,
    2. Brown JR,
    3. Newsome WT,
    4. Pauly KB
    (2013) Effective parameters for ultrasound-induced in vivo neurostimulation. Ultrasound Med Biol 39:312–331. doi:10.1016/j.ultrasmedbio.2012.09.009 pmid:23219040
    OpenUrlCrossRefPubMed
  26. ↵
    1. Krasovitski B,
    2. Frenkel V,
    3. Shoham S,
    4. 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 pmid:21300891
    OpenUrlAbstract/FREE Full Text
  27. ↵
    1. Kremkau FW,
    2. Barnes RW,
    3. McGraw CP
    (1981) Ultrasonic attenuation and propagation speed in normal human brain. J Acoust Soc Am 70:29–38. doi:10.1121/1.386578
    OpenUrlCrossRef
  28. ↵
    1. Kubanek J,
    2. Shi J,
    3. Marsh J,
    4. Chen D,
    5. Deng C,
    6. Cui J
    (2016) Ultrasound modulates ion channel currents. Sci Rep 6:24170. doi:10.1038/srep24170 pmid:27112990
    OpenUrlCrossRefPubMed
  29. ↵
    1. Lee W,
    2. Kim HC,
    3. Jung Y,
    4. Chung YA,
    5. Song IU,
    6. Lee JH,
    7. Yoo SS
    (2016) Transcranial focused ultrasound stimulation of human primary visual cortex. Sci Rep 6:34026. doi:10.1038/srep34026 pmid:27658372
    OpenUrlCrossRefPubMed
  30. ↵
    1. Li W,
    2. Feng Z,
    3. Sternberg PW,
    4. Xu XZ
    (2006) A C. elegans stretch receptor neuron revealed by a mechanosensitive TRP channel homologue. Nature 440:684–687. doi:10.1038/nature04538 pmid:16572173
    OpenUrlCrossRefPubMed
  31. ↵
    1. Menz MD,
    2. Oralkan O,
    3. Khuri-Yakub PT,
    4. Baccus SA
    (2013) Precise neural stimulation in the retina using focused ultrasound. J Neurosci 33:4550–4560. doi:10.1523/JNEUROSCI.3521-12.2013 pmid:23467371
    OpenUrlAbstract/FREE Full Text
  32. ↵
    1. Meyers R,
    2. Fry WJ,
    3. Fry FJ,
    4. Dreyer LL,
    5. Schultz DF,
    6. Noyes RF
    (1959) Early experiences with ultrasonic irradiation of the pallidofugal and nigral complexes in hyperkinetic and hypertonic disorders. J Neurosurg 16:32–54. doi:10.3171/jns.1959.16.1.0032 pmid:13621264
    OpenUrlCrossRefPubMed
  33. ↵
    1. Mihran RT,
    2. Barnes FS,
    3. Wachtel H
    (1990) Temporally-specific modification of myelinated axon excitability in vitro following a single ultrasound pulse. Ultrasound Med Biol 16:297–309. doi:10.1016/0301-5629(90)90008-Z pmid:2363236
    OpenUrlCrossRefPubMed
  34. ↵
    1. Morris P,
    2. Hurrell A,
    3. Shaw A,
    4. Zhang E,
    5. Beard P
    (2009) A Fabry-Perot fiber-optic ultrasonic hydrophone for the simultaneous measurement of temperature and acoustic pressure. J Acoust Soc Am 125:3611–3622. doi:10.1121/1.3117437 pmid:19507943
    OpenUrlCrossRefPubMed
  35. ↵
    1. Naor O,
    2. Krupa S,
    3. Shoham S
    (2016) Ultrasonic neuromodulation. J Neural Eng 13:031003. doi:10.1088/1741-2560/13/3/031003 pmid:27153566
    OpenUrlCrossRefPubMed
  36. ↵
    1. Nightingale KR,
    2. Church CC,
    3. Harris G,
    4. Wear KA,
    5. Bailey MR,
    6. Carson PL,
    7. Jiang H,
    8. Sandstrom KL,
    9. Szabo TL,
    10. Ziskin MC
    (2015) Conditionally increased acoustic pressures in nonfetal diagnostic ultrasound examinations without contrast agents: a preliminary assessment. J Ultrasound Med 34:1–41. doi:10.7863/ultra.34.7.15.13.0001 pmid:26112617
    OpenUrlAbstract/FREE Full Text
  37. ↵
    1. O'Hagan R,
    2. Chalfie M,
    3. Goodman MB
    (2005) The MEC-4 DEG/ENaC channel of Caenorhabditis elegans touch receptor neurons transduces mechanical signals. Nat Neurosci 8:43–50. doi:10.1038/nn1362 pmid:15580270
    OpenUrlCrossRefPubMed
  38. ↵
    1. Ostrow LW,
    2. Sachs F
    (2005) Mechanosensation and endothelin in astrocytes: hypothetical roles in CNS pathophysiology. Brain Res Rev 48:488–508. doi:10.1016/j.brainresrev.2004.09.005 pmid:15914254
    OpenUrlCrossRefPubMed
  39. ↵
    1. Petzold BC,
    2. Park SJ,
    3. Mazzochette EA,
    4. Goodman MB,
    5. Pruitt BL
    (2013) MEMS-based force-clamp analysis of the role of body stiffness in C. elegans touch sensation. Integr Biol (Camb) 5:853–864. doi:10.1039/c3ib20293c pmid:23598612
    OpenUrlCrossRefPubMed
  40. ↵
    1. Plaksin M,
    2. Shoham S,
    3. Kimmel E
    (2014) Intramembrane cavitation as a predictive bio-piezoelectric mechanism for ultrasonic brain stimulation. Phys Rev 4:011004.
  41. ↵
    1. Rabbitt RD,
    2. Brichta AM,
    3. Tabatabaee H,
    4. Boutros PJ,
    5. Ahn J,
    6. Della Santina CC,
    7. Poppi LA,
    8. Lim R
    (2016) Heat pulse excitability of vestibular hair cells and afferent neurons. J Neurophysiol 116:825–843. doi:10.1152/jn.00110.2016 pmid:27226448
    OpenUrlCrossRefPubMed
  42. ↵
    1. Ramot D,
    2. Johnson BE,
    3. Berry TL Jr.,
    4. Carnell L,
    5. Goodman MB
    (2008) The parallel worm tracker: a platform for measuring average speed and drug-induced paralysis in nematodes. PLoS One 3:e2208. doi:10.1371/journal.pone.0002208 pmid:18493300
    OpenUrlCrossRefPubMed
    1. Ramot D,
    2. MacInnis BL,
    3. Goodman MB
    (2008) Bidirectional temperature-sensing by a single thermosensory neuron in C. elegans. Nat Neurosci 11:908–915. doi:10.1038/nn.2157 pmid:18660808
    OpenUrlCrossRefPubMed
  43. ↵
    1. Sarvazyan AP,
    2. Rudenko OV,
    3. Nyborg WL
    (2010) Biomedical applications of radiation force of ultrasound: historical roots and physical basis. Ultrasound Med Biol 36:1379–1394. doi:10.1016/j.ultrasmedbio.2010.05.015 pmid:20800165
    OpenUrlCrossRefPubMed
  44. ↵
    1. Sassaroli E,
    2. Vykhodtseva N
    (2016) Acoustic neuromodulation from a basic science prospective. J Ther Ultrasound 4:17. doi:10.1186/s40349-016-0061-z pmid:27213044
    OpenUrlCrossRefPubMed
  45. ↵
    1. Schafer WR
    (2015) Mechanosensory molecules and circuits in C. elegans. Pflugers Arch 467:39–48. doi:10.1007/s00424-014-1574-3 pmid:25053538
    OpenUrlCrossRefPubMed
  46. ↵
    1. Schneider ER,
    2. Anderson EO,
    3. Gracheva EO,
    4. Bagriantsev SN
    (2014) Temperature sensitivity of two-pore (K2P) potassium channels. Curr Top Membr 74:113–133. doi:10.1016/B978-0-12-800181-3.00005-1 pmid:25366235
    OpenUrlCrossRefPubMed
  47. ↵
    1. Stiernagle T
    (2006) Maintenance of C. elegans. WormBook 11:1–11. doi:10.1895/wormbook.1.101.1 pmid:18050451
    OpenUrlCrossRefPubMed
  48. ↵
    1. Sukharev S,
    2. Corey DP
    (2004) Mechanosensitive channels: multiplicity of families and gating paradigms. Sci STKE 2004:re4. doi:10.1126/stke.2192004re4 pmid:14872099
    OpenUrlAbstract/FREE Full Text
  49. ↵
    1. Syeda R,
    2. Florendo MN,
    3. Cox CD,
    4. Kefauver JM,
    5. Santos JS,
    6. Martinac B,
    7. Patapoutian A
    (2016) Piezo1 channels are inherently mechanosensitive. Cell Rep 17:1739–1746. doi:10.1016/j.celrep.2016.10.033 pmid:27829145
    OpenUrlCrossRefPubMed
  50. ↵
    1. Trahey GE,
    2. Palmeri ML,
    3. Bentley RC,
    4. Nightingale KR
    (2004) Acoustic radiation force impulse imaging of the mechanical properties of arteries: in vivo and ex vivo results. Ultrasound Med Biol 30:1163–1171. doi:10.1016/j.ultrasmedbio.2004.07.022 pmid:15550320
    OpenUrlCrossRefPubMed
  51. ↵
    1. Treeby BE,
    2. Cox BT
    (2010) k-Wave: MATLAB toolbox for the simulation and reconstruction of photoacoustic wave fields. J Biomed Opt 15:021314. doi:10.1117/1.3360308 pmid:20459236
    OpenUrlCrossRefPubMed
  52. ↵
    1. Tsui PH,
    2. Wang SH,
    3. Huang CC
    (2005) In vitro effects of ultrasound with different energies on the conduction properties of neural tissue. Ultrasonics 43:560–565. doi:10.1016/j.ultras.2004.12.003 pmid:15950031
    OpenUrlCrossRefPubMed
  53. ↵
    1. Tufail Y,
    2. Matyushov A,
    3. Baldwin N,
    4. Tauchmann ML,
    5. Georges J,
    6. Yoshihiro A,
    7. Tillery SI,
    8. Tyler WJ
    (2010) Transcranial pulsed ultrasound stimulates intact brain circuits. Neuron 66:681–694. doi:10.1016/j.neuron.2010.05.008 pmid:20547127
    OpenUrlCrossRefPubMed
  54. ↵
    1. Tufail Y,
    2. Yoshihiro A,
    3. Pati S,
    4. Li MM,
    5. Tyler WJ
    (2011) Ultrasonic neuromodulation by brain stimulation with transcranial ultrasound. Nat Protoc 6:1453–1470. doi:10.1038/nprot.2011.371 pmid:21886108
    OpenUrlCrossRefPubMed
  55. ↵
    1. Tyler WJ
    (2011) Noninvasive neuromodulation with ultrasound? A continuum mechanics hypothesis. Neuroscientist 17:25–36. doi:10.1177/1073858409348066 pmid:20103504
    OpenUrlCrossRefPubMed
  56. ↵
    1. Tyler WJ
    (2012) The mechanobiology of brain function. Nat Rev Neurosci 13:867–878. doi:10.1038/nrn3383 pmid:23165263
    OpenUrlCrossRefPubMed
  57. ↵
    1. Wasserman SM,
    2. Beverly M,
    3. Bell HW,
    4. Sengupta P
    (2011) Regulation of response properties and operating range of the AFD thermosensory neurons by cGMP signaling. Curr Biol 21:353–362. doi:10.1016/j.cub.2011.01.053 pmid:21315599
    OpenUrlCrossRefPubMed
  58. ↵
    1. Way JC,
    2. Chalfie M
    (1989) The mec-3 gene of Caenorhabditis elegans requires its own product for maintained expression and is expressed in three neuronal cell types. Genes Dev 3:1823–1833. doi:10.1101/gad.3.12a.1823 pmid:2576011
    OpenUrlAbstract/FREE Full Text
  59. ↵
    1. Yan Z,
    2. Zhang W,
    3. He Y,
    4. Gorczyca D,
    5. Xiang Y,
    6. Cheng LE,
    7. Meltzer S,
    8. Jan LY,
    9. Jan YN
    (2013) Drosophila NOMPC is a mechanotransduction channel subunit for gentle-touch sensation. Nature 493:221–225. doi:10.1038/nature11685 pmid:23222543
    OpenUrlCrossRefPubMed
  60. ↵
    1. Ye PP,
    2. Brown JR,
    3. Pauly KB
    (2016) Frequency dependence of ultrasound neurostimulation in the mouse brain. Ultrasound Med Biol 42:1512–1530. doi:10.1016/j.ultrasmedbio.2016.02.012 pmid:27090861
    OpenUrlCrossRefPubMed
  61. ↵
    1. Yoo SS,
    2. Bystritsky A,
    3. Lee JH,
    4. Zhang Y,
    5. Fischer K,
    6. Min BK,
    7. McDannold NJ,
    8. Pascual-Leone A,
    9. Jolesz FA
    (2011) Focused ultrasound modulates region-specific brain activity. Neuroimage 56:1267–1275. doi:10.1016/j.neuroimage.2011.02.058 pmid:21354315
    OpenUrlCrossRefPubMed
  62. ↵
    1. Zhou W,
    2. Wang J,
    3. Wang K,
    4. Huang B,
    5. Niu L,
    6. Li F,
    7. Cai F,
    8. Chen Y,
    9. Liu X,
    10. Zhang X,
    11. Cheng H,
    12. Kang L,
    13. Meng L,
    14. Zheng H
    (2017) Ultrasound neuro-modulation chip: activation of sensory neurons in Caenorhabditis elegans by surface acoustic waves. Lab Chip 17:1725–1731. doi:10.1039/C7LC00163K pmid:28447086
    OpenUrlCrossRefPubMed
Back to top

In this issue

The Journal of Neuroscience: 38 (12)
Journal of Neuroscience
Vol. 38, Issue 12
21 Mar 2018
  • Table of Contents
  • Table of Contents (PDF)
  • About the Cover
  • Index by author
  • Advertising (PDF)
  • Ed Board (PDF)
Email

Thank you for sharing this Journal of Neuroscience article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Ultrasound Elicits Behavioral Responses through Mechanical Effects on Neurons and Ion Channels in a Simple Nervous System
(Your Name) has forwarded a page to you from Journal of Neuroscience
(Your Name) thought you would be interested in this article in Journal of Neuroscience.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Print
View Full Page PDF
Citation Tools
Ultrasound Elicits Behavioral Responses through Mechanical Effects on Neurons and Ion Channels in a Simple Nervous System
Jan Kubanek, Poojan Shukla, Alakananda Das, Stephen A. Baccus, Miriam B. Goodman
Journal of Neuroscience 21 March 2018, 38 (12) 3081-3091; DOI: 10.1523/JNEUROSCI.1458-17.2018

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Respond to this article
Request Permissions
Share
Ultrasound Elicits Behavioral Responses through Mechanical Effects on Neurons and Ion Channels in a Simple Nervous System
Jan Kubanek, Poojan Shukla, Alakananda Das, Stephen A. Baccus, Miriam B. Goodman
Journal of Neuroscience 21 March 2018, 38 (12) 3081-3091; DOI: 10.1523/JNEUROSCI.1458-17.2018
Reddit logo Twitter logo Facebook logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Jump to section

  • Article
    • Abstract
    • Introduction
    • Materials and Methods
    • Results
    • Discussion
    • Footnotes
    • References
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF

Keywords

  • Caenorhabditis elegans
  • mechanosensation
  • mechanosensitive ion channels
  • neurostimulation
  • thermosensation
  • ultrasound

Responses to this article

Respond to this article

Jump to comment:

No eLetters have been published for this article.

Related Articles

Cited By...

More in this TOC Section

Research Articles

  • Musical training facilitates exogenous temporal attention via delta phase entrainment within a sensorimotor network
  • Microglial Cytokines Mediate Plasticity Induced by 10 Hz Repetitive Magnetic Stimulation
  • Subgenual and hippocampal pathways in amygdala are set to balance affect and context processing
Show more Research Articles

Cellular/Molecular

  • Transferrin receptor is necessary for proper oligodendrocyte iron homeostasis and development
  • Drosophila tachykininergic neurons modulate the activity of two groups of receptor-expressing neurons to regulate aggressive tone
  • Coupling between the Stereocilia of Rat Sensory Inner-Hair-Cell Hair Bundles Is Weak, Shaping Their Sensitivity to Stimulation
Show more Cellular/Molecular
  • Home
  • Alerts
  • Visit Society for Neuroscience on Facebook
  • Follow Society for Neuroscience on Twitter
  • Follow Society for Neuroscience on LinkedIn
  • Visit Society for Neuroscience on Youtube
  • Follow our RSS feeds

Content

  • Early Release
  • Current Issue
  • Issue Archive
  • Collections

Information

  • For Authors
  • For Advertisers
  • For the Media
  • For Subscribers

About

  • About the Journal
  • Editorial Board
  • Privacy Policy
  • Contact
(JNeurosci logo)
(SfN logo)

Copyright © 2023 by the Society for Neuroscience.
JNeurosci Online ISSN: 1529-2401

The ideas and opinions expressed in JNeurosci do not necessarily reflect those of SfN or the JNeurosci Editorial Board. Publication of an advertisement or other product mention in JNeurosci should not be construed as an endorsement of the manufacturer’s claims. SfN does not assume any responsibility for any injury and/or damage to persons or property arising from or related to any use of any material contained in JNeurosci.