Bilateral focused ultrasound medial thalamotomies for trigeminal neuropathic pain: a randomized controlled study

OBJECTIVE

Medial thalamotomy has been shown to benefit patients with neuropathic pain, but widespread adoption of this procedure has been limited by reporting of clinical outcomes in studies without a control group. This study aimed to minimize confounders associated with medial thalamotomy for treating chronic pain by using modern MRI-guided stereotactic lesioning and a rigorous clinical design.

METHODS

This prospective, double-blinded, randomized controlled trial in 10 patients with trigeminal neuropathic pain used sham procedures as controls. Participants underwent assessments by a pain psychologist and pain management clinician, including use of the following measures: the Numeric Pain Rating Scale (NPRS); patient-reported outcome measures; and patient's impression of improvement at baseline, 1 day, 1 week, 1 month, and 3 months postprocedure. Patients in the treated group underwent bilateral focused ultrasound (FUS) medial thalamotomy targeting the central lateral nucleus. Patients in the control group underwent sham procedures with energy output disabled. The primary efficacy outcome measure was between-group differences in pain intensity (using the NPRS) at baseline and at 3 months postprocedure. Adverse events were measured for safety and included MRI analysis. Exploratory measures of connectivity and metabolism were analyzed using diffusion tensor imaging, functional MRI, and PET, respectively.

RESULTS

There were no serious complications from the FUS procedures. MRI confirmed bilateral medial thalamic ablations. There was no significant improvement in pain intensity from baseline to 3 months, either for patients undergoing FUS medial thalamotomy or for sham controls; and the between-group change in NPRS score as the primary efficacy outcome measure was not significantly different. Patient-reported outcome assessments demonstrated improvement (i.e., a decrease) only in pain interference with enjoyment of life at 3 months. There was a perception of benefit at 1 week, but only for patients treated with FUS and not for the sham cohort. Advanced neuroimaging showed that these medial thalamic lesions altered structural connectivity with the postcentral gyrus and demonstrated a trend toward hypometabolism in the insula and amygdala.

CONCLUSIONS

This randomized controlled trial of bilateral FUS medial thalamotomy did not reduce the intensity of trigeminal neuropathic pain, although it should be noted that the ability to estimate the magnitude of treatment effects is limited by the small cohort.

Autism Classification Using Topological Features and Deep Learning: A Cautionary Tale

The identification of autistic individuals using resting state functional connectivity networks can provide an objective diagnostic method for autism spectrum disorder (ASD). The present state-of-the-art machine learning model using deep learning has a classification accuracy of 70.2% on the ABIDE (Autism Brain Imaging Data Exchange) data set. In this paper, we explore the utility of topological features in the classification of ASD versus typically developing control subjects. These topological features have been shown to provide a complementary source of discriminative information in applications such as 2D object classification and social network analysis. We evaluate the performance of three different representations of topological features - persistence diagrams, persistence images, and persistence landscapes - for autism classification using neural networks, support vector machines and random forests. We also propose a hybrid approach of augmenting topological features with functional correlations, which typically outperforms the models that use functional correlations alone. With this approach, even with a simple 3-layer neural network, we are able to achieve a classification accuracy of 69.2% on the ABIDE data set. However, our experiments also show that the improvement due to topological features is not always statistically significant. Therefore, we offer a cautionary tale to the practitioners regarding the limited discriminative power of topological features derived from fMRI data for the classification of autism.

Image-based analysis and long-term clinical outcomes of deep brain stimulation for Tourette syndrome: a multisite study

BACKGROUND

Deep brain stimulation (DBS) can be an effective therapy for tics and comorbidities in select cases of severe, treatment-refractory Tourette syndrome (TS). Clinical responses remain variable across patients, which may be attributed to differences in the location of the neuroanatomical regions being stimulated. We evaluated active contact locations and regions of stimulation across a large cohort of patients with TS in an effort to guide future targeting.

METHODS

We collected retrospective clinical data and imaging from 13 international sites on 123 patients. We assessed the effects of DBS over time in 110 patients who were implanted in the centromedial (CM) thalamus (n=51), globus pallidus internus (GPi) (n=47), nucleus accumbens/anterior limb of the internal capsule (n=4) or a combination of targets (n=8). Contact locations (n=70 patients) and volumes of tissue activated (n=63 patients) were coregistered to create probabilistic stimulation atlases.

RESULTS

Tics and obsessive-compulsive behaviour (OCB) significantly improved over time (p<0.01), and there were no significant differences across brain targets (p>0.05). The median time was 13 months to reach a 40% improvement in tics, and there were no significant differences across targets (p=0.84), presence of OCB (p=0.09) or age at implantation (p=0.08). Active contacts were generally clustered near the target nuclei, with some variability that may reflect differences in targeting protocols, lead models and contact configurations. There were regions within and surrounding GPi and CM thalamus that improved tics for some patients but were ineffective for others. Regions within, superior or medial to GPi were associated with a greater improvement in OCB than regions inferior to GPi.

CONCLUSION

The results collectively indicate that DBS may improve tics and OCB, the effects may develop over several months, and stimulation locations relative to structural anatomy alone may not predict response. This study was the first to visualise and evaluate the regions of stimulation across a large cohort of patients with TS to generate new hypotheses about potential targets for improving tics and comorbidities.

Surface-Based Spatial Pyramid Matching of Cortical Regions for Analysis of Cognitive Performance

We propose a method to analyze the relationship between the shape of functional regions of the cortex and cognitive measures, such as reading ability and vocabulary knowledge. Functional regions on the cortical surface can vary not only in size and shape but also in topology and position relative to neighboring regions. Standard diffeomorphism-based shape analysis tools do not work well here because diffeomorphisms are unable to capture these topological differences, which include region splitting and merging across subjects. State-of-the-art cortical surface shape analyses compute derived regional properties (scalars), such as regional volume, cortical thickness, curvature, and gyrification index. However, these methods cannot compare the full extent of topological or shape differences in cortical regions. We propose icosahedral spatial pyramid matching (ISPM) of region borders computed on the surface of a sphere to capture this variation in regional topology, position, and shape. We then analyze how this variation corresponds to measures of cognitive performance. We compare our method to other approaches and find that it is indeed informative to consider aspects of shape beyond the standard approaches. Analysis is performed using a subset of 27 test/retest subjects from the Human Connectome Project in order to understand both the effectiveness and reproducibility of this method.

Revisiting Abnormalities in Brain Network Architecture Underlying Autism Using Topology-Inspired Statistical Inference

A large body of evidence relates autism with abnormal structural and functional brain connectivity. Structural covariance magnetic resonance imaging (scMRI) is a technique that maps brain regions with covarying gray matter densities across subjects. It provides a way to probe the anatomical structure underlying intrinsic connectivity networks (ICNs) through analysis of gray matter signal covariance. In this article, we apply topological data analysis in conjunction with scMRI to explore network-specific differences in the gray matter structure in subjects with autism versus age-, gender-, and IQ-matched controls. Specifically, we investigate topological differences in gray matter structure captured by structural correlation graphs derived from three ICNs strongly implicated in autism, namely the salience network, default mode network, and executive control network. By combining topological data analysis with statistical inference, our results provide evidence of statistically significant network-specific structural abnormalities in autism.

Riemannian Regression and Classification Models of Brain Networks Applied to Autism

Functional connectivity from resting-state functional MRI (rsfMRI) is typically represented as a symmetric positive definite (SPD) matrix. Analysis methods that exploit the Riemannian geometry of SPD matrices appropriately adhere to the positive definite constraint, unlike Euclidean methods. Recently proposed approaches for rsfMRI analysis have achieved high accuracy on public datasets, but are computationally intensive and difficult to interpret. In this paper, we show that we can get comparable results using connectivity matrices under the log-Euclidean and affine-invariant Riemannian metrics with relatively simple and interpretable models. On ABIDE Preprocessed dataset, our methods classify autism versus control subjects with 71.1% accuracy. We also show that Riemannian methods beat baseline in regressing connectome features to subject autism severity scores.

Nonparametric Aggregation of Geodesic Trends for Longitudinal Data Analysis

We propose a technique for analyzing longitudinal imaging data that models individual changes with diffeomorphic geodesic regression and aggregates these geodesics into a nonparametric group average trend. Our model is specifically tailored to the unbalanced and sparse characteristics of longitudinal imaging studies. That is, each individual has few data points measured over a short period of time, while the study population as a whole spans a wide age range. We use geodesic regression to estimate individual trends, which is an appropriate model for capturing shape changes over a short time window, as is typically found within an individual. Geodesics are also adept at handling the low sample sizes found within individuals, and can model the change between as few as two timepoints. However, geodesics are limited for modeling longer-term trends, where constant velocity may not be appropriate. Therefore, we develop a novel nonparametric regression to aggregate individual trends into an average group trend. We demonstrate the power of our method to capture non-geodesic group trends on hippocampal volume (real-valued data) and diffeomorphic registration of full 3D MRI from the longitudinal OASIS data.

Efficient Parallel Transport in the Group of Diffeomorphisms via Reduction to the Lie Algebra

This paper presents an efficient, numerically stable algorithm for parallel transport of tangent vectors in the group of diffeomorphisms. Previous approaches to parallel transport in large deformation diffeomorphic metric mapping (LDDMM) of images represent a momenta field, the dual of a tangent vector to the diffeomorphism group, as a scalar field times the image gradient. This "scalar momenta" constraint couples tangent vectors with the images being deformed and leads to computationally costly horizontal lifts in parallel transport. This paper uses the vector momenta formulation of LDDMM, which decouples the diffeomorphisms from the structures being transformed, e.g., images, point sets, etc. This decoupling leads to parallel transport expressed as a linear ODE in the Lie algebra. Solving this ODE directly is numerically stable and significantly faster than other LDDMM parallel transport methods. Results on 2D synthetic data and 3D brain MRI demonstrate that our algorithm is fast and conserves the inner products of the transported tangent vectors.

Multivariate characterization of white matter heterogeneity in autism spectrum disorder

The complexity and heterogeneity of neuroimaging findings in individuals with autism spectrum disorder has suggested that many of the underlying alterations are subtle and involve many brain regions and networks. The ability to account for multivariate brain features and identify neuroimaging measures that can be used to characterize individual variation have thus become increasingly important for interpreting and understanding the neurobiological mechanisms of autism. In the present study, we utilize the Mahalanobis distance, a multidimensional counterpart of the Euclidean distance, as an informative index to characterize individual brain variation and deviation in autism. Longitudinal diffusion tensor imaging data from 149 participants (92 diagnosed with autism spectrum disorder and 57 typically developing controls) between 3.1 and 36.83 years of age were acquired over a roughly 10-year period and used to construct the Mahalanobis distance from regional measures of white matter microstructure. Mahalanobis distances were significantly greater and more variable in the autistic individuals as compared to control participants, demonstrating increased atypicalities and variation in the group of individuals diagnosed with autism spectrum disorder. Distributions of multivariate measures were also found to provide greater discrimination and more sensitive delineation between autistic and typically developing individuals than conventional univariate measures, while also being significantly associated with observed traits of the autism group. These results help substantiate autism as a truly heterogeneous neurodevelopmental disorder, while also suggesting that collectively considering neuroimaging measures from multiple brain regions provides improved insight into the diversity of brain measures in autism that is not observed when considering the same regions separately. Distinguishing multidimensional brain relationships may thus be informative for identifying neuroimaging-based phenotypes, as well as help elucidate underlying neural mechanisms of brain variation in autism spectrum disorders.

Bayesian Covariate Selection in Mixed-Effects Models For Longitudinal Shape Analysis

The goal of longitudinal shape analysis is to understand how anatomical shape changes over time, in response to biological processes, including growth, aging, or disease. In many imaging studies, it is also critical to understand how these shape changes are affected by other factors, such as sex, disease diagnosis, IQ, etc. Current approaches to longitudinal shape analysis have focused on modeling age-related shape changes, but have not included the ability to handle covariates. In this paper, we present a novel Bayesian mixed-effects shape model that incorporates simultaneous relationships between longitudinal shape data and multiple predictors or covariates to the model. Moreover, we place an Automatic Relevance Determination (ARD) prior on the parameters, that lets us automatically select which covariates are most relevant to the model based on observed data. We evaluate our proposed model and inference procedure on a longitudinal study of Huntington's disease from PREDICT-HD. We first show the utility of the ARD prior for model selection in a univariate modeling of striatal volume, and next we apply the full high-dimensional longitudinal shape model to putamen shapes.

KERNEL PARTIAL LEAST SQUARES REGRESSION FOR RELATING FUNCTIONAL BRAIN NETWORK TOPOLOGY TO CLINICAL MEASURES OF BEHAVIOR

In this paper we present a novel method for analyzing the relationship between functional brain networks and behavioral phenotypes. Drawing from topological data analysis, we first extract topological features using persistent homology from functional brain networks that are derived from correlations in resting-state fMRI. Rather than fixing a discrete network topology by thresholding the connectivity matrix, these topological features capture the network organization across all continuous threshold values. We then propose to use a kernel partial least squares (kPLS) regression to statistically quantify the relationship between these topological features and behavior measures. The kPLS also provides an elegant way to combine multiple image features by using linear combinations of multiple kernels. In our experiments we test the ability of our proposed brain network analysis to predict autism severity from rs-fMRI. We show that combining correlations with topological features gives better prediction of autism severity than using correlations alone.

Investigating the Microstructural Correlation of White Matter in Autism Spectrum Disorder

White matter microstructure forms a complex and dynamical system that is critical for efficient and synchronized brain function. Neuroimaging findings in children with autism spectrum disorder (ASD) suggest this condition is associated with altered white matter microstructure, which may lead to atypical macroscale brain connectivity. In this study, we used diffusion tensor imaging measures to examine the extent that white matter tracts are interrelated within ASD and typical development. We assessed the strength of inter-regional white matter correlations between typically developing and ASD diagnosed individuals. Using hierarchical clustering analysis, clustering patterns of the pairwise white matter correlations were constructed and revealed to be different between the two groups. Additionally, we explored the use of graph theory analysis to examine the characteristics of the patterns formed by inter-regional white matter correlations and compared these properties between ASD and typical development. We demonstrate that the ASD sample has significantly less coherence in white matter microstructure across the brain compared to that in the typical development sample. The ASD group also presented altered topological characteristics, which may implicate less efficient brain networking in ASD. These findings highlight the potential of graph theory based network characteristics to describe the underlying networks as measured by diffusion magnetic resonance imaging and furthermore indicates that ASD may be associated with altered brain network characteristics. Our findings are consistent with those of a growing number of studies and hypotheses that have suggested disrupted brain connectivity in ASD.

Bayesian principal geodesic analysis for estimating intrinsic diffeomorphic image variability

In this paper, we present a generative Bayesian approach for estimating the low-dimensional latent space of diffeomorphic shape variability in a population of images. We develop a latent variable model for principal geodesic analysis (PGA) that provides a probabilistic framework for factor analysis in the space of diffeomorphisms. A sparsity prior in the model results in automatic selection of the number of relevant dimensions by driving unnecessary principal geodesics to zero. To infer model parameters, including the image atlas, principal geodesic deformations, and the effective dimensionality, we introduce an expectation maximization (EM) algorithm. We evaluate our proposed model on 2D synthetic data and the 3D OASIS brain database of magnetic resonance images, and show that the automatically selected latent dimensions from our model are able to reconstruct unobserved testing images with lower error than both linear principal component analysis (LPCA) in the image space and tangent space principal component analysis (TPCA) in the diffeomorphism space.

A Transformation Similarity Constraint for Groupwise Nonlinear Registration in Longitudinal Neuro Imaging Studies

Patients with Alzheimer's disease and other brain disorders often show a similar spatial distribution of volume change throughout the brain over time,^1,2 but this information is not yet used in registration algorithms to refine the quantification of change. Here, we develop a mathematical basis to incorporate that prior information into a longitudinal structural neuroimaging study. We modify the canonical minimization problem for non-linear registration to include a term that couples a collection of registrations together to enforce group similarity. More specifically, throughout the computation we maintain a group-level representation of the transformations and constrain updates to individual transformations to be similar to this representation. The derivations necessary to produce the Euler-Lagrange equations for the coupling term are presented and a gradient descent algorithm based on the formulation was implemented. We demonstrate using 57 longitudinal image pairs from the Alzheimer's Disease Neuroimaging Initiative (ADNI) that longitudinal registration with such a groupwise coupling prior is more robust to noise in estimating change, suggesting such change maps may have several important applications.

Atypical development of white matter microstructure of the corpus callosum in males with autism: a longitudinal investigation

BACKGROUND

The corpus callosum is the largest white matter structure in the brain, and it is the most consistently reported to be atypical in diffusion tensor imaging studies of autism spectrum disorder. In individuals with typical development, the corpus callosum is known to undergo a protracted development from childhood through young adulthood. However, no study has longitudinally examined the developmental trajectory of corpus callosum in autism past early childhood.

METHODS

The present study used a cohort sequential design over 9 years to examine age-related changes of the corpus callosum in 100 males with autism and 56 age-matched males with typical development from early childhood (when autism can first be reliably diagnosed) to mid-adulthood (after development of the corpus callosum has been completed) (3 to 41 years of age).

RESULTS

The group with autism demonstrated a different developmental trajectory of white matter microstructure in the anterior corpus callosum's (genu and body) fractional anisotropy, which suggests atypical brain maturation in these regions in autism. When analyses were broken down by age group, atypical developmental trajectories were present only in the youngest participants (10 years of age and younger). Significant main effects for group were found in terms of decreased fractional anisotropy across all three subregions of the corpus callosum (genu, body, and splenium) and increased mean diffusivity, radial diffusivity, and axial diffusivity in the posterior corpus callosum.

CONCLUSIONS

These longitudinal results suggest atypical early childhood development of the corpus callosum microstructure in autism that transitions into sustained group differences in adolescence and adulthood. This pattern of results provides longitudinal evidence consistent with a growing number of published studies and hypotheses regarding abnormal brain connectivity across the life span in autism.

Simultaneous Longitudinal Registration with Group-Wise Similarity Prior

Here we present an algorithm for the simultaneous registration of N longitudinal image pairs such that information acquired by each pair is used to constrain the registration of each other pair. More specifically, in the geodesic shooting setting for Large Deformation Diffeomorphic Metric Mappings (LDDMM) an average of the initial momenta characterizing the N transformations is maintained throughout and updates to individual momenta are constrained to be similar to this average. In this way, the N registrations are coupled and explore the space of diffeomorphisms as a group, the variance of which is constrained to be small. Our approach is motivated by the observation that transformations learned from images in the same diagnostic category share characteristics. The group-wise consistency prior serves to strengthen the contribution of the common signal among the N image pairs to the transformation for a specific pair, relative to features particular to that pair. We tested the algorithm on 57 longitudinal image pairs of Alzheimer's Disease patients from the Alzheimer's Disease Neuroimaging Initiative and evaluated the ability of the algorithm to produce momenta that better represent the long-term biological processes occurring in the underlying anatomy. We found that for many image pairs, momenta learned with the group-wise prior better predict a third time point image unobserved in the registration.

A Riemannian Framework for Intrinsic Comparison of Closed Genus-Zero Shapes

We present a framework for intrinsic comparison of surface metric structures and curvatures. This work parallels the work of Kurtek et al. on parameterization-invariant comparison of genus zero shapes. Here, instead of comparing the embedding of spherically parameterized surfaces in space, we focus on the first fundamental form. To ensure that the distance on spherical metric tensor fields is invariant to parameterization, we apply the conjugation-invariant metric arising from the L2 norm on symmetric positive definite matrices. As a reparameterization changes the metric tensor by a congruent Jacobian transform, this metric perfectly suits our purpose. The result is an intrinsic comparison of shape metric structure that does not depend on the specifics of a spherical mapping. Further, when restricted to tensors of fixed volume form, the manifold of metric tensor fields and its quotient of the group of unitary diffeomorphisms becomes a proper metric manifold that is geodesically complete. Exploiting this fact, and augmenting the metric with analogous metrics on curvatures, we derive a complete Riemannian framework for shape comparison and reconstruction. A by-product of our framework is a near-isometric and curvature-preserving mapping between surfaces. The correspondence is optimized using the fast spherical fluid algorithm. We validate our framework using several subcortical boundary surface models from the ADNI dataset.

Longitudinal volumetric brain changes in autism spectrum disorder ages 6-35 years

Since the impairments associated with autism spectrum disorder (ASD) tend to persist or worsen from childhood into adulthood, it is of critical importance to examine how the brain develops over this growth epoch. We report initial findings on whole and regional longitudinal brain development in 100 male participants with ASD (226 high-quality magnetic resonance imaging [MRI] scans; mean inter-scan interval 2.7 years) compared to 56 typically developing controls (TDCs) (117 high-quality scans; mean inter-scan interval 2.6 years) from childhood into adulthood, for a total of 156 participants scanned over an 8-year period. This initial analysis includes between one and three high-quality scans per participant that have been processed and segmented to date, with 21% having one scan, 27% with two scans, and 52% with three scans in the ASD sample; corresponding percentages for the TDC sample are 30%, 30%, and 40%. The proportion of participants with multiple scans (79% of ASDs and 68% of TDCs) was high in comparison to that of large longitudinal neuroimaging studies of typical development. We provide volumetric growth curves for the entire brain, total gray matter (GM), frontal GM, temporal GM, parietal GM, occipital GM, total cortical white matter (WM), corpus callosum, caudate, thalamus, total cerebellum, and total ventricles. Mean volume of cortical WM was reduced significantly. Mean ventricular volume was increased in the ASD sample relative to the TDCs across the broad age range studied. Decreases in regional mean volumes in the ASD sample most often were due to decreases during late adolescence and adulthood. The growth curve of whole brain volume over time showed increased volumes in young children with autism, and subsequently decreased during adolescence to meet the TDC curve between 10 and 15 years of age. The volume of many structures continued to decline atypically into adulthood in the ASD sample. The data suggest that ASD is a dynamic disorder with complex changes in whole and regional brain volumes that change over time from childhood into adulthood.

PARAMETRIC REGRESSION SCHEME FOR DISTRIBUTIONS: ANALYSIS OF DTI FIBER TRACT DIFFUSION CHANGES IN EARLY BRAIN DEVELOPMENT

Temporal modeling frameworks often operate on scalar variables by summarizing data at initial stages as statistical summaries of the underlying distributions. For instance, DTI analysis often employs summary statistics, like mean, for regions of interest and properties along fiber tracts for population studies and hypothesis testing. This reduction via discarding of variability information may introduce significant errors which propagate through the procedures. We propose a novel framework which uses distribution-valued variables to retain and utilize the local variability information. Classic linear regression is adapted to employ these variables for model estimation. The increased stability and reliability of our proposed method when compared with regression using single-valued statistical summaries, is demonstrated in a validation experiment with synthetic data. Our driving application is the modeling of age-related changes along DTI white matter tracts. Results are shown for the spatiotemporal population trajectory of genu tract estimated from 45 healthy infants and compared with a Krabbe's patient.

A functional network estimation method of resting-state fMRI using a hierarchical Markov random field

We propose a hierarchical Markov random field model for estimating both group and subject functional networks simultaneously. The model takes into account the within-subject spatial coherence as well as the between-subject consistency of the network label maps. The statistical dependency between group and subject networks acts as a regularization, which helps the network estimation on both layers. We use Gibbs sampling to approximate the posterior density of the network labels and Monte Carlo expectation maximization to estimate the model parameters. We compare our method with two alternative segmentation methods based on K-Means and normalized cuts, using synthetic and real fMRI data. The experimental results show that our proposed model is able to identify both group and subject functional networks with higher accuracy on synthetic data, more robustness, and inter-session consistency on the real data.

Longitudinal changes in cortical thickness in autism and typical development

The natural history of brain growth in autism spectrum disorders remains unclear. Cross-sectional studies have identified regional abnormalities in brain volume and cortical thickness in autism, although substantial discrepancies have been reported. Preliminary longitudinal studies using two time points and small samples have identified specific regional differences in cortical thickness in the disorder. To clarify age-related trajectories of cortical development, we examined longitudinal changes in cortical thickness within a large mixed cross-sectional and longitudinal sample of autistic subjects and age- and gender-matched typically developing controls. Three hundred and forty-five magnetic resonance imaging scans were examined from 97 males with autism (mean age = 16.8 years; range 3-36 years) and 60 males with typical development (mean age = 18 years; range 4-39 years), with an average interscan interval of 2.6 years. FreeSurfer image analysis software was used to parcellate the cortex into 34 regions of interest per hemisphere and to calculate mean cortical thickness for each region. Longitudinal linear mixed effects models were used to further characterize these findings and identify regions with between-group differences in longitudinal age-related trajectories. Using mean age at time of first scan as a reference (15 years), differences were observed in bilateral inferior frontal gyrus, pars opercularis and pars triangularis, right caudal middle frontal and left rostral middle frontal regions, and left frontal pole. However, group differences in cortical thickness varied by developmental stage, and were influenced by IQ. Differences in age-related trajectories emerged in bilateral parietal and occipital regions (postcentral gyrus, cuneus, lingual gyrus, pericalcarine cortex), left frontal regions (pars opercularis, rostral middle frontal and frontal pole), left supramarginal gyrus, and right transverse temporal gyrus, superior parietal lobule, and paracentral, lateral orbitofrontal, and lateral occipital regions. We suggest that abnormal cortical development in autism spectrum disorders undergoes three distinct phases: accelerated expansion in early childhood, accelerated thinning in later childhood and adolescence, and decelerated thinning in early adulthood. Moreover, cortical thickness abnormalities in autism spectrum disorders are region-specific, vary with age, and may remain dynamic well into adulthood.

4D ACTIVE CUT: AN INTERACTIVE TOOL FOR PATHOLOGICAL ANATOMY MODELING

4D pathological anatomy modeling is key to understanding complex pathological brain images. It is a challenging problem due to the difficulties in detecting multiple appearing and disappearing lesions across time points and estimating dynamic changes and deformations between them. We propose a novel semi-supervised method, called 4D active cut, for lesion recognition and deformation estimation. Existing interactive segmentation methods passively wait for user to refine the segmentations which is a difficult task in 3D images that change over time. 4D active cut instead actively selects candidate regions for querying the user, and obtains the most informative user feedback. A user simply answers 'yes' or 'no' to a candidate object without having to refine the segmentation slice by slice. Compared to single-object detection of the existing methods, our method also detects multiple lesions with spatial coherence using Markov random fields constraints. Results show improvement on the lesion detection, which subsequently improves deformation estimation.

Quantifying anatomical shape variations in neurological disorders

We develop a multivariate analysis of brain anatomy to identify the relevant shape deformation patterns and quantify the shape changes that explain corresponding variations in clinical neuropsychological measures. We use kernel Partial Least Squares (PLS) and formulate a regression model in the tangent space of the manifold of diffeomorphisms characterized by deformation momenta. The scalar deformation momenta completely encode the diffeomorphic changes in anatomical shape. In this model, the clinical measures are the response variables, while the anatomical variability is treated as the independent variable. To better understand the "shape-clinical response" relationship, we also control for demographic confounders, such as age, gender, and years of education in our regression model. We evaluate the proposed methodology on the Alzheimer's Disease Neuroimaging Initiative (ADNI) database using baseline structural MR imaging data and neuropsychological evaluation test scores. We demonstrate the ability of our model to quantify the anatomical deformations in units of clinical response. Our results also demonstrate that the proposed method is generic and generates reliable shape deformations both in terms of the extracted patterns and the amount of shape changes. We found that while the hippocampus and amygdala emerge as mainly responsible for changes in test scores for global measures of dementia and memory function, they are not a determinant factor for executive function. Another critical finding was the appearance of thalamus and putamen as most important regions that relate to executive function. These resulting anatomical regions were consistent with very high confidence irrespective of the size of the population used in the study. This data-driven global analysis of brain anatomy was able to reach similar conclusions as other studies in Alzheimer's disease based on predefined ROIs, together with the identification of other new patterns of deformation. The proposed methodology thus holds promise for discovering new patterns of shape changes in the human brain that could add to our understanding of disease progression in neurological disorders.

Abnormal lateralization of functional connectivity between language and default mode regions in autism

Background

Lateralization of brain structure and function occurs in typical development, and abnormal lateralization is present in various neuropsychiatric disorders. Autism is characterized by a lack of left lateralization in structure and function of regions involved in language, such as Broca and Wernicke areas.

Methods

Using functional connectivity magnetic resonance imaging from a large publicly available sample (n = 964), we tested whether abnormal functional lateralization in autism exists preferentially in language regions or in a more diffuse pattern across networks of lateralized brain regions.

Results

The autism group exhibited significantly reduced left lateralization in a few connections involving language regions and regions from the default mode network, but results were not significant throughout left- and right-lateralized networks. There is a trend that suggests the lack of left lateralization in a connection involving Wernicke area and the posterior cingulate cortex associates with more severe autism.

Conclusions

Abnormal language lateralization in autism may be due to abnormal language development rather than to a deficit in hemispheric specialization of the entire brain.

Diffeomorphic shape trajectories for improved longitudinal segmentation and statistics

Longitudinal imaging studies involve tracking changes in individuals by repeated image acquisition over time. The goal of these studies is to quantify biological shape variability within and across individuals, and also to distinguish between normal and disease populations. However, data variability is influenced by outside sources such as image acquisition, image calibration, human expert judgment, and limited robustness of segmentation and registration algorithms. In this paper, we propose a two-stage method for the statistical analysis of longitudinal shape. In the first stage, we estimate diffeomorphic shape trajectories for each individual that minimize inconsistencies in segmented shapes across time. This is followed by a longitudinal mixed-effects statistical model in the second stage for testing differences in shape trajectories between groups. We apply our method to a longitudinal database from PREDICT-HD and demonstrate our approach reduces unwanted variability for both shape and derived measures, such as volume. This leads to greater statistical power to distinguish differences in shape trajectory between healthy subjects and subjects with a genetic biomarker for Huntington's disease (HD).

Improved segmentation of white matter tracts with adaptive Riemannian metrics

We present a novel geodesic approach to segmentation of white matter tracts from diffusion tensor imaging (DTI). Compared to deterministic and stochastic tractography, geodesic approaches treat the geometry of the brain white matter as a manifold, often using the inverse tensor field as a Riemannian metric. The white matter pathways are then inferred from the resulting geodesics, which have the desirable property that they tend to follow the main eigenvectors of the tensors, yet still have the flexibility to deviate from these directions when it results in lower costs. While this makes such methods more robust to noise, the choice of Riemannian metric in these methods is ad hoc. A serious drawback of current geodesic methods is that geodesics tend to deviate from the major eigenvectors in high-curvature areas in order to achieve the shortest path. In this paper we propose a method for learning an adaptive Riemannian metric from the DTI data, where the resulting geodesics more closely follow the principal eigenvector of the diffusion tensors even in high-curvature regions. We also develop a way to automatically segment the white matter tracts based on the computed geodesics. We show the robustness of our method on simulated data with different noise levels. We also compare our method with tractography methods and geodesic approaches using other Riemannian metrics and demonstrate that the proposed method results in improved geodesics and segmentations using both synthetic and real DTI data.

Subject-specific prediction using nonlinear population modeling: application to early brain maturation from DTI

The term prediction implies expected outcome in the future, often based on a model and statistical inference. Longitudinal imaging studies offer the possibility to model temporal change trajectories of anatomy across populations of subjects. In the spirit of subject-specific analysis, such normative models can then be used to compare data from new subjects to the norm and to study progression of disease or to predict outcome. This paper follows a statistical inference approach and presents a framework for prediction of future observations based on past measurements and population statistics. We describe prediction in the context of nonlinear mixed effects modeling (NLME) where the full reference population's statistics (estimated fixed effects, variance-covariance of random effects, variance of noise) is used along with the individual's available observations to predict its trajectory. The proposed methodology is generic in regard to application domains. Here, we demonstrate analysis of early infant brain maturation from longitudinal DTI with up to three time points. Growth as observed in DTI-derived scalar invariants is modeled with a parametric function, its parameters being input to NLME population modeling. Trajectories of new subject's data are estimated when using no observation, only the first or the first two time points. Leave-one-out experiments result in statistics on differences between actual and predicted observations. We also simulate a clinical scenario of prediction on multiple categories, where trajectories predicted from multiple models are classified based on maximum likelihood criteria.

Hierarchical bayesian modeling, estimation, and sampling for multigroup shape analysis

This paper proposes a novel method for the analysis of anatomical shapes present in biomedical image data. Motivated by the natural organization of population data into multiple groups, this paper presents a novel hierarchical generative statistical model on shapes. The proposed method represents shapes using pointsets and defines a joint distribution on the population's (i) shape variables and (ii) object-boundary data. The proposed method solves for optimal (i) point locations, (ii) correspondences, and (iii) model-parameter values as a single optimization problem. The optimization uses expectation maximization relying on a novel Markov-chain Monte-Carlo algorithm for sampling in Kendall shape space. Results on clinical brain images demonstrate advantages over the state of the art.

SPATIOTEMPORAL MODELING OF DISCRETE-TIME DISTRIBUTION-VALUED DATA APPLIED TO DTI TRACT EVOLUTION IN INFANT NEURODEVELOPMENT

This paper proposes a novel method that extends spatiotemporal growth modeling to distribution-valued data. The method relaxes assumptions on the underlying noise models by considering the data to be represented by the complete probability distributions rather than a representative, single-valued summary statistics like the mean. When summarizing by the latter method, information on the underlying variability of data is lost early in the process and is not available at later stages of statistical analysis. The concept of 'distance' between distributions and an 'average' of distributions is employed. The framework quantifies growth trajectories for individuals and populations in terms of the complete data variability estimated along time and space. Concept is demonstrated in the context of our driving application which is modeling of age-related changes along white matter tracts in early neurodevelopment. Results are shown for a single subject with Krabbe's disease in comparison with a normative trend estimated from 15 healthy controls.

Multisite functional connectivity MRI classification of autism: ABIDE results

Background: Systematic differences in functional connectivity MRI metrics have been consistently observed in autism, with predominantly decreased cortico-cortical connectivity. Previous attempts at single subject classification in high-functioning autism using whole brain point-to-point functional connectivity have yielded about 80% accurate classification of autism vs. control subjects across a wide age range. We attempted to replicate the method and results using the Autism Brain Imaging Data Exchange (ABIDE) including resting state fMRI data obtained from 964 subjects and 16 separate international sites.

Methods: For each of 964 subjects, we obtained pairwise functional connectivity measurements from a lattice of 7266 regions of interest covering the gray matter (26.4 million “connections”) after preprocessing that included motion and slice timing correction, coregistration to an anatomic image, normalization to standard space, and voxelwise removal by regression of motion parameters, soft tissue, CSF, and white matter signals. Connections were grouped into multiple bins, and a leave-one-out classifier was evaluated on connections comprising each set of bins. Age, age-squared, gender, handedness, and site were included as covariates for the classifier.

Results: Classification accuracy significantly outperformed chance but was much lower for multisite prediction than for previous single site results. As high as 60% accuracy was obtained for whole brain classification, with the best accuracy from connections involving regions of the default mode network, parahippocampaland fusiform gyri, insula, Wernicke Area, and intraparietal sulcus. The classifier score was related to symptom severity, social function, daily living skills, and verbal IQ. Classification accuracy was significantly higher for sites with longer BOLD imaging times.

Conclusions: Multisite functional connectivity classification of autism outperformed chance using a simple leave-one-out classifier, but exhibited poorer accuracy than for single site results. Attempts to use multisite classifiers will likely require improved classification algorithms, longer BOLD imaging times, and standardized acquisition parameters for possible future clinical utility.

Neuropsychological investigation of motor impairments in autism

It is unclear how standardized neuropsychological measures of motor function relate to brain volumes of motor regions in autism spectrum disorder (ASD). An all-male sample composed of 59 ASD and 30 controls (ages 5-33 years) completed three measures of motor function: strength of grip (SOG), finger tapping test (FTT), and grooved pegboard test (GPT). Likewise, all participants underwent magnetic resonance imaging with region of interest (ROI) volumes obtained to include the following regions: motor cortex (precentral gyrus), somatosensory cortex (postcentral gyrus), thalamus, basal ganglia, cerebellum, and caudal middle frontal gyrus. These traditional neuropsychological measures of motor function are assumed to differ in motor complexity, with GPT requiring the most followed by FTT and SOG. Performance by ASD participants on the GPT and FTT differed significantly from that of controls, with the largest effect size differences observed on the more complex GPT task. Differences on the SOG task between the two groups were nonsignificant. Since more complex motor tasks tap more complex networks, poorer GPT performance by those with ASD may reflect less efficient motor networks. There was no gross pathology observed in classic motor areas of the brain in ASD, as ROI volumes did not differ, but FTT was negatively related to motor cortex volume in ASD. The results suggest a hierarchical motor disruption in ASD, with difficulties evident only in more complex tasks as well as a potential anomalous size-function relation in motor cortex in ASD.

A VECTOR MOMENTA FORMULATION OF DIFFEOMORPHISMS FOR IMPROVED GEODESIC REGRESSION AND ATLAS CONSTRUCTION

This paper presents a novel approach for diffeomorphic image regression and atlas estimation that results in improved convergence and numerical stability. We use a vector momenta representation of a diffeomorphism's initial conditions instead of the standard scalar momentum that is typically used. The corresponding variational problem results in a closed-form update for template estimation in both the geodesic regression and atlas estimation problems. While we show that the theoretical optimal solution is equivalent to the scalar momenta case, the simplification of the optimization problem leads to more stable and efficient estimation in practice. We demonstrate the effectiveness of our method for atlas estimation and geodesic regression using synthetically generated shapes and 3D MRI brain scans.

Longitudinal Heschl's gyrus growth during childhood and adolescence in typical development and autism

Heightened auditory sensitivity and atypical auditory processing are common in autism. Functional studies suggest abnormal neural response and hemispheric activation to auditory stimuli, yet the neurodevelopment underlying atypical auditory function in autism is unknown. In this study, we model longitudinal volumetric growth of Heschl's gyrus gray matter and white matter during childhood and adolescence in 40 individuals with autism and 17 typically developing participants. Up to three time points of magnetic resonance imaging data, collected on average every 2.5 years, were examined from individuals 3-12 years of age at the time of their first scan. Consistent with previous cross-sectional studies, no group differences were found in Heschl's gyrus gray matter volume or asymmetry. However, reduced longitudinal gray matter volumetric growth was found in the right Heschl's gyrus in autism. Reduced longitudinal white matter growth in the left hemisphere was found in the right-handed autism participants. Atypical Heschl's gyrus white matter volumetric growth was found bilaterally in the autism individuals with a history of delayed onset of spoken language. Heightened auditory sensitivity, obtained from the Sensory Profile, was associated with reduced volumetric gray matter growth in the right hemisphere. Our longitudinal analyses revealed dynamic gray and white matter changes in Heschl's gyrus throughout childhood and adolescence in both typical development and autism.

Voxel-wise co-analysis of macro- and microstructural brain alteration in mild cognitive impairment and Alzheimer's disease using anatomical and diffusion MRI

BACKGROUND AND PURPOSE

To determine if a voxel-wise "co-analysis" of structural and diffusion tensor magnetic resonance imaging (MRI) together reveals additional brain regions affected in mild cognitive impairment (MCI) and Alzheimer's disease (AD) than voxel-wise analysis of the individual MRI modalities alone.

METHODS

Twenty-one patients with MCI, 21 patients with AD, and 21 cognitively normal healthy elderly were studied with MRI. Maps of deformation and fractional anisotropy (FA) were computed and used as dependent variables in univariate and multivariate statistical models.

RESULTS

Univariate voxel-wise analysis of macrostructural changes in MCI showed atrophy in the right anterior temporal lobe, left posterior parietal/precuneus region, WM adjacent to the cingulate gyrus, and dorsolateral prefrontal regions, consistent with prior research. Univariate voxel-wise analysis of microstructural changes in MCI showed reduced FA in the left posterior parietal region extending into the corpus callosum, consistent with previous work. The multivariate analysis, which provides more information than univariate tests when structural and FA measures are correlated, revealed additional MCI-related changes in corpus callosum and temporal lobe.

CONCLUSION

These results suggest that in corpus callosum and temporal regions macro- and microstructural variations in MCI can be congruent, providing potentially new insight into the mechanisms of brain tissue degeneration.

MULTIVARIATE MODELING OF LONGITUDINAL MRI IN EARLY BRAIN DEVELOPMENT WITH CONFIDENCE MEASURES

The human brain undergoes rapid organization and structuring early in life. Longitudinal imaging enables the study of these changes over a developmental period within individuals through estimation of population growth trajectory and its variability. In this paper, we focus on maturation of white and gray matter depicted in structural and diffusion MRI of healthy subjects with repeated scans. We provide a framework for joint analysis of both structural MRI and DTI (Diffusion Tensor Imaging) using multivariate nonlinear mixed effect modeling of temporal changes. Our framework constructs normative growth models for all the modalities, taking into account the correlation among the modalities and individuals, along with estimation of the variability of the population trends. We apply our method to study early brain development, and to our knowledge this is the first multimodel longitudinal modeling of diffusion and signal intensity changes for this growth stage. Results show the potential of our framework to study growth trajectories, as well as neurodevelopmental disorders through comparison against the constructed normative models of multimodal 4D MRI.

A hierarchical geodesic model for diffeomorphic longitudinal shape analysis

Hierarchical linear models (HLMs) are a standard approach for analyzing data where individuals are measured repeatedly over time. However, such models are only applicable to longitudinal studies of Euclidean data. In this paper, we propose a novel hierarchical geodesic model (HGM), which generalizes HLMs to the manifold setting. Our proposed model explains the longitudinal trends in shapes represented as elements of the group of diffeomorphisms. The individual level geodesics represent the trajectory of shape changes within individuals. The group level geodesic represents the average trajectory of shape changes for the population. We derive the solution of HGMs on diffeomorphisms to estimate individual level geodesics, the group geodesic, and the residual geodesics. We demonstrate the effectiveness of HGMs for longitudinal analysis of synthetically generated shapes and 3D MRI brain scans.

Model selection and estimation of multi-compartment models in diffusion MRI with a Rician noise model

Multi-compartment models in diffusion MRI (dMRI) are used to describe complex white matter fiber architecture of the brain. In this paper, we propose a novel multi-compartment estimation method based on the ball-and-stick model, which is composed of an isotropic diffusion compartment ("ball") as well as one or more perfectly linear diffusion compartments ("sticks"). To model the noise distribution intrinsic to dMRI measurements, we introduce a Rician likelihood term and estimate the model parameters by means of an Expectation Maximization (EM) algorithm. This paper also addresses the problem of selecting the number of fiber compartments that best fit the data, by introducing a sparsity prior on the volume mixing fractions. This term provides automatic model selection and enables us to discriminate different fiber populations. When applied to simulated data, our method provides accurate estimates of the fiber orientations, diffusivities, and number of compartments, even at low SNR, and outperforms similar methods that rely on a Gaussian noise distribution assumption. We also apply our method to in vivo brain data and show that it can successfully capture complex fiber structures that match the known anatomy.

Regional characterization of longitudinal DT-MRI to study white matter maturation of the early developing brain

The human brain undergoes rapid and dynamic development early in life. Assessment of brain growth patterns relevant to neurological disorders and disease requires a normative population model of growth and variability in order to evaluate deviation from typical development. In this paper, we focus on maturation of brain white matter as shown in diffusion tensor MRI (DT-MRI), measured by fractional anisotropy (FA), mean diffusivity (MD), as well as axial and radial diffusivities (AD, RD). We present a novel methodology to model temporal changes of white matter diffusion from longitudinal DT-MRI data taken at discrete time points. Our proposed framework combines nonlinear modeling of trajectories of individual subjects, population analysis, and testing for regional differences in growth pattern. We first perform deformable mapping of longitudinal DT-MRI of healthy infants imaged at birth, 1 year, and 2 years of age, into a common unbiased atlas. An existing template of labeled white matter regions is registered to this atlas to define anatomical regions of interest. Diffusivity properties of these regions, presented over time, serve as input to the longitudinal characterization of changes. We use non-linear mixed effect (NLME) modeling where temporal change is described by the Gompertz function. The Gompertz growth function uses intuitive parameters related to delay, rate of change, and expected asymptotic value; all descriptive measures which can answer clinical questions related to quantitative analysis of growth patterns. Results suggest that our proposed framework provides descriptive and quantitative information on growth trajectories that can be interpreted by clinicians using natural language terms that describe growth. Statistical analysis of regional differences between anatomical regions which are known to mature differently demonstrates the potential of the proposed method for quantitative assessment of brain growth and differences thereof. This will eventually lead to a prediction of white matter diffusion properties and associated cognitive development at later stages given imaging data at early stages.

scMRI reveals large-scale brain network abnormalities in autism

Autism is a complex neurological condition characterized by childhood onset of dysfunction in multiple cognitive domains including socio-emotional function, speech and language, and processing of internally versus externally directed stimuli. Although gross brain anatomic differences in autism are well established, recent studies investigating regional differences in brain structure and function have yielded divergent and seemingly contradictory results. How regional abnormalities relate to the autistic phenotype remains unclear. We hypothesized that autism exhibits distinct perturbations in network-level brain architecture, and that cognitive dysfunction may be reflected by abnormal network structure. Network-level anatomic abnormalities in autism have not been previously described. We used structural covariance MRI to investigate network-level differences in gray matter structure within two large-scale networks strongly implicated in autism, the salience network and the default mode network, in autistic subjects and age-, gender-, and IQ-matched controls. We report specific perturbations in brain network architecture in the salience and default-mode networks consistent with clinical manifestations of autism. Extent and distribution of the salience network, involved in social-emotional regulation of environmental stimuli, is restricted in autism. In contrast, posterior elements of the default mode network have increased spatial distribution, suggesting a ‘posteriorization’ of this network. These findings are consistent with a network-based model of autism, and suggest a unifying interpretation of previous work. Moreover, we provide evidence of specific abnormalities in brain network architecture underlying autism that are quantifiable using standard clinical MRI.

Mixed-Effects Shape Models for Estimating Longitudinal Changes in Anatomy

In this paper, we propose a new method for longitudinal shape analysis that fits a linear mixed-effects model, while simultaneously optimizing correspondences on a set of anatomical shapes. Shape changes are modeled in a hierarchical fashion, with the global population trend as a fixed effect and individual trends as random effects. The statistical significance of the estimated trends are evaluated using specifically designed permutation tests. We also develop a permutation test based on the Hotelling T² statistic to compare the average shapes trends between two populations. We demonstrate the benefits of our method on a synthetic example of longitudinal tori and data from a developmental neuroimaging study.

Sasaki Metrics for Analysis of Longitudinal Data on Manifolds

Longitudinal data arises in many applications in which the goal is to understand changes in individual entities over time. In this paper, we present a method for analyzing longitudinal data that take values in a Riemannian manifold. A driving application is to characterize anatomical shape changes and to distinguish between trends in anatomy that are healthy versus those that are due to disease. We present a generative hierarchical model in which each individual is modeled by a geodesic trend, which in turn is considered as a perturbation of the mean geodesic trend for the population. Each geodesic in the model can be uniquely parameterized by a starting point and velocity, i.e., a point in the tangent bundle. Comparison between these parameters is achieved through the Sasaki metric, which provides a natural distance metric on the tangent bundle. We develop a statistical hypothesis test for differences between two groups of longitudinal data by generalizing the Hotelling T ² statistic to manifolds. We demonstrate the ability of these methods to distinguish differences in shape changes in a comparison of longitudinal corpus callosum data in subjects with dementia versus healthily aging controls.

STATISTICAL GROWTH MODELING OF LONGITUDINAL DT-MRI FOR REGIONAL CHARACTERIZATION OF EARLY BRAIN DEVELOPMENT

A population growth model that represents the growth trajectories of individual subjects is critical to study and understand neurodevelopment. This paper presents a framework for jointly estimating and modeling individual and population growth trajectories, and determining significant regional differences in growth pattern characteristics applied to longitudinal neuroimaging data. We use non-linear mixed effect modeling where temporal change is modeled by the Gompertz function. The Gompertz function uses intuitive parameters related to delay, rate of change, and expected asymptotic value; all descriptive measures which can answer clinical questions related to growth. Our proposed framework combines nonlinear modeling of individual trajectories, population analysis, and testing for regional differences. We apply this framework to the study of early maturation in white matter regions as measured with diffusion tensor imaging (DTI). Regional differences between anatomical regions of interest that are known to mature differently are analyzed and quantified. Experiments with image data from a large ongoing clinical study show that our framework provides descriptive, quantitative information on growth trajectories that can be directly interpreted by clinicians. To our knowledge, this is the first longitudinal analysis of growth functions to explain the trajectory of early brain maturation as it is represented in DTI.

Genetic, structural and functional imaging biomarkers for early detection of conversion from MCI to AD

With the advent of advanced imaging techniques, genotyping, and methods to assess clinical and biological progression, there is a growing need for a unified framework that could exploit information available from multiple sources to aid diagnosis and the identification of early signs of Alzheimer's disease (AD). We propose a modeling strategy using supervised feature extraction to optimally combine high-dimensional imaging modalities with several other low-dimensional disease risk factors. The motivation is to discover new imaging biomarkers and use them in conjunction with other known biomarkers for prognosis of individuals at high risk of developing AD. Our framework also has the ability to assess the relative importance of imaging modalities for predicting AD conversion. We evaluate the proposed methodology on the Alzheimer's Disease Neuroimaging Initiative (ADNI) database to predict conversion of individuals with mild cognitive impairment (MCI) to AD, only using information available at baseline.

Functional connectivity magnetic resonance imaging classification of autism

Group differences in resting state functional magnetic resonance imaging connectivity between individuals with autism and typically developing controls have been widely replicated for a small number of discrete brain regions, yet the whole-brain distribution of connectivity abnormalities in autism is not well characterized. It is also unclear whether functional connectivity is sufficiently robust to be used as a diagnostic or prognostic metric in individual patients with autism. We obtained pairwise functional connectivity measurements from a lattice of 7266 regions of interest covering the entire grey matter (26.4 million connections) in a well-characterized set of 40 male adolescents and young adults with autism and 40 age-, sex- and IQ-matched typically developing subjects. A single resting state blood oxygen level-dependent scan of 8 min was used for the classification in each subject. A leave-one-out classifier successfully distinguished autism from control subjects with 83% sensitivity and 75% specificity for a total accuracy of 79% (P = 1.1 × 10(-7)). In subjects <20 years of age, the classifier performed at 89% accuracy (P = 5.4 × 10(-7)). In a replication dataset consisting of 21 individuals from six families with both affected and unaffected siblings, the classifier performed at 71% accuracy (91% accuracy for subjects <20 years of age). Classification scores in subjects with autism were significantly correlated with the Social Responsiveness Scale (P = 0.05), verbal IQ (P = 0.02) and the Autism Diagnostic Observation Schedule-Generic's combined social and communication subscores (P = 0.05). An analysis of informative connections demonstrated that region of interest pairs with strongest correlation values were most abnormal in autism. Negatively correlated region of interest pairs showed higher correlation in autism (less anticorrelation), possibly representing weaker inhibitory connections, particularly for long connections (Euclidean distance >10 cm). Brain regions showing greatest differences included regions of the default mode network, superior parietal lobule, fusiform gyrus and anterior insula. Overall, classification accuracy was better for younger subjects, with differences between autism and control subjects diminishing after 19 years of age. Classification scores of unaffected siblings of individuals with autism were more similar to those of the control subjects than to those of the subjects with autism. These findings indicate feasibility of a functional connectivity magnetic resonance imaging diagnostic assay for autism.

Adaptive Riemannian metrics for improved geodesic tracking of white matter

We present a new geodesic approach for studying white matter connectivity from diffusion tensor imaging (DTI). Previous approaches have used the inverse diffusion tensor field as a Riemannian metric and constructed white matter tracts as geodesics on the resulting manifold. These geodesics have the desirable property that they tend to follow the main eigenvectors of the tensors, yet still have the flexibility to deviate from these directions when it results in lower costs. While this makes such methods more robust to noise, it also has the serious drawback that geodesics tend to deviate from the major eigenvectors in high-curvature areas in order to achieve the shortest path. In this paper we formulate a modification of the Riemannian metric that results in geodesics adapted to follow the principal eigendirection of the tensor even in high-curvature regions. We show that this correction can be formulated as a simple scalar field modulation of the metric and that the appropriate variational problem results in a Poisson's equation on the Riemannian manifold. We demonstrate that the proposed method results in improved geodesics using both synthetic and real DTI data.

Atypical diffusion tensor hemispheric asymmetry in autism

BACKGROUND

Biological measurements that distinguish individuals with autism from typically developing individuals and those with other developmental and neuropsychiatric disorders must demonstrate very high performance to have clinical value as potential imaging biomarkers. We hypothesized that further study of white matter microstructure (WMM) in the superior temporal gyrus (STG) and temporal stem (TS), two brain regions in the temporal lobe containing circuitry central to language, emotion, and social cognition, would identify a useful combination of classification features and further understand autism neuropathology.

METHODS

WMM measurements from the STG and TS were examined from 30 high-functioning males satisfying full criteria for idiopathic autism aged 7-28 years and 30 matched controls and a replication sample of 12 males with idiopathic autism and 7 matched controls who participated in a previous case-control diffusion tensor imaging (DTI) study. Language functioning, adaptive functioning, and psychotropic medication usage were also examined.

RESULTS

In the STG, we find reversed hemispheric asymmetry of two separable measures of directional diffusion coherence, tensor skewness, and fractional anisotropy. In autism, tensor skewness is greater on the right and fractional anisotropy is decreased on the left. We also find increased diffusion parallel to white matter fibers bilaterally. In the right not left TS, we find increased omnidirectional, parallel, and perpendicular diffusion. These six multivariate measurements possess very high ability to discriminate individuals with autism from individuals without autism with 94% sensitivity, 90% specificity, and 92% accuracy in our original and replication samples. We also report a near-significant association between the classifier and a quantitative trait index of autism and significant correlations between two classifier components and measures of language, IQ, and adaptive functioning in autism.

Multivariate statistical analysis of deformation momenta relating anatomical shape to neuropsychological measures

The purpose of this study is to characterize the neuroanatomical variations observed in neurological disorders such as dementia. We do a global statistical analysis of brain anatomy and identify the relevant shape deformation patterns that explain corresponding variations in clinical neuropsychological measures. The motivation is to model the inherent relation between anatomical shape and clinical measures and evaluate its statistical significance. We use Partial Least Squares for the multivariate statistical analysis of the deformation momenta under the Large Deformation Diffeomorphic framework. The statistical methodology extracts pertinent directions in the momenta space and the clinical response space in terms of latent variables. We report the results of this analysis on 313 subjects from the Mild Cognitive Impairment group in the Alzheimer's Disease Neuroimaging Initiative (ADNI).