Automated model-based tissue classification of MR images of the brain

Whole Brain Parcellation with Pathology: Validation on Ventriculomegaly Patients

Numerous brain disorders are associated with ventriculomegaly; normal pressure hydrocephalus (NPH) is one example. NPH presents with dementia-like symptoms and is often misdiagnosed as Alzheimer's due to its chronic nature and nonspecific presenting symptoms. However, unlike other forms of dementia NPH can be treated surgically with an over 80% success rate on appropriately selected patients. Accurate assessment of the ventricles, in particular its sub-compartments, is required to diagnose the condition. Existing segmentation algorithms fail to accurately identify the ventricles in patients with such extreme pathology. We present an improvement to a whole brain segmentation approach that accurately identifies the ventricles and parcellates them into four sub-compartments. Our work is a combination of patch-based tissue segmentation and multi-atlas registration-based labeling. We include a validation on NPH patients, demonstrating superior performance against state-of-the-art methods.

Pulmonary Lobe Segmentation With Probabilistic Segmentation of the Fissures and a Groupwise Fissure Prior

A fully automated, unsupervised lobe segmentation algorithm is presented based on a probabilistic segmentation of the fissures and the simultaneous construction of a populationmodel of the fissures. A two-class probabilistic segmentation segments the lung into candidate fissure voxels and the surrounding parenchyma. This was combined with anatomical information and a groupwise fissure prior to drive non-parametric surface fitting to obtain the final segmentation. The performance of our fissure segmentation was validated on 30 patients from the chronic obstructive pulmonary disease COPDGene cohort, achieving a high median F₁ -score of 0.90 and showed general insensitivity to filter parameters. We evaluated our lobe segmentation algorithm on the Lobe and Lung Analysis 2011 dataset, which contains 55 cases at varying levels of pathology. We achieved the highest score of 0.884 of the automated algorithms. Our method was further tested quantitatively and qualitatively on 80 patients from the COPDgene study at varying levels of functional impairment. Accurate segmentation of the lobes is shown at various degrees of fissure incompleteness for 96% of all cases. We also show the utility of including a groupwise prior in segmenting the lobes in regions of grossly incomplete fissures.

Partial volume model for brain MRI scan using MP2RAGE

MP2RAGE is a T1 weighted MRI sequence that estimates a composite image providing much reduction of the receiver bias, has a high intensity dynamic range, and provides an estimate of T1 mapping. It is, therefore, an appealing option for brain morphometry studies. However, previous studies have reported a difference in cortical thickness computed from MP2RAGE compared with widely used Multi-Echo MPRAGE. In this article, we demonstrated that using standard segmentation and partial volume estimation techniques on MP2RAGE introduces systematic errors, and we proposed a new model to estimate partial volume of the cortical gray matter. We also included in their model a local estimate of tissue intensity to take into account the natural variation of tissue intensity across the brain. A theoretical framework is provided and validated using synthetic and physical phantoms. A repeatability experiment comparing MPRAGE and MP2RAGE confirmed that MP2RAGE using our model could be considered for structural imaging in brain morphology study, with similar cortical thickness estimate than that computed with MPRAGE. Hum Brain Mapp 38:5115-5127, 2017. © 2017 Wiley Periodicals, Inc.

Automated segmentation of midbrain structures with high iron content

The substantia nigra (SN), the subthalamic nucleus (STN), and the red nucleus (RN) are midbrain structures of ample interest in many neuroimaging studies, which may benefit from the availability of automated segmentation methods. The high iron content of these structures awards them high contrast in quantitative susceptibility mapping (QSM) images. We present a novel segmentation method that leverages the information of these images to produce automated segmentations of the SN, STN, and RN. The algorithm builds a map of spatial priors for the structures by non-linearly registering a set of manually-traced training labels to the midbrain. The priors are used to inform a Gaussian mixture model of the image intensities, with smoothness constraints imposed to ensure anatomical plausibility. The method was validated on manual segmentations from a sample of 40 healthy younger and older subjects. Average Dice scores were 0.81 (0.05) for the SN, 0.66 (0.14) for the STN and 0.88 (0.04) for the RN in the left hemisphere, and similar values were obtained for the right hemisphere. In all structures, volumes of manual and automatically obtained segmentations were significantly correlated. The algorithm showed lower accuracy on R₂^* and T₂-weighted Fluid Attenuated Inversion Recovery (FLAIR) images, which are also sensitive to iron content. To illustrate an application of the method, we show that the automated segmentations were comparable to the manual ones regarding detection of age-related differences to putative iron content.

Effects of abstinence and chronic cigarette smoking on white matter microstructure in alcohol dependence: Diffusion tensor imaging at 4T

BACKGROUND

We previously reported widespread microstructural deficits of brain white matter in alcohol-dependent individuals (ALC) compared to light drinkers in a small 1.5T diffusion tensor imaging study employing tract-based spatial statistics. Using a larger dataset acquired at 4T, the present study is an extension that investigated the effects of alcohol consumption, abstinence from alcohol, and comorbid cigarette smoking on white matter microstructure.

METHODS

Tract-based spatial statistics were performed on 20 1-week-abstinent ALC, 52 1-month-abstinent ALC, and 30 controls. Regional measures of fractional anisotropy (FA) and mean diffusivity (MD) in the significant clusters were compared by Analysis of Covariance. The metrics were correlated with substance use history and behavioral measures.

RESULTS

1-week-abstinent ALC showed lower FA than controls in the corpus callosum, right cingulum, external capsule, and hippocampus. At 1 month of abstinence, only the FA in the body of the corpus callosum of ALC remained significantly different from controls. Some regional FA deficits correlated with more severe measures of drinking and smoking histories but only weakly with mood and impulsivity measures.

CONCLUSION

White matter microstructure is abnormal during early abstinence in alcohol dependent treatment seekers and recovers into the normal range within about four weeks. The compromised white matter was related to substance use severity, mood, and impulsivity. Our findings suggest that ALC may benefit from interventions that facilitate normalization of DTI metrics to maintain abstinence, via smoking cessation, cognitive-based therapy, and perhaps pharmacology to support remyelination.

Partial Discreteness: A Novel Prior for Magnetic Resonance Image Reconstruction

An important factor influencing the quality of magnetic resonance (MR) images is the reconstruction method that is employed, and specifically, the type of prior knowledge that is exploited during reconstruction. In this work, we introduce a new type of prior knowledge, partial discreteness (PD), where a small number of regions in the image are assumed to be homogeneous and can be well represented by a constant magnitude. In particular, we mathematically formalize the partial discreteness property based on a Gaussian Mixture Model (GMM) and derive a partial discreteness image representation that characterizes the salient features of partially discrete images: a constant intensity in homogeneous areas and texture in heterogeneous areas. The partial discreteness representation is then used to construct a novel prior dedicated to the reconstruction of partially discrete MR images. The strength of the proposed prior is demonstrated on various simulated and real k-space data-based experiments with partially discrete images. Results demonstrate that the PD algorithm performs competitively with state-of-the-art reconstruction methods, being flexible and easy to implement.

Brain structure in sagittal craniosynostosis

Craniosynostosis, the premature fusion of one or more cranial sutures, leads to grossly abnormal head shapes and pressure elevations within the brain caused by these deformities. To date, accepted treatments for craniosynostosis involve improving surgical skull shape aesthetics. However, the relationship between improved head shape and brain structure after surgery has not been yet established. Typically, clinical standard care involves the collection of diagnostic medical computed tomography (CT) imaging to evaluate the fused sutures and plan the surgical treatment. CT is known to provide very good reconstructions of the hard tissues in the skull but it fails to acquire good soft brain tissue contrast. This study intends to use magnetic resonance imaging to evaluate brain structure in a small dataset of sagittal craniosynostosis patients and thus quantify the effects of surgical intervention in overall brain structure. Very importantly, these effects are to be contrasted with normative shape, volume and brain structure databases. The work presented here wants to address gaps in clinical knowledge in craniosynostosis focusing on understanding the changes in brain volume and shape secondary to surgery, and compare those with normally developing children. This initial pilot study has the potential to add significant quality to the surgical care of a vulnerable patient population in whom we currently have limited understanding of brain developmental outcomes.

Machine learning in a graph framework for subcortical segmentation

Automated and reliable segmentation of subcortical structures from human brain magnetic resonance images is of great importance for volumetric and shape analyses in quantitative neuroimaging studies. However, poor boundary contrast and variable shape of these structures make the automated segmentation a tough task. We propose a 3D graph-based machine learning method, called LOGISMOS-RF, to segment the caudate and the putamen from brain MRI scans in a robust and accurate way. An atlas-based tissue classification and bias-field correction method is applied to the images to generate an initial segmentation for each structure. Then a 3D graph framework is utilized to construct a geometric graph for each initial segmentation. A locally trained random forest classifier is used to assign a cost to each graph node. The max-flow algorithm is applied to solve the segmentation problem. Evaluation was performed on a dataset of T1-weighted MRI's of 62 subjects, with 42 images used for training and 20 images for testing. For comparison, FreeSurfer and FSL approaches were also evaluated using the same dataset. Dice overlap coefficients and surface-to-surfaces distances between the automated segmentation and expert manual segmentations indicate the results of our method are statistically significantly more accurate than the other two methods, for both the caudate (Dice: 0.89 ± 0.03) and the putamen (0.89 ± 0.03).

A Rough Set Bounded Spatially Constrained Asymmetric Gaussian Mixture Model for Image Segmentation

Accurate image segmentation is an important issue in image processing, where Gaussian mixture models play an important part and have been proven effective. However, most Gaussian mixture model (GMM) based methods suffer from one or more limitations, such as limited noise robustness, over-smoothness for segmentations, and lack of flexibility to fit data. In order to address these issues, in this paper, we propose a rough set bounded asymmetric Gaussian mixture model with spatial constraint for image segmentation. First, based on our previous work where each cluster is characterized by three automatically determined rough-fuzzy regions, we partition the target image into three rough regions with two adaptively computed thresholds. Second, a new bounded indicator function is proposed to determine the bounded support regions of the observed data. The bounded indicator and posterior probability of a pixel that belongs to each sub-region is estimated with respect to the rough region where the pixel lies. Third, to further reduce over-smoothness for segmentations, two novel prior factors are proposed that incorporate the spatial information among neighborhood pixels, which are constructed based on the prior and posterior probabilities of the within- and between-clusters, and considers the spatial direction. We compare our algorithm to state-of-the-art segmentation approaches in both synthetic and real images to demonstrate the superior performance of the proposed algorithm.

Random forest regression for magnetic resonance image synthesis

By choosing different pulse sequences and their parameters, magnetic resonance imaging (MRI) can generate a large variety of tissue contrasts. This very flexibility, however, can yield inconsistencies with MRI acquisitions across datasets or scanning sessions that can in turn cause inconsistent automated image analysis. Although image synthesis of MR images has been shown to be helpful in addressing this problem, an inability to synthesize both T2-weighted brain images that include the skull and FLuid Attenuated Inversion Recovery (FLAIR) images has been reported. The method described herein, called REPLICA, addresses these limitations. REPLICA is a supervised random forest image synthesis approach that learns a nonlinear regression to predict intensities of alternate tissue contrasts given specific input tissue contrasts. Experimental results include direct image comparisons between synthetic and real images, results from image analysis tasks on both synthetic and real images, and comparison against other state-of-the-art image synthesis methods. REPLICA is computationally fast, and is shown to be comparable to other methods on tasks they are able to perform. Additionally REPLICA has the capability to synthesize both T2-weighted images of the full head and FLAIR images, and perform intensity standardization between different imaging datasets.

Two Time Point MS Lesion Segmentation in Brain MRI: An Expectation-Maximization Framework

Purpose: Lesion volume is a meaningful measure in multiple sclerosis (MS) prognosis. Manual lesion segmentation for computing volume in a single or multiple time points is time consuming and suffers from intra and inter-observer variability. Methods: In this paper, we present MSmetrix-long: a joint expectation-maximization (EM) framework for two time point white matter (WM) lesion segmentation. MSmetrix-long takes as input a 3D T1-weighted and a 3D FLAIR MR image and segments lesions in three steps: (1) cross-sectional lesion segmentation of the two time points; (2) creation of difference image, which is used to model the lesion evolution; (3) a joint EM lesion segmentation framework that uses output of step (1) and step (2) to provide the final lesion segmentation. The accuracy (Dice score) and reproducibility (absolute lesion volume difference) of MSmetrix-long is evaluated using two datasets. Results: On the first dataset, the median Dice score between MSmetrix-long and expert lesion segmentation was 0.63 and the Pearson correlation coefficient (PCC) was equal to 0.96. On the second dataset, the median absolute volume difference was 0.11 ml. Conclusions: MSmetrix-long is accurate and consistent in segmenting MS lesions. Also, MSmetrix-long compares favorably with the publicly available longitudinal MS lesion segmentation algorithm of Lesion Segmentation Toolbox.

Noise-Compensated, Bias-Corrected Diffusion Weighted Endorectal Magnetic Resonance Imaging via a Stochastically Fully-Connected Joint Conditional Random Field Model

Diffusion weighted magnetic resonance imaging (DW-MR) is a powerful tool in imaging-based prostate cancer screening and detection. Endorectal coils are commonly used in DW-MR imaging to improve the signal-to-noise ratio (SNR) of the acquisition, at the expense of significant intensity inhomogeneities (bias field) that worsens as we move away from the endorectal coil. The presence of bias field can have a significant negative impact on the accuracy of different image analysis tasks, as well as prostate tumor localization, thus leading to increased inter- and intra-observer variability. Retrospective bias correction approaches are introduced as a more efficient way of bias correction compared to the prospective methods such that they correct for both of the scanner and anatomy-related bias fields in MR imaging. Previously proposed retrospective bias field correction methods suffer from undesired noise amplification that can reduce the quality of bias-corrected DW-MR image. Here, we propose a unified data reconstruction approach that enables joint compensation of bias field as well as data noise in DW-MR imaging. The proposed noise-compensated, bias-corrected (NCBC) data reconstruction method takes advantage of a novel stochastically fully connected joint conditional random field (SFC-JCRF) model to mitigate the effects of data noise and bias field in the reconstructed MR data. The proposed NCBC reconstruction method was tested on synthetic DW-MR data, physical DW-phantom as well as real DW-MR data all acquired using endorectal MR coil. Both qualitative and quantitative analysis illustrated that the proposed NCBC method can achieve improved image quality when compared to other tested bias correction methods. As such, the proposed NCBC method may have potential as a useful retrospective approach for improving the consistency of image interpretations.

Efficient multi-scale 3D CNN with fully connected CRF for accurate brain lesion segmentation

We propose a dual pathway, 11-layers deep, three-dimensional Convolutional Neural Network for the challenging task of brain lesion segmentation. The devised architecture is the result of an in-depth analysis of the limitations of current networks proposed for similar applications. To overcome the computational burden of processing 3D medical scans, we have devised an efficient and effective dense training scheme which joins the processing of adjacent image patches into one pass through the network while automatically adapting to the inherent class imbalance present in the data. Further, we analyze the development of deeper, thus more discriminative 3D CNNs. In order to incorporate both local and larger contextual information, we employ a dual pathway architecture that processes the input images at multiple scales simultaneously. For post-processing of the network's soft segmentation, we use a 3D fully connected Conditional Random Field which effectively removes false positives. Our pipeline is extensively evaluated on three challenging tasks of lesion segmentation in multi-channel MRI patient data with traumatic brain injuries, brain tumours, and ischemic stroke. We improve on the state-of-the-art for all three applications, with top ranking performance on the public benchmarks BRATS 2015 and ISLES 2015. Our method is computationally efficient, which allows its adoption in a variety of research and clinical settings. The source code of our implementation is made publicly available.

Modeling 4D Pathological Changes by Leveraging Normative Models

With the increasing use of efficient multimodal 3D imaging, clinicians are able to access longitudinal imaging to stage pathological diseases, to monitor the efficacy of therapeutic interventions, or to assess and quantify rehabilitation efforts. Analysis of such four-dimensional (4D) image data presenting pathologies, including disappearing and newly appearing lesions, represents a significant challenge due to the presence of complex spatio-temporal changes. Image analysis methods for such 4D image data have to include not only a concept for joint segmentation of 3D datasets to account for inherent correlations of subject-specific repeated scans but also a mechanism to account for large deformations and the destruction and formation of lesions (e.g., edema, bleeding) due to underlying physiological processes associated with damage, intervention, and recovery. In this paper, we propose a novel framework that provides a joint segmentation-registration framework to tackle the inherent problem of image registration in the presence of objects not present in all images of the time series. Our methodology models 4D changes in pathological anatomy across time and and also provides an explicit mapping of a healthy normative template to a subject's image data with pathologies. Since atlas-moderated segmentation methods cannot explain appearance and locality pathological structures that are not represented in the template atlas, the new framework provides different options for initialization via a supervised learning approach, iterative semisupervised active learning, and also transfer learning, which results in a fully automatic 4D segmentation method. We demonstrate the effectiveness of our novel approach with synthetic experiments and a 4D multimodal MRI dataset of severe traumatic brain injury (TBI), including validation via comparison to expert segmentations. However, the proposed methodology is generic in regard to different clinical applications requiring quantitative analysis of 4D imaging representing spatio-temporal changes of pathologies.

Self Super-resolution for Magnetic Resonance Images

It is faster and therefore cheaper to acquire magnetic resonance images (MRI) with higher in-plane resolution than through-plane resolution. The low resolution of such acquisitions can be increased using post-processing techniques referred to as super-resolution (SR) algorithms. SR is known to be an ill-posed problem. Most state-of-the-art SR algorithms rely on the presence of external/training data to learn a transform that converts low resolution input to a higher resolution output. In this paper an SR approach is presented that is not dependent on any external training data and is only reliant on the acquired image. Patches extracted from the acquired image are used to estimate a set of new images, where each image has increased resolution along a particular direction. The final SR image is estimated by combining images in this set via the technique of Fourier Burst Accumulation. Our approach was validated on simulated low resolution MRI images, and showed significant improvement in image quality and segmentation accuracy when compared to competing SR methods. SR of FLuid Attenuated Inversion Recovery (FLAIR) images with lesions is also demonstrated.

Surface-driven registration method for the structure-informed segmentation of diffusion MR images

Current methods for processing diffusion MRI (dMRI) to map the connectivity of the human brain require precise delineations of anatomical structures. This requirement has been approached by either segmenting the data in native dMRI space or mapping the structural information from T1-weighted (T1w) images. The characteristic features of diffusion data in terms of signal-to-noise ratio, resolution, as well as the geometrical distortions caused by the inhomogeneity of magnetic susceptibility across tissues hinder both solutions. Unifying the two approaches, we propose regseg, a surface-to-volume nonlinear registration method that segments homogeneous regions within multivariate images by mapping a set of nested reference-surfaces. Accurate surfaces are extracted from a T1w image of the subject, using as target image the bivariate volume comprehending the fractional anisotropy (FA) and the apparent diffusion coefficient (ADC) maps derived from the dMRI dataset. We first verify the accuracy of regseg on a general context using digital phantoms distorted with synthetic and random deformations. Then we establish an evaluation framework using undistorted dMRI data from the Human Connectome Project (HCP) and realistic deformations derived from the inhomogeneity fieldmap corresponding to each subject. We analyze the performance of regseg computing the misregistration error of the surfaces estimated after being mapped with regseg onto 16 datasets from the HCP. The distribution of errors shows a 95% CI of 0.56-0.66mm, that is below the dMRI resolution (1.25mm, isotropic). Finally, we cross-compare the proposed tool against a nonlinear b0-to-T2w registration method, thereby obtaining a significantly lower misregistration error with regseg. The accurate mapping of structural information in dMRI space is fundamental to increase the reliability of network building in connectivity analyses, and to improve the performance of the emerging structure-informed techniques for dMRI data processing.

Fast and sequence-adaptive whole-brain segmentation using parametric Bayesian modeling

Quantitative analysis of magnetic resonance imaging (MRI) scans of the brain requires accurate automated segmentation of anatomical structures. A desirable feature for such segmentation methods is to be robust against changes in acquisition platform and imaging protocol. In this paper we validate the performance of a segmentation algorithm designed to meet these requirements, building upon generative parametric models previously used in tissue classification. The method is tested on four different datasets acquired with different scanners, field strengths and pulse sequences, demonstrating comparable accuracy to state-of-the-art methods on T1-weighted scans while being one to two orders of magnitude faster. The proposed algorithm is also shown to be robust against small training datasets, and readily handles images with different MRI contrast as well as multi-contrast data.

Surface-Based fMRI-Driven Diffusion Tractography in the Presence of Significant Brain Pathology: A Study Linking Structure and Function in Cerebral Palsy

Diffusion MRI (dMRI) tractography analyses are difficult to perform in the presence of brain pathology. Automated methods that rely on cortical parcellation for structural connectivity studies often fail, while manually defining regions is extremely time consuming and can introduce human error. Both methods also make assumptions about structure-function relationships that may not hold after cortical reorganisation. Seeding tractography with functional-MRI (fMRI) activation is an emerging method that reduces these confounds, but inherent smoothing of fMRI signal may result in the inclusion of irrelevant pathways. This paper describes a novel fMRI-seeded dMRI-analysis pipeline based on surface-meshes that reduces these issues and utilises machine-learning to generate task specific white matter pathways, minimising the requirement for manually-drawn ROIs. We directly compared this new strategy to a standard voxelwise fMRI-dMRI approach, by investigating correlations between clinical scores and dMRI metrics of thalamocortical and corticomotor tracts in 31 children with unilateral cerebral palsy. The surface-based approach successfully processed more participants (87%) than the voxel-based approach (65%), and provided significantly more-coherent tractography. Significant correlations between dMRI metrics and five clinical scores of function were found for the more superior regions of these tracts. These significant correlations were stronger and more frequently found with the surface-based method (15/20 investigated were significant; R² = 0.43–0.73) than the voxelwise analysis (2 sig. correlations; 0.38 & 0.49). More restricted fMRI signal, better-constrained tractography, and the novel track-classification method all appeared to contribute toward these differences.

Frontal Metabolite Concentration Deficits in Opiate Dependence Relate to Substance Use, Cognition, and Self-Regulation

OBJECTIVE

Proton magnetic resonance spectroscopy (1H MRS) in opiate dependence showed abnormalities in neuronal viability and glutamate concentration in the anterior cingulate cortex (ACC). Metabolite levels in dorsolateral prefrontal cortex (DLPFC) or orbitofrontal cortex (OFC) and their neuropsychological correlates have not been investigated in opiate dependence.

METHODS

Single-volume proton MRS at 4 Tesla and neuropsychological testing were conducted in 21 opiate-dependent individuals (OD) on buprenorphine maintenance therapy. Results were compared to 28 controls (CON) and 35 alcohol-dependent individuals (ALC), commonly investigated treatment-seekers providing context for OD evaluation. Metabolite concentrations were measured from ACC, DLPFC, OFC and parieto-occipital cortical (POC) regions.

RESULTS

Compared to CON, OD had lower concentrations of N-acetylaspartate (NAA), glutamate (Glu), creatine +phosphocreatine (Cr) and myo-Inositol (mI) in the DLPFC and lower NAA, Cr, and mI in the ACC. OD, ALC, and CON were equivalent on metabolite levels in the POC and γ-aminobutyric acid (GABA) concentration did not differ between groups in any region. In OD, prefrontal metabolite deficits in ACC Glu as well as DLPFC NAA and choline containing metabolites (Cho) correlated with poorer working memory, executive and visuospatial functioning; metabolite deficits in DLPFC Glu and ACC GABA and Cr correlated with substance use measures. In the OFC of OD, Glu and choline-containing metabolites were elevated and lower Cr concentration related to higher nonplanning impulsivity. Compared to 3 week abstinent ALC, OD had significant DLPFC metabolite deficits.

CONCLUSION

The anterior frontal metabolite profile of OD differed significantly from that of CON and ALC. The frontal lobe metabolite abnormalities in OD and their neuropsychological correlates may play a role in treatment outcome and could be explored as specific targets for improved OD treatment.

Segmentation of angiodysplasia lesions in WCE images using a MAP approach with Markov Random Fields

This paper deals with the segmentation of angiodysplasias in wireless capsule endoscopy images. These lesions are the cause of almost 10% of all gastrointestinal bleeding episodes, and its detection using the available software presents low sensitivity. This work proposes an automatic selection of a ROI using an image segmentation module based on the MAP approach where an accelerated version of the EM algorithm is used to iteratively estimate the model parameters. Spatial context is modeled in the prior probability density function using Markov Random Fields. The color space used was CIELab, specially the a component, which highlighted most these type of lesions. The proposed method is the first regarding this specific type of lesions, but when compared to other state-of-the-art segmentation methods, it almost doubles the results.

Bayesian longitudinal segmentation of hippocampal substructures in brain MRI using subject-specific atlases

The hippocampal formation is a complex, heterogeneous structure that consists of a number of distinct, interacting subregions. Atrophy of these subregions is implied in a variety of neurodegenerative diseases, most prominently in Alzheimer's disease (AD). Thanks to the increasing resolution of MR images and computational atlases, automatic segmentation of hippocampal subregions is becoming feasible in MRI scans. Here we introduce a generative model for dedicated longitudinal segmentation that relies on subject-specific atlases. The segmentations of the scans at the different time points are jointly computed using Bayesian inference. All time points are treated the same to avoid processing bias. We evaluate this approach using over 4700 scans from two publicly available datasets (ADNI and MIRIAD). In test-retest reliability experiments, the proposed method yielded significantly lower volume differences and significantly higher Dice overlaps than the cross-sectional approach for nearly every subregion (average across subregions: 4.5% vs. 6.5%, Dice overlap: 81.8% vs. 75.4%). The longitudinal algorithm also demonstrated increased sensitivity to group differences: in MIRIAD (69 subjects: 46 with AD and 23 controls), it found differences in atrophy rates between AD and controls that the cross sectional method could not detect in a number of subregions: right parasubiculum, left and right presubiculum, right subiculum, left dentate gyrus, left CA4, left HATA and right tail. In ADNI (836 subjects: 369 with AD, 215 with early cognitive impairment - eMCI - and 252 controls), all methods found significant differences between AD and controls, but the proposed longitudinal algorithm detected differences between controls and eMCI and differences between eMCI and AD that the cross sectional method could not find: left presubiculum, right subiculum, left and right parasubiculum, left and right HATA. Moreover, many of the differences that the cross-sectional method already found were detected with higher significance. The presented algorithm will be made available as part of the open-source neuroimaging package FreeSurfer.

Alterations in regional shape on ipsilateral and contralateral cortex contrast in children with unilateral cerebral palsy and are predictive of multiple outcomes

Congenital brain lesions result in a wide range of cerebral tissue alterations observed in children with cerebral palsy (CP) that are associated with a range of functional impairments. The relationship between injury severity and functional outcomes, however, remains poorly understood. This research investigates the differences in cortical shape between children with congenital brain lesions and typically developing children (TDC) and investigates the correlations between cortical shape and functional outcome in a large cohort of patients diagnosed with unilateral CP. Using 139 structural magnetic resonance images, including 95 patients with clinically diagnosed CP and 44 TDC, cortical segmentations were obtained using a modified expectation maximization algorithm. Three shape characteristics (cortical thickness, curvature, and sulcal depth) were computed within a number of cortical regions. Significant differences in these shape measures compared to the TDC were observed on both the injured hemisphere of children with CP (P < 0.004), as well as on the apparently uninjured hemisphere, illustrating potential compensatory mechanisms in these children. Furthermore, these shape measures were significantly correlated with several functional outcomes, including motor, cognition, vision, and communication (P < 0.012), with three out of these four models performing well on test set validation. This study highlights that cortical neuroplastic effects may be quantified using MR imaging, allowing morphological changes to be studied longitudinally, including any influence of treatment. Ultimately, such approaches could be used for the long term prediction of outcomes and the tailoring of treatment to individuals. Hum Brain Mapp 37:3588-3603, 2016. © 2016 Wiley Periodicals, Inc.

Parcellation of the Healthy Neonatal Brain into 107 Regions Using Atlas Propagation through Intermediate Time Points in Childhood

Neuroimage analysis pipelines rely on parcellated atlases generated from healthy individuals to provide anatomic context to structural and diffusion MRI data. Atlases constructed using adult data introduce bias into studies of early brain development. We aimed to create a neonatal brain atlas of healthy subjects that can be applied to multi-modal MRI data. Structural and diffusion 3T MRI scans were acquired soon after birth from 33 typically developing neonates born at term (mean postmenstrual age at birth 39⁺⁵ weeks, range 37⁺²–41⁺⁶). An adult brain atlas (SRI24/TZO) was propagated to the neonatal data using temporal registration via childhood templates with dense temporal samples (NIH Pediatric Database), with the final atlas (Edinburgh Neonatal Atlas, ENA33) constructed using the Symmetric Group Normalization (SyGN) method. After this step, the computed final transformations were applied to T2-weighted data, and fractional anisotropy, mean diffusivity, and tissue segmentations to provide a multi-modal atlas with 107 anatomical regions; a symmetric version was also created to facilitate studies of laterality. Volumes of each region of interest were measured to provide reference data from normal subjects. Because this atlas is generated from step-wise propagation of adult labels through intermediate time points in childhood, it may serve as a useful starting point for modeling brain growth during development.

BrainK for Structural Image Processing: Creating Electrical Models of the Human Head

BrainK is a set of automated procedures for characterizing the tissues of the human head from MRI, CT, and photogrammetry images. The tissue segmentation and cortical surface extraction support the primary goal of modeling the propagation of electrical currents through head tissues with a finite difference model (FDM) or finite element model (FEM) created from the BrainK geometries. The electrical head model is necessary for accurate source localization of dense array electroencephalographic (dEEG) measures from head surface electrodes. It is also necessary for accurate targeting of cerebral structures with transcranial current injection from those surface electrodes. BrainK must achieve five major tasks: image segmentation, registration of the MRI, CT, and sensor photogrammetry images, cortical surface reconstruction, dipole tessellation of the cortical surface, and Talairach transformation. We describe the approach to each task, and we compare the accuracies for the key tasks of tissue segmentation and cortical surface extraction in relation to existing research tools (FreeSurfer, FSL, SPM, and BrainVisa). BrainK achieves good accuracy with minimal or no user intervention, it deals well with poor quality MR images and tissue abnormalities, and it provides improved computational efficiency over existing research packages.

Consistent cortical reconstruction and multi-atlas brain segmentation

Whole brain segmentation and cortical surface reconstruction are two essential techniques for investigating the human brain. Spatial inconsistences, which can hinder further integrated analyses of brain structure, can result due to these two tasks typically being conducted independently of each other. FreeSurfer obtains self-consistent whole brain segmentations and cortical surfaces. It starts with subcortical segmentation, then carries out cortical surface reconstruction, and ends with cortical segmentation and labeling. However, this "segmentation to surface to parcellation" strategy has shown limitations in various cohorts such as older populations with large ventricles. In this work, we propose a novel "multi-atlas segmentation to surface" method called Multi-atlas CRUISE (MaCRUISE), which achieves self-consistent whole brain segmentations and cortical surfaces by combining multi-atlas segmentation with the cortical reconstruction method CRUISE. A modification called MaCRUISE(+) is designed to perform well when white matter lesions are present. Comparing to the benchmarks CRUISE and FreeSurfer, the surface accuracy of MaCRUISE and MaCRUISE(+) is validated using two independent datasets with expertly placed cortical landmarks. A third independent dataset with expertly delineated volumetric labels is employed to compare segmentation performance. Finally, 200MR volumetric images from an older adult sample are used to assess the robustness of MaCRUISE and FreeSurfer. The advantages of MaCRUISE are: (1) MaCRUISE constructs self-consistent voxelwise segmentations and cortical surfaces, while MaCRUISE(+) is robust to white matter pathology. (2) MaCRUISE achieves more accurate whole brain segmentations than independently conducting the multi-atlas segmentation. (3) MaCRUISE is comparable in accuracy to FreeSurfer (when FreeSurfer does not exhibit global failures) while achieving greater robustness across an older adult population. MaCRUISE has been made freely available in open source.

Development of cortical shape in the human brain from 6 to 24months of age via a novel measure of shape complexity

The quantification of local surface morphology in the human cortex is important for examining population differences as well as developmental changes in neurodegenerative or neurodevelopmental disorders. We propose a novel cortical shape measure, referred to as the 'shape complexity index' (SCI), that represents localized shape complexity as the difference between the observed distributions of local surface topology, as quantified by the shape index (SI) measure, to its best fitting simple topological model within a given neighborhood. We apply a relatively small, adaptive geodesic kernel to calculate the SCI. Due to the small size of the kernel, the proposed SCI measure captures fine differences of cortical shape. With this novel cortical feature, we aim to capture comparatively small local surface changes that capture a) the widening versus deepening of sulcal and gyral regions, as well as b) the emergence and development of secondary and tertiary sulci. Current cortical shape measures, such as the gyrification index (GI) or intrinsic curvature measures, investigate the cortical surface at a different scale and are less well suited to capture these particular cortical surface changes. In our experiments, the proposed SCI demonstrates higher complexity in the gyral/sulcal wall regions, lower complexity in wider gyral ridges and lowest complexity in wider sulcal fundus regions. In early postnatal brain development, our experiments show that SCI reveals a pattern of increased cortical shape complexity with age, as well as sexual dimorphisms in the insula, middle cingulate, parieto-occipital sulcal and Broca's regions. Overall, sex differences were greatest at 6months of age and were reduced at 24months, with the difference pattern switching from higher complexity in males at 6months to higher complexity in females at 24months. This is the first study of longitudinal, cortical complexity maturation and sex differences, in the early postnatal period from 6 to 24months of age with fine scale, cortical shape measures. These results provide information that complement previous studies of gyrification index in early brain development.

Chronic Cigarette Smoking in Healthy Middle-Aged Individuals Is Associated With Decreased Regional Brain N-acetylaspartate and Glutamate Levels

BACKGROUND

Cigarette smoking is associated with metabolite abnormalities in anterior brain regions, but it is unclear if these abnormalities are apparent in other regions. Additionally, relationships between regional brain metabolite levels and measures of decision making, risk taking, and impulsivity in smokers and nonsmokers have not been investigated.

METHODS

In young to middle-aged (predominately male) nonsmokers (n = 30) and smokers (n = 35), N-acetylaspartate (NAA), choline-containing compounds, creatine-containing compounds (Cr), myo-inositol (mI), and glutamate (Glu) levels in the anterior cingulate cortex and right dorsolateral prefrontal cortex (DLPFC) were compared via 4-tesla proton single volume magnetic resonance spectroscopy. Groups also were compared on NAA, choline-containing compounds, Cr, and mI concentrations in the gray matter and white matter of the four cerebral lobes and subcortical nuclei/regions with 1.5-tesla proton magnetic resonance spectroscopy. Associations of regional metabolite levels with neurocognitive, decision-making, risk-taking, and self-reported impulsivity measures were examined.

RESULTS

Smokers showed lower DLPFC NAA, Cr, mI and Glu concentrations and lower lenticular nuclei NAA level; smokers also demonstrated greater age-related decreases of DLPFC NAA and anterior cingulate cortex and DLPFC Glu levels. Smokers exhibited poorer decision making and greater impulsivity. Across the sample, higher NAA and Glu in the DLPFC and NAA concentrations in multiple lobar gray matter and white matter regions and subcortical nuclei were associated with better neurocognition and lower impulsivity.

CONCLUSIONS

This study provides additional novel evidence that chronic smoking in young and middle-aged individuals is associated with significant age-related neurobiological abnormalities in anterior frontal regions implicated in the development and maintenance of addictive disorders.

Differential diagnosis of neurodegenerative diseases using structural MRI data

Different neurodegenerative diseases can cause memory disorders and other cognitive impairments. The early detection and the stratification of patients according to the underlying disease are essential for an efficient approach to this healthcare challenge. This emphasizes the importance of differential diagnostics. Most studies compare patients and controls, or Alzheimer's disease with one other type of dementia. Such a bilateral comparison does not resemble clinical practice, where a clinician is faced with a number of different possible types of dementia.

Here we studied which features in structural magnetic resonance imaging (MRI) scans could best distinguish four types of dementia, Alzheimer's disease, frontotemporal dementia, vascular dementia, and dementia with Lewy bodies, and control subjects. We extracted an extensive set of features quantifying volumetric and morphometric characteristics from T1 images, and vascular characteristics from FLAIR images. Classification was performed using a multi-class classifier based on Disease State Index methodology. The classifier provided continuous probability indices for each disease to support clinical decision making.

A dataset of 504 individuals was used for evaluation. The cross-validated classification accuracy was 70.6% and balanced accuracy was 69.1% for the five disease groups using only automatically determined MRI features. Vascular dementia patients could be detected with high sensitivity (96%) using features from FLAIR images. Controls (sensitivity 82%) and Alzheimer's disease patients (sensitivity 74%) could be accurately classified using T1-based features, whereas the most difficult group was the dementia with Lewy bodies (sensitivity 32%). These results were notable better than the classification accuracies obtained with visual MRI ratings (accuracy 44.6%, balanced accuracy 51.6%). Different quantification methods provided complementary information, and consequently, the best results were obtained by utilizing several quantification methods.

The results prove that automatic quantification methods and computerized decision support methods are feasible for clinical practice and provide comprehensive information that may help clinicians in the diagnosis making.

•

Differential diagnostics of dementias was studied using structural MRI data.

•

504 patients with both T1 and FLAIR MRIs from five patient classes were evaluated.

•

Different fully automatic quantification methods were compared and combined.

•

Classification accuracy of 70.6% was obtained for 5-class classification problem.

•

Combination of several quantification methods was needed for optimal accuracy.

Patch-based augmentation of Expectation-Maximization for brain MRI tissue segmentation at arbitrary age after premature birth

Accurate automated tissue segmentation of premature neonatal magnetic resonance images is a crucial task for quantification of brain injury and its impact on early postnatal growth and later cognitive development. In such studies it is common for scans to be acquired shortly after birth or later during the hospital stay and therefore occur at arbitrary gestational ages during a period of rapid developmental change. It is important to be able to segment any of these scans with comparable accuracy. Previous work on brain tissue segmentation in premature neonates has focused on segmentation at specific ages. Here we look at solving the more general problem using adaptations of age specific atlas based methods and evaluate this using a unique manually traced database of high resolution images spanning 20 gestational weeks of development. We examine the complimentary strengths of age specific atlas-based Expectation-Maximization approaches and patch-based methods for this problem and explore the development of two new hybrid techniques, patch-based augmentation of Expectation-Maximization with weighted fusion and a spatial variability constrained patch search. The former approach seeks to combine the advantages of both atlas- and patch-based methods by learning from the performance of the two techniques across the brain anatomy at different developmental ages, while the latter technique aims to use anatomical variability maps learnt from atlas training data to locally constrain the patch-based search range. The proposed approaches were evaluated using leave-one-out cross-validation. Compared with the conventional age specific atlas-based segmentation and direct patch based segmentation, both new approaches demonstrate improved accuracy in the automated labeling of cortical gray matter, white matter, ventricles and sulcal cortical-spinal fluid regions, while maintaining comparable results in deep gray matter.

Altered white matter and cortical structure in neonates with antenatally diagnosed isolated ventriculomegaly

Ventriculomegaly (VM) is the most common central nervous system abnormality diagnosed antenatally, and is associated with developmental delay in childhood. We tested the hypothesis that antenatally diagnosed isolated VM represents a biological marker for altered white matter (WM) and cortical grey matter (GM) development in neonates.

25 controls and 21 neonates with antenatally diagnosed isolated VM had magnetic resonance imaging at 41.97(± 2.94) and 45.34(± 2.14) weeks respectively. T₂-weighted scans were segmented for volumetric analyses of the lateral ventricles, WM and cortical GM. Diffusion tensor imaging (DTI) measures were assessed using voxel-wise methods in WM and cortical GM; comparisons were made between cohorts.

Ventricular and cortical GM volumes were increased, and WM relative volume was reduced in the VM group. Regional decreases in fractional anisotropy (FA) and increases in mean diffusivity (MD) were demonstrated in WM of the VM group compared to controls. No differences in cortical DTI metrics were observed. At 2 years, neurodevelopmental delays, especially in language, were observed in 6/12 cases in the VM cohort.

WM alterations in isolated VM cases may be consistent with abnormal development of WM tracts involved in language and cognition. Alterations in WM FA and MD may represent neural correlates for later neurodevelopmental deficits.

•

This study compared brain development in neonates with isolated VM to controls.

•

Neonates with isolated VM have enlarged cortical volumes compared to controls.

•

FA was reduced and MD was increased in the WM of the VM cohort.

•

Children with antenatal isolated VM are at increased risk for language delay.

Volumetric and fiber-tracing MRI methods for gray and white matter

Magnetic resonance imaging (MRI) is capable of generating high-resolution brain images with fine anatomic detail and unique tissue contrasts that reveal structures that are not visible to the eye. Sharply defined gray- and white-matter interfaces allow for quantitative anatomic analysis that can be accurately performed with largely automated segmentation methods. In an analogous fashion, diffusion MRI in the brain provides structural information based on contrasts derived from the diffusivity of water in brain tissue, which can highlight the orientation of neuronal axons. Also using largely automated methods, diffusion MRI can be used to generate models of white-matter tracts throughout the brain, a method known as tractography, as well as characterize the microstructural integrity of neuronal axons. Tractographic analysis has helped to define connectivity in the brain that powerfully informs understanding of brain function, and, together with other diffusion metrics, is useful in evaluation of the normal and diseased brain. The quantitative methods of brain segmentation, tractography, and diffusion MRI extend MRI into a realm beyond visual inspection and provide otherwise unachievable sensitivity and specificity in the analysis of brain structure and function.

Early-Stage White Matter Lesions Detected by Multispectral MRI Segmentation Predict Progressive Cognitive Decline

White matter lesions (WML) are the main brain imaging surrogate of cerebral small-vessel disease. A new MRI tissue segmentation method, based on a discriminative clustering approach without explicit model-based added prior, detects partial WML volumes, likely representing very early-stage changes in normal-appearing brain tissue. This study investigated how the different stages of WML, from a “pre-visible” stage to fully developed lesions, predict future cognitive decline. MRI scans of 78 subjects, aged 65–84 years, from the Leukoaraiosis and Disability (LADIS) study were analyzed using a self-supervised multispectral segmentation algorithm to identify tissue types and partial WML volumes. Each lesion voxel was classified as having a small (33%), intermediate (66%), or high (100%) proportion of lesion tissue. The subjects were evaluated with detailed clinical and neuropsychological assessments at baseline and at three annual follow-up visits. We found that voxels with small partial WML predicted lower executive function compound scores at baseline, and steeper decline of executive scores in follow-up, independently of the demographics and the conventionally estimated hyperintensity volume on fluid-attenuated inversion recovery images. The intermediate and fully developed lesions were related to impairments in multiple cognitive domains including executive functions, processing speed, memory, and global cognitive function. In conclusion, early-stage partial WML, still too faint to be clearly detectable on conventional MRI, already predict executive dysfunction and progressive cognitive decline regardless of the conventionally evaluated WML load. These findings advance early recognition of small vessel disease and incipient vascular cognitive impairment.

A New MRI Masking Technique Based on Multi-Atlas Brain Segmentation in Controls and Schizophrenia: A Rapid and Viable Alternative to Manual Masking

Brain masking of MRI images separates brain from surrounding tissue and its accuracy is important for further imaging analyses. We implemented a new brain masking technique based on multi-atlas brain segmentation (MABS) and compared MABS to masks generated using FreeSurfer (FS; version 5.3), Brain Extraction Tool (BET), and Brainwash, using manually defined masks (MM) as the gold standard. We further determined the effect of different masking techniques on cortical and subcortical volumes generated by FreeSurfer.

METHODS

Images were acquired on a 3-Tesla MR Echospeed system General Electric scanner on five control and five schizophrenia subjects matched on age, sex, and IQ. Automated masks were generated from MABS, FS, BET, and Brainwash, and compared to MM using these metrics: a) volume difference from MM; b) Dice coefficients; and c) intraclass correlation coefficients.

RESULTS

Mean volume difference between MM and MABS masks was significantly less than the difference between MM and FS or BET masks. Dice coefficient between MM and MABS was significantly higher than Dice coefficients between MM and FS, BET, or Brainwash. For subcortical and left cortical regions, MABS volumes were closer to MM volumes than were BET or FS volumes. For right cortical regions, MABS volumes were closer to MM volumes than were BET volumes.

CONCLUSIONS

Brain masks generated using FreeSurfer, BET, and Brainwash are rapidly obtained, but are less accurate than manually defined masks. Masks generated using MABS, in contrast, resemble more closely the gold standard of manual masking, thereby offering a rapid and viable alternative.

Computer-aided detection of cerebral microbleeds in susceptibility-weighted imaging

Susceptibility-weighted imaging (SWI) is recognized as the preferred MRI technique for visualizing cerebral vasculature and related pathologies such as cerebral microbleeds (CMBs). Manual identification of CMBs is time-consuming, has limited reliability and reproducibility, and is prone to misinterpretation. In this paper, a novel computer-aided microbleed detection technique based on machine learning is presented: First, spherical-like objects (potential CMB candidates) with their corresponding bounding boxes were detected using a novel multi-scale Laplacian of Gaussian technique. A set of robust 3-dimensional Radon- and Hessian-based shape descriptors within each bounding box were then extracted to train a cascade of binary random forests (RF). The cascade consists of consecutive independent RF classifiers with low to high posterior probability constraints to handle imbalanced training sets (CMBs and non-CMBs), and to progressively improve detection rates. The proposed method was validated on 66 subjects whose CMBs were manually stratified into "possible" and "definite" by two medical experts. The proposed technique achieved a sensitivity of 87% and an average false detection rate of 27.1 CMBs per subject on the "possible and definite" set. A sensitivity of 93% and false detection rate of 10 CMBs per subject was also achieved on the "definite" set. The proposed automated approach outperforms state of the art methods, and promises to enhance manual expert screening. Benefits include improved reliability, minimization of intra-rater variability and a reduction in assessment time.

Regional growth and atlasing of the developing human brain

Detailed morphometric analysis of the neonatal brain is required to characterise brain development and define neuroimaging biomarkers related to impaired brain growth. Accurate automatic segmentation of neonatal brain MRI is a prerequisite to analyse large datasets. We have previously presented an accurate and robust automatic segmentation technique for parcellating the neonatal brain into multiple cortical and subcortical regions. In this study, we further extend our segmentation method to detect cortical sulci and provide a detailed delineation of the cortical ribbon. These detailed segmentations are used to build a 4-dimensional spatio-temporal structural atlas of the brain for 82 cortical and subcortical structures throughout this developmental period. We employ the algorithm to segment an extensive database of 420 MR images of the developing brain, from 27 to 45 weeks post-menstrual age at imaging. Regional volumetric and cortical surface measurements are derived and used to investigate brain growth and development during this critical period and to assess the impact of immaturity at birth. Whole brain volume, the absolute volume of all structures studied, cortical curvature and cortical surface area increased with increasing age at scan. Relative volumes of cortical grey matter, cerebellum and cerebrospinal fluid increased with age at scan, while relative volumes of white matter, ventricles, brainstem and basal ganglia and thalami decreased. Preterm infants at term had smaller whole brain volumes, reduced regional white matter and cortical and subcortical grey matter volumes, and reduced cortical surface area compared with term born controls, while ventricular volume was greater in the preterm group. Increasing prematurity at birth was associated with a reduction in total and regional white matter, cortical and subcortical grey matter volume, an increase in ventricular volume, and reduced cortical surface area.

•

A novel methodology is proposed for delineating the cortical ribbon.

•

Regional brain growth is assessed in the developing preterm brain.

•

We investigate the effect of prematurity on brain growth and cortical development.

•

A spatio-temporal neonatal atlas is constructed with 82 brain structures.

SEGMENTATION OF BRAIN MRIs BY SUPPORT VECTOR MACHINE: DETECTION AND CHARACTERIZATION OF STROKES

The aim of our work is the segmentation of healthy and pathological brains to obtain brain structures and extract strokes. We used real magnetic resonance (MR) images weighted on diffusion. The brain was isolated, and the images were filtered by an anisotropic filter, and then segmented by support vector machines (SVMs). We first applied the method on synthetic images to test the performance of the algorithm and adjust the parameters. Then, we compared our results with those obtained by a cooperative approach proposed in a previous paper.

Geodesic Information Flows: Spatially-Variant Graphs and Their Application to Segmentation and Fusion

Clinical annotations, such as voxel-wise binary or probabilistic tissue segmentations, structural parcellations, pathological regions-of-interest and anatomical landmarks are key to many clinical studies. However, due to the time consuming nature of manually generating these annotations, they tend to be scarce and limited to small subsets of data. This work explores a novel framework to propagate voxel-wise annotations between morphologically dissimilar images by diffusing and mapping the available examples through intermediate steps. A spatially-variant graph structure connecting morphologically similar subjects is introduced over a database of images, enabling the gradual diffusion of information to all the subjects, even in the presence of large-scale morphological variability. We illustrate the utility of the proposed framework on two example applications: brain parcellation using categorical labels and tissue segmentation using probabilistic features. The application of the proposed method to categorical label fusion showed highly statistically significant improvements when compared to state-of-the-art methodologies. Significant improvements were also observed when applying the proposed framework to probabilistic tissue segmentation of both synthetic and real data, mainly in the presence of large morphological variability.

Subject-Specific Sparse Dictionary Learning for Atlas-Based Brain MRI Segmentation

Quantitative measurements from segmentations of human brain magnetic resonance (MR) images provide important biomarkers for normal aging and disease progression. In this paper, we propose a patch-based tissue classification method from MR images that uses a sparse dictionary learning approach and atlas priors. Training data for the method consists of an atlas MR image, prior information maps depicting where different tissues are expected to be located, and a hard segmentation. Unlike most atlas-based classification methods that require deformable registration of the atlas priors to the subject, only affine registration is required between the subject and training atlas. A subject-specific patch dictionary is created by learning relevant patches from the atlas. Then the subject patches are modeled as sparse combinations of learned atlas patches leading to tissue memberships at each voxel. The combination of prior information in an example-based framework enables us to distinguish tissues having similar intensities but different spatial locations. We demonstrate the efficacy of the approach on the application of whole-brain tissue segmentation in subjects with healthy anatomy and normal pressure hydrocephalus, as well as lesion segmentation in multiple sclerosis patients. For each application, quantitative comparisons are made against publicly available state-of-the art approaches.

Serial longitudinal magnetic resonance imaging data indicate non-linear regional gray matter volume recovery in abstinent alcohol-dependent individuals

The trajectory of regional volume changes during the first year of sustained abstinence in those recovering from an alcohol use disorder is unclear because previous research typically employed only two assessment points. To better understand the trajectory of regional brain volume recovery in treatment-seeking alcohol-dependent individuals (ALC), regional brain volumes were measured after 1 week, 1 month and 7.5 months of sustained abstinence via magnetic resonance imaging at 1.5 T. ALC showed significant volume increases in frontal, parietal and occipital gray matter (GM) and white matter (WM), total cortical GM and total lobar WM, thalamus and cerebellum, and decreased ventricular volume over 7.5 months of abstinence. Volume increases in regional GM were significantly greater over 1 week to 1 month than from 1 month to 7.5 months of abstinence, indicating a non-linear rate of change in regional GM over 7.5 months. Overall, regional lobar WM showed linear volume increases over 7.5 months. With increasing age, smoking ALC showed lower frontal and total cortical GM volume recovery than non-smoking ALC. Despite significant volume increases, ALC showed smaller GM volumes in all regions, except the frontal cortex, than controls after 7.5 months of abstinence. ALC and controls showed no regional WM volume differences at any assessment point. In non-smoking ALC only, increasing regional GM and WM volumes were related to improving processing speed. Findings may indicate a differential rate of recovery of cell types/cellular components contributing to GM and WM volume during early abstinence, and that GM volume deficits persist after 7.5 months of sustained sobriety in this ALC cohort.