Obesity gene NEGR1 associated with white matter integrity in healthy young adults

Obesity is a crucial public health issue in developed countries, with implications for cardiovascular and brain health as we age. A number of commonly-carried genetic variants are associated with obesity. Here we aim to see whether variants in obesity-associated genes - NEGR1, FTO, MTCH2, MC4R, LRRN6C, MAP2K5, FAIM2, SEC16B, ETV5, BDNF- AS, ATXN2L, ATP2A1, KCTD15, and TNN13K - are associated with white matter microstructural properties, assessed by high angular resolution diffusion imaging (HARDI) in young healthy adults between 20 and 30. years of age from the Queensland Twin Imaging study (QTIM). We began with a multi-locus approach testing how a number of common genetic risk factors for obesity at the single nucleotide polymorphism (SNP) level may jointly influence white matter integrity throughout the brain and found a wide spread genetic effect. Risk allele rs2815752 in NEGR1 was most associated with lower white matter integrity across a substantial portion of the brain. Across the area of significance in the bilateral posterior corona radiata, each additional copy of the risk allele was associated with a 2.2% lower average FA. This is the first study to find an association between an obesity risk gene and differences in white matter integrity. As our subjects were young and healthy, our results suggest that NEGR1 has effects on brain structure independent of its effect on obesity.

Bivariate genome-wide association study of genetically correlated neuroimaging phenotypes from DTI and MRI through a seemingly unrelated regression model

Large multisite efforts (e.g., the ENIGMA Consortium), have shown that neuroimaging traits including tract integrity (from DTI fractional anisotropy, FA) and subcortical volumes (from T1-weighted scans) are highly heritable and promising phenotypes for discovering genetic variants associated with brain structure. However, genetic correlations (rg) among measures from these different modalities for mapping the human genome to the brain remain unknown. Discovering these correlations can help map genetic and neuroanatomical pathways implicated in development and inherited risk for disease. We use structural equation models and a twin design to find rg between pairs of phenotypes extracted from DTI and MRI scans. When controlling for intracranial volume, the caudate as well as related measures from the limbic system - hippocampal volume - showed high rg with the cingulum FA. Using an unrelated sample and a Seemingly Unrelated Regression model for bivariate analysis of this connection, we show that a multivariate GWAS approach may be more promising for genetic discovery than a univariate approach applied to each trait separately.

Diffusion imaging protocol effects on genetic associations

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.

Genome-wide scan of healthy human connectome discovers SPON1 gene variant influencing dementia severity

Aberrant connectivity is implicated in many neurological and psychiatric disorders, including Alzheimer's disease and schizophrenia. However, other than a few disease-associated candidate genes, we know little about the degree to which genetics play a role in the brain networks; we know even less about specific genes that influence brain connections. Twin and family-based studies can generate estimates of overall genetic influences on a trait, but genome-wide association scans (GWASs) can screen the genome for specific variants influencing the brain or risk for disease. To identify the heritability of various brain connections, we scanned healthy young adult twins with high-field, highangular resolution diffusion MRI. We adapted GWASs to screen the brain's connectivity pattern, allowing us to discover genetic variants that affect the human brain's wiring. The association of connectivity with the SPON1 variant at rs2618516 on chromosome 11 (11p15.2) reached connectome-wide, genome-wide significance after stringent statistical corrections were enforced, and it was replicated in an independent subsample. rs2618516 was shown to affect brain structure in an elderly population with varying degrees of dementia. Older people who carried the connectivity variant had significantly milder clinical dementia scores and lower risk of Alzheimer's disease. As a posthoc analysis, we conducted GWASs on several organizational and topological network measures derived from the matrices to discover variants in and around genes associated with autism (MACROD2), development (NEDD4), and mental retardation (UBE2A) significantly associated with connectivity. Connectome-wide, genome-wide screening offers substantial promise to discover genes affecting brain connectivity and risk for brain diseases.

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.

The contribution of genes to cortical thickness and volume

We analyzed brain MRI data from 372 young adult twins toidentify 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 influencewas found in frontal and parietal regions. Inaddition, 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 searchfor quantitative trait lociinfluencing brain structure.

Predicting white matter integrity from multiple common genetic variants

Several common genetic variants have recently been discovered that appear to influence white matter microstructure, as measured by diffusion tensor imaging (DTI). Each genetic variant explains only a small proportion of the variance in brain microstructure, so we set out to explore their combined effect on the white matter integrity of the corpus callosum. We measured six common candidate single-nucleotide polymorphisms (SNPs) in the COMT, NTRK1, BDNF, ErbB4, CLU, and HFE genes, and investigated their individual and aggregate effects on white matter structure in 395 healthy adult twins and siblings (age: 20-30 years). All subjects were scanned with 4-tesla 94-direction high angular resolution diffusion imaging. When combined using mixed-effects linear regression, a joint model based on five of the candidate SNPs (COMT, NTRK1, ErbB4, CLU, and HFE) explained ∼ 6% of the variance in the average fractional anisotropy (FA) of the corpus callosum. This predictive model had detectable effects on FA at 82% of the corpus callosum voxels, including the genu, body, and splenium. Predicting the brain's fiber microstructure from genotypes may ultimately help in early risk assessment, and eventually, in personalized treatment for neuropsychiatric disorders in which brain integrity and connectivity are affected.

Best individual template selection from deformation tensor minimization

We study the influence of the choice of template in tensor-based morphometry. Using 3D brain MR images from 10 monozygotic twin pairs, we defined a tensor-based distance in the log-Euclidean framework [1] between each image pair in the study. Relative to this metric, twin pairs were found to be closer to each other on average than random pairings, consistent with evidence that brain structure is under strong genetic control. We also computed the intraclass correlation and associated permutation p-value at each voxel for the determinant of the Jacobian matrix of the transformation. The cumulative distribution function (cdf) of the p-values was found at each voxel for each of the templates and compared to the null distribution. Surprisingly, there was very little difference between CDFs of statistics computed from analyses using different templates. As the brain with least log-Euclidean deformation cost, the mean template defined here avoids the blurring caused by creating a synthetic image from a population, and when selected from a large population, avoids bias by being geometrically centered, in a metric that is sensitive enough to anatomical similarity that it can even detect genetic affinity among anatomies.

A new combined surface and volume registration

3D registration of brain MRI data is vital for many medical imaging applications. However, purely intensitybased approaches for inter-subject matching of brain structure are generally inaccurate in cortical regions, due to the highly complex network of sulci and gyri, which vary widely across subjects. Here we combine a surfacebased cortical registration with a 3D fluid one for the first time, enabling precise matching of cortical folds, but allowing large deformations in the enclosed brain volume, which guarantee diffeomorphisms. This greatly improves the matching of anatomy in cortical areas. The cortices are segmented and registered with the software Freesurfer. The deformation field is initially extended to the full 3D brain volume using a 3D harmonic mapping that preserves the matching between cortical surfaces. Finally, these deformation fields are used to initialize a 3D Riemannian fluid registration algorithm, that improves the alignment of subcortical brain regions. We validate this method on an MRI dataset from 92 healthy adult twins. Results are compared to those based on volumetric registration without surface constraints; the resulting mean templates resolve consistent anatomical features both subcortically and at the cortex, suggesting that the approach is well-suited for cross-subject integration of functional and anatomic data.

Neuroimaging and genetics: Exploring, searching, and finding

This issue on the genetics of brain imaging phenotypes is a celebration of the happy marriage between two of science's highly interesting fields: neuroscience and genetics. The articles collected here are ample evidence that a good deal of synergy exists in this marriage. A wide selection of papers is presented that provide many different perspectives on how genes cause variation in brain structure and function, which in turn influence behavioral phenotypes (including psychopathology). They are examples of the many different methodologies in contemporary genetics and neuroscience research. Genetic methodology includes genome-wide association (GWA), candidate-gene association, and twin studies. Sources of data on brain phenotypes include cortical gray matter (GM) structural/volumetric measures from magnetic resonance imaging (MRI); white matter (WM) measures from diffusion tensor imaging (DTI), such as fractional anisotropy; functional- (activity-) based measures from electroencephalography (EEG), and functional MRI (fMRI). Together, they reflect a combination of scientific fields that have seen great technological advances, whether it is the single-nucleotide polymorphism (SNP) array in genetics, the increasingly high-resolution MRI imaging, or high angular resolution diffusion imaging technique for measuring WM connective properties.

Discovery and replication of dopamine-related gene effects on caudate volume in young and elderly populations (N1198) using genome-wide search

The caudate is a subcortical brain structure implicated in many common neurological and psychiatric disorders. To identify specific genes associated with variations in caudate volume, structural magnetic resonance imaging and genome-wide genotypes were acquired from two large cohorts, the Alzheimer's Disease NeuroImaging Initiative (ADNI; N=734) and the Brisbane Adolescent/Young Adult Longitudinal Twin Study (BLTS; N=464). In a preliminary analysis of heritability, around 90% of the variation in caudate volume was due to genetic factors. We then conducted genome-wide association to find common variants that contribute to this relatively high heritability. Replicated genetic association was found for the right caudate volume at single-nucleotide polymorphism rs163030 in the ADNI discovery sample (P=2.36 × 10 -6) and in the BLTS replication sample (P=0.012). This genetic variation accounted for 2.79 and 1.61% of the trait variance, respectively. The peak of association was found in and around two genes, WDR41 and PDE8B, involved in dopamine signaling and development. In addition, a previously identified mutation in PDE8B causes a rare autosomal-dominant type of striatal degeneration. Searching across both samples offers a rigorous way to screen for genes consistently influencing brain structure at different stages of life. Variants identified here may be relevant to common disorders affecting the caudate.

Mapping hippocampal and ventricular change in Alzheimer disease

We developed an anatomical mapping technique to detect hippocampal and ventricular changes in Alzheimer disease (AD). The resulting maps are sensitive to longitudinal changes in brain structure as the disease progresses. An anatomical surface modeling approach was combined with surface-based statistics to visualize the region and rate of atrophy in serial MRI scans and isolate where these changes link with cognitive decline. Fifty-two high-resolution MRI scans were acquired from 12 AD patients (age: 68.4 ± 1.9 years) and 14 matched controls (age: 71.4 ± 0.9 years), each scanned twice (2.1 ± 0.4 years apart). 3D parametric mesh models of the hippocampus and temporal horns were created in sequential scans and averaged across subjects to identify systematic patterns of atrophy. As an index of radial atrophy, 3D distance fields were generated relating each anatomical surface point to a medial curve threading down the medial axis of each structure. Hippocampal atrophic rates and ventricular expansion were assessed statistically using surface-based permutation testing and were faster in AD than in controls. Using color-coded maps and video sequences, these changes were visualized as they progressed anatomically over time. Additional maps localized regions where atrophic changes linked with cognitive decline. Temporal horn expansion maps were more sensitive to AD progression than maps of hippocampal atrophy, but both maps correlated with clinical deterioration. These quantitative, dynamic visualizations of hippocampal atrophy and ventricular expansion rates in aging and AD may provide a promising measure to track AD progression in drug trials.

Serum cholesterol and variant in cholesterol-related gene CETP predict white matter microstructure

Several common genetic variants influence cholesterol levels, which play a key role in overall health. Myelin synthesis and maintenance are highly sensitive to cholesterol concentrations, and abnormal cholesterol levels increase the risk for various brain diseases, including Alzheimer's disease. We report significant associations between higher serum cholesterol (CHOL) and high-density lipoprotein levels and higher fractional anisotropy in 403 young adults (23.8 ± 2.4years) scanned with diffusion imaging and anatomic magnetic resonance imaging at 4Tesla. By fitting a multi-locus genetic model within white matter areas associated with CHOL, we found that a set of 18 cholesterol-related, single-nucleotide polymorphisms implicated in Alzheimer's disease risk predicted fractional anisotropy. We focused on the single-nucleotide polymorphism with the largest individual effects, CETP (rs5882), and found that increased G-allele dosage was associated with higher fractional anisotropy and lower radial and mean diffusivities in voxel-wise analyses of the whole brain. A follow-up analysis detected white matter associations with rs5882 in the opposite direction in 78 older individuals (74.3 ± 7.3years). Cholesterol levels may influence white matter integrity, and cholesterol-related genes may exert age-dependent effects on the brain.

Identifying candidate gene effects by restricting search space in a multivariate genetic analysis of white matter microstructure

Several genetic variants are thought to influence white matter (WM) integrity, measured with diffusion tensor imaging (DTI). Voxel based methods can test genetic associations, but heavy multiple comparisons corrections are required to adjust for searching the whole brain and for all genetic variants analyzed. Thus, genetic associations are hard to detect even in large studies. Using a recently developed multi-SNP analysis, we examined the joint predictive power of a group of 18 cholesterol-related single nucleotide polymorphisms (SNPs) on WM integrity, measured by fractional anisotropy. To boost power, we limited the analysis to brain voxels that showed significant associations with total serum cholesterol levels. From this space, we identified two genes with effects that replicated in individual voxel-wise analyses of the whole brain. Multivariate analyses of genetic variants on a reduced anatomical search space may help to identify SNPs with strongest effects on the brain from a broad panel of genes.

Glutamatergic Modulation of Cue-Induced Drug-Seeking Behavior in the Rat

The characteristics of drug addiction include compulsive drug use despite negative consequences and re-occurring relapses, returns to drug use after a period of abstinence. Therefore, relapse prevention is one of the major challenges for the treatment of drug addiction. There are three main factors capable of inducing craving for drugs and triggering relapse long after cessation of drug use and dissipation of physical withdrawal signs: stress, re-exposure to the drug, and environmental stimuli (cues) that have been previously associated with drug use. The neurotransmitters dopamine and glutamate have been implicated in the modulation of drug-seeking behavior. The aim of this project was to examine the role of glutamatergic neurotransmission in relapse triggered by conditioned drug-associated stimuli. The focus was on clarifying whether relapse to drug seeking can be attenuated by blockade of glutamate receptors. In addition, as the nucleus accumbens has been proposed to participate in the modulation of drug-seeking behavior, the effects of glutamate receptor blockade in this brain structure on cue-induced relapse were investigated. The studies employed animals models in which rats were trained to press a lever in a test cage to obtain alcohol or intravenous cocaine. Drug availability was paired with distinct olfactory, auditory, or visual stimuli. This phase was followed by extinction training, during which lever presses did not result in the presentation of the drug or the drug-associated stimuli. Extinction training led to a gradual decrease in the number of lever presses during test sessions. Relapse was triggered by presenting the rats with the drug-associated stimuli in the absence of alcohol or cocaine. The drug-associated stimuli were alone capable of inducing resumption of lever pressing and maintaining this behavior during repeated testing. The number of lever presses during a session represented the intensity of drug-seeking and relapse behavior. The results suggest that glutamatergic neurotransmission is involved in the modulation of drug-seeking behavior. Both alcohol and cocaine relapse were attenuated by systemic pretreatment with glutamate receptor antagonists. However, differences were found in the ability of ionotropic AMPA/kainate and NMDA receptor antagonists to regulate drug-seeking behavior. The AMPA/kainate antagonists CNQX and NBQX, and L-701,324, an antagonist with affinity for the glycine site of the NMDA receptor, attenuated cue-induced drug seeking, whereas the competitive NMDA antagonist CGP39551 and the NMDA channel blocker MK-801 were without effect. MPEP, an antagonist at metabotropic mGlu5 glutamate receptors, also decreased drug seeking, but its administration was found to lead to conditioned suppression of behavior during subsequent treatment sessions, suggesting that MPEP may have undesirable side effects. The mGluR2/3 agonist LY379268 and the mGluR8 agonist (S)-3,4-DCPG decreased both cue-induced relapse to alcohol drinking and alcohol consumption. Control experiments showed however that administration of the agonists was accompanied by motor suppression limiting their usefulness. Administration of the AMPA/kainate antagonist CNQX, the NMDA antagonist D-AP5, and the mGluR5 antagonist MPEP into the nucleus accumbens resulted also in a decrease in drug-seeking behavior, suggesting that the nucleus accumbens is at least one of the anatomical sites regulating drug seeking and mediating the effects of glutamate receptor antagonists on this behavior.