Skip to main content

Main menu

  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Current Issue
    • Issue Archive
    • Blog
    • Collections
    • Podcast
  • TOPICS
    • Cognition and Behavior
    • Development
    • Disorders of the Nervous System
    • History, Teaching and Public Awareness
    • Integrative Systems
    • Neuronal Excitability
    • Novel Tools and Methods
    • Sensory and Motor Systems
  • ALERTS
  • FOR AUTHORS
  • ABOUT
    • Overview
    • Editorial Board
    • For the Media
    • Privacy Policy
    • Contact Us
    • Feedback
  • SUBMIT

User menu

Search

  • Advanced search
eNeuro

eNeuro

Advanced Search

 

  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Current Issue
    • Issue Archive
    • Blog
    • Collections
    • Podcast
  • TOPICS
    • Cognition and Behavior
    • Development
    • Disorders of the Nervous System
    • History, Teaching and Public Awareness
    • Integrative Systems
    • Neuronal Excitability
    • Novel Tools and Methods
    • Sensory and Motor Systems
  • ALERTS
  • FOR AUTHORS
  • ABOUT
    • Overview
    • Editorial Board
    • For the Media
    • Privacy Policy
    • Contact Us
    • Feedback
  • SUBMIT
PreviousNext
Research ArticleNew Research, Neuronal Excitability

Gamma Oscillations in the Basolateral Amygdala: Biophysical Mechanisms and Computational Consequences

Feng Feng, Drew B. Headley, Alon Amir, Vasiliki Kanta, Ziao Chen, Denis Paré and Satish S. Nair
eNeuro 11 January 2019, 6 (1) ENEURO.0388-18.2018; DOI: https://doi.org/10.1523/ENEURO.0388-18.2018
Feng Feng
1Department of Electrical Engineering and Computer Science, and Department of Bioengineering, University of Missouri, Columbia, MO 65211
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Drew B. Headley
2Center for Molecular and Behavioral Neuroscience, Rutgers University – Newark, Newark, NJ 07102
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Alon Amir
2Center for Molecular and Behavioral Neuroscience, Rutgers University – Newark, Newark, NJ 07102
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Vasiliki Kanta
2Center for Molecular and Behavioral Neuroscience, Rutgers University – Newark, Newark, NJ 07102
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Ziao Chen
1Department of Electrical Engineering and Computer Science, and Department of Bioengineering, University of Missouri, Columbia, MO 65211
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Denis Paré
2Center for Molecular and Behavioral Neuroscience, Rutgers University – Newark, Newark, NJ 07102
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Satish S. Nair
1Department of Electrical Engineering and Computer Science, and Department of Bioengineering, University of Missouri, Columbia, MO 65211
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • Article
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF
Loading

Abstract

The basolateral nucleus of the amygdala (BL) is thought to support numerous emotional behaviors through specific microcircuits. These are often thought to be comprised of feedforward networks of principal cells (PNs) and interneurons. Neither well-understood nor often considered are recurrent and feedback connections, which likely engender oscillatory dynamics within BL. Indeed, oscillations in the gamma frequency range (40 − 100 Hz) are known to occur in the BL, and yet their origin and effect on local circuits remains unknown. To address this, we constructed a biophysically and anatomically detailed model of the rat BL and its local field potential (LFP) based on the physiological and anatomical literature, along with in vivo and in vitro data we collected on the activities of neurons within the rat BL. Remarkably, the model produced intermittent gamma oscillations (∼50 − 70 Hz) whose properties matched those recorded in vivo, including their entrainment of spiking. BL gamma-band oscillations were generated by the intrinsic circuitry, depending upon reciprocal interactions between PNs and fast-spiking interneurons (FSIs), while connections within these cell types affected the rhythm’s frequency. The model allowed us to conduct experimentally impossible tests to characterize the synaptic and spatial properties of gamma. The entrainment of individual neurons to gamma depended on the number of afferent connections they received, and gamma bursts were spatially restricted in the BL. Importantly, the gamma rhythm synchronized PNs and mediated competition between ensembles. Together, these results indicate that the recurrent connectivity of BL expands its computational and communication repertoire.

  • amygdala
  • gamma oscillations
  • computational model
  • biophysical model
  • extracellular potential

Significance Statement

Using in vitro and in vivo data we develop the first large-scale biophysically and anatomically realistic model of the basolateral amygdala nucleus (BL), which reproduces the dynamics of the in vivo local field potential (LFP). Significantly, it predicts that BL intrinsically generates the transient gamma oscillations observed in vivo. The model permitted exploration of the poorly understood synaptic mechanisms underlying gamma genesis in BL, and the model's ability to compute LFPs at arbitrary numbers of recording sites provided insights into the characteristics of the spatial properties of gamma bursts. Furthermore, we show how gamma synchronizes principal cells (PNs) to overcome their low firing rates while simultaneously promoting competition, potentially impacting their afferent selectivity and efferent drive, and thus emotional behavior.

Introduction

The basolateral complex of the amygdala (BLA) supports emotional learning and expression (Ehrlich et al., 2009; Duvarci and Pare, 2014). Central to this function is integrating information from numerous cortical and thalamic regions (Turner and Herkenham, 1991; McDonald, 1998). It is commonly thought that these signals converge on subpopulations of principal cells (PNs), which in turn either synapse on other PNs that promote a particular behavior, or contact interneurons that inhibit PNs supporting a countervailing behavior (LeDoux, 2000; Ehrlich et al., 2009; Duvarci and Pare, 2014). However, recurrent connections also exist within the BLA (Paré et al., 1995; Samson and Paré, 2006; Beyeler et al., 2016), with fast-spiking interneurons (FSIs) synapsing on each other and forming reciprocal connections with local PNs (Smith et al., 1998; Woodruff and Sah, 2007). Yet, the importance of these connections for BLA function remains largely unknown.

One possibility is that recurrent connections support the generation of oscillations. Consistent with this possibility, computational models (Traub et al., 1996; Wang and Buzsáki, 1996) as well as in vitro (Traub et al., 1996; Sohal et al., 2009) and in vivo experiments (Penttonen et al., 1998; Cardin et al., 2009) have revealed that a dense recurrent network of PNs and FSIs produces oscillations in the gamma frequency band. The dominant model for this, known as the pyramidal-interneuron network gamma (PING) model (Whittington et al., 2000), posits that the firing of PNs excites FSIs, which in turn deliver feedback inhibition, transiently silencing PNs. As the inhibition wanes, PNs regain the ability to fire and can restart the gamma cycle.

Crucially, during affective experiences, the BLA also exhibits gamma oscillations that are especially pronounced in its basolateral nucleus (BL; Bauer et al., 2007). For instance, gamma increases when rodents regulate their anxiety level during open field exploration (Stujenske et al., 2014) or are exposed to emotionally charged stimuli (Bauer et al., 2007). Importantly, the human amygdala also produces gamma oscillations during emotionally arousing stimuli (Oya et al., 2002). Despite the prevalence of gamma oscillations in the amygdala, their cellular basis and function remain unclear.

Numerous functions have been ascribed to gamma oscillations (Wang, 2010), but two stand out in particular. First, they synchronize spiking. PNs in networks exhibiting gamma oscillations tend to fire together more often than expected by chance (Wang and Buzsáki, 1996), robustly driving downstream neurons (Salinas and Sejnowski, 2000, 2002; Zandvakili and Kohn, 2015). Second, they may mediate competitive interactions between PN ensembles (Börgers et al., 2008). During each gamma cycle, the ensemble with the strongest afferent drive will tend to recruit the local FSI network, suppressing weakly driven ensembles (de Almeida et al., 2009).

Computational models of gamma oscillations have not been reported for the amygdala; such models for other brain regions have typically used generic single cell and network configurations, with surrogate local field potential (LFP) models (e.g., Börgers et al., 2008; Palmigiano et al., 2017). To examine whether the intrinsic BL circuitry can produce the poorly-understood transient gamma oscillations and their associated functions, we created a large-scale 27,000 cell multi-compartmental biophysical model of BL, with a detailed LFP model, that recapitulates numerous features of BL activity in vivo, such as the properties of gamma oscillations in LFPs and their entrainment of spiking. Using this model, we determined the circuit elements essential for generating gamma oscillations in BL, what aspects of microcircuit architecture affect the participation of neurons in gamma, the spatiotemporal properties of BL gamma bursts, and how these oscillations support information processing.

Materials and Methods

Experimental data

To construct and validate our model, we conducted new analyses on two sets of extracellular recordings that were described in prior publications (Headley et al., 2015; Amir et al., 2018), and on one unpublished dataset. Thus, we briefly describe the methods used in these prior studies so that readers can assess the nature and quality of these data.

In vivo chronic recordings

All animal procedures were approved by the Institutional Animal Care and Use Committee at Rutgers University, in accordance with the Guide for the Care and Use of Laboratory Animals (Department of Health and Human Services). Unit activity from prefrontal (PFC), perirhinal (PR), and entorhinal (ER) cortices (Headley et al., 2015) was used to generate surrogate spike trains that simulate extrinsic afferents onto the BL model (see below, BL afferents). These data were acquired in three male Long–Evans rats weighing between 350 and 500 g that were implanted with a headcap containing microdrives loaded with tetrodes (20-µm tungsten wire, impedance <100 kΩ). Two independent drives were used to target either prefrontal or PR/ER (PFC: AP +3.0, ML +0.5, DV 3.0; PR: AP −3.0 to −8.0, ML +6.0 to +7.2, DV 5.0; ER: AP −5.4 to −8.0, ML +7.0, DV 5.5; all coordinates in mm, DV was taken with respect to the pial surface). Following recovery from surgery (>7 d), microdrives were advanced until tetrodes reached their target locations, at which point recordings began. Extracellular signals were amplified with a 96-channel system (Plexon) and digitized (National Instruments) for offline analysis. Unit activity was sorted into single units by high pass filtering wideband LFP with a moving median filter, detection of spikes with amplitude >2 SD, automatic clustering of waveforms in principal component space (KlustaKwik), and manual refinement of cluster assignment (Klusters). Validation of single unit quality and isolation can be found in our prior paper (Headley et al., 2015). Only regular spiking units (putative projection neurons), classified using k-means clustering of the negative peak to positive peak time interval of their waveform and firing rate, were used. In particular, we focused on a single epoch from each subject (subject 1: 165 s; subject 2: 140 s; subject 3: 160 s) with simultaneous recordings from multiple single units (subject 1: 36; subject 2: 14; subject 3: 20).

Unit and LFP activities from the BL recorded with silicon probes were used to calibrate and validate our model (n = 5). As detailed in Amir et al. (2018), male Sprague Dawley rats (>350 g) were implanted with either a 32 or 64 channel silicon probe (Neuronexus) mounted on a microdrive and targeting the BL nucleus of the amygdala (AP −2.2 to −3.6, ML +5 − 5.3, and DV 8.8). Recording hardware and single unit sorting was similar to that described above.

LFP spectra were also recorded from the BL of male Long–Evans rats (>350 g, n = 7). A microwire (20-µm tungsten wire, impedance <50 kΩ) attached to a fiber optic stub was stereotaxically implanted in the BLA (AP −2.5, ML +5.0, DV −7.5, from brain surface) and fixed to the skull with dental cement (Metabond and Teets Cold Cure). These subjects were also used for optogenetic experiments (AAV5-Syn-Chronos-GFP injected into BL), but these data were excluded from the present study. All recordings were performed at least 24 h following any optogenetic manipulations. Extracellular recordings from the BL were obtained with a 32-channel digitizing headstage (Intan Technologies).

All the electrophysiological data used in the present study was referenced to a screw fixed to the bone overlaying the cerebellum and was obtained while rats were allowed to behave spontaneously in a neutral plastic enclosure. We only considered data acquired in the quiet waking state (QW), which was identified by the absence of gross body movement and LFPs with relatively low power at frequencies <4 Hz.

In vitro BL slice recordings

An overdose of isoflurane was administered to deeply anesthetize male Long–Evans rats (n = 11, 53 cells). Once all reflexes had ceased, they were immediately perfused through the heart with a ice cold modified artificial CSF (aCSF) solution containing: 103 mM N-methyl-D-gluconamine, 2.5 mM KCl, 1.2 mM NaH2PO4, 30 mM NaHCO3, 10 mM MgSO4 × 7H2O, 25 mM glucose, 20 mM HEPES, 101 mM HCl, 2 mM Thiourea, 3 mM Na-Pyruvate, 12 mM N-acetyl-L-cysteine, and 0.5 mM CaCl. Following 30 s of perfusion, the brain was rapidly extracted and placed in the cutting solution, which was the same used for perfusion. Slices were cut in the coronal plane on a vibratome with a thickness of 300 − 400 µm. Following slicing sections were transferred to a chamber containing the perfusion solution at 32° C for 5 min. Slices were then placed in a holding chamber with room temperature (22°C) normal aCSF: 124 mM NaCl, 2.5 mM KCl, 1.25 mM NaH2PO4, 26 mM NaHCO3, 1 mM MgCl2, 2 mM CaCl2, and 10 mM glucose (pH 7.2 − 7.3, 305 mOsm). Slices were then transferred one at a time to a recording chamber perfused with oxygenated aCSF at 32°C.

Under visual guidance by infrared video microscopy, we obtained whole-cell patch-clamp recordings from BL neurons. We pulled patch pipettes of borosilicate glass with tip resistances of 5 − 8 MΩ and filled them with a solution of 130 mM K-gluconate, 10 mM N-2-hydroxyethylpiperazine-N’-2’-ethanesulfonic acid, 10 mM KCl, 2 mM MgCl2, 2 mM ATP-Mg, and 0.2 mM GTP-tris(hydroxy-methyl)aminomethane (pH 7.2, 280 mOsm). We did not compensate for the liquid junction potential, which is 10 mV with this solution. Current clamp recordings were obtained with a MultiClamp 700B amplifier and digitized at 20 kHz using an Axon Digidata 1550 interface.

Once whole-cell access was achieved, we characterized the neuron’s electroresponsive properties. Current pulses 500 ms in length were delivered from −200 to 360 pA in steps of 40 pA. Neurons with action potential (AP) half-widths <0.35 ms were classified as FSIs, while all others were regular spiking.

Model implementation

The single cell and network models were developed using the parallel NEURON 7.4 simulator (Carnevale and Hines, 2005), and simulations were run with a fixed time step of 50 µs. Network models were run using the NeuroScience Gateway (NSG; www.nsgportal.org) that provides free and easy access to high-performance computers (Sivagnanam et al., 2013). Model results were obtained by averaging five network model runs with different random seeds for data shown in Figures 1G, 2B, 3C, 4, 5E and 6B1-B6 and for selecting groups of neurons were used in Figure 7B–D.

Figure 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 1.

Construction of a biophysically accurate BL model. A, Our BL model was a cube structure comprised of 27,000 PNs and FSIs randomly distributed throughout. Local cell densities in 100 × 100 × 100 µm cubes (red = PN, blue = FSI) are illustrated for a typical case. For this example, PNs, the densities spanned 0–19,000/mm3 (median = 7000/mm3). FSI densities ranged from 0 to 5000/mm3 (median = 1000/mm3). B, A slice through the model shows the inhomogeneity of cell densities for PNs and FSIs. C, The dendritic processes of the neurons contained in a particular voxel are illustrated. These were randomly oriented and spanned several hundred microns as also shown in the three orthogonal views. D, (1) PNs and FSIs were distributed in space so that they spanned a volume equal to the area of the BL in the rat. Virtual current clamp electrodes (ic) could be placed into any cell, and virtual extracellular electrodes (lfp) could be placed anywhere in the model volume. PNs are indicated by red triangles and FSIs are blue circles. (2) PNs and FSIs were connected among themselves and with each other. Extrinsic glutamatergic afferents fed onto both cell types. Neither FSIs nor PNs formed autapses. E, The relative proportions of the three cell types were determined by patching neurons in BL slices prepared from adult Long–Evans rats, the same age and strain used for our LFP recordings. The cutoff for determining which PNs exhibited adaption was set to 1.5, which was between the two peaks in the distribution of adaptation ratios. F, Example recordings from neurons in the slice receiving current injection were comparable to those of our model neurons. G, The firing rate distributions for neurons in our model overlapped with the mean rates reported from the BL in vivo in a previous study (Amir et al., 2018). H, Example spike waveforms recorded with a silicon probe in vivo. Red traces are from putative PNs, while blue are FSIs. The intensity of the color is scaled to the peak of the spike wave form. I, Left, For both a model PN and FSI we delivered a suprathreshold EPSC to the a-dend and recorded the extracellular AP at different distances along the long axis of the neuron, and either at 50 or 100 µm lateral. Near the cell body the field was negative, and it decayed rapidly with distance. Positive dendritic return currents were evident as well. The FSI extracellular spike wave form was both smaller and faster than the PN’s. Right, Subthreshold stimulation resulted in a much weaker extracellular wave form (note the scale bar) reflecting the EPSC (EP) that was negative going near the stimulated dendritic branch. J, (1) Directly overlaying the extracellular APs from both cell types illustrates that those arising from PNs were much slower than those from FSIs. Amplitudes were rescaled so both spike waveforms occupy the same vertical extent of the graph. (2) AP amplitudes were much stronger than EP amplitudes for both cell types. K, For extracellular APs recorded with silicon probes in vivo, we measured how their amplitudes decayed with distance. The drop in amplitude was fit by an exponential curve (black line). The gray region is the 95% confidence bounds. Measuring the decay in our model extracellular APs along the lateral axis, we found that it fit within the in vivo distribution.

Mathematical equations for voltage-dependent ionic currents

The equation for each compartment (soma or dendrite) followed the Hodgkin-Huxley formulation (Byrne and Roberts, 2004; Kim et al., 2013) in Equation 1,Embedded Image (1)where Embedded Image are the somatic/dendritic membrane potential (mV), Embedded Image and Embedded Image are the intrinsic and synaptic currents in the soma, Embedded Image is the electrode current applied to the soma, Embedded Image is the membrane capacitance, Embedded Image is the is the conductance of leak channel, and Embedded Image is the coupling conductance between the soma and the dendrite (similar term added for other dendrites connected to the soma). The intrinsic current Embedded Image , was modeled as Embedded Image , where Embedded Image is its maximal conductance, m its activation variable (with exponent p), h its inactivation variable (with exponent q), and Embedded Image its reversal potential (a similar equation is used for the synaptic current Embedded Image but without m and h). The kinetic equation for each of the gating variables x (m or h) takes the formEmbedded Image (2)where Embedded Image is the steady state gating voltage- and/or Ca2+-dependent gating variable and Embedded Image is the voltage- and/or Ca2+-dependent time constant. The equation for the dendrite follows the same format with ‘s’ and ‘d’ switching positions in Equation 1.

PN models

To reproduce the range of spike frequency adaptation exhibited by BL PNs (Rainnie et al., 1993), we modeled two types of PNs, one with adaptation (Type-A) and one with continuous spiking (Type-C), differing solely in the magnitude of their Ca2+-dependent K+ current, either 50 or 0.2 mS/cm2, respectively (Kim et al., 2013). We conducted our own study of the relative proportion of these types because our in vivo recordings, which our model was compared with, were from adult rats (>350 g) whereas previous slice work was conducted in younger subjects, and the ionic mechanisms underlying accommodation, such as spike after-hyperpolarizations (AHPs), change with age (Disterhoft and Oh, 2007). Positive current injection into PNs (n = 56) evoked either a train of APs with increasing interspike intervals (i.e., accommodation) or at a fixed rate (i.e., continuous). This property could be quantified by the adaptation ratio:Embedded Image (3)with larger values indicating greater adaptation. We set a cutoff of 1.5 to classify PNs adapting. Our previous PN models (Li et al., 2009; Kim et al., 2013; Feng et al., 2016) were revised to incorporate the low-threshold and high-threshold oscillations (LTOs and HTOs, respectively) reported in BLA PNs (Pape et al., 1998) and modeled recently by our group (Alturki et al., 2016). The model PNs had three compartments representing a soma (diameter, 24.75 µm; length, 25 µm), where GABAergic synapses were placed, an apical dendrite (a-dend; diameter, 3 µm and length, 270 µm) where glutamatergic synapses were located, and another dendrite (p-dend; diameter, 5 µm and length, 555 µm) to match passive properties. Values of specific membrane resistance, membrane capacity and cytoplasmic (axial) resistivity were, respectively, Rm = 55 kΩ-cm2, Cm = 2.4 µF/cm2, and Ra = 150 Ω-cm. Leakage reversal potential (EL) was set to −75 mV. The resulting Vrest was −70.3 mV, input resistance (RIN) was ∼140 MΩ, and time constant (τm) was ∼30 ms, all of which were within the ranges reported in previous physiologic studies (Washburn and Moises, 1992). All compartments had the following currents: leak (IL), voltage-gated persistent muscarinic (IM), high-voltage activated Ca2+ (ICa), spike-generating sodium (INa), potassium delayed rectifier (IDR), A-type potassium (IA; Li et al., 2009; Power et al., 2011) and hyperpolarization-activated nonspecific cation (Ih) current. In addition, the soma had a slow apamin-insensitive, voltage-independent AHP current (IsAHP; Power et al., 2011; Alturki et al., 2016). See Tables 1, 2 for current equations and densities.

View this table:
  • View inline
  • View popup
Table 1.

Gating parameters of ion channels in BL PN neurons

View this table:
  • View inline
  • View popup
Table 2.

Maximal conductance densities in model BL PN neurons

Interneuron models

BL also contains local GABAergic interneurons that exhibit various firing patterns, even among neurochemically-homogeneous subgroups (Pape and Pare, 2010; Spampanato et al., 2011). However, the most prevalent are the fast-spiking parvalbumin-positive type, which has been implicated in the genesis of cortical and BL gamma (Börgers et al., 2005; Oren et al., 2006; Atallah and Scanziani, 2009; Amir et al., 2018). Accordingly, we modeled only the fast-spiking type of interneurons (FSI). The FSI model was the same as in Kim et al., 2013, with two compartments, a soma (diameter, 15 µm and length, 15 µm) and a dendrite (diameter, 10 µm and length, 150 µm). Each compartment contained a fast Na+ (INa) and a delayed rectifier K+ (IDR) current. The FSI model reproduced the short spike duration (with the spike duration at half amplitude <1 ms) that characterizes FSIs. The passive membrane properties of FSI cells were as follows: Rm = 20 kΩ-cm2, Cm = 1.0 µF/cm2, Ra = 150 Ω-cm. The FSI model also reproduced the non-adapting repetitive firing behavior of fast-spiking cells, as observed experimentally (Rainnie et al., 1993; Woodruff and Sah, 2007).

Network size and cell type proportions

Estimates of the number of neurons in rat BL vary widely (Chareyron et al., 2011), so we settled on the mean across studies, which was ∼72,000. We developed a scaled down (1:2.7) model of this region with 27,000 neurons randomly distributed in a cuboid geometry (1.4 × 1.4 × 1.4 mm), ensuring an intersoma distance >25 µm. The model included 64% PNA (n = 17,280), 26% PNC (n = 7020), and 10% FSIs. These proportions were based on in vitro results collected for this study and agreed with estimates found in the literature (McDonald and Mascagni, 2001; Muller et al., 2006). The dendrites of all neurons in the network had random orientations (McDonald, 1984).

Mathematical equations for synaptic currents

All excitatory transmission was mediated by AMPA/NMDA receptors, and inhibitory transmission by GABAA receptors. The corresponding synaptic currents were modeled by dual exponential functions (Destexhe et al., 1994; Durstewitz et al., 2000), as shown in Equations 4–6,Embedded Image Embedded Image Embedded Image (4) Embedded Image Embedded Image Embedded Image (5) Embedded Image Embedded Image Embedded Image (6)where V is the membrane potential (mV) of the compartment (dendrite or soma) where the synapse is located and w is the synaptic weight for the synapse. The synaptic reversal potentials were EAMPA = ENMDA = 0 mV and EGABAa = −75 mV (Durstewitz et al., 2000; Martina et al., 2001). The voltage-dependent variable s(V) which implements the Mg2+ block was defined as: s(V) = [1 + 0.33 exp(−0.06 V)]−1 (Zador et al., 1990). The terms ONNMDA and ONAMPA are set to 1 if the corresponding receptor is open, else to 0. The synaptic current rise and decay time constants are determined by αTmax and β (Destexhe et al., 1994). Synaptic parameter values are listed in Table 3.

View this table:
  • View inline
  • View popup
Table 3.

Parameters related to synaptic connections

Short-term presynaptic plasticity

All model AMPA and GABA synapses also exhibited short term pre-synaptic plasticity (Kim et al., 2013). Short-term depression was modeled at FSI→PN and PN→FIS connections based on experimental findings of Woodruff and Sah (2007) in BL, while between PNs, it was modeled based on results from neocortex (Silberberg et al., 2004) due to the lack of such experimental data in BL. Short term plasticity was implemented as follows (Hummos et al., 2014): for facilitation, the factor F was calculated using the equation: Embedded Image and was constrained to be ≥1. After each stimulus, F was multiplied by a constant, f (≥1) representing the amount of facilitation per pre-synaptic AP, and updated as Embedded Image . Between stimuli, F recovered exponentially back toward 1. A similar scheme was used to calculate the factor D for depression: Embedded Image and D constrained to be ≤1. After each stimulus, D was multiplied by a constant d (≤1) representing the amount of depression per pre-synaptic AP, and updated as Embedded Image . Between stimuli, D recovered exponentially back toward 1. We modeled depression using two factors d1 and d2 with d1 being fast and d2 being slow subtypes, and Embedded Image . The parameters for modeling short-term plasticity are listed in Table 4. Our model did not have long-term synaptic plasticity.

View this table:
  • View inline
  • View popup
Table 4.

Parameters related to short-term presynaptic plasticity

Intrinsic connections

Except for FSI→FSI connections with both electrical and chemical synapses, all other connections were implemented as chemical synapses. Connection probabilities have been found to be distance-dependent for PN→PN contacts in BL, and we used 3%, 2%, 1%, and 0.5% probabilities for intersoma distances of <50, 50 − 100, 100 − 200, and 200 − 600 µm, respectively (Abatis et al., 2017). For connections involving interneurons, we used data from in vitro BL reports (Woodruff and Sah, 2007), with connections limited to pairs within ∼300 µm of each other. Probabilities in the model for unidirectional connections from FSI→PN were 34% and PN→FSI were 12%. Reciprocal connections between PNs and FSIs were set to 16%. Electrical connections between FSIs were set to 8%. When a pair of FSIs were electrically coupled, they had a 50% probability of a unidirectional chemical synapse, or a 25% probability of bidirectional synaptic connectivity. FSI pairs not electrically coupled had a 19% probability of unidirectional connectivity, and a 3% probability of bidirectional probability. These connectivity numbers in our model resulted in an overall synaptic FSI→FSI connectivity of 26%, of which 20% was unidirectional and 3% bidirectional. These probabilities resulted in the following intrinsic connectivity in the model: each PN received 24.98 ± 9.5 (mean ± SD, throughout the paper unless otherwise indicated) excitatory connections from other PNs, and 42.6 ± 12.9 inhibitory connections from FSIs; each FSI received 214.8 ± 58 excitatory connections from PNs, and 21.6 ± 7.4 inhibitory connections from other FSIs. Axonal conduction delay on all connections were distance dependent. See Tables 5, 6 for details.

View this table:
  • View inline
  • View popup
Table 5.

Connection probabilities, PN-PN connections

View this table:
  • View inline
  • View popup
Table 6.

Connection probabilities, FSI-FSI and PN-FSI connections

BL afferents

PFC, PR, and ER strongly project to BL (McDonald and Mascagni, 1997; Shi and Cassell, 1999; Vertes, 2004), so we tailored our extrinsic inputs to match their spiking. To this end, we created surrogate spiking activity derived from putative PNs simultaneously recorded from the PFC, PR, and ER in vivo during the QW state (Headley et al., 2015). To do this, we created multiple surrogate spike trains where each was a random combination of the temporal firing rate profile of one neuron and the mean firing rate of another. The firing rate profile was the instantaneous firing rate of each unit, which was calculated by taking the inverse of the interspike interval in 1-ms steps. The instantaneous firing rate series was rescaled to have a mean rate of 1 Hz. Mean firing rate was the number of spikes dividing by the recording duration. A new spike train was created by multiplying one of the rescaled instantaneous rate series by a mean firing rate randomly drawn from the same population of units, and then passing each time step to a Poisson random number generator. When repeated multiple times, this creates an ensemble of spike trains that matches the distribution of firing rates, autocorrelations, and cross-correlations from the originally recorded ensemble. We generated an ensemble of 1800 Poisson spike trains to represent the activities of cortical regions upstream of BL, and each of these projected to a random 40 PNs, resulting in ∼95% of the PNs receiving at least one extrinsic input. The same extrinsic input also connected to FSIs that were within 300 µm of these 40 PNs with a probability of 0.3%; this resulted in each extrinsic input connecting to an average of 6 FSIs. Thus, PNs and FSIs (Hübner et al., 2014) received 2.97 ± 1.7 and 3.4 ± 1.9 extrinsic inputs, respectively.

To reproduce membrane potential fluctuation seen in vivo, we used a point conductance model that mimics stochastic background synaptic activity using an Ornstein–Uhlenbeck process (Destexhe et al., 2001). Specifically, stochastic background input Ibck had two independent components, excitatory and inhibitory, for both PNs and FSIs, modeled as follows:Embedded Image (7)where Embedded Image and Embedded Image are time-dependent excitatory and inhibitory conductances, respectively; Embedded Image and Embedded Image are respective reverse potentials. The two are modeled as Ornstein–Uhlenbeck processes described below:Embedded Image (8a) Embedded Image (8b)where Embedded Image and Embedded Image are average conductances, Embedded Image and Embedded Image are time constants, Embedded Image and Embedded Image are noise “diffusion” coefficients, Embedded Image and Embedded Image are Gaussian white noise of zero mean and unit SD (for these parameter values, see Table 7).

View this table:
  • View inline
  • View popup
Table 7.

Parameters related to point-conductance model

The above two stochastic differential equations can be numerically modeled by using the following update rule:Embedded Image (9a) Embedded Image (9b)where Embedded Image and Embedded Image are normal random numbers. Embedded Image and Embedded Image are amplitude coefficients with Embedded Image and Embedded Image

Table 7 lists the parameters used in the point-conductance model. For each neuron, the excitatory and inhibitory conductances were, respectively, 3.2 ± 3 and 21 ± 8 nS for PNs and 1.2 ± 0.1 and 5.7 ± 2.6 nS for FSIs. However, the background inputs to FSIs was weaker to ensure FSIs spike with input from PNs, but not with solely background inputs (Börgers and Kopell, 2003; Economo and White, 2012; Lee and Jones, 2013).

Calculation of LFP

Gamma rhythms are detected by extracellular recordings of LFPs within the brain. In contrast, biophysical models of neuronal networks have typically detected gamma rhythms using other measures, such as network spiking rates (Brunel and Hakim, 1999; Brunel and Wang, 2003; Börgers et al., 2005; Economo and White, 2012; Chalk et al., 2016; Hoseini and Wessel, 2016; Palmigiano et al., 2017) or membrane voltages (Traub et al., 2000; Bathellier et al., 2006; Neymotin et al., 2011), both of which are only indirectly linked to the LFP (Buzsáki et al., 2012; Schomburg et al., 2012). The present study modeled LFPs using a first principles approach.

We first recorded transmembrane ionic currents from each compartment of the model cells using the extracellular mechanism in NEURON (Carnevale and Hines, 2005; Parasuram et al., 2016). The extracellular potential arising from each neuronal compartment was then calculated using the line source approximation method, which provides a better approximation than point sources (Gold et al., 2006; Schomburg et al., 2012). The extracellular potential of a line compartment was estimated asEmbedded Image (10)where, I denotes the transmembrane current from just that compartment, Δs the length of the line compartment, r the radial distance from the line, h the longitudinal distance from the end of the line, and Embedded Image the distance from the start of the line (Holt, 1998; Parasuram et al., 2016). We chose a conductivity σ of the extracellular medium to be 0.3S/m (Goto et al., 2010; Einevoll et al., 2013). These individual extracellular potentials were summed linearly (Lindén et al., 2013) at 1-ms resolution, to obtain the LFP Embedded Image for an N-neuron network with n-compartment-cells using the equationEmbedded Image (11)where Ni denotes i th compartment of N th neuron in the network. All 27,000 neurons contributed to the LFP in our study, different from previous studies (Bezaire et al., 2016) where only a subset of neurons did. This permitted investigation of both the individual contribution and correlation of all modeled neurons to the LFP at any of the electrode locations.

Since our model spans the entire spatial extent of BL, but only has a fraction of the total number of neurons, the density of neurons in our model is lower than the in vivo case. This means that neurons will, on average, be farther from the simulated LFP electrode than would occur in the actual BL. To correct for this, we rescaled the LFP by a correction factor that was estimated following Lindén et al. (2011), who suggest that LFP scales as the square root of the number of neurons, N, for uncorrelated synaptic inputs. However, that scaling factor was derived from a network of pyramidal neurons positioned in a disk with uniformly oriented dendrites, and was not for density per se but actually the radius of the disk with density held constant. To determine the scaling factor in our model, we systematically varied the number of cells used to calculate the LFP from a full model run. Varying the cell count from 1000 up to 27,000, we found that the SD of the LFP scaled with density following N0.67. Since cell density in rats is reported to range from 2.5 × 104 to 2 × 105/mm3 (Tuunanen and Pitkänen, 2000; Salm et al., 2004; Pêgo et al., 2008; Rubinow and Juraska, 2009; Chareyron et al., 2011), while the model density is 9840 neurons/mm3 (27,000 neurons, 1.4 × 1.4 × 1.4 mm), the LFP correction factor would correspondingly range from 1.9 to 7.5 [Embedded Image to Embedded Image ]. We chose the average of this range, 4.7, to scale the model LFP.

Multiple extracellular electrodes

To investigate the spatial propagation of model gamma we used a 9 × 9 × 9 grid of LFP electrodes evenly spaced throughout the network. Electrode sites were spaced at 125-µm intervals and were at least 200 µm away from the edges of the model. The LFP on each electrode was computed the same as in Equation 11.

Model experiments and statistics

Spectral and cross-correlation analyses

Unless otherwise indicated, spectral decompositions were performed using Morlet wavelets ranging from 1 to 256 Hz in quarter octave steps. The width of the wavelet was seven cycles. To measure the amplitude at a particular frequency and time, we took the absolute value of the complex valued frequency domain representation of the signal. Phase was measured as the angle of the frequency domain representation.

The distribution of gamma burst properties were compared between the model and in vivo cases in Figure 2F. For both cases, wavelet power spectrograms were calculated and the mid-gamma band (64 Hz) was isolated. Peaks were detected and segmented in the time series of the mid-gamma power using a watershed function. To prevent spurious detections of peaks, before passing the data to the watershed function, we discretized gamma power into bins of half the average change in power between cycles. Once burst peaks were detected, we extracted their peak power (before discretization) and duration. Peak powers were transformed into percentile ranks and binned in steps of five percentiles, while durations were binned in steps of 5 ms.

Figure 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 2.

The model LFP exhibits gamma oscillations that are similar to those seen in vivo. A, The LFP was measured at the center of the model. B, The power spectrum of the model LFP (red) is compared with the spectra from the BL of seven rats during QW (gray). Driving the model with a homogenous Poisson input (green) still induces gamma oscillations. Shaded regions denote SD. C, The frequency where the gamma bump peaked was measured for each subject (black bars), the probability distribution was fit with a normal curve (gray). The model’s peak frequency was 64 Hz (red line), which fell within the in vivo distribution. D, Gamma oscillations occurred as intermittent bursts in vivo and in the model. There was a stronger low-frequency component in vivo (black), but filtering that out (yellow) revealed a comparable signal amplitude to our model (blue). E, The model had a similar wavelet spectrogram to that observed in vivo. F, Gamma bursts were detected in the wavelet power spectrograms and categorized based on their amplitude percentile and duration. The relationship between these features was characterized by a probability distribution. Bursts arising from the model or in vivo exhibited a similar distribution of durations when stratified by amplitude, with higher amplitude bursts tending to last longer. G, (1) Each simulated cell type was placed in current clamp and driven with a frequency modulated sinusoid. (2) Neurons responded to this input with an oscillating membrane potential that decayed with increasing frequency. (3) The impedance spectrum of the neuronal responses did not exhibit any peaks in the gamma band. H, (1) An extracellular electrode was placed near either a PN or FSI and synaptic inputs to those neurons were driven (stim). (2) The extracellular field responding arising from a single EPSC or IPSC was measured for each cell, and normalized to the strength of the PN EPSC. There was no bump in the gamma band. (3) Driving the synaptic inputs with Poisson trains at rates indicated by the x-axis did not show any obvious boost in power in the gamma range either.

For the single neuron resonance analyses (Fig. 2G,H), we calculated spectrums using the Welch Periodogram method (pwelch in MATLAB). The Hamming window taper size was 500 ms, in steps of 250 ms.

For the coherence spectrum between PN and FSI spiking (Fig. 7D), the cross-correlation and autocorrelation between these spike trains were calculated with 1-ms resolution. Then, they were converted to the frequency domain using the fast Fourier transform (FFT) and used to calculate coherence:Embedded Image (12)

Since both the PNs and FSIs exhibited gamma periodicity in their autocorrelation function, it was desirable to correct for this when estimating their cross-correlation. To do this, one can take the inverse FFT of the coherence spectrum, which yields the cross-correlation function in the time domain, but with periodicities arising from the autocorrelations factored out:Embedded Image (13)

Entrainment to LFPs

For calculating the entrainment and preferred phase of cells, we first bandpass filtered the LFPs in the frequency band of interest using a two-pole Butterworth filter implemented with the MATLAB function filtfilt, which performs forward and backward filtering to minimize phase distortion. A Hilbert transform of the resulting signal was then computed to determine the phase and amplitude at each instant (Amir et al., 2018). This was used to assign a phase to each spike from a neuron. Entrainment of spiking to the LFP was computed as the mean resultant vector length using the Circular Statistics toolbox in MATLAB (Berens, 2009). This approach matched that in Amir et al. (2018). For our analysis of entrainment when disconnecting specific connections (see Results, Circuit elements critical for gamma genesis), we instead used the pairwise phase consistency (PPC) measure, which is insensitive to changes in spike count arising from alterations in firing rate (Vinck et al., 2010).

When computing correlation between distance and entrainment (Fig. 3D), cells were selected only within 100 µm from the four longest diagonals within the modeled cuboid structure, to minimize edge effects. To directly compare the decay in entrainment between PNs and FSIs (Fig. 3D), we fit a generalized linear model with a log link function and factors distance, cell type, and their interaction.

Figure 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 3.

Spiking in the model was entrained to the gamma rhythm similarly to cells in vivo. A, Example trace of the calculated extracellular field from the center of the model and a raster plot of corresponding spiking from neurons within 300 µm (1026 PNs, 111 FSIs). Gamma bursts were evident, particularly following bandpass filtering around the peak gamma frequency. The group firing rates of FSIs and PNs exhibited modulations by the gamma rhythm, with the effect being more evident for FSIs. Spikes from individual neurons (each row is a different cell) did not spike on every gamma cycle. B, Left, Spiking from units (AP) within 300 µm of where the extracellular field was measured (LFP) were associated with the phase of the ongoing gamma oscillations. Right, Distribution of gamma phases when spikes occurred for three example PNs and FSIs, sorted by the strength of their phase locking. Dashed gray schematizes the corresponding gamma cycle. C, Left, The distribution of mean phases across each cell type. Both types of PNs preferentially spiked near the trough of the gamma oscillation, while FSIs spiked shortly thereafter during the ascending phase. Right, The distribution of phase-locking across cell types matched what we found for BL neurons in vivo. D, Left, Spiking was sampled from neurons spanning 800 µm away from the point where we measured the extracellular field, and neurons were sorted by their distance from that spot. Right, The farther a neuron was from the site were gamma was measured, the less phase-locked it was. Black lines are exponential fits. E, Left, The extracellular field arising from individual neurons was singled out to examine their distance-dependent contribution. Right, Neurons farther from the recording site contributed less to the extracellular field, and this fell off faster than the decay in the resultant vector in panel D.

Detection of gamma bursts

For gamma burst extraction, we first applied zero-phase 60 − 80Hz Butterworth bandpass filter of order 4, followed by the Hilbert transform to obtain the analytic signal of the LFP (Dvorak and Fenton, 2014). Instantaneous gamma amplitude was the modulus of the analytic signal. Gamma burst occurrence was detected when the amplitude exceeded its mean by 2 SDs. Burst duration was defined as the duration during which the amplitude was >25% of the peak value for the burst.

Propagation of gamma bursts

Using the 9 × 9 × 9 grid of LFP recording sites (see Multiple extracellular electrodes), we examined the properties of gamma bursts as spatiotemporal events. Near the edges of our model, we found that the gamma bump in the power spectrum shrunk, likely due to the decrease in the number of neurons contributing to the LFP and decrease in the absolute number of connections between PNs and FSIs. Therefore, to compensate, we multiplied the Z-score of the LFP amplitude from each recording site by a correction factor F with the form:Embedded Image (14)where GammaAmp is the area under the curve for the bump in the gamma band after removing a power law fit (Neske and Connors, 2016), TotalAmp is the total area under the curve for the power spectrum, j denotes the recording electrode, and middle indicates the value for the electrode in the center of the model. The Z-scored LFP Hilbert transformed amplitude at each site was multiplied by this factor.

We then identified gamma bursts as spatiotemporal events using a four-dimensional watershed algorithm that identifies contiguous regions of elevated gamma power that are convex. To minimize spurious detections of gamma bursts driven by small variations in the power, we discretized the Z-scored LFP Hilbert transformed amplitude into bins of 0.1. In order for a region to count as a gamma burst, its peak amplitude had to exceed two Z-scores. The borders of the burst region were delineated when power dropped by 25% of the peak. These were criteria were intentionally similar to the ones used for detecting gamma bursts on single electrodes (see above, Detection of gamma bursts).

After identification of a gamma burst, its properties were defined as follows. Gamma burst peak (Fig. 6B1) was defined as the maximum Z-scored LFP amplitude within the burst. Its time duration (Fig. 6B2) was calculated as the time duration during which at least one recording site within the watershed region was above the 0.25% of the peak Z-score. Gamma burst volume (Fig. 6B4) was defined as the number of recording sites that had reached the burst threshold with the same burst region.

To estimate propagation of gamma bursts, we identified the gamma region composed of electrodes with the same watershed label at each time step. Then, a burst center point was defined as the mean coordinate from all electrodes belonging to the same burst region. Gamma burst path length (Fig. 6B5) was calculated as the total distance that a burst center could travel in space.

To quantify instantaneous gamma burst synchrony (Fig. 6B3), we calculated the phase locking value (PLV) by using the unbiased PLV (Aydore et al., 2013) of 27 electrodes within a 3 × 3 × 3 grid located at the gamma burst peak.

Statistical analysis

Values are mean ± SD unless otherwise stated. Adjusted R 2 value were reported when data points were fit by a curve in Figure 3D,E. For Figure 5E, ANOVAs were performed on fitted general linear models and p values were reported. Two-tailed rank-sum tests were performed and p values were reported for Figure 7B,D. Watson–Williams tests by using Circular Statistics toolbox in MATLAB (Berens, 2009) were performed and p values were reported for Figure 7C. ANCOVA was tested and p values were reported for Figure 7F.

Original NEURON code for the model is available in ModelDB public database with access number 247968 (http://senselab.med.yale.edu/ModelDB/).

Results

Validation of the BL model

A biophysically realistic computational model of the rat BL was developed based on the known anatomy and cellular neurophysiology (for details and supporting empirical references, see Materials and Methods). The model was a cube with PNs and FSIs distributed randomly throughout (Fig. 1A,B). The dendritic compartments for both PNs and FSIs were oriented randomly and spanned several hundred microns from the cell body (Fig. 1C). Using this model, we tracked both the intracellular and extracellular signals associated with network events (Fig. 1D1). The model network included extrinsic afferents contacting both PNs and FSIs and intrinsic connections (Fig. 1D2), in keeping with previous anatomic (Paré et al., 1995; Muller et al., 2006) and physiologic (Washburn and Moises, 1992; Rainnie et al., 1993; Woodruff and Sah, 2007) observations.

Two types of PNs with marked (PNA) or negligible (PNC) spike frequency accommodation were used based on previous reports and our own investigations (see Materials and Methods; Fig. 1E,F). Our in vitro recordings revealed 41/56 PNA and 15/56 PNCs, so our model used similar proportions (71% PNA, 29% PNC). We also included FSIs, which are the prevalent type in BL (McDonald and Mascagni, 2001; Fishell and Rudy, 2011), well characterized biophysically and synaptically (Rainnie et al., 2006; Woodruff and Sah, 2007), and are important for gamma genesis in cortical circuits (Bartos et al., 2007). Point conductance noise levels were adjusted so that the firing rates of PNs and FSIs matched those observed in vivo (Fig. 1G).

The LFP was calculated using the relative location of cellular compartments, their transmembrane currents, and an estimate of the extracellular conductivity (see Materials and Methods). The model LFP exhibited several well-established aspects of real extracellularly recorded APs and excitatory postsynaptic potentials (EPs; Fig. 1H,I). PNs produced stronger dipoles than FSIs, resulting in greater extracellular potentials. Near the soma, the simulated extracellular AP was negative going, larger for PNs, and faster for FSIs (Fig. 1J1,2, note scale difference). A unitary EPSC evoked in a dendritic branch produced an extracellular negativity near the dendrite that was substantially smaller than that arising from an AP (Fig. 1J2). Additionally, the amplitude of extracellular APs decreased with distance from the soma similarly between the model and in vivo cases (Fig. 1H,K).

Emergence of gamma oscillations in the BL model

We next examined whether the spontaneous LFP at the center of our BL model (Fig. 2A) matched in vivo BL recordings when rats were in a QW state. A characteristic of the LFP in vivo is that its power spectrum can be fit by a 1/fa power law (Buzsáki et al., 2012). BL recordings from Long–Evans rats showed a spectrum with the 1/fa falloff, but also a bump in the gamma band between 54 and 70 Hz (Fig. 2B, gray lines). A gamma band signal was present in the model’s LFP, with similar power (101 mV2/Hz) and peak frequency (64 Hz; Fig. 2B,C, red line). However, the model’s spectrum had substantially less power below 30 Hz. This may reflect the absence of volume conducted activity or a lack of interactions with extrinsic structures. The former possibility was studied by Łęski et al. (2013), who showed that lower-frequency LFP components have a larger spatial reach and extended further outside the active population than high-frequency components. It has also been found that increasing the correlations between extrinsic afferents boosts power in the low-frequency band (Łęski et al., 2013; Hagen et al., 2016). Our model likely underestimates the contribution of correlated extrinsic inputs since our naturalistic spike trains only represented a fraction of the structures projecting to the BL, each presumably with their own correlated spiking. Consistent with this, we found that replacing the naturalistic correlated spiking processes of our extrinsic afferents with rate matched independent Poisson processes further reduced the power of low frequencies (Fig. 2B, green line). Importantly, gamma power from 50 to 80 Hz was unaffected by this change, supporting the hypothesis that it is locally generated and not imposed by upstream structures or volume conducted from adjacent regions. In addition, the gamma oscillations in our model occurred as intermittent bursts, like those found in BL in vivo (Fig. 2D,E). The probability distribution of gamma bursts with respect to their duration and amplitude was similar between the model and those recorded in vivo (Fig. 2F).

Neither the resonance properties of single cells (Fig. 2G1–3), nor the spectral content in EPSCs and IPSCs (Fig. 2H1,2) contributed to the gamma in the model. Also, there was no enhancement in the gamma-band when synaptic potentials were driven as Poisson processes with rates ranging from 1 to 200 Hz (Fig. 2H3). This leaves interactions between PNs and FSIs as the most likely candidate causing gamma.

Entrainment of neurons to gamma in the LFP

We examined whether BL neurons showed a similar entrainment (Fig. 3A) to in vivo recordings (Amir et al., 2018). For neurons within 300 µm of the point where we calculated the LFP we measured their resultant vector (Fig. 3B,C). Both types of PNs were entrained to the gamma rhythm, spiking during the trough, with a mean resultant vector of 0.19 ± 0.07 (mean ± SD) for PNAs, and 0.17 ± 0.06 for PNCs. On the other hand, FSIs tended to fire during the ascending phase, and exhibited similar mean entrainment, 0.17 ± 0.04.

PNs and FSIs exhibited decreased entrainment to gamma at longer distances (Fig. 3D; PN: 50% falloff at 409 µm, adjusted R 2 = 0.37, df = 350; FSI: 50% falloff at 355 µm, adjusted R 2 = 0.67, df = 146). A generalized linear model also indicated a distance-dependent falloff (t statistic = −15.29, df = 496, p = 1.63 × 10−43) and no difference in falloff between cell types (distance × FSI, t statistic = −0.95, df = 496, p = 0.34). To rule out the possibility that the decrease in entrainment was driven purely by the distance-dependent decay in a neuron’s contribution to the extracellular potential, we examined the contributions of individual PNs and FSIs at different distances (Fig. 3E). Arguing against this, the falloff in the contribution of both PNs and FSIs to the extracellular potential was much sharper than the decay in their entrainment (PN: 50% falloff at 125 µm, adjusted R 2 = 0.59, df = 1435; FSI: 50% falloff at 113 µm, adjusted R 2 = 0.67, df = 749).

Circuit elements critical for gamma genesis

To identify the circuit elements that support the generation of gamma oscillations, we systematically eliminated each connection type one at a time. Since the firing rates of PNs and FSIs could change dramatically in the disconnected models, we adjusted the point conductance noise of both cell types to bring them back into agreement with the original unperturbed model. For each perturbed model, we measured the power spectrum of the simulated LFP and unit entrainment to frequencies ranging from 1 to 256 Hz with the PPC measure, which is unbiased for changes in spike count (Vinck et al., 2010). To clarify the changes, we saw in the power spectrum, we calculated the cross-correlation between bursts of PN activity (time points where the firing rate of the PN population exceeded the 75th percentile) and spiking in the FSI population.

Removal of the FSI→PN connection eliminated the mid-γ bump in the LFP (Fig. 4A1), but increased power at frequencies below 20 Hz. Consistent with this, the entrainment of PNs by low frequencies was increased as well (Fig. 4A2), with a particularly prominent peak at 6 Hz (Fig. 4A2, inset). Also mirroring the LFP, PN and FSI entrainment to gamma was substantially reduced. This loss was evident in the response of the FSI network to bursts of PN activity, which now lacked a prominent gamma periodicity (Fig. 4A3). Instead, FSIs exhibited a slow undulation in firing that was capped with a transient burst of spikes in response to PN firing.

Figure 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 4.

Connections between PNs and FSIs were critical for genesis of the gamma rhythm. A, Removal of connections from FSIs to PNs. (1) Comparison of power spectrum between the original and perturbed model. Gray traces are derived from the unmodified model and are the same across all graphs in each row. The green traces are the power spectrum with the particular connection removed. (2) Entrainment to LFP oscillations for different cell types. (3) Average modulation of FSI firing rate triggered on PN population bursts. B, Removal of connections from PNs to FSIs. C, Removal of connections between FSIs. D, Removal of connections between PNs. For all graphs the solid line is the mean while the shaded region denotes ±1 SD.

PN→FSI connections were also critical to gamma generation (Fig. 4B1). Their removal reduced the entrainment of PNs and FSIs in the gamma frequency range (Fig. 4B2). Compared to the previous model, the increase in power in the low-frequency band was smaller, so too was entrainment by PNs (Fig. 4B1,2). FSIs, on the other hand, were not entrained to the low-frequency rhythm, suggesting that it originated with the PN network, from which it was now disconnected. As expected, FSIs no longer responded to bursts of PN activity (Fig. 4B3).

The elimination of gamma by removing either PN→FSI or FSI→PN connections is in agreement with a PING based mechanism for gamma genesis (Whittington et al., 2000). This does not rule out the importance of connections within the PN and FSI populations. Since these connections were also present in the model, we also examined their contribution.

An interconnected network of FSIs that receives tonic excitatory drive can produce gamma oscillations, so called interneuron network gamma (ING; Whittington et al., 2000), which may occur in the model. Eliminating FSI→FSI connections (including gap-junctions), shifted the peak frequency to 23 Hz (Fig. 4C1). Likewise, both PNs and FSIs shifted their maximal gamma entrainment to the new peak frequency (Fig. 4C2). Since the timescale of feedback inhibition affects the frequency of gamma (Wang and Buzsáki, 1996), the lower frequency likely reflects the increased duration of the PN-evoked FSI burst, which arose from the removal of feedback inhibition generated by FSIs onto themselves (Fig. 4C3).

Removing PN→PN connections had a modest effect on gamma generation, shifting the peak frequency slightly lower, which was also reflected in the entrainment of PNs to gamma (Fig. 4D1,2). Notably, there was a large drop in spectral power at lower frequencies, suggesting that it partly depends on recurrent activity within the PN network. Befitting these minor changes, FSI responses to PN bursts were only marginally affected (Fig. 4D3).

Altogether these results suggest that reciprocal interactions between PNs and FSIs support the generation of gamma oscillations in BL, and that their frequency is affected by connections within the PN and FSI populations.

Synaptic basis of gamma coordination in single cells

Since gamma in our model was similar to the in vivo case, we explored how the convergence of excitatory and inhibitory synapses affects a neuron’s entrainment to gamma. We chose 300 PNs and FSIs at random and recorded their EPSCs, IPSCs, spike times, and surrounding LFP (sampled ∼200 µm away; Fig. 5A). We limited ourselves to neurons that fired over 100 spikes in the 165 s of simulation time, so only PNcs were present in the PN sample.

Figure 5.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 5.

EPSCs and IPSCs entrained neurons in the gamma band. A, Left, For PNs and FSIs we measured the membrane voltage (ic) and their transmembrane currents while voltage clamped (vc) at the reversal potentials for either EPSCs (0 mV) or IPSCs (−75 mV). Middle, Traces from an example PN. In addition to the transmembrane currents and voltages were recorded, we also bandpass filtered the EPSC and IPSC traces around the gamma peak. Right, Examples traces from an FSI. B, The mean power spectrums of the EPSC and IPSC traces for each cell type. PNs showed greater power for IPSCs across all frequencies compared with EPSCs. FSIs had similar power levels across the spectrum for EPSCs and IPSCs. The spectrum of IPSCs from both cell types contained a bump in the gamma band. C, The coherence of the EPSCs and IPSCs for each cell type. PNs had a dip in coherence around the peak gamma frequency, while FSIs showed a prominence. Shaded regions denote SD. D, Synaptic currents were aligned to the phase of gamma cycles and centered on ones that either had, or did not have, spiking (thin blue and red lines, solid lines from cycles with spikes, dashed for cycles without spikes). The difference in the PSCs between the spike and no spike conditions is shown below (thick red and blue lines). IPSCs were modulated by gamma phase for both cell types, while only FSIs exhibited modulation of their EPSCs to gamma. E, For each cell we plotted the number of excitatory (x-axis) or inhibitory (y-axis) connections they received and their corresponding entrainment. Each scatter point is a different cell. The shaded background is the result of a linear fit that estimates the resultant vector strength based on the number of inhibitory connections, excitatory connections, and their product (interaction term). Contours in black indicate the actual resultant vector values associated with each color.

Given that synaptic inputs constitute a major determinant of neuronal firing, we first assessed the spectral content of the EPSCs and IPSCs impinging on PNs and FSIs (Fig. 5B). On average, IPSCs exhibited a pronounced bump in the mid-gamma band for both cell types. This contrasted with EPSCs, which lacked increased gamma-band power for PNs, and only exhibited a shallow increase in FSIs. The presence of a gamma bump in both synaptic current power spectra for FSIs, but not PNs, suggested that the coherence between excitatory and inhibitory synaptic currents differs between cell types. Indeed, EPSC-IPSC coherence at the gamma frequency was higher in FSIs than in PNs (Fig. 5C).

To determine whether this effect was related to spiking during gamma (Fig. 5D), we calculated the mean EPSCs (Fig. 5D1, red lines) and IPSCs (Fig. 5D2, blue) during gamma cycles with (solid) and without spikes (dashed), and their difference (Fig. 5D3). For PNs, spiking was associated with an increase in EPSCs, which started rising during the latter half of the prior gamma cycle, and peaked at the trough of the oscillation, when PNs normally fire (Fig. 3C). As was evident in the power spectrum for the PN EPSC current (Fig. 5B, left, red), no gamma periodicity was evident in the EPSC trace. On the other hand, the IPSCs tracked the ongoing LFP gamma rhythm (Fig. 5B, left, blue). While there was a slight decrease in IPSC strength during the first half of the gamma cycles containing spikes, this reversed itself during the latter half, with inhibition increasing following the phase when PN spiking normally occurs (Fig. 3C), likely reflecting feedforward activation of the FSI network.

FSIs exhibited a different pattern of postsynaptic currents during gamma cycles in which they spiked (Fig. 5D, right). The EPSCs were entrained to the ongoing gamma rhythm and elevated just before the neuron spiked, during the ascending phase of the gamma cycle. Also, unlike PNs, FSIs showed a decrease in IPSC amplitude throughout the gamma cycle in which they spiked, followed by a return to baseline. One similarity between PNs and FSIs was that their IPSCs were similarly entrained to the gamma rhythm. These results suggest that PN spiking during a gamma cycle is mostly dependent on the presence or absence of an EPSC, while FSIs are more likely to spike when both EPSC amplitudes are increased and IPSCs are decreased. Despite these differences, it is worth noting that PNs and FSIs tended to spike during the phase of the gamma cycle when EPSCs were at their maximum and IPSCs were at their minimum. Our finding that spiking during gamma cycles was accompanied by both EPSCs and disinhibition of IPSCs supports the claim that the gamma in the BL network arises through a PING mechanism, rather than the recurrent-excitation-inhibition (REI) mechanism found in models of visual cortex (Chariker et al., 2018).

Since spiking during gamma depends on the pattern of postsynaptic currents, it could be that a neuron’s entrainment to gamma depends on how many excitatory and inhibitory synapses it receives. For each neuron, we measured their entrainment to gamma using the resultant vector metric and determined how it varied as a function of their number of excitatory and inhibitory afferents (Fig. 5E). For PNs and FSIs there was a significant positive relationship between entrainment and the number of inhibitory afferents received (PN: F(1,296) = 35.8, p = 6.3 × 10−9; FSI: F(1,296) = 770.3, p = 2.3 × 10−84). However, the number of excitatory afferents did not alter entrainment in PNs (F(1,296) = 0.16, p = 0.69), but decreased it in FSIs (F(1,296) = 264.8, p = 5.7 × 10−43). Notably, FSI entrainment was enhanced when both inhibitory and excitatory afferents increased in tandem (F(1,296) = 5.1, p = 0.025). Overall, the dependence of PN and FSI entrainment on the number of inhibitory connections agrees with the relatively strong gamma periodicity of IPSCs.

Spatiotemporal properties of gamma bursts

Gamma oscillations not only persist in time but extend across space. This introduces uncertainty about the amplitude and duration of gamma bursts when recording from single sites. This is not a problem for our model, which can support an arbitrary number of recording sites arranged in any desired configuration. Thus, we characterized gamma bursts as spatiotemporal phenomena.

A 9 × 9 × 9 grid of recording sites with an interelectrode spacing of 125µm spanned the entire BL model (Fig. 6A). Spatiotemporal gamma bursts were identified in a four-dimensional space (three spatial and one temporal). We measured gamma amplitude at each site and identified spatially contiguous regions of elevated gamma power at each time step. To count as a burst power had to peak above 2 SDs of the mean and the edges of the burst were thresholded at 25% of the peak power.

Figure 6.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 6.

Spatiotemporal properties of gamma bursts. A, Left, A 9 × 9 × 9 array of electrodes was inserted into the model and the envelope of the gamma band was extracted. Right, A cube of amplitude values was generated and short lasting, spatially localized, increases in gamma power were readily evident. B, We examined the properties of these bursts, in particular their (1) amplitude, (2) duration, (3) phase synchrony at their peak, (4) number of sites they encompassed, (5) length of the path they took through the BL from initiation to termination, and (6) speed that they traversed that path.

By surveying the entirety of each gamma burst, we could unambiguously measure their properties (Fig. 6B). Gamma burst peak amplitudes hugged the cutoff threshold we had set at two Z-scores, with a mean of 2.8, and <1% above 5 (Fig. 6B1). Burst durations were 108 ms on average, enough time to encompass several gamma cycles (Fig. 6B2). We also measured the mean phase locking of the LFP in the gamma band surrounding the location of peak power in the gamma burst. Most bursts had an average PLV of 0.87 (Fig. 6B3), indicating that at the center of the burst, adjacent LFP recording sites were highly phase coherent.

Tracking across multiple recording sites revealed the spatial properties of gamma bursts. Localized increases in gamma power could start in one part of the model, spread to adjacent areas, and then dissipate. Most bursts were local, encompassing 44.3 recording sites on average out of a possible 729 total (6% of the BL volume; Fig. 6B4). If the mean volume of a gamma burst were arranged as a sphere, it would have a radius of ∼270 µm. Despite being local, the center of the burst could move throughout the BL, with bursts on average traversing 472 µm from their start to end times (Fig. 6B5). Taking the distance a burst traveled, and dividing that by the burst duration, we could measure the average speed of its center, which was 4.4 µm/ms (Fig. 6B6), or 68 µm per gamma cycle (assuming a cycle duration of ∼15 ms).

These results highlight the fact that gamma bursts, even in relatively small nuclear structures like BL, can be highly localized. However, they are not spatially confined; they can emerge at one location and terminate in another. Such propagation may mediate interactions between distant cell populations.

Computation in the BL network

To explore how the microcircuitry generating gamma affects the interaction between BL neurons, we randomly assigned PNs in the network to one of two populations (group 1 and group 2), each receiving a different set of extrinsic afferents (Fig. 7A1). The simulated LFP was recorded from the center of the model, along with the spiking activity, IPSCs, and EPSCs of nearby neurons. Each population had its extrinsic afferents driven as Poisson processes at a fixed rate (Fig. 7A2–4). For these simulations we used a reduced BL model limited to 1000 neurons, with similar connectivity rules and proportions of cells types as found in the full model.

Figure 7.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 7.

Gamma promoted synchrony between PNs responding to extrinsic afferents and the competition between weakly and strongly driven populations. A, Two separate ensembles of PNs were targeted by extrinsic afferents (group 1, orange vs group 2, green), that received different subsets of afferents. Groups could be driven independently, with (2) group 1 receiving 20-Hz stimulation while group 2 had 5 Hz; (3) group 1, 20 Hz and group 2, 20 Hz; or (4) group 1, 5 Hz and group 2, 20 Hz. B, (1) Spikes were recorded from both groups of PNs and (2–4) we calculated their cross-correlation functions. C, (1) Spikes from PNs and the mean EPSC across all FSIs were recorded. (2–4) Probability histograms of the EPSC gamma phases associated with spiking from either group 1 or group 2. The strongly driven group exhibits greater entrainment. D, (1) The coherence between PN spiking and FSI spiking for each group. (2–4) FSIs were more coherent in the gamma-band with the population that received the strongest afferent drive. Inset, PN preceded FSI spiking by 4 ms on average. E, (1) Group 1 receives 20-Hz extrinsic drive while group 2 receives a varying amount. Increasing the firing rate of extrinsic afferents onto group 2 diminished the firing rate of group 1. (2) The difference in firing rate between these groups was weakened when connections from PNs to FSIs were removed. F, Analysis of competition between groups in the full model when they either (1) fully overlap or (2) are entirely segregated in space. G, top, Gamma cycles were divided into four phases and for each phase the spike probability was measured as a function of the EPSC strength. Bottom left, For the weakly driven group there was increased sensitivity for EPSCs during the trough and ascending phases of gamma, but a dampening of responsiveness during the descending and peak phases. Bottom right, When connections from PNs to FSIs were removed, the discrepant sensitivities for EPSCs across the four phases were diminished. Graphs are either mean or mean ± SE when a shaded region is present.

We first examined how increases in afferent drive to group 1 PNs affected synchronous spiking within each group (Fig. 7B1). Randomly selected pairs of neurons had their spike cross-correlations calculated. A robust gamma frequency periodicity was evident in the cross-correlation function of both groups (Fig. 7B2–4), with the group receiving the strongest afferent drive exhibiting greater synchronous spiking [group 1 > group 2, rank-sum test z value = 11.8, p = 6.8 × 10−32 (Fig. 7B2); group 2 > group 1, z value = −13.7, p = 1.6 × 10−42 (Fig. 7B4)]. When both groups were driven with equal strength there was not significant difference in their spiking synchrony (group 1 = group 2, z value = 0.2, p = 0.81; Fig. 7B3).

Previous analyses revealed that EPSCs in FSIs had a gamma component (Fig. 5B,D), so we next determined whether the more strongly driven PN group (group 1) contributed to this effect. Spikes times were recorded in tandem with the phase of the mean EPSC current in the gamma band across all FSI neurons (Fig. 7C1). PNs tended to spike before the peak of the EPSC, with a quarter cycle lag (Fig. 7C2–4). Whichever group was driven strongly had greater entrainment to the FSI EPSCs [group 1 > group 2, group 1 spike resultant vector, RV = 0.54, group 2 RV = 0.5, Watson–Williams test, p = 1.1 × 10−16 (Fig. 7C2); group 2 > group 1, group 1 RV = 0.5, group 2 RV = 0.54, p = 4.1 × 10−13 (Fig. 7C4)]. When both groups were driven equally their entrainment with FSI EPSCs did not differ (group 1 = group 2, group 1 RV = 0.56, group 2 RV = 0.56, p = 0.17; Fig. 7C3).

Given that FSI spiking depends on EPSCs, it should be the case that their spiking will be controlled mostly by the PN population more entrained to the EPSC rhythm, i.e., the strongly driven group. We examined this by measuring the cross-correlation between PNs from each group and all spikes emitted by the FSI network (Fig. 7D1) after correcting for the autocorrelations inherent in both spike trains. This revealed a peak in the cross-correlation at a 4-ms lag (Fig. 7D2, inset), consistent with the previous analysis (Fig. 7C2–4). Casting the normalized cross-correlation function into the frequency domain returned a measure of coherence between PN and FSI spiking. This revealed that the more strongly driven PN group had higher coherence with the FSI population, particularly in the gamma band [group 1 > group 2, rank-sum test, z value = 2.6, p = 0.009 (Fig. 7D2); group 2 > group 1, rank-sum test, z value = −3.7, p = 2.3 × 10−4 (Fig. 7D4)], and there was no difference when both groups were equally driven (group 1 = group 2, z value = −0.3, p = 0.79; Fig. 7D3).

Since the more strongly activated PN group exerts a greater influence over the FSI network, it seemed likely that this would lead to suppression of the weaker PN ensemble. Indeed, we found that as group 2 was more strongly excited, it diminished the firing of group 1 (Fig. 7E1). Presumably, this effect was mediated by the PN→FSI connection, which were also critical for generating the gamma rhythm. Eliminating these connections reduced the difference in firing rates between groups 1 and 2 (Fig. 7E2), suggesting that the FSI network may mediate competition between PN ensembles.

If FSIs allow for competition between ensembles of PNs, and FSIs tend to project locally, one would expect that any competition effect will operate when the groups of PNs overlap spatially, sharing the same FSI network, and not when they are far apart from one another. To explore this possibility, we returned to the original full model with 27,000 neurons. Our spatial analyses of gamma bursts revealed that they spanned an average radius of ∼270 µm (Fig. 6B4), so we created spheroidal subgroups of PNs with a similar size that were driven by Poisson afferents. When those groups overlapped in the model, they exhibited the competition effect (Fig. 7F1). Increasing the afferent drive onto group 2 reduced the firing rate of group 1, whose drive was held constant. However, if they were completely separated then no competition was evident (for group 1 interaction between separation × group 2 afferent drive, ANCOVA, F(1,3474) = 7.8, p = 0.005; Fig. 7F2). Also in agreement with the stronger group 1 exerting downward pressure, for group 2, we found a significant main effect increase in firing rate on spatial separation (F(1,3471) = 120.5, p = 1.4 × 10−27).

The above analyses suggest that the competition between ensembles arose from the heightened gamma-band interaction between a strongly driven group of PNs and the local FSI network. Since the competition was reflected mainly as a suppression of spiking, it was probably mediated by IPSCs impinging on the PN network. Moreover, since PN spikes tended to precede FSI spikes (Fig. 7D2 , inset), the dominant PN group should, on average, control the inhibitory rhythm, which would periodically modulate neuronal responsiveness to excitatory drive, as found in sensory cortices (Cardin et al., 2009; Ni et al., 2016). To examine this, we measured the relationship between EPSC strength and spiking probability during different gamma phases. For the weakly driven group, the trough phase of the gamma cycle showed the strongest relationship between EPSC strength and spiking probability, likely because of the waning inhibition (Fig. 7G, bottom left). Once that inhibition was reinstated during the peak, it shunted EPSCs and degraded their ability to drive spiking. We sought to determine if this dependence was specific to gamma generation. Removing the PN→FSI connections eliminates gamma in the network while sparing inhibitory drive onto PNs. Doing so diminished the gamma phase dependence of sensitivity to EPSCs (Fig. 7G, bottom right).

Discussion

A biophysically and anatomically detailed 27,000-cell model of BL produced transient gamma oscillations with properties similar to those seen in vivo. This concordance extended to gamma band entrainment and preferred phase of PN and FSI spiking. BL gamma arose from reciprocal interactions between PNs and FSIs, which suggests a PING mechanism of rhythmogenesis. Having established the validity of our model, we explored questions that could not be tested experimentally. Underscoring the importance of inhibitory connectivity, both PNs and FSIs were more synchronized to gamma if they received a greater number of inhibitory afferents. Finally, this connectivity produced lateral inhibition, which mediated competition between ensembles of PNs, with the more strongly driven PN population exerting greater control over the FSI network. It is noted that computational models of transient gamma have not been reported for the amygdala, and the ones reported for other brain regions are typically not biophysically or anatomically detailed to the extent that we present here.

Emulation of BL gamma

BL has similar cell types and connectivity as found in cortical circuits (Swanson, 2003). A majority of its neurons are composed of PNs and FSIs (Rainnie et al., 1993, 2006), are reciprocally connected with each other (Muller et al., 2005; Woodruff and Sah, 2007), and like cortex (Braitenberg and Schu¨z, 1991) the BL receives a substantial portion of its input from other cortical regions (Sah et al., 2003). Also as in cortex, examples abound of gamma oscillations recorded from the amygdala (Oya et al., 2002; Bauer et al., 2007; Stujenske et al., 2014), and in vitro slice preparations can be coaxed into generating gamma through pharmacological manipulations (Sinfield and Collins, 2006; Randall et al., 2011) similar to those used in cortical slices (Whittington et al., 2000). Thus, the importance of BL gamma oscillations seems assured, and yet our understanding of how BL produces gamma is limited.

As a first step, we examined a fundamental question: are the local circuits present in that area sufficient to produce the intermittent rhythm measured in vivo? Doing this required collating the anatomic and physiologic properties of BL neurons into a model and driving it with spike trains mirroring those found in upstream areas. Then, we calculated the extracellular field that could be picked up by recording electrodes (Pesaran et al., 2018). On running our simulations, the LFP exhibited a peak in the gamma band that was in the same range of frequencies and amplitudes as observed in vivo. Moreover, the spiking of both PNs and FSIs were entrained to the gamma rhythm as found in vivo (Amir et al., 2018). This concordance suggested that the BL network, as captured by our model, is sufficient for the generation of gamma like that seen in vivo.

However, the model exhibited far less power in the low-frequency bands. Given the concordance between our model of the extracellular field for both APs and activity in the gamma band, we suspect that this discrepancy reflects the extrinsic origins of low-frequency activity. One possible source is volume conduction from adjacent cortical regions. It may also be that afferent synapses from upstream areas contribute to the low-frequency band. Supporting this argument, driving the model with spike trains derived from cortical recordings increased power in the low-frequency band compared with rate matched Poisson inputs. Besides neocortical and paleocortical afferents, BLA neurons receive direct projections from the ventral hippocampus, and in turn respond to hippocampal theta (Bienvenu et al., 2012) and intermittent sharp wave associated population bursts (Girardeau et al., 2017), both of which occupy the low-frequency range.

BL microcircuitry presented unique challenges to gamma rhythmogenesis

The generation of gamma oscillations in our BL model was not a foregone conclusion. While qualitatively the BL network contains the microcircuitry required for producing gamma oscillations, quantitative differences in the properties of this network from previous PING models may have precluded gamma generation.

For a PING rhythm to emerge PNs must fire in advance of FSIs (Whittington et al., 2000), and thus increasing the average firing rates of PNs contributes to the strength of gamma oscillations produced in PING models and in slices (Adesnik and Scanziani, 2010). PN firing rates in the BL are smaller than those in the cortical models used to simulate gamma by at least a factor of 6 (Börgers and Kopell, 2008; Hoseini and Wessel, 2016; Chariker et al., 2018), which may have diminished the network’s ability to produce a PING type gamma rhythm. Also important is the strength of inhibitory feedback in the network (Börgers et al., 2012). Two factors that contribute to this are the strength of connections between PNs and FSIs and the number of FSIs in the network. We and others have found that the BL has a lower proportion of interneurons (McDonald and Mascagni, 2001) compared with standard PING models (Börgers and Kopell, 2008; Hoseini and Wessel, 2016; Chariker et al., 2018). Moreover, the connection probability between PNs and FSIs in BL is lower (Woodruff and Sah, 2007) than the aforementioned models. Thus, a priori the BL network may not be capable of producing gamma, and so our model was essential to at least test the possibility theoretically.

Another difference between our BL model and cortical-based PING models is a lower proportion of recurrent connectivity among PNs. Brunel and Wang (2003) found that including recurrent connections between excitatory units in a PING model could downshift the frequency of the network oscillation by half. On the other hand, our BL PING model has fairly little dependence on these connections, since their elimination did not substantially affect the prominence or peak frequency of gamma (Fig. 4D). In addition, other models of gamma generation depend on the recurrent excitation among excitatory units. For instance, the REI model (Chariker et al., 2018) contains a network of reciprocally connected excitatory and inhibitory units. Unlike PING models, the gamma-periodic bursts of activity depend on the buildup of recurrent excitation in the PN network, which raises the membrane potential for all units in the network, and drives them to fire synchronously. This mechanism for producing gamma is distinctly not PING, since the excitatory and inhibitory units activate simultaneously, instead of with the inhibitory units following the excitatory. Moreover, since BL PNs and FSIs in vivo fire during different phases of the gamma cycle (Amir et al., 2018), it is likely that the REI mechanism does not operate within the BL.

Thus, our model establishes the sufficiency of the BL network to produce gamma oscillations, and for those to arise via a PING mechanism. Despite that, it was also possible that any gamma produced by the model would be of insufficient power to comprise what is observed in the LFP in vivo. To address this, our model included a realistic estimation of the LFP, which most PING models lack (Brunel and Hakim, 1999; Traub et al., 2000; Brunel and Wang, 2003; Börgers et al., 2005; Bathellier et al., 2006; Neymotin et al., 2011; Economo and White, 2012; Chalk et al., 2016; Hoseini and Wessel, 2016; Palmigiano et al., 2017; Chariker et al., 2018). Since our model LFP exhibited gamma oscillations that were close to the frequency and power of those recorded in vivo, it establishes the sufficiency of the local microcircuitry in BL to explain the gamma observed in LFP recordings in vivo.

Circuit contributions to the properties of gamma oscillations

Cortical gamma oscillations are thought to arise from an interaction between PNs and FSIs, whereby PNs strongly drive a recurrently connected network of FSIs, that in turn inhibit PNs for a short period of time, after which PNs are able to restart the cycle over again. We determined whether a similar mechanism was operating in our BL simulations, and investigated qualitative and quantitative features of the underlying microcircuits. Systematically removing each class of connections revealed that either the connections from PNs to FSIs or vice versa were crucial for generating the gamma rhythm, which is in agreement with a PING type mechanism. Although interactions between FSIs are sufficient to produce a gamma rhythm (Whittington et al., 2000), this was not the case in our model. Instead, they affected the peak frequency. This could be functionally relevant because presynaptic receptors specific to synapses between interneurons (Cossart et al., 2001) may be able to regulate the frequency of gamma rhythms.

The dynamics of EPSCs and IPSCs in our model (Fig. 5B,D) mostly aligned with those from a recent in vitro model of cortical gamma (Salkoff et al., 2015). Using paired whole-cell patch clamp recordings, synaptic currents were recorded in one cell, and correlated with the spiking of an adjacent one. The power spectrum of IPSCs in PNs and FSIs exhibited a prominence in the low gamma band, while EPSCs in PNs did not. Moreover, they observed that IPSCs weakened before PN and FSI spiking, after which they rebounded. We observed that pattern for PNs, while for FSIs the IPSCs returned to baseline after spiking. Another significant difference was that they only observed substantial increases in EPSC amplitude surrounding spikes from FSIs, but not PNs. This likely resulted from the fact that whereas FSIs receive convergent inputs from many PNs (Oswald et al., 2009), excitatory synapses are sparsely shared between PNs; in all likelihood the afferents that elicited spiking in one PN did not collateralize onto another. In contrast, we could monitor both the EPSC and spiking in the same cell and found a robust EPSC increase during PN spiking.

One unknown aspect of gamma generation that cannot be readily explored at present either in vivo or in vitro is the contribution of a neuron’s particular complement of excitatory and inhibitory afferents to its entrainment by the gamma rhythm. In addressing this, we found that the number of inhibitory afferents was particularly important for how a neuron would be affected by gamma. PNs and FSIs with more inhibitory contacts were more strongly entrained to gamma, while only FSIs showed a significant reduction in entrainment as the number of excitatory afferents increased. In agreement with these results, a previous study from our group (Amir et al., 2018) found that PNs that received a monosynaptic input from a nearby FSI showed stronger entrainment to gamma. However, our model differed from this study regarding the contribution of excitatory connections onto FSIs. In Amir et al. (2018), the presence of an excitatory monosynaptic contact onto an FSI was associated with increased entrainment, whereas in our model, it was not. Perhaps resolving this disagreement, we found that a concomitant increase in both excitatory and inhibitory contacts onto an FSI raised entrainment beyond what was expected by either connection type alone. Thus, our two studies could be reconciled if BL FSIs that receive more afferents from PNs also receive more from other FSIs. Perhaps this suggests that inhibitory and excitatory afferents need to be tuned for single neurons to be strongly entrained by the gamma rhythm. This may be reflected by experience-dependent plasticity of inhibitory networks (Donato et al., 2013).

Computations emerging from gamma generating circuits in BL

In cortical networks, several functions have been ascribed to gamma, but a particular emphasis has been placed on their ability to synchronize spiking (Tiesinga et al., 2008) and mediate competition between cell assemblies (Börgers et al., 2005; de Almeida et al., 2009; Palmigiano et al., 2017). Our model exhibited both of these phenomena. PNs synchronized their spiking by phase-locking to the gamma rhythm. This could be crucial to BL function since its neurons fire at very low rates, and so it is likely they need to coordinate their activity to drive robust postsynaptic integration. We also found that activating a subset of PNs in the network drove a suppression of those that received a weaker input, and that gamma genesis played a critical role in this effect. This competition may be important in the BL for the selective recruitment of particular PN ensembles, allowing it to drive its many downstream targets with greater specificity.

Conclusions

Gamma oscillations are a feature of cortical processing that is also expressed in BL, where they are likely generated via similar mechanisms, and perhaps perform the same functions. Moving forward, there are two reasons why having a detailed biophysical model of gamma generation in BL is important. First, we need a better understanding of the variety of situations that bring about and alter gamma oscillations, and a concrete model provides a framework in which these observations can be integrated with one another. Second, a biophysically detailed model can serve as a testbed for exploring the means and algorithms that might be used to alter these oscillations (Witt et al., 2013), and the functions of the BL itself.

Supplementary Material 1

Supplementary Code used in article. Download Supplementary Material 1, ZIP file

Footnotes

  • The authors declare no competing financial interests.

  • This work was supported by Department of Health and Human Services | National Institutes of Health Grants MH087755 and MH109122 (to S.S.N.) and MH106403 to (D.P.).

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

  1. ↵
    Abatis M, Niu R, Perin R, Giobellina G, Bito H, Markram H, Stoop R (2017) The role of local excitatory networks in the lateral amygdala in emotional memory learning. Poster presented in Society for Neuroscience. Washington, DC.
  2. ↵
    Adesnik H, Scanziani M (2010) Lateral competition for cortical space by layer-specific horizontal circuits. Nature 464:1155–1160. doi:10.1038/nature08935 pmid:20414303
    OpenUrlCrossRefPubMed
  3. ↵
    Alturki A, Feng F, Nair A, Guntu V, Nair SS (2016) Distinct current modules shape cellular dynamics in model neurons. Neuroscience 334:309–331. doi:10.1016/j.neuroscience.2016.08.016 pmid:27530698
    OpenUrlCrossRefPubMed
  4. ↵
    Amir A, Headley DB, Lee SC, Haufler D, Paré D (2018) Vigilance-associated gamma oscillations coordinate the ensemble activity of basolateral amygdala neurons. Neuron 97:656–669.e7. doi:10.1016/j.neuron.2017.12.035 pmid:29420934
    OpenUrlCrossRefPubMed
  5. ↵
    Atallah BV, Scanziani M (2009) Instantaneous modulation of gamma oscillation frequency by balancing excitation with inhibition. Neuron 62:566–577. doi:10.1016/j.neuron.2009.04.027 pmid:19477157
    OpenUrlCrossRefPubMed
  6. ↵
    Aydore S, Pantazis D, Leahy RM (2013) A note on the phase locking value and its properties. Neuroimage 74:231–244. doi:10.1016/j.neuroimage.2013.02.008 pmid:23435210
    OpenUrlCrossRefPubMed
  7. ↵
    Bartos M, Vida I, Jonas P (2007) Synaptic mechanisms of synchronized gamma oscillations in inhibitory interneuron networks. Nat Rev Neurosci 8:45–56. doi:10.1038/nrn2044 pmid:17180162
    OpenUrlCrossRefPubMed
  8. ↵
    Bathellier B, Lagier S, Faure P, Lledo PM (2006) Circuit properties generating gamma oscillations in a network model of the olfactory bulb. J Neurophysiol 95:2678–2691. doi:10.1152/jn.01141.2005 pmid:16381804
    OpenUrlCrossRefPubMed
  9. ↵
    Bauer EP, Paz R, Paré D (2007) Gamma oscillations coordinate amygdalo-rhinal interactions during learning. J Neurosci 27:9369–9379. doi:10.1523/JNEUROSCI.2153-07.2007 pmid:17728450
    OpenUrlAbstract/FREE Full Text
  10. ↵
    Berens P (2009) CircStat: a Matlab toolbox for circular statistics. J Stat Softw 31:1–21. doi:10.18637/jss.v031.i10
    OpenUrlCrossRefPubMed
  11. ↵
    Beyeler A, Namburi P, Glober GF, Simonnet C, Calhoon GG, Conyers GF, Luck R, Wildes CP, Tye KM (2016) Divergent routing of positive and negative information from the amygdala during memory retrieval. Neuron 90:348–361. doi:10.1016/j.neuron.2016.03.004 pmid:27041499
    OpenUrlCrossRefPubMed
  12. ↵
    Bezaire MJ, Raikov I, Burk K, Vyas D, Soltesz I (2016) Interneuronal mechanisms of hippocampal theta oscillations in a full-scale model of the rodent CA1 circuit. Elife 5. doi:10.7554/eLife.18566 pmid:28009257
    OpenUrlCrossRefPubMed
  13. ↵
    Bienvenu TC, Busti D, Magill PJ, Ferraguti F, Capogna M (2012) Cell-type-specific recruitment of amygdala interneurons to hippocampal theta rhythm and noxious stimuli in vivo . Neuron 74:1059–1074. doi:10.1016/j.neuron.2012.04.022 pmid:22726836
    OpenUrlCrossRefPubMed
  14. ↵
    Börgers C, Kopell N (2003) Synchronization in networks of excitatory and inhibitory neurons with sparse, random connectivity. Neural Comput 15:509–538. doi:10.1162/089976603321192059 pmid:12620157
    OpenUrlCrossRefPubMed
  15. ↵
    Börgers C, Kopell NJ (2008) Gamma oscillations and stimulus selection. Neural Comput 20:383–414. doi:10.1162/neco.2007.07-06-289 pmid:18047409
    OpenUrlCrossRefPubMed
  16. ↵
    Börgers C, Epstein S, Kopell NJ (2005) Background gamma rhythmicity and attention in cortical local circuits: a computational study. Proc Natl Acad Sci USA 102:7002–7007. doi:10.1073/pnas.0502366102 pmid:15870189
    OpenUrlAbstract/FREE Full Text
  17. ↵
    Börgers C, Epstein S, Kopell NJ (2008) Gamma oscillations mediate stimulus competition and attentional selection in a cortical network model. Proc Natl Acad Sci USA 105:18023–18028. doi:10.1073/pnas.0809511105 pmid:19004759
    OpenUrlAbstract/FREE Full Text
  18. ↵
    Börgers C, Talei Franzesi G, Lebeau FE, Boyden ES, Kopell NJ (2012) Minimal size of cell assemblies coordinated by gamma oscillations. PLoS Comput Biol 8:e1002362. doi:10.1371/journal.pcbi.1002362 pmid:22346741
    OpenUrlCrossRefPubMed
  19. ↵
    Braitenberg V, Schu¨Z A (1991) Anatomy of the cortex: statistics and geometry. Berlin; New York: Springer.
  20. ↵
    Brunel N, Hakim V (1999) Fast global oscillations in networks of integrate-and-fire neurons with low firing rates. Neural Comput 11:1621–1671. doi:10.1162/089976699300016179 pmid:10490941
    OpenUrlCrossRefPubMed
  21. ↵
    Brunel N, Wang XJ (2003) What determines the frequency of fast network oscillations with irregular neural discharges? I. Synaptic dynamics and excitation-inhibition balance. J Neurophysiol 90:415–430. doi:10.1152/jn.01095.2002 pmid:12611969
    OpenUrlCrossRefPubMed
  22. ↵
    Buzsáki G, Anastassiou CA, Koch C (2012) The origin of extracellular fields and currents–EEG, ECoG, LFP and spikes. Nat Rev Neurosci 13:407–420. doi:10.1038/nrn3241 pmid:22595786
    OpenUrlCrossRefPubMed
  23. ↵
    Byrne JH, Roberts JL (2004) From molecules to networks: an introduction to cellular and molecular neuroscience. Amsterdam; Boston: Elsevier Academic Press.
  24. Cammarota M, Losi G, Chiavegato A, Zonta M, Carmignoto G (2013) Fast spiking interneuron control of seizure propagation in a cortical slice model of focal epilepsy. J Physiol 591:807–822. doi:10.1113/jphysiol.2012.238154 pmid:23207591
    OpenUrlCrossRefPubMed
  25. ↵
    Cardin JA, Carlén M, Meletis K, Knoblich U, Zhang F, Deisseroth K, Tsai LH, Moore CI (2009) Driving fast-spiking cells induces gamma rhythm and controls sensory responses. Nature 459:663–667. doi:10.1038/nature08002 pmid:19396156
    OpenUrlCrossRefPubMed
  26. ↵
    Carnevale NT, Hines ML (2005) The NEURON book. New York, NY: Cambridge University Press.
  27. ↵
    Chalk M, Gutkin B, Deneve S (2016) Neural oscillations as a signature of efficient coding in the presence of synaptic delays. Elife 5: pii: e13824. doi:10.7554/elife.13824 pmid:27383272
    OpenUrlCrossRefPubMed
  28. ↵
    Chareyron LJ, Banta Lavenex P, Amaral DG, Lavenex P (2011) Stereological analysis of the rat and monkey amygdala. J Comp Neurol 519:3218–3239. doi:10.1002/cne.22677 pmid:21618234
    OpenUrlCrossRefPubMed
  29. ↵
    Chariker L, Shapley R, Young LS (2018) Rhythm and synchrony in a cortical network model. J Neurosci 38:8621–8634. doi:10.1523/JNEUROSCI.0675-18.2018 pmid:30120205
    OpenUrlAbstract/FREE Full Text
  30. ↵
    Cossart R, Tyzio R, Dinocourt C, Esclapez M, Hirsch JC, Ben-Ari Y, Bernard C (2001) Presynaptic kainate receptors that enhance the release of GABA on CA1 hippocampal interneurons. Neuron 29:497–508. doi:10.1016/S0896-6273(01)00221-5 pmid:11239438
    OpenUrlCrossRefPubMed
  31. ↵
    de Almeida L, Idiart M, Lisman JE (2009) A second function of gamma frequency oscillations: an E%-max winner-take-all mechanism selects which cells fire. J Neurosci 29:7497–7503. doi:10.1523/JNEUROSCI.6044-08.2009 pmid:19515917
    OpenUrlAbstract/FREE Full Text
  32. ↵
    Destexhe A, Mainen ZF, Sejnowski TJ (1994) An efficient method for computing synaptic conductances based on a kinetic model of receptor binding. Neural Comput 6:14–18. doi:10.1162/neco.1994.6.1.14
    OpenUrlCrossRef
  33. ↵
    Destexhe A, Rudolph M, Fellous JM, Sejnowski TJ (2001) Fluctuating synaptic conductances recreate in vivo-like activity in neocortical neurons. Neuroscience 107:13–24. doi:10.1016/s0306-4522(01)00344-x pmid:11744242
    OpenUrlCrossRefPubMed
  34. ↵
    Disterhoft JF, Oh MM (2007) Alterations in intrinsic neuronal excitability during normal aging. Aging Cell 6:327–336. doi:10.1111/j.1474-9726.2007.00297.x pmid:17517042
    OpenUrlCrossRefPubMed
  35. ↵
    Donato F, Rompani SB, Caroni P (2013) Parvalbumin-expressing basket-cell network plasticity induced by experience regulates adult learning. Nature 504:272–276. doi:10.1038/nature12866 pmid:24336286
    OpenUrlCrossRefPubMed
  36. ↵
    Durstewitz D, Seamans JK, Sejnowski TJ (2000) Dopamine-mediated stabilization of delay-period activity in a network model of prefrontal cortex. J Neurophysiol 83:1733–1750. doi:10.1152/jn.2000.83.3.1733 pmid:10712493
    OpenUrlCrossRefPubMed
  37. ↵
    Duvarci S, Pare D (2014) Amygdala microcircuits controlling learned fear. Neuron 82:966–980. doi:10.1016/j.neuron.2014.04.042 pmid:24908482
    OpenUrlCrossRefPubMed
  38. ↵
    Dvorak D, Fenton AA (2014) Toward a proper estimation of phase-amplitude coupling in neural oscillations. J Neurosci Methods 225:42–56. doi:10.1016/j.jneumeth.2014.01.002 pmid:24447842
    OpenUrlCrossRefPubMed
  39. ↵
    Economo MN, White JA (2012) Membrane properties and the balance between excitation and inhibition control gamma-frequency oscillations arising from feedback inhibition. PLoS Comput Biol 8:e1002354. doi:10.1371/journal.pcbi.1002354 pmid:22275859
    OpenUrlCrossRefPubMed
  40. ↵
    Ehrlich I, Humeau Y, Grenier F, Ciocchi S, Herry C, Lüthi A (2009) Amygdala inhibitory circuits and the control of fear memory. Neuron 62:757–771. doi:10.1016/j.neuron.2009.05.026 pmid:19555645
    OpenUrlCrossRefPubMed
  41. ↵
    Einevoll GT, Kayser C, Logothetis NK, Panzeri S (2013) Modelling and analysis of local field potentials for studying the function of cortical circuits. Nat Rev Neurosci 14:770–785. doi:10.1038/nrn3599 pmid:24135696
    OpenUrlCrossRefPubMed
  42. ↵
    Feng F, Samarth P, Paré D, Nair SS (2016) Mechanisms underlying the formation of the amygdalar fear memory trace: a computational perspective. Neuroscience 322:370–376. doi:10.1016/j.neuroscience.2016.02.059 pmid:26944604
    OpenUrlCrossRefPubMed
  43. ↵
    Fishell G, Rudy B (2011) Mechanisms of inhibition within the telencephalon: "where the wild things are." Annu Rev Neurosci 34:535–567. doi:10.1146/annurev-neuro-061010-113717 pmid:21469958
    OpenUrlCrossRefPubMed
  44. Galarreta M, Hestrin S (1997) Properties of GABAA receptors underlying inhibitory synaptic currents in neocortical pyramidal neurons. J Neurosci 17:7220–7227. doi:10.1523/jneurosci.17-19-07220.1997 pmid:9295368
    OpenUrlAbstract/FREE Full Text
  45. ↵
    Girardeau G, Inema I, Buzsáki G (2017) Reactivations of emotional memory in the hippocampus-amygdala system during sleep. Nat Neurosci 20:1634–1642. doi:10.1038/nn.4637 pmid:28892057
    OpenUrlCrossRefPubMed
  46. ↵
    Gold C, Henze DA, Koch C, Buzsáki G (2006) On the origin of the extracellular action potential waveform: a modeling study. J Neurophysiol 95:3113–3128. doi:10.1152/jn.00979.2005 pmid:16467426
    OpenUrlCrossRefPubMed
  47. ↵
    Goto T, Hatanaka R, Ogawa T, Sumiyoshi A, Riera J, Kawashima R (2010) An evaluation of the conductivity profile in the somatosensory barrel cortex of Wistar rats. J Neurophysiol 104:3388–3412. doi:10.1152/jn.00122.2010 pmid:20810682
    OpenUrlCrossRefPubMed
  48. Guzman SJ, Schlögl A, Frotscher M, Jonas P (2016) Synaptic mechanisms of pattern completion in the hippocampal CA3 network. Science 353:1117–1123. doi:10.1126/science.aaf1836 pmid:27609885
    OpenUrlAbstract/FREE Full Text
  49. ↵
    Hagen E, Dahmen D, Stavrinou ML, Lindén H, Tetzlaff T, van Albada SJ, Grün S, Diesmann M, Einevoll GT (2016) Hybrid scheme for modeling local field potentials from point-neuron networks. Cereb Cortex 26:4461–4496. doi:10.1093/cercor/bhw237 pmid:27797828
    OpenUrlCrossRefPubMed
  50. ↵
    Headley DB, DeLucca MV, Haufler D, Paré D (2015) Incorporating 3D-printing technology in the design of head-caps and electrode drives for recording neurons in multiple brain regions. J Neurophysiol 113:2721–2732. doi:10.1152/jn.00955.2014 pmid:25652930
    OpenUrlCrossRefPubMed
  51. ↵
    Holt GR (1998) A critical examination of some assumptions and implications of cable theory in neurobiology. PhD Thesis, California Institute of Technology.
  52. ↵
    Hoseini MS, Wessel R (2016) Coherent and intermittent ensemble oscillations emerge from networks of irregular spiking neurons. J Neurophysiol 115:457–469. doi:10.1152/jn.00578.2015 pmid:26561602
    OpenUrlCrossRefPubMed
  53. ↵
    Hübner C, Bosch D, Gall A, Lüthi A, Ehrlich I (2014) Ex vivo dissection of optogenetically activated mPFC and hippocampal inputs to neurons in the basolateral amygdala: implications for fear and emotional memory. Front Behav Neurosci 8:64. doi:10.3389/fnbeh.2014.00064 pmid:24634648
    OpenUrlCrossRefPubMed
  54. ↵
    Hummos A, Franklin CC, Nair SS (2014) Intrinsic mechanisms stabilize encoding and retrieval circuits differentially in a hippocampal network model. Hippocampus 24:1430–1448. doi:10.1002/hipo.22324 pmid:24978936
    OpenUrlCrossRefPubMed
  55. ↵
    Kim D, Paré D, Nair SS (2013) Mechanisms contributing to the induction and storage of Pavlovian fear memories in the lateral amygdala. Learn Mem 20:421–430. doi:10.1101/lm.030262.113 pmid:23864645
    OpenUrlAbstract/FREE Full Text
  56. ↵
    LeDoux JE (2000) Emotion circuits in the brain. Annu Rev Neurosci 23:155–184. doi:10.1146/annurev.neuro.23.1.155 pmid:10845062
    OpenUrlCrossRefPubMed
  57. ↵
    Lee S, Jones SR (2013) Distinguishing mechanisms of gamma frequency oscillations in human current source signals using a computational model of a laminar neocortical network. Front Hum Neurosci 7:869. doi:10.3389/fnhum.2013.00869 pmid:24385958
    OpenUrlCrossRefPubMed
  58. ↵
    Łęski S, Lindén H, Tetzlaff T, Pettersen KH, Einevoll GT (2013) Frequency dependence of signal power and spatial reach of the local field potential. PLoS Comput Biol 9:e1003137. doi:10.1371/journal.pcbi.1003137 pmid:23874180
    OpenUrlCrossRefPubMed
  59. ↵
    Li G, Nair SS, Quirk GJ (2009) A biologically realistic network model of acquisition and extinction of conditioned fear associations in lateral amygdala neurons. J Neurophysiol 101:1629–1646. doi:10.1152/jn.90765.2008 pmid:19036872
    OpenUrlCrossRefPubMed
  60. ↵
    Lindén H, Tetzlaff T, Potjans TC, Pettersen KH, Grün S, Diesmann M, Einevoll GT (2011) Modeling the spatial reach of the LFP. Neuron 72:859–872. doi:10.1016/j.neuron.2011.11.006 pmid:22153380
    OpenUrlCrossRefPubMed
  61. ↵
    Lindén H, Hagen E, Łęski S, Norheim ES, Pettersen KH, Einevoll GT (2013) LFPy: a tool for biophysical simulation of extracellular potentials generated by detailed model neurons. Front Neuroinform 7:41. doi:10.3389/fninf.2013.00041 pmid:24474916
    OpenUrlCrossRefPubMed
  62. Mahanty NK, Sah P (1998) Calcium-permeable AMPA receptors mediate long-term potentiation in interneurons in the amygdala. Nature 394:683–687. doi:10.1038/29312 pmid:9716132
    OpenUrlCrossRefPubMed
  63. ↵
    Martina M, Royer S, Paré D (2001) Cell-type-specific GABA responses and chloride homeostasis in the cortex and amygdala. J Neurophysiol 86:2887–2895. doi:10.1152/jn.2001.86.6.2887 pmid:11731545
    OpenUrlCrossRefPubMed
  64. ↵
    McDonald AJ (1984) Neuronal organization of the lateral and basolateral amygdaloid nuclei in the rat. J Comp Neurol 222:589–606. doi:10.1002/cne.902220410
    OpenUrlCrossRefPubMed
  65. ↵
    McDonald AJ (1998) Cortical pathways to the mammalian amygdala. Prog Neurobiol 55:257–332. doi:10.1016/s0301-0082(98)00003-3 pmid:9643556
    OpenUrlCrossRefPubMed
  66. ↵
    McDonald AJ, Mascagni F (1997) Projections of the lateral entorhinal cortex to the amygdala: a Phaseolus vulgaris leucoagglutinin study in the rat. Neuroscience 77:445–459. doi:10.1016/s0306-4522(96)00478-2 pmid:9472403
    OpenUrlCrossRefPubMed
  67. ↵
    McDonald AJ, Mascagni F (2001) Colocalization of calcium-binding proteins and GABA in neurons of the rat basolateral amygdala. Neuroscience 105:681–693. doi:10.1016/s0306-4522(01)00214-7 pmid:11516833
    OpenUrlCrossRefPubMed
  68. ↵
    Muller JF, Mascagni F, McDonald AJ (2005) Coupled networks of parvalbumin-immunoreactive interneurons in the rat basolateral amygdala. J Neurosci 25:7366–7376. doi:10.1523/JNEUROSCI.0899-05.2005 pmid:16093387
    OpenUrlAbstract/FREE Full Text
  69. ↵
    Muller JF, Mascagni F, McDonald AJ (2006) Pyramidal cells of the rat basolateral amygdala: synaptology and innervation by parvalbumin-immunoreactive interneurons. J Comp Neurol 494:635–650. doi:10.1002/cne.20832 pmid:16374802
    OpenUrlCrossRefPubMed
  70. ↵
    Neske GT, Connors BW (2016) Synchronized gamma-frequency inhibition in neocortex depends on excitatory-inhibitory interactions but not electrical synapses. J Neurophysiol 116:351–368. doi:10.1152/jn.00071.2016 pmid:27121576
    OpenUrlCrossRefPubMed
  71. ↵
    Neymotin SA, Lazarewicz MT, Sherif M, Contreras D, Finkel LH, Lytton WW (2011) Ketamine disrupts θ modulation of γ in a computer model of hippocampus. J Neurosci 31:11733–11743. doi:10.1523/JNEUROSCI.0501-11.2011 pmid:21832203
    OpenUrlAbstract/FREE Full Text
  72. ↵
    Ni J, Wunderle T, Lewis CM, Desimone R, Diester I, Fries P (2016) Gamma-rhythmic gain modulation. Neuron 92:240–251. doi:10.1016/j.neuron.2016.09.003 pmid:27667008
    OpenUrlCrossRefPubMed
  73. ↵
    Oren I, Mann EO, Paulsen O, Hájos N (2006) Synaptic currents in anatomically identified CA3 neurons during hippocampal gamma oscillations in vitro . J Neurosci 26:9923–9934. doi:10.1523/JNEUROSCI.1580-06.2006 pmid:17005856
    OpenUrlAbstract/FREE Full Text
  74. ↵
    Oswald AM, Doiron B, Rinzel J, Reyes AD (2009) Spatial profile and differential recruitment of GABAB modulate oscillatory activity in auditory cortex. J Neurosci 29:10321–10334. doi:10.1523/JNEUROSCI.1703-09.2009 pmid:19692606
    OpenUrlAbstract/FREE Full Text
  75. ↵
    Oya H, Kawasaki H, Howard MA 3rd., Adolphs R (2002) Electrophysiological responses in the human amygdala discriminate emotion categories of complex visual stimuli. J Neurosci 22:9502–9512. doi:10.1523/JNEUROSCI.22-21-09502.2002 pmid:12417674
    OpenUrlAbstract/FREE Full Text
  76. ↵
    Palmigiano A, Geisel T, Wolf F, Battaglia D (2017) Flexible information routing by transient synchrony. Nat Neurosci 20:1014–1022. doi:10.1038/nn.4569 pmid:28530664
    OpenUrlCrossRefPubMed
  77. ↵
    Pape HC, Pare D (2010) Plastic synaptic networks of the amygdala for the acquisition, expression, and extinction of conditioned fear. Physiol Rev 90:419–463. doi:10.1152/physrev.00037.2009 pmid:20393190
    OpenUrlCrossRefPubMed
  78. ↵
    Pape HC, Paré D, Driesang RB (1998) Two types of intrinsic oscillations in neurons of the lateral and basolateral nuclei of the amygdala. J Neurophysiol 79:205–216. doi:10.1152/jn.1998.79.1.205 pmid:9425192
    OpenUrlCrossRefPubMed
  79. ↵
    Parasuram H, Nair B, D'Angelo E, Hines M, Naldi G, Diwakar S (2016) Computational modeling of single neuron extracellular electric potentials and network local field potentials using LFPsim. Front Comput Neurosci 10:65. doi:10.3389/fncom.2016.00065 pmid:27445781
    OpenUrlCrossRefPubMed
  80. ↵
    Paré D, Smith Y, Paré JF (1995) Intra-amygdaloid projections of the basolateral and basomedial nuclei in the cat: Phaseolus vulgaris-leucoagglutinin anterograde tracing at the light and electron microscopic level. Neuroscience 69:567–583. doi:10.1016/0306-4522(95)00272-K pmid:8552250
    OpenUrlCrossRefPubMed
  81. ↵
    Pêgo JM, Morgado P, Pinto LG, Cerqueira JJ, Almeida OF, Sousa N (2008) Dissociation of the morphological correlates of stress-induced anxiety and fear. Eur J Neurosci 27:1503–1516. doi:10.1111/j.1460-9568.2008.06112.x pmid:18336570
    OpenUrlCrossRefPubMed
  82. ↵
    Penttonen M, Kamondi A, Acsády L, Buzsáki G (1998) Gamma frequency oscillation in the hippocampus of the rat: intracellular analysis in vivo . Eur J Neurosci 10:718–728. doi:10.1046/j.1460-9568.1998.00096.x pmid:9749733
    OpenUrlCrossRefPubMed
  83. ↵
    Pesaran B, Vinck M, Einevoll GT, Sirota A, Fries P, Siegel M, Truccolo W, Schroeder CE, Srinivasan R (2018) Investigating large-scale brain dynamics using field potential recordings: analysis and interpretation. Nat Neurosci 21:903–919. doi:10.1038/s41593-018-0171-8 pmid:29942039
    OpenUrlCrossRefPubMed
  84. ↵
    Power JM, Bocklisch C, Curby P, Sah P (2011) Location and function of the slow afterhyperpolarization channels in the basolateral amygdala. J Neurosci 31:526–537. doi:10.1523/JNEUROSCI.1045-10.2011 pmid:21228162
    OpenUrlAbstract/FREE Full Text
  85. ↵
    Rainnie DG, Asprodini EK, Shinnick-Gallagher P (1993) Intracellular recordings from morphologically identified neurons of the basolateral amygdala. J Neurophysiol 69:1350–1362. doi:10.1152/jn.1993.69.4.1350 pmid:8492168
    OpenUrlCrossRefPubMed
  86. ↵
    Rainnie DG, Mania I, Mascagni F, McDonald AJ (2006) Physiological and morphological characterization of parvalbumin-containing interneurons of the rat basolateral amygdala. J Comp Neurol 498:142–161. doi:10.1002/cne.21049 pmid:16856165
    OpenUrlCrossRefPubMed
  87. ↵
    Randall FE, Whittington MA, Cunningham MO (2011) Fast oscillatory activity induced by kainate receptor activation in the rat basolateral amygdala in vitro . Eur J Neurosci 33:914–922. doi:10.1111/j.1460-9568.2010.07582.x pmid:21255131
    OpenUrlCrossRefPubMed
  88. ↵
    Rubinow MJ, Juraska JM (2009) Neuron and glia numbers in the basolateral nucleus of the amygdala from preweaning through old age in male and female rats: a stereological study. J Comp Neurol 512:717–725. doi:10.1002/cne.21924 pmid:19065620
    OpenUrlCrossRefPubMed
  89. ↵
    Sah P, Faber ES, Lopez De Armentia M, Power J (2003) The amygdaloid complex: anatomy and physiology. Physiol Rev 83:803–834. doi:10.1152/physrev.00002.2003 pmid:12843409
    OpenUrlCrossRefPubMed
  90. ↵
    Salinas E, Sejnowski TJ (2000) Impact of correlated synaptic input on output firing rate and variability in simple neuronal models. J Neurosci 20:6193–6209. doi:10.1523/jneurosci.20-16-06193.2000 pmid:10934269
    OpenUrlAbstract/FREE Full Text
  91. ↵
    Salinas E, Sejnowski TJ (2002) Integrate-and-fire neurons driven by correlated stochastic input. Neural Comput 14:2111–2155. doi:10.1162/089976602320264024 pmid:12184845
    OpenUrlCrossRefPubMed
  92. ↵
    Salkoff DB, Zagha E, Yüzgeç Ö, McCormick DA (2015) Synaptic mechanisms of tight spike synchrony at gamma frequency in cerebral cortex. J Neurosci 35:10236–10251. doi:10.1523/JNEUROSCI.0828-15.2015 pmid:26180200
    OpenUrlAbstract/FREE Full Text
  93. ↵
    Salm AK, Pavelko M, Krouse EM, Webster W, Kraszpulski M, Birkle DL (2004) Lateral amygdaloid nucleus expansion in adult rats is associated with exposure to prenatal stress. Brain Res Dev Brain Res 148:159–167. doi:10.1016/j.devbrainres.2003.11.005 pmid:14766193
    OpenUrlCrossRefPubMed
  94. ↵
    Samson RD, Paré D (2006) A spatially structured network of inhibitory and excitatory connections directs impulse traffic within the lateral amygdala. Neuroscience 141:1599–1609. doi:10.1016/j.neuroscience.2006.04.077 pmid:16753264
    OpenUrlCrossRefPubMed
  95. ↵
    Schomburg EW, Anastassiou CA, Buzsáki G, Koch C (2012) The spiking component of oscillatory extracellular potentials in the rat hippocampus. J Neurosci 32:11798–11811. doi:10.1523/JNEUROSCI.0656-12.2012 pmid:22915121
    OpenUrlAbstract/FREE Full Text
  96. ↵
    Shi CJ, Cassell MD (1999) Perirhinal cortex projections to the amygdaloid complex and hippocampal formation in the rat. J Comp Neurol 406:299–328. doi:10.1002/(SICI)1096-9861(19990412)406:3<299::AID-CNE2>3.0.CO;2-9 pmid:10102498
    OpenUrlCrossRefPubMed
  97. ↵
    Silberberg G, Wu C, Markram H (2004) Synaptic dynamics control the timing of neuronal excitation in the activated neocortical microcircuit. J Physiol 556:19–27. doi:10.1113/jphysiol.2004.060962 pmid:14978208
    OpenUrlCrossRefPubMed
  98. ↵
    Sinfield JL, Collins DR (2006) Induction of synchronous oscillatory activity in the rat lateral amygdala in vitro is dependent on gap junction activity. Eur J Neurosci 24:3091–3095. doi:10.1111/j.1460-9568.2006.05202.x pmid:17156370
    OpenUrlCrossRefPubMed
  99. ↵
    Sivagnanam S, Majumdar A, Yoshimoto K, Astakhov V, Bandrowski A, Martone ME, Carnevale N (2013) Introducing the Neuroscience Gateway, IWSG, volume 993 of CEUR Workshop Proceedings, CEUR-WS.org, Zurich, Switzerland, 3-5 June, 2013.
  100. ↵
    Smith Y, Paré JF, Paré D (1998) Cat intraamygdaloid inhibitory network: ultrastructural organization of parvalbumin-immunoreactive elements. J Comp Neurol 391:164–179. doi:10.1002/(sici)1096-9861(19980209)391:2<164::aid-cne2>3.0.co;2-0 pmid:9518267
    OpenUrlCrossRefPubMed
  101. ↵
    Sohal VS, Zhang F, Yizhar O, Deisseroth K (2009) Parvalbumin neurons and gamma rhythms enhance cortical circuit performance. Nature 459:698–702. doi:10.1038/nature07991 pmid:19396159
    OpenUrlCrossRefPubMed
  102. ↵
    Spampanato J, Polepalli J, Sah P (2011) Interneurons in the basolateral amygdala. Neuropharmacology 60:765–773. doi:10.1016/j.neuropharm.2010.11.006 pmid:21093462
    OpenUrlCrossRefPubMed
  103. ↵
    Stujenske JM, Likhtik E, Topiwala MA, Gordon JA (2014) Fear and safety engage competing patterns of theta-gamma coupling in the basolateral amygdala. Neuron 83:919–933. doi:10.1016/j.neuron.2014.07.026 pmid:25144877
    OpenUrlCrossRefPubMed
  104. ↵
    Swanson LW (2003) The amygdala and its place in the cerebral hemisphere. Ann NY Acad Sci 985:174–184. doi:10.1111/j.1749-6632.2003.tb07081.x pmid:12724158
    OpenUrlCrossRefPubMed
  105. ↵
    Tiesinga P, Fellous JM, Sejnowski TJ (2008) Regulation of spike timing in visual cortical circuits. Nat Rev Neurosci 9:97–107. doi:10.1038/nrn2315 pmid:18200026
    OpenUrlCrossRefPubMed
  106. ↵
    Traub RD, Whittington MA, Colling SB, Buzsáki G, Jefferys JG (1996) Analysis of gamma rhythms in the rat hippocampus in vitro and in vivo . J Physiol 493:471–484. doi:10.1113/jphysiol.1996.sp021397 pmid:8782110
    OpenUrlCrossRefPubMed
  107. ↵
    Traub RD, Bibbig A, Fisahn A, LeBeau FE, Whittington MA, Buhl EH (2000) A model of gamma-frequency network oscillations induced in the rat CA3 region by carbachol in vitro . Eur J Neurosci 12:4093–4106. doi:10.1046/j.1460-9568.2000.00300.x pmid:11069606
    OpenUrlCrossRefPubMed
  108. ↵
    Turner BH, Herkenham M (1991) Thalamoamygdaloid projections in the rat: a test of the amygdala's role in sensory processing. J Comp Neurol 313:295–325. doi:10.1002/cne.903130208 pmid:1765584
    OpenUrlCrossRefPubMed
  109. ↵
    Tuunanen J, Pitkänen A (2000) Do seizures cause neuronal damage in rat amygdala kindling? Epilepsy Res 39:171–176. doi:10.1016/s0920-1211(99)00123-0 pmid:10759304
    OpenUrlCrossRefPubMed
  110. ↵
    Vertes RP (2004) Differential projections of the infralimbic and prelimbic cortex in the rat. Synapse 51:32–58. doi:10.1002/syn.10279 pmid:14579424
    OpenUrlCrossRefPubMed
  111. ↵
    Vinck M, van Wingerden M, Womelsdorf T, Fries P, Pennartz CM (2010) The pairwise phase consistency: a bias-free measure of rhythmic neuronal synchronization. Neuroimage 51:112–122. doi:10.1016/j.neuroimage.2010.01.073 pmid:20114076
    OpenUrlCrossRefPubMed
  112. ↵
    Wang XJ (2010) Neurophysiological and computational principles of cortical rhythms in cognition. Physiol Rev 90:1195–1268. doi:10.1152/physrev.00035.2008 pmid:20664082
    OpenUrlCrossRefPubMed
  113. ↵
    Wang XJ, Buzsáki G (1996) Gamma oscillation by synaptic inhibition in a hippocampal interneuronal network model. J Neurosci 16:6402–6413. doi:10.1523/JNEUROSCI.16-20-06402.1996
    OpenUrlAbstract/FREE Full Text
  114. ↵
    Washburn MS, Moises HC (1992) Electrophysiological and morphological properties of rat basolateral amygdaloid neurons in vitro . J Neurosci 12:4066–4079. doi:10.1523/jneurosci.12-10-04066.1992 pmid:1403101
    OpenUrlAbstract/FREE Full Text
  115. Weisskopf MG, Bauer EP, LeDoux JE (1999) L-type voltage-gated calcium channels mediate NMDA-independent associative long-term potentiation at thalamic input synapses to the amygdala. J Neurosci 19:10512–10519. doi:10.1523/JNEUROSCI.19-23-10512.1999 pmid:10575047
    OpenUrlAbstract/FREE Full Text
  116. ↵
    Whittington MA, Traub RD, Kopell N, Ermentrout B, Buhl EH (2000) Inhibition-based rhythms: experimental and mathematical observations on network dynamics. Int J Psychophysiol 38:315–336. doi:10.1016/s0167-8760(00)00173-2 pmid:11102670
    OpenUrlCrossRefPubMed
  117. ↵
    Witt A, Palmigiano A, Neef A, El Hady A, Wolf F, Battaglia D (2013) Controlling the oscillation phase through precisely timed closed-loop optogenetic stimulation: a computational study. Front Neural Circuits 7:49. doi:10.3389/fncir.2013.00049 pmid:23616748
    OpenUrlCrossRefPubMed
  118. ↵
    Woodruff AR, Sah P (2007) Networks of parvalbumin-positive interneurons in the basolateral amygdala. J Neurosci 27:553–563. doi:10.1523/JNEUROSCI.3686-06.2007 pmid:17234587
    OpenUrlAbstract/FREE Full Text
  119. ↵
    Zador A, Koch C, Brown TH (1990) Biophysical model of a Hebbian synapse. Proc Natl Acad Sci USA 87:6718–6722. doi:10.1073/pnas.87.17.6718 pmid:2168555
    OpenUrlAbstract/FREE Full Text
  120. ↵
    Zandvakili A, Kohn A (2015) Coordinated neuronal activity enhances corticocortical communication. Neuron 87:827–839. doi:10.1016/j.neuron.2015.07.026 pmid:26291164
    OpenUrlCrossRefPubMed

Synthesis

Reviewing Editor: Upinder S. Bhalla, National Centre for Biological Sciences, Tata Institute of Fundamental Research;

Decisions are customarily a result of the Reviewing Editor and the peer reviewers coming together and discussing their recommendations until a consensus is reached. When revisions are invited, a fact-based synthesis statement explaining their decision and outlining what is needed to prepare a revision will be listed below. The following reviewer(s) agreed to reveal their identity: Gaute Einevoll. Note: If this manuscript was transferred from JNeurosci and a decision was made to accept the manuscript without peer review, a brief statement to this effect will instead be what is listed below.

1. The reviewers agreed that the paper was interesting and technically

accomplished, with important findings.

2, The authors must present the anatomical organization of the modeled system,

including cellular geometries individually and with respect to the other cells

in the population.

3. The calculation for the LFP needs to be better justified in relation to

previous studies.

4. The novelty of the work should be better presented, and placed in context

of similar work especially pertaining to gamma in other brain areas.

5. The relationship between the reported gamma to other frequencies needs to

be better analyzed, both for bursts and entrainment.

6. The authors must provide their code as required by the

Journal policy: “Policy on Computer Code and Software

Computational models: New computational neuroscience computer code must be included with the initial submission as Extended Data to ensure reproducibility...”

They should really have done this at the first submission.

The detailed comments of the reviewers are appended as they provide more

details and bring up some useful additional points.

================================================================

Reviewer 1:

The paper reports a pioneering modeling study the a part of amygdala (basolateral complex). The model consists of 27000 neurons, modeled as 2-compartment or 3-compartment neurons, and interestingly and importantly the model predicts both spikes and local-field potentials (LFPs). I find the approach to be sound and the work very accomplished.

Comments to the Authors

The paper of Feng et al. reports a pioneering modeling study the a part of amygdala (basolateral complex). The model consists of 27000 neurons, modeled as 2-compartment or 3-compartment neurons, and interestingly and importantly the model predicts both spikes and local-field potentials (LFPs). I find the approach to be sound and the work very accomplished. Since I also find the questions addressed to be both important and timely, I am supportive of publishing of this work in eNeuro.

I have, however, some comments and questions that I would like the authors to address prior to publication:

1. I miss a figure of the anatomical organisation of the basolateral complex of

the amygdala. This organization, in particular on the geometric organisation of

the dendrites of the different neurons, could have a strong effect on the LFP.

2a) On l. 301 it is stated that the dendrites of the neurons are randomly oriented. Is this what is observed in the amygdala? (I believe that in other nuclei the dendrites could have a radial organisation.)

2b) If not, I would suggest that other geometrical configurations of the dendrites also should be considered. This may not only affect the overall amplitude of

the LFP, but also the shape of LFP power spectra. This follows from the observation that (i) correlated synaptic input combined with (ii) spatially aligned dendrites can give a boost of the LFP at low frequencies (see Leski et al, PLoS Comp Biol 2013)

2c) This could maybe contribute to an explanation for why the low-frequency powe of the LFP is smaller in the model than what is seen in experiments? (as mentioned on l. 564 and in the Discussion)

2d) Note that the lack of low-frequency power in predicted LFP was also observed

in models for cortical LFP, see Hagen et al, Cerebral Cortex (2016). It would be good if a connection to this previous study could be made.

3. On l. 424 it is assumed that the LFP is proportional to the density of neurons. Is there a justification for this? In, for example Linden et al, Neuron (2011), it was found that in some network regimes (low synaptic input correlation) the LFP scales like the square root of the density while with large synaptic input correlation the scaling was more linear.

Minor comments:

l. 71: anatomically detailed -&gt; anatomical realistic

l. 148: referral to “one dataset whose manuscript is in preparation”. This seems unorthodox.

l. 306-317: The used of “x” for product is unnecessary. I suggest to remove them throughout the manuscript.

l. 322: Must be stated that V is in units of mV (or whatever)

l. 332: The use of “*” for product s unnecessary. I suggest to remove them.

l. 389: Give formula for Ornstein-Uhlenbeck process

l. 454: once -&gt; one

l. 559: Can't anything be fit by a 1/f power law over a short frequency-range? What is meant?

l. 590: How can the entrainment of rates have anything to do with extracellular potentials in this model?

=============================================================

Reviewer 2

The authors constructed a biophysically and anatomically detailed model of the basolateral amygdala (BL), which reproduced several properties of the network observed experimentally. In particular, this model generated transient gamma oscillations, as reported by these authors previously. Using this model, they characterized spatial spread of gamma as well as network and synaptic properties.

I found their model comprehensive, the questions well motivated, and analysis thorough. I have some questions and comments though, as described below.

1. Given that there are several detailed ING-PING type models of gamma oscillations, the authors should specify in more detail which of their results are new in terms of modeling. I understand that parameters have been optimized for BL, which is a new thing, but some results, such as bursty nature of gamma, have been shown in other brain areas. For example, the bursty nature of gamma has been described by Bob Shapley's group, who mentioned gamma as “filtered noise”, see Burns et al., 2010 and 2011, JNS, Xing et al, 2012, JNS. Their earlier models were more abstract (based on variants of Wilson-Cowan), such as Kang et al., 2010, JCN, but more recently they also have detailed biophysical models (see work of Chariker et al., 2018, JNS). As mentioned by the authors themselves, several other groups have also used ING-PING type models to describe gamma. At the very least, the authors need to have a detailed Discussion about their model versus others, and in particular, if there is anything new that has not been shown previously in other parts of the brain.

2. Intermittent bursts of gamma (Figures 2 and 6): the model appears to have intermittent bursts of power at all frequencies, as seen in Figure 2E. Is there anything particular about the gamma band? If the authors really want to emphasize the properties of gamma bursts, they should compare gamma bursts with bursts at other frequencies. In this regard, recently it has been suggested that the bursty nature of gamma actually arises due to the variability in the spectral estimator (Subhash Chandran et al., 2018, JNP), and that previous studies have underestimated the burst duration. If the authors really want to emphasize the burst properties (Figure 6), they should ensure that their results do not suffer from signal processing artifacts (perhaps by changing the filter setting and testing if that changes the parameters) and also compare the burst statistics with nearby frequency bands (say between 40-60 and 80-100 Hz).

3. Gamma entrainment (Figure 3): I was surprised that the authors only focused on the gamma band instead of comparing the phase locking etc in the entire frequency range. For example, is the phase locking shown in Figures 3B and 3C unique to gamma band? What about the real data? The authors should show spike-field coherence or PPC for both real and model data here before dissecting the network and cellular mechanisms in later figures (in which some of these analyses are done).

4. Figure 5C shows an increase in coherence for the PN population at high frequencies. Any idea why that is the case?

5. Figure 7G - when the network does not produce gamma, the peaks/troughs etc have no meaning, right? So it is surprising that no difference in spike probability is observed?

Minor

1. Please label A1 and A2 in Figure 1.

2. Line 1307: “with” the effect, not “which”.

Author Response

Dear Dr. Bhalla and Reviewers,

We appreciate the positive review. The concerns raised have been addressed fully, and the manuscript has been updated appropriately, as we detail in the point-to-point responses below, with the review comments and comments in bold and plain font type, respectively.

The process has helped improve the manuscript, and we thank all of you for the comments.

Sincerely,

F Feng, D Headley, A Amir, V Kanta, Z Chen, D Pare, S Nair

Synthesis Statement for Author (Required):

1. The reviewers agreed that the paper was interesting and technically accomplished, with important findings.

We appreciate the comment.

2, The authors must present the anatomical organization of the modeled system, including cellular geometries individually and with respect to the other cells in the population.

We have revised Figure 1 to reflect the large scale spatial organization of our BL model and its cell densities (Figure 1A,B). We also provided an example of the random spatial orientation of the dendritic processes for our model neurons (Figure 1C). This organization fits the reports of Golgi stain studies in the rat BL (McDonald AJ 1984).

3. The calculation for the LFP needs to be better justified in relation to previous studies.

To address this, we empirically estimated the scaling factor for the LFP strength as a function of cell density. We systematically varied the cell density in the model and found how the LFP amplitude scaled as a function of it. This revealed that the LFP amplitude grows by the power of 0.67 as a function of cell density. All LFPs throughout the paper were replotted with the new scaling factor derived from this estimate. The methods section of the paper dealing with the LFP calculation has been rewritten to include the details of this calculation and make our choice of scaling factor clearer.

4. The novelty of the work should be better presented, and placed in context of similar work especially pertaining to gamma in other brain areas.

We added the following section to the discussion (starting line 922) to explain how our BL model of gamma is distinct from previous models of gamma based on cortical circuitry. We discuss how the microcircuitry of BL is not ideal for producing gamma oscillations, that our model is closer to a PING mechanism than the REI mechanism in visual cortex, and that unlike most other models we use a forward modeling approach to simulate how gamma appears in the LFP.

“BL microcircuitry presented unique challenges to gamma rhythmogenesis

The generation of gamma oscillations in our BL model was not a foregone conclusion. While qualitatively the BL network contains the microcircuitry required for producing gamma oscillations, quantitative differences in the properties of this network from previous PING models may have precluded gamma generation.

For a PING rhythm to emerge PNs must fire in advance of FSIs (Whittington et al. 2000), and thus increasing the average firing rates of PNs contributes to the strength of gamma oscillations produced in PING models and in slices (Adesnik and Scanziani 2010). PN firing rates in the BL are smaller than those in the cortical models used to simulate gamma by at least a factor of 6 (Börgers and Kopell 2008; Hoseini and Wessel 2016; Chariker et al. 2018), which may have diminished the network's ability to produce a PING type gamma rhythm. Also important is the strength of inhibitory feedback in the network (Börgers et al. 2012). Two factors that contribute to this are the strength of connections between PNs and FSIs and the number of FSIs in the network. We and others have found that the BL has a lower proportion of interneurons (McDonald and Mascagni 2001) compared with standard PING models (Börgers and Kopell 2008; Chariker et al. 2018, Hoseini and Wessel 2016). Moreover, the connection probability between PNs and FSIs in BL is lower (Woodruff and Sah 2007) than the aforementioned models. Thus, a priori the BL network may not be capable of producing gamma, and so our model was essential to at least test the possibility theoretically.

Another difference between our BL model and cortical-based PING models is a lower proportion of recurrent connectivity among PNs. Brunel and Wang (2003) found that including recurrent connections between excitatory units in a PING model could downshift the frequency of the network oscillation by half. On the other hand, our BL PING model has fairly little dependence on these connections, since their elimination did not substantially affect the prominence or peak frequency of gamma (Figure 4D). In addition, other models of gamma generation depend on the recurrent excitation among excitatory units. For instance, the REI model (Chariker et al. 2018) contains a network of reciprocally connected excitatory and inhibitory units. Unlike PING models, the gamma-periodic bursts of activity depend upon the buildup of recurrent excitation in the PN network, which raises the membrane potential for all units in the network, and drives them to fire synchronously. This mechanism for producing gamma is distinctly not PING, since the excitatory and inhibitory units activate simultaneously, instead of with the inhibitory units following the excitatory. Moreover, since BL PNs and FSIs in vivo fire during different phases of the gamma cycle (Amir et al. 2018), it is likely that the REI mechanism does not operate within the BL.

Thus, our model establishes the sufficiency of the BL network to produce gamma oscillations, and for those to arise via a PING mechanism. Despite that, it was also possible that any gamma produced by the model would be of insufficient power to comprise what is observed in the LFP in vivo. To address this, our model included a realistic estimation of the local field potential, which most PING models lack (Brunel and Hakim, 1999; Brunel and Wang, 2003; Borgers et al., 2005; Economo and White, 2012; Chalk et al., 2016; Hoseini and Wessel, 2016; Palmigiano et al., 2017; Traub et al., 2000; Bathellier et al., 2006; Neymotin et al., 2011; Chariker et al., 2018). Since our model LFP exhibited gamma oscillations that were close to the frequency and power of those recorded in vivo, it establishes the sufficiency of the local microcircuitry in BL to explain the gamma observed in LFP recordings in vivo.”

5. The relationship between the reported gamma to other frequencies needs to be better analyzed, both for bursts and entrainment.

While we agree that other frequency bands could be explored with this model, we intentionally excluded analyzing them for several reasons. First, from the outset it was clear that in vivo the power spectrum only showed a consistent prominence in the mid-gamma frequency band, with all other frequencies following the 1/f decline in power. This suggested to us that there was a distinct physiological phenomenon happening in the mid-gamma band apart from the other frequencies, and this warranted targeted investigation. Second, a priori we had little reason to investigate frequencies above the gamma band (&gt;120 Hz) in BL since there was virtually no information indicating their importance to BL function, and neurons in the BL tended not to entrain to them (Amir et al. 2018). Third, once we ran the model we found that the strength of its LFP for frequencies below gamma, &lt;30 Hz, was substantially weaker than that seen in vivo. Explanations for why this is the case and what it tells us about the origin of those low frequency signals in vivo are discussed in the manuscript (see Emergence of gamma oscillations in the BL model in the Results section and also Emulation of BL gamma in the Discussion section). However, this divergence gave us little confidence in the meaningfulness of any results arising from them. Thus, we felt it was best to focus the paper on the mid-gamma band.

That being said, the Reviewer 2 raised several important points regarding our focus on gamma and whether the mid-gamma band in particular is distinct in its relationship with network events. To address this, we have incorporated an additional analysis into the manuscript (Figure 2F) and provided analyses for the reviewer that establish the uniqueness of mid-gamma (see our responses to Reviewer 2 questions 2 and 3).

6. The authors must provide their code as required by the

Journal policy: “Policy on Computer Code and Software

Computational models: New computational neuroscience computer code must be included with the initial submission as Extended Data to ensure reproducibility...”

We have uploaded the code to the submission site as 'Extended data'

The detailed comments of the reviewers are appended as they provide more details and bring up some useful additional points.

================================================================

-

Reviewer 1:

The paper reports a pioneering modeling study the a part of amygdala (basolateral complex). The model consists of 27000 neurons, modeled as 2-compartment or 3-compartment neurons, and interestingly and importantly the model predicts both spikes and local-field potentials (LFPs). I find the approach to be sound and the work very accomplished.

Thank you.

Comments to the Authors

The paper of Feng et al. reports a pioneering modeling study the a part of amygdala (basolateral complex). The model consists of 27000 neurons, modeled as 2-compartment or 3-compartment neurons, and interestingly and importantly the model predicts both spikes and local-field potentials (LFPs). I find the approach to be sound and the work very accomplished. Since I also find the questions addressed to be both important and timely, I am supportive of publishing of this work in eNeuro.

I have, however, some comments and questions that I would like the authors to address prior to publication:

1. I miss a figure of the anatomical organisation of the basolateral complex of the amygdala. This organization, in particular on the geometric organisation of the dendrites of the different neurons, could have a strong effect on the LFP.

A new series of panels has been included (Figure 1A,B) illustrating the distribution of FSIs and PNs throughout the model, along with the random orientation of their processes and spatial extent (Figure 1C)

2a) On l. 301 it is stated that the dendrites of the neurons are randomly oriented. Is this what is observed in the amygdala? (I believe that in other nuclei the dendrites could have a radial organisation.)

Neurons in BLA do not have a layered structure as in the cortex, and the dendrites of PNs orient randomly (McDonald 1984). For this reason, we assumed random orientation for the dendrites. The relevant sentence has now been updated to include the reference (line 301):

“The dendrites of all neurons in the network had random orientations (McDonald, 1984).”

2b) If not, I would suggest that other geometrical configurations of the dendrites also should be considered. This may not only affect the overall amplitude of the LFP, but also the shape of LFP power spectra. This follows from the observation that (i) correlated synaptic input combined with (ii) spatially aligned dendrites can give a boost of the LFP at low frequencies (see Leski et al, PLoS Comp Biol 2013)

See answer to 2a) above.

2c) This could maybe contribute to an explanation for why the low-frequency powe of the LFP is smaller in the model than what is seen in experiments? (as mentioned on l. 564 and in the Discussion)

Please see responses to parts 2a) and to 2d) below.

2d) Note that the lack of low-frequency power in predicted LFP was also observed in models for cortical LFP, see Hagen et al, Cerebral Cortex (2016). It would be good if a connection to this previous study could be made.

This is again a good suggestion. We have modified the text, including both references, starting line 605, as follows:

“However, the model's spectrum had substantially less power below 30 Hz. This may reflect the absence of volume conducted activity or a lack of interactions with extrinsic structures. The former possibility was studied by Leski et al., 2013, who showed that lower-frequency LFP components have a larger spatial reach and extended further outside the active population than high-frequency components. It has also been found that increasing the correlations between extrinsic afferents boosts power in the low frequency band (Hagen et al. 2016; Leski et al. 2013). Our model likely underestimates the contribution of correlated extrinsic inputs since our naturalistic spike trains only represented a fraction of the structures projecting to the BL, each presumably with their own correlated spiking. Consistent with this, we found that replacing the naturalistic correlated spiking processes of our extrinsic afferents with rate matched independent Poisson processes further reduced the power of low frequencies (Figure 2B, green line). Importantly, gamma power from 50 to 80 Hz was unaffected by this change, supporting the hypothesis that it is locally generated and not imposed by upstream structures or volume conducted from adjacent regions.”

Leski S, Linden H, Tetzlaff T, Pettersen KH, Einevoll GT (2013) Frequency dependence of signal power and spatial reach of the local field potential. PLoS Comput Biol 9:e1003137.

Hagen E, Dahmen D, Stavrinou ML, Linden H, Tetzlaff T, van Albada SJ, Grun S, Diesmann M, Einevoll GT (2016) Hybrid Scheme for Modeling Local Field Potentials from Point-Neuron Networks. Cerebral cortex 26:4461-4496.

3. On l. 424 it is assumed that the LFP is proportional to the density of neurons. Is there a justification for this? In, for example Linden et al, Neuron (2011), it was found that in some network regimes (low synaptic input correlation) the LFP scales like the square root of the density while with large synaptic input correlation the scaling was more linear.

We agree with the reviewer and have addressed this issue in the revision. Specifically, we have revised the text, as follows (starting line 444):

“To correct for this, we rescaled the LFP by a correction factor that was estimated following Linden et al. (2011) who suggest that LFP scales as the square root of the number of neurons, N, for uncorrelated synaptic inputs. However, that scaling factor was derived from a network of pyramidal neurons positioned in a disk with uniformly oriented dendrites, and was not for density per se but actually the radius of the disk with density held constant. To determine the scaling factor in our model, we systematically varied the number of cells used to calculate the LFP from a full model run. Varying the cell count from 1000 up to 27000, we found that the standard deviation of the LFP scaled with density following N0.67. Since cell density in rats is reported to range from 2.5×104 to 2×105/mm3 (Tuunanen and Pitkanen, 2000; Salm et al., 2004; Pego et al., 2008; Rubinow and Juraska, 2009; Chareyron et al., 2011), while the model density is 9,840 neurons/mm3 (27,000 neurons, 1.4mm x 1.4mm x 1.4mm), the LFP correction factor would corresponding range from 1.9 to 7.5 [(〖2.5×〖10〗^4/9.84×〖10〗^3)〗^0.67 to (〖2×〖10〗^5/9.84×〖10〗^3)〗^0.67]. We chose the average of this range, 4.7, to scale the model LFP.”

Linden H, Tetzlaff T, Potjans TC, Pettersen KH, Grun S, Diesmann M, Einevoll GT (2011) Modeling the spatial reach of the LFP. Neuron 72:859-872.

Details for the reviewers on how to calculate the scaling factor: Linden et al. (2011) state that the scaling factor for predicting LFP amplitude is partly dependent on the correlation of synaptic inputs to cells in the network. We first calculated the correlation between the synaptic inputs using eqn. 15 in Linden et al. (2011) and found it to be 0.0012, which shows that it is close to being uncorrelated. We then used the model to estimate the dependence of LFP on the number of cells in our BL network (N). This was performed by first calculating the contribution to the LFP from all cells in the 27,000-cell network, and then iteratively dropping a random 1000 cells from the computation of the LFP, till the lower limit of N=1,000.

Such a process yielded the red circles (normalized to the value for the N=1,000 case) in the plot on the left. Blue and black dashed lines are for the linear, N or √N proportions, where N is number of cells in the network. These are again normalized to the N=1,000 case. This experiment allowed us to estimate the relationship between cell density in the model and changes in LFP amplitude.

A power law fit through the red curve shows the scaling factor to be 0.67, which is close to the scaling factor of 0.5 as reported in Linden et al. 2011. This is consistent with the low average correlation coefficient from our network that we calculated above.

Minor comments:

l. 71: anatomically detailed -&gt; anatomical realistic

Done.

l. 148: referral to “one dataset whose manuscript is in preparation”. This seems unorthodox.

We have reworded it as “one unpublished dataset”.

l. 306-317: The used of “x” for product is unnecessary. I suggest to remove them throughout the manuscript.

We have removed “x” but used “*” at some locations where the potential for ambiguity might exist (e.g., eqns. 4-6).

l. 322: Must be stated that V is in units of mV (or whatever)

Done. We state that right after equation 1.

l. 332: The use of “*” for product s unnecessary. I suggest to remove them.

Please see response two items above.

l. 389: Give formula for Ornstein-Uhlenbeck process

We have added a new section with equations and table to explain the way we modeled the background inputs (starting line 387):

“Specifically, stochastic background input Ibck had two independent components, excitatory and inhibitory, for both PNs and FSIs, modeled as follows:

I_bck=g_e (t)(V-E_e )+g_i (t)(V-E_i )

where g_e (t) and g_i (t) are time-dependent excitatory and inhibitory conductances, respectively; E_e=0 mV and E_i=-75 mV are respective reverse potentials. The two are modeled as Ornstein-Uhlenbeck processes described below:

(dg_e (t))/dt=-1/τ_e [g_e (t)-g_e0 ]+√(D_e ) χ_1 (t)

(dg_i (t))/dt=-1/τ_i [g_i (t)-g_i0 ]+√(D_i ) χ_2 (t)

where g_e0 and g_i0 are average conductances, τ_e and τ_i are time constants, D_e and D_i are noise 'diffusion' coefficients, χ_1 (t) and χ_2 (t) are Gaussian white noise of zero mean and unit standard deviation (see Table 6 for these parameter values).

The above two stochastic differential equations can be numerically modeled by using the following update rule:

g_e (t+∆t)=g_e0+[g_e (t)-g_e0 ]exp⁡(-∆t/τ_e)+A_e N_1 (0,1)

g_i (t+∆t)=g_i0+[g_i (t)-g_i0 ]exp⁡(-∆t/τ_i)+A_i N_2 (0,1)

Where N_1 (0,1) and N_2 (0,1) are normal random numbers. A_e and A_i are amplitude coefficients with

A_e=√(∂_e^2 [1-exp⁡(-2∆t/τ_e)] )

A_i=√(∂_i^2 [1-exp⁡(-2∆t/τ_i)] )

Table 6 lists the parameter values used in the point-conductance model.”

Table 6. Parameters related to point-conductance model

l. 454: once -&gt; one

Done

l. 559: Can't anything be fit by a 1/f power law over a short frequency-range? What is meant?

A power law is scale invariant, meaning that it will have the same shape and coefficients irrespective of the range of values used for its calculation.

l. 590: How can the entrainment of rates have anything to do with extracellular potentials in this model?

An issue that frequently arises in the analysis of unit entrainment to oscillations in the LFP is the degree to which a neuron's own transmembrane currents contribute to the extracellular field and cause spurious coupling. The most common example of this is the contamination of extracellular spike waveforms into the gamma band, leading to an artifactually high entrainment between spiking and LFP gamma recorded on the same electrode. To contend with this, experimenters often record spiking on one electrode and the LFP on another. One could imagine a similar phenomenon occurring in our model, where the distance dependent decay in entrainment to gamma merely reflected a decrease in the contribution of a particular cell to the LFP. In our model because we have access to the exact contribution of each neuron to the simulated LFP, we measure how fast this decays and see if it matches or decays more slowly than the change in entrainment. This analysis revealed that a neuron's contribution to the LFP decayed more than twice as fast as its entrainment to gamma, indicating that the spurious kind of spike/field coupling mentioned above was unlikely.

=============================================================

-

Reviewer 2

The authors constructed a biophysically and anatomically detailed model of the basolateral amygdala (BL), which reproduced several properties of the network observed experimentally. In particular, this model generated transient gamma oscillations, as reported by these authors previously. Using this model, they characterized spatial spread of gamma as well as network and synaptic properties.

I found their model comprehensive, the questions well motivated, and analysis thorough. I have some questions and comments though, as described below.

We appreciate the comment.

1. Given that there are several detailed ING-PING type models of gamma oscillations, the authors should specify in more detail which of their results are new in terms of modeling. I understand that parameters have been optimized for BL, which is a new thing, but some results, such as bursty nature of gamma, have been shown in other brain areas. For example, the bursty nature of gamma has been described by Bob Shapley's group, who mentioned gamma as “filtered noise”, see Burns et al., 2010 and 2011, JNS, Xing et al, 2012, JNS. Their earlier models were more abstract (based on variants of Wilson-Cowan), such as Kang et al., 2010, JCN, but more recently they also have detailed biophysical models (see work of Chariker et al., 2018, JNS). As mentioned by the authors themselves, several other groups have also used ING-PING type models to describe gamma. At the very least, the authors need to have a detailed Discussion about their model versus others, and in particular, if there is anything new that has not been shown previously in other parts of the brain.

Please see our response to Synthesis statement #4 above where we list the additional text we have added to address this concern.

2. Intermittent bursts of gamma (Figures 2 and 6): the model appears to have intermittent bursts of power at all frequencies, as seen in Figure 2E. Is there anything particular about the gamma band? If the authors really want to emphasize the properties of gamma bursts, they should compare gamma bursts with bursts at other frequencies. In this regard, recently it has been suggested that the bursty nature of gamma actually arises due to the variability in the spectral estimator (Subhash Chandran et al., 2018, JNP), and that previous studies have underestimated the burst duration. If the authors really want to emphasize the burst properties (Figure 6), they should ensure that their results do not suffer from signal processing artifacts (perhaps by changing the filter setting and testing if that changes the parameters) and also compare the burst statistics with nearby frequency bands (say between 40-60 and 80-100 Hz).

We were not aware of the recent study by Subhash Chandran et al. (2018) when we submitted our manuscript. However, we checked the matching pursuit algorithm in Chandran et al and felt it may not be applicable to our study, since it detects gamma bursts in 1D waveforms while ours focus on 3D space in Figure 6. The large number of recording sites for our 3D LFP monitoring and the extended recording durations would also have made the matching pursuit algorithm computationally prohibitive. Moreover, since Subhash Chandran et al., 2018 didn't make any claims about whether the longer burst durations they detected had any functional significance, we opted to keep our results aligned with those from most of the prior studies (e.g. Xing et al. 2012, Amir et al., 2018).

As per the reviewer's suggestion, we did re-analyze burst duration in adjacent gamma bands to see if they differed from those found at mid-gamma. The results of that analysis can be found in the following plot.

Burst duration in the mid-gamma band, which was the focus of this study, is significantly longer than the other bands we processed, perhaps because bursts were more synchronized the in mid-gamma band. Note that the burst duration we calculated in our model is based on the watershed algorithm that identifies contiguous convex regions of elevated gamma power across three spatial and one temporal dimension (see Materials and Methods). The longer burst durations in the mid-gamma band stem from the higher synchronization of gamma bursts among each individual electrode detected.

However, the Chandran et al. paper did motivate us to tackle the issue of gamma burst durations earlier in the paper. Fortuitously, we used wavelets to perform spectral analyses of our simulated and in vivo LFPs in Figure 2. In the Chrandran et al. paper, it was reported that gamma burst durations derived from wavelets were on par with those uncovered by the matching pursuit algorithm for high-amplitude gamma bursts. Thus, we calculated the probability distributions of the mid-gamma bursts in our wavelet power spectrograms with respect to their durations and amplitudes (Figure 2F). This revealed that the distribution for our model gamma bursts matched those recorded in vivo. Additionally, the mean duration for higher amplitude bursts (greater than 80th percentile in amplitude) was ~30 ms shorter than those from Figure 6. This could have arisen because gamma bursts can move throughout the BL, and so a single recording site is less likely to be both the origin and end point of a gamma burst, and thus will only 'see' the burst for a portion of its total existence.

Chandran Ks S, Seelamantula CS, Ray S (2018) Duration analysis using matching pursuit algorithm reveals longer bouts of gamma rhythm. Journal of neurophysiology 119:808-821.

Xing D, Shen Y, Burns S, Yeh CI, Shapley R, Li W (2012) Stochastic generation of gamma-band activity in primary visual cortex of awake and anesthetized monkeys. The Journal of neuroscience : the official journal of the Society for Neuroscience 32:13873-13880a.

3. Gamma entrainment (Figure 3): I was surprised that the authors only focused on the gamma band instead of comparing the phase locking etc in the entire frequency range. For example, is the phase locking shown in Figures 3B and 3C unique to gamma band? What about the real data? The authors should show spike-field coherence or PPC for both real and model data here before dissecting the network and cellular mechanisms in later figures (in which some of these analyses are done).

Please see our response to the Synthesis statement #5. Briefly, we focused on mid-gamma given its prominence in vivo and deficiencies in the interpretations of frequencies above and below the gamma band. While we strongly feel that including analyses outside the mid-gamma band would overwhelm an already packed paper and detract from its focus, we are more than willing to share with the reviewer the requested analysis.

These are the circular r values for the data from the Amir et al. 2018 study, which recorded from the BL in vivo, and several runs of our model with different random seeds. As you can see, for the most part the entrainment values predicted by the model in each frequency band fall within the distribution of values observed in vivo. Even though the results are concordant, we are still hesitant to include them because (1) we have reservations about the physiological interpretation of the low-frequency activities observed in our model and (2) frequencies above the gamma band (&gt;100 Hz) have yet to be ascribed to particular functions in the BL.

4. Figure 5C shows an increase in coherence for the PN population at high frequencies. Any idea why that is the case?

The increase in coherence above 100 Hz reflects the fact that IPSCs tended to follow EPSCs. In PNs because the EPSCs are so sparse, there is a strong lagged coordination between the onset of the EPSC and IPSC. Such coordination is much weaker in for FSIs. This can be seen in the figure below:

This is something we should have again recognized and addressed, and do so now. Since this effect is largely trivial, and the signal power was low at frequencies above 100 Hz, and so have limited the relevant plots in Fig. 5 to 100 Hz.

5. Figure 7G - when the network does not produce gamma, the peaks/troughs etc have no meaning, right? So it is surprising that no difference in spike probability is observed?

Even when the LFP lacks gamma oscillations it can still have peaks and troughs in the gamma band, since once can always filter the LFP in the frequency range of mid-gamma. Since the LFP reflects transmembrane currents, and in general neuronal depolarization will produce an LFP trough, while hyperpolarization will produce a peak, it could have been that the coordination we observed in figure 7G was merely reflecting this phenomenon irrespective of gamma. Consequently, Figure 7G was used to show that when gamma oscillations were present in BL, their troughs were the most sensitive part for generating spikes. By contrast, when BL lacked gamma oscillations, this significantly diminished the LFP phase dependent sensitivity of spikes to EPSCs (also compare Y axis on two panels), even though firing rates were unchanged. We believe that this analysis demonstrates that gamma oscillation are critical in the BL for the pacing the sensitivity of PNs to extrinsic inputs, and consequently synchronizing their spiking so that downstream targets are more effectively driven.

Minor

1. Please label A1 and A2 in Figure 1.

Done. Although now that the model geometry has been included in Figure 1, the corresponding panels are now D1 and D2.

2. Line 1307: “with” the effect, not “which”.

Done.

Back to top

In this issue

eneuro: 6 (1)
eNeuro
Vol. 6, Issue 1
January/February 2019
  • Table of Contents
  • Index by author
  • Ed Board (PDF)
Email

Thank you for sharing this eNeuro article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Gamma Oscillations in the Basolateral Amygdala: Biophysical Mechanisms and Computational Consequences
(Your Name) has forwarded a page to you from eNeuro
(Your Name) thought you would be interested in this article in eNeuro.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Print
View Full Page PDF
Citation Tools
Gamma Oscillations in the Basolateral Amygdala: Biophysical Mechanisms and Computational Consequences
Feng Feng, Drew B. Headley, Alon Amir, Vasiliki Kanta, Ziao Chen, Denis Paré, Satish S. Nair
eNeuro 11 January 2019, 6 (1) ENEURO.0388-18.2018; DOI: 10.1523/ENEURO.0388-18.2018

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Respond to this article
Share
Gamma Oscillations in the Basolateral Amygdala: Biophysical Mechanisms and Computational Consequences
Feng Feng, Drew B. Headley, Alon Amir, Vasiliki Kanta, Ziao Chen, Denis Paré, Satish S. Nair
eNeuro 11 January 2019, 6 (1) ENEURO.0388-18.2018; DOI: 10.1523/ENEURO.0388-18.2018
del.icio.us logo Digg logo Reddit logo Twitter logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Jump to section

  • Article
    • Abstract
    • Significance Statement
    • Introduction
    • Materials and Methods
    • Results
    • Discussion
    • Conclusions
    • Footnotes
    • References
    • Synthesis
    • Author Response
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF

Keywords

  • amygdala
  • gamma oscillations
  • computational model
  • biophysical model
  • extracellular potential

Responses to this article

Respond to this article

Jump to comment:

No eLetters have been published for this article.

Related Articles

Cited By...

More in this TOC Section

New Research

  • Allopregnanolone effects on inhibition in hippocampal parvalbumin interneurons
  • LINCs are vulnerable to epileptic insult and fail to provide seizure control via on-demand activation
  • Pregabalin silences oxaliplatin-activated sensory neurons to relieve cold allodynia
Show more New Research

Neuronal Excitability

  • Allopregnanolone effects on inhibition in hippocampal parvalbumin interneurons
  • Similar presynaptic action potential-calcium influx coupling in two types of large mossy fiber terminals innervating CA3 pyramidal cells and hilar mossy cells
  • Automated Image Analysis Reveals Different Localization of Synaptic Gephyrin C4 Splice Variants
Show more Neuronal Excitability

Subjects

  • Neuronal Excitability

  • Home
  • Alerts
  • Visit Society for Neuroscience on Facebook
  • Follow Society for Neuroscience on Twitter
  • Follow Society for Neuroscience on LinkedIn
  • Visit Society for Neuroscience on Youtube
  • Follow our RSS feeds

Content

  • Early Release
  • Current Issue
  • Latest Articles
  • Issue Archive
  • Blog
  • Browse by Topic

Information

  • For Authors
  • For the Media

About

  • About the Journal
  • Editorial Board
  • Privacy Policy
  • Contact
  • Feedback
(eNeuro logo)
(SfN logo)

Copyright © 2023 by the Society for Neuroscience.
eNeuro eISSN: 2373-2822

The ideas and opinions expressed in eNeuro do not necessarily reflect those of SfN or the eNeuro Editorial Board. Publication of an advertisement or other product mention in eNeuro should not be construed as an endorsement of the manufacturer’s claims. SfN does not assume any responsibility for any injury and/or damage to persons or property arising from or related to any use of any material contained in eNeuro.