Current Source Density Profiles Research Articles

Spatiotemporal dynamics across visual cortical laminae support a predictive coding framework for interpreting mismatch responses.

Context modulates neocortical processing of sensory data. Unexpected visual stimuli elicit large responses in primary visual cortex (V1)-a phenomenon known as deviance detection (DD) at the neural level, or "mismatch negativity" (MMN) when measured with EEG. It remains unclear how visual DD/MMN signals emerge across cortical layers, in temporal relation to the onset of deviant stimuli, and with respect to brain oscillations. Here we employed a visual "oddball" sequence-a classic paradigm for studying aberrant DD/MMN in neuropsychiatric populations-and recorded local field potentials in V1 of awake mice with 16-channel multielectrode arrays. Multiunit activity and current source density profiles showed that although basic adaptation to redundant stimuli was present early (50ms) in layer 4 responses, DD emerged later (150-230ms) in supragranular layers (L2/3). This DD signal coincided with increased delta/theta (2-7Hz) and high-gamma (70-80Hz) oscillations in L2/3 and decreased beta oscillations (26-36Hz) in L1. These results clarify the neocortical dynamics elicited during an oddball paradigm at a microcircuit level. They are consistent with a predictive coding framework, which posits that predictive suppression is present in cortical feed-back circuits, which synapse in L1, whereas "prediction errors" engage cortical feed-forward processing streams, which emanate from L2/3.

Urethane provides an unparalleled anaesthetic model for natural sleep: Commentary on Mondino et al., 2022.

Urethane provides an unparalleled anaesthetic model for natural sleep: Commentary on Mondino et al., 2022.

Emergence of Coordinated Activity in the Developing Entorhinal–Hippocampal Network

Correlated activity in the entorhinal–hippocampal neuronal networks, supported by oscillatory and intermittent population activity patterns is critical for learning and memory. However, when and how correlated activity emerges in these networks during development remains largely unknown. Here, we found that during the first postnatal week in non-anaesthetized head-restrained rats, activity in the superficial layers of the medial entorhinal cortex (MEC) and hippocampus was highly correlated, with intermittent population bursts in the MEC followed by early sharp waves (eSPWs) in the hippocampus. Neurons in the superficial MEC layers fired before neurons in the dentate gyrus, CA3 and CA1. eSPW current-source density profiles indicated that perforant/temporoammonic entorhinal inputs and intrinsic hippocampal connections are co-activated during entorhinal–hippocampal activity bursts. Finally, a majority of the entorhinal–hippocampal bursts were triggered by spontaneous myoclonic body movements, characteristic of the neonatal period. Thus, during the neonatal period, activity in the entorhinal cortex (EC) and hippocampus is highly synchronous, with the EC leading hippocampal activation. We propose that such correlated activity is embedded into a large-scale bottom-up circuit that processes somatosensory feedback resulting from neonatal movements, and that it is likely to instruct the development of connections between neocortex and hippocampus.

Laminar differences in response to simple and spectro-temporally complex sounds in the primary auditory cortex of ketamine-anesthetized gerbils.

In mammals, acoustic communication plays an important role during social behaviors. Despite their ethological relevance, the mechanisms by which the auditory cortex represents different communication call properties remain elusive. Recent studies have pointed out that communication-sound encoding could be based on discharge patterns of neuronal populations. Following this idea, we investigated whether the activity of local neuronal networks, such as those occurring within individual cortical columns, is sufficient for distinguishing between sounds that differed in their spectro-temporal properties. To accomplish this aim, we analyzed simple pure-tone and complex communication call elicited multi-unit activity (MUA) as well as local field potentials (LFP), and current source density (CSD) waveforms at the single-layer and columnar level from the primary auditory cortex of anesthetized Mongolian gerbils. Multi-dimensional scaling analysis was used to evaluate the degree of “call-specificity” in the evoked activity. The results showed that whole laminar profiles segregated 1.8-2.6 times better across calls than single-layer activity. Also, laminar LFP and CSD profiles segregated better than MUA profiles. Significant differences between CSD profiles evoked by different sounds were more pronounced at mid and late latencies in the granular and infragranular layers and these differences were based on the absence and/or presence of current sinks and on sink timing. The stimulus-specific activity patterns observed within cortical columns suggests that the joint activity of local cortical populations (as local as single columns) could indeed be important for encoding sounds that differ in their acoustic attributes.

An analysis of current source density profiles activated by local stimulation in the mouse auditory cortex in vitro

To examine local network properties of the mouse auditory cortex in vitro, we recorded extracellular spatiotemporal laminar profiles driven by short electric local stimulation on a planar multielectrode array substrate. The recorded local field potentials were subsequently evaluated using current source density (CSD) analysis to identify sources and sinks. Current sinks are thought to be an indicator of net synaptic current in the small volume of cortex surrounding the recording site. Thus, CSD analysis combined with multielectrode arrays enabled us to compare mean synaptic activity in response to small current stimuli on a layer-by-layer basis. We also used senescence-accelerated mice (SAM), some strains of which show earlier onset of age-related hearing loss, to examine the characteristic spatiotemporal CSD profiles stimulated by electrodes in specific cortical layers. Thus, the CSD patterns were classified into several clusters based on stimulation sites in the cortical layers. We also found some differences in CSD patterns between the two SAM strains in terms of aging according to principle component analysis with dimension reduction. For simultaneous two-site stimulation, we modeled the obtained CSD profiles as a linear superposition of the CSD profiles to individual single-site stimulation. The model analysis indicated the nonlinearity of spatiotemporal integration over stimulus-driven activity in a layer-specific manner. Finally, on the basis of these results, we discuss the auditory cortex local network properties and the effects of aging on these mouse strains.

Enhanced Excitability of the Neonatal Rat Hippocampus After Acute Exposure to Ethanol

Enhanced excitability manifested by seizures and epilepsy is one of the characteristic features of the fetal alcohol spectrum disorders (FASD). Here, we examined network excitability using a high-potassium model in the hippocampal slices prepared from the postnatal days of P7–9 rats treated with ethanol. Ethanol was administered at 6 g/kg intraperitoneally 12 h before the slice preparation. The extracellular field potential recordings from the hippocampal slices using multishank silicon probes placed along CA3-CA1 axis were performed in the interface chamber. We found that elevation of the extracellular potassium from 3.5 to 6 mM evoked seizure-like clonic or tonic-clonic discharges in 77 % of the slices from the ethanol-treated animals and only in 15 % of the slices from the control animals. Further elevation of the extracellular potassium to 8.5 mM evoked epileptiform activity in 92 and 69 % of the slices from the ethanol-treated and the control animals, respectively. The current source density profile and the multiple unit activity analysis pointed on the CA3 hippocampal region as a generator of the epileptiform activity. Thus, the hippocampal slices from the ethanol-treated neonatal rats display enhanced excitability and could serve as a FASD model to study the early epileptiform transformations following exposure to ethanol.

Texture Segregation Causes Early Figure Enhancement and Later Ground Suppression in Areas V1 and V4 of Visual Cortex.

Segregation of images into figures and background is fundamental for visual perception. Cortical neurons respond more strongly to figural image elements than to background elements, but the mechanisms of figure–ground modulation (FGM) are only partially understood. It is unclear whether FGM in early and mid-level visual cortex is caused by an enhanced response to the figure, a suppressed response to the background, or both.We studied neuronal activity in areas V1 and V4 in monkeys performing a texture segregation task. We compared texture-defined figures with homogeneous textures and found an early enhancement of the figure representation, and a later suppression of the background. Across neurons, the strength of figure enhancement was independent of the strength of background suppression.We also examined activity in the different V1 layers. Both figure enhancement and ground suppression were strongest in superficial and deep layers and weaker in layer 4. The current–source density profiles suggested that figure enhancement was caused by stronger synaptic inputs in feedback-recipient layers 1, 2, and 5 and ground suppression by weaker inputs in these layers, suggesting an important role for feedback connections from higher level areas. These results provide new insights into the mechanisms for figure–ground organization.

Spectral breadth and laminar distribution of thalamocortical inputs to A1.

The GABAergic agonist muscimol is used to inactivate brain regions in order to reveal afferent inputs in isolation. However, muscimol's use in primary auditory cortex (A1) has been questioned on the grounds that it may unintentionally suppress thalamocortical inputs. We tested whether muscimol can preferentially suppress cortical, but not thalamocortical, circuits in urethane-anesthetized mice. We recorded tone-evoked current source density profiles to determine frequency receptive fields (RFs) for three current sinks: the "layer 4" sink (fastest onset, middle-layer sink) and current sinks 100 μm above ("layer 2/3") and 300 μm below ("layer 5/6") the main input. We first determined effects of muscimol dose (0.01-1 mM) on the characteristic frequency (CF) tone-evoked layer 4 sink. An "ideal" dose (100 μM) had no effect on CF-evoked sink onset latency or initial response but reduced peak amplitude by >80%, implying inhibition of intracortical, but not thalamocortical, activity. We extended the analysis to current sinks in layers 2/3 and 5/6 and for all three sinks determined RF breadth (quarter-octave steps, 20 dB above CF threshold). Muscimol reduced RF breadth 42% in layer 2/3 (from 2.4 ± 0.14 to 1.4 ± 0.11 octaves), 14% in layer 4 (2.2 ± 0.12 to 1.9 ± 0.10 octaves), and not at all in layer 5/6 (1.8 ± 0.10 to 1.7 ± 0.12 octaves). The results provide an estimate of the laminar and spectral extent of thalamocortical projections and support the hypothesis that intracortical pathways contribute to spectral integration in A1.

Quantification of mid and late evoked sinks in laminar current source density profiles of columns in the primary auditory cortex

Current source density (CSD) analysis assesses spatiotemporal synaptic activations at somatic and/or dendritic levels in the form of depolarizing current sinks. Whereas many studies have focused on the short (<50 ms) latency sinks, associated with thalamocortical projections, sinks with longer latencies have received less attention. Here, we analyzed laminar CSD patterns for the first 600 ms after stimulus onset in the primary auditory cortex of Mongolian gerbils. By applying an algorithm for contour calculation, three distinct mid and four late evoked sinks were identified in layers I, III, Va, VIa, and VIb. Our results further showed that the patterns of intracortical information-flow remained qualitatively similar for low and for high sound pressure level stimuli at the characteristic frequency (CF) as well as for stimuli ± 1 octave from CF. There were, however, differences associated with the strength, vertical extent, onset latency, and duration of the sinks for the four stimulation paradigms used. Stimuli one octave above the most sensitive frequency evoked a new, and quite reliable, sink in layer Va whereas low level stimulation led to the disappearance of the layer VIb sink. These data indicate the presence of input sources specifically activated in response to level and/or frequency parameters. Furthermore, spectral integration above vs. below the CF of neurons is asymmetric as illustrated by CSD profiles. These results are important because synaptic feedback associated with mid and late sinks—beginning at 50 ms post stimulus latency—is likely crucial for response modulation resulting from higher order processes like memory, learning or cognitive control.

Estimation of Population-Specific Synaptic Currents from Laminar Multielectrode Local Field Potential Recordings

Multielectrode array recordings of extracellular electrical field potentials along the depth axis of the cerebral cortex is an up-and-coming approach for investigating activity of cortical neuronal circuits (Einevoll et al., J. Neurophysiol., 2007; Blomquist et al., PLoS Comp. Biol., 2009). The low-frequency band of extracellular potential, i.e., the local field potential (LFP), is assumed to reflect the synaptic activity and can be used to extract the current source density (CSD) profile. However, physiological interpretation of CSD profile is uncertain because it does not disambiguate synaptic inputs from passive return currents or identifies population-specific contributions to the signal thus obfuscating its interpretation in terms of the excitation flow in the columnar microcircuit.Here we present a novel anatomically informed model for decomposing the LFP signal into population-specific contributions and for estimating the corresponding laminar profiles of synaptic inputs. This involves a combination of 1) the linear forward model, which predicts population-specific laminar LFP in response to synaptic inputs applied to the population of cells having realistic morphologies; and 2) the linear inverse model, which reconstructs laminar profiles of synaptic inputs from laminar LFP data based on the forward prediction. Assuming spatial correlation of synaptic inputs within individual populations, the model decomposes the columnar LFP into population-specific contributions and estimates the corresponding synaptic input laminar profiles less the mean value. Constraining the solution with a priori knowledge of the spatial distribution of synaptic connectivity further allows estimating the strength of active synaptic projections from the columnar LFP profile thus fully specifying synaptic inputs.The capability of the model is demonstrated by applying it to the experimental extracellular data.

Non-homogeneous extracellular resistivity affects the current-source density profiles of up–down state oscillations

Rhythmic local field potential (LFP) oscillations observed during deep sleep are the result of synchronized electrical activities of large neuronal ensembles, which consist of alternating periods of activity and silence, termed 'up' and 'down' states, respectively. Current-source density (CSD) analysis indicates that the up states of these slow oscillations are associated with current sources in superficial cortical layers and sinks in deep layers, while the down states display the opposite pattern of source-sink distribution. We show here that a network model of up and down states displays this CSD profile only if a frequency-filtering extracellular medium is assumed. When frequency filtering was modelled as inhomogeneous conductivity, this simple model had considerably more power in slow frequencies, resulting in significant differences in LFP and CSD profiles compared with the constant-resistivity model. These results suggest that the frequency-filtering properties of extracellular media may have important consequences for the interpretation of the results of CSD analysis.

Heightened Nicotinic Regulation of Auditory Cortex during Adolescence

Adolescent smoking is associated with auditory-cognitive deficits and structural alterations to auditory thalamocortical systems, suggesting that higher auditory function is vulnerable to nicotine exposure during adolescence. Although nicotinic acetylcholine receptors (nAChRs) regulate thalamocortical processing in adults, it is not known whether they regulate processing at earlier ages since their expression pattern changes throughout postnatal development. Here we investigate nicotinic regulation of tone-evoked current source density (CSD) profiles in mouse primary auditory cortex from just after hearing onset until adulthood. At the youngest ages, systemic nicotine did not affect CSD profiles. However, beginning in early adolescence nicotine enhanced characteristic frequency (CF)-evoked responses in layers 2-4 by enhancing thalamocortical, early intracortical, and late intracortical response components. Nicotinic responsiveness developed rapidly and peaked over the course of adolescence, then declined thereafter. Generally, responsiveness in females developed more quickly, peaked earlier, and declined more abruptly and fully than in males. In contrast to the enhancement of CF-evoked responses, nicotine suppressed shorter-latency intracortical responses to spectrally distant (non-CF) stimuli while enhancing longer-latency responses. Intracortical infusion of nAChR antagonists showed that enhancement of CF-evoked intracortical processing involves α4β2*, but not α7, nAChRs, whereas both receptor subtypes regulate non-CF-evoked late intracortical responses. Notably, antagonist effects in females implied regulation by endogenous acetylcholine. Thus, nicotinic regulation of cortical processing varies with age and sex, with peak effects during adolescence that may contribute to the vulnerability of adolescents to smoking.

Dual Mechanism of Neuronal Ensemble Inhibition in Primary Auditory Cortex

Inhibition plays an essential role in shaping and refining the brain's representation of sensory stimulus attributes. In primary auditory cortex (A1), so-called "sideband" inhibition helps to sharpen the tuning of local neuronal responses. Several distinct types of anatomical circuitry could underlie sideband inhibition, including direct thalamocortical (TC) afferents, as well as indirect intracortical mechanisms. The goal of the present study was to characterize sideband inhibition in A1 and to determine its mechanism by analyzing laminar profiles of neuronal ensemble activity. Our results indicate that both lemniscal and nonlemniscal TC afferents play a role in inhibitory responses via feedforward inhibition and oscillatory phase reset, respectively. We propose that the dynamic modulation of excitability in A1 due to the phase reset of ongoing oscillations may alter the tuning of local neuronal ensembles and can be regarded as a flexible overlay on the more obligatory system of lemniscal feedforward type responses.

Current source density correlates of cerebellar Golgi and Purkinje cell responses to tactile input

The overall circuitry of the cerebellar cortex has been known for over a century, but the function of many synaptic connections remains poorly characterized in vivo. We used a one-dimensional multielectrode probe to estimate the current source density (CSD) of Crus IIa in response to perioral tactile stimuli in anesthetized rats and to correlate current sinks and sources to changes in the spike rate of corecorded Golgi and Purkinje cells. The punctate stimuli evoked two distinct early waves of excitation (at <10 and ∼ 20 ms) associated with current sinks in the granular layer. The second wave was putatively of corticopontine origin, and its associated sink was located higher in the granular layer than the first trigeminal sink. The distinctive patterns of granular-layer sinks correlated with the spike responses of corecorded Golgi cells. In general, Golgi cell spike responses could be linearly reconstructed from the CSD profile. A dip in simple-spike activity of coregistered Purkinje cells correlated with a current source deep in the molecular layer, probably generated by basket cell synapses, interspersed between sparse early sinks presumably generated by synapses from granule cells. The late (>30 ms) enhancement of simple-spike activity in Purkinje cells was characterized by the absence of simultaneous sinks in the granular layer and by the suppression of corecorded Golgi cell activity, pointing at inhibition of Golgi cells by Purkinje axon collaterals as a likely mechanism of late Purkinje cell excitation.

On the Estimation of Population-Specific Synaptic Currents from Laminar Multielectrode Recordings

Multielectrode array recordings of extracellular electrical field potentials along the depth axis of the cerebral cortex are gaining popularity as an approach for investigating the activity of cortical neuronal circuits. The low-frequency band of extracellular potential, i.e., the local field potential (LFP), is assumed to reflect synaptic activity and can be used to extract the laminar current source density (CSD) profile. However, physiological interpretation of the CSD profile is uncertain because it does not disambiguate synaptic inputs from passive return currents and does not identify population-specific contributions to the signal. These limitations prevent interpretation of the CSD in terms of synaptic functional connectivity in the columnar microcircuit. Here we present a novel anatomically informed model for decomposing the LFP signal into population-specific contributions and for estimating the corresponding activated synaptic projections. This involves a linear forward model, which predicts the population-specific laminar LFP in response to synaptic inputs applied at different positions along each population and a linear inverse model, which reconstructs laminar profiles of synaptic inputs from laminar LFP data based on the forward model. Assuming spatially smooth synaptic inputs within individual populations, the model decomposes the columnar LFP into population-specific contributions and estimates the corresponding laminar profiles of synaptic input as a function of time. It should be noted that constant synaptic currents at all positions along a neuronal population cannot be reconstructed, as this does not result in a change in extracellular potential. However, constraining the solution using a priori knowledge of the spatial distribution of synaptic connectivity provides the further advantage of estimating the strength of active synaptic projections from the columnar LFP profile thus fully specifying synaptic inputs.

An automated method for detection of layer activation order in information processing pathway of rat barrel cortex under mechanical whisker stimulation

Abstract Rodents perform object localization, texture and shape discrimination very precisely through whisking. During whisking, microcircuits in corresponding barrel columns get activated to segregate and integrate tactile information through the information processing pathway. Sensory signals are projected through the brainstem and thalamus to the corresponding ‘barrel columns’ where different cortical layers are activated during signal projection. Therefore, having precise information about the layer activation order is desirable to better understand this signal processing pathway. This work proposes an automated, computationally efficient and easy to implement method to determine the cortical layer activation from intracortically recorded local field potentials (LFPs) and derived current source density (CSD) profiles: 1. Barrel cortex LFPs are represented by a template of four subsequent events: small positive/negative (E1) → large negative (E2) → slow positive (E3)→ slow long negative (E4). The method exploits the layer specific characteristics of LFPs to obtain latencies of the individual events (E1–E4), then taking the latency of E2 for calculating the layer activation order. 2. The corresponding CSD profile is calculated from the LFPs and the first sink’s peak is considered as a reference point to calculate latencies and evaluate the layer activation order. Other reference points require manual calculation. Similar results of layer activation sequence are found using LFPs and CSDs. Extensive tests on LFPs recorded using standard borosilicate micropipettes demonstrated the method’s workability. An interpretation of layer activation order and CSD profiles on the basis of a simplified interacortical barrel column architecture is also provided.

A Comparison of Sleeplike Slow Oscillations in the Hippocampus Under Ketamine and Urethane Anesthesia

During sleep and anesthesia, a slow (<or=1 Hz) synchronized rhythmic fluctuation of the network activity in the neocortex (nCTX) is prominent. This rhythm, called the slow oscillation (SO), corresponds to sequences of neuronal activity and inactivity (UP and DOWN states) in local and extended networks. Recently, this network pattern has also been described in the hippocampus (HPC) and, interestingly, it shows a dynamic correlation with ongoing SO in the nCTX. However, because of its transient nature during both sleep and urethane anesthesia, studying the SO in the HPC is difficult. Ketamine anesthesia is known to induce a uniform and continuous SO state in the nCTX, but its effects on HPC activity and HPC-nCTX interplay are unknown. Using multisite local field potential recordings, we investigated the effects of ketamine anesthesia on HPC activity and its relation to concurrent nCTX activity. We directly compared the spectral content and spatial distribution of laminar potentials and current source density (CSD) under ketamine to urethane anesthesia. Ketamine evoked stable bouts of SO in the HPC that had a peak frequency of 1.77 +/- 0.08 Hz, slightly higher than that of urethane: 1.17 +/- 0.07 Hz. Laminar voltage and CSD profiles of SO were highly correlated across ketamine and urethane anesthesia, although a higher degree of HPC-nCTX coherence at SO frequencies appeared under ketamine as opposed to urethane. In addition, activity in the gamma bandwidth (30-40 Hz) made up 4.7% of the power spectrum under ketamine, but accounted for only 2.0% of the power spectrum in urethane. Although some differences exist, the SO under ketamine anesthesia appears highly similar to that under urethane.

Frequency-Band Coupling in Surface EEG Reflects Spiking Activity in Monkey Visual Cortex

Although the electroencephalogram (EEG) is widely used in research and clinical settings, its link to the underlying neural activity during sensory processing remains poorly understood. To investigate this, we made simultaneous recordings of surface EEG, intracortical local field potential, and multiunit activity (MUA) in the alert monkey visual cortex during presentation of natural movies. Using a general linear model, we show that in single trials, EEG power in the gamma band (30-100 Hz) and phase in delta band (2-4 Hz) are significant predictors of the MUA response. Specifically, we found that the MUA response was strongest only when increases in EEG gamma power occurred during the negative-going phase of the delta wave, thus revealing a frequency-band coupling mechanism that can be exploited to infer population spiking activity. This finding may open up a new dimension in the use and interpretation of EEG in normal and pathological conditions.

Current Source Density Profiles of Stimulus-Specific Adaptation in Rat Auditory Cortex

Neurons in primary auditory cortex (A1) are known to exhibit a phenomenon known as stimulus-specific adaptation (SSA), which means that, when tested with pure tones, they will respond more strongly to a particular frequency if it is presented as a rare, unexpected "oddball" stimulus than when the same stimulus forms part of a series of common, "standard" stimuli. Although SSA has occasionally been observed in midbrain neurons that form part of the paraleminscal auditory pathway, it is thought to be weak, rare, or nonexistent among neurons of the leminscal pathway that provide the main afferent input to A1, so that SSA seen in A1 is likely generated within A1 by local mechanisms. To study the contributions that neural processing within the different cytoarchitectonic layers of A1 may make to SSA, we recorded local field potentials in A1 of the rat in response to standard and oddball tones and subjected these to current source density analysis. Although our results show that SSA can be observed throughout all layers of A1, right from the earliest part of the response, there are nevertheless significant differences between layers, with SSA becoming significantly stronger as stimulus-related activity passes from the main thalamorecipient layers III and IV to layer V.

Spectrotemporal Analysis of Evoked and Induced Electroencephalographic Responses in Primary Auditory Cortex (A1) of the Awake Monkey

Electroencephalography is increasingly being used to probe the functional organization of auditory cortex. Modulation of the electroencephalographic (EEG) signal by tones was examined in primary auditory cortex (A1) of awake monkeys. EEG data were measured at 4 laminar depths defined by current source density profiles evoked by best frequency (BF) tones. Midlaminar multiunit activity was used to define the tuning characteristics of A1 sites. Presentation of BF tones increased EEG power across the range of frequencies examined (4-290 Hz), with maximal effects evident within the first 100 ms after stimulus onset. The largest relative increases in EEG power generally occurred at very high gamma frequency bands (130-210 Hz). Increases in EEG power for frequencies less than 70 Hz primarily represented changes in phase-locked activity, whereas increases at higher frequencies primarily represented changes in non-phase-locked activity. Power increases in higher gamma bands were better correlated with the A1 tonotopic organization than power increases in lower frequency bands. Results were similar across the 4 laminar depths examined. These findings highlight the value of examining high-frequency EEG components in exploring the functional organization of auditory cortex and may enhance interpretation of related studies in humans.