A commonly carried genetic variant in the delta opioid receptor gene, OPRD1, is associated with smaller regional brain volumes: replication in elderly and young populations.

Delta opioid receptors are implicated in a variety of psychiatric and neurological disorders. These receptors play a key role in the reinforcing properties of drugs of abuse, and polymorphisms in OPRD1 (the gene encoding delta opioid receptors) are associated with drug addiction. Delta opioid receptors are also involved in protecting neurons against hypoxic and ischemic stress. Here, we first examined a large sample of 738 elderly participants with neuroimaging and genetic data from the Alzheimer's Disease Neuroimaging Initiative. We hypothesized that common variants in OPRD1 would be associated with differences in brain structure, particularly in regions relevant to addictive and neurodegenerative disorders. One very common variant (rs678849) predicted differences in regional brain volumes. We replicated the association of this single-nucleotide polymorphism with regional tissue volumes in a large sample of young participants in the Queensland Twin Imaging study. Although the same allele was associated with reduced volumes in both cohorts, the brain regions affected differed between the two samples. In healthy elderly, exploratory analyses suggested that the genotype associated with reduced brain volumes in both cohorts may also predict cerebrospinal fluid levels of neurodegenerative biomarkers, but this requires confirmation. If opiate receptor genetic variants are related to individual differences in brain structure, genotyping of these variants may be helpful when designing clinical trials targeting delta opioid receptors to treat neurological disorders.

White matter microstructural abnormalities in bipolar disorder: A whole brain diffusion tensor imaging study.

BACKGROUND: Bipolar disorder (BD) is a chronic mental illness characterized by severe disruptions in mood and cognition. Diffusion tensor imaging (DTI) studies suggest that white matter (WM) tract abnormalities may contribute to the clinical hallmarks of the disorder. Using DTI and whole brain voxel-based analysis, we mapped the profile of WM anomalies in BD. All patients in our sample were euthymic and lithium free when scanned. METHODS: Diffusion-weighted and T1-weighted structural brain images were acquired from 23 lithium-free euthymic subjects with bipolar I disorder and 19 age- and sex-matched healthy control subjects on a 1.5 T MRI scanner. Scans were processed to provide measures of fractional anisotropy (FA) and mean and radial diffusivity (MD and RD) at each WM voxel, and processed scans were nonlinearly aligned to a customized brain imaging template for statistical group comparisons. RESULTS: Relative to controls, the bipolar group showed widespread regions of lower FA, including the corpus callosum, cortical and thalamic association fibers. MD and RD were abnormally elevated in patients in many of these same regions. CONCLUSIONS: Our findings agree with prior reports of WM abnormalities in the corpus callosum and further link a bipolar diagnosis with structural abnormalities of the tapetum, fornix and stria terminalis. Future studies assessing the diagnostic specificity and prognostic implications of these abnormalities would be of interest.

Relation between variants in the neurotrophin receptor gene, NTRK3, and white matter integrity in healthy young adults.

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.

Gene network effects on brain microstructure and intellectual performance identified in 472 twins.

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 approximately 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 nontwin siblings (mean age: 23.7 ± 2.1 SD years; 193 male/279 female). 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 computationally more 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 and the integrity of major white matter tracts, such as the callosal genu and splenium, cingulum, optic radiations, and the superior longitudinal fasciculus.

Relationship of a variant in the NTRK1 gene to white matter microstructure in young adults.

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 fractional anisotropy 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 reproducibility of results. 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.

Common Alzheimer's disease risk variant within the CLU gene affects white matter microstructure in young adults.

There is a strong genetic risk for late-onset Alzheimer's disease (AD), but so far few gene variants have been identified that reliably contribute to that risk. A newly confirmed genetic risk allele C of the clusterin (CLU) gene variant rs11136000 is carried by ∼88% of Caucasians. The C allele confers a 1.16 greater odds of developing late-onset AD than the T allele. AD patients have reductions in regional white matter integrity. We evaluated whether the CLU risk variant was similarly associated with lower white matter integrity in healthy young humans. Evidence of early brain differences would offer a target for intervention decades before symptom onset. We scanned 398 healthy young adults (mean age, 23.6 ± 2.2 years) with diffusion tensor imaging, a variation of magnetic resonance imaging sensitive to white matter integrity in the living brain. We assessed genetic associations using mixed-model regression at each point in the brain to map the profile of these associations with white matter integrity. Each C allele copy of the CLU variant was associated with lower fractional anisotropy--a widely accepted measure of white matter integrity--in multiple brain regions, including several known to degenerate in AD. These regions included the splenium of the corpus callosum, the fornix, cingulum, and superior and inferior longitudinal fasciculi in both brain hemispheres. Young healthy carriers of the CLU gene risk variant showed a distinct profile of lower white matter integrity that may increase vulnerability to developing AD later in life.

The contribution of genes to cortical thickness and volume.

We analyzed brain MRI data from 372 young adult twins to identify cortical regions in which gray matter thickness and volume are influenced by genetics. This was achieved using an A/C/E structural equation model that divides the variance of these traits, at each point on the cortex, into additive genetic (A), shared (C), and unique environmental (E) components. A strong genetic influence was found in frontal and parietal regions. In addition, we correlated cortical thickness with full-scale intelligence quotient for comparison with the A/C/E maps, and several regions where cortical structure was correlated with intelligence quotient are under genetic control. These cortical measures may be useful phenotypes to narrow the search for quantitative trait loci influencing brain structure.

BDNF gene effects on brain circuitry replicated in 455 twins.

Brain-derived neurotrophic factor (BDNF) plays a key role in learning and memory, but its effects on the fiber architecture of the living brain are unknown. We genotyped 455 healthy adult twins and their non-twin siblings (188 males/267 females; age: 23.7±2.1 years, mean±SD) and scanned them with high angular resolution diffusion tensor imaging (DTI), to assess how the BDNF Val66Met polymorphism affects white matter microstructure. By applying genetic association analysis to every 3D point in the brain images, we found that the Val-BDNF genetic variant was associated with lower white matter integrity in the splenium of the corpus callosum, left optic radiation, inferior fronto-occipital fasciculus, and superior corona radiata. Normal BDNF variation influenced the association between subjects' performance intellectual ability (as measured by Object Assembly subtest) and fiber integrity (as measured by fractional anisotropy; FA) in the callosal splenium, and pons. BDNF gene may affect the intellectual performance by modulating the white matter development. This combination of genetic association analysis and large-scale diffusion imaging directly relates a specific gene to the fiber microstructure of the living brain and to human intelligence.

Genetic influences on brain asymmetry: a DTI study of 374 twins and siblings.

Brain asymmetry, or the structural and functional specialization of each brain hemisphere, has fascinated neuroscientists for over a century. Even so, genetic and environmental factors that influence brain asymmetry are largely unknown. Diffusion tensor imaging (DTI) now allows asymmetry to be studied at a microscopic scale by examining differences in fiber characteristics across hemispheres rather than differences in structure shapes and volumes. Here we analyzed 4Tesla DTI scans from 374 healthy adults, including 60 monozygotic twin pairs, 45 same-sex dizygotic pairs, and 164 mixed-sex DZ twins and their siblings; mean age: 24.4years+/-1.9 SD). All DTI scans were nonlinearly aligned to a geometrically-symmetric, population-based image template. We computed voxel-wise maps of significant asymmetries (left/right differences) for common diffusion measures that reflect fiber integrity (fractional and geodesic anisotropy; FA, GA and mean diffusivity, MD). In quantitative genetic models computed from all same-sex twin pairs (N=210 subjects), genetic factors accounted for 33% of the variance in asymmetry for the inferior fronto-occipital fasciculus, 37% for the anterior thalamic radiation, and 20% for the forceps major and uncinate fasciculus (all L>R). Shared environmental factors accounted for around 15% of the variance in asymmetry for the cortico-spinal tract (R>L) and about 10% for the forceps minor (L>R). Sex differences in asymmetry (men>women) were significant, and were greatest in regions with prominent FA asymmetries. These maps identify heritable DTI-derived features, and may empower genome-wide searches for genetic polymorphisms that influence brain asymmetry.

A GENETIC ANALYSIS OF CORTICAL THICKNESS IN 372 TWINS.

Imaging genetics is a new field of neuroscience that blends methods from computational anatomy and quantitative genetics to identify genetic influences on brain structure and function. Here we analyzed brain MRI data from 372 young adult twins to identify cortical regions in which gray matter volume is influenced by genetic differences across subjects. Thickness maps, reconstructed from surface models of the cortical gray/white and gray/CSF interfaces, were smoothed with a 25 mm FWHM kernel and automatically parcellated into 34 regions of interest per hemisphere. In structural equation models fitted to volume values at each surface vertex, we computed components of variance due to additive genetic (A), shared (C) and unique (E) environmental factors, and tested their significance. Cortical regions in the vicinity of the perisylvian language cortex, and at the frontal and temporal poles, showed significant additive genetic variance, suggesting that volume measures from these regions may provide quantitative phenotypes to narrow the search for quantitative trait loci that influence brain structure.

STATISTICALLY ASSISTED FLUID IMAGE REGISTRATION ALGORITHM - SAFIRA.

In this paper, we develop and validate a new Statistically Assisted Fluid Registration Algorithm (SAFIRA) for brain images. A non-statistical version of this algorithm was first implemented in [2] and re-formulated using Lagrangian mechanics in [3]. Here we extend this algorithm to 3D: given 3D brain images from a population, vector fields and their corresponding deformation matrices are computed in a first round of registrations using the non-statistical implementation. Covariance matrices for both the deformation matrices and the vector fields are then obtained and incorporated (separately or jointly) in the regularizing (i.e., the non-conservative Lagrangian) terms, creating four versions of the algorithm. We evaluated the accuracy of each algorithm variant using the manually labeled LPBA40 dataset, which provides us with ground truth anatomical segmentations. We also compared the power of the different algorithms using tensor-based morphometry -a technique to analyze local volumetric differences in brain structure-applied to 46 3D brain scans from healthy monozygotic twins.

MULTIVARIATE VARIANCE-COMPONENTS ANALYSIS IN DTI.

Twin studies are a major research direction in imaging genetics, a new field, which combines algorithms from quantitative genetics and neuroimaging to assess genetic effects on the brain. In twin imaging studies, it is common to estimate the intraclass correlation (ICC), which measures the resemblance between twin pairs for a given phenotype. In this paper, we extend the commonly used Pearson correlation to a more appropriate definition, which uses restricted maximum likelihood methods (REML). We computed proportion of phenotypic variance due to additive (A) genetic factors, common (C) and unique (E) environmental factors using a new definition of the variance components in the diffusion tensor-valued signals. We applied our analysis to a dataset of Diffusion Tensor Images (DTI) from 25 identical and 25 fraternal twin pairs. Differences between the REML and Pearson estimators were plotted for different sample sizes, showing that the REML approach avoids severe biases when samples are smaller. Measures of genetic effects were computed for scalar and multivariate diffusion tensor derived measures including the geodesic anisotropy (tGA) and the full diffusion tensors (DT), revealing voxel-wise genetic contributions to brain fiber microstructure.

How does angular resolution affect diffusion imaging measures?

A key question in diffusion imaging is how many diffusion-weighted images suffice to provide adequate signal-to-noise ratio (SNR) for studies of fiber integrity. Motion, physiological effects, and scan duration all affect the achievable SNR in real brain images, making theoretical studies and simulations only partially useful. We therefore scanned 50 healthy adults with 105-gradient high-angular resolution diffusion imaging (HARDI) at 4T. From gradient image subsets of varying size (6

Mapping the regional influence of genetics on brain structure variability--a tensor-based morphometry study.

Genetic and environmental factors influence brain structure and function profoundly. The search for heritable anatomical features and their influencing genes would be accelerated with detailed 3D maps showing the degree to which brain morphometry is genetically determined. As part of an MRI study that will scan 1150 twins, we applied Tensor-Based Morphometry to compute morphometric differences in 23 pairs of identical twins and 23 pairs of same-sex fraternal twins (mean age: 23.8+/-1.8 SD years). All 92 twins' 3D brain MRI scans were nonlinearly registered to a common space using a Riemannian fluid-based warping approach to compute volumetric differences across subjects. A multi-template method was used to improve volume quantification. Vector fields driving each subject's anatomy onto the common template were analyzed to create maps of local volumetric excesses and deficits relative to the standard template. Using a new structural equation modeling method, we computed the voxelwise proportion of variance in volumes attributable to additive (A) or dominant (D) genetic factors versus shared environmental (C) or unique environmental factors (E). The method was also applied to various anatomical regions of interest (ROIs). As hypothesized, the overall volumes of the brain, basal ganglia, thalamus, and each lobe were under strong genetic control; local white matter volumes were mostly controlled by common environment. After adjusting for individual differences in overall brain scale, genetic influences were still relatively high in the corpus callosum and in early-maturing brain regions such as the occipital lobes, while environmental influences were greater in frontal brain regions that have a more protracted maturational time-course.

Genetics of brain fiber architecture and intellectual performance.

The study is the first to analyze genetic and environmental factors that affect brain fiber architecture and its genetic linkage with cognitive function. We assessed white matter integrity voxelwise using diffusion tensor imaging at high magnetic field (4 Tesla), in 92 identical and fraternal twins. White matter integrity, quantified using fractional anisotropy (FA), was used to fit structural equation models (SEM) at each point in the brain, generating three-dimensional maps of heritability. We visualized the anatomical profile of correlations between white matter integrity and full-scale, verbal, and performance intelligence quotients (FIQ, VIQ, and PIQ). White matter integrity (FA) was under strong genetic control and was highly heritable in bilateral frontal (a(2)=0.55, p=0.04, left; a(2)=0.74, p=0.006, right), bilateral parietal (a(2)=0.85, p<0.001, left; a(2)=0.84, p<0.001, right), and left occipital (a(2)=0.76, p=0.003) lobes, and was correlated with FIQ and PIQ in the cingulum, optic radiations, superior fronto-occipital fasciculus, internal capsule, callosal isthmus, and the corona radiata (p=0.04 for FIQ and p=0.01 for PIQ, corrected for multiple comparisons). In a cross-trait mapping approach, common genetic factors mediated the correlation between IQ and white matter integrity, suggesting a common physiological mechanism for both, and common genetic determination. These genetic brain maps reveal heritable aspects of white matter integrity and should expedite the discovery of single-nucleotide polymorphisms affecting fiber connectivity and cognition.

REDUCING STRUCTURAL VARIATION TO DETERMINE THE GENETICS OF WHITE MATTER INTEGRITY ACROSS HEMISPHERES - A DTI STUDY OF 100 TWINS.

Studies of cerebral asymmetry can open doors to understanding the functional specialization of each brain hemisphere, and how this is altered in disease. Here we examined hemispheric asymmetries in fiber architecture using diffusion tensor imaging (DTI) in 100 subjects, using high-dimensional fluid warping to disentangle shape differences from measures sensitive to myelination. Confounding effects of purely structural asymmetries were reduced by using co-registered structural images to fluidly warp 3D maps of fiber characteristics (fractional and geodesic anisotropy) to a structurally symmetric minimal deformation template (MDT). We performed a quantitative genetic analysis on 100 subjects to determine whether the sources of the remaining signal asymmetries were primarily genetic or environmental. A twin design was used to identify the heritable features of fiber asymmetry in various regions of interest, to further assist in the discovery of genes influencing brain micro-architecture and brain lateralization. Genetic influences and left/right asymmetries were detected in the fiber architecture of the frontal lobes, with minor differences depending on the choice of registration template.

THE MULTIVARIATE A/C/E MODEL AND THE GENETICS OF FIBER ARCHITECTURE.

We present a new algorithm to compute the voxel-wise genetic contribution to brain fiber microstructure using diffusion tensor imaging (DTI) in a dataset of 25 pairs of monozygotic (MZ) twins and 25 pairs of dizygotic (DZ) twins. First, the structural and DT scans were linearly co-registered. The structural MR scans were nonlinear mapped via a 3D fluid transformation to a geometrically centered mean template, and the deformation fields were applied to the DTI volumes. After tensor re-orientation to realign them to the anatomy, we computed several scalar and multivariate DT-derived measures including the geodesic anisotropy (GA), the tensor eigenvalues and the full diffusion tensors. A covariance-weighted distance was found between twins in the Log-Euclidean framework [2], and used as input to a maximum-likelihood based algorithm to compute the contributions from genetics (A), common environmental factors (C) and unique environmental ones (E) to fiber architecture. Quantitative genetic studies can make use of the full information in the diffusion tensor, using covariance weighted distances and statistics on the tensor manifold.

A LAGRANGIAN FORMULATION FOR STATISTICAL FLUID REGISTRATION.

We defined a new statistical fluid registration method with Lagrangian mechanics. Although several authors have suggested that empirical statistics on brain variation should be incorporated into the registration problem, few algorithms have included this information and instead use regularizers that guarantee diffeomorphic mappings. Here we combine the advantages of a large-deformation fluid matching approach with empirical statistics on population variability in anatomy. We reformulated the Riemannian fluid algorithm developed in [4], and used a Lagrangian framework to incorporate 0 th and 1 rst order statistics in the regularization process. 92 2D midline corpus callosum traces from a twin MRI database were fluidly registered using the non-statistical version of the algorithm (algorithm 0), giving initial vector fields and deformation tensors. Covariance matrices were computed for both distributions and incorporated either separately (algorithm 1 and algorithm 2) or together (algorithm 3) in the registration. We computed heritability maps and two vector and tensor-based distances to compare the power and the robustness of the algorithms.

Mapping genetic influences on ventricular structure in twins.

Despite substantial progress in measuring the anatomical and functional variability of the human brain, little is known about the genetic and environmental causes of these variations. Here we developed an automated system to visualize genetic and environmental effects on brain structure in large brain MRI databases. We applied our multi-template segmentation approach termed "Multi-Atlas Fluid Image Alignment" to fluidly propagate hand-labeled parameterized surface meshes, labeling the lateral ventricles, in 3D volumetric MRI scans of 76 identical (monozygotic, MZ) twins (38 pairs; mean age=24.6 (SD=1.7)); and 56 same-sex fraternal (dizygotic, DZ) twins (28 pairs; mean age=23.0 (SD=1.8)), scanned as part of a 5-year research study that will eventually study over 1000 subjects. Mesh surfaces were averaged within subjects to minimize segmentation error. We fitted quantitative genetic models at each of 30,000 surface points to measure the proportion of shape variance attributable to (1) genetic differences among subjects, (2) environmental influences unique to each individual, and (3) shared environmental effects. Surface-based statistical maps, derived from path analysis, revealed patterns of heritability, and their significance, in 3D. Path coefficients for the 'ACE' model that best fitted the data indicated significant contributions from genetic factors (A=7.3%), common environment (C=38.9%) and unique environment (E=53.8%) to lateral ventricular volume. Earlier-maturing occipital horn regions may also be more genetically influenced than later-maturing frontal regions. Maps visualized spatially-varying profiles of environmental versus genetic influences. The approach shows promise for automatically measuring gene-environment effects in large image databases.

Tensor-based analysis of genetic influences on brain integrity using DTI in 100 twins.

Information from the full diffusion tensor (DT) was used to compute voxel-wise genetic contributions to brain fiber microstructure. First, we designed a new multivariate intraclass correlation formula in the log-Euclidean framework. We then analyzed used the full multivariate structure of the tensor in a multivariate version of a voxel-wise maximum-likelihood structural equation model (SEM) that computes the variance contributions in the DTs from genetic (A), common environmental (C) and unique environmental (E) factors. Our algorithm was tested on DT images from 25 identical and 25 fraternal twin pairs. After linear and fluid registration to a mean template, we computed the intraclass correlation and Falconer's heritability statistic for several scalar DT-derived measures and for the full multivariate tensors. Covariance matrices were found from the DTs, and inputted into SEM. Analyzing the full DT enhanced the detection of A and C effects. This approach should empower imaging genetics studies that use DTI.