CoChat
Sign Up
Cited by 344Neuroscience and Music PerceptionHearing, Cochlea, Tinnitus, GeneticsEmbodied and Extended Cognition

The DOI matches its Crossref record; 75 of the 75 references it cites were found in a registry or an index.

Neurons along the auditory pathway exhibit a hierarchical organization of prediction error

Open access
Parras, Gloria G., Nieto-Diego, Javier, Carbajal, Guillermo V., Valdés-Baizabal, Catalina, Escera, Carles, Malmierca, Manuel S.·2017·Nature Communications
DOI 10.1038/s41467-017-02038-6PMC PMC5732270
Open this paper in CoChat — free accountOpen PDF

Abstract

Perception is characterized by a reciprocal exchange of predictions and prediction error signals between neural regions. However, the relationship between such sensory mismatch responses and hierarchical predictive processing has not yet been demonstrated at the neuronal level in the auditory pathway. We recorded single-neuron activity from different auditory centers in anaesthetized rats and awake mice while animals were played a sequence of sounds, designed to separate the responses due to prediction error from those due to adaptation effects. Here we report that prediction error is organized hierarchically along the central auditory pathway. These prediction error signals are detectable in subcortical regions and increase as the signals move towards auditory cortex, which in turn demonstrates a large-scale mismatch potential. Finally, the predictive activity of single auditory neurons underlies automatic deviance detection at subcortical levels of processing. These results demonstrate that prediction error is a fundamental component of singly auditory neuron responses.

Overview

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Integrity signals

Registry

Crossref · 10.1038/s41467-017-02038-6

Cites

75 of 75 references verified0 citation errors · 0 not checkable

Statistics

Figures

Fig. 1 Experimental design. a Sketch of experimental setup. While stimulating with sequences of pure tones, isolated neurons were recorded from three auditory nuclei of anesthetized rats: IC, MGB and AC (colored). The schematic representation of the ascending auditory pathway information flow (orange) shows how lemniscal (green) and non-lemniscal (purple) subdivisions can be distinguished in the IC, MGB and AC. b Stimulation sequences. For each neuron, 10 tones of evenly-spaced frequencies were selected to construct the stimulation sequences. Each tone fi (i = 1…10) lying inside the neuron’s receptive field could be presented in two experimental conditions (deviant and standard, in separated oddball sequences, left column), and two control conditions (cascade and many-standards, right column) for adaptation effects. Note that ascending and descending deviant tones will be compared to the control ascending or control descending sequences, respectively. They will also be compared to the many-standards sequence for both types of deviants (see Methods). c Decomposition of neuronal mismatch responses (DEV−STD) to the oddball sequence using either one of the control conditions. Under the assumption of predictive coding, CTR−STD (if positive) represents repetition suppression, and DEV−CTR (if positive) represents prediction error. d Two hypothetical scenarios according to two possible competing mechanisms accounting for the neuronal mismatch: SSA (top) and predictive coding (bottom). For SSA, there is response suppression to the standard (blue bars), which progressively increases from lower order to higher order. In addition, due to suppression of the deviant relative to control, the prediction error is increasingly negative (blue and orange bars) as one progresses to higher-order regions. For predictive coding, repetition suppression of the standard (blue bars) increases from lower to higher-order regions. Unlike SSA, responses to deviants are higher than controls, especially in higher-order regions, leading to a positive prediction error (orange bars)
Fig. 1 Experimental design. a Sketch of experimental setup. While stimulating with sequences of pure tones, isolated neurons were recorded from three auditory nuclei of anesthetized rats: IC, MGB and AC (colored). The schematic representation of the ascending auditory pathway information flow (orange) shows how lemniscal (green) and non-lemniscal (purple) subdivisions can be distinguished in the IC, MGB and AC. b Stimulation sequences. For each neuron, 10 tones of evenly-spaced frequencies were selected to construct the stimulation sequences. Each tone fi (i = 1…10) lying inside the neuron’s receptive field could be presented in two experimental conditions (deviant and standard, in separated oddball sequences, left column), and two control conditions (cascade and many-standards, right column) for adaptation effects. Note that ascending and descending deviant tones will be compared to the control ascending or control descending sequences, respectively. They will also be compared to the many-standards sequence for both types of deviants (see Methods). c Decomposition of neuronal mismatch responses (DEV−STD) to the oddball sequence using either one of the control conditions. Under the assumption of predictive coding, CTR−STD (if positive) represents repetition suppression, and DEV−CTR (if positive) represents prediction error. d Two hypothetical scenarios according to two possible competing mechanisms accounting for the neuronal mismatch: SSA (top) and predictive coding (bottom). For SSA, there is response suppression to the standard (blue bars), which progressively increases from lower order to higher order. In addition, due to suppression of the deviant relative to control, the prediction error is increasingly negative (blue and orange bars) as one progresses to higher-order regions. For predictive coding, repetition suppression of the standard (blue bars) increases from lower to higher-order regions. Unlike SSA, responses to deviants are higher than controls, especially in higher-order regions, leading to a positive prediction error (orange bars)CC BY · PubMed Central
Fig. 2 Prediction error in representative examples of neuronal responses in anaesthetized rat. a Examples of lemniscal neuronal responses in each recorded auditory station (columns). The first row contains schematics of the lemniscal subdivisions (green) within each nuclei. The second row shows the frequency-response area (representation of neuronal sensitivity to different frequency-intensity combinations) of representative lemniscal neurons from each nucleus. Ten grey dots within each frequency-response area represent the ten tones (fi) selected to build the experimental sequences (see Methods). The third row displays the measured responses of the particular neuron to each fi tone (baseline-corrected spike counts, averaged within 0–180 ms after tone onset) for all conditions tested. Note that measured conditions tend to overlap in the subcortical stations (ICL and MGBL), and only start differentiating from each other once auditory information reaches the cortex (ACL). The fourth row contains sample peri-stimulus histograms comparing the neuronal responses to each condition tested for an indicated fi tone. A thick horizontal line represents stimulus duration. A small inset within the upper right corner of each panel features the isolated spike (mean ± SEM) of that single neuron. b Examples of non-lemniscal neuronal responses in each recorded auditory nuclei, organized as in a. The first row highlights non-lemniscal divisions in purple. In the second row, note frequency-response areas tend to be more broadly tuned, as compared to lemniscal neurons. In the third row, responses to deviant conditions tend to relatively increase and distance themselves from their corresponding controls as information ascends in the auditory pathway. Also note that responses to last standards are feeble or even completely missing across all non-lemniscal stations (ICNL, MGBNL and ACNL). In the last row, the strong influence of the experimental condition over the neuronal response to the same tone can be clearly appreciated in the three nuclei
Fig. 2 Prediction error in representative examples of neuronal responses in anaesthetized rat. a Examples of lemniscal neuronal responses in each recorded auditory station (columns). The first row contains schematics of the lemniscal subdivisions (green) within each nuclei. The second row shows the frequency-response area (representation of neuronal sensitivity to different frequency-intensity combinations) of representative lemniscal neurons from each nucleus. Ten grey dots within each frequency-response area represent the ten tones (fi) selected to build the experimental sequences (see Methods). The third row displays the measured responses of the particular neuron to each fi tone (baseline-corrected spike counts, averaged within 0–180 ms after tone onset) for all conditions tested. Note that measured conditions tend to overlap in the subcortical stations (ICL and MGBL), and only start differentiating from each other once auditory information reaches the cortex (ACL). The fourth row contains sample peri-stimulus histograms comparing the neuronal responses to each condition tested for an indicated fi tone. A thick horizontal line represents stimulus duration. A small inset within the upper right corner of each panel features the isolated spike (mean ± SEM) of that single neuron. b Examples of non-lemniscal neuronal responses in each recorded auditory nuclei, organized as in a. The first row highlights non-lemniscal divisions in purple. In the second row, note frequency-response areas tend to be more broadly tuned, as compared to lemniscal neurons. In the third row, responses to deviant conditions tend to relatively increase and distance themselves from their corresponding controls as information ascends in the auditory pathway. Also note that responses to last standards are feeble or even completely missing across all non-lemniscal stations (ICNL, MGBNL and ACNL). In the last row, the strong influence of the experimental condition over the neuronal response to the same tone can be clearly appreciated in the three nucleiCC BY · PubMed Central
Fig. 3 Prediction error at population level for each station in anaesthetized rat. Distribution of normalized responses and related indices of neuronal mismatch (iMM), repetition suppression (iRS) and prediction error (iPE). a Each grey dot in these scatter plots represents the three normalized responses of a single neuron to the same tone played as deviant (DEV), as standard (STD) and as control (CTR). Indexes result from the difference between two of these normalized responses, represented in the axes surrounding the scatter plots, where the dotted black lines marks the absence of difference between conditions (index = 0). Solid black lines represent the mean of each index, corresponding their intersection to the center of gravity of the distribution of responses in the normalized space. Note how, while the intersection for lemniscal subcortical stations (ICL and MGBL) is skewed towards CTR, in their non-lemniscal counterparts (ICNL and MGBNL) as well as all over the cortex (ACL and ACNL) the center of gravity of the distribution shifts closer and closer to DEV as it moves up in the auditory pathway, increasing the iPE as auditory information reaches higher-order stations. b Histograms represent distributions within stations of the three indexes for each neuronal response. Solid black lines indicate medians. The noticeable overall tendency of the median indexes to shift towards more positive values, from IC through MGB to AC, and from lemniscal to non-lemniscal divisions, unveils a hierarchy of processing in the auditory pathway
Fig. 3 Prediction error at population level for each station in anaesthetized rat. Distribution of normalized responses and related indices of neuronal mismatch (iMM), repetition suppression (iRS) and prediction error (iPE). a Each grey dot in these scatter plots represents the three normalized responses of a single neuron to the same tone played as deviant (DEV), as standard (STD) and as control (CTR). Indexes result from the difference between two of these normalized responses, represented in the axes surrounding the scatter plots, where the dotted black lines marks the absence of difference between conditions (index = 0). Solid black lines represent the mean of each index, corresponding their intersection to the center of gravity of the distribution of responses in the normalized space. Note how, while the intersection for lemniscal subcortical stations (ICL and MGBL) is skewed towards CTR, in their non-lemniscal counterparts (ICNL and MGBNL) as well as all over the cortex (ACL and ACNL) the center of gravity of the distribution shifts closer and closer to DEV as it moves up in the auditory pathway, increasing the iPE as auditory information reaches higher-order stations. b Histograms represent distributions within stations of the three indexes for each neuronal response. Solid black lines indicate medians. The noticeable overall tendency of the median indexes to shift towards more positive values, from IC through MGB to AC, and from lemniscal to non-lemniscal divisions, unveils a hierarchy of processing in the auditory pathwayCC BY · PubMed Central
Fig. 4 Emergence of prediction error along the auditory hierarchy. a Median normalized responses (lines indicate SEM) to the deviant, standard and control within each station. b Median indices of prediction error (orange) and repetition suppression (cyan), represented with respect to the baseline set by the control. Thereby, iPE is upwards-positive while iRS is downwards-positive. Each median index corresponds to differences between normalized responses in a. Asterisks denote statistical significance of iPE against zero median (*p = 0.05, **p = 0.01, ***p = 0.001, see Table 1). c Linear model fitted for the iPE using SPL and direction (ascending / descending) as predictors. Error bars denote mean and SEM for each SPL and direction. Note the model predicts greater iPE values for ascending conditions at low intensities. d Same as in b, but only representing ascending conditions at intensities equal or lower than 40 dB SPL
Fig. 4 Emergence of prediction error along the auditory hierarchy. a Median normalized responses (lines indicate SEM) to the deviant, standard and control within each station. b Median indices of prediction error (orange) and repetition suppression (cyan), represented with respect to the baseline set by the control. Thereby, iPE is upwards-positive while iRS is downwards-positive. Each median index corresponds to differences between normalized responses in a. Asterisks denote statistical significance of iPE against zero median (*p = 0.05, **p = 0.01, ***p = 0.001, see Table 1). c Linear model fitted for the iPE using SPL and direction (ascending / descending) as predictors. Error bars denote mean and SEM for each SPL and direction. Note the model predicts greater iPE values for ascending conditions at low intensities. d Same as in b, but only representing ascending conditions at intensities equal or lower than 40 dB SPLCC BY · PubMed Central
Fig. 5 Correlation of iPE and prediction error potential (PE-LFP). a Population grand-averages for different response measures, computed for each lemniscal station (in columns, represented in first row highlighted in green). The second row shows the average LFP across all tested tones and single neurons from each station for different conditions. The third row displays the average firing rate profiles for each station as normalized spike-density functions. The fourth row contains the prediction error potentials (PE-LFP, black trace), which is the difference wave of the deviant and the control LFP. Along PE-LFP, the time course of the average iPE is plotted in orange (mean ± SEM, asterisks indicating significant iPE for the corresponding time window; Wilcoxon signed-rank test for 12 comparisons, corrected for FDR = 0.1). Next row shows an instantaneous p value (white trace) of the corresponding PE-LFP (paired t test against equal means, corrected for FDR = 0.1, critical threshold for significance set at 0.05 represented as a horizontal bar). Thick black bars of the grey panel mark time intervals for which the average PE-LFP is significant. Note that only ACL shows a significant prediction error signal. b Same as in a but computed for each non-lemniscal station (highlighted in purple in the first row). Note in the last row significant PE-LFPs appear in all three stations (ICNL, MGBNL and ACNL), and prominently in ACNL. Note also how highest iPE values are concurrent with the strongest PE-LFPs in time and location (auditory cortex, both ACL and ACNL)
Fig. 5 Correlation of iPE and prediction error potential (PE-LFP). a Population grand-averages for different response measures, computed for each lemniscal station (in columns, represented in first row highlighted in green). The second row shows the average LFP across all tested tones and single neurons from each station for different conditions. The third row displays the average firing rate profiles for each station as normalized spike-density functions. The fourth row contains the prediction error potentials (PE-LFP, black trace), which is the difference wave of the deviant and the control LFP. Along PE-LFP, the time course of the average iPE is plotted in orange (mean ± SEM, asterisks indicating significant iPE for the corresponding time window; Wilcoxon signed-rank test for 12 comparisons, corrected for FDR = 0.1). Next row shows an instantaneous p value (white trace) of the corresponding PE-LFP (paired t test against equal means, corrected for FDR = 0.1, critical threshold for significance set at 0.05 represented as a horizontal bar). Thick black bars of the grey panel mark time intervals for which the average PE-LFP is significant. Note that only ACL shows a significant prediction error signal. b Same as in a but computed for each non-lemniscal station (highlighted in purple in the first row). Note in the last row significant PE-LFPs appear in all three stations (ICNL, MGBNL and ACNL), and prominently in ACNL. Note also how highest iPE values are concurrent with the strongest PE-LFPs in time and location (auditory cortex, both ACL and ACNL)CC BY · PubMed Central
Fig. 6 Prediction error in representative examples of neuronal responses in awake mouse. a Examples of lemniscal multiunit activity recorded in two auditory nuclei (columns). The first row contains schematics of the lemniscal subdivisions (green) within each station. The second row shows a frequency-response area of each nuclei. Ten grey dots within those frequency-response area represent the ten tones (fi) selected to build the experimental sequences (see Methods). The third row displays the measured responses to each fi tone (baseline-corrected spike counts, averaged within 0–180 ms after tone onset) for all conditions tested. The fourth row contains sample peri-stimulus histagrams comparing the neuronal responses to each condition tested for an indicated fi tone. Stimulus duration is represented by a thick horizontal line. b Examples of multiunit activity recorded in non-lemniscal divisions (first row, colored purple) of each auditory nuclei, organized as in a
Fig. 6 Prediction error in representative examples of neuronal responses in awake mouse. a Examples of lemniscal multiunit activity recorded in two auditory nuclei (columns). The first row contains schematics of the lemniscal subdivisions (green) within each station. The second row shows a frequency-response area of each nuclei. Ten grey dots within those frequency-response area represent the ten tones (fi) selected to build the experimental sequences (see Methods). The third row displays the measured responses to each fi tone (baseline-corrected spike counts, averaged within 0–180 ms after tone onset) for all conditions tested. The fourth row contains sample peri-stimulus histagrams comparing the neuronal responses to each condition tested for an indicated fi tone. Stimulus duration is represented by a thick horizontal line. b Examples of multiunit activity recorded in non-lemniscal divisions (first row, colored purple) of each auditory nuclei, organized as in aCC BY · PubMed Central
Fig. 7 Population results for awake mouse. a Summary of population results for each recorded station (columns). The first row displays lemniscal (green) and non-lemniscal (purple) subdivisions of two recorded auditory nuclei of the mouse brain. The second row contains scatter plots featuring normalized responses of each multiunit recording to the same tone played as DEV, STD and CTR (grey dots) and the mean population values of each index (solid black bars). The third row contains the average LFP across all tested tones from each station for different conditions. Thick black bars at the bottom of the panels mark the time intervals were the difference between the deviant LFP and the control LFP is significant, thereby producing a prediction error potential. b Median normalized responses (bar indicate interquartile range) to the deviant, standard and control within each station. c Median indices of prediction error (orange) and repetition suppression (cyan), represented with respect to the baseline set by the control. Asterisks denote statistical significance of iPE against zero median (*p = 0.05, **p = 0.01, ***p = 0.001). Note the overall similarities with results in the anaesthetized rat (Figs. 3–5), confirming a hierarchical generation of prediction error also in awake preparations
Fig. 7 Population results for awake mouse. a Summary of population results for each recorded station (columns). The first row displays lemniscal (green) and non-lemniscal (purple) subdivisions of two recorded auditory nuclei of the mouse brain. The second row contains scatter plots featuring normalized responses of each multiunit recording to the same tone played as DEV, STD and CTR (grey dots) and the mean population values of each index (solid black bars). The third row contains the average LFP across all tested tones from each station for different conditions. Thick black bars at the bottom of the panels mark the time intervals were the difference between the deviant LFP and the control LFP is significant, thereby producing a prediction error potential. b Median normalized responses (bar indicate interquartile range) to the deviant, standard and control within each station. c Median indices of prediction error (orange) and repetition suppression (cyan), represented with respect to the baseline set by the control. Asterisks denote statistical significance of iPE against zero median (*p = 0.05, **p = 0.01, ***p = 0.001). Note the overall similarities with results in the anaesthetized rat (Figs. 3–5), confirming a hierarchical generation of prediction error also in awake preparationsCC BY · PubMed Central

Verification

Claims the paper makes

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Claims other papers make about this one

Each is checked against this paper’s own text: does it say what the citing paper says it does?

Self-check

Percentages, sums and reported differences re-computed from the paper’s own numbers.

The self-check needs the paper’s full text.

Citations

Cites

Show all 75 references

Cited by

Sorted by citations

Items 301–341 of 341

Related

See all

Metadata

Authors
Parras, G. G., Nieto-Diego, J., Carbajal, G. V., Valdés-Baizabal, C., Escera, C., & Malmierca, M. S.
Affiliations
  • Instituto de Investigación Biomédica de Salamanca (ES)
  • Institut de Recerca Sant Joan de Déu (ES)
  • Universidad de Salamanca (ES)
Corresponding author
Malmierca, Manuel S.
ORCID
Venue
Nature Communications
Publisher
Springer Science and Business Media LLC
Published
2017-12-15
Volume
8
Issue
1
Pages
2148-2148
ISSN
2041-1723 (ISSN-L)
DOI
10.1038/s41467-017-02038-6
PMC
PMC5732270
OpenAlex
W2774201964
Cited by
344
References
75
Citations by year
2026: 38 · 2025: 33 · 2024: 55 · 2023: 60 · 2022: 39 · 2021: 29 · 2020: 50 · 2019: 25 · 2018: 14
FWCI
13.96
Type
journal article
Language
en
Open access
Licence
Locations
  • Nature Communications (publishedVersion, cc-by)
  • PubMed Central (publishedVersion, cc-by)
  • PubMed (publishedVersion)
  • Gredos (University of Salamanca) (publishedVersion, cc-by)
  • Hispana (publishedVersion, cc-by)
  • RePEc: Research Papers in Economics (submittedVersion)
  • Dipòsit Digital de la Universitat de Barcelona (Universitat de Barcelona) (submittedVersion, cc-by)
  • DOAJ (DOAJ: Directory of Open Access Journals) (submittedVersion)
  • DOAJ (DOAJ: Directory of Open Access Journals) (submittedVersion, cc-by-sa)
Funders
  • Generalitat de Catalunya (SGR 2014, SGR2014-177)
  • Junta de Castilla y León
  • European Social Fund (2007–2013)
Topics
  • Life Sciences › Neuroscience › Cognitive Neuroscience › Neuroscience and Music Perception
  • Life Sciences › Neuroscience › Sensory Systems › Hearing, Cochlea, Tinnitus, Genetics
  • Life Sciences › Neuroscience › Cognitive Neuroscience › Embodied and Extended Cognition
Keywords
prediction error, auditory pathway, single-unit activity, auditory cortex, subcortical auditory processing, mismatch response, deviance detection, sensory prediction, auditory adaptation

About these checks

Registry, retraction, Cites and statistics results are computed, not judged, with one exception: a printed reference that no index lists is read once by a model before it counts as a citation error. Claims are extracted by a model and kept only when their quote appears verbatim in the text; each cited source is then checked, with a model’s help, for whether it says what the claim relies on. Those verdicts can be wrong — read the claim’s quote and the cited source before relying on one.

Statistics and arithmetic are recomputed, not judged.

Confidence bands are shown only for fields calibrated on our annotated test set. None is calibrated yet, so no band is shown.

What counts as a claim ›How we measure confidence ›How we check this ›

Are you an author, or do you see a problem with this page? Request a review