In vivo microstructural heterogeneity of white matter lesions in healthy elderly and Alzheimer's disease participants using tissue compositional analysis of diffusion MRI data

White matter hyperintensities (WMH) are regions of high signal intensity typically identified on fluid attenuated inversion recovery (FLAIR). Although commonly observed in elderly individuals, they are more prevalent in Alzheimer's disease (AD) patients. Given that WMH appear relatively homogeneous on FLAIR, they are commonly partitioned into location- or distance-based classes when investigating their relevance to disease. Since pathology indicates that such lesions are often heterogeneous, probing their microstructure in vivo may provide greater insight than relying on such arbitrary classification schemes. In this study, we investigated WMH in vivo using an advanced diffusion MRI method known as single-shell 3-tissue constrained spherical deconvolution (SS3T-CSD), which models white matter microstructure while accounting for grey matter and CSF compartments. Diffusion MRI data and FLAIR images were obtained from AD (n = 48) and healthy elderly control (n = 94) subjects. WMH were automatically segmented, and classified: (1) as either periventricular or deep; or (2) into three distance-based contours from the ventricles. The 3-tissue profile of WMH enabled their characterisation in terms of white matter-, grey matter-, and fluid-like characteristics of the diffusion signal. Our SS3T-CSD findings revealed substantial heterogeneity in the 3-tissue profile of WMH, both within lesions and across the various classes. Moreover, this heterogeneity information indicated that the use of different commonly used WMH classification schemes can result in different disease-based conclusions. We conclude that future studies of WMH in AD would benefit from inclusion of microstructural information when characterising lesions, which we demonstrate can be performed in vivo using SS3T-CSD.

Fibre-specific white matter changes in multiple sclerosis patients with optic neuritis

Long term irreversible disability in multiple sclerosis (MS) is thought to be primarily driven by axonal degeneration. Axonal degeneration leads to degenerative atrophy, therefore early markers of axonal degeneration are required to predict clinical disability and treatment efficacy. Given that additional pathologies such as inflammation, demyelination and oedema are also present in MS, it is essential to develop axonal markers that are not confounded by these processes. The present study investigated a novel method for measuring axonal degeneration in MS based on high angular resolution diffusion magnetic resonance imaging. Unlike standard methods, this novel method involved advanced acquisition and modelling for improved axonal sensitivity and specificity. Recent work has developed analytical methods, two novel axonal markers, fibre density and cross-section, that can be estimated for each fibre direction in each voxel (termed a “fixel”). This technique, termed fixel-based analysis, thus simultaneously estimates axonal density and white matter atrophy from specific white matter tracts. Diffusion-weighted imaging datasets were acquired for 17 patients with a history of acute unilateral optic neuritis (35.3 ± 10.2 years, 11 females) and 14 healthy controls (32.7 ± 4.8 years, 8 females) on a 3 T scanner. Fibre density values were compared to standard diffusion tensor imaging parameters (fractional anisotropy and mean diffusivity) in lesions and normal appearing white matter. Group comparisons were performed for each fixel to assess putative differences in fibre density and fibre cross-section. Fibre density was observed to have a comparable sensitivity to fractional anisotropy for detecting white matter pathology in MS, but was not affected by crossing axonal fibres. Whole brain fixel-based analysis revealed significant reductions in fibre density and fibre cross-section in the inferior fronto-occipital fasciculus (including the optic radiations) of patients compared to controls. We interpret this result to indicate that this fixel-based approach is able to detect early loss of fibre density and cross-section in the optic radiations in MS patients with a history of optic neuritis. Fibre-specific markers of axonal degeneration should be investigated further for use in early stage therapeutic trials, or to monitor axonal injury in early stage MS.

•

Fibre density is reduced in lesions and normal-appearing white matter in MS

•

Fibre density detects white matter pathology in regions of crossing fibres

•

Loss of fibre density and cross-section selectively evident in visual pathways of optic neuritis patients.

BIDS apps: Improving ease of use, accessibility, and reproducibility of neuroimaging data analysis methods

The rate of progress in human neurosciences is limited by the inability to easily apply a wide range of analysis methods to the plethora of different datasets acquired in labs around the world. In this work, we introduce a framework for creating, testing, versioning and archiving portable applications for analyzing neuroimaging data organized and described in compliance with the Brain Imaging Data Structure (BIDS). The portability of these applications (BIDS Apps) is achieved by using container technologies that encapsulate all binary and other dependencies in one convenient package. BIDS Apps run on all three major operating systems with no need for complex setup and configuration and thanks to the comprehensiveness of the BIDS standard they require little manual user input. Previous containerized data processing solutions were limited to single user environments and not compatible with most multi-tenant High Performance Computing systems. BIDS Apps overcome this limitation by taking advantage of the Singularity container technology. As a proof of concept, this work is accompanied by 22 ready to use BIDS Apps, packaging a diverse set of commonly used neuroimaging algorithms.

Predicting hand function after hemidisconnection

Hemidisconnections (i.e. hemispherectomies or hemispherotomies) invariably lead to contralateral hemiparesis. Many patients with a pre-existing hemiparesis, however, experience no deterioration in motor functions, and some can still grasp with their paretic hand after hemidisconnection. The scope of our study was to predict this phenomenon. Hypothesizing that preserved contralateral grasping ability after hemidisconnection can only occur in patients controlling their paretic hands via ipsilateral corticospinal projections already in the preoperative situation, we analysed the asymmetries of the brainstem (by manual magnetic resonance imaging volumetry) and of the structural connectivity of the corticospinal tracts within the brainstem (by magnetic resonance imaging diffusion tractography), assuming that marked hypoplasia or Wallerian degeneration on the lesioned side in patients who can grasp with their paretic hands indicate ipsilateral control. One hundred and two patients who underwent hemidisconnections between 0.8 and 36 years of age were included. Before the operation, contralateral hand function was normal in 3/102 patients, 47/102 patients showed hemiparetic grasping ability and 52/102 patients could not grasp with their paretic hands. After hemidisconnection, 20/102 patients showed a preserved grasping ability, and 5/102 patients began to grasp with their paretic hands only after the operation. All these 25 patients suffered from pre- or perinatal brain lesions. Thirty of 102 patients lost their grasping ability. This group included all seven patients with a post-neonatally acquired or progressive brain lesion who could grasp before the operation, and also all three patients with a preoperatively normal hand function. The remaining 52/102 patients were unable to grasp pre- and postoperatively. On magnetic resonance imaging, the patients with preserved grasping showed significantly more asymmetric brainstem volumes than the patients who lost their grasping ability. Similarly, these patients showed striking asymmetries in the structural connectivity of the corticospinal tracts. In summary, normal preoperative hand function and a post-neonatally acquired or progressive lesion predict a loss of grasping ability after hemidisconnection. A postoperatively preserved grasping ability is possible in hemiparetic patients with pre- or perinatal lesions, and this is highly likely when the brainstem is asymmetric and especially when the structural connectivity of the corticospinal tracts within the brainstem is asymmetric.

HOMOR: higher order model outlier rejection for high b-value MR diffusion data

Diffusion MR images are prone to artefacts caused by head movement and cardiac pulsation. Previous techniques for the automated voxel-wise detection of signal intensity outliers have relied on the fit of the diffusion tensor to the data (RESTORE). However, the diffusion tensor cannot appropriately model more than a single fibre population, which may lead to inaccuracies when identifying outlier voxels in crossing fibre regions, particularly when high b-values are used to obtain increased angular contrast. HOMOR (higher order model outlier rejection) was developed to overcome this limitation and is introduced in this study. HOMOR is closely related to RESTORE, but employs a higher order model capable of resolving multiple fibre populations within a voxel. Using high b-value (b=3000 s/mm2) diffusion data from a population of 90 healthy participants, as well as simulations, HOMOR was found to identify a decreased number of outlier voxels compared to RESTORE primarily within areas of crossing, bending and fanning fibres. At lower b-values, however, RESTORE and HOMOR give similar results, which is demonstrated using diffusion data acquired at b=1000 s/mm2 in a mixed cohort. This study demonstrates that, although RESTORE is suitable for low b-value data, HOMOR is better suited for high b-value data.

Reorientation of fiber orientation distributions using apodized point spread functions

Using high angular resolution diffusion-weighted images, spherical deconvolution enables multiple white matter fiber populations to be resolved within a single voxel by computing the fiber orientation distribution (FOD). Higher order information provided by FODs could improve several methods for investigating population differences in white matter, including image registration, voxel-based analysis, atlas-based segmentation and labeling, and group average fiber tractography. All of these methods require spatial normalization of FODs. In this article, a novel method to reorient the FOD is presented, which is an important step required for FOD spatial normalization. The proposed method was assessed using both qualitative and quantitative experiments, with numerical simulations and in vivo human data. Results demonstrate that the proposed method improves FOD reorientation accuracy, removes undesired artefacts, and decreases computation time compared to a previous approach. The utility of the proposed method is illustrated by nonlinear FOD spatial normalization of 10 human subjects. Accurate reorientation and normalization of FODs is a critical step toward investigating white matter tissue in the context of multiple fiber orientations.