Respiration as a dynamic modulator of sensory sampling.
The 5 matches
- [1] § Methods › Respiration preprocessing and phase extraction ↔ Code/fig1d.m, lines 13 to 64 · score 0.65 · peak prominence, breathing rates, findpeaks, troughs, respiration · Report a mismatch
- [2] § Methods › MEG preprocessing and segmentation ↔ Code/fig1d.m, lines 66 to 156 · score 0.65 · ft_preprocessing, outliers, eye, segmented, artefacts, channels · Report a mismatch
- [3] § Results › Distinct oscillatory signatures of behaviourally relevant stimulus characteristics ↔ Code/fig2.m, lines 454 to 552 · score 0.56 · Individual alpha power, posterior alpha, BF10, Violin, uncued, Figure 2 · Report a mismatch
- [4] § Methods › Respiration phase-dependent fitting of the psychometric function ↔ Code/fig1bc.m, lines 267 to 349 · score 0.56 · cumulative Gaussian, Psignifit, sigmoid, fitted, refitting, threshold · Report a mismatch
- [5] § Results › Respiration-driven spectral changes within the RMBO network mediate perception ↔ Code/fig3bcd.m, lines 100 to 205 · score 0.54 · 1–40 Hz, power spectra, respiration phase, FDR, ROIs, onset · Report a mismatch
Paper
Loaded from Europe PMC by your browser; OSCR does not store it: doi.org · Europe PMC
The paper is loaded when this pane is shown.
The authors' code
MATLAB · 166 lines · 5.6 KB · no license · 2 matches
- clearvars
- addpath(genpath('D:\RESCUE\'));
- addpath('C:\Program Files\MATLAB\Toolboxes\fieldtrip-20240111\fieldtrip-20240111\');
- ft_defaults;
- ft_warning off;
- datapath = 'D:\RESCUE\_data\';
- outpath = 'D:\RESCUE\_out\';
- ids = [2549 3051 3066 3048 3033 3050 3067 3064 3062 3063 2521 2550 2448 2314 2168 ...
- 2524 2480 2299 2749 1566 2214 2348 1822 3090 3047 3081 3053 2330 3078 3105]; % full
- %% extract respiration frequencies for rest and task
- conds = {'nasal','oral'};
- respfreqs = zeros(3,numel(ids));
- for isub = 1:numel(ids)
- % rest
- for icond = 1:2
- load([outpath 'phasevec_' num2str(ids(isub)) '_rest_' conds{icond} '.mat'],'phasevec');
- tmp = phasevec(~isnan(phasevec));
- [~, peaks] = findpeaks(tmp, 'MinPeakProminence',1);
- [~, troughs] = findpeaks(-tmp, 'MinPeakProminence',1);
- for k = 2:numel(peaks)
- diff(k-1) = peaks(k)-peaks(k-1);
- end
- respfreqs(icond,isub) = 1/(mean(diff)/600); % breathing rate (Hz)
- end
- % task
- tmpfreqs = [];
- for irun = 1:6
- load([outpath 'phasevec_' num2str(ids(isub)) '_task_' num2str(irun) '.mat'],'phasevec');
- tmp = phasevec(~isnan(phasevec));
- [~, peaks] = findpeaks(tmp, 'MinPeakProminence',1);
- [~, troughs] = findpeaks(-tmp, 'MinPeakProminence',1);
- for k = 2:numel(peaks)
- diff(k-1) = peaks(k)-peaks(k-1);
- end
- tmpfreqs(irun) = 1/(mean(diff)/600); % breathing rate (Hz)
- end
- respfreqs(3,isub) = mean(tmpfreqs);
- end
- % plot frequencies for rest (nasal, oral) and task
- addpath('D:\RespSpeech');
- cmap = {'#57315A','#8C2B59','#8C8474'};
- figure;
- v = violinplot(respfreqs');
- for k = 1:3
- v(1,k).ViolinColor = hex2rgb(cmap{k});
- v(1,k).ScatterPlot.SizeData = 25;
- v(1,k).ScatterPlot.MarkerFaceAlpha = 1;
- v(1,k).ScatterPlot.MarkerEdgeColor = [1 1 1];
- end
- ylim([0.1 0.5]);
- xlabel('condition');
- ylabel('frequency (Hz)');
- set(gca,'XTickLabel',{'nasal (rest)','oral (rest)','task mean'});
- set(gca,'YTick',[0.1:0.1:0.5]);
- title('respiration frequencies');
- % test
- [p,h,stats] = signrank(respfreqs(1,:),respfreqs(2,:)) % nasal > oral *
- [p,h,stats] = signrank(respfreqs(1,:),respfreqs(3,:)) % nasal < task **
- %% NEW 251112 - R1 Analysis for saccadic patterns across conditions
- ids = [2549 3051 3066 3048 3033 3050 3067 3064 3062 3063 2521 2550 2448 2314 2168 ...
- 2524 2480 2299 2749 1566 2214 2348 1822 3090 3047 3081 3053 2330 3078 3105];
- % create one full data set of segmented horizontal eye movements for each participant
- for isub = 1:numel(ids)
- trials = [];
- for irun = 1:6
- % load raw data sampled at 600 Hz
- idfiles = dir(fullfile([datapath 'A' num2str(ids(isub))], '*.ds'));
- fn = [datapath 'A' num2str(ids(isub)) '\' idfiles(irun+1).name];
- cfg = [];
- cfg.dataset = fn;
- cfg.continuous = 'yes';
- raw = ft_preprocessing(cfg);
- % crop to pupil
- rawpup = ft_selectdata(struct('channel',{{'UADC005'}}),raw);
- % actual trial definition
- cfg = [];
- cfg.dataset = fn;
- cfg.trialfun = 'dk_trialfun_rescue';
- cfg.trialdef.eventtype = 'frontpanel trigger';
- cfg.trialdef.eventvalue = [16 17 18 19 21 22 23 24];
- cfg.trialdef.prestim = 1.85; %1.25;
- cfg.trialdef.poststim = 0.65; %0.75;
- cfg_stimonset = ft_definetrial(cfg);
- trials = [trials; cfg_stimonset.trl(11:end,:)]; % crop first ten trials
- segrawpup = ft_redefinetrial(cfg_stimonset, rawpup);
- all{irun} = segrawpup;
- end
- cfg = [];
- cfg.keepsampleinfo = 'no';
- fullsegrawpup = ft_appenddata(cfg,all{1},all{2},all{3},all{4},all{5},all{6});
- save([outpath 'fullsegrawpup_' num2str(ids(isub)) '.mat'], 'fullsegrawpup');
- end
- % determine the number of saccades within each category
- saccnum = zeros(numel(ids),4);
- for isub = 1:numel(ids)
- isub
- load([outpath 'fullsegrawpup_' num2str(ids(isub)) '.mat'], 'fullsegrawpup');
- % crop to prestimulus time window
- cfg = [];
- cfg.latency = [-1.6 0];
- prestimrawpup = ft_selectdata(cfg, fullsegrawpup);
- % trial selection
- trlinfo = prestimrawpup.trialinfo(:,1);
- uu = find(trlinfo == 16 | trlinfo == 24); % uncued, untimed
- ut = find(trlinfo == 18 | trlinfo == 22); % uncued, timed
- cu = find(trlinfo == 17 | trlinfo == 23); % cued, untimed
- ct = find(trlinfo == 19 | trlinfo == 21); % cued, timed
- cfg = [];
- cfg.trials = uu;
- uupup = ft_selectdata(cfg, prestimrawpup);
- cfg.trials = ut;
- utpup = ft_selectdata(cfg, prestimrawpup);
- cfg.trials = cu;
- cupup = ft_selectdata(cfg, prestimrawpup);
- cfg.trials = ct;
- ctpup = ft_selectdata(cfg, prestimrawpup);
- % run the artefact detection
- cfg = [];
- cfg.artfctdef.zvalue.channel = 'all';
- cfg.artfctdef.zvalue.cutoff = 5;
- uusacc = ft_artifact_zvalue(cfg, uupup);
- utsacc = ft_artifact_zvalue(cfg, utpup);
- cusacc = ft_artifact_zvalue(cfg, cupup);
- ctsacc = ft_artifact_zvalue(cfg, ctpup);
- % save number of saccades
- saccnum(isub,1) = length(uusacc.artfctdef.zvalue.artifact);
- saccnum(isub,2) = length(utsacc.artfctdef.zvalue.artifact);
- saccnum(isub,3) = length(cusacc.artfctdef.zvalue.artifact);
- saccnum(isub,4) = length(ctsacc.artfctdef.zvalue.artifact);
- end
- % remove outlier
- saccnum(19,:) = [];
- writematrix(saccnum,[outpath 'saccnum.txt']);
- median(saccnum)
fig1d.m, no license · at the source
Overview
- Institute for Biomagnetism and Biosignal Analysis, University of Münster,Münster, Germany
- Otto Creutzfeldt Center for Cognitive and Behavioral Neuroscience, University of Münster,Münster, Germany
Abstract
Respiration dynamically modulates sensory perception by orchestrating transient states of the brain and the body. Using simultaneous recordings of high-density magneto- encephalography (MEG), respiration, and pupillometry, we show that human perceptual sensitivity to near-threshold visual stimuli was enhanced during inspiration, coinciding with respiration-modulated increases in arousal neuromodulation and cortical excitability. Participants adapted their breathing patterns to align with predictable stimulus onset, and this adaptive respiratory alignment correlated with improved performance. We further reveal that respiration-modulated changes in alpha and beta oscillations reflect distinct shifts in sensory and motor excitability, respectively. Crucially, respiration-resolved multivariate Granger causality analyses demonstrate that the breathing rhythm systematically shapes directed information flow within a widespread interoceptive network. This respiration-brain coupling was flexibly adjusted based on stimulus predictability, highlighting a mechanism for active sensing which integrates internal bodily rhythms with external sensory demands to optimise perception.
Reproduced under the paper's license (CC BY), from the paper cited above.
Repository
Its files are read in the Code ↔ Paper reader above, with 5 matches between paragraphs and lines of code.
OSF qasvp
Availability: 1 check, the latest on 29 September 2026: the link answers (HTTP 200)
- 29 September 2026: the link answers (HTTP 200)
8 files
- Code/
fig1bc.m , MATLAB, 349 lines, 1 match - Code/
fig1d.m , MATLAB, 166 lines, 2 matches - Code/
fig1ef.m , MATLAB, 144 lines - Code/
fig1gh.m , MATLAB, 557 lines - Code/
fig2.m , MATLAB, 658 lines, 1 match - Code/
fig3bcd.m , MATLAB, 357 lines, 1 match - Code/
fig3fg.m , MATLAB, 216 lines - Code/
fig4.m , MATLAB, 569 lines
Code availability
Software code is publicly available via osf.io/
Reproduced under the paper's license (CC BY), from the paper cited above.
Tracing map
Proposed by the machine: these links were found in the paper and verified at the source, without human review. The map will receive a Zenodo DOI once one of the paper's authors has validated it with their ORCID.
What the map holds:
- 1 repository of the authors' code, each at its verified commit, with its license and how the link was found in the paper;
- 8 scripts, each with its path and the digest of its content;
- 5 matches between paragraphs of the paper and lines of the code (method lexical-v1);
- neither the text of the paper nor the code itself.
Its JSON (tracing-map.json) is deposited on Zenodo with its DOI once the map is validated.
Data
No dataset and no data link were found in the paper.
Data availability
Raw data are protected by data privacy laws and can be made available by means of a formal data sharing agreement. The processed and anonymised data MEG, physiological, and behavioural data are publicly available at the Open Science Framework (OSF) via osf.io/
Reproduced under the paper's license (CC BY), from the paper cited above.
Versions
The history of this record: each version stored by the harvester or made by a correction of its authors or of the maintainers of its code, and what changed in its facts. The texts of the paper (its abstract, its availability statements) are not part of it; versions that changed only those are not listed.
Version 1, 29 September 2026: the first record
Recorded: type, language, journal, volume, issue, pages, dates, 6 authors, 2 keywords, 12 MeSH terms, 3 funders, 70 references.
Cite
This paper
Chalas, N., Saltafossi, M., Berther, T., Balestrieri, E., Abbasi, O., & Kluger, D. S. (2026). Respiration as a dynamic modulator of sensory sampling. Nature communications, 17(1), 3261. https://
BibTeX
@article{chalas2026respi
author = {Chalas, Nikos and Saltafossi, Martina and Berther, Teresa and Balestrieri, Elio and Abbasi, Omid and Kluger, Daniel S.},
title = {{Respiration as a dynamic modulator of sensory sampling}},
journal = {Nature communications},
year = {2026},
month = apr,
volume = {17},
number = {1},
pages = {3261},
publisher = {Nature Publishing Group},
issn = {2041-1723},
doi = {10.1038/
url = {https://
pmid = {41946728},
pmcid = {PMC13062087}
}
RIS
TY - JOUR
AU - Chalas, Nikos
AU - Saltafossi, Martina
AU - Berther, Teresa
AU - Balestrieri, Elio
AU - Abbasi, Omid
AU - Kluger, Daniel S.
TI - Respiration as a dynamic modulator of sensory sampling
T2 - Nature communications
J2 - Nat Commun
PY - 2026
DA - 2026/
VL - 17
IS - 1
SP - 3261
SN - 2041-1723
PB - Nature Publishing Group
DO - 10.1038/
UR - https://
LA - en
ER -
CSL-JSON
{
"id": "10.1038/
"type": "article-journal",
"title": "Respiration as a dynamic modulator of sensory sampling",
"container-title": "Nature communications",
"author": [
{
"family": "Chalas",
"given": "Nikos"
},
{
"family": "Saltafossi",
"given": "Martina"
},
{
"family": "Berther",
"given": "Teresa"
},
{
"family": "Balestrieri",
"given": "Elio"
},
{
"family": "Abbasi",
"given": "Omid"
},
{
"family": "Kluger",
"given": "Daniel S."
}
],
"container-title-short":
"volume": "17",
"issue": "1",
"page": "3261",
"DOI": "10.1038/
"PMID": "41946728",
"PMCID": "PMC13062087",
"ISSN": "2041-1723",
"publisher": "Nature Publishing Group",
"URL": "https://
"language": "en",
"issued": {
"date-parts": [
[
2026,
4,
7
]
]
}
}
The tracing map gets a citation of its own once an author has validated it and it has a DOI.
Similar papers
The papers with a page that share the most with this one: the tools found in their code, their categories, datasets, cited references and authors, the rarest counting most.
- [1] doi:10.1523/jneurosci.1173-24.2024 [code]
- A Dynamic Link between Respiration and Arousal.Journal: The Journal of neuroscience : the official journal of the Society for NeuroscienceIn common: FieldTrip, Signal Processing Toolbox, Statistics and Machine Learning Toolbox, 19 references, author Daniel S Kluger
- [2] doi:10.1038/s41467-023-40250-9 [code]
- Modulatory dynamics of periodic and aperiodic activity in respiration-brain coupling.Journal: Nature communicationsIn common: Violinplot-Matlab, CircStat, FieldTrip, 1 other tool, cognitive, 13 references, author Daniel S Kluger
- [3] doi:10.1371/journal.pcbi.1014672 [code]
- Robust circular cluster-based statistics for respiration-brain coupling.Journal: PLoS computational biologyIn common: Signal Processing Toolbox, Statistics and Machine Learning Toolbox, 12 references, author Daniel S Kluger
- [4] doi:10.1038/s43856-025-00811-z [code]
- Respiratory modulations of cortical excitability and interictal spike timing in focal epilepsy: a case report.Journal: Communications medicineIn common: CircStat, FieldTrip, Signal Processing Toolbox, 1 other tool, 8 references, author Daniel S Kluger
- [5] doi:10.1371/journal.pcbi.1013086 [code]
- The respiratory cycle modulates distinct dynamics of affective and perceptual decision-making.Journal: PLoS computational biologyIn common: cognitive, 14 references
- [6] doi:10.1113/jp287205 [code]
- The pupillary respiratory-phase response: pupil size is smallest around inhalation onset and largest during exhalation.Journal: The Journal of physiologyIn common: CircStat, Statistics and Machine Learning Toolbox, cognitive, 9 references
- [7] doi:10.1371/journal.pbio.3003982 [code]
- No evidence for modulation of the readiness potential by respiratory phase during natural breathing.Journal: PLoS biologyIn common: CircStat, FieldTrip, Signal Processing Toolbox, 1 other tool, 7 references
- [8] doi:10.1038/s41467-026-73828-0 [code]
- Human forebrain neural synchronization and entrainment to breathing during wakefulness, sleep, and external mechanical ventilation.Journal: Nature communicationsIn common: Violinplot-Matlab, Signal Processing Toolbox, Statistics and Machine Learning Toolbox, 8 references
- [9] doi:10.1098/rspb.2024.2566 [code]
- Respiration shapes response speed and accuracy with a systematic time lag.Journal: Proceedings. Biological sciencesIn common: cognitive, 11 references
- [10] doi:10.1038/s42003-024-06642-3
- Dynamic mechanisms that couple the brain and breathing to the external environment.Journal: Communications biologyIn common: 10 references
Contribute
The authors of this paper can claim it, correct its record and validate its tracing map, and the maintainers of its code (its owner, or a public member of its organization) correct what it says of their repository; anyone signed in can ask for its removal. Every request goes to OSCR's own machine, which answers it; your account page follows them.
Sign in with ORCID to claim this paper as one of its authors, correct its record or validate its tracing map: when the paper's metadata lists your ORCID iD, you are recognized at once. Maintainers of its code: sign in with GitHub, then claim the repository on your account page.
Claim this paper
Correct its record
Say what each link of this record is, remove the ones that are not the paper's, add the ones that are missing. The correction becomes a new version of the record, in its Versions section.
Validate its tracing map
You validate the map as this page shows it: 1 repository of the authors' code, each at its verified commit and with its license, 8 scripts, and 5 matches between paragraphs and code (see the Code and Map sections). It then receives a DOI on Zenodo, with you (your ORCID iD) and OSCR as its creators; the code itself is not deposited.
The map's fingerprint: sha256:72e0f85cb229a64d…
Add the badge to the code's README
If you are an author or a maintainer of this paper's code, you can add this badge. It says what OSCR knows of the paper (its code, its validated tracing map) and links the code to this page. Paste one of these snippets under the README's title. Keep the two marks around it, so that the badge can be replaced in place and is never added twice. OSCR changes nothing in your repository.
Markdown
<!-- oscr-badge:start -->
[](https://openscicode.org/badge/doi/10.1038/s41467-026-71604-8)
<!-- oscr-badge:end -->reStructuredText
.. oscr-badge:start
.. image:: https://openscicode.org/badge/doi/10.1038/s41467-026-71604-8.svg
:target: https://openscicode.org/badge/doi/10.1038/s41467-026-71604-8
:alt: OSCR
.. oscr-badge:endHTML
<!-- oscr-badge:start -->
<a href="https://openscicode.org/badge/doi/10.1038/s41467-026-71604-8"><img src="https://openscicode.org/badge/doi/10.1038/s41467-026-71604-8.svg" alt="OSCR"></a>
<!-- oscr-badge:end -->If you choose to, you can have the badge kept up to date for you, with oscr badge --write from the command line or with the GitHub Action in your repository's workflows. Every other form and option (a card, AsciiDoc, LaTeX, shields.io) is on the badge page.
Request its removal
To ask OSCR to remove this record, the copies of its authors' scripts or its tracing map, use the removal request page: signed in, you say who you are, what to remove and why, then review and confirm the request. Published rules decide every request (how).
Discussion, reproductions, activity
Discussion: code errors and mismatches between the code and the paper about this paper and its code, from readers and its authors, beside the paper and its tracing map.
Report a mismatch between the code and the paper · Report a code error (signed in)
Reproductions: reports from readers who ran the authors' code and did not get the paper's result, what they ran, with which environment, commit and data.
Report a reproduction failure (signed in)
Activity: what happens around this paper: new versions of its record, its map's validation, discussions and reproductions. It opens with sign-in.
