Application of Diffusion Tensor Imaging Parameters to Detect Change in Longitudinal Studies in Cerebral Small Vessel Disease

Eva Anna Zeestraten; Philip Benjamin; Christian Lambert; Andrew John Lawrence; Owen Alan Williams; Robin Guy Morris; Thomas Richard Barrick; Hugh Stephen Markus

doi:10.1371/journal.pone.0147836

Abstract

Cerebral small vessel disease (SVD) is the major cause of vascular cognitive impairment, resulting in significant disability and reduced quality of life. Cognitive tests have been shown to be insensitive to change in longitudinal studies and, therefore, sensitive surrogate markers are needed to monitor disease progression and assess treatment effects in clinical trials. Diffusion tensor imaging (DTI) is thought to offer great potential in this regard. Sensitivity of the various parameters that can be derived from DTI is however unknown. We aimed to evaluate the differential sensitivity of DTI markers to detect SVD progression, and to estimate sample sizes required to assess therapeutic interventions aimed at halting decline based on DTI data. We investigated 99 patients with symptomatic SVD, defined as clinical lacunar syndrome with MRI confirmation of a corresponding infarct as well as confluent white matter hyperintensities over a 3 year follow-up period. We evaluated change in DTI histogram parameters using linear mixed effect models and calculated sample size estimates. Over a three-year follow-up period we observed a decline in fractional anisotropy and increase in diffusivity in white matter tissue and most parameters changed significantly. Mean diffusivity peak height was the most sensitive marker for SVD progression as it had the smallest sample size estimate. This suggests disease progression can be monitored sensitively using DTI histogram analysis and confirms DTI’s potential as surrogate marker for SVD.

Citation: Zeestraten EA, Benjamin P, Lambert C, Lawrence AJ, Williams OA, Morris RG, et al. (2016) Application of Diffusion Tensor Imaging Parameters to Detect Change in Longitudinal Studies in Cerebral Small Vessel Disease. PLoS ONE 11(1): e0147836. https://doi.org/10.1371/journal.pone.0147836

Editor: Johannes Boltze, Fraunhofer Research Institution of Marine Biotechnology, GERMANY

Received: June 8, 2015; Accepted: January 8, 2016; Published: January 25, 2016

Copyright: © 2016 Zeestraten et al. 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 author and source are credited.

Data Availability: All relevant data are within the paper and its Supporting Information files.

Funding: The SCANS study was supported by a Wellcome Trust grant #081589, awarded to HSM (http://www.wellcome.ac.uk). HSM is supported by an English National Institute for Health Research Senior Investigator award (http://www.nihr.ac.uk). CL is supported by an Academy of Medical Sciences (GB) Clinical Lecturer Start-Up Grant (http://www.acmedsci.ac.uk). AJL is supported by a research grant from Alzheimer's Research UK: ARUK-PG2013-2 (http://www.alzheimersresearchuk.org). Recruitment was supported by the English National Institute of Health Research (NIHR) Clinical Stroke Research Network. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Competing interests: The authors have declared that no competing interests exist.

Introduction

Cerebral small vessel disease (SVD) causes about 20% of ischaemic stroke, and is the major cause of vascular cognitive impairment, including dementia. Whilst magnetic resonance imaging (MRI) techniques play a crucial role in diagnosis they are also being used increasingly to monitor disease progression and have been proposed as surrogate markers to assess the effects of therapeutic interventions. This latter use is of potential importance in detecting change given that cognitive testing has been shown to be relatively insensitive in measuring longitudinal decline [1,2].

Conventional MRI hallmarks of SVD include T2-white matter hyperintensities (WMH) and lacunes of presumed vascular origin. However, diffusion tensor imaging (DTI) has been shown to be a more sensitive measure of white matter (WM) damage in SVD detecting abnormalities not only within WMH but also in apparently normal appearing white matter (NAWM) [3,4]. In cross-sectional studies correlations between DTI parameters and levels of disability and/or cognition have been found in a number of SVD cohorts [4–9], and these correlations are stronger than those with T2 WMH volume [6,10]. To date, however, there is little longitudinal DTI data to determine the sensitivity of different DTI parameters to change over time. Determining the magnitude and variability of change over time is needed to determine the feasibility of using DTI as a biomarker, as well as to allow estimation of sample sizes needed to identify therapeutic efficacy in clinical trials.

Several quantitative measurements can be derived from DTI [11,12]. Mean diffusivity (MD) describes the overall extent of water diffusion, is orientation independent and is a marker of WM ultrastructure. Fractional anisotropy (FA) provides information on the directionality of the diffusion tensor and therefore the organisation of and damage to the ultrastructure determined by axonal fibre packing, degree of myelination and axon diameter [13]. Axial diffusivity (AD) is measured along the principal axis of the diffusion tensor and is considered a putative axonal damage marker; whereas radial diffusivity (RD) is measured perpendicular to the principal axis and thought to be sensitive to the degree of hindrance that diffusing water molecules experience due to the axonal membrane and myelin [14,15]. Analysing the various anisotropy and diffusivity measurements provides a range of information on the underlying WM tissue, but with varying reliabilities and sensitivities.

This prospective longitudinal study aimed to determine the sensitivity of different DTI parameters to WM damage over a period of 3 years in a clinically defined cohort of symptomatic SVD patients with MRI confirmation of a corresponding small subcortical infarct as well as confluent WMH. It also used the DTI data to evaluate sample sizes required for clinical trials to demonstrate treatment efficacy, and determine which DTI parameters provide the most statistical power to detect change in such studies.

Methods

2.1 Patients

Patients presenting with symptomatic SVD were recruited as part of the St George’s Cognition and Neuroimaging in Stroke (SCANS) study [6] from stroke services in three local South London hospitals (St George’s Hospital, King’s College Hospital and St Thomas’ Hospital). Symptomatic SVD was defined as a clinical lacunar stroke syndrome [16] with a neuroanatomically corresponding small subcortical infract as well as confluent areas of WMH of presumed vascular origin on MRI graded ≥2 on the semi-quantitative Fazekas scale [17]. Exclusion criteria were: cardio-embolic sources of stroke, large vessel disease with >50% stenosis in intra- or extra-cranial arteries, other central nervous system disorders, major psychiatric disorders and any cause of white matter disease other than SVD. Patients were recruited at least 3 months after their last stroke to avoid the influence of any acute ischaemic effects on neuropsychological function (the latter analysed in a separate study). Ethical approval for the study was obtained from the Wandsworth Research Ethics Committee and all patients gave written informed consent.

One hundred and twenty one patients were recruited into the SCANS study [6,18,19] and of these, 99 completed one or more follow-up scanning sessions; reasons for exit of the other patients were death (7), declining further MR scanning (6), formally withdrawing from the total study (6), and being lost to follow-up (3) (see Fig 1 and missing data section for follow-up scanning frequencies). For this study DTI analyses were performed on the 99 patients with follow-up MRI scans available.

Download:

Fig 1. SCANS longitudinal imaging cohort.

All available DTI scans without EPI warping were used in our analyses. Exit event coding: Death = patient passed away; Dementia = patient developed a primary dementia syndrome; Heart implant = patient was unable to do further MR sessions due to a heart implant; ICH = subject met exclusion criteria after intracerebral haemorrhage; Cognition only = patient withdrew from MR sessions but remained in study for cognition; Withdrawal = patient formally withdrew from the study; Loss to follow-up = patient was lost to follow-up; Major stroke = patient met exclusion criteria after major ischaemic stroke; Cardiac arrest = major cognitive impairment related to hypoxia was seen in a patient after a cardiac arrest.

https://doi.org/10.1371/journal.pone.0147836.g001

2.2 MRI acquisition protocol

MRI scans were performed using a 1.5-Tesla GE Signa HDxt system (General Electric, Milwaukee, WI, USA) with a maximum gradient amplitude of 33 mTm^-1 and a proprietary head coil. Patients were neutrally positioned in the head coil using an alignment marker at the nasal bridge to standardise head position. Foam pads and a Velcro strap across the forehead were used to minimise head movement during the scan.

The imaging protocol lasted approximately 45 minutes and included the following scan sequences which all provided whole head coverage: I. axial Fluid Attenuated Inversion Recovery (FLAIR): repetition time (TR) = 9000ms, echo time (TE) = 130ms, inversion time (TI) = 2200ms, 28 slices of 5mm thickness without slice gap, field of view (FOV) = 240x240mm², matrix = 256x192; II. coronal T1-weighted spoiled gradient recalled echo (SPGR): TR = 11.5ms, TE = 5ms, 176 slices without slice gap, FOV = 240x240mm², matrix = 256x192, flip angle = 18°, providing 1.1mm³ isotropic voxels; and III. axial single shot spin echo planar diffusion-weighted imaging: TR = 15600ms, TE = 93.4ms, 55 slices without slice gap with isotropic voxels of 2.5mm³, FOV = 240x240 mm², matrix = 96x96, 8 non-diffusion-weighted images (b = 0 smm^-2) followed by diffusion-weighted volumes with diffusion gradients applied (b = 1000 smm^-2) in 25 non-collinear directions and the negative of these.

2.3 MRI analysis

2.3.1 Structural pre-processing and tissue segmentation.

The steps used are detailed in Lambert et al [20], but are summarised below. Initially the T1-weighted and FLAIR images were rigidly aligned to MNI template space. Using SPM12 New Segment [21] the FLAIR and T1-weighted images were segmented into 1mm³ isotropic voxels, providing tissue class probability maps (TPMs) for WM, grey matter (GM), and cerebrospinal fluid (CSF). These segmentations were used to calculate deformation fields using a diffeomorphic registration algorithm (Shoot toolbox SPM12 [22]). Following skull stripping of the T1-weighted and FLAIR images, performed using the combined tissue segmentations at a threshold of 0.1, the T1-weighted and FLAIR images were warped to the group-average-template computed from this study cohort. A modified multivariate mixture of Gaussians SVD population specific TPMs [23] were generated to improve segmentation accuracy from the T1-weighted images for GM, WM and CSF. Additional TPMs for WMH were generated using the FLAIR and T1-weighted images together. SPM New Segment was subsequently re-run on the native space images using the population specific TPMs creating segmentation maps of GM, WM, CSF and WMH. To guarantee optimal correspondence of the maps with lesions as seen on the FLAIR image, WMH segmentation maps were optimised by binarisation at a manually set threshold for each patient by a blinded single rater. Manual refinement was performed where necessary to further improve accuracy. Reliability metrics were calculated as per Dunn et al [24]. The inter-rater reliability metrics were calculated using the mean values for each rater. WMH reliability metrics were checked using a subset of 20 randomly selected subjects by two raters (EAZ and CL). All potentially available time points were used for each subject in the reliability calculations. Both raters were blinded with respect to the subject identity and time point of each scan, and the scan order was randomized for each individual rater to avoid any potential bias. Intra-rater reliability metrics were: Standard error of mean = 247mm³, mean variability = 7.78% (standard deviation (SD) = 3.03%), Pearson’s intra-class correlation coefficient = 0.98. The inter-rater reliability metrics were: Standard error of mean = 187mm³, mean variability = 5.57% (SD = 3.99%), Pearson’s intraclass correlation coefficient = 0.99. Areas of lacunar damage were manually identified on native space T1-weighted images with reference to the FLAIR images, and segmented using a space-filling algorithm in ITK-SNAP (www.itksnap.org [25]). Grey scale intensity thresholds for the space-filling algorithm were set between 350 and 500 with 500 iterations. Manual refinement was required to ensure accurate representation of lacunes of presumed vascular origin adjacent to the lateral ventricles using ITK-SNAP. Reliability metrics were checked using a subset of 20 randomly selected scans by two raters (PB and CL), using the same procedure as the WMH reliability metrics. The intra-rater reliability metrics were: Standard error of mean = 3mm³, mean variability = 7.93% (SD = 4.89%), Pearson’s intra-class correlation coefficient = 0.99. The corresponding inter-rater reliability metrics were: Standard error of mean = 2mm³, mean variability = 4.32% (SD = 4.19%), Pearson’s intraclass correlation coefficient = 0.99.

The regions of lacunes of presumed vascular origin were removed from all tissue segmentation maps and classed as a fifth tissue class. All images were visually inspected to ensure consistent segmentation quality. The similar intensity profile of GM and WMH on T1-weighted images can cause tissue misclassification between GM and WM. To avoid this, the GM and WM maps were corrected using the binarised WMH and lacune of presumed vascular origin segmentations [20]. The NAWM tissue class was calculated by subtracting the WMH from the corrected WM map (Fig 2). GM, NAWM and WMH volumes were calculated in individual space by summing the corrected segmented areas. Total brain volume was defined as the sum of these three volumes.

Download:

Fig 2. A. T1-weighted image B. FLAIR image C. Corresponding tissue segmentation map.

Abbreviations: GM = grey matter; NAWM = normal appearing white matter; Lac. = lacune of presumed vascular origin; WMH = white matter hyperintensity.

https://doi.org/10.1371/journal.pone.0147836.g002

2.3.2 Diffusion pre-processing.

Diffusion-weighted images were pre-processed to remove eddy current distortions using the FSL Linear Image Registration Tool (FLIRT, FMRIB Software Library, www.fmrib.ox.ac.uk/fsl [26]). The positive and negative diffusion gradient direction images were geometrically averaged to eliminate gradient cross-terms [27]. The 8 b = 0 images were averaged to give a T2-weighted echo planar image, throughout referred to as b0 images. Affine transformations, followed by non-linear transformations, were computed between intensity inverted, skull stripped T1-weighted images and the skull stripped b0 image using FMRIB Non-linear Image Registration Tool (FNIRT) [28]. The transformations were subsequently applied to the native T1-weighted images, masks and all tissue segmentations to obtain tissue class images in b0 space. The deformation quality was assessed for each individual using the warped T1-weighted and b0 images. Masks included voxels with probabilities of GM, NAWM or WMH >0.5. Spurious CSF voxels were removed from the tissue masks by application of a diffusivity threshold (i.e. voxels with MD >0.0026 mm²s^-1 were considered to contain CSF and excluded).

2.3.3 DTI histogram analysis.

Using DTIfit in FSL a diffusion tensor model was fitted for each voxel from which tensor-derived FA, MD, RD and AD maps were subsequently calculated (Fig 3). Diffusion characteristics of NAWM, and NAWM plus WMH (All_WM) were evaluated. Additionally, to quantify the stability of DTI parameters over time, MD histogram measures were computed in CSF spaces. Normalised histograms of 1000 bins were produced for FA, MD, RD and AD maps of the WM segmentations (FA range 0–1, bin width 0.001; MD, AD and RD range 0–0.004, bin width 4×10⁻⁶). Since FA and diffusivity measures are non-normally distributed in WM, we evaluated longitudinal change of the following histogram parameters: median, peak value, peak height, skew and kurtosis. These parameters have also been shown to have associations with cognition in a healthy older adult and SVD population [6,29,30].

Download:

Fig 3. A. An axial fractional anisotropy and B. mean diffusivity map (mm²s^-1).

Images are from a 78 year old male from the SCANS cohort and representative of DTI-derived maps of symptomatic cerebral small vessel disease patients.

https://doi.org/10.1371/journal.pone.0147836.g003

2.3.4 Missing MRI data.

Out of the 99 patients included in the current study, for each of the consecutive three year follow ups, scanning sessions were performed on 98, 77 and 71 patients, with DTI respectively available for 99, 96, 74 and 68 patients (Fig 1). Disproportionate geometrical warping in diffusion-weighted images hindered accurate co-registration of four cases at baseline and two at year 1. Reasons given for withdrawal from the study are provided in Fig 1.

During follow-up three patients suffered new clinical strokes. Of these one patient had a small cortical haemorrhage which was considered a study endpoint. The other two patients had new small subcortical infarcts and were eligible to remain in the study according to the protocol, one was however unable to do so due to disability and withdrew.

2.4 Statistical analysis

Differences between baseline demographics and neuroimaging characteristics of the 99 patients included in the current study and those who did not attend any follow-up MRI were evaluated using Mann-Whitney U and χ² tests where appropriate in Statistical Package for Social Sciences version 20.0 (SPSS, Chicago, IL, USA.).

Linear mixed-effects (LME) models in MLwiN (version 2.1, Centre for Multilevel Modelling, University of Bristol [31]) were used to assess the effects of time on the change in DTI parameters in All_WM and NAWM tissue segmentations. The trajectories of the DTI parameters for each individual were modelled as a linear trend across the follow-up period as a function of time (exact time since baseline assessment in years). Initial intercept and slope of each patient’s linear trajectory were allowed to vary with both fixed and random effects. Fixed effect variation was accounted for by time, and random effect variation allowed for remaining inter-individual differences. The average fixed effect slope of time represents the average annualised change rate for a parameter. The change rates of each DTI parameter were evaluated using a Wald test. A Bonferroni corrected significance level of p ≤0.0023 was used to correct for multiple comparisons (i.e. 21 comparisons). Since sensitivity of progression markers depends on a combination of significance and magnitude of change, as well as robustness of a model, residual errors were also calculated for all LME models. For comparison reasons, the change rates and residual errors are also expressed as a percentage of the average intercept, which represent the average baseline value.

Sample size estimations were carried out for all DTI parameters using longpower statistical package in R version 3.02 (www.R-project.org) [32]. For these estimations all available data was used and minimum sample sizes per treatment arm were calculated for a hypothetical 3-year trial with measurements taken at even yearly intervals. A balanced trial design was assumed. A range of percentage treatment effects was calculated assuming 80% power with α = 0.05.

DTI analysis produces many related parameters. To determine whether a more predictive score could be developed by combining these different measures, we also performed principal component analysis (PCA) using SPSS 20.0 to reduce the number of DTI parameters into one component that reflected the underlying variability in the data. PCA was applied to baseline DTI values. Component scores for subsequent time points were calculated using the component coefficients and standardised DTI values based on baseline values. To overcome the multicollinearity issue of highly correlated variables, we chose to exclude RD and AD parameters from the PCA as they can be considered as subparameters of MD. When all candidate variables were included in the PCA we observed problems with multicollinearity (correlation matrix determinant <0.0001). When we restricted variables to those that declined significantly over time in our primary analysis this problem was ameliorated. Principal components were retained if eigenvalues were greater than 1. Promax rotation was applied. Following PCA, LME models were also created for the principal component scores.

Results

3.1 Demographics

Study cohort demographics and the volumetric measures of the various tissue class segmentations are shown in Table 1. Patients not included in the longitudinal study cohort were significantly older (p = 0.008), and had a lower score on the MMSE (p = 0.001). The mean age of the analysed longitudinal study cohort was 68.9 (SD = 10.0) years at baseline, and 75.0 (SD = 6.8) years for the patients in the non-longitudinal cohort. No differences were observed in baseline brain volume and other neuroimaging characteristics between these cohorts.

Download:

Table 1. Baseline characteristics of the longitudinal study cohort.

https://doi.org/10.1371/journal.pone.0147836.t001

3.2 DTI change over time

There was a highly significant change in most DTI parameters for both NAWM and All_WM with p-values ≤0.001 (Table 2). FA and diffusivity histogram metrics changed consistently with a decrease in anisotropy and an increase in diffusivity, indicating a general increase in WM ultrastructural damage over time. The percentage annualised change of different DTI parameters ranged from 0.01% to 4.43% in All_WM tissue, and 0.15–5.96% in NAWM.

Download:

Table 2. DTI parameter progression values calculated through linear mixed effect models.

https://doi.org/10.1371/journal.pone.0147836.t002

There was a change in the shape of the FA histograms over time due to a reduction in FA which manifested as decreased median and peak value, and an increased peak height (Fig 4A). This reduction in FA values also created a greater positive skew (i.e. longer tail to the right, All_WM: annualised change rate = 0.010, p <0.001; NAWM: annualised change rate = 0.016, p <0.001) and an increased kurtosis (i.e. a stronger peaked distribution, NAWM: annualised change rate = 0.025, p <0.001). Subtle differences over time were observed in the histogram distributions of the three diffusivity measures, MD, RD and AD (Fig 4B–4D). Increased diffusivity led to increased peak value and median, and decreased peak height values over time. Changes in skew and kurtosis were not significant for MD, but for RD and AD there was an increase in the skewness due to greater frequency of higher diffusivity voxels in the rightward tail (All_WM: annualised change rate RD skew = 0.037, p <0.001; annualised change rate AD skew = 0.061, p <0.001; NAWM: annualised change rate RD skew = 0.042, p <0.001; annualised change rate AD skew = 0.073, p <0.001). An increase in kurtosis was also found for the AD distribution representing a more peaked distribution (All_WM: annualised change rate = 0.230, p <0.001; NAWM: annualised change rate = 0.778, p <0.001).

Download:

Fig 4. Average DTI histogram distributions of the baseline longitudinal cohort and at 3 year follow-up.

A. Fractional anisotropy. B. Mean diffusivity (mm²s^-1). C. Radial diffusivity (mm²s^-1). D. Axial diffusivity (mm²s^-1).

https://doi.org/10.1371/journal.pone.0147836.g004

A significant change of most DTI parameters was also observed within regions of WMH (S1 Table). The direction of change was, however, opposite to that observed in NAWM and All_WM. The estimated progression rates of FA metrics indicated increased anisotropy over time. The FA histogram showed an increased median and peak value, and decreased peak height (median FA: annualised change rate = 0.00332, p < 0.001; peak value FA: annualised change rate = 0.00648, p = 0.002; peak height FA: annualised change rate = -0.000171, p < 0.001). Changes in histogram parameters of MD, RD and AD within WMH over time indicated reduced diffusivity, as evidenced by decreased medians (MD: annualised change rate = -8.20×10⁻⁶ mm²s^-1, p < 0.001; RD: annualised change rate = -7.64×10⁻⁶ mm²s^-1, p < 0.001) and peak values (MD: annualised change rate = -1.46×10⁻⁵ mm²s^-1, p < 0.001; RD: annualised change rate = -1.06×10⁻⁵ mm²s^-1, p < 0.001; AD: annualised change rate = -1.08×10⁻⁵ mm²s^-1, p = 0.002). Histogram peak height was reduced over time in all diffusivity measures and skew increased.

Within the CSF the MD median, peak height, peak value and kurtosis remained stable over time (see S2 Table). A significant change was found in MD skew in the CSF over time, such that the skew became less negative, indicating that the MD distribution within the CSF became more Gaussian over time (annualised change rate = 0.00267, p = 0.009).

Percentage residual error for each model was highly varied within both All_WM and NAWM. In particular kurtosis of the diffusivity measures had high residual errors: 9.515–42.658% within All_WM, and 12.235–68.440% within NAWM. Residual errors of skew models in the diffusivity DTI parameters were also relatively high. For MD and RD, these errors were larger than the annualised skew change rates.

The extracted principal component scores for All_WM and NAWM accounted respectively for 80.6% and 74.9% of the variance in the PCA. Very large annualised change rates were observed (annualised change All_WM: 714%; NAWM: 431%). Residual errors were equally large but did not exceed annualised change rates.

3.3 Sample size estimations

Table 3 shows the results of the sample size estimations. Sample sizes required for treatment trials varied greatly using the different DTI parameters. Assuming a treatment effect of 25% the minimum sample sizes needed per treatment arm for DTI parameters of All_WM ranged from 198 for MD peak height to 186043 for FA kurtosis, the range for NAWM being 356 for MD peak height to 3059376 for MD skew. The principal component scores had smaller sample size estimations than those of the majority of DTI parameters, with an estimation of 501 for the All_WM component score and 679 within NAWM if a treatment effect of 25% was observed. MD peak height however requires the smallest sample size.

Download:

Table 3. The minimum sample size estimation per arm required to test treatment effects of various percentages using DTI parameters.

https://doi.org/10.1371/journal.pone.0147836.t003

Discussion

Our results demonstrate that DTI is sensitive to change in patients with SVD over a relatively short duration of follow-up. Observed changes in DTI parameters were consistent with a progressive degradation of WM ultrastructure showing a decline in anisotropy and increase in diffusivity in WM. However our analysis demonstrates that different DTI parameters show marked differences in the sensitivity to change and errors. We have shown that the most sensitive marker for SVD progression is MD peak height. This finding indicates greater stability of MD values compared to its radial and axial counterparts, as well as sensitivity to change and a good model fit. We also provide sample size estimates, which will be useful in planning clinical trials using DTI as a surrogate marker. It follows that MD peak height has the smallest sample size estimate.

Over three year follow-up most DTI histogram parameters changed significantly. Interpretation of longitudinal DTI changes is not always straightforward. Changes in diffusivity histogram distributions are due to increased tissue ultrastructural damage that manifests as increases in diffusivity. This results in higher median and peak value, and a lower peak height as higher diffusivity values are identified with greater frequency. Increases in diffusivity values cause the rightward tail of the distribution to become broader, increasing the positive skew. Also, the related reduction of lower diffusivity values causes the diffusivity histograms to become less peaked and have lower kurtosis. This negative change is observed in MD and RD (All_WM), although it does not reach significance. For the remainder of kurtosis values, in RD (NAWM) and AD, kurtosis increases, significantly in the case of the latter. This finding most likely reflects the instability of AD as a reliable DTI measurement, indicated by the high residual errors in AD kurtosis models. Since AD is the diffusivity in the principal diffusion direction, errors in eigenvalue sorting and inaccuracies in the approximation of the principal diffusion direction using the DTI model will have a large impact on measurement of the AD parameter [33]. High residual errors in AD models therefore indicate limited value of this parameter for longitudinal assessment of SVD. MD, and to a lesser extent RD, may be more informative as they are calculated from the average of two or more eigenvalues and appear to be more stable measures.

Changes in the FA histogram distribution are due to reductions in FA that manifest as a lower median and peak value, and a higher peak height due to identification of lower FA values with greater frequency. In turn this causes a greater rightward tail of the distribution (greater positive skew) and a more pronounced peak (greater kurtosis). However, FA histogram parameters were found to have a poor longitudinal model fit with high residual error. This may be due to ambiguity of FA values. In particular, as FA is defined by a quadratic function in three related variables, similar FA values can represent different diffusion tensor (ellipsoidal) shapes [34]. The dependency of FA on the three eigenvalues of the diffusion tensor leads to reduced stability of FA compared to MD. The use of a component score combining MD and FA measurements inevitably decreases sensitivity to detection of disease progression compared to MD measurements alone. The observed sample size estimates of component scores also reflect this effect.

Samples sizes depend on the detected power of observations. Power in turn relies on error, so depends on measurement stability, effect size and in this case magnitude of change. Since scanner upgrades can affect measurement of DTI histogram parameters [35], their effect was evaluated. The results indicated that there was no significant change in MD histogram parameters in CSF spaces, apart from for MD skew. In this case the results indicated that the MD distribution was becoming more Gaussian over time, representing an improvement in MD quantification. These results assure validity of the observed rates of change in this study. The finding of lowest sample size estimates for MD peak height and median therefore demonstrates that the MD peak height benefits from being a highly stable measurement that is also sensitive to change. Whether DTI parameters are calculated for All_WM or NAWM also impacts sample size. Histologically WMH tissue shows axonal loss and glial proliferation [36]. DTI studies have shown that diffusivity is increased and anisotropy is lowered in WMHs compared to NAWM in SVD and CADASIL, a young-onset genetic variant of SVD [4,37]. This contributes to our finding of greater change in diffusivity parameters of All_WM tissue. Sample sizes estimates of All_WM diffusivity parameters are therefore reduced compared to parameter estimation in NAWM. For FA parameters, however, greatest change and lowest samples size estimates were obtained for FA measures from NAWM tissue compared with FA measures from All_WM. The severity of damage found in WMHs affects the uncertainty of the orientation of the principal diffusivities from which FA is calculated and will likely provide inconsistencies in FA [34]. This leads to greater sample size estimates for FA parameters within All_WM. Changes to diffusion characteristics of WMHs over time revealed the opposite effect to those observed in in NAWM and All_WM. The increase in anisotropy and decrease in diffusivity over time suggests improvement of tissue microstructure within WMHs; however, this is unlikely to be the case. Maillard et al., [38] have shown that voxels containing higher FA and lower diffusivity values are located at the edges of WMHs, indicating that the core of WMHs has greater microstructural damage than newly formed areas of WMH. Our observation of an apparent tissue microstructural improvement is likely to be the consequence of WMH growth (and changes in the surface area to volume ratio of the WMHs) and provides a misleading indication of improvement. Consequently, WMH diffusion characteristics provide an ambiguous and unsuitable marker for monitoring disease progression.

Despite relatively high patient dropout rates which could bias the sample to those with less decline, significant changes were detected in our cohort. The patient group who did not attend follow-up were older and had worse cognitive levels, indicating that such patients may be more vulnerable to decline and consequently an underestimation of actual change in DTI parameters in the study population. Secondly, the use of small histogram bin-widths in our analyses increased noise in measurements of peak value and peak height. A trade-off between sensitivity and noise is applicable when deciding on histogram bin size. We have however shown that annualized change rates of DTI parameter peak value and peak height are small and therefore require small histogram bin widths to enable detection of change. Also, we have generated histograms from WM tissue because SVD predominantly affects the subcortical structures, and in pilot studies we have shown that DTI parameters were sensitive to change in WM over short time periods in SVD [29]. However, it is now recognised that SVD does cause cortical GM pathology as well, with both atrophy and micro infarcts, and it is possible that a global whole brain DTI measure may be a simpler alternative. Finally, the absence of an age-matched control population is a limitation to this study, as observed changes in diffusion characteristics could be a result of SVD associated pathology and ageing. The ageing literature shows varying change rates of diffusion parameter. Most studies reveal smaller rates of change than found in our cohort. Nusbaum et al., [39] for example found age-related histogram changes in ADC, a diffusion weighted imaging equivalent to MD in SVD [35], that were approximately 8 times smaller than our MD changes. More recently, longitudinal FA and diffusivity changes in ageing have been shown to have annual percentage changes 2–3 times smaller than those in our cohort [40]. Lastly, a study directly comparing DTI parameters of CADASIL patients to age-matched controls revealed no change in MD peak value in the ageing population and an annual progression of 3.65×10⁻⁶ mm²s^-1 in CADASIL [41] which is slightly greater than the 2.85 over time indicated reduced diffusivity, as evidenced by decreased medians 10⁻⁶ mm²s^-1 change found in this study. Charlton et al., [30] on the other hand, reported slightly larger FA and MD changes in their elderly population. This discrepancy could be the result of quantifying change as the difference between two time points. Measurement noise and missing data will have a greater effect on change rates the fewer follow-up time points are acquired. Previously we have shown that analysis of fewer data points leads to greater magnitudes of change, as shown by smaller sample size estimates of hypothetical clinical trials [42]. Together these findings assure us that the progression rates as measured in our cohort represent SVD progression beyond pure age-related WM changes.

Our data suggests disease progression in SVD can be monitored sensitively using DTI histogram analysis. This is in line with earlier whole brain DTI histogram in ageing [30], sporadic SVD [29] and CADASIL [41], that have provided highly reproducible DTI markers suitable for longitudinal studies [43,44]. Since changes in cognition over the same duration in the same cohort have been small and significant for only one cognitive domain or even absent, depending on the statistical analysis applied [2,42], we suggest DTI could be attractive as a surrogate marker to follow disease progression and test the effect of therapeutic interventions in SVD. Other imaging parameters that are thought to be potential surrogate markers are brain atrophy, WMH growth and incident lacunes. Previously we have shown in the same cohort that development of new lacunes only occurred in a minority of SVD patients, there was however significant atrophy and WMH growth over a period of 3 years [42]. Additionally, in [42] we investigated progression of one DTI parameter, MD peak height of All_WM. Sample size estimates showed that WMH volume and MD peak height were most sensitive to change and most suitable as surrogate marker of disease progression. The current study provides evidence that MD peak height of All_WM is the most sensitive DTI parameter for use in longitudinal SVD studies.

Apart from being able to detect prospective change MRI markers will have to meet further requirements to fulfil all criteria of a validated surrogate marker [45]; one of the most important being, being able to predict future clinical cognitive decline. Such associations have been weak and inconsistent with regards to WMH volume [29,41,46]. Previous cross-sectional studies have on the other hand shown that DTI histogram parameters correlate strongly with cognition in SVD populations [6,8,29]. Analyses of longitudinal correlations between cognitive and DTI change, and investigating the interventional effect on DTI markers and cognition therefore now clearly represent the next phase for determining the potential of DTI as a future surrogate marker of SVD. The present study is continuing cognitive follow-up to 5 years, and plans to correlate cognitive changes over this period with MRI parameters if significant cognitive change is detectable over this longer period.

In summary, we observed significant longitudinal change in various DTI histogram parameters in WM in a SVD cohort of patients with symptomatic lacunar stroke and confluent areas of WMH of presumed vascular origin. The reliability and sensitivity of individual DTI histogram parameters as markers of disease progression is however highly variable. Our results suggest MD peak height is the most sensitive marker for detecting ultrastructural brain tissue damage in SVD. Future studies investigating progression of SVD and the use of DTI as a surrogate disease marker are needed and should focus on the relationship between the clinical outcomes of SVD and imaging markers.

Supporting Information

S1 Dataset. Minimal dataset underlying the presented findings.

https://doi.org/10.1371/journal.pone.0147836.s001

(XLS)

S1 Table. Progression of DTI parameters within white matter hyperintensities calculated through linear mixed effect models.

https://doi.org/10.1371/journal.pone.0147836.s002

(DOCX)

S2 Table. Progression of mean diffusivity values within CSF spaces.

https://doi.org/10.1371/journal.pone.0147836.s003

(DOCX)

Acknowledgments

We thank Dr Lalit Kalra and Dr Tony Rudd for help with recruitment. We also thank Dr Irina Chis Ster for statistical advice. And we thank all participants from the SCANS study.

Author Contributions

Conceived and designed the experiments: HSM TRB. Performed the experiments: EAZ PB AJL CL. Analyzed the data: EAZ PB OAW. Contributed reagents/materials/analysis tools: CL AJL TRB. Wrote the paper: EAZ PB CL AJL RGM TRB HSM.

References

1. Pearce LA, McClure LA, Anderson DC, Jacova C, Sharma M, Hart RG, et al. Effects of long-term blood pressure lowering and dual antiplatelet treatment on cognitive function in patients with recent lacunar stroke: a secondary analysis from the SPS3 randomised trial. Lancet Neurol. 2014 Dec;13(12):1177–85. pmid:25453457
- View Article
- PubMed/NCBI
- Google Scholar
2. Lawrence AJ, Brookes RL, Zeestraten EA, Barrick TR, Morris RG, Markus HS. Pattern and rate of cognitive decline in cerebral small vessel disease: a prospective study. PLoS One. 2015;10(8):e0135523. pmid:26273828
- View Article
- PubMed/NCBI
- Google Scholar
3. O’Sullivan M, Summers PE, Jones DK, Jarosz JM, Williams SCR, Markus HS. Normal-appearing white matter in ischemic leukoaraiosis: a diffusion tensor MRI study. Neurology. 2001;57(12):2301–10. pmid:11756617
- View Article
- PubMed/NCBI
- Google Scholar
4. Chabriat H, Pappata S, Poupon C, Clark CA, Vahedi K, Poupon F, et al. Clinical Severity in CADASIL Related to Ultrastructural Damage in White Matter: In Vivo Study With Diffusion Tensor MRI. Stroke. 1999 Dec 1;30(12):2637–43. pmid:10582990
- View Article
- PubMed/NCBI
- Google Scholar
5. Kim SH, Park JS, Ahn H-J, Seo SW, Lee J-M, Kim ST, et al. Voxel-based analysis of diffusion tensor imaging in patients with subcortical vascular cognitive impairment: correlates with cognitive and motor deficits. J Neuroimaging. 2011 Oct;21(4):317–24. pmid:20977535
- View Article
- PubMed/NCBI
- Google Scholar
6. Lawrence AJ, Patel B, Morris RG, MacKinnon AD, Rich PM, Barrick TR, et al. Mechanisms of cognitive impairment in cerebral small vessel disease: multimodal MRI results from the St George’s cognition and neuroimaging in stroke (SCANS) study. PLoS One. 2013 Jan;8(4):e61014. pmid:23613774
- View Article
- PubMed/NCBI
- Google Scholar
7. O’Sullivan M, Jones DK, Summers PE, Morris RG, Williams SCR, Markus HS. Evidence for cortical “disconnection” as a mechanism of age-related cognitive decline. Neurology. 2001 Aug 28;57(4):632–8. pmid:11524471
- View Article
- PubMed/NCBI
- Google Scholar
8. O’Sullivan M, Morris RG, Huckstep B, Jones DK, Williams SCR, Markus HS. Diffusion tensor MRI correlates with executive dysfunction in patients with ischaemic leukoaraiosis. J Neurol Neurosurg Psychiatry. 2004 Mar 1;75(3):441–7. pmid:14966162
- View Article
- PubMed/NCBI
- Google Scholar
9. O’Sullivan M, Singhal S, Charlton RA, Markus HS. Diffusion tensor imaging of thalamus correlates with cognition in CADASIL without dementia. Neurology. 2004 Mar 9;62(5):702–7. pmid:15007117
- View Article
- PubMed/NCBI
- Google Scholar
10. Schmidt R, Berghold A, Jokinen H, Gouw AA, van der Flier WM, Barkhof F, et al. White matter lesion progression in LADIS: frequency, clinical effects, and sample size calculations. Stroke. 2012 Oct;43(10):2643–7. pmid:22879094
- View Article
- PubMed/NCBI
- Google Scholar
11. Le Bihan D, Mangin JF, Poupon C, Clark CA, Pappata S, Molko N, et al. Diffusion tensor imaging: concepts and applications. J Magn Reson Imaging. 2001 Apr;13(4):534–46. pmid:11276097
- View Article
- PubMed/NCBI
- Google Scholar
12. Basser PJ, Jones DK. Diffusion-tensor MRI: theory, experimental design and data analysis—a technical review. NMR Biomed. 2002;15(7–8):456–67. pmid:12489095
- View Article
- PubMed/NCBI
- Google Scholar
13. Beaulieu C. The basis of anisotropic water diffusion in the nervous system—a technical review. NMR Biomed. 2002;15(7–8):435–55. pmid:12489094
- View Article
- PubMed/NCBI
- Google Scholar
14. Song S-K, Sun S-W, Ramsbottom MJ, Chang C, Russell J, Cross AH. Dysmyelination Revealed through MRI as Increased Radial (but Unchanged Axial) Diffusion of Water. Neuroimage. 2002 Nov;17(3):1429–36. pmid:12414282
- View Article
- PubMed/NCBI
- Google Scholar
15. Song S-K, Sun S-W, Ju W-K, Lin S-J, Cross AH, Neufeld AH. Diffusion tensor imaging detects and differentiates axon and myelin degeneration in mouse optic nerve after retinal ischemia. Neuroimage. 2003 Nov;20(3):1714–22. pmid:14642481
- View Article
- PubMed/NCBI
- Google Scholar
16. Bamford J, Sandercock P, Jones L, Warlow C. The natural history of lacunar infarction: the Oxfordshire Community Stroke Project. Stroke. 1987 May 1;18(3):545–51. pmid:3590244
- View Article
- PubMed/NCBI
- Google Scholar
17. Fazekas F, Chawluk JB, Alavi A, Hurtig HI, Zimmerman RA. MR signal abnormalities at 1.5 T in Alzheimer’s dementia and normal aging. Am J Roentgenol. 1987;149:351–6.
- View Article
- Google Scholar
18. Benjamin P, Lawrence AJ, Lambert C, Patel B, Chung AW, MacKinnon AD, et al. Strategic lacunes and their relationship to cognitive impairment in cerebral small vessel disease. NeuroImage Clin. 2014 Jan;4:828–37. pmid:24936433
- View Article
- PubMed/NCBI
- Google Scholar
19. Patel B, Lawrence AJ, Chung AW, Rich P, Mackinnon AD, Morris RG, et al. Cerebral microbleeds and cognition in patients with symptomatic small vessel disease. Stroke. 2013 Feb;44(2):356–61. pmid:23321452
- View Article
- PubMed/NCBI
- Google Scholar
20. Lambert C, Narean JS, Benjamin P, Zeestraten EA, Barrick TR, Markus HS. Characterising the grey matter correlates of leukoaraiosis in cerebral small vessel disease. NeuroImage Clin. 2015;9:194–205. pmid:26448913
- View Article
- PubMed/NCBI
- Google Scholar
21. Ashburner J, Friston KJ. Unified segmentation. Neuroimage. 2005 Jul 1;26(3):839–51. pmid:15955494
- View Article
- PubMed/NCBI
- Google Scholar
22. Ashburner J, Friston KJ. Diffeomorphic registration using geodesic shooting and Gauss-Newton optimisation. Neuroimage. Elsevier Inc.; 2011 Apr 1;55(3):954–67. pmid:21216294
- View Article
- PubMed/NCBI
- Google Scholar
23. Lambert C, Lutti A, Helms G, Frackowiak R, Ashburner J. Multiparametric brainstem segmentation using a modified multivariate mixture of Gaussians. NeuroImage Clin. 2013 Jan;2:684–94. pmid:24179820
- View Article
- PubMed/NCBI
- Google Scholar
24. Dunn JT, Clark-Papasavas C, Marsden P, Baker S, Cleij M, Kapur S, et al. Establishing test-retest reliability of an adapted [18F]fallypride imaging protocol in older people. J Cereb blood flow Metab. 2013 Jul;33(7):1098–103. pmid:23591647
- View Article
- PubMed/NCBI
- Google Scholar
25. Yushkevich PA, Piven J, Hazlett HC, Smith RG, Ho S, Gee JC, et al. User-guided 3D active contour segmentation of anatomical structures: significantly improved efficiency and reliability. Neuroimage. 2006 Jul 1;31(3):1116–28. pmid:16545965
- View Article
- PubMed/NCBI
- Google Scholar
26. Jenkinson M, Smith S. A global optimisation method for robust affine registration of brain images. Med Image Anal. 2001 Jun;5(2):143–56. pmid:11516708
- View Article
- PubMed/NCBI
- Google Scholar
27. Neeman M, Freyer JP, Sillerud LO. A simple method for obtaining cross-term-free images for diffusion anisotropy studies in NMR microimaging. Magn Reson Med. 1991 Sep;21(1):138–43. pmid:1943671
- View Article
- PubMed/NCBI
- Google Scholar
28. Andersson JLR, Jenkinson M, Smith S. Non-linear registration aka Spatial normalisation. FMRIB Technical Report TR07JA2. FMRIB Technical Report TR07JA2. 2007.
29. Nitkunan A, Barrick TR, Charlton RA, Clark CA, Markus HS. Multimodal MRI in cerebral small vessel disease: its relationship with cognition and sensitivity to change over time. Stroke. 2008 Jul;39(7):1999–2005. pmid:18436880
- View Article
- PubMed/NCBI
- Google Scholar
30. Charlton RA, Schiavone F, Barrick TR, Morris RG, Markus HS. Diffusion tensor imaging detects age related white matter change over a 2 year follow-up which is associated with working memory decline. J Neurol Neurosurg Psychiatry. 2010 Jan;81(1):13–9. pmid:19710050
- View Article
- PubMed/NCBI
- Google Scholar
31. Rasbash J, Charlton C, Browne WJ, Healy M, Cameron B. MLwiN Version 2.1. Centre for Multilevel Modelling, University of Bristol; 2009.
32. Liu G, Liang K-Y. Sample size calculations for studies with correlated observations. Biometrics. 1997;53(3):937–47. pmid:9290224
- View Article
- PubMed/NCBI
- Google Scholar
33. Wheeler-Kingshott CAM, Cercignani M. About “axial” and “radial” diffusivities. Magn Reson Med. 2009 May;61(5):1255–60. pmid:19253405
- View Article
- PubMed/NCBI
- Google Scholar
34. Green HAL, Pena A, Price CJ, Warburton EA, Pickard JD, Carpenter TA, et al. Increased Anisotropy in Acute Stroke: A Possible Explanation. Stroke. 2002 Jun 1;33(6):1517–21. pmid:12052984
- View Article
- PubMed/NCBI
- Google Scholar
35. Gunda B, Porcher R, Duering M, Guichard J-P, Mawet J, Jouvent E, et al. ADC Histograms from Routine DWI for Longitudinal Studies in Cerebral Small Vessel Disease: A Field Study in CADASIL. PLoS One. 2014;9(5):e97173. pmid:24819368
- View Article
- PubMed/NCBI
- Google Scholar
36. Pantoni L, Garcia JH. Pathogenesis of leukoaraiosis: a review. Stroke. 1997;28:652–9. pmid:9056627
- View Article
- PubMed/NCBI
- Google Scholar
37. Van Norden AGW, De Laat KF, Van Dijk EJ, Van Uden IWM, Van Oudheusden LJB, Gons RAR, et al. Diffusion tensor imaging and cognition in cerebral small vessel disease: the RUN DMC study. Biochim Biophys Acta. 2012 Mar;1822(3):401–17. pmid:21549191
- View Article
- PubMed/NCBI
- Google Scholar
38. Maillard P, Fletcher E, Harvey D, Carmichael O, Reed B, Mungas D, et al. White Matter Hyperintensity Penumbra. Stroke. 2011;42(7):1917–22. pmid:21636811
- View Article
- PubMed/NCBI
- Google Scholar
39. Nusbaum AO, Tang CY, Buchsbaum MS, Wei TC, Atlas SW. Regional and global changes in cerebral diffusion with normal aging. AJNR Am J Neuroradiol. 2001 Jan;22(1):136–42. pmid:11158899
- View Article
- PubMed/NCBI
- Google Scholar
40. Sexton CE, Walhovd KB, Storsve AB, Tamnes CK, Westlye LT, Johansen-Berg H, et al. Accelerated Changes in White Matter Microstructure during Aging: A Longitudinal Diffusion Tensor Imaging Study. J Neurosci. 2014;34(46):15425–36. pmid:25392509
- View Article
- PubMed/NCBI
- Google Scholar
41. Holtmannspötter M, Peters N, Opherk C, Martin D, Herzog J, Brückmann H, et al. Diffusion magnetic resonance histograms as a surrogate marker and predictor of disease progression in CADASIL: a two-year follow-up study. Stroke. 2005 Dec;36(12):2559–65. pmid:16269644
- View Article
- PubMed/NCBI
- Google Scholar
42. Benjamin P, Zeestraten EA, Lambert C, Chis Ster IC, Williams OA, Lawrence AJ, et al. Progression of MRI markers in cerebral small vessel disease: sample size considerations for clinical trials. J Cereb Blood Flow Metab. 2015 Jun 3;1–8.
- View Article
- Google Scholar
43. Cercignani M, Bammer R, Sormani MP, Fazekas F, Filippi M. Inter-sequence and inter-imaging unit variability of diffusion tensor MR imaging histogram-derived metrics of the brain in healthy volunteers. Am J Neuroradiol. 2003;24:638–43. pmid:12695195
- View Article
- PubMed/NCBI
- Google Scholar
44. Rovaris M, Filippi M, Calori G, Rodegher M, Campi A, Colombo B, et al. Intra-observer reproducibility in measuring new putative MR markers of demyelination and axonal loss in multiple sclerosis: a comparison with conventional T2-weighted images. J Neurol. 1997;244:266–70. pmid:9112597
- View Article
- PubMed/NCBI
- Google Scholar
45. Fazekas F, Ropele S, Schmidt R. Can small-vessel disease-related cerebral abnormalities be used as a surrogate marker for vascular dementia trials? J Neural Transm Suppl. 2002;(62):61–7. pmid:12456051
- View Article
- PubMed/NCBI
- Google Scholar
46. Liem MK, Oberstein SAJL, Haan J, Van der Neut IL, Ferrari MD, Van Buchem MA, et al. MRI correlates of cognitive decline in CADASIL: A 7-year follow-up study. Neurology. 2009;72(2):143–8. pmid:19139365
- View Article
- PubMed/NCBI
- Google Scholar

[ref1] 1. Pearce LA, McClure LA, Anderson DC, Jacova C, Sharma M, Hart RG, et al. Effects of long-term blood pressure lowering and dual antiplatelet treatment on cognitive function in patients with recent lacunar stroke: a secondary analysis from the SPS3 randomised trial. Lancet Neurol. 2014 Dec;13(12):1177–85. pmid:25453457
View Article
PubMed/NCBI
Google Scholar

[2] View Article

[3] PubMed/NCBI

[4] Google Scholar

[ref2] 2. Lawrence AJ, Brookes RL, Zeestraten EA, Barrick TR, Morris RG, Markus HS. Pattern and rate of cognitive decline in cerebral small vessel disease: a prospective study. PLoS One. 2015;10(8):e0135523. pmid:26273828
View Article
PubMed/NCBI
Google Scholar

[6] View Article

[7] PubMed/NCBI

[8] Google Scholar

[ref3] 3. O’Sullivan M, Summers PE, Jones DK, Jarosz JM, Williams SCR, Markus HS. Normal-appearing white matter in ischemic leukoaraiosis: a diffusion tensor MRI study. Neurology. 2001;57(12):2301–10. pmid:11756617
View Article
PubMed/NCBI
Google Scholar

[10] View Article

[11] PubMed/NCBI

[12] Google Scholar

[ref4] 4. Chabriat H, Pappata S, Poupon C, Clark CA, Vahedi K, Poupon F, et al. Clinical Severity in CADASIL Related to Ultrastructural Damage in White Matter: In Vivo Study With Diffusion Tensor MRI. Stroke. 1999 Dec 1;30(12):2637–43. pmid:10582990
View Article
PubMed/NCBI
Google Scholar

[14] View Article

[15] PubMed/NCBI

[16] Google Scholar

[ref5] 5. Kim SH, Park JS, Ahn H-J, Seo SW, Lee J-M, Kim ST, et al. Voxel-based analysis of diffusion tensor imaging in patients with subcortical vascular cognitive impairment: correlates with cognitive and motor deficits. J Neuroimaging. 2011 Oct;21(4):317–24. pmid:20977535
View Article
PubMed/NCBI
Google Scholar

[18] View Article

[19] PubMed/NCBI

[20] Google Scholar

[ref6] 6. Lawrence AJ, Patel B, Morris RG, MacKinnon AD, Rich PM, Barrick TR, et al. Mechanisms of cognitive impairment in cerebral small vessel disease: multimodal MRI results from the St George’s cognition and neuroimaging in stroke (SCANS) study. PLoS One. 2013 Jan;8(4):e61014. pmid:23613774
View Article
PubMed/NCBI
Google Scholar

[22] View Article

[23] PubMed/NCBI

[24] Google Scholar

[ref7] 7. O’Sullivan M, Jones DK, Summers PE, Morris RG, Williams SCR, Markus HS. Evidence for cortical “disconnection” as a mechanism of age-related cognitive decline. Neurology. 2001 Aug 28;57(4):632–8. pmid:11524471
View Article
PubMed/NCBI
Google Scholar

[26] View Article

[27] PubMed/NCBI

[28] Google Scholar

[ref8] 8. O’Sullivan M, Morris RG, Huckstep B, Jones DK, Williams SCR, Markus HS. Diffusion tensor MRI correlates with executive dysfunction in patients with ischaemic leukoaraiosis. J Neurol Neurosurg Psychiatry. 2004 Mar 1;75(3):441–7. pmid:14966162
View Article
PubMed/NCBI
Google Scholar

[30] View Article

[31] PubMed/NCBI

[32] Google Scholar

[ref9] 9. O’Sullivan M, Singhal S, Charlton RA, Markus HS. Diffusion tensor imaging of thalamus correlates with cognition in CADASIL without dementia. Neurology. 2004 Mar 9;62(5):702–7. pmid:15007117
View Article
PubMed/NCBI
Google Scholar

[34] View Article

[35] PubMed/NCBI

[36] Google Scholar

[ref10] 10. Schmidt R, Berghold A, Jokinen H, Gouw AA, van der Flier WM, Barkhof F, et al. White matter lesion progression in LADIS: frequency, clinical effects, and sample size calculations. Stroke. 2012 Oct;43(10):2643–7. pmid:22879094
View Article
PubMed/NCBI
Google Scholar

[38] View Article

[39] PubMed/NCBI

[40] Google Scholar

[ref11] 11. Le Bihan D, Mangin JF, Poupon C, Clark CA, Pappata S, Molko N, et al. Diffusion tensor imaging: concepts and applications. J Magn Reson Imaging. 2001 Apr;13(4):534–46. pmid:11276097
View Article
PubMed/NCBI
Google Scholar

[42] View Article

[43] PubMed/NCBI

[44] Google Scholar

[ref12] 12. Basser PJ, Jones DK. Diffusion-tensor MRI: theory, experimental design and data analysis—a technical review. NMR Biomed. 2002;15(7–8):456–67. pmid:12489095
View Article
PubMed/NCBI
Google Scholar

[46] View Article

[47] PubMed/NCBI

[48] Google Scholar

[ref13] 13. Beaulieu C. The basis of anisotropic water diffusion in the nervous system—a technical review. NMR Biomed. 2002;15(7–8):435–55. pmid:12489094
View Article
PubMed/NCBI
Google Scholar

[50] View Article

[51] PubMed/NCBI

[52] Google Scholar

[ref14] 14. Song S-K, Sun S-W, Ramsbottom MJ, Chang C, Russell J, Cross AH. Dysmyelination Revealed through MRI as Increased Radial (but Unchanged Axial) Diffusion of Water. Neuroimage. 2002 Nov;17(3):1429–36. pmid:12414282
View Article
PubMed/NCBI
Google Scholar

[54] View Article

[55] PubMed/NCBI

[56] Google Scholar

[ref15] 15. Song S-K, Sun S-W, Ju W-K, Lin S-J, Cross AH, Neufeld AH. Diffusion tensor imaging detects and differentiates axon and myelin degeneration in mouse optic nerve after retinal ischemia. Neuroimage. 2003 Nov;20(3):1714–22. pmid:14642481
View Article
PubMed/NCBI
Google Scholar

[58] View Article

[59] PubMed/NCBI

[60] Google Scholar

[ref16] 16. Bamford J, Sandercock P, Jones L, Warlow C. The natural history of lacunar infarction: the Oxfordshire Community Stroke Project. Stroke. 1987 May 1;18(3):545–51. pmid:3590244
View Article
PubMed/NCBI
Google Scholar

[62] View Article

[63] PubMed/NCBI

[64] Google Scholar

[ref17] 17. Fazekas F, Chawluk JB, Alavi A, Hurtig HI, Zimmerman RA. MR signal abnormalities at 1.5 T in Alzheimer’s dementia and normal aging. Am J Roentgenol. 1987;149:351–6.
View Article
Google Scholar

[66] View Article

[67] Google Scholar

[ref18] 18. Benjamin P, Lawrence AJ, Lambert C, Patel B, Chung AW, MacKinnon AD, et al. Strategic lacunes and their relationship to cognitive impairment in cerebral small vessel disease. NeuroImage Clin. 2014 Jan;4:828–37. pmid:24936433
View Article
PubMed/NCBI
Google Scholar

[69] View Article

[70] PubMed/NCBI

[71] Google Scholar

[ref19] 19. Patel B, Lawrence AJ, Chung AW, Rich P, Mackinnon AD, Morris RG, et al. Cerebral microbleeds and cognition in patients with symptomatic small vessel disease. Stroke. 2013 Feb;44(2):356–61. pmid:23321452
View Article
PubMed/NCBI
Google Scholar

[73] View Article

[74] PubMed/NCBI

[75] Google Scholar

[ref20] 20. Lambert C, Narean JS, Benjamin P, Zeestraten EA, Barrick TR, Markus HS. Characterising the grey matter correlates of leukoaraiosis in cerebral small vessel disease. NeuroImage Clin. 2015;9:194–205. pmid:26448913
View Article
PubMed/NCBI
Google Scholar

[77] View Article

[78] PubMed/NCBI

[79] Google Scholar

[ref21] 21. Ashburner J, Friston KJ. Unified segmentation. Neuroimage. 2005 Jul 1;26(3):839–51. pmid:15955494
View Article
PubMed/NCBI
Google Scholar

[81] View Article

[82] PubMed/NCBI

[83] Google Scholar

[ref22] 22. Ashburner J, Friston KJ. Diffeomorphic registration using geodesic shooting and Gauss-Newton optimisation. Neuroimage. Elsevier Inc.; 2011 Apr 1;55(3):954–67. pmid:21216294
View Article
PubMed/NCBI
Google Scholar

[85] View Article

[86] PubMed/NCBI

[87] Google Scholar

[ref23] 23. Lambert C, Lutti A, Helms G, Frackowiak R, Ashburner J. Multiparametric brainstem segmentation using a modified multivariate mixture of Gaussians. NeuroImage Clin. 2013 Jan;2:684–94. pmid:24179820
View Article
PubMed/NCBI
Google Scholar

[89] View Article

[90] PubMed/NCBI

[91] Google Scholar

[ref24] 24. Dunn JT, Clark-Papasavas C, Marsden P, Baker S, Cleij M, Kapur S, et al. Establishing test-retest reliability of an adapted [18F]fallypride imaging protocol in older people. J Cereb blood flow Metab. 2013 Jul;33(7):1098–103. pmid:23591647
View Article
PubMed/NCBI
Google Scholar

[93] View Article

[94] PubMed/NCBI

[95] Google Scholar

[ref25] 25. Yushkevich PA, Piven J, Hazlett HC, Smith RG, Ho S, Gee JC, et al. User-guided 3D active contour segmentation of anatomical structures: significantly improved efficiency and reliability. Neuroimage. 2006 Jul 1;31(3):1116–28. pmid:16545965
View Article
PubMed/NCBI
Google Scholar

[97] View Article

[98] PubMed/NCBI

[99] Google Scholar

[ref26] 26. Jenkinson M, Smith S. A global optimisation method for robust affine registration of brain images. Med Image Anal. 2001 Jun;5(2):143–56. pmid:11516708
View Article
PubMed/NCBI
Google Scholar

[101] View Article

[102] PubMed/NCBI

[103] Google Scholar

[ref27] 27. Neeman M, Freyer JP, Sillerud LO. A simple method for obtaining cross-term-free images for diffusion anisotropy studies in NMR microimaging. Magn Reson Med. 1991 Sep;21(1):138–43. pmid:1943671
View Article
PubMed/NCBI
Google Scholar

[105] View Article

[106] PubMed/NCBI

[107] Google Scholar

[ref28] 28. Andersson JLR, Jenkinson M, Smith S. Non-linear registration aka Spatial normalisation. FMRIB Technical Report TR07JA2. FMRIB Technical Report TR07JA2. 2007.

[ref29] 29. Nitkunan A, Barrick TR, Charlton RA, Clark CA, Markus HS. Multimodal MRI in cerebral small vessel disease: its relationship with cognition and sensitivity to change over time. Stroke. 2008 Jul;39(7):1999–2005. pmid:18436880
View Article
PubMed/NCBI
Google Scholar

[110] View Article

[111] PubMed/NCBI

[112] Google Scholar

[ref30] 30. Charlton RA, Schiavone F, Barrick TR, Morris RG, Markus HS. Diffusion tensor imaging detects age related white matter change over a 2 year follow-up which is associated with working memory decline. J Neurol Neurosurg Psychiatry. 2010 Jan;81(1):13–9. pmid:19710050
View Article
PubMed/NCBI
Google Scholar

[114] View Article

[115] PubMed/NCBI

[116] Google Scholar

[ref31] 31. Rasbash J, Charlton C, Browne WJ, Healy M, Cameron B. MLwiN Version 2.1. Centre for Multilevel Modelling, University of Bristol; 2009.

[ref32] 32. Liu G, Liang K-Y. Sample size calculations for studies with correlated observations. Biometrics. 1997;53(3):937–47. pmid:9290224
View Article
PubMed/NCBI
Google Scholar

[119] View Article

[120] PubMed/NCBI

[121] Google Scholar

[ref33] 33. Wheeler-Kingshott CAM, Cercignani M. About “axial” and “radial” diffusivities. Magn Reson Med. 2009 May;61(5):1255–60. pmid:19253405
View Article
PubMed/NCBI
Google Scholar

[123] View Article

[124] PubMed/NCBI

[125] Google Scholar

[ref34] 34. Green HAL, Pena A, Price CJ, Warburton EA, Pickard JD, Carpenter TA, et al. Increased Anisotropy in Acute Stroke: A Possible Explanation. Stroke. 2002 Jun 1;33(6):1517–21. pmid:12052984
View Article
PubMed/NCBI
Google Scholar

[127] View Article

[128] PubMed/NCBI

[129] Google Scholar

[ref35] 35. Gunda B, Porcher R, Duering M, Guichard J-P, Mawet J, Jouvent E, et al. ADC Histograms from Routine DWI for Longitudinal Studies in Cerebral Small Vessel Disease: A Field Study in CADASIL. PLoS One. 2014;9(5):e97173. pmid:24819368
View Article
PubMed/NCBI
Google Scholar

[131] View Article

[132] PubMed/NCBI

[133] Google Scholar

[ref36] 36. Pantoni L, Garcia JH. Pathogenesis of leukoaraiosis: a review. Stroke. 1997;28:652–9. pmid:9056627
View Article
PubMed/NCBI
Google Scholar

[135] View Article

[136] PubMed/NCBI

[137] Google Scholar

[ref37] 37. Van Norden AGW, De Laat KF, Van Dijk EJ, Van Uden IWM, Van Oudheusden LJB, Gons RAR, et al. Diffusion tensor imaging and cognition in cerebral small vessel disease: the RUN DMC study. Biochim Biophys Acta. 2012 Mar;1822(3):401–17. pmid:21549191
View Article
PubMed/NCBI
Google Scholar

[139] View Article

[140] PubMed/NCBI

[141] Google Scholar

[ref38] 38. Maillard P, Fletcher E, Harvey D, Carmichael O, Reed B, Mungas D, et al. White Matter Hyperintensity Penumbra. Stroke. 2011;42(7):1917–22. pmid:21636811
View Article
PubMed/NCBI
Google Scholar

[143] View Article

[144] PubMed/NCBI

[145] Google Scholar

[ref39] 39. Nusbaum AO, Tang CY, Buchsbaum MS, Wei TC, Atlas SW. Regional and global changes in cerebral diffusion with normal aging. AJNR Am J Neuroradiol. 2001 Jan;22(1):136–42. pmid:11158899
View Article
PubMed/NCBI
Google Scholar

[147] View Article

[148] PubMed/NCBI

[149] Google Scholar

[ref40] 40. Sexton CE, Walhovd KB, Storsve AB, Tamnes CK, Westlye LT, Johansen-Berg H, et al. Accelerated Changes in White Matter Microstructure during Aging: A Longitudinal Diffusion Tensor Imaging Study. J Neurosci. 2014;34(46):15425–36. pmid:25392509
View Article
PubMed/NCBI
Google Scholar

[151] View Article

[152] PubMed/NCBI

[153] Google Scholar

[ref41] 41. Holtmannspötter M, Peters N, Opherk C, Martin D, Herzog J, Brückmann H, et al. Diffusion magnetic resonance histograms as a surrogate marker and predictor of disease progression in CADASIL: a two-year follow-up study. Stroke. 2005 Dec;36(12):2559–65. pmid:16269644
View Article
PubMed/NCBI
Google Scholar

[155] View Article

[156] PubMed/NCBI

[157] Google Scholar

[ref42] 42. Benjamin P, Zeestraten EA, Lambert C, Chis Ster IC, Williams OA, Lawrence AJ, et al. Progression of MRI markers in cerebral small vessel disease: sample size considerations for clinical trials. J Cereb Blood Flow Metab. 2015 Jun 3;1–8.
View Article
Google Scholar

[159] View Article

[160] Google Scholar

[ref43] 43. Cercignani M, Bammer R, Sormani MP, Fazekas F, Filippi M. Inter-sequence and inter-imaging unit variability of diffusion tensor MR imaging histogram-derived metrics of the brain in healthy volunteers. Am J Neuroradiol. 2003;24:638–43. pmid:12695195
View Article
PubMed/NCBI
Google Scholar

[162] View Article

[163] PubMed/NCBI

[164] Google Scholar

[ref44] 44. Rovaris M, Filippi M, Calori G, Rodegher M, Campi A, Colombo B, et al. Intra-observer reproducibility in measuring new putative MR markers of demyelination and axonal loss in multiple sclerosis: a comparison with conventional T2-weighted images. J Neurol. 1997;244:266–70. pmid:9112597
View Article
PubMed/NCBI
Google Scholar

[166] View Article

[167] PubMed/NCBI

[168] Google Scholar

[ref45] 45. Fazekas F, Ropele S, Schmidt R. Can small-vessel disease-related cerebral abnormalities be used as a surrogate marker for vascular dementia trials? J Neural Transm Suppl. 2002;(62):61–7. pmid:12456051
View Article
PubMed/NCBI
Google Scholar

[170] View Article

[171] PubMed/NCBI

[172] Google Scholar

[ref46] 46. Liem MK, Oberstein SAJL, Haan J, Van der Neut IL, Ferrari MD, Van Buchem MA, et al. MRI correlates of cognitive decline in CADASIL: A 7-year follow-up study. Neurology. 2009;72(2):143–8. pmid:19139365
View Article
PubMed/NCBI
Google Scholar

[174] View Article

[175] PubMed/NCBI

[176] Google Scholar

Figures

Abstract

Introduction

Methods

2.1 Patients

2.2 MRI acquisition protocol

2.3 MRI analysis

2.3.1 Structural pre-processing and tissue segmentation.

2.3.2 Diffusion pre-processing.

2.3.3 DTI histogram analysis.

2.3.4 Missing MRI data.

2.4 Statistical analysis

Results

3.1 Demographics

3.2 DTI change over time

3.3 Sample size estimations

Discussion

Supporting Information

S1 Dataset. Minimal dataset underlying the presented findings.

S1 Table. Progression of DTI parameters within white matter hyperintensities calculated through linear mixed effect models.

S2 Table. Progression of mean diffusivity values within CSF spaces.

Acknowledgments

Author Contributions

References