Massively expedited genome-wide heritability analysis (MEGHA)

Diverging functional connectivity timescales: Capturing distinct aspects of cognitive performance in early psychosis

BACKGROUND

Psychosis spectrum disorders (PSDs) are marked by cognitive impairments, the neurobiological correlates of which remain poorly understood. Here, we investigate the entropy of time-varying functional connectivity (TVFC) patterns from resting-state functional magnetic resonance imaging (rs-fMRI) as potential biomarker for cognitive performance in PSDs. By combining our results with multimodal reference data, we hope to generate new insights into the mechanisms underlying cognitive dysfunction in PSDs. We hypothesized that low-entropy TVFC patterns (LEN) would be more behaviorally informative than high-entropy TVFC patterns (HEN), especially for tasks that require extensive integration across diverse cognitive subdomains.

METHODS

rs-fMRI and behavioral data from 97 patients in the early phases of psychosis and 53 controls were analyzed. Positron emission tomography (PET) and magnetoencephalography (MEG) data were taken from a public repository (Hansen et al., 2022). Multivariate analyses were conducted to examine relationships between TVFC patterns at multiple spatial scales and cognitive performance in patients.

RESULTS

Compared to HEN, LEN explained significantly more cognitive variance on average in PSD patients, driven by superior encoding of information on psychometrically more integrated tasks. HEN better captured information in specific subdomains of executive functioning. Nodal HEN-LEN transitions were spatially aligned with neurobiological gradients reflecting monoaminergic transporter densities and MEG beta-power. Exploratory analyses revealed a close statistical relationship between LEN and positive symptom severity in patients.

CONCLUSION

Our entropy-based analysis of TVFC patterns dissociates distinct aspects of cognition in PSDs. By linking topographies of neurotransmission and oscillatory dynamics with cognitive performance, it enhances our understanding of the mechanisms underlying cognitive deficits in PSDs.

Employing Informatics Strategies in Alzheimer's Disease Research: A Review from Genetics, Multiomics, and Biomarkers to Clinical Outcomes

Alzheimer's disease (AD) is a critical national concern, affecting 5.8 million people and costing more than $250 billion annually. However, there is no available cure. Thus, effective strategies are in urgent need to discover AD biomarkers for disease early detection and drug development. In this review, we study AD from a biomedical data scientist perspective to discuss the four fundamental components in AD research: genetics (G), molecular multiomics (M), multimodal imaging biomarkers (B), and clinical outcomes (O) (collectively referred to as the GMBO framework). We provide a comprehensive review of common statistical and informatics methodologies for each component within the GMBO framework, accompanied by the major findings from landmark AD studies. Our review highlights the potential of multimodal biobank data in addressing key challenges in AD, such as early diagnosis, disease heterogeneity, and therapeutic development. We identify major hurdles in AD research, including data scarcity and complexity, and advocate for enhanced collaboration, data harmonization, and advanced modeling techniques. This review aims to be an essential guide for understanding current biomedical data science strategies in AD research, emphasizing the need for integrated, multidisciplinary approaches to advance our understanding and management of AD.

Diverging functional connectivity timescales: Capturing distinct aspects of cognitive performance in early psychosis

Background

Psychosis spectrum disorders (PSDs) are marked by cognitive impairments, the neurobiological correlates of which remain poorly understood. Here, we investigate the entropy of time-varying functional connectivity (TVFC) patterns from resting-state fMRI (rfMRI) as potential biomarker for cognitive performance in PSDs. By combining our results with multimodal reference data, we hope to generate new insights into the mechanisms underlying cognitive dysfunction in PSDs. We hypothesized that low-entropy TVFC patterns (LEN) would be more behaviorally informative than high-entropy TVFC patterns (HEN), especially for tasks that require extensive integration across diverse cognitive subdomains.

Methods

rfMRI and behavioral data from 97 patients in the early phases of psychosis and 53 controls were analyzed. Positron-Emission Tomography (PET) and magnetoencephalography (MEG) data were taken from a public repository (Hansen et al., 2022). Multivariate analyses were conducted to examine relationships between TVFC patterns at multiple spatial scales and cognitive performance in patients.

Results

Compared to HEN, LEN explained significantly more cognitive variance on average in PSD patients, driven by superior encoding of information on psychometrically more integrated tasks. HEN better captured information in specific subdomains of executive functioning. Nodal HEN-LEN transitions were spatially aligned with neurobiological gradients reflecting monoaminergic transporter densities and MEG beta power. Exploratory analyses revealed a close statistical relationship between LEN and positive PSD symptoms.

Conclusion

Our entropy-based analysis of TVFC patterns dissociates distinct aspects of cognition in PSDs. By linking topographies of neurotransmission and oscillatory dynamics with cognitive performance, it enhances our understanding of the mechanisms underlying cognitive deficits in PSDs.

CRediT Authorship Contribution Statement

Fabian Hirsch: Conceptualization, Methodology, Software, Formal analysis, Writing - Original Draft, Writing - Review & Editing, Visualization; Ângelo Bumanglag: Methodology, Software, Formal analysis, Writing - Review & Editing; Yifei Zhang: Methodology, Software, Formal analysis, Writing - Review & Editing; Afra Wohlschlaeger: Methodology, Writing - Review & Editing, Supervision, Project administration.

SAN: Mitigating spatial covariance heterogeneity in cortical thickness data collected from multiple scanners or sites

In neuroimaging studies, combining data collected from multiple study sites or scanners is becoming common to increase the reproducibility of scientific discoveries. At the same time, unwanted variations arise by using different scanners (inter-scanner biases), which need to be corrected before downstream analyses to facilitate replicable research and prevent spurious findings. While statistical harmonization methods such as ComBat have become popular in mitigating inter-scanner biases in neuroimaging, recent methodological advances have shown that harmonizing heterogeneous covariances results in higher data quality. In vertex-level cortical thickness data, heterogeneity in spatial autocorrelation is a critical factor that affects covariance heterogeneity. Our work proposes a new statistical harmonization method called spatial autocorrelation normalization (SAN) that preserves homogeneous covariance vertex-level cortical thickness data across different scanners. We use an explicit Gaussian process to characterize scanner-invariant and scanner-specific variations to reconstruct spatially homogeneous data across scanners. SAN is computationally feasible, and it easily allows the integration of existing harmonization methods. We demonstrate the utility of the proposed method using cortical thickness data from the Social Processes Initiative in the Neurobiology of the Schizophrenia(s) (SPINS) study. SAN is publicly available as an R package.

MIR137 polygenic risk for schizophrenia and ephrin-regulated pathway: Role in lateral ventricles and corpus callosum volume

Background/Objective. Enlarged lateral ventricle (LV) volume and decreased volume in the corpus callosum (CC) are hallmarks of schizophrenia (SZ). We previously showed an inverse correlation between LV and CC volumes in SZ, with global functioning decreasing with increased LV volume. This study investigates the relationship between LV volume, CC abnormalities, and the microRNA MIR137 and its regulated genes in SZ, because of MIR137's essential role in neurodevelopment. Methods. Participants were 1224 SZ probands and 1466 unaffected controls from the GENUS Consortium. Brain MRI scans, genotype, and clinical data were harmonized across cohorts and employed in the analyses. Results. Increased LV volumes and decreased CC central, mid-anterior, and mid-posterior volumes were observed in SZ probands. The MIR137-regulated ephrin pathway was significantly associated with CC:LV ratio, explaining a significant proportion (3.42 %) of CC:LV variance, and more than for LV and CC separately. Other pathways explained variance in either CC or LV, but not both. CC:LV ratio was also positively correlated with Global Assessment of Functioning, supporting previous subsample findings. SNP-based heritability estimates were higher for CC central:LV ratio (0.79) compared to CC or LV separately. Discussion. Our results indicate that the CC:LV ratio is highly heritable, influenced in part by variation in the MIR137-regulated ephrin pathway. Findings suggest that the CC:LV ratio may be a risk indicator in SZ that correlates with global functioning.

SAN: mitigating spatial covariance heterogeneity in cortical thickness data collected from multiple scanners or sites

In neuroimaging studies, combining data collected from multiple study sites or scanners is becoming common to increase the reproducibility of scientific discoveries. At the same time, unwanted variations arise by using different scanners (inter-scanner biases), which need to be corrected before downstream analyses to facilitate replicable research and prevent spurious findings. While statistical harmonization methods such as ComBat have become popular in mitigating inter-scanner biases in neuroimaging, recent methodological advances have shown that harmonizing heterogeneous covariances results in higher data quality. In vertex-level cortical thickness data, heterogeneity in spatial autocorrelation is a critical factor that affects covariance heterogeneity. Our work proposes a new statistical harmonization method called SAN (Spatial Autocorrelation Normalization) that preserves homogeneous covariance vertex-level cortical thickness data across different scanners. We use an explicit Gaussian process to characterize scanner-invariant and scanner-specific variations to reconstruct spatially homogeneous data across scanners. SAN is computationally feasible, and it easily allows the integration of existing harmonization methods. We demonstrate the utility of the proposed method using cortical thickness data from the Social Processes Initiative in the Neurobiology of the Schizophrenia(s) (SPINS) study. SAN is publicly available as an R package.

Static and dynamic fMRI-derived functional connectomes represent largely similar information

Functional connectivity (FC) of blood oxygen level-dependent (BOLD) fMRI time series can be estimated using methods that differ in sensitivity to the temporal order of time points (static vs. dynamic) and the number of regions considered in estimating a single edge (bivariate vs. multivariate). Previous research suggests that dynamic FC explains variability in FC fluctuations and behavior beyond static FC. Our aim was to systematically compare methods on both dimensions. We compared five FC methods: Pearson's/full correlation (static, bivariate), lagged correlation (dynamic, bivariate), partial correlation (static, multivariate), and multivariate AR model with and without self-connections (dynamic, multivariate). We compared these methods by (i) assessing similarities between FC matrices, (ii) by comparing node centrality measures, and (iii) by comparing the patterns of brain-behavior associations. Although FC estimates did not differ as a function of sensitivity to temporal order, we observed differences between the multivariate and bivariate FC methods. The dynamic FC estimates were highly correlated with the static FC estimates, especially when comparing group-level FC matrices. Similarly, there were high correlations between the patterns of brain-behavior associations obtained using the dynamic and static FC methods. We conclude that the dynamic FC estimates represent information largely similar to that of the static FC.

Spatial-extent inference for testing variance components in reliability and heritability studies

Clusterwise inference is a popular approach in neuroimaging to increase sensitivity, but most existing methods are currently restricted to the General Linear Model (GLM) for testing mean parameters. Statistical methods for testing variance components, which are critical in neuroimaging studies that involve estimation of narrow-sense heritability or test-retest reliability, are underdeveloped due to methodological and computational challenges, which would potentially lead to low power. We propose a fast and powerful test for variance components called CLEAN-V (CLEAN for testing Variance components). CLEAN-V models the global spatial dependence structure of imaging data and computes a locally powerful variance component test statistic by data-adaptively pooling neighborhood information. Correction for multiple comparisons is achieved by permutations to control family-wise error rate (FWER). Through analysis of task-fMRI data from the Human Connectome Project across five tasks and comprehensive data-driven simulations, we show that CLEAN-V outperforms existing methods in detecting test-retest reliability and narrow-sense heritability with significantly improved power, with the detected areas aligning with activation maps. The computational efficiency of CLEAN-V also speaks of its practical utility, and it is available as an R package.

Static and dynamic functional connectomes represent largely similar information

Functional connectivity (FC) of blood-oxygen-level-dependent (BOLD) fMRI time series can be estimated using methods that differ in sensitivity to the temporal order of time points (static vs. dynamic) and the number of regions considered in estimating a single edge (bivariate vs. multivariate). Previous research suggests that dynamic FC explains variability in FC fluctuations and behavior beyond static FC. Our aim was to systematically compare methods on both dimensions. We compared five FC methods: Pearson's/full correlation (static, bivariate), lagged correlation (dynamic, bivariate), partial correlation (static, multivariate) and multivariate AR model with and without self-connections (dynamic, multivariate). We compared these methods by (i) assessing similarities between FC matrices, (ii) by comparing node centrality measures, and (iii) by comparing the patterns of brain-behavior associations. Although FC estimates did not differ as a function of sensitivity to temporal order, we observed differences between the multivariate and bivariate FC methods. The dynamic FC estimates were highly correlated with the static FC estimates, especially when comparing group-level FC matrices. Similarly, there were high correlations between the patterns of brain-behavior associations obtained using the dynamic and static FC methods. We conclude that the dynamic FC estimates represent information largely similar to that of the static FC.

Heritability Estimation of Cognitive Phenotypes in the ABCD Study® Using Mixed Models

Twin and family studies have historically aimed to partition phenotypic variance into components corresponding to additive genetic effects (A), common environment (C), and unique environment (E). Here we present the ACE Model and several extensions in the Adolescent Brain Cognitive Development℠ Study (ABCD Study®), employed using the new Fast Efficient Mixed Effects Analysis (FEMA) package. In the twin sub-sample (n = 924; 462 twin pairs), heritability estimates were similar to those reported by prior studies for height (twin heritability = 0.86) and cognition (twin heritability between 0.00 and 0.61), respectively. Incorporating SNP-derived genetic relatedness and using the full ABCD Study® sample (n = 9,742) led to narrower confidence intervals for all parameter estimates. By leveraging the sparse clustering method used by FEMA to handle genetic relatedness only for participants within families, we were able to take advantage of the diverse distribution of genetic relatedness within the ABCD Study® sample.

Heritability of human "directed" functional connectome

INTRODUCTION

The functional connectivity patterns in the brain are highly heritable; however, it is unclear how genetic factors influence the directionality of such "information flows." Studying the "directionality" of the brain functional connectivity and assessing how heritability modulates it can improve our understanding of the human connectome.

METHODS

Here, we investigated the heritability of "directed" functional connections using a state-space formulation of Granger causality (GC), in conjunction with blind deconvolution methods accounting for local variability in the hemodynamic response function. Such GC implementation is ideal to explore the directionality of functional interactions across a large number of networks. Resting-state functional magnetic resonance imaging data were drawn from the Human Connectome Project (total n = 898 participants). To add robustness to our findings, the dataset was randomly split into a "discovery" and a "replication" sample (each with n = 449 participants). The two cohorts were carefully matched in terms of demographic variables and other confounding factors (e.g., education). The effect of shared environment was also modeled.

RESULTS

The parieto- and prefronto-cerebellar, parieto-prefrontal, and posterior-cingulate to hippocampus connections showed the highest and most replicable heritability effects with little influence by shared environment. In contrast, shared environmental factors significantly affected the visuo-parietal and sensory-motor directed connectivity.

CONCLUSION

We suggest a robust role of heritability in influencing the directed connectivity of some cortico-subcortical circuits implicated in cognition. Further studies, for example using task-based fMRI and GC, are warranted to confirm the asymmetric effects of genetic factors on the functional connectivity within cognitive networks and their role in supporting executive functions and learning.

A review of SNP heritability estimation methods

Over the past decade, statistical methods have been developed to estimate single nucleotide polymorphism (SNP) heritability, which measures the proportion of phenotypic variance explained by all measured SNPs in the data. Estimates of SNP heritability measure the degree to which the available genetic variants influence phenotypes and improve our understanding of the genetic architecture of complex phenotypes. In this article, we review the recently developed and commonly used SNP heritability estimation methods for continuous and binary phenotypes from the perspective of model assumptions and parameter optimization. We primarily focus on their capacity to handle multiple phenotypes and longitudinal measurements, their ability for SNP heritability partition and their use of individual-level data versus summary statistics. State-of-the-art statistical methods that are scalable to the UK Biobank dataset are also elucidated in detail.

Estimating SNP heritability in presence of population substructure in biobank-scale datasets

Single nucleotide polymorphism heritability of a trait is measured as the proportion of total variance explained by the additive effects of genome-wide single nucleotide polymorphisms. Linear mixed models are routinely used to estimate single nucleotide polymorphism heritability for many complex traits, which requires estimation of a genetic relationship matrix among individuals. Heritability is usually estimated by the restricted maximum likelihood or method of moments approaches such as Haseman-Elston regression. The common practice of accounting for such population substructure is to adjust for the top few principal components of the genetic relationship matrix as covariates in the linear mixed model. This can get computationally very intensive on large biobank-scale datasets. Here, we propose a method of moments approach for estimating single nucleotide polymorphism heritability in presence of population substructure. Our proposed method is computationally scalable on biobank datasets and gives an asymptotically unbiased estimate of heritability in presence of discrete substructures. It introduces the adjustments for population stratification in a second-order estimating equation. It allows these substructures to vary in their single nucleotide polymorphism allele frequencies and in their trait distributions (means and variances) while the heritability is assumed to be the same across these substructures. Through extensive simulation studies and the application on 7 quantitative traits in the UK Biobank cohort, we demonstrate that our proposed method performs well in the presence of population substructure and much more computationally efficient than existing approaches.

A Novel Bayesian Semi-parametric Model for Learning Heritable Imaging Traits

Heritability analysis is an important research topic in brain imaging genetics. Its primary motivation is to identify highly heritable imaging quantitative traits (QTs) for subsequent in-depth imaging genetic analyses. Most existing studies perform heritability analyses on regional imaging QTs using predefined brain parcellation schemes such as the AAL atlas. However, the power to dissect genetic underpinnings under QTs defined in such an unsupervised fashion is largely deteriorate with inner partition noise and signal dilution. To bridge the gap, we propose a new semi-parametric Bayesian heritability estimation model to construct highly heritable imaging QTs. Our method leverages the aggregate of genetic signals to imaging QT construction by developing a new brain parcellation driven by voxel-level heritability. To ensure biological plausibility and clinical interpretability of the resulting brain heritability parcellations, hierarchical sparsity and smoothness, coupled with structural connectivity of the brain, are properly imposed on genetic effects to induce spatial contiguity of heritable imaging QTs. Using the ADNI imaging genetic data, we demonstrate the strength of our proposed method, in comparison with the standard GCTA method, in identifying highly heritable and biologically meaningful new imaging QTs.

Polygenic Architecture of Human Neuroanatomical Diversity

We analyzed the genomic architecture of neuroanatomical diversity using magnetic resonance imaging and single nucleotide polymorphism (SNP) data from >26 000 individuals from the UK Biobank project and 5 other projects that had previously participated in the ENIGMA (Enhancing NeuroImaging Genetics through Meta-Analysis) consortium. Our results confirm the polygenic architecture of neuroanatomical diversity, with SNPs capturing from 40% to 54% of regional brain volume variance. Chromosomal length correlated with the amount of phenotypic variance captured, r ~ 0.64 on average, suggesting that at a global scale causal variants are homogeneously distributed across the genome. At a local scale, SNPs within genes (~51%) captured ~1.5 times more genetic variance than the rest, and SNPs with low minor allele frequency (MAF) captured less variance than the rest: the 40% of SNPs with MAF <5% captured <one fourth of the genetic variance. We also observed extensive pleiotropy across regions, with an average genetic correlation of rG ~ 0.45. Genetic correlations were similar to phenotypic and environmental correlations; however, genetic correlations were often larger than phenotypic correlations for the left/right volumes of the same region. The heritability of differences in left/right volumes was generally not statistically significant, suggesting an important influence of environmental causes in the variability of brain asymmetry. Our code is available athttps://github.com/neuroanatomy/genomic-architecture.

Unpaired data empowers association tests

Motivation

There is growing interest in the biomedical research community to incorporate retrospective data, available in healthcare systems, to shed light on associations between different biomarkers. Understanding the association between various types of biomedical data, such as genetic, blood biomarkers, imaging, etc. can provide a holistic understanding of human diseases. To formally test a hypothesized association between two types of data in Electronic Health Records (EHRs), one requires a substantial sample size with both data modalities to achieve a reasonable power. Current association test methods only allow using data from individuals who have both data modalities. Hence, researchers cannot take advantage of much larger EHR samples that includes individuals with at least one of the data types, which limits the power of the association test.

Results

We present a new method called the Semi-paired Association Test (SAT) that makes use of both paired and unpaired data. In contrast to classical approaches, incorporating unpaired data allows SAT to produce better control of false discovery and to improve the power of the association test. We study the properties of the new test theoretically and empirically, through a series of simulations and by applying our method on real studies in the context of Chronic Obstructive Pulmonary Disease. We are able to identify an association between the high-dimensional characterization of Computed Tomography chest images and several blood biomarkers as well as the expression of dozens of genes involved in the immune system.

Availability and implementation

Code is available on https://github.com/batmanlab/Semi-paired-Association-Test.

Supplementary information

Supplementary data are available at Bioinformatics online.

Single-Cell RNA Sequencing in Parkinson's Disease

Single-cell and single-nucleus RNA sequencing (sc/snRNA-seq) technologies have enhanced the understanding of the molecular pathogenesis of neurodegenerative disorders, including Parkinson's disease (PD). Nonetheless, their application in PD has been limited due mainly to the technical challenges resulting from the scarcity of postmortem brain tissue and low quality associated with RNA degradation. Despite such challenges, recent advances in animals and human in vitro models that recapitulate features of PD along with sequencing assays have fueled studies aiming to obtain an unbiased and global view of cellular composition and phenotype of PD at the single-cell resolution. Here, we reviewed recent sc/snRNA-seq efforts that have successfully characterized diverse cell-type populations and identified cell type-specific disease associations in PD. We also examined how these studies have employed computational and analytical tools to analyze and interpret the rich information derived from sc/snRNA-seq. Finally, we highlighted important limitations and emerging technologies for addressing key technical challenges currently limiting the integration of new findings into clinical practice.

Heritability of individualized cortical network topography

Human cortex is patterned by a complex and interdigitated web of large-scale functional networks. Recent methodological breakthroughs reveal variation in the size, shape, and spatial topography of cortical networks across individuals. While spatial network organization emerges across development, is stable over time, and is predictive of behavior, it is not yet clear to what extent genetic factors underlie interindividual differences in network topography. Here, leveraging a nonlinear multidimensional estimation of heritability, we provide evidence that individual variability in the size and topographic organization of cortical networks are under genetic control. Using twin and family data from the Human Connectome Project (n = 1,023), we find increased variability and reduced heritability in the size of heteromodal association networks (h2 : M = 0.34, SD = 0.070), relative to unimodal sensory/motor cortex (h2 : M = 0.40, SD = 0.097). We then demonstrate that the spatial layout of cortical networks is influenced by genetics, using our multidimensional estimation of heritability (h2-multi; M = 0.14, SD = 0.015). However, topographic heritability did not differ between heteromodal and unimodal networks. Genetic factors had a regionally variable influence on brain organization, such that the heritability of network topography was greatest in prefrontal, precuneus, and posterior parietal cortex. Taken together, these data are consistent with relaxed genetic control of association cortices relative to primary sensory/motor regions and have implications for understanding population-level variability in brain functioning, guiding both individualized prediction and the interpretation of analyses that integrate genetics and neuroimaging.

Heritability of Regional Brain Volumes in Large-Scale Neuroimaging and Genetic Studies

Brain genetics is an active research area. The degree to which genetic variants impact variations in brain structure and function remains largely unknown. We examined the heritability of regional brain volumes (P ~ 100) captured by single-nucleotide polymorphisms (SNPs) in UK Biobank (n ~ 9000). We found that regional brain volumes are highly heritable in this study population and common genetic variants can explain up to 80% of their variabilities (median heritability 34.8%). We observed omnigenic impact across the genome and examined the enrichment of SNPs in active chromatin regions. Principal components derived from regional volume data are also highly heritable, but the amount of variance in brain volume explained by the component did not seem to be related to its heritability. Heritability estimates vary substantially across large-scale functional networks, exhibit a symmetric pattern across left and right hemispheres, and are consistent in females and males (correlation = 0.638). We repeated the main analysis in Alzheimer's Disease Neuroimaging Initiative (n ~ 1100), Philadelphia Neurodevelopmental Cohort (n ~ 600), and Pediatric Imaging, Neurocognition, and Genetics (n ~ 500) datasets, which demonstrated that more stable estimates can be obtained from the UK Biobank.

Bayesian sparse heritability analysis with high-dimensional neuroimaging phenotypes

Heritability analysis plays a central role in quantitative genetics to describe genetic contribution to human complex traits and prioritize downstream analyses under large-scale phenotypes. Existing works largely focus on modeling single phenotype and currently available multivariate phenotypic methods often suffer from scaling and interpretation. In this article, motivated by understanding how genetic underpinning impacts human brain variation, we develop an integrative Bayesian heritability analysis to jointly estimate heritabilities for high-dimensional neuroimaging traits. To induce sparsity and incorporate brain anatomical configuration, we impose hierarchical selection among both regional and local measurements based on brain structural network and voxel dependence. We also use a nonparametric Dirichlet process mixture model to realize grouping among single nucleotide polymorphism-associated phenotypic variations, providing biological plausibility. Through extensive simulations, we show the proposed method outperforms existing ones in heritability estimation and heritable traits selection under various scenarios. We finally apply the method to two large-scale imaging genetics datasets: the Alzheimer's Disease Neuroimaging Initiative and United Kingdom Biobank and show biologically meaningful results.

The Shared Genetic Basis of Educational Attainment and Cerebral Cortical Morphology

Individual differences in educational attainment are linked to differences in intelligence, and predict important social, economic, and health outcomes. Previous studies have found common genetic factors that influence educational achievement, cognitive performance and total brain volume (i.e., brain size). Here, in a large sample of participants from the UK Biobank, we investigate the shared genetic basis between educational attainment and fine-grained cerebral cortical morphological features, and associate this genetic variation with a related aspect of cognitive ability. Importantly, we execute novel statistical methods that enable high-dimensional genetic correlation analysis, and compute high-resolution surface maps for the genetic correlations between educational attainment and vertex-wise morphological measurements. We conduct secondary analyses, using the UK Biobank verbal-numerical reasoning score, to confirm that variation in educational attainment that is genetically correlated with cortical morphology is related to differences in cognitive performance. Our analyses relate the genetic overlap between cognitive ability and cortical thickness measurements to bilateral primary motor cortex as well as predominantly left superior temporal cortex and proximal regions. These findings extend our understanding of the neurobiology that connects genetic variation to individual differences in educational attainment and cognitive performance.

Imaging local genetic influences on cortical folding

Recent progress in deciphering mechanisms of human brain cortical folding leave unexplained whether spatially patterned genetic influences contribute to this folding. High-resolution in vivo brain MRI can be used to estimate genetic correlations (covariability due to shared genetic factors) in interregional cortical thickness, and biomechanical studies predict an influence of cortical thickness on folding patterns. However, progress has been hampered because shared genetic influences related to folding patterns likely operate at a scale that is much more local (<1 cm) than that addressed in prior imaging studies. Here, we develop methodological approaches to examine local genetic influences on cortical thickness and apply these methods to two large, independent samples. We find that such influences are markedly heterogeneous in strength, and in some cortical areas are notably stronger in specific orientations relative to gyri or sulci. The overall, phenotypic local correlation has a significant basis in shared genetic factors and is highly symmetric between left and right cortical hemispheres. Furthermore, the degree of local cortical folding relates systematically with the strength of local correlations, which tends to be higher in gyral crests and lower in sulcal fundi. The relationship between folding and local correlations is stronger in primary sensorimotor areas and weaker in association areas such as prefrontal cortex, consistent with reduced genetic constraints on the structural topology of association cortex. Collectively, our results suggest that patterned genetic influences on cortical thickness, measurable at the scale of in vivo MRI, may be a causal factor in the development of cortical folding.

Brain Imaging Genomics: Integrated Analysis and Machine Learning

Brain imaging genomics is an emerging data science field, where integrated analysis of brain imaging and genomics data, often combined with other biomarker, clinical and environmental data, is performed to gain new insights into the phenotypic, genetic and molecular characteristics of the brain as well as their impact on normal and disordered brain function and behavior. It has enormous potential to contribute significantly to biomedical discoveries in brain science. Given the increasingly important role of statistical and machine learning in biomedicine and rapidly growing literature in brain imaging genomics, we provide an up-to-date and comprehensive review of statistical and machine learning methods for brain imaging genomics, as well as a practical discussion on method selection for various biomedical applications.

Beyond heritability: improving discoverability in imaging genetics

Structural neuroimaging measures based on magnetic resonance imaging have been at the forefront of imaging genetics. Global efforts to ensure homogeneity of measurements across study sites have enabled large-scale imaging genetic projects, accumulating nearly 50K samples for genome-wide association studies (GWAS). However, not many novel genetic variants have been identified by these GWAS, despite the high heritability of structural neuroimaging measures. Here, we discuss the limitations of using heritability as a guidance for assessing statistical power of GWAS, and highlight the importance of discoverability-which is the power to detect genetic variants for a given phenotype depending on its unique genomic architecture and GWAS sample size. Further, we present newly developed methods that boost genetic discovery in imaging genetics. By redefining imaging measures independent of traditional anatomical conventions, it is possible to improve discoverability, enabling identification of more genetic effects. Moreover, by leveraging enrichment priors from genomic annotations and independent GWAS of pleiotropic traits, we can better characterize effect size distributions, and identify reliable and replicable loci associated with structural neuroimaging measures. Statistical tools leveraging novel insights into the genetic discoverability of human traits, promises to accelerate the identification of genetic underpinnings underlying brain structural variation.

Neonatal amygdalae and hippocampi are influenced by genotype and prenatal environment, and reflected in the neonatal DNA methylome

The amygdala and hippocampus undergo rapid development in early life. The relative contribution of genetic and environmental factors to the establishment of their developmental trajectories has yet to be examined. We performed imaging on neonates and examined how the observed variation in volume and microstructure of the amygdala and hippocampus varied by genotype, and compared with prenatal maternal mental health and socioeconomic status. Gene × Environment models outcompeted models containing genotype or environment only to best explain the majority of measures but some, especially of the amygdaloid microstructure, were best explained by genotype only. Models including DNA methylation measured in the neonate umbilical cords outcompeted the Gene and Gene × Environment models for the majority of amygdaloid measures and minority of hippocampal measures. This study identified brain region-specific gene networks associated with individual differences in fetal brain development. In particular, genetic and epigenetic variation within CUX1 was highlighted.

Reliable heritability estimation using sparse regularization in ultrahigh dimensional genome-wide association studies

BACKGROUND

Data from genome-wide association studies (GWASs) have been used to estimate the heritability of human complex traits in recent years. Existing methods are based on the linear mixed model, with the assumption that the genetic effects are random variables, which is opposite to the fixed effect assumption embedded in the framework of quantitative genetics theory. Moreover, heritability estimators provided by existing methods may have large standard errors, which calls for the development of reliable and accurate methods to estimate heritability.

RESULTS

In this paper, we first investigate the influences of the fixed and random effect assumption on heritability estimation, and prove that these two assumptions are equivalent under mild conditions in the theoretical aspect. Second, we propose a two-stage strategy by first performing sparse regularization via cross-validated elastic net, and then applying variance estimation methods to construct reliable heritability estimations. Results on both simulated data and real data show that our strategy achieves a considerable reduction in the standard error while reserving the accuracy.

CONCLUSIONS

The proposed strategy allows for a reliable and accurate heritability estimation using GWAS data. It shows the promising future that reliable estimations can still be obtained with even a relatively restricted sample size, and should be especially useful for large-scale heritability analyses in the genomics era.

Accelerated estimation and permutation inference for ACE modeling

There are a wealth of tools for fitting linear models at each location in the brain in neuroimaging analysis, and a wealth of genetic tools for estimating heritability for a small number of phenotypes. But there remains a need for computationally efficient neuroimaging genetic tools that can conduct analyses at the brain‐wide scale. Here we present a simple method for heritability estimation on twins that replaces a variance component model‐which requires iterative optimisation‐with a (noniterative) linear regression model, by transforming data to squared twin‐pair differences. We demonstrate that the method has comparable bias, mean squared error, false positive risk, and power to best practice maximum‐likelihood‐based methods, while requiring a small fraction of the computation time. Combined with permutation, we call this approach “Accelerated Permutation Inference for the ACE Model (APACE)” where ACE refers to the additive genetic (A) effects, and common (C), and unique (E) environmental influences on the trait. We show how the use of spatial statistics like cluster size can dramatically improve power, and illustrate the method on a heritability analysis of an fMRI working memory dataset.

Detecting heritable phenotypes without a model using fast permutation testing for heritability and set-tests

Testing for association between a set of genetic markers and a phenotype is a fundamental task in genetic studies. Standard approaches for heritability and set testing strongly rely on parametric models that make specific assumptions regarding phenotypic variability. Here, we show that resulting p-values may be inflated by up to 15 orders of magnitude, in a heritability study of methylation measurements, and in a heritability and expression quantitative trait loci analysis of gene expression profiles. We propose FEATHER, a method for fast permutation-based testing of marker sets and of heritability, which properly controls for false-positive results. FEATHER eliminated 47% of methylation sites found to be heritable by the parametric test, suggesting a substantial inflation of false-positive findings by alternative methods. Our approach can rapidly identify heritable phenotypes out of millions of phenotypes acquired via high-throughput technologies, does not suffer from model misspecification and is highly efficient.

Standard approaches for heritability and set testing in statistical genetics rely on parametric models that might not hold in reality and give inflated p-values. Here, the authors develop a fast method for permutation-based testing of marker sets and of heritability that does not suffer from model misspecification.

The use of electronic health records for psychiatric phenotyping and genomics

The widespread adoption of electronic health record (EHRs) in healthcare systems has created a vast and continuously growing resource of clinical data and provides new opportunities for population-based research. In particular, the linking of EHRs to biospecimens and genomic data in biobanks may help address what has become a rate-limiting study for genetic research: the need for large sample sizes. The principal roadblock to capitalizing on these resources is the need to establish the validity of phenotypes extracted from the EHR. For psychiatric genetic research, this represents a particular challenge given that diagnosis is based on patient reports and clinician observations that may not be well-captured in billing codes or narrative records. This review addresses the opportunities and pitfalls in EHR-based phenotyping with a focus on their application to psychiatric genetic research. A growing number of studies have demonstrated that diagnostic algorithms with high positive predictive value can be derived from EHRs, especially when structured data are supplemented by text mining approaches. Such algorithms enable semi-automated phenotyping for large-scale case-control studies. In addition, the scale and scope of EHR databases have been used successfully to identify phenotypic subgroups and derive algorithms for longitudinal risk prediction. EHR-based genomics are particularly well-suited to rapid look-up replication of putative risk genes, studies of pleiotropy (phenomewide association studies or PheWAS), investigations of genetic networks and overlap across the phenome, and pharmacogenomic research. EHR phenotyping has been relatively under-utilized in psychiatric genomic research but may become a key component of efforts to advance precision psychiatry.

Recent developments in the genetics of attention-deficit hyperactivity disorder

Attention-deficit hyperactivity disorder (ADHD) is a developmental psychiatric disorder that affects children and adults. ADHD is one of the psychiatric disorders with the strongest genetic basis according to familial, twin, and single nucleotide polymorphisms (SNP)-based epidemiological studies. In this review, we provide an update of recent insights into the genetic basis of ADHD. We discuss recent progress from genome-wide association studies (GWAS) looking at common variants as well as rare copy number variations. New analysis of gene groups, so-called functional ontologies, provide some insight into the gene networks afflicted, pointing to the role of neurodevelopmentally expressed gene networks. Bioinformatic methods, such as functional enrichment analysis and protein-protein network analysis, are used to highlight biological processes of likely relevance to the etiology of ADHD. Additionally, copy number variations seem to map on important pathways implicated in synaptic signaling and neurodevelopment. While some candidate gene associations of, for example, neurotransmitter receptors and signaling, have been replicated, they do not seem to explain significant variance in recent GWAS. We discuss insights from recent case-control SNP-GWAS that have presented the first whole-genome significant SNP in ADHD.

The Genetics of Endophenotypes of Neurofunction to Understand Schizophrenia (GENUS) consortium: A collaborative cognitive and neuroimaging genetics project

BACKGROUND

Schizophrenia has a large genetic component, and the pathways from genes to illness manifestation are beginning to be identified. The Genetics of Endophenotypes of Neurofunction to Understand Schizophrenia (GENUS) Consortium aims to clarify the role of genetic variation in brain abnormalities underlying schizophrenia. This article describes the GENUS Consortium sample collection.

METHODS

We identified existing samples collected for schizophrenia studies consisting of patients, controls, and/or individuals at familial high-risk (FHR) for schizophrenia. Samples had single nucleotide polymorphism (SNP) array data or genomic DNA, clinical and demographic data, and neuropsychological and/or brain magnetic resonance imaging (MRI) data. Data were subjected to quality control procedures at a central site.

RESULTS

Sixteen research groups contributed data from 5199 psychosis patients, 4877 controls, and 725 FHR individuals. All participants have relevant demographic data and all patients have relevant clinical data. The sex ratio is 56.5% male and 43.5% female. Significant differences exist between diagnostic groups for premorbid and current IQ (both p<1×10^-10). Data from a diversity of neuropsychological tests are available for 92% of participants, and 30% have structural MRI scans (half also have diffusion-weighted MRI scans). SNP data are available for 76% of participants. The ancestry composition is 70% European, 20% East Asian, 7% African, and 3% other.

CONCLUSIONS

The Consortium is investigating the genetic contribution to brain phenotypes in a schizophrenia sample collection of >10,000 participants. The breadth of data across clinical, genetic, neuropsychological, and MRI modalities provides an important opportunity for elucidating the genetic basis of neural processes underlying schizophrenia.

Another Round of "Clue" to Uncover the Mystery of Complex Traits

A plethora of genetic association analyses have identified several genetic risk loci. Technological and statistical advancements have now led to the identification of not only common genetic variants, but also low-frequency variants, structural variants, and environmental factors, as well as multi-omics variations that affect the phenotypic variance of complex traits in a population, thus referred to as complex trait architecture. The concept of heritability, or the proportion of phenotypic variance due to genetic inheritance, has been studied for several decades, but its application is mainly in addressing the narrow sense heritability (or additive genetic component) from Genome-Wide Association Studies (GWAS). In this commentary, we reflect on our perspective on the complexity of understanding heritability for human traits in comparison to model organisms, highlighting another round of clues beyond GWAS and an alternative approach, investigating these clues comprehensively to help in elucidating the genetic architecture of complex traits.

Sample composition alters associations between age and brain structure

Despite calls to incorporate population science into neuroimaging research, most studies recruit small, non-representative samples. Here, we examine whether sample composition influences age-related variation in global measurements of gray matter volume, thickness, and surface area. We apply sample weights to structural brain imaging data from a community-based sample of children aged 3–18 (N = 1162) to create a “weighted sample” that approximates the distribution of socioeconomic status, race/ethnicity, and sex in the U.S. Census. We compare associations between age and brain structure in this weighted sample to estimates from the original sample with no sample weights applied (i.e., unweighted). Compared to the unweighted sample, we observe earlier maturation of cortical and sub-cortical structures, and patterns of brain maturation that better reflect known developmental trajectories in the weighted sample. Our empirical demonstration of bias introduced by non-representative sampling in this neuroimaging cohort suggests that sample composition may influence understanding of fundamental neural processes.

The influence of sample composition on human neuroimaging results is unknown. Here, the authors weight a large, community-based sample to better reflect the US population and describe how applying these sample weights changes conclusions about age-related variation in brain structure.

RL-SKAT: An Exact and Efficient Score Test for Heritability and Set Tests

Testing for the existence of variance components in linear mixed models is a fundamental task in many applicative fields. In statistical genetics, the score test has recently become instrumental in the task of testing an association between a set of genetic markers and a phenotype. With few markers, this amounts to set-based variance component tests, which attempt to increase power in association studies by aggregating weak individual effects. When the entire genome is considered, it allows testing for the heritability of a phenotype, defined as the proportion of phenotypic variance explained by genetics. In the popular score-based Sequence Kernel Association Test (SKAT) method, the assumed distribution of the score test statistic is uncalibrated in small samples, with a correction being computationally expensive. This may cause severe inflation or deflation of P-values, even when the null hypothesis is true. Here, we characterize the conditions under which this discrepancy holds, and show it may occur also in large real datasets, such as a dataset from the Wellcome Trust Case Control Consortium 2 (n = 13,950) study, and, in particular, when the individuals in the sample are unrelated. In these cases, the SKAT approximation tends to be highly overconservative and therefore underpowered. To address this limitation, we suggest an efficient method to calculate exact P-values for the score test in the case of a single variance component and a continuous response vector, which can speed up the analysis by orders of magnitude. Our results enable fast and accurate application of the score test in heritability and in set-based association tests. Our method is available in http://github.com/cozygene/RL-SKAT.

FGWAS: Functional genome wide association analysis

Functional phenotypes (e.g., subcortical surface representation), which commonly arise in imaging genetic studies, have been used to detect putative genes for complexly inherited neuropsychiatric and neurodegenerative disorders. However, existing statistical methods largely ignore the functional features (e.g., functional smoothness and correlation). The aim of this paper is to develop a functional genome-wide association analysis (FGWAS) framework to efficiently carry out whole-genome analyses of functional phenotypes. FGWAS consists of three components: a multivariate varying coefficient model, a global sure independence screening procedure, and a test procedure. Compared with the standard multivariate regression model, the multivariate varying coefficient model explicitly models the functional features of functional phenotypes through the integration of smooth coefficient functions and functional principal component analysis. Statistically, compared with existing methods for genome-wide association studies (GWAS), FGWAS can substantially boost the detection power for discovering important genetic variants influencing brain structure and function. Simulation studies show that FGWAS outperforms existing GWAS methods for searching sparse signals in an extremely large search space, while controlling for the family-wise error rate. We have successfully applied FGWAS to large-scale analysis of data from the Alzheimer's Disease Neuroimaging Initiative for 708 subjects, 30,000 vertices on the left and right hippocampal surfaces, and 501,584 SNPs.

Do regional brain volumes and major depressive disorder share genetic architecture? A study of Generation Scotland (n=19 762), UK Biobank (n=24 048) and the English Longitudinal Study of Ageing (n=5766)

Major depressive disorder (MDD) is a heritable and highly debilitating condition. It is commonly associated with subcortical volumetric abnormalities, the most replicated of these being reduced hippocampal volume. Using the most recent published data from Enhancing Neuroimaging Genetics through Meta-analysis (ENIGMA) consortium's genome-wide association study of regional brain volume, we sought to test whether there is shared genetic architecture between seven subcortical brain volumes and intracranial volume (ICV) and MDD. We explored this using linkage disequilibrium score regression, polygenic risk scoring (PRS) techniques, Mendelian randomisation (MR) analysis and BUHMBOX. Utilising summary statistics from ENIGMA and Psychiatric Genomics Consortium, we demonstrated that hippocampal volume was positively genetically correlated with MDD (r_G=0.46, P=0.02), although this did not survive multiple comparison testing. None of the other six brain regions studied were genetically correlated and amygdala volume heritability was too low for analysis. Using PRS analysis, no regional volumetric PRS demonstrated a significant association with MDD or recurrent MDD. MR analysis in hippocampal volume and MDD identified no causal association, however, BUHMBOX analysis identified genetic subgrouping in GS:SFHS MDD cases only (P=0.00281). In this study, we provide some evidence that hippocampal volume and MDD may share genetic architecture in a subgroup of individuals, albeit the genetic correlation did not survive multiple testing correction and genetic subgroup heterogeneity was not replicated. In contrast, we found no evidence to support a shared genetic architecture between MDD and other regional subcortical volumes or ICV.

Heritability analysis with repeat measurements and its application to resting-state functional connectivity

Heritability, defined as the proportion of phenotypic variation attributable to genetic variation, provides important information about the genetic basis of a trait. Existing heritability analysis methods do not discriminate between stable effects (e.g., due to the subject's unique environment) and transient effects, such as measurement error. This can lead to misleading assessments, particularly when comparing the heritability of traits that exhibit different levels of reliability. Here, we present a linear mixed effects model to conduct heritability analyses that explicitly accounts for intrasubject fluctuations (e.g., due to measurement noise or biological transients) using repeat measurements. We apply the proposed strategy to the analysis of resting-state fMRI measurements-a prototypic data modality that exhibits variable levels of test-retest reliability across space. Our results reveal that the stable components of functional connectivity within and across well-established large-scale brain networks can be considerably heritable. Furthermore, we demonstrate that dissociating intra- and intersubject variation can reveal genetic influence on a phenotype that is not fully captured by conventional heritability analyses.

A peripheral epigenetic signature of immune system genes is linked to neocortical thickness and memory

Increasing age is tightly linked to decreased thickness of the human neocortex. The biological mechanisms that mediate this effect are hitherto unknown. The DNA methylome, as part of the epigenome, contributes significantly to age-related phenotypic changes. Here, we identify an epigenetic signature that is associated with cortical thickness (P=3.86 × 10⁻⁸) and memory performance in 533 healthy young adults. The epigenetic effect on cortical thickness was replicated in a sample comprising 596 participants with major depressive disorder and healthy controls. The epigenetic signature mediates partially the effect of age on cortical thickness (P<0.001). A multilocus genetic score reflecting genetic variability of this signature is associated with memory performance (P=0.0003) in 3,346 young and elderly healthy adults. The genomic location of the contributing methylation sites points to the involvement of specific immune system genes. The decomposition of blood methylome-wide patterns bears considerable potential for the study of brain-related traits.

Cortical thickness has high heritability estimates and is known to be influenced by genetic factors. Here, Freytag and colleagues show that DNA methylation patterns of peripheral blood monocytes are also correlated with cortical thickness and memory performance in human.

Phenome-wide heritability analysis of the UK Biobank

Heritability estimation provides important information about the relative contribution of genetic and environmental factors to phenotypic variation, and provides an upper bound for the utility of genetic risk prediction models. Recent technological and statistical advances have enabled the estimation of additive heritability attributable to common genetic variants (SNP heritability) across a broad phenotypic spectrum. Here, we present a computationally and memory efficient heritability estimation method that can handle large sample sizes, and report the SNP heritability for 551 complex traits derived from the interim data release (152,736 subjects) of the large-scale, population-based UK Biobank, comprising both quantitative phenotypes and disease codes. We demonstrate that common genetic variation contributes to a broad array of quantitative traits and human diseases in the UK population, and identify phenotypes whose heritability is moderated by age (e.g., a majority of physical measures including height and body mass index), sex (e.g., blood pressure related traits) and socioeconomic status (education). Our study represents the first comprehensive phenome-wide heritability analysis in the UK Biobank, and underscores the importance of considering population characteristics in interpreting heritability.

Heritability of a trait refers to the proportion of phenotypic variation that is due to genetic variation among individuals. It provides important information about the genetic basis of complex traits and indicates whether a phenotype is an appropriate target for more specific statistical and molecular genetic analyses. Recent studies have leveraged the increasingly ubiquitous genome-wide data and documented the heritability attributable to common genetic variation captured by genotyping microarrays for a wide range of human traits. However, heritability is not a fixed property of a phenotype and can vary with population-specific differences in the genetic background and environmental variation. Here, using a computationally and memory efficient heritability estimation method, we report the heritability for a large number of traits derived from the large-scale, population-based UK Biobank, and, for the first time, demonstrate the moderating effect of three major demographic variables (age, sex and socioeconomic status) on heritability estimates derived from genome-wide common genetic variation. Our study represents the first comprehensive heritability analysis across the phenotypic spectrum in the UK Biobank.

Connectome Disconnectivity and Cortical Gene Expression in Patients With Schizophrenia

BACKGROUND

Genome-wide association studies have identified several common risk loci for schizophrenia (SCZ). In parallel, neuroimaging studies have shown consistent findings of widespread white matter disconnectivity in patients with SCZ.

METHODS

We examined the role of genes in brain connectivity in patients with SCZ by combining transcriptional profiles of 43 SCZ risk genes identified by the recent genome-wide association study of the Schizophrenia Working Group of the Psychiatric Genomics Consortium with data on macroscale connectivity reductions in patients with SCZ. Expression profiles of 43 Psychiatric Genomics Consortium SCZ risk genes were extracted from the Allen Human Brain Atlas, and their average profile across the cortex was correlated to the pattern of cortical disconnectivity as derived from diffusion-weighted magnetic resonance imaging data of patients with SCZ (n = 48) and matched healthy controls (n = 43).

RESULTS

The expression profile of SCZ risk genes across cortical regions was significantly correlated with the regional macroscale disconnectivity (r = .588; p = .017). In addition, effects were found to be potentially specific to SCZ, with transcriptional profiles not related to cortical disconnectivity in patients with bipolar I disorder (diffusion-weighted magnetic resonance imaging data; 216 patients, 144 controls). Further examination of correlations across all 20,737 genes present in the Allen Human Brain Atlas showed the set of top 100 strongest correlating genes to display significant enrichment for the disorder, potentially identifying new genes involved in the pathophysiology of SCZ.

CONCLUSIONS

Our results suggest that under disease conditions, cortical areas with pronounced expression of risk genes implicated in SCZ form central areas for white matter disconnectivity.

ENIGMA and the individual: Predicting factors that affect the brain in 35 countries worldwide

In this review, we discuss recent work by the ENIGMA Consortium (http://enigma.ini.usc.edu) - a global alliance of over 500 scientists spread across 200 institutions in 35 countries collectively analyzing brain imaging, clinical, and genetic data. Initially formed to detect genetic influences on brain measures, ENIGMA has grown to over 30 working groups studying 12 major brain diseases by pooling and comparing brain data. In some of the largest neuroimaging studies to date - of schizophrenia and major depression - ENIGMA has found replicable disease effects on the brain that are consistent worldwide, as well as factors that modulate disease effects. In partnership with other consortia including ADNI, CHARGE, IMAGEN and others¹, ENIGMA's genomic screens - now numbering over 30,000 MRI scans - have revealed at least 8 genetic loci that affect brain volumes. Downstream of gene findings, ENIGMA has revealed how these individual variants - and genetic variants in general - may affect both the brain and risk for a range of diseases. The ENIGMA consortium is discovering factors that consistently affect brain structure and function that will serve as future predictors linking individual brain scans and genomic data. It is generating vast pools of normative data on brain measures - from tens of thousands of people - that may help detect deviations from normal development or aging in specific groups of subjects. We discuss challenges and opportunities in applying these predictors to individual subjects and new cohorts, as well as lessons we have learned in ENIGMA's efforts so far.

Endophenotype best practices

This review examines the current state of electrophysiological endophenotype research and recommends best practices that are based on knowledge gleaned from the last decade of molecular genetic research with complex traits. Endophenotype research is being oversold for its potential to help discover psychopathology relevant genes using the types of small samples feasible for electrophysiological research. This is largely because the genetic architecture of endophenotypes appears to be very much like that of behavioral traits and disorders: they are complex, influenced by many variants (e.g., tens of thousands) within many genes, each contributing a very small effect. Out of over 40 electrophysiological endophenotypes covered by our review, only resting heart, a measure that has received scant advocacy as an endophenotype, emerges as an electrophysiological variable with verified associations with molecular genetic variants. To move the field forward, investigations designed to discover novel variants associated with endophenotypes will need extremely large samples best obtained by forming consortia and sharing data obtained from genome wide arrays. In addition, endophenotype research can benefit from successful molecular genetic studies of psychopathology by examining the degree to which these verified psychopathology-relevant variants are also associated with an endophenotype, and by using knowledge about the functional significance of these variants to generate new endophenotypes. Even without molecular genetic associations, endophenotypes still have value in studying the development of disorders in unaffected individuals at high genetic risk, constructing animal models, and gaining insight into neural mechanisms that are relevant to clinical disorder.

Heritability of the shape of subcortical brain structures in the general population

The volumes of subcortical brain structures are highly heritable, but genetic underpinnings of their shape remain relatively obscure. Here we determine the relative contribution of genetic factors to individual variation in the shape of seven bilateral subcortical structures: the nucleus accumbens, amygdala, caudate, hippocampus, pallidum, putamen and thalamus. In 3,686 unrelated individuals aged between 45 and 98 years, brain magnetic resonance imaging and genotyping was performed. The maximal heritability of shape varies from 32.7 to 53.3% across the subcortical structures. Genetic contributions to shape extend beyond influences on intracranial volume and the gross volume of the respective structure. The regional variance in heritability was related to the reliability of the measurements, but could not be accounted for by technical factors only. These findings could be replicated in an independent sample of 1,040 twins. Differences in genetic contributions within a single region reveal the value of refined brain maps to appreciate the genetic complexity of brain structures.

Partitioning heritability analysis reveals a shared genetic basis of brain anatomy and schizophrenia

Schizophrenia is a devastating neurodevelopmental disorder with a complex genetic etiology. Widespread cortical gray matter loss has been observed in patients and prodromal samples. However, it remains unresolved whether schizophrenia-associated cortical structure variations arise due to disease etiology or secondary to the illness. Here we address this question using a partitioning-based heritability analysis of genome-wide single-nucleotide polymorphism (SNP) and neuroimaging data from 1750 healthy individuals. We find that schizophrenia-associated genetic variants explain a significantly enriched proportion of trait heritability in eight brain phenotypes (false discovery rate=10%). In particular, intracranial volume and left superior frontal gyrus thickness exhibit significant and robust associations with schizophrenia genetic risk under varying SNP selection conditions. Cross-disorder comparison suggests that the neurogenetic architecture of schizophrenia-associated brain regions is, at least in part, shared with other psychiatric disorders. Our study highlights key neuroanatomical correlates of schizophrenia genetic risk in the general population. These may provide fundamental insights into the complex pathophysiology of the illness, and a potential link to neurocognitive deficits shaping the disorder.

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.

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.

Morphometricity as a measure of the neuroanatomical signature of a trait

Complex physiological and behavioral traits, including neurological and psychiatric disorders, often associate with distributed anatomical variation. This paper introduces a global metric, called morphometricity, as a measure of the anatomical signature of different traits. Morphometricity is defined as the proportion of phenotypic variation that can be explained by macroscopic brain morphology. We estimate morphometricity via a linear mixed-effects model that uses an anatomical similarity matrix computed based on measurements derived from structural brain MRI scans. We examined over 3,800 unique MRI scans from nine large-scale studies to estimate the morphometricity of a range of phenotypes, including clinical diagnoses such as Alzheimer's disease, and nonclinical traits such as measures of cognition. Our results demonstrate that morphometricity can provide novel insights about the neuroanatomical correlates of a diverse set of traits, revealing associations that might not be detectable through traditional statistical techniques.

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.

Polygenic Epidemiology

Much of the genetic basis of complex traits is present on current genotyping products, but the individual variants that affect the traits have largely not been identified. Several traditional problems in genetic epidemiology have recently been addressed by assuming a polygenic basis for disease and treating it as a single entity. Here I briefly review some of these applications, which collectively may be termed polygenic epidemiology. Methodologies in this area include polygenic scoring, linear mixed models, and linkage disequilibrium scoring. They have been used to establish a polygenic effect, estimate genetic correlation between traits, estimate how many variants affect a trait, stratify cases into subphenotypes, predict individual disease risks, and infer causal effects using Mendelian randomization. Polygenic epidemiology will continue to yield useful applications even while much of the specific variation underlying complex traits remains undiscovered.