Mapping white matter structural and network alterations in betel quid-dependent chewers using high angular resolution diffusion imaging

BACKGROUND

To evaluate brain white matter diffusion characteristics and anatomical network alterations in betel quid dependence (BQD) chewers using high angular resolution diffusion imaging (HARDI).

METHODS

The current study recruited 53 BQD chewers and 37 healthy controls (HC) in two groups. We explored regional diffusion metrics alternations in the BQD group compared with the HC group using automated fiber quantification (AFQ). We further employed the white matter (WM) anatomical network of HARDI to explore connectivity alterations in BQD chewers using graph theory.

RESULTS

BQD chewers presented significantly lower FA values in the left and right cingulum cingulate, the left and right thalamic radiation, and the right uncinate. The BQD has a significantly higher RD value in the right uncinate fasciculus than the HC group. At the global WM anatomical network level, global network efficiency (p = 0.008) was poorer and Lp (p = 0.016) was greater in the BQD group. At the nodal WM anatomical network level, nodal efficiency (p < 0.05) was lower in the BQD group.

CONCLUSION

Our findings provide novel morphometric evidence that brain structural changes in BQD are characterized by white matter diffusivity and anatomical network connectivity among regions of the brain, potentially leading to the enhanced reward system and impaired inhibitory control.

fMRI-Targeted High-Angular Resolution Diffusion MR Tractography to Identify Functional Language Tracts in Healthy Controls and Glioma Patients

BACKGROUND

MR Tractography enables non-invasive preoperative depiction of language subcortical tracts, which is crucial for the presurgical work-up of brain tumors; however, it cannot evaluate the exact function of the fibers.

PURPOSE

A systematic pipeline was developed to combine tractography reconstruction of language fiber bundles, based on anatomical landmarks (Anatomical-T), with language fMRI cortical activations. A fMRI-targeted Tractography (fMRI-T) was thus obtained, depicting the subsets of the anatomical tracts whose endpoints are located inside a fMRI activation. We hypothesized that fMRI-T could provide additional functional information regarding the subcortical structures, better reflecting the eloquent white matter structures identified intraoperatively.

METHODS

Both Anatomical-T and fMRI-T of language fiber tracts were performed on 16 controls and preoperatively on 16 patients with left-hemisphere brain tumors, using a q-ball residual bootstrap algorithm based on High Angular Resolution Diffusion Imaging (HARDI) datasets (b = 3000 s/mm2; 60 directions); fMRI ROIs were obtained using picture naming, verbal fluency, and auditory verb generation tasks. In healthy controls, normalized MNI atlases of fMRI-T and Anatomical-T were obtained. In patients, the surgical resection of the tumor was pursued by identifying eloquent structures with intraoperative direct electrical stimulation mapping and extending surgery to the functional boundaries. Post-surgical MRI allowed to identify Anatomical-T and fMRI-T non-eloquent portions removed during the procedure.

RESULTS

MNI Atlases showed that fMRI-T is a subset of Anatomical-T, and that different task-specific fMRI-T involve both shared subsets and task-specific subsets - e.g., verbal fluency fMRI-T strongly involves dorsal frontal tracts, consistently with the phonogical-articulatory features of this task. A quantitative analysis in patients revealed that Anatomical-T removed portions of AF-SLF and IFOF were significantly greater than verbal fluency fMRI-T ones, suggesting that fMRI-T is a more specific approach. In addition, qualitative analyses showed that fMRI-T AF-SLF and IFOF predict the exact functional limits of resection with increased specificity when compared to Anatomical-T counterparts, especially the superior frontal portion of IFOF, in a subcohort of patients.

CONCLUSION

These results suggest that performing fMRI-T in addition to the 'classic' Anatomical-T may be useful in a preoperative setting to identify the 'high-risk subsets' that should be spared during the surgical procedure.

Projections of Brodmann Area 6 to the Pyramidal Tract in Humans: Quantifications Using High Angular Resolution Data

Primate studies indicate that the pyramidal tract (PyT) could originate from Brodmann area (BA) 6. However, in humans, the accurate origin of PyT from BA 6 is still uncertain owing to difficulties in visualizing anatomical features such as the fanning shape at the corona radiata and multiple crossings at the semioval centrum. High angular-resolution diffusion imaging (HARDI) could reliably replicate these anatomical features. We explored the origin of the human PyT from BA 6 using HARDI. With HARDI data of 30 adults from the Massachusetts General Hospital-Human Connectome Project (MGH-HCP) database and the HCP 1021 template (average of 1021 HCP diffusion data), we visualized the PyT at the 30-averaged group level and the 1021 large-sample level and validated the observations in each of the individuals. Endpoints of the fibers within each subregion were quantified. PyT fibers originating from the BA 6 were consistently visualized in all images. Specifically, the bilateral supplementary motor area (SMA) and dorsal premotor area (dPMA) were consistently found to contribute to the PyT. PyT fibers from BA 6 and those from BA 4 exhibited a twisting topology. The PyT contains fibers originating from the SMA and dPMA in BA 6. Infarction of these regions or aging would result in incomplete provision of information to the PyT and concomitant decreases in motor planning and coordination abilities.

Effects of Brain Parcellation on the Characterization of Topological Deterioration in Alzheimer's Disease

Alzheimer's disease (AD) causes the progressive deterioration of neural connections, disrupting structural connectivity (SC) networks within the brain. Graph-based analyses of SC networks have shown that topological properties can reveal the course of AD propagation. Different whole-brain parcellation schemes have been developed to define the nodes of these SC networks, although it remains unclear which scheme can best describe the AD-related deterioration of SC networks. In this study, four whole-brain parcellation schemes with different numbers of parcels were used to define SC network nodes. SC networks were constructed based on high angular resolution diffusion imaging (HARDI) tractography for a mixed cohort that includes 20 normal controls (NC), 20 early mild cognitive impairment (EMCI), 20 late mild cognitive impairment (LMCI), and 20 AD patients, from the Alzheimer's Disease Neuroimaging Initiative. Parcellation schemes investigated in this study include the OASIS-TRT-20 (62 regions), AAL (116 regions), HCP-MMP (180 regions), and Gordon-rsfMRI (333 regions), which have all been widely used for the construction of brain structural or functional connectivity networks. Topological characteristics of the SC networks, including the network strength, global efficiency, clustering coefficient, rich-club, characteristic path length, k-core, rich-club coefficient, and modularity, were fully investigated at the network level. Statistical analyses were performed on these metrics using Kruskal-Wallis tests to examine the group differences that were apparent at different stages of AD progression. Results suggest that the HCP-MMP scheme is the most robust and sensitive to AD progression, while the OASIS-TRT-20 scheme is sensitive to group differences in network strength, global efficiency, k-core, and rich-club coefficient at k-levels from 18 and 39. With the exception of the rich-club and modularity coefficients, AAL could not significantly identify group differences on other topological metrics. Further, the Gordon-rsfMRI atlas only significantly differentiates the groups on network strength, characteristic path length, k-core, and rich-club coefficient. Results show that the topological examination of SC networks with different parcellation schemes can provide important complementary AD-related information and thus contribute to a more accurate and earlier diagnosis of AD.

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

BACKGROUND

In Parkinson's disease cognitive impairment is an early nonmotor feature, but it is still unclear why some patients are able to maintain their cognitive performance at normal levels, as quantified by neuropsychological tests, whereas others cannot. The objectives of this study were to perform a cross-sectional study and analyze the white matter changes in the cognitive and motor bundles in patients with Parkinson's disease.

METHODS

Sixteen Parkinson's disease patients with normal cognitive performance, 19 with mild cognitive impairment (based on their performance of 1.5 standard deviations below the healthy population mean), and 16 healthy controls were compared with respect to their tractography patterns between the cortical cognitive / motor regions and subcortical structures, using high angular resolution diffusion imaging and constrained spherical deconvolution computation.

RESULTS

Motor bundles showed decreased apparent fiber density in both PD groups, associated with a significant increase in diffusivity metrics, number of reconstructed streamlines, and track volumes, compared with healthy controls. By contrast, in the cognitive bundles, decreased fiber density in both Parkinson's groups was compounded by the absence of changes in diffusivity in patients with normal cognition, whereas patients with cognitive impairment had increased diffusivity metrics, lower numbers of reconstructed streamlines, and lower track volumes.

CONCLUSIONS

Both PD groups showed similar patterns of white matter neurodegeneration in the motor bundles, whereas cognitive bundles showed a distinct pattern: Parkinson's patients with normal cognition had white matter diffusivity metrics similar to healthy controls, whereas in patients with cognitive impairment white matter showed a neurodegeneration pattern. © 2018 International Parkinson and Movement Disorder Society.

Empirical estimation of intravoxel structure with persistent angular structure and Q-ball models of diffusion weighted MRI

The diffusion tensor model is nonspecific in regions where micrometer structural patterns are inconsistent at the millimeter scale (i.e., brain regions with pathways that cross, bend, branch, fan, etc.). Numerous models have been proposed to represent crossing fibers and complex intravoxel structure from in vivo diffusion weighted magnetic resonance imaging (e.g., high angular resolution diffusion imaging-HARDI). Here, we present an empirical comparison of two HARDI approaches-persistent angular structure MRI (PAS-MRI) and Q-ball-using a newly acquired reproducibility dataset. Briefly, a single subject was scanned 11 times with 96 diffusion weighted directions and 10 reference volumes for each of two [Formula: see text] values (1000 and [Formula: see text] for a total of 2144 volumes). Empirical reproducibility of intravoxel fiber fractions (number/strength of peaks), angular orientation, and fractional anisotropy was compared with metrics from a traditional tensor analysis approach, focusing on [Formula: see text] values of 1000 and [Formula: see text]. PAS-MRI is shown to be more reproducible than Q-ball and offers advantages at low [Formula: see text] values. However, there are substantial and biologically meaningful differences between the intravoxel structures estimated both in terms of analysis method as well as by [Formula: see text] value. The two methods suggest a fundamentally different microarchitecture of the human brain; therefore, it is premature to perform meta-analysis or combine results across HARDI studies using a different analysis model or acquisition sequences.

A systematic evaluation of intraoperative white matter tract shift in pediatric epilepsy surgery using high-field MRI and probabilistic high angular resolution diffusion imaging tractography

OBJECTIVE Characterization of intraoperative white matter tract (WMT) shift has the potential to compensate for neuronavigation inaccuracies using preoperative brain imaging. This study aimed to quantify and characterize intraoperative WMT shift from the global hemispheric to the regional tract-based scale and to investigate the impact of intraoperative factors (IOFs). METHODS High angular resolution diffusion imaging (HARDI) diffusion-weighted data were acquired over 5 consecutive perioperative time points (MR₁ to MR₅) in 16 epilepsy patients (8 male; mean age 9.8 years, range 3.8-15.8 years) using diagnostic and intraoperative 3-T MRI scanners. MR₁ was the preoperative planning scan. MR₂ was the first intraoperative scan acquired with the patient's head fixed in the surgical position. MR₃ was the second intraoperative scan acquired following craniotomy and durotomy, prior to lesion resection. MR₄ was the last intraoperative scan acquired following lesion resection, prior to wound closure. MR₅ was a postoperative scan acquired at the 3-month follow-up visit. Ten association WMT/WMT segments and 1 projection WMT were generated via a probabilistic tractography algorithm from each MRI scan. Image registration was performed through pairwise MRI alignments using the skull segmentation. The MR₁ and MR₂ pairing represented the first surgical stage. The MR₂ and MR₃ pairing represented the second surgical stage. The MR₃ and MR₄ (or MR₅) pairing represented the third surgical stage. The WMT shift was quantified by measuring displacements between a pair of WMT centerlines. Linear mixed-effects regression analyses were carried out for 6 IOFs: head rotation, craniotomy size, durotomy size, resected lesion volume, presence of brain edema, and CSF loss via ventricular penetration. RESULTS The average WMT shift in the operative hemisphere was 2.37 mm (range 1.92-3.03 mm) during the first surgical stage, 2.19 mm (range 1.90-3.65 mm) during the second surgical stage, and 2.92 mm (range 2.19-4.32 mm) during the third surgical stage. Greater WMT shift occurred in the operative than the nonoperative hemisphere, in the WMTs adjacent to the surgical lesion rather than those remote to it, and in the superficial rather than the deep segment of the pyramidal tract. Durotomy size and resection size were significant, independent IOFs affecting WMT shift. The presence of brain edema was a marginally significant IOF. Craniotomy size, degree of head rotation, and ventricular penetration were not significant IOFs affecting WMT shift. CONCLUSIONS WMT shift occurs noticeably in tracts adjacent to the surgical lesions, and those motor tracts superficially placed in the operative hemisphere. Intraoperative probabilistic HARDI tractography following craniotomy, durotomy, and lesion resection may compensate for intraoperative WMT shift and improve neuronavigation accuracy.

Fast diffusion imaging with high angular resolution

PURPOSE

High angular resolution diffusion imaging (HARDI) is a well-established method to help reveal the architecture of nerve bundles, but long scan times and geometric distortions inherent to echo planar imaging (EPI) have limited its integration into clinical protocols.

METHODS

A fast imaging method is proposed here that combines accelerated multishot diffusion imaging (AMDI), multiplexed sensitivity encoding (MUSE), and crossing fiber angular resolution of intravoxel structure (CFARI) to reduce spatial distortions and reduce total scan time. A multishot EPI sequence was used to improve geometrical fidelity as compared to a single-shot EPI acquisition, and acceleration in both k-space and diffusion sampling enabled reductions in scan time. The method is regularized and self-navigated for motion correction. Seven volunteers were scanned in this study, including four with volumetric whole brain acquisitions.

RESULTS

The average similarity of microstructural orientations between undersampled datasets and their fully sampled counterparts was above 85%, with scan times below 5 min for whole-brain acquisitions. Up to 2.7-fold scan time acceleration along with four-fold distortion reduction was achieved.

CONCLUSION

The proposed imaging strategy can generate HARDI results with relatively good geometrical fidelity and low scan duration, which may help facilitate the transition of HARDI from a successful research tool to a practical clinical one. Magn Reson Med 77:696-706, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Anatomic connectivity assessed using pathway radial diffusivity is related to functional connectivity in monosynaptic pathways

This work presents a pathway-dependent anatomic and functional connectivity analysis in 19 patients with relapse-remitting multiple sclerosis (MS) and 16 age-, education-, and gender-matched controls. An MS population is used in this study as a model for anatomic connectivity, permitting us to observe relationships between anatomic and functional connectivity more easily. A combined resting-state functional magnetic resonance imaging (fMRI) and whole-brain, high angular resolution diffusion imaging analysis is performed in three independent, monosynaptic pathways. The pathways chosen were transcallosal pathway connecting the bilateral primary sensorimotor regions, right and left posterior portion of the Papez circuit, connecting the posterior cingulate cortex and hippocampus. The Papez circuit is known to be involved in memory function, one of the most frequently impacted cognitive domains in patients with MS. We show that anatomic connectivity, as measured with diffusion-weighted imaging, and functional connectivity, as measured with resting-state fMRI, are significantly reduced in patients as compared with controls for at least some of the pathways considered. In addition when all pathway measures are combined, anatomic and functional connectivity are significantly correlated in patients with MS as well as healthy controls. We suggest that anatomic and functional connectivity are related for monosynaptic pathways and that radial diffusivity, as a diffusion-tensor-based measure of white matter integrity, is a robust measure of anatomic connectivity in the general population.

Diffusivity signatures characterize trigeminal neuralgia associated with multiple sclerosis

BACKGROUND

Trigeminal neuralgia secondary to multiple sclerosis (MS-TN) is a facial neuropathic pain syndrome similar to classic trigeminal neuralgia (TN). While TN is caused by neurovascular compression of the fifth cranial nerve (CN V), how MS-related demyelination correlates with pain in MS-TN is not understood.

OBJECTIVES

We aim to examine diffusivities along CN V in MS-TN, TN, and controls in order to reveal differential neuroimaging correlates across groups.

METHODS

3T MR diffusion weighted, T1, T2 and FLAIR sequences were acquired for MS-TN, TN, and controls. Multi-tensor tractography was used to delineate CN V across cisternal, root entry zone (REZ), pontine and peri-lesional segments. Diffusion metrics including fractional anisotropy (FA), and radial (RD), axial (AD), and mean diffusivities (MD) were measured from each segment.

RESULTS

CN V segments showed distinctive diffusivity patterns. The TN group showed higher FA in the cisternal segment ipsilateral to the side of pain, and lower FA in the ipsilateral REZ segment. The MS-TN group showed lower FA in the ipsilateral peri-lesional segments, suggesting differential microstructural changes along CN V in these conditions.

CONCLUSIONS

The study demonstrates objective differences in CN V microstrucuture in TN and MS-TN using non-invasive neuroimaging. This represents a significant improvement in the methods currently available to study pain in MS.

Beyond crossing fibers: bootstrap probabilistic tractography using complex subvoxel fiber geometries

Diffusion magnetic resonance imaging fiber tractography is a powerful tool for investigating human white matter connectivity in vivo. However, it is prone to false positive and false negative results, making interpretation of the tractography result difficult. Optimal tractography must begin with an accurate description of the subvoxel white matter fiber structure, includes quantification of the uncertainty in the fiber directions obtained, and quantifies the confidence in each reconstructed fiber tract. This paper presents a novel and comprehensive pipeline for fiber tractography that meets the above requirements. The subvoxel fiber geometry is described in detail using a technique that allows not only for straight crossing fibers but for fibers that curve and splay. This technique is repeatedly performed within a residual bootstrap statistical process in order to efficiently quantify the uncertainty in the subvoxel geometries obtained. A robust connectivity index is defined to quantify the confidence in the reconstructed connections. The tractography pipeline is demonstrated in the human brain.

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

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

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

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

White matter microstructure in body dysmorphic disorder and its clinical correlates

Body dysmorphic disorder (BDD) is characterized by an often-delusional preoccupation with misperceived defects of appearance, causing significant distress and disability. Although previous studies have found functional abnormalities in visual processing, frontostriatal, and limbic systems, no study to date has investigated the microstructure of white matter connecting these systems in BDD. Participants comprised 14 medication-free individuals with BDD and 16 healthy controls who were scanned using diffusion-weighted magnetic resonance imaging (MRI). We utilized probabilistic tractography to reconstruct tracts of interest, and tract-based spatial statistics to investigate whole brain white matter. To estimate white matter microstructure, we used fractional anisotropy (FA), mean diffusivity (MD), and linear and planar anisotropy (c(l) and c(p)). We correlated diffusion measures with clinical measures of symptom severity and poor insight/delusionality. Poor insight negatively correlated with FA and c(l) and positively correlated with MD in the inferior longitudinal fasciculus (ILF) and the forceps major (FM). FA and c(l) were lower in the ILF and the inferior fronto-occipital fasciculus and higher in the FM in the BDD group, but differences were nonsignificant. This is the first diffusion-weighted MR investigation of white matter in BDD. Results suggest a relationship between impairments in insight, a clinically important phenotype, and fiber disorganization in tracts connecting visual with emotion/memory processing systems.

Effects of chronic mild traumatic brain injury on white matter integrity in Iraq and Afghanistan war veterans

Mild traumatic brain injury (TBI) is a common source of morbidity from the wars in Iraq and Afghanistan. With no overt lesions on structural MRI, diagnosis of chronic mild TBI in military veterans relies on obtaining an accurate history and assessment of behavioral symptoms that are also associated with frequent comorbid disorders, particularly posttraumatic stress disorder (PTSD) and depression. Military veterans from Iraq and Afghanistan with mild TBI (n = 30) with comorbid PTSD and depression and non-TBI participants from primary (n = 42) and confirmatory (n = 28) control groups were assessed with high angular resolution diffusion imaging (HARDI). White matter-specific registration followed by whole-brain voxelwise analysis of crossing fibers provided separate partial volume fractions reflecting the integrity of primary fibers and secondary (crossing) fibers. Loss of white matter integrity in primary fibers (P < 0.05; corrected) was associated with chronic mild TBI in a widely distributed pattern of major fiber bundles and smaller peripheral tracts including the corpus callosum (genu, body, and splenium), forceps minor, forceps major, superior and posterior corona radiata, internal capsule, superior longitudinal fasciculus, and others. Distributed loss of white matter integrity correlated with duration of loss of consciousness and most notably with "feeling dazed or confused," but not diagnosis of PTSD or depressive symptoms. This widespread spatial extent of white matter damage has typically been reported in moderate to severe TBI. The diffuse loss of white matter integrity appears consistent with systemic mechanisms of damage shared by blast- and impact-related mild TBI that involves a cascade of inflammatory and neurochemical events.

Angular versus spatial resolution trade-offs for diffusion imaging under time constraints

Diffusion weighted magnetic resonance imaging (DW-MRI) are now widely used to assess brain integrity in clinical populations. The growing interest in mapping brain connectivity has made it vital to consider what scanning parameters affect the accuracy, stability, and signal-to-noise of diffusion measures. Trade-offs between scan parameters can only be optimized if their effects on various commonly-derived measures are better understood. To explore angular versus spatial resolution trade-offs in standard tensor-derived measures, and in measures that use the full angular information in diffusion signal, we scanned eight subjects twice, 2 weeks apart, using three protocols that took the same amount of time (7 min). Scans with 3.0, 2.7, 2.5 mm isotropic voxels were collected using 48, 41, and 37 diffusion-sensitized gradients to equalize scan times. A specially designed DTI phantom was also scanned with the same protocols, and different b-values. We assessed how several diffusion measures including fractional anisotropy (FA), mean diffusivity (MD), and the full 3D orientation distribution function (ODF) depended on the spatial/angular resolution and the SNR. We also created maps of stability over time in the FA, MD, ODF, skeleton FA of 14 TBSS-derived ROIs, and an information uncertainty index derived from the tensor distribution function, which models the signal using a continuous mixture of tensors. In scans of the same duration, higher angular resolution and larger voxels boosted SNR and improved stability over time. The increased partial voluming in large voxels also led to bias in estimating FA, but this was partially addressed by using "beyond-tensor" models of diffusion.

A nonparametric Riemannian framework for processing high angular resolution diffusion images and its applications to ODF-based morphometry

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.

HARDI DATA DENOISING USING VECTORIAL TOTAL VARIATION AND LOGARITHMIC BARRIER

In this work, we wish to denoise HARDI (High Angular Resolution Diffusion Imaging) data arising in medical brain imaging. Diffusion imaging is a relatively new and powerful method to measure the three-dimensional profile of water diffusion at each point in the brain. These images can be used to reconstruct fiber directions and pathways in the living brain, providing detailed maps of fiber integrity and connectivity. HARDI data is a powerful new extension of diffusion imaging, which goes beyond the diffusion tensor imaging (DTI) model: mathematically, intensity data is given at every voxel and at any direction on the sphere. Unfortunately, HARDI data is usually highly contaminated with noise, depending on the b-value which is a tuning parameter pre-selected to collect the data. Larger b-values help to collect more accurate information in terms of measuring diffusivity, but more noise is generated by many factors as well. So large b-values are preferred, if we can satisfactorily reduce the noise without losing the data structure. Here we propose two variational methods to denoise HARDI data. The first one directly denoises the collected data S, while the second one denoises the so-called sADC (spherical Apparent Diffusion Coefficient), a field of radial functions derived from the data. These two quantities are related by an equation of the form S = S(S)exp (-b · sADC) (in the noise-free case). By applying these two different models, we will be able to determine which quantity will most accurately preserve data structure after denoising. The theoretical analysis of the proposed models is presented, together with experimental results and comparisons for denoising synthetic and real HARDI data.

The tensor distribution function

Diffusion weighted magnetic resonance imaging is a powerful tool that can be employed to study white matter microstructure by examining the 3D displacement profile of water molecules in brain tissue. By applying diffusion-sensitized gradients along a minimum of six directions, second-order tensors (represented by three-by-three positive definite matrices) can be computed to model dominant diffusion processes. However, conventional DTI is not sufficient to resolve more complicated white matter configurations, e.g., crossing fiber tracts. Recently, a number of high-angular resolution schemes with more than six gradient directions have been employed to address this issue. In this article, we introduce the tensor distribution function (TDF), a probability function defined on the space of symmetric positive definite matrices. Using the calculus of variations, we solve the TDF that optimally describes the observed data. Here, fiber crossing is modeled as an ensemble of Gaussian diffusion processes with weights specified by the TDF. Once this optimal TDF is determined, the orientation distribution function (ODF) can easily be computed by analytic integration of the resulting displacement probability function. Moreover, a tensor orientation distribution function (TOD) may also be derived from the TDF, allowing for the estimation of principal fiber directions and their corresponding eigenvalues.

Q-ball imaging of macaque white matter architecture

Diffusion-weighted magnetic resonance imaging holds substantial promise as a technique for non-invasive imaging of white matter (WM) axonal projections. For diffusion imaging to be capable of providing new insight into the connectional neuroanatomy of the human brain, it will be necessary to histologically validate the technique against established tracer methods such as horseradish peroxidase and biocytin histochemistry. The macaque monkey provides an ideal model for histological validation of the diffusion imaging method due to the phylogenetic proximity between humans and macaques, the gyrencephalic structure of the macaque cortex, the large body of knowledge on the neuroanatomic connectivity of the macaque brain and the ability to use comparable magnetic resonance acquisition protocols in both species. Recently, it has been shown that high angular resolution diffusion imaging (HARDI) can resolve multiple axon orientations within an individual imaging voxel in human WM. This capability promises to boost the accuracy of tract reconstructions from diffusion imaging. If the macaque is to serve as a model for histological validation of the diffusion tractography method, it will be necessary to show that HARDI can also resolve intravoxel architecture in macaque WM. The present study therefore sought to test whether the technique can resolve intravoxel structure in macaque WM. Using a HARDI method called q-ball imaging (QBI) it was possible to resolve composite intravoxel architecture in a number of anatomic regions. QBI resolved intravoxel structure in, for example, the dorsolateral convexity, the pontine decussation, the pulvinar and temporal subcortical WM. The paper concludes by reviewing remaining challenges for the diffusion tractography project.