INTRODUCTION

Broadband high-frequency activity (BHA; 70 to 150 Hz), also known as high gamma, a key analytic signal in human intracranial recordings, is often assumed to reflect local neural firing [multiunit activity (MUA)]. Accordingly, BHA has been used to study neuronal population responses reflecting auditory (1, 2), visual (35), language (6), mnemonic (710), and cognitive control (11, 12). BHA is arguably the electrophysiological measure best correlated with the blood oxygenation leveldependent (BOLD) signal in functional magnetic resonance imaging (fMRI) (13, 14). However, beyond the fact that BHA correlates with neuronal spiking (13, 1517), its precise physiology remains unknown; specifically, the neuronal populations and physiological processes generating BHA have not been identified. Here, we show that BHA can be dissociated from MUA in the primary visual and auditory cortex. Analyzing laminar multielectrode data recorded in monkeys, we found a bimodal distribution of stimulus-evoked BHA across the layers of cortex: an early-deep and late-superficial response. The early-deep BHA had a clear spatial (laminar) and temporal association with local MUA; however, the spatiotemporal overlap of MUA with the late-superficial BHA signal was much less clear. In many cases, particularly in V1 (70%), supragranular sites showed a strong BHA in lieu of any detectable increase in MUA. Because of volume conduction, BHAs from both the early-deep and the late-superficial generators contribute to the field potential (FP) at the pial surface, but the contribution is weighted toward the superficial BHA. Although both BHA components volume conduct to the pial surface, the strongest generators of BHA are in the superficial cortical layers. The origins of BHA include a mixture of the neuronal action potential firing and dendritic processes separable from firing. It is likely that, as typically recorded, BHA signals emphasize the latter processes to a greater extent than previously recognized.

We investigated the neural mechanisms generating BHA signal and their relationship to MUA in the neocortex by analyzing FP and concomitant MUA signals recorded with laminar multielectrodes in macaque primary visual (V1: two animals, 104 experimental sessions with whole-screen flashes, 49 sessions with free viewing) and auditory (A1: two animals, 26 experimental sessions with broadband noise, 26 sessions with best frequency tone) cortices. We used first and second derivative (current source density or CSD) analyses of the laminar FP profiles, along with a concomitant MUA to localize the neuronal generators of BHA and to determine their spatiotemporal relationship to neuronal firing. Figure 1 displays representative laminar activity profiles recorded from primary visual (Fig. 1, A, B, and E) and auditory (Fig. 1, C, D, and F) cortices. The notion that BHA directly reflects neuronal firing raises the obvious prediction that BHA and MUA should have the same spatial and temporal distributions across cortical layers. Contrary to this prediction, however, the spatiotemporal profiles of BHA and MUA were reliably different. We observed two temporally and spatially distinct BHA components, including the early-deep BHA localized to the granular and infragranular layers and the late-superficial BHA observed in the supragranular layers (Fig. 1, E and F). Early-deep BHA had clear spatiotemporal overlap with granular and infragranular MUA. However, the spatiotemporal association between late-superficial BHA and MUA was much less clear. Late-superficial BHA often had little or no MUA concomitant, particularly in V1.

Color maps show a CSD superimposed with FPs (left) and a BHA (color map) superimposed with MUA (line plots) profiles (right) in V1 (A) and A1 (C) from a representative session. (A) presents data from V1 recordings during diffuse flash stimulation. (C) shows data from A1 recordings during presentation of broadband noise (100-ms duration). x axes indicate time relative to stimulus onset. y axes represent cortical depth with increasing numbers from superficial to deep layers. Vertical lines indicate sensory events. (B and D) Line plots show the distribution of MUA-BHA cross-correlation coefficients calculated across trials for each penetration of V1 (B) and A1 (D). Box plot presents distribution of the lags obtained from the cross-correlations. Positive values indicate that MUA leads BHA. Shading in the line plot reflects SEM. Central mark and edges in box plots show median and 25th and 75th percentiles. (E and F) show maps of z statistics calculated across trials (n = 137 and 84) for individual V1 (E) and A1 (F) penetration. The statistics values are obtained from a nonparametric pairwise test on MUA (left) and BHA (middle) relative to the prestimulus baseline. Nonsignificant time points are masked (P < 0.05; controlled for multiple comparisons across channels and time with the Benjamini and Yekutieli procedure). Right panels show an overlap (yellow) between the BHA (green) and MUA (red) statistic masks.

Consistent with earlier work (13, 1517), we found that MUA and BHA were correlated at all cortical depths in both V1 and A1 (all Spearmans rho, >0.76; all P < 0.01). We used cross-correlation between BHA and MUA at supragranular, granular, and infragranular layers to test whether the relation between these two signals changes across cortical depth. The cross-correlation coefficient profiles were nonuniform across layers (Fig. 1, B and D). Specifically, the coefficients in the supragranular layers peaked later, relative to those in deeper layers [median lags in V1 supragranular, granular, and infragranular layers were 9, 1, and 2.5 ms; while median lags in A1 supragranular, granular, and infragranular layers were 3.25, 0.5, and 0.25 ms; both Kruskal-Wallis (KW) tests, P < 0.01; Fig. 1, B and D]. This suggests that the relationship between BHA and MUA might vary rather than being constant across cortical depth. To understand this divergence, we quantified the laminar/temporal distributions of both signals, which confirmed the above impressions about the changing association between BHA and MUA across cortical depth.

We tested laminar distributions of BHA and MUA across four different experiments in two cortical areas, V1 and A1. First, we used diffused light flash stimulation (V1-DF), which weakly activates the supragranular layers of V1 relative to patterned stimulation (18); this contrast has been known for decades [see, e.g., (19)]. We selected these nonoptimal stimuli because we expected the dissociation to be strongest under these conditions, where activation of the superficial layers is relatively weak. We observed that BHA and MUA distributions were both nonuniform across V1 cortical depth (both KW tests, P < 0.001, n = 104), yet they differed markedly (Fig. 2, A to D). BHA was strongest in layers where neural firing was sparse to undetected. It was more pronounced in the supragranular than both granular and infragranular layers (all Wilcoxon tests; P < 0.001). MUA displayed the opposite pattern: It was decreased in the supragranular compared to both granular and infragranular layers (both Wilcoxon tests; P < 0.001). BHA magnitude (P = 0.03) unlike MUA (P = 0.14) did differ across the granular and infragranular layers. Seventy percent of all V1 penetrations with diffused light flash stimulation showed a strong supragranular BHA in lieu of any detectable MUA increase. The fact that the BHA and MUA have a differing laminar/temporal distributions accords with the view that these signals reflect different aspects of neural activity.

(A to D) present V1 data from recordings during diffuse flash stimulation (V1-DF; n = 104 experiments in two animals). (E to H) Data from V1 recordings during free-viewing exploration of visual images [V1-FV; n = 49, same two animals as in (A to D)]. Line plots show the time course of BHA (A and E) and MUA (C and G) response across supragranular, granular, and infragranular (red, green, and blue lines) layers. x axes indicate time relative to stimulus (A, C) or fixation onset (E, G). y axes represent signal change from baseline (i.e., normalized BHA/MUA). Box plots present BHA (B and F) and MUA (D and H) distributions averaged across time after stimulus onset (B and D) and across the entire pre- and postfixation epoch (F and H). Supragranular, granular, and infragranular (S, G, and I) layers are plotted as separate box plots. Box plots indicate 25th percentile, median, and 75th percentile; whiskers extend to extreme values not considered outliers, while outliers are marked with crosses. Shading in line plots reflects SEM. Note the consistently different laminar distributions of BHA and MUA in both experiments. Despite different stimuli, levels of firing, and BHA magnitude across experiments, BHA in the supragranular layers is enhanced relative to that in the granular and infragranular layers, while MUA in supragranular is sparse compared to granular and infragranular. Note also that in the free viewing (E to H), there is substantial saccadic modulations of MUA and BHA across all layers, with a clear pattern of MUA suppression around the time of the saccade (perisaccadic suppression) and MUA increase at the end of the saccade (onset of fixation). Arrows in (A and E) indicate early (solid line) and late (dashed line) BHA components.

To test whether the observed dissociation generalizes across types of visual stimuli and task, we examined BHA-MUA dissociation in experimental data collected during free viewing of natural scenes in the same V1 penetrations (V1-FV; see Materials and Methods). It is noteworthy that natural scene viewing produces a much more effective activation of the superficial layers of the primary visual cortex (20). There are multiple differences in BHA and MUA morphology elicited by active (i.e., fixation locked) as compared to passive (stimulus locked) visual input. For example, both BHA and MUA decrease during saccade and rebound after saccade termination (Fig. 2, E to H). Despite these differences in both BHA and MUA morphology, we continue to observe the laminar dissociation between BHA and MUA. Both BHA and MUA show nonuniform distributions across layers (both KW tests, P < 0.001, n = 49), with superficial layers showing strongest BHA (all Wilcoxon tests, z > 5.7, P < 0.001) and weakest MUA (all Wilcoxon test, z > 6.0, P < 0.001).

Next, we sampled from A1 during presentation of broadband noise (A1-BBN; see Materials and Methods) to test whether these effects generalize to other areas of the sensory cortex. Despite differing temporal patterns of BHA signals in V1 and A1, key aspects of laminar dissociation between BHA and MUA generalize to A1. As in V1, event-related BHA and MUA were both observed across all three laminar compartments (Fig. 3, A to D), and both had a nonuniform laminar distribution (both KW tests, P < 0.001, n = 26). Critically, BHA was stronger in the supragranular layers compared to granular and infragranular layers (both Wilcoxon tests, P < 0.001). In contrast, supragranular MUA was weaker than granular and infragranular layers MUA (both Wilcoxon tests, P < 0.001). The poststimulus time interval during which BHA was increased from baseline was more sustained in the supragranular than in the granular and infragranular layers (73.5-, 18.7-, and 20.0-ms median supragranular, granular, infragranular, respectively; both Wilcoxon tests, P < 0.001). In contrast, the poststimulus time interval during which MUA was increased from baseline was shorter in the superficial, as compared to deeper layers (15.5-, 42.5-, and 45.5-ms median supragranular, granular, and infragranular, respectively; both Wilcoxon tests, P < 0.001). Neither BHA nor MUA durations differed between granular and infragranular layers (both Wilcoxon tests, P > 0.84). While the spatiotemporal pattern of BHA and MUA differences in A1 largely parallels with that observed in V1, there are noteworthy A1-V1 differences. BHA in V1 during both passive (stimulus locked) and active (i.e., fixation locked) visual input had two elements (early and late) with confined laminar spread; the early component was limited to deep layers, while the late component was limited to supragranular layers. In A1, we observed larger temporal separation between early and late BHA components (median of 33 ms). As in V1, the late BHA component in A1 is largely limited to the superficial layers. However, unlike the early-deep component in V1, the early BHA component in A1 extends up into the supragranular layers. We quantified this by directly comparing the BHA magnitude in A1 across layers in early (0- to 30-ms poststimulus) and late (31 to 100 ms) time window. Differential laminar distribution was noted for the magnitude of the late (KW test, P < 0.001) but not the early BHA (KW test, P = 0.76). In contrast, MUA was decreased in the supragranular layers compared to granular and infragranular layers during both time windows (all Wilcoxon tests, P < 0.001).

(A to D) presents data from A1 recordings during presentation of broadband noise (A1-BBN; duration, 100 ms; n = 26, two animals). (E to H) shows data from A1 recordings during presentation of best frequency tones [A1-BFT; 100-ms duration; n = 26, same two animals as in (A to D)]. Line plots show the time course of BHA (A and E) and MUA (C and G) response across supragranular, granular, and infragranular (red, green, and blue lines) layers. x axes indicate time relative to stimulus (A, C, E, G). y axes represent signal change from baseline (i.e., normalized BHA/MUA). Box plots present BHA (B and F) and MUA (D and H) distributions averaged across time after stimulus onset. Supragranular, granular, and infragranular (S, G, and I) layers are plotted as separate box plots. Box plots indicate 25th percentile, median, and 75th percentile; whiskers extend to extreme values not considered outliers, while outliers are marked with crosses. Shading in line plots reflects SEM. Note the consistently different laminar distributions of BHA and MUA in both experiments. Despite different stimuli (i.e., broadband noise and best frequency tones), BHA in the supragranular layers is enhanced relative to that in the granular and infragranular layers, while MUA in supragranular is sparse compared to granular and infragranular. Arrows in (A and E) indicate early (solid line) and late (dashed line) BHA components.

To test whether these results further generalize to other auditory stimuli, we analyzed data from the local best frequency tone response for each site within the same penetrations in A1 (A1-BFT; n = 26 sessions, two animals). For each individual penetration, we compared MUA responses across a set of 14 pure tones from 354 to 32 kHz (see Materials and Methods). We defined the best frequency tone as the one that elicited the largest magnitude MUA response in the layer L4 [e.g., (21, 22)]. Subsequently, we aggregated laminar BHA and MUA responses elicited by these best tone frequencies. This analysis reproduced our initial finding of nonuniform laminar distribution of BHA and MUA (KW test, both BHA and MUA, P < 0.01), with BHA being, again, strongest in the superficial layers (all Wilcoxon test, z > 3.3, P < 0.01) and MUA being weakest in the superficial layers (all Wilcoxon test, z > 2.9, P < 0.001; Fig. 3, E to H).

Next, we evaluated the extent to which these two BHA components (deep versus superficial) contribute to the BHA recorded on the pial surface [i.e., electrocorticographic (ECoG)like signal; Fig. 4, A to D]. We reasoned that a signal from the cortical depth that explains most of the variance in the pial surface BHA would indicate a major laminar generator of surface BHA. To quantify the contribution of generators from different depths to the pial surface BHA, we used a standard method of tracing continuous FP from a neural source (generator) to the cortical surface [(22, 23); see Materials and Methods]. On the basis of both the magnitude of supragranular BHA and the proximity of supragranular current generators to the pial surface of cortex, we predicted that the supragranular layers should provide the largest contribution to BHA signal at the cortical surface. To quantify this, we estimated the percentage of variance explained (adjusted R2) by four linear regression models, including the data from either individual layers (Fig. 4) or all layers together (fig. S8). We used a linear regression model because it has been shown that volume conduction itself is linear at a macroscopic scale (2325). For both V1 and A1 experiments, most of the variance was explained by supragranular BHA (see Fig. 4, E to H; KW test, P < 0.05). We also observed that BHA from deep layers could explain some significant portion of variance, particularly in V1 passive viewing. The amount of variance explained by supragranular layers reached 15, 25, 12, and 20% (for V1-DF, V1-FV, A1-BBN, and A1-BFT, respectively), while it peaked at about 5% for granular and infragranular layers in all four experiments. Thus, while the BHA from all laminar compartments contributes to the pial surface BHA, the supragranular layers appear to be the major source.

(A to D) BHA recorded from pial surface, supragranular, granular, and infragranular layers (gray, red, green, and blue, respectively) of (A) V1 during whole-screen flash stimulation (n = 31 experiments). (B) Presentation of broadband noise during A1 recordings (n = 24), (C) free viewing of visual images during V1 recordings (n = 18), and (D) presentation of best frequency pure tones during A1 recordings (n = 24). y axes represent signal change from baseline (i.e., normalized BHA). (E to H) Portions of the amount of variance (adjusted R2) in the pial surface BHA signal explained by each of three different models, containing laminar BHA from either supragranular (red), granular (green), or infragranular (blue) layers as predictors. Colored markers over panels represent cortical depth, which explains most of the variance (KW test across distributions of adjusted R2 values).

To test whether the spatial dissociation between BHA and MUA noted in monkeys could be observed in humans, we used a small dataset from similar laminar probes implanted in a to-be-resected tissue of the prefrontal cortex (PFC) of two patients with pharmacoresistant epilepsy during rest (Fig. 5). In short, spatial distributions of BHA and MUA were notably similar to those observed in monkeys with strongest BHA in the superficial electrodes and MUA in the deeper electrodes. As in our monkey data, we found nonuniform distributions across cortical depth for both BHA and MUA (both KW tests, P < 0.001). There was also a differential distribution of BHA and MUA across cortical depth that was similar to that observed in monkeys: BHA (Fig. 5, A and B) was strongest in the superficial channels as compared to middle and deep channels (both Wilcoxon tests, P < 0.001), whereas MUA (Fig. 5, C and D) was strongest in the middle channels as compared to superficial and deep (both Wilcoxon tests, P < 0.001). This is consistent with our results from nonhuman primates, showing that the main generators of BHA are localized in supragranular layers, whereas the strongest local population neuronal firing is observed in deeper layers. These data were recorded from only two patients during rest (n = 960 and 2532 segments for a duration of 1 s, resulting in a total of 16 and 42.20 min of data) and should be interpreted with caution. However, the data do suggest that BHA and MUA might dissociate in the human neocortex in a manner similar to that we observed in the nonhuman primate neocortex.

Line plots show example traces from simultaneous BHA (A and B) and MUA (C and D) recordings during rest in two patients with epilepsy implanted with linear probes in a to-be-resected tissue. Signals were assigned into three a priori depths (1, 150 to 900 m; 2, 1050 to 1800 m; and 3, 1950 to 2850 m) based on known thickness of the cortex (see Materials and Methods). The whole resting data were segmented into 1-s-long intervals (100-ms-long intervals at side of epoch were removed after filtering to avoid contamination of edge effects). The laminar distribution of BHA and MUA was quantified by testing for nonuniformity of BHA and MUA across depths (see Materials and Methods). Box plots present results aggregated across all BHA (A and B) and MUA (C and D) epochs in each patient separately. Box plots indicate 25th percentile, median, and 75th percentile; whiskers extend to extreme values not considered outliers. Outliers are shown as + signs.

Laminar activity profile recordings in rodents suggest that the most prominent current sink in the supragranular layers, which corresponds to the location generating the late-superficial BHA in our study, reflects the Ca2+-dependent spiking in the apical dendrites of pyramidal neurons (26). The synaptically evoked Ca2+ signal is largely mediated by the N-methyl-d-aspartate (NMDA) receptors (27), which regulate several processes, including neural plasticity (28) and dynamic shifts in neural excitability (29). The effect of NMDA-mediated enhancement of excitability is strongest in the supragranular layers (29), consistent with the preferential expression of NMDA receptors there (30). The laminar distribution of NMDA receptors and the time course of the NMDA-mediated depolarization (slower and primarily supragranular) suggest that the later superficial BHA signal may have a strong NMDA dependence, further dissociating BHA from neuronal firing. To address this possibility, we examined the effects of systemic administration of the noncompetitive NMDA antagonist phencyclidine (PCP) on auditory responses in A1 (n = 8 experiments, one animal). During control recordings (i.e., before PCP), stimulation elicited a sharp and transient increase in both BHA and MUA (Fig. 6). To test the overall effect of PCP on BHA and MUA, we averaged each signal across time within a 150-ms poststimulus window and compared the magnitude of these averaged responses before and after PCP administration (see Fig. 6). After PCP administration, BHA was significantly attenuated across all layers (all Wilcoxon tests, P < 0.01), while MUA showed no detectable difference between control and PCP (all Wilcoxon tests, P > 0.26) at any cortical depth. While there is no significant impact on the MUA, both early and late BHA components are attenuated by the PCP. It is also noteworthy that the early BHA component is attenuated, while the late component is abolished. Albeit based on data from only one subject, the effects of NMDA blockade are robust and provide a pharmacological dissociation between BHA and MUA. Previous studies found enhanced [e.g., (31)] and suppressed (32) BHA after NMDA-antagonist administration. These differential effects might be attributed to the complex influence that NMDA-antagonists have on the stimulus-elicited response (32). Lazarewicz et al. (32) found that an NMDA-antagonist (ketamine) simultaneously enhanced prestimulus and poststimulus gamma activity, while the relative response to stimulus was decreased. The current findings point to a pharmacological dissociation with low doses of an NMDA-antagonist affecting BHA relative response to stimulus while keeping MUA intact.

(A and B) Box plots present BHA (A) and MUA (B) averaged within the 150-ms-long poststimulus time window before (control) and after systemic administration of PCP (n = 8). Box plots indicate 25th percentile, median, and 75th percentile; whiskers extend to extreme values not considered outliers. (C to F) Line plots show the time course of BHA (C and E) recorded from supragranular, granular, and infragranular (red, green, and blue lines) layers of A1 before (C) and after (E) systemic administration of PCP. (D and F) Concomitant MUA recordings from the same supragranular, granular, and infragranular (red, green, and blue lines) layers before (D) and after (F) systemic administration of PCP. x axis indicates time relative to stimulus onset. y axis represents signal change from baseline, averaged over the 50-ms interval before the stimulus (normalized BHA/MUA). Shading reflects SEM.

BHA is a critical mesoscopic signal commonly used to bridge the gap between human electrophysiology and single-unit studies in animals. Moreover, BHA is the closest neural correlate to the BOLD fMRI signal, linking neuroimaging and neurophysiology. On the basis of influential earlier studies, a prevalent view on human electrophysiology is that BHA is a simple reflection of MUA. Across experiments in V1 and A1, we identified spatial, temporal, and pharmacological dissociations between BHA and MUA. Our results suggest that BHA and MUA index different aspects of the neural activity with divergent features across layers of the neocortex. BHA had two spatially and temporally distinct components, i.e., early deep and late superficial, observed in granular-infragranular and supragranular layers, respectively. There is a substantial spatial and temporal correspondence between MUA and the early-deep BHA component. In contrast, the MUA correlate of late-superficial BHA is much weaker and often undetectable. BHA appears to be NMDA mediated, with low doses of noncompetitive NMDA antagonist (PCP) decreasing BHA while leaving MUA intact. Our regression analyses show that, although the BHA from all laminar compartments volume conducts to the pial surface, the main generators of BHA are localized in the supragranular layers where, as shown by numerous prior studies, firing is sparse [(22, 33, 34); for review, see (35)]. This is important because it suggests that the BHA recorded at the pial surface (e.g., with an ECoG electrode) likely overrepresents signals generated in supragranular layers, while the more robust event-related spiking activity is observed in the granular and infragranular layers and often precedes activation, as indexed by BHA. The results of laminar recording experiments in rodents (26) indicate that the large, late-superficial component may be generated by Ca2+ influx during apical dendritic spiking of pyramidal cells. The vulnerability of the largest (late superficial) component of BHA signal to NMDA receptor blockade in our experiments suggests that this component may index the same process.

Several critical implications of our findings for the interpretation of the BHA signal merit further emphasis. First, BHA as recorded from the pial surface of the cortex reflects relatively modest contributions from neuronal firing. Our pharmacological findings (Fig. 6) link to those of Suzuki and Larkum (26) in pointing to an underlying NMDA-mediated process generating BHA. One possibility is that BHA originates from calcium-dependent spikes that are long-lasting (10 to 100 ms) nonsynaptic events triggered by NMDA receptormediated excitatory postsynaptic potentials (36) that have been suggested to be a mechanism for associating information carried by feedforward and feedback pathways (37).

The current finding that the supragranular BHA is the largest contributor to the pial surface signal suggests that the BHA as typically measured in ECoG may contain a substantial representation of input from cortical feedback pathways. The median of 33-ms (interquartile range of 3 ms) onset-to-onset difference between early and late supragranular BHA signals in A1 is much longer than expected for a conduction delay within a direct monosynaptic connection between granular and supragranular layers (35). One possible explanation is that the early BHA reflects feedforward signal propagating to L4 and then to extragranular layers, whereas the late-superficial BHA reflects feedback from higher auditory areas to A1, which provides a strong input to the supragranular layers. In feedback pathways, predictions and contextual information originating in extragranular layers of higher-order areas are projected to and modulate lower cortical areas. Because feedforward and feedback pathways encode different information (33), bias toward feedback circuits in ECoG-derived BHA might favor predictive and contextual information.

The idea that BHA reflects an integrative process separable from the typical action potential may help to interpret several currently unexplained observations. Niessing et al. (14) used recordings from V1 of anesthetized cats to show that BHA is a better correlate of BOLD fMRI than MUA. Both BOLD and BHA encode stimulus intensity at a finer rate than MUA. This is puzzling under the assumption of BHA being a direct consequence of spiking neurons. The current results help to understand this discrepancy by suggesting that BHA (and consequently BOLD fMRI) might reflect dendritic processes subthreshold to neural firing. This interpretation would predict that the correlation between the magnitude of BHA and firing, rather than being constant, depends on the effectiveness of a stimulus in driving the cortex. For stimuli that are more effective in eliciting action potentials, correspondence between BHA and MUA will be stronger than for stimuli that are less effective. In support of this, Nir et al. (15) showed that coupling between firing rates of individual neurons and BHA varies across time rather than being stable. Furthermore, they also showed that the level of spike-BHA coupling depends on the degree of firing rate correlations between neighboring neurons. Similarly, Smith et al. (38) observed changing correlation between BHA and MUA during epileptic seizure. In support of the subthreshold view, Rich and Wallis (17) found that the BHA in the orbitofrontal cortex, although correlated with firing, diverged from MUA on several dimensions, carried more information, and was more sensitive to spatiotemporal changes than neural firing. The current findings suggest that BHA primarily reflects a dendritic process that is separable from neural firing. Our findings also suggest that using BHA as a proxy for neural firing overestimates the onset latency and duration of firing. This outcome introduces important caveats into the interpretation of the BHA signal, but it remains clear that the signal, particularly in the deeper layers, still has a strong relationship to neuronal firing. Improved understanding of the additional neuronal contributions to BHA makes it a richer, more useful index of brain activation. While it is unlikely, it is also possible that BHA and MUA might be generated by two distinct mechanisms operating at the same time at the same cortical depths.

The current findings provide an additional step toward the understanding that in supragranular layers, neural codes could entail dense subthreshold synaptic inputs accompanied by sparse firing (35). It has been clear for some time that the firing in superficial layers is sparse compared to those in the deeper layers. For example, Sakata and Harris (34) studied the population activity in the rat auditory cortex during spontaneous and sensory-evoked stimulation. Both conditions exhibited a sparse, spatially localized activity in the layer 2/3 pyramidal cells, with a densely distributed activity in the larger layer 5 pyramidal cells and putative interneurons. Recordings from pyramidal neurons in primary somatosensory barrel cortex, both under anesthesia and in awake animals, also revealed low spontaneous and evoked firing in layer 2/3 compared to much higher firing rates in layer 5 pyramidal neurons (39, 40). Another study that looked at the barrel cortex in awake head-restrained mice during object localization found superficial layers (layer 2/3) firing was about 20 times weaker than firing rates in deeper layers (41). Similarly, in A1 of nonhuman primates, neuronal firing is typically reduced in the superficial, relative to the middle and deep layers during both spontaneous and stimulus-evoked recordings [e.g., (21, 42)]. A relative paucity of firing in superficial layers is observed in primary somatosensory cortex during stimulation (43) and in resting spontaneous activity (44). Similarly, Bastos et al. (45) observed a weaker superficial firing, as compared to layer 4 in macaque prefrontal cortex (PFC) (see the Supplementary Materials). Self et al. (46) also show weaker supragranular MUA in response to full-screen, high-contrast checkerboard stimulation and during spontaneous recordings in V1. A detailed review of possible mechanisms leading to sparser firing in the superficial layers relative to deeper layers is beyond the scope of our manuscript. However, the fact that layer 2/3 pyramidal neurons appear to require substantially more excitatory synaptic input to drive them to action potential threshold compared to L5 pyramids (35) points to one mechanistic explanation for this common observation.

Acknowledgments: Funding: B.R.A.I.N., MH111439, EYE24776, DC015780, and MH109429. Author contributions: M.L. and C.E.S. designed the study. A.B., Y.K., A.Y.F., and I.U. collected data. M.L. performed analyses. M.L. and C.E.S. wrote the manuscript. All authors contributed to the discussion and interpretation of findings and edited the manuscript. Competing interests: The authors declare that they have no competing interests. Data and materials availability: All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials. Additional data related to this paper may be requested to the authors.

Continue reading here:
Dissociation of broadband high-frequency activity and neuronal firing in the neocortex - Science Advances