MR fingerprinting (MRF) is a promising method for quantitative characterization of tissues. Often, voxel-wise measurements are made, assuming a single tissue-type per voxel. Alternatively, the... Show moreMR fingerprinting (MRF) is a promising method for quantitative characterization of tissues. Often, voxel-wise measurements are made, assuming a single tissue-type per voxel. Alternatively, the Sparsity Promoting Iterative Joint Non-negative least squares Multi-Component MRF method (SPIJN-MRF) facilitates tissue parameter estima-tion for identified components as well as partial volume segmentations. The aim of this paper was to evaluate the accuracy and repeatability of the SPIJN-MRF parameter estimations and partial volume segmentations. This was done (1) through numerical simulations based on the BrainWeb phantoms and (2) using in vivo acquired MRF data from 5 subjects that were scanned on the same week-day for 8 consecutive weeks. The partial volume segmen-tations of the SPIJN-MRF method were compared to those obtained by two conventional methods: SPM12 and FSL. SPIJN-MRF showed higher accuracy in simulations in comparison to FSL-and SPM12-based segmentations: Fuzzy Tanimoto Coefficients (FTC) comparing these segmentations and Brainweb references were higher than 0.95 for SPIJN-MRF in all the tissues and between 0.6 and 0.7 for SPM12 and FSL in white and gray matter and between 0.5 and 0.6 in CSF. For the in vivo MRF data, the estimated relaxation times were in line with literature and minimal variation was observed. Furthermore, the coefficient of variation (CoV) for estimated tissue volumes with SPIJN-MRF were 10.5% for the myelin water, 6.0% for the white matter, 5.6% for the gray matter, 4.6% for the CSF and 1.1% for the total brain volume. CoVs for CSF and total brain volume measured on the scanned data for SPIJN-MRF were in line with those obtained with SPM12 and FSL. The CoVs for white and gray mat-ter volumes were distinctively higher for SPIJN-MRF than those measured with SPM12 and FSL. In conclusion, the use of SPIJN-MRF provides accurate and precise tissue relaxation parameter estimations taking into account intrinsic partial volume effects. It facilitates obtaining tissue fraction maps of prevalent tissues including myelin water which can be relevant for evaluating diseases affecting the white matter. Show less
Results of neuroimaging datasets aggregated from multiple sites may be biased by site-specific profiles in participants' demographic and clinical characteristics, as well as MRI acquisition... Show moreResults of neuroimaging datasets aggregated from multiple sites may be biased by site-specific profiles in participants' demographic and clinical characteristics, as well as MRI acquisition protocols and scanning platforms. We compared the impact of four different harmonization methods on results obtained from analyses of cortical thickness data: (1) linear mixed-effects model (LME) that models site-specific random intercepts (LME INT), (2) LME that models both site-specific random intercepts and age-related random slopes (LME INT+ SLP), (3) ComBat, and (4) ComBat with a generalized additive model (ComBat-GAM). Our test case for comparing harmonization methods was cortical thickness data aggregated from 29 sites, which included 1,340 cases with posttraumatic stress disorder (PTSD) (6.2-81.8 years old) and 2,057 trauma-exposed controls without PTSD (6.3-85.2 years old). We found that, compared to the other data harmonization methods, data processed with ComBat-GAM was more sensitive to the detection of significant case-control differences (X-2 (3) = 63.704, p < 0.001) as well as casecontrol differences in age-related cortical thinning (X-2 (3) = 12.082, p = 0.007). Both ComBat and ComBat-GAM outperformed LME methods in detecting sex differences (X-2 (3) = 9.114, p = 0.028) in regional cortical thickness. ComBat-GAM also led to stronger estimates of age-related declines in cortical thickness (corrected p-values < 0.001), stronger estimates of case-related cortical thickness reduction (corrected p-values < 0.001), weaker estimates of age-related declines in cortical thickness in cases than controls (corrected p-values < 0.001), stronger estimates of cortical thickness reduction in females than males (corrected p-values < 0.001), and stronger estimates of cortical thickness reduction in females relative to males in cases than controls (corrected p-values < 0.001). Our results support the use of ComBat-GAM to minimize confounds and increase statistical power when harmonizing data with non-linear effects, and the use of either ComBat or ComBat-GAM for harmonizing data with linear effects. Show less
Hesse, L.S.; Aliasi, M.; Moser, F.; Haak, M.C.; Xie, W.D.; Jenkinson, M.; ... ; INTERGROWTH-21st Consortium 2022
The quantification of subcortical volume development from 3D fetal ultrasound can provide important diagnostic information during pregnancy monitoring. However, manual segmentation of subcortical... Show moreThe quantification of subcortical volume development from 3D fetal ultrasound can provide important diagnostic information during pregnancy monitoring. However, manual segmentation of subcortical structures in ultrasound volumes is time-consuming and challenging due to low soft tissue contrast, speckle and shadowing artifacts. For this reason, we developed a convolutional neural network (CNN) for the automated segmentation of the choroid plexus (CP), lateral posterior ventricle horns (LPVH), cavum septum pellucidum et vergae (CSPV), and cerebellum (CB) from 3D ultrasound. As ground-truth labels are scarce and expensive to obtain, we applied few-shot learning, in which only a small number of manual annotations (n = 9) are used to train a CNN. We compared training a CNN with only a few individually annotated volumes versus many weakly labelled volumes obtained from atlas-based segmentations. This showed that segmentation performance close to intra-observer variability can be obtained with only a handful of manual annotations. Finally, the trained models were applied to a large number (n = 278) of ultrasound image volumes of a diverse, healthy population, obtaining novel US-specific growth curves of the respective structures during the second trimester of gestation. Show less
Aims: Non-invasive measures of brain iron content would be of great benefit in neurodegeneration with brain iron accumulation (NBIA) to serve as a biomarker for disease progression and evaluation... Show moreAims: Non-invasive measures of brain iron content would be of great benefit in neurodegeneration with brain iron accumulation (NBIA) to serve as a biomarker for disease progression and evaluation of iron chelation therapy. Although magnetic resonance imaging (MRI) provides several quantitative measures of brain iron content, none of these have been validated for patients with a severely increased cerebral iron burden. We aimed to validate R 2 * as a quantitative measure of brain iron content in aceruloplasminemia, the most severely iron-loaded NBIA phenotype. Methods: Tissue samples from 50 gray-and white matter regions of a postmortem aceruloplasminemia brain and control subject were scanned at 1.5 T to obtain R 2 * , and biochemically analyzed with inductively coupled plasma mass spectrometry. For gray matter samples of the aceruloplasminemia brain, sample R 2 * values were compared with postmortem in situ MRI data that had been obtained from the same subject at 3 T - in situ R 2 * . Relationships between R 2 * and tissue iron concentration were determined by linear regression analyses. Results: Median iron concentrations throughout the whole aceruloplasminemia brain were 10 to 15 times higher than in the control subject, and R 2 * was linearly associated with iron concentration. For gray matter samples of the aceruloplasminemia subject with an iron concentration up to 1000 mg/kg, 91% of variation in R 2 * could be explained by iron, and in situ R 2 * at 3 T and sample R 2 * at 1.5 T were highly correlated. For white matter regions of the aceruloplasminemia brain, 85% of variation in R 2 * could be explained by iron. Conclusions: R 2 * is highly sensitive to variations in iron concentration in the severely iron-loaded brain, and might be used as a non-invasive measure of brain iron content in aceruloplasminemia and potentially other NBIA disorders. Show less
Kat, R.; Berg, B. van den; Perenboom, M.J.L.; Schenke, M.; Maagdenberg, A.M.J.M. van den; Bruining, H.; ... ; Kas, M.J.H. 2021
The mouse is widely used as an experimental model to study visual processing. To probe how the visual system detects changes in the environment, functional paradigms in freely behaving mice are... Show moreThe mouse is widely used as an experimental model to study visual processing. To probe how the visual system detects changes in the environment, functional paradigms in freely behaving mice are strongly needed. We developed and validated the first EEG-based method to investigate visual deviance detection in freely behaving mice. Mice with EEG implants were exposed to a visual deviant detection paradigm that involved changes in light intensity as standard and deviant stimuli. By subtracting the standard from the deviant evoked waveform, deviant detection was evident as bi-phasic negativity (starting around 70 ms) in the difference waveform. Additionally, deviance-associated evoked (beta/gamma) and induced (gamma) oscillatory responses were found. We showed that the results were stimulus-independent by applying a "flip-flop " design and the results showed good repeatability in an independent measurement. Together, we put forward a validated, easy-to-use paradigm to measure visual deviance processing in freely behaving mice. Show less
Petitclerc, L.; Hirschler, L.; Wells, J.A.; Thomas, D.L.; Walderveen, M.A.A. van; Buchem, M.A. van; Osch, M.J.P. van 2021
The study of brain clearance mechanisms is an active area of research. While we know that the cerebrospinal fluid (CSF) plays a central role in one of the main existing clearance pathways, the... Show moreThe study of brain clearance mechanisms is an active area of research. While we know that the cerebrospinal fluid (CSF) plays a central role in one of the main existing clearance pathways, the exact processes for the secretion of CSF and the removal of waste products from tissue are under debate. CSF is thought to be created by the exchange of water and ions from the blood, which is believed to mainly occur in the choroid plexus. This exchange has not been thoroughly studied in vivo.We propose a modified arterial spin labeling (ASL) MRI sequence and image analysis to track blood water as it is transported to the CSF, and to characterize its exchange from blood to CSF. We acquired six pseudo-continuous ASL sequences with varying labeling duration (LD) and post-labeling delay (PLD) and a segmented 3D-GRASE readout with a long echo train (8 echo times (TE)) which allowed separation of the very long-T-2 CSF signal. ASL signal was observed at long TEs (793 ms and higher), indicating presence of labeled water transported from blood to CSF. This signal appeared both in the CSF proximal to the choroid plexus and in the subarachnoid space surrounding the cortex. ASL signal was separated into its blood, gray matter and CSF components by fitting a triexponential function with T(2)s taken from literature. A two-compartment dynamic model was introduced to describe the exchange of water through time and TE. From this, a water exchange time from the blood to the CSF (Tbl->CSF) was mapped, with an order of magnitude of approximately 60 s. Show less
Chappell, M.A.; McConnell, F.A.K.; Golay, X.; Gunther, M.; Hernandez-Tamames, J.A.; Osch, M.J. van; Asllani, I. 2021
The mismatch in the spatial resolution of Arterial Spin Labeling (ASL) MRI perfusion images and the anatomy of functionally distinct tissues in the brain leads to a partial volume effect (PVE),... Show moreThe mismatch in the spatial resolution of Arterial Spin Labeling (ASL) MRI perfusion images and the anatomy of functionally distinct tissues in the brain leads to a partial volume effect (PVE), which in turn confounds the estimation of perfusion into a specific tissue of interest such as gray or white matter. This confound occurs because the image voxels contain a mixture of tissues with disparate perfusion properties, leading to estimated perfusion values that reflect primarily the volume proportions of tissues in the voxel rather than the perfusion of any particular tissue of interest within that volume. It is already recognized that PVE influences studies of brain perfusion, and that its effect might be even more evident in studies where changes in perfusion are co-incident with alterations in brain structure, such as studies involving a comparison between an atrophic patient population vs control subjects, or studies comparing subjects over a wide range of ages. However, the application of PVE correction (PVEc) is currently limited and the employed methodologies remain inconsistent. In this article, we outline the influence of PVE in ASL measurements of perfusion, explain the main principles of PVEc, and provide a critique of the current state of the art for the use of such methods. Furthermore, we examine the current use of PVEc in perfusion studies and whether there is evidence to support its wider adoption. We conclude that there is sound theoretical motivation for the use of PVEc alongside conventional, 'uncorrected', images, and encourage such combined reporting. Methods for PVEc are now available within standard neuroimaging toolboxes, which makes our recommendation straightforward to implement. However, there is still more work to be done to establish the value of PVEc as well as the efficacy and robustness of existing PVEc methods. Show less
Lundell, H.; Najac, C.; Bulk, M.; Kan, H.E.; Webb, A.G.; Ronen, I. 2021
Double diffusion encoding (DDE) of the water signal offers a unique ability to separate the effect of microscopic anisotropic diffusion in structural units of tissue from the overall macroscopic... Show moreDouble diffusion encoding (DDE) of the water signal offers a unique ability to separate the effect of microscopic anisotropic diffusion in structural units of tissue from the overall macroscopic orientational distribution of cells. However, the specificity in detected microscopic anisotropy is limited as the signal is averaged over different cell types and across tissue compartments. Performing side-by-side water and metabolite DDE spectroscopic (DDES) experiments provides complementary measures from which intracellular and extracellular microscopic fractional anisotropies (mu FA) and diffusivities can be estimated. Metabolites are largely confined to the intracellular space and therefore provide a benchmark for intracellular mu FA and diffusivities of specific cell types. By contrast, water DDES measurements allow examination of the separate contributions to water mu FA and diffusivity from the intraand extracellular spaces, by using a wide range of b values to gradually eliminate the extracellular contribution. Here, we aimed to estimate tissue and compartment specific human brain microstructure by combining water and metabolites DDES experiments. We performed our DDES measurements in two brain regions that contain widely different amounts of white matter (WM) and gray matter (GM): parietal white matter (PWM) and occipital gray matter (OGM) in a total of 20 healthy volunteers at 7 Tesla. Metabolite DDES measurements were performed at b = 7199 s/mm(2), while water DDES measurements were performed with a range of b values from 918 to 7199 s/mm(2). The experimental framework we employed here resulted in a set of insights pertaining to the morphology of the intracellular and extracellular spaces in both gray and white matter. Results of the metabolite DDES experiments in both PWM and OGM suggest a highly anisotropic intracellular space within neurons and glia, with the possible exception of gray matter glia. The water mu FA obtained from the DDES results at high b values in both regions converged with that of the metabolite DDES, suggesting that the signal from the extracellular space is indeed effectively suppressed at the highest b value. The mu FA measured in the OGM significantly decreased at lower b values, suggesting a considerably lower anisotropy of the extracellular space in GM compared to WM. In PWM, the water mu FA remained high even at the lowest b value, indicating a high degree of organization in the interstitial space in WM. Tortuosity values in the cytoplasm for water and tNAA, obtained with correlation analysis of microscopic parallel diffusivity with respect to GM/WM tissue fraction in the volume of interest, are remarkably similar for both molecules, while exhibiting a clear difference between gray and white matter, suggesting a more crowded cytoplasm and more complex cytomorphology of neuronal cell bodies and dendrites in GM than those found in long-range axons in WM. Show less
Background: Depression has been associated with decreased regional grey matter volume, which might partly be explained by an unhealthier lifestyle in depressed individuals which has been ignored by... Show moreBackground: Depression has been associated with decreased regional grey matter volume, which might partly be explained by an unhealthier lifestyle in depressed individuals which has been ignored by most earlier studies. Also, the longitudinal nature of depression, lifestyle and brain structure associations is largely unknown. This study investigates the relationship of depression and lifestyle with brain structure cross-sectionally and longitudinally over up to 9 years.Methods: We used longitudinal structural MRI data of persons with depression and/or anxiety disorders and controls (N-unique participants = 347, N-observations = 609). Cortical thickness of medial orbitofrontal cortex (mOFC), rostral anterior cingulate cortex (rACC) and hippocampal volume were derived using FreeSurfer. Using Generalized Estimating Equations, we investigated associations of depression and lifestyle (Body mass index (BMI), smoking, alcohol consumption, physical activity and sleep duration) with brain structure and change in brain structure over 2 (n = 179) and 9 years (n = 82).Results: Depression status (B =-.053, p = .002) and severity (B =-.002, p = .002) were negatively associated with rACC thickness. mOFC thickness was negatively associated with BMI (B =-.004, p < .001) and positively with moderate alcohol consumption (B = .030, p = .009). All associations were independent of each other. No associations were observed between (change in) depression, disease burden or lifestyle factors with brain change over time.Conclusions: Depressive symptoms and diagnosis were independently associated with thinner rACC, BMI with thinner mOFC, and moderate alcohol consumption with thicker mOFC. No longitudinal associations were observed, suggesting that regional grey matter alterations are a long-term consequence or vulnerability indicator for depression but not dynamically or progressively related to depression course trajectory. Show less
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... Show moreArterial 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. Show less
Nagtegaal, M.; Koken, P.; Amthor, T.; Bresser, J. de; Madler, B.; Vos, F.; Doneva, M. 2020
Demyelination is the key pathological process in multiple sclerosis (MS). The extent of demyelination can be quantified with magnetic resonance imaging by assessing the myelin water fraction (MWF).... Show moreDemyelination is the key pathological process in multiple sclerosis (MS). The extent of demyelination can be quantified with magnetic resonance imaging by assessing the myelin water fraction (MWF). However, long computation times and high noise sensitivity hinder the translation of MWF imaging to clinical practice. In this work, we introduce a more efficient and noise robust method to determine the MWF using a joint sparsity constraint and a pre-computed B-1(+)-T-2 dictionary.A single component analysis with this dictionary is used in an initial step to obtain a B-1(+) map. The T-2 distribution is then determined from a reduced dictionary corresponding to the estimated B-1(+) map using a combination of a non-negativity and a joint sparsity constraint.The non-negativity constraint ensures that a feasible solution with non-negative contribution of each T-2 component is obtained. The joint sparsity constraint restricts the T-2 distribution to a small set of T-2 relaxation times shared between all voxels and reduces the noise sensitivity.The applied Sparsity Promoting Iterative Joint NNLS (SPIJN) algorithm can be implemented efficiently, reducing the computation time by a factor of 50 compared to the commonly used regularized non-negative least squares algorithm. The proposed method was validated in simulations and in 8 healthy subjects with a 3D multiecho gradient- and spin echo scan at 3 T. In simulations, the absolute error in the MWF decreased from 0.031 to 0.013 compared to the regularized NNLS algorithm for SNR = 250. The in vivo results were consistent with values reported in literature and improved MWF-quantification was obtained especially in the frontal white matter. The maximum standard deviation in mean MWF in different regions of interest between subjects was smaller for the proposed method (0.0193) compared to the regularized NNLS algorithm (0.0266). In conclusion, the proposed method for MWF estimation is less computationally expensive and less susceptible to noise compared to state of the art methods. These improvements might be an important step towards clinical translation of MWF measurements. Show less
Bulk, M.; Abdelmoula, W.M.; Geut, H.; Wiarda, W.; Ronen, I.; Dijkstra, J.; Weerd, L. van der 2020
Disruption of cholinergic and serotonergic neurotransmitter systems is associated with cognitive, emotional and behavioural symptoms of Alzheimer's disease (AD). To investigate the responsiveness... Show moreDisruption of cholinergic and serotonergic neurotransmitter systems is associated with cognitive, emotional and behavioural symptoms of Alzheimer's disease (AD). To investigate the responsiveness of these systems in AD we measured the effects of a single-dose of the selective serotonin reuptake inhibitor citalopram and acetylcholinesterase inhibitor galantamine in 12 patients with AD and 12 age-matched controls on functional brain connectivity with resting state functional magnetic resonance imaging. In this randomized, double blind, placebo-controlled crossover study, functional magnetic resonance images were repeatedly obtained before and after dosing, resulting in a dataset of 432 scans. Connectivity maps of ten functional networks were extracted using a dual regression method and drug vs. placebo effects were compared between groups with a multivariate analysis with signals coming from cerebrospinal fluid and white matter as covariates at the subject level, and baseline and heart rate measurements as confound regressors in the higher-level analysis (at p < 0.05, corrected). A galantamine induced difference between groups was observed for the cerebellar network. Connectivity within the cerebellar network and between this network and the thalamus decreased after galantamine vs. placebo in AD patients, but not in controls. For citalopram, voxelwise network connectivity did not show significant group x treatment interaction effects. However, we found default mode network connectivity with the precuneus and posterior cingulate cortex to be increased in AD patients, which could not be detected within the control group. Further, in contrast to the AD patients, control subjects showed a consistent reduction in mean connectivity with all networks after administration of citalopram. Since AD has previously been characterized by reduced connectivity between the default mode network and the precuneus and posterior cingulate cortex, the effects of citalopram on the default mode network suggest a restoring potential of selective serotonin reuptake inhibitors in AD. The results of this study also confirm a change in cerebellar connections in AD, which is possibly related to cholinergic decline. Show less
Klaassens, B.L.; Gerven, J.M.A. van; Klaassen, E.S.; Grond, J. van der; Rombouts, S.A.R.B. 2019
Performing quality control to detect image artifacts and data-processing errors is crucial in structural magnetic resonance imaging, especially in developmental studies. Currently, many studies... Show morePerforming quality control to detect image artifacts and data-processing errors is crucial in structural magnetic resonance imaging, especially in developmental studies. Currently, many studies rely on visual inspection by trained raters for quality control. The subjectivity of these manual procedures lessens comparability between studies, and with growing study sizes quality control is increasingly time consuming. In addition, both inter-rater as well as intra-rater variability of manual quality control is high and may lead to inclusion of poor quality scans and exclusion of scans of usable quality. In the current study we present the Qoala-T tool, which is an easy and free to use supervised-learning model to reduce rater bias and misclassification in manual quality control procedures using FreeSurfer-processed scans. First, we manually rated quality of N = 784 FreeSurfer-processed T1-weighted scans acquired in three different waves in a longitudinal study. Different supervised-learning models were then compared to predict manual quality ratings using FreeSurfer segmented output data. Results show that the Qoala-T tool using random forests is able to predict scan quality with both high sensitivity and specificity (mean area under the curve (AUC) = 0.98). In addition, the Qoala-T tool was also able to adequately predict the quality of two novel unseen datasets (total N = 872). Finally, analyses of age effects showed that younger participants were more likely to have lower scan quality, underlining that scan quality might confound findings attributed to age effects. These outcomes indicate that this procedure could further help to reduce variability related to manual quality control, thereby benefiting the comparability of data quality between studies. Show less
