Skip to main content

Longitudinal relaxographic imaging of white matter hyperintensities in the elderly



Incidental white matter hyperintensities (WMHs) are common findings on T2-weighted magnetic resonance images of the aged brain and have been associated with cognitive decline. While a variety of pathogenic mechanisms have been proposed, the origin of WMHs and the extent to which lesions in the deep and periventricular white matter reflect distinct etiologies remains unclear. Our aim was to quantify the fractional blood volume (vb) of small WMHs in vivo using a novel magnetic resonance imaging (MRI) approach and examine the contribution of blood–brain barrier disturbances to WMH formation in the deep and periventricular white matter.


Twenty-three elderly volunteers (aged 59–82 years) underwent 7 Tesla relaxographic imaging and fluid-attenuated inversion recovery (FLAIR) MRI. Maps of longitudinal relaxation rate constant (R1) were prepared before contrast reagent (CR) injection and throughout CR washout. Voxelwise estimates of vb were determined by fitting temporal changes in R1 values to a two-site model that incorporates the effects of transendothelial water exchange. Average vb values in deep and periventricular WMHs were determined after semi-automated segmentation of FLAIR images. Ventricular permeability was estimated from the change in CSF R1 values during CR washout.


In the absence of CR, the total water fraction in both deep and periventricular WMHs was increased compared to normal appearing white matter (NAWM). The vb of deep WMHs was 1.8 ± 0.6 mL/100 g and was significantly reduced compared to NAWM (2.4 ± 0.8 mL/100 g). In contrast, the vb of periventricular WMHs was unchanged compared to NAWM, decreased with ventricular volume and showed a positive association with ventricular permeability.


Hyperintensities in the deep WM appear to be driven by vascular compromise, while those in the periventricular WM are most likely the result of a compromised ependyma in which the small vessels remain relatively intact. These findings support varying contributions of blood–brain barrier and brain-CSF interface disturbances in the pathophysiology of deep and periventricular WMHs in the aged human brain.


White matter hyperintensities (WMHs) are commonly observed in the elderly brain on T2- weighted and proton density magnetic resonance imaging (MRI). The prevalence and severity of lesions increases with age and even mild increases in WMH volume may diminish cognition and global function in later life [1, 2]. The presence of WMHs has been associated with a number of disease states in the elderly, including stroke and dementia. In patients with Alzheimer’s disease, WMH severity is often related to deficits in attention and frontal executive functions, leading to the suggestion that WMHs, in combination with amyloid deposition, may provide a ‘second hit’ to the degenerating brain, increasing the risk and/or clinical expression of the disease [3].

The etiology of incidental WMHs in the elderly is poorly understood. While neuropathological studies are informative, the myelin disruption, gliosis, and edema commonly observed in the parenchyma tend to reflect multiple inflammatory, neurodegenerative, and ischemic processes. In addition, the extent of these features varies with lesion severity, further obscuring the underlying mechanisms of their formation [4]. By and large, though, WMHs are a prominent feature of small vessel disease. Moreover, the severity of lesions is often associated with hypertension, diabetes, and hypercholesterolemia, providing additional support for the role of blood–brain barrier disturbances in their pathogenesis [5, 6]. However, accumulating evidence suggests that the contribution of vascular factors to WMH formation may differ in the deep compared to the periventricular WM. In contrast to deep lesions (dWMHs), periventricular WMHs (pWMHs) are generally independent of infarct prevalence and plasma homocysteine levels, and show inconsistent association with hypertension and systolic blood pressure, especially in their early stages [711]. Moreover, histology regularly reveals defects in the adjacent ependymal lining suggesting that cerebrospinal fluid (CSF) infiltration may play a particularly important role in pWMH formation [12].

In recent years, quantitative MRI techniques have confirmed many of the histopathological features of WMHs in vivo, including myelin disruption (decreased magnetization transfer ratio) [13] and Wallerian degeneration (reduced fractional anisotropy) [14]. In addition, arterial spin labeling [15, 16] and dynamic susceptibility contrast (DSC) MRI [17] have demonstrated WMH perfusion deficits compared to normal appearing white matter (NAWM). To extend these findings, we used T1 relaxographic imaging to measure the fractional blood volume (vb) of small, focal WMHs in the deep WM (centrum semiovale) and in mild hyperintense periventricular linings and caps. Relaxographic imaging is based on the measurement of the longitudinal relaxation rate constant of water protons (R1 ≡ 1/T1) before and after addition of a paramagnetic contrast reagent (CR). In cell suspensions and excised tissue, the technique provides information about intra- and extracellular volume and water residence times [18]. In combination with spatial encoding, the approach provides a quantitative measure of water exchange dynamics at the brain barriers. Relaxographic imaging has now been applied to a variety of human pathologies, including osteosarcoma, breast carcinoma, and multiple sclerosis [1921].

Here, we measured R1 values throughout the brain before CR injection and during washout in a group of elderly adults with minimal cerebrovascular risk factors. Absolute values of vb were then obtained by fitting the time-course of tissue R1 changes to a pharmacokinetic model that incorporates the effects of water exchange between the blood and extravascular compartments. Next, we segmented hyperintense from normal-appearing WM on fluid-attenuated inversion recovery (FLAIR) images and quantified the vb in dWMHs and pWMHs on co-registered vb maps. Finally, to investigate the potential role of degenerative factors and permeability changes at the ventricular CSF-tissue interface on pWMH formation, we examined the association of pWMH vb with ventricular volume and with the rate of CR leakage into the CSF.



Twenty-three elderly individuals in good general health were enrolled. All participants underwent complete neurological and neuropsychological evaluation and provided a non-fasting blood sample for analysis of plasma homocysteine levels by fluorescence polarization immunoassay [22]. Subjects were excluded with body mass index greater than 30 kg/m2, uncontrolled hypertension, or history of diabetes mellitus, cardiovascular disease or TIA/stroke. Ten subjects were diagnosed with mild cognitive impairment (MCI) based on research criteria [23]. These individuals satisfied all criteria for study inclusion but had lower Mini-Mental State Examination (MMSE) [24] scores than cognitively normal subjects (26 ± 2 and 29 ± 1, respectively; P = .01). Three participants showed no radiographic evidence of WMHs, leaving a total of 20 datasets for analysis. Demographic features of these participants are shown in Table 1. The study was approved by the Oregon Health & Science University Institutional Review Board and all volunteers signed written informed consent.

Table 1 Subject characteristics (N = 20) a

MRI data acquisition

All images were acquired using a 7 Tesla Magnetom scanner (Siemens Healthcare, Erlangen, Germany) with 8- or 24-channel transmit/receive head coil. High-resolution inversion recovery magnetization prepared rapid acquisition gradient echo (MPRAGE) [inversion time (TI) = 1050 ms, repetition time (TR) = 2300 ms, echo time (TE) = 2.8 ms, flip angle (FA) = 6°, field of view (FOV) = 256 × 256 mm2, reconstruction matrix = 320 × 320 pixels, and slice thickness = 0.8 mm] and FLAIR, [TE = 386 ms, TI = 2150 ms, TR = 8000 ms, FOV = 224 × 256 mm2, matrix = 280 × 320, slice thickness = 0.8 mm] images were used for volumetric measurements and WMH segmentation. The nominal spatial resolution of both sequences was 0.8 mm3. Dual T2- and proton density-weighted turbo spin echo with TE = 11 and 87 ms, TR = 10000 ms, FOV = 168 × 224 mm2, matrix = 256 × 192 mm2, slice thickness = 2 mm with 2 mm interslice gap were also obtained.

Variable inversion time MPRAGE datasets were acquired for T1-relaxography. Each dataset consisted of four axial MPRAGE images with TI = 300, 1800, or 3200 ms, or with no inversion pulse. Other sequence parameters were: TE = 2.3 ms, TR = 3500 ms, FA = 6°, FOV = 176 × 256 mm2, matrix = 176 × 256 pixels, slice thickness = 2 mm, resulting in a nominal spatial resolution of 1 × 1 × 2 mm3. All images were centered on the lateral ventricles, oriented parallel to the anterior commissure-posterior commissure line. MPRAGE datasets were acquired prior to a 0.11 mmol/kg dose of intravenous CR (gadoteridol; Bracco Diagnostics, Princeton, NJ, USA) and throughout the subsequent 50–60 minute washout period. The acquisition time for each dataset was 14.5 min.

Image processing and analysis

An overview of the image analysis procedure is shown in Figure 1. Binary masks of NAWM and CSF were prepared using FMRIB’s Software Library (FSL; Technical details of FSL image analysis tools have been described and validated previously [25]. Briefly, high-resolution MPRAGE images underwent bias correction [26], brain extraction (with eye and optic nerve cleanup) [27] and 3-class segmentation [28]. Use of default segmentation parameters resulted in exclusion of WMHs (which appear relatively hypointense on T1-weighted images) from the WM segment. NAWM masks were prepared as a subset of all WM voxels within 15 (±3) mm of the superior surface of the lateral ventricles. Ventricular CSF masks were prepared by manual removal of parenchymal pixels from the CSF binary image. CSF volume was calculated using partial volume estimates and normalized to the total intracranial volume. Binary WMH masks were prepared after linear registration (6 degrees of freedom) [29] of skull-stripped FLAIR images to high resolution MPRAGE space followed by segmentation of WMHs based on manual seeding and k-nearest-neighboring clustering (Seg3D; Individual dWMH and pWMH masks were prepared by an experienced investigator. Careful attention was paid to avoid fluid-filled boundaries to minimize partial volume effects. Voxels that did not appear on at least three consecutive slices in all three orientations were not included. In 15 subjects, masks were independently prepared by a second trained investigator. Based on a 79 and 84% pixel overlap of dWMHs and pWMHs, respectively, we estimate segmentation reproducibility to be approximately 80%.

Figure 1
figure 1

Data workflow: (Left, top) FLAIR images were linearly registered to high-resolution MPRAGE space using rigid body registration. (Left, lower) Co-registered R1 maps were prepared in high-resolution MPRAGE space after alignment of all variable-TI MPRAGE images to the TI 1800 image collected at the temporal mid-point of CR washout (see Methods for additional details). (Right) Fractional blood volume (vb) maps of pWMHs and dWMHs and corresponding FLAIR image. Insets show a magnified view of the WMHs (outlined in white) with associated color scales. The vb histograms of each WMH are also shown.

The severity of WMHs was assessed using the Fazekas rating scale [30]. WMHs were defined as WM areas with increased signal intensity on T2-spin echo and FLAIR images. Special care was taken to distinguish WMHs from lacunar infarcts, defined as parenchymal defects not extending to the cortical gray matter with a signal intensity following that of the CSF on all pulse sequences [31]. WMHs were present in all subjects; 15 participants had lesions in both deep and periventricular WM. In 17/20 subjects, lesions were seen as small scattered hyperintensities or as caps and pencil-thin linings around the ventricles (≤grade 1). In the remaining subjects, WMHs (3 periventricular, 3 deep) were more numerous and included a few patchy or more diffuse lesions. These were rated as grade 2. The total volume of pWMHs and dWMHs was 3.0 (±2.9) and 0.55 (±0.49) mL, respectively.

Pharmacokinetic analysis

R1 maps were prepared from variable inversion time MPRAGE datasets. After brain extraction, all images in the dataset collected midway through CR washout (approximately 30 min post-CR) were linearly registered to the TI 1800 MPRAGE image. Each image in the remaining datasets was then registered to one of these images by application of the TI-dependent transform. R1 maps were prepared from each dataset by nonlinear least squares curve fitting (IDL; ITT Visual Information Solutions, Boulder CO, USA) of the signal intensity at each inversion time to the MPRAGE signal equation that includes magnetization losses due to sampling during recovery [32]. After preparation, all R1 maps were linearly registered to high-resolution MPRAGE space.

Full volume vb maps were prepared by voxelwise fitting of water proton R1 values before and at each time after CR administration to a two-site (blood, extravascular) tissue model that accounts for the effects of transendothelial water exchange. The model assumes that each voxel exhibits monoexponential recovery after inversion and that CR remains confined to the intravascular space (i.e., Ktrans, the rate constant for CR extravasation, is negligible). Since tissue water molecules exchange across the capillary walls, the measured R1 value of tissue is a function of the relaxation behavior of water protons in each compartment. Due to the short range of the interactions responsible for longitudinal relaxation quenching by CR (typically on the order of several nanometers), only the R1 of blood is increased by an intravascular CR; the R1 of water protons in the extravascular space is unaffected. Therefore, the observed relaxation behavior of tissue becomes a function of the R1 value in the blood and the rate at which water exchanges across the capillary wall [33]. Under these conditions, it has been shown that the R1 behavior of water protons in tissue during CR washout from the blood, R1t(t), can be described by Eq. (1) where R1b(t) is the longitudinal relaxation rate constant of intravascular water at time t, R1t(0) and R1b(0) are the R1 values in tissue and blood before CR administration, and τb−1 is the rate constant for water extravasation [33, 34]. fw is a scaling factor that accounts for the restricted localization of CR to the blood plasma. R1e is the rate constant for extravascular water relaxation in the absence of exchange and is given by Eq. (2).


Values of vb were iteratively determined using a gradient-expansion, non-linear least-squares algorithm [35] with τb and fw held constant at reasonable values (0.3 and 0.8 s, respectively) [36, 37]. R1b(t) values were measured from a region of interest contained entirely within the sagittal sinus. Uncertainties in R1 and vb values were estimated using a Monte Carlo approach. The signal-to-noise ratio (SNR) of each MPRAGE image was calculated using the signal-background method [38]. Noise was then added by randomly sampling a Gaussian distribution with mean = 0 and standard deviation = signal intensity/SNR. Several voxels in each variable-TI dataset corresponding to WMHs were randomly selected and R1t values estimated 1000 times. R1t values were found to be robust to noise, with the variance in all voxels <0.04%. Errors in vb were estimated by randomly selecting R1(t) values within 1 standard deviation of the mean and fitting to Eq. (1), as described above. Based on 1000 simulations, we estimate vb variance to be approximately 1-2% of the mean.

Ventricular permeability

Ventricular permeability was assessed based on region of interest (ROI) analysis of the temporal changes in CSF R1 values during CR washout. Unilateral ROIs (16 ± 9 μL) were defined in the posterior horns of the lateral ventricles on R1(t) maps. To minimize partial volume errors, ROIs were selected in areas that afforded optimal conspicuity of CSF-tissue interfaces. In the majority of cases, this was in the left posterior horn. Pixels with estimated fractional CSF volumes less than 1 were excluded to further reduce partial volume contamination.

Statistical analyses

Data analysis was performed in Stata (StataCorp, College Station, TX, USA). Distributional outliers were identified by Mahalanobis and Cook’s Distance methods [39]. The mean vb of WMH voxels was determined by fitting to a Gaussian distribution. Distributions that could not be normalized by common transformations (N = 3) were fit to non-negative, bounded (gamma or Weibull) distributions using EasyFit (Mathwave Technologies, Spokane, WA, USA). Voxels with nonphysiological vb values (i.e., less than 0 or greater than 20 mL/100 g) were excluded from analysis. On average, these voxels accounted for 29 and 34% of all dWMH and pWMH voxels, respectively. Univariate associations were examined using linear regression analysis and the coefficient of determination, r2. Regional differences were examined using one-way within subject ANOVA with WM type (pWMH, dWMH, or NAWM) as the independent variable. P values < .05 were considered significant in statistical tests. Post-hoc pairwise comparisons were tested using paired t-tests. A modified Bonferroni procedure was used to adjust the level for (two-tailed) significance testing [40].


In the absence of CR, the average (±standard deviation) longitudinal relaxation time (T1 ≡ R1−1) of NAWM was 1454 (±54) ms. T1 values depend on field strength and vary somewhat with sequence and post-processing details. Nevertheless, this value corresponds well with the 1220 (±36) ms reported by Rooney and colleagues [41] and the 1500 (±100) ms reported by Li and Deoni [42]. Compared to NAWM, T1 times were increased in both pWMHs and dWMHs (1861 ± 124 ms and 1608 ± 80 ms, respectively; P < .0001). To our knowledge, T1 values of WMHs at 7Tesla have not been reported previously. Nevertheless, the 10-28% increase in WMH T1 values compared to NAWM measured here is consistent with the 5-40% increase in WMH T1 times noted at lower field strengths [43, 44]. T1 values were significantly larger in pWMHs than dWMHs (P < .0001), although the between-subject variance was also increased in the former, most likely due to increased partial volume effects. Assuming a linear dependence of R1 values on tissue water fraction [45] and that T1 differences are due solely to alterations in water content, this corresponds to roughly a 15% increase in the total water fraction of pWMHs compared to dWMHs. Associations of WMH T1 values with age, gender, memory complaint (present or absent) and homocysteine levels were not significant.

As expected, R1t values increased after CR administration in all subjects. Figure 2 shows the R1t behavior in a representative dWMH and pWMH voxel plotted against the R1b values. The best fit of the data in each voxel to Eq. (1) is also shown. The concentration of CR in plasma ([CRp]) depends on R1b[19] as:

Figure 2
figure 2

Relationship between R 1t and R 1b values in a dWMH (circle) and pWMH (square) voxel (shown at far right). Curves represent the best fit of the data in each WMH to Eq. (1). Error bars (±SD) are based on Monte Carlo simulations (see text). Upper axis represents plasma CR concentration, as calculated from Eq. (3).


where h is the microvascular hematocrit (ca. 0.40) [46] and r1 is the relaxivity of CR at 7 T (3.3 s−1 mM−1 for gadoteridol) [47] and is plotted along the upper axis. The nonlinear relationship between R1t and [CRp] observed here (especially prominent in pWMHs) is consistent with the compartmental nature of tissue and the restricted exchange of water between blood and extravascular sites.

The mean vb values in NAWM averaged 2.4 (±0.6) mL/100 g, in good agreement with the 2.0 (±0.4) mL/100 g measured by DSC MRI [48] and the 2.7 (±0.4) mL/100 g found by two-point T1 relaxography in healthy adults [49]. In dWMHs and pWMHs, vb averaged 1.8 (±0.6) and 2.4 (±0.8) mL/100 g, respectively. Mean vb values varied significantly between dWMHs, pWMHs and NAWM (F2,14 = 6.33, P = .0054). Post-hoc tests showed a significant reduction in the vb of dWMHs compared to NAWM (P < .0001). After correcting for multiple comparisons, differences between pWMHs and dWMHs were not significant (P = .033). Similarly, no significant associations of vb with age, gender, or memory complaint were observed in periventricular or deep WMHs.

In the absence of CR, CSF R1 values averaged 0.239 (±0.023) s−1, in good agreement with the 0.231- 0.237 s−1 reported previously [41, 50]. These values increased to 0.244 (± 0.022) s−1 within 13 (±5) minutes after CR administration and continued to increase with time (F3,18 = 8.71, P < .0001). Figure 3 shows the R1 histograms of ventricular CSF during CR washout. The total change in R1 values over the 50–55 minute measurement period was 0.0091 (±0.01) s−1. Assuming a linear dependence of CSF R1 values on CR concentration, this corresponds to a ventricular permeability of approximately 3.4 μM hr−1. Although not statistically significant (P = 0.09), pWMH vb showed a moderate association (r2 = 0.22) with the R1 increase observed midway through the washout period (Figure 4a). An inverse association (r2 = 0.31; P = 0.04) of pWMH vb with ventricular volume was also observed (Figure 4b).

Figure 3
figure 3

Histograms of average CSF R 1 values during CR washout. The mean (dashed line) and fitted Gaussian distribution (solid line) are also shown. Time, t, was calculated as the temporal midpoint of each variable-TI image set. Data from one subject who completed only two post-CR measurements were excluded.

Figure 4
figure 4

Linear regression plots of pWMH v b and (a) percent change in CSF R 1 values 33 (±4) minutes after CR administration, and (b) ventricular volume, as a percentage of the total intracranial volume (ICV). Error bars represent the SEM.


Our results demonstrate that the fractional blood volume of small WMHs in the deep, but not periventricular WM, is decreased compared to NAWM in elderly individuals with minimal cardiovascular risk factors. While reduced vb in the NAWM of individuals with and without WMHs has been demonstrated by 15O-positron emission tomography [51] and DSC MRI [48], the relative vb of WMHs and NAWM has been examined in only a few studies. In migraine patients with larger (3–21 mm diameter) though still focal dWMHs, the vb of lesions was found by DSC MRI to be reduced compared to NAWM, consistent with our results [44]. A similar finding has been reported in stroke subjects, although the diffuse nature of lesions in these subjects combined with the relatively poor spatial resolution makes comparisons more ambiguous [52]. However using the same technique in a group of healthy elderly individuals all between the ages of 85 and 86 years, Marstrand and colleagues were unable to find any significant reduction in the vb of WMHs compared to NAWM [17]. Whether such differences are truly absent in the very elderly brain or were obscured by the corrections to their T2*-based acquisition parameters necessary to obtain absolute vb measurements or the more limited spatial resolution of their images and subsequent partial volume effects, remains to be seen.

The absolute reduction in the vb of dWMHs compared to pWMHs approached 25% in the present study. While confidence limits were too wide to establish the statistical significance of this finding (due at least in part because of small subject numbers), the magnitude of the difference suggests that alterations in blood–brain barrier permeability and/or decreased capillary density may be an important feature of lesion pathogenesis in the deep WM. In contrast, changes in capillary permeability appear to contribute relatively little to pWMH formation. Indeed, the microvascular defects and reduced capillary density frequently observed by histology in small dWMHs are generally not seen in early pWMHs [53]. Instead, hyperintense caps and pencil-thin periventricular linings are most commonly seen in association with a disrupted ependyma, presumably driven by abnormal osmotic pressure gradients or as a compensatory mechanism for the atrophic changes and ventricular enlargement that accompany normal aging [54]. Our findings of a significant inverse association between pWMH vb and ventricular volume combined with the association of CSF [CR] and dilation of pWMH vessels are in accord with these findings. It is also possible that age-related or ischemic changes in the blood-CSF barrier may contribute to formation of pWMHs. Changes in the blood-CSF barrier are known to occur with age and cerebrovascular compromise and could explain the increased R1 values in the CSF [55]. However, since paracellular diffusion of relatively large molecular weight molecules (gadoteridol MW >500 Da) is likely to be minimal and transcellular mechanisms of CR transport across the blood CSF barrier unknown, additional experiments, ideally with increased temporal resolution in both the early and long times after CR administration, will be necessary to examine the role of the blood CSF barrier in pWMH formation more closely.

To our knowledge, the present study is the first to determine the absolute vb of WMHs using relaxographic imaging. Our use of a T1- rather than T2*- based approach minimizes the systematic errors caused by bulk magnetic susceptibility effects which can degrade vb estimates based on T2*-weighted images [56]. However, T1-based approaches typically suffer from low SNR and reduced detectability of CR. To increase the sensitivity of our measurements, all imaging was done at 7 T, where the SNR in WM is ca. 50% higher than at 4 T and the contrast to noise ratio is substantially increased [41, 57]. The former results in improved detection of small lesions while the latter provides increased CR detection sensitivity and accuracy of R1 measurements. To further improve the accuracy of vb estimates, R1 changes were fit to a pharmacokinetic model that incorporates the effects of transendothelial water exchange on the MR signal. Failure to account for exchange effects can result in large errors in vb estimates, particularly in the early stages of lesion formation when vascular disturbances are expected to be small [58]. The model also assumes negligible leakage of CR from the vessels and a fixed rate of water exchange across the capillary walls. Such assumptions are reasonable in the healthy aged WM, where vascular changes are expected to be subtle. The extent to which they remain equally valid in both deep and pWMHs, though, remains unclear. Bolus tracking (dynamic contrast enhanced) MRI studies with increased temporal resolution are currently underway to examine the validity of these assumptions in both more closely.

Our results provide in vivo evidence that permeability changes in the blood–brain barrier are likely to accompany formation of dWMHs, while extravascular changes at the CSF-tissue interface may play a more important role in pWMH formation. However, our study is not without limitations. First, to minimize secondary effects of chronic ischemia, we excluded individuals with major cardiovascular risk factors. However, cognitive impairment was not an exclusion criterion and, although mild, was present in half of the study subjects. While the vb of neither the NAWM nor WMHs was significantly different in these subjects, we did not measure cerebral blood flow. As a result, it is difficult to know to what extent our conclusions are affected by the hypoperfusion and altered cerebral autoregulation frequently associated with cognitive decline [59]. Second, the pharmacokinetic model upon which our estimates are based assumes that CR remains in the blood throughout the measurement period. If leakage into the extravascular space becomes significant, errors in vb estimates can be substantial [36]. While CR extravasation across a relatively intact blood–brain barrier is likely to be minimal in the deep WM, a similar assumption in periventricular regions, where CR in the CSF may exchange with interstitial fluid at a disrupted ependymal lining, is less clear. As a result, the pathobiological basis of vb associations in these regions should be interpreted with caution. Third, the validity of our findings depends on the accurate segmentation of WMHs. While the FLAIR images acquired at high magnetic field strength provide excellent conspicuity of small WM lesions, quantitative segmentation of WMHs remains challenging. Semiautomatic methods, like the one used here, rely on placement of a ‘seed point’ and propagation of the seeds to neighboring slices. The approach eliminates the need for manual outlining and reduces rater bias and error [60]. Nevertheless, results can be strongly dependent on the exact location of the seed point, as evidenced by an inter-rater agreement of selected pixels of only 80%. To reduce false positives, only those pixels identified as WMHs by both raters were included in our analyses. While this approach likely decreases the specificity of segmentation results, it provides an important increase in their sensitivity and the validity of WMH measurements. Fourth, the accuracy of R1 values determined from signal intensity changes in variable-TI datasets, and upon which vb estimates rely, is made difficult by the complexity of the MPRAGE sequence [32]. To minimize these errors and improve the accuracy of R1 values, we used a fitting routine that models all the RF pulses, recovery periods and the approach to steady state magnetization assuming a constant flip angle across all imaging planes. However, we did not correct for RF inhomogeneities in the images which can be substantial at high field, particularly in central brain regions [57]. While the non-selective nature of the RF pulses and 3D image encoding make the MPRAGE sequence relatively insensitive to these inhomogeneities, they remain a potential source of R1 inaccuracies in our data. Finally, 3D relaxography, like all quantitative MRI, requires a balance between spatial and temporal resolution. To achieve the spatial resolution sufficient for identification of small WMHs, each variable-TI MPRAGE dataset required a scan times of almost 15 minutes. In agreement with measurements at 3 Tesla [61], the accuracy of the R1 values obtained by this method is quite good; however, the precision of vb estimates would likely be increased with additional sampling frequency.


Our findings suggest that in the absence of major vascular risk factors, the fractional blood water in WMHs is increased compared to NAWM in the healthy aged brain. While blood–brain barrier disturbances appear to be a prominent feature of WMHs in the deep WM, disturbances in the brain-CSF interface may play a more important role in the etiology of pWMHs.

Authors’ information

VCA- PhD (Chemistry), MCR (Clinical research)


JAK- MD (Neurology)

JFQ- MD (Neurology)

PB- MD (Neurosurgery)

LPR- MD (Radiology)


WDR- PhD (Chemistry)



Contrast reagent


Contrast reagent concentration in plasma


Cerebrospinal fluid


Deep white matter hyperintensity


Dynamic susceptibility contrast


Fluid-attenuated inversion recovery


Intracranial volume


Magnetic resonance imaging


Magnetization prepared rapid acquisition gradient echo


Normal appearing white matter


Periventricular white matter hyperintensity


Blood longitudinal relaxation rate constant


Tissue longitudinal relaxation rate constant


Region of interest


Signal-to-noise ratio


Longitudinal relaxation time


Fractional blood volume.


  1. Murray M, Senjem M, Petersen R, Hollman J, Preboske G, Weignad S, Knopman D, Ferman T, Dickson D, Jack C: Functional impact of white matter hyperintensities in cognitively normal elderly subjects. Arch Neurol. 2010, 67: 1379-1385.

    Article  PubMed Central  PubMed  Google Scholar 

  2. Garde E, Mortensen L, Krabbe K, Rostrup E, Larsson H: Relation between age-related decline in intelligence and cerebral white-matter hyperintensities in healthy octogenarians: a longitudinal study. Lancet. 2000, 356: 628-634. 10.1016/S0140-6736(00)02604-0.

    Article  CAS  PubMed  Google Scholar 

  3. Provenzano F, Muraskin J, Tosto G, Narkhede A, Wasserman B, Griffith E, Guzman V, Meier I, Zimmerman M, Brickman A: White matter hyperintensities and cerebral amyloidosis. JAMA Neurol. 2013, 70: 455-461. 10.1001/jamaneurol.2013.1321.

    Article  PubMed Central  PubMed  Google Scholar 

  4. Bronge L, Wahlund L: White matter changes in dementia: does radiology matter?. Br J Radiol. 2007, 80: S115-S120. 10.1259/bjr/35265137.

    Article  PubMed  Google Scholar 

  5. Erten-Lyons D, Woltjer R, Kaye J, Mattek N, Dodge H, Green S, Tran H, Howieson D, Wild K, Silbert L: Neuropathologic basis of white matter hyperintensity accumulation with advanced age. Neurology. 2013, 81: 977-983. 10.1212/WNL.0b013e3182a43e45.

    Article  PubMed Central  PubMed  Google Scholar 

  6. Yamada S, Saiki M, Satow T, Fukuda A, Ito M, Minami S, Miyamoto S: Periventricular and deep white matter leukoaraiosis have a closer association with cerebral microbleeds than age. Eur J Neurol. 2012, 19: 98-104. 10.1111/j.1468-1331.2011.03451.x.

    Article  CAS  PubMed  Google Scholar 

  7. Sachdev P, Parslow R, Salonikas C, Lux O, Wen W, Maller J, Kumar R, Naidoo D, Christensen H, Jorm A: Homocysteine and the brain in midadult life: evidence for an increased risk of leukoaraiosis in men. Arch Neurol. 2004, 61: 1369-1376. 10.1001/archneur.61.9.1369.

    Article  PubMed  Google Scholar 

  8. Murray A, Staff R, Shenkin S, Deary I, Starr J, Whalley L: Brain white matter hyperintensities: relative importance of vascular risk factors in nondemented elderly people. Radiology. 2005, 237: 251-257. 10.1148/radiol.2371041496.

    Article  PubMed  Google Scholar 

  9. Wen W, Sachdev P, Shnier R, Brodaty H: Effect of white matter hyperintensities on cortical cerebral blood volume using perfusion MRI. Neuroimage. 2004, 21: 1350-1356. 10.1016/j.neuroimage.2003.11.015.

    Article  PubMed  Google Scholar 

  10. Spilt A, Geeraedts T, de Craen A, Westendorp R, Blauw G, van Buchem M: Age-related changes in normal-appearing brain tissue and white matter hyperintensities: more of the same or something else?. Am J Neuroradiol. 2005, 26: 725-729.

    PubMed  Google Scholar 

  11. Gunstad J, Cohen R, Tate D, Paul R, Poppas A, Hoth K, Macgregor K, Jefferson A: Blood pressure variability and white matter hyperintensities in older adults with cardiovascular disease. Blood Press. 2005, 14: 353-358. 10.1080/08037050500364117.

    Article  PubMed Central  PubMed  Google Scholar 

  12. Fazekas F, Kleinert R, Offenbacher H, Schmidt R, Kleinert G, Payer F, Radner H, Lechner H: Pathologic correlates of incidental MRI white matter signal hyperintensities. Neurology. 1993, 43: 1683-1689. 10.1212/WNL.43.9.1683.

    Article  CAS  PubMed  Google Scholar 

  13. Fazekas F, Ropele S, Enzinger C, Gorani F, Seewann A, Petrovic K, Schmidt R: MTI of white matter hyperintensities. Brain. 2005, 128: 2926-2932. 10.1093/brain/awh567.

    Article  PubMed  Google Scholar 

  14. Lee D, Fletcher E, Martinez O, Ortega M, Zozulya N, Kim J, Tran J, Buonocore M, Carmichael O, DeCarli C: Regional pattern of white matter microstructural changes in normal aging, MCI, and AD. Neurology. 2009, 73: 1722-1728. 10.1212/WNL.0b013e3181c33afb.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  15. Brickman A, Zahra A, Muraskin J, Steffener J, Holland C, Habeck C, Borogovac A, Ramos M, Brown T, Aslanni I, Stern Y: Reduction in cerebral blood flow in areas appearing as white matter hyperintensities on magnetic resonance imaging. Psychiatry Res Neuroimaging. 2009, 172: 117-120. 10.1016/j.pscychresns.2008.11.006.

    Article  PubMed  Google Scholar 

  16. Bastos-Leite A, Kuijer J, Romboutsd S, Sanz-Arigitad E, Van Straaten E, Gouw A, Van der Flier W, Scheltens P, Barkhof F: Cerebral blood flow by using pulsed arterial spin-labeling in elderly subjects with white matter hyperintensities. AJNR Am J Neuroadiol. 2008, 29: 1296-1301. 10.3174/ajnr.A1091.

    Article  CAS  Google Scholar 

  17. Marstrand J, Garde E, Rostrup E, Ring P, Rosenbaum S, Mortensen E, Larsson H: Cerebral perfusion and cerebrovascular reactivity are reduced in white matter hyperintensities. Stroke. 2002, 33: 972-976. 10.1161/01.STR.0000012808.81667.4B.

    Article  CAS  PubMed  Google Scholar 

  18. Labadie C, Lee J-H, Vetek G, Springer C: Relaxographic imaging. J Magn Reson B. 1994, 105: 99-112. 10.1006/jmrb.1994.1109.

    Article  CAS  PubMed  Google Scholar 

  19. Yankeelov T, Rooney W, Li X, Springer C: Variation of the relaxographic “shutter-speed” for transcytolemmal water exchange affects the CR bolus-tracking curve shape. Magn Reson Med. 2003, 50: 1151-1169. 10.1002/mrm.10624.

    Article  PubMed  Google Scholar 

  20. Labadie C, Lee J, Rooney W, Jarchow S, Aubert-Frecon M, Springer CJ, Moller H: Myelin water mapping by spatially regularized longitudinal relaxographic imaging at high magnetic fields. Magn Reson Med. 2014, 71: 375-387. 10.1002/mrm.24670.

    Article  PubMed  Google Scholar 

  21. Huang W, Li X, Morris E, Tudorica L, Seshan V, Rooney W, Tagge I, Wang Y, Springer C: The magnetic resonance shutter speed discriminates vascular properties of malignant and benign breast tumors in vivo. Proc Nat Acad Sci. 2008, 105: 17943-17948. 10.1073/pnas.0711226105.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  22. Hooshmand B, Polvikoski T, Kivipelto M, Tanskanen M, Myllykangas L, Erkinjuntti T, Oinas M, Paetau A, Scheltens P, Van Straaten E, Sulkava R, Solomon A: Plasma homocysteine, Alzheimer and cerebrovascular pathology: a population-based autopsy study. Brain. 2013, 136: 2707-2716. 10.1093/brain/awt206.

    Article  PubMed Central  PubMed  Google Scholar 

  23. McKhann G, Drachman D, Folstein M, Katzman R, Price D, Stadlan E: Clinical diagnosis of Alzheimer’s disease: report of the NINCDS-ADRDA Work Group under the auspices of the Department of Health and Human Services Task Force on Alzheimer’s Disease. Neurology. 1984, 34: 939-944. 10.1212/WNL.34.7.939.

    Article  CAS  PubMed  Google Scholar 

  24. Folstein M, Folstein S, McHugh P: “Mini-mental state”: a practical method for grading the cognitive state of patients for the clinician. J Psychiatr Res. 1975, 12: 189-198. 10.1016/0022-3956(75)90026-6.

    Article  CAS  PubMed  Google Scholar 

  25. Smith S, Jenkinson M, Woolrich M, Beckmann C, Behrens T, Johansen-Berg H, Bannister P, De Luca M, Drobnjak I, Flitney D, Niazy RK, Saunders J, Vickers J, Zhang Y, De Stefano N, Brady JM, Matthews PM: Advances in functional and structural MR image analysis and implementation as FSL. Neuroimage. 2004, 23: S208-S219.

    Article  PubMed  Google Scholar 

  26. Berlow Y, Sammi M, Selzer A, Bourdette D, Rooney W: Quantitative volumetrics of multiple sclerosis brain from 7T MRI. Proc Int Soc Mag Reson Med. 2013, 21: 3603-

    Google Scholar 

  27. Smith S: Fast robust automated brain extraction. Hum Brain Mapp. 2002, 17: 143-155. 10.1002/hbm.10062.

    Article  PubMed  Google Scholar 

  28. Patenaude B, Smith S, Kennedy D, Jenkinson M: A Bayesian model of shape and appearance for subcortical brain. Neuroimage. 2011, 56: 907-922. 10.1016/j.neuroimage.2011.02.046.

    Article  PubMed Central  PubMed  Google Scholar 

  29. Jenkinson M, Bannister P, Brady M, Smith S: Improved optimization for the robust and accurate linear registration and motion correction of brain images. Neuroimage. 2002, 17: 824-841.

    Article  Google Scholar 

  30. Fazekas F, Chawluk J, Alavi A, Hurtig H, Zimmerman R: MR signal abnormalities at 1.5 T in Alzherimer’s dementia and normal aging. AJR Am J Neuroradiol. 1987, 149: 351-356.

    CAS  Google Scholar 

  31. Van den Boom R, Oberstein S, Spilt A, Behloul F, Ferrari M, Haan J, Westendorp R, van Buchem M: Cerebral hemodynamics and white matter hyperintensities in CADASIL. J Cereb Blood Flow Metab. 2003, 23: 599-604.

    Article  PubMed  Google Scholar 

  32. Wright E, Mougin O, Totman J, Peters A, Brookes M, Coxon R, Morris P, Clemence M, Francis S, Bowtell R, Gowland P: Water proton T1 measurements in brain tissue at 7, 3, and 1.5 T using IR-EPI, IR-TSE, and MPRAGE: results and optimization. Magn Reson Med. 2008, 21: 121-130.

    CAS  Google Scholar 

  33. Schwarzbauer C, Morrissey S, Deichmann R, Hillenbrand C, Syha J, Adolf H, Noth U, Haase A: Quantitative magnetic resonance imaging of capillary water permeability and regional blood volume with an intravascular MR contrast agent. Magn Reson Med. 1997, 37: 769-777. 10.1002/mrm.1910370521.

    Article  CAS  PubMed  Google Scholar 

  34. Njus J, Vigeland L, Li X, Springer C, Taylor M, Telang F, Coyle P, Rooney W: Quantitative DCE 1H2O R1 measurements suggest increased fractional blood water in MS normal appearing brain tissue. Proc Int Soc Mag Reson Med. 2007, 15: 2193-

    Google Scholar 

  35. Bevington P, Robinson D: Data Reduction and Error Analysis for the Physical Sciences. 2003, New York: McGraw-Hill

    Google Scholar 

  36. Rooney W, Yankeelov T, Coyle P, Telang F, Springer C: Regional blood volumes and intravascular water lifetimes in human brain. Proc Int Soc Magn Reson Med. 2003, 11: 2188-

    Google Scholar 

  37. Landis C, Li X, Telang F, Coderre J, Micca P, Rooney W, Latour L, Vétek G, Pályka I, Springer C: Determination of the MRI contrast agent concentration time course in vivo following bolus injection: effect of equilibrium transcytolemmal water exchange. Magn Reson Med. 2000, 44: 563-574. 10.1002/1522-2594(200010)44:4<563::AID-MRM10>3.0.CO;2-#.

    Article  CAS  PubMed  Google Scholar 

  38. McRobbie D, Moore E, Graves M, Prince M: MRI: From Picture to Proton. 2007, NY: Cambridge University Press

    Google Scholar 

  39. Hodge V, Austin J: A survey of outlier detection methodologies. Artificial Intelligence Rev. 2004, 22: 85-126.

    Article  Google Scholar 

  40. Olejnik S, Li J, Supattathum S, Huberty C: Multiple testing and statistical power with modified Bonferroni procedures. J Educ Behav Stat. 1997, 22: 389-406. 10.2307/1165229.

    Article  Google Scholar 

  41. Rooney W, Johnson G, Li X, Cohen E, Kim S, Ugurbil K, Springer C: Magnetic field and tissue dependences of human brain longitudinal 1H2O relaxation in vivo. Magn Reson Med. 2007, 57: 308-318. 10.1002/mrm.21122.

    Article  CAS  PubMed  Google Scholar 

  42. Li T, Deoni C: Fast T1 mapping of the brain at 7 T with RF calibration using three point DESPOT1 method. Proc Int Soc Mag Reson Med. 2006, 14: 2643-

    Google Scholar 

  43. Gouw A, Seewann A, Vrenken H, Van der Flier W, Rozemuller J, Barkhof F, Scheltens P, Geurts J: Heterogeneity of white matter hyperintensities in Alzheimer’s disease: post-mortem quantitative MRI and neuropathology. Brain. 2008, 131: 3286-3298. 10.1093/brain/awn265.

    Article  CAS  PubMed  Google Scholar 

  44. Aradi M, Schwarcz A, Perlaki G, Orsi G, Kovacs N, Trauninger A, Kamson D, Erdelyi-Botor S, Nagy F, Nagy S, Doczi T, Komoly S, Pfund Z: Quantitative MRI studies of chronic brain white matter hyperintensities in migraine patients. Headache. 2013, 53: 752-763. 10.1111/head.12013.

    Article  PubMed  Google Scholar 

  45. Naruse S, Horikawa Y, Tanaka C, Hirakawa K, Nishikawa H, Yoshizaki K: Significance of proton relaxation time measurement in brain edema, cerebral infarction and brain tumors. Magn Reson Imag. 1986, 4: 293-304. 10.1016/0730-725X(86)91039-8.

    Article  CAS  Google Scholar 

  46. Li X, Springer CJ, Jerosch-Herold M: First-pass dynamic contrast-enhanced MRI with extravasating contrast reagent: evidence for human myocardial capillary recruitment in adenosine-induced hyperemia. NMR Biomed. 2009, 22: 148-157. 10.1002/nbm.1293.

    Article  PubMed  Google Scholar 

  47. Noebauer-Huhmann I, Kraff O, Juras V, Szomolanyi P, Maderwald S, Mlynarik V, Thyesohn J, Ladd S, Ladd M, Trattnig S: MR contrast media at 7T- preliminary study on relaxivities. Proc Int Soc Mag Reson Med. 2008, 16: 1457-

    Google Scholar 

  48. Markus H, Lythgoe D, Ostegaard L, O’Sullivan M, Williams S: Reduced cerebral blood flow in white matter in ischaemic leukoaraiosis demonstrated using quantitative exogenous contrast based perfusion MRI. J Neurol Neurosurg Psychiatry. 2000, 69: 48-53. 10.1136/jnnp.69.1.48.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  49. Shin W, Cashen T, Horowitz S, Sawlani R, Carroll T: Quantitative CBV measurement from static T1 changes in tissue and correction for intravascular water exchange. Magn Reson Med. 2006, 56: 138-145. 10.1002/mrm.20937.

    Article  PubMed  Google Scholar 

  50. Hopkins A, Yeung H, Bratton C: Multiple field strength in vivo T1 and T2 for cerebrospinal fluid protons. Magn Reson Med. 1986, 3: 303-311. 10.1002/mrm.1910030214.

    Article  CAS  PubMed  Google Scholar 

  51. Hatazawa J, Shimosegawa E, Satoh T, Toyoshima H, Okudera T: Subcortical hypoperfusion associated with asymptomatic white matter lesions on magnetic resonance imaging. Stroke. 1997, 28: 1944-1947. 10.1161/01.STR.28.10.1944.

    Article  CAS  PubMed  Google Scholar 

  52. Sachdev P, Wen W, Shnier R, Brodaty H: Cerebral blood volume in T2-weighted white matter hyperintensities using exogenous contrast based perfusion MRI. J Neuropsychiatry Clin Neurosci. 2004, 16: 83-92. 10.1176/appi.neuropsych.16.1.83.

    Article  PubMed  Google Scholar 

  53. Fazekas F, Kleinert R, Offenbacher H, Payer F, Schmidt R, Kleinert G, Radner H, Lechner H: The morphologic correlate of incidental punctate white matter hyperintensities on MR images. AJNR Am J Neuroradiol. 1991, 12: 915-921.

    CAS  PubMed  Google Scholar 

  54. Nestor S, Rupsingh R, Borrie M, Smith M, Accomazzi V, Wells J, Fogarty J, Bartha R, Initiative AsDN: Ventricular enlargement as a possible measure of Alzheimer’s disease progression validated using the Alzheimer’s disease neuroimaging initiative database. Brain. 2008, 131: 2443-2454. 10.1093/brain/awn146.

    Article  PubMed Central  PubMed  Google Scholar 

  55. Serot J, Bene M, Faure G: Choroid plexus, aging of the brain, and Alzheimer’s disease. Front Biosci. 2003, 8: s515-s521. 10.2741/1085.

    Article  PubMed  Google Scholar 

  56. Jochimsen T, Newbold R, Skare S, Clayton D, Albers G, Moseley M, Bammer R: Identifying systematic errors in quantitative dynamic-susceptibility contrast perfusion imaging by high-resolution multi-echo parallel EPI. NMR Biomed. 2007, 20: 429-438. 10.1002/nbm.1107.

    Article  PubMed Central  PubMed  Google Scholar 

  57. Vaughan J, Garwood M, Collins C, Liu W, DelaBarre L, Adriany G, Andersen P, Merkle H, Goebel R, Smith M, Ugurbil K: 7T vs. 4T: RF power, homogeneity, and signal-to-noise comparison in head images. Magn Reson Med. 2001, 46: 24-30. 10.1002/mrm.1156.

    Article  CAS  PubMed  Google Scholar 

  58. Li X, Rooney W, Springer CJ: A unified MRI pharmacokinetic theory: intravascular and extracellular contrast reagents. Magn Reson Med. 2005, 54: 1351-1359. 10.1002/mrm.20684.

    Article  CAS  PubMed  Google Scholar 

  59. Marshall R: Effects of altered cerebral hemodynamics on cognitive function. J Alzheimers Dis. 2012, 32: 633-642.

    PubMed  Google Scholar 

  60. Iorio M, Spalletta G, Chiapponi C, Luccichenti G, Caddiari C, Orfei M, Caltagirone C, Piras F: White matter hyperintensities segmentation: a new semi-automated method. Front Aging Neurosci. 2013, 5: 76-

    Article  PubMed Central  PubMed  Google Scholar 

  61. Liu J, Bock N, Silva A: Rapid high-resolution three-dimensional mapping of T1 and age-dependent variations in the non-human primate brain using magnetization-prepared rapid gradient-echo (MPRAGE) sequence. Neuroimage. 2011, 56: 1154-1163. 10.1016/j.neuroimage.2011.02.075.

    Article  PubMed Central  PubMed  Google Scholar 

Download references


This study was supported by grants from the National Institute on Aging (AG033638 and AG008017), the National Center for Research Resources (RR024140-04S1), and the Oregon Partnership for Alzheimer’s Research.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Valerie C Anderson.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

VCA, JAK, JFQ and WDR designed the study. JTO developed the software for and prepared the R1 and vb maps. LPR, PB, and DP assessed hyperintensity severity and performed lesion segmentations. VCA analyzed and interpreted the data, performed statistical analysis and drafted the manuscript. JAK, JFQ, and WDR revised the manuscript. All authors read and approved the final manuscript.

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Anderson, V.C., Obayashi, J.T., Kaye, J.A. et al. Longitudinal relaxographic imaging of white matter hyperintensities in the elderly. Fluids Barriers CNS 11, 24 (2014).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Aging
  • Blood–brain barrier
  • Blood volume
  • Relaxographic imaging
  • Periventricular
  • White matter hyperintensity
  • 7T