Author Response: Mapping Cortical Mesoscopic Networks of Single Spiking Cortical or Sub-Cortical Neurons
Dongsheng Xiao,Matthieu P. Vanni,Catalin Mitelut,Allen W. Chan,Jeffrey LeDue,Yicheng Xie,Andrew Cn Chen,Nicholas V. Swindale,Timothy H. Murphy
DOI: https://doi.org/10.7554/elife.19976.035
2017-01-01
Abstract:Article Figures and data Abstract Introduction Results Discussion Materials and methods References Decision letter Author response Article and author information Metrics Abstract Understanding the basis of brain function requires knowledge of cortical operations over wide-spatial scales, but also within the context of single neurons. In vivo, wide-field GCaMP imaging and sub-cortical/cortical cellular electrophysiology were used in mice to investigate relationships between spontaneous single neuron spiking and mesoscopic cortical activity. We make use of a rich set of cortical activity motifs that are present in spontaneous activity in anesthetized and awake animals. A mesoscale spike-triggered averaging procedure allowed the identification of motifs that are preferentially linked to individual spiking neurons by employing genetically targeted indicators of neuronal activity. Thalamic neurons predicted and reported specific cycles of wide-scale cortical inhibition/excitation. In contrast, spike-triggered maps derived from single cortical neurons yielded spatio-temporal maps expected for regional cortical consensus function. This approach can define network relationships between any point source of neuronal spiking and mesoscale cortical maps. https://doi.org/10.7554/eLife.19976.001 Introduction Neural activity ranges from the microscale of synapses to macroscale brain-wide networks. Mesoscale networks occupy an intermediate space and are well studied in cortex forming the basis of sensory and motor maps (Bohland et al., 2009). These networks are largely defined by co-activation of neurons and have been evaluated with a variety of statistical approaches that capitalize on detecting synchrony. The study of large-scale networks (meso-to macroscale) has been mostly restricted to functional magnetic resonance imaging (fMRI), or magnetoencephalography that can capture whole-brain activity patterns (de Pasquale et al., 2010; Kahn et al., 2011; Logothetis et al., 2012), but lack high spatial and temporal resolution and sensitivity. To overcome these limitations, alternative approaches including mesoscopic intrinsic signal, voltage, glutamate, or calcium sensitive indicator imaging have been employed (Kleinfeld et al., 1994; Kenet et al., 2003; Ferezou et al., 2007; Chemla and Chavane, 2010; Chen et al., 2013b; Mohajerani et al., 2013; Stroh et al., 2013; Vanni and Murphy, 2014; Carandini et al., 2015; Chan et al., 2015; Madisen et al., 2015; Wekselblatt et al., 2016; Xie et al., 2016). New preparations using large-scale craniotomies (Kim et al., 2016b) and large format imaging systems (Tsai et al., 2015; Sofroniew et al., 2016) provide the ability to link mesoscale activity patterns to individual neurons. However, these measures are restricted to superficial layers of cortex and cannot assess linkages with sub-cortical structures. Developments in fiberoptic technology allow local optical functional assessment of brain activity in sub-cortical structures (Hamel et al., 2015; Kim et al., 2016a), but cannot simultaneously resolve cortex over large fields of view. Although the evolution of imaging has revealed new aspects of cortical processing in identified neurons (Harvey et al., 2012; Chen et al., 2013a, Chen et al., 2013b; Fu et al., 2014; Guo et al., 2014), the electrically recorded action potential is still a signal of prominence given its exquisite timing and ability to reflect the output of neuronal networks (Buzsáki, 2004). We combine extracellular recordings of single units in the cortex, thalamus, and other sub-cortical sites with mesoscopic functional imaging in transgenic mice expressing the calcium indicator GCaMP (Zariwala et al., 2012; Vanni and Murphy, 2014; Silasi et al., 2016). While slower than protein-based or small molecule voltage sensors, GCaMP imaging offers a high signal-to-noise ratio and is associated with supra-threshold activity which allows a more direct comparison with spike activity. This work extends pioneering studies investigating the relationship between single neuron spiking and local neuronal population activity assessed by voltage-sensitive dye imaging. Spike-triggered averaging (STA) was used to identify the local activity profile related to the spiking activity of a single neuron within this population (Arieli et al., 1995), and it was further demonstrated that this activity profile could reveal the instantaneous spatial pattern of ongoing population activity related to a neuron’s optimal stimulus in visual cortex of anesthetized cat (Tsodyks et al., 1999). This current study extends these approaches and also exploits the main advantage of mesoscopic imaging allowing the simultaneous measurement of brain activity in multiple regions across most of cortex simultaneously and not only that of the local population of activity surrounding the recording site. This multiscale strategy has allowed us to define temporal relationships between the activity of single neurons at the microscopic scale and mesoscale cortical maps (Zingg et al., 2014; Madisen et al., 2015). Furthermore, we employ multisite, long shank, silicon probe recordings of single neuron activity that facilitates the assessment of long-distance activity relationships between remote subcortical single neuron activity and mesoscale cortical population activity. Spontaneous activity in awake and anesthetized mice was exploited as a source of diverse cortical network activity motifs (Mohajerani et al., 2010, 2013; Chan et al., 2015). Application of STA to cortical spontaneous activity linked single neurons to mesoscale networks. Single thalamic neuron spikes were found to functionally link to multiple primary sensorimotor maps, in contrast spiking cortical neurons were largely associated with consensus cortical maps. Thalamic neurons were found to both predict and report (firing before and after) specific cycles of wide-scale cortical inhibition/excitation, while cortical neuron firing was usually associated with excitation. These results are consistent with an active computational role of thalamus in sensory-motor processing (Theyel et al., 2010; Hooks et al., 2013; Petrus et al., 2014; Sheroziya and Timofeev, 2014; McCormick et al., 2015), as opposed to merely serving a relay function and is consistent with a diverse role of the thalamus in feedforward sensory processing. Thalamocortical transmission can dynamically and differentially recruit local cortical excitation and inhibition based on thalamic neuron firing patterns and where thalamocortical feedforward inhibition is a critical feature (Galarreta and Hestrin, 1998; Swadlow and Gusev, 2001; Gabernet et al., 2005; Cruikshank et al., 2007; Hu and Agmon, 2016). We expect that this spike-triggered cortical mapping technique, exploiting mesoscopic calcium imaging, can be extended to any brain location where electrodes can be placed to identify functionally linked cortical mesoscale networks. Results Linkage of individual spiking neurons to specific mesoscopic cortical maps We exploit the wide field of view of mesoscale cortical imaging using GCaMP transgenic mice (Madisen et al., 2015) in combination with cellular electrophysiology recordings to derive cortical networks that reflect activity at targeted point sources of neuronal spiking throughout the brain. Cortical and sub-cortical neuron spiking activities were recorded electrically while simultaneously imaging cortical mesoscopic activity across a 9 × 9 mm bilateral window that encompassed multiple areas of the mouse dorsal cortex including somatosensory, motor, visual, retrosplenial, parietal association and cingulate areas (Figure 1A,B). Spectral decomposition of the mesoscopic spontaneous activity using GCaMP6 revealed the presence of information below 10 Hz that was distinct from non-specific green light reflectance (Figure 1—figure supplement 1). Given the slow Ca2+ binding and unbinding kinetics of GCaMP6, we expect imaging dynamics will be prolonged compared to actual spike records. In some cases, we employed deconvolution (Pnevmatikakis et al., 2016) to improve the time course of raw calcium signals (Figure 1—figure supplement 2). While deconvolution improved the temporal dynamics of the decay of the calcium signal, spike-triggered analysis was only marginally affected and it was not used throughout. Spiking signals were initially recorded in multiple brain areas using glass electrodes (n = 8 mice) to minimize obstruction of cortical imaging and reduce potential for damage from electrode placement. Subsequently, laminar probes (16 channel with 0.1 mm contact spacing) permitted the resolution of more spiking neurons simultaneously, and facilitated the recordings in multiple sub-cortical regions (n = 16 mice). Given the invasive nature and the long duration of recordings, initial data were obtained from urethane (n = 4) or isoflurane (n = 12) anesthetized adult mice, but were later optimized to include awake recordings (n = 12, see Materials and methods). The spike-triggered average maps (STMs) obtained under both these conditions were qualitatively similar, and this observation was consistent with previous work using VSD imaging (Mohajerani et al., 2013). To perform these assessments, we identified single neuron spikes from extracellular recordings using spike sorting methods based on clustering of principal components distributions of spike signals on sets of adjacent channels (Swindale and Spacek, 2014) (Figure 1C). Figure 1 with 2 supplements see all Download asset Open asset Experimental setup and multichannel electrode recordings and spike classification. (A) Set-up for simultaneous wide-field calcium imaging and single unit recording using a glass pipette or laminar silicon probe. (Bi) Top view of wide-field transcranial window and (ii) cortical atlas adapted from the Allen Institute Brain Atlas. (C) Example of (i) cortical and (ii) subcortical pairs or spike recordings from separate channels showing the isolation in the two principal components axes. (D) The generation of a spike-triggered average map (STM) for unit located in barrel cortex. (Ei) STM generated from single neuron with 1158 spikes recorded in right barrel cortex. (ii) Red traces: Spontaneous calcium activity recorded from two different cortical areas (BCS1 and HLS1). Blue trace: spontaneous spiking activity recorded simultaneously from right BCS1. (iii) STMTD generated from average of calcium activity time-locked with each spike (red) and random spike (see Materials and methods, black, blue: subtraction of spike and random spike-evoked responses) in region-of-interest (ROI). These examples results were from mice under anesthesia. Source files for the generation of spike-triggered average map (STM) can be found at http://goo.gl/nHF29I. The folder ‘Matlab code and source data’ (see the sublink of https://github.com/catubc/sta_maps) contains calcium images (‘tif’ file), spike train and Matlab code used for the generation of STM shown in Figure 1E. The ‘tif’ file cannot not be viewed with a standard picture viewer, but must be viewed with a program, such as ‘ImageJ’. Spike times were exported as ‘txt’ file. Matlab code (named ‘STA_eLife.m’) was used for reading images and spike time files and generating STM. https://doi.org/10.7554/eLife.19976.002 To investigate how single neuron spiking activity at a cortical/sub-cortical site was related to regional cortical activity, we calculated STMs from simultaneously acquired wide-field calcium imaging (Figure 1D; see also Materials and methods). For each individual spike, we considered cortical image frames from 3s before to 3s after the spike normalized as ΔF/F0 by subtracting and dividing the average calcium activity during the 3s preceding the spike. The static STM was then defined as the peak response (in units of ΔF/F0) calculated for each pixel within a time window of ±1 s. This peak method is a better reflection of correlated activity than the average over the ±1 s interval which would smooth out highly activated – but short duration activity increases. This method revealed that the activity recorded from a single right barrel cortex neuron yielded an STM showing strong and specific GCaMP signal in the barrel and motor cortices of both hemispheres (Figure 1Ei). STMs were thus calculated by averaging calcium activity during spiking activity (Figure 1Eii) and revealed the high spatial specificity of the mapping when compared with reference region (hind limb) or random spike averaging (see Materials and methods) (Figure 1Eiii). We verified the calcium specificity of STMs (reflecting underlying neuronal activity) by imaging Thy-1 GFP-M mice (n = 6 mice) that lacked calcium-dependent neuronal fluorescent signals and failed to produce functional maps using the same procedures (Figure 2A,B). To investigate the sensitivity of the technique, the minimum number of spikes needed to make maps in GCaMP mice was measured by quantifying the similarities of pairs of STMs generated from a subset of increasing numbers of spikes. Stable STM maps were generally observed using 256 spikes (Figure 2C–E). A high stability of STMs was also confirmed by comparing the maps generated by splitting a unit's spikes into two halves, or into odd and even groups which yield similar STMs. Figure 2 Download asset Open asset Sensitivity and specificity of STMs. (A) Simultaneous calcium and spiking activity recording in GCaMP6f mouse and STM yielded from single unit recorded in barrel cortex. (B) Simultaneous GFP fluorescence and spiking activity recording in Thy-1 GFP-M mouse and STM yielded from single unit recorded in barrel cortex resulted in no clear regional map. (C) STMs generated from a subset of spikes (2–2048, on the left) randomly chosen in one experiment. Correlation coefficients (r-value on the right) between STMs were used to evaluate the consistency of mapping. In this example, STMs generated by more than 64 spikes generated a correlation >0.9 and revealed high similarity between the pairs of SPMs made using the same number of spikes. (D) Distribution of correlation values between pairs of STM for an increasing number of spikes. No significant change in r-value distribution was observed for 512 spikes in comparison to 256 spikes (Mann Whitney test, p=0.126, U = 948.5, 256 spikes group n = 58, r-value = 0.97 ± 0.01, mean ± SD; 512 spikes group n = 40, r-value = 0.98 ± 0.01). (E) STMs and profile of responses computed using spikes divided into halves or even-odd sets. These examples were performed under anesthesia. https://doi.org/10.7554/eLife.19976.005 Thalamic neurons show more diverse STMs than cortical neurons By combining wide-field calcium imaging and single unit recording, the mesoscopic network associated with any neuron of the brain can be mapped. After establishing the method, we then focused on thalamic and cortical recordings of spiking neurons and confirmed their anatomical location by labeling probes with Texas red-dextran or DiI to visualize tracks (Figure 3A, four experiments subcortical track and nucleus identification, see also Figure 2A,B). We observed that spiking cortical neurons were linked to consensus local and long-range cortical networks as defined by the stability of STMs between recordings (Figure 3B). We define consensus cortical networks as those that can be identified by assessment of correlated activity (seed pixel analysis) and reflect major mono-synaptic intra-cortical axonal projections (Mohajerani et al., 2013). Spiking barrel cortex neurons were consistently linked to regional GCaMP signal changes in the barrel and motor cortex, as well as showing signals in homotopic areas of both hemispheres, consistent with the previously observed pattern of long-distance connections (Ferezou et al., 2007; Mohajerani et al., 2013; Guo et al., 2014; Vanni and Murphy, 2014; Chan et al., 2015). To compare multiple neuron STMs, we displayed the ‘contour’ of each STM as the full-width-half-max value (Figure 3B,C-contours). Figure 3 with 2 supplements see all Download asset Open asset Topographic properties of cortical and thalamic STM. (A) Electrode track for each recording (Blue channel: DAPI, yellow: DiI). (B) STM and overlay contours of neurons recorded in barrel cortex in the same animal. Each color represents one STM contour. (C) STM and overlay contours of neurons recorded in thalamus in the electrode track presented in panels A and B. Color bar in the right side indicated the depth of each recording site. (D) Diversity of overlap of STMs between neurons on neighboring laminar electrode channels. (i) Example of overlapping STMs (red area) between two cortical neurons recorded on adjacent channels. (ii) Example of overlapping STMs for neighboring pairs of neurons recorded subcortically showing differences across depth. (iii) Average neighboring cortical neuron map overlap (blue: 93%) and neighboring sub-cortical neuron overlap (78%) show significant differences (Mann Whitney test, p<0.0001, U = 617408.0, mean percentage overlap of cortical STM pairs = 92.77% ± 0.23%, mean ± SEM, n = 966; mean percentage overlap of sub-cortical STM pairs = 78.11% ± 0.61%, mean ± SEM, n = 1936). These results are from awake mice (except Mouse #1). https://doi.org/10.7554/eLife.19976.006 We have also assessed other means of generating event triggered maps including Multi-Unit-Activity (MUA) and local field potential (LFP) frequency bands. Cortical and thalamic STMs computed from single unit activity or MUA strongly resembled one another and did not vary greatly according to laminar depth. LFP-triggered STMs for delta band activity were similar to STMs, while those associated with higher frequency bands showed more unique patterns that will be investigated in future work (Figure 3—figure supplement 1 and Figure 3—figure supplement 2). Delta band activity is where most mesoscale power functional imaging indicator power is located (Chan et al., 2015). In contrast, in the case of GCaMP6 higher frequency components are closer to hemodynamic and other noise sources making analysis more challenging. STMs for single thalamic neurons indicated not only a functional link between GCaMP maps and their consensus cortical projection areas (Hunnicutt et al., 2014; Oh et al., 2014; Zingg et al., 2014) but also showed more variability and complexity in behavior than for cortical units. Thalamic neurons were associated with both unilateral and bilateral hemispheric signals within multiple primary sensorimotor and higher order brain areas (Figure 3C). To quantify the difference of variability of STM between cortical and thalamic units, we compared the percentage overlap of static STMs (Figure 3D; red areas) for neighboring pairs of neurons (100 μm apart). We found cortical neurons in the same functional region exhibited substantial similarity, while subcortical neurons had more diversity even within the same sub-nucleus in thalamus (i.e. VPL, VPM, see Table 1 for nomenclature and compare 3 example in Figure 3C). The use of multichannel probes allowed us to obtain spiking profiles across cortical layers. However, there were only subtle topographic changes: the contours of static STMs were largely similar between and within superficial and deep layers, respectively (see the contour map in Figure 3B). To assess diversity in single neuron spiking derived cortical maps (STMs), we compared STMs derived from neighboring electrode contacts and found that sub-cortical derived maps were more varied (Figure 3Diii) Table 1 Abbreviation used to define different cortical/sub-cortical areas. https://doi.org/10.7554/eLife.19976.009 S1Primary somatosensory areaS2Supplemental or Secondary Somatosensory areaFLForelimb region of the Primary Somatosensory area (FLS1)HLHindlimb region of the Primary Somatosensory area (HLS1)BCBarrel region of the Primary Somatosensory area (BCS1)M1Primary motor areaM2Secondary motor areaMOMouth region of the Primary Somatosensory areaNONose region of the Primary Somatosensory areaTRTrunk region of the Primary Somatosensory areaUN(Unassigned) region of the Primary Somatosensory area (S1)ACAnterior Cingulate area (ACC)AAnterior or Posterior Partial Association areas: PTLp or PTAV1Primary visual cortexALAnteroLateral regions of the extrastriate visual areasAMAnteroMedial regions of the extrastriate visual areasLMLateralMedial regions of the extrastriate visual areasPLPosteroLateral regions of the extrastriate visual areasLILateralIntermediate regions of the extrastriate visual areasPMPosteroMedial regions of the extrastriate visual areasPORPostrhinal regions of the extrastriate visual areasRLRostroLateral regions of the extrastriate visual areasAUPrimary Auditory areaTEATemporal Association areaRSRetrosplenial areaPTAParietal Association areaVPMVentral posteromedial nucleus of the thalamusVPLVentral posterolateral nucleus of the thalamusPOPosterior complex of the thalamusRTReticular nucleus of the thalamusLGNLateral geniculate nucleusCPCaudoputamenHPFHippocampal formation Sub-cortical neurons are linked to cortical maps not predicted from consensus networks To determine quantitatively how cortical/sub-cortical STMs were related to intra-cortical networks, static STMs were compared using cross-correlation with a cortex-wide library of seed pixel correlation maps (SPM) (Mohajerani et al., 2013; Vanni and Murphy, 2014; Chan et al., 2015) generated iteratively for all locations from the same recording of spontaneous activity (Figure 4A,B). To create SPMs, the cross-correlation coefficient r values between the temporal profiles of one selected pixel and all the others within the field of view were calculated (Mohajerani et al., 2013; Vanni and Murphy, 2014; Chan et al., 2015) (see Materials and methods). To evaluate the similarity between the static STM and SPMs, we calculated the correlation coefficient between pixels of both types of maps for all possible SPMs in the library. We then selected the SPM that resulted in best match: highest correlation between a given STM and the library of SPMs. The library of SPMs is expected to reflect cortical consensus activity motifs (areas undergoing temporally-correlated activity) and can be largely attributed to their underlying intra-cortical axonal projections (Mohajerani et al., 2013). Single cortical neuron-derived STMs were largely predicted by the pattern of cortical connectivity using SPM correlation mapping (Figure 4A) as correlations were relatively high between these SPMs and STMs (Figure 4C). In contrast, thalamic STMs were more complex and corresponded to more unique distributions of cortical patterns and had significantly lower correlations with the cortical consensus SPM library (Figure 4B,C). For example, VPM, VPL and CP neurons can functionally link to multiple cortical areas that are not predicted by SPM (SPM made by putting the seeds in either BCS1, HLS1 or RS areas). In other words, it is possible that the cortical STMs derived from single spiking sub-cortical neurons are the super-position of 2 or more cortical networks defined by SPMs. To support this theory, we show examples comparing subcortical STMs to pairs of SPMs made from two seed locations that can lead to these potentially more complex maps (Figure 4D). For example, the VPM STM could be constructed as a combination of a BCS1 and RS SPM (seed point a and b, respectively, Figure 4D). Figure 4 Download asset Open asset STM compared with seed pixel correlation maps (SPM). (A) Cortical STM (left) and the best fitting SPM (right) according to correlation coefficient (cc) values for different electrode placements (text to the left of panel). Similarity was calculated by measuring the r-value Pearson coefficient between each pair of map pixels (in title). Group data from 12 GCaMP6f mice are reported in panel C. (B) Sub-cortical STM (left) and the most similar SPM (right). Cells #2 and #4 were from GCaMP6s mice. Cells #3, 5–7, 11 were from GCaMP3 mice. Other cells were from GCaMP6f mice. These examples were performed under anesthesia. (C) Distribution of r-values (Mann Whitney test, p<0.0001, U = 5227, sub-cortical group n = 246 r-value=0.64 ± 0.18, mean SD; cortical group n = 168 r-value=0.85 ± 0.04, mean ± SD). (D) Examples of sub-cortical STMs compared with pairs of SPMs for seed indicated by ‘a’ and ‘b’. https://doi.org/10.7554/eLife.19976.010 To better understand the underlying structural circuit basis of distinct thalamic STMs, we examined the Allen Mouse Brain Connectivity Atlas (Oh et al., 2014) as in our previous work (Figure 5) (Mohajerani et al., 2013). We processed the three-dimensional structural data and matched the composed anatomical 2D maps with our two-dimensional static STMs. As expected, we found that STMs of spiking cortical neurons correspond with underlying structural axonal projections (see BCS1 example in Figure 5, [Mohajerani et al., 2013]). However, sub-cortical STMs cannot be predicted by direct monosynaptic projections from sub-cortical to cortical areas. For example, the HPF has no strong direct structural link to RS area. Furthermore, CP was not directly linked to BCS1/HLS1 areas (Hunnicutt et al., 2014; Oh et al., 2014), indicating that sub-cortical STMs reflect apparently polysynaptic links to cortex. Figure 5 Download asset Open asset Examples of STMs and projection maps. (A) Example of STMs from neurons recorded in BCS1, VPM, VPL, HPF and CP. (B) Example of projection maps (2D surface and 3D) reconstructed from Allen Brain Atlas with injection sites (Oh et al., 2014) in the same region as our recording. For example, for a spiking neuron recorded in BCS1 we show anterograde labeling of GFP emanating from an injection site in BCS1 that extends to motor cortex and is present across cortex in the 2D surface plot of cortex. This projection pattern for BCS1 matches the STM map quite well as in previous work (Mohajerani et al., 2013). In contrast, for sub-cortical injections of GFP tracer such as in HPF there was less overlap between STMs and projection maps perhaps indicating polysynaptic pathways. Website: 2015 Allen Institute for Brain Science. Allen Mouse Brain Connectivity Atlas [Internet]. Available from: http://connectivity.brain-map.org. https://doi.org/10.7554/eLife.19976.011 Cortical and sub-cortical neuron firing is tuned to cortical network dynamics spanning millisecond to multi-second time scales By analyzing the cortical GCaMP signal, time courses corresponding to single neuron firing (firing= time 0) we observed dynamic cortical activity states (Figures 6 and 7; Videos 1–4). We determined STM Temporal Dynamics (STMTD) by identifying a region-of-interest (ROI) and tracking the time course of the maximally activated/depressed cortical pixel of the region from 3s before to 3s after spiking. As all our extracellular recordings were in right barrel cortex and predominantly right sensory thalamus, we identified the left barrel cortex as a co-activated area of interest and tracked dynamics within this ROI for spiking neurons. As observed for static STMs, the cortical calcium dynamics associated with spiking cortical neurons were relatively homogeneous with an initial peak in activity within ~100 ms – 200 ms following spiking and a return to baseline (Figure 6A). However, some cortical cells (~20%) participated in multi-second depression dynamics (see distribution of profiles in Figures 7B,C and and 8A). In contrast, STMTDs generated by thalamic cells were more varied and were dominated by depression dynamics (~80% of cells) lasting up to 3 s (Figures 6B and 7B,C and Figure 8B). Figure 6 with 6 supplements see all Download asset Open asset Montages of cortical and thalamic spatio-temporal dynamics. (A) Top: right hemisphere barrel cortex neuron spiking time montage stereotyped dynamics in left hemisphere barrel cortex region-of-interest (ROI). The maximally activated pixel in the ROI (red arrow) is tracked over time and reveals Spike-Triggered-Map Temporal Dynamics (STMTD) which rises quickly at spike time t = 0 and decays in 100–200 ms followed by 1–2 s cortical depression (red curve in right plot). Bottom: Additional examples of right hemisphere barrel cortex neuron spike-triggered montages show similar barrel-motor cortex activation pattern with peaks in cortical activation shortly following spiking and a return to baseline or prolonged depression. (B) Top: Same as in A, but for a thalamically recorded neuron (right hemisphere) which correlates strongly with motor cortex activation shortly after spiking. Bottom: Additional montage examples of thalamic neurons (also right hemisphere) reveal both the spatial diversity (i.e. different STMs) and temporal diversity (i.e. different STM dynamics). https://doi.org/10.7554/eLife.19976.012 Figure 7 with 1 supplement see all Download asset Open asset Classification of STMTD patterns. (A) Example of two STMTDs from pixels within L-BCS1 (left barrel cortex) from a single mouse recording (both cortical and subcortical neurons STMTD pixel locations shown). Each recorded cell STMTD has a slightly different maximum pixel amplitude location, but all fall within the L-BCS1 region. (B) STMTD PCA distribution from all 428 cortical and subcortical neurons recorded from all mice separated using KMEANS (k = 3). (Ci) STMTD patterns (±SD) classifications from (B). The number of neurons from cortex and thalamus used