Search

Diffusion imaging is a critical component in the pursuit of developing a better understanding of the human brain. Recent technical advances promise enabling the advancement in the quality of data that can be obtained. In this review the context for different approaches relative to the Human Connectome Project are compared. Significant new gains are anticipated from the use of high-performance head gradients. These gains can be particularly large when the high-performance gradients are employed together with ultrahigh magnetic fields. Transmit array designs are critical in realizing high accelerations in diffusion-weighted (d)MRI acquisitions, while maintaining large field of view (FOV) coverage, and several techniques for optimal signal-encoding are now available. Reconstruction and processing pipelines that precisely disentangle the acquired neuroanatomical information are established and provide the foundation for the application of deep learning in the advancement of dMRI for complex tissues. Level of Evidence: 3 Technical Efficacy Stage: Stage 3.

Purpose In this study, we sought to develop a self-navigation strategy for improving the reconstruction of diffusion weighted 3D multishot echo planar imaging (EPI). We propose a method for extracting the phase correction information from the acquisition itself, eliminating the need for a 2D navigator, further accelerating the acquisition. Methods In-vivo acquisitions at 3T with 0.9 mm and 1.5 mm isotropic resolutions were used to evaluate the performance of the self-navigation strategy. Sensitivity to motion was tested using a large difference in pitch position of the head. Using a multishell diffusion weighted acquisition, tractography results were obtained at (0.9 mm)3 to validate the quality with conventional acquisition. Results The use of 3D multislab EPI with self-navigation enables 3D diffusion-weighted spin echo EPI acquisitions that have the same efficiency as 2D single-shot acquisition. For matched acquisition time the image signal-to-noise ratio (SNR) between 3D and 2D acquisition is shown to be comparable for whole-brain coverage with (1.5 mm)3 resolution and for (0.9 mm)3 resolution the 3D acquisition has higher SNR than what can be obtained with 2D acquisitions using current state-of-art multiband techniques. The self-navigation technique was shown to be stable under inter-volume motion. In tractography analysis, the higher resolution afforded by our technique enabled clear delineation of the tapetum and posterior corona radiata. Conclusion The proposed self-navigation approach utilized a self-consistent phase in 3D diffusion weighted acquisitions. Its efficiency and stability were demonstrated for a plurality of common acquisitions. The proposed self-navigation approach allows for faster acquisition of 3D multishot EPI desirable for large field of view and/or higher resolution.

The mechanisms contributing to increased expression of motor and cognitive impairment in people with Parkinson9s disease who lack muscle atonia during rapid eye movement (REM) sleep compared to those with muscle atonia are poorly understood. This study used tract-based spatial statistics to compare diffusion measures of white matter microstructure between people with Parkinson9s disease with and without REM sleep without atonia as well as the relationships of these measures to motor and cognitive function. Thirty-eight individuals with mild-to-moderate Parkinson9s disease and twenty-one matched control subjects underwent ultra-high-field MRI (7Tesla), quantitative motor assessments of gait and bradykinesia, and neuropsychological testing. The Parkinson9s disease cohort was separated post-hoc into those with and without elevated chin or leg muscle activity during REM sleep based on polysomnography findings. Fractional anisotropy was significantly higher, and radial and mean diffusivity significantly lower, in diffuse regions of the corpus callosum, projection, and association white matter pathways in the Parkinson9s group with normal REM sleep compared to controls. In contrast, there was no significant difference in fractional anisotropy between the Parkinson9s group with elevated muscle tone and controls. Fractional anisotropy was also significantly higher in a subset of pathways in the Parkinson9s disease group with normal REM sleep muscle tone compared to those with elevated REM sleep muscle tone. The group with elevated REM sleep muscle tone had significant impairments in gait and upper arm speed compared to controls and significantly worse scores in specific cognitive domains (executive function, visuospatial memory) compared to the Parkinson9s disease group with normal REM sleep muscle tone. Regression analyses showed that gait speed and step length in the Parkinson9s disease cohort were predicted by measures of mean fractional anisotropy of the anterior corona radiata, whereas elbow flexion velocity was predicted by fractional anisotropy of the superior corona radiata. Visuospatial memory task performance was predicted by the radial diffusivity of the posterior corona radiata. These findings demonstrate that people with mild-to-moderate severity of Parkinson9s disease who have normal muscle tone during REM sleep show alterations in white matter microstructure that are associated with preserved motor and cognitive function, but these adaptations are reduced or absent in those with increased REM sleep motor tone.

Investigation of brain changes in functional connectivity and functional network topology from receiving 8-week selective serotonin reuptake inhibitor (SSRI) treatments is conducted in 12 unmedicated adolescents with major depressive disorder (MDD) by using wavelet-filtered resting-state functional magnetic resonance imaging (fMRI). Changes are observed in frontal-limbic, temporal, and default mode networks. In particular, topological analysis shows, at the global scale and in the 0.12–0.25 Hz band, that the normalized clustering coefficient and smallworldness of brain networks decreased after treatment. Regional changes in clustering coefficient and efficiency were observed in the bilateral caudal middle frontal gyrus, rostral middle frontal gyrus, superior temporal gyrus, left pars triangularis, putamen, and right superior frontal gyrus. Furthermore, changes of nodal centrality and changes of connectivity associated with these frontal and temporal regions confirm the global topological alternations. Moreover, frequency dependence is observed from FDR-controlled subnetworks for the limbic-cortical connectivity change. In the high-frequency band, the altered connections involve mostly frontal regions, while the altered connections in the low-frequency bands spread to parietal and temporal areas. Due to the limitation of small sample sizes and lack of placebo control, these preliminary findings require confirmation with future work using larger samples. Confirmation of biomarkers associated with treatment could suggest potential avenues for clinical applications such as tracking treatment response and neurobiologically informed treatment optimization.

Mild traumatic brain injury (mTBI) is a significant cause of disability, especially when symptoms become chronic. This chronicity is often linked to oculomotor dysfunction (OMD). To our knowledge, this is the first prospective study to localize aberrations in brain function between mTBI cohorts, by comparing patients with mTBI with OMD with an mTBI control group without OMD, using task and resting-state functional magnetic resonance imaging (fMRI). Ten subjects with mTBI who had OMD (OMD group) were compared with nine subjects with mTBI who had no findings of OMD (control group). These groups were determined by a developmental optometrist using objective testing for OMD. The (convergence) task fMRI data demonstrated significantly decreased brain activity, measured as decreases in the blood oxygen level dependent (BOLD) signal, in the OMD group compared with the control group in three brain regions: the left posterior lingual gyrus, the bilateral anterior lingual gyrus and cuneus, and the parahippocampal gyrus. When doing a seed-based resting state fMRI analysis in the lingual/parahippocampal region, a large cluster covering the left middle frontal gyrus and the dorsolateral pre-frontal cortex (Brodmann areas 9 and 10), with decreased functional correlation in the OMD group, was identified. Together these observations provide evidence for neural networks of interactions involving the control of eye movement for visual processing, reading comprehension, spatial localization and navigation, and spatial working memory that appear to be decreased in mTBI patients with OMD compared with mTBI patients without OMD. The clinical symptomatology associated with post-traumatic OMD correlates well with these MRI findings.

BACKGROUND AND PURPOSE Non-myelopathic degenerative cervical spinal cord compression (NMDC) frequently occurs throughout aging and may progress to potentially irreversible degenerative cervical myelopathy (DCM). Whereas standard clinical magnetic resonance imaging (MRI) and electrophysiological measures assess compression severity and neurological dysfunction, respectively, underlying microstructural deficits still have to be established in NMDC and DCM patients. The study aims to establish tract-specific diffusion MRI markers of electrophysiological deficits to predict the progression of asymptomatic NMDC to symptomatic DCM. METHODS High-resolution 3 T diffusion MRI was acquired for 103 NMDC and 21 DCM patients compared to 60 healthy controls to reveal diffusion alterations and relationships between tract-specific diffusion metrics and corresponding electrophysiological measures and compression severity. Relationship between the degree of DCM disability, assessed by the modified Japanese Orthopaedic Association scale, and tract-specific microstructural changes in DCM patients was also explored. RESULTS The study identified diffusion-derived abnormalities in the gray matter, dorsal and lateral tracts congruent with trans-synaptic degeneration and demyelination in chronic degenerative spinal cord compression with more profound alterations in DCM than NMDC. Diffusion metrics were affected in the C3-6 area as well as above the compression level at C3 with more profound rostral deficits in DCM than NMDC. Alterations in lateral motor and dorsal sensory tracts correlated with motor and sensory evoked potentials, respectively, whereas electromyography outcomes corresponded with gray matter microstructure. DCM disability corresponded with microstructure alteration in lateral columns. CONCLUSIONS Outcomes imply the necessity of high-resolution tract-specific diffusion MRI for monitoring degenerative spinal pathology in longitudinal studies.

Objective: Friedreich ataxia (FRDA) is an inherited neurological disease defined by progressive movement incoordination. We undertook a comprehensive characterization of the spatial profile and progressive evolution of structural brain abnormalities in people with FRDA. Methods: A coordinated international analysis of regional brain volume using magnetic resonance imaging data charted the whole-brain profile, interindividual variability, and temporal staging of structural brain differences in 248 individuals with FRDA and 262 healthy controls. Results: The brainstem, dentate nucleus region, and superior and inferior cerebellar peduncles showed the greatest reductions in volume relative to controls (Cohen d= 1.5–2.6). Cerebellar gray matter alterations were most pronounced in lobules I–VI (d= 0.8), whereas cerebral differences occurred most prominently in precentral gyri (d= 0.6) and corticospinal tracts (d= 1.4). Earlier onset age predicted less volume in the motor cerebellum (rmax= 0.35) and peduncles (rmax= 0.36). Disease duration and severity correlated with volume deficits in the dentate nucleus region, brainstem, and superior/inferior cerebellar peduncles (rmax= −0.49); subgrouping showed these to be robust and early features of FRDA, and strong candidates for further biomarker validation. Cerebral white matter abnormalities, particularly in corticospinal pathways, emerge as intermediate disease features. Cerebellar and cerebral gray matter loss, principally targeting motor and sensory systems, preferentially manifests later in the disease course. Interpretation: FRDA is defined by an evolving spatial profile of neuroanatomical changes beyond primary pathology in the cerebellum and spinal cord, in line with its progressive clinical course. The design, interpretation, and generalization of research studies and clinical trials must consider neuroanatomical staging and associated interindividual variability in brain measures. ANN NEUROL 2021;90:570–583.

In a companion paper by Cohen-Adad et al. we introduce the spine generic quantitative MRI protocol that provides valuable metrics for assessing spinal cord macrostructural and microstructural integrity. This protocol was used to acquire a single subject dataset across 19 centers and a multi-subject dataset across 42 centers (for a total of 260 participants), spanning the three main MRI manufacturers: GE, Philips and Siemens. Both datasets are publicly available via git-annex. Data were analysed using the Spinal Cord Toolbox to produce normative values as well as inter/intra-site and inter/intra-manufacturer statistics. Reproducibility for the spine generic protocol was high across sites and manufacturers, with an average inter-site coefficient of variation of less than 5% for all the metrics. Full documentation and results can be found at https://spine-generic.rtfd.io/ . The datasets and analysis pipeline will help pave the way towards accessible and reproducible quantitative MRI in the spinal cord.

Spinocerebellar ataxia type 1 is a progressive neurodegenerative, movement disorder. With potential therapies on the horizon, it is critical to identify biomarkers that (i) differentiate between unaffected and spinocerebellar ataxia Type 1-affected individuals; (ii) track disease progression; and (iii) are directly related to clinical changes of the patient. Magnetic resonance imaging of volumetric changes in the brain may be a suitable source of biomarkers for spinocerebellar ataxia Type 1. In a previous report on a longitudinal study of patients with spinocerebellar ataxia Type 1, we evaluated the volume and magnetic resonance spectroscopy measures of the cerebellum and pons, showing pontine volume and pontine N-acetylaspartate-to-myo-inositol ratio were sensitive to change over time. As a follow-up, the current study conducts a whole brain exploration of volumetric MRI measures with the aim to identify biomarkers for spinocerebellar ataxia Type 1 progression. We adapted a joint label fusion approach using multiple, automatically generated, morphologically matched atlases to label brain regions including cerebellar sub-regions. We adjusted regional volumes by total intracranial volume allowing for linear and power-law relationships. We then utilized Bonferroni corrected linear mixed effects models to (i) determine group differences in regional brain volume and (ii) identify change within affected patients only. We then evaluated the rate of change within each brain region to identify areas that changed most rapidly. Lastly, we used a penalized, linear mixed effects model to determine the strongest brain predictors of motor outcomes. Decrease in pontine volume and accelerating decrease in putamen volume: (i) reliably differentiated spinocerebellar ataxia Type 1-affected and -unaffected individuals; (ii) were observable in affected individuals without referencing an unaffected comparison group; (iii) were detectable within ∼6–9 months; and (iv) were associated with increased disease burden. In conclusion, volumetric change in the pons and putamen may provide powerful biomarkers to track disease progression in spinocerebellar ataxia Type 1. The methods employed here are readily translatable to current clinical settings, providing a framework for study and usage of volumetric neuroimaging biomarkers for clinical trials., There is an urgent need for biomarkers that track neurodegenerative disease progression. Pontine and putamen volume decreases are specific to affected individuals, track disease progression within 6–9 months and are directly related to clinical changes in individuals affected by spinocerebellar ataxia Type 1, making them potential biomarkers for clinical trials., Graphical Abstract

The robustness of brain structural networks, estimated from diffusion MRI data, may be relevant to cognition. We investigate whether measures of network robustness, such as Ollivier-Ricci curvature, can explain cognitive impairment in multiple sclerosis (MS). We assessed whether local (i.e., cortical area) and/or global (i.e., whole brain) robustness, differs between cognitively impaired (MSCI) and non-impaired (MSNI) MS patients. Fifty patients, with Expanded Disability Status Scale mean (m): 3.2, disease duration m: 12 years, and age m: 40 years, were enrolled. Cognitive impairment scores were estimated from the Minimal Assessment of Cognitive Function in Multiple Sclerosis. Images were obtained in a 3T MRI using a diffusion protocol with a 2 min acquisition time. Brain structural networks were created using 333 cortical areas. Local and global robustness was estimated for each individual, and comparisons were performed between MSCI and MSNI patients. 31 MSCI and 10 MSNI patients were included in the analyses. Brain structural network robustness and centrality showed significant correlations with cognitive impairment. Measures of network robustness and centrality identified specific cortical areas relevant to MS-related cognitive impairment. These measures can be obtained on clinical scanners and are succinct yet accurate potential biomarkers of cognitive impairment.

Diffusion magnetic resonance imaging (dMRI) proved promising in patients with non-myelopathic degenerative cervical cord compression (NMDCCC), i.e., without clinically manifested myelopathy. Aim of the study is to present a fast multi-shell HARDI-ZOOMit dMRI protocol and validate its usability to detect microstructural myelopathy in NMDCCC patients. In 7 young healthy volunteers, 13 age-comparable healthy controls, 18 patients with mild NMDCCC and 15 patients with severe NMDCCC, the protocol provided higher signal-to-noise ratio, enhanced visualization of white/gray matter structures in microstructural maps, improved dMRI metric reproducibility, preserved sensitivity (SE = 87.88%) and increased specificity (SP = 92.31%) of control-patient group differences when compared to DTI-RESOLVE protocol (SE = 87.88%, SP = 76.92%). Of the 56 tested microstructural parameters, HARDI-ZOOMit yielded significant patient-control differences in 19 parameters, whereas in DTI-RESOLVE data, differences were observed in 10 parameters, with mostly lower robustness. Novel marker the white-gray matter diffusivity gradient demonstrated the highest separation. HARDI-ZOOMit protocol detected larger number of crossing fibers (5–15% of voxels) with physiologically plausible orientations than DTI-RESOLVE protocol (0–8% of voxels). Crossings were detected in areas of dorsal horns and anterior white commissure. HARDI-ZOOMit protocol proved to be a sensitive and practical tool for clinical quantitative spinal cord imaging.

Development of imaging biomarkers for rare neurodegenerative diseases such as spinocerebellar ataxia (SCA) is important to non-invasively track progression of disease pathology and monitor response to interventions. Diffusion MRI (dMRI) has been shown to identify cross-sectional degeneration of white matter (WM) microstructure and connectivity between healthy controls and patients with SCAs, using various analysis methods. In this paper, we present dMRI data in SCAs type 1, 2, 3, and 6 and matched controls, including longitudinal acquisitions at 12–24-month intervals in a subset of the cohort, with up to 5 visits. The SCA1 cohort also contained 3 premanifest patients at baseline, with 2 showing ataxia symptoms at the time of the follow-up scans. We focused on two aspects: first, multimodal evaluation of the dMRI data in a cross-sectional approach, and second, longitudinal trends in dMRI data in SCAs. Three different pipelines were used to perform cross-sectional analyses in WM: region of interest (ROI), tract-based spatial statistics (TBSS), and fixel-based analysis (FBA). We further analyzed longitudinal changes in dMRI metrics throughout the brain using ROI-based analysis. Both ROI and TBSS analyses identified higher mean (MD), axial (AD), and radial (RD) diffusivity and lower fractional anisotropy (FA) in the cerebellum for all SCAs compared to controls, as well as some cerebral alterations in SCA1, 2, and 3. FBA showed lower fiber density (FD) and fiber crossing (FC) regions similar to those identified by ROI and TBSS analyses. FBA also highlighted corticospinal tract (CST) abnormalities, which was not detected by the other two pipelines. Longitudinal ROI-based analysis showed significant increase in AD in the middle cerebellar peduncle (MCP) for patients with SCA1, suggesting that the MCP may be a good candidate region to monitor disease progression. The patient who remained symptom-free throughout the study displayed no microstructural abnormalities. On the other hand, the two patients who were at the premanifest stage at baseline, and showed ataxia symptoms in their follow-up visits, displayed AD values in the MCP that were already in the range of symptomatic patients with SCA1 at their baseline visit, demonstrating that microstructural abnormalities are detectable prior to the onset of ataxia.

The success of deep brain stimulation (DBS) surgeries for the treatment of movement disorders relies on the accurate placement of an electrode within the motor portion of subcortical brain targets. However, the high number of electrodes requiring relocation indicates that today's methods do not ensure sufficient accuracy for all patients. Here, with the goal of aiding DBS targeting, we use 7 Tesla (T) MRI data to identify the functional territories and parcellate the globus pallidus pars interna (GPi) into motor, associative and limbic regions in individual subjects. 7 T MRI scans were performed in seventeen patients (prior to DBS surgery) and one healthy control. Tractography-based parcellation of each patient's GPi was performed. The cortex was divided into four masks representing motor, limbic, associative and “other” regions. Given that no direct connections between the GPi and the cortex have been shown to exist, the parcellation was carried out in two steps: 1) The thalamus was parcellated based on the cortical targets, 2) The GPi was parcellated using the thalamus parcels derived from step 1. Reproducibility, via repeated scans of a healthy subject, and validity of the findings, using different anatomical pathways for parcellation, were assessed. Lastly, post-operative imaging data was used to validate and determine the clinical relevance of the parcellation. The organization of the functional territories of the GPi observed in our individual patient population agrees with that previously reported in the literature: the motor territory was located posterolaterally, followed anteriorly by the associative region, and further antero-ventrally by the limbic territory. While this organizational pattern was observed across patients, there was considerable variability among patients. The organization of the functional territories of the GPi was remarkably reproducible in intra-subject scans. Furthermore, the organizational pattern was observed consistently by performing the parcellation of the GPi via the thalamus and via a different pathway, going through the striatum. Finally, the active therapeutic contact of the DBS electrode, identified with a combination of post-operative imaging and post-surgery DBS programming, overlapped with the high-probability “motor” region of the GPi as defined by imaging-based methods. The consistency, validity, and clinical relevance of our findings have the potential for improving DBS targeting, by increasing patient-specific knowledge of subregions of the GPi to be targeted or avoided, at the stage of surgical planning, and later, at the stage when stimulation is adjusted.

PURPOSE: To develop a fast and automated volume-of-interest (VOI) prescription pipeline (AutoVOI) for single-voxel MR spectroscopy (MRS) that removes the need for manual VOI placement, allows flexible VOI planning in any brain region, and enables high inter- and intra-subject consistency of VOI prescription. METHODS: AutoVOI was designed to transfer pre-defined VOIs from an atlas to the 3D anatomical data of the subject during the scan. The AutoVOI pipeline was optimized for consistency in VOI placement (precision), enhanced coverage of the targeted tissue (accuracy) and fast computation speed. The tool was evaluated against manual VOI placement using existing T(1)-weighted datasets and corresponding VOI prescriptions. Finally, it was implemented on two scanner platforms to acquire MRS data from clinically-relevant VOIs that span the cerebrum, cerebellum and the brainstem. RESULTS: The AutoVOI pipeline includes skull stripping, non-linear registration of the atlas to the subject’s brain, and computation of the VOI coordinates and angulations using a minimum oriented bounding box algorithm. When compared against manual prescription, AutoVOI showed higher intra- and inter-subject spatial consistency, as quantified by generalized Dice coefficients (GDC), lower intra- and inter-subject variability in tissue composition (gray matter, white matter and cerebrospinal fluid) and higher or equal accuracy, as quantified by GDC of prescribed VOI with targeted tissues. High quality spectra were obtained on Siemens and Philips 3T systems from 6 automatically prescribed VOIs by the tool. CONCLUSION: Robust automatic VOI prescription is feasible and can help facilitate clinical adoption of MRS by avoiding operator dependence of manual selection.

We present a sparse Bayesian unmixing algorithm BusineX: Bayesian Unmixing for Sparse Inference-based Estimation of Fiber Crossings (X), for estimation of white matter fiber parameters from compressed (under-sampled) diffusion MRI (dMRI) data. BusineX combines compressive sensing with linear unmixing and introduces sparsity to the previously proposed multiresolution data fusion algorithm RubiX, resulting in a method for improved reconstruction, especially from data with lower number of diffusion gradients. We formulate the estimation of fiber parameters as a sparse signal recovery problem and propose a linear unmixing framework with sparse Bayesian learning for the recovery of sparse signals, the fiber orientations and volume fractions. The data is modeled using a parametric spherical deconvolution approach and represented using a dictionary created with the exponential decay components along different possible diffusion directions. Volume fractions of fibers along these directions define the dictionary weights. The proposed sparse inference, which is based on the dictionary representation, considers the sparsity of fiber populations and exploits the spatial redundancy in data representation, thereby facilitating inference from under-sampled q-space. The algorithm improves parameter estimation from dMRI through data-dependent local learning of hyperparameters, at each voxel and for each possible fiber orientation, that moderate the strength of priors governing the parameter variances. Experimental results on synthetic and in-vivo data show improved accuracy with a lower uncertainty in fiber parameter estimates. BusineX resolves a higher number of second and third fiber crossings. For under-sampled data, the algorithm is also shown to produce more reliable estimates.

We introduce Artifact-Based Rendering (ABR), a framework of tools, algorithms, and processes that makes it possible to produce real, data-driven 3D scientific visualizations with a visual language derived entirely from colors, lines, textures, and forms created using traditional physical media or found in nature. A theory and process for ABR is presented to address three current needs: (i) designing better visualizations by making it possible for non-programmers to rapidly design and critique many alternative data-to-visual mappings; (ii) expanding the visual vocabulary used in scientific visualizations to depict increasingly complex multivariate data; (iii) bringing a more engaging, natural, and human-relatable handcrafted aesthetic to data visualization. New tools and algorithms to support ABR include front-end applets for constructing artifact-based colormaps, optimizing 3D scanned meshes for use in data visualization, and synthesizing textures from artifacts. These are complemented by an interactive rendering engine with custom algorithms and interfaces that demonstrate multiple new visual styles for depicting point, line, surface, and volume data. A within-the-research-team design study provides early evidence of the shift in visualization design processes that ABR is believed to enable when compared to traditional scientific visualization systems. Qualitative user feedback on applications to climate science and brain imaging support the utility of ABR for scientific discovery and public communication., Published in IEEE VIS 2019, 9 pages of content with 2 pages of references, 12 figures

BACKGROUND: Spinocerebellar ataxia type 1 (SCA1) causes progressive degeneration of the cerebellum and brainstem. Volumetric magnetic resonance imaging (MRI) was shown to be more sensitive to disease progression than the most sensitive clinical measure, the Scale for the Assessment and Rating of Ataxia (SARA), in longitudinal studies, and magnetic resonance spectroscopy (MRS) was shown to detect neurochemical abnormalities with high sensitivity cross‐sectionally in SCA1. OBJECTIVES: The objectives of this study were to compare the sensitivities to change of volumetric MRI, MRS, and SARA in a 3‐year longitudinal study in SCA1. METHODS: A total of 16 early‐to‐moderate stage patients with SCA1 (SARA 0‐14) and 21 matched healthy participants were scanned up to 3 times with 1.5‐year intervals. Ataxia severity was assessed with SARA. T(1)‐weighted images and magnetic resonance spectra from the cerebellar vermis, cerebellar white matter, and pons were acquired at 3T. RESULTS: The pontine total N‐acetylaspartate‐to‐myo‐inositol ratio was the most sensitive MRS measure to change (−3.9 ± 4.6%/yr in SCA1 vs. −0.3 ± 3.5%/yr in controls; P

International audience; Computational methods are crucial for the analysis of diffusion magnetic resonance imaging (MRI) of the brain. Computational diffusion MRI can provide rich information at many size scales, including local microstructure measures such as diffusion anisotropies or apparent axon diameters, whole‐brain connectivity information that describes the brain's wiring diagram and population‐based studies in health and disease. Many of the diffusion MRI analyses performed today were not possible five, ten or twenty years ago, due to the requirements for large amounts of computer memory or processor time. In addition, mathematical frameworks had to be developed or adapted from other fields to create new ways to analyze diffusion MRI data. The purpose of this review is to highlight recent computational and statistical advances in diffusion MRI and to put these advances into context by comparison with the more traditional computational methods that are in popular clinical and scientific use. We aim to provide a high‐level overview of interest to diffusion MRI researchers, with a more in‐depth treatment to illustrate selected computational advances.

The adolescent major depressive disorder is one of the top 10 debilitating psychiatric illnesses and the effectiveness of current treatment methods are constrained by the limited understanding of biological causes. In this paper, we use diffusion tensor imaging to explore changes in anatomical connectivity between the MDD group (n=37) and control group (n=27). Furthermore, along-track analysis is performed to identify locations of alterations along connections with significant connectivity change. For the connection between the hippocampus and the insular cortex in the right hemisphere, decreased connectivity in axial diffusivity (AD), mean diffusivity (MD), and radial diffusivity (RD) was discovered. Additionally, the number of tracks for the same connection is increased for the MDD group. Moreover, for the connection between the parahippocampal gyrus and the insular cortex in the left hemisphere, increased connectivity is observed in fractional anisotropy (FA), AD, MD, and RD. Furthermore, the locations of significant alterations are identified to be between 65% to 100% from the insular cortex to the hippocampus in the right hemisphere and at the 80% location from the insular cortex to the parahippocampal gyrus in the left hemisphere. The significant and consistent white matter changes at the hippocampus end of the insula-hippocampus connection suggest potential correlations to the previously reported grey matter shrinkage and functional abnormalities.

Adolescence is a period of rapid brain maturation and a critical period for the onset of Major Depressive Disorder (MDD) that usually leads to serious outcomes such as suicide. Although changes in anatomical connectivity in MDD have been reported, changes in network topology for MDD remain unclear. Additionally, whether the changes are the same for adolescent MDD and adult MDD remains unclear as well. This paper explores anatomical features including: a) anatomical connectivity defined by diffusion tensor imaging measurements between a pair of brain regions, and b) topological measurements from anatomical networks, and apply machine learning approaches to identify responsive biomarkers distinguishing MDD patients from healthy subjects. In addition to statistical tests, univariate classifiers are designed to evaluate the discriminating power of features. Furthermore, multivariate classifiers are trained for distinguishing healthy subjects from MDD patients. The best classifier achieves an accuracy of 76.56%, 81.08% sensitivity, 70.37% specificity and 78.95% precision for 64 subjects (37 MDD and 27 matched healthy control). The selected features include: 1) betweenness centrality of the right lingual gyrus of the ADC network at 12% sparsity, 2) participation coefficient of the right pars opercularis of the AD network at 16% sparsity, 3) participation coefficient of the left insular cortex of the MD network at 21% sparsity, and 4) participation coefficient of the right lateral orbitofrontal cortex in the ADC network at 10% sparsity. These features reflect changes in the topological structure of the brain anatomical network in MDD.

Due to the high resistance (35%) to the current treatment methods in adolescent Major Depressive Disorder (MDD) and its tragic outcomes, the discovery of treatmentrelated responders is critical to developing effective treatments. In this paper, the permutation test is performed to identify statistically significant changes in anatomical characteristics during pairwise comparisons among the control group (n=27), treated MDD group (n=37), and untreated MDD group (n=15). The anatomical characteristics include: 1) anatomical connectivity defined using DTI metrics between a pair of brain regions, and 2) topological measurements of anatomical networks. With the Bonferroni correction for multiple-comparison, significant alterations in community structure and local topology were identified as the p-value < 5%, which include: 1) a reduced nodal centrality (degree and strength) on right hippocampus for treated compared to untreated group, 2) an elevated clustering coefficient and local efficiency on right lateral orbitofrontal cortex for untreated compared to the combination of control and treated groups, 3) an increased participation coefficient for untreated patients on left insula cortex in the meandiffusivity network compared to the combination of control and treated groups, and 4) a degraded module degree z-score on right caudate nucleus for all the patients compared to the control group. Two connections, hippocampus-insula in the right hemisphere and parahippocampal-insula in the left hemisphere, were found significantly altered in TR, AD, and FA due to MDD.

Adolescent Major Depressive Disorder (MDD) is a common and serious mental illness that could lead to tragic outcomes including chronic adult disability and suicide. In this paper, we explore anatomical features and apply machine learning approaches to identify responsive biomarkers distinguishing MDD patients from healthy subjects. The features of interest include metrics in two categories: a) anatomical connectivity defined by diffusion tensor imaging measurements between a pair of brain regions, and b) topological measurements from anatomical networks. A combination of p-value based filtering and minimum redundancy maximum relevance method is performed to select features for optimal classification accuracy. A leave-one-out cross-validation method is used for the classification performance evaluation. The proposed methodology achieves an improved accuracy of 78%, 90.39% sensitivity, and 79.66% precision for 79 subjects. The most distinguishing features are the betweenness centrality of the right lingual gyrus of the ADC network at 12% sparsity, the participation coefficient of the right lateral occipital sulcus of the ADC network at 22% sparsity, the participation coefficient of the right pars opercularis of the AD network at 16% sparsity, and the participation coefficient of the right lateral orbitofrontal cortex in the ADC network at 10% sparsity. Those network measures reflect the change of connectivity between the regions and their associated anatomical subnetworks.

Identification of the treatment-related responders for adolescent Major Depressive Disorder (MDD) is urgently needed to develop effective treatments. In this paper, machine learning based classifiers are used to reveal anatomical features as responders for distinguishing MDD patients who have received treatment from those who never received any treatment. The features are drawn from two sets of measurements: 1) anatomical connectivity defined by diffusion tensor imaging measurements between a pair of brain regions, and 2) topological measurements from anatomical networks. Feature selection was performed based on p-value and minimum redundancy maximum relevance (mRMR) method to achieve improved classification accuracy. The classification performance is evaluated with a leave-one-out cross-validation method using 37 treated and 15 untreated subjects. The proposed methodology achieves 73% accuracy, 100% specificity, and 100% precision for 52 subjects. The most distinguishing features are the strength of the right hippocampus of the mean diffusivity (MD) network at 18% density and of the track-count (TR) network, the participation coefficient of the left middle temporal gyrus of the radial diffusivity (RD) network at 20% density, the axial diffusivity (AD) connectivity between right middle temporal gyrus and right supramarginal gyrus, the betweenness centrality of the right hippocampus of the TR network at 11% density, the apparent diffusion coefficient (ADC) connectivity between the left pars opercularis and the left rostral anterior cingulate cortex, the clustering coefficient of the middle anterior corpus callosum of the TR network at 11% density, and the AD connectivity between the left pars opercularis and the left rostral anterior cingulate cortex.

Mapping structural connectivity in healthy adults for the Human Connectome Project (HCP) benefits from high quality, high resolution, multiband (MB)-accelerated whole brain diffusion MRI (dMRI). Acquiring such data at ultrahigh fields (7 T and above) can improve intrinsic signal-to-noise ratio (SNR), but suffers from shorter T2 and T2* relaxation times, increased B1+ inhomogeneity (resulting in signal loss in cerebellar and temporal lobe regions), and increased power deposition (i.e. Specific Absorption Rate (SAR)), thereby limiting our ability to reduce the repetition time (TR). Here, we present recent developments and optimizations in 7 T image acquisitions for the HCP that allow us to efficiently obtain high-quality, high-resolution whole brain in-vivo dMRI data at 7 T. These data show spatial details typically seen only in ex-vivo studies and complement already very high quality 3 T HCP data in the same subjects. The advances are the result of intensive pilot studies aimed at mitigating the limitations of dMRI at 7 T. The data quality and methods described here are representative of the datasets that will be made freely available to the community in 2015.

We propose a sparse Bayesian learning algorithm for improved estimation of white matter fiber parameters from compressed (under-sampled q-space) multi-shell diffusion MRI data. The multi-shell data is represented in a dictionary form using a non-monoexponential decay model of diffusion, based on continuous gamma distribution of diffusivities. The fiber volume fractions with predefined orientations, which are the unknown parameters, form the dictionary weights. These unknown parameters are estimated with a linear un-mixing framework, using a sparse Bayesian learning algorithm. A localized learning of hyperparameters at each voxel and for each possible fiber orientations improves the parameter estimation. Our experiments using synthetic data from the ISBI 2012 HARDI reconstruction challenge and in-vivo data from the Human Connectome Project demonstrate the improvements.

Although brain functionality is often remarkably robust to lesions and other insults, it may be fragile when these take place in specific locations. Previous attempts to quantify robustness and fragility sought to understand how the functional connectivity of brain networks is affected by structural changes, using either model-based predictions or empirical studies of the effects of lesions. We advance a geometric viewpoint relying on a notion of network curvature, the so-called Ollivier-Ricci curvature. This approach has been proposed to assess financial market robustness and to differentiate biological networks of cancer cells from healthy ones. Here, we apply curvature-based measures to brain structural networks to identify robust and fragile brain regions in healthy subjects. We show that curvature can also be used to track changes in brain connectivity related to age and autism spectrum disorder (ASD), and we obtain results that are in agreement with previous MRI studies., The brain can often continue to function despite lesions in many areas, but damage to particular locations may have serious effects. Here, the authors use the concept of Ollivier-Ricci curvature to investigate the robustness of brain networks.

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. Magn Reson Med 72:1471–1485, 2014. © 2013 Wiley Periodicals, Inc.

Purpose: To examine the effects of the reconstruction algorithm of magnitude images from multi-channel diffusion MRI on fibre orientation estimation. Theory and Methods: It is well established that the method used to combine signals from different coil elements in multi-channel MRI can have an impact on the properties of the reconstructed magnitude image. Utilising a root-sum-of-squares (RSoS) approach results in a magnitude signal that follows an effective non-central-distribution. As a result, the noise floor, the minimum measurable in the absence of any true signal, is elevated. This is particularly relevant for diffusion-weighted MRI, where the signal attenuation is of interest. Results: In this study, we illustrate problems that such image reconstruction characteristics may cause in the estimation of fibre orientations, both for model-based and model-free approaches, when modern 32-channel coils are employed. We further propose an alternative image reconstruction method that is based on sensitivity encoding (SENSE) and preserves the Rician nature of the single-channel, magnitude MR signal. We show that for the same k-space data, RSoS can cause excessive overfitting and reduced precision in orientation estimation compared to the SENSE-based approach. Conclusion: These results highlight the importance of choosing the appropriate image reconstruction method for tractography studies that use multi-channel receiver coils for diffusion MRI acquisition.

The acquisition of high angular resolution diffusion MRI is particularly long and subject motion can become an issue. The orientation distribution function (ODF) can be reconstructed online incrementally from diffusion-weighted MRI with a Kalman filtering framework. This online reconstruction provides real-time feedback throughout the acquisition process. In this article, the Kalman filter is first adapted to the reconstruction of the ODF in constant solid angle. Then, a method called STAR (STatistical Analysis of Residuals) is presented and applied to the online detection of motion in high angular resolution diffusion images. Compared to existing techniques, this method is image based and is built on top of a Kalman filter. Therefore, it introduces no additional scan time and does not require additional hardware. The performance of STAR is tested on simulated and real data and compared to the classical generalized likelihood ratio test. Successful detection of small motion is reported (rotation under 2°) with no delay and robustness to noise., National Institutes of Health (U.S.) (NIH grant Grant P41 RR008079), National Institutes of Health (U.S.) (NIH grant P41 EB015894), National Institutes of Health (U.S.) (NIH grant P30 NS057091), National Institutes of Health (U.S.) (Human Connectome Project U54 MH091657), United States. Air Force Office of Scientific Research (NSSEFF), National Science Foundation (U.S.), United States. Army Research Office, United States. Defense Advanced Research Projects Agency, United States. National Geospatial-Intelligence Agency, France. Agence nationale de la recherche (ANR NucleiPark), Institut national de recherche en informatique et en automatique (France), National Institutes of Health (U.S.) (Human Connectome Project, Grant R01 EB008432)

Hardware constraints, scanning time limitations, patient movement, and signal-to-noise ratio (SNR) considerations, restrict the slice-selection and the in-plane resolutions of MRI differently, generally resulting in anisotropic voxels. This nonuniform sampling can be problematic, especially in image segmentation and clinical examination. To alleviate this, the acquisition is divided into (two or) three separate scans, with higher in-plane resolutions and thick slices, yet orthogonal slice-selection directions. In this work, a noniterative wavelet-based approach for combining the three orthogonal scans is adopted, and its advantages compared with other existing methods, such as Fourier techniques, are discussed, including the consideration of the actual pulse response of the MRI scanner, and its lower computational complexity. Experimental results are shown on simulated and real 7 T MRI data.

High angular resolution diffusion imaging (HARDI) has become an important technique for imaging complex oriented structures in the brain and other anatomical tissues. This has motivated the recent development of several methods for computing the orientation probability density function (PDF) at each voxel. However, much less work has been done on developing techniques for filtering, interpolation, averaging and principal geodesic analysis of orientation PDF fields. In this paper, we present a Riemannian framework for performing such operations. The proposed framework does not require that the orientation PDFs be represented by any fixed parameterization, such as a mixture of von Mises-Fisher distributions or a spherical harmonic expansion. Instead, we use a nonparametric representation of the orientation PDF. We exploit the fact that under the square-root re-parameterization, the space of orientation PDFs forms a Riemannian manifold: the positive orthant of the unit Hilbert sphere. We show that various orientation PDF processing operations, such as filtering, interpolation, averaging and principal geodesic analysis, may be posed as optimization problems on the Hilbert sphere, and can be solved using Riemannian gradient descent. We illustrate these concepts with numerous experiments on synthetic, phantom and real datasets. We show their application to studying left/right brain asymmetries.

This paper presents a coupled thermo-hydro-mechanical (THM) model enriched with a buckling-type criterion for progressive spalling. In the first part of the paper, a general fully coupled multi-phase THM model describing the behaviour of concrete at moderate and high temperatures is presented. Then the spalling criterion and its numerical implementation in the framework of the finite element method are presented. Finally, a simple 1D numerical example will illustrate the effectiveness of the implemented numerical approach.

Estimating brain wiring patterns is critical to better understand the brain organization and function. Anatomical brain connectivity models axonal pathways, while the functional brain connectivity characterizes the statistical dependencies and correlation between the activities of various brain regions. The synchronization of brain activity can be inferred through the variation of blood-oxygen-level dependent (BOLD) signal from functional MRI (fMRI) and the neural connections can be estimated using tractography from diffusion MRI (dMRI). Functional connections between brain regions are supported by anatomical connections, and the synchronization of brain activities arises through sharing of information in the form of electro-chemical signals on axon pathways. Jointly modeling fMRI and dMRI data may improve the accuracy in constructing anatomical connectivity as well as functional connectivity. Such an approach may lead to novel multimodal biomarkers potentially able to better capture functional and anatomical connectivity variations. We present a novel brain network model which jointly models the dMRI and fMRI data to improve the anatomical connectivity estimation and extract the anatomical subnetworks associated with specific functional modes by constraining the anatomical connections as structural supports to the functional connections. The key idea is similar to a multi-commodity flow optimization problem that minimizes the cost or maximizes the efficiency for flow configuration and simultaneously fulfills the supply-demand constraint for each commodity. In the proposed network, the nodes represent the grey matter (GM) regions providing brain functionality, and the links represent white matter (WM) fiber bundles connecting those regions and delivering information. The commodities can be thought of as the information corresponding to brain activity patterns as obtained for instance by independent component analysis (ICA) of fMRI data. The concept of information flow is introduced and used to model the propagation of information between GM areas through WM fiber bundles. The link capacity , i.e., ability to transfer information, is characterized by the relative strength of fiber bundles, e.g., fiber count gathered from the tractography of dMRI data. The node information demand is considered to be proportional to the correlation between neural activity at various cortical areas involved in a particular functional mode (e.g. visual, motor, etc.). These two properties lead to the link capacity and node demand constraints in the proposed model. Moreover, the information flow of a link cannot exceed the demand from either end node. This is captured by the feasibility constraints . Two different cost functions are considered in the optimization formulation in this paper. The first cost function, the reciprocal of fiber strength represents the unit cost for information passing through the link. In the second cost function, a min-max (minimizing the maximal link load) approach is used to balance the usage of each link. Optimizing the first cost function selects the pathway with strongest fiber strength for information propagation. In the second case, the optimization procedure finds all the possible propagation pathways and allocates the flow proportionally to their strength. Additionally, a penalty term is incorporated with both the cost functions to capture the possible missing and weak anatomical connections. With this set of constraints and the proposed cost functions, solving the network optimization problem recovers missing and weak anatomical connections supported by the functional information and provides the functional-associated anatomical subnetworks. Feasibility is demonstrated using realistic diffusion and functional MRI phantom data. It is shown that the proposed model recovers the maximum number of true connections, with fewest number of false connections when compared with the connectivity derived from a joint probabilistic model using the expectation-maximization (EM) algorithm presented in a prior work. We also apply the proposed method to data provided by the Human Connectome Project (HCP).

The RubiX [1] algorithm combines high SNR characteristics of low resolution data with high spacial specificity of high resolution data, to extract microstructural tissue parameters from diffusion MRI. In this paper we focus on estimating crossing fiber orientations and introduce sparsity to the RubiX algorithm, making it suitable for reconstruction from compressed (under-sampled) data. We propose a sparse Bayesian algorithm for estimation of fiber orientations and volume fractions from compressed diffusion MRI. The data at high resolution is modeled using a parametric spherical deconvolution approach and represented using a dictionary created with the exponential decay components along different possible directions. Volume fractions of fibers along these orientations define the dictionary weights. The data at low resolution is modeled using a spatial partial volume representation. The proposed dictionary representation and sparsity priors consider the dependence between fiber orientations and the spatial redundancy in data representation. Our method exploits the sparsity of fiber orientations, therefore facilitating inference from under-sampled data. Experimental results show improved accuracy and decreased uncertainty in fiber orientation estimates. For under-sampled data, the proposed method is also shown to produce more robust estimates of fiber orientations.

Quantitative investigations of fiber orientation and structural connectivity at microscopic resolution have led to great challenges for current neuroimaging techniques. Here, we present a structure tensor (ST) analysis of ex vivo rat brain images acquired by a multicontrast (MC) serial optical coherence scanner. The ST considers the gradients of images in local neighbors to generate a matrix whose eigen-decomposition can estimate the local features such as the edges, anisotropy, and orientation of tissue constituents. This computational analysis is applied on the conventional- and polarization-based contrasts of optical coherence tomography. The three-dimensional (3-D) fiber orientation maps are computed from the image stacks of sequential scans both at mesoresolution for a global view and at high-resolution for the details. The computational orientation maps demonstrate a good agreement with the optic axis orientation contrast which measures the in-plane fiber orientation. Moreover, tractography is implemented using the directional information extracted from the 3-D ST. The study provides a unique opportunity to leverage MC high-resolution information to map structural connectivity of the whole brain.

Volumetric segmentation of subcortical structures, such as the basal ganglia and thalamus, is necessary for noninvasive diagnosis and neurosurgery planning. This is a challenging problem due in part to limited boundary information between structures, similar intensity profiles across the different structures, and low contrast data. This paper presents a semiautomatic segmentation system exploiting the superior image quality of ultrahigh field (7 T) MRI. The proposed approach utilizes the complementary edge information in the multiple structural MRI modalities. It combines optimally selected two modalities from susceptibility-weighted, T2-weighted, and diffusion MRI, and introduces a tailored new edge indicator function. In addition to this, we employ prior shape and configuration knowledge of the subcortical structures in order to guide the evolution of geometric active surfaces. Neighboring structures are segmented iteratively, constraining oversegmentation at their borders with a nonoverlapping penalty. Several experiments with data acquired on a 7 T MRI scanner demonstrate the feasibility and power of the approach for the segmentation of basal ganglia components critical for neurosurgery applications such as deep brain stimulation surgery.

We introduce a framework for population analysis of white matter tracts based on diffusion-weighted images of the brain. The framework enables extraction of fibers from high angular resolution diffusion images (HARDI); clustering of the fibers based partly on prior knowledge from an atlas; representation of the fiber bundles compactly using a path following points of highest density (maximum density path; MDP); and registration of these paths together using geodesic curve matching to find local correspondences across a population. We demonstrate our method on 4-Tesla HARDI scans from 565 young adults to compute localized statistics across 50 white matter tracts based on fractional anisotropy (FA). Experimental results show increased sensitivity in the determination of genetic influences on principal fiber tracts compared to the tract-based spatial statistics (TBSS) method. Our results show that the MDP representation reveals important parts of the white matter structure and considerably reduces the dimensionality over comparable fiber matching approaches.

We consider the problem of improving the accuracy and reliability of probabilistic white matter tractography methods by improving the built-in sampling scheme, which randomly draws, from a diffusion model such as the orientation distribution function (ODF), a direction of propagation. Existing methods employing inverse transform sampling require an ad hoc thresholding step to prevent the less likely directions from being sampled. We herein propose to perform importance sampling of spherical harmonics, which redistributes an input point set on the sphere to match the ODF using hierarchical sample warping. This produces a point set that is more concentrated around the modes, allowing the subsequent inverse transform sampling to generate orientations that are in better accordance with the local fiber configuration. Integrated into a Kalman filter-based framework, our approach is evaluated through experiments on synthetic, phantom, and real datasets.

Registration of diffusion-weighted magnetic resonance images (DW-MRIs) is a key step for population studies, or construction of brain atlases, among other important tasks. Given the high dimensionality of the data, registration is usually performed by relying on scalar representative images, such as the fractional anisotropy (FA) and non diffusion-weighted (b0) images, thereby ignoring much of the directional information conveyed by DW-MR datasets itself. Alternatively, model-based registration algorithms have been proposed to exploit information on the preferred fiber orientation(s) at each voxel. Models such as the diffusion tensor or orientation distribution function (ODF) have been used for this purpose. Tensor-based registration methods rely on a model that does not completely capture the information contained in DW-MRIs, and largely depends on the accurate estimation of tensors. ODF-based approaches are more recent and computationally challenging, but also better describe complex fiber configurations thereby potentially improving the accuracy of DW-MRI registration. A new algorithm based on angular interpolation of the diffusion-weighted volumes was proposed for affine registration, and does not rely on any specific local diffusion model. In this work, we first extensively compare the performance of registration algorithms based on (i) angular interpolation, (ii) non diffusion-weighted scalar volume (b0), and (iii) diffusion tensor image (DTI). Moreover, we generalize the concept of angular interpolation (AI) to non-linear image registration, and implement it in the FMRIB Software Library (FSL). We demonstrate that AI registration of DW-MRIs is a powerful alternative to volume and tensor-based approaches. In particular, we show that AI improves the registration accuracy in many cases over existing state-of-the-art algorithms, while providing registered raw DW-MRI data, which can be used for any subsequent analysis.

We consider the problem of tracking white matter fibers in high angular resolution diffusion imaging (HARDI) data while simultaneously estimating the local fiber orientation profile. Prior work showed that an unscented Kalman filter (UKF) can be used for this problem, yet existing algorithms employ parametric mixture models to represent water diffusion and to define the state space. To address this restrictive model dependency, we propose to extend the UKF to HARDI data modeled by orientation distribution functions (ODFs), a more generic diffusion model. We consider the spherical harmonic representation of the HARDI signal as the state, enforce nonnegativity of the ODFs, and perform tractography using the directions at which the ODFs attain their peaks. In simulations, our method outperforms filtered two-tensor tractography at different levels of noise by achieving a reduction in mean Chamfer error of 0.05 to 0.27 voxels; it also produced in vivo fiber tracking that is consistent with the neuroanatomy.

The presence of noise in High Angular Resolution Diffusion Imaging (HARDI) data of the brain can limit the accuracy with which fiber pathways of the brain can be extracted. In this work, we present a variational model to denoise HARDI data corrupted by Rician noise. Numerical experiments are performed on three types of data: 2D synthetic data, 3D diffusion-weighted Magnetic Resonance Imaging (DW-MRI) data of a hardware phantom containing synthetic fibers, and 3D real HARDI brain data. Experiments show that our model is effective for denoising HARDI-type data while preserving important aspects of the fiber pathways such as fractional anisotropy and the orientation distribution functions.

Cortical connectivity is associated with cognitive and behavioral traits that are thought to vary between sexes. Using high-angular resolution diffusion imaging at 4 Tesla, we scanned 234 young adult twins and siblings (mean age: 23.4 ± 2.0 SD years) with 94 diffusion-encoding directions. We applied a novel Hough transform method to extract fiber tracts throughout the entire brain, based on fields of constant solid angle orientation distribution functions (ODFs). Cortical surfaces were generated from each subject's 3D T1-weighted structural MRI scan, and tracts were aligned to the anatomy. Network analysis revealed the proportions of fibers interconnecting 5 key subregions of the frontal cortex, including connections between hemispheres. We found significant sex differences (147 women/87 men) in the proportions of fibers connecting contralateral superior frontal cortices. Interhemispheric connectivity was greater in women, in line with long-standing theories of hemispheric specialization. These findings may be relevant for ongoing studies of the human connectome.

International audience; The orientation distribution function (ODF) can be reconstructed online incrementally from diffusion-weighted MRI with a Kalman filtering framework. This online reconstruction can provide real-time feedback to the practitioner, especially appreciated for long acquisition protocols typical in Q-ball imaging. On top of the Kalman filter, we propose a method to evaluate online the reconstruction accuracy of the estimated ODF in constant solid angle. In addition, monitoring the residuals of the Kalman filter, we design, based on statistical tests, two algorithms for online detection of subject motion. The proposed techniques, tested on real and synthetic data under various experimental conditions, can detect rotation by angle less than 3°.

High angular resolution diffusion imaging (HARDI) allows in vivo analysis of the white matter structure and connectivity. Based on orientation distribution functions (ODFs) that represent the directionality of water diffusion at each point in the brain, tractography methods can recover major axonal pathways. This enables tract-based analysis of fiber integrity and connectivity. For multi-subject comparisons, fibers may be clustered into bundles that are consistently found across subjects. To do this, we scanned 20 young adults with HARDI at 4 T. From the reconstructed ODFs, we performed whole-brain tractography with a novel Hough transform method. We then used measures of agreement between the extracted 3D curves and a co-registered probabilistic DTI atlas to select key pathways. Using median filtering and a shortest path graph search, we derived the maximum density path to compactly represent each tract in the population. With this tract-based method, we performed tract-based analysis of fractional anisotropy, and assessed how the chosen tractography algorithm influenced the results. The resulting method may expedite population-based statistical analysis of HARDI and DTI.

q-Ball imaging is a high-angular-resolution diffusion imaging technique that has been proven very successful in resolving multiple intravoxel fiber orientations in MR images. The standard computation of the orientation distribution function (the probability of diffusion in a given direction) from q-ball data uses linear radial projection, neglecting the change in the volume element along each direction. This results in spherical distributions that are different from the true orientation distribution functions. For instance, they are neither normalized nor as sharp as expected and generally require postprocessing, such as artificial sharpening. In this paper, a new technique is proposed that, by considering the solid angle factor, uses the mathematically correct definition of the orientation distribution function and results in a dimensionless and normalized orientation distribution function expression. Our model is flexible enough so that orientation distribution functions can be estimated either from single q-shell datasets or by exploiting the greater information available from multiple q-shell acquisitions. We show that the latter can be achieved by using a more accurate multiexponential model for the diffusion signal. The improved performance of the proposed method is demonstrated on artificial examples and high-angular-resolution diffusion imaging data acquired on a 7-T magnet.

International audience; The diffusion orientation distribution function (ODF) can be reconstructed from q-ball imaging (QBI) to map the complex intravoxel structure of water diffusion. As acquisition time is particularly large for high angular resolution diffusion imaging (HARDI), fast estimation algorithms have recently been proposed, as an on-line feedback on the reconstruction accuracy. Thus the acquisition could be stopped or continued on demand. We adapt these real-time algorithms to the mathematically correct definition of ODF in constant solid angle (CSA), and develop a motion detection algorithm upon this reconstruction. Results of improved fiber crossing detection by CSA ODF are shown, and motion detection was implemented and tested in vivo.

Localizing brain structures and pathways associated with intellectual performance and cognition may clarify the networks involved, their development, and their dysfunction in disease. Here we introduce a measure of fiber crossing computed from high angular resolution diffusion images (HARDI), with the goal of better characterizing complex white matter microstructure. We aimed to detect white matter features associated with intellectual performance, previously missed due to limitations of the single tensor diffusion model.