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

Machine learning-based estimation of respiratory fluctuations in a healthy adult population using resting state BOLD fMRI and head motion parameters

PURPOSE

External physiological monitoring is the primary approach to measure and remove effects of low-frequency respiratory variation from BOLD-fMRI signals. However, the acquisition of clean external respiratory data during fMRI is not always possible, so recent research has proposed using machine learning to directly estimate respiratory variation (RV), potentially obviating the need for external monitoring. In this study, we propose an extended method for reconstructing RV waveforms directly from resting state BOLD-fMRI data in healthy adult participants with the inclusion of both BOLD signals and derived head motion parameters.

METHODS

In the proposed method, 1D convolutional neural networks (1D-CNNs) used BOLD signals and head motion parameters to reconstruct the RV waveform for the whole fMRI scan time. Resting-state fMRI data and associated respiratory records from the Human Connectome Project in Young Adults (HCP-YA) dataset are used to train and test the proposed method.

RESULTS

Compared to using only BOLD-fMRI data for a CNN input, this approach yielded improvements of 14% in mean absolute error, 24% in mean square error, 14% in correlation, and 12% in dynamic time warping. When tested on independent datasets, the method demonstrated generalizability, even in data with different TRs and physiological conditions.

CONCLUSION

This study shows that the respiratory variations could be reconstructed from BOLD-fMRI data in the young adult population, and its accuracy could be improved using supportive data such as head motion parameters. The method also performed well on independent datasets with different experimental conditions.

Identifying multilayer network hub by graph representation learning

The recent advances in neuroimaging technology allow us to understand how the human brain is wired in vivo and how functional activity is synchronized across multiple regions. Growing evidence shows that the complexity of the functional connectivity is far beyond the widely used mono-layer network. Indeed, the hierarchical processing information among distinct brain regions and across multiple channels requires using a more advanced multilayer model to understand the synchronization across the brain that underlies functional brain networks. However, the principled approach for characterizing network organization in the context of multilayer topologies is largely unexplored. In this work, we present a novel multi-variate hub identification method that takes both the intra- and inter-layer network topologies into account. Specifically, we put the spotlight on the multilayer graph embeddings that allow us to separate connector hubs (connecting across network modules) with their peripheral nodes. The removal of these hub nodes breaks down the entire multilayer brain network into a set of disconnected communities. We have evaluated our novel multilayer hub identification method in task-based and resting-state functional images. Complimenting ongoing findings using mono-layer brain networks, our multilayer network analysis provides a new understanding of brain network topology that links functional connectivities with brain states and disease progression.

Quantitative mapping of cerebrovascular reactivity amplitude and delay with breath-hold BOLD fMRI when end-tidal CO2 quality is low

Cerebrovascular reactivity (CVR), the ability of cerebral blood vessels to dilate or constrict in order to regulate blood flow, is a clinically useful measure of cerebrovascular health. CVR is often measured using a breath-hold task to modulate blood CO2 levels during an fMRI scan. Measuring end-tidal CO2 (PETCO2) with a nasal cannula during the task allows CVR amplitude to be calculated in standard units (vascular response per unit change in CO2, or %BOLD/mmHg) and CVR delay to be calculated in seconds. The use of standard units allows for normative CVR ranges to be established and for CVR comparisons to be made across subjects and scan sessions. Although breath-holding can be successfully performed by diverse patient populations, obtaining accurate PETCO2 measurements requires additional task compliance; specifically, participants must breathe exclusively through their nose and exhale immediately before and after each breath hold. Meeting these requirements is challenging, even in healthy participants, and this has limited the translational potential of breath-hold fMRI for CVR mapping. Previous work has focused on using alternative regressors such as respiration volume per time (RVT), derived from respiratory belt measurements, to map CVR. Because measuring RVT does not require additional task compliance from participants, it is a more feasible measure than PETCO2. However, using RVT does not produce CVR in standard units. In this work, we explored how to achieve CVR maps, in standard units, when breath-hold task PETCO2 data quality is low. First, we evaluated whether RVT could be scaled to units of mmHg using a subset of PETCO2 data of sufficiently high quality. Second, we explored whether a PETCO2 timeseries predicted from RVT using deep learning allows for more accurate CVR measurements. Using a dense-mapping breath-hold fMRI dataset, we showed that both rescaled RVT and rescaled, predicted PETCO2 can be used to produce maps of CVR amplitude and delay in standard units with strong absolute agreement to ground-truth maps. However, the rescaled, predicted PETCO2 regressor resulted in superior accuracy for both CVR amplitude and delay. In an individual with regions of increased CVR delay due to Moyamoya disease, the predicted PETCO2 regressor also provided greater sensitivity to pathology than RVT. Ultimately, this work will increase the clinical applicability of CVR in populations exhibiting decreased task compliance.

Effects of clamping end-tidal CO2 on neurofluidic low-frequency oscillations

In recent years, low-frequency oscillations (LFOs) (0.01-0.1 Hz) have been a subject of interest in resting-state functional magnetic resonance imaging research. They are believed to have many possible driving mechanisms, from both regional and global sources. Internal fluctuations in the partial pressure of CO2 (PCO2) has long been thought of as one of these major driving forces, but its exact contributions compared with other mechanisms have yet to be fully understood. This study examined the effects of end-tidal PCO2 (PetCO2) oscillations on LF cerebral hemodynamics and cerebrospinal fluid (CSF) dynamics under "clamped PetCO2" and "free-breathing" conditions. Under clamped PetCO2, a participant's PetCO2 levels were fixed to their baseline average, whereas PetCO2 was not controlled in free breathing. Under clamped PetCO2, the fractional amplitude of hemodynamic LFOs in the occipital and sensorimotor cortex and temporal lobes were found to be significantly reduced. Additionally, the fractional amplitude of CSF LFOs, measured at the fourth ventricle, was found to be reduced by almost one-half. However, the spatiotemporal distributions of blood and CSF delay times, as measured by cross-correlation in the LF domain, were not significantly altered between conditions. This study demonstrates that, while PCO2 oscillations significantly mediate LFOs, especially those observed in the CSF, other mechanisms are able to maintain LFOs, with high correlation, even in their absence.

Automating the Human Connectome Project's Temporal ICA Pipeline

Functional magnetic resonance imaging (fMRI) data are dominated by noise and artifacts, with only a small fraction of the variance relating to neural activity. Temporal independent component analysis (tICA) is a recently developed method that enables selective denoising of fMRI artifacts related to physiology such as respiration. However, an automated and easy to use pipeline for tICA has not previously been available; instead, two manual steps have been necessary: 1) setting the group spatial ICA dimensionality after MELODIC's Incremental Group-PCA (MIGP) and 2) labeling tICA components as artifacts versus signals. Moreover, guidance has been lacking as to how many subjects and timepoints are needed to adequately re-estimate the temporal ICA decomposition and what alternatives are available for smaller groups or even individual subjects. Here, we introduce a nine-step fully automated tICA pipeline which removes global artifacts from fMRI dense timeseries after sICA+FIX cleaning and MSMAll alignment driven by functionally relevant areal features. Additionally, we have developed an automated "reclean" Pipeline for improved spatial ICA (sICA) artifact removal. Two major automated components of the pipeline are 1) an automatic group spatial ICA (sICA) dimensionality selection for MIGP data enabled by fitting multiple Wishart distributions; 2) a hierarchical classifier to distinguish group tICA signal components from artifactual components, equipped with a combination of handcrafted features from domain expert knowledge and latent features obtained via self-supervised learning on spatial maps. We demonstrate that the dimensionality estimated for the MIGP data from HCP Young Adult 3T and 7T datasets is comparable to previous manual tICA estimates, and that the group sICA decomposition is highly reproducible. We also show that the tICA classifier achieved over 0.98 Precision-Recall Area Under Curve (PR-AUC) and that the correctly classified components account for over 95% of the tICA-represented variance on multiple held-out evaluation datasets including the HCP-Young Adult, HCP-Aging and HCP-Development datasets under various settings. Our automated tICA pipeline is now available as part of the HCP pipelines, providing a powerful and user-friendly tool for the neuroimaging community.

Denoising task-correlated head motion from motor-task fMRI data with multi-echo ICA

Motor-task functional magnetic resonance imaging (fMRI) is crucial in the study of several clinical conditions, including stroke and Parkinson's disease. However, motor-task fMRI is complicated by task-correlated head motion, which can be magnified in clinical populations and confounds motor activation results. One method that may mitigate this issue is multi-echo independent component analysis (ME-ICA), which has been shown to separate the effects of head motion from the desired blood oxygenation level dependent (BOLD) signal but has not been tested in motor-task datasets with high amounts of motion. In this study, we collected an fMRI dataset from a healthy population who performed a hand grasp task with and without task-correlated amplified head motion to simulate a motor-impaired population. We analyzed these data using three models: single-echo (SE), multi-echo optimally combined (ME-OC), and ME-ICA. We compared the models' performance in mitigating the effects of head motion on the subject level and group level. On the subject level, ME-ICA better dissociated the effects of head motion from the BOLD signal and reduced noise. Both ME models led to increased t-statistics in brain motor regions. In scans with high levels of motion, ME-ICA additionally mitigated artifacts and increased stability of beta coefficient estimates, compared to SE. On the group level, all three models produced activation clusters in expected motor areas in scans with both low and high motion, indicating that group-level averaging may also sufficiently resolve motion artifacts that vary by subject. These findings demonstrate that ME-ICA is a useful tool for subject-level analysis of motor-task data with high levels of task-correlated head motion. The improvements afforded by ME-ICA are critical to improve reliability of subject-level activation maps for clinical populations in which group-level analysis may not be feasible or appropriate, for example, in a chronic stroke cohort with varying stroke location and degree of tissue damage.

Improved cerebrovascular reactivity mapping using coherence weighted general linear model in the frequency domain

Cerebrovascular reactivity (CVR) is a prognostic indicator of cerebrovascular health. Estimating CVR from endogenous end-tidal carbon dioxide (CO2) fluctuation and MRI signal recorded under resting state can be difficult due to the poor signal-to-noise ratio (SNR) of signals. Thus, we aimed to improve the method of estimating CVR from end-tidal CO2 and MRI signals. We proposed a coherence weighted general linear model (CW-GLM) to estimate CVR from the Fourier coefficients weighted by the signal coherence in frequency domain, which confers two advantages. First, it requires no signal alignment in time domain, which simplifies experimental methods. Second, it limits the GLM analysis within the frequency band where CO2 and MRI signals are highly correlated, which automatically suppresses noise and nuisance signals. We compared the performance of our method with time-domain GLM (TD-GLM) and frequency-domain GLM (FD-GLM) in both synthetic and in-vivo data; wherein we calculated CVR from signals recorded under both resting state and sinusoidal stimulus. In synthetic data, CW-GLM has a remarkable performance on CVR estimation from narrow band signals with a mean-absolute error of 0.7 % (gray matter) and 1.2 % (white matter), which was lower than all the other methods. Meanwhile, CW-GLM maintains a comparable performance on CVR estimation from resting signals, with a mean-absolute error of 4.1 % (gray matter) and 8 % (white matter). The superior performance was maintained across the 36 in-vivo measurements, with CW-GLM exhibiting limits of agreement of -16.7 % - 9.5 % between CVR calculated from the resting and sinusoidal CO2 paradigms which was 12 % - 209 % better than current time-domain methods. Evaluating of the cross-coherence spectrum revealed highest signal coherence within the frequency band from 0.01 Hz to 0.05 Hz, which overlaps with previously recommended frequency band (0.02 Hz to 0.04 Hz) for CVR analysis. Our data demonstrates that CW-GLM can work as a self-adaptive band-pass filter to improve CVR robustness, while also avoiding the need for signal temporal alignment.

Denoising task-correlated head motion from motor-task fMRI data with multi-echo ICA

Motor-task functional magnetic resonance imaging (fMRI) is crucial in the study of several clinical conditions, including stroke and Parkinson's disease. However, motor-task fMRI is complicated by task-correlated head motion, which can be magnified in clinical populations and confounds motor activation results. One method that may mitigate this issue is multi-echo independent component analysis (ME-ICA), which has been shown to separate the effects of head motion from the desired BOLD signal but has not been tested in motor-task datasets with high amounts of motion. In this study, we collected an fMRI dataset from a healthy population who performed a hand grasp task with and without task-correlated amplified head motion to simulate a motor-impaired population. We analyzed these data using three models: single-echo (SE), multi-echo optimally combined (ME-OC), and ME-ICA. We compared the models' performance in mitigating the effects of head motion on the subject level and group level. On the subject level, ME-ICA better dissociated the effects of head motion from the BOLD signal and reduced noise. Both ME models led to increased t-statistics in brain motor regions. In scans with high levels of motion, ME-ICA additionally mitigated artifacts and increased stability of beta coefficient estimates, compared to SE. On the group level, all three models produced activation clusters in expected motor areas in scans with both low and high motion, indicating that group-level averaging may also sufficiently resolve motion artifacts that vary by subject. These findings demonstrate that ME-ICA is a useful tool for subject-level analysis of motor-task data with high levels of task-correlated head motion. The improvements afforded by ME-ICA are critical to improve reliability of subject-level activation maps for clinical populations in which group-level analysis may not be feasible or appropriate, for example in a chronic stroke cohort with varying stroke location and degree of tissue damage.

Neurofluid coupling during sleep and wake states

BACKGROUND

In clinical populations, the movement of cerebrospinal fluid (CSF) during sleep is a growing area of research with potential mechanistic connections in both neurodegenerative (e.g., Alzheimer's Disease) and neurodevelopmental disorders. However, we know relatively little about the processes that influence CSF movement. To inform clinical intervention targets this study assesses the coupling between (a) real-time CSF movement, (b) neuronal-driven movement, and (c) non-neuronal systemic physiology driven movement.

METHODS

This study included eight young, healthy volunteers, with concurrently acquired neurofluid dynamics using functional Magnetic Resonance Imaging (MRI), neural activity using Electroencephalography (EEG), and non-neuronal systemic physiology with peripheral functional Near-Infrared Spectroscopy (fNIRS). Neuronal and non-neuronal drivers were assessed temporally; wherein, EEG measured slow wave activity that preceded CSF movement was considered neuronally driven. Similarly, slow wave oscillations (assessed via fNIRS) that coupled with CSF movement were considered non-neuronal systemic physiology driven.

RESULTS AND CONCLUSIONS

Our results document neural contributions to CSF movement were only present during light NREM sleep but low-frequency non-neuronal oscillations were strongly coupled with CSF movement in all assessed states - awake, NREM-1, NREM-2. The clinical/research implications of these findings are two-fold. First, neuronal-driven oscillations contribute to CSF movement outside of deep sleep (NREM-3); therefore, interventions aimed at increasing CSF movement may yield meaningful increases with the promotion of NREM sleep more generally - a focus on NREM S3 may not be needed. Second, non-neuronal systemic oscillations contribute across wake and sleep stages; therefore, interventions may increase CSF movement by manipulating systemic physiology.

From Lung to Brain: Respiration Modulates Neural and Mental Activity

Respiration protocols have been developed to manipulate mental states, including their use for therapeutic purposes. In this systematic review, we discuss evidence that respiration may play a fundamental role in coordinating neural activity, behavior, and emotion. The main findings are: (1) respiration affects the neural activity of a wide variety of regions in the brain; (2) respiration modulates different frequency ranges in the brain's dynamics; (3) different respiration protocols (spontaneous, hyperventilation, slow or resonance respiration) yield different neural and mental effects; and (4) the effects of respiration on the brain are related to concurrent modulation of biochemical (oxygen delivery, pH) and physiological (cerebral blood flow, heart rate variability) variables. We conclude that respiration may be an integral rhythm of the brain's neural activity. This provides an intimate connection of respiration with neuro-mental features like emotion. A respiratory-neuro-mental connection holds the promise for a brain-based therapeutic usage of respiration in mental disorders.

Neural signals-based respiratory motion tracking: a proof-of-concept study

Objective.Respiratory motion tracking techniques can provide optimal treatment accuracy for thoracoabdominal radiotherapy and robotic surgery. However, conventional imaging-based respiratory motion tracking techniques are time-lagged owing to the system latency of medical linear accelerators and surgical robots. This study aims to investigate the precursor time of respiratory-related neural signals and analyze the potential of neural signals-based respiratory motion tracking.Approach.The neural signals and respiratory motion from eighteen healthy volunteers were acquired simultaneously using a 256-channel scalp electroencephalography (EEG) system. The neural signals were preprocessed using the MNE python package to extract respiratory-related EEG neural signals. Cross-correlation analysis was performed to assess the precursor time and cross-correlation coefficient between respiratory-related EEG neural signals and respiratory motion.Main results.Respiratory-related neural signals that precede the emergence of respiratory motion are detectable via non-invasive EEG. On average, the precursor time of respiratory-related EEG neural signals was 0.68 s. The representative cross-correlation coefficients between EEG neural signals and respiratory motion of the eighteen healthy subjects varied from 0.22 to 0.87.Significance.Our findings suggest that neural signals have the potential to compensate for the system latency of medical linear accelerators and surgical robots. This indicates that neural signals-based respiratory motion tracking is a potential promising solution to respiratory motion and could be useful in thoracoabdominal radiotherapy and robotic surgery.

Identifying the distinct spectral dynamics of laminar-specific interhemispheric connectivity with bilateral line-scanning fMRI

Despite extensive efforts to identify interhemispheric functional connectivity (FC) with resting-state (rs-) fMRI, correlated low-frequency rs-fMRI signal fluctuation across homotopic cortices originates from multiple sources. It remains challenging to differentiate circuit-specific FC from global regulation. Here, we developed a bilateral line-scanning fMRI method to detect laminar-specific rs-fMRI signals from homologous forepaw somatosensory cortices with high spatial and temporal resolution in rat brains. Based on spectral coherence analysis, two distinct bilateral fluctuation spectral features were identified: ultra-slow fluctuation (<0.04 Hz) across all cortical laminae versus Layer (L) 2/3-specific evoked BOLD at 0.05 Hz based on 4 s on/16 s off block design and resting-state fluctuations at 0.08-0.1 Hz. Based on the measurements of evoked BOLD signal at corpus callosum (CC), this L2/3-specific 0.05 Hz signal is likely associated with neuronal circuit-specific activity driven by the callosal projection, which dampened ultra-slow oscillation less than 0.04 Hz. Also, the rs-fMRI power variability clustering analysis showed that the appearance of L2/3-specific 0.08-0.1 Hz signal fluctuation is independent of the ultra-slow oscillation across different trials. Thus, distinct laminar-specific bilateral FC patterns at different frequency ranges can be identified by the bilateral line-scanning fMRI method.

Interoceptive Awareness of the Breath Preserves Attention and Language Networks amidst Widespread Cortical Deactivation: A Within-Participant Neuroimaging Study

Interoception, the representation of the body's internal state, serves as a foundation for emotion, motivation, and wellbeing. Yet despite its centrality in human experience, the neural mechanisms of interoceptive attention are poorly understood. The Interoceptive/Exteroceptive Attention Task (IEAT) is a novel neuroimaging paradigm that compares behavioral tracking of the respiratory cycle (Active Interoception) to tracking of a visual stimulus (Active Exteroception). Twenty-two healthy participants completed the IEAT during two separate scanning sessions (N = 44) as part of a randomized control trial of mindful awareness in body-oriented therapy (MABT). Compared with Active Exteroception, Active Interoception deactivated somatomotor and prefrontal regions. Greater self-reported interoceptive sensibility (MAIA scale) predicted sparing from deactivation within the anterior cingulate cortex (ACC) and left-lateralized language regions. The right insula, typically described as a primary interoceptive cortex, was only specifically implicated by its deactivation during an exogenously paced respiration condition (Active Matching) relative to self-paced Active Interoception. Psychophysiological interaction (PPI) analysis characterized Active Interoception as promoting greater ACC connectivity with lateral prefrontal and parietal regions commonly referred to as the dorsal attention network (DAN). In contrast to evidence relating accurate detection of liminal interoceptive signals such as the heartbeat to anterior insula activity, interoceptive attention toward salient signals such as the respiratory cycle may involve reduced cortical activity but greater ACC-DAN connectivity, with greater sensibility linked to reduced deactivation within the ACC and language-processing regions.

Direct machine learning reconstruction of respiratory variation waveforms from resting state fMRI data in a pediatric population

In many functional magnetic resonance imaging (fMRI) studies, respiratory signals are unavailable or do not have acceptable quality due to issues with subject compliance, equipment failure or signal error. In large databases, such as the Human Connectome Projects, over half of the respiratory recordings may be unusable. As a result, the direct removal of low frequency respiratory variations from the blood oxygen level-dependent (BOLD) signal time series is not possible. This study proposes a deep learning-based method for reconstruction of respiratory variation (RV) waveforms directly from BOLD fMRI data in pediatric participants (aged 5 to 21 years old), and does not require any respiratory measurement device. To do this, the Lifespan Human Connectome Project in Development (HCP-D) dataset, which includes respiratory measurements, was used to both train a convolutional neural network (CNN) and evaluate its performance. Results show that a CNN can capture informative features from the BOLD signal time course and reconstruct accurate RV timeseries, especially when the subject has a prominent respiratory event. This work advances the use of direct estimation of physiological parameters from fMRI, which will eventually lead to reduced complexity and decrease the burden on participants because they may not be required to wear a respiratory bellows.

Comparing end-tidal CO2, respiration volume per time (RVT), and average gray matter signal for mapping cerebrovascular reactivity amplitude and delay with breath-hold task BOLD fMRI

Cerebrovascular reactivity (CVR), defined as the cerebral blood flow response to a vasoactive stimulus, is an imaging biomarker with demonstrated utility in a range of diseases and in typical development and aging processes. A robust and widely implemented method to map CVR involves using a breath-hold task during a BOLD fMRI scan. Recording end-tidal CO₂ (P_ETCO₂) changes during the breath-hold task is recommended to be used as a reference signal for modeling CVR amplitude in standard units (%BOLD/mmHg) and CVR delay in seconds. However, obtaining reliable P_ETCO₂ recordings requires equipment and task compliance that may not be achievable in all settings. To address this challenge, we investigated two alternative reference signals to map CVR amplitude and delay in a lagged general linear model (lagged-GLM) framework: respiration volume per time (RVT) and average gray matter BOLD response (GM-BOLD). In 8 healthy adults with multiple scan sessions, we compare spatial agreement of CVR maps from RVT and GM-BOLD to those generated with P_ETCO₂. We define a threshold to determine whether a P_ETCO₂ recording has "sufficient" quality for CVR mapping and perform these comparisons in 16 datasets with sufficient P_ETCO₂ and 6 datasets with insufficient P_ETCO₂. When P_ETCO₂ quality is sufficient, both RVT and GM-BOLD produce CVR amplitude maps that are nearly identical to those from P_ETCO₂ (after accounting for differences in scale), with the caveat they are not in standard units to facilitate between-group comparisons. CVR delays are comparable to P_ETCO₂ with an RVT regressor but may be underestimated with the average GM-BOLD regressor. Importantly, when P_ETCO₂ quality is insufficient, RVT and GM-BOLD CVR recover reasonable CVR amplitude and delay maps, provided the participant attempted the breath-hold task. Therefore, our framework offers a solution for achieving high quality CVR maps in both retrospective and prospective studies where sufficient P_ETCO₂ recordings are not available and especially in populations where obtaining reliable measurements is a known challenge (e.g., children). Our results have the potential to improve the accessibility of CVR mapping and to increase the prevalence of this promising metric of vascular health.

Generating dynamic carbon-dioxide traces from respiration-belt recordings: Feasibility using neural networks and application in functional magnetic resonance imaging

INTRODUCTION

In the context of functional magnetic resonance imaging (fMRI), carbon dioxide (CO2) is a well-known vasodilator that has been widely used to monitor and interrogate vascular physiology. Moreover, spontaneous fluctuations in end-tidal carbon dioxide (PETCO2) reflects changes in arterial CO2 and has been demonstrated as the largest physiological noise source for denoising the low-frequency range of the resting-state fMRI (rs-fMRI) signal. However, the majority of rs-fMRI studies do not involve CO2 recordings, and most often only heart rate and respiration are recorded. While the intrinsic link between these latter metrics and CO2 led to suggested possible analytical models, they have not been widely applied.

METHODS

In this proof-of-concept study, we propose a deep-learning (DL) approach to reconstruct CO2 and PETCO2 data from respiration waveforms in the resting state.

RESULTS

We demonstrate that the one-to-one mapping between respiration and CO2 recordings can be well predicted using fully convolutional networks (FCNs), achieving a Pearson correlation coefficient (r) of 0.946 ± 0.056 with the ground truth CO2. Moreover, dynamic PETCO2 can be successfully derived from the predicted CO2, achieving r of 0.512 ± 0.269 with the ground truth. Importantly, the FCN-based methods outperform previously proposed analytical methods. In addition, we provide guidelines for quality assurance of respiration recordings for the purposes of CO2 prediction.

DISCUSSION

Our results demonstrate that dynamic CO2 can be obtained from respiration-volume using neural networks, complementing the still few reports in DL of physiological fMRI signals, and paving the way for further research in DL based bio-signal processing.

Effects of remifentanil on brain responses to noxious stimuli during deep propofol sedation

BACKGROUND

The safety of anaesthesia has improved as a result of better control of anaesthetic depth. However, conventional monitoring does not inform on the nature of nociceptive processes during unconsciousness. A means of inferring the quality of potentially painful experiences could derive from analysis of brain activity using neuroimaging. We have evaluated the dose effects of remifentanil on brain response to noxious stimuli during deep sedation and spontaneous breathing.

METHODS

Optimal data were obtained in 26 healthy subjects. Pressure stimulation that proved to be moderately painful before the experiment was applied to the thumbnail. Functional MRI was acquired in 4-min periods at low (0.5 ng ml^-1), medium (1 ng ml^-1), and high (1.5 ng ml^-1) target plasma concentrations of remifentanil at a stable background infusion of propofol adjusted to induce a state of light unconsciousness.

RESULTS

At low remifentanil doses, we observed partial activation in brain areas processing sensory-discriminative and emotional-affective aspects of pain. At medium doses, relevant changes were identified in structures highly sensitive to general brain arousal, including the brainstem, cerebellum, thalamus, auditory and visual cortices, and the frontal lobe. At high doses, no significant activation was observed.

CONCLUSIONS

The response to moderately intense focal pressure in pain-related brain networks is effectively eliminated with safe remifentanil doses. However, the safety margin in deep sedation-analgesia would be narrowed in minimising not only nociceptive responses, but also arousal-related biological stress.

Using carpet plots to analyze blood transit times in the brain during hypercapnic challenge magnetic resonance imaging

Blood arrival time and blood transit time are useful metrics in characterizing hemodynamic behaviors in the brain. Functional magnetic resonance imaging in combination with a hypercapnic challenge has been proposed as a non-invasive imaging tool to determine blood arrival time and replace dynamic susceptibility contrast (DSC) magnetic resonance imaging, a current gold-standard imaging tool with the downsides of invasiveness and limited repeatability. Using a hypercapnic challenge, blood arrival times can be computed by cross-correlating the administered CO₂ signal with the fMRI signal, which increases during elevated CO₂ due to vasodilation. However, whole-brain transit times derived from this method can be significantly longer than the known cerebral transit time for healthy subjects (nearing 20 s vs. the expected 5-6 s). To address this unrealistic measurement, we here propose a novel carpet plot-based method to compute improved blood transit times derived from hypercapnic blood oxygen level dependent fMRI, demonstrating that the method reduces estimated blood transit times to an average of 5.32 s. We also investigate the use of hypercapnic fMRI with cross-correlation to compute the venous blood arrival times in healthy subjects and compare the computed delay maps with DSC-MRI time to peak maps using the structural similarity index measure (SSIM). The strongest delay differences between the two methods, indicated by low structural similarity index measure, were found in areas of deep white matter and the periventricular region. SSIM measures throughout the remainder of the brain reflected a similar arrival sequence derived from the two methods despite the exaggerated spread of voxel delays computed using CO₂ fMRI.

Assessing the Repeatability of Multi-Frequency Multi-Layer Brain Network Topologies Across Alternative Researcher's Choice Paths

There is a growing interest in the neuroscience community on the advantages of multilayer functional brain networks. Researchers usually treated different frequencies separately at distinct functional brain networks. However, there is strong evidence that these networks share complementary information while their interdependencies could reveal novel findings. For this purpose, neuroscientists adopt multilayer networks, which can be described mathematically as an extension of trivial single-layer networks. Multilayer networks have become popular in neuroscience due to their advantage to integrate different sources of information. Here, Ι will focus on the multi-frequency multilayer functional connectivity analysis on resting-state fMRI (rs-fMRI) recordings. However, constructing a multilayer network depends on selecting multiple pre-processing steps that can affect the final network topology. Here, I analyzed the rs-fMRI dataset from a single human performing scanning over a period of 18 months (84 scans in total), and the rs-fMRI dataset containing 25 subjects with 3 repeat scans. I focused on assessing the reproducibility of multi-frequency multilayer topologies exploring the effect of two filtering methods for extracting frequencies from BOLD activity, three connectivity estimators, with or without a topological filtering scheme, and two spatial scales. Finally, I untangled specific combinations of researchers' choices that yield consistently brain networks with repeatable topologies, giving me the chance to recommend best practices over consistent topologies.

Modeling the dynamics of cerebrovascular reactivity to carbon dioxide in fMRI under task and resting-state conditions

Conventionally, cerebrovascular reactivity (CVR) is estimated as the amplitude of the hemodynamic response to vascular stimuli, most commonly carbon dioxide (CO₂). While the CVR amplitude has established clinical utility, the temporal characteristics of CVR (dCVR) have been increasingly explored and may yield even more pathology-sensitive parameters. This work is motivated by the current need to evaluate the feasibility of dCVR modeling in various experimental conditions. In this work, we present a comparison of several recently published/utilized model-based deconvolution (response estimation) approaches for estimating the CO₂ response function h(t), including maximum a posteriori likelihood (MAP), inverse logit (IL), canonical correlation analysis (CCA), and basis expansion (using Gamma and Laguerre basis sets). To aid the comparison, we devised a novel simulation framework that incorporates a wide range of SNRs, ranging from 10 to -7 dB, representative of both task and resting-state CO₂ changes. In addition, we built ground-truth h(t) into our simulation framework, overcoming the conventional limitation that the true h(t) is unknown. Moreover, to best represent realistic noise found in fMRI scans, we extracted noise from in-vivo resting-state scans. Furthermore, we introduce a simple optimization of the CCA method (CCA_opt) and compare its performance to these existing methods. Our findings suggest that model-based methods can accurately estimate dCVR even amidst high noise (i.e. resting-state), and in a manner that is largely independent of the underlying model assumptions for each method. We also provide a quantitative basis for making methodological choices, based on the desired dCVR parameters, the estimation accuracy and computation time. The BEL method provided the highest accuracy and robustness, followed by the CCA_opt and IL methods. Of the three, the CCA_opt method has the lowest computational requirements. These findings lay the foundation for wider adoption of dCVR estimation in CVR mapping.

fMRI studies evaluating central respiratory control in humans

A plethora of neural centers in the central nervous system control the fundamental respiratory pattern. This control is ensured by neurons that act as pacemakers, modulating activity through chemical control driven by changes in the O2/CO2 balance. Most of the respiratory neural centers are located in the brainstem, but difficult to localize on magnetic resonance imaging (MRI) due to their small size, lack of visually-detectable borders with neighboring areas, and significant physiological noise hampering detection of its activity with functional MRI (fMRI). Yet, several approaches make it possible to study the normal response to different abnormal stimuli or conditions such as CO2 inhalation, induced hypercapnia, volitional apnea, induced hypoxia etc. This review provides a comprehensive overview of the majority of available studies on central respiratory control in humans.

Separating neuronal activity and systemic low-frequency oscillation related BOLD responses at nodes of the default mode network during resting-state fMRI with multiband excitation echo-planar imaging

Functional magnetic resonance imaging (fMRI) evaluates brain activity using blood oxygenation level-dependent (BOLD) contrast. Resting-state fMRI (rsfMRI) examines spontaneous brain function using BOLD in the absence of a task, and the default mode network (DMN) has been identified from that. The DMN is a set of nodes within the brain that appear to be active and in communication when the subject is in an awake resting state. In addition to signal changes related to neural activity, it is thought that the BOLD signal may be affected by systemic low-frequency oscillations (SysLFOs) that are non-neuronal in source and likely propagate throughout the brain to arrive at different regions at different times. However, it may be difficult to distinguish between the response due to neuronal activity and the arrival of a SysLFO in specific regions. Conventional single-shot EPI (Conv) acquisition requires a longish repetition time, but faster image acquisition has recently become possible with multiband excitation EPI (MB). In this study, we evaluated the time-lag between nodes of the DMN using both Conv and MB protocols to determine whether it is possible to distinguish between neuronal activity and SysLFO related responses during rsfMRI. While the Conv protocol data suggested that SysLFOs substantially influence the apparent time-lag of neuronal activity, the MB protocol data implied that the effects of SysLFOs and neuronal activity on the BOLD response may be separated. Using a higher time-resolution acquisition for rsfMRI might help to distinguish neuronal activity induced changes to the BOLD response from those induced by non-neuronal sources.

Human CSF movement influenced by vascular low frequency oscillations and respiration

Cerebrospinal fluid (CSF) movement through the pathways within the central nervous system is of high significance for maintaining normal brain health and function. Low frequency hemodynamics and respiration have been shown to drive CSF in humans independently. Here, we hypothesize that CSF movement may be driven simultaneously (and in synchrony) by both mechanisms and study their independent and coupled effects on CSF movement using novel neck fMRI scans. Caudad CSF movement at the fourth ventricle and hemodynamics of the major neck blood vessels (internal carotid arteries and internal jugular veins) was measured from 11 young, healthy volunteers using novel neck fMRI scans with simultaneous measurement of respiration. Two distinct models of CSF movement (1. Low-frequency hemodynamics and 2. Respiration) and possible coupling between them were investigated. We show that the dynamics of brain fluids can be assessed from the neck by studying the interrelationships between major neck blood vessels and the CSF movement in the fourth ventricle. We also demonstrate that there exists a cross-frequency coupling between these two separable mechanisms. The human CSF system can respond to multiple coupled physiological forces at the same time. This information may help inform the pathological mechanisms behind CSF movement-related disorders.

Systemic physiology augmented functional near-infrared spectroscopy: a powerful approach to study the embodied human brain

In this Outlook paper, we explain why an accurate physiological interpretation of functional near-infrared spectroscopy (fNIRS) neuroimaging signals is facilitated when systemic physiological activity (e.g., cardiorespiratory and autonomic activity) is measured simultaneously by employing systemic physiology augmented functional near-infrared spectroscopy (SPA-fNIRS). The rationale for SPA-fNIRS is twofold: (i) SPA-fNIRS enables a more complete interpretation and understanding of the fNIRS signals measured at the head since they contain components originating from neurovascular coupling and from systemic physiological sources. The systemic physiology signals measured with SPA-fNIRS can be used for regressing out physiological confounding components in fNIRS signals. Misinterpretations can thus be minimized. (ii) SPA-fNIRS enables to study the embodied brain by linking the brain with the physiological state of the entire body, allowing novel insights into their complex interplay. We envisage the SPA-fNIRS approach will become increasingly important in the future.

Performance of Temporal and Spatial Independent Component Analysis in Identifying and Removing Low-Frequency Physiological and Motion Effects in Resting-State fMRI

Effective separation of signal from noise (including physiological processes and head motion) is one of the chief challenges for improving the sensitivity and specificity of resting-state fMRI (rs-fMRI) measurements and has a profound impact when these noise sources vary between populations. Independent component analysis (ICA) is an approach for addressing these challenges. Conventionally, due to the lower amount of temporal than spatial information in rs-fMRI data, spatial ICA (sICA) is the method of choice. However, with recent developments in accelerated fMRI acquisitions, the temporal information is becoming enriched to the point that the temporal ICA (tICA) has become more feasible. This is particularly relevant as physiological processes and motion exhibit very different spatial and temporal characteristics when it comes to rs-fMRI applications, leading us to conduct a comparison of the performance of sICA and tICA in addressing these types of noise. In this study, we embrace the novel practice of using theory (simulations) to guide our interpretation of empirical data. We find empirically that sICA can identify more noise-related signal components than tICA. However, on the merit of functional-connectivity results, we find that while sICA is more adept at reducing whole-brain motion effects, tICA performs better in dealing with physiological effects. These interpretations are corroborated by our simulation results. The overall message of this study is that if ICA denoising is to be used for rs-fMRI, there is merit in considering a hybrid approach in which physiological and motion-related noise are each corrected for using their respective best-suited ICA approach.

Resting-state functional magnetic resonance imaging signal variations in aging: The role of neural activity

Resting‐state functional magnetic resonance imaging (rs‐fMRI) has been extensively used to study brain aging, but the age effect on the frequency content of the rs‐fMRI signal has scarcely been examined. Moreover, the neuronal implications of such age effects and age–sex interaction remain unclear. In this study, we examined the effects of age and sex on the rs‐fMRI signal frequency using the Leipzig mind–brain–body data set. Over a frequency band of up to 0.3 Hz, we found that the rs‐fMRI fluctuation frequency is higher in the older adults, although the fluctuation amplitude is lower. The rs‐fMRI signal frequency is also higher in men than in women. Both age and sex effects on fMRI frequency vary with the frequency band examined but are not found in the frequency of physiological‐noise components. This higher rs‐fMRI frequency in older adults is not mediated by the electroencephalograph (EEG)‐frequency increase but a likely link between fMRI signal frequency and EEG entropy, which vary with age and sex. Additionally, in different rs‐fMRI frequency bands, the fMRI‐EEG amplitude ratio is higher in young adults. This is the first study to investigate the neuronal contribution to age and sex effects in the frequency dimension of the rs‐fMRI signal and may lead to the development of new, frequency‐based rs‐fMRI metrics. Our study demonstrates that Fourier analysis of the fMRI signal can reveal novel information about aging. Furthermore, fMRI and EEG signals reflect different aspects of age‐ and sex‐related brain differences, but the signal frequency and complexity, instead of amplitude, may hold their link.

Structure can predict function in the human brain: a graph neural network deep learning model of functional connectivity and centrality based on structural connectivity

Although functional connectivity and associated graph theory measures (e.g., centrality; how centrally important to the network a region is) are widely used in brain research, the full extent to which these functional measures are related to the underlying structural connectivity is not yet fully understood. Graph neural network deep learning methods have not yet been applied for this purpose, and offer an ideal model architecture for working with connectivity data given their ability to capture and maintain inherent network structure. Here, we applied this model to predict functional connectivity from structural connectivity in a sample of 998 participants from the Human Connectome Project. Our results showed that the graph neural network accounted for 89% of the variance in mean functional connectivity, 56% of the variance in individual-level functional connectivity, 99% of the variance in mean functional centrality, and 81% of the variance in individual-level functional centrality. These results represent an important finding that functional centrality can be robustly predicted from structural connectivity. Regions of particular importance to the model's performance as determined through lesioning are discussed, whereby regions with higher centrality have a higher impact on model performance. Future research on models of patient, demographic, or behavioural data can also benefit from this graph neural network method as it is ideally-suited for depicting connectivity and centrality in brain networks. These results have set a new benchmark for prediction of functional connectivity from structural connectivity, and models like this may ultimately lead to a way to predict functional connectivity in individuals who are unable to do fMRI tasks (e.g., non-responsive patients).

Lost Dynamics and the Dynamics of Loss: Longitudinal Compression of Brain Signal Variability is Coupled with Declines in Functional Integration and Cognitive Performance

Reduced moment-to-moment blood oxygen level-dependent (BOLD) signal variability has been consistently linked to advanced age and poorer cognitive performance, showing potential as a functional marker of brain aging. To date, however, this promise has rested exclusively on cross-sectional comparisons. In a sample of 74 healthy adults, we provide the first longitudinal evidence linking individual differences in BOLD variability, age, and performance across multiple cognitive domains over an average period of 2.5 years. As expected, those expressing greater loss of BOLD variability also exhibited greater decline in cognition. The fronto-striato-thalamic system emerged as a core neural substrate for these change-change associations. Preservation of signal variability within regions of the fronto-striato-thalamic system also cohered with preservation of functional integration across regions of this system, suggesting that longitudinal maintenance of "local" dynamics may require across-region communication. We therefore propose this neural system as a primary target in future longitudinal studies on the neural substrates of cognitive aging. Given that longitudinal change-change associations between brain and cognition are notoriously difficult to detect, the presence of such an association within a relatively short follow-up period bolsters the promise of brain signal variability as a viable, experimentally sensitive probe for studying individual differences in human cognitive aging.

Effects of Thresholding on Voxel-Wise Correspondence of Breath-Hold and Resting-State Maps of Cerebrovascular Reactivity

Functional magnetic resonance imaging for presurgical brain mapping enables neurosurgeons to identify viable tissue near a site of operable pathology which might be at risk of surgery-induced damage. However, focal brain pathology (e.g., tumors) may selectively disrupt neurovascular coupling while leaving the underlying neurons functionally intact. Such neurovascular uncoupling can result in false negatives on brain activation maps thereby compromising their use for surgical planning. One way to detect potential neurovascular uncoupling is to map cerebrovascular reactivity using either an active breath-hold challenge or a passive resting-state scan. The equivalence of these two methods has yet to be fully established, especially at a voxel level of resolution. To quantitatively compare breath-hold and resting-state maps of cerebrovascular reactivity, we first identified threshold settings that optimized coverage of gray matter while minimizing false responses in white matter. When so optimized, the resting-state metric had moderately better gray matter coverage and specificity. We then assessed the spatial correspondence between the two metrics within cortical gray matter, again, across a wide range of thresholds. Optimal spatial correspondence was strongly dependent on threshold settings which if improperly set tended to produce statistically biased maps. When optimized, the two CVR maps did have moderately good correspondence with each other (mean accuracy of 73.6%). Our results show that while the breath-hold and resting-state maps may appear qualitatively similar they are not quantitatively identical at a voxel level of resolution.

Physiological noise modeling in fMRI based on the pulsatile component of photoplethysmograph

The blood oxygenation level-dependent (BOLD) contrast mechanism allows the noninvasive monitoring of changes in deoxyhemoglobin content. As such, it is commonly used in functional magnetic resonance imaging (fMRI) to study brain activity since levels of deoxyhemoglobin are indirectly related to local neuronal activity through neurovascular coupling mechanisms. However, the BOLD signal is severely affected by physiological processes as well as motion. Due to this, several noise correction techniques have been developed to correct for the associated confounds. The present study focuses on cardiac pulsatility fMRI confounds, aiming to refine model-based techniques that utilize the photoplethysmograph (PPG) signal. Specifically, we propose a new technique based on convolution filtering, termed cardiac pulsatility model (CPM) and compare its performance with the cardiac-related RETROICOR (Card-RETROICOR), which is a technique commonly used to model fMRI fluctuations due to cardiac pulsatility. Further, we investigate whether variations in the amplitude of the PPG pulses (PPG-Amp) covary with variations in amplitude of pulse-related fMRI fluctuations, as well as with the systemic low frequency oscillations (SLFOs) component of the fMRI global signal (GS - defined as the mean signal across all gray matter voxels). Capitalizing on 3T fMRI data from the Human Connectome Project, CPM was found to explain a significantly larger fraction of the fMRI signal variance compared to Card-RETROICOR, particularly for subjects with larger heart rate variability during the scan. The amplitude of the fMRI pulse-related fluctuations did not covary with PPG-Amp; however, PPG-Amp explained significant variance in the GS that was not attributed to variations in heart rate or breathing patterns. Our results suggest that the proposed approach can model high-frequency fluctuations due to pulsation as well as low-frequency physiological fluctuations more accurately compared to model-based techniques commonly employed in fMRI studies.

Physiological and motion signatures in static and time-varying functional connectivity and their subject identifiability

Human brain connectivity yields significant potential as a noninvasive biomarker. Several studies have used fMRI-based connectivity fingerprinting to characterize individual patterns of brain activity. However, it is not clear whether these patterns mainly reflect neural activity or the effect of physiological and motion processes. To answer this question, we capitalize on a large data sample from the Human Connectome Project and rigorously investigate the contribution of the aforementioned processes on functional connectivity (FC) and time-varying FC, as well as their contribution to subject identifiability. We find that head motion, as well as heart rate and breathing fluctuations, induce artifactual connectivity within distinct resting-state networks and that they correlate with recurrent patterns in time-varying FC. Even though the spatiotemporal signatures of these processes yield above-chance levels in subject identifiability, removing their effects at the preprocessing stage improves identifiability, suggesting a neural component underpinning the inter-individual differences in connectivity.

A novel method of quantifying hemodynamic delays to improve hemodynamic response, and CVR estimates in CO2 challenge fMRI

Elevated carbon dioxide (CO2) in breathing air is widely used as a vasoactive stimulus to assess cerebrovascular functions under hypercapnia (i.e., "stress test" for the brain). Blood-oxygen-level-dependent (BOLD) is a contrast mechanism used in functional magnetic resonance imaging (fMRI). BOLD is used to study CO2-induced cerebrovascular reactivity (CVR), which is defined as the voxel-wise percentage BOLD signal change per mmHg change in the arterial partial pressure of CO2 (PaCO2). Besides the CVR, two additional important parameters reflecting the cerebrovascular functions are the arrival time of arterial CO2 at each voxel, and the waveform of the local BOLD signal. In this study, we developed a novel analytical method to accurately calculate the arrival time of elevated CO2 at each voxel using the systemic low frequency oscillations (sLFO: 0.01-0.1 Hz) extracted from the CO2 challenge data. In addition, 26 candidate hemodynamic response functions (HRF) were used to quantitatively describe the temporal brain reactions to a CO2 stimulus. We demonstrated that our approach improved the traditional method by allowing us to accurately map three perfusion-related parameters: the relative arrival time of blood, the hemodynamic response function, and CVR during a CO2 challenge.

A novel approach for assessing hypoperfusion in stroke using spatial independent component analysis of resting-state fMRI

Individualized treatment of acute stroke depends on the timely detection of ischemia and potentially salvageable tissue in the brain. Using functional MRI (fMRI), it is possible to characterize cerebral blood flow from blood‐oxygen‐level‐dependent (BOLD) signals without the administration of exogenous contrast agents. In this study, we applied spatial independent component analysis to resting‐state fMRI data of 37 stroke patients scanned within 24 hr of symptom onset, 17 of whom received follow‐up scans the next day. Our analysis revealed “Hypoperfusion spatially‐Independent Components” (HICs) whose spatial patterns of BOLD signal resembled regions of delayed perfusion depicted by dynamic susceptibility contrast MRI. These HICs were detected even in the presence of excessive patient motion, and disappeared following successful tissue reperfusion. The unique spatial and temporal features of HICs allowed them to be distinguished with high accuracy from other components in a user‐independent manner (area under the curve = 0.93, balanced accuracy = 0.90, sensitivity = 1.00, and specificity = 0.85). Our study therefore presents a new, noninvasive method for assessing blood flow in acute stroke that minimizes interpretative subjectivity and is robust to severe patient motion.

The effects of age on resting-state BOLD signal variability is explained by cardiovascular and cerebrovascular factors

Accurate identification of brain function is necessary to understand neurocognitive aging, and thereby promote health and well-being. Many studies of neurocognitive aging have investigated brain function with the blood-oxygen level-dependent (BOLD) signal measured by functional magnetic resonance imaging. However, the BOLD signal is a composite of neural and vascular signals, which are differentially affected by aging. It is, therefore, essential to distinguish the age effects on vascular versus neural function. The BOLD signal variability at rest (known as resting state fluctuation amplitude, RSFA), is a safe, scalable, and robust means to calibrate vascular responsivity, as an alternative to breath-holding and hypercapnia. However, the use of RSFA for normalization of BOLD imaging assumes that age differences in RSFA reflecting only vascular factors, rather than age-related differences in neural function (activity) or neuronal loss (atrophy). Previous studies indicate that two vascular factors, cardiovascular health (CVH) and cerebrovascular function, are insufficient when used alone to fully explain age-related differences in RSFA. It remains possible that their joint consideration is required to fully capture age differences in RSFA. We tested the hypothesis that RSFA no longer varies with age after adjusting for a combination of cardiovascular and cerebrovascular measures. We also tested the hypothesis that RSFA variation with age is not associated with atrophy. We used data from the population-based, lifespan Cam-CAN cohort. After controlling for cardiovascular and cerebrovascular estimates alone, the residual variance in RSFA across individuals was significantly associated with age. However, when controlling for both cardiovascular and cerebrovascular estimates, the variance in RSFA was no longer associated with age. Grey matter volumes did not explain age differences in RSFA, after controlling for CVH. The results were consistent between voxel-level analysis and independent component analysis. Our findings indicate that cardiovascular and cerebrovascular signals are together sufficient predictors of age differences in RSFA. We suggest that RSFA can be used to separate vascular from neuronal factors, to characterize neurocognitive aging. We discuss the implications and make recommendations for the use of RSFA in the research of aging.

A practical modification to a resting state fMRI protocol for improved characterization of cerebrovascular function

Cerebrovascular reactivity (CVR), defined here as the Blood Oxygenation Level Dependent (BOLD) response to a CO₂ pressure change, is a useful metric of cerebrovascular function. Both the amplitude and the timing (hemodynamic lag) of the CVR response can bring insight into the nature of a cerebrovascular pathology and aid in understanding noise confounds when using functional Magnetic Resonance Imaging (fMRI) to study neural activity. This research assessed a practical modification to a typical resting-state fMRI protocol, to improve the characterization of cerebrovascular function. In 9 healthy subjects, we modelled CVR and lag in three resting-state data segments, and in data segments which added a 2–3 minute breathing task to the start of a resting-state segment. Two different breathing tasks were used to induce fluctuations in arterial CO₂ pressure: a breath-hold task to induce hypercapnia (CO₂ increase) and a cued deep breathing task to induce hypocapnia (CO₂ decrease). Our analysis produced voxel-wise estimates of the amplitude (CVR) and timing (lag) of the BOLD-fMRI response to CO₂ by systematically shifting the CO₂ regressor in time to optimize the model fit. This optimization inherently increases gray matter CVR values and fit statistics. The inclusion of a simple breathing task, compared to a resting-state scan only, increases the number of voxels in the brain that have a significant relationship between CO₂ and BOLD-fMRI signals, and improves our confidence in the plausibility of voxel-wise CVR and hemodynamic lag estimates. We demonstrate the clinical utility and feasibility of this protocol in an incidental finding of Moyamoya disease, and explore the possibilities and challenges of using this protocol in younger populations. This hybrid protocol has direct applications for CVR mapping in both research and clinical settings and wider applications for fMRI denoising and interpretation.

The Effect of CO2 on Resting-State Functional Connectivity: Isocapnia vs. Poikilocapnia

The normal variability in breath size and frequency results in breath-to-breath variability of end-tidal PCO₂ (P_ETCO₂), the measured variable, and arterial partial pressure of carbon dioxide (PaCO₂), the independent variable affecting cerebral blood flow (CBF). This study examines the effect of variability in PaCO₂ on the pattern of resting-state functional MRI (rs-fMRI) connectivity. A region of interest (ROI)-to-ROI and Seed-to-Voxel first-level bivariate correlation, hemodynamic response function (hrf)-weighted analysis for measuring rs-fMRI connectivity was performed during two resting-state conditions: (a) normal breathing associated with breath-to-breath variation in PaCO₂ (poikilocapnia), and (b) normal breathing with breath-to-breath variability of P_ETCO₂ dampened using sequential rebreathing (isocapnia). End-tidal PCO₂ (P_ETCO₂) was used as a measurable surrogate for fluctuations of PaCO₂. During poikilocapnia, enhanced functional connections were found between the cerebellum and inferior frontal and supramarginal gyrus (SG), visual cortex and occipital fusiform gyrus; and between the primary visual network (PVN) and the hippocampal formation. During isocapnia, these associations were not seen, rather enhanced functional connections were identified in the corticostriatal pathway between the putamen and intracalacarine cortex, supracalcarine cortex (SCC), and precuneus cortex. We conclude that vascular responses to variations in P_ETCO₂, account for at least some of the observed resting state synchronization of blood oxygenation level-dependent (BOLD) signals.

Respiratory, cardiac, EEG, BOLD signals and functional connectivity over multiple microsleep episodes

Falling asleep is common in fMRI studies. By using long eyelid closures to detect microsleep onset, we showed that the onset and termination of short sleep episodes invokes a systematic sequence of BOLD signal changes that are large, widespread, and consistent across different microsleep durations. The signal changes are intimately intertwined with shifts in respiration and heart rate, indicating that autonomic contributions are integral to the brain physiology evaluated using fMRI and cannot be simply treated as nuisance signals. Additionally, resting state functional connectivity (RSFC) was altered in accord with the frequency of falling asleep and in a manner that global signal regression does not eliminate. Our findings point to the need to develop a consensus among neuroscientists using fMRI on how to deal with microsleep intrusions. SIGNIFICANCE STATEMENT: Sleep, breathing and cardiac action are influenced by common brainstem nuclei. We show that falling asleep and awakening are associated with a sequence of BOLD signal changes that are large, widespread and consistent across varied durations of sleep onset and awakening. These signal changes follow closely those associated with deceleration and acceleration of respiration and heart rate, calling into question the separation of the latter signals as 'noise' when the frequency of falling asleep, which is commonplace in RSFC studies, correlates with the extent of RSFC perturbation. Autonomic and central nervous system contributions to BOLD signal have to be jointly considered when interpreting fMRI and RSFC studies.

Inhibition of emotions in healthy aging: age-related differences in brain network connectivity

INTRODUCTION

Successful inhibition of distracting emotions is important for preserving well-being and daily functioning. There is conflicting evidence regarding the impact of healthy aging on emotional inhibition, and possible age-related alterations in the neuronal underpinnings of emotional interference processing are unexplored.

METHODS

Thirty younger (mean age 26 years; 15 women) and 30 older (mean age 71 years; 13 women) healthy adults performed a face-word emotional Stroop task while undergoing functional magnetic resonance imaging of the brain. A resting-state scan was acquired for calculating the amplitude of low-frequency fluctuations as an estimate of vascular reactivity. Comparisons of brain activation during the task were assessed in a whole-brain, voxel-wise analysis, contrasting congruent, and incongruent conditions. The canonical regions of the frontoparietal, salience, dorsal attention, and default mode networks were used as seed regions for assessing functional connectivity within and between large-scale brain networks. Task performance was evaluated using response accuracy and response time.

RESULTS

The older adults had longer response times and lower task accuracy than the younger adults, but the emotional interference effect was not significantly different between the groups. Whole-brain analysis revealed no significant age-related differences in brain activation patterns. Rescaling the data for estimated variability in vascular reactivity did not affect the results. In older adults, there was relatively stronger functional connectivity with the default mode network, the sensorimotor network, and the dorsal attention network for the frontoparietal and salience network seeds during the task. Conversely, younger adults had relatively stronger connections within and between the frontoparietal and salience networks.

CONCLUSION

In this first fMRI study of emotional Stroop interference in older and younger adults, we found that the emotional interference effect was unchanged in healthy aging and replicated the finding from non-emotional task studies that older adults have greater between-network and less within-network connectivity compared to younger adults.

Vascular origins of low-frequency oscillations in the cerebrospinal fluid signal in resting-state fMRI: Interpretation using photoplethysmography

In vivo mapping of cerebrovascular oscillations in the 0.05–0.15 Hz remains difficult. Oscillations in the cerebrospinal fluid (CSF) represent a possible avenue for noninvasively tracking these oscillations using resting‐state functional MRI (rs‐fMRI), and have been used to correct for vascular oscillations in rs‐fMRI functional connectivity. However, the relationship between low‐frequency CSF and vascular oscillations remains unclear. In this study, we investigate this relationship using fast simultaneous rs‐fMRI and photoplethysmogram (PPG), examining the 0.1 Hz PPG signal, heart‐rate variability (HRV), pulse‐intensity ratio (PIR), and the second derivative of the PPG (SDPPG). The main findings of this study are: (a) signals in different CSF regions are not equivalent in their associations with vascular and tissue rs‐fMRI signals; (b) the PPG signal is maximally coherent with the arterial and CSF signals at the cardiac frequency, but coherent with brain tissue at ~0.2 Hz; (c) PIR is maximally coherent with the CSF signal near 0.03 Hz; and (d) PPG‐related vascular oscillations only contribute to ~15% of the CSF (and arterial) signal in rs‐fMRI. These findings caution against averaging all CSF regions when extracting physiological nuisance regressors in rs‐fMRI applications, and indicate the drivers of the CSF signal are more than simply cardiac. Our study is an initial attempt at the refinement and standardization of how the CSF signal in rs‐fMRI can be used and interpreted. It also paves the way for using rs‐fMRI in the CSF as a potential tool for tracking cerebrovascular health through, for instance, the potential relationship between PIR and the CSF signal.

The neuronal associations of respiratory-volume variability in the resting state

The desire to enhance the sensitivity and specificity of resting-state (rs-fMRI) measures has prompted substantial recent research into removing noise components. Chief among contributions to noise in rs-fMRI are physiological processes, and the neuronal implications of respiratory-volume variability (RVT), a main rs-fMRI-relevant physiological process, is incompletely understood. The potential implications of RVT in modulating and being modulated by autonomic nervous regulation, has yet to be fully understood by the rs-fMRI community. In this work, we use high-density electroencephalography (EEG) along with simultaneously acquired RVT recordings to help address this question. We hypothesize that (1) there is a significant relationship between EEG and RVT in multiple EEG bands, and (2) that this relationship varies by brain region. Our results confirm our first hypothesis, although all brain regions are shown to be equally implicated in RVT-related EEG-signal fluctuations. The lag between RVT and EEG is consistent with previously reported values. However, an interesting finding is related to the polarity of the correlation between RVT and EEG. Our results reveal potentially two main regimes of EEG-RVT association, one in which EEG leads RVT with a positive association between the two, and one in which RVT leads EEG but with a negative association between the two. We propose that these two patterns can be interpreted differently in terms of the involvement of higher cognition. These results further suggest that treating RVT simply as noise is likely a questionable practice, and that more work is needed to avoid discarding cognitively relevant information when performing physiological correction rs-fMRI.

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.

Separating vascular and neuronal effects of age on fMRI BOLD signals

Accurate identification of brain function is necessary to understand the neurobiology of cognitive ageing, and thereby promote well-being across the lifespan. A common tool used to investigate neurocognitive ageing is functional magnetic resonance imaging (fMRI). However, although fMRI data are often interpreted in terms of neuronal activity, the blood oxygenation level-dependent (BOLD) signal measured by fMRI includes contributions of both vascular and neuronal factors, which change differentially with age. While some studies investigate vascular ageing factors, the results of these studies are not well known within the field of neurocognitive ageing and therefore vascular confounds in neurocognitive fMRI studies are common. Despite over 10 000 BOLD-fMRI papers on ageing, fewer than 20 have applied techniques to correct for vascular effects. However, neurovascular ageing is not only a confound in fMRI, but an important feature in its own right, to be assessed alongside measures of neuronal ageing. We review current approaches to dissociate neuronal and vascular components of BOLD-fMRI of regional activity and functional connectivity. We highlight emerging evidence that vascular mechanisms in the brain do not simply control blood flow to support the metabolic needs of neurons, but form complex neurovascular interactions that influence neuronal function in health and disease. This article is part of the theme issue 'Key relationships between non-invasive functional neuroimaging and the underlying neuronal activity'.

Impact of physiological noise in characterizing the functional MRI default-mode network in Alzheimer's disease

The functional connectivity of the default-mode network (DMN) monitored by functional magnetic resonance imaging (fMRI) in Alzheimer's disease (AD) patients has been found weaker than that in healthy participants. Since breathing and heart beating can cause fluctuations in the fMRI signal, these physiological activities may affect the fMRI data differently between AD patients and healthy participants. We collected resting-state fMRI data from AD patients and age-matched healthy participants. With concurrent cardiac and respiratory recordings, we estimated both physiological responses phase-locked and non-phase-locked to heart beating and breathing. We found that the cardiac and respiratory physiological responses in AD patients were 3.00 ± 0.51 s and 3.96 ± 0.52 s later (both p < 0.0001) than those in healthy participants, respectively. After correcting the physiological noise in the resting-state fMRI data by population-specific physiological response functions, the DMN estimated by seed-correlation was more localized to the seed region. The DMN difference between AD patients and healthy controls became insignificant after suppressing physiological noise. Our results indicate the importance of controlling physiological noise in the resting-state fMRI analysis to obtain clinically related characterizations in AD.

Reconstruction of respiratory variation signals from fMRI data

Functional MRI signals can be heavily influenced by systemic physiological processes in addition to local neural activity. For example, widespread hemodynamic fluctuations across the brain have been found to correlate with natural, low-frequency variations in the depth and rate of breathing over time. Acquiring peripheral measures of respiration during fMRI scanning not only allows for modeling such effects in fMRI analysis, but also provides valuable information for interrogating brain-body physiology. However, physiological recordings are frequently unavailable or have insufficient quality. Here, we propose a computational technique for reconstructing continuous low-frequency respiration volume (RV) fluctuations from fMRI data alone. We evaluate the performance of this approach across different fMRI preprocessing strategies. Further, we demonstrate that the predicted RV signals can account for similar patterns of temporal variation in resting-state fMRI data compared to measured RV fluctuations. These findings indicate that fluctuations in respiration volume can be extracted from fMRI alone, in the common scenario of missing or corrupted respiration recordings. The results have implications for enriching a large volume of existing fMRI datasets through retrospective addition of respiratory variations information.

Dynamic brain-body coupling of breath-by-breath O2-CO2 exchange ratio with resting state cerebral hemodynamic fluctuations

BACKGROUND

The origin of low frequency cerebral hemodynamic fluctuations (CHF) in the resting state remains unknown. Breath-by breath O2-CO2 exchange ratio (bER) has been reported to correlate with the cerebrovascular response to brief breath hold challenge at the frequency range of 0.008-0.03Hz in healthy adults. bER is defined as the ratio of the change in the partial pressure of oxygen (ΔPO2) to that of carbon dioxide (ΔPCO2) between end inspiration and end expiration. In this study, we aimed to investigate the contribution of respiratory gas exchange (RGE) metrics (bER, ΔPO2 and ΔPCO2) to low frequency CHF during spontaneous breathing.

METHODS

Twenty-two healthy adults were included. We used transcranial Doppler sonography to evaluate CHF by measuring the changes in cerebral blood flow velocity (ΔCBFv) in bilateral middle cerebral arteries. The regional CHF were mapped with blood oxygenation level dependent (ΔBOLD) signal changes using functional magnetic resonance imaging. Temporal features and frequency characteristics of RGE metrics during spontaneous breathing were examined, and the simultaneous measurements of RGE metrics and CHF (ΔCBFv and ΔBOLD) were studied for their correlation.

RESULTS

We found that the time courses of ΔPO2 and ΔPCO2 were interdependent but not redundant. The oscillations of RGE metrics were coherent with resting state CHF at the frequency range of 0.008-0.03Hz. Both bER and ΔPO2 were superior to ΔPCO2 in association with CHF while CHF could correlate more strongly with bER than with ΔPO2 in some brain regions. Brain regions with the strongest coupling between bER and ΔBOLD overlapped with many areas of default mode network including precuneus and posterior cingulate.

CONCLUSION

Although the physiological mechanisms underlying the strong correlation between bER and CHF are unclear, our findings suggest the contribution of bER to low frequency resting state CHF, providing a novel insight of brain-body interaction via CHF and oscillations of RGE metrics.

Heart Rate and Respiration Affect the Functional Connectivity of Default Mode Network in Resting-State Functional Magnetic Resonance Imaging

A growing number of brain imaging studies show functional connectivity (FC) between regions during emotional and cognitive tasks in humans. However, emotions are accompanied by changes in physiological parameters such as heart rate and respiration. These changes may affect blood oxygen level-dependent signals, as well as connectivity between brain areas. This study aimed to clarify the effects of physiological noise on the connectivity between areas related to the default mode network using resting-state functional magnetic resonance imaging (rs-fMRI). Healthy adult volunteers (age range: 19–51 years, mean age: 26.9 ± 9.1 years, 8 males and 8 females) underwent rs-fMRI for 10 min using a clinical 3T scanner (MAGNETOM Trio A Tim System, Siemens) with simultaneously recorded respiration and cardiac output. Physiological noise signals were subsequently removed from the acquired fMRI data using the DRIFTER toolbox. Image processing and analysis of the FC between areas related to the default mode network were performed using DPARSF. Network-Based Statistic (NBS) analysis of the functional connectome of the DMN and DMN-related area was used to perform three groups of comparison: without physiological noise correction, with cardiac noise correction, and with cardiac and respiratory noise correction. NBS analysis identified 36 networks with significant differences in three conditions in FC matrices. Post hoc comparison showed no differences between the three conditions, indicating that all three had the same networks. Among the 36 networks, strength of FC of 8 networks was modified under physiological noise correction. Connectivity between left and right anterior medial frontal regions increased strength of connectivity. These areas are located on the medial cerebral hemisphere, close to the sagittal sinus and arteries in the cerebral hemispheres, suggesting that medial frontal areas may be sensitive to cardiac rhythm close to arteries. The other networks observed temporal regions and showed a decrease in their connectivity strength by removing physiological noise, indicating that physiological noise, especially respiration, may be sensitive to BOLD signal in the temporal regions during resting state. Temporal lobe was highly correlated with anxiety-related respiration changes (Masaoka and Homma, 2000), speech processing, and respiratory sensation. These factors may affect the rs-fMRI signaling sensitivity.

Resting-state "physiological networks"

Slow changes in systemic brain physiology can elicit large fluctuations in fMRI time series, which manifest as structured spatial patterns of temporal correlations between distant brain regions. Here, we investigated whether such "physiological networks"-sets of segregated brain regions that exhibit similar responses following slow changes in systemic physiology-resemble patterns associated with large-scale networks typically attributed to remotely synchronized neuronal activity. By analyzing a large group of subjects from the 3T Human Connectome Project (HCP) database, we demonstrate brain-wide and noticeably heterogenous dynamics tightly coupled to either respiratory variation or heart rate changes. We show, using synthesized data generated from physiological recordings across subjects, that these physiologically-coupled fluctuations alone can produce networks that strongly resemble previously reported resting-state networks, suggesting that, in some cases, the "physiological networks" seem to mimic the neuronal networks. Further, we show that such physiologically-relevant connectivity estimates appear to dominate the overall connectivity observations in multiple HCP subjects, and that this apparent "physiological connectivity" cannot be removed by the use of a single nuisance regressor for the entire brain (such as global signal regression) due to the clear regional heterogeneity of the physiologically-coupled responses. Our results challenge previous notions that physiological confounds are either localized to large veins or globally coherent across the cortex, therefore emphasizing the necessity to consider potential physiological contributions in fMRI-based functional connectivity studies. The rich spatiotemporal patterns carried by such "physiological" dynamics also suggest great potential for clinical biomarkers that are complementary to large-scale neuronal networks.

Cerebrovascular reactivity mapping using intermittent breath modulation

Cerebrovascular reactivity (CVR), an index of brain vessel's dilatory capacity, is typically measured using hypercapnic gas inhalation or breath-holding as a vasoactive challenge. However, these methods require considerable subject cooperation and could be challenging in clinical studies. More recently, there have been attempts to use resting-state BOLD data to map CVR by utilizing spontaneous changes in breathing pattern. However, in subjects who have small fluctuations in their spontaneous breathing pattern, the CVR results could be noisy and unreliable. In this study, we aim to develop a new method for CVR mapping that does not require gas-inhalation yet provides substantially higher sensitivity than resting-state CVR mapping. This new method is largely based on resting-state scan, but introduces intermittent modulation of breathing pattern in the subject to enhance fluctuations in their end-tidal CO2 (EtCO2) level. Here we examined the comfort level, sensitivity, and accuracy of this method in two studies. First, in 8 healthy young subjects, we developed the intermittent breath-modulation method using two different modulation frequencies, 6 ​s per breath and 12 ​s per breath, respectively, and compared the results to three existing CVR methods, specifically hypercapnic gas inhalation, breath-holding, and resting-state. Our results showed that the comfort level of the 6-s breath-modulation method was significantly higher than breath-holding (p ​= ​0.007) and CO2-inhalation (p ​= ​0.015) methods, while not different from the resting-state, i.e. free breathing method (p ​= ​0.52). When comparing the sensitivity of CVR methods, the breath-modulation methods revealed higher Z-statistics compared to the resting-state scan (p ​< ​0.008) and was comparable to breath-holding results. Next, we tested the feasibility of breath-modulation CVR mapping (6 ​s per breath) in 21 cognitively normal elderly participants and compared quantitative CVR values to that obtained with the CO2-inhalation method. Whole-brain CVR was found to be 0.150 ​± ​0.055 and 0.154 ​± ​0.032 %ΔBOLD/mmHg for the breath-modulation and CO2-inhalation method, respectively, with a significant correlation between them (y ​= ​0.97x, p ​= ​0.007). CVR mapping with intermittent breath modulation may be a useful method that combines the advantages of resting-state and CO2-inhalation based approaches.

Frontal cerebral oxygenation asymmetry: intersubject variability and dependence on systemic physiology, season, and time of day

Significance: Our study reveals that frontal cerebral oxygenation asymmetry (FCOA), i.e. a difference in the oxygenation between the right and left prefrontal cortex (PFC), is a real phenomenon in healthy human subjects at rest. Aim: To investigate FCOA, we performed a study with 134 healthy right-handed subjects with the systemic physiology augmented functional near infrared spectroscopy (SPA-fNIRS) approach. Approach: Subjects were measured 2 to 4 times on different days resulting in an unprecedented number of 518 single measurements of the absolute values of tissue oxygen saturation (StO 2 ) and total hemoglobin concentration ([tHb]) of the right and left PFC. Measurements were performed with frequency-domain functional near-infrared spectroscopy. In addition, the cardiorespiratory parameters were measured simultaneously. Results: We found that (i) subjects showed an FCOA (higherStO 2 on the right PFC), but not for tHb; (ii) intrasubject variability was excellent for bothStO 2 and tHb, and fair for FCOA; (iii) StO 2 correlated significantly with bloodCO 2 concentration, [tHb] with heart rate, respiration rate (RR), and the pulse-respiration quotient (PRQ), and FCOA with RR and PRQ; (iv) FCOA andStO 2 were dependent on season and time of day, respectively; (v) FCOA was negatively correlated with the room temperature; and (vi) StO 2 and tHb were not correlated with the subjects mood but with their chronotype, whereas FCOA was not dependent on the chronotype. Conclusion: Our study demonstrates that FCOA is real, and it provides unique insights into this remarkable phenomenon.