ORIGINAL RESEARCH article
Age-Related Changes of Peak Width Skeletonized Mean Diffusivity (PSMD) Across the Adult Lifespan: A Multi-Cohort Study
- 1Institute of Neurodegenerative Diseases (IMN), CNRS, CEA, Bordeaux, France
- 2Institute of Neurodegenerative Diseases (IMN), University of Bordeaux, Bordeaux, France
- 3Bordeaux Population Health Research Center, Inserm, Bordeaux, France
- 4Institute of Neuroscience and Medicine (INM-1), Research Centre Juelich, Juelich, Germany
- 5Institute for Anatomy I, Medical Faculty, Heinrich Heine University Dusseldorf, Dusseldorf, Germany
- 6Research Institute, The Hospital for Sick Children, Toronto, ON, Canada
- 7Department of Physiology and Nutritional Sciences, University of Toronto, Toronto, ON, Canada
- 8Bloorview Research Institute, Holland Bloorview Kids Rehabilitation Hospital, Toronto, ON, Canada
- 9Departments of Psychology and Psychiatry, University of Toronto, Toronto, ON, Canada
- 10Department of Neurology, Medical University of Graz, Graz, Austria
- 11Centre for Healthy Brain Ageing (CHeBA), School of Psychiatry, UNSW Medicine, University of New South Wales, Sydney, NSW, Australia
- 12Neuropsychiatric Institute, Neuropsychiatric Institute Prince of Wales Hospital, Randwick, NSW, Australia
- 13Mathematics and Statistics, Murdoch University, Perth, WA, Australia
- 14Centre for Cognitive Ageing and Cognitive Epidemiology, Department of Psychology, University of Edinburgh, Edinburgh, United Kingdom
- 15Brain Research Imaging Centre, Neuroimaging Sciences, The University of Edinburgh, Edinburgh, United Kingdom
- 16Departmet of Neurology, Max Planck Institute for Human Cognitive and Brain Sciences, Leipzig, Germany
- 17Institute for Stroke and Dementia Research (ISD), University Hospital, LMU Munich, Munich, Germany
- 18Department of Neurology, Bordeaux University Hospital, Bordeaux, France
Parameters of water diffusion in white matter derived from diffusion-weighted imaging (DWI), such as fractional anisotropy (FA), mean, axial, and radial diffusivity (MD, AD, and RD), and more recently, peak width of skeletonized mean diffusivity (PSMD), have been proposed as potential markers of normal and pathological brain ageing. However, their relative evolution over the entire adult lifespan in healthy individuals remains partly unknown during early and late adulthood, and particularly for the PSMD index. Here, we gathered and analyzed cross-sectional diffusion tensor imaging (DTI) data from 10 population-based cohort studies in order to establish the time course of white matter water diffusion phenotypes from post-adolescence to late adulthood. DTI data were obtained from a total of 20,005 individuals aged 18.1 to 92.6 years and analyzed with the same pipeline for computing skeletonized DTI metrics from DTI maps. For each individual, MD, AD, RD, and FA mean values were computed over their FA volume skeleton, PSMD being calculated as the 90% peak width of the MD values distribution across the FA skeleton. Mean values of each DTI metric were found to strongly vary across cohorts, most likely due to major differences in DWI acquisition protocols as well as pre-processing and DTI model fitting. However, age effects on each DTI metric were found to be highly consistent across cohorts. RD, MD, and AD variations with age exhibited the same U-shape pattern, first slowly decreasing during post-adolescence until the age of 30, 40, and 50 years, respectively, then progressively increasing until late life. FA showed a reverse profile, initially increasing then continuously decreasing, slowly until the 70s, then sharply declining thereafter. By contrast, PSMD constantly increased, first slowly until the 60s, then more sharply. These results demonstrate that, in the general population, age affects PSMD in a manner different from that of other DTI metrics. The constant increase in PSMD throughout the entire adult life, including during post-adolescence, indicates that PSMD could be an early marker of the ageing process.
Parameters of water diffusion in white matter derived from diffusion-weighted imaging (DWI), such as fractional anisotropy (FA), mean, axial, and radial diffusivity (MD, AD, and RD) are well-established markers of normal brain maturation (1–5) and ageing (6–12) and have been proposed as potential tools for the investigation of various brain disorders (13–19).
More recently, peak width of skeletonized mean diffusivity (PSMD) (20), a new phenotype of white matter microstructure that can be derived from DWI, has been proposed as an imaging biomarker of small vessel disease (SVD) (20, 21) and a correlate of cognitive impairment, particularly processing speed (20–22). Recall that PSMD is the difference between the 5th and the 95th percentiles of the distribution of the voxel MD value across a skeleton of the brain white matter. Note that PSMD is a dispersion statistic, as opposed to other diffusion tensor imaging (DTI) metric, such as AD, MD, RD, and FA, that are central tendency statistics. So far, our knowledge of the PSMD distribution in healthy individuals has been limited to these three previously mentioned studies that all included people aged over 50 years. In addition, none of these studies addressed the issue of changes in PSMD across lifespan, which is critical for establishing whether PSMD could be used as an imaging marker of brain aging as well as an early predictor of age-related disorders or to serve as a tool to monitor outcomes in clinical trials. Here, we gathered and analyzed cross-sectional DTI data from 10 population-based cohort studies to establish the time course, from post-adolescence to late adulthood, of the PSMD distribution and compare it with that of more commonly used white matter water diffusion phenotypes in white matter.
Materials and Methods
Ten independent data sets coming from cross-sectional cohort studies were gathered in the present study, namely, MRi-Share, BIL&GIN, SYS, LIFE-Adult, 1000 BRAINS, UKBiobank, ASPSF, OATS, LBC1936, MAS (see acronym definition in Table 1 caption). All but three (LIFE-Adult, 1000BRAINS, and UKBiobank) were part of the BRIDGET Consortium (BRain Imaging, cognition, Dementia, and next generation GEnomics: a Transdisciplinary approach to search for risk and protective factors of neurodegenerative disease), supported by EU-JPND (European Union Joint Programme—Neurodegenerative Disease Research). The 10 data sets included a total of 20,005 individuals (age range, 18.1 to 92.6 years; 10,807 women and 9,198 men). Table 1 detail sample size and age distribution for the 10 cohorts that were all of cross.
Diffusion-Weighted Image Acquisition and Preprocessing
Table 2 summarize the main characteristics of the DWI acquisition and preprocessing for the 10 cohorts. Overall, there was considerable variability between studies regarding almost all acquisition parameters, including scanner manufacturer, field strength, gradient strength, diffusion pulse sequence, resolution, and number of directions. For the present work, it was not possible to access raw DWI data at different sites in order to harmonize processing from the initial DICOM data. For this reason, DWI data sets were pre-processed with procedures specific to each site, including exclusion of data upon visual detection of major artifacts due to eddy current distortions or head motion. AD, RD, MD, and FA maps were computed by fitting the DTI model parameters in each voxel from these preprocessed DWI volumes. Additional details on DWI preprocessing and DTI parameter map computation for each data set are provided in the Supplementary Material section.
Derivation of DTI Metrics
Various metrics were derived from the DTI data using a script developed by Baykara et al. (http://www.psmd-marker.com) (20). This original script was designed to extract PSMD, an index of the dispersion of MD values across the white matter skeleton. Briefly, the computation included two steps: 1) WM skeletonizing using the FA map, and 2) analyzing the voxel value distribution histogram in the MD volume masked by the WM-FA skeleton. The 3D brain image of FA of each individual was skeletonized using the FSL-TBSS software, part of the FMRIB Software Library (FSL) (23, 24), using the FMRIB 1 mm FA template and applying a 0.2 threshold on FA maps. Then the MD volume of the same individual was masked, keeping only voxels within the FA skeleton. Furthermore, in order to reduce contamination of the skeleton by CSF voxels, the FA-masked MD volumes were further masked by both a standard FA skeleton with a threshold of 0.3 on FA values and a custom mask (provided with the PSMD software tool) designed so as to exclude regions adjacent to the ventricles, such as the fornix. Finally, PSMD was computed as the difference between the 95th and 5th percentiles of the so-masked MD volume voxel value distribution. Here, we extended this script in order to obtain, in addition to PSMD values, estimates of the mean values of axial, radial, and mean diffusivity (AD, RD, MD, respectively) as well as of FA over the same customized skeleton. All 10 cohorts were processed separately with this customized script, and the results were sent to the Bordeaux site where they were combined for further statistical analysis.
Age Category Definition
Due to previously reported non-linear effects of age on DTI metrics (1, 3, 8), we divided each cohort sample into subsamples of 10-year age range starting at 18 years of age, the last subsample (i.e. [78 to 98]) including all subjects older than 78 years, as there was only a small number of individuals older than 88 years. Table 1 and Figure 1 detail the contribution of each cohort to each age category. Because we planned running analyses at the age category by cohort by sex level, we discarded subsamples of small sizes, namely, containing less than 30 individuals.
Figure 1 Box-plots and histogram distribution of the age of the participants of the 10 contributing data sets.
Assessing Age-Related Changes of PSMD and Other DTI Metrics
For each of the five DTI metrics and each age category, we performed an analysis of variance including “age” as the main effect, and “sex,” total intracranial volume (TIV), and “Cohort” as confounding factors. The Cohort effect was included in order to account for apparent large differences in DTI metric average values across cohorts contributing to the same age category data set (see Figure 2). “Sex” and “TIV” were included as covariates since mixed results have been reported regarding the impact of sex and TIV on DTI measures [(12, 25–28), see review in (3)]. In order to document differences of age effects between cohorts contributing to the same age category, we also performed an analysis of variance for each age category and each cohort, including “age,” “sex,” and “TIV” as effects. Moreover, we analyzed the effects of age on DTI metrics in men and women separately.
Figure 2 Mean and standard deviation (bars) of the five DTI metrics for each age subcategory and each data set. The dotted lines connect values for the different age categories of the same data set.
Assessing the Effects of Sex and TIV on PSMD and Other DTI Metrics
For each of the five DTI metrics and each age category, we also performed an analysis of variance including “Sex” and “TIV” as main factors and “Cohort” as confounding factors.
All statistical analyses were performed using the JMP Pro Software (version 14.3.0, SAS Institute Inc.).
Tables 3A–E provide basic statistics across age categories and cohorts for PSMD and the four other DTI metrics, while Figure 2 illustrates their respective profiles across the adult lifespan. From both Table 3 and Figure 2, it is apparent that there is considerable variability in all five DTI parameter values, and that within a given age category the variability across cohorts is larger than the variability between individuals of the same cohort (see the extreme case of AD and FA values for the LBC1936 study performed at 1.5T, for example).
Table 3A Basic statistics of axial diffusivity (AD, in mm2 s−1 × 10−4) average across a white matter skeleton for each age category and each contributing data set.
Table 3B Basic statistics of radial diffusivity (RD, in mm2 s−1 × 10−4) average across a white matter skeleton for each age category and each contributing data set.
Table 3C Basic statistics of mean diffusivity (MD, in mm2 s−1 × 10−4) average across a white matter skeleton for each age category and each contributing data set.
Table 3D Basic statistics of peak width skeletonized mean diffusivity (PSMD, in mm2 s−1 × 10−4) for each age category and each contributing data set.
Table 3E Basic statistics of fractional anisotropy (FA, unitless) average across individual FA skeleton for each age range category and each contributing data set.
Figure 3 compares the inter-individual variability of PSMD (in terms of its coefficient of variation, CV in %) within each cohort to those of the other DTI metrics, again for each age category and each cohort, revealing that PSMD CV's is in the order of 10% to 15% (with values as high as 20% for later ages) while those for AD, RD, MD, and FA are in the order of 2% to 5%. Note also that the CVs of all DTI metrics increase with age, more for PSMD than for the other parameters.
Figure 3 Coefficients of variation of the five DTI metrics for each age subcategory and each data set.
Effects of Age on PSMD and Other DTI Metrics
The evolution of PSMD across the adult life is different from that of the other metrics (Figure 2). Indeed, PSMD seems to increase monotonically with age, whereas AD, RD, and MD exhibit a similar J-shape profile, initially slightly decreasing during post-adolescence before later increasing during adulthood. As for FA, it shows a reverse profile to that of AD, RD, and MD, with an initial small increase followed by a later decrease.
This apparent specific lifespan profile of PSMD was confirmed by the quantitative estimates of the effects of age on PSMD and other DTI metrics provided by the between-cohort ANOVA (see Table 4) and their profiles of evolution across age categories as shown in Figure 4. Estimates of the age effect on PSMD were indeed positive for all age categories and significant for all but the 28- to 38-year age period. This increase in PSMD accelerated during late life periods, its value being multiplied by a factor of 3 between the 58- to 68-year and the 78- to 98-year periods. By contrast, AD, RD, and MD age variation profiles were characterized by an initial small but significant decrease (negative age effect), followed by a stable period (non-significant age effect) and a final significant increase (positive age effect). What distinguished AD, RD, and MD profiles was their respective timings, the initial decrease in RD and MD being significant only for the 18 to 28 years subsample, while it extended over the 38- to 48-year period for AD. Note that since the stable period covered two periods of 10 years for all three metrics, the increase in AD was delayed by 10 years compared with RD and MD. At the same time, FA exhibited the reverse contrast consisting of an initial increase during the 18- to 28-year period followed by a stable period (28 to 38 years) before an accelerated decrease during the rest of the adult life course. The age effects on DTI metrics were not significantly different between men and women as revealed by the separate sex-specific analyses of age effects.
Table 4 ANOVA effect of age (estimate (standard error) and significance p-value) on five DTI metrics (AD, RD, MD, PSMD, and FA evaluated across individual FA skeletons).
Figure 4 ANOVA estimates and 95% confidence intervals of the effects of age on the five DT metrics for each age subcategory. Units are mm2/s/year × 10−3 for AD, RD, and MD, mm2/s/year × 10−4 for PSMD, year × 10−3 for FA.
Effects of Sex and TIV on DTI
Amplitude of sex effects on DTI metric average values were found to be quite variable across the various cohorts for the different age categories (see Figure 5). When pooling all data sets, we found that sex had significant effects on all DTI parameters except on FA. Women had higher mean AD and MD values than men, who conversely had higher PSMD values than women (see Table 5). As for TIV, we found that it had positive effects on all DTI parameters (see Table 5), including PSMD, these effects being very significant (p < 10−4) in all cases but RD (p = 0.53).
Figure 5 ANOVA estimates and 95% confidence intervals of the effects of sex (females minus males) on the five DT metrics for each age subcategory and each data set.
Table 5 ANOVA effect of sex and TIV (estimate (standard error) and significance p value) on five DTI metrics (AD, RD, MD, PSMD, and FA evaluated across individual FA skeletons).
We will first discuss methodological issues and potential limitations in interpreting results of our study. Second, we will discuss what the present study adds to the already existing literature on age effects on classical DTI parameters. In the third part, we will focus on the original findings regarding PSMD distribution and evolution over the adult life.
Methodological Issues and Potential Limitations of the Present Study
In order to study the effects of age over the entire adult lifespan, we gathered DTI data from 20,005 individuals scanned at 10 different sites with the major objective of maximizing statistical power. By so doing, we were aware that the variability of DTI parameters derived from the entire data set would be much larger than for a single-site study because of between-site differences in data acquisition protocols, such as scanner manufacturer, field strength, number and strength of diffusion-encoding gradients, image voxel size, post-processing of raw DWI data, and so on (23, 24, 29–32). Accordingly, a key issue was to reduce as much as possible the variability due to sources that were controllable so that we could maximize benefits of the very large sample size.
Concerns with DTI multi-site studies have been emphasized by several authors (33, 34) and strategies have been proposed for harmonizing DWI data before proceeding to statistical analysis of DTI parameter maps. However, it has also been reported that DTI parameters estimated over the whole white matter or large regions of interest actually exhibit high intra- and inter-scanner reproducibility, making them suitable for multisite studies without extensive harmonization (17, 35, 36). This is the approach that we implemented here as it was not possible to access raw DWI data at different sites in order to harmonize processing from the very earliest stages (see Supplementary Materials). In our study, harmonization was limited to the post-processing of DTI maps pre-computed by each site using the same script for generating individual skeletons and computing individual DTI parameter values across these skeletons. However, our findings indicate that it is possible to draw meaningful conclusions from such minimal harmonization by focusing on the effects of interest (age in our case) rather than the absolute values of the measurements, an approach previously used by others (37).
Multi-site studies do have some advantages as compared with single-site investigations, namely, statistical power due to very large sample size, and the ability to recruit a sufficient number of individuals over the entire life or adult lifespan. As emphasized by others, these are important features in the context of clinical studies that are often multisite in nature (17, 35). In fact, there have been several previous reports of multisite studies on DTI parameters in adults but they dealt with issues other than age effects such as methodology (17, 30, 33) and genetic effects (38, 39), for example.
In this study, we focused on the mean DTI metrics over a white matter skeleton, rather than over the entire white-matter compartment, since our primary metric of interest was PSMD, which is defined over the FSL-TBSS white matter skeleton derived from FA data. This approach should in theory minimize any partial volume effects by limiting the measurement over the core of white matter tracts. Indeed, when we compared the mean DTI metrics over the white matter skeleton with those measured over the whole white-matter compartment in two cohorts of the study (MRi-Share and BIL&GIN), mean FA values were higher and diffusivity values were lower when averaged over the white matter skeleton than when averaged over the entire white-matter compartment. However, the cross-subject variability measured as the CV of the mean DTI metrics decreased only marginally (a fraction of a percent for CVs ranging from 2% to 4%) when mean values were computed over the white matter skeleton rather than the whole white matter. In contrast, the CV of PSMD decreased markedly (on average, from 24% to 8%) when its computation was performed on the white matter skeleton rather than on than the global white-matter mask. This demonstrates the importance of choosing a measure of MD dispersion values over a white matter skeleton for controlling between subject variability.
In the present work, we restricted the analysis to classical DTI metrics as only two of the contributing cohorts had high angular resolution and/or multi-shell acquisition schemes that could be used for estimating advanced white-matter microstructural parameters with more sophisticated models (29). Here, DWI data processing was solely based on the classical DTI model. The DTI model limitations are well known (24), and it has been shown, for example, that correction for free water has a major impact on classical DTI parameter values (40, 41). However, although investigating advanced white-matter microstructural parameters is highly desirable, it was beyond the scope of our study: it would require additional data sets with multi-shell acquisition, especially for individuals aged 30 to 50 years or over 70 years, in order to supplement existing data (8) on the adult lifespan trajectory of these microstructural parameters.
Mixed results have been reported regarding the impact of sex and TIV on DTI measures [(12, 25–28), see review in (3)]. Here, we also observed mixed results across the different cohorts, although very significant sex effects on all DTI parameters, except FA, were uncovered when combining the entire data set. Note, however, these sex effects were of very small size (for PSMD, for example, omega2 = 6.7 × 10−3 for the sex effect to be compared with 7.6 × 10−2 for the age effect), which could explain the mixed findings in the literature, and suggests further investigations are required in order to understand their biological origins. TIV effects on DTI parameters are not well established in the literature. In our study, we found that TIV was positively correlated with all DTI parameters except for RD. Similar to sex, TIV effects when significant were very small (for PSMD again, omega2 = 1.5 × 10−3 for the TIV effect). Here again additional investigations are needed to understand the origins of these effects.
Finally, and importantly, it should be stressed that interpretation of the results of the present study should be taken with caution because of the cross-sectional nature of the data that we analyzed. Numerous reports have indeed pointed out the caveats of cross-sectional design for assessing effects of age and demonstrated how such design may lead to spurious findings when compared to those obtained with longitudinal data (12, 42, 43). In the present work, there was no attempt to use a single model to describe the variation of DTI parameters with age over the entire adulthood period. Rather, we selected a piecewise linear model to examine/compare age-related changes in 10-year duration consecutive time bins, thereby minimizing the generation bias between cohorts of nearby categories. Understandably, such an approach does not eradicate the intrinsic limits of our cross-sectional study. But it should be reminded that a fully longitudinal design is quasi impossible to implement in the context of lifespan research, since, in practice, measures in an individual can be repeated only a few times and at short duration intervals. As a consequence, such longitudinal studies suffer from some of the limitations of cross-sectional ones. This may explain why the results of the present study are compatible with those a previously published longitudinal study (12) in which individuals aged between 20 and 84 years were observed twice 3 years apart.
Adult Lifespan Profiles of Variation of Classical DTI Parameters AD, RD, MD, and FA
Effects of age on white-matter microstructure assessed with DTI have been intensively investigated over the past decade from a developmental perspective (4) as well as in a lifespan/ageing framework (1, 3, 5, 8, 11, 12, 27, 40). Briefly, and considering only DTI metrics estimated at the global level, AD, RD, and thus, MD were reported to follow similar U- or better J-shape age variation patterns, initially decreasing during childhood and adolescence (see (4) for review) then exhibiting an accelerated increase during the adult life (8, 12), while FA followed a reverse profile. Our own findings agree with this body of results during the adult life course. Raw data plots show J-shape profiles for AD, RD, and MD, and the reverse profile for FA, as well as acceleration of these changes during late life.
Maximum global FA values and minimum global MD, RD, and AD values have been reported to occur before the age of 40 years (1, 5, 11, 12, 27), although large variations were found when considering individual tracts (5, 26, 44, 45). Here, we found extreme values for RD, MD, and FA occurring between 28 and 38 years, well in line with these previous findings. In addition, we found that the decrease of AD in the post-adolescence period extended into adulthood by about 10 years more than for RD and MD, thereby uncovering heterochrony of AD and RD variations during adulthood. Such a heterochrony during adulthood was not detected in a previous longitudinal study (12) possibly due to an insufficient sample size and large DTI metric variability between individuals (as can be seen in Figure 7 of the mentioned report). Note that two recent studies (46, 47) have reported opposite age effects for AD (decrease) and RD (increase) with stable MD during the 18 to 55 year age period; however, as both studies used simple linear modeling due to small sample sizes, no age at extreme value could be observed. Rather, our findings are compatible with the AD-RD variation heterochrony that has been noticed earlier during childhood and adolescence at the individual tract level with stronger decrease for RD than for AD (34, 44). According to these and our findings, the AD decrease/RD increase profile (6) would occur only during mid-adulthood.
PSMD Is a Diffusion Imaging Phenotype With a Profile of Variation Across the Adult Lifespan That Differs From That of Other DTI Parameters
The distribution of PSMD values observed for the different cohorts and age categories of the present study are consistent with the few comparable data reported in the literature for older participants (no data are available in young adults). For example, Baykara et al. reported in their pioneering article a PSMD median value around 3.0 (in mm2 s−1 × 10−4, range [2.5, 4.9]) in a sample of healthy individuals aged 60 to 80 years drawn from the ASPF cohort [see Table 2 of Baykara et al. (20)], values that are comparable to those reported in Table 3D of our study in subsamples of other cohorts of similar age category. Similarly, Wei et al. (21) recently reported a PSMD average value of 2.4 × 10−4 mm2 s−1 in a sample of healthy controls aged around 60 years. In both studies, PSMD CVs were close to 10%, a value again similar to those observed in our own study. That the CV of PSMD is two to three times larger than the CVs of other DTI metrics could be expected since PSMD is a dispersion rather than a central tendency statistic. Moreover, the larger increase in PSMD CV as age advances (as compared to the other DTI metrics) indicates that this phenotype should be used with caution especially during the late life period. However, it is important to note that the CV of PSMD was found to be quite stable across cohorts with similar age ranges.
The main goal of the present study was to document the profile of PSMD evolution across age bands during adulthood. In this respect, and the proviso that the data we gathered were not longitudinal, our results show that PSMD increases continuously from post-adolescence to late adult life, that this increase is accelerating at later ages, and that this acceleration is larger than for the other DTI metrics. As there are no available data of PSMD in childhood and adolescence, it is not possible to decide whether the lifetime PSMD evolution profile is similar to those of AD, RD, MD, i.e. with a decrease during childhood that reaches the minimum value before adulthood, or if it shows continuous increase throughout the lifespan. Nevertheless, it remains the case that the continuous and accelerating increase of PSMD during adulthood is an indication that it is an adequate and potentially valuable marker of white matter ageing. In particular, it is notable that PSMD increases during early adulthood when the other DTI metrics variations appear to be still undergoing late maturational processes.
The biological mechanisms of the origin of PSMD evolution with age are at present unknown, but one can think of several reasons why PSMD may be more prone to increase with age as compared with the other metrics. First, it is important to remember that PSMD is a measure of MD values dispersion across a skeleton of white matter. As such, it will be directly affected by differences across MD values of the individual tracts. Consequently, regional heterogeneity as well as heterochrony in MD values of the fiber tracts will result in higher PSMD values more than in average MD values. Second, MD itself is a weighted average of AD (1/3) and RD (2/3) values, and thus MD value dispersion will also be affected by heterochrony in AD and RD variations with age. Overall, what possibly makes PSMD an early and sensitive imaging marker of ageing is that it captures multiple sources of heterogeneity in white matter water diffusion parameters. With this regard, it would be interesting to investigate variations in the pattern of MD dispersion at the regional level using tract-based DTI metrics since it is well established that heterochrony is a major feature of the development and aging of the different fiber tracts [see for example (5, 8, 26, 44)]. Accordingly, variations of PSMD value provide only a gross and possibly biased estimate of the white matter-microstructure dynamics. We did not implement regional analysis as our study focused on PSMD that is by definition a dispersion statistic over the entire white matter skeleton. Nevertheless, a regional approach would certainly be interesting and feasible since peak width of MD values could be measured on a white-matter skeleton at the tract level in the same manner as it has been done for other DTI metrics (see (34, 46) for example).
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation, to any qualified researcher.
The 10 cohort studies involving human participants were reviewed and approved by the following ethic committees: Comité de Protection des Personnes Sud-Ouest (MRi-Share), Comité d’Ethique de Basse-Normandie (BIL&GIN), Research Ethics Committee of the Chicoutimi Hospital (SYS), Ethics committee of the University of Duisburg-Essen (1000BRAINS), Ethics board of the Medical Faculty of the University of Leipzig (LIFE), North West Multi-center Research Ethics Committee (UKBiobank), Ethics Committee of the Medical University Graz (ASPS), Ethics committees of the Australian Twin Registry, University of New South Wales, University of Melbourne, Queensland Institute of Medical Research and the South-Eastern Sydney and Illawarra Area Health Service (OATS), Multi-Centre Research Ethics Committee for Scotland and the Lothian Research Ethics Committee (LBC1936), Ethics Committee of the University of New South Wales (MAS). The participants provided written informed consent to participate in these studies.
Study conception: BM and SD. Data collection: GB, CT, SC, ZP, TP, RS, PS, HB, NK, JT, ID, AW, AV, and BM. Data analysis: GB, LPe, SC, JS, YP, LPi, PS, WW, NA, MB, SM, AW, MD, and BM. Drafting: GB, AT, and BM. Revising the manuscript: LPe, CT, SC, JS, ZP, YP, TP, RS, LPi, PS, WW, NA, ID, MB, JW, SM, AW, AV, MD, SD, and BM).
GB has been supported by an EU-ERC starting grant (SEGWAY, PI S Debette, with funding from the European Union's Horizon 2020 research and innovation program under grant agreement No 640643). AT is supported by a grant from the Fondation pour la Recherche Médicale (DIC202161236446). SC was supported by the Initiative and Networking Fund of the Helmholtz Association and the European Union's Horizon 2020 Research and Innovation Program under Grant Agreement 785907 (Human Brain Project SGA2). The BRIDGET project is supported by the Fondation Leducq (Transatlantic Network of Excellence on the Pathogenesis of SVD of the Brain) and is an EU Joint Program -Neurodegenerative Disease Research (JPND) project. The project is supported through the following funding organizations under the aegis of JPND—www.jpnd.eu: Australia, National Health and Medical Research Council, Austria, Federal Ministry of Science, Research and Economy; Canada, Canadian Institutes of Health Research; France, French National Research Agency; Germany, Federal Ministry of Education and Research; Netherlands, The Netherlands Organization for Health Research and Development; United Kingdom, Medical Research Council. This project has received funding from the European Union's Horizon 2020 research and innovation program under grant agreement No 643417. This project has also received funding from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation program under grant agreement No 640643.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpsyt.2020.00342/full#supplementary-material
1. Hasan KM, Sankar A, Halphen C, Kramer LA, Brandt ME, Juranek J, et al. Development and organization of the human brain tissue compartments across the lifespan using diffusion tensor imaging. Neuroreport (2007) 18:1735–9. doi: 10.1097/WNR.0b013e3282f0d40c
2. Giorgio A, Watkins KE, Chadwick M, James S, Winmill L, Douaud G, et al. Longitudinal changes in grey and white matter during adolescence. NeuroImage (2010) 49:94–103. doi: 10.1016/j.neuroimage.2009.08.003
3. Lebel C, Gee M, Camicioli R, Wieler M, Martin W, Beaulieu C. Diffusion tensor imaging of white matter tract evolution over the lifespan. NeuroImage (2012) 60:340–52. doi: 10.1016/j.neuroimage.2011.11.094
6. Bennett IJ, Madden DJ, Vaidya CJ, Howard DV, Howard JH. Age-related differences in multiple measures of white matter integrity: A diffusion tensor imaging study of healthy aging. Hum Brain Mapp (2010) 31:378–90. doi: 10.1002/hbm.20872
7. Bartzokis G, Lu PH, Heydari P, Couvrette A, Lee GJ, Kalashyan G, et al. Multimodal Magnetic Resonance Imaging Assessment of White Matter Aging Trajectories Over the Lifespan of Healthy Individuals. Biol Psychiatry (2012) 72:1026–34. doi: 10.1016/j.biopsych.2012.07.010
8. Cox SR, Ritchie SJ, Tucker-Drob EM, Liewald DC, Hagenaars SP, Davies G, et al. Ageing and brain white matter structure in 3,513 UK Biobank participants. Nat Commun (2016) 7:13629. doi: 10.1038/ncomms13629
10. Rathee R, Rallabandi VPS, Roy PK. Age-related differences in white matter integrity in healthy human brain: evidence from structural MRI and diffusion tensor imaging. Magn Reson Insights (2016) 9:9–20. doi: 10.4137/MRI.S39666
11. Westlye LT, Walhovd KB, Dale AM, Bjornerud A, Due-Tonnessen P, Engvig A, et al. Life-Span Changes of the Human Brain White Matter: Diffusion Tensor Imaging (DTI) and Volumetry. Cereb Cortex (2010) 20:2055–68. doi: 10.1093/cercor/bhp280
12. 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:15425–36. doi: 10.1523/JNEUROSCI.0203-14.2014
13. for the Alzheimer's Disease Neuroimaging Initiative, Chandra A, Dervenoulas G, Politis M. Magnetic resonance imaging in Alzheimer's disease and mild cognitive impairment. J Neurol (2019) 266:1293–302. doi: 10.1007/s00415-018-9016-3
15. Shen X, Adams MJ, Ritakari TE, Cox SR, McIntosh AM, Whalley HC. White matter microstructure and its relation to longitudinal measures of depressive symptoms in mid- and late life. Biol Psychiatry (2019) 86:759–68. doi: 10.1016/j.biopsych.2019.06.011
16. Suzuki H, Gao H, Bai W, Evangelou E, Glocker B, O'Regan DP, et al. Abnormal brain white matter microstructure is associated with both pre-hypertension and hypertension. PloS One (2017) 12:e0187600. doi: 10.1371/journal.pone.0187600
17. Croall ID, Lohner V, Moynihan B, Khan U, Hassan A, O'Brien JT, et al. Using DTI to assess white matter microstructure in cerebral small vessel disease (SVD) in multicentre studies. Clin Sci (2017) 131:1361–73. doi: 10.1042/CS20170146
19. Zavaliangos-Petropulu A, Nir TM, Thomopoulos SI, Reid RI, Bernstein MA, Borowski B, et al. Diffusion MRI Indices and their relation to cognitive impairment in brain aging: the updated multi-protocol approach in ADNI3. Front Neuroinform (2019) 13:2. doi: 10.3389/fninf.2019.00002
20. Baykara E, Gesierich B, Adam R, Tuladhar AM, Biesbroek JM, Koek HL, et al. a novel imaging marker for small vessel disease based on skeletonization of white matter tracts and diffusion histograms. Ann Neurol (2016) 80:581–92. doi: 10.1002/ana.24758
21. Wei N, Deng Y, Yao L, Jia W, Wang J, Shi Q, et al. A Neuroimaging marker based on diffusion tensor imaging and cognitive impairment due to cerebral white matter lesions. Front Neurol (2019) 10:81. doi: 10.3389/fneur.2019.00081
22. Deary IJ, Ritchie SJ, Muñoz Maniega S, Cox SR, Valdés Hernández MC, Luciano M, et al. Brain Peak Width of Skeletonized Mean Diffusivity (PSMD) and cognitive function in later life. Front Psychiatry (2019) 10:524. doi: 10.3389/fpsyt.2019.00524
23. Barrio-Arranz G, de Luis-García R, Tristán-Vega A, Martín-Fernández M, Aja-Fernández S. Impact of MR acquisition parameters on DTI scalar indexes: a tractography based approach. PloS One (2015) 10:e0137905. doi: 10.1371/journal.pone.0137905
24. Curran KM, Emsell L, Leemans A. Quantitative DTI Measures. In: Van Hecke W, Emsell L, Sunaert S, editors. Diffusion Tensor Imaging. New York, NY: Springer New York (2016). p. 65–87. doi: 10.1007/978-1-4939-3118-7_5
25. Sullivan EV, Adalsteinsson E, Hedehus M, Ju C, Moseley M, Lim KO, et al. Equivalent disruption of regional white matter microstructure in ageing healthy men and women. Neuroreport (2001) 12:99–104. doi: 10.1097/00001756-200101220-00027
26. Hasan KM, Kamali A, Abid H, Kramer LA, Fletcher JM, Ewing-Cobbs L. Quantification of the spatiotemporal microstructural organization of the human brain association, projection and commissural pathways across the lifespan using diffusion tensor tractography. Brain Struct Funct (2010) 214:361–73. doi: 10.1007/s00429-009-0238-0
27. Kochunov P, Glahn DC, Lancaster J, Thompson PM, Kochunov V, Rogers B, et al. Fractional anisotropy of cerebral white matter and thickness of cortical gray matter across the lifespan. NeuroImage (2011) 58:41–9. doi: 10.1016/j.neuroimage.2011.05.050
28. Hsu J-L, Leemans A, Bai C-H, Lee C-H, Tsai Y-F, Chiu H-C, et al. Gender differences and age-related white matter changes of the human brain: A diffusion tensor imaging study. NeuroImage (2008) 39:566–77. doi: 10.1016/j.neuroimage.2007.09.017
30. Helmer KG, Chou M-C, Preciado RI, Gimi B, Rollins NK, Song A, et al. Multi-site study of diffusion metric variability: characterizing the effects of site, vendor, field strength, and echo time using the histogram distance. In: Gimi B, Krol A, editors. Proc SPIE Int Soc Opt Eng. San Diego, California, United States: Int Soc Opt. Eng. (2016). p. 9788. doi: 10.1117/12.2217449
31. Qin W, Shui Yu C, Zhang F, Du XY, Jiang H, Xia Yan Y, et al. Effects of echo time on diffusion quantification of brain white matter at 1.5T and 3.0T. Magn Reson Med (2009) 61:755–60. doi: 10.1002/mrm.21920
32. Hutchinson EB, Avram AV, Irfanoglu MO, Koay CG, Barnett AS, Komlosh ME, et al. Analysis of the effects of noise, DWI sampling, and value of assumed parameters in diffusion MRI models: Cross-Model Analysis of Noise and DWI Sampling. Magn Reson Med (2017) 78:1767–80. doi: 10.1002/mrm.26575
33. Fortin J-P, Parker D, Tunç B, Watanabe T, Elliott MA, Ruparel K, et al. Harmonization of multi-site diffusion tensor imaging data. NeuroImage (2017) 161:149–70. doi: 10.1016/j.neuroimage.2017.08.047
34. Pohl KM, Sullivan EV, Rohlfing T, Chu W, Kwon D, Nichols BN, et al. Harmonizing DTI measurements across scanners to examine the development of white matter microstructure in 803 adolescents of the NCANDA study. NeuroImage (2016) 130:194–213. doi: 10.1016/j.neuroimage.2016.01.061
35. Grech-Sollars M, Hales PW, Miyazaki K, Raschke F, Rodriguez D, Wilson M, et al. Multi-centre reproducibility of diffusion MRI parameters for clinical sequences in the brain: Multi-centre reproducibility of diffusion MRI using clinical sequences. NMR BioMed (2015) 28:468–85. doi: 10.1002/nbm.3269
36. Prohl AK, Scherrer B, Tomas-Fernandez X, Filip-Dhima R, Kapur K, Velasco-Annis C, et al. Reproducibility of Structural and Diffusion Tensor Imaging in the TACERN Multi-Center Study. Front Integr Neurosci (2019) 13:24. doi: 10.3389/fnint.2019.00024
37. Jockwitz C, Mérillat S, Liem F, Oschwald J, Amunts K, Caspers S, et al. Generalizing age effects on brain structure and cognition: A two-study comparison approach. Hum Brain Mapp (2019) 40:2305–19. doi: 10.1002/hbm.24524
38. Jahanshad N, Kochunov PV, Sprooten E, Mandl RC, Nichols TE, Almasy L, et al. Multi-site genetic analysis of diffusion images and voxelwise heritability analysis: A pilot project of the ENIGMA–DTI working group. NeuroImage (2013) 81:455–69. doi: 10.1016/j.neuroimage.2013.04.061
39. Kochunov P, Jahanshad N, Sprooten E, Nichols TE, Mandl RC, Almasy L, et al. Multi-site study of additive genetic effects on fractional anisotropy of cerebral white matter: Comparing meta and megaanalytical approaches for data pooling. NeuroImage (2014) 95:136–50. doi: 10.1016/j.neuroimage.2014.03.033
40. Chad JA, Pasternak O, Salat DH, Chen JJ. Re-examining age-related differences in white matter microstructure with free-water corrected diffusion tensor imaging. Neurobiol Aging (2018) 71:161–70. doi: 10.1016/j.neurobiolaging.2018.07.018
41. Wu Y-C, Field AS, Whalen PJ, Alexander AL. Age- and gender-related changes in the normal human brain using hybrid diffusion imaging (HYDI). NeuroImage (2011) 54:1840–53. doi: 10.1016/j.neuroimage.2010.09.067
42. Fjell Anders M, Walhovd KB, Westlye LT, Østby Y, Tamnes CK, Jernigan TL, et al. When does brain aging accelerate? Dangers of quadratic fits in cross-sectional studies. NeuroImage (2010) 50:1376–83. doi: 10.1016/j.neuroimage.2010.01.061
43. Pfefferbaum A, Sullivan EV. Cross-sectional versus longitudinal estimates of age-related changes in the adult brain: overlaps and discrepancies. Neurobiol Aging (2015) 36:2563–7. doi: 10.1016/j.neurobiolaging.2015.05.005
44. Lebel C, Walker L, Leemans A, Phillips L, Beaulieu C. Microstructural maturation of the human brain from childhood to adulthood. NeuroImage (2008) 40:1044–55. doi: 10.1016/j.neuroimage.2007.12.053
46. Kodiweera C, Alexander AL, Harezlak J, McAllister TW, Wu Y-C. Age effects and sex differences in human brain white matter of young to middle-aged adults: A DTI, NODDI, and q-space study. NeuroImage (2016) 128:180–92. doi: 10.1016/j.neuroimage.2015.12.033
Keywords: ageing, white matter, neurodegeneration, MRI, diffusion, PSMD
Citation: Beaudet G, Tsuchida A, Petit L, Tzourio C, Caspers S, Schreiber J, Pausova Z, Patel Y, Paus T, Schmidt R, Pirpamer L, Sachdev PS, Brodaty H, Kochan N, Trollor J, Wen W, Armstrong NJ, Deary IJ, Bastin ME, Wardlaw JM, Munõz Maniega S, Witte AV, Villringer A, Duering M, Debette S and Mazoyer B (2020) Age-Related Changes of Peak Width Skeletonized Mean Diffusivity (PSMD) Across the Adult Lifespan: A Multi-Cohort Study. Front. Psychiatry 11:342. doi: 10.3389/fpsyt.2020.00342
Received: 06 January 2020; Accepted: 06 April 2020;
Published: 04 May 2020.
Edited by:Elisa Ambrosi, Hospital of the Brothers of St. John of God, Austria
Reviewed by:Arvind Caprihan, Mind Research Network (MRN), United States
Fabrizio Piras, Santa Lucia Foundation (IRCCS), Italy
Copyright © 2020 Beaudet, Tsuchida, Petit, Tzourio, Caspers, Schreiber, Pausova, Patel, Paus, Schmidt, Pirpamer, Sachdev, Brodaty, Kochan, Trollor, Wen, Armstrong, Deary, Bastin, Wardlaw, Munõz Maniega, Witte, Villringer, Duering, Debette and Mazoyer. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Bernard Mazoyer, email@example.com