M. Mallar Chakravarty

An anatomically comprehensive atlas of the adult human brain transcriptome

Neuroanatomically precise, genome-wide maps of transcript distributions are critical resources to complement genomic sequence data and to correlate functional and genetic brain architecture. Here we describe the generation and analysis of a transcriptional atlas of the adult human brain, comprising extensive histological analysis and comprehensive microarray profiling of ∼900 neuroanatomically precise subdivisions in two individuals. Transcriptional regulation varies enormously by anatomical location, with different regions and their constituent cell types displaying robust molecular signatures that are highly conserved between individuals. Analysis of differential gene expression and gene co-expression relationships demonstrates that brain-wide variation strongly reflects the distributions of major cell classes such as neurons, oligodendrocytes, astrocytes and microglia. Local neighbourhood relationships between fine anatomical subdivisions are associated with discrete neuronal subtypes and genes involved with synaptic transmission. The neocortex displays a relatively homogeneous transcriptional pattern, but with distinct features associated selectively with primary sensorimotor cortices and with enriched frontal lobe expression. Notably, the spatial topography of the neocortex is strongly reflected in its molecular topography—the closer two cortical regions, the more similar their transcriptomes. This freely accessible online data resource forms a high-resolution transcriptional baseline for neurogenetic studies of normal and abnormal human brain function.

Identification of common variants associated with human hippocampal and intracranial volumes

Identifying genetic variants influencing human brain structures may reveal new biological mechanisms underlying cognition and neuropsychiatric illness. The volume of the hippocampus is a biomarker of incipient Alzheimers disease and is reduced in schizophrenia, major depression and mesial temporal lobe epilepsy. Whereas many brain imaging phenotypes are highly heritable, identifying and replicating genetic influences has been difficult, as small effects and the high costs of magnetic resonance imaging (MRI) have led to underpowered studies. Here we report genome-wide association meta-analyses and replication for mean bilateral hippocampal, total brain and intracranial volumes from a large multinational consortium. The intergenic variant rs7294919 was associated with hippocampal volume (12q24.22; N = 21,151; P = 6.70 × 10−16) and the expression levels of the positional candidate gene TESC in brain tissue. Additionally, rs10784502, located within HMGA2, was associated with intracranial volume (12q14.3; N = 15,782; P = 1.12 × 10−12). We also identified a suggestive association with total brain volume at rs10494373 within DDR2 (1q23.3; N = 6,500; P = 5.81 × 10−7).

Fronto-striatal connections in the human brain: A probabilistic diffusion tractography study

Anatomical studies in animals have described multiple striatal circuits and suggested that sub-components of the striatum, although functionally related, project to distinct cortical areas. To date, anatomical investigations in humans have been limited by methodological constraints such that most of our knowledge of fronto-striatal networks relies on nonhuman primate studies. To better identify the fronto-striatal pathways in the human brain, we used Diffusion Tensor Imaging (DTI) tractography to reconstruct neural connections between the frontal cortex and the caudate nucleus and putamen in vivo. We demonstrate that the human caudate nucleus is interconnected with the prefrontal cortex, inferior and middle temporal gyrus, frontal eye fields, cerebellum and thalamus; the putamen is interconnected with the prefrontal cortex, primary motor area, primary somatosensory cortex, supplementary motor area, premotor area, cerebellum and thalamus. A connectivity-based seed classification analysis identified connections between the dorsolateral prefrontal areas (DLPFC) and the dorsal-posterior caudate nucleus and between the ventrolateral prefrontal areas (VLPFC) and the ventral-anterior caudate nucleus. For the putamen, connections exist between the supplementary motor area (SMA) and dorsal-posterior putamen while the premotor area projects to medial putamen, and the primary motor area to the lateral putamen. Identifying the anatomical organization of the fronto-striatal network has important implications for understanding basal ganglia function and associated disease processes.

Resting-state networks link invasive and noninvasive brain stimulation across diverse psychiatric and neurological diseases.

Significance Brain stimulation is a powerful treatment for an increasing number of psychiatric and neurological diseases, but it is unclear why certain stimulation sites work or where in the brain is the best place to stimulate to treat a given patient or disease. We found that although different types of brain stimulation are applied in different locations, targets used to treat the same disease most often are nodes in the same brain network. These results suggest that brain networks might be used to understand why brain stimulation works and to improve therapy by identifying the best places to stimulate the brain. Brain stimulation, a therapy increasingly used for neurological and psychiatric disease, traditionally is divided into invasive approaches, such as deep brain stimulation (DBS), and noninvasive approaches, such as transcranial magnetic stimulation. The relationship between these approaches is unknown, therapeutic mechanisms remain unclear, and the ideal stimulation site for a given technique is often ambiguous, limiting optimization of the stimulation and its application in further disorders. In this article, we identify diseases treated with both types of stimulation, list the stimulation sites thought to be most effective in each disease, and test the hypothesis that these sites are different nodes within the same brain network as defined by resting-state functional-connectivity MRI. Sites where DBS was effective were functionally connected to sites where noninvasive brain stimulation was effective across diseases including depression, Parkinsons disease, obsessive-compulsive disorder, essential tremor, addiction, pain, minimally conscious states, and Alzheimer’s disease. A lack of functional connectivity identified sites where stimulation was ineffective, and the sign of the correlation related to whether excitatory or inhibitory noninvasive stimulation was found clinically effective. These results suggest that resting-state functional connectivity may be useful for translating therapy between stimulation modalities, optimizing treatment, and identifying new stimulation targets. More broadly, this work supports a network perspective toward understanding and treating neuropsychiatric disease, highlighting the therapeutic potential of targeted brain network modulation.

A Genome-Wide Association Study Identifies Five Loci Influencing Facial Morphology in Europeans

Inter-individual variation in facial shape is one of the most noticeable phenotypes in humans, and it is clearly under genetic regulation; however, almost nothing is known about the genetic basis of normal human facial morphology. We therefore conducted a genome-wide association study for facial shape phenotypes in multiple discovery and replication cohorts, considering almost ten thousand individuals of European descent from several countries. Phenotyping of facial shape features was based on landmark data obtained from three-dimensional head magnetic resonance images (MRIs) and two-dimensional portrait images. We identified five independent genetic loci associated with different facial phenotypes, suggesting the involvement of five candidate genes—PRDM16, PAX3, TP63, C5orf50, and COL17A1—in the determination of the human face. Three of them have been implicated previously in vertebrate craniofacial development and disease, and the remaining two genes potentially represent novel players in the molecular networks governing facial development. Our finding at PAX3 influencing the position of the nasion replicates a recent GWAS of facial features. In addition to the reported GWA findings, we established links between common DNA variants previously associated with NSCL/P at 2p21, 8q24, 13q31, and 17q22 and normal facial-shape variations based on a candidate gene approach. Overall our study implies that DNA variants in genes essential for craniofacial development contribute with relatively small effect size to the spectrum of normal variation in human facial morphology. This observation has important consequences for future studies aiming to identify more genes involved in the human facial morphology, as well as for potential applications of DNA prediction of facial shape such as in future forensic applications.

Correlated gene expression supports synchronous activity in brain networks

Cooperating brain regions express similar genes When the brain is at rest, a number of distinct areas are functionally connected. They tend to be organized in networks. Richiardi et al. compared brain imaging and gene expression data to build computational models of these networks. These functional networks are underpinned by the correlated expression of a core set of 161 genes. In this set, genes coding for ion channels and other synaptic functions such as neurotransmitter release dominate. Science, this issue p. 1241 Gene expression is more similar than expected by chance in brain regions that are functionally connected. During rest, brain activity is synchronized between different regions widely distributed throughout the brain, forming functional networks. However, the molecular mechanisms supporting functional connectivity remain undefined. We show that functional brain networks defined with resting-state functional magnetic resonance imaging can be recapitulated by using measures of correlated gene expression in a post mortem brain tissue data set. The set of 136 genes we identify is significantly enriched for ion channels. Polymorphisms in this set of genes significantly affect resting-state functional connectivity in a large sample of healthy adolescents. Expression levels of these genes are also significantly associated with axonal connectivity in the mouse. The results provide convergent, multimodal evidence that resting-state functional networks correlate with the orchestrated activity of dozens of genes linked to ion channel activity and synaptic function.

Performing label-fusion-based segmentation using multiple automatically generated templates.

Classically, model‐based segmentation procedures match magnetic resonance imaging (MRI) volumes to an expertly labeled atlas using nonlinear registration. The accuracy of these techniques are limited due to atlas biases, misregistration, and resampling error. Multi‐atlas‐based approaches are used as a remedy and involve matching each subject to a number of manually labeled templates. This approach yields numerous independent segmentations that are fused using a voxel‐by‐voxel label‐voting procedure. In this article, we demonstrate how the multi‐atlas approach can be extended to work with input atlases that are unique and extremely time consuming to construct by generating a library of multiple automatically generated templates of different brains (MAGeT Brain). We demonstrate the efficacy of our method for the mouse and human using two different nonlinear registration algorithms (ANIMAL and ANTs). The input atlases consist a high‐resolution mouse brain atlas and an atlas of the human basal ganglia and thalamus derived from serial histological data. MAGeT Brain segmentation improves the identification of the mouse anterior commissure (mean Dice Kappa values (κ = 0.801), but may be encountering a ceiling effect for hippocampal segmentations. Applying MAGeT Brain to human subcortical structures improves segmentation accuracy for all structures compared to regular model‐based techniques (κ = 0.845, 0.752, and 0.861 for the striatum, globus pallidus, and thalamus, respectively). Experiments performed with three manually derived input templates suggest that MAGeT Brain can approach or exceed the accuracy of multi‐atlas label‐fusion segmentation (κ = 0.894, 0.815, and 0.895 for the striatum, globus pallidus, and thalamus, respectively). Hum Brain Mapp 34:2635–2654, 2013.

Longitudinal four-dimensional mapping of subcortical anatomy in human development

Significance Our spatiotemporal understanding of subcortical development in humans lags far behind that of the cortical sheet. This disparity ignores that developmental refinements and disruptions of complex behavior involve systems spanning both components of the brain. We begin redressing this imbalance by applying new techniques for striatal, pallidal and thalamic morphometry to large-scale longitudinal neuroimaging data extending from childhood through early adulthood. This work (i) establishes the curvilinear, sexual dimorphic and often protracted nature of global volume change within each structure, (ii) reveals profound spatiotemporal complexities in striatal, pallidal and thalamic maturation that are organized by the known topography of primate cortico-subcortical connectivity, and (iii) identifies focal sex differences in subcortical maturation that strike regions implicated in psychopathologies with an adolescent-emergent sex-bias. Growing access to large-scale longitudinal structural neuroimaging data has fundamentally altered our understanding of cortical development en route to human adulthood, with consequences for basic science, medicine, and public policy. In striking contrast, basic anatomical development of subcortical structures such as the striatum, pallidum, and thalamus has remained poorly described—despite these evolutionarily ancient structures being both intimate working partners of the cortical sheet and critical to diverse developmentally emergent skills and disorders. Here, to begin addressing this disparity, we apply methods for the measurement of subcortical volume and shape to 1,171 longitudinally acquired structural magnetic resonance imaging brain scans from 618 typically developing males and females aged 5–25 y. We show that the striatum, pallidum, and thalamus each follow curvilinear trajectories of volume change, which, for the striatum and thalamus, peak after cortical volume has already begun to decline and show a relative delay in males. Four-dimensional mapping of subcortical shape reveals that (i) striatal, pallidal, and thalamic domains linked to specific fronto-parietal association cortices contract with age whereas other subcortical territories expand, and (ii) each structure harbors hotspots of sexually dimorphic change over adolescence—with relevance for sex-biased mental disorders emerging in youth. By establishing the developmental dynamism, spatial heterochonicity, and sexual dimorphism of human subcortical maturation, these data bring our spatiotemporal understanding of subcortical development closer to that of the cortex—allowing evolutionary, basic, and clinical neuroscience to be conducted within a more comprehensive developmental framework.

Neurite density from magnetic resonance diffusion measurements at ultrahigh field: Comparison with light microscopy and electron microscopy

Due to its unique sensitivity to tissue microstructure, diffusion-weighted magnetic resonance imaging (MRI) has found many applications in clinical and fundamental science. With few exceptions, a more precise correspondence between physiological or biophysical properties and the obtained diffusion parameters remain uncertain due to lack of specificity. In this work, we address this problem by comparing diffusion parameters of a recently introduced model for water diffusion in brain matter to light microscopy and quantitative electron microscopy. Specifically, we compare diffusion model predictions of neurite density in rats to optical myelin staining intensity and stereological estimation of neurite volume fraction using electron microscopy. We find that the diffusion model describes data better and that its parameters show stronger correlation with optical and electron microscopy, and thus reflect myelinated neurite density better than the more frequently used diffusion tensor imaging (DTI) and cumulant expansion methods. Furthermore, the estimated neurite orientations capture dendritic architecture more faithfully than DTI diffusion ellipsoids.

Multi-atlas segmentation of the whole hippocampus and subfields using multiple automatically generated templates

INTRODUCTION Advances in image segmentation of magnetic resonance images (MRI) have demonstrated that multi-atlas approaches improve segmentation over regular atlas-based approaches. These approaches often rely on a large number of manually segmented atlases (e.g. 30-80) that take significant time and expertise to produce. We present an algorithm, MAGeT-Brain (Multiple Automatically Generated Templates), for the automatic segmentation of the hippocampus that minimises the number of atlases needed whilst still achieving similar agreement to multi-atlas approaches. Thus, our method acts as a reliable multi-atlas approach when using special or hard-to-define atlases that are laborious to construct. METHOD MAGeT-Brain works by propagating atlas segmentations to a template library, formed from a subset of target images, via transformations estimated by nonlinear image registration. The resulting segmentations are then propagated to each target image and fused using a label fusion method. We conduct two separate Monte Carlo cross-validation experiments comparing MAGeT-Brain and basic multi-atlas whole hippocampal segmentation using differing atlas and template library sizes, and registration and label fusion methods. The first experiment is a 10-fold validation (per parameter setting) over 60 subjects taken from the Alzheimers Disease Neuroimaging Database (ADNI), and the second is a five-fold validation over 81 subjects having had a first episode of psychosis. In both cases, automated segmentations are compared with manual segmentations following the Pruessner-protocol. Using the best settings found from these experiments, we segment 246 images of the ADNI1:Complete 1Yr 1.5 T dataset and compare these with segmentations from existing automated and semi-automated methods: FSL FIRST, FreeSurfer, MAPER, and SNT. Finally, we conduct a leave-one-out cross-validation of hippocampal subfield segmentation in standard 3T T1-weighted images, using five high-resolution manually segmented atlases (Winterburn et al., 2013). RESULTS In the ADNI cross-validation, using 9 atlases MAGeT-Brain achieves a mean Dices Similarity Coefficient (DSC) score of 0.869 with respect to manual whole hippocampus segmentations, and also exhibits significantly lower variability in DSC scores than multi-atlas segmentation. In the younger, psychosis dataset, MAGeT-Brain achieves a mean DSC score of 0.892 and produces volumes which agree with manual segmentation volumes better than those produced by the FreeSurfer and FSL FIRST methods (mean difference in volume: 80 mm(3), 1600 mm(3), and 800 mm(3), respectively). Similarly, in the ADNI1:Complete 1Yr 1.5 T dataset, MAGeT-Brain produces hippocampal segmentations well correlated (r>0.85) with SNT semi-automated reference volumes within disease categories, and shows a conservative bias and a mean difference in volume of 250 mm(3) across the entire dataset, compared with FreeSurfer and FSL FIRST which both overestimate volume differences by 2600 mm(3) and 2800 mm(3) on average, respectively. Finally, MAGeT-Brain segments the CA1, CA4/DG and subiculum subfields on standard 3T T1-weighted resolution images with DSC overlap scores of 0.56, 0.65, and 0.58, respectively, relative to manual segmentations. CONCLUSION We demonstrate that MAGeT-Brain produces consistent whole hippocampal segmentations using only 9 atlases, or fewer, with various hippocampal definitions, disease populations, and image acquisition types. Additionally, we show that MAGeT-Brain identifies hippocampal subfields in standard 3T T1-weighted images with overlap scores comparable to competing methods.