Abstract
White matter lesions (WML) have been linked to cognitive decline in aging as well as in Alzheimer’s disease. While hypoperfusion is frequently considered a cause of WMLs due to the resulting reduction in oxygen availability to brain tissue, such reductions could also be caused by impaired oxygen exchange. Here, we tested the hypothesis that venous hyperintense signal (VHS) in arterial spin labeling (ASL) magnetic resonance imaging (MRI) may represent a marker of impaired oxygen extraction in aging older adults. In participants aged 60–80 years (n = 30), we measured cerebral blood flow and VHS with arterial spin labeling, maximum oxygen extraction fraction (OEFmax) with dynamic susceptibility contrast, and WML volume with T1-weighted MRI. We found a significant interaction between OEFmax and VHS presence on WML volume (p = 0.02), where lower OEFmax was associated with higher WML volume in participants with VHS, and higher OEFmax was associated with higher WML volume in participants without VHS. These results indicate that VHS in perfusion-weighted ASL data may represent a distinct cerebrovascular aging pattern involving oxygen extraction inefficiency as well as hypoperfusion.
Introduction
White matter lesions (WMLs) are a common manifestation of cerebral small vessel disease, 1 and the presence of these lesions has been linked to declining cognitive function in aging as well as in patients with Alzheimer’s disease and related dementias.2–8 While the etiology of these lesions is often presumed to be vascular in nature,9,10 microvascular hemodynamic mechanisms associated with the formation and progression of these lesions are highly complex and remain incompletely understood. Thus, there is an urgent need to improve our understanding of these mechanisms so that therapies aimed at preventing or delaying WML-related cognitive decline may be developed and titrated to high-risk individuals.
Noninvasive magnetic resonance imaging (MRI) techniques including arterial spin labeling (ASL) 11 have frequently been used to study how cerebral blood flow (CBF) changes with aging. 12–17 In pseudo-continuous ASL, 18 the magnetization of arterial blood water in the neck is inverted (i.e., labeled) using a long chain of radiofrequency pulses. After a post-labeling delay to allow this labeled blood to reach and perfuse tissue, MR images of a volume containing the brain are acquired. These ‘label’ images are subtracted from ‘control’ images acquired without labeling to yield perfusion-weighted images that are converted to CBF using a kinetic model. 19
While ASL has proven difficult for measuring perfusion in white matter, gray matter CBF has been extensively studied as a biomarker in various diseases. Using ASL or other MRI approaches, overwhelming evidence has been presented that cerebral perfusion decreases with increasing age and represents a marker of aging-related hemodynamic impairment.13,20,21 Therefore, WML formation and progression has been frequently attributed to aging-related hypoperfusion or impairment in the regulation of CBF.
22
However, perfusion represents only one component of the equation that governs the relationship between cerebral oxygen supply and demand:
Recently, venous hyperintense signal (VHS) in perfusion-weighted ASL images has been characterized as a marker of capillary shunting and reduced OEF in young adults with sickle cell anemia.25,26 In the absence of pathology affecting microvascular flow, it is expected that labeled arterial blood water arrives in the capillary bed and exchanges sufficiently with tissue water, resulting in outflow into veins that contains a good mixture of labeled and unlabeled spins. Therefore, little to no ASL difference signal is expected to be present in large veins. However, in the presence of disturbed microvascular flow, it is possible that labeled blood water passes through the capillaries and into the venous system too quickly to allow for good mixing with tissue water. In this case, the presence of labeled arterial blood causes bright signal to appear inside of large draining veins (Figure 1). Since the efficiency of oxygen exchange in the capillaries is dependent in part on these microvascular transit times, we posited that the presence of VHS could be indicative of oxygen exchange efficiency.

Venous hyperintense signal: cerebral blood flow (CBF)-weighted images derived with arterial spin labeling (ASL) were assessed for venous hyperintense signal (VHS). Representative axial, sagittal, and coronal CBF images (from left to right) are shown here for individuals with VHS present (a; VHS+; 79 y.o. female) and VHS absent (b; VHS−; 64 y.o. female). Participants were classified as VHS+ if bright signals were contiguously observed at locations where large draining veins, including the superior sagittal sinus and straight sinus, are expected to be present (green arrows). Images where no signal enhancement was observed relative to surroundings in such regions (green circle) were classified as VHS−.
Prior studies in young anemic adults have shown that the presence of VHS was indicative of reduced OEF and could be associated with elevated risk for silent lesions. 26 Interestingly, such relationships were also observed, but at a lower prevalence, in young adults who were not anemic but exhibited radiological evidence of WMLs. 26 In the present study, we extended this work to investigate whether VHS may also be present in older adults and whether it may represent a more general marker of hemodynamic impairment associated with oxygen extraction efficiency in an aging population. Based on prior findings in young adults, we hypothesized that older adults with VHS would exhibit impaired OEFmax and elevated burden of WMLs compared to those without VHS.
Methods
Study design and participants
All components of this study were performed in compliance with the ethical standards instituted by the Mass General Brigham Research Compliance office and as set forth by the Belmont Report. Older adult volunteers (60–80 years of age) were recruited through the Brain Aging and Dementia Laboratory at Massachusetts General Hospital and provided informed written consent as part of this Mass General Brigham Institutional Review Board-approved study (IRB #2013P001639). Exclusion criteria included major neurological or psychiatric disorders, including dementia or clinical stroke, as well as the presence of substantial systemic illnesses likely to confound the study. Participants underwent two imaging sessions: a non-contrast MRI scan at 3 Tesla (3T; Siemens Biograph mMR) and a contrast-enhanced MRI scan at 7 Tesla (7T; Siemens Magnetom). A battery of neuropsychological testing, including the Montreal Cognitive Assessment (MoCA), was administered to measure overall cognitive health as previously described, 27 and participants were classified as exhibiting mild cognitive impairment (MCI) based on performance falling one standard deviation or more below published normative values on at least two cognitive domains. 28 Of note, only data from participants who underwent both imaging sessions were included in this retrospective study.
Image acquisition
T1-weighted MR images were acquired at 3T using a multi-echo magnetization-prepared rapid gradient echo (ME-MPRAGE) sequence (repetition time [TR] = 2530 ms; inversion time [TI] = 1100 ms; echo times [TE] = 1.69 ms, 3.55 ms, 5.41 ms, and 7.27 ms; flip angle = 7°; spatial resolution = 1.0 × 1.0 × 1.0 mm3). The individual echoes were combined by using a root mean square to generate a single T1-weighted volume for tissue segmentation and image registration. Fluid-attenuated-inversion-recovery (FLAIR) T2-weighted MR images were also acquired at 3T (TR = 5000 ms; TE = 400 ms; flip angle = 120°; spatial resolution = 0.5 × 0.5 × 1.0 mm3).
ASL MR images were acquired at 3T with pseudo-continuous labeling (total labeling duration = 1500 ms and post-labeling delay = 1200 ms) and without background or vascular suppression using a 2D echo planar imaging readout (TR = 4000 ms; TE = 12 ms; spatial resolution = 3.5 × 3.5 × 5.0 mm3; control/label pairs = 40).
Dynamic susceptibility contrast (DSC) MRI data were acquired at 7T with gadolinium contrast. A bolus of 0.2 ml/kg (0.1 mmol/kg) of gadopentetate dimeglumine (MAGNEVIST, Bayer Healthcare) was administered intravenously at a constant rate between 1.0–2.0 ml/s, followed by a saline flush at the same rate. A 2D single-shot gradient-echo echo planar imaging sequence was used for the readout (TR = 1500 ms; TE = 22 ms; flip angle = 75°; spatial resolution = 1.5 ×1.5 × 1.5 mm3; nominal echo spacing = 0.69 ms; Generalized Auto-calibrating Partially Parallel Acquisition factor = 3; multi-band factor = 2). A total of 100 frames were collected.
Structural image processing
T1-weighted images were automatically processed to reconstruct cortical surfaces and to segment volume regions-of-interest using the Freesurfer ‘recon-all’ procedure.29,30 Segmentation of these images using a watershed/surface deformation procedure and automatically registered into a common surface template using a surface-based averaging technique by considering cortical folding patterns. 31 A pre-defined gray matter mask in the standard space was then inversely mapped to each participant’s native space using a high-dimensional spherical morphing procedure to define gray matter masks in each participant’s images. WML volumes were calculated using T1-weighted images with the FreeSurfer-based procedure for measuring white matter hypointensities as described previously. 2
Physiological image processing
ASL images from each participant were motion-corrected
32
and pair-wise subtracted.
33
CBF was then computed using a two-compartment kinetic model as previously described:
34
DSC-MRI data were motion-corrected 32 and processed using the Penguin/pgui software to convert the T2* signal time curves of each voxel to contrast concentration time curves. 35 An arterial input function (AIF) from middle cerebral artery branches was selected by the software based on cluster analysis of voxel-wise dynamic curves. 36 The AIF was then used in a vascular model developed by Ostergaard et al.24,35,37 to derive perfusion parametric maps, including the OEFmax images. DSC-MRI data were aligned within-subject to corresponding T1-weighted images using a boundary-based registration algorithm, 31 and this transform was applied to corresponding OEFmax images to enable quantitation in a similar gray matter region-of-interest, as with CBF. DSC data were also used to derive other parameters, including mean transit time, CBF, and cerebral blood volume. Analyses with these parameters as a function of VHS presence are presented as supplementary information.
Statistics
Continuous variables are described with the corresponding median and interquartile range, and prevalence measures are described with percentages. All statistical analysis was performed using two-tailed tests when applicable and at a 0.05 level of significance. Continuous variables were assessed for normality using Shapiro-Wilk tests.
The primary statistical objective of this study was to test the hypothesis that older age is associated with lower OEFmax and elevated WML volume in VHS+ versus VHS− participants. First, we assessed associations between CBF and age and between OEFmax and age using Spearman’s correlation analysis. To study the effect of VHS presence, we then performed linear regressions with the inclusion of a term modeling the interaction between VHS and age on each of CBF and OEFmax separately. Second, we assessed whether WML volume was different between VHS+ and VHS− participants using a two-tailed Wilcoxon rank-sum test. We then assessed between WML volume and age using Spearman’s correlation analysis followed by a linear regression modeling the interaction between VHS and age on WML volume. Finally, we assessed the relationships between CBF and WML volume and between OEFmax and WML volume using Spearman’s correlation analysis. We also performed linear regressions modeling the interaction between VHS and each of CBF and OEFmax on WML volume separately.
Results
Demographics and prevalence of venous hyperintense signal
A total of 30 individuals satisfied the inclusion criteria for this study. Demographic information, cognitive data, and vascular risk factor prevalence from these participants are shown in Table 1. Overall, ten participants (33%) exhibited neuropsychologically-defined criteria for MCI.
Participant demographics, cognitive scores, and vascular risk factor prevalence.
MoCA: Montreal Cognitive Assessment; IQR: interquartile range; VHS: venous hyperintense signal. p-values are for comparisons between groups with versus without VHS.
Overall, VHS was present in 30.0% of participants (9/30). Participants with VHS (median age = 70 years) were older than those without VHS (median age = 65 years), but this difference was not statistically significant (p = 0.18). There was no significant difference in the male-to-female ratio in the groups with VHS (33% male) and without VHS (52% male; p = 0.34).
Of the continuous variables presented below, age (p = 0.02) and WML volume (p < 0.01) were normally distributed while CBF (p = 0.22) and OEFmax (p = 0.67) were not.
Cerebrovascular physiology
Age was not correlated with ASL-derived gray matter CBF (

Age and hemodynamic physiology: There was no significant interaction between age and venous hyperintense signal (VHS) presence (VHS+) or absence (VHS−) on the association between cerebral blood flow (CBF) and age (a) and There was a significant interaction between age and VHS on maximum oxygen extraction fraction (OEFmax), with older age associated with higher OEFmax in VHS− participants and lower OEFmax in VHS+ participants (b).
There was no relationship between OEFmax and age overall (
White matter lesion volume
Overall, participants with VHS (median volume = 4146 mm3) did not exhibit a significantly greater WML burden compared to those without VHS (median volume = 2548 mm3; p = 0.18; Figure 3(a)). However, older age is significantly correlated with higher WML volume (

Age and white matter lesion volume: Overall, participants with VHS (VHS+) did not exhibit a significantly higher white matter lesion (WML) volume compared to participants without VHS (VHS−) (a) and However, there was a significant interaction between age and VHS on WML volume, with older age associated with a higher WML volume in VHS+ compared to VHS− participants (b).
WML volume was not significantly correlated with CBF (

White matter lesion volume and hemodynamic physiology. There was no significant interaction between cerebral blood flow (CBF) and venous hyperintense signal (VHS) presence (VHS+) or absence (VHS−) on white matter lesion (WML) volume (a). However, there was a significant interaction between maximum oxygen extraction fraction (OEFmax) and VHS on WML volume, with lower OEFmax associated with higher WML volume in VHS+ participants and lower WML volume in VHS− participants (b).
Outputs from regression analyses reported above are summarized in Table 2. Representative examples of participants with and without VHS (Figure 5) demonstrate that these group-level findings can also be observed at the individual level. Analyses conducted with other DSC-based measures, including mean transit time and cerebral blood volume, are shown in Supplementary Figure 2.
Coefficients and p-values of terms from the multivariate regression analyses.
*Represents statistically significant regression parameters. CBF: cerebral blood flow; OEFmax: maximum oxygen extraction fraction; VHS: venous hyperintense signal; WML: white matter lesions.

Representative examples: Venous hyperintense signal (VHS) assessment, cerebral blood flow (CBF), maximum oxygen extraction fraction (OEFmax), T1-weighted, and T2-weighted fluid-attenuated inversion recovery (FLAIR) MR images are shown for two representative participants matched for age and sex but differing by the presence of VHS. CBF exhibited a small difference between the participant with (a) versus without VHS (b), but OEFmax was markedly lower in the participant with VHS, and white matter lesion (WML) volume was much higher in the participant with VHS. As a note, T1-weighted images were used for WML volume quantification, while T2-FLAIR images are shown here for visualization of the lesions.
Discussion
In this study, we characterized a novel imaging marker of oxygen extraction inefficiency and its relationship with age, cerebral hemodynamic physiology, and white matter lesion burden in a cohort of older adults between 60–80 years of age. We used ASL to assess for CBF and for venous signal indicative of capillary shunting and DSC to compute OEFmax as a measure of oxygen exchange efficiency. Our main findings were that VHS was prevalent in 30% of our sample of older adults, and participants with VHS were older (median age = 70 years) compared to participants without VHS (median age = 65 years). While not significantly related to CBF, the relationship between age and OEFmax was dependent on VHS presence, with lower OEFmax corresponding to older age when VHS was present. Finally, higher WML volumes were correlated with lower OEFmax in participants with VHS but not in those without VHS. These results indicate that the presence of VHS on perfusion-weighted ASL data could serve as a marker of a distinct cerebrovascular aging pattern that potentially involves physiological mechanisms involving oxygen extraction inefficiency.
The relationship between oxygen demand and oxygen supply in the brain can be described using Fick’s first law, where the cerebral metabolic rate of oxygen consumption is supported by the product of CBF and OEF. 38 Prior studies that have included a wide age range of participants have shown decreases in CBF with aging,39–41 but relatively few have included participants older than 60 years. Those that include individuals between 60–80 years suggest a weaker association with age and higher variability during these two decades.40,42–44 Additionally, lower CBF was shown to be associated with greater severity of periventricular hyperintensities, which is also suggested by our comparisons in individuals with and without VHS. However, the relationship between OEF and severity was not as clear, with the highest severity group displaying lower OEF than the moderate severity group.45 From our understanding of ischemic disease, it is thought that OEF should increase as CBF decreases to maintain CMRO2.46,47 In the presence of disturbed capillary flow patterns, however, it may no longer be possible for this relationship to persist due to resulting impairment in oxygen extraction efficiency. 48 It has been previously shown that OEFmax is limited by capillary flow patterns, 23 and faster transit across the capillary bed with aging or disease could result in a decrease in OEFmax. While such a gradual decrease may first result in compensatory increases in CBF, sufficient oxygen extraction may only be maintained after a certain threshold only if CBF is suppressed to prolong capillary transit times. 49 This mismatch may in turn result in a chronic imbalance between oxygen demand and supply and elevated levels of hypoxic insult presenting as higher WML volumes. However, the degree to which this phenomenon is present and associated with white matter lesions is not yet fully understood. In this study, we implemented an approach that uses DSC MRI to calculate OEFmax. While these methods were developed for use at 3T, they assume a linear relationship between contrast agent concentration and changes in tissue and blood water relaxation rates and should be equally applicable at 7T.
Our experimental results closely match these theoretical expectations regarding the relationship between aging, OEFmax, and WML volumes. The lack of an expected relationship between CBF and age could be related to the narrow age range studied coupled with relatively small sample size. Meanwhile, we found that OEFmax decreases with increasing age in participants with VHS but increases with older age in participants without VHS. Prior studies in healthy adults have described a steady decrease in venous oxygenation with advancing age, implying that OEF increases with age, across the adult lifespan. 41 However, our results indicate that there is heterogeneity in the directionality of the relationship between age and OEFmax, which may be described by the presence or absence of VHS physiology. Our analyses examining the relationship between hemodynamics and WML volumes provided a very similar result. Higher WML volumes were associated with lower OEFmax in VHS+ participants and with higher OEFmax in VHS− participants. Taken together, these could indicate that VHS may be indicative of impairment in oxygen exchange, which could lead to greater WML burden with aging.
Venous signals in perfusion-weighted MRI acquired with ASL were first characterized as markers of oxygen extraction efficiency in young adults (aged 16−39 years) with sickle cell anemia. 25 In this patient population, OEF was shown to be lower when VHS was present in patients with a history of silent infarcts or white matter lesions. 26 The rationale for capillary flow disruptions leading to tissue hypoxia in this disease where high flow has been demonstrated has been previously described. 50 While the labeling durations (1000 ms vs. 1500 ms) and post-labeling delays (1900 ms vs. 1200 ms) of the ASL sequences used in the two studies are different, the time since the start of labeling (i.e., the sum of the labeling duration and the post-labeling delay) is similar (2900 ms vs. 2700 ms). For VHS-related considerations, this is the quantity of importance since the difference between start of labeling and the time of imaging places a limit on how short the sum of the arterial transit time and capillary transit time must be for labeled blood water to be detected in the veins. We have recently shown that arterial transit times in adults between 60–80 years of age are ∼1.5 s, 14 and typical capillary transit times are thought to be ∼2 s. Therefore, in the absence of pathology, we expect that the shortest imaging time at which labeled signal should be detectable is at least ∼3.5 s. This estimate does not include the time necessary for the labeled blood to travel the venous system to the draining veins, but including this time would only increase this lower limit. However, venous signals in this study were observed at just 2.7 s following the start of labeling, indicating that either the arterial transit time and/or the capillary transit time must be reduced in participants with VHS. Our prior work showing strong increases in arterial transit times with increasing age 14 leads us to therefore posit that VHS is likely due to reduced capillary transit times. Supplementary analyses indeed indicate that capillary mean transit time may also decrease with age in VHS+ individuals (Supplementary Figure 2), but further studies are necessary to confirm this hypothesis. Importantly, only venous signal apparent when using asymmetric labeling schemes (versus symmetric or non-selective labeling) may be interpreted as microvascular transit related since the spin labeling occurs in cervical arteries of the neck and not in the brain. In addition, this phenomenon is distinct from arteriovenous malformations, which result in physical shunts and may also play a role in white matter lesion pathophysiology. 51
Limitations
Our results should be considered in the context of the following limitations. First, our sample size (n = 30 participants) is modest, limiting our ability to generalize the findings from the data. Second, while the timing of our ASL acquisition may be particularly well-suited for studying capillary flow disturbances, the post-labeling delay (=1200 ms) is shorter than what has been recently recommended for using in older adults 19 due to the onset of this study in 2014 prior to wide dissemination of these recommendations. The use of a shorter delay time could introduce arterial transit time as a potential confounding factor into our CBF measures as transit times are known to increase with age. However, this is somewhat mitigated by our use of a 2D echo planar readout and an ascending slice progression for the ASL acquisition, since the effective post-labeling increases for each slice and reaches ∼1500 ms for the most superior slices. While DSC can also provide a measure of CBF, this measure was not used as the primary CBF parameter in this study to avoid the potential for correlated effects with OEFmax. However, as illustrated in Supplementary Figure 3, CBF in WM also indicates differences in hemodynamic relationships based on VHS presence. Third, OEFmax represents only the theoretical upper bound and does not provide a measure of the true OEF of brain tissue. While there are elegant MRI-based approaches available for measuring OEFmax, this parameter was not of primary interest in the study for which these data were acquired. However, ongoing work aims to utilize such approaches to address this issue. Fourth, while our structural marker of interest (WML volume) is a quantity belonging to white matter, our physiological measurements were made in gray matter. ASL in white matter is inherently more difficult due to much lower blood volumes and much longer transit times compared to gray matter. While we have recently shown that these limitations may be overcome with state-of-the-art sequences and updated processing schemes, 14 ASL data used in this study were acquired prior to the implementation of these more recently developed innovations. Fifth, our measure of VHS is currently a qualitative one, but ongoing work involves the development of a multiple post-labeling delay approach for automatically detecting and quantifying venous ASL signal. Such an approach would also allow for estimation of arterial transit times as well as a more robust computation of CBF. Finally, while large vessel disease has been shown to influence microvascular hemodynamics, magnetic resonance angiography data were not acquired as part of this study. Therefore, it was not possible to determine the degree to which these large vessel effects could be present. However, none of the participants included in this study exhibited a history of stroke.
Conclusions
In conclusion, we have characterized the presence of labeled arterial blood water in large draining veins (i.e., capillary shunting) as a potential noninvasive mechanism contributing to reduced oxygen availability to the brain with aging. Using independent measures of VHS (with ASL), OEFmax (with DSC), and WML volume (with T1-weighted imaging), we have shown that the presence of venous signals in perfusion-weighted ASL MRI may be indicative of oxygen extraction inefficiency-related impairment with aging. However, as the physiology underlying these lesions is highly complex, future work is required to understand how VHS and capillary hemodynamics may interact with different risk factors for WMLs and AD to influence progression of these lesions, how they may be associated with measures of cognitive function, and potential regional differences in physiology based on vessel-specific presence of VHS.
Supplemental Material
sj-pdf-1-jcb-10.1177_0271678X221105986 - Supplemental material for Oxygen extraction efficiency and white matter lesion burden in older adults exhibiting radiological evidence of capillary shunting
Supplemental material, sj-pdf-1-jcb-10.1177_0271678X221105986 for Oxygen extraction efficiency and white matter lesion burden in older adults exhibiting radiological evidence of capillary shunting by Meher R Juttukonda, Kimberly A Stephens, Yi-Fen Yen, Casey M Howard, Jonathan R Polimeni, Bruce R Rosen and David H Salat in Journal of Cerebral Blood Flow & Metabolism
Footnotes
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was supported by the National Institutes of Health (R01NR010827) and the American Heart Association (19CDA34790002).
Acknowledgements
We would like to thank Dr. Leif Ostergaard for helpful discussions regarding the applicability of 7T dynamic susceptibility contrast approaches for quantifying hemodynamic parameters relevant to this study.
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
Meher R. Juttukonda: Involved in study conception and design; analysis and interpretation of data; manuscript drafting and revision. Kimberly A. Stephens: Involved in data acquisition; analysis and interpretation of data; manuscript revision; Yi-Fen Yen: Involved in analysis and interpretation of data; manuscript revision. Casey M. Howard: Involved in data acquisition; manuscript revision. Jonathan R. Polimeni: Involved in data acquisition; manuscript revision; Bruce R. Rosen: Involved in analysis and interpretation of data; manuscript revision; David H. Salat: Involved in study conception and design; analysis and interpretation of data; manuscript revision.
Supplemental 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.
