Search

The recent advancements in large-scale activity imaging of neuronal ensembles offer valuable opportunities to comprehend the process involved in generating brain activity patterns and understanding how information is transmitted between neurons or neuronal ensembles. However, existing methodologies for extracting the underlying properties that generate overall dynamics are still limited. In this study, we applied previously unexplored methodologies to analyze time-lapse 3D imaging (4D imaging) data of head neurons of the nematode Caenorhabditis elegans. By combining time-delay embedding with the independent component analysis, we successfully decomposed whole-brain activities into a small number of component dynamics. Through the integration of results from multiple samples, we extracted common dynamics from neuronal activities that exhibit apparent divergence across different animals. Notably, while several components show common cooperativity across samples, some component pairs exhibited distinct relationships between individual samples. We further developed time series prediction models of synaptic communications. By combining dimension reduction using the general framework, gradient kernel dimension reduction, and probabilistic modeling, the overall relationships of neural activities were incorporated. By this approach, the stochastic but coordinated dynamics were reproduced in the simulated whole-brain neural network. We found that noise in the nervous system is crucial for generating realistic whole-brain dynamics. Furthermore, by evaluating synaptic interaction properties in the models, strong interactions within the core neural circuit, variable sensory transmission and importance of gap junctions were inferred. Virtual optogenetics can be also performed using the model. These analyses provide a solid foundation for understanding information flow in real neural networks.

Abstract Background In the field of neuroscience, neural modules and circuits that control biological functions have been found throughout entire neural networks. Correlations in neural activity can be used to identify such neural modules. Recent technological advances enable us to measure whole-brain neural activity with single-cell resolution in several species including $$Caenorhabditis\ elegans$$ C a e n o r h a b d i t i s e l e g a n s . Because current neural activity data in C. elegans contain many missing data points, it is necessary to merge results from as many animals as possible to obtain more reliable functional modules. Results In this work, we developed a new time-series clustering method, WormTensor, to identify functional modules using whole-brain activity data from C. elegans. WormTensor uses a distance measure, modified shape-based distance to account for the lags and the mutual inhibition of cell–cell interactions and applies the tensor decomposition algorithm multi-view clustering based on matrix integration using the higher orthogonal iteration of tensors (HOOI) algorithm (MC-MI-HOOI), which can estimate both the weight to account for the reliability of data from each animal and the clusters that are common across animals. Conclusion We applied the method to 24 individual C. elegans and successfully found some known functional modules. Compared with a widely used consensus clustering method to aggregate multiple clustering results, WormTensor showed higher silhouette coefficients. Our simulation also showed that WormTensor is robust to contamination from noisy data. WormTensor is freely available as an R/CRAN package https://cran.r-project.org/web/packages/WormTensor .

Abstract Background Annotation of cell identity is an essential process in neuroscience that allows comparison of cells, including that of neural activities across different animals. In Caenorhabditis elegans, although unique identities have been assigned to all neurons, the number of annotatable neurons in an intact animal has been limited due to the lack of quantitative information on the location and identity of neurons. Results Here, we present a dataset that facilitates the annotation of neuronal identities, and demonstrate its application in a comprehensive analysis of whole-brain imaging. We systematically identified neurons in the head region of 311 adult worms using 35 cell-specific promoters and created a dataset of the expression patterns and the positions of the neurons. We found large positional variations that illustrated the difficulty of the annotation task. We investigated multiple combinations of cell-specific promoters driving distinct fluorescence and generated optimal strains for the annotation of most head neurons in an animal. We also developed an automatic annotation method with human interaction functionality that facilitates annotations needed for whole-brain imaging. Conclusion Our neuron ID dataset and optimal fluorescent strains enable the annotation of most neurons in the head region of adult C. elegans, both in full-automated fashion and a semi-automated version that includes human interaction functionalities. Our method can potentially be applied to model species used in research other than C. elegans, where the number of available cell-type-specific promoters and their variety will be an important consideration.

Cellular homeostasis is regulated by signals through multiple molecular networks that include protein phosphorylation and metabolites. However, where and when the signal flows through a network and regulates homeostasis has not been explored. We have developed a reconstruction method for the signal flow based on time-course phosphoproteome and metabolome data, using multiple databases, and have applied it to acute action of insulin, an important hormone for metabolic homeostasis. An insulin signal flows through a network, through signaling pathways that involve 13 protein kinases, 26 phosphorylated metabolic enzymes, and 35 allosteric effectors, resulting in quantitative changes in 44 metabolites. Analysis of the network reveals that insulin induces phosphorylation and activation of liver-type phosphofructokinase 1, thereby controlling a key reaction in glycolysis. We thus provide a versatile method of reconstruction of signal flow through the network using phosphoproteome and metabolome data.

Abstract Insulin governs systemic glucose metabolism, including glycolysis, gluconeogenesis and glycogenesis, through temporal change and absolute concentration. However, how insulin‐signalling pathway selectively regulates glycolysis, gluconeogenesis and glycogenesis remains to be elucidated. To address this issue, we experimentally measured metabolites in glucose metabolism in response to insulin. Step stimulation of insulin induced transient response of glycolysis and glycogenesis, and sustained response of gluconeogenesis and extracellular glucose concentration (GLCex). Based on the experimental results, we constructed a simple computational model that characterises response of insulin‐signalling‐dependent glucose metabolism. The model revealed that the network motifs of glycolysis and glycogenesis pathways constitute a feedforward (FF) with substrate depletion and incoherent feedforward loop (iFFL), respectively, enabling glycolysis and glycogenesis responsive to temporal changes of insulin rather than its absolute concentration. In contrast, the network motifs of gluconeogenesis pathway constituted a FF inhibition, enabling gluconeogenesis responsive to absolute concentration of insulin regardless of its temporal patterns. GLCexwas regulated by gluconeogenesis and glycolysis. These results demonstrate the selective control mechanism of glucose metabolism by temporal patterns of insulin.

To measure the activity of neurons using whole-brain activity imaging, precise detection of each neuron or its nucleus is required. In the head region of the nematode C. elegans, the neuronal cell bodies are distributed densely in three-dimensional (3D) space. However, no existing computational methods of image analysis can separate them with sufficient accuracy. Here we propose a highly accurate segmentation method based on the curvatures of the iso-intensity surfaces. To obtain accurate positions of nuclei, we also developed a new procedure for least squares fitting with a Gaussian mixture model. Combining these methods enables accurate detection of densely distributed cell nuclei in a 3D space. The proposed method was implemented as a graphical user interface program that allows visualization and correction of the results of automatic detection. Additionally, the proposed method was applied to time-lapse 3D calcium imaging data, and most of the nuclei in the images were successfully tracked and measured.

Homeostatic control of blood glucose is regulated by a complex feedback loop between glucose and insulin, of which failure leads to diabetes mellitus. However, physiological and pathological nature of the feedback loop is not fully understood. We made a mathematical model of the feedback loop between glucose and insulin using time course of blood glucose and insulin during consecutive hyperglycemic and hyperinsulinemic-euglycemic clamps in 113 subjects with variety of glucose tolerance including normal glucose tolerance (NGT), impaired glucose tolerance (IGT) and type 2 diabetes mellitus (T2DM). We analyzed the correlation of the parameters in the model with the progression of glucose intolerance and the conserved relationship between parameters. The model parameters of insulin sensitivity and insulin secretion significantly declined from NGT to IGT, and from IGT to T2DM, respectively, consistent with previous clinical observations. Importantly, insulin clearance, an insulin degradation rate, significantly declined from NGT, IGT to T2DM along the progression of glucose intolerance in the mathematical model. Insulin clearance was positively correlated with a product of insulin sensitivity and secretion assessed by the clamp analysis or determined with the mathematical model. Insulin clearance was correlated negatively with postprandial glucose at 2h after oral glucose tolerance test. We also inferred a square-law between the rate constant of insulin clearance and a product of rate constants of insulin sensitivity and secretion in the model, which is also conserved among NGT, IGT and T2DM subjects. Insulin clearance shows a conserved relationship with the capacity of glucose disposal among the NGT, IGT and T2DM subjects. The decrease of insulin clearance predicts the progression of glucose intolerance.

Animals navigate their environment by manipulating their movements and adjusting their trajectory which requires a sophisticated integration of sensory data with their current motor status. Here, we utilize the nematode Caenorhabditis elegans to explore the neural mechanisms of processing the sensory and motor information for navigation. We developed a microfluidic device which allows animals to freely move their heads while receiving temporal NaCl stimuli. We found that C. elegans regulates neck bending direction in response to temporal NaCl concentration changes in a way which is consistent with a C. elegans' navigational strategy which regulates traveling direction toward preferred NaCl concentrations. Our analysis also revealed that the activity of a neck motor neuron is significantly correlated with neck bending and activated by the decrease in NaCl concentration in a phase-dependent manner. By combining the analysis of behavioral and neural response to NaCl stimuli and optogenetic perturbation experiments, we revealed that NaCl decrease during ventral bending activates the neck motor neuron which counteracts ipsilateral bending. Simulations further suggest that this phase-dependent response of neck motor neurons can facilitate curving toward preferred salt concentrations. [ABSTRACT FROM AUTHOR]

Animals demonstrate flexible behaviors through associative learning based on their experiences. Deciphering the neural mechanisms for sensing and integrating multiple types of sensory information is critical for understanding such behavioral controls. The soil nematode Caenorhabditis elegans avoids salt concentrations it has previously experienced under starvation conditions. Here, we identify a pair of sensory neurons, the ASG neuron pair, which in cooperation with the ASER salt-sensing neuron generate starvation-dependent salt avoidance. Animals whose sensory input is restricted to only ASER failed to show learned avoidance due to inappropriately directed navigation behaviors. However, their navigation through a salt concentration gradient was improved by allowing sensory inputs to ASG in addition to ASER. Detailed behavioral analyses of these animals revealed that input from ASG neurons is required not only for controlling the frequency of initiating a set of sharp turns (called pirouettes) based on detected ambient salt concentrations but also adjusting the migration direction during pirouettes. Optogenetic activation of ASER by ChR2 induced turning behaviors in a salt concentration-dependent manner where presence of intact ASG was important for the starvation-dependent responses. Calcium imaging of the activity of ASG neurons in freely moving worms revealed that ASG is activated upon turning behavior. Our results indicate that ASG neurons cooperate with the ASER neuron to generate destination-directed reorientation in starvationassociated salt concentration avoidance. [ABSTRACT FROM AUTHOR]

Sensory receptor neurons match their dynamic range to ecologically relevant stimulus intensities. How this tuning is achieved is poorly understood in most receptors. The roundworm Caenorhabditis elegans avoids 21% O2 and hypoxia and prefers intermediate O2 concentrations. We show how this O2 preference is sculpted by the antagonistic action of a neuroglobin and an O2-binding soluble guanylate cyclase. These putative molecular O2 sensors confer a sigmoidal O2 response curve in the URX neurons that has highest slope between 15 and 19% O2 and approaches saturation when O2 reaches 21%. In the absence of the neuroglobin, the response curve is shifted to lower O2 values and approaches saturation at 14% O2. In behavioral terms, neuroglobin signaling broadens the O2 preference of Caenorhabditis elegans while maintaining avoidance of 21% O2. A computational model of aerotaxis suggests the relationship between GLB-5-modulated URX responses and reversal behavior is sufficient to broaden O2 preference. In summary, we show that a neuroglobin can shift neural information coding leading to altered behavior. Antagonistically acting molecular sensors may represent a common mechanism to sharpen tuning of sensory neurons. [ABSTRACT FROM AUTHOR]