Computational models of O-LM cells are recruited by low or high theta frequency inputs depending on h-channel distributions

Synaptic and dendritic architecture of different types of hippocampal somatostatin interneurons

GABAergic inhibitory neurons fundamentally shape the activity and plasticity of cortical circuits. A major subset of these neurons contains somatostatin (SOM); these cells play crucial roles in neuroplasticity, learning, and memory in many brain areas including the hippocampus, and are implicated in several neuropsychiatric diseases and neurodegenerative disorders. Two main types of SOM-containing cells in area CA1 of the hippocampus are oriens-lacunosum-moleculare (OLM) cells and hippocampo-septal (HS) cells. These cell types show many similarities in their soma-dendritic architecture, but they have different axonal targets, display different activity patterns in vivo, and are thought to have distinct network functions. However, a complete understanding of the functional roles of these interneurons requires a precise description of their intrinsic computational properties and their synaptic interactions. In the current study we generated, analyzed, and make available several key data sets that enable a quantitative comparison of various anatomical and physiological properties of OLM and HS cells in mouse. The data set includes detailed scanning electron microscopy (SEM)-based 3D reconstructions of OLM and HS cells along with their excitatory and inhibitory synaptic inputs. Combining this core data set with other anatomical data, patch-clamp electrophysiology, and compartmental modeling, we examined the precise morphological structure, inputs, outputs, and basic physiological properties of these cells. Our results highlight key differences between OLM and HS cells, particularly regarding the density and distribution of their synaptic inputs and mitochondria. For example, we estimated that an OLM cell receives about 8,400, whereas an HS cell about 15,600 synaptic inputs, about 16% of which are GABAergic. Our data and models provide insight into the possible basis of the different functionality of OLM and HS cell types and supply essential information for more detailed functional models of these neurons and the hippocampal network.

The alpha2 nicotinic acetylcholine receptor, a subunit with unique and selective expression in inhibitory interneurons associated with principal cells

Nicotinic acetylcholine receptors (nAChRs) play crucial roles in various human disorders, with the α7, α4, α6, and α3-containing nAChR subtypes extensively studied in relation to conditions such as Alzheimer's disease, Parkinson's disease, nicotine dependence, mood disorders, and stress disorders. In contrast, the α2-nAChR subunit has received less attention due to its more restricted expression and the scarcity of specific agonists and antagonists for studying its function. Nevertheless, recent research has shed light on the unique expression pattern of the Chrna2 gene, which encodes the α2-nAChR subunit, and its involvement in distinct populations of inhibitory interneurons. This review highlights the structure, pharmacology, localization, function, and disease associations of α2-containing nAChRs and points to the unique expression pattern of the Chrna2 gene and its role in different inhibitory interneuron populations. These populations, including the oriens lacunosum moleculare (OLM) cells in the hippocampus, Martinotti cells in the neocortex, and Renshaw cells in the spinal cord, share common features and contribute to recurrent inhibitory microcircuits. Thus, the α2-nAChR subunit's unique expression pattern in specific interneuron populations and its role in recurrent inhibitory microcircuits highlight its importance in various physiological processes. Further research is necessary to uncover the comprehensive functionality of α2-containing nAChRs, delineate their specific contributions to neuronal circuits, and investigate their potential as therapeutic targets for related disorders.

Theta-gamma phase amplitude coupling in a hippocampal CA1 microcircuit

Phase amplitude coupling (PAC) between slow and fast oscillations is found throughout the brain and plays important functional roles. Its neural origin remains unclear. Experimental findings are often puzzling and sometimes contradictory. Most computational models rely on pairs of pacemaker neurons or neural populations tuned at different frequencies to produce PAC. Here, using a data-driven model of a hippocampal microcircuit, we demonstrate that PAC can naturally emerge from a single feedback mechanism involving an inhibitory and excitatory neuron population, which interplay to generate theta frequency periodic bursts of higher frequency gamma. The model suggests the conditions under which a CA1 microcircuit can operate to elicit theta-gamma PAC, and highlights the modulatory role of OLM and PVBC cells, recurrent connectivity, and short term synaptic plasticity. Surprisingly, the results suggest the experimentally testable prediction that the generation of the slow population oscillation requires the fast one and cannot occur without it.

Reduced oriens-lacunosum/moleculare cell model identifies biophysical current balances for in vivo theta frequency spiking resonance

Conductance-based models have played an important role in the development of modern neuroscience. These mathematical models are powerful "tools" that enable theoretical explorations in experimentally untenable situations, and can lead to the development of novel hypotheses and predictions. With advances in cell imaging and computational power, multi-compartment models with morphological accuracy are becoming common practice. However, as more biological details are added, they make extensive explorations and analyses more challenging largely due to their huge computational expense. Here, we focus on oriens-lacunosum/moleculare (OLM) cell models. OLM cells can contribute to functionally relevant theta rhythms in the hippocampus by virtue of their ability to express spiking resonance at theta frequencies, but what characteristics underlie this is far from clear. We converted a previously developed detailed multi-compartment OLM cell model into a reduced single compartment model that retained biophysical fidelity with its underlying ion currents. We showed that the reduced OLM cell model can capture complex output that includes spiking resonance in in vivo-like scenarios as previously obtained with the multi-compartment model. Using the reduced model, we were able to greatly expand our in vivo-like scenarios. Applying spike-triggered average analyses, we were able to to determine that it is a combination of hyperpolarization-activated cation and muscarinic type potassium currents that specifically allow OLM cells to exhibit spiking resonance at theta frequencies. Further, we developed a robust Kalman Filtering (KF) method to estimate parameters of the reduced model in real-time. We showed that it may be possible to directly estimate conductance parameters from experiments since this KF method can reliably extract parameter values from model voltage recordings. Overall, our work showcases how the contribution of cellular biophysical current details could be determined and assessed for spiking resonance. As well, our work shows that it may be possible to directly extract these parameters from current clamp voltage recordings.

Deciphering how interneuron specific 3 cells control oriens lacunosum-moleculare cells to contribute to circuit function

The wide diversity of inhibitory cells across the brain makes them suitable to contribute to network dynamics in specialized fashions. However, the contributions of a particular inhibitory cell type in a behaving animal are challenging to untangle as one needs to both record cellular activities and identify the cell type being recorded. Thus, using computational modeling and theory to predict and hypothesize cell-specific contributions is desirable. Here, we examine potential contributions of interneuron-specific 3 (I-S3) cells - an inhibitory interneuron found in CA1 hippocampus that only targets other inhibitory interneurons - during simulated theta rhythms. We use previously developed multi-compartment models of oriens lacunosum-moleculare (OLM) cells, the main target of I-S3 cells, and explore how I-S3 cell inputs during in vitro and in vivo scenarios contribute to theta. We find that I-S3 cells suppress OLM cell spiking, rather than engender its spiking via post-inhibitory rebound mechanisms, and contribute to theta frequency spike resonance during simulated in vivo scenarios. To elicit recruitment similar to in vitro experiments, inclusion of disinhibited pyramidal cell inputs is necessary, implying that I-S3 cell firing broadens the window for pyramidal cell disinhibition. Using in vivo virtual networks, we show that I-S3 cells contribute to a sharpening of OLM cell recruitment at theta frequencies. Further, shifting the timing of I-S3 cell spiking due to external modulation shifts the timing of the OLM cell firing and thus disinhibitory windows. We propose a specialized contribution of I-S3 cells to create temporally precise coordination of modulation pathways.

Effects of Several Classes of Voltage-Gated Ion Channel Conductances on Gamma and Theta Oscillations in a Hippocampal Microcircuit Model

Gamma and theta oscillations have been functionally associated with cognitive processes, such as learning and memory. Synaptic conductances play an important role in the generation of intrinsic network rhythmicity, but few studies have examined the effects of voltage-gated ion channels (VGICs) on these rhythms. In this report, we have used a pyramidal-interneuron-gamma (PING) network consisting of excitatory pyramidal cells and two types of inhibitory interneurons. We have constructed a conductance-based neural network incorporating a persistent sodium current (I_NaP), a delayed rectifier potassium current (I_KDR), a inactivating potassium current (I_A) and a hyperpolarization-activated current (I_H). We have investigated the effects of several conductances on network theta and gamma frequency oscillations. Variation of all conductances of interest changed network rhythmicity. Theta power was altered by all conductances tested. Gamma rhythmogenesis was dependent on I_A and I_H. The I_KDR currents in excitatory pyramidal cells as well as both types of inhibitory interneurons were essential for theta rhythmogenesis and altered gamma rhythm properties. Increasing I_NaP suppressed both gamma and theta rhythms. Addition of noise did not alter these patterns. Our findings suggest that VGICs strongly affect brain network rhythms. Further investigations in vivo will be of great interest, including potential effects on neural function and cognition.

Septohippocampal transmission from parvalbumin-positive neurons features rapid recovery from synaptic depression

Parvalbumin-containing projection neurons of the medial-septum-diagonal band of Broca ([Formula: see text]) are essential for hippocampal rhythms and learning operations yet are poorly understood at cellular and synaptic levels. We combined electrophysiological, optogenetic, and modeling approaches to investigate [Formula: see text] neuronal properties. [Formula: see text] neurons had intrinsic membrane properties distinct from acetylcholine- and somatostatin-containing MS-DBB subtypes. Viral expression of the fast-kinetic channelrhodopsin ChETA-YFP elicited action potentials to brief (1-2 ms) 470 nm light pulses. To investigate [Formula: see text] transmission, light pulses at 5-50 Hz frequencies generated trains of inhibitory postsynaptic currents (IPSCs) in CA1 stratum oriens interneurons. Using a similar approach, optogenetic activation of local hippocampal PV ([Formula: see text]) neurons generated trains of [Formula: see text]-mediated IPSCs in CA1 pyramidal neurons. Both synapse types exhibited short-term depression (STD) of IPSCs. However, relative to [Formula: see text] synapses, [Formula: see text] synapses possessed lower initial release probability, transiently resisted STD at gamma (20-50 Hz) frequencies, and recovered more rapidly from synaptic depression. Experimentally-constrained mathematical synapse models explored mechanistic differences. Relative to the [Formula: see text] model, the [Formula: see text] model exhibited higher sensitivity to calcium accumulation, permitting a faster rate of calcium-dependent recovery from STD. In conclusion, resistance of [Formula: see text] synapses to STD during short gamma bursts enables robust long-range GABAergic transmission from MS-DBB to hippocampus.

Modeling Reveals Human-Rodent Differences in H-Current Kinetics Influencing Resonance in Cortical Layer 5 Neurons

While our understanding of human neurons is often inferred from rodent data, inter-species differences between neurons can be captured by building cellular models specifically from human data. This includes understanding differences at the level of ion channels and their implications for human brain function. Thus, we here present a full spiking, biophysically detailed multi-compartment model of a human layer 5 (L5) cortical pyramidal cell. Model development was primarily based on morphological and electrophysiological data from the same human L5 neuron, avoiding confounds of experimental variability. Focus was placed on describing the behavior of the hyperpolarization-activated cation (h-) channel, given increasing interest in this channel due to its role in pacemaking and differentiating cell types. We ensured that the model exhibited post-inhibitory rebound spiking considering its relationship with the h-current, along with other general spiking characteristics. The model was validated against data not used in its development, which highlighted distinctly slower kinetics of the human h-current relative to the rodent setting. We linked the lack of subthreshold resonance observed in human L5 neurons to these human-specific h-current kinetics. This work shows that it is possible and necessary to build human-specific biophysical neuron models in order to understand human brain dynamics.

Integration of Within-Cell Experimental Data With Multi-Compartmental Modeling Predicts H-Channel Densities and Distributions in Hippocampal OLM Cells

Determining biophysical details of spatially extended neurons is a challenge that needs to be overcome if we are to understand the dynamics of brain function from cellular perspectives. Moreover, we now know that we should not average across recordings from many cells of a given cell type to obtain quantitative measures such as conductance since measures can vary multiple-fold for a given cell type. In this work we examine whether a tight combination of experimental and computational work can address this challenge. The oriens-lacunosum/moleculare (OLM) interneuron operates as a “gate” that controls incoming sensory and ongoing contextual information in the CA1 of the hippocampus, making it essential to understand how its biophysical properties contribute to memory function. OLM cells fire phase-locked to the prominent hippocampal theta rhythms, and we previously used computational models to show that OLM cells exhibit high or low theta spiking resonance frequencies that depend respectively on whether their dendrites have hyperpolarization-activated cation channels (h-channels) or not. However, whether OLM cells actually possess dendritic h-channels is unknown at present. We performed a set of whole-cell recordings of OLM cells from mouse hippocampus and constructed three multi-compartment models using morphological and electrophysiological parameters extracted from the same OLM cell, including per-cell pharmacologically isolated h-channel currents. We found that the models best matched experiments when h-channels were present in the dendrites of each of the three model cells created. This strongly suggests that h-channels must be present in OLM cell dendrites and are not localized to their somata. Importantly, this work shows that a tight integration of model and experiment can help tackle the challenge of characterizing biophysical details and distributions in spatially extended neurons. Full spiking models were built for two of the OLM cells, matching their current clamp cell-specific electrophysiological recordings. Overall, our work presents a technical advancement in modeling OLM cells. Our models are available to the community to use to gain insight into cellular dynamics underlying hippocampal function.

Computationally going where experiments cannot: a dynamical assessment of dendritic ion channel currents during in vivo-like states

Background: Despite technological advances, how specific cell types are involved in brain function remains shrouded in mystery. Further, little is known about the contribution of different ion channel currents to cell excitability across different neuronal subtypes and their dendritic compartments in vivo. The picture that we do have is largely based on somatic recordings performed in vitro. Uncovering dendritic ion channel current contributions in neuron subtypes that represent a minority of the neuronal population is not currently a feasible task using purely experimental means.

Methods: We employ two morphologically-detailed multi-compartment models of a specific type of inhibitory interneuron, the oriens lacunosum moleculare (OLM) cell. The OLM cell is a well-studied cell type in CA1 hippocampus that is important in gating sensory and contextual information. We create in vivo-like states for these cellular models by including levels of synaptic bombardment that would occur in vivo. Using visualization tools and analyses we assess the ion channel current contribution profile across the different somatic and dendritic compartments of the models.

Results: We identify changes in dendritic excitability, ion channel current contributions and co-activation patterns between in vitro and in vivo-like states. Primarily, we find that the relative timing between ion channel currents are mostly invariant between states, but exhibit changes in magnitudes and decreased propagation across dendritic compartments. We also find enhanced dendritic hyperpolarization-activated cyclic nucleotide-gated channel (h-channel) activation during in vivo-like states, which suggests that dendritically located h-channels are functionally important in altering signal propagation in the behaving animal.

Conclusions: Overall, we have demonstrated, using computational modelling, the dynamical changes that can occur to ion channel mechanisms governing neuronal spiking. Simultaneous access to dendritic compartments during simulated in vivo states shows that the magnitudes of some ion channel current contributions are differentially altered during in vivo-like states relative to in vitro.

Presynaptic GABAB receptors functionally uncouple somatostatin interneurons from the active hippocampal network

Information processing in cortical neuronal networks relies on properly balanced excitatory and inhibitory neurotransmission. A ubiquitous motif for maintaining this balance is the somatostatin interneuron (SOM-IN) feedback microcircuit. Here, we investigated the modulation of this microcircuit by presynaptic GABA_B receptors (GABA_BRs) in the rodent hippocampus. Whole-cell recordings from SOM-INs revealed that both excitatory and inhibitory synaptic inputs are strongly inhibited by GABA_BRs, while optogenetic activation of the interneurons shows that their inhibitory output is also strongly suppressed. Electron microscopic analysis of immunogold-labelled freeze-fracture replicas confirms that GABA_BRs are highly expressed presynaptically at both input and output synapses of SOM-INs. Activation of GABA_BRs selectively suppresses the recruitment of SOM-INs during gamma oscillations induced in vitro. Thus, axonal GABA_BRs are positioned to efficiently control the input and output synapses of SOM-INs and can functionally uncouple them from local network with implications for rhythmogenesis and the balance of entorhinal versus intrahippocampal afferents.

Using computational models to predict in vivo synaptic inputs to interneuron specific 3 (IS3) cells of CA1 hippocampus that also allow their recruitment during rhythmic states

Brain coding strategies are enabled by the balance of synaptic inputs that individual neurons receive as determined by the networks in which they reside. Inhibitory cell types contribute to brain function in distinct ways but recording from specific, inhibitory cell types during behaviour to determine their contributions is highly challenging. In particular, the in vivo activities of vasoactive intestinal peptide-expressing interneuron specific 3 (IS3) cells in the hippocampus that only target other inhibitory cells are unknown at present. We perform a massive, computational exploration of possible synaptic inputs to IS3 cells using multi-compartment models and optimized synaptic parameters. We find that asynchronous, in vivo-like states that are sensitive to additional theta-timed inputs (8 Hz) exist when excitatory and inhibitory synaptic conductances are approximately equally balanced and with low numbers of activated synapses receiving correlated inputs. Specifically, under these balanced conditions, the input resistance is larger with higher mean spike firing rates relative to other activated synaptic conditions investigated. Incoming theta-timed inputs result in strongly increased spectral power relative to baseline. Thus, using a generally applicable computational approach we predict the existence and features of background, balanced states in hippocampal circuits.

Input-Specific Synaptic Location and Function of the α5 GABAA Receptor Subunit in the Mouse CA1 Hippocampal Neurons

Hippocampus-dependent learning processes are coordinated via a large diversity of GABAergic inhibitory mechanisms. The α5 subunit-containing GABA_A receptor (α5-GABA_AR) is abundantly expressed in the hippocampus populating primarily the extrasynaptic domain of CA1 pyramidal cells, where it mediates tonic inhibitory conductance and may cause functional deficits in synaptic plasticity and hippocampus-dependent memory. However, little is known about synaptic expression of the α5-GABA_AR and, accordingly, its location site-specific function. We examined the cell- and synapse-specific distribution of the α5-GABA_AR in the CA1 stratum oriens/alveus (O/A) using a combination of immunohistochemistry, whole-cell patch-clamp recordings and optogenetic stimulation in hippocampal slices obtained from mice of either sex. In addition, the input-specific role of the α5-GABA_AR in spatial learning and anxiety-related behavior was studied using behavioral testing and chemogenetic manipulations. We demonstrate that α5-GABA_AR is preferentially targeted to the inhibitory synapses made by the vasoactive intestinal peptide (VIP)- and calretinin-positive terminals onto dendrites of somatostatin-expressing interneurons. In contrast, synapses made by the parvalbumin-positive inhibitory inputs to O/A interneurons showed no or little α5-GABA_AR. Inhibiting the α5-GABA_AR in control mice in vivo improved spatial learning but also induced anxiety-like behavior. Inhibiting the α5-GABA_AR in mice with inactivated CA1 VIP input could still improve spatial learning and was not associated with anxiety. Together, these data indicate that the α5-GABA_AR-mediated phasic inhibition via VIP input to interneurons plays a predominant role in the regulation of anxiety while the α5-GABA_AR tonic inhibition via this subunit may control spatial learning.SIGNIFICANCE STATEMENT The α5-GABA_AR subunit exhibits high expression in the hippocampus, and regulates the induction of synaptic plasticity and the hippocampus-dependent mnemonic processes. In CA1 principal cells, this subunit occupies mostly extrasynaptic sites and mediates tonic inhibition. Here, we provide evidence that, in CA1 somatostatin-expressing interneurons, the α5-GABA_AR subunit is targeted to synapses formed by the VIP- and calretinin-expressing inputs, and plays a specific role in the regulation of anxiety-like behavior.

Ventral hippocampal OLM cells control type 2 theta oscillations and response to predator odor

Dorsal and ventral hippocampus regions exert cognition and emotion-related functions, respectively. Since both regions display rhythmic activity, specific neural oscillatory pacemakers may underlie their functional dichotomy. Type 1 theta oscillations are independent of cholinergic transmission and are observed in the dorsal hippocampus during movement and exploration. In contrast, type 2 theta depends on acetylcholine and appears when animals are exposed to emotionally laden contexts such as a predator presence. Despite its involvement in emotions, type 2 theta has not been associated with the ventral hippocampus. Here, we show that optogenetic activation of oriens-lacunosum moleculare (OLM) interneurons in the ventral hippocampus drives type 2 theta. Moreover, we found that type 2 theta generation is associated with increased risk-taking behavior in response to predator odor. These results demonstrate that two theta oscillations subtypes originate in the two hippocampal regions that predominantly underlie either cognitive or emotion-related functions.

Combining Theory, Model, and Experiment to Explain How Intrinsic Theta Rhythms Are Generated in an In Vitro Whole Hippocampus Preparation without Oscillatory Inputs

Scientists have observed local field potential theta rhythms (3-12 Hz) in the hippocampus for decades, but understanding the mechanisms underlying their generation is complicated by their diversity in pharmacological and frequency profiles. In addition, interactions with other brain structures and oscillatory drives to the hippocampus during distinct brain states has made it difficult to identify hippocampus-specific properties directly involved in theta generation. To overcome this, we develop cellular-based network models using a whole hippocampus in vitro preparation that spontaneously generates theta rhythms. Building on theoretical and computational analyses, we find that spike frequency adaptation and postinhibitory rebound constitute a basis for theta generation in large, minimally connected CA1 pyramidal (PYR) cell network models with fast-firing parvalbumin-positive (PV⁺) inhibitory cells. Sparse firing of PYR cells and large excitatory currents onto PV⁺ cells are present as in experiments. The particular theta frequency is more controlled by PYR-to-PV⁺ cell interactions rather than PV⁺-to-PYR cell interactions. We identify two scenarios by which theta rhythms can emerge, and they can be differentiated by the ratio of excitatory to inhibitory currents to PV⁺ cells, but not to PYR cells. Only one of the scenarios is consistent with data from the whole hippocampus preparation, which leads to the prediction that the connection probability from PV⁺ to PYR cells needs to be larger than from PYR to PV⁺ cells. Our models can serve as a platform on which to build and develop an understanding of in vivo theta generation.