Sage Journals: Discover world-class research

Abstract

We investigate whether common circle of Willis (CoW) variants relate to cerebral blood flow (CBF) characteristics among aging adults. Vanderbilt Memory and Aging Project participants free of clinical stroke (n = 390, 71 ± 8 years, 55% male) completed magnetic resonance angiography to assess cerebral artery variants and pseudocontinuous arterial spin labeling to quantify resting CBF and spatial coefficient of variation (sCoV), a marker of perfusion heterogeneity or stability. Linear regression models related categorical variants to regional CBF and sCoV, adjusting for age, sex, race/ethnicity, education, cognitive status, Framingham Stroke Risk Profile, and APOE-ε4 status. Absence of at least one posterior communicating artery (PcoA) related to lower occipital lobe CBF (p = 0.001) and higher global (p = 0.003), temporal (p = 0.001), and occipital lobe sCoV (p = 0.01). Bilateral PcoA absence related to lower global (p = 0.01), parietal (p = 0.01), and occipital lobe CBF (p = 0.0001) and higher global (p < 0.0001), temporal (p < 0.0001), and occipital lobe sCoV (p = 0.0003). Hypoplastic and fetal-type variants did not relate to outcomes (p > 0.22). Findings demonstrate that missing PcoAs influence both perfusion magnitude and efficiency, suggesting that common CoW variants may produce regionally coherent hemodynamic effects in heart-healthy aging adults and represent markers of early cerebrovascular vulnerability.

Keywords

Circle of Willis cerebral arteries cerebral blood flow hemodynamics magnetic resonance imaging

Introduction

The circle of Willis (CoW) plays critical roles in maintaining stable perfusion, including collateral circulation,^1–3 dynamic flow regulation,⁴ hemodynamic stress buffering,⁵ and long-term adaptation to chronic vascular changes.^6,7 Structural variations are observed in approximately 70% of adults.⁸ Most prior work has focused on individuals with clinically evident cerebrovascular disease, revealing that posterior (PcoA) and anterior (AcoA) communicating artery variants increase risk, particularly under conditions requiring compensation (e.g. stenosis, occlusion).^9–11 There has been comparatively less emphasis on how variants affect dynamic cerebrovascular indicators. Emerging work suggests posterior CoW variants may contribute to perfusion asymmetries,^12,13 and missing multiple segments may contribute to cognitive deficits,¹⁴ suggesting threshold effects. Investigating hemodynamic correlates of CoW variation in a healthy aging cohort may clarify how subtle differences in arterial anatomy influence cerebral perfusion prior to the development of overt vascular disease.

Common CoW variants may exacerbate vulnerability to cerebrovascular dysregulation, even in otherwise healthy adults. While research on cerebral perfusion has traditionally emphasized the microcirculation and local autoregulation mechanisms, experimental evidence also highlights the influence of large cerebral arteries,^15–17 which play increasingly important roles in aging and cerebrovascular disease.¹⁸ The spatial coefficient of variation (sCoV)¹⁹ offers a valuable complement to cerebral blood flow (CBF) magnetic resonance imaging (MRI) measures, capturing clinically significant perfusion heterogeneity²⁰ that can signal early or subtle hemodynamic stress. Although single post-labeling delay (PLD) arterial spin labeling (ASL) MRI remains sensitive to both true physiological variation and arterial transit time (ATT) effects,^12,13 integrated CBF and sCoV analyses can improve differentiation of physiological and measurement contributions. Novel signatures of emerging hemodynamic insufficiency can be identified, including preserved CBF alongside elevated sCoV (reflecting active compensation and early perfusion instability), reduced CBF alongside preserved sCoV (reflecting uniform hypoperfusion or moderate risk of ATT bias), or reduced CBF alongside elevated sCoV (reflecting heterogenous vascular compromise or high risk of ATT bias).

The Vanderbilt Memory and Aging Project provides a unique opportunity to isolate the effects of common CoW variants on early hemodynamic changes in a cardiovascular-healthy and stroke-free cohort, minimizing confounding from overt disease. We hypothesize that absent or hypoplastic CoW segments, particularly PcoAs, may lead to localized hemodynamic disturbances in anatomically vulnerable regions, providing insight into subclinical pathways linking vascular anatomy to cerebrovascular and cognitive risk.

Methods

Study cohort

The Vanderbilt Memory and Aging Project²¹ (VMAP) is a longitudinal study investigating vascular health and brain aging. Between 2012 and 2014, participants were recruited into the “Legacy Cohort” and required to be ⩾60 years of age, have English proficiency, have adequate auditory and visual acuity, have a reliable study partner, and be cognitively unimpaired (CU) or meet diagnostic criteria for early mild cognitive impairment (MCI)²² or MCI.²³ Starting in 2021, the cohort was expanded (“Expansion Cohort”) to include participants ⩾50 years of age who were required to have English proficiency, have adequate auditory and visual acuity, have a reliable study partner, and be CU. Exclusion criteria for all participants included MRI contraindication; a history of neurological disease (e.g. clinical dementia, clinical stroke, head injury with loss of consciousness >5 min); heart failure; major psychiatric illness; or a systemic or terminal illness affecting follow-up examination participation. At enrollment, participants completed a comprehensive evaluation, including (but not limited to) fasting blood draw, physical examination, clinical interview, medication review, neuropsychological assessment, and multimodal brain MRI. The current study excluded participants for missing predictor, covariate, or outcome data. The protocol, guided by the Belmont Report, was approved by the Vanderbilt University Medical Center Institutional Review Board, and written informed consent was obtained prior to data collection. Data, analytic methods, and study materials can be obtained by contacting the corresponding author.

Brain MRI—CoW variants

Between 2012 and 2017, participants were scanned at the Vanderbilt University Institute of Imaging Science on a 3 T MRI system (Philips Achieva, Best, the Netherlands) with an eight-channel phased-array sensitivity encoding (SENSE) receiver head coil. In 2017, the system was upgraded to a 32-channel dStream head coil. Differences in scanner head coils were evaluated and accounted for as needed (see harmonization methods below). Time-of-flight (TOF) magnetic resonance angiography (MRA) images (TR = 6.8 ms, TE = 3.8 ms, spatial resolution = 0.35 × 0.35 × 0.50 mm³) were acquired as part of the larger multimodal neuroimaging protocol. 2D source MRA images were used to assess the patency of the CoW due to their ability to better preserve signal intensity, vessel continuity, and small or low-flow vessel visibility compared to 3D reconstructions.²⁴ All images were reviewed on two separate occasions by a board-certified neuroradiologist (LTD) blinded to clinical information. CoW communicating artery (AcoA, left PcoA, right PcoA), P1 (left P1, right P1), and A1 (left A1, right A1) segments were coded based on manual measurements of single image slices as normal (⩾0.8 mm), hypoplastic (<0.8 mm), or aplastic (invisible on MRA). All CoW segments that were coded differently between the two measurements were further assessed by a separate reviewer (CWB) using 3D MRA reconstructed images in OsiriX (OsiriX, Geneva, Switzerland) for the Legacy Cohort, which offered complementary spatial information for secondary evaluations of how vessels connect or branch. The 3D reconstructed images were used to make a final decision on whether segments with discrepancies were normal, hypoplastic, or aplastic (i.e. missing). The entire CoW was then classified as normal across all segments, hypoplastic in at least one segment, or missing in at least one segment for descriptive purposes. Six CoW variant categories were created for analytical purposes, including missing communicating arteries (AcoA missing, unilateral or bilateral PcoAs missing, AcoA and unilateral or bilateral PcoAs missing) and the reference group of otherwise patent communicating arteries (AcoA and bilateral PcoAs normal or hypoplastic); hypoplastic communicating arteries (AcoA hypoplastic, unilateral or bilateral PcoA hypoplastic, AcoA and unilateral or bilateral PcoA hypoplastic) and the reference group of otherwise normal communicating arteries (AcoA and bilateral PcoAs normal); full fetal-type posterior cerebral artery (FTP) variant (unilateral or bilateral P1 missing) and the reference group of otherwise patent P1s (bilateral P1s normal or hypoplastic); and partial FTP (unilateral or bilateral P1 hypoplasia and normal ipsilateral PcoA) and the reference group of otherwise normal P1s (bilateral P1s normal). Variants were further subdivided by unilateral and bilateral characteristics (i.e. “lateralized subgroups”) when sample sizes were sufficient (e.g. missing unilateral PcoA only, missing bilateral PcoA only, missing AcoA + unilateral PcoA, missing AcoA + bilateral PcoA).

Brain MRI—CBF and CBF sCoV

In the VMAP Legacy Cohort, T1-weighted MPRAGE images were acquired using the following parameters: TR = 8.9 ms, TE = 4.6 ms, flip angle = 5°, and spatial resolution = 1 × 1 × 1 mm³. The inversion time (TI) was session-specific, with the “shortest” viable TI selected, ranging from 916.44–978.81 ms (mean = 937.04 ms). In the VMAP Expansion Cohort, T1-weighted MPRAGE images were acquired using the following parameters: TR = 6.5 ms, TE = 2.9 ms, TI = 900 ms, FA = 9°, and spatial resolution = 1 × 1 × 1 mm³. Differences in the T1-weighted MPRAGE sequence were evaluated and accounted for as described in the harmonization methods below. T1-weighted MPRAGE images were post-processed using a NiChart_DLMUSE pipeline (DLMUSE) with parcellation of five gray matter regions of interest (ROI), including whole brain and frontal, temporal, parietal, and occipital lobes.

Standard pseudo-continuous ASL MRI with a 2D EPI readout (label duration = 1.65 s; post-labeling delay = 1.525 s; spatial resolution = 3 × 3 × 7.5 mm³; TR = 3900 ms, TE = 13 ms) assessed CBF (ml of blood/100 g tissue/min) at rest using a reproducible protocol,^21,25 enabling reliable measures of CBF and sCoV and enhanced sensitivity to subtle perfusion differences between groups. An equilibrium magnetization (M₀) image was acquired separately with identical geometry and TR = 20,000 ms. Images were corrected for motion using the 3dvolreg tool in AFNI version 23.1.10.²⁶ Skull stripping was performed using SynthStrip.²⁷ Distortion correction of M₀ and ASL images was performed using a field map-less method.²⁸ Images were slice-time-corrected, normalized by the equilibrium magnetization (M₀), and converted to absolute CBF units using recommended guidelines and participant-specific T1b (calculated from participant-specific hematocrit).²⁹ Individual control-label pairs with poor quality, including large ATT-related signal fluctuations, were excluded using ENhancement of Automated Blood fLow Estimates (ENABLE), an automated method that down-weights or removes image pairs exhibiting excessive motion, signal dropout, temporal variance, and slice-wise artifacts.³⁰ For our PLD, which fell within the range of the predefined ENABLE PLD index, relative weights for these quality measures were linearly interpolated. Sessions with mean framewise-displacement above 1 mm were excluded due to excessive motion.³¹ Final CBF maps were visually inspected and excluded for motion or labeling error artifacts, and vascular artifacts were flagged for exclusion in sensitivity analyses. CBF maps were co-registered to the anatomic T1-weighted map and standard Montreal Neurological Institute template.³² Transformation matrices were applied to CBF maps. Mean resting CBF was calculated within the gray matter ROIs described above.

Following absolute CBF calculations, the CBF spatial coefficient of variation (sCoV) was assessed as a measure of perfusion heterogeneity or stability.¹⁹ sCoV was calculated by dividing the standard deviation of the CBF signal within a given gray matter ROI by the corresponding mean CBF signal (i.e. σ (CBF_ROI)/μ (CBF_ROI) × 100). We additionally calculated a nuisance covariate, pseudo sCoV, within each ROI using a pseudo CBF image created from smoothed gray matter and white matter segmentations to adjust for partial volume-related signal variability, as previously described.³³ Partial volume corrected (PVC) CBF and pseudo sCoV maps were generated using DLMUSE segmentations to account for tissue-specific signal variability.

Brain MRI—harmonization

The present study combined baseline data across the VMAP Legacy Cohort and Expansion Cohort. To account for batch variable effects (e.g. differences in scanner coil, scanner software, T1-sequence parameters), ASL-based outputs (CBF, sCoV, PVC–CBF, pseudo sCoV) were harmonized using longitudinal ComBat³⁴ across the combined VMAP Legacy and Expansion Cohorts. T1 sequence was included as a batch variable in CBF harmonization to account for scanner-specific differences, while MPRAGE images themselves were not directly harmonized prior to ASL analysis.

Harmonization outputs were evaluated primarily via visual inspection of spaghetti plots to check for outliers, epoch-wise issues, and general patterns. This approach ensures consistent population-level adjustment for scanner and sequence effects while preserving physiologically interpretable measures for CBF and sCoV analyses. Following longitudinal harmonization, data were subset for the present cross-sectional analysis.

Analytical plan

Covariates were identified a priori for their potential to confound the analytical models, including age,^35–37 sex,³⁸ education,³⁹ race/ethnicity,^40,41 cognitive status,⁴² Framingham Stroke Risk Profile (FSRP) excluding points assigned for age,^43,44 and apolipoprotein E ε4 (APOE-ε4) carrier status (defined as positive (ε2/ε4, ε3/ε4, ε4/ε4) or negative (ε2/ε2, ε2/ε3, ε3/ε3)).⁴⁵ FSRP assigned points by sex for age, systolic blood pressure (accounting for antihypertensive medication use), diabetes mellitus, cigarette smoking, left ventricular hypertrophy, cardiovascular disease, and atrial fibrillation.^46,47

Linear regression models using ordinary least square estimates related categorical CoW variants (i.e. missing communicating arteries; missing communicating arteries by lateralized subgroups; hypoplastic communicating arteries; hypoplastic communicating arteries by lateralized subgroups; full FTP; partial FTP) to region-specific mean resting CBF (one test per model) or log-transformed sCoV (one test per model). Each model yielded a global F-test from analysis of variance (ANOVA) to assess whether CoW variant status collectively explained a significant proportion of variance in CBF or sCoV, as well as level-specific estimates comparing each CoW variant to the reference category. Models were adjusted for age, sex, race/ethnicity, education, cognitive status, FSRP (excluding points assigned for age), and APOE-ε4 status. Sensitivity analyses were performed excluding statistical outliers (defined as ⩾4 standard deviations from the sample mean) to test if these conditions accounted for significant results. Additional sensitivity analyses assessed whether partial-volume related signal variability (e.g. tissue atrophy) confounded significant results by re-running CBF models using PVC CBF outcomes, and separately, adding pseudo sCoV covariates to sCoV models. To account for multiple comparisons, the Benjamini-Hochberg procedure was used to control the false discovery rate (FDR) at a threshold of Q = 0.05.⁴⁸ FDR was chosen for its ability to balance error control and sensitivity in datasets with expected correlation among ROI measures. The number of comparisons was set to m = 10 outcomes per predictor, consistent with hypotheses (i.e. variant effect profiles across regional perfusion outcomes). This overall approach supports hypothesis-driven inference while accounting for dependence among related measures. Final sensitivity analyses applied a global FDR correction across all tests (m = 60) to assess whether observed results reflect stable, reproducible effects. All analyses were conducted using R version 3.2.3.⁴⁷

Results

Participant characteristics

The cohort was composed of 390 adults (71 ± 8 years, 55% male, 88% self-identified as non-Hispanic White), including 245 CU participants, 121 participants with MCI, and 24 with early MCI. Mean CBF was 48 ± 10 ml/100 g/min and mean sCoV was 41% ± 9% across whole brain gray matter. Sex differences existed across missing communicating artery subgroups (p = 0.002). See Table 1 for more details.

Table 1.

Participant characteristics by communicating artery variant.

Characteristic	Total (n = 390)	AcoA and PcoAs Present (n = 167)	Missing AcoA (n = 36)	Missing PcoA (n = 149)	Missing AcoA and PcoA (n = 38)	p value
Demographic
Age (years)	71 ± 8	71 ± 8	71 ± 9	72 ± 8	73 ± 7	0.34
Sex (%), male	55	50	36	66	50	0.002
Race (%), non-Hispanic White	88	87	75	90	92	0.08
Education (years)	16 ± 3	16 ± 3	16 ± 3	16 ± 3	16 ± 2	0.66
Cognitive status (%), MCI	31	34	39	27	29	0.55
APOE-ε4 (%), carrier	34	39	28	29	42	0.15
Body mass index (kg/m²)	27.4 ± 5	27.7 ± 5	26.6 ± 4	27.4 ± 5	27.1 ± 4	0.66
Framingham Stroke Risk Profile (total)	11.5 ± 4.5	11.4 ± 4.6	12.0 ± 5.2	11.4 ± 4.4	11.9 ± 3.7	0.85
Systolic blood pressure (mmHg)	138 ± 18	140 ± 18	140 ± 19	137 ± 17	133 ± 15	0.24
Anti-hypertensive medication usage (%)	50	48	47	53	53	0.79
Diabetes (%)	16	14	22	17	13	0.61
Cigarette smoking (%), current	2	2	3	1	5	0.48
Prevalent CVD (%)	4	3	0	6	8	0.21
Atrial fibrillation (%)	6	5	8	6	3	0.75
Left ventricular hypertrophy (%)	6	6	8	5	11	0.66
Intracranial volume (cm³)	1385 ± 140	1390 ± 140	1305 ± 145	1402 ± 134	1375 ± 135	0.001
Cerebral blood flow (ml/100 g/min)
Whole brain	48 ± 10	50 ± 10	48 ± 11	47 ± 10	50 ± 11	0.14
Frontal lobes	54 ± 11	55 ± 11	54 ± 11	52 ± 12	56 ± 12	0.35
Temporal lobes	40 ± 8	40 ± 8	39 ± 8	39 ± 8	41 ± 9	0.22
Parietal lobes	55 ± 13	57 ± 13	55 ± 13	53 ± 13	58 ± 15	0.11
Occipital lobes	44 ± 12	46 ± 11	45 ± 12	41 ± 12	45 ± 13	0.003
Cerebral blood flow sCoV (%)
Whole brain	41 ± 9	39 ± 8	41 ± 11	43 ± 10	41 ± 8	0.001
Frontal lobes	43 ± 7	42 ± 6	43 ± 8	44 ± 8	42 ± 6	0.21
Temporal lobes	55 ± 12	53 ± 11	55 ± 12	58 ± 12	55 ± 11	<0.001
Parietal lobes	45 ± 12	44 ± 10	45 ± 14	47 ± 13	44 ± 12	0.22
Occipital lobes	51 ± 18	48 ± 15	50 ± 19	55 ± 20	50 ± 16	0.003

AcoA: anterior communicating artery; APOE: apolipoprotein E; CVD: cardiovascular disease; MCI: mild cognitive impairment; ml/100 g/min: milliliters of blood per 100 grams of tissue per minute; PcoA: posterior communicating artery; sCoV: spatial coefficient of variation.

Descriptive statistics were calculated using mean ± standard deviation for continuous variables and frequencies for categorical variables. Between-group characteristics were statistically compared using rank-based Kruskal–Wallis test for continuous variables and Pearson’s chi-square test for categorical variables. Bolded values indicate p < 0.05. Framingham Stroke Risk Profile (minus age points) score is 6.1 ± 3.2 for the total cohort.

CoW variant prevalence

Variant CoW structures were observed among 96.2% of participants at baseline. Variant communicating arteries were observed among 91.3% of participants, including 34.1% with at least one hypoplastic communicating artery and 57.2% with at least one missing communicating artery. Among individual arteries, A1 variants were observed among 9.7% of participants (including 6.9% hypoplastic and 2.8% missing), AcoA variants among 59.3% of participants (including 40.3% hypoplastic and 19.0% missing), P1 variants among 22.1% of participants (including 14.9% hypoplastic and 8.5% missing), and PcoA variants among 77.4% of participants (including 38.2% hypoplastic and 47.9% missing). Unilateral variants were observed among 9.7% of participants for A1 segments (including 6.9% hypoplastic and 2.8% missing), 21.3% of participants for P1 segments (including 13.1% hypoplastic and 8.2% missing), and 51.5% of participants for PcoA segments (including 25.9% hypoplastic and 25.6% missing). Bilateral variants were observed among 2.1% of participants for P1 segments (including 1.8% hypoplastic and 0.3% missing) and 34.6% of participants for PcoA segments (including 12.3% hypoplastic and 22.3% missing); no bilateral A1 variants were observed. The prevalence of full FTP variants was 8.5% and partial FTP variants was 13.6%. See Table 2 for further details.

Table 2.

Prevalence of CoW variants.

Segment(s)	Typical	Hypoplastic	Missing
CoW^a	3.8% (15)	34.4% (134)	61.8% (241)
Communicating arteries^b	8.7% (34)	34.1 (133)	57.2% (223)
Individual artery
A1	90.3% (352)	6.9% (27)	2.8% (11)
Unilateral A1	–	6.9% (27)	2.8% (11)
Bilateral A1	–	0% (0)	0% (0)
AcoA	40.7% (159)	40.3% (157)	19.0% (74)
P1	77.9% (304)	14.9% (58)	8.5% (33)
Unilateral P1	–	13.1% (51)	8.2% (32)
Bilateral P1	–	1.8% (7)	0.3% (1)
PcoA	22.6% (90)	38.2% (149)	47.9% (187)
Unilateral PcoA	–	25.9% (101)	25.6% (100)
Bilateral PcoA	–	12.3% (48)	22.3% (87)
FTP
Full FTP variant	–	–	8.5% (33)
Unilateral full FTP	–	–	8.2% (32)
Bilateral full FTP	–	–	0.3% (1)
Partial FTP variant	–	13.6% (53)	–
Unilateral partial FTP	–	11.8% (46)	–
Bilateral partial FTP	–	1.8% (7)	–

A1: A1 segment of the anterior cerebral artery; AcoA: anterior communicating artery; CoW: circle of Willis; FTP: fetal-type posterior cerebral artery variant, including missing (full FTP) or hypoplastic (partial FTP) P1 segment of the posterior cerebral artery relative to a typical ipsilateral PcoA; P1: P1 segment of poster cerebral artery; PcoA: posterior communicating artery.

CoW combined structure was measured across seven cerebral arteries, including AcoA and bilateral A1, P1, and PcoA segments. Typical group indicates presence of all normal segments, hypoplastic group indicates presence of at least one hypoplastic segment (but no missing segments), and missing group indicates presence of at least one aplastic segment.

Communicating arteries of the CoW include the AcoA, left PcoA, and right PcoA.

CoW variants and CBF

Missing communicating arteries related to differences in occipital lobe CBF (p = 0.01). The association survived outlier exclusion (p = 0.01), partial volume correction (p = 0.03), and vascular artifact exclusion (p = 0.02) but was attenuated with FDR correction (p = 0.07). Level-specific estimates indicated that missing PcoAs related to lower occipital lobe CBF (β = −4.15 ml/100 g/min, p = 0.001), and other variant levels were not related (p-values > 0.43). Missing communicating arteries by lateralized subgroup related to differences in whole brain (p = 0.04), parietal (p = 0.02), and occipital lobe CBF (p = 0.001). Associations survived outlier exclusion (p-values < 0.04) and vascular artifact exclusion (p-values < 0.04); while parietal and occipital lobe CBF associations survived FDR correction (FDR-adjusted p-values < 0.05) and partial volume correction (p-values < 0.04), whole brain CBF associations were attenuated with FDR correction (p = 0.07) and partial volume correction (p = 0.09). Level-specific estimates indicated that missing bilateral PcoAs related to lower whole brain (β = −3.93 ml/100 g/min, p = 0.01), parietal (β = −4.87 ml/100 g/min, p = 0.01), and occipital lobe CBF (β = −6.54 ml/100 g/min, p = 0.0001); missing AcoA and bilateral PcoAs related to higher parietal lobe CBF (β = 6.47 ml/100 g/min, p < 0.05); and other variant levels were not related (p-values > 0.21). Neither hypoplastic communicating arteries (p-values > 0.40), hypoplastic communicating arteries by lateralized subgroup (p-values > 0.32), full FTP (p-values > 0.70), nor partial FTP variants (p-values > 0.22) related to any of the CBF ROIs. In sensitivity analyses using a global FDR correction, missing communicating arteries by lateralized subgroup remained related to differences in occipital lobe CBF (p = 0.03). See Table 3 for details and Figure 1 for illustrations.

Figure 1.

Resting CBF across missing communicating artery variants. Occipital lobe CBF differed significantly across missing communicating artery variants (a), such that missing PcoAs related to lower occipital lobe CBF compared to the patent AcoA and PcoA reference level, indicating hypoperfusion. Whole brain, parietal, and occipital lobe CBF differed significantly across missing communicating artery variants when examined by unilateral or bilateral characteristics (b), such that missing bilateral PcoAs related to lower CBF levels compared to the patent AcoA and PcoA reference level. In addition, within parietal lobe CBF models, concurrently missing AcoA and bilateral PcoAs marginally related to higher CBF compared to the patent AcoA and PcoA reference level.

Table 3.

CoW variants and resting CBF associations.

Categorical predictor	n	Whole brain CBF			Frontal lobe CBF			Parietal lobe CBF			Temporal lobe CBF			Occipital lobe CBF
Categorical predictor	n	β	CI	p value	β	CI	p value	β	CI	p value	β	CI	p value	β	CI	p value
Missing AcoA or PcoA	390			0.08			0.16			0.06			0.14			0.01 ^a,b
Normal or hypoplastic AcoA + PcoAs (ref.)	167
Missing AcoA only	36	−1.01	−4.63 to 2.61	0.58	−0.48	−4.49 to 3.53	0.81	−1.82	−6.38 to 2.75	0.43	−1.02	3.82 to 1.79	0.48	−1.64	−5.75 to 2.46	0.43
Missing PcoA only	149	−2.53	−4.78 to −0.30	0.03	−2.07	−4.55 to 0.41	0.10	−3.16	−5.98 to −0.33	0.03	−1.73	−3.47 to 0.002	0.05	−4.15	−6.69 to −1.61	0.001
Missing AcoA + PcoA	38	1.26	−2.27 to 4.79	0.48	2.01	−1.89 to 5.92	0.31	2.07	−2.38 to 6.52	0.36	0.90	−1.84 to 3.63	0.52	0.15	−3.85 to 4.15	0.94
Missing AcoA or PcoA, lateralized	390			0.04 ^a			0.13			0.02 ^a,b,c			0.10			0.001 ^a,b,c
Normal or hypoplastic AcoA + PcoAs (ref.)	167
Missing AcoA only	36	−1.03	−4.64 to 2.57	0.57	−0.50	−4.50 to 3.50	0.81	−1.84	−6.38 to 2.70	0.43	−1.03	−3.83 to 1.76	0.47	−1.68	−5.75 to 2.39	0.42
Missing unilateral PcoA only	77	−1.27	−3.98 to 1.45	0.36	−0.97	−3.98 to 2.04	0.53	−1.61	−5.02 to 1.81	0.36	−0.77	−2.88 to 1.33	0.47	−1.97	−5.03 to 1.09	0.21
Missing bilateral PcoA only	72	−3.93	−6.72 to −1.15	0.01	−3.29	−6.38 to −0.20	0.04	−4.87	−8.38 to −1.36	0.01	−2.78	−4.95 to −0.62	0.01	−6.54	−9.69 to −3.40	0.0001
Missing AcoA + unilateral PcoA	23	−0.64	−5.00 to 3.71	0.77	0.13	−4.70 to 4.96	0.96	−1.00	−6.48 to 4.48	0.72	−0.15	−3.53 to 3.22	0.93	−2.18	−7.10 to 2.74	0.38
Missing AcoA + bilateral PcoA	15	4.16	−1.10 to 9.41	0.12	4.87	−0.96 to 10.70	0.10	6.74	0.12–13.36	<0.05	2.49	−1.58 to 6.57	0.23	3.70	−2.24 to 9.63	0.22
Hypoplastic AcoA or PcoA	167			0.62			0.60			0.54			0.77			0.40
Normal AcoA + PcoAs (ref.)	34
Hypoplastic AcoA only	36	−2.93	−7.52 to 1.67	0.21	−3.36	−8.37 to 1.65	0.19	−4.33	−10.16 to 1.50	0.14	−1.73	−5.32 to 1.86	0.34	−3.15	−8.35 to 2.04	0.23
Hypoplastic PcoA only	48	−2.18	−6.48 to 2.12	0.32	−2.11	−6.80 to 2.57	0.37	−2.30	−7.75 to 3.15	0.41	−1.36	−4.72 to 2.00	0.43	−3.68	−8.53 to 1.18	0.14
Hypoplastic AcoA + PcoA	49	−1.47	−5.81 to 2.86	0.50	−1.30	−6.03 to 3.43	0.59	−2.10	−7.60 to 3.40	0.45	−0.59	−3.97 to 2.80	0.73	−3.80	−8.70 to 1.10	0.13
Hypoplastic AcoA or PcoA, lateralized	167			0.40			0.32			0.34			0.58			0.52
Normal AcoA + PcoAs (ref.)	34
Hypoplastic AcoA only	36	−2.83	−7.41 to 1.75	0.22	−3.24	−8.22 to 1.74	0.20	−4.20	−10.01 to 1.60	0.15	−1.66	−5.24 to 1.93	0.36	−3.08	−8.29 to 2.13	0.24
Hypoplastic unilateral PcoA only	25	−4.00	−9.10 to 1.10	0.12	−4.38	−9.93 to 1.16	0.12	−4.55	−11.01 to 1.91	0.17	−2.62	−6.61 to 1.37	0.20	−4.83	−10.63 to 0.97	0.10
Hypoplastic bilateral PcoA only	23	−0.23	−5.37 to 4.92	0.93	0.32	−5.29 to 5.91	0.91	0.12	−6.40 to 6.64	0.97	0.003	−4.03 to 4.03	0.99	−2.42	−8.28 to 3.43	0.42
Hypoplastic AcoA + unilateral PcoA	29	−0.06	−4.93 to 4.81	0.98	0.28	−5.02 to 5.58	0.92	−0.17	−6.34 to 6.00	0.96	0.40	−3.41 to 4.21	0.84	−2.72	−8.26 to 2.82	0.33
Hypoplastic AcoA + bilateral PcoA	20	−3.63	−9.06 to 1.81	0.19	−3.71	−9.62 to 2.20	0.22	−5.02	−11.90 to 1.86	0.15	−2.08	−6.34 to 2.17	0.33	−5.43	−11.62 to 0.75	0.08
Full FTP	390			0.89			0.70			0.77			0.80			0.85
Normal or hypoplastic P1s (ref.)	357
Full FTP	33	0.26	−3.38 to 3.89	0.89	0.79	−3.23 to 4.80	0.70	−0.69	−5.27 to 3.89	0.77	0.36	−2.45 to 3.17	0.80	0.41	−3.74 to 4.56	0.85
Partial FTP	357			0.40			0.22			0.32			0.25			0.64
Normal P1s (ref.)	304
Partial FTP	53	−1.25	−4.19 to 1.69	0.40	−2.02	−5.25 to 1.22	0.22	−1.89	−5.60 to 1.82	0.32	−1.33	−3.59 to 0.93	0.25	0.80	−2.56 to 4.15	0.64

AcoA: anterior communicating artery; CBF: cerebral blood flow; CI: 95% confidence interval; CoW: circle of Willis; FTP: fetal-type posterior cerebral artery variant; P1s: P1 segments of poster cerebral artery (left and right); PcoAs: posterior communicating arteries (left and right).

Bolded values indicate p < 0.05.

Result survived exclusion for statistical outliers.

Result survived partial volume correction.

Result survived FDR correction for multiple hypothesis testing.

CoW variants and CBF sCoV

Missing communicating arteries related to differences in whole brain (p = 0.03), temporal (p = 0.01), and occipital lobe sCoV (p = 0.04). Associations survived outlier exclusion (p-values < 0.05) but were attenuated with FDR correction (p-values > 0.06); temporal and occipital lobe sCoV associations survived partial volume correction (p-values < 0.05) but whole brain sCoV did not (p = 0.05); and temporal lobe associations survived vascular artifact exclusion (p = 0.03) but whole brain and occipital lobe sCoV did not (p-values > 0.07). Level-specific estimates indicated that missing PcoAs related to higher whole brain (p = 0.002), temporal (p = 0.001), and occipital lobe sCoV (p = 0.01), and other variant levels were not related (p-values > 0.17). Missing communicating arteries by lateralized subgroup related to differences in whole brain (p = 0.002), temporal (p = 0.001), and occipital lobe sCoV (p = 0.01). Associations survived FDR correction (p-values < 0.01), outlier exclusion (p-values < 0.01), partial volume correction (p-values < 0.005), and vascular artifact exclusion (p-values < 0.01). Level-specific estimates indicated that missing bilateral PcoAs related to higher whole brain (p < 0.0001), temporal (p < 0.0001), and occipital lobe sCoV (p = 0.0003). Neither hypoplastic communicating arteries (p-values > 0.22), hypoplastic communicating arteries by lateralized subgroup (p-values > 0.32), full FTP (p-values > 0.60), nor partial FTP variants (p-values > 0.82) related to sCoV ROIs. In sensitivity analyses using a global FDR correction, missing communicating arteries by lateralized subgroup remained related to differences in whole brain (p = 0.03) and temporal lobe sCoV (p = 0.03). See Table 4 for details and Figure 2 for illustrations.

Table 4.

CoW variants and CBF sCoV associations.

Categorical predictor	n	Whole brain sCoV			Frontal lobe sCoV			Parietal lobe sCoV			Temporal lobe sCoV			Occipital lobe sCoV
Categorical predictor	n	β	CI	p value	β	CI	p value	β	CI	p value	β	CI	p value	β	CI	p value
Missing AcoA or PcoA	390			0.03 ^a			0.54			0.36			0.01 ^a,b			0.04 ^a,b
Normal or hypoplastic AcoA + PcoAs (ref.)	167
Missing AcoA only	36	0.03	−0.03 to 0.09	0.39	0.01	−0.03 to 0.06	0.58	0.02	−0.06 to 0.10	0.60	0.04	−0.02 to 0.11	0.17	0.04	−0.06 to 0.14	0.48
Missing PcoA only	149	0.05	0.01–0.10	0.003	0.02	−0.01 to 0.05	0.27	0.03	−0.02 to 0.07	0.26	0.07	0.03–0.10	0.001	0.09	0.03–0.15	0.01
Missing AcoA + PcoA	38	0.02	−0.04 to 0.09	0.46	−0.01	−0.06 to 0.03	0.60	−0.04	−0.11 to 0.04	0.34	0.03	−0.03 to 0.09	0.35	0.02	−0.08 to 0.11	0.73
Missing AcoA or PcoA, lateralized	390			0.002 ^a,b,c			0.34			0.05			0.001 ^a,b,c			0.01 ^a,b,c
Normal or hypoplastic AcoA + PcoAs (ref.)	167
Missing AcoA only	36	0.03	−0.03 to 0.09	0.37	0.01	−0.03 to 0.06	0.57	0.02	−0.05 to 0.10	0.59	0.04	−0.02 to 0.11	0.16	0.04	−0.06 to 0.14	0.46
Missing unilateral PcoA only	77	0.02	−0.03 to 0.07	0.41	−0.001	−0.04 to 0.04	0.96	−0.001	−0.06 to 0.10	0.59	0.03	−0.02 to 0.08	0.21	0.04	−0.04 to 0.11	0.30
Missing bilateral PcoA only	72	0.10	0.05–0.15	<0.0001	0.04	−0.001 to 0.07	0.06	0.06	−0.002 to 0.12	0.05	0.11	0.05–0.15	<0.0001	0.14	0.07–0.22	0.0003
Missing AcoA + unilateral PcoA	23	0.06	−0.02 to 0.13	0.13	0.002	−0.06 to 0.06	0.94	0.02	−0.07 to 0.12	0.60	0.08	0.005–0.15	<0.05	0.08	−0.04 to 0.20	0.18
Missing AcoA + bilateral PcoA	15	−0.03	−0.12 to 0.06	0.51	−0.04	−0.11 to 0.04	0.33	−0.13	−0.24 to −0.02	0.02	−0.04	−0.13 to 0.05	0.38	−0.08	−0.22 to 0.06	0.27
Hypoplastic AcoA or PcoA	167			0.59			0.22			0.23			0.22			0.23
Normal AcoA + PcoAs (ref.)	34
Hypoplastic AcoA only	36	0.03	−0.04 to 0.10	0.44	0.01	−0.05 to 0.06	0.79	0.02	−0.07 to 0.10	0.73	0.03	−0.05 to 0.11	0.38	0.02	−0.09 to 0.14	0.68
Hypoplastic PcoA only	48	0.04	−0.03 to 0.10	0.29	0.01	−0.05 to 0.06	0.85	−0.001	−0.08 to 0.08	0.98	0.04	−0.03 to 0.11	0.30	0.06	−0.05 to 0.16	0.28
Hypoplastic AcoA + PcoA	49	0.05	−0.02 to 0.11	0.18	0.01	−0.05 to 0.06	0.81	0.03	−0.06 to 0.11	0.55	0.04	−0.03 to 0.12	0.25	0.09	−0.01 to 0.20	0.08
Hypoplastic AcoA or PcoA, lateralized	167			0.79			0.32			0.32			0.32			0.32
Normal AcoA + PcoAs (ref.)	34
Hypoplastic AcoA only	36	0.03	−0.04 to 0.10	0.45	0.01	−0.05 to 0.06	0.80	0.01	−0.07 to 0.10	0.75	0.03	−0.04 to 0.11	0.40	0.02	−0.09 to 0.14	0.69
Hypoplastic unilateral PcoA only	25	0.05	−0.03 to 0.12	0.23	0.01	−0.04 to 0.08	0.53	0.01	−0.09 to 0.11	0.86	0.04	−0.04 to 0.13	0.32	0.06	−0.06 to 0.19	0.32
Hypoplastic bilateral PcoA only	23	0.02	−0.06 to 0.10	0.58	−0.01	−0.07 to 0.05	0.75	−0.01	−0.11 to 0.09	0.81	0.03	−0.06 to 0.12	0.46	0.05	−0.08 to 0.18	0.43
Hypoplastic AcoA + unilateral PcoA	29	0.04	−0.04 to 0.12	0.31	0.01	−0.05 to 0.07	0.84	−0.01	−0.09 to 0.10	0.89	0.03	−0.05 to 0.11	0.47	0.08	−0.04 to 0.20	0.19
Hypoplastic AcoA + bilateral PcoA	20	0.06	−0.03 to 0.14	0.21	0.01	−0.06 to 0.07	0.82	0.05	−0.05 to 0.16	0.32	0.06	−0.03 to 0.15	0.20	0.12	−0.02 to 0.25	0.09
Full FTP	390			0.78			0.63			0.94			0.66			0.84
Normal or hypoplastic P1s (ref.)	357
Full FTP	33	−0.01	−0.07 to 0.05	0.78	0.01	−0.04 to 0.06	0.63	0.003	−0.07 to 0.08	0.94	0.01	−0.05 to 0.08	0.66	0.01	−0.09 to 0.11	0.84
Partial FTP	357			0.53			0.97			0.96			0.36			0.27
Normal P1s (ref.)	304
Partial FTP	53	−0.02	−0.07 to 0.04	0.53	0.001	−0.04 to 0.04	0.97	0.002	−0.06 to 0.06	0.96	−0.02	−0.08 to 0.03	0.36	−0.05	−0.13 to 0.04	0.27

Bolded values indicate p < 0.05.

Result survived exclusion for statistical outliers.

Result survived partial volume correction.

Result survived FDR correction for multiple hypothesis testing.

Figure 2.

CBF sCoV across missing communicating artery variants. Whole brain, temporal, and occipital lobe sCoV differed significantly across missing communicating artery variants (a), such that the missing PcoAs related to higher sCoV compared to the patent AcoA and PcoA reference level, indicating higher perfusion instability. The same sCoV regions differed significantly across missing communicating artery variants when examined by unilateral or bilateral characteristics (b), such that missing bilateral PcoAs related to higher whole brain, temporal, and occipital lobe sCoV compared to the patent AcoA and PcoA reference level. In addition, within temporal lobe sCoV models, concurrently missing AcoA and unilateral PcoA marginally related to higher temporal lobe sCoV compared to the patent AcoA and PcoA reference level.

Discussion

Among community-dwelling middle-aged and older adults free of clinical stroke at baseline, structural variants of the CoW communicating arteries related to greater cerebral perfusion variability and modest perfusion reductions. Bilateral absence of the PcoAs was linked to significantly higher whole brain sCoV (most notably in the temporal lobes), lower occipital lobe CBF alongside higher occipital sCoV, and moderately lower parietal lobe CBF. Results reveal a directionally consistent pattern in which the loss of posterior collateral pathways appears to reduce both perfusion stability and magnitude in anatomically vulnerable regions. It is noteworthy that these effects were statistically independent of vascular risk factors (e.g. hypertension, diabetes, smoking), indicating that even among relatively healthy aging adults with minimal cardiovascular disease, prevalent PcoA absence may still relate to compromised cerebral hemodynamics or emerging risk.

The association of PcoA absence with elevated sCoV may reflect posterior perfusion instability and the potential compensatory limits of collaterals. Missing PcoAs limit posterior-to-anterior compensation capacity, and occipital and temporal lobes may be particularly vulnerable compared to anterior regions.^49,50 Loss of PcoA support may shift perfusion toward longer, indirect, or unstable collateral routes, which may promote preserved but unstable blood flow (as observed in temporal lobes), compromised and unstable flow (as observed in occipital lobes), or overall prolonged ATT. In temporal lobes, mixed arterial inputs may enable compensatory adjustments when posterior collaterals are compromised, increasing blood flow heterogeneity despite its overall preservation. In contrast, sCoV may be more tightly coupled to CBF in occipital lobe regions with distal arterial positioning and limited collateral support, which increase compensation failure risk. Elevated sCoV patterns may thus serve as sensitive markers of posterior circulation stress, reflecting subtle compensation adjustments or emerging vascular insufficiency.¹⁹ Importantly, the link between PcoA absence and higher sCoV may also reflect ATT-related measurement bias, particularly with our single-PLD ASL methodology. Definitive separation of ATT effects from true perfusion deficits would require multi-PLD ASL or PET. Findings suggest that the preservation of regional CBF may depend on the presence and redundancy of collateral pathways, with temporal lobes demonstrating compensatory capacity and occipital lobes susceptible to compromised flow.

Observed CBF reductions were modest (e.g. −9% in parietal lobe) compared with larger sCoV increases (e.g. +20% in temporal lobes), consistent with mild hypoperfusion in vulnerable territories. Parietal regions, largely supplied by the middle cerebral artery, may offer higher resilience to posterior circulation variants and the capacity to redistribute flow uniformly. Although missing PcoAs may affect perfusion to parietal regions with minor posterior inputs (e.g. watershed areas), most flow likely still arrives via well-established anterior arteries, maintaining uniform distribution. While ATT-related CBF underestimation remains possible even with stable sCoV, parietal lobe patterns of modest CBF reductions with preserved sCoV represent strong candidates for true hypoperfusion and stable flow redistribution rather than perfusion instability or collateral failure. Taken together, results suggest region-dependent hemodynamic consequences of PcoA absence, wherein parietal lobes maintain effective compensation with modestly reduced but uniform flow, temporal lobes demonstrate greater flow variability and likely reliance on indirect collateral routes, and occipital lobes appear the most hemodynamically vulnerable.

As expected, most CoW variants were observed among communicating arteries (91.3%), including missing PcoAs (47.9%), hypoplastic AcoA (40.3%), and hypoplastic PcoAs (38.2%). The high prevalence of PcoA variants, especially absence, likely reflects this segment’s vulnerability to regression, collapse, or stenosis. Most existing evidence suggests CoW configurations are established early,⁵¹ but age-related changes (e.g. arterial stiffening, hypertension, atherosclerosis, microvascular rarefaction, reduced autoregulatory capacity) may also promote variant emergence, particularly since small caliber vessels are inherently fragile.^52–54 Given our cohort’s healthy cardiovascular profile, variants likely represent baseline collateral anatomy, although age-related modification is still possible. If congenital, CoW variants may reflect inherent vascular configurations that have shaped lifelong resilience to hemodynamic stress, whereas variants acquired from cumulative degenerative processes may relate to emerging signatures of flow reorganization. Notably, while our TOF MRA imaging parameters offer relatively high-resolution for small artery detection, this method may still overestimate AcoA absence compared to previous studies.⁵⁵ PcoA or AcoA absence may thus arise from any combination of congenital differences, cumulative pathology, and methodological limitations.

Several subthreshold findings demonstrated notable effect sizes with unexpected regional or directional trends, likely reflecting small subgroup sizes, regional measurement variability, or physiological heterogeneity. These associations nonetheless provide insight into regional vulnerability and physiologically plausible compensation patterns (e.g. reduced anterior or watershed perfusion capacity when anterior-to-posterior pathways are eliminated⁹; a shift in perfusion reliance to distal pathways when central CoW redundancy is lost^2,56), which could offer hypotheses for future studies with larger samples or higher-sensitivity imaging. Given limited statistical power and unstable estimates, these findings should be interpreted cautiously and do not provide definitive evidence of physiological effects.

Our study has several limitations. Although TOF MRA likely captured most hemodynamically relevant segments,^57,58 very small segments still capable of hemodynamic function⁵⁹ may still be missed or misclassified, creating heterogeneous subgroups that mask true effects. Combined ACoA/PCoA modeling may also reduce anatomical specificity and conflate distinct effects. Single-PLD ASL with a one-compartment model may underestimate CBF in regions with delayed ATT due to late blood arrival or the label decaying in transit. We are thus limited in our ability to distinguish true hypoperfusion from ATT artifact. While our rigorous quality-control procedures (e.g. motion exclusion, ENABLE filtering, vascular artifact screening), region-specific comparative analyses, and usage of a heart-healthy cohort reduce the risk of ATT-related bias, residual confounding cannot be ruled out. Although statistically significant, most observed associations were small relative to natural variability, limiting predictive or diagnostic value at the individual level. Nevertheless, findings may still be meaningful in the context of cardiovascular or demographic risk factors and offer new insights into early regional vulnerability and collateral support limits. Small effects are more susceptible to noise, measurement error, or unmeasured confounders, and residual confounding remains possible; subclinical cerebrovascular disease or physiological fluctuations (e.g. end-tidal CO₂) may contribute variability. Finally, the cross-sectional design and small subgroup sizes for rarer variants limit causal inference. These factors may attenuate or obscure subtle effects, though observed patterns remain consistent and biologically plausible.

Despite these limitations, this study has several strengths. Our extensive MRA protocol, including double review by a board-certified neuroradiologist (LTD) and reconciliation by a secondary reviewer (CWB), reduces CoW measurement error. Classifications using quantitative thresholds (0.8 mm) increase sensitivity to subtle vessel variations and facilitate reproducibility. Comprehensive variant coding, capturing unilateral/bilateral, single-segment, and co-occurring variations (e.g. FTP) provides novel, fine-grained arterial specificity over global CoW patency approaches. Studying healthy older adults with low cardiovascular risk offers insight into early hemodynamic changes before overt pathology, minimizing confounding by major stenosis, stroke, or collateral failure, and enables isolation of anatomical variant effects. Concurrent CBF and sCoV analyses provide complementary information on perfusion magnitude and efficiency, revealing early-stage hemodynamic dysregulation patterns. While findings are broadly relevant (e.g. community-dwelling adults, high prevalence variants⁶⁰), our cohort may not fully reflect the general population, as participants were predominantly White (non-Hispanic), formally educated, and relatively healthy; associations may be stronger in higher-risk groups. Lastly, our adjustments for multiple comparisons and potential confounds (e.g. statistical outliers, cardiovascular risk factors, partial volume effects, vascular artifact) increase the robustness of findings. Future research is warranted to replicate findings, expand CoW characterizations, clarify mechanisms of early stage hemodynamic insufficiency (e.g. VEASL MRI), and assess the clinical significance of our findings.

In conclusion, this study demonstrates that structural variants of the CoW, particularly bilateral PcoA absence, exert regionally coherent effects on cerebral perfusion dynamics in community-dwelling, heart-healthy aging adults. Bilateral variants produced larger effects than unilateral, suggesting a dose-like relationship. Our findings expand frameworks for interpreting ASL perfusion heterogeneity in healthy aging. Combined sCoV and CBF profiles revealed region-specific patterns of effective flow redistribution, emerging hemodynamic stress, and potential measurement artifact risk. Over time, the subtle perfusion changes observed may increase vulnerability to age-related cerebrovascular dysfunction and cognitive decline.⁶¹ However, future studies are needed to disentangle true perfusion deficits from transit time related confounds. These results provide novel insight into the link between common large artery variants and hemodynamic dysregulation, highlighting early hemodynamic vulnerabilities that may precede overt cerebrovascular disease and provide a framework for assessing subclinical vascular risk in healthy aging.

Footnotes

Acknowledgements

The authors would like to thank the dedicated Vanderbilt Memory and Aging Project participants, their loved ones, and our devoted staff and trainees who contributed to the recruitment, screening, and enrollment of the baseline cohort.

Non-standard abbreviations and acronyms

AcoA = anterior communicating artery, APOE = apolipoprotein E, ASL = arterial spin labeling, ATT = arterial transit time, CBF = cerebral blood flow, CMR = cardiac magnetic resonance imaging, CVD = cerebrovascular disease, CoW = circle of Willis, CI = 95% confidence interval, CU = cognitively unimpaired, DLMUSE = NiChart_DLMUSE pipeline; ENABLE = ENhancement of Automated Blood fLow Estimates, FSRP = Framingham Stroke Risk Profile, FTP = fetal-type posterior cerebral artery variant, PLD = post-labeling delay, PVC = partial volume corrected, VMAP = Vanderbilt Memory and Aging Project, MCI = mild cognitive impairment, MRA = magnetic resonance angiography, MRI = magnetic resonance imaging, PcoA = posterior communicating artery, ROI = region of interest.

Author contributions

FEC and ALJ conceived and designed the study. CWB, LTD, FEC, WHR, KRP, NS, TBJ, and TJH, curated and managed the imaging and clinical data. PZ, DL, and FEC performed the statistical analyses. FEC interpreted the results and drafted the manuscript. PZ created visual representations of results. ALJ supervised research activities and acquired funding and resources. All authors reviewed, provided critical revisions, and approved the final manuscript.

Funding

The authors disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This research was supported by T32-AG058524 (FEC), Howard Hughes Medical Institute James H Gilliam Fellowship for Advanced Study (FEC), F31-AG059345 (FEC), Alzheimer’s Association IIRG-08-88733 (ALJ), R01-AG034962 (ALJ), R01-NS100980 (ALJ), K24-AG046373 (ALJ), F31-AG066358 (CWB), T32-AG058524 (CWB), UL1-TR000445 and UL1-TR002243 (Vanderbilt Clinical Translational Science Award), S10-OD023680 (Vanderbilt’s High-Performance Computer Cluster for Biomedical Research), S10-OD021771 (VUIIS Center for Human Imaging), P20-AG068082 (Vanderbilt Alzheimer’s Disease Research Center), Richard Eugene Hickman Alzheimer’s Disease Research Endowment, and the Vanderbilt Memory and Alzheimer’s Center.

Declaration of conflicting interests

The authors declared the following potential conflicts of interest with respect to the research, authorship, and/or publication of this article: TJH serves on the Scientific Advisory Board for Vivid Genomics, Deputy Editor for Alzheimer’s and Dementia: TRCI, and Senior Associate Editor for Alzheimer’s and Dementia. LTD serves as a consultant for Nashville Biosciences, LLC. ALJ has served as an advisor for Lantheus — Diagnostic and Therapeutic Innovations. DL is a member of the Data Monitoring Committee for a drug study sponsored by Regeneron Pharmaceuticals. All remaining authors have nothing to disclose.

ORCID iD

Angela L Jefferson

References

Hartkamp

van Der Grond

van Everdingen

, et al. Circle of Willis collateral flow investigated by magnetic resonance angiography. Stroke 1999; 30: 2671–2678.

Liu

Zhu

Hou

, et al. Study on the correlation between the circle of Willis structure and collateral circulation in bilateral carotid artery occlusion. Neurol Sci 2021; 42: 5335–5342.

Romero

Pikula

Nguyen

, et al. Cerebral collateral circulation in carotid artery disease. Curr Cardiol Rev 2009; 5: 279–288.

Guo

Z-N

Sun

Liu

, et al. The impact of variational primary collaterals on cerebral autoregulation. Front Physiol 2018; 9: 759.

Vrselja

Brkic

Mrdenovic

, et al. Function of circle of Willis. J Cereb Blood Flow Metab 2014; 34: 578–584.

van Hespen

Mackaaij

Waas

ISE

, et al. Arterial remodeling of the intracranial arteries in patients with hypertension and controls: a postmortem study. Hypertension 2021; 77: 135–146.

Ramirez-Velandia

Mensah

Salih

, et al. Patterns of dynamic adaptability of the circle of Willis in response to major branch artery coverage with a flow diverter. World Neurosurg 2025; 193: 1065–1075.

Jones

Castanho

Bazira

, et al. Anatomical variations of the circle of Willis and their prevalence, with a focus on the posterior communicating artery: a literature review and meta-analysis. Clin Anat 2021; 34: 978–990.

Schomer

Marks

Steinberg

, et al. The anatomy of the posterior communicating artery as a risk factor for ischemic cerebral infarction. N Engl J Med 1994; 330: 1565–1570.

10.

Chuang

Y-M

Liu

C-Y

Pan

P-J

, et al. Posterior communicating artery hypoplasia as a risk factor for acute ischemic stroke in the absence of carotid artery occlusion. J Clin Neurosci 2008; 15: 1376–1381.

11.

Zhang

X-J

Gao

B-L

Hao

W-L

, et al. Presence of anterior communicating artery aneurysm is associated with age, bifurcation angle, and vessel diameter. Stroke 2018; 49: 341–347.

12.

Barkeij Wolf

JJH

Foster-Dingley

Moonen

JEF

, et al. Unilateral fetal-type circle of Willis anatomy causes right-left asymmetry in cerebral blood flow with pseudo-continuous arterial spin labeling: a limitation of arterial spin labeling-based cerebral blood flow measurements? J Cereb Blood Flow Metab 2016; 36: 1570–1578.

13.

Noorbakhsh

Farid

Bolar

DS.

Apparent posterior cerebral artery territory perfusion asymmetry on arterial spin labeling MRI is a common non-pathologic finding in patients with a unilateral fetal posterior cerebral artery. Neuroradiology 2022; 64: 513–520.

14.

Okada

Tabara

Kohara

, et al. Morphologic variation of the circle of Willis is associated with cognitive impairment in a general population: the J-SHIPP study. J Stroke Cerebrovasc Dis 2025; 34: 108287.

15.

Kontos

Wei

Navari

, et al. Responses of cerebral arteries and arterioles to acute hypotension and hypertension. Am J Physiol 1978; 234: H371–H383.

16.

Harper

Bohlen

Rubin

MJ.

Arterial and microvascular contributions to cerebral cortical autoregulation in rats. Am J Physiol 1984; 246: H17–H24.

17.

Szabo

Rosengarten

Juhasz

, et al. Effect of non-steroid anti-inflammatory drugs on neurovascular coupling in humans. J Neurol Sci 2014; 336: 227–231.

18.

Faraci

Heistad

DD.

Regulation of large cerebral arteries and cerebral microvascular pressure. Circ Res 1990; 66: 8–17.

19.

Mutsaerts

Petr

Václavů

, et al. The spatial coefficient of variation in arterial spin labeling cerebral blood flow images. J Cereb Blood Flow Metab 2017; 37: 3184–3192.

20.

Morgan

Melzer

Roberts

, et al. Spatial variation of perfusion MRI reflects cognitive decline in mild cognitive impairment and early dementia. Sci Rep 2021; 11: 23325.

21.

Jefferson

Cambronero

Liu

, et al. Higher aortic stiffness is related to lower cerebral blood flow and preserved cerebrovascular reactivity in older adults. Circulation 2018; 138: 1951–1962.

22.

Aisen

Petersen

Donohue

, et al. Clinical core of the Alzheimer’s disease neuroimaging initiative: progress and plans. Alzheimers Dement 2010; 6: 239–246.

23.

Albert

DeKosky

Dickson

, et al. The diagnosis of mild cognitive impairment due to Alzheimer’s disease: recommendations from the National Institute on Aging–Alzheimer’s Association workgroups on diagnostic guidelines for Alzheimer’s disease. Alzheimers Dement 2011; 7: 270–279.

24.

Korogi

Takahashi

Nakagawa

, et al. Intracranial vascular stenosis and occlusion: MR angiographic findings. AJNR Am J Neuroradiol 1997; 18: 135–143.

25.

Jefferson

Liu

Gupta

, et al. Lower cardiac index levels relate to lower cerebral blood flow in older adults. Neurology 2017; 89: 2327–2334.

26.

Cox

RW.

AFNI: software for analysis and visualization of functional magnetic resonance neuroimages. Comput Biomed Res 1996; 29: 162–173.

27.

Hoopes

Mora

Dalca

, et al. SynthStrip: skull-stripping for any brain image. Neuroimage 2022; 260: 119474.

28.

Schilling

Blaber

Huo

, et al. Synthesized b0 for diffusion distortion correction (Synb0-DisCo). Magn Reson Imaging 2019; 64: 62–70.

29.

Alsop

Detre

Golay

, et al. Recommended implementation of arterial spin-labeled perfusion MRI for clinical applications: a consensus of the ISMRM perfusion study group and the European consortium for ASL in dementia. Magn Reson Med 2015; 73: 102–116.

30.

Shirzadi

Stefanovic

Chappell

, et al. Enhancement of automated blood flow estimates (ENABLE) from arterial spin-labeled MRI. J Magn Reson Imaging 2018; 47: 647–655.

31.

Adebimpe

Bertolero

Dolui

, et al. ASLPrep: a platform for processing of arterial spin labeled MRI and quantification of regional brain perfusion. Nat Methods 2022; 19: 683–686.

32.

Jenkinson

Beckmann

Behrens

TEJ

, et al. FSL. Neuroimage 2012; 62: 782–790.

33.

Chappell

Kirk

Craig

, et al. BASIL: a toolbox for perfusion quantification using arterial spin labelling. Imaging Neurosci 2023; 1: 1–16.

34.

Beer

Tustison

Cook

, et al. Longitudinal ComBat: a method for harmonizing longitudinal multi-scanner imaging data. Neuroimage 2020; 220: 117129.

35.

Zhang

Gordon

Goldberg

TE.

Cerebral blood flow measured by arterial spin labeling MRI at resting state in normal aging and Alzheimer’s disease. Neurosci Biobehav Rev 2017; 72: 168–175.

36.

Juttukonda

Almaktoum

, et al. Characterizing cerebral hemodynamics across the adult lifespan with arterial spin labeling MRI data from the Human Connectome Project-Aging. Neuroimage 2021; 230: 117807.

37.

Tomoto

Khan

, et al. Cerebral blood flow and cerebrovascular resistance across the adult lifespan: a multimodality approach. J Cereb Blood Flow Metab 2023; 43: 962–976.

38.

Aanerud

Borghammer

Rodell

, et al. Sex differences of human cortical blood flow and energy metabolism. J Cereb Blood Flow Metab 2017; 37: 2433–2440.

39.

Chiu

N-T

Lee

B-F

Hsiao

, et al. Educational level influences regional cerebral blood flow in patients with Alzheimer’s disease. J Nucl Med 2004; 45: 1860–1863.

40.

Moonen

JEF

Nasrallah

Detre

, et al. Race, sex, and mid-life changes in brain health: cardia MRI substudy. Alzheimers Dement 2022; 18: 2428–2437.

41.

Kurian

Cardarelli

KM.

Racial and ethnic differences in cardiovascular disease risk factors: a systematic review. Ethn Dis 2007; 17: 143–152.

42.

Binnewijzend

MAA

Kuijer

JPA

Benedictus

, et al. Cerebral blood flow measured with 3D pseudocontinuous arterial spin-labeling MR imaging in Alzheimer disease and mild cognitive impairment: a marker for disease severity. Radiology 2013; 267: 221–230.

43.

Muller

van der Graaf

Visseren

, et al. Hypertension and longitudinal changes in cerebral blood flow: the SMART-MR study. Ann Neurol 2012; 71: 825–833.

44.

Suri

Topiwala

Chappell

, et al. Association of midlife cardiovascular risk profiles with cerebral perfusion at older ages. JAMA Netw Open 2019; 2: e195776.

45.

Wierenga

Clark

Dev

, et al. Interaction of age and APOE genotype on cerebral blood flow at rest. J Alzheimers Dis 2013; 34: 921–935.

46.

Wolf

D’Agostino

Belanger

, et al. Probability of stroke: a risk profile from the Framingham study. Stroke 1991; 22: 312–318.

47.

D’Agostino

Wolf

Belanger

, et al. Stroke risk profile: adjustment for antihypertensive medication. The Framingham study. Stroke 1994; 25: 40–43.

48.

Benjamini

Hochberg

Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc B Stat Methodol 1995; 57: 289–300.

49.

Zhang

Association of cerebral hypoperfusion and poor collaterals with cognitive impairment in patients with severe vertebrobasilar artery stenosis. J Alzheimers Dis Rep 2024; 8: 999–1007.

50.

Sheth

Liebeskind

DS.

Imaging evaluation of collaterals in the brain: physiology and clinical translation. Curr Radiol Rep 2014; 2: 29.

51.

Iqbal

A comprehensive study of the anatomical variations of the circle of Willis in adult human brains. J Clin Diagn Res 2013; 7: 2423–2427.

52.

Zaninovich

Ramey

Walter

, et al. Completion of the circle of Willis varies by gender, age, and indication for computed tomography angiography. World Neurosurg 2017; 106: 953–963.

53.

Eaton

Shah

Dornbos

, et al. Demographic age-related variation in circle of Willis completeness assessed by digital subtraction angiography. Brain Circ 2020; 6: 31–37.

54.

Chopra

Neuberger

Prince

, et al. Age-related changes in the completeness of the circle of Willis in children. Childs Nerv Syst 2022; 38: 1181–1184.

55.

Triantafyllou

Tudose

Tsiouris

, et al. The anterior communicating artery variants: a meta-analysis with a proposed classification system. Surg Radiol Anat 2024; 46: 697–716.

56.

Binder

El Amki

Glück

, et al. Leptomeningeal collaterals regulate reperfusion in ischemic stroke and rescue the brain from futile recanalization. Neuron 2024; 112: 1456–1472.e6.

57.

Hindenes

Håberg

Johnsen

, et al. Variations in the circle of Willis in a large population sample using 3D TOF angiography: the Tromsø study. PLoS One 2020; 15: e0241373.

58.

Wang

Wei

X-E

, et al. Examination of structural variations of the circle of Willis by 3D time-of-flight magnetic resonance angiography. Front Neurosci 2020; 14: 71.

59.

Cassot

Vergeur

Bossuet

, et al. Effects of anterior communicating artery diameter on cerebral hemodynamics in internal carotid artery disease: a model study. Circulation 1995; 92: 3122–3131.

60.

Forgo

Tarnoki

, et al. Are the variants of the circle of Willis determined by genetic or environmental factors? Results of a twin study and review of the literature. Twin Res Hum Genet 2018; 21: 384–393.

61.

Chen

X-Y

Soo

, et al. Benign oligemia in subacute stage is associated with borderzone infarction in stroke patients caused by intracranial large artery disease. Eur Neurol 2017; 77: 80–86.

Posterior communicating artery variation is linked to compromised cerebral hemodynamics in middle-aged and older adults

Abstract

Keywords

Introduction

Methods

Study cohort

Brain MRI—CoW variants

Brain MRI—CBF and CBF sCoV

Brain MRI—harmonization

Analytical plan

Results

Participant characteristics

CoW variant prevalence

CoW variants and CBF

CoW variants and CBF sCoV

Discussion

Footnotes

Acknowledgements

Non-standard abbreviations and acronyms

Author contributions

Funding

Declaration of conflicting interests

ORCID iD

References