We aimed to improve on the single-atlas ventricular segmentation method of (Carmichael, O.T., Thompson, P.M., Dutton, R.A., Lu, A., Lee, S.E., Lee, J.Y., Kuller, L.H., Lopez, O.L., Aizenstein, H.J., Meltzer, C.C., Liu, Y., Toga, A.W., Becker, J.T., 2006. Mapping ventricular changes related to dementia and mild cognitive impairment in a large community-based cohort. IEEE ISBI. 315-318) by using multi-atlas segmentation, which has been shown to lead to more accurate segmentations (Chou, Y., Leporé, N., de Zubicaray, G., Carmichael, O., Becker, J., Toga, A., Thompson, P., 2008. Automated ventricular mapping with multi-atlas fluid image alignment reveals genetic effects in Alzheimer's disease, NeuroImage 40(2): 615-630); with this method, we calculated minimal numbers of subjects needed to detect correlations between clinical scores and ventricular maps. We also assessed correlations between emerging CSF biomarkers of Alzheimer's disease pathology and localizable deficits in the brain, in 80 AD, 80 mild cognitive impairment (MCI), and 80 healthy controls from the Alzheimer's Disease Neuroimaging Initiative. Six expertly segmented images and their embedded parametric mesh surfaces were fluidly registered to each brain; segmentations were averaged within subjects to reduce errors. Surface-based statistical maps revealed powerful correlations between surface morphology and 4 variables: (1) diagnosis, (2) depression severity, (3) cognitive function at baseline, and (4) future cognitive decline over the following year. Cognitive function was assessed using the mini-mental state exam (MMSE), global and sum-of-boxes clinical dementia rating (CDR) scores, at baseline and 1-year follow-up. Lower CSF Abeta(1-42) protein levels, a biomarker of AD pathology assessed in 138 of the 240 subjects, were correlated with lateral ventricular expansion. Using false discovery rate (FDR) methods, 40 and 120 subjects, respectively, were needed to discriminate AD and MCI from normal groups. 120 subjects were required to detect correlations between ventricular enlargement and MMSE, global CDR, sum-of-boxes CDR and clinical depression scores. Ventricular expansion maps correlate with pathological and cognitive measures in AD, and may be useful in future imaging-based clinical trials.

Imaging traits provide a powerful and biologically relevant substrate to examine the influence of genetics on the brain. Interest in genome-wide, brain-wide search for influential genetic variants is growing, but has mainly focused on univariate, SNP-based association tests. Moving to gene-based multivariate statistics, we can test the combined effect of multiple genetic variants in a single test statistic. Multivariate models can reduce the number of statistical tests in gene-wide or genome-wide scans and may discover gene effects undetectable with SNP-based methods. Here we present a gene-based method for associating the joint effect of single nucleotide polymorphisms (SNPs) in 18,044 genes across 31,662 voxels of the whole brain in 731 elderly subjects (mean age: 75.56±6.82SD years; 430 males) from the Alzheimer's Disease Neuroimaging Initiative (ADNI). Structural MRI scans were analyzed using tensor-based morphometry (TBM) to compute 3D maps of regional brain volume differences compared to an average template image based on healthy elderly subjects. Using the voxel-level volume difference values as the phenotype, we selected the most significantly associated gene (out of 18,044) at each voxel across the brain. No genes identified were significant after correction for multiple comparisons, but several known candidates were re-identified, as were other genes highly relevant to brain function. GAB2, which has been previously associated with late-onset AD, was identified as the top gene in this study, suggesting the validity of the approach. This multivariate, gene-based voxelwise association study offers a novel framework to detect genetic influences on the brain.

The aim of this study was to investigate the longitudinal positron emission tomography (PET) metabolic changes in the elderly.

One key issue that must be addressed during the development of image segmentation algorithms is the accuracy of the results they produce. Algorithm developers require this so they can see where methods need to be improved and see how new developments compare with existing ones. Users of algorithms also need to understand the characteristics of algorithms when they select and apply them to their neuroimaging analysis applications. Many metrics have been proposed to characterize error and success rates in segmentation, and several datasets have also been made public for evaluation. Still, the methodologies used in analyzing and reporting these results vary from study to study, so even when studies use the same metrics their numerical results may not necessarily be directly comparable. To address this problem, we developed a web-based resource for evaluating the performance of skull-stripping in T1-weighted MRI. The resource provides both the data to be segmented and an online application that performs a validation study on the data. Users may download the test dataset, segment it using whichever method they wish to assess, and upload their segmentation results to the server. The server computes a series of metrics, displays a detailed report of the validation results, and archives these for future browsing and analysis. We applied this framework to the evaluation of 3 popular skull-stripping algorithms--the Brain Extraction Tool [Smith, S.M., 2002. Fast robust automated brain extraction. Hum. Brain Mapp. 17 (3),143-155 (Nov)], the Hybrid Watershed Algorithm [Ségonne, F., Dale, A.M., Busa, E., Glessner, M., Salat, D., Hahn, H.K., Fischl, B., 2004. A hybrid approach to the skull stripping problem in MRI. NeuroImage 22 (3), 1060-1075 (Jul)], and the Brain Surface Extractor [Shattuck, D.W., Sandor-Leahy, S.R., Schaper, K.A., Rottenberg, D.A., Leahy, R.M., 2001. Magnetic resonance image tissue classification using a partial volume model. NeuroImage 13 (5), 856-876 (May) under several different program settings. Our results show that with proper parameter selection, all 3 algorithms can achieve satisfactory skull-stripping on the test data.

Contemporary informatics and genomics research require efficient, flexible and robust management of large heterogeneous data, advanced computational tools, powerful visualization, reliable hardware infrastructure, interoperability of computational resources, and detailed data and analysis-protocol provenance. The Pipeline is a client-server distributed computational environment that facilitates the visual graphical construction, execution, monitoring, validation and dissemination of advanced data analysis protocols.

Supernumerary sex chromosome aneuploidies (sSCA) are characterized by the presence of one or more additional sex chromosomes in an individual's karyotype; they affect around 1 in 400 individuals. Although there is high variability, each sSCA subtype has a characteristic set of cognitive and physical phenotypes. Here, we investigated the differences in the morphometry of the human corpus callosum (CC) between sex-matched controls 46,XY (N =99), 46,XX (N =93), and six unique sSCA karyotypes: 47,XYY (N =29), 47,XXY (N =58), 48,XXYY (N =20), 47,XXX (N =30), 48,XXXY (N =5), and 49,XXXXY (N =6).

Correcting the effect of multiple testing is important in statistical parametric mapping. If the threshold is too liberal, then spurious claims may flood in; if it is too conservative, then true hints may be overlooked. It is highly desirable to combine random field theory and the false discovery rate (FDR) to achieve more powerful detection under gauged topological errors. However, the current FDR method based on peak height does not fully meet this expectation, and sometimes is more conservative than the traditional family-wise error rate method, for unexplained reasons. In this paper, we introduce a new topological FDR method based on signal height. As analyzed in theory and validated with extensive experiments, it controls error rates much more accurately than the peak FDR method does, and substantially gains detection power. In addition, we discover reasons behind the peak FDR method's under-performance, and formulate equations to predict the two methods' behavior.

Whole-genome and exome sequencing have already proven to be essential and powerful methods to identify genes responsible for simple Mendelian inherited disorders. These methods can be applied to complex disorders as well, and have been adopted as one of the current mainstream approaches in population genetics. These achievements have been made possible by next generation sequencing (NGS) technologies, which require substantial bioinformatics resources to analyze the dense and complex sequence data. The huge analytical burden of data from genome sequencing might be seen as a bottleneck slowing the publication of NGS papers at this time, especially in psychiatric genetics. We review the existing methods for processing NGS data, to place into context the rationale for the design of a computational resource. We describe our method, the Graphical Pipeline for Computational Genomics (GPCG), to perform the computational steps required to analyze NGS data. The GPCG implements flexible workflows for basic sequence alignment, sequence data quality control, single nucleotide polymorphism analysis, copy number variant identification, annotation, and visualization of results. These workflows cover all the analytical steps required for NGS data, from processing the raw reads to variant calling and annotation. The current version of the pipeline is freely available at http://pipeline.loni.ucla.edu. These applications of NGS analysis may gain clinical utility in the near future (e.g., identifying miRNA signatures in diseases) when the bioinformatics approach is made feasible. Taken together, the annotation tools and strategies that have been developed to retrieve information and test hypotheses about the functional role of variants present in the human genome will help to pinpoint the genetic risk factors for psychiatric disorders.

Characterization of the complex branching architecture of cerebral arteries across a representative sample of the human population is important for diagnosing, analyzing, and predicting pathological states. Brain arterial vasculature can be visualized by magnetic resonance angiography (MRA). However, most MRA studies are limited to qualitative assessments, partial morphometric analyses, individual (or small numbers of) subjects, proprietary datasets, or combinations of the above limitations. Neuroinformatics tools, developed for neuronal arbor analysis, were used to quantify vascular morphology from 3T time-of-flight MRA high-resolution (620 μm isotropic) images collected in 61 healthy volunteers (36/25 F/M, average age=31.2 ± 10.7, range=19-64 years). We present in-depth morphometric analyses of the global and local anatomical features of these arbors. The overall structure and size of the vasculature did not significantly differ across genders, ages, or hemispheres. The total length of the three major arterial trees stemming from the circle of Willis (from smallest to largest: the posterior, anterior, and middle cerebral arteries; or PCAs, ACAs, and MCAs, respectively) followed an approximate 1:2:4 proportion. Arterial size co-varied across individuals: subjects with one artery longer than average tended to have all other arteries also longer than average. There was no net right-left difference across the population in any of the individual arteries, but ACAs were more lateralized than MCAs. MCAs, ACAs, and PCAs had similar branch-level properties such as bifurcation angles. Throughout the arterial vasculature, there were considerable differences between branch types: bifurcating branches were significantly shorter and straighter than terminating branches. Furthermore, the length and meandering of bifurcating branches increased with age and with path distance from the circle of Willis. All reconstructions are freely distributed through a public database to enable additional analyses and modeling (cng.gmu.edu/brava).

The NTRK3 gene (also known as TRKC) encodes a high affinity receptor for the neurotrophin 3'-nucleotidase (NT3), which is implicated in oligodendrocyte and myelin development. We previously found that white matter integrity in young adults is related to common variants in genes encoding neurotrophins and their receptors. This underscores the importance of neurotrophins for white matter development. NTRK3 variants are putative risk factors for schizophrenia, bipolar disorder, and obsessive-compulsive disorder hoarding, suggesting that some NTRK3 variants may affect the brain. To test this, we scanned 392 healthy adult twins and their siblings (mean age, 23.6 ± 2.2 years; range: 20-29 years) with 105-gradient 4-Tesla diffusion tensor imaging (DTI). We identified 18 single nucleotide polymorphisms (SNPs) in the NTRK3 gene that have been associated with neuropsychiatric disorders. We used a multi-SNP model, adjusting for family relatedness, age, and sex, to relate these variants to voxelwise fractional anisotropy (FA) - a DTI measure of white matter integrity. FA was optimally predicted (based on the highest false discovery rate critical p), by five SNPs (rs1017412, rs2114252, rs16941261, rs3784406, and rs7176429; overall FDR critical p=0.028). Gene effects were widespread and included the corpus callosum genu and inferior longitudinal fasciculus - regions implicated in several neuropsychiatric disorders and previously associated with other neurotrophin-related genetic variants in an overlapping sample of subjects. NTRK3 genetic variants, and neurotrophins more generally, may influence white matter integrity in brain regions implicated in neuropsychiatric disorders.

The insula, hidden deep within the Sylvian fissures, has proven difficult to study from a connectivity perspective. Most of our current information on the anatomical connectivity of the insula comes from studies of nonhuman primates and post mortem human dissections. To date, only two neuroimaging studies have successfully examined the connectivity of the insula. Here we examine how the connectivity of the insula develops between ages 12 and 30, in 307 young adolescent and adult subjects scanned with 4-Tesla high angular resolution diffusion imaging (HARDI). The density of fiber connections between the insula and the frontal and parietal cortex decreased with age, but the connection density between the insula and the temporal cortex generally increased with age. This trajectory is in line with well-known patterns of cortical development in these regions. In addition, males and females showed different developmental trajectories for the connection between the left insula and the left precentral gyrus. The insula plays many different roles, some of them affected in neuropsychiatric disorders; this information on the insula's connectivity may help efforts to elucidate mechanisms of brain disorders in which it is implicated.

Available approaches to the investigation of traumatic brain injury (TBI) are frequently hampered, to some extent, by the unsatisfactory abilities of existing methodologies to efficiently define and represent affected structural connectivity and functional mechanisms underlying TBI-related pathology. In this paper, we describe a patient-tailored framework which allows mapping and characterization of TBI-related structural damage to the brain via multimodal neuroimaging and personalized connectomics. Specifically, we introduce a graphically driven approach for the assessment of trauma-related atrophy of white matter connections between cortical structures, with relevance to the quantification of TBI chronic case evolution. This approach allows one to inform the formulation of graphical neurophysiological and neuropsychological TBI profiles based on the particular structural deficits of the affected patient. In addition, it allows one to relate the findings supplied by our workflow to the existing body of research that focuses on the functional roles of the cortical structures being targeted. A graphical means for representing patient TBI status is relevant to the emerging field of personalized medicine and to the investigation of neural atrophy.

Diffusion weighted magnetic resonance imaging (DW-MRI) are now widely used to assess brain integrity in clinical populations. The growing interest in mapping brain connectivity has made it vital to consider what scanning parameters affect the accuracy, stability, and signal-to-noise of diffusion measures. Trade-offs between scan parameters can only be optimized if their effects on various commonly-derived measures are better understood. To explore angular versus spatial resolution trade-offs in standard tensor-derived measures, and in measures that use the full angular information in diffusion signal, we scanned eight subjects twice, 2 weeks apart, using three protocols that took the same amount of time (7 min). Scans with 3.0, 2.7, 2.5 mm isotropic voxels were collected using 48, 41, and 37 diffusion-sensitized gradients to equalize scan times. A specially designed DTI phantom was also scanned with the same protocols, and different b-values. We assessed how several diffusion measures including fractional anisotropy (FA), mean diffusivity (MD), and the full 3D orientation distribution function (ODF) depended on the spatial/angular resolution and the SNR. We also created maps of stability over time in the FA, MD, ODF, skeleton FA of 14 TBSS-derived ROIs, and an information uncertainty index derived from the tensor distribution function, which models the signal using a continuous mixture of tensors. In scans of the same duration, higher angular resolution and larger voxels boosted SNR and improved stability over time. The increased partial voluming in large voxels also led to bias in estimating FA, but this was partially addressed by using "beyond-tensor" models of diffusion.

Several cortical regions are reported to vary in meditation practitioners. However, prior analyses have focused primarily on examining gray matter or cortical thickness. Thus, additional effects with respect to other cortical features might have remained undetected. Gyrification (the pattern and degree of cortical folding) is an important cerebral characteristic related to the geometry of the brain's surface. Thus, exploring cortical gyrification in long-term meditators may provide additional clues with respect to the underlying anatomical correlates of meditation. This study examined cortical gyrification in a large sample (n = 100) of meditators and controls, carefully matched for sex and age. Cortical gyrification was established by calculating mean curvature across thousands of vertices on individual cortical surface models. Pronounced group differences indicating larger gyrification in meditators were evident within the left precentral gyrus, right fusiform gyrus, right cuneus, as well as left and right anterior dorsal insula (the latter representing the global significance maximum). Positive correlations between gyrification and the number of meditation years were similarly pronounced in the right anterior dorsal insula. Although the exact functional implications of larger cortical gyrification remain to be established, these findings suggest the insula to be a key structure involved in aspects of meditation. For example, variations in insular complexity could affect the regulation of well-known distractions in the process of meditation, such as daydreaming, mind-wandering, and projections into past or future. Moreover, given that meditators are masters in introspection, awareness, and emotional control, increased insular gyrification may reflect an integration of autonomic, affective, and cognitive processes. Due to the cross-sectional nature of this study, further research is necessary to determine the relative contribution of nature and nurture to links between cortical gyrification and meditation.

Autism spectrum disorder is a heterogeneous disorder of brain development with wide ranging cognitive deficits. Typically diagnosed before age 3, autism spectrum disorder is behaviorally defined but patients are thought to have protracted alterations in brain maturation. With longitudinal magnetic resonance imaging (MRI), we mapped an anomalous developmental trajectory of the brains of autistic compared with those of typically developing children and adolescents. Using tensor-based morphometry, we created 3D maps visualizing regional tissue growth rates based on longitudinal brain MRI scans of 13 autistic and seven typically developing boys (mean age/interscan interval: autism 12.0 ± 2.3 years/2.9 ± 0.9 years; control 12.3 ± 2.4/2.8 ± 0.8). The typically developing boys demonstrated strong whole brain white matter growth during this period, but the autistic boys showed abnormally slowed white matter development (P = 0.03, corrected), especially in the parietal (P = 0.008), temporal (P = 0.03), and occipital lobes (P = 0.02). We also visualized abnormal overgrowth in autism in gray matter structures such as the putamen and anterior cingulate cortex. Our findings reveal aberrant growth rates in brain regions implicated in social impairment, communication deficits and repetitive behaviors in autism, suggesting that growth rate abnormalities persist into adolescence. Tensor-based morphometry revealed persisting growth rate anomalies long after diagnosis, which has implications for evaluation of therapeutic effects.

Large-archives of neuroimaging data present many opportunities for re-analysis and mining that can lead to new findings of use in basic research or in the characterization of clinical syndromes. However, interaction with such archives tends to be driven textually, based on subject or image volume meta-data, not the actual neuroanatomical morphology itself, for which the imaging was performed to measure. What is needed is a content-driven approach for examining not only the image content itself but to explore brains that are anatomically similar, and identifying patterns embedded within entire sets of neuroimaging data. With the aim of visual navigation of large- scale neurodatabases, we introduce the concept of brain meta-spaces. The meta-space encodes pair-wise dissimilarities between all individuals in a population and shows the relationships between brains as a navigable framework for exploration. We employ multidimensional scaling (MDS) to implement meta-space processing for a new coordinate system that distributes all data points (brain surfaces) in a common frame-of-reference, with anatomically similar brain data located near each other. To navigate within this derived meta-space, we have developed a fully interactive 3D visualization environment that allows users to examine hundreds of brains simultaneously, visualize clusters of brains with similar characteristics, zoom in on particular instances, and examine the surface topology of an individual brain's surface in detail. The visualization environment not only displays the dissimilarities between brains, but also renders complete surface representations of individual brain structures, allowing an instant 3D view of the anatomies, as well as their differences. The data processing is implemented in a grid-based setting using the LONI Pipeline workflow environment. Additionally users can specify a range of baseline brain atlas spaces as the underlying scale for comparative analyses. The novelty in our approach lies in the user ability to simultaneously view and interact with many brains at once but doing so in a vast meta-space that encodes (dis) similarity in morphometry. We believe that the concept of brain meta-spaces has important implications for the future of how users interact with large-scale archives of primary neuroimaging data.

Measures of brain change can be computed from sequential MRI scans, providing valuable information on disease progression, e.g., for patient monitoring and drug trials. Tensor-based morphometry (TBM) creates maps of these brain changes, visualizing the 3D profile and rates of tissue growth or atrophy, but its sensitivity depends on the contrast and geometric stability of the images. As part of the Alzheimer's Disease Neuroimaging Initiative (ADNI), 17 normal elderly subjects were scanned twice (at a 2-week interval) with several 3D 1.5 T MRI pulse sequences: high and low flip angle SPGR/FLASH (from which Synthetic T1 images were generated), MP-RAGE, IR-SPGR (N = 10) and MEDIC (N = 7) scans. For each subject and scan type, a 3D deformation map aligned baseline and follow-up scans, computed with a nonlinear, inverse-consistent elastic registration algorithm. Voxelwise statistics, in ICBM stereotaxic space, visualized the profile of mean absolute change and its cross-subject variance; these maps were then compared using permutation testing. Image stability depended on: (1) the pulse sequence; (2) the transmit/receive coil type (birdcage versus phased array); (3) spatial distortion corrections (using MEDIC sequence information); (4) B1-field intensity inhomogeneity correction (using N3). SPGR/FLASH images acquired using a birdcage coil had least overall deviation. N3 correction reduced coil type and pulse sequence differences and improved scan reproducibility, except for Synthetic T1 images (which were intrinsically corrected for B1-inhomogeneity). No strong evidence favored B0 correction. Although SPGR/FLASH images showed least deviation here, pulse sequence selection for the ADNI project was based on multiple additional image analyses, to be reported elsewhere.

The Parkinson's Progression Markers Initiative (PPMI) is an observational, international study designed to establish biomarker-defined cohorts and identify clinical, imaging, genetic, and biospecimen Parkinson's disease (PD) progression markers to accelerate disease-modifying therapeutic trials.

The objective of this study was to assess longitudinal change in clinical and dopamine transporter imaging outcomes in early, untreated PD.

A long-standing question is how to best use brain morphometric and genetic data to distinguish Alzheimer's disease (AD) patients from cognitively normal (CN) subjects and to predict those who will progress from mild cognitive impairment (MCI) to AD. Here, we use a neural network (NN) framework on both magnetic resonance imaging-derived quantitative structural brain measures and genetic data to address this question. We tested the effectiveness of NN models in classifying and predicting AD. We further performed a novel analysis of the NN model to gain insight into the most predictive imaging and genetics features and to identify possible interactions between features that affect AD risk. Data were obtained from the AD Neuroimaging Initiative cohort and included baseline structural MRI data and single nucleotide polymorphism (SNP) data for 138 AD patients, 225 CN subjects, and 358 MCI patients. We found that NN models with both brain and SNP features as predictors perform significantly better than models with either alone in classifying AD and CN subjects, with an area under the receiver operating characteristic curve (AUC) of 0.992, and in predicting the progression from MCI to AD (AUC=0.835). The most important predictors in the NN model were the left middle temporal gyrus volume, the left hippocampus volume, the right entorhinal cortex volume, and the APOE (a gene that encodes apolipoprotein E) ɛ4 risk allele. Furthermore, we identified interactions between the right parahippocampal gyrus and the right lateral occipital gyrus, the right banks of the superior temporal sulcus and the left posterior cingulate, and SNP rs10838725 and the left lateral occipital gyrus. Our work shows the ability of NN models to not only classify and predict AD occurrence but also to identify important AD risk factors and interactions among them.