Mechanisms Underlying the Recruitment of Inhibitory Interneurons in Fictive Swimming in Developing Xenopus laevis Tadpoles.

Developing spinal circuits generate patterned motor outputs while many neurons with high membrane resistances are still maturing. In the spinal cord of hatchling frog tadpoles of unknown sex, we found that the firing reliability in swimming of inhibitory interneurons with commissural and ipsilateral ascending axons was negatively correlated with their cellular membrane resistance. Further analyses showed that neurons with higher resistances had outward rectifying properties, low firing thresholds, and little delay in firing evoked by current injections. Input synaptic currents these neurons received during swimming, either compound, unitary current amplitudes, or unitary synaptic current numbers, were scaled with their membrane resistances, but their own synaptic outputs were correlated with membrane resistances of their postsynaptic partners. Analyses of neuronal dendritic and axonal lengths and their activities in swimming and cellular input resistances did not reveal a clear correlation pattern. Incorporating these electrical and synaptic properties into a computer swimming model produced robust swimming rhythms, whereas randomizing input synaptic strengths led to the breakdown of swimming rhythms, coupled with less synchronized spiking in the inhibitory interneurons. We conclude that the recruitment of these developing interneurons in swimming can be predicted by cellular input resistances, but the order is opposite to the motor-strength-based recruitment scheme depicted by Henneman's size principle. This form of recruitment/integration order in development before the emergence of refined motor control is progressive potentially with neuronal acquisition of mature electrical and synaptic properties, among which the scaling of input synaptic strengths with cellular input resistance plays a critical role.SIGNIFICANCE STATEMENT The mechanisms on how interneurons are recruited to participate in circuit function in developing neuronal systems are rarely investigated. In 2-d-old frog tadpole spinal cord, we found the recruitment of inhibitory interneurons in swimming is inversely correlated with cellular input resistances, opposite to the motor-strength-based recruitment order depicted by Henneman's size principle. Further analyses showed the amplitude of synaptic inputs that neurons received during swimming was inversely correlated with cellular input resistances. Randomizing/reversing the relation between input synaptic strengths and membrane resistances in modeling broke down swimming rhythms. Therefore, the recruitment or integration of these interneurons is conditional on the acquisition of several electrical and synaptic properties including the scaling of input synaptic strengths with cellular input resistances.

From decision to action: Detailed modelling of frog tadpoles reveals neuronal mechanisms of decision-making and reproduces unpredictable swimming movements in response to sensory signals.

How does the brain process sensory stimuli, and decide whether to initiate locomotor behaviour? To investigate this question we develop two whole body computer models of a tadpole. The "Central Nervous System" (CNS) model uses evidence from whole-cell recording to define 2300 neurons in 12 classes to study how sensory signals from the skin initiate and stop swimming. In response to skin stimulation, it generates realistic sensory pathway spiking and shows how hindbrain sensory memory populations on each side can compete to initiate reticulospinal neuron firing and start swimming. The 3-D "Virtual Tadpole" (VT) biomechanical model with realistic muscle innervation, body flexion, body-water interaction, and movement is then used to evaluate if motor nerve outputs from the CNS model can produce swimming-like movements in a volume of "water". We find that the whole tadpole VT model generates reliable and realistic swimming. Combining these two models opens new perspectives for experiments.

A computational model of familiarity detection for natural pictures, abstract images, and random patterns: Combination of deep learning and anti-Hebbian training.

We present a neural network model for familiarity recognition of different types of images in the perirhinal cortex (the FaRe model). The model is designed as a two-stage system. At the first stage, the parameters of an image are extracted by a pretrained deep learning convolutional neural network. At the second stage, a two-layer feed forward neural network with anti-Hebbian learning is used to make the decision about the familiarity of the image. FaRe model simulations demonstrate high capacity of familiarity recognition memory for natural pictures and low capacity for both abstract images and random patterns. These findings are in agreement with psychological experiments.

A simple method defines 3D morphology and axon projections of filled neurons in a small CNS volume: Steps toward understanding functional network circuitry.

BACKGROUND: Fundamental to understanding neuronal network function is defining neuron morphology, location, properties, and synaptic connectivity in the nervous system. A significant challenge is to reconstruct individual neuron morphology and connections at a whole CNS scale and bring together functional and anatomical data to understand the whole network. NEW METHOD: We used a PC controlled micropositioner to hold a fixed whole mount of Xenopus tadpole CNS and replace the stage on a standard microscope. This allowed direct recording in 3D coordinates of features and axon projections of one or two neurons dye-filled during whole-cell recording to study synaptic connections. Neuron reconstructions were normalised relative to the ventral longitudinal axis of the nervous system. Coordinate data were stored as simple text files. RESULTS: Reconstructions were at 1 µm resolution, capturing axon lengths in mm. The output files were converted to SWC format and visualised in 3D reconstruction software NeuRomantic. Coordinate data are tractable, allowing correction for histological artefacts. Through normalisation across multiple specimens we could infer features of network connectivity of mapped neurons of different types. COMPARISON WITH EXISTING METHODS: Unlike other methods using fluorescent markers and utilising large-scale imaging, our method allows direct acquisition of 3D data on neurons whose properties and synaptic connections have been studied using whole-cell recording. CONCLUSIONS: This method can be used to reconstruct neuron 3D morphology and follow axon projections in the CNS. After normalisation to a common CNS framework, inferences on network connectivity at a whole nervous system scale contribute to network modelling to understand CNS function.

A Population Model of Deep Brain Stimulation in Movement Disorders From Circuits to Cells.

For more than 30 years, deep brain stimulation (DBS) has been used to target the symptoms of a number of neurological disorders and in particular movement disorders such as Parkinson's disease (PD) and essential tremor (ET). It is known that the loss of dopaminergic neurons in the substantia nigra leads to PD, while the exact impact of this on the brain dynamics is not fully understood, the presence of beta-band oscillatory activity is thought to be pathological. The cause of ET, however, remains uncertain, however pathological oscillations in the thalamocortical-cerebellar network have been linked to tremor. Both of these movement disorders are treated with DBS, which entails the surgical implantation of electrodes into a patient's brain. While DBS leads to an improvement in symptoms for many patients, the mechanisms underlying this improvement is not clearly understood, and computational modeling has been used extensively to improve this. Many of the models used to study DBS and its effect on the human brain have mainly utilized single neuron and single axon biophysical models. We have previously shown in separate models however, that the use of population models can shed much light on the mechanisms of the underlying pathological neural activity in PD and ET in turn, and on the mechanisms underlying DBS. Together, this work suggested that the dynamics of the cerebellar-basal ganglia thalamocortical network support oscillations at frequency range relevant to movement disorders. Here, we propose a new combined model of this network and present new results that demonstrate that both Parkinsonian oscillations in the beta band and oscillations in the tremor frequency range arise from the dynamics of such a network. We find regions in the parameter space demonstrating the different dynamics and go on to examine the transition from one oscillatory regime to another as well as the impact of DBS on these different types of pathological activity. This work will allow us to better understand the changes in brain activity induced by DBS, and allow us to optimize this clinical therapy, particularly in terms of target selection and parameter setting.

The decision to move: response times, neuronal circuits and sensory memory in a simple vertebrate.

All animals use sensory systems to monitor external events and have to decide whether to move. Response times are long and variable compared to reflexes, and fast escape movements. The complexity of adult vertebrate brains makes it difficult to trace the neuronal circuits underlying basic decisions to move. To simplify the problem, we investigate the nervous system and responses of hatchling frog tadpoles which swim when their skin is stimulated. Studying the neuron-by-neuron pathway from sensory to hindbrain neurons, where the decision to swim is made, has revealed two simple pathways generating excitation which sums to threshold in these neurons to initiate swimming. The direct pathway leads to short, and reliable delays like an escape response. The other includes a population of sensory processing neurons which extend firing to introduce noise and delay into responses. These neurons provide a brief, sensory memory of the stimulus, that allows tadpoles to integrate stimuli occurring within a second or so of each other. We relate these findings to other studies and conclude that sensory memory makes a fundamental contribution to simple decisions and is present in the brainstem of a basic vertebrate at a surprisingly early stage in development.

A simple decision to move in response to touch reveals basic sensory memory and mechanisms for variable response times.

KEY POINTS: Short-term working memory and decision-making are usually studied in the cerebral cortex; in many models of simple decision making, sensory signals build slowly and noisily to threshold to initiate a motor response after long, variable delays. When touched, hatchling frog tadpoles decide whether to swim; we define the long and variable delays to swimming and use whole-cell recordings to uncover the neurons and processes responsible. Firing in sensory and sensory pathway neurons is short latency, and too brief and invariant to explain these delays, while recordings from hindbrain reticulospinal neurons controlling swimming reveal a prolonged and variable build-up of synaptic excitation which can reach firing threshold and initiate swimming. We propose this excitation provides a sensory memory of the stimulus and may be generated by small reverberatory hindbrain networks. Our results uncover fundamental network mechanisms that allow animals to remember brief sensory stimuli and delay simple motor decisions. ABSTRACT: Many motor responses to sensory input, like locomotion or eye movements, are much slower than reflexes. Can simpler animals provide fundamental answers about the cellular mechanisms for motor decisions? Can we observe the 'accumulation' of excitation to threshold proposed to underlie decision making elsewhere? We explore how somatosensory touch stimulation leads to the decision to swim in hatchling Xenopus tadpoles. Delays measured to swimming in behaving and immobilised tadpoles are long and variable. Activity in their extensively studied sensory and sensory pathway neurons is too short-lived to explain these response delays. Instead, whole-cell recordings from the hindbrain reticulospinal neurons that drive swimming show that these receive prolonged, variable synaptic excitation lasting for nearly a second following a brief stimulus. They fire and initiate swimming when this excitation reaches threshold. Analysis of the summation of excitation requires us to propose extended firing in currently undefined presynaptic hindbrain neurons. Simple models show that a small excitatory recurrent-network inserted in the sensory pathway can mimic this process. We suggest that such a network may generate slow, variable summation of excitation to threshold. This excitation provides a simple memory of the sensory stimulus. It allows temporal and spatial integration of sensory inputs and explains the long, variable delays to swimming. The process resembles the 'accumulation' of excitation proposed for cortical circuits in mammals. We conclude that fundamental elements of sensory memory and decision making are present in the brainstem at a surprisingly early stage in development.

Bifurcations of Limit Cycles in a Reduced Model of the Xenopus Tadpole Central Pattern Generator.

We present the study of a minimal microcircuit controlling locomotion in two-day-old Xenopus tadpoles. During swimming, neurons in the spinal central pattern generator (CPG) generate anti-phase oscillations between left and right half-centres. Experimental recordings show that the same CPG neurons can also generate transient bouts of long-lasting in-phase oscillations between left-right centres. These synchronous episodes are rarely recorded and have no identified behavioural purpose. However, metamorphosing tadpoles require both anti-phase and in-phase oscillations for swimming locomotion. Previous models have shown the ability to generate biologically realistic patterns of synchrony and swimming oscillations in tadpoles, but a mathematical description of how these oscillations appear is still missing. We define a simplified model that incorporates the key operating principles of tadpole locomotion. The model generates the various outputs seen in experimental recordings, including swimming and synchrony. To study the model, we perform detailed one- and two-parameter bifurcation analysis. This reveals the critical boundaries that separate different dynamical regimes and demonstrates the existence of parameter regions of bi-stable swimming and synchrony. We show that swimming is stable in a significantly larger range of parameters, and can be initiated more robustly, than synchrony. Our results can explain the appearance of long-lasting synchrony bouts seen in experiments at the start of a swimming episode.

Structural and functional properties of a probabilistic model of neuronal connectivity in a simple locomotor network.

Although, in most animals, brain connectivity varies between individuals, behaviour is often similar across a species. What fundamental structural properties are shared across individual networks that define this behaviour? We describe a probabilistic model of connectivity in the hatchling Xenopus tadpole spinal cord which, when combined with a spiking model, reliably produces rhythmic activity corresponding to swimming. The probabilistic model allows calculation of structural characteristics that reflect common network properties, independent of individual network realisations. We use the structural characteristics to study examples of neuronal dynamics, in the complete network and various sub-networks, and this allows us to explain the basis for key experimental findings, and make predictions for experiments. We also study how structural and functional features differ between detailed anatomical connectomes and those generated by our new, simpler, model (meta-model).

Winner-take-all in a phase oscillator system with adaptation.

We consider a system of generalized phase oscillators with a central element and radial connections. In contrast to conventional phase oscillators of the Kuramoto type, the dynamic variables in our system include not only the phase of each oscillator but also the natural frequency of the central oscillator, and the connection strengths from the peripheral oscillators to the central oscillator. With appropriate parameter values the system demonstrates winner-take-all behavior in terms of the competition between peripheral oscillators for the synchronization with the central oscillator. Conditions for the winner-take-all regime are derived for stationary and non-stationary types of system dynamics. Bifurcation analysis of the transition from stationary to non-stationary winner-take-all dynamics is presented. A new bifurcation type called a Saddle Node on Invariant Torus (SNIT) bifurcation was observed and is described in detail. Computer simulations of the system allow an optimal choice of parameters for winner-take-all implementation.

Studying the role of axon fasciculation during development in a computational model of the Xenopus tadpole spinal cord.

During nervous system development growing axons can interact with each other, for example by adhering together in order to produce bundles (fasciculation). How does such axon-axon interaction affect the resulting axonal trajectories, and what are the possible benefits of this process in terms of network function? In this paper we study these questions by adapting an existing computational model of the development of neurons in the Xenopus tadpole spinal cord to include interactions between axons. We demonstrate that even relatively weak attraction causes bundles to appear, while if axons weakly repulse each other their trajectories diverge such that they fill the available space. We show how fasciculation can help to ensure axons grow in the correct location for proper network formation when normal growth barriers contain gaps, and use a functional spiking model to show that fasciculation allows the network to generate reliable swimming behaviour even when overall synapse counts are artificially lowered. Although we study fasciculation in one particular organism, our approach to modelling axon growth is general and can be widely applied to study other nervous systems.

To swim or not to swim: A population-level model of Xenopus tadpole decision making and locomotor behaviour.

We present a detailed computational model of interacting neuronal populations that mimic the hatchling Xenopus tadpole nervous system. The model includes four sensory pathways, integrators of sensory information, and a central pattern generator (CPG) network. Sensory pathways of different modalities receive inputs from an "environment"; these inputs are then processed and integrated to select the most appropriate locomotor action. The CPG populations execute the selected action, generating output in motor neuron populations. Thus, the model describes a detailed and biologically plausible chain of information processing from external signals to sensors, sensory pathways, integration and decision-making, action selection and execution and finally, generation of appropriate motor activity and behaviour. We show how the model produces appropriate behaviours in response to a selected scenario, which consists of a sequence of "environmental" signals. These behaviours might be relatively complex due to noisy sensory pathways and the possibility of spontaneous actions.

Advanced correlation grid: Analysis and visualisation of functional connectivity among multiple spike trains.

BACKGROUND: This study analyses multiple spike trains (MST) data, defines its functional connectivity and subsequently visualises an accurate diagram of connections. This is a challenging problem. For example, it is difficult to distinguish the common input and the direct functional connection of two spike trains. NEW METHOD: The new method presented in this paper is based on the traditional pairwise cross-correlation function (CCF) and a new combination of statistical techniques. First, the CCF is used to create the Advanced Correlation Grid (ACG) correlation where both the significant peak of the CCF and the corresponding time delay are used for detailed analysis of connectivity. Second, these two features of functional connectivity are used to classify connections. Finally, the visualization technique is used to represent the topology of functional connections. RESULTS: Examples are presented in the paper to demonstrate the new Advanced Correlation Grid method and to show how it enables discrimination between (i) influence from one spike train to another through an intermediate spike train and (ii) influence from one common spike train to another pair of analysed spike trains. COMPARISON WITH EXISTING METHODS: The ACG method enables scientists to automatically distinguish between direct connections from spurious connections such as common source connection and indirect connection whereas existing methods require in-depth analysis to identify such connections. CONCLUSIONS: The ACG is a new and effective method for studying functional connectivity of multiple spike trains. This method can identify accurately all the direct connections and can distinguish common source and indirect connections automatically.

Reaction times in visual search can be explained by a simple model of neural synchronization.

We present an oscillatory neural network model that can account for reaction times in visual search experiments. The model consists of a central oscillator that represents the central executive of the attention system and a number of peripheral oscillators that represent objects in the display. The oscillators are described as generalized Kuramoto type oscillators with adapted parameters. An object is considered as being included in the focus of attention if the oscillator associated with this object is in-phase with the central oscillator. The probability for an object to be included in the focus of attention is determined by its saliency that is described in formal terms as the strength of the connection from the peripheral oscillator to the central oscillator. By computer simulations it is shown that the model can reproduce reaction times in visual search tasks of various complexities. The dependence of the reaction time on the number of items in the display is represented by linear functions of different steepness which is in agreement with biological evidence.

A Network Model of Local Field Potential Activity in Essential Tremor and the Impact of Deep Brain Stimulation.

Essential tremor (ET), a movement disorder characterised by an uncontrollable shaking of the affected body part, is often professed to be the most common movement disorder, affecting up to one percent of adults over 40 years of age. The precise cause of ET is unknown, however pathological oscillations of a network of a number of brain regions are implicated in leading to the disorder. Deep brain stimulation (DBS) is a clinical therapy used to alleviate the symptoms of a number of movement disorders. DBS involves the surgical implantation of electrodes into specific nuclei in the brain. For ET the targeted region is the ventralis intermedius (Vim) nucleus of the thalamus. Though DBS is effective for treating ET, the mechanism through which the therapeutic effect is obtained is not understood. To elucidate the mechanism underlying the pathological network activity and the effect of DBS on such activity, we take a computational modelling approach combined with electrophysiological data. The pathological brain activity was recorded intra-operatively via implanted DBS electrodes, whilst simultaneously recording muscle activity of the affected limbs. We modelled the network hypothesised to underlie ET using the Wilson-Cowan approach. The modelled network exhibited oscillatory behaviour within the tremor frequency range, as did our electrophysiological data. By applying a DBS-like input we suppressed these oscillations. This study shows that the dynamics of the ET network support oscillations at the tremor frequency and the application of a DBS-like input disrupts this activity, which could be one mechanism underlying the therapeutic benefit.

Bifurcation study of phase oscillator systems with attractive and repulsive interaction.

We study a model of globally coupled phase oscillators that contains two groups of oscillators with positive (synchronizing) and negative (desynchronizing) incoming connections for the first and second groups, respectively. This model was previously studied by Hong and Strogatz (the Hong-Strogatz model) in the case of a large number of oscillators. We consider a generalized Hong-Strogatz model with a constant phase shift in coupling. Our approach is based on the study of invariant manifolds and bifurcation analysis of the system. In the case of zero phase shift, various invariant manifolds are analytically described and a new dynamical mode is found. In the case of a nonzero phase shift we obtained a set of bifurcation diagrams for various systems with three or four oscillators. It is shown that in these cases system dynamics can be complex enough and include multistability and chaotic oscillations.

The generation of antiphase oscillations and synchrony by a rebound-based vertebrate central pattern generator.

Many neural circuits are capable of generating multiple stereotyped outputs after different sensory inputs or neuromodulation. We have previously identified the central pattern generator (CPG) for Xenopus tadpole swimming that involves antiphase oscillations of activity between the left and right sides. Here we analyze the cellular basis for spontaneous left-right motor synchrony characterized by simultaneous bursting on both sides at twice the swimming frequency. Spontaneous synchrony bouts are rare in most tadpoles, and they instantly emerge from and switch back to swimming, most frequently within the first second after skin stimulation. Analyses show that only neurons that are active during swimming fire action potentials in synchrony, suggesting both output patterns derive from the same neural circuit. The firing of excitatory descending interneurons (dINs) leads that of other types of neurons in synchrony as it does in swimming. During synchrony, the time window between phasic excitation and inhibition is 7.9 ± 1 ms, shorter than that in swimming (41 ± 2.3 ms). The occasional, extra midcycle firing of dINs during swimming may initiate synchrony, and mismatches of timing in the left and right activity can switch synchrony back to swimming. Computer modeling supports these findings by showing that the same neural network, in which reciprocal inhibition mediates rebound firing, can generate both swimming and synchrony without circuit reconfiguration. Modeling also shows that lengthening the time window between phasic excitation and inhibition by increasing dIN synaptic/conduction delay can improve the stability of synchrony.

A developmental approach to predicting neuronal connectivity from small biological datasets: a gradient-based neuron growth model.

Relating structure and function of neuronal circuits is a challenging problem. It requires demonstrating how dynamical patterns of spiking activity lead to functions like cognitive behaviour and identifying the neurons and connections that lead to appropriate activity of a circuit. We apply a "developmental approach" to define the connectome of a simple nervous system, where connections between neurons are not prescribed but appear as a result of neuron growth. A gradient based mathematical model of two-dimensional axon growth from rows of undifferentiated neurons is derived for the different types of neurons in the brainstem and spinal cord of young tadpoles of the frog Xenopus. Model parameters define a two-dimensional CNS growth environment with three gradient cues and the specific responsiveness of the axons of each neuron type to these cues. The model is described by a nonlinear system of three difference equations; it includes a random variable, and takes specific neuron characteristics into account. Anatomical measurements are first used to position cell bodies in rows and define axon origins. Then a generalization procedure allows information on the axons of individual neurons from small anatomical datasets to be used to generate larger artificial datasets. To specify parameters in the axon growth model we use a stochastic optimization procedure, derive a cost function and find the optimal parameters for each type of neuron. Our biologically realistic model of axon growth starts from axon outgrowth from the cell body and generates multiple axons for each different neuron type with statistical properties matching those of real axons. We illustrate how the axon growth model works for neurons with axons which grow to the same and the opposite side of the CNS. We then show how, by adding a simple specification for dendrite morphology, our model "developmental approach" allows us to generate biologically-realistic connectomes.

Can simple rules control development of a pioneer vertebrate neuronal network generating behavior?

How do the pioneer networks in the axial core of the vertebrate nervous system first develop? Fundamental to understanding any full-scale neuronal network is knowledge of the constituent neurons, their properties, synaptic interconnections, and normal activity. Our novel strategy uses basic developmental rules to generate model networks that retain individual neuron and synapse resolution and are capable of reproducing correct, whole animal responses. We apply our developmental strategy to young Xenopus tadpoles, whose brainstem and spinal cord share a core vertebrate plan, but at a tractable complexity. Following detailed anatomical and physiological measurements to complete a descriptive library of each type of spinal neuron, we build models of their axon growth controlled by simple chemical gradients and physical barriers. By adding dendrites and allowing probabilistic formation of synaptic connections, we reconstruct network connectivity among up to 2000 neurons. When the resulting "network" is populated by model neurons and synapses, with properties based on physiology, it can respond to sensory stimulation by mimicking tadpole swimming behavior. This functioning model represents the most complete reconstruction of a vertebrate neuronal network that can reproduce the complex, rhythmic behavior of a whole animal. The findings validate our novel developmental strategy for generating realistic networks with individual neuron- and synapse-level resolution. We use it to demonstrate how early functional neuronal connectivity and behavior may in life result from simple developmental "rules," which lay out a scaffold for the vertebrate CNS without specific neuron-to-neuron recognition.

The emergence of two anti-phase oscillatory neural populations in a computational model of the Parkinsonian globus pallidus.

Experiments in rodent models of Parkinson's disease have demonstrated a prominent increase of oscillatory firing patterns in neurons within the Parkinsonian globus pallidus (GP) which may underlie some of the motor symptoms of the disease. There are two main pathways from the cortex to GP: via the striatum and via the subthalamic nucleus (STN), but it is not known how these inputs sculpt the pathological pallidal firing patterns. To study this we developed a novel neural network model of conductance-based spiking pallidal neurons with cortex-modulated input from STN neurons. Our results support the hypothesis that entrainment occurs primarily via the subthalamic pathway. We find that as a result of the interplay between excitatory input from the STN and mutual inhibitory coupling between GP neurons, a homogeneous population of GP neurons demonstrates a self-organizing dynamical behavior where two groups of neurons emerge: one spiking in-phase with the cortical rhythm and the other in anti-phase. This finding mirrors what is seen in recordings from the GP of rodents that have had Parkinsonism induced via brain lesions. Our model also includes downregulation of Hyperpolarization-activated Cyclic Nucleotide-gated (HCN) channels in response to burst firing of GP neurons, since this has been suggested as a possible mechanism for the emergence of Parkinsonian activity. We found that the downregulation of HCN channels provides even better correspondence with experimental data but that it is not essential in order for the two groups of oscillatory neurons to appear. We discuss how the influence of inhibitory striatal input will strengthen our results.