Spherical Deconvolution Research Articles

Personalized biomechanical tongue models based on diffusion-weighted MRI and validated using optical tracking of range of motion

For advanced tongue cancer, the choice between surgery and organ-sparing treatment is often dependent on the expected loss of tongue functionality after treatment. Biomechanical models might assist in this choice by simulating the post-treatment function loss. However, this function loss varies between patients and should, therefore, be predicted for each patient individually. In the present study, the goal was to better predict the postoperative range of motion (ROM) of the tongue by personalizing biomechanical models using diffusion-weighted MRI and constrained spherical deconvolution reconstructions of tongue muscle architecture. Diffusion-weighted MRI scans of ten healthy volunteers were obtained to reconstruct their tongue musculature, which were subsequently registered to a previously described population average or atlas. Using the displacement fields obtained from the registration, the segmented muscle fiber tracks from the atlas were morphed back to create personalized muscle fiber tracks. Finite element models were created from the fiber tracks of the atlas and those of the individual tongues. Via inverse simulation of a protruding, downward, left and right movement, the ROM of the tongue was predicted. This prediction was compared to the ROM measured with a 3D camera. It was demonstrated that biomechanical models with personalized muscles bundles are better in approaching the measured ROM than a generic model. However, to achieve this result a correction factor was needed to compensate for the small magnitude of motion of the model. Future versions of these models may have the potential to improve the estimation of function loss after treatment for advanced tongue cancer.

From Neurosurgical Planning to Histopathological Brain Tumor Characterization: Potentialities of Arcuate Fasciculus Along-Tract Diffusion Tensor Imaging Tractography Measures.

Background: Tractography has been widely adopted to improve brain gliomas' surgical planning and guide their resection. This study aimed to evaluate state-of-the-art of arcuate fasciculus (AF) tractography for surgical planning and explore the role of along-tract analyses in vivo for characterizing tumor histopathology.Methods: High angular resolution diffusion imaging (HARDI) images were acquired for nine patients with tumors located in or near language areas (age: 41 ± 14 years, mean ± standard deviation; five males) and 32 healthy volunteers (age: 39 ± 16 years; 16 males). Phonemic fluency task fMRI was acquired preoperatively for patients. AF tractography was performed using constrained spherical deconvolution diffusivity modeling and probabilistic fiber tracking. Along-tract analyses were performed, dividing the AF into 15 segments along the length of the tract defined using the Laplacian operator. For each AF segment, diffusion tensor imaging (DTI) measures were compared with those obtained in healthy controls (HCs). The hemispheric laterality index (LI) was calculated from language task fMRI activations in the frontal, parietal, and temporal lobe parcellations. Tumors were grouped into low/high grade (LG/HG).Results: Four tumors were LG gliomas (one dysembryoplastic neuroepithelial tumor and three glioma grade II) and five HG gliomas (two grade III and three grade IV). For LG tumors, gross total removal was achieved in all but one case, for HG in two patients. Tractography identified the AF trajectory in all cases. Four along-tract DTI measures potentially discriminated LG and HG tumor patients (false discovery rate < 0.1): the number of abnormal MD and RD segments, median AD, and MD measures. Both a higher number of abnormal AF segments and a higher AD and MD measures were associated with HG tumor patients. Moreover, correlations (unadjusted p < 0.05) were found between the parietal lobe LI and the DTI measures, which discriminated between LG and HG tumor patients. In particular, a more rightward parietal lobe activation (LI < 0) correlated with a higher number of abnormal MD segments (R = −0.732) and RD segments (R = −0.724).Conclusions: AF tractography allows to detect the course of the tract, favoring the safer-as-possible tumor resection. Our preliminary study shows that along-tract DTI metrics can provide useful information for differentiating LG and HG tumors during pre-surgical tumor characterization.

Detecting Corticospinal Tract Impairment in Tumor Patients With Fiber Density and Tensor-Based Metrics.

Tumors infiltrating the motor system lead to significant disability, often caused by corticospinal tract injury. The delineation of the healthy-pathological white matter (WM) interface area, for which diffusion magnetic resonance imaging (dMRI) has shown promising potential, may improve treatment outcome. However, up to 90% of white matter (WM) voxels include multiple fiber populations, which cannot be correctly described with traditional metrics such as fractional anisotropy (FA) or apparent diffusion coefficient (ADC). Here, we used a novel fixel-based along-tract analysis consisting of constrained spherical deconvolution (CSD)-based probabilistic tractography and fixel-based apparent fiber density (FD), capable of identifying fiber orientation specific microstructural metrics. We addressed this novel methodology’s capability to detect corticospinal tract impairment. We measured and compared tractogram-related FD and traditional microstructural metrics bihemispherically in 65 patients with WHO grade III and IV gliomas infiltrating the motor system. The cortical tractogram seeds were based on motor maps derived by transcranial magnetic stimulation. We extracted 100 equally distributed cross-sections along each streamline of corticospinal tract (CST) for along-tract statistical analysis. Cross-sections were then analyzed to detect differences between healthy and pathological hemispheres. All metrics showed significant differences between healthy and pathologic hemispheres over the entire tract and between peritumoral segments. Peritumoral values were lower for FA and FD, but higher for ADC within the entire cohort. FD was more specific to tumor-induced changes in CST than ADC or FA, whereas ADC and FA showed higher sensitivity. The bihemispheric along-tract analysis provides an approach to detect subject-specific structural changes in healthy and pathological WM. In the current clinical dataset, the more complex FD metrics did not outperform FA and ADC in terms of describing corticospinal tract impairment.

Equivariant Spherical Deconvolution: Learning Sparse Orientation Distribution Functions from Spherical Data.

We present a rotation-equivariant self-supervised learning framework for the sparse deconvolution of non-negative scalar fields on the unit sphere. Spherical signals with multiple peaks naturally arise in Diffusion MRI (dMRI), where each voxel consists of one or more signal sources corresponding to anisotropic tissue structure such as white matter. Due to spatial and spectral partial voluming, clinically-feasible dMRI struggles to resolve crossing-fiber white matter configurations, leading to extensive development in spherical deconvolution methodology to recover underlying fiber directions. However, these methods are typically linear and struggle with small crossing-angles and partial volume fraction estimation. In this work, we improve on current methodologies by nonlinearly estimating fiber structures via self-supervised spherical convolutional networks with guaranteed equivariance to spherical rotation. We perform validation via extensive single and multi-shell synthetic benchmarks demonstrating competitive performance against common base-lines. We further show improved downstream performance on fiber tractography measures on the Tractometer benchmark dataset. Finally, we show downstream improvements in terms of tractography and partial volume estimation on a multi-shell dataset of human subjects.

Tracking white-matter brain modifications in chronic non-bothersome acoustic trauma tinnitus

Subjective tinnitus is a symptom characterized by the perception of sound with no external acoustic source, most often accompanied by co-morbidities. To date, the specific role of white matter abnormalities related to tinnitus reaches no consensus in the literature. The goal of this study was to explore the structural connectivity related to tinnitus percept per se, thus focusing on a specific population presenting chronic non-bothersome tinnitus of similar etiology (noise induced) without co-morbidities. We acquired diffusion-weighted images with high angular resolution in a homogeneous group of mildly impacted tinnitus participants (n = 19) and their matched controls (n = 19). We focused the study on two subsets of fiber bundles of interest: on one hand, we extracted the acoustic radiation and further included any intersecting fiber bundles; on the other hand, we explored the tracts related to the limbic system. We modeled the diffusion signal using constrained spherical deconvolution. We conducted a deep-learning based tractography segmentation and mapped Apparent Fiber Density (AFD) on the bundles of interest. C, as well as Fractional Anisotropy (FA) and FOD peak amplitude for comparison. Between group statistical comparison was performed along the 27 tracts of interest controlling for confounding hearing loss, tinnitus severity, and duration since onset. We tested a potential correlation with hearing loss, tinnitus duration and tinnitus handicap score along these tracts. In the tinnitus group, we observed increased AFD related to chronic tinnitus percept after acoustic trauma in two main white matter regions. First, in the right hemisphere, in the isthmus between inferior temporal and inferior frontal cortices, in the uncinate fasciculus (UF), and in the inferior fronto-occipital bundle (IFO). Second, in the left hemisphere, underneath the superior parietal region in the thalamo parietal tract and parieto-occipital pontine tract. Between-group differences in the acoustic radiations were not significant with AFD but were with FA. Furthermore, significant correlations with hearing loss were found in the left hemisphere in the inferior longitudinal fasciculus and in the fronto-pontine tract. No additional correlation was found with tinnitus duration nor with tinnitus handicap, as reflected by THI scores. The regions that displayed tinnitus related increased AFD also displayed increased FA. The isthmus of the UF and IFO in the right hemisphere appear to be involved with a number of neuropsychiatric and traumatic disorders confirming the involvement of the limbic system even in chronic non-bothersome tinnitus subjects, potentially suggesting a common pathway between these pathologies. White matter changes underneath the superior parietal cortex found here in tinnitus participants supports the implication of an auditory-somatosensory pathway in tinnitus perception.

Tractographic and Microstructural Analysis of the Dentato-Rubro-Thalamo-Cortical Tracts in Children Using Diffusion MRI.

The dentato-rubro-thalamo-cortical tract (DRTC) is the main outflow pathway of the cerebellum, contributing to a finely balanced corticocerebellar loop involved in cognitive and sensorimotor functions. Damage to the DRTC has been implicated in cerebellar mutism syndrome seen in up to 25% of children after cerebellar tumor resection. Multi-shell diffusion MRI (dMRI) combined with quantitative constrained spherical deconvolution tractography and multi-compartment spherical mean technique modeling was used to explore the frontocerebellar connections and microstructural signature of the DRTC in 30 healthy children. The highest density of DRTC connections were to the precentral (M1) and superior frontal gyri (F1), and from cerebellar lobules I-IV and IX. The first evidence of a topographic organization of anterograde projections to the frontal cortex at the level of the superior cerebellar peduncle (SCP) is demonstrated, with streamlines terminating in F1 lying dorsomedially in the SCP compared to those terminating in M1. The orientation dispersion entropy of DRTC regions appears to exhibit greater contrast than that shown by fractional anisotropy. Analysis of a separate reproducibility cohort demonstrates good consistency in the dMRI metrics described. These novel anatomical insights into this well-studied pathway may prove to be of clinical relevance in the surgical resection of cerebellar tumors.

Diffusion MRI tractography filtering techniques change the topology of structural connectomes

Objective. The use of non-invasive techniques for the estimation of structural brain networks (i.e. connectomes) opened the door to large-scale investigations on the functioning and the architecture of the brain, unveiling the link between neurological disorders and topological changes of the brain network. This study aims at assessing if and how the topology of structural connectomes estimated non-invasively with diffusion MRI is affected by the employment of tractography filtering techniques in structural connectomic pipelines. Additionally, this work investigates the robustness of topological descriptors of filtered connectomes to the common practice of density-based thresholding. Approach. We investigate the changes in global efficiency, characteristic path length, modularity and clustering coefficient on filtered connectomes obtained with the spherical deconvolution informed filtering of tractograms and using the convex optimization modelling for microstructure informed tractography. The analysis is performed on both healthy subjects and patients affected by traumatic brain injury and with an assessment of the robustness of the computed graph-theoretical measures with respect to density-based thresholding of the connectome. Main results. Our results demonstrate that tractography filtering techniques change the topology of brain networks, and thus alter network metrics both in the pathological and the healthy cases. Moreover, the measures are shown to be robust to density-based thresholding. Significance. The present work highlights how the inclusion of tractography filtering techniques in connectomic pipelines requires extra caution as they systematically change the network topology both in healthy subjects and patients affected by traumatic brain injury. Finally, the practice of low-to-moderate density-based thresholding of the connectomes is confirmed to have negligible effects on the topological analysis.

Constrained spherical deconvolution of nonspherically sampled diffusion MRI data.

Constrained spherical deconvolution (CSD) of diffusion‐weighted MRI (DW‐MRI) is a popular analysis method that extracts the full white matter (WM) fiber orientation density function (fODF) in the living human brain, noninvasively. It assumes that the DW‐MRI signal on the sphere can be represented as the spherical convolution of a single‐fiber response function (RF) and the fODF, and recovers the fODF through the inverse operation. CSD approaches typically require that the DW‐MRI data is sampled shell‐wise, and estimate the RF in a purely spherical manner using spherical basis functions, such as spherical harmonics (SH), disregarding any radial dependencies. This precludes analysis of data acquired with nonspherical sampling schemes, for example, Cartesian sampling. Additionally, nonspherical sampling can also arise due to technical issues, for example, gradient nonlinearities, resulting in a spatially dependent bias of the apparent tissue densities and connectivity information. Here, we adopt a compact model for the RFs that also describes their radial dependency. We demonstrate that the proposed model can accurately predict the tissue response for a wide range of b‐values. On shell‐wise data, our approach provides fODFs and tissue densities indistinguishable from those estimated using SH. On Cartesian data, fODF estimates and apparent tissue densities are on par with those obtained from shell‐wise data, significantly broadening the range of data sets that can be analyzed using CSD. In addition, gradient nonlinearities can be accounted for using the proposed model, resulting in much more accurate apparent tissue densities and connectivity metrics.

Characterisation of white matter asymmetries in the healthy human brain using diffusion MRI fixel-based analysis

The diffusion tensor model for diffusion MRI has been used extensively to study asymmetry in the human brain white matter. However, given the limitations of the tensor model, the nature of any underlying asymmetries remains uncertain, particularly in crossing fibre regions. Here, we provide a more robust characterisation of human brain white matter asymmetries based on fibre-specific diffusion MRI metrics and a whole-brain data-driven approach. We used high-quality diffusion MRI data (n = 100) from the Human Connectome Project, the spherical deconvolution model for fibre orientation distribution estimation, and the Fixel-Based Analysis framework to utilise crossing fibre information in registration, data smoothing and statistical inference. We found many significant asymmetries, widespread throughout the brain white matter, with both left>right and right>left dominances observed in different pathways. No influences of sex, age, or handedness on asymmetry were found. We also report on the relative contributions of microstructural and morphological white matter properties toward the asymmetry findings. Our findings should provide important information to future studies focussing on how these asymmetries are affected by disease, development/ageing, or how they correlate to functional/cognitive measures.

Assessment of the structural complexity of diffusion MRI voxels using 3D electron microscopy in the rat brain

Validation and interpretation of diffusion magnetic resonance imaging (dMRI) requires detailed understanding of the actual microstructure restricting the diffusion of water molecules. In this study, we used serial block-face scanning electron microscopy (SBEM), a three-dimensional electron microscopy (3D-EM) technique, to image seven white and grey matter volumes in the rat brain. SBEM shows excellent contrast of cellular membranes, which are the major components restricting the diffusion of water in tissue. Additionally, we performed 3D structure tensor (3D-ST) analysis on the SBEM volumes and parameterised the resulting orientation distributions using Watson and angular central Gaussian (ACG) probability distributions as well as spherical harmonic (SH) decomposition. We analysed how these parameterisations described the underlying orientation distributions and compared their orientation and dispersion with corresponding parameters from two dMRI methods, neurite orientation dispersion and density imaging (NODDI) and constrained spherical deconvolution (CSD). Watson and ACG parameterisations and SH decomposition captured well the 3D-ST orientation distributions, but ACG and SH better represented the distributions due to its ability to model asymmetric dispersion. The dMRI parameters corresponded well with the 3D-ST parameters in the white matter volumes, but the correspondence was less evident in the more complex grey matter. SBEM imaging and 3D-ST analysis also revealed that the orientation distributions were often not axially symmetric, a property neatly captured by the ACG distribution. Overall, the ability of SBEM to image diffusion barriers in intricate detail, combined with 3D-ST analysis and parameterisation, provides a step forward toward interpreting and validating the dMRI signals in complex brain tissue microstructure.

Fiber-Specific Changes in White Matter Microstructure in Individuals With X-Linked Auditory Neuropathy.

Auditory neuropathy (AN) is the term used to describe a group of hearing disorders, in which the hearing impairment occurs as a result of abnormal auditory nerve function. While our understanding of this condition has advanced significantly over recent years, the ability to determine the site of lesion and the extent of dysfunction in affected individuals remains a challenge. To this end, we investigated potential axonal degeneration in the white matter tracts of the brainstem in individuals with X-linked AN. We hypothesized that individuals with X-linked AN would show focal degeneration within the VIII nerve and/or auditory brainstem tracts, and the degree of degeneration would correlate with the extent of auditory perceptual impairment. This was achieved using a higher-order diffusion magnetic resonance imaging (dMRI)-based quantitative measure called apparent fiber density as obtained from a technique called single-shell 3-tissue constrained spherical deconvolution and analyzed with the fixel-based analysis framework. Eleven subjects with genetically confirmed X-linked AN and 11 controls with normal hearing were assessed using behavioral and objective auditory measures. dMRI data were also collected for each participant. Fixel-based analysis of the brainstem region showed that subjects with X-linked AN had significantly lower apparent fiber density in the VIII nerve compared with controls, consistent with axonal degeneration in this region. Subsequent analysis of the auditory brainstem tracts specifically showed that degeneration was also significant in these structures overall. The apparent fiber density findings were supported by objective measures of auditory function, such as auditory brainstem responses, electrocochleography, and otoacoustic emissions, which showed VIII nerve activity was severely disrupted in X-linked AN subjects while cochlear sensory hair cell function was relatively unaffected. Moreover, apparent fiber density results were significantly correlated with temporal processing ability (gap detection task) in affected subjects, suggesting that the degree of VIII nerve degeneration may impact the ability to resolve temporal aspects of an acoustic signal. Auditory assessments of sound detection, speech perception, and the processing of binaural cues were also significantly poorer in the X-linked AN group compared with the controls with normal hearing. The results of this study suggest that the dMRI-based measure of apparent fiber density may provide a useful adjunct to existing auditory assessments in the characterization of the site of lesion and extent of dysfunction in individuals with AN. Additionally, the ability to determine the degree of degeneration has the potential to guide rehabilitation strategies in the future.

Diffusion magnetic resonance imaging assessment of regional white matter maturation in preterm neonates

PurposeDiffusion magnetic resonance imaging (dMRI) studies report altered white matter (WM) development in preterm infants. Neurite orientation dispersion and density imaging (NODDI) metrics provide more realistic estimations of neurite architecture in vivo compared with standard diffusion tensor imaging (DTI) metrics. This study investigated microstructural maturation of WM in preterm neonates scanned between 25 and 45 weeks postmenstrual age (PMA) with normal neurodevelopmental outcomes at 2 years using DTI and NODDI metrics.MethodsThirty-one neonates (n = 17 male) with median (range) gestational age (GA) 32+1 weeks (24+2–36+4) underwent 3 T brain MRI at median (range) post menstrual age (PMA) 35+2 weeks (25+3–43+1). WM tracts (cingulum, fornix, corticospinal tract (CST), inferior longitudinal fasciculus (ILF), optic radiations) were delineated using constrained spherical deconvolution and probabilistic tractography in MRtrix3. DTI and NODDI metrics were extracted for the whole tract and cross-sections along each tract to assess regional development.ResultsPMA at scan positively correlated with fractional anisotropy (FA) in the CST, fornix and optic radiations and neurite density index (NDI) in the cingulum, CST and fornix and negatively correlated with mean diffusivity (MD) in all tracts. A multilinear regression model demonstrated PMA at scan influenced all diffusion measures, GA and GAxPMA at scan influenced FA, MD and NDI and gender affected NDI. Cross-sectional analyses revealed asynchronous WM maturation within and between WM tracts.).ConclusionWe describe normal WM maturation in preterm neonates with normal neurodevelopmental outcomes. NODDI can enhance our understanding of WM maturation compared with standard DTI metrics alone.

The effect of gradient nonlinearities on fiber orientation estimates from spherical deconvolution of diffusion magnetic resonance imaging data.

Gradient nonlinearities in magnetic resonance imaging (MRI) cause spatially varying mismatches between the imposed and the effective gradients and can cause significant biases in rotationally invariant diffusion MRI measures derived from, for example, diffusion tensor imaging. The estimation of the orientational organization of fibrous tissue, which is nowadays frequently performed with spherical deconvolution techniques ideally using higher diffusion weightings, can likewise be biased by gradient nonlinearities. We explore the sensitivity of two established spherical deconvolution approaches to gradient nonlinearities, namely constrained spherical deconvolution (CSD) and damped Richardson‐Lucy (dRL). Additionally, we propose an extension of dRL to take into account gradient imperfections, without the need of data interpolation. Simulations show that using the effective b‐matrix can improve dRL fiber orientation estimation and reduces angular deviations, while CSD can be more robust to gradient nonlinearity depending on the implementation. Angular errors depend on a complex interplay of many factors, including the direction and magnitude of gradient deviations, underlying microstructure, SNR, anisotropy of the effective response function, and diffusion weighting. Notably, angular deviations can also be observed at lower b‐values in contrast to the perhaps common assumption that only high b‐value data are affected. In in vivo Human Connectome Project data and acquisitions from an ultrastrong gradient (300 mT/m) scanner, angular differences are observed between applying and not applying the effective gradients in dRL estimation. As even small angular differences can lead to error propagation during tractography and as such impact connectivity analyses, incorporating gradient deviations into the estimation of fiber orientations should make such analyses more reliable.

Preterm birth leads to impaired rich-club organization and fronto-paralimbic/limbic structural connectivity in newborns

Prematurity disrupts brain development during a critical period of brain growth and organization and is known to be associated with an increased risk of neurodevelopmental impairments. Investigating whole-brain structural connectivity alterations accompanying preterm birth may provide a better comprehension of the neurobiological mechanisms related to the later neurocognitive deficits observed in this population.Using a connectome approach, we aimed to study the impact of prematurity on neonatal whole-brain structural network organization at term-equivalent age. In this cohort study, twenty-four very preterm infants at term-equivalent age (VPT-TEA) and fourteen full-term (FT) newborns underwent a brain MRI exam at term age, comprising T2-weighted imaging and diffusion MRI, used to reconstruct brain connectomes by applying probabilistic constrained spherical deconvolution whole-brain tractography. The topological properties of brain networks were quantified through a graph-theoretical approach. Furthermore, edge-wise connectivity strength was compared between groups.Overall, VPT-TEA infants’ brain networks evidenced increased segregation and decreased integration capacity, revealed by an increased clustering coefficient, increased modularity, increased characteristic path length, decreased global efficiency and diminished rich-club coefficient. Furthermore, in comparison to FT, VPT-TEA infants had decreased connectivity strength in various cortico-cortical, cortico-subcortical and intra-subcortical networks, the majority of them being intra-hemispheric fronto-paralimbic and fronto-limbic. Inter-hemispheric connectivity was also decreased in VPT-TEA infants, namely through connections linking to the left precuneus or left dorsal cingulate gyrus – two regions that were found to be hubs in FT but not in VPT-TEA infants. Moreover, posterior regions from Default-Mode-Network (DMN), namely precuneus and posterior cingulate gyrus, had decreased structural connectivity in VPT-TEA group.Our finding that VPT-TEA infants’ brain networks displayed increased modularity, weakened rich-club connectivity and diminished global efficiency compared to FT infants suggests a delayed transition from a local architecture, focused on short-range connections, to a more distributed architecture with efficient long-range connections in those infants. The disruption of connectivity in fronto-paralimbic/limbic and posterior DMN regions might underlie the behavioral and social cognition difficulties previously reported in the preterm population.

Accelerating in vivo fast spin echo high angular resolution diffusion imaging with an isotropic resolution in mice through compressed sensing.

Echo planar imaging (EPI) is commonly used to acquire the many volumes needed for high angular resolution diffusion Imaging (HARDI), posing a higher risk for artifacts, such as distortion and deformation. An alternative to EPI is fast spin echo (FSE) imaging, which has fewer artifacts but is inherently slower. The aim is to accelerate FSE such that a HARDI data set can be acquired in a time comparable to EPI using compressed sensing. Compressed sensing was applied in either q-space or simultaneously in k-space and q-space, by undersampling the k-space in the phase-encoding direction or retrospectively eliminating diffusion directions for different degrees of undersampling. To test the replicability of the acquisition and reconstruction, brain data were acquired from six mice, and a numerical phantom experiment was performed. All HARDI data were analyzed individually using constrained spherical deconvolution, and the apparent fiber density and complexity metric were evaluated, together with whole-brain tractography. The apparent fiber density and complexity metric showed relatively minor differences when only q-space undersampling was used, but deteriorate when k-space undersampling was applied. Likewise, the tract density weighted image showed good results when only q-space undersampling was applied using 15 directions or more, but information was lost when fewer volumes or k-space undersampling were used. It was found that acquiring 15 to 20 diffusion directions with a full k-space and reconstructed using compressed sensing could suffice for a replicable measurement of quantitative measures in mice, where areas near the sinuses and ear cavities are untainted by signal loss.

No Association Between Structural Properties of Corpus Callosum and Handedness: Evidence from the Constrained Spherical Deconvolution Approach

Handedness is the most prominent trait of functional asymmetry in humans, associated with lateralized cognitive functions and considered in relation to mental disorders. However, the neuroanatomical correlates of handedness are still unclear. It has been hypothesized that the structural properties of sub-regions of the corpus callosum (CC) are linked to handedness. Nevertheless, tractography studies of the relation between directly measured structural properties of CC subregions and handedness are lacking. The Constrained Spherical Deconvolution (CSD) approach enables full reconstruction of the sub-regions of the CC. The current study aimed to investigate the relation between the structural properties of the CC, such as volume and the CSD metric, referred to as hindrance modulated orientational anisotropy (HMOA), and handedness. Handedness was considered in two dimensions: direction (right-handed, ambidextrous, left-handed) and degree (the absolute values of Handedness quotient). We found no association between 1) volume or HMOA as a proxy of microstructural properties, namely the axonal diameter and fiber dispersion, of each sub-region and 2) either the direction or the degree of handedness. These findings suggest the absence of a direct relation between sub-regions of the CC and handedness, demonstrating the necessity of future tractography studies.

Spherical deconvolution with tissue-specific response functions and multi-shell diffusion MRI to estimate multiple fiber orientation distributions (mFODs)

In diffusion MRI, spherical deconvolution approaches can estimate local white matter (WM) fiber orientation distributions (FOD) which can be used to produce fiber tractography reconstructions. The applicability of spherical deconvolution to gray matter (GM), however, is still limited, despite its critical role as start/endpoint of WM fiber pathways. The advent of multi-shell diffusion MRI data offers additional contrast to model the GM signal but, to date, only isotropic models have been applied to GM. Evidence from both histology and high-resolution diffusion MRI studies suggests a marked anisotropic character of the diffusion process in GM, which could be exploited to improve the description of the cortical organization. In this study, we investigated whether performing spherical deconvolution with tissue specific models of both WM and GM can improve the characterization of the latter while retaining state-of-the-art performances in WM. To this end, we developed a framework able to simultaneously accommodate multiple anisotropic response functions to estimate multiple, tissue-specific, fiber orientation distributions (mFODs). As proof of principle, we used the diffusion kurtosis imaging model to represent the WM signal, and the neurite orientation dispersion and density imaging (NODDI) model to represent the GM signal. The feasibility of the proposed approach is shown with numerical simulations and with data from the Human Connectome Project (HCP). The performance of our method is compared to the current state of the art, multi-shell constrained spherical deconvolution (MSCSD). The simulations show that with our new method we can accurately estimate a mixture of two FODs at SNR≥50. With HCP data, the proposed method was able to reconstruct both tangentially and radially oriented FODs in GM, and performed comparably well to MSCSD in computing FODs in WM. When performing fiber tractography, the trajectories reconstructed with mFODs reached the cortex with more spatial continuity and for a longer distance as compared to MSCSD and allowed to reconstruct short trajectories tangential to the cortical folding. In conclusion, we demonstrated that our proposed method allows to perform spherical deconvolution of multiple anisotropic response functions, specifically improving the performances of spherical deconvolution in GM tissue.

Characterizing White Matter Tract Organization in Polymicrogyria and Lissencephaly: A Multifiber Diffusion MRI Modeling and Tractography Study.

Polymicrogyria and lissencephaly may be associated with abnormal organization of the undelying white matter tracts that have been rarely investigated so far. Our aim was to characterize white matter tract organization in polymicrogyria and lissencephaly using constrained spherical deconvolution, a multifiber diffusion MR imaging modeling technique for white matter tractography reconstruction. We retrospectively reviewed 50 patients (mean age, 8.3 ± 5.4 years; range, 1.4-21.2 years; 27 males) with different polymicrogyria (n = 42) and lissencephaly (n = 8) subtypes. The fiber direction-encoded color maps and 6 different white matter tracts reconstructed from each patient were visually compared with corresponding images reconstructed from 7 age-matched, healthy control WM templates. Each white matter tract was assessed by 2 experienced pediatric neuroradiologists and scored in consensus on the basis of the severity of the structural abnormality, ranging from the white matter tracts being absent to thickened. The results were summarized by different polymicrogyria and lissencephaly subgroups. More abnormal-appearing white matter tracts were identified in patients with lissencephaly compared with those with polymicrogyria (79.2% versus 37.3%). In lissencephaly, structural abnormalities were identified in all studied white matter tracts. In polymicrogyria, the more frequently affected white matter tracts were the cingulum, superior longitudinal fasciculus, inferior longitudinal fasciculus, and optic radiation-posterior corona radiata. The severity of superior longitudinal fasciculus and cingulum abnormalities was associated with the polymicrogyria distribution and extent. A thickened superior fronto-occipital fasciculus was demonstrated in 3 patients. We demonstrated a range of white matter tract structural abnormalities in patients with polymicrogyria and lissencephaly. The patterns of white matter tract involvement are related to polymicrogyria and lissencephaly subgroups, distribution, and, possibly, their underlying etiologies.

Optimized rectification of fiber orientation density function.

To demonstrate an optimized rectification strategy for fiber orientation density functions (fODFs). In white matter, fODFs can be estimated with diffusion MRI. However, because of signal noise, imaging artifacts and other factors, experimentally determined fODFs may take on unphysical negative values in some directions. Here, we show how to rectify such fODFs to eliminate all negative values while minimizing the mean square difference between the original and rectified fODFs. The method is demonstrated for a mathematical model and for fODFs estimated from experimental human data using both constrained spherical deconvolution and fiber ball imaging. Comparison with an alternative nonoptimized rectification approach is also provided. For the mathematical model, it is found that the optimized rectification procedure removes negative fODF values while at the same time reducing the mean square error. Relative to the alternative rectification approach, the optimized fODFs are substantially more accurate. For the experimental data, the optimized fODFs have a lower average fractional anisotropy axonal and often fewer small peaks than the original, unrectified fODFs. The calculation of optimized fODFs is straightforward where the main step is the finding of the root to an equation in one variable, as may be efficiently accomplished with the bisection method. Unphysical negative fODF values can be easily eliminated in a manner that minimizes the mean square difference between the original and rectified fODFs. Optimized fODF rectification may be useful in applications for which negative values are problematic.

Maturation and interhemispheric asymmetry in neurite density and orientation dispersion in early childhood

BackgroundThe brain's white matter undergoes profound changes during early childhood, which are believed to underlie the rapid development of cognitive and behavioral skills during this period. Neurite density, and complexity of axonal projections, have been shown to change across the life span, though changes during early childhood are poorly characterized. Here, we utilize neurite orientation dispersion and density imaging (NODDI) to investigate maturational changes in tract-wise neurite density index (NDI) and orientation dispersion index (ODI) during early childhood. Additionally, we assess hemispheric asymmetry of tract-wise NDI and ODI values, and longitudinal changes. MethodsTwo sets of diffusion weighted images with different diffusion-weighting were collected from 125 typically developing children scanned at baseline (N = 125; age range = 4.14–7.29; F/M = 73/52), 6-month (N = 8; F/M = 8/0), and 12-month (N = 52; F/M = 39/13) timepoints. NODDI and template-based tractography using constrained spherical deconvolution were utilized to calculate NDI and ODI values for major white matter tracts. Mixed-effects models controlling for sex, handedness, and in-scanner head motion were utilized to assess developmental changes in tract-wise NDI and ODI. Additional mixed-effects models were used to assess interhemispheric differences in tract-wise NDI and ODI values and hemispheric asymmetries in tract-wise development. ResultsMaturational increases in NDI were seen in all major white matter tracts, though we did not observe the expected tract-wise pattern of maturational rates (e.g. fast commissural/projection and slow frontal/temporal tract change). ODI did not change significantly with age in any tract. We observed greater NDI and ODI values in the right as compared to the left hemisphere for most tracts, but no hemispheric asymmetry for rate of change with age. ConclusionsThese findings suggest that neurite density, but not orientation dispersion, increases with age during early childhood. In relation to NDI growth trends reported in infancy and late-childhood, our results suggest that early childhood may be a transitional period for neurite density maturation wherein commissural/projection fibers are approaching maturity, maturation in long range association fibers is increasing, and changes in limbic/frontal fibers remain modest. Rightward asymmetry in NDI and ODI values, but no asymmetry in developmental changes, suggests that rightward asymmetry of neurite density and orientation dispersion is established prior to age 4.