Multi-atlas tool for automated segmentation of brain gray matter nuclei and quantification of their magnetic susceptibility.

Quantification of tissue magnetic susceptibility using MRI offers a non-invasive measure of important tissue components in the brain, such as iron and myelin, potentially providing valuable information about normal and pathological conditions during aging. Despite many advances made in recent years on imaging techniques of quantitative susceptibility mapping (QSM), accurate and robust automated segmentation tools for QSM images that can help generate universal and sharable susceptibility measures in a biologically meaningful set of structures are still not widely available. In the present study, we developed an automated process to segment brain nuclei and quantify tissue susceptibility in these regions based on a susceptibility multi-atlas library, consisting of 10 atlases with T1-weighted images, gradient echo (GRE) magnitude images and QSM images of brains with different anatomic patterns. For each atlas in this library, 10 regions of interest in iron-rich deep gray matter structures that are better defined by QSM contrast were manually labeled, including caudate, putamen, globus pallidus internal/external, thalamus, pulvinar, subthalamic nucleus, substantia nigra, red nucleus and dentate nucleus in both left and right hemispheres. We then tested different pipelines using different combinations of contrast channels to bring the set of labels from the multi-atlases to each target brain and compared them with the gold standard manual delineation. The results showed that the segmentation accuracy using dual contrasts QSM/T1 pipeline outperformed other dual-contrast or single-contrast pipelines. The dice values of 0.77 ±â€¯0.09 using the QSM/T1 multi-atlas pipeline rivaled with the segmentation reliability obtained from multiple evaluators with dice values of 0.79 ±â€¯0.07 and gave comparable or superior performance in segmenting subcortical nuclei in comparison with standard FSL FIRST or recent multi-atlas package of volBrain. The segmentation performance of the QSM/T1 multi-atlas was further tested on QSM images acquired using different acquisition protocols and platforms and showed good reliability and reproducibility with average dice of 0.79 ±â€¯0.08 to manual labels and 0.89 ±â€¯0.04 in an inter-protocol manner. The extracted quantitative magnetic susceptibility values in the deep gray matter nuclei also correlated well between different protocols with inter-protocol correlation constants all larger than 0.97. Such reliability and performance was ultimately validated in an external dataset acquired at another study site with consistent susceptibility measures obtained using the QSM/T1 multi-atlas approach in comparison to those using manual delineation. In summary, we designed a susceptibility multi-atlas tool for automated and reliable segmentation of QSM images and for quantification of magnetic susceptibilities. It is publicly available through our cloud-based platform (www.mricloud.org). Further improvement on the performance of this multi-atlas tool is expected by increasing the number of atlases in the future.

ASL-MRICloud: An online tool for the processing of ASL MRI data.

Arterial spin labeling (ASL) MRI is increasingly used in research and clinical settings. The purpose of this work is to develop a cloud-based tool for ASL data processing, referred to as ASL-MRICloud, which may be useful to the MRI community. In contrast to existing ASL toolboxes, which are based on software installation on the user's local computer, ASL-MRICloud uses a web browser for data upload and results download, and the computation is performed on the remote server. As such, this tool is independent of the user's operating system, software version, and CPU speed. The ASL-MRICloud tool was implemented to be compatible with data acquired by scanners from all major MRI manufacturers, is capable of processing several common forms of ASL, including pseudo-continuous ASL and pulsed ASL, and can process single-delay and multi-delay ASL data. The outputs of ASL-MRICloud include absolute and relative values of cerebral blood flow, arterial transit time, voxel-wise masks indicating regions with potential hyper-perfusion and hypo-perfusion, and an image quality index. The ASL tool is also integrated with a T1 -based brain segmentation and normalization tool in MRICloud to allow generation of parametric maps in standard brain space as well as region-of-interest values. The tool was tested on a large data set containing 309 ASL scans as well as on publicly available ASL data from the Alzheimer's Disease Neuroimaging Initiative (ADNI) study.

Resource atlases for multi-atlas brain segmentations with multiple ontology levels based on T1-weighted MRI.

Technologies for multi-atlas brain segmentation of T1-weighted MRI images have rapidly progressed in recent years, with highly promising results. This approach, however, relies on a large number of atlases with accurate and consistent structural identifications. Here, we introduce our atlas inventories (n=90), which cover ages 4-82years with unique hierarchical structural definitions (286 structures at the finest level). This multi-atlas library resource provides the flexibility to choose appropriate atlases for various studies with different age ranges and structure-definition criteria. In this paper, we describe the details of the atlas resources and demonstrate the improved accuracy achievable with a dynamic age-matching approach, in which atlases that most closely match the subject's age are dynamically selected. The advanced atlas creation strategy, together with atlas pre-selection principles, is expected to support the further development of multi-atlas image segmentation.

A Bayesian approach to the creation of a study-customized neonatal brain atlas.

Atlas-based image analysis (ABA), in which an anatomical "parcellation map" is used for parcel-by-parcel image quantification, is widely used to analyze anatomical and functional changes related to brain development, aging, and various diseases. The parcellation maps are often created based on common MRI templates, which allow users to transform the template to target images, or vice versa, to perform parcel-by-parcel statistics, and report the scientific findings based on common anatomical parcels. The use of a study-specific template, which represents the anatomical features of the study population better than common templates, is preferable for accurate anatomical labeling; however, the creation of a parcellation map for a study-specific template is extremely labor intensive, and the definitions of anatomical boundaries are not necessarily compatible with those of the common template. In this study, we employed a volume-based template estimation (VTE) method to create a neonatal brain template customized to a study population, while keeping the anatomical parcellation identical to that of a common MRI atlas. The VTE was used to morph the standardized parcellation map of the JHU-neonate-SS atlas to capture the anatomical features of a study population. The resultant "study-customized" T1-weighted and diffusion tensor imaging (DTI) template, with three-dimensional anatomical parcellation that defined 122 brain regions, was compared with the JHU-neonate-SS atlas, in terms of the registration accuracy. A pronounced increase in the accuracy of cortical parcellation and superior tensor alignment were observed when the customized template was used. With the customized atlas-based analysis, the fractional anisotropy (FA) detected closely approximated the manual measurements. This tool provides a solution for achieving normalization-based measurements with increased accuracy, while reporting scientific findings in a consistent framework.

Tools for multiple granularity analysis of brain MRI data for individualized image analysis.

Voxel-based analysis is widely used for quantitative analysis of brain MRI. While this type of analysis provides the highest granularity level of spatial information (i.e., each voxel), the sheer number of voxels and noisy information from each voxel often lead to low sensitivity for detection of abnormalities. To ameliorate this issue, granularity reduction is commonly performed by applying isotropic spatial filtering. This study proposes a systematic reduction of the spatial information using ontology-based hierarchical structural relationships. The 254 brain structures were first defined in multiple (n=29) geriatric atlases. The multiple atlases were then applied to T1-weighted MR images of each subject's data for automated brain parcellation and five levels of ontological relationships were established, which further reduced the spatial dimension to as few as 11 structures. At each ontology level, the amount of atrophy was evaluated, providing a unique view of low-granularity analysis. This reduction of spatial information allowed us to investigate the anatomical features of each patient, demonstrated in an Alzheimer's disease group.

Diffeomorphic brain mapping based on T1-weighted images: improvement of registration accuracy by multichannel mapping.

PURPOSE: To improve image registration accuracy in neurodegenerative populations. MATERIALS AND METHODS: This study used primary progressive aphasia, aged control, and young control T1-weighted images. Mapping to a template image was performed using single-channel Large Deformation Diffeomorphic Metric Mapping (LDDMM), a dual-channel method with ventricular anatomy in the second channel, and a dual-channel with appendage method, which utilized a priori knowledge of template ventricular anatomy in the deformable atlas. RESULTS: Our results indicated substantial improvement in the registration accuracy over single-contrast-based brain mapping, mainly in the lateral ventricles and regions surrounding them. Dual-channel mapping significantly (P < 0.001) reduced the number of misclassified lateral ventricle voxels (based on a manually defined reference) over single-channel mapping. The dual-channel (w/appendage) method further reduced (P < 0.001) misclassification over the dual-channel method, indicating that the appendage provides more accurate anatomical correspondence for deformation. CONCLUSION: Brain anatomical mapping by shape normalization is widely used for quantitative anatomical analysis. However, in many geriatric and neurodegenerative disorders, severe tissue atrophy poses a unique challenge for accurate mapping of voxels, especially around the lateral ventricles. In this study we demonstrate our ability to improve mapping accuracy by incorporating ventricular anatomy in LDDMM and by utilizing a priori knowledge of ventricular anatomy in the deformable atlas.

Early amygdala and ERC atrophy linked to 3D reconstruction of rostral neurofibrillary tau tangle pathology in Alzheimer's disease.

Previous research has emphasized the unique impact of Alzheimer's Disease (AD) pathology on the medial temporal lobe (MTL), a reflection that tau pathology is particularly striking in the entorhinal and transentorhinal cortex (ERC, TEC) early in the course of disease. However, other brain regions are affected by AD pathology during its early phases. Here, we use longitudinal diffeomorphometry to measure the atrophy rate from MRI of the amygdala compared with that in the ERC and TEC in cognitively unimpaired (CU) controls, CU individuals who progressed to mild cognitive impairment (MCI), and individuals with MCI who progressed to dementia of the AD type (DAT), using a dataset from the Alzheimer's Disease Neuroimaging Initiative (ADNI). Our results show significantly higher atrophy rates of the amygdala in both groups of 'converters' (CU→MCI, MCI→DAT) compared to controls, with rates of volume loss comparable to rates of thickness loss in the ERC and TEC. We localize atrophy within the amygdala within each of these groups using fixed effects modeling. Controlling for the familywise error rate highlights the medial regions of the amygdala as those with significantly higher atrophy in both groups of converters than in controls. Using our recently developed method, referred to as Projective LDDMM, we map measures of neurofibrillary tau tangles (NFTs) from digital pathology to MRI atlases and reconstruct dense 3D spatial distributions of NFT density within regions of the MTL. The distribution of NFTs is consistent with the spatial distribution of MR measured atrophy rates, revealing high densities (and atrophy) in the amygdala (particularly medial), ERC, and rostral third of the MTL. The similarity of the location of NFTs in AD and shape changes in a well-defined clinical population suggests that amygdalar atrophy rate, as measured through MRI may be a viable biomarker for AD.

Fornix integrity and hippocampal volume predict memory decline and progression to Alzheimer's disease.

BACKGROUND: The fornix is the predominant outflow tract of the hippocampus, a brain region known to be affected early in the course of Alzheimer's disease (AD). The aims of the present study were to: (1) examine the cross-sectional relationship between fornix diffusion tensor imaging (DTI) measurements (fractional anisotropy [FA], mean diffusivity [MD], axial diffusivity, and radial diffusivity), hippocampal volume, and memory performance, and (2) compare fornix DTI measures with hippocampal volumes as predictors of progression and transition from amnestic mild cognitive impairment to AD dementia. METHODS: Twenty-three mild cognitive impairment participants for whom hippocampal volumetry and DTI were conducted at baseline received detailed evaluations at baseline; 3, 6, and 12 months; and 2.5 years. Six participants converted to AD over the follow-up period. Fornix and posterior cingulum DTI measurements and hippocampal volumes were ascertained using manual measures. Random effects models assessed each of the neuroimaging measures as predictors of decline on the Mini-Mental State Examination, Clinical Dementia Rating-sum of boxes, and memory z scores; receiver operating characteristic analyses examined the predictive value for conversion to AD. RESULTS: There was a significant correlation between fornix FA and hippocampal volumes. However, only the fornix measurements (FA, MD, radial diffusivity, and axial diffusivity) were cross-sectionally correlated with memory z scores. Both fornix FA and hippocampal volumes were predictive of memory decline. Individually, fornix FA and MD and hippocampal volumes were very good predictors of progression, with likelihood ratios >83, and better than 90% accuracy. CONCLUSION: Fornix FA both cross-sectionally correlated with and longitudinally predicted memory decline and progression to AD. Manually drawn region of interest within the fornix shows promise comparable with hippocampal volume as a predictive biomarker of progression, and this finding warrants replication in a larger study.

Brain-derived neurotrophic factor Val66Met variant on brain volumes in infants.

The brain-derived neurotrophic factor (BDNF) has many important roles in neurogenesis and neuronal health. BDNF is also involved in learning and memory. Individuals with BDNF-Val66Met variant (Met +) are at higher risk for neuropsychiatric disorders and have smaller hippocampi and amgydalae compared to those without this variant (Met -). Whether these smaller brain volumes are already present at birth is unknown and were evaluated. 66 newborn infants were genotyped for BDNF-rs6265 and had brain MRI scans. The T1-weighted images were automatically parcellated for hippocampus and amygdala, as well as the intracranial volume (ICV), total brain volume, total gray and white matter, using a multi-atlas label fusion method implemented in the MRICloud ( https://braingps.anatomyworks.org ). The segmented brain volumes were normalized to the ICV for group comparisons. The two infant groups were not different in their demographics and birth characteristics. However, compared to Met - infants, the Met + infants had smaller hippocampi (p = 0.013), smaller amygdalae (p = 0.041), and less steep age-related declines in total brain volume and % white matter volume. The smaller relative hippocampal and amygdala volumes in Met + infants suggest that the Met + genotype affected prenatal developmental processes. In addition, the slower age-dependent declines in the relative total brain and white matter volumes of the Met + group in this cross-sectional dataset suggest the BDNF-Val66Met variant might have an ongoing negative influence on the postnatal developmental processes.

A 7 Tesla Amygdalar-Hippocampal Shape Analysis of Lithium Response in Bipolar Disorder.

Research to discover clinically useful predictors of lithium response in patients with bipolar disorder has largely found them to be elusive. We demonstrate here that detailed neuroimaging may have the potential to fill this important gap in mood disorder therapeutics. Lithium treatment and bipolar disorder have both been shown to affect anatomy of the hippocampi and amygdalae but there is no consensus on the nature of their effects. We aimed to investigate structural surface anatomy changes in amygdala and hippocampus correlated with treatment response in bipolar disorder. Patients with bipolar disorder (N = 14) underwent lithium treatment, were classified by response status at acute and long-term time points, and scanned with 7 Tesla structural MRI. Large Deformation Diffeomorphic Metric Mapping was applied to detect local differences in hippocampal and amygdalar anatomy between lithium responders and non-responders. Anatomy was also compared to 21 healthy comparison participants. A patch of the ventral surface of the left hippocampus was found to be significantly atrophied in non-responders as compared to responders at the acute time point and was associated at a trend-level with long-term response status. We did not detect an association between response status and surface anatomy of the right hippocampus or amygdala. To the best of our knowledge, this is the first shape analysis of hippocampus and amygdala in bipolar disorder using 7 Tesla MRI. These results can inform future work investigating possible neuroimaging predictors of lithium response in bipolar disorder.

Plasma ceramides are altered in mild cognitive impairment and predict cognitive decline and hippocampal volume loss.

BACKGROUND: A blood-based biomarker of Alzheimer's disease (AD) would be superior to cerebrospinal fluid (CSF) and neuroimaging measures in terms of cost, invasiveness, and feasibility for repeated measures. We previously reported that blood ceramides varied in relation to timing of memory impairment in a population-based study. The present objective was to examine whether plasma ceramides varied by AD severity in a well-characterized clinic sample and were associated with cognitive decline and hippocampal volume loss over 1 year. METHODS: Participants included 25 normal controls (NC), 17 amnestic Mild Cognitive Impairment (MCI), and 21 early probable AD. A thorough neuropsychological battery and neuroimaging with hippocampal volume determination were conducted at baseline and 1 year later. Plasma ceramides were assayed at baseline using high performance liquid chromatography coupled electrospray ionization tandem mass spectrometry. RESULTS: Although all saturated ceramides were lower in MCI compared with AD at baseline, ceramides C22:0 and C24:0 were significantly lower in the MCI group compared with both NC and AD groups (P < .01). Ceramide levels did not differ (P > .05) in AD versus NC. There were no cross-sectional associations between ceramides C22:0 and C24:0 and either cognitive performance or hippocampal volume among any group. However, among the MCI group, higher baseline ceramide C22:0 and C24:0 levels were predictive of cognitive decline and hippocampal volume loss 1 year later. CONCLUSION: Results suggest that very long-chain plasma ceramides C22:0 and C24:0 are altered in MCI and predict memory loss and right hippocampal volume loss among subjects with MCI. These plasma ceramides may be early indicators of AD progression.

Multi-contrast large deformation diffeomorphic metric mapping for diffusion tensor imaging.

Diffusion tensor imaging (DTI) can reveal detailed white matter anatomy and has the potential to detect abnormalities in specific white matter structures. Such detection and quantification are, however, not straightforward. The voxel-based analysis after image normalization is one of the most widely used methods for quantitative image analyses. To apply this approach to DTI, it is important to examine if structures in the white matter are well registered among subjects, which would be highly dependent on employed algorithms for normalization. In this paper, we evaluate the accuracy of normalization of DTI data using a highly elastic transformation algorithm, called large deformation diffeomorphic metric mapping. After simulation-based validation of the algorithm, DTI data from normal subjects were used to measure the registration accuracy. To examine the impact of morphological abnormalities on the accuracy, the algorithm was also tested using data from Alzheimer's disease (AD) patients with severe brain atrophy. The accuracy level was measured by using manual landmark-based white matter matching and surface-based brain and ventricle matching as gold standard. To improve the accuracy level, cascading and multi-contrast approaches were developed. The accuracy level for the white matter was 1.88+/-0.55 and 2.19+/-0.84 mm for the measured locations in the controls and patients, respectively.

A Multi-Atlas Label Fusion Tool for Neonatal Brain MRI Parcellation and Quantification.

Structure-by-structure analysis, in which the brain magnetic resonance imaging (MRI) is parcellated based on its anatomical units, is widely used to investigate chronological changes in morphology or signal intensity during normal development, as well as to identify the alterations seen in various diseases or conditions. The multi-atlas label fusion (MALF) method is considered a highly accurate parcellation approach, and anticipated for clinical application to quantitatively evaluate early developmental processes. However, the current MALF methods, which are designed for neonatal brain segmentations, are not widely available. In this study, we developed a T1-weighted, neonatal, multi-atlas repository and integrated it into the MALF-based brain segmentation tools in the cloud-based platform, MRICloud. The cloud platform ensures users instant access to the advanced MALF tool for neonatal brains, with no software or installation requirements for the client. The Web platform by braingps.mricloud.org will eliminate the dependence on a particular operating system (eg, Windows, Macintosh, or Linux) and the requirement for high computational performance of the user's computers. The MALF-based, fully automated, image parcellation could achieve excellent agreement with manual parcellation, and the whole and regional brain volumes quantified through this method demonstrated developmental trajectories comparable to those from a previous publication. This solution will make the latest MALF tools readily available to users, with minimum barriers, and will expedite and accelerate advancements in developmental neuroscience research, neonatology, and pediatric neuroradiology.

Test-retest reproducibility of a multi-atlas automated segmentation tool on multimodality brain MRI.

INTRODUCTION: The increasing use of large sample sizes for population and personalized medicine requires high-throughput tools for imaging processing that can handle large amounts of data with diverse image modalities, perform a biologically meaningful information reduction, and result in comprehensive quantification. Exploring the reproducibility of these tools reveals the specific strengths and weaknesses that heavily influence the interpretation of results, contributing to transparence in science. METHODS: We tested-retested the reproducibility of MRICloud, a free automated method for whole-brain, multimodal MRI segmentation and quantification, on two public, independent datasets of healthy adults. RESULTS: The reproducibility was extremely high for T1-volumetric analysis, high for diffusion tensor images (DTI) (however, regionally variable), and low for resting-state fMRI. CONCLUSION: In general, the reproducibility of the different modalities was slightly superior to that of widely used software. This analysis serves as a normative reference for planning samples and for the interpretation of structure-based MRI studies.

Correction of B0 susceptibility induced distortion in diffusion-weighted images using large-deformation diffeomorphic metric mapping.

Geometric distortion caused by B0 inhomogeneity is one of the most important problems for diffusion-weighted images (DWI) using single-shot, echo planar imaging (SS-EPI). In this study, large-deformation, diffeomorphic metric mapping (LDDMM) algorithm has been tested for the correction of geometric distortion in diffusion tensor images (DTI). Based on data from nine normal subjects, the amount of distortion caused by B0 susceptibility in the 3-T magnet was characterized. The distortion quality was validated by manually placing landmarks in the target and DTI images before and after distortion correction. The distortion was found to be up to 15 mm in the population-averaged map and could be more than 20 mm in individual images. Both qualitative demonstration and quantitative statistical results suggest that the highly elastic geometric distortion caused by spatial inhomogeneity of the B0 field in DTI using SS-EPI can be effectively corrected by LDDMM. This postprocessing method is especially useful for correcting existent DTI data without phase maps.

Atlas pre-selection strategies to enhance the efficiency and accuracy of multi-atlas brain segmentation tools.

Multi-atlas brain segmentation of human brain MR images allows quantification research in structural neuroimaging. To achieve high accuracy and computational efficiency of segmentation relies on a custom subset of atlases for each target subject. However, the criterion for atlas pre-selection remains an open question. In this study, two atlas pre-selection approaches based on location-based feature matching were proposed and compared to random and mutual information-based methods using a database of 47 atlases. A varying number of atlases ranked top with hierarchical structural granularity were compared using Dice overlap. The results indicated that the proposed 4L approach consistently led to the highest level of accuracy at a given number of employed atlases in both adult and geriatric populations. In addition, the proposed two methods (4L and LV) can reduce 20 times computational time compared with the stereotypical mutual information-based method. Our pre-selection strategy would provide better segmentation performance in terms of both accuracy and efficiency. The proposed atlas pre-selection will be further implemented into our online automatic brain image segmentation system (www.mricloud.org).

Large deformation diffeomorphism and momentum based hippocampal shape discrimination in dementia of the Alzheimer type.

In large-deformation diffeomorphic metric mapping (LDDMM), the diffeomorphic matching of images are modeled as evolution in time, or a flow, of an associated smooth velocity vector field v controlling the evolution. The initial momentum parameterizes the whole geodesic and encodes the shape and form of the target image. Thus, methods such as principal component analysis (PCA) of the initial momentum leads to analysis of anatomical shape and form in target images without being restricted to small-deformation assumption in the analysis of linear displacements. We apply this approach to a study of dementia of the Alzheimer type (DAT). The left hippocampus in the DAT group shows significant shape abnormality while the right hippocampus shows similar pattern of abnormality. Further, PCA of the initial momentum leads to correct classification of 12 out of 18 DAT subjects and 22 out of 26 control subjects.

Direct estimation of patient attributes from anatomical MRI based on multi-atlas voting.

MRI brain atlases are widely used for automated image segmentation, and in particular, recent developments in multi-atlas techniques have shown highly accurate segmentation results. In this study, we extended the role of the atlas library from mere anatomical reference to a comprehensive knowledge database with various patient attributes, such as demographic, functional, and diagnostic information. In addition to using the selected (heavily-weighted) atlases to achieve high segmentation accuracy, we tested whether the non-anatomical attributes of the selected atlases could be used to estimate patient attributes. This can be considered a context-based image retrieval (CBIR) approach, embedded in the multi-atlas framework. We first developed an image similarity measurement to weigh the atlases on a structure-by-structure basis, and then, the attributes of the multiple atlases were weighted to estimate the patient attributes. We tested this concept first by estimating age in a normal population; we then performed functional and diagnostic estimations in Alzheimer's disease patients. The accuracy of the estimated patient attributes was measured against the actual clinical data, and the performance was compared to conventional volumetric analysis. The proposed CBIR framework by multi-atlas voting would be the first step toward a knowledge-based support system for quantitative radiological image reading and diagnosis.

Evaluation of Cross-Protocol Stability of a Fully Automated Brain Multi-Atlas Parcellation Tool.

Brain parcellation tools based on multiple-atlas algorithms have recently emerged as a promising method with which to accurately define brain structures. When dealing with data from various sources, it is crucial that these tools are robust for many different imaging protocols. In this study, we tested the robustness of a multiple-atlas, likelihood fusion algorithm using Alzheimer's Disease Neuroimaging Initiative (ADNI) data with six different protocols, comprising three manufacturers and two magnetic field strengths. The entire brain was parceled into five different levels of granularity. In each level, which defines a set of brain structures, ranging from eight to 286 regions, we evaluated the variability of brain volumes related to the protocol, age, and diagnosis (healthy or Alzheimer's disease). Our results indicated that, with proper pre-processing steps, the impact of different protocols is minor compared to biological effects, such as age and pathology. A precise knowledge of the sources of data variation enables sufficient statistical power and ensures the reliability of an anatomical analysis when using this automated brain parcellation tool on datasets from various imaging protocols, such as clinical databases.

Segmentation of brain magnetic resonance images based on multi-atlas likelihood fusion: testing using data with a broad range of anatomical and photometric profiles.

We propose a hierarchical pipeline for skull-stripping and segmentation of anatomical structures of interest from T1-weighted images of the human brain. The pipeline is constructed based on a two-level Bayesian parameter estimation algorithm called multi-atlas likelihood fusion (MALF). In MALF, estimation of the parameter of interest is performed via maximum a posteriori estimation using the expectation-maximization (EM) algorithm. The likelihoods of multiple atlases are fused in the E-step while the optimal estimator, a single maximizer of the fused likelihoods, is then obtained in the M-step. There are two stages in the proposed pipeline; first the input T1-weighted image is automatically skull-stripped via a fast MALF, then internal brain structures of interest are automatically extracted using a regular MALF. We assess the performance of each of the two modules in the pipeline based on two sets of images with markedly different anatomical and photometric contrasts; 3T MPRAGE scans of pediatric subjects with developmental disorders vs. 1.5T SPGR scans of elderly subjects with dementia. Evaluation is performed quantitatively using the Dice overlap as well as qualitatively via visual inspections. As a result, we demonstrate subject-level differences in the performance of the proposed pipeline, which may be accounted for by age, diagnosis, or the imaging parameters (particularly the field strength). For the subcortical and ventricular structures of the two datasets, the hierarchical pipeline is capable of producing automated segmentations with Dice overlaps ranging from 0.8 to 0.964 when compared with the gold standard. Comparisons with other representative segmentation algorithms are presented, relative to which the proposed hierarchical pipeline demonstrates comparative or superior accuracy.