Search

Background: The “central vein sign” (CVS), a linear hypointensity on T2*-weighted imaging corresponding to a central vein/venule, is associated with multiple sclerosis (MS) lesions. The effect of lesion-size exclusion criteria on MS diagnostic accuracy has not been extensively studied. Objective: Investigate the optimal lesion-size exclusion criteria for CVS use in MS diagnosis. Methods: Cross-sectional study of 163 MS and 51 non-MS, and radiological/histopathological correlation of 5 MS and 1 control autopsy cases. The effects of lesion-size exclusion on MS diagnosis using the CVS, and intralesional vein detection on histopathology were evaluated. Results: CVS+ lesions were larger compared to CVS− lesions, with effect modification by MS diagnosis (mean difference +7.7 mm3, p = 0.004). CVS percentage-based criteria with no lesion-size exclusion showed the highest diagnostic accuracy in differentiating MS cases. However, a simple count of three or more CVS+ lesions greater than 3.5 mm is highly accurate and can be rapidly implemented (sensitivity 93%; specificity 88%). On magnetic resonance imaging (MRI)-histopathological correlation, the CVS had high specificity for identifying intralesional veins (0/7 false positives). Conclusion: Lesion-size measures add important information when using CVS+ lesion counts for MS diagnosis. The CVS is a specific biomarker corresponding to intralesional veins on histopathology.

Computational models of the brain and its biomechanical response to skull accelerations are important tools for understanding and predicting traumatic brain injuries (TBIs). However, most models have been developed using experimental data collected on animal models and cadaveric specimens, both of which differ from the living human brain. Here we describe efforts to noninvasively measure the biomechanical response of the human brain with MRI—at non-injurious strain levels—and generate data that can be used to develop, calibrate, and evaluate computational brain biomechanics models. Specifically, this paper reports on a project supported by the National Institute of Neurological Disorders and Stroke to comprehensively image brain anatomy and geometry, mechanical properties, and brain deformations that arise from impulsive and harmonic skull loadings. The outcome of this work will be a publicly available dataset ( http://www.nitrc.org/projects/bbir ) that includes measurements on both males and females across an age range from adolescence to older adulthood. This article describes the rationale and approach for this study, the data available, and how these data may be used to develop new computational models and augment existing approaches; it will serve as a reference to researchers interested in using these data.

Susceptibility-weighted imaging (SWI) provides superior image contrast of cerebral microhemorrhages (CMBs). It is based on a three-dimensional (3D) gradient echo (GRE) sequence with a relatively long imaging time.To evaluate whether an accelerated 3D segmented echo planar imaging SWI is comparable to GRE SWI in detecting CMBs in traumatic brain injury (TBI).Prospective.Four healthy volunteers and 46 consecutive subjects (38.0 ± 14.4 years, 16 females; 12 mild, 13 moderate, and 7 severe TBI).A 3 T scanner/3D gradient echo and 3D segmented echo planar imaging (segEPI).Brain images were acquired using GRE and segEPI in a single session (imaging time = 9 minutes 47 seconds and 1 minute 30 seconds, respectively). The signal-to-noise ratio (SNR) calculated from healthy volunteer thalamus and centrum semiovale were compared. CMBs were counted by three raters blinded to diagnostic information.A t-test was used to assess SNR difference. Pearson correlation and Wilcoxon signed-rank test were performed using CMB counts. The intermethod agreement was evaluated using Bland-Altman method. Intermethod and interrater reliabilities of image-based diffuse axonal injury (DAI) diagnoses were evaluated using Cohen's kappa and percent agreement. P ≤ 0.05 was considered statistically significant.Thalamus SNRs were 16.9 ± 2.2 and 16.5 ± 3 for GRE and segEPI (P = 0.84), respectively. Centrum semiovale SNRs were 25.8 ± 4.6 and 21.1 ± 2.7 (P = 0.13). The correlation coefficient of CMBs was 0.93, and differences were not significant (P = 0.56-0.85). For DAI diagnoses, Cohen's kappa was 0.62-0.84 and percent agreement was 85%-94%.CMB counts on segEPI and GRE were highly correlated, and DAI diagnosis was made equally effectively. segEPI SWI can potentially replace GRE SWI in detecting TBI CMBs, especially when time constraints are critical.1 TECHNICAL EFFICACY: Stage 2.

The Sørensen-Dice index (SDI) is a widely used measure for evaluating medical image segmentation algorithms. It offers a standardized measure of segmentation accuracy which has proven useful. However, it offers diminishing insight when the number of objects is unknown, such as in white matter lesion segmentation of multiple sclerosis (MS) patients. We present a refinement for finer grained parsing of SDI results in situations where the number of objects is unknown. We explore these ideas with two case studies showing what can be learned from our two presented studies. Our first study explores an inter-rater comparison, showing that smaller lesions cannot be reliably identified. In our second case study, we demonstrate fusing multiple MS lesion segmentation algorithms based on the insights into the algorithms provided by our analysis to generate a segmentation that exhibits improved performance. This work demonstrates the wealth of information that can be learned from refined analysis of medical image segmentations.

Background and ObjectivesThe central vein sign (CVS), a central linear hypointensity within lesions on T2*-weighted imaging, has been established as a sensitive and specific biomarker for the diagnosis of multiple sclerosis (MS). However, the CVS has not yet been comprehensively studied in newly developing MS lesions. We aimed to identify the CVS profiles of new white matter lesions in patients with MS followed over time and investigate demographic and clinical risk factors associated with new CVS+ or CVS− lesion development.MethodsIn this retrospective longitudinal cohort study, adults from the NIH MS Natural History Study were considered for inclusion. Participants with new T2 or enhancing lesions were identified through review of the radiology report and/or longitudinal subtraction imaging. Each new lesion was evaluated for the CVS. Clinical characteristics were identified through chart review.ResultsA total of 153 adults (95 relapsing-remitting MS, 27 secondary progressive MS, 16 primary progressive MS, 5 clinically isolated syndrome, and 10 healthy; 67% female) were included. Of this cohort, 96 had at least 1 new T2 or contrast-enhancing lesion during median 3.1 years (Q1–Q3: 0.7–6.3) of follow-up; lesions eligible for CVS evaluation were found in 62 (65%). Of 233 new CVS-eligible lesions, 159 (68%) were CVS+, with 30 (48%) individuals having only CVS+, 12 (19%) only CVS−, and 20 (32%) both CVS+ and CVS− lesions. In gadolinium-enhancing (Gd+) lesions, the CVS+ percentage increased from 102/152 (67%) at the first time point where the lesion was observed, to 92/114 (82%) after a median follow-up of 2.8 years. Younger age (OR = 0.5 per 10-year increase, 95% CI = 0.3–0.8) and higher CVS+ percentage at baseline (OR = 1.4 per 10% increase, 95% CI = 1.1–1.9) were associated with increased likelihood of new CVS+ lesion development.DiscussionIn a cohort of adults with MS followed over a median duration of 3 years, most newly developing T2 or enhancing lesions were CVS+ (68%), and nearly half (48%) developed new CVS+ lesions only. Importantly, the effects of edema and T2 signal changes can obscure small veins in Gd+ lesions; therefore, caution and follow-up is necessary when determining their CVS status.Trial Registration InformationClinical trial registration number NCT00001248.Classification of EvidenceThis study provides Class III evidence that younger age and higher CVS+ percentage at baseline are associated with new CVS+ lesion development.

Objective: To obtain dense spatiotemporal measurements of brain deformation from two distinct but complementary head motion experiments: linear and rotational accelerations. Methods: This study introduces a strategy for integrating harmonic phase analysis of tagged magnetic resonance imaging (MRI) and finite-element models to extract mechanically representative deformation measurements. The method was calibrated using simulated as well as experimental data, demonstrated in a phantom including data with image artifacts, and used to measure brain deformation in human volunteers undergoing rotational and linear acceleration. Results: Evaluation methods yielded a displacement error of 1.1 mm compared to human observers and strain errors between ${\text{0.1}}\pm {\text{0.2}}{\% \,(\text{mean}\pm \text{std}.\,\text{dev.)}}$ for linear acceleration and ${\text{0.7}}\pm {\text{0.3}}\% $ for rotational acceleration. This study also demonstrates an approach that can reduce error by 86% in the presence of corrupted data. Analysis of results shows consistency with 2-D motion estimation, agreement with external sensors, and the expected physical behavior of the brain. Conclusion: Mechanical regularization is useful for obtaining dense spatiotemporal measurements of in vivo brain deformation under different loading regimes. Significance: The measurements suggest that the brain's 3-D response to mild accelerations includes distinct patterns observable using practical MRI resolutions. This type of measurement can provide validation data for computer models for the study of traumatic brain injury.

Magnetic resonance imaging (MRI) is a powerful tool for visualizing traumatic brain injury(TBI)-related lesions. Trauma-induced encephalomalacia is frequently identified by its hyperintense appearance on fluid-attenuated inversion recovery (FLAIR) sequences. In addition to parenchymal lesions, TBI commonly results in cerebral microvascular injury, but its anatomical relationship to parenchymal encephalomalacia is not well characterized. The current study utilized a multi-modal MRI protocol to assess microstructural tissue integrity (by mean diffusivity [MD] and fractional aniosotropy [FA]) and altered vascular function (by cerebral blood flow [CBF] and cerebral vascular reactivity [CVR]) within regions of visible encephalomalacia and normal appearing tissue in 27 chronic TBI (minimum 6 months post-injury) subjects. Fifteen subjects had visible encephalomalacias whereas 12 did not have evident lesions on MRI. Imaging from 14 age-matched healthy volunteers were used as controls. CBF was assessed by arterial spin labeling (ASL) and CVR by measuring the change in blood-oxygen-level–dependent (BOLD) MRI during a hypercapnia challenge. There was a significant reduction in FA, CBF, and CVR with a complementary increase in MD within regions of FLAIR-visible encephalomalacia (p

Improvements in resolution, pulse sequences, and access to scanners allow for increasingly subtle changes to be detected in magnetic resonance images, while increasing the volume of data to be interpreted by radiologists. The burden of reviewing this vast amount of information is further exacerbated by a lack of spatial alignment between images as well as the potential presence of multifocal pathologies. In this work, we describe an informatics approach for incorporating difference images into clinical workflows for highlighting 1) changes in anatomy over time, or 2) contrast enhanced regions following administration of gadolinium. These difference images enhance features of interest while suppressing stable features, providing efficient detection of changes over time, whether due to treatment, age, medical intervention, or progression of disease. Our pipeline is initiated by the physician sending images to a DICOM receiver, which then applies spatial transformation, bias-field correction, and intensity normalization before returning the newly generated difference images to the clinical PACS for review. This pipeline has been used to process 2683 neuroimaging data sets for evaluating disease progression and contrast enhancement in patients with brain lesions. We demonstrate that difference images quantitatively improve contrast-to-noise ratio (CNR) of new lesions while reducing clutter, which should translate to improved accuracy and efficiency in radiological interpretations.

Natural modes and frequencies of three-dimensional (3D) deformation of the human brain were identified from in vivo tagged magnetic resonance images (MRI) acquired dynamically during transient mild acceleration of the head. Twenty 3D strain fields, estimated from tagged MRI image volumes in 19 adult subjects, were analyzed using dynamic mode decomposition (DMD). These strain fields represented dynamic, 3D brain deformations during constrained head accelerations, either involving rotation about the vertical axis of the neck or neck extension. DMD results reveal fundamental oscillatory modes of deformation at damped frequencies near 7 Hz (in neck rotation) and 11 Hz (in neck extension). Modes at these frequencies were found consistently among all subjects. These characteristic features of 3D human brain deformation are important for understanding the response of the brain in head impacts and provide valuable quantitative criteria for the evaluation and use of computer models of brain mechanics.

Traumatic Brain Injury (TBI) is associated with both diffuse axonal injury (DAI) and diffuse vascular injury (DVI), which result from inertial shearing forces. These terms are often used interchangeably, but the spatial relationships between DAI and DVI have not been carefully studied. Multimodal magnetic resonance imaging (MRI) can help distinguish these injury mechanisms: diffusion tensor imaging (DTI) provides information about axonal integrity, while arterial spin labeling (ASL) can be used to measure cerebral blood flow (CBF), and the reactivity of the Blood Oxygen Level Dependent (BOLD) signal to a hypercapnia challenge reflects cerebrovascular reactivity (CVR). Subjects with chronic TBI (n = 27) and healthy controls (n = 14) were studied with multimodal MRI. Mean values of mean diffusivity (MD), fractional anisotropy (FA), CBF, and CVR were extracted for pre-determined regions of interest (ROIs). Normalized z-score maps were generated from the pool of healthy controls. Abnormal ROIs in one modality were not predictive of abnormalities in another. Approximately 9-10% of abnormal voxels for CVR and CBF also showed an abnormal voxel value for MD, while only 1% of abnormal CVR and CBF voxels show a concomitant abnormal FA value. These data indicate that DAI and DVI represent two distinct TBI endophenotypes that are spatially independent.

Measurements of how the brain deforms in response to head impact are often obtained in human cadavers because loading conditions can be investigated that are on the order of those that could lead to injury. How the deformation response of the postmortem human brain differs from the live human brain remains unanswered. In this study, we used tagged magnetic resonance imaging (MRI) to measure brain deformation in two human cadavers in response to non-injurious head impact. We then compared these strain fields to those obtained in 10 healthy volunteers (HV) during the same neck rotation under similar loading conditions. Our results showed a number of similarities, such as a similar magnitude of the largest maximum principal strain values and high values within the cortex, and some differences, such as high strain concentrations near the falx, lower median strain values at each time point, and an increased frequency of oscillation in cadavers. These initial results provide insight into potential differences in the deformation response of the brain in cadavers compared to HVs, though additional experiments are needed to further investigate the effect of age and atrophy on the deformation response.

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/.

High resolution magnetic resonance (MR) images are desired in many clinical and research applications. Acquiring such images with high signal-to-noise (SNR), however, can require a long scan duration, which is difficult for patient comfort, is more costly, and makes the images susceptible to motion artifacts. A very common practical compromise for both 2D and 3D MR imaging protocols is to acquire volumetric MR images with high in-plane resolution, but lower through-plane resolution. In addition to having poor resolution in one orientation, 2D MRI acquisitions will also have aliasing artifacts, which further degrade the appearance of these images. This paper presents an approach SMORE 1 based on convolutional neural networks (CNNs) that restores image quality by improving resolution and reducing aliasing in MR images. 2 This approach is self-supervised, which requires no external training data because the high-resolution and low-resolution data that are present in the image itself are used for training. For 3D MRI, the method consists of only one self-supervised super-resolution (SSR) deep CNN that is trained from the volumetric image data. For 2D MRI, there is a self-supervised anti-aliasing (SAA) deep CNN that precedes the SSR CNN, also trained from the volumetric image data. Both methods were evaluated on a broad collection of MR data, including filtered and downsampled images so that quantitative metrics could be computed and compared, and actual acquired low resolution images for which visual and sharpness measures could be computed and compared. The super-resolution method is shown to be visually and quantitatively superior to previously reported methods.

Highlights • PML has characteristic dynamic changes in brain and lesion volume on MRI. • JCnet is an automated method for brain atrophy and lesion segmentation in PML. • JCnet improves PML lesion segmentation accuracy compared to conventional methods. • JCnet can accurately track PML lesion changes over time., 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 meninges serve as a functional barrier surrounding the brain, critical to the immune response, and can be compromised following head trauma. Meningeal enhancement can be detected on contrast-enhanced MRI in patients presenting with acute traumatic brain injury, even when head CT is negative. Following head trauma, gadolinium-based contrast appears to extravasate from the vasculature, enhancing the dura within minutes, and later permeates the subarachnoid space. The aims of this study were to characterize the initial kinetics of the uptake of contrast agent after injury and the delayed redistribution of contrast enhancement in the subarachnoid space in hyperacute patients. Neuroimaging was obtained prospectively in two large ongoing observational studies of patients aged 18 years or older presenting to the emergency department with suspected acute head injury. Dynamic contrast-enhanced MRI studies in a cohort of consecutively enrolling patients with mild traumatic brain injury (n = 36) determined that the kinetic half-life of dural-related meningeal enhancement was 1.3 ± 0.6 min (95% enhancement within 6 min). The extravasation of contrast into the subarachnoid space was investigated in a cohort of CT negative mild traumatic brain injury patients initially imaged within 6 h of injury (hyperacute) who subsequently underwent a delayed MRI, with no additional contrast administration, several hours after the initial MRI. Of the 32 patients with delayed post-contrast imaging, 18 (56%) had conspicuous expansion of the contrast enhancement into the subarachnoid space, predominantly along the falx and superior sagittal sinus. Patients negative for traumatic meningeal enhancement on initial hyperacute MRI continued to have no evidence of meningeal enhancement on the delayed MRI. These studies demonstrate that (i) the initial enhancement of the traumatically injured meninges occurs within minutes of contrast injection, suggesting highly permeable meningeal vasculature, and that (ii) contrast in the meninges redistributes within the subarachnoid space over the period of hours, suggesting a compromise in the blood–brain and/or blood-cerebrospinal barriers. Data from the parent study indicate that up to one in two patients with mild traumatic brain injury have traumatic brain injury on acute (

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.

Image synthesis and superresolution have a long history of research and use in medical imaging. Recently, with the emergence of powerful machine learning methods, both processes have risen to new levels of capability. This chapter presents some historical and technical background in each of these areas and describes some of the new approaches. Magnetic resonance imaging is the primary modality considered.

Background and purpose Magnetic resonance imaging (MRI) is crucial for in vivo detection and characterization of white matter lesions (WMLs) in multiple sclerosis. While WMLs have been studied for over two decades using MRI, automated segmentation remains challenging. Although the majority of statistical techniques for the automated segmentation of WMLs are based on single imaging modalities, recent advances have used multimodal techniques for identifying WMLs. Complementary modalities emphasize different tissue properties, which help identify interrelated features of lesions. Methods Method for Inter-Modal Segmentation Analysis (MIMoSA), a fully automatic lesion segmentation algorithm that utilizes novel covariance features from intermodal coupling regression in addition to mean structure to model the probability lesion is contained in each voxel, is proposed. MIMoSA was validated by comparison with both expert manual and other automated segmentation methods in two datasets. The first included 98 subjects imaged at Johns Hopkins Hospital in which bootstrap cross-validation was used to compare the performance of MIMoSA against OASIS and LesionTOADS, two popular automatic segmentation approaches. For a secondary validation, a publicly available data from a segmentation challenge were used for performance benchmarking. Results In the Johns Hopkins study, MIMoSA yielded average Sorensen-Dice coefficient (DSC) of .57 and partial AUC of .68 calculated with false positive rates up to 1%. This was superior to performance using OASIS and LesionTOADS. The proposed method also performed competitively in the segmentation challenge dataset. Conclusion MIMoSA resulted in statistically significant improvements in lesion segmentation performance compared with LesionTOADS and OASIS, and performed competitively in an additional validation study.

Being able to focus on a complex task and inhibit unwanted actions or interfering information (i.e., inhibitory control) are essential human cognitive abilities. However, it remains unknown the extent to which mild traumatic brain injury (mTBI) may impact these critical functions. In this study, seventeen patients and age-matched healthy controls (HC) performed a variant of the Stroop task and attention-demanding 4-choice response tasks (4CRT) with identical stimuli but two contexts: one required only routine responses and the other with occasional response conflicts. The results showed that mTBI patients performed equally well as the HC when the 4CRT required only routine responses. However, when the task conditions included occasional response conflicts, mTBI patients with even a single concussion showed a significant slow-down in all responses and higher error rates relative to the HC. Results from event-related functional magnetic resonance imaging (efMRI) revealed altered neural activity in the mTBI patients in the cerebellum-thalamo-cortical and the fronto-basal-ganglia networks regulating inhibitory control. These results suggest that even without apparent difficulties in performing complex attention-demanding but routine tasks, patients with mTBI may experience long-lasting deficits in regulating inhibitory control when situations call for rapid conflict resolutions.

Diffuse axonal injury (DAI) is a debilitating consequence of traumatic brain injury (TBI) attributed to abnormal stretching of axons caused by blunt head trauma or acceleration of the head. We developed an anatomically accurate, subject-specific, three-dimensional (3D) computational model of the human brain, and used it to study the dynamic deformations in the substructures of the brain when the head is subjected to rotational accelerations. The computational head models use anatomy and morphology of the white matter fibers obtained using MRI. Subject-specific full-field shearing motions in live human brains obtained through a recently developed tagged MRI imaging technique are then used to validate the models by comparing the measured and predicted heterogeneous dynamic mechanical response of the brain. These results are used to elucidate the dynamics of local shearing deformations in the brain substructures caused by rotational acceleration of the head. Our work demonstrates that the rotational dynamics of the brain has a timescale of ∼100 ms as determined by the shearing wave speeds, and thus the injuries associated with rotational accelerations likely occur over these time scales. After subject-specific validation using the live human subject data, a representative subject-specific head model is used to simulate a real life scenario that resulted in a concussive injury. Results suggest that regions of the brain, in the form of a toroid, encompassing the white matter, the cortical gray matter, and outer parts of the limbic system have a higher susceptibility to injury under axial rotations of the head.

Purpose To explore (i) the variability of upper cervical cord area (UCCA) measurements from volumetric brain 3D T1-weighted scans related to gradient nonlinearity (GNL) and subject positioning; (ii) the effect of vendor-implemented GNL corrections; and (iii) easily applicable methods that can be used to retrospectively correct data. Methods A multiple sclerosis patient was scanned at seven sites using 3T MRI scanners with the same 3D T1-weighted protocol without GNL-distortion correction. Two healthy subjects and a phantom were additionally scanned at a single site with varying table positions. The 2D and 3D vendor-implemented GNL-correction algorithms and retrospective methods based on (i) phantom data fit, (ii) normalization with C2 vertebral body diameters, and (iii) the Jacobian determinant of nonlinear registrations to a template were tested. Results Depending on the positioning of the subject, GNL introduced up to 15% variability in UCCA measurements from volumetric brain T1-weighted scans when no distortion corrections were used. The 3D vendor-implemented correction methods and the three proposed methods reduced this variability to less than 3%. Conclusions Our results raise awareness of the significant impact that GNL can have on quantitative UCCA studies, and point the way to prospectively and retrospectively managing GNL distortions in a variety of settings, including clinical environments. Magn Reson Med, 2017. © 2017 International Society for Magnetic Resonance in Medicine.

We assessed the utility of a brief MRI protocol, appropriate for the acute setting, to detect acute traumatic brain injury (TBI) in patients with suspected mild TBI (mTBI) and distinguish traumatic from nontraumatic brain injury by comparing trauma with nontrauma patients. Twenty-two patients with suspected mTBI were included in this exploratory study over a period of 9 months. Median time from injury to MR scanning was 5.4 h (interquartile range 3.6-15.3). To determine the specificity of certain findings for TBI, 61 patients presenting with suspected minor acute stroke were included as a comparative group using the same MRI methods. A selected series of MRI sequences (diffusion-weighted imaging, fluid attenuated inversion recovery [FLAIR], and T2* weighted) were independently evaluated by two neuroradiologists blinded to clinical diagnosis, for presence of specific findings. In a separate session, all cases in which at least one MRI sequence above was positive were classified as TBI, stroke, or indeterminate. Intracranial MRI abnormalities were observed in 47 (57%) of the 83 studied patients. Based on findings on MRI, 12 (55%) of 22 suspected mTBI patients were classified as having traumatic injury. Nine (47%) of the 19 suspected mTBI patients with a negative CT had findings on MRI. Abnormalities on MRI consistent with trauma were observed most frequently on postcontrast FLAIR (83%) and T2*-weighted (58%) sequences. We demonstrated the ability of a fast MRI protocol to identify trauma-related abnormalities not seen on CT, and differentiate acute trauma from nonspecific chronic disease in a blinded cohort of mTBI patients.

Head impact can cause traumatic brain injury (TBI) through axonal overstretch or subsequent inflammation and understanding the biomechanics of the impact event is useful for TBI prevention research. Tagged magnetic resonance imaging (MRI) acquired during a mild-acceleration impact has enabled measurement and visualization of brain deformation in vivo. However, measurements using MRI are subject to error, and having independent validation while imaging in vivo is very difficult. Thus, characterizing the accuracy of these measurements needs to be done in a separate experiment using a phantom where a gold standard is available. This study describes a method for error quantification using a calibration phantom compatible with MRI and high-speed video (the gold standard). During linear acceleration, the maximum shear strain (MSS) in the phantom ranged from 0 to 12%, which is similar to in vivo brain deformation at a similar acceleration. The mean displacement error against video was 0.3 ± 0.3 mm, and the MSS error was 1.4 ± 0.3%. To match resolutions, video data was filtered temporally using an averaging filter. Compared to the unfiltered results, resolution matching improved the agreement between MRI and video results by 15%. In conclusion, tagged MRI analysis compares well to video data provided that resolutions are matched—a finding that is also applicable when using MRI to validate simulations.

Total brain white matter lesion (WML) volume is the most widely established magnetic resonance imaging (MRI) outcome measure in studies of multiple sclerosis (MS). To estimate WML volume, there are a number of automatic segmentation methods available, yet manual delineation remains the gold standard approach. Automatic approaches often yield a probability map to which a threshold is applied to create lesion segmentation masks. Unfortunately, few approaches systematically determine the threshold employed; many methods use a manually selected threshold, thus introducing human error and bias into the automated procedure. In this study, we propose and validate an automatic thresholding algorithm, Thresholding Approach for Probability Map Automatic Segmentation in Multiple Sclerosis (TAPAS), to obtain subject-specific threshold estimates for probability map automatic segmentation of T2-weighted (T2) hyperintense WMLs. Using multimodal MRI, the proposed method applies an automatic segmentation algorithm to obtain probability maps. We obtain the true subject-specific threshold that maximizes the Sorensen-Dice similarity coefficient (DSC). Then the subject-specific thresholds are modeled on a naive estimate of volume using a generalized additive model. Applying this model, we predict a subject-specific threshold in data not used for training. We ran a Monte Carlo-resampled split-sample cross-validation (100 validation sets) using two data sets: the first obtained from the Johns Hopkins Hospital (JHH) on a Philips 3 Tesla (3T) scanner (n = 94) and a second collected at the Brigham and Women's Hospital (BWH) using a Siemens 3T scanner (n = 40). By means of the proposed automated technique, in the JHH data we found an average reduction in subject-level absolute error of 0.1 mL per one mL increase in manual volume. Using Bland-Altman analysis, we found that volumetric bias associated with group-level thresholding was mitigated when applying TAPAS. The BWH data showed similar absolute error estimates using group-level thresholding or TAPAS likely since Bland-Altman analyses indicated no systematic biases associated with group or TAPAS volume estimates. The current study presents the first validated fully automated method for subject-specific threshold prediction to segment brain lesions.

Total brain white matter lesion (WML) volume is the most widely established magnetic resonance imaging (MRI) outcome measure in studies of multiple sclerosis (MS). To estimate WML volume, there are a number of automatic segmentation methods, yet, manual delineation remains the gold standard approach. These approaches often yield a probability map to which a threshold is applied to create lesion segmentation masks. Unfortunately, few approaches systematically determine the threshold employed; many methods use a manually selected threshold, thus introducing human error and bias into the automated procedure. In this study, we propose and validate an automatic thresholding algorithm, Thresholding Approach for Probability Map Automatic Segmentation in Multiple Sclerosis (TAPAS), to obtain subject-specific threshold estimates for probability map automatic segmentation of T2-weighted (T2) hyperintense WMLs. Using multimodal MRI, the proposed method applies an automatic segmentation algorithm to obtain probability maps. We obtain the true subject-specific threshold that maximizes Sørensen-Dice Similarity Coefficient (DSC). Then the subject-specific thresholds are modeled on a naive estimate of volume using a general additive model. Applying this model, we predict a subject-specific threshold in data not used for training. We ran a Monte Carlo-resampled split-sample cross-validation (100 validation sets) using two data sets: the first obtained from the Johns Hopkins Hospital (JHH) on a Philips 3 Tesla (3T) scanner (n = 94) and a second collected at the Brigham and Women’s Hospital (BWH) using a Siemens 3T scanner (n = 40). By means of the proposed automated technique, in the JHH data, we found an average reduction in subject-level absolute error of 0.1 mL per one mL increase in manual volume. Using Bland-Altman analysis, we found that volumetric bias associated with group-level thresholding is mitigated when applying TAPAS. The BWH data showed similar absolute error estimates using group-level thresholding or TAPAS likely since Bland-Altman analyses indicate no systematic biases associated with group or TAPAS volume estimates. The current study presents the first validated fully automated method for subject-specific threshold prediction to segment brain lesions.

Chordoma is a rare type of tumor that usually appears in the bone near the spinal cord and skull base. Due to their location in the skull base and diverse appearance in size and shape, automatic segmentation of chordoma tumors from magnetic resonance images (MRI) is a challenging task. In addition, similar MR intensity distributions of different anatomical regions, specifically sinuses, make the segmentation task from MRI more challenging. In comparison, most of the state-of-the-art lesion segmentation methods are designed to segment pathologies inside the brain. In this work, we propose an automatic chordoma segmentation framework using two cascaded 3D convolutional neural networks (CNN) via an auto-context model. While the first network learns to detect all potential tumor voxels, the second network fine-tunes the classifier to distinguish true tumor voxels from the false positives detected by the first network. The proposed method is evaluated using multi-contrast MR images of 22 longitudinal scans from 8 patients. Preliminary results showed a linear correlation of 0.71 between the detected and manually outlined tumor volumes, compared to 0.40 for a random forest (RF) based method. Furthermore, the response of tumor growth over time, i.e. increasing, decreasing, or stable, is evaluated according to the response evaluation criteria in solid tumors with an outcome of 0.26 kappa coefficient, compared to 0.13 for the RF based method.

Background: The effects of disease-modifying therapies (DMTs) on region-specific brain atrophy in multiple sclerosis (MS) are unclear. Objective: To determine the effects of higher versus lower efficacy DMTs on rates of brain substructure atrophy in MS. Methods: A non-randomized, observational cohort of people with MS followed with annual brain magnetic resonance imaging (MRI) was evaluated retrospectively. Whole brain, subcortical gray matter (GM), cortical GM, and cerebral white matter (WM) volume fractions were obtained. DMTs were categorized as higher (DMT-H: natalizumab and rituximab) or lower (DMT-L: interferon-beta and glatiramer acetate) efficacy. Follow-up epochs were analyzed if participants had been on a DMT for ⩾6 months prior to baseline and had at least one follow-up MRI while on DMTs in the same category. Results: A total of 86 DMT epochs (DMT-H: n = 32; DMT-L: n = 54) from 78 participants fulfilled the study inclusion criteria. Mean follow-up was 2.4 years. Annualized rates of thalamic (−0.15% vs −0.81%; p = 0.001) and putaminal (−0.27% vs −0.73%; p = 0.001) atrophy were slower during DMT-H compared to DMT-L epochs. These results remained significant in multivariate analyses including demographics, clinical characteristics, and T2 lesion volume. Conclusion: DMT-H treatment may be associated with slower rates of subcortical GM atrophy, especially of the thalamus and putamen. Thalamic and putaminal volumes are promising imaging biomarkers in MS.

Magnetic resonance imaging (MRI) is a flexible medical imaging modality that often lacks reproducibility between protocols and scanners. It has been shown that even when care is taken to standardize acquisitions, any changes in hardware, software, or protocol design can lead to differences in quantitative results. This greatly impacts the quantitative utility of MRI in multi-site or long-term studies, where consistency is often valued over image quality. We propose a method of contrast harmonization, called DeepHarmony, which uses a U-Net-based deep learning architecture to produce images with consistent contrast. To provide training data, a small overlap cohort (n = 8) was scanned using two different protocols. Images harmonized with DeepHarmony showed significant improvement in consistency of volume quantification between scanning protocols. A longitudinal MRI dataset of patients with multiple sclerosis was also used to evaluate the effect of a protocol change on atrophy calculations in a clinical research setting. The results show that atrophy calculations were substantially and significantly affected by protocol change, whereas such changes have a less significant effect and substantially reduced overall difference when using DeepHarmony. This establishes that DeepHarmony can be used with an overlap cohort to reduce inconsistencies in segmentation caused by changes in scanner protocol, allowing for modernization of hardware and protocol design in long-term studies without invalidating previously acquired data.

Empirical studies of taste function in multiple sclerosis (MS) are rare. Moreover, a detailed assessment of whether quantitative measures of taste function correlate with the punctate and patchy myelin-related lesions found throughout the CNS of MS patients has not been made. We administered a 96-trial test of sweet (sucrose), sour (citric acid), bitter (caffeine) and salty (NaCl) taste perception to the left and right anterior (CN VII) and posterior (CN IX) tongue regions of 73 MS patients and 73 matched controls. The number and volume of lesions were assessed using quantitative MRI in 52 brain regions of 63 of the MS patients. Taste identification scores were significantly lower in the MS patients for sucrose (p = 0.0002), citric acid (p = 0.0001), caffeine (p = 0.0372) and NaCl (p = 0.0004) and were present in both anterior and posterior tongue regions. The percent of MS patients with identification scores falling below the 5th percentile of controls was 15.07 % for caffeine, 21.9 % for citric acid, 24.66 % for sucrose, and 31.50 % for NaCl. Such scores were inversely correlated with lesion volumes in the temporal, medial frontal, and superior frontal lobes, and with the number of lesions in the left and right superior frontal lobes, right anterior cingulate gyrus, and left parietal operculum. Regardless of the subject group, women outperformed men on the taste measures. These findings indicate that a sizable number of MS patients exhibit taste deficits that are associated with MS-related lesions throughout the brain., U.S. National Institutes of Health ( RO1 DC 02974 -- R01 NS 37172 -- R01 NS070906)

Longitudinal analysis of magnetic resonance images of the human brain provides knowledge of brain changes during both normal aging as well as the progression of many diseases. Previous longitudinal segmentation methods have either ignored temporal information or have incorporated temporal consistency constraints within the algorithm. In this work, we assume that some anatomical brain changes can be explained by temporal transitions in image intensities. Once the images are aligned in the same space, the intensities of each scan at the same voxel constitute a temporal (or 4D) intensity trend at that voxel. Temporal intensity variations due to noise or other artifacts are corrected by a 4D intensity-based filter that smooths the intensity values where appropriate, while preserving real anatomical changes such as atrophy. Here smoothing refers to removal of sudden changes or discontinuities in intensities. Images processed with the 4D filter can be used as a pre-processing step to any segmentation method. We show that such a longitudinal pre-processing step produces robust and consistent longitudinal segmentation results, even when applying 3D segmentation algorithms. We compare with state-of-the-art 4D segmentation algorithms. Specifically, we experimented on three longitudinal datasets containing 4–12 time-points, and showed that the 4D temporal filter is more robust and has more power in distinguishing between healthy subjects and those with dementia, mild cognitive impairment, as well as different phenotypes of multiple sclerosis., Highlights • Longitudinal analysis of MR images of the human brain provides knowledge about brain changes during aging and diseases. • Unlike previous methods, some temporal anatomical changes in the brain has been modeled by MR intensity change. • A 4D intensity based filter is proposed as a pre-processing step to 3D segmentation methods. • The 4D filter ensures consistency in segmentation while retaining sensitivity to true anatomical changes. • Experiments show improved atrophy detection in presence of dementia, cognitive impairment and different phenotypes of multiple sclerosis.

On average, African Americans with multiple sclerosis demonstrate higher inflammatory disease activity, faster disability accumulation, greater visual dysfunction, more pronounced brain tissue damage and higher lesion volume loads compared to Caucasian Americans with multiple sclerosis. Neurodegeneration is an important component of multiple sclerosis, which in part accounts for the clinical heterogeneity of the disease. Brain atrophy appears to be widespread, although it is becoming increasingly recognized that regional substructure atrophy may be of greater clinical relevance. Patient race (within the limitations of self-identified ancestry) is regarded as an important contributing factor. However, there is a paucity of studies examining differences in neurodegeneration and brain substructure volumes over time in African Americans relative to Caucasian American patients. Optical coherence tomography is a non-invasive and reliable tool for measuring structural retinal changes. Recent studies support its utility for tracking neurodegeneration and disease progression in vivo in multiple sclerosis. Relative to Caucasian Americans, African American patients have been found to have greater retinal structural injury in the inner retinal layers. Increased thickness of the inner nuclear layer and the presence of microcystoid macular pathology at baseline predict clinical and radiological inflammatory activity, although whether race plays a role in these changes has not been investigated. Similarly, assessment of outer retinal changes according to race in multiple sclerosis remains incompletely characterized. Twenty-two African Americans and 60 matched Caucasian Americans with multiple sclerosis were evaluated with brain MRI, and 116 African Americans and 116 matched Caucasian Americans with multiple sclerosis were monitored with optical coherence tomography over a mean duration of 4.5 years. Mixed-effects linear regression models were used in statistical analyses. Grey matter (-0.9%/year versus -0.5%: P =0.02), white matter (-0.7%/year versus -0.3%: P =0.04) and nuclear thalamic (-1.5%/year versus -0.7%/year: P =0.02) atrophy rates were approximately twice as fast in African Americans. African Americans also exhibited higher proportions of microcystoid macular pathology (12.1% versus 0.9%, P =0.001). Retinal nerve fibre layer (-1.1% versus -0.8%: P =0.02) and ganglion cell+ inner plexiform layer (-0.7%/year versus -0.4%/year: P =0.01) atrophy rates were faster in African versus Caucasian Americans. African Americans on average exhibited more rapid neurodegeneration than Caucasian Americans and had significantly faster brain and retinal tissue loss. These results corroborate the more rapid clinical progression reported to occur, in general, in African Americans with multiple sclerosis and support the need for future studies involving African Americans in order to identify individual differences in treatment responses in multiple sclerosis.

Automatic brain extraction or skull stripping from magnetic resonance images (MRI) is an important pre-processing step in many image processing pipelines. Most skull stripping methods are optimized for normal brains and applicable to single T 1 -w MR images. However, other contrasts, such as T 2 , can provide complementary information about the boundary. This is especially true in the presence of traumatic brain injury (TBI) and other diseases, where lesions can confound boundary definitions. In this paper, we propose a deep learning based framework to extract intracranial tissues from multi-contrast MR images in the presence of TBI. Our approach is based on state-of-the-art convolutional neural network architecture to learn a transformation from multi-contrast atlas MR images to their stripping masks without using any deformable registration. An advantage of our framework is that it can be applied to different species. We applied our approach to 19 human patients with mild to severe TBI, as well as 16 normal mice images, and another 10 mice brains with TBI. We compared the approach with 3 separate state-of-the-art human and rodent brain extraction methods. Using only a few manually delineated atlases, we showed significant improvement in brain extraction accuracy in both healthy and pathological human and rodent images.

Dizziness, postural instability, and ataxia are among the most debilitating symptoms of multiple sclerosis (MS), reflecting, in large part, dysfunctional integration of visual, somatosensory, and vestibular sensory cues. However, the role of MS-related supratentorial lesions in producing such symptoms is poorly understood. In this study, motor control test (MCT) and dynamic sensory organization test (SOT) scores of 58 MS patients were compared to those of 72 healthy controls; correlations were determined between the MS scores of 49 patients and lesion volumes within 26 brain regions. Depending upon platform excursion direction and magnitude, MCT latencies, which were longer in MS patients than controls (p

Magnetic Resonance (MR) imaging allows the acquisition of images with different contrast properties depending on the acquisition protocol and the magnetic properties of tissues. Many MR brain image processing techniques, such as tissue segmentation, require multiple MR contrasts as inputs, and each contrast is treated differently. Thus it is advantageous to automate the identification of image contrasts for various purposes, such as facilitating image processing pipelines, and managing and maintaining large databases via content-based image retrieval (CBIR). Most automated CBIR techniques focus on a two-step process: extracting features from data and classifying the image based on these features. We present a novel 3D deep convolutional neural network (CNN)-based method for MR image contrast classification. The proposed CNN automatically identifies the MR contrast of an input brain image volume. Specifically, we explored three classification problems: (1) identify T1-weighted (T1-w), T2-weighted (T2-w), and fluid-attenuated inversion recovery (FLAIR) contrasts, (2) identify pre vs post-contrast T1, (3) identify pre vs post-contrast FLAIR. A total of 3418 image volumes acquired from multiple sites and multiple scanners were used. To evaluate each task, the proposed model was trained on 2137 images and tested on the remaining 1281 images. Results showed that image volumes were correctly classified with 97.57% accuracy., Comment: Github: https://github.com/sremedios/phinet

High resolution magnetic resonance (MR) images are desired in many clinical applications, yet acquiring such data with an adequate signal-to-noise ratio requires a long time, making them costly and susceptible to motion artifacts. A common way to partly achieve this goal is to acquire MR images with good in-plane resolution and poor through-plane resolution (i.e., large slice thickness). For such 2D imaging protocols, aliasing is also introduced in the through-plane direction, and these high-frequency artifacts cannot be removed by conventional interpolation. Super-resolution (SR) algorithms which can reduce aliasing artifacts and improve spatial resolution have previously been reported. State-of-the-art SR methods are mostly learning-based and require external training data consisting of paired low resolution (LR) and high resolution (HR) MR images. However, due to scanner limitations, such training data are often unavailable. This paper presents an anti-aliasing (AA) and self super-resolution (SSR) algorithm that needs no external training data. It takes advantage of the fact that the in-plane slices of those MR images contain high frequency information. Our algorithm consists of three steps: (1) We build a self AA (SAA) deep network followed by (2) an SSR deep network, both of which can be applied along different orientations within the original images, and (3) recombine the multiple orientations output from Steps 1 and 2 using Fourier burst accumulation. We perform our SAA+SSR algorithm on a diverse collection of MR data without modification or preprocessing other than N4 inhomogeneity correction, and demonstrate significant improvement compared to competing SSR methods.

Objective: Pathology in both cortex and deep gray matter contribute to disability in multiple sclerosis (MS). We used the increased signal-to-noise ratio of 7-tesla (7T) MRI to visualize small lesions within the thalamus and to relate this to clinical information and cortical lesions. Methods: We obtained 7T MRI scans on 34 MS cases and 15 healthy volunteers. Thalamic lesion number and volume were related to demographic data, clinical disability measures, and lesions in cortical gray matter. Results: Thalamic lesions were found in 24/34 of MS cases. Two lesion subtypes were noted: discrete, ovoid lesions, and more diffuse lesional areas lining the periventricular surface. The number of thalamic lesions was greater in progressive MS compared to relapsing–remitting (mean ±SD, 10.7 ±0.7 vs. 3.0 ±0.7, respectively, p < 0.001). Thalamic lesion burden (count and volume) correlated with EDSS score and measures of cortical lesion burden, but not with white matter lesion burden or white matter volume. Conclusions: Using 7T MRI allows identification of thalamic lesions in MS, which are associated with disability, progressive disease, and cortical lesions. Thalamic lesion analysis may be a simpler, more rapid estimate of overall gray matter lesion burden in MS.

Integrated PET/MR systems are becoming increasingly popular in clinical and research applications. Quantitative PET reconstruction requires correction for γ-photon attenuations using an attenuation coefficient map (μ map) that is a measure of the electron density. One challenge of PET/MR, in contrast to PET/CT, lies in the accurate computation of μ maps. Unlike CT, MR imaging measures physical properties not directly related to electron density. Previous approaches have computed the attenuation coefficients using a segmentation of MR images or using deformable registration of atlas CT images to the space of the subject MR images. Methods: In this work, we propose a patch-based method to generate whole-head μ maps from ultrashort echo-time (UTE) MR imaging sequences. UTE images are preferred to other MR sequences because of the increased signal from bone. To generate a synthetic CT image, we use patches from a reference dataset, which consists of dual-echo UTE images and a coregistered CT scan from the same subject. Matching of patches between the reference and target images allows corresponding patches from the reference CT scan to be combined via a Bayesian framework. No registration or segmentation is required. Results: For evaluation, UTE, CT, and PET data acquired from 5 patients under an institutional review board–approved protocol were used. Another patient (with UTE and CT data only) was selected to be the reference to generate synthetic CT images for these 5 patients. PET reconstructions were attenuation-corrected using the original CT, our synthetic CT, Siemens Dixon-based μ maps, Siemens UTE-based μ maps, and deformable registration-based CT. Our synthetic CT–based PET reconstruction showed higher correlation (average ρ = 0.996, R2 = 0.991) to the original CT-based PET, as compared with the segmentation- and registration-based methods. Synthetic CT–based reconstruction had minimal bias (regression slope, 0.990), as compared with the segmentation-based methods (regression slope, 0.905). A peak signal-to-noise ratio of 35.98 dB in the reconstructed PET activity was observed, compared with 29.767, 29.34, and 27.43 dB for the Siemens Dixon-, UTE-, and registration-based μ maps. Conclusion: A patch-matching approach to synthesize CT images from dual-echo UTE images leads to significantly improved accuracy of PET reconstruction as compared with actual CT scans. The PET reconstruction is improved over segmentation- (Dixon and Siemens UTE) and registration-based methods, even in subjects with pathologic findings.

The falx cerebri and tentorium cerebelli are dural structures found in the brain. Due to the roles both structures play in constraining brain motion, the falx and tentorium must be identified and included in finite element models of the head to accurately predict brain dynamics during injury events. To date there has been very little research work on automatically segmenting these two structures, which is understandable given that their 1) thin structure challenges the resolution limits of in vivo 3D imaging, and 2) contrast with respect to surrounding tissue is low in standard magnetic resonance imaging. An automatic segmentation algorithm to find the falx and tentorium which uses the results of a multi-atlas segmentation and cortical reconstruction algorithm is proposed. Gray matter labels are used to find the location of the falx and tentorium. The proposed algorithm is applied to five datasets with manual delineations. 3D visualizations of the final results are provided, and Hausdorff distance (HD) and mean surface distance (MSD) is calculated to quantify the accuracy of the proposed method. For the falx, the mean HD is 43.84 voxels and the mean MSD is 2.78 voxels, with the largest errors occurring at the frontal inferior falx boundary. For the tentorium, the mean HD is 14.50 voxels and mean MSD is 1.38 voxels.

Understanding of in vivo brain biomechanical behavior is critical in the study of traumatic brain injury (TBI) mechanisms and prevention. Using tagged magnetic resonance imaging, we measured spatiotemporal brain deformations in 34 healthy human volunteers under mild angular accelerations of the head. Two-dimensional (2D) Lagrangian strains were examined throughout the brain in each subject. Strain metrics peaked shortly after contact with a padded stop, corresponding to the inertial response of the brain after head deceleration. Maximum shear strain of at least 3% was experienced at peak deformation by an area fraction (median±standard error) of 23.5±1.8% of cortical gray matter, 15.9±1.4% of white matter, and 4.0±1.5% of deep gray matter. Cortical gray matter strains were greater in the temporal cortex on the side of the initial contact with the padded stop and also in the contralateral temporal, frontal, and parietal cortex. These tissue-level deformations from a population of healthy volunteers provide the first in vivo measurements of full-volume brain deformation in response to known kinematics. Although strains differed in different tissue type and cortical lobes, no significant differences between male and female head accelerations or strain metrics were found. These cumulative results highlight important kinematic features of the brain's mechanical response and can be used to facilitate the evaluation of computational simulations of TBI.

The falx cerebri is a meningeal projection of dura in the brain, separating the cerebral hemispheres. It has stiffer mechanical properties than surrounding tissue and must be accurately segmented for building computational models of traumatic brain injury. In this work, we propose a method to segment the falx using T1-weighted magnetic resonance images (MRI) and susceptibility-weighted MRI (SWI). Multi-atlas whole brain segmentation is performed using the T1-weighted MRI and the gray matter cerebrum labels are extended into the longitudinal fissure using fast marching to find an initial estimate of the falx. To correct the falx boundaries, we register and then deform a set of SWI with manually delineated falx boundaries into the subject space. The continuous-STAPLE algorithm fuses sets of corresponding points to produce an estimate of the corrected falx boundary. Correspondence between points on the deformed falx boundaries is obtained using coherent point drift. We compare our method to manual ground truth, a multi-atlas approach without correction, and single-atlas approaches.

Tagged magnetic resonance imaging has been used for decades to observe and quantify motion and strain of deforming tissue. It is challenging to obtain three-dimensional (3D) motion estimates due to a tradeoff between image slice density and acquisition time. Typically, interpolation methods are used either to combine two-dimensional motion extracted from sparse slice acquisitions into 3D motion or to construct a dense volume from sparse acquisitions before image registration methods are applied. This paper proposes a new phase-based 3D motion estimation technique that first computes harmonic phase volumes from interpolated tagged slices and then matches them using an image registration framework. The approach uses several concepts from diffeomorphic image registration with a key novelty that defines a symmetric similarity metric on harmonic phase volumes from multiple orientations. The material property of harmonic phase solves the aperture problem of optical flow and intensity-based methods and is robust to tag fading. A harmonic magnitude volume is used in enforcing incompressibility in the tissue regions. The estimated motion fields are dense, incompressible, diffeomorphic, and inverse-consistent at a 3D voxel level. The method was evaluated using simulated phantoms, human brain data in mild head accelerations, human tongue data during speech, and an open cardiac dataset. The method shows comparable accuracy to three existing methods while demonstrating low computation time and robustness to tag fading and noise.

Magnetic resonance imaging (MRI) is crucial for in vivo detection and characterization of white matter lesions (WML) in multiple sclerosis. While these lesions have been studied for over two decades using MRI technology, automated segmentation remains challenging. Although the majority of statistical techniques for the automated segmentation of WML are based on a single imaging modality, recent advances have used multimodal techniques for identifying WML. Complementary imaging modalities emphasize different tissue properties, which can help identify and characterize interrelated features of lesions. However, prior work has ignored relationships between imaging modalities, which may be informative in this clinical context. To harness the coherent changes in these measurements, we utilized inter-modal coupling regression (IMCo) to estimate the covariance structure across modalities. We then used a local logistic regression, MIMoSA, which leverages new covariance features from IMCo regression as well as the mean structure of each imaging modality in order to model the probability that any voxel is part of a lesion. Finally, we introduced a novel thresholding algorithm to fully automate the estimation of the probability maps to generate fully automated segmentations masks for 94 subjects. To evaluate the performance of the automated segmentations generated using MIMoSA we compared results with gold standard manual segmentations. We demonstrate the superiority of MIMoSA to other automated segmentation techniques by comparing it to the OASIS algorithm as well as LesionTOADS. MIMoSA resulted in statistically significant improvement in lesion segmentation.

The simulation of magnetic resonance (MR) images plays an important role in the validation of image analysis algorithms such as image segmentation, due to lack of sufficient ground truth in real MR images. Previous work on MRI simulation has focused on explicitly modeling the MR image formation process. However, because of the overwhelming complexity of MR acquisition these simulations must involve simplifications and approximations that can result in visually unrealistic simulated images. In this work, we describe an example-based simulation framework, which uses an “atlas” consisting of an MR image and its anatomical models derived from the hard segmentation. The relationships between the MR image intensities and its anatomical models are learned using a patch-based regression that implicitly models the physics of the MR image formation. Given the anatomical models of a new brain, a new MR image can be simulated using the learned regression. This approach has been extended to also simulate intensity inhomogeneity artifacts based on the statistical model of training data. Results show that the example based MRI simulation method is capable of simulating different image contrasts and is robust to different choices of atlas. The simulated images resemble real MR images more than simulations produced by a physics-based model.

Different magnetic resonance imaging pulse sequences are used to generate image contrasts based on physical properties of tissues, which provide different and often complementary information about them. Therefore multiple image contrasts are useful for multimodal analysis of medical images. Often, medical image processing algorithms are optimized for particular image contrasts. If a desirable contrast is unavailable, contrast synthesis (or modality synthesis) methods try to “synthesize” the unavailable constrasts from the available ones. Most of the recent image synthesis methods generate synthetic brain images, while whole head magnetic resonance (MR) images can also be useful for many applications. We propose an atlas based patch matching algorithm to synthesize \(T_2-\)w whole head (including brain, skull, eyes etc.) images from \(T_1-\)w images for the purpose of distortion correction of diffusion weighted MR images. The geometric distortion in diffusion MR images due to inhomogeneous \(B_0\) magnetic field are often corrected by non-linearly registering the corresponding \(b=0\) image with zero diffusion gradient to an undistorted \(T_2-\)w image. We show that our synthetic \(T_2-\)w images can be used as a template in absence of a real \(T_2-\)w image. Our patch based method requires multiple atlases with \(T_1\) and \(T_2\) to be registered to a given target \(T_1\). Then for every patch on the target, multiple similar looking matching patches are found on the atlas \(T_1\) images and corresponding patches on the atlas \(T_2\) images are combined to generate a synthetic \(T_2\) of the target. We experimented on image data obtained from 44 patients with traumatic brain injury (TBI), and showed that our synthesized \(T_2\) images produce more accurate distortion correction than a state-of-the-art registration based image synthesis method.

PURPOSE: To explore: (1) the variability of upper cervical cord area (UCCA) measurements from volumetric brain 3D T1-weighted scans related to gradient nonlinearity (GNL) and subject positioning; (2) the effect of vendor-implemented GNL-corrections; and (3) easily applicable methods that can be used to retrospectively correct data. METHODS: A multiple sclerosis patient was scanned at 7 sites using 3T MRI scanners with the same 3D T1-weighted protocol without GNL-distortion correction. Two healthy subjects and a phantom were additionally scanned at a single site with varying table positions. 2D and 3D vendor-implemented GNL-correction algorithms and retrospective methods based on: (1) phantom data fit; (2) normalization with C2 vertebral body diameters; and (3) the Jacobian determinant of non-linear registrations to a template, were tested. RESULTS: Depending on positioning of the subject, GNL introduced up to 15% variability in UCCA measurements from volumetric brain T1-weighted scans when no distortion corrections were utilized. 3D vendor-implemented correction methods and the three proposed methods reduced this variability to less than 3%. CONCLUSION: Our results raise awareness of the significant impact that GNL can have on quantitative UCCA studies and point the way to prospectively and retrospectively managing GNL-distortions in a variety of settings, including clinical environments.

Noninvasive measurements of tissue deformation provide biomechanical insights of an organ, which can be used as clinical functional biomarkers or experimental data for validating computational simulations. However, acquisition of 3D displacement information is susceptible to experimental inconsistency and limited scan time. In this research, we describe the process of tracking tagged magnetic resonance imaging (MRI) as enforcing harmonic phase conservation in finite-element (FE) models. This concept is demonstrated as a tool for motion estimation in an experimental brain phantom, and images from the human heart and tongue. Our results demonstrate that the new methodology offers robustness to edge and large-displacement artifacts, and that it can be seamlessly coupled with numerical simulations for estimating fiber stretch in residually stressed tissue, or for inverse identification of muscle activation.

Mild traumatic brain injury (mTBI) and posttraumatic stress disorder (PTSD) may share common symptom and neuropsychological profiles in military service members (SMs) following deployment; while a connection between the two conditions is plausible, the relationship between them has been difficult to discern. The intent of this report is to enhance our understanding of the relationship between findings on structural and functional brain imaging and symptoms of PTSD. Within a cohort of SMs who did not meet criteria for PTSD but were willing to complete a comprehensive assessment within 2 months of their return from combat deployment, we conducted a nested case-control analysis comparing those with combat-related mTBI to age/gender-matched controls with diffusion tensor imaging, resting state functional magnetic resonance imaging and a range of psychological measures. We report degraded white matter integrity in those with a history of combat mTBI, and a positive correlation between the white matter microstructure and default mode network (DMN) connectivity. Higher clinician-administered and self-reported subthreshold PTSD symptoms were reported in those with combat mTBI. Our findings offer a potential mechanism through which mTBI may alter brain function, and in turn, contribute to PTSD symptoms.