Cited by 198Neural dynamics and brain functionMemory and Neural MechanismsNeuroscience and Neuropharmacology Research

The DOI matches its Crossref record; 160 of the 160 references it cites were found in a registry or an index; none of the 87 sources its claims cite is checked for support yet.

Interneuronal mechanisms of hippocampal theta oscillations in a full-scale model of the rodent CA1 circuit

Open access

Published version of a preprint

Abstract

The hippocampal theta rhythm plays important roles in information processing; however, the mechanisms of its generation are not well understood. We developed a data-driven, supercomputer-based, full-scale (1:1) model of the rodent CA1 area and studied its interneurons during theta oscillations. Theta rhythm with phase-locked gamma oscillations and phase-preferential discharges of distinct interneuronal types spontaneously emerged from the isolated CA1 circuit without rhythmic inputs. Perturbation experiments identified parvalbumin-expressing interneurons and neurogliaform cells, as well as interneuronal diversity itself, as important factors in theta generation. These simulations reveal new insights into the spatiotemporal organization of the CA1 circuit during theta oscillations.

Overview

AI snapshot

Written from the paper’s full text.

TL;DR

A data-driven, full-scale model of the isolated rodent CA1 spontaneously generated theta (~7.8 Hz) and phase-locked gamma oscillations without rhythmic external inputs. Perturbation experiments implicated PV+ interneurons, neurogliaform cells, pyramidal recurrent connections, and interneuronal diversity itself in theta generation.

Key findings

  • The model contained 338,740 cells (311,500 pyramidal cells and 27,240 interneurons), 454,700 artificial stimulating afferents, and 5.19 billion synapses; it spontaneously produced stable theta at 7.8 Hz and gamma oscillations (25–80 Hz), with gamma amplitude phase-locked to the theta trough.
  • Distinct interneuron types showed phase-preferential firing: PV+ basket, bistratified, and O-LM cells fired near the theta trough, whereas ivy, neurogliaform, CCK+ basket, and axo-axonic cells fired least around the trough; 71% (5/7) of interneuron types with experimental data matched anesthetized in vivo phase maxima.
  • Muting pyramidal cell output or CA1 pyramidal-to-pyramidal connections collapsed theta; muting PV-family, CCK+ basket, or neurogliaform cell output strongly reduced theta, while muting O-LM, ivy, or SC-A cells had no effect. Partial PV+ muting disrupted theta more than partial SOM+ muting.
  • Reducing interneuron diversity abolished theta: giving all interneurons one electrophysiological profile, or morphing them toward PV+ basket cells, did not restore theta. Removing the GABAB component of neurogliaform synapses reduced theta power, but replacing it with a scaled-up fast GABAA conductance with equivalent charge transfer restored theta.

Methods

  • The CA1 network was implemented in parallel NEURON 7.4 with one multicompartmental pyramidal cell type and eight interneuron types; connectivity and cell numbers were taken from Bezaire and Soltesz (2013), and simulations used fixed time steps of 0.01–0.025 ms for 2000 or 4000 ms on supercomputers.
  • CA3 and entorhinal afferents were modeled as independent Poisson spike trains; control theta emerged at a mean afferent firing rate of 0.65 Hz, and analyses included spike density functions, Welch periodograms, spectrograms, an LFP analog, Rayleigh tests for phase modulation, and ANOVA for GABAB-related simulations.
  • Perturbation experiments muted each cell type’s output, altered pyramidal-to-pyramidal connection numbers, reduced interneuron diversity, and manipulated neurogliaform GABAB/GABAA signaling.

Limitations

  • Stated: about 30% of known CA1 interneurons were omitted due to paucity of quantitative data, likely making inhibition weaker than biological; pyramidal cells fired more than in vivo, and although scaled-up inhibition produced only 1% active pyramidal cells firing at 1.8 Hz, it was not used because the source of the additional inhibition was unidentified.
  • Stated: the model did not include neuromodulators, gap junctions, short-term plasticity, glial dynamics, pyramidal cell diversity, or cell-to-cell variability; the emergent theta likely resembled the atropine-resistant form, some single-cell electrophysiology came from mouse while the network was rat, and no isolated-CA1 data were available for interneuron phase preferences, so comparisons used anesthetized in vivo data.
  • Evident: theta appeared only within particular afferent excitation levels (around 0.65 Hz; lower or higher levels disrupted it), gamma was not examined further, and the authors note they could not rule out a hypothetical theta regime somewhere in the low-diversity parameter space.

Model-written summary. Check it against the paper before relying on it.

Integrity signals

Registry

Crossref · 10.7554/elife.18566

Cites

160 of 160 references verified0 citation errors · 0 not checkable

Statistics

Nothing to recompute.

Figures

Figure 1. CA1 network connectivity. (A) The model network is arranged in a layered prism with the lengths of each dimension similar to the actual dimensions of the CA1 region and its layers. (B) The model cell somata within a small chunk of stratum pyramidale (as depicted in A) are plotted to show the regular distribution of model cells throughout the layer in which they are found. (C) Each pyramidal cell in the network has detailed morphology with realistic incoming synapse placement along the dendrites and soma. (D,E) Diagrams illustrate connectivity between types of cells. (D) The network includes one principal cell type (pyramidal cells) and eight interneuron types. Cell types that may connect are linked by a line colored according to the presynaptic cell type. Most cell types can connect to most other cell types. Total number of cells of each type are displayed, as are the number of local output synapses (boutons) from all cells of each type. (E) The number, position, and cell types of each connection are biologically constrained, as are the numbers and positions of the cells. See Figure 1—figure supplement 1) for details about the convergence onto each cell type. Also see Table 1 and Figure 1—figure supplement 2 for information about the cell-type combinations of the 5 billion connections and the axonal distributions followed by each cell type, as well as detailed connectivity results at http://doi.org/10.6080/K05H7D60. DOI: http://dx.doi.org/10.7554/eLife.18566.002
Figure 1. CA1 network connectivity. (A) The model network is arranged in a layered prism with the lengths of each dimension similar to the actual dimensions of the CA1 region and its layers. (B) The model cell somata within a small chunk of stratum pyramidale (as depicted in A) are plotted to show the regular distribution of model cells throughout the layer in which they are found. (C) Each pyramidal cell in the network has detailed morphology with realistic incoming synapse placement along the dendrites and soma. (D,E) Diagrams illustrate connectivity between types of cells. (D) The network includes one principal cell type (pyramidal cells) and eight interneuron types. Cell types that may connect are linked by a line colored according to the presynaptic cell type. Most cell types can connect to most other cell types. Total number of cells of each type are displayed, as are the number of local output synapses (boutons) from all cells of each type. (E) The number, position, and cell types of each connection are biologically constrained, as are the numbers and positions of the cells. See Figure 1—figure supplement 1) for details about the convergence onto each cell type. Also see Table 1 and Figure 1—figure supplement 2 for information about the cell-type combinations of the 5 billion connections and the axonal distributions followed by each cell type, as well as detailed connectivity results at http://doi.org/10.6080/K05H7D60. DOI: http://dx.doi.org/10.7554/eLife.18566.002CC BY · PubMed Central
Figure 1—figure supplement 1. Quantitative network connectivity. The average number of incoming synapses per postsynaptic cell of the given type are shown for (A) all inputs to the cells, (B) all excitatory inputs to the cells and (C) all inhibitory inputs to the cells. DOI: http://dx.doi.org/10.7554/eLife.18566.003
Figure 1—figure supplement 1. Quantitative network connectivity. The average number of incoming synapses per postsynaptic cell of the given type are shown for (A) all inputs to the cells, (B) all excitatory inputs to the cells and (C) all inhibitory inputs to the cells. DOI: http://dx.doi.org/10.7554/eLife.18566.003CC BY · PubMed Central
Figure 1—figure supplement 2. Anatomically constrained connectivity. The axonal distributions are shown per presynaptic cell type. The distribution of boutons is plotted as a function of distance from the presynaptic cell’s soma. Boutons connecting to all possible types of postsynaptic cells are included in the plot. The colors correspond to each presynaptic cell type using the same color code as previous figures. DOI: http://dx.doi.org/10.7554/eLife.18566.004
Figure 1—figure supplement 2. Anatomically constrained connectivity. The axonal distributions are shown per presynaptic cell type. The distribution of boutons is plotted as a function of distance from the presynaptic cell’s soma. Boutons connecting to all possible types of postsynaptic cells are included in the plot. The colors correspond to each presynaptic cell type using the same color code as previous figures. DOI: http://dx.doi.org/10.7554/eLife.18566.004CC BY · PubMed Central
Figure 2. Electrophysiology of the model network components. (A) Ion channel densities vary as a function of location (top) in the morphologically detailed pyramidal cell model (bottom; adapted from Poolos et al., 2002). Scale bar: 100 μm and 0.01 μF/cm². (B–C) The sodium channel found in the pyramidal cell soma is characterized in terms of (B) the activation/inactivation curves and (C) the current-voltage relation at peak (transient) current and steady state. (D–G) Current sweeps are shown for four model cell types: (D) PV+ basket cell, (E) CCK+ basket cell, (F) O-LM cell, and (G) neurogliaform cell. Scale bar: 100 ms and 20 mV. (H–J) Electrophysiological properties for each cell type, including (H) input resistance, (I) membrane time constant, and (J) action potential threshold. (K–L) Pyramidal cell synaptic connections are characterized as post-synaptic currents with the postsynaptic cell voltage clamped at −50 mV; (K) synapses made onto the pyramidal cell from all other cell types and (L) synapses made by the pyramidal cell onto all network cell types. Cells represented by same colors as in Figure 1. Source Data available for electrophysiological characterizations shown here. Additional details available in the Methods, Table 3, and the Appendix. DOI: http://dx.doi.org/10.7554/eLife.18566.005
Figure 2. Electrophysiology of the model network components. (A) Ion channel densities vary as a function of location (top) in the morphologically detailed pyramidal cell model (bottom; adapted from Poolos et al., 2002). Scale bar: 100 μm and 0.01 μF/cm². (B–C) The sodium channel found in the pyramidal cell soma is characterized in terms of (B) the activation/inactivation curves and (C) the current-voltage relation at peak (transient) current and steady state. (D–G) Current sweeps are shown for four model cell types: (D) PV+ basket cell, (E) CCK+ basket cell, (F) O-LM cell, and (G) neurogliaform cell. Scale bar: 100 ms and 20 mV. (H–J) Electrophysiological properties for each cell type, including (H) input resistance, (I) membrane time constant, and (J) action potential threshold. (K–L) Pyramidal cell synaptic connections are characterized as post-synaptic currents with the postsynaptic cell voltage clamped at −50 mV; (K) synapses made onto the pyramidal cell from all other cell types and (L) synapses made by the pyramidal cell onto all network cell types. Cells represented by same colors as in Figure 1. Source Data available for electrophysiological characterizations shown here. Additional details available in the Methods, Table 3, and the Appendix. DOI: http://dx.doi.org/10.7554/eLife.18566.005CC BY · PubMed Central
Figure 3. Detailed network activity. (A–D) One second of network activity is shown. (A–B) The LFP analog, filtered at (A) the theta range of 5–10 Hz and (B) the low gamma range of 25–40 Hz, shows consistent theta and gamma signals. Scale bar represents 100 ms and 0.2 mV (theta) or 0.27 mV (gamma) for filtered LFP traces. (C) Raster of all spikes from cells within 100 μm of the reference electrode point. (D) Representative intracellular somatic membrane potential traces from cells near the reference electrode point. Scale bar represents 100 ms and 50 mV for the intracellular traces. DOI: http://dx.doi.org/10.7554/eLife.18566.037
Figure 3. Detailed network activity. (A–D) One second of network activity is shown. (A–B) The LFP analog, filtered at (A) the theta range of 5–10 Hz and (B) the low gamma range of 25–40 Hz, shows consistent theta and gamma signals. Scale bar represents 100 ms and 0.2 mV (theta) or 0.27 mV (gamma) for filtered LFP traces. (C) Raster of all spikes from cells within 100 μm of the reference electrode point. (D) Representative intracellular somatic membrane potential traces from cells near the reference electrode point. Scale bar represents 100 ms and 50 mV for the intracellular traces. DOI: http://dx.doi.org/10.7554/eLife.18566.037CC BY · PubMed Central
Figure 4. Spectral analysis of model activity. (A) A spectrogram of the local pyramidal-layer LFP analog (including contributions from all pyramidal cells within 100 μm of the reference electrode and 10% of pyramidal cells outside that radius) shows the stability and strength of the theta oscillation over time. The oscillation also featured strong harmonics at multiples of the theta frequency of 7.8 Hz. (B,D) Welch’s periodogram of the spike density function for each cell type, normalized by cell type and by displayed frequency range, shows the dominant network frequencies of (B) theta (7.8 Hz) and (D) gamma (71 Hz). Power is normalized to the peak power displayed in the power spectrum for each cell type. (C) Cross-frequency coupling between theta and gamma components of the LFP analog shows that the gamma oscillation is theta modulated. The gamma envelope is a function of the theta phase with the largest amplitude gamma oscillations occurring at the trough of the theta oscillation. Following convention, the theta trough was designated 0°/360°; see e.g., Varga et al. (2012). A graphical explanation of the relation between a spike train and its spike density function is shown in Figure 4—figure supplement 1. DOI: http://dx.doi.org/10.7554/eLife.18566.041
Figure 4. Spectral analysis of model activity. (A) A spectrogram of the local pyramidal-layer LFP analog (including contributions from all pyramidal cells within 100 μm of the reference electrode and 10% of pyramidal cells outside that radius) shows the stability and strength of the theta oscillation over time. The oscillation also featured strong harmonics at multiples of the theta frequency of 7.8 Hz. (B,D) Welch’s periodogram of the spike density function for each cell type, normalized by cell type and by displayed frequency range, shows the dominant network frequencies of (B) theta (7.8 Hz) and (D) gamma (71 Hz). Power is normalized to the peak power displayed in the power spectrum for each cell type. (C) Cross-frequency coupling between theta and gamma components of the LFP analog shows that the gamma oscillation is theta modulated. The gamma envelope is a function of the theta phase with the largest amplitude gamma oscillations occurring at the trough of the theta oscillation. Following convention, the theta trough was designated 0°/360°; see e.g., Varga et al. (2012). A graphical explanation of the relation between a spike train and its spike density function is shown in Figure 4—figure supplement 1. DOI: http://dx.doi.org/10.7554/eLife.18566.041CC BY · PubMed Central
Figure 4—figure supplement 1. Different views of cell activity. Several ways of characterizing model cell activity per cell type are shown using the spikes from the ivy cells as an example. (A) The spike times of each ivy cell are plotted as a function of time and ivy cell number. A subset of ivy cells positioned within 100 μm of the reference electrode location (whose spikes are shown in black) are then carried forward in the remaining calculations. (B) The spikes of the local ivy cells are binned into 1 ms windows to give spike counts per window. (C) A continuous representation of the ivy cell spikes as a function of time is given in the spike density function (SDF) computed from the ivy cell spike times. (D) A Welch’s Periodogram is computed, which summarizes the power of each oscillation frequency in the ivy cell SDF Although only a part of the simulation is shown, the full simulation length (except the first 50 ms) was used in the spectral analysis. DOI: http://dx.doi.org/10.7554/eLife.18566.044
Figure 4—figure supplement 1. Different views of cell activity. Several ways of characterizing model cell activity per cell type are shown using the spikes from the ivy cells as an example. (A) The spike times of each ivy cell are plotted as a function of time and ivy cell number. A subset of ivy cells positioned within 100 μm of the reference electrode location (whose spikes are shown in black) are then carried forward in the remaining calculations. (B) The spikes of the local ivy cells are binned into 1 ms windows to give spike counts per window. (C) A continuous representation of the ivy cell spikes as a function of time is given in the spike density function (SDF) computed from the ivy cell spike times. (D) A Welch’s Periodogram is computed, which summarizes the power of each oscillation frequency in the ivy cell SDF Although only a part of the simulation is shown, the full simulation length (except the first 50 ms) was used in the spectral analysis. DOI: http://dx.doi.org/10.7554/eLife.18566.044CC BY · PubMed Central
Figure 5. Model and experimental cell theta phases. All model results are based on the spiking of the cells within 100 μm of the reference electrode. (A–B) Firing probability by cell type as a function of theta phase for (A) model and (B) experimental cells under anesthesia (histograms adapted with permission from Figure 2, Figure 5B left, and Figure 6F respectively from Klausberger and Somogyi, 2008; Fuentealba et al., 2008; Fuentealba et al., 2010). The model histograms are normalized; see Figure 5—figure supplement 1 for firing rates. (C) Theta phase preference and theta modulation level were correlated; better modulated cell types spiked closer to the LFP analog trough near the phase preference of pyramidal cells. (D) Theta phase preference plotted on an idealized LFP wave for model data (base of arrow signifies the model phase preference and head of the arrow shows the distance to anesthetized, experimental phase preference). DOI: http://dx.doi.org/10.7554/eLife.18566.045
Figure 5. Model and experimental cell theta phases. All model results are based on the spiking of the cells within 100 μm of the reference electrode. (A–B) Firing probability by cell type as a function of theta phase for (A) model and (B) experimental cells under anesthesia (histograms adapted with permission from Figure 2, Figure 5B left, and Figure 6F respectively from Klausberger and Somogyi, 2008; Fuentealba et al., 2008; Fuentealba et al., 2010). The model histograms are normalized; see Figure 5—figure supplement 1 for firing rates. (C) Theta phase preference and theta modulation level were correlated; better modulated cell types spiked closer to the LFP analog trough near the phase preference of pyramidal cells. (D) Theta phase preference plotted on an idealized LFP wave for model data (base of arrow signifies the model phase preference and head of the arrow shows the distance to anesthetized, experimental phase preference). DOI: http://dx.doi.org/10.7554/eLife.18566.045CC BY · PubMed Central
Figure 5—figure supplement 1. Firing rates of model and experimental cells of each type. For experimental cells, firing rates in both the anesthetized and awake states were included where available. See Table 6 for sources of experimental data. DOI: http://dx.doi.org/10.7554/eLife.18566.057
Figure 5—figure supplement 1. Firing rates of model and experimental cells of each type. For experimental cells, firing rates in both the anesthetized and awake states were included where available. See Table 6 for sources of experimental data. DOI: http://dx.doi.org/10.7554/eLife.18566.057CC BY · PubMed Central
Figure 5—figure supplement 2. Theta phase-specific firing preferences of various biological hippocampal cell types as reported in the literature. The trough of the pyramidal-layer LFP is designated as 0/360 and the peak as 180. There is variation in phase preference for given cell types as a function of the experimental preparation. Shown are (A) anesthetized and (B) awake experimental conditions. Reference subscripts correspond to: 1: Klausberger et al. (2003), 2: Klausberger et al. (2004), 3: Klausberger et al. (2005), 4: Lapray et al. (2012), 5: Varga et al. (2012), 6: Fuentealba et al. (2008), 7: Fuentealba et al. (2010), 8: Varga et al. (2014). See Table 6 for further details. DOI: http://dx.doi.org/10.7554/eLife.18566.058
Figure 5—figure supplement 2. Theta phase-specific firing preferences of various biological hippocampal cell types as reported in the literature. The trough of the pyramidal-layer LFP is designated as 0/360 and the peak as 180. There is variation in phase preference for given cell types as a function of the experimental preparation. Shown are (A) anesthetized and (B) awake experimental conditions. Reference subscripts correspond to: 1: Klausberger et al. (2003), 2: Klausberger et al. (2004), 3: Klausberger et al. (2005), 4: Lapray et al. (2012), 5: Varga et al. (2012), 6: Fuentealba et al. (2008), 7: Fuentealba et al. (2010), 8: Varga et al. (2014). See Table 6 for further details. DOI: http://dx.doi.org/10.7554/eLife.18566.058CC BY · PubMed Central
Figure 6. Altered network configurations. Oscillation power (in mV2²/Hz) of the spike density function (SDF) for pyramidal cells within 100 μm of the reference electrode, at the peak frequency within theta range (5–10 Hz) in altered network configurations. For corresponding peak frequencies, see Figure 6—figure supplement 1. (A) Theta is present at some excitation levels. (B) Muting each cell type’s output caused a range of effects. (C) The stability and frequency of spontaneous theta in the network was sensitive to the presence and number of recurrent connections between CA1 pyramidal cells. (D) Partially muting the broad classes of PV+ or SOM+ cells by 50% showed that PV+ muting disrupted the network more than SOM+ muting. (E) Theta falls apart when all interneurons are given the same electrophysiological profile, whether it be of a PV+ basket, CCK+ basket, neurogliaform, or O-LM cell. (F) Gradually setting all interneuron properties to those of PV+ basket cells did not restore theta. From left to right: control network; PV+ basket cell electrophysiology; also weights of incoming synapses; also numbers of incoming synapses; then all interneurons being PV+ basket cells (with the addition of the output synapse numbers, weights, and kinetics); then variable RMP (normal distribution with standard deviation of 8 mV). (G) A wide range in excitation was unable to produce theta in the PV+ B. network. (H) Removing the GABAB component from the neurogliaform synapses onto other neurogliaform cells and pyramidal cells showed a significant drop in theta power. Massively increasing the weight of the GABAA component to produce a similar amount of charge transfer restored theta power (compare the IPSCs corresponding to each condition in Figure 6—figure supplement 2). Standard deviations (n = 3) shown; significance (p=1.8e-05). DOI: http://dx.doi.org/10.7554/eLife.18566.061
Figure 6. Altered network configurations. Oscillation power (in mV2²/Hz) of the spike density function (SDF) for pyramidal cells within 100 μm of the reference electrode, at the peak frequency within theta range (5–10 Hz) in altered network configurations. For corresponding peak frequencies, see Figure 6—figure supplement 1. (A) Theta is present at some excitation levels. (B) Muting each cell type’s output caused a range of effects. (C) The stability and frequency of spontaneous theta in the network was sensitive to the presence and number of recurrent connections between CA1 pyramidal cells. (D) Partially muting the broad classes of PV+ or SOM+ cells by 50% showed that PV+ muting disrupted the network more than SOM+ muting. (E) Theta falls apart when all interneurons are given the same electrophysiological profile, whether it be of a PV+ basket, CCK+ basket, neurogliaform, or O-LM cell. (F) Gradually setting all interneuron properties to those of PV+ basket cells did not restore theta. From left to right: control network; PV+ basket cell electrophysiology; also weights of incoming synapses; also numbers of incoming synapses; then all interneurons being PV+ basket cells (with the addition of the output synapse numbers, weights, and kinetics); then variable RMP (normal distribution with standard deviation of 8 mV). (G) A wide range in excitation was unable to produce theta in the PV+ B. network. (H) Removing the GABAB component from the neurogliaform synapses onto other neurogliaform cells and pyramidal cells showed a significant drop in theta power. Massively increasing the weight of the GABAA component to produce a similar amount of charge transfer restored theta power (compare the IPSCs corresponding to each condition in Figure 6—figure supplement 2). Standard deviations (n = 3) shown; significance (p=1.8e-05). DOI: http://dx.doi.org/10.7554/eLife.18566.061CC BY · PubMed Central
Figure 6—figure supplement 1. Peak frequencies of oscillations in altered networks. Peak theta frequency (within 5–10 Hz) of the spike density function (SDF) for all pyramidal cells within 100 μm of the reference electrode in each altered network configuration. For networks where no pyramidal cells spiked, resulting in zero power within the spectral analysis of the pyramidal cell spike density function, their peak frequencies are listed as ‘not available’ or ‘n/a’. (A) Spontaneous theta oscillation accelerated out of theta range with more excitation. (B) Muting each cell type shifted the oscillation out of range (neurogliaform, CCK+ basket, and axo-axonic cells), disrupted theta but not gamma (not shown; pyramidal, PV+ basket, and bistratified cells), or had little effect (S.C.-A., O-LM, and ivy cells). (C) Doubling the connections between CA1 pyramidal cells increased the theta frequency, while networks with half the number or no recurrent collaterals lost the slow oscillation but kept gamma. (D) Removing 50% of PV+ cell inhibition (PV+ basket, bistratified, and axo-axonic cells) or 50% of SOM+ cell inhibition (bistratified or O-LM cells) shifted the oscillation out of theta range or lost the slow oscillation entirely but kept gamma. (E) Peak oscillation shifted out of theta range when all interneurons had the same electrophysiological profile, regardless of the profile used. (F) Converging all properties to PV+ basket cells, gamma was restored (not shown) but not theta (left to right: control; network with 1: diverse interneurons with same electrophysiology; 2: also with same weights of incoming synapses; 3: also with same numbers of incoming synapses; 4: complete conversion to PV+ basket cells; 5: added variability in resting membrane potential (normal distribution with st. dev.=8 mV)). (G) In the all-PV+ basket cell network, a wide range of excitation levels could not produce a spontaneous theta rhythm. (H) Removing GABAB increased the oscillation frequency. DOI: http://dx.doi.org/10.7554/eLife.18566.064
Figure 6—figure supplement 1. Peak frequencies of oscillations in altered networks. Peak theta frequency (within 5–10 Hz) of the spike density function (SDF) for all pyramidal cells within 100 μm of the reference electrode in each altered network configuration. For networks where no pyramidal cells spiked, resulting in zero power within the spectral analysis of the pyramidal cell spike density function, their peak frequencies are listed as ‘not available’ or ‘n/a’. (A) Spontaneous theta oscillation accelerated out of theta range with more excitation. (B) Muting each cell type shifted the oscillation out of range (neurogliaform, CCK+ basket, and axo-axonic cells), disrupted theta but not gamma (not shown; pyramidal, PV+ basket, and bistratified cells), or had little effect (S.C.-A., O-LM, and ivy cells). (C) Doubling the connections between CA1 pyramidal cells increased the theta frequency, while networks with half the number or no recurrent collaterals lost the slow oscillation but kept gamma. (D) Removing 50% of PV+ cell inhibition (PV+ basket, bistratified, and axo-axonic cells) or 50% of SOM+ cell inhibition (bistratified or O-LM cells) shifted the oscillation out of theta range or lost the slow oscillation entirely but kept gamma. (E) Peak oscillation shifted out of theta range when all interneurons had the same electrophysiological profile, regardless of the profile used. (F) Converging all properties to PV+ basket cells, gamma was restored (not shown) but not theta (left to right: control; network with 1: diverse interneurons with same electrophysiology; 2: also with same weights of incoming synapses; 3: also with same numbers of incoming synapses; 4: complete conversion to PV+ basket cells; 5: added variability in resting membrane potential (normal distribution with st. dev.=8 mV)). (G) In the all-PV+ basket cell network, a wide range of excitation levels could not produce a spontaneous theta rhythm. (H) Removing GABAB increased the oscillation frequency. DOI: http://dx.doi.org/10.7554/eLife.18566.064CC BY · PubMed Central
Figure 6—figure supplement 2. IPSCs from the neurogliaform to pyramidal cell synapse corresponding to the different conditions in Figure 6H. These traces are from pyramidal cells clamped at −50 mV during a paired recording from a presynaptic neurogliaform cell with a GABAA reversal potential of −60 mV and a GABAB reversal potential of −90 mV. The currents shown are averages from 10 recordings. Scale bar = 100 ms and 5 pA. DOI: http://dx.doi.org/10.7554/eLife.18566.065
Figure 6—figure supplement 2. IPSCs from the neurogliaform to pyramidal cell synapse corresponding to the different conditions in Figure 6H. These traces are from pyramidal cells clamped at −50 mV during a paired recording from a presynaptic neurogliaform cell with a GABAA reversal potential of −60 mV and a GABAB reversal potential of −90 mV. The currents shown are averages from 10 recordings. Scale bar = 100 ms and 5 pA. DOI: http://dx.doi.org/10.7554/eLife.18566.065CC BY · PubMed Central
Appendix 1—figure 1. Firing Rates of Experimental Cells. Rebound spiking, which occurs in some O-LM cells at hyperpolarized current injection levels, is not shown in this graph. DOI: http://dx.doi.org/10.7554/eLife.18566.067
Appendix 1—figure 1. Firing Rates of Experimental Cells. Rebound spiking, which occurs in some O-LM cells at hyperpolarized current injection levels, is not shown in this graph. DOI: http://dx.doi.org/10.7554/eLife.18566.067CC BY · PubMed Central
Appendix 1—figure 2. Physiological properties of experimental and model cells. Experimental data are shown with closed markers for the mean and error bars for cell types where n > 1. The model cell properties are plotted as open circles. Calculation of properties is explained in the text. (A) resting membrane potential, (B) threshold, and (C) spike amplitude. DOI: http://dx.doi.org/10.7554/eLife.18566.070
Appendix 1—figure 2. Physiological properties of experimental and model cells. Experimental data are shown with closed markers for the mean and error bars for cell types where n > 1. The model cell properties are plotted as open circles. Calculation of properties is explained in the text. (A) resting membrane potential, (B) threshold, and (C) spike amplitude. DOI: http://dx.doi.org/10.7554/eLife.18566.070CC BY · PubMed Central
Appendix 1—figure 3. Physiological properties, continued. (A) sag time constant, (B) sag amplitude, and (C) amplitude of afterhyperpolarization (AHP). DOI: http://dx.doi.org/10.7554/eLife.18566.071
Appendix 1—figure 3. Physiological properties, continued. (A) sag time constant, (B) sag amplitude, and (C) amplitude of afterhyperpolarization (AHP). DOI: http://dx.doi.org/10.7554/eLife.18566.071CC BY · PubMed Central
Appendix 1—figure 4. Physiological properties, continued. (A) rheobase, (B) membrane time constant, (C) interspike interval (ISI), and (D) input resistance. DOI: http://dx.doi.org/10.7554/eLife.18566.072
Appendix 1—figure 4. Physiological properties, continued. (A) rheobase, (B) membrane time constant, (C) interspike interval (ISI), and (D) input resistance. DOI: http://dx.doi.org/10.7554/eLife.18566.072CC BY · PubMed Central
Appendix 1—figure 5. Pyramidal (A) model and (B) experimental current sweep. (C) Firing rates of model and experimental cells. DOI: http://dx.doi.org/10.7554/eLife.18566.073
Appendix 1—figure 5. Pyramidal (A) model and (B) experimental current sweep. (C) Firing rates of model and experimental cells. DOI: http://dx.doi.org/10.7554/eLife.18566.073CC BY · PubMed Central
Appendix 1—figure 6. Connections onto (A) and (B) from model Pyramidal cells, under voltage clamp at −50 mV with physiological reversal potentials. DOI: http://dx.doi.org/10.7554/eLife.18566.081
Appendix 1—figure 6. Connections onto (A) and (B) from model Pyramidal cells, under voltage clamp at −50 mV with physiological reversal potentials. DOI: http://dx.doi.org/10.7554/eLife.18566.081CC BY · PubMed Central
Appendix 1—figure 7. Axo-axonic (A) model and (B) experimental current sweep. (C) Firing rates of model and experimental cells. DOI: http://dx.doi.org/10.7554/eLife.18566.082
Appendix 1—figure 7. Axo-axonic (A) model and (B) experimental current sweep. (C) Firing rates of model and experimental cells. DOI: http://dx.doi.org/10.7554/eLife.18566.082CC BY · PubMed Central
Appendix 1—figure 8. Connections onto (A) and (B) from model Axo-axonic cells, under voltage clamp at −50 mV with physiological reversal potentials. DOI: http://dx.doi.org/10.7554/eLife.18566.089
Appendix 1—figure 8. Connections onto (A) and (B) from model Axo-axonic cells, under voltage clamp at −50 mV with physiological reversal potentials. DOI: http://dx.doi.org/10.7554/eLife.18566.089CC BY · PubMed Central
Appendix 1—figure 9. Bistratified (A) model and (B) experimental current sweep. (C) Firing rates of model and experimental cells. DOI: http://dx.doi.org/10.7554/eLife.18566.090
Appendix 1—figure 9. Bistratified (A) model and (B) experimental current sweep. (C) Firing rates of model and experimental cells. DOI: http://dx.doi.org/10.7554/eLife.18566.090CC BY · PubMed Central
Appendix 1—figure 10. Connections onto (A) and (B) from model Bistratified cells, under voltage clamp at −50 mV with physiological reversal potentials. DOI: http://dx.doi.org/10.7554/eLife.18566.098
Appendix 1—figure 10. Connections onto (A) and (B) from model Bistratified cells, under voltage clamp at −50 mV with physiological reversal potentials. DOI: http://dx.doi.org/10.7554/eLife.18566.098CC BY · PubMed Central
Appendix 1—figure 11. CCK+ Basket (A) model and (B) experimental current sweep. (C) Firing rates of model and experimental cells. DOI: http://dx.doi.org/10.7554/eLife.18566.099
Appendix 1—figure 11. CCK+ Basket (A) model and (B) experimental current sweep. (C) Firing rates of model and experimental cells. DOI: http://dx.doi.org/10.7554/eLife.18566.099CC BY · PubMed Central
Appendix 1—figure 12. Connections onto (A) and (B) from model CCK+ Basket cells, under voltage clamp at −50 mV with physiological reversal potentials. DOI: http://dx.doi.org/10.7554/eLife.18566.106
Appendix 1—figure 12. Connections onto (A) and (B) from model CCK+ Basket cells, under voltage clamp at −50 mV with physiological reversal potentials. DOI: http://dx.doi.org/10.7554/eLife.18566.106CC BY · PubMed Central
Appendix 1—figure 13. Ivy (A) model and (B) experimental current sweep. (fig:ivypage:firing) Firing rates of model and experimental cells. DOI: http://dx.doi.org/10.7554/eLife.18566.107
Appendix 1—figure 13. Ivy (A) model and (B) experimental current sweep. (fig:ivypage:firing) Firing rates of model and experimental cells. DOI: http://dx.doi.org/10.7554/eLife.18566.107CC BY · PubMed Central
Appendix 1—figure 14. Connections onto (A) and (B) from model Ivy cells, under voltage clamp at −50 mV with physiological reversal potentials. DOI: http://dx.doi.org/10.7554/eLife.18566.115
Appendix 1—figure 14. Connections onto (A) and (B) from model Ivy cells, under voltage clamp at −50 mV with physiological reversal potentials. DOI: http://dx.doi.org/10.7554/eLife.18566.115CC BY · PubMed Central
Appendix 1—figure 15. Neurogliaform (A) model and (B) experimental current sweep. (C) Firing rates of model and experimental cells. DOI: http://dx.doi.org/10.7554/eLife.18566.116
Appendix 1—figure 15. Neurogliaform (A) model and (B) experimental current sweep. (C) Firing rates of model and experimental cells. DOI: http://dx.doi.org/10.7554/eLife.18566.116CC BY · PubMed Central
Appendix 1—figure 16. Connections onto (A) and (B) from model Neurogliaform cells, under voltage clamp at −50 mV with physiological reversal potentials. DOI: http://dx.doi.org/10.7554/eLife.18566.124
Appendix 1—figure 16. Connections onto (A) and (B) from model Neurogliaform cells, under voltage clamp at −50 mV with physiological reversal potentials. DOI: http://dx.doi.org/10.7554/eLife.18566.124CC BY · PubMed Central
Appendix 1—figure 17. O-LM (A) model and (B) experimental current sweep. (C) Firing rates of model and experimental cells. DOI: http://dx.doi.org/10.7554/eLife.18566.125
Appendix 1—figure 17. O-LM (A) model and (B) experimental current sweep. (C) Firing rates of model and experimental cells. DOI: http://dx.doi.org/10.7554/eLife.18566.125CC BY · PubMed Central
Appendix 1—figure 18. Connections onto (A) and (B) from model O-LM cells, under voltage clamp at −50 mV with physiological reversal potentials. DOI: http://dx.doi.org/10.7554/eLife.18566.132
Appendix 1—figure 18. Connections onto (A) and (B) from model O-LM cells, under voltage clamp at −50 mV with physiological reversal potentials. DOI: http://dx.doi.org/10.7554/eLife.18566.132CC BY · PubMed Central
Appendix 1—figure 19. PV+ Basket (A) model and (B) experimental current sweep. (C) Firing rates of model and experimental cells. DOI: http://dx.doi.org/10.7554/eLife.18566.133
Appendix 1—figure 19. PV+ Basket (A) model and (B) experimental current sweep. (C) Firing rates of model and experimental cells. DOI: http://dx.doi.org/10.7554/eLife.18566.133CC BY · PubMed Central
Appendix 1—figure 20. Connections onto (A) and (B) from model PV+ Basket cells, under voltage clamp at −50 mV with physiological reversal potentials. DOI: http://dx.doi.org/10.7554/eLife.18566.141
Appendix 1—figure 20. Connections onto (A) and (B) from model PV+ Basket cells, under voltage clamp at −50 mV with physiological reversal potentials. DOI: http://dx.doi.org/10.7554/eLife.18566.141CC BY · PubMed Central
Appendix 1—figure 21. Schaffer Collateral-Associated (A) model and (B) experimental current sweep. (C) Firing rates of model and experimental cells. DOI: http://dx.doi.org/10.7554/eLife.18566.142
Appendix 1—figure 21. Schaffer Collateral-Associated (A) model and (B) experimental current sweep. (C) Firing rates of model and experimental cells. DOI: http://dx.doi.org/10.7554/eLife.18566.142CC BY · PubMed Central
Appendix 1—figure 22. Connections onto (A) and (B) from model Schaffer Collateral-Associated cells, under voltage clamp at −50 mV with physiological reversal potentials. DOI: http://dx.doi.org/10.7554/eLife.18566.150
Appendix 1—figure 22. Connections onto (A) and (B) from model Schaffer Collateral-Associated cells, under voltage clamp at −50 mV with physiological reversal potentials. DOI: http://dx.doi.org/10.7554/eLife.18566.150CC BY · PubMed Central
Appendix 1—figure 23. Calcium channel currents. DOI: http://dx.doi.org/10.7554/eLife.18566.156
Appendix 1—figure 23. Calcium channel currents. DOI: http://dx.doi.org/10.7554/eLife.18566.156CC BY · PubMed Central
Appendix 1—figure 24. HCN channel currents. DOI: http://dx.doi.org/10.7554/eLife.18566.157
Appendix 1—figure 24. HCN channel currents. DOI: http://dx.doi.org/10.7554/eLife.18566.157CC BY · PubMed Central
Appendix 1—figure 25. Delayed rectifier potassium channel currents. DOI: http://dx.doi.org/10.7554/eLife.18566.158
Appendix 1—figure 25. Delayed rectifier potassium channel currents. DOI: http://dx.doi.org/10.7554/eLife.18566.158CC BY · PubMed Central
Appendix 1—figure 26. A-type potassium channel currents. DOI: http://dx.doi.org/10.7554/eLife.18566.159
Appendix 1—figure 26. A-type potassium channel currents. DOI: http://dx.doi.org/10.7554/eLife.18566.159CC BY · PubMed Central
Appendix 1—figure 27. Other potassium channel currents. Because they didn’t have a voltage-sensitive inactivation component, only the activation curve, which is equivalent to the IV Peak curve, need be shown here. DOI: http://dx.doi.org/10.7554/eLife.18566.160
Appendix 1—figure 27. Other potassium channel currents. Because they didn’t have a voltage-sensitive inactivation component, only the activation curve, which is equivalent to the IV Peak curve, need be shown here. DOI: http://dx.doi.org/10.7554/eLife.18566.160CC BY · PubMed Central
Appendix 1—figure 28. Calcium-dependent potassium channel dependence on calcium concentration. (a) The normalized conductance of the channels are plotted as a function of test voltage step and calcium concentration. (b) and (c) The current-voltage relation is shown at several calcium concentrations for (b) KvCaB channel and (c) KCaS channel. Note that the KCaS channel is only active at the highest calcium concentration and is not dependent on voltage (although the voltage continues to set the driving force) when it is active. DOI: http://dx.doi.org/10.7554/eLife.18566.161
Appendix 1—figure 28. Calcium-dependent potassium channel dependence on calcium concentration. (a) The normalized conductance of the channels are plotted as a function of test voltage step and calcium concentration. (b) and (c) The current-voltage relation is shown at several calcium concentrations for (b) KvCaB channel and (c) KCaS channel. Note that the KCaS channel is only active at the highest calcium concentration and is not dependent on voltage (although the voltage continues to set the driving force) when it is active. DOI: http://dx.doi.org/10.7554/eLife.18566.161CC BY · PubMed Central
Appendix 1—figure 29. Sodium channel voltage dependence. The normalized conductance of the sodium channel is plotted as (a) a function of test voltage step to show activation and (b) as a function of holding voltage prior to the test step to show inactivation. DOI: http://dx.doi.org/10.7554/eLife.18566.162
Appendix 1—figure 29. Sodium channel voltage dependence. The normalized conductance of the sodium channel is plotted as (a) a function of test voltage step to show activation and (b) as a function of holding voltage prior to the test step to show inactivation. DOI: http://dx.doi.org/10.7554/eLife.18566.162CC BY · PubMed Central

Verification

Claims the paper makes

53 Not checked yet · 46 Own result · 29 Own interpretation · 9 Own method · 6 No basis given · 1 Sources not found

Verify all 86 references — free account

Show all 91 claims

  • It may state two separate claims as one. 23
  • Probably not a substantial claim. 37
53 of 53
  • The paper relies on this without support.

    The CA1 area is one of the most extensively studied brain regions, with abundant anatomical and electrophysiological data about its organization.

    “Because the CA1 area is one of the most extensively studied brain regions, there are abundant anatomical and electrophysiological data about its organization, making it a logical choice for the development of a full-scale model.”
  • The paper relies on this without support.

    The model parameters that could not be constrained by experimental data do not critically affect the model's conclusions.

    “In addition, the remaining parameters that could not be constrained by experimental data were documented, with the assumptions used to arrive at them explicitly listed in Table 2 of Bezaire and Soltesz (2013) and additional parameter calculations described in the Appendix of the present paper, section 'Inhibitory connectivity'.”
  • The paper relies on this without support.

    The missing 30% of interneurons do not qualitatively alter the generation of theta oscillations in the model.

    “The remaining 30% of the interneurons were not included in the model due to paucity of quantitative data (Bezaire and Soltesz, 2013).”
  • The paper relies on this without support.

    The experimental phase preferences from anesthetized rats are comparable to the model's conditions.

    “Firing probability by cell type as a function of theta phase for (A) model and (B) experimental cells under anesthesia (histograms adapted with permission from Figure 2, Figure 5B left, and Figure 6F respectively from Klausberger and Somogyi, 2008; Fuentealba et al., 2008; Fuentealba et al., 2010).”
  • The paper relies on this without support.

    The phase preferences of interneurons in anesthetized animals are comparable to those in the model's conditions for the purpose of validating the model.

    “because no data are available on the phase preferential firing of morphologically identified interneurons from the isolated CA1 preparation, comparison is made here with results from anesthetized animals, from which the most complete data sets are available”
  • The paper relies on this without support.

    The parameter space explored for low-diversity configurations is sufficient to conclude that interneuronal diversity is necessary for spontaneous theta oscillations.

    “While we could not rule out a hypothetical theta regime somewhere in the parameter space of such low-diversity configurations, any theta solution space would likely be smaller and more elusive than we were able to determine in the control configuration (Figure 6A).”
  • The paper relies on this without support.

    The phase of LFP theta recorded in the stratum radiatum is approximately 180 degrees out of phase with the pyramidal cell layer theta used in this paper.

    “note that Ferguson and colleagues used LFP theta recorded in the stratum radiatum as reference, which is approximately 180 degrees out of phase with the pyramidal cell layer theta used in this paper”
  • The paper relies on this without support.

    Model neurons sometimes behave much differently than expected when subjected to current sweep protocols or synaptic inputs that are outside the range of the original protocols used to construct the model.

    “Model neurons sometimes behave much differently than expected when subjected to current sweep protocols or synaptic inputs that are outside the range of the original protocols used to construct the model”
  • The authors state this as an assumption.

    All available incoming inhibitory synapses onto interneurons in a layer have an equal chance of being innervated by the available inhibitory boutons targeting interneurons in that layer.

    “All available incoming inhibitory synapses onto interneurons in a layer have an equal chance of being innervated by the available inhibitory boutons targeting interneurons in that layer”
  • The paper relies on this without support.

    There is no reported evidence for major, systematic differences in key parameters such as the phase specific firing of rat and mouse interneurons in vivo.

    “there is no reported evidence for major, systematic differences in key parameters such as the phase specific firing of rat and mouse interneurons in vivo”
10 of 10

Claims other papers make about this one

Each is checked against this paper’s own text: does it say what the citing paper says it does?

Citations

Cites

Show all 160 references

Cited by

Sorted by citations

Metadata

Authors
Bezaire, M. J., Raikov, I., Burk, K., Vyas, D., & Soltesz, I.
Affiliations
  • University of California, Irvine (US)
  • Department of Neurosurgery, Stanford University, Stanford, United States
Venue
eLife
Published
2016-12-23
Volume
5
ISSN
2050-084X (ISSN-L)
OpenAlex
W2561068888
Cited by
198
References
160
Citations by year
2026: 8 · 2025: 16 · 2024: 27 · 2023: 26 · 2022: 32 · 2021: 25 · 2020: 15 · 2019: 28 · 2018: 15 · 2017: 5 · 2016: 1
FWCI
10.39
Type
journal article
Language
en
Open access
Locations
  • eLife (publishedVersion, cc-by)
  • bioRxiv (Cold Spring Harbor Laboratory) (acceptedVersion, cc-by)
  • Europe PMC (PubMed Central) (publishedVersion, cc-by)
  • PubMed Central (publishedVersion, cc-by)
  • PubMed (publishedVersion)
  • DOAJ (DOAJ: Directory of Open Access Journals) (submittedVersion, cc-by-sa)
Funders
  • National Institutes of Health (F32NS090753, NS35915, NS090583)
  • National Science Foundation (DGE-0808392, XSEDE Allocations - TG-IBN140007, XSEDE Allocations - TG-IBN130022, XSEDE Allocations - TG-IBN100011, IOS-1310378)
Topics
  • Life Sciences › Neuroscience › Cognitive Neuroscience › Neural dynamics and brain function
  • Life Sciences › Neuroscience › Cognitive Neuroscience › Memory and Neural Mechanisms
  • Life Sciences › Neuroscience › Cellular and Molecular Neuroscience › Neuroscience and Neuropharmacology Research
Keywords
hippocampal theta rhythm, interneurons, parvalbumin interneurons, neurogliaform cells, gamma oscillations, phase locking, computational neuroscience, spontaneous oscillations
MeSH terms
Animals, Interneurons, Neuroglia, Theta Rhythm, Neural Networks, Computer, Mice, Rats, CA1 Region, Hippocampal, Spatio-Temporal Analysis

About these checks

Registry, retraction, Cites and statistics results are computed, not judged, with one exception: a printed reference that no index lists is read once by a model before it counts as a citation error. Claims are extracted by a model and kept only when their quote appears verbatim in the text; each cited source is then checked, with a model’s help, for whether it says what the claim relies on. Those verdicts can be wrong — read the claim’s quote and the cited source before relying on one.

A language model drew the claims from the paper on 2026-10-09, under version 3 of the claim definition (model deepseek-v4.1-flash, prompt extract-v3).

A second model labelled each claim on 2026-10-09 (model jev-1.13-20260917, definition revision gold-guideline-v0.2).

Statistics and arithmetic are recomputed, not judged.

Confidence bands are shown only for fields calibrated on our annotated test set. None is calibrated yet, so no band is shown.

What counts as a claim ›How we measure confidence ›How we check this ›

Are you an author, or do you see a problem with this page? Request a review