Multi-manifold diffeomorphic metric mapping for aligning cortical hemispheric surfaces

Amyloid-β Accumulation in Relation to Functional Connectivity in Aging: a Longitudinal Study

The brain undergoes many changes at pathological and functional levels in healthy aging. This study employed a longitudinal and multimodal imaging dataset from the OASIS-3 study (n=300) and explored possible relationships between amyloid beta (Aβ) accumulation and functional brain organization over time in healthy aging. We used positron emission tomography (PET) with Pittsburgh compound-B (PIB) to quantify the Aβ accumulation in the brain and resting-state functional MRI (rs-fMRI) to measure functional connectivity (FC) among brain regions. Each participant had at least 2 to 3 follow-up visits. A linear mixed-effect model was used to examine longitudinal changes of Aβ accumulation and FC throughout the whole brain. We found that the limbic and frontoparietal networks had a greater annual Aβ accumulation and a slower decline in FC in aging. Additionally, the amount of the Aβ deposition in the amygdala network at baseline slowed down the decline in its FC in aging. Furthermore, the functional connectivity of the limbic, default mode network (DMN), and frontoparietal networks accelerated the Aβ propagation across their functionally highly connected regions. The functional connectivity of the somatomotor and visual networks accelerated the Aβ propagation across the brain regions in the limbic, frontoparietal, and DMN networks. These findings suggested that the slower decline in the functional connectivity of the functional hubs may compensate for their greater Aβ accumulation in aging. The Aβ propagation from one brain region to the other may depend on their functional connectivity strength.

Predicting diagnosis 4 years prior to Alzheimer's disease incident

This study employed a deep learning longitudinal model, graph convolutional and recurrent neural network (graph-CNN-RNN), on a series of brain structural MRI scans for AD prognosis. It characterized whole-brain morphology via incorporating longitudinal cortical and subcortical morphology and defined a probabilistic risk for the prediction of AD as a function of age prior to clinical diagnosis. The graph-CNN-RNN model was trained on half of the Alzheimer's Disease Neuroimaging Initiative dataset (ADNI, n = 1559) and validated on the other half of the ADNI dataset and the Open Access Series of Imaging Studies-3 (OASIS-3, n = 930). Our findings demonstrated that the graph-CNN-RNN can reliably and robustly diagnose AD at the accuracy rate of 85% and above across all the time points for both datasets. The graph-CNN-RNN predicted the AD conversion from 0 to 4 years before the AD onset at ∼80% of accuracy. The AD probabilistic risk was associated with clinical traits, cognition, and amyloid burden assessed using [18F]-Florbetapir (AV45) positron emission tomography (PET) across all the time points. The graph-CNN-RNN provided the quantitative trajectory of brain morphology from prognosis to overt stages of AD. Such a deep learning tool and the AD probabilistic risk have great potential in clinical applications for AD prognosis.

Interethnic differences in neuroimaging markers and cognition in Asians, a population-based study

We examined interethnic differences in the prevalence of neuroimaging markers of cerebrovascular and neurodegenerative disease in 3 major Asian ethnicities (Chinese, Malays, and Indians), as well as their role in cognitive impairment. 3T MRI brain scans were acquired from 792 subjects (mean age: 70.0 ± 6.5years, 52.1% women) in the multi-ethnic Epidemiology of Dementia In Singapore study. Markers of cerebrovascular disease and neurodegeneration were identified. Cognitive performance was evaluated using Mini Mental State Examination (MMSE), Montreal Cognitive Assessment (MoCA), and a neuropsychological assessment. Compared to Chinese, Malays had a higher burden of intracranial stenosis (OR: 2.28. 95%CI: 1.23-4.20) and cortical atrophy (β: -0.60. 95%CI: -0.78, -0.41), while Indians had a higher burden of subcortical atrophy (β: -0.23. 95%CI: -0.40, -0.06). Moreover, Malay and Indian ethnicities were likely to be cognitively impaired (OR for Malays: 3.79. 95%CI: 2.29-6.26; OR for Indians: 2.87. 95%CI: 1.74-4.74) and showed worse performance in global cognition (β for Malays: -0.51. 95%CI: -0.66, -0.37; and Indians: -0.32. 95%CI: -0.47, -0.17). A higher burden of cerebrovascular and neurodegenerative markers were found in Malays and Indians when compared to Chinese. Further research is required to fully elucidate the factors and pathways that contribute to these observed differences.

Neural Transcription Correlates of Multimodal Cortical Phenotypes during Development

During development, cellular events such as cell proliferation, migration, and synaptogenesis determine the structural organization of the brain. These processes are driven in part by spatiotemporally regulated gene expression. We investigated how the genetic signatures of specific neural cell types shape cortical organization of the human brain throughout infancy and childhood. Using a transcriptional atlas and in vivo magnetic resonance imaging (MRI) data, we demonstrated time-dependent associations between the expression levels of neuronal and glial genes and cortical macro- and microstructure. Neonatal cortical phenotypes were associated with prenatal glial but not neuronal gene expression. These associations reflect cell migration and proliferation during fetal development. Childhood cortical phenotypes were associated with neuronal and astrocyte gene expression related to synaptic signaling processes, reflecting the refinement of cortical connections. These findings indicate that sequential developmental stages contribute to distinct MRI measures at different time points. This helps to bridge the gap between the genetic mechanisms driving cellular changes and widely used neuroimaging techniques.

Long-term Influences of Prenatal Maternal Depressive Symptoms on the Amygdala-Prefrontal Circuitry of the Offspring From Birth to Early Childhood

BACKGROUND

Prenatal maternal depression may have long-term impacts on amygdala-cortical development. This study explored associations of prenatal maternal depressive symptoms on the amygdala-cortical structural covariance of the offspring from birth to early childhood, derived from a longitudinal birth cohort.

METHODS

Structural magnetic resonance imaging was performed to obtain the amygdala volume and cortical thickness at each time point. Prenatal maternal depressive symptoms were measured using the Edinburgh Postnatal Depression Scale at 26 weeks of pregnancy. Regression analysis was used to examine the effects of the Edinburgh Postnatal Depression Scale on a structural coupling between the amygdala volume and cortical thickness at birth (n = 167) and 4.5 years of age (n = 199).

RESULTS

Girls whose mothers had high prenatal maternal depressive symptoms showed a positive coupling between the amygdala volume and insula thickness at birth (β = .617, p = .001) but showed a negative coupling between the amygdala volume and inferior frontal thickness at 4.5 years of age (β = -.369, p = .008). No findings were revealed in boys at any time point.

CONCLUSIONS

The development of the amygdala-prefrontal circuitry is vulnerable to environmental factors related to depression. Such a vulnerability might be sex dependent.

The Use of Physiological Signals in Brainstem/Midbrain fMRI

Brainstem and midbrain nuclei are closely linked to cognitive performance and autonomic function. To advance the localization in this area, precise functional imaging is fundamental. In this study, we used a sophisticated fMRI technique as well as physiological recordings to investigate the involvement of brainstem/midbrain nuclei in cognitive control during a Stroop task. The temporal signal-to-noise ratio (tSNR) increased due to physiological noise correction (PNC) especially in regions adjacent to arteries and cerebrospinal fluid. Within the brainstem/cerebellum template an average tSNR of 68 ± 16 was achieved after the simultaneous application of a high-resolution fMRI, specialized co-registration, and PNC. The analysis of PNC data revealed an activation of the substantia nigra in the Stroop interference contrast whereas no significant results were obtained in the midbrain or brainstem when analyzing uncorrected data. Additionally, we found that pupil size indicated the level of cognitive effort. The Stroop interference effect on pupillary responses was correlated to the effect on reaction times (R 2 = 0.464, p < 0.05). When Stroop stimuli were modulated by pupillary responses, we observed a significant activation of the LC in the Stroop interference contrast. Thus, we demonstrated the beneficial effect of PNC on data quality and statistical results when analyzing neuronal responses to a cognitive task. Parametric modulation of task events with pupillary responses improved the model of LC BOLD activations in the Stroop interference contrast.

Effects of Antenatal Maternal Depressive Symptoms and Socio-Economic Status on Neonatal Brain Development are Modulated by Genetic Risk

This study included 168 and 85 mother–infant dyads from Asian and United States of America cohorts to examine whether a genomic profile risk score for major depressive disorder (GPRS_MDD) moderates the association between antenatal maternal depressive symptoms (or socio-economic status, SES) and fetal neurodevelopment, and to identify candidate biological processes underlying such association. Both cohorts showed a significant interaction between antenatal maternal depressive symptoms and infant GPRS_MDD on the right amygdala volume. The Asian cohort also showed such interaction on the right hippocampal volume and shape, thickness of the orbitofrontal and ventromedial prefrontal cortex. Likewise, a significant interaction between SES and infant GPRS_MDD was on the right amygdala and hippocampal volumes and shapes. After controlling for each other, the interaction effect of antenatal maternal depressive symptoms and GPRS_MDD was mainly shown on the right amygdala, while the interaction effect of SES and GPRS_MDD was mainly shown on the right hippocampus. Bioinformatic analyses suggested neurotransmitter/neurotrophic signaling, SNAp REceptor complex, and glutamate receptor activity as common biological processes underlying the influence of antenatal maternal depressive symptoms on fetal cortico-limbic development. These findings suggest gene–environment interdependence in the fetal development of brain regions implicated in cognitive–emotional function. Candidate biological mechanisms involve a range of brain region-specific signaling pathways that converge on common processes of synaptic development.

Psychiatric polygenic risk associates with cortical morphology and functional organization in aging

Common brain abnormalities in cortical morphology and functional organization are observed in psychiatric disorders and aging, reflecting shared genetic influences. This preliminary study aimed to examine the contribution of a polygenetic risk for psychiatric disorders (PRS_cross) to aging brain and to identify molecular mechanisms through the use of multimodal brain images, genotypes, and transcriptome data. We showed age-related cortical thinning in bilateral inferior frontal cortex (IFC) and superior temporal gyrus and alterations in the functional connectivity between bilateral IFC and between right IFC and right inferior parietal lobe as a function of PRS_cross. Interestingly, the genes in PRS_cross, that contributed most to aging neurodegeneration, were expressed in the functioanlly connected cortical regions. Especially, genes identified through the genotype-functional connectivity association analysis were commonly expressed in both cortical regions and formed strong gene networks with biological processes related to neural plasticity and synaptogenesis, regulated by glutamatergic and GABAergic transmission, neurotrophin signaling, and metabolism. This study suggested integrating genotype and transcriptome with neuroimage data sheds new light on the mechanisms of aging brain.

Minimum action principle and shape dynamics

In this paper, we propose a new method for computing a distance between two shapes embedded in three-dimensional space. Instead of comparing directly the geometric properties of the two shapes, we measure the cost of deforming one of the two shapes into the other. The deformation is computed as the geodesic between the two shapes in the space of shapes. The geodesic is found as a minimizer of the Onsager-Machlup action, based on an elastic energy for shapes that we define. Its length is set to be the integral of the action along that path; it defines an intrinsic quasi-metric on the space of shapes. We illustrate applications of our method to geometric morphometrics using three datasets representing bones and teeth of primates. Experiments on these datasets show that the variational quasi-metric we have introduced performs remarkably well both in shape recognition and in identifying evolutionary patterns, with success rates similar to, and in some cases better than, those obtained by expert observers.

Large Deformation Multiresolution Diffeomorphic Metric Mapping for Multiresolution Cortical Surfaces: A Coarse-to-Fine Approach

Brain surface registration is an important tool for characterizing cortical anatomical variations and understanding their roles in normal cortical development and psychiatric diseases. However, surface registration remains challenging due to complicated cortical anatomy and its large differences across individuals. In this paper, we propose a fast coarse-to-fine algorithm for surface registration by adapting the large diffeomorphic deformation metric mapping (LDDMM) framework for surface mapping and show improvements in speed and accuracy via a multiresolution analysis of surface meshes and the construction of multiresolution diffeomorphic transformations. The proposed method constructs a family of multiresolution meshes that are used as natural sparse priors of the cortical morphology. At varying resolutions, these meshes act as anchor points where the parameterization of multiresolution deformation vector fields can be supported, allowing the construction of a bundle of multiresolution deformation fields, each originating from a different resolution. Using a coarse-to-fine approach, we show a potential reduction in computation cost along with improvements in sulcal alignment when compared with LDDMM surface mapping.

Modulative effects of COMT haplotype on age-related associations with brain morphology

Catechol-O-methyltransferase (COMT), located on chromosome 22q11.2, encodes an enzyme critical for dopamine flux in the prefrontal cortex. Genetic variants of COMT have been suggested to functionally manipulate prefrontal morphology and function in healthy adults. This study aims to investigate modulative roles of individuals COMT SNPs (rs737865, val158met, rs165599) and its haplotypes in age-related brain morphology using an Asian sample with 174 adults aged from 21 to 80 years. We showed an age-related decline in cortical thickness of the dorsal visual pathway, including the left dorsolateral prefrontal cortex, bilateral angular gyrus, right superior frontal cortex, and age-related shape compression in the basal ganglia as a function of the genotypes of the individual COMT SNPs, especially COMT val158met. Using haplotype trend regression analysis, COMT haplotype probabilities were estimated and further revealed an age-related decline in cortical thickness in the default mode network (DMN), including the posterior cingulate, precuneus, supramarginal and paracentral cortex, and the ventral visual system, including the occipital cortex and left inferior temporal cortex, as a function of the COMT haplotype. Our results provided new evidence on an antagonistic pleiotropic effect in COMT, suggesting that genetically programmed neural benefits in early life may have a potential bearing towards neural susceptibility in later life. Hum Brain Mapp 37:2068-2082, 2016. © 2016 Wiley Periodicals, Inc.

Optimal mass transport for shape matching and comparison

Surface based 3D shape analysis plays a fundamental role in computer vision and medical imaging. This work proposes to use optimal mass transport map for shape matching and comparison, focusing on two important applications including surface registration and shape space. The computation of the optimal mass transport map is based on Monge-Brenier theory, in comparison to the conventional method based on Monge-Kantorovich theory, this method significantly improves the efficiency by reducing computational complexity from O(n(2)) to O(n) . For surface registration problem, one commonly used approach is to use conformal map to convert the shapes into some canonical space. Although conformal mappings have small angle distortions, they may introduce large area distortions which are likely to cause numerical instability thus resulting failures of shape analysis. This work proposes to compose the conformal map with the optimal mass transport map to get the unique area-preserving map, which is intrinsic to the Riemannian metric, unique, and diffeomorphic. For shape space study, this work introduces a novel Riemannian framework, Conformal Wasserstein Shape Space, by combing conformal geometry and optimal mass transport theory. In our work, all metric surfaces with the disk topology are mapped to the unit planar disk by a conformal mapping, which pushes the area element on the surface to a probability measure on the disk. The optimal mass transport provides a map from the shape space of all topological disks with metrics to the Wasserstein space of the disk and the pullback Wasserstein metric equips the shape space with a Riemannian metric. We validate our work by numerous experiments and comparisons with prior approaches and the experimental results demonstrate the efficiency and efficacy of our proposed approach.

Multiresolution Diffeomorphic Mapping for Cortical Surfaces

Due to the convoluted folding pattern of the cerebral cortex, accurate alignment of cortical surfaces remains challenging. In this paper, we present a multiresolution diffeomorphic surface mapping algorithm under the framework of large deformation diffeomorphic metric mapping (LDDMM). Our algorithm takes advantage of multiresolution analysis (MRA) for surfaces and constructs cortical surfaces at multiresolution. This family of multiresolution surfaces are used as natural sparse priors of the cortical anatomy and provide the anchor points where the parametrization of deformation vector fields is supported. This naturally constructs tangent bundles of diffeomorphisms at different resolution levels and hence generates multiresolution diffeomorphic transformation. We show that our construction of multiresolution LDDMM surface mapping can potentially reduce computational cost and improves the mapping accuracy of cortical surfaces.

Unified heat kernel regression for diffusion, kernel smoothing and wavelets on manifolds and its application to mandible growth modeling in CT images

We present a novel kernel regression framework for smoothing scalar surface data using the Laplace-Beltrami eigenfunctions. Starting with the heat kernel constructed from the eigenfunctions, we formulate a new bivariate kernel regression framework as a weighted eigenfunction expansion with the heat kernel as the weights. The new kernel method is mathematically equivalent to isotropic heat diffusion, kernel smoothing and recently popular diffusion wavelets. The numerical implementation is validated on a unit sphere using spherical harmonics. As an illustration, the method is applied to characterize the localized growth pattern of mandible surfaces obtained in CT images between ages 0 and 20 by regressing the length of displacement vectors with respect to a surface template.

Developmental synchrony of thalamocortical circuits in the neonatal brain

The thalamus is a deep gray matter structure and consists of axonal fibers projecting to the entire cortex, which provide the anatomical support for its sensorimotor and higher-level cognitive functions. There is limited in vivo evidence on the normal thalamocortical development, especially in early life. In this study, we aimed to investigate the developmental patterns of the cerebral cortex, the thalamic substructures, and their connectivity with the cortex in the first few weeks of the postnatal brain. We hypothesized that there is developmental synchrony of the thalamus, its cortical projections, and corresponding target cortical structures. We employed diffusion tensor imaging (DTI) and divided the thalamus into five substructures respectively connecting to the frontal, precentral, postcentral, temporal, and parietal and occipital cortex. T2-weighted magnetic resonance imaging (MRI) was used to measure cortical thickness. We found age-related increases in cortical thickness of bilateral frontal cortex and left temporal cortex in the early postnatal brain. We also found that the development of the thalamic substructures was synchronized with that of their respective thalamocortical connectivity in the first few weeks of the postnatal life. In particular, the right thalamo-frontal substructure had the fastest growth in the early postnatal brain. Our study suggests that the distinct growth patterns of the thalamic substructures are in synchrony with those of the cortex in early life, which may be critical for the development of the cortical and subcortical functional specialization.

Individualized diffeomorphic mapping of brains with large cortical infarcts

Whole brain mapping of stroke patients with large cortical infarcts is not trivial due to the complexity of infarcts' anatomical location and appearance in magnetic resonance image. In this study, we proposed an individualized diffeomorphic mapping framework for solving this problem. This framework is based on our recent work of large deformation diffeomorphic metric mapping (LDDMM) in Du et al. (2011) and incorporates anatomical features, such as sulcal/gyral curves, cortical surfaces, brain intensity image, and masks of infarcted regions, in order to align a normal brain to the brain of stroke patients. We applied this framework to synthetic data and data of stroke patients and validated the mapping accuracy in terms of the alignment of gyral/sulcal curves, sulcal regions, and brain segmentation. Our results revealed that this framework provided comparable mapping results for stroke patients and healthy controls, suggesting the importance of incorporating individualized anatomical features in whole brain mapping of brains with large cortical infarcts.

Gauss-Newton inspired preconditioned optimization in large deformation diffeomorphic metric mapping

In this work, we propose a novel preconditioned optimization method in the paradigm of Large Deformation Diffeomorphic Metric Mapping (LDDMM). The preconditioned update scheme is formulated for the non-stationary and the stationary parameterizations of diffeomorphisms, yielding three different LDDMM methods. The preconditioning matrices are inspired in the Hessian approximation used in Gauss-Newton method. The derivatives are computed using Frechet differentials. Thus, optimization is performed in a Sobolev space, in contrast to optimization in L(2) commonly used in non-rigid registration literature. The proposed LDDMM methods have been evaluated and compared with their respective implementations of gradient descent optimization. Evaluation has been performed using real and simulated images from the Non-rigid Image Registration Evaluation Project (NIREP). The experiments conducted in this work reported that our preconditioned LDDMM methods achieved a performance similar or superior to well-established-in-literature gradient descent non-stationary LDDMM in the great majority of cases. Moreover, preconditioned optimization showed a substantial reduction in the execution time with an affordable increase of the memory usage per iteration. Additional experiments reported that optimization using Frechet differentials should be preferable to optimization using L(2) differentials.

Three-Class Differential Diagnosis among Alzheimer Disease, Frontotemporal Dementia, and Controls

Biomarkers derived from brain magnetic resonance (MR) imaging have promise in being able to assist in the clinical diagnosis of brain pathologies. These have been used in many studies in which the goal has been to distinguish between pathologies such as Alzheimer's disease and healthy aging. However, other dementias, in particular, frontotemporal dementia, also present overlapping pathological brain morphometry patterns. Hence, a classifier that can discriminate morphometric features from a brain MRI from the three classes of normal aging, Alzheimer's disease (AD), and frontotemporal dementia (FTD) would offer considerable utility in aiding in correct group identification. Compared to the conventional use of multiple pair-wise binary classifiers that learn to discriminate between two classes at each stage, we propose a single three-way classification system that can discriminate between three classes at the same time. We present a novel classifier that is able to perform a three-class discrimination test for discriminating among AD, FTD, and normal controls (NC) using volumes, shape invariants, and local displacements (three features) of hippocampi and lateral ventricles (two structures times two hemispheres individually) obtained from brain MR images. In order to quantify its utility in correct discrimination, we optimize the three-class classifier on a training set and evaluate its performance using a separate test set. This is a novel, first-of-its-kind comparative study of multiple individual biomarkers in a three-class setting. Our results demonstrate that local atrophy features in lateral ventricles offer the potential to be a biomarker in discriminating among AD, FTD, and NC in a three-class setting for individual patient classification.

Effects of the neurogranin variant rs12807809 on thalamocortical morphology in schizophrenia

Although the genome wide supported psychosis susceptibility neurogranin (NRGN) gene is expressed in human brains, it is unclear how it impacts brain morphology in schizophrenia. We investigated the influence of NRGN rs12807809 on cortical thickness, subcortical volumes and shapes in patients with schizophrenia. One hundred and fifty six subjects (91 patients with schizophrenia and 65 healthy controls) underwent structural MRI scans and their blood samples were genotyped. A brain mapping algorithm, large deformation diffeomorphic metric mapping, was used to perform group analysis of subcortical shapes and cortical thickness. Patients with risk TT genotype were associated with widespread cortical thinning involving frontal, parietal and temporal cortices compared with controls with TT genotype. No volumetric difference in subcortical structures (hippocampus, thalamus, amygdala, basal ganglia) was observed between risk TT genotype in patients and controls. However, patients with risk TT genotype were associated with thalamic shape abnormalities involving regions related to pulvinar and medial dorsal nuclei. Our results revealed the influence of the NRGN gene on thalamocortical morphology in schizophrenia involving widespread cortical thinning and thalamic shape abnormalities. These findings help to clarify underlying NRGN mediated pathophysiological mechanisms involving cortical-subcortical brain networks in schizophrenia.

Association of silent lacunar infarct with brain atrophy and cognitive impairment

OBJECTIVE

Silent lacunar infarct (SLI) is associated with cognitive decline and linked to an increased risk of stroke and dementia. We examined the association of SLI with MRI measures of cortical thickness, subcortical and lateral ventricular shapes and cognition in 285 ethnic Chinese elderly.

METHODS

SLI, cortical thickness, shapes of subcortical and ventricular structures were quantified using MRI. The cognitive performance was assessed using comprehensive neuropsychological tests. Linear regression was used to examine associations among SLI, brain measures and cognition.

RESULTS

SLI was associated with atrophy in multiple subcortical structures, ventricular enlargement and widespread cortical thinning. Both SLI and atrophy were independently related to poorer performance in attention, memory and language domains. Only SLI was associated with visuomotor speed and executive function, while atrophy mediated the association between SLI and visuoconstruction.

CONCLUSIONS

Our findings support a vascular contribution to neurodegeneration and cognitive impairment.

Evolution of hippocampal shapes across the human lifespan

Aberrant hippocampal morphology plays an important role in the pathophysiology of aging. Volumetric analysis of the hippocampus has been performed in aging studies; however, the shape morphometry--which is potentially more informative in terms of related cognition--has yet to be examined. In this paper, we employed an advanced brain mapping technique, large deformation diffeomorphic metric mapping (LDDMM), and a dimensionality reduction approach, locally linear diffeomorphic metric embedding (LLDME), to explore age-related changes in hippocampal shape as delineated from magnetic resonance (MR) images of 302 healthy adults aged from 18 to 94 years. Compared with the hippocampal volumes, the hippocampal shapes clearly showed the nonlinear trajectory of biological aging across the human lifespan, where the variation of hippocampal shapes by age was characterized by a cubic polynomial. By integrating of LDDMM and LLDME, we were also able to illustrate the average hippocampal shapes in each individual decade. In addition, LDDMM and LLDME facilitated the identification of 63 years as a threshold beyond which hippocampal morphological changes were accelerated. Adults over 63 years of age showed the inward-deformation bilaterally in the head of the hippocampi and the left subiculum regardless of hippocampal volume reduction when compared to adults younger than 63. Hence, we demonstrated that the shape of anatomical structures added another dimension of structural morphological quantification beyond the volume in understanding aging.

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.

Ricci Flow-based Spherical Parameterization and Surface Registration

This paper presents an improved Euclidean Ricci flow method for spherical parameterization. We subsequently invent a scale space processing built upon Ricci energy to extract robust surface features for accurate surface registration. Since our method is based on the proposed Euclidean Ricci flow, it inherits the properties of Ricci flow such as conformality, robustness and intrinsicalness, facilitating efficient and effective surface mapping. Compared with other surface registration methods using curvature or sulci pattern, our method demonstrates a significant improvement for surface registration. In addition, Ricci energy can capture local differences for surface analysis as shown in the experiments and applications.

Abnormalities of cortical thickness, subcortical shapes, and white matter integrity in subcortical vascular cognitive impairment

Subcortical vascular cognitive impairment (sVCI) is caused by lacunar infarcts or extensive and/or diffuse lesions in the white matter that may disrupt the white matter circuitry connecting cortical and subcortical regions and result in the degeneration of neurons in these regions. This study used structural magnetic resonance imaging (MRI) and high angular resolution diffusion imaging (HARDI) techniques to examine cortical thickness, subcortical shapes, and white matter integrity in mild vascular cognitive impairment no dementia (VCIND Mild) and moderate-to-severe VCI (MSVCI). Our study found that compared to controls (n = 25), VCIND Mild (n = 25), and MSVCI (n = 30) showed thinner cortex predominantly in the frontal cortex. The cortex in MSVCI was thinner in the parietal and lateral temporal cortices than that in VCIND Mild. Moreover, compared to controls, VCIND Mild and MSVCI showed smaller shapes (i.e., volume reduction) in the thalamus, putamen, and globus pallidus and ventricular enlargement. Finally, compared to controls, VCIND Mild, and MSVCI showed an increased mean diffusivity in the white matter, while decreased generalized fractional anisotropy was only found in the MSVCI subjects. The major axonal bundles involved in the white matter abnormalities were mainly toward the frontal regions, including the internal capsule/corona radiata, uncinate fasciculus, and anterior section of the inferior fronto-occipital fasciculus, and were anatomically connected to the affected cortical and subcortical structures. Our findings suggest that abnormalities in cortical, subcortical, and white matter morphology in sVCI occur in anatomically connected structures, and that abnormalities progress along a similar trajectory from the mild to moderate and severe conditions.

Amygdala-hippocampal shape and cortical thickness abnormalities in first-episode schizophrenia and mania

BACKGROUND

Abnormalities in cortical thickness and subcortical structures have been studied in schizophrenia but little is known about corresponding changes in mania and brain structural differences between these two psychiatric conditions, especially early in the stage of the illness. In this study we aimed to compare cortical thickness and shape of the amygdala-hippocampal complex in first-episode schizophrenia (FES) and mania (FEM). Method Structural magnetic resonance imaging (MRI) was performed on 28 FES patients, 28 FEM patients and 28 healthy control subjects who were matched for age, gender and handedness.

RESULTS

Overall, the shape of the amygdala was deformed in both patient groups, relative to controls. Compared to FEM patients, FES patients had significant inward shape deformation in the left hippocampal tail, right hippocampal body and a small region in the right amygdala. Cortical thinning was more widespread in FES patients, with significant differences found in the temporal brain regions when compared with FEM and controls.

CONCLUSIONS

Significant differences were observed between the two groups of patients with FES and FEM in terms of the hippocampal shape and cortical thickness in the temporal region, highlighting that distinguishable brain structural changes are present early in the course of schizophrenia and mania.

Setup and data analysis for functional magnetic resonance imaging of awake cat visual cortex

Functional magnetic resonance imaging (fMRI) is one of the most commonly used methods in cognitive neuroscience on humans. In recent decades, fMRI has also been used in the awake monkey experiments to localize functional brain areas and to compare the functional differences between human and monkey brains. Several procedures and paradigms have been developed to maintain proper head fixation and to perform motion control training. In this study, we extended the application of fMRI to awake cats without training, receiving a flickering checkerboard visual stimulus projected to a screen in front of them in a block-design paradigm. We found that body movement-induced non-rigid motion introduced artifacts into the functional scans, especially those around the eye and neck. To correct for these artifacts, we developed two methods: one for general experimental design, and the other for studies of whether a checkerboard task could be used as a localizer to optimize the motion-correction parameters. The results demonstrated that, with proper animal fixation and motion correction procedures, it is possible to perform fMRI experiments with untrained awake cats.

Functional networks in parallel with cortical development associate with executive functions in children

Children begin performing similarly to adults on tasks requiring executive functions in late childhood, a transition that is probably due to neuroanatomical fine-tuning processes, including myelination and synaptic pruning. In parallel to such structural changes in neuroanatomical organization, development of functional organization may also be associated with cognitive behaviors in children. We examined 6- to 10-year-old children's cortical thickness, functional organization, and cognitive performance. We used structural magnetic resonance imaging (MRI) to identify areas with cortical thinning, resting-state fMRI to identify functional organization in parallel to cortical development, and working memory/response inhibition tasks to assess executive functioning. We found that neuroanatomical changes in the form of cortical thinning spread over bilateral frontal, parietal, and occipital regions. These regions were engaged in 3 functional networks: sensorimotor and auditory, executive control, and default mode network. Furthermore, we found that working memory and response inhibition only associated with regional functional connectivity, but not topological organization (i.e., local and global efficiency of information transfer) of these functional networks. Interestingly, functional connections associated with "bottom-up" as opposed to "top-down" processing were more clearly related to children's performance on working memory and response inhibition, implying an important role for brain systems involved in late childhood.

Applying tensor-based morphometry to parametric surfaces can improve MRI-based disease diagnosis

Many methods have been proposed for computer-assisted diagnostic classification. Full tensor information and machine learning with 3D maps derived from brain images may help detect subtle differences or classify subjects into different groups. Here we develop a new approach to apply tensor-based morphometry to parametric surface models for diagnostic classification. We use this approach to identify cortical surface features for use in diagnostic classifiers. First, with holomorphic 1-forms, we compute an efficient and accurate conformal mapping from a multiply connected mesh to the so-called slit domain. Next, the surface parameterization approach provides a natural way to register anatomical surfaces across subjects using a constrained harmonic map. To analyze anatomical differences, we then analyze the full Riemannian surface metric tensors, which retain multivariate information on local surface geometry. As the number of voxels in a 3D image is large, sparse learning is a promising method to select a subset of imaging features and to improve classification accuracy. Focusing on vertices with greatest effect sizes, we train a diagnostic classifier using the surface features selected by an L1-norm based sparse learning method. Stability selection is applied to validate the selected feature sets. We tested the algorithm on MRI-derived cortical surfaces from 42 subjects with genetically confirmed Williams syndrome and 40 age-matched controls, multivariate statistics on the local tensors gave greater effect sizes for detecting group differences relative to other TBM-based statistics including analysis of the Jacobian determinant and the largest eigenvalue of the surface metric. Our method also gave reasonable classification results relative to the Jacobian determinant, the pair of eigenvalues of the Jacobian matrix and volume features. This analysis pipeline may boost the power of morphometry studies, and may assist with image-based classification.

Inattention and hyperactivity predict alterations in specific neural circuits among 6-year-old boys

OBJECTIVE

Assessment of inattention and hyperactivity in preschoolers is highly dependent upon parental reports. Such reports are compromised by parental attitudes and mental health. Our study aimed to examine associations of inattention and hyperactivity/impulsivity from maternal reports on the Conners' Parent Rating Scale (CPRS) with brain morphology assessed using structural magnetic resonance imaging (MRI) and diffusion tensor imaging (DTI) in 6-year-old boys.

METHOD

Large deformation diffeomorphic metric brain mapping was used to assess brain morphology on MRI and DTI in 96 six-year-old boys, including cortical thickness, subcortical shapes, and fractional anisotropy (FA) of deep white matter tracts (DWMTs). Linear regression examined associations between these measures of brain structures and mothers' CPRS ratings of their child's inattention and hyperactivity/impulsivity.

RESULTS

Our results revealed that temporal and parietal cortices, as well as posterior white matter and callosal tracts are associated with inattention and hyperactivity/impulsivity symptoms among six-year-old boys. Inattention and hyperactivity/impulsivity symptoms share common neural circuits, but hyperactivity/impulsivity ratings associate with more extensive cortical areas, such as frontal regions, and with white matter tracts emphasizing executive control. There were no associations detected between inattention (or hyperactivity/impulsivity) and the shape of subcortical structures.

CONCLUSIONS

Our results suggested specific rather than widespread neural circuits involved in inattention and hyperactivity/impulsivity in young children, which is congruent with existing findings in older children and adolescents, and in adults with attention-deficit/hyperactivity disorder (ADHD). Hence, our study supported the dimensional view of ADHD, that is, that symptoms of inattention and hyperactivity/impulsivity lie on a continuum.

Emotional expressions in voice and music: same code, same effect?

Scholars have documented similarities in the way voice and music convey emotions. By using functional magnetic resonance imaging (fMRI) we explored whether these similarities imply overlapping processing substrates. We asked participants to trace changes in either the emotion or pitch of vocalizations and music using a joystick. Compared to music, vocalizations more strongly activated superior and middle temporal cortex, cuneus, and precuneus. However, despite these differences, overlapping rather than differing regions emerged when comparing emotion with pitch tracing for music and vocalizations, respectively. Relative to pitch tracing, emotion tracing activated medial superior frontal and anterior cingulate cortex regardless of stimulus type. Additionally, we observed emotion specific effects in primary and secondary auditory cortex as well as in medial frontal cortex that were comparable for voice and music. Together these results indicate that similar mechanisms support emotional inferences from vocalizations and music and that these mechanisms tap on a general system involved in social cognition.

Brain surface conformal parameterization with the Ricci flow

In brain mapping research, parameterized 3-D surface models are of great interest for statistical comparisons of anatomy, surface-based registration, and signal processing. Here, we introduce the theories of continuous and discrete surface Ricci flow, which can create Riemannian metrics on surfaces with arbitrary topologies with user-defined Gaussian curvatures. The resulting conformal parameterizations have no singularities and they are intrinsic and stable. First, we convert a cortical surface model into a multiple boundary surface by cutting along selected anatomical landmark curves. Secondly, we conformally parameterize each cortical surface to a parameter domain with a user-designed Gaussian curvature arrangement. In the parameter domain, a shape index based on conformal invariants is computed, and inter-subject cortical surface matching is performed by solving a constrained harmonic map. We illustrate various target curvature arrangements and demonstrate the stability of the method using longitudinal data. To map statistical differences in cortical morphometry, we studied brain asymmetry in 14 healthy control subjects. We used a manifold version of Hotelling's T(2) test, applied to the Jacobian matrices of the surface parameterizations. A permutation test, along with the cumulative distribution of p-values, were used to estimate the overall statistical significance of differences. The results show our algorithm's power to detect subtle group differences in cortical surfaces.

Principal component based diffeomorphic surface mapping

We present a new diffeomorphic surface mapping algorithm under the framework of large deformation diffeomorphic metric mapping (LDDMM). Unlike existing LDDMM approaches, this new algorithm reduces the complexity of the estimation of diffeomorphic transformations by incorporating a shape prior in which a nonlinear diffeomorphic shape space is represented by a linear space of initial momenta of diffeomorphic geodesic flows from a fixed template. In addition, for the first time, the diffeomorphic mapping is formulated within a decision-theoretic scheme based on Bayesian modeling in which an empirical shape prior is characterized by a low dimensional Gaussian distribution on initial momentum. This is achieved using principal component analysis (PCA) to construct the eigenspace of the initial momentum. A likelihood function is formulated as the conditional probability of observing surfaces given any particular value of the initial momentum, which is modeled as a random field of vector-valued measures characterizing the geometry of surfaces. We define the diffeomorphic mapping as a problem that maximizes a posterior distribution of the initial momentum given observable surfaces over the eigenspace of the initial momentum. We demonstrate the stability of the initial momentum eigenspace when altering training samples using a bootstrapping method. We then validate the mapping accuracy and show robustness to outliers whose shape variation is not incorporated into the shape prior.

Approximations of the diffeomorphic metric and their applications in shape learning

In neuroimaging studies based on anatomical shapes, it is well-known that the dimensionality of the shape information is much higher than the number of subjects available. A major challenge in shape analysis is to develop a dimensionality reduction approach that is able to efficiently characterize anatomical variations in a low-dimensional space. For this, there is a need to characterize shape variations among individuals for N given subjects. Therefore, one would need to calculate (2(N)) mappings between any two shapes and obtain their distance matrix. In this paper, we propose a method that reduces the computational burden to N mappings. This is made possible by making use of the first- and second-order approximations of the metric distance between two brain structural shapes in a diffeomorphic metric space. We directly derive these approximations based on the so-called conservation law of momentum, i.e., the diffeomorphic transformation acting on anatomical shapes along the geodesic is completely determined by its velocity at the origin of a fixed template. This allows for estimating morphological variation of two shapes through the first- and second-order approximations of the initial velocity in the tangent space of the diffeomorphisms at the template. We also introduce an alternative representation of these approximations through the initial momentum, i.e., a linear transformation of the initial velocity, and provide a simple computational algorithm for the matrix of the diffeomorphic metric. We employ this algorithm to compute the distance matrix of hippocampal shapes among an aging population used in a dimensionality reduction analysis, namely, ISOMAP. Our results demonstrate that the first- and second-order approximations are sufficient to characterize shape variations when compared to the diffeomorphic metric constructed through (2(N)) mappings in ISOMAP analysis.

Diffeomorphic brain registration under exhaustive sulcal constraints

The alignment and normalization of individual brain structures is a prerequisite for group-level analyses of structural and functional neuroimaging data. The techniques currently available are either based on volume and/or surface attributes, with limited insight regarding the consistent alignment of anatomical landmarks across individuals. This article details a global, geometric approach that performs the alignment of the exhaustive sulcal imprints (cortical folding patterns) across individuals. This DIffeomorphic Sulcal-based COrtical (DISCO) technique proceeds to the automatic extraction, identification and simplification of sulcal features from T1-weighted Magnetic Resonance Image (MRI) series. These features are then used as control measures for fully-3-D diffeomorphic deformations. Quantitative and qualitative evaluations show that DISCO correctly aligns the sulcal folds and gray and white matter volumes across individuals. The comparison with a recent, iconic diffeomorphic approach (DARTEL) highlights how the absence of explicit cortical landmarks may lead to the misalignment of cortical sulci. We also feature DISCO in the automatic design of an empirical sulcal template from group data. We also demonstrate how DISCO can efficiently be combined with an image-based deformation (DARTEL) to further improve the consistency and accuracy of alignment performances. Finally, we illustrate how the optimized alignment of cortical folds across subjects improves sensitivity in the detection of functional activations in a group-level analysis of neuroimaging data.

Consistent sulcal parcellation of longitudinal cortical surfaces

Automated accurate and consistent sulcal parcellation of longitudinal cortical surfaces is of great importance in studying longitudinal morphological and functional changes of human brains, since longitudinal cortical changes are normally very subtle, especially in aging brains. However, applying the existing methods (which were typically developed for cortical sulcal parcellation of a single cortical surface) independently to longitudinal cortical surfaces might generate longitudinally-inconsistent results. To overcome this limitation, this paper presents a novel energy function based method for accurate and consistent sulcal parcellation of longitudinal cortical surfaces. Specifically, both spatial and temporal smoothness are imposed in the energy function to obtain consistent longitudinal sulcal parcellation results. The energy function is efficiently minimized by a graph cut method. The proposed method has been successfully applied to sulcal parcellation of both real and simulated longitudinal inner cortical surfaces of human brain MR images. Both qualitative and quantitative evaluation results demonstrate the validity of the proposed method.

Whole brain diffeomorphic metric mapping via integration of sulcal and gyral curves, cortical surfaces, and images

This paper introduces a novel large deformation diffeomorphic metric mapping algorithm for whole brain registration where sulcal and gyral curves, cortical surfaces, and intensity images are simultaneously carried from one subject to another through a flow of diffeomorphisms. To the best of our knowledge, this is the first time that the diffeomorphic metric from one brain to another is derived in a shape space of intensity images and point sets (such as curves and surfaces) in a unified manner. We describe the Euler-Lagrange equation associated with this algorithm with respect to momentum, a linear transformation of the velocity vector field of the diffeomorphic flow. The numerical implementation for solving this variational problem, which involves large-scale kernel convolution in an irregular grid, is made feasible by introducing a class of computationally friendly kernels. We apply this algorithm to align magnetic resonance brain data. Our whole brain mapping results show that our algorithm outperforms the image-based LDDMM algorithm in terms of the mapping accuracy of gyral/sulcal curves, sulcal regions, and cortical and subcortical segmentation. Moreover, our algorithm provides better whole brain alignment than combined volumetric and surface registration (Postelnicu et al., 2009) and hierarchical attribute matching mechanism for elastic registration (HAMMER) (Shen and Davatzikos, 2002) in terms of cortical and subcortical volume segmentation.

Locally Linear Diffeomorphic Metric Embedding (LLDME) for surface-based anatomical shape modeling

This paper presents the algorithm, Locally Linear Diffeomorphic Metric Embedding (LLDME), for constructing efficient and compact representations of surface-based brain shapes whose variations are characterized using Large Deformation Diffeomorphic Metric Mapping (LDDMM). Our hypothesis is that the shape variations in the infinite-dimensional diffeomorphic metric space can be captured by a low-dimensional space. To do so, traditional Locally Linear Embedding (LLE) that reconstructs a data point from its neighbors in Euclidean space is extended to LLDME that requires interpolating a shape from its neighbors in the infinite-dimensional diffeomorphic metric space. This is made possible through the conservation law of momentum derived from LDDMM. It indicates that initial momentum, a linear transformation of the initial velocity of diffeomorphic flows, at a fixed template shape determines the geodesic connecting the template to a subject's shape in the diffeomorphic metric space and becomes the shape signature of an individual subject. This leads to the compact linear representation of the nonlinear diffeomorphisms in terms of the initial momentum. Since the initial momentum is in a linear space, a shape can be approximated by a linear combination of its neighbors in the diffeomorphic metric space. In addition, we provide efficient computations for the metric distance between two shapes through the first order approximation of the geodesic using the initial momentum as well as for the reconstruction of a shape given its low-dimensional Euclidean coordinates using the geodesic shooting with the initial momentum as the initial condition. Experiments are performed on the hippocampal shapes of 302 normal subjects across the whole life span (18-94years). Compared with Principal Component Analysis and ISOMAP, LLDME provides the most compact and efficient representation of the age-related hippocampal shapes. Even though the hippocampal volumes among young adults are as variable as those in older adults, LLDME disentangles the hippocampal local shape variation from the hippocampal size and thus reveals the nonlinear relationship of the hippocampal morphometry with age.

Quantitative evaluation of LDDMM, FreeSurfer, and CARET for cortical surface mapping

Cortical surface mapping has been widely used to compensate for individual variability of cortical shape and topology in anatomical and functional studies. While many surface mapping methods were proposed based on landmarks, curves, spherical or native cortical coordinates, few studies have extensively and quantitatively evaluated surface mapping methods across different methodologies. In this study we compared five cortical surface mapping algorithms, including large deformation diffeomorphic metric mapping (LDDMM) for curves (LDDMM-curve), for surfaces (LDDMM-surface), multi-manifold LDDMM (MM-LDDMM), FreeSurfer, and CARET, using 40 MRI scans and 10 simulated datasets. We computed curve variation errors and surface alignment consistency for assessing the mapping accuracy of local cortical features (e.g., gyral/sulcal curves and sulcal regions) and the curvature correlation for measuring the mapping accuracy in terms of overall cortical shape. In addition, the simulated datasets facilitated the investigation of mapping error distribution over the cortical surface when the MM-LDDMM, FreeSurfer, and CARET mapping algorithms were applied. Our results revealed that the LDDMM-curve, MM-LDDMM, and CARET approaches best aligned the local curve features with their own curves. The MM-LDDMM approach was also found to be the best in aligning the local regions and cortical folding patterns (e.g., curvature) as compared to the other mapping approaches. The simulation experiment showed that the MM-LDDMM mapping yielded less local and global deformation errors than the CARET and FreeSurfer mappings.