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

Axonal injury, sleep disturbances, and memory following traumatic brain injury

OBJECTIVES

Traumatic brain injury (TBI) is associated with sleep deficits, but it is not clear why some report sleep disturbances and others do not. The objective of this study was to assess the associations between axonal injury, sleep, and memory in chronic and acute TBI.

METHODS

Data were acquired from two independent datasets which included 156 older adult veterans (69.8 years) from the Alzheimer's Disease Neuroimaging Initiative (ADNI) with prior moderate-to-severe TBIs and 90 (69.2 years) controls and 374 (39.6 years) from Transforming Research and Clinical Knowledge in TBI (TRACK-TBI) with a recent mild TBI (mTBI) and 87 controls (39.6 years), all who completed an MRI, memory assessment, and sleep questionnaire.

RESULTS

Older adults with a prior TBI had a significant association between axonal injury and sleep disturbances [β = 9.52, 95% CI (4.1, 14.9), p = 0.01]. Axonal injury predicted changes in memory over 1-year in TBI [β = -8.72, 95% CI (-18, -2.7), p = 0.03]. We externally validated those findings in TRACK-TBI where axonal injury within 2 weeks after mTBI was significantly associated with higher sleep disturbances in the TBI group at 2 weeks[β = -7.2, 95% CI (-14, -0.50), p = 0.04], 6 months [β = -16, 95% CI (-24, -7.6), p ≤ 0.01], and 12 months post-injury [β = -11, 95% CI (-19, -0.85), p = 0.03]. These associations were not significant in controls.

INTERPRETATIONS

Axonal injury, specifically to the left anterior internal capsule is robustly associated with sleep disturbances in multiple TBI populations. Early assessment of axonal injury following mTBI could identify those at risk for persistent sleep disturbances following injury.

Reduced cross-scanner variability using vendor-agnostic sequences for single-shell diffusion MRI

PURPOSE

To reduce the inter-scanner variability of diffusion MRI (dMRI) measures between scanners from different vendors by developing a vendor-neutral dMRI pulse sequence using the open-source vendor-agnostic Pulseq platform.

METHODS

We implemented a standard EPI based dMRI sequence in Pulseq. We tested it on two clinical scanners from different vendors (Siemens Prisma and GE Premier), systematically evaluating and comparing the within- and inter-scanner variability across the vendors, using both the vendor-provided and Pulseq dMRI sequences. Assessments covered both a diffusion phantom and three human subjects, using standard error (SE) and Lin's concordance correlation to measure the repeatability and reproducibility of standard DTI metrics including fractional anisotropy (FA) and mean diffusivity (MD).

RESULTS

Identical dMRI sequences were executed on both scanners using Pulseq. On the phantom, the Pulseq sequence showed more than a 2.5× reduction in SE (variability) across Siemens and GE scanners. Furthermore, Pulseq sequences exhibited markedly reduced SE in-vivo, maintaining scan-rescan repeatability while delivering lower variability in FA and MD (more than 50% reduction in cortical/subcortical regions) compared to vendor-provided sequences.

CONCLUSION

The Pulseq diffusion sequence reduces the cross-scanner variability for both phantom and in-vivo data, which will benefit multi-center neuroimaging studies and improve the reproducibility of neuroimaging studies.

Diffusion-weighted MR spectroscopy: Consensus, recommendations, and resources from acquisition to modeling

Brain cell structure and function reflect neurodevelopment, plasticity, and aging; and changes can help flag pathological processes such as neurodegeneration and neuroinflammation. Accurate and quantitative methods to noninvasively disentangle cellular structural features are needed and are a substantial focus of brain research. Diffusion-weighted MRS (dMRS) gives access to diffusion properties of endogenous intracellular brain metabolites that are preferentially located inside specific brain cell populations. Despite its great potential, dMRS remains a challenging technique on all levels: from the data acquisition to the analysis, quantification, modeling, and interpretation of results. These challenges were the motivation behind the organization of the Lorentz Center workshop on "Best Practices & Tools for Diffusion MR Spectroscopy" held in Leiden, the Netherlands, in September 2021. During the workshop, the dMRS community established a set of recommendations to execute robust dMRS studies. This paper provides a description of the steps needed for acquiring, processing, fitting, and modeling dMRS data, and provides links to useful resources.

Emotional Resilience Predicts Preserved White Matter Microstructure Following Mild Traumatic Brain Injury

BACKGROUND

Adult patients with mild traumatic brain injury (mTBI) exhibit distinct phenotypes of emotional and cognitive functioning identified by latent profile analysis of clinical neuropsychological assessments. When discerned early after injury, these latent clinical profiles have been found to improve prediction of long-term outcomes from mTBI. The present study hypothesized that white matter (WM) microstructure is better preserved in an emotionally resilient mTBI phenotype compared with a neuropsychiatrically distressed mTBI phenotype.

METHODS

The present study used diffusion magnetic resonance imaging to investigate and compare WM microstructure in major association, projection, and commissural tracts between the two phenotypes and over time. Diffusion magnetic resonance images from 172 patients with mTBI were analyzed to compute individual diffusion tensor imaging maps at 2 weeks and 6 months after injury.

RESULTS

By comparing the diffusion tensor imaging parameters between the two phenotypes at global, regional, and voxel levels, emotionally resilient patients were shown to have higher axial diffusivity compared with neuropsychiatrically distressed patients early after mTBI. Longitudinal analysis revealed greater compromise of WM microstructure in neuropsychiatrically distressed patients, with greater decrease of global axial diffusivity and more widespread decrease of regional axial diffusivity during the first 6 months after injury compared with emotionally resilient patients.

CONCLUSIONS

These results provide neuroimaging evidence of WM microstructural differences underpinning mTBI phenotypes identified from neuropsychological assessments and show differing longitudinal trajectories of these biological effects. These findings suggest that diffusion magnetic resonance imaging can provide short- and long-term imaging biomarkers of resilience.

Practical considerations of diffusion-weighted MRS with ultra-strong diffusion gradients

Introduction

Diffusion-weighted magnetic resonance spectroscopy (DW-MRS) offers improved cellular specificity to microstructure-compared to water-based methods alone-but spatial resolution and SNR is severely reduced and slow-diffusing metabolites necessitate higher b-values to accurately characterize their diffusion properties. Ultra-strong gradients allow access to higher b-values per-unit time, higher SNR for a given b-value, and shorter diffusion times, but introduce additional challenges such as eddy-current artefacts, gradient non-uniformity, and mechanical vibrations.

Methods

In this work, we present initial DW-MRS data acquired on a 3T Siemens Connectom scanner equipped with ultra-strong (300 mT/m) gradients. We explore the practical issues associated with this manner of acquisition, the steps that may be taken to mitigate their impact on the data, and the potential benefits of ultra-strong gradients for DW-MRS. An in-house DW-PRESS sequence and data processing pipeline were developed to mitigate the impact of these confounds. The interaction of TE, b-value, and maximum gradient amplitude was investigated using simulations and pilot data, whereby maximum gradient amplitude was restricted. Furthermore, two DW-MRS voxels in grey and white matter were acquired using ultra-strong gradients and high b-values.

Results

Simulations suggest T2-based SNR gains that are experimentally confirmed. Ultra-strong gradient acquisitions exhibit similar artefact profiles to those of lower gradient amplitude, suggesting adequate performance of artefact mitigation strategies. Gradient field non-uniformity influenced ADC estimates by up to 4% when left uncorrected. ADC and Kurtosis estimates for tNAA, tCho, and tCr align with previously published literature.

Discussion

In conclusion, we successfully implemented acquisition and data processing strategies for ultra-strong gradient DW-MRS and results indicate that confounding effects of the strong gradient system can be ameliorated, while achieving shorter diffusion times and improved metabolite SNR.

Accurate, repeatable, and geometrically precise diffusion-weighted imaging on a 0.35 T magnetic resonance imaging-guided linear accelerator

Background and purpose

Diffusion weighted imaging (DWI) allows for the interrogation of tissue cellularity, which is a surrogate for cellular proliferation. Previous attempts to incorporate DWI into the workflow of a 0.35 T MR-linac (MRL) have lacked quantitative accuracy. In this study, accuracy, repeatability, and geometric precision of apparent diffusion coefficient (ADC) maps produced using an echo planar imaging (EPI)-based DWI protocol on the MRL system is illustrated, and in vivo potential for longitudinal patient imaging is demonstrated.

Materials and methods

Accuracy and repeatability were assessed by measuring ADC values in a diffusion phantom at three timepoints and comparing to reference ADC values. System-dependent geometric distortion was quantified by measuring the distance between 93 pairs of phantom features on ADC maps acquired on a 0.35 T MRL and a 3.0 T diagnostic scanner and comparing to spatially precise CT images. Additionally, for five sarcoma patients receiving radiotherapy on the MRL, same-day in vivo ADC maps were acquired on both systems, one of which at multiple timepoints.

Results

Phantom ADC quantification was accurate on the 0.35 T MRL with significant discrepancies only seen at high ADC. Average geometric distortions were 0.35 (±0.02) mm and 0.85 (±0.02) mm in the central slice and 0.66 (±0.04) mm and 2.14 (±0.07) mm at 5.4 cm off-center for the MRL and diagnostic system, respectively. In the sarcoma patients, a mean pretreatment ADC of 910x10-6 (±100x10-6) mm2/s was measured on the MRL.

Conclusions

The acquisition of accurate, repeatable, and geometrically precise ADC maps is possible at 0.35 T with an EPI approach.

PREdiction and Diagnosis using Imaging and Clinical biomarkers Trial in Traumatic Brain Injury (PREDICT-TBI) study protocol: an observational, prospective, multicentre cohort study for the prediction of outcome in moderate-to-severe TBI

INTRODUCTION

Traumatic brain injury (TBI) is a heterogeneous condition with a broad spectrum of injury severity, pathophysiological processes and variable outcomes. For moderate-to-severe TBI survivors, recovery is often protracted and outcomes can range from total dependence to full recovery. Despite advances in medical treatment options, prognosis remains largely unchanged. The objective of this study is to develop a machine learning predictive model for neurological outcomes at 6 months in patients with a moderate-to-severe TBI, incorporating longitudinal clinical, multimodal neuroimaging and blood biomarker predictor variables.

METHODS AND ANALYSIS

A prospective, observational, cohort study will enrol 300 patients with moderate-to-severe TBI from seven Australian hospitals over 3 years. Candidate predictors including demographic and general health variables, and longitudinal clinical, neuroimaging (CT and MRI), blood biomarker and patient-reported outcome measures will be collected at multiple time points within the acute phase of injury. The predictor variables will populate novel machine learning models to predict the Glasgow Outcome Scale Extended 6 months after injury. The study will also expand on current prognostic models by including novel blood biomarkers (circulating cell-free DNA), and the results of quantitative neuroimaging such as Quantitative Susceptibility Mapping and Dynamic Contrast Enhanced MRI as predictor variables.

ETHICS AND DISSEMINATION

Ethical approval has been obtained by the Royal Brisbane and Women's Hospital Human Research Ethics Committee, Queensland. Participants or their substitute decision-maker/s will receive oral and written information about the study before providing written informed consent. Study findings will be disseminated by peer-review publications and presented at national and international conferences and clinical networks.

TRIAL REGISTRATION NUMBER

ACTRN12620001360909.

Machine Learning for Brain MRI Data Harmonisation: A Systematic Review

BACKGROUND

Magnetic Resonance Imaging (MRI) data collected from multiple centres can be heterogeneous due to factors such as the scanner used and the site location. To reduce this heterogeneity, the data needs to be harmonised. In recent years, machine learning (ML) has been used to solve different types of problems related to MRI data, showing great promise.

OBJECTIVE

This study explores how well various ML algorithms perform in harmonising MRI data, both implicitly and explicitly, by summarising the findings in relevant peer-reviewed articles. Furthermore, it provides guidelines for the use of current methods and identifies potential future research directions.

METHOD

This review covers articles published through PubMed, Web of Science, and IEEE databases through June 2022. Data from studies were analysed based on the criteria of Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA). Quality assessment questions were derived to assess the quality of the included publications.

RESULTS

a total of 41 articles published between 2015 and 2022 were identified and analysed. In the review, MRI data has been found to be harmonised either in an implicit (n = 21) or an explicit (n = 20) way. Three MRI modalities were identified: structural MRI (n = 28), diffusion MRI (n = 7) and functional MRI (n = 6).

CONCLUSION

Various ML techniques have been employed to harmonise different types of MRI data. There is currently a lack of consistent evaluation methods and metrics used across studies, and it is recommended that the issue be addressed in future studies. Harmonisation of MRI data using ML shows promises in improving performance for ML downstream tasks, while caution should be exercised when using ML-harmonised data for direct interpretation.

Magnetic resonance biomarker assessment software (MR-BIAS): an automated open-source tool for the ISMRM/NIST system phantom

Objective.To provide an open-source software for repeatable and efficient quantification ofT₁andT₂relaxation times with the ISMRM/NIST system phantom. Quantitative magnetic resonance imaging (qMRI) biomarkers have the potential to improve disease detection, staging and monitoring of treatment response. Reference objects, such as the system phantom, play a major role in translating qMRI methods into the clinic. The currently available open-source software for ISMRM/NIST system phantom analysis, Phantom Viewer (PV), includes manual steps that are subject to variability.Approach.We developed the Magnetic Resonance BIomarker Assessment Software (MR-BIAS) to automatically extract system phantom relaxation times. The inter-observer variability (IOV) and time efficiency of MR-BIAS and PV was observed in six volunteers analysing three phantom datasets. The IOV was measured with the coefficient of variation (CV) of percent bias (%bias) inT₁andT₂with respect to NMR reference values. The accuracy of MR-BIAS was compared to a custom script from a published study of twelve phantom datasets. This included comparison of overall bias and %bias for variable inversion recovery (T₁^VIR), variable flip angle (T₁^VFA) and multiple spin-echo (T₂^MSE) relaxation models.Main results.MR-BIAS had a lower mean CV withT₁^VIR(0.03%) andT₂^MSE(0.05%) in comparison to PV withT₁^VIR(1.28%) andT₂^MSE(4.55%). The mean analysis duration was 9.7 times faster for MR-BIAS (0.8 min) than PV (7.6 min). There was no statistically significant difference in the overall bias, or the %bias for the majority of ROIs, as calculated by MR-BIAS or the custom script for all models.Significance.MR-BIAS has demonstrated repeatable and efficient analysis of the ISMRM/NIST system phantom, with comparable accuracy to previous studies. The software is freely available to the MRI community, providing a framework to automate required analysis tasks, with the flexibility to explore open questions and accelerate biomarker research.

Diffusion-Weighted Imaging in Mild Traumatic Brain Injury: A Systematic Review of the Literature

There is evidence that diffusion-weighted imaging (DWI) is able to detect tissue alterations following mild traumatic brain injury (mTBI) that may not be observed on conventional neuroimaging; however, findings are often inconsistent between studies. This systematic review assesses patterns of differences in DWI metrics between those with and without a history of mTBI. A PubMed literature search was performed using relevant indexing terms for articles published prior to May 14, 2020. Findings were limited to human studies using DWI in mTBI. Articles were excluded if they were not full-length, did not contain original data, if they were case studies, pertained to military populations, had inadequate injury severity classification, or did not report post-injury interval. Findings were reported independently for four subgroups: acute/subacute pediatric mTBI, acute/subacute adult mTBI, chronic adult mTBI, and sport-related concussion, and all DWI acquisition and analysis methods used were included. Patterns of findings between studies were reported, along with strengths and weaknesses of the current state of the literature. Although heterogeneity of sample characteristics and study methods limited the consistency of findings, alterations in DWI metrics were most commonly reported in the corpus callosum, corona radiata, internal capsule, and long association pathways. Many acute/subacute pediatric studies reported higher FA and lower ADC or MD in various regions. In contrast, acute/subacute adult studies most commonly indicate lower FA within the context of higher MD and RD. In the chronic phase of recovery, FA may remain low, possibly indicating overall demyelination or Wallerian degeneration over time. Longitudinal studies, though limited, generally indicate at least a partial normalization of DWI metrics over time, which is often associated with functional improvement. We conclude that DWI is able to detect structural mTBI-related abnormalities that may persist over time, although future DWI research will benefit from larger samples, improved data analysis methods, standardized reporting, and increasing transparency.

Evaluation of Apparent Diffusion Coefficient Repeatability and Reproducibility for Preclinical MRIs Using Standardized Procedures and a Diffusion-Weighted Imaging Phantom

Relevant to co-clinical trials, the goal of this work was to assess repeatability, reproducibility, and bias of the apparent diffusion coefficient (ADC) for preclinical MRIs using standardized procedures for comparison to performance of clinical MRIs. A temperature-controlled phantom provided an absolute reference standard to measure spatial uniformity of these performance metrics. Seven institutions participated in the study, wherein diffusion-weighted imaging (DWI) data were acquired over multiple days on 10 preclinical scanners, from 3 vendors, at 6 field strengths. Centralized versus site-based analysis was compared to illustrate incremental variance due to processing workflow. At magnet isocenter, short-term (intra-exam) and long-term (multiday) repeatability were excellent at within-system coefficient of variance, wCV [±CI] = 0.73% [0.54%, 1.12%] and 1.26% [0.94%, 1.89%], respectively. The cross-system reproducibility coefficient, RDC [±CI] = 0.188 [0.129, 0.343] µm2/ms, corresponded to 17% [12%, 31%] relative to the reference standard. Absolute bias at isocenter was low (within 4%) for 8 of 10 systems, whereas two high-bias (>10%) scanners were primary contributors to the relatively high RDC. Significant additional variance (>2%) due to site-specific analysis was observed for 2 of 10 systems. Base-level technical bias, repeatability, reproducibility, and spatial uniformity patterns were consistent with human MRIs (scaled for bore size). Well-calibrated preclinical MRI systems are capable of highly repeatable and reproducible ADC measurements.

Validation of cross-sectional and longitudinal ComBat harmonization methods for magnetic resonance imaging data on a travelling subject cohort

Background

The growth in multi-center neuroimaging studies generated a need for methods that mitigate the differences in hardware and acquisition protocols across sites i.e., scanner effects. ComBat harmonization methods have shown promise but have not yet been tested on all the data types commonly studied with magnetic resonance imaging (MRI). This study aimed to validate neuroCombat, longCombat and gamCombat on both structural and diffusion metrics in both cross-sectional and longitudinal data.

Methods

We used a travelling subject design whereby 73 healthy volunteers contributed 161 scans across two sites and four machines using one T1 and five diffusion MRI protocols. Scanner was defined as a composite of site, machine and protocol. A common pipeline extracted two structural metrics (volumes and cortical thickness) and two diffusion tensor imaging metrics (mean diffusivity and fractional anisotropy) for seven regions of interest including gray and (except for cortical thickness) white matter regions.

Results

Structural data exhibited no significant scanner effect and therefore did not benefit from harmonization in our particular cohort. Indeed, attempting harmonization obscured the true biological effect for some regions of interest. Diffusion data contained marked scanner effects and was successfully harmonized by all methods, resulting in smaller scanner effects and better detection of true biological effects. LongCombat less effectively reduced the scanner effect for cross-sectional white matter data but had a slightly lower probability of incorrectly finding group differences in simulations, compared to neuroCombat and gamCombat. False positive rates for all methods and all metrics did not significantly exceed 5%.

Conclusions

Statistical harmonization of structural data is not always necessary and harmonization in the absence of a scanner effect may be harmful. Harmonization of diffusion MRI data is highly recommended with neuroCombat, longCombat and gamCombat performing well in cross-sectional and longitudinal settings.

Optimization of quasi-diffusion magnetic resonance imaging for quantitative accuracy and time-efficient acquisition

PURPOSE

Quasi-diffusion MRI (QDI) is a novel quantitative technique based on the continuous time random walk model of diffusion dynamics. QDI provides estimates of the diffusion coefficient, <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mspace/> <mml:msub><mml:mi>D</mml:mi> <mml:mrow><mml:mn>1</mml:mn> <mml:mo>,</mml:mo> <mml:mn>2</mml:mn></mml:mrow> </mml:msub> </mml:mrow> <mml:annotation>$$ {D}_{1,2} $$</mml:annotation></mml:semantics> </mml:math> in mm²  s^-1 and a fractional exponent, <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mi>α</mml:mi></mml:mrow> <mml:annotation>$$ \upalpha $$</mml:annotation></mml:semantics> </mml:math> , defining the non-Gaussianity of the diffusion signal decay. Here, the b-value selection for rapid clinical acquisition of QDI tensor imaging (QDTI) data is optimized.

METHODS

Clinically appropriate QDTI acquisitions were optimized in healthy volunteers with respect to a multi-b-value reference (MbR) dataset comprising 29 diffusion-sensitized images arrayed between <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mi>b</mml:mi> <mml:mo>=</mml:mo> <mml:mn>0</mml:mn></mml:mrow> <mml:annotation>$$ b=0 $$</mml:annotation></mml:semantics> </mml:math> and 5000 s mm^-2 . The effects of varying maximum b-value ( <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow><mml:msub><mml:mi>b</mml:mi> <mml:mi>max</mml:mi></mml:msub> </mml:mrow> <mml:annotation>$$ {b}_{\mathrm{max}} $$</mml:annotation></mml:semantics> </mml:math> ), number of b-value shells, and the effects of Rician noise were investigated.

RESULTS

QDTI measures showed <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow><mml:msub><mml:mi>b</mml:mi> <mml:mi>max</mml:mi></mml:msub> </mml:mrow> <mml:annotation>$$ {b}_{\mathrm{max}} $$</mml:annotation></mml:semantics> </mml:math> dependence, most significantly for <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mi>α</mml:mi></mml:mrow> <mml:annotation>$$ \upalpha $$</mml:annotation></mml:semantics> </mml:math> in white matter, which monotonically decreased with higher <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow><mml:msub><mml:mi>b</mml:mi> <mml:mi>max</mml:mi></mml:msub> </mml:mrow> <mml:annotation>$$ {b}_{\mathrm{max}} $$</mml:annotation></mml:semantics> </mml:math> leading to improved tissue contrast. Optimized 2 b-value shell acquisitions showed small systematic differences in QDTI measures relative to MbR values, with overestimation of <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mspace/> <mml:mspace/> <mml:msub><mml:mi>D</mml:mi> <mml:mrow><mml:mn>1</mml:mn> <mml:mo>,</mml:mo> <mml:mn>2</mml:mn></mml:mrow> </mml:msub> </mml:mrow> <mml:annotation>$$ \kern0.50em {D}_{1,2} $$</mml:annotation></mml:semantics> </mml:math> and underestimation of <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mi>α</mml:mi></mml:mrow> <mml:annotation>$$ \upalpha $$</mml:annotation></mml:semantics> </mml:math> in white matter, and overestimation of <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow><mml:msub><mml:mi>D</mml:mi> <mml:mrow><mml:mn>1</mml:mn> <mml:mo>,</mml:mo> <mml:mn>2</mml:mn></mml:mrow> </mml:msub> </mml:mrow> <mml:annotation>$$ {D}_{1,2} $$</mml:annotation></mml:semantics> </mml:math> and <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mi>α</mml:mi></mml:mrow> <mml:annotation>$$ \upalpha $$</mml:annotation></mml:semantics> </mml:math> anisotropies in gray and white matter. Additional shells improved the accuracy, precision, and reliability of QDTI estimates with 3 and 4 shells at <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow><mml:msub><mml:mi>b</mml:mi> <mml:mi>max</mml:mi></mml:msub> <mml:mo>=</mml:mo> <mml:mn>5000</mml:mn></mml:mrow> <mml:annotation>$$ {b}_{\mathrm{max}}=5000 $$</mml:annotation></mml:semantics> </mml:math>  s mm^-2 , and 4 b-value shells at <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow><mml:msub><mml:mi>b</mml:mi> <mml:mi>max</mml:mi></mml:msub> <mml:mo>=</mml:mo> <mml:mn>3960</mml:mn></mml:mrow> <mml:annotation>$$ {b}_{\mathrm{max}}=3960 $$</mml:annotation></mml:semantics> </mml:math>  s mm^-2 , providing minimal bias in <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow><mml:msub><mml:mi>D</mml:mi> <mml:mrow><mml:mn>1</mml:mn> <mml:mo>,</mml:mo> <mml:mn>2</mml:mn></mml:mrow> </mml:msub> </mml:mrow> <mml:annotation>$$ {D}_{1,2} $$</mml:annotation></mml:semantics> </mml:math> and <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mi>α</mml:mi></mml:mrow> <mml:annotation>$$ \upalpha $$</mml:annotation></mml:semantics> </mml:math> compared to the MbR.

CONCLUSION

A highly detailed optimization of non-Gaussian dMRI for in vivo brain imaging was performed. QDI provided robust parameterization of non-Gaussian diffusion signal decay in clinically feasible imaging times with high reliability, accuracy, and precision of QDTI measures.

Diffusion Tensor Imaging Reveals Elevated Diffusivity of White Matter Microstructure that Is Independently Associated with Long-Term Outcome after Mild Traumatic Brain Injury: A TRACK-TBI Study

Diffusion tensor imaging (DTI) literature on single-center studies contains conflicting results regarding acute effects of mild traumatic brain injury (mTBI) on white matter (WM) microstructure and the prognostic significance. This larger-scale multi-center DTI study aimed to determine how acute mTBI affects WM microstructure over time and how early WM changes affect long-term outcome. From Transforming Research and Clinical Knowledge in Traumatic Brain Injury (TRACK-TBI), a cohort study at 11 United States level 1 trauma centers, a total of 391 patients with acute mTBI ages 17 to 60 years were included and studied at two weeks and six months post-injury. Demographically matched friends or family of the participants were the control group (n = 148). Axial diffusivity (AD), fractional anisotropy (FA), mean diffusivity (MD), and radial diffusivity (RD) were the measures of WM microstructure. The primary outcome was the Glasgow Outcome Scale Extended (GOSE) score of injury-related functional limitations across broad life domains at six months post-injury. The AD, MD, and RD were higher and FA was lower in mTBI versus friend control (FC) at both two weeks and six months post-injury throughout most major WM tracts of the cerebral hemispheres. In the mTBI group, AD and, to a lesser extent, MD decreased in WM from two weeks to six months post-injury. At two weeks post-injury, global WM AD and MD were both independently associated with six-month incomplete recovery (GOSE <8 vs = 8) even after accounting for demographic, clinical, and other imaging factors. DTI provides reliable imaging biomarkers of dynamic WM microstructural changes after mTBI that have utility for patient selection and treatment response in clinical trials. Continued technological advances in the sensitivity, specificity, and precision of diffusion magnetic resonance imaging hold promise for routine clinical application in mTBI.

Validating the accuracy of multispectral metal artifact suppressed diffusion-weighted imaging

BACKGROUND

Diffusion-weighted imaging (DWI) provides quantitative measurement of random water displacement in tissue as calculated by the apparent diffusion coefficient (ADC). While heavily utilized in stroke and oncology applications, DWI is a promising tool to map microstructural changes in musculoskeletal applications including evaluation of synovial reactions resulting from total hip arthroplasty (THA). One major challenge facing the application of DWI in THA is the significant artifacts related to the conventional echo-planar imaging (EPI) readout used. Multispectral imaging (MSI) techniques, including the multiacquisition with variable resonance image combination (MAVRIC), have been shown to effectively reduce metallic susceptibility artifacts around total joint replacements to render clinically useful images. Recently, a 2D periodically rotated overlapping parallel line with enhanced reconstruction (PROPELLER) FSE acquisition that incorporates a diffusion preparation pulse with 2D-MAVRIC has been developed to mitigate both distortion and dropout artifacts. While there have been some preliminary assessments of DWI-MAVRIC, the repeatability of DWI-MAVRIC and the effects of key parameters, such as the number of spectral bins, are unknown.

PURPOSE

To evaluate the quantitative accuracy of DWI-MAVRIC as compared to conventional diffusion sequences.

METHODS

A diffusion phantom with different reference diffusivities (ADC = 113-1123 μm2 /s) was used. Scans were performed on two 1.5T MRI scanners. DWI-EPI and DWI-MAVRIC were acquired in both the axial and coronal planes. Three spatial offsets (0 cm, 10 cm left, and 10 cm right off iso-center) were used to evaluate effects of off-isocenter positioning. To assess intraday and interday repeatability, DWI-EPI and DWI-MAVRIC acquisitions were repeated on one scanner at same-day and 9-month intervals. To assess inter-scanner repeatability, DWI-EPI and DWI-MAVRIC acquisitions were compared between two scanners. ADC maps were generated with and without gradient nonlinearity correction (GNC). Linear regression, correlation, and error statistics were determined between calculated and reference ADC values. Bland-Altman plots were generated to evaluate intraday, interday, and interscanner repeatability.

RESULTS

DWI-MAVRIC had excellent correlation to reference values but at reduced linearity (r = 1.00, slope = 0.91-0.94) as compared to DWI-EPI (r = 1.00, slope = 0.99-1.01). A greater than 5% ADC bias was observed at the lowest ADC values, predominantly in the DWI-MAVRIC scans. ADC values did not vary with DWI-MAVRIC parameters. DWI-EPI acquisitions had intraday, interday, and interscanner repeatability of 3.18 μm2 /s, 19.2 μm2 /s, and 20.2 μm2 /s, respectively. DWI-MAVRIC acquisitions had inferior intraday, interday, and interscanner repeatability of 13.3 μm2 /s, 44.7 μm2 /s, 110 μm2 /s, respectively. Lower ADC errors were found at isocenter, as compared to the left and right positions. GNC reduced the absolute error by 0.31% ± 0.89%, 3.6% ± 1.4%, 0.65% ± 2.4% for the center, left, and right positions, respectively.

CONCLUSIONS

DWI-MAVRIC provides good linearity with respect to reference values and good intra- and interday repeatability.

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.

Clinical translation of quantitative magnetic resonance imaging biomarkers - An overview and gap analysis of current practice

PURPOSE

This overview of the current landscape of quantitative magnetic resonance imaging biomarkers (qMR IBs) aims to support the standardisation of academic IBs to assist their translation to clinical practice.

METHODS

We used three complementary approaches to investigate qMR IB use and quality management practices within the UK: 1) a literature search of qMR and quality management terms during 2011-2015 and 2016-2020; 2) a database search for clinical research studies using qMR IBs during 2016-2020; and 3) a survey to ascertain the current availability and quality management practices for clinical MRI scanners and associated equipment at research institutions across the UK.

RESULTS

The analysis showed increased use of all qMR methods between the periods 2011-2015 and 2016-2020 and diffusion-tensor MRI and volumetry to be popular methods. However, the "translation ratio" of journal articles to clinical research studies was higher for qMR methods that have evidence of clinical translation via a commercial route, such as fat fraction and T2 mapping. The number of journal articles citing quality management terms doubled between the periods 2011-2015 and 2016-2020; although, its proportion relative to all journal articles only increased by 3.0%. The survey suggested that quality assurance (QA) and quality control (QC) of data acquisition procedures are under-reported in the literature and that QA/QC of acquired data/data analysis are under-developed and lack consistency between institutions.

CONCLUSIONS

We summarise current attempts to standardise and translate qMR IBs, and conclude by outlining the ideal quality management practices and providing a gap analysis between current practice and a metrological standard.

Reproducibility and repeatability of magnetic resonance imaging in dementia

PURPOSE

Individualised predictive models of cognitive decline require disease-monitoring markers that are repeatable. For wide-spread adoption, such markers also need to be reproducible at different locations. This study assessed the repeatability and reproducibility of MRI markers derived from a dementia protocol.

METHODS

Six participants were scanned at three different sites with a 3T MRI scanner. The protocol employed: T1-weighted (T1w) imaging, resting state functional MRI (rsfMRI), arterial spin labelling (ASL), diffusion-weighted imaging (DWI), T2-weighted fluid attenuation inversion recovery (FLAIR), T2-weighted (T2w) imaging, and susceptibility weighted imaging (SWI). Participants were scanned repeatedly, up to six times over a maximum period of five years. One participant was also scanned a further three times on sequential days on one scanner. Fifteen derived metrics were computed from the seven different modalities.

RESULTS

Reproducibility (coefficient of variation; CoV, across sites) was best for T1w derived grey matter, white matter and hippocampal volume (CoV < 1.5%), compared to rsfMRI and SWI derived metrics (CoV, 19% and 21%). For a given metric, long-term repeatability (CoV across time) was comparable to reproducibility, with short-term repeatability considerably better.

CONCLUSIONS

Reproducibility and repeatability were assessed for a suite of markers calculated from a dementia MRI protocol. In general, structural markers were less variable than functional MRI markers. Variability over time on the same scanner was comparable to variability measured across different scanners. Overall, the results support the viability of multi-site longitudinal studies for monitoring cognitive decline.

WA, Boyle J, Coll‐Font J, Dall'Armellina E, Ennis DB, Ferreira PF, Kalra P, Kolipaka A, Kozerke S, Lohr D, Mongeon F, Moulin K, Nguyen C, Nielles‐Vallespin S, Raterman B, Schreiber LM, Scott AD, Sosnovik DE, Stoeck CT, Tous C, Tunnicliffe EM, Weng AM, Croisille P, Viallon M, Schneider JE. Validation of cardiac diffusion tensor imaging sequences: A multicentre test-retest phantom study

Cardiac diffusion tensor imaging (DTI) is an emerging technique for the in vivo characterisation of myocardial microstructure, and there is a growing need for its validation and standardisation. We sought to establish the accuracy, precision, repeatability and reproducibility of state-of-the-art pulse sequences for cardiac DTI among 10 centres internationally. Phantoms comprising 0%-20% polyvinylpyrrolidone (PVP) were scanned with DTI using a product pulsed gradient spin echo (PGSE; N = 10 sites) sequence, and a custom motion-compensated spin echo (SE; N = 5) or stimulated echo acquisition mode (STEAM; N = 5) sequence suitable for cardiac DTI in vivo. A second identical scan was performed 1-9 days later, and the data were analysed centrally. The average mean diffusivities (MDs) in 0% PVP were (1.124, 1.130, 1.113) x 10^-3  mm² /s for PGSE, SE and STEAM, respectively, and accurate to within 1.5% of reference data from the literature. The coefficients of variation in MDs across sites were 2.6%, 3.1% and 2.1% for PGSE, SE and STEAM, respectively, and were similar to previous studies using only PGSE. Reproducibility in MD was excellent, with mean differences in PGSE, SE and STEAM of (0.3 ± 2.3, 0.24 ± 0.95, 0.52 ± 0.58) x 10^-5  mm² /s (mean ± 1.96 SD). We show that custom sequences for cardiac DTI provide accurate, precise, repeatable and reproducible measurements. Further work in anisotropic and/or deforming phantoms is warranted.

MAGNIMS recommendations for harmonization of MRI data in MS multicenter studies

•

Sharing data from cooperative studies is essential to develop new biomarkers in MS.

•

Differences in MRI acquisition, analysis, storage represent a substantial constraint.

•

We review the state of the art and developments in the harmonization of MRI.

•

We provide recommendations to harmonize large MRI datasets in the MS field.

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.

The Case for Optimized Edge-Centric Tractography at Scale

The anatomic validity of structural connectomes remains a significant uncertainty in neuroimaging. Edge-centric tractography reconstructs streamlines in bundles between each pair of cortical or subcortical regions. Although edge bundles provides a stronger anatomic embedding than traditional connectomes, calculating them for each region-pair requires exponentially greater computation. We observe that major speedup can be achieved by reducing the number of streamlines used by probabilistic tractography algorithms. To ensure this does not degrade connectome quality, we calculate the identifiability of edge-centric connectomes between test and re-test sessions as a proxy for information content. We find that running PROBTRACKX2 with as few as 1 streamline per voxel per region-pair has no significant impact on identifiability. Variation in identifiability caused by streamline count is overshadowed by variation due to subject demographics. This finding even holds true in an entirely different tractography algorithm using MRTrix. Incidentally, we observe that Jaccard similarity is more effective than Pearson correlation in calculating identifiability for our subject population.

Quantitative MRI Harmonization to Maximize Clinical Impact: The RIN-Neuroimaging Network

Neuroimaging studies often lack reproducibility, one of the cardinal features of the scientific method. Multisite collaboration initiatives increase sample size and limit methodological flexibility, therefore providing the foundation for increased statistical power and generalizable results. However, multisite collaborative initiatives are inherently limited by hardware, software, and pulse and sequence design heterogeneities of both clinical and preclinical MRI scanners and the lack of benchmark for acquisition protocols, data analysis, and data sharing. We present the overarching vision that yielded to the constitution of RIN-Neuroimaging Network, a national consortium dedicated to identifying disease and subject-specific in-vivo neuroimaging biomarkers of diverse neurological and neuropsychiatric conditions. This ambitious goal needs efforts toward increasing the diagnostic and prognostic power of advanced MRI data. To this aim, 23 Italian Scientific Institutes of Hospitalization and Care (IRCCS), with technological and clinical specialization in the neurological and neuroimaging field, have gathered together. Each IRCCS is equipped with high- or ultra-high field MRI scanners (i.e., ≥3T) for clinical or preclinical research or has established expertise in MRI data analysis and infrastructure. The actions of this Network were defined across several work packages (WP). A clinical work package (WP1) defined the guidelines for a minimum standard clinical qualitative MRI assessment for the main neurological diseases. Two neuroimaging technical work packages (WP2 and WP3, for clinical and preclinical scanners) established Standard Operative Procedures for quality controls on phantoms as well as advanced harmonized quantitative MRI protocols for studying the brain of healthy human participants and wild type mice. Under FAIR principles, a web-based e-infrastructure to store and share data across sites was also implemented (WP4). Finally, the RIN translated all these efforts into a large-scale multimodal data collection in patients and animal models with dementia (i.e., case study). The RIN-Neuroimaging Network can maximize the impact of public investments in research and clinical practice acquiring data across institutes and pathologies with high-quality and highly-consistent acquisition protocols, optimizing the analysis pipeline and data sharing procedures.

Conformance of a 3T Radiotherapy MRI Scanner to the QIBA Diffusion Profile

Purpose

To assess the technical performance of the apparent diffusion coefficient (ADC) on a dedicated 3T radiotherapy scanner, using a standardized phantom and sequences. Investigations into factors that could impact the technical performance of ADC in the clinic were also completed, including changing the slice‐encoded imaging direction and the reference sample ADC value.

Methods

ADC acquisitions were performed monthly on an isotropic diffusion phantom over 1 year. Measurements of ADC %bias, coefficients of variation for short‐/long‐term repeatability and precision (CV_ST/CV_LT and CV_P), and b‐value dependency (Dep_b) were calculated. The measurements were then assessed according to the Quantitative Imaging Biomarker Alliance (QIBA) Diffusion Profile specifications.

Results

The average of all measurements over the year was within Profile recommended ranges. This included when testing was performed in different imaging directions, and on samples that had different ADC reference values (0.4–1.1 μm²/ms). Results in the axial plane for the central water vial included a bias of +0.05%, CV_ST /CV_LT/CV_P = 0.1%/ 0.9%/0.4% and Dep_b = 0.4%.

Conclusions

The technical performance of ADC on a radiotherapy dedicated MRI scanner over the course of 12 months was considered conformant to the QIBA Profile. Quantifying these metrics and factors that may affect the performance is essential in progressing the use of ADC clinically: ensuring that the observed change of ADC in a tissue is due to a physiological response and not measurement variability.

Longitudinal assessment of quality assurance measurements in a 1.5 T MR-linac: Part II-Magnetic resonance imaging

PURPOSE

To describe and report longitudinal quality assurance (QA) measurements for the magnetic resonance imaging (MRI) component of the Elekta Unity MR-linac during the first year of clinical use in our institution.

MATERIALS AND METHODS

The performance of the MRI component of Unity was evaluated with daily, weekly, monthly, and annual QA testing. The measurements monitor image uniformity, signal-to-noise ratio (SNR), resolution/detectability, slice position/thickness, linearity, central frequency, and geometric accuracy. In anticipation of routine use of quantitative imaging (qMRI), we characterize B0/B1 uniformity and the bias/reproducibility of longitudinal/transverse relaxation times (T1/T2) and apparent diffusion coefficient (ADC). Tolerance levels for QA measurements of qMRI biomarkers are derived from weekly monitoring of T1, T2, and ADC.

RESULTS

The 1-year assessment of QA measurements shows that daily variations in each MR quality metric are well below the threshold for failure. Routine testing procedures can reproducibly identify machine issues. The longitudinal three-dimensional (3D) geometric analysis reveals that the maximum distortion in a diameter of spherical volume (DSV) of 20, 30, 40, and 50 cm is 0.4, 0.6, 1.0, and 3.1 mm, respectively. The main source of distortion is gradient nonlinearity. Maximum peak-to-peak B0 inhomogeneity is 3.05 ppm, with gantry induced B0 inhomogeneities an order of magnitude smaller. The average deviation from the nominal B1 is within 2%, with minimal dependence on gantry angle. Mean ADC, T1, and T2 values are measured with high reproducibility. The median coefficient of variation for ADC, T1, and T2 is 1.3%, 1.1%, and 0.5%, respectively. The median bias for ADC, T1, and T2 is -0.8%, -0.1%, and 3.9%, respectively.

CONCLUSION

The MRI component of Unity operates within the guidelines and recommendations for scanner performance and stability. Our findings support the recently published guidance in establishing clinically acceptable tolerance levels for image quality. Highly reproducible qMRI measurements are feasible in Unity.

Technical Note: Temperature and concentration dependence of water diffusion in polyvinylpyrrolidone solutions

Objective

The goal of this work is to provide temperature and concentration calibration of water diffusivity in polyvinylpyrrolidone (PVP) solutions used in phantoms to assess system bias and linearity in apparent diffusion coefficient (ADC) measurements.

Method

ADC measurements were performed for 40 kDa (K40) PVP of six concentrations (0%, 10%, 20%, 30%, 40%, and 50% by weight) at three temperatures (19.5°C, 22.5°C, and 26.4°C), with internal phantom temperature monitored by optical thermometer (±0.2°C). To achieve ADC measurement and fit accuracy of better than 0.5%, three orthogonal diffusion gradients were calibrated using known water diffusivity at 0°C and system gradient nonlinearity maps. Noise‐floor fit bias was also controlled by limiting the maximum b‐value used for ADC calculation of each sample. The ADC temperature dependence was modeled by Arrhenius functions of each PVP concentration. The concentration dependence was modeled by quadratic function for ADC normalized by the theoretical water diffusion values. Calibration coefficients were obtained from linear regression model fits.

Results

Measured phantom ADC values increased with temperature and decreasing PVP concentration, [PVP]. The derived Arrhenius model parameters for [PVP] between 0% and 50%, are reported and can be used for K40 ADC temperature calibration with absolute ADC error within ±0.016 μm²/ms. Arrhenius model fit parameters normalized to water value scaled with [PVP] between 10% and 40%, and proportional change in activation energy increased faster than collision frequency. ADC normalization by water diffusivity, D_W, from the Speedy–Angell relation accounted for the bulk of temperature dependence (±0.035 μm²/ms) and yielded quadratic calibration for ADC_PVP/D_W = (12.5 ± 0.7) ·10⁻⁵·[PVP]² − (23.2 ± 0.3)·10⁻³·[PVP]+1, nearly independent of PVP molecular weight and temperature.

Conclusion

The study provides ground‐truth ADC values for K40 PVP solutions commonly used in diffusion phantoms for scanning at ambient room temperature. The described procedures and the reported calibration can be used for quality control and standardization of measured ADC values of PVP at different concentrations and temperatures.

Examining brain white matter after pediatric mild traumatic brain injury using neurite orientation dispersion and density imaging: An A-CAP study

•

We examined white matter microstructure after pediatric mTBI using NODDI and DTI.

•

Children with mTBI did not significantly differ from those with OI on any metrics.

•

Minor alterations, if any, may be present in children at the post-acute stage after mTBI.

•

Large longitudinal studies are needed to understand long-term brain changes post injury.

Background

Pediatric mild traumatic brain injury (mTBI) affects millions of children annually. Diffusion tensor imaging (DTI) is sensitive to axonal injuries and white matter microstructure and has been used to characterize the brain changes associated with mild traumatic brain injury (mTBI). Neurite orientation dispersion and density imaging (NODDI) is a diffusion model that can provide additional insight beyond traditional DTI metrics, but has not been examined in pediatric mTBI. The goal of this study was to employ DTI and NODDI to gain added insight into white matter alterations in children with mTBI compared to children with mild orthopedic injury (OI).

Methods

Children (mTBI n = 320, OI n = 176) aged 8–16.99 years (12.39 ± 2.32 years) were recruited from emergency departments at five hospitals across Canada and underwent 3 T MRI on average 11 days post-injury. DTI and NODDI metrics were calculated for seven major white matter tracts and compared between groups using univariate analysis of covariance controlling for age, sex, and scanner type. False discovery rate (FDR) was used to correct for multiple comparisons.

Results

Univariate analysis revealed no significant group main effects or interactions in DTI or NODDI metrics. Fractional anisotropy and neurite density index in all tracts exhibited a significant positive association with age and mean diffusivity in all tracts exhibited a significant negative association with age in the whole sample.

Conclusions

Overall, there were no significant differences between mTBI and OI groups in brain white matter microstructure from either DTI or NODDI in the seven tracts. This indicates that mTBI is associated with relatively minor white matter differences, if any, at the post-acute stage. Brain differences may evolve at later stages of injury, so longitudinal studies with long-term follow-up are needed.

MASiVar: Multisite, multiscanner, and multisubject acquisitions for studying variability in diffusion weighted MRI

PURPOSE

Diffusion-weighted imaging allows investigators to identify structural, microstructural, and connectivity-based differences between subjects, but variability due to session and scanner biases is a challenge.

METHODS

To investigate DWI variability, we present MASiVar, a multisite data set consisting of 319 diffusion scans acquired at 3 T from b = 1000 to 3000 s/mm2 across 14 healthy adults, 83 healthy children (5 to 8 years), three sites, and four scanners as a publicly available, preprocessed, and de-identified data set. With the adult data, we demonstrate the capacity of MASiVar to simultaneously quantify the intrasession, intersession, interscanner, and intersubject variability of four common DWI processing approaches: (1) a tensor signal representation, (2) a multi-compartment neurite orientation dispersion and density model, (3) white-matter bundle segmentation, and (4) structural connectomics. Respectively, we evaluate region-wise fractional anisotropy, mean diffusivity, and principal eigenvector; region-wise CSF volume fraction, intracellular volume fraction, and orientation dispersion index; bundle-wise shape, volume, fractional anisotropy, and length; and whole connectome correlation and maximized modularity, global efficiency, and characteristic path length.

RESULTS

We plot the variability in these measures at each level and find that it consistently increases with intrasession to intersession to interscanner to intersubject effects across all processing approaches and that sometimes interscanner variability can approach intersubject variability.

CONCLUSIONS

This study demonstrates the potential of MASiVar to more globally investigate DWI variability across multiple levels and processing approaches simultaneously and suggests harmonization between scanners for multisite analyses should be considered before inference of group differences on subjects.

Diffusion MRI harmonization enables joint-analysis of multicentre data of patients with cerebral small vessel disease

OBJECTIVES

Acquisition-related differences in diffusion magnetic resonance imaging (dMRI) hamper pooling of multicentre data to achieve large sample sizes. A promising solution is to harmonize the raw diffusion signal using rotation invariant spherical harmonic (RISH) features, but this has not been tested in elderly subjects. Here we aimed to establish if RISH harmonization effectively removes acquisition-related differences in multicentre dMRI of elderly subjects with cerebral small vessel disease (SVD), while preserving sensitivity to disease effects.

METHODS

Five cohorts of patients with SVD (N = 397) and elderly controls (N = 175) with 3 Tesla MRI on different systems were included. First, to establish effectiveness of harmonization, the RISH method was trained with data of 13 to 15 age and sex-matched controls from each site. Fractional anisotropy (FA) and mean diffusivity (MD) were compared in matched controls between sites using tract-based spatial statistics (TBSS) and voxel-wise analysis, before and after harmonization. Second, to assess sensitivity to disease effects, we examined whether the contrast (effect sizes of FA, MD and peak width of skeletonized MD - PSMD) between patients and controls within each site remained unaffected by harmonization. Finally, we evaluated the association between white matter hyperintensity (WMH) burden, FA, MD and PSMD using linear regression analyses both within individual cohorts as well as with pooled scans from multiple sites, before and after harmonization.

RESULTS

Before harmonization, significant differences in FA and MD were observed between matched controls of different sites (p < 0.05). After harmonization these site-differences were removed. Within each site, RISH harmonization did not alter the effect sizes of FA, MD and PSMD between patients and controls (relative change in Cohen's d = 4 %) nor the strength of association with WMH volume (relative change in R2 = 2.8 %). After harmonization, patient data of all sites could be aggregated in a single analysis to infer the association between WMH volume and FA (R2 = 0.62), MD (R2 = 0.64), and PSMD (R2 = 0.60).

CONCLUSIONS

We showed that RISH harmonization effectively removes acquisition-related differences in dMRI of elderly subjects while preserving sensitivity to SVD-related effects. This study provides proof of concept for future multicentre SVD studies with pooled datasets.

Same Brain, Different Look?-The Impact of Scanner, Sequence and Preprocessing on Diffusion Imaging Outcome Parameters

In clinical diagnostics and longitudinal studies, the reproducibility of MRI assessments is of high importance in order to detect pathological changes, but developments in MRI hard- and software often outrun extended periods of data acquisition and analysis. This could potentially introduce artefactual changes or mask pathological alterations. However, if and how changes of MRI hardware, scanning protocols or preprocessing software affect complex neuroimaging outcomes from, e.g., diffusion weighted imaging (DWI) remains largely understudied. We therefore compared DWI outcomes and artefact severity of 121 healthy participants (age range 19–54 years) who underwent two matched DWI protocols (Siemens product and Center for Magnetic Resonance Research sequence) at two sites (Siemens 3T Magnetom Verio and Skyra^fit). After different preprocessing steps, fractional anisotropy (FA) and mean diffusivity (MD) maps, obtained by tensor fitting, were processed with tract-based spatial statistics (TBSS). Inter-scanner and inter-sequence variability of skeletonised FA values reached up to 5% and differed largely in magnitude and direction across the brain. Skeletonised MD values differed up to 14% between scanners. We here demonstrate that DTI outcome measures strongly depend on imaging site and software, and that these biases vary between brain regions. These regionally inhomogeneous biases may exceed and considerably confound physiological effects such as ageing, highlighting the need to harmonise data acquisition and analysis. Future studies thus need to implement novel strategies to augment neuroimaging data reliability and replicability.

Axonal marker neurofilament light predicts long-term outcomes and progressive neurodegeneration after traumatic brain injury

[Figure: see text].

Biomarkers for Traumatic Brain Injury: Data Standards and Statistical Considerations

Recent biomarker innovations hold potential for transforming diagnosis, prognostic modeling, and precision therapeutic targeting of traumatic brain injury (TBI). However, many biomarkers, including brain imaging, genomics, and proteomics, involve vast quantities of high-throughput and high-content data. Management, curation, analysis, and evidence synthesis of these data are not trivial tasks. In this review, we discuss data management concepts and statistical and data sharing strategies when dealing with biomarker data in the context of TBI research. We propose that application of biomarkers involves three distinct steps-discovery, evaluation, and evidence synthesis. First, complex/big data has to be reduced to useful data elements at the stage of biomarker discovery. Second, inferential statistical approaches must be applied to these biomarker data elements for assessment of biomarker clinical utility and validity. Last, synthesis of relevant research is required to support practice guidelines and enable health decisions informed by the highest quality, up-to-date evidence available. We focus our discussion around recent experiences from the International Traumatic Brain Injury Research (InTBIR) initiative, with a specific focus on four major clinical projects (Transforming Research and Clinical Knowledge in TBI, Collaborative European NeuroTrauma Effectiveness Research in TBI, Collaborative Research on Acute Traumatic Brain Injury in Intensive Care Medicine in Europe, and Approaches and Decisions in Acute Pediatric TBI Trial), which are currently enrolling subjects in North America and Europe. We discuss common data elements, data collection efforts, data-sharing opportunities, and challenges, as well as examine the statistical techniques required to realize successful adoption and use of biomarkers in the clinic as a foundation for precision medicine in TBI.

Empirical validation of gradient field models for an accurate ADC measured on clinical 3T MR systems in body oncologic applications

PURPOSE

To empirically corroborate vendor-provided gradient nonlinearity (GNL) characteristics and demonstrate efficient GNL bias correction for human brain apparent diffusion coefficient (ADC) across 3T MR systems and spatial locations.

METHODS

Spatial distortion vector fields (DVF) were mapped in 3D using a surface fiducial array phantom for individual gradient channels on three 3T MR platforms from different vendors. Measured DVF were converted into empirical 3D GNL tensors and compared with their theoretical counterparts derived from vendor-provided spherical harmonic (SPH) coefficients. To illustrate spatial impact of GNL on ADC, diffusion weighted imaging using three orthogonal gradient directions was performed on a volunteer brain positioned at isocenter (as a reference) and offset superiorly by 10-17 cm (>10% predicted GNL bias). The SPH tensor-based GNL correction was applied to individual DWI gradient directions, and derived ADC was compared with low-bias reference for human brain white matter (WM) ROIs.

RESULTS

Empiric and predicted GNL errors were comparable for all three studied 3T MR systems, with <1.0% differences in the median and width of spatial histograms for individual GNL tensor elements. Median (±width) of ADC (10-3mm2/s) histograms measured at isocenter in WM reference ROIs from three MR systems were: 0.73 ± 0.11, 0.71 ± 0.14, 0.74 ± 0.17, and at off-isocenters (before versus after GNL correction) were respectively 0.63 ± 0.14 versus 0.72 ± 0.11, 0.53 ± 0.16 versus 0.74 ± 0.18, and 0.65 ± 0.16 versus 0.76 ± 0.18.

CONCLUSION

The phantom-based spatial distortion measurements validated vendor-provided gradient fields, and accurate WM ADC was recovered regardless of spatial locations and clinical MR platforms using system-specific tensor-based GNL correction for routine DWI.

In-depth characterization of a mouse model of post-traumatic epilepsy for biomarker and drug discovery

Post-traumatic epilepsy (PTE) accounts for 5% of all epilepsies and 10–20% of the acquired forms. The latency between traumatic brain injury (TBI) and epilepsy onset in high-risk patients offers a therapeutic window for intervention to prevent or improve the disease course. However, progress towards effective treatments has been hampered by the lack of sensitive prognostic biomarkers of PTE, and of therapeutic targets. There is therefore a pressing clinical need for preclinical PTE models suitable for biomarker discovery and drug testing. We characterized in-depth a model of severe TBI induced by controlled cortical impact evolving into PTE in CD1 adult male mice. To identify sensitive measures predictive of PTE development and severity, TBI mice were longitudinally monitored by video-electrocorticography (ECoG), examined by MRI, and tested for sensorimotor and cognitive deficits and locomotor activity. At the end of the video-ECoG recording mice were killed for brain histological analysis. PTE occurred in 58% of mice with frequent motor seizures (one seizure every other day), as determined up to 5 months post-TBI. The weight loss of PTE mice in 1 week after TBI correlated with the number of spontaneous seizures at 5 months. Moreover, the recovery rate of the sensorimotor deficit detected by the SNAP test before the predicted time of epilepsy onset was significantly lower in PTE mice than in those without epilepsy. Neuroscore, beam walk and cognitive deficit were similar in all TBI mice. The increase in the contusion volume, the volume of forebrain regions contralateral to the lesioned hemisphere and white matter changes over time assessed by MRI were similar in PTE and no-PTE mice. However, brain histology showed a more pronounced neuronal cell loss in the cortex and hippocampus contralateral to the injured hemisphere in PTE than in no-PTE mice. The extensive functional and neuropathological characterization of this TBI model, provides a tool to identify sensitive measures of epilepsy development and severity clinically useful for increasing PTE prediction in high-risk TBI patients. The high PTE incidence and spontaneous seizures frequency in mice provide an ideal model for biomarker discovery and for testing new drugs.

MRI diffusion phantoms: ADC and relaxometric measurement comparisons between polyacrylamide and agarose gels

PURPOSE

The aim of this study is to compare polyacrylamide and agarose gels, as components of a simple MRI phantom, for the measurements of Apparent Diffusion Coefficient (ADC), T1 and T2 relaxation times.

MATERIALS AND METHODS

Five (5) test tubes with polyacrylamide gels of different monomer concentrations and six (6) test tubes of different agarose gel concentrations were used as a phantom for ADC, T1 and T2 measurements, which were expressed as 2D color parametric maps, on a 1.5 T clinical MRI system. ADC and T2 maps were calculated utilizing a Weighted Linear (WL) regression fitting algorithm. T1 maps were calculated utilizing a standard non-linear fitting algorithm.

RESULTS

In agarose gels, ADC measurements are independent of the agarose concentration, whereas the T1 and T2 relaxation times decrease with increasing agarose concentration. On the contrary, in polyacrylamide gels, ADC measurements decrease quadratically while increasing the monomer concentration, whereas the T1 and T2 relaxation times reveal a linear decrease with increasing monomer concentration.

CONCLUSION

Polyacrylamide gels can serve as a better means for simulating ADC values, as compared with the agarose gels used in this study.

Molecular and Diffusion Tensor Imaging Biomarkers of Traumatic Brain Injury: Principles for Investigation and Integration

The last 20 years have seen the advent of new technologies that enhance the diagnosis and prognosis of traumatic brain injury (TBI). There is recognition that TBI affects the brain beyond initial injury, in some cases inciting a progressive neuropathology that leads to chronic impairments. Medical researchers are now searching for biomarkers to detect and monitor this condition. Perhaps the most promising developments are in the biomolecular and neuroimaging domains. Molecular assays can identify proteins indicative of neuronal injury and/or degeneration. Diffusion imaging now allows sensitive evaluations of the brain's cellular microstructure. As the pace of discovery accelerates, it is important to survey the research landscape and identify promising avenues of investigation. In this review, we discuss the potential of molecular and diffusion tensor imaging (DTI) biomarkers in TBI research. Integration of these technologies could advance models of disease prognosis, ultimately improving care. To date, however, few studies have explored relationships between molecular and DTI variables in patients with TBI. Here, we provide a short primer on each technology, review the latest research, and discuss how these biomarkers may be incorporated in future studies.

Quantitative MRI: Defining repeatability, reproducibility and accuracy for prostate cancer imaging biomarker development

INTRODUCTION

Quantitative MRI (qMRI) parameters have been increasingly used to develop predictive models to accurately monitor treatment response in prostate cancer after radiotherapy. To reliably detect changes in signal due to treatment response, predictive models require qMRI parameters with high repeatability and reproducibility. The purpose of this study was to measure qMRI parameter uncertainties in both commercial and in-house developed phantoms to guide the development of robust predictive models for monitoring treatment response.

MATERIALS AND METHODS

ADC, T1, and R2* values were acquired across three 3 T scanners with a prostate-specific qMRI protocol using the NIST/ISMRM system phantom, RSNA/NIST diffusion phantom, and an in-house phantom. A B1 field map was acquired to correct for flip angle inhomogeneity in T1 maps. All sequences were repeated in each scan to assess within-session repeatability. Weekly scans were acquired on one scanner for three months with the in-house phantom. Between-session repeatability was measured with test-retest scans 6-months apart on all scanners with all phantoms. Accuracy, defined as percentage deviation from reference value for ADC and T1, was evaluated using the system and diffusion phantoms. Repeatability and reproducibility coefficients of variation (%CV) were calculated for all qMRI parameters on all phantoms.

RESULTS

Overall, repeatability CV of ADC was <2.40%, reproducibility CV was <3.98%, and accuracy ranged between -8.0% to 2.7% across all scanners. Applying B1 correction on T1 measurements significantly improved the repeatability and reproducibility (p<0.05) but increased error in accuracy (p<0.001). Repeatability and reproducibility of R2* was <4.5% and <7.3% respectively in the system phantom across all scanners.

CONCLUSION

Repeatability, reproducibility, and accuracy in qMRI parameters from a prostate-specific protocol was estimated using both commercial and in-house phantoms. Results from this work will be used to identify robust qMRI parameters for use in the development of predictive models to longitudinally monitor treatment response for prostate cancer in current and future clinical trials.

Ficoll as testing material for diffusion weighted imaging-quality assurance phantoms

PURPOSE

The aim of this work is to test the use of aqueous solutions of Ficoll®**, a highly branched polymer displaying crowding properties, to build a phantom suitable for Diffusion Weighted Imaging (DWI) in Magnetic Resonance Imaging (MRI).

METHODS

We developed a test object made of a cylindrical plastic container with a precise geometrical arrangement suitable for measuring several samples at the same time. The container was designed to host single vials with variable geometry and number, and to fit inside common commercial head coils for MRI scanners. In our experiments, vials were filled with 8 aqueous solutions of Ficoll 70 and Ficoll 400 spanning a range of polymer concentration from 5 to 30% by weight. Vials containing ultra-pure water were also used as reference. Experiments were performed on both 1.5 and 3 T clinical scanners (GE, Philips and Siemens), under the conditions of a standard clinical examination.

RESULTS

The geometry of the phantom provided reduced imaging artifacts, especially image distortions at magnetic interfaces. We found that the Apparent Diffusion Coefficient (ADC) varied in the range of 0.00125-0.00223 mm²/s and decreased with Ficoll concentration. ADC vs Ficoll concentration exhibited a linear trend. Results were consistent over time and among different MRI clinical scanners, showing an average variability of 3% at 1.5 T and of 7.5% at 3 T. Moreover, no substantial difference was found between Ficoll 70 and 400. By varying Ficoll concentration, ADC can be modulated to approach tissue-mimicking values. Preliminary results for relaxation measurements proved that both T₁ and T₂ decreased with Ficoll concentration in the ranges 1.3-2.4 s and 150-800 ms respectively.

CONCLUSIONS

In this work, we propose a 3D phantom design based on the widespread crowding agent Ficoll, which is suitable for DWI quality assurance purposes in MRI acquisitions. Aqueous Ficoll solutions provide good performance in terms of stability, ease of preparation, and safety.

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 Structural MRI in Neurologically Healthy Adults

BACKGROUND

Structural brain MRI measures are frequently examined in both healthy and clinical groups, so an understanding of how these measures vary over time is desirable.

PURPOSE

To test the stability of structural brain MRI measures over time.

POPULATION

In all, 112 healthy volunteers across four sites.

STUDY TYPE

Retrospective analysis of prospectively acquired data.

FIELD STRENGTH/SEQUENCE

3 T, magnetization prepared - rapid gradient echo, and single-shell diffusion sequence.

ASSESSMENT

Diffusion, cortical thickness, and volume data from the sensorimotor network were assessed for stability over time across 3 years. Two sites used a Siemens MRI scanner, two sites a Philips scanner.

STATISTICAL TESTS

The stability of structural measures across timepoints was assessed using intraclass correlation coefficients (ICC) for absolute agreement, cutoff ≥0.80, indicating high reliability. Mixed-factorial analysis of variance (ANOVA) was used to examine between-site and between-scanner type differences in individuals over time.

RESULTS

All cortical thickness and gray matter volume measures in the sensorimotor network, plus all diffusivity measures (fractional anisotropy plus mean, axial and radial diffusivities) for primary and premotor cortices, primary somatosensory thalamic connections, and the cortico-spinal tract met ICC. The majority of measures differed significantly between scanners, with a trend for sites using Siemens scanners to produce larger values for connectivity, cortical thickness, and volume measures than sites using Philips scanners.

DATA CONCLUSION

Levels of reliability over time for all tested structural MRI measures were generally high, indicating that any differences between measurements over time likely reflect underlying biological differences rather than inherent methodological variability.

LEVEL OF EVIDENCE

4.

TECHNICAL EFFICACY STAGE

1.

Harmonization of Brain Diffusion MRI: Concepts and Methods

MRI diffusion data suffers from significant inter- and intra-site variability, which hinders multi-site and/or longitudinal diffusion studies. This variability may arise from a range of factors, such as hardware, reconstruction algorithms and acquisition settings. To allow a reliable comparison and joint analysis of diffusion data across sites and over time, there is a clear need for robust data harmonization methods. This review article provides a comprehensive overview of diffusion data harmonization concepts and methods, and their limitations. Overall, the methods for the harmonization of multi-site diffusion images can be categorized in two main groups: diffusion parametric map harmonization (DPMH) and diffusion weighted image harmonization (DWIH). Whereas DPMH harmonizes the diffusion parametric maps (e.g., FA, MD, and MK), DWIH harmonizes the diffusion-weighted images. Defining a gold standard harmonization technique for dMRI data is still an ongoing challenge. Nevertheless, in this paper we provide two classification tools, namely a feature table and a flowchart, which aim to guide the readers in selecting an appropriate harmonization method for their study.

Image acquisition and quality assurance in the Boston Adolescent Neuroimaging of Depression and Anxiety study

The Connectomes Related to Human Diseases (CRHD) initiative was developed with the Human Connectome Project (HCP) to provide high-resolution, open-access, multi-modal MRI data to better understand the neural correlates of human disease. Here, we present an introduction to a CRHD project, the Boston Adolescent Neuroimaging of Depression and Anxiety (BANDA) study, which is collecting multimodal neuroimaging, clinical, and neuropsychological data from 225 adolescents (ages 14-17), 150 of whom are expected to have a diagnosis of depression and/or anxiety. Our transdiagnostic recruitment approach samples the full spectrum of depressed/anxious symptoms and their comorbidity, consistent with NIMH Research Domain Criteria (RDoC). We focused on an age range that is critical for brain development and for the onset of mental illness. This project sought to harmonize imaging sequences, hardware, and functional tasks with other HCP studies, although some changes were made to canonical HCP methods to accommodate our study population and questions. We present a thorough overview of our imaging sequences, hardware, and scanning protocol. We detail similarities and differences between this study and other HCP studies. We evaluate structural-, diffusion-, and functional-image-quality measures that may be influenced by clinical factors (e.g., disorder, symptomatology). Signal-to-noise and motion estimates from the first 140 adolescents suggest minimal influence of clinical factors on image quality. We anticipate enrollment of an additional 85 participants, most of whom are expected to have a diagnosis of anxiety and/or depression. Clinical and neuropsychological data from the first 140 participants are currently freely available through the National Institute of Mental Health Data Archive (NDA).

Test-retest reliability and sample size estimates after MRI scanner relocation

OBJECTIVE

Many factors can contribute to the reliability and robustness of MRI-derived metrics. In this study, we assessed the reliability and reproducibility of three MRI modalities after an MRI scanner was relocated to a new hospital facility.

METHODS

Twenty healthy volunteers (12 females, mean age (standard deviation) ​= ​41 (11) years, age range [25-66]) completed three MRI sessions. The first session (S1) was one week prior to the 3T GE HDxt scanner relocation. The second (S2) occurred nine weeks after S1 and at the new location; a third session (S3) was acquired 4 weeks after S2. At each session, we acquired structural T1-weighted, pseudo-continuous arterial spin labelled, and diffusion tensor imaging sequences. We used longitudinal processing streams to create 12 summary MRI metrics, including total gray matter (GM), cortical GM, subcortical GM, white matter (WM), and lateral ventricle volume; mean cortical thickness; total surface area; average gray matter perfusion, and average diffusion tensor metrics along principal white matter pathways. We compared mean MRI values and variance at the old scanner location to multiple sessions at the new location using Bayesian multi-level regression models. K-fold cross validation allowed identification of important predictors. Whole-brain analyses were used to investigate any regional differences. Furthermore, we calculated within-subject coefficient of variation (wsCV), intraclass correlation coefficient (ICC), and dice similarity index (SI) of cortical segmentations across scanner relocation and within-site. Additionally, we estimated sample sizes required to robustly detect a 4% difference between two groups across MRI metrics.

RESULTS

All global MRI metrics exhibited little mean difference and small variability (bar cortical gray matter perfusion) both across scanner relocation and within-site repeat. T1- and DTI-derived tissue metrics showed ​< ​|0.3|% mean difference and <1.2% variance across scanner location and <|0.4|% mean difference and <0.8% variance within the new location, with between-site intraclass correlation coefficient (ICC) ​> ​0.80 and within-subject coefficient of variation (wsCV) ​< ​1.4%. Mean cortical gray matter perfusion had the highest between-session variability (6.7% [0.3, 16.7], estimate [95% uncertainty interval]), and hence the smallest ICC (0.71 [0.44,0.92]) and largest wsCV (13.4% [5.4, 18.1]). No global metric exhibited evidence of a meaningful mean difference between scanner locations. However, surface area showed evidence of a mean difference within-site repeat (between S2 and S3). Whole-brain analyses revealed no significant areas of difference between scanner relocation or within-site. For all metrics, we found no support for a systematic difference in variance across relocation sites compared to within-site test-retest reliability. Necessary sample sizes to detect a 4% difference between two independent groups varied from a maximum of n ​= ​362 per group (cortical gray matter perfusion), to total gray matter volume (n ​= ​114), average fractional anisotropy (n ​= ​23), total gray matter volume normalized by intracranial volume (n ​= ​19), and axial diffusivity (n ​= ​3 per group).

CONCLUSION

Cortical gray matter perfusion was the most variable metric investigated (necessitating large sample sizes to identify group differences), with other metrics showing substantially less variability. Scanner relocation appeared to have a negligible effect on variability of the global MRI metrics tested. This manuscript reports within-site test-retest variability to act as a tool for calculating sample size in future investigations. Our results suggest that when all other parameters are held constant (e.g., sequence parameters and MRI processing), the effect of scanner relocation is indistinguishable from within-site variability, but may need to be considered depending on the question being investigated.

Repeatability of Quantitative Diffusion-Weighted Imaging Metrics in Phantoms, Head-and-Neck and Thyroid Cancers: Preliminary Findings

The aim of this study was to establish the repeatability measures of quantitative Gaussian and non-Gaussian diffusion metrics using diffusion-weighted imaging (DWI) data from phantoms and patients with head-and-neck and papillary thyroid cancers. The Quantitative Imaging Biomarker Alliance (QIBA) DWI phantom and a novel isotropic diffusion kurtosis imaging phantom were scanned at 3 different sites, on 1.5T and 3T magnetic resonance imaging systems, using standardized multiple b-value DWI acquisition protocol. In the clinical component of this study, a total of 60 multiple b-value DWI data sets were analyzed for test–retest, obtained from 14 patients (9 head-and-neck squamous cell carcinoma and 5 papillary thyroid cancers). Repeatability of quantitative DWI measurements was assessed by within-subject coefficient of variation (wCV%) and Bland–Altman analysis. In isotropic diffusion kurtosis imaging phantom vial with 2% ceteryl alcohol and behentrimonium chloride solution, the mean apparent diffusion (D_app × 10⁻³ mm²/s) and kurtosis (K_app, unitless) coefficient values were 1.02 and 1.68 respectively, capturing in vivo tumor cellularity and tissue microstructure. For the same vial, D_app and K_app mean wCVs (%) were ≤1.41% and ≤0.43% for 1.5T and 3T across 3 sites. For pretreatment head-and-neck squamous cell carcinoma, apparent diffusion coefficient, D, D*, K, and f mean wCVs (%) were 2.38%, 3.55%, 3.88%, 8.0%, and 9.92%, respectively; wCVs exhibited a higher trend for papillary thyroid cancers. Knowledge of technical precision and bias of quantitative imaging metrics enables investigators to properly design and power clinical trials and better discern between measurement variability versus biological change.

Scan-rescan and inter-vendor reproducibility of neurite orientation dispersion and density imaging metrics

Purpose

The reproducibility of neurite orientation dispersion and density imaging (NODDI) metrics in the human brain has not been explored across different magnetic resonance (MR) scanners from different vendors. This study aimed to evaluate the scan–rescan and inter-vendor reproducibility of NODDI metrics in white and gray matter of healthy subjects using two 3-T MR scanners from two vendors.

Methods

Ten healthy subjects (7 males; mean age 30 ± 7 years, range 23–37 years) were included in the study. Whole-brain diffusion-weighted imaging was performed with b-values of 1000 and 2000 s/mm² using two 3-T MR scanners from two different vendors. Automatic extraction of the region of interest was performed to obtain NODDI metrics for whole and localized areas of white and gray matter. The coefficient of variation (CoV) and intraclass correlation coefficient (ICC) were calculated to assess the scan–rescan and inter-vendor reproducibilities of NODDI metrics.

Results

The scan–rescan and inter-vendor reproducibility of NODDI metrics (intracellular volume fraction and orientation dispersion index) were comparable with those of diffusion tensor imaging (DTI) metrics. However, the inter-vendor reproducibilities of NODDI (CoV = 2.3–14%) were lower than the scan–rescan reproducibility (CoV: scanner A = 0.8–3.8%; scanner B = 0.8–2.6%). Compared with the finding of DTI metrics, the reproducibility of NODDI metrics was lower in white matter and higher in gray matter.

Conclusion

The lower inter-vendor reproducibility of NODDI in some brain regions indicates that data acquired from different MRI scanners should be carefully interpreted.

Replication and generalization in applied neuroimaging

There is much interest in translating neuroimaging findings into meaningful clinical diagnostics. The goal of scientific discoveries differs from clinical diagnostics. Scientific discoveries must replicate under a specific set of conditions; to translate to the clinic we must show that findings using purpose-built scientific instruments will be observable in clinical populations and instruments. Here we describe and evaluate data and computational methods designed to translate a scientific observation to a clinical setting. Using diffusion weighted imaging (DWI), Wahl et al. (2010) observed that across subjects the mean fractional anisotropy (FA) of homologous pairs of tracts is highly correlated. We hypothesize that this is a fundamental biological trait that should be present in most healthy participants, and deviations from this assessment may be a useful diagnostic metric. Using this metric as an illustration of our methods, we analyzed six pairs of homologous white matter tracts in nine different DWI datasets with 44 subjects each. Considering the original FA measurement as a baseline, we show that the new metric is between 2 and 4 times more precise when used in a clinical context. Our framework to translate research findings into clinical practice can be applied, in principle, to other neuroimaging results.

Associations between modifiable risk factors and white matter of the aging brain: insights from diffusion tensor imaging studies

There is increasing interest in factors that may modulate white matter (WM) breakdown and, consequentially, age-related cognitive and behavioral deficits. Recent diffusion tensor imaging studies have examined the relationship of such factors with WM microstructure. This review summarizes the evidence regarding the relationship between WM microstructure and recognized modifiable factors, including hearing loss, hypertension, diabetes, obesity, smoking, depressive symptoms, physical (in) activity, and social isolation, as well as sleep disturbances, diet, cognitive training, and meditation. Current cross-sectional evidence suggests a clear link between loss of WM integrity (lower fractional anisotropy and higher mean diffusivity) and hypertension, obesity, diabetes, and smoking; a relationship that seems to hold for hearing loss, social isolation, depressive symptoms, and sleep disturbances. Physical activity, cognitive training, diet, and meditation, on the other hand, may protect WM with aging. Preliminary evidence from cross-sectional studies of treated risk factors suggests that modification of factors could slow down negative effects on WM microstructure. Careful intervention studies are needed for this literature to contribute to public health initiatives going forward.

Reproducibility of Structural and Diffusion Tensor Imaging in the TACERN Multi-Center Study

BACKGROUND

Multi-site MRI studies are often necessary for recruiting sufficiently sized samples when studying rare conditions. However, they require pooling data from multiple scanners into a single data set, and therefore it is critical to evaluate the variability of quantitative MRI measures within and across scanners used in multi-site studies. The aim of this study was to evaluate the reproducibility of structural and diffusion weighted (DW) MRI measurements acquired on seven scanners at five medical centers as part of the Tuberous Sclerosis Complex Autism Center of Excellence Research Network (TACERN) multisite study.

METHODS

The American College of Radiology (ACR) phantom was imaged monthly to measure reproducibility of signal intensity and uniformity within and across seven 3T scanners from General Electric, Philips, and Siemens vendors. One healthy adult male volunteer was imaged repeatedly on all seven scanners under the TACERN structural and DW protocol (5 b = 0 s/mm2 and 30 b = 1000 s/mm2) over a period of 5 years (age 22-27 years). Reproducibility of inter- and intra-scanner brain segmentation volumes and diffusion tensor imaging metrics fractional anisotropy (FA) and mean diffusivity (MD) within white matter regions was quantified with coefficient of variation.

RESULTS

The American College of Radiology Phantom signal intensity and uniformity were similar across scanners and changed little over time, with a mean intra-scanner coefficient of variation of 3.6 and 1.8%, respectively. The mean inter- and intra-scanner coefficients of variation of brain structure volumes derived from T1-weighted (T1w) images of the human phantom were 3.3 and 1.1%, respectively. The mean inter- and intra-scanner coefficients of variation of FA in white matter regions were 4.5 and 2.5%, while the mean inter- and intra-scanner coefficients of variation of MD in white matter regions were 5.4 and 1.5%.

CONCLUSION

Our results suggest that volumetric and diffusion tensor imaging (DTI) measurements are highly reproducible between and within scanners and provide typical variation amplitudes that can be used as references to interpret future findings in the TACERN network.

White Matter Indices of Medication Response in Major Depression: A Diffusion Tensor Imaging Study

BACKGROUND

While response to antidepressants in major depressive disorder is difficult to predict, characterizing the organization and integrity of white matter in the brain with diffusion tensor imaging (DTI) may provide the means to distinguish between antidepressant responders and nonresponders.

METHODS

DTI data were collected at 6 sites (Canadian Biomarker Integration Network in Depression-1 [CAN-BIND-1 study]) from 200 (127 women) depressed and 112 (71 women) healthy participants at 3 time points: at baseline, 2 weeks, and 8 weeks following initiation of selective serotonin reuptake inhibitor treatment. Therapeutic response was established by a 50% reduction of symptoms at 8 weeks. Analysis on responders, nonresponders, and control subjects yielded 4 scalar metrics: fractional anisotropy and mean, axial, and radial diffusivity. Region-of-interest analysis was carried out on 40 white matter regions using a skeletonization approach. Mixed-effects regression was incorporated to test temporal trends.

RESULTS

The data acquired at baseline showed that axial diffusivity in the external capsule, which overlaps the superior longitudinal fasciculus, was significantly associated with medication response. Regression analysis revealed further baseline differences of responders compared with nonresponders in the cingulum regions, sagittal stratum, and corona radiata. Additional group differences relative to control subjects were seen in the internal capsule, posterior thalamic radiation, and uncinate fasciculus. Most effect sizes were moderate (near 0.5), with a maximum of 0.76 in the cingulum-hippocampus region. No temporal changes in DTI metrics were observed over the 8-week study period.

CONCLUSIONS

Several DTI measures of altered white matter specifically distinguished medication responders and nonresponders at baseline and show promise for predicting treatment response in depression.