Distinct superficial and deep laminar domains of activity in the visual cortex during rest and stimulation

Understanding music and aging through the lens of Bayesian inference

Bayesian inference has recently gained momentum in explaining music perception and aging. A fundamental mechanism underlying Bayesian inference is the notion of prediction. This framework could explain how predictions pertaining to musical (melodic, rhythmic, harmonic) structures engender action, emotion, and learning, expanding related concepts of music research, such as musical expectancies, groove, pleasure, and tension. Moreover, a Bayesian perspective of music perception may shed new insights on the beneficial effects of music in aging. Aging could be framed as an optimization process of Bayesian inference. As predictive inferences refine over time, the reliance on consolidated priors increases, while the updating of prior models through Bayesian inference attenuates. This may affect the ability of older adults to estimate uncertainties in their environment, limiting their cognitive and behavioral repertoire. With Bayesian inference as an overarching framework, this review synthesizes the literature on predictive inferences in music and aging, and details how music could be a promising tool in preventive and rehabilitative interventions for older adults through the lens of Bayesian inference.

The laminar organization of cell types in macaque cortex and its relationship to neuronal oscillations

The canonical microcircuit (CMC) has been hypothesized to be the fundamental unit of information processing in cortex. Each CMC unit is thought to be an interconnected column of neurons with specific connections between excitatory and inhibitory neurons across layers. Recently, we identified a conserved spectrolaminar motif of oscillatory activity across the primate cortex that may be the physiological consequence of the CMC. The spectrolaminar motif consists of local field potential (LFP) gamma-band power (40-150 Hz) peaking in superficial layers 2 and 3 and alpha/beta-band power (8-30 Hz) peaking in deep layers 5 and 6. Here, we investigate whether specific conserved cell types may produce the spectrolaminar motif. We collected laminar histological and electrophysiological data in 11 distinct cortical areas spanning the visual hierarchy: V1, V2, V3, V4, TEO, MT, MST, LIP, 8A/FEF, PMD, and LPFC (area 46), and anatomical data in DP and 7A. We stained representative slices for the three main inhibitory subtypes, Parvalbumin (PV), Calbindin (CB), and Calretinin (CR) positive neurons, as well as pyramidal cells marked with Neurogranin (NRGN). We found a conserved laminar structure of PV, CB, CR, and pyramidal cells. We also found a consistent relationship between the laminar distribution of inhibitory subtypes with power in the local field potential. PV interneuron density positively correlated with gamma (40-150 Hz) power. CR and CB density negatively correlated with alpha (8-12 Hz) and beta (13-30 Hz) oscillations. The conserved, layer-specific pattern of inhibition and excitation across layers is therefore likely the anatomical substrate of the spectrolaminar motif.

Significance Statement

Neuronal oscillations emerge as an interplay between excitatory and inhibitory neurons and underlie cognitive functions and conscious states. These oscillations have distinct expression patterns across cortical layers. Does cellular anatomy enable these oscillations to emerge in specific cortical layers? We present a comprehensive analysis of the laminar distribution of the three main inhibitory cell types in primate cortex (Parvalbumin, Calbindin, and Calretinin positive) and excitatory pyramidal cells. We found a canonical relationship between the laminar anatomy and electrophysiology in 11 distinct primate areas spanning from primary visual to prefrontal cortex. The laminar anatomy explained the expression patterns of neuronal oscillations in different frequencies. Our work provides insight into the cortex-wide cellular mechanisms that generate neuronal oscillations in primates.

A ubiquitous spectrolaminar motif of local field potential power across the primate cortex

The mammalian cerebral cortex is anatomically organized into a six-layer motif. It is currently unknown whether a corresponding laminar motif of neuronal activity patterns exists across the cortex. Here we report such a motif in the power of local field potentials (LFPs). Using laminar probes, we recorded LFPs from 14 cortical areas across the cortical hierarchy in five macaque monkeys. The laminar locations of recordings were histologically identified by electrolytic lesions. Across all areas, we found a ubiquitous spectrolaminar pattern characterized by an increasing deep-to-superficial layer gradient of high-frequency power peaking in layers 2/3 and an increasing superficial-to-deep gradient of alpha-beta power peaking in layers 5/6. Laminar recordings from additional species showed that the spectrolaminar pattern is highly preserved among primates-macaque, marmoset and human-but more dissimilar in mouse. Our results suggest the existence of a canonical layer-based and frequency-based mechanism for cortical computation.

High-Resolution Laminar Identification in Macaque Primary Visual Cortex Using Neuropixels Probes

Laminar electrode arrays allow simultaneous recording of activity of many cortical neurons and assignment to correct layers using current source density (CSD) analyses. Electrode arrays with 100-micron contact spacing can estimate borders between layer 4 versus superficial or deep layers, but in macaque primary visual cortex (V1) there are far more layers, such as 4A which is only 50-100 microns thick. Neuropixels electrode arrays have 20-micron spacing, and thus could potentially discern thinner layers and more precisely identify laminar borders. Here we show that CSD signals lack the spatial resolution required to take advantage of high density Neuropixels arrays and describe the development of approaches based on higher resolution electrical signals and analyses, including spike waveforms and spatial spread, unit density, high-frequency action potential (AP) power spectrum, temporal power change, and coherence spectrum, that afford far higher resolution of laminar distinctions, including the ability to precisely detect the borders of even the thinnest layers of V1.

Neural Oscillations and Multisensory Processing

Neural oscillations play a role in sensory processing by coordinating synchronized neuronal activity. Synchronization of gamma oscillations is engaged in local computation of feedforward signals and synchronization of alpha-beta oscillations is engaged in feedback processing over long-range areas. These spatially and spectrally segregated bi-directional signals may be integrated by a mechanism of cross-frequency coupling. Synchronization of neural oscillations has also been proposed as a mechanism for information integration across multiple sensory modalities. A transient stimulus or rhythmic stimulus from one modality may lead to phase alignment of ongoing neural oscillations in multiple sensory cortices, through a mechanism of cross-modal phase reset or cross-modal neural entrainment. Synchronized activities in multiple sensory cortices are more likely to boost stronger activities in downstream areas. Compared to synchronized oscillations, asynchronized oscillations may impede signal processing, and may contribute to sensory selection by setting the oscillations in the target-related cortex and the oscillations in the distractor-related cortex to opposite phases.

Enabling brain-wide mapping of directed functional connectivity at 3T via layer-dependent fMRI with draining-vein suppression

Layer-dependent functional magnetic resonance imaging (fMRI) offers a compelling avenue for investigating directed functional connectivity (FC). To construct a comprehensive map of brain-wide directed FC, several technical criteria must be met, including sub-mm spatial resolution, adequate temporal resolution, functional sensitivity, global brain coverage, and high spatial specificity. Although gradient echo (GE)-based echo planar imaging (EPI) is commonly used for rapid fMRI acquisition, it faces significant challenges due to the draining-vein effect, particularly when utilizing blood oxygen level-dependent (BOLD) contrast. In this study, we mitigated this effect by incorporating velocity-nulling (VN) gradients into a GE-BOLD fMRI sequence, opting for a 3T magnetic field strength over 7T. We also integrated several advanced techniques, such as simultaneous multi-slice (SMS) acceleration and NORDIC denoising, to enhance temporal resolution, spatial coverage, and signal sensitivity. Collectively, the VN fMRI method exhibited notable spatial specificity, as evidenced by the identification of double-peak activation patterns within the primary motor cortex (M1) during a finger-tapping task. Additionally, the technique demonstrated BOLD sensitivity in the lateral geniculate nucleus (LGN). Furthermore, our VN fMRI technique displayed superior robustness when compared to conventional fMRI approaches across participants. Our findings of directed FC elucidate several layer-specific functional relationships between different brain regions and align closely with existing literature. Given the widespread availability of 3T scanners, this technical advancement has the potential for significant impact across multiple domains of neuroscience research.

Feedforward attentional selection in sensory cortex

Salient objects grab attention because they stand out from their surroundings. Whether this phenomenon is accomplished by bottom-up sensory processing or requires top-down guidance is debated. We tested these alternative hypotheses by measuring how early and in which cortical layer(s) neural spiking distinguished a target from a distractor. We measured synaptic and spiking activity across cortical columns in mid-level area V4 of male macaque monkeys performing visual search for a color singleton. A neural signature of attentional capture was observed in the earliest response in the input layer 4. The magnitude of this response predicted response time and accuracy. Errant behavior followed errant selection. Because this response preceded top-down influences and arose in the cortical layer not targeted by top-down connections, these findings demonstrate that feedforward activation of sensory cortex can underlie attentional priority.

Cortical recurrence supports resilience to sensory variance in the primary visual cortex

Our daily endeavors occur in a complex visual environment, whose intrinsic variability challenges the way we integrate information to make decisions. By processing myriads of parallel sensory inputs, our brain is theoretically able to compute the variance of its environment, a cue known to guide our behavior. Yet, the neurobiological and computational basis of such variance computations are still poorly understood. Here, we quantify the dynamics of sensory variance modulations of cat primary visual cortex neurons. We report two archetypal neuronal responses, one of which is resilient to changes in variance and co-encodes the sensory feature and its variance, improving the population encoding of orientation. The existence of these variance-specific responses can be accounted for by a model of intracortical recurrent connectivity. We thus propose that local recurrent circuits process uncertainty as a generic computation, advancing our understanding of how the brain handles naturalistic inputs.

Near-field potentials index local neural computations more accurately than population spiking

Local field potentials (LFP) are low-frequency extracellular voltage fluctuations thought to primarily arise from synaptic activity. However, unlike highly localized neuronal spiking, LFP is spatially less specific. LFP measured at one location is not entirely generated there due to far-field contributions that are passively conducted across volumes of neural tissue. We sought to quantify how much information within the locally generated, near-field low-frequency activity (nfLFP) is masked by volume-conducted far-field signals. To do so, we measured laminar neural activity in primary visual cortex (V1) of monkeys viewing sequences of multifeatured stimuli. We compared information content of regular LFP and nfLFP that was mathematically stripped of volume-conducted far-field contributions. Information content was estimated by decoding stimulus properties from neural responses via spatiotemporal multivariate pattern analysis. Volume-conducted information differed from locally generated information in two important ways: (1) for stimulus features relevant to V1 processing (orientation and eye-of-origin), nfLFP contained more information. (2) in contrast, the volume-conducted signal was more informative regarding temporal context (relative stimulus position in a sequence), a signal likely to be coming from elsewhere. Moreover, LFP and nfLFP differed both spectrally as well as spatially, urging caution regarding the interpretations of individual frequency bands and/or laminar patterns of LFP. Most importantly, we found that population spiking of local neurons was less informative than either the LFP or nfLFP, with nfLFP containing most of the relevant information regarding local stimulus processing. These findings suggest that the optimal way to read out local computational processing from neural activity is to decode the local contributions to LFP, with significant information loss hampering both regular LFP and local spiking.

Author’s Contributions

Conceptualization, D.A.T., J.A.W, and A.M.; Data Collection, J.A.W., M.A.C., K.D.; Formal Analysis, D.A.T. and J.A.W.; Data Visualization, D.A.T. and J.A.W.; Original Draft, D.A.T., J.A.W., and A.M.; Revisions and Final Draft, D.A.T., J.A.W., M.A.C., K.D., M.T.W., A.M.B., and A.M.

Competing Interests

The authors declare no conflicts of interest.

Sleep/wake changes in perturbational complexity in rats and mice

In humans, the level of consciousness is assessed by quantifying the spatiotemporal complexity of cortical responses using Perturbational Complexity Index (PCI) and related PCI^st (st, state transitions). Here we validate PCI^st in freely moving rats and mice by showing that it is lower in NREM sleep and slow wave anesthesia than in wake or REM sleep, as in humans. We then show that (1) low PCI^st is associated with the occurrence of an OFF period of neuronal silence; (2) stimulation of deep, but not superficial, cortical layers leads to reliable PCI^st changes across sleep/wake and anesthesia; (3) consistent PCI^st changes are independent of which single area is being stimulated or recorded, except for recordings in mouse prefrontal cortex. These experiments show that PCI^st can reliably measure vigilance states in unresponsive animals and support the hypothesis that it is low when an OFF period disrupts causal interactions in cortical networks.

Decreased coherence in the model of the dorsal visual pathway associated with Alzheimer's disease

Decreased coherence in electroencephalogram (EEG) has been reported in Alzheimer's disease (AD) experimentally, which could be considered as a typical electrophysiological characteristic in AD. This work aimed to investigate the effect of AD on coherence in the dorsal visual pathway by the technique of neurocomputation. Firstly, according to the hierarchical organization of the cerebral cortex and the information flows of the dorsal visual pathway, a more physiologically plausible neural mass model including cortical areas v1, v2, and v5 was established in the dorsal visual pathway. The three interconnected cortical areas were connected by ascending and descending projections. Next, the pathological condition of loss of long synaptic projections in AD was simulated by reducing the parameters of long synaptic projections in the model. Then, the loss of long synaptic projections on coherence among different visual cortex areas was explored by means of power spectral analysis and coherence function. The results demonstrate that the coherence between these interconnected cortical areas showed an obvious decline with the gradual decrease of long synaptic projections, i.e. decrease in descending projections from area v2 to v1 and v5 to v2 and ascending projection from area v2 to v5. Hopefully, the results of this study could provide theoretical guidance for understanding the dynamical mechanism of AD.

LFP polarity changes across cortical and eccentricity in primary visual cortex

Local field potentials (LFPs) can evaluate neural population activity in the cortex and their interaction with other cortical areas. Analyzing current source density (CSD) rather than LFPs is very significant due to the reduction of volume conduction effects. Current sinks are construed as net inward transmembrane currents, while current sources are net outward ones. Despite extensive studies of LFPs and CSDs, their morphology in different cortical layers and eccentricities are still largely unknown. Because LFP polarity changes provide a measure of neural activity, they can be useful in implanting brain-computer interface (BCI) chips and effectively communicating the BCI devices to the brain. We hypothesize that sinks and sources analyses could be a way to quantitatively achieve their characteristics in response to changes in stimulus size and layer-dependent differences with increasing eccentricities. In this study, we show that stimulus properties play a crucial role in determining the flow. The present work focusses on the primary visual cortex (V1). In this study, we investigate a map of the LFP-CSD in V1 area by presenting different stimulus properties (e.g., size and type) in the visual field area of Macaque monkeys. Our aim is to use the morphology of sinks and sources to measure the input and output information in different layers as well as different eccentricities. According to the value of CSDs, the results show that the stimuli smaller than RF's size had lower strength than the others and the larger RF's stimulus size showed smaller strength than the optimized stimulus size, which indicated the suppression phenomenon. Additionally, with the increased eccentricity, CSD's strengths were increased across cortical layers.

Weakly Correlated Local Cortical State Switches under Anesthesia Lead to Strongly Correlated Global States

During recovery from anesthesia, brain activity switches abruptly between a small set of discrete states. Surprisingly, this switching also occurs under constant doses of anesthesia, even in the absence of stimuli. These metastable states and the transitions between them are thought to form a "scaffold" that ultimately guides the brain back to wakefulness. The processes that constrain cortical activity patterns to these states and govern how states are coordinated between different cortical regions are unknown. If state transitions were driven by subcortical modulation, different cortical sites should exhibit near-synchronous state transitions. Conversely, spatiotemporal heterogeneity would suggest that state transitions are coordinated through corticocortical interactions. To differentiate between these hypotheses, we quantified synchrony of brain states in male rats exposed to a fixed isoflurane concentration. States were defined from spectra of local field potentials recorded across layers of visual and motor cortices. A transition synchrony measure shows that most state transitions are highly localized. Furthermore, while most pairs of cortical sites exhibit statistically significant coupling of both states and state transition times, coupling strength is typically weak. States and state transitions in the thalamic input layer (L4) are particularly decoupled from those in supragranular and infragranular layers. This suggests that state transitions are not imposed on the cortex by broadly projecting modulatory systems. Although each pairwise interaction is typically weak, we show that the multitude of such weak interactions is sufficient to confine global activity to a small number of discrete states.SIGNIFICANCE STATEMENT The brain consistently recovers to wakefulness after anesthesia, but this process is poorly understood. Previous work revealed that, during recovery from anesthesia, corticothalamic activity falls into one of several discrete patterns. The neuronal mechanisms constraining the cortex to just a few discrete states remain unknown. Global states could be coordinated by fluctuations in subcortical nuclei that project broadly to the cortex. Alternatively, these states may emerge from interactions within the cortex itself. Here, we provide evidence for the latter possibility by demonstrating that most pairs of cortical sites exhibit weak coupling. We thereby lay groundwork for future investigations of the specific cellular and network mechanisms of corticocortical activity state coupling.

Simultaneous Observation of Mouse Cortical and Hippocampal Neural Dynamics under Anesthesia through a Cranial Microprism Window

The fluorescence microscope has been widely used to explore dynamic processes in vivo in mouse brains, with advantages of a large field-of-view and high spatiotemporal resolution. However, owing to background light and tissue scattering, the single-photon wide-field microscope fails to record dynamic neural activities in the deep brain. To achieve simultaneous imaging of deep-brain regions and the superficial cortex, we combined the extended-field-of-view microscopy previously proposed with a novel prism-based cranial window to provide a longitudinal view. As well as a right-angle microprism for imaging above 1 mm, we also designed a new rectangular-trapezoidal microprism cranial window to extend the depth of observation to 1.5 mm and to reduce brain injury. We validated our method with structural imaging of microglia cells in the superficial cortex and deep-brain regions. We also recorded neuronal activity from the mouse brains in awake and anesthesitized states. The results highlight the great potential of our methods for simultaneous dynamic imaging in the superficial and deep layers of mouse brains.

Anesthetic modulations dissociate neuroelectric characteristics between sensory-evoked and spontaneous activities across bilateral rat somatosensory cortical laminae

Spontaneous neural activity has been widely adopted to construct functional connectivity (FC) amongst distant brain regions. Although informative, the functional role and signaling mechanism of the resting state FC are not intuitive as those in stimulus/task-evoked activity. In order to bridge the gap, we investigated anesthetic modulation of both resting-state and sensory-evoked activities. We used two well-studied GABAergic anesthetics of varying dose (isoflurane: 0.5–2.0% and α-chloralose: 30 and 60 mg/kg∙h) and recorded changes in electrophysiology using a pair of laminar electrode arrays that encompass the entire depth of the bilateral somatosensory cortices (S1fl) in rats. Specifically, the study focused to describe how varying anesthesia conditions affect the resting state activities and resultant FC between bilateral hemispheres in comparison to those obtained by evoked responses. As results, isoflurane decreased the amplitude of evoked responses in a dose-dependent manner mostly due to the habituation of repetitive responses. However, α-chloralose rather intensified the amplitude without exhibiting habituation. No such diverging trend was observed for the spontaneous activity, in which both anesthetics increased the signal power. For α-chloralose, overall FC was similar to that obtained with the lowest dose of isoflurane at 0.5% while higher doses of isoflurane displayed increased FC. Interestingly, only α-chloralose elicited relatively much greater increases in the ipsi-stimulus evoked response (i.e., in S1fl ipsilateral to the stimulated forelimb) than those associated with the contra-stimulus response, suggesting enhanced neuronal excitability. Taken together, the findings demonstrate modulation of the FC profiles by anesthesia is highly non-linear, possibly with a distinct underlying mechanism that affects either resting state or evoked activities differently. Further, the current study warrants thorough investigation of the basal neuronal states prior to the interpretation of resting state FC and evoked activities for accurate understanding of neural signal processing and circuitry.

On-scalp magnetocorticography with optically pumped magnetometers: Simulated performance in resolving simultaneous sources

Currently, the gold standard for high-resolution mapping of cortical electrophysiological activity is invasive electrocorticography (ECoG), a procedure that carries with it the risk of serious morbidity and mortality. Due to these risks, the use of ECoG is largely limited to pre-surgical mapping in intractable epilepsy. Nevertheless, many seminal studies in neuroscience have utilized ECoG to explore domains such as visual perception, attention, auditory processing, and sensorimotor behavior. Studies such as these, occurring in patients with epilepsy rather than healthy controls, may lack generalizability, and are limited by the placement of the electrode arrays over the presumed seizure focus. This manuscript explores the use of optically pumped magnetometers (OPMs) to create a non-invasive alternative to ECoG, which we refer to as magnetocorticography. Because prior ECoG studies reveal that most cognitive processes are driven by multiple, simultaneous independent neuronal assemblies, we characterize the ability of a theoretical 56-channel dense OPM array to resolve simultaneous independent sources, and compare it to currently available SQUID devices, as well as OPM arrays with inter-sensor spacings more typical of other systems in development. Our evaluation of this theoretical system assesses many potential sources of error, including errors of sensor calibration and position. In addition, we investigate the influence of geometrical and anatomical factors on array performance. Our simulations reveal the potential of high-density, on-scalp OPM MEG devices to localize electrophysiological brain responses at unprecedented resolution for a non-invasive device.

Modelling neuron fiber interaction and coupling in non-myelinated bundled fiber

Understanding the local dynamics of a neural network relies heavily on local field potential and cell-field interaction. But it is still unclear how local the local potential is and what kinds of consequences the trans-membrane current flow and produced electric field have on the local neural fiber. Mimicking signal transmission in neighboring nerve fiber, a simulation model is built to analyze local behavior due to trans-membrane current, cell-field interactions, and their repercussions on the bundled fiber system. Simulation studies reveal that depending on the coupling parameters, activity in one fiber can depolarize or hyper-polarize adjacent fibers. The suggested cell-field interaction model was tested using an orientation-selective coupled retinal ganglion cell network, which was compared to its uncoupled counterpart. The proposed work has been used to model and simulate local signal dynamics in a bundled fiber system of an orientation-selective RGC network due to cell-field interaction, as well as to gain insight into the possible significance of dendritic fiber coupling in orientation selectivity bandwidth adjustment.

1024-channel electrophysiological recordings in macaque V1 and V4 during resting state

Co-variations in resting state activity are thought to arise from a variety of correlated inputs to neurons, such as bottom-up activity from lower areas, feedback from higher areas, recurrent processing in local circuits, and fluctuations in neuromodulatory systems. Most studies have examined resting state activity throughout the brain using MRI scans, or observed local co-variations in activity by recording from a small number of electrodes. We carried out electrophysiological recordings from over a thousand chronically implanted electrodes in the visual cortex of non-human primates, yielding a resting state dataset with unprecedentedly high channel counts and spatiotemporal resolution. Such signals could be used to observe brain waves across larger regions of cortex, offering a temporally detailed picture of brain activity. In this paper, we provide the dataset, describe the raw and processed data formats and data acquisition methods, and indicate how the data can be used to yield new insights into the ‘background’ activity that influences the processing of visual information in our brain.

Measurement(s)	brain activity measurement • eye movement
Technology Type(s)	extracellular electrophysiology recording • multi-electrode array system • infrared camera
Factor Type(s)	visual stimulation
Sample Characteristic - Organism	Macaca mulatta

Laminar microcircuitry of visual cortex producing attention-associated electric fields

Cognitive operations are widely studied by measuring electric fields through EEG and ECoG. However, despite their widespread use, the neural circuitry giving rise to these signals remains unknown because the functional architecture of cortical columns producing attention-associated electric fields has not been explored. Here, we detail the laminar cortical circuitry underlying an attention-associated electric field measured over posterior regions of the brain in humans and monkeys. First, we identified visual cortical area V4 as one plausible contributor to this attention-associated electric field through inverse modeling of cranial EEG in macaque monkeys performing a visual attention task. Next, we performed laminar neurophysiological recordings on the prelunate gyrus and identified the electric-field-producing dipoles as synaptic activity in distinct cortical layers of area V4. Specifically, activation in the extragranular layers of cortex resulted in the generation of the attention-associated dipole. Feature selectivity of a given cortical column determined the overall contribution to this electric field. Columns selective for the attended feature contributed more to the electric field than columns selective for a different feature. Last, the laminar profile of synaptic activity generated by V4 was sufficient to produce an attention-associated signal measurable outside of the column. These findings suggest that the top-down recipient cortical layers produce an attention-associated electric field that can be measured extracortically with the relative contribution of each column depending upon the underlying functional architecture.

Stimulus dependence of directed information exchange between cortical layers in macaque V1

Perception and cognition require the integration of feedforward sensory information with feedback signals. Using different sized stimuli, we isolate spectral signatures of feedforward and feedback signals, and their effect on communication between layers in primary visual cortex of male macaque monkeys. Small stimuli elicited gamma frequency oscillations predominantly in the superficial layers. These Granger-causally originated in upper layer 4 and lower supragranular layers. Unexpectedly, large stimuli generated strong narrow band gamma oscillatory activity across cortical layers. They Granger-causally arose in layer 5, were conveyed through layer six to superficial layers, and violated existing models of feedback spectral signatures. Equally surprising, with large stimuli, alpha band oscillatory activity arose predominantly in granular and supragranular layers and communicated in a feedforward direction. Thus, oscillations in specific frequency bands are dynamically modulated to serve feedback and feedforward communication and are not restricted to specific cortical layers in V1.

Operating principles of the cerebral cortex as a six-layered network in primates: beyond the classic canonical circuit model

The cerebral cortex performs its computations with many six-layered fundamental units, collectively spreading along the cortical sheet. What is the local network structure and the operating dynamics of such a fundamental unit? Previous investigations of primary sensory areas revealed a classic "canonical" circuit model, leading to an expectation of similar circuit organization and dynamics throughout the cortex. This review clarifies the different circuit dynamics at play in the higher association cortex of primates that implements computation for high-level cognition such as memory and attention. Instead of feedforward processing of response selectivity through Layers 4 to 2/3 that the classic canonical circuit stipulates, memory recall in primates occurs in Layer 5/6 with local backward projection to Layer 2/3, after which the retrieved information is sent back from Layer 6 to lower-level cortical areas for further retrieval of nested associations of target attributes. In this review, a novel "dynamic multimode module (D3M)" in the primate association cortex is proposed, as a new "canonical" circuit model performing this operation.

Pop-out search instigates beta-gated feature selectivity enhancement across V4 layers

Visual search is a workhorse for investigating how attention interacts with processing of sensory information. Attentional selection has been linked to altered cortical sensory responses and feature preferences (i.e., tuning). However, attentional modulation of feature selectivity during search is largely unexplored. Here we map the spatiotemporal profile of feature selectivity during singleton search. Monkeys performed a search where a pop-out feature determined the target of attention. We recorded laminar neural responses from visual area V4. We first identified "feature columns" which showed preference for individual colors. In the unattended condition, feature columns were significantly more selective in superficial relative to middle and deep layers. Attending a stimulus increased selectivity in all layers but not equally. Feature selectivity increased most in the deep layers, leading to higher selectivity in extragranular layers as compared to the middle layer. This attention-induced enhancement was rhythmically gated in phase with the beta-band local field potential. Beta power dominated both extragranular laminar compartments, but current source density analysis pointed to an origin in superficial layers, specifically. While beta-band power was present regardless of attentional state, feature selectivity was only gated by beta in the attended condition. Neither the beta oscillation nor its gating of feature selectivity varied with microsaccade production. Importantly, beta modulation of neural activity predicted response times, suggesting a direct link between attentional gating and behavioral output. Together, these findings suggest beta-range synaptic activation in V4's superficial layers rhythmically gates attentional enhancement of feature tuning in a way that affects the speed of attentional selection.

Inferring entire spiking activity from local field potentials

Extracellular recordings are typically analysed by separating them into two distinct signals: local field potentials (LFPs) and spikes. Previous studies have shown that spikes, in the form of single-unit activity (SUA) or multiunit activity (MUA), can be inferred solely from LFPs with moderately good accuracy. SUA and MUA are typically extracted via threshold-based technique which may not be reliable when the recordings exhibit a low signal-to-noise ratio (SNR). Another type of spiking activity, referred to as entire spiking activity (ESA), can be extracted by a threshold-less, fast, and automated technique and has led to better performance in several tasks. However, its relationship with the LFPs has not been investigated. In this study, we aim to address this issue by inferring ESA from LFPs intracortically recorded from the motor cortex area of three monkeys performing different tasks. Results from long-term recording sessions and across subjects revealed that ESA can be inferred from LFPs with good accuracy. On average, the inference performance of ESA was consistently and significantly higher than those of SUA and MUA. In addition, local motor potential (LMP) was found to be the most predictive feature. The overall results indicate that LFPs contain substantial information about spiking activity, particularly ESA. This could be useful for understanding LFP-spike relationship and for the development of LFP-based BMIs.

Cross Laminar Traveling Components of Field Potentials due to Volume Conduction of Non-Traveling Neuronal Activity in Macaque Sensory Cortices

Field potentials (FPs) reflect neuronal activities in the brain, and often exhibit traveling peaks across recording sites. While traveling FPs are interpreted as propagation of neuronal activity, not all studies directly reveal such propagating patterns of neuronal activation. Neuronal activity is associated with transmembrane currents that form dipoles and produce negative and positive fields. Thereby, FP components reverse polarity between those fields and have minimal amplitudes at the center of dipoles. Although their amplitudes could be smaller, FPs are never flat even around these reversals. What occurs around the reversal has not been addressed explicitly, although those are rationally in the middle of active neurons. We show that sensory FPs around the reversal appeared with peaks traveling across cortical laminae in macaque sensory cortices. Interestingly, analyses of current source density did not depict traveling patterns but lamina-delimited current sinks and sources. We simulated FPs produced by volume conduction of a simplified 2 dipoles' model mimicking sensory cortical laminar current source density components. While FPs generated by single dipoles followed the temporal patterns of the dipole moments without traveling peaks, FPs generated by concurrently active dipole moments appeared with traveling components in the vicinity of dipoles by superimposition of individually non-traveling FPs generated by single dipoles. These results indicate that not all traveling FP are generated by traveling neuronal activity, and that recording positions need to be taken into account to describe FP peak components around active neuronal populations.SIGNIFICANCE STATEMENT Field potentials (FPs) generated by neuronal activity in the brain occur with fields of opposite polarity. Likewise, in the cerebral cortices, they have mirror-imaged waveforms in upper and lower layers. We show that FPs appear like traveling across the cortical layers. Interestingly, the traveling FPs occur without traveling components of current source density, which represents transmembrane currents associated with neuronal activity. These seemingly odd findings are explained using current source density models of multiple dipoles. Concurrently active, non-traveling dipoles produce FPs as mixtures of FPs produced by individual dipoles, and result in traveling FP waveforms as the mixing ratio depends on the distances from those dipoles. The results suggest that not all traveling FP components are associated with propagating neuronal activity.

Directed information exchange between cortical layers in macaque V1 and V4 and its modulation by selective attention

Achieving behavioral goals requires integration of sensory and cognitive information across cortical laminae and cortical regions. How this computation is performed remains unknown. Using local field potential recordings and spectrally resolved conditional Granger causality (cGC) analysis, we mapped visual information flow, and its attentional modulation, between cortical layers within and between macaque brain areas V1 and V4. Stimulus-induced interlaminar information flow within V1 dominated upwardly, channeling information toward supragranular corticocortical output layers. Within V4, information flow dominated from granular to supragranular layers, but interactions between supragranular and infragranular layers dominated downwardly. Low-frequency across-area communication was stronger from V4 to V1, with little layer specificity. Gamma-band communication was stronger in the feedforward V1-to-V4 direction. Attention to the receptive field of V1 decreased communication between all V1 layers, except for granular-to-supragranular layer interactions. Communication within V4, and from V1 to V4, increased with attention across all frequencies. While communication from V4 to V1 was stronger in lower-frequency bands (4 to 25 Hz), attention modulated cGCs from V4 to V1 across all investigated frequencies. Our data show that top-down cognitive processes result in reduced communication within cortical areas, increased feedforward communication across all frequency bands, and increased gamma-band feedback communication.

Theta, but Not Gamma Oscillations in Area V4 Depend on Input from Primary Visual Cortex

Theta (3-9 Hz) and gamma (30-100 Hz) oscillations have been observed at different levels along the hierarchy of cortical areas and across a wide set of cognitive tasks. In the visual system, the emergence of both rhythms in primary visual cortex (V1) and mid-level cortical areas V4 has been linked with variations in perceptual reaction times.1-5 Based on analytical methods to infer causality in neural activation patterns, it was concluded that gamma and theta oscillations might both reflect feedforward sensory processing from V1 to V4.6-10 Here, we report on experiments in macaque monkeys in which we experimentally assessed the presence of both oscillations in the neural activity recorded from multi-electrode arrays in V1 and V4 before and after a permanent V1 lesion. With intact cortex, theta and gamma oscillations could be reliably elicited in V1 and V4 when monkeys viewed a visual contour illusion and showed phase-to-amplitude coupling. Laminar analysis in V1 revealed that both theta and gamma oscillations occurred primarily in the supragranular layers, the cortical output compartment of V1. However, there was a clear dissociation between the two rhythms in V4 that became apparent when the major feedforward input to V4 was removed by lesioning V1: although V1 lesioning eliminated V4 theta, it had little effect on V4 gamma power except for delaying its emergence by >100 ms. These findings suggest that theta is more tightly associated with feedforward processing than gamma and pose limits on the proposed role of gamma as a feedforward mechanism.

Decision Signals in the Local Field Potentials of Early and Mid-Level Macaque Visual Cortex

The neural basis of perceptual decision making has typically been studied using measurements of single neuron activity, though decisions are likely based on the activity of large neuronal ensembles. Local field potentials (LFPs) may, in some cases, serve as a useful proxy for population activity and thus be useful for understanding the neural basis of perceptual decision making. However, little is known about whether LFPs in sensory areas include decision-related signals. We therefore analyzed LFPs recorded using two 48-electrode arrays implanted in primary visual cortex (V1) and area V4 of macaque monkeys trained to perform a fine orientation discrimination task. We found significant choice information in low (0-30 Hz) and higher (70-500 Hz) frequency components of the LFP, but little information in gamma frequencies (30-70 Hz). Choice information was more robust in V4 than V1 and stronger in LFPs than in simultaneously measured spiking activity. LFP-based choice information included a global component, common across electrodes within an area. Our findings reveal the presence of robust choice-related signals in the LFPs recorded in V1 and V4 and suggest that LFPs may be a useful complement to spike-based analyses of decision making.

Stimulus Feature-Specific Information Flow Along the Columnar Cortical Microcircuit Revealed by Multivariate Laminar Spiking Analysis

Most of the mammalian neocortex is comprised of a highly similar anatomical structure, consisting of a granular cell layer between superficial and deep layers. Even so, different cortical areas process different information. Taken together, this suggests that cortex features a canonical functional microcircuit that supports region-specific information processing. For example, the primate primary visual cortex (V1) combines the two eyes' signals, extracts stimulus orientation, and integrates contextual information such as visual stimulation history. These processes co-occur during the same laminar stimulation sequence that is triggered by the onset of visual stimuli. Yet, we still know little regarding the laminar processing differences that are specific to each of these types of stimulus information. Univariate analysis techniques have provided great insight by examining one electrode at a time or by studying average responses across multiple electrodes. Here we focus on multivariate statistics to examine response patterns across electrodes instead. Specifically, we applied multivariate pattern analysis (MVPA) to linear multielectrode array recordings of laminar spiking responses to decode information regarding the eye-of-origin, stimulus orientation, and stimulus repetition. MVPA differs from conventional univariate approaches in that it examines patterns of neural activity across simultaneously recorded electrode sites. We were curious whether this added dimensionality could reveal neural processes on the population level that are challenging to detect when measuring brain activity without the context of neighboring recording sites. We found that eye-of-origin information was decodable for the entire duration of stimulus presentation, but diminished in the deepest layers of V1. Conversely, orientation information was transient and equally pronounced along all layers. More importantly, using time-resolved MVPA, we were able to evaluate laminar response properties beyond those yielded by univariate analyses. Specifically, we performed a time generalization analysis by training a classifier at one point of the neural response and testing its performance throughout the remaining period of stimulation. Using this technique, we demonstrate repeating (reverberating) patterns of neural activity that have not previously been observed using standard univariate approaches.

Layer and rhythm specificity for predictive routing

An established theoretical model, predictive coding, states that the brain is constantly building models (signifying changing predictions) of the environment. The brain does this by forming predictions and signaling sensory inputs which deviate from predictions (“prediction errors”). Various hypotheses exist about how predictive coding could be implemented in the brain. We recorded neural spiking and oscillations with laminar resolution in a network of cortical areas as monkeys performed a working memory task with changing stimulus predictability. Predictability modulated the patterns of feedforward/feedback flow, cortical layers, and oscillations used to process a visual stimulus. These data support the theory of predictive coding but suggest an alternate model for its neural implementation: predictive routing.

In predictive coding, experience generates predictions that attenuate the feeding forward of predicted stimuli while passing forward unpredicted “errors.” Different models have suggested distinct cortical layers, and rhythms implement predictive coding. We recorded spikes and local field potentials from laminar electrodes in five cortical areas (visual area 4 [V4], lateral intraparietal [LIP], posterior parietal area 7A, frontal eye field [FEF], and prefrontal cortex [PFC]) while monkeys performed a task that modulated visual stimulus predictability. During predictable blocks, there was enhanced alpha (8 to 14 Hz) or beta (15 to 30 Hz) power in all areas during stimulus processing and prestimulus beta (15 to 30 Hz) functional connectivity in deep layers of PFC to the other areas. Unpredictable stimuli were associated with increases in spiking and in gamma-band (40 to 90 Hz) power/connectivity that fed forward up the cortical hierarchy via superficial-layer cortex. Power and spiking modulation by predictability was stimulus specific. Alpha/beta power in LIP, FEF, and PFC inhibited spiking in deep layers of V4. Area 7A uniquely showed increases in high-beta (∼22 to 28 Hz) power/connectivity to unpredictable stimuli. These results motivate a conceptual model, predictive routing. It suggests that predictive coding may be implemented via lower-frequency alpha/beta rhythms that “prepare” pathways processing-predicted inputs by inhibiting feedforward gamma rhythms and associated spiking.

Dissociation of broadband high-frequency activity and neuronal firing in the neocortex

Broadband high-frequency activity (BHA; 70 to 150 Hz), also known as "high gamma," a key analytic signal in human intracranial (electrocorticographic) recordings, is often assumed to reflect local neural firing [multiunit activity (MUA)]. As the precise physiological substrates of BHA are unknown, this assumption remains controversial. Our analysis of laminar multielectrode data from V1 and A1 in monkeys outlines two components of stimulus-evoked BHA distributed across the cortical layers: an "early-deep" and "late-superficial" response. Early-deep BHA has a clear spatial and temporal overlap with MUA. Late-superficial BHA was more prominent and accounted for more of the BHA signal measured near the cortical pial surface. However, its association with local MUA is weak and often undetectable, consistent with the view that it reflects dendritic processes separable from local neuronal firing.

Cuprizone-induced oligodendrocyte loss and demyelination impairs recording performance of chronically implanted neural interfaces

Biological inflammation induced during penetrating cortical injury can disrupt functional neuronal and glial activity within the cortex, resulting in potential recording failure of chronically implanted neural interfaces. Oligodendrocytes provide critical support for neuronal health and function through direct contact with neuronal soma and axons within the cortex. Given their fundamental role to regulate neuronal activity via myelin, coupled with their heightened vulnerability to metabolic brain injury due to high energetic demands, oligodendrocytes are hypothesized as a possible source of biological failure in declining recording performances of intracortical microelectrode devices. To determine the extent of their contribution to neuronal activity and function, a cuprizone-inducible model of oligodendrocyte depletion and demyelination in mice was performed prior to microelectrode implantation. At 5 weeks of cuprizone exposure, mice demonstrated significantly reduced cortical oligodendrocyte density and myelin expression. Mice were then implanted with functional recording microelectrodes in the visual cortex and neuronal activity was evaluated up to 7 weeks alongside continued cuprizone administration. Cuprizone-induced oligodendrocyte loss and demyelination was associated with significantly reduced recording performances at the onset of implantation, which remained relatively stable over time. In contast, recording performances for mice on a normal diet were intially elevated before decreasing over time to the recording level of tcuprizone-treated mice. Further electrophysiological analysis revealed deficits in multi-unit firing rates, frequency-dependent disruptions in neuronal oscillations, and altered laminar communication within the cortex of cuprizone-treated mice. Post-mortem immunohistochemistry revealed robust depletion of oligodendrocytes around implanted microelectrode arrays alongside comparable neuronal densities to control mice, suggesting that oligodendrocyte loss was a possible contributor to chronically impaired device performances. This study highlights potentially significant contributions from the oligodendrocyte lineage population concerning the biological integration and long-term functional performance of neural interfacing technology.

Estimates of cortical column orientation improve MEG source inversion

Determining the anatomical source of brain activity non-invasively measured from EEG or MEG sensors is challenging. In order to simplify the source localization problem, many techniques introduce the assumption that current sources lie on the cortical surface. Another common assumption is that this current flow is orthogonal to the cortical surface, thereby approximating the orientation of cortical columns. However, it is not clear which cortical surface to use to define the current source locations, and normal vectors computed from a single cortical surface may not be the best approximation to the orientation of cortical columns. We compared three different surface location priors and five different approaches for estimating dipole vector orientation, both in simulations and visual and motor evoked MEG responses. We show that models with source locations on the white matter surface and using methods based on establishing correspondences between white matter and pial cortical surfaces dramatically outperform models with source locations on the pial or combined pial/white surfaces and which use methods based on the geometry of a single cortical surface in fitting evoked visual and motor responses. These methods can be easily implemented and adopted in most M/EEG analysis pipelines, with the potential to significantly improve source localization of evoked responses.

Cortical Responses to Input From Distant Areas are Modulated by Local Spontaneous Alpha/Beta Oscillations

Any given area in human cortex may receive input from multiple, functionally heterogeneous areas, potentially representing different processing threads. Alpha (8-13 Hz) and beta oscillations (13-20 Hz) have been hypothesized by other investigators to gate local cortical processing, but their influence on cortical responses to input from other cortical areas is unknown. To study this, we measured the effect of local oscillatory power and phase on cortical responses elicited by single-pulse electrical stimulation (SPES) at distant cortical sites, in awake human subjects implanted with intracranial electrodes for epilepsy surgery. In 4 out of 5 subjects, the amplitudes of corticocortical evoked potentials (CCEPs) elicited by distant SPES were reproducibly modulated by the power, but not the phase, of local oscillations in alpha and beta frequencies. Specifically, CCEP amplitudes were higher when average oscillatory power just before distant SPES (-110 to -10 ms) was high. This effect was observed in only a subset (0-33%) of sites with CCEPs and, like the CCEPs themselves, varied with stimulation at different distant sites. Our results suggest that although alpha and beta oscillations may gate local processing, they may also enhance the responsiveness of cortex to input from distant cortical sites.

Brain-wave equation incorporating axodendritic connectivity

We introduce an integral model of a two-dimensional neural field that includes a third dimension representing space along a dendritic tree that can incorporate realistic patterns of axodendritic connectivity. For natural choices of this connectivity we show how to construct an equivalent brain-wave partial differential equation that allows for efficient numerical simulation of the model. This is used to highlight the effects that passive dendritic properties can have on the speed and shape of large scale traveling cortical waves.

Temporal dynamics of binocular integration in primary visual cortex

Whenever we open our eyes, our brain quickly integrates the two eyes' perspectives into a combined view. This process of binocular integration happens so rapidly that even incompatible stimuli are briefly fused before one eye's view is suppressed in favor of the other (binocular rivalry). The neuronal basis for this brief period of fusion during incompatible binocular stimulation is unclear. Neuroanatomically, the eyes provide two largely separate streams of information that are integrated into a binocular response by the primary visual cortex (V1). However, the temporal dynamics underlying the formation of this binocular response are largely unknown. To address this question, we examined the temporal profile of binocular responses in V1 of fixating monkeys. We found that V1 processes binocular stimuli in a dynamic sequence that comprises at least two distinct temporal phases. An initial transient phase is characterized by enhanced spiking responses for both compatible and incompatible binocular stimuli compared to monocular stimulation. This transient is followed by a sustained response that differed markedly between congruent and incongruent binocular stimulation. Specifically, incompatible binocular stimulation resulted in overall response reduction relative to monocular stimulation (binocular suppression). In contrast, responses to compatible stimuli were either suppressed or enhanced (binocular facilitation) depending on the neurons' ocularity (selectivity for one eye over the other) and laminar location. These results suggest that binocular integration in V1 occurs in at least two sequential steps that comprise initial additive combination of the two eyes' signals followed by widespread differentiation between binocular concordance and discordance.

A Sound-Sensitive Source of Alpha Oscillations in Human Non-Primary Auditory Cortex

The functional organization of human auditory cortex can be probed by characterizing responses to various classes of sound at different anatomical locations. Along with histological studies this approach has revealed a primary field in posteromedial Heschl's gyrus (HG) with pronounced induced high-frequency (70–150 Hz) activity and short-latency responses that phase-lock to rapid transient sounds. Low-frequency neural oscillations are also relevant to stimulus processing and information flow, however, their distribution within auditory cortex has not been established. Alpha activity (7–14 Hz) in particular has been associated with processes that may differentially engage earlier versus later levels of the cortical hierarchy, including functional inhibition and the communication of sensory predictions. These theories derive largely from the study of occipitoparietal sources readily detectable in scalp electroencephalography. To characterize the anatomical basis and functional significance of less accessible temporal-lobe alpha activity we analyzed responses to sentences in seven human adults (4 female) with epilepsy who had been implanted with electrodes in superior temporal cortex. In contrast to primary cortex in posteromedial HG, a non-primary field in anterolateral HG was characterized by high spontaneous alpha activity that was strongly suppressed during auditory stimulation. Alpha-power suppression decreased with distance from anterolateral HG throughout superior temporal cortex, and was more pronounced for clear compared to degraded speech. This suppression could not be accounted for solely by a change in the slope of the power spectrum. The differential manifestation and stimulus-sensitivity of alpha oscillations across auditory fields should be accounted for in theories of their generation and function.

SIGNIFICANCE STATEMENT To understand how auditory cortex is organized in support of perception, we recorded from patients implanted with electrodes for clinical reasons. This allowed measurement of activity in brain regions at different levels of sensory processing. Oscillations in the alpha range (7–14 Hz) have been associated with functions including sensory prediction and inhibition of regions handling irrelevant information, but their distribution within auditory cortex is not known. A key finding was that these oscillations dominated in one particular non-primary field, anterolateral Heschl's gyrus, and were suppressed when subjects listened to sentences. These results build on our knowledge of the functional organization of auditory cortex and provide anatomical constraints on theories of the generation and function of alpha oscillations.

V1 microcircuit dynamics: altered signal propagation suggests intracortical origins for adaptation in response to visual repetition

Repetitive visual stimulation profoundly changes sensory processing in the primary visual cortex (V1). We show how the associated adaptive changes are linked to an altered flow of synaptic activation across the V1 laminar microcircuit. Using repeated visual stimulation, we recorded layer-specific responses in V1 of two fixating monkeys. We found that repetition-related spiking suppression was most pronounced outside granular V1 layers that receive the main retinogeniculate input. This repetition-related response suppression was robust to alternating stimuli between the eyes, in line with the notion that repetition-related adaptation is predominantly of cortical origin. Most importantly, current source density (CSD) analysis, which provides an estimate of local net depolarization, revealed that synaptic processing during repeated stimulation was most profoundly affected within supragranular layers, which harbor the bulk of cortico-cortical connections. Direct comparison of the temporal evolution of laminar CSD and spiking activity showed that stimulus repetition first affected supragranular synaptic currents, which translated into a reduction of stimulus-evoked spiking across layers. Together, these results suggest that repetition induces an altered state of intracortical processing that underpins visual adaptation. NEW & NOTEWORTHY Our survival depends on our brains rapidly adapting to ever changing environments. A well-studied form of adaptation occurs whenever we encounter the same or similar stimuli repeatedly. We show that this repetition-related adaptation is supported by systematic changes in the flow of sensory activation across the laminar cortical microcircuitry of primary visual cortex. These results demonstrate how adaptation impacts neuronal interactions across cortical circuits.

Phase shifts in high-beta- and low-gamma-band local field potentials predict the focus of visual spatial attention

The local field potential (LFP) contains rich information about activity in local neuronal populations. However, it has been challenging to establish direct links between LFP modulations and task-relevant behavior or cognitive processes, such as attention. We sought to determine whether LFP amplitude or phase modulations are predictive of the allocation of visual spatial attention. LFPs were recorded simultaneously in multiple early visual brain structures of alert macaque monkeys performing attention-demanding detection and discrimination tasks. Attention directed toward the receptive field of recorded neurons generated systematically larger phase shifts in high-beta- and low-gamma-frequency LFPs compared with LFP phase shifts on trials in which attention was directed away from the receptive field. This attention-mediated temporal advance corresponded to ~10 ms. LFP phase shifts also correlated with reaction times when monkeys were engaged in the tasks. Importantly, attentional modulation of LFP phase was consistent across monkeys, tasks, visual brain structures, and cortical layers. In contrast, attentional modulation of LFP amplitude varied across frequency bands, visual structures/layers, and tasks. Because LFP phase shifts were robust, consistent, and predictive of spatial attention, they could serve as a reliable marker for attention signals in the brain. NEW & NOTEWORTHY Local field potentials (LFPs) reflect the activity of spatially localized populations of neurons. Whether alterations in LFP activity are indicative of cognitive processes, such as attention, is unclear. We found that shifts in the phase of LFPs measured in multiple visual brain areas reliably predicted the focus of spatial attention. LFP phase shifts could therefore serve as a marker for behaviorally relevant attention signals in the brain.

Toward a Neurobiologically Plausible Model of Language-Related, Negative Event-Related Potentials

Language-related event-related potential (ERP) components such as the N400 have traditionally been associated with linguistic or cognitive functional interpretations. By contrast, it has been considerably more difficult to relate these components to neurobiologically grounded accounts of language. Here, we propose a theoretical framework based on a predictive coding architecture, within which negative language-related ERP components such as the N400 can be accounted for in a neurobiologically plausible manner. Specifically, we posit that the amplitude of negative language-related ERP components reflects precision-weighted prediction error signals, i.e., prediction errors weighted by the relevance of the information source leading to the error. From this perspective, precision has a direct link to cue validity in a particular language and, thereby, to relevance of individual linguistic features for internal model updating. We view components such as the N400 and LAN as members of a family with similar functional characteristics and suggest that latency and topography differences between these components reflect the locus of prediction errors and model updating within a hierarchically organized cortical predictive coding architecture. This account has the potential to unify findings from the full range of the N400 literature, including word-level, sentence-, and discourse-level results as well as cross-linguistic differences.

Correlation Structure in Micro-ECoG Recordings is Described by Spatially Coherent Components

Electrocorticography (ECoG) is becoming more prevalent due to improvements in fabrication and recording technology as well as its ease of implantation compared to intracortical electrophysiology, larger cortical coverage, and potential advantages for use in long term chronic implantation. Given the flexibility in the design of ECoG grids, which is only increasing, it remains an open question what geometry of the electrodes is optimal for an application. Conductive polymer, PEDOT:PSS, coated microelectrodes have an advantage that they can be made very small without losing low impedance. This makes them suitable for evaluating the required granularity of ECoG recording in humans and experimental animals. We used two-dimensional (2D) micro-ECoG grids to record intra-operatively in humans and during acute implantations in mouse with separation distance between neighboring electrodes (i.e., pitch) of 0.4 mm and 0.2/0.25 mm respectively. To assess the spatial properties of the signals, we used the average correlation between electrodes as a function of the pitch. In agreement with prior studies, we find a strong frequency dependence in the spatial scale of correlation. By applying independent component analysis (ICA), we find that the spatial pattern of correlation is largely due to contributions from multiple spatially extended, time-locked sources present at any given time. Our analysis indicates the presence of spatially structured activity down to the sub-millimeter spatial scale in ECoG despite the effects of volume conduction, justifying the use of dense micro-ECoG grids.

Electrocorticography (ECoG) is a type of electrophysiological monitoring that uses electrodes placed directly on the exposed surface of the brain. ECoG is a promising technique for studying the brain, and EcoG signals can be used to control brain-computer interfaces. Advances have made it possible to record simultaneously with an increasing number of smaller, and more closely spaced electrodes. However, a property of electrical recording from outside the brain is that common signals appear on different electrodes at different locations, and this affects decisions about how to best distribute a limited number of electrodes to maximize the information that can be gathered. Large spacing of electrodes around one centimeter apart on the brain’s surface has proven useful for clinical and research use, but how much benefit there is to recording from more locations in a smaller area remains to be answered. We found that we can explain the commonality between the different locations as the combination of different patterns of brain activity that are present at multiple electrode locations, and that signals recorded from very closely spaced electrodes, around a millimeter or less apart, are able to identify patterns that are at this small scale.

Layer-Specific Physiological Features and Interlaminar Interactions in the Primary Visual Cortex of the Mouse

The relationship between mesoscopic local field potentials (LFPs) and single-neuron firing in the multi-layered neocortex is poorly understood. Simultaneous recordings from all layers in the primary visual cortex (V1) of the behaving mouse revealed functionally defined layers in V1. The depth of maximum spike power and sink-source distributions of LFPs provided consistent laminar landmarks across animals. Coherence of gamma oscillations (30-100 Hz) and spike-LFP coupling identified six physiological layers and further sublayers. Firing rates, burstiness, and other electrophysiological features of neurons displayed unique layer and brain state dependence. Spike transmission strength from layer 2/3 cells to layer 5 pyramidal cells and interneurons was stronger during waking compared with non-REM sleep but stronger during non-REM sleep among deep-layer excitatory neurons. A subset of deep-layer neurons was active exclusively in the DOWN state of non-REM sleep. These results bridge mesoscopic LFPs and single-neuron interactions with laminar structure in V1.

Binocular Modulation of Monocular V1 Neurons

In humans and other primates, sensory signals from each eye remain separated until they arrive in the primary visual cortex (V1), but their exact meeting point is unknown. In V1, some neurons respond to stimulation of only one eye (monocular neurons), while most neurons respond to stimulation of either eye (binocular neurons). The main input layers of V1 contain most of the monocular neurons while binocular neurons dominate the layers above and below. This observation has given rise to the idea that the two eyes' signals remain separate until they converge outside V1's input layers. Here, we show that, despite responding to only one eye, monocular neurons in all layers, including the input layers, of V1 discriminate between stimulation of their driving eye alone and stimulation of both eyes. Some monocular V1 neurons' responses were significantly enhanced, or facilitated, when both eyes were stimulated. Binocular facilitation within V1's input layers tended to occur at the onset of the visual response, which could be explained by converging thalamocortical inputs. However, most V1 monocular neurons were significantly reduced, or suppressed, to binocular stimulation. In contrast to facilitation, binocular suppression occurred several milliseconds following the onset of the visual response, suggesting that the bulk of binocular modulation involves cortical inhibition. These findings, combined, suggest that binocular signals arise at an earlier processing stage than previously appreciated, as even so-called monocular neurons in V1's input layers encode what is shown to both eyes.

Alpha Oscillations Modulate Preparatory Activity in Marmoset Area 8Ad

Cognitive control often requires suppression of prepotent stimulus-driven responses in favor of less potent alternatives. Suppression of prepotent saccades has been shown to require proactive inhibition in the frontoparietal saccade network. Electrophysiological evidence in macaque monkeys has revealed neural correlates of such inhibition in this network; however, the interlaminar instantiation of inhibitory processes remains poorly understood because these areas lie deep within sulci in macaques, rendering them inaccessible to laminar recordings. Here, we addressed this gap by exploiting the mostly lissencephalic cortex of the common marmoset (Callithrix jacchus). We inserted linear electrode arrays into areas 8Ad-the putative marmoset frontal eye field-and the lateral intraparietal area of two male marmosets and recorded neural activity during performance of a task comprised of alternating blocks of trials requiring a saccade either toward a large, high-luminance stimulus or the inhibition of this prepotent response in favor of a saccade toward a small, low-luminance stimulus. We observed prominent task-dependent activity in both alpha/gamma bands of the LFP and discharge rates of single neurons in area 8Ad during a prestimulus task epoch in which the animals had been instructed which of these two tasks to perform but before peripheral stimulus onset. These data are consistent with a model in which rhythmic alpha-band activity in deeper layers inhibits spiking in upper layers to support proactive inhibitory saccade control.SIGNIFICANCE STATEMENT Failures to inhibit automatic saccadic responses are a hallmark of many neuropsychiatric disorders, but how this process is implemented across the cortical layers in the frontoparietal saccade network remains unknown because many of the areas are inaccessible to laminar recordings in macaques. Here, we investigated laminar neural activity in marmoset monkeys, which have a smooth cortex. Monkeys were required either to generate or inhibit a prepotent saccade response. In area 8Ad, the putative frontal eye field in marmosets, rhythmic alpha-band activity (9-14 Hz) was higher in deeper layers and spiking activity was lower in upper layers when the animals were instructed to suppress a saccade toward a peripheral stimulus. Reduced alpha power during task preparation may be the underlying common neural basis of a saccade suppression deficit.

Functional connectivity with cortical depth assessed by resting state fMRI of subregions of S1 in squirrel monkeys

Whereas resting state blood oxygenation-level dependent (BOLD) functional MRI has been widely used to assess functional connectivity between cortical regions, the laminar specificity of such measures is poorly understood. This study aims to determine: (a) whether the resting state functional connectivity (rsFC) between two functionally related cortical regions varies with cortical depth, (b) the relationship between layer-resolved tactile stimulus-evoked activation pattern and interlayer rsFC pattern between two functionally distinct but related somatosensory areas 3b and 1, and (c) the effects of spatial resolution on rsFC measures. We examined the interlayer rsFC between areas 3b and 1 of squirrel monkeys under anesthesia using tactile stimulus-driven and resting state BOLD acquisitions at submillimeter resolution. Consistent with previous observations in the areas 3b and 1, we detected robust stimulus-evoked BOLD activations with foci were confined mainly to the upper layers (centered at 21% of the cortical depth). By carefully placing seeds in upper, middle, and lower layers of areas 3b and 1, we observed strong rsFC between upper and middle layers of these two areas. The layer-resolved activation patterns in areas 3b and 1 agree with their interlayer rsFC patterns, and are consistent with the known anatomical connections between layers. In summary, using BOLD rsFC pattern, we identified an interlayer interareal microcircuit that shows strong intrinsic functional connections between upper and middle layer areas 3b and 1. RsFC can be used as a robust invasive tool to probe interlayer corticocortical microcircuits.

Spiking Suppression Precedes Cued Attentional Enhancement of Neural Responses in Primary Visual Cortex

Attending to a visual stimulus increases its detectability, even if gaze is directed elsewhere. This covert attentional selection is known to enhance spiking across many brain areas, including the primary visual cortex (V1). Here we investigate the temporal dynamics of attention-related spiking changes in V1 of macaques performing a task that separates attentional selection from the onset of visual stimulation. We found that preceding attentional enhancement there was a sharp, transient decline in spiking following presentation of an attention-guiding cue. This disruption of V1 spiking was not observed in a task-naïve subject that passively observed the same stimulus sequence, suggesting that sensory activation is insufficient to cause suppression. Following this suppression, attended stimuli evoked more spiking than unattended stimuli, matching previous reports of attention-related activity in V1. Laminar analyses revealed a distinct pattern of activation in feedback-associated layers during both the cue-induced suppression and subsequent attentional enhancement. These findings suggest that top-down modulation of V1 spiking can be bidirectional and result in either suppression or enhancement of spiking responses.

Multimodal Modeling of Neural Network Activity: Computing LFP, ECoG, EEG, and MEG Signals With LFPy 2.0

Recordings of extracellular electrical, and later also magnetic, brain signals have been the dominant technique for measuring brain activity for decades. The interpretation of such signals is however nontrivial, as the measured signals result from both local and distant neuronal activity. In volume-conductor theory the extracellular potentials can be calculated from a distance-weighted sum of contributions from transmembrane currents of neurons. Given the same transmembrane currents, the contributions to the magnetic field recorded both inside and outside the brain can also be computed. This allows for the development of computational tools implementing forward models grounded in the biophysics underlying electrical and magnetic measurement modalities. LFPy (LFPy.readthedocs.io) incorporated a well-established scheme for predicting extracellular potentials of individual neurons with arbitrary levels of biological detail. It relies on NEURON (neuron.yale.edu) to compute transmembrane currents of multicompartment neurons which is then used in combination with an electrostatic forward model. Its functionality is now extended to allow for modeling of networks of multicompartment neurons with concurrent calculations of extracellular potentials and current dipole moments. The current dipole moments are then, in combination with suitable volume-conductor head models, used to compute non-invasive measures of neuronal activity, like scalp potentials (electroencephalographic recordings; EEG) and magnetic fields outside the head (magnetoencephalographic recordings; MEG). One such built-in head model is the four-sphere head model incorporating the different electric conductivities of brain, cerebrospinal fluid, skull and scalp. We demonstrate the new functionality of the software by constructing a network of biophysically detailed multicompartment neuron models from the Neocortical Microcircuit Collaboration (NMC) Portal (bbp.epfl.ch/nmc-portal) with corresponding statistics of connections and synapses, and compute in vivo-like extracellular potentials (local field potentials, LFP; electrocorticographical signals, ECoG) and corresponding current dipole moments. From the current dipole moments we estimate corresponding EEG and MEG signals using the four-sphere head model. We also show strong scaling performance of LFPy with different numbers of message-passing interface (MPI) processes, and for different network sizes with different density of connections. The open-source software LFPy is equally suitable for execution on laptops and in parallel on high-performance computing (HPC) facilities and is publicly available on GitHub.com.

Isoflurane and ketamine differentially influence spontaneous and evoked laminar electrophysiology in mouse V1

General anesthesia is ubiquitous in research and medicine, yet although the molecular mechanisms of anesthetics are well characterized, their ultimate influence on cortical electrophysiology remains unclear. Moreover, the influence that different anesthetics have on sensory cortexes at neuronal and ensemble scales is mostly unknown and represents an important gap in knowledge that has widespread relevance for neural sciences. To address this knowledge gap, this work explored the effects of isoflurane and ketamine/xylazine, two widely used anesthetic paradigms, on electrophysiological behavior in mouse primary visual cortex. First, multiunit activity and local field potentials were examined to understand how each anesthetic influences spontaneous activity. Then, the interlaminar relationships between populations of neurons at different cortical depths were studied to assess whether anesthetics influenced resting-state functional connectivity. Lastly, the spatiotemporal dynamics of visually evoked multiunit and local field potentials were examined to determine how each anesthetic alters communication of visual information. We found that isoflurane enhanced the rhythmicity of spontaneous ensemble activity at 10-40 Hz, which coincided with large increases in coherence between layer IV with superficial and deep layers. Ketamine preferentially increased local field potential power from 2 to 4 Hz, and the largest increases in coherence were observed between superficial and deep layers. Visually evoked responses across layers were diminished under isoflurane, and enhanced under ketamine anesthesia. These findings demonstrate that isoflurane and ketamine anesthesia differentially impact sensory processing in V1. NEW & NOTEWORTHY We directly compared electrophysiological responses in awake and anesthetized (isoflurane or ketamine) mice. We also proposed a method for quantifying and visualizing highly variable, evoked multiunit activity. Lastly, we observed distinct oscillatory responses to stimulus onset and offset in awake and isoflurane-anesthetized mice.

Working Memory 2.0

Working memory is the fundamental function by which we break free from reflexive input-output reactions to gain control over our own thoughts. It has two types of mechanisms: online maintenance of information and its volitional or executive control. Classic models proposed persistent spiking for maintenance but have not explicitly addressed executive control. We review recent theoretical and empirical studies that suggest updates and additions to the classic model. Synaptic weight changes between sparse bursts of spiking strengthen working memory maintenance. Executive control acts via interplay between network oscillations in gamma (30-100 Hz) in superficial cortical layers (layers 2 and 3) and alpha and beta (10-30 Hz) in deep cortical layers (layers 5 and 6). Deep-layer alpha and beta are associated with top-down information and inhibition. It regulates the flow of bottom-up sensory information associated with superficial layer gamma. We propose that interactions between different rhythms in distinct cortical layers underlie working memory maintenance and its volitional control.

Lamina-specific cortical dynamics in human visual and sensorimotor cortices

Distinct anatomical and spectral channels are thought to play specialized roles in the communication within cortical networks. While activity in the alpha and beta frequency range (7 - 40 Hz) is thought to predominantly originate from infragranular cortical layers conveying feedback-related information, activity in the gamma range (>40 Hz) dominates in supragranular layers communicating feedforward signals. We leveraged high precision MEG to test this proposal, directly and non-invasively, in human participants performing visually cued actions. We found that visual alpha mapped onto deep cortical laminae, whereas visual gamma predominantly occurred more superficially. This lamina-specificity was echoed in movement-related sensorimotor beta and gamma activity. These lamina-specific pre- and post- movement changes in sensorimotor beta and gamma activity suggest a more complex functional role than the proposed feedback and feedforward communication in sensory cortex. Distinct frequency channels thus operate in a lamina-specific manner across cortex, but may fulfill distinct functional roles in sensory and motor processes.

The Effect of Common Signals on Power, Coherence and Granger Causality: Theoretical Review, Simulations, and Empirical Analysis of Fruit Fly LFPs Data

When analyzing neural data it is important to consider the limitations of the particular experimental setup. An enduring issue in the context of electrophysiology is the presence of common signals. For example a non-silent reference electrode adds a common signal across all recorded data and this adversely affects functional and effective connectivity analysis. To address the common signals problem, a number of methods have been proposed, but relatively few detailed investigations have been carried out. As a result, our understanding of how common signals affect neural connectivity estimation is incomplete. For example, little is known about recording preparations involving high spatial-resolution electrodes, used in linear array recordings. We address this gap through a combination of theoretical review, simulations, and empirical analysis of local field potentials recorded from the brains of fruit flies. We demonstrate how a framework that jointly analyzes power, coherence, and quantities based on Granger causality reveals the presence of common signals. We further show that subtracting spatially adjacent signals (bipolar derivations) largely removes the effects of the common signals. However, in some special cases this operation itself introduces a common signal. We also show that Granger causality is adversely affected by common signals and that a quantity referred to as “instantaneous interaction” is increased in the presence of common signals. The theoretical review, simulation, and empirical analysis we present can readily be adapted by others to investigate the nature of the common signals in their data. Our contributions improve our understanding of how common signals affect power, coherence, and Granger causality and will help reduce the misinterpretation of functional and effective connectivity analysis.