Generalized diffusion spectrum magnetic resonance imaging (GDSI) for model-free reconstruction of the ensemble average propagator

Recent Progress in Brain Network Models for Medical Applications: A Review

Importance: Pathological perturbations of the brain often spread via connectome to fundamentally alter functional consequences. By integrating multimodal neuroimaging data with mathematical neural mass modeling, brain network models (BNMs) enable to quantitatively characterize aberrant network dynamics underlying multiple neurological and psychiatric disorders. We delved into the advancements of BNM-based medical applications, discussed the prevalent challenges within this field, and provided possible solutions and future directions. Highlights: This paper reviewed the theoretical foundations and current medical applications of computational BNMs. Composed of neural mass models, the BNM framework allows to investigate large-scale brain dynamics behind brain diseases by linking the simulated functional signals to the empirical neurophysiological data, and has shown promise in exploring neuropathological mechanisms, elucidating therapeutic effects, and predicting disease outcome. Despite that several limitations existed, one promising trend of this research field is to precisely guide clinical neuromodulation treatment based on individual BNM simulation. Conclusion: BNM carries the potential to help understand the mechanism underlying how neuropathology affects brain network dynamics, further contributing to decision-making in clinical diagnosis and treatment. Several constraints must be addressed and surmounted to pave the way for its utilization in the clinic.

Q-space imaging based on Gaussian radial basis function with Laplace regularization

PURPOSE

To introduce the diffusion signal characteristics presented by spherical harmonics (SH) basis into the q-space imaging method based on Gaussian radial basis function (GRBF) to robustly reconstruct ensemble average diffusion propagator (EAP) in diffusion MRI (dMRI).

METHODS

We introduced the Laplacian regularization of the signal into the dMRI imaging method based on GRBF, and derived the relevant indicators of microstructure imaging and the orientation distribution function (ODF) providing fiber bundle direction information based on EAP. In addition, this method is combined with a multi-compartment model to calculate the diameter of fiber bundle axons. The evaluation of the results included qualitative comparisons and quantitative assessments of the signal fitting.

RESULTS

The results show that the proposed method achieves the more significant accuracy improvement in reconstructing signal. Meanwhile, ODFs estimated by the proposed method show the sharper profiles and less spurious peaks, even under the sparse and noisy conditions. In the 36 sets of axon diameter estimation experiments, 34 and 30 sets of results showed that the proposed method reduced the mean and SD of axon diameter estimates, respectively. Moreover, compared with the current state-of-the-art method, the mean and SD of axon diameter estimated by the proposed method are mostly lower, with 32 and 29 of 36 groups.

CONCLUSION

The proposed method outperforms the GRBF regarding signal fitting and the estimation of the EAP and ODF with multi-shell sparse samples. Moreover, it shows the potential to recover important features of microstructures with less uncertainty by using proposed method together with multi-compartment models.

Research Progress in Diffusion Spectrum Imaging

Studies have demonstrated that many regions in the human brain include multidirectional fiber tracts, in which the diffusion of water molecules within image voxels does not follow a Gaussian distribution. Therefore, the conventional diffusion tensor imaging (DTI) that hypothesizes a single fiber orientation within a voxel is intrinsically incapable of revealing the complex microstructures of brain tissues. Diffusion spectrum imaging (DSI) employs a pulse sequence with different b-values along multiple gradient directions to sample the diffusion information of water molecules in the entire q-space and then quantitatively estimates the diffusion profile using a probability density function with a high angular resolution. Studies have suggested that DSI can reliably observe the multidirectional fibers within each voxel and allow fiber tracking along different directions, which can improve fiber reconstruction reflecting the true but complicated brain structures that were not observed in the previous DTI studies. Moreover, with increasing angular resolution, DSI is able to reveal new neuroimaging biomarkers used for disease diagnosis and the prediction of disorder progression. However, so far, this method has not been used widely in clinical studies, due to its overly long scanning time and difficult post-processing. Within this context, the current paper aims to conduct a comprehensive review of DSI research, including the fundamental principles, methodology, and application progress of DSI tractography. By summarizing the DSI studies in recent years, we propose potential solutions towards the existing problem in the methodology and applications of DSI technology as follows: (1) using compressed sensing to undersample data and to reconstruct the diffusion signal may be an efficient and promising method for reducing scanning time; (2) the probability density function includes more information than the orientation distribution function, and it should be extended in application studies; and (3) large-sample study is encouraged to confirm the reliability and reproducibility of findings in clinical diseases. These findings may help deepen the understanding of the DSI method and promote its development in clinical applications.

Diffusion MRI data analysis assisted by deep learning synthesized anatomical images (DeepAnat)

Diffusion MRI is a useful neuroimaging tool for non-invasive mapping of human brain microstructure and structural connections. The analysis of diffusion MRI data often requires brain segmentation, including volumetric segmentation and cerebral cortical surfaces, from additional high-resolution T1-weighted (T1w) anatomical MRI data, which may be unacquired, corrupted by subject motion or hardware failure, or cannot be accurately co-registered to the diffusion data that are not corrected for susceptibility-induced geometric distortion. To address these challenges, this study proposes to synthesize high-quality T1w anatomical images directly from diffusion data using convolutional neural networks (CNNs) (entitled "DeepAnat"), including a U-Net and a hybrid generative adversarial network (GAN), and perform brain segmentation on synthesized T1w images or assist the co-registration using synthesized T1w images. The quantitative and systematic evaluations using data of 60 young subjects provided by the Human Connectome Project (HCP) show that the synthesized T1w images and results for brain segmentation and comprehensive diffusion analysis tasks are highly similar to those from native T1w data. The brain segmentation accuracy is slightly higher for the U-Net than the GAN. The efficacy of DeepAnat is further validated on a larger dataset of 300 more elderly subjects provided by the UK Biobank. Moreover, the U-Nets trained and validated on the HCP and UK Biobank data are shown to be highly generalizable to the diffusion data from Massachusetts General Hospital Connectome Diffusion Microstructure Dataset (MGH CDMD) acquired with different hardware systems and imaging protocols and therefore can be used directly without retraining or with fine-tuning for further improved performance. Finally, it is quantitatively demonstrated that the alignment between native T1w images and diffusion images uncorrected for geometric distortion assisted by synthesized T1w images substantially improves upon that by directly co-registering the diffusion and T1w images using the data of 20 subjects from MGH CDMD. In summary, our study demonstrates the benefits and practical feasibility of DeepAnat for assisting various diffusion MRI data analyses and supports its use in neuroscientific applications.

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

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

Population-Average Brain Templates and Application to Automated Voxel-Wise Analysis Pipelines for Cynomolgus Macaque

The growing number of neuroimaging studies of cynomolgus macaques require extending existing templates to facilitate species-specific application of voxel-wise neuroimaging methodologies. This study aimed to create population-averaged structural magnetic resonance imaging (MRI) and diffusion tensor imaging (DTI) templates for the cynomolgus macaques and apply the templates in fully automated voxel-wise analyses. We presented the development of symmetric and asymmetric MRI and DTI templates from a sample of 63 young male cynomolgus monkeys with the use of optimized template creation approaches. We also generated the associated average tissue probability maps and Diffeomorphic Anatomical Registration using Exponentiated Lie Algebra templates for use with the Statistical Parametric Mapping (SPM), as well as the average fractional anisotropy/skeleton targets for incorporation into tract-based spatial statistics (TBSS) framework. Both asymmetric and symmetric templates in a standardized coordinate space demonstrated low bias and high contrast. Fully automated processing using SPM was accomplished for all native MRI datasets and demonstrated outstanding performance regarding skull-stripping, segmentation, and normalization when using the MRI templates. Automated normalization to the DTI template was excellently achieved for all native DTI images using the TBSS pipeline. The cynomolgus MRI and DTI templates are anticipated to provide a common platform for precise single-subject data analysis and facilitate comparison of neuroimaging findings in cynomolgus monkeys across studies and sites. It is also hoped that the procedures of template creation and fully-automated voxel-wise frameworks will provide a straightforward avenue for investigating brain function, development, and neuro-psychopathological disorders in non-human primate models.

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

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

Comprehensive diffusion MRI dataset for in vivo human brain microstructure mapping using 300 mT/m gradients

Strong gradient systems can improve the signal-to-noise ratio of diffusion MRI measurements and enable a wider range of acquisition parameters that are beneficial for microstructural imaging. We present a comprehensive diffusion MRI dataset of 26 healthy participants acquired on the MGH-USC 3 T Connectome scanner equipped with 300 mT/m maximum gradient strength and a custom-built 64-channel head coil. For each participant, the one-hour long acquisition systematically sampled the accessible diffusion measurement space, including two diffusion times (19 and 49 ms), eight gradient strengths linearly spaced between 30 mT/m and 290 mT/m for each diffusion time, and 32 or 64 uniformly distributed directions. The diffusion MRI data were preprocessed to correct for gradient nonlinearity, eddy currents, and susceptibility induced distortions. In addition, scan/rescan data from a subset of seven individuals were also acquired and provided. The MGH Connectome Diffusion Microstructure Dataset (CDMD) may serve as a test bed for the development of new data analysis methods, such as fiber orientation estimation, tractography and microstructural modelling.

DeepDTI: High-fidelity six-direction diffusion tensor imaging using deep learning

Diffusion tensor magnetic resonance imaging (DTI) is unsurpassed in its ability to map tissue microstructure and structural connectivity in the living human brain. Nonetheless, the angular sampling requirement for DTI leads to long scan times and poses a critical barrier to performing high-quality DTI in routine clinical practice and large-scale research studies. In this work we present a new processing framework for DTI entitled DeepDTI that minimizes the data requirement of DTI to six diffusion-weighted images (DWIs) required by conventional voxel-wise fitting methods for deriving the six unique unknowns in a diffusion tensor using data-driven supervised deep learning. DeepDTI maps the input non-diffusion-weighted (b ​= ​0) image and six DWI volumes sampled along optimized diffusion-encoding directions, along with T1-weighted and T2-weighted image volumes, to the residuals between the input and high-quality output b = 0 image and DWI volumes using a 10-layer three-dimensional convolutional neural network (CNN). The inputs and outputs of DeepDTI are uniquely formulated, which not only enables residual learning to boost CNN performance but also enables tensor fitting of resultant high-quality DWIs to generate orientational DTI metrics for tractography. The very deep CNN used by DeepDTI leverages the redundancy in local and non-local spatial information and across diffusion-encoding directions and image contrasts in the data. The performance of DeepDTI was systematically quantified in terms of the quality of the output images, DTI metrics, DTI-based tractography and tract-specific analysis results. We demonstrate rotationally-invariant and robust estimation of DTI metrics from DeepDTI that are comparable to those obtained with two b ​= ​0 images and 21 DWIs for the primary eigenvector derived from DTI and two b ​= ​0 images and 26-30 DWIs for various scalar metrics derived from DTI, achieving 3.3-4.6 × ​acceleration, and twice as good as those of a state-of-the-art denoising algorithm at the group level. The twenty major white-matter tracts can be accurately identified from the tractography of DeepDTI results. The mean distance between the core of the major white-matter tracts identified from DeepDTI results and those from the ground-truth results using 18 ​b ​= ​0 images and 90 DWIs measures around 1-1.5 ​mm. DeepDTI leverages domain knowledge of diffusion MRI physics and power of deep learning to render DTI, DTI-based tractography, major white-matter tracts identification and tract-specific analysis more feasible for a wider range of neuroscientific and clinical studies.

In Vivo Imaging of Gray Matter Microstructure in Major Psychiatric Disorders: Opportunities for Clinical Translation

Postmortem studies reveal that individuals with major neuropsychiatric disorders such as schizophrenia and autism spectrum disorder have gray matter microstructural abnormalities. These include abnormalities in neuropil organization, expression of proteins supporting neuritic and synaptic integrity, and myelination. Genetic and postmortem studies suggest that these changes may be causally linked to the pathogenesis of these disorders. Advances in diffusion-weighted magnetic resonance image (dMRI) acquisition techniques and biophysical modeling allow for the quantification of gray matter microstructure in vivo. While several biophysical models for imaging microstructural properties are available, one in particular, neurite orientation dispersion and density imaging (NODDI), holds great promise for clinical applications. NODDI can be applied to both gray and white matter and requires only a single extra shell beyond a standard dMRI acquisition. Since its development only a few years ago, the NODDI algorithm has been used to characterize gray matter microstructure in schizophrenia, Alzheimer's disease, healthy aging, and development. These investigations have shown that microstructural findings in vivo, using NODDI, align with postmortem findings. Not only do NODDI and other advanced dMRI-based modeling methods provide a window into the brain previously only available postmortem, but they may be more sensitive to certain brain changes than conventional magnetic resonance imaging approaches. This opens up exciting new possibilities for clinicians to more rapidly detect disease signatures and allows earlier intervention in the course of the disease. Given that neurites and gray matter microstructure have the capacity to rapidly remodel, these novel dMRI-based methods represent an opportunity to noninvasively monitor neuroplastic changes posttherapy within much shorter time scales.