Neurological soft signs in adolescents are associated with brain structure.

Neurological soft signs (NSS) are minor deviations in motor performance. During childhood and adolescence, NSS are examined for functional motor phenotyping to describe development, to screen for comorbidities, and to identify developmental vulnerabilities. Here, we investigate underlying brain structure alterations in association with NSS in physically trained adolescents. Male adolescent athletes (n = 136, 13-16 years) underwent a standardized neurological examination including 28 tests grouped into 6 functional clusters. Non-optimal performance in at least 1 cluster was rated as NSS (NSS+ group). Participants underwent T1- and diffusion-weighted magnetic resonance imaging. Cortical volume, thickness, and local gyrification were calculated using Freesurfer. Measures of white matter microstructure (Free-water (FW), FW-corrected fractional anisotropy (FAt), axial and radial diffusivity (ADt, RDt)) were calculated using tract-based spatial statistics. General linear models with age and handedness as covariates were applied to assess differences between NSS+ and NSS- group. We found higher gyrification in a large cluster spanning the left superior frontal and parietal areas, and widespread lower FAt and higher RDt compared with the NSS- group. This study shows that NSS in adolescents are associated with brain structure alterations. Underlying mechanisms may include alterations in synaptic pruning and axon myelination, which are hallmark processes of brain maturation.

Mean kurtosis-Curve (MK-Curve) correction improves the test-retest reproducibility of diffusion kurtosis imaging at 3 T.

Diffusion kurtosis imaging (DKI) is applied to gain insights into the microstructural organization of brain tissues. However, the reproducibility of DKI outside brain white matter, particularly in combination with advanced estimation to remedy its noise sensitivity, remains poorly characterized. Therefore, in this study, we investigated the variability and reliability of DKI metrics while correcting implausible values with a fit method called mean kurtosis (MK)-Curve. A total of 10 volunteers (four women; age: 41.4 ± 9.6 years) were included and underwent two MRI examinations of the brain. The images were acquired on a clinical 3-T scanner and included a T1-weighted image and a diffusion sequence with multiple diffusion weightings suitable for DKI. Region of interest analysis of common kurtosis and tensor metrics derived with the MK-Curve DKI fit was performed, including intraclass correlation (ICC) and Bland-Altman (BA) plot statistics. A p value of less than 0.05 was considered statistically significant. The analyses showed good to excellent agreement of both kurtosis tensor- and diffusion tensor-derived MK-Curve-corrected metrics (ICC values: 0.77-0.98 and 0.87-0.98, respectively), with the exception of two DKI-derived metrics (axial kurtosis in the cortex: ICC = 0.68, and radial kurtosis in deep gray matter: ICC = 0.544). Non-MK-Curve-corrected kurtosis tensor-derived metrics ranged from 0.01 to 0.52 and diffusion tensor-derived metrics from 0.06 to 0.66, indicating poor to moderate reliability. No structural bias was observed in the BA plots for any of the diffusion metrics. In conclusion, MK-Curve-corrected DKI metrics of the human brain can be reliably acquired in white and gray matter at 3 T and DKI metrics have good to excellent agreement in a test-retest setting.

Anatomically informed multi-level fiber tractography for targeted virtual dissection.

OBJECTIVES: Diffusion-weighted MRI can assist preoperative planning by reconstructing the trajectory of eloquent fiber pathways, such as the corticospinal tract (CST). However, accurate reconstruction of the full extent of the CST remains challenging with existing tractography methods. We suggest a novel tractography algorithm exploiting unused fiber orientations to produce more complete and reliable results. METHODS: Our novel approach, referred to as multi-level fiber tractography (MLFT), reconstructs fiber pathways by progressively considering previously unused fiber orientations at multiple levels of tract propagation. Anatomical priors are used to minimize the number of false-positive pathways. The MLFT method was evaluated on synthetic data and in vivo data by reconstructing the CST while compared to conventional tractography approaches. RESULTS: The radial extent of MLFT reconstructions is comparable to that of probabilistic reconstruction: [Formula: see text] for the left and [Formula: see text] for the right hemisphere according to Wilcoxon test, while achieving significantly higher topography preservation compared to probabilistic tractography: [Formula: see text]. DISCUSSION: MLFT provides a novel way to reconstruct fiber pathways by adding the capability of including branching pathways in fiber tractography. Thanks to its robustness, feasible reconstruction extent and topography preservation, our approach may assist in clinical practice as well as in virtual dissection studies.

Quantitative diffusion-weighted MRI response assessment in rhabdomyosarcoma: an international retrospective study on behalf of the European paediatric Soft tissue sarcoma Study Group Imaging Committee.

OBJECTIVE: To investigate the feasibility of diffusion-weighted magnetic resonance imaging (DW-MRI) as a predictive imaging marker after neoadjuvant chemotherapy in patients with rhabdomyosarcoma. MATERIAL AND METHODS: We performed a multicenter retrospective study including pediatric, adolescent and young adult patients with rhabdomyosarcoma, Intergroup Rhabdomyosarcoma Study group III/IV, treated according to the European paediatric Soft tissue sarcoma Study Group (EpSSG) RMS2005 or MTS2008 studies. DW-MRI was performed according to institutional protocols. We performed two-dimensional single-slice tumor delineation. Areas of necrosis or hemorrhage were delineated to be excluded in the primary analysis. Mean, median and 5th and 95th apparent diffusion coefficient (ADC) were extracted. RESULTS: Of 134 included patients, 82 had measurable tumor at diagnosis and response and DW-MRI scans of adequate quality and were included in the analysis. Technical heterogeneity in scan acquisition protocols and scanners was observed. Mean ADC at diagnosis was 1.1 (95% confidence interval [CI]: 1.1-1.2) (all ADC expressed in * 10-3 mm2/s), versus 1.6 (1.5-1.6) at response assessment. The 5th percentile ADC was 0.8 (0.7-0.9) at diagnosis and 1.1 (1.0-1.2) at response. Absolute change in mean ADC after neoadjuvant chemotherapy was 0.4 (0.3-0.5). Exploratory analyses for association between ADC and clinical parameters showed a significant difference in mean ADC at diagnosis for alveolar versus embryonal histology. Landmark analysis at nine weeks after the date of diagnosis showed no significant association (hazard ratio 1.3 [0.6-3.2]) between the mean ADC change and event-free survival. CONCLUSION: A significant change in the 5th percentile and the mean ADC after chemotherapy was observed. Strong heterogeneity was identified in DW-MRI acquisition protocols between centers and in individual patients.

Cross-site harmonization of multi-shell diffusion MRI measures based on rotational invariant spherical harmonics (RISH).

Quantification methods based on the acquisition of diffusion magnetic resonance imaging (dMRI) with multiple diffusion weightings (e.g., multi-shell) are becoming increasingly applied to study the in-vivo brain. Compared to single-shell data for diffusion tensor imaging (DTI), multi-shell data allows to apply more complex models such as diffusion kurtosis imaging (DKI), which attempts to capture both diffusion hindrance and restriction effects, or biophysical models such as NODDI, which attempt to increase specificity by separating biophysical components. Because of the strong dependence of the dMRI signal on the measurement hardware, DKI and NODDI metrics show scanner and site differences, much like other dMRI metrics. These effects limit the implementation of multi-shell approaches in multicenter studies, which are needed to collect large sample sizes for robust analyses. Recently, a post-processing technique based on rotation invariant spherical harmonics (RISH) features was introduced to mitigate cross-scanner differences in DTI metrics. Unlike statistical harmonization methods, which require repeated application to every dMRI metric of choice, RISH harmonization is applied once on the raw data, and can be followed by any analysis. RISH features harmonization has been tested on DTI features but not its generalizability to harmonize multi-shell dMRI. In this work, we investigated whether performing the RISH features harmonization of multi-shell dMRI data removes cross-site differences in DKI and NODDI metrics while retaining longitudinal effects. To this end, 46 subjects underwent a longitudinal (up to 3 time points) two-shell dMRI protocol at 3 imaging sites. DKI and NODDI metrics were derived before and after harmonization and compared both at the whole brain level and at the voxel level. Then, the harmonization effects on cross-sectional and on longitudinal group differences were evaluated. RISH features averaged for each of the 3 sites exhibited prominent between-site differences in the frontal and posterior part of the brain. Statistically significant differences in fractional anisotropy, mean diffusivity and mean kurtosis were observed both at the whole brain and voxel level between all the acquisition sites before harmonization, but not after. The RISH method also proved effective to harmonize NODDI metrics, particularly in white matter. The RISH based harmonization maintained the magnitude and variance of longitudinal changes as compared to the non-harmonized data of all considered metrics. In conclusion, the application of RISH feature based harmonization to multi-shell dMRI data can be used to remove cross-site differences in DKI metrics and NODDI analyses, while retaining inherent relations between longitudinal acquisitions.

Insights from the IronTract challenge: Optimal methods for mapping brain pathways from multi-shell diffusion MRI.

Limitations in the accuracy of brain pathways reconstructed by diffusion MRI (dMRI) tractography have received considerable attention. While the technical advances spearheaded by the Human Connectome Project (HCP) led to significant improvements in dMRI data quality, it remains unclear how these data should be analyzed to maximize tractography accuracy. Over a period of two years, we have engaged the dMRI community in the IronTract Challenge, which aims to answer this question by leveraging a unique dataset. Macaque brains that have received both tracer injections and ex vivo dMRI at high spatial and angular resolution allow a comprehensive, quantitative assessment of tractography accuracy on state-of-the-art dMRI acquisition schemes. We find that, when analysis methods are carefully optimized, the HCP scheme can achieve similar accuracy as a more time-consuming, Cartesian-grid scheme. Importantly, we show that simple pre- and post-processing strategies can improve the accuracy and robustness of many tractography methods. Finally, we find that fiber configurations that go beyond crossing (e.g., fanning, branching) are the most challenging for tractography. The IronTract Challenge remains open and we hope that it can serve as a valuable validation tool for both users and developers of dMRI analysis methods.

Association between Motor Planning and the Frontoparietal Network in Children: An Exploratory Multimodal Study.

OBJECTIVE: Evidence from adult literature shows the involvement of cortical grey matter areas of the frontoparietal lobe and the white matter bundle, the superior longitudinal fasciculus (SLF) in motor planning. This is yet to be confirmed in children. METHOD: A multimodal study was designed to probe the neurostructural basis of childhood motor planning. Behavioural (motor planning), magnetic resonance imaging (MRI) and diffusion weighted imaging (DWI) data were acquired from 19 boys aged 8-11 years. Motor planning was assessed using the one and two colour sequences of the octagon task. The MRI data were preprocessed and analysed using FreeSurfer 6.0. Cortical thickness and cortical surface area were extracted from the caudal middle frontal gyrus (MFG), superior frontal gyrus (SFG), precentral gyrus (PcG), supramarginal gyrus (SMG), superior parietal lobe (SPL) and the inferior parietal lobe (IPL) using the Desikan-Killiany atlas. The DWI data were preprocessed and analysed using ExploreDTI 4.8.6 and the white matter tract, the SLF was reconstructed. RESULTS: Motor planning of the two colour sequence was associated with cortical thickness of the bilateral MFG and left SFG, PcG, IPL and SPL. The right SLF was related to motor planning for the two colour sequence as well as with the left cortical thickness of the SFG. CONCLUSION: Altogether, morphology within frontodorsal circuity, and the white matter bundles that support communication between them, may be associated with individual differences in childhood motor planning.

On the generalizability of diffusion MRI signal representations across acquisition parameters, sequences and tissue types: Chronicles of the MEMENTO challenge.

Diffusion MRI (dMRI) has become an invaluable tool to assess the microstructural organization of brain tissue. Depending on the specific acquisition settings, the dMRI signal encodes specific properties of the underlying diffusion process. In the last two decades, several signal representations have been proposed to fit the dMRI signal and decode such properties. Most methods, however, are tested and developed on a limited amount of data, and their applicability to other acquisition schemes remains unknown. With this work, we aimed to shed light on the generalizability of existing dMRI signal representations to different diffusion encoding parameters and brain tissue types. To this end, we organized a community challenge - named MEMENTO, making available the same datasets for fair comparisons across algorithms and techniques. We considered two state-of-the-art diffusion datasets, including single-diffusion-encoding (SDE) spin-echo data from a human brain with over 3820 unique diffusion weightings (the MASSIVE dataset), and double (oscillating) diffusion encoding data (DDE/DODE) of a mouse brain including over 2520 unique data points. A subset of the data sampled in 5 different voxels was openly distributed, and the challenge participants were asked to predict the remaining part of the data. After one year, eight participant teams submitted a total of 80 signal fits. For each submission, we evaluated the mean squared error, the variance of the prediction error and the Bayesian information criteria. The received submissions predicted either multi-shell SDE data (37%) or DODE data (22%), followed by cartesian SDE data (19%) and DDE (18%). Most submissions predicted the signals measured with SDE remarkably well, with the exception of low and very strong diffusion weightings. The prediction of DDE and DODE data seemed more challenging, likely because none of the submissions explicitly accounted for diffusion time and frequency. Next to the choice of the model, decisions on fit procedure and hyperparameters play a major role in the prediction performance, highlighting the importance of optimizing and reporting such choices. This work is a community effort to highlight strength and limitations of the field at representing dMRI acquired with trending encoding schemes, gaining insights into how different models generalize to different tissue types and fiber configurations over a large range of diffusion encodings.

Tractography dissection variability: What happens when 42 groups dissect 14 white matter bundles on the same dataset?

White matter bundle segmentation using diffusion MRI fiber tractography has become the method of choice to identify white matter fiber pathways in vivo in human brains. However, like other analyses of complex data, there is considerable variability in segmentation protocols and techniques. This can result in different reconstructions of the same intended white matter pathways, which directly affects tractography results, quantification, and interpretation. In this study, we aim to evaluate and quantify the variability that arises from different protocols for bundle segmentation. Through an open call to users of fiber tractography, including anatomists, clinicians, and algorithm developers, 42 independent teams were given processed sets of human whole-brain streamlines and asked to segment 14 white matter fascicles on six subjects. In total, we received 57 different bundle segmentation protocols, which enabled detailed volume-based and streamline-based analyses of agreement and disagreement among protocols for each fiber pathway. Results show that even when given the exact same sets of underlying streamlines, the variability across protocols for bundle segmentation is greater than all other sources of variability in the virtual dissection process, including variability within protocols and variability across subjects. In order to foster the use of tractography bundle dissection in routine clinical settings, and as a fundamental analytical tool, future endeavors must aim to resolve and reduce this heterogeneity. Although external validation is needed to verify the anatomical accuracy of bundle dissections, reducing heterogeneity is a step towards reproducible research and may be achieved through the use of standard nomenclature and definitions of white matter bundles and well-chosen constraints and decisions in the dissection process.

White matter microstructure of the neural emotion regulation circuitry in mild traumatic brain injury.

Emotion regulation is related to recovery after mild traumatic brain injury (mTBI). This longitudinal tractography study examined white matter tracts subserving emotion regulation across the spectrum of mTBI, with a focus on persistent symptoms. Four groups were examined: (a) symptomatic (n = 33) and (b) asymptomatic (n = 20) patients with uncomplicated mTBI (i.e., no lesions on computed tomography [CT]), (c) patients with CT-lesions in the frontal areas (n = 14), and (d) healthy controls (HC) (n = 20). Diffusion and conventional MRI were performed approximately 1- and 3-months post-injury. Whole-brain deterministic tractography followed by region of interest analyses was used to identify forceps minor (FM), uncinate fasciculus (UF), and cingulum bundle as tracts of interest. An adjusted version of the ExploreDTI Atlas Based Tractography method was used to obtain reliable tracts for every subject. Mean fractional anisotropy (FA), mean, radial and axial diffusivity (MD, RD, AD), and number of streamlines were studied per tract. Linear mixed models showed lower FA, and higher MD, and RD of the right UF in asymptomatic patients with uncomplicated mTBI relative to symptomatic patients and HC. Diffusion alterations were most pronounced in the group with frontal lesions on CT, particularly in the FM and UF; these effects increased over time. Within the group of patients with uncomplicated mTBI, there were no associations of diffusion measures with the number of symptoms nor with lesions on conventional MRI. In conclusion, mTBI can cause microstructural changes in emotion regulation tracts, however, no explanation was found for the presence of symptoms.

The effect of gradient nonlinearities on fiber orientation estimates from spherical deconvolution of diffusion magnetic resonance imaging data.

Gradient nonlinearities in magnetic resonance imaging (MRI) cause spatially varying mismatches between the imposed and the effective gradients and can cause significant biases in rotationally invariant diffusion MRI measures derived from, for example, diffusion tensor imaging. The estimation of the orientational organization of fibrous tissue, which is nowadays frequently performed with spherical deconvolution techniques ideally using higher diffusion weightings, can likewise be biased by gradient nonlinearities. We explore the sensitivity of two established spherical deconvolution approaches to gradient nonlinearities, namely constrained spherical deconvolution (CSD) and damped Richardson-Lucy (dRL). Additionally, we propose an extension of dRL to take into account gradient imperfections, without the need of data interpolation. Simulations show that using the effective b-matrix can improve dRL fiber orientation estimation and reduces angular deviations, while CSD can be more robust to gradient nonlinearity depending on the implementation. Angular errors depend on a complex interplay of many factors, including the direction and magnitude of gradient deviations, underlying microstructure, SNR, anisotropy of the effective response function, and diffusion weighting. Notably, angular deviations can also be observed at lower b-values in contrast to the perhaps common assumption that only high b-value data are affected. In in vivo Human Connectome Project data and acquisitions from an ultrastrong gradient (300 mT/m) scanner, angular differences are observed between applying and not applying the effective gradients in dRL estimation. As even small angular differences can lead to error propagation during tractography and as such impact connectivity analyses, incorporating gradient deviations into the estimation of fiber orientations should make such analyses more reliable.

Spherical deconvolution with tissue-specific response functions and multi-shell diffusion MRI to estimate multiple fiber orientation distributions (mFODs).

In diffusion MRI, spherical deconvolution approaches can estimate local white matter (WM) fiber orientation distributions (FOD) which can be used to produce fiber tractography reconstructions. The applicability of spherical deconvolution to gray matter (GM), however, is still limited, despite its critical role as start/endpoint of WM fiber pathways. The advent of multi-shell diffusion MRI data offers additional contrast to model the GM signal but, to date, only isotropic models have been applied to GM. Evidence from both histology and high-resolution diffusion MRI studies suggests a marked anisotropic character of the diffusion process in GM, which could be exploited to improve the description of the cortical organization. In this study, we investigated whether performing spherical deconvolution with tissue specific models of both WM and GM can improve the characterization of the latter while retaining state-of-the-art performances in WM. To this end, we developed a framework able to simultaneously accommodate multiple anisotropic response functions to estimate multiple, tissue-specific, fiber orientation distributions (mFODs). As proof of principle, we used the diffusion kurtosis imaging model to represent the WM signal, and the neurite orientation dispersion and density imaging (NODDI) model to represent the GM signal. The feasibility of the proposed approach is shown with numerical simulations and with data from the Human Connectome Project (HCP). The performance of our method is compared to the current state of the art, multi-shell constrained spherical deconvolution (MSCSD). The simulations show that with our new method we can accurately estimate a mixture of two FODs at SNR≥50. With HCP data, the proposed method was able to reconstruct both tangentially and radially oriented FODs in GM, and performed comparably well to MSCSD in computing FODs in WM. When performing fiber tractography, the trajectories reconstructed with mFODs reached the cortex with more spatial continuity and for a longer distance as compared to MSCSD and allowed to reconstruct short trajectories tangential to the cortical folding. In conclusion, we demonstrated that our proposed method allows to perform spherical deconvolution of multiple anisotropic response functions, specifically improving the performances of spherical deconvolution in GM tissue.

The adverse effect of gradient nonlinearities on diffusion MRI: From voxels to group studies.

Nonlinearities of gradient magnetic fields in diffusion MRI (dMRI) can introduce systematic errors in estimates of diffusion measures. While there are correction methods that can compensate for these errors, as presented in the Human Connectome Project, such nonlinear effects are often assumed to be negligible for typical applications, and as a result, gradient nonlinearities are mostly left uncorrected. In this work, we perform a systematic analysis to investigate the effect of gradient nonlinearities on dMRI studies, from voxel-wise estimates to group study outcomes. We present a novel framework to quantify and visualize these effects by decomposing them into their magnitude and angle components. Mean magnitude deviation and fractional gradient anisotropy are introduced to quantify the distortions in the size and shape of gradient vector distributions. By means of Monte-Carlo simulations and real data from the Human Connectome Project, the errors on dMRI measures derived from the diffusion tensor imaging and diffusional kurtosis imaging are highlighted. We perform a group study to showcase the alteration in the significance and effect size due to ignoring gradient nonlinearity correction. Our results indicate that the effect of gradient field nonlinearities on dMRI studies can be significant and may complicate the interpretation of the results and conclusions.

Generalized Richardson-Lucy (GRL) for analyzing multi-shell diffusion MRI data.

Spherical deconvolution is a widely used approach to quantify the fiber orientation distribution (FOD) from diffusion MRI data of the brain. The damped Richardson-Lucy (dRL) is an algorithm developed to perform robust spherical deconvolution on single-shell diffusion MRI data while suppressing spurious FOD peaks due to noise or partial volume effects. Due to recent progress in acquisition hardware and scanning protocols, it is becoming increasingly common to acquire multi-shell diffusion MRI data, which allows for the modelling of multiple tissue types beyond white matter. While the dRL algorithm could, in theory, be directly applied to multi-shell data, it is not optimised to exploit its information content to model the signal from multiple tissue types. In this work, we introduce a new framework based on dRL - dubbed generalized Richardson-Lucy (GRL) - that uses multi-shell data in combination with user-chosen tissue models to disentangle partial volume effects and increase the accuracy in FOD estimation. Further, GRL estimates signal fraction maps associated to each user-selected model, which can be used during fiber tractography to dissect and terminate the tracking without need for additional structural data. The optimal weighting of multi-shell data in the fit and the robustness to noise and to partial volume effects of GRL was studied with synthetic data. Subsequently, we investigated the performance of GRL in comparison to dRL and to multi-shell constrained spherical deconvolution (MSCSD) on a high-resolution diffusion MRI dataset from the Human Connectome Project and on an MRI dataset acquired at 3T on a clinical scanner. In line with previous studies, we described the signal of the cerebrospinal-fluid and of the grey matter with isotropic diffusion models, whereas four diffusion models were considered to describe the white matter. With a third dataset including small diffusion weightings, we studied the feasibility of including intra-voxel incoherent motion effects due to blood pseudo-diffusion in the modelling. Further, the reliability of GRL was demonstrated with a test-retest scan of a subject acquired at 3T. Results of simulations show that GRL can robustly disentangle different tissue types at SNR above 20 with respect to the non-weighted image, and that it improves the angular accuracy of the FOD estimation as compared to dRL. On real data, GRL provides signal fraction maps that are physiologically plausible and consistent with those obtained with MSCSD, with correlation coefficients between the two methods up to 0.96. When considering IVIM effects, a high blood pseudo-diffusion fraction is observed in the medial temporal lobe and in the sagittal sinus. In comparison to dRL and MSCSD, GRL provided sharper FODs and less spurious peaks in presence of partial volume effects, but the FOD reconstructions are also highly dependent on the chosen modelling of white matter. When performing fiber tractography, GRL allows to terminate fiber tractography using the signal fraction maps, which results in a better tract termination at the grey-white matter interface or at the outer cortical surface. In terms of inter-scan reliability, GRL performed similarly to or better than compared methods. In conclusion, GRL offers a new modular and flexible framework to perform spherical deconvolution of multi-shell data.

Harmonization of diffusion MRI data sets with adaptive dictionary learning.

Diffusion magnetic resonance imaging can indirectly infer the microstructure of tissues and provide metrics subject to normal variability in a population. Potentially abnormal values may yield essential information to support analysis of controls and patients cohorts, but subtle confounds could be mistaken for purely biologically driven variations amongst subjects. In this work, we propose a new harmonization algorithm based on adaptive dictionary learning to mitigate the unwanted variability caused by different scanner hardware while preserving the natural biological variability of the data. Our harmonization algorithm does not require paired training data sets, nor spatial registration or matching spatial resolution. Overcomplete dictionaries are learned iteratively from all data sets at the same time with an adaptive regularization criterion, removing variability attributable to the scanners in the process. The obtained mapping is applied directly in the native space of each subject toward a scanner-space. The method is evaluated with a public database which consists of two different protocols acquired on three different scanners. Results show that the effect size of the four studied diffusion metrics is preserved while removing variability attributable to the scanner. Experiments with alterations using a free water compartment, which is not simulated in the training data, shows that the modifications applied to the diffusion weighted images are preserved in the diffusion metrics after harmonization, while still reducing global variability at the same time. The algorithm could help multicenter studies pooling their data by removing scanner specific confounds, and increase statistical power in the process.

Tractography reproducibility challenge with empirical data (TraCED): The 2017 ISMRM diffusion study group challenge.

BACKGROUND: Fiber tracking with diffusion-weighted MRI has become an essential tool for estimating in vivo brain white matter architecture. Fiber tracking results are sensitive to the choice of processing method and tracking criteria. PURPOSE: To assess the variability for an algorithm in group studies reproducibility is of critical context. However, reproducibility does not assess the validity of the brain connections. Phantom studies provide concrete quantitative comparisons of methods relative to absolute ground truths, yet do no capture variabilities because of in vivo physiological factors. The ISMRM 2017 TraCED challenge was created to fulfill the gap. STUDY TYPE: A systematic review of algorithms and tract reproducibility studies. SUBJECTS: Single healthy volunteers. FIELD STRENGTH/SEQUENCE: 3.0T, two different scanners by the same manufacturer. The multishell acquisition included b-values of 1000, 2000, and 3000 s/mm2 with 20, 45, and 64 diffusion gradient directions per shell, respectively. ASSESSMENT: Nine international groups submitted 46 tractography algorithm entries each consisting 16 tracts per scan. The algorithms were assessed using intraclass correlation (ICC) and the Dice similarity measure. STATISTICAL TESTS: Containment analysis was performed to assess if the submitted algorithms had containment within tracts of larger volume submissions. This also serves the purpose to detect if spurious submissions had been made. RESULTS: The top five submissions had high ICC and Dice >0.88. Reproducibility was high within the top five submissions when assessed across sessions or across scanners: 0.87-0.97. Containment analysis shows that the top five submissions are contained within larger volume submissions. From the total of 16 tracts as an outcome relatively the number of tracts with high, moderate, and low reproducibility were 8, 4, and 4. DATA CONCLUSION: The different methods clearly result in fundamentally different tract structures at the more conservative specificity choices. Data and challenge infrastructure remain available for continued analysis and provide a platform for comparison. LEVEL OF EVIDENCE: 5 Technical Efficacy Stage: 1 J. Magn. Reson. Imaging 2020;51:234-249.

Early prediction of unilateral cerebral palsy in infants at risk: MRI versus the hand assessment for infants.

BACKGROUND: Neonates with unilateral perinatal brain injury (UPBI) are at risk for developing unilateral spastic cerebral palsy (USCP). This study compares several predictors for USCP later in life. METHODS: Twenty-one preterm and 24 term born infants with UPBI were included, with an MRI scan including diffusion tensor imaging (DTI) performed at term equivalent age or around 3 months after birth, respectively. T2-weighted images and DTI-based tractography were used to measure the surface area, diameter, and fractional anisotropy (FA) of both corticospinal tracts (CSTs). The hand assessment for infants (HAI) was performed before 5, between 5 and 8 and between 8 and 12 months of (corrected) age. Asymmetry indices were derived from all techniques and related to USCP at ≥2 years of age. RESULTS: MRI measures and HAI scores were significantly lower for the affected compared to the unaffected side. Before 5 months of age, FA asymmetry on DTI yielded the highest area under the curve compared to conventional MRI and HAI. CONCLUSIONS: Prediction of USCP after UPBI is reliable using asymmetry of the CST on MRI, as well as clinical hand assessment. Before 5 months of age, DTI tractography provides strongest predictive information, while HAI specifically aids to prognosis of USCP at later age points.

Indirect frontocingulate structural connectivity predicts clinical response to accelerated rTMS in major depressive disorder

Background: Repetitive transcranial magnetic stimulation (rTMS) is an established treatment for major depressive disorder (MDD), but its clinical efficacy remains rather modest. One reason for this could be that the propagation of rTMS effects via structural connections from the stimulated area to deeper brain structures (such as the cingulate cortices) is suboptimal. Methods: We investigated whether structural connectivity ­ derived from diffusion MRI data ­ could serve as a biomarker to predict treatment response. We hypothesized that stronger structural connections between the patient-specific stimulation position in the left dorsolateral prefrontal cortex (dlPFC) and the cingulate cortices would predict better clinical outcomes. We applied accelerated intermittent theta burst stimulation (aiTBS) to the left dlPFC in 40 patients with MDD. We correlated baseline structural connectivity, quantified using various metrics (fractional anisotropy, mean diffusivity, tract density, tract volume and number of tracts), with changes in depression severity scores after aiTBS. Results: Exploratory results (p < 0.05) showed that structural connectivity between the patient-specific stimulation site and the caudal and posterior parts of the cingulate cortex had predictive potential for clinical response to aiTBS. Limitations: We used the diffusion tensor to perform tractography. A main limitation was that multiple fibre directions within voxels could not be resolved, which might have led to missing connections in some patients. Conclusion: Stronger structural frontocingular connections may be of essence to optimally benefit from left dlPFC rTMS treatment in MDD. Even though the results are promising, further investigation with larger numbers of patients, more advanced tractography algorithms and classic daily rTMS treatment paradigms is warranted. Clinical trial registration: http://clinicaltrials.gov/show/NCT01832805

Reducing variability in along-tract analysis with diffusion profile realignment.

Diffusion weighted magnetic resonance imaging (dMRI) provides a non invasive virtual reconstruction of the brain's white matter structures through tractography. Analyzing dMRI measures along the trajectory of white matter bundles can provide a more specific investigation than considering a region of interest or tract-averaged measurements. However, performing group analyses with this along-tract strategy requires correspondence between points of tract pathways across subjects. This is usually achieved by creating a new common space where the representative streamlines from every subject are resampled to the same number of points. If the underlying anatomy of some subjects was altered due to, e.g., disease or developmental changes, such information might be lost by resampling to a fixed number of points. In this work, we propose to address the issue of possible misalignment, which might be present even after resampling, by realigning the representative streamline of each subject in this 1D space with a new method, coined diffusion profile realignment (DPR). Experiments on synthetic datasets show that DPR reduces the coefficient of variation for the mean diffusivity, fractional anisotropy and apparent fiber density when compared to the unaligned case. Using 100 in vivo datasets from the human connectome project, we simulated changes in mean diffusivity, fractional anisotropy and apparent fiber density. Independent Student's t-tests between these altered subjects and the original subjects indicate that regional changes are identified after realignment with the DPR algorithm, while preserving differences previously detected in the unaligned case. This new correction strategy contributes to revealing effects of interest which might be hidden by misalignment and has the potential to improve the specificity in longitudinal population studies beyond the traditional region of interest based analysis and along-tract analysis workflows.

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

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