Abstract
Synchronous activity of neuronal networks is found in many brain areas and correlates with cognition and behavior. Gamma synchrony is particularly strong in the dentate gyrus, which is thought to process contextual information in the hippocampus. Several network mechanisms for synchrony generation have been proposed and studied computationally. One such mechanism relies solely on recurrent inhibitory interneuron connectivity, but it requires a large enough number of synapses. Here, we incorporate previously published connectivity data of the dentate gyrus from mice of either sex into a biophysical computational model to test its ability to generate synchronous activity. We find that recurrent interneuron connectivity is insufficient to induce synchronous activity. This applies to an interneuron ring network and the broader dentate gyrus circuitry. Despite asynchronous input, recurrent interneuron connectivity can have small synchronizing effects but can also desynchronize the network for some types of synaptic input. Our results suggest that biologically plausible recurrent inhibitory connectivity alone is likely insufficient to synchronize the dentate gyrus.
Significance Statement
Neurons in the brain do not activate randomly but show synchronous activity with other neurons during states of high activity. In the hippocampus, a brain area responsible for memory storage, synchronous activity is well known and the recurrent inhibitory connections have been proposed to synchronize neurons. Here, we simulated a detailed model of a hippocampal brain area called the dentate gyrus with modern connectivity estimates. We found that the number of connections is insufficient to synchronize the network. This means that other models are more likely to explain synchronous activity in the dentate gyrus. Furthermore, we predict that future experiments interfering with recurrent connectivity will not affect synchronous activity.
Introduction
Synchronous neuronal activity is important for a range of brain functions. Synchrony can occur at different frequencies, and the gamma range (30 − 100 Hz) has been implied by correlation and experiment with cognition and behavior (Fries, 2009; Girardeau et al., 2009; Wang, 2010; Kitanishi et al., 2015). Therefore, synchronous activity plays an important role in brain disorders such as epilepsies and Parkinson’s disease (Uhlhaas and Singer, 2006; Hammond et al., 2007). In the dentate gyrus (DG), gamma oscillations are particularly strong (Buzsáki and Vanderwolf, 1983). Many mechanisms come together to create these gamma oscillations. However, the postsynaptic currents of synchronously active interneurons (INs) are a crucial component (Buzsáki and Wang, 2012).
Several network mechanisms that could generate synchrony have been proposed, the interneuron gamma generator (ING) being one of them. Experiments show that inhibitory INs have a central role in the generation of gamma oscillations (Bartos et al., 2007; Wulff et al., 2009; Buzsáki and Wang, 2012; Allen and Monyer, 2015; Antonoudiou et al., 2020). Two major circuit models could allow INs to generate synchrony: the ING model and the pyramidal-interneuron gamma (PING) model (Tiesinga and Sejnowski, 2009). In the ING model, INs synchronize themselves during asynchronous input solely through recurrent inhibition (Wang and Buzsáki, 1996; Bartos et al., 2002). In the PING model, on the other hand, the network becomes synchronized through the reciprocal connections of principal cells and INs. The ING and the PING models both work well in-silico and have shown merit in explaining experimental results in some brain areas. However, the ING mechanism requires dense recurrent inhibitory connections, where the density depends on intrinsic and synaptic properties (Wang and Buzsáki, 1996; Golomb and Hansel, 2000). Therefore, determining whether the ING mechanism is biologically plausible in a given area requires a careful analysis of the interplay between connectivity and intrinsic properties. We use a biophysical conductance-based model with simplified dendrites to test the ING hypothesis under biologically detailed conditions.
The DG is a hippocampal brain region that receives input primarily from the entorhinal cortex (EC) and sends its output to CA3 (Amaral et al., 2007). It is essential for pattern separation (Yassa and Stark, 2011) and has the largest gamma power among the hippocampal regions (Buzsáki and Vanderwolf, 1983). The DG has been proposed to act as a gamma generator. However, unlike the CA3 generator, it is not independent of EC input (Csicsvari et al., 2003). Since the EC is already gamma frequency modulated, it is unclear to what extent the DG generates its own gamma rhythm or whether it merely relays EC gamma. At least at theta frequencies, the EC strongly drives the DG (Mizuseki et al., 2009). To restore healthy brain rhythms in brain disorders, the distinction between internally and externally generated oscillations is crucial.
IN types and circuitry of the DG have been extensively studied, and parvalbumin-positive (PV+) INs in particular, have revealed extensive recurrent connections (both electrical and chemical) with other PV+ INs (Espinoza et al., 2018). This connectivity, together with their intrinsic properties and their fast signaling (Bartos et al., 2002), in theory, makes them an excellent substrate for an ING mechanism. However, optogenetic experiments have not been performed on DG PV+ INs while measuring synchronous activity and such experiments would not distinguish between ING and PING mechanisms, unless the results could be integrated into a computational model (Tiesinga, 2012). In practice, it is therefore unclear whether PV+ IN connectivity synchronizes the network. We have thus incorporated the PV+ IN connectivity as measured by Espinoza et al. (2018) into a biophysical DG model to study its ability to generate synchronous activity.
Here, we use a biophysical computational network model of the DG (Braganza et al., 2020; Müller-Komorowska et al., 2023) with connectivity data of the PV+ INs in the DG of mice (Espinoza et al., 2018) and find that an ING mechanism is unlikely to support synchronous activity. Our results predict that experimental inhibition of recurrent PV+ IN connectivity would have little to no effects on synchronous activity in the DG.
Methods
PV+ interneuron ring model
We used a biophysical parvalbumin-positive (PV+) interneuron (IN) model as well as the inhibitory PV+-PV+ synapse parameters that are described in Braganza et al. (2020). However, our previous model did not have continuous distance-dependent connection probabilities for chemical synapses and did not model gap junctions (GJ). To model connection probabilities we fit sigmoid functions to data in Espinoza et al. (2018) for chemical synapses as well as GJs (Fig. 1C). That data was acquired from mice of either sex. To model the distance between neurons, we distributed them evenly on a ring with a circumference of 4 mm and calculated intersomatic distance on the ring. The intersomatic distance matrix is then converted to a connection matrix by making a connection with the probability given in Figure 1C. Figure 1D shows example connectivity matrices. For our main simulations, we chose to model 120 PV+ INs, which is roughly the number of PV+ INs in a 300 μm coronal slice of rat hippocampus. The entire unilateral rat dentate gyrus (DG) contains about 4,000 PV+ INs (Huusko et al., 2015) and it is about 1 cm in length from septal to temporal tip. We chose 120 PV+ INs for most simulations because our DG model is the model of a dorsal DG slice and upscaling the model to accommodate 4,000 PV+ INs would have made it extremely slow to simulate. Self-connectivity was excluded. Random realizations of this connectivity are shown in Figure 1D. To model GJs we used the NEURON (Carnevale and Hines, 2006) mechanism from Saraga et al. (2006). The resistance parameter was set manually to result in a plausible coupling coefficient of about 0.01, measured from soma to soma (Fig. 1B). The GJ mechanisms were placed in the middle of a randomly chosen proximal dendritic segment, about 32 μm from the middle of the soma. To bring the network into an asynchronous state, each neuron was injected at the soma with a constant current randomly chosen from a normal distribution with Iμ = 300 pA and Iσ = 50 pA. Iμ was gradually increased in Figure 3 to increase the overall network activity. For simulations of PV+ INs with Ih we used a NEURON mechanisms first described in Saraga et al. (2003). We added this mechanism to the soma of our PV+ IN model and adjusted the channel density parameter to result in a sag amplitude similar to that of CA3 PV+ INs (Papp et al., 2013). Extended Data Figure 2-2A shows the properties of the PV+ IN with Ih.
Dense recurrent connectivity synchronizes the network but biologically plausible connectivity does not. A, Intrinsic properties of the biophysical PV+ IN model. Note that it does not have Ih, as seen in the hyperpolarizing current step. Extended Data Figure 1-2 shows the models response to a chirp current injection B, Gap junction coupling implementation. The current injected in the red neuron increases the voltage in the blue neuron. The coupling coefficient is measured at the steady-state soma-to-soma. Gap junction resistance was hand-tuned to lead to the coupling coefficient. The gap junction was placed in the middle of the proximal dendritic segment (32 μm from the center of the soma). Extended Data Figure 1-1 shows the transmission of a granule cell EPSP through an interneuron gap juction. C, The biologically plausible sigmoid connection probability functions for chemical synapses and gap junctions. These were designed to fit the data in Espinoza et al. (2018). D, Random realization from the sigmoid functions fit to the measured connection probabilities from Espinoza et al. (2018) and the distances of the ring network (see methods). E, Spike raster plots of 120 PV+ INs with biologically plausible connectivity on the left in blue (D shows a biologically plausible connectivity example) and dense (nearly full) chemical connectivity on the right in orange. Nsyn is the number of chemical synapses received by the average neuron in this example run. F, The convolved output is calculated from the spike rasters above by convolving with a synaptic kernel. Shows the partially synchronous state for the full connectivity on the right. G, shows the coherence measure k(τ) for different values of τ (see methods) for the biologically plausible on the left (blue) and the dense connectivity on the right (orange). The dashed black line shows the line between (0, 0) and
Figure 1-1
Granule cell EPSP is measurable in a coupled PV+ IN. Gap junction resistance and placement as in Figure 1 B. Download Figure 1-1, TIF file.
Figure 1-2
The PV+ IN model does not resonate below threshold. The bottom shows a chirp current that is injected into the neuron. The frequency of the chirp current rises linearly from 1 Hz to 100 Hz over 5 s. The top shows the voltage response to the current as measured at the soma. On the left, the current was injected into the soma. On the right, the current was injected into a randomly chosen proximal dendritic segment. Download Figure 1-2, TIF file.
Figure 1-3
Doubling the strength of the PV-PV synapse does not induce synchrony with biologically plausible connectivity. A shows the strenght of the PV-PV synapse as it is used throughout the paper. The postsynaptic neuron is voltage clamped to 0 mV. The peak amplitude of the IPSC is 229 pA. That is larger than most IPSCs measured by Espinoza et al. (2018). B shows the spiking activity in the PV+ ring network with doubled synaptic conductance (15.2 nS) and biologically plausible connectivity. Download Figure 1-3, TIF file.
Figure 1-4
Keeping the number of synaptic inputs homogeneous does not change measured synchrony. In these simulations, each postynaptic neuron randomly chooses Nsyn presynaptic neurons to receive input from. Therefore, each neuron receives the same number of synapses. A,B&C show the same information as in main Figure 1, but for the simulations with equal synaptic inputs. D,E&F show the synchrony measures. In black are the same data as shown in main Figure 1. In green are the simulations with equal synapse numbers. Download Figure 1-4, TIF file.
Figure 1-5
Identical neurons with identical number of incoming synapses are exhibit perfectly synchronous activity. Each neuron receives input from 8 randomly chosen other neurons in the ring network. Each neuron receives a constant somatic current injection of 300 pA. Under these conditions all neurons fire at exactly the same time, there their correlation coefficient is 1, their
Figure 1-6
Varying the reversal potential of the inhibitory synapses in the ring model does not affect synchrony measures. For reference, the equilibrium potential in the baseline model is -70 mV. Download Figure 1-6, TIF file.
Figure 1-7
A network with low variance input current is initially synchronous but desynchronizes over time. For reference, the sigma of the input current in the baseline model is 0.3 nA. Download Figure 1-7, TIF file.
Mechanisms and parameters of the dentate gyrus model
The DG model is based on a model by Santhakumar et al. (2005) and we use largely the same ion channels but replace their double exponential synapse with Tsodyks-Markram type synapses (see below). The DG model used the following ion channel models. All neuron models used a Borg-Graham generic A-type potassium channel (borgka), a voltage dependent calcium activated potassium channel (cagk), a non-voltage-dependent calcium-activated potassium channel (gskch), a mechanism combining Hodgkin-Huxley style sodium and potassium conductances (ichan2) and an L-type calcium channel (lca). Furthermore, all neurons modeled intracellular calcium accumulation with exponential decay to baseline (ccanl). An N-type calcium channel (nca) was used in all cell types except for hilar-perforant path associated (HIPP) cells. T-Type calcium channels (cat) were used only in granule cells (GCs). Persistently modified h-channels (hyperde3) were used in HIPP and mossy cells (MCs). All conductances in all dendritic compartments and other intrinsic parameters are detailed in Extended Data Tables 4-1–4-5.
Chemical synaptic connections were all modeled with Tsodyks-Markram synapses (Tsodyks et al., 1998). The input spike trains are modeled as Poisson processes. The input innvervates GCs and PV+ cells. GCs excite MCs, PV+ cells and HIPP cells. MCs excite MCs, PV+ cells and HIPP cells. PV+ cells inhibit GCs, MCs and other PV+ cells. HIPP cells inhibit GCs, MCs and PV+ cells. The PV+ to PV+ connection is modeled with the distance dependent connectivity described in detail above. The other connections have a target pool of nearest postsynaptic cells. The size of the target pool differs from connection to connection. Within the target pool, the cells are chosen with uniform probability. Synaptic parameters are detailed in Extended Data Table 4-6.
Large scale PV+ IN plane model
To check whether the effects of IN connectivity on synchrony generation differ depending on the size of the network, we simulated 4,000 PV+ INs, since that represents roughly their total number in an entire unilateral DG. We also change from a ring network topology to a plane. The 2D plane was 1 cm in length and 0.4 cm in width. Cells were distributed uniformly on it. Distance-dependent connectivity was maintained as above through sigmoid functions. The results are shown in Extended Data Figure 2-1.
The PV+ IN ring network in the larger DG model
We integrated the PV+ IN ring model with distance-dependent chemical and GJ connectivity into our larger DG model (Braganza et al., 2020; Müller-Komorowska et al., 2023). To accommodate 120 PV+ INs we scaled up the model by a factor of five to 10,000 GCs, 300 MCs, 120 PV+ INs, and 120 HIPP cells. We removed the input current that generated asynchrony in the above versions of the ring model and instead added Poisson input processes that model input from 120 entorhinal cortex cells. For the homogeneous Poisson input (Fig. 4), the average input rate was 15 Hz. For the gamma oscillation conditions (Fig. 5), an inhomogeneous Poisson process was used with a sinusoidal waveform at 30 Hz (slow gamma) or 80 Hz (fast gamma). For Figure 6, we additionally recorded the somatic membrane voltage of all GCs.
Network synchrony measures
We mainly use the k(τ) measure (Wang and Buzsáki, 1996) to quantify the overall network synchrony. k(τ) is calculated as the cross correlation between two binned spike trains. The spike trains can be binned at different bin sizes τ. Given the total length of the simulation T and bin size τ we have
Finally, to quantify the power at different frequencies (Fig. 6) we calculated the mean membrane voltage of all GCs. We then estimated the power spectral density (PSD). For statistical testing, the PSD at theta (4 Hz–12 Hz) and gamma (30 Hz–100 Hz) was summed up to get a single number for power in each frequency band.
Software and simulation
Simulations were done with NEURON 8.0.0 (Carnevale and Hines, 2006) in Python 3.9.12 (Hines et al., 2009). Plotting was done with matplotlib 3.5.3 (Hunter, 2007) and seaborn 0.11.2 (Waskom, 2021). Tabular data was organized with pandas 1.4.3 (Reback et al., 2022). The statistical test to assess the effect of chemical synapses and GJs was the two-way ANOVA as implemented in the statsmodels 0.13.2 Python package (Seabold and Perktold, 2010). To calculate the PSD, we used the periodogram function from the scipy signal toolbox (Virtanen et al., 2020). p values were considered significant if p < 0.05. p and F values were rounded to the fourth significant decimal. For each condition, ten different networks were sampled using different random seeds. The Δt of all simulations was 0.0001 s. The duration of all simulations was 2 s and the entire duration was used in the analysis.
Results
Synchrony emerges from dense recurrent inhibition
To investigate the ability of biologically plausible connectivity to generate synchrony we started with a ring network consisting solely of parvalbumin-positive (PV+) interneurons (INs). Without connectivity, the network is asynchronous because we inject a steady current that is chosen from a normal distribution into the soma of each neuron. We used the biophysical model of a dentate gyrus (DG) PV+ IN that we used previously as part of a DG network model (Braganza et al., 2020; Müller-Komorowska et al., 2023). This model exhibits continuous firing in response to somatic step current injections and does not contain Ih channels as seen from the lack of slow depolarization in the hyperpolarizing current step (Fig. 1A. See also Extended Data Fig. 1-2, which shows no subthreshold resonance). To model the biologically plausible connectivity, we used paired patch clamp data from Espinoza et al. (2018). For gap junctions, they measured no distance dependence of the coupling coefficient. Therefore, we gave all gap junctions the same resistance to result in a coupling coefficient of about 0.1 measured from soma to soma (Fig. 1B). For lack of DG specific anatomical data, gap junctions were placed in the middle of the proximal dendritic segment, about 32.5 μm from the center of the somatic segment. We also ensured that a granule cell EPSP (as modeled in Braganza et al., 2020) is measurable at the soma of a coupled PV+ IN (Extended Data Fig. 1-1), a phenomenon measured by Espinoza et al. (2018).
To model the connection probability we hand-fit sigmoid functions to resemble the somatic distance-connection probability functions reported in Espinoza et al. (2018) for both gap junctions and chemical synapses (Fig. 1C). To model the overall number of PV+ INs and their somatic distance, we assumed that a unilateral DG contains about 4,000 PV+ INs (Huusko et al., 2015). Furthermore, to represent a 300 μm coronal slice we simulate 120 PV+ INs assuming that the DG is about 1 cm from temporal to medial tip in the rat. To simplify anatomy, we assumed a ring where neurons are distributed evenly, resulting in 33 μm between neurons. Figure 1D shows an example of a random realization of the sigmoidal connection probabilities and the neuronal distance. We took chemical synapses’ strength and decay time constant from our previous model implementation.
With this ring model in place, we moved to investigate synchrony in the spiking output of the network. To induce asynchronous activity, each neuron receives a current injection randomly drawn from a normal distribution. We found that for the biologically plausible connectivity, where the average PV+ IN receives 6.8 (± 2.31 SD) chemical synapses, the network remains asynchronous. However, when the sigmoid function is adjusted for nearly full chemical connectivity (118.9 synaptic inputs per neuron), the network is synchronous (Fig. 1E,F). To quantify synchrony, we calculate the measure k(τ) on the binned spike trains (Fig. 1G, see methods). Furthermore, we also calculate the average correlation coefficient of each neuronal pair in the network (Fig. 1H). All synchrony measures show that there is a discontinuity in network state at about 60–70 synapses per neuron. At sparser connectivity, synchrony increases only slightly with increasing synapse number. At about 60–70 synapses, synchrony changes slope and then increases linearly. This suggests that at dense connectivity with about 60–70 synapses, the network enters a partially synchronous state. At sparse biological connectivity (6.8 ± 2.31 SD synapses) on the other hand network activity is asynchronous.
Because in a complex neuronal network model synchrony depends not only on the average number of synapses, we performed various controls while holding the recurrent connectivity at the biologically plausible level. The synaptic conductance in our model induces a 229 pA peak current, which is on the higher end of what has been measured in vitro. But to make sure that synaptic strength does not preclude synchrony, we doubled the synaptic conductance (from 7.6 nS to 15.2 nS) and found no synchrony with biological connectivity (Extended Data Fig. 1-3). Besides, the connection strength, the synaptic equilibrium potential also affects synchrony (Jeong and Gutkin, 2007). We therefore changed the equilibrium potential but in our model we found no systematic effect of equilibrium potential on synchrony (Extended Data Fig. 1-6). For inhibition to induce synchrony, it has been shown that the synpatic rise-time should be at least as long as the action potential duration (Van et al., 1994). While in our model, the conductance rises instantaneously, the voltage change is filtered by the membrane capacitance, and therefore, the IPSP rise-time is indeed longer than the action potential. We also tested the influence of the variance of the injected current that desynchronizes neurons in our model. We found that extremely low variance has some effect on inital synchrony but even with variance as low as 0.0001 nA, the network becomes asynchronous after about 200 ms (Extended Data Fig. 1-7). Only in a network where each neuron receives identical current and has the same number of incoming synapses all neurons remain synchronous indefinitely (Extended Data Fig. 1-5). But as expected, keeping the number of incoming synapses identical does not induce synchrony, because of the variable current (Extended Data Fig. 1-4). Next, we selectively turned the biologically plausible gap junction and chemical connectivity on, to investigate their effects on synchrony separately.
Biologically plausible recurrent connectivity is insufficient for synchrony
Neither gap junctions, nor chemical synapses nor their interaction, are sufficient to induce the partially synchronous state (Fig. 2A). However, both have significant effects on the average firing rate of the network (Fig. 2B). Statistical analysis of the three synchrony measures (Fig. 2C–E) showed a significant increase only of
Adding PV+ IN connectivity to the ring network does not consistently increase network synchrony. A Spike raster plots (top) and convolved output (bottom) for three different conditions of the PV+ ring network. GJ-/Sy- has neither gap junctions nor chemical synapses. GJ+/Sy- has only gap junctions, GJ-/Sy+ has only chemical synapses and GJ+/Sy+ has gap junctions and chemical synapses. GJ+/Sy+ is a different random realization of the connectivity shown in Figure 1D. The convolved output does not show a partially synchronous state. B, Average frequency including all neurons of the network. Both gap junctions and chemical synapses significantly decrease the average activity of the ring network with a much larger effect of chemical synapses. Two-way ANOVA: Interaction: F = 0.0133, p = 0.9084, Main effects: GJ, F = 8.7578, p < 0.01; Sy, F = 5206.6925, p < 0.001. C–E, The three synchrony measures. Two-way ANOVA only showed significant effects for
Figure 2-1
Large PV + IN plane network. 4000 neurons were randomly distributed on a rectangular plane. A The spike raster plots show only every 40th neuron. Raster plot and convolved output show that the state is asynchronous. B Average frequency of all neurons in the network. Two-way ANOVA: Interaction: F = 1.893, p = 0.1743, Main effects: GJ, F = 14.0517, p < 0.001; Sy, F = 164105.4068, p < 0.001. C-D The synchrony measures. Two-way ANOVA for
Figure 2-2
Ring network simulations with PV+ INs containing Ih current. A Voltage response of the modified PV+ IN with Ih added. Top: voltage response to depolarizing current step. Middle: response to hyperpolarizing current step with the slow depolarization characteristic for Ih. Bottom: The current injections that induce the voltage responses above. B Spike raster plot and convolved output show an asynchronous state. C Average frequency of all neurons in the network. Twoway ANOVA: Interaction: F = 0.0926, p = 0.7626, Main effects: GJ, F = 1.8648, p = 0.1805; Sy, F = 5884.6182, p < 0.001. D-F The synchrony measures. Two-way ANOVA for
Figure 2-3
Ring network simulation with strong input current resulting in high average frequency. A Spike raster plot and convolved output show that some synchronous oscillations are larger when chemical synapses are added. However, this synchronous state is different from the one in Figure 1, where cells become silent during a cycle and furthermore does not reach the magnitude. B Average frequency of all neurons in the network. Two-way ANOVA: Interaction: F = 0.1127, p = 0.7391, Main effects: GJ, F = 5.3695, p < 0.05; Sy, F = 7474.6760, p < 0.001. C-E The synchrony measures. Two-way ANOVA for
Increasing the input drive increases network synchrony but does not reach the magnitude found with dense network connectivity. A, Spike raster plots for the PV+ ring network with biologically plausible connectivity for the standard model on the left (blue, input strength as in Figs. 1, 2), and a model with strong input drive on the right. The mean frequency fμ is calculated from all neurons in the network. B, The convolved output of the above spiking activity. C, The coherence measure k(τ) for different values of τ (see methods) for the baseline activity condition on the left (blue) and the high activity on the right (orange). The dashed black line shows the line between (0, 0) and
The PV+ IN ring network connectivity has no effect on GC synchrony and inconsistent effects on the PV+ ring network itself in the biophysical DG model. A,G, Schematic of the DG model. Green highlights the cells that are analyzed. B-F, GCs and H-L PV+ cells. Only the connectivity of the ring network was changed in the different conditions. B, Spike raster plots (top) and corresponding convolved output of EC neurons (black) and GCs. Note that the time constant of the kernel to calculate the convolved output was kept constant to the fast 1.8 ms of the PV+ neuron synapse for all cell types. C, Average frequency of the network calculated from all neurons. Two-way ANOVA showed no significant effects. D–F, The synchrony measures. Two-way ANOVA showed no significant effects. Note that in E values are large and negative. This is likely a failure of the supralinear k measure that occurs when the average frequency is slow. H, Same as B but with spike raster plots and convolved output for the PV+ INs instead of the GCs. I, Average frequency of the network calculated from all neurons. Two-way ANOVA: Interaction: F = 0.0881, p = 0.7683, Main effects: GJ, F = 1.8815, p = 0.1787; Sy, F = 30.9370, p < 0.001. J–L, The synchrony measures as in D-F but for the PV+ neurons. Two-way ANOVA showed no significant results in J. Two-way ANOVA for Supralinear K: Interaction: F = 0.9901, p = 0.3264, Main effects: GJ, F = 5.7611, p < 0.05; Sy, F = 0.0054, p = 0.9421. Two-way ANOVA for Correlation: Interaction: F = 0.2310, p = 0.6337, Main effects: GJ, F = 0.7404, p = 0.3952; Sy, F = 5.5301, p < 0.05. Extended Data Figure 4-1 shows the activity of HIPP cells and mossy cells. Extended Data Tables 4-1–4-5 contain the intrinsic parameters of the DG model. Extended Data Table 4-6 contains the synaptic parameters of the DG model.
Figure 4-1
Results for the HC (top) and MC (bottom), which were simulated for Figure 4 but not shown. Two-way ANOVA was done for all of the box plots but showed no significant result for any of them. Download Figure 4-1, TIF file.
Table 4-1
Parameters that are equal in all four neuron types. ccanl, intracellular calcium accumulation with exponential decay to baseline. Download Table 4-1, DOCX file.
Table 4-2
The intrinsic parameters of the granule cell model. borgka, Borg-Grahamgeneric A-type potassium channel; cagk, a voltage dependent calcium activated potassium channel; gskch, a nonvoltage-dependent calcium-activated potassium channel; ichan2, a mechanism combining Hodgkin-Huxley style sodium and potassium conductances; lca, an L-type calcium channel;ore, An N-type calcium channel (nca) was used in all cell types except for HIPP cells. T-Type calcium channels (cat) were used only in granule cells. Persistently modified h-channels (hyperde3) were used in HIPP and Mossy cells. Download Table 4-2, DOCX file.
Table 4-3
The intrinsic parameters of the basket cell model. Download Table 4-3, DOCX file.
Table 4-4
The intrinsic parameters of the mossy cell model. Download Table 4-4, DOCX file.
Table 4-5
The intrinsic parameters of the HIPP cell model. Download Table 4-5, DOCX file.
Table 4-6
The parameters of the synaptic connections. Columns: Pre, the presynaptic population; Post, the postsynaptic populations; N Target, number of nearby neurons that a neuron can connect to; Dendrite, dendrite of the postsynaptic neuron that is targeted by the synapse; Divergence, number of postsynaptic neurons that each presynaptic neuronc chooses randomly from the population defined by N Target; τ1, Decay time constant; τfacil, facilitation time constant of the Tsodyks-Markam dynamics. Download Table 4-6, DOCX file.
Recurrent connectivity desynchronizes PV+ INs in the DG network model (same as in Fig. 4A,C) when the input is 30 Hz frequency modulated but not when the input is 80 Hz modulated. This figure shows results from the biophysical DG model. Results for the other cell types are shown in Extended Data Figure s 5-1 and 5-2 A, Spike raster plots (top) and corresponding convolved output of EC neurons (black) and GCs in the slow gamma (30 Hz) modulated condition. Note that the time constant of the kernel to calculate the convolved output was kept constant to the fast 1.8 ms of the PV+ neuron synapse for all cell types. B, Average frequency of the network calculated from all neurons. Two-way ANOVA: Interaction: F = 0.0215, p = 0.8842, Main effects: GJ, F = 0.0256, p = 0.8737; Sy, F = 80.9939, p < 0.001. C–E, The synchrony measures for the PV+ neurons. Two-way ANOVA for
Figure 5-1
Cell types that were simulated in the slow gamma condition in Figure 5 but not shown. Boxplots were omitted because none of them was significant. Download Figure 5-1, TIF file.
Figure 5-2
Cell types that were simulated in the fast gamma condition in Figure 5 but not shown. Boxplots were omitted because none of them was significant with the sole exception of a significant increase of average HC frequency in the fast gamma condition when synaptic connectivity is added in the PV+ ring network: F = 5.322589, p < 0.05. Download Figure 5-2, TIF file.
Inhibitory recurrence does not change theta or gamma power measured through the GC membrane voltage. A, Single run examples of the average membrane potential across all granule cells. The color labeling for the conditions is between A and B and applies to both panels A & B. B, The power spectral density of the average membrane potential averaged across ten different simulation runs per condition. C, Boxplots of the summed PSD in the theta range. Two-way ANOVA showed no significant main effects or interactions for any of the conditions. D, Boxplots of the summed PSD in the gamma range. Two-way ANOVA showed no significant main effects or interactions for any of the conditions.
So far we simulated networks of 120 PV+ INs that represent a small part of the entire DG PV+ population. To ensure that we do not miss the effects of network size, we simulated networks with 4,000 PV+ INs with somatic positions randomly distributed on a 2D plane. This also changes the even distribution of the cells in the ring model. In this large network, we again found no synchronous state, despite clear effects of both gap junctions and chemical synapses on the average rate (Extended Data Fig. 2-1).
Although no strong Ih dependent sag potential has been measured in DG PV+ INs (Lübke et al., 1998), could Ih be necessary in other areas to enable the synchronizing effects of connectivity? Ih could be important in our model, because the slow current is conducted well by gap junctions. Furthermore, Ih is important for for synchrony in the inferior olive (Bal and McCormick, 1997; Lüthi and McCormick, 1998). We therefore added Ih to the soma of the PV+ model. The result was the same: an asynchronous network. However, the significant effect of gap junctions on the average firing rate was abolished by Ih (Extended Data Fig. 2-2). Within the asynchronous regime, synaptic connections caused small but significant increases in supralinear k and the correlation measure. The input drive of the neurons was the same in the models thus far. Because overall activity can influence synchrony, we next simulated different input strengths.
To increase the overall network activity, we increased the mean of the distribution (
In the DG the GCs are the principal cell type and outnumber basket cells by a factor of 1:100 in the suprapyramidal and 1:180 in the infrapyramidal blade (Amaral et al., 2007). We therefore first investigated whether connectivity in the PV+ ring network affects GC activity (Fig. 4A). As expected, biologically plausible IN connectivity did not synchronize the GC population (Fig. 4B). Even the average frequency of GCs is not significantly affected by IN connectivity (Fig. 4C) and neither are any of the three synchrony measures (Fig. 4D–F). This is in contrast to the PV+ neurons (Fig. 4G). While they did not reach a partially synchronous state (Fig. 4H), their average frequency was significantly decreased by inhibitory synaptic connectivity (Fig. 4I). The correlation measure of PV+ neurons on the other hand was significantly decreased by chemical synapses (Fig. 4L). However, like the upward trends in the ring network, the effect was small in magnitude and importantly, does not affect GC synchrony. The data for mossy cells and HC cells is shown in Extended Data Fig. 4-1 but none of the parameters showed statistically significant effects for them. Overall, this data does not give any indication that PV+ neurons can generate synchrony given homogeneous Poisson input from entorhinal cortex (EC) neurons. But biologically plausible input is unlikely to be homogeneous. It is rather expected to have some oscillatory structure. We therefore went on to simulate synaptic input that is modulated at slow gamma (30 Hz) or high gamma (80 Hz) frequency.
Because GCs showed no significant effects with the homogeneous Poisson input, we focused this analysis on the PV+ INs. For slow gamma (30 Hz) modulated input the PV+ INs did not enter a partially synchronous state (Fig. 5A), although their average frequency was significantly affected by their connectivity (Fig. 5B). The synchrony measures all showed small but significant decreases in synchrony (Fig. 5C–E). No significant effects of gap junctions or the interaction were found. For the fast gamma (80 Hz) modulated inputs (Fig. 5F), the result was overall similar, however
So far, we have used generic synchrony measures of the output spike train. During experiments in biological tissue, synchrony is measured through field potential and evaluated at specific frequencies. We therefore finished our analysis by calculating the power spectrum of the average granule cell membrane potentials, because they broadly integrate direct as well as recurrent input currents. We ran simulations for the DG model and the different conditions as shown in the above figures, but recorded the somatic membrane voltage of all 10,000 GCs. Figure 6A shows the average membrane potential of all GCs of a representative model run for the homogeneous Poisson input, slow gamma (30 Hz) modulated and fast gamma (80 Hz) modulated. Regardless of the input modulation, gap junctions and synapses do not have a clear effect on the oscillatory frequencies. This is also true when looking at the power spectral density (PSD) in Figure 6B. This shows the average PSD over 10 model runs in each condition. For statistical testing, we calculated the summed PSD at theta (4 Hz-12 Hz) and gamma (30 Hz-100 Hz). No statistically significant main effects or interactions for any of the conditions were found. This highlights our conclusion that IN recurrence alone does not synchronize the DG network.
Discussion
We have performed a variety of simulations, all intending to find synchronous neuronal activity and all with the same result: biologically plausible connectivity does not generate synchronous neuronal activity. We started by varying connectivity from sparse and biologically plausible to fully connected. As reported before (Wang and Buzsáki, 1996; Bartos et al., 2002), we found that a partially synchronous state is possible in the ring network if the number of recurrent synapses is large enough. This is consistent with the interneuron gamma generator (ING) model of synchrony. However, more synapses are required than we expect from distance-based connectivity probabilities (Espinoza et al., 2018). This is also true when we combine chemical and electrical synapses, which have been shown to have complementary roles in synchronizing neuronal firing (Kopell and Ermentrout, 2004). Comparing simulations with biologically plausible recurrent connectivity to those without such connectivity showed small but significant effects of connectivity on some synchrony measures. Effects were particularly clear when Ih current was added to the neuronal dynamics or the average frequency was high. However, the degree of synchrony did not reach the same level as we have observed in the partially synchronous state and was varied between the three synchrony measures. Furthermore, we observed neither synchrony nor increased theta or gamma power in the larger dentate gyrus (DG) model due to recurrent IN connectivity. On the contrary, PV+ INs showed very slight desynchronization for frequency-modulated inputs. Therefore, we conclude that in the DG, recurrent interneuron connectivity is too sparse to cause synchronous activity.
Previous papers most likely found synchrony in ring networks of recurrently connected INs (Wang and Buzsáki, 1996; Bartos et al., 2002) because they have assumed more synapses than we have in our model. Theoretical work has well established that in a ring network, synchrony requires a sufficiently large number of inhibitory synapses (Golomb and Hansel, 2000). Below that number, the network is asynchronous. Wang and Buzsáki (1996) also investigated the influence of synapse number and found that at least 60 synaptic inputs per neuron are required for synchrony. Our estimate of 70 synapses is not far off, although we did not attempt to replicate the exact neuron model, which is known to be relevant to synchronization (Nomura et al., 2003). Wang and Buzsáki (1996) consider 60 a sparse and plausible number based on an estimation of synaptic contacts made by a biocytin filled PV+ basket cell in CA1 (Sik et al., 1995) and anatomical data regarding PV+ neuron density in CA1 (Aika et al., 1994). Bartos et al. (2002) extrapolate these CA1 estimates to the DG and thereby likely overstate the number of recurrent DG synapses. The primary difference between ours and previous studies is that we use DG connectivity estimates (Espinoza et al., 2018) instead of CA1 estimates.
The circuit mechanisms of synchrony can also be studied experimentally. Unfortunately, experiments that inhibit PV+ INs are inadequate for conclusions about recurrent inhibition as they perturb the influence of all inputs, not only inhibitory recurrence. A study by Wulff et al. (2009) has come close to directly testing the effect of inhibitory recurrence by deleting the GABAA γ2 subunit in PV+ INs. Measuring in CA1, they found an effect on gamma oscillations. We predict the same result for the DG. However, Wulff et al. (2009) find effects on theta power, which our model does not predict for the DG. Moreover, the genetic deletion approach has two important limitations: it affects all PV+ neurons of the brain and affects mice throughout their development. Therefore, compensatory effects could obscure the role of recurrent inhibition in gamma generation. An opto-/chemogenetic approach targeting recurrent synapses specifically could overcome these limitations, but it is currently infeasible. Alternatively, an inducible genetic deletion could achieve the goal. However, the genetic deletion results from Wulff et al. (2009) support our prediction that recurrent inhibition has no major role in gamma generation.
Besides recurrent inhibition, feedback between pyramidal cells and interneurons can also generate synchrony. This is called the pyramidal-interneuron gamma generator (PING; Tiesinga and Sejnowski, 2009) model. Since hippocampal networks feature IN to IN and PC to IN connectivity, both mechanisms are assumed to cooperate in gamma generation (Buzsáki and Wang, 2012). When we turn recurrent connections between INs off, we find no effect on theta or gamma power, suggesting that the two mechanisms do not meaningfully interact. However, while our model has granule cell (GC) to IN feedback connectivity, it does not show gamma power. Therefore, our model does not act as a gamma generator and may not be a good representative of a PING model. For example, our GC model does not exhibit bursting, a known GC dynamic in-vitro and in-vivo (Pan and Stringer, 1996; Pernia-Andrade and Jonas, 2014; Neubrandt et al., 2018). GC bursting could be a necessary component of a DG PING model. In its current implementation, however, recurrent inhibition affects neither theta nor gamma power.
Although our results cannot exclude the possible existence of a parameter region in which recurrent inhibitory connectivity is sufficient to synchronize DG granule cells, they indicate that such a region, if any exists, is likely to be small. We therefore predict that experimental inhibition of recurrent PV+ IN connectivity would have little to no effects on synchronous activity in the DG.
Footnotes
The authors have no conflicts of interest to declare.
We thank Joanna Komorowska-Müller for her helpful comments on the manuscript. This work was partly supported by JSPS KAKENHI no. JP23H05476.
This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license, which permits unrestricted use, distribution and reproduction in any medium provided that the original work is properly attributed.
References
Data and Code Availability
The code will be available as a release in the pydentate repository. The data is available on Zenodo: https://zenodo.org/records/15186260.