Genomic Structural Equation Modeling Reveals Latent Phenotypes in the Human Cortex with Distinct Genetic Architecture

Abstract

ABSTRACT Genetic contributions to human cortical structure manifest pervasive pleiotropy. This pleiotropy may be harnessed to identify unique genetically-informed parcellations of the cortex that are neurobiologically distinct from anatomical, functional, cytoarchitectural, or other cortical parcellation schemes. We investigated genetic pleiotropy by applying genomic structural equation modeling (SEM) to model the genetic architecture of cortical surface area (SA) and cortical thickness (CT) of 34 brain regions recently reported in the ENIGMA cortical GWAS. Genomic SEM uses the empirical genetic covariance estimated from GWAS summary statistics with LD score regression (LDSC) to discover factors underlying genetic covariance. Genomic SEM can fit a multivariate GWAS from summary statistics, which can subsequently be used for LD score regression (LDSC). We found the best-fitting model of cortical SA was explained by 6 latent factors and CT was explained by 4 latent factors. The multivariate GWAS of these latent factors identified 74 genome-wide significant (GWS) loci (p<5×10 −8 ), including many previously implicated in neuroimaging phenotypes, behavioral traits, and psychiatric conditions. LDSC of latent factor GWAS results found that SA-derived factors had a positive genetic correlation with bipolar disorder (BPD), and major depressive disorder (MDD), and a negative genetic correlation with attention deficit hyperactivity disorder (ADHD), MDD, and insomnia, while CT factors displayed a negative genetic correlation with alcohol dependence. Jointly modeling the genetic architecture of complex traits and investigating multivariate genetic links across phenotypes offers a new vantage point for mapping genetically informed cortical networks. HIGHLIGHTS Genomic SEM can examine genetic correlation across cortical regions. We inferred regional genetic networks of cortical thickness and surface area. Network-associated variants have been implicated in multiple traits. These networks are genetically correlated with several psychiatric disorders including MDD, bipolar, ADHD, and alcohol dependence.