Identical, but not the same: intra-site and inter-site reproducibility of fractional anisotropy measures on two 3.0T scanners

MidRISH: Unbiased harmonization of rotationally invariant harmonics of the diffusion signal

Data harmonization is necessary for removing confounding effects in multi-site diffusion image analysis. One such harmonization method, LinearRISH, scales rotationally invariant spherical harmonic (RISH) features from one site ("target") to the second ("reference") to reduce confounding scanner effects. However, reference and target site designations are not arbitrary and resultant diffusion metrics (fractional anisotropy, mean diffusivity) are biased by this choice. In this work we propose MidRISH: rather than scaling reference RISH features to target RISH features, we project both sites to a mid-space. We validate MidRISH with the following experiments: harmonizing scanner differences from 37 matched patients free of cognitive impairment, and harmonizing acquisition and study differences on 117 matched patients free of cognitive impairment. We find that MidRISH reduces bias of reference selection while preserving harmonization efficacy of LinearRISH. Users should be cautious when performing LinearRISH harmonization. To select a reference site is to choose diffusion metric effect-size. Our proposed method eliminates the bias-inducing site selection step.

ComBating inter-site differences in field strength: harmonizing preclinical traumatic brain injury MRI data

Integrating datasets from multiple sites and scanners can increase statistical power for neuroimaging studies but can also introduce significant inter-site confounds. We evaluated the effectiveness of ComBat, an empirical Bayes approach, to combine longitudinal preclinical MRI data acquired at 4.7 or 9.4 T at two different sites in Australia. Male Sprague Dawley rats underwent MRI on Days 2, 9, 28, and 150 following moderate/severe traumatic brain injury (TBI) or sham injury as part of Project 1 of the NIH/NINDS-funded Centre Without Walls EpiBioS4Rx project. Diffusion-weighted and multiple-gradient-echo images were acquired, and outcomes included QSM, FA, and ADC. Acute injury measures including apnea and self-righting reflex were consistent between sites. Mixed-effect analysis of ipsilateral and contralateral corpus callosum (CC) summary values revealed a significant effect of site on FA and ADC values, which was removed following ComBat harmonization. Bland-Altman plots for each metric showed reduced variability across sites following ComBat harmonization, including for QSM, despite appearing to be largely unaffected by inter-site differences and no effect of site observed. Following harmonization, the combined inter-site data revealed significant differences in the imaging metrics consistent with previously reported outcomes. TBI resulted in significantly reduced FA and increased susceptibility in the ipsilateral CC, and significantly reduced FA in the contralateral CC compared with sham-injured rats. Additionally, TBI rats also exhibited a reversal in ipsilateral CC ADC values over time with significantly reduced ADC at Day 9, followed by increased ADC 150 days after injury. Our findings demonstrate the need for harmonizing multi-site preclinical MRI data and show that this can be successfully achieved using ComBat while preserving phenotypical changes due to TBI.

A diffusion MRI tractography atlas for concurrent white matter mapping across Eastern and Western populations

The study of brain differences across Eastern and Western populations provides vital insights for understanding potential cultural and genetic influences on cognition and mental health. Diffusion MRI (dMRI) tractography is an important tool in assessing white matter (WM) connectivity and brain tissue microstructure across different populations. However, a comprehensive investigation into WM fiber tracts between Eastern and Western populations is challenged due to the lack of a cross-population WM atlas and the large site-specific variability of dMRI data. This study presents a dMRI tractography atlas, namely the East-West WM Atlas, for concurrent WM mapping between Eastern and Western populations and creates a large, harmonized dMRI dataset (n=306) based on the Human Connectome Project and the Chinese Human Connectome Project. The curated WM atlas, as well as subject-specific data including the harmonized dMRI data, the whole brain tractography data, and parcellated WM fiber tracts and their diffusion measures, are publicly released. This resource is a valuable addition to facilitating the exploration of brain commonalities and differences across diverse cultural backgrounds.

Tractometry of Human Visual White Matter Pathways in Health and Disease

Diffusion-weighted MRI (dMRI) provides a unique non-invasive view of human brain tissue properties. The present review article focuses on tractometry analysis methods that use dMRI to assess the properties of brain tissue within the long-range connections comprising brain networks. We focus specifically on the major white matter tracts that convey visual information. These connections are particularly important because vision provides rich information from the environment that supports a large range of daily life activities. Many of the diseases of the visual system are associated with advanced aging, and tractometry of the visual system is particularly important in the modern aging society. We provide an overview of the tractometry analysis pipeline, which includes a primer on dMRI data acquisition, voxelwise model fitting, tractography, recognition of white matter tracts, and calculation of tract tissue property profiles. We then review dMRI-based methods for analyzing visual white matter tracts: the optic nerve, optic tract, optic radiation, forceps major, and vertical occipital fasciculus. For each tract, we review background anatomical knowledge together with recent findings in tractometry studies on these tracts and their properties in relation to visual function and disease. Overall, we find that measurements of the brain's visual white matter are sensitive to a range of disorders and correlate with perceptual abilities. We highlight new and promising analysis methods, as well as some of the current barriers to progress toward integration of these methods into clinical practice. These barriers, such as variability in measurements between protocols and instruments, are targets for future development.

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.

Implementing ABCD studyⓇ MRI sequences for multi-site cohort studies: Practical guide to necessary steps, preprocessing methods, and challenges

Large multi-site studies that combine magnetic resonance imaging (MRI) data across research sites present exceptional opportunities to advance neuroscience research. However, scanner or site variability and non-standardised image acquisition protocols, data processing and analysis pipelines can adversely affect the reliability and repeatability of MRI derived brain measures. We implemented a standardised MRI protocol based on that used in the Adolescent Brain Cognition Development (ABCD)Ⓡ study in two sites, and across four MRI scanners. Twice repeated measurements of a single healthy volunteer were obtained in two sites and in four 3T MRI scanners (vendors: Siemens, Philips, and GE). Imaging data included anatomical scans (T1 weighted, T2 weighted), diffusion weighted imaging (DWI) and resting state functional MRI (rs-fMRI). Standardised containerized pipelines were utilised to pre-process the data and different image quality metrics and test-retest variability of different brain metrics were evaluated. The implementation of the MRI protocols was possible with minor adjustments in acquisition (e.g. repetition time (TR), higher b-values) and exporting (DICOM formats) of images due to different technical performance of the scanners. This study provides practical insights into the implementation of standardised sequences and data processing for multisite studies, showcase the benefits of containerised preprocessing tools, and highlights the need for careful optimisation of multisite image acquisition.

Diffusion MRI harmonization via personalized template mapping

One fundamental challenge in diffusion magnetic resonance imaging (dMRI) harmonization is to disentangle the contributions of scanner-related effects from the variable brain anatomy for the observed imaging signals. Conventional harmonization methods rely on establishing an atlas space to resolve anatomical variability and generate a unified inter-site mapping function. However, this approach is limited in accounting for the misalignment of neuroanatomy that still widely persists even after registration, especially in regions close to cortical boundaries. To overcome this challenge, we propose a personalized framework in this paper to more effectively address the confounding from the misalignment of neuroanatomy in dMRI harmonization. Instead of using a common template representing site-effects for all subjects, the main novelty of our method is the adaptive computation of personalized templates for both source and target scanning sites to estimate the inter-site mapping function. We integrate our method with the rotation invariant spherical harmonics (RISH) features to achieve the harmonization of dMRI signals between sites. In our experiments, the proposed approach is applied to harmonize the dMRI data acquired from two scanning platforms: Siemens Prisma and GE MR750 from the Adolescent Brain Cognitive Development dataset and compared with a state-of-the-art method based on RISH features. Our results indicate that the proposed harmonization framework achieves superior performance not only in reducing inter-site variations due to scanner differences but also in preserving sex-related biological variability in original cohorts. Moreover, we assess the impact of harmonization on the estimation of fiber orientation distributions and show the robustness of the personalized harmonization procedure in preserving the fiber orientation of original dMRI signals.

Lifespan reference curves for harmonizing multi-site regional brain white matter metrics from diffusion MRI

Age-related white matter (WM) microstructure maturation and decline occur throughout the human lifespan, complementing the process of gray matter development and degeneration. Here, we create normative lifespan reference curves for global and regional WM microstructure by harmonizing diffusion MRI (dMRI)-derived data from ten public datasets (N = 40,898 subjects; age: 3-95 years; 47.6% male). We tested three harmonization methods on regional diffusion tensor imaging (DTI) based fractional anisotropy (FA), a metric of WM microstructure, extracted using the ENIGMA-DTI pipeline. ComBat-GAM harmonization provided multi-study trajectories most consistent with known WM maturation peaks. Lifespan FA reference curves were validated with test-retest data and used to assess the effect of the ApoE4 risk factor for dementia in WM across the lifespan. We found significant associations between ApoE4 and FA in WM regions associated with neurodegenerative disease even in healthy individuals across the lifespan, with regional age-by-genotype interactions. Our lifespan reference curves and tools to harmonize new dMRI data to the curves are publicly available as eHarmonize (https://github.com/ahzhu/eharmonize).

Harmonized diffusion MRI data and white matter measures from the Adolescent Brain Cognitive Development Study

The Adolescent Brain Cognitive Development (ABCD) Study® has collected data from over 10,000 children across 21 sites, providing insights into adolescent brain development. However, site-specific scanner variability has made it challenging to use diffusion MRI (dMRI) data from this study. To address this, a dataset of harmonized and processed ABCD dMRI data (from release 3) has been created, comprising quality-controlled imaging data from 9,345 subjects, focusing exclusively on the baseline session, i.e., the first time point of the study. This resource required substantial computational time (approx. 50,000 CPU hours) for harmonization, whole-brain tractography, and white matter parcellation. The dataset includes harmonized dMRI data, 800 white matter clusters, 73 anatomically labeled white matter tracts in full and low resolution, and 804 different dMRI-derived measures per subject (72.3 TB total size). Accessible via the NIMH Data Archive, it offers a large-scale dMRI dataset for studying structural connectivity in child and adolescent neurodevelopment. Additionally, several post-harmonization experiments were conducted to demonstrate the success of the harmonization process on the ABCD dataset.

Traveling Subject-Informed Harmonization Increases Reliability of Brain Diffusion Tensor and Neurite Mapping

Diffusion-weighted magnetic resonance imaging (dMRI) of brain has helped elucidate the microstructural changes of psychiatric and neurodegenerative disorders. Inconsistency between MRI models has hampered clinical application of dMRI-based metrics. Using harmonized dMRI data of 300 scans from 69 traveling subjects (TS) scanning the same individuals at multiple conditions with 13 MRI models and 2 protocols, the widely-used metrics such as diffusion tensor imaging (DTI) and neurite orientation dispersion and density imaging (NODDI) were evaluated before and after harmonization with a combined association test (ComBat) or TS-based general linear model (TS-GLM). Results showed that both ComBat and TS-GLM significantly reduced the effects of the MRI site, model, and protocol for diffusion metrics while maintaining the intersubject biological effects. The harmonization power of TS-GLM based on TS data model is more powerful than that of ComBat. In conclusion, our research demonstrated that although ComBat and TS-GLM harmonization approaches were effective at reducing the scanner effects of the site, model, and protocol for DTI and NODDI metrics in WM, they exhibited high retainability of biological effects. Therefore, we suggest that, after harmonizing DTI and NODDI metrics, a multisite study with large cohorts can accurately detect small pathological changes by retaining pathological effects.

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

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

MidRISH: Unbiased harmonization of rotationally invariant harmonics of the diffusion signal

Objective

Data harmonization is necessary for removing confounding effects in multi-site diffusion image analysis. One such harmonization method, LinearRISH, scales rotationally invariant spherical harmonic (RISH) features from one site ("target") to the second ("reference") to reduce confounding scanner effects. However, reference and target site designations are not arbitrary and resultant diffusion metrics (fractional anisotropy, mean diffusivity) are biased by this choice. In this work we propose MidRISH: rather than scaling reference RISH features to target RISH features, we project both sites to a mid-space.

Methods

We validate MidRISH with the following experiments: harmonizing scanner differences from 37 matched patients free of cognitive impairment, and harmonizing acquisition and study differences on 117 matched patients free of cognitive impairment.

Conclusion

MidRISH reduces bias of reference selection while preserving harmonization efficacy of LinearRISH.

Significance

Users should be cautious when performing LinearRISH harmonization. To select a reference site is to choose diffusion metric effect-size. Our proposed method eliminates the bias-inducing site selection step.

versaFlow: a versatile pipeline for resolution adapted diffusion MRI processing and its application to studying the variability of the PRIME-DE database

The lack of "gold standards" in Diffusion Weighted Imaging (DWI) makes validation cumbersome. To tackle this task, studies use translational analysis where results in humans are benchmarked against findings in other species. Non-Human Primates (NHP) are particularly interesting for this, as their cytoarchitecture is closely related to humans. However, tools used for processing and analysis must be adapted and finely tuned to work well on NHP images. Here, we propose versaFlow, a modular pipeline implemented in Nextflow, designed for robustness and scalability. The pipeline is tailored to in vivo NHP DWI at any spatial resolution; it allows for maintainability and customization. Processes and workflows are implemented using cutting-edge and state-of-the-art Magnetic Resonance Imaging (MRI) processing technologies and diffusion modeling algorithms, namely Diffusion Tensor Imaging (DTI), Constrained Spherical Deconvolution (CSD), and DIstribution of Anisotropic MicrOstructural eNvironments in Diffusion-compartment imaging (DIAMOND). Using versaFlow, we provide an in-depth study of the variability of diffusion metrics computed on 32 subjects from 3 sites of the Primate Data Exchange (PRIME-DE), which contains anatomical T1-weighted (T1w) and T2-weighted (T2w) images, functional MRI (fMRI), and DWI of NHP brains. This dataset includes images acquired over a range of resolutions, using single and multi-shell gradient samplings, on multiple scanner vendors. We perform a reproducibility study of the processing of versaFlow using the Aix-Marseilles site's data, to ensure that our implementation has minimal impact on the variability observed in subsequent analyses. We report very high reproducibility for the majority of metrics; only gamma distribution parameters of DIAMOND display less reproducible behaviors, due to the absence of a mechanism to enforce a random number seed in the software we used. This should be taken into consideration when future applications are performed. We show that the PRIME-DE diffusion data exhibits a great level of variability, similar or greater than results obtained in human studies. Its usage should be done carefully to prevent instilling uncertainty in statistical analyses. This hints at a need for sufficient harmonization in acquisition protocols and for the development of robust algorithms capable of managing the variability induced in imaging due to differences in scanner models and/or vendors.

Anesthetic modulation of water diffusion: Insights from a diffusion tensor imaging study

Diffusion tensor imaging (DTI) in animal models are essential for translational neuroscience studies. A critical step in animal studies is the use of anesthetics. Understanding the influence of specific anesthesia regimes on DTI-derived parameters, such as fractional anisotropy (FA) and mean diffusivity (MD), is imperative when comparing results between animal studies using different anesthetics. Here, the quantification of FA and MD under different anesthetic regimes, alpha-chloralose and isoflurane, is discussed. We also used a range of b-values to determine whether the anesthetic effect was b-value dependent. The first group of rats (n = 6) was anesthetized with alpha-chloralose (80 mg/kg), whereas the second group of rats (n = 7) was anesthetized with isoflurane (1.5%). DTI was performed with b-values of 500, 1500, and 1500s/mm2, and the MD and FA were assessed individually. Anesthesia-specific differences in MD were apparent, as manifested by the higher estimated MD under isoflurane anesthesia than that under alpha-chloralose anesthesia (P < 0.001). MD values increased with decreasing b-value in all regions studied, and the degree of increase when rats were anesthetized with isoflurane was more pronounced than that associated with alpha-chloralose (P < 0.05). FA quantitation was also influenced by anesthesia regimens to varying extents, depending on the brain regions and b-values. In conclusion, both scanning parameters and the anesthesia regimens significantly impacted the quantification of DTI indices.

Regional microstructural differences in ADHD youth with and without a family history of bipolar I disorder

BACKGROUND

Having a first-degree relative with bipolar I disorder (BD) in conjunction with prodromal attention deficit/hyperactivity disorder (ADHD) may represent a unique phenotype that confers greater risk for developing BD than ADHD alone. However, underlying neuropathoetiological mechanisms remain poorly understood. This cross-sectional study compared regional microstructure in psychostimulant-free ADHD youth with ('high-risk', HR) and without ('low-risk', LR) a first-degree relative with BD, and healthy controls (HC).

METHODS

A total of 140 (high-risk, n = 44; low-risk, n = 49; and HC, n = 47) youth (mean age: 14.1 ± 2.5 years, 65 % male) were included in the analysis. Diffusion tensor images were collected and fractional anisotropy (FA) and mean diffusivity (MD) maps were calculated. Both tract-based and voxel-based analyses were performed. Correlations between clinical ratings and microstructural metrics that differed among groups were examined.

RESULTS

No significant group differences in major long-distance fiber tracts were observed. The high-risk ADHD group exhibited predominantly higher FA and lower MD in frontal, limbic, and striatal subregions compared with the low-risk ADHD group. Both low-risk and high-risk ADHD groups exhibited higher FA in unique and overlapping regions compared with HC subjects. Significant correlations between regional microstructural metrics and clinical ratings were observed in ADHD groups.

LIMITATIONS

Prospective longitudinal studies will be required to determine the relevance of these findings to BD risk progression.

CONCLUSIONS

Psychostimulant-free ADHD youth with a BD family history exhibit different microstructure alterations in frontal, limbic, and striatal regions compared with ADHD youth without a BD family history, and may therefore represent a unique phenotype relevant to BD risk progression.

RELIEF: A structured multivariate approach for removal of latent inter-scanner effects

Combining data collected from multiple study sites is becoming common and is advantageous to researchers to increase the generalizability and replicability of scientific discoveries. However, at the same time, unwanted inter-scanner biases are commonly observed across neuroimaging data collected from multiple study sites or scanners, rendering difficulties in integrating such data to obtain reliable findings. While several methods for handling such unwanted variations have been proposed, most of them use univariate approaches that could be too simple to capture all sources of scanner-specific variations. To address these challenges, we propose a novel multivariate harmonization method called RELIEF (REmoval of Latent Inter-scanner Effects through Factorization) for estimating and removing both explicit and latent scanner effects. Our method is the first approach to introduce the simultaneous dimension reduction and factorization of interlinked matrices to a data harmonization context, which provides a new direction in methodological research for correcting inter-scanner biases. Analyzing diffusion tensor imaging (DTI) data from the Social Processes Initiative in Neurobiology of the Schizophrenia (SPINS) study and conducting extensive simulation studies, we show that RELIEF outperforms existing harmonization methods in mitigating inter-scanner biases and retaining biological associations of interest to increase statistical power. RELIEF is publicly available as an R package.

Reproducibility of diffusion tensor imaging-derived parameters: implications for the streptozotocin-induced type 1 diabetic rats

OBJECTIVE

Diffusion tensor imaging (DTI) is a useful approach for studying neuronal integrity in animals. However, the test-retest reproducibility of DTI techniques in animals has not been discussed. Therefore, the first part of this work was to systematically elucidate the reliability of DTI-derived parameters in an animal study. Subsequently, we applied the DTI approach to an animal model of diabetes in a longitudinal manner.

MATERIALS AND METHODS

In Study 1, nine rats underwent two DTI sessions using the same scanner and protocols, with a gap of 4 weeks. The reliability of the DTI-derived parameters was evaluated in terms of sessions and raters. In Study 2, nine rats received a single intraperitoneal injection of 70 mg/kg streptozotocin (STZ) to develop diabetes. Longitudinal DTI scans were used to assess brain alterations before and 4 weeks after STZ administration.

RESULTS

In the test-retest evaluation, the inter-scan coefficient of variation (CoV) ranged from 3.04 to 3.73% and 2.12-2.59% for fractional anisotropy (FA) and mean diffusivity (MD), respectively, in different brain regions, suggesting excellent reproducibility. Moreover, rater-dependence had minimal effects on FA and MD quantification, with all inter-rater CoV values less than 4%. Following the onset of diabetes, FA in striatum and cortex were noted to be significantly lower relative to the period where they had not developed diabetes (both P < 0.05). However, when compared to the control group, a significant change in FA caused by diabetes was detected only in the striatum (P < 0.05), but not in the cortex.

CONCLUSION

These results demonstrate good inter-rater and inter-scan reliability of DTI in animal studies, and the longitudinal setting has a beneficial effect on detecting small changes in the brain due to diseases.

A comprehensive evaluation of multicentric reliability of single-subject cortical morphological networks on traveling subjects

Despite the prevalence of research on single-subject cerebral morphological networks in recent years, whether they can offer a reliable way for multicentric studies remains largely unknown. Using two multicentric datasets of traveling subjects, this work systematically examined the inter-site test-retest (TRT) reliabilities of single-subject cerebral morphological networks, and further evaluated the effects of several key factors. We found that most graph-based network measures exhibited fair to excellent reliabilities regardless of different analytical pipelines. Nevertheless, the reliabilities were affected by choices of morphological index (fractal dimension > sulcal depth > gyrification index > cortical thickness), brain parcellation (high-resolution > low-resolution), thresholding method (proportional > absolute), and network type (binarized > weighted). For the factor of similarity measure, its effects depended on the thresholding method used (absolute: Kullback-Leibler divergence > Jensen-Shannon divergence; proportional: Jensen-Shannon divergence > Kullback-Leibler divergence). Furthermore, longer data acquisition intervals and different scanner software versions significantly reduced the reliabilities. Finally, we showed that inter-site reliabilities were significantly lower than intra-site reliabilities for single-subject cerebral morphological networks. Altogether, our findings propose single-subject cerebral morphological networks as a promising approach for multicentric human connectome studies, and offer recommendations on how to determine analytical pipelines and scanning protocols for obtaining reliable results.

Harmonized diffusion MRI data and white matter measures from the Adolescent Brain Cognitive Development Study

The Adolescent Brain Cognitive Development (ABCD) study has collected data from over 10,000 children across 21 sites, providing valuable insights into adolescent brain development. However, site-specific scanner variability has made it challenging to use diffusion MRI (dMRI) data from this study. To address this, a database of harmonized and processed ABCD dMRI data has been created, comprising quality-controlled imaging data from 9345 subjects. This resource required significant computational effort, taking ~50,000 CPU hours to harmonize the data, perform white matter parcellation, and run whole brain tractography. The database includes harmonized dMRI data, 800 white matter clusters, 73 anatomically labeled white matter tracts both in full-resolution (for analysis) and low-resolution (for visualization), and 804 different dMRI-derived measures per subject. It is available via the NIMH Data Archive and offers tremendous potential for scientific discoveries in structural connectivity studies of neurodevelopment in children and adolescents. Additionally, several post-harmonization experiments were conducted to demonstrate the success of the harmonization process on the ABCD dataset.

Meta-analysis of the normal diffusion tensor imaging values of the peripheral nerves in the upper limb

Peripheral neuropathy affects 1 in 10 adults over the age of 40 years. Given the absence of a reliable diagnostic test for peripheral neuropathy, there has been a surge of research into diffusion tensor imaging (DTI) because it characterises nerve microstructure and provides reproducible proxy measures of myelination, axon diameter, fibre density and organisation. Before researchers and clinicians can reliably use diffusion tensor imaging to assess the 'health' of the major nerves of the upper limb, we must understand the "normal" range of values and how they vary with experimental conditions. We searched PubMed, Embase, medRxiv and bioRxiv for studies which reported the findings of DTI of the upper limb in healthy adults. Four review authors independently triple extracted data. Using the meta suite of Stata 17, we estimated the normal fractional anisotropy (FA) and diffusivity (mean, MD; radial, RD; axial AD) values of the median, radial and ulnar nerve in the arm, elbow and forearm. Using meta-regression, we explored how DTI metrics varied with age and experimental conditions. We included 20 studies reporting data from 391 limbs, belonging to 346 adults (189 males and 154 females, ~ 1.2 M:1F) of mean age 34 years (median 31, range 20-80). In the arm, there was no difference in the FA (pooled mean 0.59 mm²/s [95% CI 0.57, 0.62]; I² 98%) or MD (pooled mean 1.13 × 10^-3 mm²/s [95% CI 1.08, 1.18]; I² 99%) of the median, radial and ulnar nerves. Around the elbow, the ulnar nerve had a 12% lower FA than the median and radial nerves (95% CI - 0.25, 0.00) and significantly higher MD, RD and AD. In the forearm, the FA (pooled mean 0.55 [95% CI 0.59, 0.64]; I² 96%) and MD (pooled mean 1.03 × 10^-3 mm²/s [95% CI 0.94, 1.12]; I² 99%) of the three nerves were similar. Multivariable meta regression showed that the b-value, TE, TR, spatial resolution and age of the subject were clinically important moderators of DTI parameters in peripheral nerves. We show that subject age, as well as the b-value, TE, TR and spatial resolution are important moderators of DTI metrics from healthy nerves in the adult upper limb. The normal ranges shown here may inform future clinical and research studies.

Inter-vendor and inter-observer reliability of diffusion tensor imaging in the musculoskeletal system: a multiscanner MR study

BACKGROUND

To evaluate the inter-observer and inter-vendor reliability of diffusion tensor imaging parameters in the musculoskeletal system.

METHODS

This prospective study included six healthy volunteers three men (mean age: 42; range: 31-52 years) and three women (mean age: 36; range: 30-44 years). Each subject was scanned using different 3 Tesla magnetic resonance scanners from three different vendors at three different sites bilaterally. First, the intra-class correlation coefficient was used to determine between-observers agreement for overall measurements and clinical sites. Next, between-group comparisons were made through the nonparametric Friedman's test. Finally, the Bland-Altman method was used to determine agreement among the three scanner measurements, comparing them two by two.

RESULTS

A total of 792 measurement were calculated. ICC reported high levels of agreement between the two observers. ICC related to MD, FA, and RD measurements ranged from 0.88 (95% CI 0.85-0.90) to 0.95 (95% CI 0.94-0.96), from 0.85 (95% CI 0.81-0.88) to 0.95 (95% CI 0.93-0.96), and from 0.89 (0.85-0.90) to 0.92 (0.90-0.94). No statistically significant inter-vendor differences were observed. The Bland-Altmann method confirmed a high correlation between parameter values.

CONCLUSION

An excellent inter-observer and inter-vendor reliability was found in our study.

The comparison of diffusion tensor imaging in human hearts between 1.5 T and 3.0 T

BACKGROUND

The aim was to compare the diffusion tensor imaging (DTI) indices derived from human hearts between 1.5 T and 3.0 T scanners. Additionally, the reproducibility of DTI indices was assessed between 1.5 T and 3.0 T scanners.

METHODS

A total of 18 ex-vivo hearts were derived from patients who underwent heart transplantation. The DTI schemes were performed at 1.5 T and 3.0 T, respectively. Then, the same slices from each ex-vivo heart were selected for image analysis. The student's t-test or Wilcoxon-rank test was used to compare the statistical differences. The agreement of DTI indices was mainly reported as the interclass correlation coefficient (ICC).

RESULTS

No significant differences (all P > 0.05) were found in the DTI indices between 1.5 T and 3.0 T scanners. Interestingly, the ICC of all DTI indices was relatively lower with a low b-value. The reproducibility of the helix angle (HA) was relatively lower when compared to the other DTI indices.

CONCLUSION

The DTI indices of ex-vivo human hearts between 1.5 T and 3.0 T scanners had no significant differences. The consistency of DTI indices needed caution using a low b-value with different field strengths, and the relatively low reproducibility of HA should be considered.

Free water corrected diffusion tensor imaging discriminates between good and poor outcomes of comatose patients after cardiac arrest

OBJECTIVES

Approximately 50% of comatose patients after cardiac arrest never regain consciousness. Cerebral ischaemia may lead to cytotoxic and/or vasogenic oedema, which can be detected by diffusion tensor imaging (DTI). Here, we evaluate the potential value of free water corrected mean diffusivity (MD) and fractional anisotropy (FA) based on DTI, for the prediction of neurological recovery of comatose patients after cardiac arrest.

METHODS

A total of 50 patients after cardiac arrest were included in this prospective cohort study in two Dutch hospitals. DTI was obtained 2-4 days after cardiac arrest. Outcome was assessed at 6 months, dichotomised as poor (cerebral performance category 3-5; n = 20) or good (n = 30) neurological outcome. We calculated the whole brain mean MD and FA and compared between patients with good and poor outcomes. In addition, we compared a preliminary prediction model based on clinical parameters with or without the addition of MD and FA.

RESULTS

We found significant differences between patients with good and poor outcome of mean MD (good: 726 [702-740] × 10^-6 mm²/s vs. poor: 663 [575-736] × 10^-6 mm²/s; p = 0.01) and mean FA (0.30 ± 0.03 vs. 0.28 ± 0.03; p = 0.03). An exploratory prediction model combining clinical parameters, MD and FA increased the sensitivity for reliable prediction of poor outcome from 60 to 85%, compared to the model containing clinical parameters only, but confidence intervals are overlapping.

CONCLUSIONS

Free water-corrected MD and FA discriminate between patients with good and poor outcomes after cardiac arrest and hold the potential to add to multimodal outcome prediction.

KEY POINTS

• Whole brain mean MD and FA differ between patients with good and poor outcome after cardiac arrest. • Free water-corrected MD can better discriminate between patients with good and poor outcome than uncorrected MD. • A combination of free water-corrected MD (sensitive to grey matter abnormalities) and FA (sensitive to white matter abnormalities) holds potential to add to the prediction of outcome.

Aberrant claustrum structure in preterm-born neonates: an MRI study

The human claustrum is a gray matter structure in the white matter between insula and striatum. Previous analysis found altered claustrum microstructure in very preterm-born adults associated with lower cognitive performance. As the claustrum development is related to hypoxia-ischemia sensitive transient cell populations being at-risk in premature birth, we hypothesized that claustrum structure is already altered in preterm-born neonates. We studied anatomical and diffusion-weighted MRIs of 83 preterm- and 83 term-born neonates at term-equivalent age. Additionally, claustrum development was analyzed both in a spectrum of 377 term-born neonates and longitudinally in 53 preterm-born subjects. Data was provided by the developing Human Connectome Project. Claustrum development showed increasing volume, increasing fractional anisotropy (FA), and decreasing mean diffusivity (MD) around term both across term- and preterm-born neonates. Relative to term-born ones, preterm-born neonates had (i) increased absolute and relative claustrum volumes, both indicating increased cellular and/or extracellular matter and being in contrast to other subcortical gray matter regions of decreased volumes such as thalamus; (ii) lower claustrum FA and higher claustrum MD, pointing at increased extracellular matrix and impaired axonal integrity; and (iii) aberrant covariance between claustrum FA and MD, respectively, and that of distributed gray matter regions, hinting at relatively altered claustrum microstructure. Results together demonstrate specifically aberrant claustrum structure in preterm-born neonates, suggesting altered claustrum development in prematurity, potentially relevant for later cognitive performance.

Higher reliability and validity of Wavelet-ALFF of resting-state fMRI: From multicenter database and application to rTMS modulation

Amplitude of low-frequency fluctuation (ALFF) has been widely used for localization of abnormal activity at the single-voxel level in resting-state fMRI (RS-fMRI) studies. However, previous ALFF studies were based on fast Fourier transform (FFT-ALFF). Our recent study found that ALFF based on wavelet transform (Wavelet-ALFF) showed better sensitivity and reproducibility than FFT-ALFF. The current study aimed to test the reliability and validity of Wavelet-ALFF, and apply Wavelet-ALFF to investigate the modulation effect of repetitive transcranial magnetic stimulation (rTMS). The reliability and validity were assessed on multicenter RS-fMRI datasets under eyes closed (EC) and eyes open (EO) conditions (248 healthy participants in total). We then detected the sensitivity of Wavelet-ALFF using a rTMS modulation dataset (24 healthy participants). For each dataset, Wavelet-ALFF based on five mother wavelets (i.e., db2, bior4.4, morl, meyr and sym3) and FFT-ALFF were calculated in the conventional band and five frequency sub-bands. The results showed that the reliability of both inter-scanner and intra-scanner was higher with Wavelet-ALFF than with FFT-ALFF across multiple frequency bands, especially db2-ALFF in the higher frequency band slow-2 (0.1992-0.25 Hz). In terms of validity, the multicenter ECEO datasets showed that the effect sizes of Wavelet-ALFF with all mother wavelets (especially for db2-ALFF) were larger than those of FFT-ALFF across multiple frequency bands. Furthermore, Wavelet-ALFF detected a larger modulation effect than FFT-ALFF. Collectively, Wavelet db2-ALFF showed the best reliability and validity, suggesting that db2-ALFF may offer a powerful metric for inspecting regional spontaneous brain activities in future studies.

Contrastive semi-supervised harmonization of single-shell to multi-shell diffusion MRI

Diffusion weighted MRI (DW-MRI) harmonization is necessary for multi-site or multi-acquisition studies. Current statistical methods address the need to harmonize from one site to another, but do not simultaneously consider the use of multiple datasets which are comprised of multiple sites, acquisitions protocols, and age demographics. This work explores deep learning methods which can generalize across these variations through semi-supervised and unsupervised learning while also learning to estimate multi-shell data from single-shell data using the Multi-shell Diffusion MRI Harmonization Challenge (MUSHAC) and Baltimore Longitudinal Study on Aging (BLSA) datasets. We compare disentanglement harmonization models, which seek to encode anatomy and acquisition in separate latent spaces, and a CycleGAN harmonization model, which uses generative adversarial networks (GAN) to perform style transfer between sites, to the baseline preprocessing and to SHORE interpolation. We find that the disentanglement models achieve superior performance in harmonizing all data while at the same transforming the input data to a single target space across several diffusion metrics (fractional anisotropy, mean diffusivity, mean kurtosis, primary eigenvector).

Improved sensitivity and precision in multicentre diffusion MRI network analysis using thresholding and harmonization

PURPOSE

To investigate if network thresholding and raw data harmonization improve consistency of diffusion MRI (dMRI)-based brain networks while also increasing precision and sensitivity to detect disease effects in multicentre datasets.

METHODS

Brain networks were reconstructed from dMRI of five samples with cerebral small vessel disease (SVD; 629 patients, 166 controls), as a clinically relevant exemplar condition for studies on network integrity. We evaluated consistency of network architecture in age-matched controls, by calculating cross-site differences in connection probability and fractional anisotropy (FA). Subsequently we evaluated precision and sensitivity to disease effects by identifying connections with low FA in sporadic SVD patients relative to controls, using more severely affected patients with a pure form of genetically defined SVD as reference.

RESULTS

In controls, thresholding and harmonization improved consistency of network architecture, minimizing cross-site differences in connection probability and FA. In patients relative to controls, thresholding improved precision to detect disrupted connections by removing false positive connections (precision, before: 0.09-0.19; after: 0.38-0.70). Before harmonization, sensitivity was low within individual sites, with few connections surviving multiple testing correction (k = 0-25 connections). Harmonization and pooling improved sensitivity (k = 38), while also achieving higher precision when combined with thresholding (0.97).

CONCLUSION

We demonstrated that network consistency, precision and sensitivity to detect disease effects in SVD are improved by thresholding and harmonization. We recommend introducing these techniques to leverage large existing multicentre datasets to better understand the impact of disease on brain networks.

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.

Personalized dMRI Harmonization on Cortical Surface

The inter-site variability of diffusion magnetic resonance imaging (dMRI) hinders the aggregation of dMRI data from multiple centers. This necessitates dMRI harmonization for removing non-biological site-effects. Recently, the emergence of high-resolution dMRI data across various connectome imaging studies allows the large-scale analysis of cortical micro-structure. Existing harmonization methods, however, perform poorly in the harmonization of dMRI data in cortical areas because they rely on image registration methods to factor out anatomical variations, which have known difficulty in aligning cortical folding patterns. To overcome this fundamental challenge in dMRI harmonization, we propose a framework of personalized dMRI harmonization on the cortical surface to improve the dMRI harmonization of gray matter by adaptively estimating the inter-site harmonization mappings. In our experiments, we demonstrate the effectiveness of the proposed method by applying it to harmonize dMRI across the Human Connectome Project (HCP) and the Lifespan Human Connectome Projects in Development (HCPD) studies and achieved much better performance in comparison with conventional methods based on image registration.

Does Loss of Integrity of the Cingulum Bundle Link Amyloid-β Accumulation and Neurodegeneration in Alzheimer’s Disease?

Background: Alzheimer’s disease is characterized by the accumulation of amyloid-β (Aβ) into plaques, aggregation of tau into neurofibrillary tangles, and neurodegenerative processes including atrophy. However, there is a poorly understood spatial discordance between initial Aβ deposition and local neurodegeneration. Objective: Here, we test the hypothesis that the cingulum bundle links Aβ deposition in the cingulate cortex to medial temporal lobe (MTL) atrophy. Methods: 21 participants with mild cognitive impairment (MCI) from the UMC Utrecht memory clinic (UMCU, discovery sample) and 37 participants with MCI from Alzheimer’s Disease Neuroimaging Initiative (ADNI, replication sample) with available Aβ-PET scan, T1-weighted and diffusion-weighted MRI were included. Aβ load of the cingulate cortex was measured by the standardized uptake value ratio (SUVR), white matter integrity of the cingulum bundle was assessed by mean diffusivity and atrophy of the MTL by normalized MTL volume. Relationships were tested with linear mixed models, to accommodate multiple measures for each participant. Results: We found at most a weak association between cingulate Aβ and MTL volume (added R2 <0.06), primarily for the posterior hippocampus. In neither sample, white matter integrity of the cingulum bundle was associated with cingulate Aβ or MTL volume (added R2 <0.01). Various sensitivity analyses (Aβ-positive individuals only, posterior cingulate SUVR, MTL sub region volume) provided similar results. Conclusion: These findings, consistent in two independent cohorts, do not support our hypothesis that loss of white matter integrity of the cingulum is a connecting factor between cingulate gyrus Aβ deposition and MTL atrophy.

A multivariate Bayesian classification algorithm for cerebral stage prediction by diffusion tensor imaging in amyotrophic lateral sclerosis

•

Novel DTI-based classification of ALS disease stages by a Bayesian approach is applied.

•

Bayesian classification algorithm improves threshold-based staging method.

•

Step forward in MRI-based patient stratification in ALS in vivo.

Background and Objective

Diffusion tensor imaging (DTI) can be used to tract-wise map correlates of the sequential disease progression and, therefore, to assess disease stages of amyotrophic lateral sclerosis (ALS) in vivo. According to a threshold-based sequential scheme, a classification of ALS patients into disease stages is possible, however, several patients cannot be staged for methodological reasons. This study aims to implement a multivariate Bayesian classification algorithm for disease stage prediction at an individual ALS patient level based on DTI metrics of involved tract systems to improve disease stage mapping.

Methods

The analysis of fiber tracts involved in each stage of ALS was performed in 325 ALS patients and 130 age- and gender-matched healthy controls. Based on Bayes’ theorem and in accordance with the sequential disease progression, a multistage classifier was implemented. Patients were categorized into in vivo DTI stages using the threshold-based method and the Bayesian algorithm. By the margin of confidence, the reliability of the Bayesian categorizations was accessible.

Results

Based on the Bayesian multistage classifier, 88% of all ALS patients could be assigned into an ALS stage compared to 77% using the threshold-based staging scheme. Additionally, the confidence of all classifications could be estimated.

Conclusions

By the application of the multi-stage Bayesian classifier, an individualized in vivo cerebral staging of ALS patients was possible based on the sequentially involved tract systems and, furthermore, the reliability of the respective classifications could be determined. The Bayesian classification algorithm is an improvement of the threshold-based staging method and could provide a framework for extending the DTI-based in vivo cerebral staging in ALS.

Parameters from site classification to harmonize MRI clinical studies: Application to a multi-site Parkinson's disease dataset

Multi‐site MRI datasets are crucial for big data research. However, neuroimaging studies must face the batch effect. Here, we propose an approach that uses the predictive probabilities provided by Gaussian processes (GPs) to harmonize clinical‐based studies. A multi‐site dataset of 216 Parkinson's disease (PD) patients and 87 healthy subjects (HS) was used. We performed a site GP classification using MRI data. The outcomes estimated from this classification, redefined like Weighted HARMonization PArameters (WHARMPA), were used as regressors in two different clinical studies: A PD versus HS machine learning classification using GP, and a VBM comparison (FWE‐p < .05, k = 100). Same studies were also conducted using conventional Boolean site covariates, and without information about site belonging. The results from site GP classification provided high scores, balanced accuracy (BAC) was 98.39% for grey matter images. PD versus HS classification performed better when the WHARMPA were used to harmonize (BAC = 78.60%; AUC = 0.90) than when using the Boolean site information (BAC = 56.31%; AUC = 0.71) and without it (BAC = 57.22%; AUC = 0.73). The VBM analysis harmonized using WHARMPA provided larger and more statistically robust clusters in regions previously reported in PD than when the Boolean site covariates or no corrections were added to the model. In conclusion, WHARMPA might encode global site‐effects quantitatively and allow the harmonization of data. This method is user‐friendly and provides a powerful solution, without complex implementations, to clean the analyses by removing variability associated with the differences between sites.

Brain Imaging Biomarkers for Chronic Pain

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

Potential Pitfalls of Using Fractional Anisotropy, Axial Diffusivity, and Radial Diffusivity as Biomarkers of Cerebral White Matter Microstructure

Fractional anisotropy (FA), axial diffusivity (AD), and radial diffusivity (RD) are commonly used as MRI biomarkers of white matter microstructure in diffusion MRI studies of neurodevelopment, brain aging, and neurologic injury/disease. Some of the more frequent practices include performing voxel-wise or region-based analyses of these measures to cross-sectionally compare individuals or groups, longitudinally assess individuals or groups, and/or correlate with demographic, behavioral or clinical variables. However, it is now widely recognized that the majority of cerebral white matter voxels contain multiple fiber populations with different trajectories, which renders these metrics highly sensitive to the relative volume fractions of the various fiber populations, the microstructural integrity of each constituent fiber population, and the interaction between these factors. Many diffusion imaging experts are aware of these limitations and now generally avoid using FA, AD or RD (at least in isolation) to draw strong reverse inferences about white matter microstructure, but based on the continued application and interpretation of these metrics in the broader biomedical/neuroscience literature, it appears that this has perhaps not yet become common knowledge among diffusion imaging end-users. Therefore, this paper will briefly discuss the complex biophysical underpinnings of these measures in the context of crossing fibers, provide some intuitive “thought experiments” to highlight how conventional interpretations can lead to incorrect conclusions, and suggest that future studies refrain from using (over-interpreting) FA, AD, and RD values as standalone biomarkers of cerebral white matter microstructure.

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.

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.

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.

Reproducibility of Longitudinal Changes in Cortical Thickness Determined by Surface-Based Morphometry Between Non-Accelerated and Accelerated MR Imaging

BACKGROUND

Scan acceleration such as parallel imaging reduces scan time, but shorter scan time may reduce the signal-to-noise ratio and affect image quality. The reproducibility of longitudinal changes in the brain structure between non-accelerated and accelerated imaging by surface-based analysis is unclear.

PURPOSE

To determine the reproducibility of longitudinal changes in cortical thickness, measured by surface-based morphometry, between non-accelerated and accelerated structural T₁ -weighted imaging in the healthy elderly and those with mild cognitive impairment (MCI) and Alzheimer's disease (AD).

STUDY TYPE

Retrospective.

SUBJECTS

Fifty healthy elderly subjects (age = 73 ± 5 years, 29 females, 21 males), 54 MCI patients (age = 71 ± 7 years, 23 females, 31 males), and 8 AD patients (age = 78 ± 6 years, 6 females, 2 males).

FIELD STRENGTH/SEQUENCE

3 T, magnetization-prepared rapid gradient-echo.

ASSESSMENT

Longitudinal changes in cortical thickness estimated by the longitudinal stream in FreeSurfer from 2-year interval data, and visual assessment of image quality by three radiologists.

STATISTICAL TESTS

Intraclass correlation coefficient (ICC) and Kruskal-Wallis test. A P value <0.05 was considered significant.

RESULTS

Healthy elderly subjects, MCI patients, and AD patients showed different patterns in the ICC maps. For the smoothing of 20 mm full width at half maximum, the mean ICC was 0.45 overall (healthy elderly, 0.33; MCI patients, 0.49; AD patients, 0.31). The within-subject SDs of the symmetrized percent changes were similar between healthy elderly subjects (mean, 1.3%/year) and MCI patients (mean, 1.3%/year) but larger in AD patients (mean, 1.7%/year). Image quality did not significantly differ per group (P = 0.18).

DATA CONCLUSION

The results of this study indicate the reproducibility of longitudinal changes in cortical thickness measured by surface-based morphometry between non-accelerated and accelerated imaging, and that the reproducibility varies by disease and region.

LEVEL OF EVIDENCE

3 TECHNICAL EFFICACY: Stage 1.

Reproducibility of Brain Volume Changes in Longitudinal Voxel-Based Morphometry Between Non-Accelerated and Accelerated Magnetic Resonance Imaging

Background: Scan acceleration techniques, such as parallel imaging, can reduce scan times, but reliability is essential to implement these techniques in neuroimaging. Objective: To evaluate the reproducibility of the longitudinal changes in brain morphology determined by longitudinal voxel-based morphometry (VBM) between non-accelerated and accelerated magnetic resonance images (MRI) in normal aging, mild cognitive impairment (MCI), and Alzheimer’s disease (AD). Methods: Using data from the Alzheimer’s Disease Neuroimaging Initiative (ADNI) 2 database, comprising subjects who underwent non-accelerated and accelerated structural T1-weighted MRI at screening and at a 2-year follow-up on 3.0 T Philips scanners, we examined the reproducibility of longitudinal gray matter volume changes determined by longitudinal VBM processing between non-accelerated and accelerated imaging in 50 healthy elderly subjects, 54 MCI patients, and eight AD patients. Results: The intraclass correlation coefficient (ICC) maps differed among the three groups. The mean ICC was 0.72 overall (healthy elderly, 0.63; MCI, 0.75; AD, 0.63), and the ICC was good to excellent (0.6–1.0) for 81.4%of voxels (healthy elderly, 64.8%; MCI, 85.0%; AD, 65.0%). The differences in image quality (head motion) were not significant (Kruskal–Wallis test, p = 0.18) and the within-subject standard deviations of longitudinal gray matter volume changes were similar among the groups. Conclusion: The results indicate that the reproducibility of longitudinal gray matter volume changes determined by VBM between non-accelerated and accelerated MRI is good to excellent for many regions but may vary between diseases and regions.

Bi-parametric magnetic resonance imaging based radiomics for the identification of benign and malignant prostate lesions: cross-vendor validation

The purpose of this study was to develop Bi-parametric Magnetic Resonance Imaging (BP-MRI) based radiomics models for differentiation between benign and malignant prostate lesions, and to cross-vendor validate the generalization ability of the models. The prebiopsy BP-MRI data (T2-Weighted Image [T2WI] and the Apparent Diffusion Coefficient [ADC]) of 459 patients with clinical suspicion of prostate cancer were acquired using two scanners from different vendors. The prostate biopsies are the reference standard for diagnosing benign and malignant prostate lesions. The training set was 168 patients' data from Siemens (Vendor 1), and the inner test set was 70 patients' data from the same vendor. The external test set was 221 patients' data from GE (Vendor 2). The lesion Region of Interest (ROI) was manually delineated by experienced radiologists. A total of 851 radiomics features including shape, first-order statistical, texture, and wavelet features were extracted from ROI in T2WI and ADC, respectively. Two feature-ranking methods (Minimum Redundancy Maximum Relevance [MRMR] and Wilcoxon Rank-Sum Test [WRST]) and three classifiers (Random Forest [RF], Support Vector Machine [SVM], and the Least Absolute Shrinkage and Selection Operator [LASSO] regression) were investigated for their efficacy in building single-parametric radiomics signatures. A biparametric radiomics model was built by combining the optimal single-parametric radiomics signatures. A comprehensive diagnosis model was built by combining the biparametric radiomics model with age and Prostate Specific Antigen (PSA) value using multivariable logistic regression. All models were built in the training set and independently validated in the inner and external test sets, and the performances of models in the diagnosis of benign and malignant prostate lesions were quantified by the Area Under the Receiver Operating Characteristic Curve (AUC). The mean AUCs of the inner and external test sets were calculated for each model. The non-inferiority test was used to test if the AUC of model in external test was not inferior to the AUC of model in inner test. Combining MRMR and LASSO produced the best-performing single-parametric radiomics signatures with the highest mean AUC of 0.673 for T2WI (inner test AUC = 0.729 vs. external test AUC = 0.616, p = 0.569) and the highest mean AUC of 0.810 for ADC (inner test AUC = 0.822 vs. external test AUC = 0.797, p = 0.102). The biparametric radiomics model produced a mean AUC of 0.833 (inner test AUC = 0.867 vs. external test AUC = 0.798, p = 0.051). The comprehensive diagnosis model had an improved mean AUC of 0.911 (inner test AUC = 0.935 vs. external test AUC = 0.886, p = 0.010). The comprehensive diagnosis model for differentiating benign from malignant prostate lesions was accurate and generalizable.

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

Importance

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

Objective

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

Design, Setting, and Participants

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

Exposures

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

Main Outcomes and Measures

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

Results

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

Conclusions and Relevance

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

The role of diffusion tensor imaging in idiopathic normal pressure hydrocephalus: A literature review

Idiopathic normal pressure hydrocephalus (iNPH) is a syndrome that comprises a triad of gait disturbance, dementia and urinary incontinence, associated with ventriculomegaly in the absence of elevated intraventricular cerebrospinal fluid (CSF) pressure. It is important to identify patients with iNPH because some of its clinical features may be reversed by the insertion of a CSF shunt. The diagnosis is based on clinical history, physical examination and brain imaging, especially magnetic resonance imaging (MRI). Recently, some papers have investigated the role of diffusion tensor imaging (DTI) in evaluating white matter alterations in patients with iNPH. DTI analysis in specific anatomical regions seems to be a promising MR biomarker of iNPH and could also be used in the differential diagnosis from other dementias. However, there is a substantial lack of structured reviews on this topic. Thus, we performed a literature search and analyzed the most recent and pivotal articles that investigated the role of DTI in iNPH in order to provide an up-to-date overview of the application of DTI in this setting. We reviewed studies published between January 2000 and June 2020. Thirty-eight studies and four reviews were included. Despite heterogeneity in analysis approaches, the majority of studies reported significant correlations between DTI and clinical symptoms in iNPH patients, as well as different DTI patterns in patients with iNPH compared to those with Alzheimer or Parkinson diseases. It remains to be determined whether DTI could predict the success after CSF shunting.

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.

Longitudinal Reproducibility of Neurite Orientation Dispersion and Density Imaging (NODDI) Derived Metrics in the White Matter

Diffusion-weighted magnetic resonance imaging (DWI) is undergoing constant evolution with the ambitious goal of developing in-vivo histology of the brain. A recent methodological advancement is Neurite Orientation Dispersion and Density Imaging (NODDI), a histologically validated multi-compartment model to yield microstructural features of brain tissue such as geometric complexity and neurite packing density, which are especially useful in imaging the white matter. Since NODDI is increasingly popular in clinical research and fields such as developmental neuroscience and neuroplasticity, it is of vast importance to characterize its reproducibility (or reliability). We acquired multi-shell DWI data in 29 healthy young subjects twice over a rescan interval of 4 weeks to assess the within-subject coefficient of variation (CV_WS), between-subject coefficient of variation (CV_BS) and the intraclass correlation coefficient (ICC), respectively. Using these metrics, we compared regional and voxel-by-voxel reproducibility of the most common image analysis approaches (tract-based spatial statistics [TBSS], voxel-based analysis with different extents of smoothing ["VBM-style"], ROI-based analysis). We observed high test-retest reproducibility for the orientation dispersion index (ODI) and slightly worse results for the neurite density index (NDI). Our findings also suggest that the choice of analysis approach might have significant consequences for the results of a study. Collectively, the voxel-based approach with Gaussian smoothing kernels of ≥4 mm FWHM and ROI-averaging yielded the highest reproducibility across NDI and ODI maps (CV_WS mostly ≤3%, ICC mostly ≥0.8), respectively, whilst smaller kernels and TBSS performed consistently worse. Furthermore, we demonstrate that image quality (signal-to-noise ratio [SNR]) is an important determinant of NODDI metric reproducibility. We discuss the implications of these results for longitudinal and cross-sectional research designs commonly employed in the neuroimaging field.

MICRA: Microstructural image compilation with repeated acquisitions

We provide a rich multi-contrast microstructural MRI dataset acquired on an ultra-strong gradient 3T Connectom MRI scanner comprising 5 repeated sets of MRI microstructural contrasts in 6 healthy human participants. The availability of data sets that support comprehensive simultaneous assessment of test-retest reliability of multiple microstructural contrasts (i.e., those derived from advanced diffusion, multi-component relaxometry and quantitative magnetisation transfer MRI) in the same population is extremely limited. This unique dataset is offered to the imaging community as a test-bed resource for conducting specialised analyses that may assist and inform their current and future research. The Microstructural Image Compilation with Repeated Acquisitions (MICRA) dataset includes raw data and computed microstructure maps derived from multi-shell and multi-direction encoded diffusion, multi-component relaxometry and quantitative magnetisation transfer acquisition protocols. Our data demonstrate high reproducibility of several microstructural MRI measures across scan sessions as shown by intra-class correlation coefficients and coefficients of variation. To illustrate a potential use of the MICRA dataset, we computed sample sizes required to provide sufficient statistical power a priori across different white matter pathways and microstructure measures for different statistical comparisons. We also demonstrate whole brain white matter voxel-wise repeatability in several microstructural maps. The MICRA dataset will be of benefit to researchers wishing to conduct similar reliability tests, power estimations or to evaluate the robustness of their own analysis pipelines.

Tractography Alterations in the Arcuate and Uncinate Fasciculi in Post-Stroke Aphasia

Fiber tractography based on diffuse tensor imaging (DTI) can reveal three-dimensional white matter connectivity of the human brain. Tractography is a non-invasive method of visualizing cerebral white matter structures in vivo, including neural pathways surrounding the ischemic area. DTI may be useful for elucidating alterations in brain connectivity resulting from neuroplasticity after stroke. We present a case of a male patient who developed significant mixed aphasia following ischemic stroke. The patient had been treated by mechanical thrombectomy followed by an early rehabilitation, in conjunction with transcranial direct current stimulation (tDCS). DTI was used to examine the arcuate fasciculus and uncinate fasciculus upon admission and again at three months post-stroke. Results showed an improvement in the patient’s symptoms of aphasia, which was associated with changes in the volume and numbers of tracts in the uncinate fasciculus and the arcuate fasciculus.

Cerebrovascular pathology in Alzheimer's disease: Hopes and gaps

Alzheimer's disease (AD) is recognized as multifactorial and heterogeneous disease with multiple contributors to its pathophysiology, including vascular dysfunction. Given that a revision of the AT(N) classification is expected in the near future, the present work supports the importance to add an additional vascular (V) category to the framework. In particular, we attempt to shed light on the vascular markers and risk factors that are currently ready-to-be-added to the framework: i) lacunes, ii) white matter hyperintensities and iii) microbleeds seen in Flair, T2* weighted imaging and susceptibility images (SWI). Next, we discuss the added value of other types of imaging, such as diffusion-based metrics and advanced perfusion sequences to encompass more subtle vascular dysfunction. Finally, we highlight the importance to add information about the following cardiovascular risk factors to the framework: history of hypertension, obesity, and diabetes. We believe that adding a V category to the AT(N) framework will improve AD classification and foster efforts to apply the right drug(s) at the right time in the right AD subgroups. Brief communication The present work supports the importance to add an additional vascular (V) category to the AT(N) framework and shed light on the vascular MRI markers and risk factors that are currently ready-to-be-added to the framework.