MICRA: Microstructural image compilation with repeated acquisitions

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.

Neuroimaging to Facilitate Clinical Trials in Huntington's Disease: Current Opinion from the EHDN Imaging Working Group

 Neuroimaging is increasingly being included in clinical trials of Huntington's disease (HD) for a wide range of purposes from participant selection and safety monitoring, through to demonstration of disease modification. Selection of the appropriate modality and associated analysis tools requires careful consideration. On behalf of the EHDN Imaging Working Group, we present current opinion on the utility and future prospects for inclusion of neuroimaging in HD trials. Covering the key imaging modalities of structural-, functional- and diffusion- MRI, perfusion imaging, positron emission tomography, magnetic resonance spectroscopy, and magnetoencephalography, we address how neuroimaging can be used in HD trials to: 1) Aid patient selection, enrichment, stratification, and safety monitoring; 2) Demonstrate biodistribution, target engagement, and pharmacodynamics; 3) Provide evidence for disease modification; and 4) Understand brain re-organization following therapy. We also present the challenges of translating research methodology into clinical trial settings, including equipment requirements and cost, standardization of acquisition and analysis, patient burden and invasiveness, and interpretation of results. We conclude, that with appropriate consideration of modality, study design and analysis, imaging has huge potential to facilitate effective clinical trials in HD.

Mapping tissue microstructure across the human brain on a clinical scanner with soma and neurite density image metrics

Soma and neurite density image (SANDI) is an advanced diffusion magnetic resonance imaging biophysical signal model devised to probe in vivo microstructural information in the gray matter (GM). This model requires acquisitions that include b values that are at least six times higher than those used in clinical practice. Such high b values are required to disentangle the signal contribution of water diffusing in soma from that diffusing in neurites and extracellular space, while keeping the diffusion time as short as possible to minimize potential bias due to water exchange. These requirements have limited the use of SANDI only to preclinical or cutting-edge human scanners. Here, we investigate the potential impact of neglecting water exchange in the SANDI model and present a 10-min acquisition protocol that enables to characterize both GM and white matter (WM) on 3 T scanners. We implemented analytical simulations to (i) evaluate the stability of the fitting of SANDI parameters when diminishing the number of shells; (ii) estimate the bias due to potential exchange between neurites and extracellular space in such reduced acquisition scheme, comparing it with the bias due to experimental noise. Then, we demonstrated the feasibility and assessed the repeatability and reproducibility of our approach by computing microstructural metrics of SANDI with AMICO toolbox and other state-of-the-art models on five healthy subjects. Finally, we applied our protocol to five multiple sclerosis patients. Results suggest that SANDI is a practical method to characterize WM and GM tissues in vivo on performant clinical scanners.

HYDI-DSI revisited: Constrained non-parametric EAP imaging without q-space re-gridding

Hybrid Diffusion Imaging (HYDI) was one of the first attempts to use multi-shell samplings of the q-space to infer diffusion properties beyond Diffusion Tensor Imaging (DTI) or High Angular Resolution Diffusion Imaging (HARDI). HYDI was intended as a flexible protocol embedding both DTI (for lower b-values) and HARDI (for higher b-values) processing, as well as Diffusion Spectrum Imaging (DSI) when the entire data set was exploited. In the latter case, the spherical sampling of the q-space is re-gridded by interpolation to a Cartesian lattice whose extent covers the range of acquired b-values, hence being acquisition-dependent. The Discrete Fourier Transform (DFT) is afterwards used to compute the corresponding Cartesian sampling of the Ensemble Average Propagator (EAP) in an entirely non-parametric way. From this lattice, diffusion markers such as the Return To Origin Probability (RTOP) or the Mean Squared Displacement (MSD) can be numerically estimated. We aim at re-formulating this scheme by means of a Fourier Transform encoding matrix that eliminates the need for q-space re-gridding at the same time it preserves the non-parametric nature of HYDI-DSI. The encoding matrix is adaptively designed at each voxel according to the underlying DTI approximation, so that an optimal sampling of the EAP can be pursued without being conditioned by the particular acquisition protocol. The estimation of the EAP is afterwards carried out as a regularized Quadratic Programming (QP) problem, which allows to impose positivity constraints that cannot be trivially embedded within the conventional HYDI-DSI. We demonstrate that the definition of the encoding matrix in the adaptive space allows to analytically (as opposed to numerically) compute several popular descriptors of diffusion with the unique source of error being the cropping of high frequency harmonics in the Fourier analysis of the attenuation signal. They include not only RTOP and MSD, but also Return to Axis/Plane Probabilities (RTAP/RTPP), which are defined in terms of specific spatial directions and are not available with the former HYDI-DSI. We report extensive experiments that suggest the benefits of our proposal in terms of accuracy, robustness and computational efficiency, especially when only standard, non-dedicated q-space samplings are available.

Efficient estimation of propagator anisotropy and non-Gaussianity in multishell diffusion MRI with micro-structure adaptive convolution kernels and dual Fourier integral transforms

PURPOSE

We seek to reformulate the so-called Propagator Anisotropy (PA) and Non-Gaussianity (NG), originally conceived for the Mean Apparent Propagator diffusion MRI (MAP-MRI), to the Micro-Structure adaptive convolution kernels and dual Fourier Integral Transforms (MiSFIT). These measures describe relevant normalized features of the Ensemble Average Propagator (EAP).

THEORY AND METHODS

First, the indices, which are defined as the EAP's dissimilarity from an isotropic (PA) or a Gaussian (NG) one, are analytically reformulated within the MiSFIT framework. Then a comparison between the resulting maps is drawn by means of a visual analysis, a quantitative assessment via numerical simulations, a test-retest study across the MICRA dataset (6 subjects scanned five times) and, finally, a computational time evaluation.

RESULTS

Findings illustrate the visual similarity between the indices computed with either technique. Evaluation against synthetic ground truth data, however, demonstrates MiSFIT's improved accuracy. In addition, the test-retest study reveals MiSFIT's higher degree of reliability in most of white matter regions. Finally, the computational time evaluation shows MiSFIT's time reduction up to two orders of magnitude.

CONCLUSIONS

Despite being a direct development on the MAP-MRI representation, the PA and the NG can be reliably and efficiently computed within MiSFIT's framework. This, together with the previous findings in the original MiSFIT's article, could mean the difference that definitely qualifies diffusion MRI to be incorporated into regular clinical settings.

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.

Mapping the Human Connectome using Diffusion MRI at 300 mT/m Gradient Strength: Methodological Advances and Scientific Impact

Tremendous efforts have been made in the last decade to advance cutting-edge MRI technology in pursuit of mapping structural connectivity in the living human brain with unprecedented sensitivity and speed. The first Connectom 3T MRI scanner equipped with a 300 mT/m whole-body gradient system was installed at the Massachusetts General Hospital in 2011 and was specifically constructed as part of the Human Connectome Project. Since that time, numerous technological advances have been made to enable the broader use of the Connectom high gradient system for diffusion tractography and tissue microstructure studies and leverage its unique advantages and sensitivity to resolving macroscopic and microscopic structural information in neural tissue for clinical and neuroscientific studies. The goal of this review article is to summarize the technical developments that have emerged in the last decade to support and promote large-scale and scientific studies of the human brain using the Connectom scanner. We provide a brief historical perspective on the development of Connectom gradient technology and the efforts that led to the installation of three other Connectom 3T MRI scanners worldwide - one in the United Kingdom in Cardiff, Wales, another in Continental Europe in Leipzig, Germany, and the latest in Asia in Shanghai, China. We summarize the key developments in gradient hardware and image acquisition technology that have formed the backbone of Connectom-related research efforts, including the rich array of high-sensitivity receiver coils, pulse sequences, image artifact correction strategies and data preprocessing methods needed to optimize the quality of high-gradient strength dMRI data for subsequent analyses. Finally, we review the scientific impact of the Connectom MRI scanner, including advances in diffusion tractography, tissue microstructural imaging, ex vivo validation, and clinical investigations that have been enabled by Connectom technology. We conclude with brief insights into the unique value of strong gradients for dMRI and where the field is headed in the coming years.

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.

Fiber tractography bundle segmentation depends on scanner effects, vendor effects, acquisition resolution, diffusion sampling scheme, diffusion sensitization, and bundle segmentation workflow

When investigating connectivity and microstructure of white matter pathways of the brain using diffusion tractography bundle segmentation, it is important to understand potential confounds and sources of variation in the process. While cross-scanner and cross-protocol effects on diffusion microstructure measures are well described (in particular fractional anisotropy and mean diffusivity), it is unknown how potential sources of variation effect bundle segmentation results, which features of the bundle are most affected, where variability occurs, nor how these sources of variation depend upon the method used to reconstruct and segment bundles. In this study, we investigate six potential sources of variation, or confounds, for bundle segmentation: variation (1) across scan repeats, (2) across scanners, (3) across vendors (4) across acquisition resolution, (5) across diffusion schemes, and (6) across diffusion sensitization. We employ four different bundle segmentation workflows on two benchmark multi-subject cross-scanner and cross-protocol databases, and investigate reproducibility and biases in volume overlap, shape geometry features of fiber pathways, and microstructure features within the pathways. We find that the effects of acquisition protocol, in particular acquisition resolution, result in the lowest reproducibility of tractography and largest variation of features, followed by vendor-effects, scanner-effects, and finally diffusion scheme and b-value effects which had similar reproducibility as scan-rescan variation. However, confounds varied both across pathways and across segmentation workflows, with some bundle segmentation workflows more (or less) robust to sources of variation. Despite variability, bundle dissection is consistently able to recover the same location of pathways in the deep white matter, with variation at the gray matter/ white matter interface. Next, we show that differences due to the choice of bundle segmentation workflows are larger than any other studied confound, with low-to-moderate overlap of the same intended pathway when segmented using different methods. Finally, quantifying microstructure features within a pathway, we show that tractography adds variability over-and-above that which exists due to noise, scanner effects, and acquisition effects. Overall, these confounds need to be considered when harmonizing diffusion datasets, interpreting or combining data across sites, and when attempting to understand the successes and limitations of different methodologies in the design and development of new tractography or bundle segmentation methods.

Changes in white matter microstructure and MRI-derived cerebral blood flow after 1-week of exercise training

Exercise is beneficial for brain health, inducing neuroplasticity and vascular plasticity in the hippocampus, which is possibly mediated by brain-derived neurotrophic factor (BDNF) levels. Here we investigated the short-term effects of exercise, to determine if a 1-week intervention is sufficient to induce brain changes. Fifteen healthy young males completed five supervised exercise training sessions over seven days. This was preceded and followed by a multi-modal magnetic resonance imaging (MRI) scan (diffusion-weighted MRI, perfusion-weighted MRI, dual-calibrated functional MRI) acquired 1 week apart, and blood sampling for BDNF. A diffusion tractography analysis showed, after exercise, a significant reduction relative to baseline in restricted fraction-an axon-specific metric-in the corpus callosum, uncinate fasciculus, and parahippocampal cingulum. A voxel-based approach found an increase in fractional anisotropy and reduction in radial diffusivity symmetrically, in voxels predominantly localised in the corpus callosum. A selective increase in hippocampal blood flow was found following exercise, with no change in vascular reactivity. BDNF levels were not altered. Thus, we demonstrate that 1 week of exercise is sufficient to induce microstructural and vascular brain changes on a group level, independent of BDNF, providing new insight into the temporal dynamics of plasticity, necessary to exploit the therapeutic potential of exercise.

Detecting microstructural deviations in individuals with deep diffusion MRI tractometry

Most diffusion magnetic resonance imaging studies of disease rely on statistical comparisons between large groups of patients and healthy participants to infer altered tissue states in the brain; however, clinical heterogeneity can greatly challenge their discriminative power. There is currently an unmet need to move away from the current approach of group-wise comparisons to methods with the sensitivity to detect altered tissue states at the individual level. This would ultimately enable the early detection and interpretation of microstructural abnormalities in individual patients, an important step towards personalized medicine in translational imaging. To this end, Detect was developed to advance diffusion magnetic resonance imaging tractometry towards single-patient analysis. By operating on the manifold of white-matter pathways and learning normative microstructural features, our framework captures idiosyncrasies in patterns along white-matter pathways. Our approach paves the way from traditional group-based comparisons to true personalized radiology, taking microstructural imaging from the bench to the bedside.

Temporal Dynamics of Brain White Matter Plasticity in Sighted Subjects during Tactile Braille Learning: A Longitudinal Diffusion Tensor Imaging Study

The white matter (WM) architecture of the human brain changes in response to training, though fine-grained temporal characteristics of training-induced white matter plasticity remain unexplored. We investigated white matter microstructural changes using diffusion tensor imaging at five different time points in 26 sighted female adults during 8 months of training on tactile braille reading. Our results show that training-induced white matter plasticity occurs both within and beyond the trained sensory modality, as reflected by fractional anisotropy (FA) increases in somatosensory and visual cortex, respectively. The observed changes followed distinct time courses, with gradual linear FA increase along the training in the somatosensory cortex and sudden visual cortex cross-modal plasticity occurring after braille input became linguistically meaningful. WM changes observed in these areas returned to baseline after the cessation of learning in line with the supply-demand model of plasticity. These results also indicate that the temporal dynamics of microstructural plasticity in different cortical regions might be modulated by the nature of computational demands. We provide additional evidence that observed FA training-induced changes are behaviorally relevant to tactile reading. Together, these results demonstrate that WM plasticity is a highly dynamic process modulated by the introduction of novel experiences.SIGNIFICANCE STATEMENT Throughout the lifetime the human brain is shaped by various experiences. Training-induced reorganization in white matter (WM) microstructure has been reported, but we know little about its temporal dynamics. To fill this gap, we scanned sighted subjects five times during tactile braille reading training. We observed different dynamics of WM plasticity in the somatosensory and visual cortices implicated in braille reading. The former showed a continuous increase in WM tissue anisotropy along with tactile training, while microstructural changes in the latter were observed only after the participants learned to read braille words. Our results confirm the supply-demand model of brain plasticity and provide evidence that WM reorganization depends on distinct computational demands and functional roles of regions involved in the trained skill.

Scan-rescan repeatability of axonal imaging metrics using high-gradient diffusion MRI and statistical implications for study design

Axon diameter mapping using diffusion MRI in the living human brain has attracted growing interests with the increasing availability of high gradient strength MRI systems. A systematic assessment of the consistency of axon diameter estimates within and between individuals is needed to gain a comprehensive understanding of how such methods extend to quantifying differences in axon diameter index between groups and facilitate the design of neurobiological studies using such measures. We examined the scan-rescan repeatability of axon diameter index estimation based on the spherical mean technique (SMT) approach using diffusion MRI data acquired with gradient strengths up to 300 mT/m on a 3T Connectom system in 7 healthy volunteers. We performed statistical power analyses using data acquired with the same protocol in a larger cohort consisting of 15 healthy adults to investigate the implications for study design. Results revealed a high degree of repeatability in voxel-wise restricted volume fraction estimates and tract-wise estimates of axon diameter index derived from high-gradient diffusion MRI data. On the region of interest (ROI) level, across white matter tracts in the whole brain, the Pearson’s correlation coefficient of the axon diameter index estimated between scan and rescan experiments was r = 0.72 with an absolute deviation of 0.18 μm. For an anticipated 10% effect size in studies of axon diameter index, most white matter regions required a sample size of less than 15 people to observe a measurable difference between groups using an ROI-based approach. To facilitate the use of high-gradient strength diffusion MRI data for neuroscientific studies of axonal microstructure, the comprehensive multi-gradient strength, multi-diffusion time data used in this work will be made publicly available, in support of open science and increasing the accessibility of such data to the greater scientific community.