4.1.

Large Spatial Extent of the Point Spread: Further Characterization

It is important to note that for the vast majority of optical imaging studies, the PS and, in particular, its entire spatial extent was rarely the intended subject of research. Indeed, it has been fortuitous for researchers interested in the PS that some popular sensory stimuli such as a single whisker or a pure tone also constitute point stimulations. In many functional imaging reports, there was no intention to quantify the spatial extent of the PS. Consequently, in some reports, the information about the spatial extent of the PS could be gleaned only from the scale bar that accompanied the image of the PS. In other reports, the size of the PS spatial extent could be inferred from the superposition of the PS image on the underlying map of cortical architecture, such as cytochrome oxidase maps. Moreover, even if quantified, there was no standard way to compare the spatial extent of PSs across different sensory cortical areas within the same species or across different species [see Fig. 1(a)]; for example, should comparisons be made on an absolute scale in units of millimeters, or on a relative scale in units of percentage relative to either the entire area of a given primary sensory cortex or perhaps the area of the entire cortical hemisphere? Compounding this issue is the difficulty in determining how similar is the stimulus between the different point stimulations being used because: (i) typically, the description of the point stimulation is less than complete (e.g., a single whisker was stimulated, but what was the angle of deflection?) and (ii) there is ambiguity in identifying the equivalent stimulation parameters across the different sensory modalities (e.g., a particular lumens for visual stimulation versus a particular dB for auditory stimulation versus a particular angular displacement for somatosensory stimulation). Finally, it is not clear what criteria should be used to decide where evoked activity is still significant over spontaneous activity, criteria that can strongly influence the size of the reported spatial extent of the PS. Nevertheless, despite all the above-mentioned issues, a picture emerges regarding the large spatial extent of the PS and its ubiquity in different sensory cortices of different species.

By employing a CCD camera system that permitted the direct (i.e., without a microscope) wide-field imaging of intrinsic signals in the anesthetized rat primary somatosensory cortex through the thinned skull, Masino et al.²¹ achieved the first optical imaging of a PS; more specifically, a PS evoked by stimulating a single whisker. In order to create an image of a $3\times 3$ matrix of neighboring whisker PSs (each stimulated and imaged separately), Masino et al.²¹ needed to apply a high activity threshold of 2% to 5% of peak evoked optical activity of each whisker PS to exclude areas of overlap with the other single whisker PSs [Fig. 2(a), top right panel]. Without thresholding, however, the total spatial extent of a single whisker PS can be seen and its spatial extent is rather large²² [see Fig. 2(a), left panel], $\sim 15{\mathrm{mm}}^2$ and thus, 2 orders of magnitude larger than the $\sim 0.15{\mathrm{mm}}^2$ area for a single whisker barrel²³ [compare left panel and bottom right panel of Fig. 2(a); see also Fig. 1(a), left panel, for approximate size of rat SI]. The initial imaging studies of single whisker PSs strongly suggested that a single whisker PS has an even larger spatial extent than indicated by prior optical recording experiments, a PS that extends far beyond its own barrel and neighboring barrels²² resulting in highly overlapping whisker PSs.²⁴ Since then, as long as a quantification threshold was employed that permitted a more comprehensive inclusion of the spatial extent²⁵ the Frostig lab continued to find a large PS evoked by a single whisker;^10,23,26,27 for mouse, see Ref. 28.

Fig. 2

Cortical PSs in various primary sensory cortices of the anesthetized rat as imaged with ISOI. (a) The large spatial extent of a PS in SI as evoked by a single whisker. Upper right: nine whiskers (B1, B2, B3, C1, C2, C3, Dl, D2, and D3) were individually stimulated and imaged through the thinned skull. The region of strongest activation was determined for each whisker (plotted in hot colors) and compiled into the composite image as shown here. Lower right: Photomontage of the cytochrome oxidase staining of the anatomical representations in layer IV cortex for the same nine whiskers from the same rat that underwent imaging. Note the similarity in the relative size and location of the whisker representations as determined by imaging versus anatomy. From Ref. 21. Left: ISOI data from a single whisker (C2) is chosen and a transformation from a 2-D view of the data (specifically containing only peak activity) to a 3-D view of the data from the entire recorded area is shown here. The white dotted horizontal line is used to illustrate the level of high thresholding needed in order to isolate only the peak activity, as shown in the top right panel. Note that a high thresholding can lead to an incomplete view of the PS by excluding most of the PS (barrel of C2 $\sim 0.15{\mathrm{mm}}^2$, imaged PS $\sim 15{\mathrm{mm}}^2$). Also note the large spatial extent of a single whisker PS relative to the size of the entire SI [see Fig. 1(a), left panel]. Modified from Chen-Bee and Frostig.²² (b) The large spatial extent of a PS for other examples of point stimulation. Left: PS in SI evoked by a somatosensory point stimulus delivered to the fur located on the knuckle of the 4th forepaw digit. The stimulus consisted of five pulses delivered at 5 Hz, each $\text{pulse duration}=10\mathrm{ms}$. Middle: PS in VI evoked by a visual point stimulus delivered to the eye. The stimulus consisted of a white LED point ($0.35\mathrm{deg}\times 0.35\mathrm{deg}$) delivered in five 10-ms pulses at 5 Hz. The “break” in the smooth decline of the left side of the PS originates from a blood vessel artifact. Right: PS in AI evoked by an auditory point stimulus delivered to the ear. The stimulus consisted of five pure tone pips (8 KHz, 87 dB, 100 ms on 100 ms off) delivered at 5 Hz. Note that large PSs were observed in all three primary sensory regions [compare to Fig. 1(a), left panel] and also note their similarity. Chen-Bee and Frostig unpublished data.

Employing ISOI, Das and Gilbert²⁹ were the first to image a large PS (diameter ranging between 3.2 to 5.2 mm) in the anesthetized cat primary visual cortex following a point stimulus (line segment 0.5 deg). Sharon et al.² were the first to employ VSDOI to image the PS from the anesthetized cat visual cortex and they also found a large PS ($\sim 7\mathrm{mm}\times 5\mathrm{mm}$)—following a 4-deg diameter of moving gratings stimulation—peaking over the expected retinotopic cortical area. In the first application of ISOI in the auditory cortex of a rat in response to a pure tone, Bakin et al.³⁰ imaged a large auditory PS ($4\mathrm{mm}\times 2\mathrm{mm}$) following 45-dB pure tone (23 KHz) stimulation [see Fig. 1(a), left panel for approximate size of rat A1], findings which were later replicated by Tsytsarev and Tanaka.³¹ Additionally, large PSs were imaged with ISOI or VSDOI from the auditory cortex of ferrets,³² cats,³³ and guinea pigs.³⁴

The ferret ISOI study also demonstrated that the PS spatial extent could depend on the magnitude of the stimulation: the stronger the stimulation magnitude of the PS the larger the areal extent of the PS; at 70 and 80 dB, large areas covering the entire auditory cortex could be activated, as also shown in the cat auditory cortex³³ and using VSDOI imaging in the guinea pig auditory cortex.³⁴ Areal extent and amplitude of the whisker PS were also progressively dependent on stimulation strength. Petersen et al.³⁵ demonstrated that a very weak (0.6 deg) whisker stimulation evoked a low amplitude PS that lasted for a short duration and spread over a relatively small area of the cortex (just above the appropriate barrel area). Stronger stimulation amplitudes (2 and 6 deg) evoked a progressively stronger amplitude, longer duration, and larger spread, respectively. The 6-deg stimulation—still a relatively weak whisker stimulation for the awake, behaving rodent³⁶—resulted in a large PS over the entire PMBSF;³⁵ see also Berger et al.³⁷ and Deneux and Grinvald³⁸ for similar findings regarding the spatial extent of the PS following different whisker stimulation amplitudes. However, the amplitude of the stimulation is unlikely a major variable that explains the large spatial extent, despite differences in anesthetic agent, stimulus type (e.g., mechanical, air-puffs, magnetic), stimulation parameters and image analysis techniques, similar large-scale (most or the entire PMBSF, or even larger than PMBSF) whisker PSs were reported employing both imaging techniques when the whisker stimulation was at least of a moderate amplitude as can be inferred from rat studies by Blood et al.;³⁹ Jones et al.;⁴⁰ Derdikman et al.;⁴¹ Devor et al.;⁴² Nemoto et al.;⁴³ Civillico and Contreras;^44,45 Lippert et al.;⁴⁶ Devonshire et al.;⁴⁷ Lustig et al.;⁴⁸ and Ollerenshaw et al.;⁴⁹ and in the mouse Prakash et al.;²⁸ Petersen et al.;³⁵ Ferezou et al.;^20,50 Berger et al.;³⁷ Lim et al.;⁵¹ and Mohajerani et al.⁵²

In most studies, the PSs were imaged in or near the center of the cortical area of interest and typically, PSs peaking near borders with other cortical areas were rare. Such preferences could bias interpretations regarding general characteristics of different PSs. Case in point are whisker PSs in the PMBSF. In most studies, single-whisker PSs were recorded or imaged following stimulation of whiskers located at the center of the whisker array on the snout and therefore, due to the somatotopic mapping of these whiskers, at the center of the PMBSF [e.g., whisker C2 barrel shown in Fig. 3(a)], whereas whisker PSs whose barrels are localized near the border of the PMBSF (border whisker PSs) remained rarely imaged. For example, whisker A2 barrel is located at the border of PMBSF, which also constitutes the border of SI, whereas whisker E2 barrel is located at the border of PMBSF and the trunk representation that is part of the body representation within SI [see Fig. 3(a) for locations]. Employing ISOI in the anesthetized rat, Brett-Green et al.²³ compared the PS for the central whisker C2 and the border whiskers A2 and E2, and found that irrespective of their location within PMBSF, the PSs (which in all cases co-registered with the location of their appropriate barrels) for central versus border whiskers shared a similar large spatial extent and smooth decline of evoked activity in all directions away from the peak.

Fig. 3

Electrophysiological confirmation of the large spatial extent of the PS in rat SI. Representative case of single units (SU) and LFP responses evoked by whisker A2 stimulation as recorded using an eight-electrode array arranged in a straight line. (a) Schematic of cortical cytoarchitecture based on flattened layer IV CO-stained brain slice. Barrels A2, C2, and E2 are highlighted in colors. MCx, motor cortex; VCx, visual cortex; ACx, auditory cortex; SII, secondary somatosensory cortex; DZ, dysgranular zone cortex; PV, parietal ventral cortex; IC, insular cortex. (b) Lesions produced by the eight-electrode array. Note electrode #1 lesion is located at the A2 barrel while electrodes 2 to 3 are located in the dysgranular zone, and electrodes 4 to 8 span almost the entire auditory cortex (ACx). (c) Evoked SUs (averaged from 128 trials) decay over cortical distance and completely disappear after electrode 4. (d) Evoked LFPs, recorded from the same microelectrodes, also decay over distance but are still present at the last electrode. (e) Evoked LFPs are abolished in electrodes 7 to 8 (red traces) after lidocaine injection between electrodes 7 and 8, with full recovery after 45 min (f). (g) Control saline injection had no effect. Bottom, magnification of traces boxed in (d)–(f). From Ref. 58.

An important development has been the ability to image the PS in awake animals such as awake monkeys trained to fixate on a point stimulus.^53,54 Imaging V1, Slovin et al.⁵³ found PS dimensions of $9\mathrm{mm}\times 6.5\mathrm{mm}$ following a tiny $0.1\mathrm{deg}\times 0.1\mathrm{deg}$ visual stimulus, and Palmer et al.⁵⁴ found a PS of $6\mathrm{mm}\times 8\mathrm{mm}$. Notably, Palmer et al.⁵⁴ found that for a given stimulation amplitude, the PS was almost constant (invariant) in spatial extent at different (yet limited range of) eccentricities despite changes in magnification and population receptive fields scatter at those eccentricities. This finding suggests a potential for generalization regarding the spatial extent of PSs at different parts of a cortical area as this finding resembles the similarity in PS spatial extent for center versus border whiskers, as described above.²³ Indeed, whisker PSs in the highly magnified PMBSF subarea are similar in their spatial extent to other PSs in nonmagnified subareas of SI [Fig. 2(b), left panel] and in different sensory cortices of the rat [Fig. 2(b), middle and right panels], suggesting further support for similarity of PSs irrespective of differences in cortical mapping characteristics such as magnification.

As with the progress of PS mapping in the awake monkey, technological progress enabled whisker PS imaging and mapping in the awake, behaving rodent. Specifically, using fiber optics and VSDOI, this was accomplished in freely moving mice by the Petersen group.⁵⁰ Notably, these important experiments demonstrated that passive single whisker stimulation evoked a PS that is larger and lasts longer than passive single whisker stimulation in the same mouse when anesthetized, results that were recently reproduced by Mohajerani et al.⁵² Together, these findings demonstrate that the awake animal PS could exhibit an even larger spatial extent and stronger amplitude as compared to the anesthetized cortex.

4.2.

Electrophysiological Studies of the Point Spread

Thus far, PS results have been described for wide-field optical recordings and imaging studies. Compared to the findings obtained from the microelectrode studies used to establish topography as a fundamental principle of cortical organization, these imaging PS results were congruent with the expectation that the peak evoked activity should be registered with a specific topographical location but were disparate in that the large spatial extent of the PS contrasted with the expectation for the spread of activity not to radiate too far away from the peak. Possibly, this disparity was due to the differences in what was being recorded by the different techniques: whereas microelectrode recordings measure supra- and subthreshold activity directly from neurons, VSDOI directly images population-evoked activity mostly originating from the membranes of dendritic trees, and ISOI directly measures hemodynamic-based intrinsic signals as an indirect measure of neuronal activity. Therefore, the implications of the functional imaging results obtained from employing either optical technique would benefit from an explicit confirmation by microelectrode recordings. The optimal approach would be to first image the PS and then to perform microelectrode recordings postimaging within the same animal using the same stimulation regimen, with the imaging results guiding the placement of microelectrodes at different locations inside and outside the imaged PS.

Already during the early study by Grinvald et al.¹⁸ mentioned above, postoptical recordings results were investigated by using a single microelectrode to record sequentially at different locations of the imaged PS in the primary visual cortex of the anesthetized monkey. Using the same stimulation as for their VSDOI imaging, they found that the evoked action potentials comprised only a small area centered with the location of the imaged PS peak activity; in particular, microelectrode recordings located at the periphery of the imaged PS did not detect any evoked action potentials despite a clear presence of evoked VSD activity. These findings introduced the notion that most of the large spread is composed of evoked subthreshold activity. Indeed, microelectrode recordings at the periphery of the optically evoked PS did not detect any evoked spikes despite clear presence of evoked VSD activity. These findings suggested that there exists in primary visual cortex of monkeys a small cortical zone surrounding peak location, where evoked suprathreshold activity can be recorded and a larger zone surrounding that peak zone characterized by evoked subthreshold activity. Further support emerged from subsequent studies in the cat primary visual cortex by Das and Gilbert,²⁹ who followed ISOI imaging of the PS with postimaging microelectrode mapping of spiking neurons and estimated that the spiking PS area represented only 5% of the area of evoked PS, and by Sharon et al.,² who employed single unit recordings of evoked spiking responses post-VSDOI imaging and found spiking activity to be localized only at a small cortical zone around the peak activity. Together, these findings suggested that a large evoked subthreshold area surrounding peak activity location constitutes most of the imaged PS area in the visual cortex of cats and monkeys. These findings were congruent with direct intracellular electrophysiology recordings in the cat primary visual cortex describing a large subthreshold activation area following a small visual stimulation;⁵⁵ although a potential contribution of astrocytes to ISOI imaging in the visual cortex has also been suggested.⁵⁶ In parallel, electrophysiological results in the primary auditory cortex of the rat following pure tone stimulation also demonstrated a large spread of evoked local field potentials (LFPs representing population evoked synaptic activity) in the rat auditory cortex following a pure tone stimulation.⁵⁷

To better understand and further characterize the relationship between these suprathreshold and subthreshold zones of the PS, Frostig et al.⁵⁸ conducted a study in which imaging a single whisker PS with ISOI was followed by detailed neuronal recordings utilizing an array of eight independently controlled microelectrodes organized in a straight line to span across different locations of the PS with increasing distances away from peak activity [for example, see Fig. 3(b)] and enabling recordings at different cortical depths. As in the previous study,²³ both central (C2) and border whiskers (A2 and E2) neuronal PSs [see Fig. 3(a) for locations] were investigated.

The obtained results of Frostig et al. in the rat PMBSF confirmed and extended previous results obtained from the visual cortex, namely by explicitly demonstrating the simultaneous existence of two activity zones comprising the whisker PSs: a suprathreshold and a subthreshold—both of which registered with the peak, with the latter being much larger than and surrounding the former [Figs. 3(c) and 3(d)]. Evoked suprathreshold whisker PS amplitude progressively declined over cortical distance away from peak activity as far as 1.5 mm away in $>50\%$ of rats [see PSTHs in Fig. 3(c)], which would be equivalent to a circular areal extent of $7.1{\mathrm{mm}}^2$ and as far as 2.5 mm maximal distance, which would be equivalent to a circular areal extent of $19.6{\mathrm{mm}}^2$. The evoked subthreshold (LFP) whisker PS also showed a progressive decline of its amplitude over cortical distance away from peak activity [Fig. 3(d)], but unlike the evoked suprathreshold case, evoked subthreshold PS always reached the last microelectrode (3.5 mm away from peak PS, a radial distance equivalent to a circular areal extent of $38.5{\mathrm{mm}}^2$ at the last microelectrode). Notably, the LFP amplitude at the last microelectrode was still, on average, 11% of the peak amplitude, suggesting that the evoked LFP may radiate even farther than 3.5 mm from the peak. These findings on the extent of the subthreshold PS were obtained irrespective of recording array direction away from peak location, its cortical depth, or whether recordings were obtained from central or border whisker PSs. As with previous ISOI results in the PMBSF, both suprathreshold and subthreshold PSs were symmetrically spreading away from peak location, and the only difference among different whisker PSs was the location of their peak. Notably, targeted injections of sodium channel blocker lidocaine into the territory of the last microelectrodes (3.0 to 3.5 mm away from peak) completely flattened evoked LFPs at these distant locations, clearly demonstrating that the large subthreshold PS is not a case of “volume conductance” artifact but rather a case of actual long-range subthreshold neuronal activation [Figs. 3(e)–3(g)]. In addition, thin gray-matter transections between the array’s middle microelectrodes resulted in flattening of the evoked LFP amplitude to noise level for microelectrodes located beyond the transection line, directly demonstrating for the first time that evoked single whisker subthreshold activity spreads horizontally through cortical gray matter (Fig. 4).

Fig. 4

Evoked LFP following single whisker C2 stimulation spreads horizontally within cortical gray matter. An eight-electrode array is used to record LFPs from the anesthetized rat cortex, with electrode #1 coregistering with whisker C2 barrel. Layers 2 to 3 LFP amplitudes declining over cortical distance (denoted on the $x$-axis as electrode number, with 0.5 mm between electrodes) as evoked by whisker C2 stimulation (black line plot) were almost abolished at electrodes 5 to 8 (green line plot) after thin gray matter transection between electrodes 4 and 5 (black arrows), an effect not observed after an earlier transection parallel to the electrode array (purple line plot). Inset: Nissl-based micrograph of the gray matter transection (arrows); $\text{scale bar}=500\mu \mathrm{m}$. From Ref. 58.

4.3.

Point Spread and the Underlying Cortical Cytoarchitecture

Optical imaging and electrophysiological results summarized in the previous two subsections support the notion that an evoked PS can be very large, with an activity amplitude smoothly decaying to far distances in all directions away from the peak. Because the peak for any given PS is registered at the topographically correct cortical location, the findings of a large and relatively symmetrical spatial extent even for PSs whose peaks are located at the border of a sensory cortical region introduced the possibility of a PS to extend outside its associated sensory cortical region and perhaps even into regions of other sensory modalities. In the Frostig et al. study⁵⁸ already reviewed in the previous section, it should be noted that immediately following the single unit and LFP recordings acquired from the microelectrode array, the authors performed localized lesions to mark the location of the recording microelectrodes in order to explicitly compare the spread of activity to the underlying cortical cytoarchitecture as stained with cytochrome oxidase (CO). In doing so, they were able to directly demonstrate that stimulating border whisker A2 resulted in an evoked subthreshold PS that not only spanned the dysgranular area surrounding SI but also, and more surprisingly, continued and spanned the entire auditory cortex [Fig. 3(b)]. Indeed, as seen in [Fig. 3(b)], only the first microelectrode of the array was located within the PMBSF (above the A2 barrel) and, therefore, still within the SI territory, whereas the remaining seven microelectrodes in the array were located outside the SI territory. Notably, microelectrodes 2 to 3 were located in the dysgranular zone and microelectrodes 4 to 8 spanned almost the entire auditory cortex. Furthermore, even stimulating the central whisker C2, whose barrel is located at the center of the PMBSF and, therefore, farther away from the SI border, still evoked a large subthreshold PS that crossed dysgranular areas to invade parts of the visual and auditory cortices. Similarly, stimulating whisker E2 on the opposite side of C2 [refer to locations in Fig. 3(a)] resulted in a subthreshold PS that spanned the trunk representation of SI and reached the border of primary motor cortex MI. These findings suggested that in order to keep a single whisker PS relatively symmetric and similar in spatial extent (see previous two subsections) irrespective of which whisker was stimulated, it was necessary for a given PS to cross cytoarchitectural cortical borders. Furthermore, these findings also suggested that as the difference between whisker PSs is only in the location of their evoked peak, the closer the peak activity location of a whisker PS is to a cytoarchitectonic cortical border, the deeper the whisker PS invades into other cortical areas that are located closest to the border.

Further supporting evidence was obtained in the mouse PMBSF by the Petersen lab employing VSDOI. When the entire mouse hemisphere was imaged in awake, head-fixed mice, two important aspects were revealed: the whisker PS did not stop at the borders of PMBSF, and a second PS was imaged in the motor cortex that was found to be spreading away from the stimulated whisker anatomical motor representation. Specifically, within 40 ms of active whisker stimulation (hitting an object with single whisker), the entire surface of the imaged hemisphere was activated as a result of the merging of both somatosensory and motor PSs as the evoked activity radiated away from their respective peak activity locations;²⁰ see Fig. 5. Clearly, the whisker PS in this case was spreading far beyond borders of their respective PMBSF and motor cortex and thus, trespassing deeply into other cortical areas, although the authors did not identify the different cortical areas activated by these PSs. Similar VSDOI imaging results supporting the spread beyond the border of PMBSF by a whisker PS were obtained by the Murphy lab⁵² and recently reproduced by the Petersen lab.⁵⁹ In addition, similar results about the spread beyond borders can also be seen in the mouse when imaged with genetically encoded indicators of voltage.⁶⁰ Together, these findings demonstrate that there are two additional characteristics shared between imaged PS and subthreshold PS, which lend further support to the notion that subthreshold activity underlies imaging PS: (1) both types of PSs cross borders and (2) both invade into other cortical areas.

Fig. 5

The large PS following active whisker C2 stimulation as imaged with VSDOI. A single-trial example of active touch imaged in an awake mouse. The C2 whisker actively touches an object (upper row image sequence), evoking a PS first in the somatosensory cortex and then in the motor cortex (lower row image sequence). Note after 40 ms activity is detected over the entire hemisphere. From Ref. 20; reprinted from Neuron, 56(5), Ferezou et al., Spatiotemporal dynamics of cortical sensorimotor integration in behaving mice, 907-923, Copyright (2007), with permission from Elsevier.

4.4.

Relationship between Functional and Anatomical Point Spread

The functional results obtained with optical imaging and electrophysiology recordings summarized so far raise an important question: what could be the underlying anatomical system responsible for supporting the large horizontal spread of PSs? The earlier imaging and electrophysiological studies by Grinvald et al.¹⁸ in the monkey visual cortex and Das and Gilbert²⁹ in the cat visual cortex proposed that cortical long-range horizontal projection within gray matter could constitute the underlying system supporting the PS spread of activity. Several studies were conducted to address this possibility in the PMBSF, but the challenge in this case was the need to demonstrate whether such an anatomical system could support the functional findings, i.e., a case of a system of long-range horizontal projections that are capable of crossing cytoarchitectural borders and, therefore, spatially match with the imaged and electrophysiological PSs, although the option of multisynaptic projections could also be relevant.

In the Frostig et al. electrophysiological recording study,⁵⁸ localized (point) injections of the anterograde tract-tracer biotinylated dextran amine (BDA) were deposited into imaged peak activity locations as a means to characterize the underlying anatomy of the PS for a single whisker. Findings showed clear spread of axons or axon collateral segments diffusely projecting to all directions away from the PMBSF injection site. The point injection and the resulting spread of projections constitute the cortical anatomical PS. The density of the anatomical PS labeled axon segments declined over cortical distance, indicating that the anatomical PS is composed of short, medium, and long-range axonal projections. Congruent with results obtained by optical imaging techniques and microelectrode array recordings, such axonal projections also clearly ignored cytoarchitectonic borders by trespassing into the auditory and visual cortices as identified by layer IV flattened cortex cytochrome oxidase staining of the same rat. After crossing the dysgranular area separating PMBSF and visual cortex, traced axonal projection reached 2.5 mm in length spanning into visual cortex. In addition, the authors found dense staining locations of known area-to-area projections through white matter from the PMBSF into expected targets such as secondary somatosensory cortex (SII), dysgranular cortex (ParP, posterior parietal), motor cortex, and insular cortex.

The preliminary anatomical results demonstrating the existence of horizontal long-range, border crossing, diffuse projection system within the gray matter were followed up and extended by a detailed study of this projection system in the PMBSF by Stehberg et al.⁶¹ Using anterograde tract-tracer BDA, injections and detailed mappings of axonal projections at different locations within PMBSF, a clear pattern emerged that supported imaging and electrophysiological PS mapping findings. These findings suggested that irrespective of BDA injection location within PMBSF, two major systems can be clearly described: (1) the traditional area-to-areas system, where axons project through white matter into specific areas (e.g., SII, PVT) and (2) a diffuse, apparently nonspecific projection system within gray matter that projects to all directions and includes very long, border crossing horizontal projections, see Figs. 6(a)–6(c) for a schematic view of these systems.

Fig. 6

Proposed model of anatomical projections for rat PMBSF and comparison between anatomical versus functional PS as imaged with ISOI. (a–c) Schematic of the proposed distinction between a specific versus a diffuse system of anatomical projections for the rat PMBSF. Schematics depict relevant cortical areas (a) to be compared to the proposed diffuse system of long-range border crossing projections (b) and the more familiar, specific system of main outputs projections (c). VI, primary visual cortex; AI, primary auditory cortex; SI, primary somatosensory cortex; SII, secondary somatosensory cortex; dys, dysgranular cortex; PMBSF, posteromedial barrel subfield of SI; PVT, parietoventral cortex. From Ref. 61. (d) and (e) The similarity between functional PS and anatomical PS of a single whisker. (d) Functional PS. Evoked activity spreads horizontally for long distances following single whisker stimulation. The intrinsic signal optical imaging response following stimulation of the C2 whisker (the first 500 ms containing the maximal areal extent of the initial dip activity) was averaged across 37 rats, as described by Chen-Bee et al.⁶⁶ and was plotted as a false-color image of fractional change relative to prestimulus values. The “outer circle” has a diameter of 7 mm and represents an extrapolation of the farthest electrode used in the recordings of Frostig et al.,⁵⁸ at which an evoked field potential response could be detected in 100% of animals. Black outlines show locations of whisker barrels and cytoarchitectonic areas that were detected by cytochrome oxidase staining in a representative animal. Trunk, trunk region of primary somatosensory cortex; VCx, visual cortex; ACx, auditory cortex. (e) Anatomical PS. Average density of projection pattern across all brains and section depths for supragranular injections of AAV-CMV-GFP. After averaging across four section depths, data arrays corresponding to the 7.2-mm diameter analysis region were re-expressed for each of the nine brains injected with AAV-CMV-GFP in supragranular layers. The re-expressed arrays were then averaged across the nine brains to reveal a largely symmetrical axonal radiation in which axonal density declines with distance from the injection site. Scale: white denotes no detectable projections, the transition from green to yellow denotes the mean density across the analyzed region, and red denotes the maximum density across the region. From Ref. 63; Brain Structure and Function, Long, intrinsic horizontal axons radiating through and beyond rat barrel cortex have spatial distributions similar to horizontal spreads of activity evoked by whisker stimulation, 221, 2016, 3617-3639, Johnson and Frostig, © Springer-Verlag Berlin Heidelberg 2015, With permission of Springer.

Johnson and Frostig^62,63 further extended the BDA findings by employing detailed qualitative and quantitative analysis of the anatomical PS within PMBSF using modern tract-tracing techniques. Specifically, discrete deposits of adeno-associated virus (AAV) were injected into the supragranular or infragranular layers of the PMBSF. The AAV vectors used either the direct expression of enhanced green fluorescent protein (GFP) under cytome-galovitus (CMV) or enhanced yellow fluorescent protein under a calcium/calmodulin-dependent protein kinase IIα (CaMKIIα), with the former labeling all types of neurons and the latter labeling only excitatory neurons. The detailed qualitative and quantitative analysis of the anatomical PS for both types of vectors confirmed and further extended previous results of Frostig et al.⁵⁸ and Stehberg et al.⁶¹ Analysis revealed anatomic PS of labeled axons diffusely radiating in all directions for distances $>3.5\mathrm{mm}$ originating both from supragranular and infragranular injections with declining density over cortical distance—patterns that were common to both types of vectors. Detailed reconstruction of single axons originating from each injection site demonstrated how projections radiated away from the injection site and across the PMBSF, branched and sometimes crossed into other sensory cortices as identified by the underlying layer IV cytochrome oxidase staining. The spatial similarity between functional and anatomical PSs is shown in Figs. 6(d) and 6(e). It should be noted that the anatomical mapping results of the diffuse system constitute a conservative estimate of the horizontal projections density, because this mapping was achieved using $20\times$ magnification and while more projections could be seen with $100\times$ magnification, the higher magnification was used only for full-length reconstruction of single axons. Despite the conservative estimate, the anatomical PS was still spatially matching the imaged and electrophysiological PSs by spanning distances in the order of several millimeters.

The combined BDA and AAV findings indicated that at least for a single whisker in the rat PMBSF, the anatomical PS [Fig. 6(e)] shared many characteristics with the imaged PS [Fig. 6(d)] or subthreshold electrophysiological recordings (Fig. 3). These include large and relatively symmetrical spatial extent, ability for border crossing and ability for trespassing into other cortical areas, diffuse and smooth decline over cortical distance. The BDA and AAV experiments can be combined together with the transection (Fig. 4) and lidocaine (Fig. 3) electrophysiology experiments of Frostig et al.⁵⁸ to implicate the diffuse, horizontal system—composed of short, medium, and/or long-range monosynaptic and multiple synaptic projections—as the underlying anatomical PS supporting the functional PS as imaged by optical imaging techniques and confirmed by electrode array recordings. These findings also suggested that a large-scale, structure–function spatial similarity exists involving the functional and anatomical PSs at the mesoscopic level of cortex.

4.5.

Point Spread is a Dynamic Entity Capable of Contextual Modulation and Plasticity

As described above, the results from PS imaging and postimaging neuronal recording in the visual cortex and PMBSF of anesthetized and awake animals suggested that the PS spatial extent is relatively constant (invariant) to its location within a given cortical area (i.e., different eccentricities in visual cortex, different barrels in PMBSF) for a given stimulus amplitude. However, the PS is a dynamic entity that can exhibit contextual modulation and plasticity. The dynamic nature of the PS was first demonstrated in the anesthetized monkey in the Grinvald et al. study using VSDOI.¹⁸ The visual cortex PS amplitude was always attenuated by various patterns of flanking stimulations irrespective of their size and orientation, although the degree of attenuation was dependent on these variables. Similar findings were observed in the awake monkey using ISOI. The PS amplitude was attenuated following the presentation of identical, two flanking PSs.⁶⁴ Similar to the monkey, contextual interactions among PSs were observed in the PMBSF of the anesthetized rat using ISOI. The amplitude of response to simultaneously stimulating an entire array of whiskers was attenuated as compared to what might be expected based on a linear summation of single-whisker PSs associated with each of the whiskers in the array. In general, this amplitude attenuation was dependent on the number of whiskers being stimulated: the more whiskers stimulated, the stronger the response attenuation (see more details below in Sec. 5.3).^65,66 Similar findings were obtained using VSDOI for pairs of stimulated whiskers in the anesthetized mouse and rat.^44,48,67,68 Together, these studies suggested that attenuation by contextual modulation of the PS by flanking stimulations could be a general rule irrespective of cortical area imaged, anesthesia versus awake preparation, and imaging technique, but more data are needed to further substantiate this potential rule.

A whisker PS can exhibit use-dependent plasticity of its spatial extent and amplitude when the rat’s environment changes. Repeated imaging supported by single unit recordings demonstrated that a spared whisker PS contracts its spatial extent and its amplitude diminishes when rats have the opportunity for whisker-dependent scanning of an environment outside their cage.^27,69 Conversely, its PS expands when the same single spared whisker is only used inside its standard cage, and no changes are found for control, nonspared whisker rats. Allowing the rest of the previously deprived whiskers to regrow to their original length results in the return of the single PS back to its baseline values (Fig. 7).⁶⁹ Despite this strong case of plasticity, the basic characteristics of the PS, such as peaking at the expected topographical location and the symmetrical declining away from the peak, remained intact. Moreover, even without deprivation of whiskers—just allowing rats to live for 28 days in a “naturalistic habitat” that promotes tunnel digging, interactions with other rats, and foraging activities—results in PS contraction and reduction in its amplitude as compared to controls that are housed in standard cages residing within the same room.²⁷ While the findings were obtained after 28 days, cellular and molecular studies comparing rats living in a naturalistic habitat to matched standard cage controls in the same room showed that plasticity can start as early as 2 days after transfer to the naturalistic environment,⁷⁰ therefore, raising the possibility that future imaging studies of whisker PS could detect such plasticity at short periods following transfer to the naturalistic habitat. However, extremely long, continuous living in the naturalistic habitat ($\sim 13$ months) showed no difference in PS spatial extent and amplitude compared to standard cage controls within the same room (Kvasnak and Frostig, unpublished results). Taken together, these findings suggest that while the whisker PS can exhibit significant experience dependent plastic changes, there seems to be a spatial extent and amplitude “set point” that it is returning to, implying a possibility for homeostatic control on the spatial extent and amplitude of the PS.

Fig. 7

Imaging PS plasticity using ISOI. Plasticity of a spared whisker’s functional PS is reversible upon restoration of normal sensory input. (a) Schematic of the experimental design. (b) Examples of ISOI images for the spared whisker PS obtained from a sensory-deprived animal that remained in its home cage (top row) or was given an opportunity for spared whisker guided exploration outside its home cage (bottom row) are taken before deprivation, after 28 days of deprivation, and after 28 days of whisker regrowth. Ratio values are converted to grayscale values in which the prestimulus baseline is shown as gray, and the black and white values on the grayscale bar are set to a decrease or increase of $2.5\times {10}^{-4}$ from baseline values, respectively. For detailed analysis of the plasticity at all activity levels of the spared whisker PS, see Ref. 69.

Additional supportive evidence regarding the behaviorally relevant dynamic nature of the whisker PS originated from the studies of the PMBSF in alert or behaving mice demonstrating that the areal extent and amplitude of a whisker PS were dependent on a wakefulness state and the behavioral context. Passive whisker stimulation during quiet wakefulness resulted in a larger, stronger PS during wakefulness as compared to the PS in the same mouse when anesthetized.⁵⁰ Similarly, Mohajerani et al.⁵² replicated these findings and extended them by demonstrating that the large PS imaged in auditory cortex in response to pure tone during quiet awake state has stronger amplitude and lasts longer compared to the PS evoked by an identical stimulus in the same mouse when anesthetized. Furthermore, whisker PS became smaller and exhibited weaker amplitude during active whisking in the air, yet became large and strong again (similar to quiet wakefulness) following an active whisker touch.⁵⁰ Recently, in the Petersen lab, VSDOI imaging was used to image the cortex of mice trained in a simple, goal directed learning task: to lick a water reward spout after a 1-ms deflection of a single whisker. Hit trials, in which the mouse licked after a whisker stimulus, were accompanied by a larger whisker PS compared to miss trials. Prestimulus whisking decreased behavioral performance by increasing the fraction of miss trials, and these miss trials had attenuated cortical responses, a strong example of behavioral modulation of the PS.⁵⁹

5.1.

Point Spread as a Potential Structure–Function Neocortical Motif

The findings reviewed thus far collectively introduce the possibility that a large PS, whether it is functional or anatomical, may be a motif repeated across neocortical gray matter. For example, as reviewed in Sec. 4.1, large functional PSs have been observed for different species (rat, mouse, cat, monkey, ferret, guinea pig) and cortical areas of different sensory modalities (visual, somatosensory, auditory). As reviewed in Sec. 4.4, large anatomical PSs have also been described and quantified specifically for the rat somatosensory cortex, in particular, the PMBSF region of somatosensory cortex. Since then, the Frostig research team have pursued additional experiments to explicitly address whether large functional and anatomical PSs are observed across a myriad of cortical areas including even nonsensory regions.

Both imaging and anatomical experiments were pursued for various sensory areas of the rat cortex. While employing the same technical (e.g., imaging data collection and data processing) and surgical (e.g., anesthesia) protocols as those applied for imaging whisker PSs [Fig. 2(a)], the Frostig research team imaged the functional PSs for various sensory cortical areas in order to address potential variability typically encountered when comparing results obtained from different labs. As seen in [Fig. 2(b)], point tactile stimulation to the skin at the knuckle forepaw area, point visual stimulation to the eye, and point (pure tone) stimulation to the auditory system all resulted in a PS that is similarly large as to the PS obtained in the rat PMBSF, characterized by one peak above the appropriate topographical map location and declining symmetrically over cortical distance. Therefore, experiments were also conducted to determine whether a large anatomical PS occurs for cortical areas other than the rat PMBSF that could explain the similarity of such imaged PSs. The anatomical PSs for various primary sensory cortical areas were mapped by the anterograde tract-tracer BDA.⁶¹ BDA was injected into the supragranular layers of the somatosensory, visual and auditory cortex. Similar to the imaged PSs obtained for these primary sensory areas [Fig. 2(b)], the anatomical PSs were also found to be large and were characterized by a progressive reduction of projections density with cortical distance away from peak density. Furthermore, the anatomical PSs for these other sensory areas shared other features with that for PMBSF: diffuse long-range projections spreading in all directions and spanning into other primary cortices; the spatial extent of PS projections into other unimodal cortices was dependent on the location of the injection site (i.e., the closer to the cytoarchitectonic border the deeper the axonal projections into other cortices); and axonal labeling in specific known targets were also detected presumably connected by white matter projections. When combined with the findings previously described earlier in the review, collectively these imaging and anatomical findings suggest that large functional and anatomical PSs may constitute a motif shared across the various sensory cortices in the rat. Furthermore, they suggest that beyond the traditional area-to-area projections to specific targets through white matter, there is a system of diffuse, long-range horizontal projections through gray matter that could directly connect primary cortices. In other words, the “two-system” schematic of projections [Figs. 6(a)–6c)] that was based on results obtained in the PMBSF could also be applicable for other sensory areas. The potential implications of these findings to cross-modal interactions and plasticity have been discussed in Ref. 61.

Still, other experiments have been pursued to determine whether a large PS can be observed even for nonsensory cortical regions. Employing AAV virus vectors and the same qualitative and quantitative analyses as described above for the PMBSF, Johnson and Frostig⁷¹ have further demonstrated that injections in the supragranular layers of all known types of neocortical tissue including granular (trunk area of somatosensory cortex), dysgranular (several dysgranular areas surrounding PMBSF), and agranular cortex (motor cortex) reveal diffuse anatomical PSs that radiate horizontally through gray matter in all directions and trespass cytoarchitectonic borders into other areas surrounding them. Indeed, the quantified images of the anatomical PS mapping in all these cortical areas were quite similar to the ones obtained in the PMBSF. That the findings obtained for nonsensory cortical areas are similar to those obtained for SI, AI, and VI suggest that large PSs could constitute a unifying motif of organization for the entire neocortex. In other words, a possibility exists for every point in gray matter to be the source (and recipient) of the diffuse horizontal system in a distance-dependent fashion.

The ubiquity of large PSs observed across different species from mice to monkeys and across different cortical areas from primary sensory to dysgranular and agranular cortex suggests the possibility that a large PS may be an evolutionary conserved general motif of the mammalian neocortex. Moreover, the fact that functional and anatomical PSs trespass borders strengthen the evidence described above that the neocortex can be viewed more of a continuum rather than a parceled entity as already suggested by Frostig et al.⁵⁸ in their 2008 study. The ubiquitous presence of large PSs, which require precious metabolic resources for their development and especially their maintenance, begs the question: what is the relevance for a PS to be large? The next two subsections provide additional findings about PSs that point toward possible relevance for their existence.

5.2.

Relative Spatiotemporal Invariance of a Point Spread

As mentioned earlier in Sec. 1, the PS is composed of many thousands of neurons engaged in coordinated activity and thus, can be considered a classic case of Cajal’s “neuronal avalanche” ensemble. A follow-up study of the PS spatiotemporal characteristics provided the opportunity to discover emergent properties not readily attainable simply by studying cortical activity at the single neuron level. As already reviewed in an earlier section, for a given stimulus amplitude, the spatial extent of the PS was similar irrespective of where its peak location resided within a given sensory cortex, and thus, this characteristic of the PS could be thought of as being invariant with respect to its location in the cortex. Findings from the follow-up study suggest that another characteristic of the PS, namely its relative spatiotemporal profile in millisecond resolution, can also be invariant to stimulation magnitude.

Following ISOI of a whisker PS, a 32-microelectrode array (eight microelectrodes organized in a straight line in layer 1, eight in layers 2 to 3, eight in layer 4, and eight in layer 5) was employed to simultaneously record from all microelectrodes both suprathreshold and subthreshold neuronal activity within the PMBSF. Jacobs et al.⁷² characterized the spatiotemporal profile of a single whisker PS in response to a variety of stimulation amplitudes delivered at 5 Hz: four angular displacements ranging from 0.035 to 7.5 deg that differed in a logarithmic (base 6) scale ranging from barely perceptible to the eye (0.035 deg) up to moderate strength (7.5 deg). The acquired data were analyzed based on 1-ms frames of spatiotemporal movies obtained from each microelectrode. Detailed spatiotemporal movies of evoked cortical activity were obtained for each of the stimulation amplitudes. These movies of subthreshold activity (LFP) and suprathreshold activity (multiunit potentials) recorded from the same microelectrodes clearly demonstrated that the PS exhibits a stronger evoked amplitude and larger spatial extent when a stronger stimulation amplitude is delivered, which is congruent with results of previous studies that also employed various amplitudes of sensory stimulation, as previously reviewed in Sec. 4.1. However, when these spatiotemporal movies were normalized to the largest amplitude of evoked response per stimulus strength, the normalized movies demonstrated that the spatiotemporal characteristics of the PS are the same irrespective of stimulation amplitude, despite the major difference between stimulation amplitudes ($>200$ fold difference between the weakest and strongest stimulation amplitude). Figure 8 shows an example of the raw and normalized LFP spatiotemporal movies created for all stimuli strengths for the first deflection of the whisker. Similar results were obtained for consecutive deflections of the 5-Hz stimulation for both LFP and multiunit potentials. In other words, while the absolute spatiotemporal profile of a PS can depend on stimulation amplitude, its relative profile was invariant to major changes in stimulus strength. These findings suggest that the relative spatiotemporal profile can be added alongside the spatial extent size as being another characteristic of the PS that emerges at the ensemble level as being invariant. Such findings may be of interest to those studying underlying physiological mechanisms of sensory coding for perceptual invariance.

Fig. 8

Invariance of the relative spatiotemporal response profile of a whisker PS. (a) $7\times 4$ microelectrode arrays were used to record LFP responses evoked by a single whisker (five deflections delivered at 5 Hz) across a mesoscopic section of SI extending beyond the boundaries of barrel cortex and penetrating through most cortical layers in six rats (left and middle panels). A color scale can be applied to the LFP response magnitudes collected from the $7\times 4$ electrode arrays as a means to provide “snapshots” of mesoscopic activity at a given time point (example provided in right panel). (b) 1-ms frame movies of LFP responses evoked by the first deflection of a single whisker for four stimulus amplitudes (0.035, 0.2, 1.25, and 7.5 deg). Note the laminar and lateral spread of LFP response across time as evoked by a single whisker, and that the amplitude of response is dependent on stimulation amplitude. (c) The relative spatial profile of evoked LFP spread for each stimulus amplitude can be obtained by normalizing each movie to its own maximum value across all pixels and time points. Note that after normalization, the relative spatiotemporal profile is similar across all stimulus amplitudes. From Ref. 72.

5.3.

Point Spread could Serve as a Building Block for Integrated Cortical Activation

While a whisker PS is quite useful as an animal model for cortical function, the stimulation of a single whisker is a rather artificial situation when compared to what actually happens in an awake, behaving rodent. The study of simultaneously stimulating the entire whisker array on the snout is important because it represents a more naturalistic type of stimulation that would be actually encountered by the rat. Due to their large spatial spread, neighboring whisker PSs highly overlap with each other, which begs the question: what would be the integrated cortical response when multiple PSs are simultaneously activated? To address this question, a study was conducted that employed ISOI and microelectrode array recordings, along with modeling and pharmacological manipulations, to characterize the cortical activity spread in response to multiwhisker stimulation.⁶⁶

To formulate predictions, the PS for whisker C2 as averaged across 37 rats was obtained [Fig. 6(d)] and used as a template for any whisker PS, justified by our previously described findings that PSs of different whiskers are similar.²³ A model of the possible cortical response to multiwhisker stimulation (i.e., simultaneous activation of multiple PSs) was generated simply by linearly summating all individual PSs associated with the stimulated whiskers. Specifically, a copy of the PS template was used for each stimulated whisker such that its peak activity coregistered with the location of that whisker’s barrel within PMBSF, and then all appropriately aligned PSs were spatially summated. The generated model for all 24 large facial whiskers on the snout is illustrated in [Fig. 9(a), left]. This model provided predictions about various aspects of the cortical response, such as the peak magnitude, peak location within the PMBSF, and overall 3-D shape of the response magnitude as it decayed away from the peak. Interestingly, despite the stimulation of many whiskers, the model predicted a single peak at the center of the activated area, as opposed to 24 local peaks as possibly could be expected based on the one-to-one principle of organization for the whisker-to-barrel system. Furthermore, the model predicted that the response magnitude would decay relatively symmetrically in all directions away from the peak activity (in a fashion reminiscent of a single whisker PS) as though the activity collectively evoked by stimulating multiple whiskers was being funneled toward a single peak in its center (response to multiwhisker stimulation could also be referred to as a “funneled” response). In other words, the area of maximal overlap between the spatially summated PSs evolves into the peak of the funneled response, and at progressive distances away from the peak, the overlap among PSs decreases resulting in the weakening of the evoked amplitude with cortical distance away from the area with maximal overlap. To generalize this simulation result, a model was generated in a similar manner for a $2\times 2$ matrix, hence 4 whiskers, located on a corner of the snout [Fig. 9(a), right]. Again, the model predicted a single, funneled peak at the center of the evoked area that weakens relatively symmetrically in all directions away from peak activity, but it also predicted that the peak would shift to the appropriate location within the PBMSF—i.e., the location of peak evoked activity at the center of the barrels associated with the $2\times 2$ matrix of whiskers; except for the peak location, however, the four-whiskers model shared many features with the 24-whiskers model.

Fig. 9

Modeling versus in vivo (ISOI imaging and electrophysiological confirmation) results for stimulating all the large whiskers (left column) or the $2\times 2$ corner matrix of whiskers (right column) of the whisker array. (a) The predicted model of imaging cortical response to stimulating 24 or 4 whiskers. Modeling is based on using whisker C2 functional PS averaged across 37 rats as a building block [see Fig. 6(d) for a partial view of the building block when using a 3.5 mm radius from peak activity] by serving as a representative PS for any individual whisker and performing linear spatial summation of the PSs associated with all the individual whiskers intended for stimulation. Provided here are the 2-D (insets; $\text{scale bar}=1\mathrm{mm}$; imaging data superimposed on barrels map; rostral toward the right; lateral toward the bottom) and 3-D images of the predicted response for the same $6\times 6\mathrm{mm}$ cortical area to stimulating either all 24 (left panel) or 4 (right panel; corner whiskers D3, D4, E3, and E4) large whiskers. The same color scale is applied to both the 2-D and 3-D images. Note the differences in the overall size, shape, and location of the predicted response to stimulating 24 (left panel) versus 4 (right panel) whiskers when compared to each other as well as to single whisker C2 [see Figs. 2(b) and 6(d)]. In particular, the modeling predicts that the cortical response to stimulating all 24 whiskers is spatially larger compared to a single whisker PS and has a single activity peak that is centrally located within the barrels map, while the cortical response to stimulating corner whiskers D3, D4, E3, and E4 is predicted to exhibit a spatial extent larger than for a single whisker PS, but smaller than for 24 whiskers, and with a single peak whose location has shifted in the anterior direction away from the center of the barrels map to the center of the $4\times 4$ activation. (b) The in vivo imaging data averaged across rats. Note the similarity of the in vivo data compared to the predicted models. (c) The rostrocaudal line plot for the model versus in vivo cortical response. Magenta shading indicates SE. Note the ability of the model to predict many aspects of the cortical response obtained in vivo, including the overall shape of the activity spread, the single peak of activity, and the location of peak activity (center of $x$-axis co-registers with center of barrel map). The models explain 80% and 79% of the in vivo variability for the 24 and 4 whiskers cases, respectively. (d) LFP electrophysiology confirmation using an eight-electrode array arranged in a straight line within supragranular layers, with the third electrode located centrally within the barrels map. Average response to stimulating either 24 or 4 whiskers (magenta) is compared to the average response to stimulating whisker C2 (black) within the same rats. 0 mm on $x$-axis co-registers with center of barrel map. (e) Same as for (d) except for multiunit recordings from the same electrodes within the same rats. From Ref. 66.

Comparison between the simulation results and the functional imaging (ISOI) results obtained in vivo while simultaneously stimulating multiple whiskers is summarized in Figs. 9(a)–9(c), along with neuronal recordings (multiunits and LFP) confirmation [Figs. 9(d) and 9(e)]. The models succeeded in predicting many aspects of the obtained in vivo funneled responses, including a single peak, the location of that peak within the PMBSF, and the decay in magnitude as the response radiated away from the peak activity. Interestingly, the model grossly overestimated the peak amplitude of the experimental data [Figs. 7(e) and 7(f) in Ref. 66], suggesting that some form of mutual cortical activity attenuation occurs when simultaneously stimulating multiple whiskers. Collectively, these findings suggest that a PS-based rule (specifically, the linear spatial summation of individual PSs that are topographically aligned) can be used to model various aspects of the cortical response to simultaneous stimulation of multiple whiskers.

Results from additional experiments [Figs. (10) and (11)] further extended these findings. When the temporary sodium channel blocker lidocaine was locally injected into one side of the funneled response evoked by stimulating 24 whiskers [Fig. 10(a)], rather than simply silencing the portion of the funneled response receiving the lidocaine [Fig. 10(b)], the resultant response appeared as a new funneled response with its peak shifted away from the site of lidocaine injection [Figs. 10(c)–10(e)]. This result reconfirmed the observation that the evoked funneled peak response of multiwhisker stimulation is always found toward the center of the response area. In addition, the result suggested that such a change in peak location originates within the cortex [Fig. 10(c)], because if the single funneled peak location of multiwhisker activation was only dependent on thalamocortical input, then it would not relocate following the lidocaine injection [Fig. 10(b)]. A different set of experiments was conducted to determine whether results could be replicated for a cortical area other than somatosensory cortex. As illustrated in Fig. 11, modeling based on the spatial summation of individual PSs can also be used to predict various features of the response in visual cortex to simultaneous stimulation of two visual point stimuli. Last, in the previous section, it was already mentioned that using subthreshold (LFPs) and suprathreshold (multiunit potential) neuronal recordings, Jacobs et al.⁷² found the relative spatiotemporal profile of a single whisker PS to be invariant to major changes in whisker stimulation amplitude. The PS-based modeling of the cortical response to multiwhisker stimulation, where each whisker’s relative spatiotemporal profile is invariant to stimulus amplitude, would predict that the relative spatiotemporal profile of the funneled response evoked by multiwhisker stimulation should also be invariant to stimulus strength. Indeed, Jacobs et al.⁷² obtained subthreshold and suprathreshold neuronal recordings results evoked by multiwhisker stimulation that were found to be relatively invariant to stimulation amplitude. These findings further support the suggestion that a single whisker PS is potentially a building block of integrated activity following multipoint stimulations.

Fig. 10

Cortical response evoked by the 24-whiskers array can be predicted by modeling using single whisker PSs as building blocks for activity interactions occurring at the cortical level. (a) Eight electrodes spaced 0.5-mm apart recorded the funneled response for 24 whiskers (middle electrodes aimed at the funneled peak activity location) before and after a local lidocaine injection (green) deposited distal (1.5 mm) to the middle electrodes, at 300- to $450\text{-}\mu \mathrm{m}$ cortical depth. (b) If subcortical activity interactions were the sole contributors to the funneled response, then no shifting in the funneled peak activity location would be expected to occur outside the lidocaine site. (c) If activity interactions occurring at the cortical level (double-headed arrows) contributed to the funneled response, then the local silencing of cortical activity should lead to a new spatial overlap of the remaining unsilenced single whisker PSs that in turn should lead to a shift in the funneled peak location away from the lidocaine site and toward the center of the new activated area. Results from (d) sub- and (e) suprathreshold neuronal recordings initiated before versus a few minutes after a targeted lidocaine injection are congruent with predictions of panel (c). Recordings initiated 1 h after the lidocaine injection revealed a recovery of response almost to preinjection levels (insets). From Ref. 66.

Fig. 11

Separate and simultaneous activation of two functional PSs at different locations in rat visual cortex as imaged with ISOI; model and experimental results. Labeling in left column applies to remaining columns. Rostral is toward the left and lateral is toward the bottom for all panels in top two rows. In the same animal, ISOI was used to image the cortical response evoked by two visual point stimulations, whether delivered individually (50 deg $\text{azimuth}=\text{first column}$; 80 deg $\text{azimuth}=\text{second column}$; both at elevation of $+15\mathrm{deg}$) or simultaneously (third column). The fourth column provides the predicted model based on linear spatial summation of the individual PSs as seen in columns 1 and 2. Upper row: 2-D maps scaled in fractional change units. Middle row: topographical map displaying isoactivity lines. Bottom row: 3-D images scaled in fractional change units. From Chen-Bee and Frostig, unpublished data.

Taken together, the results (Figs. 8–11) collectively suggest that PSs can be used as a building block for the cortical response to simultaneous multipoint stimulation. Further research is needed to test the robustness of such PS-based modeling for other combinations of multipoint stimulation as well as across various cortical areas and mammalian species. If PS-based integration, based on our modeling, is found to be robust, it would be interesting to see whether it may be useful for addressing questions related to integrative evoked cortical activity patterns associated with perception and cognition.