Recent discovery of approximately 270 common genetic variants associated with schizophrenia has enabled polygenic risk scores (PRS) to be measured in the population. We hypothesized that normal variation in PRS would be associated with magnetic resonance imaging (MRI) phenotypes of brain morphometry and tissue composition. We used the largest extant genome-wide association dataset (N = 69,369 cases and N = 236,642 healthy controls) to measure PRS for schizophrenia in a large sample of adults from the UK Biobank (Nmax = 29,878) who had multiple micro- and macrostructural MRI metrics measured at each of 180 cortical areas, seven subcortical structures, and 15 major white matter tracts. Linear mixed-effect models were used to investigate associations between PRS and brain structure at global and regional scales, controlled for multiple comparisons. Polygenic risk was significantly associated with reduced neurite density index (NDI) at global brain scale, at 149 cortical regions, five subcortical structures, and 14 white matter tracts. Other microstructural parameters, e.g., fractional anisotropy, that were correlated with NDI were also significantly associated with PRS. Genetic effects on multiple MRI phenotypes were co-located in temporal, cingulate, and prefrontal cortical areas, insula, and hippocampus. Post-hoc bidirectional Mendelian randomization analyses provided preliminary evidence in support of a causal relationship between (reduced) thalamic NDI and (increased) risk of schizophrenia. Risk-related reduction in NDI is plausibly indicative of reduced density of myelinated axons and dendritic arborization in large-scale cortico-subcortical networks. Cortical, subcortical, and white matter microstructure may be linked to the genetic mechanisms of schizophrenia.

Our understanding of the genetic architecture of the human cerebral cortex is limited both in terms of the diversity of brain structural phenotypes and the anatomical granularity of their associations with genetic variants. Here, we conducted genome-wide association meta-analysis of 13 structural and diffusion magnetic resonance imaging derived cortical phenotypes, measured globally and at 180 bilaterally averaged regions in 36,843 individuals from the UK Biobank and the ABCD cohorts. These phenotypes include cortical thickness, surface area, grey matter volume, and measures of folding, neurite density, and water diffusion. We identified 4,349 experiment-wide significant loci associated with global and regional phenotypes. Multiple lines of analyses identified four genetic latent structures and causal relationships between surface area and some measures of cortical folding. These latent structures partly relate to different underlying gene expression trajectories during development and are enriched for different cell types. We also identified differential enrichment for neurodevelopmental and constrained genes and demonstrate that common genetic variants associated with surface area and volume specifically are associated with cephalic disorders. Finally, we identified complex inter-phenotype and inter-regional genetic relationships among the 13 phenotypes which reflect developmental differences among them. These analyses help refine the role of common genetic variants in human cortical development and organisation.

Background Recent discovery of hundreds of common gene variants associated with schizophrenia has enabled polygenic risk scores (PRS) to be measured in the population. It is hypothesized that normal variation in genetic risk of schizophrenia should be associated with MRI changes in brain morphometry and tissue composition. Methods We used the largest extant genome-wide association dataset (N = 69,369 cases and N = 236,642 healthy controls) to measure PRS for schizophrenia in a large sample of adults from the UK Biobank (Nmax = 29,878) who had multiple micro- and macro-structural MRI metrics measured at each of 180 cortical areas and seven subcortical structures. Linear mixed effect models were used to investigate associations between schizophrenia PRS and brain structure at global and regional scales, controlled for multiple comparisons. Results Micro-structural phenotypes were more robustly associated with schizophrenia PRS than macro-structural phenotypes. Polygenic risk was significantly associated with reduced neurite density index (NDI) at global brain scale, at 149 cortical regions, and five subcortical structures. Other micro-structural parameters, e.g., fractional anisotropy, that were correlated with NDI were also significantly associated with schizophrenia PRS. Genetic effects on multiple MRI phenotypes were co-located in temporal, cingulate and prefrontal cortical areas, insula, and hippocampus. Conclusions We show widespread cortical and subcortical grey matter micro-structure associations with schizophrenia PRS. Across all investigated phenotypes NDI, a measure of the density of myelinated axons and dendrites, showed the most robust associations with schizophrenia PRS. We interpret these results as indicative of reduced density of myelinated axons and dendritic arborization in large-scale cortico-subcortical networks mediating the genetic risk for schizophrenia.

Recent studies suggest shared genetic effects on both schizophrenia and brain structure, but it has been challenging to specify which genes mediate this pleiotropic association. We accessed genome-wide association data on schizophrenia (N=69,369 cases; 236,642 controls), and on three magnetic resonance imaging (MRI) metrics (surface area, cortical thickness, neurite density index) measured at 180 cortical areas (N=36,843). Using Hi-C-coupled MAGMA, we identified 61 genes that were significantly associated with both schizophrenia and one or more MRI metrics. Whole genome analysis demonstrated significant genetic covariation between schizophrenia and area or thickness of most cortical regions. Genetic similarity between cortical areas was strongly coupled to covariance of their MRI metrics, and genetic covariation between schizophrenia and cortical regional phenotypes was greatest in the hubs of the corresponding structural covariance network. Three genomic regions, on chromosomes 3p21, 17q21 and 11p11, were enriched for neurodevelopmental processes and consistently implicated in these pleiotropic associations between schizophrenia and cortical network organization.

No abstract