Arrow Research search

Author name cluster

Frederik Barkhof

Possible papers associated with this exact author name in Arrow. This page groups case-insensitive exact name matches and is not a full identity disambiguation profile.

71 papers
1 author row

Possible papers

71

AAAI Conference 2026 Conference Paper

A Stage-Aware Mixture of Experts Framework for Neurodegenerative Disease Progression Modelling

  • Tiantian He
  • Keyue Jiang
  • An Zhao
  • Anna Schroder
  • Elinor Thompson
  • Sonja Soskic
  • Frederik Barkhof
  • Daniel C. Alexander

The long-term progression of neurodegenerative diseases is commonly conceptualized as a spatiotemporal diffusion process that consists of a graph diffusion process across the structural brain connectome and a localized reaction process within brain regions. However, modeling this progression remains challenging due to 1) the scarcity of longitudinal data obtained through irregular and infrequent subject visits and 2) the complex interplay of pathological mechanisms across brain regions and disease stages, where traditional models assume fixed mechanisms throughout disease progression. To address these limitations, we propose a novel stage-aware Mixture of Experts (MoE) framework that explicitly models how different contributing mechanisms dominate at different disease stages through time-dependent expert weighting. This architecture is a key innovation designed to maximize the utility of small datasets and provide interpretable insights into disease etiology. Data-wise, we utilize an iterative dual optimization method to properly estimate the temporal position of individual observations, constructing a cohort-level progression trajectory from irregular snapshots. Model-wise, we enhance the spatial component with an inhomogeneous graph neural diffusion model (IGND) that allows diffusivity to vary based on node states and time, providing more flexible representations of brain networks. We also introduce a localized neural reaction module to capture complex dynamics beyond standard processes.The resulting IGND-MoE model dynamically integrates these components across temporal states, offering a principled way to understand how stage-specific pathological mechanisms contribute to progression. When used to model tau pathology propagation in human brains, IGND-MoE outperforms purely pathophysiological and purely neural baselines in long-term prediction accuracy. Moreover, its stage-wise weights yield novel clinical insights that align with literature, suggesting that graph-related processes are more influential at early stages, while other unknown physical processes become dominant later on. Our findings highlight the necessity of designing hybrid and expert-constrained models that account for the evolving nature of neurodegenerative processes.

YNICL Journal 2026 Journal Article

Blood-brain barrier water exchange in relation to amyloid, cognition and cerebrovascular burden

  • Beatriz E. Padrela
  • Sandra Tecelão
  • Bjørn-Eivind Kirsebom
  • Oliver Geier
  • Mario Tranfa
  • Federico Masserini
  • Markus H. Sneve
  • Maksim Slivka

Blood-brain barrier (BBB) water exchange may serve as a sensitive early biomarker for Alzheimer's disease and age-related cognitive decline. This study applied a non-invasive multi-echo arterial spin labeling (ASL) technique to measure BBB water exchange time (Tex), cerebral blood flow (CBF), and arterial transit time (ATT) in 160 adults aged 50 years and older. Participants were classified as cognitively normal (CN), having subjective cognitive decline (SCD), or mild cognitive impairment (MCI). They were assessed for amyloid status and cerebrovascular burden. Compared to CN participants, Tex was significantly lower in both SCD (-9.5 %) and MCI (-14.5 %) groups, suggesting that reductions in BBB water exchange emerge early in the course of cognitive decline. In contrast, CBF was reduced only in MCI participants (-20.8 % compared to CN), and ATT was significantly increased only in individuals with severe cerebrovascular burden (Fazekas score 3). Notably, Tex showed a stepwise decrease with increasing Fazekas scores (1-2), supporting its sensitivity to moderate small vessel disease. No associations were found between Tex and amyloid positivity after adjusting for age and sex. These findings indicate that Tex alterations may precede changes in traditional perfusion markers and are more closely related to vascular and early cognitive changes than to amyloid pathology. BBB water exchange mapping may therefore provide a promising, non-invasive tool to detect early neurovascular dysfunction that contributes to cognitive decline in aging populations, potentially offering a useful biomarker for early intervention trials targeting vascular contributions to dementia.

YNIMG Journal 2026 Journal Article

Decreased functional connectivity in post-COVID syndrome patients with high neuroinflammatory activity

  • Denise Visser
  • Leo Pieperhoff
  • Luigi Lorenzini
  • Sandeep S.V. Golla
  • Sander C.J. Verfaillie
  • Anouk Verveen
  • Elsmarieke van de Giessen
  • Marijke E. den Hollander

INTRODUCTION: Previous research suggests some post-COVID patients with neurocognitive complaints (NCC) show neuroinflammation. Like in inflammatory diseases, this may affect network connectivity. This study aimed to compare resting-state functional connectivity in individuals with and without i) neuroinflammation and ii) persistent post-COVID NCC. METHODS: F]DPA-714 PET scan with arterial sampling for neuroinflammation, and 3T MRI scan for resting-state functional connectivity. Twenty resting-state networks (RSNs) were identified using independent component analysis. Group differences in within-RSN connectivity were analyzed using general linear models. Differences in subcortico-cortical between-RSN connectivity-between brainstem or thalamus and cortical RSNs-were assessed using interaction models. RESULTS: F]DPA-714 binding in participants with NCC. Lastly, thalamic-somatomotor and brainstem-control network connectivity (between-RSN) was altered in both individuals with persistent NCC and those with neuroinflammation, with thalamic-somatomotor changes mainly driven by NCC and brainstem-control changes by neuroinflammation. CONCLUSION: Our results suggest that neuroinflammation in individuals with persistent NCC after SARS-CoV-2 infection is linked to altered functional connectivity in RSNs central to higher-order cognitive functions.

YNICL Journal 2026 Journal Article

Longitudinal MRI-based changes in intracranial volume and skull thickness observed in both metachromatic leukodystrophy and multiple sclerosis

  • Guus H.J. Vorst
  • Nicole I. Wolf
  • David R. van Nederpelt
  • Frederik Barkhof
  • Marjo S. van der Knaap
  • Menno M. Schoonheim
  • Eva M.M. Strijbis
  • Petra J.W. Pouwels

Intracranial volume (ICV) is often used as normalization factor in volumetrics and considered to be stable in young adults. We noticed thick skulls on MRI scans of leukodystrophy patients, suggesting potentially changing skull morphology. In this study we aimed to quantify skull thickness and ICV in metachromatic leukodystrophy (MLD) patients and people with multiple sclerosis (pwMS). We retrospectively analyzed single-center cross-sectional and longitudinal MRI scans. Skull thickness and ICV were determined using automated segmentation techniques. Participants included MLD (n = 32, 11 male, scans = 136, median age first scan = 14.1 [IQR 7.9-25.7] years), MS (n = 232, 78 male, median age first scan = 47.3 [IQR 39.6-55.4] years, scans = 431), and controls (n = 139, 67 male, median age first scan = 30.7 [IQR 10.7-48.9] years, scans = 283). Both ICV and skull thickness showed natural growth in young controls. In young MLD participants, ICV decreased (-18.8 ± 22.4 mL/year, p < 0.001). Above age 20, ICV and skull thickness remained stable in controls. In comparison to controls, we observed ICV loss in MLD participants (-4.01 ± 8.29 mL/year, p < 0.001) and in pwMS (-2.99 ± 2.69 mL/year, p < 0.001), as well as skull thickening (MLD: 0.16 ± 0.14 mm/year, p < 0.001, pwMS: 0.04 ± 0.09 mm/year, p = 0.009). In adult patient groups, negative correlations were found between ICV and skull thickness (MLD: -19.24 mL/mm, p < 0.001, pwMS: -11.56 mL/mm, p < 0.001), but not in controls (p = 0.11). Despite limitations due to scanner variations, segmentation reliability and absence of validation against ground truth, these findings demonstrate a reduction in ICV in pathologies, already in young adulthood. Although the observed changes are small, they may lead to underestimations of atrophy when using ICV as a normalization factor.

YNIMG Journal 2025 Journal Article

A large-scale multi-centre study characterising atrophy heterogeneity in Alzheimer’s disease

  • Vikram Venkatraghavan
  • Damiano Archetti
  • Pierrick Bourgeat
  • Chenyang Jiang
  • Mara ten Kate
  • Anna C. van Loenhoud
  • Rik Ossenkoppele
  • Charlotte E. Teunissen

Previous studies identified atrophy-based Alzheimer's disease(AD) subtypes linked to distinct clinical symptoms, but their consistency across subtyping approaches remains unclear. This large-scale study evaluates subtype concordance using two data-driven approaches. In this work, we analyzed data from n=10,011 patients across 10 AD cohorts spanning Europe, the US, and Australia, extracting regional volumes using Freesurfer. To characterize atrophy heterogeneity in the AD continuum, we developed a two-step approach, Snowphlake (Staging NeurOdegeneration With PHenotype informed progression timeLine of biomarKErs), to identify subtypes and atrophy-event sequences within each subtype. Results were compared with SuStaIn (Subtype and Stage Inference), which jointly estimates subtypes and staging, using similar training and validation. Training included Aβ+ participants (n=1,195) and Aβ- cognitively unimpaired controls (n=1,692). We validated model-staging in a held-out clinical dataset (n=6,362) and an independent dataset (n=762), and assessed clinical significance in Aβ+ subsets(n=1,796 held-out; n=159 external). Concordance analysis evaluated consistency between methods. In the AD dementia(AD-D) training data, both Snowphlake and SuStaIn identified four subtypes. In the validation datasets, staging with both methods correlated with Mini-Mental State Examination(MMSE) scores. The Snowphlake subtypes assigned in Aβ+ validation datasets were associated with alterations in specific cognitive domains(Cohen's f:[0.15-0.33]). Similarly, the SuStaIn subtypes were also associated specific cognitive domains(Cohen's f:[0.17-0.34]). However, we observed low concordance between Snowphlake and SuStaIn, with 39.7% of AD-D patients grouped in concordant subtypes by both methods. In conclusion, Snowphlake and SuStaIn identified four atrophy-based subtypes that linked to distinct symptom profiles. While this highlights that the neuro-anatomically defined subtypes also meaningfully associate with different cognitive impairments at a group level, the low concordance between methods suggests that future research is needed to better understand the biological and methodological factors contributing to the observed variability.

YNICL Journal 2025 Journal Article

Amyloid PET predicts atrophy in older adults without dementia: Results from the AMYPAD Prognostic & Natural History study

  • Leonard Pieperhoff
  • Luigi Lorenzini
  • Sophie Mastenbroek
  • Mario Tranfa
  • Mahnaz Shekari
  • Alle Meije Wink
  • Robin Wolz
  • Sylke Grootoonk

Highlights • Aβ-PET predicts atrophy in key brain regions in older adults without dementia.• Fusiform volumetric loss is linked to Aβ independent of CSF tau levels.• Temporal atrophy is stronger in women with higher Aβ burden.• APOE-ε4 carriers show larger Aβ-driven frontal & hippocampal atrophy.

YNICL Journal 2025 Journal Article

Individualized lesion-symptom mapping using explainable artificial intelligence for the cognitive impact of white matter hyperintensities

  • Ryanne Offenberg
  • Alberto de Luca
  • Geert Jan Biessels
  • Frederik Barkhof
  • Wiesje M. van der Flier
  • Argonde C. van Harten
  • Ewoud van der Lelij
  • Josien Pluim

Lesion-symptom mapping methods assess the relationship between lesions caused by cerebral small vessel disease and cognition, but current technology like support vector regression (SVR)) primarily provide group-level results. We propose a novel lesion-symptom mapping approach that can indicate how lesion patterns contribute to cognitive impairment on an individual level. A convolutional neural network (CNN) predicts cognitive scores and is combined with explainable artificial intelligence (XAI) to map the relation between cognition and vascular lesions. This method was evaluated primarily using real white matter hyperintensity maps of 821 memory clinic patients and simulated cognitive data, with weighted lesions and noise levels. Simulated data provided ground truth locations to assess predictive performance of the CNN and accuracy of strategic lesion identification by XAI, using an established lesion-symptom mapping method, SVR, and a simple fully connected neural network (FNN) as benchmarks. Real cognitive scores were used in a final proof-of-principle analysis. Predictive performance in simulation experiments was high for the CNN (R2 = 0. 964), SVR (R2 = 0. 875), and FNN (R2 = 0. 863). CNN with XAI provided patient-specific attribution maps that highlighted the ground truth locations. All methods showed similar sensitivity to noise. Using real cognitive scores, SVR (R2 = 0. 291) obtained a somewhat higher predictive performance than the CNN (R2 = 0. 216), although both methods substantially exceeded the predictive performance of total WMH volume alone (R2 = 0. 013). The FNN performed worse on real data (R2 = 0. 020). To conclude, results show that CNNs combined with XAI can perform lesion-symptom mapping and generate individual attribution maps, which could be a valuable feature with further method development.

YNICL Journal 2025 Journal Article

Regional cerebellar atrophy related to disability and cognitive progression in multiple sclerosis

  • Myrte Strik
  • Iris Dekker
  • Aurélie Ruet
  • Hanneke E. Hulst
  • Mike P. Wattjes
  • Frederik Barkhof
  • Bernard M.J. Uitdehaag
  • Joep Killestein

OBJECTIVE: The implications of cerebellar pathology on clinical disease progression in multiple sclerosis (MS) remain unclear. This study investigated regional cerebellar atrophy related to physical disability and cognitive impairment progression. METHODS: We included 331 MS patients and 95 controls (Amsterdam MS Cohort, 229 patients and 58 controls re-evaluated after 5 years). Assessments included baseline MRI, and disability and cognition at baseline and follow up. Cerebellar (sub)cortex was parcellated, volumetric data were determined and related to baseline disability and cognition. Longitudinal progression was explored only for regions with significant baseline correlations. RESULTS: At baseline, patients had mild disability (median EDSS 3.0) and 46% showing mild-to-severe cognitive impairment. At follow-up, 34.5% showed EDSS progression and 26.6% cognitive decline. All global and most regional volumes showed atrophy. Cross-sectionally, atrophy of several regions encompassing both anterior and posterior lobes correlated with both disability and cognition, while some correlated with EDSS only. Additionally, cerebellar nuclei only correlated with cognition. Cerebellar volumes were mainly related to information processing speed, working and verbal memory. Longitudinally, atrophy in the posterior lobe, lobule VI and VIIIb, and vermis VI, correlated with cognitive decline, while no variables correlated with disability progression. CONCLUSION: Regional cerebellar atrophy in both anterior and posterior lobes correlated with disability and cognitive impairment. Posterior regional atrophy was correlated with longitudinal cognitive decline, but none correlated with disability progression. Further research is required to elucidate these relationships.

YNIMG Journal 2024 Journal Article

Deep learning enables accurate brain tissue microstructure analysis based on clinically feasible diffusion magnetic resonance imaging

  • Yuxing Li
  • Zhizheng Zhuo
  • Chenghao Liu
  • Yunyun Duan
  • Yulu Shi
  • Tingting Wang
  • Runzhi Li
  • Yanli Wang

Diffusion magnetic resonance imaging (dMRI) allows non-invasive assessment of brain tissue microstructure. Current model-based tissue microstructure reconstruction techniques require a large number of diffusion gradients, which is not clinically feasible due to imaging time constraints, and this has limited the use of tissue microstructure information in clinical settings. Recently, approaches based on deep learning (DL) have achieved promising tissue microstructure reconstruction results using clinically feasible dMRI. However, it remains unclear whether the subtle tissue changes associated with disease or age are properly preserved with DL approaches and whether DL reconstruction results can benefit clinical applications. Here, we provide the first evidence that DL approaches to tissue microstructure reconstruction yield reliable brain tissue microstructure analysis based on clinically feasible dMRI scans. Specifically, we reconstructed tissue microstructure from four different brain dMRI datasets with only 12 diffusion gradients, a clinically feasible protocol, and the neurite orientation dispersion and density imaging (NODDI) and spherical mean technique (SMT) models were considered. With these results we show that disease-related and age-dependent alterations of brain tissue were accurately identified. These findings demonstrate that DL tissue microstructure reconstruction can accurately quantify microstructural alterations in the brain based on clinically feasible dMRI.

YNICL Journal 2024 Journal Article

Scanner-specific optimisation of automated lesion segmentation in MS

  • David R. van Nederpelt
  • Giuseppe Pontillo
  • Mar Barrantes-Cepas
  • Iman Brouwer
  • Eva M.M. Strijbis
  • Menno M. Schoonheim
  • Bastiaan Moraal
  • Bas Jasperse

BACKGROUND & OBJECTIVE: Automatic lesion segmentation techniques on MRI scans of people with multiple sclerosis (pwMS) could support lesion detection and segmentation in trials and clinical practice. However, knowledge on their reliability across scanners is limited, hampering clinical implementation. The aim of this study was to investigate the within-scanner repeatability and between-scanner reproducibility of lesion segmentation tools in pwMS across three different scanners and examine their accuracy compared to manual segmentations with and without optimization. METHODS: 30 pwMS underwent a scan and rescan on three MRI scanners. GE Discovery MR750 (3.0 T), Siemens Sola (1.5 T) and Siemens Vida (3.0 T)). 3D-FLuid Attenuated Inversion Recovery (3D-FLAIR) and 3D T1-weighted scans were acquired on each scanner. Lesion segmentation involved preprocessing and automatic segmentation using the Lesion Segmentation Toolbox (LST) and nicMSlesions (nicMS) as well as manual segmentation. Both automated segmentation techniques were used with default settings, and with settings optimized to match manual segmentations for each scanner specifically and combined for the three scanners. LST settings were optimized by adjusting the threshold to improve the Dice similarity coefficient (DSC) for each scanner separately and a combined threshold for all scanners. For nicMS the last layers were retrained, once with the multi-scanner data to represent a combined optimization and once separately for each scanner for scanner specific optimization. Volumes and counts were extracted. DSC was calculated for accuracy, and reliability was assessed using intra-class correlation coefficients (ICC). Differences in DSC between software was tested with a repeated measures ANOVA and when appropriate post-hoc paired t-tests using Bonferroni correction. RESULTS: Scanner-specific optimization significantly improved DSC for LST compared to default and combined settings, except for the GE scanner. NicMS showed significantly higher DSC for both the scanner-specific and combined optimization than default. Within-scanner repeatability was excellent (ICC>0.9) for volume and counts. Between-scanner ICC for volume between Vida and Sola was higher (0.94-0.99) than between GE MR750 and Vida or Sola (0.18-0.93), with improved ICCs for nicMS scanner-specific (0.87-0.93) compared to others (0.18-0.79). This was not present for Sola vs. Vida where all ICCs were excellent (>0.94). CONCLUSION: Scanner-specific optimization strategies proved effective in mitigating inter-scanner variability, addressing the issue of insufficient reproducibility and accuracy found with default settings.

YNIMG Journal 2023 Journal Article

Evaluation of novel data-driven metrics of amyloid β deposition for longitudinal PET studies

  • Ariane Bollack
  • Pawel J Markiewicz
  • Alle Meije Wink
  • Lloyd Prosser
  • Johan Lilja
  • Pierrick Bourgeat
  • Jonathan M Schott
  • William Coath

PURPOSE: Positron emission tomography (PET) provides in vivo quantification of amyloid-β (Aβ) pathology. Established methods for assessing Aβ burden can be affected by physiological and technical factors. Novel, data-driven metrics have been developed to account for these sources of variability. We aimed to evaluate the performance of four of these amyloid PET metrics against conventional techniques, using a common set of criteria. METHODS: and CL, variability of the rate of change and sample size estimates to detect a 25% slowing in Aβ accumulation. RESULTS: and Aβ load compared to the CL. CONCLUSION: and Aβ load could offer a more precise alternative to CL, although further studies in larger cohorts should be conducted.

YNICL Journal 2023 Journal Article

Quantitative MRI in leukodystrophies

  • Menno D. Stellingwerff
  • Petra J.W. Pouwels
  • Stefan D. Roosendaal
  • Frederik Barkhof
  • Marjo S. van der Knaap

Leukodystrophies constitute a large and heterogeneous group of genetic diseases primarily affecting the white matter of the central nervous system. Different disorders target different white matter structural components. Leukodystrophies are most often progressive and fatal. In recent years, novel therapies are emerging and for an increasing number of leukodystrophies trials are being developed. Objective and quantitative metrics are needed to serve as outcome measures in trials. Quantitative MRI yields information on microstructural properties, such as myelin or axonal content and condition, and on the chemical composition of white matter, in a noninvasive fashion. By providing information on white matter microstructural involvement, quantitative MRI may contribute to the evaluation and monitoring of leukodystrophies. Many distinct MR techniques are available at different stages of development. While some are already clinically applicable, others are less far developed and have only or mainly been applied in healthy subjects. In this review, we explore the background, current status, potential and challenges of available quantitative MR techniques in the context of leukodystrophies.

YNICL Journal 2023 Journal Article

Spatial distributions of white matter hyperintensities on brain MRI: A pooled analysis of individual participant data from 11 memory clinic cohorts

  • Mirthe Coenen
  • Geert Jan Biessels
  • Charles DeCarli
  • Evan F. Fletcher
  • Pauline M. Maillard
  • Frederik Barkhof
  • Josephine Barnes
  • Thomas Benke

INTRODUCTION: The spatial distribution of white matter hyperintensities (WMH) on MRI is often considered in the diagnostic evaluation of patients with cognitive problems. In some patients, clinicians may classify WMH patterns as "unusual", but this is largely based on expert opinion, because detailed quantitative information about WMH distribution frequencies in a memory clinic setting is lacking. Here we report voxel wise 3D WMH distribution frequencies in a large multicenter dataset and also aimed to identify individuals with unusual WMH patterns. METHODS: Individual participant data (N = 3525, including 777 participants with subjective cognitive decline, 1389 participants with mild cognitive impairment and 1359 patients with dementia) from eleven memory clinic cohorts, recruited through the Meta VCI Map Consortium, were used. WMH segmentations were provided by participating centers or performed in Utrecht and registered to the Montreal Neurological Institute (MNI)-152 brain template for spatial normalization. To determine WMH distribution frequencies, we calculated WMH probability maps at voxel level. To identify individuals with unusual WMH patterns, region-of-interest (ROI) based WMH probability maps, rule-based scores, and a machine learning method (Local Outlier Factor (LOF)), were implemented. RESULTS: WMH occurred in 82% of voxels from the white matter template with large variation between subjects. Only a small proportion of the white matter (1.7%), mainly in the periventricular areas, was affected by WMH in at least 20% of participants. A large portion of the total white matter was affected infrequently. Nevertheless, 93.8% of individual participants had lesions in voxels that were affected in less than 2% of the population, mainly located in subcortical areas. Only the machine learning method effectively identified individuals with unusual patterns, in particular subjects with asymmetric WMH distribution or with WMH at relatively rarely affected locations despite common locations not being affected. DISCUSSION: Aggregating data from several memory clinic cohorts, we provide a detailed 3D map of WMH lesion distribution frequencies, that informs on common as well as rare localizations. The use of data-driven analysis with LOF can be used to identify unusual patterns, which might serve as an alert that rare causes of WMH should be considered.

YNICL Journal 2023 Journal Article

The spatio-temporal relationship between concurrent lesion and brain atrophy changes in early multiple sclerosis: A post-hoc analysis of the REFLEXION study

  • Giordano Gentile
  • Rozemarijn M. Mattiesing
  • Iman Brouwer
  • Ronald A. van Schijndel
  • Bernard M.J. Uitdehaag
  • Jos W.R. Twisk
  • Ludwig Kappos
  • Mark S. Freedman

BACKGROUND: White matter (WM) lesions and brain atrophy are present early in multiple sclerosis (MS). However, their spatio-temporal relationship remains unclear. METHODS: Yearly magnetic resonance images were analysed in 387 patients with a first clinical demyelinating event (FCDE) from the 5-year REFLEXION study. Patients received early (from baseline; N = 258; ET) or delayed treatment (from month-24; N = 129; DT) with subcutaneous interferon beta-1a. FSL-SIENA/VIENA were used to provide yearly percentage volume change of brain (PBVC) and ventricles (PVVC). Yearly total lesion volume change (TLVC) was determined by a semi-automated method. Using linear mixed models and voxel-wise analyses, we firstly investigated the overall relationship between TLVC and PBVC and between TLVC and PVVC in the same follow-up period. Analyses were then separately performed for: the untreated period of DT patients (first two years), the first year of treatment (year 1 for ET and year 3 for DT), and a period where patients had received at least 1 year of treatment (stable treatment; ET: years 2, 3, 4, and 5; DT: years 4 and 5). RESULTS: Whole brain: across the whole study period, lower TLVC was related to faster atrophy (PBVC: B = 0.046, SE = 0.013, p < 0.001; PVVC: B = -0.466, SE = 0.118, p < 0.001). Within the untreated period of DT patients, lower TLVC was related to faster atrophy (PBVC: B = 0.072, SE = 0.029, p = 0.013; PVVC: B = -0.917, SE = 0.306, p = 0.003). A similar relationship was found within the first year of treatment of ET patients (PBVC: B = 0.081, SE = 0.027, p = 0.003; PVVC: B = -1.08, SE = 0.284, p < 0.001), consistent with resolving oedema and pseudo-atrophy. Voxel-wise: overall, higher TLVC was related to faster ventricular enlargement. Lower TLVC was related to faster widespread atrophy in year 1 in both ET (first year of treatment) and DT (untreated) patients. In the second untreated year of DT patients and within the stable treatment period of ET patients (year 4), faster periventricular and occipital lobe atrophy was associated with higher TLVC. CONCLUSIONS: WM lesion changes and atrophy occurred simultaneously in early MS. Spatio-temporal correspondence of these two processes involved mostly the periventricular area. Within the first year of the study, in both treatment groups, faster atrophy was linked to lower lesion volume changes, consistent with higher shrinking and disappearing lesion activity. This might reflect the pseudo-atrophy phenomenon that is probably related to the therapy driven (only in ET patients, as they received treatment from baseline) and "natural" (both ET and DT patients entered the study after a FCDE) resolution of oedema. In an untreated period and later on during stable treatment, (real) atrophy was related to higher lesion volume changes, consistent with increased new and enlarging lesion activity.

YNIMG Journal 2023 Journal Article

Transferability of Alzheimer's disease progression subtypes to an independent population cohort

  • Hanyi Chen
  • Alexandra Young
  • Neil P. Oxtoby
  • Frederik Barkhof
  • Daniel C. Alexander
  • Andre Altmann

In the past, methods to subtype or biotype patients using brain imaging data have been developed. However, it is unclear whether and how these trained machine learning models can be successfully applied to population cohorts to study the genetic and lifestyle factors underpinning these subtypes. This work, using the Subtype and Stage Inference (SuStaIn) algorithm, examines the generalisability of data-driven Alzheimer's disease (AD) progression models. We first compared SuStaIn models trained separately on Alzheimer's disease neuroimaging initiative (ADNI) data and an AD-at-risk population constructed from the UK Biobank dataset. We further applied data harmonization techniques to remove cohort effects. Next, we built SuStaIn models on the harmonized datasets, which were then used to subtype and stage subjects in the other harmonized dataset. The first key finding is that three consistent atrophy subtypes were found in both datasets, which match the previously identified subtype progression patterns in AD: 'typical', 'cortical' and 'subcortical'. Next, the subtype agreement was further supported by high consistency in individuals' subtypes and stage assignment based on the different models: more than 92% of the subjects, with reliable subtype assignment in both ADNI and UK Biobank dataset, were assigned to an identical subtype under the model built on the different datasets. The successful transferability of AD atrophy progression subtypes across cohorts capturing different phases of disease development enabled further investigations of associations between AD atrophy subtypes and risk factors. Our study showed that (1) the average age is highest in the typical subtype and lowest in the subcortical subtype; (2) the typical subtype is associated with statistically more-AD-like cerebrospinal fluid biomarkers values in comparison to the other two subtypes; and (3) in comparison to the subcortical subtype, the cortical subtype subjects are more likely to associate with prescription of cholesterol and high blood pressure medications. In summary, we presented cross-cohort consistent recovery of AD atrophy subtypes, showing how the same subtypes arise even in cohorts capturing substantially different disease phases. Our study opened opportunities for future detailed investigations of atrophy subtypes with a broad range of early risk factors, which will potentially lead to a better understanding of the disease aetiology and the role of lifestyle and behaviour on AD.

YNICL Journal 2022 Journal Article

MAGNIMS recommendations for harmonization of MRI data in MS multicenter studies

  • Nicola De Stefano
  • Marco Battaglini
  • Deborah Pareto
  • Rosa Cortese
  • Jian Zhang
  • Niels Oesingmann
  • Ferran Prados
  • Maria A. Rocca

There is an increasing need of sharing harmonized data from large, cooperative studies as this is essential to develop new diagnostic and prognostic biomarkers. In the field of multiple sclerosis (MS), the issue has become of paramount importance due to the need to translate into the clinical setting some of the most recent MRI achievements. However, differences in MRI acquisition parameters, image analysis and data storage across sites, with their potential bias, represent a substantial constraint. This review focuses on the state of the art, recent technical advances, and desirable future developments of the harmonization of acquisition, analysis and storage of large-scale multicentre MRI data of MS cohorts. Huge efforts are currently being made to achieve all the requirements needed to provide harmonized MRI datasets in the MS field, as proper management of large imaging datasets is one of our greatest opportunities and challenges in the coming years. Recommendations based on these achievements will be provided here. Despite the advances that have been made, the complexity of these tasks requires further research by specialized academical centres, with dedicated technical and human resources. Such collective efforts involving different professional figures are of crucial importance to offer to MS patients a personalised management while minimizing consumption of resources.

YNICL Journal 2022 Journal Article

Slowly expanding lesions relate to persisting black-holes and clinical outcomes in relapse-onset multiple sclerosis

  • Alberto Calvi
  • Carmen Tur
  • Declan Chard
  • Jonathan Stutters
  • Olga Ciccarelli
  • Rosa Cortese
  • Marco Battaglini
  • Anna Pietroboni

BACKGROUND: Slowly expanding lesions (SELs) are MRI markers of chronic active lesions in multiple sclerosis (MS). T1-hypointense black holes, and reductions in magnetization transfer ratio (MTR) are pathologically correlated with myelin and axonal loss. While all associated with progressive MS, the relationship between these lesion's metrics and clinical outcomes in relapse-onset MS has not been widely investigated. OBJECTIVES: To explore the relationship of SELs with T1-hypointense black holes, and longitudinal T1 intensity contrast ratio and MTR, their correlation to brain volume, and their contribution to MS disability in relapse-onset patients. METHODS: 135 patients with relapsing-remitting MS (RRMS) were studied with clinical assessments and brain MRI (T2/FLAIR and T1-weighted scans at 1.5/3 T) at baseline and two subsequent follow-ups; a subset of 83 patients also had MTR acquisitions. Early-onset patients were defined when the baseline disease duration was ≤ 5 years (n = 85). SELs were identified using deformation field maps from the manually segmented baseline T2 lesions and differentiated from the non-SELs. Persisting black holes (PBHs) were defined as a subset of T2 lesions with a signal below a patient-specific grey matter T1 intensity in a semi-quantitative manner. SELs, PBH counts, and brain volume were computed, and their associations were assessed through Spearman and Pearson correlation. Clusters of patients according to low (up to 2), intermediate (3 to 10), or high (more than 10) SEL counts were determined with a Gaussian generalised mixture model. Mixed-effects and logistic regression models assessed volumes, T1 and MTR within SELs, and their correlation with Expanded Disability Status Scale (EDSS) and confirmed disability progression (CDP). RESULTS: Mean age at study onset was 35.5 years (73% female), disease duration 5.5 years and mean time to last follow-up 6.5 years (range 1 to 12.5); median baseline EDSS 1.5 (range 0 to 5.5) and a mean EDSS change of 0.31 units at final follow-up. Among 4007 T2 lesions, 27% were classified as SELs and 10% as PBHs. Most patients (n = 65) belonged to the cluster with an intermediate SEL count (3 to 10 SELs). The percentage of PBHs was higher in SELs than non-SELs (up to 61% vs 44%, p < 0.001) and within-patient SEL volumes positively correlated with PBH volumes (r = 0.53, p < 0.001). SELs showed a decrease in T1 intensity over time (beta = -0.004, 95%CI -0.005 to -0.003, p < 0.001), accompanied by lower cross-sectional baseline and follow-up MTR. In mixed-effects models, EDSS worsening was predicted by the SEL log-volumes increase over time (beta = 0.11, 95%CI 0.03 to 0.20, p = 0.01), which was confirmed in the sub-cohort of patients with early onset MS (beta = 0.14, 95%CI 0.04 to 0.25, p = 0.008). In logistic regressions, a higher risk for CDP was associated with SEL volumes (OR = 5.15, 95%CI 1.60 to 16.60, p = 0.006). CONCLUSIONS: SELs are associated with accumulation of more destructive pathology as indicated by an association with PBH volume, longitudinal reduction in T1 intensity and MTR. Higher SEL volumes are associated with clinical progression, while lower ones are associated with stability in relapse-onset MS.

YNICL Journal 2022 Journal Article

The Open-Access European Prevention of Alzheimer’s Dementia (EPAD) MRI dataset and processing workflow

  • Luigi Lorenzini
  • Silvia Ingala
  • Alle Meije Wink
  • Joost P.A. Kuijer
  • Viktor Wottschel
  • Mathijs Dijsselhof
  • Carole H. Sudre
  • Sven Haller

The European Prevention of Alzheimer Dementia (EPAD) is a multi-center study that aims to characterize the preclinical and prodromal stages of Alzheimer's Disease. The EPAD imaging dataset includes core (3D T1w, 3D FLAIR) and advanced (ASL, diffusion MRI, and resting-state fMRI) MRI sequences. Here, we give an overview of the semi-automatic multimodal and multisite pipeline that we developed to curate, preprocess, quality control (QC), and compute image-derived phenotypes (IDPs) from the EPAD MRI dataset. This pipeline harmonizes DICOM data structure across sites and performs standardized MRI preprocessing steps. A semi-automated MRI QC procedure was implemented to visualize and flag MRI images next to site-specific distributions of QC features - i.e. metrics that represent image quality. The value of each of these QC features was evaluated through comparison with visual assessment and step-wise parameter selection based on logistic regression. IDPs were computed from 5 different MRI modalities and their sanity and potential clinical relevance were ascertained by assessing their relationship with biological markers of aging and dementia. The EPAD v1500.0 data release encompassed core structural scans from 1356 participants 842 fMRI, 831 dMRI, and 858 ASL scans. From 1356 3D T1w images, we identified 17 images with poor quality and 61 with moderate quality. Five QC features - Signal to Noise Ratio (SNR), Contrast to Noise Ratio (CNR), Coefficient of Joint Variation (CJV), Foreground-Background energy Ratio (FBER), and Image Quality Rate (IQR) - were selected as the most informative on image quality by comparison with visual assessment. The multimodal IDPs showed greater impairment in associations with age and dementia biomarkers, demonstrating the potential of the dataset for future clinical analyses.

YNICL Journal 2022 Journal Article

The spatio-temporal relationship between white matter lesion volume changes and brain atrophy in clinically isolated syndrome and early multiple sclerosis

  • Rozemarijn M. Mattiesing
  • Giordano Gentile
  • Iman Brouwer
  • Ronald A. van Schijndel
  • Bernard M.J. Uitdehaag
  • Jos W.R. Twisk
  • Ludwig Kappos
  • Mark S. Freedman

BACKGROUND: White matter lesions and brain atrophy are both present early in multiple sclerosis. However, the spatio-temporal relationship between atrophy and lesion processes remains unclear. METHODS: Yearly magnetic resonance images were analyzed in 392 patients with clinically isolated syndrome from the 5-year REFLEX/REFLEXION studies. Patients received early treatment (from baseline; N = 262) or delayed treatment (from month-24; N = 130) with subcutaneous interferon beta-1a. Global and central atrophy were assessed using FSL-SIENA to provide yearly percentage volume change of brain and ventricles, respectively. Yearly total lesion volume change was calculated by subtracting the sum of the negative lesion volume change (disappearing + shrinking) from the positive lesion volume change (new + enlarging) for each yearly interval, as determined by an in-house developed semi-automated method. Using linear mixed models, during the period where patients had received ≥1 year of treatment, we investigated whether total lesion volume change was associated with percentage brain volume change or percentage ventricular volume change in the next year, and vice versa. RESULTS: Higher total lesion volume change was related to significantly faster global atrophy (percentage brain volume change) in the next year (B = - 0.113, SE = 0.022, p < 0.001). In patients receiving early treatment only, total lesion volume change was also associated with percentage ventricular volume change in the next year (B = 1.348, SE = 0.181, p < 0.001). Voxel-wise analyses showed that in patients receiving early treatment, higher total lesion volume change in years 2, 3, and 4 was related to faster atrophy in the next year, and in year 4 this relationship was stronger in patients receiving delayed treatment. Interestingly, faster atrophy was related to higher total lesion volume change in the next year (percentage brain volume change: B = - 0.136, SE = 0.062, p = 0.028; percentage ventricular volume change: B = 0.028, SE = 0.008, p < 0.001). CONCLUSIONS: Higher lesion volume changes were associated with faster atrophy in the next year. Interestingly, there was also an association between faster atrophy and higher lesion volume changes in the next year.

YNICL Journal 2021 Journal Article

Development and evaluation of a manual segmentation protocol for deep grey matter in multiple sclerosis: Towards accelerated semi-automated references

  • Alexandra de Sitter
  • Jessica Burggraaff
  • Fabian Bartel
  • Miklos Palotai
  • Yaou Liu
  • Jorge Simoes
  • Serena Ruggieri
  • Katharina Schregel

BACKGROUND: Deep grey matter (dGM) structures, particularly the thalamus, are clinically relevant in multiple sclerosis (MS). However, segmentation of dGM in MS is challenging; labeled MS-specific reference sets are needed for objective evaluation and training of new methods. OBJECTIVES: This study aimed to (i) create a standardized protocol for manual delineations of dGM; (ii) evaluate the reliability of the protocol with multiple raters; and (iii) evaluate the accuracy of a fast-semi-automated segmentation approach (FASTSURF). METHODS: A standardized manual segmentation protocol for caudate nucleus, putamen, and thalamus was created, and applied by three raters on multi-center 3D T1-weighted MRI scans of 23 MS patients and 12 controls. Intra- and inter-rater agreement was assessed through intra-class correlation coefficient (ICC); spatial overlap through Jaccard Index (JI) and generalized conformity index (CIgen). From sparse delineations, FASTSURF reconstructed full segmentations; accuracy was assessed both volumetrically and spatially. RESULTS: All structures showed excellent agreement on expert manual outlines: intra-rater JI > 0.83; inter-rater ICC ≥ 0.76 and CIgen ≥ 0.74. FASTSURF reproduced manual references excellently, with ICC ≥ 0.97 and JI ≥ 0.92. CONCLUSIONS: The manual dGM segmentation protocol showed excellent reproducibility within and between raters. Moreover, combined with FASTSURF a reliable reference set of dGM segmentations can be produced with lower workload.

YNICL Journal 2021 Journal Article

Differential patterns of gray matter volumes and associated gene expression profiles in cognitively-defined Alzheimer’s disease subgroups

  • Colin Groot
  • Michel J. Grothe
  • Shubhabrata Mukherjee
  • Irina Jelistratova
  • Iris Jansen
  • Anna Catharina van Loenhoud
  • Shannon L. Risacher
  • Andrew J. Saykin

The clinical presentation of Alzheimer's disease (AD) varies widely across individuals but the neurobiological mechanisms underlying this heterogeneity are largely unknown. Here, we compared regional gray matter (GM) volumes and associated gene expression profiles between cognitively-defined subgroups of amyloid-β positive individuals clinically diagnosed with AD dementia (age: 66 ± 7, 47% male, MMSE: 21 ± 5). All participants underwent neuropsychological assessment with tests covering memory, executive-functioning, language and visuospatial-functioning domains. Subgroup classification was achieved using a psychometric framework that assesses which cognitive domain shows substantial relative impairment compared to the intra-individual average across domains, which yielded the following subgroups in our sample; AD-Memory (n = 41), AD-Executive (n = 117), AD-Language (n = 33), AD-Visuospatial (n = 171). We performed voxel-wise contrasts of GM volumes derived from 3Tesla structural MRI between subgroups and controls (n = 127, age 58 ± 9, 42% male, MMSE 29 ± 1), and observed that differences in regional GM volumes compared to controls closely matched the respective cognitive profiles. Specifically, we detected lower medial temporal lobe GM volumes in AD-Memory, lower fronto-parietal GM volumes in AD-Executive, asymmetric GM volumes in the temporal lobe (left < right) in AD-Language, and lower GM volumes in posterior areas in AD-Visuospatial. In order to examine possible biological drivers of these differences in regional GM volumes, we correlated subgroup-specific regional GM volumes to brain-wide gene expression profiles based on a stereotactic characterization of the transcriptional architecture of the human brain as provided by the Allen human brain atlas. Gene-set enrichment analyses revealed that variations in regional expression of genes involved in processes like mitochondrial respiration and metabolism of proteins were associated with patterns of regional GM volume across multiple subgroups. Other gene expression vs GM volume-associations were only detected in particular subgroups, e.g., genes involved in the cell cycle for AD-Memory, specific sets of genes related to protein metabolism in AD-Language, and genes associated with modification of gene expression in AD-Visuospatial. We conclude that cognitively-defined AD subgroups show neurobiological differences, and distinct biological pathways may be involved in the emergence of these differences.

YNICL Journal 2021 Journal Article

Differential trajectories of hypometabolism across cognitively-defined Alzheimer’s disease subgroups

  • Colin Groot
  • Shannon L. Risacher
  • J.Q. Alida Chen
  • Ellen Dicks
  • Andrew J. Saykin
  • Christine L. Mac Donald
  • Jesse Mez
  • Emily H. Trittschuh

F]FDG-PET standardized uptake value ratios (SUVRs) to map hypometabolism across cognitively-defined AD subgroups. Participants were 384 amyloid-positive individuals with an AD dementia diagnosis from ADNI who had a total of 1028 FDG-scans (mean time between first and last scan: 1.6 ± 1.8 years). These participants were categorized into subgroups on the basis of substantial impairment at time of dementia diagnosis in a specific cognitive domain relative to the average across domains. This approach resulted in groups of AD-Memory (n = 135), AD-Executive (n = 8), AD-Language (n = 22), AD-Visuospatial (n = 44), AD-Multiple Domains (n = 15) and AD-No Domains (for whom no domain showed substantial relative impairment; n = 160). Voxelwise contrasts against controls revealed that all AD-subgroups showed progressive hypometabolism compared to controls across temporoparietal regions at time of AD diagnosis. Voxelwise and regions-of-interest (ROI)-based linear mixed model analyses revealed there were also subgroup-specific hypometabolism patterns and trajectories. The AD-Memory group had more pronounced hypometabolism compared to all other groups in the medial temporal lobe and posterior cingulate, and faster decline in metabolism in the medial temporal lobe compared to AD-Visuospatial. The AD-Language group had pronounced lateral temporal hypometabolism compared to all other groups, and the pattern of metabolism was also more asymmetrical (left < right) than all other groups. The AD-Visuospatial group had faster decline in metabolism in parietal regions compared to all other groups, as well as faster decline in the precuneus compared to AD-Memory and AD-No Domains. Taken together, in addition to a common pattern, cognitively-defined subgroups of people with AD dementia show subgroup-specific hypometabolism patterns, as well as differences in trajectories of metabolism over time. These findings provide support to the notion that cognitively-defined subgroups are biologically distinct.

YNICL Journal 2021 Journal Article

Manual and automated tissue segmentation confirm the impact of thalamus atrophy on cognition in multiple sclerosis: A multicenter study

  • Jessica Burggraaff
  • Yao Liu
  • Juan C. Prieto
  • Jorge Simoes
  • Alexandra de Sitter
  • Serena Ruggieri
  • Iman Brouwer
  • Birgit I. Lissenberg-Witte

BACKGROUND AND RATIONALE: Thalamus atrophy has been linked to cognitive decline in multiple sclerosis (MS) using various segmentation methods. We investigated the consistency of the association between thalamus volume and cognition in MS for two common automated segmentation approaches, as well as fully manual outlining. METHODS: Standardized neuropsychological assessment and 3-Tesla 3D-T1-weighted brain MRI were collected (multi-center) from 57 MS patients and 17 healthy controls. Thalamus segmentations were generated manually and using five automated methods. Agreement between the algorithms and manual outlines was assessed with Bland-Altman plots; linear regression assessed the presence of proportional bias. The effect of segmentation method on the separation of cognitively impaired (CI) and preserved (CP) patients was investigated through Generalized Estimating Equations; associations with cognitive measures were investigated using linear mixed models, for each method and vendor. RESULTS: In smaller thalami, automated methods systematically overestimated volumes compared to manual segmentations [ρ=(-0.42)-(-0.76); p-values < 0.001). All methods significantly distinguished CI from CP MS patients, except manual outlines of the left thalamus (p = 0.23). Poorer global neuropsychological test performance was significantly associated with smaller thalamus volumes bilaterally using all methods. Vendor significantly affected the findings. CONCLUSION: Automated and manual thalamus segmentation consistently demonstrated an association between thalamus atrophy and cognitive impairment in MS. However, a proportional bias in smaller thalami and choice of MRI acquisition system might impact the effect size of these findings.

YNIMG Journal 2021 Journal Article

Parametric imaging of dual-time window [18F]flutemetamol and [18F]florbetaben studies

  • Fiona Heeman
  • Maqsood Yaqub
  • Janine Hendriks
  • Ilona Bader
  • Frederik Barkhof
  • Juan Domingo Gispert
  • Bart N.M. van Berckel
  • Isadora Lopes Alves

Optimal pharmacokinetic models for quantifying amyloid beta (Aβ) burden using both [18F]flutemetamol and [18F]florbetaben scans have previously been identified at a region of interest (ROI) level. The purpose of this study was to determine optimal quantitative methods for parametric analyses of [18F]flutemetamol and [18F]florbetaben scans. Forty-six participants were scanned on a PET/MR scanner using a dual-time window protocol and either [18F]flutemetamol (N=24) or [18F]florbetaben (N=22). The following parametric approaches were used to derive DVR estimates: reference Logan (RLogan), receptor parametric mapping (RPM), two-step simplified reference tissue model (SRTM2) and multilinear reference tissue models (MRTM0, MRTM1, MRTM2), all with cerebellar grey matter as reference tissue. In addition, a standardized uptake value ratio (SUVR) was calculated for the 90–110 min post injection interval. All parametric images were assessed visually. Regional outcome measures were compared with those from a validated ROI method, i. e. DVR derived using RLogan. Visually, RPM, and SRTM2 performed best across tracers and, in addition to SUVR, provided highest AUC values for differentiating between Aβ-positive vs Aβ-negative scans ([18F]flutemetamol: range AUC=0. 96–0. 97 [18F]florbetaben: range AUC=0. 83–0. 85). Outcome parameters of most methods were highly correlated with the reference method (R 2≥0. 87), while lowest correlation were observed for MRTM2 (R 2=0. 71–0. 80). Furthermore, bias was low (≤5%) and independent of underlying amyloid burden for MRTM0 and MRTM1. The optimal parametric method differed per evaluated aspect; however, the best compromise across aspects was found for MRTM0 followed by SRTM2, for both tracers. SRTM2 is the preferred method for parametric imaging because, in addition to its good performance, it has the advantage of providing a measure of relative perfusion (R 1), which is useful for measuring disease progression.

YNICL Journal 2021 Journal Article

The sequence of structural, functional and cognitive changes in multiple sclerosis

  • Iris Dekker
  • Menno M. Schoonheim
  • Vikram Venkatraghavan
  • Anand J.C. Eijlers
  • Iman Brouwer
  • Esther E. Bron
  • Stefan Klein
  • Mike P. Wattjes

BACKGROUND: As disease progression remains poorly understood in multiple sclerosis (MS), we aim to investigate the sequence in which different disease milestones occur using a novel data-driven approach. METHODS: We analysed a cohort of 295 relapse-onset MS patients and 96 healthy controls, and considered 28 features, capturing information on T2-lesion load, regional brain and spinal cord volumes, resting-state functional centrality ("hubness"), microstructural tissue integrity of major white matter (WM) tracts and performance on multiple cognitive tests. We used a discriminative event-based model to estimate the sequence of biomarker abnormality in MS progression in general, as well as specific models for worsening physical disability and cognitive impairment. RESULTS: We demonstrated that grey matter (GM) atrophy of the cerebellum, thalamus, and changes in corticospinal tracts are early events in MS pathology, whereas other WM tracts as well as the cognitive domains of working memory, attention, and executive function are consistently late events. The models for disability and cognition show early functional changes of the default-mode network and earlier changes in spinal cord volume compared to the general MS population. Overall, GM atrophy seems crucial due to its early involvement in the disease course, whereas WM tract integrity appears to be affected relatively late despite the early onset of WM lesions. CONCLUSION: Data-driven modelling revealed the relative occurrence of both imaging and non-imaging events as MS progresses, providing insights into disease propagation mechanisms, and allowing fine-grained staging of patients for monitoring purposes.

YNIMG Journal 2021 Journal Article

Uncertainty analysis of MR-PET image registration for precision neuro-PET imaging

  • Pawel J. Markiewicz
  • Julian C. Matthews
  • John Ashburner
  • David M. Cash
  • David L. Thomas
  • Enrico De Vita
  • Anna Barnes
  • M. Jorge Cardoso

Accurate regional brain quantitative PET measurements, particularly when using partial volume correction, rely on robust image registration between PET and MR images. We argue here that the precision, and hence the uncertainty, of MR-PET image registration is mainly driven by the registration implementation and the quality of PET images due to their lower resolution and higher noise compared to the structural MR images. We propose a dedicated uncertainty analysis for quantifying the precision of MR-PET registration, centred around the bootstrap resampling of PET list-mode events to generate multiple PET image realisations with different noise (count) levels. The effects of PET image reconstruction parameters, such as the use of attenuation and scatter corrections and different number of iterations, on the precision and accuracy of MR-PET registration were investigated. In addition, the performance of four software packages with their default settings for rigid inter-modality image registration were considered: NiftyReg, Vinci, FSL and SPM. Four distinct PET image distributions made of two early time frames (similar to cortical FDG) and two late frames using two amyloid PET dynamic acquisitions of one amyloid positive and one amyloid negative participants were investigated. For the investigated four PET frames, the biggest impact on the uncertainty was observed between registration software packages (up to 10-fold difference in precision) followed by the reconstruction parameters. On average, the lowest uncertainty for different PET frames and brain regions was observed with SPM and two iterations of fully quantitative image reconstruction. The observed uncertainty for the varying PET count-level (from 5% to 60%) was slightly lower than for the reconstruction parameters. We also observed that the registration uncertainty in quantitative PET analysis depends on amyloid status of the considered PET frames, with increased uncertainty (up to three times) when using post-reconstruction partial volume correction. This analysis is applicable for PET data obtained from either PET/MR or PET/CT scanners.

YNIMG Journal 2020 Journal Article

Amygdalar nuclei and hippocampal subfields on MRI: Test-retest reliability of automated volumetry across different MRI sites and vendors

  • Giulia Quattrini
  • Michela Pievani
  • Jorge Jovicich
  • Marco Aiello
  • Núria Bargalló
  • Frederik Barkhof
  • David Bartres-Faz
  • Alberto Beltramello

BACKGROUND: The amygdala and the hippocampus are two limbic structures that play a critical role in cognition and behavior, however their manual segmentation and that of their smaller nuclei/subfields in multicenter datasets is time consuming and difficult due to the low contrast of standard MRI. Here, we assessed the reliability of the automated segmentation of amygdalar nuclei and hippocampal subfields across sites and vendors using FreeSurfer in two independent cohorts of older and younger healthy adults. METHODS: Sixty-five healthy older (cohort 1) and 68 younger subjects (cohort 2), from the PharmaCog and CoRR consortia, underwent repeated 3D-T1 MRI (interval 1-90 days). Segmentation was performed using FreeSurfer v6.0. Reliability was assessed using volume reproducibility error (ε) and spatial overlapping coefficient (DICE) between test and retest session. RESULTS: (for hippocampal subfields, except for molecular layer) had the best test-retest reproducibility (ε ​ ​0.80). CONCLUSION: Our results support the use of volumetric measures of larger amygdalar nuclei and hippocampal subfields in multisite MRI studies. These measures could be useful for disease tracking and assessment of efficacy in drug trials.

NeurIPS Conference 2020 Conference Paper

Disentangling Human Error from Ground Truth in Segmentation of Medical Images

  • Le Zhang
  • Ryutaro Tanno
  • Mou-Cheng Xu
  • Chen Jin
  • Joseph Jacob
  • Olga Cicarrelli
  • Frederik Barkhof
  • Daniel Alexander

Recent years have seen increasing use of supervised learning methods for segmentation tasks. However, the predictive performance of these algorithms depends on the quality of labels. This problem is particularly pertinent in the medical image domain, where both the annotation cost and inter-observer variability are high. In a typical label acquisition process, different human experts provide their estimates of the ``true'' segmentation labels under the influence of their own biases and competence levels. Treating these noisy labels blindly as the ground truth limits the performance that automatic segmentation algorithms can achieve. In this work, we present a method for jointly learning, from purely noisy observations alone, the reliability of individual annotators and the true segmentation label distributions, using two coupled CNNs. The separation of the two is achieved by encouraging the estimated annotators to be maximally unreliable while achieving high fidelity with the noisy training data. We first define a toy segmentation dataset based on MNIST and study the properties of the proposed algorithm. We then demonstrate the utility of the method on three public medical imaging segmentation datasets with simulated (when necessary) and real diverse annotations: 1) MSLSC (multiple-sclerosis lesions); 2) BraTS (brain tumours); 3) LIDC-IDRI (lung abnormalities). In all cases, our method outperforms competing methods and relevant baselines particularly in cases where the number of annotations is small and the amount of disagreement is large. The experiments also show strong ability to capture the complex spatial characteristics of annotators' mistakes. Our code is available at \url{https: //github. com/moucheng2017/Learn Noisy Labels Medical Images}.

YNIMG Journal 2020 Journal Article

ExploreASL: An image processing pipeline for multi-center ASL perfusion MRI studies

  • Henk J.M.M. Mutsaerts
  • Jan Petr
  • Paul Groot
  • Pieter Vandemaele
  • Silvia Ingala
  • Andrew D. Robertson
  • Lena Václavů
  • Inge Groote

Arterial spin labeling (ASL) has undergone significant development since its inception, with a focus on improving standardization and reproducibility of its acquisition and quantification. In a community-wide effort towards robust and reproducible clinical ASL image processing, we developed the software package ExploreASL, allowing standardized analyses across centers and scanners. The procedures used in ExploreASL capitalize on published image processing advancements and address the challenges of multi-center datasets with scanner-specific processing and artifact reduction to limit patient exclusion. ExploreASL is self-contained, written in MATLAB and based on Statistical Parameter Mapping (SPM) and runs on multiple operating systems. To facilitate collaboration and data-exchange, the toolbox follows several standards and recommendations for data structure, provenance, and best analysis practice. ExploreASL was iteratively refined and tested in the analysis of >10, 000 ASL scans using different pulse-sequences in a variety of clinical populations, resulting in four processing modules: Import, Structural, ASL, and Population that perform tasks, respectively, for data curation, structural and ASL image processing and quality control, and finally preparing the results for statistical analyses on both single-subject and group level. We illustrate ExploreASL processing results from three cohorts: perinatally HIV-infected children, healthy adults, and elderly at risk for neurodegenerative disease. We show the reproducibility for each cohort when processed at different centers with different operating systems and MATLAB versions, and its effects on the quantification of gray matter cerebral blood flow. ExploreASL facilitates the standardization of image processing and quality control, allowing the pooling of cohorts which may increase statistical power and discover between-group perfusion differences. Ultimately, this workflow may advance ASL for wider adoption in clinical studies, trials, and practice.

YNIMG Journal 2020 Journal Article

Generalised boundary shift integral for longitudinal assessment of spinal cord atrophy

  • Ferran Prados
  • Marcello Moccia
  • Aubrey Johnson
  • Marios Yiannakas
  • Francesco Grussu
  • Manuel Jorge Cardoso
  • Olga Ciccarelli
  • Sebastien Ourselin

Spinal cord atrophy measurements obtained from structural magnetic resonance imaging (MRI) are associated with disability in many neurological diseases and serve as in vivo biomarkers of neurodegeneration. Longitudinal spinal cord atrophy rate is commonly determined from the numerical difference between two volumes (based on 3D surface fitting) or two cross-sectional areas (CSA, based on 2D edge detection) obtained at different time-points. Being an indirect measure, atrophy rates are susceptible to variable segmentation errors at the edge of the spinal cord. To overcome those limitations, we developed a new registration-based pipeline that measures atrophy rates directly. We based our approach on the generalised boundary shift integral (GBSI) method, which registers 2 scans and uses a probabilistic XOR mask over the edge of the spinal cord, thereby measuring atrophy more accurately than segmentation-based techniques. Using a large cohort of longitudinal spinal cord images (610 subjects with multiple sclerosis from a multi-centre trial and 52 healthy controls), we demonstrated that GBSI is a sensitive, quantitative and objective measure of longitudinal spinal cord volume change. The GBSI pipeline is repeatable, reproducible, and provides more precise measurements of longitudinal spinal cord atrophy than segmentation-based methods in longitudinal spinal cord atrophy studies.

YNICL Journal 2020 Journal Article

Spinal cord atrophy in a primary progressive multiple sclerosis trial: Improved sample size using GBSI

  • Marcello Moccia
  • Nicola Valsecchi
  • Olga Ciccarelli
  • Ronald Van Schijndel
  • Frederik Barkhof
  • Ferran Prados

BACKGROUND: We aimed to evaluate the implications for clinical trial design of the generalised boundary-shift integral (GBSI) for spinal cord atrophy measurement. METHODS: ), acquired separately. We obtained segmentation-based cross-sectional spinal cord area (CSA) at C1-2 (from both brain and spinal cord MRI) and C2-5 levels (from spinal cord MRI) using DeepSeg, and, then, we computed corresponding GBSI. RESULTS: Depending on the spinal cord segment, we included 67.4-98.1% patients for CSA measurements, and 66.9-84.2% for GBSI. Spinal cord atrophy measurements obtained with GBSI had lower measurement variability, than corresponding CSA. Looking at the image noise floor, the lowest median standard deviation of the MRI signal within the cerebrospinal fluid surrounding the spinal cord was found on brain MRI at the C1-2 level. Spinal cord atrophy derived from brain MRI was related to the corresponding measures from dedicated spinal cord MRI, more strongly for GBSI than CSA. Spinal cord atrophy measurements using GBSI, but not CSA, were associated with upper and lower limb motor progression. DISCUSSION: Notwithstanding the reduced measurement variability, the clinical correlates, and the possibility of using brain acquisitions, spinal cord atrophy using GBSI should remain a secondary outcome measure in MS studies, until further advancements increase the quality of acquisition and reliability of processing.

YNICL Journal 2019 Journal Article

Automated fiber tract reconstruction for surgery planning: Extensive validation in language-related white matter tracts

  • Matteo Mancini
  • Sjoerd B. Vos
  • Vejay N. Vakharia
  • Aidan G. O'Keeffe
  • Karin Trimmel
  • Frederik Barkhof
  • Christian Dorfer
  • Salil Soman

Diffusion MRI and tractography hold great potential for surgery planning, especially to preserve eloquent white matter during resections. However, fiber tract reconstruction requires an expert with detailed understanding of neuroanatomy. Several automated approaches have been proposed, using different strategies to reconstruct the white matter tracts in a supervised fashion. However, validation is often limited to comparison with manual delineation by overlap-based measures, which is limited in characterizing morphological and topological differences. In this work, we set up a fully automated pipeline based on anatomical criteria that does not require manual intervention, taking advantage of atlas-based criteria and advanced acquisition protocols available on clinical-grade MRI scanners. Then, we extensively validated it on epilepsy patients with specific focus on language-related bundles. The validation procedure encompasses different approaches, including simple overlap with manual segmentations from two experts, feasibility ratings from external multiple clinical raters and relation with task-based functional MRI. Overall, our results demonstrate good quantitative agreement between automated and manual segmentation, in most cases better performances of the proposed method in qualitative terms, and meaningful relationships with task-based fMRI. In addition, we observed significant differences between experts in terms of both manual segmentation and external ratings. These results offer important insights on how different levels of validation complement each other, supporting the idea that overlap-based measures, although quantitative, do not offer a full perspective on the similarities and differences between automated and manual methods.

YNICL Journal 2019 Journal Article

Detecting frontotemporal dementia syndromes using MRI biomarkers

  • Marie Bruun
  • Juha Koikkalainen
  • Hanneke F.M. Rhodius-Meester
  • Marta Baroni
  • Le Gjerum
  • Mark van Gils
  • Hilkka Soininen
  • Anne M. Remes

BACKGROUND: Diagnosing frontotemporal dementia may be challenging. New methods for analysis of regional brain atrophy patterns on magnetic resonance imaging (MRI) could add to the diagnostic assessment. Therefore, we aimed to develop automated imaging biomarkers for differentiating frontotemporal dementia subtypes from other diagnostic groups, and from one another. METHODS: In this retrospective multicenter cohort study, we included 1213 patients (age 67 ± 9, 48% females) from two memory clinic cohorts: 116 frontotemporal dementia, 341 Alzheimer's disease, 66 Dementia with Lewy bodies, 40 vascular dementia, 104 other dementias, 229 mild cognitive impairment, and 317 subjective cognitive decline. Three MRI atrophy biomarkers were derived from the normalized volumes of automatically segmented cortical regions: 1) the anterior vs. posterior index, 2) the asymmetry index, and 3) the temporal pole left index. We used the following performance metrics: area under the receiver operating characteristic curve (AUC), sensitivity, and specificity. To account for the low prevalence of frontotemporal dementia we pursued a high specificity of 95%. Cross-validation was used in assessing the performance. The generalizability was assessed in an independent cohort (n = 200). RESULTS: The anterior vs. posterior index performed with an AUC of 83% for differentiation of frontotemporal dementia from all other diagnostic groups (Sensitivity = 59%, Specificity = 95%, positive likelihood ratio = 11.8, negative likelihood ratio = 0.4). The asymmetry index showed highest performance for separation of primary progressive aphasia and behavioral variant frontotemporal dementia (AUC = 85%, Sensitivity = 79%, Specificity = 92%, positive likelihood ratio = 9.9, negative likelihood ratio = 0.2), whereas the temporal pole left index was specific for detection of semantic variant primary progressive aphasia (AUC = 85%, Sensitivity = 82%, Specificity = 80%, positive likelihood ratio = 4.1, negative likelihood ratio = 0.2). The validation cohort provided corresponding results for the anterior vs. posterior index and temporal pole left index. CONCLUSION: This study presents three quantitative MRI biomarkers, which could provide additional information to the diagnostic assessment and assist clinicians in diagnosing frontotemporal dementia.

YNICL Journal 2019 Journal Article

Medial temporal lobe atrophy and posterior atrophy scales normative values

  • Matteo Cotta Ramusino
  • Daniele Altomare
  • Ruggero Bacchin
  • Silvia Ingala
  • Claudio Bnà
  • Matteo Bonetti
  • Alfredo Costa
  • Frederik Barkhof

OBJECTIVES: The medial temporal lobe atrophy (MTA) and the posterior atrophy (PA) scales allow to assess the degree hippocampal and parietal atrophy from magnetic resonance imaging (MRI) scans. Despite reliable, easy and widespread employment, appropriate normative values are still missing. We aim to provide norms for the Italian population. METHODS: Two independent raters assigned the highest MTA and PA score between hemispheres, based on 3D T1-weighted MRI of 936 Italian Brain Normative Archive subjects (age: mean ± SD: 50.2 ± 14.7, range: 20-84; MMSE>26 or CDR = 0). The inter-rater agreement was assessed with the absolute intraclass correlation coefficient (aICC). We assessed the association between MTA and PA scores and sociodemographic features and APOE status, and normative data were established by age decade based on percentile distributions. RESULTS: Raters agreed in 90% of cases for MTA (aICC = 0.86; 95% CI = 0.69-0.98) and in 86% for PA (aICC = 0.82; 95% CI = 0.58-0.98). For both rating scales, score distribution was skewed, with MTA = 0 in 38% of the population and PA = 0 in 52%, while a score ≥ 2 was only observed in 12% for MTA and in 10% for PA. Median denoted overall hippocampal (MTA: median = 1, IQR = 0-1) and parietal (PA: median = 0, IQR = 0-1) integrity. The 90th percentile of the age-specific distributions increased from 1 (at age 20-59) for both scales, to 2 for PA over age 60, and up to 4 for MTA over age 80. Gender, education and APOE status did not significantly affect the percentile distributions in the whole sample, nor in the subset over age 60. CONCLUSIONS: Our normative data for the MTA and PA scales are consistent with previous studies and overcome their main limitations (in particular uneven representation of ages and missing percentile distributions), defining the age-specific norms to be considered for proper brain atrophy assessment.

YNICL Journal 2019 Journal Article

Modeling grey matter atrophy as a function of time, aging or cognitive decline show different anatomical patterns in Alzheimer's disease

  • Ellen Dicks
  • Lisa Vermunt
  • Wiesje M. van der Flier
  • Pieter Jelle Visser
  • Frederik Barkhof
  • Philip Scheltens
  • Betty M. Tijms

BACKGROUND: Grey matter (GM) atrophy in Alzheimer's disease (AD) is most commonly modeled as a function of time. However, this approach does not take into account inter-individual differences in initial disease severity or changes due to aging. Here, we modeled GM atrophy within individuals across the AD clinical spectrum as a function of time, aging and MMSE, as a proxy for disease severity, and investigated how these models influence estimates of GM atrophy. METHODS: We selected 523 individuals from ADNI (100 preclinical AD, 288 prodromal AD, 135 AD dementia) with abnormal baseline amyloid PET/CSF and ≥1 year of MRI follow-up. We calculated total and 90 regional GM volumes for 2281 MRI scans (median [IQR]; 4 [3-5] scans per individual over 2 [1.6-4] years) and used linear mixed models to investigate atrophy as a function of time, aging and decline on MMSE. Analyses included clinical stage as interaction with the predictor and were corrected for baseline age, sex, education, field strength and total intracranial volume. We repeated analyses for a sample of participants with normal amyloid (n = 387) to assess whether associations were specific for amyloid pathology. RESULTS: < 0.05). Associations with time, aging and MMSE remained when accounting for these effects in the other models, suggesting that all measures explain part of the variance in GM atrophy. Repeating analyses in amyloid normal individuals, effects for time and aging showed similar widespread anatomical patterns, while associations with MMSE were largely reduced. CONCLUSION: Effects of time, aging and MMSE all explained variance in GM atrophy slopes within individuals. Associations with MMSE were weaker than those for time or age, but specific for amyloid pathology. This suggests that at least some of the atrophy observed in time or age models may not be specific to AD.

YNICL Journal 2019 Journal Article

Multi-study validation of data-driven disease progression models to characterize evolution of biomarkers in Alzheimer's disease

  • Damiano Archetti
  • Silvia Ingala
  • Vikram Venkatraghavan
  • Viktor Wottschel
  • Alexandra L. Young
  • Maura Bellio
  • Esther E. Bron
  • Stefan Klein

Understanding the sequence of biological and clinical events along the course of Alzheimer's disease provides insights into dementia pathophysiology and can help participant selection in clinical trials. Our objective is to train two data-driven computational models for sequencing these events, the Event Based Model (EBM) and discriminative-EBM (DEBM), on the basis of well-characterized research data, then validate the trained models on subjects from clinical cohorts characterized by less-structured data-acquisition protocols. Seven independent data cohorts were considered totalling 2389 cognitively normal (CN), 1424 mild cognitive impairment (MCI) and 743 Alzheimer's disease (AD) patients. The Alzheimer's Disease Neuroimaging Initiative (ADNI) data set was used as training set for the constriction of disease models while a collection of multi-centric data cohorts was used as test set for validation. Cross-sectional information related to clinical, cognitive, imaging and cerebrospinal fluid (CSF) biomarkers was used. Event sequences obtained with EBM and DEBM showed differences in the ordering of single biomarkers but according to both the first biomarkers to become abnormal were those related to CSF, followed by cognitive scores, while structural imaging showed significant volumetric decreases at later stages of the disease progression. Staging of test set subjects based on sequences obtained with both models showed good linear correlation with the Mini Mental State Examination score (R2 EBM = 0. 866; R2 DEBM = 0. 906). In discriminant analyses, significant differences (p-value ≤ 0. 05) between the staging of subjects from training and test sets were observed in both models. No significant difference between the staging of subjects from the training and test was observed (p-value > 0. 05) when considering a subset composed by 562 subjects for which all biomarker families (cognitive, imaging and CSF) are available. Event sequence obtained with DEBM recapitulates the heuristic models in a data-driven fashion and is clinically plausible. We demonstrated inter-cohort transferability of two disease progression models and their robustness in detecting AD phases. This is an important step towards the adoption of data-driven statistical models into clinical domain.

YNICL Journal 2019 Journal Article

Normal Aging Brain Collection Amsterdam (NABCA): A comprehensive collection of postmortem high-field imaging, neuropathological and morphometric datasets of non-neurological controls

  • Laura E. Jonkman
  • Yvon Galis-de Graaf
  • Marjolein Bulk
  • Eliane Kaaij
  • Petra J.W. Pouwels
  • Frederik Barkhof
  • Annemieke J.M. Rozemuller
  • Louise van der Weerd

Well-characterized, high-quality brain tissue of non-neurological control subjects is a prerequisite to study the healthy aging brain, and can serve as a control for the study of neurological disorders. The Normal Aging Brain Collection Amsterdam (NABCA) provides a comprehensive collection of post-mortem (ultra-)high-field MRI (3Tesla and 7 Tesla) and neuropathological datasets of non-neurological controls. By providing MRI within the pipeline, NABCA uniquely stimulates translational neurosciences; from molecular and morphometric tissue studies to the clinical setting. We describe our pipeline, including a description of our on-call autopsy team, donor selection, in situ and ex vivo post-mortem MRI protocols, brain dissection and neuropathological diagnosis. A demographic, radiological and pathological overview of five selected cases on all these aspects is provided. Additionally, information is given on data management, data and tissue application procedures, including review by a scientific advisory board, and setting up a material transfer agreement before distribution of tissue. Finally, we focus on future prospects, which includes laying the foundation for a unique platform for neuroanatomical, histopathological and neuro-radiological education, of professionals, students and the general (lay) audience.

YNICL Journal 2019 Journal Article

Novel imaging phantom for accurate and robust measurement of brain atrophy rates using clinical MRI

  • Houshang Amiri
  • Iman Brouwer
  • Joost P.A. Kuijer
  • Jan C. de Munck
  • Frederik Barkhof
  • Hugo Vrenken

Brain volume loss, or atrophy, has been proven to be an important characteristic of neurological diseases such as Alzheimer's disease and multiple sclerosis. To use atrophy rate as a reliable clinical biomarker and to increase statistical power in clinical treatment trials, measurement variability needs to be minimized. Among other sources, systematic differences between different MR scanners are suspected to contribute to this variability. In this study we developed and performed initial validation tests of an MR-compatible phantom and analysis software for robust and reliable evaluation of the brain volume loss. The phantom contained three inflatable models of brain structures, i.e. cerebral hemisphere, putamen, and caudate nucleus. Software to reliably quantify volumes form the phantom images was also developed. To validate the method, the phantom was imaged using 3D T1-weighted protocols at three clinical 3T MR scanners from different vendors. Calculated volume change from MRI was compared with the known applied volume change using ICC and mean absolute difference. As assessed by the ICC, the agreement between our developed software and the applied volume change for different structures ranged from 0.999-1 for hemisphere, 0.976-0.998 for putamen, and 0.985-0.999 for caudate nucleus. The mean absolute differences between measured and applied volume change were 109-332 μL for hemisphere, 2.9-11.9 μL for putamen, and 2.2-10.1 μL for caudate nucleus. This method offers a reliable and robust measurement of volume change using MR images and could potentially be used to standardize clinical measurement of atrophy rates.

YNICL Journal 2019 Journal Article

Prion disease diagnosis using subject-specific imaging biomarkers within a multi-kernel Gaussian process

  • Liane S. Canas
  • Carole H. Sudre
  • Enrico De Vita
  • Akin Nihat
  • Tze How Mok
  • Catherine F. Slattery
  • Ross W. Paterson
  • Alexander J.M. Foulkes

Prion diseases are a group of rare neurodegenerative conditions characterised by a high rate of progression and highly heterogeneous phenotypes. Whilst the most common form of prion disease occurs sporadically (sporadic Creutzfeldt-Jakob disease, sCJD), other forms are caused by prion protein gene mutations, or exposure to prions in the diet or by medical procedures, such us surgeries. To date, there are no accurate quantitative imaging biomarkers that can be used to predict the future clinical diagnosis of a healthy subject, or to quantify the progression of symptoms over time. Besides, CJD is commonly mistaken for other forms of dementia. Due to the heterogeneity of phenotypes and the lack of a consistent geometrical pattern of disease progression, the approaches used to study other types of neurodegenerative diseases are not satisfactory to capture the progression of human form of prion disease. In this paper, using a tailored framework, we aim to classify and stratify patients with prion disease, according to the severity of their illness. The framework is initialised with the extraction of subject-specific imaging biomarkers. The extracted biomakers are then combined with genetic and demographic information within a Gaussian Process classifier, used to calculate the probability of a subject to be diagnosed with prion disease in the next year. We evaluate the effectiveness of the proposed method in a cohort of patients with inherited and sporadic forms of prion disease. The model has shown to be effective in the prediction of both inherited CJD (92% of accuracy) and sporadic CJD (95% of accuracy). However the model has shown to be less effective when used to stratify the different stages of the disease, in which the average accuracy is 85%, whilst the recall is 59%. Finally, our framework was extended as a differential diagnosis tool to identify both forms of CJD among another neurodegenerative disease. In summary we have developed a novel method for prion disease diagnosis and prediction of clinical onset using multiple sources of features, which may have use in other disorders with heterogeneous imaging features.

YNICL Journal 2019 Journal Article

SVM recursive feature elimination analyses of structural brain MRI predicts near-term relapses in patients with clinically isolated syndromes suggestive of multiple sclerosis

  • Viktor Wottschel
  • Declan T. Chard
  • Christian Enzinger
  • Massimo Filippi
  • Jette L. Frederiksen
  • Claudio Gasperini
  • Antonio Giorgio
  • Maria A. Rocca

Machine learning classification is an attractive approach to automatically differentiate patients from healthy subjects, and to predict future disease outcomes. A clinically isolated syndrome (CIS) is often the first presentation of multiple sclerosis (MS), but it is difficult at onset to predict who will have a second relapse and hence convert to clinically definite MS. In this study, we thus aimed to distinguish CIS converters from non-converters at onset of a CIS, using recursive feature elimination and weight averaging with support vector machines. We also sought to assess the influence of cohort size and cross-validation methods on the accuracy estimate of the classification. We retrospectively collected 400 patients with CIS from six European MAGNIMS MS centres. Patients underwent brain MRI at onset of a CIS according to local standard-of-care protocols. The diagnosis of clinically definite MS at one-year follow-up was the standard against which the accuracy of the model was tested. For each patient, we derived MRI-based features, such as grey matter probability, white matter lesion load, cortical thickness, and volume of specific cortical and white matter regions. Features with little contribution to the classification model were removed iteratively through an interleaved sample bootstrapping and feature averaging approach. Classification of CIS outcome at one-year follow-up was performed with 2-fold, 5-fold, 10-fold and leave-one-out cross-validation for each centre cohort independently and in all patients together. The estimated classification accuracy across centres ranged from 64.9% to 88.1% using 2-fold cross-validation and from 73% to 92.9% using leave-one-out cross-validation. The classification accuracy estimate was higher in single-centre, smaller data sets than in combinations of data sets, being the lowest when all patients were merged together. Regional MRI features such as WM lesions, grey matter probability in the thalamus and the precuneus or cortical thickness in the cuneus and inferior temporal gyrus predicted the occurrence of a second relapse in patients at onset of a CIS using support vector machines. The increased accuracy estimate of the classification achieved with smaller and single-centre samples may indicate a model bias (overfitting) when data points were limited, but also more homogeneous. We provide an overview of classifier performance from a range of cross-validation schemes to give insight into the variability across schemes. The proposed recursive feature elimination approach with weight averaging can be used both in single- and multi-centre data sets in order to bridge the gap between group-level comparisons and making predictions for individual patients.

YNICL Journal 2018 Journal Article

Hypometabolism of the posterior cingulate cortex is not restricted to Alzheimer's disease

  • Nienke M.E. Scheltens
  • Kars van der Weijden
  • Sofie M. Adriaanse
  • Danielle van Assema
  • Priscilla P. Oomen
  • Welmoed A. Krudop
  • Adriaan A. Lammertsma
  • Frederik Barkhof

F]FDG standard uptake value ratio (SUVr) cut-off, which was derived from the receiver operating characteristic (ROC) separating probable AD from CN, resulting in 9-14/33 bvFTD patients having PCC hypometabolism, depending on the reference tissue used. In conclusion, PCC hypometabolism is not restricted to AD.

YNICL Journal 2018 Journal Article

Progressive brain rich-club network disruption from clinically isolated syndrome towards multiple sclerosis

  • Ni Shu
  • Yunyun Duan
  • Jing Huang
  • Zhuoqiong Ren
  • Zheng Liu
  • Huiqing Dong
  • Frederik Barkhof
  • Kuncheng Li

Objective: To investigate the rich-club organization in clinically isolated syndrome (CIS) and multiple sclerosis (MS), and to characterize its relationships with physical disabilities and cognitive impairments. Methods: We constructed high-resolution white matter (WM) structural networks in 41 CIS, 32 MS and 35 healthy controls (HCs) using diffusion MRI and deterministic tractography. Group differences in rich-club organization, global and local network metrics were investigated. The relationship between the altered network metrics, brain lesions and clinical variables including EDSS, MMSE, PASAT, disease duration were calculated. Additionally, reproducibility analysis was performed using different parcellation schemes. Results: = 0.0004). The reduced strength of rich-club and feeder connections was correlated with cognitive impairments in the MS group. These results were independent of lesion distribution and reproducible across different brain parcellation schemes. Conclusion: The rich-club organization was disrupted in MS patients and relatively preserved in CIS. The disrupted rich-club connectivity was correlated with cognitive impairment in MS. These findings suggest that impaired rich-club connectivity is an essential feature of progressive structural network disruption, heralding the development of clinical disability in MS.

YNICL Journal 2018 Journal Article

Urgent challenges in quantification and interpretation of brain grey matter atrophy in individual MS patients using MRI

  • Houshang Amiri
  • Alexandra de Sitter
  • Kerstin Bendfeldt
  • Marco Battaglini
  • Claudia A.M. Gandini Wheeler-Kingshott
  • Massimiliano Calabrese
  • Jeroen J.G. Geurts
  • Maria A. Rocca

Atrophy of the brain grey matter (GM) is an accepted and important feature of multiple sclerosis (MS). However, its accurate measurement is hampered by various technical, pathological and physiological factors. As a consequence, it is challenging to investigate the role of GM atrophy in the disease process as well as the effect of treatments that aim to reduce neurodegeneration. In this paper we discuss the most important challenges currently hampering the measurement and interpretation of GM atrophy in MS. The focus is on measurements that are obtained in individual patients rather than on group analysis methods, because of their importance in clinical trials and ultimately in clinical care. We discuss the sources and possible solutions of the current challenges, and provide recommendations to achieve reliable measurement and interpretation of brain GM atrophy in MS.

YNICL Journal 2017 Journal Article

Agreement of MSmetrix with established methods for measuring cross-sectional and longitudinal brain atrophy

  • Martijn D. Steenwijk
  • Houshang Amiri
  • Menno M. Schoonheim
  • Alexandra de Sitter
  • Frederik Barkhof
  • Petra J.W. Pouwels
  • Hugo Vrenken

INTRODUCTION: MRI data. METHODS: data from 100 MS patients and 20 matched healthy controls. In fifty of the MS patients a second time point was available. In this subgroup, we additionally analyzed the whole-brain and GM volume change using the longitudinal pipeline of MSmetrix and compared the results with those of FreeSurfer (longitudinal pipeline) and SIENA. RESULTS: . Among the cross-sectional methods, Deming regression analyses revealed proportional errors particularly in MSmetrix and SPM. The mean difference percentage brain volume change (PBVC) was lowest between longitudinal MSmetrix and SIENA (+ 0.16 ± 0.91%). A strong proportional error was present between longitudinal percentage gray matter volume change (PGVC) measures of MSmetrix and FreeSurfer (slope = 2.48). All longitudinal methods were sensitive to the MRI hardware upgrade that occurred during the time of the study. CONCLUSION: MSmetrix, FreeSurfer, FSL and SPM show differences in atrophy measurements, even at the whole-brain level, that are large compared to typical atrophy rates observed in MS. Especially striking are the proportional errors between methods. Cross-sectional MSmetrix behaved similarly to SPM, both in terms of mean volume difference as well as proportional error. Longitudinal MSmetrix behaved most similar to SIENA. Our results indicate that brain volume measurement and normalization from T1-weighted images remains an unsolved problem that requires much more attention.

YNICL Journal 2017 Journal Article

Connectomic profile and clinical phenotype in newly diagnosed glioma patients

  • Jolanda Derks
  • Anne R. Dirkson
  • Philip C. de Witt Hamer
  • Quinten van Geest
  • Hanneke E. Hulst
  • Frederik Barkhof
  • Petra J.W. Pouwels
  • Jeroen J.G. Geurts

Gliomas are primary brain tumors, originating from the glial cells in the brain. In contrast to the more traditional view of glioma as a localized disease, it is becoming clear that global brain functioning is impacted, even with respect to functional communication between brain regions remote from the tumor itself. However, a thorough investigation of glioma-related functional connectomic profiles is lacking. Therefore, we constructed functional brain networks using functional MR scans of 71 glioma patients and 19 matched healthy controls using the automated anatomical labelling (AAL) atlas and interregional Pearson correlation coefficients. The frequency distributions across connectivity values were calculated to depict overall connectomic profiles and quantitative features of these distributions (full-width half maximum (FWHM), peak position, peak height) were calculated. Next, we investigated the spatial distribution of the connectomic profile. We defined hub locations based on the literature and determined connectivity (1) between hubs, (2) between hubs and non-hubs, and (3) between non-hubs. Results show that patients had broader and flatter connectivity distributions compared to controls. Spatially, glioma patients particularly showed increased connectivity between non-hubs and hubs. Furthermore, connectivity distributions and hub-non-hub connectivity differed within the patient group according to tumor grade, while relating to Karnofsky performance status and progression-free survival. In conclusion, newly diagnosed glioma patients have globally altered functional connectomic profiles, which mainly affect hub connectivity and relate to clinical phenotypes. These findings underscore the promise of using connectomics as a future biomarker in this patient population.

YNICL Journal 2017 Journal Article

Five-class differential diagnostics of neurodegenerative diseases using random undersampling boosting

  • Tong Tong
  • Christian Ledig
  • Ricardo Guerrero
  • Andreas Schuh
  • Juha Koikkalainen
  • Antti Tolonen
  • Hanneke Rhodius
  • Frederik Barkhof

Differentiating between different types of neurodegenerative diseases is not only crucial in clinical practice when treatment decisions have to be made, but also has a significant potential for the enrichment of clinical trials. The purpose of this study is to develop a classification framework for distinguishing the four most common neurodegenerative diseases, including Alzheimer's disease, frontotemporal lobe degeneration, Dementia with Lewy bodies and vascular dementia, as well as patients with subjective memory complaints. Different biomarkers including features from images (volume features, region-wise grading features) and non-imaging features (CSF measures) were extracted for each subject. In clinical practice, the prevalence of different dementia types is imbalanced, posing challenges for learning an effective classification model. Therefore, we propose the use of the RUSBoost algorithm in order to train classifiers and to handle the class imbalance training problem. Furthermore, a multi-class feature selection method based on sparsity is integrated into the proposed framework to improve the classification performance. It also provides a way for investigating the importance of different features and regions. Using a dataset of 500 subjects, the proposed framework achieved a high accuracy of 75.2% with a balanced accuracy of 69.3% for the five-class classification using ten-fold cross validation, which is significantly better than the results using support vector machine or random forest, demonstrating the feasibility of the proposed framework to support clinical decision making.

YNIMG Journal 2017 Journal Article

Performance of five research-domain automated WM lesion segmentation methods in a multi-center MS study

  • Alexandra de Sitter
  • Martijn D. Steenwijk
  • Aurélie Ruet
  • Adriaan Versteeg
  • Yaou Liu
  • Ronald A. van Schijndel
  • Petra J.W. Pouwels
  • Iris D. Kilsdonk

Background and Purpose In vivoidentification of white matter lesions plays a key-role in evaluation of patients with multiple sclerosis (MS). Automated lesion segmentation methods have been developed to substitute manual outlining, but evidence of their performance in multi-center investigations is lacking. In this work, five research-domain automated segmentation methods were evaluated using a multi-center MS dataset. Methods 70 MS patients (median EDSS of 2. 0 [range 0. 0–6. 5]) were included from a six-center dataset of the MAGNIMS Study Group (www. magnims. eu) which included 2D FLAIR and 3D T1 images with manual lesion segmentation as a reference. Automated lesion segmentations were produced using five algorithms: Cascade; Lesion Segmentation Toolbox (LST) with both the Lesion growth algorithm (LGA) and the Lesion prediction algorithm (LPA); Lesion-Topology preserving Anatomical Segmentation (Lesion-TOADS); and k-Nearest Neighbor with Tissue Type Priors (kNN-TTP). Main software parameters were optimized using a training set (N = 18), and formal testing was performed on the remaining patients (N = 52). To evaluate volumetric agreement with the reference segmentations, intraclass correlation coefficient (ICC) as well as mean difference in lesion volumes between the automated and reference segmentations were calculated. The Similarity Index (SI), False Positive (FP) volumes and False Negative (FN) volumes were used to examine spatial agreement. All analyses were repeated using a leave-one-center-out design to exclude the center of interest from the training phase to evaluate the performance of the method on ‘unseen’ center. Results Compared to the reference mean lesion volume (4. 85 ± 7. 29 mL), the methods displayed a mean difference of 1. 60 ± 4. 83 (Cascade), 2. 31 ± 7. 66 (LGA), 0. 44 ± 4. 68 (LPA), 1. 76 ± 4. 17 (Lesion-TOADS) and −1. 39 ± 4. 10 mL (kNN-TTP). The ICCs were 0. 755, 0. 713, 0. 851, 0. 806 and 0. 723, respectively. Spatial agreement with reference segmentations was higher for LPA (SI = 0. 37 ± 0. 23), Lesion-TOADS (SI = 0. 35 ± 0. 18) and kNN-TTP (SI = 0. 44 ± 0. 14) than for Cascade (SI = 0. 26 ± 0. 17) or LGA (SI = 0. 31 ± 0. 23). All methods showed highly similar results when used on data from a center not used in software parameter optimization. Conclusion The performance of the methods in this multi-center MS dataset was moderate, but appeared to be robust even with new datasets from centers not included in training the automated methods.

YNIMG Journal 2017 Journal Article

The European DTI Study on Dementia — A multicenter DTI and MRI study on Alzheimer's disease and Mild Cognitive Impairment

  • Katharina Brueggen
  • Michel J. Grothe
  • Martin Dyrba
  • Andreas Fellgiebel
  • Florian Fischer
  • Massimo Filippi
  • Federica Agosta
  • Peter Nestor

The European DTI Study on Dementia (EDSD) is a multicenter framework created to study the diagnostic accuracy and inter-site variability of DTI-derived markers in patients with manifest and prodromal Alzheimer's disease (AD). The dynamically growing database presently includes 493 DTI, 512 T1-weighted MRI, and 300 FLAIR scans from patients with AD dementia, patients with Mild Cognitive Impairment (MCI) and matched Healthy Controls, acquired on 13 different scanner platforms. The imaging data is publicly available, along with the subjects' demographic and clinical characterization. Detailed neuropsychological information, cerebrospinal fluid information on biomarkers and clinical follow-up diagnoses are included for a subset of subjects. This paper describes the rationale and structure of the EDSD, summarizes the available data, and explains how to gain access to the database. The EDSD is a useful database for researchers seeking to investigate the contribution of DTI to dementia diagnostics.

YNICL Journal 2016 Journal Article

Differential diagnosis of neurodegenerative diseases using structural MRI data

  • Juha Koikkalainen
  • Hanneke Rhodius-Meester
  • Antti Tolonen
  • Frederik Barkhof
  • Betty Tijms
  • Afina W. Lemstra
  • Tong Tong
  • Ricardo Guerrero

Different neurodegenerative diseases can cause memory disorders and other cognitive impairments. The early detection and the stratification of patients according to the underlying disease are essential for an efficient approach to this healthcare challenge. This emphasizes the importance of differential diagnostics. Most studies compare patients and controls, or Alzheimer's disease with one other type of dementia. Such a bilateral comparison does not resemble clinical practice, where a clinician is faced with a number of different possible types of dementia. Here we studied which features in structural magnetic resonance imaging (MRI) scans could best distinguish four types of dementia, Alzheimer's disease, frontotemporal dementia, vascular dementia, and dementia with Lewy bodies, and control subjects. We extracted an extensive set of features quantifying volumetric and morphometric characteristics from T1 images, and vascular characteristics from FLAIR images. Classification was performed using a multi-class classifier based on Disease State Index methodology. The classifier provided continuous probability indices for each disease to support clinical decision making. A dataset of 504 individuals was used for evaluation. The cross-validated classification accuracy was 70.6% and balanced accuracy was 69.1% for the five disease groups using only automatically determined MRI features. Vascular dementia patients could be detected with high sensitivity (96%) using features from FLAIR images. Controls (sensitivity 82%) and Alzheimer's disease patients (sensitivity 74%) could be accurately classified using T1-based features, whereas the most difficult group was the dementia with Lewy bodies (sensitivity 32%). These results were notable better than the classification accuracies obtained with visual MRI ratings (accuracy 44.6%, balanced accuracy 51.6%). Different quantification methods provided complementary information, and consequently, the best results were obtained by utilizing several quantification methods. The results prove that automatic quantification methods and computerized decision support methods are feasible for clinical practice and provide comprehensive information that may help clinicians in the diagnosis making.

YNIMG Journal 2016 Journal Article

Longitudinal reproducibility of default-mode network connectivity in healthy elderly participants: A multicentric resting-state fMRI study

  • Jorge Jovicich
  • Ludovico Minati
  • Moira Marizzoni
  • Rocco Marchitelli
  • Roser Sala-Llonch
  • David Bartrés-Faz
  • Jennifer Arnold
  • Jens Benninghoff

To date, limited data are available regarding the inter-site consistency of test–retest reproducibility of functional connectivity measurements, in particular with regard to integrity of the Default Mode Network (DMN) in elderly participants. We implemented a harmonized resting-state fMRI protocol on 13 clinical scanners at 3. 0T using vendor-provided sequences. Each site scanned a group of 5 healthy elderly participants twice, at least a week apart. We evaluated inter-site differences and test–retest reproducibility of both temporal signal-to-noise ratio (tSNR) and functional connectivity measurements derived from: i) seed-based analysis (SBA) with seed in the posterior cingulate cortex (PCC), ii) group independent component analysis (ICA) separately for each site (site ICA), and iii) consortium ICA, with group ICA across the whole consortium. Despite protocol harmonization, significant and quantitatively important inter-site differences remained in the tSNR of resting-state fMRI data; these were plausibly driven by hardware and pulse sequence differences across scanners which could not be harmonized. Nevertheless, the tSNR test–retest reproducibility in the consortium was high (ICC=0. 81). The DMN was consistently extracted across all sites and analysis methods. While significant inter-site differences in connectivity scores were found, there were no differences in the associated test–retest error. Overall, ICA measurements were more reliable than PCC-SBA, with site ICA showing higher reproducibility than consortium ICA. Across the DMN nodes, the PCC yielded the most reliable measurements (≈4% test–retest error, ICC=0. 85), the medial frontal cortex the least reliable (≈12%, ICC=0. 82) and the lateral parietal cortices were in between (site ICA). Altogether these findings support usage of harmonized multisite studies of resting-state functional connectivity to characterize longitudinal effects in studies that assess disease progression and treatment response.

YNICL Journal 2015 Journal Article

Automatic segmentation and volumetry of multiple sclerosis brain lesions from MR images

  • Saurabh Jain
  • Diana M. Sima
  • Annemie Ribbens
  • Melissa Cambron
  • Anke Maertens
  • Wim Van Hecke
  • Johan De Mey
  • Frederik Barkhof

The location and extent of white matter lesions on magnetic resonance imaging (MRI) are important criteria for diagnosis, follow-up and prognosis of multiple sclerosis (MS). Clinical trials have shown that quantitative values, such as lesion volumes, are meaningful in MS prognosis. Manual lesion delineation for the segmentation of lesions is, however, time-consuming and suffers from observer variability. In this paper, we propose MSmetrix, an accurate and reliable automatic method for lesion segmentation based on MRI, independent of scanner or acquisition protocol and without requiring any training data. In MSmetrix, 3D T1-weighted and FLAIR MR images are used in a probabilistic model to detect white matter (WM) lesions as an outlier to normal brain while segmenting the brain tissue into grey matter, WM and cerebrospinal fluid. The actual lesion segmentation is performed based on prior knowledge about the location (within WM) and the appearance (hyperintense on FLAIR) of lesions. The accuracy of MSmetrix is evaluated by comparing its output with expert reference segmentations of 20 MRI datasets of MS patients. Spatial overlap (Dice) between the MSmetrix and the expert lesion segmentation is 0.67 ± 0.11. The intraclass correlation coefficient (ICC) equals 0.8 indicating a good volumetric agreement between the MSmetrix and expert labelling. The reproducibility of MSmetrix' lesion volumes is evaluated based on 10 MS patients, scanned twice with a short interval on three different scanners. The agreement between the first and the second scan on each scanner is evaluated through the spatial overlap and absolute lesion volume difference between them. The spatial overlap was 0.69 ± 0.14 and absolute total lesion volume difference between the two scans was 0.54 ± 0.58 ml. Finally, the accuracy and reproducibility of MSmetrix compare favourably with other publicly available MS lesion segmentation algorithms, applied on the same data using default parameter settings.

YNICL Journal 2015 Journal Article

Joint assessment of white matter integrity, cortical and subcortical atrophy to distinguish AD from behavioral variant FTD: A two-center study

  • Christiane Möller
  • Anne Hafkemeijer
  • Yolande A.L. Pijnenburg
  • Serge A.R.B. Rombouts
  • Jeroen van der Grond
  • Elise Dopper
  • John van Swieten
  • Adriaan Versteeg

We investigated the ability of cortical and subcortical gray matter (GM) atrophy in combination with white matter (WM) integrity to distinguish behavioral variant frontotemporal dementia (bvFTD) from Alzheimer's disease (AD) and from controls using voxel-based morphometry, subcortical structure segmentation, and tract-based spatial statistics. To determine which combination of MR markers differentiated the three groups with the highest accuracy, we conducted discriminant function analyses. Adjusted for age, sex and center, both types of dementia had more GM atrophy, lower fractional anisotropy (FA) and higher mean (MD), axial (L1) and radial diffusivity (L23) values than controls. BvFTD patients had more GM atrophy in orbitofrontal and inferior frontal areas than AD patients. In addition, caudate nucleus and nucleus accumbens were smaller in bvFTD than in AD. FA values were lower; MD, L1 and L23 values were higher, especially in frontal areas of the brain for bvFTD compared to AD patients. The combination of cortical GM, hippocampal volume and WM integrity measurements, classified 97-100% of controls, 81-100% of AD and 67-75% of bvFTD patients correctly. Our results suggest that WM integrity measures add complementary information to measures of GM atrophy, thereby improving the classification between AD and bvFTD.

YNIMG Journal 2015 Journal Article

Standardized evaluation of algorithms for computer-aided diagnosis of dementia based on structural MRI: The CADDementia challenge

  • Esther E. Bron
  • Marion Smits
  • Wiesje M. van der Flier
  • Hugo Vrenken
  • Frederik Barkhof
  • Philip Scheltens
  • Janne M. Papma
  • Rebecca M.E. Steketee

Algorithms for computer-aided diagnosis of dementia based on structural MRI have demonstrated high performance in the literature, but are difficult to compare as different data sets and methodology were used for evaluation. In addition, it is unclear how the algorithms would perform on previously unseen data, and thus, how they would perform in clinical practice when there is no real opportunity to adapt the algorithm to the data at hand. To address these comparability, generalizability and clinical applicability issues, we organized a grand challenge that aimed to objectively compare algorithms based on a clinically representative multi-center data set. Using clinical practice as the starting point, the goal was to reproduce the clinical diagnosis. Therefore, we evaluated algorithms for multi-class classification of three diagnostic groups: patients with probable Alzheimer's disease, patients with mild cognitive impairment and healthy controls. The diagnosis based on clinical criteria was used as reference standard, as it was the best available reference despite its known limitations. For evaluation, a previously unseen test set was used consisting of 354 T1-weighted MRI scans with the diagnoses blinded. Fifteen research teams participated with a total of 29 algorithms. The algorithms were trained on a small training set (n=30) and optionally on data from other sources (e. g. , the Alzheimer's Disease Neuroimaging Initiative, the Australian Imaging Biomarkers and Lifestyle flagship study of aging). The best performing algorithm yielded an accuracy of 63. 0% and an area under the receiver-operating-characteristic curve (AUC) of 78. 8%. In general, the best performances were achieved using feature extraction based on voxel-based morphometry or a combination of features that included volume, cortical thickness, shape and intensity. The challenge is open for new submissions via the web-based framework: http: //caddementia. grand-challenge. org.

YNICL Journal 2014 Journal Article

Changes in MEG resting-state networks are related to cognitive decline in type 1 diabetes mellitus patients

  • Matteo Demuru
  • Eelco van Duinkerken
  • Matteo Fraschini
  • Francesco Marrosu
  • Frank J. Snoek
  • Frederik Barkhof
  • Martin Klein
  • Michaela Diamant

OBJECTIVE: Integrity of resting-state functional brain networks (RSNs) is important for proper cognitive functioning. In type 1 diabetes mellitus (T1DM) cognitive decrements are commonly observed, possibly due to alterations in RSNs, which may vary according to microvascular complication status. Thus, we tested the hypothesis that functional connectivity in RSNs differs according to clinical status and correlates with cognition in T1DM patients, using an unbiased approach with high spatio-temporal resolution functional network. METHODS: Resting-state magnetoencephalographic (MEG) data for T1DM patients with (n = 42) and without (n = 41) microvascular complications and 33 healthy participants were recorded. MEG time-series at source level were reconstructed using a recently developed atlas-based beamformer. Functional connectivity within classical frequency bands, estimated by the phase lag index (PLI), was calculated within eight commonly found RSNs. Neuropsychological tests were used to assess cognitive performance, and the relation with RSNs was evaluated. RESULTS: Significant differences in terms of RSN functional connectivity between the three groups were observed in the lower alpha band, in the default-mode (DMN), executive control (ECN) and sensorimotor (SMN) RSNs. T1DM patients with microvascular complications showed the weakest functional connectivity in these networks relative to the other groups. For DMN, functional connectivity was higher in patients without microangiopathy relative to controls (all p < 0.05). General cognitive performance for both patient groups was worse compared with healthy controls. Lower DMN alpha band functional connectivity correlated with poorer general cognitive ability in patients with microvascular complications. DISCUSSION: Altered RSN functional connectivity was found in T1DM patients depending on clinical status. Lower DMN functional connectivity was related to poorer cognitive functioning. These results indicate that functional connectivity may play a key role in T1DM-related cognitive dysfunction.

YNIMG Journal 2014 Journal Article

Hippocampal volume change measurement: Quantitative assessment of the reproducibility of expert manual outlining and the automated methods FreeSurfer and FIRST

  • Emma R. Mulder
  • Remko A. de Jong
  • Dirk L. Knol
  • Ronald A. van Schijndel
  • Keith S. Cover
  • Pieter J. Visser
  • Frederik Barkhof
  • Hugo Vrenken

Background To measure hippocampal volume change in Alzheimer's disease (AD) or mild cognitive impairment (MCI), expert manual delineation is often used because of its supposed accuracy. It has been suggested that expert outlining yields poorer reproducibility as compared to automated methods, but this has not been investigated. Aim To determine the reproducibilities of expert manual outlining and two common automated methods for measuring hippocampal atrophy rates in healthy aging, MCI and AD. Methods From the Alzheimer's Disease Neuroimaging Initiative (ADNI), 80 subjects were selected: 20 patients with AD, 40 patients with mild cognitive impairment (MCI) and 20 healthy controls (HCs). Left and right hippocampal volume change between baseline and month-12 visit was assessed by using expert manual delineation, and by the automated software packages FreeSurfer (longitudinal processing stream) and FIRST. To assess reproducibility of the measured hippocampal volume change, both back-to-back (BTB) MPRAGE scans available for each visit were analyzed. Hippocampal volume change was expressed in μL, and as a percentage of baseline volume. Reproducibility of the 1-year hippocampal volume change was estimated from the BTB measurements by using linear mixed model to calculate the limits of agreement (LoA) of each method, reflecting its measurement uncertainty. Using the delta method, approximate p-values were calculated for the pairwise comparisons between methods. Statistical analyses were performed both with inclusion and exclusion of visibly incorrect segmentations. Results Visibly incorrect automated segmentation in either one or both scans of a longitudinal scan pair occurred in 7. 5% of the hippocampi for FreeSurfer and in 6. 9% of the hippocampi for FIRST. After excluding these failed cases, reproducibility analysis for 1-year percentage volume change yielded LoA of ±7. 2% for FreeSurfer, ±9. 7% for expert manual delineation, and ±10. 0% for FIRST. Methods ranked the same for reproducibility of 1-yearμL volume change, with LoA of ±218μL for FreeSurfer, ±319μL for expert manual delineation, and ±333μL for FIRST. Approximate p-values indicated that reproducibility was better for FreeSurfer than for manual or FIRST, and that manual and FIRST did not differ. Inclusion of failed automated segmentations led to worsening of reproducibility of both automated methods for 1-year raw and percentage volume change. Conclusion Quantitative reproducibility values of 1-year microliter and percentage hippocampal volume change were roughly similar between expert manual outlining, FIRST and FreeSurfer, but FreeSurfer reproducibility was statistically significantly superior to both manual outlining and FIRST after exclusion of failed segmentations.

YNICL Journal 2013 Journal Article

Accurate white matter lesion segmentation by k nearest neighbor classification with tissue type priors (kNN-TTPs)

  • Martijn D. Steenwijk
  • Petra J.W. Pouwels
  • Marita Daams
  • Jan Willem van Dalen
  • Matthan W.A. Caan
  • Edo Richard
  • Frederik Barkhof
  • Hugo Vrenken

INTRODUCTION: The segmentation and volumetric quantification of white matter (WM) lesions play an important role in monitoring and studying neurological diseases such as multiple sclerosis (MS) or cerebrovascular disease. This is often interactively done using 2D magnetic resonance images. Recent developments in acquisition techniques allow for 3D imaging with much thinner sections, but the large number of images per subject makes manual lesion outlining infeasible. This warrants the need for a reliable automated approach. Here we aimed to improve k nearest neighbor (kNN) classification of WM lesions by optimizing intensity normalization and using spatial tissue type priors (TTPs). METHODS: The kNN-TTP method used kNN classification with 3.0 T 3DFLAIR and 3DT1 intensities as well as MNI-normalized spatial coordinates as features. Additionally, TTPs were computed by nonlinear registration of data from healthy controls. Intensity features were normalized using variance scaling, robust range normalization or histogram matching. The algorithm was then trained and evaluated using a leave-one-out experiment among 20 patients with MS against a reference segmentation that was created completely manually. The performance of each normalization method was evaluated both with and without TTPs in the feature set. Volumetric agreement was evaluated using intra-class coefficient (ICC), and voxelwise spatial agreement was evaluated using Dice similarity index (SI). Finally, the robustness of the method across different scanners and patient populations was evaluated using an independent sample of elderly subjects with hypertension. RESULTS: The intensity normalization method had a large influence on the segmentation performance, with average SI values ranging from 0.66 to 0.72 when no TTPs were used. Independent of the normalization method, the inclusion of TTPs as features increased performance particularly by reducing the lesion detection error. Best performance was achieved using variance scaled intensity features and including TTPs in the feature set: this yielded ICC = 0.93 and average SI = 0.75 ± 0.08. Validation of the method in an independent sample of elderly subjects with hypertension, yielded even higher ICC = 0.96 and SI = 0.84 ± 0.14. CONCLUSION: Adding TTPs increases the performance of kNN based MS lesion segmentation methods. Best performance was achieved using variance scaling for intensity normalization and including TTPs in the feature set, showing excellent agreement with the reference segmentations across a wide range of lesion severity, irrespective of the scanner used or the pathological substrate of the lesions.

YNIMG Journal 2013 Journal Article

Functional segmentation of the hippocampus in the healthy human brain and in Alzheimer's disease

  • Mojtaba Zarei
  • Christian F. Beckmann
  • Maja A.A. Binnewijzend
  • Menno M. Schoonheim
  • Mohammad Ali Oghabian
  • Ernesto J. Sanz-Arigita
  • Philip Scheltens
  • Paul M. Matthews

In this study we segment the hippocampus according to functional connectivity assessed from resting state functional magnetic resonance images in healthy subjects and in patients with Alzheimer's disease (AD). We recorded the resting FMRI signal from 16 patients and 22 controls. We used seed-based functional correlation analyses to calculate partial correlations of all voxels in the hippocampus relative to characteristic regional signal changes in the thalamus, the prefrontal cortex (PFC) and the posterior cingulate cortex (PCC), while controlling for ventricular CSF and white matter signals. Group comparisons were carried out controlling for age, gender, hippocampal volume and brain volume. The strength of functional connectivity in each region also was correlated with neuropsychological measures. We found that the hippocampus can be segmented into three distinct functional subregions (head, body, and tail), according to the relative connectivity with PFC, PCC and thalamus, respectively. The AD group showed stronger hippocampus–PFC and weaker hippocampus–PCC functional connectivity, the magnitudes of which correlated with MMSE in both cases. The results are consistent with an adaptive role of the PFC in the context of progression of dysfunction in PCC during earlier stages of AD. Extension of our approach could integrate regional volume measures for the hippocampus with their functional connectivity patterns in ways that should increase sensitivity for assessment of AD onset and progression.

YNIMG Journal 2010 Journal Article

Combining shape and connectivity analysis: An MRI study of thalamic degeneration in Alzheimer's disease

  • Mojtaba Zarei
  • Brian Patenaude
  • Jessica Damoiseaux
  • Ciro Morgese
  • Steve Smith
  • Paul M. Matthews
  • Frederik Barkhof
  • Serge Rombouts

Alzheimer's disease (AD) is associated with neuronal loss not only in the hippocampus and amygdala but also in the thalamus. Anterodorsal, centromedial, and pulvinar nuclei are the main sites of degeneration in AD. Here we combined shape analysis and diffusion tensor imaging (DTI) tractography to study degeneration in AD in the thalamus and its connections. Structural and diffusion tensor MRI scans were obtained from 16 AD patients and 22 demographically similar healthy volunteers. The thalamus, hippocampus, and amygdala were automatically segmented using our locally developed algorithm, and group comparisons were carried out for each surface vertex. We also employed probabilistic diffusion tractography to obtain connectivity measures between individual thalamic voxels and hippocampus/amygdala voxels and to segment the internal medullary lamina (IML). Shape analysis showed significant bilateral regional atrophy in the dorsal–medial part of the thalamus in AD patients compared to controls. Probabilistic tractography demonstrated that these regions are mainly connected with the hippocampus, temporal, and prefrontal cortex. Intrathalamic FA comparisons showed reductions in the anterodorsal region of thalamus. Intrathalamic tractography from this region revealed that the IML was significantly smaller in AD patients than in controls. We suggest that these changes can be attributed to the degeneration of the anterodorsal and intralaminar nuclei, respectively. In addition, based on previous neuropathological reports, ventral and dorsal–medial shape change in the thalamus in AD patients is likely to be driven by IML atrophy. This combined shape and connectivity analysis provides MRI evidence of regional thalamic degeneration in AD.

YNIMG Journal 2008 Journal Article

Segmentation of age-related white matter changes in a clinical multi-center study

  • Tim B. Dyrby
  • Egill Rostrup
  • William F.C. Baaré
  • Elisabeth C.W. van Straaten
  • Frederik Barkhof
  • Hugo Vrenken
  • Stefan Ropele
  • Reinhold Schmidt

Age-related white matter changes (WMC) are thought to be a marker of vascular pathology, and have been associated with motor and cognitive deficits. In the present study, an optimized artificial neural network was used as an automatic segmentation method to produce probabilistic maps of WMC in a clinical multi-center study. The neural network uses information from T1- and T2-weighted and fluid attenuation inversion recovery (FLAIR) magnetic resonance (MR) scans, neighboring voxels and spatial location. Generalizability of the neural network was optimized by including the Optimal Brain Damage (OBD) pruning method in the training stage. Six optimized neural networks were produced to investigate the impact of different input information on WMC segmentation. The automatic segmentation method was applied to MR scans of 362 non-demented elderly subjects from 11 centers in the European multi-center study Leukoaraiosis And Disability (LADIS). Semi-manually delineated WMC were used for validating the segmentation produced by the neural networks. The neural network segmentation demonstrated high consistency between subjects and centers, making it a promising technique for large studies. For WMC volumes less than 10 ml, an increasing discrepancy between semi-manual and neural network segmentation was observed using the similarity index (SI) measure. The use of all three image modalities significantly improved cross-center generalizability compared to neural networks using the FLAIR image only. Expert knowledge not available to the neural networks was a minor source of discrepancy, while variation in MR scan quality constituted the largest source of error.

YNIMG Journal 2007 Journal Article

Regional brain atrophy development is related to specific aspects of clinical dysfunction in multiple sclerosis

  • Bas Jasperse
  • Hugo Vrenken
  • Ernesto Sanz-Arigita
  • Vincent de Groot
  • Stephen M. Smith
  • Chris H. Polman
  • Frederik Barkhof

Brain atrophy in multiple sclerosis (MS) is thought to reflect irreversible tissue damage leading to persistent clinical deficit. Little is known about the rate of atrophy in specific brain regions in relation to specific clinical deficits. We determined the displacement of the brain surface between two T1-weighted MRI images obtained at baseline and after a median follow-up time of 2. 2 years for 79 recently diagnosed, mildly disabled MS patients. Voxel- and cluster-wise permutation-based statistics were used to identify brain regions in which atrophy development was significantly related to Expanded Disability Status Scale (EDSS), Timed Walk Test (TWT), Paced Auditory Serial Addition Test (PASAT) and 9-Hole Peg Test (HPT). Clusters were considered significant at a corrected cluster-wise p-value of 0. 05. Worse EDSS change-score and worse follow-up EDSS were related to atrophy development of periventricular and brainstem regions and right-sided parietal, occipital and temporal regions. Worse PASAT at follow-up was significantly related to atrophy of the ventricles. A worse TWT change-score and worse follow-up TWT were exclusively related to atrophy around the ventricles and of the brainstem. Worse HPT change-score and worse follow-up HPT of either arm were significantly related to the atrophy of widely distributed peripheral regions, as well as atrophy of periventricular and brainstem regions. Our findings suggest that decline in ambulatory function is related to atrophy of central brain regions exclusively, whereas decline in neurologically more complex tasks for coordinated hand function is related to atrophy of both central and peripheral brain regions.

YNIMG Journal 2006 Journal Article

Automated segmentation of multiple sclerosis lesion subtypes with multichannel MRI

  • Ying Wu
  • Simon K. Warfield
  • I. Leng Tan
  • William M. Wells
  • Dominik S. Meier
  • Ronald A. van Schijndel
  • Frederik Barkhof
  • Charles R.G. Guttmann

Purpose. To automatically segment multiple sclerosis (MS) lesions into three subtypes (i. e. , enhancing lesions, T1 “black holes”, T2 hyperintense lesions). Materials and methods. Proton density-, T2- and contrast-enhanced T1-weighted brain images of 12 MR scans were pre-processed through intracranial cavity (IC) extraction, inhomogeneity correction and intensity normalization. Intensity-based statistical k-nearest neighbor (k-NN) classification was combined with template-driven segmentation and partial volume artifact correction (TDS+) for segmentation of MS lesions subtypes and brain tissue compartments. Operator-supervised tissue sampling and parameter calibration were performed on 2 randomly selected scans and were applied automatically to the remaining 10 scans. Results from this three-channel TDS+ (3ch-TDS+) were compared to those from a previously validated two-channel TDS+ (2ch-TDS+) method. The results of both the 3ch-TDS+ and 2ch-TDS+ were also compared to manual segmentation performed by experts. Results. Intra-class correlation coefficients (ICC) of 3ch-TDS+ for all three subtypes of lesions were higher (ICC between 0. 95 and 0. 96) than that of 2ch-TDS+ for T2 lesions (ICC = 0. 82). The 3ch-TDS+ also identified the three lesion subtypes with high specificity (98. 7–99. 9%) and accuracy (98. 5–99. 9%). Sensitivity of 3ch-TDS+ for T2 lesions was 16% higher than with 2ch-TDS+. Enhancing lesions were segmented with the best sensitivity (81. 9%). “Black holes” were segmented with the least sensitivity (62. 3%). Conclusion. 3ch-TDS+ is a promising method for automated segmentation of MS lesion subtypes.

YNIMG Journal 2005 Journal Article

Changes in brain electrical activity during extended continuous word recognition

  • Jan W. Van Strien
  • Rogier E. Hagenbeek
  • Cornelis J. Stam
  • Serge A.R.B. Rombouts
  • Frederik Barkhof

Twenty healthy subjects (10 men, 10 women) participated in an EEG study with an extended continuous recognition memory task, in which each of 30 words was randomly shown 10 times and subjects were required to make old vs. new decisions. Both event-related brain potentials (ERPs) and induced band power (IBP) were investigated. We hypothesized that repeated presentations affect recollection rather than familiarity. For the 300- to 500-ms time window, an ‘old/new’ ERP effect was found for the first vs. second word presentations. The correct recognition of an ‘old’ word was associated with a more positive waveform than the correct identification of a new word. The old/new effect was most pronounced at and around the midline parietal electrode position. For the 500- to 800-ms time window, a linear repetition effect was found for multiple word repetitions. Correct recognition after an increasing number of repetitions was associated with increasing positivity. The multiple repetitions effect was most pronounced at the midline central (Cz) and fronto-central (FCz) electrode positions and reflects a graded recollection process: the stronger the memory trace grows, the more positive the ERP in the 500- to 800-ms time window. The ERP results support a dual-processing model, with familiarity being discernable from a more graded recollection state that depends on memory strengths. For IBP, we found ‘old/new’ effects for the lower-2 alpha, theta, and delta bands, with higher bandpower during ‘old’ words. The lower-2 alpha ‘old/new’ effect most probably reflects attentional processes, whereas the theta and delta effects reflect encoding and retrieval processes. Upon repeated word presentations, the magnitude of induced delta power in the 375- to 750-ms time window diminished linearly. Correlation analysis suggests that decreased delta power is moderately associated with faster decision speed and higher accuracy.

YNIMG Journal 2005 Journal Article

Delayed rather than decreased BOLD response as a marker for early Alzheimer's disease

  • Serge A.R.B. Rombouts
  • Rutger Goekoop
  • Cornelis J. Stam
  • Frederik Barkhof
  • Philip Scheltens

Functional MRI (fMRI) in established Alzheimer's disease (AD) shows regionally altered blood oxygenation level dependent (BOLD) responses. Mild cognitive impairment (MCI) is thought to represent an intermediate state between health and early Alzheimer's disease. To study this probable early dementia stage pathology, we studied in detail the BOLD response in MCI during visual encoding. 28 MCI patients, 18 AD patients, and 41 healthy elderly controls performed a face encoding task during fMRI scanning. Data were analyzed using orthogonal regressors, each representing different phases of the BOLD response (from slow to fast). Using a mixed effects model, regressor × group interactions were analyzed applying P < 0. 05, corrected. In occipital regions, MCI patients could be distinguished significantly better from controls and AD patients with a regressor of the early phase of the (fast) BOLD response than with the regressor of the late (slow) BOLD phase. Occipitally, the early phase BOLD response was significantly diminished in MCI patients compared to controls, and significantly increased when compared to AD. AD patients showed diminished early phase activation in widespread regions throughout the brain when compared to controls. There were no differences in the late (slow) phase of the BOLD response. This study stresses the importance of analyzing early phase BOLD responses and not only using one model of the BOLD response in neurodegenerative diseases. The increasing delay of the BOLD response from controls to MCI to AD may be consistent with the idea that MCI is a transitional state between healthy aging and dementia. Analyzing differences in different phases of the BOLD response introduces new opportunities to understand changes in regional brain dynamics in MCI and how well this may serve as an early marker of AD pathology.

YNIMG Journal 2005 Journal Article

Noradrenaline mediates amygdala activation in men and women during encoding of emotional material

  • Anda H. van Stegeren
  • Rutger Goekoop
  • Walter Everaerd
  • Philip Scheltens
  • Frederik Barkhof
  • Joost P.A. Kuijer
  • Serge A.R.B. Rombouts

The amygdala is a pivotal structure in humans for encoding of emotional information, as shown by recent imaging studies. It is unknown which neurotransmitters are specifically involved in the human amygdala, although in animal studies noradrenaline was shown to be essential. In our study, participants received the betablocker propranolol (which blocks the noradrenergic response) or placebo when watching neutral to highly negative arousing pictures. Amygdala activation, monitored with functional magnetic resonance imaging (fMRI), increased with emotional intensity of the pictures under placebo condition. Betablockade selectively decreased amygdala activation for emotional pictures of the second highest category, but not for the highest or lower (neutral) category pictures. Two findings add to the existing knowledge in this area. First, the activation pattern in the amygdala under placebo condition shows a nonlinearity related to the emotional categories of the pictures. Second, propranolol disturbs this activation pattern in the amygdala. Explorations with respect to gender show a similar effect of betablockade on amygdala activation in both men and women, but a difference in its effect on long-term memory for emotional pictures. This study supports the hypothesis that the neurotransmitter noradrenaline also mediates amygdala activity in humans when processing emotional stimuli and that betablockers can disrupt the normal activation pattern in the amygdala.

YNIMG Journal 2004 Journal Article

Challenging the cholinergic system in mild cognitive impairment: a pharmacological fMRI study

  • Rutger Goekoop
  • Serge A.R.B. Rombouts
  • Cees Jonker
  • Anneke Hibbel
  • Dirk L. Knol
  • Luc Truyen
  • Frederik Barkhof
  • Philip Scheltens

Mild cognitive impairment (MCI) often represents an early form of Alzheimer disease (AD). In both MCI and AD, characteristic cholinergic changes may occur. Functional magnetic resonance imaging (fMRI) may help to examine neurochemical changes in early disease by studying signal reactivity to pharmacological challenge. In this study, MCI patients [n = 28; mean age 73. 6 ± 7. 5; mini mental state examination (MMSE) 27. 0 ± 1. 2] were scanned during task performance in a randomized trial under three different medication regimes: at baseline [BL; no galantamine (GAL)], after a single oral dose of GAL (SD), and after prolonged exposure (steady state: SS). Memory tasks included an episodic face-encoding task and a parametric n-letter back working memory (WM) task. Alterations in brain activation patterns before and after treatment were analyzed for both tasks using multilevel statistical analysis. Significant increases in brain activation from BL were observed after prolonged exposure only. For face encoding (n = 28), these involved left prefrontal areas, the anterior cingulate gyrus, left occipital areas, and left posterior hippocampus. For working memory (n = 28), increased activation was found in right precuneus and right middle frontal gyrus, coinciding with increased accuracy scores after GAL treatment. In conclusion, cholinergic challenge produces alterations in brain activation patterns in elderly MCI patients that can be detected with fMRI. This should encourage further functional imaging studies to examine the status of neurotransmitter systems in disease.

YNIMG Journal 2003 Journal Article

Frontostriatal system in planning complexity: a parametric functional magnetic resonance version of tower of london task

  • Odile A van den Heuvel
  • Henk J Groenewegen
  • Frederik Barkhof
  • Richard H.C Lazeron
  • Richard van Dyck
  • Dick J Veltman

In the present study, we sought to investigate which brain structures are recruited in planning tasks of increasing complexity. For this purpose, a parametric self-paced pseudo-randomized event-related functional MRI version of the Tower of London task was designed. We tested 22 healthy subjects, enabling assessment of imaging results at a second (random effects) level of analysis. Compared with baseline, planning activity was correlated with increased blood oxygenation level-dependent (BOLD) signal in the dorsolateral prefrontal cortex, striatum, premotor cortex, supplementary motor area, and visuospatial system (precuneus and inferior parietal cortex). Task load was associated with increased activity in these same regions. In addition, increasing task complexity was correlated with activity in the left anterior prefrontal cortex, a region supposed to be specifically involved in third-order higher cognitive functioning.

v2026.09.13