Acetylcholine modulates the precision of prediction error in the auditory cortex
eLife assessment
The findings of this study are valuable as they provide new insights into the role of acetylcholine in modulating sensory processing in the auditory cortex. This paper reports a systematic measurement of cell activity in the auditory cortex before and after the microiontophoretic application of Ach during an oddball and cascade sequence of auditory stimuli. The evidence presented is convincing, as the study used a rigorous experimental design and statistical analysis. The manuscript will interest researchers in auditory neuroscience and neuromodulation, as well as clinicians and individuals with auditory processing disorders.
https://doi.org/10.7554/eLife.91475.3.sa0Valuable: Findings that have theoretical or practical implications for a subfield
- Landmark
- Fundamental
- Important
- Valuable
- Useful
Convincing: Appropriate and validated methodology in line with current state-of-the-art
- Exceptional
- Compelling
- Convincing
- Solid
- Incomplete
- Inadequate
During the peer-review process the editor and reviewers write an eLife Assessment that summarises the significance of the findings reported in the article (on a scale ranging from landmark to useful) and the strength of the evidence (on a scale ranging from exceptional to inadequate). Learn more about eLife Assessments
Abstract
A fundamental property of sensory systems is their ability to detect novel stimuli in the ambient environment. The auditory brain contains neurons that decrease their response to repetitive sounds but increase their firing rate to novel or deviant stimuli; the difference between both responses is known as stimulus-specific adaptation or neuronal mismatch (nMM). Here, we tested the effect of microiontophoretic applications of ACh on the neuronal responses in the auditory cortex (AC) of anesthetized rats during an auditory oddball paradigm, including cascade controls. Results indicate that ACh modulates the nMM, affecting prediction error responses but not repetition suppression, and this effect is manifested predominantly in infragranular cortical layers. The differential effect of ACh on responses to standards, relative to deviants (in terms of averages and variances), was consistent with the representational sharpening that accompanies an increase in the precision of prediction errors. These findings suggest that ACh plays an important role in modulating prediction error signaling in the AC and gating the access of these signals to higher cognitive levels.
Introduction
Neuromodulation strongly impacts sensory processing by influencing neuronal excitability or synaptic processes in neuronal circuits (Hasselmo and Giocomo, 2006; Lucas-Meunier et al., 2003; Metherate, 2011; Muñoz and Rudy, 2014; Picciotto et al., 2012). Acetylcholine (ACh) is a widely distributed neuromodulator throughout the brain, including the AC, and modulates processes such as attention, learning, memory, arousal, sleep, and/or cognitive reinforcement (Batista-Brito et al., 2018; Dalley et al., 2004; Franklin and Frank, 2015). The main source of ACh to the AC is the basal forebrain (Chavez and Zaborszky, 2017; Mesulam, 2013; Zaborszky et al., 2008). In the auditory system, cholinergic modulation is known to alter frequency response areas generating changes across frequency tuning, decreasing the acoustic threshold at the characteristic frequency, and changing the encoding of the spectral representation of many auditory neurons (Ma and Suga, 2005; Metherate, 2011). Thus, ACh promotes neuronal and synaptic plasticity at different temporal scales (Kamke et al., 2005; Kilgard and Merzenich, 1998).
Here, we characterized the effect of ACh on stimulus-specific adaptation (SSA), a type of neuronal adaptation found in the AC. SSA is elicited by an oddball paradigm, which comprises a pattern of repeating sounds (standards), interrupted by a low-probability and unexpected sound (deviant). The deviant usually differs in frequency, but could differ on any other physical dimension from the standard sound, or otherwise violate a pattern of regularity established by the standard (Chen et al., 2015; Nieto-Diego and Malmierca, 2016; Parras et al., 2017; Ulanovsky et al., 2003; von der Behrens et al., 2009).
Thus, neurons that exhibit SSA adapt specifically to the standard stimulus but resume their firing when a deviant stimulus appears. SSA has been proposed to be a neuronal correlate of ‘mismatch negativity (MMN)’, an evoked potential obtained in human and animal electroencephalographic studies using the oddball paradigm (Nieto-Diego and Malmierca, 2016; Ulanovsky et al., 2003). Hence, SSA has been also referred to as nMM in animal models (Casado-Román et al., 2020; Lao-Rodríguez et al., 2023; Parras et al., 2020; Parras et al., 2017; Malmierca et al., 2019; Carbajal and Malmierca, 2020), a nomenclature that highlights deviance detection, as well as other underlying processes other than neuronal adaptation. Because nMM is considered a form of short-term plasticity (Ogawa and Oka, 2015) and ACh has been shown to play a role in this type of neural plasticity (Marshall et al., 2016; Moran et al., 2013; Parr and Friston, 2017; Vossel et al., 2012) it is plausible that ACh may be involved in the generation and/or modulation of nMM. Furthermore, it has been shown that ACh differentially modulates the neural response to the standard stimulus in neurons of the IC (Ayala and Malmierca, 2015).
Currently, MMN and nMM are best explained by the predictive coding theory (Friston and Sporns, 2008; Friston, 2005). According to the predictive coding framework, higher-level cortical areas generate predictions about sensory streams that are conveyed in a top-down manner to lower hierarchical levels, to suppress the ascending neuronal activity evoked by sensory events that can be anticipated. However, when current predictions do not match the sensory inputs, the lower levels then send forward bottom-up prediction errors to higher hierarchical levels (Friston and Kiebel, 2009).
In the setting of an oddball paradigm, every evoked response is considered to be a prediction error response. When the stimulus does not change, sensory learning enables better predictions to progressively attenuate the response to standards until an oddball or deviant stimulus is encountered — and a large prediction error is seen. The prediction error response to oddball stimuli may be further augmented by increasing the precision or gain of prediction errors: predictive coding models weigh sensory prediction errors by their precision, which is the inverse of the variance. Precision weighting of prediction errors, therefore, encodes the confidence in the accuracy of the information (i.e. prediction errors) in terms of their signal-to-noise ratio (Parr and Friston, 2017). In neurobiological terms, precision has been suggested to be mediated by synaptic gain modulation; namely, the sensitivity of postsynaptic responses to neural inputs (Adams et al., 2022), and is likely implemented by cholinergic neuromodulation (Moran et al., 2013; Yu and Dayan, 2005; Yu and Dayan, 2002). Increasing the gain of certain neuronal populations — e.g., superficial pyramidal cells reporting prediction errors – affords them privileged access to higher levels of processing and ensuing Bayesian belief updating. This can be likened to attentional gain in psychological terms – or increasing the Kalman gain in terms of Bayesian filtering.
According to the predictive coding model, there are, therefore, two mechanisms underlying the MMN/nMM (Auksztulewicz and Friston, 2016; Carbajal and Malmierca, 2018; Friston and Sporns, 2008; Harms et al., 2021; Parras et al., 2017). First, nMM could reflect the repetition suppression (RS) of the response to the predictable stimuli (standards) due to sensory learning – that resolves sensory prediction errors as predictions become more accurate. In addition, nMM could also reflect an increase in the precision or gain afforded prediction errors, when the (oddball) stimulus cannot be predicted; enabling the ascending prediction errors to revise higher-level representations about the auditory stream. In short, oddball stimuli elicit an enhanced neural response when an unexpected (deviant) stimulus is presented. Repetition suppression and prediction error modulation can be distinguished using control sequences (Ruhnau et al., 2012), and there is now evidence in both humans and rodents that MMN/nMM receives contributions from both prediction error (modulation) and repetition suppression (sensory learning) at various levels of the auditory system (Ishishita et al., 2019; Parras et al., 2017).
The main goal of the present study was to determine the role that ACh plays in the modulation of prediction errors and repetition suppression in the rat AC; utilizing control sequences to disambiguate these mechanisms (Nieto-Diego and Malmierca, 2016; Parras et al., 2017). We used microiontophoretic injections of ACh to determine how ACh neurotransmission affects the responses of AC neurons that exhibit nMM, and to establish the role of ACh on prediction error modulation and repetition suppression, in both primary and secondary AC areas.
Results
To explore the influence of cholinergic modulation on nMM, we recorded a total of 113 units in the AC before, during and after the microiontophoretic injection of acetylcholine in the rat. Recording depths ranged from 120–1191 μm including neurons from all layers.
To allocate each recorded neuron to a specific field in the AC, we recorded the frequency response area (FRA) and analyzed the topographical distribution of CF (characteristic frequency) for each unit. Each recording was assigned to a dorsoventral and rostrocaudal coordinate system relative to bregma as in previous studies (Nieto-Diego and Malmierca, 2016; Parras et al., 2017; Polley et al., 2007). This analysis allowed us to pool the data from all animals (Figure 1a) and construct a synthetic map of the CF across the entire rat auditory cortex (Nieto-Diego and Malmierca, 2016; Parras et al., 2017). Similar to these previous studies, we found a high-frequency reversal zone between ventral auditory field (VAF, caudally) and anterior auditory field (AAF, rostrally), a low-frequency reversal zone between A1 and posterior auditory field (PAF, dorso-caudally), and a high-frequency reversal between VAF and supra-rhinal auditory field (SRAF, ventrally). Thus, we could reliably define the lemniscal (A1, AAF, and VAF) and non-lemniscal (SRAF, PAF) auditory cortical fields as shown in Figure 1b.
ACh modulates the magnitude of neuronal responses and nMM
After the FRA was identified, we selected a pair of pure tones (10–30 dB above the minimum threshold) within the FRA at each recording site to test the adaptation of the neuronal response under the oddball paradigm. To do so we quantified nMM using the common SSA index (CSI). CSI values may range from –1 to +1, with negative values for units that respond more to standard stimuli and positive values for units that respond more to deviant stimuli (see Methods). In the following, we describe the effect of ACh on nMM (Figure 2; Figure 3; Figure 4; Figure 5; Figure 6), in terms of prediction error modulation and repetition suppression (Figure 7). Figure 2 illustrates representative responses of two units from the lemniscal auditory pathway — during and after the injection of ACh — that showed a change in the CSI during ACh release (Figure 2a and b, decreased and increased CSI, respectively). Similarly, the non-lemniscal, divisions (high-order regions) of the AC (Figure 3) also show decreased and increased CSI during ACh application (Figure 3a and b; respectively). In both Figures 2 and 3, peri-stimulus time histograms (PSTH) are also depicted for standard- and deviant stimuli (blue and red lines, respectively). On average, the firing rate decreased during the application of ACh (mean control: 5.8±7.0 spikes/s, mean ACh: 4.7±10.3 spikes/s; Wilcoxon signed-rank test, p<0.001). Most neurons recovered their basal firing rates after 60–90 min after the ACh injection (recovery mean: 6.5±8.1 spikes/s, Wilcoxon signed-rank test, p=0.4, control vs. recovery).
Given the changes in the FRAs observed in individual units during the application of ACh (Figures 2 and 3), we checked for bandwidth changes at the population level. For that, we measured the FRA bandwidth in octaves at 10 and 30 dB above the minimum threshold (BW10 and BW30, respectively). Numerous individual units showed either broadening or sharpening of their FRA bandwidth during ACh application, which were averaged out at the population level (Figure 4a, b). Indeed, no significant changes occur in neither BW10 (control: 0.76±0.57 oc taves; effect: 0.88±0.75 oc taves; Wilcoxon signed-rank test, p=0.248) nor BW30 (control: 1.33±0.73 oc taves; effect: 1.43±0.80 oc taves; Wilcoxon signed-rank test, p=0.664) (Figure 4c).
Next, we analyzed the effect of ACh on the CSI at the population level (n=85; comprising the neurons in our sample with significant responses to both f1 and f2 auditory stimuli). Figures 2 and 3 illustrate four examples of units from lemniscal and non-lemniscal subdivisions before, during and after the application of ACh. Modulation of CSI levels was observed in all individual examples. ACh produced a significant increase of the CSI value in 26% (Figure 5a, purple dots; 22/85, bootstrapping, 95% c.i.) and a significant decrease in 31% of the neurons recorded (Figure 5a, orange dots; 26/85, bootstrapping, 95% c.i.). Thirty-seven units (43%) were unaffected by ACh. The average CSI during ACh injection (CSI = 0.54±0.27) was similar to that of the control condition (CSI = 0.57±0.23; Wilcoxon signed-rank test, p=0.359; Figure 5b).
We analyzed the effect of ACh on the response to the deviant- or standard stimuli, separately. We observed that ACh produced a significant decrease in the mean firing rate in response to the deviant tones (Figure 5c; 11.75±17.46 Hz), as compared to the control condition (16.51±16.19 Hz; Wilcoxon signed-rank test, p<0.001), as well as in the mean firing rate in response to standard tones (Figure 5c; 4.58±6.38 vs. 3.97±9.90 Hz, control vs ACh; Wilcoxon signed-rank test, p=0.001).
In order to establish how ACh affects discriminability of deviant- and standard stimuli — i.e., how well the stimulus identity can be decoded from the neuronal firing rate – we calculated the area under the curve (AUC) of the receiver operating characteristic (ROC) for each tone (Figure 5d, left panel) and the Mutual Information to determine how much the response distributions for deviant- and standard stimuli overlap before and after ACh injections and how well the stimulus identity can be decoded from the neuronal firing rate. AUC indicates how different the neuronal responses to deviant- and standard stimuli are, and thus, how well a neuron is able to discriminate between both types of stimuli, considering that the neuron would respond differently to them (in terms of firing rate). Completely overlapping firing rate distributions (in response to deviants and standards) would yield an AUC = 0.5, while separate distributions would result in AUC = 1. On the other hand, Mutual Information is defined as the reduction in the uncertainty of one variable (in this case the neuronal response) due to the knowledge of another (in this case, the type of stimulus; Timme and Lapish, 2018). In other words, this measure reflects how much information a neural response (firing rate) contains about the type of stimulus that is played. In practical terms, high Mutual Information implies that a neuron is able to discriminate between stimuli, and thus respond differently to them. This analysis shows that the mean AUC decreased significantly during the application of ACh (0.74±0.13 vs. 0.67±0.15, control vs. effect; Wilcoxon signed-rank test, p<0.001). The same happened for Mutual Information (Figure 5d, right panel; 0.09±0.07 vs. 0.07±0.08, control vs. effect; Wilcoxon signed-rank test, p=0.003). These AUC and Mutual Information analyses indicate that under the influence of ACh, neurons are less able to discriminate deviant stimuli from standard stimuli.
Different mechanisms for SI increase or decrease during ACh application
As mentioned above, the CSI is an indicator of the average level of nMM for a single unit, since it takes into account the responses of two frequencies (f1 and f2) that alternatively play the role of deviant and standard. However, in order to better understand the effect of ACh on the response to each stimulus class — depending on whether it was presented as a deviant or as a standard — the frequency-specific SSA index (SI, see Methods) is a better indicator. We calculated the SI for 205 individual frequencies, from 113 units that elicited a significant response to that particular frequency from our sample. This analysis shows that ACh produced a significant increase of the SI in 27% of the frequencies recorded (Figure 6a, purple dots; 55/205, bootstrapping, 95% c.i.), 31% (64/205; Figure 6a, orange dots) showed a significant decrease of their SI and 42% frequencies (86/205) were unaffected by ACh. The mean SI (SI = 0.55±0.31, control condition) was reduced during ACh injection (SI = 0.40±0.50) although this decrease was not significant (Wilcoxon signed-rank test, p=0.086; Figure 6b).
Since we found that the application of ACh had divergent effects on the SI (increase or decrease), we further analyzed separately those two groups to identify different mechanisms. We examined the firing rates in response to deviant- and standard stimuli, separately for those frequencies showing a SI increase or decrease during the application of ACh. In both cases, responses to the deviant- and standard stimuli were affected differently, when measuring the firing rate change ratio (amount of change caused by ACh, relative to the baseline firing rate). In frequencies showing SI decrease, the effect of ACh was significantly different for the deviant- and standard stimuli (Figure 6c, top panel; p<0.001, Wilcoxon signed-rank test). Furthermore, this ratio was clearly negative for deviant sounds (p<0.001, Wilcoxon signed-rank test) and not significant for standards (p=0.256, Wilcoxon signed-rank test), indicating that ACh caused a reduction of the firing rate differentially during deviant stimuli. On the other hand, in frequencies showing SI increase the effect of ACh also had a distinct effect on the deviant- and standard stimuli (Figure 6d, top panel; p<0.001, Wilcoxon signed-rank test), but in this case, the effects were reverse: the change ratio was clearly negative for standard sounds (p=0.001, Wilcoxon signed-rank test) and not significant for deviant ones (p=0.379, Wilcoxon signed-rank test), hence ACh caused a reduction of the firing rate primarily during standard stimuli. In other words, while the effect of ACh was a reduction of the firing rate, in some frequencies it only affected responses to the deviant stimuli, thus decreasing SI (Figure 6c, top panel), while in others it only affected responses to the standard stimuli, resulting in an increment of the SI (Figure 6d, top panel).
Interestingly, stimuli discriminability was affected differently by ACh in both groups. In frequencies showing SI decrease (Figure 6c, bottom panel) both ROC AUC and Mutual Information decreased significantly during ACh application (p<0.001 in both cases, Wilcoxon signed-rank test), in line with the results obtained for the whole sample (compare to Figure 5d). However, in frequencies showing SI increase (Figure 6d, bottom panel) ROC AUC did not change significantly (p=0.339), while Mutual Information increased slightly (p=0.008). These results suggest that frequencies showing SI decrease are the ones mediating the main effects of ACh in the population.
ACh affects neuronal mismatch, modulating prediction error responses but not repetition suppression
According to the predictive coding framework, nMM can be split into two functionally distinct components: prediction error modulation (i.e. an increased response to an oddball caused by the violation of a regularity) and repetition suppression (i.e. reduction in the response to a standard caused by sensory learning of a repeated stimulus). These differential effects can be quantified by comparing the responses to the oddball paradigm with a suitable control sequence (Figure 7—figure supplement 1). Therefore, to test if ACh modulates prediction error signaling in AC, we also recorded the cascade control sequence (for the corresponding oddball sequences) in 66 neurons of the sample, using an approach previously published (Carbajal and Malmierca, 2018; Parras et al., 2017; Ruhnau et al., 2012; Valdés-Baizabal et al., 2019). First, we calculated the normalized spike counts in response to the cascade control, deviant, and standard stimuli (CAS, DEV, and STD, respectively; see Methods). ACh led to a significant decrease in the DEV responses in those units where SI decreased (Figure 7a; Wilcoxon Signed Rank test, p<0.001) and to a significant increase in those units where SI increased (Figure 7b; p=0.011). The opposite effect occurs on STD stimuli, increasing in units where SI decreased (Figure 7a; p<0.001) and decreasing in units where SI increased (Figure 7b; p<0.001). While this trend is very clear, there were no significant changes in CAS (p=0.091 and p=0.072, for SI decrease and increase, respectively).
Using the normalized CAS, DEV, and STD responses, we computed the indexes iMM, iPE, and iRS before and after the ACh injection in SI decrease units (Figure 7c, iMM: 0.63±0.19 vs. 0.13±42; iPE: 0.29±0.27 to -0.05±0.38; iRS: 0.34±0.20 vs. 0.17±0.49; control vs. ACh, respectively) and SI increase units (Figure 7d, iMM: 0.47±0.31 vs. 0.67±0.22; iPE: 0.15±0.36 vs. 0.33±43; iRS: 0.32±0.30 vs. 0.34±0.31). iMM is equivalent to the CSI, but using normalized response values, and since it includes information about the cascade control, it can be divided into prediction error and repetition suppression components (iPE and iRS, respectively; see Methods and Figure 7—figure supplement 1c). Significant effects of ACh administration were found for iMM and iPE, both in SI decrease units (iPE: p<0.001; iMM: p<0.001; Wilcoxon signed-rank test; Figure 7c) and SI increase units (iPE: p=0.023; iMM: p<0.001; Wilcoxon signed-rank test; Figure 7d), while iRS remained unchanged (SI decrease units iRS: p=0.121; SI increase units iRS: p=0.675; Wilcoxon signed-rank test).
Under the predictive coding framework, the precision of prediction errors encodes the confidence or reliability afforded such errors (Friston and Sporns, 2008), enabling precise prediction errors to access higher levels of the auditory hierarchy. To test whether neuronal responses show differential effects of precision weighting — depending on their level of deviance sensitivity – we split our sample into two similar-sized groups (of frequencies), using the median SI as cutoff (0.6143: SI before the application of ACh): low SI (n=65 frequencies with cascade control) and high SI (n=66 frequencies with cascade control) (Figure 8a–c). Given this grouping, we found that ACh decreased significantly iMM (0.72±0.13 vs. 0.52±0.36, control vs. effect; p<0.001) and iPE (0.33±0.33 vs. 0.145±0.48, control vs. effect; p=0.004) in high SI frequencies, while only decreasing iRS significantly in low SI frequencies (0.28±0.25 vs. 0.18±0.36, control vs. effect; p<0.037) (Figure 8a–c, top row).
Given the divergent defects of ACh on putative prediction error modulation, we hypothesized that acetylcholine would increase the variance of deviance-sensitive responses over groups of cells. We calculated the variance of the above distributions of predictive coding indexes (iMM, iPE, and iRS) as an estimate of the (inverse of) their precision, as well as its 95% confidence interval using a 1000-sample bootstrap strategy, and then used an F-test to check for differences in the variances before and during ACh application (Figure 8a–c, bottom row). These results demonstrate that ACh application significantly increases the variance of iMM and iRS for both low SI and high SI groups, but interestingly ACh also significantly increases the variance of iPE only for the high SI group.
Then, we repeated the same variance analysis grouping the frequencies based on the anatomical layer location (depth) of the recording sites (Figure 8d–f). In this case, we split the frequencies into deep and superficial, with a depth cutoff of 650 µm, which is near the boundary between layers IV and V in the rat AC (Games and Winer, 1988; Malmierca, 2015). In this case, we found no significant effect of ACh on iMM, iPE, nor iRS for none of the depth groups (Figure 8d–f, top row). However, the general trend was an increase in the variance for all indexes and depth groups during the application of ACh (except for iPE in superficial units), which was significant for all indexes in the case of deep units only (Figure 8d–f, top row).
Thus, our results indicate that ACh produces a general increase in the variance for iMM, iPE, and iRS. This general increase occurs regardless of the baseline SI value, but is more noticeable for neurons located in the infragranular layers.
ACh effect on nMM as a function of topographic distribution
A remaining question is whether the effect of ACh is uniformly distributed across AC neurons in different fields and layers. Figure 9 shows the CSI (first row) and spike count levels for the standard and deviant tones (middle and bottom rows, respectively) before and after ACh application, as well as the difference between the ACh and control conditions. The highest levels of CSI in the control condition (left column) were found in SRAF. ACh application (middle column) produced an overall increment of CSI in AAF and a decrement in the rest of the fields, however, such changes did not reach statistical significance (p>0.79 in all fields, Wilcoxon signed-rank test with Holm-Bonferroni correction). Regarding the responses to the standard stimuli (middle row), ACh application caused an increase of the firing rate in A1 and PAF, and a decrease in AAF, VAF, and SRAF, but such changes were only significant in AAF (p=0.001, Wilcoxon signed-rank test with Holm-Bonferroni correction). By contrast, the firing rates in response to the deviant sounds (bottom row) decreased in all fields during ACh application, with significant changes in AAF, VAF, and SRAF (p=0.011, p=0.026, and p=0.009, respectively; Wilcoxon signed-rank test with Holm-Bonferroni correction).
In order to better understand the effect of ACh on the responses to the individual frequencies of the oddball paradigm (f1 and f2), we calculated the SSA index (SI) independently for each of the paradigm frequencies. We found that ACh was more likely to cause an SI decrease in all AC fields except for AAF, where ACh was more likely to cause an SI increase (Figure 10a; Figure 10b) shows the distribution of recording depths within the AC, for those units that showed a significant SI decrease (orange) or increase (purple) during the application of ACh. The shaded background shows the probability density function for those distributions. We found that frequencies that experienced an SI decrease during ACh application were most likely located in the granular layer (layer IV) or the upper infragranular regions (dorsal parts of layer V). In contrast, frequencies that experienced an SI decrease during ACh application were more likely located in deeper infragranular regions (layers V and VI).
Discussion
We recorded neuronal responses from the AC under an auditory oddball paradigm while performing microiontophoretic applications of ACh to investigate the role of the cholinergic system in the processing of unexpected or surprising events. Furthermore, we employed cascade sequences to control for repetition suppression and assess cholinergic effects on prediction error relative to repetition suppression. Our results show that (1) ACh modulates nMM, evidenced by decrements or increments of the neuronal mismatch indexes (iMM, iPE, and iRS). The decreased SI is due to a diminished response to the deviant tone while increased SI is the consequence of a diminished response to the standard sounds (Figure 6, and examples in Figures 2–3); (2) ACh exerts a global and distributed effect on AC neurons that is layer-specific but field unspecific (Figures 9 and 10); (3) The decreased- and increased SI are driven by the effects of ACh on iPE but not iRS (Figure 7), and (4) ACh increases the variance of iMM, iPE, and iRS. Together, these findings reveal how ACh modulates prediction error signals and repetition suppression differentially in the AC, likely gating these signals, as they ascend beyond the auditory cortex to frontal, higher cognitive regions.
ACh increases the amplitude of AC neuronal responses to sound stimulation, decreases the auditory threshold, and sharpens the receptive field in the AC. (Edeline, 2012; Irvine, 2018a; Irvine, 2018b; Kilgard and Merzenich, 1998; Ma and Suga, 2005; Metherate, 2011; Metherate et al., 1992; Puckett et al., 2007). These ACh-mediated effects are produced by a rapid disinhibition of neuronal responses, modifying synaptic strength, enhancing excitatory-inhibitory balance, and reorganizing cortical circuits promoting cortical plasticity (Froemke et al., 2007; Irvine, 2018a; Irvine, 2018b; Metherate, 2011; Picciotto et al., 2012). But ACh also produces direct inhibitory effects leading to a decrease of cortical responses (Colangelo et al., 2019; Dasgupta et al., 2018; Eggermann and Feldmeyer, 2009; Qi and Feldmeyer, 2022). Our results confirm these mixed effects (excitatory and inhibitory, during the application of ACh) at the neuronal level. However, when looking at the population level, the overall effect in our sample was a differential reduction of the responses to either deviant or standard stimuli. Apart from the direct hyperpolarizing effect of ACh on some cell types, it is also possible that some of the recorded neurons received direct inhibitory inputs, and those inhibitory neurons were the ones affected by ACh. By increasing the firing rates of such inhibitory interneurons, the observed overall effect of ACh could lead to a reduction of the activity of the recorded neurons. An important functional motif — for state modulation in the cortex — is through disinhibition consisting of parvalbumin-, somatostatin-, and vasoactive intestinal peptide- expressing interneurons. When the latter receives cholinergic projections from the basal forebrain, they remove the inhibition on layer 2/3 pyramidal neurons exerted by somatostatin interneurons (Fu et al., 2014). The connectivity relationships among inhibitory neurons seem to be critical. Moreover, it has been recently described how each of those cell types responds differently to deviants and standards in the auditory cortex (Yarden et al., 2022). ACh plays important roles in arousal, attention, and sensory learning (Hasselmo, 1999; Hasselmo and Sarter, 2011; Metherate, 2011; Picciotto et al., 2012; Sarter et al., 2001; Weinberger, 2004). Thus, it is not surprising that ACh may have a critical role in shaping nMM, which has many properties in common with behavioral habituation to a repeated stimulus and can be considered a simple form of learning (Irvine, 2018a; Nelken, 2014; Netser et al., 2011). Furthermore, a recent study in rats has shown that the administration of muscarinic ACh receptor agonists and antagonists affected MMN amplitude in A1 and PAF (Schöbi et al., 2021), suggesting that the cholinergic system modulates the connectivity between those areas during the occurrence of deviant sounds. But to the best of our knowledge, the present study demonstrates a specific effect of ACh on prediction error responses.
In the IC, ACh selectively increased the response to the repetitious standard stimulus (Ayala and Malmierca, 2015), while we show here that ACh in the AC not only increases the responses to the unexpected deviant sounds leading to an increase of SI (Figure 6d) but also decreases responses to the repeated standard that result in a diminished CSI index (Figure 6c). This divergent effect of ACh on subcortical and cortical nMM may be related to the different origin of ACh and/or the unique organization of neuronal circuitries in IC and AC (Figure 11). While the sources of ACh to IC emerge from the pedunculopontine and laterodorsal tegmental regions in the brainstem (Motts and Schofield, 2009), those to AC originate in the basal forebrain (Chavez and Zaborszky, 2017).
According to the predictive coding formulations, cortical areas send predictions to lower hierarchical areas to inhibit neuronal populations encoding prediction errors. When top-down predictions match bottom-up sensory input, the prediction error response is attenuated. Conversely, those lower areas send error information to higher hierarchical centers when event predictions fail to drive higher-level representations and improve the top-down predictions. If a stimulus is not predicted, the difference between predicted and received inputs yields large prediction errors, while if a stimulus is constantly repeated, sensory learning enables more accurate predictions and a decrease in prediction error responses – a phenomenon referred to as repetition suppression (Auksztulewicz and Friston, 2016). Thus, nMM can be explained under the predictive coding framework (Carbajal and Malmierca, 2018; Carbajal and Malmierca, 2020; Parras et al., 2017; Malmierca et al., 2019), and recent studies that found LFP and spiking responses to omitted tones in the rat auditory cortex (Auksztulewicz et al., 2023; Lao-Rodríguez et al., 2023), among others, strongly support this interpretation. However, predictions and prediction errors may also be modulated in a context-dependent manner, meaning that responses can be modulated and given precedence depending on the context in which the stimulus is perceived (Keller and Mrsic-Flogel, 2018). The source of such a modulating or gating signal may arise from neuromodulatory inputs (Fu et al., 2014; Moran et al., 2013; Pinto et al., 2013; Polack et al., 2013; Thiele and Bellgrove, 2018) that can not only gate plasticity (Kilgard and Merzenich, 1998; Martins and Froemke, 2015; Weinberger, 2004), but also change the balance of top-down versus bottom-up influence (Yu and Dayan, 2002).
Within the predictive coding framework, classical neuromodulators such as ACh are often thought to increase the precision of prediction error signaling (Auksztulewicz et al., 2018; Auksztulewicz and Friston, 2016; Feldman and Friston, 2010; Moran et al., 2013). Indeed, Yu and Dayan, 2005; Yu and Dayan, 2002 proposed that ACh levels reflect the ‘expected uncertainty’ (when there is less confidence in the prediction) associated with top-down information and modulate the interaction between top-down and bottom-up processing in determining the appropriate neural representations for inputs (Auksztulewicz et al., 2017; Yu and Dayan, 2005; Yu and Dayan, 2002). Accordingly, the precision of prediction errors would encode the confidence or reliability afforded by such errors (Friston, 2005), increasing the influence of precise prediction errors on higher-level processing. Thus, the canonical view – within the predictive coding framework — is that precision would affect postsynaptic gain through the action of neuromodulators such as ACh, enhancing the responses of units broadcasting prediction errors. Using computational neuronal modeling and EEG recordings under oddball paradigms in subjects manipulated pharmacologically with galantamine a competitive inhibitor of acetylcholinesterase, Moran et al., 2013 have shown that ACh enhances the precision of bottom-up synaptic transmission in cortical hierarchies by optimizing the gain of supragranular pyramidal neurons, amounting to an increase in the precision of prediction error signaling. Thus, our results and those from Moran and colleagues (2013) suggest that ACh mediates the representation of precision.
Our results are consistent with this predictive coding picture; particularly when noting that an increase in precision leads to a representational sharpening (Kok et al., 2012a), in which the majority of (prediction error) neurons are suppressed and a few neurons encoding precise prediction errors are enhanced (Friston, 2018; Kok et al., 2012a; Shipp, 2016). This is consistent with the effect of acetylcholine on average over neurons studied and the differential effect on standard and oddball responses; dependent upon their SI.
ACh increases the variance of the prediction error component of neuronal responses in those units with higher SI (Figure 8b, top), especially in units with strong deviance detection properties (Figure 8b, bottom). This finding suggests that ACh plays a key role in encoding precision or uncertainty. Interestingly, ACh not only increases the variance of prediction error responses, but also the variance of repetition suppression (Figure 8c).
A change in the precision of prediction errors may denote — or alert to — a variable environment, where it is difficult to extract rules and only generate weak (and less useful) predictions. Under this scenario, some fluctuations in the sensory input can be expected, assigning less precision to prediction errors (Yon and Frith, 2021). If so, ACh would adjust the predictions’ confidence, so we would not be surprised by small changes in the environment, since they are likely, expected noisy, fluctuations. Otherwise, an excessive confidence on predictions in a volatile environment may override the incoming evidence, leading to the false inference of expected by absent events (Yon and Frith, 2021). Attention is thought to boost the precision of prediction errors (Kok et al., 2012b), and our experiments carried out under anesthesia clearly lack attention, where baseline cholinergic activity is low, making it difficult to interpret how the results generalize to more ecologically valid situations (Minces et al., 2017). Interestingly, there is a wealth of data suggesting that minimal impairments to the encoding or defective signaling of precision or uncertainty are sufficient to explain psychopathic conditions such as schizophrenia or autism (Friston, 2023). Furthermore, while we are not aware of any known direct connection between auditory processing disorders and ACh, individuals with auditory processing disorders do have difficulties with auditory selective attention, so perhaps one could speculate that ACh, by modulating SSA/prediction error, could impact encoding of salient events, and if disrupted could lead to problems with selective attention. Moore, 2006; Moore, 2012 speculated that auditory processing disorders may arise from unbalanced processing in bottom-up and top-down contributions, i.e., there might be some aspects of predictive processing altered. A recent study (Felix et al., 2019) analyzed altered temporal processing at the level of the brainstem in α7-subunit of the nicotinic acetylcholine receptor (α7-nAChR)-deficient mice. After studying α7-nAChR knockout mice of both sexes and wild-type colony controls, they concluded that the malfunction of the CHRNA7 gene that encodes the α7-nAChR may contribute to degraded spike timing in the midbrain, which may underlie the observed timing delay in the ABR signals. These authors (Felix et al., 2019) proposed that their findings are consistent with a role for the α7-nAChR in types of neurodevelopmental and auditory processing disorders. There is also evidence of cholinergic system disfunction being related to the pathophysiology of Alzheimer’s disease (Pérez-González et al., 2022). For instance, disfunction of the synapses of cholinergic neurons in the hippocampus and nucleus basalis of Meynert, as well as decreased choline acetyltransferase activity, is associated with memory disorders in Alzheimer’s disease (Hampel et al., 2018). Also, Alzheimer’s disease patients show reduced amounts of the vesicular ACh transporter in some brain areas (Aghourian et al., 2017). Finally, cholinesterase inhibitors seem to have some favorable effects in the treatment of Alzheimer’s disease patients (Sharma, 2019).
It is unclear to what degree the amount of agonist liberated used iontophoretic application of cholinergic agonists emulates a relevant physiological condition (Minces et al., 2017). The application of exogenous ACh, as we used here, could potentially surpass in physiological concentrations, making the results hard to extrapolate to a behaving animal. Furthermore, future studies should investigate the effects of different doses of ACh on auditory processing to determine if there is a dose-response relationship as well as the effect of muscarinic- and nicotinic receptors blockade to theoretically observe more akin, analogous physiological effects. Similarly, forthcoming studies using optogenetic techniques to activate only basal forebrain cholinergic neurons while also recording from auditory cortex neurons will allow a fine-grained investigation of the role of cortical acetylcholine on sensory representation. Moreover, because our study was carried out in anesthetized animals, it precludes generalization of the effects of ACh on behaving states (Quintela-Vega et al., 2023). Future studies using an awake preparation will confirm a definitive corroboration of the effect of ACh on precision and determine the behavioral relevance of these findings. The available neurobiological evidence suggests that precision is encoded by the synaptic gain of superficial pyramidal cells encoding prediction errors, which are controlled by neuromodulatory systems (e.g. dopaminergic, cholinergic) and/or synchronized neural activity (Feldman and Friston, 2010; Kok et al., 2012b). However, we do not find changes in prediction error precision in the superficial layers; our analysis suggests instead that the effects of ACh were more limited to the deep AC layers (Figure 8e). We observed a distinct and differential distribution of units throughout the cortical layers, contingent upon the influence of ACh on their selectivity index (SI). Units exhibiting an SI decrease during ACh application were predominantly located and biased in layers III-IV and its adjacent regions, while those demonstrating an SI increment during ACh application spanned from the deeper layer IV to layer VI (Figure 10b). In the initial conceptualizations of predictive coding, it was posited that prediction units reside in the infragranular layers (layers V and VI), with error units situated in supragranular layers (layers I – III) (Bastos et al., 2012; Friston, 2005). Alternatively, some models suggest that error units would be localized in the granular layer (layer IV), while prediction units would be detected in both supra- and infragranular layers biased competition model of predictive coding (Spratling, 2008). In our recordings, units displaying an SI reduction during ACh application were also those with a higher SI under baseline conditions (i.e. larger response to deviant than to standard stimuli), thereby potentially qualifying as error units. These units were identified in layer IV and neighboring areas, which appears to support the model proposed by Spratling, 2008. Furthermore, according to Eggermann and Feldmeyer, 2009, the application of ACh would filter weak sensory inputs to layer IV and amplify signals at superficial and deep layers, effectively increasing the signal-to-noise ratio. In contrast, units demonstrating an SI increase during ACh application could be considered potential prediction units; such units were detected in granular and deep layers but were absent in superficial layers. This observation aligns with both predictive coding models, as each postulates that deep layers contain prediction units. A very interesting study would have been to estimate whether the neurons we recorded were excitatory or inhibitory based on the spike shape. Indeed, we have previously shown that both types of responses show prediction error (Pérez-González et al., 2021) but ACh may affect them differentially and this potential divergence may explain some of the variance of the data. This would involve implementing a similar analysis as done in our previous study (Pérez-González et al., 2021), however, that study used single-unit recordings and here we used multi-unit recordings. In the present data, the recordings captured occasional spikes from other neurons which introduced alterations in the average spike shape, and thus precluded an accurate categorization. Therefore, we decided to discard this analysis for the sake of accuracy and appropriate scrutiny awaits future studies.
In conclusion, we have shown that ACh plays a multifold role in modulating nMM in AC, affecting the precision of prediction error signaling and gating prediction errors to hierarchically higher processing levels.
Materials and methods
Subjects and surgical procedures
Request a detailed protocolThe experimental protocols were approved conforming to the University of Salamanca Animal Care Committee standards and the European Union (Directive 2010/63/EU) for the use of animals in neuroscience research. Experiments were performed on 37 adult female Long-Evans rats with body weights within 180–250 g. We used female rats as our previous work is based on this sex but no attempt to study gender dimension was aimed, which would have required a larger number of animals. Surgical anesthesia was induced and maintained with urethane (1.5 g/kg, intraperitoneal), with supplementary doses (0.5 g/kg, intraperitoneal) given as needed. Dexamethasone (0.25 mg/kg) and atropine sulfate (0.1 mg/kg) were administered at the beginning of the surgery to reduce brain edema and the viscosity of bronchial secretions, respectively. Normal hearing was verified with auditory brainstem responses (ABR) recorded with subcutaneous needle electrodes, using a RZ6 Multi I/O Processor (Tucker-Davis Technologies, TDT) and processed with BioSig software (TDT). ABR stimuli consisted of 100 µs clicks at a 21 /s rate, delivered monaurally to the right ear in 10 dB steps, from 10 to 90 decibels of sound pressure level (dB SPL), using a closed-field speaker. After the animal reached a surgical plane of anesthesia, the trachea was cannulated for artificial ventilation and a cisternal drain was introduced to prevent brain hernia and brain edema. Isotonic glucosaline solution was administered periodically (5–10 ml every 7 hr, subcutaneous) throughout the experiment to prevent dehydration. Body temperature was monitored with a rectal probe thermometer and maintained between 37 and 38°C with a homeothermic blanket system.
The auditory cortex surgery was described in previous works (Nieto-Diego and Malmierca, 2016; Parras et al., 2017). The temporal bone was exposed and the auditory cortex was located using stereotactic coordinates (Murone et al., 1997). A craniotomy was performed over the auditory cortex, the dura was removed carefully, and the exposed area was filled with a layer of agar to prevent desiccation and to stabilize the recordings. Before applying the agar, a magnified picture (25x) of the exposed cortex was taken with a digital camera coupled to the surgical microscope (Zeiss) through a lens adapter (TTI Medical). The picture included a pair of reference points previously marked on the dorsal ridge of the temporal bone, indicating the absolute scale and position of the image with respect to bregma (the reference point). This picture was displayed on a computer screen and was overlapped with a micrometric grid, over which the placement of the multibarrel for every recording was marked. The micrometric grid allowed the generation of coordinates in a two-dimensional axis, from the superimposed image of the auditory cortex of each animal. Once the coordinates of each of the recorded units of all the animals were obtained, we used the functions ‘griddata’ and ‘contourf’ of MATLAB to generate topographic maps (through a graduation of colors) of the characteristic frequency (CF), CSI levels, and deviant and standard tone responses in all auditory cortical fields.
Electrophysiological recording and microiontophoresis
Request a detailed protocolA tungsten electrode (1–3 MΩ) was used to record multiunit neuronal activity. It was attached to a five barrel multibarrel borosilicate glass pipette that carried the drug solution to be delivered in the vicinity of the recorded neuron. The multibarrel’s tip was cut to a diameter of 20–30 µm. The center barrel was filled with saline solution for current compensation (165 mM NaCl) whereas the other barrels were filled with 1 M acetylcholine chloride (Sigma, catalog no. A6625). The pH was adjusted between 4.0–4.2 (Ayala and Malmierca, 2015). Drugs were retained by applying a –15 nA current, and were ejected when required, typically using 30–40 nA currents (Neurophore BH-2 System, Harvard Apparatus). The duration of the drug ejection usually lasted 8–10 min and the recording protocols were extended until the effect of the drug had disappeared (60–90 min approx.). The multibarrel assembly was positioned over the pial surface of the auditory cortex, forming a 30° angle with the horizontal plane towards the rostral direction, and advanced using a piezoelectric micromanipulator (Sensapex) until we observed a strong spiking activity synchronized with the train of search stimuli made of white noise. Analog signals were digitalized with a RZ6 Multi I/O Processor, a RA16PA Medusa Preamplifier, and a ZC16 headstage (TDT) at a 12 kHz sampling rate and amplified 251 x. Neurophysiological signals for multiunit activity were band-pass filtered between 0.5 and 4.5 kHz.
Experimental design and stimulation paradigms
Request a detailed protocolSound stimuli were generated using the RZ6 Multi I/O Processor (TDT) and custom software programmed with the OpenEx Suite (TDT) and MATLAB. Sounds were presented monaurally through a speaker, in a closed-field condition to the ear contralateral to the left auditory cortex. We calibrated the speaker to ensure a flat spectrum up to ~75 dB SPL between 0.5 and 44 kHz; the second and third harmonics were at least 40 dB lower than the fundamental at the loudest output level for all the frequencies. The experimental stimuli were pure tones in the range 0.5–44 kHz, with a duration of 75 ms, including 5 ms rise/fall ramps presented at a rate of 4 stimuli/s. Once a suitable neuron was found, the frequency response area (FRA) of the cell, i.e., the combination of frequencies and intensities capable of evoking a suprathreshold response, was obtained automatically using a randomized paradigm that presented tones between 0.5–44 kHz in 25 logarithmic steps, with intensities spaced by 10 dB steps, from 0 to 70 dB SPL, at a repetition rate of 4 /s. Based on this information, we selected a pair of frequencies evoking similar responses at 10–30 dB above the threshold. We used pure tones at these frequencies as the stimuli in the oddball paradigm. Oddball sequences consisted of frequently repeating stimuli (standard tones) which were pseudo-randomly interleaved with rare events (deviant tones). Two oddball sequences with fixed parameters (400 trials each, 75 ms stimulus duration, 0.5 octaves frequency separation, 10% deviant probability, 250 ms onset to onset, and a minimum of three standard tones before a deviant) were presented for every pair of stimuli thus selected. In one of the sequences, the low frequency (f1) was the ‘standard’ and the high frequency (f2) was the ‘deviant,’ and in the other sequence, their roles were inverted. The order of presentation of these two sequences was randomized across sites. In some cases, one or more extra pairs of stimuli were selected, and the oddball sequences were repeated with the new stimuli.
According to the predictive coding framework, mismatch responses like those obtained during an oddball paradigm can be divided into two components: repetition suppression, a reduction in the response caused by a repeated stimulus, and prediction error, an increased response caused by the violation of a regularity (Carbajal and Malmierca, 2018). In a subset of the experiments, we used control sequences to evaluate the separate contribution of repetition suppression and prediction error. Control sequences consisted of 10 tones evenly spaced by 0.5 octaves (same as in the oddball sequences), including the tones used in the oddball paradigm, and all stimuli at the same previously chosen sound level. Each control sequence lasted 400 trials, the duration of all stimuli was 75 ms, and the presentation rate 4 /s. We used two different control sequences, namely the many-standards and cascade sequences (Figure 7—figure supplement 1). The many-standards control is the consecutive presentation of blocks of 10 tones randomly ordered within the block, each tone with a 10% probability of occurrence (Schröger and Wolff, 1996). In this sequence the tones are unpredictable, and it is not possible to establish a rule which could be used to predict the following tones. On the other hand, the cascade control consists of the regular presentation of the same 10 tones in ascending or descending frequency succession (Ruhnau et al., 2012). This sequence also avoids the effects of the repetition of a single stimulus, but in this case, it maintains a predictable context. Cascade sequences are considered a more rigorous control than many-standard sequences because, like the oddball sequence, a regularity is established (Carbajal and Malmierca, 2018).
Data analysis
Request a detailed protocolWe calculated the bandwidth of the FRAs (Figure 4) for each recorded unit by measuring the width of the response area in octaves (relative to the low-frequency border), at 10 and 30 dB above the minimum threshold (BW10 and BW30, respectively).
The degree of nMM was quantified by the CSI, reported previously (Ayala et al., 2015; Duque et al., 2016; Malmierca et al., 2009; Ulanovsky et al., 2004; Ulanovsky et al., 2003). The CSI reflects the difference between the neural responses to the deviant and standard stimuli, normalized to the total of responses to both stimuli, and is defined as:
where and are responses to each frequency or when they were the deviant (d) or standard (s) stimulus in the oddball paradigm, respectively.
Similarly, a frequency-specific SSA index (SI) was calculated for each of the frequencies presented in an oddball paradigm, as:
Where d(f) and s(f) are responses to frequency f when it was the deviant (d) or standard (s) stimulus in the oddball paradigm, respectively.
For the subset of experiments where we recorded the many-standards and cascade controls, we compared the responses to the same physical stimulus when it took the role of a standard, a deviant (ascending or descending, depending on whether the preceding standard was of lower or higher frequency, respectively), or as part of a cascade sequence (matching the ascending or descending condition of the corresponding deviant). Alongside the oddball paradigm, we recorded responses of neurons to two cascade sequences (ascending and descending), which consisted of 10 tones selected within the FRA presented in a predictable succession of increasing or decreasing frequencies. We did not include the many-standards control in these analyses because these experiments are time-consuming and the need of holding the recording neuron for long enough before, during and after the drug injection. However, we have previously demonstrated that the many-standards and cascade control responses were comparable, and the latter is considered to be more rigorous. (Casado-Román et al., 2019; Parras et al., 2017; Ruhnau et al., 2012; Valdés-Baizabal et al., 2019).
The responses of each neuron were normalized for each tested frequency as
where,
is the Euclidean norm of the vector defined by the deviant, standard, and cascade responses, so that the normalized responses take values in the range of 0–1.
With these normalized baseline-corrected spike counts, we next computed the indices of neuronal mismatch (iMM), repetition suppression (iRS), and prediction error (iPE) as:
Index values ranged between –1 and 1 and facilitated the quantitative decomposition of neuronal mismatch into repetition suppression and prediction error since
Which is largely comparable to the CSI calculated for the oddball paradigm (Parras et al., 2017). In order to determine a significant effect of the drugs on the CSI (Figure 5a), we calculated an empirical distribution of CSI values by performing 2000 bootstraps of the responses to each standard or deviant stimulus for each neuron in the control condition, from which we determined a 95% confidence interval (vertical lines in Figure 5a). An effect was considered significant at α=0.05 if the CSI during the drug condition fell beyond the limits of the control CSI confidence interval (red dots in Figure 5a). An equivalent approach was followed for the SI values in Figure 6a.
The area under the ROC curve (AUC) was calculated using the firing rates in response to deviant or standard stimuli as predictors (Figures 5d, 6c and d). Mutual Information on the firing rates in response to deviant or standard stimuli was calculated using the Information Breakdown Toolbox (Magri et al., 2009).
The confidence intervals for the variance of the different index distributions (Figure 8), were obtained by performing 1000 bootstraps.
Since parameters did not follow a normal distribution, non-parametric statistics were used: Mann-Whitney (for independent data) or Wilcoxon signed-rank test (for paired data). Multiple comparisons were corrected using the Holm-Bonferroni approach. Unless stated otherwise, all data are reported as mean ± SD. All of the data analyses were performed with Sigma Plot v12.5 and MATLAB software, using the built-in functions, the Statistics and Machine Learning toolbox, or custom scripts and functions developed in our laboratory. Custom software used in this paper is available upon reasonable request to the authors.
Data availability
Source data in the manuscript contains the numerical data used to generate the figures.
References
-
Repetition suppression and its contextual determinants in predictive codingCortex; a Journal Devoted to the Study of the Nervous System and Behavior 80:125–140.https://doi.org/10.1016/j.cortex.2015.11.024
-
Cholinergic modulation of stimulus-specific adaptation in the inferior colliculusThe Journal of Neuroscience 35:12261–12272.https://doi.org/10.1523/JNEUROSCI.0909-15.2015
-
Modulation of cortical circuits by top-down processing and arousal state in health and diseaseCurrent Opinion in Neurobiology 52:172–181.https://doi.org/10.1016/j.conb.2018.06.008
-
The Senses: A Comprehensive ReferenceNovelty processing in the auditory system: detection, adaptation or expectation, The Senses: A Comprehensive Reference, Elsevier, 10.1016/B978-0-12-809324-5.24154-0.
-
Cellular, synaptic and network effects of acetylcholine in the neocortexFrontiers in Neural Circuits 13:24.https://doi.org/10.3389/fncir.2019.00024
-
Prefrontal executive and cognitive functions in rodents: neural and neurochemical substratesNeuroscience and Biobehavioral Reviews 28:771–784.https://doi.org/10.1016/j.neubiorev.2004.09.006
-
Beyond traditional approaches to understanding the functional role of neuromodulators in sensory corticesFrontiers in Behavioral Neuroscience 6:45.https://doi.org/10.3389/fnbeh.2012.00045
-
Attention, uncertainty, and free-energyFrontiers in Human Neuroscience 4:215.https://doi.org/10.3389/fnhum.2010.00215
-
Nicotinic acetylcholine receptor subunit α7-knockout mice exhibit degraded auditory temporal processingJournal of Neurophysiology 122:451–465.https://doi.org/10.1152/jn.00170.2019
-
A theory of cortical responsesPhilosophical Transactions of the Royal Society of London. Series B, Biological Sciences 360:815–836.https://doi.org/10.1098/rstb.2005.1622
-
Hierarchical models in the brainPLOS Computational Biology 4:e1000211.https://doi.org/10.1371/journal.pcbi.1000211
-
Cortical circuits for perceptual inferenceNeural Networks 22:1093–1104.https://doi.org/10.1016/j.neunet.2009.07.023
-
Does predictive coding have a future?Nature Neuroscience 21:1019–1021.https://doi.org/10.1038/s41593-018-0200-7
-
Computational psychiatry: from synapses to sentienceMolecular Psychiatry 28:256–268.https://doi.org/10.1038/s41380-022-01743-z
-
Revisiting the cholinergic hypothesis in alzheimer’s disease: Emerging evidence from translational and clinical researchThe Journal Of Prevention of Alzheimer’s Disease 6:1–14.https://doi.org/10.14283/jpad.2018.43
-
Neuromodulation: acetylcholine and memory consolidationTrends in Cognitive Sciences 3:351–359.https://doi.org/10.1016/s1364-6613(99)01365-0
-
Cholinergic modulation of cortical functionJournal of Molecular Neuroscience 30:133–135.https://doi.org/10.1385/JMN:30:1:133
-
Modes and models of forebrain cholinergic neuromodulation of cognitionNeuropsychopharmacology 36:52–73.https://doi.org/10.1038/npp.2010.104
-
Plasticity in the auditory systemHearing Research 362:61–73.https://doi.org/10.1016/j.heares.2017.10.011
-
Attention reverses the effect of prediction in silencing sensory signalsCerebral Cortex 22:2197–2206.https://doi.org/10.1093/cercor/bhr310
-
Cholinergic modulation of the cortical neuronal networkPflugers Archiv 446:17–29.https://doi.org/10.1007/s00424-002-0999-2
-
Stimulus-specific adaptation in the inferior colliculus of the anesthetized ratThe Journal of Neuroscience 29:5483–5493.https://doi.org/10.1523/JNEUROSCI.4153-08.2009
-
BookAuditory systemIn: Paxinos G, editors. The Rat Nervous System. Elsevier - Academic Press. pp. 865–946.
-
BookDeviance Detection and Encoding Acoustic Regularity in the Auditory MidbrainIn: Kandler Karl, editors. The Oxford Handbook of the Auditory Brainstem. Oxford Academic. pp. 707–740.https://doi.org/10.1093/oxfordhb/9780190849061.001.0001
-
Coordinated forms of noradrenergic plasticity in the locus coeruleus and primary auditory cortexNature Neuroscience 18:1483–1492.https://doi.org/10.1038/nn.4090
-
Cholinergic circuitry of the human nucleus basalis and its fate in Alzheimer’s diseaseThe Journal of Comparative Neurology 521:4124–4144.https://doi.org/10.1002/cne.23415
-
Functional connectivity and cholinergic modulation in auditory cortexNeuroscience and Biobehavioral Reviews 35:2058–2063.https://doi.org/10.1016/j.neubiorev.2010.11.010
-
Listening difficulties in children: bottom-up and top-down contributionsJournal of Communication Disorders 45:411–418.https://doi.org/10.1016/j.jcomdis.2012.06.006
-
Free energy, precision and learning: the role of cholinergic neuromodulationThe Journal of Neuroscience 33:8227–8236.https://doi.org/10.1523/JNEUROSCI.4255-12.2013
-
Spatiotemporal specificity in cholinergic control of neocortical functionCurrent Opinion in Neurobiology 26:149–160.https://doi.org/10.1016/j.conb.2014.02.015
-
Distribution of bradykinin B2 receptors in sheep brain and spinal cord visualized by in vitro autoradiographyThe Journal of Comparative Neurology 381:203–218.https://doi.org/10.1002/(sici)1096-9861(19970505)381:2<203::aid-cne7>3.0.co;2-7
-
Stimulus-specific adaptation: can it be a neural correlate of behavioral habituation?The Journal of Neuroscience 31:17811–17820.https://doi.org/10.1523/JNEUROSCI.4790-11.2011
-
Direction-specific adaptation in neuronal and behavioral responses of an insect mechanosensory systemThe Journal of Neuroscience 35:11644–11655.https://doi.org/10.1523/JNEUROSCI.1378-15.2015
-
Uncertainty, epistemics and active inferenceJournal of the Royal Society, Interface 14:20170376.https://doi.org/10.1098/rsif.2017.0376
-
Alzheimer’s disease, hearing loss, and deviance detectionFrontiers in Neuroscience 16:879480.https://doi.org/10.3389/fnins.2022.879480
-
Fast modulation of visual perception by basal forebrain cholinergic neuronsNature Neuroscience 16:1857–1863.https://doi.org/10.1038/nn.3552
-
Cellular mechanisms of brain state-dependent gain modulation in visual cortexNature Neuroscience 16:1331–1339.https://doi.org/10.1038/nn.3464
-
Multiparametric auditory receptive field organization across five cortical fields in the albino ratJournal of Neurophysiology 97:3621–3638.https://doi.org/10.1152/jn.01298.2006
-
Plasticity in the rat posterior auditory field following nucleus basalis stimulationJournal of Neurophysiology 98:253–265.https://doi.org/10.1152/jn.01309.2006
-
Novelty detection in an auditory oddball task on freely moving ratsCommunications Biology 6:1063.https://doi.org/10.1038/s42003-023-05403-y
-
Finding the right control: the mismatch negativity under investigationClinical Neurophysiology 123:507–512.https://doi.org/10.1016/j.clinph.2011.07.035
-
The cognitive neuroscience of sustained attention: where top-down meets bottom-upBrain Research. Brain Research Reviews 35:146–160.https://doi.org/10.1016/s0165-0173(01)00044-3
-
Cholinesterase inhibitors as Alzheimer’s therapeutics (Review)Molecular Medicine Reports 20:1479–1487.https://doi.org/10.3892/mmr.2019.10374
-
Neural elements for predictive codingFrontiers in Psychology 7:1792.https://doi.org/10.3389/fpsyg.2016.01792
-
Reconciling predictive coding and biased competition models of cortical functionFrontiers in Computational Neuroscience 2:4.https://doi.org/10.3389/neuro.10.004.2008
-
A tutorial for information theory in neuroscienceeNeuro 5:ENEURO.0052-18.2018.https://doi.org/10.1523/ENEURO.0052-18.2018
-
Processing of low-probability sounds by cortical neuronsNature Neuroscience 6:391–398.https://doi.org/10.1038/nn1032
-
Multiple time scales of adaptation in auditory cortex neuronsThe Journal of Neuroscience 24:10440–10453.https://doi.org/10.1523/JNEUROSCI.1905-04.2004
-
Correlating stimulus-specific adaptation of cortical neurons and local field potentials in the awake ratThe Journal of Neuroscience 29:13837–13849.https://doi.org/10.1523/JNEUROSCI.3475-09.2009
-
Deconstructing the architecture of dorsal and ventral attention systems with dynamic causal modelingThe Journal of Neuroscience 32:10637–10648.https://doi.org/10.1523/JNEUROSCI.0414-12.2012
-
Specific long-term memory traces in primary auditory cortexNature Reviews. Neuroscience 5:279–290.https://doi.org/10.1038/nrn1366
-
Context-dependent inhibitory control of stimulus-specific adaptationThe Journal of Neuroscience 42:4629–4651.https://doi.org/10.1523/JNEUROSCI.0988-21.2022
-
Precision and the Bayesian brainCurrent Biology 31:R1026–R1032.https://doi.org/10.1016/j.cub.2021.07.044
-
Acetylcholine in cortical inferenceNeural Networks 15:719–730.https://doi.org/10.1016/s0893-6080(02)00058-8
Article and author information
Author details
Funding
Agencia Estatal de Investigación (PID2019-104570RB-I00)
- Manuel S Malmierca
- Ana Belén Lao-Rodríguez
Fundación Ramón Areces (CIVP20A6616)
- Manuel S Malmierca
- David Pérez-González
Horizon 2020 Framework Programme (No 952378-BrainTwin)
- Manuel S Malmierca
- Ana Belén Lao-Rodríguez
- David Pérez-González
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank Drs. Francisco Garcia-Rosales and Julio Hechavarria for their advice in the mutual information analysis and Drs. Ryszard Auksztulewicz, Yaneri A Ayala, Donald Caspary, Diego Contreras, Karl Friston, Lauren Harms, and Marina Picciotto for their constructive criticisms of a previous version.
Ethics
All experimental procedures were carried out at the University of Salamanca and all methodological procedures were approved by the Bioethics Committee for Animal Care of the University of Salamanca (USAL-ID-574 and 683) and performed in compliance with the standards of the European Convention ETS 123, the European Union Directive 2010/63/EU, and the Spanish Royal Decree 53/2013 for the use of animals in scientific research.
Version history
- Sent for peer review:
- Preprint posted:
- Reviewed Preprint version 1:
- Reviewed Preprint version 2:
- Version of Record published:
- Version of Record updated:
Cite all versions
You can cite all versions using the DOI https://doi.org/10.7554/eLife.91475. This DOI represents all versions, and will always resolve to the latest one.
Copyright
© 2023, Pérez-González, Lao-Rodríguez et al.
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
-
- 1,542
- views
-
- 148
- downloads
-
- 4
- 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
-
- Neuroscience
Each sensory modality has its own primary and secondary thalamic nuclei. While the primary thalamic nuclei are well understood to relay sensory information from the periphery to the cortex, the role of secondary sensory nuclei is elusive. We trained head-fixed mice to attend to one sensory modality while ignoring a second modality, namely to attend to touch and ignore vision, or vice versa. Arrays were used to record simultaneously from the secondary somatosensory thalamus (POm) and secondary visual thalamus (LP). In mice trained to respond to tactile stimuli and ignore visual stimuli, POm was robustly activated by touch and largely unresponsive to visual stimuli. A different pattern was observed when mice were trained to respond to visual stimuli and ignore touch, with POm now more robustly activated during visual trials. This POm activity was not explained by differences in movements (i.e. whisking, licking, pupil dilation) resulting from the two tasks. Post hoc histological reconstruction of array tracks through POm revealed that subregions varied in their degree of plasticity. LP exhibited similar phenomena. We conclude that behavioral training reshapes activity in secondary thalamic nuclei. Secondary nuclei respond to the same behaviorally relevant, reward-predicting stimuli regardless of stimulus modality.
-
- Neuroscience
The relation between neural activity and behaviorally relevant variables is at the heart of neuroscience research. When strong, this relation is termed a neural representation. There is increasing evidence, however, for partial dissociations between activity in an area and relevant external variables. While many explanations have been proposed, a theoretical framework for the relationship between external and internal variables is lacking. Here, we utilize recurrent neural networks (RNNs) to explore the question of when and how neural dynamics and the network’s output are related from a geometrical point of view. We find that training RNNs can lead to two dynamical regimes: dynamics can either be aligned with the directions that generate output variables, or oblique to them. We show that the choice of readout weight magnitude before training can serve as a control knob between the regimes, similar to recent findings in feedforward networks. These regimes are functionally distinct. Oblique networks are more heterogeneous and suppress noise in their output directions. They are furthermore more robust to perturbations along the output directions. Crucially, the oblique regime is specific to recurrent (but not feedforward) networks, arising from dynamical stability considerations. Finally, we show that tendencies toward the aligned or the oblique regime can be dissociated in neural recordings. Altogether, our results open a new perspective for interpreting neural activity by relating network dynamics and their output.