Sleep waves in a large‐scale corticothalamic model constrained by activities intrinsic to neocortical networks and single thalamic neurons

Abstract Aim Many biophysical and non‐biophysical models have been able to reproduce the corticothalamic activities underlying different EEG sleep rhythms but none of them included the known ability of neocortical networks and single thalamic neurons to generate some of these waves intrinsically. Methods We built a large‐scale corticothalamic model with a high fidelity in anatomical connectivity consisting of a single cortical column and first‐ and higher‐order thalamic nuclei. The model is constrained by different neocortical excitatory and inhibitory neuronal populations eliciting slow (<1 Hz) oscillations and by thalamic neurons generating sleep waves when isolated from the neocortex. Results Our model faithfully reproduces all EEG sleep waves and the transition from a desynchronized EEG to spindles, slow (<1 Hz) oscillations, and delta waves by progressively increasing neuronal membrane hyperpolarization as it occurs in the intact brain. Moreover, our model shows that slow (<1 Hz) waves most often start in a small assembly of thalamocortical neurons though they can also originate in cortical layer 5. Moreover, the input of thalamocortical neurons increases the frequency of EEG slow (<1 Hz) waves compared to those generated by isolated cortical networks. Conclusion Our simulations challenge current mechanistic understanding of the temporal dynamics of sleep wave generation and suggest testable predictions.

Using current knowledge of anatomical connectivity, [27][28][29][30][31][32][33][34][35][36][37] we built a large-scale corticothalamic model that contains a single cortical column (with fast spiking, FS, inhibitory interneurons, and different excitatory neuronal populations) that is reciprocally connected to first-order and higher-order thalamic nuclei and its respective γ-aminobutyric acid (GABA)ergic neurons of the nucleus reticularis thalami (NRT).[25] Our corticothalamic model accurately reproduces all sleep EEG rhythms and the transition from a desynchronized EEG to various sleep rhythms.Our results show that the EEG slow (<1 Hz) waves are most often initiated by a small assembly of thalamocortical neurons though they can also originate in cortical layer 5 (L5) and that their frequency is decreased in the absence of thalamic afferent activity.Moreover, our simulations suggest testable predictions.

| Model network connectivity
The model network connectivity and the connection weights are provided in (Figure S1 and Tables S1,S2).Connections were organized topographically with sources and targets located in matching regions of their corresponding structures.A neuron did not synapse onto itself and could only form a single synapse on its target neuron.The number of contacts that a source neuron could form in a target structure was defined by the parameter P (a projection radius) (Table S1).
Other key connectivity parameters (e.g., connection weight, postsynaptic potential shape, synaptic transmission latency, synaptic receptors) are listed in Table S2 and Supplementary Methods.
We made several simplifying assumptions.In certain instances, we scaled down differences in anatomical projection radii.This was true for intra-cortical connections simplifying network calibration.
Moreover, we scaled down the physiological N-methyl D-aspartate (NMDA) component in cortical synapses (i.e., 12-14 times smaller rather than three to five times as the AMPA component) in order to eliminate a slow periodicity component in paroxysmal simulations which had not been observed experimentally.NRT cells were connected by gap junctions, 42 that is, the first-and second-degree neighbors formed single 3 GΩ and 4.5 GΩ junctions, respectively.

| Model neurons
TC and NRT neurons were single-compartment Hodgkin-Huxley models whereas cortical neurons were Hodgkin-Huxley models with separate axosomatic and dendritic compartments.The equations for both thalamic and cortical neuron models are listed in the Supplementary Methods, that also describe the intrinsic and synaptic membrane currents and their dependence on intracellular ion concentration (where applicable) (Tables S3-S7).

| Simulations
All simulations were carried out in NEURON 43

| Data analysis
Simulation data were analyzed and visualized with the help of custom-written Matlab (MathWorks Inc routines).The scalp EEG signal produced by the simulations was estimated as described by Bédard et al (2004)  45 (see Supplementary Methods).

| Slow (<1 Hz) oscillations in the isolated neocortical network
We were able to reproduce the typical firing patterns observed at different levels of membrane polarization in RS and EF neurons, the pyramidal cell type that shows early firing during slow (<1 Hz) oscillations 22 (Figure S2), as well as in IB and FS neurons 22 (Figure S3).
We also developed models of RIB neurons and of the ND pyramidal neurons that "drive" slow (<1 Hz) oscillations in neocortical slices. 22e firing patterns of RIB and ND neurons well matched those observed experimentally 22 (Figure S4A), including their input-output curves (Figure S4B).When all these neocortical excitatory and inhibitory model neurons were synaptically connected in a cortical column (Figure S1), a pattern of non-rhythmic, low-amplitude activity was observed in the EEG for low K + leak conductance (g KL ) values (Table S8) in all component neurons (Figure 1L,M).When g KL was increased (Table S8), however, the neocortical network generated clear EEG slow (<1 Hz) waves (Figure 1C-E) that could also give rise to faster waves within the delta frequency band (Figure 1F-H).Thus, all neurons in all cortical layers oscillated at a frequency of 0.44 Hz (Figure 1C-E), that is, within the slow (<1 Hz) oscillation frequency, whereas at lower g KL values the oscillation reached the frequency of 1.57 Hz, that is, within the δ frequency range (Figure 1F-H): indeed, as g KL was changed the neocortical oscillation frequency smoothly changed accordingly (not shown).Finally, further small progressive decreases in g KL in all cortical cells (except EF neurons) (Table S8) resulted in the break-down of the regular pattern of oscillations (Figure 1I-K), and eventually led to a desynchronized, low-amplitude EEG similar to that of the wake state (Figure 1L,M) where single neurons fired randomly.
The slow (<1 Hz) oscillation elicited by the isolated neocortical network (Figure 1C-E) was initiated in L5 pyramidal neurons since the firing appeared there earlier than in other layers, followed by pyramidal neurons in L6, L2/3, and L4 (Figure S5A,B).In each layer, the inhibitory neuron firing occurred after that of the pyramidal neurons with little difference among layers (Figure S5A,B).The ND neurons in L5 are the cells that initiate the oscillation, followed by the EF neurons in L/6 and L2/3 (Figure S5C).Thus, the onset timing of the simulated neocortical slow (<1 Hz) oscillation is highly consistent with experimental data obtained in cortical slices. 22,38Indeed, the duration of the up-and down-states of the simulated slow (<1 Hz) oscillations in the isolated neocortex was similar to that observed in vitro (~250 ms and ~840 ms, respectively) 22 (Figure S5D,E).
The two main contributors to the maintenance of the up-state were, in order of importance, the persistent Na + current (I Na(P) ) and the excitatory postsynaptic potential (EPSP) barrage in IB (Figure S6G).Blocking I K[Ca] in all model neurons except ND increased the up-state duration confirming its causal role in the termination (Figure S7I).Similarly, block of I M increased the duration and irregularity of the up-states confirming its role in curtailing upstates (Figure S7J).On the other hand, the block of I K[Na] resulted in an almost-continuous up-state, that is, the up-state never properly terminated (Figure S7K) since EF neurons fired continuously (not shown).
The contribution to the initiation, maintenance, and termination of up-states by the T-type Ca 2+ current (I T ), the high-threshold Ca 2+ current (I HVA ), and the hyperpolarization-activated, cyclic nucleotidegated current (I h ) was negligible in IB and RS model cells (Figure S6H).Blocking I T , I HVA , and I h (except for ND) had no significant effect on the duration or shape of the up-states (Figure S7L-N).However, I h and I HVA , as the key pacemaker currents, initiated the up-states in ND neurons (not shown).
Blocking GABA A Rs resulted in a continuous paroxysmal oscillation (Figure S7D) as it was observed experimentally in cortical slices after application of bicuculline. 48Hence, the effect of GABA A Rs is to reduce the intensity of cortical up-states and, in turn, to increase their duration by limiting the build-up of I K[Ca] and I K[Na] , an observation consistent with the experimental data. 49The role of GABA B Rs was to shorten the up-states since their blockade resulted in upstates longer than 1 s (Figure S6E), an effect opposite to that of GABA A R and consistent with in vitro experiments. 49ving established the dynamics of initiation and maintenance of the EEG slow (<1 Hz) oscillation in the isolated cortical network, we next investigated whether the membrane potential waveforms of the different cortical populations during the up-and down-states dynamics were similar to those observed in neocortical slices. 22,38e excitatory RS, EF, and ND neurons faithfully reproduced the up-and down-state dynamics, the firing patterns, and the bimodal distribution of the membrane potential characteristic of neocortical slow (<1 Hz) oscillations in vitro 22 (Figure 2A-H, M-P).In particular, ND neurons exhibited the diverse activity observed in this pyramidal neuron type at an early and late stages of the slow (<1 Hz) oscillation as shown following its pharmacological induction in cortical slices 22 (Figure 2O,P).Moreover, the membrane potential waveform of IB and RIB pyramidal neuron models during the simulated neocortical slow (<1 Hz) oscillations well matched that reported in cortical slices 22 (Figure S8).A slow (<1 Hz) oscillation was also present in the inhibitory FS neurons, although a clear bimodal distribution was evident in their firing but not in their membrane potential distribution plots as observed experimentally 22 (Figure 2I-L).Finally, a faithful reproduction of the slow (<1 Hz) oscillation and its membrane potential distribution plots was also found for IB and RIB neurons 22 (Figure S8).We then simulated the intrinsic generation of delta and slow (<1 Hz) oscillations observed in TC FO neurons in vitro following pharmacological activation of metabotropic glutamate receptors 1a (mGluR1a), [23][24][25] which was mimicked by reducing g KL (Figure 3).The TC FO neurons closely reproduced the hyperpolarization-dependent transition from tonic firing to slow (<1 Hz) oscillations at increasing frequency (due to shorter up states) and then to delta waves (Figure 3B,D  NRT FO , first-order NRT neurons; NRT HO , higher-order NRT neurons; TC FO , first-order TC neurons; TC HO , higher-order TC neurons.

| Sleep waves in the full corticothalamic model
Experimentally it was shown that whereas the neocortex can generate the slow (<1 Hz) oscillation when isolated from the thalamus, its up-states are less rhythmic and frequent in the absence of the thalamus, 52 a finding that was faithfully reproduced by our full corticothalamic model (Figure S12).Connecting cortex and thalamus produced stronger up-state firing in both structures: as a result, K +   4B), the up-state-linked firing appeared initially in a small number of TC FO neurons that, notably, was different from one up-state to the next (Figure 4B,C).This was followed by firing in NRT FO and NRT HO neurons and then in L4, L5, L6, and L2/3, with the start of the EEG slow wave coinciding with the start of the L4 neuron firing (Figure 4C).Notably, when ND cells were depolarized, up-states could also start in L5 (Figure S13).[25] The membrane potential waveforms of different neuronal types during the slow (<1 Hz) oscillation were similar to those observed experimentally in vivo 1,46,52 (Figure 4B).The most active up-states (i.e., containing APs) among the cortical neurons were in ND, EF, and FS neurons, whereas RS, IB, and RIB displayed on average a sparser firing (Figure 4B).The firing in thalamic neurons was more regular than in cortex and had a clear intrinsic aspect to it (Figure 4B).The onset of regularly shaped up-states in all thalamic neurons followed a down-state of an almost fixed length and invariably terminated with an I T -mediated burst of action potentials (Figures 3D and 4A), as observed in thalamic neurons in vitro [23][24][25][26] : indeed, thalamic activity was essentially an intrinsic oscillation that was shaped by, and synchronized with, the neocortex, a result supported by experimental findings. 26,53tably, when g KL of TC and NRT neurons was increased (Table S9) delta waves started to occur during the down-states of the slow (<1 Hz) oscillations of these thalamic neurons, as observed experimentally [23][24][25] : their reflection in the cortical territory resulted in a speeding up of the slow (<1 Hz) oscillation into the delta frequency range (Figure 5B) as seen in the EEG (Figure 5A).Whereas during slow (<1 Hz) waves firing was observed first in a small number of TC FO neurons, during delta waves there was little time lag across the entire TC FO neuron population and the delay between thalamic and L4 neuron firing in each cycle was much smaller (<100 msec) Finally, by adjusting g KL and increasing the fast component of K + after hyperpolarizing current (I AHP ) in NRT cells the corticothalamic model was able to reproduce EEG sleep spindles (Figures 6 and S14).
Spindle waves appeared regularly in the EEG and had the "classical" waveform observed in vivo 1,46 (Figures 6A,B and S14A,B).During a spindle wave, firing started in NRT FO neurons, followed in a few milliseconds by TC HO and NRT HO neurons, and then TC FO and cortical neurons (Figure 6B-E).Notably, the first firing in TC FO neurons resulted from IPSPs summation leading to a rebound burst, partly explaining the firing initiation in these thalamic neurons (Figure 6B,E).
However, the first firing of a spindle cycle could also originate in TC HO neurons (Figure S14).Notably, in both cases it is evident that the spindle wave builds up in the thalamus before being reflected into the cortical territory (Figures 6C and S14B).

| DISCUSS ION
The main finding of this study, that used a bottom-up modeling strategy to constrain the corticothalamic model with parameters obtained from in vivo and in vitro experiments, is the ability to replicate in the EEG and in cortical and thalamic neuronal populations the slow (<1 Hz), delta, and spindle waves and the transition between these

| Model limitations
The ND neuron model has limitations, including the down-states being more depolarized and having a larger variance due to a pronounced I h -mediated sag than ND neurons recorded in vitro. 22It may be possible that I T is the main pacemaker in the real ND neurons as opposed to I HVA in the model cells.In the absence of experimental data, initial attempts to model ND neurons used the I T of TC neurons.However, either because I T parameters are not appropriate for the cortical axosomatic and dendritic compartments or the model ND neuron morphology was inadequate, I T did not support the intrinsic slow (<1 Hz) oscillation.Having I T as the actual pacemaker current and not I HVA could drastically stabilize slow (<1 Hz) oscillations by increasing the membrane potential difference between upand down-states in ND neurons and by increasing the AP frequency within a burst leading to more densely packed EPSPs in the postsynaptic neurons and thus helping the transition to a new up-state.
The frequency of slow/delta intrinsic oscillations of ND cells increases as they are more depolarized, a result similar to that observed in vitro 22 This, in turn, is reflected in the frequency of modeled cortical network slow/delta oscillations, although this is the opposite of the results in cortical slices. 22As indicated above, this is most likely the consequence of the lack of any data on the relative contribution of different intrinsic currents to the excitability of these neurons.
The model does not include the tight regulation of the intrinsic oscillations of thalamic neurons by changes in the function of their metabotropic glutamate receptors as shown experimentally, [23][24][25] that is, the ability of these neurons to behave as conditional oscillators. 26This could be achieved by linking activation of these modulatory receptors to the firing of the thalamic projecting L5-6 axons, as shown with electrical stimulation of the corticothalamic afferents in thalamic slices. 54

| Model strengths
Up-states in the cortical network model were initiated by ND and EF neurons, and, to the best of our knowledge, this study is the first that successfully simulate EEG slow (<1 Hz) waves with these intrinsic mechanisms.However, having these intrinsic initiation processes does not exclude the involvement of other mechanisms of initiation and maintenance of this sleep rhythm, in particular the accumulation of spontaneous synaptic activity 55 and the synchronization of already active neuronal assemblies. 56Notably, however, in our model simulated excitatory synaptic conductances were slightly higher than the inhibitory ones, in contrast to experimental studies that reported the opposite 57 or an equal influence by the two, 58 although these studies were performed in vitro.Thus, our model confirms and enlarges our previous simulations obtained in a thalamocortical model with a limited number of cortical neuron types and fewer thalamic neuron conductances. 19Notably, our simulations are also in agreement with the experimental data of David et al. ( 2013) 51 as they show that connecting thalamus to cortex increases the frequency and rhythmicity of EEG slow (<1 Hz) waves.
Whereas in the isolated neocortex, the appearance of the slow (<1 Hz) oscillations matches the temporal profile reported in cortical slices, 22,38 that is, L5, and L4 and L6, in the full corticothalamic model the first cortical firing is preferentially observed in L4 followed by L5 and then L6.This is not surprising because of the prevailing input to L4 from the TC FO neurons that are the first to start firing in the full corticothalamic model.However, there are cases in the corticothalamic model where firing starts in L5, as in neocortical slices. 22,38us, whereas the full corticothalamic model maintains the inherent on a desktop computer or one of the following computing clusters: the Neuroscience Gateway (NSG) Portal for Computational Neuroscience 44 or the Cardiff School of Biosciences Biocomputing Hub HPC/Cloud infrastructure.

F I G U R E 1
Slow (<1 Hz) and delta oscillations in the isolated neocortical model.EEG (top trace in A, C, F, I, L) and corresponding colorcoded membrane potential graphs for the indicated neuronal populations (bottom traces in A, C, F, I, L), EEG autocorrelograms (B, D, G, J, M) and EEG power graphs (E, H, K) of the indicated activity.(A, B), with a large g KL the cortical model does not elicit any activity.(C-E), Decreasing g KL of ND and EF cells results in spontaneous slow (<1 Hz) oscillations in the entire cortical network.(F-H), Further reducing g KL of ND cells leads to delta waves in the isolated cortical model.(I-K), When g KL of other cortical neurons (except ND and EF neurons) is slightly decreased a non-regular activity pattern is observed, that is, there is a break-down of the oscillations.(L, M), The network shows a desynchronized EEG and tonic firing in all neuronal populations when g KL of all cell types (except ND cells) is further decreased.L4 PY, pyramidal neurons in cortical layer 4; L4 IN, interneurons in cortical layer 4; L5 PY, pyramidal neurons in cortical layer 5; L5 IN, interneurons in cortical layer 5; L6 PY, pyramidal neurons in cortical layer 6; L6 IN, interneurons in cortical layer 6; L2/3 PY, pyramidal neurons in cortical layers 2 and 3; L2/3 IN, interneurons in cortical layers 2 and 3.
) and RS (not shown) neurons.Reducing α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA)-mediated EPSPs or I Na(P) abolished up-states (FigureS7B,F) confirming their causal roles in the initiation and maintenance of the slow (<1 Hz) oscillation, whereas blockade of NMDA receptors had little impact on up-states, indicating that NMDA receptor function is not necessary for the maintenance of cortical up-states (FigureS7C).Other significant contributors in maintaining the up-and downstate dynamics were the A-type K + current (I A ) and the inhibitory postsynaptic potential (IPSP) barrage (FiguresS6F and S7G).Their blockade resulted in over-excitation and the transformation of the slow (<1 Hz) oscillation into a paroxysmal oscillation (Figures S7D) (implying a balancing role for GABA A receptors) or in the disappearance of the down-states as the remaining K + currents could no longer terminate up-states (FigureS7G) (implying a role for I A in the termination of the slow (<1 Hz) oscillation).The total intrinsic hyperpolarizing and depolarizing currents were similar as were the excitatory and inhibitory synaptic currents (FigureS6D,E,I).The termination of the up-state occurred because of the gradual accumulation of the Ca 2+ -activated K + current (I K[Ca] ) but also the M-type K + current (I M ) and the Na + -activated K + current (I K[Na] ) Our model TC FO , TC HO , NRT FO , and NRT HO neurons reproduced the 'classical' activity patterns of these cells, that is, tonic single action potential (AP) and T-type Ca 2+ channel-mediated burst firing, in particular the typical burst signature of decelerando and F I G U R E 2 Experimental and simulated membrane potential dynamics in cortical neurons during slow (<1 Hz) oscillations in the isolated cortical network.(A, B), Simultaneous local field potential and membrane potential dynamics of an RS neuron recorded in vitro (experiment) during slow (<1 Hz) oscillations and its simulated activity in the isolated cortical network (simulation).The left-hand traces show the initial stage of the oscillation while the right-hand traces show the oscillation at a later stage (increased neuromodulatory drive).(C, D), Normalized membrane potential distribution plots for the experimental and simulated slow (<1 Hz) oscillations of an RS neuron with the typical peaks of the up-and down-states.(E, F), Same as (C, D) but for an EF neuron.(G, H), Same as (A, B) for an EF neuron.(I, J), Same as (A, B) but for an FS neuron.(K, L), Same as (C, D) but for an FS neuron.Note the lack of a clear bimodal membrane potential distribution in the plots of both the experimental and simulated data of the FS neurons.(M, N), Same as (C, D) but for an ND neuron.(O, P), Same as (A, B) but for an ND neuron.Experimental data are reproduced with permission from Lorincz et al. (2015). 22accelerando-decelerando burst firing patterns of TC and NRT neurons, 50 respectively (not shown).
).Moreover, the full spectrum of delta and slow (<1 Hz) oscillations observed at different levels of membrane polarization were elicited by the NRT FO neuron model (Figures S9).

25 F I G U R E 4
Notably, I Na(P)and the Ca 2+ -activated non-specific cation current (I CAN ) played a major role in the duration of the up-state of NRT FO neurons since their reduction led to a progressive shortening, and eventually to the abolition of, the up-states, that is, the slow (<1 Hz) oscillation was gradually transformed into a delta wave (FigureS10B1-C2).Notably, NRT FO neurons were capable of eliciting both slow and fast sleep spindle waves (6.5-16 Hz) in isolation, that is, when disconnected from TC FO neurons (FigureS11).Finally, TC HO and NRT HO neurons reproduced slow (<1 Hz) and delta oscillations similar to those of TC FO and NRT FO neurons, respectively (not shown).F I G U R E 3 Experimental and simulated membrane potential dynamics of TC FO neurons during intrinsically generated slow (<1 Hz) oscillations.(A), Tonic firing (top trace), quiescence (middle trace), and delta oscillations (bottom trace) are observed at different level of membrane polarization (injected current at the bottom of each trace) in a cat ventrobasal TC neuron recorded in a thalamic slice maintained in a standard recording solution.(B), After application of 100 μM trans-ACPD (a metabotropic glutamate receptor, mGluR, agonist) the same neuron shows slow (<1 Hz) oscillations between the quiescent state and the delta oscillations (bottom trace).(C), Simulation in the TC FO neuron model of the experimental activity shown in (A).(D), Simulation of experimental activity shown in (B) was obtained after reducing g KL of the TC FO neuron model to mimic the effect of mGluR activation in vitro.Dashed lines on the left of each trace indicate −60 mV.(A) and (B) are reproduced with permission from Zhu et al. (2006).Slow (<1 Hz) oscillations in the full corticothalamic model.(A), EEG showing the rhythmic pattern of slow (<1 Hz) oscillations.(B), EEG (top trace) and color-coded membrane potential plots of the indicated cortical and thalamic neuronal populations during the five cycles of the slow (<1 Hz) oscillation highlighted in (A) (note the two separate color-scales for the cortical and thalamic neurons).Below are the corresponding membrane potential waveforms of the two neurons indicated by the red arrow on the left of the corresponding colorcoded plots.(C), EEG (top trace), AP rastergrams of the firing in each neuronal population for the slow (<1 Hz) oscillation cycle highlighted in (B).Red dashed vertical line represents the first AP of the up-state in each population.The latency (indicated in red below each rastergram) is measured relatively to the first AP of the cycle in the TC FO neuron that fires first (time zero).Below the rastergrams are the corresponding membrane potential waveforms of that cycle for the indicated neuron.(D), Cross-correlations of EEG and membrane potential (black trace) and EEG and APs (color trace) for the indicated neuronal populations, calculated over a 485 s-long simulation.Shaded regions are 95% confidence intervals.(E), Enlargement of the highlighted sections of the membrane potential color plot of the TC FO neurons.L4 FS, FS neuron in cortical layer 4; L4 IB, IB neuron in cortical layer 4; L4 IN, interneurons in cortical layer 4; L4 PY, pyramidal neurons in cortical layer 4; L5 FS, FS neuron in cortical layer 5; L5 IN, interneurons in cortical layer 5; L5 PY, pyramidal neurons in cortical layer 5; L5 RS, RS neuron in cortical layer 5; currents accumulated faster and the up-states terminated earlier in both structures (not shown).Having shorter thalamic up-states shortened the whole oscillation cycle because the duration of thalamic down-states changed little since they are mostly controlled by intrinsic currents.As a result, the next global oscillation cycle was brought forward by an early onset of I T -mediated burst firing of TC FO neurons which initiated the up-states in the cortical neurons (Figure S12).

Figure
Figure4Ashows simulated slow (<1 Hz) oscillation in the full corticothalamic model with up-states occurring with a regular periodicity across all neuronal populations.Most commonly (see the first four up-states in Figure4B), the up-state-linked firing appeared

F I G U R E 5
Delta waves in the full corticothalamic model.(A), EEG showing a period of rhythmic delta waves.(B), EEG (top trace) and color-coded membrane potential plots of the indicated cortical and thalamic neuronal populations during the section of delta waves highlighted in (A) (note the two separate color-scales for the cortical and thalamic neurons).Below are the corresponding membrane potential waveforms of the two neurons indicated by the red arrow on the left of the corresponding color-coded plots.(C), EEG (top trace), AP rastergrams of the firing of each neuronal population for the delta wave cycle highlighted in (B).Red dashed vertical line represents the first AP of the up-state in each population.The latency (indicated in red below each rastergram) is measured relatively to the first AP of the cycle in the TC FO neuron that the fires first (time zero).Below the rastergrams are the corresponding membrane potential waveforms of that cycle for the indicated neuron.(D), Cross-correlations of EEG and membrane potential (black trace) and EEG and APs (color trace) for the indicated neuronal populations, calculated over a 485 s-long simulation.Shaded regions are 95% confidence intervals.Dashed vertical line indicates zero lag.(E), Distribution of the first AP in a delta wave cycle with respect to the EEG for all APs of the indicated neuronal populations.Shaded regions are 95% confidence intervals.Dashed vertical line indicates zero lag.L4 FS, FS neuron in cortical layer 4; L4 IN, interneurons in cortical layer 4; L4 PY, pyramidal neurons in cortical layer 4; L4 RS, RS neuron in cortical layer 4; L5 FS, FS neuron in cortical layer 5; L5 IB, IB neuron in cortical layer 5; L5 IN, interneurons in cortical layer 5; L5 PY, pyramidal neurons in cortical layer 5; NRT FO , firstorder NRT neurons; NRT HO , higher-order NRT neurons; TC FO , first-order TC neurons; TC HO , higher-order TC neurons.sleep states by mimicking the neuronal input resistance changes induced by the hypothalamic and brain stem neuromodulatory drives.Indeed, our simulation results remarkably well reproduced the neuronal firing patterns and their sequence of occurrence during different sleep waves, the membrane potential bistability underlying the slow (<1 Hz) oscillation in single cortical and thalamic neurons, the onset and duration dynamics of up-and down-states, and the increased regularity and frequency imposed by the thalamus on the slow (<1 Hz) oscillation of the isolated neocortical network.
ability of cortical networks to start slow (<1 Hz) oscillations in L5, the thalamocortical input most often prevails when the reciprocal F I G U R E 6 Sleep spindles in the full corticothalamic model.(A), EEG showing the rhythmic pattern of sleep spindles.(B), EEG (top trace) and color-coded membrane potential plots of the indicated cortical and thalamic neuron populations during the sleep spindle highlighted in (A) (note the two separate color-coded scales for the cortex and the thalamus).Below are the corresponding membrane potential waveforms of the two neurons indicated by the red arrow on the left of the corresponding color-coded plots.(C), EEG (top trace), AP rastergrams of the firing of the first AP in each neuronal population for the onset of sleep spindles highlighted in (B).Red dashed vertical line represents the first AP of the up-state in each population.The latency (indicated below each rastergram) is measured relatively to the first AP of the cycle in an NRT FO neuron (time zero).Below the rastergrams are the corresponding membrane potential waveforms of that cycle for the indicated neuron.(D), Cross-correlations of EEG and membrane potential (black trace) and EEG and APs (color trace) for the indicated neuronal populations, calculated over a 485 s-long simulation.Shaded regions are 95% confidence intervals.Dashed vertical line indicates the zero lag.(E), AP distribution with respect to the EEG for all APs of the indicated neuronal populations.Shaded regions are 95% confidence intervals.Dashed vertical line indicates the zero lag.L4 FS, FS neuron in cortical layer 4; L4 IB, IB neuron in cortical layer 4; L4 IN, interneurons in cortical layer 4; L4 PY, pyramidal neurons in cortical layer 4; L5 FS, FS neuron in cortical layer 5; L5 IN, interneurons in cortical layer 5; L5 PY, pyramidal neurons in cortical layer 5; L5 RS, RS neuron in cortical layer 5; NRT FO , first-order NRT neurons; NRT HO , higher-order NRT neurons; TC FO , first-order TC neurons; TC HO , higher-order TC neurons.