The amygdala is composed of multiple nuclei with unique functions and connections in the limbic system and to the rest of the brain. However, standard in vivo neuroimaging tools to automatically ...delineate the amygdala into its multiple nuclei are still rare. By scanning postmortem specimens at high resolution (100–150µm) at 7T field strength (n = 10), we were able to visualize and label nine amygdala nuclei (anterior amygdaloid, cortico-amygdaloid transition area; basal, lateral, accessory basal, central, cortical medial, paralaminar nuclei). We created an atlas from these labels using a recently developed atlas building algorithm based on Bayesian inference. This atlas, which will be released as part of FreeSurfer, can be used to automatically segment nine amygdala nuclei from a standard resolution structural MR image. We applied this atlas to two publicly available datasets (ADNI and ABIDE) with standard resolution T1 data, used individual volumetric data of the amygdala nuclei as the measure and found that our atlas i) discriminates between Alzheimer's disease participants and age-matched control participants with 84% accuracy (AUC=0.915), and ii) discriminates between individuals with autism and age-, sex- and IQ-matched neurotypically developed control participants with 59.5% accuracy (AUC=0.59). For both datasets, the new ex vivo atlas significantly outperformed (all p < .05) estimations of the whole amygdala derived from the segmentation in FreeSurfer 5.1 (ADNI: 75%, ABIDE: 54% accuracy), as well as classification based on whole amygdala volume (using the sum of all amygdala nuclei volumes; ADNI: 81%, ABIDE: 55% accuracy). This new atlas and the segmentation tools that utilize it will provide neuroimaging researchers with the ability to explore the function and connectivity of the human amygdala nuclei with unprecedented detail in healthy adults as well as those with neurodevelopmental and neurodegenerative disorders.
•We visualized 9 nuclei boundaries (anterior amygdaloid area, cortico-amygdaloid transition area; basal, lateral, accessory basal, central, cortical medial, paralaminar nuclei) using ultra-high-resolution ex vivo imaging.•Nuclei were consistent across cases and raters.•We built a segmentation atlas of the amygdala nuclei, which will be distributed with FreeSurfer.•Atlas was applied to 2 datasets and showed higher discriminability of Alzheimer's & autism than previously possible.•The atlas will provide neuroimaging researchers with the ability to test nucleus function with greater spatial specificity.
In this paper, we focus on the retrospective topology correction of surfaces. We propose a technique to accurately correct the spherical topology of cortical surfaces. Specifically, we construct a ...mapping from the original surface onto the sphere to detect topological defects as minimal nonhomeomorphic regions. The topology of each defect is then corrected by opening and sealing the surface along a set of nonseparating loops that are selected in a Bayesian framework. The proposed method is a wholly self-contained topology correction algorithm, which determines geometrically accurate, topologically correct solutions based on the magnetic resonance imaging (MRI) intensity profile and the expected local curvature. Applied to real data, our method provides topological corrections similar to those made by a trained operator
Cross-sectional magnetic resonance imaging (MRI) studies of cortical thickness and volume have shown age effects on large areas, but there are substantial discrepancies across studies regarding the ...localization and magnitude of effects. These discrepancies hinder understanding of effects of aging on brain morphometry, and limit the potential usefulness of MR in research on healthy and pathological age-related brain changes. The present study was undertaken to overcome this problem by assessing the consistency of age effects on cortical thickness across 6 different samples with a total of 883 participants. A surface-based segmentation procedure (FreeSurfer) was used to calculate cortical thickness continuously across the brain surface. The results showed consistent age effects across samples in the superior, middle, and inferior frontal gyri, superior and middle temporal gyri, precuneus, inferior and superior parietal cortices, fusiform and lingual gyri, and the temporo-parietal junction. The strongest effects were seen in the superior and inferior frontal gyri, as well as superior parts of the temporal lobe. The inferior temporal lobe and anterior cingulate cortices were relatively less affected by age. The results are discussed in relation to leading theories of cognitive aging.
Volume deficits of the hippocampus in schizophrenia have been consistently reported. However, the hippocampus is anatomically heterogeneous; it remains unclear whether certain portions of the ...hippocampus are affected more than others in schizophrenia. In this study, we aimed to determine whether volume deficits in schizophrenia are confined to specific subfields of the hippocampus and to measure the subfield volume trajectories over the course of the illness. Magnetic resonance imaging scans were obtained from Data set 1: 155 patients with schizophrenia (mean duration of illness of 7 years) and 79 healthy controls, and Data set 2: an independent cohort of 46 schizophrenia patients (mean duration of illness of 18 years) and 46 healthy controls. In addition, follow-up scans were collected for a subset of Data set 1. A novel, automated method based on an atlas constructed from ultra-high resolution, post-mortem hippocampal tissue was used to label seven hippocampal subfields. Significant cross-sectional volume deficits in the CA1, but not of the other subfields, were found in the schizophrenia patients of Data set 1. However, diffuse cross-sectional volume deficits across all subfields were found in the more chronic and ill schizophrenia patients of Data set 2. Consistent with this pattern, the longitudinal analysis of Data set 1 revealed progressive illness-related volume loss (~2-6% per year) that extended beyond CA1 to all of the other subfields. This decline in volume correlated with symptomatic worsening. Overall, these findings provide converging evidence for early atrophy of CA1 in schizophrenia, with extension to other hippocampal subfields and accompanying clinical sequelae over time.
Atlas-based approaches have demonstrated the ability to automatically identify detailed brain structures from 3-D magnetic resonance (MR) brain images. Unfortunately, the accuracy of this type of ...method often degrades when processing data acquired on a different scanner platform or pulse sequence than the data used for the atlas training. In this paper, we improve the performance of an atlas-based whole brain segmentation method by introducing an intensity renormalization procedure that automatically adjusts the prior atlas intensity model to new input data. Validation using manually labeled test datasets has shown that the new procedure improves the segmentation accuracy (as measured by the Dice coefficient) by 10% or more for several structures including hippocampus, amygdala, caudate, and pallidum. The results verify that this new procedure reduces the sensitivity of the whole brain segmentation method to changes in scanner platforms and improves its accuracy and robustness, which can thus facilitate multicenter or multisite neuroanatomical imaging studies
We present a novel skull-stripping algorithm based on a hybrid approach that combines watershed algorithms and deformable surface models. Our method takes advantage of the robustness of the former as ...well as the surface information available to the latter. The algorithm first localizes a single white matter voxel in a T1-weighted MRI image, and uses it to create a global minimum in the white matter before applying a watershed algorithm with a preflooding height. The watershed algorithm builds an initial estimate of the brain volume based on the three-dimensional connectivity of the white matter. This first step is robust, and performs well in the presence of intensity nonuniformities and noise, but may erode parts of the cortex that abut bright nonbrain structures such as the eye sockets, or may remove parts of the cerebellum. To correct these inaccuracies, a surface deformation process fits a smooth surface to the masked volume, allowing the incorporation of geometric constraints into the skull-stripping procedure. A statistical atlas, generated from a set of accurately segmented brains, is used to validate and potentially correct the segmentation, and the MRI intensity values are locally re-estimated at the boundary of the brain. Finally, a high-resolution surface deformation is performed that accurately matches the outer boundary of the brain, resulting in a robust and automated procedure. Studies by our group and others outperform other publicly available skull-stripping tools.
We present a technique for automatically assigning a neuroanatomical label to each location on a cortical surface model based on probabilistic information estimated from a manually labeled training ...set. This procedure incorporates both geometric information derived from the cortical model, and neuroanatomical convention, as found in the training set. The result is a complete labeling of cortical sulci and gyri. Examples are given from two different training sets generated using different neuroanatomical conventions, illustrating the flexibility of the algorithm. The technique is shown to be comparable in accuracy to manual labeling.
The human cerebral cortex is made up of a mosaic of structural areas, frequently referred to as Brodmann areas (BAs). Despite the widespread use of cortical folding patterns to perform ad hoc ...estimations of the locations of the BAs, little is understood regarding 1) how variable the position of a given BA is with respect to the folds, 2) whether the location of some BAs is more variable than others, and 3) whether the variability is related to the level of a BA in a putative cortical hierarchy. We use whole-brain histology of 10 postmortem human brains and surface-based analysis to test how well the folds predict the locations of the BAs. We show that higher order cortical areas exhibit more variability than primary and secondary areas and that the folds are much better predictors of the BAs than had been previously thought. These results further highlight the significance of cortical folding patterns and suggest a common mechanism for the development of the folds and the cytoarchitectonic fields.
We present the Spherical Demons algorithm for registering two spherical images. By exploiting spherical vector spline interpolation theory, we show that a large class of regularizors for the modified ...Demons objective function can be efficiently approximated on the sphere using iterative smoothing. Based on one parameter subgroups of diffeomorphisms, the resulting registration is diffeomorphic and fast. The Spherical Demons algorithm can also be modified to register a given spherical image to a probabilistic atlas. We demonstrate two variants of the algorithm corresponding to warping the atlas or warping the subject. Registration of a cortical surface mesh to an atlas mesh, both with more than 160 k nodes requires less than 5 min when warping the atlas and less than 3 min when warping the subject on a Xeon 3.2 GHz single processor machine. This is comparable to the fastest nondiffeomorphic landmark-free surface registration algorithms. Furthermore, the accuracy of our method compares favorably to the popular FreeSurfer registration algorithm. We validate the technique in two different applications that use registration to transfer segmentation labels onto a new image (1) parcellation of in vivo cortical surfaces and (2) Brodmann area localization in ex vivo cortical surfaces.
Highly accurate surface models of the cerebral cortex are becoming increasingly important as tools in the investigation of the functional organization of the human brain. The construction of such ...models is difficult using current neuroimaging technology due to the high degree of cortical folding. Even single voxel mis-classifications can result in erroneous connections being created between adjacent banks of a sulcus, resulting in a topologically inaccurate model. These topological defects cause the cortical model to no longer be homeomorphic to a sheet, preventing the accurate inflation, flattening, or spherical morphing of the reconstructed cortex. Surface deformation techniques can guarantee the topological correctness of a model, but are time-consuming and may result in geometrically inaccurate models. In order to address this need the authors have developed a technique for taking a model of the cortex, detecting and fixing the topological defects while leaving that majority of the model intact, resulting in a surface that is both geometrically accurate and topologically correct.