Feedback Inhibition and Throughput Properties of an Integrate-and-Fire-or-Burst Network Model of Retinogeniculate Transmission

Computational modeling has played an important role in the dissection of the biophysical basis of rhythmic oscillations in thalamus that are associated with sleep and certain forms of epilepsy. In contrast, the dynamic filter properties of thalamic relay nuclei during states of arousal are not well understood. Here we present a modeling and simulation study of the throughput properties of the visually driven dorsal lateral geniculate nucleus (dLGN) in the presence of feedback inhibition from the perigeniculate nucleus (PGN). We employ thalamocortical (TC) and thalamic reticular (RE) versions of a minimal integrate-and-fire-or-burst type model and a one-dimensional, two-layered network architecture. Potassium leakage conductances control the neuromodulatory state of the network and eliminate rhythmic bursting in the presence of spontaneous input (i.e., wake up the network). The aroused dLGN/PGN network model is subsequently stimulated by spatially homogeneous spontaneous retinal input or spatio-temporally patterned input consistent with the activity of X-type retinal ganglion cells during full-field or drifting grating visual stimulation. The throughput properties of this visually-driven dLGN/PGN network model are characterized and quantified as a function of stimulus parameters such as contrast, temporal frequency, and spatial frequency. During low-frequency oscillatory full-field stimulation, feedback inhibition from RE neurons often leads to TC neuron burst responses, while at high frequency tonic responses dominate. Depending on the average rate of stimulation, contrast level, and temporal frequency of modulation, the TC and RE cell bursts may or may not be phase-locked to the visual stimulus. During drifting-grating stimulation, phase-locked bursts often occur for sufficiently high contrast so long as the spatial period of the grating is not small compared to the synaptic footprint length, i.e., the spatial scale of the network connectivity.     

A multi-compartment model for interneurons in the dorsal lateral geniculate nucleus

GABAergic interneurons (INs) in the dorsal lateral geniculate nucleus (dLGN) shape the information flow from retina to cortex, presumably by controlling the number of visually evoked spikes in geniculate thalamocortical (TC) neurons, and refining their receptive field. The INs exhibit a rich variety of firing patterns: Depolarizing current injections to the soma may induce tonic firing, periodic bursting or an initial burst followed by tonic spiking, sometimes with prominent spike-time adaptation. When released from hyperpolarization, some INs elicit rebound bursts, while others return more passively to the resting potential. A full mechanistic understanding that explains the function of the dLGN on the basis of neuronal morphology, physiology and circuitry is currently lacking. One way to approach such an understanding is by developing a detailed mathematical model of the involved cells and their interactions. Limitations of the previous models for the INs of the dLGN region prevent an accurate representation of the conceptual framework needed to understand the computational properties of this region. We here present a detailed compartmental model of INs using, for the first time, a morphological reconstruction and a set of active dendritic conductances constrained by experimental somatic recordings from INs under several different current-clamp conditions. The model makes a number of experimentally testable predictions about the role of specific mechanisms for the firing properties observed in these neurons. In addition to accounting for the significant features of all experimental traces, it quantitatively reproduces the experimental recordings of the action-potential- firing frequency as a function of injected current. We show how and why relative differences in conductance values, rather than differences in ion channel composition, could account for the distinct differences between the responses observed in two different neurons, suggesting that INs may be individually tuned to optimize network operation under different input conditions.     

Active dendritic conductances dynamically regulate GABA release from thalamic interneurons

Inhibitory interneurons in the dorsal lateral geniculate nucleus (dLGN) process visual information by precisely controlling spike timing and by refining the receptive fields of thalamocortical (TC) neurons. Previous studies indicate that dLGN interneurons inhibit TC neurons by releasing GABA from both axons and dendrites. However, the mechanisms controlling GABA release are poorly understood. Here, using simultaneous whole-cell recordings from interneurons and TC neurons and two-photon calcium imaging, we find that synchronous activation of multiple retinal ganglion cells (RGCs) triggers sodium spikes that propagate throughout interneuron axons and dendrites, and calcium spikes that invade dendrites but not axons. These distinct modes of interneuron firing can trigger both a rapid and a sustained component of inhibition onto TC neurons. Our studies suggest that active conductances make LGN interneurons flexible circuit-elements that can shift their spatial and temporal properties of GABA release in response to coincident activation of functionally related subsets of RGCs.     

Dendritic development in the dorsal lateral geniculate nucleus of ferrets in the postnatal absence of retinal input: A golgi study

In order to determine the ongoing role of retinal fibers in the development of dorsal lateral geniculate nucleus (dLGN) neurons during postnatal development, the development of dLGN neurons in the postnatal absence of retinal input was studied in pigmented ferrets using the Golgi-Hortega technique. The development of four dLGN cell classes, defined on the basis of somatic and dendritic morphology, was described previously in normal ferrets (Sutton and Brunso-Bechtold, 1991, J. Comp. Neurol. 309 : 71–85). The present results indicate that the morphological development of dLGN neurons is strikingly similar in normal and experimental ferrets. The exuberant dendritic appendages that appear after eye opening in normal ferrets are overproduced and eliminated in the postnatal absence of retinal input; however, the final reduction of these transient appendages is delayed. Because exuberant appendages develop in the absence of retinal input, their production cannot depend upon visual experience. Differences in cell body size between normal and experimental ferrets are apparent only after neurons can be classified at the end of the first postnatal month. Cell body size is markedly reduced for class 1 neurons; class 2 cells also are reduced in size but to a far lesser extent. As there is a general trend for class 1 neurons to have the functional properties of Y-cells, it is likely that the dLGN neurons most affected by the absence of retinal input also are Y-cells. © 1993 John Wiley & Sons, Inc.     

Bursting as an Effective Relay Mode in a Minimal Thalamic Model

In recent years, accumulating evidence indicates that thalamic bursts are present during wakefulness and participate in information transmission as an effective relay mode with distinctive properties from the tonic activity. Thalamic bursts originate from activation of the low threshold calcium cannels via a local feedback inhibition, exerted by the thalamic reticular neurons upon the relay neurons. This article, examines if this simple mechanism is sufficient to explain the distinctive properties of thalamic bursting as an effective relay mode. A minimal model of thalamic circuit composed of a retinal spike train, a relay neuron and a reticular neuron is simulated to generate the tonic and burst firing modes. The integrate-and-fire-or-burst model is used to simulate the neurons. After discriminating the burst events with criteria based on inter-spike-intervals, statistical indices show that the bursts of the minimal model are stereotypic events. The relation between the rate of bursts and the parameters of the input spike train demonstrates marked nonlinearities. Burst response is shown to be selective to spike-silence-spike sequences in the input spike train. Moreover, burst events represent the input more reliably than the tonic spike in a considerable range of the parameters of the model. In conclusion, many of the distinctive properties of thalamic bursts such as stereotypy, nonlinear dependence on the sensory stimulus, feature selectivity and reliability are reproducible in the minimal model. Furthermore, the minimal model predicts that while the bursts are more frequent in the spike train of the off-center X relay neurons (corresponding to off-center X retinal ganglion cells), they are more reliable when generated by the on-center ones (corresponding to on-center X ganglion cells).     

A computational model of how an interaction between the thalamocortical and thalamic reticular neurons transforms the low-frequency oscillations of the globus pallidus

In Parkinson’s disease, neurons of the internal segment of the globus pallidus (GPi) display the low-frequency tremor-related oscillations. These oscillatory activities are transmitted to the thalamic relay nuclei. Computer models of the interacting thalamocortical (TC) and thalamic reticular (RE) neurons were used to explore how the TC-RE network processes the low-frequency oscillations of the GPi neurons. The simulation results show that, by an interaction between the TC and RE neurons, the TC-RE network transforms a low-frequency oscillatory activity of the GPi neurons to a higher frequency of oscillatory activity of the TC neurons (the superharmonic frequency transformation). In addition to the interaction between the TC and RE neurons, the low-threshold calcium current in the RE and TC neurons and the hyperpolarization-activated cation current (I _h) in the TC neurons have significant roles in the superharmonic frequency transformation property of the TC-RE network. The external globus pallidus (GPe) oscillatory activity, which is directly transmitted to the RE nucleus also displays a significant modulatory effect on the superharmonic frequency transformation property of the TC-RE network. Action Editor: John Rinzel     

Synchronization of Isolated Downstates (K-Complexes) May Be Caused by Cortically-Induced Disruption of Thalamic Spindling

Sleep spindles and K-complexes (KCs) define stage 2 NREM sleep (N2) in humans. We recently showed that KCs are isolated downstates characterized by widespread cortical silence. We demonstrate here that KCs can be quasi-synchronous across scalp EEG and across much of the cortex using electrocorticography (ECOG) and localized transcortical recordings (bipolar SEEG). We examine the mechanism of synchronous KC production by creating the first conductance based thalamocortical network model of N2 sleep to generate both spontaneous spindles and KCs. Spontaneous KCs are only observed when the model includes diffuse projections from restricted prefrontal areas to the thalamic reticular nucleus (RE), consistent with recent anatomical findings in rhesus monkeys. Modeled KCs begin with a spontaneous focal depolarization of the prefrontal neurons, followed by depolarization of the RE. Surprisingly, the RE depolarization leads to decreased firing due to disrupted spindling, which in turn is due to depolarization-induced inactivation of the low-threshold Ca²⁺ current (I_T). Further, although the RE inhibits thalamocortical (TC) neurons, decreased RE firing causes decreased TC cell firing, again because of disrupted spindling. The resulting abrupt removal of excitatory input to cortical pyramidal neurons then leads to the downstate. Empirically, KCs may also be evoked by sensory stimuli while maintaining sleep. We reproduce this phenomenon in the model by depolarization of either the RE or the widely-projecting prefrontal neurons. Again, disruption of thalamic spindling plays a key role. Higher levels of RE stimulation also cause downstates, but by directly inhibiting the TC neurons. SEEG recordings from the thalamus and cortex in a single patient demonstrated the model prediction that thalamic spindling significantly decreases before KC onset. In conclusion, we show empirically that KCs can be widespread quasi-synchronous cortical downstates, and demonstrate with the first model of stage 2 NREM sleep a possible mechanism whereby this widespread synchrony may arise.     

大鼠海马癫痫电网络重建中爆发式放电神经元的活动

本文探讨双侧海马(hippoeampus,HPC)神经网络中爆发式放电神经元(bursting-firing neurons,BFN)的活动规律及其与海马癫痫网络重建的关系。实验用雄性SD大鼠140只(150-250 g),急性强直电刺激(60 Hz,2 s,0.4-0.6 mA)右后背HPC CAl区(acute tetanization of the posterior dorsal hippocampus,ATPDH),同步记录同侧或对侧前背HPC单位放电和深部电图;强直电刺激右前背HPC(acute tetanization of the anterior dorsal hippocampus,AT-ADH),同步记录双侧前背HPC单位放电。实验共记录了13.8％(19/138)双侧前背HPC的BFN,其中13个为刺激诱发性BFN,6个为自发性BFN。强直电刺激引起的诱发反应包括:(1)ATPDH明显调制同侧前背HPC的BFN,产生规则的节律性爆发式放电,刺激后串内动作电位间期(bursting interspike interval,BISI)减小(P<0.001);(2)AT-PDH引起对侧前背HPC的BFN出现抑制后轻度调制效应,刺激后动作电位间期(interspike interval,ISI)增大(P<0.001);(3)ATADH后易化对侧前背HPC的自发性BFN节律,增加ISI(P<0.001)和IBI(P=0.01);(4)ATPDH诱导双侧前背HPC的BFN产生规则的节律性爆发式放电,伴有同步或非同步性网络癫痫的形成。上述实验结果提示,ATPDH沿同侧HPC长轴,跨大脑半球诱发前背HPC单个BFN的形成,其节律性爆     

Intrinsic bursting activity in the pre-Bötzinger Complex: Role of persistent sodium and potassium currents

Computational models of single pacemaker neuron and neural population in the pre-Bötzinger Complex (pBC) were developed based on the previous models by Butera et al. (1999a,b). Our modeling study focused on the conditions that could define endogenous bursting vs. tonic activity in single pacemaker neurons and population bursting vs. asynchronous firing in populations of pacemaker neurons. We show that both bursting activity in single pacemaker neurons and population bursting activity may be released or suppressed depending on the expression of persistent sodium (I_NaP) and delayed-rectifier potassium (I_K) currents. Specifically, a transition from asynchronous firing to population bursting could be induced by a reduction of I_K via a direct suppression of the potassium conductance or through an elevation of extracellular potassium concentration. Similar population bursting activity could be triggered by an augmentation of I_NaP. These findings are discussed in the context of the possible role of population bursting activity in the pBC in the respiratory rhythm generation in vivo vs. in vitro and during normal breathing in vivo vs. gasping.     

The influence of the corticothalamic projection on responses in thalamus and cortex

We review results on the in vivo properties of neurons in the dorsal lateral geniculate nucleus (dLGN) that receives its afferent input from the retina and projects to the visual cortex. In addition, the dLGN receives input from the brain stem and from a rather strong corticothalamic back-projection, which originates in layer 6 of the visual cortex. We compare the behaviour of dLGN cells during spontaneous changes of the frequency contents of the electroencephalograph (EEG) (which are mainly related to a changing brain stem influence), with those that are obtained when experimentally silencing the corticothalamic feedback. The spatial and temporal response properties of dLGN cells are compared during these two conditions, and we report that the neurons behave similarly during a synchronized EEG state and during inactive corticothalamic feedback. In both situations, dLGN cells are rather phasic and their remaining tonic activity is temporally dispersed, indicating a hyperpolarizing effect. By means of a novel method, we were able to chronically eliminate a large proportion of the corticothalamic projection neurons from the otherwise intact cortex. In this condition, we found that cortical cells also lose their EEG specific response differences but, in this instance, probably due to a facilitatory (depolarizing) plasticity reaction of the remaining network.     

Synchronous and asynchronous bursting states: role of intrinsic neural dynamics

Brain signals such as local field potentials often display gamma-band oscillations (30-70 Hz) in a variety of cognitive tasks. These oscillatory activities possibly reflect synchronization of cell assemblies that are engaged in a cognitive function. A type of pyramidal neurons, i.e., chattering neurons, show fast rhythmic bursting (FRB) in the gamma frequency range, and may play an active role in generating the gamma-band oscillations in the cerebral cortex. Our previous phase response analyses have revealed that the synchronization between the coupled bursting neurons significantly depends on the bursting mode that is defined as the number of spikes in each burst. Namely, a network of neurons bursting through a Ca(2+)-dependent mechanism exhibited sharp transitions between synchronous and asynchronous firing states when the neurons exchanged the bursting mode between singlet, doublet and so on. However, whether a broad class of bursting neuron models commonly show such a network behavior remains unclear. Here, we analyze the mechanism underlying this network behavior using a mathematically tractable neuron model. Then we extend our results to a multi-compartment version of the NaP current-based neuron model and prove a similar tight relationship between the bursting mode changes and the network state changes in this model. Thus, the synchronization behavior couples tightly to the bursting mode in a wide class of networks of bursting neurons.     

Ca2+-induced Ca2+ release supports the relay mode of activity in thalamocortical cells

Ca2+ ions play an important role during rhythmic bursting of thalamocortical neurons within sleep. The function of Ca2+ during the tonic relay mode of these neurons during wakefulness is less clear. Here, we report that tonic activity in thalamocortical cells results in an increase in the intracellular Ca2+ concentration and subsequent release of Ca2+ from intracellular stores mediated via ryanodine receptors (RyRs). Blockade of Ca2+ release shifted the regular firing of single action potentials toward the generation of spike clusters. Regular spike firing and intracellular Ca2+ release thus appear to be functionally coupled in a positive feedback manner, thereby supporting the relay mode of thalamocortical cells during wakefulness. Regulatory influences may be coupled to this system via the cyclic ADP ribose pathway.     

Changing excitation and inhibition in simulated neural networks: effects on induced bursting behavior

 The development of synchronous bursting in neuronal ensembles represents an important change in network behavior. To determine the influences on development of such synchronous bursting behavior we study the dynamics of small networks of sparsely connected excitatory and inhibitory neurons using numerical simulations. The synchronized bursting activities in networks evoked by background spikes are investigated. Specifically, patterns of bursting activity are examined when the balance between excitation and inhibition on neuronal inputs is varied and the fraction of inhibitory neurons in the network is changed. For quantitative comparison of bursting activities in networks, measures of the degree of synchrony are used. We demonstrate how changes in the strength of excitation on inputs of neurons can be compensated by changes in the strength of inhibition without changing the degree of synchrony in the network. The effects of changing several network parameters on the network activity are analyzed and discussed. These changes may underlie the transition of network activity from normal to potentially pathologic (e.g., epileptic) states. Received: 21 May 2002 / Accepted in revised form: 3 December 2002 / Published online: 7 March 2003 Correspondence to: P. Kudela (e-mail: pkudela@jhmi.edu) Acknowledgements. This research was supported by NIH grant NS 38958.     

Patterns of correlated spontaneous bursting activity in the developing mammalian retina

The adult visual system is highly organized in its patterns of connectivity. Connections between the retina and its central target, the dorsal lateral geniculate nucleus (dLGN), are remodeled during development as inappropriate synaptic inputs are eliminated by a process that requires retinal activity. Multineuronal recordings of the neonatal ferret retina reveal that during the refinement period, retinal ganglion cells spontaneously display rhythmic bursting activity in which the bursts of neighboring cells are correlated by propagating excitatory waves. These spontaneous retinal waves have temporal and spatial properties that appear instructive for the refinement of the early patterns of retinogeniculate connections prior to visual stimulation.     

Rhythmic activities of hypothalamic magnocellular neurons: Autocontrol mechanisms

Electrophysiological recordings in lactating rats show that oxytocin (OT) and vasopressin (AVP) neurons exhibit specific patterns of activities in relation to peripheral stimuli: periodic bursting firing for OT neurons during suckling, phasic firing for AVP neurons during hyperosmolarity (systemic injection of hypertonic saline). These activities are autocontrolled by OT and AVP released somato-dentritically within the hypothalamic magnocellular nuclei. In vivo, OT enhances the amplitude and frequency of bursts, an effect accompanied with an increase in basal firing rate. However, the characteristics of firing change as facilitation proceeds: the spike patterns become very irregular with clusters of spikes spaced by long silences; the firing rate is highly variable and clearly oscillates before facilitated bursts. This unstable behaviour dramatically decreases during intense tonic activation which temporarily interrupts bursting, and could therefore be a prerequisite for bursting. In vivo, the effects of AVP depend on the initial firing pattern of AVP neurons: AVP excites weakly active neurons (increasing duration of active periods and decreasing silences), inhibits highly active neurons, and does not affect neurons with intermediate phasic activity. AVP brings the entire population of AVP neurons to discharge with a medium phasic activity characterised by periods of firing and silence lasting 20–40 s, a pattern shown to optimise the release of AVP from the neurohypophysis. Each of the peptides (OT or AVP) induces an increase in intracellular Ca²⁺ concentration, specifically in the neurons containing either OT or AVP respectively. OT evokes the release of Ca²⁺ from IP3-sensitive intracellular stores. AVP induces an influx of Ca²⁺ through voltage-dependent Ca²⁺ channels of T-, L- and N-types. We postulate that the facilitatory autocontrol of OT and AVP neurons could be mediated by Ca²⁺ known to play a key role in the control of the patterns of phasic neurons.     

Encoding whisker deflection velocity within the rodent barrel cortex using phase-delayed inhibition

The primary sensory feature represented within the rodent barrel cortex is the velocity with which a whisker has been deflected. Whisker deflection velocity is encoded within the thalamus via population synchrony (higher deflection velocities entail greater synchrony among the corresponding thalamic population). Thalamic (TC) cells project to regular spiking (RS) cells within the barrel cortex, as well as to inhibitory cortical fast-spiking (FS) neurons, which in turn project to RS cells. Thus, TC spikes result in EPSPs followed, with a small time lag, by IPSPs within an RS cell, and hence the RS cell decodes TC population synchrony by employing a phase-delayed inhibition synchrony detection scheme. As whisker deflection velocity is increased, the probability that an RS cell spikes rises, while jitter in the timing of RS cell spikes remains constant. Furthermore, repeated whisker deflections with fixed velocity lead to system adaptation – TC →RS, TC →FS, and FS →RS synapses all weaken substantially, leading to a smaller probability of spiking of the RS cell and increased jitter in the timing of RS cell spikes. Interestingly, RS cell activity is better able to distinguish among different whisker deflection velocities after adaptation. In this work, we construct a biophysical model of a basic ‘building block’ of barrel cortex – the feedforward circuit consisting of TC cells, FS cells, and a single RS cell – and we examine the ability of the purely feedforward circuit to explain the experimental data on RS cell spiking probability, jitter, adaptation, and deflection velocity discrimination. Moreover, we study the contribution of the phase-delayed inhibition network structure to the ability of an RS cell to decode whisker deflection velocity encoded via TC population synchrony.     

Single Versus Repetitive Spiking to the Current Stimulus of A-β Mechanosensitive Neurons in the Crotaline Snake Trigeminal Ganglion

1. Intrasomal recordings of potentials produced by current stimulation in vivo were made from 24 (A-) touch and 19 vibrotactile neurons in the trigeminal ganglion of 29 crotaline snakes, Trimeresurus flavoviridis. 2. Usually touch neurons responded with a single action potential at the beginning of a prolonged depolarizing pulse, whereas all vibrotactile neurons responded with multiple spikes.3. The electrophysiological parameters examined were membrane potential, threshold current, input resistance and capacitance, time constant, rebound latency, and its threshold current. Touch neurons had higher input resistance (and lower input capacitance) than vibrotactile neurons.4. In conclusion, current injection, which elicits a single or multiple spiking, seems a useful way to separate touch neurons from vibrotactile neurons without confirming the receptor response, and some membrane properties are also specific to the sensory modality.     

Optogenetic Stimulation of the Corticothalamic Pathway Affects Relay Cells and GABAergic Neurons Differently in the Mouse Visual Thalamus

The dorsal lateral geniculate nucleus (dLGN) serves as the primary conduit of retinal information to visual cortex. In addition to retinal input, dLGN receives a large feedback projection from layer VI of visual cortex. Such input modulates thalamic signal transmission in different ways that range from gain control to synchronizing network activity in a stimulus-specific manner. However, the mechanisms underlying such modulation have been difficult to study, in part because of the complex circuitry and diverse cell types this pathway innervates. To address this and overcome some of the technical limitations inherent in studying the corticothalamic (CT) pathway, we adopted a slice preparation in which we were able to stimulate CT terminal arbors in the visual thalamus of the mouse with blue light by using an adeno-associated virus to express the light-gated ion channel, ChIEF, in layer VI neurons. To examine the postsynaptic responses evoked by repetitive CT stimulation, we recorded from identified relay cells in dLGN, as well as GFP expressing GABAergic neurons in the thalamic reticular nucleus (TRN) and intrinsic interneurons of dLGN. Relay neurons exhibited large glutamatergic responses that continued to increase in amplitude with each successive stimulus pulse. While excitatory responses were apparent at postnatal day 10, the strong facilitation noted in adult was not observed until postnatal day 21. GABAergic neurons in TRN exhibited large initial excitatory responses that quickly plateaued during repetitive stimulation, indicating that the degree of facilitation was much larger for relay cells than for TRN neurons. The responses of intrinsic interneurons were smaller and took the form of a slow depolarization. These differences in the pattern of excitation for different thalamic cell types should help provide a framework for understanding how CT feedback alters the activity of visual thalamic circuitry during sensory processing as well as different behavioral or pathophysiological states.     

A thalamo-cortical neural mass model for the simulation of brain rhythms during sleep

Cortico-thalamic interactions are known to play a pivotal role in many brain phenomena, including sleep, attention, memory consolidation and rhythm generation. Hence, simple mathematical models that can simulate the dialogue between the cortex and the thalamus, at a mesoscopic level, have a great cognitive value. In the present work we describe a neural mass model of a cortico-thalamic module, based on neurophysiological mechanisms. The model includes two thalamic populations (a thalamo-cortical relay cell population, TCR, and its related thalamic reticular nucleus, TRN), and a cortical column consisting of four connected populations (pyramidal neurons, excitatory interneurons, inhibitory interneurons with slow and fast kinetics). Moreover, thalamic neurons exhibit two firing modes: bursting and tonic. Finally, cortical synapses among pyramidal neurons incorporate a disfacilitation mechanism following prolonged activity. Simulations show that the model is able to mimic the different patterns of rhythmic activity in cortical and thalamic neurons (beta and alpha waves, spindles, delta waves, K-complexes, slow sleep waves) and their progressive changes from wakefulness to deep sleep, by just acting on modulatory inputs. Moreover, simulations performed by providing short sensory inputs to the TCR show that brain rhythms during sleep preserve the cortex from external perturbations, still allowing a high cortical activity necessary to drive synaptic plasticity and memory consolidation. In perspective, the present model may be used within larger cortico-thalamic networks, to gain a deeper understanding of mechanisms beneath synaptic changes during sleep, to investigate the specific role of brain rhythms, and to explore cortical synchronization achieved via thalamic influences.     

Two types of independent bursting mechanisms in inspiratory neurons: an integrative model

The network of coupled neurons in the pre-Bötzinger complex (pBC) of the medulla generates a bursting rhythm, which underlies the inspiratory phase of respiration. In some of these neurons, bursting persists even when synaptic coupling in the network is blocked and respiratory rhythmic discharge stops. Bursting in inspiratory neurons has been extensively studied, and two classes of bursting neurons have been identified, with bursting mechanism depends on either persistent sodium current or changes in intracellular Ca²⁺, respectively. Motivated by experimental evidence from these intrinsically bursting neurons, we present a two-compartment mathematical model of an isolated pBC neuron with two independent bursting mechanisms. Bursting in the somatic compartment is modeled via inactivation of a persistent sodium current, whereas bursting in the dendritic compartment relies on Ca²⁺ oscillations, which are determined by the neuromodulatory tone. The model explains a number of conflicting experimental results and is able to generate a robust bursting rhythm, over a large range of parameters, with a frequency adjusted by neuromodulators.