Topology-preserving tissue classification of magnetic resonance brain images

Age-related differences in human cortical microstructure depend on the distance to the nearest vein

Age-related differences in cortical microstructure are used to understand the neuronal mechanisms that underlie human brain ageing. The cerebral vasculature contributes to cortical ageing, but its precise interaction with cortical microstructure is poorly understood. In a cross-sectional study, we combine venous imaging with vessel distance mapping to investigate the interaction between venous distances and age-related differences in the microstructural architecture of the primary somatosensory cortex, the primary motor cortex and additional areas in the frontal cortex as non-sensorimotor control regions. We scanned 18 younger adults and 17 older adults using 7 Tesla MRI to measure age-related changes in longitudinal relaxation time (T1) and quantitative susceptibility mapping (QSM) values at 0.5 mm isotropic resolution. We modelled different cortical depths using an equi-volume approach and assessed the distance of each voxel to its nearest vein using vessel distance mapping. Our data reveal a dependence of cortical quantitative T1 values and positive QSM values on venous distance. In addition, there is an interaction between venous distance and age on quantitative T1 values, driven by lower quantitative T1 values in older compared to younger adults in voxels that are closer to a vein. Together, our data show that the local venous architecture explains a significant amount of variance in standard measures of cortical microstructure and should be considered in neurobiological models of human brain organisation and cortical ageing.

Comparative analysis of alignment algorithms for macular optical coherence tomography imaging

BACKGROUND

Optical coherence tomography (OCT) is the most important and commonly utilized imaging modality in ophthalmology and is especially crucial for the diagnosis and management of macular diseases. Each OCT volume is typically only available as a series of cross-sectional images (B-scans) that are accessible through proprietary software programs which accompany the OCT machines. To maximize the potential of OCT imaging for machine learning purposes, each OCT image should be analyzed en bloc as a 3D volume, which requires aligning all the cross-sectional images within a particular volume.

METHODS

A dataset of OCT B-scans obtained from 48 age-related macular degeneration (AMD) patients and 50 normal controls was used to evaluate five registration algorithms. After alignment of B-scans from each patient, an en face surface map was created to measure the registration quality, based on an automatically generated Laplace difference of the surface map-the smoother the surface map, the smaller the average Laplace difference. To demonstrate the usefulness of B-scan alignment, we trained a 3D convolutional neural network (CNN) to detect age-related macular degeneration (AMD) on OCT images and compared the performance of the model with and without B-scan alignment.

RESULTS

The mean Laplace difference of the surface map before registration was 27 ± 4.2 pixels for the AMD group and 26.6 ± 4 pixels for the control group. After alignment, the smoothness of the surface map was improved, with a mean Laplace difference of 5.5 ± 2.7 pixels for Advanced Normalization Tools Symmetric image Normalization (ANTs-SyN) registration algorithm in the AMD group and a mean Laplace difference of 4.3 ± 1.4.2 pixels for ANTs in the control group. Our 3D CNN achieved superior performance in detecting AMD, when aligned OCT B-scans were used (AUC 0.95 aligned vs. 0.89 unaligned).

CONCLUSIONS

We introduced a novel metric to quantify OCT B-scan alignment and compared the effectiveness of five alignment algorithms. We confirmed that alignment could be improved in a statistically significant manner with readily available alignment algorithms that are available to the public, and the ANTs algorithm provided the most robust performance overall. We further demonstrated that alignment of OCT B-scans will likely be useful for training 3D CNN models.

The 3D Structural Architecture of the Human Hand Area Is Nontopographic

The functional topography of the human primary somatosensory cortex hand area is a widely studied model system to understand sensory organization and plasticity. It is so far unclear whether the underlying 3D structural architecture also shows a topographic organization. We used 7 Tesla (7T) magnetic resonance imaging (MRI) data to quantify layer-specific myelin, iron, and mineralization in relation to population receptive field maps of individual finger representations in Brodman area 3b (BA 3b) of human S1 in female and male younger adults. This 3D description allowed us to identify a characteristic profile of layer-specific myelin and iron deposition in the BA 3b hand area, but revealed an absence of structural differences, an absence of low-myelin borders, and high similarity of 3D microstructure profiles between individual fingers. However, structural differences and borders were detected between the hand and face areas. We conclude that the 3D structural architecture of the human hand area is nontopographic, unlike in some monkey species, which suggests a high degree of flexibility for functional finger organization and a new perspective on human topographic plasticity.SIGNIFICANCE STATEMENT Using ultra-high-field MRI, we provide the first comprehensive in vivo description of the 3D structural architecture of the human BA 3b hand area in relation to functional population receptive field maps. High similarity of precise finger-specific 3D profiles, together with an absence of structural differences and an absence of low-myelin borders between individual fingers, reveals the 3D structural architecture of the human hand area to be nontopographic. This suggests reduced structural limitations to cortical plasticity and reorganization and allows for shared representational features across fingers.

Coordinate Translator for Learning Deformable Medical Image Registration

The majority of deep learning (DL) based deformable image registration methods use convolutional neural networks (CNNs) to estimate displacement fields from pairs of moving and fixed images. This, however, requires the convolutional kernels in the CNN to not only extract intensity features from the inputs but also understand image coordinate systems. We argue that the latter task is challenging for traditional CNNs, limiting their performance in registration tasks. To tackle this problem, we first introduce Coordinate Translator, a differentiable module that identifies matched features between the fixed and moving image and outputs their coordinate correspondences without the need for training. It unloads the burden of understanding image coordinate systems for CNNs, allowing them to focus on feature extraction. We then propose a novel deformable registration network, im2grid, that uses multiple Coordinate Translator's with the hierarchical features extracted from a CNN encoder and outputs a deformation field in a coarse-to-fine fashion. We compared im2grid with the state-of-the-art DL and non-DL methods for unsupervised 3D magnetic resonance image registration. Our experiments show that im2grid outperforms these methods both qualitatively and quantitatively.

Initial evidence of abnormal brain plasticity in anorexia nervosa: an ultra-high field study

Anorexia Nervosa has been associated with white matter abnormalities implicating subcortical abnormal myelination. Extending these findings to intracortical myelin has been challenging but ultra-high field neuroimaging offers new methodological opportunities. To test the integrity of intracortical myelin in AN we used 7 T neuroimaging to acquire T1-weighted images optimized for intracortical myelin from seven females with AN (age range: 18-33) and 11 healthy females (age range: 23-32). Intracortical T₁ values (inverse index of myelin concentration) were extracted from 148 cortical regions at ten depth-levels across the cortical ribbon. Across all cortical regions, these levels were averaged to generate estimates of total intracortical myelin concentration and were clustered using principal component analyses into two clusters; the outer cluster comprised T₁ values across depth-levels ranging from the CSF boundary to the middle of the cortical regions and the inner cluster comprised T₁ values across depth-levels ranging from the middle of the cortical regions to the gray/white matter boundary. Individuals with AN exhibited higher T₁ values (i.e., decreased intracortical myelin concentration) in all three metrics. It remains to be established if these abnormalities result from undernutrition or specific lipid nutritional imbalances, or are trait markers; and whether they may contribute to neurobiological deficits seen in AN.

BISON: Brain tissue segmentation pipeline using T1 -weighted magnetic resonance images and a random forest classifier

PURPOSE

Tissue segmentation from T₁ -weighted (T1W) MRI is a critical requirement in many neuroscience and clinical applications. However, accurate tissue segmentation is challenging because of the variabilities in tissue intensity profiles caused by differences in scanner models, acquisition protocols, and age. In addition, many methods assume healthy anatomy and fail in the presence of pathology such as white matter hyperintensities (WMHs). We present BISON (Brain tISsue segmentatiON), a new pipeline for tissue segmentation using a random forest classifier and a set of intensity and location priors based on T1W MRI.

METHODS

BISON was developed and cross-validated using multiscanner manual labels of 72 subjects aged 5 to 96 years. We also assessed the test-retest reliability of BISON on two data sets: 20 subjects with scan/rescan MR images and manual segmentations and 90 scans from a single individual. The results were compared against Atropos, a state-of-the-art commonly used tissue classification method from advanced normalization tools (ANTs).

RESULTS

BISON cross-validation dice kappa values against manual segmentations of 72 MRI volumes yielded κ_GM = 0.88, κ_WM = 0.85, κ_CSF = 0.77, outperforming Atropos (κ_GM = 0.79, κ_WM = 0.84, κ_CSF = 0.64), test-retest values on 20 subjects of κ_GM = 0.94, κ_WM = 0.92, κ_CSF = 0.77 outperforming both manual (κ_GM = 0.92, κ_WM = 0.91, κ_CSF =0.74) and Atropos (κ_GM = 0.87, κ_WM = 0.92, κ_CSF = 0.79). Finally, BISON outperformed Atropos, FAST (fast automated segmentation tool) from the FMRIB (Functional Magnetic Resonance Imaging of the Brain) Software Library, and SPM12 (statistical parametric mapping 12) in the presence of WMHs.

CONCLUSION

BISON can provide accurate and robust segmentations in data from various age ranges and scanner models, making it ideal for performing tissue classification in large multicenter and multiscanner databases.

Automatic segmentation of brain MRI using a novel patch-wise U-net deep architecture

Accurate segmentation of brain magnetic resonance imaging (MRI) is an essential step in quantifying the changes in brain structure. Deep learning in recent years has been extensively used for brain image segmentation with highly promising performance. In particular, the U-net architecture has been widely used for segmentation in various biomedical related fields. In this paper, we propose a patch-wise U-net architecture for the automatic segmentation of brain structures in structural MRI. In the proposed brain segmentation method, the non-overlapping patch-wise U-net is used to overcome the drawbacks of conventional U-net with more retention of local information. In our proposed method, the slices from an MRI scan are divided into non-overlapping patches that are fed into the U-net model along with their corresponding patches of ground truth so as to train the network. The experimental results show that the proposed patch-wise U-net model achieves a Dice similarity coefficient (DSC) score of 0.93 in average and outperforms the conventional U-net and the SegNet-based methods by 3% and 10%, respectively, for on Open Access Series of Imaging Studies (OASIS) and Internet Brain Segmentation Repository (IBSR) dataset.

Triple visual hemifield maps in a case of optic chiasm hypoplasia

In humans, each hemisphere comprises an overlay of two visuotopic maps of the contralateral visual field, one from each eye. Is the capacity of the visual cortex limited to these two maps or are plastic mechanisms available to host more maps? We determined the cortical organization of the visual field maps in a rare individual with chiasma hypoplasia, where visual cortex plasticity is challenged to accommodate three hemifield maps. Using high-resolution fMRI at 7T and diffusion-weighted MRI at 3T, we found three hemiretinal inputs, instead of the normal two, to converge onto the left hemisphere. fMRI-based population receptive field mapping of the left V1-V3 at 3T revealed three superimposed hemifield representations in the left visual cortex, i.e. two representations of opposing visual hemifields from the left eye and one right hemifield representation from the right eye. We conclude that developmental plasticity including the re-wiring of local intra- and cortico-cortical connections is pivotal to support the coexistence and functioning of three hemifield maps within one hemisphere.

Differential Changes in Functional Connectivity of Striatum-Prefrontal and Striatum-Motor Circuits in Premanifest Huntington's Disease

BACKGROUND

Huntington's disease (HD) is a progressive neurodegenerative disorder. The striatum is one of the first brain regions that show detectable atrophy in HD. Previous studies using functional magnetic resonance imaging (fMRI) at 3 tesla (3 T) revealed reduced functional connectivity between striatum and motor cortex in the prodromal period of HD. Neuroanatomical and neurophysiological studies have suggested segregated corticostriatal pathways with distinct loops involving different cortical regions, which may be investigated using fMRI at an ultra-high field (7 T) with enhanced sensitivity compared to lower fields.

OBJECTIVES

We performed fMRI at 7 T to assess functional connectivity between the striatum and several chosen cortical areas including the motor and prefrontal cortex, in order to better understand brain changes in the striatum-cortical pathways.

METHOD

13 manifest subjects (age 51 ± 13 years, cytosine-adenine-guanine [CAG] repeat 45 ± 5, Unified Huntington's Disease Rating Scale [UHDRS] motor score 32 ± 17), 8 subjects in the close-to-onset premanifest period (age 38 ± 10 years, CAG repeat 44 ± 2, UHDRS motor score 8 ± 2), 11 subjects in the far-from-onset premanifest period (age 38 ± 11 years, CAG repeat 42 ± 2, UHDRS motor score 1 ± 2), and 16 healthy controls (age 44 ± 15 years) were studied. The functional connectivity between the striatum and several cortical areas was measured by resting state fMRI at 7 T and analyzed in all participants.

RESULTS

Compared to controls, functional connectivity between striatum and premotor area, supplementary motor area, inferior frontal as well as middle frontal regions was altered in HD (all p values <0.001). Specifically, decreased striatum-motor connectivity but increased striatum-prefrontal connectivity were found in premanifest HD subjects. Altered functional connectivity correlated consistently with genetic burden, but not with clinical scores.

CONCLUSIONS

Differential changes in functional connectivity of striatum-prefrontal and striatum-motor circuits can be found in early and premanifest HD. This may imply a compensatory mechanism, where additional cortical regions are recruited to subserve functions that have been impaired due to HD pathology. Our results suggest the potential value of functional connectivity as a marker for future clinical trials in HD.

Brain and retinal atrophy in African-Americans versus Caucasian-Americans with multiple sclerosis: a longitudinal study

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.

Iterative spatial fuzzy clustering for 3D brain magnetic resonance image supervoxel segmentation

BACKGROUND

Although supervoxel segmentation methods have been employed for brain Magnetic Resonance Image (MRI) processing and analysis, due to the specific features of the brain, including complex-shaped internal structures and partial volume effect, their performance remains unsatisfactory.

NEW METHODS

To address these issues, this paper presents a novel iterative spatial fuzzy clustering (ISFC) algorithm to generate 3D supervoxels for brain MRI volume based on prior knowledge. This work makes use of the common topology among the human brains to obtain a set of seed templates from a population-based brain template MRI image. After selecting the number of supervoxels, the corresponding seed template is projected onto the considered individual brain for generating reliable seeds. Then, to deal with the influence of partial volume effect, an efficient iterative spatial fuzzy clustering algorithm is proposed to allocate voxels to the seeds and to generate the supervoxels for the overall brain MRI volume.

RESULTS

The performance of the proposed algorithm is evaluated on two widely used public brain MRI datasets and compared with three other up-to-date methods.

CONCLUSIONS

The proposed algorithm can be utilized for several brain MRI processing and analysis, including tissue segmentation, tumor detection and segmentation, functional parcellation and registration.

Depth-dependent intracortical myelin organization in the living human brain determined by in vivo ultra-high field magnetic resonance imaging

BACKGROUND

Intracortical myelin is a key determinant of neuronal synchrony and plasticity that underpin optimal brain function. Magnetic resonance imaging (MRI) facilitates the examination of intracortical myelin but presents with methodological challenges. Here we describe a whole-brain approach for the in vivo investigation of intracortical myelin in the human brain using ultra-high field MRI.

METHODS

Twenty-five healthy adults were imaged in a 7 Tesla MRI scanner using diffusion-weighted imaging and a T₁-weighted sequence optimized for intracortical myelin contrast. Using an automated pipeline, T₁ values were extracted at 20 depth-levels from each of 148 cortical regions. In each cortical region, T₁ values were used to infer myelin concentration and to construct a non-linearity index as a measure the spatial distribution of myelin across the cortical ribbon. The relationship of myelin concentration and the non-linearity index with other neuroanatomical properties were investigated. Five patients with multiple sclerosis were also assessed using the same protocol as positive controls.

RESULTS

Intracortical T₁ values decreased between the outer brain surface and the gray-white matter boundary following a slope that showed a slight leveling between 50% and 75% of cortical depth. Higher-order regions in the prefrontal, cingulate and insular cortices, displayed higher non-linearity indices than sensorimotor regions. Across all regions, there was a positive association between T₁ values and non-linearity indices (P < 10^-5). Both T₁ values (P < 10^-5) and non-linearity indices (P < 10^-15) were associated with cortical thickness. Higher myelin concentration but only in the deepest cortical levels was associated with increased subcortical fractional anisotropy (P = 0.05).

CONCLUSIONS

We demonstrate the usefulness of an automatic, whole-brain method to perform depth-dependent examination of intracortical myelin organization. The extracted metrics, T₁ values and the non-linearity index, have characteristic patterns across cortical regions, and are associated with thickness and underlying white matter microstructure.

Urgent challenges in quantification and interpretation of brain grey matter atrophy in individual MS patients using MRI

Atrophy of the brain grey matter (GM) is an accepted and important feature of multiple sclerosis (MS). However, its accurate measurement is hampered by various technical, pathological and physiological factors. As a consequence, it is challenging to investigate the role of GM atrophy in the disease process as well as the effect of treatments that aim to reduce neurodegeneration. In this paper we discuss the most important challenges currently hampering the measurement and interpretation of GM atrophy in MS. The focus is on measurements that are obtained in individual patients rather than on group analysis methods, because of their importance in clinical trials and ultimately in clinical care. We discuss the sources and possible solutions of the current challenges, and provide recommendations to achieve reliable measurement and interpretation of brain GM atrophy in MS.

Multi-atlas segmentation of the hydrocephalus brain using an adaptive ventricle atlas

Normal pressure hydrocephalus (NPH) is a brain disorder caused by disruption of the flow of cerebrospinal fluid (CSF). The dementia-like symptoms of NPH are often mistakenly attributed to Alzheimer's disease. However, if correctly diagnosed, NPH patients can potentially be treated and their symptoms reversed through surgery. Observing the dilated ventricles through magnetic resonance imaging (MRI) is one element in diagnosing NPH. Diagnostic accuracy therefore benefits from accurate, automatic parcellation of the ventricular system into its sub-compartments. We present an improvement to a whole brain segmentation approach designed for subjects with enlarged and deformed ventricles. Our method incorporates an adaptive ventricle atlas from an NPH-atlas-based segmentation as a prior and uses a more robust relaxation scheme for the multi-atlas label fusion approach that accurately labels the four sub-compartments of the ventricular system. We validated our method on NPH patients, demonstrating improvement over state-of-the-art segmentation techniques.

Cross contrast multi-channel image registration using image synthesis for MR brain images

Multi-modal deformable registration is important for many medical image analysis tasks such as atlas alignment, image fusion, and distortion correction. Whereas a conventional method would register images with different modalities using modality independent features or information theoretic metrics such as mutual information, this paper presents a new framework that addresses the problem using a two-channel registration algorithm capable of using mono-modal similarity measures such as sum of squared differences or cross-correlation. To make it possible to use these same-modality measures, image synthesis is used to create proxy images for the opposite modality as well as intensity-normalized images from each of the two available images. The new deformable registration framework was evaluated by performing intra-subject deformation recovery, intra-subject boundary alignment, and inter-subject label transfer experiments using multi-contrast magnetic resonance brain imaging data. Three different multi-channel registration algorithms were evaluated, revealing that the framework is robust to the multi-channel deformable registration algorithm that is used. With a single exception, all results demonstrated improvements when compared against single channel registrations using the same algorithm with mutual information.

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.

Segmentation and labeling of the ventricular system in normal pressure hydrocephalus using patch-based tissue classification and multi-atlas labeling

Normal pressure hydrocephalus (NPH) affects older adults and is thought to be caused by obstruction of the normal flow of cerebrospinal fluid (CSF). NPH typically presents with cognitive impairment, gait dysfunction, and urinary incontinence, and may account for more than five percent of all cases of dementia. Unlike most other causes of dementia, NPH can potentially be treated and the neurological dysfunction reversed by shunt surgery or endoscopic third ventriculostomy (ETV), which drain excess CSF. However, a major diagnostic challenge remains to robustly identify shunt-responsive NPH patients from patients with enlarged ventricles due to other neurodegenerative diseases. Currently, radiologists grade the severity of NPH by detailed examination and measurement of the ventricles based on stacks of 2D magnetic resonance images (MRIs). Here we propose a new method to automatically segment and label different compartments of the ventricles in NPH patients from MRIs. While this task has been achieved in healthy subjects, the ventricles in NPH are both enlarged and deformed, causing current algorithms to fail. Here we combine a patch-based tissue classification method with a registration-based multi-atlas labeling method to generate a novel algorithm that labels the lateral, third, and fourth ventricles in subjects with ventriculomegaly. The method is also applicable to other neurodegenerative diseases such as Alzheimer's disease; a condition considered in the differential diagnosis of NPH. Comparison with state of the art segmentation techniques demonstrate substantial improvements in labeling the enlarged ventricles, indicating that this strategy may be a viable option for the diagnosis and characterization of NPH.

Automatic Contrast Enhancement of Brain MR Images Using Hierarchical Correlation Histogram Analysis

Parkinson’s disease is a progressive neurodegenerative disorder that has a higher probability of occurrence in middle-aged and older adults than in the young. With the use of a computer-aided diagnosis (CAD) system, abnormal cell regions can be identified, and this identification can help medical personnel to evaluate the chance of disease. This study proposes a hierarchical correlation histogram analysis based on the grayscale distribution degree of pixel intensity by constructing a correlation histogram, that can improves the adaptive contrast enhancement for specific objects. The proposed method produces significant results during contrast enhancement preprocessing and facilitates subsequent CAD processes, thereby reducing recognition time and improving accuracy. The experimental results show that the proposed method is superior to existing methods by using two estimation image quantitative methods of PSNR and average gradient values. Furthermore, the edge information pertaining to specific cells can effectively increase the accuracy of the results.

Realistic analytical polyhedral MRI phantoms

PURPOSE

Analytical phantoms have closed form Fourier transform expressions and are used to simulate MRI acquisitions. Existing three-dimensional (3D) analytical phantoms are unable to accurately model shapes of biomedical interest. The goal of this study was to demonstrate that polyhedral analytical phantoms have closed form Fourier transform expressions and can accurately represent 3D biomedical shapes.

METHODS

The Fourier transform of a polyhedron was implemented and its accuracy in representing faceted and smooth surfaces was characterized. Realistic anthropomorphic polyhedral brain and torso phantoms were constructed and their use in simulated 3D and two-dimensional (2D) MRI acquisitions was described.

RESULTS

Using polyhedra, the Fourier transform of faceted shapes can be computed to within machine precision. Smooth surfaces can be approximated with increasing accuracy by increasing the number of facets in the polyhedron; the additional accumulated numerical imprecision of the Fourier transform of polyhedra with many faces remained small. Simulations of 3D and 2D brain and 2D torso cine acquisitions produced realistic reconstructions free of high frequency edge aliasing compared with equivalent voxelized/rasterized phantoms.

CONCLUSION

Analytical polyhedral phantoms are easy to construct and can accurately simulate shapes of biomedical interest. Magn Reson Med 76:663-678, 2016. © 2015 Wiley Periodicals, Inc.

Optical coherence tomography reflects brain atrophy in multiple sclerosis: A four-year study

OBJECTIVE

The aim of this work was to determine whether atrophy of specific retinal layers and brain substructures are associated over time, in order to further validate the utility of optical coherence tomography (OCT) as an indicator of neuronal tissue damage in patients with multiple sclerosis (MS).

METHODS

Cirrus high-definition OCT (including automated macular segmentation) was performed in 107 MS patients biannually (median follow-up: 46 months). Three-Tesla magnetic resonance imaging brain scans (including brain-substructure volumetrics) were performed annually. Individual-specific rates of change in retinal and brain measures (estimated with linear regression) were correlated, adjusting for age, sex, disease duration, and optic neuritis (ON) history.

RESULTS

Rates of ganglion cell + inner plexiform layer (GCIP) and whole-brain (r = 0.45; p < 0.001), gray matter (GM; r = 0.37; p < 0.001), white matter (WM; r = 0.28; p = 0.007), and thalamic (r = 0.38; p < 0.001) atrophy were associated. GCIP and whole-brain (as well as GM and WM) atrophy rates were more strongly associated in progressive MS (r = 0.67; p < 0.001) than relapsing-remitting MS (RRMS; r = 0.33; p = 0.007). However, correlation between rates of GCIP and whole-brain (and additionally GM and WM) atrophy in RRMS increased incrementally with step-wise refinement to exclude ON effects; excluding eyes and then patients (to account for a phenotype effect), the correlation increased to 0.45 and 0.60, respectively, consistent with effect modification. In RRMS, lesion accumulation rate was associated with GCIP (r = -0.30; p = 0.02) and inner nuclear layer (r = -0.25; p = 0.04) atrophy rates.

INTERPRETATION

Over time GCIP atrophy appears to mirror whole-brain, and particularly GM, atrophy, especially in progressive MS, thereby reflecting underlying disease progression. Our findings support OCT for clinical monitoring and as an outcome in investigative trials.

Digital Topology and Geometry in Medical Imaging: A Survey

Digital topology and geometry refers to the use of topologic and geometric properties and features for images defined in digital grids. Such methods have been widely used in many medical imaging applications, including image segmentation, visualization, manipulation, interpolation, registration, surface-tracking, object representation, correction, quantitative morphometry etc. Digital topology and geometry play important roles in medical imaging research by enriching the scope of target outcomes and by adding strong theoretical foundations with enhanced stability, fidelity, and efficiency. This paper presents a comprehensive yet compact survey on results, principles, and insights of methods related to digital topology and geometry with strong emphasis on understanding their roles in various medical imaging applications. Specifically, this paper reviews methods related to distance analysis and path propagation, connectivity, surface-tracking, image segmentation, boundary and centerline detection, topology preservation and local topological properties, skeletonization, and object representation, correction, and quantitative morphometry. A common thread among the topics reviewed in this paper is that their theory and algorithms use the principle of digital path connectivity, path propagation, and neighborhood analysis.

MR image synthesis by contrast learning on neighborhood ensembles

Automatic processing of magnetic resonance images is a vital part of neuroscience research. Yet even the best and most widely used medical image processing methods will not produce consistent results when their input images are acquired with different pulse sequences. Although intensity standardization and image synthesis methods have been introduced to address this problem, their performance remains dependent on knowledge and consistency of the pulse sequences used to acquire the images. In this paper, an image synthesis approach that first estimates the pulse sequence parameters of the subject image is presented. The estimated parameters are then used with a collection of atlas or training images to generate a new atlas image having the same contrast as the subject image. This additional image provides an ideal source from which to synthesize any other target pulse sequence image contained in the atlas. In particular, a nonlinear regression intensity mapping is trained from the new atlas image to the target atlas image and then applied to the subject image to yield the particular target pulse sequence within the atlas. Both intensity standardization and synthesis of missing tissue contrasts can be achieved using this framework. The approach was evaluated on both simulated and real data, and shown to be superior in both intensity standardization and synthesis to other established methods.

Using image synthesis for multi-channel registration of different image modalities

This paper presents a multi-channel approach for performing registration between magnetic resonance (MR) images with different modalities. In general, a multi-channel registration cannot be used when the moving and target images do not have analogous modalities. In this work, we address this limitation by using a random forest regression technique to synthesize the missing modalities from the available ones. This allows a single channel registration between two different modalities to be converted into a multi-channel registration with two mono-modal channels. To validate our approach, two openly available registration algorithms and five cost functions were used to compare the label transfer accuracy of the registration with (and without) our multi-channel synthesis approach. Our results show that the proposed method produced statistically significant improvements in registration accuracy (at an α level of 0.001) for both algorithms and all cost functions when compared to a standard multi-modal registration using the same algorithms with mutual information.

Volumetric brain analysis in neurosurgery: Part 1. Particle filter segmentation of brain and cerebrospinal fluid growth dynamics from MRI and CT images

OBJECT

Accurate edge tracing segmentation remains an incompletely solved problem in brain image analysis. The authors propose a novel algorithm using a particle filter to follow the boundary of the brain in the style often used in autonomous air and ground vehicle navigation. Their goals were to create a versatile tool to segment brain and fluid in MRI and CT images of the developing brain, lay the foundation for an intelligent automated edge tracker that is modality independent, and segment normative data from MRI that can be applied to both MRI and CT.

METHODS

Simulated MRI data sets were used to train and evaluate the particle filter segmentation algorithm. The method was then applied to produce normative growth curves for children and adolescents from 0 to 18 years of age for brain and fluid from MR images from the National Institutes of Health pediatric database and these data were compared to historical results. The authors further adapted this method for use with CT images of pediatric hydrocephalus and compared the results with hand-segmented data.

RESULTS

Segmentation of simulated MRI data with varied levels of noise (0%-9%) and spatial inhomogeneity (0%-40%) resulted in percent errors ranging from 0.06% to 5.38% for brain volume and 2.45% to 22.3% for fluid volume. The authors used this tool to create normal brain and CSF growth curves from MR images. The calculated growth curves showed excellent consistency with historical data. Additionally, compared with manual segmentation the particle filter accurately segmented brain and fluid volumes from CT scans of 5 pediatric patients with hydrocephalus (p<0.001).

CONCLUSIONS

The authors have produced the first normative brain and CSF growth curves for children and adolescents 0-18 years of age. In addition, this study includes the first use of a particle filter as an edge tracker in image segmentation and offers a semiautomatic method to segment both pediatric and adult brain data from MR and CT images. The particle filter has the potential to be further automated toward a clinical rather than research tool with both of these modalities. Because of its modality independence, it has the capability to allow CT to be a more effective diagnostic tool for neurological disorders, a task of substantial importance in emergency settings and in developing countries where CT is often the only available method of brain imaging.

Statistical normalization techniques for magnetic resonance imaging

While computed tomography and other imaging techniques are measured in absolute units with physical meaning, magnetic resonance images are expressed in arbitrary units that are difficult to interpret and differ between study visits and subjects. Much work in the image processing literature on intensity normalization has focused on histogram matching and other histogram mapping techniques, with little emphasis on normalizing images to have biologically interpretable units. Furthermore, there are no formalized principles or goals for the crucial comparability of image intensities within and across subjects. To address this, we propose a set of criteria necessary for the normalization of images. We further propose simple and robust biologically motivated normalization techniques for multisequence brain imaging that have the same interpretation across acquisitions and satisfy the proposed criteria. We compare the performance of different normalization methods in thousands of images of patients with Alzheimer's disease, hundreds of patients with multiple sclerosis, and hundreds of healthy subjects obtained in several different studies at dozens of imaging centers.

Example Based Lesion Segmentation

Automatic and accurate detection of white matter lesions is a significant step toward understanding the progression of many diseases, like Alzheimer's disease or multiple sclerosis. Multi-modal MR images are often used to segment T₂ white matter lesions that can represent regions of demyelination or ischemia. Some automated lesion segmentation methods describe the lesion intensities using generative models, and then classify the lesions with some combination of heuristics and cost minimization. In contrast, we propose a patch-based method, in which lesions are found using examples from an atlas containing multi-modal MR images and corresponding manual delineations of lesions. Patches from subject MR images are matched to patches from the atlas and lesion memberships are found based on patch similarity weights. We experiment on 43 subjects with MS, whose scans show various levels of lesion-load. We demonstrate significant improvement in Dice coefficient and total lesion volume compared to a state of the art model-based lesion segmentation method, indicating more accurate delineation of lesions.

Topology-Preserving Rigid Transformation of 2D Digital Images

We provide conditions under which 2D digital images preserve their topological properties under rigid transformations. We consider the two most common digital topology models, namely dual adjacency and well-composedness. This paper leads to the proposal of optimal preprocessing strategies that ensure the topological invariance of images under arbitrary rigid transformations. These results and methods are proved to be valid for various kinds of images (binary, gray-level, label), thus providing generic and efficient tools, which can be used in particular in the context of image registration and warping.

Comparison of 10 brain tissue segmentation methods using revisited IBSR annotations

PURPOSE

Ground-truth annotations from the well-known Internet Brain Segmentation Repository (IBSR) datasets consider Sulcal cerebrospinal fluid (SCSF) voxels as gray matter. This can lead to bias when evaluating the performance of tissue segmentation methods. In this work we compare the accuracy of 10 brain tissue segmentation methods analyzing the effects of SCSF ground-truth voxels on accuracy estimations.

MATERIALS AND METHODS

The set of methods is composed by FAST, SPM5, SPM8, GAMIXTURE, ANN, FCM, KNN, SVPASEG, FANTASM, and PVC. Methods are evaluated using original IBSR ground-truth and ranked by means of their performance on pairwise comparisons using permutation tests. Afterward, the evaluation is repeated using IBSR ground-truth without considering SCSF.

RESULTS

The Dice coefficient of all methods is affected by changes in SCSF annotations, especially on SPM5, SPM8 and FAST. When not considering SCSF voxels, SVPASEG (0.90 ± 0.01) and SPM8 (0.91 ± 0.01) are the methods from our study that appear more suitable for gray matter tissue segmentation, while FAST (0.89 ± 0.02) is the best tool for segmenting white matter tissue.

CONCLUSION

The performance and the accuracy of methods on IBSR images vary notably when not considering SCSF voxels. The fact that three of the most common methods (FAST, SPM5, and SPM8) report an important change in their accuracy suggest to consider these differences in labeling for new comparative studies.

Reconstruction of the human cerebral cortex robust to white matter lesions: method and validation

Cortical atrophy has been reported in a number of diseases, such as multiple sclerosis and Alzheimer's disease, that are also associated with white matter (WM) lesions. However, most cortical reconstruction techniques do not account for these pathologies, thereby requiring additional processing to correct for the effect of WM lesions. In this work, we introduce CRUISE(+), an automated process for cortical reconstruction from magnetic resonance brain images with WM lesions. The process extends previously well validated methods to allow for multichannel input images and to accommodate for the presence of WM lesions. We provide new validation data and tools for measuring the accuracy of cortical reconstruction methods on healthy brains as well as brains with multiple sclerosis lesions. Using this data, we validate the accuracy of CRUISE(+) and compare it to another state-of-the-art cortical reconstruction tool. Our results demonstrate that CRUISE(+) has superior performance in the cortical regions near WM lesions, and similar performance in other regions.

Automatic iceball segmentation with adapted shape priors for MRI-guided cryoablation

PURPOSE

To develop and evaluate an automatic segmentation method that extracts the 3D configuration of the ablation zone, the iceball, from images acquired during the freezing phase of MRI-guided cryoablation.

MATERIALS AND METHODS

Intraprocedural images at 63 timepoints from 13 kidney tumor cryoablation procedures were examined retrospectively. The images were obtained using a 3 Tesla wide-bore MRI scanner and axial HASTE sequence. Initialized with semiautomatically localized cryoprobes, the iceball was segmented automatically at each timepoint using the graph cut (GC) technique with adapted shape priors.

RESULTS

The average Dice Similarity Coefficients (DSC), compared with manual segmentations, were 0.88, 0.92, 0.92, 0.93, and 0.93 at 3, 6, 9, 12, and 15 min timepoints, respectively, and the average DSC of the total 63 segmentations was 0.92 ± 0.03. The proposed method improved the accuracy significantly compared with the approach without shape prior adaptation (P = 0.026). The number of probes involved in the procedure had no apparent influence on the segmentation results using our technique. The average computation time was 20 s, which was compatible with an intraprocedural setting.

CONCLUSION

Our automatic iceball segmentation method demonstrated high accuracy and robustness for practical use in monitoring the progress of MRI-guided cryoablation.

Magnetic Resonance Image Example-Based Contrast Synthesis

The performance of image analysis algorithms applied to magnetic resonance images is strongly influenced by the pulse sequences used to acquire the images. Algorithms are typically optimized for a targeted tissue contrast obtained from a particular implementation of a pulse sequence on a specific scanner. There are many practical situations, including multi-institution trials, rapid emergency scans, and scientific use of historical data, where the images are not acquired according to an optimal protocol or the desired tissue contrast is entirely missing. This paper introduces an image restoration technique that recovers images with both the desired tissue contrast and a normalized intensity profile. This is done using patches in the acquired images and an atlas containing patches of the acquired and desired tissue contrasts. The method is an example-based approach relying on sparse reconstruction from image patches. Its performance in demonstrated using several examples, including image intensity normalization, missing tissue contrast recovery, automatic segmentation, and multimodal registration. These examples demonstrate potential practical uses and also illustrate limitations of our approach.

Self-reported sleep and β-amyloid deposition in community-dwelling older adults

IMPORTANCE

Older adults commonly report disturbed sleep, and recent studies in humans and animals suggest links between sleep and Alzheimer disease biomarkers. Studies are needed that evaluate whether sleep variables are associated with neuroimaging evidence of β-amyloid (Aβ) deposition. OBJECTIVE To determine the association between self-reported sleep variables and Aβ deposition in community-dwelling older adults.

DESIGN, SETTING, AND PARTICIPANTS

Cross-sectional study of 70 adults (mean age, 76 [range, 53-91] years) from the neuroimaging substudy of the Baltimore Longitudinal Study of Aging, a normative aging study. EXPOSURE Self-reported sleep variables. MAIN OUTCOMES AND MEASURES β-Amyloid burden, measured by carbon 11-labeled Pittsburgh compound B positron emission tomography distribution volume ratios (DVRs).

RESULTS

After adjustment for potential confounders, reports of shorter sleep duration were associated with greater Aβ burden, measured by mean cortical DVR (B = 0.08 [95% CI, 0.03-0.14]; P = .005) and precuneus DVR (B = 0.11 [0.03-0.18]; P = .007). Reports of lower sleep quality were associated with greater Aβ burden measured by precuneus DVR (B = 0.08 [0.01-0.15]; P = .03).

CONCLUSIONS AND RELEVANCE

Among community-dwelling older adults, reports of shorter sleep duration and poorer sleep quality are associated with greater Aβ burden. Additional studies with objective sleep measures are needed to determine whether sleep disturbance causes or accelerates Alzheimer disease.

Deformable templates guided discriminative models for robust 3D brain MRI segmentation

Automatically segmenting anatomical structures from 3D brain MRI images is an important task in neuroimaging. One major challenge is to design and learn effective image models accounting for the large variability in anatomy and data acquisition protocols. A deformable template is a type of generative model that attempts to explicitly match an input image with a template (atlas), and thus, they are robust against global intensity changes. On the other hand, discriminative models combine local image features to capture complex image patterns. In this paper, we propose a robust brain image segmentation algorithm that fuses together deformable templates and informative features. It takes advantage of the adaptation capability of the generative model and the classification power of the discriminative models. The proposed algorithm achieves both robustness and efficiency, and can be used to segment brain MRI images with large anatomical variations. We perform an extensive experimental study on four datasets of T1-weighted brain MRI data from different sources (1,082 MRI scans in total) and observe consistent improvement over the state-of-the-art systems.

Elevated arteriolar cerebral blood volume in prodromal Huntington's disease

BACKGROUND

Neurovascular alterations have been implicated in the pathophysiology of Huntington's disease (HD). Because arterioles are most responsive to metabolic alterations, arteriolar cerebral blood volume (CBVa) is an important indicator of cerebrovascular regulation. The objective of this pilot study was to investigate potential neurovascular (CBVa ) abnormality in prodromal-HD patients and compare it with the widely used imaging marker: brain atrophy.

METHODS

CBVa and brain volumes were measured with ultra-high-field (7.0-Telsa) magnetic resonance imaging in seven prodromal-HD patients and nine age-matched controls.

RESULTS

Cortical CBVa was elevated significantly in prodromal-HD patients compared with controls (relative difference, 38.5%; effect size, 1.48). Significant correlations were found between CBVa in the frontal cortex and genetic measures. By contrast, no significant brain atrophy was detected in the prodromal-HD patients.

CONCLUSIONS

CBVa may be abnormal in prodromal-HD, even before substantial brain atrophy occurs. Further investigation with a larger cohort and longitudinal follow-up is merited to determine whether CBVa could be used as a potential biomarker for clinical trials.

Skull stripping of neonatal brain MRI: using prior shape information with graph cuts

In this paper, we propose a novel technique for skull stripping of infant (neonatal) brain magnetic resonance images using prior shape information within a graph cut framework. Skull stripping plays an important role in brain image analysis and is a major challenge for neonatal brain images. Popular methods like the brain surface extractor (BSE) and brain extraction tool (BET) do not produce satisfactory results for neonatal images due to poor tissue contrast, weak boundaries between brain and non-brain regions, and low spatial resolution. Inclusion of prior shape information helps in accurate identification of brain and non-brain tissues. Prior shape information is obtained from a set of labeled training images. The probability of a pixel belonging to the brain is obtained from the prior shape mask and included in the penalty term of the cost function. An extra smoothness term is based on gradient information that helps identify the weak boundaries between the brain and non-brain region. Experimental results on real neonatal brain images show that compared to BET, BSE, and other methods, our method achieves superior segmentation performance for neonatal brain images and comparable performance for adult brain images.

LONGITUDINAL INTENSITY NORMALIZATION IN THE PRESENCE OF MULTIPLE SCLEROSIS LESIONS

This paper proposes a longitudinal intensity normalization algorithm for T₁-weighted magnetic resonance images of human brains in the presence of multiple sclerosis lesions, aiming towards stable and consistent longitudinal segmentations. Unlike previous longitudinal segmentation methods, we propose a 4D intensity normalization that can be used as a preprocessing step to any segmentation method. The variability in intensities arising from the relapsing and remitting nature of the multiple sclerosis lesions is modeled into an otherwise smooth intensity transform based on first order autoregressive models, resulting in smooth changes in segmentation statistics of normal tissues, while keeping the lesion information unaffected. We validated our method on both simulated and real longitudinal normal subjects and on multiple sclerosis subjects.

Cortical surface reconstruction via unified Reeb analysis of geometric and topological outliers in magnetic resonance images

In this paper we present a novel system for the automated reconstruction of cortical surfaces from T1-weighted magnetic resonance images. At the core of our system is a unified Reeb analysis framework for the detection and removal of geometric and topological outliers on tissue boundaries. Using intrinsic Reeb analysis, our system can pinpoint the location of spurious branches and topological outliers, and correct them with localized filtering using information from both image intensity distributions and geometric regularity. In this system, we have also developed enhanced tissue classification with Hessian features for improved robustness to image inhomogeneity, and adaptive interpolation to achieve sub-voxel accuracy in reconstructed surfaces. By integrating these novel developments, we have a system that can automatically reconstruct cortical surfaces with improved quality and dramatically reduced computational cost as compared with the popular FreeSurfer software. In our experiments, we demonstrate on 40 simulated MR images and the MR images of 200 subjects from two databases: the Alzheimer’s Disease Neuroimaging Initiative (ADNI) and International Consortium of Brain Mapping (ICBM), the robustness of our method in large scale studies. In comparisons with FreeSurfer, we show that our system is able to generate surfaces that better represent cortical anatomy and produce thickness features with higher statistical power in population studies.

Multi-object geodesic active contours (MOGAC)

An emerging topic is to build image segmentation systems that can segment hundreds to thousands of objects (i.e. cell segmentation\tracking, full brain parcellation, full body segmentation, etc.). Multi-object Level Set Methods (MLSM) perform this task with the benefit of sub-pixel precision. However, current implementations of MLSM are not as computationally or memory efficient as their region growing and graph cut counterparts which lack sub-pixel precision. To address this performance gap, we present a novel parallel implementation of MLSM that leverages the sparse properties of the algorithm to minimize its memory footprint for multiple objects. The new method, Multi-Object Geodesic Active Contours (MOGAC), can represent N objects with just two functions: a label mask image and unsigned distance field. The time complexity of the algorithm is shown to be O((M (power)d)/P) for M (power)d pixels and P processing units in dimension d = {2,3}, independent of the number of objects. Results are presented for 2D and 3D image segmentation problems.

System for integrated neuroimaging analysis and processing of structure

Mapping brain structure in relation to neurological development, function, plasticity, and disease is widely considered to be one of the most essential challenges for opening new lines of neuro-scientific inquiry. Recent developments with MRI analysis of structural connectivity, anatomical brain segmentation, cortical surface parcellation, and functional imaging have yielded fantastic advances in our ability to probe the neurological structure-function relationship in vivo. To date, the image analysis efforts in each of these areas have typically focused on a single modality. Here, we extend the cortical reconstruction using implicit surface evolution (CRUISE) methodology to perform efficient, consistent, and topologically correct analyses in a natively multi-parametric manner. This effort combines and extends state-of-the-art techniques to simultaneously consider and analyze structural and diffusion information alongside quantitative and functional imaging data. Robust and consistent estimates of the cortical surface extraction, cortical labeling, diffusion-inferred contrasts, diffusion tractography, and subcortical parcellation are demonstrated in a scan-rescan paradigm. Accompanying this demonstration, we present a fully automated software system complete with validation data.

Automated segmentation of the cerebellar lobules using boundary specific classification and evolution

The cerebellum is instrumental in coordinating many vital functions ranging from speech and balance to eye movement. The effect of cerebellar pathology on these functions is frequently examined using volumetric studies that depend on consistent and accurate delineation, however, no existing automated methods adequately delineate the cerebellar lobules. In this work, we describe a method we call the Automatic Classification of Cerebellar Lobules Algorithm using Implicit Multi-boundary evolution (ACCLAIM). A multiple object geometric deformable model (MGDM) enables each boundary surface of each individual lobule to be evolved under different level set speeds. An important innovation described in this work is that the speed for each lobule boundary is derived from a classifier trained specifically to identify that boundary. We compared our method to segmentations obtained using the atlas-based and multi-atlas fusion techniques, and demonstrate ACCLAIM's superior performance.

An open source multivariate framework for n-tissue segmentation with evaluation on public data

We introduce Atropos, an ITK-based multivariate n-class open source segmentation algorithm distributed with ANTs ( http://www.picsl.upenn.edu/ANTs). The Bayesian formulation of the segmentation problem is solved using the Expectation Maximization (EM) algorithm with the modeling of the class intensities based on either parametric or non-parametric finite mixtures. Atropos is capable of incorporating spatial prior probability maps (sparse), prior label maps and/or Markov Random Field (MRF) modeling. Atropos has also been efficiently implemented to handle large quantities of possible labelings (in the experimental section, we use up to 69 classes) with a minimal memory footprint. This work describes the technical and implementation aspects of Atropos and evaluates its performance on two different ground-truth datasets. First, we use the BrainWeb dataset from Montreal Neurological Institute to evaluate three-tissue segmentation performance via (1) K-means segmentation without use of template data; (2) MRF segmentation with initialization by prior probability maps derived from a group template; (3) Prior-based segmentation with use of spatial prior probability maps derived from a group template. We also evaluate Atropos performance by using spatial priors to drive a 69-class EM segmentation problem derived from the Hammers atlas from University College London. These evaluation studies, combined with illustrative examples that exercise Atropos options, demonstrate both performance and wide applicability of this new platform-independent open source segmentation tool.

Brain metabolite alterations and cognitive dysfunction in early Huntington's disease

Huntington's disease (HD) is a neurodegenerative disorder characterized by early cognitive decline that progresses at later stages to dementia and severe movement disorder. HD is caused by a cytosine-adenine-guanine triplet-repeat expansion mutation in the Huntingtin gene, allowing early diagnosis by genetic testing. This study aimed to identify the relationship of N-acetylaspartate and other brain metabolites to cognitive function in HD-mutation carriers by using high-field-strength magnetic resonance spectroscopy (MRS) at 7 Tesla. Twelve individuals with the HD mutation in premanifest or early-stage disease versus 12 healthy controls underwent (1)H magnetic resonance spectroscopy (7.2 mL voxel in the posterior cingulate cortex) at 7 Tesla, and also T1-weighted structural magnetic resonance imaging. All participants received standardized tests of cognitive functioning including the Montreal Cognitive Assessment and standardized quantified neurological examination within an hour before scanning. Individuals with the HD mutation had significantly lower posterior cingulate cortex N-acetylaspartate (-9.6%, P = .02) and glutamate (-10.1%, P = .02) levels than did controls. In contrast, in this small group, measures of brain morphology including striatal and ventricle volumes did not differ significantly. Linear regression with Montreal Cognitive Assessment scores revealed significant correlations with N-acetylaspartate (r(2) = 0.50, P = .01) and glutamate (NAA) (r(2) = 0.64, P = .002) in HD subjects. Our data suggest a relationship between reduced N-acetylaspartate and glutamate levels in the posterior cingulate cortex with cognitive decline in the early stages of HD. N-acetylaspartate and glutamate magnetic resonance spectroscopy signals of the posterior cingulate cortex region may serve as potential biomarkers of disease progression or treatment outcome in HD and other neurodegenerative disorders with early cognitive dysfunction, when structural brain changes are still minor.

Automatic Sulcal Curve Extraction with MRF Based Shape Prior

Extracting and labeling sulcal curves on the human cerebral cortex is important for many neuroscience studies, however manually annotating the sulcal curves is a time-consuming task. In this paper, we present an automatic sulcal curve extraction method by registering a set of dense landmark points representing the sulcal curves to the subject cortical surface. A Markov random field is used to model the prior distribution of these landmark points, with short edges in the graph preserving the curve structure and long edges modeling the global context of the curves. Our approach is validated using a leave-one-out strategy of training and evaluation on fifteen cortical surfaces, and a quantitative error analysis on the extracted major sulcal curves.

Impaired cortico-striatal functional connectivity in prodromal Huntington's Disease

Huntington's Disease (HD) is a neurodegenerative disease caused by a CAG triplet-repeat expansion-mutation in the Huntingtin gene. Subjects at risk for HD can be identified by genetic testing in the prodromal phase. Structural changes of basal-ganglia nuclei such as the caudate nucleus are well-replicated findings observable early in prodromal-HD subjects and may be preceded by distinct functional alterations of cortico-striatal circuits. This study aims to assess functional integrity of the motor system as a cortico-striatal circuit with particular clinical relevance in HD. Ten subjects in the prodromal phase of HD and ten matched controls were administered blood oxygen level dependent (BOLD) functional magnetic resonance imaging (fMRI) at rest (3T). Functional connectivity was measured as synchrony of BOLD activity between the caudate nucleus and thirteen cortical brain regions (seeds). Basal-ganglia volumes were assessed as established markers of disease progression in prodromal-HD. Linear regression analysis was performed to test for a relationship between structural changes and group differences in functional connectivity. Prodromal-HD subjects showed reduced BOLD synchrony between two seeds in the premotor cortex (BA6) and the caudate nucleus. While similar effect sizes could be observed for reduced basal-ganglia volumes and differences in functional connectivity, coefficients of determination indicate a moderate relationship between functional connectivity and striatal atrophy. Our data show reduced cortico-striatal functional connectivity at rest in prodromal-HD and suggest a relation to early structural brain changes. Additional longitudinal studies are necessary to elucidate the temporal relationship between functional alterations and earliest structural brain changes in prodromal-HD.

Simultaneous Cortical Surface Labeling and Sulcal Curve Extraction

Automatic labeling of the gyri and sulci on the cortical surface is important for studying cortical morphology and brain functions within populations. A method to simultaneously label gyral regions and extract sulcal curves is proposed. Assuming that the gyral regions parcellate the whole cortical surface into contiguous regions with certain fixed topology, the proposed method labels the subject cortical surface by deformably registering a network of curves that form the boundary of gyral regions to the subject cortical surface. In the registration process, the curves are encouraged to follow the fine details of the sulcal geometry and to observe the shape statistics learned from training data. Using the framework of probabilistic point set registration methods, the proposed algorithm finds the sulcal curve network that maximizes the posterior probability by Expectation-Maximization (EM). The automatic labeling method was evaluated on 15 cortical surfaces using a leave-one-out strategy. Quantitative error analysis is carried out on both labeled regions and major sulcal curves.

Consistent segmentation using a Rician classifier

Several popular classification algorithms used to segment magnetic resonance brain images assume that the image intensities, or log-transformed intensities, satisfy a finite Gaussian mixture model. In these methods, the parameters of the mixture model are estimated and the posterior probabilities for each tissue class are used directly as soft segmentations or combined to form a hard segmentation. It is suggested and shown in this paper that a Rician mixture model fits the observed data better than a Gaussian model. Accordingly, a Rician mixture model is formulated and used within an expectation maximization (EM) framework to yield a new tissue classification algorithm called Rician Classifier using EM (RiCE). It is shown using both simulated and real data that RiCE yields comparable or better performance to that of algorithms based on the finite Gaussian mixture model. As well, we show that RiCE yields more consistent segmentation results when used on images of the same individual acquired with different T1-weighted pulse sequences. Therefore, RiCE has the potential to stabilize segmentation results in brain studies involving heterogeneous acquisition sources as is typically found in both multi-center and longitudinal studies.

Robust statistical fusion of image labels

Image labeling and parcellation (i.e., assigning structure to a collection of voxels) are critical tasks for the assessment of volumetric and morphometric features in medical imaging data. The process of image labeling is inherently error prone as images are corrupted by noise and artifacts. Even expert interpretations are subject to subjectivity and the precision of the individual raters. Hence, all labels must be considered imperfect with some degree of inherent variability. One may seek multiple independent assessments to both reduce this variability and quantify the degree of uncertainty. Existing techniques have exploited maximum a posteriori statistics to combine data from multiple raters and simultaneously estimate rater reliabilities. Although quite successful, wide-scale application has been hampered by unstable estimation with practical datasets, for example, with label sets with small or thin objects to be labeled or with partial or limited datasets. As well, these approaches have required each rater to generate a complete dataset, which is often impossible given both human foibles and the typical turnover rate of raters in a research or clinical environment. Herein, we propose a robust approach to improve estimation performance with small anatomical structures, allow for missing data, account for repeated label sets, and utilize training/catch trial data. With this approach, numerous raters can label small, overlapping portions of a large dataset, and rater heterogeneity can be robustly controlled while simultaneously estimating a single, reliable label set and characterizing uncertainty. The proposed approach enables many individuals to collaborate in the construction of large datasets for labeling tasks (e.g., human parallel processing) and reduces the otherwise detrimental impact of rater unavailability.