A computational model for motor pattern switching between taste-induced ingestion and rejection oromotor behaviors

Impact of precisely-timed inhibition of gustatory cortex on taste behavior depends on single-trial ensemble dynamics

Sensation and action are necessarily coupled during stimulus perception - while tasting, for instance, perception happens while an animal decides to expel or swallow the substance in the mouth (the former via a behavior known as 'gaping'). Taste responses in the rodent gustatory cortex (GC) span this sensorimotor divide, progressing through firing-rate epochs that culminate in the emergence of action-related firing. Population analyses reveal this emergence to be a sudden, coherent and variably-timed ensemble transition that reliably precedes gaping onset by 0.2-0.3s. Here, we tested whether this transition drives gaping, by delivering 0.5s GC perturbations in tasting trials. Perturbations significantly delayed gaping, but only when they preceded the action-related transition - thus, the same perturbation impacted behavior or not, depending on the transition latency in that particular trial. Our results suggest a distributed attractor network model of taste processing, and a dynamical role for cortex in driving motor behavior.

Propensity for Bistability of Bursting and Silence in the Leech Heart Interneuron

The coexistence of neuronal activity regimes has been reported under normal and pathological conditions. Such multistability could enhance the flexibility of the nervous system and has many implications for motor control, memory, and decision making. Multistability is commonly promoted by neuromodulation targeting specific membrane ionic currents. Here, we investigated how modulation of different ionic currents could affect the neuronal propensity for bistability. We considered a leech heart interneuron model. It exhibits bistability of bursting and silence in a narrow range of the leak current parameters, conductance (g_leak) and reversal potential (E_leak). We assessed the propensity for bistability of the model by using bifurcation diagrams. On the diagram (g_leak, E_leak), we mapped bursting and silent regimes. For the canonical value of E_leak we determined the range of g_leak which supported the bistability. We use this range as an index of propensity for bistability. We investigated how this index was affected by alterations of ionic currents. We systematically changed their conductances, one at a time, and built corresponding bifurcation diagrams in parameter planes of the maximal conductance of a given current and the leak conductance. We found that conductance of only one current substantially affected the index of propensity; the increase of the maximal conductance of the hyperpolarization-activated cationic current increased the propensity index. The second conductance with the strongest effect was the conductance of the low-threshold fast Ca²⁺ current; its reduction increased the propensity index although the effect was about two times smaller in magnitude. Analyzing the model with both changes applied simultaneously, we found that the diagram (g_leak, E_leak) showed a progressively expanded area of bistability of bursting and silence.

Sensory Cortical Activity Is Related to the Selection of a Rhythmic Motor Action Pattern

Rats produce robust, highly distinctive orofacial rhythms in response to taste stimuli-responses that aid in the consumption of palatable tastes and the ejection of aversive tastes, and that are sourced in a multifunctional brainstem central pattern generator. Several pieces of indirect evidence suggest that primary gustatory cortex (GC) may be a part of a distributed forebrain circuit involved in the selection of particular consumption-related rhythms, although not in the production of individual mouth movements per se. Here, we performed a series of tests of this hypothesis. We first examined the temporal relationship between GC activity and orofacial behaviors by performing paired single-neuron and electromyographic recordings in awake rats. Using a trial-by-trial analysis, we found that a subset of GC neurons shows a burst of activity beginning before the transition between nondistinct and taste-specific (i.e., consumption-related) orofacial rhythms. We further showed that shifting the latency of consumption-related behavior by selective cueing has an analogous impact on the timing of GC activity. Finally, we showed the complementary result, demonstrating that optogenetic perturbation of GC activity has a modest but significant impact on the probability that a specific rhythm will be produced in response to a strongly aversive taste. GC appears to be a part of a distributed circuit that governs the selection of taste-induced orofacial rhythms.

SIGNIFICANCE STATEMENT

In many well studied (typically invertebrate) sensorimotor systems, top-down modulation helps motor-control regions "select" movement patterns. Here, we provide evidence that gustatory cortex (GC) may be part of the forebrain circuit that performs this function in relation to oral behaviors ("gapes") whereby a substance in the mouth is rejected as unpalatable. We show that GC palatability coding is well timed to play this role, and that the latency of these codes changes as the latency of gaping shifts with learning. We go on to show that by silencing these neurons, we can change the likelihood of gaping. These data help to break down the sensory/motor divide by showing a role for sensory cortex in the selection of motor behavior.

High prevalence of multistability of rest states and bursting in a database of a model neuron

Flexibility in neuronal circuits has its roots in the dynamical richness of their neurons. Depending on their membrane properties single neurons can produce a plethora of activity regimes including silence, spiking and bursting. What is less appreciated is that these regimes can coexist with each other so that a transient stimulus can cause persistent change in the activity of a given neuron. Such multistability of the neuronal dynamics has been shown in a variety of neurons under different modulatory conditions. It can play either a functional role or present a substrate for dynamical diseases. We considered a database of an isolated leech heart interneuron model that can display silent, tonic spiking and bursting regimes. We analyzed only the cases of endogenous bursters producing functional half-center oscillators (HCOs). Using a one parameter (the leak conductance ()) bifurcation analysis, we extended the database to include silent regimes (stationary states) and systematically classified cases for the coexistence of silent and bursting regimes. We showed that different cases could exhibit two stable depolarized stationary states and two hyperpolarized stationary states in addition to various spiking and bursting regimes. We analyzed all cases of endogenous bursters and found that 18% of the cases were multistable, exhibiting coexistences of stationary states and bursting. Moreover, 91% of the cases exhibited multistability in some range of . We also explored HCOs built of multistable neuron cases with coexisting stationary states and a bursting regime. In 96% of cases analyzed, the HCOs resumed normal alternating bursting after one of the neurons was reset to a stationary state, proving themselves robust against this perturbation.

It is often not appreciated that different activity regimes can coexist with each other in a given neuron so that a transient stimulus can cause a persistent change of activity. Such multistability of the neuronal dynamics has in fact been shown in a variety of neurons and can play either a functional role or present a substrate for neurological diseases. We explored the propensity for multistability in a database of a leech heart interneuron model, testing each case (parameter set) in a database for multistability. We found a large proportion of multistable cases, especially the coexistence of silent and bursting regimes. This was a surprising result, since these cells pace the heartbeat of the leech, and the coexistence of silence and bursting could disrupt the functional pattern, threatening the viability of the leech. Analysis of networks of mutually inhibitory multistable neurons, however, showed robustness in maintaining functional activity, suggesting that the mutually inhibitory coupling can act as a protective mechanism against failures induced by multistability.

Bistability of bursting and silence regimes in a model of a leech heart interneuron

Bursting is one of the primary activity regimes of neurons. Our study is focused on determining a generic biophysical mechanism underlying the coexistence of the bursting and silent regimes observed in a neuron model. We show that the main ingredient for this mechanism is a saddle periodic orbit. The stable manifold of the orbit sets a threshold between the regimes of activity. Thus, the range of the controlling parameters, where the coexistence is observed, is limited by the bifurcations' values at which the saddle orbit appears and disappears. We show that it appears through the subcritical Andronov-Hopf bifurcation, where the equilibrium representing the silent regime loses stability, and disappears at the homoclinic bifurcation. Correspondingly, the bursting regime disappears in close proximity to the homoclinic bifurcation.

Six types of multistability in a neuronal model based on slow calcium current

BACKGROUND

Multistability of oscillatory and silent regimes is a ubiquitous phenomenon exhibited by excitable systems such as neurons and cardiac cells. Multistability can play functional roles in short-term memory and maintaining posture. It seems to pose an evolutionary advantage for neurons which are part of multifunctional Central Pattern Generators to possess multistability. The mechanisms supporting multistability of bursting regimes are not well understood or classified.

METHODOLOGY/PRINCIPAL FINDINGS

Our study is focused on determining the bio-physical mechanisms underlying different types of co-existence of the oscillatory and silent regimes observed in a neuronal model. We develop a low-dimensional model typifying the dynamics of a single leech heart interneuron. We carry out a bifurcation analysis of the model and show that it possesses six different types of multistability of dynamical regimes. These types are the co-existence of 1) bursting and silence, 2) tonic spiking and silence, 3) tonic spiking and subthreshold oscillations, 4) bursting and subthreshold oscillations, 5) bursting, subthreshold oscillations and silence, and 6) bursting and tonic spiking. These first five types of multistability occur due to the presence of a separating regime that is either a saddle periodic orbit or a saddle equilibrium. We found that the parameter range wherein multistability is observed is limited by the parameter values at which the separating regimes emerge and terminate.

CONCLUSIONS

We developed a neuronal model which exhibits a rich variety of different types of multistability. We described a novel mechanism supporting the bistability of bursting and silence. This neuronal model provides a unique opportunity to study the dynamics of networks with neurons possessing different types of multistability.

Modulation of inhibitory strength and kinetics facilitates regulation of persistent inward currents and motoneuron excitability following spinal cord injury

Spasticity is commonly observed after chronic spinal cord injury (SCI) and many other central nervous system disorders (e.g., multiple sclerosis, stroke). SCI-induced spasticity has been associated with motoneuron hyperexcitability partly due to enhanced activation of intrinsic persistent inward currents (PICs). Disrupted spinal inhibitory mechanisms also have been implicated. Altered inhibition can result from complex changes in the strength, kinetics, and reversal potential (E(Cl(-))) of γ-aminobutyric acid A (GABA(A)) and glycine receptor currents. Development of optimal therapeutic strategies requires an understanding of the impact of these interacting factors on motoneuron excitability. We employed computational methods to study the effects of conductance, kinetics, and E(Cl(-)) of a dendritic inhibition on PIC activation and motoneuron discharge. A two-compartment motoneuron with enhanced PICs characteristic of SCI and receiving recurrent inhibition from Renshaw cells was utilized in these simulations. This dendritic inhibition regulated PIC onset and offset and exerted its strongest effects at motoneuron recruitment and in the secondary range of the current-frequency relationship during PIC activation. Increasing inhibitory conductance compensated for moderate depolarizing shifts in E(Cl(-)) by limiting PIC activation and self-sustained firing. Furthermore, GABA(A) currents exerted greater control on PIC activation than glycinergic currents, an effect attributable to their slower kinetics. These results suggest that modulation of the strength and kinetics of GABA(A) currents could provide treatment strategies for uncontrollable spasms.

Intrinsic membrane properties of pre-oromotor neurons in the intermediate zone of the medullary reticular formation

Neurons in the lower brainstem that control consummatory behavior are widely distributed in the reticular formation (RF) of the pons and medulla. The intrinsic membrane properties of neurons within this distributed system shape complex excitatory and inhibitory inputs from both orosensory and central structures implicated in homeostatic control to produce coordinated oromotor patterns. The current study explored the intrinsic membrane properties of neurons in the intermediate subdivision of the medullary reticular formation (IRt). Neurons in the IRt receive input from the overlying (gustatory) nucleus of the solitary tract and project to the oromotor nuclei. Recent behavioral pharmacology studies as well as computational modeling suggest that inhibition in the IRt plays an important role in the transition from a taste-initiated oromotor pattern of ingestion to one of rejection. The present study explored the impact of hyperpolarization on membrane properties. In response to depolarization, neurons responded with either a tonic discharge, an irregular/burst pattern or were spike-adaptive. A hyperpolarizing pre-pulse modulated the excitability of most (82%) IRt neurons to subsequent depolarization. Instances of both increased (30%) and decreased (52%) excitability were observed. Currents induced by the hyperpolarization included an outward 4-aminopyridine (4-AP) sensitive K+ current that suppressed excitability and an inward cation current that increased excitability. These currents are also present in other subpopulations of RF neurons that influence the oromotor nuclei and we discuss how these currents could alter firing characteristics to impact pattern generation.

Taste-specific cell assemblies in a biologically informed model of the nucleus of the solitary tract

Although the cellular organization of many primary sensory nuclei has been well characterized, questions remain about the functional architecture of the first central relay for gustation, the rostral nucleus of the solitary tract (NTS). Here we used electrophysiological data recorded from single cells in the NTS to inform a network model of taste processing. Previous studies showed that electrical stimulation of the chorda tympani (CT) nerve initiates two types of inhibitory influences with different time courses in separate groups of NTS cells. Each type of inhibition targeted cells with distinct taste response properties. Further analyses of these data identified three NTS cell types differentiated by their latency of evoked response, time course of CT evoked inhibition, and degree of selectivity across taste qualities. Based on these results, we designed a model of the NTS consisting of discrete, reciprocally connected, stimulus-specific "cell" assemblies. Input to the network of integrate-and-fire model neurons was based on electrophysiological recordings from the CT nerve. Following successful simulation of paired-pulse CT stimulation, the network was tested for its ability to discriminate between two "taste" stimuli. Network dynamics of the model produced biologically plausible responses from each unit type and enhanced discrimination between taste qualities. We propose that an interactive network of taste quality specific cell assemblies, similar to our model, may account for the coherence in across-neuron patterns of NTS responses between similar tastants.

Cortical networks produce three distinct 7-12 Hz rhythms during single sensory responses in the awake rat

Cortical rhythms in the alpha/mu frequency range (7-12 Hz) have been variously related to "idling," anticipation, seizure, and short-term or working memory. This overabundance of interpretations suggests that sensory cortex may be able to produce more than one (and even more than two) distinct alpha/mu rhythms. Here we describe simultaneous local field potential and single-neuron recordings made from primary sensory (gustatory) cortex of awake rats and reveal three distinct 7-12 Hz de novo network rhythms within single sessions: an "early," taste-induced approximately 11 Hz rhythm, the first peak of which was a short-latency gustatory evoked potential; a "late," significantly lower-frequency (approximately 7 Hz) rhythm that replaced this first rhythm at approximately 750-850 ms after stimulus onset (consistently timed with a previously described shift in taste temporal codes); and a "spontaneous" spike-and-wave rhythm of intermediate peak frequency (approximately 9 Hz) that appeared late in the session, as part of a oft-described reduction in arousal/attention. These rhythms proved dissociable on many grounds: in addition to having different peak frequencies, amplitudes, and shapes and appearing at different time points (although often within single 3 s snippets of activity), the early and late rhythms proved to have completely uncorrelated session-to-session variability, and the spontaneous rhythm affected the early rhythm only (having no impact on the late rhythm). Analysis of spike-to-wave coupling suggested that the early and late rhythms are a unified part of discriminative taste process: the identity of phase-coupled single-neuron ensembles differed from taste to taste, and coupling typically lasted across the change in frequency. These data reveal that even rhythms confined to a narrow frequency band may still have distinct properties.

Expression of Fos during sham sucrose intake in rats with central gustatory lesions

For humans and rodents, ingesting sucrose is rewarding. This experiment tested the prediction that the neural activity produced by sapid sucrose reaches reward systems via projections from the pons through the limbic system. Gastric cannulas drained ingested fluid before absorption. For 10 days, the rats alternated an hour of this sham ingestion between sucrose and water. On the final test day, half of them sham drank water and the other half 0.6 M sucrose. Thirty minutes later, the rats were killed and their brains immunohistochemically stained for Fos. The groups consisted of controls and rats with excitotoxic lesions in the gustatory thalamus (TTA), the medial (gustatory) parabrachial nucleus (PBN), or the lateral (visceral afferent) parabrachial nucleus. In controls, compared with water, sham ingesting sucrose produced significantly more Fos-positive neurons in the nucleus of the solitary tract, PBN, TTA, and gustatory cortex (GC). In the ventral forebrain, sucrose sham licking increased Fos in the bed nucleus of the stria terminalis, central nucleus of amygdala, and the shell of nucleus accumbens. Thalamic lesions blocked the sucrose effect in GC but not in the ventral forebrain. After lateral PBN lesions, the Fos distributions produced by distilled H(2)O or sucrose intake did not differ from controls. Bilateral medial PBN damage, however, eliminated the sucrose-induced Fos increase not only in the TTA and GC but also in the ventral forebrain. Thus ventral forebrain areas associated with affective responses appear to be activated directly by PBN gustatory neurons rather than via the thalamocortical taste system.

Local circuit input to the medullary reticular formation from the rostral nucleus of the solitary tract

The intermediate reticular formation (IRt) subjacent to the rostral (gustatory) nucleus of the solitary tract (rNST) receives projections from the rNST and appears essential to the expression of taste-elicited ingestion and rejection responses. We used whole cell patch-clamp recording and calcium imaging to characterize responses from an identified population of prehypoglossal neurons in the IRt to electrical stimulation of the rNST in a neonatal rat pup slice preparation. The calcium imaging studies indicated that IRt neurons could be activated by rNST stimulation and that many neurons were under tonic inhibition. Whole cell patch-clamp recording revealed mono- and polysynaptic projections from the rNST to identified prehypoglossal neurons. The projection was primarily excitatory and glutamatergic; however, there were some inhibitory GABAergic projections, and many neurons received excitatory and inhibitory inputs. There was also evidence of disinhibition. Overall, bath application of GABA(A) antagonists increased the amplitude of excitatory currents, and, in several neurons, stimulation of the rNST systematically decreased inhibitory currents. We have hypothesized that the transition from licks to gapes by natural stimuli, such as quinine monohydrochloride, could occur via such disinhibition. We present an updated dynamic model that summarizes the complex synaptic interface between the rNST and the IRt and demonstrates how inhibition could contribute to the transition from ingestion to rejection.