Spherical Deconvolution Research Articles

VRfRNet: Volumetric ROI fODF reconstruction network for estimation of multi-tissue constrained spherical deconvolution with only single shell dMRI

Diffusion MRI (dMRI) is one of the most popular techniques for studying the brain structure, mainly the white matter region. Among several sampling methods in dMRI, the high angular resolution diffusion imaging (HARDI) technique has attracted researchers due to its more accurate fiber orientation estimation. However, the current single-shell HARDI makes the intravoxel structure challenging to estimate accurately. While multi-shell acquisition can address this problem, it takes a longer scanning time, restricting its use in clinical applications. In addition, most existing dMRI scanners with low gradient-strengths often acquire single-shell up to b=1000s/mm2 because of signal-to-noise ratio issues and severe image artefacts. Hence, we propose a novel generative adversarial network, VRfRNet, for the reconstruction of multi-shell multi-tissue fiber orientation distribution function from single-shell HARDI volumes. Such a transformation learning is performed in the spherical harmonics (SH) space, as raw input HARDI volume is transformed to SH coefficients to soften gradient directions. The proposed VRfRNet consists of several modules, such as multi-context feature enrichment module, feature level attention, and softmax level attention. In addition, three loss functions have been used to optimize network learning, including L1, adversarial, and total variation. The network is trained and tested using standard qualitative and quantitative performance metrics on the publicly available HCP data-set.

Topological Alteration Is Associated with Non-dependent Alcohol Use in Bipolar Disorder.

Introduction: Structural alterations in cortical thickness and the microstructural organization of white matter are independently associated with non-dependent alcohol consumption and bipolar disorder (BD). Identifying their interactive and network-level effects on brain topology may identify the impact of alcohol on reward and emotion circuitry, and its contribution to relapse in BD. Materials and Methods: Thirty-four BD-I (DSM-IV-TR) and 38 healthy controls (HC) underwent T1 and diffusion-weighted magnetic resonance imaging scanning, and the Alcohol Use Disorders Identification Test-Consumption to assess alcohol use. Connectomes comprising 34 cortical and 9 subcortical nodes bilaterally (Freesurfer v5.3) connected by fractional anisotropy-weighted edges derived from non-tensor based deterministic constrained spherical deconvolution tractography (ExploreDTI v4.8.6) underwent permutation-based topological analysis (NBS v1.2) and were examined for the effects of alcohol use and diagnosis-by-alcohol use accounting for age, sex, and diagnosis. Results: Alcohol was significantly related to a subnetwork, encompassing connections between fronto-limbic, basal ganglia, and temporal nodes (Frange = 5-8.4, p = 0.031) and it was not detected to have an effect on global brain integration or segregation. A portion of this network (18%), involving cortico-limbic and basal ganglia connections, was differentially impacted by alcohol in the BD relative to the control group (Frange = 5-8.8, p = 0.033), despite the groups' consuming similar amounts of alcohol (BD: mean ± standard deviation 4.95 ± 3.0; HC 3.62 ± 3.0, T = 1.88, p = 0.06). Discussion: Non-dependent alcohol use impacts brain architectural organization and connectivity within salience, reward, and affective circuitry. The relationship between alcohol use and topology of the network in BD suggests an interactive effect between specific biological vulnerability and alcohol use, which may explain the susceptibility to an increased risk of relapse in the disorder. Impact statement The association between non-dependent alcohol use and neural architecture in bipolar disorder (BD) is unknown, despite the poor clinical trajectory and increased likelihood of relapse associated with alcohol use in BD. We demonstrate that together alcohol and a diagnosis of BD is associated with a subnetwork involving nodes of the cortico-limbic and reward networks. This subnetwork, demonstrated in BD and absent in controls, differentially involves nodes that are specific to reward and emotion processes. This suggests a diagnosis-specific biological vulnerability for alcohol use and may be consistent with known mood lability and thus relapse associated with alcohol use in BD.

Accuracy and reliability of diffusion imaging models

Diffusion imaging aims to non-invasively characterize the anatomy and integrity of the brain's white matter fibers. We evaluated the accuracy and reliability of commonly used diffusion imaging methods as a function of data quantity and analysis method, using both simulations and highly sampled individual-specific data (927–1442 diffusion weighted images [DWIs] per individual). Diffusion imaging methods that allow for crossing fibers (FSL's BedpostX [BPX], DSI Studio's Constant Solid Angle Q-Ball Imaging [CSA-QBI], MRtrix3’s Constrained Spherical Deconvolution [CSD]) estimated excess fibers when insufficient data were present and/or when the data did not match the model priors. To reduce such overfitting, we developed a novel Bayesian Multi-tensor Model-selection (BaMM) method and applied it to the popular ball-and-stick model used in BedpostX within the FSL software package. BaMM was robust to overfitting and showed high reliability and the relatively best crossing-fiber accuracy with increasing amounts of diffusion data. Thus, sufficient data and an overfitting resistant analysis method enhance precision diffusion imaging. For potential clinical applications of diffusion imaging, such as neurosurgical planning and deep brain stimulation (DBS), the quantities of data required to achieve diffusion imaging reliability are lower than those needed for functional MRI.

Longitudinal tracking of axonal loss using diffusion magnetic resonance imaging in multiple sclerosis.

Axonal loss in the CNS is a key driver of progressive neurological impairments in people with multiple sclerosis. Currently, there are no established methods for tracking axonal loss clinically. This study aimed to determine the sensitivity of longitudinal diffusion MRI-derived fibre-specific measures of axonal loss in people with multiple sclerosis. Fibre measures were derived from diffusion MRI acquired as part of a standard radiological MRI protocol and were compared (i) to establish measures of neuro-axonal degeneration: brain parenchymal fraction and retinal nerve fibre layer thickness and (ii) between different disease stages: clinically isolated syndrome and early/late relapsing–remitting multiple sclerosis. Retrospectively identified data from 59 people with multiple sclerosis (18 clinically isolated syndrome, 22 early and 19 late relapsing–remitting) who underwent diffusion MRI as part of their routine clinical monitoring were collated and analysed. Twenty-six patients had 1-year and 14 patients had a 2-year follow-up. Brain parenchymal fraction was calculated from 3D MRI scans, and fibre-specific measures were calculated from diffusion MRI using multi-tissue constrained spherical deconvolution. At each study visit, patients underwent optical coherence tomography to determine retinal nerve fibre layer thickness, and standard neurological assessment expanded the disability status scale. We found a significant annual fibre-specific neuro-axonal degeneration (mean ± SD = −3.49 ± 3.32%, P < 0.001) that was ∼7 times larger than the annual change of brain parenchymal fraction (−0.53 ± 0.95%, P < 0.001), and more than four times larger than annual retinal nerve fibre layer thinning (−0.75 ± 2.50% P = 0.036). Only fibre-specific measures showed a significant difference in annual degeneration between the disease stages (P = 0.029). Reduced brain parenchymal fraction, retinal nerve fibre layer thickness and fibre-specific measures were moderately related to higher expanded disability status scale (rho = −0.368, rho = −0.408 and rho = −0.365, respectively). Fibre-specific measures can be measured from data collected within a standard radiological multiple sclerosis study and are substantially more sensitive to longitudinal change compared with brain atrophy and retinal nerve fibre layer thinning.

An atlas of white matter anatomy, its variability, and reproducibility based on constrained spherical deconvolution of diffusion MRI

Virtual dissection of white matter (WM) using diffusion MRI tractography is confounded by its poor reproducibility. Despite the increased adoption of advanced reconstruction models, early region-of-interest driven protocols based on diffusion tensor imaging (DTI) remain the dominant reference for virtual dissection protocols. Here we bridge this gap by providing a comprehensive description of typical WM anatomy reconstructed using a reproducible automated subject-specific parcellation-based approach based on probabilistic constrained-spherical deconvolution (CSD) tractography. We complement this with a WM template in MNI space comprising 68 bundles, including all associated anatomical tract selection labels and associated automated workflows. Additionally, we demonstrate bundle inter- and intra-subject variability using 40 (20 test-retest) datasets from the human connectome project (HCP) and 5 sessions with varying b-values and number of b-shells from the single-subject Multiple Acquisitions for Standardization of Structural Imaging Validation and Evaluation (MASSIVE) dataset. The most reliably reconstructed bundles were the whole pyramidal tracts, primary corticospinal tracts, whole superior longitudinal fasciculi, frontal, parietal and occipital segments of the corpus callosum and middle cerebellar peduncles. More variability was found in less dense bundles, e.g., the fornix, dentato-rubro-thalamic tract (DRTT), and premotor pyramidal tract. Using the DRTT as an example, we show that this variability can be reduced by using a higher number of seeding attempts. Overall inter-session similarity was high for HCP test-retest data (median weighted-dice = 0.963, stdev = 0.201 and IQR = 0.099). Compared to the HCP-template bundles there was a high level of agreement for the HCP test-retest data (median weighted-dice = 0.747, stdev = 0.220 and IQR = 0.277) and for the MASSIVE data (median weighted-dice = 0.767, stdev = 0.255 and IQR = 0.338). In summary, this WM atlas provides an overview of the capabilities and limitations of automated subject-specific probabilistic CSD tractography for mapping white matter fasciculi in healthy adults. It will be most useful in applications requiring a reproducible parcellation-based dissection protocol, and as an educational resource for applied neuroimaging and clinical professionals.

Avoiding data loss: Synthetic MRIs generated from diffusion imaging can replace corrupted structural acquisitions for freesurfer-seeded tractography.

Magnetic Resonance Imaging (MRI) motion artefacts frequently complicate structural and diffusion MRI analyses. While diffusion imaging is easily 'scrubbed' of motion affected volumes, the same is not true for T1w or T2w 'structural' images. Structural images are critical to most diffusion-imaging pipelines thus their corruption can lead to disproportionate data loss. To enable diffusion-image processing when structural images are missing or have been corrupted, we propose a means by which synthetic structural images can be generated from diffusion MRI. This technique combines multi-tissue constrained spherical deconvolution, which is central to many existing diffusion analyses, with the Bloch equations that allow simulation of MRI intensities for given scanner parameters and magnetic resonance (MR) tissue properties. We applied this technique to 32 scans, including those acquired on different scanners, with different protocols and with pathology present. The resulting synthetic T1w and T2w images were visually convincing and exhibited similar tissue contrast to acquired structural images. These were also of sufficient quality to drive a Freesurfer-based tractographic analysis. In this analysis, probabilistic tractography connecting the thalamus to the primary sensorimotor cortex was delineated with Freesurfer, using either real or synthetic structural images. Tractography for real and synthetic conditions was largely identical in terms of both voxels encountered (Dice 0.88-0.95) and mean fractional anisotropy (intrasubject absolute difference 0.00-0.02). We provide executables for the proposed technique in the hope that these may aid the community in analysing datasets where structural image corruption is common, such as studies of children or cognitively impaired persons.

Automatic oculomotor nerve identification based on data-driven fiber clustering.

The oculomotor nerve (OCN) is the main motor nerve innervating eye muscles and can be involved in multiple flammatory, compressive, or pathologies. The diffusion magnetic resonance imaging (dMRI) tractography is now widely used to describe the trajectory of the OCN. However, the complex cranial structure leads to difficulties in fiber orientation distribution (FOD) modeling, fiber tracking, and region of interest (ROI) selection. Currently, the identification of OCN relies on expert manual operation, resulting in challenges, such as the carries high clinical, time‐consuming, and labor costs. Thus, we propose a method that can automatically identify OCN from dMRI tractography. First, we choose the multi‐shell multi‐tissue constraint spherical deconvolution (MSMT‐CSD) FOD estimation model and deterministic tractography to describe the 3D trajectory of the OCN. Then, we rely on the well‐established computational pipeline and anatomical expertise to create a data‐driven OCN tractography atlas from 40 HCP data. We identify six clusters belonging to the OCN from the atlas, including the structures of three kinds of positional relationships (pass between, pass through, and go around) with the red nuclei and two kinds of positional relationships with medial longitudinal fasciculus. Finally, we apply the proposed OCN atlas to identify the OCN automatically from 40 new HCP subjects and two patients with brainstem cavernous malformation. In terms of spatial overlap and visualization, experiment results show that the automatically and manually identified OCN fibers are consistent. Our proposed OCN atlas provides an effective tool for identifying OCN by avoiding the traditional selection strategy of ROIs.

Brain microstructural antecedents of visual difficulties in infants born very preterm

Infants born very preterm (VPT) are at risk of later visual problems. Although neonatal screening can identify ophthalmologic abnormalities, subtle perinatal brain injury and/or delayed brain maturation may be significant contributors to complex visual-behavioral problems. Our aim was to assess the micro and macrostructural antecedents of early visual-behavioral difficulties in VPT infants by using diffusion MRI (dMRI) at term-equivalent age.We prospectively recruited a cohort of 262 VPT infants (≤32 weeks gestational age [GA]) from five neonatal intensive care units. We obtained structural and diffusion MRI at term-equivalent age and administered the Preverbal Visual Assessment (PreViAs) questionnaire to parents at 3–4 months corrected age. We used constrained spherical deconvolution to reconstruct nine white matter tracts of the visual pathways with high reliability and performed fixel-based analysis to derive fiber density (FD), fiber-bundle cross-section (FC), and combined fiber density and cross-section (FDC). In multiple logistic regression analyses, we related these tract metrics to visual-behavioral function.Of 262 infants, 191 had both high-quality dMRI and completed PreViAs, constituting the final cohort: mean (SD) GA was 29.3 (2.4) weeks, 90 (47.1%) were males, and postmenstrual age (PMA) at MRI was 42.8 (1.3) weeks. FD and FC of several tracts were altered in infants with (N = 59) versus those without retinopathy of prematurity (N = 132). FDC of the left posterior thalamic radiations (PTR), left inferior longitudinal fasciculus (ILF), right superior longitudinal fasciculus (SLF), and left inferior fronto-occipital fasciculus (IFOF) were significantly associated with visual attention scores, prior to adjusting for confounders. After adjustment for PMA at MRI, GA, severe retinopathy of prematurity, and total brain volume, FDC of the left PTR, left ILF, and left IFOF remained significantly associated with visual attention.Early visual-behavioral difficulties in VPT infants are preceded by micro and macrostructural abnormalities in several major visual pathways at term-equivalent age.

Insular gliomas and tractographic visualization of the connectome.

In this video, the authors present a connectome-guided surgical resection of an insular glioma in a 39-year-old woman. Preoperative study with constrained spherical deconvolution (CSD)–based tractography revealed the surrounding brain connectome architecture around the tumor relevant for safe surgical resection. Connectomic information provided detailed maps of the surrounding language and salience networks, including eloquent white matter fibers and cortical regions, which were visualized intraoperatively with image guidance and artificial intelligence (AI)–based brain mapping software. Microsurgical dissection is presented with detailed discussion of the safe boundaries and angles of resection when entering the insular operculum defined by connectomic information.The video can be found here: https://stream.cadmore.media/r10.3171/2021.10.FOCVID21194

Structural disconnection of the posterior medial frontal cortex reduces speech error monitoring

Optimal performance in any task relies on the ability to detect and correct errors. The anterior cingulate cortex and the broader posterior medial frontal cortex (pMFC) are active during error processing. However, it is unclear whether damage to the pMFC impairs error monitoring. We hypothesized that successful error monitoring critically relies on connections between the pMFC and broader cortical networks involved in executive functions and the task being monitored. We tested this hypothesis in the context of speech error monitoring in people with post-stroke aphasia. Diffusion weighted images were collected in 51 adults with chronic left-hemisphere stroke and 37 age-matched control participants. Whole-brain connectomes were derived using constrained spherical deconvolution and anatomically-constrained probabilistic tractography. Support vector regressions identified white matter connections in which lost integrity in stroke survivors related to reduced error detection during confrontation naming. Lesioned connections to the bilateral pMFC were related to reduce error monitoring, including many connections to regions associated with speech production and executive function. We conclude that connections to the pMFC support error monitoring. Error monitoring in speech production is supported by the structural connectivity between the pMFC and regions involved in speech production, comprehension, and executive function. Interactions between pMFC and other task-relevant processors may similarly be critical for error monitoring in other task contexts.

A robust framework for characterising diffusion metrics of the median and ulnar nerves: Exploiting state-of-the-art tracking methods.

Diffusion-weighted imaging has been used to quantify peripheral nerve properties; however, traditional post-processing techniques have several limitations. Advanced neuroimaging techniques, which overcome many of these limitations, have not been applied to peripheral nerves. Here, we use state-of-the-art diffusion analysis tools to reconstruct the median and ulnar nerves and quantify their diffusion properties. Diffusion-weighted MRI scans were obtained from eight healthy adult subjects. Constrained spherical deconvolution was combined with probabilistic fibre tracking to compute track-weighted fibre orientation distribution (TW-FOD). The tensor was computed and used along with the tracks to estimate TW apparent diffusion coefficient (TW-ADC), TW fractional anisotropy (TW-FA), TW axial diffusivity (TW-AD), and TW radial diffusivity (TW-RD). Variability of TW measurements was used to estimate power size information. The population intersession mean (± SD) measurements for the median nerve were TW-FOD 1.30 (±0.17), TW-ADC 1.16 (±0.13) × 10-3 mm2 /s, TW-FA 0.60 (±0.05), TW-AD 2.05 (±0.16) × 10-3 mm2 /s, and TW-RD 0.72 (±0.12) × 10-3 mm2 /s. The corresponding measurements for the ulnar nerve were TW-FOD 1.25 (±0.14), TW-ADC 1.13 (±0.10) × 10-3 mm2 /s, TW-FA 0.56 (±0.06), TW-AD 1.93 (±0.01) × 10-3 mm2 /s, and TW-RD 0.74 (±0.12) × 10-3 mm2 /s. Based on these measurements, a sample size of 37 would be sufficient to detect a 10% difference in any of the measured TW metrics. A sample size of 20 would be large enough to detect within-subject differences as small as 2.9% (TW-AD, ulnar nerve) and between-subject differences as small as 3.8% (TW-AD, ulnar nerve).

Alterations in white matter fiber in Parkinson disease across different cognitive stages

Cognitive impairment in Parkinson disease (PD) leads to substantial disability. Unlike external manifestations such as tremor, the decay of cognitive function is often an underlying process, and its neuroanatomic substrates are not yet fully elucidated. Knowledge regarding cognitive-related alterations in white matter (WM) pathways helps us understand the mechanisms of cognitive decline in patients with PD. Previous voxel-based analyses with Diffusion tensor imaging (DTI) metrics, such as fractional anisotropy (FA) and mean diffusivity (MD) have uncovered white matter differences in groupwise, but the conclusions were inconsistent. That was partially due to white matter fibers that are known to affect cognition, such as the corpus callosum (CC) and superior longitudinal fasciculus that cross in voxel, and are hard to interpreted by the abovementioned metrics. Furthermore, cognitive decay is a continuous process, it is difficult to reflect the continuous changes of white matter fibers between groups comparison. In the present work, we chose the constrained spherical deconvolution (CSD) and the fixel model, which avoided the effect of crossing fibers. To compare the white matter fiber in different cognitive stages of patients with PD, the results found that the CC, the cingulum bundle (CB), and the corticospinal tract (CST) showed the same trend in the decline of cognitive function, and this change may lead to the impairment of cognitive function. Our findings can help physicians determine the cognitive stage of PD from the perspective of white matter fiber and provide a reference for clinical trials and predictions.

Tracking the Corticospinal Tract in Patients With High-Grade Glioma: Clinical Evaluation of Multi-Level Fiber Tracking and Comparison to Conventional Deterministic Approaches.

While the diagnosis of high-grade glioma (HGG) is still associated with a considerably poor prognosis, neurosurgical tumor resection provides an opportunity for prolonged survival and improved quality of life for affected patients. However, successful tumor resection is dependent on a proper surgical planning to avoid surgery-induced functional deficits whilst achieving a maximum extent of resection (EOR). With diffusion magnetic resonance imaging (MRI) providing insight into individual white matter neuroanatomy, the challenge remains to disentangle that information as correctly and as completely as possible. In particular, due to the lack of sensitivity and accuracy, the clinical value of widely used diffusion tensor imaging (DTI)-based tractography is increasingly questioned. We evaluated whether the recently developed multi-level fiber tracking (MLFT) technique can improve tractography of the corticospinal tract (CST) in patients with motor-eloquent HGGs. Forty patients with therapy-naïve HGGs (mean age: 62.6 ± 13.4 years, 57.5% males) and preoperative diffusion MRI [repetition time (TR)/echo time (TE): 5000/78 ms, voxel size: 2x2x2 mm3, one volume at b=0 s/mm2, 32 volumes at b=1000 s/mm2] underwent reconstruction of the CST of the tumor-affected and unaffected hemispheres using MLFT in addition to deterministic DTI-based and deterministic constrained spherical deconvolution (CSD)-based fiber tractography. The brain stem was used as a seeding region, with a motor cortex mask serving as a target region for MLFT and a region of interest (ROI) for the other two algorithms. Application of the MLFT method substantially improved bundle reconstruction, leading to CST bundles with higher radial extent compared to the two other algorithms (delineation of CST fanning with a wider range; median radial extent for tumor-affected vs. unaffected hemisphere – DTI: 19.46° vs. 18.99°, p=0.8931; CSD: 30.54° vs. 27.63°, p=0.0546; MLFT: 81.17° vs. 74.59°, p=0.0134). In addition, reconstructions by MLFT and CSD-based tractography nearly completely included respective bundles derived from DTI-based tractography, which was however favorable for MLFT compared to CSD-based tractography (median coverage of the DTI-based CST for affected vs. unaffected hemispheres – CSD: 68.16% vs. 77.59%, p=0.0075; MLFT: 93.09% vs. 95.49%; p=0.0046). Thus, a more complete picture of the CST in patients with motor-eloquent HGGs might be achieved based on routinely acquired diffusion MRI data using MLFT.

Structural Brain Network Reproducibility: Influence of Different Diffusion Acquisition and Tractography Reconstruction Schemes on Graph Metrics.

Background: Graph metrics of structural brain networks demonstrate to be a powerful tool for investigating brain topology at a large scale. However, the variability of the results related to applying different magnetic resonance acquisition schemes and tractography reconstruction techniques is not fully characterized. Materials and Methods: The present work aims to evaluate the influence of different combinations of diffusion acquisition schemes (single and multishell), diffusion models (tensor and spherical deconvolution), and tractography reconstruction approaches (deterministic and probabilistic) on the reproducibility of graph metrics derived from structural connectome on test/retest (TRT) data released by the Human Connectome Project. From each implemented experimental setup, both global and local graph metrics were evaluated and their reproducibility was estimated by the intraclass correlation coefficient (ICC). Moreover, the percentage relative standard deviation (pRSD) from the ICC values of local graph metrics was calculated to quantify how much the reproducibility varied across nodes within each experimental setup. Results: The presented results show that different combinations of diffusion acquisition schemes, diffusion models, and tractography algorithms can strongly affect the reproducibility of global and local graph metrics. The combination of constrained spherical deconvolution (CSD) and deterministic tractography gave generally high reproducibility (ICCs >0.75) and lowest pRSD for the considered graph metrics, meanwhile probabilistic CSD with a high b-value returned the highest reproducibility. Notably, the introduction of streamline selection filters on CSD can substantially affect the reproducibility. Discussion: This work demonstrates that the TRT reproducibility of graph metrics is generally high but can vary substantially with different combinations of acquisition and reconstruction schemes. Impact statement This work demonstrates the influence of different diffusion acquisition schemes, diffusion models, and tractography reconstruction approaches on the reproducibility of graph metrics derived from structural connectome. The presented findings impact on the choice of both acquisition protocol and processing pipeline for topological analyses to produce reproducible measurements for brain network studies.

DETERMINATION OF FIBER ORIENTATION IN HIGH ANGULAR RESOLUTION DIFFUSION IMAGING USING SPHERICAL HARMONICS AND PARTICLE SWARM OPTIMIZATION

The fiber directions in High Angular Resolution Diffusion Imaging (HARDI) with low fractional anisotropy or low Signal to Noise Ratio (SNR) cannot be estimated accurately. In this paper, the fiber directions are estimated using Particle Swarm Optimization and Spherical Deconvolution (PSO-SD). Fiber orientation is modeled as a Dirac delta function in [Formula: see text]. The Spherical Harmonic Coefficients (SHC) of the Dirac delta function in the [Formula: see text] direction are obtained using the rotational harmonic matrix and the SHC of the Dirac delta function in the [Formula: see text]-axis. The PSO-SD method is used to determine ([Formula: see text]). We generated noise-free synthetic data for isotropic regions (FA varied from 0.1 to 0.8) and synthetic data with two crossing fibers for anisotropic regions with SNRs of 20, 15, 10 and 5 (FA [Formula: see text] 0.78). In the noise-free signal (FA [Formula: see text] 0.3), the Success Ratio (SR) and Mean Difference Angle (MDA) of the PSO-SD method were 1∘ and 9.48∘, respectively. In the noisy signal (FA [Formula: see text] 0.78, SNR [Formula: see text] 10, crossing angle [Formula: see text] 40), the SR and MDA of PSO-SD (with [Formula: see text]) were 0.46∘ and 12.3∘, respectively. The PSO-SD method can estimate fiber directions in HARDI with low fractional anisotropy and low SNR. Moreover, it has a higher SR and lower MDA in comparison with those of the super-CSD method.

Multi-tissue spherical deconvolution of tensor-valued diffusion MRI

Multi-tissue constrained spherical deconvolution (MT-CSD) leverages the characteristic b-value dependency of each tissue type to estimate both the apparent tissue densities and the white matter fiber orientation distribution function from diffusion MRI data. In this work, we generalize MT-CSD to tensor-valued diffusion encoding with arbitrary b-tensor shapes. This enables the use of data encoded with mixed b-tensors, rather than being limited to the subset of linear (conventional) b-tensors. Using the complete set of data, including all b-tensor shapes, provides a categorical improvement in the estimation of apparent tissue densities, fiber ODF, and resulting tractography. Furthermore, we demonstrate that including multiple b-tensor shapes in the analysis provides improved contrast between tissue types, in particular between gray matter and white matter. We also show that our approach provides high-quality apparent tissue density maps and high-quality fiber tracking from data, even with sparse sampling across b-tensors that yield whole-brain coverage at 2 mm isotropic resolution in approximately 5:15 min.

Mapping the Microstructure and Striae of the Human Olfactory Tract with Diffusion MRI.

The human sense of smell plays an important role in appetite and food intake, detecting environmental threats, social interactions, and memory processing. However, little is known about the neural circuity supporting its function. The olfactory tracts project from the olfactory bulb along the base of the frontal cortex, branching into several striae to meet diverse cortical regions. Historically, using diffusion magnetic resonance imaging (dMRI) to reconstruct the human olfactory tracts has been prevented by susceptibility and motion artifacts. Here, we used a dMRI method with readout segmentation of long variable echo-trains (RESOLVE) to minimize image distortions and characterize the human olfactory tracts in vivo. We collected high-resolution dMRI data from 25 healthy human participants (12 male and 13 female) and performed probabilistic tractography using constrained spherical deconvolution (CSD). At the individual subject level, we identified the lateral, medial, and intermediate striae with their respective cortical connections to the piriform cortex and amygdala (AMY), olfactory tubercle (OT), and anterior olfactory nucleus (AON). We combined individual results across subjects to create a normalized, probabilistic atlas of the olfactory tracts. We then investigated the relationship between olfactory perceptual scores and measures of white matter integrity, including mean diffusivity (MD). Importantly, we found that olfactory tract MD negatively correlated with odor discrimination performance. In summary, our results provide a detailed characterization of the connectivity of the human olfactory tracts and demonstrate an association between their structural integrity and olfactory perceptual function.SIGNIFICANCE STATEMENT This study provides the first detailed in vivo description of the cortical connectivity of the three olfactory tract striae in the human brain, using diffusion magnetic resonance imaging (dMRI). Additionally, we show that tract microstructure correlates with performance on an odor discrimination task, suggesting a link between the structural integrity of the olfactory tracts and odor perception. Lastly, we generated a normalized probabilistic atlas of the olfactory tracts that may be used in future research to study its integrity in health and disease.

Comparison of Diffusion Signal Models for Fiber Tractography in Eloquent Glioma Surgery–Determination of Accuracy Under Awake Craniotomy Conditions

Fiber tractography (FT) has become an important noninvasive tool to ensure maximal safe tumor resection in eloquent glioma surgery. Intraoperatively applied FT is still predominantly based on diffusion tensor imaging (DTI). However, reconstruction schemes of high angular resolution diffusion imaging data for high-resolution FT (HRFT) are gaining increasing attention. The aim of this prospective study was to compare the accuracy of sophisticated HRFT models compared with DTI-FT. Ten patients with eloquent gliomas underwent surgery under awake craniotomy conditions. The localization of acquisition points, representing deteriorations during intraoperative electrostimulation (IOM) and neuropsychological mapping, were documented. The offsets of acquisition points to the respective fiber bundle were calculated. Probabilistic Q-ball imaging (QBI) and constrained spherical deconvolution (CSD)-FT were compared with DTI-FT for the major language-associated fiber bundles (superior longitudinal fasciculus [SLF] II-IV, inferior fronto-occipital fasciculus, and inferior longitudinal fasciculus/medial longitudinal fasciculus). Among 186 offset values, 46% were located closer than 10 mm to the estimated fiber bundle (CSD, 36%; DTI, 40% and QBI, 60%). Moreover, only 10 offsets were further away than 30 mm (5%). Lowest mean minimum offsets (SLF, 7.7 ± 7.9 mm; inferior fronto-occipital fasciculus, 12.7 ± 8.3 mm; inferior longitudinal fasciculus/medial longitudinal fasciculus, 17.7 ± 6.7 mm) were found for QBI, indicating a significant advantage compared with CSD or DTI (P < 0.001), respectively. No significant differences were found between CSD-FT and DTI-FT offsets (P= 0.105), albeit for the compound SLF exclusively (P < 0.001). Comparing HRFT techniques QBI and CSD with DTI, QBI delivered significantly better results with lowest offsets and good correlation to IOM results. Besides, QBI-FT was feasible for neurosurgical preoperative and intraoperative applications. Our findings suggest that a combined approach of QBI-FT and IOM under awake craniotomy is considerable for best preservation of neurological function in the presented setting. Overall, the implementation of selected HRFT models into neuronavigation systems seems to be a promising tool in glioma surgery.

Whole-brain white matter network reorganization in HIV.

The human immunodeficiency virus (HIV) causes an infectious disease with a high viral tropism toward CD4 T-lymphocytes and macrophage. Since the advent of combined antiretroviral therapy (CART), the number of opportunistic infectious disease has diminished, turning HIV into a chronic condition. Nevertheless, HIV-infected patients suffer from several life-long symptoms, including the HIV-associated neurocognitive disorder (HAND), whose biological substrates remain unclear. HAND includes a range of cognitive impairments which have a huge impact on daily patient life. The aim of this study was to examine putative structural brain network changes in HIV-infected patient to test whether diffusion-imaging-related biomarkers could be used to discover and characterize subtle neurological alterations in HIV infection. To this end, we employed multi-shell, multi-tissue constrained spherical deconvolution in conjunction with probabilistic tractography and graph-theoretical analyses. We found several statistically significant effects in both local (right postcentral gyrus, right precuneus, right inferior parietal lobule, right transverse temporal gyrus, right inferior temporal gyrus, right putamen and right pallidum) and global graph-theoretical measures (global clustering coefficient, global efficiency and transitivity). Our study highlights a global and local reorganization of the structural connectome which support the possible application of graph theory to detect subtle alteration of brain regions in HIV patients.Clinical Relevance-Brain measures able to detect subtle alteration in HIV patients could also be used in e.g. evaluating therapeutic responses, hence empowering clinical trials.

The human dorsal hippocampal commissure: Delineating connections across the midline using multi-modal neuroimaging in major depressive disorder

Reclining on the lyra of the fornix, the dorsal hippocampal commissure (DHC) facilitates communication between the left and right hippocampi. Despite being a key commissural fiber bundle of the hippocampus, the DHC remains relatively underexplored in humans. Using high-resolution T1 and T2 magnetic resonance imaging in conjunction with High Angular Resolution Diffusion Imaging and Constrained Spherical Deconvolution tractography, we investigated if there are differences between the DHC of a sample of 54 healthy controls and 38 patients diagnosed with Major Depressive Disorder (MDD). After isolating the DHC using an anatomically derived protocol, segmented volumes of the left and right temporal pole, amygdala, hippocampus and parahippocampal gyrus extracted. In both controls and MDD patients, DHC fibers were spread across each of these four regions, with the parahippocampal gyri contributing the most to the volume of this white matter tract. Furthermore, diffusion metrics showed no effect of aging or sex on the integrity of the DHC in either group. This study is the first to investigate the DHC in MDD and the first to demonstrate that the DHC connects with the amygdala and temporal pole on both sides of the brain. Our data shows that the human DHC, rather than simply connecting the left and right hippocampi, allows connections between each hippocampus and the contralateral temporal pole, amygdala and parahippocampal gyrus. This suggests that the human temporal lobe, and the parahippocampal gyrus, may exert more influence on the contralateral limbic system via the DHC than previously thought.