Search

Remyelination is crucial to recover from inflammatory demyelination in multiple sclerosis (MS). Investigating remyelination in vivo using magnetic resonance imaging (MRI) is difficult in MS, where collecting serial short-interval scans is challenging. Using experimental autoimmune encephalomyelitis (EAE) in common marmosets, a model of MS that recapitulates focal cerebral inflammatory demyelinating lesions, we investigated whether MRI is sensitive to, and can characterize, remyelination. In six animals followed with multisequence 7 T MRI, 31 focal lesions, predicted to be demyelinated or remyelinated based on signal intensity on proton density-weighted images, were subsequently assessed with histopathology. Remyelination occurred in four of six marmosets and 45% of lesions. Radiological-pathological comparison showed that MRI had high statistical sensitivity (100%) and specificity (90%) for detecting remyelination. This study demonstrates the prevalence of spontaneous remyelination in marmoset EAE and the ability of in vivo MRI to detect it, with implications for preclinical testing of pro-remyelinating agents.

Regional analysis of volumes examined in normalized space (RAVENS) are transformation images used in the study of brain morphometry. In this paper, RAVENS images are analyzed using a longitudinal variant of voxel-based morphometry (VBM) and longitudinal functional principal component analysis (LFPCA) for high-dimensional images. We demonstrate that the latter overcomes the limitations of standard longitudinal VBM analyses, which does not separate registration errors from other longitudinal changes and baseline patterns. This is especially important in contexts where longitudinal changes are only a small fraction of the overall observed variability, which is typical in normal aging and many chronic diseases. Our simulation study shows that LFPCA effectively separates registration error from baseline and longitudinal signals of interest by decomposing RAVENS images measured at multiple visits into three components: a subject-specific imaging random intercept that quantifies the cross-sectional variability, a subject-specific imaging slope that quantifies the irreversible changes over multiple visits, and a subject-visit specific imaging deviation. We describe strategies to identify baseline/longitudinal variation and registration errors combined with covariates of interest. Our analysis suggests that specific regional brain atrophy and ventricular enlargement are associated with multiple sclerosis (MS) disease progression.

Brain atrophy is a well-accepted imaging biomarker of multiple sclerosis (MS) that partially correlates with both physical disability and cognitive impairment.Based on MRI scans of 60 MS cases and 37 healthy volunteers, we measured the volumes of white matter (WM) lesions, cortical gray matter (GM), cerebral WM, caudate nucleus, putamen, thalamus, ventricles, and brainstem using a validated and completely automated segmentation method. We correlated these volumes with the Expanded Disability Status Scale (EDSS), MS Severity Scale (MSSS), MS Functional Composite (MSFC), and quantitative measures of ankle strength and toe sensation. Normalized volumes of both cortical and subcortical GM structures were abnormally low in the MS group, whereas no abnormality was found in the volume of the cerebral WM. High physical disability was associated with low cerebral WM, thalamus, and brainstem volumes (partial correlation coefficients ~0.3-0.4) but not with low cortical GM volume. Thalamus volumes were inversely correlated with lesion load (r = -0.36, p

Quantitative Susceptibility Mapping (QSM) is an MRI tool with the potential to reveal pathological changes from magnetic susceptibility measurements. Before phase data can be used to recover susceptibility (Δχ), the QSM process begins with two steps: data acquisition and phase estimation. We assess the performance of these steps, when applied without user intervention, on several variations of a phantom imaging task. We used a rotating-tube phantom with five tubes ranging from Δχ=0.05 ppm to Δχ=0.336 ppm. MRI data was acquired at nine angles of rotation for four different pulse sequences. The images were processed by 10 phase estimation algorithms including Laplacian, region-growing, branch-cut, temporal unwrapping, and maximum-likelihood methods, resulting in approximately 90 different combinations of data acquisition and phase estimation methods. We analyzed errors between measured and expected phases using the probability mass function and Cumulative Distribution Function. Repeatable acquisition and estimation methods were identified based on the probability of relative phase errors. For single-echo GRE and segmented EPI sequences, a region-growing method was most reliable with Pr (relative error

Advances in brain imaging and computational methods have facilitated the creation of subject-specific computational brain models that aid researchers in investigating brain trauma using simulated impacts. The emergence of magnetic resonance elastography (MRE) as a non-invasive mechanical neuroimaging tool has enabled in vivo estimation of material properties at low-strain, harmonic loading. An open question in the field has been how this data can be integrated into computational models. The goals of this study were to use a novel MRI dataset acquired in human volunteers to generate models with subject-specific anatomy and material properties, and then to compare simulated brain deformations to subject-specific brain deformation data under non-injurious loading. Models of five subjects were simulated with linear viscoelastic (LVE) material properties estimated directly from MRE data. Model predictions were compared to experimental brain deformation acquired in the same subjects using tagged MRI. Outcomes from the models matched the spatial distribution and magnitude of the measured peak strain components as well as the 95th percentile in-plane peak strains within 0.005 mm/mm and maximum principal strain within 0.012 mm/mm. Sensitivity to material heterogeneity was also investigated. Simulated brain deformations from a model with homogenous brain properties and a model with brain properties discretized with up to ten regions were very similar (a mean absolute difference less than 0.0015 mm/mm in peak strains). Incorporating material properties directly from MRE into a biofidelic subject-specific model is an important step toward future investigations of higher-order model features and simulations under more severe loading conditions. Statement of Significance: The study presents a method to calibrate and evaluate subject-specific finite element brain models using a combination of advanced magnetic resonance imaging (MRI) data. The imaging data is acquired in human volunteers and includes anatomical MRI, magnetic resonance elastography (MRE), and tagged MRI to generate subject-specific geometry, calibrate subject-specific material properties, and evaluate simulation response using subject-specific brain deformation. This dataset of MRE and tagged MRI allows for a unique evaluation of whether material properties from MRE can be used to create biofidelic computational models of the human brain. The study develops a calibration procedure to readily calculate linear viscoelastic material parameters from MRE data and then provides a sensitivity study of the effect of mechanical heterogeneity of the brain on simulation response. The calibrated computational models are used to simulate each subject's tagged MRI experiment; the results show good agreement between the simulated and experimental strain fields. The presented study and results will be informative in guiding the calibration of subject-specific computational brain model from experimental MRE data. The processed MRI, MRE, and tagged MRI data are publicly available at https://www.nitrc.org/projects/bbir/.

The rapid deformation of brain tissue in response to head impact can lead to traumatic brain injury. In vivo measurements of brain deformation during non-injurious head impacts are necessary to understand the underlying mechanisms of traumatic brain injury and compare to computational models of brain biomechanics. Using tagged magnetic resonance imaging (MRI), we obtained measurements of three-dimensional strain tensors that resulted from a mild head impact after neck rotation or neck extension. Measurements of maximum principal strain (MPS) peaked shortly after impact, with maximal values of 0.019–0.053 that correlated strongly with peak angular velocity. Subject-specific patterns of MPS were spatially heterogeneous and consistent across subjects for the same motion, though regions of high deformation differed between motions. The largest MPS values were seen in the cortical gray matter and cerebral white matter for neck rotation and the brainstem and cerebellum for neck extension. Axonal fiber strain (Ef) was estimated by combining the strain tensor with diffusion tensor imaging data. As with MPS, patterns of Ef varied spatially within subjects, were similar across subjects within each motion, and showed group differences between motions. Values were highest and most strongly correlated with peak angular velocity in the corpus callosum for neck rotation and in the brainstem for neck extension. The different patterns of brain deformation between head motions highlight potential areas of greater risk of injury between motions at higher loading conditions. Additionally, these experimental measurements can be directly compared to predictions of generic or subject-specific computational models of traumatic brain injury. Statement of Significance: Traumatic brain injury can result from the rapid acceleration of the skull, leading to deformation of brain tissue and elongation of axonal fibers. Because treatment options and prognostic models for patients are lacking, a better understanding of injury mechanisms is needed. Here, we use tagged magnetic resonance imaging to measure deformation throughout the live, human brain during non-injurious head accelerations. We present the first in vivo measurements of axonal stretch and compare MPS and axonal stretch experienced during neck rotation and neck extension. These results are important to elucidate brain regions at risk for injury. Additionally, they can be directly used to evaluate computational models of brain injury, which are used to predict risk of concussion during head impacts and design protective equipment.

Progressive multifocal leukoencephalopathy (PML) is a rare opportunistic brain infection caused by the JC virus and associated with substantial morbidity and mortality. Accurate MRI assessment of PML lesion burden and brain parenchymal atrophy is of decisive value in monitoring the disease course and response to therapy. However, there are currently no validated automatic methods for quantification of PML lesion burden or associated parenchymal volume loss. Furthermore, manual brain or lesion delineations can be tedious, require the use of valuable time resources by radiologists or trained experts, and are often subjective. In this work, we introduce JCnet (named after the causative viral agent), an end-to-end, fully automated method for brain parenchymal and lesion segmentation in PML using consecutive 3D patch-based convolutional neural networks. The network architecture consists of multi-view feature pyramid networks with hierarchical residual learning blocks containing embedded batch normalization and nonlinear activation functions. The feature maps across the bottom-up and top-down pathways of the feature pyramids are merged, and an output probability membership generated through convolutional pathways, thus rendering the method fully convolutional. Our results show that this approach outperforms and improves longitudinal consistency compared to conventional, state-of-the-art methods of healthy brain and multiple sclerosis lesion segmentation, utilized here as comparators given the lack of available methods validated for use in PML. The ability to produce robust and accurate automated measures of brain atrophy and lesion segmentation in PML is not only valuable clinically but holds promise toward including standardized quantitative MRI measures in clinical trials of targeted therapies. Code is available at: https://github.com/omarallouz/JCnet.

The data presented in this article is related to the research article entitled “Longitudinal multiple sclerosis lesion segmentation: Resource and challenge” (Carass et al., 2017) [1]. In conjunction with the 2015 International Symposium on Biomedical Imaging, we organized a longitudinal multiple sclerosis (MS) lesion segmentation challenge providing training and test data to registered participants. The training data consists of five subjects with a mean of 4.4 (±0.55) time-points, and test data of fourteen subjects with a mean of 4.4 (±0.67) time-points. All 82 data sets had the white matter lesions associated with multiple sclerosis delineated by two human expert raters. The training data including multi-modal scans and manually delineated lesion masks is available for download. In addition, the testing data is also being made available in conjunction with a website for evaluating the automated analysis of the testing data.