Cortical surface segmentation and mapping

The 3D Structural Architecture of the Human Hand Area Is Nontopographic

The functional topography of the human primary somatosensory cortex hand area is a widely studied model system to understand sensory organization and plasticity. It is so far unclear whether the underlying 3D structural architecture also shows a topographic organization. We used 7 Tesla (7T) magnetic resonance imaging (MRI) data to quantify layer-specific myelin, iron, and mineralization in relation to population receptive field maps of individual finger representations in Brodman area 3b (BA 3b) of human S1 in female and male younger adults. This 3D description allowed us to identify a characteristic profile of layer-specific myelin and iron deposition in the BA 3b hand area, but revealed an absence of structural differences, an absence of low-myelin borders, and high similarity of 3D microstructure profiles between individual fingers. However, structural differences and borders were detected between the hand and face areas. We conclude that the 3D structural architecture of the human hand area is nontopographic, unlike in some monkey species, which suggests a high degree of flexibility for functional finger organization and a new perspective on human topographic plasticity.SIGNIFICANCE STATEMENT Using ultra-high-field MRI, we provide the first comprehensive in vivo description of the 3D structural architecture of the human BA 3b hand area in relation to functional population receptive field maps. High similarity of precise finger-specific 3D profiles, together with an absence of structural differences and an absence of low-myelin borders between individual fingers, reveals the 3D structural architecture of the human hand area to be nontopographic. This suggests reduced structural limitations to cortical plasticity and reorganization and allows for shared representational features across fingers.

Structural brain abnormalities in epilepsy with myoclonic atonic seizures

OBJECTIVE

Epilepsy with myoclonic atonic seizure (EMAS) occurs in young children with previously normal to subnormal development. The outcome ranges from seizure freedom with preserved cognitive abilities to refractory epilepsy with intellectual disability (ID). Routine brain imaging typically shows no abnormalities. We aimed to compare the brain morphometry of EMAS patients with healthy subjects several years after epilepsy onset, and to correlate it to epilepsy severity and cognitive findings.

METHODS

Fourteen EMAS patients (4 females, 5-14 years) and 14 matched healthy controls were included. Patients were classified into three outcome groups (good, intermediate, poor) according to seizure control and cognitive and behavioral functioning. Individual anatomical data (T1-weighted sequence) were processed using the FreeSurfer pipeline. Cortical volume (CV), cortical thickness (CT), local gyrification index (LGI), and subcortical volumes were used for group-comparison and linear regression analyses.

RESULTS

Morphometric comparison between EMAS patients and healthy controls revealed that patients have 1) reduced CV in frontal, temporal and parietal lobes (p = <.001; 0.009 and 0.024 respectively); 2) reduced CT and LGI in frontal lobes (p = 0.036 and 0.032 respectively); and 3) a neat cerebellar volume reduction (p = 0.011). Neither the number of anti-seizure medication nor the duration of epilepsy was related to cerebellar volume (both p > 0.62). Poor outcome group was associated with lower LGI. Patients in good and intermediate outcome groups had a comparable LGI to their matched healthy controls (p > 0.27 for all lobes).

CONCLUSIONS

Structural brain differences were detectable in our sample of children with EMAS, mainly located in the frontal lobes and cerebellum. These findings are similar to those found in patients with genetic/idiopathic generalized epilepsies. Outcome groups correlated best with LGI. Whether these anatomical changes reflect genetically determined abnormal neuronal networks or a consequence of sustained epilepsy remains to be solved with prospective longitudinal studies.

Feedback contribution to surface motion perception in the human early visual cortex

Human visual surface perception has neural correlates in early visual cortex, but the role of feedback during surface segmentation in human early visual cortex remains unknown. Feedback projections preferentially enter superficial and deep anatomical layers, which provides a hypothesis for the cortical depth distribution of fMRI activity related to feedback. Using ultra-high field fMRI, we report a depth distribution of activation in line with feedback during the (illusory) perception of surface motion. Our results fit with a signal re-entering in superficial depths of V1, followed by a feedforward sweep of the re-entered information through V2 and V3. The magnitude and sign of the BOLD response strongly depended on the presence of texture in the background, and was additionally modulated by the presence of illusory motion perception compatible with feedback. In summary, the present study demonstrates the potential of depth-resolved fMRI in tackling biomechanical questions on perception.

Cortical depth profiles of luminance contrast responses in human V1 and V2 using 7 T fMRI

Neural activity in early visual cortex is modulated by luminance contrast. Cortical depth (i.e., laminar) contrast responses have been studied in monkey early visual cortex, but not in humans. In addition to the high spatial resolution needed and the ensuing low signal-to-noise ratio, laminar studies in humans using fMRI are hampered by the strong venous vascular weighting of the fMRI signal. In this study, we measured luminance contrast responses in human V1 and V2 with high-resolution fMRI at 7 T. To account for the effect of intracortical ascending veins, we applied a novel spatial deconvolution model to the fMRI depth profiles. Before spatial deconvolution, the contrast response in V1 showed a slight local maximum at mid cortical depth, whereas V2 exhibited a monotonic signal increase toward the cortical surface. After applying the deconvolution, both V1 and V2 showed a pronounced local maximum at mid cortical depth, with an additional peak in deep grey matter, especially in V1. Moreover, we found a difference in contrast sensitivity between V1 and V2, but no evidence for variations in contrast sensitivity as a function of cortical depth. These findings are in agreement with results obtained in nonhuman primates, but further research will be needed to validate the spatial deconvolution approach.

Analysis of cortical morphometric variability using labeled cortical distance maps

Morphometric (i.e., shape and size) differences in the anatomy of cortical structures are associated with neurodevelopmental and neuropsychiatric disorders. Such differences can be quantized and detected by a powerful tool called Labeled Cortical Distance Map (LCDM). The LCDM method provides distances of labeled gray matter (GM) voxels from the GM/white matter (WM) surface for specific cortical structures (or tissues). Here we describe a method to analyze morphometric variability in the particular tissue using LCDM distances. To extract more of the information provided by LCDM distances, we perform pooling and censoring of LCDM distances. In particular, we employ Brown-Forsythe (BF) test of homogeneity of variance (HOV) on the LCDM distances. HOV analysis of pooled distances provides an overall analysis of morphometric variability of the LCDMs due to the disease in question, while the HOV analysis of censored distances suggests the location(s) of significant variation in these differences (i.e., at which distance from the GM/WM surface the morphometric variability starts to be significant). We also check for the influence of assumption violations on the HOV analysis of LCDM distances. In particular, we demonstrate that BF HOV test is robust to assumption violations such as the non-normality and within sample dependence of the residuals from the median for pooled and censored distances and are robust to data aggregation which occurs in analysis of censored distances. We recommend HOV analysis as a complementary tool to the analysis of distribution/location differences. We also apply the methodology on simulated normal and exponential data sets and assess the performance of the methods when more of the underlying assumptions are satisfied. We illustrate the methodology on a real data example, namely, LCDM distances of GM voxels in ventral medial prefrontal cortices (VMPFCs) to see the effects of depression or being of high risk to depression on the morphometry of VMPFCs. The methodology used here is also valid for morphometric analysis of other cortical structures.

Volumetric relief map for intracranial cerebrospinal fluid distribution analysis

Cerebrospinal fluid imaging plays a significant role in the clinical diagnosis of brain disorders, such as hydrocephalus and Alzheimer's disease. While three-dimensional images of cerebrospinal fluid are very detailed, the complex structures they contain can be time-consuming and laborious to interpret. This paper presents a simple technique that represents the intracranial cerebrospinal fluid distribution as a two-dimensional image in such a way that the total fluid volume is preserved. We call this a volumetric relief map, and show its effectiveness in a characterization and analysis of fluid distributions and networks in hydrocephalus patients and healthy adults.

Reconstruction of the human cerebral cortex robust to white matter lesions: method and validation

Cortical atrophy has been reported in a number of diseases, such as multiple sclerosis and Alzheimer's disease, that are also associated with white matter (WM) lesions. However, most cortical reconstruction techniques do not account for these pathologies, thereby requiring additional processing to correct for the effect of WM lesions. In this work, we introduce CRUISE(+), an automated process for cortical reconstruction from magnetic resonance brain images with WM lesions. The process extends previously well validated methods to allow for multichannel input images and to accommodate for the presence of WM lesions. We provide new validation data and tools for measuring the accuracy of cortical reconstruction methods on healthy brains as well as brains with multiple sclerosis lesions. Using this data, we validate the accuracy of CRUISE(+) and compare it to another state-of-the-art cortical reconstruction tool. Our results demonstrate that CRUISE(+) has superior performance in the cortical regions near WM lesions, and similar performance in other regions.

Quantitative analysis of structural neuroimaging of mesial temporal lobe epilepsy

Mesial temporal lobe epilepsy (MTLE) is the most common of the surgically remediable drug-resistant epilepsies. MRI is the primary diagnostic tool to detect anatomical abnormalities and, when combined with EEG, can more accurately identify an epileptogenic lesion, which is often hippocampal sclerosis in cases of MTLE. As structural imaging technology has advanced the surgical treatment of MTLE and other lesional epilepsies, so too have the analysis techniques that are used to measure different structural attributes of the brain. These techniques, which are reviewed here and have been used chiefly in basic research of epilepsy and in studies of MTLE, have identified different types and the extent of anatomical abnormalities that can extend beyond the affected hippocampus. These results suggest that structural imaging and sophisticated imaging analysis could provide important information to identify networks capable of generating spontaneous seizures and ultimately help guide surgical therapy that improves postsurgical seizure-freedom outcomes.

Diffeomorphic sulcal shape analysis on the cortex

We present a diffeomorphic approach for constructing intrinsic shape atlases of sulci on the human cortex. Sulci are represented as square-root velocity functions of continuous open curves in R³, and their shapes are studied as functional representations of an infinite-dimensional sphere. This spherical manifold has some advantageous properties--it is equipped with a Riemannian L² metric on the tangent space and facilitates computational analyses and correspondences between sulcal shapes. Sulcal shape mapping is achieved by computing geodesics in the quotient space of shapes modulo scales, translations, rigid rotations, and reparameterizations. The resulting sulcal shape atlas preserves important local geometry inherently present in the sample population. The sulcal shape atlas is integrated in a cortical registration framework and exhibits better geometric matching compared to the conventional euclidean method. We demonstrate experimental results for sulcal shape mapping, cortical surface registration, and sulcal classification for two different surface extraction protocols for separate subject populations.

Collaborative Labeling of Malignant Glioma with WebMILL: A First Look

Malignant gliomas are the most common form of primary neoplasm in the central nervous system, and one of the most rapidly fatal of all human malignancies. They are treated by maximal surgical resection followed by radiation and chemotherapy. Herein, we seek to improve the methods available to quantify the extent of tumors using newly presented, collaborative labeling techniques on magnetic resonance imaging. Traditionally, labeling medical images has entailed that expert raters operate on one image at a time, which is resource intensive and not practical for very large datasets. Using many, minimally trained raters to label images has the possibility of minimizing laboratory requirements and allowing high degrees of parallelism. A successful effort also has the possibility of reducing overall cost. This potentially transformative technology presents a new set of problems, because one must pose the labeling challenge in a manner accessible to people with little or no background in labeling medical images and raters cannot be expected to read detailed instructions. Hence, a different training method has to be employed. The training must appeal to all types of learners and have the same concepts presented in multiple ways to ensure that all the subjects understand the basics of labeling. Our overall objective is to demonstrate the feasibility of studying malignant glioma morphometry through statistical analysis of the collaborative efforts of many, minimally-trained raters. This study presents preliminary results on optimization of the WebMILL framework for neoplasm labeling and investigates the initial contributions of 78 raters labeling 98 whole-brain datasets.

Cortical surface mapping using topology correction, partial flattening and 3D shape context-based non-rigid registration for use in quantifying atrophy in Alzheimer's disease

Magnetic resonance (MR) provides a non-invasive way to investigate changes in the brain resulting from aging or neurodegenerative disorders such as Alzheimer's disease (AD). Performing accurate analysis for population studies is challenging because of the interindividual anatomical variability. A large set of tools is found to perform studies of brain anatomy and population analysis (FreeSurfer, SPM, FSL). In this paper we present a newly developed surface-based processing pipeline (MILXCTE) that allows accurate vertex-wise statistical comparisons of brain modifications, such as cortical thickness (CTE). The brain is first segmented into the three main tissues: white matter, gray matter and cerebrospinal fluid, after CTE is computed, a topology corrected mesh is generated. Partial inflation and non-rigid registration of cortical surfaces to a common space using shape context are then performed. Each of the steps was firstly validated using MR images from the OASIS database. We then applied the pipeline to a sample of individuals randomly selected from the AIBL study on AD and compared with FreeSurfer. For a population of 50 individuals we found correlation of cortical thickness in all the regions of the brain (average r=0.62 left and r=0.64 right hemispheres). We finally computed changes in atrophy in 32 AD patients and 81 healthy elderly individuals. Significant differences were found in regions known to be affected in AD. We demonstrated the validity of the method for use in clinical studies which provides an alternative to well established techniques to compare different imaging biomarkers for the study of neurodegenerative diseases.

Pediatric epileptology

Challenges facing children with epilepsy are understanding the neurobiology of pharmacoresistance of epileptic encephalopathies and the development of effective surgical treatment options for those with "non-lesional" epilepsy. Although, understanding the genetics of childhood epilepsy has advanced, an effective treatment intervention has not occurred. Recently, understanding the neurobiology of hamartin and tuberin in the development of epilepsy and cognitive impairment associated with tuberous sclerosis complex allowed the development of sirolimus and everolimus to be used in human clinical trials. In spite of these breakthroughs a large number of children are likely to be outside the scope of interventional therapies. For such patients the burden of seizures is onerous and psycho-social consequences debilitating. Surgical resective options are often limited by the lack of a well defined epileptic lesion. Co-registered synthesis of advanced functional, structural and electrographic seizure onset allows identification of a focus in patients thought to have "non-lesional" epilepsy. Developments of a Pipeline for prospective data sharing are likely to increase understanding and validation of the epileptogenic zone and offer the hope of seizure freedom. Two outstanding young investigators provide a review of their exciting research and its implications in pediatric epilepsy.

Magnetic resonance relaxation and quantitative measurement in the brain

Underlying the exquisite soft tissue contrast provided by magnetic resonance imaging are the inherent biophysical processes of relaxation. Through the intricate relationships between tissue microstructure and biochemistry and the longitudinal and transverse relaxation rates, quantitative measurement of these relaxation parameters is informative of tissue change associated with disease, neural plasticity, and other biological processes. Quantitative imaging studies can further facilitate more detailed characterizations of tissue, providing a more direct link between modern MR imaging and classic histochemical and histological studies. In this chapter, we briefly review the biophysical basis of relaxation, introducing and focusing specifically on the T(1), T(2), and T(2)(*) relaxation times and detail some of the more widely used and clinically feasible techniques for their in vivo measurement. Methods for analyzing relaxation data are covered, and a summary of significant results from reported neuroimaging studies is provided. Finally, the combination of relaxation time data with other quantitative imaging data, including diffusion tensor and magnetization transfer, is examined, with the aim of providing more thorough characterization of brain tissue.

Multiple cortical surface correspondence using pairwise shape similarity

Accurately corresponding a population of human cortical surfaces provides important shape information for the diagnosis of many brain diseases. This problem is very challenging due to the highly convoluted nature of cortical surfaces. Pairwise methods using a fixed template may not handle well the case when a target cortical surface is substantially different from the template. In this paper, we develop a new method to organize the population of cortical surfaces into pairs with high shape similarity and only correspond such similar pairs to achieve a higher accuracy. In particular, we use the geometric information to identify co-located gyri and sulci for defining a new measure of shape similarity. We conduct experiments on 40 instances of the cortical surface, resulting in an improved performance over several existing shape-correspondence methods.

Volumetric topological analysis: a novel approach for trabecular bone classification on the continuum between plates and rods

Trabecular bone (TB) is a complex quasi-random network of interconnected plates and rods. TB constantly remodels to adapt to the stresses to which it is subjected (Wolff's Law). In osteoporosis, this dynamic equilibrium between bone formation and resorption is perturbed, leading to bone loss and structural deterioration. Both bone loss and structural deterioration increase fracture risk. Bone's mechanical behavior can only be partially explained by variations in bone mineral density, which led to the notion of bone structural quality. Previously, we developed digital topological analysis (DTA) which classifies plates, rods, profiles, edges, and junctions in a TB skeletal representation. Although the method has become quite popular, a major limitation of DTA is that it provides only hard classifications of different topological entities, failing to distinguish between narrow and wide plates. Here, we present a new method called volumetric topological analysis (VTA) for regional quantification of TB topology. At each TB location, the method uniquely classifies its topology on the continuum between perfect plates and perfect rods, facilitating early detections of TB alterations from plates to rods according to the known etiology of osteoporotic bone loss. Several new ideas, including manifold distance transform, manifold scale, and feature propagation have been introduced here and combined with existing DTA and distance transform methods, leading to the new VTA technology. This method has been applied to multidetector computed tomography (CT) and micro-computed tomography ( μCT) images of four cadaveric distal tibia and five distal radius specimens. Both intra- and inter-modality reproducibility of the method has been examined using repeat CT and μCT scans of distal tibia specimens. Also, the method's ability to predict experimental biomechanical properties of TB via CT imaging under in vivo conditions has been quantitatively examined and the results found are very encouraging.

Automatic parcellation of human cortical gyri and sulci using standard anatomical nomenclature

Precise localization of sulco-gyral structures of the human cerebral cortex is important for the interpretation of morpho-functional data, but requires anatomical expertise and is time consuming because of the brain's geometric complexity. Software developed to automatically identify sulco-gyral structures has improved substantially as a result of techniques providing topologically correct reconstructions permitting inflated views of the human brain. Here we describe a complete parcellation of the cortical surface using standard internationally accepted nomenclature and criteria. This parcellation is available in the FreeSurfer package. First, a computer-assisted hand parcellation classified each vertex as sulcal or gyral, and these were then subparcellated into 74 labels per hemisphere. Twelve datasets were used to develop rules and algorithms (reported here) that produced labels consistent with anatomical rules as well as automated computational parcellation. The final parcellation was used to build an atlas for automatically labeling the whole cerebral cortex. This atlas was used to label an additional 12 datasets, which were found to have good concordance with manual labels. This paper presents a precisely defined method for automatically labeling the cortical surface in standard terminology.

Atlas-based improved prediction of magnetic field inhomogeneity for distortion correction of EPI data

We describe a method for atlas-based segmentation of structural MRI for calculation of magnetic fieldmaps. CT data sets are used to construct a probabilistic atlas of the head and corresponding MR is used to train a classifier that segments soft tissue, air, and bone. Subject-specific fieldmaps are computed from the segmentations using a perturbation field model. Previous work has shown that distortion in echo-planar images can be corrected using predicted fieldmaps. We obtain results that agree well with acquired fieldmaps: 90% of voxel shifts from predicted fieldmaps show subvoxel disagreement with those computed from acquired fieldmaps. In addition, our fieldmap predictions show statistically significant improvement following inclusion of the atlas.

Patterns of cortical degeneration in an elderly cohort with cerebral small vessel disease

Emerging noninvasive neuroimaging techniques allow for the morphometric analysis of patterns of gray and white matter degeneration in vivo, which may help explain and predict the occurrence of cognitive impairment and Alzheimer's disease. A single center prospective follow-up study (Radboud University Nijmegen Diffusion tensor and Magnetic resonance imaging Cohort study (RUN DMC)) was performed involving 503 nondemented elderly individuals (50-85 years) with a history of symptomatic cerebral small vessel disease (SVD). Age was associated with a global reduction in cortical thickness, and this relationship was strongest for ventrolateral prefrontal cortex, auditory cortex, Wernicke's area, superior temporal lobe, and primary visual cortex. Right and left hemispheres differed in the thickness of language-related areas. White matter (WM) lesions were generally negatively correlated with cortical thickness, primarily in individuals over the age of 60, with the notable exception of Brodmann areas 4 and 5, which were positively correlated in age groups 50-60 and 60-70, respectively. The observed pattern of age-related decline may explain problems in memory and executive functions, which are already well documented in individuals with SVD. The additional gray matter loss affecting visual and auditory cortex, and specifically the head region of primary motor cortex, may indicate morphological correlates of impaired sensory and motor functions. The paradoxical positive relationship between WM lesion volume and cortical thickness in some areas may reflect early compensatory hypertrophy. This study raises a further interest in the mechanisms underlying cerebral gray and white matter degeneration in association with SVD, which will require further investigation with diffusion weighted and longitudinal MR studies.

Near-isometric flattening of brain surfaces

Flattened representations of brain surfaces are often used to visualize and analyze spatial patterns of structural organization and functional activity. Here, we present a set of rigorous criteria and accompanying test cases with which to evaluate flattening algorithms that attempt to preserve shortest-path distances on the original surface. We also introduce a novel flattening algorithm that is the first to satisfy all of these criteria and demonstrate its ability to produce accurate flat maps of human and macaque visual cortex. Using this algorithm, we have recently obtained results showing a remarkable, unexpected degree of consistency in the shape and topographic structure of visual cortical areas within humans and macaques, as well as between these two species.

A topology-preserving approach to the segmentation of brain images with multiple sclerosis lesions

We describe a new fully automatic method for the segmentation of brain images that contain multiple sclerosis white matter lesions. Multichannel magnetic resonance images are used to delineate multiple sclerosis lesions while segmenting the brain into its major structures. The method is an atlas-based segmentation technique employing a topological atlas as well as a statistical atlas. An advantage of this approach is that all segmented structures are topologically constrained, thereby allowing subsequent processing such as cortical unfolding or diffeomorphic shape analysis techniques. Evaluation with both simulated and real data sets demonstrates that the method has an accuracy competitive with state-of-the-art MS lesion segmentation methods, while simultaneously segmenting the whole brain.

Automatic extraction of sulcal lines on cortical surfaces based on anisotropic geodesic distance

Analyzing cortical sulci is important for studying cortical morphology and brain functions. Although sulcal lines on cortical surfaces can be defined in various ways, it is critical in a neuroimaging study to define a sulcal line along the valley of a cortical surface with a high curvature within a sulcus. To extract the sulcal lines automatically, we present a new geometric algorithm based on the computation of anisotropic skeletons of sulcal regions. Because anisotropic skeletons are highly adaptive to the anisotropic nature of the surface shape, the resulting sulcal lines lie accurately on the valleys of the sulcal areas. Our sulcal lines remain unchanged under local shape variabilities in different human brains. Through experiments, we show that the errors of the sulcal lines for both synthetic data and real cortical surfaces were nearly as constant as the function of random noise. By measuring the changes in sulcal shape in Alzheimer's disease (AD) patients, we further investigated the effectiveness of the accuracy of our sulcal lines using a large sample of MRI data. This study involved 70 normal controls (n [men/women]: 29/41, age [mean+/-SD]: 71.7+/-4.9 years), and 100 AD subjects (37/63, 72.3+/-5.5). We observe significantly lower absolute average mean curvature and shallower sulcal depth in AD subjects, where the group difference becomes more significant if we measure the quantities along the sulcal lines rather than over the entire sulcal area. The most remarkable difference in the AD patients was the average sulcal depth (control: 11.70 and AD: 11.34).

The emerging discipline of Computational Functional Anatomy

Computational Functional Anatomy (CFA) is the study of functional and physiological response variables in anatomical coordinates. For this we focus on two things: (i) the construction of bijections (via diffeomorphisms) between the coordinatized manifolds of human anatomy, and (ii) the transfer (group action and parallel transport) of functional information into anatomical atlases via these bijections. We review advances in the unification of the bijective comparison of anatomical submanifolds via point-sets including points, curves and surface triangulations as well as dense imagery. We examine the transfer via these bijections of functional response variables into anatomical coordinates via group action on scalars and matrices in DTI as well as parallel transport of metric information across multiple templates which preserves the inner product.

Depth potential function for folding pattern representation, registration and analysis

Some surfaces present folding patterns formed by juxtapositions of ridges and valleys as, for example, the cortical surface of the human brain. The fundamental problem with ridges is to find a correspondence among and analyze the variability among them. Many techniques to achieve these goals exist but use scalar functions. Depth maps are used to efficiently project the geometry of folds into a scalar function in the case where a natural projection plane exists. However, in most cases of curved surfaces, there is no natural projection plane to represent folding patterns. This paper studies the problem of shape matching and analysis of folding patterns by extending the notion of depth maps when no natural projection plane exists. The novel depth measure is called a depth potential function. The depth potential function integrates the information known from the curvature of the surface into a point-of-view invariant representation. The main advantage of the depth potential function is that it is computed by solving a time independent Poisson equation. The Poisson equation endows our surface representation with a significant computational advantage that makes it orders of magnitude faster to compute compared with other available surface representations. The method described in this paper was validated using both synthetic surfaces and cortical surfaces of human brain acquired by magnetic resonance imaging. On average, the improvement in shape matching when using the depth potential was of 11%, which is considerable.

Region-of-interest-based analysis with application of cortical thickness variation of left planum temporale in schizophrenia and psychotic bipolar disorder

In neuroimaging studies, spatial normalization and multivariate testing are central problems in characterizing group variation of functions (e.g., cortical thickness, curvature, functional response) in an atlas coordinate system across clinical populations. We present a region-of-interest (ROI)-based analysis framework for detecting such a group variation. This framework includes two main techniques: ROI-based registration via large deformation diffeomorphic metric surface mapping and a multivariate testing using a Gaussian random field (GRF) model on the cortical surface constructed by the eigenfunctions of the Laplace-Beltramioperator. We compared our GRF statistical model with a pointwise hypothesis testing approach, whose P-value is corrected using false discovery rate or random field theory at several smoothness scales. As an illustration, we applied this framework to a clinical study of the cortical thickness of the left planum temporale (PT) in subjects with psychotic bipolar disorder, schizophrenia, and healthy comparison controls. Our results show that the anterior portion of the left PT is thinner in the psychotic bipolar and schizophrenic groups than in the healthy control group, and the posterior portion of the left PT shows the reversal finding. Moreover, there may be a greater thickness variation in the left PT in psychotic bipolar patients when compared with that in schizophrenic patients.

Anisotropy induced by macroscopic boundaries: surface-normal mapping using diffusion-weighted imaging

In MRI, macroscopic boundaries lead to a diffusion-related increase in signal intensity near them--an effect commonly referred to as edge-enhancement. In diffusion-weighted imaging protocols where the signal attenuation due to diffusion results predominantly from the application of magnetic field gradients, edge-enhancement will depend on the orientation of these diffusion gradients. The resulting diffusion anisotropy can be exploited to map the direction normal to the macroscopic boundary. Simulations suggest that the hypothesized anisotropy may be within observable limits even when the voxel contains no boundary itself--hence, the name remote-anisotropy. Moreover, for certain experimental parameters there may be significant phase cancellations within the voxel that may lead to an edge detraction effect. When this is avoided, the eigenvector corresponding to the smallest eigenvalue of the diffusion tensor obtained from diffusion-tensor imaging can be used to create surface-normal maps conveniently. Experiments performed on simple geometric constructs as well as real tissue demonstrate the feasibility of using the edge-enhancement mechanism to map orientations orthogonal to macroscopic surfaces, which may be used to assess the integrity of tissue and organ boundaries noninvasively.

Anisotropy induced by macroscopic boundaries: surface-normal mapping using diffusion-weighted imaging

In MRI, macroscopic boundaries lead to a diffusion-related increase in signal intensity near them--an effect commonly referred to as edge-enhancement. In diffusion-weighted imaging protocols where the signal attenuation due to diffusion results predominantly from the application of magnetic field gradients, edge-enhancement will depend on the orientation of these diffusion gradients. The resulting diffusion anisotropy can be exploited to map the direction normal to the macroscopic boundary. Simulations suggest that the hypothesized anisotropy may be within observable limits even when the voxel contains no boundary itself--hence, the name remote-anisotropy. Moreover, for certain experimental parameters there may be significant phase cancellations within the voxel that may lead to an edge detraction effect. When this is avoided, the eigenvector corresponding to the smallest eigenvalue of the diffusion tensor obtained from diffusion-tensor imaging can be used to create surface-normal maps conveniently. Experiments performed on simple geometric constructs as well as real tissue demonstrate the feasibility of using the edge-enhancement mechanism to map orientations orthogonal to macroscopic surfaces, which may be used to assess the integrity of tissue and organ boundaries noninvasively.

A hierarchical algorithm for MR brain image parcellation

We introduce an algorithm for segmenting brain magnetic resonance (MR) images into anatomical compartments such as the major tissue classes and neuro-anatomical structures of the gray matter. The algorithm is guided by prior information represented within a tree structure. The tree mirrors the hierarchy of anatomical structures and the subtrees correspond to limited segmentation problems. The solution to each problem is estimated via a conventional classifier. Our algorithm can be adapted to a wide range of segmentation problems by modifying the tree structure or replacing the classifier. We evaluate the performance of our new segmentation approach by revisiting a previously published statistical group comparison between first-episode schizophrenia patients, first-episode affective psychosis patients, and comparison subjects. The original study is based on 50 MR volumes in which an expert identified the brain tissue classes as well as the superior temporal gyrus, amygdala, and hippocampus. We generate analogous segmentations using our new method and repeat the statistical group comparison. The results of our analysis are similar to the original findings, except for one structure (the left superior temporal gyrus) in which a trend-level statistical significance (p = 0.07) was observed instead of statistical significance.

Combining anatomical manifold information via diffeomorphic metric mappings for studying cortical thinning of the cingulate gyrus in schizophrenia

Spatial normalization is a crucial step in assessing patterns of neuroanatomical structure and function associated with health and disease. Errors that occur during spatial normalization can influence hypothesis testing due to the dimensionalities of mapping algorithms and anatomical manifolds (landmarks, curves, surfaces, volumes) used to drive the mapping algorithms. The primary aim of this paper is to improve statistical inference using multiple anatomical manifolds and large deformation diffeomorphic metric mapping (LDDMM) algorithms. We propose that combining information generated by the various manifolds and algorithms improves the reliability of hypothesis testing. We used this unified approach to assess variation in the thickness of the cingulate gyrus in subjects with schizophrenia and healthy comparison subjects. Three different LDDMM algorithms for mapping landmarks, curves and triangulated meshes were used to transform thickness maps of the cingulate surfaces into an atlas coordinate system. We then tested for group differences by combining the information from the three types of anatomical manifolds and LDDMM mapping algorithms. The unified approach provided reliable statistical results and eliminated ambiguous results due to surface mismatches. Subjects with schizophrenia had non-uniform cortical thinning over the left and right cingulate gyri, especially in the anterior portion, as compared to healthy comparison subjects.

Brain functional localization: a survey of image registration techniques

Functional localization is a concept which involves the application of a sequence of geometrical and statistical image processing operations in order to define the location of brain activity or to produce functional/parametric maps with respect to the brain structure or anatomy. Considering that functional brain images do not normally convey detailed structural information and, thus, do not present an anatomically specific localization of functional activity, various image registration techniques are introduced in the literature for the purpose of mapping functional activity into an anatomical image or a brain atlas. The problems addressed by these techniques differ depending on the application and the type of analysis, i.e., single-subject versus group analysis. Functional to anatomical brain image registration is the core part of functional localization in most applications and is accompanied by intersubject and subject-to-atlas registration for group analysis studies. Cortical surface registration and automatic brain labeling are some of the other tools towards establishing a fully automatic functional localization procedure. While several previous survey papers have reviewed and classified general-purpose medical image registration techniques, this paper provides an overview of brain functional localization along with a survey and classification of the image registration techniques related to this problem.

Diffeomorphic metric surface mapping in subregion of the superior temporal gyrus

This paper describes the application of large deformation diffeomorphic metric mapping to cortical surfaces based on the shape and geometric properties of subregions of the superior temporal gyrus in the human brain. The anatomical surfaces of the cortex are represented as triangulated meshes. The diffeomorphic matching algorithm is implemented by defining a norm between the triangulated meshes, based on the algorithms of Vaillant and Glaunès. The diffeomorphic correspondence is defined as a flow of the extrinsic three dimensional coordinates containing the cortical surface that registers the initial and target geometry by minimizing the norm. The methods are demonstrated in 40 high-resolution MRI cortical surfaces of planum temporale (PT) constructed from subsets of the superior temporal gyrus (STG). The effectiveness of the algorithm is demonstrated via the Euclidean positional distance, distance of normal vectors, and curvature before and after the surface matching as well as the comparison with a landmark matching algorithm. The results demonstrate that both the positional and shape variability of the anatomical configurations are being represented by the diffeomorphic maps.

MR imaging anatomy in neurodegeneration: a robust volumetric parcellation method of the frontal lobe gyri with quantitative validation in patients with dementia

BACKGROUND AND PURPOSE

Brain volumetry is widely used for evaluating tissue degeneration; however, the parcellation methods are rarely validated and use arbitrary planes to mark boundaries of brain regions. The goal of this study was to develop, validate, and apply an MR imaging tracing method for the parcellation of 3 major gyri of the frontal lobe, which uses only local landmarks intrinsic to the structures of interest, without the need for global reorientation or the use of dividing planes or lines.

METHODS

Studies were performed on 25 subjects--healthy controls and subjects diagnosed with Lewy body dementia and Alzheimer disease--with significant variation in the underlying gyral anatomy and state of atrophy. The protocol was evaluated by using multiple observers tracing scans of subjects diagnosed with neurodegenerative disease and those aging normally, and the results were compared by spatial overlap agreement. To confirm the results, observers marked the same locations in different brains. We illustrated the variabilities of the key boundaries that pose the greatest challenge to defining consistent parcellations across subjects.

RESULTS

The resulting gyral volumes were evaluated, and their consistency across raters was used as an additional assessment of the validity of our marking method. The agreement on a scale of 0-1 was found to be 0.83 spatial and 0.90 volumetric for the same rater and 0.85 spatial and 0.90 volumetric for 2 different raters. The results revealed that the protocol remained consistent across different neurodegenerative conditions.

CONCLUSION

Our method provides a simple and reliable way for the volumetric evaluation of frontal lobe neurodegeneration and can be used as a resource for larger comparative studies as well as a validation procedure of automated algorithms.

Mapping heritability and molecular genetic associations with cortical features using probabilistic brain atlases: methods and applications to schizophrenia

There is an urgent need to decipher the complex nature of genotype-phenotype relationships within the multiple dimensions of brain structure and function that are compromised in neuropsychiatric syndromes such as schizophrenia. Doing so requires sophisticated methodologies to represent population variability in neural traits and to probe their heritable and molecular genetic bases. We have recently developed and applied computational algorithms to map the heritability of, as well as genetic linkage and association to, neural features encoded using brain imaging in the context of three-dimensional (3D), populationbased, statistical brain atlases. One set of algorithms builds on our prior work using classical twin study methods to estimate heritability by fitting biometrical models for additive genetic, unique, and common environmental influences. Another set of algorithms performs regression-based (Haseman-Elston) identical-bydescent linkage analysis and genetic association analysis of DNA polymorphisms in relation to neural traits of interest in the same 3D population-based brain atlas format. We demonstrate these approaches using samples of healthy monozygotic (MZ) and dizygotic (DZ) twin pairs, as well as MZ and DZ twin pairs discordant for schizophrenia, but the methods can be generalized to other classes of relatives and to other diseases. The results confirm prior evidence of genetic influences on gray matter density in frontal brain regions. They also provide converging evidence that the chromosome 1q42 region is relevant to schizophrenia by demonstrating linkage and association of markers of the Transelin-Associated-Factor-X and Disrupted-In- Schizophrenia-1 genes with prefrontal cortical gray matter deficits in twins discordant for schizophrenia.

Computational anatomy: shape, growth, and atrophy comparison via diffeomorphisms

Computational anatomy (CA) is the mathematical study of anatomy I in I = I(alpha) o G, an orbit under groups of diffeomorphisms (i.e., smooth invertible mappings) g in G of anatomical exemplars I(alpha) in I. The observable images are the output of medical imaging devices. There are three components that CA examines: (i) constructions of the anatomical submanifolds, (ii) comparison of the anatomical manifolds via estimation of the underlying diffeomorphisms g in G defining the shape or geometry of the anatomical manifolds, and (iii) generation of probability laws of anatomical variation P(.) on the images I for inference and disease testing within anatomical models. This paper reviews recent advances in these three areas applied to shape, growth, and atrophy.

Validation of semiautomated methods for quantifying cingulate cortical metrics in schizophrenia

This paper validates semiautomated methods for reconstructing cortical surfaces of the cingulate gyrus from high-resolution magnetic resonance (MR) images. Bayesian segmentation was used to delineate the image voxels into five tissue types: cerebrospinal fluid (CSF), gray matter (GM), white matter (WM), and partial volumes of CSF/GM and GM/WM; the tissues were then recalibrated as CSF, GM, and WM via the Neyman-Pearson Likelihood Ratio Test. To generate cortical surfaces at the interface of GM and WM, the thresholds between the tissue types were first used to reassign partial volume voxels to CSF, GM, and WM with minimum error (that varied from 0.06 to 0.15 for the 10 subjects). Next, topology-correct cortical surfaces were generated and validated with almost all surface vertices lying within one voxel (0.5 mm) of hand contours. Dynamic programming was used to delineate and extract the cingulate gyrus from the cortical surfaces based on its gyral and sulcal boundaries. The intraclass correlation coefficient for surface area obtained by two raters for all 10 surfaces was 0.82. In addition, by repeating the entire procedure three times in one subject, we obtained a coefficient of variation of 0.0438 for surface area.