Abstract
Cerebrovascular reactivity (CVR) and cerebral pulsatility (CP) are important indicators of cerebrovascular health, which are associated with physical activity (PA). While sex differences influence the impact of PA on cerebrovascular health, sex-specific effects of PA intensity and dose on CP and CVR remains unknown. This study aimed to evaluate the sex-specific effects of self-reported PA dose and intensity on CVR and CP. The Human Connectome Project - Aging dataset was used, including 626 participants (350 females, 276 males) aged 36–85. The effect of menopausal status was also assessed. Resting state fMRI data was used to estimate both CVR and CP. Weekly self-reported PA was quantified as metabolic equivalent of task. Females presented a unique non-linear relationship between relative CVR and total PA in the cerebral cortex. Females and menopausal subgroups revealed negative linear relationships with total and walking PA in occipital and cingulate regions. Males exhibited negative linear relationships between total and vigorous PA and CVR in parietal and cingulate regions. Postmenopausal females showed greater reductions across more regions in CP than other groups. Overall, males and females appear to benefit from different amounts and intensities of PA, with menopause status influencing the effect of PA on cerebrovascular health.
Keywords
Introduction
Physical inactivity is associated with an increased risk of developing cardiovascular 1 and neurological diseases. 2 More specifically, physical inactivity is associated with greater arterial stiffness, 3 reduced endothelial function 4 and decreased vasodilatory function. 5 Conversely, participating in physical activity (PA) has been shown to enhance vascular health, thereby delaying or preventing disease onset and progression. 6 In the brain, PA has been shown to lead to changes in cerebral blood flow (CBF) and cerebrovascular reactivity (CVR). 7 Importantly, PA may slow or even reverse cognitive decline associated with aging.8,9 Furthermore, improvements in cognition resulting from PA engagement correlate with changes in cerebral hemodynamics, including higher CBF 10 and lower cerebral pulsatility (CP). 11 Thus, PA may promote healthy aging in part through maintenance of cerebrovascular function.
The current PA recommendations from the World Health Organisation are identical for adult males and females 12 across the lifespan. However, sex differences have been observed in the effects of PA on health, including in the brain. For example, sex differences have been identified in the relationship between grey and white matter volumes and PA among older adults. 13 Moreover, older aged males and females exhibit differences in the association between cardiorespiratory fitness and functional connectivity. 14 Although the physiological underpinnings of these sex differences remain unclear, PA appears to offer unique protective mechanisms for each sex. 15 Therefore, it is likely that males and females benefit in terms of cerebrovascular health from different intensity and dose of PA. Assessing the dose and intensity of PA needed to derive benefit for each sex could furthermore help refine future recommendations in the context of a personalized medicine approach to disease prevention.
Sex differences in the prevalence of diseases, including Alzheimer’s disease where two-thirds of patients are female, highlight the need to better understand sex-specifc mechanisms. 16 However, lifelong sex differences alone cannot fully explain the vascular and metabolic changes observed in an aging population; rather, the menopausal transition in females influences these physiological alterations. 17 Menopause has been identified as a significant risk factor for the development of cardiovascular diseases. 18 Nevertheless, PA, independent of the menopausal status, has been shown to reduce cardiovascular risk factors. 19 Additionally, menopause is associated with a decline in physical performance, especially strength and power. 20 However, a higher dose of PA can significantly mitigate these reductions in physical performance. 21 Hence, to understand the effects of PA on female cerebrovascular health, it is crucial to stratify according to the menopausal stage.
Magnetic resonance imaging (MRI) can be used to assess several cerebrovascular biomarkers, including CBF, CVR and CP. Furthermore, CVR and CP are critical markers of vascular brain health, which are affected in small vessel disease. 22 CVR captures the ability of blood vessels to dilate in response to a vasodilatory stimulus, such as carbon dioxide (CO2). CVR is a measure of vascular reserve and has been shown to be more sensitive than CBF to the effects of aging.23,24 The effects of fitness and PA participation on CVR is unclear however, with the extant literature showing both positive,25 –27 negative linear7,28 and inverted U-shaped relationships 29 in older adults. While some of these effects likely arise from methodological differences, it is currently unknown if some of these contradictory effects may be influenced by sex-specific effects of PA on CVR, or in differential effects of menopause. Here, we will address this gap in the literature. Several methods can be used to measure CVR. Functional MRI (fMRI) blood oxygen level dependent (BOLD) CVR is typically measured as the percent change in BOLD signal in response to a gas manipulation challenge or injection of acetazolamide. 30 Recently, a novel technique 31 was introduced to acquire CVR exclusively from the BOLD signal during rest, without any respiratory challenge. This method is based on the ability to detect natural end-tidal CO2 fluctuations from the lower frequencies of the BOLD signal.
CP measures the extent to which blood vessels expand and recoil in response to fluctuations in blood pressure during the cardiac cycle. CP primarily reflects arterial stiffness through the measurement of pulse pressure waves. 32 It therefore represents complementary information to CVR, which represents the vasodilatory capacity. 30 Higher CP is associated with a decline in microvascular integrity, resulting from increased arterial stiffness. 33 Sex differences in the effects of PA on CP have been investigated using transcranial doppler (TCD) and 4D flow MRI. Females were shown to exhibit lower CP in the basilar artery with higher levels of PA, whereas males displayed the opposite relationship in the middle cerebral artery. Thus, there may be important sexual dimorphisms in the relationship between PA and CP. 34 MRI studies of CP have typically used 4D flow in large vessels. 35 CP at the levels of smaller vessels is also an important marker of cerebrovascular health, though few techniques currently exist to measure it. A technique using near-infrared spectroscopy exists, 11 but few MRI-based techniques. Critically sampled fMRI data enables the estimation of cardiorespiratory pulsation with high accuracy. 36 However, many large publicly available datasets only provide lower-sampled fMRI data. Recently, a technique based on hypersampling of the cardiac cycle and the resting state fMRI signal was developed, 37 as well as a tool to derive the cardiac cycle from the BOLD signal alone in non critically sampled data. 38 This opens the door to measurement of CP in any dataset that includes a resting state fMRI acquisition. Here, we will assess the effects of PA on CP in sex-disaggregated data, to further understand sex-specific patterns in this relationship.
In summary, the primary objective of this study was to evaluate the sex-specific effects of intensity and amount of time spent in PA on CVR and CP in middle-aged and older adults. We hypothesized that PA would be negatively associated with CVR7,28,29 and CP. 34 The secondary objective was to assess the effect of menopause status in the effects of intensity and dose of PA on CVR and CP. We hypothesized that postmenopausal females would exhibit a stronger response to the effects of PA on CVR and CP compared to premenopausal females. 39 The exploration of these objectives will allow us to better understand how PA participation can be optimized to benefit cerebrovascular health in male and female aging.
Material and methods
Participants
The study included 725 adults (406 females & 319 males) older than 36 years old from the Human Connectome Project – Aging (HCP-A) dataset. The objective of HCP-A is to explore ‘normal’ aging, including vascular risk factors. Participants were screened and excluded for any history of major psychiatric or neurological diseases. More details on the screening process of the HCP-A can be found here. 40
From the initial sample, 25 participants (11 females, 14 males) were excluded due to PA outliers with abnormal PA values (example: 29 hours of PA per day), 25 participants (14 females, 11 males) due to missing systolic blood pressure information and 3 participants (3 females) due to missing body mass index information. Individuals older than 85 years old were also excluded due to the variability induced from the superagers profile (25 females, 17 males). 41 Following these exclusions, 626 participants (350 females, 276 males) were included in the overall statistical analysis. Additionally, the HCP-A dataset provided demographic information on menopausal status. 42 In follow-up analyses on the impact of menopause status, 7 female participants were excluded due to missing menopausal status data. Consequently, 126 premenopausal females and 217 postmenopausal females were included in the statistical analysis.
All participants of the HCP-A dataset provided signed informed consent in accordance with the declaration of Helsinki. Secondary use of this publicly-available data was approved by Concordia University’s Human Research Ethics Committee.
Physical activity measurements
The international physical activity questionnaire (IPAQ) was used to quantify the weekly intensity and dose of PA. The IPAQ was designed to provide a valid measurement of self-reported PA that is reproducible across countries. 43 The activities from the IPAQ were categorized into metabolic equivalent of task (METs) from the 2011 Compendium of Physical Activities. The PA intensities were classified as walking, moderate (3.0 to 6.0 METs) and vigorous (>6.0 METs) PA. The total dose of PA was measured as a combination of all PA levels.
The weekly dose of PA in METs were further dissociated into participation levels. Participants were divided into four groups according to the amount of weekly PA participation including: no (0 METs), low (0.01 to 500 METs), medium (500 to 1500 METs) and high (>1500 METs) PA. These ranges of levels of PA were based on work by. 44
MRI acquisition
All MRI data were acquired using a Siemens 3T Prisma machine across four sites. The BOLD sequence was obtained using a 2D multiband (MB) gradient echo echo-planar imaging (EPI), with 2 mm isotropic whole-brain coverage (resolution = 2.0 mm isotropic, repetition time = 0.8 s, echo time = 0.037 s, MP factor = 8, flip angle = 52°, Volumes = 488). The HCP-A dataset included 4 resting state fMRI runs across two sessions, with two runs per session, each 6.5 minutes. For this study, a randomly chosen single run from the first session, utilizing a posterior-to-anterior phase encoding polarity, was selected. Although the HCP-A dataset included concatenated runs of preprocessed fMRI data, concatenated data was not compatible with the CP measurement so we elected to compute relative CVR and CP from the same fMRI session and encoding polarity. The T1-weighted structural data was acquired using a multi-echo magnetization prepared rapid gradient echo imaging sequence (resolution = 0.8 mm isotropic, echo time = 1.8/3.6/5.4/7.2 ms, repetition time = 2.5 s, flip angle = 8°) More details on the HCP-A resting state fMRI acquisition can be found here. 45
Image processing
The HCP-A dataset included minimally preprocessed resting state fMRI data, with key preprocessing steps. Freesurfer 46 and FSL 47 were employed for functional and structural data processing. 48 Gradient distortion correction was applied to the BOLD images to account for B0 field inhomogeneities. Subsequently, the BOLD time series were aligned to a single band reference to correct for motion artifacts. The BOLD data were linearly demeaned and detrended to remove signal drifts. Echo planar imaging distortion correction was performed using a spin echo field map, further refining spatial accuracy. 49 From the preprocessed structural data, the T1-weighted structural image was subsequently registered to the single-band reference. 50 The resting state fMRI data was non-linearly registered to MNI space using a one-step spline resampling, followed by intensity normalization and brain masking using Freesurfer. Bandpass filtering was not applied, as the lower BOLD frequencies were necessary for the relative CVR maps. The HCP-A preprocessing pipeline, described in detail by, 48 was used to ensure consistent and high-quality data processing.
Relative CVR analysis
The relative CVR maps were computed from the preprocessed resting state fMRI data. The first step of the relative CVR processing pipeline (Figure 1) was to perform smoothing of BOLD images (FWHM = 8 mm). A time domain-filtering (0–0.1164 Hz) of the BOLD data was conducted with scipy to extract a proxy of the end-tidal CO2 (ETCO2). This frequency range was chosen based on previous work to isolate the relevant physiological signal associated with ETCO2. 51 This processing procedure precluded the necessity for a measurement of the ETCO2, which is not included in the HCP-A data. Additionally, the filtering of the BOLD signal generates an approximation of resting ETCO2 fluctuations, providing relative rather than absolute values. Thus, this CVR method bypasses the need for temporal alignment of the ETCO2 data, as the proxy is directly estimated from BOLD signal changes. This contrasts with traditional CVR gas manipulation techniques, where ETCO2 is recorded externally at the nose or mouth. The filtered BOLD signal, acting as a proxy of ETCO2, was included as a regressor in the general linear model to measure the CVR indices. The CVR index map was normalised to a whole brain CVR average to obtain the relative CVR map. The resting state CVR processing pipeline was based on previous work by. 31 Relative CVR was averaged over bilateral regions of interests (ROIs) in grey matter of the LBPA40 atlas 52 normalised to MNI-152 space.

Relative CVR pipeline from preprocessed BOLD data.
Cerebral pulsatility analysis
A proxy of cerebral pulsatility was extracted from the unprocessed BOLD signal using the Hypersampling by Analytic Phase Projection - Yay (Happy) tool, 38 a part of the Rapidtide package (doi:10.5281/zenodo.814990). Happy is a Python based tool that can be used to estimate the cardiac waveform from the highest signal-to-noise ratio voxels in unprocessed resting state fMRI data. This signal is then cleaned using a deep learning filter. The resulting cardiac trace is used to estimate voxel-wise pulsatility from the BOLD signal timecourse using Voss and colleagues’ hypersampling technique. 37 This method harnesses the cardiac contamination within the BOLD signal to extract the cardiac cycle and estimate pulsatility. Happy provides a proxy for cerebral pulsatility measured as the BOLD amplitude variations over a cardiac cycle. Because veins are thought to be passively compliant, cardiac variations in the BOLD signal arise as veins passively dilate in response to changes in pressure throughout the cycle. Therefore, while this is not a direct measure of pulsatility of arteries, it is here used as a proxy for local pulsatility.
The unprocessed resting-state fMRI data from the HCP-A was chosen from the same session and phase encoding polarity as the relative CVR analysis to maintain consistency. Cerebral pulsatility (CP) proxy maps were obtained from averaged amplitude BOLD variations over a cardiac cycle using Happy. A brain mask of the fMRI data was extracted using the brain extraction (BET) tool from FSL and was included within the Happy analysis. The CP maps were registered to MNI space using FNIRT from FSL. Additionally, the Jacobian files from the HCP-A were added to the registration of the CP maps to correct from the significant EPI distortions. 48 Lastly, a rigid registration was applied with ANTs to improve the alignment of the CP maps to the MNI space. More information on the methodology of Happy can be found here. 38 Values were averaged over the bilateral regions of the LBPA40 atlas for subsequent statistical analysis.
Statistical analysis
Analyses were performed initially in all subjects with sex as a covariate to investigate sex differences. Then, the analyses were repeated in each sex separately to assess the sex-specific relationships between cerebrovascular properties and PA. This approach is recommended to comprehensively assess the presence of sex-dependent effects.53,54 Analyses controlled for age, systolic blood pressure and body mass index. Linear and quadratic multiple regressions were conducted to assess the relationship between the intensity and dose of PA and CVR using ordinary least squares from the statsmodels python program. A Tukey’s Honest Significant Difference Analysis was performed to assess the amount PA participation within a specific intensity or dose of PA on CVR using the pairwise post hoc test from the python library statsmodels. All statistical analyses were repeated to evaluate the relationship between PA and CP. Multiple comparisons across all analyses were corrected with false discovery rate (FDR). Lastly, a corrected p-value threshold of less than 0.05 was used to determine statistical significance.
Results
The demographic data of all groups are shown in Table 1. In all participants, we observed significant positive associations between relative CVR and age (p < 0.001; R = 0.201) (Supplementary Figure 1), as well as a significant negative correlation between CP and age (p < 0.001; R = 0.332) (Supplementary Figure 2) in the whole gray matter.
Group differences between males and females in demographic, exercise and brain-related variables.
SBP: systolic blood pressure; DBP: diastolic blood pressure; MAP: mean arterial pressure; BMI: body mass index; METs: metabolic equivalent of tasks; CVR: cerebrovascular reactivity; CP: cerebral pulsatility.
Sex-specific effects in self-reported PA and relative CVR
In all participants, a significant negative linear relationship between relative CVR and total PA was observed in the whole gray matter (p = 0.010; R = 0.241) (Supplementary Figure 3). There was no significant effect of sex on the model. The effects of the different doses of PA on relative CVR in both sexes separately are demonstrated in Figure 2. In females, a medium dose of PA (500 to 1500 METs) was associated with higher relative CVR in the cerebral cortex than low (p = 0.041; Cohen’s d = −0.45) or high (p = 0.007; Cohen’s d = 0.40) PA doses. In males, different doses of PA did not lead to significant differences in cerebral cortex CVR.

Effects of total PA dose on CVR. Box plots from Tukey’s tests showing the relationship between relative CVR and total PA dose in males and females. Green) Relative cerebral cortex CVR and total PA dose for males (No: n = 4; Low: n = 28; Medium: n = 81; High: n = 163); Red) Relative cerebral cortex CVR and total PA dose for females (No: n = 9; Low: n = 48; Medium: n = 105; High: n = 188). Significance (FDR corrected p < 0.05) is indicated by *.
Males
The effects of total PA, as well as dose of specific PA intensities on relative CVR was assessed. Males showed significant negative linear relationships between relative CVR and total PA in the bilateral superior parietal gyri (p = 0.013; R = 0.195) and the insula (p = 0.026; R = 0.158), while a significant positive relationship was observed in the cerebellum (p = 0.035; R = 0.195). Additionally, negative linear relationships were identified between relative CVR and walking PA in the bilateral superior parietal gyri (p = 0.043; R = 0.212) and between relative CVR and vigorous PA in the bilateral cingulate gyri (p = 0.046; R = 0.164). These correlations are displayed in Figure 3.

Effects of different PA intensity on CVR in males. Multiple linear regressions displaying the relationship between relative CVR and different doses of PA at different intensities in males. Linear relations between relative CVR and (a) total PA in the cerebellum; (b) total PA in the superior parietal gyrus; (c) total PA in the insula; (d) walking PA in the superior parietal gyrus and (e) vigorous PA in the cingulate gyrus. P-values shown were FDR corrected.
Females
Females showed a single negative linear relationship between relative CVR and walking PA in the bilateral cingulate gyri (p = 0.033; R = 0.141). Further analyses were conducted with females separated according to menopausal status into premenopausal and postmenopausal groups. Premenopausal females revealed significant negative linear relationships between relative CVR and total PA in the superior occipital gyri (p = 0.026; R = 0.319). We also observed relationships between walking PA in the bilateral superior occipital gyri (p = 0.008; R = 0.361), as well as in the middle occipital gyri (p = 0.011; R = 0.286). In postmenopausal females, we only observed a negative linear relationship with walking PA in the bilateral inferior occipital gyri (p = 0.002; R = 0.237). These results are presented in Figure 4.

Effects of different PA intensity on CVR in females. Multiple linear regressions exhibiting the relationship between relative CVR and dose of PA at different intensities in all, premenopausal and postmenopausal females. Linear relations between relative CVR and (a) total PA in the superior occipital gyrus; (b) walking PA in the superior occipital gyrus; (c) walking PA in the middle occipital gyrus in premenopausal females; (d) walking PA in the inferior occipital gyrus in postmenopausal females and (e) walking PA in the cingulate gyrus in all females. P-values shown were FDR corrected.
Sex-specfic effects in self-reported PA and CP
There were no significant interactions between PA, sex, and CP when all participants were pooled. The sex-specific effects of total PA on CP are presented in Figure 5. Females showed significantly lower CP in medium compared to low PA participation (p = 0.016; Cohen’s d = 0.49) and from high to low PA participation (p = 0.027; Cohen’s d = 0.44) in the cerebral cortex. However, males did not demonstrate any significant effects of the overall dose of PA on CP in the cerebral cortex.

Effects of PA dose on cerebral cortex CP. Box plots from Tukey’s tests presenting the relationship between CP and total PA in males and females. Green) CP and total PA in the cerebral cortex of males (No: n = 4; Low: n = 27; Medium: n = 77; High: n = 159); Red) CP and total PA in the cerebral cortex of females (No: n = 9; Low: n = 48; Medium: n = 103; High: n = 177). Significance (FDR corrected p < 0.05) is indicated by *.
No effects of PA dose at any intensity was associated with CP in males or when all females were combined in the same group.
The effects of the total PA participation levels on regional CP in all groups are presented in Figure 6. In males, a significantly lower CP was observed in the bilateral middle temporal gyri from no to low PA participation (p = 0.049; Cohen’s d = 1.28). In females, the lowest CP was consistently observed at medium total PA dose (p < 0.05). This included a significantly lower CP in medium as compared to low PA participation in all females in the superior temporal gyri (p = 0.039; Cohen’s d = 0.43). Postmenopausal females were the group where this effect was the most pronounced, with lower CP in the no to medium PA participation in the bilateral inferior temporal gyri (p = 0.032; Cohen’s d = 1.23), cingulate gyri (p = 0.039; Cohen’s d = 1.37) and the cerebellum (p = 0.027; Cohen’s d = 1.08), as well as from no to low PA participation in that same region (p = 0.028; Cohen’s d = 1.05). Conversely, we observed a higher CP from medium to high PA participation in both the bilateral cingulate gyri (p = 0.044; Cohen’s d = −0.43) and the hippocampi (p = 0.025; Cohen’s d = −0.45).

Effects of total PA on regional CP. Box plots from Tukey’s tests displaying the relationship between CP and dose of PA in males, all females and postmenopausal females. (a) Males CP in the middle temporal gyrus (No: n = 4; Low: n = 27; Medium: n = 77; High: n = 159); (b) Females CP in the superior temporal gyrus (No: n = 9; Low: n = 48; Medium: n = 103; High: n = 177); Postmenopausal females CP in (c) Inferior Temporal Gyrus; (d) Cingulate Gyrus; (e) Cerebellum and (f) Hippocampus (No: n = 6; Low: n = 30; Medium: n = 59; High: n = 113). Significance (FDR corrected p < 0.05) is indicated by *.
Finally, different PA intensities were not associated with differential effects on CP in males or in the combined group of females. However, some associations were identified once females were segregated into postmenopausal and premenopausal groups. These correlations can be found in (Supplementary Figure 4).
To confirm the validity of using SBP to control for hypertension, all analyses were repeated with mean arterial pressure (MAP) Supplementary Tables 1 & 2, resulting in almost identical findings.
Discussion
This study assessed sex-specifc effects of dose of PA at different intensities on proxies of cerebrovascular reactivity and cerebral pulsatility. In males, negative linear relationships between relative CVR and PA were shown in several cortical gray matter ROIs particularly for total PA, while a positive linear relationship was identified in the cerebellum. The effects of PA on CP were limited to a single cortical gray matter ROI in males, showing lower pulsatility from no to a low dose of total PA in the middle temporal gyri. In females, relative CVR and CP both showed a pattern whereby the greatest effects of overall dose of PA across all intensities was for a medium PA dose. This medium dose corresponded to the highest relative CVR and the lowest CP values. However, PA dose at specific intensities showed a more complex picture for the relationship between PA and relative CVR, where higher dose of total or walking PA were associated with lower relative CVR, especially in occipital regions as well as the cingulate. This pattern was present in a greater number of regions pre-menopause. In terms of CP, we observed a consistently lower CP for medium level of participation, as compared to no, low and high participation, and this effect was present over more regions in postmenopausal females than other groups. Together, these findings suggest that males and females’ cerebral vasculature responds differently to PA, with menopausal status having a significant effect.
Sex-specific effects of PA on relative CVR
The results of the CVR-PA relationships indicated predominantly negative linear patterns, with a greater reduction in relative CVR from higher doses of PA. Notably, these findings were both sex-specific and present across several gray matter regions. These regional inverse relationships are consistent with MRI finding showing lower CVR with higher cardiorespiratory fitness or in Master athletes.7,28 The relationship may be more complex however and this inverse relationship was mostly identified in healthier-than average samples. In contrast, Dubose and colleagues 29 discovered a non-linear relationship between CVR and cardiorespiratory fitness whereby individuals with lower cardiorespiratory fitness showed higher CVR with greater fitness, while individuals at the more highly fit end of the spectrum exhibited a negative relationship. This inverse-U nonlinear effect echoes our findings in the cerebral cortex region in females, where we observed the highest relative CVR in females with a medium dose of total PA, as compared both low or high PA dose. Furthermore, this reduced regional relative CVR in individuals with high PA dose is not unique to females since we observed a similar negative effect in males.
A relationship between PA and relative CVR was also observed for different PA intensities. More specifically, males were responsive to all intensities of PA except for moderate PA, whereas all three female groups were primarily sensitive to walking PA. Furthermore, the bilateral cingulate gyri was the only common significant ROI between both sexes. Therefore, the regional variations in the effects of PA on relative CVR between males and females may stem in part from differences in regional perfusion. Consistent with this, sex differences in regional and whole-brain perfusion have been previously reported. 55 Thus, the effects of PA on relative CVR could be influenced by these sex-specifc regional differences in cerebrovascular changes.
In all groups, the relative CVR-PA relationship was predominantly observed in cortical gray matter ROIs. However, females' unique relative CVR-PA relationship in the overall region of the cerebral cortex could indicate a greater whole cortical brain sensitivity to relative CVR. It is possible that this relationship is due to the fact that females have higher CBF in the whole-brain, including gray and white matter and greater perfusion differences between gray and white matter 56 than males.
Our results highlight an important difference in the PA intensities that are associated with cerebrovascular adaptations. While males displayed relative CVR relationships predominantly with total or vigorous PA, all three female groups (all, premenopausal and postmenopausal females) displayed significant negative linear relative CVR-PA relationships with walking PA. This is in line with studies showing that females benefit more from lower intensities of PA compared to males. 57 For instance, walking PA has been shown to greatly reduce the prevalence of cardiovascular diseases in females. 58 Additionally, females appear to experience improved mental health from lighter PA, while males may benefit more from higher intensities of PA. 59 In postmenopausal females, the high PA participation group showed higher CP compared to the medium group in the cingulate gyrus and hippocampus. While negative associations are typically documented, Reed and colleagues (2024) 60 observed a positive association following 12-weeks of training in middle-aged adults with elevated blood pressure or stage-1 hypertension, and speculated that this could arise from an increased stroke volume prior to vascular adaptation that dampen pulsatile flow. However, mechanistic studies are needed to confirm and explain these effects.
Sex-specific effects of PA on CP
The findings of this study revealed an overall lower CP with higher PA participation across all groups. Since cerebral pulsatility is increased with arterial stiffness, this direction of change is consistent with a beneficial effect of PA on vascular health. Because we did not find any relationship between specific PA intensities and CP, the benefits of PA on CP appear to be primarily influenced by total PA participation rather than a specific PA intensity. These results align with those reported by Mohammadi and colleagues, where a 12-month PA intervention resulted in a significant reduction in CP measured using near-infrared spectroscopy in older adults. 11 The improvements in CP from PA are likely due to reduced arterial stiffness which decreases the amplitude of the cardiac pressure wave in the brain. 61 Notably, lower CP measured by 4D flow MRI, is strongly correlated with lower arterial stiffness in many large cerebral arteries. 62 Moreover, greater cardiorespiratory fitness has been inversely associated with 4D flow MRI CP in the internal carotid artery and basilar artery. 63 Therefore, total PA, regardless of any specific intensity, likely reduces CP by improving arterial stiffness.
Males showed limited benefits of PA on CP, in contrast to the stronger relationship observed with relative CVR in this group. Therefore, these findings suggest that in males, cerebrovascular health may be more strongly influenced by factors related to vascular reactivity than by protective effects of PA on resting state fMRI-based CP measurement. Consistent with this, the effects of PA on CVR are well documented in males,7,28,29 while its impact on CP is mainly observed in females. 64 These sex-specific effects are likely driven by the earlier onset of cerebral artery stiffness in females, resulting from the rapid decline in estrogen levels after menopause, which amplifies cranial pulse waves. 65 Consequently, the effects of PA on CP may be smaller in other groups as compared to postmenopausal females.
Our results in all and postmenopausal females are consistent with previous studies that observed a negative relationship between CP and cardiorespiratory fitness only in females.66,67 Reduced estrogen after menopause is associated with decreased nitric oxide production and impaired endothelial function, 39 leading to greater arterial stiffness. This, in turn, reduces the arteries’ ability to absorb the pulsatility, resulting in increased pulsatility reaching to the brain. 68 Furthermore, older females have been identified as being at higher risk for elevated CP, 69 highlighting the importance of PA at midlife during the transition to menopause and after. 70
Potential mechanisms
A few potential mechanisms underlying the negative CVR-PA relationship have been previously proposed, including changes in chemosensitivity,7,28 cerebral autoregulation and neuroimaging modality differences.7,29 However, Dubose and colleagues 29 discussed the lack of evidence supporting an attenuation of chemosensitivity due to PA. Additionally, several studies did not report delayed cerebral autoregulation from PA.71 –73 Conversely, in individuals with mild traumatic brain injury, PA has been recognised as a potential intervention to improve cerebral autoregulation. 74 In contrast, differences in methodology likely explain some contradictory results in the literature. Especially, MRI and TCD measure very different aspects of the vasculature, with TCD being sensitive to changes in blood velocity in large arteries, and MRI measurements being predominantly sensitive to smaller vessels. In fact, MRI and TCD CVR measures have been shown to lack correlation. 75
One aspect that may contribute to this negative relationship between PA and CVR is arteriole structural remodelling. The concept of the athlete’s artery was first introduced by Green and colleagues (2012). Endurance trained individuals were shown to have increased lumen size and reduced vessel wall thickness. Evidence from the literature suggests that initially, in untrained individuals, PA may elicit beneficial adaptation to the function of the blood vessel through endothelial improvements, increasing nitric oxide expression 76 and decreasing oxidative stress. 77 After a few months of PA participation, repetitive bouts of shear stress are believed to lead to permanent vascular structural remodelling. 76 This structural remodelling allows increases in blood flow. Notably, a 16 week aerobic training intervention in young adults has been shown to enhance CBF in association with greater cerebral artery lumen size. 78 Therefore, vascular remodelling could result in lower CVR by enhancing the capacity to wash out CO2 from a greater lumen size. Conversely, Folklow and colleagues discovered that a greater wall-to-lumen ratio induces exaggerated responses to vasodilatory stimuli. 79 Thus, low levels of PA may cause an excessive vasodilation response, leading to abnormally high CVR, whereas trained individuals are likely to exhibit a normal CVR response. Our findings revealed, in agreement with Dubose and colleagues, 29 that the CVR-PA relationship may not be solely negative, as low to medium levels of PA participation could enhance CVR, at least in the cerebral cortex region in females. Thus, PA may initially improve the function of the cerebral arterioles, increasing CVR, while later structural changes could potentially decrease CVR.
Myogenic tone may also contribute to this negative relationship. Myogenic tone controls how blood vessels contract and dilate in response to changes in blood pressure. 80 Pre-clinical studies have shown that exercise training induces greater myogenic reactivity, 76 meaning more vasoconstriction with increased pressure. In humans, post-exercise recovery has been associated with increased myogenic activity, which improves microcirculation perfusion and facilitates the supply of nutrients to support tissue. 81 Because CO2 increases cardiac output and therefore pressure, higher myogenic tone may lead to a more attenuated response to CO2 and a lower CVR. Therefore, PA could potentially enhance the myogenic tone of arterioles, minimising CVR.
Finally, because the relationship between ETCO2 and CVR is non-linear, higher resting ETCO2 has been linked to reduced CVR.82,83 Individuals with lower respiratory rates have been shown to have elevated ETCO2. Moreover, individuals with a higher baseline ETCO2 levels have an enhanced cardiorespiratory fitness. 84 Higher body weight increases CO2 production which further contributes to greater resting ETCO2 levels. 85 Additionally, males exhibit higher resting ETCO2 levels compared to females regardless of breathing pattern, which may partly explain sex-specifc effects in the CVR-PA relationship.86,87 Thus, lower CVR associated with PA may be influenced by resting ETCO2 levels. Overall, these potential mechanisms offer promising insights into beneficial vascular physiological adaptations.
Limitations
There are a few limitations to this study as the methods to the cerebrovascular measurements rely on resting state fMRI fluctuations. The measurement of relative CVR relies on natural ETCO2 fluctuations. For instance, a steady breathing rate diminishes the sensitivity to CVR since the ETCO2 changes are small. Therefore, hypoperfused regions may suffer from lower CVR values from a lack of detectible vasodilation. Conversely, individuals with erratic breathing due factors like stress could lead to artificially elevated CVR. Future work should validate these exploratory analyses by using quantitative CVR measurements with ETCO2 recordings. The estimation of CP through BOLD originates from venous blood which is not ideal to capture changes in arteries and arterioles. Furthermore, the sensitivity of this CP method is highly reliant on the accuracy of the reconstruction of the heart pulse waveform. Despite these limitations, these techniques are highly valuable, since fMRI BOLD acquisitions are simple to measure and common in both research and clinical routines, making these measurements accessible in already acquired datasets. The fact that these resting state measurements are sensitive not only to the effects of PA, but also to sex specific effects of PA demonstrates their promise as biomarkers of cerebrovascular health.
A second limitation is that relative CVR and CP were measured using non-critically sampled fMRI data (TR = 800 ms). While the techniques used here provide estimates of vasoreactivity and pulsatility, critically sampled fMRI data would provide more accurate measurements of relative CVR and CP. An example of a subject’s estimated cardiac cycle derived from the Happy technique alongside plethysmogram data, (Supplementary Figure 5) as well as the estimated CO2 trace from our CVR pipeline compared to respiratory belt measurements (Supplementary Figure 6) are provided to show data quality.
A third limitation is that PA was quantified from a self-reported questionnaire. Although the IPAQ has been validated across many countries, the results must be interpreted with caution as they are self-reported PA levels. Furthermore, the IPAQ did not assess the participants’ lifetime exposure to PA. Thus, these factors inherently introduced variability in the measurement of PA intensities and doses.
Lastly, our findings demonstrated sex- and menopause status-specific relationships between PA, relative CVR, and CP. However, the mechanisms underlying these relationships remain to be elucidated and future mechanistic studies are needed to understand the biological underpinnings of these observations.
Conclusions
In summary, this study showed that males and females may benefit differently from varying intensities and doses of PA in relation to cerebral hemodynamics. Females, particularly postmenopausal females, presented more significant effects of total PA on CP compared to males. Because the mechanisms that lead to CVR, CP and their measurement are complex, some of the sex-specifc effects of PA on CVR and CP may be due to different dominant mechanisms of vascular decline and remodelling in males and females. The observed negative relationship between relative CVR and PA may be influenced by physiological mechanisms such as vascular structural remodeling, myogenic reactivity and resting ETCO2 levels. Moving forward, research employing gold standard techniques in large, well-characterized cohorts is essential to unravel the complex physiological processes involved and to optimize PA interventions for promoting cerebral vascular health.
Supplemental Material
sj-pdf-1-jcb-10.1177_0271678X251325399 - Supplemental material for Sex-specific effects of intensity and dose of physical activity on BOLD-fMRI cerebrovascular reactivity and cerebral pulsatility
Supplemental material, sj-pdf-1-jcb-10.1177_0271678X251325399 for Sex-specific effects of intensity and dose of physical activity on BOLD-fMRI cerebrovascular reactivity and cerebral pulsatility by Zacharie Potvin-Jutras, Brittany Intzandt, Hanieh Mohammadi, Peiying Liu, Jean J Chen and Claudine J Gauthier in Journal of Cerebral Blood Flow & Metabolism
Footnotes
Ethics approval and informed consent statements
All participants of the HCP-A dataset provided signed informed consent in accordance with the declaration of Helsinki. Secondary use of this publicly-available data was approved by Concordia University’s Human Research Ethics Committee.
Data availability statement
The Human Connectome Project - Aging dataset supporting this study is openly available at (https://www.humanconnectome.org/study/hcp-lifespan-aging). The method to the processing of the cerebral pulsatility is openly available at (
).
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: Research reported in this publication was supported by the National Institute On Aging of the National Institutes of Health under Award Number U01AG052564 and by funds provided by the McDonnell Center for Systems Neuroscience at Washington University in St. Louis. Funding for this project is provided by the Heart and Stroke Foundation of Canada and Brain Canada (to Zacharie Potvin-Jutras). Hanieh Mohammadi is supported by a grant from Canadian Heart and Stroke Foundation. This study was supported by the Canadian Natural Sciences and Engineering Research Council (RGPIN-2024-06455, to Claudine J. Gauthier), the Michal and Renata Hornstein Chair in Cardiovascular Imaging (to Claudine J. Gauthier) and Canadian Institutes of Health Research (468740, to Claudine J. Gauthier)
Acknowledgements
Declaration of conflicting interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Authors’ contributions
Zacharie Potvin-Jutras contributed to the study conception and design; data curation; methodology; data analysis; writing of the manuscript; reviewing and editing of manuscript. Brittany Intzandt contributed to data analysis; reviewing and editing of the manuscript. Hanieh Mohammadi contributed to reviewing and editing of the manuscript. Peiying Liu contributed to methodology; reviewing and editing of the manuscript. Jean J. Chen contributed to reviewing and editing of the manuscript. Claudine J. Gauthier contributed to study conception and design; methodology; reviewing and editing of manuscript; supervision.
Supplementary material
Supplemental material for this article is available online.
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
