The Center for Computational Biology (CCB) is a multidisciplinary program where biomedical scientists, engineers, and clinicians work jointly to combine modern mathematical and computational techniques, to perform phenotypic and genotypic studies of biological structure, function, and physiology in health and disease. CCB has developed a computational framework built around the Manifold Atlas, an integrated biomedical computing environment that enables statistical inference on biological manifolds. These manifolds model biological structures, features, shapes, and flows, and support sophisticated morphometric and statistical analyses. The Manifold Atlas includes tools, workflows, and services for multimodal population-based modeling and analysis of biological manifolds. The broad spectrum of biomedical topics explored by CCB investigators include the study of normal and pathological brain development, maturation and aging, discovery of associations between neuroimaging and genetic biomarkers, and the modeling, analysis, and visualization of biological shape, form, and size. CCB supports a wide range of short-term and long-term collaborations with outside investigators, which drive the center's computational developments and focus the validation and dissemination of CCB resources to new areas and scientific domains.

The link between brain structure and intelligence is a well-investigated topic, but existing analyses have mainly focused on adult samples. Studies in healthy children and adolescents are rare, and normative data specifically addressing the association between corpus callosum morphology and intellectual abilities is quite limited. To advance this field of research, we mapped the correlations between standardized intelligence measures and callosal thickness based on high-resolution magnetic resonance imaging (MRI) data. Our large and well-matched sample included 200 normally developing subjects (100 males, 100 females) ranging from 6 to 17 years of age. Although the strongest correlations were negative and confined to the splenium, the strength and the direction of intelligence-callosal thickness associations varied considerably with respect to age and sex. While significant correlations in females were mainly positive, significant correlations in males were exclusively negative. However, only the negative correlations in the overall sample (i.e., males and females combined) remained significant when controlling for multiple comparisons. The observed negative correlations between callosal thickness and intelligence in children and adolescents contrast with the positive correlations typically reported in adult samples. However, negative correlations are in line with reports from other pediatric studies relating cognitive measures to other brain attributes such as cortical thickness, gray matter volume, and gray matter density. Altogether, these findings suggest that relationships between callosal morphology and cognition are highly dynamic during brain maturation. Sex effects on links between callosal thickness and intelligence during childhood and adolescence are present but appear rather weak in general.

The white matter of the brain consists of fiber tracts that connect different regions of the brain. Among these tracts, the intrahemispheric cortico-cortical connections are called association fibers. The U-fibers are short association fibers that connect adjacent gyri. These fibers were thought to work as part of the cortico-cortical networks to execute associative brain functions. However, their anatomy and functions have not been documented in detail for the human brain. In past studies, U-fibers have been characterized in the human brain with diffusion tensor imaging (DTI). However, the validity of such findings remains unclear. In this study, DTI of the macaque brain was performed, and the anatomy of U-fibers was compared with that of the human brain reported in a previous study. The macaque brain was chosen because it is the most commonly used animal model for exploring cognitive functions and the U-fibers of the macaque brain have been already identified by axonal tracing studies, which makes it an ideal system for confirming the DTI findings. Ten U-fibers found in the macaque brain were also identified in the human brain, with a similar organization and topology. The delineation of these species-conserved white matter structures may provide new options for understanding brain anatomy and function.

Brain structural abnormalities and neurocognitive dysfunction have been observed in individuals with fetal alcohol spectrum disorders (FASDs). Little is known about how white matter integrity is related to these functional and morphological deficits. We used a combination of diffusion tensor and T1-weighted magnetic resonance imaging to evaluate white matter integrity in individuals with FASDs and related these findings to neurocognitive deficits. Seventeen children and adolescents with FASDs were compared with 19 typically developing age-and gender-matched controls. Lower fractional anisotropy (FA) was observed in individuals with FASDs relative to controls in the right lateral temporal lobe and bilaterally in the lateral aspects of the splenium of the corpus callosum. White matter density was also lower in some, but not all regions in which FA was lower. FA abnormalities were confirmed to be in areas of white matter in post hoc region of interest analyses, further supporting that less myelin or disorganized fiber tracts are associated with heavy prenatal alcohol exposure. Significant correlations between performance on a test of visuomotor integration and FA in bilateral splenium, but not temporal regions were observed within the FASD group. Correlations between the visuomotor task and FA within the splenium were not significant with in the control group, and were not significant for measures of reading ability. This suggests that this region of white matter is particularly susceptible to damage from prenatal alcohol exposure and that disruption of splenial fibers in this group is associated with poorer visuomotor integration.

Rapidly evolving neuroimaging techniques are producing unprecedented quantities of digital data at the same time that many research studies are evolving into global, multi-disciplinary collaborations between geographically distributed scientists. While networked computers have made it almost trivial to transmit data across long distances, collecting and analyzing this data requires extensive metadata if the data is to be maximally shared. Though it is typically straightforward to encode text and numerical values into files and send content between different locations, it is often difficult to attach context and implicit assumptions to the content. As the number of and geographic separation between data contributors grows to national and global scales, the heterogeneity of the collected metadata increases and conformance to a single standardization becomes implausible. Neuroimaging data repositories must then not only accumulate data but must also consolidate disparate metadata into an integrated view. In this article, using specific examples from our experiences, we demonstrate how standardization alone cannot achieve full integration of neuroimaging data from multiple heterogeneous sources and why a fundamental change in the architecture of neuroimaging data repositories is needed instead.

The corpus callosum changes structurally throughout life, but most dramatically during childhood and adolescence. Even so, existing studies of callosal development tend to use parcellation schemes that may not capture the complex spatial profile of anatomical changes. Thus, more detailed mapping of callosal growth processes is desirable to create a normative reference. This will help to relate and interpret other structural, functional, and behavioral measurements, both from healthy subjects and pediatric patients. We applied computational surface-based mesh-modeling methods to analyze callosal morphology at extremely high spatial resolution. We mapped callosal development and explored sex differences in a large and well-matched sample of healthy children and adolescents (n=190) aged 5 to 18 years. Except for the rostrum in females, callosal thickness increased across the whole surface, with sex- and region-specific rates of growth, and also shrinkage at times. The temporally distinct changes in callosal thickness are likely to be a consequence of varying degrees of axonal myelination, redirection, and pruning. Alternating phases of callosal growth and shrinkage may reflect a permanent adjustment and fine-tuning of fibers connecting homologous cortical areas during childhood and adolescence. Our findings emphasize the importance of taking into account sex differences in future studies, as existing developmental effects might remain disguised (or biased towards the effect of the dominant sex in unbalanced statistical designs) when pooling male and female samples.

Summary

Naturally occurring mutants and genetically manipulated strains of mice are widely used to model a variety of human diseases. Atlases are an invaluable aid in understanding the impact of such manipulations by providing a standard for comparison and to facilitate the integration of anatomic, genetic, and physiologic observations from multiple subjects and experiments. We have developed digital atlases of the C57BL/6J mouse brain (adult and neonate) as comprehensive frameworks for storing and accessing the myriad types of information about the mouse brain. Along with raw and annotated images, these contain database management systems and a set of tools for comparing information from different techniques and different animals. Each atlas establishes a canonical representation of the mouse brain and provides the tools for the manipulation and analysis of new data. We describe both these atlases and discuss how they may be put to use in organizing and analyzing data from mouse models of epilepsy.

Although the systematic study of meditation is still in its infancy, research has provided evidence for meditation-induced improvements in psychological and physiological well-being. Moreover, meditation practice has been shown not only to benefit higher-order cognitive functions but also to alter brain activity. Nevertheless, little is known about possible links to brain structure. Using high-resolution MRI data of 44 subjects, we set out to examine the underlying anatomical correlates of long-term meditation with different regional specificity (i.e., global, regional, and local). For this purpose, we applied voxel-based morphometry in association with a recently validated automated parcellation approach. We detected significantly larger gray matter volumes in meditators in the right orbito-frontal cortex (as well as in the right thalamus and left inferior temporal gyrus when co-varying for age and/or lowering applied statistical thresholds). In addition, meditators showed significantly larger volumes of the right hippocampus. Both orbito-frontal and hippocampal regions have been implicated in emotional regulation and response control. Thus, larger volumes in these regions might account for meditators’ singular abilities and habits to cultivate positive emotions, retain emotional stability, and engage in mindful behavior. We further suggest that these regional alterations in brain structures constitute part of the underlying neurological correlate of long-term meditation independent of a specific style and practice. Future longitudinal analyses are necessary to establish the presence and direction of a causal link between meditation practice and brain anatomy.

Alzheimer's disease (AD) is the most common type of dementia worldwide. Hippocampal atrophy and ventricular enlargement have been associated with AD but also with normal aging. We analyzed 1.5T brain MRI data from 46 cognitively normal elderly (NC), 33 mild cognitive impairment (MCI) and 43 AD subjects. Hippocampal and ventricular analyses were conducted with two novel semi-automated segmentation approaches followed by the radial distance mapping technique. Multiple linear regression was used to assess effects of age and diagnosis on hippocampal and ventricular volumes and radial distance. Additional 3D map correction for multiple comparisons was conducted with permutation testing. As expected, most significant hippocampal atrophy and ventricular enlargement were seen in the AD vs. NC comparison. MCI subjects showed intermediate levels of hippocampal atrophy and ventricular enlargement. Significant effects of age on hippocampal volume and radial distance were seen in the pooled sample as well as in the NC and AD groups considered separately. Age-associated differences were detected in all hippocampal subfields and the frontal and body/occipital horn portions of the lateral ventricles. Aging affects both the hippocampus and lateral ventricles independent of AD pathology and should be included as covariate in all structural hippocampal and ventricular analyses when possible.

Diffusion tensor imaging (DTI) is sensitive to the directionally- constrained flow of water, which diffuses preferentially along axons. Tractography programs may be used to infer matrices of connectivity (anatomical networks) between pairs of brain regions. Little is known about how these computed connectivity measures depend on the scans’ spatial and angular resolutions. To determine this, we scanned 8 young adults with DTI at 2.5 and 3 mm resolutions, and an additional subject at 4 resolutions between 2–4 mm. We computed 70×70 connectivity matrices, using whole-brain tractography to measure fiber density between all pairs of 70 cortical and subcortical regions. Spatial and angular resolution affected the computed connectivity for narrower tracts (internal capsule and cerebellum), but also for the corticospinal tract. Data resolution affected the apparent role of some key structures in cortical anatomic networks. Care is needed when comparing network data across studies, and interpreting apparent disagreements among findings.

Penalized or sparse regression methods are gaining increasing attention in imaging genomics, as they can select optimal regressors from a large set of predictors whose individual effects are small or mostly zero. We applied a multivariate approach, based on L1-L2-regularized regression (elastic net) to predict a magnetic resonance imaging (MRI) tensor-based morphometry-derived measure of temporal lobe volume from a genome-wide scan in 740 Alzheimer’s Disease Neuroimaging Initiative (ADNI) subjects. We tuned the elastic net model’s parameters using internal crossvalidation and evaluated the model on independent test sets. Compared to 100,000 permutations performed with randomized imaging measures, the predictions were found to be statistically significant (p ~ 0.001). The rs9933137 variant in the RBFOX1 gene was a highly contributory genotype, along with rs10845840 in GRIN2B and rs2456930, discovered previously in a univariate genomewide search.

Alzheimer’s Disease (AD) has long been considered a cortical degenerative disease, but impaired brain connectivity, due to white matter injury, may exacerbate cognitive problems. Predicting brain changes is critically important for early treatment. In a longitudinal diffusion tensor imaging study, we investigated white matter fiber integrity in 19 patients (mean age: 74.7 +/− 8.4 yrs at baseline) displaying early signs of mild cognitive impairment (eMCI). We first examined whether baseline average fractional anisotropy (FA) measures in the corpus callosum (CC) predicted changes in white matter integrity over the following 6 months. We then examined whether “small world” architecture measures - calculated from baseline connectivity maps - predicted white matter changes over the next 6 months. While average CC FA measures at baseline were not associated with future changes in FA, network measures were a sensitive biomarker for predicting white matter changes during this critical time before AD strikes.

Large multi-site image-analysis studies have successfully discovered genetic variants that affect brain structure in tens of thousands of subjects scanned worldwide. Candidate genes have also associated with brain integrity, measured using fractional anisotropy in diffusion tensor images (DTI). To evaluate the heritability and robustness of DTI measures as a target for genetic analysis, we compared 417 twins and siblings scanned on the same day on the same high field scanner (4-Tesla) with two protocols: (1) 94-directions; 2mm-thick slices, (2) 27-directions; 5mm-thickness. Using mean FA in white matter ROIs and FA ‘skeletons’ derived using FSL, we (1) examined differences in voxelwise means, variances, and correlations among the measures; and (2) assessed heritability with structural equation models, using the classical twin design. FA measures from the genu of the corpus callosum were highly heritable, regardless of protocol. Genome-wide analysis of the genu mean FA revealed differences across protocols in the top associations.

Graph theory can be applied to matrices that represent the brain’s anatomical connections, to better understand global properties of anatomical networks, such as their clustering, efficiency and “small-world” topology. Network analysis is popular in adult studies of connectivity, but only one study – in just 30 subjects – has examined how network measures change as the brain develops over this period. Here we assessed the developmental trajectory of graph theory metrics of structural brain connectivity in a cross-sectional study of 467 subjects, aged 12 to 30. We computed network measures from 70×70 connectivity matrices of fiber density generated using whole-brain tractography in 4-Tesla 105-gradient high angular resolution diffusion images (HARDI). We assessed global efficiency and modularity, and both age and age2 effects were identified. HARDI-based connectivity maps are sensitive to the remodeling and refinement of structural brain connections as the human brain develops.

Human brain connectivity is disrupted in a wide range of disorders – from Alzheimer’s disease to autism – but little is known about which specific genes affect it. Here we conducted a genome-wide association for connectivity matrices that capture information on the density of fiber connections between 70 brain regions. We scanned a large twin cohort (N=366) with 4-Tesla high angular resolution diffusion imaging (105-gradient HARDI). Using whole brain HARDI tractography, we extracted a relatively sparse 70×70 matrix representing fiber density between all pairs of cortical regions automatically labeled in co-registered anatomical scans. Additive genetic factors accounted for 1–58% of the variance in connectivity between 90 (of 122) tested nodes. We discovered genome-wide significant associations between variants and connectivity. GWAS permutations at various levels of heritability, and split-sample replication, validated our genetic findings. The resulting genes may offer new leads for mechanisms influencing aberrant connectivity and neurodegeneration.

A major challenge in neuroscience is finding which genes affect brain integrity, connectivity, and intellectual function. Discovering influential genes holds vast promise for neuroscience, but typical genome-wide searches assess around one million genetic variants one-by-one, leading to intractable false positive rates, even with vast samples of subjects. Even more intractable is the question of which genes interact and how they work together to affect brain connectivity. Here we report a novel approach that discovers which genes contribute to brain wiring and fiber integrity at all pairs of points in a brain scan. We studied genetic correlations between thousands of points in human brain images from 472 twins and their non-twin siblings (mean age: 23.7±2.1 SD years; 193 M/279 F). We combined clustering with genome-wide scanning to find brain systems with common genetic determination. We then filtered the image in a new way to boost power to find causal genes. Using network analysis, we found a network of genes that affect brain wiring in healthy young adults. Our new strategy makes it more computationally tractable to discover genes that affect brain integrity. The gene network showed small-world and scale-free topologies, suggesting efficiency in genetic interactions, and resilience to network disruption. Genetic variants at hubs of the network influence intellectual performance by modulating associations between performance intelligence quotient (IQ) and the integrity of major white matter tracts, such as the callosal genu and splenium, cingulum, optic radiations, and the superior longitudinal fasciculus.

The different brain anatomy of men and women is both a classic and continuing topic of major interest. Among the most replicated and robust sex differences are larger overall brain dimensions in men, and relative increases of global and regional gray matter (GM) in women. However, the question remains whether sex-typical differences in brain size (i.e., larger male and smaller female brains) or biological sex itself account for the observed sex effects on tissue amount and distribution. Exploring cerebral structures in men and women with similar brain size may clarify the true contribution of biological sex. We thus examined a sample of 24 male and 24 female subjects with brains identical in size, in addition to 24 male and 24 female subjects with considerable brain size differences. Using this large set of brains (n=96), we applied a well-validated and automated voxel-based approach to examine regional volumes of GM. While we revealed significant main effects of sex, there were no significant effects of brain size (and no significant interactions between sex and brain size). When conducting post hoc tests, we revealed a number of regions where women had larger GM volumes compared to men. Importantly, these sex effects remained evident when comparing men and women with the same brain size. Altogether, our findings suggest that the observed increased regional GM volumes in female brains constitute sex-dependent redistributions of tissue volume, rather than individual adjustments attributable to brain size.

The arcuate fasciculus (AF) connects cortical regions important in language processing, but how fiber coherence and organization relates to gray matter macrostructure remains uncharacterized. We used high-resolution structural and 30-direction diffusion imaging data from 36 healthy adults (24 male/12 female; mean age: 30.5±9.8 years) to establish the relationships between AF microstructure and regional variations in cortical gray matter within language networks. Cortical pattern matching algorithms were used to measure gray matter thickness at high spatial density, and a validated diffusion tractography method was used to reconstruct the AF in the left and right hemisphere of each subject. Relationships between imaging measures and neuropsychological scores of verbal fluency were additionally assessed. Results revealed positive and highly topographical associations between arcuate fractional anisotropy and cortical thickness within anterior and posterior language regions and surrounding cortices, more prominently in the left hemisphere. These regional cortical thickness/fractional anisotropy relationships were primarily attributable to variations in radial diffusivity. Associations between cortical thickness and verbal fluency were observed in perisylvian language-related regions. Language scores associated with left hemisphere AF axial diffusivity, but not with AF fractional anisotropy or radial diffusivity. These findings thus suggest that particular components of white matter microstructure and regional increases in cortical thickness benefit aspects of language processing. Further, the topographical relationships between independent measures of white matter and gray matter integrity suggest that rich developmental or environmental interactions influence brain structure and function where the presence and strength of such associations may elucidate pathophysiological processes influencing language systems.

Abstract

Although neuroimaging is essential for prompt and proper management of traumatic brain injury (TBI), there is a regrettable and acute lack of robust methods for the visualization and assessment of TBI pathophysiology, especially for of the purpose of improving clinical outcome metrics. Until now, the application of automatic segmentation algorithms to TBI in a clinical setting has remained an elusive goal because existing methods have, for the most part, been insufficiently robust to faithfully capture TBI-related changes in brain anatomy. This article introduces and illustrates the combined use of multimodal TBI segmentation and time point comparison using 3D Slicer, a widely-used software environment whose TBI data processing solutions are openly available. For three representative TBI cases, semi-automatic tissue classification and 3D model generation are performed to perform intra-patient time point comparison of TBI using multimodal volumetrics and clinical atrophy measures. Identification and quantitative assessment of extra- and intra-cortical bleeding, lesions, edema, and diffuse axonal injury are demonstrated. The proposed tools allow cross-correlation of multimodal metrics from structural imaging (e.g., structural volume, atrophy measurements) with clinical outcome variables and other potential factors predictive of recovery. In addition, the workflows described are suitable for TBI clinical practice and patient monitoring, particularly for assessing damage extent and for the measurement of neuroanatomical change over time. With knowledge of general location, extent, and degree of change, such metrics can be associated with clinical measures and subsequently used to suggest viable treatment options.

The NTRK1 gene (also known as TRKA) encodes a high affinity receptor for NGF, a neurotrophin involved in nervous system development and myelination. NTRK1 has been implicated in neurological function via links between the T allele at rs6336 (NTRK1-T) and schizophrenia risk. A variant in the neurotrophin gene, BDNF, was previously associated with white matter integrity in young adults, highlighting the importance of neurotrophins to white matter development. We hypothesized that NTRK1-T would relate to lower FA in healthy adults.

We scanned 391 healthy adult human twins and their siblings (mean age: 23.6 ± 2.2 years; 31 NTRK1-T carriers, 360 non-carriers) using 105-gradient diffusion tensor imaging at 4 Tesla. We evaluated in brain white matter how NTRK1-T and NTRK1 rs4661063 allele A (rs4661063-A, which is in moderate linkage disequilibrium with rs6336) related to voxelwise fractional anisotropy – a common diffusion tensor imaging measure of white matter microstructure. We used mixed-model regression to control for family relatedness, age, and sex. The sample was split in half to test results reproducibility. The false discovery rate method corrected for voxelwise multiple comparisons.

NTRK1-T and rs4661063-A correlated with lower white matter fractional anisotropy, independent of age and sex (multiple comparisons corrected: false discovery rate critical p = 0.038 for NTRK1-T and 0.013 for rs4661063-A). In each half-sample, the NTRK1-T effect was replicated in the cingulum, corpus callosum, superior and inferior longitudinal fasciculi, inferior fronto-occipital fasciculus, superior corona radiata, and uncinate fasciculus. Our results suggest that NTRK1-T is important for developing white matter microstructure.

Structural deficits in the frontotemporal network have been shown in individuals with psychopathy and are posited to contribute to neuropsychological impairments such as response perseveration. However, no study to date has examined structural correlates of response perseveration in individuals with psychopathy. In this structural MRI study, the authors found higher correlations between increased response perseveration and reduced cortical thickness in the orbitofrontal and anterior temporal regions in individuals with psychopathy than in healthy-comparison subjects. The findings provide preliminary evidence suggesting potential contributions of frontotemporal structural deficits in neurocognitive impairment with perseveration in individuals with psychopathy.

With the advancement of image acquisition and analysis methods in recent decades, unique opportunities have emerged to study the neuroanatomical correlates of intelligence. Traditional approaches examining global measures have been complemented by insights from more regional analyses based on pre-defined areas. Newer state-of-the-art approaches have further enhanced our ability to localize the presence of correlations between cerebral characteristics and intelligence with high anatomic precision. These in vivo assessments have confirmed mainly positive correlations, suggesting that optimally increased brain regions are associated with better cognitive performance. Findings further suggest that the models proposed to explain the anatomical substrates of intelligence should address contributions from not only (pre)frontal regions, but also widely distributed networks throughout the whole brain.

Normal human aging is accompanied by progressive brain tissue loss and cognitive decline; however, several factors are thought to influence brain aging. We applied tensor-based morphometry to high-resolution brain MRI scans to determine whether educational level or physical activity was associated with brain tissue volumes in the elderly, particularly in regions susceptible to age-related atrophy. We mapped the 3D profile of brain volume differences in 226 healthy elderly subjects (130F/96M; 77.9 ± 3.6 SD years) from the Cardiovascular Health Study-Cognition Study. Statistical maps revealed the 3D profile of brain regions whose volumes were associated with educational level and physical activity (based on leisure-time energy expenditure). After controlling for age, sex, and physical activity, higher educational levels were associated with ~2–3% greater tissue volumes, on average, in the temporal lobe gray matter. After controlling for age, sex, and education, greater physical activity was associated with ~2–2.5% greater average tissue volumes in the white matter of the corona radiata extending into the parietal-occipital junction. Body mass index (BMI) was highly correlated with both education and physical activity, so we examined BMI as a contributing factor by including physical activity, education, and BMI in the same model; only BMI effects remained significant. This is one of the largest MRI studies of factors influencing structural brain aging, and BMI may be a key factor explaining the observed relationship between education, physical activity, and brain structure. Independent contributions to brain structure could not be teased apart as all these factors were highly correlated with one another.

Recently, carriers of a common variant in the autism risk gene, CNTNAP2, were found to have altered functional brain connectivity using functional MRI. Here we scanned 328 young adults with high-field (4-Tesla) diffusion imaging, to test the hypothesis that carriers of this gene variant would have altered structural brain connectivity. All participants (209 females, 119 males, age: 23.4 +/−2.17 SD years) were scanned with 105-gradient high angular diffusion imaging (HARDI) at 4 Tesla. After performing a whole-brain fiber tractography using the full angular resolution of the diffusion scans, 70 cortical surface-based regions of interest were created from each individual’s co-registered anatomical data to compute graph metrics for all pairs of cortical regions. In graph theory analyses, subjects homozygous for the risk allele (CC) had lower characteristic path length, greater small-worldness and global efficiency in whole brain analyses, as well as greater eccentricity (maximum path length) in 60 of 70 nodes in regional analyses. These results were not reducible to differences in more commonly studied traits such as fiber density or fractional anisotropy. This is the first study to link graph theory metrics of brain structural connectivity to a common genetic variant linked with autism and will help us understand the neurobiology of circuits implicated in risk for autism.