Validation of Alternating Kernel Mixture Method: Application to Tissue Segmentation of Cortical and Subcortical Structures

This paper describes the application of the alternating Kernel mixture (AKM) segmentation algorithm to high resolution MRI subvolumes acquired from a 1.5T scanner (hippocampus, n = 10 and prefrontal cortex, n = 9) and a 3T scanner (hippocampus, n = 10 and occipital lobe, n = 10). Segmentation of the subvolumes into cerebrospinal fluid, gray matter, and white matter tissue is validated by comparison with manual segmentation. When compared with other segmentation methods that use traditional Bayesian segmentation, AKM yields smaller errors (P < .005, exact Wilcoxon signed rank test) demonstrating the robustness and wide applicability of AKM across different structures. By generating multiple mixtures for each tissue compartment, AKM mimics the increased variation of manual segmentation in partial volumes due to the highly folded tissues. AKM's superior performance makes it useful for tissue segmentation of subcortical and cortical structures in large-scale neuroimaging studies

3D Normal Coordinate Systems for Cortical Areas

A surface-based diffeomorphic algorithm to generate 3D coordinate grids in the cortical ribbon is described. In the grid, normal coordinate lines are generated by the diffeomorphic evolution from the grey/white (inner) surface to the grey/csf (outer) surface. Specifically, the cortical ribbon is described by two triangulated surfaces with open boundaries. Conceptually, the inner surface sits on top of the white matter structure and the outer on top of the gray matter. It is assumed that the cortical ribbon consists of cortical columns which are orthogonal to the white matter surface. This might be viewed as a consequence of the development of the columns in the embryo. It is also assumed that the columns are orthogonal to the outer surface so that the resultant vector field is orthogonal to the evolving surface. Then the distance of the normal lines from the vector field such that the inner surface evolves diffeomorphically towards the outer one can be construed as a measure of thickness. Applications are described for the auditory cortices in human adults and cats with normal hearing or hearing loss. The approach offers great potential for cortical morphometry

Large Deformation Diffeomorphic Metric Mapping Registration of Reconstructed 3D Histological Section Images and in vivo MR Images

Our current understanding of neuroanatomical abnormalities in neuropsychiatric diseases is based largely on magnetic resonance imaging (MRI) and post mortem histological analyses of the brain. Further advances in elucidating altered brain structure in these human conditions might emerge from combining MRI and histological methods. We propose a multistage method for registering 3D volumes reconstructed from histological sections to corresponding in vivo MRI volumes from the same subjects: (1) manual segmentation of white matter (WM), gray matter (GM) and cerebrospinal fluid (CSF) compartments in histological sections, (2) alignment of consecutive histological sections using 2D rigid transformation to construct a 3D histological image volume from the aligned sections, (3) registration of reconstructed 3D histological volumes to the corresponding 3D MRI volumes using 3D affine transformation, (4) intensity normalization of images via histogram matching, and (5) registration of the volumes via intensity based large deformation diffeomorphic metric (LDDMM) image matching algorithm. Here we demonstrate the utility of our method in the transfer of cytoarchitectonic information from histological sections to identify regions of interest in MRI scans of nine adult macaque brains for morphometric analyses. LDDMM improved the accuracy of the registration via decreased distances between GM/CSF surfaces after LDDMM (0.39 ± 0.13 mm) compared to distances after affine registration (0.76 ± 0.41 mm). Similarly, WM/GM distances decreased to 0.28 ± 0.16 mm after LDDMM compared to 0.54 ± 0.39 mm after affine registration. The multistage registration method may find broad application for mapping histologically based information, for example, receptor distributions, gene expression, onto MRI volumes

Patient Specific Dosimetry Phantoms Using Multichannel LDDMM of the Whole Body

This paper describes an automated procedure for creating detailed patient-specific pediatric dosimetry phantoms from a small set of segmented organs in a child's CT scan. The algorithm involves full body mappings from adult template to pediatric images using multichannel large deformation diffeomorphic metric mapping (MC-LDDMM). The parallel implementation and performance of MC-LDDMM for this application is studied here for a sample of 4 pediatric patients, and from 1 to 24 processors. 93.84% of computation time is parallelized, and the efficiency of parallelization remains high until more than 8 processors are used. The performance of the algorithm was validated on a set of 24 male and 18 female pediatric patients. It was found to be accurate typically to within 1-2 voxels (2–4 mm) and robust across this large and variable data set

Regional subcortical shape analysis in premanifest Huntington’s disease

Huntington’s disease (HD) involves preferential and progressive degeneration of striatum and other subcortical regions as well as regional cortical atrophy. It is caused by a CAG repeat expansion in the Huntingtin gene, and the longer the expansion the earlier the age of onset. Atrophy begins prior to manifest clinical signs and symptoms, and brain atrophy in premanifest expansion carriers can be studied. We employed a diffeomorphometric pipeline to contrast subcortical structures’ morphological properties in a control group with three disease groups representing different phases of premanifest HD (far, intermediate, and near to onset) as defined by the length of the CAG expansion and the participant’s age (CAG‐Age‐Product). A total of 1,428 magnetic resonance image scans from 694 participants from the PREDICT‐HD cohort were used. We found significant region‐specific atrophies in all subcortical structures studied, with the estimated abnormality onset time varying from structure to structure. Heterogeneous shape abnormalities of caudate nuclei were present in premanifest HD participants estimated furthest from onset and putaminal shape abnormalities were present in participants intermediate to onset. Thalamic, hippocampal, and amygdalar shape abnormalities were present in participants nearest to onset. We assessed whether the estimated progression of subcortical pathology in premanifest HD tracked specific pathways. This is plausible for changes in basal ganglia circuits but probably not for changes in hippocampus and amygdala. The regional shape analyses conducted in this study provide useful insights into the effects of HD pathology in subcortical structures.Peer Reviewedhttps://deepblue.lib.umich.edu/bitstream/2027.42/148249/1/hbm24456.pdfhttps://deepblue.lib.umich.edu/bitstream/2027.42/148249/2/hbm24456_am.pd

Longitudinal imaging highlights preferential basal ganglia circuit atrophy in Huntington's disease

Huntington's disease is caused by a CAG repeat expansion in the Huntingtin gene (HTT), coding for polyglutamine in the Huntingtin protein, with longer CAG repeats causing earlier age of onset. The variable 'Age' × ('CAG'-L), where 'Age' is the current age of the individual, 'CAG' is the repeat length and L is a constant (reflecting an approximation of the threshold), termed the 'CAG Age Product' (CAP) enables the consideration of many individuals with different CAG repeat expansions at the same time for analysis of any variable and graphing using the CAG Age Product score as the X axis. Structural MRI studies have showed that progressive striatal atrophy begins many years prior to the onset of diagnosable motor Huntington's disease, confirmed by longitudinal multicentre studies on three continents, including PREDICT-HD, TRACK-HD and IMAGE-HD. However, previous studies have not clarified the relationship between striatal atrophy, atrophy of other basal ganglia structures, and atrophy of other brain regions. The present study has analysed all three longitudinal datasets together using a single image segmentation algorithm and combining data from a large number of subjects across a range of CAG Age Product score. In addition, we have used a strategy of normalizing regional atrophy to atrophy of the whole brain, in order to determine which regions may undergo preferential degeneration. This made possible the detailed characterization of regional brain atrophy in relation to CAG Age Product score. There is dramatic selective atrophy of regions involved in the basal ganglia circuit-caudate, putamen, nucleus accumbens, globus pallidus and substantia nigra. Most other regions of the brain appear to have slower but steady degeneration. These results support (but certainly do not prove) the hypothesis of circuit-based spread of pathology in Huntington's disease, possibly due to spread of mutant Htt protein, though other connection-based mechanisms are possible. Therapeutic targets related to prion-like spread of pathology or other mechanisms may be suggested. In addition, they have implications for current neurosurgical therapeutic approaches, since delivery of therapeutic agents solely to the caudate and putamen may miss other structures affected early, such as nucleus accumbens and output nuclei of the striatum, the substantia nigra and the globus pallidus

Soliton Dynamics in Computational Anatomy

Computational anatomy (CA) has introduced the idea of anatomical structures being transformed by geodesic deformations on groups of diffeomorphisms. Among these geometric structures, landmarks and image outlines in CA are shown to be singular solutions of a partial differential equation that is called the geodesic EPDiff equation. A recently discovered momentum map for singular solutions of EPDiff yields their canonical Hamiltonian formulation, which in turn provides a complete parameterization of the landmarks by their canonical positions and momenta. The momentum map provides an isomorphism between landmarks (and outlines) for images and singular soliton solutions of the EPDiff equation. This isomorphism suggests a new dynamical paradigm for CA, as well as new data representation.Comment: published in NeuroImag

Normal and equivolumetric coordinate systems for cortical areas

We describe coordinate systems adapted for the space between two surfaces, such as those delineating the highly folded cortex in mammalian brains. These systems are estimated in order to satisfy geometric priors, including streamline normality or equivolumetric conditions on layers. We give a precise mathematical formulation of these problems, and present numerical simulations based on diffeomorphic registration methods, comparing them with recent approaches.Our method involves • Diffeomorphic registration of inner and outer folded folded surfaces. • Followed by equivolumetric reparametrization of layers to yield coordinate system