This is the concluding article in a series of 3 studies that investigate the anatomical determinants of thalamocortical (TC) input to excitatory neurons in a cortical column of rat primary somatosensory cortex (S1). We used viral synaptophysin-enhanced green fluorescent protein expression in thalamic neurons and reconstructions of biocytin-labeled cortical neurons in TC slices to quantify the number and distribution of boutons from the ventral posterior medial (VPM) and posteromedial (POm) nuclei potentially innervating dendritic arbors of excitatory neurons located in layers (L)2–6 of a cortical column in rat somatosensory cortex. We found that 1) all types of excitatory neurons potentially receive substantial TC input (90–580 boutons per neuron); 2) pyramidal neurons in L3–L6 receive dual TC input from both VPM and POm that is potentially of equal magnitude for thick-tufted L5 pyramidal neurons (ca. 300 boutons each from VPM and POm); 3) L3, L4, and L5 pyramidal neurons have multiple (2–4) subcellular TC innervation domains that match the dendritic compartments of pyramidal cells; and 4) a subtype of thick-tufted L5 pyramidal neurons has an additional VPM innervation domain in L4. The multiple subcellular TC innervation domains of L5 pyramidal neurons may partly explain their specific action potential patterns observed in vivo. We conclude that the substantial potential TC innervation of all excitatory neuron types in a cortical column constitutes an anatomical basis for the initial near-simultaneous representation of a sensory stimulus in different neuron types.

Highlights d The consequences of activating single pyramidal cells were studied in turtle cortex d Two spikes from single neurons activated tens of neurons sequentially and reliably d The activated neurons were both interneurons and other pyramidal cells d The evoked sequences lasted up to 200 ms, spread over 1 mm 2 , and never ignited cortex

Understanding circuit computation in the nervous system requires sampling activity over large neural populations and maximizing the number of features that can be extracted. By combining planar arrays of extracellular electrodes with the three-layered cortex of turtles, we show that synaptic signals induced along individual axons as well as action potentials can be easily captured. Two types of information can be extracted from these signals, the neuronal subtype (inhibitory or excitatory)-whose identification is more reliable than with traditional measures such as action potential width-and a (partial) spatial map of functional axonal projections from individual neurons. Because our approach is algorithmic, it can be carried out in parallel on hundreds of simultaneously recorded neurons. Combining our approach with soma triangulation, we reveal an axonal projection bias among a population of pyramidal neurons in turtle cortex and confirm this bias through anatomical reconstructions.

Optimizing spike-sorting algorithms is difficult because sorted clusters can rarely be checked against independently obtained “ground truth” data. In most spike-sorting algorithms in use today, the optimality of a clustering solution is assessed relative to some assumption on the distribution of the spike shapes associated with a particular single unit (e.g., Gaussianity) and by visual inspection of the clustering solution followed by manual validation. When the spatiotemporal waveforms of spikes from different cells overlap, the decision as to whether two spikes should be assigned to the same source can be quite subjective, if it is not based on reliable quantitative measures. We propose a new approach, whereby spike clusters are identified from the most consensual partition across an ensemble of clustering solutions. Using the variability of the clustering solutions across successive iterations of the same clustering algorithm (template matching based on K-means clusters), we estimate the probability of spikes being clustered together and identify groups of spikes that are not statistically distinguishable from one another. Thus, we identify spikes that are most likely to be clustered together and therefore correspond to consistent spike clusters. This method has the potential advantage that it does not rely on any model of the spike shapes. It also provides estimates of the proportion of misclassified spikes for each of the identified clusters. We tested our algorithm on several datasets for which there exists a ground truth (simultaneous intracellular data), and show that it performs close to the optimum reached by a support vector machine trained on the ground truth. We also show that the estimated rate of misclassification matches the proportion of misclassified spikes measured from the ground truth data.

No abstract