Diffusion MRI fiber tractography of the brain

Towards linking diffusion MRI based macro- and microstructure measures with cortico-cortical transmission in brain tumor patients

We aimed to link macro- and microstructure measures of brain white matter obtained from diffusion MRI with effective connectivity measures based on a propagation of cortico-cortical evoked potentials induced with intrasurgical direct electrical stimulation. For this, we compared streamline lengths and log-transformed ratios of streamlines computed from presurgical diffusion-weighted images, and the delays and amplitudes of N1 peaks recorded intrasurgically with electrocorticography electrodes in a pilot study of 9 brain tumor patients. Our results showed positive correlation between these two modalities in the vicinity of the stimulation sites (Pearson coefficient 0.54±0.13 for N1 delays, and 0.47±0.23 for N1 amplitudes), which could correspond to the neural propagation via U-fibers. In addition, we reached high sensitivities (0.78±0.07) and very high specificities (0.93±0.03) in a binary variant of our comparison. Finally, we used the structural connectivity measures to predict the effective connectivity using a multiple linear regression model, and showed a significant role of brain microstructure-related indices in this relation.

Structural connectome and connectivity lateralization of the multimodal vestibular cortical network

Unlike other sensory systems, the structural connectivity patterns of the human vestibular cortex remain a matter of debate. Based on their functional properties and hypothesized centrality within the vestibular network, the ‘core’ cortical regions of this network are thought to be areas in the posterior peri-sylvian cortex, in particular the retro-insula (previously named the posterior insular cortex-PIC), and the subregion OP2 of the parietal operculum.

To study the vestibular network, structural connectivity matrices from n=974 healthy individuals drawn from the public Human Connectome Project (HCP) repository were estimated using multi-shell diffusion-weighted data followed by probabilistic tractography and spherical-deconvolution informed filtering of tractograms in combination with subject-specific grey-matter parcellations. Weighted graph-theoretical measures, modularity, and ‘hubness’ of the multimodal vestibular network were then estimated, and a structural lateralization index was defined in order to assess the difference in fiber density of homonym regions in the right and left hemisphere. Differences in connectivity patterns between OP2 and PIC were also estimated.

We found that the bilateral intraparietal sulcus, PIC, and to a lesser degree OP2, are key ‘hub’ regions within the multimodal vestibular network. PIC and OP2 structural connectivity patterns were lateralized to the left hemisphere, while structural connectivity patterns of the posterior peri-sylvian supramarginal and superior temporal gyri were lateralized to the right hemisphere. These lateralization patterns were independent of handedness.

We also found that the structural connectivity pattern of PIC is consistent with a key role of PIC in visuo-vestibular processing and that the structural connectivity pattern of OP2 is consistent with integration of mainly vestibular somato-sensory and motor information. These results suggest an analogy between PIC and the simian visual posterior sylvian (VPS) area and OP2 and the simian parieto-insular vestibular cortex (PIVC).

Overall, these findings may provide novel insights to the current models of vestibular function, as well as to the understanding of the complexity and lateralized signs of vestibular syndromes.

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.

Entrack: Probabilistic Spherical Regression with Entropy Regularization for Fiber Tractography

White matter tractography, based on diffusion-weighted magnetic resonance images, is currently the only available in vivo method to gather information on the structural brain connectivity. The low resolution of diffusion MRI data suggests to employ probabilistic methods for streamline reconstruction, i.e., for fiber crossings. We propose a general probabilistic model for spherical regression based on the Fisher-von-Mises distribution, which efficiently estimates maximum entropy posteriors of local streamline directions with machine learning methods. The optimal precision of posteriors for streamlines is determined by an information-theoretic technique, the expected log-posterior agreement concept. It relies on the requirement that the posterior distributions of streamlines, inferred on retest measurements of the same subject, should yield stable results within the precision determined by the noise level of the data source.

Reproducing macaque lateral grasping and oculomotor networks using resting state functional connectivity and diffusion tractography

Cortico-cortical networks involved in motor control have been well defined in the macaque using a range of invasive techniques. The advent of neuroimaging has enabled non-invasive study of these large-scale functionally specialized networks in the human brain; however, assessing its accuracy in reproducing genuine anatomy is more challenging. We set out to assess the similarities and differences between connections of macaque motor control networks defined using axonal tracing and those reproduced using structural and functional connectivity techniques. We processed a cohort of macaques scanned in vivo that were made available by the open access PRIME-DE resource, to evaluate connectivity using diffusion imaging tractography and resting state functional connectivity (rs-FC). Sectors of the lateral grasping and exploratory oculomotor networks were defined anatomically on structural images, and connections were reproduced using different structural and functional approaches (probabilistic and deterministic whole-brain and seed-based tractography; group template and native space functional connectivity analysis). The results showed that parieto-frontal connections were best reproduced using both structural and functional connectivity techniques. Tractography showed lower sensitivity but better specificity in reproducing connections identified by tracer data. Functional connectivity analysis performed in native space had higher sensitivity but lower specificity and was better at identifying connections between intrasulcal ROIs than group-level analysis. Connections of AIP were most consistently reproduced, although those connected with prefrontal sectors were not identified. We finally compared diffusion MR modelling with histology based on an injection in AIP and speculate on anatomical bases for the observed false negatives. Our results highlight the utility of precise ex vivo techniques to support the accuracy of neuroimaging in reproducing connections, which is relevant also for human studies.

Opportunities of connectomic neuromodulation

The process of altering neural activity - neuromodulation - has long been used to treat patients with brain disorders and answer scientific questions. Deep brain stimulation in particular has provided clinical benefit to over 150,000 patients. However, our understanding of how neuromodulation impacts the brain is evolving. Instead of focusing on the local impact at the stimulation site itself, we are considering the remote impact on brain regions connected to the stimulation site. Brain connectivity information derived from advanced magnetic resonance imaging data can be used to identify these connections and better understand clinical and behavioral effects of neuromodulation. In this article, we review studies combining neuromodulation and brain connectomics, highlighting opportunities where this approach may prove particularly valuable. We focus on deep brain stimulation, but show that the same principles can be applied to other forms of neuromodulation, such as transcranial magnetic stimulation and MRI-guided focused ultrasound. We outline future perspectives and provide testable hypotheses for future work.

Fine-Grained Topography and Modularity of the Macaque Frontal Pole Cortex Revealed by Anatomical Connectivity Profiles

The frontal pole cortex (FPC) plays key roles in various higher-order functions and is highly developed in non-human primates. An essential missing piece of information is the detailed anatomical connections for finer parcellation of the macaque FPC than provided by the previous tracer results. This is important for understanding the functional architecture of the cerebral cortex. Here, combining cross-validation and principal component analysis, we formed a tractography-based parcellation scheme that applied a machine learning algorithm to divide the macaque FPC (2 males and 6 females) into eight subareas using high-resolution diffusion magnetic resonance imaging with the 9.4T Bruker system, and then revealed their subregional connections. Furthermore, we applied improved hierarchical clustering to the obtained parcels to probe the modular structure of the subregions, and found that the dorsolateral FPC, which contains an extension to the medial FPC, was mainly connected to regions of the default-mode network. The ventral FPC was mainly involved in the social-interaction network and the dorsal FPC in the metacognitive network. These results enhance our understanding of the anatomy and circuitry of the macaque brain, and contribute to FPC-related clinical research.

Tractography-Based Analysis of Morphological and Anatomical Characteristics of the Uncinate Fasciculus in Human Brains

(1) Background: The uncinate fasciculus (UF) is a white matter bundle connecting the prefrontal cortex and temporal lobe. The functional role of the uncinate fasciculus is still uncertain. The role of the UF is attributed to the emotional empathy network. The present study aimed to more accurately the describe anatomical variability of the UF by focusing on the volume of fibers and testing for correlations with sex and age. (2) Material and Methods: Magnetic resonance imaging of adult patients with diffusion tensor imaging (DTI) was performed on 34 patients. The total number of fibers, volume of UF, and number of tracts were processed using DSI studio software. The DSI studio allows for mapping of different nerve pathways and visualizing of the obtained results using spatial graphics. (3) Results: The total number of UF tracts was significantly higher in the right hemisphere compared to the left hemisphere (right M ± SD = 52 ± 24; left: 39 ± 25, p < 0.05). A hook-shaped UF was the most common variant (91.7%). The UF volumes were larger in men (1410 ± 150.7 mm³) as compared to women (1325 ± 133.2 mm³) (p < 0.05). The mean fractional anisotropy (FA) values of the UF were significantly larger on the left side 0.597, while the right UF had an average of 0.346 (p < 0.05). Patients older than 50 years old had a significantly higher value of mean diffusivity (MD) (p = 0.034). In 73.5% of patients, a greater number of fibers terminated in the inferior part of the inferior frontal gyrus. (4) Conclusions: The morphological characteristics of the UF, unlike the shape, are associated with sex and are characterized by hemispheric dominance. These findings confirm the results of the previous studies. Future research should examine the potential correlation among the UF volume, number of fibers, and total brain volume in both sexes and patient psychological state.

Brain Metabolism and Structure in Chronic Migraine

PURPOSE OF REVIEW

The purpose of this paper is to review and synthesize current literature in which neurochemical and structural brain imaging were used to investigate chronic migraine (CM) pathophysiology and to further discuss the clinical implications.

RECENT FINDINGS

Spectroscopic and structural MRI studies have shown the presence of both impaired metabolism and structural alterations in the brain of CM patients. Metabolic changes in key brain regions support the notion of altered energetics and homeostasis as part of CM pathophysiology. Furthermore, CM, like other chronic pain disorders, may undergo structural reorganization in pain-related brain regions following near persistent endogenous painful input. Finally, both imaging techniques may provide potential biomarkers of disease state and progression and may help guide novel therapeutic interventions or strategies. Spectroscopic and structural MRI have revealed novel aspects of CM pathophysiology. Findings from the former support the metabolic theory of migraine pathogenesis.

Topographic organization of the human subcortex unveiled with functional connectivity gradients

Brain atlases are fundamental to understanding the topographic organization of the human brain, yet many contemporary human atlases cover only the cerebral cortex, leaving the subcortex a terra incognita. We use functional MRI (fMRI) to map the complex topographic organization of the human subcortex, revealing large-scale connectivity gradients and new areal boundaries. We unveil four scales of subcortical organization that recapitulate well-known anatomical nuclei at the coarsest scale and delineate 27 new bilateral regions at the finest. Ultrahigh field strength fMRI corroborates and extends this organizational structure, enabling the delineation of finer subdivisions of the hippocampus and the amygdala, while task-evoked fMRI reveals a subtle subcortical reorganization in response to changing cognitive demands. A new subcortical atlas is delineated, personalized to represent individual differences and used to uncover reproducible brain-behavior relationships. Linking cortical networks to subcortical regions recapitulates a task-positive to task-negative axis. This new atlas enables holistic connectome mapping and characterization of cortico-subcortical connectivity.

Classifyber, a robust streamline-based linear classifier for white matter bundle segmentation

Virtual delineation of white matter bundles in the human brain is of paramount importance for multiple applications, such as pre-surgical planning and connectomics. A substantial body of literature is related to methods that automatically segment bundles from diffusion Magnetic Resonance Imaging (dMRI) data indirectly, by exploiting either the idea of connectivity between regions or the geometry of fiber paths obtained with tractography techniques, or, directly, through the information in volumetric data. Despite the remarkable improvement in automatic segmentation methods over the years, their segmentation quality is not yet satisfactory, especially when dealing with datasets with very diverse characteristics, such as different tracking methods, bundle sizes or data quality. In this work, we propose a novel, supervised streamline-based segmentation method, called Classifyber, which combines information from atlases, connectivity patterns, and the geometry of fiber paths into a simple linear model. With a wide range of experiments on multiple datasets that span from research to clinical domains, we show that Classifyber substantially improves the quality of segmentation as compared to other state-of-the-art methods and, more importantly, that it is robust across very diverse settings. We provide an implementation of the proposed method as open source code, as well as web service.

Structure-function abnormalities in cortical sensory projections in embouchure dystonia

BACKGROUND

Embouchure dystonia (ED) is a task-specific focal dystonia in professional brass players leading to abnormal orofacial muscle posturing/spasms during performance. Previous studies have outlined abnormal cortical sensorimotor function during sensory/motor tasks and in the resting state as well as abnormal cortical sensorimotor structure. Yet, potentially underlying white-matter tract abnormalities in this network disease are unknown.

OBJECTIVE

To delineate structure-function abnormalities within cerebral sensorimotor trajectories in ED.

METHOD

Probabilistic tractography and seed-based functional connectivity analysis were performed in 16/16 ED patients/healthy brass players within a simple literature-informed network model of cortical sensorimotor processing encompassing supplementary motor, superior parietal, primary somatosensory and motor cortex as well as the putamen. Post-hoc grey matter volumetry was performed within cortices of abnormal trajectories.

RESULTS

ED patients showed average axial diffusivity reduction within projections between the primary somatosensory cortex and putamen, with converse increases within projections between supplementary motor and superior parietal cortex in both hemispheres. Increase in the mode of anisotropy in patients was accompanying the latter left-hemispheric projection, as well as in the supplementary motor area's projection to the left primary motor cortex. Patient's left primary somatosensory functional connectivity with the putamen was abnormally reduced and significantly associated with the axial diffusivity reduction. Left primary somatosensory grey matter volume was increased in patients.

CONCLUSION

Correlates of abnormal tract integrity within primary somatosensory cortico-subcortical projections and higher-order sensorimotor projections support the key role of dysfunctional sensory information propagation in ED pathophysiology. Differential directionality of cortico-cortical and cortico-subcortical abnormalities hints at non-uniform sensory system changes.

Structural Network Analysis Using Diffusion MRI Tractography in Parkinson's Disease and Correlations With Motor Impairment

Functional impairment of spatially distributed brain regions in Parkinson's disease (PD) suggests changes in integrative and segregative network characteristics, for which novel analysis methods are available. To assess underlying structural network differences between PD patients and controls, we employed MRI T1 gray matter segmentation and diffusion MRI tractography to construct connectivity matrices to compare patients and controls with data originating from two different centers. In the Dutch dataset (Data-NL), 14 PD patients, and 15 healthy controls were analyzed, while 19 patients and 18 controls were included in the Canadian dataset (Data-CA). All subjects underwent T1 and diffusion-weighted MRI. Patients were assessed with Part 3 of the Unified Parkinson's Disease Rating Scale (UPDRS). T1 images were segmented using FreeSurfer, while tractography was performed using ExploreDTI. The regions of interest from the FreeSurfer segmentation were combined with the white matter streamline sets resulting from the tractography, to construct connectivity matrices. From these matrices, both global and local efficiencies were calculated, which were compared between the PD and control groups and related to the UPDRS motor scores. The connectivity matrices showed consistent patterns among the four groups, without significant differences between PD patients and control subjects, either in Data-NL or in Data-CA. In Data-NL, however, global and local efficiencies correlated negatively with UPDRS scores at both the whole-brain and the nodal levels [false discovery rate (FDR) 0.05]. At the nodal level, particularly, the posterior parietal cortex showed a negative correlation between UPDRS and local efficiency, while global efficiency correlated negatively with the UPDRS in the sensorimotor cortex. The spatial patterns of negative correlations between UPDRS and parameters for network efficiency seen in Data-NL suggest subtle structural differences in PD that were below sensitivity thresholds in Data-CA. These correlations are in line with previously described functional differences. The methodological approaches to detect such differences are discussed.

New Cognitive Neurotechnology Facilitates Studies of Cortical-Subcortical Interactions

Most of the studies employing neuroimaging have focused on cortical and subcortical signals individually to obtain neurophysiological signatures of cognitive functions. However, understanding the dynamic communication between the cortex and subcortical structures is essential for unraveling the neural correlates of cognition. In this quest, magnetoencephalography (MEG) and electroencephalography (EEG) are the methods of choice because they are noninvasive electrophysiological recording techniques with high temporal resolution. Sophisticated MEG/EEG source estimation techniques and network analysis methods, developed recently, can provide a more comprehensive understanding of the neurophysiological mechanisms of fundamental cognitive processes. Used together with noninvasive modulation of cortical-subcortical communication, these approaches may open up new possibilities for expanding the repertoire of noninvasive cognitive neurotechnology.

Recent advances of deep learning in psychiatric disorders

Deep learning (DL) is a recently proposed subset of machine learning methods that has gained extensive attention in the academic world, breaking benchmark records in areas such as visual recognition and natural language processing. Different from conventional machine learning algorithm, DL is able to learn useful representations and features directly from raw data through hierarchical nonlinear transformations. Because of its ability to detect abstract and complex patterns, DL has been used in neuroimaging studies of psychiatric disorders, which are characterized by subtle and diffuse alterations. Here, we provide a brief review of recent advances and associated challenges in neuroimaging studies of DL applied to psychiatric disorders. The results of these studies indicate that DL could be a powerful tool in assisting the diagnosis of psychiatric diseases. We conclude our review by clarifying the main promises and challenges of DL application in psychiatric disorders, and possible directions for future research.

Brain Characteristics Noted Prior to and Following Cranial Orthotic Treatment

Objective:

This case report aims to assess a potential association between cranial asymmetry, brain deformation, and associated developmental delay.

Study Design:

Two infants born at ≥37 weeks pursuing cranial orthotic treatment for severe Deformational Plagiocephaly (DP) (cranial vault asymmetry index >8.75%) underwent developmental assessment using Mullen Scales of Early Learning (MSEL) and non-sedated brain structural and diffusion magnetic resonance imaging (MRI) prior to and following cranial orthotic treatment.

Results:

In both infants with DP, tractography results revealed alterations in the white matter pathways of the brain. Both infants also had low to low/normal visual receptivity and fine motor skills. After cranial orthotic treatment, cranial asymmetry improved but did not completely resolve, tractography demonstrated a change toward normalized white matter pathways, and visual receptivity and fine motor skills improved.

Conclusions:

These preliminary findings suggest a potential link between DP, altered brain structures, and developmental assessment. Further investigation with a larger sample is warranted.

NODDI in clinical research

Diffusion MRI (dMRI) has proven to be a useful imaging approach for both clinical diagnosis and research investigating the microstructures of nervous tissues, and it has helped us to better understand the neurophysiological mechanisms of many diseases. Though diffusion tensor imaging (DTI) has long been the default tool to analyze dMRI data in clinical research, acquisition with stronger diffusion weightings beyond the DTI regimen is now possible with modern clinical scanners, potentially enabling even more detailed characterization of tissue microstructures. To take advantage of such data, neurite orientation dispersion and density imaging (NODDI) has been proposed as a way to relate the dMRI signal to tissue features via biophysically inspired modeling. The number of reports demonstrating the potential clinical utility of NODDI is rapidly increasing. At the same time, the pitfalls and limitations of NODDI, and general challenges in microstructure modeling, are becoming increasingly recognized by clinicians. dMRI microstructure modeling is a rapidly evolving field with great promise, where people from different scientific backgrounds, such as physics, medicine, biology, neuroscience, and statistics, are collaborating to build novel tools that contribute to improving human healthcare. Here, we review the applications of NODDI in clinical research and discuss future perspectives for investigations toward the implementation of dMRI microstructure imaging in clinical practice.

The efficacy of different preprocessing steps in reducing motion-related confounds in diffusion MRI connectomics

Head motion is a major confounding factor in neuroimaging studies. While numerous studies have investigated how motion impacts estimates of functional connectivity, the effects of motion on structural connectivity measured using diffusion MRI have not received the same level of attention, despite the fact that, like functional MRI, diffusion MRI relies on elaborate preprocessing pipelines that require multiple choices at each step. Here, we report a comprehensive analysis of how these choices influence motion-related contamination of structural connectivity estimates. Using a healthy adult sample (N = 294), we evaluated 240 different preprocessing pipelines, devised using plausible combinations of different choices related to explicit head motion correction, tractography propagation algorithms, track seeding methods, track termination constraints, quantitative metrics derived for each connectome edge, and parcellations. We found that an approach to motion correction that includes outlier replacement and within-slice volume correction led to a dramatic reduction in cross-subject correlations between head motion and structural connectivity strength, and that motion contamination is more severe when quantifying connectivity strength using mean tract fractional anisotropy rather than streamline count. We also show that the choice of preprocessing strategy can significantly influence subsequent inferences about network organization, with the location of network hubs varying considerably depending on the specific preprocessing steps applied. Our findings indicate that the impact of motion on structural connectivity can be successfully mitigated using recent motion-correction algorithms that include outlier replacement and within-slice motion correction.

Using structural connectivity to augment community structure in EEG functional connectivity

Recently, EEG recording techniques and source analysis have improved, making it feasible to tap into fast network dynamics. Yet, analyzing whole-cortex EEG signals in source space is not standard, partly because EEG suffers from volume conduction: Functional connectivity (FC) reflecting genuine functional relationships is impossible to disentangle from spurious FC introduced by volume conduction. Here, we investigate the relationship between white matter structural connectivity (SC) and large-scale network structure encoded in EEG-FC. We start by confirming that FC (power envelope correlations) is predicted by SC beyond the impact of Euclidean distance, in line with the assumption that SC mediates genuine FC. We then use information from white matter structural connectivity in order to smooth the EEG signal in the space spanned by graphs derived from SC. Thereby, FC between nearby, structurally connected brain regions increases while FC between nonconnected regions remains unchanged, resulting in an increase in genuine, SC-mediated FC. We analyze the induced changes in FC, assessing the resemblance between EEG-FC and volume-conduction- free fMRI-FC, and find that smoothing increases resemblance in terms of overall correlation and community structure. This result suggests that our method boosts genuine FC, an outcome that is of interest for many EEG network neuroscience questions.

The Importance of Cerebellar Connectivity on Simulated Brain Dynamics

The brain shows a complex multiscale organization that prevents a direct understanding of how structure, function and dynamics are correlated. To date, advances in neural modeling offer a unique opportunity for simulating global brain dynamics by embedding empirical data on different scales in a mathematical framework. The Virtual Brain (TVB) is an advanced data-driven model allowing to simulate brain dynamics starting from individual subjects' structural and functional connectivity obtained, for example, from magnetic resonance imaging (MRI). The use of TVB has been limited so far to cerebral connectivity but here, for the first time, we have introduced cerebellar nodes and interconnecting tracts to demonstrate the impact of cerebro-cerebellar loops on brain dynamics. Indeed, the matching between the empirical and simulated functional connectome was significantly improved when including the cerebro-cerebellar loops. This positive result should be considered as a first step, since issues remain open about the best strategy to reconstruct effective structural connectivity and the nature of the neural mass or mean-field models generating local activity in the nodes. For example, signal processing is known to differ remarkably between cortical and cerebellar microcircuits. Tackling these challenges is expected to further improve the predictive power of functional brain activity simulations, using TVB or other similar tools, in explaining not just global brain dynamics but also the role of cerebellum in determining brain states in physiological conditions and in the numerous pathologies affecting the cerebro-cerebellar loops.

Imaging Alzheimer's genetic risk using diffusion MRI: A systematic review

Diffusion magnetic resonance imaging (dMRI) is an imaging technique which probes the random motion of water molecules in tissues and has been widely applied to investigate changes in white matter microstructure in Alzheimer's Disease. This paper aims to systematically review studies that examined the effect of Alzheimer's risk genes on white matter microstructure. We assimilated findings from 37 studies and reviewed their diffusion pre-processing and analysis methods. Most studies estimate the diffusion tensor (DT) and compare derived quantitative measures such as fractional anisotropy and mean diffusivity between groups. Those with increased AD genetic risk are associated with reduced anisotropy and increased diffusivity across the brain, most notably the temporal and frontal lobes, cingulum and corpus callosum. Structural abnormalities are most evident amongst those with established Alzheimer's Disease. Recent studies employ signal representations and analysis frameworks beyond DT MRI but show that dMRI overall lacks specificity to disease pathology. However, as the field advances, these techniques may prove useful in pre-symptomatic diagnosis or staging of Alzheimer's disease.

Non-invasive inference of information flow using diffusion MRI, functional MRI, and MEG

OBJECTIVE

To infer information flow in the white matter of the brain and recover cortical activity using functional MRI, diffusion MRI, and MEG without a manual selection of the white matter connections of interest.

APPROACH

A Bayesian network which encodes the priors knowledge of possible brain states is built from imaging data. Diffusion MRI is used to enumerate all possible connections between cortical regions. Functional MRI is used to prune connections without manual intervention and increase the likelihood of specific regions being active. MEG data is used as evidence into this network to obtain a posterior distribution on cortical regions and connections.

MAIN RESULTS

We show that our proposed method is able to identify connections associated with the a sensory-motor task. This allows us to build the Bayesian network with no manual selection of connections of interest. Using sensory-motor MEG evoked response as evidence into this network, our method identified areas known to be involved in a visuomotor task. In addition, information flow along white matter fiber bundles connecting those regions was also recovered.

SIGNIFICANCE

Current methods to estimate white matter information flow are extremely invasive, therefore limiting our understanding of the interaction between cortical regions. The proposed method makes use of functional MRI, diffusion MRI, and M/EEG to infer communication between cortical regions, therefore opening the door to the non-invasive exploration of information flow in the white matter.

Analytical and fast Fiber Orientation Distribution reconstruction in 3D-Polarized Light Imaging

Three dimensional Polarized Light Imaging (3D-PLI) is an optical technique which allows mapping the spatial fiber architecture of fibrous postmortem tissues, at sub-millimeter resolutions. Here, we propose an analytical and fast approach to compute the fiber orientation distribution (FOD) from high-resolution vector data provided by 3D-PLI. The FOD is modeled as a sum of K orientations/Diracs on the unit sphere, described on a spherical harmonics basis and analytically computed using the spherical Fourier transform. Experiments are performed on rich synthetic data which simulate the geometry of the neuronal fibers and on human brain data. Results indicate the analytical FOD is computationally efficient and very fast, and has high angular precision and angular resolution. Furthermore, investigations on the right occipital lobe illustrate that our strategy of FOD computation enables the bridging of spatial scales from microscopic 3D-PLI information to macro- or mesoscopic dimensions of diffusion Magnetic Resonance Imaging (MRI), while being a means to evaluate prospective resolution limits for diffusion MRI to reconstruct region-specific white matter tracts. These results demonstrate the interest and great potential of our analytical approach.

Connections of the Human Orbitofrontal Cortex and Inferior Frontal Gyrus

The direct connections of the orbitofrontal cortex (OFC) were traced with diffusion tractography imaging and statistical analysis in 50 humans, to help understand better its roles in emotion and its disorders. The medial OFC and ventromedial prefrontal cortex have direct connections with the pregenual and subgenual parts of the anterior cingulate cortex; all of which are reward-related areas. The lateral OFC (OFClat) and its closely connected right inferior frontal gyrus (rIFG) have direct connections with the supracallosal anterior cingulate cortex; all of which are punishment or nonreward-related areas. The OFClat and rIFG also have direct connections with the right supramarginal gyrus and inferior parietal cortex, and with some premotor cortical areas, which may provide outputs for the OFClat and rIFG. Another key finding is that the ventromedial prefrontal cortex shares with the medial OFC especially strong outputs to the nucleus accumbens and olfactory tubercle, which comprise the ventral striatum, whereas the other regions have more widespread outputs to the striatum. Direct connections of the OFC and IFG were with especially the temporal pole part of the temporal lobe. The left IFG, which includes Broca's area, has direct connections with the left angular and supramarginal gyri.

A half-century of innovation in technology-preparing MRI for the 21st century

MRI developed during the last half-century from a very basic concept to an indispensable non-ionising medical imaging technique that has found broad application in diagnostics, therapy control and far beyond. Due to its excellent soft-tissue contrast and the huge variety of accessible tissue- and physiological-parameters, MRI is often preferred to other existing modalities. In the course of its development, MRI underwent many substantial transformations. From the beginning, starting as a proof of concept, much effort was expended to develop the appropriate basic scanning technology and methodology, and to establish the many clinical contrasts (e.g., T₁, T₂, flow, diffusion, water/fat, etc.) that MRI is famous for today. Beyond that, additional prominent innovations to the field have been parallel imaging and compressed sensing, leading to significant scanning time reductions, and the move towards higher static magnetic field strengths, which led to increased sensitivity and improved image quality. Improvements in workflow and the use of artificial intelligence are among many current trends seen in this field, paving the way for a broad use of MRI. The 125th anniversary of the BJR is a good point to reflect on all these changes and developments and to offer some slightly speculative ideas as to what the future may bring.

Tractostorm: The what, why, and how of tractography dissection reproducibility

Investigative studies of white matter (WM) brain structures using diffusion MRI (dMRI) tractography frequently require manual WM bundle segmentation, often called "virtual dissection." Human errors and personal decisions make these manual segmentations hard to reproduce, which have not yet been quantified by the dMRI community. It is our opinion that if the field of dMRI tractography wants to be taken seriously as a widespread clinical tool, it is imperative to harmonize WM bundle segmentations and develop protocols aimed to be used in clinical settings. The EADC-ADNI Harmonized Hippocampal Protocol achieved such standardization through a series of steps that must be reproduced for every WM bundle. This article is an observation of the problematic. A specific bundle segmentation protocol was used in order to provide a real-life example, but the contribution of this article is to discuss the need for reproducibility and standardized protocol, as for any measurement tool. This study required the participation of 11 experts and 13 nonexperts in neuroanatomy and "virtual dissection" across various laboratories and hospitals. Intra-rater agreement (Dice score) was approximately 0.77, while inter-rater was approximately 0.65. The protocol provided to participants was not necessarily optimal, but its design mimics, in essence, what will be required in future protocols. Reporting tractometry results such as average fractional anisotropy, volume or streamline count of a particular bundle without a sufficient reproducibility score could make the analysis and interpretations more difficult. Coordinated efforts by the diffusion MRI tractography community are needed to quantify and account for reproducibility of WM bundle extraction protocols in this era of open and collaborative science.

Asymmetric fiber trajectory distribution estimation using streamline differential equation

Fiber orientation distribution estimation with diffusion magnetic resonance imaging (dMRI) is critical in white matter fiber tractography which is most commonly based on symmetric crossing structures. Ambiguous spatial correspondence between estimated diffusion directions and fiber geometry, such as asymmetric crossing, bending, fanning, or kissing, makes tractography challenging. Consequently, numerous tracts suggest intertwined connections in unexpected regions of the white matter or actually stop prematurely in the white matter. In this work, we propose a novel asymmetric fiber trajectory distribution (FTD) function defined on neighboring voxels based on a streamline differential equation from fluid mechanics. The spatial consistency with intra- and inter-voxel constraints is derived for FTD estimation by introducing the concept of divergence. At a local level, the FTD is a series of curve flows that minimize the energy function, characterizes the relations between fibers and joint fiber fragments within the same fiber bundle. Experiments are performed on FiberCup phantom, ISMRM 2015 Tractography challenge data, and in vivo brain dMRI data for qualitative and quantitative evaluations. Results show that our approach can reveal continuous asymmetric FTD details that are potentially useful for robust tractography.

Mapping Short Association Fibers in the Early Cortical Visual Processing Stream Using In Vivo Diffusion Tractography

Short association fibers (U-fibers) connect proximal cortical areas and constitute the majority of white matter connections in the human brain. U-fibers play an important role in brain development, function, and pathology but are underrepresented in current descriptions of the human brain connectome, primarily due to methodological challenges in diffusion magnetic resonance imaging (dMRI) of these fibers. High spatial resolution and dedicated fiber and tractography models are required to reliably map the U-fibers. Moreover, limited quantitative knowledge of their geometry and distribution makes validation of U-fiber tractography challenging. Submillimeter resolution diffusion MRI—facilitated by a cutting-edge MRI scanner with 300 mT/m maximum gradient amplitude—was used to map U-fiber connectivity between primary and secondary visual cortical areas (V1 and V2, respectively) in vivo. V1 and V2 retinotopic maps were obtained using functional MRI at 7T. The mapped V1–V2 connectivity was retinotopically organized, demonstrating higher connectivity for retinotopically corresponding areas in V1 and V2 as expected. The results were highly reproducible, as demonstrated by repeated measurements in the same participants and by an independent replication group study. This study demonstrates a robust U-fiber connectivity mapping in vivo and is an important step toward construction of a more complete human brain connectome.

Effect of Fixed-Density Thresholding on Structural Brain Networks: A Demonstration in Cerebral Small Vessel Disease

A popular solution to control for edge density variability in structural brain network analysis is to threshold the networks to a fixed density across all subjects. However, it remains unclear how this type of thresholding affects the basic network architecture in terms of edge weights, hub location, and hub connectivity and, especially, how it affects the sensitivity to detect disease-related abnormalities. We investigated these two questions in a cohort of patients with cerebral small vessel disease and age-matched controls. Brain networks were reconstructed from diffusion magnetic resonance imaging data using deterministic fiber tractography. Networks were thresholded to a fixed density by removing edges with the lowest number of streamlines. We compared edge length (mm), fractional anisotropy (FA), proportion of hub connections, and hub location between the unthresholded and the thresholded networks of each subject. Moreover, we compared weighted graph measures of global and local connectivity obtained from the (un)thresholded networks between patients and controls. We performed these analyses over a range of densities (2-20%). Results indicate that fixed-density thresholding disproportionally removes edges composed of long streamlines, but is independent of FA. The edges removed were not preferentially connected to hub or nonhub nodes. Over half of the original hubs were reproducible when networks were thresholded to a density ≥10%. Furthermore, the between-group differences in graph measures observed in the unthresholded network remained present after thresholding, irrespective of the chosen density. We therefore conclude that moderate fixed-density thresholds can successfully be applied to control for the effects of density in structural brain network analysis.

The role of diffusion tractography in refining glial tumor resection

Primary brain tumors are notoriously hard to resect surgically. Due to their infiltrative nature, finding the optimal resection boundary without damaging healthy tissue can be challenging. One potential tool to help make this decision is diffusion-weighted magnetic resonance imaging (dMRI) tractography. dMRI exploits the diffusion of water molecule along axons to generate a 3D modelization of the white matter bundles in the brain. This feature is particularly useful to visualize how a tumor affects its surrounding white matter and plan a surgical path. This paper reviews the different ways in which dMRI can be used to improve brain tumor resection, its benefits and also its limitations. We expose surgical tools that can be paired with dMRI to improve its impact on surgical outcome, such as loading the 3D tractography in the neuronavigation system and direct electrical stimulation to validate the position of the white matter bundles of interest. We also review articles validating dMRI findings using other anatomical investigation techniques, such as postmortem dissections, manganese-enhanced MRI, electrophysiological stimulations, and phantom studies with known ground truth. We will be discussing the areas of the brain where dMRI performs well and where the future challenges are. We will conclude this review with suggestions and take home messages for neurosurgeons, tractographers, and vendors for advancing the field and on how to benefit from tractography's use in clinical practice.

Diffusion tensor imaging of the corpus callosum in healthy aging: Investigating higher order polynomial regression modelling

Previous diffusion tensor imaging (DTI) studies confirmed the vulnerability of corpus callosum (CC) fibers to aging. However, most studies employed lower order regressions to study the relationship between age and white matter microstructure. The present study investigated whether higher order polynomial regression modelling can better describe the relationship between age and CC DTI metrics compared to lower order models in 140 healthy participants (ages 18-85). The CC was found to be non-uniformly affected by aging, with accelerated and earlier degradation occurring in anterior portion; callosal volume, fiber count, fiber length, mean fibers per voxel, and FA decreased with age while mean, axial, and radial diffusivities increased. Half of the parameters studied also displayed significant age-sex interaction or intracranial volume effects. Higher order models were chosen as the best fit, based on Bayesian Information Criterion minimization, in 16 out of 23 significant cases when describing the relationship between DTI measurements and age. Higher order model fits provided different estimations of aging trajectory peaks and decline onsets than lower order models; however, a likelihood ratio test found that higher order regressions generally did not fit the data significantly better than lower order polynomial or linear models. The results contrast the modelling approaches and highlight the importance of using higher order polynomial regression modelling when investigating associations between age and CC white matter microstructure.

Mitigating gyral bias in cortical tractography via asymmetric fiber orientation distributions

Diffusion tractography in brain connectomics often involves tracing axonal trajectories across gray-white matter boundaries in gyral blades of complex cortical convolutions. To date, gyral bias is observed in most tractography algorithms with streamlines predominantly terminating at gyral crowns instead of sulcal banks. This work demonstrates that asymmetric fiber orientation distribution functions (AFODFs), computed via a multi-tissue global estimation framework, can mitigate the effects of gyral bias, enabling fiber streamlines at gyral blades to make sharper turns into the cortical gray matter. We use ex-vivo data of an adult rhesus macaque and in-vivo data from the Human Connectome Project (HCP) to show that the fiber streamlines given by AFODFs bend more naturally into the cortex than the conventional symmetric FODFs in typical gyral blades. We demonstrate that AFODF tractography improves cortico-cortical connectivity and provides highly consistent outcomes between two different field strengths (3T and 7T).

Tractography reproducibility challenge with empirical data (TraCED): The 2017 ISMRM diffusion study group challenge

BACKGROUND

Fiber tracking with diffusion-weighted MRI has become an essential tool for estimating in vivo brain white matter architecture. Fiber tracking results are sensitive to the choice of processing method and tracking criteria.

PURPOSE

To assess the variability for an algorithm in group studies reproducibility is of critical context. However, reproducibility does not assess the validity of the brain connections. Phantom studies provide concrete quantitative comparisons of methods relative to absolute ground truths, yet do no capture variabilities because of in vivo physiological factors. The ISMRM 2017 TraCED challenge was created to fulfill the gap.

STUDY TYPE

A systematic review of algorithms and tract reproducibility studies.

SUBJECTS

Single healthy volunteers.

FIELD STRENGTH/SEQUENCE

3.0T, two different scanners by the same manufacturer. The multishell acquisition included b-values of 1000, 2000, and 3000 s/mm2 with 20, 45, and 64 diffusion gradient directions per shell, respectively.

ASSESSMENT

Nine international groups submitted 46 tractography algorithm entries each consisting 16 tracts per scan. The algorithms were assessed using intraclass correlation (ICC) and the Dice similarity measure.

STATISTICAL TESTS

Containment analysis was performed to assess if the submitted algorithms had containment within tracts of larger volume submissions. This also serves the purpose to detect if spurious submissions had been made.

RESULTS

The top five submissions had high ICC and Dice >0.88. Reproducibility was high within the top five submissions when assessed across sessions or across scanners: 0.87-0.97. Containment analysis shows that the top five submissions are contained within larger volume submissions. From the total of 16 tracts as an outcome relatively the number of tracts with high, moderate, and low reproducibility were 8, 4, and 4.

DATA CONCLUSION

The different methods clearly result in fundamentally different tract structures at the more conservative specificity choices. Data and challenge infrastructure remain available for continued analysis and provide a platform for comparison.

LEVEL OF EVIDENCE

5 Technical Efficacy Stage: 1 J. Magn. Reson. Imaging 2020;51:234-249.

Regions of white matter abnormalities in the arcuate fasciculus in veterans with anger and aggression problems

Aggression after military deployment is a common occurrence in veterans. Neurobiological research has shown that aggression is associated with a dysfunction in a network connecting brain regions implicated in threat processing and emotion regulation. However, aggression may also be related to deficits in networks underlying communication and social cognition. The uncinate and arcuate fasciculi are integral to these networks, thus studying potential abnormalities in these white matter connections can further our understanding of anger and aggression problems in military veterans. Here, we use diffusion tensor imaging tractography to investigate white matter microstructural properties of the uncinate fasciculus and the arcuate fasciculus in veterans with and without anger and aggression problems. A control tract, the parahippocampal cingulum was also included in the analyses. More specifically, fractional anisotropy (FA) estimates are derived along the trajectory from all fiber pathways and compared between both groups. No between-group FA differences are observed for the uncinate fasciculus and the cingulum, however parts of the arcuate fasciculus show a significantly lower FA in the group of veterans with aggression and anger problems. Our data suggest that abnormalities in arcuate fasciculus white matter connectivity that are related to self-regulation may play an important role in the etiology of anger and aggression in military veterans.

Spatial versus angular resolution for tractography-assisted planning of deep brain stimulation

•

Deep brain stimulation (DBS) benefits from precise targeting of white matter tracts.

•

Better to increase spatial vs. angular resolution for separating parallel tracts.

•

Scanning time trade-off between angular & spatial resolution depends on local anatomy.

•

We recommend increased spatial resolution dMRI for tract-guided internal capsule DBS.

Given the restricted total scanning time for clinical neuroimaging, it is unclear whether clinical diffusion MRI protocols would benefit more from higher spatial resolution or higher angular resolution. In this work, we investigated the relative benefit of improving spatial or angular resolution in diffusion MRI to separate two parallel running white matter tracts that are targets for deep brain stimulation: the anterior thalamic radiation and the supero-lateral branch of the medial forebrain bundle. Both these tracts are situated in the ventral anterior limb of the internal capsule, and recent studies suggest that targeting a specific tract could improve treatment efficacy. Therefore, we scanned 19 healthy volunteers at 3T and 7T according to three diffusion MRI protocols with respectively standard clinical settings, increased spatial resolution of 1.4 mm, and increased angular resolution (64 additional gradient directions at b = 2200s/mm²). We performed probabilistic tractography for all protocols and quantified the separability of both tracts. The higher spatial resolution protocol improved separability by 41% with respect to the clinical standard, presumably due to decreased partial voluming. The higher angular resolution protocol resulted in increased apparent tract volumes and overlap, which is disadvantageous for application in precise treatment planning. We thus recommend to increase the spatial resolution for deep brain stimulation planning to 1.4 mm while maintaining angular resolution. This recommendation complements the general advice to aim for high angular resolution to resolve crossing fibers, confirming that the specific application and anatomical considerations are leading in clinical diffusion MRI protocol optimization.

Differentiation of multiple system atrophy from Parkinson's disease by structural connectivity derived from probabilistic tractography

Recent studies combining diffusion tensor-derived metrics and machine learning have shown promising results in the discrimination of multiple system atrophy (MSA) and Parkinson’s disease (PD) patients. This approach has not been tested using more complex methodologies such as probabilistic tractography. The aim of this work is assessing whether the strength of structural connectivity between subcortical structures, measured as the number of streamlines (NOS) derived from tractography, can be used to classify MSA and PD patients at the single-patient level. The classification performance of subcortical FA and MD was also evaluated to compare the discriminant ability between diffusion tensor-derived metrics and NOS. Using diffusion-weighted images acquired in a 3 T MRI scanner and probabilistic tractography, we reconstructed the white matter tracts between 18 subcortical structures from a sample of 54 healthy controls, 31 MSA patients and 65 PD patients. NOS between subcortical structures were compared between groups and entered as features into a machine learning algorithm. Reduced NOS in MSA compared with controls and PD were found in connections between the putamen, pallidum, ventral diencephalon, thalamus, and cerebellum, in both right and left hemispheres. The classification procedure achieved an overall accuracy of 78%, with 71% of the MSA subjects and 86% of the PD patients correctly classified. NOS features outperformed the discrimination performance obtained with FA and MD. Our findings suggest that structural connectivity derived from tractography has the potential to correctly distinguish between MSA and PD patients. Furthermore, NOS measures obtained from tractography might be more useful than diffusion tensor-derived metrics for the detection of MSA.

Negative Affect Shared with Siblings is Associated with Structural Brain Network Efficiency and Loneliness in Adolescents

Loneliness has a strong neurobiological basis reflected by its specific relationships with structural brain connectivity. Critically, affect traits are highly related to loneliness, which shows close association with the onset and severity of major depressive disorder. This diffusion imaging study was conducted on a sample of adolescent siblings to examine whether positive and negative affect traits were related to loneliness, with brain network efficiency playing a mediating role. The findings of this study confirmed that both global and average local efficiency negatively mediated the association between low positive affect and high negative affect and loneliness, and the mediation was more sensitive to sibling-shared affect traits. The findings have important implications for interventions targeted at reducing the detrimental impact of familiar negative emotional experiences and loneliness.

Age-Related Differences in Functional and Structural Connectivity in the Spatial Navigation Brain Network

Spatial navigation involves multiple cognitive processes including multisensory integration, visuospatial coding, memory, and decision-making. These functions are mediated by the interplay of cerebral structures that can be broadly separated into a posterior network (subserving visual and spatial processing) and an anterior network (dedicated to memory and navigation planning). Within these networks, areas such as the hippocampus (HC) are known to be affected by aging and to be associated with cognitive decline and navigation impairments. However, age-related changes in brain connectivity within the spatial navigation network remain to be investigated. For this purpose, we performed a neuroimaging study combining functional and structural connectivity analyses between cerebral regions involved in spatial navigation. Nineteen young (μ = 27 years, σ = 4.3; 10 F) and 22 older (μ = 73 years, σ = 4.1; 10 F) participants were examined in this study. Our analyses focused on the parahippocampal place area (PPA), the retrosplenial cortex (RSC), the occipital place area (OPA), and the projections into the visual cortex of central and peripheral visual fields, delineated from independent functional localizers. In addition, we segmented the HC and the medial prefrontal cortex (mPFC) from anatomical images. Our results show an age-related decrease in functional connectivity between low-visual areas and the HC, associated with an increase in functional connectivity between OPA and PPA in older participants compared to young subjects. Concerning the structural connectivity, we found age-related differences in white matter integrity within the navigation brain network, with the exception of the OPA. The OPA is known to be involved in egocentric navigation, as opposed to allocentric strategies which are more related to the hippocampal region. The increase in functional connectivity between the OPA and PPA may thus reflect a compensatory mechanism for the age-related alterations around the HC, favoring the use of the preserved structural network mediating egocentric navigation. Overall, these findings on age-related differences of functional and structural connectivity may help to elucidate the cerebral bases of spatial navigation deficits in healthy and pathological aging.

Brain structural connectivity network alterations in insomnia disorder reveal a central role of the right angular gyrus

•

People with insomnia show widespread brain structural hyperconnectivity.

•

The right angular gyrus is central to the structural connectivity alterations.

•

Connectivity of this angular gyrus subnetwork correlates with reactive hyperarousal.

•

Brain structural hyperconnectivity may mark vulnerability to insomnia.

Insomnia Disorder (ID) is a prevalent and persistent condition, yet its neural substrate is not well understood. The cognitive, emotional, and behavioral characteristics of ID suggest that vulnerability involves distributed brain networks rather than a single brain area or connection. The present study utilized probabilistic diffusion tractography to compare the whole-brain structural connectivity networks of people with ID and those of matched controls without sleep complaints.

Diffusion-weighted images and T1-weighed images were acquired in 51 people diagnosed with ID (21–69 years of age, 37 female) and 48 matched controls without sleep complaints (22–70 years of age, 31 female). Probabilistic tractography was performed to construct the whole-brain structural connectivity network of each participant. Case–control differences in connectivity strength and network efficiency were evaluated with permutation tests.

People with ID showed structural hyperconnectivity within a subnetwork that spread over frontal, parietal, temporal, and subcortical regions and was anchored at the right angular gyrus. The result was robust across different edge-weighting strategies. Moreover, converging support was given by the finding of heightened right angular gyrus nodal efficiency (harmonic centrality) across varying graph density in people with ID. Follow-up correlation analyses revealed that subnetwork connectivity was associated with self-reported reactive hyperarousal.

The findings demonstrate that the right angular gyrus is a hub of enhanced structural connectivity in ID. Hyperconnectivity within the identified subnetwork may contribute to increased reactivity to stimuli and may signify vulnerability to ID.

The Seven Deadly Sins of Measuring Brain Structural Connectivity Using Diffusion MRI Streamlines Fibre-Tracking

There is great interest in the study of brain structural connectivity, as white matter abnormalities have been implicated in many disease states. Diffusion magnetic resonance imaging (MRI) provides a powerful means to characterise structural connectivity non-invasively, by using a fibre-tracking algorithm. The most widely used fibre-tracking strategy is based on the step-wise generation of streamlines. Despite their popularity and widespread use, there are a number of practical considerations that must be taken into account in order to increase the robustness of streamlines tracking results, particularly when these methods are used to study brain structural connectivity, and the connectome. This review article describes what we consider the ‘seven deadly sins’ of mapping structural connections using diffusion MRI streamlines fibre-tracking, with particular emphasis on ‘sins’ that can be practically avoided and they can have an important impact in the results. It is shown that there are important ‘deadly sins’ to be avoided at every step of the pipeline, such as during data acquisition, during data modelling to estimate local fibre architecture, during the fibre-tracking process itself, and during quantification of the tracking results. The recommendations here are intended to inform users on potential important shortcomings of their current tracking protocols, as well as to guide future users on some of the key issues and decisions that must be faced when designing their processing pipelines.

Early Detection of Cerebral Palsy Using Sensorimotor Tract Biomarkers in Very Preterm Infants

BACKGROUND

Our objectives were to evaluate the brain's sensorimotor network microstructure using diffusion magnetic resonance imaging (MRI) at term-corrected age and test the ability of sensorimotor microstructural parameters to accurately predict cerebral palsy in extremely-low-birth-weight infants.

METHODS

We enrolled a prospective pilot cohort of extremely-low-birth-weight preterm infants (birth weight ≤ 1000 g) before neonatal intensive care unit discharge and studied them with structural and diffusion MRI at term-corrected age. Six sensorimotor tracts were segmented, and microstructural parameters from these tracts were evaluated for their ability to predict later development of cerebral palsy, diagnosed at 18 to 22 months corrected age.

RESULTS

We found significant differences in multiple diffusion MRI parameters from five of the six sensorimotor tracts in infants who developed cerebral palsy (n = 5) versus those who did not (n = 36). When compared with structural MRI or individual diffusion MRI biomarkers, the combination of two individual biomarkers-fractional anisotropy of superior thalamic radiations (sensory component) and radial diffusivity of the corticospinal tract-exhibited the highest sensitivity (80%), specificity (97%), and positive likelihood ratio (28.0) for prediction of cerebral palsy. This combination of diffusion MRI biomarkers accurately classified 95% of the study infants.

CONCLUSIONS

Development of cerebral palsy in very preterm infants is preceded by early brain injury or immaturity to one or more sensorimotor tracts. A larger study is warranted to evaluate if a combination of sensorimotor microstructural biomarkers could accurately facilitate early diagnosis of cerebral palsy.

In vivo magnetic resonance imaging and spectroscopy. Technological advances and opportunities for applications continue to abound

Over the past decades, the field of in vivo magnetic resonance (MR) has built up an impressive repertoire of data acquisition and analysis technologies for anatomical, functional, physiological, and molecular imaging, the description of which requires many book volumes. As such it is impossible for a few authors to have an authoritative overview of the field and for a brief article to be inclusive. We will therefore focus mainly on data acquisition and attempt to give some insight into the principles underlying current advanced methods in the field and the potential for further innovation. In our view, the foreseeable future is expected to show continued rapid progress, for instance in imaging of microscopic tissue properties in vivo, assessment of functional and anatomical connectivity, higher resolution physiologic and metabolic imaging, and even imaging of receptor binding. In addition, acquisition speed and information content will continue to increase due to the continuous development of approaches for parallel imaging (including simultaneous multi-slice imaging), compressed sensing, and MRI fingerprinting. Finally, artificial intelligence approaches are becoming more realistic and will have a tremendous effect on both acquisition and analysis strategies. Together, these developments will continue to provide opportunity for scientific discovery and, in combination with large data sets from other fields such as genomics, allow the ultimate realization of precision medicine in the clinic.

A Web-Based Atlas Combining MRI and Histology of the Squirrel Monkey Brain

The squirrel monkey (Saimiri sciureus) is a commonly-used surrogate for humans in biomedical research. In the neuroimaging community, MRI and histological atlases serve as valuable resources for anatomical, physiological, and functional studies of the brain; however, no digital MRI/histology atlas is currently available for the squirrel monkey. This paper describes the construction of a web-based multi-modal atlas of the squirrel monkey brain. The MRI-derived information includes anatomical MRI contrast (i.e., T2-weighted and proton-density-weighted) and diffusion MRI metrics (i.e., fractional anisotropy and mean diffusivity) from data acquired both in vivo and ex vivo on a 9.4 Tesla scanner. The histological images include Nissl and myelin stains, co-registered to the corresponding MRI, allowing identification of cyto- and myelo-architecture. In addition, a bidirectional neuronal tracer, biotinylated dextran amine (BDA) was injected into the primary motor cortex, enabling highly specific identification of regions connected to the injection location. The atlas integrates the results of common image analysis methods including diffusion tensor imaging glyphs, labels of 57 white-matter tracts identified using DTI-tractography, and 18 cortical regions of interest identified from Nissl-revealed cyto-architecture. All data are presented in a common space, and all image types are accessible through a web-based atlas viewer, which allows visualization and interaction of user-selectable contrasts and varying resolutions. By providing an easy to use reference system of anatomical information, our web-accessible multi-contrast atlas forms a rich and convenient resource for comparisons of brain findings across subjects or modalities. The atlas is called the Combined Histology-MRI Integrated Atlas of the Squirrel Monkey (CHIASM). All images are accessible through our web-based viewer ( https://chiasm.vuse.vanderbilt.edu /), and data are available for download at ( https://www.nitrc.org/projects/smatlas/ ).

Cognitive Training in Young Patients With Traumatic Brain Injury: A Fixel-Based Analysis

Background. Traumatic brain injury (TBI) is associated with altered white matter organization and impaired cognitive functioning. Objective. We aimed to investigate changes in white matter and cognitive functioning following computerized cognitive training. Methods. Sixteen adolescents with moderate-to-severe TBI (age 15.6 ± 1.8 years, 1.2-4.6 years postinjury) completed the 8-week BrainGames program and diffusion weighted imaging (DWI) and cognitive assessment at time point 1 (before training) and time point 2 (after training). Sixteen healthy controls (HC) (age 15.6 ± 1.8 years) completed DWI assessment at time point 1 and cognitive assessment at time point 1 and 2. Fixel-based analyses were used to examine fractional anisotropy (FA), mean diffusivity (MD), and fiber cross-section (FC) on a whole brain level and in tracts of interest. Results. Patients with TBI showed cognitive impairments and extensive areas with decreased FA and increased MD together with an increase in FC in the body of the corpus callosum and left superior longitudinal fasciculus (SLF) at time point 1. Patients improved significantly on the inhibition measure at time point 2, whereas the HC group remained unchanged. No training-induced changes were observed on the group level in diffusion metrics. Exploratory correlations were found between improvements on verbal working memory and reduced MD of the left SLF and between increased performance on an information processing speed task and increased FA of the right precentral gyrus. Conclusions. Results are indicative of positive effects of BrainGames on cognitive functioning and provide preliminary evidence for neuroplasticity associated with cognitive improvements following cognitive intervention in TBI.

Optimization of data acquisition and analysis for fiber ball imaging

The inverse Funk transform of high angular resolution diffusion imaging (HARDI) data provides an estimate for the fiber orientation density function (fODF) in white matter (WM). Since the inverse Funk transform is a straightforward linear transformation, this technique, referred to as fiber ball imaging (FBI), offers a practical means of calculating the fODF that avoids the need for a response function or nonlinear numerical fitting. Nevertheless, the accuracy of FBI depends on both the choice of b-value and the number of diffusion-encoding directions used to acquire the HARDI data. To inform the design of optimal scan protocols for its implementation, FBI predictions are investigated here with in vivo data from healthy adult volunteers acquired at 3 T for b-values spanning 1000 to 10,000 s/mm², for diffusion-encoding directions varying in number from 30 to 256 and for TE ranging from 90 to 120 ms. Our results suggest b-values above 4000 s/mm² with at least 64 diffusion-encoding directions are adequate to achieve reasonable accuracy with FBI for calculating axon-specific diffusion measures and for performing WM fiber tractography (WMFT).

Disrupted structural connectivity of fronto-deep gray matter pathways in progressive supranuclear palsy

Background

Structural connectivity is a promising methodology to detect patterns of neural network dysfunction in neurodegenerative diseases. This approach has not been tested in progressive supranuclear palsy (PSP).

Objectives

The aim of this study is reconstructing the structural connectome to characterize and detect the pathways of degeneration in PSP patients compared with healthy controls and their correlation with clinical features. The second objective is to assess the potential of structural connectivity measures to distinguish between PSP patients and healthy controls at the single-subject level.

Methods

Twenty healthy controls and 19 PSP patients underwent diffusion-weighted MRI with a 3T scanner. Structural connectivity, represented by number of streamlines, was derived from probabilistic tractography. Global and local network metrics were calculated based on graph theory.

Results

Reduced numbers of streamlines were predominantly found in connections between frontal areas and deep gray matter (DGM) structures in PSP compared with controls. Significant changes in structural connectivity correlated with clinical features in PSP patients. An abnormal small-world architecture was detected in the subnetwork comprising the frontal lobe and DGM structures in PSP patients. The classification procedure achieved an overall accuracy of 82.23% with 94.74% sensitivity and 70% specificity.

Conclusion

Our findings suggest that modelling the brain as a structural connectome is a useful method to detect changes in the organization and topology of white matter tracts in PSP patients. Secondly, measures of structural connectivity have the potential to correctly discriminate between PSP patients and healthy controls.

•

Reduced structural connectivity in PSP patients compared with healthy controls

•

Connectivity reductions in fronto-DGM tracts correlate with PSPRS and FAB scores

•

PSP patients present abnormal small-world architecture in the fronto-DGM network.

Reducing variability in along-tract analysis with diffusion profile realignment

Diffusion weighted magnetic resonance imaging (dMRI) provides a non invasive virtual reconstruction of the brain's white matter structures through tractography. Analyzing dMRI measures along the trajectory of white matter bundles can provide a more specific investigation than considering a region of interest or tract-averaged measurements. However, performing group analyses with this along-tract strategy requires correspondence between points of tract pathways across subjects. This is usually achieved by creating a new common space where the representative streamlines from every subject are resampled to the same number of points. If the underlying anatomy of some subjects was altered due to, e.g., disease or developmental changes, such information might be lost by resampling to a fixed number of points. In this work, we propose to address the issue of possible misalignment, which might be present even after resampling, by realigning the representative streamline of each subject in this 1D space with a new method, coined diffusion profile realignment (DPR). Experiments on synthetic datasets show that DPR reduces the coefficient of variation for the mean diffusivity, fractional anisotropy and apparent fiber density when compared to the unaligned case. Using 100 in vivo datasets from the human connectome project, we simulated changes in mean diffusivity, fractional anisotropy and apparent fiber density. Independent Student's t-tests between these altered subjects and the original subjects indicate that regional changes are identified after realignment with the DPR algorithm, while preserving differences previously detected in the unaligned case. This new correction strategy contributes to revealing effects of interest which might be hidden by misalignment and has the potential to improve the specificity in longitudinal population studies beyond the traditional region of interest based analysis and along-tract analysis workflows.