TWISP: A Transgenic Worm for Interrogating Signal Propagation in C. elegans.

Genetically encoded optical indicators and actuators of neural activity allow for all-optical investigations of signaling in the nervous system. But commonly used indicators, actuators and expression strategies are poorly suited for systematic measurements of signal propagation at brain scale and cellular resolution. Large scale measurements of the brain require indicators and actuators with compatible excitation spectra to avoid optical crosstalk. They must be highly expressed in every neuron but at the same time avoid lethality and permit the animal to reach adulthood. Their expression must also be compatible with additional fluorescent labels to locate and identify neurons, such as those in the NeuroPAL cell identification system. We present TWISP, a Transgenic Worm for Interrogating Signal Propagation, that addresses these needs and enables optical measurements of evoked calcium activity at brain scale and cellular resolution in the nervous system of the nematode Caenorhabditis elegans. We express in every neuron a non-conventional optical actuator, the gustatory receptor homolog GUR-3 + PRDX-2 under the control of a drug-inducible system QF + hGR, and calcium indicator GCAMP6 s, in a background with additional fluorophores of the NeuroPAL cell ID system. We show that this combination, but not others tested, avoids optical-crosstalk, creates strong expression in the adult, and generates stable transgenic lines for systematic measurements of signal propagation in the worm brain.

An inhibitory acetylcholine receptor gates context dependent mechanosensory processing in C. elegans.

An animal's current behavior influences its response to sensory stimuli, but the molecular and circuit-level mechanisms of this context-dependent decision-making is not well understood. In the nematode C. elegans, inhibitory feedback from turning associated neurons alter downstream mechanosensory processing to gate the animal's response to stimuli depending on whether the animal is turning or moving forward [1-3]. Until now, the specific neurons and receptors that mediate this inhibitory feedback were not known. We use genetic manipulations, single-cell rescue experiments and high-throughput closed-loop optogenetic perturbations during behavior to reveal the specific neuron and receptor responsible for receiving inhibition and altering sensorimotor processing. An inhibitory acetylcholine gated chloride channel comprised of lgc-47 and acc-1 expressed in neuron RIM receives inhibitory signals from turning neurons and performs the gating that disrupts the worm's mechanosensory evoked reversal response.

Light evokes stereotyped global brain dynamics in Caenorhabditis elegans.

Stereotyped oscillations in population neural activity recordings from immobilized Caenorhabditis elegans have garnered interest for their striking low dimensionality and their evocative state-space trajectories or manifolds. Previously these oscillations have been interpreted as intrinsically driven global motor commands. Here we test whether these oscillations are intrinsic. We show that similar oscillations are evoked by high-intensity blue light commonly used for calcium imaging. Oscillations are reduced or absent and have a lower frequency when a longer imaging wavelength is used. Under the original blue light illumination, oscillations are reduced or have a lower frequency in animals that lack GUR-3, an endogenous light- and hydrogen-peroxide-sensitive gustatory receptor. Additional experiments with hydrogen peroxide are consistent with GUR-3's involvement. We therefore propose that blue light evokes global oscillations in part through the creation of reactive oxygen species that activate the hydrogen-peroxide-sensing receptor GUR-3.

Olfactory learning alters navigation strategies and behavioral variability in C. elegans.

Animals adjust their behavioral response to sensory input adaptively depending on past experiences. The flexible brain computation is crucial for survival and is of great interest in neuroscience. The nematode C. elegans modulates its navigation behavior depending on the association of odor butanone with food (appetitive training) or starvation (aversive training), and will then climb up the butanone gradient or ignore it, respectively. However, the exact change in navigation strategy in response to learning is still unknown. Here we study the learned odor navigation in worms by combining precise experimental measurement and a novel descriptive model of navigation. Our model consists of two known navigation strategies in worms: biased random walk and weathervaning. We infer weights on these strategies by applying the model to worm navigation trajectories and the exact odor concentration it experiences. Compared to naive worms, appetitive trained worms up-regulate the biased random walk strategy, and aversive trained worms down-regulate the weathervaning strategy. The statistical model provides prediction with $>90 \%$ accuracy of the past training condition given navigation data, which outperforms the classical chemotaxis metric. We find that the behavioral variability is altered by learning, such that worms are less variable after training compared to naive ones. The model further predicts the learning-dependent response and variability under optogenetic perturbation of the olfactory neuron AWC$^\mathrm{ON}$. Lastly, we investigate neural circuits downstream from AWC$^\mathrm{ON}$ that are differentially recruited for learned odor-guided navigation. Together, we provide a new paradigm to quantify flexible navigation algorithms and pinpoint the underlying neural substrates.

Neural signal propagation atlas of Caenorhabditis elegans.

Establishing how neural function emerges from network properties is a fundamental problem in neuroscience1. Here, to better understand the relationship between the structure and the function of a nervous system, we systematically measure signal propagation in 23,433 pairs of neurons across the head of the nematode Caenorhabditis elegans by direct optogenetic activation and simultaneous whole-brain calcium imaging. We measure the sign (excitatory or inhibitory), strength, temporal properties and causal direction of signal propagation between these neurons to create a functional atlas. We find that signal propagation differs from model predictions that are based on anatomy. Using mutants, we show that extrasynaptic signalling not visible from anatomy contributes to this difference. We identify many instances of dense-core-vesicle-dependent signalling, including on timescales of less than a second, that evoke acute calcium transients-often where no direct wired connection exists but where relevant neuropeptides and receptors are expressed. We propose that, in such cases, extrasynaptically released neuropeptides serve a similar function to that of classical neurotransmitters. Finally, our measured signal propagation atlas better predicts the neural dynamics of spontaneous activity than do models based on anatomy. We conclude that both synaptic and extrasynaptic signalling drive neural dynamics on short timescales, and that measurements of evoked signal propagation are crucial for interpreting neural function.

Inhibitory feedback from the motor circuit gates mechanosensory processing in Caenorhabditis elegans.

Animals must integrate sensory cues with their current behavioral context to generate a suitable response. How this integration occurs is poorly understood. Previously, we developed high-throughput methods to probe neural activity in populations of Caenorhabditis elegans and discovered that the animal's mechanosensory processing is rapidly modulated by the animal's locomotion. Specifically, we found that when the worm turns it suppresses its mechanosensory-evoked reversal response. Here, we report that C. elegans use inhibitory feedback from turning-associated neurons to provide this rapid modulation of mechanosensory processing. By performing high-throughput optogenetic perturbations triggered on behavior, we show that turning-associated neurons SAA, RIV, and/or SMB suppress mechanosensory-evoked reversals during turns. We find that activation of the gentle-touch mechanosensory neurons or of any of the interneurons AIZ, RIM, AIB, and AVE during a turn is less likely to evoke a reversal than activation during forward movement. Inhibiting neurons SAA, RIV, and SMB during a turn restores the likelihood with which mechanosensory activation evokes reversals. Separately, activation of premotor interneuron AVA evokes reversals regardless of whether the animal is turning or moving forward. We therefore propose that inhibitory signals from SAA, RIV, and/or SMB gate mechanosensory signals upstream of neuron AVA. We conclude that C. elegans rely on inhibitory feedback from the motor circuit to modulate its response to sensory stimuli on fast timescales. This need for motor signals in sensory processing may explain the ubiquity in many organisms of motor-related neural activity patterns seen across the brain, including in sensory processing areas.

TWISP: A Transgenic Worm for Interrogating Signal Propagation in C. elegans.

Genetically encoded optical indicators and actuators of neural activity allow for all-optical investigations of signaling in the nervous system. But commonly used indicators, actuators and expression strategies are poorly suited for systematic measurements of signal propagation at brain scale and cellular resolution. Large scale measurements of the brain require indicators and actuators with compatible excitation spectra to avoid optical crosstalk. They must be highly expressed in every neuron but at the same time avoid lethality and permit the animal to reach adulthood. And finally, their expression must be compatible with additional fluorescent labels to locate and identify neurons, such as those in the NeuroPAL cell identification system. We present TWISP, a Transgenic Worm for Interrogating Signal Propagation, that address these needs and enables optical measurements of evoked calcium activity at brain scale and cellular resolution in the nervous system of the nematode Caenorhabditis elegans. We express in every neuron a non-conventional optical actuator, the gustatory receptor homolog GUR-3+PRDX-2 under the control of a drug-inducible system QF+hGR, and calcium indicator GCAMP6s, in a background with additional fluorophores of the NeuroPAL cell ID system. We show that this combination, but not others tested, avoids optical-crosstalk, creates strong expression in the adult, and generates stable transgenic lines for systematic measurements of signal propagation in the worm brain.

Continuous odor profile monitoring to study olfactory navigation in small animals.

Olfactory navigation is observed across species and plays a crucial role in locating resources for survival. In the laboratory, understanding the behavioral strategies and neural circuits underlying odor-taxis requires a detailed understanding of the animal's sensory environment. For small model organisms like Caenorhabditis elegans and larval Drosophila melanogaster, controlling and measuring the odor environment experienced by the animal can be challenging, especially for airborne odors, which are subject to subtle effects from airflow, temperature variation, and from the odor's adhesion, adsorption, or reemission. Here, we present a method to control and measure airborne odor concentration in an arena compatible with an agar substrate. Our method allows continuous controlling and monitoring of the odor profile while imaging animal behavior. We construct stationary chemical landscapes in an odor flow chamber through spatially patterned odorized air. The odor concentration is measured with a spatially distributed array of digital gas sensors. Careful placement of the sensors allows the odor concentration across the arena to be continuously inferred in space and monitored through time. We use this approach to measure the odor concentration that each animal experiences as it undergoes chemotaxis behavior and report chemotaxis strategies for C. elegans and D. melanogaster larvae populations as they navigate spatial odor landscapes.

Correcting motion induced fluorescence artifacts in two-channel neural imaging.

Imaging neural activity in a behaving animal presents unique challenges in part because motion from an animal's movement creates artifacts in fluorescence intensity time-series that are difficult to distinguish from neural signals of interest. One approach to mitigating these artifacts is to image two channels simultaneously: one that captures an activity-dependent fluorophore, such as GCaMP, and another that captures an activity-independent fluorophore such as RFP. Because the activity-independent channel contains the same motion artifacts as the activity-dependent channel, but no neural signals, the two together can be used to identify and remove the artifacts. However, existing approaches for this correction, such as taking the ratio of the two channels, do not account for channel-independent noise in the measured fluorescence. Here, we present Two-channel Motion Artifact Correction (TMAC), a method which seeks to remove artifacts by specifying a generative model of the two channel fluorescence that incorporates motion artifact, neural activity, and noise. We use Bayesian inference to infer latent neural activity under this model, thus reducing the motion artifact present in the measured fluorescence traces. We further present a novel method for evaluating ground-truth performance of motion correction algorithms by comparing the decodability of behavior from two types of neural recordings; a recording that had both an activity-dependent fluorophore and an activity-independent fluorophore (GCaMP and RFP) and a recording where both fluorophores were activity-independent (GFP and RFP). A successful motion correction method should decode behavior from the first type of recording, but not the second. We use this metric to systematically compare five models for removing motion artifacts from fluorescent time traces. We decode locomotion from a GCaMP expressing animal 20x more accurately on average than from control when using TMAC inferred activity and outperforms all other methods of motion correction tested, the best of which were ~8x more accurate than control.

Large-scale neural recordings call for new insights to link brain and behavior.

Neuroscientists today can measure activity from more neurons than ever before, and are facing the challenge of connecting these brain-wide neural recordings to computation and behavior. In the present review, we first describe emerging tools and technologies being used to probe large-scale brain activity and new approaches to characterize behavior in the context of such measurements. We next highlight insights obtained from large-scale neural recordings in diverse model systems, and argue that some of these pose a challenge to traditional theoretical frameworks. Finally, we elaborate on existing modeling frameworks to interpret these data, and argue that the interpretation of brain-wide neural recordings calls for new theoretical approaches that may depend on the desired level of understanding. These advances in both neural recordings and theory development will pave the way for critical advances in our understanding of the brain.

A high-throughput method to deliver targeted optogenetic stimulation to moving C. elegans populations.

We present a high-throughput optogenetic illumination system capable of simultaneous closed-loop light delivery to specified targets in populations of moving Caenorhabditis elegans. The instrument addresses three technical challenges: It delivers targeted illumination to specified regions of the animal's body such as its head or tail; it automatically delivers stimuli triggered upon the animal's behavior; and it achieves high throughput by targeting many animals simultaneously. The instrument was used to optogenetically probe the animal's behavioral response to competing mechanosensory stimuli in the the anterior and posterior gentle touch receptor neurons. Responses to more than 43,418 stimulus events from a range of anterior-posterior intensity combinations were measured. The animal's probability of sprinting forward in response to a mechanosensory stimulus depended on both the anterior and posterior stimulation intensity, while the probability of reversing depended primarily on the anterior stimulation intensity. We also probed the animal's response to mechanosensory stimulation during the onset of turning, a relatively rare behavioral event, by delivering stimuli automatically when the animal began to turn. Using this closed-loop approach, over 9,700 stimulus events were delivered during turning onset at a rate of 9.2 events per worm hour, a greater than 25-fold increase in throughput compared to previous investigations. These measurements validate with greater statistical power previous findings that turning acts to gate mechanosensory evoked reversals. Compared to previous approaches, the current system offers targeted optogenetic stimulation to specific body regions or behaviors with many fold increases in throughput to better constrain quantitative models of sensorimotor processing.

Fast deep neural correspondence for tracking and identifying neurons in C. elegans using semi-synthetic training.

We present an automated method to track and identify neurons in C. elegans, called 'fast Deep Neural Correspondence' or fDNC, based on the transformer network architecture. The model is trained once on empirically derived semi-synthetic data and then predicts neural correspondence across held-out real animals. The same pre-trained model both tracks neurons across time and identifies corresponding neurons across individuals. Performance is evaluated against hand-annotated datasets, including NeuroPAL (Yemini et al., 2021). Using only position information, the method achieves 79.1% accuracy at tracking neurons within an individual and 64.1% accuracy at identifying neurons across individuals. Accuracy at identifying neurons across individuals is even higher (78.2%) when the model is applied to a dataset published by another group (Chaudhary et al., 2021). Accuracy reaches 74.7% on our dataset when using color information from NeuroPAL. Unlike previous methods, fDNC does not require straightening or transforming the animal into a canonical coordinate system. The method is fast and predicts correspondence in 10 ms making it suitable for future real-time applications.

Understanding the intricacies of the brain often requires spotting and tracking specific neurons over time and across different individuals. For instance, scientists may need to precisely monitor the activity of one neuron even as the brain moves and deforms; or they may want to find universal patterns by comparing signals from the same neuron across different individuals. Both tasks require matching which neuron is which in different images and amongst a constellation of cells. This is theoretically possible in certain 'model' animals where every single neuron is known and carefully mapped out. Still, it remains challenging: neurons move relative to one another as the animal changes posture, and the position of a cell is also slightly different between individuals. Sophisticated computer algorithms are increasingly used to tackle this problem, but they are far too slow to track neural signals as real-time experiments unfold. To address this issue, Yu et al. designed a new algorithm based on the Transformer, an artificial neural network originally used to spot relationships between words in sentences. To learn relationships between neurons, the algorithm was fed hundreds of thousands of 'semi-synthetic' examples of constellations of neurons. Instead of painfully collated actual experimental data, these datasets were created by a simulator based on a few simple measurements. Testing the new algorithm on the tiny worm Caenorhabditis elegans revealed that it was faster and more accurate, finding corresponding neurons in about 10ms. The work by Yu et al. demonstrates the power of using simulations rather than experimental data to train artificial networks. The resulting algorithm can be used immediately to help study how the brain of C. elegans makes decisions or controls movements. Ultimately, this research could allow brain-machine interfaces to be developed.

Decoding locomotion from population neural activity in moving C. elegans.

We investigated the neural representation of locomotion in the nematode C. elegans by recording population calcium activity during movement. We report that population activity more accurately decodes locomotion than any single neuron. Relevant signals are distributed across neurons with diverse tunings to locomotion. Two largely distinct subpopulations are informative for decoding velocity and curvature, and different neurons' activities contribute features relevant for different aspects of a behavior or different instances of a behavioral motif. To validate our measurements, we labeled neurons AVAL and AVAR and found that their activity exhibited expected transients during backward locomotion. Finally, we compared population activity during movement and immobilization. Immobilization alters the correlation structure of neural activity and its dynamics. Some neurons positively correlated with AVA during movement become negatively correlated during immobilization and vice versa. This work provides needed experimental measurements that inform and constrain ongoing efforts to understand population dynamics underlying locomotion in C. elegans.

Nonequilibrium Green's Functions for Functional Connectivity in the Brain.

A theoretical framework describing the set of interactions between neurons in the brain, or functional connectivity, should include dynamical functions representing the propagation of signal from one neuron to another. Green's functions and response functions are natural candidates for this but, while they are conceptually very useful, they are usually defined only for linear time-translationally invariant systems. The brain, instead, behaves nonlinearly and in a time-dependent way. Here, we use nonequilibrium Green's functions to describe the time-dependent functional connectivity of a continuous-variable network of neurons. We show how the connectivity is related to the measurable response functions, and provide two illustrative examples via numerical calculations, inspired from Caenorhabditis elegans.

Measuring and modeling whole-brain neural dynamics in Caenorhabditis elegans.

The compact nervous system of the nematode Caenorhabditis elegans makes it a powerful playground to study how neural dynamics constrained by neuroanatomy generate neural function and behavior. The ability to record neural activity from the whole brain simultaneously in this worm has opened several research avenues and is providing insights into brain-wide neural coding of locomotion, sleep, and other behaviors. We review these findings and the development of new methods, including new microscopes, new genetic tools, and new modeling approaches. We conclude with a discussion of the role of theory in interpreting or driving new experiments in C. elegans and potential paths forward.

Searching for collective behavior in a small brain.

In large neuronal networks, it is believed that functions emerge through the collective behavior of many interconnected neurons. Recently, the development of experimental techniques that allow simultaneous recording of calcium concentration from a large fraction of all neurons in Caenorhabditis elegans-a nematode with 302 neurons-creates the opportunity to ask whether such emergence is universal, reaching down to even the smallest brains. Here, we measure the activity of 50+ neurons in C. elegans, and analyze the data by building the maximum entropy model that matches the mean activity and pairwise correlations among these neurons. To capture the graded nature of the cells' responses, we assign each cell multiple states. These models, which are equivalent to a family of Potts glasses, successfully predict higher statistical structure in the network. In addition, these models exhibit signatures of collective behavior: the state of single cells can be predicted from the state of the rest of the network; the network, despite being sparse in a way similar to the structural connectome, distributes its response globally when locally perturbed; the distribution over network states has multiple local maxima, as in models of memory; and the parameters that describe the real network are close to a critical surface in this family of models.

Temporal processing and context dependency in Caenorhabditis elegans response to mechanosensation.

A quantitative understanding of how sensory signals are transformed into motor outputs places useful constraints on brain function and helps to reveal the brain's underlying computations. We investigate how the nematode Caenorhabditis elegans responds to time-varying mechanosensory signals using a high-throughput optogenetic assay and automated behavior quantification. We find that the behavioral response is tuned to temporal properties of mechanosensory signals, such as their integral and derivative, that extend over many seconds. Mechanosensory signals, even in the same neurons, can be tailored to elicit different behavioral responses. Moreover, we find that the animal's response also depends on its behavioral context. Most dramatically, the animal ignores all tested mechanosensory stimuli during turns. Finally, we present a linear-nonlinear model that predicts the animal's behavioral response to stimulus.

Automatically tracking neurons in a moving and deforming brain.

Advances in optical neuroimaging techniques now allow neural activity to be recorded with cellular resolution in awake and behaving animals. Brain motion in these recordings pose a unique challenge. The location of individual neurons must be tracked in 3D over time to accurately extract single neuron activity traces. Recordings from small invertebrates like C. elegans are especially challenging because they undergo very large brain motion and deformation during animal movement. Here we present an automated computer vision pipeline to reliably track populations of neurons with single neuron resolution in the brain of a freely moving C. elegans undergoing large motion and deformation. 3D volumetric fluorescent images of the animal's brain are straightened, aligned and registered, and the locations of neurons in the images are found via segmentation. Each neuron is then assigned an identity using a new time-independent machine-learning approach we call Neuron Registration Vector Encoding. In this approach, non-rigid point-set registration is used to match each segmented neuron in each volume with a set of reference volumes taken from throughout the recording. The way each neuron matches with the references defines a feature vector which is clustered to assign an identity to each neuron in each volume. Finally, thin-plate spline interpolation is used to correct errors in segmentation and check consistency of assigned identities. The Neuron Registration Vector Encoding approach proposed here is uniquely well suited for tracking neurons in brains undergoing large deformations. When applied to whole-brain calcium imaging recordings in freely moving C. elegans, this analysis pipeline located 156 neurons for the duration of an 8 minute recording and consistently found more neurons more quickly than manual or semi-automated approaches.

Whole-brain calcium imaging with cellular resolution in freely behaving Caenorhabditis elegans.

The ability to acquire large-scale recordings of neuronal activity in awake and unrestrained animals is needed to provide new insights into how populations of neurons generate animal behavior. We present an instrument capable of recording intracellular calcium transients from the majority of neurons in the head of a freely behaving Caenorhabditis elegans with cellular resolution while simultaneously recording the animal's position, posture, and locomotion. This instrument provides whole-brain imaging with cellular resolution in an unrestrained and behaving animal. We use spinning-disk confocal microscopy to capture 3D volumetric fluorescent images of neurons expressing the calcium indicator GCaMP6s at 6 head-volumes/s. A suite of three cameras monitor neuronal fluorescence and the animal's position and orientation. Custom software tracks the 3D position of the animal's head in real time and two feedback loops adjust a motorized stage and objective to keep the animal's head within the field of view as the animal roams freely. We observe calcium transients from up to 77 neurons for over 4 min and correlate this activity with the animal's behavior. We characterize noise in the system due to animal motion and show that, across worms, multiple neurons show significant correlations with modes of behavior corresponding to forward, backward, and turning locomotion.

Simultaneous optogenetic manipulation and calcium imaging in freely moving C. elegans.

Understanding how an organism's nervous system transforms sensory input into behavioral outputs requires recording and manipulating its neural activity during unrestrained behavior. Here we present an instrument to simultaneously monitor and manipulate neural activity while observing behavior in a freely moving animal, the nematode Caenorhabditis elegans. Neural activity is recorded optically from cells expressing a calcium indicator, GCaMP3. Neural activity is manipulated optically by illuminating targeted neurons expressing the optogenetic protein Channelrhodopsin. Real-time computer vision software tracks the animal's behavior and identifies the location of targeted neurons in the nematode as it crawls. Patterned illumination from a DMD is used to selectively illuminate subsets of neurons for either calcium imaging or optogenetic stimulation. Real-time computer vision software constantly updates the illumination pattern in response to the worm's movement and thereby allows for independent optical recording or activation of different neurons in the worm as it moves freely. We use the instrument to directly observe the relationship between sensory neuron activation, interneuron dynamics and locomotion in the worm's mechanosensory circuit. We record and compare calcium transients in the backward locomotion command interneurons AVA, in response to optical activation of the anterior mechanosensory neurons ALM, AVM or both.