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.     

Propofol suppresses synaptic responsiveness of somatosensory relay neurons to excitatory input by potentiating GABAA receptor chloride channels

Propofol is a widely used intravenous general anesthetic. Propofol-induced unconsciousness in humans is associated with inhibition of thalamic activity evoked by somatosensory stimuli. However, the cellular mechanisms underlying the effects of propofol in thalamic circuits are largely unknown. We investigated the influence of propofol on synaptic responsiveness of thalamocortical relay neurons in the ventrobasal complex (VB) to excitatory input in mouse brain slices, using both current- and voltage-clamp recording techniques. Excitatory responses including EPSP temporal summation and action potential firing were evoked in VB neurons by electrical stimulation of corticothalamic fibers or pharmacological activation of glutamate receptors. Propofol (0.6 – 3 μM) suppressed temporal summation and spike firing in a concentration-dependent manner. The thalamocortical suppression was accompanied by a marked decrease in both EPSP amplitude and input resistance, indicating that a shunting mechanism was involved. The propofol-mediated thalamocortical suppression could be blocked by a GABA_A receptor antagonist or chloride channel blocker, suggesting that postsynaptic GABA_A receptors in VB neurons were involved in the shunting inhibition. GABA_A receptor-mediated inhibitory postsynaptic currents (IPSCs) were evoked in VB neurons by electrical stimulation of the reticular thalamic nucleus. Propofol markedly increased amplitude, decay time, and charge transfer of GABA_A IPSCs. The results demonstrated that shunting inhibition of thalamic somatosensory relay neurons by propofol at clinically relevant concentrations is primarily mediated through the potentiation of the GABA_A receptor chloride channel-mediated conductance, and such inhibition may contribute to the impaired thalamic responses to sensory stimuli seen during propofol-induced anesthesia.     

Hypocretin (orexin) induces calcium transients in single spines postsynaptic to identified thalamocortical boutons in prefrontal slice

In vivo, thalamocortical axons are susceptible to the generation of terminal spikes which antidromically promote bursting in the thalamus. Although neurotransmitters could elicit such ectopic action potentials at thalamocortical boutons, this hypothesis has never been confirmed. Prefrontal cortex is the cortical area most implicated in arousal and is innervated by thalamic neurons that are unusual since they burst rhythmically during waking. We show that a neurotransmitter critical for alertness, hypocretin (orexin), directly excites prefrontal thalamocortical synapses in acute slice. This TTX-sensitive activation of thalamic axons was demonstrated electrophysiologically and by two-photon sampling of calcium transients at single spines in apposition to thalamic boutons anterogradely labeled in vivo. Spines receiving these long-range projections constituted a unique population in terms of the presynaptic excitatory action of hypocretin. By this mechanism, the hypocretin projection to prefrontal cortex may play a larger role in prefrontal or "executive" aspects of alertness and attention than previously anticipated.     

A multivariate population density model of the dLGN/PGN relay

Using a population density approach we study the dynamics of two interacting collections of integrate-and-fire-or-burst (IFB) neurons representing thalamocortical (TC) cells from the dorsal lateral geniculate nucleus (dLGN) and thalamic reticular (RE) cells from the perigeniculate nucleus (PGN). Each population of neurons is described by a multivariate probability density function that satisfies a conservation equation with appropriately defined probability fluxes and boundary conditions. The state variables of each neuron are the membrane potential and the inactivation gating variable of the low-threshold Ca²⁺ current I_T. The synaptic coupling of the populations and external excitatory drive are modeled by instantaneous jumps in the membrane potential of postsynaptic neurons. The population density model is validated by comparing its response to time-varying retinal input to Monte Carlo simulations of the corresponding IFB network composed of 100 to 1000 cells per population. In the absence of retinal input, the population density model exhibits rhythmic bursting similar to the 7 to 14 Hz oscillations associated with slow wave sleep that require feedback inhibition from RE to TC cells. When the TC and RE cell potassium leakage conductances are adjusted to represent cholingergic neuromodulation and arousal of the network, rhythmic bursting of the probability density model may either persists or be eliminated depending on the number of excitatory (TC to RE) or inhibitory (RE to TC) connections made by each presynaptic cell. When the probability density model is stimulated with constant retinal input (10–100 spikes/sec), a wide range of responses are observed depending on cellular parameters and network connectivity. These include asynchronous burst and tonic spikes, sleep spindle-like rhythmic bursting, and oscillations in population firing rate that are distinguishable from sleep spindles due to their amplitude, frequency, or the presence of tonic spikes. In this context of dLGN/PGN network modeling, we find the population density approach using 2,500 mesh points and resolving membrane voltage to 0.7 mV is over 30 times more efficient than 1000-cell Monte Carlo simulations. Action Editor: David Golomb     

Self-sustained asynchronous irregular states and Up–Down states in thalamic,cortical and thalamocortical networks of nonlinear integrate-and-fire neurons

Randomly-connected networks of integrate-and-fire (IF) neurons are known to display asynchronous irregular (AI) activity states, which resemble the discharge activity recorded in the cerebral cortex of awake animals. However, it is not clear whether such activity states are specific to simple IF models, or if they also exist in networks where neurons are endowed with complex intrinsic properties similar to electrophysiological measurements. Here, we investigate the occurrence of AI states in networks of nonlinear IF neurons, such as the adaptive exponential IF (Brette-Gerstner-Izhikevich) model. This model can display intrinsic properties such as low-threshold spike (LTS), regular spiking (RS) or fast-spiking (FS). We successively investigate the oscillatory and AI dynamics of thalamic, cortical and thalamocortical networks using such models. AI states can be found in each case, sometimes with surprisingly small network size of the order of a few tens of neurons. We show that the presence of LTS neurons in cortex or in thalamus, explains the robust emergence of AI states for relatively small network sizes. Finally, we investigate the role of spike-frequency adaptation (SFA). In cortical networks with strong SFA in RS cells, the AI state is transient, but when SFA is reduced, AI states can be self-sustained for long times. In thalamocortical networks, AI states are found when the cortex is itself in an AI state, but with strong SFA, the thalamocortical network displays Up and Down state transitions, similar to intracellular recordings during slow-wave sleep or anesthesia. Self-sustained Up and Down states could also be generated by two-layer cortical networks with LTS cells. These models suggest that intrinsic properties such as adaptation and low-threshold bursting activity are crucial for the genesis and control of AI states in thalamocortical networks.     

Inhibition modifies the effects of slow calcium-activated potassium channels on epileptiform activity in a neuronal network model

Generation of epileptiform activity typically results from a change in the balance between network excitation and inhibition. Experimental evidence indicates that alterations of either synaptic activity or intrinsic membrane properties can produce increased network excitation. The slow Ca²⁺-activated K⁺ currents (sI_AHP) are important modulators of neuronal firing rate and excitability and have important established and potential roles in epileptogenesis. While the effects of changes in sI_AHP on individual neuronal excitability are readily studied and well established, the effects of such changes on network behavior are less well known. The experiments here utilize a defined small network model of multicompartment pyramidal cells and an inhibitory interneuron to study the effects of changes in sI_AHP on network behavior. The benefits of this model system include the ability to observe activity in all cells in a network and the effects of interactions of multiple simultaneous influences. In the model with no inhibitory interneuron, increasing sI_AHP results in progressively decreasing burst activity. Adding an inhibitory interneuron changes the observed effects; at modest inhibitory strengths, increasing sI_AHP in all network neurons actually results in increased network bursting (except at very high values). The duration of the burst activity is influenced by the length of delay in a feedback loop, with longer loops resulting in more prolonged bursting. These observations illustrate that the study of potential antiepileptogenic membrane effects must be extended to realistic networks. Network inhibition can dramatically alter the observations seen in pure excitatory networks.     

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     

A mathematical model of adult GnRH neurons in mouse brain and its bifurcation analysis

GnRH neurons are hypothalamic neurons that secrete gonadotropin-releasing hormone (GnRH) which stimulates the release of gonadotropins, one of the crucial hormones for sexual development, fertility and maturation. A mathematical model was built to help elucidate the mechanisms underlying electrical bursting and synchronous [Ca²⁺] transients in GnRH neurons (Lee et al., 2010). The model predicted that bursting in GnRH neurons (at least of the short-bursting type) requires the existence of a [Ca²⁺]-dependent slow after-hyperpolarisation current (sI_AHP-UCL), and this predicted current was found experimentally. GnRH behaviour under a wide range of conditions (inhibition of Na⁺ channels, IP₃ receptors, [Ca²⁺]-dependent K⁺ channels, or Ca²⁺ pumps, or in the presence of zero extracellular [Ca²⁺]) is successfully reproduced by the model. In this paper, a simplified version of the previous model, with the same qualitative behaviour, is constructed and studied using timescale separation techniques and bifurcation analysis.     

Somatosensory integration controlled by dynamic thalamocortical feed-forward inhibition

The temporal features of tactile stimuli are faithfully represented by the activity of neurons in the somatosensory cortex. However, the cellular mechanisms that enable cortical neurons to report accurate temporal information are not known. Here, we show that in the rodent barrel cortex, the temporal window for integration of thalamic inputs is under the control of thalamocortical feed-forward inhibition and can vary from 1 to 10 ms. A single thalamic fiber can trigger feed-forward inhibition and contacts both excitatory and inhibitory cortical neurons. The dynamics of feed-forward inhibition exceed those of each individual synapse in the circuit and are captured by a simple disynaptic model of the thalamocortical projection. The variations in the integration window produce changes in the temporal precision of cortical responses to whisker stimulation. Hence, feed-forward inhibitory circuits, classically known to sharpen spatial contrast of tactile inputs, also increase the temporal resolution in the somatosensory cortex.     

Perisomatic GABA release and thalamocortical integration onto neocortical excitatory cells are regulated by neuromodulators

Neuromodulators such as acetylcholine, serotonin, and noradrenaline are powerful regulators of neocortical activity. Although it is well established that cortical inhibition is the target of these modulations, little is known about their effects on GABA release from specific interneuron types. This knowledge is necessary to gain a mechanistic understanding of the actions of neuromodulators because different interneuron classes control specific aspects of excitatory cell function. Here, we report that GABA release from fast-spiking (FS) cells, the most prevalent interneuron subtype in neocortex, is robustly inhibited following activation of muscarinic, serotonin, adenosine, and GABA(B) receptors--an effect that regulates FS cell control of excitatory neuron firing. The potent muscarinic inhibition of GABA release from FS cells suppresses thalamocortical feedforward inhibition. This is supplemented by the muscarinic-mediated depolarization of thalamo-recipient excitatory neurons and the nicotinic enhancement of thalamic input onto these neurons to promote thalamocortical excitation.     

Mechanisms of Generation of Pacemaker Activity by Identified Helix Neurons

We studied the mechanisms of generation of pacemaker activity in identified neurons of Helix pomatia. For this purpose, we isolated the PPa2 and PPa7 neurons generating spontaneous rhythmic monomodal activity and PPa1 neuron with bursting activity. It was demonstrated that isolated PPa2 and PPa7 cells produce endogenous rhythmic activity that was not considerably modified by external application of 1 mM CdCl₂. Sometimes, only low-amplitude dendritic action potentials (AP) were observed instead of generation of full-amplitude somatic AP. In contrast, isolation of the PPa1 neuron eliminated its bursting activity, but subsequent application of oxytocin on this neuron recovered such activity. This finding shows that the bursting activity of the PPa1 neuron is of an exogenous nature. Application of 1 mM CdCl₂ suppressed this bursting activity, but when Cd²⁺ was applied against the background of superfusion of the neuron with Ringer solution containing a bursting activity-initiating neuropeptide obtained from the molluscan CNS, this blocker was incapable of suppressing the bursting activity. A blocker of the hyperpolarization-activated current (I _h , H current), Cs⁺ (10 mM) exerted no noticeable effect on the activity of the studied neurons. Our findings allow us to conclude that the pacemaker activity is initiated within the dendritic tree of a cell and is then electrotonically spread to the soma, where full-amplitude AP are generated. It seems probable that Ca²⁺ ions and H current are not directly involved in generation of the pacemaker activity in the studied snail neurons.     

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.     

Evidence that descending cortical axons are essential for thalamocortical axons to cross the pallial-subpallial boundary in the embryonic forebrain

Developing thalamocortical axons traverse the subpallium to reach the cortex located in the pallium. We tested the hypothesis that descending corticofugal axons are important for guiding thalamocortical axons across the pallial-subpallial boundary, using conditional mutagenesis to assess the effects of blocking corticofugal axonal development without disrupting thalamus, subpallium or the pallial-subpallial boundary. We found that thalamic axons still traversed the subpallium in topographic order but did not cross the pallial-subpallial boundary. Co-culture experiments indicated that the inability of thalamic axons to cross the boundary was not explained by mutant cortex developing a long-range chemorepulsive action on thalamic axons. On the contrary, cortex from conditional mutants retained its thalamic axonal growth-promoting activity and continued to express Nrg-1, which is responsible for this stimulatory effect. When mutant cortex was replaced with control cortex, corticofugal efferents were restored and thalamic axons from conditional mutants associated with them and crossed the pallial-subpallial boundary. Our study provides the most compelling evidence to date that cortical efferents are required to guide thalamocortical axons across the pallial-subpallial boundary, which is otherwise hostile to thalamic axons. These results support the hypothesis that thalamic axons grow from subpallium to cortex guided by cortical efferents, with stimulation from diffusible cortical growth-promoting factors.     

Modeling shifts in the rate and pattern of subthalamopallidal network activity during deep brain stimulation

Deep brain stimulation (DBS) of the subthlamic nucleus (STN) represents an effective treatment for medically refractory Parkinson’s disease; however, understanding of its effects on basal ganglia network activity remains limited. We constructed a computational model of the subthalamopallidal network, trained it to fit in vivo recordings from parkinsonian monkeys, and evaluated its response to STN DBS. The network model was created with synaptically connected single compartment biophysical models of STN and pallidal neurons, and stochastically defined inputs driven by cortical beta rhythms. A least mean square error training algorithm was developed to parameterize network connections and minimize error when compared to experimental spike and burst rates in the parkinsonian condition. The output of the trained network was then compared to experimental data not used in the training process. We found that reducing the influence of the cortical beta input on the model generated activity that agreed well with recordings from normal monkeys. Further, during STN DBS in the parkinsonian condition the simulations reproduced the reduction in GPi bursting found in existing experimental data. The model also provided the opportunity to greatly expand analysis of GPi bursting activity, generating three major predictions. First, its reduction was proportional to the volume of STN activated by DBS. Second, GPi bursting decreased in a stimulation frequency dependent manner, saturating at values consistent with clinically therapeutic DBS. And third, ablating STN neurons, reported to generate similar therapeutic outcomes as STN DBS, also reduced GPi bursting. Our theoretical analysis of stimulation induced network activity suggests that regularization of GPi firing is dependent on the volume of STN tissue activated and a threshold level of burst reduction may be necessary for therapeutic effect.     

Convergence between projections from different thalamic nuclei to columns of the rat somatosensory cortex

Primuline fluorochrome retrograde transport technique was used to investigate sources of thalamocortical projections to a single rat somatosensory cortex column connected with the projection of the C₃ vibrissa. Labeled cells were identified in eight different thalamic nuclei: two specific, five nonspecific, and one association nucleus. Labeled neurons differed in the degree of stain accumulated as well as cell numbers and density of distribution from one nucleus to another, indicative of the different arborization patterns of their axons within the cortex. Highest numbers of heavily stained cells as well as highest density of distribution were observed in the ventral thalamic nucleus. The convergence seen between different thalamocortical inputs on to a single somatosensory cortex column explains the functional differences observed between neurons belonging to the same column and makes the formation of functionally distinct neuronal groupings appear possible on this structural basis.Neurocybernetics Research Institute, Rostov-on-Don. Translated from Neirofiziologiya, Vol. 21, No. 2, pp. 168–174, March–April, 1989.     

The effects of DBS patterns on basal ganglia activity and thalamic relay

Thalamic neurons receive inputs from cortex and their responses are modulated by the basal ganglia (BG). This modulation is necessary to properly relay cortical inputs back to cortex and downstream to the brain stem when movements are planned. In Parkinson's disease (PD), the BG input to thalamus becomes pathological and relay of motor-related cortical inputs is compromised, thereby impairing movements. However, high frequency (HF) deep brain stimulation (DBS) may be used to restore relay reliability, thereby restoring movements in PD patients. Although therapeutic, HF stimulation consumes significant power forcing surgical battery replacements, and may cause adverse side effects. Here, we used a biophysical-based model of the BG-Thalamus motor loop in both healthy and PD conditions to assess whether low frequency stimulation can suppress pathological activity in PD and enable the thalamus to reliably relay movement-related cortical inputs. We administered periodic pulse train DBS waveforms to the sub-thalamic nucleus (STN) with frequencies ranging from 0-140 Hz, and computed statistics that quantified pathological bursting, oscillations, and synchronization in the BG as well as thalamic relay of cortical inputs. We found that none of the frequencies suppressed all pathological activity in BG, though the HF waveforms recovered thalamic reliability. Our rigorous study, however, led us to a novel DBS strategy involving low frequency multi-input phase-shifted DBS, which successfully suppressed pathological symptoms in all BG nuclei and enabled reliable thalamic relay. The neural restoration remained robust to changes in the model parameters characterizing early to late PD stages.     

A unified model for two modes of bursting in GnRH neurons

Gonadotropin-releasing hormone (GnRH) neurons exhibit at least two intrinsic modes of action potential burst firing, referred to as parabolic and irregular bursting. Parabolic bursting is characterized by a slow wave in membrane potential that can underlie periodic clusters of action potentials with increased interspike interval at the beginning and at the end of each cluster. Irregular bursting is characterized by clusters of action potentials that are separated by varying durations of interburst intervals and a relatively stable baseline potential. Based on recent studies of isolated ionic currents, a stochastic Hodgkin-Huxley (HH)-like model for the GnRH neuron is developed to reproduce each mode of burst firing with an appropriate set of conductances. Model outcomes for bursting are in agreement with the experimental recordings in terms of interburst interval, interspike interval, active phase duration, and other quantitative properties specific to each mode of bursting. The model also shows similar outcomes in membrane potential to those seen experimentally when tetrodotoxin (TTX) is used to block action potentials during bursting, and when estradiol transitions cells exhibiting slow oscillations to irregular bursting mode in vitro. Based on the parameter values used to reproduce each mode of bursting, the model suggests that GnRH neurons can switch between the two through changes in the maximum conductance of certain ionic currents, notably the slow inward Ca²⁺ current I _s, and the Ca²⁺ -activated K⁺ current I _KCa. Bifurcation analysis of the model shows that both modes of bursting are similar from a dynamical systems perspective despite differences in burst characteristics.     

Mechanisms of membrane potential oscillation in bursting neurons on the snail,Helix pomatia

• 1.1. The mechanism of generation of membrane potential (MP) oscillations was studied in identified bursting neurons from the snail Helix pomatia.
• 2.2. Long-lasting stimulation of an identified peptidergic interneuron produced a persistent bursting activity in a non-active burster.
• 3.3. External application of calcium channel blockers (1 mM Cd²⁺ or 5 mM La²⁺) resulted in a transient increase in the slow-wave amplitude and subsequent prevention of pacemaker activity generation in bursting neurons. Application of these blockers together with endogenous neuropeptide initiating bursting activity generation, increased MP wave amplitude without prevention of bursting activity generation.
• 4.4. Replacement of all NaCl in normal Ringer's solution with isoosmotic CaCl₂, glucose or Tris-HCl produced a reversible block of bursting activity generation. Stationary current-voltage relation (CVR) of bursting neuron membrane has a region of negative resistance (NRR) and does not intersect the potential axis in threshold region for action potential (AP) generation in normal Ringer's solution. In Na-free solution stationary CVR is linear and intersects the potential axis near — 52 mV.
• 5.5. Novel potential- and time-dependent outward (E_rev = − 58 mV) current, I_B, activated by hyperpolarization was found in the bursting neuron membrane. Having achieved a maximal value, this current decayed with a time constant of about 1 sec. Hyperpolarization inactivated maximal conductance, g_B, responsible for I_B, and depolarization abolished inactivation of g_B.
• 6.6. Short-lasting (0.01 sec) hyperpolarization of the bursting neuron membrane by inward current pulse induced the development of prolonged hyperpolarization wave lasting up to 10 sec.
• 7.7. These results suggest that: (a) persistent bursting activity of RPal neuron in the snail Helix pomatia is not endogenous but is due to a constant activation of peptidergic synaptic inputs of these neurons; (b) Ca²⁺ ions do not play a pivotal role in the ionic mechanism of MP oscillations but play a determining role in the process of secretion of a peptide initiating bursting activity by the interneuron presynaptic terminal; (c) depolarizing phase of the MP wave is due to specific properties of stationary CVR and hyperpolarization phase is due to regenerative properties of hyperpolarization-activated outward current I_B. The minimal mathematical version of MP oscillations based on the experimental data is presented.

     

Convergent neuromodulation onto a network neuron can have divergent effects at the network level

Different neuromodulators often target the same ion channel. When such modulators act on different neuron types, this convergent action can enable a rhythmic network to produce distinct outputs. Less clear are the functional consequences when two neuromodulators influence the same ion channel in the same neuron. We examine the consequences of this seeming redundancy using a mathematical model of the crab gastric mill (chewing) network. This network is activated in vitro by the projection neuron MCN1, which elicits a half-center bursting oscillation between the reciprocally-inhibitory neurons LG and Int1. We focus on two neuropeptides which modulate this network, including a MCN1 neurotransmitter and the hormone crustacean cardioactive peptide (CCAP). Both activate the same voltage-gated current (I _MI ) in the LG neuron. However, I _MI-MCN1 , resulting from MCN1 released neuropeptide, has phasic dynamics in its maximal conductance due to LG presynaptic inhibition of MCN1, while I _MI-CCAP retains the same maximal conductance in both phases of the gastric mill rhythm. Separation of time scales allows us to produce a 2D model from which phase plane analysis shows that, as in the biological system, I _MI-MCN1 and I _MI-CCAP primarily influence the durations of opposing phases of this rhythm. Furthermore, I _MI-MCN1 influences the rhythmic output in a manner similar to the Int1-to-LG synapse, whereas I _MI-CCAP has an influence similar to the LG-to-Int1 synapse. These results show that distinct neuromodulators which target the same voltage-gated ion channel in the same network neuron can nevertheless produce distinct effects at the network level, providing divergent neuromodulator actions on network activity.     

The early topography of thalamocortical projections is shifted in Ebf1 and Dlx1/2 mutant mice

The prevailing model to explain the formation of topographic projections in the nervous system stipulates that this process is governed by information located within the projecting and targeted structures. In mammals, different thalamic nuclei establish highly ordered projections with specific neocortical domains and the mechanisms controlling the initial topography of these projections remain to be characterized. To address this issue, we examined Ebf1(-/-) embryos in which a subset of thalamic axons does not reach the neocortex. We show that the projections that do form between thalamic nuclei and neocortical domains have a shifted topography, in the absence of regionalization defects in the thalamus or neocortex. This shift is first detected inside the basal ganglia, a structure on the path of thalamic axons, and which develops abnormally in Ebf1(-/-) embryos. A similar shift in the topography of thalamocortical axons inside the basal ganglia and neocortex was observed in Dlx1/2(-/-) embryos, which also have an abnormal basal ganglia development. Furthermore, Dlx1 and Dlx2 are not expressed in the dorsal thalamus or in cortical projections neurons. Thus, our study shows that: (1) different thalamic nuclei do not establish projections independently of each other; (2) a shift in thalamocortical topography can occur in the absence of major regionalization defects in the dorsal thalamus and neocortex; and (3) the basal ganglia may contain decision points for thalamic axons' pathfinding and topographic organization. These observations suggest that the topography of thalamocortical projections is not strictly determined by cues located within the neocortex and may be regulated by the relative positioning of thalamic axons inside the basal ganglia.