Context-dependent relationships between locus coeruleus firing patterns and coordinated neural activity in the anterior cingulate cortex
Abstract
Ascending neuromodulatory projections from the locus coeruleus (LC) affect cortical neural networks via the release of norepinephrine (NE). However, the exact nature of these neuromodulatory effects on neural activity patterns in vivo is not well understood. Here, we show that in awake monkeys, LC activation is associated with changes in coordinated activity patterns in the anterior cingulate cortex (ACC). These relationships, which are largely independent of changes in firing rates of individual ACC neurons, depend on the type of LC activation: ACC pairwise correlations tend to be reduced when ongoing (baseline) LC activity increases but enhanced when external events evoke transient LC responses. Both relationships covary with pupil changes that reflect LC activation and arousal. These results suggest that modulations of information processing that reflect changes in coordinated activity patterns in cortical networks can result partly from ongoing, context-dependent, arousal-related changes in activation of the LC-NE system.
Editor's evaluation
This is a timely and important study that systematically assesses the relationships between neuronal activity in the locus coeruleus (LC) and the anterior cingulate cortex (ACC) in non-human primates. The LC is a major source of cortical norepinephrine that has reciprocal connectivity with the ACC, and the authors have convincingly shown that LC spiking is associated with changes in ACC spike correlations. Further, these changes have consistent phase relationships with pupil size. This is a rare data set that is technically challenging to acquire, and the results are an important advance toward understanding a circuit that is likely to play a role in regulating brain states such as arousal or attention.
https://doi.org/10.7554/eLife.63490.sa0Introduction
Changes in brain state are associated with different levels of arousal, attention, motivation, surprise, and other factors that can affect the activity patterns of large populations of cortical neurons (McAdams and Maunsell, 1999; Purcell et al., 2012; Chang et al., 2012; Falkner et al., 2013; Downer et al., 2015; Ecker et al., 2016; Thiele et al., 2016). These changes are thought to result, in part, from the widespread release of neuromodulators (Aoki et al., 1987; Devauges and Sara, 1990; Schultz et al., 1993; Dalley et al., 2001; Bouret and Sara, 2005; Aston-Jones and Cohen, 2005; McGaughy et al., 2008; Salamone et al., 2009; Pinto et al., 2013; Varazzani et al., 2015; Khani and Rainer, 2016; Minces et al., 2017). Different neuromodulatory systems have different anatomical and physiological properties (Kupfermann, 1979; Aoki et al., 1987; Xiang et al., 1998; Flores-Hernandez et al., 2000; Fernández-Pastor and Meana, 2002; Ding and Perkel, 2002; Salgado et al., 2011; Herrero et al., 2013; Sugihara et al., 2016; Doyle and Meeks, 2017) that are thought to support their different roles in neural information processing (Abbott and Dayan, 1999; Nirenberg and Latham, 2003; Averbeck et al., 2006; Silver, 2010; Beck et al., 2011; Moreno-Bote et al., 2014; Kanitscheider et al., 2015; Lin et al., 2015; Rodenkirch et al., 2019). One prominent example is the locus coeruleus (LC)-norepinephrine (NE) system, whose diffuse projections throughout the brain, close relationship to the sleep-wake cycle, and relationship to electroencephalography (EEG) and pupillometry have led to several theories of its role in arousal-related modulations of cortical activity and function (Aston-Jones and Cohen, 2005; Aston-Jones and Bloom, 1981; Nieuwenhuis et al., 2011; Gilzenrat et al., 2010; Einhäuser et al., 2010; Murphy et al., 2011; Murphy et al., 2014; Joshi et al., 2016; Joshi and Gold, 2020). However, only a small number of studies have shown direct relationships between (artificial) LC activity and cortical activity in vivo (e.g., Devilbiss and Waterhouse, 2011), limiting our understanding of the exact nature of these relationships.
Our aim was to test if and how endogenous, ongoing activity and sensory-driven, evoked responses in the LC relate to changes in neural activity patterns in the anterior cingulate cortex (ACC) of the primate brain. We targeted the ACC because it has strong reciprocal connectivity, both structural and functional, with the LC (Lewis et al., 1979; Morrison et al., 1979; Porrino and Goldman-Rakic, 1982; Jones and Olpe, 1984; Fernández-Pastor et al., 2005; Gompf et al., 2010; Chandler and Waterhouse, 2012; Chandler et al., 2013; Tervo et al., 2014; Köhler et al., 2016; De Gee et al., 2017; Koga et al., 2020). Moreover, ACC neural activity can encode computations that underlie adaptive, goal-directed behaviors, including those that are also associated with indirect measures of LC-linked arousal, such as changes in pupil size and the P300 component of the event-related potential (ERP; Shima and Tanji, 1998; Gehring and Willoughby, 2002; Holroyd and Coles, 2002; Critchley et al., 2005; Matsumoto et al., 2007; Hayden et al., 2011; Shenhav et al., 2013; Ebitz and Platt, 2015; Sarafyazd and Jazayeri, 2019). Thus, interactions between LC and ACC neural activity patterns are likely to have broad behavioral relevance.
We recorded neural activity simultaneously in LC and ACC (Figure 1) and measured the pupil size of alert monkeys under two behavioral conditions: (1) performing a fixation task and (2) performing a fixation task with randomly presented sounds. Both of these conditions have been associated with variations in arousal that covary with pupil size and can affect cognition and behavior (Aston-Jones and Bloom, 1981; Nieuwenhuis et al., 2011; Einhäuser et al., 2010; Gilzenrat et al., 2010; Murphy et al., 2011; Murphy et al., 2014; Varazzani et al., 2015; Joshi et al., 2016). Conditioned on the firing (or lack of firing) of LC neurons, we measured activity patterns of individual neurons and coordinated activity between pairs of neurons in ACC, both of which govern the information processing capacities of neural networks (Britten et al., 1992; Parker and Newsome, 1998; Cohen and Kohn, 2011; Moreno-Bote et al., 2014; Lin et al., 2015; Kohn et al., 2016). We were particularly interested in the timescales over which these features of neural activity in LC and ACC are related, which can provide insights into putative underlying mechanisms. For example, neuromodulatory effects of the LC-NE system on ACC might be expected to have a relatively long time course compared with the time course of typical synaptic events that are mediated by faster, glutamatergic neurotransmission (Feldman RS et al., 1997; McCormick and Prince, 1988; Wang and McCormick, 1993; Schmidt et al., 2013; Timmons et al., 2004). Below we show that such relatively long timescale relationships are evident between LC and ACC neural activity patterns, in particular changes in coordinated activity between pairs of ACC neurons that tend to decrease or increase in relation to spontaneous or evoked LC spikes, respectively.
Results
We analyzed neural activity from simultaneously measured sites in the LC (n = 84 single units, including 78 simultaneously recorded pairs from 35 sites in monkey Sp; 84/80/31 in monkey Ci) and ACC (372/4336/35 in monkey Sp; 275/2875/31 in monkey Ci) and from LC-only recordings (71/12/73 in monkey Ci and 36/8/34 in monkey Oz) while the monkeys maintained fixation on a visual spot. Our analyses focused on whether and how mean spiking activity and individual and pairwise neuronal variability measured in one brain region (ACC or LC) related to neuronal activity in the other brain region (LC or ACC, respectively). To understand the temporal dynamics of these potential relationships, we systematically tested for effects across a broad range of time windows. These windows ranged in duration from 100 ms, which is often used to study neural synchrony, to 1 s, which is consistent with typical timescales of slow neuromodulatory influences (Feldman RS et al., 1997; McCormick and Prince, 1988; Wang and McCormick, 1993; Schmidt et al., 2013; Timmons et al., 2004; Doiron et al., 2016). Because LC and ACC are connected reciprocally, we tested for relationships in both directions. Reliable relationships in either direction are included in the main figures, and the remaining results are shown in figure supplements.
Single-neuron activity during passive fixation
During passive fixation, LC neurons were weakly active, with no spikes measured on more than a third of all trials and otherwise a median (interquartile range [IQR]) firing rate of 1.8 [0.9–3.6] sp/s measured during 1.1 s of stable fixation (Figure 2A). For each session, we divided trials into six groups based on the given LC neuron’s firing rate: trials with firing rate = 0 (LCzero), 1, 2, 3, and ≥4 sp/s formed five groups, and all trials with firing rate >0 formed a sixth group (LCnon-zero). We then assessed ACC neural activity measured at the same time, conditioned on the LC group.
The mean and variance of ACC spike counts, along with their Fano factor (the ratio of variance over mean), increased steadily as a function of the duration of the counting window. These trends showed an apparent, small relationship with simultaneously measured LC spiking activity, such that the mean, variance, and Fano factor measured in the ACC were, on average, slightly higher for LCzero versus LCnon-zero trials (Figure 2B–D). However, these relationships were not statistically reliable when considering data from both monkeys combined together and only in one case (Figure 2G) when considering data separately for each monkey (sign-rank test for H0: median difference between LCzero and LCnon-zero conditions computed per ACC unit, or ANOVA test for an effect of LC firing rate group, p > 0.05 in all cases; Figure 2B–G). Across ACC units, differences in Fano factor for LCzero versus LCnon-zero trials were associated strongly with differences in spike-count variance and more weakly with differences in spike counts (Figure 2—figure supplement 1).
In addition, the mean, variance, and Fano factor of LC spike counts also increased steadily as a function of the duration of the counting window. However, none of these trends showed a statistically significant relationship to simultaneously measured ACC spiking activity (sign-rank test and ANOVA, p > 0.05 in all cases; Figure 2—figure supplement 2).
Coordinated activity during passive fixation
During passive fixation, pairs of neurons in ACC had coordinated spiking activity that varied considerably in sign and magnitude across pairs but also, as has been reported previously, depended systematically on the size of the time window used to count spikes (de la Rocha et al., 2007). Specifically, pairwise spike-count correlations (rsc) in ACC had median [IQR] values of 0.03 [–0.01 0.07] using 200 ms counting windows and systematically increased and became more variable across pairs with larger counting windows up to 1 s (0.05 [–0.03 0.14]; ANOVA test for effect of window size, p = 0.0243).
These ACC spike-count correlations also depended on the spiking activity of LC neurons measured at the same time. Specifically, the correlations were largest when (1) they were measured using relatively large time windows; (2) for ACC pairs that exhibited relatively large, positive correlations independent of LC firing; and (3) when the simultaneously recorded LC neuron was not active (Figures 3 and 4). For the example ACC pair illustrated in Figure 3A–C, rsc was near zero when computed using relatively small time windows and then increased steadily with increasing bin sizes. The magnitude of these increases in rsc with bin size was larger on trials in which the simultaneously measured LC neuron was not spiking versus spiking (Figure 3C).
The population of ACC pairs exhibited similar trends, with rsc values tending to be smaller on LCnon-zero versus LCzero trials, particularly for larger time windows and for ACC pairs with non-negative LC-independent rsc values. To visualize these effects, we first divided ACC pairs into terciles, according to their LC-independent rsc values (Figure 4). Under these conditions, ACC rsc values from the upper two terciles of LC-independent rsc values were reduced by up to ~50% when computed on trials with at least one or more LC spike versus trials with no LC spikes (Figure 4D–F). The magnitude of these reductions was not related in a consistent manner to the magnitude of LCnon-zero activity (ANOVA test for group effect of LC firing rates > 0). These relationships were not evident on shuffled trials, supporting the idea that the trial-by-trial relationships that we identified were not spurious reflections of the spiking statistics of each region considered separately (Figure 4—figure supplement 1). Furthermore, these LC-linked differences in ACC rsc values did not result simply from changes in ACC neuron firing rates because the two measures were unrelated (Figure 4—figure supplement 2A). We found a weak but statistically reliable relationship between LC-linked changes in ACC rsc and changes in ACC Fano factor (Figure 4—figure supplement 2B).
In contrast, spike-count correlations (rsc) from simultaneously measured pairs of LC neurons did not depend systematically on concurrently measured ACC activity. In general, ACC-independent LC rsc values were roughly similar to LC-independent ACC rsc values, increasing steadily and becoming more variable as a function of the counting-window size; for example, median [IQR] values were 0.02 [0 0.09] for 200 ms windows and 0.06 [–0.02 0.19] for 1 s windows. However, these values did not show reliable differences when compared on trials with high versus low ACC firing, analyzed in the same way as the LC-linked ACC rsc values (Figure 4—figure supplement 3). Moreover, the distributions of these ACC-linked LC rsc values did not appear to come from the same (shifted) distribution as the LC-linked ACC rsc values for all five time bins (Kolmogorov–Smirnov test for H0: both sets of values come from the same distribution, p < 0.0189 in all five cases). Thus, differences in single-unit LC activity were associated with differences in coordinated activity patterns of pairs of neurons in ACC, but differences in single-unit ACC activity were not associated with differences in coordinated activity patterns of pairs of neurons in LC.
Relationship between pupil diameter, LC spiking, and ACC coordinated activity during passive fixation
To further examine relationships between coordinated activity in ACC and activation of LC-linked arousal systems, we examined both relative to changes in pupil diameter. During passive fixation, there is a quasi-periodic fluctuation of the pupil (Lowenstein and Loewenfeld, 1969; Pong and Fuchs, 2000; Joshi et al., 2016). Consistent with our previous findings (Joshi et al., 2016), the timing of LC spiking was related to the phase of these ongoing pupil fluctuations, such that LC spiking tended to be higher preceding dilation versus constriction. These modulations showed a systematic precession with dilation phase that corresponded to a delay of ~270 ms from the maximum modulation of LC activity to the relevant pupil change (Figure 5A; median of intercepts from linear regression of per-time-bin peaks of individual LC PETHs; bootstrapped 95% confidence interval from regression for each LC neuron = [177–353 ms]; for examples of pupil cycles and measurement epochs, see Figure 5—figure supplements 1 and 2, respectively).
Combined with our current findings that pairwise correlations in ACC tended to be smaller when LC is active during passive fixation, these results predicted that pairwise correlations in ACC should also vary systematically with the phase of ongoing fluctuations in pupil diameter. We found such a relationship between pupil phase and ACC rsc values. Specifically, over a wide range of pupil phase, decreases in ACC rsc followed increases in LC spiking (the fitted line in panel B is shifted to the right relative to the fitted line in A; note also that increases in ACC rsc, seen as the bright yellow band in the lower-left corner of panel B, occurred even earlier and are examined in more detail in Figure 9). These modulations showed a systematic progression with pupil phase over roughly the same time frame as LC firing, implying that modulations of both LC firing and ACC rsc had relatively fixed temporal relationships to pupil fluctuations and therefore to each other, as well (Figure 5B; Wilcoxon rank-sum test for H0: equal LC PETH and ACC rsc slopes, p = 0.23).
Single-neuron activity in response to startling events
We further examined relationships between LC and ACC neural activity in the context of external events that can cause a startle response. We showed previously that a brief, loud sound played on randomly selected fixation trials (‘beep trials’) can elicit a transient pupil dilation as well as transient responses from individual LC neurons (Joshi et al., 2016; evoked responses from an example LC neuron are shown in Figure 6A, and the LC population average from three monkeys is shown in Figure 6B). The beep stimulus also elicited a consistent, albeit weaker, response in the ACC (evoked responses from an example ACC neuron are shown in Figure 6D, and the ACC population average from two monkeys is shown in Figure 6E). The averaged ACC response included a slight excitation and then inhibition relative to baseline within 200 ms of stimulus presentation. Although we found some individual cases with a biphasic response (e.g., Figure 6D), the average response likely resulted from the range of heterogeneous responses of individual ACC neurons (Figure 6—figure supplement 1).
The startle events also caused transient reductions in the trial-to-trial variability (measured as the Fano factor) of spiking activity of individual neurons in both LC and ACC, with a slightly more sustained effect in the ACC (Figure 6C and F). This reduction in neuronal variability following the onset of a stimulus, also referred to as quenching, has been reported previously for neurons in a range of other cortical regions and for a range of task conditions (Churchland et al., 2010).
Coordinated activity in response to startling events
In the ACC, we found no evidence that the startling sound caused systematic changes in coordinated activity when measured independently of the LC response. Specifically, spike-count correlations were highly variable across ACC neuron pairs, with no statistically reliable, systematic differences when compared before versus after the beep (comparing ACC rsc before versus after beep for spikes measured in 1 s bins; Mann–Whitney U-test, p > 0.05 for both monkeys).
In contrast, we found systematic relationships between the beep-induced changes in coordinated activity in the ACC and the simultaneously measured LC response. Specifically, LC neurons did not show a characteristic transient response on every presentation of a startle stimulus: the median fraction of trials with a response that included a transient increase and then decrease in spiking activity was 54% for monkey Sp and 43% for monkey Ci (Figure 7—figure supplement 1). Therefore, we tested if and how arousal-linked changes in ACC correlated activity depended on whether or not the simultaneously recorded LC neuron responded transiently to the startling sound. Specifically, we computed ACC rsc in time bins of different sizes, separately for 1 s preceding and 1 s following the sound. ACC rsc tended to be larger after versus before the beep stimulus, but only for the subset of trials in which the beep also elicited the characteristic LC response (Figure 6A and B), for both monkeys (Figure 7A). These modulations did not reflect simply transient differences in LC firing rates but rather specific differences between the presence and absence of the stimulus-evoked LC response (Figure 7B). These modulations of ACC pairwise correlations also did not reflect LC-driven differences in the stimulus-driven quenching of response variability in ACC, which was similar on trials with and without LC-evoked responses (Figure 7—figure supplement 2).
Relationship between pupil diameter, LC spiking, and ACC coordinated activity in response to startling events
We showed previously that the magnitude of beep-evoked LC responses is correlated positively with the size of the simultaneously evoked change in pupil diameter (Joshi et al., 2016). We reproduced those results in our current data set (Figure 8—figure supplement 1). We extended those results by showing that these pupil-LC relationships also relate systematically to beep-evoked changes in ACC rsc. Specifically, we measured beep-evoked changes in ACC rsc separately for trials grouped by LC and pupil responses, focusing on how large or small pupil responses related to the evoked-LC versus no-evoked-LC differences in ACC rsc shown in Figure 7. In general, the beep-evoked changes in ACC rsc that required an LC-evoked response were larger on trials in which the evoked LC response was accompanied by a large versus small pupil dilation (Figure 8A). These effects were not evident when considering baseline, not evoked, pupil size (Figure 8B).
Relative timing of neural activity patterns in LC and ACC
To better understand the relative timing of firing patterns (both individual neuron firing rate changes and paired-neuron rsc changes) in LC and ACC, we measured their time courses relative to two external events with well-defined timing: onset of stable fixation on no-beep trials and onset of the sound on beep trials (Figure 9). For no-beep trials, LC firing rates increased after the fixation point turned on and before stable fixation was acquired (Figure 9A). This epoch also included elevated ACC rsc (particularly when the simultaneously recorded LC neuron was not active, consistent with the results shown in Figure 4) that preceded the elevation in LC firing (Figure 9B, gray lines). Both LC firing and ACC rsc tended to stabilize after stable fixation was attained. These results are consistent with the possibility that LC activation prior to fixation onset drives (or at least coincides with) a lasting decrease in ACC rsc during the remainder of the trial.
Following a beep, there was a different pattern of relationships between LC firing and ACC rsc. Specifically, we divided trials that did (Figure 9C and D) or did not (Figure 9E and F) have a beep-evoked LC response (like in Figures 7 and 8). In general, the beep caused a transient increase in ACC rsc, and this increase tended to occur more rapidly on trials with (Figure 9D, gold shaded area) versus without (Figure 9F, gold shaded area) a concomitant LC response (Mann–Whitney U-test, for H0: the linear rate of increase of ACC rsc in this period was larger when LC neuron did versus did not respond, p = 0.0297). Note that the LC transient response (Figure 9C) peaked at the beginning of the epoch in which the ACC rsc started to increase, which is consistent with the idea that the LC-evoked response drives (or at least coincides with) a faster rate of ACC rsc increase (steeper slope in the gold region in Figure 9D compared with F). The evoked LC response was also associated with a more reliable and larger increase in ACC rsc later in the trial (compare Figure 9D and F; Mann–Whitney U-test for H0: ACC rsc 500–1000 ms after the beep was the same on trials with versus without a beep-evoked LC response, p = 0.0350 for both monkeys’ data pooled together) but was not associated with differences in ACC rsc in the time just before the beep (i.e., ACC rsc did not anticipate whether or not the beep would evoke an LC response; Mann–Whitney U-test for H0: ACC rsc –500–0 ms before the beep was the same on trials with versus without a beep-evoked LC response, p = 0.0693 for both monkeys’ data pooled together).
Discussion
We measured relationships between pupil-linked activation of the LC-NE system and changes in cortical neural activity patterns. Previously, we showed that changes in pupil size during passive, near fixation and when driven by external startling events covary with the timing of spiking activity in LC and parts of brainstem and cortex, including the ACC (Joshi et al., 2016). Other studies have shown that pupil size and other indirect measures of LC activation can also correspond to changes in coordinated activity in sensory cortex under certain conditions (Reimer et al., 2014; Reimer et al., 2016; Vinck et al., 2015; McGinley et al., 2015). These previous studies did not explore the relationships between ongoing versus evoked LC activation (or NE release) and concurrent changes in both coordinated cortical activity and pupil size. We extended those findings, including one study that elicited changes in coordinated activity in cortex via manipulation of LC activity (Devilbiss and Waterhouse, 2011), by showing that LC neuronal activity is reliably associated with changes in correlated spiking activity (spike-count correlations, or rsc) between pairs of neurons in ACC. We further showed that these LC-linked changes in ACC rsc depend on the nature of the LC activation. During passive fixation and in the absence of controlled external inputs, ongoing LC firing was associated with a reduction in ACC rsc. In contrast, startling sounds drove transient increases in LC firing that were associated with an increase in ACC rsc. Under both conditions, LC and pupil-linked changes in ACC rsc were most pronounced over relatively long time windows (>500 ms) that are consistent with neuromodulatory timescales and thus might involve the effects of LC-mediated NE release in the ACC (Feldman RS et al., 1997; McCormick and Prince, 1988; Wang and McCormick, 1993; Schmidt et al., 2013; Timmons et al., 2004).
The context dependence of these effects implies that spontaneous versus sensory-evoked activation of the LC-NE system can have different effects on coordinated patterns of activity in cortical networks. These effects seem likely to involve the different LC firing patterns under these conditions, which for sensory-evoked activity include not just a large, transient increase in spike rate but also a reduction in spike-count variability and heterogenous relationships to pupil-linked changes in arousal (Joshi et al., 2016). In principle, different temporal patterns of NE release could mediate a range of synaptic effects via distinct NE receptor subtypes with different affinities and different spatial distributions in different parts of the brain (Aoki et al., 1987; Nicholas et al., 1993; Arnsten et al., 1998; Arnsten, 2000; Berridge et al., 2012; Berridge and Spencer, 2016). These synaptic effects can influence both excitatory and inhibitory neurons as well as astrocytes, resulting in changes in network dynamics in thalamus, cortex, and elsewhere (Segal and Bloom, 1976; Dillier et al., 1978; Waterhouse et al., 1980; McCormick and Prince, 1988; Wang and McCormick, 1993; McLean and Waterhouse, 1994; Fernández-Pastor and Meana, 2002; Salgado et al., 2011; Salgado et al., 2016; Paukert et al., 2014; Guan et al., 2015; Schiemann et al., 2015; Sherpa et al., 2016; Garcia-Junco-Clemente et al., 2019; Aston-Jones and Cohen, 2005; Ohshima et al., 2017; Rodenkirch et al., 2019). Exactly how these diverse mechanisms support context-dependent changes in coordinated activity patterns merits further study, perhaps by selectively controlling the magnitude and timing of NE release in cortex while assessing changes in coordinated activity and using biophysically inspired models to test hypotheses about the underlying mechanisms (de la Rocha et al., 2007; Doiron et al., 2016).
In contrast to the relationships that we identified between LC activation and coordinated ACC activity patterns, we did not find similarly reliable relationships between the spiking activity of individual neurons measured simultaneously in each of the two brain regions. For example, we found slight increases and decreases in mean firing rates and a marked reduction in the trial-to-trial variability of firing rates of individual units in the ACC in response to external stimuli, as has been reported previously for other cortical regions (Churchland et al., 2010). However, these changes in single-unit ACC responses were not related reliably to properties of the concurrently measured LC response that we could measure and thus may involve mechanisms other than the LC-NE system. These findings do not appear to be consistent with previous work showing improvements in signal-to-noise ratios in cortex in response to LC-NE activation or during states of higher versus lower global arousal (Kolta et al., 1987; McLean and Waterhouse, 1994; Lee et al., 2018; Lombardo et al., 2018). This difference could reflect differences in the specific neurons that we targeted in ACC, the different task conditions we tested relative to previous studies, analyses based on ongoing versus evoked LC activity, or some combination of these factors, and merits further study.
What might be the function of LC- and arousal-linked changes in ACC correlations? In general, correlations in spiking activity can be useful or detrimental depending on a number of factors, many of which remain unexplored, particularly outside of sensory cortex (Averbeck et al., 2006; Cohen and Kohn, 2011; Kohn et al., 2016). For example, an increase in correlations is considered a possible mechanism for connecting neural populations over a range of spatial and temporal scales (Gray et al., 1989; Singer, 1999; Riehle et al., 1997). Some computational schemes can benefit from increases in correlated variability (Abbott and Dayan, 1999; Singer and Gray, 1995; Gray, 1999; Kohn et al., 2016; Valente et al., 2021). However, an increase in correlations can also negatively impact the information coding capacity of a large neural population, particularly over longer integration windows (Zohary et al., 1994; Bair et al., 2001; Averbeck et al., 2006; Renart et al., 2010; but also see Nirenberg and Latham, 2003 Moreno-Bote et al., 2014, for alternative interpretations). Global states such as arousal and attention (traditionally linked with LC-NE and cholinergic systems, respectively) can modulate cortical correlations (Cohen and Maunsell, 2009; Mitchell et al., 2009; Herrero et al., 2013; Schmidt et al., 2013). Some work has also suggested that correlations induced by common inputs must be actively decorrelated by the action of local recurrent excitation and inhibition to preserve information fidelity (Ecker et al., 2010; Renart et al., 2010).
Our results suggest that both increases and decreases in coordinated neural activity patterns in cortex may be under context-dependent, moment-by-moment control of the LC-NE system. Event-driven transient activation of LC could help to synchronize populations of cortical neurons by shifting them to a more correlated state. Conversely, during nonevoked, ongoing (event-independent) firing, the NE release could enhance information processing and signal-to-noise ratios by reducing cortical correlations. Further work is needed to identify if and how NE-mediated changes in network dynamics subserve these functions, particularly in the context of goal-directed behaviors that involve multiple brain regions, including LC and ACC, as well as other neuromodulator systems (Hayden et al., 2011; Varazzani et al., 2015; Ebitz and Platt, 2015; Alamia et al., 2019).
Materials and methods
Three adult male rhesus monkeys (Macaca mulatta) were used for this study (monkeys Oz, Ci, and Sp). All training, surgery, and experimental procedures were performed in accordance with the NIH’s Guide for the Care of Use of Laboratory Animals and were approved by the University of Pennsylvania Institutional Animal Care and Use Committee (protocol 806027). The behavioral task and pupillometry recording and analysis techniques were identical to those we used previously (Joshi et al., 2016). Briefly, fixation trials were of variable length (1–5 s, uniformly distributed). The monkey was rewarded with a drop of water or diluted Kool-Aid for maintaining fixation until the end of the trial. On a subset of randomly chosen trials (~25%), after 1–1.5 s of fixation a sound (1 kHz, 0.5 s) was played over a speaker in the experimental booth (‘beep trials’). The monkey was required to maintain fixation through the presentation of the sound, until the fixation point was turned off.
Electrophysiology
Request a detailed protocolEach monkey was implanted with a recording cylinder that provided access to LC+ (the LC and adjacent, NE-containing subcoeruleus nucleus; Sharma et al., 2010; Paxinos et al., 2008; Kalwani et al., 2014), inferior colliculus (IC), and superior colliculus (SC). The detailed methodology for targeting and surgically implanting the recording cylinder and then targeting, identifying, and confirming recording sites in these three brain regions is described in detail elsewhere (Kalwani et al., 2014; Joshi et al., 2016). Briefly, the LC was targeted initially using custom (Kalwani et al., 2009) and/or commercial software (Brainsight). Tracks were then refined using electrophysiological recordings and microstimulation in brain regions dorsal to LC. Neurons in the intermediate layers of SC (SCi) exhibited spatial tuning on a visually guided saccade task and could elicit saccades via electrical microstimulation (Robinson, 1972; Sparks and Nelson, 1987). IC units exhibited clear responses to auditory stimuli. Activity in the trigeminal mesencephalic tract (me5), located immediately dorsal to the LC, showed distinct activity related to orofacial movements such as sipping. LC+ units had relatively long action potential waveforms, were sensitive to arousing external stimuli (e.g., door knocking), and decreased firing when the monkey was drowsy (e.g., eyelids drooped; Aston-Jones et al., 1994; Bouret and Sara, 2004; Bouret and Richmond, 2009). Nonoptimal tracks also helped with mapping; for example, a more medial track could miss me5 but lead to the trochlear decussation with characteristic ramp-and-hold activity related to downward saccades. Likewise, tracks that encountered IC but not SCi were likely too lateral and often missed me5 and LC. Sites were verified using MRI and assessing the effects of systemic injection of clonidine on LC+ responses in monkeys Oz and Ci and by histology with electrolytic lesions and electrode-tract reconstruction in monkey Oz (Kalwani et al., 2014). Recording and microstimulation in these brainstem-targeting tracks were conducted using custom-made electrodes (made from quartz-coated platinum-tungsten stock wire from Thomas Recording) and a Multichannel Acquisition Processor (Plexon, Inc).
ACC cylinders were placed at Horsley–Clarke coordinates 33 mm anterior-posterior (AP), 8 mm lateral (L), in the left hemisphere for monkey Sp and in the right hemisphere for monkey Ci. For ACC recordings, we targeted the dorsal bank of the anterior cingulate sulcus, ~4–6 mm below the cortical surface. ACC tracks were planned and refined using MRI and Brainsight software, as well as by listening for characteristic patterns of white and gray matter during recordings in an initial series of mapping experiments. Recordings were conducted using either custom-made single electrodes or multicontact linear electrode arrays (8- and 16-channel V-probe, Plexon).
For each brain region, we recorded and analyzed data from all stable, well-isolated units that we encountered. Neural recordings were filtered between 100 Hz and 8 kHz for spikes (Plexon MAP). Spikes were sorted offline (Plexon offline sorter). Electrical microstimulation in SCi consisted of biphasic (negative-positive) pulses, 0.3 ms long, 100 ms in duration, and delivered at 300 Hz via a Grass S-88 stimulator through a pair of constant-current stimulus isolation units (Grass SIU6) that were linked together to generate the biphasic pulse.
Data analysis
Request a detailed protocolFor each recorded neuron, we considered spiking activity only during stable fixation, defined as a 1.1 s window that began 1 s after attaining fixation in which the monkey’s gaze remained within a square window 0.2° per side centered on the fixation point. To assess ACC spiking activity patterns conditioned on LC activation, we divided trials within each session based on whether the single LC neuron from which recordings were being made either did (LCnon-zero) or did not (LCzero) produce at least one action potential during stable fixation. We also divided session trials into four additional groups of trials in which the LC neuron fired 1, 2, 3, and ≥4 spikes.
We used 10 bin sizes ranging from 100 ms to 1 s, spaced logarithmically, to count spikes. The mean spike count, variance of spike counts, and the Fano factor (the ratio of the variance to the mean of spike counts) for each neuron were calculated across trials for each bin size. Pairwise spike-count correlations (rsc) were calculated for each bin size as follows. First, trial spike counts from each neuron in the pair were z-scored, and trials on which the response of either neuron was >3 standard deviations different from its mean were removed to avoid effects of outlier responses (Kohn and Smith, 2005; Smith and Kohn, 2008). Then, the MATLAB function corrcoef was used to obtain the Pearson correlation coefficient. Shuffled estimates were made by calculating rsc from pairs of neurons with spike-count vectors generated from randomly selected trials for each neuron.
ACC rsc values were calculated for all trials for each pair of well-isolated neurons, independently of spiking in LC and also separately for trials divided into groups depending on spiking in LC (LCzero and the five LC > 0 groups). The LC-independent rsc values collected from all pairs and measured using ACC spikes counted using each time bin were divided into terciles that corresponded broadly to pairs that were negatively correlated (tercile 1), uncorrelated (tercile 2), or positively correlated (tercile 3). This analysis allowed us to assess whether changes in rsc in one region associated with spiking in the other region depended on the pairs being correlated to begin with or not.
Data Availability Statement
Request a detailed protocolData and Matlab code for all figures in this manuscript are available at here copy archived at swh:1:rev:22b8c94380c057201fb1fcc56c6037b70c56021c; Joshi and Gold, 2021.
Data availability
Data and Matlab code for all figures in this manuscript are available at: https://github.com/TheGoldLab/LC_ACC_paper_Joshi_Gold_2021 (copy archived at swh:1:rev:22b8c94380c057201fb1fcc56c6037b70c56021c).
References
-
Pupil-Linked Arousal Responds to Unconscious SurprisalThe Journal of Neuroscience 39:5369–5376.https://doi.org/10.1523/JNEUROSCI.3010-18.2019
-
Norepinephrine-containing locus coeruleus neurons in behaving rats exhibit pronounced responses to non-noxious environmental stimuliThe Journal of Neuroscience 1:887–900.
-
Locus coeruleus neurons in monkey are selectively activated by attended cues in a vigilance taskThe Journal of Neuroscience 14:4467–4480.https://doi.org/10.1523/JNEUROSCI.14-07-04467.1994
-
An integrative theory of locus coeruleus-norepinephrine function: adaptive gain and optimal performanceAnnual Review of Neuroscience 28:403–450.https://doi.org/10.1146/annurev.neuro.28.061604.135709
-
Neural correlations, population coding and computationNature Reviews. Neuroscience 7:358–366.https://doi.org/10.1038/nrn1888
-
Correlated firing in macaque visual area MT: time scales and relationship to behaviorThe Journal of Neuroscience 21:1676–1697.https://doi.org/10.1523/JNEUROSCI.21-05-01676.2001
-
Reward expectation, orientation of attention and locus coeruleus-medial frontal cortex interplay during learningThe European Journal of Neuroscience 20:791–802.https://doi.org/10.1111/j.1460-9568.2004.03526.x
-
Network reset: a simplified overarching theory of locus coeruleus noradrenaline functionTrends in Neurosciences 28:574–582.https://doi.org/10.1016/j.tins.2005.09.002
-
Relation of locus coeruleus neurons in monkeys to Pavlovian and operant behaviorsJournal of Neurophysiology 101:898–911.https://doi.org/10.1152/jn.91048.2008
-
The analysis of visual motion: a comparison of neuronal and psychophysical performanceThe Journal of Neuroscience 12:4745–4765.https://doi.org/10.1523/JNEUROSCI.12-12-04745.1992
-
Stimulus onset quenches neural variability: a widespread cortical phenomenonNature Neuroscience 13:369–378.https://doi.org/10.1038/nn.2501
-
Attention improves performance primarily by reducing interneuronal correlationsNature Neuroscience 12:1594–1600.https://doi.org/10.1038/nn.2439
-
Measuring and interpreting neuronal correlationsNature Neuroscience 14:811–819.https://doi.org/10.1038/nn.2842
-
Distinct changes in cortical acetylcholine and noradrenaline efflux during contingent and noncontingent performance of a visual attentional taskThe Journal of Neuroscience 21:4908–4914.
-
Activation of the noradrenergic system facilitates an attentional shift in the ratBehavioural Brain Research 39:19–28.https://doi.org/10.1016/0166-4328(90)90118-x
-
Dopamine modulates excitability of spiny neurons in the avian basal gangliaThe Journal of Neuroscience 22:5210–5218.https://doi.org/10.1523/JNEUROSCI.22-12-05210.2002
-
The mechanics of state-dependent neural correlationsNature Neuroscience 19:383–393.https://doi.org/10.1038/nn.4242
-
Task engagement selectively modulates neural correlations in primary auditory cortexThe Journal of Neuroscience 35:7565–7574.https://doi.org/10.1523/JNEUROSCI.4094-14.2015
-
Heterogeneous effects of norepinephrine on spontaneous and stimulus-driven activity in the male accessory olfactory bulbJournal of Neurophysiology 117:1342–1351.https://doi.org/10.1152/jn.00871.2016
-
On the Structure of Neuronal Population Activity under Fluctuations in Attentional StateThe Journal of Neuroscience 36:1775–1789.https://doi.org/10.1523/JNEUROSCI.2044-15.2016
-
Pupil dilation betrays the timing of decisionsFrontiers in Human Neuroscience 4:18.https://doi.org/10.3389/fnhum.2010.00018
-
Spatial representation and cognitive modulation of response variability in the lateral intraparietal area priority mapThe Journal of Neuroscience 33:16117–16130.https://doi.org/10.1523/JNEUROSCI.5269-12.2013
-
In vivo tonic modulation of the noradrenaline release in the rat cortex by locus coeruleus somatodendritic alpha(2)-adrenoceptorsEuropean Journal of Pharmacology 442:225–229.https://doi.org/10.1016/s0014-2999(02)01543-1
-
Pupil diameter tracks changes in control state predicted by the adaptive gain theory of locus coeruleus functionCognitive, Affective & Behavioral Neuroscience 10:252–269.https://doi.org/10.3758/CABN.10.2.252
-
Locus ceruleus and anterior cingulate cortex sustain wakefulness in a novel environmentThe Journal of Neuroscience 30:14543–14551.https://doi.org/10.1523/JNEUROSCI.3037-10.2010
-
dependence and differential modulation by norepinephrineJournal of Neurophysiology 113:2014–2032.https://doi.org/10.1152/jn.00524.2014
-
Pupil Size as a Window on Neural Substrates of CognitionTrends in Cognitive Sciences 24:466–480.https://doi.org/10.1016/j.tics.2020.03.005
-
A method for localizing microelectrode trajectories in the macaque brain using MRIJournal of Neuroscience Methods 176:104–111.https://doi.org/10.1016/j.jneumeth.2008.08.034
-
Measuring Fisher information accurately in correlated neural populationsPLOS Computational Biology 11:e1004218.https://doi.org/10.1371/journal.pcbi.1004218
-
Neural and neurochemical basis of reinforcement-guided decision makingJournal of Neurophysiology 116:724–741.https://doi.org/10.1152/jn.01113.2015
-
Stimulus dependence of neuronal correlation in primary visual cortex of the macaqueThe Journal of Neuroscience 25:3661–3673.https://doi.org/10.1523/JNEUROSCI.5106-04.2005
-
Correlations and Neuronal Population InformationAnnual Review of Neuroscience 39:237–256.https://doi.org/10.1146/annurev-neuro-070815-013851
-
Modulatory actions of neurotransmittersAnnual Review of Neuroscience 2:447–465.https://doi.org/10.1146/annurev.ne.02.030179.002311
-
Norepinephrine enhances song responsiveness and encoding in the auditory forebrain of male zebra finchesJournal of Neurophysiology 119:209–220.https://doi.org/10.1152/jn.00251.2017
-
State dependence of stimulus-induced variability tuning in macaque MTPLOS Computational Biology 14:e1006527.https://doi.org/10.1371/journal.pcbi.1006527
-
Medial prefrontal cell activity signaling prediction errors of action valuesNature Neuroscience 10:647–656.https://doi.org/10.1038/nn1890
-
Effects of attention on orientation-tuning functions of single neurons in macaque cortical area V4The Journal of Neuroscience 19:431–441.https://doi.org/10.1523/JNEUROSCI.19-01-00431.1999
-
Noradrenergic modulation of firing pattern in guinea pig and cat thalamic neurons, in vitroJournal of Neurophysiology 59:978–996.https://doi.org/10.1152/jn.1988.59.3.978
-
Pupil diameter covaries with BOLD activity in human locus coeruleusHuman Brain Mapping 35:4140–4154.https://doi.org/10.1002/hbm.22466
-
Distributions of mRNAs for alpha-2 adrenergic receptor subtypes in rat brain: an in situ hybridization studyThe Journal of Comparative Neurology 328:575–594.https://doi.org/10.1002/cne.903280409
-
Sense and the single neuron: probing the physiology of perceptionAnnual Review of Neuroscience 21:227–277.https://doi.org/10.1146/annurev.neuro.21.1.227
-
BookThe Rhesus Monkey Brain in Stereotaxic CoordinatesSan Diego, California, USA: Academic Press.
-
Fast modulation of visual perception by basal forebrain cholinergic neuronsNature Neuroscience 16:1857–1863.https://doi.org/10.1038/nn.3552
-
Characteristics of the pupillary light reflex in the macaque monkey: metricsJournal of Neurophysiology 84:953–963.https://doi.org/10.1152/jn.2000.84.2.953
-
Brainstem innervation of prefrontal and anterior cingulate cortex in the rhesus monkey revealed by retrograde transport of HRPThe Journal of Comparative Neurology 205:63–76.https://doi.org/10.1002/cne.902050107
-
Eye movements evoked by collicular stimulation in the alert monkeyVision Research 12:1795–1808.https://doi.org/10.1016/0042-6989(72)90070-3
-
Dopamine, behavioral economics, and effortFrontiers in Behavioral Neuroscience 3:13.https://doi.org/10.3389/neuro.08.013.2009
-
Comparative anatomy of the locus coeruleus in humans and nonhuman primatesThe Journal of Comparative Neurology 518:963–971.https://doi.org/10.1002/cne.22249
-
Visual feature integration and the temporal correlation hypothesisAnnual Review of Neuroscience 18:555–586.https://doi.org/10.1146/annurev.ne.18.030195.003011
-
Spatial and temporal scales of neuronal correlation in primary visual cortexThe Journal of Neuroscience 28:12591–12603.https://doi.org/10.1523/JNEUROSCI.2929-08.2008
-
Sensory and motor maps in the mammalian superior colliculusTrends in Neurosciences 10:312–317.https://doi.org/10.1016/0166-2236(87)90085-3
-
Cell-specific modulation of plasticity and cortical state by cholinergic inputs to the visual cortexJournal of Physiology, Paris 110:37–43.https://doi.org/10.1016/j.jphysparis.2016.11.004
-
Attention Induced Gain Stabilization in Broad and Narrow-Spiking Cells in the Frontal Eye-Field of Macaque MonkeysThe Journal of Neuroscience 36:7601–7612.https://doi.org/10.1523/JNEUROSCI.0872-16.2016
Article and author information
Author details
Funding
National Institutes of Health (R21 MH107001)
- Joshua I Gold
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank Long Ding, Adrian Radillo, Alice Dallstream, Kyra Schapiro, David Kleinfeld, and Takahiro Doi for valuable comments, Rishi Kalwani for piloting the LC/fixation studies, and Jean Zweigle for expert animal care and training. This study was funded by R21 MH107001.
Ethics
All animal training, surgery and experimental procedures were performed in accordance with the NIH's Guide for the Care and Use of Laboratory Animals and were approved by the University of Pennsylvania Institutional Animal Care and Use Committee (protocol 806027).
Copyright
© 2022, Joshi and Gold
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 2,333
- views
-
- 357
- downloads
-
- 37
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Medicine
- Neuroscience
The advent of midazolam holds profound implications for modern clinical practice. The hypnotic and sedative effects of midazolam afford it broad clinical applicability. However, the specific mechanisms underlying the modulation of altered consciousness by midazolam remain elusive. Herein, using pharmacology, optogenetics, chemogenetics, fiber photometry, and gene knockdown, this in vivo research revealed the role of locus coeruleus (LC)-ventrolateral preoptic nucleus noradrenergic neural circuit in regulating midazolam-induced altered consciousness. This effect was mediated by α1 adrenergic receptors. Moreover, gamma-aminobutyric acid receptor type A (GABAA-R) represents a mechanistically crucial binding site in the LC for midazolam. These findings will provide novel insights into the neural circuit mechanisms underlying the recovery of consciousness after midazolam administration and will help guide the timing of clinical dosing and propose effective intervention targets for timely recovery from midazolam-induced loss of consciousness.
-
- Neuroscience
Gamma oscillations in brain activity (30–150 Hz) have been studied for over 80 years. Although in the past three decades significant progress has been made to try to understand their functional role, a definitive answer regarding their causal implication in perception, cognition, and behavior still lies ahead of us. Here, we first review the basic neural mechanisms that give rise to gamma oscillations and then focus on two main pillars of exploration. The first pillar examines the major theories regarding their functional role in information processing in the brain, also highlighting critical viewpoints. The second pillar reviews a novel research direction that proposes a therapeutic role for gamma oscillations, namely the gamma entrainment using sensory stimulation (GENUS). We extensively discuss both the positive findings and the issues regarding reproducibility of GENUS. Going beyond the functional and therapeutic role of gamma, we propose a third pillar of exploration, where gamma, generated endogenously by cortical circuits, is essential for maintenance of healthy circuit function. We propose that four classes of interneurons, namely those expressing parvalbumin (PV), vasointestinal peptide (VIP), somatostatin (SST), and nitric oxide synthase (NOS) take advantage of endogenous gamma to perform active vasomotor control that maintains homeostasis in the neuronal tissue. According to this hypothesis, which we call GAMER (GAmma MEdiated ciRcuit maintenance), gamma oscillations act as a ‘servicing’ rhythm that enables efficient translation of neural activity into vascular responses that are essential for optimal neurometabolic processes. GAMER is an extension of GENUS, where endogenous rather than entrained gamma plays a fundamental role. Finally, we propose several critical experiments to test the GAMER hypothesis.