High-angular resolution diffusion imaging generation using 3d u-net.

PURPOSE: To investigate the effects on tractography of artificial intelligence-based prediction of motion-probing gradients (MPGs) in diffusion-weighted imaging (DWI). METHODS: The 251 participants in this study were patients with brain tumors or epileptic seizures who underwent MRI to depict tractography. DWI was performed with 64 MPG directions and b = 0 s/mm2 images. The dataset was divided into a training set of 191 (mean age 45.7 [± 19.1] years), a validation set of 30 (mean age 41.6 [± 19.1] years), and a test set of 30 (mean age 49.6 [± 18.3] years) patients. Supervised training of a convolutional neural network was performed using b = 0 images and the first 32 axes of MPG images as the input data and the second 32 axes as the reference data. The trained model was applied to the test data, and tractography was performed using (a) input data only; (b) input plus prediction data; and (c) b = 0 images and the 64 MPG data (as a reference). RESULTS: In Q-ball imaging tractography, the average dice similarity coefficient (DSC) of the input plus prediction data was 0.715 (± 0.064), which was significantly higher than that of the input data alone (0.697 [± 0.070]) (p < 0.05). In generalized q-sampling imaging tractography, the average DSC of the input plus prediction data was 0.769 (± 0.091), which was also significantly higher than that of the input data alone (0.738 [± 0.118]) (p < 0.01). CONCLUSION: Diffusion tractography is improved by adding predicted MPG images generated by an artificial intelligence model.

General orientation transform for the estimation of fiber orientations in white matter tissues.

PURPOSE: The orientation distribution function (ODF), which is obtained from the radial integral of the probability density function weighted by rn$$ {r}^n $$ ( r$$ r $$ is the radial length), has been used to estimate fiber orientations of white matter tissues. Currently, there is no general expression of the ODF that is suitable for any n value in the HARDI methods. THEORY AND METHODS: A novel methodology is proposed to calculate the ODF for any n>-1$$ n>-1 $$ through the Taylor series expansion and a generalized expression for -1

Delineating the Decussating Dentato-rubro-thalamic Tract and Its Connections in Humans Using Diffusion Spectrum Imaging Techniques.

The objective of this study was to identify the decussating dentato-rubro-thalamic tract (d-DRTT) and its afferent and efferent connections in healthy humans using diffusion spectrum imaging (DSI) techniques. In the present study, the trajectory and lateralization of the d-DRTT was explored using data from subjects in the Massachusetts General Hospital-Human Connectome Project adult diffusion dataset. The afferent and efferent networks that compose the cerebello-thalamo-cerebral pathways were also reconstructed. Correlation analysis was performed to identify interrelationships between subdivisions of the cerebello-dentato-rubro-thalamic and thalamo-cerebral connections. The d-DRTT was visualized bilaterally in 28 subjects. According to a normalized quantitative anisotropy and lateralization index evaluation, the left and right d-DRTT were relatively symmetric. Afferent regions were found mainly in the posterior cerebellum, especially the entire lobule VII (crus I, II and VIIb). Efferent fibers mainly are projected to the contralateral frontal cortex, including the motor and nonmotor regions. Correlations between cerebello-thalamic connections and thalamo-cerebral connections were positive, including the lobule VIIa (crus I and II) to the medial prefrontal cortex (MPFC) and the dorsolateral prefrontal cortex and lobules VI, VIIb, VIII, and IX, to the MPFC and motor and premotor areas. These results provide DSI-based tratographic evidence showing segregated and parallel cerebellar outputs to cerebral regions. The posterior cerebellum may play an important role in supporting and handling cognitive activities through d-DRTT. Future studies will allow for a more comprehensive understanding of cerebello-cerebral connections.

Identifying and characterizing projections from the subthalamic nucleus to the cerebellum in humans.

A connection between the subthalamic nucleus (STN) and the cerebellum which has been shown to exist in non-human primates, was recently identified in humans. However, its anatomical features, network properties and function have yet to be elucidated in humans. In the present study, we quantified the STN-cerebellum pathway in humans and explored its function based on structural observations. Anatomical features and asymmetry index (AI) were explored using high definition fiber tractography data of 30 individuals from the Massachusetts General Hospital - Human Connectome Project adult diffusion database. Pearson's correlation analysis was performed to determine the interrelationship between the subdivisions of the STN-cerebellum and the global cortical-STN connections. The pathway was visualized bilaterally in all the subjects. Typically, after setting out from the STN, the STN-cerebellum projections incorporated into the nearby corticopontine tracts, passing through the cerebral peduncle, mediated by the pontine nucleus and then connecting in two opposite directions to join the bilateral middle cerebellar peduncle. On the group averaged level, 78.03% and 62.54% of fibers from the right and left STN respectively, distributed to Crus I in the cerebellum, part of the remaining fibers projected to Crus II, with most of the fibers crossing contralaterally. According to the AI evaluation, 60% of the participants were right STN dominant, 23% were left STN dominant, and 17% were relatively symmetric. Pearson's correlation analysis further indicated that the number of pathways from mesial Brodmann area 8 to the STN (hyperdirect pathway associated with decision making) was positively correlated with the number of fibers from the right STN to Crus I. The insertion and termination, the right-side dominance, and the positive correlation with the hyperdirect pathway all suggest that the STN-cerebellum pathway might be involved in decision-making processes.

Characterizing the role of the structural connectome in seizure dynamics.

How does the human brain's structural scaffold give rise to its intricate functional dynamics? This is a central question in translational neuroscience that is particularly relevant to epilepsy, a disorder affecting over 50 million subjects worldwide. Treatment for medication-resistant focal epilepsy is often structural-through surgery or laser ablation-but structural targets, particularly in patients without clear lesions, are largely based on functional mapping via intracranial EEG. Unfortunately, the relationship between structural and functional connectivity in the seizing brain is poorly understood. In this study, we quantify structure-function coupling, specifically between white matter connections and intracranial EEG, across pre-ictal and ictal periods in 45 seizures from nine patients with unilateral drug-resistant focal epilepsy. We use high angular resolution diffusion imaging (HARDI) tractography to construct structural connectivity networks and correlate these networks with time-varying broadband and frequency-specific functional networks derived from coregistered intracranial EEG. Across all frequency bands, we find significant increases in structure-function coupling from pre-ictal to ictal periods. We demonstrate that short-range structural connections are primarily responsible for this increase in coupling. Finally, we find that spatiotemporal patterns of structure-function coupling are highly stereotyped for each patient. These results suggest that seizures harness the underlying structural connectome as they propagate. Mapping the relationship between structural and functional connectivity in epilepsy may inform new therapies to halt seizure spread, and pave the way for targeted patient-specific interventions.

Tractography and the connectome in neurosurgical treatment of gliomas: the premise, the progress, and the potential.

The ability of diffusion tensor MRI to detect the preferential diffusion of water in cerebral white matter tracts enables neurosurgeons to noninvasively visualize the relationship of lesions to functional neural pathways. Although viewed as a research tool in its infancy, diffusion tractography has evolved into a neurosurgical tool with applications in glioma surgery that are enhanced by evolutions in crossing fiber visualization, edema correction, and automated tract identification. In this paper the current literature supporting the use of tractography in brain tumor surgery is summarized, highlighting important clinical studies on the application of diffusion tensor imaging (DTI) for preoperative planning of glioma resection, and risk assessment to analyze postoperative outcomes. The key methods of tractography in current practice and crucial white matter fiber bundles are summarized. After a review of the physical basis of DTI and post-DTI tractography, the authors discuss the methodologies with which to adapt DT image processing for surgical planning, as well as the potential of connectomic imaging to facilitate a network approach to oncofunctional optimization in glioma surgery.

White matter degeneration profile in the cognitive cortico-subcortical tracts in Parkinson's disease.

BACKGROUND: In Parkinson's disease cognitive impairment is an early nonmotor feature, but it is still unclear why some patients are able to maintain their cognitive performance at normal levels, as quantified by neuropsychological tests, whereas others cannot. The objectives of this study were to perform a cross-sectional study and analyze the white matter changes in the cognitive and motor bundles in patients with Parkinson's disease. METHODS: Sixteen Parkinson's disease patients with normal cognitive performance, 19 with mild cognitive impairment (based on their performance of 1.5 standard deviations below the healthy population mean), and 16 healthy controls were compared with respect to their tractography patterns between the cortical cognitive / motor regions and subcortical structures, using high angular resolution diffusion imaging and constrained spherical deconvolution computation. RESULTS: Motor bundles showed decreased apparent fiber density in both PD groups, associated with a significant increase in diffusivity metrics, number of reconstructed streamlines, and track volumes, compared with healthy controls. By contrast, in the cognitive bundles, decreased fiber density in both Parkinson's groups was compounded by the absence of changes in diffusivity in patients with normal cognition, whereas patients with cognitive impairment had increased diffusivity metrics, lower numbers of reconstructed streamlines, and lower track volumes. CONCLUSIONS: Both PD groups showed similar patterns of white matter neurodegeneration in the motor bundles, whereas cognitive bundles showed a distinct pattern: Parkinson's patients with normal cognition had white matter diffusivity metrics similar to healthy controls, whereas in patients with cognitive impairment white matter showed a neurodegeneration pattern. © 2018 International Parkinson and Movement Disorder Society.

Prospective motion correction in diffusion-weighted imaging using intermediate pseudo-trace-weighted images.

Subject head motion is a major challenge in diffusion-weighted imaging, which requires a precise alignment of images from different time points to allow a reliable quantification of diffusion parameters within each voxel. The technique requires long measurement times, making it highly sensitive to long-term subject motion, even when head restraint is used. Current methods of data analysis rely on retrospective motion correction, but there are potential benefits to using prospective motion correction, in which motion is tracked and compensated for during data acquisition. This technique is regularly used to enhance image quality in blood-oxygen-level dependent (BOLD) imaging, but its application to diffusion-weighted imaging has been limited by the contrast variation between images acquired with different diffusion-gradient directions. This paper describes a novel approach to this topic that exploits the rotational invariance of the trace of the diffusion tensor to reduce the effect of this contrast variation, making it possible to perform a fast image registration using a least-squares cost function. This results in an image-based motion detection algorithm that can be applied in real time during data acquisition to adapt the slice position and orientation in response to subject motion. The motion detection capabilities of the technique were evaluated in a study of ten subjects with b-values up to 3000s/mm². The resulting motion-parameter estimates were in close agreement with reference values provided by interleaved low-b-value images with a correlation coefficient of R=0.9634 for the voxel displacements measured across all subjects and b-values. The technique was also used to perform prospective motion correction on a standard clinical MRI system with b-values up to 2000s/mm². The correction was evaluated in 3 subjects using interleaved low-b-value images, retrospective image registration using the AFNI processing package and mean diffusivity histogram analysis. Compared to acquisitions without motion correction, prospective motion correction based on pseudo-trace-weighted images was found to provide a robust method for substantially reducing the level of misregistration between volumes. In most cases, misregistrations were reduced to less than 0.2mm of translation and 0.2° of rotation for an isotropic voxel size of 2mm, yielding high-quality diffusion parameter maps even in the absence of head restraint and post-acquisition image registration.

Fast diffusion imaging with high angular resolution.

PURPOSE: High angular resolution diffusion imaging (HARDI) is a well-established method to help reveal the architecture of nerve bundles, but long scan times and geometric distortions inherent to echo planar imaging (EPI) have limited its integration into clinical protocols. METHODS: A fast imaging method is proposed here that combines accelerated multishot diffusion imaging (AMDI), multiplexed sensitivity encoding (MUSE), and crossing fiber angular resolution of intravoxel structure (CFARI) to reduce spatial distortions and reduce total scan time. A multishot EPI sequence was used to improve geometrical fidelity as compared to a single-shot EPI acquisition, and acceleration in both k-space and diffusion sampling enabled reductions in scan time. The method is regularized and self-navigated for motion correction. Seven volunteers were scanned in this study, including four with volumetric whole brain acquisitions. RESULTS: The average similarity of microstructural orientations between undersampled datasets and their fully sampled counterparts was above 85%, with scan times below 5 min for whole-brain acquisitions. Up to 2.7-fold scan time acceleration along with four-fold distortion reduction was achieved. CONCLUSION: The proposed imaging strategy can generate HARDI results with relatively good geometrical fidelity and low scan duration, which may help facilitate the transition of HARDI from a successful research tool to a practical clinical one. Magn Reson Med 77:696-706, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

The importance of correcting for signal drift in diffusion MRI.

PURPOSE: To investigate previously unreported effects of signal drift as a result of temporal scanner instability on diffusion MRI data analysis and to propose a method to correct this signal drift. METHODS: We investigated the signal magnitude of non-diffusion-weighted EPI volumes in a series of diffusion-weighted imaging experiments to determine whether signal magnitude changes over time. Different scan protocols and scanners from multiple vendors were used to verify this on phantom data, and the effects on diffusion kurtosis tensor estimation in phantom and in vivo data were quantified. Scalar metrics (eigenvalues, fractional anisotropy, mean diffusivity, mean kurtosis) and directional information (first eigenvectors and tractography) were investigated. RESULTS: Signal drift, a global signal decrease with subsequently acquired images in the scan, was observed in phantom data on all three scanners, with varying magnitudes up to 5% in a 15-min scan. The signal drift has a noticeable effect on the estimation of diffusion parameters. All investigated quantitative parameters as well as tractography were affected by this artifactual signal decrease during the scan. CONCLUSION: By interspersing the non-diffusion-weighted images throughout the session, the signal decrease can be estimated and compensated for before data analysis; minimizing the detrimental effects on subsequent MRI analyses. Magn Reson Med 77:285-299, 2017. © 2016 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine.

Fiber ball imaging.

By modeling axons as thin cylinders, it is shown that the inverse Funk transform of the diffusion MRI (dMRI) signal intensity obtained on a spherical shell in q-space gives an estimate for a fiber orientation density function (fODF), where the accuracy improves with increasing b-value provided the signal-to-noise ratio is sufficient. The method is similar to q-ball imaging, except that the Funk transform of q-ball imaging is replaced by its inverse. We call this new approach fiber ball imaging. The fiber ball method is demonstrated for healthy human brain, and fODF estimates are compared to diffusion orientation distribution function (dODF) approximations obtained with q-ball imaging. The fODFs are seen to have sharper features than the dODFs, reflecting an enhancement of the higher degree angular frequencies. The inverse Funk transform of the dMRI signal intensity data provides a simple and direct method of estimating a fODF. In addition, fiber ball imaging leads to an estimate for the ratio of the fraction of MRI visible water confined to the intra-axonal space divided by the square root of the intra-axonal diffusivity. This technique may be useful for white matter fiber tractography, as well as other types of microstructural modeling of brain tissue.

Diffusivity signatures characterize trigeminal neuralgia associated with multiple sclerosis.

BACKGROUND: Trigeminal neuralgia secondary to multiple sclerosis (MS-TN) is a facial neuropathic pain syndrome similar to classic trigeminal neuralgia (TN). While TN is caused by neurovascular compression of the fifth cranial nerve (CN V), how MS-related demyelination correlates with pain in MS-TN is not understood. OBJECTIVES: We aim to examine diffusivities along CN V in MS-TN, TN, and controls in order to reveal differential neuroimaging correlates across groups. METHODS: 3T MR diffusion weighted, T1, T2 and FLAIR sequences were acquired for MS-TN, TN, and controls. Multi-tensor tractography was used to delineate CN V across cisternal, root entry zone (REZ), pontine and peri-lesional segments. Diffusion metrics including fractional anisotropy (FA), and radial (RD), axial (AD), and mean diffusivities (MD) were measured from each segment. RESULTS: CN V segments showed distinctive diffusivity patterns. The TN group showed higher FA in the cisternal segment ipsilateral to the side of pain, and lower FA in the ipsilateral REZ segment. The MS-TN group showed lower FA in the ipsilateral peri-lesional segments, suggesting differential microstructural changes along CN V in these conditions. CONCLUSIONS: The study demonstrates objective differences in CN V microstrucuture in TN and MS-TN using non-invasive neuroimaging. This represents a significant improvement in the methods currently available to study pain in MS.

Advanced diffusion MRI fiber tracking in neurosurgical and neurodegenerative disorders and neuroanatomical studies: A review.

Diffusion MRI enabled in vivo microstructural imaging of the fiber tracts in the brain resulting in its application in a wide range of settings, including in neurological and neurosurgical disorders. Conventional approaches such as diffusion tensor imaging (DTI) have been shown to have limited applications due to the crossing fiber problem and the susceptibility of their quantitative indices to partial volume effects. To overcome these limitations, the recent focus has shifted to the advanced acquisition methods and their related analytical approaches. Advanced white matter imaging techniques provide superior qualitative data in terms of demonstration of multiple crossing fibers in their spatial orientation in a three dimensional manner in the brain. In this review paper, we discuss the advancements in diffusion MRI and introduce their roles. Using examples, we demonstrate the role of advanced diffusion MRI-based fiber tracking in neuroanatomical studies. Results from its preliminary application in the evaluation of intracranial space occupying lesions, including with respect to future directions for prognostication, are also presented. Building upon the previous DTI studies assessing white matter disease in Huntington's disease and Amyotrophic lateral sclerosis; we also discuss approaches which have led to encouraging preliminary results towards developing an imaging biomarker for these conditions.

A new compression format for fiber tracking datasets.

A single diffusion MRI streamline fiber tracking dataset may contain hundreds of thousands, and often millions of streamlines and can take up to several gigabytes of memory. This amount of data is not only heavy to compute, but also difficult to visualize and hard to store on disk (especially when dealing with a collection of brains). These problems call for a fiber-specific compression format that simplifies its manipulation. As of today, no fiber compression format has yet been adopted and the need for it is now becoming an issue for future connectomics research. In this work, we propose a new compression format, .zfib, for streamline tractography datasets reconstructed from diffusion magnetic resonance imaging (dMRI). Tracts contain a large amount of redundant information and are relatively smooth. Hence, they are highly compressible. The proposed method is a processing pipeline containing a linearization, a quantization and an encoding step. Our pipeline is tested and validated under a wide range of DTI and HARDI tractography configurations (step size, streamline number, deterministic and probabilistic tracking) and compression options. Similar to JPEG, the user has one parameter to select: a worst-case maximum tolerance error in millimeter (mm). Overall, we find a compression factor of more than 96% for a maximum error of 0.1mm without any perceptual change or change of diffusion statistics (mean fractional anisotropy and mean diffusivity) along bundles. This opens new opportunities for connectomics and tractometry applications.

Acceleration of high angular and spatial resolution diffusion imaging using compressed sensing with multichannel spiral data.

PURPOSE: To accelerate the acquisition of simultaneously high spatial and angular resolution diffusion imaging. METHODS: Accelerated imaging is achieved by recovering the diffusion signal at all voxels simultaneously from under-sampled k-q space data using a compressed sensing algorithm. The diffusion signal at each voxel is modeled as a sparse complex Gaussian mixture model. The joint recovery scheme enables incoherent under-sampling of the 5-D k-q space, obtained by randomly skipping interleaves of a multishot variable density spiral trajectory. This sampling and reconstruction strategy is observed to provide considerably improved reconstructions than classical k-q under-sampling and reconstruction schemes. The complex model enables to account for the noise statistics without compromising the computational efficiency and theoretical convergence guarantees. The reconstruction framework also incorporates compensation of motion induced phase errors that result from the multishot acquisition. RESULTS: Reconstructions of the diffusion signal from under-sampled data using the proposed method yields accurate results with errors less that 5% for different accelerations and b-values. The proposed method is also shown to perform better than standard k-q acceleration schemes. CONCLUSIONS: The proposed scheme can significantly accelerate the acquisition of high spatial and angular resolution diffusion imaging by accurately reconstructing crossing fiber architectures from under-sampled data.

Development of a high angular resolution diffusion imaging human brain template.

Brain diffusion templates contain rich information about the microstructure of the brain, and are used as references in spatial normalization or in the development of brain atlases. The accuracy of diffusion templates constructed based on the diffusion tensor (DT) model is limited in regions with complex neuronal micro-architecture. High angular resolution diffusion imaging (HARDI) overcomes limitations of the DT model and is capable of resolving intravoxel heterogeneity. However, when HARDI is combined with multiple-shot sequences to minimize image artifacts, the scan time becomes inappropriate for human brain imaging. In this work, an artifact-free HARDI template of the human brain was developed from low angular resolution multiple-shot diffusion data. The resulting HARDI template was produced in ICBM-152 space based on Turboprop diffusion data, was shown to resolve complex neuronal micro-architecture in regions with intravoxel heterogeneity, and contained fiber orientation information consistent with known human brain anatomy.

Quantification of track-weighted imaging (TWI): characterisation of within-subject reproducibility and between-subject variability.

Recently several novel image contrasts derived from whole-brain fibre tracking-data (tractograms) have been introduced. The novel contrasts of these track-weighted imaging (TWI) methods may provide important information for clinical neuroscience studies. However, before they can be used reliably to generate quantitative measures, it is important to characterise their within-subject reproducibility, and between-subject variability. In this work we compute the within-subject reproducibility (intra-scan, intra-session and inter-session), and between-subject variability of TWI for a number of different TWI contrasts across multiple subjects. The results are used in simple voxel-wise power calculations within illustrative regions of interest to provide guidelines for required sample sizes and observable effect sizes for individual subjects and between groups. It was found that the required sample sizes and observable effect sizes varied considerably between different TWI maps and for different ROIs. For some TWI contrast and ROI combinations, the power calculations yielded clinically practical values. These results provide important information concerning the potential usefulness and sensitivity of TWI maps for individual diagnosis, longitudinal studies and group comparisons, as well as for study designs.

Abnormalities of cortical thickness, subcortical shapes, and white matter integrity in subcortical vascular cognitive impairment.

Subcortical vascular cognitive impairment (sVCI) is caused by lacunar infarcts or extensive and/or diffuse lesions in the white matter that may disrupt the white matter circuitry connecting cortical and subcortical regions and result in the degeneration of neurons in these regions. This study used structural magnetic resonance imaging (MRI) and high angular resolution diffusion imaging (HARDI) techniques to examine cortical thickness, subcortical shapes, and white matter integrity in mild vascular cognitive impairment no dementia (VCIND Mild) and moderate-to-severe VCI (MSVCI). Our study found that compared to controls (n = 25), VCIND Mild (n = 25), and MSVCI (n = 30) showed thinner cortex predominantly in the frontal cortex. The cortex in MSVCI was thinner in the parietal and lateral temporal cortices than that in VCIND Mild. Moreover, compared to controls, VCIND Mild and MSVCI showed smaller shapes (i.e., volume reduction) in the thalamus, putamen, and globus pallidus and ventricular enlargement. Finally, compared to controls, VCIND Mild, and MSVCI showed an increased mean diffusivity in the white matter, while decreased generalized fractional anisotropy was only found in the MSVCI subjects. The major axonal bundles involved in the white matter abnormalities were mainly toward the frontal regions, including the internal capsule/corona radiata, uncinate fasciculus, and anterior section of the inferior fronto-occipital fasciculus, and were anatomically connected to the affected cortical and subcortical structures. Our findings suggest that abnormalities in cortical, subcortical, and white matter morphology in sVCI occur in anatomically connected structures, and that abnormalities progress along a similar trajectory from the mild to moderate and severe conditions.

Structural and functional connectivity of the human brain in autism spectrum disorders and attention-deficit/hyperactivity disorder: A rich club-organization study.

Attention-deficit/hyperactive disorder (ADHD) and autism spectrum disorders (ASD) are two of the most common and vexing neurodevelopmental disorders among children. Although the two disorders share many behavioral and neuropsychological characteristics, most MRI studies examine only one of the disorders at a time. Using graph theory combined with structural and functional connectivity, we examined the large-scale network organization among three groups of children: a group with ADHD (8-12 years, n = 20), a group with ASD (7-13 years, n = 16), and typically developing controls (TD) (8-12 years, n = 20). We apply the concept of the rich-club organization, whereby central, highly connected hub regions are also highly connected to themselves. We examine the brain into two different network domains: (1) inside a rich-club network phenomena and (2) outside a rich-club network phenomena. The ASD and ADHD groups had markedly different patterns of rich club and non rich-club connections in both functional and structural data. The ASD group exhibited higher connectivity in structural and functional networks but only inside the rich-club networks. These findings were replicated using the autism brain imaging data exchange dataset with ASD (n = 85) and TD (n = 101). The ADHD group exhibited a lower generalized fractional anisotropy and functional connectivity inside the rich-club networks, but a higher number of axonal fibers and correlation coefficient values outside the rich club. Despite some shared biological features and frequent comorbity, these data suggest ADHD and ASD exhibit distinct large-scale connectivity patterns in middle childhood.

Estimation of the CSA-ODF using Bayesian compressed sensing of multi-shell HARDI.

PURPOSE: Diffusion MRI provides important information about the brain white matter structures and has opened new avenues for neuroscience and translational research. However, acquisition time needed for advanced applications can still be a challenge in clinical settings. There is consequently a need to accelerate diffusion MRI acquisitions. METHODS: A multi-task Bayesian compressive sensing (MT-BCS) framework is proposed to directly estimate the constant solid angle orientation distribution function (CSA-ODF) from under-sampled (i.e., accelerated image acquisition) multi-shell high angular resolution diffusion imaging (HARDI) datasets, and accurately recover HARDI data at higher resolution in q-space. The proposed MT-BCS approach exploits the spatial redundancy of the data by modeling the statistical relationships within groups (clusters) of diffusion signal. This framework also provides uncertainty estimates of the computed CSA-ODF and diffusion signal, directly computed from the compressive measurements. Experiments validating the proposed framework are performed using realistic multi-shell synthetic images and in vivo multi-shell high angular resolution HARDI datasets. RESULTS: Results indicate a practical reduction in the number of required diffusion volumes (q-space samples) by at least a factor of four to estimate the CSA-ODF from multi-shell data. CONCLUSION: This work presents, for the first time, a multi-task Bayesian compressive sensing approach to simultaneously estimate the full posterior of the CSA-ODF and diffusion-weighted volumes from multi-shell HARDI acquisitions. It demonstrates improvement of the quality of acquired datasets by means of CS de-noising, and accurate estimation of the CSA-ODF, as well as enables a reduction in the acquisition time by a factor of two to four, especially when "staggered" q-space sampling schemes are used. The proposed MT-BCS framework can naturally be combined with parallel MR imaging to further accelerate HARDI acquisitions.