# Beyond the frontiers of neuronal types

Demian Battaglia, Anastassios Karagiannis, Thierry Gallopin, Harold W. Gutch, Bruno Cauli

► To cite this version:

Demian Battaglia, Anastassios Karagiannis, Thierry Gallopin, Harold W. Gutch, Bruno Cauli. Beyond the frontiers of neuronal types. Frontiers in Neural Circuits, 2013, 7, pp.13. 10.3389/fncir.2013.00013. hal-01537215

HAL Id: hal-01537215

https://hal.sorbonne-universite.fr/hal-01537215

Submitted on 12 Jun 2017

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
<figure><img src="image_2.png" /></figure>

# Beyond the frontiers of neuronal types

Demian Battaglia<fnref n="1" />,<fnref n="2" />*, Anastassios Karagiannis<fnref n="3" />,<fnref n="4" />, Thierry Gallopin<fnref n="5" />, Harold W. Gutch<fnref n="1" />,<fnref n="2" />,<fnref n="6" /> and Bruno Cauli<fnref n="3" />*

**Edited by:**

Yoshiyuki Kubota, National Institute for Physiological Sciences, Japan

**Reviewed by:**

Ariel Agmon, West Virginia University, USA
Rosa Cossart, Institut National de la Recherche Médicale, France

**\*Correspondence:**

Demian Battaglia, Max Planck Institute for Dynamics and Self-Organization, and Bernstein Center for Computational Neuroscience, Am Fößberg 17, D-37077 Göttingen, Germany.
e-mail: demian@nld ds.mpg.de

Bruno Cauli, CNRS UMR7102, Laboratoire de Neurobiologie des processus adaptatifs, Université Pierre et Marie Curie, 9 Quai St Bernard, F-75005 Paris, France.
e-mail: bruno.cauli@snv.jussieu.fr

Cortical neurons and, particularly, inhibitory interneurons display a large diversity of morphological, synaptic, electrophysiological, and molecular properties, as well as diverse embryonic origins. Various authors have proposed alternative classification schemes that rely on the concomitant observation of several multimodal features. However, a broad variability is generally observed even among cells that are grouped into a same class. Furthermore, the attribution of specific neurons to a single defined class is often difficult, because individual properties vary in a highly graded fashion, suggestive of continua of features between types. Going beyond the description of representative traits of distinct classes, we focus here on the analysis of atypical cells. We introduce a novel paradigm for neuronal type classification, assuming explicitly the existence of a structured continuum of diversity. Our approach, grounded on the theory of fuzzy sets, identifies a small optimal number of model archetypes. At the same time, it quantifies the degree of similarity between these archetypes and each considered neuron. This allows highlighting archetypal cells, which bear a clear similarity to a single model archetype, and edge cells, which manifest a convergence of traits from multiple archetypes.

**Keywords:** neuronal diversity, interneuron diversity, barrel cortex, petilla terminology, atypical cells, fuzzy sets, unsupervised clustering

## INTRODUCTION

Cortical physiology and function depend on a delicate interplay between excitatory glutamatergic neurons and diverse inhibitory GABAergic interneurons in a temporally, spatially and cell-type-specific manner (Gupta et al., 2000; Somogyi and Klausberger, 2005; Sohal et al., 2009; Cauli and Hamel, 2010; Kätzel et al., 2011; Gentet, 2012). A pertinent understanding of information processing requires thus a thorough description of this neocortical neuronal diversity.

Historically, inhibitory interneurons have been subdivided according to a large repertoire of morphological, electrophysiological, and molecular properties (Ascoli et al., 2008). More recently, classification schemes based on developmental criteria largely confirmed and further supported (Butt et al., 2005; Miyoshi et al., 2007; Batista-Brito and Fishell, 2009; Vucurovic et al., 2010) previous classifications (McCormick et al., 1985; Kawaguchi and Kubota, 1993, 1996; Cauli et al., 1997). This emphasizes the concept that neuronal types can and must be defined by the convergence of common features and not solely by a limited set of prescribed properties (Tyner, 1975).

Unsupervised clustering methods have so far established themselves as the state-of-the-art approach to identify neuronal types based on the simultaneous consideration of many features (Tamás et al., 1997; Cauli et al., 2000; Toledo-Rodriguez et al.,

2004; Gallopin et al., 2006; Halabisky et al., 2006; Dumitriu et al., 2007; Helmstaedter et al., 2009; Karagiannis et al., 2009; McGarry et al., 2010; Suzuki and Bekkers, 2010; Perrenoud et al., 2012a,b). Although useful, these classification schemes are hampered by at least two shortcomings. First, a distinctive trait can be shared by multiple cell types. For instance, parvalbumin (PV) expression is common to both some basket cells and chandelier cells (Kawaguchi and Kubota, 1993; Toledo-Rodriguez et al., 2004; Woodruff et al., 2010). Second, a feature discriminative for certain neuronal types can be irrelevant and highly variable for other types. Examples are spike width and input resistance, respectively, very thin and small for Fast Spiking (FS)-PV cells, but highly heterogeneous for both Somatostatin-(SOM) and Vasoactive Intestinal Polypeptide (VIP)-expressing cells (Kawaguchi and Kubota, 1996; Cauli et al., 2000; Ma et al., 2006).

Because of these shortcomings, the boundaries between different types are blurred. Yet, despite this, all clustering methods used so far for neuronal classification inherently exclude the existence of a graded separation between neuronal classes (Tyner, 1975; Parra et al., 1998). As an alternative, we propose here a novel approach, based on the theory of Fuzzy Sets (Zadeh, 1965, 2008), in which the absence of a strict dividing line between neuronal types is actually assumed as an explicit categorizing hypothesis.

<footnote n="1">¹ Department of Nonlinear Dynamics, Max Planck Institute for Dynamics and Self-Organization (MPIDS), Göttingen, Germany</footnote>

<footnote n="2">² Bernstein Center for Computational Neuroscience, Göttingen, Germany</footnote>

<footnote n="3">³ CNRS UMR7102, Laboratoire de Neurobiologie des processus adaptatifs, Université Pierre et Marie Curie, Paris, France</footnote>

<footnote n="4">⁴ Neuroscience, Physiology and Pharmacology, Division of Biosciences, University College London, London, UK</footnote>

<footnote n="5">⁵ CNRS UMR 7637, Laboratoire de Neurobiologie et Diversité Cellulaire, ESPCI ParisTech, Paris, France</footnote>

<footnote n="6">⁶ Department of Mathematics, Technische Universität München, Garching, Germany</footnote>
By applying this approach to a well-characterized multi-parametric dataset (Karagiannis et al., 2009), we confirm the identification of four main types of inhibitory neurons (canonic archetypes). Beyond this, our method provides a way to quantify how individual cells are similar to these archetypes, opening the way to a systematic study of deviations from typical tendencies. In particular, we unmask atypical cells displaying properties simultaneously reminiscent of multiple archetypes. These edge cells might contribute to robustness and flexibility of performance in neural circuits by enriching, with their polyvalent phenotypes, the reservoir of available functional resources.

# MATERIALS AND METHODS

## THE DATASET

We performed here an alternative analysis of a dataset presented in Karagiannis et al. (2009). This previous study focused on typical properties of interneuron classes, while we emphasize here deviations from these typical aspects. The sample consisted of 200 neurons from slices from juvenile male Wistar rat somatosensory cortex (postnatal days 19 ± 2, SD, Charles River, L'Arbresle, France). Overall 43 features were measured for each neuron, including digitized value of laminar location of the soma (with non-integer values indicating cells at the border between different layers, e.g., 1.5 meaning between layers I and II), 32 electrophysiological features (5 sub-threshold, 4 just-above threshold, 5 firing, and 18 action potential properties) and expression of 10 well-established molecular markers of neuronal diversity (Ascoli et al., 2008). Morphology, reconstructed only for a limited subset of cells, was not considered in the current analysis. A complete list of the used features is reported in Tables 1–4. The reader is invited to refer to (Karagiannis et al., 2009) for all details on feature definitions and on the used experimental procedures.

## FUZZY SETS AND PARTITIONS

In Fuzzy Set Theory (Zadeh, 1965; Jang et al., 1997; Xu et al., 2008), a data-point *i* can belong to a fuzzy class *α* with different degrees of membership, quantified by a membership value *m*<sub>*iα*</sub>. A membership value *m*<sub>*iα*</sub> = 1.0 indicates that the data-point *i* displays all the defining attributes of the class *α* and therefore fully belongs to it. A membership value *m*<sub>*iα*</sub> = 0.0 indicates that the data-point *i* does not display any of the defining attributes of the class *α* and therefore it does not belong to it at all. In addition to these “black and white” cases, and extending ordinary set theory, intermediate membership values 0 < *m*<sub>*iα*</sub> < 1 are admitted (Figure 1A), corresponding to a continuum of “gray” cases in which the data-point *i* displays only some of the defining attributes of class *α*. Fuzzy Set Theory thus provides, in the words of its initiator,

<table><caption>Table 1 | Localization and subthreshold features of different archetypes.</caption><thead><tr><th></th><th>Glutamatergic<br>(n = 49)</th><th>FS-PV<br>(n = 33)</th><th>Adapting SOM<br>(n = 21)</th><th>Adapting VIP<br>(n = 31)</th><th>Adapting NPY<br>(n = 33)</th><th>UFO<br>(n = 33)</th></tr></thead><tbody><tr><td>(1) Digitized laminar location (Layer)</td><td>3.4 ± 0.8</td><td>3.1 ± 0.8</td><td>2.8 ± 0.7</td><td>2.7 ± 0.7</td><td>2.3 ± 0.6</td><td>2.2 ± 0.7</td></tr><tr><td></td><td colspan="6">UFO, Adapt. NPY &lt; Adapt. VIP, Adapt. SOM &lt;&lt; FS-PV, Glutamatergic</td></tr><tr><td>(2) Resting potential (V<sub>m</sub>, mV)</td><td>-74.3 ± 4.7</td><td>-72.7 ± 4.0</td><td>-65.0 ± 3.9</td><td>-71.0 ± 5.1</td><td>-72.2 ± 4.5</td><td>-70.8 ± 5.5</td></tr><tr><td></td><td colspan="6">Glutamatergic, FS-PV, Adapt. NPY, "UFO," Adapt. VIP &lt;&lt; Adapt. SOM<br>Glutamatergic &lt;&lt; UFO; Glutamatergic &lt; Adapt. VIP</td></tr><tr><td>(3) Input resistance (R<sub>m</sub>, MΩ)</td><td>379 ± 138</td><td>204 ± 77</td><td>272 ± 76</td><td>514 ± 166</td><td>314 ± 116</td><td>552 ± 266</td></tr><tr><td></td><td colspan="6">FS-PV &lt;&lt; Adapt. SOM &lt;&lt; Glutamatergic &lt;&lt; Adapt. VIP, UFO<br>FS-PV &lt;&lt;&lt; Adapt. NPY &lt;&lt;&lt; Adapt. VIP, "UFO"</td></tr><tr><td>(4) Time constant (τ<sub>m</sub>, ms)</td><td>36.1 ± 10.2</td><td>15.6 ± 6.0</td><td>24.2 ± 8.7</td><td>24.6 ± 10.0</td><td>22.8 ± 8.3</td><td>29.0 ± 12.0</td></tr><tr><td></td><td colspan="6">FS-PV &lt;&lt;&lt; Adapt. NPY, Adapt. VIP &lt; UFO &lt;&lt; Glutamatergic<br>FS-PV &lt;&lt;&lt; Adapt. SOM &lt; UFO</td></tr><tr><td>(5) Membrane capacitance (C<sub>m</sub>, pF)</td><td>104.0 ± 32.3</td><td>79.1 ± 21.7</td><td>91.3 ± 23.4</td><td>48.5 ± 14.3</td><td>76.8 ± 22.7</td><td>56.0 ± 20.3</td></tr><tr><td></td><td colspan="6">Adapt. VIP, UFO &lt;&lt;&lt; Adapt. NPY, FS-PV &lt;&lt; Glutamatergic<br>Adapt. VIP, "UFO" &lt;&lt;&lt; Adapt. SOM</td></tr><tr><td>(6) Sag index (Sag, %)</td><td>19.7 ± 9.2</td><td>9.4 ± 5.3</td><td>27.6 ± 12.0</td><td>8.2 ± 3.7</td><td>9.0 ± 4.8</td><td>6.4 ± 3.2</td></tr><tr><td></td><td colspan="6">UFO &lt; Adapt. VIP, FS-PV &lt;&lt;&lt; Glutamatergic &lt; Adapt. SOM<br>UFO &lt;&lt; Adapt. NPY &lt;&lt;&lt; Glutamatergic &lt; Adapt. SOM</td></tr></tbody></table>

Values are weighted means ± SD; n, number of cells; < significantly smaller with $P \le 0.05$; << significantly smaller with $P \le 0.01$; &lt;&lt; significantly smaller with $P \le 0.001$.
<table><caption>Table 2 | Above threshold features of different archetypes.</caption><thead><tr><th></th><th>Glutamatergic<br>(n = 49)</th><th>FS-PV<br>(n = 33)</th><th>Adapting SOM<br>(n = 21)</th><th>Adapting VIP<br>(n = 31)</th><th>Adapting NPY<br>(n = 33)</th><th>UFO<br>(n = 33)</th></tr></thead><tbody><tr><td>(7) Rheobase (I<sub>rheo</sub>, pA)</td><td>35.8 ± 26.2</td><td><strong>103.1 ± 46.4</strong></td><td><strong>-4.3 ± 30.2</strong></td><td>16.3 ± 12.7</td><td>53.0 ± 28.5</td><td>32.7 ± 22.0</td></tr><tr><td></td><td colspan="6"><strong>Adapt. SOM</strong> &lt;&lt;&lt; Adapt. VIP &lt;&lt; "UFO," Glutamatergic &lt;&lt; Adapt. NPY &lt;&lt;&lt; <strong>FS-PV</strong></td></tr><tr><td>(8) First spike latency<br>(τ<sub>1st</sub>, ms)</td><td>127.1 ± 46.5</td><td>370.7 ± 245.4</td><td>110.8 ± 88.3</td><td>125.2 ± 113.7</td><td>218.5 ± 209.9</td><td><strong>106.6 ± 112.2</strong></td></tr><tr><td></td><td colspan="6"><strong>"UFO"</strong> &lt; Adapt. SOM &lt; Adapt. VIP, Glutamatergic, Adapt. NPY<br>Glutamatergic, Adapt. SOM &lt;&lt; FS-PV</td></tr><tr><td>(9) Just-above threshold<br>adaptation (Ad<sub>thr</sub>, Hz/s)</td><td><strong>-46.8 ± 74.9</strong></td><td><strong>3.3 ± 20.4</strong></td><td>-15.6 ± 25.9</td><td>-6.9 ± 18.4</td><td><strong>-0.4 ± 4.1</strong></td><td><strong>-32.6 ± 38.6</strong></td></tr><tr><td></td><td colspan="6"><strong>Glutamatergic, "UFO"</strong> &lt; Adapt. VIP &lt;&lt; Adapt. NPY, FS-PV<br>Adapt. SOM &lt;&lt; Adapt. NPY, FS-PV</td></tr><tr><td>(10) Minimal steady state<br>frequency (F<sub>min</sub>, Hz)</td><td>30.9 ± 53.8</td><td>14.9 ± 13.3</td><td>13.7 ± 12.3</td><td>13.1 ± 10.1</td><td><strong>6.1 ± 3.0</strong></td><td><strong>31.7 ± 27.1</strong></td></tr><tr><td></td><td colspan="6"><strong>Adapt. NPY</strong> &lt;&lt;&lt; Adapt. VIP, Adapt. SOM, FS-PV, <strong>"UFO"</strong>; Adapt. NPY &lt;&lt; Glutamatergic<br>Adapt. VIP, Adapt. SOM &lt; <strong>"UFO"</strong></td></tr><tr><td>(11) Amplitude<br>accommodation (A<sub>hump</sub>, mV)</td><td><strong>24.3 ± 8.8</strong></td><td><strong>1.1 ± 1.2</strong></td><td>3.8 ± 4.3</td><td>5.2 ± 4.8</td><td>10.6 ± 7.4</td><td>7.1 ± 5.6</td></tr><tr><td></td><td colspan="6"><strong>FS-PV</strong> &lt;&lt;&lt; Adapt. VIP, "UFO" &lt; Adapt. NPY &lt;&lt;&lt; <strong>Glutamatergic</strong><br>FS-PV &lt;&lt; Adapt. SOM &lt; "UFO"; Adapt. SOM &lt;&lt;&lt; Adapt. NPY</td></tr><tr><td>(12) Amplitude of early<br>adaptation (A<sub>sat</sub>, Hz)</td><td><strong>153.1 ± 59.4</strong></td><td><strong>52.7 ± 24.8</strong></td><td>93.9 ± 24.3</td><td>110.7 ± 41.3</td><td>127.4 ± 37.0</td><td>109.8 ± 37.1</td></tr><tr><td></td><td colspan="6"><strong>FS-PV</strong> &lt;&lt;&lt; Adapt. SOM, "UFO," Adapt. VIP &lt;&lt; <strong>Glutamatergic</strong><br>Adapt. SOM &lt;&lt; Adapt. NPY &lt; <strong>Glutamatergic</strong></td></tr><tr><td>(13) Time constant of early<br>adaptation (τ<sub>sat</sub>, ms)</td><td>25.7 ± 11.6</td><td><strong>18.4 ± 15.9</strong></td><td><strong>37.8 ± 7.8</strong></td><td>23.9 ± 8.8</td><td>25.8 ± 6.1</td><td>27.1 ± 8.8</td></tr><tr><td></td><td colspan="6"><strong>FS-PV</strong> &lt; Adapt. VIP, Adapt. NPY, Glutamatergic, "UFO" &lt;&lt;&lt; <strong>Adapt. SOM</strong></td></tr><tr><td>(14) Late adaptation<br>(Ad<sub>sat</sub>, Hz/s)</td><td><strong>-9.6 ± 7.3</strong></td><td>-27.6 ± 13.1</td><td>-20.6 ± 9.2</td><td><strong>-35.6 ± 11.3</strong></td><td>-21.3 ± 10.0</td><td>26.1 ± 14.9</td></tr><tr><td></td><td colspan="6"><strong>Adapt. VIP</strong> &lt; FS-PV, "UFO," Adapt. NPY, Adapt. SOM &lt;&lt;&lt; <strong>Glutamatergic</strong><br>Adapt. VIP &lt;&lt; "UFO"; Adapt. VIP &lt;&lt;&lt; Adapt. NPY, Adapt. SOM</td></tr><tr><td>(15) Maximal steady state<br>frequency (F<sub>max</sub>, Hz)</td><td><strong>29.9 ± 8.3</strong></td><td><strong>140.4 ± 30.4</strong></td><td>69.1 ± 19.9</td><td>87.0 ± 29.2</td><td>61.1 ± 13.9</td><td>64.3 ± 22.0</td></tr><tr><td></td><td colspan="6"><strong>Glutamatergic</strong> &lt;&lt;&lt; Adapt. NPY, "UFO" &lt;&lt; Adapt. VIP &lt;&lt;&lt; <strong>FS-PV</strong><br>Glutamatergic &lt;&lt;&lt; Adapt. SOM &lt;&lt;&lt; <strong>FS-PV</strong></td></tr></tbody></table>

Values are weighted means ± SD; n, number of cells; < significantly smaller with $P \le 0.05$; &lt; significantly smaller with $P \le 0.01$; &lt;&lt; significantly smaller with $P \le 0.001$.

"a precise language to describe imprecise similarities" (Zadeh, 2008).

A dataset can be partitioned into *c* fuzzy classes. In a fuzzy partition each data-point can belong to multiple classes with different degrees of membership. Formally, a fuzzy partition of a dataset with elements *i* = 1... *N* into *c* fuzzy clusters is uniquely determined by a set of *N* membership vectors *m*<sub>*i*</sub> with *c* components each. The *α*-th entry *m*<sub>*i*</sub>*α* of the vector *m*<sub>*i*</sub> gives the membership of the data-point *i* to the *α*-th class in the fuzzy partition. The memberships vectors are normalized in such a way that $\sum_{\alpha=1}^{c} m_{i\alpha} = 1$.

## UNSUPERVISED FUZZY CLUSTERING

Each cell *i* (*i* = 1... *N*) is represented as a vector of features *f*<sub>*i*</sub> with entries *f*<sub>*i*</sub> where the index *ℓ* = 1... 43 runs over the considered features (see rows of Tables 1–4). Each entry corresponds to the centered and reduced value of the corresponding feature (Cauli et al., 2000; Karagiannis et al., 2009). An unsupervised clustering algorithm was used to determine an optimized fuzzy partition of the dataset. We resorted to the fuzzy *c*-means algorithm (Dunn, 1973; Jang et al., 1997; Xu et al., 2008), a fuzzy analog of the crisp *k*-means algorithm previously used for unsupervised classification of interneurons (Karagiannis et al., 2009).
<table><caption>Table 3 | Action potentials features of different archetypes.</caption><thead><tr><th></th><th>Glutamatergic<br>(n = 49)</th><th>FS-PV<br>(n = 33)</th><th>Adapting SOM<br>(n = 21)</th><th>Adapting VIP<br>(n = 31)</th><th>Adapting NPY<br>(n = 33)</th><th>"UFO"<br>(n = 33)</th></tr></thead><tbody><tr><td>(16) First spike amplitude<br>(A<sub>1</sub>, mV)</td><td>93.0 ± 7.2</td><td><strong>81.9 ± 7.5</strong></td><td>94.9 ± 7.6</td><td><strong>98.0 ± 6.5</strong></td><td>93.4 ± 5.7</td><td>91.5 ± 13.0</td></tr><tr><td></td><td colspan="6"><strong>FS-PV</strong> &lt;&lt;&lt; Glutamatergic, Adapt. NPY, Adapt. SOM &lt;&lt; <strong>Adapt. VIP</strong><br><strong>FS-PV</strong> &lt;&lt; "UFO" &lt; <strong>Adapt. VIP</strong></td></tr><tr><td>(17) Second spike amplitude<br>(A<sub>2</sub>, mV)</td><td>85.2 ± 11.7</td><td><strong>79.1 ± 8.0</strong></td><td>90.7 ± 7.7</td><td>92.7 ± 7.5</td><td>90.0 ± 5.7</td><td><strong>79.9 ± 9.3</strong></td></tr><tr><td></td><td colspan="6"><strong>FS-PV</strong>, "UFO" &lt;&lt;&lt; Adapt. NPY, Adapt. SOM, Adapt. VIP<br><strong>FS-PV</strong>, "UFO" &lt; Glutamatergic &lt;&lt; Adapt. VIP</td></tr><tr><td>(18) First spike duration<br>(D<sub>1</sub>, ms)</td><td><strong>1.4 ± 0.2</strong></td><td><strong>0.6 ± 0.2</strong></td><td>0.9 ± 0.2</td><td>0.8 ± 0.2</td><td>1.0 ± 0.2</td><td>1.0 ± 0.3</td></tr><tr><td></td><td colspan="6"><strong>FS-PV</strong> &lt;&lt;&lt; Adapt. VIP &lt;&lt;&lt; Adapt. NPY, "UFO" &lt;&lt;&lt; <strong>Glutamatergic</strong><br>Adapt. VIP &lt; Adapt. SOM &lt;&lt;&lt; <strong>Glutamatergic</strong></td></tr><tr><td>(19) Second spike duration<br>(D<sub>2</sub>, ms)</td><td><strong>1.6 ± 0.3</strong></td><td><strong>0.6 ± 0.2</strong></td><td>1.0 ± 0.2</td><td>0.8 ± 0.2</td><td>1.1 ± 0.2</td><td>1.1 ± 0.3</td></tr><tr><td></td><td colspan="6"><strong>FS-PV</strong> &lt;&lt;&lt; Adapt. VIP &lt;&lt;&lt; Adapt. NPY, "UFO" &lt;&lt;&lt; <strong>Glutamatergic</strong><br>Adapt. VIP &lt;&lt; Adapt. SOM &lt;&lt;&lt; <strong>Glutamatergic</strong></td></tr><tr><td>(20) Amplitude Reduction<br>(Δ(Amp), %)</td><td>-8.4 ± 9.7</td><td><strong>-1.2 ± 6.5</strong></td><td>-4.8 ± 3.7</td><td>-5.3 ± 4.9</td><td>-3.4 ± 2.7</td><td>-11.8 ± 9.7</td></tr><tr><td></td><td colspan="6"><strong>FS-PV</strong> &lt;&lt;&lt; Adapt. VIP &lt; "UFO"; <strong>FS-PV</strong> &lt;&lt;&lt; Glutamatergic<br><strong>FS-PV</strong> &lt;&lt; Adapt. NPY &lt;&lt;&lt; "UFO"; Adapt. NPY &lt; Glutamatergic<br><strong>FS-PV</strong> &lt;&lt; Adapt. SOM &lt; "UFO"</td></tr><tr><td>(21) Duration Increase<br>(Δ(Dur), %)</td><td>10.7 ± 12.0</td><td><strong>1.3 ± 5.3</strong></td><td>7.4 ± 4.6</td><td>2.8 ± 4.3</td><td>9.0 ± 7.4</td><td><strong>12.9 ± 10.1</strong></td></tr><tr><td></td><td colspan="6"><strong>FS-PV</strong> &lt;&lt; Adapt. VIP &lt;&lt;&lt; Adapt. SOM, Adapt. NPY, "UFO"<br><strong>FS-PV</strong> &lt;&lt;&lt; Glutamatergic; Adapt. VIP &lt;&lt; Glutamatergic</td></tr><tr><td>(22) First spike, first<br>component AHP (AHP<sub>f</sub>, mV)</td><td>-6.7 ± 3.7</td><td><strong>-23.2 ± 3.1</strong></td><td>-11.4 ± 3.8</td><td>-14.1 ± 3.8</td><td>-13.7 ± 4.4</td><td>-13.1 ± 4.6</td></tr><tr><td></td><td colspan="6"><strong>FS-PV</strong> &lt;&lt;&lt; Adapt. VIP, Adapt. NPY, "UFO," Adapt. SOM &lt;&lt;&lt; <strong>Glutamatergic</strong></td></tr><tr><td>(23) First spike, second<br>component AHP (AHP<sub>s</sub>, mV)</td><td>-8.0 ± 6.9</td><td>-2.1 ± 6.3</td><td>-6.5 ± 4.2</td><td>-8.6 ± 3.2</td><td><strong>-15.3 ± 3.8</strong></td><td><strong>-1.9 ± 4.2</strong></td></tr><tr><td></td><td colspan="6"><strong>Adapt. NPY</strong> &lt;&lt;&lt; Adapt. VIP &lt; Adapt. SOM &lt;&lt; FS-PV &lt;&lt;&lt; "UFO"<br><strong>Adapt. NPY</strong> &lt;&lt;&lt; Glutamatergic &lt; FS-PV</td></tr><tr><td>(24) Second spike, first<br>component AHP (AHP<sub>f,2</sub>, mV)</td><td><strong>-8.8 ± 4.1</strong></td><td><strong>-23.3 ± 3.4</strong></td><td>-10.3 ± 6.6</td><td>-15.5 ± 4.1</td><td>-15.3 ± 4.0</td><td>-15.0 ± 5.0</td></tr><tr><td></td><td colspan="6"><strong>FS-PV</strong> &lt;&lt;&lt; Adapt. NPY &lt;&lt; Adapt. SOM &lt; <strong>Glutamatergic</strong><br><strong>FS-PV</strong> &lt;&lt;&lt; Adapt. VIP, "UFO" &lt; Adapt. SOM</td></tr><tr><td>(25) Second spike, second<br>component AHP (AHP<sub>s,2</sub>, mV)</td><td><strong>-14.7 ± 8.8</strong></td><td><strong>-2.1 ± 5.7</strong></td><td>-7.5 ± 3.9</td><td>-10.7 ± 2.6</td><td><strong>-13.2 ± 8.3</strong></td><td><strong>-2.8 ± 5.3</strong></td></tr><tr><td></td><td colspan="6"><strong>Glutamatergic</strong> &lt;&lt;&lt; Adapt. VIP &lt;&lt;&lt; "UFO," <strong>FS-PV</strong><br><strong>Adapt. NPY</strong> &lt;&lt; Adapt. VIP &lt;&lt; Adapt. SOM &lt;&lt; "UFO," <strong>FS-PV</strong></td></tr><tr><td>(26) First spike, first AHP<br>component latency (τ<sub>AHPf</sub>, ms)</td><td>14.2 ± 20.7</td><td><strong>2.9 ± 1.3</strong></td><td>3.2 ± 1.1</td><td><strong>2.7 ± 0.7</strong></td><td>4.7 ± 1.7</td><td>4.2 ± 2.2</td></tr><tr><td></td><td colspan="6"><strong>Adapt. VIP, FS-PV</strong> &lt;&lt;&lt; "UFO" &lt;&lt;&lt; <strong>Glutamatergic</strong><br><strong>Adapt. VIP, FS-PV</strong> &lt;&lt;&lt; Adapt. NPY &lt;&lt; <strong>Glutamatergic</strong><br><strong>FS-PV</strong> &lt; Adapt. SOM &lt;&lt;&lt; Adapt. NPY</td></tr><tr><td>(27) First spike, second AHP<br>component latency (τ<sub>AHPs</sub>, ms)</td><td>46.0 ± 42.3</td><td><strong>1.1 ± 3.1</strong></td><td>27.5 ± 16.8</td><td>34.0 ± 17.2</td><td>20.1 ± 6.8</td><td><strong>5.1 ± 12.1</strong></td></tr><tr><td></td><td colspan="6"><strong>FS-PV</strong>, "UFO" &lt;&lt;&lt; Adapt. NPY, Adapt. SOM, Adapt. VIP, Glutamatergic<br>Adapt. NPY &lt;&lt;&lt; Adapt. VIP</td></tr></tbody></table>

(Continued)
<table><caption>**Table 3 | Continued**</caption><thead><tr><th></th><th>Glutamatergic<br>(n = 49)</th><th>FS-PV<br>(n = 33)</th><th>Adapting SOM<br>(n = 21)</th><th>Adapting VIP<br>(n = 31)</th><th>Adapting NPY<br>(n = 33)</th><th>"UFO"<br>(n = 33)</th></tr></thead><tbody><tr><td>(28) Second spike, first AHP<br>component latency (τ<sub>AHPf.2</sub>, ms)</td><td>18.9 ± 23.5</td><td>2.9 ± 1.3</td><td>3.4 ± 1.0</td><td>2.8 ± 0.7</td><td>7.0 ± 4.4</td><td>5.7 ± 4.2</td></tr><tr><td></td><td colspan="6">Adapt. VIP, FS-PV &lt;&lt;&lt; "UFO" &lt;&lt;&lt; Glutamatergic<br>Adapt. VIP, FS-PV, Adapt. SOM &lt;&lt;&lt; Adapt. NPY &lt;&lt; Glutamatergic<br>FS-PV &lt;&lt; Adapt. SOM &lt;&lt; "UFO"</td></tr><tr><td>(29) Second spike, second AHP<br>component latency (τ<sub>AHPs.2</sub>, ms)</td><td>59.3 ± 38.6</td><td>1.1 ± 3.1</td><td>30.1 ± 16.9</td><td>31.4 ± 14.0</td><td>16.4 ± 9.6</td><td>10.3 ± 18.6</td></tr><tr><td></td><td colspan="6">FS-PV &lt;&lt;&lt; Adapt. NPY &lt;&lt;&lt; Adapt. VIP &lt;&lt; Glutamatergic<br>FS-PV &lt;&lt;&lt; Adapt. SOM &lt;&lt;&lt; Glutamatergic; Adapt. NPY &lt; Adapt. SOM<br>"UFO" &lt;&lt;&lt; Adapt. VIP, Glutamatergic; "UFO" &lt;&lt; Adapt. SOM; "UFO" &lt; Adapt. NPY</td></tr><tr><td>(30) First spike ADP (ADP, mV)</td><td>1.3 ± 1.8</td><td>0.5 ± 1.5</td><td>4.7 ± 2.5</td><td>7.9 ± 3.4</td><td>0.6 ± 0.7</td><td>0.3 ± 0.8</td></tr><tr><td></td><td colspan="6">"UFO," FS-PV, Adapt. NPY, Glutamatergic &lt;&lt;&lt; Adapt. SOM, Adapt. VIP<br>"UFO," FS-PV &lt;&lt; Adapt. NPY; "UFO" &lt;&lt; Glutamatergic</td></tr><tr><td>(31) Second spike ADP<br>(ADP<sub>2</sub>, mV)</td><td>0.3 ± 0.6</td><td>0.5 ± 1.6</td><td>3.7 ± 2.1</td><td>7.0 ± 2.8</td><td>0.2 ± 0.4</td><td>0.7 ± 2.0</td></tr><tr><td></td><td colspan="6">Adapt. NPY, Glutamatergic, FS-PV, "UFO" &lt;&lt;&lt; Adapt. SOM &lt;&lt;&lt; Adapt. VIP</td></tr><tr><td>(32) First spike ADP latency<br>(τ<sub>ADP</sub>, ms)</td><td>5.6 ± 6.1</td><td>0.7 ± 1.9</td><td>10.4 ± 5.7</td><td>9.8 ± 2.7</td><td>5.1 ± 3.5</td><td>1.5 ± 3.9</td></tr><tr><td></td><td colspan="6">FS-PV, "UFO" &lt;&lt; Adapt. NPY &lt;&lt;&lt; Adapt. VIP, Adapt. SOM;<br>FS-PV, "UFO" &lt;&lt; Glutamatergic<br>Glutamatergic &lt;&lt; Adapt. VIP; Glutamatergic &lt; Adapt. SOM</td></tr><tr><td>(33) Second spike ADP latency<br>(τ<sub>ADP2</sub>, ms)</td><td>3.6 ± 4.2</td><td>0.6 ± 1.7</td><td>8.5 ± 4.2</td><td>9.8 ± 2.1</td><td>3.0 ± 3.6</td><td>1.8 ± 4.1</td></tr><tr><td></td><td colspan="6">FS-PV &lt;&lt; Glutamatergic &lt;&lt;&lt; Adapt. SOM, Adapt. VIP<br>"UFO," Adapt. NPY &lt;&lt;&lt; Adapt. SOM, Adapt. VIP; FS-PV &lt;&lt; Adapt. NPY</td></tr></tbody></table>

Values are weighted means ± SD; n, number of cells; < significantly smaller with $P \le 0.05$; &lt;&lt; significantly smaller with $P \le 0.01$; &lt;&lt;&lt; significantly smaller with $P \le 0.001$.

This algorithm has two adjustable parameters. The first parameter is the maximally allowed number of classes, c. Analogously to k-means, the algorithm attempts to partition the dataset into c different fuzzy clusters. However, unlike in k-means, some fuzzy clusters can actually “compenetrate” and finally coalesce, leading to a fuzzy partition with a smaller effective number of clusters. We assumed here a quite large potential number of archetypes, c = 20. The second parameter is the fuzziness parameter $\mu > 1$, allowing to control the level of fuzziness of the obtained partition. In the limit of $\mu$ approaching unit from higher values, the fuzzy c-means algorithm is equivalent to k-means for $k = c$. Larger values of this parameter lead to increasingly fuzzier partitions, with individual data-points sharing their membership between a smaller number of effective fuzzy classes. The fuzziness parameter was varied in the range $1 < \mu < 2$ for the analysis of Figure 2 and set to an optimal value of $\mu = 1.349$ for the rest of the analyses in Figures 3–6 (see below for criteria guiding the choice of $\mu$).

The fuzzy c-means algorithm builds a fuzzy partition of the neuronal dataset through an iterative optimization process. At a given iteration t, c cluster centroids are given by vectors $\mathbf{u}_{\alpha}^{(t)}$

($\alpha = 1...c$) with components $u_{\alpha\ell}^{(t)}$. Associated membership vectors $\mathbf{m}_{i}^{(t)}$ are computed as:

$$ \frac{1}{m_{i\alpha}^{(t)}} = \sum_{\lambda=1}^{c} \left( \frac{d_{i\lambda}^{(t)}}{d_{i\lambda}^{(t)}} \right)^{\frac{2}{\mu-1}} \quad (1) $$

where $d_{i\lambda}^{(t)}$ is the Euclidean distance between the data-point $f_i$ and the centroid $\mathbf{u}_{\lambda}^{(t)}$.

These membership vectors are used in turn to compute a new set of cluster centroids $\mathbf{u}_{\alpha}^{(t+1)}$ with coordinates:

$$ u_{\alpha\ell}^{(t+1)} = \frac{\sum_{i=1}^{N} (m_{i\alpha}^{(t)})^{\mu} f_{i\ell}}{\sum_{i=1}^{N} (m_{i\alpha}^{(t)})^{\mu}} \quad (2) $$

This procedure is designed to minimize a specific cost function (Dunn, 1973; Jang et al., 1997), namely the sum of the squared
<table><caption>Table 4 | Occurrence of molecular markers in different archetypes.</caption><thead><tr><th></th><th>Glutamatergic<br>(n = 49)</th><th>FS-PV<br>(n = 33)</th><th>Adapting SOM<br>(n = 21)</th><th>Adapting VIP<br>(n = 31)</th><th>Adapting NPY<br>(n = 33)</th><th>"UFO"<br>(n = 33)</th></tr></thead><tbody><tr><td>(34) VGlut1</td><td><strong>99%</strong></td><td>26%</td><td>22%</td><td>7%</td><td>34%</td><td>21%</td></tr><tr><td></td><td colspan="6">Adapt. VIP, "UFO," Adapt. SOM, FS-PV, Adapt. NPY &lt;&lt;&lt; <strong>Glutamatergic</strong></td></tr><tr><td>(35) GAD</td><td>3%</td><td><strong>100%</strong></td><td><strong>100%</strong></td><td><strong>100%</strong></td><td><strong>100%</strong></td><td><strong>100%</strong></td></tr><tr><td></td><td colspan="6">Glutamatergic &lt;&lt;&lt; <strong>FS-PV</strong>, Adapt. SOM, Adapt. VIP, "UFO," Adapt. NPY</td></tr><tr><td>(36) NOS-1</td><td>0%</td><td>8%</td><td>3%</td><td>2%</td><td>23%</td><td>5%</td></tr><tr><td></td><td colspan="6">Glutamatergic &lt; Adapt. NPY</td></tr><tr><td></td><td>42%</td><td><strong>58%</strong></td><td><strong>77%</strong></td><td>8%</td><td>4%</td><td>6%</td></tr><tr><td>(37) CB</td><td colspan="6">Adapt. NPY, "UFO," Adapt. VIP &lt;&lt;&lt; <strong>FS-PV</strong>, Adapt. SOM<br>Adapt. NPY, "UFO" &lt;&lt;&lt; Glutamatergic &lt; <strong>Adapt. SOM</strong><br>Adapt. VIP &lt;&lt; Glutamatergic</td></tr><tr><td></td><td>32%</td><td><strong>100%</strong></td><td>27%</td><td>12%</td><td>26%</td><td>16%</td></tr><tr><td>(38) PV</td><td colspan="6">Adapt. VIP, "UFO," Adapt. NPY, Adapt. SOM, Glutamatergic &lt;&lt;&lt; <strong>FS-PV</strong></td></tr><tr><td></td><td>2%</td><td>2%</td><td>19%</td><td>22%</td><td>11%</td><td>24%</td></tr><tr><td>(39) CR</td><td colspan="6">Glutamatergic, FS-PV &lt; Adapt. VIP, "UFO"</td></tr><tr><td></td><td>2%</td><td>29%</td><td>63%</td><td>14%</td><td><strong>83%</strong></td><td>24%</td></tr><tr><td>(40) NPY</td><td colspan="6">Glutamatergic, Adapt. VIP, "UFO," FS-PV &lt;&lt;&lt; <strong>Adapt. NPY</strong><br>Glutamatergic &lt;&lt; "UFO," FS-PV, Adapt. SOM<br>Adapt. VIP &lt;&lt; Adapt. SOM</td></tr><tr><td></td><td>3%</td><td>2%</td><td>5%</td><td><strong>92%</strong></td><td>4%</td><td><strong>55%</strong></td></tr><tr><td>(41) VIP</td><td colspan="6">FS-PV, Glutamatergic, Adapt. NPY, Adapt. SOM &lt;&lt;&lt; <strong>"UFO"</strong> &lt; <strong>Adapt. VIP</strong></td></tr><tr><td></td><td>2%</td><td>6%</td><td><strong>88%</strong></td><td>0%</td><td>0%</td><td>9%</td></tr><tr><td>(42) SOM</td><td colspan="6">Adapt. VIP, Adapt. NPY, Glutamatergic, FS-PV, "UFO" &lt;&lt;&lt; <strong>Adapt. SOM</strong></td></tr><tr><td>(43) CCK</td><td>7%</td><td>0%</td><td>2%</td><td>10%</td><td>2%</td><td>12%</td></tr></tbody></table>

Values are weighted probabilities of expression; n, number of cells; > significantly larger with $P \le 0.05$; >> significantly larger with $P \le 0.01$; >>> significantly larger with $P \le 0.001$.

distances of the data-points from the different centroids, weighted by the relative fuzzy memberships:

$$J^{(t)} = \sum_{i=1}^{N} \sum_{\lambda=1}^{c} \left(m_{i\lambda}^{(t)}\right)^{\mu} \cdot \left(d_{i\lambda}^{(t)}\right)^{2} \quad (3)$$

In practice, as a first step, we randomly initialized a collection of c cluster centroids $u_{\alpha}^{(0)}$ in the feature space, by selecting c arbitrary data-points $f_i$. Initial membership vectors $m_i^{(0)}$ were then computed using Equation (1). Equations (1) and (2) were then iterated until the positions of the c centroids converged to a fixed point (with a prescribed tolerance) or until a fixed maximal number of iterations was reached. The final set of centroids was then inspected to identify potential coalescences and drop redundant centroids. Whenever,

the Euclidean distances between different centroids was smaller than a tolerance threshold (conventionally set to $\epsilon = 0.001$), the associated fuzzy classes were merged, and the membership vectors of data-points correspondingly shrunk to a length $c^* < c$, by adding up memberships of the merged classes. Thus, given a dataset and a maximum number of allowed clusters c, the effective number $c^*$ of clusters in the final fuzzy partition depended on $\mu$. The larger $\mu$ was, the more coalescences occurred (Figure 2E).

Note that fuzzy c-means, as classic k-means, can converge to different partitions for different initial conditions. We therefore used a seeding pre-initialization strategy (see Basu et al., 2002), by selecting 7 out of 20 initial centroids to match data-points close to the centroids of the 7 k-means clusters analyzed in Karagiannis et al. (2009), to obtain fuzzy partitions maximally correlated with our previous crisp analyses.
<figure><img src="image_3.png" /></figure>

## TYPICALITY QUANTIFICATION

We denoted as archetypal classes of diversity (or, briefly, archetypes) the c\* fuzzy classes in an optimized fuzzy partition. Given a set of archetypes and membership vectors for a prescribed fuzziness parameter μ, the main type of a cell i was the archetype α with the largest membership component $m_{iα} = m_i^{(1st)}$. Denoting then as $m_i^{(2nd)}$ the second largest component of $m_i$ (associated to a secondary type), we defined the typicality coefficient of i as

$$d(i) = m_i^{(1st)} - m_i^{(2nd)} \quad (4)$$

The typicality coefficient d is bounded in the interval $0 \le d \le 1$.

The sample average $\bar{d}$ of typicality coefficients and their dispersion, quantified robustly by the halved mean absolute deviation $\Delta = d - d/2$ (Sachs, 1984), were quantified over the entire sample. Based on the observed distribution of typicality coefficients (Figure 6B), cells were then considered as archetypal if $d > d_{\text{Archetypal}} = \bar{d} + \Delta$ and as atypical if $d < d_{\text{Atypical}} = \bar{d} - \Delta$.

Finally, we denoted as edge cells atypical neurons with a very small typicality coefficient such that $d \le \Delta$, reflecting very similar memberships to the main and to another secondary type.

## CHARACTERIZATION OF ARCHETYPAL PROPERTIES

Properties of different archetypes (cf. Tables 1–4) were computed as averages over the features of cells with a given main type, weighted by memberships toward this main type. Standard deviations were also evaluated using weighted expressions (Taylor, 1997). Defining $I_α$ as the subset of cells having archetype α as main type, we computed, for each feature $ℓ$

and each archetype α:

$$\bar{f}_{\ell}(\alpha) = \frac{\sum_{i \in I_{\alpha}} m_{i\alpha} f_{i\ell}}{\sum_{i \in I_{\alpha}} m_{i\alpha}} \quad (5)$$

$$\sigma_{\ell}(\alpha) = \sqrt{\frac{\left(\sum_{i \in I_{\alpha}} m_{i\alpha}\right) \cdot \sum_{i \in I_{\alpha}} m_{i\alpha} \left(f_{i\ell} - \bar{f}_{\ell}(\alpha)\right)^2}{\left(\sum_{i \in I_{\alpha}} m_{i\alpha}\right)^2 - \sum_{i \in I_{\alpha}} m_{i\alpha}^2}} \quad (6)$$

These weighted averages and standard deviations reflected dominantly feature values of archetypal cells and damped the contribution of atypical cells. The typical range of a feature $ℓ$ for an archetype α was defined accordingly as the interval $\bar{f}_{\ell}(\alpha) - \sigma_{\ell}(\alpha) \le f \le \bar{f}_{\ell}(\alpha) + \sigma_{\ell}(\alpha)$ and such convention was used to determine the coloring of table entries in the panel (C) of Figures 3–5.

Properties were then compared between every pair of archetypes α and β, by looking for differences between the distributions of features over the subsamples $I_α$ and $I_β$. Significance of pairwise comparisons was tested using the (two-tailed) Mann–Whitney U non-parametric test for continuous-valued features (i.e., lamination and electrophysiological properties) and Fisher’s exact test for Boolean features (i.e., expression of molecular markers). All significant comparisons are listed in Tables 1–4.

## RELEVANCE OF DIFFERENT PROPERTIES FOR CLASSIFICATION

To analyze the impact of different properties on the quality of the extracted fuzzy partition of our dataset we compared the cost of our fuzzy partition (i.e., the value of the cost function $J^{(t_{\text{end}})}$ as given by Equation (3) evaluated at the final iteration $t_{\text{end}}$ of the clustering algorithm) with the costs of partitions extracted from partially randomized datasets.

To perform the randomization, the dataset was represented as a matrix with 43 columns, corresponding to different features, and 200 rows, corresponding to different neurons. We built randomized datasets by permuting randomly and independently the order of entries within selected columns. Such a random scrambling does not alter means and standard deviations of the randomized properties but it does disrupt their correlations with other features. The randomized datasets were then clustered with fuzzy c-means, using the same parameters as for the original dataset. A loss in the quality of the extracted fuzzy partition was quantified by measuring the increase of the associated cost function $J^{(t_{\text{end}})}$. The larger the increase in residual cost after randomization of a given subset of properties, the more this property subset is considered to contribute to the quality of the reference classification.

To estimate the increase $\Delta J$ of the cost function due to randomization, the mean, and standard deviation of the achieved cost variation were taken over 1000 independent randomizations of each considered property subset.

This analysis is analogous to the one introduced in Karagiannis et al. (2009). In the present study, however, we quantify losses
A $\mu = 1.732$

<figure><img src="image_4.png" /></figure>

B $\mu = 1.627$

<figure><img src="image_5.png" /></figure>

<figure><img src="image_6.png" /></figure>

C $\mu = 1.389$

<figure><img src="image_7.png" /></figure>

D $\mu = 1.349$

<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /></figure>

**FIGURE 2** | **Emergence of archetypes.** Fuzzy partitions with decreasing μ values are visualized as membership matrices (**A–D**, upper panels). Rows correspond to different fuzzy clusters and columns to individual neurons. Membership values of single neurons to each class are color coded (right bars). Schematic representation of the interrelations between archetypes (lower panels). Colored discs depict different archetypes and their overlaps denote cells with shared memberships. The Glutamatergic archetype is the first to emerge (black, **A**), followed by FS-PV interneurons (Red, **B**), and Adapting VIP interneurons (blue, **C**). Note that the glutamatergic archetype splits into three sub-groups. Adapting SOM (green) and Adapting NPY (orange) archetypes and a residual UFO phenotype of highly atypical inhibitory interneurons are the last

groups to singularize (D). (E) Effective number of clusters generated by different fuzziness parameters. The number of archetypes included in the partition is indicated to the left of the corresponding range of μ. Partitions with more than eight archetypes can be only retrieved within very narrow ranges of low μ values (red). The range leading to classification with the largest number of robust archetypes is marked in green. (F) Metaphoric example illustrating the impact of fuzziness on the relevance of partitions and numbers of archetypes. Fuzzy partitions with too few archetypes (large μ, bottom) convey a too blurred image of the dataset. Conversely, fuzzy partitions with too many archetypes (small μ, top) are scarcely representative being strongly affected by outliers. Such issue is graphically represented as impulse noise on the image.

in classification quality by monitoring increases of fuzzy cost function, rather than decreases of average silhouette (Rousseuw, 1987). This is due to the fact that standard silhouette analysis does not take into account the membership profiles of different cells.

Different features were ranked in order of decreasing ΔJ. Fuzzy clustering was then performed taking into account only the top K highest ranked features, to extract reduced classifications based on a smaller number of relevant properties. The same c and the same seeding were always used for pre-initialization at every number of included properties K in order to guarantee the extraction of comparable fuzzy clusters. The fuzzy clusters in the reduced fuzzy partitions were mapped to the archetypes of the reference classification. More specifically, a fuzzy cluster

$\alpha_{\text{reduced}}$ in a reduced classification was matched to the phenotype $\alpha_{\text{ref}}$ in the original classification with the largest mutual overlap. For each phenotype, a classification matching fraction was then defined as the percentage of cells whose original main type $\alpha_{\text{ref}}$ coincided with its main type $\alpha_{\text{reduced}}$ in the considered reduced classification.

# RESULTS

## EMERGING ARCHETYPES

Following analogous approaches in numerical taxonomy and ecology (Bezdek, 1974; Equiha, 1990; Salski, 2007), we used an established unsupervised algorithm (Dunn, 1973) to extract from our sample—based on the same features considered in Karagiannis et al. (2009)— a small number of fuzzy clusters, or
<figure><img src="image_10.png" /></figure>

**FIGURE 3** | Walking outward from the FS-PV archetype. (A) Tetrahedral representation linking FS-PV (red), Adapting NPY (orange), and Adapting SOM archetypes (green). Individual neurons are depicted as dots in the 3D space of memberships. The X, Y, and Z axes correspond to memberships to the FS-PV, Adapting NPY, and Adapting SOM archetypes, respectively. Vertical stems denoting the projection of selected cells on the bidimensional base plane are introduced as depth clues. (B) Current-clamp recordings of illustrative transition neurons (cells a to o colored in panel A) in response to rheobase current and to a 100 pA hyperpolarizing current pulse (black traces, scale bars 50 mV, 400 ms). Insets: details of the repolarization phase of the

first spike (red traces, scale bars 5 mV, 20 ms). (C) Table summarizing 10 electrophysiological and 3 molecular discriminative properties of transition cells between the Adapting NPY, the Adapting SOM, and the FS-PV archetypes. Orange, green, and red backgrounds indicate distinctive values for the Adapting NPY, the Adapting SOM, or the FS-PV archetype, respectively. Gradient backgrounds indicate values falling in a range typical for multiple. Bold colored entries indicate extreme values for an archetypal trend. Thick contours highlight columns corresponding to archetypal cells. Atypical cells display a heterogeneous mixture of property values which are not compatible with a single archetype or which fall in transition ranges.

archetypes. Archetypes—like usual clusters (called “crisp” in the context of fuzzy set theory)— are characterized by canonic properties, manifesting common phenotypic trends. In addition, as detailed in the section “Materials and Methods,” our analysis provided lists of *memberships* (Zadeh, 1965; Jang et al., 1997; Xu et al., 2008), quantifying the similarity of each neuron with different archetypes. Cells were considered as *archetypal*, if they had a large membership toward a unique archetype and as *atypical*, if they had comparable memberships toward more archetypes (Figure 1B).

Our clustering algorithm depends on few parameters, most notably on the fuzziness parameter μ (see Materials and Methods), the tuning of which affected the number of identified

archetypes. For a sufficiently large fuzziness parameter, all cells were grouped indistinctly into a unique all-embracing fuzzy class. Decreasing the fuzziness, Glutamatergic neurons detached then first at μ = 1.81 and remained the only singularized archetype down to μ = 1.73 (Figures 2A,E). At this fuzziness level and down to μ = 1.63 a second class emerged (Figures 2B,E). Its GABAergic interneurons were PV-positive with electrophysiological properties reminiscent of FS interneurons and could be therefore identified as an archetype of FS-PV neurons (McCormick et al., 1985; Kawaguchi and Kubota, 1993).

At the fuzziness value of μ = 1.49 and down to μ = 1.37, a new GABAergic archetype differentiated (Figures 2C,E). It was mainly constituted by VIP-positive interneurons with an adapting
<figure><img src="image_11.png" /><figcaption>FIGURE 4 | Walking between Adapting GABAergic archetypes.</figcaption></figure>

(A) Tetrahedral representation linking Adapting VIP (blue), Adapting SOM (green), and Adapting NPY (orange) archetypes. Individual neurons are depicted as dots in the 3D space of memberships. The X, Y, and Z axes correspond to memberships to the Adapting VIP, Adapting NPY, and Adapting SOM archetypes, respectively. Vertical stems denoting the projection of selected cells on the bidimensional base plane are introduced as depth clues. (B) Current-clamp recordings of illustrative transition neurons (cells *g* to *r* colored in panel A) in response to rheobase current and to a 100 pA hyperpolarizing current pulse (black traces, scale bars 50 mV, 400 ms). Insets: details of the repolarization

phase of the first spike (red traces, scale bars 5 mV, 20 ms). (C) Table summarizing 8 electrophysiological and 3 molecular discriminative properties of transitions cells between the Adapting VIP, the Adapting SOM, and the Adapting NPY archetypes. Blue, green, and orange backgrounds indicate distinctive values for the Adapting VIP, the Adapting SOM, or the Adapting NPY prototype, respectively. Gradient backgrounds indicate values falling in a range typical for multiple. Bold colored entries indicate extreme values for an archetypal trend. Thick contours highlight columns corresponding to archetypal cells. Atypical cells display a heterogeneous mixture of property values which are not compatible with a single prototype or which fall in transition ranges.

firing pattern (Kawaguchi and Kubota, 1996; Cauli et al., 1997, 2000) and thus corresponded to the Adapting VIP prototype. Furthermore, within this same range of μ, the Glutamatergic type started differentiating into further subtypes.

In the range of 1.37 ≥ μ ≥ 1.21 two further archetypes of mostly SOM- or Neuropeptide Y (NPY)-positive interneurons progressively emerged (Figures 2D,E). The GABAergic

archetypes found in that range (Figure 2E) reproduced to a good extent the GABAergic classes obtained with the two independent crisp clustering algorithms we previously used in Karagiannis et al. (2009), confirming their reliability. In the obtained soft hierarchy of fuzzy archetypes, a fuzzy class encompassing all neurons with a yet uncharacterized main type was found at any μ and was denoted as the class of “Unidentified Firing Objects” (UFOs).
A

<figure><img src="image_12.png" /></figure>

B

<figure><img src="image_13.png" /></figure>

C

<figure><img src="image_14.png" /></figure>

**FIGURE 5** | **Visiting the UFOs.** (A) Tridimensional representation joining three heterogeneous UFOs. Adapting VIP (blue), Adapting NPY (orange), and UFOs (pink) archetypes. Individual neurons are depicted as dots in a 3D space of memberships. The X, Y, and Z axes correspond to memberships to the Adapting VIP, Adapting NPY, and UFOs archetypes, respectively. Vertical stems denoting the projection of selected cells on the bidimensional base plane are introduced as depth clues. (B) Current-clamp recordings of illustrative transition neurons (cells *s* to *u* colored in panel A) in response to

rhебос current and to a 100 pA hyperpolarizing current pulse (black traces, scale bars 50 mV, 400 ms). Insets: details of the repolarization phase of the first spike (red traces, scale bars 5 mV, 20 ms). (C) Table summarizing 17 electrophysiological and 5 molecular properties reminiscent of FS-PV (Red), Adapting VIP (blue), Adapting SOM (Green), or Adapting NPY (orange) archetypes in three different UFOs. Colored backgrounds indicate values falling in ranges typical for an archetype. Gradient backgrounds indicate values falling in a range typical for multiple archetypes.

Finally, classifications with a larger number of archetypes could not be robustly extracted from our dataset, since these partitions were valid only for very narrow ranges of μ (Figures 2E,F). In the rest of this study, we focus on the specific classification obtained for μ = 1.349, leading to four GABAergic archetypes, together with one (out of three merged) Glutamatergic archetype and a residual UFO class (Figures 2D,E).

Weighting features by their phenotype memberships limits the confounding contribution of atypical cells, so our fuzzy method allowed for a better quantitative characterization of the typical features of these archetypes (summarized in Tables 1–4).

### THREE GLUTAMATERGIC SUBCLASSES

The whole population of Glutamatergic neurons was distinctly characterized by relatively deeply located somata, hyperpolarized resting membrane potentials, and slow time constants due to large membrane capacitances (Table 1). Their action potentials were of long duration with slowly developing afterhyperpolarization (AHP) of small amplitude (Table 3). When strongly depolarized these neurons fired action potentials with a marked amplitude accommodation, a pronounced frequency adaptation, and were unable to sustain high firing rates (Table 2).

As previously described for glutamatergic neurons (Kubota et al., 1994; Andjelic et al., 2009), they expressed vGluT1 and only infrequently interneuronal markers, except for CB (Table 4). Although these neurons were the first to singularize, they subsequently split up into three subclasses at μ values required for the emergence of Adapting VIP, SOM, and NPY archetypes (Figures 2C–E). These classes of glutamatergic neurons only differed by the laminar location of their soma and by electrophysiological properties, as reported in Table 5.

Briefly, class A glutamatergic neurons whose somata were deeply located (Table 5) displayed a pronounced voltage sag at hyperpolarized potentials (Table 5). They distinctly differed from the other glutamatergic neurons by a first action potential exhibiting a very slowly developing and deep secondary AHP. Compared with other glutamatergic neurons the duration of their spikes was fairly constant (Table 5).

Class B glutamatergic neurons were the most superficial group of glutamatergic neurons and exhibited modest voltage sag when hyperpolarized (Table 5). The profile of their spike amplitude accommodation was of small amplitude for a glutamatergic neuron (Table 5). Their action potentials were of long duration and
<figure><img src="image_15.png" /></figure>

**FIGURE 6** | Archetype segregations. (A) Pairwise comparisons of phenotype segregations. Two-dimensional projections of memberships of neurons belonging to the Glutamatergic (black), FS-PV (red), Adapting SOM (green), Adapting VIP (blue), Adapting NPY (orange), and UFO (magenta) archetypes. Dashed lines represent identical memberships and gray zones the mean absolute deviation of typicality. Neurons falling within the gray zone correspond to "edge cells." Note the absence of "edge cells" between,

Glutamatergic, FS-PV, and Adapting VIP neurons. (B) Overall distribution of typicality coefficients. The right and left peaks correspond to archetypal and atypical cells, respectively. The grayed background denotes the range of typicalities associated to edge cells. In the stacked histogram, sections with different colors indicate cells with different main type. Archetypal and atypical cells are unequally distributed across archetypes. The bimodal distribution indicates that archetypes tend to separate, but only imperfectly.

their first component AHP was slowly developing and of relatively large amplitude (Table 5).

Similarly to class A glutamatergic neurons, *class C glutamatergic neurons* had deeply located somata and a very pronounced voltage sag (Table 5). In contrast to the other glutamatergic neurons, none of them were able to fire below a minimal frequency and displayed strong adaptation when depolarized just above threshold (Table 5), two firing features reminiscent of a bursting phenotype (Hodgkin, 1948). Their first, but not their second action potential displayed a monophasic repolarization. The duration and the amplitude of their second action potential were respectively longer and smaller than those of their first spike. Also consistent with a bursting phenotype they displayed a very pronounced amplitude accommodation when strongly depolarized. Altogether these features indicate that class C glutamatergic neurons correspond to intrinsically bursting neurons (McCormick et al., 1985; Connors and Gutnick, 1990).

From the different properties available for all neurons in the dataset, it was not possible to associate other distinctive properties to these classes of glutamatergic neurons. They included molecular (Table 5) and basic morphological somatic features

(data not shown). The absence of relevant markers for these neurons is attributable to the low occurrence of interneuron markers in glutamatergic neurons (Andjelic et al., 2009). In addition, the difficulty to correlate unambiguously the morphology of spiny granular and supragranular neurons with a firing pattern (McCormick et al., 1985; Connors and Gutnick, 1990; Staiger et al., 2004) further supports the idea that the three classes of glutamatergic neurons do not differ by their somatic features. Nevertheless, it is possible that these neurons exhibit distinctive patterns of dendritic and/or axonal arborizations, as it has been observed for bursting and non-bursting Layer V pyramidal cells (Chagnac-Amitai et al., 1990; Christophe et al., 2005). Since such properties are not available in our dataset, glutamatergic neurons will be considered as a single unified archetype in the rest of this study.

## FOUR GABAergic ARCHETYPES

Archetypal *FS-PV* interneurons had a very low input resistance, a very fast membrane time constant, and displayed the lowest electrical excitability (Tables 1 and 3, see also cell *a* in Figure 3 as an illustrative example). Their spikes were short with fast
<table><caption>Table 5 | Significant comparisons between glutamatergic subtypes.</caption><thead><tr><th></th><th>Glutamatergic class A<br>(n = 25)</th><th>Glutamatergic class B<br>(n = 12)</th><th>Glutamatergic class C<br>(n = 9)</th></tr></thead><tbody><tr><td>(1) Digitized laminar location (Layer)</td><td>3.7 ± 0.6</td><td><strong>2.5 ± 0.8</strong></td><td>4.0 ± 0.2</td></tr><tr><td></td><td colspan="3"><strong>Class B</strong> &lt;&lt; Classes A, C</td></tr><tr><td>(4) Time constant (τ<sub>m</sub>, ms)</td><td>33.4 ± 8.1</td><td><strong>42.6 ± 11.7</strong></td><td>34.2 ± 9.2</td></tr><tr><td></td><td colspan="3"><strong>Class A</strong> &lt; <strong>Class B</strong></td></tr><tr><td>(6) Sag index (Sag, %)</td><td>21.5 ± 7.5</td><td><strong>11.3 ± 7.7</strong></td><td>25.8 ± 7.6</td></tr><tr><td></td><td colspan="3"><strong>Class B</strong> &lt;&lt; Classes A, C</td></tr><tr><td>(9) Just-above threshold adaptation (Ad<sub>thr</sub>, Hz/s)</td><td>-18 ± 35</td><td>-4.7 ± 3.0</td><td><strong>-163 ± 76</strong></td></tr><tr><td></td><td colspan="3"><strong>Class C</strong> &lt;&lt;&lt; Classes A, B</td></tr><tr><td>(10) Minimal steady state frequency (F<sub>min</sub>, Hz)</td><td>14 ± 18</td><td><strong>7.9 ± 2.8</strong></td><td><strong>77 ± 15</strong></td></tr><tr><td></td><td colspan="3"><strong>Classes A, B</strong> &lt;&lt;&lt; <strong>Class C</strong></td></tr><tr><td>(11) Amplitude accommodation (A<sub>hump</sub>, mV)</td><td>25.2 ± 6.2</td><td><strong>15.6 ± 6.9</strong></td><td><strong>33.2 ± 8.0</strong></td></tr><tr><td></td><td colspan="3"><strong>Class B</strong> &lt;&lt; Class A &lt; <strong>Class C</strong></td></tr><tr><td>(13) Time constant of early adaptation (τ<sub>sat</sub>, ms)</td><td>23.6 ± 4.8</td><td>20.1 ± 3.5</td><td><strong>38 ± 18</strong></td></tr><tr><td></td><td colspan="3"><strong>Classes A, B</strong> &lt;&lt;&lt; <strong>Class C</strong></td></tr><tr><td>(14) Late adaptation (Ad<sub>sat</sub>, Hz/s)</td><td>-9.7 ± 6.6</td><td>-14.5 ± 8.4</td><td><strong>-3.7 ± 3.8</strong></td></tr><tr><td></td><td colspan="3"><strong>Classes A, B</strong> &lt;&lt; <strong>Class C</strong></td></tr><tr><td>(15) Maximal steady state frequency (F<sub>max</sub>, Hz)</td><td>30.5 ± 8.6</td><td>32.6 ± 6.8</td><td><strong>23.2 ± 6.4</strong></td></tr><tr><td></td><td colspan="3"><strong>Class C</strong> &lt;&lt; Class B</td></tr><tr><td>(17) Second spike amplitude (A<sub>2</sub>, mV)</td><td>90.4 ± 7.8</td><td>89.1 ± 9.3</td><td><strong>69.5 ± 5.2</strong></td></tr><tr><td></td><td colspan="3"><strong>Class C</strong> &lt;&lt;&lt; Classes A, B</td></tr><tr><td>(18) First spike duration (D<sub>1</sub>, ms)</td><td>1.40 ± 0.17</td><td><strong>1.60 ± 0.18</strong></td><td>1.41 ± 0.17</td></tr><tr><td></td><td colspan="3"><strong>Class A</strong> &lt;&lt;&lt; <strong>Class B</strong></td></tr><tr><td>(20) Amplitude Reduction (Δ(Amp), %)</td><td>-3.7 ± 2.0</td><td>-2.9 ± 1.9</td><td><strong>-24.8 ± 3.9</strong></td></tr><tr><td></td><td colspan="3"><strong>Classes A, B</strong> &lt;&lt;&lt; <strong>Class C</strong></td></tr><tr><td>(21) Duration Increase (Δ(Dur), %)</td><td>3.4 ± 2.8</td><td>8.3 ± 4.1</td><td><strong>31.3 ± 9.9</strong></td></tr><tr><td></td><td colspan="3"><strong>Class A</strong> &lt;&lt; <strong>Class B</strong> &lt;&lt;&lt; <strong>Class C</strong></td></tr><tr><td>(22) First spike, first component AHP (AHP<sub>f</sub>, mV)</td><td>-5.5 ± 2.6</td><td><strong>-10.1 ± 4.5</strong></td><td>-5.4 ± 2.5</td></tr><tr><td></td><td colspan="3"><strong>Class B</strong> &lt; Classes A, C</td></tr><tr><td>(23) First spike, second component AHP (AHP<sub>s</sub>, mV)</td><td><strong>-13.6 ± 2.8</strong></td><td>-3.5 ± 6.0</td><td>0.0 ± 0.0</td></tr><tr><td></td><td colspan="3"><strong>Class A</strong> &lt;&lt; Classes B, C</td></tr><tr><td>(24) Second spike, first component AHP (AHP<sub>f,2</sub>, mV)</td><td>-6.9 ± 2.6</td><td>-13.6 ± 4.2</td><td>-7.4 ± 2.8</td></tr><tr><td></td><td colspan="3"><strong>Class B</strong> &lt; Classes A, C</td></tr><tr><td>(25) Second spike, second component AHP (AHP<sub>s,2</sub>, mV)</td><td>-18.2 ± 2.2</td><td><strong>-0.7 ± 3.7</strong></td><td><strong>-23.0 ± 2.0</strong></td></tr><tr><td></td><td colspan="3"><strong>Class C</strong> &lt;&lt;&lt; <strong>Class A</strong> &lt;&lt;&lt; <strong>Class B</strong></td></tr></tbody></table>

(Continued)
<table><caption>**Table 5 | Continued**</caption><thead><tr><th></th><th>Glutamatergic class A<br>(n = 25)</th><th>Glutamatergic class B<br>(n = 12)</th><th>Glutamatergic class C<br>(n = 9)</th></tr></thead><tbody><tr><td>(26) First spike, first AHP component latency (τ<sub>AHPf</sub>, ms)</td><td>5.3 ± 1.5</td><td><strong>42.6 ± 25.5</strong></td><td>4.5 ± 0.6</td></tr><tr><td></td><td colspan="3">Classes A, C &lt;&lt;&lt; <strong>Class B</strong></td></tr><tr><td>(27) First spike, second AHP component latency (τ<sub>AHPs</sub>, ms)</td><td><strong>78.2 ± 23.7</strong></td><td>18.8 ± 33.7</td><td>0.0 ± 0.0</td></tr><tr><td></td><td colspan="3">Classes B, C &lt;&lt;&lt; <strong>Class A</strong></td></tr><tr><td>(28) Second spike, first AHP component latency (τ<sub>AHPf,2</sub>, ms)</td><td>6.1 ± 1.6</td><td><strong>57.5 ± 15.1</strong></td><td>7.7 ± 1.7</td></tr><tr><td></td><td colspan="3">Class A &lt; Class C &lt;&lt;&lt; <strong>Class B</strong></td></tr><tr><td>(29) Second spike, second AHP component latency (τ<sub>AHPs,2</sub>, ms)</td><td>73.3 ± 17.0</td><td><strong>2.2 ± 12.5</strong></td><td>95.3 ± 23.1</td></tr><tr><td></td><td colspan="3"><strong>Class B</strong> &lt;&lt;&lt; Class A &lt; Class C</td></tr><tr><td>(30) first spike ADP (ADP, mV)</td><td><strong>2.5 ± 1.8</strong></td><td>0.1 ± 0.3</td><td>0.0 ± 0.0</td></tr><tr><td></td><td colspan="3">Classes B, C &lt;&lt;&lt; <strong>Class A</strong></td></tr><tr><td>(31) Second spike ADP (ADP<sub>2</sub>, mV)</td><td>0.4 ± 0.5</td><td><strong>0.01 ± 0.07</strong></td><td>0.3 ± 0.5</td></tr><tr><td></td><td colspan="3"><strong>Class B</strong> &lt; Classes A, C</td></tr><tr><td>(32) First spike ADP latency (τ<sub>ADP</sub>, ms)</td><td><strong>10.5 ± 4.3</strong></td><td>0.6 ± 2.5</td><td>0.0 ± 0.0</td></tr><tr><td></td><td colspan="3">Classes B, C &lt;&lt;&lt; <strong>Class A</strong></td></tr><tr><td>(33) Second spike ADP latency (τ<sub>ADP2</sub>, ms)</td><td>4.6 ± 3.9</td><td>0.2 ± 1.4</td><td>5.1 ± 5.3</td></tr><tr><td></td><td colspan="3"><strong>Class B</strong> &lt; Classes A, C</td></tr></tbody></table>

Values are weighted means ± SD; n, number of cells; < significantly smaller with $P \le 0.05$; << significantly smaller with $P \le 0.01$; &lt;&lt;&lt; significantly smaller with $P \le 0.001$.

and deep monophasic AHP (Table 3) and they were able to sustain very high firing frequencies (Table 2). They expressed the GABA synthesizing enzymes GAD65 and/or 67, (GADs) and PV, frequently Calbindin (CB) and, to a lesser extent, NPY (Table 4).

Archetypal *Adapting VIP interneurons*, whose emergence followed that of FS-PV interneurons (Figure 2D), displayed a high resistance and a very small membrane capacitance (Table 1, see also cell g in Figure 4). Their spike waveform was complex with a fast AHP, followed by a marked afterdepolarization (ADP) and a secondary AHP (Table 3). At high stimulation, they showed a fast early adaptation followed by a marked late adaptation (Table 2). They expressed GADs and VIP and, relatively frequently, Calretinin (CR) (Table 4).

Archetypal *Adapting SOM interneurons* displayed a marked voltage sag in response to hyperpolarization (Table 1, see also cell k in Figures 3, 4). They were the most excitable neuronal prototype with a depolarized resting potential and a low rheobase (Table 2). Similarly to *Adapting VIP interneurons* their complex spike waveform was characterized by a prominent ADP (Table 3). At strong stimulation intensities, they very slowly adapted their firing frequency (Table 2). They expressed GADs and SOM and, often, CB and NPY (Table 4) (Kawaguchi and Kubota, 1996; Cauli et al., 2000; Karagiannis et al., 2009).

Archetypal *Adapting NPY interneurons* were weakly excitable, although less than the FS-PV prototype (Table 2, see also cell o in

Figures 3, 4), and characterized by the virtual absence of adaptation near threshold (Table 2). They fired rather wide spikes with a very low minimal firing frequency (Table 3). Repolarization behavior consisted in a biphasic AHP with a slow first component and virtually no ADP (Table 3). They expressed GADs, often NPY and, relatively frequently, the neuronal isoform of nitric oxide synthase (NOS-1) (Table 4).

## WALKS THROUGH DIVERSITY

Neurons can be represented as points in a high-dimensional space with coordinates given by their phenotype memberships $m_q$. For the purpose of visualization, positions of individual cells were plotted in different 3-dimensional projections of the full membership space (Figures 3–5). Such geometric views allowed following cell sequences with properties gradually varying between archetypes. Archetypal cells had, by definition, a single dominant membership and were thus located nearby an apex. Conversely, atypical cells were lying in the tetrahedral volume defined by three different apices.

FS-PV neurons were the first phenotype of interneurons to emerge when decreasing the fuzziness parameter (Figure 2). We first describe a sequence linking these neurons with other archetypes. Since no marked shared membership between the FS-PV and Adapting VIP archetypes could be found (as visible from Figures 2 and 6A), we only considered linking sequences
with the Adapting SOM and NPY archetypes (Figure 3A). A closed loop linking cells of the Adapting VIP, SOM, and NPY archetypes was examined in another sequence (Figure 4A).

The electrophysiological behaviors of interneurons sitting on these paths of membership space are depicted in Figures 3, 4B and their selected features are reported in tables of Figures 3, 4C. The background color of table entries encodes different ranges of feature values, which are typical for different archetypes (Tables 1–4). Accordingly columns in Figures 3, 4C corresponding to the features of cells a ($m_{FS-PV}$ = 97%), g ($m_{AdVIP}$ = 94%), k ($m_{AdSOM}$ = 79%), and o ($m_{AdNPY}$ = 83%), archetypal for FS-PV, Adapting VIP, SOM, and NPY archetypes respectively are dominated by a single color. Archetypal cells could also exhibit extreme distinctive values of certain hallmark features. For instance, cell g had a very small voltage sag and cell o a very low excitability. Conversely and indicative of transition phenotypes (Figures 3, 4C), the features of cells progressively diverging from one archetype and converging to another, gradually lose the color code of their starting archetype to adopt the color code of the target one.

## ATYPICAL CELLS RELATED WITH THE FS-PV ARCHETYPE

Figure 3 highlights transition cells on sequences linking the FS-PV to the Adapting SOM and Adapting NPY archetypes, of which neurons a, k, and o are representative examples.

Cell c is an example of a transition between FS-PV ($m_{FS-PV}$ = 29%) and Adapting SOM ($m_{AdSOM}$ = 21%) archetypes (Table 6). This neuron displayed a relatively small resistance, a high firing rate with narrow spikes and PV expression characteristic

of FS-PV cells. Unusual for FS-PV cells, it exhibited an action potential waveform with a fast and deep AHP followed by a small ADP. Also reminiscent of its secondary Adapting SOM membership, it had a slow time constant and a prominent voltage sag in response to hyperpolarizing currents. Neuron b had also transition phenotypes between the FS-PV and the Adapting SOM archetypes.

Cell e illustrates a transition between Adapting NPY ($m_{AdNPY}$ = 38%) and FS-PV ($m_{FS-PV}$ = 22%) archetypes (Table 6). Its maximal firing frequency at high stimulation intensities was moderate and it expressed NPY, consistent with its main type. Reminiscent of its secondary type it exhibited a small resistance, a low electrical excitability, fired narrow spikes and expressed PV. Most other features assumed values in ranges compliant with both the Adapting NPY and the FS-PV archetypes except for its complex repolarization phase displaying an uncommon small ADP. Similarly, neuron d had transition phenotypes between the FS-PV and the Adapting NPY archetypes.

Finally, cell f illustrates a neuron sharing multiple memberships between Adapting SOM ($m_{AdSOM}$ = 31%), FS-PV, and Adapting NPY ($m_{FS-PV}$ = $m_{AdNPY}$ = 20%) archetypes (Table 6). It had mostly intermediate properties but displayed a depolarized resting potential and voltage sag characteristic of Adapting SOM neurons, a short time constant, and deep monophasic repolarization reminiscent of its FS-PV membership. However, its AHP latency was rather long for a FS-PV neuron and was closer to that of Adapting NPY cells. Notably, it expressed the three distinctive molecular markers, PV, SOM, and NPY of its parent archetypes.

<table><caption>Table 6 | Membership values of cells in the sequences of Figures 2, 3, and 4.</caption><thead><tr><th>Cell</th><th>Glutama-tergic</th><th>FS-PV</th><th>Adapt. SOM</th><th>Adapt. VIP</th><th>Adapt. NPY</th><th>UFO</th><th>Comments</th></tr></thead><tbody><tr><td>a</td><td>0.00</td><td>0.97</td><td>0.00</td><td>0.00</td><td>0.01</td><td>0.01</td><td>Archetypal (FS-PV)</td></tr><tr><td>b</td><td>0.05</td><td>0.39</td><td>0.16</td><td>0.13</td><td>0.14</td><td>0.12</td><td>Atypical (FS-PV)</td></tr><tr><td>c</td><td>0.05</td><td>0.29</td><td>0.21</td><td>0.14</td><td>0.13</td><td>0.17</td><td>Edge cell (FS-PV—Adapting SOM)</td></tr><tr><td>d</td><td>0.05</td><td>0.40</td><td>0.15</td><td>0.07</td><td>0.21</td><td>0.12</td><td>Atypical cell (FS-PV)</td></tr><tr><td>e</td><td>0.08</td><td>0.22</td><td>0.12</td><td>0.07</td><td>0.38</td><td>0.13</td><td>Atypical cell (Adapting NPY)</td></tr><tr><td>f</td><td>0.10</td><td>0.20</td><td>0.31</td><td>0.05</td><td>0.20</td><td>0.13</td><td>Edge cell (Ad. SOM—FS-PV—Ad. NPY)</td></tr><tr><td>g</td><td>0.01</td><td>0.00</td><td>0.02</td><td>0.94</td><td>0.02</td><td>0.01</td><td>Archetypal (Adapting VIP)</td></tr><tr><td>h</td><td>0.04</td><td>0.02</td><td>0.23</td><td>0.49</td><td>0.09</td><td>0.13</td><td>Atypical (Adapting VIP)</td></tr><tr><td>i</td><td>0.06</td><td>0.02</td><td>0.40</td><td>0.26</td><td>0.11</td><td>0.15</td><td>Edge cell (Adapt. SOM—Adapt. VIP)</td></tr><tr><td>j</td><td>0.04</td><td>0.03</td><td>0.48</td><td>0.18</td><td>0.14</td><td>0.12</td><td>Atypical (Adapting SOM)</td></tr><tr><td>k</td><td>0.05</td><td>0.01</td><td>0.79</td><td>0.06</td><td>0.05</td><td>0.04</td><td>Archetypal (Adapting SOM)</td></tr><tr><td>l</td><td>0.01</td><td>0.03</td><td>0.07</td><td>0.49</td><td>0.24</td><td>0.16</td><td>Atypical (Adapting VIP)</td></tr><tr><td>m</td><td>0.06</td><td>0.05</td><td>0.14</td><td>0.34</td><td>0.22</td><td>0.19</td><td>Edge cell (Adapt. VIP—Adapt. NPY)</td></tr><tr><td>n</td><td>0.13</td><td>0.02</td><td>0.09</td><td>0.28</td><td>0.34</td><td>0.14</td><td>Edge cell (Adapt. NPY—Adapt. VIP)</td></tr><tr><td>o</td><td>0.01</td><td>0.01</td><td>0.04</td><td>0.01</td><td>0.83</td><td>0.10</td><td>Archetypal (Adapting NPY)</td></tr><tr><td>p</td><td>0.08</td><td>0.04</td><td>0.22</td><td>0.10</td><td>0.42</td><td>0.14</td><td>Atypical (Adapting NPY)</td></tr><tr><td>q</td><td>0.04</td><td>0.09</td><td>0.24</td><td>0.13</td><td>0.35</td><td>0.14</td><td>Edge cell (Adapt. NPY—Adapt. SOM)</td></tr><tr><td>r</td><td>0.10</td><td>0.13</td><td>0.33</td><td>0.09</td><td>0.22</td><td>0.13</td><td>Edge cell (Adapt. SOM—Adapt. NPY)</td></tr><tr><td>s</td><td>0.06</td><td>0.15</td><td>0.04</td><td>0.23</td><td>0.15</td><td>0.38</td><td>“UFO” main type (Burst. VIP)</td></tr><tr><td>t</td><td>0.20</td><td>0.05</td><td>0.15</td><td>0.07</td><td>0.16</td><td>0.37</td><td>“UFO” main type</td></tr><tr><td>u</td><td>0.04</td><td>0.09</td><td>0.06</td><td>0.06</td><td>0.34</td><td>0.41</td><td>“UFO” main type</td></tr></tbody></table>

Memberships evaluated at $\mu$ = 1.349, rounded at two decimals (note that, due to rounding, membership sums can deviate slightly from unit).
## ATYPICAL CELLS RELATED WITH THE ADAPTING VIP ARCHETYPE

**Figure 4** depicts another loop sequence linking the Adapting VIP, SOM, and NPY archetypes, of which neurons *g*, *k*, and *o* are representative examples.

Cell *i* is an illustrative example of a transition between Adapting VIP (*m*<sub>AdVIP</sub> = 26%) and Adapting SOM archetypes (*m*<sub>AdSOM</sub> = 40%, **Table 6**). Consistent with its Adapting SOM main type, it was fairly depolarized at rest, displayed a marked voltage sag, and slowly adapted at saturating frequencies (τ<sub>sat</sub> = 34.0 ms, see **Table 2**). Indicative of its secondary membership, its membrane capacitance was rather small and it did not express SOM but VIP. Cells *h* and *j* had also mixed Adapting VIP and SOM traits.

Neuron *m* depicts a transition cell between Adapting VIP (*m*<sub>AdVIP</sub> = 34%) and Adapting NPY (*m*<sub>AdNPY</sub> = 22%) archetypes (**Table 6**). Reminiscent of its main type, it displayed a noticeable ADP, followed by a second component of AHP. It was also able to discharge at a large firing frequency (*F*<sub>max</sub> = 120.9 Hz, cf. **Table 2**). Compatible with its secondary, but not with its main type, it exhibited a small resistance, a large rheobase and a long first spike latency (*t*<sub>1st</sub> = 478.4 ms, see **Table 2**). Cells *l* and *n* had also transition phenotypes between the Adapting VIP and NPY archetypes.

## ATYPICAL CELLS BETWEEN THE ADAPTING SOM AND NPY ARCHETYPES

In addition to the atypical cells with noticeable Adapting SOM or NPY memberships already described in the previous two sections and related as well with the FS-PV or Adapting VIP archetypes, **Figure 4** shows other atypical cells lying between the Adapting SOM and NPY archetypes.

As an example, cell *q* had Adapting NPY main type (*m*<sub>AdNPY</sub> = 35%) and a second important Adapting SOM membership (*m*<sub>AdSOM</sub> = 24%, **Table 6**). Consistent with its main type it was weakly excitable and expressed NPY. It displayed a complex spike waveform with uncommon features such as a slow ADP, suggestive of an Adapting SOM phenotype, combined with a marked and fast second component of AHP, characteristic of Adapting NPY cells. In agreement with its secondary but, intriguingly, not with its main type (**Tables 1–4**), cell *q* had an intermediate minimal firing frequency and expressed CB. Cells *p* and *r* also exhibited mixed Adapting NPY and SOM traits.

## OTHER ATYPICAL CELLS

Finally, some highly atypical cells displayed heterogeneous phenotypes that could not be simply described as transitional between defined archetypes. All these cells had a large membership toward the residual UFO class. Neurons *s*, *t*, and *u* are illustrative example UFOs illustrated in **Figure 5**.

Despite being GABAergic, Cell *t* (*m*<sub>UFO</sub> = 37%) displayed anomalous features, e.g., broad spikes and an inability to discharge at high firing rates, which are reminiscent of the Glutamatergic archetype (*m*<sub>Glut</sub> = 20%, **Table 6**). It shared as well properties with the FS-PV (small spike amplitude), the Adapting SOM (depolarized resting membrane potential and intermediate first AHP amplitude), and the Adapting NPY

archetypes (slow first AHP component) and it had a membrane resistance extreme even for the Adapting VIP archetype.

Few cells with UFO main type (*n* = 18 out of 33) displayed-marked bursting phenotype, defined by the inability to fire below a minimal frequency (Hodgkin, 1948) and characterized by strong near-threshold adaptation and spike amplitude reduction at the onset of firing. This hinted at the potential emergence of a further Bursting VIP archetype (Karagiannis et al., 2009) out of the heterogeneous and yet indistinct collection of cells constituting the UFO class. Cell *s* represents here a characteristic example of this tendency (**Figure 5**) and was indeed classified as belonging to the Bursting VIP type in the previous crisp analysis (Karagiannis et al., 2009).

## GRADED AND STRUCTURED SEPARATIONS BETWEEN ARCHETYPES

These partial lists of atypical cells suggest that, with the exception of Glutamatergic, FS-PV, and Adapting VIP archetypes (**Figure 2**), a neat separation between archetypes may be lacking. **Figure 6A** allows visualizing in a systematic way, through a series of 2-dimensional cross-sections, how neurons were scattered in the membership space. The two axes of every triangular panel in **Figure 6A** correspond to memberships toward a different pair of archetypes for each panel. The black dashed diagonal line corresponds to equal membership toward the two considered archetypes and it is surrounded by a shaded gray region corresponding to the halved mean absolute deviation. Cells falling within this region had therefore almost identical memberships to the two archetypes and were lying precisely at the geometric frontier between them. Consequently, we denoted them as *edge cells* (see Materials and Methods for definition).

There were no Glutamatergic cells at the edge with well-defined archetypes (i.e., any archetype apart from the UFO class), corresponding to a complete separation of the Glutamatergic archetype from all the other GABAergic archetypes. The separation between the FS-PV and the Adapting VIP archetypes, which were the first two clusters of GABAergic neurons emerging with decreasing value of μ (Figures 2B,C,E), was also very distinct and there was no edge cell between them.

Edge cells were found however at the fringe of other Adapting Archetypes. Considering Adapting SOM neurons, there were edge cells between the Adapting SOM and the FS-PV archetypes (*n* = 3) and between the Adapting SOM and the Adapting VIP archetypes (*n* = 5). Regarding Adapting NPY neurons, there were edge cells between the Adapting NPY and the FS-PV archetypes (*n* = 2) and between the Adapting NPY and the Adapting VIP archetypes (*n* = 4). There was also a conspicuous overlap between the two Adapting SOM and NPY archetypes themselves (*n* = 8 edge cells). Finally, every archetype had cells at the edge with the UFO class, but these cells could not easily be described as transitional, since there is no clear phenotype associated to the UFO class.

To achieve an accurate and synthetic description of cell distribution in a high-dimensional space, we quantified for each cell its distance from the boundary between its main and secondary archetypes by computing its *typicality coefficient d*. As defined in the *Materials and Methods*, cells with typicalities larger or smaller than specific thresholds, related to the halved mean
absolute deviation of the typicality distribution, were considered respectively as archetypal or atypical. The chosen thresholds captured well the bimodal structure of the sample distribution of *d* at μ = 1.349, in which peaks associated to strictly archetypal and atypical cells were clearly visible, as shown in **Figure 6B**.

According to the used criteria, 40% of the cells (n = 80 of 200) resulted as archetypal and 42% (n = 83 out of 200) as atypical. However, the distribution of archetypal and atypical cells was not homogeneous across different archetypes. The Glutamatergic, the FS-PV and, to a lesser extent, the Adapting VIP archetypes, were well-detached from the other archetypes, as indicated by their large average typicality coefficients and fractions of strictly archetypal cells (**Figure 6**). Other GABAergic archetypes displayed larger strictly atypical fractions. UFOs were an exceptional case, since none of them were archetypal (see **Figure 6B** and **Table 6**).

We note that the actual distribution of typicality coefficients depended on the chosen value of the fuzziness parameter. However, since the relative degrees of typicality of different cells were very well-preserved over a wide range of parameters, the fractions of cells labeled as atypical or archetypal were robustly conserved (not shown).

## THE RELEVANCE OF DIFFERENT FEATURES FOR CLASSIFICATION

The archetypes and the typicality spectrum observed in our sample reflected the existence of structured correlations between the values of different features in different cells. A first natural question is therefore establishing which of the measured features are the most decisive for the determination of the extracted fuzzy partition. A second question is whether it would be possible to use fewer features and, yet, obtain an equivalent classification.

Generalizing an approach we previously introduced (Karagiannis et al., 2009), we assessed the relevance for classification of different properties or groups of properties by means of partial randomization analyses. The values of a subset of properties were “scrambled” by permuting them randomly between different cells. Such a manipulation preserved by construction the sample average and standard deviation of the scrambled features, but destroyed their correlation with other features. Scrambling properties of higher significance to the observed fuzzy partition will lead to a lower classification quality. We quantified such variations in quality of fuzzy clustering by comparing the *cost function* of the extracted fuzzy partitions (see Materials and Methods) before and after randomization.

The increase in classification cost determined by randomization of individual features varied in the 2–4% range. Properties are ranked in **Table 7** in descending order of cost increase after individual randomization. GAD expression (feature 35, **Table 4**), discriminating excitatory and inhibitory neurons, was the property associated to the largest cost increase. It was, followed by the maximum steady state frequency $F_{max}$ (feature 15, **Table 3**), distinctly high for FS-PV neurons and very low for Glutamatergic neurons (**Table 2**). Note that the first two archetypes to emerge when decreasing the fuzziness μ were the Glutamatergic and the FS-PV archetypes.

We next sought to determine whether it was possible to reproduce the reference classification with a parsimonious list

<table><caption>**Table 7 | Relevance of properties for the classification of our sample.**</caption><thead><tr><th>Rank</th><th>Feature</th><th>ΔJ (%)</th></tr></thead><tbody><tr><td>#1</td><td>(35) GAD</td><td>4.80 ± 0.60</td></tr><tr><td>#2</td><td>(15) Maximal steady state frequency</td><td>4.46 ± 0.70</td></tr><tr><td>#3</td><td>(22) First spike, first component AHP</td><td>4.45 ± 0.75</td></tr><tr><td>#4</td><td>(18) First spike duration</td><td>4.41 ± 0.72</td></tr><tr><td>#5</td><td>(19) Second spike duration</td><td>4.41 ± 0.65</td></tr><tr><td>#6</td><td>(11) Amplitude accommodation</td><td>4.21 ± 0.68</td></tr><tr><td>#7</td><td>(24) Second spike, first component AHP</td><td>4.01 ± 0.76</td></tr><tr><td>#8</td><td>(27) First spike, second AHP component latency</td><td>3.96 ± 0.70</td></tr><tr><td>#9</td><td>(32) First spike ADP latency</td><td>3.69 ± 0.70</td></tr><tr><td>#10</td><td>(31) Second spike ADP</td><td>3.65 ± 0.81</td></tr><tr><td>#11</td><td>(34) VGluT1</td><td>3.63 ± 0.77</td></tr><tr><td>#12</td><td>(29) Second spike, second AHP component latency</td><td>3.59 ± 0.74</td></tr><tr><td>#13</td><td>(30) first spike ADP</td><td>3.55 ± 0.80</td></tr><tr><td>#14</td><td>(4) Time constant</td><td>3.54 ± 0.78</td></tr><tr><td>#15</td><td>(14) Late adaptation</td><td>3.46 ± 0.76</td></tr><tr><td>#16</td><td>(41) VIP</td><td>3.44 ± 0.76</td></tr><tr><td>#17</td><td>(38) PV</td><td>3.37 ± 0.82</td></tr><tr><td>#18</td><td>(33) Second spike ADP latency</td><td>3.35 ± 0.86</td></tr><tr><td>#19</td><td>(21) Duration Increase</td><td>3.35 ± 0.67</td></tr><tr><td>#20</td><td>(23) First spike, second component AHP</td><td>3.33 ± 0.72</td></tr><tr><td>#21</td><td>(12) Amplitude of early adaptation</td><td>3.33 ± 0.84</td></tr><tr><td>#22</td><td>(7) Rheobase</td><td>3.31 ± 0.75</td></tr><tr><td>#23</td><td>(16) First spike amplitude</td><td>3.27 ± 0.80</td></tr><tr><td>#24</td><td>(25) Second spike, second component AHP</td><td>3.20 ± 0.82</td></tr><tr><td>#25</td><td>(5) Membrane capacitance</td><td>3.19 ± 0.78</td></tr><tr><td>#26</td><td>(17) Second spike amplitude</td><td>3.18 ± 0.71</td></tr><tr><td>#27</td><td>(20) Amplitude reduction</td><td>3.16 ± 0.66</td></tr><tr><td>#28</td><td>(3) Input Resistance</td><td>3.10 ± 0.85</td></tr><tr><td>#29</td><td>(26) First spike, first component AHP latency</td><td>3.07 ± 0.65</td></tr><tr><td>#30</td><td>(40) NPY</td><td>3.06 ± 0.79</td></tr><tr><td>#31</td><td>(9) Just-above threshold adaptation</td><td>3.06 ± 0.63</td></tr><tr><td>#32</td><td>(6) Sag index</td><td>3.03 ± 0.76</td></tr><tr><td>#33</td><td>(8) First spike latency</td><td>2.94 ± 0.76</td></tr><tr><td>#34</td><td>(28) Second spike, first component AHP latency</td><td>2.89 ± 0.98</td></tr><tr><td>#35</td><td>(42) SOM</td><td>2.88 ± 0.76</td></tr><tr><td>#36</td><td>(2) Resting Potential</td><td>2.85 ± 0.80</td></tr><tr><td>#37</td><td>(1) Digitized laminar location</td><td>2.85 ± 0.70</td></tr><tr><td>#38</td><td>(10) Minimal steady state frequency</td><td>2.76 ± 0.69</td></tr><tr><td>#39</td><td>(43) CCK</td><td>2.73 ± 0.83</td></tr><tr><td>#40</td><td>(13) Time constant of early adaptation</td><td>2.69 ± 0.75</td></tr><tr><td>#41</td><td>(39) CR</td><td>2.69 ± 0.95</td></tr><tr><td>#42</td><td>(37) CB</td><td>2.58 ± 0.78</td></tr><tr><td>#43</td><td>(36) NOS1</td><td>2.45 ± 0.91</td></tr></tbody></table>

Increase of the residual cost of the fuzzy partition estimated via single-property randomization. Mean values and sample standard deviations are taken over 1000 independent randomizations.

of features (**Figure 7**). We, thus, generated different classifications, starting with the most relevant property (i.e., the expression of GAD according to the ranking in **Table 7**) and then gradually adding new properties in decreasing order of relevance (see Materials and Methods). We found that while certain archetypes could be well-discriminated by a small number of top-ranked
<figure><img src="image_16.png" /><figcaption>**FIGURE 7** | Relevance of features for classification. Matching between classifications based on a reduced number of top-ranked properties with the reference classification based on all the 43 features (see **Table 7** for relevance ranking of the different properties). Classification matching is analyzed separately for every well-defined archetype and is measured by the fraction of cells with a given main type, matching in both the reference, and a reduced classification. Matching classification fraction for all archetypes confounded is also shown for comparison. Note the order of correctly classified archetypes corresponding to the historically characterized neuronal types.</figcaption></figure>

features, the correct discrimination of other archetypes required the consideration of a much larger number of properties. About half of the neurons with Glutamatergic or FS-PV main type could be identified by considering only the top three-ranked properties, i.e., GAD expression, maximum steady state frequency $F_{max}$ and first AHP component of the first spike $AHP_f$ (feature 22, **Table 3**). Considering also the duration $D_1$ of the first and $D_2$ of the second spikes (features 18–19, **Table 3**), allowed the correct classification of over half of the neurons with Adapting VIP main type and raised the fraction of correctly identified FS-PV neurons to over 70%. At least nine features were required for the identification of more than half of the cells with Adapting SOM main type, and at least 13 features for the identification of more than half of the cells with Adapting NPY main type. Correct identification of over 90% of all cells in the dataset was only possible by considering at least 35 features, the 30th and the 35th of which being respectively the molecular markers NPY and SOM.

Interestingly, this observed order of correctly classified archetypes is similar to the order of emergence of archetypes obtained by decreasing the fuzziness parameter (Figure 2) as well as with their historical identification (Mountcastle et al., 1969; McCormick et al., 1985; Kawaguchi and Kubota, 1993, 1996; Karagiannis et al., 2009).

It should be however noted that the results of these property relevance analyses are strongly dependent on the fine structure of our sample. Furthermore, the effects of randomizing more than one property simultaneously cannot be predicted trivially

by the knowledge of the effects of randomization of individual properties. For instance, individual randomization of SOM and NPY expression have a small overall impact on classification, but their joint randomization can strongly affect the correct detection of Adapting SOM and NPY archetypes (not shown). More systematic analyses over multiple independent samples would thus be required to make general statements on the impact on classification of different properties or property groups.

# DISCUSSION

## A STRUCTURED CONTINUUM

Neuronal features varied in a remarkably graded way. However, this landscape of diversity was far from being a continuum without structure, since definite trends and trend changes were evident along archetype transitions (Figures 3–4).

Our approach allowed describing this structured continuum of phenotypes, by mediating between contrasting philosophies of neuronal classification (Tyner, 1975; Parra et al., 1998). On the one hand, we avoided the provocative view of the existence of a separate class for each cell (Parra et al., 1998). On the other hand, we could preserve information about the observed phenotypic diversity, which subdivisions into sharply distinct classes (Markram et al., 2004; Burkhalter, 2008) would have curtailed.

Unlike possibly more powerful but unavoidably more abstract dimensionality reduction techniques (Bishop, 2006; Hinton and Salakhutdinov, 2006), fuzzy memberships had a straight biological relevance, representing directly and compactly the similarity between single cells and canonic template archetypes. Because of an intuitive spatial representation of sample diversity, the analysis of memberships allowed thus the systematic detection of atypical and archetypal cells within a dataset already too large to be parsed without neuroinformatic assistance.

## THE "RIGHT" NUMBER OF ARCHETYPES

The four reference GABAergic archetypes and the typicality spectrum which we found at $\mu = 1.349$ were robust in a wide parameter range. They were also optimal, in the sense that any attempt to reliably extract more archetypes from our sample was unsuccessful. Indeed, for $\mu < 1.204$, even tiny changes in $\mu$ shattered entirely the extracted fuzzy partition (Figure 2E). Introducing a metaphor, one might say that adopting a too small $\mu$ would be like looking at a poor quality printed image through a macro lens. This would reveal most likely the low resolution of the halftone screen, rather than additional relevant details. Conversely, using a too large $\mu$ would be like taking an out-of-focus snapshot (Figure 2F). Therefore, our present list of archetypes should be considered, citing again Tyner (1975, p. 91), as “a temporary convenience, a rather crude data storage system,” reflecting only our present state of knowledge (and ignorance).

The failure to extract more archetypes was likely due to the lack of supplemental properties allowing subdivision of neuronal populations (Tyner, 1975). Such features might correspond to morphology and connectivity patterns (Markram et al., 2004; Christophe et al., 2005; Krimer et al., 2005; Perin et al., 2011), synaptic activity (Dumitriu et al., 2007), temporal dynamics of activity (Klausberger and Somogyi, 2008), and/or mechanisms of synaptic plasticity (Szabo et al., 2012). The analysis
summarized by **Figure 7** revealed that a smaller number of features was required to correctly discriminate the archetypes that emerged earlier at a high value of μ. The consideration of additional features in our classification would therefore most likely allow the discrimination of additional neuronal types, not yet emerged out of the residual UFO class. The simple results summarized by **Figure 7** and **Table 7** constitute only the first steps of a systematic study of the impact of different properties on neuronal classification, which goes beyond the limits of the present study.

Finally, since specific selection criteria (i.e., superficial layers and prominent radial processes) were used to increase the probability of collecting NPY neurons (Karagiannis et al., 2009), our results are not necessarily representative of the entire population of cortical interneurons. Nevertheless, in a complementary study aimed at characterizing deep layers neurons not necessarily exhibiting radial processes (Perrenoud et al., 2013), very similar classes of interneurons were obtained. Therefore, we believe that our approach provides at least a comprehensive view of interneuronal types, even if it does not claim to give an absolute and definitive classification of them.

## FROM DEVELOPMENT TO GRADED DIVERSITY

Neuronal phenotypes are largely determined by developmental processes, which depend not only on the embryonic origin and date of birth (Butt et al., 2005; Miyoshi et al., 2007) but also on the signaling environment (Batista-Brito and Fishell, 2009) or the ongoing network activity (Cossart, 2011).

Glutamatergic neurons have embryonic origins, dates of birth, and migration patterns very different from that of interneurons (Wonders and Anderson, 2006). They were the first archetype to singularize and had no edge cell with the four well-defined interneuron archetypes (Figures 2, 6). Similarly, the FS-PV and the Adapting VIP archetypes exhibiting also clearly distinct embryonic origins (Butt et al., 2005; Miyoshi et al., 2007; Lee et al., 2010; Vucurovic et al., 2010) were the two first interneuron archetypes to emerge and no edge cell could be found between them. These observations suggest that neuronal archetypes that are primarily specified by clearly distinct spatial and temporal origins generally display very different phenotypic traits (Wonders and Anderson, 2006). Such multiple distinctive features result in an early singularization and therefore in a limited number of edge cells.

Neurons initially sharing similar origins progressively specify during development to form distinct neuronal classes. FS-PV and Adapting SOM neurons, both derived from the medial ganglionic eminence (Wonders and Anderson, 2006; Miyoshi et al., 2007), provide an example of such a developmental divergence. The observation of edge cells between these archetypes could speculatively reflect a vestige of their shared spatial embryonic origin. Alternatively, since neurons were collected from juvenile rats (P19) the existence of edge cells between FS-PV and Adapting SOM neurons (e.g., cell b) which exhibited a relatively high input resistance could also correspond to maturating FS-PV neurons. Indeed, these neurons progressively reduce their input resistance by acquiring two-pore K⁺ leak conductances (Okaty et al., 2009).

Conversely, neurons of distinct embryonic origins can also converge to a similar functional class. Adapting NPY neurons provide such an example of developmental convergence. Although they display common electrophysiological, molecular, and morphological features (Povysheva et al., 2008; Gelman et al., 2009; Karagiannis et al., 2009; Zaitsev et al., 2009; Lee et al., 2010) they remarkably originate either from the medial or caudal ganglionic eminences, or the preoptic area (Gelman et al., 2009; Tricoire et al., 2010; Vucurovic et al., 2010).

## FROM LOCAL CUES TO GRADED DIVERSITY

The overall converging phenotypic homogeneity of Adapting NPY neurons is however mitigated by the differential expression of NOS-1 (Cauli et al., 2004; Karagiannis et al., 2009; Kubota et al., 2011; Perrenoud et al., 2012a), which in the hippocampus (Tricoire et al., 2010; Jaglin et al., 2012), but not in the neocortex (Magno et al., 2012), is dependent on the embryonic origins.

This hints at other potential determinants of graded diversity. The density of cortical NOS-1 interneurons is enriched in the vicinity of blood vessels (Cauli et al., 2004; Tricoire et al., 2010; Perrenoud et al., 2013; Rockland and Nayyar, 2012). Given the shared local cues governing neuritic branching and vascular sprouting (Stubbs et al., 2009; Adams and Eichmann, 2010), it is likely that the local vascular environment might also orientate phenotypic traits.

More generally, an appealing hypothesis deserving future investigation is that certain atypical traits could emerge as a phenotypic adaptation to context-dependent signaling. The acquisition of secondary functional specializations would contribute to blur inter-archetype boundaries.

## A RESERVOIR OF PHENOTYPES FOR FUNCTIONAL DEGENERACY

An important motivation underlying many studies of neuronal diversity is to establish a correspondence between neuronal types and their respective functions. Under this perspective, the diversity of neuronal types would correspond tightly to the multiplicity of possible functional specializations (Burkhalter, 2008). Such a view is however oversimplified under many aspects.

First, relatively homogeneous groups of neurons can exert multiple functions. For instance, the overall population of Adapting NPY/neurogliaform neurons (Karagiannis et al., 2009; Zaitsev et al., 2009) is responsible for slow GABAergic inhibition (Tamás et al., 2003; Szabadics et al., 2007; Oláh et al., 2009). Since NPY and NO are respectively potent vasoconstrictor and vasodilator of cortical blood vessels (Cauli et al., 2004; Cauli and Hamel, 2010), the Adapting NPY archetype might also be specialized in a complex control of blood perfusion (Estrada and DeFelipe, 1998).

Second, different neuronal populations can contribute and/or cooperate to a same function. The neurovascular coupling is indeed mediated by multiple vasoactive messengers produced by different cell types (Attwell et al., 2010; Cauli and Hamel, 2010). Several neuronal types, including NOS-1 and VIP-expressing interneurons (Cauli et al., 2004; Perrenoud et al., 2013) as well as COX-2 expressing pyramidal cells (Niwa et al., 2000; Wang et al., 2005; Lecrux et al., 2011) are the primary sources of these vasodilatory compounds. Such a partial functional overlap
between different neuronal types confers robustness to the neu-
rovascular response (Leithner et al., 2010), essential for normal
brain functions (Iadecola, 2004).

Therefore, the will to maintain a one-to-one matching between
neuronal types and functions would ineluctably lead to the defi-
nition of a combinatorial number of neuronal types (Parra et al.,
1998).

The multiplicity of cell types devoted to the neurovascul-
lar response constitutes an example of functional degeneracy
(Edelman, 1987; Tononi et al., 1999), defined as the ability of
heterogeneous elements to perform the same function. Beyond
redundancy, occurring when a given function is achieved by repli-
cating homogeneous elements, degeneracy confers higher robust-
ness through adaptability. Indeed, heterogeneous elements can
react differently in different contexts providing a considerable
margin of safety over a wide spectrum of conditions.

A RESERVOIR OF PHENOTYPES FOR FUNCTIONAL COMPENSATION

Another interesting question raised by our systematic exploration
is how specific functions could be reliably achieved by networks
whose cellular components are heterogeneous in a structured and
graded manner.

A same circuit function can arise from very diverse com-
binations of neuronal properties and interaction mechanisms
(Prinz et al., 2004; Goaillard et al., 2009). Indeed, heterogeneity
of cellular properties does not necessarily translate to different

performance at the network level for at least two reasons: First
many of these properties might be irrelevant for the function of
the network, and, in absence of specific constraints, they would
freely assume values over broad ranges. Second, even variations
of certain critical properties in a cell might be compensated by
covariations in other cells in the network. These two aspects
might contribute to the intra-archetype diversity emphasized
by our analysis. In addition, compensation between heteroge-
neous elements might provide a more flexible and robust solution
than optimization of individual elements to achieve network
functions (Marder and Goaillard, 2006). In this context the exis-
tence of a structured continuum of diversity might be viewed
as a prerequisite for the emergence of reliable operation at the
system level.

ACKNOWLEDGMENTS

The authors thank Bertrand Lambolez, Bâlint Lasztòczi, Jochen Staiger, Philippe Faure, and Ludovic Tricoire for useful discussions. We acknowledge financial support by the German Federal Ministry of Education and Research (BMBF) via the Bernstein Center for Computational Neuroscience—Göttingen (01GQ1005B), by the Human Frontier Science Program (RGY0070/2007) and by the Agence Nationale pour la Recherche (ANR 2011 MALZ 003 01). Anastassios Karagiannis was supported by a Fondation pour la Recherche Médicale fellowship (FDT20100920106).

REFERENCES

Adams, R. H., and Eichmann, A. (2010). Axon guidance molecules in vascular patterning. *Cold Spring Harb. Perspect. Biol.* 2:a001875. doi: 10.1101/cshperspect.a001875

Andjelic, S., Gallopin, T., Cauli, B., Hill, E. L., Roux, L., Badr, S., et al. (2009). Glutamatergic nonpyramidal neurons from neocortical layer VI and their comparison with pyramidal and spiny stellate neurons. *J. Neurophysiol.* 101, 641–654.

Ascoli, G. A., Alonso-Nanclares, L., Anderson, S. A., Barrionuevo, G., Benavides-Piccione, R., Burkhalter, A., et al. (2008). Petilla Terminology, Nomenclature of features of GABAergic interneurons of the cerebral cortex. *Nat. Rev. Neurosci.* 9, 557–568.

Attwell, D., Buchan, A. M., Charpak, S., Lauritzen, M., MacVicar, B. A., and Newman, E. A. (2010). Glial and neuronal control of brain blood flow. *Nature* 468, 232–243.

Basu, S., Banerjee, A., and Mooney, R. J. (2002). “Semi-supervised clustering by seeding,” in *Proceedings of the Nineteenth International Conference on Machine Learning*, eds C. Sammut and A. G. Hoffmann (San Francisco, CA: Morgan Kaufmann), 27–34.

Batista-Brito, R., and Fishell, G. (2009).
The developmental integration of

cortical interneurons into a func-
tional network. *Curr. Top. Dev. Biol.*
87, 81–118.

Bezdek, J. C. (1974). Numerical taxonomy with fuzzy sets. *J. Math. Biol.* 1, 57–71.

Bishop, C. M. (2006). *Pattern Recognition and Machine Learning*. New York, NY: Springer.

Burkhalter, A. (2008). Many specialists for suppressing cortical excitation. *Front. Neurosci.* 2:2. doi: 10.3389/neuro.01.026.2008

Butt, S. J., Fuccillo, M., Nery, S., Noctor, S., Kriegstein, A., Corbin, J. G., et al. (2005). The temporal and spatial origins of cortical interneurons predict their physiological subtype. *Neuron* 48, 591–604.

Cauli, B., Audinat, E., Lambolez, B., Angulo, M. C., Ropert, N., Tsuzuki, K., et al. (1997). Molecular and physiological diversity of cortical nonpyramidal cells. *J. Neurosci.* 17, 3894–3906.

Cauli, B., and Hamel, E. (2010). Revisiting the role of neurons in neurovascular coupling. *Front. Neuroenerg.* 2:9. doi: 10.3389/fnene.2010.00009

Cauli, B., Porter, J. T., Tsuzuki, K., Lambolez, B., Rossier, J., Quenet, B., et al. (2000). Classification of fusiform neocortical interneurons based on unsupervised clustering.

Proc. Natl. Acad. Sci. U.S.A. 97, 6144-6149.

Cauli, B., Tong, X. K., Rancillac, A., Serluca, N., Lambolez, B., Rossier, J., et al. (2004). Cortical GABA interneurons in neurovascular coupling: relays for subcortical vasoactive pathways. *J. Neurosci.* 24, 8940–8949.

Chagnac-Amitai, Y., Luhmann, H. J., and Prince, D. A. (1990). Burst generating and regular spiking layer 5 pyramidal neurons of rat neocortex have different morphological features. *J. Comp. Neurol.* 296, 598–613.

Christophe, E., Doerflinger, N., Lavery, D. J., Molnar, Z., Charpak, S., and Audinat, E. (2005). Two populations of layer V pyramidal cells of the mouse neocortex: development and sensitivity to anesthetics. *J. Neurophysiol.* 94, 3357–3367.

Connors, B. W., and Gutnick, M. J. (1990). Intrinsic firing patterns of diverse neocortical neurons. *Trends Neurosci.* 13, 99–104.

Cossart, R. (2011). The maturation of cortical interneuron diversity: how multiple developmental journeys shape the emergence of proper network function. *Curr. Opin. Neurobiol.* 21, 160–168.

Dumitriu, D., Cossart, R., Huang, J., and Yuste, R. (2007). Correlation between axonal morphologies and

synaptic input kinetics of interneu-
rons from mouse visual cortex.
Cereb. Cortex 17, 81–91.

Dunn, J. C. (1973). A fuzzy relative of the ISODATA process and its use in detecting compact well-separated clusters. *J. Cybern.* 3, 32–57.

Edelman, G. M. (1987). *Neural Darwinism: The Theory of Neuronal Group Selection*. New York, NY: Basic Books.

Equiha, M. (1990). Fuzzy clustering of ecological data. *J. Ecol.* 78, 519–534.

Estrada, C., and DeFelipe, J. (1998). Nitric oxide-producing neurons in the neocortex: morphological and functional relationship with intraparenchymal microvasculature. *Cereb. Cortex* 8, 193–203.

Gallopin, T., Geoffroy, H., Rossier, J., and Lambolez, B. (2006). Cortical sources of CRF, NKB, and CCK and their effects on pyramidal cells in the neocortex. *Cereb. Cortex* 16, 1440–1452.

Gelman, D. M., Martini, F. J., Nóbrega-Pereira, S., Pierani, A., Kessaris, N., and Marin, O. (2009). The embryonic preoptic area is a novel source of cortical GABAergic interneurons. *J. Neurosci.* 29, 9380–9389.

Gentet, L. J. (2012). Functional diversity of supragranular GABAergic neurons in the barrel cortex. *Front. Neural Circuits* 6:52. doi: 10.3389/fncir.2012.00052
Goaillard, J. M., Taylor, A. L., Schulz, D. J., and Marder, E. (2009). Functional consequences of animal-to-animal variation in circuit parameters. *Nat. Neurosci.* 12, 1424–1430.

Gupta, A., Wang, Y., and Markram, H. (2000). Organizing principles for a diversity of GABAergic interneurons and synapses in the neocortex. *Science* 287, 273–278.

Halabisky, B., Shen, F., Huguenard, J. R., and Prince, D. A. (2006). Electrophysiological classification of somatostatin-positive interneurons in mouse sensorimotor cortex. *J. Neurophysiol.* 96, 834–845.

Helmstaedter, M., Sakmann, B., and Feldmeyer, D. (2009). The relation between dendritic geometry, electrical excitability, and axonal projections of L2/3 interneurons in rat barrel cortex. *Cereb. Cortex* 19, 938–950.

Hinton, G. E., and Salakhutdinov, R. R. (2006). Reducing the dimensionality of data with neural networks. *Science* 313, 504–507.

Hodgkin, A. L. (1948). The local electric changes associated with repetitive action in a non-medullated axon. *J. Physiol.* 107, 165–181.

Iadecola, C. (2004). Neurovascular regulation in the normal brain and in Alzheimer's disease. *Nat. Rev. Neurosci.* 5, 347–360.

Jaglin, X. H., Hjerling-Leffler, J., Fishell, G., and Batista-Brito, R. (2012). The origin of neocortical nitric oxide synthase-expressing inhibitory neurons. *Front. Neural Circuits* 6:44. doi: 10.3389/fncir.2012.00004

Jang, J. S. R., Sun, C. T., and Mizutani, E. (1997). *Neuro-Fuzzy and Soft Computing*. Upper Saddle River, NJ: Prentice Hall.

Karagiannis, A., Gallopin, T., Dávid, C., Battaglia, D., Geoffroy, H., Rossier, J., et al. (2009). Classification of NPY-expressing neocortical interneurons. *J. Neurosci.* 29, 3642–3659.

Kätzel, D., Zemelman, B. V., Buetfering, C., Wölfel, M., and Miesenböck, G. (2011). The columnar and laminar organization of inhibitory connections to neocortical excitatory cells. *Nat. Neurosci.* 14, 100–107.

Kawaguchi, Y., and Kubota, Y. (1993). Correlation of physiological subgroupings of nonpyramidal cells with parvalbumin- and calbindinD28k-immunoreactive neurons in layer V of rat frontal cortex. *J. Neurophysiol.* 70, 387–396.

Kawaguchi, Y., and Kubota, Y. (1996). Physiological and morphological identification of somatostatin- or

vasoactive intestinal polypeptide-containing cells among GABAergic cell subtypes in rat frontal cortex. *J. Neurosci.* 16, 2701–2715.

Klausberger, T., and Somogyi, P. (2008). Neuronal diversity and temporal dynamics: the unity of hippocampal circuit operations. *Science* 321, 53–57.

Krimer, L. S., Zaitsev, A. V., Czanner, G., Kröner, S., González-Burgos, G., Povysheva, N. V., et al. (2005). Cluster analysis-based physiological classification and morphological properties of inhibitory neurons in layers 2-3 of monkey dorsolateral prefrontal cortex. *J. Neurophysiol.* 94, 3009–3022.

Kubota, Y., Hattori, R., and Yui, Y. (1994). Three distinct subpopulations of GABAergic neurons in rat frontal agranular cortex. *Brain Res.* 649, 159–173.

Kubota, Y., Shigematsu, N., Karube, F., Sekigawa, A., Kato, S., Yamaguchi, N., et al. (2011). Selective coexpression of multiple chemical markers defines discrete populations of neocortical GABAergic neurons. *Cereb. Cortex* 21, 1803–1817.

Lecrux, C., Toussay, X., Kocharyan, A., Fernandes, P., Neupane, S., Lévesque, M., et al. (2011). Pyramidal neurons are “neurogenic hubs” in the neurovascular coupling response to whisker stimulation. *J. Neurosci.* 31, 9836–9847.

Lee, S., Hjerling-Leffler, J., Zagha, E., Fishell, G., and Rudy, B. (2010). The largest group of superficial neocortical GABAergic interneurons expresses ionotropic serotonin receptors. *J. Neurosci.* 30, 16796–16808.

Leithner, C., Royl, G., Offenhauser, N., Füchtemeier, M., Kohl-Bareis, M., Villringer, A., et al. (2010). Pharmacological uncoupling of activation induced increases in CBF and CMRO2. *J. Cereb. Blood Flow Metab.* 30, 311–322.

Ma, Y., Hu, H., Berrebi, A. S., Mathers, P. H., and Agmon, A. (2006). Distinct subtypes of somatostatin-containing neocortical interneurons revealed in transgenic mice. *J. Neurosci.* 26, 5069–5082.

Magno, L., Oliveira, M. G., Mucha, M., Rubin, A. N., and Kessaris, N. (2012). Multiple embryonic origins of nitric oxide synthase-expressing GABAergic neurons of the neocortex. *Front. Neural Circuits* 6:65. doi: 10.3389/fncir.2012.00065

Marder, E., and Goaillard, J. M. (2006). Variability, compensation and homeostasis in neuron and network function. *Nat. Rev. Neurosci.* 7, 563–574.

Markram, H., Toledo-Rodriguez, M., Wang, Y., Gupta, A., Silberberg, G., and Wu, C. (2004). Interneurons of the neocortical inhibitory system. *Nat. Rev. Neurosci.* 5, 793–807.

McCormick, D. A., Connors, B. W., Lighthall, J. W., and Prince, D. A. (1985). Comparative electrophysiology of pyramidal and sparsely spiny stellate neurons of the neocortex. *J. Neurophysiol.* 54, 782–806.

McGarry, L. M., Packer, A. M., Fino, E., Nikolenko, V., Sippy, T., and Yuste, R. (2010). Quantitative classification of somatostatin-positive neocortical interneurons identifies three interneuron subtypes. *Front. Neural Circuits* 4:12. doi: 10.3389/fncir.2010.00012

Miyoshi, G., Butt, S. J., Takebayashi, H., and Fishell, G. (2007). Physiologically distinct temporal cohorts of cortical interneurons arise from telencephalic Olig2-expressing precursors. *J. Neurosci.* 27, 7786–7798.

Mountcastle, V. B., Talbot, W. H., Sakata, H., and Hyvarinen, J. (1969). Cortical neuronal mechanisms in flutter-vibration studied in unanesthetized monkeys. Neuronal periodicity and frequency discrimination. *J. Neurophysiol.* 32, 452–484.

Niwa, K., Araki, E., Morham, S. G., Ross, M. E., and Iadecola, C. (2000). Cyclooxygenase-2 contributes to functional hyperemia in whisker-barrel cortex. *J. Neurosci.* 20, 763–770.

Okaty, B. W., Miller, M. N., Sugino, K., Hempel, C. M., and Nelson, S. B. (2009). Transcriptional and electrophysiological maturation of neocortical fast-spiking GABAergic interneurons. *J. Neurosci.* 29, 7040–7052.

Oláh, S., Fule, M., Komlósi, G., Varga, C., Báldi, R., Barzó, P., et al. (2009). Regulation of cortical microcircuits by unitary GABA-mediated volume transmission. *Nature* 461, 1278–1281.

Parra, P., Gulyás, A. I., and Miles, R. (1998). How many subtypes of inhibitory cells in the hippocampus? *Neuron* 20, 983–993.

Perin, R., Berger, T. K., and Markram, H. (2011). A synaptic organizing principle for cortical neuronal groups. *Proc. Natl. Acad. Sci. U.S.A.* 108, 5419–5424.

Perrenoud, Q., Geoffroy, H., Gauthier, B., Rancillac, A., Alfonsi, F., Kessaris, N., et al. (2012a). Characterization of type I and type II nNOS-expressing interneurons in the barrel cortex of mouse.

Front. Neural Circuits 6:36. doi: 10.3389/fncir.2012.00036

Perrenoud, Q., Rossier, J., Férézou, I., Geoffroy, H., Gallopin, T., Vitalis, T., et al. (2012b). Activation of cortical 5-HT(3) receptor-expressing interneurons induces NO mediated vasodilatations and NPY mediated vasoconstrictions. *Front. Neural Circuits* 6:50. doi: 10.3389/fncir.2012.00050

Perrenoud, Q., Rossier, J., Geoffroy, H., Vitalis, T., and Gallopin, T. (2013). Diversity of GABAergic interneurons in layer V1a and V1b of mouse barrel cortex. *Cereb. Cortex.* 23, 423–441.

Povysheva, N. V., Zaitsev, A. V., Rotaru, D. C., González-Burgos, G., Lewis, D. A., and Krimer, L. S. (2008). Parvalbumin-positive basket interneurons in monkey and rat prefrontal cortex. *J. Neurophysiol.* 100, 2348–2360.

Prinz, A. A., Bucher, D., and Marder, E. (2004). Similar network activity from disparate circuit parameters. *Nat. Neurosci.* 7, 1345–1352.

Rockland, K. S., and Nayyar, N. (2012). Association of type I neurons positive for NADPH-dihydrolase with blood vessels in the adult monkey corpus callosum. *Front. Neural Circuits* 6:4. doi: 10.3389/fncir.2012.00004

Rousseeuw, P. J. (1987). Silhouettes: a graphical aid to the interpretation and validation of cluster analysis. *J. Comput. Appl. Math.* 20, 53–65.

Sachs, L. (1984). *Applied Statistics: A Handbook of Techniques*. New York, NY: Springer.

Salski, A. (2007). Fuzzy clustering of fuzzy ecological data. *Ecol. Inform.* 3, 262–269.

Sohal, V. S., Zhang, F., Yizhar, O., and Deisseroth, K. (2009). Parvalbumin neurons and gamma rhythms enhance cortical circuit performance. *Nature* 459, 698–702.

Somogyi, P., and Klausberger, T. (2005). Defined types of cortical interneurons structure space and spike timing in the hippocampus. *J. Physiol.* 562, 9–26.

Staiger, J. F., Flagmeyer, I., Schubert, D., Zilles, K., Kotter, R., and Luhmann, H. J. (2004). Functional diversity of layer IV spiny neurons in rat somatosensory cortex: quantitative morphology of electrophysiology characterised and biocytin labeled cells. *Cereb. Cortex* 14, 690–701.

Stubbs, D., DeProto, J., Nie, K., Englund, C., Mahmud, I., Hevner, R., et al. (2009). Neurovascular congruence during cerebral
cortical development. *Cereb. Cortex* 19(Suppl. 1), i32–i41.

Suzuki, N., and Bekkers, J. M. (2010). Distinctive classes of GABAergic interneurons provide layer-specific phasic inhibition in the anterior piriform cortex. *Cereb. Cortex* 20, 2971–2984.

Szabadics, J., Tamás, G., and Soltesz, I. (2007). Different transmitter transients underlie presynaptic cell type specificity of GABAA, slow and GABAA, fast. *Proc. Natl. Acad. Sci. U.S.A.* 104, 14831–14836.

Szabo, A., Somogyi, J., Cauli, B., Lambolez, B., Somogyi, P., and Lamsa, K. P. (2012). Calcium-permeable AMPA receptors provide a common mechanism for LTP in glutamatergic synapses of distinct hippocampal interneuron types. *J. Neurosci.* 32, 6511–6516.

Tamás, G., Buhl, E. H., and Somogyi, P. (1997). Fast IPSPs elicited via multiple synaptic release sites by different types of GABAergic neurone in the cat visual cortex. *J. Physiol.* 500, 715–738.

Tamás, G., Lorincz, A., Simon, A., and Szabadics, J. (2003). Identified sources and targets of

slow inhibition in the neocortex. *Science* 299, 1902–1905.

Taylor, J. R. (1997). *An Introduction to Error Analysis, 2nd Edn*. Sausalito, CA: University Science Books.

Toledo-Rodriguez, M., Blumenfeld, B., Wu, C., Luo, J., Attali, B., Goodman, P., et al. (2004). Correlation maps allow neuronal electrical properties to be predicted from single-cell gene expression profiles in rat neocortex. *Cereb. Cortex* 14, 1310–1327.

Tononi, G., Sporns, O., and Edelman, G. M. (1999). Measures of degeneracy and redundancy in biological networks. *Proc. Natl. Acad. Sci. U.S.A.* 96, 3257–3262.

Tricoire, L., Pelkey, K. A., Daw, M. I., Sousa, V. H., Miyoshi, G., Jeffries, B., et al. (2010). Common origins of hippocampal Ivy and nitric oxide synthase expressing neurogliaform cells. *J. Neurosci.* 30, 2165–2176.

Tyner, C. F. (1975). The naming of neurons: applications of taxonomic theory to the study of cellular populations. *Brain Behav. Evol.* 12, 75–96.

Vucurovic, K., Gallopin, T., Ferezou, I., Rancillac, A., Chameau, P., van

Hooft, J. A., et al. (2010). Serotonin 3A receptor subtype as an early and protracted marker of cortical interneuron subpopulations. *Cereb. Cortex* 20, 2333–2347.

Wang, H., Hitron, I. M., Iadecola, C., and Pickel, V. M. (2005). Synaptic and vascular associations of neurons containing cyclooxygenase-2 and nitric oxide synthase in rat somatosensory cortex. *Cereb. Cortex* 15, 1250–1260.

Wonders, C. P., and Anderson, S. A. (2006). The origin and specification of cortical interneurons. *Nat. Rev. Neurosci.* 7, 687–696.

Woodruff, A. R., Anderson, S. A., and Yuste, R. (2010). The enigmatic function of chandelier cells. *Front. Neurosci.* 4:201. doi: 10.3389/fnins.2010.00201

Xu, D., Keller, J. M., Popescu, M., and Bondugula, R. (2008). *Applications of Fuzzy Logic in Bioinformatics*. London: Imperial College Press.

Zadeh, L. A. (1965). Fuzzy sets. *Inform. Control* 8, 338–353.

Zadeh, L. A. (2008). Is there a need for fuzzy logic? *Inf. Sci.* 178, 2751–2779.

Zaitsev, A. V., Povysheva, N. V., Gonzalez-Burgos, G., Rotaru, D., Fish, K. N., Krimer, L. S., et al.

(2009). Interneuron diversity in layers 2-3 of monkey prefrontal cortex. *Cereb. Cortex* 19, 1597–1615.

**Conflict of Interest Statement:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Received: 19 November 2012; paper pending published: 09 December 2012; accepted: 21 January 2013; published online: 07 February 2013.

Citation: Battaglia D, Karagiannis A, Gallopin T, Gutch HW and Cauli B (2013) Beyond the frontiers of neuronal types. *Front. Neural Circuits* 7:13. doi: 10.3389/fncir.2013.00013

Copyright © 2013 Battaglia, Karagiannis, Gallopin, Gutch and Cauli. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in other forums, provided the original authors and source are credited and subject to any copyright notices concerning any third-party graphics etc.