Analysis of the electroencephalographic activity associated with thalamic tumors

A physiologically based model of corticothalamic dynamics is used to investigate the electroencephalographic (EEG) activity associated with tumors of the thalamus. Tumor activity is modeled by introducing localized two-dimensional spatial non-uniformities into the model parameters, and calculating the resulting activity via the coupling of spatial eigenmodes. The model is able to reproduce various qualitative features typical of waking eyes-closed EEGs in the presence of a thalamic tumor, such as the appearance of abnormal peaks at theta ( approximately 3Hz) and spindle ( approximately 12Hz) frequencies, the attenuation of normal eyes-closed background rhythms, and the onset of epileptic activity, as well as the relatively normal EEGs often observed. The results indicate that the abnormal activity at theta and spindle frequencies arises when a small portion of the brain is forced into an over-inhibited state due to the tumor, in which there is an increase in the firing of (inhibitory) thalamic reticular neurons. The effect is heightened when there is a concurrent decrease in the firing of (excitatory) thalamic relay neurons, which are in any case inhibited by the reticular ones. This is likely due to a decrease in the responsiveness of the peritumoral region to cholinergic inputs from the brainstem, and a corresponding depolarization of thalamic reticular neurons, and hyperpolarization of thalamic relay neurons, similar to the mechanism active during slow-wave sleep. The results indicate that disruption of normal thalamic activity is essential to generate these spectral peaks. Furthermore, the present work indicates that high-voltage and epileptiform EEGs are caused by a tumor-induced local over-excitation of the thalamus, which propagates to the cortex. Experimental findings relating to local over-inhibition and over-excitation are discussed. It is also confirmed that increasing the size of the tumor leads to greater abnormalities in the observable EEG. The usefulness of EEG for localizing the tumor is investigated.     

Thalamic circuitry and thalamocortical synchrony

The corticothalamic system has an important role in synchronizing the activities of thalamic and cortical neurons. Numerically, its synapses dominate the inputs to relay cells and to the gamma-amino butyric acid (GABA)ergic cells of the reticular nucleus (RTN). The capacity of relay neurons to operate in different voltage-dependent functional modes determines that the inputs from the cortex have the capacity directly to excite the relay cells, or indirectly to inhibit them via the RTN, serving to synchronize high- or low-frequency oscillatory activity respectively in the thalamocorticothalamic network. Differences in the alpha-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA) subunit composition of receptors at synapses formed by branches of the same corticothalamic axon in the RTN and dorsal thalamus are an important element in the capacity of the cortex to synchronize low-frequency oscillations in the network. Interactions of focused corticothalamic axons arising from layer VI cortical cells and diffuse corticothalamic axons arising from layer V cortical cells, with the specifically projecting core relay cells and diffusely projecting matrix cells of the dorsal thalamus, form a substrate for synchronization of widespread populations of cortical and thalamic cells during high-frequency oscillations that underlie discrete conscious events.     

Selective GABAergic control of higher-order thalamic relays

GABAergic signaling is central to the function of the thalamus and has been traditionally attributed primarily to the nucleus reticularis thalami (nRT). Here we present a GABAergic pathway, distinct from the nRT, that exerts a powerful inhibitory effect selectively in higher-order thalamic relays of the rat. Axons originating in the anterior pretectal nucleus (APT) innervated the proximal dendrites of relay cells via large GABAergic terminals with multiple release sites. Stimulation of the APT in an in vitro slice preparation revealed a GABA(A) receptor-mediated, monosynaptic IPSC in relay cells. Activation of presumed single APT fibers induced rebound burst firing in relay cells. Different APT neurons recorded in vivo displayed fast bursting, tonic, or rhythmic firing. Our data suggest that selective extrareticular GABAergic control of relay cell activity will result in effective, state-dependent gating of thalamocortical information transfer in higher-order but not in first-order relays.     

Modelling corticothalamic feedback and the gating of the thalamus by the cerebral cortex

Morphological studies have shown that excitatory synapses from the cortex constitute the major source of synapses in the thalamus. However, the effect of these corticothalamic synapses on the function of the thalamus is not well understood because thalamic neurones have complex intrinsic firing properties and interact through multiple types of synaptic receptors. Here we investigate these complex interactions using computational models. We show first, using models of reconstructed thalamic relay neurones, that the effect of corticothalamic synapses on relay cells can be similar to that of afferent synapses, in amplitude, kinetics and timing, although these synapses are located in different regions of the dendrites. This suggests that cortical EPSPs may complement (or predict) the afferent information. Second, using models of reconstructed thalamic reticular neurones, we show that high densities of the low-threshold Ca2+ current in dendrites can give these cells an exquisite sensitivity to cortical EPSPs, but only if their dendrites are hyperpolarized. This property has consequences at the level of thalamic circuits, where corticothalamic EPSPs evoke bursts in reticular neurones and recruit relay cells predominantly through feedforward inhibition. On the other hand, with depolarized dendrites, thalamic reticular neurones do not generate bursts and the cortical influence on relay cells is mostly excitatory. Models therefore suggest that the cortical influence can either promote or antagonize the relay of information, depending on the state of the dendrites of reticular neurones. The control of these dendrites may therefore be a determinant of attentional mechanisms. We also review the effect of corticothalamic feedback at the network level, and show how the cortical control over the thalamus is essential in co-ordinating widespread, coherent oscillations. We suggest mechanisms by which different modes of corticothalamic interaction would allow oscillations of very different spatiotemporal coherence to coexist in the thalamocortical system.     

Synchronized oscillations at alpha and theta frequencies in the lateral geniculate nucleus

In relaxed wakefulness, the EEG exhibits robust rhythms in the alpha band (8-13 Hz), which decelerate to theta (approximately 2-7 Hz) frequencies during early sleep. In animal models, these rhythms occur coherently with synchronized activity in the thalamus. However, the mechanisms of this thalamic activity are unknown. Here we show that, in slices of the lateral geniculate nucleus maintained in vitro, activation of the metabotropic glutamate receptor (mGluR) mGluR1a induces synchronized oscillations at alpha and theta frequencies that share similarities with thalamic alpha and theta rhythms recorded in vivo. These in vitro oscillations are driven by an unusual form of burst firing that is present in a subset of thalamocortical neurons and are synchronized by gap junctions. We propose that mGluR1a-induced oscillations are a potential mechanism whereby the thalamus promotes EEG alpha and theta rhythms in the intact brain.     

The initiation of bursts in thalamic neurons and the cortical control of thalamic sensitivity

Thalamic neurons generate high-frequency bursts of action potentials when a low-threshold (T-type) calcium current, located in soma and dendrites, becomes activated. Computational models were used to investigate the bursting properties of thalamic relay and reticular neurons. These two types of thalamic cells differ fundamentally in their ability to generate bursts following either excitatory or inhibitory events. Bursts generated with excitatory inputs in relay cells required a high degree of convergence from excitatory inputs, whereas moderate excitation drove burst discharges in reticular neurons from hyperpolarized levels. The opposite holds for inhibitory rebound bursts, which are more difficult to evoke in reticular neurons than in relay cells. The differences between the reticular neurons and thalamocortical neurons were due to different kinetics of the T-current, different electrotonic properties and different distribution patterns of the T-current in the two cell types. These properties enable the cortex to control the sensitivity of the thalamus to inputs and are also important for understanding states such as absence seizures.     

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.     

Relating Alpha Power and Phase to Population Firing and Hemodynamic Activity Using a Thalamo-cortical Neural Mass Model

Oscillations are ubiquitous phenomena in the animal and human brain. Among them, the alpha rhythm in human EEG is one of the most prominent examples. However, its precise mechanisms of generation are still poorly understood. It was mainly this lack of knowledge that motivated a number of simultaneous electroencephalography (EEG) – functional magnetic resonance imaging (fMRI) studies. This approach revealed how oscillatory neuronal signatures such as the alpha rhythm are paralleled by changes of the blood oxygenation level dependent (BOLD) signal. Several such studies revealed a negative correlation between the alpha rhythm and the hemodynamic BOLD signal in visual cortex and a positive correlation in the thalamus. In this study we explore the potential generative mechanisms that lead to those observations. We use a bursting capable Stefanescu-Jirsa 3D (SJ3D) neural-mass model that reproduces a wide repertoire of prominent features of local neuronal-population dynamics. We construct a thalamo-cortical network of coupled SJ3D nodes considering excitatory and inhibitory directed connections. The model suggests that an inverse correlation between cortical multi-unit activity, i.e. the firing of neuronal populations, and narrow band local field potential oscillations in the alpha band underlies the empirically observed negative correlation between alpha-rhythm power and fMRI signal in visual cortex. Furthermore the model suggests that the interplay between tonic and bursting mode in thalamus and cortex is critical for this relation. This demonstrates how biophysically meaningful modelling can generate precise and testable hypotheses about the underpinnings of large-scale neuroimaging signals.     

Cholinergic mechanisms in the pathogenesis of genetically-caused absence epilepsy

Frontoparietal cortex and the thalamocortical circuit comprising reticular thalamic nucleus (RTN) and relay nuclei of the ventrolateral thalamus (VLT) are critical structures in the generation of spike-wave discharges (SWD) during absence seizures. The activity of these nuclei is under the control of the ascending cholinergic projections of nucleus basalis of Meynert. The aim of our study is to make an attempt to change the pattern of SWD in WAG/Rij rats by injecting of cholinotoxine AF64A to the area of RTN. Spontaneous SWD were registered in cortex of WAG/Rij rats with genetically determined absences. The spectral content of SWD was analyzed by means of the Fast Fourier Transformation (FFT) procedure. Unilateral injections of AF64A (1 nmol) to RTN led the decrease in duration and number of SWD comparing to the basal EEG recordings 2 days after the lesion. The FFT analysis showed the disappearance of 17-18 Hz spike on the side of the lesion compared with the intact side. The immunohistochemical study for acetylcholinetransferase (ChaT)-containing neurons showed the loss of ChaT-positive cells in the nucleus basalis area on the side of the lesion. The removal of cholinergic afferentation of RTN and cortex from nucleus basalis inhibits the SWD developing most likely due to the decrease of cortical excitability. Moreover, possibly cholinergic transmission is involved in the transforation of the synchronized phenomena (SWD) to another with close mechanism of generation.     

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     

Simulation of GABA B -receptor-mediated K+ current in thalamocortical relay neurons: tonic firing,bursting, and oscillations

Until recently, the presence of γ-aminobutyric acid (GABA) in the thalamus has usually been associated with the "classical" GABA _A Cl⁻-dependent receptor. However, the discovery of a slower, long-lasting, K⁺-dependent inhibitory postsynaptic potential (IPSP) mediated by GABA _B receptors in projection cells of the dorsal lateral geniculate nucleus has led researchers to reconsider its role in modulating the behavior of these cell groups (Crunelli et al. 1988; Crunelli and Leresche 1991). Of particular interest is the role of this K⁺ current in the activation of the low-threshold Ca²⁺ current, I _T , of thalamocortical relay (TCR) neurons responsible for bursting activity (Jahnsen and Llinás 1984 a, b). Considering the time scale on which the GABA _B -receptor-activated K⁺ current operates, it is ideally suited to foster sustained rhythmicity in TCR cells reciprocally connected to neurons of the nucleus reticularis thalami (NRT) as well as interneurons at frequencies observed in vivo (Steriade and Llinás 1988). In this study we show that small changes in the duration and amplitude of the K⁺-dependent IPSPs can have marked effects on TCR cell groups including a shift from single-spike firing (tonic) to bursting behavior. We further show that a single GABA _B -mediated IPSP is sufficient to activate the low-threshold Ca²⁺ response and that sustained oscillations are possible given the presence of excitatory TCR connections to GABAergic NRT cells or interneurons of the dorsal lateral thalamus. These combined effects are examined with regard to their role in generating the well known 7 – 14 Hz spindle rhythm as well as slower 6 – 8 Hz oscillations observed in TCR cells in vivo (Steriade and Llinás 1988). Received: 13 October 1993 / Accepted in revised form: 24 February 1994     

The site of RanGTP generation can act as an organizational cue for mitotic microtubules

Background information. RanGTP, which is generated on chromosomes during mitosis, is required for microtubule spindle assembly. Due to its restricted spatial generation within the cell it has been suggested that RanGTP acts as a spatial cue to organize site‐specific spindle assembly within the cell. However, the absence of a detectable sharp gradient of RanGTP in somatic cells has led to suggestions that it may only act as a spatial cue in large cells and that it may operate as a general activator of the mitotic cytosol in somatic cells. Results. We report that ectopic generation of RanGTP at the plasma membrane stimulates the formation of organized arrays of microtubules at the plasma membrane. Conclusions. These results suggest that the site of RanGTP generation in a mitotic somatic cell can generate critical spatial information that specifies where microtubules grow towards and where microtubules are organized. As RanGTP is normally generated on chromosomes, these results suggest that RanGTP may play an important role in specifying that spindle assembly occurs around chromosomes.     

Stimulus-induced transitions between spike-wave discharges and spindles with the modulation of thalamic reticular nucleus

It is believed that thalamic reticular nucleus (TRN) controls spindles and spike-wave discharges (SWD) in seizure or sleeping processes. The dynamical mechanisms of spatiotemporal evolutions between these two types of activity, however, are not well understood. In light of this, we first use a single-compartment thalamocortical neural field model to investigate the effects of TRN on occurrence of SWD and its transition. Results show that the increasing inhibition from TRN to specific relay nuclei (SRN) can lead to the transition of system from SWD to slow-wave oscillation. Specially, it is shown that stimulations applied in the cortical neuronal populations can also initiate the SWD and slow-wave oscillation from the resting states under the typical inhibitory intensity from TRN to SRN. Then, we expand into a 3-compartment coupled thalamocortical model network in linear and circular structures, respectively, to explore the spatiotemporal evolutions of wave states in different compartments. The main results are: (i) for the open-ended model network, SWD induced by stimulus in the first compartment can be transformed into sleep-like slow UP-DOWN and spindle states as it propagates into the downstream compartments; (ii) for the close-ended model network, weak stimulations performed in the first compartment can result in the consistent experimentally observed spindle oscillations in all three compartments; in contrast, stronger periodic single-pulse stimulations applied in the first compartment can induce periodic transitions between SWD and spindle oscillations. Detailed investigations reveal that multi-attractor coexistence mechanism composed of SWD, spindles and background state underlies these state evolutions. What’s more, in order to demonstrate the state evolution stability with respect to the topological structures of neural network, we further expand the 3-compartment coupled network into 10-compartment coupled one, with linear and circular structures, and nearest-neighbor (NN) coupled network as well as its realization of small-world (SW) topology via random rewiring, respectively. Interestingly, for the cases of linear and circular connetivities, qualitatively similar results were obtained in addition to the more irregularity of firings. However, SWD can be eventually transformed into the consistent low-amplitude oscillations for both NN and SW networks. In particular, SWD evolves into the slow spindling oscillations and background tonic oscillations within the NN and SW network, respectively. Our modeling and simulation studies highlight the effect of network topology in the evolutions of SWD and spindling oscillations, which provides new insights into the mechanisms of cortical seizures development.     

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.     

Modeling the contribution of lamina 5 neuronal and network dynamics to low frequency EEG phenomena

The Electroencephalogram (EEG) is an important clinical and research tool in neurophysiology. With the advent of recording techniques, new evidence is emerging on the neuronal populations and wiring in the neocortex. A main challenge is to relate the EEG generation mechanisms to the underlying circuitry of the neocortex. In this paper, we look at the principal intrinsic properties of neocortical cells in layer 5 and their network behavior in simplified simulation models to explain the emergence of several important EEG phenomena such as the alpha rhythms, slow-wave sleep oscillations, and a form of cortical seizure. The models also predict the ability of layer 5 cells to produce a resonance-like neuronal recruitment known as the augmenting response. While previous models point to deeper brain structures, such as the thalamus, as the origin of many EEG rhythms (spindles), the current model suggests that the cortical circuitry itself has intrinsic oscillatory dynamics which could account for a wide variety of EEG phenomena.Electronic Supplementary Material Supplementary material is available for this article at Sensorimotor Control Project- MIT Harvard NeuroEngineering Research Collaborative.     

The role of the thalamus in the flow of information to the cortex

The lateral geniculate nucleus is the best understood thalamic relay and serves as a model for all thalamic relays. Only 5-10% of the input to geniculate relay cells derives from the retina, which is the driving input. The rest is modulatory and derives from local inhibitory inputs, descending inputs from layer 6 of the visual cortex, and ascending inputs from the brainstem. These modulatory inputs control many features of retinogeniculate transmission. One such feature is the response mode, burst or tonic, of relay cells, which relates to the attentional demands at the moment. This response mode depends on membrane potential, which is controlled effectively by the modulator inputs. The lateral geniculate nucleus is a first-order relay, because it relays subcortical (i.e. retinal) information to the cortex for the first time. By contrast, the other main thalamic relay of visual information, the pulvinar region, is largely a higher-order relay, since much of it relays information from layer 5 of one cortical area to another. All thalamic relays receive a layer-6 modulatory input from cortex, but higher-order relays in addition receive a layer-5 driver input. Corticocortical processing may involve these corticothalamocortical 're-entry' routes to a far greater extent than previously appreciated. If so, the thalamus sits at an indispensable position for the modulation of messages involved in corticocortical processing.     

Neocortical Dynamics: Implications for Understanding the Role of Neurofeedback and Related Techniques for the Enhancement of Attention

For nearly 25 years, EEG biofeedback (neurofeedback) has been utilized in research and clinical settings for the treatment and investigation of a number of disorders ranging from attention deficit hyperactivity disorder to seizure disorders as well as many other established and investigational applications. Until recently, mechanisms underlying the generation and origins of EEG have been poorly understood but now are beginning to become much more clarified. Now it is important to combine the information gathered on the genesis of EEG and neocortical dynamics with the findings from neurofeedback investigations. This will help us to develop models of how neurofeedback might operate in producing the changes in EEG and in clinical symptomatology. We know that the cortex operates in terms of resonant loops between neocortical columns of cells known as local, regional, and global resonances. These resonances determine the specific EEG frequencies and are often activated by groups of cells in the thalamus known as pacemakers. There are complex excitatory and inhibitory interactions within the cortex and between the cortex and the thalamus that allow these loops to operate and provide the basis for learning. Neurofeedback is a technique for modifying these resonant loops, and hence, modifying the neurophysiological and neurological basis for learning and for the management of a number of neurologically based disorders. This paper provides an introduction to understanding EEG and neocortical dynamics and how these concepts can be used to explain the results of neurofeedback training and other interventions particularly in the context of understanding attentive mechanisms and for the management of attention deficit/hyperactivity disorders.     

基于生理解剖知识的入睡机制神经元群网络模型研究

以生理解剖知识为基础,在已有的丘脑网状核细胞和丘脑皮质细胞间组成的入睡机制的两细胞环路模型［１］和由此两细胞环路组成的网络模型［２］的基础上,提出了增加皮层细胞在内的三种细胞组成的环路模型和网络模型,以使模型更符合近来认为睡眠机制是皮层和丘脑环路中出现特定的同步振荡的看法［３］。并能使模型的仿真结果可以和规定人体睡眠分期的脑电特征波相对应。这一网络模型的仿真结果,在一定条件下,确能在皮层细胞处出现符合睡眠分期中规定的标志入睡的纺锤波,这一初步结果,启示我们用模型仿真方法来进一步探讨睡眠机制和用模型仿真方法来进一步探讨人脑的微观神经元的电活动是如何通过同步振荡整合到宏观功能状态的某些信息处理过程的可能性。     

Correlations of neuronal spike discharges of VL neurons during spontaneous firing and during the activity evoked by peripheral stimulation

The temporal relations between simultaneously recorded neurons of the nucleus ventralis lateralis (VL) of cat thalamus were studied. The interaction and the functional connections between individual VL neurons are described. This was achieved with an application of cross correlation techniques. The response patterns of different individual neurons to somatic sensory and photic stimuli were also analyzed. For the purpose of classifying neurons as thalamocortical relay cells (T-C) and non relay cells (N-C) which do not project to the motor sensory cortex antidromic cortical stimulation was used. This stimulation was also used as conditioning one when proceeded peripheral stimuli. To analyze the nonspecific specific interactions upon single neurons conditioning photic stimuli were applied. The results show that T-C neurons are antidromically excited from a wide cortical areas and that the functional interaction between T-C neurons is mediated by a shared input from common sources. It is further postulated that N-C cells interposed between relay neurons subserve the functions of gating units modifying the neuronal network of lateral ventral nucleus of the thalamus.     

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