A recurrent neural network model of C. elegans responses to aversive stimuli

Revealing neural dynamical structure of C. elegans with deep learning

Caenorhabditis elegans serves as a common model for investigating neural dynamics and functions of biological neural networks. Data-driven approaches have been employed in reconstructing neural dynamics. However, challenges remain regarding the curse of high-dimensionality and stochasticity in realistic systems. In this study, we develop a deep neural network (DNN) approach to reconstruct the neural dynamics of C. elegans and study neural mechanisms for locomotion. Our model identifies two limit cycles in the neural activity space: one underpins basic pirouette behavior, essential for navigation, and the other introduces extra Ω turns. The combination of two limit cycles elucidates predominant locomotion patterns in neural imaging data. The corresponding energy landscape explains the switching strategies between two limit cycles, quantitatively, and provides testable predictions on neural functions and circuit roles. Our work provides a general approach to study neural dynamics by combining imaging data and stochastic modeling.

Biophysical modeling of the whole-cell dynamics of C. elegans motor and interneurons families

The nematode Caenorhabditis elegans is a widely used model organism for neuroscience. Although its nervous system has been fully reconstructed, the physiological bases of single-neuron functioning are still poorly explored. Recently, many efforts have been dedicated to measuring signals from C. elegans neurons, revealing a rich repertoire of dynamics, including bistable responses, graded responses, and action potentials. Still, biophysical models able to reproduce such a broad range of electrical responses lack. Realistic electrophysiological descriptions started to be developed only recently, merging gene expression data with electrophysiological recordings, but with a large variety of cells yet to be modeled. In this work, we contribute to filling this gap by providing biophysically accurate models of six classes of C. elegans neurons, the AIY, RIM, and AVA interneurons, and the VA, VB, and VD motor neurons. We test our models by comparing computational and experimental time series and simulate knockout neurons, to identify the biophysical mechanisms at the basis of inter and motor neuron functioning. Our models represent a step forward toward the modeling of C. elegans neuronal networks and virtual experiments on the nematode nervous system.

See Elegans: Simple-to-use, accurate, and automatic 3D detection of neural activity from densely packed neurons

In the emerging field of whole-brain imaging at single-cell resolution, which represents one of the new frontiers to investigate the link between brain activity and behavior, the nematode Caenorhabditis elegans offers one of the most characterized models for systems neuroscience. Whole-brain recordings consist of 3D time series of volumes that need to be processed to obtain neuronal traces. Current solutions for this task are either computationally demanding or limited to specific acquisition setups. Here, we propose See Elegans, a direct programming algorithm that combines different techniques for automatic neuron segmentation and tracking without the need for the RFP channel, and we compare it with other available algorithms. While outperforming them in most cases, our solution offers a novel method to guide the identification of a subset of head neurons based on position and activity. The built-in interface allows the user to follow and manually curate each of the processing steps. See Elegans is thus a simple-to-use interface aimed at speeding up the post-processing of volumetric calcium imaging recordings while maintaining a high level of accuracy and low computational demands. (Contact: enrico.lanza@iit.it).

A recurrent Hopfield network for estimating meso-scale effective connectivity in MEG

The architecture of communication within the brain, represented by the human connectome, has gained a paramount role in the neuroscience community. Several features of this communication, e.g., the frequency content, spatial topology, and temporal dynamics are currently well established. However, identifying generative models providing the underlying patterns of inhibition/excitation is very challenging. To address this issue, we present a novel generative model to estimate large-scale effective connectivity from MEG. The dynamic evolution of this model is determined by a recurrent Hopfield neural network with asymmetric connections, and thus denoted Recurrent Hopfield Mass Model (RHoMM). Since RHoMM must be applied to binary neurons, it is suitable for analyzing Band Limited Power (BLP) dynamics following a binarization process. We trained RHoMM to predict the MEG dynamics through a gradient descent minimization and we validated it in two steps. First, we showed a significant agreement between the similarity of the effective connectivity patterns and that of the interregional BLP correlation, demonstrating RHoMM's ability to capture individual variability of BLP dynamics. Second, we showed that the simulated BLP correlation connectomes, obtained from RHoMM evolutions of BLP, preserved some important topological features, e.g, the centrality of the real data, assuring the reliability of RHoMM. Compared to other biophysical models, RHoMM is based on recurrent Hopfield neural networks, thus, it has the advantage of being data-driven, less demanding in terms of hyperparameters and scalable to encompass large-scale system interactions. These features are promising for investigating the dynamics of inhibition/excitation at different spatial scales.

Learning the dynamics of realistic models of C. elegans nervous system with recurrent neural networks

Given the inherent complexity of the human nervous system, insight into the dynamics of brain activity can be gained from studying smaller and simpler organisms. While some of the potential target organisms are simple enough that their behavioural and structural biology might be well-known and understood, others might still lead to computationally intractable models that require extensive resources to simulate. Since such organisms are frequently only acting as proxies to further our understanding of underlying phenomena or functionality, often one is not interested in the detailed evolution of every single neuron in the system. Instead, it is sufficient to observe the subset of neurons that capture the effect that the profound nonlinearities of the neuronal system have in response to different stimuli. In this paper, we consider the well-known nematode Caenorhabditis elegans and seek to investigate the possibility of generating lower complexity models that capture the system's dynamics with low error using only measured or simulated input-output information. Such models are often termed black-box models. We show how the nervous system of C. elegans can be modelled and simulated with data-driven models using different neural network architectures. Specifically, we target the use of state-of-the-art recurrent neural network architectures such as Long Short-Term Memory and Gated Recurrent Units and compare these architectures in terms of their properties and their accuracy (Root Mean Square Error), as well as the complexity of the resulting models. We show that Gated Recurrent Unit models with a hidden layer size of 4 are able to accurately reproduce the system response to very different stimuli. We furthermore explore the relative importance of their inputs as well as scalability to more scenarios.

Systematic generation of biophysically detailed models with generalization capability for non-spiking neurons

Unlike spiking neurons which compress continuous inputs into digital signals for transmitting information via action potentials, non-spiking neurons modulate analog signals through graded potential responses. Such neurons have been found in a large variety of nervous tissues in both vertebrate and invertebrate species, and have been proven to play a central role in neuronal information processing. If general and vast efforts have been made for many years to model spiking neurons using conductance-based models (CBMs), very few methods have been developed for non-spiking neurons. When a CBM is built to characterize the neuron behavior, it should be endowed with generalization capabilities (i.e. the ability to predict acceptable neuronal responses to different novel stimuli not used during the model’s building). Yet, since CBMs contain a large number of parameters, they may typically suffer from a lack of such a capability. In this paper, we propose a new systematic approach based on multi-objective optimization which builds general non-spiking models with generalization capabilities. The proposed approach only requires macroscopic experimental data from which all the model parameters are simultaneously determined without compromise. Such an approach is applied on three non-spiking neurons of the nematode Caenorhabditis elegans (C. elegans), a well-known model organism in neuroscience that predominantly transmits information through non-spiking signals. These three neurons, arbitrarily labeled by convention as RIM, AIY and AFD, represent, to date, the three possible forms of non-spiking neuronal responses of C. elegans.

Neural model generating klinotaxis behavior accompanied by a random walk based on C. elegans connectome

Klinotaxis is a strategy of chemotaxis behavior in Caenorhabditis elegans (C. elegans), and random walking is evident during its locomotion. As yet, the understanding of the neural mechanisms underlying these behaviors has remained limited. In this study, we present a connectome-based simulation model of C. elegans to concurrently realize realistic klinotaxis and random walk behaviors and explore their neural mechanisms. First, input to the model is derived from an ASE sensory neuron model in which the all-or-none depolarization characteristic of ASEL neuron is incorporated for the first time. Then, the neural network is evolved by an evolutionary algorithm; klinotaxis emerged spontaneously. We identify a plausible mechanism of klinotaxis in this model. Next, we propose the liquid synapse according to the stochastic nature of biological synapses and introduce it into the model. Adopting this, the random walk is generated autonomously by the neural network, providing a new hypothesis as to the neural mechanism underlying the random walk. Finally, simulated ablation results are fairly consistent with the biological conclusion, suggesting the similarity between our model and the biological network. Our study is a useful step forward in behavioral simulation and understanding the neural mechanisms of behaviors in C. elegans.

C. elegans-based chemosensation strategy for the early detection of cancer metabolites in urine samples

Chemosensory receptors play a crucial role in distinguishing the wide range of volatile/soluble molecules by binding them with high accuracy. Chemosensation is the main sensory modality in organisms lacking long-range sensory mechanisms like vision/hearing. Despite its low number of sensory neurons, the nematode Caenorhabditis elegans possesses several chemosensory receptors, allowing it to detect about as many odorants as mammals. Here, we show that C. elegans displays attraction towards urine samples of women with breast cancer, avoiding control ones. Behavioral assays on animals lacking AWC sensory neurons demonstrate the relevance of these neurons in sensing cancer odorants: calcium imaging on AWC increases the accuracy of the discrimination (97.22%). Also, chemotaxis assays on animals lacking GPCRs expressed in AWC allow to identify receptors involved in binding cancer metabolites, suggesting that an alteration of a few metabolites is sufficient for the cancer discriminating behavior of C. elegans, which may help identify a fundamental fingerprint of breast cancer.

A Shearless Microfluidic Device Detects a Role in Mechanosensitivity for AWCON Neuron in Caenorhabditis elegans

AWC olfactory neurons are fundamental for chemotaxis toward volatile attractants in Caenorhabditis elegans. Here, it is shown that AWC^ON responds not only to chemicals but also to mechanical stimuli caused by fluid flow changes in a microfluidic device. The dynamics of calcium events are correlated with the stimulus amplitude. It is further shown that the mechanosensitivity of AWC^ON neurons has an intrinsic nature rather than a synaptic origin, and the calcium transient response is mediated by TAX-4 cGMP-gated cation channel, suggesting the involvement of one or more "odorant" receptors in AWC^ON mechano-transduction. In many cases, the responses show plateau properties resembling bistable calcium dynamics where neurons can switch from one stable state to the other. To investigate the unprecedentedly observed mechanosensitivity of AWC^ON neurons, a novel microfluidic device is designed to minimize the fluid shear flow in the arena hosting the nematodes. Animals in this device show reduced neuronal activation of AWC^ON neurons. The results observed indicate that the tangential component of the mechanical stress is the main contributor to the mechanosensitivity of AWC^ON . Furthermore, the microfluidic platform, integrating shearless perfusion and calcium imaging, provides a novel and more controlled solution for in vivo analysis both in micro-organisms and cultured cells.