Zuhong Lu

Predictive models of autism spectrum disorder based on brain regional cortical thickness

Autism spectrum disorder (ASD) is a neurodevelopmental disorder with a wide phenotypic range, often affecting personality and communication. Previous voxel-based morphometry (VBM) studies of ASD have identified both gray- and white-matter volume changes. However, the cerebral cortex is a 2-D sheet with a highly folded and curved geometry, which VBM cannot directly measure. Surface-based morphometry (SBM) has the advantage of being able to measure cortical surface features, such as thickness. The goals of this study were twofold: to construct diagnostic models for ASD, based on regional thickness measurements extracted from SBM, and to compare these models to diagnostic models based on volumetric morphometry. Our study included 22 subjects with ASD (mean age 9.2+/-2.1 years) and 16 volunteer controls (mean age 10.0+/-1.9 years). Using SBM, we obtained regional cortical thicknesses for 66 brain structures for each subject. In addition, we obtained volumes for the same 66 structures for these subjects. To generate diagnostic models, we employed four machine-learning techniques: support vector machines (SVMs), multilayer perceptrons (MLPs), functional trees (FTs), and logistic model trees (LMTs). We found that thickness-based diagnostic models were superior to those based on regional volumes. For thickness-based classification, LMT achieved the best classification performance, with accuracy=87%, area under the receiver operating characteristic (ROC) curve (AUC)=0.93, sensitivity=95%, and specificity=75%. For volume-based classification, LMT achieved the highest accuracy, with accuracy=74%, AUC=0.77, sensitivity=77%, and specificity=69%. The thickness-based diagnostic model generated by LMT included 7 structures. Relative to controls, children with ASD had decreased cortical thickness in the left and right pars triangularis, left medial orbitofrontal gyrus, left parahippocampal gyrus, and left frontal pole, and increased cortical thickness in the left caudal anterior cingulate and left precuneus. Overall, thickness-based classification outperformed volume-based classification across a variety of classification methods.

White matter impairments in autism, evidence from voxel-based morphometry and diffusion tensor imaging.

This study explored white matter abnormalities in a group of Chinese children with high functioning autism (HFA). Twelve male children with HFA and ten matched typically developing children underwent diffusion tensor imaging (DTI) as well three-dimensional T1-weighted MRI for voxel-based morphometry (VBM). We found a significant decrease of the white matter density in the right frontal lobe, left parietal lobe and right anterior cingulate and a significant increase in the right frontal lobe, left parietal lobe and left cingulate gyrus in the HFA group compared with the control group. The HFA group also had decreased FA in the frontal lobe and left temporal lobe. By combining DT-MRI FA and MRI volumetric analyses based on the VBM model, the results showed consistent white matter abnormalities in a group of Chinese children with HFA.

Analyzing brain networks with PCA and conditional Granger causality.

Identifying directional influences in anatomical and functional circuits presents one of the greatest challenges for understanding neural computations in the brain. Granger causality mapping (GCM) derived from vector autoregressive models of data has been employed for this purpose, revealing complex temporal and spatial dynamics underlying cognitive processes. However, the traditional GCM methods are computationally expensive, as signals from thousands of voxels within selected regions of interest (ROIs) are individually processed, and being based on pairwise Granger causality, they lack the ability to distinguish direct from indirect connectivity among brain regions. In this work a new algorithm called PCA based conditional GCM is proposed to overcome these problems. The algorithm implements the following two procedures: (i) dimensionality reduction in ROIs of interest with principle component analysis (PCA), and (ii) estimation of the direct causal influences in local brain networks, using conditional Granger causality. Our results show that the proposed method achieves greater accuracy in detecting network connectivity than the commonly used pairwise Granger causality method. Furthermore, the use of PCA components in conjunction with conditional GCM greatly reduces the computational cost relative to the use of individual voxel time series. Hum Brain Mapp, 2009.

Voxel-based morphometry study on brain structure in children with high-functioning autism.

Earlier studies have suggested abnormal brain volumes in autism, but inconsistencies exist. Using voxel-based morphometry, we compared global and regional brain volumes in 17 high-functioning autistic children with 15 matched controls. We identified significant reduction in left white matter volume and white/gray matter ratio in autism. Regional brain volume reductions were detected for right anterior cingulate, left superior parietal lobule white matter volumes, and right parahippocampal gyrus gray matter volume, whereas enlargements in bilateral supramarginal gyrus, right postcentral gyrus, right medial frontal gyrus, and right posterior lobe of cerebellum gray matter in autism. Our findings showed global and regional brain volumes abnormality in high-functioning autism.

Detecting abnormalities of corpus callosum connectivity in autism using magnetic resonance imaging and diffusion tensor tractography

The corpus callosum (CC) has emerged as one of the primary targets of autism research. To detect aberrant CC interhemispheric connectivity in autism, we performed T1-weighted magnetic resonance imaging (MRI) and diffusion tensor imaging (DTI)-based tractography in 18 children with high functioning autism (HFA) and 16 well-matched typically developing (TD) children. We compared global and regional T1 measures (CC volume, and CC density), and the DTI measures [fractional anisotropy (FA), apparent diffusion coefficient (ADC), average fiber length (AFL), and fiber number (FN)] of transcallosal fibers, between the two groups. We also evaluated the relationships between scores on the Childhood Autism Rating Scale (CARS) and CC T1 or DTI measurements. Significantly less white matter density in the anterior third of the CC, and higher ADC and lower FN values of the anterior third transcallosal fiber tracts were found in HFA patients compared to TD children. These results suggested that the anterior third CC density and transcallosal fiber connectivity were affected in HFA children.

Gender difference in hemodynamic responses of prefrontal area to emotional stress by near-infrared spectroscopy.

Presentation of negative pictures was used as emotional stress to assess gender differences in prefrontal area activation in a functional near-infrared spectroscopy (NIRS) study. Compared with neutral condition, the response of oxy-HB for men yielded no significant difference during stress period, but the response induced by stress pictures for women showed significant enhancement. It was indicated that it is crucial to take gender difference into account when negative stimuli are used in functional brain imaging.

Detecting directional influence in fMRI connectivity analysis using PCA based Granger causality

An fMRI connectivity analysis approach combining both principal component analysis (PCA) and Granger causality method (GCM) is proposed to study directional influence between functional brain regions. Both simulated data and human fMRI data obtained during behavioral tasks were used to validate this method. If PCA is first used to reduce number of fMRI time series, then more energy and information features in the signal can be preserved than using averaged values from brain regions of interest. Subsequently, GCM can be applied to principal components extracted in order to further investigate effective connectivity. The simulation demonstrated that by using GCM with PCA, between-region causalities were better represented than using GCM with average values. Furthermore, after localizing an emotion task-induced activation in the anterior cingulate cortex, inferior frontal sulcus and amygdala, the directional influences among these brain regions were resolved using our new approach. These results indicate that using PCA may improve upon application of existing GCMs in study of human brain effective connectivity.

A conditional Granger causality model approach for group analysis in functional magnetic resonance imaging

Granger causality model (GCM) derived from multivariate vector autoregressive models of data has been employed to identify effective connectivity in the human brain with functional magnetic resonance imaging (fMRI) and to reveal complex temporal and spatial dynamics underlying a variety of cognitive processes. In the most recent fMRI effective connectivity measures, pair-wise GCM has commonly been applied based on single-voxel values or average values from special brain areas at the group level. Although a few novel conditional GCM methods have been proposed to quantify the connections between brain areas, our study is the first to propose a viable standardized approach for group analysis of fMRI data with GCM. To compare the effectiveness of our approach with traditional pair-wise GCM models, we applied a well-established conditional GCM to preselected time series of brain regions resulting from general linear model (GLM) and group spatial kernel independent component analysis of an fMRI data set in the temporal domain. Data sets consisting of one task-related and one resting-state fMRI were used to investigate connections among brain areas with the conditional GCM method. With the GLM-detected brain activation regions in the emotion-related cortex during the block design paradigm, the conditional GCM method was proposed to study the causality of the habituation between the left amygdala and pregenual cingulate cortex during emotion processing. For the resting-state data set, it is possible to calculate not only the effective connectivity between networks but also the heterogeneity within a single network. Our results have further shown a particular interacting pattern of default mode network that can be characterized as both afferent and efferent influences on the medial prefrontal cortex and posterior cingulate cortex. These results suggest that the conditional GCM approach based on a linear multivariate vector autoregressive model can achieve greater accuracy in detecting network connectivity than the widely used pair-wise GCM, and this group analysis methodology can be quite useful to extend the information obtainable in fMRI.

Investigating univariate temporal patterns for intrinsic connectivity networks based on complexity and low-frequency oscillation: A test–retest reliability study

Intrinsic connectivity networks (ICNs) are composed of spatial components and time courses. The spatial components of ICNs were discovered with moderate-to-high reliability. So far as we know, few studies focused on the reliability of the temporal patterns for ICNs based their individual time courses. The goals of this study were twofold: to investigate the test-retest reliability of temporal patterns for ICNs, and to analyze these informative univariate metrics. Additionally, a correlation analysis was performed to enhance interpretability. Our study included three datasets: (a) short- and long-term scans, (b) multi-band echo-planar imaging (mEPI), and (c) eyes open or closed. Using dual regression, we obtained the time courses of ICNs for each subject. To produce temporal patterns for ICNs, we applied two categories of univariate metrics: network-wise complexity and network-wise low-frequency oscillation. Furthermore, we validated the test-retest reliability for each metric. The network-wise temporal patterns for most ICNs (especially for default mode network, DMN) exhibited moderate-to-high reliability and reproducibility under different scan conditions. Network-wise complexity for DMN exhibited fair reliability (ICC<0.5) based on eyes-closed sessions. Specially, our results supported that mEPI could be a useful method with high reliability and reproducibility. In addition, these temporal patterns were with physiological meanings, and certain temporal patterns were correlated to the node strength of the corresponding ICN. Overall, network-wise temporal patterns of ICNs were reliable and informative and could be complementary to spatial patterns of ICNs for further study.

Predictive models for subtypes of autism spectrum disorder based on single-nucleotide polymorphisms and magnetic resonance imaging.

PURPOSE Autism spectrum disorder (ASD) is a neurodevelopmental disorder, of which Asperger syndrome and high-functioning autism are subtypes. Our goal is: 1) to determine whether a diagnostic model based on single-nucleotide polymorphisms (SNPs), brain regional thickness measurements, or brain regional volume measurements can distinguish Asperger syndrome from high-functioning autism; and 2) to compare the SNP, thickness, and volume-based diagnostic models. MATERIAL AND METHODS Our study included 18 children with ASD: 13 subjects with high-functioning autism and 5 subjects with Asperger syndrome. For each child, we obtained 25 SNPs for 8 ASD-related genes; we also computed regional cortical thicknesses and volumes for 66 brain structures, based on structural magnetic resonance (MR) examination. To generate diagnostic models, we employed five machine-learning techniques: decision stump, alternating decision trees, multi-class alternating decision trees, logistic model trees, and support vector machines. RESULTS For SNP-based classification, three decision-tree-based models performed better than the other two machine-learning models. The performance metrics for three decision-tree-based models were similar: decision stump was modestly better than the other two methods, with accuracy = 90%, sensitivity = 0.95 and specificity = 0.75. All thickness and volume-based diagnostic models performed poorly. The SNP-based diagnostic models were superior to those based on thickness and volume. For SNP-based classification, rs878960 in GABRB3 (gamma-aminobutyric acid A receptor, beta 3) was selected by all tree-based models. CONCLUSION Our analysis demonstrated that SNP-based classification was more accurate than morphometry-based classification in ASD subtype classification. Also, we found that one SNP--rs878960 in GABRB3--distinguishes Asperger syndrome from high-functioning autism.