Touch deforms, or strains, the skin beyond the immediate point of contact. The spatiotemporal nature of the touch-induced strain fields depend on the mechanical properties of the skin and the tissues below. Somatosensory neurons that sense touch branch out within the skin and rely on a set of mechano-electrical transduction channels distributed within their dendrites to detect mechanical stimuli. Here, we sought to understand how tissue mechanics shape touch-induced mechanical strain across the skin over time and how individual channels located in different regions of the strain field contribute to the overall touch response. We leveraged Caenorhabditis elegans’ touch receptor neurons as a simple model amenable to in vivo whole-cell patch-clamp recording and an integrated experimental-computational approach to dissect the mechanisms underlying the spatial and temporal dynamics we observed. Consistent with the idea that strain is produced at a distance, we show that delivering strong stimuli outside the anatomical extent of the neuron is sufficient to evoke MRCs. The amplitude and kinetics of the MRCs depended on both stimulus displacement and speed. Finally, we found that the main factor responsible for touch sensitivity is the recruitment of progressively more distant channels by stronger stimuli, rather than modulation of channel open probability. This principle may generalize to somatosensory neurons with more complex morphologies.

## Introduction

Sensory receptor neurons are classified by their receptive fields—the region of their sensory space in which a stimulus elicits a response, such as a spot in the visual field or a frequency of sound. The receptive field of a somatosensory neuron is a combined function of its dendritic arbor and the mechanical strain, or deformation, induced by touch. Mechanical strain decreases with distance from the stimulus; the extent of the strain field and the speed of propagation depend on the intensity of the stimulus and the material properties of the neurons and skin being touched. With light, a single photon can activate only a single receptor protein at a time, even if many are in close proximity. With touch, a single indentation will simultaneously affect all eligible receptors that fall within the strain field, to a degree dependent on their distance from the point of stimulation. In the case of somatosensation, this means that multiple channels at spatially distinct sites within a neuron, or even multiple neurons, are liable to be activated to varying degrees by the same stimulus.

Intricate anatomical structures dedicated to mechanosensation are present in worms, flies, mice, cats, and humans. From the chordotonal organs and campaniform receptors of insects to Merkel cell touch domes and low-threshold mechanoreceptors arrayed around hair follicles in mammals, these structures place mechanosensitive proteins in specific arrangements around features of the skin or other mechanosensory organs (Katta et al., 2015). Increasing anatomical and functional data have provided insights into how these structures work. In many cases, the complexity of skin has made it difficult to understand the strain field. Although useful models have been created (Lesniak and Gerling, 2009; Quindlen et al., 2015; Quindlen-Hotek and Barocas, 2018; Sanzeni et al., 2019), few studies build on these models to generate an integrated understanding of how mechanical strain is detected by mechanosensitive proteins, including ion channels, distributed within somatosensory neurons.

Here, we probed how mechano-electrical transduction (MeT) channel distribution interacts with the strain field of a mechanical stimulus in a Caenorhabditis elegans touch receptor neuron (TRN). This system is a useful model because the factors governing the strain field are straightforward and the arrangement of channels is simple and well characterized. Furthermore, we can directly measure the mechanoreceptor currents (MRCs). The six C. elegans TRNs (ALMR/L, PLMR/L, AVM, and PVM) extend long, unusually straight neurites (Krieg et al., 2017) that are embedded in the epidermal cells that form the worm’s skin or cuticle (Chalfie and Thomson, 1979; Chalfie and Sulston, 1981). External mechanical loads evoke MRCs that are carried by MEC-4-dependent MeT channels and are activated at both the application and removal of mechanical stimuli (O’Hagan et al., 2005). The MEC-4 protein belongs to a large superfamily of non–voltage-gated ion channels conserved in animals but absent from microbes and plants (Katta et al., 2015). MEC-4 proteins localize to discrete puncta that are arrayed along the entire length of TRN neurites (Chelur et al., 2002; Emtage et al., 2004; Cueva et al., 2007). These anatomical properties make it easier to interpret how the geometry of the strain field depends on the size and speed of a given stimulus and to determine which subset of MeT channels is likely to be affected by each stimulus.

Because MRC amplitude depends more on indentation than it does on applied force (Eastwood et al., 2015), we developed and deployed a stimulator system enabling fast indentation and concurrent optical monitoring of probe movement. We used this system to measure MRCs evoked by mechanical stimuli delivered inside and outside of the anatomical receptive field of the ALM neuron and as a function of indentation depth and speed. By combining these experimental results with simulations (Sanzeni et al., 2019), we identified a biophysical mechanism linking indentation to MRCs. Specifically, we show that MRC size and time course depend on both the modulation of open probability and the recruitment of distal channels in a model somatosensory neuron.

## Materials and methods

### Nematode strains

We used three strains of transgenic C. elegans nematodes: TU2769 uIs31[mec-17p::gfp] III (O’Hagan et al., 2005) for electrophysiology, GN865 uIs31[mec-17p::GFP] III; kaIs12[col-19::GFP] for imaging cuticle indentation, and GN753 pgSi116[mec-17p::mNeonGreen::3xFLAG::MEC-4::tbb-2 3′UTR] II for visualizing MEC-4 puncta. The uIs31 transgene is a TRN-specific GFP marker that enables us to conduct in vivo recordings in TRNs. The kaIs12 transgene encodes a GFP fusion to a collagen that labels cuticular annuli (Thein et al., 2003). The pgSi116 transgene expresses MEC-4 tagged with mNeonGreen in the TRNs to visualize MEC-4 puncta. We grew animals on Escherichia coli strain OP50 and used well-fed subjects as late-L4 larvae or young adults. On average, the animals used for electrophysiology had a body length of 1,298 ± 93 µm (mean ± SD, n = 53; range of 1,064–1,547) and diameter of 39 ± 4 µm (range of 32–50) at the terminal bulb of the pharynx.

### Imaging cuticle deformation

We immobilized GN865 worms using 2% agarose pads with WormGlu (GluStitch); subjects were either left intact or dissected as described previously (O’Hagan et al., 2005; Eastwood et al., 2015). We used an Orca Flash 4.0LTv2 camera (Hamamatsu) controlled by μManager (Edelstein et al., 2010) on an upright microscope (E600FN; Nikon) under a 60× objective for imaging.

### Electrophysiology

We recorded from the ALMR and ALML TRNs in TU2769 worms. Due to geometric constraints imposed by our stimulator system, we recorded from ALMR when stimulating anterior to the cell body and from ALML when stimulating posterior to the cell body. These two neurons are bilaterally symmetric and exhibit no detectable variation in voltage- or touch-evoked currents (not shown). The extracellular solution contained (in mM): NaCl (145), KCl (5), MgCl2 (5), CaCl2 (1), and Na-HEPES (10), adjusted to a pH of 7.2 with NaOH. Before using the solution, we added 20 mM of d-glucose, which brought the osmolarity to ∼325 mOsm. The intracellular solution contained (in mM): K-gluconate (125), KCl (18), NaCl (4), MgCl2 (1), CaCl2 (0.6), K-HEPES (10), and K2EGTA (10), adjusted to a pH of 7.2 with KOH. Before use, we added 1 mM of sulforhodamine 101 (Invitrogen) to help visualize whether the neuron was being successfully recorded.

We used an EPC-10 USB amplifier controlled by Patchmaster software (version 2x90.1; HEKA/Harvard Biosciences) to set the membrane voltage and control the mechanical stimulator. We corrected the membrane voltage for the liquid-junction potential (−14 mV) between the extracellular and intracellular solutions; we also corrected for errors that may have resulted from an imperfect space clamp based on stimulus distance, as described previously (Goodman et al., 1998; O’Hagan et al., 2005). We filtered the analogue data at 2.9 kHz and digitized it at 10 kHz.

### Mechanical stimulation

Some of the previous electrophysiological studies on C. elegans mechanoreceptor neurons have used open-loop systems with a piezoelectric bimorph (O’Hagan et al., 2005; Bounoutas et al., 2009; Arnadóttir et al., 2011; Geffeney et al., 2011; Chen and Chalfie, 2015) or a piezoelectric stack with no independent measurement of stimulator motion (Kang et al., 2010). We previously employed a slow closed-loop system driven by a piezoelectric stack with a piezoresistive cantilever for force detection (Eastwood et al., 2015). Here, we used an open-loop system adapted from the piezoelectric stack system with a photodiode motion detector described by Peng and Ricci (2016) in their study on hair cells. This system enables faster stimulation than the force-clamp system and allows us to measure the time course of stimulation.

Using marine epoxy (Loctite), we attached an open-loop piezoelectric stack actuator (PAS-005; ThorLabs, 20 µm travel) to a tungsten rod (length = 8 in.) to damp vibration. We mounted the rod-mounted actuator on a micromanipulator for positioning (MP-225; Sutter) and controlled the piezoelectric stack via an analogue voltage output on a HEKA EPC-10 PLUS amplifier and Patchmaster software. We filtered this analogue signal at 2.5 kHz for steps or 5 kHz for ramps and sines on an 8-pole Bessel filter (LPF-8; Warner Instruments) and used a high-voltage, high-current amplifier (Peng and Ricci, 2016) to achieve a signal between 0 and 75 V. The stack was biased with a starting offset of 3–4 µm, and the largest displacement we used was 3–4 µm short of the 20-µm total travel limit. This protocol ensured that stack motion was linearly related to the analogue voltage signal.

Beads glued to force-clamp cantilevers create a defined and reproducible contact surface that is easy to model (Petzold et al., 2013; Eastwood et al., 2015). We adapted this technique for a system relying on stiff glass, rather than silicon probes, and for axial stimulus motion. We pulled and polished borosilicate glass pipettes (Sutter; BF150-86-10) to a tip diameter that was half or three-fourths of the diameter of the bead that we intended to use. We attached black polyethylene beads (20–24-µm diameter; BKPMS-1.2; Cospheric) to the glass pipettes using UV-curable glue (Loctite 352; Henkel). The bead-attached pipettes were then waxed in place to a 3-D-printed acrylic pipette holder (custom design available from https://3dprint.nih.gov as model ID 3DPX-010770) and epoxied to a steel tip (PAA001; ThorLabs) that was mounted on the piezo stack.

Step protocols had a 250- or 300-ms hold at the commanded indentation and used an interstimulus interval (ISI) of 1 s. Trapezoidal protocols had a 300-ms hold and 2-s ISI. Sinusoidal stimulus profiles consisted of a 1-s sinusoidal stimulus epoch flanked by 300-ms steps and 2-s ISI.

### Motion detection

To monitor the motion of the stimulus probe, we modified the system reported by Peng and Ricci (2016) to detect larger movements by using the SPOT-2D segmented photodiode (OSI Optoelectronics) and a higher resistance in the differential amplifier circuit. The photodiode was mounted in an XY translator on top of a rotation stage (ST1XY-D, LCP02R; ThorLabs) such that the photodiode gap was aligned perpendicular to the direction of the probe motion. This apparatus was fixed above a secondary camera port on the microscope (Eclipse E600FN; Nikon) with no additional magnification.

Before obtaining an on-cell, high-resistance seal for patch clamp, we aligned the front edge of the bead under the 60× objective with the highest clearly visible edge of the worm’s cuticle. After achieving whole-cell access, we moved the microscope to visually align the image of the bead with the outer edge of the photodiode. As the bead moved, the output from the photodiode circuit was read directly by Patchmaster, and we used this signal to check that resonance during steps and relative attenuation of sine amplitude at high frequencies were small.

At the end of each recording, we calibrated the image motion to actual motion. We drove the XY translator in steps of known distances in a direction that was opposite to the probe motion using a stepper motor actuator (ZFS-06; ThorLabs) controlled by Kinesis software (ThorLabs). We took these steps once with the probe at the edge of the worm and once with only the worm. Distances were measured relative to this starting position. The worm-only trace was subtracted from the probe + worm trace to obtain the relationship between the motion of the probe image and the known movement of the photodiode. We corrected stimulus traces from the recording using a piecewise cubic polynomial interpolation method (pchip; Matlab R2015b; MathWorks) of the mean voltage during calibration steps versus the commanded photodiode motion at each step.

### Data analysis

We measured whole-cell capacitance and series resistance as previously described (Goodman et al., 1998). We conducted data analysis using Matlab (data import and analysis functions are available online at https://github.com/wormsenseLab/Matlab-PatchMaster/tree/vJGP and Igor Pro 6.37 (Wavemetrics).

Our analysis included only recordings with a holding current less than −10 pA at −60 mV and a series resistance 250 MΩ. As these experiments were conducted at constant voltage, the voltage errors due to uncompensated series resistance were negligible, and we did not correct the reported data for these errors. Where noted, we corrected for voltage attenuation due to space-clamp errors as previously described (O’Hagan et al., 2005). Briefly, we computed the theoretical capacitance for each neuron based on body size and compared with to the measured capacitance to estimate the ratio of membrane resistance to axial resistance for ALM. This allowed us to calculate the length constant for each neuron and estimate the actual voltage at the center of the stimulus, which we used as the average channel location for the purpose of accounting for space-clamp error. We calculated total charge, or the area under the MRC curve, from stimulus start to 150 ms after stimulus end and subtracted the charge at baseline for an equivalent time period.

### MEC-4 puncta imaging and data analysis

We mounted young adult GN753 pgSi116 mNeonGreen::3xFLAG::MEC-4 transgenic animals on grooved 5% agarose pads dissolved in M9 buffer. Grooved pads were prepared by casting agarose on the surface of a vinyl record (Rivera Gomez and Schvarzstein, 2018). Animals were immobilized using levamisole (5 mM) dissolved in M9 buffer. Micrographs were acquired with an automated inverted epifluorescence microscope (Keyence BZ-X800) equipped with software and hardware to generate stitched images of entire ALM neurons with a 40× objective (Nikon plan apo 40×/1.0 oil).

The stitched images were preprocessed using ImageJ (National Institutes of Health), where we manually traced a line along the neurite from the cell body to the distal tip to create a 20-pixel-wide straightened image of the neurite. We then used these straightened neurite images to analyze MEC-4 puncta distribution using a custom Python script (https://github.com/wormsenseLab/Puncta_analysis/tree/vJGP%2Blicense). Briefly, we averaged the pixel intensities in the middle rows to calculate the raw neurite fluorescence, subtracted the average background fluorescence at the top and bottom, divided by the background to normalize across animals, and smoothed with a 2-pixel-maximum filter to reduce noise. Peaks in the fluorescence traces were detected automatically (find_peaks, scipy library). Interpunctum intervals (IPIs) were calculated by measuring the distance between adjacent peaks along the neurite; the location of each IPI is calculated as the distance of the midpoint between the adjacent peaks from the cell body.

### Simulations

For simulations, we used the model described in Sanzeni et al. (2019) and model parameters derived by fitting the dataset from Eastwood et al. (2015) and the Falcon closed-loop stimulation system (Table 1). We took this approach for the following reasons. First, the previous dataset directly linked force, indentation, and current, enabling the inference of model parameters. By contrast, the current dataset was collected using a mechanical stimulation system that enables high-speed indentation at the expense of this direct link. Second, single-channel parameters are, to a first approximation, preserved across worms and independent of internal pressure (Sanzeni et al., 2019). Consequently, we use single-channel parameters derived from the dataset of Eastwood et al. (2015), which allows us to infer internal pressure, and show model predictions obtained using either a typical value of 4-kPa internal pressure (Sanzeni et al., 2019) or different values of internal pressure (e.g., Fig. 3). The stimulator was modeled as a rigid sphere with a 20-µm diameter, in accord with the stimulation apparatus used in this study.

As in our recent computational study (Sanzeni et al., 2019), we built a finite element mechanical model of the worm as a pressurized, thin-shelled cylinder. A three-state (closed, subconducting, open) gating model of the channels allowed us to estimate the probabilities for individual channels to open. We opted to include the subconducting state because single-channel recordings of heterologously expressed MEC-4-dependent channels demonstrate a subconducting state (Brown et al., 2007). Additionally, including the subconducting state improved the match between simulated and experimental data (Eastwood et al., 2015). The total current flowing along the TRN is the sum of the currents flowing through the channels arrayed along the neuron that spans the length of the worm. In simulations, we used a fixed IPI of 1.4 µm and assumed that each punctum contained a single functional MeT channel. These choices are supported by empirical studies of endogenous channel distribution and estimates of the number of functional channels per punctum (Cueva et al., 2007).

The strain energy density (SED) needed to activate a single channel was estimated as follows. We considered a single channel located 1 µm below the cuticle surface and 1 µm from the stimulus center and computed its current response to step stimuli of different displacements. For any given displacement w, we also computed the SED produced at the channel position x (Audoly and Pomeau, 2010) as

$SED=∫0wσij(x)δεij(x),$

where $σij$ and $δεij$ are the stress tensor and the change in strain tensor, respectively. Using this approach and assuming that reliable neuronal activation would be achieved following stimuli sufficient to achieve a half-maximal MRC amplitude, we estimated the SED needed to activate a TRN.

### Online supplementary material

The model developed by Sanzeni et al. (2019) predicts the ionic current flowing through MeT channels distributed along the neurite. The amplitude of this current depends nonlinearly on indentation depth and channel distance from the center of the stimulating bead. Fig. S1 shows an example of this analysis.

## Results

MRCs in C. elegans TRNs are carried by MeT channels distributed in puncta along the length of each sensory neurite (Chelur et al., 2002; Emtage et al., 2004; Cueva et al., 2007). For many other types of channels, an increase in macroscopic current is thought to reflect an increase in the open probability of a population of channels that receive the same stimulus at the same time. For MeT channels, however, an increase in the macroscopic current cannot be interpreted solely as an increase in open probability. As local indentation deforms tissue at distal sites (Fig. 1, A–C) and in a manner that depends on indentation depth (Elmi et al., 2017; Sanzeni et al., 2019), stimulus intensity governs how many channels are reachable by the stimulus and the probability that a given channel will open (Fig. 1 A). Similarly, we expect that stimulus speed influences the evolution of the strain field, thereby introducing an additional factor that affects the synchrony or phase of MeT channel activation. As such, the spatial distribution of these MeT channels can affect how they respond to stimuli that vary across space and time.

### MRC amplitude depends on stimulus probe size, displacement, and speed

To enable a thorough investigation of how indentation depth and speed influence activation of MeT channels in their native environment, we built a high-speed mechanical stimulator (Materials and methods). While recording MRCs, we positioned this stimulator perpendicular to the worm’s body and used the stimulator to push a stiff glass probe carrying a large (∼22-µm diameter) polyethylene bead into living, wild-type worms (Fig. 1, B and C). We applied displacements up to 12 µm, which is ∼30% of the average body width (39 ± 4 µm, mean ± SD, n = 53). Unexpectedly, the peak MRCs we report here are approximately fivefold larger than those previously reported by us (Fig. 1 E). Here, we recorded from ALM, rather than PLM. We pushed the bead perpendicularly, rather than tangentially into the worm (O’Hagan et al., 2005; Arnadóttir et al., 2011; Chen and Chalfie, 2015; Eastwood et al., 2015) and used larger beads to deliver the mechanical stimulus (Eastwood et al., 2015). We also used high-speed stimulation with minimal resonant vibration.

Intrinsically larger responses in ALM versus PLM might account for some of the difference (Chen and Chalfie, 2015), but ALM currents measured with a force-clamped stimulus (Eastwood et al., 2015) were still significantly smaller than what we measured here (Fig. 1 E). Thus, the discrepancy in peak MRC size is likely to reflect differences in the stimulus paradigm or worm preparation between studies: the use of a larger bead to contact the animal’s skin, high-speed stimulation, variations in body mechanics, or a combination of these factors. For a given indentation, our model (Sanzeni et al., 2019) predicts that larger beads and stiffer animals will generate larger strain fields, which would in turn increase the number of channels available for activation. In addition to these factors, the model predicts that current amplitude increases with stimulation speed (Sanzeni et al., 2019). The indentation steps we applied here are at least an order of magnitude faster than those of the Falcon system (Eastwood et al., 2015; Sanzeni et al., 2019). Although dissection of cell bodies for recording creates variation in body stiffness that affects peak current amplitude (Sanzeni et al., 2019), this is unlikely to account for the 10-fold change in mean MRC amplitude (Fig. 1 E).

Another factor affecting MRC amplitude is the expression of MEC-4, a key pore-forming subunit of the native MeT channel (O’Hagan et al., 2005). Whether visualized by immunofluorescence labeling of fixed samples (Cueva et al., 2007; Chen and Chalfie, 2015) or by direct observation of transgenic animals expressing a fluorescent protein fused to MEC-4 (Chelur et al., 2002; Emtage et al., 2004; Arnadóttir et al., 2011; Petzold et al., 2013; Vásquez et al., 2014), MEC-4 localizes to discrete puncta along TRN neurites. In prior studies (Cueva et al., 2007; Arnadóttir et al., 2011; Petzold et al., 2013; Vásquez et al., 2014; Chen and Chalfie, 2015), adjacent puncta were separated by 1.3–3.8 µm, on average, and these average values were obtained by analyzing arbitrary neurite fragments. To gain a comprehensive and global view of the distribution of MEC-4 protein, we tagged MEC-4 with mNeonGreen, an exceptionally bright fluorescent protein suitable for in vivo applications (Hostettler et al., 2017), and developed an image processing workflow to analyze IPIs in full-length TRNs (Materials and methods). With this approach, we estimate that a single ALM neuron harbors 113 puncta (5,783 puncta detected in 51 neurons), on average. Additionally, the distance between neighboring mNeonGreen::MEC-4 puncta increases with distance from the ALM cell body (Fig. 2).

Pooling IPIs without regard to position along the neurite yields a log-normal distribution (Fig. 2 C), reminiscent of the distribution first reported by Cueva et al. (2007). This finding suggests an explanation for the variation in prior reports of average IPIs: namely, that datasets sampling fragments closer to the cell body would yield smaller IPIs than those sampling fragments farther away. We tested this idea computationally, by sampling 100-µm fragments centered at distances of 50 to 100 µm from the ALM cell body, matching the delivery of mechanical stimuli in this study (Fig. 2 D, green). As expected, this resulted in a narrower distribution with a shorter average IPI than when we sampled 100-µm fragments at any distance (Fig. 2 D, blue). Thus, mechanical stimuli delivered closer to the cell body might be expected to evoke larger MRCs than those delivered farther away, as seen in the PLM neuron (Chen and Chalfie, 2015). This factor is unlikely to account for the 10-fold change in mean MRC amplitude (Fig. 1 E), however, since the present and prior (O’Hagan et al., 2005; Eastwood et al., 2015) studies also relied on stimuli delivered ∼100 µm away from the cell body. Based on these considerations, we infer that the probe size and stimulus speed are the dominant factors that account for the larger currents.

### Peak MRCs increase with stimulus size but do not saturate

Indentation depth also affects peak current amplitude and total charge transferred (Fig. 3, A and B, left). Increasing the indentation depth has two effects, as demonstrated in our computational model linking tissue mechanics to MeT-channel activation in living animals (Sanzeni et al., 2019). The first effect is an expansion of the strain field, which increases the number of channels available for activation. The second is an increase in the probability that a given MeT channel is activated. This model also shows that the channels closest to the point of stimulation, where the deformation is greatest, are more likely to open than distant channels (Fig. S1). These two effects combine to determine MRC peak amplitude as well as the total charge transferred and result in nonlinear, nonsaturating current–indentation curves (Fig. 3, A and B, right).

Although we expect variations in body stiffness and internal pressure due to natural variation and variation in the dissection procedure (Eastwood et al., 2015), these parameters are largely outside of experimental control. To circumvent this technical limitation, we used simulations based on the model of Sanzeni et al. (2019) to systematically explore the influence of variations in pressure (Fig. 3, A and B, right). With the exception of internal pressure, the free parameters (Table 1) of the simulation were derived by fitting a dataset collected with a different stimulation system (Eastwood et al., 2015). Higher internal pressures correspond to stiffer animals and lower pressures correspond to softer ones. The simulated current amplitude for a 10-µm indentation is higher in stiff animals compared with softer ones (Fig. 3 A). Henceforth, the simulations rely on the parameters in Table 1 and an internal pressure of 4 kPa, which reflects an intermediate stiffness.

To test the effect of body stiffness in vivo, we obtained two high-quality recordings under the minimally dissected (stiff) regimen described in Eastwood et al. (2015). While these currents were larger than the average, they fell within the observed range for dissected preparations in this study (Fig. 3, A and B, dark gray). This may be due to inadvertent selection wherein we were able to successfully isolate cell bodies only from the least pressurized animals, but we cannot distinguish between these possibilities in the present dataset. Alternatively, the factors that we believe led to larger currents in this study (larger bead size and faster stimulation) may have minimized the effect of variations in stiffness.

In simulations, we used a fixed IPI of 1.4 µm, a value drawn from immunofluorescence labeling of native protein (Cueva et al., 2007). To determine how this choice affected our results, we repeated simulations with values between 1.4 and 4 µm. For each value, we determined single-channel model parameters from the best fits to recorded neural responses as discussed in Materials and methods (Table 1). We found that the peak current decreased slightly with the IPI (Fig. 3 C, left), but that the shape of the relationship between displacement and normalized peak current was unaffected (Fig. 3 C, right). These simulations show that, while the quantitative aspect of our results depends on the IPI chosen, the scaling relations found are robust to the specific choice. In the following plots, we report simulations performed with an IPI of 1.4 µm; by using the smallest value within the physiologically reported range, we avoided position-related artifacts in the responses to small displacements.

In all prior studies of MRCs in C. elegans TRNs, current–indentation curves for MeT channels have been treated as a saturating function of stimulus size. These experimental and computational results indicate that such a treatment neglects the influence of tissue mechanics as well as the delocalized nature of MeT channels. Furthermore, our results show that the shape of the current–indentation curve is independent of IPI but varies with other factors, such as body stiffness, that are extrinsic to the MeT channels themselves.

Using the Sanzeni et al. (2019) model, we further examined the effect of displacement on current by separating the effect of indentation on open probability from the effect of indentation on increased channel recruitment via an expansion of the strain field. In this model, the strength of the stimulus received by each channel, and hence its predicted open probability, depends on its location relative to the stimulus. Based on single-channel recording that revealed a robust subconductance state in MEC-4 channels (Brown et al., 2008), we included a subconductance state between closed and open states in the model and derived the subconductance current is relative to the open-channel current io from the fit (Table 1). Subconductance and fully open states carry 0.35 and 1.6 pA of inward current, respectively, with the primary parameter set. (Subconductance current ranged from 0.35 to 0.84 pA when fitted with different IPIs [Table 1].) These values are consistent with both single-channel recordings of channels expressed in heterologous cells (Brown et al., 2008) and estimates derived from nonstationary noise analysis (O’Hagan et al., 2005).

As a concise method of describing the single-channel response, we computed the mean of the peak current carried by each channel, regardless of whether we considered only channels that reach the open state (O) or also included channels in a subconductance state (O+SC). Based on this measurement, we observed that the mean current per channel rose with small displacements but then held steady (Fig. 3 D). The inclusion of subconducting channels more closely reflects the behavior of both current and charge. As such, all subsequent plots of this type reflect the O+SC condition. In contrast, the number of open channels continued to increase with increased displacement, behaving similarly to the current (Fig. 3 E). The discontinuity in the number of channels that occurs near 5-µm displacement is due to a nonlinear recruitment of a population of distant channels with low open probability (Fig. S1), which in turn decreases the average current per channel (Fig. 3 D). Collectively, these results suggest that channel recruitment dominates the overall current–indentation relationship, while increases in open probability and transitions from subconducting to open states steepen the curve, especially for small indentations.

### MeT channels are activated by distant stimuli

Our results show that indentation at one point affects multiple channels along the neurite. Thus, we expect that MRC amplitude depends on the overlap between the strain field and the neurite (Fig. 4 A). We explored this concept experimentally by holding the stimulus size constant and varying the distance between the ALM cell body and the stimulator in both the anterior and posterior directions (Fig. 4 B). For recordings obtained while stimulating anterior to the cell body, we corrected for errors arising from the fact that the TRNs are not isopotential (O’Hagan et al., 2005; Materials and methods). Although we did not detect a systematic relationship between current amplitude and distance from the cell body in the anterior direction (Fig. 4 B, blue), the high variation in current amplitudes across individual recordings could mask such a relationship. Additional experiments involving measurements at several positions during a single recording may be required to resolve this uncertainty.

There are no channels from ALM posterior to its cell body. Despite this fact, stimuli delivered posterior to ALM can evoke calcium transients in the ALM cell body (Suzuki et al., 2003), which are likely driven by MRCs (Fig. 4 B). On average, these currents had smaller amplitudes than currents elicited by stimuli of an equivalent magnitude anterior to the cell body (Figs. 4 B and 3 C). Although the variation in current amplitude was lower than the variation we observed for anterior stimulation, there was no obvious relationship between current amplitude and distance within 100 µm of the cell body in the posterior direction. For more distant stimuli, however, we could not detect currents even for large displacements (n = 4).

Consistent with the reduced overlap between the strain field and the neurite for posterior stimulation, larger displacements were necessary for eliciting currents under this condition (Figs. 4 C and 5, A and B). Peak current–displacement curves were shifted to the right relative to curves observed for anterior stimulation, such that the posterior stimulation evoked consistently detectable MRCs only for indentations 4 µm. Simulations captured this shift in peak current–indentation curves (Fig. 5 A, right) and suggested that the size of the shift is proportional to the distance from the cell body. We also plotted the total charge (the area under the MRC curve) against displacement to account for differences in the timing of channel opening. This reduced the difference in operating range in our experimental data (Fig. 5 B, left), but not in the simulations (Fig. 5 B, right). This suggests that channels activated during posterior stimulation are activated in a less synchronous manner than channels activated during anterior stimulation. For both measures, larger stimuli are necessary for activating channels in the simulations.

Given the discrepancy between peak current and total charge, we looked at how the activation and decay rates varied with displacement in anterior and posterior recordings. Passive conduction over a long distance can broaden the time course of electrical events measured at the cell body (Rall, 1967; Bekkers and Stevens, 1996). Therefore, for both anterior and posterior recordings, we limited this analysis to stimuli within 90 µm of the cell body, less than the average length constant λ = 96 ± 3 µm (mean ± SEM, n = 83) calculated for our recordings.

We estimated activation rates by measuring the time required for the current to rise to half the peak value from the start of the stimulus, and decay rates from the time to peak to when the current dropped to half the peak value. Activation rate increased with stimulus size and was similar for anterior and posterior stimulation (Fig. 5 C, left). We failed to detect any difference in latency. The decay rate for anterior, but not posterior, stimulation increased slowly with displacement (Fig. 5 D, left). This observation is consistent with space clamp error for channels beyond the point of stimulation (the point at which we corrected for voltage attenuation): larger stimuli affect channels that are much farther away from the cell body, and the contribution of these channels to the total current measured at the cell body is slowed and decreased by passive conduction. Indeed, when we included recordings in which the stimulator was placed even farther anterior to the cell body, we found that both the activation and decay were slower (not depicted). Channels activated by posterior stimulation are located near the cell body, where this error is low, but their position at the edge of the strain field means they are activated only weakly. The simulation does not consider space clamp error, but predicts faster decay with posterior stimulation because tissue experiencing low strain needs to move only a short distance to relax. Although the smaller shift of total charge is not entirely explained by this finding, simulations of activation and decay rates agree qualitatively with the experimental results in at least one other way: the farther posterior the stimulator is placed, the larger the displacement required to reach any channels in the neurite (Fig. 5, C and D, right). The simulation predicts slower activation times for posterior stimulation because strain takes longer to reach distant channels than proximal channels. However, we could not clearly detect such a difference in our experiments with this fast step protocol.

As noted above, the macroscopic currents we record depend on (a) stimulus-dependent changes in a channel’s open probability and (b) the number of channels exposed to each stimulus. Using average peak single-channel current as a proxy for open probability, our simulations predict that channel open probability saturates at smaller displacements for anterior stimuli than for posterior stimuli (Fig. 5 E). Also, the fact that larger stimuli reach more channels than smaller stimuli is the dominant factor controlling peak current–displacement curves for both anterior and posterior stimuli (Fig. 5 A).

### MRC size and activation rate depend on stimulus speed and direction

The peak amplitude of MRCs in C. elegans TRNs increases with the speed of indentation (Figs. 6 A and 7 A, left). This could reflect speed-dependent activation of individual channels. It could also reflect a speed-dependent increase in synchrony, or in-phase activation, of a population of channels. To differentiate between these possibilities, we analyzed peak current and charge across a wide range of speeds and with fine resolution. Unlike peak MRC size, charge accounts for all channel openings regardless of timing. It follows, therefore, that if the speed dependence of peak MRC amplitude was entirely caused by variations in the phase of channel opening, then the charge would be independent of speed. Contrary to this prediction, we found that charge was speed dependent in both experiments and simulations (Fig. 7 B). Collectively, these findings indicate that variations in the phase of channel activation are not sufficient to account for the speed dependence of MRC generation.

In contrast to displacement sensitivity, measured current and charge saturate at high speed, and simulated current and charge clearly saturate at speeds 1,000 µm/s (Fig. 7, A and B, left). This finding suggests that indentation depth governs the number of channels eligible for recruitment. To estimate the relative contributions of channel recruitment and gating to speed dependence, we turned to simulations, which predicted a similar dependence on speed as with displacement (Fig. 7, A–D, right). We found that both peak channel current and the number of open channels saturate at high speeds (Fig. 7, E and F) and that the number of open channels increases with stimulus speed (Fig. 7 F).

The saturating number of channels depends on the value of $g0$, a model parameter controlling the effect of linker elongation on channel opening (Fig. 7 F). Smaller parameter values render the channels more sensitive to linker elongation, allowing recruitment of more channels. Thus, for a fixed final displacement, there seems to be a maximal number of channels that can be opened, and this number is expected to be inversely proportional to the force needed to activate single channels. Even with $g0$ and the final displacement held constant, channel recruitment increases with stimulus speed; this results from the larger forces generated by quicker changes in strain. Thus, MRC amplitude increases with speed because more distant channels are activated in response to faster stimuli.

Both “on” and “off” MRCs increase in size with stimulus speed (Fig. 6, A–C, left; and Fig. 7, A and B, left), but to different extents. At low speeds, off MRCs have a larger peak current than on MRCs, but the ratio reverses at high speeds (Fig. 6 D). This observation was not recapitulated in simulations, suggesting that its origin lies in an element that is absent from the present computational model or an assumption that may not hold. For instance, the model assumes that tissue indentation follows the stimulator speed precisely. It is possible, however, that the stimulator bead detaches from the body during fast off stimuli, creating asymmetries in the mechanical loads delivered at high speed. Above, we argued that synchrony or phase differences in channel activation do not fully account for speed dependence generally. When we explored the possibility that phase differences contributed to the on/off asymmetry observed in peak current amplitudes by plotting the ratio of off to on charge against stimulus speed, we found that off responses consistently carried less charge than on responses. This ratio was not obviously speed dependent (Fig. 6 E).

On currents activate slowly at slow speeds, while off currents activate rapidly and vary less with stimulus speed in both experimental and simulated data (Fig. 7 C). The simulation predicts that recruiting more channels by decreasing $g0$ would increase the asymmetry in activation rate (Fig. 7 C, right). This is not true for the decay rate, however. Both on and off MRCs seem to decay at a relatively constant rate regardless of stimulus speed (Fig. 7 D). This suggests that activation depends on stimulus speed, while decay depends on tissue relaxation. Intuitively, we might expect channels to be more likely to be opened in-phase during off stimuli because all channels affected by the deformation will start moving once the stimulator moves. In contrast, during slow on stimuli, channels closer to the stimulator will be affected before the deformation reaches more distant channels, resulting in out-of-phase opening and slower MRC activation. This is supported by the observation that activation of on MRCs appears to be faster at high speeds (Figs. 6 B and 7 C). The activation rate of off currents is much less dependent on stimulus speed (Figs. 6 C and 7 C), but the amplitude of off currents still increases with speed despite this greater synchrony (Fig. 7 A), reinforcing the idea of another mechanism for speed dependence.

### Sinusoidal stimuli evoke steady MRCs in a frequency-dependent manner

The stimuli a worm encounters in the wild or in the laboratory (Nekimken et al., 2017b) involve complex signals with a variety of frequencies. Prior work shows that the TRNs seem poised to respond to variable signals: calcium responses to sinusoidal “buzz” stimuli are much stronger than those evoked by simple steps (Suzuki et al., 2003; Nekimken et al., 2017a; Fehlauer et al., 2018). Therefore, we systematically examined the frequency response to sinusoidal mechanical stimuli at frequencies ranging from 10 to 500 Hz. The representative traces in Fig. 8 A show both the signal from the photodiode measuring bead motion and the currents evoked by this stimulation. We see little evidence of adaptation following the application of a buzz-like stimulus: MRCs evoked by a step stimulus before and after a 1-s sinusoid show no more than 10–15% attenuation. This result starkly contrasts with studies on Piezo1 and Piezo2 channels, which adapt following sinusoidal stimulation (Lewis et al., 2017).

In agreement with our previously reported experimental and simulation results (Eastwood et al., 2015; Sanzeni et al., 2019), the expanded traces show that inward currents occur at twice the frequency of the stimulus (Fig. 8 B), especially at frequencies 500 Hz. This suggests an alternation between on and off responses. We next examined the power spectral density of the stimulus and the MRC response. We found that MRCs have roughly equal power at 1× and 2× the stimulus frequency (Fig. 8 C), whereas the spectrum of the stimulator shows a larger peak at the stimulus frequency and a smaller one at its second harmonic. The relative power of the second harmonic peak in the MRC response decreases slightly at higher frequencies, which is reflected in the decreasing asymmetry between subsequent (on vs. off) peaks (Fig. 8 B). This finding demonstrates that the neuron’s response to stimulation is nonlinear.

The average steady-state current continues to increase with frequency (Fig. 8, A and D). Although the experimental data are shifted to the right relative to the predictions of the model with the previously set parameters (τ = 2.9 ms, dashed line), both approaches show the tissue acting as a high-pass filter (Fig. 8 D). This is consistent with TRNs being relatively insensitive to slowly varying stimuli and thus ignoring self-motion. By contrast, the variance in the current, as measured by the RMS of the response at steady state, shows bandpass behavior (Fig. 8 E). Again, the peak of the experimental curve is shifted to higher frequencies than the model predicts with the previous parameters (dashed line). The bandpass behavior in Fig. 8 E can be understood as being due to the fact that channels do not have time to transition to lower-conductance states during high-frequency stimulation. Because the channels are held in an open or subconducting state, their contribution to the RMS noise declines with frequency.

We hypothesized that the simulated results were shifted to the left because channels were not able to close in phase with the sinusoidal stimulus. To test this possibility, we shortened the time required for the elastic filament to relax to baseline (τ = 0.6), which allowed us to better match the peak RMS (Fig. 8 E, dotted line). However, this change also shifted the steady-state current to the right. Collectively, these observations suggest that additional modifications to model parameters will be needed to fine-tune the mechanical filter that links touch sensation to MeT channel activation.

### Single-channel activation is a function of position and stimulus intensity

Our model allows us to explore spatial and temporal dynamics not only at the level of macroscopic currents, but also at the individual channel level. We selected four channels at varying distances from one realization of the simulation and calculated the expected on and off responses to various step displacements, variable speed trapezoidal profiles, and sine frequencies. Our goal was to examine how MRC dynamics vary along the length of the neurite and as a function of stimulus type.

The simulated channel current reflects the weighted probability that a channel at this location is in an open or subconducting state. The maximal current possible is −1.6 pA, which is the value calculated from the measured single-channel conductance (Brown et al., 2007) and our recording conditions. Channels that are more distant from the center of the stimulator respond only to large displacements (Fig. 9 A, highlighted). Channels directly under the stimulator (Fig. 9 A, rightmost column) reach a fully open state even with small displacements; with larger displacements, the open probability remains high for longer periods of time even when the step takes the same amount of time. Because the model parameters were fitted based on a previous dataset, the total current predicted by the simulation (Fig. 9 A, left panel) is smaller than what we observed at high speeds. This discrepancy could be due to uncertainty in some assumptions used in the model. For instance, channels in the model were distributed at regular 1.4-µm spacing, which is shorter than the average 2.95-µm spacing we observed in the region we were stimulating (Fig. 2). However, increasing the IPI would lead to smaller simulated currents. Another possibility is that individual puncta contain more than a single functional channel. Additional investigations of the response to small (≤1-µm) displacements with small stimulators, superresolution studies of channel protein position, or both are needed to gain further insight into this question.

When we held displacement constant and varied the speed, we found that on and off currents activate faster at higher speeds through a combination of synchrony and increasing open probability. Channels distant from the stimulus take longer to receive and respond to stimulus onset at slow speeds (Fig. 9 B, highlighted). Channels at the point of stimulation (Fig. 9 B, rightmost column) behave in the opposite manner: although the off current initially rises like the on current, it then slows and takes longer to reach the peak. Intermediate channels (Fig. 9 B, second to right column) open both faster and more fully with increased speed. This dependence on channel distance is inherited from the underlying strain dynamics and explains the predicted slow offset activation times (Fig. 7 C), which depend more strongly on the most proximal channels, and reflects the fact that the time course of strain rate depends on distance from the stimulus.

The model nicely captures the increase in steady-state current that occurs when we increase the frequency of sinusoidal stimuli. However, the fluctuations in the current (measured earlier as the RMS) are larger experimentally (Fig. 9 C, left panels). In both the observed and simulated data, the relative amplitudes of the peaks at 2× the stimulus frequency also decrease with frequency. Thus, the off response seems weaker. The open probability for channels near the stimulus remains more constant at higher frequencies (Fig. 9 C, highlighted column). The intensity of the stimulus reaching the channel is attenuated in distal channels, as seen by their lack of response to smaller or slower indentations (Fig. 9, A and B), which contributes to the decrease in RMS.

### Single-channel sensitivity also determines how many channels are opened

Having established the relationship between stimulus profile and the activation of channels proximal and distal to the stimulator, we next sought to derive an estimate of the energy required to activate single channels. Following Elmi et al. (2017), we calculated the SED, a measure of the energy required to deform the tissue at a given channel location; this quantity has been used as a proxy for the stimulus experienced by individual mechanoreceptors in various models of somatosensation (Grigg and Hoffman, 1984; Dandekar et al., 2003; Lesniak et al., 2014; Elmi et al., 2017). Values for SED were between 0.87 and 7.69 kPa, depending on the value for the IPI used in the model (Table 1). Our estimates of SED were two to three orders of magnitude larger than those reported by Elmi et al. (2017). There are two main factors contributing to this difference. First, SED is proportional to the Young’s modulus, which is 40 times larger in our simulations. Although the origin of this discrepancy is not clear, it may arise from differences in the finite element models used to infer the Young’s modulus. Second, different assumptions were used to deduce the minimum SED for single-channel activation. Rather than assuming that maximal currents are produced by full gating of all the available channels (Elmi et al., 2017), we accounted for activation of both nearby and distant channels, many of which remained in a subconducting state (Fig. S1).

Because the SED estimate applies only to the responses to displacement steps and involves several uncertainties, we turned to our model linking body indentation to activation of single channels via elongation of a hypothetical gating filament. In our model (Sanzeni et al., 2019), the filament extension needed to activate a single channel is of the order $g0/g1$; such an elongation is generated by a deformation of constant rate $g0/g1/τ$ if it is applied for a time of order, $τ$. Using values reported in Table 1, the model predicts that a single channel is able to sense a deformation on the order of 0.003 times the length of the linker generated over a time of 3 ms.

## Discussion

Using experiments and simulations, we advance understanding of how touch is transformed into depolarizing currents in TRNs. The picture emerging from this work is that touch sensitivity and filtering depend on a combination of factors: tissue mechanics and the distribution of ion channels within sensory neurons. These factors come together to generate mechanical forces transferred to individual channels whose activity accounts for the spatiotemporal dynamics of MRCs.

Additionally, our results show that the depth and speed of indentation affect the amplitude and time course of MRCs. Current–displacement curves do not saturate and instead display inflections that vary with internal pressure. By contrast, when we hold displacement constant, both the peak MRC–speed and charge–speed curves saturate at high speeds. The increase in peak current with speed occurs partially due to increased synchrony of channel openings with faster stimuli. Yet the total charge, which reflects all channel openings during the MRC, also increases with stimulus speed, suggesting that even without the effect of synchrony, MRCs depend on stimulus speed. This idea is also supported by the fact that off responses, which have more synchronous activation, show a similar dependence on stimulus speed. Collectively, these observations indicate that activating MRCs in their native environment depends on tissue mechanics and dynamics as well as on the intrinsic mechanosensitivity of the MeT channels that carry the currents.

### Two mechanisms govern MRC amplitude by defining the strain field generated during touch

Displacement (indentation) and speed affect MRCs via two mechanisms: (a) recruiting individual channels and (b) modulating the open probability for each channel. In general, the macroscopic current ($I$) is given by the product of the current carried by each individual channel ($i$), the number of channels available ($N$) and their probability of being open ($Po$). In many systems, we justifiably assume that $N$ is independent of stimulus intensity. We can also assume that while Po varies with stimulus intensity, it remains consistent across channels. For touch sensation, neither of these assumptions hold true. Our experimental and simulated data are consistent with the idea that larger displacements recruit increasingly distant channels by expanding the strain field. We previously showed that a displacement of a given magnitude elicits a larger response when the displacement begins from a preindented position (Sanzeni et al., 2019); this may occur due to a combination of wider recruitment and increased synchrony of activation. When we positioned the stimulus posterior to all channels, the current–displacement curve was shifted to the right (larger displacements) because we began to observe responses only once the displacement was large enough for the mechanical strain to reach the nearest channels (Figs. 4, 5, and 9). Our model also shows that both larger displacements and faster speeds increase the number of channels that fall within the indentation-induced strain field (Figs. 3 D and 7 F). Increasing either displacement or speed increases the probability that a channel will be in a subconducting or open state in a manner that depends on the channel’s distance from the stimulus (Fig. 9). The macroscopic current must take into account the different stimulus experienced by each channel as a result of tissue mechanics.

Further evidence that the joint effects of distal channel recruitment and modulation of open probability govern MRC amplitude comes from the finding that neither experimental nor simulated current–speed curves were well fitted by the saturating Boltzmann function expected if open probability were the dominant factor affecting current size. The total charge, in fact, appeared to have two or more distinct steps that may reflect differences in the relative contributions of increased open probability and recruitment (Fig. 7 B). Although we held the final displacement constant across speeds, the rate of strain, which determines the magnitude of the force acting on the channels in the model, increased with faster stimulation. Even if $N$ remains constant across all speeds, the distribution of open probabilities likely changes with curvature and strain near the stimulator. Curvature and strain, in turn, depend on body mechanics. Our model predicts that channels distant from the stimulator position are more likely to be in a subconducting rather than a fully open state. Moreover, increasing stimulator speed might increase the proportion of fully open versus subconducting channels in a nonlinear manner. Collectively, these two mechanisms work together to expand the TRN operating range, enabling this neuron to detect stimuli that are small and fast as well as those that are slow and large.

How might nematodes leverage these two mechanisms of TRN activation? One possibility is to escape from predatory fungi that ensnare nematodes by means of constricting traps. These traps inflate within ∼100 ms of contact (Higgins and Pramer, 1967), constricting to produce an initial indentation of ∼15% of the worm’s body width (Barron, 1979; Nordbring-Hertz et al., 2011). This initial rapid, but modest, stimulation may be followed by a slower but more complete closure of the constricting ring (Barron, 1979). Thus, the initial indentation would be expected to activate the TRNs swiftly, while a larger indentation could recruit distal channels to ensure a robust escape if the initial response failed. Consistent with this idea, C. elegans that lack functional TRNs are more likely to be trapped by such fungi in the laboratory (Maguire et al., 2011).

### MRC kinetics and on/off asymmetry depend on channel recruitment

The model supports the idea that channels located farther away from the center of stimulation are recruited and opened later. However, while the rate at which the strain changes is as important for channel activation as the magnitude of the strain, the difference in timing is noticeable only at slower speeds (Fig. 9 B). As seen in hair cell stereocilia (Ó Maoiléidigh and Ricci, 2019), distant channels experience increases in strain with a delayed time course because displacement along the worm is not instantaneous. The rate at which strain changes for distal channels can be high only when the displacement is large enough to reach them, which occurs late in an on ramp or early in an off ramp. Proximal channels, conversely, experience the greatest changes in strain at small displacements, where the curvature changes most drastically. As the proportion of distal to proximal channels increases with greater recruitment, the overall time course shifts toward the distal temporal profile, with slow on activation and fast off activation, creating asymmetry in MRC kinetics.

Similar to intact mammalian Pacinian corpuscles (Mendelson and Lowenstein, 1964), the DEG/ENaC/ASIC channels in C. elegans TRNs respond with depolarizing currents at both the onset and offset of an indentation step (O’Hagan et al., 2005; Arnadóttir et al., 2011; Eastwood et al., 2015). In this study, we also observed nearly identical responses at stimulus speeds of ∼1 mm/s, but varying the speed varied the asymmetry between peak on and off currents. The significant increase in activation time and the relative invariance of total charge across speeds suggest that the high off/on ratio at slow speeds is due to a change in synchrony (Fig. 7). Many MeT channels show asymmetrical responses to the onset and offset of a stimulus, due to inactivation (Lewis et al., 2017) or because the channels are part of asymmetrical structures that provide directional mechanical stimuli (Katta et al., 2015). Our current and previous models assume that individual channels receive symmetrical forces at the onset and offset of stimuli (Eastwood et al., 2015; Sanzeni et al., 2019). Here, we show that asymmetry of the stimulus is dynamic and is produced at the population level by the spatiotemporal recruitment of channels at different locations.

### Implications for understanding touch in vivo

The sensory neurons innervating Pacinian corpuscles show rapidly adapting, phasic responses similar to those of C. elegans mechanoreceptor neurons (Geffeney and Goodman, 2012). Adaptation of Pacinian corpuscle afferents slows following delamination of the corpuscle, indicating that the high-pass filter depends on the mechanics of surrounding tissues (Loewenstein and Mendelson, 1965; Loewenstein and Skalak, 1966). In addition, Pacinian corpuscles seem to require the lamellae to sustain the off response. Although C. elegans TRNs differ in their anatomic structure, the mechanics of the body as a fluid-filled shell creates a similar mechanical filter that leads to adaptation. Furthermore, the congruence between model dynamics and experimental dynamics lends credence to the model. We can pursue this line of research by directly measuring the time course and distribution of strain throughout the body and along the TRN neurite. How would these MeT channels respond if we stimulated them in cultured neurons without body mechanics filtering the force? Would the channels still show the same rapidly adapting behavior and responses to both stimulus onset and offset? The results presented here and in Sanzeni et al. (2019) indicate that tissue mechanics play a crucial role in shaping the response to touch. Thus, we would not expect to recapitulate the dependence on indentation depth and speed in vitro that is found in vivo.

Although the model reflects and explains many of the dynamics we see experimentally without being tuned for this dataset, there are still differences that hint at biological mechanisms as yet unaccounted for by the model. While many of the quantitative discrepancies can be resolved by better matching existing parameters to the recordings, a few qualitative observations may require more consideration. For large displacements, the total charge for on responses is consistently higher than the off response at all speeds, which might suggest asymmetry in how force is transferred from the bead to the tissue during indentation versus stimulus removal, inactivation at the level of the channel, or an irreversible dismantling of the molecules involved in transmitting force with repeated indentation. The presence of responses to distal posterior stimulation might be accounted for with simulation parameters tuned to this dataset, but it might also reflect our new observations about channel distribution. Perhaps the greater concentration of channels near the cell body allows the neuron to respond even when the strain field falls only partially on the neuron. Alternatively, extracellular collagen networks or cytoskeletal structures such as microtubules could render the neuron stiffer than the surrounding tissue and direct force down the neurite. Luckily, C. elegans also provides many genetic tools for perturbing these elements and assessing their contributions to the spatial and temporal dynamics of mechanoreceptor responses. We can now explore the effects of mutations in skin and neuronal proteins on mechano-electrical transduction within the context of a larger, overarching biophysical model.

## Acknowledgments

Richard W. Aldrich served as editor.

We thank Z. Liao for superb technical support and animal husbandry; V. Vásquez, A. Peng, F. Loizeau and A.R. Ricci for assistance with the design, fabrication, and testing of the mechanical stimulator and optical displacement sensor; and D. Ó’Maoiléidigh and members of the Goodman laboratory for comments on the manuscript.

Research supported by a Ruth L. Kirschstein fellowship (F31-NS-093825) to S. Katta and National Institutes of Health grants (R01-NS-047715 and R35-NS-105092) to M.B. Goodman.

The authors declare no competing financial interests.

All authors contributed to conceptualization, methodology, and revision of the final manuscript. Collection and analysis of electrophysiological data and indentation imaging data were performed by S. Katta, while A. Das was responsible for analyzing the position of MEC-4 puncta, both with supervision by M.B. Goodman. Simulations were carried out by A. Sanzeni, with supervision by M. Vergassola. S. Katta created the initial draft and visualizations.

## References

References
,
J.
,
R.
O’Hagan
,
Y.
Chen
,
M.B.
Goodman
, and
M.
Chalfie
.
2011
.
The DEG/ENaC protein MEC-10 regulates the transduction channel complex in Caenorhabditis elegans touch receptor neurons
.
J. Neurosci.
31
:
12695
12704
.
Audoly
,
B.
, and
Y.
Pomeau
.
2010
.
Elasticity and Geometry: From hair curls to the nonlinear response of shells.
Oxford University Press
,
Oxford, New York
.
Barron
,
G.L.
1979
.
Observations on predatory fungi
.
Can. J. Bot.
57
:
187
193
.
Bekkers
,
J.M.
, and
C.F.
Stevens
.
1996
.
Cable properties of cultured hippocampal neurons determined from sucrose-evoked miniature EPSCs
.
J. Neurophysiol.
75
:
1250
1255
.
Bounoutas
,
A.
,
R.
O’Hagan
, and
M.
Chalfie
.
2009
.
The multipurpose 15-protofilament microtubules in C. elegans have specific roles in mechanosensation
.
Curr. Biol.
19
:
1362
1367
.
Brown
,
A.L.
,
S.M.
Fernandez-Illescas
,
Z.
Liao
, and
M.B.
Goodman
.
2007
.
Gain-of-function mutations in the MEC-4 DEG/ENaC sensory mechanotransduction channel alter gating and drug blockade
.
J. Gen. Physiol.
129
:
161
173
.
Brown
,
A.L.
,
Z.
Liao
, and
M.B.
Goodman
.
2008
.
MEC-2 and MEC-6 in the Caenorhabditis elegans sensory mechanotransduction complex: auxiliary subunits that enable channel activity
.
J. Gen. Physiol.
131
:
605
616
.
Chalfie
,
M.
, and
J.
Sulston
.
1981
.
Developmental genetics of the mechanosensory neurons of Caenorhabditis elegans
.
Dev. Biol.
82
:
358
370
.
Chalfie
,
M.
, and
J.N.
Thomson
.
1979
.
Organization of neuronal microtubules in the nematode Caenorhabditis elegans
.
J. Cell Biol.
82
:
278
289
.
Chelur
,
D.S.
,
G.G.
Ernstrom
,
M.B.
Goodman
,
C.A.
Yao
,
L.
Chen
,
R.
O’ Hagan
, and
M.
Chalfie
.
2002
.
The mechanosensory protein MEC-6 is a subunit of the C. elegans touch-cell degenerin channel
.
Nature.
420
:
669
673
.
Chen
,
X.
, and
M.
Chalfie
.
2015
.
Regulation of mechanosensation in C. elegans through ubiquitination of the MEC-4 mechanotransduction channel
.
J. Neurosci.
35
:
2200
2212
.
Cueva
,
J.G.
,
A.
Mulholland
, and
M.B.
Goodman
.
2007
.
Nanoscale organization of the MEC-4 DEG/ENaC sensory mechanotransduction channel in Caenorhabditis elegans touch receptor neurons
.
J. Neurosci.
27
:
14089
14098
.
Dandekar
,
K.
,
B.I.
Raju
, and
M.A.
Srinivasan
.
2003
.
3-D finite-element models of human and monkey fingertips to investigate the mechanics of tactile sense
.
J. Biomech. Eng.
125
:
682
691
.
Eastwood
,
A.L.
,
A.
Sanzeni
,
B.C.
Petzold
,
S.-J.
Park
,
M.
Vergassola
,
B.L.
Pruitt
, and
M.B.
Goodman
.
2015
.
Tissue mechanics govern the rapidly adapting and symmetrical response to touch
.
112
:
E6955
E6963
.
Edelstein
,
A.
,
Amodaj
,
N.
,
Hoover
,
K.
,
Vale
,
R.
,
Stuurman
,
N.
,
2010
. Computer Control of Microscopes Using µManager. Curr. Protoc. Mol. Biol. 92:14.20.1–14.20.17.
Elmi
,
M.
,
V.M.
Pawar
,
M.
Shaw
,
D.
Wong
,
H.
Zhan
, and
M.A.
Srinivasan
.
2017
.
Determining the biomechanics of touch sensation in C. elegans
.
Sci. Rep.
7
:
12329
.
Emtage
,
L.
,
G.
Gu
,
E.
Hartwieg
, and
M.
Chalfie
.
2004
.
Extracellular proteins organize the mechanosensory channel complex in C. elegans touch receptor neurons
.
Neuron.
44
:
795
807
.
Fehlauer
,
H.
,
A.L.
Nekimken
,
A.A.
Kim
,
B.L.
Pruitt
,
M.B.
Goodman
, and
M.
Krieg
.
2018
.
Using a Microfluidics Device for Mechanical Stimulation and High Resolution Imaging of C. elegans
.
J. Vis. Exp.
(
132
):
56530
.
Geffeney
,
S.L.
, and
M.B.
Goodman
.
2012
.
How we feel: ion channel partnerships that detect mechanical inputs and give rise to touch and pain perception
.
Neuron.
74
:
609
619
.
Geffeney
,
S.L.
,
J.G.
Cueva
,
D.A.
Glauser
,
J.C.
Doll
,
T.H.-C.
Lee
,
M.
Montoya
,
S.
Karania
,
A.M.
Garakani
,
B.L.
Pruitt
, and
M.B.
Goodman
.
2011
.
DEG/ENaC but not TRP channels are the major mechanoelectrical transduction channels in a C. elegans nociceptor
.
Neuron.
71
:
845
857
.
Goodman
,
M.B.
,
D.H.
Hall
,
L.
Avery
, and
S.R.
Lockery
.
1998
.
Active currents regulate sensitivity and dynamic range in C. elegans neurons
.
Neuron.
20
:
763
772
.
Grigg
,
P.
, and
A.H.
Hoffman
.
1984
.
Ruffini mechanoreceptors in isolated joint capsule: responses correlated with strain energy density
.
Somatosens. Res.
2
:
149
162
.
Higgins
,
M.L.
, and
D.
Pramer
.
1967
.
Fungal Morphogenesis: Ring Formation and Closure by Arthrobotrys dactyloides
.
Science.
155
:
345
346
.
Hostettler
,
L.
,
Grundy
,
L.
,
Käser-Pébernard
,
S.
,
Wicky
,
C.
,
Schafer
,
W.R.
,
Glauser
,
D.A.
,
2017
. The Bright Fluorescent Protein mNeonGreen Facilitates Protein Expression Analysis In Vivo. G3 Genes Genomes Genet. 7:607–615.
Kang
,
L.
,
J.
Gao
,
W.R.
Schafer
,
Z.
Xie
, and
X.Z.S.
Xu
.
2010
.
C. elegans TRP family protein TRP-4 is a pore-forming subunit of a native mechanotransduction channel
.
Neuron.
67
:
381
391
.
Katta
,
S.
,
M.
Krieg
, and
M.B.
Goodman
.
2015
.
Feeling force: physical and physiological principles enabling sensory mechanotransduction
.
Annu. Rev. Cell Dev. Biol.
31
:
347
371
.
Krieg
,
M.
,
J.
Stühmer
,
J.G.
Cueva
,
R.
Fetter
,
K.
Spilker
,
D.
Cremers
,
K.
Shen
,
A.R.
Dunn
, and
M.B.
Goodman
.
2017
.
Genetic defects in β-spectrin and tau sensitize C. elegans axons to movement-induced damage via torque-tension coupling
.
eLife.
6
:e20172.
Lesniak
,
D.R.
, and
G.J.
Gerling
.
2009
.
Predicting SA-I mechanoreceptor spike times with a skin-neuron model
.
Math. Biosci.
220
:
15
23
.
Lesniak
,
D.R.
,
K.L.
Marshall
,
S.A.
Wellnitz
,
B.A.
Jenkins
,
Y.
Baba
,
M.N.
Rasband
,
G.J.
Gerling
, and
E.A.
Lumpkin
.
2014
.
Computation identifies structural features that govern neuronal firing properties in slowly adapting touch receptors
.
eLife.
3
:e01488.
Lewis
,
A.H.
,
A.F.
Cui
,
M.F.
McDonald
, and
J.
Grandl
.
2017
.
Transduction of Repetitive Mechanical Stimuli by Piezo1 and Piezo2 Ion Channels
.
Cell Reports.
19
:
2572
2585
.
Loewenstein
,
W.R.
, and
M.
Mendelson
.
1965
.
Components of receptor adaptation in a Pacinian corpuscle
.
J. Physiol.
177
:
377
397
.
Loewenstein
,
W.R.
, and
R.
Skalak
.
1966
.
Mechanical transmission in a Pacinian corpuscle. An analysis and a theory
.
J. Physiol.
182
:
346
378
.
Maguire
,
S.M.
,
C.M.
Clark
,
J.
Nunnari
,
J.K.
Pirri
, and
M.J.
Alkema
.
2011
.
The C. elegans touch response facilitates escape from predacious fungi
.
Curr. Biol.
21
:
1326
1330
.
Mendelson
,
M.
, and
W.R.
Lowenstein
.
1964
.
.
Science.
144
:
554
555
.
Nekimken
,
A.L.
,
H.
Fehlauer
,
A.A.
Kim
,
S.N.
Manosalvas-Kjono
,
P.
,
F.
Memon
,
D.
Gopisetty
,
V.
Sanchez
,
M.B.
Goodman
,
B.L.
Pruitt
, and
M.
Krieg
.
2017
a
.
Pneumatic stimulation of C. elegans mechanoreceptor neurons in a microfluidic trap
.
Lab Chip.
17
:
1116
1127
.
Nekimken
,
A.L.
,
E.A.
Mazzochette
,
M.B.
Goodman
, and
B.L.
Pruitt
.
2017
b
.
Forces applied during classical touch assays for Caenorhabditis elegans
.
PLoS One.
12
:e0178080.
Nordbring-Hertz
,
B.
,
H.-B.
Jansson
, and
A.
Tunlid
.
2011
. Nematophagous Fungi. In
eLS.
John Wiley Sons
,
Chichester, UK
.
O’Hagan
,
R.
,
M.
Chalfie
, and
M.B.
Goodman
.
2005
.
The MEC-4 DEG/ENaC channel of Caenorhabditis elegans touch receptor neurons transduces mechanical signals
.
Nat. Neurosci.
8
:
43
50
.
Ó Maoiléidigh
,
D.
, and
A.J.
Ricci
.
2019
.
A Bundle of Mechanisms: Inner-Ear Hair-Cell Mechanotransduction
.
Trends Neurosci.
42
:
221
236
.
Peng
,
A.W.
, and
A.J.
Ricci
.
2016
.
Glass Probe Stimulation of Hair Cell Stereocilia
. In
B.
Sokolowski
, editor.
Auditory and Vestibular Research: Methods and Protocols.
Springer
,
New York
;
487
500
.
Petzold
,
B.C.
,
S.-J.
Park
,
E.A.
Mazzochette
,
M.B.
Goodman
, and
B.L.
Pruitt
.
2013
.
MEMS-based force-clamp analysis of the role of body stiffness in C. elegans touch sensation
.
Integr. Biol.
5
:
853
864
.
Quindlen
,
J.C.
,
V.K.
Lai
, and
V.H.
Barocas
.
2015
.
Multiscale Mechanical Model of the Pacinian Corpuscle Shows Depth and Anisotropy Contribute to the Receptor’s Characteristic Response to Indentation
.
PLOS Comput. Biol.
11
:e1004370.
Quindlen-Hotek
,
J.C.
, and
V.H.
Barocas
.
2018
.
A finite-element model of mechanosensation by a Pacinian corpuscle cluster in human skin
.
Biomech. Model. Mechanobiol.
17
:
1053
1067
.
Rall
,
W.
1967
.
Distinguishing theoretical synaptic potentials computed for different soma-dendritic distributions of synaptic input
.
J. Neurophysiol.
30
:
1138
1168
.
Rivera Gomez
,
K.A.
, and
M.
Schvarzstein
.
2018
.
Immobilization of nematodes for live imaging using an agarose pad produced with a vinyl record.
Sanzeni
,
A.
,
S.
Katta
,
B.
Petzold
,
B.L.
Pruitt
,
M.B.
Goodman
, and
M.
Vergassola
.
2019
.
Somatosensory neurons integrate the geometry of skin deformation and mechanotransduction channels to shape touch sensing
.
eLife.
8
:e43226.
Suzuki
,
H.
,
R.
Kerr
,
L.
Bianchi
,
C.
Frøkjaer-Jensen
,
D.
Slone
,
J.
Xue
,
B.
Gerstbrein
,
M.
Driscoll
, and
W.R.
Schafer
.
2003
.
In vivo imaging of C. elegans mechanosensory neurons demonstrates a specific role for the MEC-4 channel in the process of gentle touch sensation
.
Neuron.
39
:
1005
1017
.
Thein
,
M.C.
,
G.
McCormack
,
A.D.
Winter
,
I.L.
Johnstone
,
C.B.
Shoemaker
, and
A.P.
Page
.
2003
.
Caenorhabditis elegans exoskeleton collagen COL-19: an adult-specific marker for collagen modification and assembly, and the analysis of organismal morphology
.
Dev. Dyn.
226
:
523
539
.
Vásquez
,
V.
,
M.
Krieg
,
D.
, and
M.B.
Goodman
.
2014
.
Phospholipids that contain polyunsaturated fatty acids enhance neuronal cell mechanics and touch sensation
.
Cell Reports.
6
:
70
80
.

## Author notes

A. Sanzeni’s present address is Dept. of Neurobiology, Duke University, Durham, NC.

S. Katta’s present address is American Society of Human Genetics, Rockville, MD.