Behavioral Forgetting of Olfactory Learning Is Mediated by Interneuron-Regulated Network Plasticity in Caenorhabditis elegans.

Forgetting is important for animals to manage acquired memories to enable adaptation to changing environments; however, the neural network in mechanisms of forgetting is not fully understood. To understand the mechanisms underlying forgetting, we examined olfactory adaptation, a form of associative learning, in Caenorhabditis elegans The forgetting of diacetyl olfactory adaptation in C. elegans is regulated by secreted signals from AWC sensory neurons via the TIR-1/JNK-1 pathway. These signals cause a decline of the sensory memory trace in AWA neurons, where diacetyl is mainly sensed. To further understand the neural network that regulates this forgetting, we investigated the function of interneurons downstream of AWA and AWC neurons. We found that a pair of interneurons, AIA, is indispensable for the proper regulation of behavioral forgetting of diacetyl olfactory adaptation. Loss or inactivation of AIA caused the impairment of the chemotaxis recovery after adaptation without causing severe chemotaxis defects in the naive animal. AWA Ca2+ imaging analyses suggested that loss or inactivation of AIA interneurons did not affect the decline of the sensory memory trace after the recovery. Furthermore, AIA responses to diacetyl were observed in naive animals and after the recovery, but not just after the conditioning, suggesting that AIA responses after the recovery are required for the chemotaxis to diacetyl. We propose that the functional neuronal circuit for attractive chemotaxis to diacetyl is changed temporally at the recovery phase so that AIA interneurons are required for chemotaxis, although AIAs are dispensable for attractive chemotaxis to diacetyl in naive animals.

Cohesive and anisotropic vascular endothelial cell motility driving angiogenic morphogenesis.

Vascular endothelial cells (ECs) in angiogenesis exhibit inhomogeneous collective migration called "cell mixing", in which cells change their relative positions by overtaking each other. However, how such complex EC dynamics lead to the formation of highly ordered branching structures remains largely unknown. To uncover hidden laws of integration driving angiogenic morphogenesis, we analyzed EC behaviors in an in vitro angiogenic sprouting assay using mouse aortic explants in combination with mathematical modeling. Time-lapse imaging of sprouts extended from EC sheets around tissue explants showed directional cohesive EC movements with frequent U-turns, which often coupled with tip cell overtaking. Imaging of isolated branches deprived of basal cell sheets revealed a requirement of a constant supply of immigrating cells for ECs to branch forward. Anisotropic attractive forces between neighboring cells passing each other were likely to underlie these EC motility patterns, as evidenced by an experimentally validated mathematical model. These results suggest that cohesive movements with anisotropic cell-to-cell interactions characterize the EC motility, which may drive branch elongation depending on a constant cell supply. The present findings provide novel insights into a cell motility-based understanding of angiogenic morphogenesis.

SPF-CellTracker: Tracking Multiple Cells with Strongly-Correlated Moves Using a Spatial Particle Filter.

Tracking many cells in time-lapse 3D image sequences is an important challenging task of bioimage informatics. Motivated by a study of brain-wide 4D imaging of neural activity in C. elegans, we present a new method of multi-cell tracking. Data types to which the method is applicable are characterized as follows: (i) cells are imaged as globular-like objects, (ii) it is difficult to distinguish cells on the basis of shape and size only, (iii) the number of imaged cells in the several-hundred range, (iv) movements of nearly-located cells are strongly correlated, and (v) cells do not divide. We developed a tracking software suite that we call SPF-CellTracker. Incorporating dependency on the cells' movements into the prediction model is the key for reducing the tracking errors: the cell switching and the coalescence of the tracked positions. We model the target cells' correlated movements as a Markov random field and we also derive a fast computation algorithm, which we call spatial particle filter. With the live-imaging data of the nuclei of C. elegans neurons in which approximately 120 nuclei of neurons were imaged, the proposed method demonstrated improved accuracy compared to the standard particle filter and the method developed by Tokunaga et al. (2014).

Accurate Automatic Detection of Densely Distributed Cell Nuclei in 3D Space.

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.

Automated detection and tracking of many cells by using 4D live-cell imaging data.

MOTIVATION: Automated fluorescence microscopes produce massive amounts of images observing cells, often in four dimensions of space and time. This study addresses two tasks of time-lapse imaging analyses; detection and tracking of the many imaged cells, and it is especially intended for 4D live-cell imaging of neuronal nuclei of Caenorhabditis elegans. The cells of interest appear as slightly deformed ellipsoidal forms. They are densely distributed, and move rapidly in a series of 3D images. Thus, existing tracking methods often fail because more than one tracker will follow the same target or a tracker transits from one to other of different targets during rapid moves. RESULTS: The present method begins by performing the kernel density estimation in order to convert each 3D image into a smooth, continuous function. The cell bodies in the image are assumed to lie in the regions near the multiple local maxima of the density function. The tasks of detecting and tracking the cells are then addressed with two hill-climbing algorithms. The positions of the trackers are initialized by applying the cell-detection method to an image in the first frame. The tracking method keeps attacking them to near the local maxima in each subsequent image. To prevent the tracker from following multiple cells, we use a Markov random field (MRF) to model the spatial and temporal covariation of the cells and to maximize the image forces and the MRF-induced constraint on the trackers. The tracking procedure is demonstrated with dynamic 3D images that each contain >100 neurons of C.elegans. AVAILABILITY: http://daweb.ism.ac.jp/yoshidalab/crest/ismb2014 SUPPLEMENTARY INFORMATION: Supplementary data are available at http://daweb.ism.ac.jp/yoshidalab/crest/ismb2014

Separation of stationary and non-stationary sources with a generalized eigenvalue problem.

Non-stationary effects are ubiquitous in real world data. In many settings, the observed signals are a mixture of underlying stationary and non-stationary sources that cannot be measured directly. For example, in EEG analysis, electrodes on the scalp record the activity from several sources located inside the brain, which one could only measure invasively. Discerning stationary and non-stationary contributions is an important step towards uncovering the mechanisms of the data generating system. To that end, in Stationary Subspace Analysis (SSA), the observed signal is modeled as a linear superposition of stationary and non-stationary sources, where the aim is to separate the two groups in the mixture. In this paper, we propose the first SSA algorithm that has a closed form solution. The novel method, Analytic SSA (ASSA), is more than 100 times faster than the state-of-the-art, numerically stable, and guaranteed to be optimal when the covariance between stationary and non-stationary sources is time-constant. In numerical simulations on wide range of settings, we show that our method yields superior results, even for signals with time-varying group-wise covariance. In an application to geophysical data analysis, ASSA extracts meaningful components that shed new light on the Pi 2 pulsations of the geomagnetic field.