Uncovering circuit mechanisms of current sinks and sources with biophysical simulations of primary visual cortex

Uncovering population contributions to the extracellular potential in the mouse visual system using Laminar Population Analysis

The local field potential (LFP), the low-frequency part of the extracellular potential, reflects transmembrane currents in the vicinity of the recording electrode. Thought mainly to stem from currents caused by synaptic input, it provides information about neural activity complementary to that of spikes, the output of neurons. However, the many neural sources contributing to the LFP, and likewise the derived current source density (CSD), can often make it challenging to interpret. Efforts to improve its interpretability have included the application of statistical decomposition tools like principal component analysis (PCA) and independent component analysis (ICA) to disentangle the contributions from different neural sources. However, their underlying assumptions of, respectively, orthogonality and statistical independence are not always valid for the various processes or pathways generating LFP. Here, we expand upon and validate a decomposition algorithm named Laminar Population Analysis (LPA), which is based on physiological rather than statistical assumptions. LPA utilizes the multiunit activity (MUA) and LFP jointly to uncover the contributions of different populations to the LFP. To perform the validation of LPA, we used data simulated with the large-scale, biophysically detailed model of mouse V1 developed by the Allen Institute. We find that LPA can identify laminar positions within V1 and the temporal profiles of laminar population firing rates from the MUA. We also find that LPA can estimate the salient current sinks and sources generated by feedforward input from the lateral geniculate nucleus (LGN), recurrent activity in V1, and feedback input from the lateromedial (LM) area of visual cortex. LPA identifies and distinguishes these contributions with a greater accuracy than the alternative statistical decomposition methods, PCA and ICA. The contributions from different cortical layers within V1 could however not be robustly separated and identified with LPA. This is likely due to substantial synchrony in population firing rates across layers, which may be reduced with other stimulus protocols in the future. Lastly, we also demonstrate the application of LPA on experimentally recorded MUA and LFP from 24 animals in the publicly available Visual Coding dataset. Our results suggest that LPA can be used both as a method to estimate positions of laminar populations and to uncover salient features in LFP/CSD contributions from different populations.

Layer-specific dynamics of local field potentials in monkey V1 during electrical stimulation

The mammalian neocortex, organized into six cellular layers or laminae, forms a cortical network within layers. Layer specific computations are crucial for sensory processing of visual stimuli within primary visual cortex. Laminar recordings of local field potentials (LFPs) are a powerful tool to study neural activity within cortical layers. Electric brain stimulation is widely used in basic neuroscience and in a large range of clinical applications. However, the layer-specific effects of electric stimulation on LFPs remain unclear. To address this gap, we conducted laminar LFP recordings of the primary visual cortex in monkeys while presenting a flash visual stimulus. Simultaneously, we applied a low frequency sinusoidal current to the occipital lobe with offset frequency to the flash stimulus repetition rate. We analyzed the modulation of visual-evoked potentials with respect to the applied phase of the electric stimulation. Our results reveal that only the deeper layers, but not the superficial layers, show phase-dependent changes in LFP components with respect to the applied current. Employing a cortical column model, we demonstrate that these in vivo observations can be explained by phase-dependent changes in the driving force within neurons of deeper layers. Our findings offer crucial insight into the selective modulation of cortical layers through electric stimulation, thus advancing approaches for more targeted neuromodulation.

Reconciliation of weak pairwise spike-train correlations and highly coherent local field potentials across space

Multi-electrode arrays covering several square millimeters of neural tissue provide simultaneous access to population signals such as extracellular potentials and spiking activity of one hundred or more individual neurons. The interpretation of the recorded data calls for multiscale computational models with corresponding spatial dimensions and signal predictions. Multi-layer spiking neuron network models of local cortical circuits covering about $1\,{\text{mm}^{2}}$ have been developed, integrating experimentally obtained neuron-type-specific connectivity data and reproducing features of observed in-vivo spiking statistics. Local field potentials can be computed from the simulated spiking activity. We here extend a local network and local field potential model to an area of $4\times 4\,{\text{mm}^{2}}$, preserving the neuron density and introducing distance-dependent connection probabilities and conduction delays. We find that the upscaling procedure preserves the overall spiking statistics of the original model and reproduces asynchronous irregular spiking across populations and weak pairwise spike-train correlations in agreement with experimental recordings from sensory cortex. Also compatible with experimental observations, the correlation of local field potential signals is strong and decays over a distance of several hundred micrometers. Enhanced spatial coherence in the low-gamma band around $50\,\text{Hz}$ may explain the recent report of an apparent band-pass filter effect in the spatial reach of the local field potential.

Large-Scale Mechanistic Models of Brain Circuits with Biophysically and Morphologically Detailed Neurons

Understanding the brain requires studying its multiscale interactions from molecules to networks. The increasing availability of large-scale datasets detailing brain circuit composition, connectivity, and activity is transforming neuroscience. However, integrating and interpreting this data remains challenging. Concurrently, advances in supercomputing and sophisticated modeling tools now enable the development of highly detailed, large-scale biophysical circuit models. These mechanistic multiscale models offer a method to systematically integrate experimental data, facilitating investigations into brain structure, function, and disease. This review, based on a Society for Neuroscience 2024 MiniSymposium, aims to disseminate recent advances in large-scale mechanistic modeling to the broader community. It highlights (1) examples of current models for various brain regions developed through experimental data integration; (2) their predictive capabilities regarding cellular and circuit mechanisms underlying experimental recordings (e.g., membrane voltage, spikes, local-field potential, electroencephalography/magnetoencephalography) and brain function; and (3) their use in simulating biomarkers for brain diseases like epilepsy, depression, schizophrenia, and Parkinson's, aiding in understanding their biophysical underpinnings and developing novel treatments. The review showcases state-of-the-art models covering hippocampus, somatosensory, visual, motor, auditory cortical, and thalamic circuits across species. These models predict neural activity at multiple scales and provide insights into the biophysical mechanisms underlying sensation, motor behavior, brain signals, neural coding, disease, pharmacological interventions, and neural stimulation. Collaboration with experimental neuroscientists and clinicians is essential for the development and validation of these models, particularly as datasets grow. Hence, this review aims to foster interest in detailed brain circuit models, leading to cross-disciplinary collaborations that accelerate brain research.

Is there a ubiquitous spectrolaminar motif of local field potential power across primate neocortex?

Mendoza-Halliday, Major et al., 2024 ("The Paper")1 advocates a local field potential (LFP)-based approach to functional identification of cortical layers during "laminar" (simultaneous recordings from all cortical layers) multielectrode recordings in nonhuman primates (NHPs). The Paper describes a "ubiquitous spectrolaminar motif" in the primate neocortex: 1) 75-150 Hz power peaks in the supragranular layers, 2) 10-19 Hz power peaks in the infragranular layers and 3) the crossing point of their laminar power gradients identifies Layer 4 (L4). Identification of L4 is critical in general, but especially for The Paper as the "motif" discovery is couched within a framework whose central hypothesis is that gamma activity originates in the supragranular layers and reflects feedforward activity, while alpha-beta activity originates in the infragranular layers and reflects feedback activity. In an impressive scientific effort, The Paper analyzed laminar data from 14 cortical areas in 2 prior macaque studies and compared them to marmoset, mouse, and human data to further bolster the canonical nature of the motif. Identification of such canonical principles of brain operation is clearly a topic of broad scientific interest. Similarly, a reliable online method for L4 identification would be of broad scientific value for the rapidly increasing use of laminar recordings using numerous evolving technologies. Despite The Paper's strengths, and its potential for scientific impact, a series of concerns that are fundamental to the analysis and interpretation of laminar activity profile data in general, and local field potential (LFP) signals in particular, led us to question its conclusions. We thus evaluated the generality of The Paper's methods and findings using new sets of data comprised of stimulus-evoked laminar response profiles from primary and higher-order auditory cortices (A1 and belt cortex), and primary visual cortex (V1). The rationale for using these areas as a test bed for new methods is that their laminar anatomy and physiology have already been extensively characterized by prior studies, and there is general agreement across laboratories on key matters like L4 identification. Our analyses indicate that The Paper's findings do not generalize well to any of these cortical areas. In particular, we find The Paper's methods for L4 identification to be unreliable. Moreover, both methodological and statistical concerns, outlined below and in the supplement, question the stated prevalence of the motif in The Paper's published dataset. After summarizing our findings and related broader concerns, we briefly critique the evidence from biophysical modeling studies cited to support The Paper's conclusions. While our findings are at odds with the proposition of a ubiquitous spectrolaminar motif in the primate neocortex, The Paper already has, and will continue to spark debate and further experimentation. Hopefully this countervailing presentation will lead to robust collegial efforts to define optimal strategies for applying laminar recording methods in future studies.

State-dependent complexity of the local field potential in the primary visual cortex

The local field potential (LFP) is as a measure of the combined activity of neurons within a region of brain tissue. While biophysical modeling schemes for LFP in cortical circuits are well established, there is a paramount lack of understanding regarding the LFP properties along the states assumed in cortical circuits over long periods. Here we use a symbolic information approach to determine the statistical complexity based on Jensen disequilibrium measure and Shannon entropy of LFP data recorded from the primary visual cortex (V1) of urethane-anesthetized rats and freely moving mice. Using these information quantifiers, we find consistent relations between LFP recordings and measures of cortical states at the neuronal level. More specifically, we show that LFP's statistical complexity is sensitive to cortical state (characterized by spiking variability), as well as to cortical layer. In addition, we apply these quantifiers to characterize behavioral states of freely moving mice, where we find indirect relations between such states and spiking variability.

Multitask learning of a biophysically-detailed neuron model

The human brain operates at multiple levels, from molecules to circuits, and understanding these complex processes requires integrated research efforts. Simulating biophysically-detailed neuron models is a computationally expensive but effective method for studying local neural circuits. Recent innovations have shown that artificial neural networks (ANNs) can accurately predict the behavior of these detailed models in terms of spikes, electrical potentials, and optical readouts. While these methods have the potential to accelerate large network simulations by several orders of magnitude compared to conventional differential equation based modelling, they currently only predict voltage outputs for the soma or a select few neuron compartments. Our novel approach, based on enhanced state-of-the-art architectures for multitask learning (MTL), allows for the simultaneous prediction of membrane potentials in each compartment of a neuron model, at a speed of up to two orders of magnitude faster than classical simulation methods. By predicting all membrane potentials together, our approach not only allows for comparison of model output with a wider range of experimental recordings (patch-electrode, voltage-sensitive dye imaging), it also provides the first stepping stone towards predicting local field potentials (LFPs), electroencephalogram (EEG) signals, and magnetoencephalography (MEG) signals from ANN-based simulations. While LFP and EEG are an important downstream application, the main focus of this paper lies in predicting dendritic voltages within each compartment to capture the entire electrophysiology of a biophysically-detailed neuron model. It further presents a challenging benchmark for MTL architectures due to the large amount of data involved, the presence of correlations between neighbouring compartments, and the non-Gaussian distribution of membrane potentials.

Single neuron contributions to the auditory brainstem EEG

The auditory brainstem response (ABR) is an acoustically evoked EEG potential that is an important diagnostic tool for hearing loss, especially in newborns. The ABR originates from the response sequence of auditory brainstem nuclei, and a click-evoked ABR typically shows three positive peaks ('waves') within the first six milliseconds. However, an assignment of the waves of the ABR to specific sources is difficult, and a quantification of contributions to the ABR waves is not available. Here, we exploit the large size and physical separation of the barn owl first-order cochlear nucleus magnocellularis (NM) to estimate single-cell contributions to the ABR. We simultaneously recorded NM neurons' spikes and the EEG, and found that ≳ 5, 000 spontaneous single-cell spikes are necessary to isolate a significant spike-triggered average response at the EEG electrode. An average single-neuron contribution to the ABR was predicted by convolving the spike-triggered average with the cell's peri-stimulus time histogram. Amplitudes of predicted contributions of single NM cells typically reached 32.9 ± 1.1 nV (mean ± SE, range: 2.5 - 162.7 nV), or 0.07 ± 0.02% (median ± SE range: 0.01 - 4.0%) of the ABR amplitude. The time of the predicted peak coincided best with the peak of the ABR wave II, and this coincidence was independent of the click sound level. Our results suggest that wave II of the ABR is shaped by a small fraction of NM units.

Uncovering population contributions to the extracellular potential in the mouse visual system using Laminar Population Analysis

The local field potential (LFP), the low-frequency part of the extracellular potential, reflects transmembrane currents in the vicinity of the recording electrode. Thought mainly to stem from currents caused by synaptic input, it provides information about neural activity complementary to that of spikes, the output of neurons. However, the many neural sources contributing to the LFP, and likewise the derived current source density (CSD), can often make it challenging to interpret. Efforts to improve its interpretability have included the application of statistical decomposition tools like principal component analysis (PCA) and independent component analysis (ICA) to disentangle the contributions from different neural sources. However, their underlying assumptions of, respectively, orthogonality and statistical independence are not always valid for the various processes or pathways generating LFP. Here, we expand upon and validate a decomposition algorithm named Laminar Population Analysis (LPA), which is based on physiological rather than statistical assumptions. LPA utilizes the multiunit activity (MUA) and LFP jointly to uncover the contributions of different populations to the LFP. To perform the validation of LPA, we used data simulated with the large-scale, biophysically detailed model of mouse V1 developed by the Allen Institute. We find that LPA can identify laminar positions within V1 and the temporal profiles of laminar population firing rates from the MUA. We also find that LPA can estimate the salient current sinks and sources generated by feedforward input from the lateral geniculate nucleus (LGN), recurrent activity in V1, and feedback input from the lateromedial (LM) area of visual cortex. LPA identifies and distinguishes these contributions with a greater accuracy than the alternative statistical decomposition methods, PCA and ICA. Lastly, we also demonstrate the application of LPA on experimentally recorded MUA and LFP from 24 animals in the publicly available Visual Coding dataset. Our results suggest that LPA can be used both as a method to estimate positions of laminar populations and to uncover salient features in LFP/CSD contributions from different populations.

Data-driven multiscale model of macaque auditory thalamocortical circuits reproduces in vivo dynamics

We developed a detailed model of macaque auditory thalamocortical circuits, including primary auditory cortex (A1), medial geniculate body (MGB), and thalamic reticular nucleus, utilizing the NEURON simulator and NetPyNE tool. The A1 model simulates a cortical column with over 12,000 neurons and 25 million synapses, incorporating data on cell-type-specific neuron densities, morphology, and connectivity across six cortical layers. It is reciprocally connected to the MGB thalamus, which includes interneurons and core and matrix-layer-specific projections to A1. The model simulates multiscale measures, including physiological firing rates, local field potentials (LFPs), current source densities (CSDs), and electroencephalography (EEG) signals. Laminar CSD patterns, during spontaneous activity and in response to broadband noise stimulus trains, mirror experimental findings. Physiological oscillations emerge spontaneously across frequency bands comparable to those recorded in vivo. We elucidate population-specific contributions to observed oscillation events and relate them to firing and presynaptic input patterns. The model offers a quantitative theoretical framework to integrate and interpret experimental data and predict its underlying cellular and circuit mechanisms.

Cortical origin of theta error signals

A multi-scale approach elucidated the origin of the error-related-negativity (ERN), with its associated theta-rhythm, and the post-error-positivity (Pe) in macaque supplementary eye field (SEF). Using biophysical modeling, synaptic inputs to a subpopulation of layer-3 (L3) and layer-5 (L5) pyramidal cells (PCs) were optimized to reproduce error-related spiking modulation and inter-spike intervals. The intrinsic dynamics of dendrites in L5 but not L3 error PCs generate theta rhythmicity with random phases. Saccades synchronized the phases of the theta-rhythm, which was magnified on errors. Contributions from error PCs to the laminar current source density (CSD) observed in SEF were negligible and could not explain the observed association between error-related spiking modulation in L3 PCs and scalp-EEG. CSD from recorded laminar field potentials in SEF was comprised of multipolar components, with monopoles indicating strong electro-diffusion, dendritic/axonal electrotonic current leakage outside SEF, or violations of the model assumptions. Our results also demonstrate the involvement of secondary cortical regions, in addition to SEF, particularly for the later Pe component. The dipolar component from the observed CSD paralleled the ERN dynamics, while the quadrupolar component paralleled the Pe. These results provide the most advanced explanation to date of the cellular mechanisms generating the ERN.

Diversity of cortical activity changes beyond depression during Spreading Depolarizations

Spreading depolarizations (SDs) are classically thought to be associated with spreading depression of cortical activity. Here, we found that SDs in patients with subarachnoid hemorrhage produce variable, ranging from depression to booming, changes in electrocorticographic activity, especially in the delta frequency band. In rats, depression of activity was characteristic of high-potassium-induced full SDs, whereas partial superficial SDs caused either little change or a boom of activity at the cortical vertex, supported by volume conduction of signals from spared delta generators in the deep cortical layers. Partial SDs also caused moderate neuronal depolarization and sustained excitation, organized in gamma oscillations in a narrow sub-SD zone. Thus, our study challenges the concept of homology between spreading depolarization and spreading depression by showing that SDs produce variable, from depression to booming, changes in activity at the cortical surface and in different cortical layers depending on the depth of SD penetration.

Uncovering circuit mechanisms of current sinks and sources with biophysical simulations of primary visual cortex

Local field potential (LFP) recordings reflect the dynamics of the current source density (CSD) in brain tissue. The synaptic, cellular, and circuit contributions to current sinks and sources are ill-understood. We investigated these in mouse primary visual cortex using public Neuropixels recordings and a detailed circuit model based on simulating the Hodgkin-Huxley dynamics of >50,000 neurons belonging to 17 cell types. The model simultaneously captured spiking and CSD responses and demonstrated a two-way dissociation: firing rates are altered with minor effects on the CSD pattern by adjusting synaptic weights, and CSD is altered with minor effects on firing rates by adjusting synaptic placement on the dendrites. We describe how thalamocortical inputs and recurrent connections sculpt specific sinks and sources early in the visual response, whereas cortical feedback crucially alters them in later stages. These results establish quantitative links between macroscopic brain measurements (LFP/CSD) and microscopic biophysics-based understanding of neuron dynamics and show that CSD analysis provides powerful constraints for modeling beyond those from considering spikes.