The effects of respiratory CO2 fluctuations in the resting-state BOLD signal differ between eyes open and eyes closed

Agreement in cerebrovascular reactivity assessed with diffuse correlation spectroscopy across experimental paradigms improves with short separation regression

Significance

Cerebrovascular reactivity (CVR), i.e., the ability of cerebral vasculature to dilate or constrict in response to vasoactive stimuli, is a biomarker of vascular health. Exogenous administration of inhaled carbon dioxide, i.e., hypercapnia (HC), remains the "gold-standard" intervention to assess CVR. More tolerable paradigms that enable CVR quantification when HC is difficult/contraindicated have been proposed. However, because these paradigms feature mechanistic differences in action, an assessment of agreement of these more tolerable paradigms to HC is needed.

Aim

We aim to determine the agreement of CVR assessed during HC, breath-hold (BH), and resting state (RS) paradigms.

Approach

Healthy adults were subject to HC, BH, and RS paradigms. End tidal carbon dioxide (EtCO2) and cerebral blood flow (CBF, assessed with diffuse correlation spectroscopy) were monitored continuously. CVR (%/mmHg) was quantified via linear regression of CBF versus EtCO2 or via a general linear model (GLM) that was used to minimize the influence of systemic and extracerebral signal contributions.

Results

Strong agreement ( CCC ≥ 0.69 ; R ≥ 0.76 ) among CVR paradigms was demonstrated when utilizing a GLM to regress out systemic/extracerebral signal contributions. Linear regression alone showed poor agreement across paradigms ( CCC ≤ 0.35 ; R ≤ 0.45 ).

Conclusions

More tolerable experimental paradigms coupled with regression of systemic/extracerebral signal contributions may offer a viable alternative to HC for assessing CVR.

Study Protocol: The Heart and Brain Study

BACKGROUND

It is well-established that what is good for the heart is good for the brain. Vascular factors such as hypertension, diabetes, and high cholesterol, and genetic factors such as the apolipoprotein E4 allele increase the risk of developing both cardiovascular disease and dementia. However, the mechanisms underlying the heart-brain association remain unclear. Recent evidence suggests that impairments in vascular phenotypes and cerebrovascular reactivity (CVR) may play an important role in cognitive decline. The Heart and Brain Study combines state-of-the-art vascular ultrasound, cerebrovascular magnetic resonance imaging (MRI) and cognitive testing in participants of the long-running Whitehall II Imaging cohort to examine these processes together. This paper describes the study protocol, data pre-processing and overarching objectives.

METHODS AND DESIGN

The 775 participants of the Whitehall II Imaging cohort, aged 65 years or older in 2019, have received clinical and vascular risk assessments at 5-year-intervals since 1985, as well as a 3T brain MRI scan and neuropsychological tests between 2012 and 2016 (Whitehall II Wave MRI-1). Approximately 25% of this cohort are selected for the Heart and Brain Study, which involves a single testing session at the University of Oxford (Wave MRI-2). Between 2019 and 2023, participants will undergo ultrasound scans of the ascending aorta and common carotid arteries, measures of central and peripheral blood pressure, and 3T MRI scans to measure CVR in response to 5% carbon dioxide in air, vessel-selective cerebral blood flow (CBF), and cerebrovascular lesions. The structural and diffusion MRI scans and neuropsychological battery conducted at Wave MRI-1 will also be repeated. Using this extensive life-course data, the Heart and Brain Study will examine how 30-year trajectories of vascular risk throughout midlife (40-70 years) affect vascular phenotypes, cerebrovascular health, longitudinal brain atrophy and cognitive decline at older ages.

DISCUSSION

The study will generate one of the most comprehensive datasets to examine the longitudinal determinants of the heart-brain association. It will evaluate novel physiological processes in order to describe the optimal window for managing vascular risk in order to delay cognitive decline. Ultimately, the Heart and Brain Study will inform strategies to identify at-risk individuals for targeted interventions to prevent or delay dementia.

Cerebrovascular Reactivity Mapping Without Gas Challenges: A Methodological Guide

Cerebrovascular reactivity (CVR) is defined as the ability of vessels to alter their caliber in response to vasoactive factors, by means of dilating or constricting, in order to increase or decrease regional cerebral blood flow (CBF). Importantly, CVR may provide a sensitive biomarker for pathologies where vasculature is compromised. Furthermore, the spatiotemporal dynamics of CVR observed in healthy subjects, reflecting regional differences in cerebral vascular tone and response, may also be important in functional MRI studies based on neurovascular coupling mechanisms. Assessment of CVR is usually based on the use of a vasoactive stimulus combined with a CBF measurement technique. Although transcranial Doppler ultrasound has been frequently used to obtain global flow velocity measurements, MRI techniques are being increasingly employed for obtaining CBF maps. For the vasoactive stimulus, vasodilatory hypercapnia is usually induced through the manipulation of respiratory gases, including the inhalation of increased concentrations of carbon dioxide. However, most of these methods require an additional apparatus and complex setups, which not only may not be well-tolerated by some populations but are also not widely available. For these reasons, strategies based on voluntary breathing fluctuations without the need for external gas challenges have been proposed. These include the task-based methodologies of breath holding and paced deep breathing, as well as a new generation of methods based on spontaneous breathing fluctuations during resting-state. Despite the multitude of alternatives to gas challenges, existing literature lacks definitive conclusions regarding the best practices for the vasoactive modulation and associated analysis protocols. In this work, we perform an extensive review of CVR mapping techniques based on MRI and CO2 variations without gas challenges, focusing on the methodological aspects of the breathing protocols and corresponding data analysis. Finally, we outline a set of practical guidelines based on generally accepted practices and available data, extending previous reports and encouraging the wider application of CVR mapping methodologies in both clinical and academic MRI settings.

Assessing low-frequency oscillations in cerebrovascular diseases and related conditions with near-infrared spectroscopy: a plausible method for evaluating cerebral autoregulation?

BACKGROUND

Cerebral autoregulation (CA) is the brain's ability to always maintain an adequate and relatively constant blood supply, which is often impaired in cerebrovascular diseases. Near-infrared spectroscopy (NIRS) examines oxygenated hemoglobin (OxyHb) in the cerebral cortex. Low- and very low-frequency oscillations ( LFOs ≈ 0.1    Hz and VLFOs ≈ 0.05 to 0.01 Hz) in OxyHb have been proposed to reflect CA.

AIM

To systematically review published results on OxyHb LFOs and VLFOs in cerebrovascular diseases and related conditions measured with NIRS.

APPROACH

A systematic search was performed in the MEDLINE database, which generated 36 studies relevant for inclusion.

RESULTS

Healthy people have relatively stable LFOs. LFO amplitude seems to reflect myogenic CA being decreased by vasomotor paralysis in stroke, by smooth muscle damage or as compensatory action in other conditions but can also be influenced by the sympathetic tone. VLFO amplitude is believed to reflect neurogenic and metabolic CA and is lower in stroke, atherosclerosis, and with aging. Both LFO and VLFO synchronizations appear disturbed in stroke, while the former is also altered in internal carotid stenosis and hypertension.

CONCLUSION

We conclude that amplitudes of LFOs and VLFOs are relatively robust measures for evaluating mechanisms of CA and synchronization analyses can show temporal disruption of CA. Further research and more coherent methodologies are needed.

Reprint of 'Noise contributions to the fMRI signal: An Overview'

The ability to discriminate signal from noise plays a key role in the analysis and interpretation of functional magnetic resonance imaging (fMRI) measures of brain activity. Over the past two decades, a number of major sources of noise have been identified, including system-related instabilities, subject motion, and physiological fluctuations. This article reviews the characteristics of the various noise sources as well as the mechanisms through which they affect the fMRI signal. Approaches for distinguishing signal from noise and the associated challenges are also reviewed. These challenges reflect the fact that some noise sources, such as respiratory activity, are generated by the same underlying brain networks that give rise to functional signals that are of interest.

The Effect of Low-Frequency Physiological Correction on the Reproducibility and Specificity of Resting-State fMRI Metrics: Functional Connectivity, ALFF, and ReHo

The resting-state fMRI (rs-fMRI) signal is affected by a variety of low-frequency physiological phenomena, including variations in cardiac-rate (CRV), respiratory-volume (RVT), and end-tidal CO₂ (PETCO₂). While these effects have become better understood in recent years, the impact that their correction has on the quality of rs-fMRI measurements has yet to be clarified. The objective of this paper is to investigate the effect of correcting for CRV, RVT and PETCO₂ on the rs-fMRI measurements. Nine healthy subjects underwent a test-retest rs-fMRI acquisition using repetition times (TRs) of 2 s (long-TR) and 0.323 s (short-TR), and the data were processed using eight different physiological correction strategies. Subsequently, regional homogeneity (ReHo), amplitude of low-frequency fluctuation (ALFF), and resting-state connectivity of the motor and default-mode networks are calculated for each strategy. Reproducibility is calculated using intra-class correlation and the Dice Coefficient, while the accuracy of functional-connectivity measures is assessed through network separability, sensitivity and specificity. We found that: (1) the reproducibility of the rs-fMRI measures improved significantly after correction for PETCO₂; (2) separability of functional networks increased after PETCO₂ correction but was not affected by RVT and CRV correction; (3) the effect of physiological correction does not depend on the data sampling-rate; (4) the effect of physiological processes and correction strategies is network-specific. Our findings highlight limitations in our understanding of rs-fMRI quality measures, and underscore the importance of using multiple quality measures to determine the optimal physiological correction strategy.

Biomarkers, designs, and interpretations of resting-state fMRI in translational pharmacological research: A review of state-of-the-Art, challenges, and opportunities for studying brain chemistry

A decade of research and development in resting-state functional MRI (RSfMRI) has opened new translational and clinical research frontiers. This review aims to bridge between technical and clinical researchers who seek reliable neuroimaging biomarkers for studying drug interactions with the brain. About 85 pharma-RSfMRI studies using BOLD signal (75% of all) or arterial spin labeling (ASL) were surveyed to investigate the acute effects of psychoactive drugs. Experimental designs and objectives include drug fingerprinting dose-response evaluation, biomarker validation and calibration, and translational studies. Common biomarkers in these studies include functional connectivity, graph metrics, cerebral blood flow and the amplitude and spectrum of BOLD fluctuations. Overall, RSfMRI-derived biomarkers seem to be sensitive to spatiotemporal dynamics of drug interactions with the brain. However, drugs cause both central and peripheral effects, thus exacerbate difficulties related to biological confounds, structured noise from motion and physiological confounds, as well as modeling and inference testing. Currently, these issues are not well explored, and heterogeneities in experimental design, data acquisition and preprocessing make comparative or meta-analysis of existing reports impossible. A unifying collaborative framework for data-sharing and data-mining is thus necessary for investigating the commonalities and differences in biomarker sensitivity and specificity, and establishing guidelines. Multimodal datasets including sham-placebo or active control sessions and repeated measurements of various psychometric, physiological, metabolic and neuroimaging phenotypes are essential for pharmacokinetic/pharmacodynamic modeling and interpretation of the findings. We provide a list of basic minimum and advanced options that can be considered in design and analyses of future pharma-RSfMRI studies. Hum Brain Mapp 38:2276-2325, 2017. © 2017 Wiley Periodicals, Inc.

Noise contributions to the fMRI signal: An overview

The ability to discriminate signal from noise plays a key role in the analysis and interpretation of functional magnetic resonance imaging (fMRI) measures of brain activity. Over the past two decades, a number of major sources of noise have been identified, including system-related instabilities, subject motion, and physiological fluctuations. This article reviews the characteristics of the various noise sources as well as the mechanisms through which they affect the fMRI signal. Approaches for distinguishing signal from noise and the associated challenges are also reviewed. These challenges reflect the fact that some noise sources, such as respiratory activity, are generated by the same underlying brain networks that give rise to functional signals that are of interest.

Quantitative mapping of cerebrovascular reactivity using resting-state BOLD fMRI: Validation in healthy adults

In conventional neuroimaging, cerebrovascular reactivity (CVR) is quantified primarily using the blood-oxygenation level-dependent (BOLD) functional MRI (fMRI) signal, specifically, as the BOLD response to intravascular carbon dioxide (CO2) modulations, in units of [%ΔBOLD/mmHg]. While this method has achieved wide appeal and clinical translation, the tolerability of CO2-related tasks amongst patients and the elderly remains a challenge in more routine and large-scale applications. In this work, we propose an improved method to quantify CVR by exploiting intrinsic fluctuations in CO2 and corresponding changes in the resting-state BOLD signal (rs-qCVR). Our rs-qCVR approach requires simultaneous monitoring of PETCO2, cardiac pulsation and respiratory volume. In 16 healthy adults, we compare our quantitative CVR estimation technique to the prospective CO2-targeting based CVR quantification approach (qCVR, the "standard"). We also compare our rs-CVR to non-quantitative alternatives including the resting-state fluctuation amplitude (RSFA), amplitude of low-frequency fluctuation (ALFF) and global-signal regression. When all subjects were pooled, only RSFA and ALFF were significantly associated with qCVR. However, for characterizing regional CVR variations within each subject, only the PETCO2-based rs-qCVR measure is strongly associated with standard qCVR in 100% of the subjects (p≤0.1). In contrast, for the more qualitative CVR measures, significant within-subject association with qCVR was only achieved in 50-70% of the subjects. Our work establishes the feasibility of extracting quantitative CVR maps using rs-fMRI, opening the possibility of mapping functional connectivity and qCVR simultaneously.

Functional connectivity and information flow of the respiratory neural network in chronic obstructive pulmonary disease

Breathing involves a complex interplay between the brainstem automatic network and cortical voluntary command. How these brain regions communicate at rest or during inspiratory loading is unknown. This issue is crucial for several reasons: (i) increased respiratory loading is a major feature of several respiratory diseases, (ii) failure of the voluntary motor and cortical sensory processing drives is among the mechanisms that precede acute respiratory failure, (iii) several cerebral structures involved in responding to inspiratory loading participate in the perception of dyspnea, a distressing symptom in many disease. We studied functional connectivity and Granger causality of the respiratory network in controls and patients with chronic obstructive pulmonary disease (COPD), at rest and during inspiratory loading. Compared with those of controls, the motor cortex area of patients exhibited decreased connectivity with their contralateral counterparts and no connectivity with the brainstem. In the patients, the information flow was reversed at rest with the source of the network shifted from the medulla towards the motor cortex. During inspiratory loading, the system was overwhelmed and the motor cortex became the sink of the network. This major finding may help to understand why some patients with COPD are prone to acute respiratory failure. Network connectivity and causality were related to lung function and illness severity. We validated our connectivity and causality results with a mathematical model of neural network. Our findings suggest a new therapeutic strategy involving the modulation of brain activity to increase motor cortex functional connectivity and improve respiratory muscles performance in patients. Hum Brain Mapp 37:2736–2754, 2016. © 2016 The Authors Human Brain Mapping Published by Wiley Periodicals, Inc.

The association between cerebrovascular reactivity and resting-state fMRI functional connectivity in healthy adults: The influence of basal carbon dioxide

Although widely used in resting-state fMRI (fMRI) functional connectivity measurement (fcMRI), the BOLD signal is only an indirect measure of neuronal activity, and is inherently modulated by both neuronal activity and vascular physiology. For instance, cerebrovascular reactivity (CVR) varies widely across individuals irrespective of neuronal function, but the implications for fcMRI are currently unknown. This knowledge gap compromises our ability to correctly interpret fcMRI measurements. In this work, we investigate the relationship between CVR and resting fcMRI measurements in healthy young adults, in both the motor and the executive-control networks. We modulate CVR within each individual by subtly increasing and decreasing resting vascular tension through baseline end-tidal CO2 (PETCO2), and measure fcMRI during these hypercapnic, hypocapnic and normocapnic states. Furthermore, we assess the association between CVR and fcMRI within and across individuals. Within individuals, resting PETCO2 is found to significantly influence both CVR and resting fcMRI values. In addition, we find resting fcMRI to be significantly and positively associated with CVR across the group in both networks. This relationship is potentially mediated by concomitant alterations in BOLD signal fluctuation amplitude. This work clearly demonstrates and quantifies a major vascular modulator of resting fcMRI, one that is also subject and regional dependent. We suggest that individualized correction for CVR effects in fcMRI measurements is essential for fcMRI studies of healthy brains, and can be even more important in studying diseased brains.

Frequency-Dependent Modulation of Regional Synchrony in the Human Brain by Eyes Open and Eyes Closed Resting-States

The eyes-open (EO) and eyes-closed (EC) states have differential effects on BOLD-fMRI signal dynamics, affecting both the BOLD oscillation frequency of a single voxel and the regional homogeneity (ReHo) of several neighboring voxels. To explore how the two resting-states modulate the local synchrony through different frequency bands, we decomposed the time series of each voxel into several components that fell into distinct frequency bands. The ReHo in each of the bands was calculated and compared between the EO and EC conditions. The cross-voxel correlations between the mean frequency and the overall ReHo of each voxel’s original BOLD series in different brain areas were also calculated and compared between the two states. Compared with the EC state, ReHo decreased with EO in a wide frequency band of 0.01–0.25 Hz in the bilateral thalamus, sensorimotor network, and superior temporal gyrus, while ReHo increased significantly in the band of 0–0.01 Hz in the primary visual cortex, and in a higher frequency band of 0.02–0.1 Hz in the higher order visual areas. The cross-voxel correlations between the frequency and overall ReHo were negative in all the brain areas but varied from region to region. These correlations were stronger with EO in the visual network and the default mode network. Our results suggested that different frequency bands of ReHo showed different sensitivity to the modulation of EO-EC states. The better spatial consistency between the frequency and overall ReHo maps indicated that the brain might adopt a stricter frequency-dependent configuration with EO than with EC.

Characterizing Resting-State Brain Function Using Arterial Spin Labeling

Arterial spin labeling (ASL) is an increasingly established magnetic resonance imaging (MRI) technique that is finding broader applications in studying the healthy and diseased brain. This review addresses the use of ASL to assess brain function in the resting state. Following a brief technical description, we discuss the use of ASL in the following main categories: (1) resting-state functional connectivity (FC) measurement: the use of ASL-based cerebral blood flow (CBF) measurements as an alternative to the blood oxygen level-dependent (BOLD) technique to assess resting-state FC; (2) the link between network CBF and FC measurements: the use of network CBF as a surrogate of the metabolic activity within corresponding networks; and (3) the study of resting-state dynamic CBF-BOLD coupling and cerebral metabolism: the use of dynamic CBF information obtained using ASL to assess dynamic CBF-BOLD coupling and oxidative metabolism in the resting state. In addition, we summarize some future challenges and interesting research directions for ASL, including slice-accelerated (multiband) imaging as well as the effects of motion and other physiological confounds on perfusion-based FC measurement. In summary, this work reviews the state-of-the-art of ASL and establishes it as an increasingly viable MRI technique with high translational value in studying resting-state brain function.

Agreement and repeatability of vascular reactivity estimates based on a breath-hold task and a resting state scan

FMRI BOLD responses to changes in neural activity are influenced by the reactivity of the vasculature. By complementing a task-related BOLD acquisition with a vascular reactivity measure obtained through breath-holding or hypercapnia, this unwanted variance can be statistically reduced in the BOLD responses of interest. Recently, it has been suggested that vascular reactivity can also be estimated using a resting state scan.

This study aimed to compare three breath-hold based analysis approaches (block design, sine–cosine regressor and CO₂ regressor) and a resting state approach (CO₂ regressor) to measure vascular reactivity. We tested BOLD variance explained by the model and repeatability of the measures. Fifteen healthy participants underwent a breath-hold task and a resting state scan with end-tidal CO₂ being recorded during both. Vascular reactivity was defined as CO₂-related BOLD percent signal change/mm Hg change in CO₂.

Maps and regional vascular reactivity estimates showed high repeatability when the breath-hold task was used. Repeatability and variance explained by the CO₂ trace regressor were lower for the resting state data based approach, which resulted in highly variable measures of vascular reactivity.

We conclude that breath-hold based vascular reactivity estimations are more repeatable than resting-based estimates, and that there are limitations with replacing breath-hold scans by resting state scans for vascular reactivity assessment.

Mapping the end-tidal CO2 response function in the resting-state BOLD fMRI signal: spatial specificity, test-retest reliability and effect of fMRI sampling rate

The blood oxygenation level dependent (BOLD) signal measures brain function indirectly through physiological processes and hence is susceptible to global physiological changes. Specifically, fluctuations in end-tidal CO2 (PETCO2), in addition to cardiac rate variation (CRV), and respiratory volume per time (RVT) variations, have been known to confound the resting-state fMRI (rs-fMRI) signal. Previous studies addressed the resting-state fMRI response function to CRV and RVT, but no attempt has been made to directly estimate the voxel-wise response function to PETCO2. Moreover, the potential interactions among PETCO2, CRV, and RVT necessitate their simultaneous inclusion in a multi-regression model to estimate the PETCO2 response. In this study, we use such a model to estimate the voxel-wise PETCO2 response functions directly from rs-fMRI data of nine healthy subjects. We also characterized the effect of sampling rate (TR=2seconds vs. 323ms) on the temporal and spatial variability of the PETCO2 response function in addition to that of CRV and RVT. In addition, we assess the test-retest reproducibility of the response functions to PETCO2, CRV and RVT. We found that despite overlaps across their spatial patterns, PETCO2 explains a unique portion of the rs-fMRI signal variance compared to RVT and CRV. We also found the shapes of the estimated responses are very similar between long- and short-TR data, although responses estimated from short-TR data have higher reproducibility.

Response of end tidal CO2 pressure to impulse exercise

The purpose of the present study was to examine how end tidal CO(2) pressure (PETCO(2)) is controlled in impulse exercise. After pre-exercise at 25 watts for 5 min, impulse exercise for 10 sec with 200 watts followed by post exercise at 25 watts was performed. Ventilation (VE) significantly increased until the end of impulse exercise and significantly re-increased after a sudden decrease. Heart rate (HR) significantly increased until the end of impulse exercise and then decreased to the pre-exercise level. PETCO(2) remained constant during impulse exercise. PETCO(2) significantly increased momentarily after impulse exercise and then significantly decreased to the pre-exercise level. PETCO(2) showed oscillation. The average peak frequency of power spectral density in PETCO(2) appeared at 0.0078 Hz. Cross correlations were obtained after impulse exercise. The peak cross correlations between VE and PETCO(2), HR and PETCO(2), and VE and HR were 0.834 with a time delay of -7 sec, 0.813 with a time delay of 7 sec and 0.701 with a time delay of -15 sec, respectively. We demonstrated that PETCO(2) homeodynamics was interactively maintained by PETCO(2) itself, CO(2) transportation (product of cardiac output and mixed venous CO(2) content) into the lungs by heart pumping and CO(2) elimination by ventilation, and it oscillates as a result of their interactions.

Dynamic functional connectivity: promise, issues, and interpretations

The brain must dynamically integrate, coordinate, and respond to internal and external stimuli across multiple time scales. Non-invasive measurements of brain activity with fMRI have greatly advanced our understanding of the large-scale functional organization supporting these fundamental features of brain function. Conclusions from previous resting-state fMRI investigations were based upon static descriptions of functional connectivity (FC), and only recently studies have begun to capitalize on the wealth of information contained within the temporal features of spontaneous BOLD FC. Emerging evidence suggests that dynamic FC metrics may index changes in macroscopic neural activity patterns underlying critical aspects of cognition and behavior, though limitations with regard to analysis and interpretation remain. Here, we review recent findings, methodological considerations, neural and behavioral correlates, and future directions in the emerging field of dynamic FC investigations.

Neural mechanisms underlying breathing complexity

Breathing is maintained and controlled by a network of automatic neurons in the brainstem that generate respiratory rhythm and receive regulatory inputs. Breathing complexity therefore arises from respiratory central pattern generators modulated by peripheral and supra-spinal inputs. Very little is known on the brainstem neural substrates underlying breathing complexity in humans. We used both experimental and theoretical approaches to decipher these mechanisms in healthy humans and patients with chronic obstructive pulmonary disease (COPD). COPD is the most frequent chronic lung disease in the general population mainly due to tobacco smoke. In patients, airflow obstruction associated with hyperinflation and respiratory muscles weakness are key factors contributing to load-capacity imbalance and hence increased respiratory drive. Unexpectedly, we found that the patients breathed with a higher level of complexity during inspiration and expiration than controls. Using functional magnetic resonance imaging (fMRI), we scanned the brain of the participants to analyze the activity of two small regions involved in respiratory rhythmogenesis, the rostral ventro-lateral (VL) medulla (pre-Bötzinger complex) and the caudal VL pons (parafacial group). fMRI revealed in controls higher activity of the VL medulla suggesting active inspiration, while in patients higher activity of the VL pons suggesting active expiration. COPD patients reactivate the parafacial to sustain ventilation. These findings may be involved in the onset of respiratory failure when the neural network becomes overwhelmed by respiratory overload We show that central neural activity correlates with airflow complexity in healthy subjects and COPD patients, at rest and during inspiratory loading. We finally used a theoretical approach of respiratory rhythmogenesis that reproduces the kernel activity of neurons involved in the automatic breathing. The model reveals how a chaotic activity in neurons can contribute to chaos in airflow and reproduces key experimental fMRI findings.