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

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.

Estimation of cerebrovascular reactivity amplitude and lag using breath-holding fMRI and the global BOLD signal: Application in diabetes and hypertension

In this work, we demonstrate a data-driven approach for estimating cerebrovascular reactivity (CVR) amplitude and lag from breathhold (BH) fMRI data alone. Our approach employs a frequency-domain approach that is independent of external recordings. CVR amplitude is estimated from the BOLD frequency spectrum and CVR lag is estimated from the Fourier phase using the global-mean BOLD signal as reference. Unlike referencing to external recordings, these lags are specific to the brain. We demonstrated our method in detecting regional CVR amplitude and lag differences across healthy (CTL), hypertensive (HT) and hypertension-plus-type-2-diabetes (HT + DM) groups of similar ages and sex ratios, with a total N of 49. We found CVR amplitude to be significantly higher in CTL compared to HT + DM, with minimal difference between CTL and HT. Also, voxelwise CVR lag estimated in the Fourier domain is a more sensitive marker of vascular dysfunction than CVR amplitude. CVR lag in HT is significantly shorter than in CTL, with minimal difference between CTL and HT + DM. Our results support the importance of joint CVR amplitude and lag assessments in clinical applications.

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.

Macrovascular blood flow and microvascular cerebrovascular reactivity are regionally coupled in adolescence

Cerebrovascular imaging assessments are particularly challenging in adolescent cohorts, where not all modalities are appropriate, and rapid brain maturation alters hemodynamics at both macro- and microvascular scales. In a preliminary sample of healthy adolescents (n = 12, 8-25 years), we investigated relationships between 4D flow MRI-derived blood velocity and blood flow in bilateral anterior, middle, and posterior cerebral arteries and BOLD cerebrovascular reactivity (CVR) in associated vascular territories. As hypothesized, higher velocities in large arteries are associated with an earlier response to a vasodilatory stimulus (cerebrovascular reactivity delay) in the downstream territory. Higher blood flow through these arteries is associated with a larger BOLD response to a vasodilatory stimulus (cerebrovascular reactivity amplitude) in the associated territory. These trends are consistent in a case study of adult moyamoya disease. In our small adolescent cohort, macrovascular-microvascular relationships for velocity/delay and flow/CVR change with age, though underlying mechanisms are unclear. Our work emphasizes the need to better characterize this key stage of human brain development, when cerebrovascular hemodynamics are changing, and standard imaging methods offer limited insight into these processes. We provide important normative data for future comparisons in pathology, where combining macro- and microvascular assessments may better help us prevent, stratify, and treat cerebrovascular disease.

Macrovascular blood flow and microvascular cerebrovascular reactivity are regionally coupled in adolescence

Cerebrovascular imaging assessments are particularly challenging in adolescent cohorts, where not all modalities are appropriate, and rapid brain maturation alters hemodynamics at both macro- and microvascular scales. In a preliminary sample of healthy adolescents (n=12, 8-25 years), we investigated relationships between 4D flow MRI-derived blood velocity and blood flow in bilateral anterior, middle, and posterior cerebral arteries and BOLD cerebrovascular reactivity in associated vascular territories. As hypothesized, higher velocities in large arteries are associated with an earlier response to a vasodilatory stimulus (cerebrovascular reactivity delay) in the downstream territory. Higher blood flow through these arteries is associated with a larger BOLD response to a vasodilatory stimulus (cerebrovascular reactivity amplitude) in the associated territory. These trends are consistent in a case study of adult moyamoya disease. In our small adolescent cohort, macrovascular-microvascular relationships for velocity/delay and flow/CVR change with age, though underlying mechanisms are unclear. Our work emphasizes the need to better characterize this key stage of human brain development, when cerebrovascular hemodynamics are changing, and standard imaging methods offer limited insight into these processes. We provide important normative data for future comparisons in pathology, where combining macro- and microvascular assessments may better help us prevent, stratify, and treat cerebrovascular disease.

Alteration in early resting‑state functional MRI activity in comatose survivors of cardiac arrest: a prospective cohort study

BACKGROUND

This study aimed to explore the characteristics of abnormal regional resting-state functional magnetic resonance imaging (rs-fMRI) activity in comatose patients in the early period after cardiac arrest (CA), and to investigate their relationships with neurological outcomes. We also explored the correlations between jugular venous oxygen saturation (SjvO2) and rs-fMRI activity in resuscitated comatose patients. We also examined the relationship between the amplitude of the N20-baseline and the rs-fMRI activity within the intracranial conduction pathway of somatosensory evoked potentials (SSEPs).

METHODS

Between January 2021 and January 2024, eligible post-resuscitated patients were screened to undergo fMRI examination. The amplitude of low-frequency fluctuation (ALFF), fractional ALFF (fALFF), and regional homogeneity (ReHo) of rs-fMRI blood oxygenation level-dependent (BOLD) signals were used to characterize regional neural activity. Neurological outcomes were evaluated using the Glasgow-Pittsburgh cerebral performance category (CPC) scale at 3 months after CA.

RESULTS

In total, 20 healthy controls and 31 post-resuscitated patients were enrolled in this study. The rs-fMRI activity of resuscitated patients revealed complex changes, characterized by increased activity in some local brain regions and reduced activity in others compared to healthy controls (P < 0.05). However, the mean ALFF values of the whole brain were significantly greater in CA patients (P = 0.011). Among the clusters of abnormal rs-fMRI activity, the cluster values of ALFF in the left middle temporal gyrus and inferior temporal gyrus and the cluster values of ReHo in the right precentral gyrus, superior frontal gyrus and middle frontal gyrus were strongly correlated with the CPC score (P < 0.001). There was a strong correlation between the mean ALFF and SjvO2 in CA patients (r = 0.910, P < 0.001). The SSEP N20-baseline amplitudes in CA patients were negatively correlated with thalamic rs-fMRI activity (all P < 0.001).

CONCLUSIONS

This study revealed that abnormal rs-fMRI BOLD signals in resuscitated patients showed complex changes, characterized by increased activity in some local brain regions and reduced activity in others. Abnormal BOLD signals were associated with neurological outcomes in resuscitated patients. The mean ALFF values of the whole brain were closely related to SjvO2 levels, and changes in the thalamic BOLD signals correlated with the N20-baseline amplitudes of SSEP responses.

TRIAL REGISTRATION

NCT05966389 (Registered July 27, 2023).

Delayed cerebrovascular reactivity in individuals with spinal cord injury in the right inferior parietal lobe: a breath-hold functional near-infrared spectroscopy study

Cerebrovascular reactivity (CVR) reflects the ability of blood vessels to dilate or constrict in response to a vasoactive stimulus, and allows researchers to assess the brain's vascular health. Individuals with spinal cord injury (SCI) are at an increased risk for autonomic dysfunction in addition to cognitive impairments, which have been linked to a decline in CVR; however, there is currently a lack of brain-imaging studies that investigate how CVR is altered after SCI. In this study, we used a breath-holding hypercapnic stimulus and functional near-infrared spectroscopy (fNIRS) to investigate CVR alterations in individuals with SCI (n = 20, 14M, 6F, mean age = 46.3 ± 10.2 years) as compared to age- and sex-matched able-bodied (AB) controls (n = 25, 19M, 6F, mean age = 43.2 ± 12.28 years). CVR was evaluated by its amplitude and delay components separately by using principal component analysis and cross-correlation analysis, respectively. We observed significantly delayed CVR in the right inferior parietal lobe in individuals with SCI compared to AB controls (linear mixed-effects model, fixed-effects estimate = 6.565, Satterthwaite's t-test, t = 2.663, p = 0.008), while the amplitude of CVR was not significantly different. The average CVR delay in the SCI group in the right inferior parietal lobe was 14.21 s (sd: 6.60 s), and for the AB group, the average delay in the right inferior parietal lobe was 7.08 s (sd: 7.39 s). CVR delays were also associated with the duration since injury in individuals with SCI, in which a longer duration since injury was associated with a shortened delay in CVR in the right inferior parietal region (Pearson's r-correlation, r = -0.59, p = 0.04). This study shows that fNIRS can be used to quantify changes in CVR in individuals with SCI, and may be further used in rehabilitative settings to monitor the cerebrovascular health of individuals with SCI.

MRI mapping of hemodynamics in the human spinal cord

Impaired spinal cord vascular function contributes to numerous neurological pathologies, making it important to be able to noninvasively characterize these changes. Here, we propose a functional magnetic resonance imaging (fMRI)-based method to map spinal cord vascular reactivity (SCVR). We used a hypercapnic breath-holding task, monitored with end-tidal CO2 (PETCO2), to evoke a systemic vasodilatory response during concurrent blood oxygenation level-dependent (BOLD) fMRI. SCVR amplitude and hemodynamic delay were mapped at the group level in 27 healthy participants as proof-of-concept of the approach, and then in two highly-sampled participants to probe feasibility/stability of individual SCVR mapping. Across the group and the highly-sampled individuals, a strong ventral SCVR amplitude was initially observed without accounting for local regional variation in the timing of the vasodilatory response. Shifted breathing traces (PETCO2) were used to account for temporal differences in the vasodilatory response across the spinal cord, producing maps of SCVR delay. These delay maps reveal an earlier ventral and later dorsal response and demonstrate distinct gray matter regions concordant with territories of arterial supply. The SCVR fMRI methods described here enable robust mapping of spatiotemporal hemodynamic properties of the human spinal cord. This noninvasive approach has exciting potential to provide early insight into pathology-driven vascular changes in the cord, which may precede and predict future irreversible tissue damage and guide the treatment of several neurological pathologies involving the spine.

The significance and limited influence of cerebrovascular reactivity on age and sex effects in task- and resting-state brain activity

Functional MRI measures the blood-oxygen-level dependent signals, which provide an indirect measure of neural activity mediated by neurovascular responses. Cerebrovascular reactivity affects both task-induced and resting-state blood-oxygen-level dependent activity and may confound inter-individual effects, such as those related to aging and biological sex. We examined a large dataset containing breath-holding, checkerboard, and resting-state tasks. We used the breath-holding task to measure cerebrovascular reactivity, used the checkerboard task to obtain task-based activations, and quantified resting-state activity with amplitude of low-frequency fluctuations and regional homogeneity. We hypothesized that cerebrovascular reactivity would be correlated with blood-oxygen-level dependent measures and that accounting for these correlations would result in better estimates of age and sex effects. We found that cerebrovascular reactivity was correlated with checkerboard task activations in the visual cortex and with amplitude of low-frequency fluctuations and regional homogeneity in widespread fronto-parietal regions, as well as regions with large vessels. We also found significant age and sex effects in cerebrovascular reactivity, some of which overlapped with those observed in amplitude of low-frequency fluctuations and regional homogeneity. However, correcting for the effects of cerebrovascular reactivity had very limited influence on the estimates of age and sex. Our results highlight the limitations of accounting for cerebrovascular reactivity with the current breath-holding task.

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.

Alteration of the large-scale white-matter functional networks in autism spectrum disorder

Autism spectrum disorder is a neurodevelopmental disorder whose core deficit is social dysfunction. Previous studies have indicated that structural changes in white matter are associated with autism spectrum disorder. However, few studies have explored the alteration of the large-scale white-matter functional networks in autism spectrum disorder. Here, we identified ten white-matter functional networks on resting-state functional magnetic resonance imaging data using the K-means clustering algorithm. Compared with the white matter and white-matter functional network connectivity of the healthy controls group, we found significantly decreased white matter and white-matter functional network connectivity mainly located within the Occipital network, Middle temporo-frontal network, and Deep network in autism spectrum disorder. Compared with healthy controls, findings from white-matter gray-matter functional network connectivity showed the decreased white-matter gray-matter functional network connectivity mainly distributing in the Occipital network and Deep network. Moreover, we compared the spontaneous activity of white-matter functional networks between the two groups. We found that the spontaneous activity of Middle temporo-frontal and Deep network was significantly decreased in autism spectrum disorder. Finally, the correlation analysis showed that the white matter and white-matter functional network connectivity between the Middle temporo-frontal network and others networks and the spontaneous activity of the Deep network were significantly correlated with the Social Responsiveness Scale scores of autism spectrum disorder. Together, our findings indicate that changes in the white-matter functional networks are associated behavioral deficits in autism spectrum disorder.

The Significance and Limited Influence of Cerebrovascular Reactivity on Age and Sex Effects in Task- and Resting-State Brain Activity

Functional MRI (fMRI) measures the blood-oxygen-level dependent (BOLD) signals, which provide an indirect measure of neural activity mediated by neurovascular responses. Cerebrovascular reactivity affects both task-induced and resting-state BOLD activity and may confound inter-individual effects observed in BOLD-based measures, such as those related to aging and biological sex. To investigate this, we examined a large open-access fMRI dataset containing a breath-holding task, checkerboard task, and resting-state scans. We used the breath-holding task to measure cerebrovascular reactivity, used the checkerboard task to obtain task-based activations, and from the resting-state data, we quantified the resting-state amplitude of low-frequency fluctuations (ALFF), and resting-state regional homogeneity (ReHo). We hypothesized that cerebrovascular reactivity would be correlated with BOLD measures and that accounting for these correlations would result in better estimates of age and sex effects. Our analysis showed that cerebrovascular reactivity was correlated with checkerboard task activations in the visual cortex and with ALFF and ReHo in widespread fronto-parietal regions, as well as regions with large vessels. We also found significant age and sex effects in cerebrovascular reactivity, some of which overlapped with those observed in ALFF and ReHo scores. Finally, we demonstrated that correcting for the effects of cerebrovascular reactivity had very limited influence on the estimates of age and sex. Our results highlight the limitations of accounting for cerebrovascular reactivity with the current breath-holding task.