Synaptic Conductance Estimates of the Connection Between Local Inhibitor Interneurons and Pyramidal Neurons in Layer 2/3 of a Cortical Column

Recovery kinetics of short-term depression of GABAergic and glutamatergic synapses at layer 2/3 pyramidal cells in the mouse barrel cortex

Introduction

Short-term synaptic plasticity (STP) is a widespread mechanism underlying activity-dependent modifications of cortical networks.

Methods

To investigate how STP influences excitatory and inhibitory synapses in layer 2/3 of mouse barrel cortex, we combined whole-cell patch-clamp recordings from visually identified pyramidal neurons (PyrN) and parvalbumin-positive interneurons (PV-IN) of cortical layer 2/3 in acute slices with electrical stimulation of afferent fibers in layer 4 and optogenetic activation of PV-IN.

Results

These experiments revealed that electrical burst stimulation (10 pulses at 10 Hz) of layer 4 afferents to layer 2/3 neurons induced comparable short-term depression (STD) of glutamatergic postsynaptic currents (PSCs) in PyrN and in PV-IN, while disynaptic GABAergic PSCs in PyrN showed a stronger depression. Burst-induced depression of glutamatergic PSCs decayed within <4 s, while the decay of GABAergic PSCs required >11 s. Optogenetically-induced GABAergic PSCs in PyrN also demonstrated STD after burst stimulation, with a decay of >11 s. Excitatory postsynaptic potentials (EPSPs) in PyrN were unaffected after electrical burst stimulation, while a selective optogenetic STD of GABAergic synapses caused a transient increase of electrically evoked EPSPs in PyrN.

Discussion

In summary, these results demonstrate substantial short-term plasticity at all synapses investigated and suggest that the prominent STD observed in GABAergic synapses can moderate the functional efficacy of glutamatergic STD after repetitive synaptic stimulations. This mechanism may contribute to a reliable information flow toward the integrative layer 2/3 for complex time-varying sensory stimuli.

The role of superficial and deep layers in the generation of high frequency oscillations and interictal epileptiform discharges in the human cortex

Describing intracortical laminar organization of interictal epileptiform discharges (IED) and high frequency oscillations (HFOs), also known as ripples. Defining the frequency limits of slow and fast ripples. We recorded potential gradients with laminar multielectrode arrays (LME) for current source density (CSD) and multi-unit activity (MUA) analysis of interictal epileptiform discharges IEDs and HFOs in the neocortex and mesial temporal lobe of focal epilepsy patients. IEDs were observed in 20/29, while ripples only in 9/29 patients. Ripples were all detected within the seizure onset zone (SOZ). Compared to hippocampal HFOs, neocortical ripples proved to be longer, lower in frequency and amplitude, and presented non-uniform cycles. A subset of ripples (≈ 50%) co-occurred with IEDs, while IEDs were shown to contain variable high-frequency activity, even below HFO detection threshold. The limit between slow and fast ripples was defined at 150 Hz, while IEDs' high frequency components form clusters separated at 185 Hz. CSD analysis of IEDs and ripples revealed an alternating sink-source pair in the supragranular cortical layers, although fast ripple CSD appeared lower and engaged a wider cortical domain than slow ripples MUA analysis suggested a possible role of infragranularly located neural populations in ripple and IED generation. Laminar distribution of peak frequencies derived from HFOs and IEDs, respectively, showed that supragranular layers were dominated by slower (< 150 Hz) components. Our findings suggest that cortical slow ripples are generated primarily in upper layers while fast ripples and associated MUA in deeper layers. The dissociation of macro- and microdomains suggests that microelectrode recordings may be more selective for SOZ-linked ripples. We found a complex interplay between neural activity in the neocortical laminae during ripple and IED formation. We observed a potential leading role of cortical neurons in deeper layers, suggesting a refined utilization of LMEs in SOZ localization.

Unitary synaptic responses of parvalbumin interneurons evoked by excitatory neurons in the mouse barrel cortex

The mammalian cortex integrates and processes information to transform sensory inputs into perceptions and motor outputs. These operations are performed by networks of excitatory and inhibitory neurons distributed through the cortical layers. Parvalbumin interneurons (PVIs) are the most abundant type of inhibitory cortical neuron. With axons projecting within and between layers, PVIs supply feedforward and feedback inhibition to control and modulate circuit function. Distinct populations of excitatory neurons recruit different PVI populations, but the specializations of these synapses are poorly understood. Here, we targeted a genetically encoded hybrid voltage sensor to PVIs and used fluorescence imaging in mouse somatosensory cortex slices to record their voltage changes. Stimulating a single visually identified excitatory neuron with small-tipped theta-glass electrodes depolarized multiple PVIs, and a common threshold suggested that stimulation elicited unitary synaptic potentials in response to a single excitatory neuron. Excitatory neurons depolarized PVIs in multiple layers, with the most residing in the layer of the stimulated neuron. Spiny stellate cells depolarized PVIs more strongly than pyramidal cells by up to 77%, suggesting a greater role for stellate cells in recruiting PVI inhibition and controlling cortical computations. Response half-width also varied between different excitatory inputs. These results demonstrate functional differences between excitatory synapses on PVIs.

Disinhibitory circuit mediated by connections from vasoactive intestinal polypeptide to somatostatin interneurons underlies the paradoxical decrease in spike synchrony with increased border ownership selective neuron firing rate

The activity of border ownership selective (BOS) neurons in intermediate-level visual areas indicates which side of a contour owns a border relative to its classical receptive field and provides a fundamental component of figure-ground segregation. A physiological study reported that selective attention facilitates the activity of BOS neurons with a consistent border ownership preference, defined as two neurons tuned to respond to the same visual object. However, spike synchrony between this pair is significantly suppressed by selective attention. These neurophysiological findings are derived from a biologically-plausible microcircuit model consisting of spiking neurons including two subtypes of inhibitory interneurons, somatostatin (SOM) and vasoactive intestinal polypeptide (VIP) interneurons, and excitatory BOS model neurons. In our proposed model, BOS neurons and SOM interneurons cooperate and interact with each other. VIP interneurons not only suppress SOM interneuron responses but also are activated by feedback signals mediating selective attention, which leads to disinhibition of BOS neurons when they are directing selective attention toward an object. Our results suggest that disinhibition arising from the synaptic connections from VIP to SOM interneurons plays a critical role in attentional modulation of neurons in intermediate-level visual areas.

A microcircuit model involving parvalbumin, somatostatin, and vasoactive intestinal polypeptide inhibitory interneurons for the modulation of neuronal oscillation during visual processing

Various subtypes of inhibitory interneurons contact one another to organize cortical networks. Most cortical inhibitory interneurons express 1 of 3 genes: parvalbumin (PV), somatostatin (SOM), or vasoactive intestinal polypeptide (VIP). This diversity of inhibition allows the flexible regulation of neuronal responses within and between cortical areas. However, the exact roles of these interneuron subtypes and of excitatory pyramidal (Pyr) neurons in regulating neuronal network activity and establishing perception (via interactions between feedforward sensory and feedback attentional signals) remain largely unknown. To explore the regulatory roles of distinct neuronal types in cortical computation, we developed a computational microcircuit model with biologically plausible visual cortex layers 2/3 that combined Pyr neurons and the 3 inhibitory interneuron subtypes to generate network activity. In simulations with our model, inhibitory signals from PV and SOM neurons preferentially induced neuronal firing at gamma (30-80 Hz) and beta (20-30 Hz) frequencies, respectively, in agreement with observed physiological results. Furthermore, our model indicated that rapid inhibition from VIP to SOM subtypes underlies marked attentional modulation for low-gamma frequency (30-50 Hz) in Pyr neuron responses. Our results suggest the distinct but cooperative roles of inhibitory interneuron subtypes in the establishment of visual perception.

Unique properties of dually innervated dendritic spines in pyramidal neurons of the somatosensory cortex uncovered by 3D correlative light and electron microscopy

Pyramidal neurons (PNs) are covered by thousands of dendritic spines receiving excitatory synaptic inputs. The ultrastructure of dendritic spines shapes signal compartmentalization, but ultrastructural diversity is rarely taken into account in computational models of synaptic integration. Here, we developed a 3D correlative light-electron microscopy (3D-CLEM) approach allowing the analysis of specific populations of synapses in genetically defined neuronal types in intact brain circuits. We used it to reconstruct segments of basal dendrites of layer 2/3 PNs of adult mouse somatosensory cortex and quantify spine ultrastructural diversity. We found that 10% of spines were dually innervated and 38% of inhibitory synapses localized to spines. Using our morphometric data to constrain a model of synaptic signal compartmentalization, we assessed the impact of spinous versus dendritic shaft inhibition. Our results indicate that spinous inhibition is locally more efficient than shaft inhibition and that it can decouple voltage and calcium signaling, potentially impacting synaptic plasticity.

Persistent sodium current blockers can suppress seizures caused by loss of low-threshold D-type potassium currents: Predictions from an in silico study of Kv1 channel disorders

Objective

Ion channels belonging to subfamily A of voltage‐gated potassium channels (K_v1) are highly expressed on axons, where they play a key role in determining resting membrane potential, in shaping action potentials, and in modulating action potential frequency during repetitive neuronal firing. We aimed to study the genesis of seizures caused by mutations affecting K_v1 channels and searched for potential therapeutic targets.

Methods

We used a novel in silico model, the laminar cortex model (LCM), to examine changes in neuronal excitability and network dynamics associated with loss‐of‐function mutations in K_v1 channels. The LCM simulates the activities of a network of tens of thousands of interconnected neurons and incorporates the kinetics of 11 types of ion channel and three classes of neurotransmitter receptor. Changes in two types of potassium currents conducted by K_v1 channels were examined: slowly inactivating D‐type currents and rapidly inactivating A‐type currents. Effects on neuronal firing rate, action potential shape, and neuronal oscillation state were evaluated. A systematic parameter scan was performed to identify parameter changes that can reverse the effects of the changes.

Results

Reduced axonal D‐type currents led to lower firing threshold and widened action potentials, both lowering the seizure threshold. Two potential therapeutic targets for treating seizures caused by loss‐of‐function changes in K_v1 channels were identified: persistent sodium channels and NMDA receptors. Blocking persistent sodium channels restored the firing threshold and reduced action potential width. NMDA receptor antagonists reduced excitatory postsynaptic currents from excessive glutamate release related to widened action potentials.

Significance

Riluzole reduces persistent sodium currents and excitatory postsynaptic currents from NMDA receptor activation. Our results suggest that this FDA‐approved drug can be repurposed to treat epilepsies caused by mutations affecting axonal K_v1 channels.

Combined Positive and Negative Feedback Allows Modulation of Neuronal Oscillation Frequency during Sensory Processing

A key step in sensory information processing involves modulation and integration of neuronal oscillations in disparate frequency bands, a poorly understood process. Here, we investigate how top-down input causes frequency changes in slow oscillations during sensory processing and, in turn, how the slow oscillations are combined with fast oscillations (which encode sensory input). Using experimental connectivity patterns and strengths of interneurons, we develop a system-level model of a neuronal circuit controlling these oscillatory behaviors, allowing us to understand the mechanisms responsible for the observed oscillatory behaviors. Our analysis discovers a circuit capable of producing the observed oscillatory behaviors and finds that a detailed balance in the strength of synaptic connections is the critical determinant to produce such oscillatory behaviors. We not only uncover how disparate frequency bands are modulated and combined but also give insights into the causes of abnormal neuronal activities present in brain disorders.

Anatomical Correlates of Local, Translaminar, and Transcolumnar Inhibition by Layer 6 GABAergic Interneurons in Somatosensory Cortex

In the vibrissal area of rodent somatosensory cortex, information on whisker stimulation is processed by neuronal networks in a corresponding cortical column. To understand how sensory stimuli are represented in a column, it is essential to identify cell types constituting these networks. Layer 6 (L6) comprises 25% of all neurons in a column. In rats, 430 of these are inhibitory interneurons (INs). Little is known about the axon projection of L6 INs with reference to columnar and laminar organization. We quantified axonal projections of L6 INs (n = 68) with reference to columns and layers in somatosensory cortex of rats. We found distinct projection types differentially targeting layers of a cortical column. The majority of L6 INs did not show a column-specific innervation, densely projecting to neighboring columns as well as the home column. However, a small fraction targeted granular and supragranular layers, where axon projections were confined to the home column. We also quantified putative innervation of pyramidal cells as a functional correlate of axonal distribution. Electrophysiological properties were not correlated to axon projection. The quantitative data on axonal projections and electrophysiological properties of L6 INs can guide future studies investigating cortical processing of sensory information at the single cell level.

Reconstructing neuronal circuitry from parallel spike trains

State-of-the-art techniques allow researchers to record large numbers of spike trains in parallel for many hours. With enough such data, we should be able to infer the connectivity among neurons. Here we develop a method for reconstructing neuronal circuitry by applying a generalized linear model (GLM) to spike cross-correlations. Our method estimates connections between neurons in units of postsynaptic potentials and the amount of spike recordings needed to verify connections. The performance of inference is optimized by counting the estimation errors using synthetic data. This method is superior to other established methods in correctly estimating connectivity. By applying our method to rat hippocampal data, we show that the types of estimated connections match the results inferred from other physiological cues. Thus our method provides the means to build a circuit diagram from recorded spike trains, thereby providing a basis for elucidating the differences in information processing in different brain regions.

Current techniques have enabled the simultaneous collection of spike train data from large numbers of neurons. Here, the authors report a method to infer the underlying neural circuit connectivity diagram based on a generalized linear model applied to spike cross-correlations between neurons.

Quantifying How Staining Methods Bias Measurements of Neuron Morphologies

The process through which neurons are labeled is a key methodological choice in measuring neuron morphology. However, little is known about how this choice may bias measurements. To quantify this bias we compare the extracted morphology of neurons collected from the same rodent species, experimental condition, gender distribution, age distribution, brain region and putative cell type, but obtained with 19 distinct staining methods. We found strong biases on measured features of morphology. These were largest in features related to the coverage of the dendritic tree (e.g., the total dendritic tree length). Understanding measurement biases is crucial for interpreting morphological data.

Leveraging heterogeneity for neural computation with fading memory in layer 2/3 cortical microcircuits

Complexity and heterogeneity are intrinsic to neurobiological systems, manifest in every process, at every scale, and are inextricably linked to the systems’ emergent collective behaviours and function. However, the majority of studies addressing the dynamics and computational properties of biologically inspired cortical microcircuits tend to assume (often for the sake of analytical tractability) a great degree of homogeneity in both neuronal and synaptic/connectivity parameters. While simplification and reductionism are necessary to understand the brain’s functional principles, disregarding the existence of the multiple heterogeneities in the cortical composition, which may be at the core of its computational proficiency, will inevitably fail to account for important phenomena and limit the scope and generalizability of cortical models. We address these issues by studying the individual and composite functional roles of heterogeneities in neuronal, synaptic and structural properties in a biophysically plausible layer 2/3 microcircuit model, built and constrained by multiple sources of empirical data. This approach was made possible by the emergence of large-scale, well curated databases, as well as the substantial improvements in experimental methodologies achieved over the last few years. Our results show that variability in single neuron parameters is the dominant source of functional specialization, leading to highly proficient microcircuits with much higher computational power than their homogeneous counterparts. We further show that fully heterogeneous circuits, which are closest to the biophysical reality, owe their response properties to the differential contribution of different sources of heterogeneity.

Cortical microcircuits are highly inhomogeneous dynamical systems whose information processing capacity is determined by the characteristics of its heterogeneous components and their complex interactions. The high degree of variability that characterizes macroscopic population dynamics, both during ongoing, spontaneous activity and active processing states reflects the underlying complexity and heterogeneity which has the potential to dramatically constrain the space of functions that any given circuit can compute, leading to richer and more expressive information processing systems. In this study, we identify different tentative sources of heterogeneity and assess their differential and cumulative contribution to the microcircuit’s dynamics and information processing capacity. We study these properties in a generic Layer 2/3 cortical microcircuit model, built and constrained by multiple sources of experimental data, and demonstrate that heterogeneity in neuronal properties and microconnectivity structure are important sources of functional specialization, greatly improving the circuit’s processing capacity, while capturing various important features of cortical physiology.

Brief activation of GABAergic interneurons initiates the transition to ictal events through post-inhibitory rebound excitation

Activation of γ-aminobutyric acid (GABA_A) receptors have been associated with the onset of epileptiform events. To investigate if a causal relationship exists between GABA_A receptor activation and ictal event onset, we activated inhibitory GABAergic networks in the superficial layer (2/3) of the somatosensory cortex during hyperexcitable conditions using optogenetic techniques in mice expressing channelrhodopsin-2 in all GABAergic interneurons. We found that a brief 30ms light pulse reliably triggered either an interictal-like event (IIE) or ictal-like ("ictal") event in the in vitro cortical 4-Aminopyridine (4-AP) slice model. The link between light pulse and epileptiform event onset was lost following blockade of GABA_A receptors with bicuculline methiodide. Additionally, recording the chronological sequence of events following a light pulse in a variety of configurations (whole-cell, gramicidin-perforated patch, and multi-electrode array) demonstrated an initial hyperpolarization followed by post-inhibitory rebound spiking and a subsequent slow depolarization at the transition to epileptiform activity. Furthermore, the light-triggered ictal events were independent of the duration or intensity of the initiating light pulse, suggesting an underlying regenerative mechanism. Moreover, we demonstrated that brief GABA_A receptor activation can initiate ictal events in the in vivo 4-AP mouse model, in another common in vitro model of epileptiform activity, and in neocortical tissue resected from epilepsy patients. Our findings reveal that the synchronous activation of GABAergic interneurons is a robust trigger for ictal event onset in hyperexcitable cortical networks.

SynEM, automated synapse detection for connectomics

Nerve tissue contains a high density of chemical synapses, about 1 per µm3 in the mammalian cerebral cortex. Thus, even for small blocks of nerve tissue, dense connectomic mapping requires the identification of millions to billions of synapses. While the focus of connectomic data analysis has been on neurite reconstruction, synapse detection becomes limiting when datasets grow in size and dense mapping is required. Here, we report SynEM, a method for automated detection of synapses from conventionally en-bloc stained 3D electron microscopy image stacks. The approach is based on a segmentation of the image data and focuses on classifying borders between neuronal processes as synaptic or non-synaptic. SynEM yields 97% precision and recall in binary cortical connectomes with no user interaction. It scales to large volumes of cortical neuropil, plausibly even whole-brain datasets. SynEM removes the burden of manual synapse annotation for large densely mapped connectomes.

Robustness of sensory-evoked excitation is increased by inhibitory inputs to distal apical tuft dendrites

Cortical inhibitory interneurons (INs) are subdivided into a variety of morphologically and functionally specialized cell types. How the respective specific properties translate into mechanisms that regulate sensory-evoked responses of pyramidal neurons (PNs) remains unknown. Here, we investigated how INs located in cortical layer 1 (L1) of rat barrel cortex affect whisker-evoked responses of L2 PNs. To do so we combined in vivo electrophysiology and morphological reconstructions with computational modeling. We show that whisker-evoked membrane depolarization in L2 PNs arises from highly specialized spatiotemporal synaptic input patterns. Temporally L1 INs and L2-5 PNs provide near synchronous synaptic input. Spatially synaptic contacts from L1 INs target distal apical tuft dendrites, whereas PNs primarily innervate basal and proximal apical dendrites. Simulations of such constrained synaptic input patterns predicted that inactivation of L1 INs increases trial-to-trial variability of whisker-evoked responses in L2 PNs. The in silico predictions were confirmed in vivo by L1-specific pharmacological manipulations. We present a mechanism-consistent with the theory of distal dendritic shunting-that can regulate the robustness of sensory-evoked responses in PNs without affecting response amplitude or latency.