Distinct gamma oscillations in the distal dendritic fields of the dentate gyrus and the CA1 area of mouse hippocampus

The molecular layer of the dentate gyrus and the anatomically adjacent stratum lacunosum-moleculare of CA1 area, represent afferent areas at distinct levels of the hippocampal trisynaptic loop. Afferents to the dentate gyrus and CA1 area originate from different cell populations, including projection cells in entorhinal cortex layers two and three, respectively. To determine the organization of oscillatory activities along these terminal fields, we recorded local field potentials from multiple sites in the dentate gyrus and CA1 area of the awake mice, and localized gamma frequency (30–150 Hz) oscillations in different layers by means of current source density analysis. During theta oscillations, we observed different temporal and spectral organization of gamma oscillations in the dendritic layers of the dentate gyrus and CA1 area, with a sharp transition across the hippocampal fissure. In CA1 stratum lacunosum-moleculare, transient mid-frequency gamma oscillations (CA1-gammaM; 80 Hz) occurred on theta cycle peaks, while in the dentate gyrus, fast (DG-gammaF; 110 Hz), and slow (DG-gammaS; 40 Hz) gamma oscillations preferentially occurred on troughs of theta waves. Units in dentate gyrus, in contrast to units in CA1 pyramidal layer, phase-coupled to DG-gammaF, which was largely independent from CA1 fast gamma oscillations (CA1-gammaF) of similar frequency and timing. Spike timing of units recorded in either CA1 area or dentate gyrus were modulated by CA1-gammaM. Our experiments disclosed a set of gamma oscillations that differentially regulate neuronal activity in the dentate gyrus and CA1 area, and may allow flexible segregation and integration of information across different levels of hippocampal circuitry.


Introduction
The hippocampus receives inputs from associative cortical areas, and is a key structure for spatial navigation and episodic memory (O'Keefe and Nadel 1978;Buzsaki and Moser 2013). The majority of its external afferents originate in the entorhinal cortex (EC), and reach the hippocampal formation via the temporoammonic (TA), and perforant pathways (PP), that terminate in anatomically adjacent areas, separated by the hippocampal fissure (Witter 2012). The PP is involved in contextual memory , and connects layer two (L2) of the EC to the dentate gyrus (DG; van Groen et al. 2003;Witter 2012). This pathway forms the very initial stage of information processing in hippocampal circuitry. By contrast, the TA, which is indispensable in temporal association memory (Suh et al. 2011), originates from the more upstream layer three (L3) of the EC. The TA innervates a more downstream level of the hippocampal trisynaptic loop via terminals located in the stratum lacunosum-moleculare of the CA1 area (van Groen et al. 2003;Witter 2012).
Based on their distinct theta phase preference, and spectral and spatial distributions, gamma frequency  oscillations in the CA1 area have been further classified as slow (30-60 Hz; CA1-gamma S ), mid-frequency (50-100 Hz; CA1-gamma M ), and fast (90-150 Hz; CA1-gamma F ) oscillations (Belluscio et al. 2012;Scheffer-Teixeira et al. 2012;Schomburg et al. 2014;Lasztoczi and Klausberger 2016). In the CA1 area, different gamma oscillations are associated with pathways selectively terminating in different layers (Colgin et al. 2009;Scheffer-Teixeira et al. 2012;Schomburg et al. 2014;Lasztoczi and Klausberger 2014). For example, transient CA1-gamma M oscillations appear strongest in the stratum lacunosummoleculare, and preferentially occur on the peaks of theta oscillations (Scheffer-Teixeira et al. 2012;Schomburg et al. 2014;Lasztoczi and Klausberger 2016). Simultaneously, a similar oscillation in L3 of the medial EC (mEC) occurs, and pyramidal cells in L3 discharge (Schomburg et al. 2014;Mizuseki et al. 2009). Thus, CA1-gamma M oscillations may result from concerted synaptic activity in the TA pathway (Schomburg et al. 2014;Yamamoto et al. 2014;Colgin et al. 2009;Lasztoczi and Klausberger 2014) and regulate the communication from EC L3 to area CA1 of the hippocampus (Schomburg et al. 2014;Lasztoczi and Klausberger 2016). Afferents projecting to the DG arise from a different set of parent cells, and DG gamma oscillations appear different from CA1-gamma M (Scheffer-Teixeira et al. 2012). Axons of the PP originate from calbindin immunonegative, EC L2 stellate cells (Ray et al. 2014;Kitamura et al. 2014). Stellate cells, together with GABAergic L2 basket cells, form a microcircuit capable of autonomously generating powerful theta-nested gamma oscillations independent of the L3, in vitro (Pastoll et al. 2013;Middleton et al. 2008;Couey et al. 2013). Interestingly, most mEC L2 projection cells, and hilar mossy cells that give rise to a second major excitatory pathway to the DG (Witter 2012;Scharfman, 2016), fire counter-phase to L3 cells, on the trough of theta waves in vivo (Mizuseki et al. 2009;Senzai and Buzsaki 2017). To test if CA1 and DG network operations are segregated by different temporal organizations of gamma oscillations, we recorded local field potentials (LFP) and calculated instantaneous current source density (CSD) in the dendritic layers of DG and CA1. In addition, we recorded spiking activity of CA1 and DG units and investigated their phase coupling to different gamma oscillations in head-restrained mice during movement.

Materials and methods
All animal procedures were carried out under licences approved by the Austrian Ministry of Science and in accordance with the relevant regulations of the Medical University of Vienna. Adult male C57BL/6J mice were implanted with a plastic head plate under isoflurane anaesthesia (3-4% for induction and 1.5-2% for maintenance). After recovery (1-2 days), the animals were habituated to head restraint (additional 1-2 days), were water restricted (1 ml water/day), and trained to perform unidirectional runs in a 4 m long linear virtual reality maze (Phenosys), for a small water reward. Animals controlled the maze by rotating an air-supported styrofoam ball in all directions, of which only rotations along the long axis of the maze were registered. For craniotomy (9-129 days after head plate implantation), animals were briefly reanaesthetized, a small cranial window was drilled above their right dorsal hippocampus (1.3 mm lateral and 1.9 mm caudal from the Bregma), the dura was removed, and the brain surface was sealed with silicone (Kwik-Sil, World Precision Instruments). At least 4 h were given as recovery before the first recordings were performed.
On recording days, mice were head-fixed in the apparatus, and the recording electrodes were inserted. To record LFP, we used a linear silicon probe with 16 recording sites at 50 lm spacing (Neuronexus), inserted 1.3 mm lateral and 1.7-2.0 mm caudal from the Bregma, with the dorsoventral positioning guided by the profile of theta oscillations, sharp waves, and ripple oscillations. Spiking activity was recorded from the CA1 stratum pyramidale or stratum granulosum of the DG, using another silicon probe inserted in an 8°-10°angle, 300-600 lm away (Neuronexus; four shanks spaced at 150 lm with four contacts/shank in tetrode arrangement, or two shanks spaced at 200 lm distance with eight staggered contacts/shank; in final position the shank tips were 100-300 lm from the other silicon probe). After the recording locations were reached, the brain surface was covered (saline, mineral oil or wax), and additional 15-20 min was allowed for the electrode positions to stabilize. The exact position of the recording electrodes within the hippocampal formation was inferred post hoc, by comparing electrode tracks in histological analysis to the electrophysiological activity profiles Klausberger 2014, 2016). Only experiments where the linear silicon probe spanned all the layers from CA1 stratum pyramidale to DG molecular layer, and at least three contacts were ventral from the hippocampal fissure, were included. Units were analysed only from shanks unequivocally positioned in CA1 stratum pyramidale or DG stratum granulosum. From nine experiments (in five animals), data on CA1 place cells and their coupling to CA1 gamma oscillations during place field traversals have been reported in an earlier publication (Lasztoczi and Klausberger 2016).
Signals from the silicon probes were pre-amplified (19; RA16AC, Tucker-Davis), amplified (10009), band-pass filtered (0.1-475 Hz for LFP, and 0.1-6000 Hz for unit recordings; Lynx-8 signal conditioners; Neuralynx), and digitized (2 kHz for LFP, and 20 kHz for unit recordings; Power1401mkII controlled by Spike2 software, Cambridge Electronic Design). All analyses were limited to periods with theta oscillations, which almost exclusively corresponded to periods of movement, despite no behavioural variables were considered during their definition. A single contact was selected (typically from CA1 stratum oriens or pyramidale), and theta periods were semi-automatically detected when the theta (5-12 Hz) to delta (2-4 Hz) power ratio exceeded 4 in the LFP signal, with period boundaries checked and adjusted manually if necessary (Lasztoczi and Klausberger 2016;Lapray et al. 2012). In all analyses, instantaneous theta phase was taken from stratum pyramidale or equivalent (stratum oriens), by linear interpolation between peaks (180°) and troughs (0°and 360°) detected in theta-filtered (5-12 Hz), and down-sampled (400 Hz) LFP trace.
Instantaneous CSD traces on contact n were derived from LFP traces of the linear silicon probe, by estimating the second spatial derivative at every sampled time point t, using the equation: where LFP n,t , LFP n-1,t , and LFP n?1,t are the LFP samples at time t, from contacts n and the two neighbouring contacts, respectively, and Dz is the spacing (50 lm; Lasztoczi and Klausberger 2016;Mitzdorf 1985;Bragin et al. 1995).
In CSD traces sinks/sources are presented downwards/upwards. To analyse gamma oscillations, CSD traces were subjected to wavelet transformation (Lasztoczi and Klausberger 2016), with a complex Morlet wavelet (20-150 Hz; 53 logarithmically equidistant frequencies; wavelet parameters of 1 and 1.5; maximum sink at 0°and 360°and maximum source at 180°). Using complex wavelets allowed us to extract the instantaneous phase and amplitude values for each wavelet scale, at any time point.
To analyse cross-frequency phase-amplitude coupling, amplitudes of CSD gamma oscillations were Z-scored over theta periods of the entire recording, averaged within 18°t heta phase bins, and averaged across theta cycles. Modulation index spectra were generated by subtracting the minima from the maxima of these averages at each frequency (wavelet scale). For each gamma oscillation a frequency range was defined in a selected contact by identifying a local maximum of the modulation index, and extending the range until it decreased to either 80%, or a local minimum. Selected contacts were typically from stratum lacunosum-moleculare (CA1-gamma M ), stratum radiatum (CA1-gamma S ), stratum pyramidale (CA1gamma F ), and the DG molecular layer (DG-gamma F and DG-gamma S ). To quantify the theta modulation of gamma oscillations, average Z-scored amplitudes from the selected contact were averaged over the frequency ranges. Theta ranges for each gamma oscillation were defined as theta phases when these averages were positive. To calculate theta phase dependent phase coherence between CSD signals in contact pairs, in each tenth of a theta cycle, a time sample was drawn, and the corresponding instantaneous phase spectra were extracted in all contacts, and for the contact pairs phase difference spectra were generated. By repeating this procedure for all theta cycles in a recording, for each contact pair, each tenth of theta cycle and each frequency, we derived a phase-difference distribution and the corresponding statistics, such as the mean phase difference, the phase locking value (PLV; the mean vector length of angles; Lachaux et al. 1999), and an estimated Rayleigh's P value. To quantify gamma phase coherence between two spatially distinct gamma oscillations, we averaged the PLV values from the phase coherence matrix of the contact pair relevant for those two oscillations, over the overlap in frequency and theta phase ranges. Spikes were extracted and clustered from silicon probe recordings using standard procedures (Csicsvari et al. 2003a;Lasztoczi and Klausberger 2016;Royer et al. 2012;Rossant et al. 2016). Extracted spikes were sorted into clusters putatively originating from single units by first running an automatic clustering algorithm, KlustaKwik (Rossant et al. 2016;Harris et al. 2000), followed by manual and automatic refinement of the clusters in Klusters or KlustaViewa software (Hazan et al. 2006;Rossant et al. 2016). In experiments with silicon probe shanks with eight staggered contacts, we used the masked EM algorithm of the new version of KlustaKwik for both spike extraction and clustering (Rossant et al. 2016). Units were classified as putative principal cells (or putative GABAergic cells), if they had an overall firing rate below (above) 3 Hz, a spike width at 90% of the peak amplitude above (below) 0.5 ms, and an event autocorrelogram value below 10 ms (above 15 ms). Event autocorrelogram value was calculated by taking the count-weighted average of offset times at which event autocorrelogram counts exceeded the mean count calculated over a 50 ms window. Other units remained unclassified. Modulation of spike timing by theta oscillations was tested in units with at least 20 spikes during theta oscillations. To analyse the modulation of spike timing by gamma oscillations, instantaneous wavelet transform spectra corresponding to unit spikes were extracted, and summarized with phase coupling statistics calculated in a spectral manner. Phase coupling for a particular gamma oscillation in a particular recording contact was considered significant, if we observed non-uniform (P \ 0.05, Rayleigh test) phase distribution within the frequency range of gamma oscillation (individually defined for each experiment, see above) for at least half of the frequencies (scales). Phase coupling to gamma oscillations was tested only if at least 200 spikes occurred during theta oscillations. Phase coupling to a particular gamma oscillation in general was considered significant if it was significant in Brain Struct Funct (2017) 222:3355-3365 3357 the contact used to define the frequency range of that oscillation (see above), and in the case of DG-gamma F , DG-gamma S , CA1-gamma M , and CA1-gamma S , in one additional designated contact, typically in the DG molecular layer and the CA1 stratum radiatum. To quantify coupling strength and mean phase to a particular gamma oscillation, r values for the most relevant contact were averaged over the respective frequency range and the mean phase was calculated over the scales with significant coupling within this range. Unless stated otherwise, data are presented as mean ± SD, or mean angle ± circular SD, as appropriate.

Results
We recorded LFP, and calculated instantaneous CSD from multiple sites, in different layers of the CA1 area and the suprapyramidal blade of the DG, in the right dorsal hippocampus of head-fixed awake mice (N = 17 recording experiments from 9 animals). During theta oscillations, we observed strong and diverse oscillatory activity in the gamma frequency range (30-150 Hz), in CSD traces from all layers of the CA1 and DG (Fig. 1). Gamma oscillations were transient, and their occurrence depended on the phase of ongoing theta oscillation, and appeared strikingly different in the stratum lacunosum-moleculare of the CA1 area, and the molecular layer of the DG. In the CA1 stratum lacunosum-moleculare, high amplitude gamma oscillation transients were often observed on the peak of theta cycles, measured in the LFP of the CA1 pyramidal cell layer. These oscillations corresponded to CA1-gamma M (Schomburg et al. 2014;Lasztoczi and Klausberger 2016), reached maximum amplitude at 193°± 8°, and occupied a theta phase range between 117°and 279°, and a frequency range between 57 and 86 Hz (Figs. 1, 2a). By contrast, CSD traces in the DG displayed short bouts of faster oscillations, and longer transients of slower oscillations at the trough of theta cycles, counter-phase to CA1-gamma M (Fig. 1). These oscillations have been mentioned before (Scheffer-Teixeira et al. 2012), here we characterize them in detail and we term them DG fast, and DG slow gamma oscillations (DG-gamma F and DG-gamma S ), respectively. DG-gamma F and DG-gamma S occupied a narrow range of the theta cycle around its trough (-63°to 63°and -63°to 81°, respectively; on average maximum amplitudes were at 8°± 13°and 6°± 21°, respectively; Figs. 1 and 2a).
Despite the similar theta phase preference of DG-gamma F and DG-gamma S , the modulation index spectra in DG showed biphasic frequency distributions, with the two oscillations occurring at frequency ranges 75-150, and 31-46 Hz, respectively (Figs. 1b, 2a). In the frequency domain, DG-gamma S did not overlap with CA1-gamma M but DG-gamma F did (Fig. 2a). Within the theta cycle, however, DG-gamma F and DG-gamma S did not overlap with CA1-gamma M (Figs. 1b, 2a, b).
In line with previous data (Schomburg et al. 2014;Lasztoczi and Klausberger 2016), in CA1 stratum radiatum and stratum pyramidale, we observed CA1-gamma S (32-39 Hz) and CA1-gamma F (92-150 Hz) with frequency ranges similar to DG-gamma S and DG-gamma F (Figs. 1,  2a). Both CA1-gamma S and CA1-gamma F were widely distributed (from -135°to 45°, and from -175°to 45°, respectively) on the descending phase of the theta waves recorded in the CA1 pyramidal cell layer (maximal amplitudes at 316°± 17°, and 301°± 11°, respectively; Figs. 1b, 2a), and thus overlapped substantially with the DG oscillations (Fig. 2b). This prompted us to analyse, if oscillatory phase of DG-gamma S -CA1-gamma S , and DGgamma F -CA1-gamma F oscillation pairs correlated in time. For quantification, we used the phase locking value (PLV; Lachaux et al. 1999; see Materials and Methods). Although often reaching significance thresholds (at a = 0.05, Rayleigh test), even when multiple comparisons were accounted for, the PLV values for DG-gamma F -CA1gamma F pairs were small (0.061 ± 0.028, range 0.018-0.129), indicating that these two gamma oscillations were largely independent (Fig. 2c). By contrast, pairs of DG-gamma s -CA1-gamma s showed stronger phase coherence (0.297 ± 0.066, range 0.134-0.386), suggesting that the two oscillations may partly reflect the same underlying oscillatory process (Fig. 2c). Unlike the DG oscillations, CA1-gamma F and CA1-gamma S showed substantial overlap with CA1-gamma M , during the early descending phase of the theta cycle (Figs. 1, 2b).
Along with gamma oscillations in the hippocampal formation, we recorded spikes of units from the granule cell layer of DG (N = 122 units, 44 putative principal cells, 21 putative interneurons, in 7 recording experiments from 3 mice), and the pyramidal layer of CA1 (N = 294 units, 161 putative pyramidal cells, 23 putative interneurons, in 10 recording experiments from 6 mice). Most DG (N = 112, 96%) and CA1 (N = 278, 95%) units were significantly (P \ 0.05 with Rayleigh test) modulated by theta oscillations, and coupled with variable strength to the late descending (DG units) or descending (CA1 units) phase of the CA1 theta cycle (Fig. 3a, c).
Spectra of spike timing modulation by oscillation phase in the gamma frequency range showed a high degree of variability in different contacts, and across different units, even from the same hippocampal subfield. Importantly, DG-gamma F entrained 37% (N = 37 of 99) of DG units with mean r of 0.072 ± 0.05 (range 0.018-0.261), but had only weak influence on the spike-timing (r 0.030 ± 0.012, range 0.014-0.042; P = 0.012 compared to DG units with Mann-Whitney U test) in a small subpopulation of CA1 units (N = 5 of 284, 1.8%; P = 1.6 9 10 -22 , v 2 = 95, compared to DG units with v 2 test; Fig. 3b, d-f). However, CA1-gamma M phase-modulated the spike timing of CA1 and DG units to similar extent (r in CA1: 0.053 ± 0.025, range 0.017-0.148; r in DG: 0.048 ± 0.025, range 0.016-0.139; P = 0.3, Mann-Whitney U test), and in similar proportion of units (CA1: N = 77 of 284, 27%, DG: N = 35 of 99; 35%; P = 0.12; v 2 = 2.4; v 2 test; Fig. 3b, d-f). Despite CA1-gamma F and DG-gamma F showing little coherence, and CA1 units firing independent of DG-gamma F , DG units were modulated by CA1gamma F in similar proportion (N = 26 of 99, 26%; P = 0.21; v 2 = 1.6; v 2 test), yet substantially weaker (r 0.032 ± 0.015, range 0.010-0.073; P = 1.8 9 10 -6 , Mann-Whitney U test) than CA1 units (N = 94 of 284, 33%; r 0.078 ± 0.063, range 0.019-0.334; Fig. 3e, f). Phase of slow gamma oscillations, recorded in either DG or CA1, modulated the firing of more DG units than CA1 units (DG-gamma S , N = 65 of 99 DG units, 65%, and 54 of 284 CA1 units, 19%; P = 5.9 9 10 -18 ; v 2 = 75; v 2 Fig. 1 Gamma oscillations in the dentate gyrus and CA1 area of the hippocampus. a High-pass filtered (at 25 Hz) CSD traces calculated for silicon probe contacts located in different hippocampal input layers (source is upwards) and the corresponding LFP trace recorded from the CA1 pyramidal layer (bottom trace), during theta oscillations. Layers are indicated on the left. Theta troughs from the CA1 pyramidal layer LFP are marked by vertical dotted lines, for reference. Coloured arrowheads indicate instances of distinct gamma oscillations, as indicated. Dentate gyrus, and CA1 area are indicated by red and green background colours, respectively. b Left, mean amplitude of gamma oscillations (Z-score of CSD wavelet amplitude), plotted for each contact as a function of theta phase in CA1 pyramidal layer (18°bins, the theta cycle is duplicated for visualization), and gamma frequency (53 logarithmically equidistant wavelet scales between 20 and 150 Hz). Right, mean amplitude spectra (grey), and phase-amplitude modulation index spectra (black), of CSD from individual contacts. Coloured arrowheads indicate peak modulation index positions for various gamma oscillations in the most relevant contacts. Note the markedly different oscillatory dynamics across the fissure. s. gran. granule cell layer, s.mol. molecular layer, s. l-m. stratum lacunosum-moleculare, s.rad. stratum radiatum, s. pyr. stratum pyramidale, CSD current source density, LFP local field potential, DG dentate gyrus, CA1 cornu ammonis area 1 Fig. 2 Co-occurrence and coherence of gamma oscillations in the dentate gyrus and CA1 area. a Frequency distribution (upper plot) and theta phase-amplitude modulation of different gamma oscillations. Theta phase here, and throughout the manuscript, was measured in the pyramidal cell layer of CA1 area. b Overlap in the occurrence (defined as a positive mean Z score) of different gamma oscillation pairs during the theta cycle. Note that the DG gamma F -CA1 gamma F and DG gamma S -CA1 gamma S oscillation pairs show substantial overlap in their frequency and theta phase distributions. c Phase coherence (measured as phase locking value; PLV) between oscillatory activities in contact pairs, displayed as a function of theta phase, and frequency. In the top row, and the left column, mean CSD oscillation amplitude Z-scores are plotted for selected contacts in stratum moleculare, stratum lacunosum-moleculare, stratum radiatum, and stratum pyramidale, as indicated. At the intercept of these, phase coherence between pairs of contacts (displayed as colour-coded phase locking value) is plotted as a function of theta phase (18°bins, the theta cycle is duplicated for visualization), and gamma frequency (53 logarithmically equidistant wavelet scales between 20 and 150 Hz). In plots at the intercept of a contact with itself, the theta phase and frequency ranges of different gamma oscillations prominent in the particular contact are displayed, for reference (also plotted as black dotted lines on the left column amplitude plots). White dotted lines indicate these same ranges on coherence plots, to indicate overlaps in frequency and theta phase ranges of different gamma oscillations test; CA1-gamma S , N = 73 of 99 DG units, 74%, and 155 of 284 CA1 units, 55%; P = 0.00028; v 2 = 13; v 2 test). Also, DG units displayed stronger modulation than the CA1 units (DG-gamma S , r 0.098 ± 0.056 vs. 0.053 ± 0.025, ranges 0.033-0.327 and 0.014-0.132; P = 1.3 9 10 -8 , Mann-Whitney U test; CA1-gamma S , r 0.095 ± 0.038 vs. 0.074 ± 0.031, ranges 0.031-0.196 and 0.017-0.164; P = 6*10 -5 , Mann-Whitney U test; Fig. 3e, f). These data indicate that DG-gamma F represents an oscillation independent of CA1-gamma F , while DGgamma S and CA1-gamma S may, at least partly, reflect the same underlying oscillatory process.

Discussion
Here we report and characterise a class of gamma oscillations expressed in the molecular layer of the DG. Using CSD analysis in different layers of the hippocampus, in awake head-fixed mice, we observed short transients of around 110 Hz oscillations localised to the recording sites positioned within the DG, and restricted to a narrow temporal window around the trough of theta oscillations. We termed these DG fast gamma oscillations (DG-gamma F ). In the rodent hippocampus, multiple gamma oscillations have been described (Bragin et al. 1995;Csicsvari et al. 2003b;Colgin et al. 2009;Scheffer-Teixeira et al. 2012;Schomburg et al. 2014;Klausberger 2014, 2016). However, its high frequency, and preferred theta phase, distinguishes DG-gamma F from CA1-gamma S (frequency around 35 Hz) and CA1-gamma M (occurring on theta peaks). Moreover, DG-gamma F and CA1-gamma F display little phase coherence and unlike CA1-gamma F , DGgamma F influence the spike-timing of DG units, but not CA1 units. These two lines of evidence suggest that DGgamma F are different also from CA1-gamma F , and thus represent a distinct hippocampal gamma oscillation, likely corresponding to the wide-band, fast gamma oscillations observed in LFP recordings from the rat DG (Scheffer-Teixeira et al. 2012). Extracellular potentials resulting from neuronal spikes have significant power in the gamma F frequency range (Ray and Maunsell 2011;Belluscio et al. 2012). This raises the possibility that uncoordinated multiunit activity (axonal or somatic) may have contributed to both DG-gamma F and CA1-gamma F (Belluscio et al. 2012;Scheffer-Teixeira et al. 2013), and caused spurious phase modulation of spike-timing (Zanos et al. 2011;Scheffer-Teixeira et al. 2013). However, the amplitude of unitary extracellular spikes decreases by more than an order of magnitude over 100 lm Schomburg et al. 2012), and we detected phase coupling of units to DG-gamma F and CA1-gamma F in electrodes typically [100 lm apart. Moreover, we observed coherent DG-gamma F oscillations at multiple non-neighbouring, dendritic recording sites. Thus, gamma F , and in particular DGgamma F , result from spiking and synaptic activities temporally coordinated over spatially extended populations of CA1 and DG cells (Schomburg et al. 2012). Spiking activity in CA1 stratum pyramidale and-more surprisingly-in DG as well, was phase-modulated by CA1gamma F . Since DG-gamma F and CA1-gamma F are not coherent, this weak modulation in the DG may occur on descending phase of the theta cycle, when CA1-gamma F , but not DG-gamma F are present, and may be mediated by a yet unknown common oscillatory input, or GABAergic back-projections from CA1 to the DG (Katona et al. 2016;Klausberger et al. 2005;Fuentealba et al. 2010).
The temporal organization and frequency distribution of gamma oscillations is markedly different in the molecular layer of the DG and the stratum lacunosum-moleculare of the CA1 area, suggesting that network operations are segregated in these two structures. Glutamatergic innervation to the CA1 stratum lacunosum-moleculare arise from the pyramidal cells of mEC L3 (van Groen et al. 2003;Yamamoto et al. 2014;Witter 2012;Suh et al. 2011), where the frequency and theta phase preference of gamma oscillations mirror CA1-gamma M , suggesting that coordinated firing of L3 pyramidal cells on theta cycle peaks (Mizuseki et al. 2009) are transmitted via the TA, and generate CA1-gamma M (Schomburg et al. 2014;Colgin et al. 2009;Yamamoto et al. 2014;Sun et al. 2015;Klausberger 2014, 2016). Units in the lateral EC (lEC) display weaker theta modulation (Deshmukh et al. 2010) and their contribution to CA1-gamma M is poorly understood. On the other hand, the origins of DG gamma oscillations remain more elusive. Most glutamatergic terminals in the molecular layer of DG, where DG-gamma F and DG-gamma S were recorded, arise from the reelin-expressing, calbindin immunonegative principal cells of the mEC (stellate cells; Ray et al. 2014;Kitamura et al. 2014) and the lEC (Leitner et al. 2016), and the recurrent collaterals of hilar mossy cells (Scharfman 2016). Most L2 projection cells in the mEC (Mizuseki et al. 2009;Quilichini et al. 2010), and hilar mossy cells (Senzai and Buzsaki 2017), discharge in a narrow time window around the trough of theta oscillations, coincident with gamma oscillations in the DG (both fast and slow) and counterphase to CA1 gamma M . Indeed, when rhythmically excited at theta frequency in vitro, the L2 microcircuit of the mEC (including stellate cells and mutually connected GABAergic basket cells) generates transient synchronous gamma oscillations at every theta cycle, even when its L3 afferents are severed (Pastoll et al. 2013;Middleton et al. 2008;Couey et al. 2013), pointing to the PP as a potential source of DG-gamma F . Alternative explanations also exist (e.g. GABAergic networks within the DG generating gamma oscillations; Bragin et al. 1995;Pernia-Andrade and Jonas 2014), and further experimental work is needed to discriminate between these possibilities. Despite the lack of DG cell innervation by L3 (van Groen et al. 2003;Witter 2012), the spike timing of DG units was detectably correlated to the phase of CA1-gamma M . Potential mechanism for such entrainment include the GABAergic back-projection from CA1 to the DG (Katona et al. 2016;Klausberger et al. 2005;Fuentealba et al. 2010) or the entrainment of L2 network by L3 cells at theta peaks in the absence of its intrinsic rhythm as observed in vitro (Middleton et al. 2008).
The phase of DG-gamma S and CA1-gamma S were correlated during theta oscillations, and both slow gamma oscillations modulated the spike timing of units recorded from either the DG or the CA1 area. This indicates that DG-gamma S is not independent from CA1-gamma S , and slow gamma oscillations may rather represent manifestations of the same underlying oscillatory process, that regulates information processing across the entire hippocampal formation, and perhaps beyond. Slow gamma oscillations have been linked to the communication between CA3 and CA1, and have been suggested to arise in the CA3 microcircuit (Csicsvari et al. 2003b;Colgin et al. 2009). In this case GABAergic (Lasztoczi et al. 2011) and glutamatergic (Scharfman 2007) back-projections from CA3 to DG may mediate the entrainment of DG units by CA1-gamma S . Alternatively, slow gamma oscillations may originate upstream from CA3 in the DG or hilus, as originally proposed (Bragin et al. 1995). This latter scenario also offers a simple explanation to our finding that DG units couple stronger than CA1 units to both DG-gamma S and CA1-gamma S . Indeed, it has been recently demonstrated that slow gamma oscillations persist when CA3 communication to CA1 area is blocked (Middleton and McHugh 2016) and Granger causality analysis indicated that slow gamma oscillations are imposed on CA3 pyramidal cells, by upstream microcircuits in DG (Hsiao et al. 2016). The preferred firing phase of both DG granule cells (Mizuseki et al. 2009;Senzai and Buzsaki 2017) and hilar mossy cells (Senzai and Buzsaki 2017) is consistent with their roles in generating DG-gamma S and CA1-gamma S .
Multiple gamma oscillations are instrumental in regulating the communication along different extrinsic and intrinsic connections of the hippocampus (Colgin et al. 2009;Schomburg et al. 2014); yet experimental support on the actual roles of distinct gamma oscillations is sparse. Synchronization of CA1 pyramidal cells to CA1-gamma M contributes to hippocampal network operations (Lasztoczi and Klausberger 2016;Yamamoto et al. 2014;Schomburg et al. 2014), and different gamma oscillations have been linked to different navigation strategies in rodents (Cabral et al. 2014). Our data indicates that the complexity and flexibility of animal cognition and behaviour may be supported by a diversity of gamma oscillations exceeding that previously thought, and that synaptic communication to different levels of hippocampal information processing is regulated by a set of distinct gamma oscillations. b Fig. 3 Phase modulation of spike-timing of units in the dentate gyrus and area CA1 by theta and different gamma oscillations. a Left, theta phase histogram of firing of one unclassified unit, recorded from the dentate gyrus. Right, mean theta phase of firing (abscissa), and modulation strength (mean r value, ordinate) of all significantly modulated putative principal cells (green), putative interneurons (red), and unclassified units (black), recorded from the dentate gyrus. b Spike-timing modulation by gamma oscillations during theta oscillations of one unit, recorded from the dentate gyrus (same unit as in a). Left, amplitude modulation of gamma oscillations by theta oscillations, in CSD traces from selected contacts in stratum moleculare of the dentate gyrus, and stratum lacunosum-moleculare, stratum radiatum, and stratum pyramidale of the CA1 area, as indicated. Middle, spike density (greyscale-coded) of the unit, plotted as a function of CSD gamma phase and frequency for the same contacts. The mean phase values for the frequencies significantly phase-modulating the spike timing are plotted in black. Right, modulation strength spectra (mean r values) for the unit, relative to gamma oscillations in different contacts (values at frequencies significantly modulating the spike-timing are plotted in black). Frequency ranges of the different gamma oscillations are displayed by colour-coded bars at the plots of most relevant contacts. c Same as in A for units recorded from the CA1 area. d Same as in b, for one unit recorded in the CA1 area (the same putative interneuron as in c). e Spike-timing modulation strength (mean r values) of units recorded in the dentate gyrus (DG) and the CA1 area, by different gamma oscillations. Only units with significantly non-uniform phase distribution are shown. Putative principal cells, interneurons, and unclassified units are plotted in green, red, and black, respectively. Asterisks denote significant differences by Mann-Whitney U test with Holm-Bonferroni correction procedure for multiple comparisons at a general a = 0.05. Note that CA1 units do not phase couple to DG-gamma F and DG units are only weakly phase modulated by CA1-gamma F . f Percentage of DG and CA1 units with spike-timing significantly modulated by different gamma oscillations. Asterisks denote significant differences by the v 2 test with Holm-Bonferroni correction procedure for multiple comparisons at a general a = 0.05