TOWARDS OPTIMAL LINEAR ESTIMATION OF ORIENTATION DISTRIBUTION FUNCTIONS WITH ARBITRARILY SAMPLED DIFFUSION MRI DATA

Diffusion MRI with Machine Learning

Diffusion-weighted magnetic resonance imaging (dMRI) of the brain offers unique capabilities including noninvasive probing of tissue microstructure and structural connectivity. It is widely used for clinical assessment of disease and injury, and for neuroscience research. Analyzing the dMRI data to extract useful information for medical and scientific purposes can be challenging. The dMRI measurements may suffer from strong noise and artifacts, and may exhibit high inter-session and inter-scanner variability in the data, as well as inter-subject heterogeneity in brain structure. Moreover, the relationship between measurements and the phenomena of interest can be highly complex. Recent years have witnessed increasing use of machine learning methods for dMRI analysis. This manuscript aims to assess these efforts, with a focus on methods that have addressed data preprocessing and harmonization, microstructure mapping, tractography, and white matter tract analysis. We study the main findings, strengths, and weaknesses of the existing methods and suggest topics for future research. We find that machine learning may be exceptionally suited to tackle some of the difficult tasks in dMRI analysis. However, for this to happen, several shortcomings of existing methods and critical unresolved issues need to be addressed. There is a pressing need to improve evaluation practices, to increase the availability of rich training datasets and validation benchmarks, as well as model generalizability, reliability, and explainability concerns.

BrainSuite BIDS App: Containerized Workflows for MRI Analysis

There has been a concerted effort by the neuroimaging community to establish standards for computational methods for data analysis that promote reproducibility and portability. In particular, the Brain Imaging Data Structure (BIDS) specifies a standard for storing imaging data, and the related BIDS App methodology provides a standard for implementing containerized processing environments that include all necessary dependencies to process BIDS datasets using image processing workflows. We present the BrainSuite BIDS App, which encapsulates the core MRI processing functionality of BrainSuite within the BIDS App framework. Specifically, the BrainSuite BIDS App implements a participant-level workflow comprising three pipelines and a corresponding set of group-level analysis workflows for processing the participant-level outputs. The Anatomical Pipeline extracts cortical surface models from a T1-weighted (T1w) MRI. It then performs surface-constrained volumetric registration to align the T1w MRI to a labeled anatomical atlas, which is used to delineate anatomical regions of interest in the MRI brain volume and on the cortical surface models. The Diffusion Pipeline processes diffusion-weighted imaging (DWI) data, with steps that include coregistering the DWI data to the T1w scan, correcting for susceptibility-induced geometric image distortion, and fitting diffusion models to the DWI data. The Functional Pipeline performs fMRI processing using a combination of FSL, AFNI, and BrainSuite tools. It coregisters the fMRI data to the T1w image, then transforms the data to the anatomical atlas space and to the Human Connectome Project's grayordinate space. The outputs of each pipeline can then be processed during group-level analysis. The outputs of the Anatomical Pipeline and the Diffusion Pipeline are analyzed using the BrainSuite Statistics Toolbox in R (bstr), which provides functionality for hypothesis testing and statistical modeling. The outputs of the Functional Pipeline can be analyzed using atlas-based or atlas-free statistical methods during group-level processing. These analyses include the application of BrainSync, which synchronizes the time-series data temporally and enables comparison of resting-state or task-based fMRI data across scans. We also present the BrainSuite Dashboard quality control system, which provides a browser-based interface for reviewing the outputs of individual modules of the participant-level pipelines across a study in real-time as they are generated. BrainSuite Dashboard facilitates rapid review of intermediate results, enabling users to identify processing errors and make adjustments to processing parameters if necessary. The comprehensive functionality included in the BrainSuite BIDS App provides a mechanism for rapidly deploying the BrainSuite workflows into new environments to perform large-scale studies. We demonstrate the capabilities of the BrainSuite BIDS App using structural, diffusion, and functional MRI data from the Amsterdam Open MRI Collection's Population Imaging of Psychology dataset.

Geometric deep learning for diffusion MRI signal reconstruction with continuous samplings (DISCUS)

Diffusion-weighted magnetic resonance imaging (dMRI) permits a detailed in-vivo analysis of neuroanatomical microstructure, invaluable for clinical and population studies. However, many measurements with different diffusion-encoding directions and possibly b-values are necessary to infer the underlying tissue microstructure within different imaging voxels accurately. Two challenges particularly limit the utility of dMRI: long acquisition times limit feasible scans to only a few directional measurements, and the heterogeneity of acquisition schemes across studies makes it difficult to combine datasets. Left unaddressed by previous learning-based methods that only accept dMRI data adhering to the specific acquisition scheme used for training, there is a need for methods that accept and predict signals for arbitrary diffusion encodings. Addressing these challenges, we describe the first geometric deep learning method for continuous dMRI signal reconstruction for arbitrary diffusion sampling schemes for both the input and output. Our method combines the reconstruction accuracy and robustness of previous learning-based methods with the flexibility of model-based methods, for example, spherical harmonics or SHORE. We demonstrate that our method outperforms model-based methods and performs on par with discrete learning-based methods on single-, multi-shell, and grid-based diffusion MRI datasets. Relevant for dMRI-derived analyses, we show that our reconstruction translates to higher-quality estimates of frequently used microstructure models compared to other reconstruction methods, enabling high-quality analyses even from very short dMRI acquisitions.

Chronic anemia: The effects on the connectivity of white matter

Chronic anemia is commonly observed in patients with hemoglobinopathies, mainly represented by disorders of altered hemoglobin (Hb) structure (sickle cell disease, SCD) and impaired Hb synthesis (e.g. thalassemia syndromes, non-SCD anemia). Both hemoglobinopathies have been associated with white matter (WM) alterations. Novel structural MRI research in our laboratory demonstrated that WM volume was diffusely lower in deep, watershed areas proportional to anemia severity. Furthermore, diffusion tensor imaging analysis has provided evidence that WM microstructure is disrupted proportionally to Hb level and oxygen saturation. SCD patients have been widely studied and demonstrate lower fractional anisotropy (FA) in the corticospinal tract and cerebellum across the internal capsule and corpus callosum. In the present study, we compared 19 SCD and 15 non-SCD anemia patients with a wide range of Hb values allowing the characterization of the effects of chronic anemia in isolation of sickle Hb. We performed a tensor analysis to quantify FA changes in WM connectivity in chronic anemic patients. We calculated the volumetric mean of FA along the pathway of tracks connecting two regions of interest defined by BrainSuite's BCI-DNI atlas. In general, we found lower FA values in anemic patients; indicating the loss of coherence in the main diffusion direction that potentially indicates WM injury. We saw a positive correlation between FA and hemoglobin in these same regions, suggesting that decreased WM microstructural integrity FA is highly driven by chronic hypoxia. The only connection that did not follow this pattern was the connectivity within the left middle-inferior temporal gyrus. Interestingly, more reductions in FA were observed in non-SCD patients (mainly along with intrahemispheric WM bundles and watershed areas) than the SCD patients (mainly interhemispheric).