Region of interest correction factors improve reliability of diffusion imaging measures within and across scanners and field strengths

Characterizing patterns of DTI variance in aging brains

Background

As large analyses merge data across sites, a deeper understanding of variance in statistical assessment across the sources of data becomes critical for valid analyses. Diffusion tensor imaging (DTI) exhibits spatially varying and correlated noise, so care must be taken with distributional assumptions.

Purpose

We characterize the role of physiology, subject compliance, and the interaction of subject with the scanner in the understanding of DTI variability, as modeled in spatial variance of derived metrics in homogeneous regions.

Methods

We analyze DTI data from 1035 subjects in the Baltimore Longitudinal Study of Aging (BLSA), with ages ranging from 22.4 to 103 years old. For each subject, up to 12 longitudinal sessions were conducted. We assess variance of DTI scalars within regions of interest (ROIs) defined by four segmentation methods and investigate the relationships between the variance and covariates, including baseline age, time from the baseline (referred to as "interval"), motion, sex, and whether it is the first scan or the second scan in the session.

Results

Covariate effects are heterogeneous and bilaterally symmetric across ROIs. Inter-session interval is positively related (p ≪ 0.001) to FA variance in the cuneus and occipital gyrus, but negatively (p ≪ 0.001) in the caudate nucleus. Males show significantly (p ≪ 0.001) higher FA variance in the right putamen, thalamus, body of the corpus callosum, and cingulate gyrus. In 62 out of 176 ROIs defined by the Eve type-1 atlas, an increase in motion is associated (p < 0.05) with a decrease in FA variance. Head motion increases during the rescan of DTI (Δμ = 0.045 millimeters per volume).

Conclusions

The effects of each covariate on DTI variance, and their relationships across ROIs are complex. Ultimately, we encourage researchers to include estimates of variance when sharing data and consider models of heteroscedasticity in analysis. This work provides a foundation for study planning to account for regional variations in metric variance.

Removal of site effects and enhancement of signal using dual projection independent component analysis for pooling multi-site MRI data

Combining magnetic resonance imaging (MRI) data from multi-site studies is a popular approach for constructing larger datasets to greatly enhance the reliability and reproducibility of neuroscience research. However, the scanner/site variability is a significant confound that complicates the interpretation of the results, so effective and complete removal of the scanner/site variability is necessary to realise the full advantages of pooling multi-site datasets. Independent component analysis (ICA) and general linear model (GLM) based harmonisation methods are the two primary methods used to eliminate scanner/site effects. Unfortunately, there are challenges with both ICA-based and GLM-based harmonisation methods to remove site effects completely when the signals of interest and scanner/site effects-related variables are correlated, which may occur in neuroscience studies. In this study, we propose an effective and powerful harmonisation strategy that implements dual projection (DP) theory based on ICA to remove the scanner/site effects more completely. This method can separate the signal effects correlated with site variables from the identified site effects for removal without losing signals of interest. Both simulations and vivo structural MRI datasets, including a dataset from Autism Brain Imaging Data Exchange II and a travelling subject dataset from the Strategic Research Program for Brain Sciences, were used to test the performance of a DP-based ICA harmonisation method. Results show that DP-based ICA harmonisation has superior performance for removing site effects and enhancing the sensitivity to detect signals of interest as compared with GLM-based and conventional ICA harmonisation methods.

A Brain Morphometry Study with Across-Site Harmonization Using a ComBat-Generalized Additive Model in Children and Adolescents

Regional anatomical structures of the brain are intimately connected to functions corresponding to specific regions and the temporospatial pattern of genetic expression and their functions from the fetal period to old age. Therefore, quantitative brain morphometry has often been employed in neuroscience investigations, while controlling for the scanner effect of the scanner is a critical issue for ensuring accuracy in brain morphometric studies of rare orphan diseases due to the lack of normal reference values available for multicenter studies. This study aimed to provide across-site normal reference values of global and regional brain volumes for each sex and age group in children and adolescents. We collected magnetic resonance imaging (MRI) examinations of 846 neurotypical participants aged 6.0-17.9 years (339 male and 507 female participants) from 5 institutions comprising healthy volunteers or neurotypical patients without neurological disorders, neuropsychological disorders, or epilepsy. Regional-based analysis using the CIVET 2.1.0. pipeline provided regional brain volumes, and the measurements were across-site combined using ComBat-GAM harmonization. The normal reference values of global and regional brain volumes and lateral indices in our study could be helpful for evaluating the characteristics of the brain morphology of each individual in a clinical setting and investigating the brain morphology of ultra-rare diseases.

ComBat Harmonization: Empirical Bayes versus fully Bayes approaches

Studying small effects or subtle neuroanatomical variation requires large-scale sample size data. As a result, combining neuroimaging data from multiple datasets is necessary. Variation in acquisition protocols, magnetic field strength, scanner build, and many other non-biologically related factors can introduce undesirable bias into studies. Hence, harmonization is required to remove the bias-inducing factors from the data. ComBat is one of the most common methods applied to features from structural images. ComBat models the data using a hierarchical Bayesian model and uses the empirical Bayes approach to infer the distribution of the unknown factors. The empirical Bayes harmonization method is computationally efficient and provides valid point estimates. However, it tends to underestimate uncertainty. This paper investigates a new approach, fully Bayesian ComBat, where Monte Carlo sampling is used for statistical inference. When comparing fully Bayesian and empirical Bayesian ComBat, we found Empirical Bayesian ComBat more effectively removed scanner strength information and was much more computationally efficient. Conversely, fully Bayesian ComBat better preserved biological disease and age-related information while performing more accurate harmonization on traveling subjects. The fully Bayesian approach generates a rich posterior distribution, which is useful for generating simulated imaging features for improving classifier performance in a limited data setting. We show the generative capacity of our model for augmenting and improving the detection of patients with Alzheimer's disease. Posterior distributions for harmonized imaging measures can also be used for brain-wide uncertainty comparison and more principled downstream statistical analysis.Code for our new fully Bayesian ComBat extension is available at https://github.com/batmanlab/BayesComBat.

Altered cerebellar-motor loop in benign adult familial myoclonic epilepsy type 1: The structural basis of cortical tremor

OBJECTIVE

Cortical tremor/myoclonus is the hallmark feature of benign adult familial myoclonic epilepsy (BAFME), the mechanism of which remains elusive. A hypothesis is that a defective control in the preexisting cerebellar-motor loop drives cortical tremor. Meanwhile, the basal ganglia system might also participate in BAFME. This study aimed to discover the structural basis of cortical tremor/myoclonus in BAFME.

METHODS

Nineteen patients with BAFME type 1 (BAFME1) and 30 matched healthy controls underwent T1-weighted and diffusion tensor imaging scans. FreeSurfer and spatially unbiased infratentorial template (SUIT) toolboxes were utilized to assess the motor cortex and the cerebellum. Probabilistic tractography was generated for two fibers to test the hypothesis: the dentato-thalamo-(M1) (primary motor cortex) and globus pallidus internus (GPi)-thalamic projections. Average fractional anisotropy (FA), axial diffusivity (AD), mean diffusivity (MD), and radial diffusivity (RD) of each tract were extracted.

RESULTS

Cerebellar atrophy and dentate nucleus alteration were observed in the patients. In addition, patients with BAFME1 exhibited reduced AD and FA in the left and right dentato-thalamo-M1 nondecussating fibers, respectively false discovery rate (FDR) correction q < .05. Cerebellar projections showed negative correlations with somatosensory-evoked potential P25-N33 amplitude and were independent of disease duration and medication. BAFME1 patients also had increased FA and decreased MD in the left GPi-thalamic projection. Higher FA and lower RD in the right GPi-thalamic projection were also observed (FDR q < .05).

SIGNIFICANCE

The present findings support the hypothesis that the cerebello-thalamo-M1 loop might be the structural basis of cortical tremor in BAFME1. The basal ganglia system also participates in BAFME1 and probably serves a regulatory role.

Longitudinal monitoring of amyotrophic lateral sclerosis by diffusion tensor imaging: Power calculations for group studies

Introduction

Diffusion tensor imaging (DTI) can be used to map disease progression in amyotrophic lateral sclerosis (ALS) and therefore is a promising candidate for a biomarker in ALS. To this end, longitudinal study protocols need to be optimized and validated regarding group sizes and time intervals between visits. The objective of this study was to assess the influences of sample size, the schedule of follow-up measurements, and measurement uncertainties on the statistical power to optimize longitudinal DTI study protocols in ALS.

Patients and methods

To estimate the measurement uncertainty of a tract-of–interest-based DTI approach, longitudinal test-retest measurements were applied first to a normal data set. Then, DTI data sets of 80 patients with ALS and 50 healthy participants were analyzed in the simulation of longitudinal trajectories, that is, longitudinal fractional anisotropy (FA) values for follow-up sessions were simulated for synthetic patient and control groups with different rates of FA decrease in the corticospinal tract. Monte Carlo simulations of synthetic longitudinal study groups were used to estimate the statistical power and thus the potentially needed sample sizes for a various number of scans at one visit, different time intervals between baseline and follow-up measurements, and measurement uncertainties.

Results

From the simulation for different longitudinal FA decrease rates, it was found that two scans per session increased the statistical power in the investigated settings unless sample sizes were sufficiently large and time intervals were appropriately long. The positive effect of a second scan per session on the statistical power was particularly pronounced for FA values with high measurement uncertainty, for which the third scan per session increased the statistical power even further.

Conclusion

With more than one scan per session, the statistical power of longitudinal DTI studies can be increased in patients with ALS. Consequently, sufficient statistical power can be achieved even with limited sample sizes. An improved longitudinal DTI study protocol contributes to the detection of small changes in diffusion metrics and thereby supports DTI as an applicable and reliable non-invasive biomarker in ALS.

Longitudinal corpus callosum microstructural decline in early-stage Parkinson’s disease in association with akinetic-rigid symptom severity

Previous diffusion tensor imaging (DTI) studies of Parkinson’s disease (PD) show reduced microstructural integrity of the corpus callosum (CC) relative to controls, although the characteristics of such callosal degradation remain poorly understood. Here, we utilized a longitudinal approach to identify microstructural decline in the entire volume of the CC and its functional subdivisions over 2 years and related the callosal changes to motor symptoms in early-stage PD. The study sample included 61 PD subjects (N = 61, aged 45–82, 38 M & 23 F, H&Y ≤ 2) from the Parkinson’s Progressive Markers Initiative database (PPMI). Whole-brain voxel-wise results revealed significant fractional anisotropy (FA) and mean diffusivity (MD) changes in the CC, especially in the genu and splenium. Using individually drawn CC regions of interest (ROI), our analysis further revealed that almost all subdivisions of the CC show significant decline in FA to certain extents over the two-year timeframe. Additionally, FA seemed lower in the right hemisphere of the CC at both time-points, and callosal FA decline was associated with FA and MD decline in widespread cortical and subcortical areas. Notably, multiple regression analysis revealed that across-subject akinetic-rigid severity was negatively associated with callosal FA at baseline and 24 months follow-up, and the effect was strongest in the anterior portion of the CC. These results suggest that callosal microstructure alterations in the anterior CC may serve as a viable biomarker for akinetic-rigid symptomology and disease progression, even in early PD.

Accelerated decline in white matter microstructure in subsequently impaired older adults and its relationship with cognitive decline

Little is known about a longitudinal decline in white matter microstructure and its associations with cognition in preclinical dementia. Longitudinal diffusion tensor imaging and neuropsychological testing were performed in 50 older adults who subsequently developed mild cognitive impairment or dementia (subsequently impaired) and 200 cognitively normal controls. Rates of white matter microstructural decline were compared between groups using voxel-wise linear mixed-effects models. Associations between change in white matter microstructure and cognition were examined. Subsequently impaired individuals had a faster decline in fractional anisotropy in the right inferior fronto-occipital fasciculus and bilateral splenium of the corpus callosum. A decline in right inferior fronto-occipital fasciculus fractional anisotropy was related to a decline in verbal memory, visuospatial ability, processing speed and mini-mental state examination. A decline in bilateral splenium fractional anisotropy was related to a decline in verbal fluency, processing speed and mini-mental state examination. Accelerated regional white matter microstructural decline is evident during the preclinical phase of mild cognitive impairment/dementia and is related to domain-specific cognitive decline.

Brain Imaging Biomarkers for Chronic Pain

The prevalence of chronic pain has reached epidemic levels. In addition to personal suffering chronic pain is associated with psychiatric and medical co-morbidities, notably substance misuse, and a huge a societal cost amounting to hundreds of billions of dollars annually in medical cost, lost wages, and productivity. Chronic pain does not have a cure or quantitative diagnostic or prognostic tools. In this manuscript we provide evidence that this situation is about to change. We first start by summarizing our current understanding of the role of the brain in the pathogenesis of chronic pain. We particularly focus on the concept of learning in the emergence of chronic pain, and the implication of the limbic brain circuitry and dopaminergic signaling, which underly emotional learning and decision making, in this process. Next, we summarize data from our labs and from other groups on the latest brain imaging findings in different chronic pain conditions focusing on results with significant potential for translation into clinical applications. The gaps in the study of chronic pain and brain imaging are highlighted in throughout the overview. Finally, we conclude by discussing the costs and benefits of using brain biomarkers of chronic pain and compare to other potential markers.

Batch effect reduction of microarray data with dependent samples using an empirical Bayes approach (BRIDGE)

Batch-effects present challenges in the analysis of high-throughput molecular data and are particularly problematic in longitudinal studies when interest lies in identifying genes/features whose expression changes over time, but time is confounded with batch. While many methods to correct for batch-effects exist, most assume independence across samples; an assumption that is unlikely to hold in longitudinal microarray studies. We propose Batch effect Reduction of mIcroarray data with Dependent samples usinGEmpirical Bayes (BRIDGE), a three-step parametric empirical Bayes approach that leverages technical replicate samples profiled at multiple timepoints/batches, so-called "bridge samples", to inform batch-effect reduction/attenuation in longitudinal microarray studies. Extensive simulation studies and an analysis of a real biological data set were conducted to benchmark the performance of BRIDGE against both ComBat and longitudinalComBat. Our results demonstrate that while all methods perform well in facilitating accurate estimates of time effects, BRIDGE outperforms both ComBat and longitudinal ComBat in the removal of batch-effects in data sets with bridging samples, and perhaps as a result, was observed to have improved statistical power for detecting genes with a time effect. BRIDGE demonstrated competitive performance in batch effect reduction of confounded longitudinal microarray studies, both in simulated and a real data sets, and may serve as a useful preprocessing method for researchers conducting longitudinal microarray studies that include bridging samples.

Age affects white matter microstructure and episodic memory across the older adult lifespan

Diffusion imaging studies have observed age-related degradation of white matter that contributes to cognitive deficits separately in younger-old (ages 65-89) and oldest-old (ages 90+) adults. But it remains unclear whether these age effects are magnified in advanced age groups, which may reflect disease-related pathology. Here, we tested whether age-related differences in white matter microstructure followed linear or nonlinear patterns across the entire older adult lifespan (65-98 years), these patterns were influenced by oldest-old adults at increased risk of dementia (cognitive impairment no dementia, CIND), and they explained age effects on episodic memory. Results revealed nonlinear microstructure declines across fiber classes (medial temporal, callosal, association, projection and/or thalamic) that were largest for medial temporal fibers. These patterns remained after excluding oldest-old participants with CIND, indicating that aging of white matter microstructure cannot solely be explained by pathology associated with early cognitive impairment. Moreover, finding that the effect of age on episodic memory was mediated by medial temporal fiber microstructure suggests it is essential for facilitating memory-related neural signals across the older adult lifespan.

Association of Poorer Hearing With Longitudinal Change in Cerebral White Matter Microstructure

Importance

There is a dearth of studies that examine the association between poorer hearing and change in cerebral white matter (WM) microstructure.

Objective

To examine the association of poorer hearing with baseline and change in WM microstructure among older adults.

Design, Setting, and Participants

This was a prospective cohort study that evaluated speech-in-noise, pure-tone audiometry, and WM microstructure, as measured by mean diffusivity (MD) and fractional anisotropy (FA), both of which were evaluated by diffusion tensor imaging (DTI) in 17 WM regions. Data were collected between October 2012 and December 2018 and analyzed between March 2019 and August 2019 with a mean follow-up time of 1.7 years. The study evaluated responses to the Baltimore Longitudinal Study of Aging among 356 cognitively normal adults who had at least 1 hearing assessment and DTI session. Excluded were those with baseline cognitive impairment, stroke, head injuries, Parkinson disease, and/or bipolar disorder.

Exposures

Peripheral auditory function was measured by pure-tone average in the better-hearing ear. Central auditory function was measured by signal-to-noise ratio score from a speech-in-noise task and adjusted by pure-tone average.

Main Outcomes and Measures

Linear mixed-effects models with random intercepts and slopes were used to examine the association of poorer peripheral and central auditory function with baseline and longitudinal DTI metrics in 17 WM regions, adjusting for baseline characteristics (age, sex, race, hypertension, elevated total cholesterol, and obesity).

Results

Of 356 cognitively normal adults included in the study, the mean (SD) age was 73.5 (8.8) years, and 204 (57.3%) were women. There were no baseline associations between hearing and DTI measures. Longitudinally, poorer peripheral hearing was associated with increases in MD in the inferior fronto-occipital fasciculus (β = 0.025; 95% CI, 0.008-0.042) and the body (β = 0.050; 95% CI, 0.015-0.085) of the corpus callosum, but there were no associations of peripheral hearing with FA changes in these tracts. Poorer central auditory function was associated with longitudinal MD increases (β = 0.031; 95% CI, 0.010-0.052) and FA declines (β = -1.624; 95% CI, -2.511 to -0.738) in the uncinate fasciculus.

Conclusions and Relevance

Findings of this cohort study suggest that poorer hearing is related to change in integrity of specific WM regions involved with auditory processing.

Interval-censored survival analysis of mild traumatic brain injury with outcome based neuroimaging clinical applications

Objective

The purpose of this study was to assess the relationship between MRI findings and clinical presentation and outcomes in patients following mild traumatic brain injury (mTBI). We hypothesize that imaging findings other than hemorrhages and contusions may be used to predict symptom presentation and longevity following mTBI.

Methods

Patients (n = 250) diagnosed with mTBI and in litigation for brain injury underwent 3T magnetic resonance imaging (MRI). A retrospective chart review was performed to assess symptom presentation and improvement/resolution. To account for variable times of clinical presentation, nonuniform follow-up, and uncertainty in the dates of symptom resolution, a right censored, interval censored statistical analysis was performed. Incidence and resolution of headache, balance, cognitive deficit, fatigue, anxiety, depression, and emotional lability were compared among patients. Image findings analyzed included white matter hyperintensities (WMH), Diffusion Tensor Imaging (DTI) fractional anisotropy (FA) values, MR perfusion, auditory functional MRI (fMRI) activation, hippocampal atrophy (HA) and hippocampal asymmetry as defined by NeuroQuant ® volumetric software.

Results

Patients who reported LOC were significantly more likely to present with balance problems (p < 0.001), cognitive deficits (p = 0.010), fatigue (p = 0.025), depression (p = 0.002), and emotional lability (p = 0.002). Patients with LOC also demonstrated significantly slower recovery of cognitive function than those who did not lose consciousness (p = 0.044). Patients over the age of 40 had significantly higher odds of presenting with balance problems (p = 0.006). Additionally, these older patients were slower to recover cognitive function (p = 0.001) and less likely to experience improvement of headaches (p = 0.007). Abnormal MRI did not correlate significantly with symptom presentation, but was a strong indicator of symptom progression, with slower recovery of balance (p = 0.009) and cognitive deficits (p < 0.001).

Conclusion

This analysis demonstrates the utility of clinical data analysis using interval-censored survival statistical technique in head trauma patients. Strong statistical associations between neuroimaging findings and aggregate clinical outcomes were identified in patients with mTBI.

A deep learning-based method for improving reliability of multicenter diffusion kurtosis imaging with varied acquisition protocols

Multicenter magnetic resonance imaging is gaining more popularity in large-sample projects. Since both varying hardware and software across different centers cause unavoidable data heterogeneity across centers, its impact on reliability in study outcomes has also drawn much attention recently. One fundamental issue arises in how to derive model parameters reliably from image data of varying quality. This issue is even more challenging for advanced diffusion methods such as diffusion kurtosis imaging (DKI). Recently, deep learning-based methods have been demonstrated with their potential for robust and efficient computation of diffusion-derived measures. Inspired by these approaches, the current study specifically designed a framework based on a three-dimensional hierarchical convolutional neural network, to jointly reconstruct and harmonize DKI measures from multicenter acquisition to reformulate these to a state-of-the-art hardware using data from traveling subjects. The results from the harmonized data acquired with different protocols show that: 1) the inter-scanner variation of DKI measures within white matter was reduced by 51.5% in mean kurtosis, 65.9% in axial kurtosis, 53.7% in radial kurtosis, and 61.5% in kurtosis fractional anisotropy, respectively; 2) data reliability of each single scanner was enhanced and brought to the level of the reference scanner; and 3) the harmonization network was able to reconstruct reliable DKI values from high data variability. Overall the results demonstrate the feasibility of the proposed deep learning-based method for DKI harmonization and help to simplify the protocol setup procedure for multicenter scanners with different hardware and software configurations.

Longitudinal ComBat: A method for harmonizing longitudinal multi-scanner imaging data

While aggregation of neuroimaging datasets from multiple sites and scanners can yield increased statistical power, it also presents challenges due to systematic scanner effects. This unwanted technical variability can introduce noise and bias into estimation of biological variability of interest. We propose a method for harmonizing longitudinal multi-scanner imaging data based on ComBat, a method originally developed for genomics and later adapted to cross-sectional neuroimaging data. Using longitudinal cortical thickness measurements from 663 participants in the Alzheimer's Disease Neuroimaging Initiative (ADNI) study, we demonstrate the presence of additive and multiplicative scanner effects in various brain regions. We compare estimates of the association between diagnosis and change in cortical thickness over time using three versions of the ADNI data: unharmonized data, data harmonized using cross-sectional ComBat, and data harmonized using longitudinal ComBat. In simulation studies, we show that longitudinal ComBat is more powerful for detecting longitudinal change than cross-sectional ComBat and controls the type I error rate better than unharmonized data with scanner included as a covariate. The proposed method would be useful for other types of longitudinal data requiring harmonization, such as genomic data, or neuroimaging studies of neurodevelopment, psychiatric disorders, or other neurological diseases.

Reproducibility, reliability and variability of FA and MD in the older healthy population: A test-retest multiparametric analysis

•

In older healthy subjects, FA and MD show overall good test-retest reliability & reproducibility.

•

MD is sistematically more reproducible than FA across the entire brain anatomy.

•

FA is more reliable than MD in subcortical white matter regions.

•

In high reliability & low reproducibility regions, variability between subjects is high and statistical power is low.

•

In low reliability & high reproducibility regions, variability between subjects is low and statistical power is high.

Cross-scanner and cross-protocol diffusion MRI data harmonisation: A benchmark database and evaluation of algorithms

Diffusion MRI is being used increasingly in studies of the brain and other parts of the body for its ability to provide quantitative measures that are sensitive to changes in tissue microstructure. However, inter-scanner and inter-protocol differences are known to induce significant measurement variability, which in turn jeopardises the ability to obtain 'truly quantitative measures' and challenges the reliable combination of different datasets. Combining datasets from different scanners and/or acquired at different time points could dramatically increase the statistical power of clinical studies, and facilitate multi-centre research. Even though careful harmonisation of acquisition parameters can reduce variability, inter-protocol differences become almost inevitable with improvements in hardware and sequence design over time, even within a site. In this work, we present a benchmark diffusion MRI database of the same subjects acquired on three distinct scanners with different maximum gradient strength (40, 80, and 300 mT/m), and with 'standard' and 'state-of-the-art' protocols, where the latter have higher spatial and angular resolution. The dataset serves as a useful testbed for method development in cross-scanner/cross-protocol diffusion MRI harmonisation and quality enhancement. Using the database, we compare the performance of five different methods for estimating mappings between the scanners and protocols. The results show that cross-scanner harmonisation of single-shell diffusion data sets can reduce the variability between scanners, and highlight the promises and shortcomings of today's data harmonisation techniques.

Bias Introduced by Multiple Head Coils in MRI Research: An 8 Channel and 32 Channel Coil Comparison

Neuroimaging MRI data in scientific research is increasingly pooled, but the reliability of such studies may be hampered by the use of different hardware elements. This might introduce bias, for example when cross-sectional studies pool data acquired with different head coils, or when longitudinal clinical studies change head coils halfway. In the present study, we aimed to estimate this possible bias introduced by using different head coils to create awareness and to avoid misinterpretation of results. We acquired, with both an 8 channel and 32 channel head coil, T1-weighted, diffusion tensor imaging and resting state fMRI images at 3T MRI (Philips Achieva) with stable acquisition parameters in a large group of cognitively healthy participants (n = 77). Standard analysis methods, i.e., voxel-based morphometry, tract-based spatial statistics and resting state functional network analyses, were used in a within-subject design to compare 8 and 32 channel head coil data. Signal-to-noise ratios (SNR) for both head coils showed similar ranges, although the 32 channel SNR profile was more homogeneous. Our data demonstrates specific patterns of gray and white matter volume differences between head coils (relative volume change of 6 to 9%), related to altered image contrast and therefore, altered tissue segmentation. White matter connectivity (fractional anisotropy and diffusivity measures) showed hemispherical dependent differences between head coils (relative connectivity change of 4 to 6%), and functional connectivity in resting state networks was higher using the 32 channel head coil in posterior cortical areas (relative change up to 27.5%). This study shows that, even when acquisition protocols are harmonized, the results of standardized analysis models can be severely affected by the use of different head coils. Researchers should be aware of this when combining multiple neuroimaging MRI datasets, to prevent coil-related bias and avoid misinterpretation of their findings.

Synthesized b0 for diffusion distortion correction (Synb0-DisCo)

Diffusion magnetic resonance images typically suffer from spatial distortions due to susceptibility induced off-resonance fields, which may affect the geometric fidelity of the reconstructed volume and cause mismatches with anatomical images. State-of-the art susceptibility correction (for example, FSL's TOPUP algorithm) typically requires data acquired twice with reverse phase encoding directions, referred to as blip-up blip-down acquisitions, in order to estimate an undistorted volume. Unfortunately, not all imaging protocols include a blip-up blip-down acquisition, and cannot take advantage of the state-of-the art susceptibility and motion correction capabilities. In this study, we aim to enable TOPUP-like processing with historical and/or limited diffusion imaging data that include only a structural image and single blip diffusion image. We utilize deep learning to synthesize an undistorted non-diffusion weighted image from the structural image, and use the non-distorted synthetic image as an anatomical target for distortion correction. We evaluate the efficacy of this approach (named Synb0-DisCo) and show that our distortion correction process results in better matching of the geometry of undistorted anatomical images, reduces variation in diffusion modeling, and is practically equivalent to having both blip-up and blip-down non-diffusion weighted images.

Vascular burden and APOE ε4 are associated with white matter microstructural decline in cognitively normal older adults

White matter microstructure can be measured with diffusion tensor imaging (DTI). While increasing age is a predictor of white matter (WM) microstructure changes, roles of other possible modifiers, such as cardiovascular risk factors, APOE ε4 allele status and biological sex have not been clarified. We investigated 665 cognitively normal participants from the Baltimore Longitudinal Study of Aging (age 50-95, 56.7% female) with a total of 1384 DTI scans. WM microstructure was assessed by fractional anisotropy (FA) and mean diffusivity (MD). A vascular burden score was defined as the sum of five risk factors (hypertension, obesity, elevated cholesterol, diabetes and smoking status). Linear mixed effects models assessed the association of baseline vascular burden on baseline and on rates of change of FA and MD over a mean follow-up of 3.6 years, while controlling for age, race, and scanner type. We also compared DTI trajectories in APOE ε4 carriers vs. non-carriers and men vs. women. At baseline, higher vascular burden was associated with lower FA and higher MD in many WM structures including association, commissural, and projection fibers. Higher baseline vascular burden was also associated with greater longitudinal decline in FA in the hippocampal part of the cingulum and the fornix (crus)/stria terminalis and splenium of the corpus callosum, and with greater increases in MD in the splenium of the corpus callosum. APOE ε4 carriers did not differ from non-carriers in baseline DTI metrics but had greater decline in FA in the genu and splenium of the corpus callosum. Men had higher FA and lower MD in multiple WM regions at baseline but showed greater increase in MD in the genu of the corpus callosum. Women showed greater decreases over time in FA in the gyrus part of the cingulum, compared to men. Our findings show that modifiable vascular risk factors (1) have a negative impact on white matter microstructure and (2) are associated with faster microstructural deterioration of temporal WM regions and the splenium of the corpus callosum in cognitively normal adults. Reducing vascular burden in aging could modify the rate of WM deterioration and could decrease age-related cognitive decline and impairment.

DTI Profiles for Rapid Description of Cohorts at the Clinical-Research Interface

Normal pressure hydrocephalus (NPH) is a syndrome comprising gait disturbance, cognitive decline and urinary incontinence that is an unique model of reversible brain injury, but it presents as a challenging spectrum of disease cohorts. Diffusion Tensor Imaging (DTI), with its ability to interrogate structural white matter patterns at a microarchitectural level, is a potentially useful tool for the confirmation and characterization of disease cohorts at the clinical-research interface. However, obstacles to its widespread use involve the need for consistent DTI analysis and interpretation tools across collaborator sites. We present the use of DTI profiles, a simplistic methodology to interpret white matter injury patterns based on the morphology of diffusivity parameters. We examined 13 patients with complex NPH, i.e., patients with NPH and overlay from multiple comorbidities, including vascular risk burden and neurodegenerative disease, undergoing extended CSF drainage, clinical assessments, and multi-modal MR imaging. Following appropriate exclusions, we compared the morphology of DTI profiles in such complex NPH patients (n = 12, comprising 4 responders and 8 non-responders) to exemplar DTI profiles from a cohort of classic NPH patients (n = 16) demonstrating responsiveness of white matter injury to ventriculo-peritoneal shunting. In the cohort of complex NPH patients, mean age was 71.3 ± 7.6 years (10 males, 2 females) with a mean MMSE score of 21.1. There were 5 age-matched healthy controls, mean age was 73.4 ± 7.2 years (1 male, 4 females) and mean MMSE score was 26.8. In the exemplar cohort of classic NPH patients, mean age was 74.7 ± 5.9 years (10 males, 6 females) and mean MMSE score was 24.1. There were 9 age-matched healthy controls, mean age was 69.4 ± 9.7 years (4 males, 5 females) and mean MMSE score was 28.6. We found that, despite the challenges of acquiring DTI metrics from differing scanners across collaborator sites and NPH patients presenting as differing cohorts along the spectrum of disease, DTI profiles for responsiveness to interventions were comparable. Distinct DTI characteristics were demonstrated for complex NPH responders vs. non-responders. The morphology of DTI profiles for complex NPH responders mimicked DTI patterns found in predominantly shunt-responsive patients undergoing intervention for classic NPH. However, DTI profiles for complex NPH non-responders was suggestive of atrophy. Our findings suggest that it is possible to use DTI profiles to provide a methodology for rapid description of differing cohorts of disease at the clinical-research interface. By describing DTI measures morphologically, it was possible to consistently compare white matter injury patterns across international collaborator datasets.

Retrospective harmonization of multi-site diffusion MRI data acquired with different acquisition parameters

A joint and integrated analysis of multi-site diffusion MRI (dMRI) datasets can dramatically increase the statistical power of neuroimaging studies and enable comparative studies pertaining to several brain disorders. However, dMRI data sets acquired on multiple scanners cannot be naively pooled for joint analysis due to scanner specific nonlinear effects as well as differences in acquisition parameters. Consequently, for joint analysis, the dMRI data has to be harmonized, which involves removing scanner-specific differences from the raw dMRI signal. In this work, we propose a dMRI harmonization method that is capable of removing scanner-specific effects, while accounting for minor differences in acquisition parameters such as b-value, spatial resolution and number of gradient directions. We validate our algorithm on dMRI data acquired from two sites: Philadelphia Neurodevelopmental Cohort (PNC) with 800 healthy adolescents (ages 8-22 years) and Brigham and Women's Hospital (BWH) with 70 healthy subjects (ages 14-54 years). In particular, we show that gender and age-related maturation differences in different age groups are preserved after harmonization, as measured using effect sizes (small, medium and large), irrespective of the test sample size. Since we use matched control subjects from different scanners to estimate scanner-specific effects, our goal in this work is also to determine the minimum number of well-matched subjects needed from each site to achieve best harmonization results. Our results indicate that at-least 16 to 18 well-matched healthy controls from each site are needed to reliably capture scanner related differences. The proposed method can thus be used for retrospective harmonization of raw dMRI data across sites despite differences in acquisition parameters, while preserving inter-subject anatomical variability.

Scalar diffusion-MRI measures invariant to acquisition parameters: A first step towards imaging biomarkers

An imaging biomarker is a biologic feature in an image that is relevant to a patient's diagnosis or prognosis. In order to qualify as a biomarker, a measure must be robust and reproducible. However, the usual scalar measures derived from diffusion tensor imaging are known to be highly dependent on the variation of the acquisition parameters, which prevents their possible use as biomarkers. In this work, we propose a new set of quantitative measures based on diffusion magnetic resonance imaging from single-shell acquisitions that are designed to be robust to the variations of several acquisition parameters (number of gradient directions, b-value and SNR) while keeping a high discrimination power on differences in the diffusion characteristics of the tissue. These new scalar measures are analytically obtained from a generic diffusion function that does not require the calculation of a diffusion tensor. This way, on one hand, we avoid the use of a specific diffusion model and, on the other hand, we make easier the statistical characterization of the measures. Accordingly, the analysis of the measures bias is carried out and it is used to minimize their dependency with respect to the acquisition noise for different SNRs. The robustness and discrimination power of the measures are tested for different number of gradients, b-values and SNRs using a realistic phantom and three real datasets: (1) 13 control subjects and different acquisition parameters; (2) a public data set from a single subject acquired using multiple shells and (3) 32 schizophrenia patients and 32 age and sex-matched healthy controls with a varying number of gradient directions. The proposed quantitative measures exhibit low variability to the changes of the acquisition parameters, while at the same time they preserve a discrimination power that is able to detect significant changes in the anisotropy of the diffusion.

Multi-site harmonization of diffusion MRI data in a registration framework

Diffusion MRI (dMRI) data acquired on different scanners varies significantly in its content throughout the brain even if the acquisition parameters are nearly identical. Thus, proper harmonization of such data sets is necessary to increase the sample size and thereby the statistical power of neuroimaging studies. In this paper, we present a novel approach to harmonize dMRI data (the raw signal, instead of dMRI derived measures such as fractional anisotropy) using rotation invariant spherical harmonic (RISH) features embedded within a multi-modal image registration framework. All dMRI data sets from all sites are registered to a common template and voxel-wise differences in RISH features between sites at a group level are used to harmonize the signal in a subject-specific manner. We validate our method on diffusion data acquired from seven different sites (two GE, three Philips, and two Siemens scanners) on a group of age-matched healthy subjects. We demonstrate the efficacy of our method by statistically comparing diffusion measures such as fractional anisotropy, mean diffusivity and generalized fractional anisotropy across these sites before and after data harmonization. Validation was also done on a group oftest subjects, which were not used to "learn" the harmonization parameters. We also show results using TBSS before and after harmonization for independent validation of the proposed methodology. Using synthetic data, we show that any abnormality in diffusion measures due to disease is preserved during the harmonization process. Our experimental results demonstrate that, for nearly identical acquisition protocol across sites, scanner-specific differences in the signal can be removed using the proposed method in a model independent manner.

Assessment of quantitative magnetic resonance imaging metrics in the brain through the use of a novel phantom

OBJECTIVE

Multisite and longitudinal neuroimaging studies are important in uncovering trajectories of recovery and neurodegeneration following traumatic brain injury (TBI) and concussion through the use of diffusion tensor imaging (DTI) and other imaging modalities. This study assessed differences in anisotropic diffusion measurement across four scanners using a human and a novel phantom developed in conjunction with the Chronic Effects of Neurotrauma Consortium.

METHOD

Human scans provided measurement within biological tissue, and the novel physical phantom provided measures of anisotropic intra-tubular diffusion to serve as a model for intra-axonal water diffusion. Intra- and inter-scanner measurement variances were compared, and the impact on effect size was calculated.

RESULTS

Intra-scanner test-retest reliability estimates for fractional anisotropy (FA) demonstrated relative stability over testing intervals. The human tissue and phantom showed similar FA ranges, high linearity and large within-device effect sizes. However, inter-scanner measures of FA indicated substantial differences, some of which exceeded typical DTI effect sizes in mild TBI.

CONCLUSION

The diffusion phantom may be used to better elucidate inter-scanner variability in DTI-based measurement and provides an opportunity to better calibrate results obtained from scanners used in multisite and longitudinal studies. Novel solutions are being evaluated to understand and potentially overcome these differences.

White matter lesions: Spatial heterogeneity, links to risk factors, cognition, genetics, and atrophy

OBJECTIVES

To investigate spatial heterogeneity of white matter lesions or hyperintensities (WMH).

METHODS

MRI scans of 1,836 participants (median age 52.2 ± 13.16 years) encompassing a wide age range (22-84 years) from the cross-sectional Study of Health in Pomerania (Germany) were included as discovery set identifying spatially distinct components of WMH using a structural covariance approach. Scans of 307 participants (median age 73.8 ± 10.2 years, with 747 observations) from the Baltimore Longitudinal Study of Aging (United States) were included to examine differences in longitudinal progression of these components. The associations of these components with vascular risk factors, cortical atrophy, Alzheimer disease (AD) genetics, and cognition were then investigated using linear regression.

RESULTS

WMH were found to occur nonuniformly, with higher frequency within spatially heterogeneous patterns encoded by 4 components, which were consistent with common categorizations of deep and periventricular WMH, while further dividing the latter into posterior, frontal, and dorsal components. Temporal trends of the components differed both cross-sectionally and longitudinally. Frontal periventricular WMH were most distinctive as they appeared in the fifth decade of life, whereas the other components appeared later in life during the sixth decade. Furthermore, frontal WMH were associated with systolic blood pressure and with pronounced atrophy including AD-related regions. AD polygenic risk score was associated with the dorsal periventricular component in the elderly. Cognitive decline was associated with the dorsal component.

CONCLUSIONS

These results support the hypothesis that the appearance of WMH follows age and disease-dependent regional distribution patterns, potentially influenced by differential underlying pathophysiologic mechanisms, and possibly with a differential link to vascular and neurodegenerative changes.

Scan-rescan repeatability and cross-scanner comparability of DTI metrics in healthy subjects in the SPRINT-MS multicenter trial

PURPOSE

To assess intrascanner repeatability and cross-scanner comparability for diffusion tensor imaging (DTI) metrics in a multicenter clinical trial.

METHODS

DTI metrics (including longitudinal diffusivity [LD], fractional anisotropy [FA], mean diffusivity [MD], and transverse diffusivity [TD]) from pyramidal tracts for healthy controls were calculated from images acquired on twenty-seven 3T MR scanners (Siemens and GE) with 6 different scanner models and 7 different software versions as part of the NN102/SPRINT-MS clinical trial. Each volunteer underwent two scanning sessions on the same scanner. Signal-to-noise ratio (SNR) and signal-to-noise floor ratio (SNFR) were also assessed.

RESULTS

DTI metrics showed good scan-rescan repeatability. There were no significant differences between scans and rescans in LD, FA, MD, or TD values. Although the cross-scanner coefficient of variation (CV) values for all DTI metrics were <5.7%, significant differences were observed for LD (p < 3.3e-5) and FA (p < 0.0024) when GE scanners were compared with Siemens scanners. Significant differences were also observed for SNR when comparing GE scanners and Siemens Skyra scanners (p < 1.4e-7) and when comparing Siemens Skyra scanners and TIM Trio scanners (p < 1.0e-10). Analysis of background signal also demonstrated differences between GE and Siemens scanners in terms of signal statistics. The measured signal intensity from a background noise region of interest was significantly higher for GE scanners than for Siemens scanners (p < 1.2e-12). Significant differences were also observed for SNFR when comparing GE scanners and Siemens Skyra scanners (p < 2.5e-11), GE scanners and Siemens Trio scanners (p < 7.5e-11), and Siemens Skyra scanners and TIM Trio scanners (p < 2.5e-9).

CONCLUSIONS

The good repeatability of the DTI metrics among the 27 scanners used in this study confirms the feasibility of combining DTI data from multiple centers using high angular resolution sequences. Our observations support the feasibility of longitudinal multicenter clinical trials using DTI outcome measures. The noise floor level and SNFR are important parameters that must be assessed when comparing studies that used different scanner models.

Toward Precision and Reproducibility of Diffusion Tensor Imaging: A Multicenter Diffusion Phantom and Traveling Volunteer Study

BACKGROUND AND PURPOSE

Precision medicine is an approach to disease diagnosis, treatment, and prevention that relies on quantitative biomarkers that minimize the variability of individual patient measurements. The aim of this study was to assess the intersite variability after harmonization of a high-angular-resolution 3T diffusion tensor imaging protocol across 13 scanners at the 11 academic medical centers participating in the Transforming Research and Clinical Knowledge in Traumatic Brain Injury multisite study.

MATERIALS AND METHODS

Diffusion MR imaging was acquired from a novel isotropic diffusion phantom developed at the National Institute of Standards and Technology and from the brain of a traveling volunteer on thirteen 3T MR imaging scanners representing 3 major vendors (GE Healthcare, Philips Healthcare, and Siemens). Means of the DTI parameters and their coefficients of variation across scanners were calculated for each DTI metric and white matter tract.

RESULTS

For the National Institute of Standards and Technology diffusion phantom, the coefficients of variation of the apparent diffusion coefficient across the 13 scanners was <3.8% for a range of diffusivities from 0.4 to 1.1 × 10^-6 mm²/s. For the volunteer, the coefficients of variations across scanners of the 4 primary DTI metrics, each averaged over the entire white matter skeleton, were all <5%. In individual white matter tracts, large central pathways showed good reproducibility with the coefficients of variation consistently below 5%. However, smaller tracts showed more variability, with the coefficients of variation of some DTI metrics reaching 10%.

CONCLUSIONS

The results suggest the feasibility of standardizing DTI across 3T scanners from different MR imaging vendors in a large-scale neuroimaging research study.

Brain atrophy in Alzheimer's Disease and aging

Thanks to its safety and accessibility, magnetic resonance imaging (MRI) is extensively used in clinical routine and research field, largely contributing to our understanding of the pathophysiology of neurodegenerative disorders such as Alzheimer's disease (AD). This review aims to provide a comprehensive overview of the main findings in AD and normal aging over the past twenty years, focusing on the patterns of gray and white matter changes assessed in vivo using MRI. Major progresses in the field concern the segmentation of the hippocampus with novel manual and automatic segmentation approaches, which might soon enable to assess also hippocampal subfields. Advancements in quantification of hippocampal volumetry might pave the way to its broader use as outcome marker in AD clinical trials. Patterns of cortical atrophy have been shown to accurately track disease progression and seem promising in distinguishing among AD subtypes. Disease progression has also been associated with changes in white matter tracts. Recent studies have investigated two areas often overlooked in AD, such as the striatum and basal forebrain, reporting significant atrophy, although the impact of these changes on cognition is still unclear. Future integration of different MRI modalities may further advance the field by providing more powerful biomarkers of disease onset and progression.

Inter-site and inter-scanner diffusion MRI data harmonization

We propose a novel method to harmonize diffusion MRI data acquired from multiple sites and scanners, which is imperative for joint analysis of the data to significantly increase sample size and statistical power of neuroimaging studies. Our method incorporates the following main novelties: i) we take into account the scanner-dependent spatial variability of the diffusion signal in different parts of the brain; ii) our method is independent of compartmental modeling of diffusion (e.g., tensor, and intra/extra cellular compartments) and the acquired signal itself is corrected for scanner related differences; and iii) inter-subject variability as measured by the coefficient of variation is maintained at each site. We represent the signal in a basis of spherical harmonics and compute several rotation invariant spherical harmonic features to estimate a region and tissue specific linear mapping between the signal from different sites (and scanners). We validate our method on diffusion data acquired from seven different sites (including two GE, three Philips, and two Siemens scanners) on a group of age-matched healthy subjects. Since the extracted rotation invariant spherical harmonic features depend on the accuracy of the brain parcellation provided by Freesurfer, we propose a feature based refinement of the original parcellation such that it better characterizes the anatomy and provides robust linear mappings to harmonize the dMRI data. We demonstrate the efficacy of our method by statistically comparing diffusion measures such as fractional anisotropy, mean diffusivity and generalized fractional anisotropy across multiple sites before and after data harmonization. We also show results using tract-based spatial statistics before and after harmonization for independent validation of the proposed methodology. Our experimental results demonstrate that, for nearly identical acquisition protocol across sites, scanner-specific differences can be accurately removed using the proposed method.

Peripheral sphingolipids are associated with variation in white matter microstructure in older adults

Sphingolipids serve important structural and functional roles in cellular membranes and myelin sheaths. Plasma sphingolipids have been shown to predict cognitive decline and Alzheimer's disease. However, the association between plasma sphingolipid levels and brain white matter (WM) microstructure has not been examined. We investigated whether plasma sphingolipids (ceramides and sphingomyelins) were associated with magnetic resonance imaging-based diffusion measures, fractional anisotropy (FA), and mean diffusivity, 10.5 years later in 17 WM regions of 150 cognitively normal adults (mean age 67.2). Elevated ceramide species (C20:0, C22:0, C22:1, and C24:1) were associated with lower FA in multiple WM regions, including total cerebral WM, anterior corona radiata, and the cingulum of the cingulate gyrus. Higher sphingomyelins (C18:1 and C20:1) were associated with lower FA in regions such as the anterior corona radiata and body of the corpus callosum. Furthermore, lower sphingomyelin to ceramide ratios (C22:0, C24:0, and C24:1) were associated with lower FA or higher mean diffusivity in regions including the superior and posterior corona radiata. However, although these associations were significant at the a priori p < 0.05, only associations with some regional diffusion measures for ceramide C22:0 and sphingomyelin C18:1 survived correction for multiple comparisons. These findings suggest plasma sphingolipids are associated with variation in WM microstructure in cognitively normal aging.