Generalized tensor-based morphometry of HIV/AIDS using multivariate statistics on deformation tensors

Surface-based morphometry of the corpus callosum in young children of ages 1-5

The corpus callosum (CC) is a large white matter fiber bundle in the brain and is involved in various cognitive, sensory, and motor processes. While implicated in various developmental and psychiatric disorders, much is yet to be uncovered about the normal development of this structure, especially in young children. Additionally, while sexual dimorphism has been reported in prior literature, observations have not necessarily been consistent. In this study, we use morphometric measures including surface tensor-based morphometry (TBM) to investigate local changes in the shape of the CC in children between the ages of 12 and 60 months, in intervals of 12 months. We also analyze sex differences in each of these age groups. We observed larger significant clusters in the earlier ages between 12 v 24 m and between 48 v 60 m and localized differences in the anterior region of the body of the CC. Sex differences were most pronounced in the 12 m group. This study adds to the growing literature of work aiming to understand the developing brain and emphasizes the utility of surface TBM as a useful tool for analyzing regional differences in neuroanatomical morphometry.

Classification of Alzheimer's disease based on hippocampal multivariate morphometry statistics

BACKGROUND

Alzheimer's disease (AD) is a neurodegenerative disease characterized by progressive cognitive decline, and mild cognitive impairment (MCI) is associated with a high risk of developing AD. Hippocampal morphometry analysis is believed to be the most robust magnetic resonance imaging (MRI) markers for AD and MCI. Multivariate morphometry statistics (MMS), a quantitative method of surface deformations analysis, is confirmed to have strong statistical power for evaluating hippocampus.

AIMS

We aimed to test whether surface deformation features in hippocampus can be employed for early classification of AD, MCI, and healthy controls (HC).

METHODS

We first explored the differences in hippocampus surface deformation among these three groups by using MMS analysis. Additionally, the hippocampal MMS features of selective patches and support vector machine (SVM) were used for the binary classification and triple classification.

RESULTS

By the results, we identified significant hippocampal deformation among the three groups, especially in hippocampal CA1. In addition, the binary classification of AD/HC, MCI/HC, AD/MCI showed good performances, and area under curve (AUC) of triple-classification model achieved 0.85. Finally, positive correlations were found between the hippocampus MMS features and cognitive performances.

CONCLUSIONS

The study revealed significant hippocampal deformation among AD, MCI, and HC. Additionally, we confirmed that hippocampal MMS can be used as a sensitive imaging biomarker for the early diagnosis of AD at the individual level.

Spatiotemporal Atlas of the Fetal Brain Depicts Cortical Developmental Gradient

The fetal brains experience rapid and complex development in utero during the second and third trimesters. In utero MRI of the fetal brain in this period enables us to quantify normal fetal brain development in the spatiotemporal domain. In this study, we established a high-quality spatiotemporal atlas between 23 and 38 weeks gestational age (GA) from 90 healthy Chinese human fetuses of both sexes using a pairwise and groupwise registration pipeline. We quantified the fetal cortical morphology indices and characterized their spatiotemporal developmental pattern. The cortical thickness exhibited a biphasic pattern that first increased and then decreased; the curvature fitted well into the Gompertz growth model; sulcal depth increased linearly, while surface area expanded exponentially. The cortical thickness and curvature trajectories consistently pointed to a characteristic time point around GA of 31 weeks. The characteristic GA and growth rate obtained from individual cortical regions suggested a central-to-peripheral developmental gradient, with the earliest development in the parietal lobe, and we also observed a superior-to-inferior gradient within the temporal lobe. These findings may be linked to biophysical events, such as dendritic arborization and thalamocortical fibers ingrowth. The proposed atlas was also compared with an existing fetal atlas from a white/mixed population. Finally, we examined the structural asymmetry of the fetal brains and found extensive asymmetry that dynamically changed with development. The current study depicted a comprehensive profile of fetal cortical development, and the established atlas could be used as a normative reference for neurodevelopmental and diagnostic purposes, especially in the Chinese population.SIGNIFICANCE STATEMENT We generated a high-quality 4D spatiotemporal atlas of the normal fetal brain development from 23 to 38 gestational weeks in a Chinese population and characterized the spatiotemporal developmental pattern of cortical morphology. According to the cortical development trajectories, the fetal cerebral cortex development follows a central-to-peripheral developmental gradient that may be related to the underlying cellular events. The majority of cortical regions already exhibit significant asymmetry during the fetal period.

Bladder surface dose modeling in prostate cancer radiotherapy: An analysis of motion-induced variations and the cumulative dose across the treatment

PURPOSE

To introduce a novel surface-based dose mapping method to improve quantitative bladder dosimetric assessment in prostate cancer (PC) radiotherapy.

METHODS

Based on the planning and daily pre and postfraction MRIs of 12 PC patients, bladder surface models (SMs) were generated on manually delineated contours and regionally aligned via surface-based registration. Subsequently, bladder surface dose models (SDMs) were created using face-wise dose sampling. To determine the bladder intrafractional and interfractional motion and dose variation, we performed a pose analysis between pre and postfraction bladder SMs, as well as surface mapping for fractional SMs. Discrepancies between the received dose, accumulated from daily SDMs, and the planned dose were then assessed on the corresponding SDMs. Complementary to the surface dose mapping, dose surface histogram (DSH)-based comparisons were also performed.

RESULTS

The intrafraction pose analysis revealed a significant (p < 0.05) bladder expansion, as well as an anterior/superior drift during the treatment. The intrafraction motion substantially altered dose to mid-bladder body, but not the bladder surface areas distal to or contiguous with the target. A similar pattern of dose variations was also detected by interfraction comparisons. With surface registration to the common SM, the cumulative bladder dose significantly differs from the planned dose. The discrepancy is evident in the mid-posterior range that corresponds to a mid- to high-dose region. The received DSH significantly differs from the planned DSH after permutation correction (p = 0.0122), while the overall surface-based comparison after multiple comparison correction is nonsignificant (p = 0.0800).

CONCLUSIONS

We developed a novel surface-based intra and interdose mapping framework applied to a unique daily MR dataset for image-guided radiotherapy. The framework identified significant intrafraction bladder positional changes, localized the intra and interfraction variations, and quantified planned versus received dose differences on the bladder surface. The result indicates the importance of adopting the motion-integrated bladder SDM for bladder dose management.

An Image Registration-Based Morphing Technique for Generating Subject-Specific Brain Finite Element Models

Finite element (FE) models of the brain are crucial for investigating the mechanisms of traumatic brain injury (TBI). However, FE brain models are often limited to a single neuroanatomy because the manual development of subject-specific models is time consuming. The objective of this study was to develop a pipeline to automatically generate subject-specific FE brain models using previously developed nonlinear image registration techniques, preserving both external and internal neuroanatomical characteristics. To verify the morphing-induced mesh distortions did not influence the brain deformation response, strain distributions predicted using the morphed model were compared to those from manually created voxel models of the same subject. Morphed and voxel models were generated for 44 subjects ranging in age, and simulated using head kinematics from a football concussion case. For each subject, brain strain distributions predicted by each model type were consistent, and differences in strain prediction was less than 4% between model type. This automated technique, taking approximately 2 h to generate a subject-specific model, will facilitate interdisciplinary research between the biomechanics and neuroimaging fields and could enable future use of biomechanical models in the clinical setting as a tool for improving diagnosis.

Applying surface-based morphometry to study ventricular abnormalities of cognitively unimpaired subjects prior to clinically significant memory decline

•

A completely automated surface-based ventricular morphometry system.

•

Generate a whole connected 3D ventricular shape model.

•

Test-retest the system in two independent CU subject cohorts.

•

Subregional ventricular abnormalities prior to clinically memory decline.

Ventricular volume (VV) is a widely used structural magnetic resonance imaging (MRI) biomarker in Alzheimer’s disease (AD) research. Abnormal enlargements of VV can be detected before clinically significant memory decline. However, VV does not pinpoint the details of subregional ventricular expansions. Here we introduce a ventricular morphometry analysis system (VMAS) that generates a whole connected 3D ventricular shape model and encodes a great deal of ventricular surface deformation information that is inaccessible by VV. VMAS contains an automated segmentation approach and surface-based multivariate morphometry statistics. We applied VMAS to two independent datasets of cognitively unimpaired (CU) groups. To our knowledge, it is the first work to detect ventricular abnormalities that distinguish normal aging subjects from those who imminently progress to clinically significant memory decline. Significant bilateral ventricular morphometric differences were first shown in 38 members of the Arizona APOE cohort, which included 18 CU participants subsequently progressing to the clinically significant memory decline within 2 years after baseline visits (progressors), and 20 matched CU participants with at least 4 years of post-baseline cognitive stability (non-progressors). VMAS also detected significant differences in bilateral ventricular morphometry in 44 Alzheimer’s Disease Neuroimaging Initiative (ADNI) subjects (18 CU progressors vs. 26 CU non-progressors) with the same inclusion criterion. Experimental results demonstrated that the ventricular anterior horn regions were affected bilaterally in CU progressors, and more so on the left. VMAS may track disease progression at subregional levels and measure the effects of pharmacological intervention at a preclinical stage.

Brain structural correlates of trajectories to cognitive impairment in men with and without HIV disease

There are distinct trajectories to cognitive impairment among participants in the Multicenter AIDS Cohort Study (MACS). Here we analyzed the relationship between regional brain volumes and the individual trajectories to impairment in a subsample (n = 302) of the cohort. 302 (167 HIV-infected; mean age = 55.7 yrs.; mean education: 16.2 yrs.) of the men enrolled in the MACS MRI study contributed data to this analysis. We used voxel-based morphometry (VBM) to segment the brain images to analyze gray and white matter volume at the voxel-level. A Mixed Membership Trajectory Model had previously identified three distinct profiles, and each study participant had a membership weight for each of these three trajectories. We estimated VBM model parameters for 100 imputations, manually performed the post-hoc contrasts, and pooled the results. We examined the associations between brain volume at the voxel level and the MMTM membership weights for two profiles: one considered "unhealthy" and the other considered "Premature aging." The unhealthy profile was linked to the volume of the posterior cingulate gyrus/precuneus, the inferior frontal cortex, and the insula, whereas the premature aging profile was independently associated with the integrity of a portion of the precuneus. Trajectories to cognitive impairment are the result, in part, of atrophy in cortical regions linked to normal and pathological aging. These data suggest the possibility of predicting cognitive morbidity based on patterns of CNS atrophy.

Morphological changes in subregions of hippocampus and amygdala in major depressive disorder patients

Despite many neuroimaging studies in the past years, the neuroanatomical substrates of major depressive disorder (MDD) subcortical structures are still not well understood. Since hippocampus and amygdala are the two vital subcortical structures that most susceptible to MDD, finding the evidence of morphological changes in their subregions may bring some new insights for MDD research. Combining structural magnetic resonance imaging (MRI) with novel morphometry analysis methods, we recruited 25 MDD patients and 28 healthy controls (HC), and investigated their volume and morphological differences in hippocampus and amygdala. Relative to volumetric method, our methods detected more significant global morphological atrophies (p<0.05). More precisely, subiculum and cornu ammonis (CA) 1 subregions of bilateral hippocampus, lateral (LA) and basolateral ventromedial (BLVM) of left amygdala and LA, BLVM, central (CE), amygdalostriatal transition area (ASTR), anterior cortical (ACO) and anterior amygdaloid area (AAA) of right amygdala were demonstrated prone to atrophy. Correlation analyses between each subject's surface eigenvalues and Hamilton Depression Scale (HAMD) were then performed. Correlation results showed that atrophy areas in hippocampus and amygdala have slight tendencies of expanding into other subregions with the development of MDD. Finally, we performed group morphometric analysis and drew the atrophy and expansion areas between MDD-Medicated group (only 19 medicated subjects in MDD group were included) and HC group, found some preliminary evidence about subregional morphological resilience of hippocampus and amygdala. These findings revealed new pathophysiologic patterns in the subregions of hippocampus and amygdala, which can help with subsequent smaller-scale MDD research.

Applying surface-based hippocampal morphometry to study APOE-E4 allele dose effects in cognitively unimpaired subjects

Apolipoprotein E (APOE) e4 is the major genetic risk factor for late-onset Alzheimer's disease (AD). The dose-dependent impact of this allele on hippocampal volumes has been documented, but its influence on general hippocampal morphology in cognitively unimpaired individuals is still elusive. Capitalizing on the study of a large number of cognitively unimpaired late middle aged and older adults with two, one and no APOE-e4 alleles, the current study aims to characterize the ability of our automated surface-based hippocampal morphometry algorithm to distinguish between these three levels of genetic risk for AD and demonstrate its superiority to a commonly used hippocampal volume measurement. We examined the APOE-e4 dose effect on cross-sectional hippocampal morphology analysis in a magnetic resonance imaging (MRI) database of 117 cognitively unimpaired subjects aged between 50 and 85 years (mean = 57.4, SD = 6.3), including 36 heterozygotes (e3/e4), 37 homozygotes (e4/e4) and 44 non-carriers (e3/e3). The proposed automated framework includes hippocampal surface segmentation and reconstruction, higher-order hippocampal surface correspondence computation, and hippocampal surface deformation analysis with multivariate statistics. In our experiments, the surface-based method identified APOE-e4 dose effects on the left hippocampal morphology. Compared to the widely-used hippocampal volume measure, our hippocampal morphometry statistics showed greater statistical power by distinguishing cognitively unimpaired subjects with two, one, and no APOE-e4 alleles. Our findings mirrored previous studies showing that APOE-e4 has a dose effect on the acceleration of brain structure deformities. The results indicated that the proposed surface-based hippocampal morphometry measure is a potential preclinical AD imaging biomarker for cognitively unimpaired individuals.

•

Applied surface-based hippocampal morphometry on cognitively unimpaired subjects.

•

Our study identified APOE-e4 dose effects on cognitively unimpaired subjects.

•

Surface-based hippocampal morphometry outperformed the hippocampal volume measure.

•

Surface-based hippocampal morphometry may be a potential preclinical AD biomarker.

MIDAS: Regionally linear multivariate discriminative statistical mapping

Statistical parametric maps formed via voxel-wise mass-univariate tests, such as the general linear model, are commonly used to test hypotheses about regionally specific effects in neuroimaging cross-sectional studies where each subject is represented by a single image. Despite being informative, these techniques remain limited as they ignore multivariate relationships in the data. Most importantly, the commonly employed local Gaussian smoothing, which is important for accounting for registration errors and making the data follow Gaussian distributions, is usually chosen in an ad hoc fashion. Thus, it is often suboptimal for the task of detecting group differences and correlations with non-imaging variables. Information mapping techniques, such as searchlight, which use pattern classifiers to exploit multivariate information and obtain more powerful statistical maps, have become increasingly popular in recent years. However, existing methods may lead to important interpretation errors in practice (i.e., misidentifying a cluster as informative, or failing to detect truly informative voxels), while often being computationally expensive. To address these issues, we introduce a novel efficient multivariate statistical framework for cross-sectional studies, termed MIDAS, seeking highly sensitive and specific voxel-wise brain maps, while leveraging the power of regional discriminant analysis. In MIDAS, locally linear discriminative learning is applied to estimate the pattern that best discriminates between two groups, or predicts a variable of interest. This pattern is equivalent to local filtering by an optimal kernel whose coefficients are the weights of the linear discriminant. By composing information from all neighborhoods that contain a given voxel, MIDAS produces a statistic that collectively reflects the contribution of the voxel to the regional classifiers as well as the discriminative power of the classifiers. Critically, MIDAS efficiently assesses the statistical significance of the derived statistic by analytically approximating its null distribution without the need for computationally expensive permutation tests. The proposed framework was extensively validated using simulated atrophy in structural magnetic resonance imaging (MRI) and further tested using data from a task-based functional MRI study as well as a structural MRI study of cognitive performance. The performance of the proposed framework was evaluated against standard voxel-wise general linear models and other information mapping methods. The experimental results showed that MIDAS achieves relatively higher sensitivity and specificity in detecting group differences. Together, our results demonstrate the potential of the proposed approach to efficiently map effects of interest in both structural and functional data.

Imaging studies of the HIV-infected brain

Human immunodeficiency virus (HIV) enters the brain early after infecting humans and may remain in the central nervous system despite successful antiretroviral treatment. Many neuroimaging techniques were used to study HIV+ patients with or without opportunistic infections. These techniques assessed abnormalities in brain structures (using computed tomography, structural magnetic resonance imaging (MRI), diffusion MRI) and function (using functional MRI at rest or during a task, and perfusion MRI with or without a contrast agent). In addition, single-photon emission computed tomography with various tracers (e.g., thallium-201, Tc99-HMPAO) and positron emission tomography with various agents (e.g., [18F]-dexoyglucose, [11C]-PiB, and [11C]-TSPO tracers), were applied to study opportunistic infections or HIV-associated neurocognitive disorders. Neuroimaging provides diagnoses and biomarkers to quantitate the severity of brain injury or to monitor treatment effects, and may yield insights into the pathophysiology of HIV infection. As the majority of antiretroviral-stable HIV+ patients are living longer, age-related comorbid disorders (e.g., additional neuroinflammation, cerebrovascular disorders, or other dementias) will need to be considered. Other highly prevalent conditions, such as substance use disorders, psychiatric illnesses, and the long-term effects of combined antiretroviral therapy, all may lead to additional brain injury. Neuroimaging studies could provide knowledge regarding how these comorbid conditions impact the HIV-infected brain. Lastly, specific molecular imaging agents may be needed to assess the central nervous system viral reservoir.

A novel registration-based methodology for prediction of trabecular bone fabric from clinical QCT: A comprehensive analysis

Osteoporosis leads to hip fractures in aging populations and is diagnosed by modern medical imaging techniques such as quantitative computed tomography (QCT). Hip fracture sites involve trabecular bone, whose strength is determined by volume fraction and orientation, known as fabric. However, bone fabric cannot be reliably assessed in clinical QCT images of proximal femur. Accordingly, we propose a novel registration-based estimation of bone fabric designed to preserve tensor properties of bone fabric and to map bone fabric by a global and local decomposition of the gradient of a non-rigid image registration transformation. Furthermore, no comprehensive analysis on the critical components of this methodology has been previously conducted. Hence, the aim of this work was to identify the best registration-based strategy to assign bone fabric to the QCT image of a patient’s proximal femur. The normalized correlation coefficient and curvature-based regularization were used for image-based registration and the Frobenius norm of the stretch tensor of the local gradient was selected to quantify the distance among the proximal femora in the population. Based on this distance, closest, farthest and mean femora with a distinction of sex were chosen as alternative atlases to evaluate their influence on bone fabric prediction. Second, we analyzed different tensor mapping schemes for bone fabric prediction: identity, rotation-only, rotation and stretch tensor. Third, we investigated the use of a population average fabric atlas. A leave one out (LOO) evaluation study was performed with a dual QCT and HR-pQCT database of 36 pairs of human femora. The quality of the fabric prediction was assessed with three metrics, the tensor norm (TN) error, the degree of anisotropy (DA) error and the angular deviation of the principal tensor direction (PTD). The closest femur atlas (CTP) with a full rotation (CR) for fabric mapping delivered the best results with a TN error of 7.3 ± 0.9%, a DA error of 6.6 ± 1.3% and a PTD error of 25 ± 2°. The closest to the population mean femur atlas (MTP) using the same mapping scheme yielded only slightly higher errors than CTP for substantially less computing efforts. The population average fabric atlas yielded substantially higher errors than the MTP with the CR mapping scheme. Accounting for sex did not bring any significant improvements. The identified fabric mapping methodology will be exploited in patient-specific QCT-based finite element analysis of the proximal femur to improve the prediction of hip fracture risk.

Contribution to speech development of the right anterior putamen revealed with multivariate tensor-based morphometry

In our previous study1, we suggested that the difference between tensor-based metrics in the anterior part of the right putamen between 21 and 18 months age groups associated with speech development during this ages. Here we used a correlational analysis between verbal scores and determinant of the Jacobian matrix to confirm our hypothesis. Significant correlations in anterior part of the right putamen between verbal scores and surface metric were revealed in the 18 and 21 age groups.

Feature selective temporal prediction of Alzheimer's disease progression using hippocampus surface morphometry

INTRODUCTION

Prediction of Alzheimer's disease (AD) progression based on baseline measures allows us to understand disease progression and has implications in decisions concerning treatment strategy. To this end, we combine a predictive multi-task machine learning method (cFSGL) with a novel MR-based multivariate morphometric surface map of the hippocampus (mTBM) to predict future cognitive scores of patients.

METHODS

Previous work has shown that a multi-task learning framework that performs prediction of all future time points simultaneously (cFSGL) can be used to encode both sparsity as well as temporal smoothness. The authors showed that this method is able to predict cognitive outcomes of ADNI subjects using FreeSurfer-based baseline MRI features, MMSE score demographic information and ApoE status. Whilst volumetric information may hold generalized information on brain status, we hypothesized that hippocampus specific information may be more useful in predictive modeling of AD. To this end, we applied a multivariate tensor-based parametric surface analysis method (mTBM) to extract features from the hippocampal surfaces.

RESULTS

We combined mTBM features with traditional surface features such as middle axis distance, the Jacobian determinant as well as 2 of the Jacobian principal eigenvalues to yield 7 normalized hippocampal surface maps of 300 points each. By combining these 7 × 300 = 2100 features together with the previous ~350 features, we illustrate how this type of sparsifying method can be applied to an entire surface map of the hippocampus that yields a feature space that is 2 orders of magnitude larger than what was previously attempted.

CONCLUSIONS

By combining the power of the cFSGL multi-task machine learning framework with the addition of AD sensitive mTBM feature maps of the hippocampus surface, we are able to improve the predictive performance of ADAS cognitive scores 6, 12, 24, 36 and 48 months from baseline.

Riemannian Nonlinear Mixed Effects Models: Analyzing Longitudinal Deformations in Neuroimaging

Statistical machine learning models that operate on manifold-valued data are being extensively studied in vision, motivated by applications in activity recognition, feature tracking and medical imaging. While non-parametric methods have been relatively well studied in the literature, efficient formulations for parametric models (which may offer benefits in small sample size regimes) have only emerged recently. So far, manifold-valued regression models (such as geodesic regression) are restricted to the analysis of cross-sectional data, i.e., the so-called "fixed effects" in statistics. But in most "longitudinal analysis" (e.g., when a participant provides multiple measurements, over time) the application of fixed effects models is problematic. In an effort to answer this need, this paper generalizes non-linear mixed effects model to the regime where the response variable is manifold-valued, i.e., f : Rd → ℳ. We derive the underlying model and estimation schemes and demonstrate the immediate benefits such a model can provide - both for group level and individual level analysis - on longitudinal brain imaging data. The direct consequence of our results is that longitudinal analysis of manifold-valued measurements (especially, the symmetric positive definite manifold) can be conducted in a computationally tractable manner.

Topographies of Cortical and Subcortical Volume Loss in HIV and Aging in the cART Era

OBJECTIVES

Studies of HIV-associated brain atrophy often focus on a priori brain regions of interest, which can introduce bias. A data-driven, minimally biased approach was used to analyze changes in brain volumetrics associated with HIV and their relationship to aging, viral factors, combination antiretroviral therapy (cART), and gender, and smoking.

DESIGN

A cross-sectional study of 51 HIV-uninfected (HIV-) and 146 HIV-infected (HIV+) participants.

METHODS

Structural MRI of participants was analyzed using principal component analysis (PCA) to reduce dimensionality and determine topographies of volumetric changes. Neuropsychological (NP) assessment was examined using global and domain-specific scores. The effects of HIV disease factors (eg, viral load, CD4, etc.) on brain volumes and neuropsychological were investigated using penalized regression (LASSO).

RESULTS

Two components of interest were visualized using principal component analysis. An aging effect predominated for both components. The first component, a cortically weighted topography, accounted for a majority of variance across participants (43.5% of variance) and showed independent effects of HIV and smoking. A secondary, subcortically weighted topography (4.6%) showed HIV-status accentuated age-related volume loss. In HIV+ patients, the cortical topography correlated with global neuropsychological scores and nadir CD4, whereas subcortical volume loss was associated with recent viral load.

CONCLUSIONS

Cortical regions showed the most prominent volumetric changes because of aging and HIV. Within HIV+ participants, cortical volumes were associated with immune history, whereas subcortical changes correlated with current immune function. Cognitive function was primarily associated with cortical volume changes. Observed volumetric changes in chronic HIV+ patients may reflect both past infection history and current viral status.

Ventricular shape and relative position abnormalities in preterm neonates

Recent neuroimaging findings have highlighted the impact of premature birth on subcortical development and morphological changes in the deep grey nuclei and ventricular system. To help characterize subcortical microstructural changes in preterm neonates, we recently implemented a multivariate tensor-based method (mTBM). This method allows to precisely measure local surface deformation of brain structures in infants. Here, we investigated ventricular abnormalities and their spatial relationships with surrounding subcortical structures in preterm neonates. We performed regional group comparisons on the surface morphometry and relative position of the lateral ventricles between 19 full-term and 17 preterm born neonates at term-equivalent age. Furthermore, a relative pose analysis was used to detect individual differences in translation, rotation, and scale of a given brain structure with respect to an average. Our mTBM results revealed broad areas of alterations on the frontal horn and body of the left ventricle, and narrower areas of differences on the temporal horn of the right ventricle. A significant shift in the rotation of the left ventricle was also found in preterm neonates. Furthermore, we located significant correlations between morphology and pose parameters of the lateral ventricles and that of the putamen and thalamus. These results show that regional abnormalities on the surface and pose of the ventricles are also associated with alterations on the putamen and thalamus. The complementarity of the information provided by the surface and pose analysis may help to identify abnormal white and grey matter growth, hinting toward a pattern of neural and cellular dysmaturation.

Complement C5aR1 Signaling Promotes Polarization and Proliferation of Embryonic Neural Progenitor Cells through PKCζ

The complement system, typically associated with innate immunity, is emerging as a key controller of nonimmune systems including in development, with recent studies linking complement mutations with neurodevelopmental disease. A key effector of the complement response is the activation fragment C5a, which, through its receptor C5aR1, is a potent driver of inflammation. Surprisingly, C5aR1 is also expressed during early mammalian embryogenesis; however, no clearly defined function is ascribed to C5aR1 in development. Here we demonstrate polarized expression of C5aR1 on the apical surface of mouse embryonic neural progenitor cells in vivo and on human embryonic stem cell-derived neural progenitors. We also show that signaling of endogenous C5a during mouse embryogenesis drives proliferation of neural progenitor cells within the ventricular zone and is required for normal brain histogenesis. C5aR1 signaling in neural progenitors was dependent on atypical protein kinase C ζ, a mediator of stem cell polarity, with C5aR1 inhibition reducing proliferation and symmetric division of apical neural progenitors in human and mouse models. C5aR1 signaling was shown to promote the maintenance of cell polarity, with exogenous C5a increasing the retention of polarized rosette architecture in human neural progenitors after physical or chemical disruption. Transient inhibition of C5aR1 during neurogenesis in developing mice led to behavioral abnormalities in both sexes and MRI-detected brain microstructural alterations, in studied males, demonstrating a requirement of C5aR1 signaling for appropriate brain development. This study thus identifies a functional role for C5a-C5aR1 signaling in mammalian neurogenesis and provides mechanistic insight into recently identified complement gene mutations and brain disorders.SIGNIFICANCE STATEMENT The complement system, traditionally known as a controller of innate immunity, now stands as a multifaceted signaling family with a broad range of physiological actions. These include roles in the brain, where complement activation is associated with diseases, including epilepsy and schizophrenia. This study has explored complement regulation of neurogenesis, identifying a novel relationship between the complement activation peptide C5a and the neural progenitor proliferation underpinning formation of the mammalian brain. C5a was identified as a regulator of cell polarity, with inhibition of C5a receptors during embryogenesis leading to abnormal brain development and behavioral deficits. This work demonstrates mechanisms through which dysregulation of complement causes developmental disease and highlights the potential risk of complement inhibition for therapeutic purposes in pregnancy.

Mapping the basal ganglia alterations in children chronically exposed to manganese

Chronic manganese (Mn) exposure is associated with neuromotor and neurocognitive deficits, but the exact mechanism of Mn neurotoxicity is still unclear. With the advent of magnetic resonance imaging (MRI), in-vivo analysis of brain structures has become possible. Among different sub-cortical structures, the basal ganglia (BG) has been investigated as a putative anatomical biomarker in MR-based studies of Mn toxicity. However, previous investigations have yielded inconsistent results in terms of regional MR signal intensity changes. These discrepancies may be due to the subtlety of brain alterations caused by Mn toxicity, coupled to analysis techniques that lack the requisite detection power. Here, based on brain MRI, we apply a 3D surface-based morphometry method on 3 bilateral basal ganglia structures in school-age children chronically exposed to Mn through drinking water to investigate the effect of Mn exposure on brain anatomy. Our method successfully pinpointed significant enlargement of many areas of the basal ganglia structures, preferentially affecting the putamen. Moreover, these areas showed significant correlations with fine motor performance, indicating a possible link between altered basal ganglia neurodevelopment and declined motor performance in high Mn exposed children.

Altered structural brain changes and neurocognitive performance in pediatric HIV

Pediatric HIV patients often suffer with neurodevelopmental delay and subsequently cognitive impairment. While tissue injury in cortical and subcortical regions in the brain of adult HIV patients has been well reported there is sparse knowledge about these changes in perinatally HIV infected pediatric patients. We analyzed cortical thickness, subcortical volume, structural connectivity, and neurocognitive functions in pediatric HIV patients and compared with those of pediatric healthy controls. With informed consent, 34 perinatally infected pediatric HIV patients and 32 age and gender matched pediatric healthy controls underwent neurocognitive assessment and brain magnetic resonance imaging (MRI) on a 3 T clinical scanner. Altered cortical thickness, subcortical volumes, and abnormal neuropsychological test scores were observed in pediatric HIV patients. The structural network connectivity analysis depicted lower connection strengths, lower clustering coefficients, and higher path length in pediatric HIV patients than healthy controls. The network betweenness and network hubs in cortico-limbic regions were distorted in pediatric HIV patients. The findings suggest that altered cortical and subcortical structures and regional brain connectivity in pediatric HIV patients may contribute to deficits in their neurocognitive functions. Further, longitudinal studies are required for better understanding of the effect of HIV pathogenesis on brain structural changes throughout the brain development process under standard ART treatment.

•

Structural brain MRI and cognitive assessments were performed in pediatric HIV.

•

Pediatric HIV showed altered cortical thickness and subcortical volumes.

•

Disrupted structural connectivity was observed in pediatric HIV.

•

Altered brain structures and connectivity contribute to deficits in neurocognition.

A T1 and DTI fused 3D corpus callosum analysis in MCI subjects with high and low cardiovascular risk profile

Understanding the extent to which vascular disease and its risk factors are associated with prodromal dementia, notably Alzheimer's disease (AD), may enhance predictive accuracy as well as guide early interventions. One promising avenue to determine this relationship consists of looking for reliable and sensitive in-vivo imaging methods capable of characterizing the subtle brain alterations before the clinical manifestations. However, little is known from the imaging perspective about how risk factors such as vascular disease influence AD progression. Here, for the first time, we apply an innovative T1 and DTI fusion analysis of 3D corpus callosum (CC) on mild cognitive impairment (MCI) populations with different levels of vascular profile, aiming to de-couple the vascular factor in the prodromal AD stage. Our new fusion method successfully increases the detection power for differentiating MCI subjects with high from low vascular risk profiles, as well as from healthy controls. MCI subjects with high and low vascular risk profiles showed differed alteration patterns in the anterior CC, which may help to elucidate the inter-wired relationship between MCI and vascular risk factors.

Biomechanical Analysis of Normal Brain Development during the First Year of Life Using Finite Strain Theory

The first year of life is the most critical time period for structural and functional development of the human brain. Combining longitudinal MR imaging and finite strain theory, this study aimed to provide new insights into normal brain development through a biomechanical framework. Thirty-three normal infants were longitudinally imaged using MRI from 2 weeks to 1 year of age. Voxel-wise Jacobian determinant was estimated to elucidate volumetric changes while Lagrange strains (both normal and shear strains) were measured to reveal directional growth information every 3 months during the first year of life. Directional normal strain maps revealed that, during the first 6 months, the growth pattern of gray matter is anisotropic and spatially inhomogeneous with higher left-right stretch around the temporal lobe and interhemispheric fissure, anterior-posterior stretch in the frontal and occipital lobes, and superior-inferior stretch in right inferior occipital and right inferior temporal gyri. In contrast, anterior lateral ventricles and insula showed an isotropic stretch pattern. Volumetric and directional growth rates were linearly decreased with age for most of the cortical regions. Our results revealed anisotropic and inhomogeneous brain growth patterns of the human brain during the first year of life using longitudinal MRI and a biomechanical framework.

Multivariate surface-based analysis of corpus callosum in patients with sickle cell disease

Sickle cell disease (SCD) is a genetic hematological disease in which the hemoglobin molecule in red blood cells is abnormal. It is closely associated with many symptoms, including pain, anemia, chest syndrome and neurocognitive impairment. One of the most debilitating symptoms is elevated risk for cerebro-vascular accidents. The corpus callosum (CC), as the largest and most prominent white matter (WM) structure in the brain, can reflect the chronic cerebrovascular damage resulting from silent strokes or infarctions in asymptomatic SCD patients. While a lot of studies have reported WM alterations in this cohort, little is known about the shape deformation of the CC. Here we perform the first surface morphometry analysis of the CC in SCD patients using four different shape metrics on T1-weighted magnetic resonance images. We detect regional surface morphological differences in the CC between 11 patients and 10 healthy control subjects. Differences are located in the genu, posterior midbody and splenium, potentially casting light on the anatomical substrates underlying neuropsychological test differences between the SCD and control groups.

Lognormal Distributions and Geometric Averages of Symmetric Positive Definite Matrices

This article gives a formal definition of a lognormal family of probability distributions on the set of symmetric positive definite (SPD) matrices, seen as a matrix-variate extension of the univariate lognormal family of distributions. Two forms of this distribution are obtained as the large sample limiting distribution via the central limit theorem of two types of geometric averages of i.i.d. SPD matrices: the log-Euclidean average and the canonical geometric average. These averages correspond to two different geometries imposed on the set of SPD matrices. The limiting distributions of these averages are used to provide large-sample confidence regions and two-sample tests for the corresponding population means. The methods are illustrated on a voxelwise analysis of diffusion tensor imaging data, permitting a comparison between the various average types from the point of view of their sampling variability.

Putamen Development in Children 12 to 21 Months Old

We studied the developmental trajectory of the putamen in 13-21 months old children using multivariate surface tensor-based morphometry. Our results indicate surface changes between 12 and 15 months' age groups in the middle superior part the left putamen. The growth of the left putamen at earlier ages slows down after 15 months. The most important surface changes were detected in the right putamen between 18 and 21 months and were located in the anterior part of the structure. Our results demonstrate the heterochronic growth of the right and left putamen related to different functional subregions within putamen. Our results are compatible with previous studies devoted to total putamen volume changes during normal development.

Influence of APOE Genotype on Hippocampal Atrophy over Time - An N=1925 Surface-Based ADNI Study

The apolipoprotein E (APOE) e4 genotype is a powerful risk factor for late-onset Alzheimer’s disease (AD). In the Alzheimer’s Disease Neuroimaging Initiative (ADNI) cohort, we previously reported significant baseline structural differences in APOE e4 carriers relative to non-carriers, involving the left hippocampus more than the right—a difference more pronounced in e4 homozygotes than heterozygotes. We now examine the longitudinal effects of APOE genotype on hippocampal morphometry at 6-, 12- and 24-months, in the ADNI cohort. We employed a new automated surface registration system based on conformal geometry and tensor-based morphometry. Among different hippocampal surfaces, we computed high-order correspondences, using a novel inverse-consistent surface-based fluid registration method and multivariate statistics consisting of multivariate tensor-based morphometry (mTBM) and radial distance. At each time point, using Hotelling’s T² test, we found significant morphological deformation in APOE e4 carriers relative to non-carriers in the full cohort as well as in the non-demented (pooled MCI and control) subjects at each follow-up interval. In the complete ADNI cohort, we found greater atrophy of the left hippocampus than the right, and this asymmetry was more pronounced in e4 homozygotes than heterozygotes. These findings, combined with our earlier investigations, demonstrate an e4 dose effect on accelerated hippocampal atrophy, and support the enrichment of prevention trial cohorts with e4 carriers.

Deformation-based atrophy computation by surface propagation and its application to Alzheimer's disease

Obtaining regional volume changes from a deformation field is more precise when using simplex counting (SC) compared with Jacobian integration (JI) due to the numerics involved in the latter. Although SC has been proposed before, numerical properties underpinning the method and a thorough evaluation of the method against JI is missing in the literature. The contributions of this paper are: (a) we propose surface propagation (SP)-a simplification to SC that significantly reduces its computational complexity; (b) we will derive the orders of approximation of SP which can also be extended to SC. In the experiments, we will begin by empirically showing that SP is indeed nearly identical to SC, and that both methods are more stable than JI in presence of moderate to large deformation noise. Since SC and SP are identical, we consider SP as a representative of both the methods for a practical evaluation against JI. In a real application on Alzheimer's disease neuroimaging initiative data, we show the following: (a) SP produces whole brain and medial temporal lobe atrophy numbers that are significantly better than JI at separating between normal controls and Alzheimer's disease patients; (b) SP produces disease group atrophy differences comparable to or better than those obtained using FreeSurfer, demonstrating the validity of the obtained clinical results. Finally, in a reproducibility study, we show that the voxel-wise application of SP yields significantly lower variance when compared to JI.

Impact of Early and Late Visual Deprivation on the Structure of the Corpus Callosum: A Study Combining Thickness Profile with Surface Tensor-Based Morphometry

Blindness represents a unique model to study how visual experience may shape the development of brain organization. Exploring how the structure of the corpus callosum (CC) reorganizes ensuing visual deprivation is of particular interest due to its important functional implication in vision (e.g., via the splenium of the CC). Moreover, comparing early versus late visually deprived individuals has the potential to unravel the existence of a sensitive period for reshaping the CC structure. Here, we develop a novel framework to capture a complete set of shape differences in the CC between congenitally blind (CB), late blind (LB) and sighted control (SC) groups. The CCs were manually segmented from T1-weighted brain MRI and modeled by 3D tetrahedral meshes. We statistically compared the combination of local area and thickness at each point between subject groups. Differences in area are found using surface tensor-based morphometry; thickness is estimated by tracing the streamlines in the volumetric harmonic field. Group differences were assessed on this combined measure using Hotelling's T(2) test. Interestingly, we observed that the total callosal volume did not differ between the groups. However, our fine-grained analysis reveals significant differences mostly localized around the splenium areas between both blind groups and the sighted group (general effects of blindness) and, importantly, specific dissimilarities between the LB and CB groups, illustrating the existence of a sensitive period for reorganization. The new multivariate statistics also gave better effect sizes for detecting morphometric differences, relative to other statistics. They may boost statistical power for CC morphometric analyses.

Structural Brain Changes in Early-Onset Alzheimer's Disease Subjects Using the LONI Pipeline Environment

BACKGROUND AND PURPOSE

This study investigates 36 subjects aged 55-65 from the Alzheimer's Disease Neuroimaging Initiative (ADNI) database to expand our knowledge of early-onset (EO) Alzheimer's Disease (EO-AD) using neuroimaging biomarkers.

METHODS

Nine of the subjects had EO-AD, and 27 had EO mild cognitive impairment (EO-MCI). The structural ADNI data were parcellated using BrainParser, and the 15 most discriminating neuroimaging markers between the two cohorts were extracted using the Global Shape Analysis (GSA) Pipeline workflow. Then the Local Shape Analysis (LSA) Pipeline workflow was used to conduct local (per-vertex) post-hoc statistical analyses of the shape differences based on the participants' diagnoses (EO-MCI+EO-AD). Tensor-based Morphometry (TBM) and multivariate regression models were used to identify the significance of the structural brain differences based on the participants' diagnoses.

RESULTS

The significant between-group regional differences using GSA were found in 15 neuroimaging markers. The results of the LSA analysis workflow were based on the subject diagnosis, age, years of education, apolipoprotein E (ε4), Mini-Mental State Examination, visiting times, and logical memory as regressors. All the variables had significant effects on the regional shape measures. Some of these effects survived the false discovery rate (FDR) correction. Similarly, the TBM analysis showed significant effects on the Jacobian displacement vector fields, but these effects were reduced after FDR correction.

CONCLUSIONS

These results may explain some of the differences between EO-AD and EO-MCI, and some of the characteristics of the EO cognitive impairment subjects.

Novel Neuroimaging Methods to Understand How HIV Affects the Brain

In much of the developed world, the HIV epidemic has largely been controlled by antiretroviral treatment. Even so, there is growing concern that HIV-infected individuals may be at risk for accelerated brain aging and a range of cognitive impairments. What promotes or resists these changes is largely unknown. There is also interest in discovering factors that promote resilience to HIV and combat its adverse effects in children. Here, we review recent developments in brain imaging that reveal how the virus affects the brain. We relate these brain changes to changes in blood markers, cognitive function, and other patient outcomes or symptoms, such as apathy or neuropathic pain. We focus on new and emerging techniques, including new variants of brain MRI. Diffusion tensor imaging, for example, can map the brain's structural connections, while fMRI can uncover functional connections. Finally, we suggest how large-scale global research alliances, such as ENIGMA, may resolve controversies over effects where evidence is now lacking. These efforts pool scans from tens of thousands of individuals and offer a source of power not previously imaginable for brain imaging studies.

Shape analysis of the human brain: a brief survey

The survey outlines and compares popular computational techniques for quantitative description of shapes of major structural parts of the human brain, including medial axis and skeletal analysis, geodesic distances, Procrustes analysis, deformable models, spherical harmonics, and deformation morphometry, as well as other less widely used techniques. Their advantages, drawbacks, and emerging trends, as well as results of applications, in particular, for computer-aided diagnostics, are discussed.

Diffusion MRI abnormalities detection with orientation distribution functions: a multiple sclerosis longitudinal study

We propose a new algorithm for the voxelwise analysis of orientation distribution functions between one image and a group of reference images. It relies on a generic framework for the comparison of diffusion probabilities on the sphere, sampled from the underlying models. We demonstrate that this method, combined to dimensionality reduction through a principal component analysis, allows for more robust detection of lesions on simulated data when compared to classical tensor-based analysis. We then demonstrate the efficiency of this pipeline on the longitudinal comparison of multiple sclerosis patients at an early stage of the disease: right after their first clinically isolated syndrome (CIS) and three months later. We demonstrate the predictive value of ODF-based scores for the early detection of lesions that will appear or heal.

A T1 and DTI fused 3D Corpus Callosum analysis in pre- vs. post-season contact sports players

Sports related traumatic brain injury (TBI) is a worldwide public health issue, and damage to the corpus callosum (CC) has been considered as an important indicator of TBI. However, contact sports players suffer repeated hits to the head during the course of a season even in the absence of diagnosed concussion, and less is known about their effect on callosal anatomy. In addition, T1-weighted and diffusion tensor brain magnetic resonance images (DTI) have been analyzed separately, but a joint analysis of both types of data may increase statistical power and give a more complete understanding of anatomical correlates of subclinical concussions in these athletes. Here, for the first time, we fuse T1 surface-based morphometry and a new DTI analysis on 3D surface representations of the CCs into a single statistical analysis on these subjects. Our new combined method successfully increases detection power in detecting differences between pre- vs. post-season contact sports players. Alterations are found in the ventral genu, isthmus, and splenium of CC. Our findings may inform future health assessments in contact sports players. The new method here is also the first truly multimodal diffusion and T1-weighted analysis of the CC in TBI, and may be useful to detect anatomical changes in the corpus callosum in other multimodal datasets.

Regional gray matter atrophy in relapsing remitting multiple sclerosis: baseline analysis of multi-center data

Regional gray matter (GM) atrophy in multiple sclerosis (MS) at disease onset and its temporal variation can provide objective information regarding disease evolution. An automated pipeline for estimating atrophy of various GM structures was developed using tensor based morphometry (TBM) and implemented on a multi-center sub-cohort of 1008 relapsing remitting MS (RRMS) patients enrolled in a Phase 3 clinical trial. Four hundred age and gender matched healthy controls were used for comparison. Using the analysis of covariance, atrophy differences between MS patients and healthy controls were assessed on a voxel-by-voxel analysis. Regional GM atrophy was observed in a number of deep GM structures that included thalamus, caudate nucleus, putamen, and cortical GM regions. General linear regression analysis was performed to analyze the effects of age, gender, and scanner field strength, and imaging sequence on the regional atrophy. Correlations between regional GM volumes and expanded disability status scale (EDSS) scores, disease duration (DD), T2 lesion load (T2 LL), T1 lesion load (T1 LL), and normalized cerebrospinal fluid (nCSF) were analyzed using Pearson׳s correlation coefficient. Thalamic atrophy observed in MS patients compared to healthy controls remained consistent within subgroups based on gender and scanner field strength. Weak correlations between thalamic volume and EDSS (r=-0.133; p<0.001) and DD (r=-0.098; p=0.003) were observed. Of all the structures, thalamic volume moderately correlated with T2 LL (r=-0.492; P-value<0.001), T1 LL (r=-0.473; P-value<0.001) and nCSF (r=-0.367; P-value<0.001).

Thalamic alterations in preterm neonates and their relation to ventral striatum disturbances revealed by a combined shape and pose analysis

Finding the neuroanatomical correlates of prematurity is vital to understanding which structures are affected, and to designing efficient prevention and treatment strategies. Converging results reveal that thalamic abnormalities are important indicators of prematurity. However, little is known about the localization of the abnormalities within the subnuclei of the thalamus, or on the association of altered thalamic development with other deep gray matter disturbances. Here, we aim to investigate the effect of prematurity on the thalamus and the putamen in the neonatal brain, and further investigate the associated abnormalities between these two structures. Using brain structural magnetic resonance imaging, we perform a novel combined shape and pose analysis of the thalamus and putamen between 17 preterm (41.12 ± 5.08 weeks) and 19 term-born (45.51 ± 5.40 weeks) neonates at term equivalent age. We also perform a set of correlation analyses between the thalamus and the putamen, based on the surface and pose results. We locate significant alterations on specific surface regions such as the anterior and ventral anterior (VA) thalamic nuclei, and significant relative pose changes of the left thalamus and the right putamen. In addition, we detect significant association between the thalamus and the putamen for both surface and pose parameters. The regions that are significantly associated include the VA, and the anterior and inferior putamen. We detect statistically significant surface deformations and pose changes on the thalamus and putamen, and for the first time, demonstrate the feasibility of using relative pose parameters as indicators for prematurity in neonates. Our methods show that regional abnormalities of the thalamus are associated with alterations of the putamen, possibly due to disturbed development of shared pre-frontal connectivity. More specifically, the significantly correlated regions in these two structures point to frontal-subcortical pathways including the dorsolateral prefrontal-subcortical circuit, the lateral orbitofrontal-subcortical circuit, the motor circuit, and the oculomotor circuit. These findings reveal new insight into potential subcortical structural covariates for poor neurodevelopmental outcomes in the preterm population.

A new fast accurate nonlinear medical image registration program including surface preserving regularization

Recently inexpensive graphical processing units (GPUs) have become established as a viable alternative to traditional CPUs for many medical image processing applications. GPUs offer the potential of very significant improvements in performance at low cost and with low power consumption. One way in which GPU programs differ from traditional CPU programs is that increasingly elaborate calculations per voxel may not impact of the overall processing time because memory accesses can dominate execution time. This paper presents a new GPU based elastic image registration program named Ezys. The Ezys image registration algorithm belongs to the wide class of diffeomorphic demons but uses surface preserving image smoothing and regularization filters designed for a GPU that would be computationally expensive on a CPU. We describe the methods used in Ezys and present results from two important neuroscience applications. Firstly inter-subject registration for transfer of anatomical labels and secondly longitudinal intra-subject registration to quantify atrophy in individual subjects. Both experiments showed that Ezys registration compares favorably with other popular elastic image registration programs. We believe Ezys is a useful tool for neuroscience and other applications, and also demonstrates the value of developing of novel image processing filters specifically designed for GPUs.

Plasma cytokine levels are related to brain volumes in HIV-infected individuals

HIV-infected individuals frequently exhibit brain dysfunction despite antiretroviral treatment. The neuropathological mechanisms underlying these abnormalities remain unclear, pointing to the importance of identifying biomarkers sensitive to brain dysfunction. We examined 74 medically stable HIV-infected individuals using T1-weighted MRI. Volumes of the cortical grey matter (GM), white matter (WM), caudate, putamen, globus pallidus, thalamus, hippocampus, amygdala, and ventricles were derived using automated parcellation. A panel of plasma cytokines was measured using multiplexed bead array immunoassay. A model selection algorithm was used to select the combination of clinical and cytokine markers that best predicted each brain volumetric measure in a series of linear regression models. Higher CD4 nadir, shorter HIV infection duration, and antiretroviral treatment were significantly related to higher volumes of the putamen, thalamus, hippocampus, and WM. Older age was related to lower volumes in most brain regions and higher ventricular volume. Higher IFN-γ, MCP-1, and TNF-α were related to higher volumes of the putamen, pallidum, amygdala, GM, and WM. Higher IL-1β, IL-6, IL-16, IL-18, IP-10, MIP-1β, and SDF-1α were related to lower volumes of the putamen, pallidum, thalamus, hippocampus, amygdala, GM, and WM; and higher ventricular volume. The current findings provide evidence linking smaller brain volumes to HIV disease history, antiretroviral treatment, and advanced age. Cytokine markers, especially IL-6 and IL-16, showed robust association with brain volumes even after accounting for other clinical variables, demonstrating their utility in examining the mechanisms of HIV-associated brain abnormalities.

Genetic influence of apolipoprotein E4 genotype on hippocampal morphometry: An N = 725 surface-based Alzheimer's disease neuroimaging initiative study

The apolipoprotein E (APOE) e4 allele is the most prevalent genetic risk factor for Alzheimer's disease (AD). Hippocampal volumes are generally smaller in AD patients carrying the e4 allele compared to e4 noncarriers. Here we examined the effect of APOE e4 on hippocampal morphometry in a large imaging database-the Alzheimer's Disease Neuroimaging Initiative (ADNI). We automatically segmented and constructed hippocampal surfaces from the baseline MR images of 725 subjects with known APOE genotype information including 167 with AD, 354 with mild cognitive impairment (MCI), and 204 normal controls. High-order correspondences between hippocampal surfaces were enforced across subjects with a novel inverse consistent surface fluid registration method. Multivariate statistics consisting of multivariate tensor-based morphometry (mTBM) and radial distance were computed for surface deformation analysis. Using Hotelling's T(2) test, we found significant morphological deformation in APOE e4 carriers relative to noncarriers in the entire cohort as well as in the nondemented (pooled MCI and control) subjects, affecting the left hippocampus more than the right, and this effect was more pronounced in e4 homozygotes than heterozygotes. Our findings are consistent with previous studies that showed e4 carriers exhibit accelerated hippocampal atrophy; we extend these findings to a novel measure of hippocampal morphometry. Hippocampal morphometry has significant potential as an imaging biomarker of early stage AD.

Common variants in psychiatric risk genes predict brain structure at birth

Studies in adolescents and adults have demonstrated that polymorphisms in putative psychiatric risk genes are associated with differences in brain structure, but cannot address when in development these relationships arise. To determine if common genetic variants in disrupted-in-schizophrenia-1 (DISC1; rs821616 and rs6675281), catechol-O-methyltransferase (COMT; rs4680), neuregulin 1 (NRG1; rs35753505 and rs6994992), apolipoprotein E (APOE; ε3ε4 vs. ε3ε3), estrogen receptor alpha (ESR1; rs9340799 and rs2234693), brain-derived neurotrophic factor (BDNF; rs6265), and glutamate decarboxylase 1 (GAD1; rs2270335) are associated with individual differences in brain tissue volumes in neonates, we applied both automated region-of-interest volumetry and tensor-based morphometry to a sample of 272 neonates who had received high-resolution magnetic resonance imaging scans. ESR1 (rs9340799) predicted intracranial volume. Local variation in gray matter (GM) volume was significantly associated with polymorphisms in DISC1 (rs821616), COMT, NRG1, APOE, ESR1 (rs9340799), and BDNF. No associations were identified for DISC1 (rs6675281), ESR1 (rs2234693), or GAD1. Of note, neonates homozygous for the DISC1 (rs821616) serine allele exhibited numerous large clusters of reduced GM in the frontal lobes, and neonates homozygous for the COMT valine allele exhibited reduced GM in the temporal cortex and hippocampus, mirroring findings in adults. The results highlight the importance of prenatal brain development in mediating psychiatric risk.

Automatic clustering and population analysis of white matter tracts using maximum density paths

We introduce a framework for population analysis of white matter tracts based on diffusion-weighted images of the brain. The framework enables extraction of fibers from high angular resolution diffusion images (HARDI); clustering of the fibers based partly on prior knowledge from an atlas; representation of the fiber bundles compactly using a path following points of highest density (maximum density path; MDP); and registration of these paths together using geodesic curve matching to find local correspondences across a population. We demonstrate our method on 4-Tesla HARDI scans from 565 young adults to compute localized statistics across 50 white matter tracts based on fractional anisotropy (FA). Experimental results show increased sensitivity in the determination of genetic influences on principal fiber tracts compared to the tract-based spatial statistics (TBSS) method. Our results show that the MDP representation reveals important parts of the white matter structure and considerably reduces the dimensionality over comparable fiber matching approaches.

Joint statistics on cardiac shape and fiber architecture

Cardiac fiber architecture plays an important role in electrophysiological and mechanical functions of the heart. Yet, its inter-subject variability and more particularly, its relationship to the shape of the myocardium, is not fully understood. In this paper, we extend the statistical analysis of cardiac fiber architecture beyond its description with a fixed average geometry. We study the co-variation of fiber architecture with either shape or strain-based information by exploring their principal modes of joint variations. We apply our general framework to a dataset of 8 ex vivo canine hearts, and find that strain-based information appears to correlate best with the fiber architecture. Furthermore, compared to current approaches that warp an average atlas to the patient geometry, our preliminary results show that joint statistics improves fiber synthesis from shape by 8.0%, with cases up to 25.9%. Our experiments also reveal evidence on a possible relation between architectural variability and myocardial thickness.

Mapping white matter integrity in elderly people with HIV

People with HIV are living longer as combination antiretroviral therapy (cART) becomes more widely available. However, even when plasma viral load is reduced to untraceable levels, chronic HIV infection is associated with neurological deficits and brain atrophy beyond that of normal aging. HIV is often marked by cortical and subcortical atrophy, but the integrity of the brain's white matter (WM) pathways also progressively declines. Few studies focus on older cohorts where normal aging may be compounded with HIV infection to influence deficit patterns. In this relatively large diffusion tensor imaging (DTI) study, we investigated abnormalities in WM fiber integrity in 56 HIV+ adults with access to cART (mean age: 63.9 ± 3.7 years), compared to 31 matched healthy controls (65.4 ± 2.2 years). Statistical 3D maps revealed the independent effects of HIV diagnosis and age on fractional anisotropy (FA) and diffusivity, but we did not find any evidence for an age by diagnosis interaction in our current sample. Compared to healthy controls, HIV patients showed pervasive FA decreases and diffusivity increases throughout WM. We also assessed neuropsychological (NP) summary z-score associations. In both patients and controls, fiber integrity measures were associated with NP summary scores. The greatest differences were detected in the corpus callosum and in the projection fibers of the corona radiata. These deficits are consistent with published NP deficits and cortical atrophy patterns in elderly people with HIV.

Voxelwise single-subject analysis of imaging metabolic response to therapy in neuro-oncology

F-18-Fluorodeoxyglucose positron emission tomography (FDG-PET) has been used to evaluate the metabolic response of metastatic brain tumors to treatment by comparing their tumor glucose metabolism before and after treatment. The standard analysis based on regions-of-interest has the advantage of simplicity. However, it is by definition restricted to those regions and is subject to observer variability. In addition, the observed changes in tumor metabolism are often confounded by normal changes in the tissue background, which can be heterogenous. We propose an analysis pipeline for automatically detecting the change at each voxel in the entire brain of a single subject, while adjusting for changes in the background. The complete analysis includes image registration, segmentation, a hierarchical model for background adjustment and voxelwise statistical comparisons. We demonstrate the method's ability to identify areas of tumor response and/or progression in two subjects enrolled in a clinical trial using FDG-PET to evaluate lapatinib for the treatment of brain metastases in breast cancer patients.

The adaptive FEM elastic model for medical image registration

This paper proposes an adaptive mesh refinement strategy for the finite element method (FEM) based elastic registration model. The signature matrix for mesh refinement takes into account the regional intensity variance and the local deformation displacement. The regional intensity variance reflects detailed information for improving registration accuracy and the deformation displacement fine-tunes the mesh refinement for a more efficient algorithm. The gradient flows of two different similarity metrics, the sum of the squared difference and the spatially encoded mutual information for the mono-modal and multi-modal registrations, are used to derive external forces to drive the model to the equilibrium state. We compared our approach to three other models: (1) the conventional multi-resolution FEM registration algorithm; (2) the FEM elastic method that uses variation information for mesh refinement; and (3) the robust block matching based registration. Comparisons among different methods in a dataset with 20 CT image pairs upon artificial deformation demonstrate that our registration method achieved significant improvement in accuracies. Experimental results in another dataset of 40 real medical image pairs for both mono-modal and multi-modal registrations also show that our model outperforms the other three models in its accuracy.

Structural patterns of the proximal femur in relation to age and hip fracture risk in women

Fractures of the proximal femur are the most devastating outcome of osteoporosis. It is generally understood that age-related changes in hip structure confer increased risk, but there have been few explicit comparisons of such changes in healthy subjects to those with hip fracture. In this study, we used quantitative computed tomography and tensor-based morphometry (TBM) to identify three-dimensional internal structural patterns of the proximal femur associated with age and with incident hip fracture. A population-based cohort of 349 women representing a broad age range (21-97years) was included in this study, along with a cohort of 222 older women (mean age 79±7years) with (n=74) and without (n=148) incident hip fracture. Images were spatially normalized to a standardized space, and age- and fracture-specific morphometric features were identified based on statistical maps of shape features described as local changes of bone volume. Morphometric features were visualized as maps of local contractions and expansions, and significance was displayed as Student's t-test statistical maps. Significant age-related changes included local expansions of regions low in volumetric bone mineral density (vBMD) and local contractions of regions high in vBMD. Some significant fracture-related features resembled an accentuated aging process, including local expansion of the superior aspect of the trabecular bone compartment in the femoral neck, with contraction of the adjoining cortical bone. However, other features were observed only in the comparison of hip fracture subjects with age-matched controls including focal contractions of the cortical bone at the superior aspect of the femoral neck, the lateral cortical bone just inferior to the greater trochanter, and the anterior intertrochanteric region. Results of this study support the idea that the spatial distribution of morphometric features is relevant to age-related changes in bone and independent to fracture risk. In women, the identification by TBM of fracture-specific morphometric alterations of the proximal femur, in conjunction with vBMD and clinical risk factors, may improve hip fracture prediction.

Surface fluid registration of conformal representation: application to detect disease burden and genetic influence on hippocampus

In this paper, we develop a new automated surface registration system based on surface conformal parameterization by holomorphic 1-forms, inverse consistent surface fluid registration, and multivariate tensor-based morphometry (mTBM). First, we conformally map a surface onto a planar rectangle space with holomorphic 1-forms. Second, we compute surface conformal representation by combining its local conformal factor and mean curvature and linearly scale the dynamic range of the conformal representation to form the feature image of the surface. Third, we align the feature image with a chosen template image via the fluid image registration algorithm, which has been extended into the curvilinear coordinates to adjust for the distortion introduced by surface parameterization. The inverse consistent image registration algorithm is also incorporated in the system to jointly estimate the forward and inverse transformations between the study and template images. This alignment induces a corresponding deformation on the surface. We tested the system on Alzheimer's Disease Neuroimaging Initiative (ADNI) baseline dataset to study AD symptoms on hippocampus. In our system, by modeling a hippocampus as a 3D parametric surface, we nonlinearly registered each surface with a selected template surface. Then we used mTBM to analyze the morphometry difference between diagnostic groups. Experimental results show that the new system has better performance than two publicly available subcortical surface registration tools: FIRST and SPHARM. We also analyzed the genetic influence of the Apolipoprotein E[element of]4 allele (ApoE4), which is considered as the most prevalent risk factor for AD. Our work successfully detected statistically significant difference between ApoE4 carriers and non-carriers in both patients of mild cognitive impairment (MCI) and healthy control subjects. The results show evidence that the ApoE genotype may be associated with accelerated brain atrophy so that our work provides a new MRI analysis tool that may help presymptomatic AD research.

Disrupted cerebral metabolite levels and lower nadir CD4 + counts are linked to brain volume deficits in 210 HIV-infected patients on stable treatment

Cognitive impairment and brain injury are common in people with HIV/AIDS, even when viral replication is effectively suppressed with combined antiretroviral therapies (cART). Metabolic and structural abnormalities may promote cognitive decline, but we know little about how these measures relate in people on stable cART. Here we used tensor-based morphometry (TBM) to reveal the 3D profile of regional brain volume variations in 210 HIV + patients scanned with whole-brain MRI at 1.5 T (mean age: 48.6 ± 8.4 years; all receiving cART). We identified brain regions where the degree of atrophy was related to HIV clinical measures and cerebral metabolite levels assessed with magnetic resonance spectroscopy (MRS). Regional brain volume reduction was linked to lower nadir CD4 + count, with a 1–2% white matter volume reduction for each 25-point reduction in nadir CD4 +. Even so, brain volume measured by TBM showed no detectable association with current CD4 + count, AIDS Dementia Complex (ADC) stage, HIV RNA load in plasma or cerebrospinal fluid (CSF), duration of HIV infection, antiretroviral CNS penetration-effectiveness (CPE) scores, or years on cART, after controlling for demographic factors, and for multiple comparisons. Elevated glutamate and glutamine (Glx) and lower N-acetylaspartate (NAA) in the frontal white matter, basal ganglia, and mid frontal cortex — were associated with lower white matter, putamen and thalamus volumes, and ventricular and CSF space expansion. Reductions in brain volumes in the setting of chronic and stable disease are strongly linked to a history of immunosuppression, suggesting that delays in initiating cART may result in imminent and irreversible brain damage.

•

We mapped the 3D pattern of brain abnormalities in 210 HIV patients on stable cART.

•

Brain atrophy was linked to MRS metabolite disturbances reflecting neuronal injury.

•

Lower nadir CD4 + count was associated with greater white matter atrophy.

A multivariate surface-based analysis of the putamen in premature newborns: regional differences within the ventral striatum

Many children born preterm exhibit frontal executive dysfunction, behavioral problems including attentional deficit/hyperactivity disorder and attention related learning disabilities. Anomalies in regional specificity of cortico-striato-thalamo-cortical circuits may underlie deficits in these disorders. Nonspecific volumetric deficits of striatal structures have been documented in these subjects, but little is known about surface deformation in these structures. For the first time, here we found regional surface morphological differences in the preterm neonatal ventral striatum. We performed regional group comparisons of the surface anatomy of the striatum (putamen and globus pallidus) between 17 preterm and 19 term-born neonates at term-equivalent age. We reconstructed striatal surfaces from manually segmented brain magnetic resonance images and analyzed them using our in-house conformal mapping program. All surfaces were registered to a template with a new surface fluid registration method. Vertex-based statistical comparisons between the two groups were performed via four methods: univariate and multivariate tensor-based morphometry, the commonly used medial axis distance, and a combination of the last two statistics. We found statistically significant differences in regional morphology between the two groups that are consistent across statistics, but more extensive for multivariate measures. Differences were localized to the ventral aspect of the striatum. In particular, we found abnormalities in the preterm anterior/inferior putamen, which is interconnected with the medial orbital/prefrontal cortex and the midline thalamic nuclei including the medial dorsal nucleus and pulvinar. These findings support the hypothesis that the ventral striatum is vulnerable, within the cortico-stiato-thalamo-cortical neural circuitry, which may underlie the risk for long-term development of frontal executive dysfunction, attention deficit hyperactivity disorder and attention-related learning disabilities in preterm neonates.

Spatial-temporal atlas of human fetal brain development during the early second trimester

During the second trimester, the human fetal brain undergoes numerous changes that lead to substantial variation in the neonatal in terms of its morphology and tissue types. As fetal MRI is more and more widely used for studying the human brain development during this period, a spatiotemporal atlas becomes necessary for characterizing the dynamic structural changes. In this study, 34 postmortem human fetal brains with gestational ages ranging from 15 to 22 weeks were scanned using 7.0 T MR. We used automated morphometrics, tensor-based morphometry and surface modeling techniques to analyze the data. Spatiotemporal atlases of each week and the overall atlas covering the whole period with high resolution and contrast were created. These atlases were used for the analysis of age-specific shape changes during this period, including development of the cerebral wall, lateral ventricles, Sylvian fissure, and growth direction based on local surface measurements. Our findings indicate that growth of the subplate zone is especially striking and is the main cause for the lamination pattern changes. Changes in the cortex around Sylvian fissure demonstrate that cortical growth may be one of the mechanisms for gyration. Surface deformation mapping, revealed by local shape analysis, indicates that there is global anterior-posterior growth pattern, with frontal and temporal lobes developing relatively quickly during this period. Our results are valuable for understanding the normal brain development trajectories and anatomical characteristics. These week-by-week fetal brain atlases can be used as reference in in vivo studies, and may facilitate the quantification of fetal brain development across space and time.