Scanned optogenetic control of mammalian somatosensory input to map input-specific behavioral outputs
Abstract
Somatosensory stimuli guide and shape behavior, from immediate protective reflexes to longer-term learning and higher-order processes related to pain and touch. However, somatosensory inputs are challenging to control in awake mammals due to the diversity and nature of contact stimuli. Application of cutaneous stimuli is currently limited to relatively imprecise methods as well as subjective behavioral measures. The strategy we present here overcomes these difficulties, achieving ‘remote touch’ with spatiotemporally precise and dynamic optogenetic stimulation by projecting light to a small defined area of skin. We mapped behavioral responses in freely behaving mice with specific nociceptor and low-threshold mechanoreceptor inputs. In nociceptors, sparse recruitment of single-action potentials shapes rapid protective pain-related behaviors, including coordinated head orientation and body repositioning that depend on the initial body pose. In contrast, activation of low-threshold mechanoreceptors elicited slow-onset behaviors and more subtle whole-body behaviors. The strategy can be used to define specific behavioral repertoires, examine the timing and nature of reflexes, and dissect sensory, motor, cognitive, and motivational processes guiding behavior.
eLife digest
To safely navigate their world, animals need to be able to tell apart a gentle touch from an eye-watering pinch, detect cold water or sense the throbbing pain stemming from an infected cut. These ‘somatic’ sensations are relayed through thousands of nerve endings embedded in the skin and other tissues. Yet the neurological mechanisms that underpin these abilities are complex and still poorly understood.
Indeed, these nerve endings can be stimulated by extreme temperatures, harmful chemicals, friction or even internal signals such as inflammation. One event can also recruit many different types of endings: a cut for example, will involve responses to mechanical pressure, tissue damage and local immune response. To disentangle these different actors and how they affect behavior, scientists need to develop approaches that allow them to deliver specific stimuli with increased precision, and to monitor the impact on an animal.
To achieve this goal, Schorscher-Petcu et al. used mice in which blue light could trigger specific types of nerve endings. For instance, depending on the genetic background of the animals, a laser could either activate nerve endings involved in pain or gentle touch. Crucially, this could be done from a distance by beaming light with exquisite precision onto the paws of the mice without physically touching or disturbing the animals.
How the mice responded could then be observed without any interference. Their behavior was analyzed using a combination of high-speed videos, computer-driven recording systems, and machine learning. This revealed subtle changes in behavior that had not been detected before, spotting microscopic movements of the stimulated paw and mapping simultaneous whole-body movements such as changes in posture or head orientation. The approach therefore allows scientists to assess the impact of touch, pain or temperature sensation in freely behaving mice. It could also be harnessed to develop much needed treatments against chronic pain.
Introduction
The survival of an organism depends on its ability to detect and respond appropriately to its environment. Afferent neurons innervating the skin provide sensory information to guide and refine behavior (Seymour, 2019; Zimmerman et al., 2014). Cutaneous stimuli are used to study a wide range of neurobiological mechanisms since neurons densely innervating skin function to provide diverse information as the body interfaces with its immediate environment. These afferents maintain the integrity of the body by recruiting rapid sensorimotor responses, optimize movement through feedback loops, provide teaching signals that drive learning, and update internal models of the environment through higher-order perceptual and cognitive processes (Barik et al., 2018; Brecht, 2017; Corder et al., 2019; de Haan and Dijkerman, 2020; Haggard et al., 2013; Huang et al., 2019; Petersen, 2019; Seymour, 2019). Damaging stimuli, for example, evoke rapid motor responses to minimize immediate harm and generate pain that motivates longer-term behavioral changes.
Compared to visual, olfactory, and auditory stimuli, somatosensory inputs are challenging to deliver in awake unrestrained mammals. This is due to the nature of stimuli that require contact and the diversity of stimulus features encoded by afferents that innervate skin. Cutaneous afferent neurons are functionally and genetically heterogeneous, displaying differential tuning, spike thresholds, adaptation rates, and conduction velocities (Abraira and Ginty, 2013; Dubin and Patapoutian, 2010; Gatto et al., 2019; Häring et al., 2018). The arborization of their peripheral terminals can delineate spatial and temporal dimensions of the stimulus (Pruszynski and Johansson, 2014), particularly once many inputs are integrated by the central nervous system (Prescott et al., 2014). Cutaneous stimulation in freely moving mice often requires the experimenter to manually touch or approach the skin. This results in inaccurate timing, duration, and localization of stimuli. The close proximity of the experimenter can cause observer-induced changes in animal behavior (Sorge et al., 2014). Stimuli also activate a mixture of sensory neuron populations. For example, intense stimuli can co-activate fast-conducting low-threshold afferents that encode innocuous stimuli simultaneously with more slowly conducting high-threshold afferents (Wang et al., 2018). The latter are nociceptors that trigger fast protective behaviors and pain. Consequently, mixed cutaneous inputs recruit cells, circuits, and behaviors that are not specific to the neural mechanism under study. A way to control genetically defined afferent populations is to introduce opsins into these afferents and optogenetically stimulate them through the skin (Abdo et al., 2019; Arcourt et al., 2017; Barik et al., 2018; Beaudry et al., 2017; Browne et al., 2017; Daou et al., 2013; Iyer et al., 2014). However, these methods in their current form do not fully exploit the properties of light.
The behaviors that are evoked by cutaneous stimuli are also typically measured with limited and often subjective means. Manual scoring introduces unnecessary experimenter bias and omits key features of behavior. Behavioral assays have traditionally focused on a snapshot of the stimulated body part rather than dynamics of behavior involving the body as a whole (Gatto et al., 2019). Recent advances in machine vision and markerless pose estimation have enabled the dissection of animal behavioral sequences (Mathis et al., 2018; Pereira et al., 2019; Wiltschko et al., 2015). However, these have not been adapted to study behavioral outputs relating to specific cutaneous inputs.
Here we developed an approach to project precise optogenetic stimuli onto the skin of freely behaving mice (Figure 1A). The strategy elicits time-locked individual action potentials in genetically targeted afferents innervating a small stimulation field targeted to the skin. Stimuli can be delivered remotely as predefined microscale patterns, lines, or moving points. The utility of the system was demonstrated by precisely stimulating nociceptors, or Aβ low threshold mechanoreceptors (LTMRs), in freely behaving mice to map behavioral outputs at high speed. We provide an analysis toolkit that quantifies the millisecond-timescale dynamics of behavioral responses using machine vision methods. We dissect discrete behavioral components of local paw responses, head orienting and body repositioning behaviors, and determine how these specific behavioral components relate to precise somatosensory inputs.

Remote and precise somatosensory input and analysis of behavior.
(A) Afferent neurons expressing ChR2 are controlled remotely in freely behaving mice by projecting laser light with sub-millimeter precision to the skin. This enables precise non-contact stimulation with microscale patterns, lines, and points using scanned transdermal optogenetics. Time-locked triggering of single-action potential volleys is achieved through high temporal control of the laser. Behavioral responses can be automatically recorded and analyzed using a combination of computational methods. (B) Schematic of the stimulation laser (in blue) and infrared imaging (in red) paths. Mirrors (M1 and M2) direct the laser beam through a set of lenses (L1–L3), which allow the beam to be focused manually to pre-calibrated spot sizes. A dichroic mirror (DM) guides the laser beam into a pair of galvanometer mirrors, which are remotely controlled to enable precise targeting of the beam onto the glass platform. Near-infrared frustrated total internal reflection (NIR-FTIR) signal from the glass platform is descanned through the galvanometers and imaged using a high-speed infrared camera. A second wide-field camera is used to concomitantly record a below view of the entire glass platform. (C) Rendering of the assembled components. A Solidworks assembly is available at https://github.com/browne-lab/throwinglight.
Results
Design of the optical approach
The design of the optical strategy had eight criteria: (1) that somatosensory stimuli are delivered non-invasively without touching or approaching the mice; (2) localization of stimuli are spatially precise and accurate (<10 μm); (3) freely moving mice can be targeted anywhere within a relatively large (400 cm2) arena; (4) stimuli can be controlled with a computer interface from outside the behavior room; (5) stimulation patterns, lines, and points are generated by rapidly scanning the stimuli between predefined locations; (6) stimulation size can be controlled down to ≥150 μm diameter; (7) stimuli are temporally precise to control individual action potentials using sub-millisecond time-locked pulses; and (8) behavioral responses are recorded at high speed at the stimulated site and across the whole body simultaneously. An optical system was assembled to meet these specific criteria (Figure 1B and C).
The stimulation path uses two mirror galvanometers to remotely target the laser stimulation to any location on a large glass stimulation floor. A series of lenses expands the beam and then focuses it down to 0.018 mm2 (150 μm beam diameter) at the surface of this floor. This was defocused to provide a range of calibrated stimulation spot sizes up to 2.307 mm2, with separable increments that were stable over long periods of time (Figure 1—figure supplement 1A). The optical power density could be kept equal between these different stimulation spot sizes. The glass floor was far (400 mm) from the galvanometers, resulting in a maximum focal length variability of <1.5% (see Materials and methods). This design yielded a spatial targeting resolution of 6.2 μm while minimizing variability in laser stimulation spot sizes across the large stimulation plane (coefficient of variation ≤0.1, Figure 1—figure supplement 1B). The beam ellipticity was 74.3% ± 14.3% (median± MAD, range of 36–99%) for all spot sizes. The optical power was uniform across the stimulation plane (Figure 1—figure supplement 1C). The galvanometers allow rapid small angle step (300 µs) responses to scan the laser beam between adjacent positions and shape stimulation patterns using brief laser pulses (diode laser rise and fall time: 2.5 ns). Custom software (see Materials and methods) was developed to remotely control the laser stimulation position, trigger laser pulses, synchronize galvanometer jumps, and trigger the camera acquisition (Figure 1—figure supplement 2).
The camera acquisition path was used to manually target the location of the laser stimulation pulse(s); the path was descanned through the galvanometers so that the alignment between the laser and camera is fixed (Figure 1B). The camera feed is displayed in the user interface and enables the operator to use this image to target the laser to the desired location. High signal-to-noise recordings were obtained using near-infrared frustrated total internal reflection (NIR-FTIR) in the glass stimulation floor (Roberson, D. P. et al., manuscript submitted). If a medium (skin, hair, tail, etc.) is within a few hundred microns of the glass, it causes reflection of the evanescent wave and this signal decreases non-linearly with distance from the glass such that very minor movements of the paw can be detected. The acquisition camera acquired the NIR-FTIR signal in high-speed (up to 1000 frames/s) with a pixel size of 110 μm. A second camera was used to record the entire arena and capture behaviors involving the whole body before and after stimulation. Offline quantification was carried out using custom analysis code combined with markerless tracking tools (Mathis et al., 2018).
Mapping high-speed local responses to nociceptive input
To validate the strategy, we first crossed Trpv1-IRES-Cre (TRPV1Cre) and R26-CAG-LSL-ChR2-tdTomato mice to obtain a line (TRPV1Cre::ChR2) in which ChR2 is selectively expressed in a broad class of nociceptors innervating glabrous skin (Browne et al., 2017). These mice were allowed to freely explore individual chambers placed on the stimulation plane. When mice were idle (still and awake), a time-locked laser pulse was targeted to the hind paw. Stimuli could be controlled remotely from outside the behavior room. We recorded paw withdrawal dynamics with millisecond resolution. For example, a single, small 1 ms laser pulse initiated a behavioral response at 29 ms, progressing to complete removal of the hind paw from the glass floor just 5 ms later (Figure 2A and Figure 2—video 1). The stimulus used for this protocol was S6, 0.577 mm2 in area, which corresponds to less than 1% of the glabrous paw area and highlights the sensitivity of the nociceptive system. Motion energy, individual pixel latencies, and response dynamics could be extracted from these high-speed recordings (Figure 2B and C).

Scanned optogenetic stimuli reveal relationships with local behaviors.
(A) Millisecond-timescale changes in hind paw near-infrared frustrated total internal reflection (NIR-FTIR) signal in response to a single 1 ms laser pulse (laser spot size S6 = 0.577 mm2) recorded at 1000 frames/s. (B) Motion energy analysis (left) and response latencies calculated for each pixel (right) for the same trial as in (A). (C) Example traces of the NIR-FTIR signal time course as measured within a circular region of interest centered on the stimulation site. Six traces from two animals are depicted (1 ms pulse, spot size S6 = 0.577 mm2). The red trace corresponds to the example trial illustrated in (A) and (B). (D) Paw response probability increases as a function of laser pulse duration when stimulation size is constant (spot size S6 = 0.577 mm2; 37–42 trials for each pulse duration from eight mice, mean probability ± SEM). Light pulses 10 ms or less with the same intensity and wavelength have been shown to generate just a single-action potential in each nociceptor activated in the TRPV1Cre::ChR2 line (Browne et al., 2017). Note that a 30 ms might generate more than one action potential but the response already plateaus at 10 ms duration, suggesting one action potential per nociceptor shapes the response. (E) Paw response probability increases as a function of laser stimulation spot size when pulse duration is constant. Data are 34–45 trials for each spot size per pulse duration from 7 to 8 mice, shown as mean probability ± SEM. The dataset for (D) and (E) is provided in Figure 2—source data 1. (F) Stimulation patterning shows that the absolute size, rather than the geometric shape, of the nociceptive stimulus determines the withdrawal probability (Friedman’s non-parametric test for within subject repeated measures S(5) = 22.35, p=0.0004). Paw response probabilities in response to a single large laser spot (S7 = 1.15 mm2), a single small spot (S4 = 0.176 mm2; p=0.018 compared to S7 and p=0.013 compared to the line pattern), a 10 ms train of seven small 1 ms spots targeting the same site (p=0.039, compared to S7 and p=0.030 compared to the line pattern) or spatially translated to produce different patterns. Note that the cumulative area of the seven small spots approximates the area of the large spot, and no statistically significant difference was detected between any of their response probabilities. Data shown as mean probability ± SEM are from n = 6 mice, with each 6–10 trials per pattern. The dataset for (F) is provided in Figure 2—source data 2.
-
Figure 2—source data 1
Time courses of paw movement recorded at 1000 frames/s with stimuli that vary in duration and size.
Stimuli (40 mW/mm2) were delivered at 0 ms. Data are from TRPV1Cre::ChR2 mice and littermate controls.
- https://cdn.elifesciences.org/articles/62026/elife-62026-fig2-data1-v2.csv
-
Figure 2—source data 2
Time courses of paw movement recorded at 500 frames/s with single point and patterned stimuli.
Stimuli (1 ms, 40 mW/mm2) were delivered at 0 ms. Data are from TRPV1Cre::ChR2 mice and littermate controls.
- https://cdn.elifesciences.org/articles/62026/elife-62026-fig2-data2-v2.csv
We probed multiple sites across the plantar surface and digits and found that the hind paw heel gave the most robust responses (Figure 2—figure supplement 1). This region was targeted in all subsequent experiments. Littermates that did not express the Cre recombinase allele confirmed that the laser stimulation did not produce non-specific responses. These mice did not show any behavioral responses, even with the largest stimuli (spot size S8, 30 ms pulse, Figure 2—figure supplement 2). We next provide some examples of the utility of the strategy by examining the relationship between nociceptive input and protective behaviors.
Probabilistic nociceptor recruitment determines the nature, timing, and extent of behavior
Fast protective withdrawal behaviors can be triggered by the first action potential arriving at the spinal cord from cutaneous nociceptors. A brief optogenetic stimulus generates just a single-action potential in each nociceptor activated (Browne et al., 2017). This is due to the rapid closing rate of ChR2 relative to the longer minimal interspike interval of nociceptors. The same transient optogenetic stimulus (Browne et al., 2017), or a pinprick stimulus (Arcourt et al., 2017), initiates behavior before a second action potential would have time to arrive at the spinal cord. That the first action potential can drive protective behaviors places constraints on how stimulus intensity can be encoded, suggesting that the total population of nociceptors firing a single-action potential can provide information as a "Boolean array." The consequences of this have not been investigated previously as precise control of specific nociceptive input had not been possible. We predicted that the relative number of nociceptors firing a single-action potential determines the features of the behavioral response.
Varying the pulse duration with nanosecond precision influences the probability of each nociceptor generating a single-action potential within the stimulation site. A pulse as short as 300 μs elicited behavioral responses but with relatively low probability (Figure 2D). This probability increased with pulse duration until it approached unity, closely matching the on-kinetics of the ChR2 used (τ = 1.9 ms; Lin, 2011). We next controlled the spatial, rather than temporal, properties of the stimulation in two further experiments. Firstly, we find that the total area of stimulated skin determines the behavioral response probability, such that the larger the nociceptive input the larger the response probability (Figure 2E). Secondly, we generated different stimulation patterns. We find that sub-threshold stimulations are additive (Figure 2F). Specifically, seven spatially displaced small sub-threshold stimulations could reproduce the response probability of a single large stimulation that was approximately seven times their size. This could not be achieved by repeated application of the small stimulations to the same site (Figure 2F).
Time-locking the stimulus enabled us to examine the hind paw responses with high temporal resolution. The nociceptive input size influenced the behavioral response latency: for example, a 3 ms pulse resulted in response latencies of 27 ± 1 ms, 30 ± 2 ms, 33 ± 5 ms, and 112 ± 46 ms for spot sizes S8, S7, S6, and S5, respectively (Figure 3A and B). The shorter latencies are consistent with medium-conduction velocity Aδ-fibers that arrive at the spinal cord before slower C-fiber action potentials (>35 ms) (Browne et al., 2017). The rank order of response latencies follows the nociceptive input size for both pulse durations, and they fit well with log-log regressions (3 ms pulse R2 = 0.87, 1 ms pulse R2 = 0.90). Once a hind limb motor response was initiated, it developed rapidly, lifting from the glass with rise times that show the vigor of the motor response was also dependent on nociceptive input size (Figure 3C). These responses, in >65% of cases, proceeded to full withdrawal. However, in a fraction of trials the paw moved but did not withdraw (Figure 3D), highlighting the sensitivity of the acquisition system. Even the smallest of nociceptive inputs still produced a large fraction of full withdrawal responses, despite decreases in response probability (Figure 3E). The fraction of full withdrawal responses increased with the size of nociceptive input. The onset latency of both full and partial responses decreased as nociceptive input increased (Figure 3F).

Paw response latency and magnitude are influenced by the sparse recruitment of nociceptors.
(A) Raster plots of hind paw responses for five different 3 ms laser stimulation sizes, sorted by response latency. The paw response latency is indicated in red. (B) Paw response latencies to trials with single 3 ms (blue, top) and 1 ms (green, bottom) stimulations at different spot sizes, sorted by latency. (C) Response vigor (hind paw rise time, 20–80%) to single 3 ms (blue, top) or 1 ms (green, bottom) pulses with a range of stimulation spot sizes. Rise times to a 3 ms pulse were 4 ± 1 ms, 4 ± 1 ms, 4 ± 1 ms, and 9 ± 5 ms for spot sizes S8, S7, S6, and S5, respectively, and to a 1 ms pulse were 4 ± 1 ms, 5 ± 2 ms, and 6 ± 3 ms for spot sizes S8, S7, and S6, respectively. (D) Extent of responses (%NIR-FTIR signal decrease). The threshold for a full response and partial response is 75% of baseline signal (red line). (E) The probability of responses to reach completion (full response) as a function of the probability of response for four stimulation spot sizes and two pulse durations (green 1 ms; blue 3 ms). (F) Response latency distributions for trials that reach completion (full response) shown with Gaussian kernel density estimation of data (left). Rug plot inset representing individual response latencies for each color-coded spot size. No correlation was observed between response latency and extent for partial responses when stimulation duration was 3 ms. Data from 7 to 8 mice with 39–44 trials per spot size for 1 ms pulse duration and 34–44 per spot size for 3 ms pulse duration. The dataset is provided in Figure 2—source data 1. NIR-FTIR: near-infrared frustrated total internal reflection.
Whole-body behavioral responses to remote and precise nociceptive input
Pain-related responses are not limited to the affected limb but involve simultaneous movement of other parts of the body (Blivis et al., 2017; Browne et al., 2017). These non-local behaviors theoretically serve several protective purposes: to investigate and identify the potential source of danger, move the entire body away from this danger, attend to the affected area of the body (Huang et al., 2019) and to maintain balance (Sherrington, 1910). Whole-body movements were quantified as motion energy (Figure 4—figure supplement 1A) and high-speed recordings show this initiated with a mean response latency of 30 ± 1 ms, with the first movement bout displaying a mean duration of 136 ± 14 ms (80 trials from 10 mice) (Figure 4—figure supplement 2). The magnitude of whole-body movement increased with the stimulation spot size (Figure 4—figure supplement 1B). Peak motion energy had a lognormal relationship with nociceptive input size (R2 = 0.99). This indicates that global behaviors are also proportional to the relative size of the nociceptive input; the recruited nociceptors firing a single-action potential (Figure 4—figure supplement 1B).
Sparse nociceptor stimulation triggers coordinated postural adjustments
Most behaviors arise from the complex coordination of discrete body parts, which can be tracked individually. To dissect specific components of these behaviors, we implemented DeepLabCut (Mathis et al., 2018) by training a network using frames from the high-speed (400 frames/s) videos to track 18 user-defined body parts across the mouse (for details, refer to Materials and methods, Global behaviors during optogenetic stimulation). The high-speed video recordings of stimulation trials were analyzed using this network. Specific nociceptive input at the hind paw (S8, 2.307 mm2, 10 ms pulse) causes behavior that initiates simultaneously across the body. Inspection of the movements of each body part relative to the baseline pose (Figure 4A) shows fast outward movement of the stimulated and contralateral hind paws, and concomitant initiation of head orientation (two example responses in Figure 4B). Based on these observations, we examined the behavioral trajectories in the first 115 ms across the population of 80 trials. The first three principal components (PCs) were fit using six body part x and y values at 115 ms after the stimulus onset. These PCs explain 88.8% of the variance (50.4, 26.5, and 11.9% for PC1, PC2, and PC3, respectively). PC1 is dominated by hind paw translation, PC2 by head and body movement, and PC3 by head orientation (Figure 4C). Projecting the entire time course onto these same PCs can explain 78.1% of the variance (37.1, 24.3, and 16.7% for PC1, PC2, and PC3, respectively). The response trajectories revealed that movements occur largely in same direction within PC space with a circular standard deviation of 52.9° (Figure 4D and E). Shuffling body parts on each trial gave non-directional trajectories with a circular standard deviation of 126.8° (Figure 4—figure supplement 3). Behavioral trajectories also show that the response magnitude in PC space can be partly explained by initial PC1 and PC2 values (Figure 4F and G). This suggests that the initial pose influences these fast behavioral responses.

Mapping whole-body behavioral repertoires to precise nociceptive input.
(A) Example spatiotemporal structure of a noxious stimulus response superimposed on the baseline image taken immediately before stimulus. The color indicates the timing of nose and hind paw trajectories. In this example, the left hind paw of the mouse was stimulated, which is the right hind paw as viewed in the image. For ease, we refer to the stimulation side as viewed in the image, rather than the side with respect to the mouse. (B) Example graphical representation showing the sequence of postural adjustment following nociceptive stimulus in two trials. Left: the left (as viewed) hind paw was stimulated. Right: the right (as viewed) hind paw was stimulated. (C) Principal component analysis of the x and y values for six body parts – nose, left hind paw digits, left hind paw heel, right hind paw digits, right hind paw heel, and tail base – across all 80 trials. Coordinates were egocentrically aligned by the baseline pose, setting the tail base as origin and the stimulated paw on the right. This allowed the reconstruction of these locations using the first three principal components (PCs). Using the mean values of PC1, PC2, and PC3 with the stimulated hind paw indicated in blue (top); the mean values of PC2 and PC3, while varying PC1 either side of its mean by one standard deviation (middle-top); the mean values of PC1 and PC3, while varying PC2 (middle-bottom); and the mean values of PC1 and PC2, varying PC3 (bottom). (D) Behavioral trajectories of the 80 trials in PC space, showing 35–115 ms after stimulation. Only the first two PCs are shown for clarity. (E) PC vectors based on (D) show that trajectories are largely in the same direction. (F) The response magnitude (shown by colors that represent shift in PC2) varies as a function the initial pose, reduced to the first two PCs. (G) The initial PC values correlate with the shift in PC2 (left three plots). The initial PC3 value also correlates with the shift in PC3 (right). Least-squares linear fits are shown in blue and r values are Pearson’s correlation coefficients. (H) Raster plots of the distances that each tracked body part moves relative to baseline in 80 trials from 10 mice. All raster plots are sorted by maximum distances achieved by the stimulated paw within 300 ms of the stimulation. (I) Six representative traces showing the Euclidean distance between the stimulated paw and nose. (J) This expansion and shortening of Euclidean distance between the stimulated paw and the nose are shown up to 300 ms post-stimulus for all 80 trials by plotting the maximum distances as a function of the minimum distance. Corresponding rug plots (orange ticks) and a kernel density estimate (gray lines) are shown. (K) Traces showing the angle of the nose normalized to mean baseline angle between the nose and tail base. The tail base reflects the origin in these calculations. 80 trials are shown, with stimulation on the left hind paw and right hind paw (top). Average traces are shown in blue and red for left and right hind paw stimulations, respectively. Polar histograms for mean nose yaw during 300 ms post-stimulus, corresponding to the traces directly above (below). The dataset is provided in Figure 4—source data 2.
-
Figure 4—source data 1
Whole-body motion energy recorded at 40 frames/s with different size stimuli.
Stimuli (3 ms, 40 mW/mm2) were delivered at 0 ms. Data are from TRPV1Cre::ChR2 mice.
- https://cdn.elifesciences.org/articles/62026/elife-62026-fig4-data1-v2.csv
-
Figure 4—source data 2
Time courses for coordinates of six tracked body parts recorded at 400 frames/s.
Stimuli (10 ms, 40 mW/mm2, laser spot size S8 = 2.307 mm2) were delivered at 0 ms and body parts tracked with DeepLabCut. Data are from TRPV1Cre::ChR2 mice.
- https://cdn.elifesciences.org/articles/62026/elife-62026-fig4-data2-v2.csv
Examining specific features of these behaviors over a slightly longer period (300 ms) provides further insights. Displacement of each body part relative to their baseline position reveals the response timing, extent, and coordination (Figure 4H). The stimulated paw started moving at 29 ± 1 ms, the contralateral hind paw at 34 ± 4 ms, and the nose at 33 ± 2 ms (80 trials from 10 mice). With this intense stimulus, only in 6% of trials did the hind paws or single body parts move alone, although the magnitude of the head movement varied between trials. The distance traveled by the nose positively correlates with the distance for the stimulated paw (Pearson’s r = 0.64, n = 80 trials from 10 mice). Examining the relative distance between the nose and stimulated hind paw shows a reliably short latency (Figure 4I), indicating that these responses are driven by Aδ-nociceptor input rather than more slowly conducting C-fibers. A diversity of responses was observed: the head and stimulated paw move closer together in some trials and in others moved further apart (Figure 4I and J). This could result from the head moving towards or away from the stimulated paw but also the stimulated paw moving backwards as the body rotates. Indeed, consistent with initial observations (Figure 4A and B) and principal component analysis (PCA; Figure 4C–G), we find that the head selectively and rapidly orients to the stimulated side (Figure 4K). The presence of head orientation suggests that a brief nociceptive input can rapidly generate a coordinated spatially organized behavioral response. This is likely integral to protective pain-related behaviors and might function to gather sensory information about the stimulus or its consequences, and potentially provides coping strategies. Protective behaviors can be statistically categorized (Abdus-Saboor et al., 2019) and computational discrimination of high-speed hind paw responses used as a score of pain (Jones et al., 2020). We have shown that the analysis can easily be customized to incorporate computational tools that facilitate quantification and reveal insights into complex behavioral responses.
Behavioral responses to precise LTMR input
The vesicular glutamate transporter-1 (Vglut1) is a known marker of Aβ-LTMRs (Alvarez et al., 2004; Brumovsky et al., 2007). To demonstrate the utility of the system in the broader context of somatosensation, we crossed Slc17a7-IRES2-Cre-D (Vglut1Cre) mice with R26-CAG-LSL-ChR2-tdTomato mice to generate a line (Vglut1Cre::ChR2) that express ChR2 in LTMRs (Harris et al., 2014). A recent detailed anatomical and physiological characterization of Vglut1Cre::ChR2 mice further confirmed that in DRG neurons, ChR2 is restricted to broad class of myelinated Aβ-LTMRs (Chamessian et al., 2019). Here, we find that a single 3 ms stimulus (S7 = 1.155 mm2) precisely delivered to the hind paw of these mice rarely elicited hind paw responses (mean paw withdrawal probability = 0.10 ± 0.03 SEM, 99 trials from n = 11 mice), with the earliest response occurring at 206 ms after stimulation (Figure 5A and B), which is an order of magnitude slower than we observed in TRPV1Cre::ChR2 mice (fastest response: 19 ms). Trains of five pulses, however, frequently elicited responses, showing mean paw withdrawal probabilities of 0.31 ± 0.09 (SEM, 108 trials from n = 12 mice) for 5 Hz and 0.40 ± 0.10 (SEM, 117 trials from n = 12 mice) for 10 Hz trains (Figure 5C). Increasing stimulation frequency to 20 Hz did not result in higher withdrawal probabilities, which may reflect ChR2 desensitization, rather than a physiological process (Lin, 2011). While the responses at first seem to be frequency-dependent (Figure 5D, left), inspection of recordings indicated that these occurred after the second or third pulse in most trials, regardless of stimulation frequency (Figure 5A). We find that the response distributions superimpose when withdrawal latencies are normalized to the interstimulus interval (pulse-matched latencies in Figure 5D, right). This observation suggests that response probability is likely driven by pulse summation, rather than by stimulation frequency. Indeed, we find that the probabilities and latencies can be explained by the probability sum rule, using the values for a single pulse to predict the values for five pulses (Figure 5C and D).

Scanned transdermal optogenetic activation of Aβ-LTMRs triggers slow-onset responses.
(A) Example traces of the near-infrared frustrated total internal reflection signal time course for three different stimulation protocols in Vglut1Cre::ChR2 mice: single pulse, five pulses at 5 Hz, and five pulses at 10 Hz (pulse duration 3 ms, spot size S7 = 1.155 mm2). (B) Corresponding raster plots of hind paw responses sorted by latency. The paw response latency is indicated in red (99–103 trials/protocol from n = 11–12 mice) and the 3 ms laser stimuli shown with blue carets. (C) Paw response probability peaks at 10 Hz stimulation frequency in Vglut1Cre::ChR2 mice (pulse duration 3 ms, spot size S7 = 1.155 mm2; 99–103 trials/protocol from n = 11–12 mice, mean probability ± SEM). (D) Left panel: paw response latencies in trials with a single 3 ms stimulation or with trains of five 3 ms stimuli at 5 Hz or at 10 Hz. Right panel: paw response latencies normalized to the interstimulus interval. The estimated probability in (C) and (D) (dashed gray lines) was calculated using P(X ≥ 1) = 1–(1–p)n, where p is the probability of a response on a single pulse (0.096) and n is the number of pulses (5). The dataset is provided in Figure 5—source data 1.
-
Figure 5—source data 1
Time courses of paw movement recorded at 1000 frames/s with stimuli that vary in frequency.
Stimuli (3 ms, laser spot size S7 = 1.181 mm2, 40 mW/mm2) were delivered at 0 ms. Data are from Vglut1Cre::ChR2 mice and littermate controls.
- https://cdn.elifesciences.org/articles/62026/elife-62026-fig5-data1-v2.csv
-
Figure 5—source data 2
Whole-body motion energy recorded at 400 frames/s with stimuli that vary in frequency.
Stimuli (3 ms, laser spot size S7 = 1.181 mm2, 40 mW/mm2) were delivered at 0 ms. Data are from Vglut1Cre::ChR2 mice.
- https://cdn.elifesciences.org/articles/62026/elife-62026-fig5-data2-v2.csv
The magnitude of whole-body motion was not altered by increasing frequencies (Figure 5—figure supplement 1). In contrast to the TRPV1Cre::ChR2 line, whole-body behaviors in response to optogenetic stimulation of Vglut1Cre::ChR2 mice were subtle: visual inspection of high-speed whole-body behavior videos revealed that responses were mostly limited to small hind paw lifts or shifts towards the center of the body in cases where the stimulated paw was initially further away from the body. In most instances, these movements did not disturb balance or alter the animal’s posture. Interestingly, we observed that whisking and, to a lesser extent, circular movements of the upheld forepaws would precede hind paw responses and initiate as early as the first pulse, even in trials that would not proceed to withdrawal. We speculate that mice may perceive the stimulation early on, but only act on this after a delay.
Discussion
We describe a strategy for remote, precise, dynamic somatosensory input and behavioral mapping in awake unrestrained mice. The approach can remotely deliver spatiotemporally accurate optogenetic stimuli to the skin with predefined size, geometry, duration, timing, and location, while simultaneously monitoring behavior in the millisecond timescale. Microscale optogenetic stimulation can be used to simulate patterns, edges, and moving points on the skin. Responses to these precisely defined points and patterns can be mapped using machine vision approaches. The design is modular, for example, additional lasers for multicolor optogenetic control or naturalistic infrared stimuli can be added and complementary machine vision analysis approaches readily implemented. As an example, we combine this with DeepLabCut (Mathis et al., 2018), for markerless tracking of individual body parts to further dissect specific components of whole-body responses.
We validated the system in two transgenic mouse lines, providing optical control of broad-class Aδ and C-nociceptors, and Aβ-LTMRs. Advances in transcriptional profiling have identified a vast array of genetically defined primary afferent neuron populations involved in specific aspects of temperature, mechanical, and itch sensation (Usoskin et al., 2015). Selective activation of these populations is expected to recruit a specific combination of downstream cells and circuits depending on their function. For example, nociceptive input generates immediate sensorimotor responses and also pain that acts as a teaching signal. This strategy can be thus combined with techniques to modify genes, manipulate cells and neural circuits, and record neural activity in freely behaving mice to probe these mechanisms (Boyden et al., 2005; Kim et al., 2017). We provide approaches to map behavioral responses to defined afferent inputs across the spectrum of somatosensory modalities (Browne et al., 2017; Huang et al., 2019).
We find that the probabilistic recruitment of nociceptors determines the behavioral response probability, latency, and magnitude. We propose that the aggregate number of first action potentials arriving from nociceptors to the spinal cord can be utilized to optimize the timing and extent of rapid protective responses. These first action potentials could be summated by spinal neurons so that appropriate behaviors are selected based on thresholds. Resultant fast behaviors are diverse but include coordinated head orientation and body repositioning that depends on the initial pose. In contrast, responses to optogenetic activation of Aβ-LTMRs occurred with slower onset, lower probability, and resulted in more subtle whole-body movements. Using a fixed number of pulses, we find that responses from multiple Aβ-LTMR inputs can be explained by the sum rule of probabilities rather than frequency-dependence (Chamessian et al., 2019). This does not, however, rule out the tuning of responses to more spatially or temporally complex stimuli. We used broad-class Cre driver lines to selectively stimulate either nociceptors or Aβ-LTMRs, and it is possible that their respective subpopulations exploit a diversity of coding strategies. This optical approach can reveal how such subpopulation and their specific downstream circuits guide behavior.
In summary, we have developed a strategy to precisely control afferents in the skin without touching or approaching them by projecting light to optogenetically generate somatosensory input in patterns, lines, or points. This is carried out non-invasively in awake freely behaving mice in a way that is remote yet precise. Remote control of temporally and spatially precise input addresses the many limitations of manually applied contact stimuli. The timing, extent, directionality, and coordination of resultant millisecond-timescale behavioral responses can be investigated computationally with specific sensory inputs. This provides a way to map behavioral responses, circuits, and cells recruited by defined afferent inputs and dissect the neural basis of processes associated with pain and touch. This strategy thus enables the investigation of sensorimotor, perceptual, cognitive, and motivational processes that guide and shape behavior in health and disease.
Materials and methods
Reagent type (species) or resource | Designation | Source or reference | Identifiers | Additional information |
---|---|---|---|---|
Genetic reagent (Mus musculus) | R26-CAG-LSL-hChR2(H134R)-tdTomato (Ai27D) | Jackson Laboratory | Stock #: 012567RRID: IMSR_JAX:012567 | PMID:22446880 |
Genetic reagent (M. musculus) | Trpv1-IRES-Cre (TRPV1Cre) | Jackson Laboratory | Stock #: 017769RRID: IMSR_JAX:017769 | PMID:21752988 |
Genetic reagent (M. musculus) | Slc17a7-IRES2-Cre-D(Vglut1Cre) | Jackson Laboratory | Stock #: 023527RRID: IMSR_JAX:023527 | PMID:21752988 |
Software, algorithm | RStudio | RStudio http://www.rstudio.com/ | RRID:SCR_000432 | Version 1.2.5019 |
Software, algorithm | Python | Python http://www.python.org/ | RRID:SCR_008394 | Version3.6.8 |
Software, algorithm | Fiji | Fiji http://fiji.sc | RRID:SCR_002285 | Version2.0.0 |
Software, algorithm | Prism 7 | GraphPad Prism http://www.graphpad.com/ | RRID:SCR_002798 | Version 7 |
Software, algorithm | Seaborn | Seaborn http://www.seaborn.pydata.org | RRID:SCR_018132 | |
Software, algorithm | Adobe Illustrator | Adobe http://www.adobe.com | RRID:SCR_010279 | Version 24.0 |
Optical system design, components, and assembly
Request a detailed protocolOptical elements, optomechanical components, mirror galvanometers, the diode laser, LEDs, controllers, machine vision cameras, and structural parts for the optical platform are listed in the table in Supplementary file 1. These components were assembled on an aluminum breadboard as shown in the Solidworks rendering in Figure 1C. The laser was aligned to the center of all lenses and exiting the midpoint of the mirror galvanometer housing aperture when the mirrors were set to the center of their working range. A series of lenses (L1–L3) expanded the beam before focusing it on to the glass stimulation plane, on which mice are placed during experiments. The glass stimulation platform was constructed of 5-mm-thick borosilicate glass framed by aluminum extrusions. NIR-FTIR was achieved by embedding an infrared LED ribbon inside the aluminum frame adjacent to the glass edges (Roberson, D. P. et al., manuscript submitted). The non-rotating L1 lens housing was calibrated to obtain eight defined laser spot sizes, ranging from 0.0185 mm2 to 2.307 mm2, by translating this lens along the beam path at set points to defocus the laser spot at the 200 mm × 200 mm stimulation plane. The beam size can be altered manually using this rotating lens tube per design, but this is modular and could be altered by the user. To ensure a relatively flat field in the stimulation plane, the galvanometer housing aperture was placed at a distance of 400 mm from its center. In this configuration, the corners of the stimulation plane were at a distance of 424 mm from the galvanometer housing aperture and variability of the focal length was below 1.5%.
Optical power density was kept constant by altering the laser power according to the laser spot area. Neutral density (ND) filters were used so that the power at the laser aperture was above a minimum working value (≥8 mW) and to minimize potential changes in the beam profile at the stimulation plane. The laser and mirror galvanometers were controlled through a multifunction DAQ (National Instruments, USB-6211) using custom software written in LabVIEW. The software displays the NIR-FTIR camera feed, whose path through the mirror galvanometers is shared with the laser beam, so that they are always in alignment with one another. Computationally adjusting mirror galvanometer angles causes identical shifts in both the descanned NIR-FTIR image field of view and intended laser stimulation site, so that the laser can be targeted to user-identified locations. Shaped stimulation patterns were achieved by programmatically scaling the mirror galvanometer angles to the glass stimulation plane using a calibration grid array (Thorlabs, R1L3S3P). The timings of laser pulse trains were synchronized with the mirror galvanometers to computationally implement predefined shapes and lines using small angle steps that could be as short as 300 µs. The custom software also synchronized image acquisition from the two cameras, so that time-locked high-speed local paw responses were recorded (camera 1: 160 pixels × 160 pixels, 250–1000 frames/s depending on the experiment). Time-locked global whole-body responses were recorded above video-frame rate (camera 2: 664 pixels × 660 pixels, 40 frames/s) or at high speed (camera 2: 560 pixels × 540 pixels, 400 frames/s) across the entire stimulation platform.
Technical calibration and characterization of the optical system
Request a detailed protocolTo calibrate the L1 lens housing and ensure consistency of laser spot sizes across the glass stimulation platform, we designed a 13.90 ± 0.05 mm thick aluminum alignment mask. This flat aluminum mask was used to replace the glass stimulation platform and was combined with custom acrylic plates that align the aperture of a rotating scanning-slit optical beam profiler (Thorlabs, BP209-VIS/M) to nine defined coordinates at different locations covering the stimulation plane. The laser power was set to a value that approximates powers used in behavioral experiments (40 mW). The laser power was then attenuated with an ND filter to match the operating range of the beam profiler. Using Thorlabs Beam Software, Gaussian fits were used to determine x-axis and y-axis 1/e2 diameters and ellipticities for each laser spot size over three replicates at all nine coordinates. The averages of replicates were used to calculate the area of the eight different laser spot sizes that were measured in each of the nine coordinates (Figure 1—figure supplement 1A) and then fitted with a two-dimensional polynomial equation in MATLAB to create heatmaps (Figure 1—figure supplement 1 B).
The average values over the nine coordinates were defined for each laser spot size: S1 = 0.0185 mm2, S2 = 0.0416 mm2, S3 = 0.0898 mm2, S4 = 0.176 mm2, S5 = 0.308 mm2, S6 = 0.577 mm2, S7 = 1.155 mm2, S8 = 2.307 mm2. These measurements were repeated 6 months after extensive use of the optical system to ensure stability over time (Figure 1—figure supplement 1A). In addition, the uniformity of laser power was assessed by measuring optical power at five positions of the experimental platform with a power meter (Thorlabs, PM100D) (Figure 1—figure supplement 1C).
Experimental animals
Request a detailed protocolExperiments were performed using mice on a C57BL/6j background. Targeted expression of ChR2-tdTomato in broad-class cutaneous nociceptors was achieved by breeding mice homozygous for Cre-dependent ChR2(H134R)-tdTomato at the Rosa26 locus (RRID:IMSR_JAX:012567, R26-CAG-LSL-hChR2(H134R)-tdTomato, Ai27D; Madisen et al., 2012) with mice that have Cre recombinase inserted downstream of the Trpv1 gene in one allele (RRID:IMSR_JAX:017769, Trpv1-IRES-Cre, TRPV1Cre; Cavanaugh et al., 2011). Aβ-LTMRs were selectively stimulated by breeding homozygous Ai27D mice with mice in which Cre recombinase is targeted to cells expressing the vesicular glutamate transporter 1 (RRID:IMSR_JAX: 023527, Slc17a7-IRES2-Cre-D, Vglut1Cre; Harris et al., 2014). Resultant mice were heterozygous for both transgenes and were housed with control littermates that do not encode Cre recombinase but do encode Cre-dependent ChR2-tdTomato. Adult (2–4 months old) male and female mice were used in experiments. Mice were given ad libitum access to food and water and were housed in 21°C ± 2°C, 55% relative humidity and a 12 hr light:12 hr dark cycle. Experiments were carried out on at least two separate cohorts of mice, each cohort contained 4–6 mice. Experiments were spaced by at least one day in the case where the same cohort of mice was used in different experiments. All animal procedures were approved by University College London ethical review committees and conformed to UK Home Office regulations.
Optogenetic stimulation and resultant behaviors
Request a detailed protocolPrior to the first experimental day, mice underwent two habituation sessions during which each mouse was individually placed in a plexiglass chamber (100 mm × 100 mm, 130 mm tall) on a mesh wire floor for 1 hr, then on a glass platform for another hour. On the experimental day, mice were again placed on the mesh floor for 1 hr, then up to six mice were transferred to six enclosures (95 mm × 60 mm, 75 mm tall) positioned on the 200 mm × 200 mm glass stimulation platform. Mice were allowed to settle down and care was taken to stimulate mice that were calm, still, and awake in an ‘idle’ state. The laser was remotely targeted to the hind paw glabrous skin using the descanned NIR-FTIR image feed. The laser spot size was manually set using the calibrated L1 housing, while laser power and neutral density filters were used to achieve a power density of 40 mW/mm2 regardless of spot size. The software was then employed to trigger a laser pulse of defined duration (between 100 μs and 30 ms) and simultaneously acquire high-speed (1000, 500, or 250 frames/s depending on experiment) NIR-FTIR recordings of the stimulated paw, as well as a global view of the mice with a second camera (400 frames/s or 40 frames/s) (Figure 1C). Recordings of stimulations of TRPV1Cre::ChR2 mice were 1500 ms in duration, with the laser pulse initiated at 500 ms. For each stimulation protocol, six pulses, three on each hind paw, spaced by at least 1 min were delivered to eight mice, split into two cohorts. For experiments involving Vglut1Cre::ChR2 mice, we used a single stimulation spot size (S7 = 1.155 mm2) and duration (3 ms). In addition to the single-pulse stimulation, these mice received a train of five pulses applied at 5, 10, or 20 Hz. The recording time for each trial was extended to 2000 ms to accommodate for the longer stimulation period. For each protocol, Vglut1Cre::ChR2 mice were stimulated in 10 trials, split equally between the two hind paws. Data was collected from 12 Vglut1Cre::ChR2 mice and 8 littermate controls lacking Cre recombinase split into five cohorts. In all experiments, the behavioral withdrawal of the stimulated hind paw was also manually recorded by the experimenter.
Patterned stimulation protocols
Request a detailed protocolTRPV1Cre::ChR2 mice were stimulated on the heel of the hind paw with each of the following protocols: (1) a single 1 ms pulse with spot size S7 (1.155 mm2); (2) a single 1 ms pulse with spot size S4 (0.176 mm2); (3) seven 1 ms pulses with spot size S4, superimposed on the same stimulation site and spaced by 500 μs intervals; (4) seven 1 ms pulses with spot size S4, spaced by 500 μs intervals and spatially displacing stimuli with 0.3791 mm jumps such as to draw a small hexagon; (5) seven 1 ms pulses with spot size S4, spaced by 500 μs intervals and spatially displacing stimuli with 0.5687 mm jumps such as to draw a hexagon expanded by 50% compared to the previous shape; and (6) seven 1 ms pulses with spot size S4, spaced by 500 μs intervals and spatially displacing stimuli with 0.3791 mm jumps such as to draw a straight line. The power density of the stimulations was kept constant at 40 mW/mm2 as before. Seven mice, split into two cohorts, received 10 stimulations per protocol (five on each hind paw) after a baseline epoch of 500 ms. An additional cohort of four littermates lacking Cre recombinase were stimulated in the same way and served as negative controls. Finally, three TRPV1Cre::ChR2 mice were stimulated (spot size S8, 10 ms pulse duration) with a single pulse adjacent to the hind paw, five times on each side, in order to control for potential off-target effects. The NIR-FTIR signal was recorded at 500 frames/s.
Whole-body behaviors during optogenetic stimulation
Request a detailed protocolTo obtain recordings optimized for markerless tracking with DeepLabCut, a single acrylic chamber (100 mm × 100 mm, 150 mm tall) was centered on the glass stimulation platform of the system. Rapid movements were recorded at 400 frames/s using a below-view camera (FLIR, BFS-U3-04S2M-CS). Two white and two infrared LED panels illuminated the sides of the behavioral chamber in order to optimize lighting for these short exposure times and achieve high contrast images. NIR-FTIR was not used in this configuration. TRPV1Cre::ChR2 mice received between 10 and 20 single-shot laser pulse stimulations of 10 ms each, at least 1 min apart and equally split between right and left hind paw and using spot size S8 (2.31 mm2). The first 10 trials that exceeded quality control were used (see Markerless tracking of millisecond-timescale global behaviors, Data processing). Each trial consisted of a 500 ms baseline and 4000 ms after-stimulus recording epoch.
Automated analysis of optogenetically evoked local withdrawal events
Request a detailed protocolHigh-speed NIR-FTIR recordings were saved as uncompressed AVI files. A Python script was implemented in Fiji to verify the integrity of the high-speed NIR-FTIR recordings and extract average 8-bit intensity values from all frames within a circular region of interest on the stimulation site (60 pixels diameter). This output was then fed into RStudio to calculate the average intensity and associated standard deviation of the baseline recording (first 500 ms). A hind paw response was defined as a drop of intensity equal to or below the mean of the baseline minus five times its standard deviation. Paw response latency was defined as time between the start of the pulse and the time at which a hind paw response was first detected. For purposes of quality control, only recordings with a baseline NIR-FTIR intensity mean ≥ 3 and a standard deviation/mean of the baseline ratio ≥23 were retained for analysis. Another criterion was that response latencies are not 10 ms or shorter since this would be too short to be generated by the stimulus itself. Only one trial out of 2369 trials did not meet this criterion (spot size S6, 1 ms pulse, 8 ms response latency). In addition to this two-step workflow using Fiji/Python to process AVI files and then RStudio to analyze the resulting output, alternative code was written in Python 3, which combines both steps and also computes individual pixel latencies and motion energy using NumPy and Pandas packages. A median filter (radius = 2 pixels) was applied to the NIR-FTIR recordings used to create the representative time series in Figure 2A and Figure 2—video 1. For raster plots of hind paw response dynamics in Figure 4A, NIR-FTIR intensity values were normalized to the average baseline value. For the patterned stimulation experiments in Figure 2F and Vglut1Cre::ChR2 experiments in Figure 5A–D, trials were analyzed as stated to compute local response probabilities, but an additional rule was introduced to further minimize the risk of false positives. A response required the signal to fall by 20% and exceed a threshold of four times the standard deviation of baseline. Compared to the performance of an experimenter manually processing the videos with Fiji, the automated analysis pipeline was substantially faster for similar accuracy. For example, it took an experimenter two working days to analyze 127 videos, whereas the Fiji/Python pipeline generated the identical output within 90 s.
Automated analysis of whole-body protective behavior
Request a detailed protocolVideos of the entire stimulation platform were cropped into individual mouse chambers (200 × 315 pixels) and then analyzed using RStudio to quantify the amount of whole-body movements, including those stemming from the response of the stimulated limb, herein referred to as global behavior (GB). GB was approximated as the binarized motion energy: the summed number of pixels changing by more than five 8-bit values between two subsequent frames (Pixel Change). Briefly, for each pixeln (n = 63,000 pixels/frame), the 8-bit value at a given frame (Fn) was subtracted from the corresponding pixeln at the previous frame (Fn-1). If the resulting absolute value was ≤5, 0 would be assigned to the pixel. If the absolute resulting value was >5, 1 would be assigned to the pixel. The threshold was chosen to discard background noise from the recording. The pixel binary values were then summed for each frame pair to obtain binarized motion energy. Normalized binarized motion energy was calculated by subtracting each post-stimulus frame binarized motion energy from the average baseline binarized motion energy. As an alternative to this analysis strategy, we have developed code in Python that processes the video files and calculates motion energy. The peak normalized binarized motion energy was determined and only trials displaying a peak response ≥5 standard deviations of the baseline mean were retained for further analysis and plotting. For TRPV1Cre::ChR2 mice, the analysis was restricted to a time window of 100 ms after stimulus onset (first three frame pairs proceeding the stimulus frame) to enable time-locking to the stimulus. Between 41 and 47 videos from eight mice were analyzed per spot size. For experiments with Vglut1Cre::ChR2 mice, the peak normalized binary motion energy exceeding five standard deviations of the baseline mean was determined for the entire 1.5 s recording epoch proceeding stimulus onset. Between 51 and 80 trials from 11 to 12 mice were analyzed per stimulation frequency.
Markerless tracking of millisecond-timescale global behaviors
DeepLabCut installation
Request a detailed protocolDeepLabCut (version 2.0.1) was installed on a computer (Intel-Core-i7-7800 × 3.5 GHz CPU, NVIDIA GTX GeForce 1080 Ti GPU, quad-core 64 GB RAM, Windows 10, manufactured by PC Specialist Ltd.) with an Anaconda virtual environment and was coupled to Tensorflow-GPU (v.1.8.0, with CUDA v.9.01 and cUdNN v. 5.4).
Data compression
Request a detailed protocolAll recordings were automatically cropped with Python MoviePy package and compressed with standard compression using the H.264 format, then saved in mp4 format. This compression method was previously shown to result in robust improvement of processing rate with minimal compromise on detection error.
Training the network
Request a detailed protocolDeepLabCut was used with default network and training settings. Pilot stimulation trials were collected for initial training with 1,030,000 iterations from 253 labeled images from 50 videos. The videos were selected to represent the whole range of behavioral responses and conditions (25 videos of males and 25 videos of females from six different recording sessions). Out of the 25 videos, 15 were selected from the most vigorous responses, 5 were selected from less vigorous responses, and 5 from control mice. Ground truth images were selected manually, aiming to include the most variable images from each video (up to 14 frames per video). 18 body parts were labeled, namely the nose, approximate center of the mouse, two points on each sides of the torso and one point at each side of the neck, the fore paws, distal and proximal points on the hind paw, between the hind limbs, and three points on the tail. While most of these labels were not used in subsequent analysis, labeling more body parts on the image enhanced performance. The resulting network output was visually assessed. Erroneously labeled frames were manually corrected and used to retrain the network while also adding new recordings. Four sequential retraining sessions with 1,030,000 iterations each were conducted adding a total of 109 frames from 38 videos. This resulted in a reduction in the pixel root mean square error (RMSE) from 4.97 down to 2.66 on the test set, which is comparable to human ground truth variability quantified elsewhere.
Data processing
Request a detailed protocolOnly labels of interest were used for analysis. These were ipsilateral and contralateral hind paws (distal), the tail base, and the nose labels. To minimize error, points were removed if they (1) were labeled with less than 0.95 p-cutoff confidence by DeepLabCut, (2) jumped at least 10 pixels in one single frame compared to the previous frame, (3) had not returned on the subsequent frame, and (4) were from the five stimulation frames. Code for data processing was written in Python using the NumPy and Pandas packages. Additional post-hoc quality control was performed on the network output to identify and remove poorly labeled trials. To this end, heatmaps of distances between labels were created and inspected for dropped labels and sudden changes in distance. Trials identified in this manner were then manually inspected and removed if more than 10% of labels were missing or more than 10 frames were mislabeled. In total, 4.7% of trials were discarded. Only the first eight trials for each of the 10 mice that met this video quality control were used in analysis.
Automated detection of the stimulated limb
Request a detailed protocolDisabling NIR-FTIR illumination reduces the baseline saturation and thus allowed us to automate stimulated paw detection using pixel saturation from the stimulation laser. To determine which of the left or right paw had been stimulated in a given trial, the number of saturated pixels within a 60 × 60 pixels window close to the hind paw label was compared 7.5 ms prior and 5 ms after stimulus onset.
Detection of movement latency of discrete body parts
Request a detailed protocolMovement latencies of hind paws and head (nose) were computed based on significant changes from the baseline position. Baseline positions were calculated as the average x and y values from 10 consecutive frames prior to stimulus onset. A post-stimulus response was considered to be meaningful if the position of the label changed by at least 0.5 pixels (~0.16 mm) compared to baseline and continued moving at a rate of at least 0.5 pixel/frame for the subsequent 10 frames.
Dimensionality reduction
Request a detailed protocolWe carried out dimensionality reduction on x and y values for six body parts (nose, left hind paw digits, left hind paw heel, right hind paw digits, right hind paw heel, and tail base) determined at a single time point. These were egocentrically aligned using the tail base as the origin, and the stimulated paw always on the right. PCA was carried out by extracting the first three PCs using these 12 features at 115 ms after stimulus onset. The PCA was cross-validated by pseudo-randomly splitting the 80 trials into training and test datasets (80:20). The training dataset showed 49.5, 27.4, and 12.3% variance was explained by PC1, PC2, and PC3, respectively. The same PCs explained 53.5, 23.2, and 10.1% variance in the test dataset. PCA of these 80 trials together (at 115 ms) gave explained variance values 50.4% (PC1), 26.5% (PC2), and 11.9% (PC3). Projecting the time courses onto these same PCs resulted in explained variance values 37.1% (PC1), 24.3% (PC2), and 16.7% (PC3). In all cases, the shifts seen in PC1–3 were similar to that shown in Figure 4C.
Motion energy calculations in millisecond-timescale global behaviors
Request a detailed protocolGB was analyzed within a 1 ms time frame following stimulation by computing the binarized motion energy relative to a baseline reference frame 5 ms prior to stimulation as described above. Here, the threshold for pixel change was set to seven 8-bit values. The binarized motion energy (sum of pixel binaries) of a given frame was normalized to the total number of pixels within that frame after removing those frames that had been affected by the stimulation laser pulse. The global response latency of movement initiation was determined as the time when binarized motion energy was greater than 10 times the standard deviation at baseline. Termination of movement was determined as the time point when binarized motion energy returned below 10 times standard deviation from baseline following the first movement bout.
Statistical analysis
Request a detailed protocolData was analyzed in RStudio 1.2.5019, Python 3.6.8, ImageJ/Fiji 2.0.0 and Prism 7 and visualized using Seaborn, Prism 7, and Adobe Illustrator 24.0. In all experiments, repeated measurements were taken from multiple mice. Paw responses to patterned stimulation were reported as mean probabilities ± standard error of the mean (SEM) and analyzed using Friedman’s non-parametric test for within-subject repeated measures followed by Dunn’s signed-rank test for multiple comparisons (Figure 2F). In this experiment, one of the seven TRPV1Cre::ChR2 mice was removed from the dataset because it displayed saturating responses to Protocol 3 preventing comparison of values across a dynamic range. Response latencies, response rise times, and response durations were computed using a hierarchical bootstrap procedure (Saravanan et al., 2020) modified to acquire bootstrap estimates of the median with balanced resampling. Briefly, mice are sampled with replacement for the number of times that there are mice. For each mouse within this sample, its trials were sampled with replacement, but the number of selected trials was balanced, ensuring each mouse contributes equally to the number of trials in the sample. The median was taken for this resampled population and this entire process was repeated 10,000 times. Bootstrap estimates from 1000 simulated experiments show that an additional 1.6–3.1% of values fall within 1% of the population median for seven mice with between 2 and 6 responses. Values provided are the mean bootstrap estimate of the median ± the standard error of this estimate. The median bias was small due to the resampled population size from hierarchically nested data and only moderate distribution skew. Global peak motion energy (Figure 4B) was examined in a similar way, except the mean of resampled populations was used as it represents a better estimator of the population mean. In this case, we report the mean bootstrap estimate of the mean ± the standard error of this estimate. Pearson’s correlation coefficients were determined to compare maximum distances moved from baseline for each body part (Figure 4F). Experimental units and n values are indicated in the figure legends.
Data availability
All components necessary to assemble the optical system are listed in Supplementary file 1. A Solidworks assembly, the optical system control and acquisition software and behavioral analysis toolkit are available at https://github.com/browne-lab/throwinglight. The data that support the findings of this study are provided as source data files.
References
-
Vesicular glutamate transporters in the spinal cord, with special reference to sensory primary afferent synapsesThe Journal of Comparative Neurology 472:257–280.https://doi.org/10.1002/cne.20012
-
Millisecond-timescale, genetically targeted optical control of neural activityNature Neuroscience 8:1263–1268.https://doi.org/10.1038/nn1525
-
Remote optogenetic activation and sensitization of pain pathways in freely moving miceThe Journal of Neuroscience 33:18631–18640.https://doi.org/10.1523/JNEUROSCI.2424-13.2013
-
Somatosensation in the Brain: A Theoretical Re-evaluation and a New ModelTrends in Cognitive Sciences 24:529–541.https://doi.org/10.1016/j.tics.2020.04.003
-
Nociceptors: The sensors of the pain pathwayThe Journal of Clinical Investigation 120:3760–3772.https://doi.org/10.1172/JCI42843
-
Neuronal diversity in the somatosensory system: bridging the gap between cell type and functionCurrent Opinion in Neurobiology 56:167–174.https://doi.org/10.1016/j.conb.2019.03.002
-
Spatial sensory organization and body representation in pain perceptionCurrent Biology 23:R164–R176.https://doi.org/10.1016/j.cub.2013.01.047
-
Anatomical characterization of cre driver mice for neural circuit mapping and manipulationFrontiers in Neural Circuits 8:76.https://doi.org/10.3389/fncir.2014.00076
-
Integration of optogenetics with complementary methodologies in systems neuroscienceNature Reviews. Neuroscience 18:222–235.https://doi.org/10.1038/nrn.2017.15
-
Deeplabcut: Markerless pose estimation of user-defined body parts with deep learningNature Neuroscience 21:1281–1289.https://doi.org/10.1038/s41593-018-0209-y
-
Fast animal pose estimation using deep neural networksNature Methods 16:117–125.https://doi.org/10.1038/s41592-018-0234-5
-
Sensorimotor processing in the rodent barrel cortexNature Reviews. Neuroscience 20:533–546.https://doi.org/10.1038/s41583-019-0200-y
-
Normal and abnormal coding of somatosensory stimuli causing painNature Neuroscience 17:183–191.https://doi.org/10.1038/nn.3629
-
Edge-orientation processing in first-order tactile neuronsNature Neuroscience 17:1404–1409.https://doi.org/10.1038/nn.3804
-
Application of the hierarchical bootstrap to multi-level data in neuroscienceNeurons, Behavior, Data Analysis and Theory 3:819334.https://doi.org/10.1101/819334
-
Flexion-reflex of the limb, crossed extension-reflex, and reflex stepping and standingThe Journal of Physiology 40:28–121.https://doi.org/10.1113/jphysiol.1910.sp001362
-
Unbiased classification of sensory neuron types by large-scale single-cell RNA sequencingNature Neuroscience 18:145–153.https://doi.org/10.1038/nn.3881
Decision letter
-
Ishmail Abdus-SaboorReviewing Editor; University of Pennsylvania, United States
-
Andrew J KingSenior Editor; University of Oxford, United Kingdom
In the interests of transparency, eLife publishes the most substantive revision requests and the accompanying author responses.
Acceptance summary:
We acknowledge the effort and time the authors took to address all of the comments made by reviewers which have improved this paper. We strongly believe that this approach will be widely adopted in the somatosensory research community to deliver "remote optogenetic touch" to freely behaving animals in a highly spatial and temporal fashion. The addition of a non-nociceptive line to this paper, with distinct behavioral outputs from the nociceptive line, is an added major strength of this work.
Decision letter after peer review:
Thank you for submitting your article "Scanned optogenetic control of mammalian somatosensory input to map input-specific behavioral outputs" for consideration by eLife. Your article has been reviewed by 3 peer reviewers, one of whom is a member of our board of Reviewing Editors and the evaluation has been overseen by Andrew King as the Senior Editor. The reviewers have opted to remain anonymous.
The reviewers have discussed the reviews with one another and the Reviewing Editor has drafted this decision to help you prepare a revised submission.
The manuscript by Schorscher-Petcu is a very innovative study approaching an important problem in pain and somatosensory neuroscience – precise and remote delivery of sensory stimuli. This work is timely, and there are many clear applications to understanding peripheral somatosensory encoding using this strategy. The rationale for such a tool developed here is widely agreed upon in the field, and if others can easily adopt this strategy, this could become the standard for peripheral optogenetic stimulation of the hind paw. Nonetheless, the reviewers would like the authors to address the points listed here before consideration of publication in eLife.
1. The strength of this paper was in the technique and not the new biology uncovered. We would like the authors to remove any language of a "sparse code" proposition because the data presented here do not support such a claim.
2. More detail on how to use this system so that new users can use this off-the-shelf. In particular: we had a hard time evaluating the hierarchical bootstrap procedure, which references a pre-print. Is this method really ensuring that the results are more rigorous? How do the authors define a withdraw? More detail and commentary on how this approach interfaces with Deep Lab Cut. In general, focus more on the technique and less on biology.
3. For widespread applicability and to determine the range, strengths, and weaknesses of this new tool to the pain field, the authors should extend their behavioral analyses. The reviewers preferred the authors to do as they mention in the discussion, which is to add an additional somatosensory line (perhaps a non-pain line) and see how their platform performs in comparison to Trpv1-ChR2. The less preferred option if the authors are not able to breed new somatosensory lines in reasonable time, is to try the Trpv1-ChR2 line in different contexts (inflammatory and/or neuropathic pain). In either case, at baseline or during chronic pain states in the Trpv1-ChR2 line, the authors should use an analgesic and show that their tool is modular and can detect decreases in pain-related signatures. The authors should take care to have N numbers closer to 10 animals per group, as the N of 4 in their studies is on the lower side.
https://doi.org/10.7554/eLife.62026.sa1Author response
1. The strength of this paper was in the technique and not the new biology uncovered. We would like the authors to remove any language of a "sparse code" proposition because the data presented here do not support such a claim.
We thank the reviewers and agree that as a technical report we should focus on the technical aspects. In demonstrating the utility of the system, we reveal biological insights that should be described clearly and fully. The language used is critical, so based on the reviewers comments we have removed the term "sparse code". We agree that "code" unintentionally implied that we show noxious stimuli are encoded naturally in this manner in this combination of broad-class nociceptors. More precisely, we are reporting that the relative size of a single action potential nociceptor volley can determine the nature, timing and extent of resultant protective behaviors. This relationship between nociceptive input and behavior has not been demonstrated before and so it is important to describe this. We do however retain the use of the term "sparse" as this summarizes the optogenetic stimulus for the reasons below.
First, the input is spatially sparse – stimuli activate a small percentage of nociceptors in the total population. In Figure 2D, we show even an area as small as 0.6 mm² can tune the behavioral response. We now emphasize this point in the manuscript “The stimulus used for this protocol was S6, 0.577 mm² in area, which corresponds to less than 1% of the glabrous paw area and highlights the sensitivity of the nociceptive system”. Second, the input is temporally sparse – the stimuli only generate a single action potential in each nociceptor and this is enough to generate coordinated behaviors. That this is temporally sparse is in line with common usage of the term (e.g. Hahnloser et al., Nature 2002). This builds on our previous work that showed that a single pulse only generates a single action potential in each nociceptor.
We confirm that the subtle difference in laser rise times have no implications on generalizing the previous findings. Indeed, the TRPV1-lineage used here have long refractory periods (minimum interspike intervals, ~20 ms), as measured by patch clamp recordings (Browne et al., Cell Reports 2017). We had previously reported that even a 10 ms optogenetic stimulation of this mouse line does not generate more than a single action potential (see Browne et al., Cell Reports 2017 Supplementary Information, in vitro electrophysiology), but this is perhaps not prominent enough in the previous paper. In the system we describe here, the diode laser delivers the same wavelength and intensity of light as before. The light reaches a stable intensity in less than 1 microsecond, which allowed us to precisely compare the effects of short pulses as the ChR2 kinetics (tau=1.9 ms, Lin, 2011) and are what limits the rate of depolarization, not these lasers. Figure 2E uses 1 ms and 3 ms pulses only, and Figure 2F uses 1 ms pulses. We show that paw response probability increases from 0.3 ms up to 10 ms (Figure 2D), which is not long enough to generate more than one action potential. We now point out alongside Figure 2D in the legend that the 30 ms pulse duration may result in more than one action potential but it is important to note this does not alter the interpretation as the probability already plateaued at 10 ms.
Sparse recruitment of broad-class nociceptors was achieved by gaining precise control over the pulse width and spatial properties (size and shape) of the stimuli. It is not known whether certain pulse durations are more selective for the fibers they activate, so we have used a range of approaches to control stimulation size. We find that the larger the stimulus, irrespective of increasing the spot size or increasing the duration of the pulse, the shorter the response latencies. As illumination area increases so does the probability that Aδ and C fibres are activated. Our analysis of hind paw latencies was focused on Aδ responses. This is possible due to the stimuli being time-locked and the differential conduction velocities of these fibres – the responses typically initiate before 35 ms when the C-fiber action potentials (1.5 m/s) will only just be arriving at the spinal cord. We have previously shown using electrophysiology that TRPV1-Cre::ChR2 drives expression of ChR2 in broad-class nociceptors including Aδ and C fibers (Browne et al., Cell Reports, 2017). This is also in line with evidence from Usoskin et al., (Nature Neuroscience, 2015) showing the PEP2 population expresses TRPV1, and discussed by Caterina and Julius (Annual Review of Neuroscience, 2001) and Mitchell et al., (Mol Pain 2010).
2. More detail on how to use this system so that new users can use this off-the-shelf. In particular: we had a hard time evaluating the hierarchical bootstrap procedure, which references a pre-print. Is this method really ensuring that the results are more rigorous? How do the authors define a withdraw? More detail and commentary on how this approach interfaces with Deep Lab Cut. In general, focus more on the technique and less on biology.
We thank the reviewers and have made changes throughout. Specific responses and actions are given below.
Bootstrap. The bootstrap procedure allows statistical error to be calculated and provides better estimates of central tendency in skewed distributions. Here, we use hierarchical bootstrap, which as the now peer-reviewed article from Sober and colleagues points out (Saravanan et al., Neuron Behav Data Anal Theory, 2020), is an extensively used statistical method for nested data. We modified this to ensure that each mouse contributes equally to the number of trials in the sample, to avoid bias. By running our own analysis of this procedure, we confirm that the results are more rigorous; for 7 mice with between 2 and 6 responses, an additional 1.6-3.1% of bootstrap estimates from 1000 simulated experiments (with skewed distributions) fall within 1% of the population median. We have added this to the manuscript.
Targeting stimuli. The laser is fixed to the center of view in the targeting camera, which is displayed in real time in the GUI. This is due to the "descanned" configuration of the laser and targeting camera both using the same path through the galvanometers. Moving the target location of the laser moves the field of view that is displayed. The experimenter chooses what exact location of the mouse to stimulate by aiming the laser with this camera feed. We used a computer mouse cursor to rapidly target stimuli. On triggering a trial, the same camera records 0.5 seconds before the stimulus and 1.5 second after stimulus onset. These high-speed (up to 1000 fps) videos are used to verify if any movement occurred prior to the stimulation. There are occasions when the mouse moves before the stimulation is delivered, which is unavoidable in freely moving mice but the pre-stimulus baseline from all recordings are examined for movement using our code. In practical terms, we assess whether there has been movement occurring during the 0.5 seconds preceding the stimulus by calculating the mean signal intensity for this epoch and dividing it by its standard deviation. If the animal is still, the standard deviation will be small and the incurring ratio large. If the mouse moved, the standard deviation will be high and the ratio small. We empirically determined a threshold value, which we found to reliably indicate whether the animal is not moving up until the delivery of the stimulation. If there is any movement, the trial will not be used thus removing the possibility of recordings with unintended stimulus delivery.
To provide clarification on the targeting of stimuli, we have added further description to the Results section Design and assembly of the optical stimulation approach, which now states:
“The camera-feed is displayed in the user interface and enables the operator to use this image to target the laser to the desired location”.
This makes it clearer how the laser is remotely targeted. The reader may also refer to the Materials and methods section for a more detailed description:
“The software displays the NIR-FTIR camera feed, whose path through the mirror galvanometers is shared with the laser beam, so that they are always in alignment with one another. Computationally adjusting mirror galvanometer angles causes identical shifts in both the descanned NIR-FTIR image field of view and intended laser stimulation site, so that the laser can be targeted to user-identified locations.”
Paw responses. We define a hind paw response as a decrease in the intensity equal to or below the mean of the baseline minus five times its standard deviation. This is stated in Materials and methods section Automated analysis of optogenetically evoked local withdrawal events. By this definition, a paw response does not require the paw to lift / withdraw. We examine the extent of the response and can categorize partial and full responses. For example, if only part of the paw moves without lifting from the floor this is not a full paw withdrawal but still considered a response and it contains information that can be included in the analysis. The responses are detected using frustrated total internal reflection of infrared light at 1,000 frames per second. These resultant videos reveal nuance that would not be detected with the human eye. The videos are analyzed with code written in R and Python. We extract each of the 25,600 pixels for each frame and statistically analyze changes in each pixel or the intensity across a region of interest. The code is freely available. The videos are inspected by eye to validate findings. We have made changes in the manuscript that make this clearer.
Open science. We provide all information required to build the system: a solidworks assembly, part list, software and analysis code, along with technical specifications of the system in the Results section and Materials and methods sections Optical system design, components and assembly and Patterned stimulation protocols.
DeepLabCut. We have added additional details and discussion regarding the implementation of DeepLabCut to the Results section Sparse nociceptor stimulation triggers coordinated postural adjustments, and in the Discussion.
3. For widespread applicability and to determine the range, strengths, and weaknesses of this new tool to the pain field, the authors should extend their behavioral analyses. The reviewers preferred the authors to do as they mention in the discussion, which is to add an additional somatosensory line (perhaps a non-pain line) and see how their platform performs in comparison to Trpv1-ChR2. The less preferred option if the authors are not able to breed new somatosensory lines in reasonable time, is to try the Trpv1-ChR2 line in different contexts (inflammatory and/or neuropathic pain). In either case, at baseline or during chronic pain states in the Trpv1-ChR2 line, the authors should use an analgesic and show that their tool is modular and can detect decreases in pain-related signatures. The authors should take care to have N numbers closer to 10 animals per group, as the N of 4 in their studies is on the lower side.
This is an excellent suggestion. We have included an additional somatosensory Cre driver line, which took time to obtain and breed but substantially improves the manuscript. The somatosensory line we use is the Vglut1-Cre::ChR2-tdTomato line, which was originally characterized by Ru-Rong Ji and colleagues (Chamessian et al., J. Neurosci. 2019) in a relevant context. Here, we show that behavioral responses are specific to this stimulus. We therefore demonstrate the system is sensitive enough to detect behavioral responses associated with specific non-nociceptive and nociceptive input, and these responses can be readily discriminated.
Centrally acting analgesics (morphine, for example) may or may not suppress single-shot optogenetic stimuli; the interpretation of such an experiment does not address if pain is or is not experienced by such a stimulus, as they typically have other effects on arousal and locomotion that confound responses to mild inputs. A grooming or sleeping mouse, for example, does not respond to noxious stimuli as often as when it is at rest (Browne et al., Cell Reports 2017). We have previously shown that local lidocaine injection can block the optogenetic-evoked responses in these TRPV1-Cre::ChR2 mice (Browne et al., Cell Reports 2017). This blocks the input from the periphery. Further, the editor has previously developed a mouse pain scale and showed clearly that optogenetic stimulation causes pain-like responses (Abdus-Saboor et al., Cell Reports 2019). Therefore, experiments with analgesics would not contribute further to the existing understanding here. Computational classification of pain vs non-pain behaviors is an exciting area of research but is outside of the purpose of this paper. Our goal is to describe an approach for "remote touch" in freely behaving mice and provide some clear examples of its utility.
We have substantially extended our behavioral analysis of the TRPV1-Cre::ChR2 mice not to discriminate between pain vs non-pain, but based on the excellent suggestion to examine whole-body movement “…using a vector instead of a scalar”. We have extended the DeepLabCut analysis and reworked Figure 4, providing in depth computational analysis of the behavioral trajectories, demonstrating coordinated head orienting and body reposition. This further demonstrates the analysis and insights that can be achieved using this technical advance.
Regarding n numbers, all experiments used between 7 and 12 mice with one exception; the patterned stimulation experiment in Figure 2F used 6 mice, and 3 littermate controls and 4 off-target controls to confirm that these stimuli were indeed specific. These controls together show no responses, supporting what is already known with optogenetic stimulation of this mouse line. We have made the numbers of mice used clearer throughout the manuscript.
https://doi.org/10.7554/eLife.62026.sa2Article and author information
Author details
Funding
Wellcome Trust (Sir Henry Dale Fellowship 109372/Z/15/Z)
- Liam E Browne
Royal Society (Sir Henry Dale Fellowship 109372/Z/15/Z)
- Liam E Browne
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We are grateful to Dr Mehmet Fisek and Dr Adam M Packer for initial advice on the optical system, and thank Dr David P Roberson and the Woolf lab for sharing the NIR-FTIR technology. We gratefully acknowledge feedback on the manuscript from Dr Adam M Packer and Professor John N Wood. This work was supported by a Sir Henry Dale Fellowship jointly funded by the Wellcome Trust and the Royal Society (109372/Z/15/Z).
Ethics
All animal procedures were approved by University College London ethical review committees and conformed to UK Home Office regulations.
Senior Editor
- Andrew J King, University of Oxford, United Kingdom
Reviewing Editor
- Ishmail Abdus-Saboor, University of Pennsylvania, United States
Publication history
- Preprint posted: August 10, 2020 (view preprint)
- Received: August 11, 2020
- Accepted: July 28, 2021
- Accepted Manuscript published: July 29, 2021 (version 1)
- Version of Record published: September 2, 2021 (version 2)
Copyright
© 2021, Schorscher-Petcu et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 1,967
- Page views
-
- 274
- Downloads
-
- 2
- Citations
Article citation count generated by polling the highest count across the following sources: Crossref, PubMed Central, Scopus.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Neuroscience
One signature of the human brain is its ability to derive knowledge from language inputs, in addition to nonlinguistic sensory channels such as vision and touch. How does human language experience modulate the mechanism by which semantic knowledge is stored in the human brain? We investigated this question using a unique human model with varying amounts and qualities of early language exposure: early deaf adults who were born to hearing parents and had reduced early exposure and delayed acquisition of any natural human language (speech or sign), with early deaf adults who acquired sign language from birth as the control group that matches on nonlinguistic sensory experiences. Neural responses in a semantic judgment task with 90 written words that were familiar to both groups were measured using fMRI. The deaf group with reduced early language exposure, compared with the deaf control group, showed reduced semantic sensitivity, in both multivariate pattern (semantic structure encoding) and univariate (abstractness effect) analyses, in the left dorsal anterior temporal lobe (dATL). These results provide positive, causal evidence that language experience drives the neural semantic representation in the dATL, highlighting the roles of language in forming human neural semantic structures beyond nonverbal sensory experiences.
-
- Neuroscience
Across phyla, males often produce species-specific vocalizations to attract females. Although understanding the neural mechanisms underlying behavior has been challenging in vertebrates, we previously identified two anatomically distinct central pattern generators (CPGs) that drive the fast and slow clicks of male Xenopus laevis, using an ex vivo preparation that produces fictive vocalizations. Here, we extended this approach to four additional species, X. amieti, X. cliivi, X. petersii, and X. tropicalis, by developing ex vivo brain preparation from which fictive vocalizations are elicited in response to a chemical or electrical stimulus. We found that even though the courtship calls are species-specific, the CPGs used to generate clicks are conserved across species. The fast CPGs, which critically rely on reciprocal connections between the parabrachial nucleus and the nucleus ambiguus, are conserved among fast-click species, and slow CPGs are shared among slow-click species. In addition, our results suggest that testosterone plays a role in organizing fast CPGs in fast-click species, but not in slow-click species. Moreover, fast CPGs are not inherited by all species but monopolized by fast-click species. The results suggest that species-specific calls of the genus Xenopus have evolved by utilizing conserved slow and/or fast CPGs inherited by each species.