-
PDF
- Split View
-
Views
-
Cite
Cite
Christian Wozny, Stephen R. Williams, Specificity of Synaptic Connectivity between Layer 1 Inhibitory Interneurons and Layer 2/3 Pyramidal Neurons in the Rat Neocortex, Cerebral Cortex, Volume 21, Issue 8, August 2011, Pages 1818–1826, https://doi.org/10.1093/cercor/bhq257
- Share Icon Share
Abstract
Understanding the structure and function of the neocortical microcircuit requires a description of the synaptic connectivity between identified neuronal populations. Here, we investigate the electrophysiological properties of layer 1 (L1) neurons of the rat somatosensory neocortex (postnatal day 24–36) and their synaptic connectivity with supragranular pyramidal neurons. The active and passive properties of visually identified L1 neurons (n = 266) suggested division into 4 groups according to the Petilla classification scheme with characteristics of neurogliaform cells (NGFCs) (n = 72), classical-accommodating (n = 137), fast-spiking (n = 23), and burst-spiking neurons (n = 34). Anatomical reconstructions of L1 neurons supported the existence of 4 major neuronal groups. Multiparameter unsupervised cluster analysis confirmed the existence of 4 groups, revealing a high degree of similarity with the Petilla scheme. Simultaneous recordings between synaptically connected L1 neurons and L2/3 pyramidal neurons (n = 384) demonstrated neuronal class specificity in both excitatory and inhibitory connectivity and the properties of synaptic potentials. Notably, all groups of L1 neurons received monosynaptic excitatory input from L2/3 pyramidal neurons (n = 33), with the exception of NGFCs (n = 68 pairs tested). In contrast, NGFCs strongly inhibited L2/3 pyramidal neurons (n = 12 out 27 pairs tested). These data reveal a high specificity of excitatory and inhibitory connections in the superficial layers of the neocortex.
Introduction
The neocortex is a 6-layered structure which is functionally organized into columns (Mountcastle 1997). Understanding the wiring diagram of a cortical column has received much attention in recent years (Douglas and Martin 2004). The excitatory synaptic circuitry of the cortical column has been well described (Lubke and Feldmeyer 2007; Schubert et al. 2007). The neocortex, however, is not only composed of excitatory neurons but also γ-aminobutyric acid (GABA)ergic inhibitory interneurons, which make up ∼20% of the neocortical neuronal population. Inhibitory interneurons can be classified by morphological, electrophysiological, neurochemical features (Ascoli et al. 2008), and their transcriptional regulation during development (Butt et al. 2007; Batista-Brito and Fishell 2009). The distribution of interneuron classes, however, differs across neocortical layers (Markram et al. 2004; Ascoli et al. 2008).
Layer 1 (L1) of the neocortex is a cell-sparse synaptic and axon dense zone. In contrast to other neocortical layers, L1 is composed of more than 90% of GABAergic neurons (Winer and Larue 1989; Prieto et al. 1994). In the primary somatosensory cortex, whisker-evoked sensory information is rapidly relayed to L1 neurons (Zhu Y and Zhu JJ 2004), which, in turn, act to powerfully inhibit sensory-evoked responses in L2/3 (Shlosberg et al. 2006). Surprisingly, the interneuronal composition and synaptic connectivity of L1 neurons have not been explored in detail (Hestrin and Armstrong 1996; Zhou and Hablitz 1996a, 1996b; Chu et al. 2003; Zhu Y and Zhu JJ 2004). Layer 1 constitutes an important site of integration in the neocortical column as it contains a dense plexus of the apical dendrites of layer 2/3 and layer 5 pyramidal neurons and axons originating from a variety of cortical areas (Cauller et al. 1998; Mitchell and Cauller 2001; Petreanu et al. 2009) and the thalamus (Rubio-Garrido et al. 2009). Long-range intracortical excitatory axons, conveyed in L1, synapse with both the apical dendrites of pyramidal neurons (Petreanu et al. 2009) and local GABAergic interneurons (Anderson and Martin 2006). Thus, L1 inhibitory interneurons are ideally placed to control dendritic synaptic integration in pyramidal neurons (Williams and Stuart 2002; Waters et al. 2003; Williams 2004; Larkum et al. 2009).
In this study, we have characterized the electrophysiological properties of L1 neurons of the rat neocortex and investigated the synaptic connectivity between L1 inhibitory interneurons and L2/3 pyramidal neurons. Our results provide new insights into the wiring diagram of the superficial layers of the neocortical column.
Materials and Methods
Slices and Recordings
Wistar rats (P24 to P36) were decapitated under deep isoflurane anesthesia following UK Home office and institutional guidelines. Brains were quickly removed and placed into ice-cold solution containing (mM): NaCl 125, NaHCO3 25, KCl 3, NaH2PO4 1.25, CaCl2 1, MgCl2 6, sodium pyruvate 3, and glucose 25, oxygenated at 95% O2 and 5% CO2. After cutting coronal sections (300–350 μm) of the somatosensory cortex (slices were collected corresponding to the coordinates of Bregma ca. –3.8 mm to Bregma ca. –1.8 mm Paxinos and Watson 1998), slices were incubated at 34 °C for 30 min and then stored at room temperature. For recordings, a single brain slice was placed in a chamber perfused with a solution of composition (mM): NaCl 125, NaHCO3 25, KCl 3, NaH2PO4 1.25, CaCl2 2, MgCl2 1, sodium pyruvate 3, and glucose 25 at 34–36 °C. Double and triple whole-cell recordings from the somata of layer 1 (L1) and layer 2/3 (L2/3) somatosensory neocortical neurons were made with identical current clamp amplifiers (BVC 700A; Dagan). The L1–L2 border was defined by the abrupt change in the density of neurons visualized under infrared differential interference contrast video microscopy. Pipettes were filled with a solution containing (mM): potassium gluconate 135, NaCl 7, Hepes 10, Na2-ATP 2, Na-GTP 0.3, MgCl2 2, and occasionally with Alexa Fluor 568 0.01–0.04 (Molecular Probes) for visualization after recording; pH was adjusted to 7.2–7.3 with KOH.
Intrinsic Properties of L1 and L2/3 Neurons
The action potential (AP) firing pattern of neurons was investigated by the injection of a series of 600 ms positive current steps through the recording pipette (100 pA steps up to a maximum of 1 nA). Input resistance (Rin) was calculated from the steady-state voltage change evoked by the injection of a 100 pA negative current step. Resting membrane potential (RMP) was measured at the beginning of each recording and was not corrected for the liquid junction potential. The accommodation of a neuron's AP firing frequency was assessed by calculating an accommodation index, which is defined as the ratio of the first and last interspike interval at threshold for repetitive AP firing. We did not estimate the rheobase of neurons as we injected current through the patch pipette in 100 pA steps and therefore the accommodation index was calculated at near-threshold voltages.
Cluster Analysis
To cluster L1 neurons into groups, unsupervised clustering was performed using Ward's method (WinSTAT) (for details, see Ward 1963; Cauli et al. 2000). The results of the cluster analysis were compared with the classification of neurons according to the Petilla terminology (see Results; Ascoli et al. 2008).
Synaptic Connectivity
Excitatory and inhibitory synaptic transmission was examined between simultaneously recorded L2/3 pyramidal neurons and L1 neurons. Synaptic connectivity was tested for by the repeated generation of a pair of APs (evoked by 2 ms test pulses; 0.5–5 nA; separated by 50 ms; delivered at 0.33 Hz). If neurons were not synaptically connected, one of the recording electrodes was withdrawn and either a new pyramidal neuron or interneuron recorded. Voltage and current signals were low-pass filtered at 10 kHz and acquired at 25–50 kHz using an ITC-18 interface (HEKA) controlled by an Apple computer. Data were analyzed and curve fitting performed using Axograph X. Unitary postsynaptic potential (uPSP) onset time, rise time, and half-width was measured from either single episodes or from averages (usually 50 responses). Onset was determined at 5% of peak response, rise time was measured as the time to reach between 10 and 90% of the maximal amplitude of the uPSP, and the decay time constants were determined by fitting the decay with a single exponential function. To accurately measure the amplitude of the second uPSP following a paired-pulse stimulus, the decay of first uPSP was fitted with an exponential function, which was subsequently subtracted from the trace. The amplitude of the second uPSP was measured over a window centered around peak amplitude relative to the subtracted trace. Statistical analysis was performed with Student's t-test (Excel, Microsoft), chi-squared test, Fisher's Exact test, or with an analysis of variance (ANOVA) test (all Sigma-Plot).
Reconstructions
Neurons were filled with biocytin (0.5%, added to the pipette recording solution) and slices fixed overnight in 4% paraformaldehyde dissolved in 0.1 M sodium phosphate-buffered saline (PBS, pH 7.4). Subsequently, slices were incubated for 12–18 h in PBS supplemented with 0.5–2% Triton X-100 and 0.2% streptavidin, Alexa Fluor 488 conjugate (Invitrogen) at 4–6 °C. After washing with PBS, slices were mounted on slides, covered with Vectashield and coverslipped (Vector Laboratories). Confocal images were taken with Zeiss 510 or Zeiss 710 confocal microscope, respectively. Neuronal morpholology was reconstructed from confocal image stacks using Neurolucida software (MBF Bioscience).
Results
Electrophysiological Properties of Layer 1 Neurons
We recorded from 271 layer 1 (L1) neurons located at least 20 μm from the L1–L2 border, selected irrespective of their morphology under infrared differential interference contrast video microscopy, using whole-cell current clamp techniques. Despite this selection criteria, 5 of 271 neurons had electrophysiological properties similar to those of simultaneously recorded L2/3 pyramidal neurons (n = 384; low input resistance and prominent rectification of the voltage–current relationship, see Supplementary Fig. 1). We first classified the remaining 266 L1 neurons according to their repetitive AP firing patterns in response to suprathreshold injection of positive current steps (Fig. 1). We observed key differences between groups of neurons in the initial AP firing frequency, the half-width of the first AP, and accommodation index (Fig. 1; see Materials and Methods). According to the Petilla classification scheme (Ascoli et al. 2008), we observed regular- (n = 209), fast- (n = 23), and burst-spiking (BS) neurons (n = 34). Regular-spiking behavior was observed in a group of neurons termed classical-accommodating neurons (c-AC, n = 137) and in a group resembling neurogliaform cells (NGFCs, n = 72) (Tamas et al. 2003; Olah et al. 2007; Szabadics et al. 2007). Most of the fast-spiking (FS) neurons were observed at the L1-L2 boarder (see anatomical features below). BS neurons showed high initial (burst) frequency (Fig. 1A,B). We divided neuronal classes according to their accommodation index (see Materials and Methods), which separated c-AC from NGFC cells (Fig. 1C) and the half-width of the first AP, which was narrowest in FS cells (Fig. 1D, Table 1). These quantitative measures allowed the clear division of L1 neurons into 4 groups.
RMP (mV) | Rin (MΩ) | Firing Freqinitial.(Hz) | AP1 half-width (ms) | Accommodation index (t(AP2 − AP1)/t(APn−1 − APn)) | |
Petilla | |||||
NGFC (72) | −73.4 ± 0.5 | 107 ± 3 | 23.0 ± 1.6 | 0.48 ± 0.01 | 1.12 ± 0.04 |
FS (23) | −72.0 ± 1.1 | 87 ± 7 | 93.0 ± 12.5 | 0.28 ± 0.01 | 0.93 ± 0.1 |
c-AC (137) | 67.3 ± 0.4 | 166 ± 5 | 44.5 ± 4.2 | 0.5 ± 0.01 | 0.50 ± 0.02 |
BS (34) | −69.2 ± 0.9 | 201 ± 15 | 285.9 ± 11.3 | 0.43 ± 0.01 | 0.09 ± 0.01 |
Clustering | |||||
Cluster 1 (81) (81% NGFC) | −74.0 ± 0.4 (−74.1 ± 0.5) | 108 ± 3 (106 ± 3) | 23.8 ± 2.3 (21.7 ± 1.5) | 0.49 ± 0.1 (0.48 ± 0.01) | 1.06 ± 0.04 (1.14 ± 0.04) |
Cluster 2 (24) (96% FS) | −72.0 ± 1.1 (−72.0 ± 1.1) | 88 ± 7 (87 ± 8) | 92.8 ± 12.0 (93.0 ± 12.5) | 0.28 ± 0.01 (0.28 ± 0.01) | 0.91 ± 0.1 (0.93 ± 0.1) |
Cluster 3 (125) (94% c-AC) | −66.2 ± 0.4 (−66.4 ± 0.4) | 167 ± 4 (169 ± 4) | 41.2 ± 3.2 (39.7 ± 3.2) | 0.50 ± 0.01 (0.50 ± 0.01) | 0.51 ± 0.02 (0.49 ± 0.02) |
Cluster 4 (36) (89% BS) | −69.7 ± 0.7 (−69.8 ± 0.8) | 209 ± 17 (202 ± 16) | 286.3 ± 11.3 (295.2 ± 9.9) | 0.43 ± 0.01 (0.43 ± 0.01) | 0.09 ± 0.01 (0.08 ± 0.01) |
RMP (mV) | Rin (MΩ) | Firing Freqinitial.(Hz) | AP1 half-width (ms) | Accommodation index (t(AP2 − AP1)/t(APn−1 − APn)) | |
Petilla | |||||
NGFC (72) | −73.4 ± 0.5 | 107 ± 3 | 23.0 ± 1.6 | 0.48 ± 0.01 | 1.12 ± 0.04 |
FS (23) | −72.0 ± 1.1 | 87 ± 7 | 93.0 ± 12.5 | 0.28 ± 0.01 | 0.93 ± 0.1 |
c-AC (137) | 67.3 ± 0.4 | 166 ± 5 | 44.5 ± 4.2 | 0.5 ± 0.01 | 0.50 ± 0.02 |
BS (34) | −69.2 ± 0.9 | 201 ± 15 | 285.9 ± 11.3 | 0.43 ± 0.01 | 0.09 ± 0.01 |
Clustering | |||||
Cluster 1 (81) (81% NGFC) | −74.0 ± 0.4 (−74.1 ± 0.5) | 108 ± 3 (106 ± 3) | 23.8 ± 2.3 (21.7 ± 1.5) | 0.49 ± 0.1 (0.48 ± 0.01) | 1.06 ± 0.04 (1.14 ± 0.04) |
Cluster 2 (24) (96% FS) | −72.0 ± 1.1 (−72.0 ± 1.1) | 88 ± 7 (87 ± 8) | 92.8 ± 12.0 (93.0 ± 12.5) | 0.28 ± 0.01 (0.28 ± 0.01) | 0.91 ± 0.1 (0.93 ± 0.1) |
Cluster 3 (125) (94% c-AC) | −66.2 ± 0.4 (−66.4 ± 0.4) | 167 ± 4 (169 ± 4) | 41.2 ± 3.2 (39.7 ± 3.2) | 0.50 ± 0.01 (0.50 ± 0.01) | 0.51 ± 0.02 (0.49 ± 0.02) |
Cluster 4 (36) (89% BS) | −69.7 ± 0.7 (−69.8 ± 0.8) | 209 ± 17 (202 ± 16) | 286.3 ± 11.3 (295.2 ± 9.9) | 0.43 ± 0.01 (0.43 ± 0.01) | 0.09 ± 0.01 (0.08 ± 0.01) |
Petilla: The active and passive properties of L1 neurons grouped according to the Petilla classification scheme (c-AC, NGFC, FS, and BS). Number of neurons are given in brackets. Clustering: unsupervised clustering (hierarchical clustering using Ward's method) was performed and a cut off criteria was determined to reveal 4 clusters. Number of neurons in each cluster are given in brackets. The primary composition of neurons in each cluster when identified by the Petilla classification scheme are indicated: 81% of neurons in cluster 1 were NGFCs, 96% of neurons in cluster 2 FS neurons, 94% of neurons in cluster 3 c-AC neurons, and 89% of neurons in cluster 4 BS neurons. The average values for the active and passive properties for the primary composition of neurons in each cluster. All values are given as mean ± standard error of the mean.
RMP (mV) | Rin (MΩ) | Firing Freqinitial.(Hz) | AP1 half-width (ms) | Accommodation index (t(AP2 − AP1)/t(APn−1 − APn)) | |
Petilla | |||||
NGFC (72) | −73.4 ± 0.5 | 107 ± 3 | 23.0 ± 1.6 | 0.48 ± 0.01 | 1.12 ± 0.04 |
FS (23) | −72.0 ± 1.1 | 87 ± 7 | 93.0 ± 12.5 | 0.28 ± 0.01 | 0.93 ± 0.1 |
c-AC (137) | 67.3 ± 0.4 | 166 ± 5 | 44.5 ± 4.2 | 0.5 ± 0.01 | 0.50 ± 0.02 |
BS (34) | −69.2 ± 0.9 | 201 ± 15 | 285.9 ± 11.3 | 0.43 ± 0.01 | 0.09 ± 0.01 |
Clustering | |||||
Cluster 1 (81) (81% NGFC) | −74.0 ± 0.4 (−74.1 ± 0.5) | 108 ± 3 (106 ± 3) | 23.8 ± 2.3 (21.7 ± 1.5) | 0.49 ± 0.1 (0.48 ± 0.01) | 1.06 ± 0.04 (1.14 ± 0.04) |
Cluster 2 (24) (96% FS) | −72.0 ± 1.1 (−72.0 ± 1.1) | 88 ± 7 (87 ± 8) | 92.8 ± 12.0 (93.0 ± 12.5) | 0.28 ± 0.01 (0.28 ± 0.01) | 0.91 ± 0.1 (0.93 ± 0.1) |
Cluster 3 (125) (94% c-AC) | −66.2 ± 0.4 (−66.4 ± 0.4) | 167 ± 4 (169 ± 4) | 41.2 ± 3.2 (39.7 ± 3.2) | 0.50 ± 0.01 (0.50 ± 0.01) | 0.51 ± 0.02 (0.49 ± 0.02) |
Cluster 4 (36) (89% BS) | −69.7 ± 0.7 (−69.8 ± 0.8) | 209 ± 17 (202 ± 16) | 286.3 ± 11.3 (295.2 ± 9.9) | 0.43 ± 0.01 (0.43 ± 0.01) | 0.09 ± 0.01 (0.08 ± 0.01) |
RMP (mV) | Rin (MΩ) | Firing Freqinitial.(Hz) | AP1 half-width (ms) | Accommodation index (t(AP2 − AP1)/t(APn−1 − APn)) | |
Petilla | |||||
NGFC (72) | −73.4 ± 0.5 | 107 ± 3 | 23.0 ± 1.6 | 0.48 ± 0.01 | 1.12 ± 0.04 |
FS (23) | −72.0 ± 1.1 | 87 ± 7 | 93.0 ± 12.5 | 0.28 ± 0.01 | 0.93 ± 0.1 |
c-AC (137) | 67.3 ± 0.4 | 166 ± 5 | 44.5 ± 4.2 | 0.5 ± 0.01 | 0.50 ± 0.02 |
BS (34) | −69.2 ± 0.9 | 201 ± 15 | 285.9 ± 11.3 | 0.43 ± 0.01 | 0.09 ± 0.01 |
Clustering | |||||
Cluster 1 (81) (81% NGFC) | −74.0 ± 0.4 (−74.1 ± 0.5) | 108 ± 3 (106 ± 3) | 23.8 ± 2.3 (21.7 ± 1.5) | 0.49 ± 0.1 (0.48 ± 0.01) | 1.06 ± 0.04 (1.14 ± 0.04) |
Cluster 2 (24) (96% FS) | −72.0 ± 1.1 (−72.0 ± 1.1) | 88 ± 7 (87 ± 8) | 92.8 ± 12.0 (93.0 ± 12.5) | 0.28 ± 0.01 (0.28 ± 0.01) | 0.91 ± 0.1 (0.93 ± 0.1) |
Cluster 3 (125) (94% c-AC) | −66.2 ± 0.4 (−66.4 ± 0.4) | 167 ± 4 (169 ± 4) | 41.2 ± 3.2 (39.7 ± 3.2) | 0.50 ± 0.01 (0.50 ± 0.01) | 0.51 ± 0.02 (0.49 ± 0.02) |
Cluster 4 (36) (89% BS) | −69.7 ± 0.7 (−69.8 ± 0.8) | 209 ± 17 (202 ± 16) | 286.3 ± 11.3 (295.2 ± 9.9) | 0.43 ± 0.01 (0.43 ± 0.01) | 0.09 ± 0.01 (0.08 ± 0.01) |
Petilla: The active and passive properties of L1 neurons grouped according to the Petilla classification scheme (c-AC, NGFC, FS, and BS). Number of neurons are given in brackets. Clustering: unsupervised clustering (hierarchical clustering using Ward's method) was performed and a cut off criteria was determined to reveal 4 clusters. Number of neurons in each cluster are given in brackets. The primary composition of neurons in each cluster when identified by the Petilla classification scheme are indicated: 81% of neurons in cluster 1 were NGFCs, 96% of neurons in cluster 2 FS neurons, 94% of neurons in cluster 3 c-AC neurons, and 89% of neurons in cluster 4 BS neurons. The average values for the active and passive properties for the primary composition of neurons in each cluster. All values are given as mean ± standard error of the mean.

Electrophysiological characterization of L1 neurons. (A) Representative AP firing patterns of L1 neurons (upper traces), generated in response to the injection of positive current steps (lower traces). L1 neurons were classified according to their firing properties with respect to the Petilla classification scheme: c-AC, NGFCs, FS, and BS neurons. (B–D) The cumulative probability distributions of initial firing frequency, accommodation index, and half-width of the first AP1 are shown for each type of L1 neurons (color code: c-AC: black, n = 137; NGFC: green, n = 72; FS: orange, n = 23; BS: blue, n = 34).
When we examined the passive properties of neurons grouped by their AP properties, we found that the input resistance (Rin) and RMP of each group possessed surprisingly narrow boundaries. c-AC cells had an average input resistance of 166 ± 5 MΩ (n = 137), NGFC cells 107 ± 3 MΩ (n = 72), FS cells 87 ± 7 MΩ (n = 23), and BS cells 201 ± 15 MΩ (n = 34) (ANOVA on ranks, Kruskal–Wallis followed by a Dunn's test: P < 0.05 for all pairwise multiple comparisons, except for c-AC vs. BS and NGFC vs. FS). Seventy-nine percent of NGFC cells and 83% of FS cells had Rin of less than 130 MΩ, whereas 80% of either c-AC cells and 82% BS cells had Rin of more than 130 MΩ (Fig. 2A,B). The RMP of 81% of NGFC cells and 78% of FS cells had an RMP depolarized to −70 mV, whereas the RMP of 23% of c-AC cells were hyperpolarized to −70 mV (Fig. 2C). BS cells showed a heterogeneous RMP distribution (Fig. 2C). These results are summarized in Table 1 (RMP: ANOVA on ranks, Kruskal–Wallis followed by a Dunn's test: P < 0.05 for all pairwise multiple comparisons, except for c-AC vs. BS and NGFC vs. FS).

Passive properties of c-AC, NGFCs, FS, and BS neurons. (A) Example traces of voltage responses (upper traces) evoked the injection of a 100 pA negative current step (lower trace) in each cell type. The RMP is indicated. (B) Distribution of apparent input resistance (Rin, measured at steady state, indicated as ΔV by the arrow in A) and (C) RMP of each cell type. Data are presented as cumulative probability plots (color code: c-AC: black, n = 137; NGFC: green, n = 72; FS: orange, n = 23; BS: blue, n = 34).
Anatomical Features of L1 Neurons
Neurons were regularly filled with biocytin and subsequently processed for morphological reconstruction (n = 71). In 40 of 71 neurons, the axonal arborization was sufficiently stained for further analysis. Typical reconstructions for each neuronal type are shown in Figure 3. The neurites of 13 of 21 c-AC neurons were confined to L1 (62%, Fig. 3A), 5 had 1 or 2 axonal branches descending to deeper layers (Supplementary Fig. 2B), and a further 2 neurons had more than 2 descending axonal branches. One c-AC neuron showed a divergent morphology. All successfully reconstructed NGFCs (n = 12) had characteristic anatomical features, typified by a large and dense axonal field and short dendrites (Fig. 4B) (Tamas et al. 2003; Olah et al. 2007; Szabadics et al. 2007). The axonal arborization of 5 NGFCs was confined to layer 1, while 7 neurons had multiple descending axons (Fig. 3B). Located close to the L1-L2/3 boarder 5 FS cells were reconstructed; 1 had multipolar dendrites and an axonal arbor in deep L1 and superficial part of L2 (Fig. 3C) and 4 were found to be chandelier cells, which had the well-described feature of several hundreds of vertically oriented axonal segments targeting the axon initial segments of pyramidal neurons (Supplementary Fig. 3) (Somogyi 1977; Howard et al. 2005). The active properties of FS neurons that were anatomically identified as chandelier cells had unique electrophysiological properties, exhibiting a left-shifted relationship between injected current and AP firing frequency (Supplementary Fig. 4A) and transient stuttering AP firing upon threshold current injection (Supplementary Fig. 4C1,4C2), which was not be observed in other FS cells (Supplementary Fig. 4B) (Woodruff et al. 2009). Two BS neurons were also found to have multipolar dendrites and axonal branches ending in layer 1 and layer 2/3 (Fig. 3D).

Anatomical reconstructions of L1 neurons. (A1) Typical example of a c-AC neuron with neurites confined to L1, dendrites are shown in black, and axon in red. (B1) Representative example of an NGFC with anatomical characteristics of a large and dense axonal field and short dendrites. (C1) Reconstruction of an FS neuron located at the boarder of L1 and L2. (D1) Example of a BS neuron with multipolar dendritic and axonal arbor in L1 and L2/3. Laminar borders are indicated with dashed lines. (A2–D2) AP firing pattern of the reconstructed L1 neurons (upper traces) in response to current injection (lower traces).

Cluster analysis of the electrophysiological properties of L1 neurons. (A) Tree diagram obtained using Ward's method. Five parameters were used to cluster neurons: initial firing rate, accommodation index, half-width of the first AP, Rin, and RMP. The dashed line indicates the chosen cut off. Clusters are numbered from 1 to 4. (B) Pie chart illustrating the composition of each cluster with reference to the Petilla classification scheme. The first cluster comprises 81 neurons, 66 of which had been classified as NGFCs and 15 as c-ACs. The second cluster is composed of 24 neurons, 23 FS and 1 c-AC neurons, the third cluster of 117 c-AC, 6 NGFCs, 2 BS neurons, and the fourth cluster of 32 BS and 4 c-AC neurons (color code: c-AC: black, NGFC: green, FS: orange, BS: blue).
Unsupervised Clustering
We used unsupervised cluster analysis (using Ward's method) to independently classify L1 neurons. Five electrophysiological parameters were used in this analysis: initial firing rate, accommodation index, half-width of the first AP, Rin, and RMP. Figure 4A shows the tree diagram of clusters based on these parameters; when a threshold of 4 clusters was selected, we found that neuronal grouping was similar to the Petilla classification. The first cluster was composed of 66 NGFCs and 15 c-AC neurons. The second cluster contained 96% FS cells. The majority of L1 neurons were grouped in the third cluster: 117 c-AC neurons, 6 NGFCs, and 2 BS neurons. The fourth cluster was composed of 89% BS cells (Fig. 4). Thus, unsupervised cluster analysis confirmed the presence of 4 groups: NGFCs, FS, c-AC, and BS neurons. For all subsequent analysis, we excluded all neurons that did not unequivocally belong to the same group using both classification schemes. We therefore excluded 8% of NGFC, 15% of c-AC, and 6% of BS cells from our initial groups based on the Petilla scheme (Fig. 4B). The electrophysiological properties of neurons grouped according to the Petilla classification scheme and following refinement based on cluster analysis are summarized in Table 1.
Synaptic Interaction between L2/3 Pyramidal Neurons and L1 Interneurons
Next, we investigated how efficiently L2/3 pyramidal neurons excited L1 interneurons (vertical distance up to ∼250 μm, horizontal distance less than 100 μm between L2/3 Pyr–L1 pairs). The properties of unitary excitatory postsynaptic potentials (uEPSPs) were found to be dependent on the identity of the postsynaptic neuron. We analyzed the latency, the rise time, the amplitude, and the half-width of the uEPSPs (Fig. 5). L2/3-FS uEPSPs had a latency of 1.4 ± 0.2 ms, a rise time of 1.1 ± 0.1 ms, amplitude of 1.1 ± 0.3 mV, a half-width of 5.7 ± 0.7 ms, and a decay time constant of 7.4 ± 1.1 ms (n = 7, Supplementary Fig. 5). L2/3-c-AC uEPSPs had a latency of 1.8 ± 0.1 ms (P = 0.06), a slower mean rise time of 2.2 ± 0.2 ms (P < 0.001), a similar average mean amplitude of 1.1 ± 0.2 mV, a longer half-width of 17.2 ± 1.2 ms (P < 0.001), and a greater decay time constant of 25.0 ± 4.2 (P < 0.001, n = 21, Supplementary Fig. 5). L2/3-BS connections had a latency similar to the L2/3-c-AC and L2/3-FS connections (1.4 ± 0.3 ms, P > 0.05), a rise time of 2.4 ± 0.1 ms, an amplitude of 1.1 ± 0.4 mV, a half-width of 14.5 ± 2.0 ms, and a decay time constant of 21.5 ± 3.5 ms (n = 5, Supplementary Fig. 5). The latter properties were not significantly different from the properties of L2/3-c-AC connections but distinct from L2/3-FS connections (rise time, P < 0.001; half-width, P < 0.05; decay time constant, P < 0.05). The connectivity probability between identified L2/3 pyramidal neurons and L1 interneurons was calculated as the ratio of the number of synaptic connected neurons to the number of connections tested (Fig. 5). Surprisingly, none of the NGFCs tested received excitatory input from L2/3 pyramidal neurons (68 connections tested), whereas L2/3 pyramidal neurons had a high probability of connectivity with FS cells (0.35, 7/20 connections tested, significantly different from L2/3 Pyr-NGFC ratio, P < 0.001, Fisher's Exact test, and from L2/3 Pyr-c-AC and L2/3 Pyr-BS ratio, P < 0.05, chi-squared test). c-AC and BS neurons were found to receive excitatory input from L2/3 pyramidal neurons with a probability of 0.16 (21/134) and 0.12 (5/42), respectively (both significantly different from L2/3-NGFC ratio, Fisher's Exact test, P < 0.05). These data indicate that L2/3 pyramidal neurons target L1 interneurons in a class-specific manner.

Properties of unitary EPSPs recorded in different types of L1 interneurons evoked by single AP firing in L2/3 pyramidal neurons. (A) Overlain traces of uEPSPs recorded from a c-AC (black), FS (orange), and BS(green) neuron. Average uEPSPs are shown as thick black traces. Pre, presynaptic L2/3 pyramidal neuron, post, postsynaptic L1 neuron. Inset shows probability of connectivity for each type of L1 neuron. (B–E) Distributions of the peak amplitude (B), half-width (C), 10–90% rise time (D), and latency (E) of uEPSPs recorded from 21 c-AC (n = 792 uEPSPs, black), 7 FS (n = 182 uEPSPs, orange), and 5 BS neurons (n = 111 uEPSPs, blue). Data are presented as cumulative probability plots.
We next investigated if the use-dependent properties (use-dependent depression or facilitation) of the excitatory synaptic input to classes of L1 interneurons were target cell dependent. To test this, we evoked 2 APs in L2/3 pyramidal neurons at an interval of 50 ms. Notably, use-dependent plasticity of L2/3 Pyr-FS and L2/3 Pyr-c-AC uEPSPs was distinct, L2/3 Pyr-FS uEPSPs exhibited paired-pulse facilitation (1.39 ± 0.09; n = 7), while L2/3 Pyr-c-AC connections displayed, on average, paired-pulse depression (0.92 ± 0.07; n = 21; P < 0.01; Fig. 6). L2/3 Pyr-BS connections, however, showed both facilitation and depression upon paired-pulse stimulation (1.1 ± 0.2, n = 5, data not shown).

Target-specific short-term plasticity of L2/3 to L1 unitary EPSPs. (A1, B1) AP firing pattern of a c-AC (left) and an FS neuron (right, upper traces) in response to current injection (lower traces). (A2, B2) Paired-pulse stimulation in synaptically coupled L2/3 pyramidal–L1 neuron pairs (same neurons as shown in A1, B1). Top trace (pre, presynaptic L2/3 pyramidal neuron) and 10 consecutive uEPSPs (post, postsynaptic L1 neuron) that were on average (lower trace) depressing in a c-AC and facilitating in an FS neuron. (C) Bar plot of the paired-pulse ratio (EPSP2/EPSP1) for L2/3 Pyr-c-AC (n = 21) and L2/3 Pyr-FS (n = 7) connections (Student's t-test, **P < 0.01).
Inhibitory Synaptic Connectivity between L1 Interneurons and L2/3 Pyramidal Neurons
The properties of unitary inhibitory postsynaptic potentials (uIPSPs) recorded from L2/3 pyramidal neurons were found to be dependent on the identity of the presynaptic L1 neuron. In contrast to excitatory connectivity, NGFCs were the most prominent source of inhibition to L2/3 pyramidal neurons, with a connectivity probability of 0.44 (12 connected out 27 connections tested). NGFC-L2/3 Pyr uIPSPs were of long latency, slowly rising, and large amplitude (latency 4.1 ± 0.9 ms, rise time 53.2 ± 10.8 ms, amplitude 0.58 ± 0.1 mV, and half-width 100 ± 19 ms, n = 12). In support of recent observations (Tamas et al. 2003), we observed that a single AP in NGFCs evoked uIPSPs mediated by the activation of both GABAA and GABAB receptors (GABAA receptor antagonist SR 95531, 2 μM and GABAB receptor antagonist CGP 52432, 10 μM, n = 10, Fig. 7E,F). In contrast, the other class of regular-spiking L1 neuron, c-AC cells, generated uIPSPs with a significantly shorter latency (1.8 ± 0.3 ms, P < 0.05), smaller amplitude (0.28 ± 0.04 mV, P < 0.001), and faster kinetics (rise time 6.5 ± 0.5 ms, P < 0.001; half-width 29 ± 8 ms, P < 0.01; n = 4; Supplementary Fig. 5; decay time constant 39.6 ± 6.9 ms [data not shown]). Notably, the probability of connectivity between cAC and L2/3 pyramidal neurons was low (0.06, 4 out of 67 connections tested). Furthermore, c-AC neurons inhibited L2/3 pyramidal neurons through the activation of GABAA receptors alone as uIPSPs were completely blocked by SR 95531 (2 μM, Fig. 7B,7D, c-AC: n = 3).

Properties of unitary IPSPs recorded from L2/3 pyramidal neurons evoked by single AP firing in L1 neurons. (A) Typical example of the firing pattern of a c-AC neuron (upper trace) in response to current injection (lower traces). (B) Average uIPSPs recorded from an L2/3 pyramidal neuron (average of 10 uIPSPs, Post) evoked by AP firing in a c-AC L1 neuron (Pre, same neuron as in A). uIPSPs were blocked by the GABAA receptor antagonist SR 95531 (2 μM). (C) Typical example of the firing pattern of an FS neuron. (D) Average uIPSPs in this FS-L2/3 pair were blocked by SR 95531 (2 μM, Pre, same neuron as in C). (E) Firing pattern of an NGFC. (F) Average uIPSPs in an NGFC-L2/3 connection (Pre, same neuron as in E). Note the longer time base. The slow component of the uIPSPs was blocked by the GABAB receptor antagonist CGP 52432 (10 μM, shaded area in lower panel) and the fast component by SR 95531 (2 μM). B, D, F, Insets show probability of connectivity for each type of L1 neuron.
FS neurons are known to act as fast signaling devices in other neuronal circuits (Jonas et al. 2004). Consistent with this role in the superficial cortical layers, FS cells evoked uIPSPs in L2/3 pyramidal neurons with a shorter latency than NGFCs and a faster rise time than other classes of L1 interneurons (uIPSP latency 0.9 ± 0.1 ms, FS vs. NGFC P < 0.01, FS vs. c-AC P = 0.08; rise time 4 ± 0.4 ms, FS vs. NGFC P < 0.001, FS vs. c-AC P <0.05; amplitude 0.27 ± 0.04 mV, FS vs. NGFC P < 0.001, FS vs. c-AC P = 0.77; half-width 27 ± 1 ms, FS vs. NGFC P < 0.001, FS vs. c-AC P = 0.85; n = 3; Supplementary Fig. 5; and decay time constant 31.7 ± 9.5 ms [data not shown]). Moreover, a high connection probability was found between FS and L2/3 pyramidal neurons (0.33, 3/9 connections tested). One chandelier cell was found to be synaptically connected with a layer 2/3 pyramidal neuron and following reconstruction, putative synaptic contacts were found decorating the axon initial segment of the postsynaptic neuron (Supplementary Fig. 7). In line with recent studies (Szabadics et al. 2006; Woodruff et al. 2006, 2009); but see (Glickfeld et al. 2009), we observed that the chandelier cell evoked depolarizing synaptic responses at a membrane potential of −60 mV (Supplementary Fig. 7), in contrast to hyperpolarizing uIPSPs evoked in L2/3 pyramidal neurons by all other classes of L1 interneurons at this membrane potential (see Fig. 7). In common with hyperpolarizing uIPSPs evoked by FS neurons, depolarizing Chandlier-L2/3 Pyr uIPSPs were blocked by the GABAA receptor antagonist SR 95531 (2 μM, FS: n = 2; Supplementary Fig. 6). Interestingly, none of the BS neurons we recorded were synaptically connected with L2/3 pyramidal neurons (20 connections tested); however, one BS-BS connection confirmed the GABAergic nature of this connection (Supplementary Fig. 8).
Discussion
In this study, we have shown that neocortical L1 neurons can be divided into classes according to the Petilla classification scheme (Ascoli et al. 2008). Hierarchical unsupervised cluster analysis confirmed the division of neurons into 4 groups. L2/3 pyramidal neurons were found to excite L1 neurons in a cell class–specific manner. Surprisingly, we found no excitatory connection between L2/3 pyramidal neurons and NGFCs. In contrast, NGFCs provided strong inhibition of L2/3 pyramidal neurons.
How Many Classes of Interneuron Exist in L1 of the Neocortex?
Previous studies have distinguished, at least, 2 classes of L1 neurons based on electrophysiological and anatomical criteria in brain slices prepared from early postnatal (P7–24) rats (Hestrin and Armstrong 1996; Zhou and Hablitz 1996a; Chu et al. 2003). Local circuit neurons were found to have axonal and dendritic arbors restricted to L1, while deeper layer–projecting neurons possessed descending axons that innervated deeper neocortical layers (Hestrin and Armstrong 1996; Zhou and Hablitz 1996a; Chu et al. 2003). Chu et al. (2003) identified delayed AP firing neurons as NGFCs, that were restricted to L1, while regular firing neurons had axonal arborizations that invaded the deeper layer of the neocortex. In this study, we investigated the properties of L1 neurons maintained in brain slices prepared from P24–36 rats. Over this postnatal age range, the electrophysiological properties of pyramidal (Zhu 2000; Atkinson and Williams 2009) and interneurons (Doischer et al. 2008) as well as the use-dependent dynamics of unitary excitatory and inhibitory synaptic transmission reach adult levels (Reyes and Sakmann 1999; Angulo et al. 1999; Frick et al. 2007; Feldmeyer and Radnikow 2009). In the mature neocortex, we have identified 4 classes of L1 neurons using unsupervised hierarchical cluster analysis of active and passive electrophysiological properties, a classification scheme supported by morphological analysis.
The role of inhibitory interneurons in neuronal computation is largely determined by the subcellular sites at which inhibitory synaptic contacts are made on pyramidal neurons (Somogyi et al. 1998). For example, classes of cortical interneurons target the axon initial segment, perisomatic regions, or the distal dendrites of pyramidal neurons (Markram et al. 2004; Pouille and Scanziani 2004; Ascoli et al. 2008; Klausberger and Somogyi 2008). Here, we find that the axons of 60% of c-AC neurons and more than 40% of NGFCs were restricted to L1, suggesting that the majority of L1 neurons inhibit pyramidal neurons at the level of the apical dendritic arbor. There are only few reports on FS neurons described in L1 (Foehring et al. 2002). L1 is relatively void of parvalbumin-positive cells, a calcium-binding protein expressed in FS neurons of other neocortical layers (Xu et al. 2010). We have, however, also identified a sparse population (23 of 266) of neurons in L1 with electrophysiological properties of FS neurons that possessed anatomical characteristics of axon or perisomatic-targeting interneurons. We cannot rule out that a small fraction of these neurons have been assigned to the wrong cortical layer because of the proximity of their somata with the L1-L2/3 border. However, based on the finding that only 5 of 271 neurons recorded from layer 1 had properties reminiscent of layer 2/3 pyramidal neurons, we estimate this error to be less than 2%. In support of the view that FS L1 neurons target proximal sites of L2/3 pyramidal neurons, we observed that uIPSPs recorded between FS L1 neurons and L2/3 pyramidal neurons had significantly faster rise times than those generated by dendritic targeting c-AC neurons, presumably because of the influence of dendritic filtering (Miles et al. 1996). As the rise time of uIPSPs also depends on the subunit composition of postsynaptic GABAA receptors, it is, however, possible that GABAA receptor composition may influence these findings (Farrant and Nusser 2005).
The minority of L1 neurons possessed axonal projections to deeper layers of the neocortex and were electrophysiologically classified as BS neurons, NGFCs, and chandelier cells. As we failed to observe any inhibitory synaptic input from BS neurons to L2/3 pyramidal neurons, we speculate that such a pathway is sparse and that BS neurons might belong to a group of interneurons which have previously been described to preferentially target other interneurons (Supplementary Fig. 8; (Gulyas et al. 1996; Caputi et al. 2009).
Neuronal Microcircuits in Superficial Layers of the Neocortex
It is widely believed that a canonical microcircuit exists in various cortical areas. In sensory cortices, the flow of information in excitatory neurons of the microcircuit is well described: Thalamocortical projection neurons mainly target spiny stellate cells in L4, which, in turn, innervate L2/3 pyramidal neurons (Feldmeyer et al. 2002; Shepherd and Svoboda 2005). L2/3 pyramidal neurons target neighboring L2/3 (Feldmeyer et al. 2006) and L5 pyramidal neurons (Thomson and Bannister 1998; Schubert et al. 2001; Williams and Atkinson 2007). However, L2/3 and L5 pyramidal neurons also show extensive axonal arborization within L1 (Feldmeyer et al. 2006; Brown and Hestrin 2009). We have found that various interneuron types of L1 receive monosynaptic excitatory input from L2/3 pyramidal neurons. We find the highest connectivity probability in L2/3-FS cell pairs, followed by L2/3-c-AC and L2/3-BS connections. As FS and c-AC interneurons were found to inhibit L2/3 pyramidal neurons, we propose that these cell types may function to provide feedback inhibition.
The use-dependent properties of excitatory synapses have been suggested to determine their role in network function, by shaping the flow of information in neuronal circuits during repetitive AP firing (Abbott and Regehr 2004). Target-specific use-dependent modifications of uEPSPs have been reported in several excitatory synaptic connections in the neocortex (Tsodyks and Markram 1997; Reyes et al. 1998; Markram et al. 1998; Gupta et al. 2000; Rozov et al. 2001; Koester and Johnston 2005; West et al. 2006; Williams and Atkinson 2007). Indeed, the use-dependent properties of uEPSPs evoked from a single presynaptic pyramidal neuron has been found to be dependent upon the identity of the postsynaptic partner (Reyes et al. 1998; Markram et al. 1998; Koester and Johnston 2005). In line with this, we found that the use-dependent properties of excitatory synaptic transmission was dependent upon the postsynaptic target in L1, with L2/3-FS uEPSPs showing paired-pulse facilitation, and L2/3-c-AC depression. Paired recordings from adult rat brain slices have shown that L5 and L6 pyramidal neurons generate facilitating uEPSPs in FS neurons (Angulo et al. 1999; West et al. 2006), whereas synaptic depression has been reported for L2/3-FS synapses in younger animals (P14) (Reyes et al. 1998) and for L4 regular-spiking neurons onto FS neurons (P14–21) (Beierlein et al. 2003). A developmental study, however, has demonstrated that in motor cortices from Wistar rats excitation onto FS undergo marked changes during maturation switching from use-dependent depression at P14–P20 to facilitation at P27–P36 (Angulo et al. 1999).
L1 neurons have been shown to powerfully control the excitability of L2/3 pyramidal neurons in vivo (Shlosberg et al. 2006). However, the functional role of different classes of L1 inhibitory neurons is unknown. Here, we have revealed the existence of 2 inhibitory feedback circuits between L2/3 pyramidal and L1 neurons. We suggest that the contrasting use-dependent properties observed in the excitatory arm of these circuits suggest that c-AC neurons, which target the apical dendritic arbor of L2/3 pyramidal neurons and perisomatically targeting FS neurons will be differentially engaged in a frequency-dependent manner. Both cell types provide feedback inhibition when L2/3 neurons fire APs at low frequencies, however, use-dependent depression at L2/3-c-AC synapses will ensure that dendritic inhibition is dampened during higher frequencies of AP firing. This frequency-dependent routing of feed-back inhibition is in contrast to that found in area CA1 of the hippocampus, where pyramidal to interneuron synapses show use-dependent depression for interneurons that innervate perisomatic sites but facilitation for interneuron that innervate apical dendritic sites of pyramidal neurons (Pouille and Scanziani 2004). In contrast to other classes of L1 neuron, we failed to find excitatory synaptic drive to NGFCs from L2/3 pyramidal neurons. We, however, found that NGFCs generate combined GABAA- and GABAB-mediated inhibitory responses to powerfully inhibit L2/3 pyramidal neurons. Notably, in line with the long latency of NGFCs evoked IPSPs, NGFCs have been shown to release GABA into the extracellular space providing widespread inhibitory input to neighboring neurons (Olah et al. 2009). Because of the large and dense axonal arborization of NGFCs in layer 1, these neurons are ideally positioned in the columnar circuit to control synaptic integration in distal apical dendrites of L2/3 and L5 pyramidal neurons (Waters et al. 2003; Williams 2004; Larkum et al. 2007; Larkum et al. 2009). Indeed, the activation of GABAB receptors has been shown to powerfully inhibit apical dendritic spike generation in layer 5 pyramidal neurons (Perez-Garci et al. 2006). Similar to hippocampal feed-forward inhibition (Elfant et al. 2008), we suggest that NGFCs may not receive their excitatory input from the columnar circuitry but from long-range intracortical excitatory pathways, conveyed in L1, that synapses with both the apical dendrites of pyramidal neurons (Petreanu et al. 2009) and L1 interneurons (Anderson and Martin 2006). We speculate therefore that NGFCs function to provide powerful feed-forward inhibition of long-range intracortical excitatory input.
Funding
Intramural funding from the Medical Research Council to S.R.W.; EU Marie-Curie Fellowship to C.W.
We would like to thank the present and previous members of Stephen Williams' laboratory for helpful discussions. We are grateful to Rick Livesey from the Gurdon Institute, Cambridge, UK and the members of the Gurdon Imaging Facility for initial help with confocal images. Conflict of Interest: None declared.