Spike-wave discharges in absence epilepsy: segregation of electrographic components reveals distinct pathways of seizure activity.

KEY POINTS: The major electrophysiological hallmarks of absence seizures are spike and wave discharges (SWDs), consisting of a sharp spike component and a slow wave component. In a widely accepted scheme, these components are functionally coupled and reflect an iterative progression of neuronal excitation during the spike and post-excitatory silence during the wave. In a genetic rat model of absence epilepsy, local pharmacological inhibition of the centromedian thalamus (CM) selectively suppressed the spike component, leaving self-contained waves in epidural recordings. Thalamic inputs induced activity in cortical microcircuits underlying the spike component, while intracortical oscillations generated the wave component. Based on these findings, we propose a model in which oscillatory waves provide adequate time windows for integration of thalamocortical inputs and feedback responses during generation of a synchronized SWD. ABSTRACT: Spike and wave discharges (SWDs) are the electrographic hallmark of absence seizures and the major diagnostic criterion for childhood absence epilepsy (CAE). In a widely accepted scheme, the alternating sequence of spikes and waves reflects an iterative progression of neuronal excitation during the spike component and post-excitatory silence during the wave component. Here we challenge this view by showing that these two components are not necessarily coupled. In a genetic rat model of CAE, self-contained waves occurred in motor cortex in synchrony with SWDs in the somatosensory system during blockade of afferent input from the thalamus. Current-source density analyses of multi-site local field potentials (LFPs) revealed layer-specific activity, in which thalamic inputs induced a sequence of cellular-synaptic events underlying the spike component, while intracortical oscillations generated the wave component. These findings indicate novel principles of SWDs, where oscillatory cortical waves provide adequate time windows for integration of thalamocortical inputs and feedback responses during generation of seizure activity.

Laminar analysis of the slow wave activity in the somatosensory cortex of anesthetized rats.

Rhythmic slow waves characterize brain electrical activity during natural deep sleep and under anesthesia, reflecting the synchronous membrane potential fluctuations of neurons in the thalamocortical network. Strong evidence indicates that the neocortex plays an important role in the generation of slow wave activity (SWA), however, contributions of individual cortical layers to the SWA generation are still unclear. The anatomically correct laminar profiles of SWA were revealed under ketamine/xylazine anesthesia, with combined local field potential recordings, multiple-unit activity (MUA), current source density (CSD) and time-frequency analyses precisely co-registered with histology. The up-state related negative field potential wave showed the largest amplitude in layer IV, the CSD was largest in layers I and III, whereas MUA was maximal in layer V, suggesting spatially dissociated firing and synaptic/transmembrane processes in the rat somatosensory cortex. Up-state related firing could start in virtually any layers (III-VI) of the cortex, but were most frequently initiated in layer V. However, in a subset of experiments, layer IV was considerably active in initiating up-state related MUA even in the absence of somatosensory stimulation. Somatosensory stimulation further strengthened up-state initiation in layer IV. Our results confirm that cortical layer V firing may have a major contribution to the up-state generation of ketamine/xylazine-induced SWA, however, thalamic influence through the thalamorecipient layer IV can also play an initiating role, even in the absence of sensory stimulation.

Invariance in current dipole moment density across brain structures and species: physiological constraint for neuroimaging.

Although anatomical constraints have been shown to be effective for MEG and EEG inverse solutions, there are still no effective physiological constraints. Strength of the current generator is normally described by the moment of an equivalent current dipole Q. This value is quite variable since it depends on size of active tissue. In contrast, the current dipole moment density q, defined as Q per surface area of active cortex, is independent of size of active tissue. Here we studied whether the value of q has a maximum in physiological conditions across brain structures and species. We determined the value due to the primary neuronal current (q primary) alone, correcting for distortions due to measurement conditions and secondary current sources at boundaries separating regions of differing electrical conductivities. The values were in the same range for turtle cerebellum (0.56-1.48 nAm/mm(2)), guinea pig hippocampus (0.30-1.34 nAm/mm(2)), and swine neocortex (0.18-1.63 nAm/mm(2)), rat neocortex (~2.2 nAm/mm(2)), monkey neocortex (~0.40 nAm/mm(2)) and human neocortex (0.16-0.77 nAm/mm(2)). Thus, there appears to be a maximum value across the brain structures and species (1-2 nAm/mm(2)). The empirical values closely matched the theoretical values obtained with our independently validated neural network model (1.6-2.8 nAm/mm(2) for initial spike and 0.7-3.1 nAm/mm(2) for burst), indicating that the apparent invariance is not coincidental. Our model study shows that a single maximum value may exist across a wide range of brain structures and species, varying in neuron density, due to fundamental electrical properties of neurons. The maximum value of q primary may serve as an effective physiological constraint for MEG/EEG inverse solutions.

Midazolam and atropine alter theta oscillations in the hippocampal CA1 region by modulating both the somatic and distal dendritic dipoles.

Theta (4-12 Hz) oscillations in the hippocampus play an important role in learning and memory. They are altered by a wide variety of drugs that impair memory, and these effects may underlie or contribute to drug-induced amnesia. However, the network mechanisms linking drug actions with changes in memory formation remain poorly defined. Here, we used a multisite linear electrode array to measure local field potentials simultaneously across the CA1 layers of the hippocampus during active exploration, and employed current source density analysis and computational modeling to investigate how midazolam and atropine-two amnestic drugs that are used clinically and experimentally-change the relative timing and strength of the drivers of θ-oscillations. We found that two dipoles are present, with active inputs that are centered at the soma and the distal apical dendrite and passive return pathways that overlap in the mid-apical dendrite. Both drugs shifted the position of the phase reversal in the local field potential that occurred in the mid-apical dendritic region, but in opposite directions, by changing the strength of the dendritic pole, without altering the somatic pole or relative timing. Computational modeling showed that this constellation of changes, as well as an additional effect on a variably present mid-apical pole, could be produced by simultaneous changes in the active somatic and distal dendritic inputs. These network-level changes, produced by two amnestic drugs that target different types of receptors, may thus serve as a common basis for impaired memory encoding.

In vitro intrinsic optical imaging can be used for source determination in cortical slices.

In the last decades intrinsic optical imaging has become a widely used technique for monitoring activity in vivo and in vitro. It is assumed that in brain slices the source of intrinsic optical signals (IOSs) is the change in light scattering caused by cell swelling or shrinkage. The aim of the present study was to find a correlation between electrical activity and parallel optical characteristics, elicited by 4-aminopyridine-containing or Mg(2+) -free medium in rat cortical brain slices. Electrophysiological signals and reflected light alterations were recorded during spontaneous seizure activity. Current source density (CSD) analysis was performed on the electrophysiological records. Direct correlation analysis of IOS to CSD was made, and source distribution provided by IOS and CSD methods was compared by determining Matthews correlation coefficient. The gradual development of seizure-like activity elicited the reduction of light reflectance. The main findings of our experiments are that long-term epileptiform activity resulted in persistent alteration in IOSs of brain slices. The observed IOS pattern remained stable after 1 h incubation in convulsants. The pattern of IOS shows good correlation with the data obtained from the CSD analysis. Persistent IOS changes provide information about the area-specific changes of basic excitability, which can serve as a background for ictal and interictal-like epileptiform activity. We can conclude that changes in IOSs correlate well with electrophysiological recordings under different conditions. Our experiments provide evidence that underlying synchronised neuronal processes produce parallel alterations in IOSs and electrophysiological activity.

Current Source Density Analysis of Electroantennogram Recordings: A Tool for Mapping the Olfactory Response in an Insect Antenna.

The set of chemosensory receptors expressed by the olfactory receptor neurons lying in an insect's antennae and maxillary palps define the ability of this insect to perceive the volatile chemicals of its environment. The main two electrophysiological methods of antennal recordings for studying the range of chemicals that activate chemosensory receptors have limitations. Single-sensillum recording (SSR) samples a subset of olfactory receptor neurons and therefore does not reveal the full capacity of an insect to perceive an odor. Electroantennography (EAG), even if less resolutive than SSRs, is sometimes preferred since it samples the activity of a large number of the olfactory receptor neurons. But, at least in flies, the amplitude of the EAG signal is not directly correlated with the degree of sensitivity of the insect to the olfactory compound. Such dual methodology was also used to study mammalian brains, and the current source density (CSD) analysis was developed to bridge the gap between the cellular and the population recordings. This paper details the use of a similar approach adapted to the study of olfactory responses within insects with bulbous antennae. The EAG was recorded at multiple antennal positions and the CSD that generates the EAG potentials were estimated. The method measures the activation of olfactory receptor neurons (ORNs) across the antennae and thus it quantifies the olfactory sensitivity of the insect. It allows a rapid mapping of olfactory responses and thus can be used to guide further SSRs or to determine that two chemicals are detected by independent ORNs. This study further explored biases resulting from a limited number of recording positions or from an approximation of the antennal geometry that should be considered for interpreting the CSD maps. It also shows that the CSD analysis of EAGs is compatible with a gas chromatograph stimulator for analyzing the response to complex odors. Finally, I discuss the origin of the EAG signal in light of the CSD theory.

waveCSD: A method for estimating transmembrane currents originated from propagating neuronal activity in the neocortex: Application to study cortical spreading depression.

BACKGROUND: Recent years have witnessed an upsurge in the development of methods for estimating current source densities (CSDs) in the neocortical tissue from their recorded local field potential (LFP) reflections using microelectrode arrays. Among these, methods utilizing linear arrays work under the assumption that CSDs vary as a function of cortical depth; whereas they are constant in the tangential direction, infinitely or in a confined cylinder. This assumption is violated in the analysis of neuronal activity propagating along the neocortical sheet, e.g. propagation of alpha waves or cortical spreading depression. NEW METHOD: Here, we developed a novel mathematical method (waveCSD) for CSD analysis of LFPs associated with a planar wave of neocortical neuronal activity propagating at a constant velocity towards a linear probe. RESULTS: Results show that the algorithm is robust to the presence of noise in LFP data and uncertainties in knowledge of propagation velocity. Also, results show high level of accuracy of the method in a wide range of electrode resolutions. Using in vivo experimental recordings from the rat neocortex, we employed waveCSD to characterize transmembrane currents associated with cortical spreading depressions. COMPARISON WITH EXISTING METHOD(S): Simulation results indicate that waveCSD has a significantly higher reconstruction accuracy compared to the widely-used inverse CSD method (iCSD), and the regularized kernel CSD method (kCSD), in the analysis of CSDs originating from propagating neuronal activity. CONCLUSIONS: The waveCSD method provides a theoretical platform for estimation of transmembrane currents from their LFPs in experimental paradigms involving wave propagation.

Interhemispheric Transfer Time Asymmetry of Visual Information Depends on Eye Dominance: An Electrophysiological Study.

The interhemispheric transfer of information is a fundamental process in the human brain. When a visual stimulus appears eccentrically in one visual-hemifield, it will first activate the contralateral hemisphere but also the ipsilateral one with a slight delay due to the interhemispheric transfer. This interhemispheric transfer of visual information is believed to be faster from the right to the left hemisphere in right-handers. Such an asymmetry is considered as a relevant fact in the context of the lateralization of the human brain. We show here using current source density (CSD) analyses of visually evoked potential (VEP) that, in right-handers and, to a lesser extent in left-handers, this asymmetry is in fact dependent on the sighting eye dominance, the tendency we have to prefer one eye for monocular tasks. Indeed, in right-handers, a faster interhemispheric transfer of visual information from the right to left hemisphere was observed only in participants with a right dominant eye (DE). Right-handers with a left DE showed the opposite pattern, with a faster transfer from the left to the right hemisphere. In left-handers, albeit a smaller number of participants has been tested and hence confirmation is required, only those with a right DE showed an asymmetrical interhemispheric transfer with a faster transfer from the right to the left hemisphere. As a whole these results demonstrate that eye dominance is a fundamental determinant of asymmetries in interhemispheric transfer of visual information and suggest that it is an important factor of brain lateralization.

Spatial Localization of Sources in the Rat Subthalamic Motor Region Using an Inverse Current Source Density Method.

Objective: In this study we introduce the use of the current source density (CSD) method as a way to visualize the spatial organization of evoked responses in the rat subthalamic nucleus (STN) at fixed time stamps resulting from motor cortex stimulation. This method offers opportunities to visualize neuronal input and study the relation between the synaptic input and the neural output of neural populations. Approach: Motor cortex evoked local field potentials and unit activity were measured in the subthalamic region, with a 3D measurement grid consisting of 320 measurement points and high spatial resolution. This allowed us to visualize the evoked synaptic input by estimating the current source density (CSD) from the measured local field potentials, using the inverse CSD method. At the same time, the neuronal output of the cells within the grid is assessed by calculating post stimulus time histograms. Main results: The CSD method resulted in clear and distinguishable sources and sinks of the neuronal input activity in the STN after motor cortex stimulation. We showed that the center of the synaptic input of the STN from the motor cortex is located dorsal to the input from globus pallidus. Significance: For the first time we have performed CSD analysis on motor cortex stimulation evoked LFP responses in the rat STN as a proof of principle. Our results suggest that the CSD method can be used to gain new insights into the spatial extent of synaptic pathways in brain structures.

Spatiotemporal patterns of current source density in the prefrontal cortex of a behaving monkey.

One of the fundamental missions of neuroscience is to explore the input and output properties of neuronal networks to reveal their functional significance. However, it is technically difficult to examine synaptic inputs into neuronal circuits in behaving animals. Here, we conducted current source density (CSD) analysis on local field potentials (LFPs) recorded simultaneously using a multi-contact electrode in the prefrontal cortex (PFC) of a behaving monkey. We observed current sink task-dependent spatiotemporal patterns considered to reflect the synaptic input to neurons adjacent to the recording site. Specifically, the inferior convex current sink in the PFC was dominant during the delay period, whereas the current sink was prominent in the principal sulcus during the sample cue and test cue periods. Surprisingly, sulcus current sink patterns were spatially periodic, which corresponds to the columnar structure suggested by previous anatomical studies. The approaches used in the current study will help to elucidate how the PFC network performs executive functions according to its synaptic input.

Early exposure to urethane anesthesia: Effects on neuronal activity in the piriform cortex of the developing brain.

Exposure to urethane anesthesia reportedly produces selective neuronal cell loss in the piriform cortex of young brains; however, resulting functional deficits have not been investigated. The present study found abnormalities in piriform cortex activity of isolated brains in vitro that were harvested from guinea pigs exposed to urethane anesthesia at 14 days of age. Current source density (CSD) analysis and voltage-sensitive dye (VSD) imaging experiments were conducted 48h after urethane injection. We applied paired-pulse stimulation to the lateral olfactory tract (LOT) and assessed short-interval intra-cortical inhibition in the piriform cortex. CSD analysis revealed that a current sink in layer Ib remained active in response to successive stimuli, with an inter-stimulus interval of 30-60 ms, which was typically strongly inhibited. VSD imaging demonstrated stronger and extended neural activity in the urethane-treated piriform cortex, even in response to a second stimulus delivered in short succession. We identified gamma-aminobutyric acid (GABA) ergic neurons in the piriform cortex of sham and urethane-treated animals and found a decrease in GABA-immunoreactive cell density in the urethane group. These results suggest that urethane exposure induces loss of GABAergic interneurons and a subsequent reduction in paired-pulse inhibition in the immature piriform cortex.

Spatiotemporal characteristics and pharmacological modulation of multiple gamma oscillations in the CA1 region of the hippocampus.

Multiple components of "γ-oscillations" between 30-170 Hz in the CA1 region of the hippocampus have been described, based on their coherence with oscillations in other brain regions and on their cross-frequency coupling with local θ-oscillations. However, it remains unclear whether the different sub-bands are generated by a single broadband oscillator coupled to multiple external inputs, or by separate oscillators that incorporate distinct circuit elements. To distinguish between these possibilities, we used high-density linear array recording electrodes in awake behaving mice to examine the spatiotemporal characteristics of γ-oscillations and their responses to midazolam and atropine. We characterized oscillations using current source density (CSD) analysis, and measured θ-γ phase-amplitude coupling by cross frequency coupling (CFC) analysis. Prominent peaks were present in the CSD signal in the mid- and distal apical dendritic layers at all frequencies, and at stratum pyramidale for γ(slow) (30-45 Hz) and γ(mid) (50-90 Hz), but not γ(fast) (90-170 Hz) oscillations. Differences in the strength and timing of θ-γ(slow) and θ-γ(mid) cross frequency coupling, and a lack of coupling at the soma and mid-apical region for γ(fast) oscillations, indicated that separate circuit components generate the three sub-bands. Midazolam altered CSD amplitudes and cross-frequency coupling in a lamina- and frequency specific manner, providing further evidence for separate generator circuits. Atropine altered CSD amplitudes and θ-γ CFC uniformly at all locations. Simulations using a detailed compartmental model were consistent with γ(slow) and γ(mid) oscillations driven primarily by inputs at the mid-apical dendrites, and γ(fast) at the distal apical dendrite. Our results indicate that multiple distinct local circuits generate γ-oscillations in the CA1 region of the hippocampus, and provide detailed information about their spatiotemporal characteristics.