Abstract
Gamma frequency (30–100 Hz) oscillations in the mature cortex underlie higher cognitive functions. Fast signaling in GABAergic interneuron networks plays a key role in the generation of these oscillations. During development of the rodent brain, gamma activity appears at the end of the first postnatal week, but frequency and synchrony reach adult levels only by the fourth week. However, the mechanisms underlying the maturation of gamma activity are unclear. Here we demonstrate that hippocampal basket cells (BCs), the proposed cellular substrate of gamma oscillations, undergo marked changes in their morphological, intrinsic, and synaptic properties between postnatal day 6 (P6) and P25. During maturation, action potential duration, propagation time, duration of the release period, and decay time constant of IPSCs decreases by ∼30–60%. Thus, postnatal development converts BCs from slow into fast signaling devices. Computational analysis reveals that BC networks with young intrinsic and synaptic properties as well as reduced connectivity generate oscillations with moderate coherence in the lower gamma frequency range. In contrast, BC networks with mature properties and increased connectivity generate highly coherent activity in the upper gamma frequency band. Thus, late postnatal maturation of BCs enhances coherence in neuronal networks and will thereby contribute to the development of cognitive brain functions.
Introduction
Fast rhythmic network activity patterns in the gamma frequency band (30–100 Hz) are thought to be critically important for higher brain functions. Gamma oscillations have been proposed to represent reference signals for temporal encoding of information, sensory binding of features into a coherent percept, and storage and recall of information (Gray and Singer, 1989; Ribary et al., 1991; Buzsáki, 2006). In the hippocampus, amplitude and frequency of gamma oscillations change markedly during postnatal maturation (Lahtinen et al., 2002; Leinekugel et al., 2002). In rodents, gamma oscillations are first detected between postnatal day 6 (P6) and P10, and both coherence and frequency increase to reach adult levels at the fourth postnatal week (Lahtinen et al., 2002; Leinekugel et al., 2002). It is believed that these changes underlie the maturation of higher cognitive functions both in rodents and humans (Ben-Ari et al., 2007; Yurgelun-Todd, 2007). However, the mechanisms involved in the maturation process of gamma activity in neuronal networks are insufficiently understood.
Networks of GABAergic interneurons play a key role in the generation of gamma oscillations. Among the various types of GABAergic cells (Buhl et al., 1994; Freund and Buzsáki, 1996; Miles et al., 1996; Parra et al., 1998; McBain and Fisahn, 2001; Somogyi and Klausberger, 2005), fast-spiking parvalbumin (PV)-expressing basket cells (BCs) are thought to be particularly important for network function (for review, see Bartos et al., 2007). Previous studies suggested that the fast signaling properties of BCs are necessary to optimize synchronization and to provide frequency tuning of network oscillations to the gamma band (Bartos et al., 2001, 2002; Galarreta and Hestrin, 2001; Buzsáki and Draguhn, 2004; Vida et al., 2006; Fuchs et al., 2007). The key role of BCs in the generation of gamma activity raises the intriguing possibility that developmental changes in functional properties and connectivity of BC–BC synapses are directly related to alterations in power and frequency of gamma oscillations. Consistent with this idea, during development the abundance of perisomatic inhibitory synapses increases (Tyzio et al., 1999), and the time course of miniature IPSCs in hippocampal principal cells becomes faster (Hollrigel and Soltesz, 1997; Cohen et al., 2000). Both factors will have a strong influence on the coherence and frequency of gamma activity (Traub et al., 1996; Wang and Buzsáki, 1996; Bartos et al., 2001, 2002). However, a systematic analysis of functional properties of BC–BC synapses in different phases of postnatal development has not been performed.
In the present study, we investigated the maturation of intrinsic and synaptic characteristics of BCs in the dentate gyrus of glutamate decarboxylase 67–green fluorescent protein (GAD67–GFP)-expressing mice between the first and the fourth postnatal week by applying a combined neuroanatomical, electrophysiological, and computational strategy. We find that intrinsic and synaptic differentiation converts BCs from slow to fast signaling devices at multiple levels of cellular function. Increased axonal conduction velocity, faster kinetics of unitary IPSCs at BC output synapses, and higher connectivity in the network will support the developmental increase in coherence and frequency of gamma oscillations observed in vivo (Lahtinen et al., 2002).
Materials and Methods
Patch-clamp recording from BCs.
To identify interneurons at early developmental stages, heterozygous GAD67–GFP (Δneo) mice expressing GFP under the control of the endogenous GAD67 gene promoter were used (Tamamaki et al., 2003). Transgenic mice were crossed with wild-type mice (C57BL/6), and GFP expression in the offspring was monitored by illuminating the heads with UV light (<400 nm). Transverse hippocampal slices (300 μm thickness) were cut from brains of 6- to 25-day-old mice using a Vibratome (DTK-1000; Dosaka; or a home-made slicer). Animals were killed by decapitation, in accordance with national and institutional guidelines. Patch pipettes were pulled from borosilicate glass tubing (outer diameter, 2 mm; inner diameter, 1 mm). When filled with intracellular solution, the resistance was 2.2–3.5 MΩ. GFP-expressing interneurons were identified during the experiment using epifluorescence illumination. Recordings were obtained from neurons in the dentate gyrus under visual control using infrared differential interference contrast video microscopy (Stuart et al., 1993; Koh et al., 1995). BCs chosen for recording had somata located at the granule cell (GC) layer–hilus border. GCs examined had cell bodies located in the outer third of the GC layer. The recording temperature was 33–34°C.
Two Axopatch 200B amplifiers (Molecular Devices) were used for paired recordings. Presynaptic cells were held in the I-Clamp normal mode. For the identification of passive and active membrane properties, BC recordings were performed in the I-Clamp fast mode. The holding potential was −70 mV and maintained by hyperpolarizing current injections (<100 pA). The series resistance (Rs) in the presynaptic neuron (5–18 MΩ) was fully compensated. For paired recordings, the presynaptic neuron was stimulated at a frequency of 0.2 Hz. Action potentials were initiated by brief current pulses (duration, 2 ms; amplitude, 0.2–1.4 nA). The postsynaptic cell was held in the V-Clamp mode at −70 mV with Rs compensation (∼65–90%; lag, 10–35 μs; Rs before compensation, 5–12 MΩ). The stationarity of Rs in the postsynaptic neuron was assessed from the amplitude of the capacitive current in response to a 10 mV pulse, and the compensation was readjusted during the experiment when necessary. Action potentials and IPSCs were filtered at 5–10 kHz using the built-in four-pole, low-pass Bessel filters of the amplifiers and were digitized at 20–40 kHz using a 1401plus laboratory interface (Cambridge Electronic Design) connected to a Pentium personal computer. Igor programs (F-Pulse, U. Fröbe Physiologisches Institut, Freiburg, Germany) were used for stimulus generation and data acquisition.
Antidromic action potentials were evoked in BCs by extracellular stimulation in the GC layer using glass pipettes filled with Na+-rich HEPES-buffered solution (2–3 MΩ; stimulus duration, 0.2 ms; stimulus amplitude, <30 V; repetition frequency, 0.2 Hz). In these experiments, 10 μm 6-cyano-7-nitroquinoxaline-2,3-dione (CNQX), 20 μm bicuculline methiodide (BIC), and 100 μm d-2-amino-5-phosphonopentanoic acid were added to the bath solution.
Solutions.
The physiological extracellular solution contained (in mm) 125 NaCl, 25 NaHCO3, 25 glucose, 2.5 KCl, 1.25 NaH2PO4, 2 CaCl2, and 1 MgCl2 (equilibrated with a 95% O2/5% CO2 gas mixture). In a subset of experiments, slices were stored in a solution containing (in mm) 87 NaCl, 25 NaHCO3, 2.5 KCl, 1.25 NaH2PO4, 0.5 CaCl2, 7 MgCl2, 25 glucose, and 75 sucrose. Intrinsic membrane properties from BCs were recorded using an intracellular solution that contained (in mm) 135 K-gluconate, 20 KCl, 0.1 EGTA, 2 MgCl2, 2 Na2ATP, and 10 HEPES. Paired recordings and recordings during bath application of Zolpidem were performed with an intracellular solution containing (in mm) 110 K-gluconate, 40 KCl, 0.1 EGTA, 2 MgCl2, 2 Na2ATP, and 10 HEPES; the pH was adjusted to 7.2 with KOH, and the osmolarity was 310–315 mOsm. Paired recordings were performed in the presence of 10 μm CNQX in the extracellular bath solution to block EPSCs. In some experiments (five pairs), 20 μm BIC was additionally added to block unitary IPSCs. BIC and tetrodotoxin (TTX) were purchased from Sigma, CNQX and Zolpidem were obtained from Tocris, and other chemicals were from Merck, Sigma, Riedel-de Haën, or Gerbu.
Data analysis.
Resting membrane potential was determined immediately after transition into the whole-cell configuration. Total cell capacitance (cm) and input resistance (Rin) were measured under voltage-clamp conditions from the current during a 1 s, 10 mV voltage pulse. To determine the membrane time constant (τm), short hyperpolarizing pulses (0.2 ms, 0.5–1 nA) were applied in the current-clamp mode in the presence of 10 μm CNQX, 20 μm BIC, and 0.5 μm TTX; 300–500 voltage traces were averaged; 100% to 10% of the decay was logarithmically transformed, and the late component was analyzed by linear regression. Single action potential parameters (peak amplitude, half-duration, maximal rate of rise and decay) were measured from the action potential threshold determined as the first point in the voltage trajectory in which the slope exceeded 20 V s−1 (Bekkers and Delaney, 2001). Membrane potentials reported in the text were not corrected for junction potentials.
Functional properties of unitary IPSCs were determined from averages of 10–50 traces including failures. IPSCs were aligned to the steepest point in the rise of the presynaptic action potential. The synaptic latency was determined as the time interval between the steepest point in the rise of the action potential and the onset of the IPSC; the onset point was determined from the intersection of a line through the 20 and 80% points with the baseline. The rise time of evoked IPSCs was determined as the time interval between the points corresponding to 20 and 80% of the peak amplitude. The peak current was determined as the maximum within a window of 1–4 ms duration after the presynaptic action potential. The decay phase of the IPSCs was fitted with the sum of two exponentials [A exp(−t/τ1) + B exp(−t/τ2)], using a nonlinear least-squares fit algorithm; time constants are reported as an amplitude-weighted mean [τw = (Aτ1 + Bτ2)/(A+ B)]. A trace was classified as failure when the peak amplitude was less than three times the SD of the baseline. Coefficients of variation (CVs; SD/mean) of synaptic latencies of unitary IPSCs were calculated from 30 to 50 traces during stationary periods.
To determine the velocity of axonal action potential propagation, antidromic action potentials were evoked in BC axons at four to six different stimulation sites in the GC layer. The distance between stimulation and recording site was measured on the video screen. The conduction latency was obtained as the time interval between the beginning of the stimulus artifact and the steepest point in rise of the somatically recorded action potential. Average conduction latencies were measured for 20–30 consecutive action potentials evoked at a given stimulation site. To determine the conduction velocity, average latency was plotted against distance, and data were analyzed by linear regression.
Home-made programs (Stimfit; C. Schmidt-Hieber and P. Jonas, Physiologisches Institut I, Freiburg, Germany) and Mathematica 4.1.2 (Wolfram Research) were used for data analysis. All values are given as mean ± SEM. Significance of differences was assessed by a nonparametric two-tailed Mann–Whitney test, and correlations between two parameters were examined by Spearman rank correlation analysis; the corresponding significance levels are indicated as p values.
Morphological reconstruction and analysis.
Morphological analysis of BCs was performed as reported previously (Bartos et al., 2001, 2002). Neurons were filled with biocytin (0.1–0.2%) during the recording. After withdrawal of the pipettes, slices were fixed in 2.5% paraformaldehyde, 1.25% glutaraldehyde, and 15% saturated picric acid in 100 mm phosphate buffer (PB), pH 7.4. Filled cells were visualized using avidin-biotinylated peroxidase complex and 3,3′-diaminobenzidine tetrahydrochloride (DAB) as chromogen.
Axons and dendrites of all recorded neurons were examined by light microscopy using a 63× oil-immersion objective [numerical aperture (NA) 1.4]. From all morphologically identified BCs, a subset of four young BCs (P6–P10) and five mature BCs (P18–P21) was selected for detailed morphological reconstruction using Neurolucida 7 (MicroBrightField). Morphometric data analysis was performed using NeuroExplorer (MicroBrightField).
Immunohistochemistry.
For immunohistochemical double-labeling studies, biocytin-filled cells were fixed in 4% paraformaldehyde. After washing in PBS (0.025 m, pH 7.3), slices were incubated in a monoclonal antibody against PV (mouse; 1:10,000; Swant) in PBS containing 5% goat serum and 0.3% Triton X-100 for 24 h at 22°C. The secondary antibody (goat anti-mouse conjugated with Alexa Fluor 568; 1:5000; Invitrogen) was applied together with streptavidin conjugated with Alexa Fluor 647 (1:500; Invitrogen) in PBS and 0.3% Triton X-100 for 6–12 h at 22°C. Slices were washed in 0.1 m PB and embedded in Prolong Antifade (Invitrogen). Labeled neurons were examined with a confocal microscope (LSM 510; Zeiss).
For a quantitative analysis on the colocalization of GAD67–GFP and PV, animals were deeply anesthetized with Narkodorm-n (180 mg/kg, i.p.; Alvetra), and the hearts were surgically exposed for perfusion fixation. First, the vascular system was flushed by circulating 0.9% saline for 1 min. This was followed by transcardial perfusion with a freshly prepared fixative solution containing 4% paraformaldehyde in 0.1 m PB for 13 min. After perfusion, brains were removed, and tissue blocks containing the hippocampus were cut in 50-μm-thick slices for PV immunohistochemistry.
To examine changes in the number and density of GAD67–GFP-labeled cells during development, confocal image stacks were taken from slice preparations of the dentate gyrus from two young (P9, 13 slices, 50 μm thickness) and two mature (P20, 16 slices) mice (40× oil-immersion objective, NA 1.3; Plan Neofluar; Zeiss). GFP-labeled cells were counted in a randomly selected 100 × 100 μm region of interest, which was centered on the GC layer (P9, total 102 cells; P20, total 98 cells). Shrinkage of slices after mounting was independent of the age (∼50%), indicating that similar tissue volumes were compared.
Network simulations.
Simulations were performed using a previously developed inhibitory interneuron network model (Vida et al., 2006) in the software environment of Neuron 6.1 (Hines and Carnevale, 1997). Fast-spiking interneurons were represented as single compartments endowed with voltage-gated Na+ and K+ conductances (Wang and Buzsáki, 1996). To reproduce the longer τm and larger Rin in young interneurons, the specific leak conductance was assumed as 0.13 to 0.05 mS cm−2, and the membrane surface area was taken as 100 and 60 μm2 for mature and young BCs, respectively. To mimic the action potential waveform of young BCs, we reduced the density of voltage-gated Na+ conductances from 0.035 to 0.025 mS cm−2. This change resulted in action potentials with 80% peak amplitude and 175% duration compared with those in the mature modeled interneurons.
The interneuron network models were assembled from 200 neurons arranged on a virtual ring with 50 μm distance between adjacent cells (Bartos et al., 2002). Each neuron was randomly connected to a subset of its 100 nearest neighbors by inhibitory synapses with a mean probability of 0.4 in young networks (this study) and of 0.6 in mature networks (Sik et al., 1995) to account for the experimentally observed differences in connection probability. Additionally, each neuron was randomly connected to four of its eight nearest neighbors by electrical synapses (transcellular conductance, 0.01 mS cm−2) (Bartos et al., 2002).
Properties of synapses in the model were based on experimental data. In the mature network, the 20–80% rise time of the conductance was set to 0.3 ms, and the weighted decay time constant was 1.9 ms (Bartos et al., 2002). To model the slower conductance in the young network, the rise time was set to 0.5 ms, and the decay time constant was 8 ms (present data). The synaptic latency of the postsynaptic conductance consisted of a constant synaptic delay (0.5 ms) and a variable conduction delay, calculated from the distance between presynaptic and postsynaptic cells along the circumference of the ring. Conduction velocity was set to 0.25 m s−1 in the mature network and 0.13 m s−1 in the young network, and the cell-to-cell distance was reduced from 50 to 40 μm to account for the reduced size of the young dentate gyrus.
Neurons were initialized to their resting membrane potential, and activity was initiated by applying tonic excitatory currents (I) to each neuron with amplitudes randomly chosen from a normal distribution (with mean Iμ and SD Iσ). The heterogeneity of the excitatory drive (Iσ /Iμ) was 10%. To examine the emergence of synchronization from a disordered initial condition, driving currents were applied to individual neurons with randomized onset times taken from a uniform distribution in the range −150 ms ≤ t < −100 ms. Chemical and electrical synapses were inactive at −150 ms ≤ t < 0 and enabled at t = 0. In all simulations, the time step was 10 μs. The amplitude of the tonic excitatory current and the synaptic conductance were calculated from density values (Iμ = 0.5–2.5 μA cm−2 and gsyn = 0.001–1 mS cm−2, respectively) and the surface area of mature neurons (100 μm2). To account for the higher input resistance of young cells, the same current and conductance values were applied to a smaller surface of cells in networks with young intrinsic properties. The three-dimensional plots in Figure 7 indicate the density values for the mature cells. Coherence (κ) and frequency (f) values shown in Figure 7 are means of five simulations
Results
Morphological identification of young BCs
To obtain recordings from identified BCs in acute hippocampal slices, we had previously used transgenic mice expressing the enhanced GFP under the control of the PV promoter (Bartos et al., 2002; Meyer et al., 2002). However, this approach cannot be used early in postnatal development because PV expression begins at approximately P10–P14 (Du et al., 1996; Meyer et al., 2002; Tansey et al., 2002). To identify BCs at earlier developmental stages, we applied the following alternative strategy.
To identify BCs, we used a GAD67–GFP knock-in mouse line (Fig. 1) in which GABAergic cells are labeled early in development (Tamamaki et al., 2003). In acute hippocampal slices, several GFP-labeled cells were observed in the dentate gyrus in both young and mature mice (Fig. 1A,B). The density of GFP-labeled cells in the GC layer did not change significantly during postnatal development (P9, 7.9 ± 0.5 cells/5 × 105 μm3; P20, 6.1 ± 0.7 cells/5 × 105 μm3; p > 0.05; see Materials and Methods), consistent with previous data showing that cell death is almost absent in the dentate gyrus during postnatal development (Dupuy and Houser, 1997). We first confirmed that PV-positive interneurons were labeled in this mouse line. In mature mice, double labeling revealed that 14% of GFP-expressing cells colocalized PV in the dentate gyrus (39 from 274 GFP-labeled cells in 14 slices), in agreement with ∼15% of PV-immunopositive neurons within the GABAergic population in this hippocampal subregion (Kosaka and Hama, 1985). Furthermore, 43% of GFP-positive cells with somata located in the GC layer coexpressed PV. The remaining GFP-positive cells comprised interneuron types expressing other Ca2+-binding proteins, such as calbindin and calretinin, and cells colocalizing neuropeptides, such as cholecystokinin (J. A. Hosp, I. Vida, P. Jonas, and M. Bartos, unpublished observations). These results indicate that GAD67–GFP expression is suitable for identification of PV-positive BCs in the hippocampus. We then combined GFP labeling with light-microscopical analysis of the axonal localization of neurons filled with biocytin during recording. Neurons with an axonal arborization confined to the GC layer (>80%), a hallmark of BCs (Buhl et al., 1994; Freund and Buzsáki, 1996; Somogyi and Klausberger, 2005), were first observed around P6 (Fig. 1C,D). Thus, in early stages of postnatal development, BCs can be unequivocally identified on the basis of GAD67–GFP expression and location of the axon in the GC layer. At later stages of postnatal development, BCs were additionally identified immunocytochemically using an antibody against PV (Fig. 1E–G). At P12, five of five recorded BCs were PV positive. In summary, all BCs included in the present study fulfilled the criteria of GFP labeling and perisomatic location of the axon, and a subset of BCs additionally fulfilled the criterion of PV immunoreactivity.
Developmental changes in morphological and intrinsic membrane properties of BCs
We first examined the morphological properties of BCs at different developmental stages (Fig. 2). Quantitative morphological analysis of young (P6–P10) and mature (P18–P21) BCs after reconstruction using a camera lucida revealed that dendrites and axon continue to grow during postnatal maturation, resulting in a significant increase of 354 and 608%, respectively, in their total length (dendrites: 1521 ± 298 μm in four young cells vs 5387 ± 1015 μm in five mature cells, p < 0.05; axon: 3871 ± 1137 μm vs 23,522 ± 2971 μm, p < 0.001) (Fig. 2A–C). At the same time, the diameter of the dendrites (e.g., second-order apical dendrites; P6–P10, 3.4 ± 0.5 μm vs P18–P21, 4.1 ± 0.4 μm) and the axon (P6–P10, 1.3 ± 0.2 μm; P18–P21, 1.3 ± 0.3 μm) at a distance of 15–60 μm from the soma remained constant (p > 0.05). Axon collaterals of young BCs were primarily located (86 ± 7%; four cells) in the GC layer, comparable to those of the mature BCs (81 ± 7%; five cells). The growth of the axon was paralleled by a fivefold increase in the total number of branch points (p < 0.01) (supplemental Table 1, available at www.jneurosci.org as supplemental material) and increasing length of axonal segments (Fig. 2D). Thus, both the dendritic tree and the axonal arborization of BCs show substantial extension during development.
We next examined whether the observed morphological changes were accompanied by alterations in the passive membrane properties (Fig. 3). Consistent with an increase in membrane surface area (P6–P10, 4907.5 ± 634.2 μm2 vs P18–P21, 14,286.6 ± 3689 μm2; p < 0.05) (supplemental Table 1, available at www.jneurosci.org as supplemental material), the total membrane capacitance (cm) increased as a function of postnatal age from 46 ± 4 pF in young BCs (P6–P10, 21 cells) to 70 ± 4 pF in mature BCs (P18–P25, 16 cells; p < 0.001) (Fig. 3A). Furthermore, the input resistance (Rin) decreased substantially from 323 ± 38 MΩ (P6–P10) to 86 ± 6 MΩ (P18–P25; p < 0.001) (Fig. 3B). Finally, the membrane time constant (τm) depended on postnatal age. In young cells, τm was 21.9 ± 4.3 ms (P6–P10, 11 cells), significantly slower than in mature neurons with 7.3 ± 0.5 ms (P18–P25, 13 cells; p < 0.001) (Fig. 3D), suggesting changes in specific membrane properties. All passive membrane parameters changed approximately exponentially as a function of postnatal age.
A hallmark of mature BCs is the generation of brief action potentials and the ability to discharge at high frequency (>200 Hz at near-physiological temperature) in response to long depolarizing current pulses (Rudy and McBain, 2001). To test whether this spiking phenotype is also expressed in young BCs, we first analyzed the shape of somatically evoked single action potentials at various postnatal ages (Fig. 3E–H). The maximal rate of both rise and decay (dV/dtmax) increased significantly, resulting in a twofold reduction of the half-duration of single action potentials from 744 ± 56 μs in young BCs to 398 ± 36 μs in mature BCs (p < 0.001). Furthermore, the peak amplitude of action potentials increased as a function of postnatal age from 68 ± 2.6 mV at P6–P10 (14 cells) to 85.6 ± 1.9 mV at P18–P21 (9 cells; p < 0.05). Changes in action potential waveform were paralleled by a significant increase in current threshold for the initiation of single action potentials from 134.6 ± 18.5 pA in young BCs (13 cells) to 600 ± 44 pA in mature BCs (8 cells; p < 0.0005). Finally, young BCs fired high-frequency trains of action potentials when stimulated with long depolarizing current injections (Fig. 3I). The frequency–current relationship was steep for both young and mature cells (data not shown); however, maximal discharge frequencies were significantly lower in young neurons (P6–P10, 127.7 ± 8 Hz, 24 cells vs P18–P21, 263.3 ± 26 Hz, 9 cells; p < 0.0005) (Fig. 3J). To determine whether trains of action potentials evoked by long-lasting depolarizing current injections change their pattern during development, the ratio of the last and the first interspike interval at maximal discharge activity was calculated. The ratio did not change significantly during development (P6–P10, 1.4 ± 0.04, 9 cells; P18–P21, 1.2 ± 0.07, 9 cells; p > 0.1). In summary, both passive and active membrane properties of BCs are developmentally regulated, favoring slow signaling in young cells but rapid signaling in mature neurons.
Velocity and reliability of action potential propagation in BC axons
The observed developmental change in the maximal rate of rise of action potentials may lead to differences in the velocity of axonal action potential propagation between young and mature BCs. We tested this hypothesis by evoking single action potentials in BC axons by extracellular stimulation in the GC layer (Fig. 4). Antidromic action potentials were recorded at the soma under current-clamp conditions at a holding potential (Vhold) of −70 mV (Fig. 4A). In young BCs, the average conduction velocity (see Materials and Methods) was low, with a mean value of 0.14 m s−1 (four cells) (Fig. 4B). In contrast, in mature BCs, the conduction velocity was almost two times higher, with a mean value of 0.25 m s−1 (four cells). Additionally, a pronounced jitter in the onset of somatically recorded action potentials was evident in young BCs (Fig. 4C, left). In contrast, in mature BCs, antidromic action potentials were precisely timed (Fig. 4C, right). To quantify temporal precision, we used the CV of the conduction delay between the beginning of the stimulus artifact and the time corresponding to the maximal rate of rise of the somatically recorded action potential (Fig. 4D). On average, the CV of the conduction delay was four times higher in young (P6–P10) than in mature (P18–P25) BCs (0.08 ± 0.01 vs 0.02 ± 0.002, respectively; p < 0.01). Furthermore, there was a difference in the reliability of action potential propagation. To exclude stimulation failures, suprathreshold stimuli (∼1.3 times threshold) were used in these experiments. In young BCs, only 83.6 ± 0.9% of stimuli resulted in an action potential propagating to the soma, whereas in mature BCs, the reliability was 100% (Fig. 4D). Although we could not directly control for precise action potential initiation times in the axon, our data indicate that velocity, temporal precision, and reliability of action potential propagation in BC axons increased during postnatal development.
Precision and reliability of transmitter release at BC output synapses
To characterize developmental changes in GABAergic transmission at BC output synapses, we performed paired recordings from synaptically coupled BC–BC and BC–GC pairs at different postnatal ages (Fig. 5).
Young and mature BC output synapses differed in several properties of transmitter release. First, synaptic latency at BC–BC and BC–GC synapses decreased as a function of postnatal age (Fig. 5B). The latency–age relationship could be fitted exponentially for both BC–BC and BC–GC output synapses. Because the results were not significantly different between BC–BC and BC–GC pairs, data were pooled (Fig. 5E). Although the distance between presynaptic and postsynaptic recording sites was comparable (44 ± 3.8 μm), the mean synaptic latency was almost twofold longer at young (P6–P10) than at mature (P18–P25) BC output synapses (1.5 ± 0.2 ms, 15 pairs vs 0.9 ± 0.1 ms, 17 pairs, respectively; p < 0.005) (supplemental Table 2, available at www.jneurosci.org as supplemental material). This reduction in latency was in quantitative agreement with the increase in velocity of axonal action potential conduction (Fig. 4B), implying that the change in latency at BC–BC and BC–GC synapses was primarily due to the faster action potential propagation along the axon.
Second, at young BC–BC and BC–GC pairs, we observed a large jitter in synaptic latency of unitary IPSCs, suggesting imprecise GABA release (Fig. 5A). In contrast, at mature BC output synapses, latency jitter was substantially smaller, indicating that the temporal precision of GABA release increased during postnatal development. This change was most evident in the CV of the latency distribution that was used as a measure for temporal precision in transmitter release. The CV of latency decreased substantially at BC output synapses during development (p < 0.01, respectively) (Fig. 5C) and was almost three times higher at young than at mature BC synapses (0.28 ± 0.05, 12 pairs vs 0.1 ± 0.01, 10 pairs) (Fig. 5E). This is also reflected in a reduced half-duration of the latency distribution from 1.1 ms in young BC synapses (P6–P10, 12 pairs) to 0.31 ms in mature BC synapses (P18–P25, 10 pairs; BC–BC and BC–GC synapses, pooled data) (see Fig. 8). Because the CV of synaptic latency is larger for unitary IPSCs (Fig. 5C,E) than for antidromically propagated action potentials (Fig. 4D), an increase in temporal precision of both GABA release and action potential propagation contributes to the improved timing of inhibitory signaling at BC output synapses.
Finally, probability of failures in GABAergic transmission decreased markedly during postnatal development at BC–BC synapses (young, 32 ± 11% vs mature, 4.5 ± 3%; p < 0.05) (supplemental Table 2, available at www.jneurosci.org as supplemental material) and BC–GC synapses (young, 42 ± 5% vs mature, 2 ± 2%; p < 0.005) (Fig. 5D,E). Failures at output synapses of young BCs can be largely explained by the stochastic properties of transmitter release and partially (∼16%) by the unreliable action potential propagation along the axon.
In contrast to the marked changes in chemical transmission during development, we observed significant changes neither in the incidence of electrical coupling among BCs (P6–P14, 9.5%; P18–P25, 11%) nor in the ratio of electrical coupling at BC–BC pairs (P6–P14: 0.013 ± 0.01, four pairs; P18–P25: 0.011 ± 0.01, five pairs), consistent with previous reports in the neocortex (Long et al., 2005). In summary, BC output synapses undergo marked developmental changes in the properties of transmitter release, leading to a developmental shift from imprecise and unreliable to precise and reliable signaling devices.
Increased peak amplitude and accelerated decay time course of unitary IPSCs
Light-microscopic analysis of intracellularly labeled BC–BC and BC–GC pairs showed that putative output synapses of both young (Fig. 1) and mature BCs (Geiger et al., 1997; Bartos et al., 2001) are located close to the soma of their target cells. This allowed us to measure the postsynaptic inhibitory conductance change as IPSCs in the voltage-clamp configuration at various postnatal ages (Fig. 6).
During maturation, a fourfold increase in the peak amplitude of unitary IPSCs (Fig. 6B) was observed at both BC–BC (P6–P10, 27 ± 5 pA, 4 pairs vs P18–P25, 104 ± 24 pA, 7 pairs; p < 0.01) and BC–GC (P6–P10, 37 ± 19 pA, 11 pairs vs P18–P25, 148 ± 34 pA, 10 pairs; p < 0.01) synapses (supplemental Table 2, available at www.jneurosci.org as supplemental material). This increase in synaptic strength was paralleled by a growing probability of obtaining paired BC–BC recordings from 7% at P6–P8 to 15% at P18–P25 (Fig. 6C). Similarly, the connection probability for BC–GC synapses increased from 10% at P6–P8 to 55% at P18–P25. These data together indicate that both the number of target cells and the number of contact sites per target cell increase during maturation, consistent with the observed extension of the axonal arborization. Paired recording analysis also revealed a marked decrease in the amplitude-weighted decay time constant (τw) of unitary IPSCs at both BC–BC (P6–P10, 7.5 ± 1.5 ms vs P18–P25, 2.0 ± 0.2 ms; p < 0.0005) and BC–GC (P6–P10, 8.8 ± 1.4 ms vs P18–P25, 4.4 ± 0.5 ms; p < 0.005) synapses (Fig. 6B). Intriguingly, the relative differences between BC–BC and BC–GC synapses were small in the young dentate gyrus but became greatly accentuated during development (p < 0.005 at P18–P25). Thus, the maturation of the fast kinetics of inhibition is target-cell specific.
It has been suggested that the α1 subunit is a key molecular determinant of rapid kinetics of recombinant GABAA receptors (Lavoie et al., 1997; Cope et al., 2005) and that the postsynaptic sites of BC–BC synapses are highly enriched in α1 subunit immunoreactivity (Fritschy and Möhler, 1995; Nusser et al., 1995; Klausberger et al., 2002). To probe the α1 content of BC–BC synapses in young and mature hippocampus, we tested the effects of the GABAA receptor modulator Zolpidem (Fig. 6D). At low concentrations, Zolpidem selectively affects GABAA receptors assembled from α1βγ2 subunits (Thomson et al., 2000; Cope et al., 2005). The addition of 5 μm Zolpidem to the bath solution significantly increased τw of compound IPSCs evoked by extracellular stimulation in the GC layer in mature cells (159.6 ± 36%; six cells), indicating the presence of α1 subunits. In contrast, no significant effect was observed on τw of IPSCs recorded in young neurons (8.5 ± 16%; five cells), suggesting that their synapses are devoid of α1 subunits. Thus, the developmental acceleration of the unitary IPSC decay time course at BC–BC synapses can be explained by a change in the subunit composition, especially a higher α1 content of postsynaptic GABAA receptors.
Generation of gamma oscillations in young interneuron network models
The functional properties of mature BCs appear to be optimized for the generation of synchronous gamma frequency oscillations in interneuron networks (Bartos et al., 2001, 2002, 2007; Vida et al., 2006). However, the present results indicate that the properties of young BCs are very different from those of mature cells. Can interneuron networks with young intrinsic and synaptic properties support synchronization at gamma frequencies?
To answer this question, we modified stepwise our previously developed interneuron network model (Vida et al., 2006) to reproduce salient experimentally determined properties of young BCs and examined the effects of these changes on coherence and frequency of network oscillations (Fig. 7). As shown previously, the mature interneuron network with fast action potential conduction (0.25 m s−1), rapid IPSCs (τw = 2 ms), and high connectivity (mean number of postsynaptic neurons, Msyn = 60) (Sik et al., 1995) generated highly coherent oscillations in the upper gamma frequency band (maximal κ = 0.74 at a network frequency of 56 Hz) when activated by a low, moderately heterogeneous tonic excitatory drive (Iμ = 0.5–2.5 μA cm−2; Iσ /Iμ = 0.1) (Fig. 7A,B) [Vida et al. (2006), their Fig. 3B]. Introduction of young intrinsic properties to the mature network, including slow membrane time constant, broad action potential waveform, and low conduction velocity (0.13 m s−1), resulted in a reduction in both network coherence and frequency (maximal κ = 0.7 at 43 Hz, Fig. 7C). Similarly, when the time course of IPSCs was slowed to match the measured values in young BCs (τw = 8 ms), both peak coherence and network frequency were reduced (maximal κ = 0.68 at 39 Hz) (Fig. 7D). When connectivity was reduced in the mature network (Msyn = 40), the peak coherence dropped markedly (maximal κ = 0.57), whereas the corresponding network frequency remained high (56 Hz) (Fig. 7E). Interestingly, when slow intrinsic properties and slow IPSCs were combined in the networks, gamma oscillations with high coherence but low frequency emerged (maximal κ = 0.89 at 34 Hz) (Fig. 7F). In contrast, when low connectivity was combined to the slow intrinsic properties or slow IPSCs, coherence was further reduced (0.39 at 55 Hz and 0.56 at 44 Hz, respectively) (Fig 7G,H). Finally, when all three properties were introduced together, to reproduce an immature interneuron network, oscillations were generated with low coherence in the lower gamma frequency band (maximal κ = 0.49 at 34 Hz) (Fig. 7I,J).
In summary, our data show that intrinsic, synaptic, and structural properties of interneuron networks are important parameters that jointly define coherence and frequency of network oscillations. Changes in the intrinsic properties and the kinetics of ISPCs act synergistically to promote oscillations, first in the lower and later in the upper gamma frequency band during postnatal maturation. In contrast, increasing network connectivity promotes synchrony in the neuron network and thereby underlies the increasing power of oscillations as observed in the developing brain (Lahtinen et al., 2002).
Discussion
Dentate gyrus BCs undergo marked changes in their morphological, intrinsic, and synaptic properties at P6–P25 (Fig. 8). First, membrane time constant, action potential waveform, and axonal conduction velocity become faster. Second, precision, strength, and reliability of transmitter release increase at BC output synapses. Finally, kinetics of unitary IPSCs become rapid. Thus, concerted changes, presumably driven by a complex genetic program, convert BCs from slow to fast signaling devices during the late postnatal period. Our computational analysis also reveals that combined changes in cellular and synaptic properties of BCs can support the developmental increase in frequency and coherence of gamma oscillations observed in vivo (Lahtinen et al., 2002).
Developmental changes in BC morphology
Size and functional properties of BC dendrites change substantially during development. The total dendritic length increases 3.5-fold at P6–P21 (Fig. 2). In correlation with the associated enlargement in surface area, we found an increase in the apparent membrane capacitance (cm) and a decrease in input resistance (Rin) and membrane time constant (τm) (Fig. 3A–D). Because the capacitance of lipid membranes is relatively constant (∼1 μF cm−2) (Hille, 2001), the latter changes indicate a decrease in specific membrane resistance (Rm). Our results therefore suggest major changes in the extent and the integrative properties of BC dendrites.
The developmental reduction in Rm could be explained by increases in densities of leak channels (e.g., two-pore domain K+ channels) (Taverna et al., 2005; Torborg et al., 2006). Changes in morphology and membrane properties jointly redefine somatodendritic integration of EPSPs in BCs. In young BCs, the size of the dendritic tree is smaller, Ri is higher, and τm is slower than in mature BCs. The combination of these properties favors spatial and temporal integration of synaptic inputs in young BCs but coincidence detection in mature BCs.
Substantial changes during development also occur in BC axons. The total length of the axon increases sixfold at P6–P21 (Fig. 2). In parallel, the probability of obtaining paired recordings increases approximately twofold (Fig. 6C). Thus, axon maturation will lead to the innervation of a higher number of target cells in a larger area. Within the same developmental window, velocity and reliability of axonal action potential propagation also increases (Fig. 4B,D), supporting rapid and reliable synaptic transmission, a prerequisite for the function of mature BC synapses (Bartos et al., 2001, 2002). Both structural and functional changes could explain the increase in conduction velocity and synchrony of transmitter release (Williams et al., 1998). Faster action potential propagation may be caused by a larger axon diameter; however, we found comparable axon diameters in young and mature cells. Alternatively, an enhanced density of voltage-dependent Na+ channels in BC axons may be responsible (Jack et al., 1975). Consistent with this hypothesis, we found that maximal rate of rise of axonal action potentials is twofold faster in mature than in young BCs (Fig. 3G). An increasing Na+ channel density in developing BCs not only explains the observed change in conduction velocity but also the increased reliability of action potential propagation. This effect may further be supported by increases in myelination of BC axons. However, this hypothesis needs additional experimental proof. Collectively, developmental changes in BC axons lead to a situation in which BCs innervate a higher number of target cells in a larger area with shorter delays and increased reliability.
Developmental changes at BC output synapses
GABAergic transmission at BC output synapses is characterized by a threefold higher synchrony at mature than at young synapses (Fig. 5B,C,E). The CV of synaptic latency is much larger than that of action potential conduction velocity in young BCs, indicating that the lower synchrony can only partly be explained by variable and slow spike propagation.
Several mechanisms may contribute to the enhanced synchronization of transmitter release. First, a reduction in presynaptic action potential duration may occur. Indeed, we found a substantial developmental shortening of somatic action potentials in BCs (Fig. 3H). Similar changes at presynaptic sites would explain increases in synchrony of transmitter release. Consistent with this hypothesis, the density of Kv3-type potassium channels, a major determinant of action potential repolarization in interneurons (Rudy and McBain, 2001), increases during development (Du et al., 1996; Tansey et al., 2002). Second, properties of presynaptic Ca2+ channels may change developmentally. Indeed, transmitter release at BC output synapses depends on N- and P/Q-type channels in young BCs (M. Bartos, unpublished observations) but relies exclusively on P/Q-type channels in mature cells (Hefft and Jonas, 2005). Coupling of these channels to the release machinery appears to be different. P/Q type are coupled more tightly to the release machinery than N-type channels, leading to sharper Ca2+ transients at the Ca2+ sensor (Bucurenciu et al., 2008). Finally, it is possible that properties of the exocytotic machinery change. For example, the expression of synaptotagmins, the putative Ca2+ sensors, is developmentally regulated (Berton et al., 1997). Thus, different presynaptic changes may converge to synchronize transmitter release. Indeed, combined presynaptic alternations underlying synchronization of transmitter release during development were reported at the calyx of Held (Taschenberger and von Gersdorff, 2000; Yang and Wang, 2006).
Additional developmental changes take place at postsynaptic sites of BC output synapses. The decay time constant of unitary IPSCs is approximately fourfold faster at mature than at young synapses (Fig. 6). How can we explain these changes? In situ hybridization and immunohistochemical data suggest changes in the subunit composition of GABAA receptors during development. Although α2 subunits are highly expressed at early stages, α1 subunits are upregulated at later times (Laurie et al., 1992; Fritschy and Möhler, 1995). Immunocytochemical analysis revealed that mature BC output synapses are highly enriched in α1 subunits (Klausberger et al., 2002). GABAA receptors containing α1 deactivate faster than those with α2 (Lavoie et al., 1997). Thus, a developmental switch in subunit composition of postsynaptic GABAA receptors could explain changes in IPSC kinetics. This hypothesis is consistent with our result that IPSCs are sensitive to the α1-selective benzodiazepine agonist Zolpidem at mature but not at young BC synapses (Fig. 6D).
Implications for network function
Developmental alterations in BC properties have several implications for neuronal network function. First, we predict that young and mature BCs are differentially involved in feedforward and feedback inhibition. Young BCs integrate asynchronous inputs, generate low-frequency trains of action potentials (this study; Daw et al., 2007), and elicit weak and slow IPSCs in their target neurons. Thus, young BCs generate a slow form of inhibition. In contrast, mature BCs detect coincident excitatory inputs, generate high-frequency trains of action potentials (present results; Geiger et al., 1997; Galarreta and Hestrin, 2001; Penttonen et al., 1998), and generate a fast inhibitory synaptic output, providing phasic inhibition with high temporal precision to their targets (Hefft and Jonas, 2005). Thus, fast signaling properties of mature BCs contribute to speed and temporal precision of feedforward and feedback inhibition (Pouille and Scanziani, 2001, 2004; Galarreta and Hestrin, 2002). These developmental changes may account for the sharpening of place fields of GCs and the emergence of critical periods of synaptic plasticity during development (Martin and Berthoz, 2002; Hensch, 2005).
Developmental changes in BC signaling are relevant for the generation of network oscillations. Experimental and theoretical evidence indicate that mature BCs play a key role in the generation of robust gamma oscillations in vitro and in vivo (Bragin et al., 1995; Whittington et al., 1995; Traub et al., 1996, 2004; Wang and Buzsáki, 1996; Fisahn et al., 1998; Towers et al., 2002; Csicsvari et al., 2003). BCs form networks via mutual inhibitory synapses, and the decay time constant of IPSCs at BC synapses has been identified as key parameter setting frequency and coherence of network oscillations (M. Bartos, unpublished observations; Traub et al., 1996; Wang and Buzsáki, 1996; Cope et al., 2005). During postnatal development, gamma activity emerges at the end of the first postnatal week (Lahtinen et al., 2002; Leinekugel et al., 2002). However, frequency and coherence of the oscillations are low at this stage and reach adult levels only by the fourth postnatal week. This time course corresponds well to the appearance and subsequent functional maturation of BC-mediated inhibition observed in this study. Our simulations reveal that neuronal networks with immature properties generate gamma oscillations (Fig. 7), but with markedly lower frequency and coherence than those with mature characteristics (Bartos et al., 2002; Vida et al., 2006). Whereas changes in intrinsic and synaptic properties of BCs set the frequency of the oscillations, changes in connectivity can explain the developmental increase in coherence of gamma activity observed in vivo (Lahtinen et al., 2002). Although our network model reproduces adequately gamma oscillations in young (Lahtinen et al., 2002) and mature (Bragin et al., 1995; Penttonen et al., 1998) animals in vivo, the next challenge is to determine how interactions among various neuron types contribute to developmental changes in network dynamics.
In summary, our results show that, during late postnatal development, BCs differentiate functionally from slow to fast signaling devices. Changes in morphological, intrinsic and synaptic properties of BCs associated with this process are crucial for the emergence of robust and highly coherent gamma oscillations and may thereby contribute to the maturation of higher cognitive brain functions in both rodents and humans (Ben-Ari et al., 2007; Yurgelun-Todd, 2007).
Footnotes
-
This work was supported by the Deutsche Forschungsgemeinschaft (SFB 505/C6, SFB 780/A5), the Bundesministerium für Bildung und Forschung (OI GQ 0420), the Volkswagenstiftung (I/78563-564), the Royal Society (M.B.), and by Grants-in-Aid for Scientific Research from Ministry of Education, Culture, Sports, Science, and Technology of Japan (Y.Y.). We thank W. Wisden for critically reading this manuscript. We thank J. Sauer for the quantification of GAD67/GFP-labeled cells in the dentate gyrus. We also thank K. Winterhalter, S. Becherer, M. Northemann, and U. Nöller for technical assistance.
- Corresponding should be addressed to Dr. Marlene Bartos, Institute of Medical Sciences, University of Aberdeen, Foresterhill, Aberdeen AB25 2ZD, UK. m.bartos{at}abdn.ac.uk