Common genetic variants influence human subcortical brain structures

Multidimensional heritability analysis of neuroanatomical shape

In the dawning era of large-scale biomedical data, multidimensional phenotype vectors will play an increasing role in examining the genetic underpinnings of brain features, behaviour and disease. For example, shape measurements derived from brain MRI scans are multidimensional geometric descriptions of brain structure and provide an alternate class of phenotypes that remains largely unexplored in genetic studies. Here we extend the concept of heritability to multidimensional traits, and present the first comprehensive analysis of the heritability of neuroanatomical shape measurements across an ensemble of brain structures based on genome-wide SNP and MRI data from 1,320 unrelated, young and healthy individuals. We replicate our findings in an extended twin sample from the Human Connectome Project (HCP). Our results demonstrate that neuroanatomical shape can be significantly heritable, above and beyond volume, and can serve as a complementary phenotype to study the genetic determinants and clinical relevance of brain structure.

HASE: Framework for efficient high-dimensional association analyses

High-throughput technology can now provide rich information on a person’s biological makeup and environmental surroundings. Important discoveries have been made by relating these data to various health outcomes in fields such as genomics, proteomics, and medical imaging. However, cross-investigations between several high-throughput technologies remain impractical due to demanding computational requirements (hundreds of years of computing resources) and unsuitability for collaborative settings (terabytes of data to share). Here we introduce the HASE framework that overcomes both of these issues. Our approach dramatically reduces computational time from years to only hours and also requires several gigabytes to be exchanged between collaborators. We implemented a novel meta-analytical method that yields identical power as pooled analyses without the need of sharing individual participant data. The efficiency of the framework is illustrated by associating 9 million genetic variants with 1.5 million brain imaging voxels in three cohorts (total N = 4,034) followed by meta-analysis, on a standard computational infrastructure. These experiments indicate that HASE facilitates high-dimensional association studies enabling large multicenter association studies for future discoveries.

Association between polygenic risk for schizophrenia, neurocognition and social cognition across development

Breakthroughs in genomics have begun to unravel the genetic architecture of schizophrenia risk, providing methods for quantifying schizophrenia polygenic risk based on common genetic variants. Our objective in the current study was to understand the relationship between schizophrenia genetic risk variants and neurocognitive development in healthy individuals. We first used combined genomic and neurocognitive data from the Philadelphia Neurodevelopmental Cohort (4303 participants ages 8-21 years) to screen 26 neurocognitive phenotypes for their association with schizophrenia polygenic risk. Schizophrenia polygenic risk was estimated for each participant based on summary statistics from the most recent schizophrenia genome-wide association analysis (Psychiatric Genomics Consortium 2014). After correction for multiple comparisons, greater schizophrenia polygenic risk was significantly associated with reduced speed of emotion identification and verbal reasoning. These associations were significant by age 9 years and there was no evidence of interaction between schizophrenia polygenic risk and age on neurocognitive performance. We then looked at the association between schizophrenia polygenic risk and emotion identification speed in the Harvard/MGH Brain Genomics Superstruct Project sample (695 participants ages 18-35 years), where we replicated the association between schizophrenia polygenic risk and emotion identification speed. These analyses provide evidence for a replicable association between polygenic risk for schizophrenia and a specific aspect of social cognition. Our findings indicate that individual differences in genetic risk for schizophrenia are linked with the development of aspects of social cognition and potentially verbal reasoning, and that these associations emerge relatively early in development.

Abnormal asymmetries in subcortical brain volume in schizophrenia

Subcortical structures, which include the basal ganglia and parts of the limbic system, have key roles in learning, motor control and emotion, but also contribute to higher-order executive functions. Prior studies have reported volumetric alterations in subcortical regions in schizophrenia. Reported results have sometimes been heterogeneous, and few large-scale investigations have been conducted. Moreover, few large-scale studies have assessed asymmetries of subcortical volumes in schizophrenia. Here, as a work completely independent of a study performed by the ENIGMA consortium, we conducted a large-scale multisite study of subcortical volumetric differences between patients with schizophrenia and controls. We also explored the laterality of subcortical regions to identify characteristic similarities and differences between them. T1-weighted images from 1680 healthy individuals and 884 patients with schizophrenia, obtained with 15 imaging protocols at 11 sites, were processed with FreeSurfer. Group differences were calculated for each protocol and meta-analyzed. Compared with controls, patients with schizophrenia demonstrated smaller bilateral hippocampus, amygdala, thalamus and accumbens volumes as well as intracranial volume, but larger bilateral caudate, putamen, pallidum and lateral ventricle volumes. We replicated the rank order of effect sizes for subcortical volumetric changes in schizophrenia reported by the ENIGMA consortium. Further, we revealed leftward asymmetry for thalamus, lateral ventricle, caudate and putamen volumes, and rightward asymmetry for amygdala and hippocampal volumes in both controls and patients with schizophrenia. Also, we demonstrated a schizophrenia-specific leftward asymmetry for pallidum volume. These findings suggest the possibility of aberrant laterality in neural pathways and connectivity patterns related to the pallidum in schizophrenia.

The domesticated brain: genetics of brain mass and brain structure in an avian species

As brain size usually increases with body size it has been assumed that the two are tightly constrained and evolutionary studies have therefore often been based on relative brain size (i.e. brain size proportional to body size) rather than absolute brain size. The process of domestication offers an excellent opportunity to disentangle the linkage between body and brain mass due to the extreme selection for increased body mass that has occurred. By breeding an intercross between domestic chicken and their wild progenitor, we address this relationship by simultaneously mapping the genes that control inter-population variation in brain mass and body mass. Loci controlling variation in brain mass and body mass have separate genetic architectures and are therefore not directly constrained. Genetic mapping of brain regions indicates that domestication has led to a larger body mass and to a lesser extent a larger absolute brain mass in chickens, mainly due to enlargement of the cerebellum. Domestication has traditionally been linked to brain mass regression, based on measurements of relative brain mass, which confounds the large body mass augmentation due to domestication. Our results refute this concept in the chicken.

Alterations in the hippocampus and thalamus in individuals at high risk for psychosis

Reduction in hippocampal volume is a hallmark of schizophrenia and already present in the clinical high-risk state. Nevertheless, other subcortical structures, such as the thalamus, amygdala and pallidum can differentiate schizophrenia patients from controls. We studied the role of hippocampal and subcortical structures in clinical high-risk individuals from two cohorts. High-resolution T₁-weighted structural MRI brain scans of a total of 91 clinical high-risk individuals and 64 healthy controls were collected in two centers. The bilateral volume of the hippocampus, the thalamus, the caudate, the putamen, the pallidum, the amygdala, and the accumbens were automatically segmented using FSL-FIRST. A linear mixed-effects model and a prospective meta-analysis were applied to assess group-related volumetric differences. We report reduced hippocampal and thalamic volumes in clinical high-risk individuals compared to healthy controls. No volumetric alterations were detected for the caudate, the putamen, the pallidum, the amygdala, or the accumbens. Moreover, we found comparable medium effect sizes for group-related comparison of the thalamus in the two analytical methods. These findings underline the relevance of specific alterations in the hippocampal and subcortical volumes in the high-risk state. Further analyses may allow hippocampal and thalamic volumes to be used as biomarkers to predict psychosis.

A Population-Based Study of Four Genes Associated with Heroin Addiction in Han Chinese

Recent studies have shown that variants in FAT atypical cadherin 3 (FAT3), kinectin 1 (KTN1), discs large homolog2 (DLG2) and deleted in colorectal cancer (DCC) genes influence the structure of the human mesolimbic reward system. We conducted a systematic analysis of the potential functional single nucleotide polymorphisms (SNPs) in these genes associated with heroin addiction. We scanned the functional regions of these genes and identified 20 SNPs for genotyping by using the SNaPshot method. A total of 1080 samples, comprising 523 cases and 557 controls, were analyzed. We observed that DCC rs16956878, rs12607853, and rs2292043 were associated with heroin addiction. The T alleles of rs16956878 (p = 0.0004) and rs12607853 (p = 0.002) were significantly enriched in the case group compared with the controls. A lower incidence of the C allele of rs2292043 (p = 0.002) was observed in the case group. In block 2 of DCC (rs2292043-rs12607853-rs16956878), the frequency of the T-T-T haplotype was significantly higher in the case group than in the control group (p = 0.024), and fewer C-C-C haplotypes (p = 0.006) were detected in the case group. DCC may be an important candidate gene in heroin addiction, and rs16956878, rs12607853, and rs2292043 may be risk factors, thereby providing a basis for further genetic and biological research.

Generalized reduced rank latent factor regression for high dimensional tensor fields, and neuroimaging-genetic applications

We propose a generalized reduced rank latent factor regression model (GRRLF) for the analysis of tensor field responses and high dimensional covariates. The model is motivated by the need from imaging-genetic studies to identify genetic variants that are associated with brain imaging phenotypes, often in the form of high dimensional tensor fields. GRRLF identifies from the structure in the data the effective dimensionality of the data, and then jointly performs dimension reduction of the covariates, dynamic identification of latent factors, and nonparametric estimation of both covariate and latent response fields. After accounting for the latent and covariate effects, GRLLF performs a nonparametric test on the remaining factor of interest. GRRLF provides a better factorization of the signals compared with common solutions, and is less susceptible to overfitting because it exploits the effective dimensionality. The generality and the flexibility of GRRLF also allow various statistical models to be handled in a unified framework and solutions can be efficiently computed. Within the field of neuroimaging, it improves the sensitivity for weak signals and is a promising alternative to existing approaches. The operation of the framework is demonstrated with both synthetic datasets and a real-world neuroimaging example in which the effects of a set of genes on the structure of the brain at the voxel level were measured, and the results compared favorably with those from existing approaches.

Population neuroscience

Population neuroscience endeavors to identify influences shaping the human brain from conception onwards, thus generating knowledge relevant for building and maintaining brain health throughout the life span. This can be achieved by studying large samples of participants drawn from the general population and evaluated with state-of-the-art tools for assessing (a) genes and their regulation; (b) external and internal environments; and (c) brain properties. This chapter reviews the three elements of population neuroscience (principles, tools, innovations, limitations), and discusses future directions in this field.

Quantitative magnetic resonance imaging traits as endophenotypes for genetic mapping in epilepsy

Over the last decade, the field of imaging genomics has combined high-throughput genotype data with quantitative magnetic resonance imaging (QMRI) measures to identify genes associated with brain structure, cognition, and several brain-related disorders. Despite its successful application in different psychiatric and neurological disorders, the field has yet to be advanced in epilepsy. In this article we examine the relevance of imaging genomics for future genetic studies in epilepsy from three perspectives. First, we discuss prior genome-wide genetic mapping efforts in epilepsy, considering the possibility that some studies may have been constrained by inherent theoretical and methodological limitations of the genome-wide association study (GWAS) method. Second, we offer a brief overview of the imaging genomics paradigm, from its original inception, to its role in the discovery of important risk genes in a number of brain-related disorders, and its successful application in large-scale multinational research networks. Third, we provide a comprehensive review of past studies that have explored the eligibility of brain QMRI traits as endophenotypes for epilepsy. While the breadth of studies exploring QMRI-derived endophenotypes in epilepsy remains narrow, robust syndrome-specific neuroanatomical QMRI traits have the potential to serve as accessible and relevant intermediate phenotypes for future genetic mapping efforts in epilepsy.

•

QMRI traits have the potential to serve as robust intermediate phenotypes for brain-related disorders.

•

Hippocampal volume is the most promising neuroimaging endophenotype for MTLE + HS.

•

Imaging genomics holds great promise in advancing epilepsy genetic research.

•

Studies are encouraged to explore the validity of QMRI traits as endophenotypes for epilepsy.

Distinct Genetic Influences on Cortical and Subcortical Brain Structures

This study examined the heritability of brain grey matter structures in a subsample of older adult twins (93 MZ and 68 DZ twin pairs; mean age 70 years) from the Older Australian Twins Study. The heritability estimates of subcortical regions ranged from 0.41 (amygdala) to 0.73 (hippocampus), and of cortical regions, from 0.55 (parietal lobe) to 0.78 (frontal lobe). Corresponding structures in the two hemispheres were influenced by the same genetic factors and high genetic correlations were observed between the two hemispheric regions. There were three genetically correlated clusters, comprising (i) the cortical lobes (frontal, temporal, parietal and occipital lobes); (ii) the basal ganglia (caudate, putamen and pallidum) with weak genetic correlations with cortical lobes, and (iii) the amygdala, hippocampus, thalamus and nucleus accumbens grouped together, which genetically correlated with both basal ganglia and cortical lobes, albeit relatively weakly. Our study demonstrates a complex but patterned and clustered genetic architecture of the human brain, with divergent genetic determinants of cortical and subcortical structures, in particular the basal ganglia.

Genetic Variation in Schizophrenia Liability is Shared With Intellectual Ability and Brain Structure

BACKGROUND

Alterations in intellectual ability and brain structure are important genetic markers for schizophrenia liability. How variations in these phenotypes interact with variance in schizophrenia liability due to genetic or environmental factors is an area of active investigation. Studying these genetic markers using a multivariate twin modeling approach can provide novel leads for (genetic) pathways of schizophrenia development.

METHODS

In a sample of 70 twins discordant for schizophrenia and 130 healthy control twins, structural equation modeling was applied to quantify unique contributions of genetic and environmental factors on human brain structure (cortical thickness, cortical surface and global white matter fractional anisotropy [FA]), intellectual ability and schizophrenia liability.

RESULTS

In total, up to 28.1% of the genetic variance (22.8% of total variance) in schizophrenia liability was shared with intelligence quotient (IQ), global-FA, cortical thickness, and cortical surface. The strongest contributor was IQ, sharing on average 16.4% of the genetic variance in schizophrenia liability, followed by cortical thickness (6.3%), global-FA (4.7%) and cortical surface (0.5%). Furthermore, we found that up to 57.4% of the variation due to environmental factors (4.6% of total variance) in schizophrenia was shared with IQ (34.2%) and cortical surface (13.4%).

CONCLUSIONS

Intellectual ability, FA and cortical thickness show significant and independent shared genetic variance with schizophrenia liability. This suggests that measuring brain-imaging phenotypes helps explain genetic variance in schizophrenia liability that is not captured by variation in IQ.

Analysis of 23andMe antidepressant efficacy survey data: implication of circadian rhythm and neuroplasticity in bupropion response

Genetic predisposition may contribute to the differences in drug-specific, class-specific or antidepressant-wide treatment resistance. Clinical studies with the genetic data are often limited in sample sizes. Drug response obtained from self-reports may offer an alternative approach to conduct a study with much larger sample size. Using the phenotype data collected from 23andMe 'Antidepressant Efficacy and Side Effects' survey and genotype data from 23andMe's research participants, we conducted genome-wide association study (GWAS) on subjects of European ancestry using four groups of phenotypes (a) non-treatment-resistant depression (n=7795) vs treatment-resistant depression (TRD, n=1311), (b) selective serotonin reuptake inhibitors (SSRI) responders (n=6348) vs non-responders (n=3340), (c) citalopram/escitalopram responders (n=2963) vs non-responders (n=2005), and (d) norepinephrine-dopamine reuptake inhibitor (NDRI, bupropion) responders (n=2675) vs non-responders (n=1861). Each of these subgroups was also compared with controls (n ~ 190 000). The most significant association was from bupropion responders vs non-responders analysis. Variant rs1908557 (P=2.6 × 10(-8), OR=1.35) passed the conventional genome-wide significance threshold (P=5 × 10(-8)) and was located within the intron of human spliced expressed sequence tags in chromosome 4. Gene sets associated with long-term depression, circadian rhythm and vascular endothelial growth factor (VEGF) pathway were enriched in the bupropion analysis. No single-nucleotide polymorphism passed genome-wide significance threshold in other analyses. The heritability estimates for each response group compared with controls were between 0.15 and 0.25, consistent with the known heritability for major depressive disorder.

Phenome-wide analysis of genome-wide polygenic scores

Genome-wide polygenic scores (GPS), which aggregate the effects of thousands of DNA variants from genome-wide association studies (GWAS), have the potential to make genetic predictions for individuals. We conducted a systematic investigation of associations between GPS and many behavioral traits, the behavioral phenome. For 3152 unrelated 16-year-old individuals representative of the United Kingdom, we created 13 GPS from the largest GWAS for psychiatric disorders (for example, schizophrenia, depression and dementia) and cognitive traits (for example, intelligence, educational attainment and intracranial volume). The behavioral phenome included 50 traits from the domains of psychopathology, personality, cognitive abilities and educational achievement. We examined phenome-wide profiles of associations for the entire distribution of each GPS and for the extremes of the GPS distributions. The cognitive GPS yielded stronger predictive power than the psychiatric GPS in our UK-representative sample of adolescents. For example, education GPS explained variation in adolescents' behavior problems (~0.6%) and in educational achievement (~2%) but psychiatric GPS were associated with neither. Despite the modest effect sizes of current GPS, quantile analyses illustrate the ability to stratify individuals by GPS and opportunities for research. For example, the highest and lowest septiles for the education GPS yielded a 0.5 s.d. difference in mean math grade and a 0.25 s.d. difference in mean behavior problems. We discuss the usefulness and limitations of GPS based on adult GWAS to predict genetic propensities earlier in development.

Cross-Phenotype Polygenic Risk Score Analysis of Persistent Post-Concussive Symptoms in U.S. Army Soldiers with Deployment-Acquired Traumatic Brain Injury

Traumatic brain injury (TBI) contributes to the increased rates of suicide and post-traumatic stress disorder in military personnel and veterans, and it is also associated with the risk for neurodegenerative and psychiatric disorders. A cross-phenotype high-resolution polygenic risk score (PRS) analysis of persistent post-concussive symptoms (PCS) was conducted in 845 U.S. Army soldiers who sustained TBI during their deployment. We used a prospective longitudinal survey of three brigade combat teams to assess deployment-acquired TBI and persistent physical, cognitive, and emotional PCS. PRS was derived from summary statistics of large genome-wide association studies of Alzheimer's disease, Parkinson's disease, schizophrenia, bipolar disorder, and major depressive disorder (MDD); and for years of schooling, college completion, childhood intelligence, infant head circumference (IHC), and adult intracranial volume. Although our study had more than 95% of statistical power to detect moderate-to-large effect sizes, no association was observed with neurodegenerative and psychiatric disorders, suggesting that persistent PCS does not share genetic components with these traits to a moderate-to-large degree. We observed a significant finding: subjects with high IHC PRS recovered better from cognitive/emotional persistent PCS than the other individuals (R² = 1.11%; p = 3.37 × 10^-3). Enrichment analysis identified two significant Gene Ontology (GO) terms related to this result: GO:0050839∼Cell adhesion molecule binding (p = 8.9 × 10^-6) and GO:0050905∼Neuromuscular process (p = 9.8 × 10^-5). In summary, our study indicated that the genetic predisposition to persistent PCS after TBI does not have substantial overlap with neurodegenerative and psychiatric diseases, but mechanisms related to early brain growth may be involved.

COINSTAC: A Privacy Enabled Model and Prototype for Leveraging and Processing Decentralized Brain Imaging Data

The field of neuroimaging has embraced the need for sharing and collaboration. Data sharing mandates from public funding agencies and major journal publishers have spurred the development of data repositories and neuroinformatics consortia. However, efficient and effective data sharing still faces several hurdles. For example, open data sharing is on the rise but is not suitable for sensitive data that are not easily shared, such as genetics. Current approaches can be cumbersome (such as negotiating multiple data sharing agreements). There are also significant data transfer, organization and computational challenges. Centralized repositories only partially address the issues. We propose a dynamic, decentralized platform for large scale analyses called the Collaborative Informatics and Neuroimaging Suite Toolkit for Anonymous Computation (COINSTAC). The COINSTAC solution can include data missing from central repositories, allows pooling of both open and "closed" repositories by developing privacy-preserving versions of widely-used algorithms, and incorporates the tools within an easy-to-use platform enabling distributed computation. We present an initial prototype system which we demonstrate on two multi-site data sets, without aggregating the data. In addition, by iterating across sites, the COINSTAC model enables meta-analytic solutions to converge to "pooled-data" solutions (i.e., as if the entire data were in hand). More advanced approaches such as feature generation, matrix factorization models, and preprocessing can be incorporated into such a model. In sum, COINSTAC enables access to the many currently unavailable data sets, a user friendly privacy enabled interface for decentralized analysis, and a powerful solution that complements existing data sharing solutions.

A study of genetic and environmental contributions to structural brain changes over time in twins concordant and discordant for bipolar disorder

This is the first longitudinal twin study examining genetic and environmental contributions to the association between liability to bipolar disorder (BD) and changes over time in global brain volumes, and global and regional measures of cortical surface area, cortical thickness and cortical volume. A total of 50 twins from pairs discordant or concordant for BD (monozygotic: 8 discordant and 3 concordant pairs, and 1 patient and 3 co-twins from incomplete pairs; dizygotic: 6 discordant and 2 concordant pairs, and 1 patient and 7 co-twins from incomplete pairs) underwent magnetic resonance imaging twice. In addition, 57 twins from healthy twin pairs (15 monozygotic and 10 dizygotic pairs, and 4 monozygotic and 3 dizygotic subjects from incomplete pairs) were also scanned twice. Mean follow-up duration for all twins was 7.5 years (standard deviation: 1.5 years). Data were analyzed using structural equation modeling software OpenMx. The liability to BD was not associated with global or regional structural brain changes over time. Although we observed a subtle increase in cerebral white matter in BD patients, this effect disappeared after correction for multiple comparisons. Heritability of brain changes over time was generally low to moderate. Structural brain changes appear to follow similar trajectories in BD patients and healthy controls. Existing brain abnormalities in BD do not appear to progressively change over time, but this requires additional confirmation. Further study with large cohorts is recommended to assess genetic and environmental influences on structural brain abnormalities in BD, while taking into account the influence of lithium on the brain.

Impact of polygenic schizophrenia-related risk and hippocampal volumes on the onset of psychosis

Alterations in hippocampal volume are a known marker for first-episode psychosis (FEP) as well as for the clinical high-risk state. The Polygenic Schizophrenia-related Risk Score (PSRS), derived from a large case-control study, indicates the polygenic predisposition for schizophrenia in our clinical sample. A total of 65 at-risk mental state (ARMS) and FEP patients underwent structural magnetic resonance imaging. We used automatic segmentation of hippocampal volumes using the FSL-FIRST software and an odds-ratio-weighted PSRS based on the publicly available top single-nucleotide polymorphisms from the Psychiatric Genomics Consortium genome-wide association study (GWAS). We observed a negative association between the PSRS and hippocampal volumes (β=-0.42, P=0.01, 95% confidence interval (CI)=(-0.72 to -0.12)) across FEP and ARMS patients. Moreover, a higher PSRS was significantly associated with a higher probability of an individual being assigned to the FEP group relative to the ARMS group (β=0.64, P=0.03, 95% CI=(0.08-1.29)). These findings provide evidence that a subset of schizophrenia risk variants is negatively associated with hippocampal volumes, and higher values of this PSRS are significantly associated with FEP compared with the ARMS. This implies that FEP patients have a higher genetic risk for schizophrenia than the total cohort of ARMS patients. The identification of associations between genetic risk variants and structural brain alterations will increase our understanding of the neurobiology underlying the transition to psychosis.

Conservation of Distinct Genetically-Mediated Human Cortical Pattern

The many subcomponents of the human cortex are known to follow an anatomical pattern and functional relationship that appears to be highly conserved between individuals. This suggests that this pattern and the relationship among cortical regions are important for cortical function and likely shaped by genetic factors, although the degree to which genetic factors contribute to this pattern is unknown. We assessed the genetic relationships among 12 cortical surface areas using brain images and genotype information on 2,364 unrelated individuals, brain images on 466 twin pairs, and transcriptome data on 6 postmortem brains in order to determine whether a consistent and biologically meaningful pattern could be identified from these very different data sets. We find that the patterns revealed by each data set are highly consistent (p<10⁻³), and are biologically meaningful on several fronts. For example, close genetic relationships are seen in cortical regions within the same lobes and, the frontal lobe, a region showing great evolutionary expansion and functional complexity, has the most distant genetic relationship with other lobes. The frontal lobe also exhibits the most distinct expression pattern relative to the other regions, implicating a number of genes with known functions mediating immune and related processes. Our analyses reflect one of the first attempts to provide an assessment of the biological consistency of a genetic phenomenon involving the brain that leverages very different types of data, and therefore is not just statistical replication which purposefully use very similar data sets.

Although functional and anatomical connections among cortical regions have been intensively explored, genetically-mediated relationships between cortical regions have not been pursued to the same degree. Identifying genetic factors that mediate these relationships among different brain subcomponents can provide insight into how the human brain is organized and functions. We have assessed the genetic relationships among cortical regions using an integrated approach that considers twin data, genotype information among a large set of unrelated individuals, and gene expression measurements from postmortem neural tissues. We looked for evidence that subsets of cortical brain regions are under common or unique genetic control. We found that the patterns of genetic relationships are highly consistent across three independent data sets and multiple lines of evidence, suggesting that the patterning of cortical surface area is strongly mediated by genetic factors and, furthermore, likely reflects underlying anatomical and possibly functional relationships among cortical brain regions.

Polygenic risk of Alzheimer disease is associated with early- and late-life processes

OBJECTIVE

To examine associations between aggregate genetic risk and Alzheimer disease (AD) markers in stages preceding the clinical symptoms of dementia using data from 2 large observational cohort studies.

METHODS

We computed polygenic risk scores (PGRS) using summary statistics from the International Genomics of Alzheimer's Project genome-wide association study of AD. Associations between PGRS and AD markers (cognitive decline, clinical progression, hippocampus volume, and β-amyloid) were assessed within older participants with dementia. Associations between PGRS and hippocampus volume were additionally examined within healthy younger participants (age 18-35 years).

RESULTS

Within participants without dementia, elevated PGRS was associated with worse memory (p = 0.002) and smaller hippocampus (p = 0.002) at baseline, as well as greater longitudinal cognitive decline (memory: p = 0.0005, executive function: p = 0.01) and clinical progression (p < 0.00001). High PGRS was associated with AD-like levels of β-amyloid burden as measured with florbetapir PET (p = 0.03) but did not reach statistical significance for CSF β-amyloid (p = 0.11). Within the younger group, higher PGRS was associated with smaller hippocampus volume (p = 0.05). This pattern was evident when examining a PGRS that included many loci below the genome-wide association study (GWAS)-level significance threshold (16,123 single nucleotide polymorphisms), but not when PGRS was restricted to GWAS-level significant loci (18 single nucleotide polymorphisms).

CONCLUSIONS

Effects related to common genetic risk loci distributed throughout the genome are detectable among individuals without dementia. The influence of this genetic risk may begin in early life and make an individual more susceptible to cognitive impairment in late life. Future refinement of polygenic risk scores may help identify individuals at risk for AD dementia.

Limited Evidence for Association of Genome-Wide Schizophrenia Risk Variants on Cortical Neuroimaging Phenotypes

BACKGROUND

There are now over 100 established genetic risk variants for schizophrenia; however, their influence on brain structure and circuitry across the human lifespan are not known.

METHODS

We examined healthy individuals 8-86 years of age, from the Centre for Addiction and Mental Health, the Zucker Hillside Hospital, and the Philadelphia Neurodevelopmental Cohort. Following thorough quality control procedures, we investigated associations of established genetic risk variants with heritable neuroimaging phenotypes relevant to schizophrenia, namely thickness of frontal and temporal cortical regions (n = 565) and frontotemporal and interhemispheric white matter tract fractional anisotropy (FA) (n = 530).

RESULTS

There was little evidence for association of risk variants with imaging phenotypes. No association with cortical thickness of any region was present. Only rs12148337, near a long noncoding RNA region, was associated with white matter FA (splenium of corpus callosum) following multiple comparison correction (corrected p = .012); this single nucleotide polymorphism was also associated with genu FA and superior longitudinal fasciculus FA at p <.005 (uncorrected). There was no association of polygenic risk score with white matter FA or cortical thickness.

CONCLUSIONS

In sum, our findings provide limited evidence for association of schizophrenia risk variants with cortical thickness or diffusion imaging white matter phenotypes. When taken with recent lack of association of these variants with subcortical brain volumes, our results either suggest that structural neuroimaging approaches at current resolution are not sufficiently sensitive to detect effects of these risk variants or that multiple comparison correction in correlated phenotypes is too stringent, potentially "eliminating" biologically important signals.

Additive sex-specific influence of common non-synonymous DISC1 variants on amygdala, basal ganglia, and white cortical surface area in healthy young adults

The disrupted-in-schizophrenia-1 (DISC1) gene is known for its role in the development of mental disorders. It is also involved in neurodevelopment, cognition, and memory. To investigate the association between DISC1 variants and brain morphology, we analyzed the influence of the three common non-synonymous polymorphisms in DISC1 on specific brain structures in healthy young adults. The volumes of brain regions were determined in 145 subjects by magnetic resonance imaging and automated analysis using FreeSurfer. Genotyping was performed by high resolution melting of amplified products. In an additive genetic model, rs6675281 (Leu607Phe), rs3738401 (Arg264Gln), and rs821616 (Ser704Cys) significantly explained the volume variance of the amygdala (p = 0.007) and the pallidum (p = 0.004). A higher cumulative portion of minor alleles was associated with larger volumes of the amygdala (p = 0.005), the pallidum (p = 0.001), the caudate (p = 0.024), and the putamen (p = 0.007). Sex-stratified analysis revealed a strong genetic effect of rs6675281 on putamen and pallidum in females but not in males and an opposite influence of rs3738401 on the white cortical surface in females compared to males. The strongest single association was found for rs821616 and the amygdala volume in male subjects (p < 0.001). No effect was detected for the nucleus accumbens. We report-to our knowledge-for the first time a significant and sex-specific influence of common DISC1 variants on volumes of the basal ganglia, the amygdala and on the cortical surface area. Our results demonstrate that the additive model of all three polymorphisms outperforms their single analysis.

Genetic variants in Alzheimer disease - molecular and brain network approaches

Genetic studies in late-onset Alzheimer disease (LOAD) are aimed at identifying core disease mechanisms and providing potential biomarkers and drug candidates to improve clinical care of AD. However, owing to the complexity of LOAD, including pathological heterogeneity and disease polygenicity, extraction of actionable guidance from LOAD genetics has been challenging. Past attempts to summarize the effects of LOAD-associated genetic variants have used pathway analysis and collections of small-scale experiments to hypothesize functional convergence across several variants. In this Review, we discuss how the study of molecular, cellular and brain networks provides additional information on the effects of LOAD-associated genetic variants. We then discuss emerging combinations of these omic data sets into multiscale models, which provide a more comprehensive representation of the effects of LOAD-associated genetic variants at multiple biophysical scales. Furthermore, we highlight the clinical potential of mechanistically coupling genetic variants and disease phenotypes with multiscale brain models.

Impact of PICALM and CLU on hippocampal degeneration

PICALM and CLU are two major risk genes of late-onset Alzheimer's disease (LOAD), and there is strong molecular evidence suggesting their interaction on amyloid-beta deposition, hence finding functional dependency between their risk genotypes may lead to better understanding of their roles in LOAD development and greater clinical utility. In this study, we mainly investigated interaction effects of risk loci PICALM rs3581179 and CLU rs11136000 on hippocampal degeneration in both young and elderly adults in order to understand their neural mechanism on aging process, which may help identify robust biomarkers for early diagnosis and intervention. Besides volume we also assessed hippocampal shape phenotypes derived from diffeomorphic metric mapping and nonlinear dimensionality reduction. In elderly individuals (75.6 ± 6.7 years) significant interaction effects existed on hippocampal volume (P < 0.001), whereas in young healthy adults (19.4 ± 1.1 years) such effects existed on a shape phenotype (P = 0.01) indicating significant variation at hippocampal head and tail that mirror most AD vulnerable regions. Voxel-wise analysis also pointed to the same regions but lacked statistical power. In both cohorts, PICALM protective genotype AA only exhibited protective effects on hippocampal degeneration and cognitive performance when combined with CLU protective T allele, but adverse effects with CLU risk CC. This study revealed novel PICALM and CLU interaction effects on hippocampal degeneration along aging, and validated effectiveness of diffeomorphometry in imaging genetics study. Hum Brain Mapp 37:2419-2430, 2016. © 2016 Wiley Periodicals, Inc.

Gene-set and multivariate genome-wide association analysis of oppositional defiant behavior subtypes in attention-deficit/hyperactivity disorder

Oppositional defiant disorder (ODD) is a frequent psychiatric disorder seen in children and adolescents with attention-deficit-hyperactivity disorder (ADHD). ODD is also a common antecedent to both affective disorders and aggressive behaviors. Although the heritability of ODD has been estimated to be around 0.60, there has been little research into the molecular genetics of ODD. The present study examined the association of irritable and defiant/vindictive dimensions and categorical subtypes of ODD (based on latent class analyses) with previously described specific polymorphisms (DRD4 exon3 VNTR, 5-HTTLPR, and seven OXTR SNPs) as well as with dopamine, serotonin, and oxytocin genes and pathways in a clinical sample of children and adolescents with ADHD. In addition, we performed a multivariate genome-wide association study (GWAS) of the aforementioned ODD dimensions and subtypes. Apart from adjusting the analyses for age and sex, we controlled for "parental ability to cope with disruptive behavior." None of the hypothesis-driven analyses revealed a significant association with ODD dimensions and subtypes. Inadequate parenting behavior was significantly associated with all ODD dimensions and subtypes, most strongly with defiant/vindictive behaviors. In addition, the GWAS did not result in genome-wide significant findings but bioinformatics and literature analyses revealed that the proteins encoded by 28 of the 53 top-ranked genes functionally interact in a molecular landscape centered around Beta-catenin signaling and involved in the regulation of neurite outgrowth. Our findings provide new insights into the molecular basis of ODD and inform future genetic studies of oppositional behavior. © 2015 The Authors. American Journal of Medical Genetics Part B: Neuropsychiatric Genetics Published by Wiley Periodicals, Inc.

Heritability of cortical thickness changes over time in twin pairs discordant for schizophrenia

BACKGROUND

Cortical thickness and surface area changes have repeatedly been found in schizophrenia. Whether progressive loss in cortical thickness and surface area are mediated by genetic or disease related factors is unknown. Here we investigate to what extent genetic and/or environmental factors contribute to the association between change in cortical thickness and surface area and liability to develop schizophrenia.

METHOD

Longitudinal magnetic resonance imaging study over a 5-year interval. Monozygotic (MZ) and dizygotic (DZ) twin pairs discordant for schizophrenia were compared with healthy control twin pairs using repeated measures analysis of variance (RM-ANOVA) and structural equation modeling (SEM). Twins discordant for schizophrenia and healthy control twins were recruited from the twin cohort at the University Medical Centre Utrecht, The Netherlands. A total of 90 individuals from 46 same sex twin pairs were included: 9 MZ and 10 DZ discordant for schizophrenia and 14 MZ and 13 (11 complete and 2 incomplete) DZ healthy twin-pairs. Age varied between 19 and 57years.

RESULTS

Higher genetic liability for schizophrenia was associated with progressive global thinning of the cortex, particularly of the left superior temporal cortex. Higher environmental liability for schizophrenia was associated with global attenuated thinning of the cortex, and including of the left superior temporal cortex. Cortical surface area change was heritable, but not significantly associated with higher genetic or environmental liability for schizophrenia.

CONCLUSIONS

Excessive cortical thinning, particularly of the left superior temporal cortex, may represent a genetic risk marker for schizophrenia.

Genetic and Environmental Factors Associated with Cannabis Involvement

Approximately 50-70% of the variation in cannabis use and use disorders can be attributed to heritable factors. For cannabis use, the remaining variance can be parsed in to familial and person-specific environmental factors while for use disorders, only the latter contribute. While numerous candidate gene studies have identified the role of common variation influencing liability to cannabis involvement, replication has been elusive. To date, no genomewide association study has been sufficiently powered to identify significant loci. Despite this, studies adopting polygenic techniques and integrating genetic variation with neural phenotypes and measures of environmental risk, such as childhood adversity, are providing promising new leads. It is likely that the small effect sizes associated with variants related to cannabis involvement will only be robustly identified in substantially larger samples. Results of such large-scale efforts will provide valuable single variant targets for translational research in neurogenetic, pharmacogenetic and non-human animal models as well as polygenic risk indices that can be used to explore a host of other genetic hypotheses related to cannabis use and misuse.

Genome-wide association study of cognitive functions and educational attainment in UK Biobank (N=112 151)

People's differences in cognitive functions are partly heritable and are associated with important life outcomes. Previous genome-wide association (GWA) studies of cognitive functions have found evidence for polygenic effects yet, to date, there are few replicated genetic associations. Here we use data from the UK Biobank sample to investigate the genetic contributions to variation in tests of three cognitive functions and in educational attainment. GWA analyses were performed for verbal-numerical reasoning (N=36 035), memory (N=112 067), reaction time (N=111 483) and for the attainment of a college or a university degree (N=111 114). We report genome-wide significant single-nucleotide polymorphism (SNP)-based associations in 20 genomic regions, and significant gene-based findings in 46 regions. These include findings in the ATXN2, CYP2DG, APBA1 and CADM2 genes. We report replication of these hits in published GWA studies of cognitive function, educational attainment and childhood intelligence. There is also replication, in UK Biobank, of SNP hits reported previously in GWA studies of educational attainment and cognitive function. GCTA-GREML analyses, using common SNPs (minor allele frequency>0.01), indicated significant SNP-based heritabilities of 31% (s.e.m.=1.8%) for verbal-numerical reasoning, 5% (s.e.m.=0.6%) for memory, 11% (s.e.m.=0.6%) for reaction time and 21% (s.e.m.=0.6%) for educational attainment. Polygenic score analyses indicate that up to 5% of the variance in cognitive test scores can be predicted in an independent cohort. The genomic regions identified include several novel loci, some of which have been associated with intracranial volume, neurodegeneration, Alzheimer's disease and schizophrenia.

Genome research in pre-dementia stages of Alzheimer's disease

Genetic characterization of individuals at risk of Alzheimer's disease (AD), i.e. people having amyloid deposits in the brain without symptoms, people suffering from subjective cognitive decline (SCD) or mild cognitive impairment (MCI), has spurred the interests of researchers. However, their pre-dementia genetic profile remains mostly unexplored. In this study, we reviewed the loci related to phenotypes of AD, MCI and SCD from literature and performed the first meta-analyses evaluating the role of apolipoprotein E (APOE) in the risk of conversion from a healthy status to MCI and SCD. For AD dementia risk, an increased number of loci have been identified; to date, 28 genes have been associated with Late Onset AD. In MCI syndrome, APOE is confirmed as a pheno-conversion factor leading from MCI to AD, and clusterin is a promising candidate. Additionally, our meta-analyses revealed APOE as genetic risk factor to convert from a healthy status to MCI [OR = 1.849 (1.587-2.153); P = 2.80  × 10-15] and to a lesser extent from healthy status to SCD [OR = 1.151 (1.015-1.304); P = 0.028]. Thus, we believe that genetic studies in longitudinal SCD and MCI series may provide new therapeutic targets and improve the existing knowledge of AD. This type of studies must be completed on healthy subjects to better understand the natural disease resistance to brain insults and neurodegeneration.

Longitudinal Evidence for Smaller Hippocampus Volume as a Vulnerability Factor for Perceived Stress

Hippocampal volume has been found to be smaller in individuals with stress-related disorders, but it remains unclear whether smaller volume is a consequence of stress or rather a vulnerability factor. Here, we examined this issue by relating stress levels to hippocampal volumes in healthy participants examined every 5 years in a longitudinal population-based study. Based on scores of 25- to 60-year–old participants on the perceived stress questionnaire, we defined moderately to high (n = 35) and low (n = 76) stress groups. The groups were re-examined after 5 years (at the 6th study wave). Historical data on subjective stress were available up to 10 years prior to Wave 5. At the first MRI session, the moderately to high stress group had a significantly smaller hippocampal volume, as measured by FreeSurfer (version 5.3), compared with the low-stress group. At follow-up, group differences in stress levels and hippocampal volume remained unchanged. In retrospective analyses of subjective stress, the observed group difference in stress was found to be stable. The long-term stability of group differences in perceived stress and hippocampal volume suggests that a small hippocampal volume may be a vulnerability factor for stress-related disorders.

Developmentally Sensitive Interaction Effects of Genes and the Social Environment on Total and Subcortical Brain Volumes

Smaller total brain and subcortical volumes have been linked to psychopathology including attention-deficit/hyperactivity disorder (ADHD). Identifying mechanisms underlying these alterations, therefore, is of great importance. We investigated the role of gene-environment interactions (GxE) in interindividual variability of total gray matter (GM), caudate, and putamen volumes. Brain volumes were derived from structural magnetic resonance imaging scans in participants with (N = 312) and without ADHD (N = 437) from N = 402 families (age M = 17.00, SD = 3.60). GxE effects between DAT1, 5-HTT, and DRD4 and social environments (maternal expressed warmth and criticism; positive and deviant peer affiliation) as well as the possible moderating effect of age were examined using linear mixed modeling. We also tested whether findings depended on ADHD severity. Deviant peer affiliation was associated with lower caudate volume. Participants with low deviant peer affiliations had larger total GM volumes with increasing age. Likewise, developmentally sensitive GxE effects were found on total GM and putamen volume. For total GM, differential age effects were found for DAT1 9-repeat and HTTLPR L/L genotypes, depending on the amount of positive peer affiliation. For putamen volume, DRD4 7-repeat carriers and DAT1 10/10 homozygotes showed opposite age relations depending on positive peer affiliation and maternal criticism, respectively. All results were independent of ADHD severity. The presence of differential age-dependent GxE effects might explain the diverse and sometimes opposing results of environmental and genetic effects on brain volumes observed so far.

Gene networks and developmental context: the importance of understanding complex gene expression patterns in evolution

Animal development is the product of distinct components and interactions-genes, regulatory networks, and cells-and it exhibits emergent properties that cannot be inferred from the components in isolation. Often the focus is on the genotype-to-phenotype map, overlooking the process of development that turns one into the other. We propose a move toward micro-evolutionary analysis of development, incorporating new tools that enable cell type resolution and single-cell microscopy. Using the sex determination pathway in Drosophila to illustrate potential avenues of research, we highlight some of the questions that these emerging technologies can address. For example, they provide an unprecedented opportunity to study heterogeneity within cell populations, and the potential to add the dimension of time to gene regulatory network analysis. Challenges still remain in developing methods to analyze this data and to increase the throughput. However this line of research has the potential to bridge the gaps between previously more disparate fields, such as population genetics and development, opening up new avenues of research.

Corticolimbic anatomical characteristics predetermine risk for chronic pain

SEE TRACEY DOI101093/BRAIN/AWW147 FOR A SCIENTIFIC COMMENTARY ON THIS ARTICLE: Mechanisms of chronic pain remain poorly understood. We tracked brain properties in subacute back pain patients longitudinally for 3 years as they either recovered from or transitioned to chronic pain. Whole-brain comparisons indicated corticolimbic, but not pain-related circuitry, white matter connections predisposed patients to chronic pain. Intra-corticolimbic white matter connectivity analysis identified three segregated communities: dorsal medial prefrontal cortex-amygdala-accumbens, ventral medial prefrontal cortex-amygdala, and orbitofrontal cortex-amygdala-hippocampus. Higher incidence of white matter and functional connections within the dorsal medial prefrontal cortex-amygdala-accumbens circuit, as well as smaller amygdala volume, represented independent risk factors, together accounting for 60% of the variance for pain persistence. Opioid gene polymorphisms and negative mood contributed indirectly through corticolimbic anatomical factors, to risk for chronic pain. Our results imply that persistence of chronic pain is predetermined by corticolimbic neuroanatomical factors.

Structural Brain Abnormalities in Youth With Psychosis Spectrum Symptoms

IMPORTANCE

Structural brain abnormalities are prominent in psychotic disorders, including schizophrenia. However, it is unclear when aberrations emerge in the disease process and if such deficits are present in association with less severe psychosis spectrum (PS) symptoms in youth.

OBJECTIVE

To investigate the presence of structural brain abnormalities in youth with PS symptoms.

DESIGN, SETTING, AND PARTICIPANTS

The Philadelphia Neurodevelopmental Cohort is a prospectively accrued, community-based sample of 9498 youth who received a structured psychiatric evaluation. A subsample of 1601 individuals underwent neuroimaging, including structural magnetic resonance imaging, at an academic and children's hospital health care network between November 1, 2009, and November 30, 2011.

MAIN OUTCOMES AND MEASURES

Measures of brain volume derived from T1-weighted structural neuroimaging at 3 T. Analyses were conducted at global, regional, and voxelwise levels. Regional volumes were estimated with an advanced multiatlas regional segmentation procedure, and voxelwise volumetric analyses were conducted as well. Nonlinear developmental patterns were examined using penalized splines within a general additive model. Psychosis spectrum (PS) symptom severity was summarized using factor analysis and evaluated dimensionally.

RESULTS

Following exclusions due to comorbidity and image quality assurance, the final sample included 791 participants aged youth 8 to 22 years. Fifty percent (n = 393) were female. After structured interviews, 391 participants were identified as having PS features (PS group) and 400 participants were identified as typically developing comparison individuals without significant psychopathology (TD group). Compared with the TD group, the PS group had diminished whole-brain gray matter volume (P = 1.8 × 10-10) and expanded white matter volume (P = 2.8 × 10-11). Voxelwise analyses revealed significantly lower gray matter volume in the medial temporal lobe (maximum z score = 5.2 and cluster size of 1225 for the right and maximum z score = 4.5 and cluster size of 310 for the left) as well as in frontal, temporal, and parietal cortex. Volumetric reduction in the medial temporal lobe was correlated with PS symptom severity.

CONCLUSIONS AND RELEVANCE

Structural brain abnormalities that have been commonly reported in adults with psychosis are present early in life in youth with PS symptoms and are not due to medication effects. Future longitudinal studies could use the presence of such abnormalities in conjunction with clinical presentation, cognitive profile, and genomics to predict risk and aid in stratification to guide early interventions.

The CHRM3 gene is implicated in abnormal thalamo-orbital frontal cortex functional connectivity in first-episode treatment-naive patients with schizophrenia

BACKGROUND

The genetic influences in human brain structure and function and impaired functional connectivities are the hallmarks of the schizophrenic brain. To explore how common genetic variants affect the connectivities in schizophrenia, we applied genome-wide association studies assaying the abnormal neural connectivities in schizophrenia as quantitative traits.

METHOD

We recruited 161 first-onset and treatment-naive patients with schizophrenia and 150 healthy controls. All the participants underwent scanning with a 3 T-magnetic resonance imaging scanner to acquire structural and functional imaging data and genotyping using the HumanOmniZhongHua-8 BeadChip. The brain-wide association study approach was employed to account for the inherent modular nature of brain connectivities.

RESULTS

We found differences in four abnormal functional connectivities [left rectus to left thalamus (REC.L-THA.L), left rectus to right thalamus (REC.L-THA.R), left superior orbital cortex to left thalamus (ORBsup.L-THA.L) and left superior orbital cortex to right thalamus (ORBsup.L-THA.R)] between the two groups. Univariate single nucleotide polymorphism (SNP)-based association revealed that the SNP rs6800381, located nearest to the CHRM3 (cholinergic receptor, muscarinic 3) gene, reached genomic significance (p = 1.768 × 10-8) using REC.L-THA.R as the phenotype. Multivariate gene-based association revealed that the FAM12A (family with sequence similarity 12, member A) gene nearly reached genomic significance (nominal p = 2.22 × 10-6, corrected p = 0.05).

CONCLUSIONS

Overall, we identified the first evidence that the CHRM3 gene plays a role in abnormal thalamo-orbital frontal cortex functional connectivity in first-episode treatment-naive patients with schizophrenia. Identification of these genetic variants using neuroimaging genetics provides insights into the causes of variability in human brain development, and may help us determine the mechanisms of dysfunction in schizophrenia.

Magnetic resonance imaging in Alzheimer's Disease Neuroimaging Initiative 2

INTRODUCTION

Alzheimer's Disease Neuroimaging Initiative (ADNI) is now in its 10th year. The primary objective of the magnetic resonance imaging (MRI) core of ADNI has been to improve methods for clinical trials in Alzheimer's disease (AD) and related disorders.

METHODS

We review the contributions of the MRI core from present and past cycles of ADNI (ADNI-1, -Grand Opportunity and -2). We also review plans for the future-ADNI-3.

RESULTS

Contributions of the MRI core include creating standardized acquisition protocols and quality control methods; examining the effect of technical features of image acquisition and analysis on outcome metrics; deriving sample size estimates for future trials based on those outcomes; and piloting the potential utility of MR perfusion, diffusion, and functional connectivity measures in multicenter clinical trials.

DISCUSSION

Over the past decade the MRI core of ADNI has fulfilled its mandate of improving methods for clinical trials in AD and will continue to do so in the future.

An Association Study Between Genetic Polymorphisms in Functional Regions of Five Genes and the Risk of Schizophrenia

Schizophrenia is a severe mental disorder that is likely to be strongly determined by genetic factors. To identify markers of disks, large homolog 2 (DLG2), FAT atypical cadherin 3 (FAT3), kinectin1 (KTN1), deleted in colorectal carcinoma (DCC), and glycogen synthase kinase-3β (GSK3β) that contribute to the genetic susceptibility to schizophrenia, we systematically screened for polymorphisms in the functional regions of these genes. A total of 22 functional single-nucleotide polymorphisms (SNPs) in 940 Chinese subjects were genotyped using SNaPshot. The results first suggested that the allelic and genotypic frequencies of the DCC polymorphism rs2229080 were nominally associated with schizophrenia. The patients were significantly less likely to be CC homozygous (P = 0.005, odds ratio [OR] = 0.635, 95 % confidence interval [95 % CI] = 0.462-0.873), and the schizophrenia subjects exhibited lower frequency of the C allele (P = 0.024, OR = 0.811, 95 % CI = 0.676-0.972). Regarding GSK3β, there was a significant difference in genotype distribution of rs3755557 between schizophrenia and healthy control subjects (P = 0.009). The patients exhibited a significantly lower frequency of the T allele of rs3755557 (P = 0.002, OR = 0.654, 95 % CI = 0.498-0.860). Our results point to the polymorphisms of DCC and GSK3β as contributors to the genetic basis of individual differences in the susceptibility to schizophrenia.

Gene Prioritization for Imaging Genetics Studies Using Gene Ontology and a Stratified False Discovery Rate Approach

Imaging genetics is an emerging field in which the association between genes and neuroimaging-based quantitative phenotypes are used to explore the functional role of genes in neuroanatomy and neurophysiology in the context of healthy function and neuropsychiatric disorders. The main obstacle for researchers in the field is the high dimensionality of the data in both the imaging phenotypes and the genetic variants commonly typed. In this article, we develop a novel method that utilizes Gene Ontology, an online database, to select and prioritize certain genes, employing a stratified false discovery rate (sFDR) approach to investigate their associations with imaging phenotypes. sFDR has the potential to increase power in genome wide association studies (GWAS), and is quickly gaining traction as a method for multiple testing correction. Our novel approach addresses both the pressing need in genetic research to move beyond candidate gene studies, while not being overburdened with a loss of power due to multiple testing. As an example of our methodology, we perform a GWAS of hippocampal volume using both the Enhancing NeuroImaging Genetics through Meta-Analysis (ENIGMA2) and the Alzheimer's Disease Neuroimaging Initiative datasets. The analysis of ENIGMA2 data yielded a set of SNPs with sFDR values between 10 and 20%. Our approach demonstrates a potential method to prioritize genes based on biological systems impaired in a disease.

Toward an integrative science of the developing human mind and brain: Focus on the developing cortex

Based on the Huttenlocher lecture, this article describes the need for a more integrative scientific paradigm for addressing important questions raised by key observations made over 2 decades ago. Among these are the early descriptions by Huttenlocher of variability in synaptic density in cortex of postmortem brains of children of different ages and the almost simultaneous reports of cortical volume reductions on MR imaging in children and adolescents. In spite of much progress in developmental neurobiology, developmental cognitive neuroscience, and behavioral and imaging genetics, we still do not know how these early observations relate to each other. It is argued that large scale, collaborative research programs are needed to establish the associations between behavioral differences among children and imaging biomarkers, and to link the latter to cellular changes in the developing brain. Examples of progress and challenges remaining are illustrated with data from the Pediatric Imaging, Neurocognition, and Genetics Project (PING).

Sex-specific mediation effect of the right fusiform face area volume on the association between variants in repeat length of AVPR1A RS3 and altruistic behavior in healthy adults

Microsatellite variants in the arginine vasopressin receptor 1A gene (AVPR1A) RS3 have been associated with normal social behaviors variation and autism spectrum disorders (ASDs) in a sex-specific manner. However, neural mechanisms underlying these associations remain largely unknown. We hypothesized that AVPR1A RS3 variants affect altruistic behavior by modulating the gray matter volume (GMV) of specific brain regions in a sex-specific manner. We investigated 278 young healthy adults using the Dictator Game to assess altruistic behavior. All subjects were genotyped and main effect of AVPR1A RS3 repeat polymorphisms and interaction of genotype-by-sex on the GMV were assessed in a voxel-wise manner. We observed that male subjects with relatively short repeats allocated less money to others and exhibited a significantly smaller GMV in the right fusiform face area (FFA) compared with male long homozygotes. In male subjects, the GMV of the right FFA exhibited a significant positive correlation with altruistic behavior. A mixed mediation and moderation analysis further revealed both a significant mediation effect of the GMV of the right FFA on the association between AVPR1A RS3 repeat polymorphisms and allocation sums and a significant moderation effect of sex (only in males) on the mediation effect. Post hoc analysis showed that the GMV of the right FFA was significantly smaller in male subjects carrying allele 426 than in non-426 carriers. These results suggest that the GMV of the right FFA may be a potential mediator whereby the genetic variants in AVPR1A RS3 affect altruistic behavior in healthy male subjects. Hum Brain Mapp 37:2700-2709, 2016. © 2016 Wiley Periodicals, Inc.

"RADIOTRANSCRIPTOMICS": A synergy of imaging and transcriptomics in clinical assessment

Recent advances in quantitative imaging and "omics" technology have generated a wealth of mineable biological "big data". With the push towards a P4 "predictive, preventive, personalized, and participatory" approach to medicine, researchers began integrating complementary tools to further tune existing diagnostic and therapeutic models. The field of radiogenomics has long pioneered such multidisciplinary investigations in neuroscience and oncology, correlating genotypic and phenotypic signatures to study structural and functional changes in relation to altered molecular behavior. Given the innate dynamic nature of complex disorders and the role of environmental and epigenetic factors in pathogenesis, the transcriptome can further elucidate serial modifications undetected at the genome level. We therefore propose "radiotranscriptomics" as a new member of the P4 medicine initiative, combining transcriptome information, including gene expression and isoform variation, and quantitative image annotations.

Genetic influences on schizophrenia and subcortical brain volumes: large-scale proof of concept

Schizophrenia is a devastating psychiatric illness with high heritability. Brain structure and function differ, on average, between people with schizophrenia and healthy individuals. As common genetic associations are emerging for both schizophrenia and brain imaging phenotypes, we can now use genome-wide data to investigate genetic overlap. Here we integrated results from common variant studies of schizophrenia (33,636 cases, 43,008 controls) and volumes of several (mainly subcortical) brain structures (11,840 subjects). We did not find evidence of genetic overlap between schizophrenia risk and subcortical volume measures either at the level of common variant genetic architecture or for single genetic markers. These results provide a proof of concept (albeit based on a limited set of structural brain measures) and define a roadmap for future studies investigating the genetic covariance between structural or functional brain phenotypes and risk for psychiatric disorders.

Heritability and reliability of automatically segmented human hippocampal formation subregions

The human hippocampal formation can be divided into a set of cytoarchitecturally and functionally distinct subregions, involved in different aspects of memory formation. Neuroanatomical disruptions within these subregions are associated with several debilitating brain disorders including Alzheimer's disease, major depression, schizophrenia, and bipolar disorder. Multi-center brain imaging consortia, such as the Enhancing Neuro Imaging Genetics through Meta-Analysis (ENIGMA) consortium, are interested in studying disease effects on these subregions, and in the genetic factors that affect them. For large-scale studies, automated extraction and subsequent genomic association studies of these hippocampal subregion measures may provide additional insight. Here, we evaluated the test-retest reliability and transplatform reliability (1.5T versus 3T) of the subregion segmentation module in the FreeSurfer software package using three independent cohorts of healthy adults, one young (Queensland Twins Imaging Study, N=39), another elderly (Alzheimer's Disease Neuroimaging Initiative, ADNI-2, N=163) and another mixed cohort of healthy and depressed participants (Max Planck Institute, MPIP, N=598). We also investigated agreement between the most recent version of this algorithm (v6.0) and an older version (v5.3), again using the ADNI-2 and MPIP cohorts in addition to a sample from the Netherlands Study for Depression and Anxiety (NESDA) (N=221). Finally, we estimated the heritability (h(2)) of the segmented subregion volumes using the full sample of young, healthy QTIM twins (N=728). Test-retest reliability was high for all twelve subregions in the 3T ADNI-2 sample (intraclass correlation coefficient (ICC)=0.70-0.97) and moderate-to-high in the 4T QTIM sample (ICC=0.5-0.89). Transplatform reliability was strong for eleven of the twelve subregions (ICC=0.66-0.96); however, the hippocampal fissure was not consistently reconstructed across 1.5T and 3T field strengths (ICC=0.47-0.57). Between-version agreement was moderate for the hippocampal tail, subiculum and presubiculum (ICC=0.78-0.84; Dice Similarity Coefficient (DSC)=0.55-0.70), and poor for all other subregions (ICC=0.34-0.81; DSC=0.28-0.51). All hippocampal subregion volumes were highly heritable (h(2)=0.67-0.91). Our findings indicate that eleven of the twelve human hippocampal subregions segmented using FreeSurfer version 6.0 may serve as reliable and informative quantitative phenotypes for future multi-site imaging genetics initiatives such as those of the ENIGMA consortium.

Genetics of Bipolar Disorder: Recent Update and Future Directions

Although genetic studies of Bipolar Disorder have been pursued for decades, it has only been in the last several years that clearly replicated findings have emerged. These findings, typically of modest effects, point to a polygenic genetic architecture consisting of multiple common and rare susceptibility variants. While larger genome-wide association studies are ongoing, the advent of whole exome and genome sequencing should lead to the identification of rare, and potentially more penetrant, variants. Progress along both fronts will provide novel insights into the biology of Bipolar Disorder and help usher in a new era of personalized medicine and improved treatments.

Annual Research Review: Discovery science strategies in studies of the pathophysiology of child and adolescent psychiatric disorders--promises and limitations

BACKGROUND

Psychiatric science remains descriptive, with a categorical nosology intended to enhance interobserver reliability. Increased awareness of the mismatch between categorical classifications and the complexity of biological systems drives the search for novel frameworks including discovery science in Big Data. In this review, we provide an overview of incipient approaches, primarily focused on classically categorical diagnoses such as schizophrenia (SZ), autism spectrum disorder (ASD), and attention-deficit/hyperactivity disorder (ADHD), but also reference convincing, if focal, advances in cancer biology, to describe the challenges of Big Data and discovery science, and outline approaches being formulated to overcome existing obstacles.

FINDINGS

A paradigm shift from categorical diagnoses to a domain/structure-based nosology and from linear causal chains to complex causal network models of brain-behavior relationship is ongoing. This (r)evolution involves appreciating the complexity, dimensionality, and heterogeneity of neuropsychiatric data collected from multiple sources ('broad' data) along with data obtained at multiple levels of analysis, ranging from genes to molecules, cells, circuits, and behaviors ('deep' data). Both of these types of Big Data landscapes require the use and development of robust and powerful informatics and statistical approaches. Thus, we describe Big Data analysis pipelines and the promise and potential limitations in using Big Data approaches to study psychiatric disorders.

CONCLUSIONS

We highlight key resources available for psychopathological studies and call for the application and development of Big Data approaches to dissect the causes and mechanisms of neuropsychiatric disorders and identify corresponding biomarkers for early diagnosis.

Modulative effects of COMT haplotype on age-related associations with brain morphology

Catechol-O-methyltransferase (COMT), located on chromosome 22q11.2, encodes an enzyme critical for dopamine flux in the prefrontal cortex. Genetic variants of COMT have been suggested to functionally manipulate prefrontal morphology and function in healthy adults. This study aims to investigate modulative roles of individuals COMT SNPs (rs737865, val158met, rs165599) and its haplotypes in age-related brain morphology using an Asian sample with 174 adults aged from 21 to 80 years. We showed an age-related decline in cortical thickness of the dorsal visual pathway, including the left dorsolateral prefrontal cortex, bilateral angular gyrus, right superior frontal cortex, and age-related shape compression in the basal ganglia as a function of the genotypes of the individual COMT SNPs, especially COMT val158met. Using haplotype trend regression analysis, COMT haplotype probabilities were estimated and further revealed an age-related decline in cortical thickness in the default mode network (DMN), including the posterior cingulate, precuneus, supramarginal and paracentral cortex, and the ventral visual system, including the occipital cortex and left inferior temporal cortex, as a function of the COMT haplotype. Our results provided new evidence on an antagonistic pleiotropic effect in COMT, suggesting that genetically programmed neural benefits in early life may have a potential bearing towards neural susceptibility in later life. Hum Brain Mapp 37:2068-2082, 2016. © 2016 Wiley Periodicals, Inc.

Early developmental gene enhancers affect subcortical volumes in the adult human brain

Genome-wide association screens aim to identify common genetic variants contributing to the phenotypic variability of complex traits, such as human height or brain morphology. The identified genetic variants are mostly within noncoding genomic regions and the biology of the genotype-phenotype association typically remains unclear. In this article, we propose a complementary targeted strategy to reveal the genetic underpinnings of variability in subcortical brain volumes, by specifically selecting genomic loci that are experimentally validated forebrain enhancers, active in early embryonic development. We hypothesized that genetic variation within these enhancers may affect the development and ultimately the structure of subcortical brain regions in adults. We tested whether variants in forebrain enhancer regions showed an overall enrichment of association with volumetric variation in subcortical structures of >13,000 healthy adults. We observed significant enrichment of genomic loci that affect the volume of the hippocampus within forebrain enhancers (empirical P = 0.0015), a finding which robustly passed the adjusted threshold for testing of multiple brain phenotypes (cutoff of P < 0.0083 at an alpha of 0.05). In analyses of individual single nucleotide polymorphisms (SNPs), we identified an association upstream of the ID2 gene with rs7588305 and variation in hippocampal volume. This SNP-based association survived multiple-testing correction for the number of SNPs analyzed but not for the number of subcortical structures. Targeting known regulatory regions offers a way to understand the underlying biology that connects genotypes to phenotypes, particularly in the context of neuroimaging genetics. This biology-driven approach generates testable hypotheses regarding the functional biology of identified associations. Hum Brain Mapp 37:1788-1800, 2016. © 2016 Wiley Periodicals, Inc.

Estimating Effect Sizes and Expected Replication Probabilities from GWAS Summary Statistics

Genome-wide Association Studies (GWAS) result in millions of summary statistics (“z-scores”) for single nucleotide polymorphism (SNP) associations with phenotypes. These rich datasets afford deep insights into the nature and extent of genetic contributions to complex phenotypes such as psychiatric disorders, which are understood to have substantial genetic components that arise from very large numbers of SNPs. The complexity of the datasets, however, poses a significant challenge to maximizing their utility. This is reflected in a need for better understanding the landscape of z-scores, as such knowledge would enhance causal SNP and gene discovery, help elucidate mechanistic pathways, and inform future study design. Here we present a parsimonious methodology for modeling effect sizes and replication probabilities, relying only on summary statistics from GWAS substudies, and a scheme allowing for direct empirical validation. We show that modeling z-scores as a mixture of Gaussians is conceptually appropriate, in particular taking into account ubiquitous non-null effects that are likely in the datasets due to weak linkage disequilibrium with causal SNPs. The four-parameter model allows for estimating the degree of polygenicity of the phenotype and predicting the proportion of chip heritability explainable by genome-wide significant SNPs in future studies with larger sample sizes. We apply the model to recent GWAS of schizophrenia (N = 82,315) and putamen volume (N = 12,596), with approximately 9.3 million SNP z-scores in both cases. We show that, over a broad range of z-scores and sample sizes, the model accurately predicts expectation estimates of true effect sizes and replication probabilities in multistage GWAS designs. We assess the degree to which effect sizes are over-estimated when based on linear-regression association coefficients. We estimate the polygenicity of schizophrenia to be 0.037 and the putamen to be 0.001, while the respective sample sizes required to approach fully explaining the chip heritability are 10⁶ and 10⁵. The model can be extended to incorporate prior knowledge such as pleiotropy and SNP annotation. The current findings suggest that the model is applicable to a broad array of complex phenotypes and will enhance understanding of their genetic architectures.

Neurobiological origin of spurious brain morphological changes: A quantitative MRI study

The high gray‐white matter contrast and spatial resolution provided by T1‐weighted magnetic resonance imaging (MRI) has made it a widely used imaging protocol for computational anatomy studies of the brain. While the image intensity in T1‐weighted images is predominantly driven by T1, other MRI parameters affect the image contrast, and hence brain morphological measures derived from the data. Because MRI parameters are correlates of different histological properties of brain tissue, this mixed contribution hampers the neurobiological interpretation of morphometry findings, an issue which remains largely ignored in the community. We acquired quantitative maps of the MRI parameters that determine signal intensities in T1‐weighted images (R₁ (=1/T1), R₂*, and PD) in a large cohort of healthy subjects (n = 120, aged 18–87 years). Synthetic T1‐weighted images were calculated from these quantitative maps and used to extract morphometry features—gray matter volume and cortical thickness. We observed significant variations in morphometry measures obtained from synthetic images derived from different subsets of MRI parameters. We also detected a modulation of these variations by age. Our findings highlight the impact of microstructural properties of brain tissue—myelination, iron, and water content—on automated measures of brain morphology and show that microstructural tissue changes might lead to the detection of spurious morphological changes in computational anatomy studies. They motivate a review of previous morphological results obtained from standard anatomical MRI images and highlight the value of quantitative MRI data for the inference of microscopic tissue changes in the healthy and diseased brain. Hum Brain Mapp 37:1801–1815, 2016. © 2016 The Authors. Human Brain Mapping Published by Wiley Periodicals, Inc.