ANTsX neuroimaging-derived structural phenotypes of UK Biobank

UK Biobank is a large-scale epidemiological resource for investigating prospective correlations between various lifestyle, environmental, and genetic factors with health and disease progression. In addition to individual subject information obtained through surveys and physical examinations, a comprehensive neuroimaging battery consisting of multiple modalities provides imaging-derived phenotypes (IDPs) that can serve as biomarkers in neuroscience research. In this study, we augment the existing set of UK Biobank neuroimaging structural IDPs, obtained from well-established software libraries such as FSL and FreeSurfer, with related measurements acquired through the Advanced Normalization Tools Ecosystem. This includes previously established cortical and subcortical measurements defined, in part, based on the Desikan-Killiany-Tourville atlas. Also included are morphological measurements from two recent developments: medial temporal lobe parcellation of hippocampal and extra-hippocampal regions in addition to cerebellum parcellation and thickness based on the Schmahmann anatomical labeling. Through predictive modeling, we assess the clinical utility of these IDP measurements, individually and in combination, using commonly studied phenotypic correlates including age, fluid intelligence, numeric memory, and several other sociodemographic variables. The predictive accuracy of these IDP-based models, in terms of root-mean-squared-error or area-under-the-curve for continuous and categorical variables, respectively, provides comparative insights between software libraries as well as potential clinical interpretability. Results demonstrate varied performance between package-based IDP sets and their combination, emphasizing the need for careful consideration in their selection and utilization.

Neurological Effects of Repeated Blast Exposure in Special Operations Personnel

Exposure to blast overpressure has been a pervasive feature of combat-related injuries. Studies exploring the neurological correlates of repeated low-level blast exposure in career "breachers" demonstrated higher levels of tumor necrosis factor alpha (TNFα) and interleukin (IL)-6 and decreases in IL-10 within brain-derived extracellular vesicles (BDEVs). The current pilot study was initiated in partnership with the U.S. Special Operations Command (USSOCOM) to explore whether neuroinflammation is seen within special operators with prior blast exposure. Data were analyzed from 18 service members (SMs), inclusive of 9 blast-exposed special operators with an extensive career history of repeated blast exposures and 9 controls matched by age and duration of service. Neuroinflammation was assessed utilizing positron emission tomography (PET) imaging with [18F]DPA-714. Serum was acquired to assess inflammatory biomarkers within whole serum and BDEVs. The Blast Exposure Threshold Survey (BETS) was acquired to determine blast history. Both self-report and neurocognitive measures were acquired to assess cognition. Similarity-driven Multi-view Linear Reconstruction (SiMLR) was used for joint analysis of acquired data. Analysis of BDEVs indicated significant positive associations with a generalized blast exposure value (GBEV) derived from the BETS. SiMLR-based analyses of neuroimaging demonstrated exposure-related relationships between GBEV, PET-neuroinflammation, cortical thickness, and volume loss within special operators. Affected brain networks included regions associated with memory retrieval and executive functioning, as well as visual and heteromodal processing. Post hoc assessments of cognitive measures failed to demonstrate significant associations with GBEV. This emerging evidence suggests neuroinflammation may be a key feature of the brain response to blast exposure over a career in operational personnel. The common thread of neuroinflammation observed in blast-exposed populations requires further study.

ANTsX neuroimaging-derived structural phenotypes of UK Biobank

UK Biobank is a large-scale epidemiological resource for investigating prospective correlations between various lifestyle, environmental, and genetic factors with health and disease progression. In addition to individual subject information obtained through surveys and physical examinations, a comprehensive neuroimaging battery consisting of multiple modalities provides imaging-derived phenotypes (IDPs) that can serve as biomarkers in neuroscience research. In this study, we augment the existing set of UK Biobank neuroimaging structural IDPs, obtained from well-established software libraries such as FSL and FreeSurfer, with related measurements acquired through the Advanced Normalization Tools Ecosystem. This includes previously established cortical and subcortical measurements defined, in part, based on the Desikan-Killiany-Tourville atlas. Also included are morphological measurements from two recent developments: medial temporal lobe parcellation of hippocampal and extra-hippocampal regions in addition to cerebellum parcellation and thickness based on the Schmahmann anatomical labeling. Through predictive modeling, we assess the clinical utility of these IDP measurements, individually and in combination, using commonly studied phenotypic correlates including age, fluid intelligence, numeric memory, and several other sociodemographic variables. The predictive accuracy of these IDP-based models, in terms of root-mean-squared-error or area-under-the-curve for continuous and categorical variables, respectively, provides comparative insights between software libraries as well as potential clinical interpretability. Results demonstrate varied performance between package-based IDP sets and their combination, emphasizing the need for careful consideration in their selection and utilization.

Image- versus histogram-based considerations in semantic segmentation of pulmonary hyperpolarized gas images

PURPOSE

To characterize the differences between histogram-based and image-based algorithms for segmentation of hyperpolarized gas lung images.

METHODS

Four previously published histogram-based segmentation algorithms (ie, linear binning, hierarchical k-means, fuzzy spatial c-means, and a Gaussian mixture model with a Markov random field prior) and an image-based convolutional neural network were used to segment 2 simulated data sets derived from a public (n = 29 subjects) and a retrospective collection (n = 51 subjects) of hyperpolarized 129Xe gas lung images transformed by common MRI artifacts (noise and nonlinear intensity distortion). The resulting ventilation-based segmentations were used to assess algorithmic performance and characterize optimization domain differences in terms of measurement bias and precision.

RESULTS

Although facilitating computational processing and providing discriminating clinically relevant measures of interest, histogram-based segmentation methods discard important contextual spatial information and are consequently less robust in terms of measurement precision in the presence of common MRI artifacts relative to the image-based convolutional neural network.

CONCLUSIONS

Direct optimization within the image domain using convolutional neural networks leverages spatial information, which mitigates problematic issues associated with histogram-based approaches and suggests a preferred future research direction. Further, the entire processing and evaluation framework, including the newly reported deep learning functionality, is available as open source through the well-known Advanced Normalization Tools ecosystem.

Similarity-driven multi-view embeddings from high-dimensional biomedical data

Diverse, high-dimensional modalities collected in large cohorts present new opportunities for the formulation and testing of integrative scientific hypotheses. Similarity-driven multi-view linear reconstruction (SiMLR) is an algorithm that exploits inter-modality relationships to transform large scientific datasets into smaller, more well-powered and interpretable low-dimensional spaces. SiMLR contributes an objective function for identifying joint signal, regularization based on sparse matrices representing prior within-modality relationships and an implementation that permits application to joint reduction of large data matrices. We demonstrate that SiMLR outperforms closely related methods on supervised learning problems in simulation data, a multi-omics cancer survival prediction dataset and multiple modality neuroimaging datasets. Taken together, this collection of results shows that SiMLR may be applied to joint signal estimation from disparate modalities and may yield practically useful results in a variety of application domains.

Longitudinal Mapping of Cortical Thickness Measurements: An Alzheimer's Disease Neuroimaging Initiative-Based Evaluation Study

Longitudinal studies of development and disease in the human brain have motivated the acquisition of large neuroimaging data sets and the concomitant development of robust methodological and statistical tools for quantifying neurostructural changes. Longitudinal-specific strategies for acquisition and processing have potentially significant benefits including more consistent estimates of intra-subject measurements while retaining predictive power. Using the first phase of the Alzheimer's Disease Neuroimaging Initiative (ADNI-1) data, comprising over 600 subjects with multiple time points from baseline to 36 months, we evaluate the utility of longitudinal FreeSurfer and Advanced Normalization Tools (ANTs) surrogate thickness values in the context of a linear mixed-effects (LME) modeling strategy. Specifically, we estimate the residual variability and between-subject variability associated with each processing stream as it is known from the statistical literature that minimizing the former while simultaneously maximizing the latter leads to greater scientific interpretability in terms of tighter confidence intervals in calculated mean trends, smaller prediction intervals, and narrower confidence intervals for determining cross-sectional effects. This strategy is evaluated over the entire cortex, as defined by the Desikan-Killiany-Tourville labeling protocol, where comparisons are made with the cross-sectional and longitudinal FreeSurfer processing streams. Subsequent linear mixed effects modeling for identifying diagnostic groupings within the ADNI cohort is provided as supporting evidence for the utility of the proposed ANTs longitudinal framework which provides unbiased structural neuroimage processing and competitive to superior power for longitudinal structural change detection.

Functional and Structural Neuroimaging Correlates of Repetitive Low-Level Blast Exposure in Career Breachers

Combat military and civilian law enforcement personnel may be exposed to repetitive low-intensity blast events during training and operations. Persons who use explosives to gain entry (i.e., breach) into buildings are known as “breachers” or dynamic entry personnel. Breachers operate under the guidance of established safety protocols, but despite these precautions, breachers who are exposed to low-level blast throughout their careers frequently report performance deficits and symptoms to healthcare providers. Although little is known about the etiology linking blast exposure to clinical symptoms in humans, animal studies demonstrate network-level changes in brain function, alterations in brain morphology, vascular and inflammatory changes, hearing loss, and even alterations in gene expression after repeated blast exposure. To explore whether similar effects occur in humans, we collected a comprehensive data battery from 20 experienced breachers exposed to blast throughout their careers and 14 military and law enforcement controls. This battery included neuropsychological assessments, blood biomarkers, and magnetic resonance imaging measures, including cortical thickness, diffusion tensor imaging of white matter, functional connectivity, and perfusion. To better understand the relationship between repetitive low-level blast exposure and behavioral and imaging differences in humans, we analyzed the data using similarity-driven multi-view linear reconstruction (SiMLR). SiMLR is specifically designed for multiple modality statistical integration using dimensionality-reduction techniques for studies with high-dimensional, yet sparse, data (i.e., low number of subjects and many data per subject). We identify significant group effects in these data spanning brain structure, function, and blood biomarkers.

Small Animal Multivariate Brain Analysis (SAMBA) - a High Throughput Pipeline with a Validation Framework

While many neuroscience questions aim to understand the human brain, much current knowledge has been gained using animal models, which replicate genetic, structural, and connectivity aspects of the human brain. While voxel-based analysis (VBA) of preclinical magnetic resonance images is widely-used, a thorough examination of the statistical robustness, stability, and error rates is hindered by high computational demands of processing large arrays, and the many parameters involved therein. Thus, workflows are often based on intuition or experience, while preclinical validation studies remain scarce. To increase throughput and reproducibility of quantitative small animal brain studies, we have developed a publicly shared, high throughput VBA pipeline in a high-performance computing environment, called SAMBA. The increased computational efficiency allowed large multidimensional arrays to be processed in 1–3 days—a task that previously took ~1 month. To quantify the variability and reliability of preclinical VBA in rodent models, we propose a validation framework consisting of morphological phantoms, and four metrics. This addresses several sources that impact VBA results, including registration and template construction strategies. We have used this framework to inform the VBA workflow parameters in a VBA study for a mouse model of epilepsy. We also present initial efforts towards standardizing small animal neuroimaging data in a similar fashion with human neuroimaging. We conclude that verifying the accuracy of VBA merits attention, and should be the focus of a broader effort within the community. The proposed framework promotes consistent quality assurance of VBA in preclinical neuroimaging, thus facilitating the creation and communication of robust results.

Multivariate MR biomarkers better predict cognitive dysfunction in mouse models of Alzheimer's disease

To understand multifactorial conditions such as Alzheimer's disease (AD) we need brain signatures that predict the impact of multiple pathologies and their interactions. To help uncover the relationships between pathology affected brain circuits and cognitive markers we have used mouse models that represent, at least in part, the complex interactions altered in AD, while being raised in uniform environments and with known genotype alterations. In particular, we aimed to understand the relationship between vulnerable brain circuits and memory deficits measured in the Morris water maze, and we tested several predictive modeling approaches. We used in vivo manganese enhanced MRI traditional voxel based analyses to reveal regional differences in volume (morphometry), signal intensity (activity), and magnetic susceptibility (iron deposition, demyelination). These regions included hippocampus, olfactory areas, entorhinal cortex and cerebellum, as well as the frontal association area. The properties of these regions, extracted from each of the imaging markers, were used to predict spatial memory. We next used eigenanatomy, which reduces dimensionality to produce sets of regions that explain the variance in the data. For each imaging marker, eigenanatomy revealed networks underpinning a range of cognitive functions including memory, motor function, and associative learning, allowing the detection of associations between context, location, and responses. Finally, the integration of multivariate markers in a supervised sparse canonical correlation approach outperformed single predictor models and had significant correlates to spatial memory. Among a priori selected regions, expected to play a role in memory dysfunction, the fornix also provided good predictors, raising the possibility of investigating how disease propagation within brain networks leads to cognitive deterioration. Our cross-sectional results support that modeling approaches integrating multivariate imaging markers provide sensitive predictors of AD-like behaviors. Such strategies for mapping brain circuits responsible for behaviors may help in the future predict disease progression, or response to interventions.

Convolutional Neural Networks with Template-Based Data Augmentation for Functional Lung Image Quantification

RATIONALE AND OBJECTIVES

We propose an automated segmentation pipeline based on deep learning for proton lung MRI segmentation and ventilation-based quantification which improves on our previously reported methodologies in terms of computational efficiency while demonstrating accuracy and robustness. The large data requirement for the proposed framework is made possible by a novel template-based data augmentation strategy. Supporting this work is the open-source ANTsRNet-a growing repository of well-known deep learning architectures first introduced here.

MATERIALS AND METHODS

Deep convolutional neural network (CNN) models were constructed and trained using a custom multilabel Dice metric loss function and a novel template-based data augmentation strategy. Training (including template generation and data augmentation) employed 205 proton MR images and 73 functional lung MRI. Evaluation was performed using data sets of size 63 and 40 images, respectively.

RESULTS

Accuracy for CNN-based proton lung MRI segmentation (in terms of Dice overlap) was left lung: 0.93 ± 0.03, right lung: 0.94 ± 0.02, and whole lung: 0.94 ± 0.02. Although slightly less accurate than our previously reported joint label fusion approach (left lung: 0.95 ± 0.02, right lung: 0.96 ± 0.01, and whole lung: 0.96 ± 0.01), processing time is <1 second per subject for the proposed approach versus ∼30 minutes per subject using joint label fusion. Accuracy for quantifying ventilation defects was determined based on a consensus labeling where average accuracy (Dice multilabel overlap of ventilation defect regions plus normal region) was 0.94 for the CNN method; 0.92 for our previously reported method; and 0.90, 0.92, and 0.94 for expert readers.

CONCLUSION

The proposed framework yields accurate automated quantification in near real time. CNNs drastically reduce processing time after offline model construction and demonstrate significant future potential for facilitating quantitative analysis of functional lung MRI.

Amyloid beta-positive subjects exhibit longitudinal network-specific reductions in spontaneous brain activity

Amyloid beta (Aβ) deposition and cognitive decline are key features of Alzheimer's disease. The relationship between Aβ status and changes in neuronal function over time, however, remains unclear. We evaluated the effect of baseline Aβ status on reference region spontaneous brain activity (SBA-rr) using resting-state functional magnetic resonance imaging and fluorodeoxyglucose positron emission tomography in patients with mild cognitive impairment. Patients (N = 62, [43 Aβ-positive]) from the Alzheimer's Disease Neuroimaging Initiative were divided into Aβ-positive and Aβ-negative groups via prespecified cerebrospinal fluid Aβ42 or 18F-florbetapir positron emission tomography standardized uptake value ratio cutoffs measured at baseline. We analyzed interaction of biomarker-confirmed Aβ status with SBA-rr change over a 2-year period using mixed-effects modeling. SBA-rr differences between Aβ-positive and Aβ-negative subjects increased significantly over time within subsystems of the default and visual networks. Changes exhibit an interaction with memory performance over time but were independent of glucose metabolism. Results reinforce the value of resting-state functional magnetic resonance imaging in evaluating Alzheimer''s disease progression and suggest spontaneous neuronal activity changes are concomitant with cognitive decline.

A population level atlas of Mus musculus craniofacial skeleton and automated image-based shape analysis

Laboratory mice are staples for evo/devo and genetics studies. Inbred strains provide a uniform genetic background to manipulate and understand gene-environment interactions, while their crosses have been instrumental in studies of genetic architecture, integration and modularity, and mapping of complex biological traits. Recently, there have been multiple large-scale studies of laboratory mice to further our understanding of the developmental basis, evolution, and genetic control of shape variation in the craniofacial skeleton (i.e. skull and mandible). These experiments typically use micro-computed tomography (micro-CT) to capture the craniofacial phenotype in 3D and rely on manually annotated anatomical landmarks to conduct statistical shape analysis. Although the common choice for imaging modality and phenotyping provides the potential for collaborative research for even larger studies with more statistical power, the investigator (or lab-specific) nature of the data collection hampers these efforts. Investigators are rightly concerned that subtle differences in how anatomical landmarks were recorded will create systematic bias between studies that will eventually influence scientific findings. Even if researchers are willing to repeat landmark annotation on a combined dataset, different lab practices and software choices may create obstacles for standardization beyond the underlying imaging data. Here, we propose a freely available analysis system that could assist in the standardization of micro-CT studies in the mouse. Our proposal uses best practices developed in biomedical imaging and takes advantage of existing open-source software and imaging formats. Our first contribution is the creation of a synthetic template for the adult mouse craniofacial skeleton from 25 inbred strains and five F1 crosses that are widely used in biological research. The template contains a fully segmented cranium, left and right hemi-mandibles, endocranial space, and the first few cervical vertebrae. We have been using this template in our lab to segment and isolate cranial structures in an automated fashion from a mixed population of mice, including craniofacial mutants, aged 4-12.5 weeks. As a secondary contribution, we demonstrate an application of nearly automated shape analysis, using symmetric diffeomorphic image registration. This approach, which we call diGPA, closely approximates the popular generalized Procrustes analysis (GPA) but negates the collection of anatomical landmarks. We achieve our goals by using the open-source advanced normalization tools (ANT) image quantification library, as well as its associated R library (ANTsR) for statistical image analysis. Finally, we make a plea to investigators to commit to using open imaging standards and software in their labs to the extent possible to increase the potential for data exchange and improve the reproducibility of findings. Future work will incorporate more anatomical detail (such as individual cranial bones, turbinals, dentition, middle ear ossicles) and more diversity into the template.

Multi-template analysis of human perirhinal cortex in brain MRI: Explicitly accounting for anatomical variability

RATIONAL

The human perirhinal cortex (PRC) plays critical roles in episodic and semantic memory and visual perception. The PRC consists of Brodmann areas 35 and 36 (BA35, BA36). In Alzheimer's disease (AD), BA35 is the first cortical site affected by neurofibrillary tangle pathology, which is closely linked to neural injury in AD. Large anatomical variability, manifested in the form of different cortical folding and branching patterns, makes it difficult to segment the PRC in MRI scans. Pathology studies have found that in ~97% of specimens, the PRC falls into one of three discrete anatomical variants. However, current methods for PRC segmentation and morphometry in MRI are based on single-template approaches, which may not be able to accurately model these discrete variants METHODS: A multi-template analysis pipeline that explicitly accounts for anatomical variability is used to automatically label the PRC and measure its thickness in T2-weighted MRI scans. The pipeline uses multi-atlas segmentation to automatically label medial temporal lobe cortices including entorhinal cortex, PRC and the parahippocampal cortex. Pairwise registration between label maps and clustering based on residual dissimilarity after registration are used to construct separate templates for the anatomical variants of the PRC. An optimal path of deformations linking these templates is used to establish correspondences between all the subjects. Experimental evaluation focuses on the ability of single-template and multi-template analyses to detect differences in the thickness of medial temporal lobe cortices between patients with amnestic mild cognitive impairment (aMCI, n=41) and age-matched controls (n=44).

RESULTS

The proposed technique is able to generate templates that recover the three dominant discrete variants of PRC and establish more meaningful correspondences between subjects than a single-template approach. The largest reduction in thickness associated with aMCI, in absolute terms, was found in left BA35 using both regional and summary thickness measures. Further, statistical maps of regional thickness difference between aMCI and controls revealed different patterns for the three anatomical variants.

Arterial spin labeling perfusion predicts longitudinal decline in semantic variant primary progressive aphasia

The objective of the study was to evaluate the prognostic value of regional cerebral blood flow (CBF) measured by arterial spin labeled (ASL) perfusion MRI in patients with semantic variant primary progressive aphasia (svPPA). We acquired pseudo-continuous ASL (pCASL) MRI and whole-brain T1-weighted structural MRI in svPPA patients (N = 13) with cerebrospinal fluid biomarkers consistent with frontotemporal lobar degeneration pathology. Follow-up T1-weighted MRI was available in a subset of patients (N = 8). We performed whole-brain comparisons of partial volume-corrected CBF and cortical thickness between svPPA and controls, and compared baseline and follow-up cortical thickness in regions of significant hypoperfusion and hyperperfusion. Patients with svPPA showed partial volume-corrected hypoperfusion relative to controls in left temporal lobe and insula. svPPA patients also had typical cortical thinning in anterior temporal, insula, and inferior frontal regions at baseline. Volume-corrected hypoperfusion was seen in areas of significant cortical thinning such as the left temporal lobe and insula. Additional regions of hypoperfusion corresponded to areas without cortical thinning. We also observed regions of hyperperfusion, some associated with cortical thinning and others without cortical thinning, including right superior temporal, inferior parietal, and orbitofrontal cortices. Regions of hypoperfusion and hyperperfusion near cortical thinning at baseline had significant longitudinal thinning between baseline and follow-up scans, but perfusion changes in distant areas did not show progressive thinning. Our findings suggest ASL MRI may be sensitive to functional changes not readily apparent in structural MRI, and specific changes in perfusion may be prognostic markers of disease progression in a manner consistent with cell-to-cell spreading pathology.

Neuropsychological Testing Predicts Cerebrospinal Fluid Amyloid-β in Mild Cognitive Impairment

BACKGROUND

Psychometric tests predict conversion of mild cognitive impairment (MCI) to probable Alzheimer's disease (AD). Because the definition of clinical AD relies on those same psychometric tests, the ability of these tests to identify underlying AD pathology remains unclear.

OBJECTIVE

To determine the degree to which psychometric testing predicts molecular evidence of AD amyloid pathology, as indicated by cerebrospinal fluid (CSF) amyloid-β (Aβ)1 - 42, in patients with MCI, as compared to neuroimaging biomarkers.

METHODS

We identified 408 MCI subjects with CSF Aβ levels, psychometric test data, FDG-PET scans, and acceptable volumetric MR scans from the Alzheimer's Disease Neuroimaging Initiative (ADNI). We used psychometric tests and imaging biomarkers in univariate and multivariate models to predict Aβ status.

RESULTS

The 30-min delayed recall score of the Rey Auditory Verbal Learning Test was the best predictor of Aβ status among the psychometric tests, achieving an AUC of 0.67 ± 0.02 and odds ratio of 2.5 ± 0.4. FDG-PET was the best imaging-based biomarker (AUC 0.67 ± 0.03, OR 3.2 ± 1.2), followed by hippocampal volume (AUC 0.64 ± 0.02, OR 2.4 ± 0.3). A multivariate analysis based on the psychometric tests improved on the univariate predictors, achieving an AUC of 0.68 ± 0.03 (OR 3.38 ± 1.2). Adding imaging biomarkers to the multivariate analysis did not improve the AUC.

CONCLUSION

Psychometric tests perform as well as imaging biomarkers to predict presence of molecular markers of AD pathology in MCI patients and should be considered in the determination of the likelihood that MCI is due to AD.

White matter hyperintensities are more highly associated with preclinical Alzheimer's disease than imaging and cognitive markers of neurodegeneration

INTRODUCTION

Cognitive tests and nonamyloid imaging biomarkers do not consistently identify preclinical AD. The objective of this study was to evaluate whether white matter hyperintensity (WMH) volume, a cerebrovascular disease marker, is more associated with preclinical AD than conventional AD biomarkers and cognitive tests.

METHODS

Elderly controls enrolled in the Alzheimer's Disease Neuroimaging Initiative (ADNI, n = 158) underwent florbetapir-PET scans, psychometric testing, neuroimaging with MRI and PET, and APOE genetic testing. Elderly controls the Parkinson's progression markers initiative (PPMI, n = 58) had WMH volume, cerebrospinal fluid (CSF) Aβ1-42, and APOE status measured.

RESULTS

In the ADNI cohort, only WMH volume and APOE ε4 status were associated with cerebral Aβ (standardized β = 0.44 and 1.25, P = .03 and .002). The association between WMH volume and APOE ε4 status with cerebral Aβ (standardized β = 1.12 and 0.26, P = .048 and .045) was confirmed in the PPMI cohort.

DISCUSSION

WMH volume is more highly associated with preclinical AD than other AD biomarkers.

Optimal Symmetric Multimodal Templates and Concatenated Random Forests for Supervised Brain Tumor Segmentation (Simplified) with ANTsR

Segmenting and quantifying gliomas from MRI is an important task for diagnosis, planning intervention, and for tracking tumor changes over time. However, this task is complicated by the lack of prior knowledge concerning tumor location, spatial extent, shape, possible displacement of normal tissue, and intensity signature. To accommodate such complications, we introduce a framework for supervised segmentation based on multiple modality intensity, geometry, and asymmetry feature sets. These features drive a supervised whole-brain and tumor segmentation approach based on random forest-derived probabilities. The asymmetry-related features (based on optimal symmetric multimodal templates) demonstrate excellent discriminative properties within this framework. We also gain performance by generating probability maps from random forest models and using these maps for a refining Markov random field regularized probabilistic segmentation. This strategy allows us to interface the supervised learning capabilities of the random forest model with regularized probabilistic segmentation using the recently developed ANTsR package--a comprehensive statistical and visualization interface between the popular Advanced Normalization Tools (ANTs) and the R statistical project. The reported algorithmic framework was the top-performing entry in the MICCAI 2013 Multimodal Brain Tumor Segmentation challenge. The challenge data were widely varying consisting of both high-grade and low-grade glioma tumor four-modality MRI from five different institutions. Average Dice overlap measures for the final algorithmic assessment were 0.87, 0.78, and 0.74 for "complete", "core", and "enhanced" tumor components, respectively.

The Multimodal Brain Tumor Image Segmentation Benchmark (BRATS)

In this paper we report the set-up and results of the Multimodal Brain Tumor Image Segmentation Benchmark (BRATS) organized in conjunction with the MICCAI 2012 and 2013 conferences. Twenty state-of-the-art tumor segmentation algorithms were applied to a set of 65 multi-contrast MR scans of low- and high-grade glioma patients-manually annotated by up to four raters-and to 65 comparable scans generated using tumor image simulation software. Quantitative evaluations revealed considerable disagreement between the human raters in segmenting various tumor sub-regions (Dice scores in the range 74%-85%), illustrating the difficulty of this task. We found that different algorithms worked best for different sub-regions (reaching performance comparable to human inter-rater variability), but that no single algorithm ranked in the top for all sub-regions simultaneously. Fusing several good algorithms using a hierarchical majority vote yielded segmentations that consistently ranked above all individual algorithms, indicating remaining opportunities for further methodological improvements. The BRATS image data and manual annotations continue to be publicly available through an online evaluation system as an ongoing benchmarking resource.

Assessment of MRI-Based Automated Fetal Cerebral Cortical Folding Measures in Prediction of Gestational Age in the Third Trimester

BACKGROUND AND PURPOSE

Traditional methods of dating a pregnancy based on history or sonographic assessment have a large variation in the third trimester. We aimed to assess the ability of various quantitative measures of brain cortical folding on MR imaging in determining fetal gestational age in the third trimester.

MATERIALS AND METHODS

We evaluated 8 different quantitative cortical folding measures to predict gestational age in 33 healthy fetuses by using T2-weighted fetal MR imaging. We compared the accuracy of the prediction of gestational age by these cortical folding measures with the accuracy of prediction by brain volume measurement and by a previously reported semiquantitative visual scale of brain maturity. Regression models were constructed, and measurement biases and variances were determined via a cross-validation procedure.

RESULTS

The cortical folding measures are accurate in the estimation and prediction of gestational age (mean of the absolute error, 0.43 ± 0.45 weeks) and perform better than (P = .024) brain volume (mean of the absolute error, 0.72 ± 0.61 weeks) or sonography measures (SDs approximately 1.5 weeks, as reported in literature). Prediction accuracy is comparable with that of the semiquantitative visual assessment score (mean, 0.57 ± 0.41 weeks).

CONCLUSIONS

Quantitative cortical folding measures such as global average curvedness can be an accurate and reliable estimator of gestational age and brain maturity for healthy fetuses in the third trimester and have the potential to be an indicator of brain-growth delays for at-risk fetuses and preterm neonates.

The pediatric template of brain perfusion

Magnetic resonance imaging (MRI) captures the dynamics of brain development with multiple modalities that quantify both structure and function. These measurements may yield valuable insights into the neural patterns that mark healthy maturation or that identify early risk for psychiatric disorder. The Pediatric Template of Brain Perfusion (PTBP) is a free and public neuroimaging resource that will help accelerate the understanding of childhood brain development as seen through the lens of multiple modality neuroimaging and in relation to cognitive and environmental factors. The PTBP uses cross-sectional and longitudinal MRI to quantify cortex, white matter, resting state functional connectivity and brain perfusion, as measured by Arterial Spin Labeling (ASL), in 120 children 7-18 years of age. We describe the PTBP and show, as a demonstration of validity, that global summary measurements capture the trajectories that demarcate critical turning points in brain maturation. This novel resource will allow a more detailed understanding of the network-level, structural and functional landmarks that are obtained during normal adolescent brain development.

Automatic clustering and thickness measurement of anatomical variants of the human perirhinal cortex

The entorhinal cortex (ERC) and the perirhinal cortex (PRC) are subregions of the medial temporal lobe (MTL) that play important roles in episodic memory representations, as well as serving as a conduit between other neocortical areas and the hippocampus. They are also the sites where neuronal damage first occurs in Alzheimer's disease (AD). The ability to automatically quantify the volume and thickness of the ERC and PRC is desirable because these localized measures can potentially serve as better imaging biomarkers for AD and other neurodegenerative diseases. However, large anatomical variation in the PRC makes it a challenging area for analysis. In order to address this problem, we propose an automatic segmentation, clustering, and thickness measurement approach that explicitly accounts for anatomical variation. The approach is targeted to highly anisotropic (0.4x0.4x2.0mm3 ) T2-weighted MRI scans that are preferred by many authors for detailed imaging of the MTL, but which pose challenges for segmentation and shape analysis. After automatically labeling MTL substructures using multi-atlas segmentation, our method clusters subjects into groups based on the shape of the PRC, constructs unbiased population templates for each group, and uses the smooth surface representations obtained during template construction to extract regional thickness measurements in the space of each subject. The proposed thickness measures are evaluated in the context of discrimination between patients with Mild Cognitive Impairment (MCI) and normal controls (NC).

Decomposing cerebral blood flow MRI into functional and structural components: a non-local approach based on prediction

We present RIPMMARC (Rotation Invariant Patch-based Multi-Modality Analysis aRChitecture), a flexible and widely applicable method for extracting information unique to a given modality from a multi-modal data set. We use RIPMMARC to improve the interpretation of arterial spin labeling (ASL) perfusion images by removing the component of perfusion that is predicted by the underlying anatomy. Using patch-based, rotation invariant descriptors derived from the anatomical image, we learn a predictive relationship between local neuroanatomical structure and the corresponding perfusion image. This relation allows us to produce an image of perfusion that would be predicted given only the underlying anatomy and a residual image that represents perfusion information that cannot be predicted by anatomical features. Our learned structural features are significantly better at predicting brain perfusion than tissue probability maps, which are the input to standard partial volume correction techniques. Studies in test-retest data show that both the anatomically predicted and residual perfusion signals are highly replicable for a given subject. In a pediatric population, both the raw perfusion and structurally predicted images are tightly linked to age throughout adolescence throughout the brain. Interestingly, the residual perfusion also shows a strong correlation with age in selected regions including the hippocampi (corr = 0.38, p-value <10(-6)), precuneus (corr = -0.44, p < 10(-5)), and combined default mode network regions (corr = -0.45, p < 10(-8)) that is independent of global anatomy-perfusion trends. This finding suggests that there is a regionally heterogeneous pattern of functional specialization that is distinct from that of cortical structural development.

Large-scale evaluation of ANTs and FreeSurfer cortical thickness measurements

Many studies of the human brain have explored the relationship between cortical thickness and cognition, phenotype, or disease. Due to the subjectivity and time requirements in manual measurement of cortical thickness, scientists have relied on robust software tools for automation which facilitate the testing and refinement of neuroscientific hypotheses. The most widely used tool for cortical thickness studies is the publicly available, surface-based FreeSurfer package. Critical to the adoption of such tools is a demonstration of their reproducibility, validity, and the documentation of specific implementations that are robust across large, diverse imaging datasets. To this end, we have developed the automated, volume-based Advanced Normalization Tools (ANTs) cortical thickness pipeline comprising well-vetted components such as SyGN (multivariate template construction), SyN (image registration), N4 (bias correction), Atropos (n-tissue segmentation), and DiReCT (cortical thickness estimation). In this work, we have conducted the largest evaluation of automated cortical thickness measures in publicly available data, comparing FreeSurfer and ANTs measures computed on 1205 images from four open data sets (IXI, MMRR, NKI, and OASIS), with parcellation based on the recently proposed Desikan-Killiany-Tourville (DKT) cortical labeling protocol. We found good scan-rescan repeatability with both FreeSurfer and ANTs measures. Given that such assessments of precision do not necessarily reflect accuracy or an ability to make statistical inferences, we further tested the neurobiological validity of these approaches by evaluating thickness-based prediction of age and gender. ANTs is shown to have a higher predictive performance than FreeSurfer for both of these measures. In promotion of open science, we make all of our scripts, data, and results publicly available which complements the use of open image data sets and the open source availability of the proposed ANTs cortical thickness pipeline.

Subject-specific functional parcellation via prior based eigenanatomy

We present a new framework for prior-constrained sparse decomposition of matrices derived from the neuroimaging data and apply this method to functional network analysis of a clinically relevant population. Matrix decomposition methods are powerful dimensionality reduction tools that have found widespread use in neuroimaging. However, the unconstrained nature of these totally data-driven techniques makes it difficult to interpret the results in a domain where network-specific hypotheses may exist. We propose a novel approach, Prior Based Eigenanatomy (p-Eigen), which seeks to identify a data-driven matrix decomposition but at the same time constrains the individual components by spatial anatomical priors (probabilistic ROIs). We formulate our novel solution in terms of prior-constrained ℓ1 penalized (sparse) principal component analysis. p-Eigen starts with a common functional parcellation for all the subjects and refines it with subject-specific information. This enables modeling of the inter-subject variability in the functional parcel boundaries and allows us to construct subject-specific networks with reduced sensitivity to ROI placement. We show that while still maintaining correspondence across subjects, p-Eigen extracts biologically-relevant and patient-specific functional parcels that facilitate hypothesis-driven network analysis. We construct default mode network (DMN) connectivity graphs using p-Eigen refined ROIs and use them in a classification paradigm. Our results show that the functional connectivity graphs derived from p-Eigen significantly aid classification of mild cognitive impairment (MCI) as well as the prediction of scores in a Delayed Recall memory task when compared to graph metrics derived from 1) standard registration-based seed ROI definitions, 2) totally data-driven ROIs, 3) a model based on standard demographics plus hippocampal volume as covariates, and 4) Ward Clustering based data-driven ROIs. In summary, p-Eigen incarnates a new class of prior-constrained dimensionality reduction tools that may improve our understanding of the relationship between MCI and functional connectivity.

Relating brain anatomy and cognitive ability using a multivariate multimodal framework

Linking structural neuroimaging data from multiple modalities to cognitive performance is an important challenge for cognitive neuroscience. In this study we examined the relationship between verbal fluency performance and neuroanatomy in 54 patients with frontotemporal degeneration (FTD) and 15 age-matched controls, all of whom had T1- and diffusion-weighted imaging. Our goal was to incorporate measures of both gray matter (voxel-based cortical thickness) and white matter (fractional anisotropy) into a single statistical model that relates to behavioral performance. We first used eigenanatomy to define data-driven regions of interest (DD-ROIs) for both gray matter and white matter. Eigenanatomy is a multivariate dimensionality reduction approach that identifies spatially smooth, unsigned principal components that explain the maximal amount of variance across subjects. We then used a statistical model selection procedure to see which of these DD-ROIs best modeled performance on verbal fluency tasks hypothesized to rely on distinct components of a large-scale neural network that support language: category fluency requires a semantic-guided search and is hypothesized to rely primarily on temporal cortices that support lexical-semantic representations; letter-guided fluency requires a strategic mental search and is hypothesized to require executive resources to support a more demanding search process, which depends on prefrontal cortex in addition to temporal network components that support lexical representations. We observed that both types of verbal fluency performance are best described by a network that includes a combination of gray matter and white matter. For category fluency, the identified regions included bilateral temporal cortex and a white matter region including left inferior longitudinal fasciculus and frontal-occipital fasciculus. For letter fluency, a left temporal lobe region was also selected, and also regions of frontal cortex. These results are consistent with our hypothesized neuroanatomical models of language processing and its breakdown in FTD. We conclude that clustering the data with eigenanatomy before performing linear regression is a promising tool for multimodal data analysis.

The Insight ToolKit image registration framework

Publicly available scientific resources help establish evaluation standards, provide a platform for teaching and improve reproducibility. Version 4 of the Insight ToolKit (ITK⁴) seeks to establish new standards in publicly available image registration methodology. ITK⁴ makes several advances in comparison to previous versions of ITK. ITK⁴ supports both multivariate images and objective functions; it also unifies high-dimensional (deformation field) and low-dimensional (affine) transformations with metrics that are reusable across transform types and with composite transforms that allow arbitrary series of geometric mappings to be chained together seamlessly. Metrics and optimizers take advantage of multi-core resources, when available. Furthermore, ITK⁴ reduces the parameter optimization burden via principled heuristics that automatically set scaling across disparate parameter types (rotations vs. translations). A related approach also constrains steps sizes for gradient-based optimizers. The result is that tuning for different metrics and/or image pairs is rarely necessary allowing the researcher to more easily focus on design/comparison of registration strategies. In total, the ITK⁴ contribution is intended as a structure to support reproducible research practices, will provide a more extensive foundation against which to evaluate new work in image registration and also enable application level programmers a broad suite of tools on which to build. Finally, we contextualize this work with a reference registration evaluation study with application to pediatric brain labeling.¹

The power of neuroimaging biomarkers for screening frontotemporal dementia

Frontotemporal dementia (FTD) is a clinically and pathologically heterogeneous neurodegenerative disease that can result from either frontotemporal lobar degeneration (FTLD) or Alzheimer's disease (AD) pathology. It is critical to establish statistically powerful biomarkers that can achieve substantial cost-savings and increase the feasibility of clinical trials. We assessed three broad categories of neuroimaging methods to screen underlying FTLD and AD pathology in a clinical FTD series: global measures (e.g., ventricular volume), anatomical volumes of interest (VOIs) (e.g., hippocampus) using a standard atlas, and data-driven VOIs using Eigenanatomy. We evaluated clinical FTD patients (N = 93) with cerebrospinal fluid, gray matter (GM) magnetic resonance imaging (MRI), and diffusion tensor imaging (DTI) to assess whether they had underlying FTLD or AD pathology. Linear regression was performed to identify the optimal VOIs for each method in a training dataset and then we evaluated classification sensitivity and specificity in an independent test cohort. Power was evaluated by calculating minimum sample sizes required in the test classification analyses for each model. The data-driven VOI analysis using a multimodal combination of GM MRI and DTI achieved the greatest classification accuracy (89% sensitive and 89% specific) and required a lower minimum sample size (N = 26) relative to anatomical VOI and global measures. We conclude that a data-driven VOI approach using Eigenanatomy provides more accurate classification, benefits from increased statistical power in unseen datasets, and therefore provides a robust method for screening underlying pathology in FTD patients for entry into clinical trials.

Single-subject structural networks with closed-form rotation invariant matching mprove power in developmental studies of the cortex

Although much attention has recently been focused on single-subject functional networks, using methods such as resting-state functional MRI, methods for constructing single-subject structural networks are in their infancy. Single-subject cortical networks aim to describe the self-similarity across the cortical structure, possibly signifying convergent developmental pathways. Previous methods for constructing single-subject cortical networks have used patch-based correlations and distance metrics based on curvature and thickness. We present here a method for constructing similarity-based cortical structural networks that utilizes a rotation-invariant representation of structure. The resulting graph metrics are closely linked to age and indicate an increasing degree of closeness throughout development in nearly all brain regions, perhaps corresponding to a more regular structure as the brain matures. The derived graph metrics demonstrate a four-fold increase in power for detecting age as compared to cortical thickness. This proof of concept study indicates that the proposed metric may be useful in identifying biologically relevant cortical patterns.

Explicit B-spline regularization in diffeomorphic image registration

Diffeomorphic mappings are central to image registration due largely to their topological properties and success in providing biologically plausible solutions to deformation and morphological estimation problems. Popular diffeomorphic image registration algorithms include those characterized by time-varying and constant velocity fields, and symmetrical considerations. Prior information in the form of regularization is used to enforce transform plausibility taking the form of physics-based constraints or through some approximation thereof, e.g., Gaussian smoothing of the vector fields [a la Thirion's Demons (Thirion, 1998)]. In the context of the original Demons' framework, the so-called directly manipulated free-form deformation (DMFFD) (Tustison et al., 2009) can be viewed as a smoothing alternative in which explicit regularization is achieved through fast B-spline approximation. This characterization can be used to provide B-spline “flavored” diffeomorphic image registration solutions with several advantages. Implementation is open source and available through the Insight Toolkit and our Advanced Normalization Tools (ANTs) repository. A thorough comparative evaluation with the well-known SyN algorithm (Avants et al., 2008), implemented within the same framework, and its B-spline analog is performed using open labeled brain data and open source evaluation tools.

Genetic and neuroanatomic associations in sporadic frontotemporal lobar degeneration

Genome-wide association studies have identified single nucleotide polymorphisms (SNPs) that are sensitive for tau or TDP-43 pathology in frontotemporal lobar degeneration (FTLD). Neuroimaging analyses have revealed distinct distributions of disease in FTLD patients with genetic mutations. However, genetic influences on neuroanatomic structure in sporadic FTLD have not been assessed. In this report, we use novel multivariate tools, Eigenanatomy, and sparse canonical correlation analysis to identify associations between SNPs and neuroanatomic structure in sporadic FTLD. Magnetic resonance imaging analyses revealed that rs8070723 (MAPT) was associated with gray matter variance in the temporal cortex. Diffusion tensor imaging analyses revealed that rs1768208 (MOBP), rs646776 (near SORT1), and rs5848 (PGRN) were associated with white matter variance in the midbrain and superior longitudinal fasciculus. In an independent autopsy series, we observed that rs8070723 and rs1768208 conferred significant risk of tau pathology relative to TDP-43, and rs646776 conferred increased risk of TDP-43 pathology relative to tau. Identified brain regions and SNPs may help provide an in vivo screen for underlying pathology in FTLD and contribute to our understanding of sporadic FTLD.

Histology-derived volumetric annotation of the human hippocampal subfields in postmortem MRI

Recently, there has been a growing effort to analyze the morphometry of hippocampal subfields using both in vivo and postmortem magnetic resonance imaging (MRI). However, given that boundaries between subregions of the hippocampal formation (HF) are conventionally defined on the basis of microscopic features that often lack discernible signature in MRI, subfield delineation in MRI literature has largely relied on heuristic geometric rules, the validity of which with respect to the underlying anatomy is largely unknown. The development and evaluation of such rules are challenged by the limited availability of data linking MRI appearance to microscopic hippocampal anatomy, particularly in three dimensions (3D). The present paper, for the first time, demonstrates the feasibility of labeling hippocampal subfields in a high resolution volumetric MRI dataset based directly on microscopic features extracted from histology. It uses a combination of computational techniques and manual post-processing to map subfield boundaries from a stack of histology images (obtained with 200μm spacing and 5μm slice thickness; stained using the Kluver-Barrera method) onto a postmortem 9.4Tesla MRI scan of the intact, whole hippocampal formation acquired with 160μm isotropic resolution. The histology reconstruction procedure consists of sequential application of a graph-theoretic slice stacking algorithm that mitigates the effects of distorted slices, followed by iterative affine and diffeomorphic co-registration to postmortem MRI scans of approximately 1cm-thick tissue sub-blocks acquired with 200μm isotropic resolution. These 1cm blocks are subsequently co-registered to the MRI of the whole HF. Reconstruction accuracy is evaluated as the average displacement error between boundaries manually delineated in both the histology and MRI following the sequential stages of reconstruction. The methods presented and evaluated in this single-subject study can potentially be applied to multiple hippocampal tissue samples in order to construct a histologically informed MRI atlas of the hippocampal formation.

Associations between children's socioeconomic status and prefrontal cortical thickness

Childhood socioeconomic status (SES) predicts executive function performance and measures of prefrontal cortical function, but little is known about its anatomical correlates. Structural MRI and demographic data from a sample of 283 healthy children from the NIH MRI Study of Normal Brain Development were used to investigate the relationship between SES and prefrontal cortical thickness. Specifically, we assessed the association between two principal measures of childhood SES, family income and parental education, and gray matter thickness in specific subregions of prefrontal cortex and on the asymmetry of these areas. After correcting for multiple comparisons and controlling for potentially confounding variables, parental education significantly predicted cortical thickness in the right anterior cingulate gyrus and left superior frontal gyrus. These results suggest that brain structure in frontal regions may provide a meaningful link between SES and cognitive function among healthy, typically developing children.

White matter imaging helps dissociate tau from TDP-43 in frontotemporal lobar degeneration

BACKGROUND

Frontotemporal lobar degeneration (FTLD) is most commonly associated with TAR-DNA binding protein (TDP-43) or tau pathology at autopsy, but there are no in vivo biomarkers reliably discriminating between sporadic cases. As disease-modifying treatments emerge, it is critical to accurately identify underlying pathology in living patients so that they can be entered into appropriate etiology-directed clinical trials. Patients with tau inclusions (FTLD-TAU) appear to have relatively greater white matter (WM) disease at autopsy than those patients with TDP-43 (FTLD-TDP). In this paper, we investigate the ability of white matter (WM) imaging to help discriminate between FTLD-TAU and FTLD-TDP during life using diffusion tensor imaging (DTI).

METHODS

Patients with autopsy-confirmed disease or a genetic mutation consistent with FTLD-TDP or FTLD-TAU underwent multimodal T1 volumetric MRI and diffusion weighted imaging scans. We quantified cortical thickness in GM and fractional anisotropy (FA) in WM. We performed Eigenanatomy, a statistically robust dimensionality reduction algorithm, and used leave-one-out cross-validation to predict underlying pathology. Neuropathological assessment of GM and WM disease burden was performed in the autopsy-cases to confirm our findings of an ante-mortem GM and WM dissociation in the neuroimaging cohort.

RESULTS

ROC curve analyses evaluated classification accuracy in individual patients and revealed 96% sensitivity and 100% specificity for WM analyses. FTLD-TAU had significantly more WM degeneration and inclusion severity at autopsy relative to FTLD-TDP.

CONCLUSIONS

These neuroimaging and neuropathological investigations provide converging evidence for greater WM burden associated with FTLD-TAU, and emphasise the role of WM neuroimaging for in vivo discrimination between FTLD-TAU and FTLD-TDP.

A digital atlas of the dog brain

There is a long history and a growing interest in the canine as a subject of study in neuroscience research and in translational neurology. In the last few years, anatomical and functional magnetic resonance imaging (MRI) studies of awake and anesthetized dogs have been reported. Such efforts can be enhanced by a population atlas of canine brain anatomy to implement group analyses. Here we present a canine brain atlas derived as the diffeomorphic average of a population of fifteen mesaticephalic dogs. The atlas includes: 1) A brain template derived from in-vivo, T1-weighted imaging at 1 mm isotropic resolution at 3 Tesla (with and without the soft tissues of the head); 2) A co-registered, high-resolution (0.33 mm isotropic) template created from imaging of ex-vivo brains at 7 Tesla; 3) A surface representation of the gray matter/white matter boundary of the high-resolution atlas (including labeling of gyral and sulcal features). The properties of the atlas are considered in relation to historical nomenclature and the evolutionary taxonomy of the Canini tribe. The atlas is available for download (https://cfn.upenn.edu/aguirre/wiki/public:data_plosone_2012_datta).

Robust Automated Amygdala Segmentation via Multi-Atlas Diffeomorphic Registration

Here, we describe a novel method for volumetric segmentation of the amygdala from MRI images collected from 35 human subjects. This approach is adapted from open-source techniques employed previously with the hippocampus (Suh et al., 2011; Wang et al., 2011a,b). Using multi-atlas segmentation and machine learning-based correction, we were able to produce automated amygdala segments with high Dice (Mean = 0.918 for the left amygdala; 0.916 for the right amygdala) and Jaccard coefficients (Mean = 0.850 for the left; 0.846 for the right) compared to rigorously hand-traced volumes. This automated routine also produced amygdala segments with high intra-class correlations (consistency = 0.830, absolute agreement = 0.819 for the left; consistency = 0.786, absolute agreement = 0.783 for the right) and bivariate (r = 0.831 for the left; r = 0.797 for the right) compared to hand-drawn amygdala. Our results are discussed in relation to other cutting-edge segmentation techniques, as well as commonly available approaches to amygdala segmentation (e.g., Freesurfer). We believe this new technique has broad application to research with large sample sizes for which amygdala quantification might be needed.

Logical circularity in voxel-based analysis: normalization strategy may induce statistical bias

Recent discussions within the neuroimaging community have highlighted the problematic presence of selection bias in experimental design. Although initially centering on the selection of voxels during the course of fMRI studies, we demonstrate how this bias can potentially corrupt voxel-based analyses. For such studies, template-based registration plays a critical role in which a representative template serves as the normalized space for group alignment. A standard approach maps each subject's image to a representative template before performing statistical comparisons between different groups. We analytically demonstrate that in these scenarios the popular sum of squared difference (SSD) intensity metric, implicitly surrogating as a quantification of anatomical alignment, instead explicitly maximizes effect size--an experimental design flaw referred to as "circularity bias." We illustrate how this selection bias varies in strength with the similarity metric used during registration under the hypothesis that while SSD-related metrics, such as Demons, will manifest similar effects, other metrics which are not formulated based on absolute intensity differences will produce less of an effect. Consequently, given the variability in voxel-based analysis outcomes with similarity metric choice, we caution researchers specifically in the use of SSD and SSD-related measures where normalization and statistical analysis involve the same image set. Instead, we advocate a more cautious approach where normalization of the individual subject images to the reference space occurs through corresponding image sets which are independent of statistical testing. Alternatively, one can use similarity terms that are less sensitive to this bias.

Quantitative mouse brain phenotyping based on single and multispectral MR protocols

Sophisticated image analysis methods have been developed for the human brain, but such tools still need to be adapted and optimized for quantitative small animal imaging. We propose a framework for quantitative anatomical phenotyping in mouse models of neurological and psychiatric conditions. The framework encompasses an atlas space, image acquisition protocols, and software tools to register images into this space. We show that a suite of segmentation tools (Avants, Epstein et al., 2008) designed for human neuroimaging can be incorporated into a pipeline for segmenting mouse brain images acquired with multispectral magnetic resonance imaging (MR) protocols. We present a flexible approach for segmenting such hyperimages, optimizing registration, and identifying optimal combinations of image channels for particular structures. Brain imaging with T1, T2* and T2 contrasts yielded accuracy in the range of 83% for hippocampus and caudate putamen (Hc and CPu), but only 54% in white matter tracts, and 44% for the ventricles. The addition of diffusion tensor parameter images improved accuracy for large gray matter structures (by >5%), white matter (10%), and ventricles (15%). The use of Markov random field segmentation further improved overall accuracy in the C57BL/6 strain by 6%; so Dice coefficients for Hc and CPu reached 93%, for white matter 79%, for ventricles 68%, and for substantia nigra 80%. We demonstrate the segmentation pipeline for the widely used C57BL/6 strain, and two test strains (BXD29, APP/TTA). This approach appears promising for characterizing temporal changes in mouse models of human neurological and psychiatric conditions, and may provide anatomical constraints for other preclinical imaging, e.g. fMRI and molecular imaging. This is the first demonstration that multiple MR imaging modalities combined with multivariate segmentation methods lead to significant improvements in anatomical segmentation in the mouse brain.

From label fusion to correspondence fusion: a new approach to unbiased groupwise registration

Label fusion strategies are used in multi-atlas image segmentation approaches to compute a consensus segmentation of an image, given a set of candidate segmentations produced by registering the image to a set of atlases [19, 11, 8]. Effective label fusion strategies, such as local similarity-weighted voting [1, 13] substantially reduce segmentation errors compared to single-atlas segmentation. This paper extends the label fusion idea to the problem of finding correspondences across a set of images. Instead of computing a consensus segmentation, weighted voting is used to estimate a consensus coordinate map between a target image and a reference space. Two variants of the problem are considered: (1) where correspondences between a set of atlases are known and are propagated to the target image; (2) where correspondences are estimated across a set of images without prior knowledge. Evaluation in synthetic data shows that correspondences recovered by fusion methods are more accurate than those based on registration to a population template. In a 2D example in real MRI data, fusion methods result in more consistent mappings between manual segmentations of the hippocampus.

AGREEMENT BETWEEN THE WHITE MATTER CONNECTIVITY BASED ON THE TENSOR-BASED MORPHOMETRY AND THE VOLUMETRIC WHITE MATTER PARCELLATIONS BASED ON DIFFUSION TENSOR IMAGING

We are interested in investigating white matter connectivity using a novel computational framework that does not use diffusion tensor imaging (DTI) but only uses T1-weighted magnetic resonance imaging. The proposed method relies on correlating Jacobian determinants across different voxels based on the tensor-based morphometry (TBM) framework. In this paper, we show agreement between the TBM-based white matter connectivity and the DTI-based white matter atlas. As an application, altered white matter connectivity in a clinical population is determined.

RECONSTRUCTION OF THE HUMAN HIPPOCAMPUS IN 3D FROM HISTOLOGY AND HIGH-RESOLUTION EX-VIVO MRI

In this paper, we present methods for the reconstruction of 3D histological volumes of the human hippocampal formation from histology slices. Inter-slice alignment is guided by a graph-theoretic approach that minimizes the impact of badly distorted slices. The reconstruction is refined by iterative affine and deformable co-registration with a high-resolution MRI of the postmortem tissue sample. We present an evaluation of reconstruction accuracy that is based on measures of similarity between boundaries drawn on both histology and MRI. Our methodology is currently being applied to an MRI atlas of the human hippocampal formation, in which atlas anatomical labels are derived from segmentation of reconstructed histology.

Measuring longitudinal change in the hippocampal formation from in vivo high-resolution T2-weighted MRI

The hippocampal formation (HF) is a brain structure of great interest because of its central role in learning and memory, and its associated vulnerability to several neurological disorders. In vivo oblique coronal T2-weighted MRI with high in-plane resolution (~0.5 mm × 0.5 mm), thick slices (~2.0 mm), and a field of view tailored to imaging the hippocampal formation (denoted HF-MRI in this paper) has been advanced as a useful imaging modality for detailed hippocampal morphometry. Cross-sectional analysis of volume measurements derived from HF-MRI has shown the modality's promise to yield sensitive imaging-based biomarker for neurological disorders such as Alzheimer's disease. However, the utility of this modality for making measurements of longitudinal change has not yet been demonstrated. In this paper, using an unbiased deformation-based morphometry (DBM) pipeline, we examine the suitability of HF-MRI for estimating longitudinal change by comparing atrophy rates measured in the whole hippocampus from this modality with those measured from more common isotropic (~1 mm³) T1-weighted MRI in the same set of individuals, in a cohort of healthy controls and patients with cognitive impairment. While measurements obtained from HF-MRI were largely consistent with those obtained from T1-MRI, HF-MRI yielded slightly larger group effect of greater atrophy rates in patients than in controls. The estimated minimum sample size required for detecting a 25% change in patients' atrophy rate in the hippocampus compared to the control group with a statistical power β=0.8 was N=269. For T1-MRI, the equivalent sample size was N=325. Using a dataset of test-retest scans, we show that the measurements were free of additive bias. We also demonstrate that these results were not a confound of certain methodological choices made in the DBM pipeline to address the challenges of making longitudinal measurements from HF-MRI, using a region of interest (ROI) around the HF to globally align serial images, followed by slice-by-slice deformable registration to measure local volume change. Additionally, we present a preliminary study of atrophy rate measurements within hippocampal subfields using HF-MRI. Cross-sectional differences in atrophy rates were detected in several subfields.

Ventilation-based segmentation of the lungs using hyperpolarized (3)He MRI

PURPOSE

To develop an automated segmentation method to differentiate the ventilated lung volume on (3) He magnetic resonance imaging (MRI).

MATERIALS AND METHODS

Computational processing (CP) for each subject consisted of the following three essential steps: 1) inhomogeneity bias correction, 2) whole lung segmentation, and 3) subdivision of the lung segmentation into regions of similar ventilation. Evaluation consisted of two comparative analyses: i) comparison of the number of defects scored by two human readers in 43 subjects, and ii) simultaneous truth and performance level estimation (STAPLE) in 18 subjects in which the ventilation defects were manually segmented by four human readers.

RESULTS

There was excellent correlation between the number of ventilation defects tabulated by CP and reader #1 (intraclass correlation coefficient [ICC] = 0.86), CP and reader #2 (ICC = 0.85), and between the two readers (ICC = 0.97). The STAPLE results from the second analysis yielded the following sensitivity/specificity numbers: CP (0.898/0.905), radiologist #1 (0.743/0.897), radiologist #2 (0.501/0.985), radiologist #3 (0.898/0.848), and the first author (0.600/0.984).

CONCLUSION

We developed and evaluated an automated method for quantifying the ventilated lung volume on (3) He MRI. The findings strongly indicate that our proposed algorithmic processing may be a reliable, automatic method for quantitating ventilation defects.

Topological well-composedness and glamorous glue: a digital gluing algorithm for topologically constrained front propagation

We propose a new approach to front propagation algorithms based on a topological variant of well-composedness which contrasts with previous methods based on simple point detection. This provides for a theoretical justification, based on the digital Jordan separation theorem, for digitally "gluing" evolved well-composed objects separated by well-composed curves or surfaces. Additionally, our framework can be extended to more relaxed topologically constrained algorithms based on multisimple points. For both methods this framework has the additional benefit of obviating the requirement for both a user-specified connectivity and a topologically-consistent marching cubes/squares algorithm in meshing the resulting segmentation.

Estimation of perfusion and arterial transit time in myocardium using free-breathing myocardial arterial spin labeling with navigator-echo

Arterial spin labeling (ASL) provides noninvasive measurement of tissue blood flow, but sensitivity to motion has limited its application to imaging of myocardial blood flow. Although different cardiac phases can be synchronized using electrocardiography triggering, breath holding is generally required to minimize effects of respiratory motion during ASL scanning, which may be challenging in clinical populations. Here a free-breathing myocardial ASL technique with the potential for reliable clinical application is presented, by combining ASL with a navigator-gated, electrocardiography-triggered TrueFISP readout sequence. Dynamic myocardial perfusion signals were measured at multiple delay times that allowed simultaneous fitting of myocardial blood flow and arterial transit time. With the assist of a nonrigid motion correction program, the estimated mean myocardial blood flow was 1.00 ± 0.55 mL/g/min with a mean transit time of ∼ 400 msec. The intraclass correlation coefficient of repeated scans was 0.89 with a mean within subject coefficient of variation of 22%. Perfusion response during mild to moderate stress was further measured. The capability for noninvasive, free-breathing assessment of myocardial blood flow using ASL may offer an alternative approach to first-pass perfusion MRI for clinical evaluation of patients with coronary artery disease.

STRUCTURAL CONNECTIVITY VIA THE TENSOR-BASED MORPHOMETRY

The tensor-based morphometry (TBM) has been widely used in characterizing tissue volume difference between populations at voxel level. We present a novel computational framework for investigating the white matter connectivity using TBM. Unlike other diffusion tensor imaging (DTI) based white matter connectivity studies, we do not use DTI but only T1-weighted magnetic resonance imaging (MRI). To construct brain network graphs, we have developed a new data-driven approach called the ε-neighbor method that does not need any predetermined parcellation. The proposed pipeline is applied in detecting the topological alteration of the white matter connectivity in maltreated children.