Large-scale GWAS reveals genetic architecture of brain white matter microstructure and genetic overlap with cognitive and mental health traits (n = 17,706)

Principled distillation of UK Biobank phenotype data reveals underlying structure in human variation

Data within biobanks capture broad yet detailed indices of human variation, but biobank-wide insights can be difficult to extract due to complexity and scale. Here, using large-scale factor analysis, we distill hundreds of variables (diagnoses, assessments and survey items) into 35 latent constructs, using data from unrelated individuals with predominantly estimated European genetic ancestry in UK Biobank. These factors recapitulate known disease classifications, disentangle elements of socioeconomic status, highlight the relevance of psychiatric constructs to health and improve measurement of pro-health behaviours. We go on to demonstrate the power of this approach to clarify genetic signal, enhance discovery and identify associations between underlying phenotypic structure and health outcomes. In building a deeper understanding of ways in which constructs such as socioeconomic status, trauma, or physical activity are structured in the dataset, we emphasize the importance of considering the interwoven nature of the human phenome when evaluating public health patterns.

A multivariate to multivariate approach for voxel-wise genome-wide association analysis

The joint analysis of imaging-genetics data facilitates the systematic investigation of genetic effects on brain structures and functions with spatial specificity. We focus on voxel-wise genome-wide association analysis, which may involve trillions of single nucleotide polymorphism (SNP)-voxel pairs. We attempt to identify underlying organized association patterns of SNP-voxel pairs and understand the polygenic and pleiotropic networks on brain imaging traits. We propose a bi-clique graph structure (ie, a set of SNPs highly correlated with a cluster of voxels) for the systematic association pattern. Next, we develop computational strategies to detect latent SNP-voxel bi-cliques and an inference model for statistical testing. We further provide theoretical results to guarantee the accuracy of our computational algorithms and statistical inference. We validate our method by extensive simulation studies, and then apply it to the whole genome genetic and voxel-level white matter integrity data collected from 1052 participants of the human connectome project. The results demonstrate multiple genetic loci influencing white matter integrity measures on splenium and genu of the corpus callosum.

Multi-organ imaging-derived polygenic indexes for brain and body health

The UK Biobank (UKB) imaging project is a crucial resource for biomedical research, but is limited to 100,000 participants due to cost and accessibility barriers. Here we used genetic data to predict heritable imaging-derived phenotypes (IDPs) for a larger cohort. We developed and evaluated 4,375 IDP genetic scores (IGS) derived from UKB brain and body images. When applied to UKB participants who were not imaged, IGS revealed links to numerous phenotypes and stratified participants at increased risk for both brain and somatic diseases. For example, IGS identified individuals at higher risk for Alzheimer's disease and multiple sclerosis, offering additional insights beyond traditional polygenic risk scores of these diseases. When applied to independent external cohorts, IGS also stratified those at high disease risk in the All of Us Research Program and the Alzheimer's Disease Neuroimaging Initiative study. Our results demonstrate that, while the UKB imaging cohort is largely healthy and may not be the most enriched for disease risk management, it holds immense potential for stratifying the risk of various brain and body diseases in broader external genetic cohorts.

An adaptive and robust method for multi-trait analysis of genome-wide association studies using summary statistics

Genome-wide association studies (GWAS) have identified thousands of genetic variants associated with human traits or diseases in the past decade. Nevertheless, much of the heritability of many traits is still unaccounted for. Commonly used single-trait analysis methods are conservative, while multi-trait methods improve statistical power by integrating association evidence across multiple traits. In contrast to individual-level data, GWAS summary statistics are usually publicly available, and thus methods using only summary statistics have greater usage. Although many methods have been developed for joint analysis of multiple traits using summary statistics, there are many issues, including inconsistent performance, computational inefficiency, and numerical problems when considering lots of traits. To address these challenges, we propose a multi-trait adaptive Fisher method for summary statistics (MTAFS), a computationally efficient method with robust power performance. We applied MTAFS to two sets of brain imaging derived phenotypes (IDPs) from the UK Biobank, including a set of 58 Volumetric IDPs and a set of 212 Area IDPs. Through annotation analysis, the underlying genes of the SNPs identified by MTAFS were found to exhibit higher expression and are significantly enriched in brain-related tissues. Together with results from a simulation study, MTAFS shows its advantage over existing multi-trait methods, with robust performance across a range of underlying settings. It controls type 1 error well and can efficiently handle a large number of traits.

Diffusion imaging genomics provides novel insight into early mechanisms of cerebral small vessel disease

Cerebral small vessel disease (cSVD) is a leading cause of stroke and dementia. Genetic risk loci for white matter hyperintensities (WMH), the most common MRI-marker of cSVD in older age, were recently shown to be significantly associated with white matter (WM) microstructure on diffusion tensor imaging (signal-based) in young adults. To provide new insights into these early changes in WM microstructure and their relation with cSVD, we sought to explore the genetic underpinnings of cutting-edge tissue-based diffusion imaging markers across the adult lifespan. We conducted a genome-wide association study of neurite orientation dispersion and density imaging (NODDI) markers in young adults (i-Share study: N = 1 758, (mean[range]) 22.1[18-35] years), with follow-up in young middle-aged (Rhineland Study: N = 714, 35.2[30-40] years) and late middle-aged to older individuals (UK Biobank: N = 33 224, 64.3[45-82] years). We identified 21 loci associated with NODDI markers across brain regions in young adults. The most robust association, replicated in both follow-up cohorts, was with Neurite Density Index (NDI) at chr5q14.3, a known WMH locus in VCAN. Two additional loci were replicated in UK Biobank, at chr17q21.2 with NDI, and chr19q13.12 with Orientation Dispersion Index (ODI). Transcriptome-wide association studies showed associations of STAT3 expression in arterial and adipose tissue (chr17q21.2) with NDI, and of several genes at chr19q13.12 with ODI. Genetic susceptibility to larger WMH volume, but not to vascular risk factors, was significantly associated with decreased NDI in young adults, especially in regions known to harbor WMH in older age. Individually, seven of 25 known WMH risk loci were associated with NDI in young adults. In conclusion, we identified multiple novel genetic risk loci associated with NODDI markers, particularly NDI, in early adulthood. These point to possible early-life mechanisms underlying cSVD and to processes involving remyelination, neurodevelopment and neurodegeneration, with a potential for novel approaches to prevention.

The mechanisms of white matter injury and immune system crosstalk in promoting the progression of Parkinson's disease: a narrative review

Parkinson's disease (PD) is neurodegenerative disease in middle-aged and elderly people with some pathological mechanisms including immune disorder, neuroinflammation, white matter injury and abnormal aggregation of alpha-synuclein, etc. New research suggests that white matter injury may be important in the development of PD, but how inflammation, the immune system, and white matter damage interact to harm dopamine neurons is not yet understood. Therefore, it is particularly important to delve into the crosstalk between immune cells in the central and peripheral nervous system based on the study of white matter damage in PD. This crosstalk could not only exacerbate the pathological process of PD but may also reveal new therapeutic targets. By understanding how immune cells penetrate through the blood-brain barrier and activate inflammatory responses within the central nervous system, we can better grasp the impact of structural destruction of white matter in PD and explore how this process can be modulated to mitigate or combat disease progression. Microglia, astrocytes, oligodendrocytes and peripheral immune cells (especially T cells) play a central role in its pathological process where these immune cells produce and respond to pro-inflammatory cytokines such as tumor necrosis factor (TNF-α), interleukin-1β(IL-1β) and interleukin-6(IL-6), and white matter injury causes microglia to become pro-inflammatory and release inflammatory mediators, which attract more immune cells to the damaged area, increasing the inflammatory response. Moreover, white matter damage also causes dysfunction of blood-brain barrier, allows peripheral immune cells and inflammatory factors to invade the brain further, and enhances microglia activation forming a vicious circle that intensifies neuroinflammation. And these factors collectively promote the neuroinflammatory environment and neurodegeneration changes of PD. Overall, these findings not only deepen our understanding of the complexity of PD, but also provide new targets for the development of therapeutic strategies focused on inflammation and immune regulation mechanisms. In summary, this review provided the theoretical basis for clarifying the pathogenesis of PD, summarized the association between white matter damage and the immune cells in the central and peripheral nervous systems, and then emphasized their potential specific mechanisms of achieving crosstalk with further aggravating the pathological process of PD.

Polygenic architecture of brain structure and function, behaviors, and psychopathologies in children

The human brain undergoes structural and functional changes during childhood, a critical period in cognitive and behavioral development. Understanding the genetic architecture of the brain development in children can offer valuable insights into the development of the brain, cognition, and behaviors. Here, we integrated brain imaging-genetic-phenotype data from over 8,600 preadolescent children of diverse ethnic backgrounds using multivariate statistical techniques. We found a low-to-moderate level of SNP-based heritability in most IDPs, which is lower compared to the adult brain. Using sparse generalized canonical correlation analysis (SGCCA), we identified several covariation patterns among genome-wide polygenic scores (GPSs) of 29 traits, 7 different modalities of brain imaging-derived phenotypes (IDPs), and 266 cognitive and psychological phenotype data. In structural MRI, significant positive associations were observed between total grey matter volume, left ventral diencephalon volume, surface area of right accumbens and the GPSs of cognition-related traits. Conversely, negative associations were found with the GPSs of ADHD, depression and neuroticism. Additionally, we identified a significant positive association between educational attainment GPS and regional brain activation during the N-back task. The BMI GPS showed a positive association with fractional anisotropy (FA) of connectivity between the cerebellum cortex and amygdala in diffusion MRI, while the GPSs for educational attainment and cannabis use were negatively associated with the same IDPs. Our GPS-based prediction models revealed substantial genetic contributions to cognitive variability, while the genetic basis for many mental and behavioral phenotypes remained elusive. This study delivers a comprehensive map of the relationships between genetic profiles, neuroanatomical diversity, and the spectrum of cognitive and behavioral traits in preadolescence.

Genetic and Clinical Correlates of AI-Based Brain Aging Patterns in Cognitively Unimpaired Individuals

Importance

Brain aging elicits complex neuroanatomical changes influenced by multiple age-related pathologies. Understanding the heterogeneity of structural brain changes in aging may provide insights into preclinical stages of neurodegenerative diseases.

Objective

To derive subgroups with common patterns of variation in participants without diagnosed cognitive impairment (WODCI) in a data-driven manner and relate them to genetics, biomedical measures, and cognitive decline trajectories.

Design, Setting, and Participants

Data acquisition for this cohort study was performed from 1999 to 2020. Data consolidation and harmonization were conducted from July 2017 to July 2021. Age-specific subgroups of structural brain measures were modeled in 4 decade-long intervals spanning ages 45 to 85 years using a deep learning, semisupervised clustering method leveraging generative adversarial networks. Data were analyzed from July 2021 to February 2023 and were drawn from the Imaging-Based Coordinate System for Aging and Neurodegenerative Diseases (iSTAGING) international consortium. Individuals WODCI at baseline spanning ages 45 to 85 years were included, with greater than 50 000 data time points.

Exposures

Individuals WODCI at baseline scan.

Main Outcomes and Measures

Three subgroups, consistent across decades, were identified within the WODCI population. Associations with genetics, cardiovascular risk factors (CVRFs), amyloid β (Aβ), and future cognitive decline were assessed.

Results

In a sample of 27 402 individuals (mean [SD] age, 63.0 [8.3] years; 15 146 female [55%]) WODCI, 3 subgroups were identified in contrast with the reference group: a typical aging subgroup, A1, with a specific pattern of modest atrophy and white matter hyperintensity (WMH) load, and 2 accelerated aging subgroups, A2 and A3, with characteristics that were more distinct at age 65 years and older. A2 was associated with hypertension, WMH, and vascular disease-related genetic variants and was enriched for Aβ positivity (ages ≥65 years) and apolipoprotein E (APOE) ε4 carriers. A3 showed severe, widespread atrophy, moderate presence of CVRFs, and greater cognitive decline. Genetic variants associated with A1 were protective for WMH (rs7209235: mean [SD] B = -0.07 [0.01]; P value = 2.31 × 10-9) and Alzheimer disease (rs72932727: mean [SD] B = 0.1 [0.02]; P value = 6.49 × 10-9), whereas the converse was observed for A2 (rs7209235: mean [SD] B = 0.1 [0.01]; P value = 1.73 × 10-15 and rs72932727: mean [SD] B = -0.09 [0.02]; P value = 4.05 × 10-7, respectively); variants in A3 were associated with regional atrophy (rs167684: mean [SD] B = 0.08 [0.01]; P value = 7.22 × 10-12) and white matter integrity measures (rs1636250: mean [SD] B = 0.06 [0.01]; P value = 4.90 × 10-7).

Conclusions and Relevance

The 3 subgroups showed distinct associations with CVRFs, genetics, and subsequent cognitive decline. These subgroups likely reflect multiple underlying neuropathologic processes and affect susceptibility to Alzheimer disease, paving pathways toward patient stratification at early asymptomatic stages and promoting precision medicine in clinical trials and health care.

Unsupervised deep representation learning enables phenotype discovery for genetic association studies of brain imaging

Understanding the genetic architecture of brain structure is challenging, partly due to difficulties in designing robust, non-biased descriptors of brain morphology. Until recently, brain measures for genome-wide association studies (GWAS) consisted of traditionally expert-defined or software-derived image-derived phenotypes (IDPs) that are often based on theoretical preconceptions or computed from limited amounts of data. Here, we present an approach to derive brain imaging phenotypes using unsupervised deep representation learning. We train a 3-D convolutional autoencoder model with reconstruction loss on 6130 UK Biobank (UKBB) participants' T1 or T2-FLAIR (T2) brain MRIs to create a 128-dimensional representation known as Unsupervised Deep learning derived Imaging Phenotypes (UDIPs). GWAS of these UDIPs in held-out UKBB subjects (n = 22,880 discovery and n = 12,359/11,265 replication cohorts for T1/T2) identified 9457 significant SNPs organized into 97 independent genetic loci of which 60 loci were replicated. Twenty-six loci were not reported in earlier T1 and T2 IDP-based UK Biobank GWAS. We developed a perturbation-based decoder interpretation approach to show that these loci are associated with UDIPs mapped to multiple relevant brain regions. Our results established unsupervised deep learning can derive robust, unbiased, heritable, and interpretable brain imaging phenotypes.

Genetic architecture of the structural connectome

Myelinated axons form long-range connections that enable rapid communication between distant brain regions, but how genetics governs the strength and organization of these connections remains unclear. We perform genome-wide association studies of 206 structural connectivity measures derived from diffusion magnetic resonance imaging tractography of 26,333 UK Biobank participants, each representing the density of myelinated connections within or between a pair of cortical networks, subcortical structures or cortical hemispheres. We identify 30 independent genome-wide significant variants after Bonferroni correction for the number of measures studied (126 variants at nominal genome-wide significance) implicating genes involved in myelination (SEMA3A), neurite elongation and guidance (NUAK1, STRN, DPYSL2, EPHA3, SEMA3A, HGF, SHTN1), neural cell proliferation and differentiation (GMNC, CELF4, HGF), neuronal migration (CCDC88C), cytoskeletal organization (CTTNBP2, MAPT, DAAM1, MYO16, PLEC), and brain metal transport (SLC39A8). These variants have four broad patterns of spatial association with structural connectivity: some have disproportionately strong associations with corticothalamic connectivity, interhemispheric connectivity, or both, while others are more spatially diffuse. Structural connectivity measures are highly polygenic, with a median of 9.1 percent of common variants estimated to have non-zero effects on each measure, and exhibited signatures of negative selection. Structural connectivity measures have significant genetic correlations with a variety of neuropsychiatric and cognitive traits, indicating that connectivity-altering variants tend to influence brain health and cognitive function. Heritability is enriched in regions with increased chromatin accessibility in adult oligodendrocytes (as well as microglia, inhibitory neurons and astrocytes) and multiple fetal cell types, suggesting that genetic control of structural connectivity is partially mediated by effects on myelination and early brain development. Our results indicate pervasive, pleiotropic, and spatially structured genetic control of white-matter structural connectivity via diverse neurodevelopmental pathways, and support the relevance of this genetic control to healthy brain function.

Neuroimaging genomics as a window into the evolution of human sulcal organization

Primate brain evolution has involved prominent expansions of the cerebral cortex, with largest effects observed in the human lineage. Such expansions were accompanied by fine-grained anatomical alterations, including increased cortical folding. However, the molecular bases of evolutionary alterations in human sulcal organization are not yet well understood. Here, we integrated data from recently completed large-scale neuroimaging genetic analyses with annotations of the human genome relevant to various periods and events in our evolutionary history. These analyses identified single-nucleotide polymorphism (SNP) heritability enrichments in fetal brain human-gained enhancer (HGE) elements for a number of sulcal structures, including the central sulcus, which is implicated in human hand dexterity. We zeroed in on a genomic region that harbors DNA variants associated with left central sulcus shape, an HGE element, and genetic loci involved in neurogenesis including ZIC4, to illustrate the value of this approach for probing the complex factors contributing to human sulcal evolution.

Transcriptomic Similarity Informs Neuromorphic Deviations in Depression Biotypes

BACKGROUND

Major depressive disorder (MDD) is complicated by population heterogeneity, motivating the investigation of biotypes through imaging-derived phenotypes. However, neuromorphic heterogeneity in MDD remains unclear, and how the correlated gene expression (CGE) connectome constrains these neuromorphic anomalies in MDD biotypes has not yet been studied.

METHODS

Here, we related cortical thickness deviations in MDD biotypes to a pattern of CGE connectome. Cortical thickness was estimated from 3-dimensional T1-weighted magnetic resonance images in 2 independent cohorts (discovery cohort: N = 425; replication cohort: N = 217). The transcriptional activity was measured according to Allen Human Brain Atlas. A density peak-based clustering algorithm was used to identify MDD biotypes.

RESULTS

We found that patients with MDD were clustered into 2 replicated biotypes based on single-patient regional deviations from healthy control participants across 2 datasets. Biotype 1 mainly exhibited cortical thinning across the brain, whereas biotype 2 mainly showed cortical thickening in the brain. Using brainwide gene expression data, we found that deviations of transcriptionally connected neighbors predicted regional deviation for both biotypes. Furthermore, putative CGE-informed epicenters of biotype 1 were concentrated on the cognitive control circuit, whereas biotype 2 epicenters were located in the social perception circuit. The patterns of epicenter likelihood were separately associated with depression- and anxiety-response maps, suggesting that epicenters of MDD biotypes may be associated with clinical efficacies.

CONCLUSIONS

Our findings linked the CGE connectome and neuromorphic deviations to identify distinct epicenters in MDD biotypes, providing insight into how microscale gene expressions informed MDD biotypes.

Genome-Wide Association Studies of 3 Distinct Recovery Phenotypes in Mild Ischemic Stroke

BACKGROUND AND OBJECTIVES

Stroke genetic research has made substantial progress in the past decade. Its recovery application, however, remains behind, in part due to its reliance on the modified Rankin Scale (mRS) score as a measure of poststroke outcome. The mRS does not map well to biological processes because numerous psychosocial factors drive much of what the mRS captures. Second, the mRS contains multiple disparate biological events into a single measure further limiting its use for biological discovery. This led us to investigate the effect of distinct stroke recovery phenotypes on genetic variation associations with Genome-Wide Association Studies (GWASs) by repurposing the NIH Stroke Scale (NIHSS) and its subscores.

METHODS

In the Vitamin Intervention for Stroke Prevention cohort, we estimated changes in cognition, motor, and global impairments over 2 years using specific measures. We included genotyped participants with a total NIHSS score greater than zero at randomization and excluded those with recurrent stroke during the trial. A GWAS linear mixed-effects model predicted score changes, with participant as a random effect, and included initial score, age, sex, treatment group, and the first 5 ancestry principal components.

RESULTS

In total, 1,270 participants (64% male) were included with a median NIHSS score of 2 (interquartile range [IQR] 1-3) and median age 68 (IQR 59-75) years. At randomization, 20% had cognitive deficits (NIHSS Cog-4 score >0) and 70% had ≥1 motor deficits (impairment score >1). At 2 years, these percentages improved to 7.2% with cognitive deficits and 30% with motor deficits. GWAS identified novel suggestive gene-impairment associations (p < 5e-6) for cognition (CAMK2D, EVX2, LINC0143, PTPRM, SGMS1, and SMAD2), motor (ACBD6, KDM4B, MARK4, PTPRS, ROBO1, and ROBO2), and global (MSR1 and ROBO2) impairments.

DISCUSSION

Defining domain-specific stroke recovery phenotypes and using longitudinal clinical trial designs can help detect novel genes associated with chronic recovery. These data support the use of granular endpoints to identify genetic associations related to stroke recovery.

Identification and validation of supervariants reveal novel loci associated with human white matter microstructure

As an essential part of the central nervous system, white matter coordinates communications between different brain regions and is related to a wide range of neurodegenerative and neuropsychiatric disorders. Previous genome-wide association studies (GWASs) have uncovered loci associated with white matter microstructure. However, GWASs suffer from limited reproducibility and difficulties in detecting multi-single-nucleotide polymorphism (multi-SNP) and epistatic effects. In this study, we adopt the concept of supervariants, a combination of alleles in multiple loci, to account for potential multi-SNP effects. We perform supervariant identification and validation to identify loci associated with 22 white matter fractional anisotropy phenotypes derived from diffusion tensor imaging. To increase reproducibility, we use United Kingdom (UK) Biobank White British (n = 30,842) data for discovery and internal validation, and UK Biobank White but non-British (n = 1927) data, Europeans from the Adolescent Brain Cognitive Development study (n = 4399) data, and Europeans from the Human Connectome Project (n = 319) data for external validation. We identify 23 novel loci on the discovery set that have not been reported in the previous GWASs on white matter microstructure. Among them, three supervariants on genomic regions 5q35.1, 8p21.2, and 19q13.32 have P-values lower than 0.05 in the meta-analysis of the three independent validation data sets. These supervariants contain genetic variants located in genes that have been related to brain structures, cognitive functions, and neuropsychiatric diseases. Our findings provide a better understanding of the genetic architecture underlying white matter microstructure.

Vascular risk factors, white matter microstructure, and depressive symptoms: a longitudinal analysis in the UK Biobank

BACKGROUND

Cumulative burden from vascular risk factors (VRFs) has been associated with an increased risk of depressive symptoms in mid- and later life. It has been hypothesised that this association arises because VRFs disconnect fronto-subcortical white matter tracts involved in mood regulation, which puts older adults at higher risk of developing depressive symptoms. However, evidence for the hypothesis that disconnection of white matter tracts underlies the association between VRF burden and depressive symptoms from longitudinal studies is scarce.

METHODS

This preregistered study analysed longitudinal data from 6,964 middle-aged and older adults from the UK Biobank who participated in consecutive assessments of VRFs, brain imaging, and depressive symptoms. Using mediation modelling, we directly tested to what extend white matter microstructure mediates the longitudinal association between VRF burden and depressive symptoms.

RESULTS

VRF burden showed a small association with depressive symptoms at follow-up. However, there was no evidence that fractional anisotropy (FA) of white matter tracts mediated this association. Additional analyses also yielded no mediating effects using alternative operationalisations of VRF burden, mean diffusivity (MD) of single tracts, or overall average of tract-based white matter microstructure (global FA, global MD, white matter hyperintensity volume).

CONCLUSIONS

Our results lend no support to the hypothesis that disconnection of white matter tracts underlies the association between VRF burden and depressive symptoms, while highlighting the relevance of using longitudinal data to directly test pathways linking vascular and mental health.

Longitudinal development of the human white matter structural connectome and its association with brain transcriptomic and cellular architecture

From childhood to adolescence, the spatiotemporal development pattern of the human brain white matter connectome and its underlying transcriptomic and cellular mechanisms remain largely unknown. With a longitudinal diffusion MRI cohort of 604 participants, we map the developmental trajectory of the white matter connectome from global to regional levels and identify that most brain network properties followed a linear developmental trajectory. Importantly, connectome-transcriptomic analysis reveals that the spatial development pattern of white matter connectome is potentially regulated by the transcriptomic architecture, with positively correlated genes involve in ion transport- and development-related pathways expressed in excitatory and inhibitory neurons, and negatively correlated genes enriches in synapse- and development-related pathways expressed in astrocytes, inhibitory neurons and microglia. Additionally, the macroscale developmental pattern is also associated with myelin content and thicknesses of specific laminas. These findings offer insights into the underlying genetics and neural mechanisms of macroscale white matter connectome development from childhood to adolescence.

An augmented Mendelian randomization approach provides causality of brain imaging features on complex traits in a single biobank-scale dataset

Mendelian randomization (MR) is an effective approach for revealing causal risk factors that underpin complex traits and diseases. While MR has been more widely applied under two-sample settings, it is more promising to be used in one single large cohort given the rise of biobank-scale datasets that simultaneously contain genotype data, brain imaging data, and matched complex traits from the same individual. However, most existing multivariable MR methods have been developed for two-sample setting or a small number of exposures. In this study, we introduce a one-sample multivariable MR method based on partial least squares and Lasso regression (MR-PL). MR-PL is capable of considering the correlation among exposures (e.g., brain imaging features) when the number of exposures is extremely upscaled, while also correcting for winner's curse bias. We performed extensive and systematic simulations, and demonstrated the robustness and reliability of our method. Comprehensive simulations confirmed that MR-PL can generate more precise causal estimates with lower false positive rates than alternative approaches. Finally, we applied MR-PL to the datasets from UK Biobank to reveal the causal effects of 36 white matter tracts on 180 complex traits, and showed putative white matter tracts that are implicated in smoking, blood vascular function-related traits, and eating behaviors.

The genetic relationships between brain structure and schizophrenia

Genetic risks for schizophrenia are theoretically mediated by genetic effects on brain structure but it has been unclear which genes are associated with both schizophrenia and cortical phenotypes. We accessed genome-wide association studies (GWAS) of schizophrenia (N = 69,369 cases; 236,642 controls), and of three magnetic resonance imaging (MRI) metrics (surface area, cortical thickness, neurite density index) measured at 180 cortical areas (N = 36,843, UK Biobank). Using Hi-C-coupled MAGMA, 61 genes were significantly associated with both schizophrenia and one or more MRI metrics. Whole genome analysis with partial least squares demonstrated significant genetic covariation between schizophrenia and area or thickness of most cortical regions. Genetic similarity between cortical areas was strongly coupled to their phenotypic covariance, and genetic covariation between schizophrenia and brain phenotypes was strongest in the hubs of structural covariance networks. Pleiotropically associated genes were enriched for neurodevelopmental processes and positionally concentrated in chromosomes 3p21, 17q21 and 11p11. Mendelian randomization analysis indicated that genetically determined variation in a posterior cingulate cortical area could be causal for schizophrenia. Parallel analyses of GWAS on bipolar disorder, Alzheimer's disease and height showed that pleiotropic association with MRI metrics was stronger for schizophrenia compared to other disorders.

Psychiatric disorders and brain white matter exhibit genetic overlap implicating developmental and neural cell biology

Improved understanding of the shared genetic architecture between psychiatric disorders and brain white matter may provide mechanistic insights for observed phenotypic associations. Our objective is to characterize the shared genetic architecture of bipolar disorder (BD), major depression (MD), and schizophrenia (SZ) with white matter fractional anisotropy (FA) and identify shared genetic loci to uncover biological underpinnings. We used genome-wide association study (GWAS) summary statistics for BD (n = 413,466), MD (n = 420,359), SZ (n = 320,404), and white matter FA (n = 33,292) to uncover the genetic architecture (i.e., polygenicity and discoverability) of each phenotype and their genetic overlap (i.e., genetic correlations, overlapping trait-influencing variants, and shared loci). This revealed that BD, MD, and SZ are at least 7-times more polygenic and less genetically discoverable than average FA. Even in the presence of weak genetic correlations (range = -0.05 to -0.09), average FA shared an estimated 42.5%, 43.0%, and 90.7% of trait-influencing variants as well as 12, 4, and 28 shared loci with BD, MD, and SZ, respectively. Shared variants were mapped to genes and tested for enrichment among gene-sets which implicated neurodevelopmental expression, neural cell types, myelin, and cell adhesion molecules. For BD and SZ, case vs control tract-level differences in FA associated with genetic correlations between those same tracts and the respective disorder (rBD = 0.83, p = 4.99e-7 and rSZ = 0.65, p = 5.79e-4). Genetic overlap at the tract-level was consistent with average FA results. Overall, these findings suggest a genetic basis for the involvement of brain white matter aberrations in the pathophysiology of psychiatric disorders.

Genome-wide association study of hippocampal blood-oxygen-level-dependent-cerebral blood flow correlation in Chinese Han population

Correlation between blood-oxygen-level-dependent (BOLD) and cerebral blood flow (CBF) has been used as an index of neurovascular coupling. Hippocampal BOLD-CBF correlation is associated with neurocognition, and the reduced correlation is associated with neuropsychiatric disorders. We conducted the first genome-wide association study of the hippocampal BOLD-CBF correlation in 4,832 Chinese Han subjects. The hippocampal BOLD-CBF correlation had an estimated heritability of 16.2-23.9% and showed reliable genome-wide significant association with a locus at 3q28, in which many variants have been linked to neuroimaging and cerebrospinal fluid markers of Alzheimer's disease. Gene-based association analyses showed four significant genes (GMNC, CRTC2, DENND4B, and GATAD2B) and revealed enrichment for mast cell calcium mobilization, microglial cell proliferation, and ubiquitin-related proteolysis pathways that regulate different cellular components of the neurovascular unit. This is the first unbiased identification of the association of hippocampal BOLD-CBF correlation, providing fresh insights into the genetic architecture of hippocampal neurovascular coupling.

Assessing efficiency of fine-mapping obesity-associated variants through leveraging ancestry architecture and functional annotation using PAGE and UKBB cohorts

Inadequate representation of non-European ancestry populations in genome-wide association studies (GWAS) has limited opportunities to isolate functional variants. Fine-mapping in multi-ancestry populations should improve the efficiency of prioritizing variants for functional interrogation. To evaluate this hypothesis, we leveraged ancestry architecture to perform comparative GWAS and fine-mapping of obesity-related phenotypes in European ancestry populations from the UK Biobank (UKBB) and multi-ancestry samples from the Population Architecture for Genetic Epidemiology (PAGE) consortium with comparable sample sizes. In the investigated regions with genome-wide significant associations for obesity-related traits, fine-mapping in our ancestrally diverse sample led to 95% and 99% credible sets (CS) with fewer variants than in the European ancestry sample. Lead fine-mapped variants in PAGE regions had higher average coding scores, and higher average posterior probabilities for causality compared to UKBB. Importantly, 99% CS in PAGE loci contained strong expression quantitative trait loci (eQTLs) in adipose tissues or harbored more variants in tighter linkage disequilibrium (LD) with eQTLs. Leveraging ancestrally diverse populations with heterogeneous ancestry architectures, coupled with functional annotation, increased fine-mapping efficiency and performance, and reduced the set of candidate variants for consideration for future functional studies. Significant overlap in genetic causal variants across populations suggests generalizability of genetic mechanisms underpinning obesity-related traits across populations.

Deciphering the causal relationship between blood pressure and regional white matter integrity: A two-sample Mendelian randomization study

Elevated arterial blood pressure (BP) is a common risk factor for cerebrovascular and cardiovascular diseases, but no causal relationship has been established between BP and cerebral white matter (WM) integrity. In this study, we performed a two-sample Mendelian randomization (MR) analysis with individual-level data by defining two nonoverlapping sets of European ancestry individuals (genetics-exposure set: N = 203,111; mean age = 56.71 years, genetics-outcome set: N = 16,156; mean age = 54.61 years) from UK Biobank to evaluate the causal effects of BP on regional WM integrity, measured by fractional anisotropy of diffusion tensor imaging. Two BP traits: systolic and diastolic blood pressure were used as exposures. Genetic variant was carefully selected as instrumental variable (IV) under the MR analysis assumptions. We existing large-scale genome-wide association study summary data for validation. The main method used was a generalized version of inverse-variance weight method while other MR methods were also applied for consistent findings. Two additional MR analyses were performed to exclude the possibility of reverse causality. We found significantly negative causal effects (FDR-adjusted p < .05; every 10 mmHg increase in BP leads to a decrease in FA value by .4% ~ 2%) of BP traits on a union set of 17 WM tracts, including brain regions related to cognitive function and memory. Our study extended the previous findings of association to causation for regional WM integrity, providing insights into the pathological processes of elevated BP that might chronically alter the brain microstructure in different regions.

Structural neuroimaging phenotypes and associated molecular and genomic underpinnings in autism: a review

Autism has been associated with differences in the developmental trajectories of multiple neuroanatomical features, including cortical thickness, surface area, cortical volume, measures of gyrification, and the gray-white matter tissue contrast. These neuroimaging features have been proposed as intermediate phenotypes on the gradient from genomic variation to behavioral symptoms. Hence, examining what these proxy markers represent, i.e., disentangling their associated molecular and genomic underpinnings, could provide crucial insights into the etiology and pathophysiology of autism. In line with this, an increasing number of studies are exploring the association between neuroanatomical, cellular/molecular, and (epi)genetic variation in autism, both indirectly and directly in vivo and across age. In this review, we aim to summarize the existing literature in autism (and neurotypicals) to chart a putative pathway from (i) imaging-derived neuroanatomical cortical phenotypes to (ii) underlying (neuropathological) biological processes, and (iii) associated genomic variation.

Genetic analyses identify brain structures related to cognitive impairment associated with elevated blood pressure

BACKGROUND AND AIMS

Observational studies have linked elevated blood pressure (BP) to impaired cognitive function. However, the functional and structural changes in the brain that mediate the relationship between BP elevation and cognitive impairment remain unknown. Using observational and genetic data from large consortia, this study aimed to identify brain structures potentially associated with BP values and cognitive function.

METHODS AND RESULTS

Data on BP were integrated with 3935 brain magnetic resonance imaging-derived phenotypes (IDPs) and cognitive function defined by fluid intelligence score. Observational analyses were performed in the UK Biobank and a prospective validation cohort. Mendelian randomisation (MR) analyses used genetic data derived from the UK Biobank, International Consortium for Blood Pressure, and COGENT consortium. Mendelian randomisation analysis identified a potentially adverse causal effect of higher systolic BP on cognitive function [-0.044 standard deviation (SD); 95% confidence interval (CI) -0.066, -0.021] with the MR estimate strengthening (-0.087 SD; 95% CI -0.132, -0.042), when further adjusted for diastolic BP. Mendelian randomisation analysis found 242, 168, and 68 IDPs showing significant (false discovery rate P < 0.05) association with systolic BP, diastolic BP, and pulse pressure, respectively. Most of these IDPs were inversely associated with cognitive function in observational analysis in the UK Biobank and showed concordant effects in the validation cohort. Mendelian randomisation analysis identified relationships between cognitive function and the nine of the systolic BP-associated IDPs, including the anterior thalamic radiation, anterior corona radiata, or external capsule.

CONCLUSION

Complementary MR and observational analyses identify brain structures associated with BP, which may be responsible for the adverse effects of hypertension on cognitive performance.

Bio-psycho-social factors' associations with brain age: a large-scale UK Biobank diffusion study of 35,749 participants

Brain age refers to age predicted by brain features. Brain age has previously been associated with various health and disease outcomes and suggested as a potential biomarker of general health. Few previous studies have systematically assessed brain age variability derived from single and multi-shell diffusion magnetic resonance imaging data. Here, we present multivariate models of brain age derived from various diffusion approaches and how they relate to bio-psycho-social variables within the domains of sociodemographic, cognitive, life-satisfaction, as well as health and lifestyle factors in midlife to old age (N = 35,749, 44.6-82.8 years of age). Bio-psycho-social factors could uniquely explain a small proportion of the brain age variance, in a similar pattern across diffusion approaches: cognitive scores, life satisfaction, health and lifestyle factors adding to the variance explained, but not socio-demographics. Consistent brain age associations across models were found for waist-to-hip ratio, diabetes, hypertension, smoking, matrix puzzles solving, and job and health satisfaction and perception. Furthermore, we found large variability in sex and ethnicity group differences in brain age. Our results show that brain age cannot be sufficiently explained by bio-psycho-social variables alone. However, the observed associations suggest to adjust for sex, ethnicity, cognitive factors, as well as health and lifestyle factors, and to observe bio-psycho-social factor interactions' influence on brain age in future studies.

Proteogenomic analysis of human cerebrospinal fluid identifies neurologically relevant regulation and informs causal proteins for Alzheimer's disease

The integration of quantitative trait loci (QTL) with disease genome-wide association studies (GWAS) has proven successful at prioritizing candidate genes at disease-associated loci. QTL mapping has mainly been focused on multi-tissue expression QTL or plasma protein QTL (pQTL). Here we generated the largest-to-date cerebrospinal fluid (CSF) pQTL atlas by analyzing 7,028 proteins in 3,107 samples. We identified 3,373 independent study-wide associations for 1,961 proteins, including 2,448 novel pQTLs of which 1,585 are unique to CSF, demonstrating unique genetic regulation of the CSF proteome. In addition to the established chr6p22.2-21.32 HLA region, we identified pleiotropic regions on chr3q28 near OSTN and chr19q13.32 near APOE that were enriched for neuron-specificity and neurological development. We also integrated this pQTL atlas with the latest Alzheimer's disease (AD) GWAS through PWAS, colocalization and Mendelian Randomization and identified 42 putative causal proteins for AD, 15 of which have drugs available. Finally, we developed a proteomics-based risk score for AD that outperforms genetics-based polygenic risk scores. These findings will be instrumental to further understand the biology and identify causal and druggable proteins for brain and neurological traits.

Structural architecture and brain network efficiency link polygenic scores to intelligence

Intelligence is highly heritable. Genome-wide association studies (GWAS) have shown that thousands of alleles contribute to variation in intelligence with small effect sizes. Polygenic scores (PGS), which combine these effects into one genetic summary measure, are increasingly used to investigate polygenic effects in independent samples. Whereas PGS explain a considerable amount of variance in intelligence, it is largely unknown how brain structure and function mediate this relationship. Here, we show that individuals with higher PGS for educational attainment and intelligence had higher scores on cognitive tests, larger surface area, and more efficient fiber connectivity derived by graph theory. Fiber network efficiency as well as the surface of brain areas partly located in parieto-frontal regions were found to mediate the relationship between PGS and cognitive performance. These findings are a crucial step forward in decoding the neurogenetic underpinnings of intelligence, as they identify specific regional networks that link polygenic predisposition to intelligence.

The genetic architecture of fornix white matter microstructure and their involvement in neuropsychiatric disorders

The fornix is a white matter bundle located in the center of the hippocampaldiencephalic limbic circuit that controls memory and executive functions, yet its genetic architectures and involvement in brain disorders remain largely unknown. We carried out a genome-wide association analysis of 30,832 UK Biobank individuals of the six fornix diffusion magnetic resonance imaging (dMRI) traits. The post-GWAS analysis allowed us to identify causal genetic variants in phenotypes at the single nucleotide polymorphisms (SNP), locus, and gene levels, as well as genetic overlap with brain health-related traits. We further generalized our GWAS in adolescent brain cognitive development (ABCD) cohort. The GWAS identified 63 independent significant variants within 20 genomic loci associated (P < 8.33 × 10^-9) with the six fornix dMRI traits. Geminin coiled-coil domain containing (GMNC) and NUAK family SNF1-like kinase 1 (NUAK1) gene were highlighted, which were found in UKB and replicated in ABCD. The heritability of the six traits ranged from 10% to 27%. Gene mapping strategies identified 213 genes, where 11 were supported by all of four methods. Gene-based analyses revealed pathways relating to cell development and differentiation, with astrocytes found to be significantly enriched. Pleiotropy analyses with eight neurological and psychiatric disorders revealed shared variants, especially with schizophrenia under the conjFDR threshold of 0.05. These findings advance our understanding of the complex genetic architectures of fornix and their relevance in neurological and psychiatric disorders.

The genetic architecture of the corpus callosum and its genetic overlap with common neuropsychiatric diseases

BACKGROUND

The corpus callosum (CC) is the main structure transferring information between the cerebral hemispheres. Although previous large-scale genome-wide association study (GWAS) has illustrated the genetic architecture of white matter integrity of CC, CC volume is less stressed.

METHODS

Using MRI data from 33,861 individuals in UK Biobank, we conducted univariate and multivariate GWAS for CC fractional anisotropy (FA) and volume with PLINK 2.0 and MOSTest. All discovered SNPs in the multivariate framework were functionally annotated in FUMA v1.3.8. In the meanwhile, a series of gene property analyses was conducted simultaneously. In addition, we estimated genetic relationship between CC metrics and other neuropsychiatric traits and diseases.

RESULTS

We identified a total of 36 and 82 significant genomic loci for CC FA and volume (P < 5 × 10^-8). And 53 and 27 genes were respectively mapped by four mapping strategies. For CC volume, gene-set analysis revealed pathways mainly relating to cell migration; cell-type analysis found the top enrichment in neuroglia while for CC FA in GABAergic neurons. Furthermore, we found a lot of genetic overlap and shared loci between CC FA and volume and common neuropsychiatric diseases.

DISCUSSION

Collectively, this study helps to better understand the genetic architecture of whole CC and CC subregions. However, the way to divide CC FA and volume in our study restricts the interpretations of our results. Future work will be needed to pay attention to the genetic structure of white matter volume, and an appropriate division of CC may help to better understand CC structure.

A global view of the genetic basis of Alzheimer disease

The risk of Alzheimer disease (AD) increases with age, family history and informative genetic variants. Sadly, there is still no cure or means of prevention. As in other complex diseases, uncovering genetic causes of AD could identify underlying pathological mechanisms and lead to potential treatments. Rare, autosomal dominant forms of AD occur in middle age as a result of highly penetrant genetic mutations, but the most common form of AD occurs later in life. Large-scale, genome-wide analyses indicate that 70 or more genes or loci contribute to AD. One of the major factors limiting progress is that most genetic data have been obtained from non-Hispanic white individuals in Europe and North America, preventing the development of personalized approaches to AD in individuals of other ethnicities. Fortunately, emerging genetic data from other regions - including Africa, Asia, India and South America - are now providing information on the disease from a broader range of ethnicities. Here, we summarize the current knowledge on AD genetics in populations across the world. We predominantly focus on replicated genetic discoveries but also include studies in ethnic groups where replication might not be feasible. We attempt to identify gaps that need to be addressed to achieve a complete picture of the genetic and molecular factors that drive AD in individuals across the globe.

Refining patterns of MEF2C effects in white matter microstructure and psychiatric features

Several GWAS reported Myocyte Enhancer Factor 2 C (MEF2C) gene associations with white matter microstructure and psychiatric disorders, and MEF2C involvement in pathways related to neuronal development suggests a common biological factor underlying these phenotypes. We aim to refine the MEF2C effects in the brain relying on an integrated analysis of white matter and psychiatric phenotypes in an extensively characterized sample. This study included 870 Brazilian adults (47% from an attention-deficit/hyperactivity disorder outpatient clinic) assessed through standardized psychiatric interviews, 139 of which underwent a magnetic resonance imaging scan. We evaluated variants in the MEF2C region using two approaches: 1) a gene-wide analysis, which uses the sum of polymorphism effects, and 2) SNP analyses, restricted to the independent variants within the gene. The outcomes included psychiatric phenotypes and fractional anisotropy for brain images. Results: The gene-wide analyses pointed to a nominal association between MEF2C and the Temporal Portion of the Superior Longitudinal Fasciculus (SLFTEMP). The SNP analysis identified four independent variants significantly associated with SLFTEMP and one (rs4218438) with Substance Use Disorder. Our findings showing specific associations of MEF2C variants with temporal-frontal circuitry components may help to elucidate how the MEF2C gene underlies a broad range of psychiatric phenotypes since these regions are relevant to executive and cognitive functions.

Genome-wide association study of cerebellar white matter microstructure and genetic overlap with common brain disorders

BACKGROUND

The cerebellum is recognized as being involved in neurocognitive and motor functions with communication with extra-cerebellar regions relying on the white matter integrity of the cerebellar peduncles. However, the genetic determinants of cerebellar white matter integrity remain largely unknown.

METHODS

We conducted a genome-wide association analysis of cerebellar white matter microstructure using diffusion tensor imaging data from 25,415 individuals from UK Biobank. The integrity of cerebellar white matter microstructure was measured as fractional anisotropy (FA) and mean diffusivity (MD). Identification of independent genomic loci, functional annotation, and tissue and cell-type analysis were conducted with FUMA. The linkage disequilibrium score regression (LDSC) was used to calculate genetic correlations between cerebellar white matter microstructure and regional brain volumes and brain-related traits. Furthermore, the conditional/conjunctional false discovery rate (condFDR/conjFDR) framework was employed to identify the shared genetic basis between cerebellar white matter microstructure and common brain disorders.

RESULTS

We identified 11 genetic loci (P < 8.3 × 10^-9) and 86 genes associated with cerebellar white matter microstructure. Further functional enrichment analysis implicated the involvement of GABAergic neurons and cholinergic pathways. Significant polygenetic overlap between cerebellar white matter tracts and their anatomically connected or adjacent brain regions was detected. In addition, we report the overall genetic correlation and specific loci shared between cerebellar white matter microstructural integrity and brain-related traits, including movement, cognitive, psychiatric, and cerebrovascular categories.

CONCLUSIONS

Collectively, this study represents a step forward in understanding the genetics of cerebellar white matter microstructure and its shared genetic etiology with common brain disorders.

Adolescent neurodevelopment and psychopathology: The interplay between adversity exposure and genetic risk for accelerated brain ageing

In adulthood, stress exposure and genetic risk heighten psychological vulnerability by accelerating neurobiological senescence. To investigate whether molecular and brain network maturation processes play a similar role in adolescence, we analysed genetic, as well as longitudinal task neuroimaging (inhibitory control, incentive processing) and early life adversity (i.e., material deprivation, violence) data from the Adolescent Brain and Cognitive Development study (N = 980, age range: 9-13 years). Genetic risk was estimated separately for Major Depressive Disorder (MDD) and Alzheimer's Disease (AD), two pathologies linked to stress exposure and allegedly sharing a causal connection (MDD-to-AD). Adversity and genetic risk for MDD/AD jointly predicted functional network segregation patterns suggestive of accelerated (GABA-linked) visual/attentional, but delayed (dopamine [D2]/glutamate [GLU5R]-linked) somatomotor/association system development. A positive relationship between brain maturation and psychopathology emerged only among the less vulnerable adolescents, thereby implying that normatively maladaptive neurodevelopmental alterations could foster adjustment among the more exposed and genetically more stress susceptible youths. Transcriptomic analyses suggested that sensitivity to stress may underpin the joint neurodevelopmental effect of adversity and genetic risk for MDD/AD, in line with the proposed role of negative emotionality as a precursor to AD, likely to account for the alleged causal impact of MDD on dementia onset.

Brain expression quantitative trait locus and network analyses reveal downstream effects and putative drivers for brain-related diseases

Identification of therapeutic targets from genome-wide association studies (GWAS) requires insights into downstream functional consequences. We harmonized 8,613 RNA-sequencing samples from 14 brain datasets to create the MetaBrain resource and performed cis- and trans-expression quantitative trait locus (eQTL) meta-analyses in multiple brain region- and ancestry-specific datasets (n ≤ 2,759). Many of the 16,169 cortex cis-eQTLs were tissue-dependent when compared with blood cis-eQTLs. We inferred brain cell types for 3,549 cis-eQTLs by interaction analysis. We prioritized 186 cis-eQTLs for 31 brain-related traits using Mendelian randomization and co-localization including 40 cis-eQTLs with an inferred cell type, such as a neuron-specific cis-eQTL (CYP24A1) for multiple sclerosis. We further describe 737 trans-eQTLs for 526 unique variants and 108 unique genes. We used brain-specific gene-co-regulation networks to link GWAS loci and prioritize additional genes for five central nervous system diseases. This study represents a valuable resource for post-GWAS research on central nervous system diseases.

Genetic architecture of the white matter connectome of the human brain

White matter tracts form the structural basis of large-scale brain networks. We applied brain-wide tractography to diffusion images from 30,810 adults (U.K. Biobank) and found significant heritability for 90 node-level and 851 edge-level network connectivity measures. Multivariate genome-wide association analyses identified 325 genetic loci, of which 80% had not been previously associated with brain metrics. Enrichment analyses implicated neurodevelopmental processes including neurogenesis, neural differentiation, neural migration, neural projection guidance, and axon development, as well as prenatal brain expression especially in stem cells, astrocytes, microglia, and neurons. The multivariate association profiles implicated 31 loci in connectivity between core regions of the left-hemisphere language network. Polygenic scores for psychiatric, neurological, and behavioral traits also showed significant multivariate associations with structural connectivity, each implicating distinct sets of brain regions with trait-relevant functional profiles. This large-scale mapping study revealed common genetic contributions to variation in the structural connectome of the human brain.

Epigenetic and integrative cross-omics analyses of cerebral white matter hyperintensities on MRI

Cerebral white matter hyperintensities on MRI are markers of cerebral small vessel disease, a major risk factor for dementia and stroke. Despite the successful identification of multiple genetic variants associated with this highly heritable condition, its genetic architecture remains incompletely understood. More specifically, the role of DNA methylation has received little attention. We investigated the association between white matter hyperintensity burden and DNA methylation in blood at ∼450 000 cytosine-phosphate-guanine (CpG) sites in 9732 middle-aged to older adults from 14 community-based studies. Single CpG and region-based association analyses were carried out. Functional annotation and integrative cross-omics analyses were performed to identify novel genes underlying the relationship between DNA methylation and white matter hyperintensities. We identified 12 single CpG and 46 region-based DNA methylation associations with white matter hyperintensity burden. Our top discovery single CpG, cg24202936 (P = 7.6 × 10-8), was associated with F2 expression in blood (P = 6.4 × 10-5) and co-localized with FOLH1 expression in brain (posterior probability = 0.75). Our top differentially methylated regions were in PRMT1 and in CCDC144NL-AS1, which were also represented in single CpG associations (cg17417856 and cg06809326, respectively). Through Mendelian randomization analyses cg06809326 was putatively associated with white matter hyperintensity burden (P = 0.03) and expression of CCDC144NL-AS1 possibly mediated this association. Differentially methylated region analysis, joint epigenetic association analysis and multi-omics co-localization analysis consistently identified a role of DNA methylation near SH3PXD2A, a locus previously identified in genome-wide association studies of white matter hyperintensities. Gene set enrichment analyses revealed functions of the identified DNA methylation loci in the blood-brain barrier and in the immune response. Integrative cross-omics analysis identified 19 key regulatory genes in two networks related to extracellular matrix organization, and lipid and lipoprotein metabolism. A drug-repositioning analysis indicated antihyperlipidaemic agents, more specifically peroxisome proliferator-activated receptor-alpha, as possible target drugs for white matter hyperintensities. Our epigenome-wide association study and integrative cross-omics analyses implicate novel genes influencing white matter hyperintensity burden, which converged on pathways related to the immune response and to a compromised blood-brain barrier possibly due to disrupted cell-cell and cell-extracellular matrix interactions. The results also suggest that antihyperlipidaemic therapy may contribute to lowering risk for white matter hyperintensities possibly through protection against blood-brain barrier disruption.

Genetic underpinnings of brain structural connectome for young adults

With distinct advantages in power over behavioral phenotypes, brain imaging traits have become emerging endophenotypes to dissect molecular contributions to behaviors and neuropsychiatric illnesses. Among different imaging features, brain structural connectivity (i.e., structural connectome) which summarizes the anatomical connections between different brain regions is one of the most cutting edge while under-investigated traits; and the genetic influence on the structural connectome variation remains highly elusive. Relying on a landmark imaging genetics study for young adults, we develop a biologically plausible brain network response shrinkage model to comprehensively characterize the relationship between high dimensional genetic variants and the structural connectome phenotype. Under a unified Bayesian framework, we accommodate the topology of brain network and biological architecture within the genome; and eventually establish a mechanistic mapping between genetic biomarkers and the associated brain sub-network units. An efficient expectation-maximization algorithm is developed to estimate the model and ensure computing feasibility. In the application to the Human Connectome Project Young Adult (HCP-YA) data, we establish the genetic underpinnings which are highly interpretable under functional annotation and brain tissue eQTL analysis, for the brain white matter tracts connecting the hippocampus and two cerebral hemispheres. We also show the superiority of our method in extensive simulations.

Genetic and Environmental Variation in Continuous Phenotypes in the ABCD Study®

Twin studies yield valuable insights into the sources of variation, covariation and causation in human traits. The ABCD Study® (abcdstudy.org) was designed to take advantage of four universities known for their twin research, neuroimaging, population-based sampling, and expertise in genetic epidemiology so that representative twin studies could be performed. In this paper we use the twin data to: (i) provide initial estimates of heritability for the wide range of phenotypes assessed in the ABCD Study using a consistent direct variance estimation approach, assuring that both data and methodology are sound; and (ii) provide an online resource for researchers that can serve as a reference point for future behavior genetic studies of this publicly available dataset. Data were analyzed from 772 pairs of twins aged 9-10 years at study inception, with zygosity determined using genotypic data, recruited and assessed at four twin hub sites. The online tool provides twin correlations and both standardized and unstandardized estimates of additive genetic, and environmental variation for 14,500 continuously distributed phenotypic features, including: structural and functional neuroimaging, neurocognition, personality, psychopathology, substance use propensity, physical, and environmental trait variables. The estimates were obtained using an unconstrained variance approach, so they can be incorporated directly into meta-analyses without upwardly biasing aggregate estimates. The results indicated broad consistency with prior literature where available and provided novel estimates for phenotypes without prior twin studies or those assessed at different ages. Effects of site, self-identified race/ethnicity, age and sex were statistically controlled. Results from genetic modeling of all 53,172 continuous variables, including 38,672 functional MRI variables, will be accessible via the user-friendly open-access web interface we have established, and will be updated as new data are released from the ABCD Study. This paper provides an overview of the initial results from the twin study embedded within the ABCD Study, an introduction to the primary research domains in the ABCD study and twin methodology, and an evaluation of the initial findings with a focus on data quality and suitability for future behavior genetic studies using the ABCD dataset. The broad introductory material is provided in recognition of the multidisciplinary appeal of the ABCD Study. While this paper focuses on univariate analyses, we emphasize the opportunities for multivariate, developmental and causal analyses, as well as those evaluating heterogeneity by key moderators such as sex, demographic factors and genetic background.

Evaluating the efficacy and mechanism of metformin targets on reducing Alzheimer's disease risk in the general population: a Mendelian randomisation study

AIMS/HYPOTHESIS

Metformin use has been associated with reduced incidence of dementia in diabetic individuals in observational studies. However, the causality between the two in the general population is unclear. This study uses Mendelian randomisation (MR) to investigate the causal effect of metformin targets on Alzheimer's disease and potential causal mechanisms in the brain linking the two.

METHODS

Genetic proxies for the effects of metformin drug targets were identified as variants in the gene for the corresponding target that associated with HbA1c level (N=344,182) and expression level of the corresponding gene (N≤31,684). The cognitive outcomes were derived from genome-wide association studies comprising 527,138 middle-aged Europeans, including 71,880 with Alzheimer's disease or Alzheimer's disease-by-proxy. MR estimates representing lifelong metformin use on Alzheimer's disease and cognitive function in the general population were generated. Effect of expression level of 22 metformin-related genes in brain cortex (N=6601 donors) on Alzheimer's disease was further estimated.

RESULTS

Genetically proxied metformin use, equivalent to a 6.75 mmol/mol (1.09%) reduction on HbA1c, was associated with 4% lower odds of Alzheimer's disease (OR 0.96 [95% CI 0.95, 0.98], p=1.06×10-4) in non-diabetic individuals. One metformin target, mitochondrial complex 1 (MCI), showed a robust effect on Alzheimer's disease (OR 0.88, p=4.73×10-4) that was independent of AMP-activated protein kinase. MR of expression in brain cortex tissue showed that decreased MCI-related gene (NDUFA2) expression was associated with lower Alzheimer's disease risk (OR 0.95, p=4.64×10-4) and favourable cognitive function.

CONCLUSIONS/INTERPRETATION

Metformin use may cause reduced Alzheimer's disease risk in the general population. Mitochondrial function and the NDUFA2 gene are plausible mechanisms of action in dementia protection.

Implicating Causal Brain Magnetic Resonance Imaging in Glaucoma Using Mendelian Randomization

Background

Glaucoma is hypothesized to originate in the brain but manifests as an eye disease as it possesses the common features of neurodegeneration diseases. But there is no evidence to demonstrate the primary brain changes in glaucoma patients. In the present study, we have used Mendelian randomization (MR) to understand the causal effect of brain alterations on glaucoma.

Methods

Our MR study was carried out using summary statistics from genome-wide associations for 110 diffusion tensor imaging (DTI) measurements of white matter (WM) tracts (17,706 individuals), 101 brain region-of-interest (ROI) volumes (19,629 individuals), and glaucoma (8,591 cases, 210,201 control subjects). The causal relationship was evaluated by multiplicative random effects inverse variance weighted (IVW) method and verified by two other MR methods, including MR Egger, weighted median, and extensive sensitivity analyses.

Results

Genetic liability to fornix fractional anisotropy (FX.FA) (OR = 0.71, 95%CI = 0.56–0.88, P = 2.44 × 10^–3), and uncinate fasciculus UNC.FA (OR = 0.65, 95%CI = 0.48–0.88, P = 5.57 × 10^–3) was associated with a low risk of glaucoma. Besides, the right ventral diencephalon (OR = 1.72, 95%CI = 1.17–2.52, P = 5.64 × 10^–3) and brain stem (OR = 1.35, 95%CI = 1.08–1.69, P = 8.94 × 10^–3) were associated with the increased risk of glaucoma. No heterogeneity and pleiotropy were detected.

Conclusion

Our study suggests that the fornix and uncinate fasciculus degenerations and injures of the right ventral diencephalon and brain stem potentially increase the occurrence of glaucoma and reveal the existence of the brain-eye axis.

Phenotypic and genetic associations of quantitative magnetic susceptibility in UK Biobank brain imaging

A key aim in epidemiological neuroscience is identification of markers to assess brain health and monitor therapeutic interventions. Quantitative susceptibility mapping (QSM) is an emerging magnetic resonance imaging technique that measures tissue magnetic susceptibility and has been shown to detect pathological changes in tissue iron, myelin and calcification. We present an open resource of QSM-based imaging measures of multiple brain structures in 35,273 individuals from the UK Biobank prospective epidemiological study. We identify statistically significant associations of 251 phenotypes with magnetic susceptibility that include body iron, disease, diet and alcohol consumption. Genome-wide associations relate magnetic susceptibility to 76 replicating clusters of genetic variants with biological functions involving iron, calcium, myelin and extracellular matrix. These patterns of associations include relationships that are unique to QSM, in particular being complementary to T2* signal decay time measures. These new imaging phenotypes are being integrated into the core UK Biobank measures provided to researchers worldwide, creating the potential to discover new, non-invasive markers of brain health.

Imaging Genetics in Epilepsy: Current Knowledge and New Perspectives

Epilepsy is a neurological network disease with genetics playing a much greater role than was previously appreciated. Unfortunately, the relationship between genetic basis and imaging phenotype is by no means simple. Imaging genetics integrates multidimensional datasets within a unified framework, providing a unique opportunity to pursue a global vision for epilepsy. This review delineates the current knowledge of underlying genetic mechanisms for brain networks in different epilepsy syndromes, particularly from a neural developmental perspective. Further, endophenotypes and their potential value are discussed. Finally, we highlight current challenges and provide perspectives for the future development of imaging genetics in epilepsy.

Diverse functions associate with non-coding polymorphisms shared between humans and chimpanzees

Background

Long-term balancing selection (LTBS) can maintain allelic variation at a locus over millions of years and through speciation events. Variants shared between species in the state of identity-by-descent, hereafter “trans-species polymorphisms”, can result from LTBS, often due to host–pathogen interactions. For instance, the major histocompatibility complex (MHC) locus contains TSPs present across primates. Several hundred candidate LTBS regions have been identified in humans and chimpanzees; however, because many are in non-protein-coding regions of the genome, the functions and potential adaptive roles for most remain unknown.

Results

We integrated diverse genomic annotations to explore the functions of 60 previously identified regions with multiple shared polymorphisms (SPs) between humans and chimpanzees, including 19 with strong evidence of LTBS. We analyzed genome-wide functional assays, expression quantitative trait loci (eQTL), genome-wide association studies (GWAS), and phenome-wide association studies (PheWAS) for all the regions. We identify functional annotations for 59 regions, including 58 with evidence of gene regulatory function from GTEx or functional genomics data and 19 with evidence of trait association from GWAS or PheWAS. As expected, the SPs associate in humans with many immune system phenotypes, including response to pathogens, but we also find associations with a range of other phenotypes, including body size, alcohol intake, cognitive performance, risk-taking behavior, and urate levels.

Conclusions

The diversity of traits associated with non-coding regions with multiple SPs support previous hypotheses that functions beyond the immune system are likely subject to LTBS. Furthermore, several of these trait associations provide support and candidate genetic loci for previous hypothesis about behavioral diversity in human and chimpanzee populations, such as the importance of variation in risk sensitivity.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-02020-x.

Multivariate genome-wide association study on tissue-sensitive diffusion metrics highlights pathways that shape the human brain

The molecular determinants of tissue composition of the human brain remain largely unknown. Recent genome-wide association studies (GWAS) on this topic have had limited success due to methodological constraints. Here, we apply advanced whole-brain analyses on multi-shell diffusion imaging data and multivariate GWAS to two large scale imaging genetic datasets (UK Biobank and the Adolescent Brain Cognitive Development study) to identify and validate genetic association signals. We discover 503 unique genetic loci that have impact on multiple regions of human brain. Among them, more than 79% are validated in either of two large-scale independent imaging datasets. Key molecular pathways involved in axonal growth, astrocyte-mediated neuroinflammation, and synaptogenesis during development are found to significantly impact the measured variations in tissue-specific imaging features. Our results shed new light on the biological determinants of brain tissue composition and their potential overlap with the genetic basis of neuropsychiatric disorders.

Genome-Wide Association Study of Alzheimer's Disease Brain Imaging Biomarkers and Neuropsychological Phenotypes in the European Medical Information Framework for Alzheimer's Disease Multimodal Biomarker Discovery Dataset

Alzheimer's disease (AD) is the most frequent neurodegenerative disease with an increasing prevalence in industrialized, aging populations. AD susceptibility has an established genetic basis which has been the focus of a large number of genome-wide association studies (GWAS) published over the last decade. Most of these GWAS used dichotomized clinical diagnostic status, i.e., case vs. control classification, as outcome phenotypes, without the use of biomarkers. An alternative and potentially more powerful study design is afforded by using quantitative AD-related phenotypes as GWAS outcome traits, an analysis paradigm that we followed in this work. Specifically, we utilized genotype and phenotype data from n = 931 individuals collected under the auspices of the European Medical Information Framework for Alzheimer's Disease Multimodal Biomarker Discovery (EMIF-AD MBD) study to perform a total of 19 separate GWAS analyses. As outcomes we used five magnetic resonance imaging (MRI) traits and seven cognitive performance traits. For the latter, longitudinal data from at least two timepoints were available in addition to cross-sectional assessments at baseline. Our GWAS analyses revealed several genome-wide significant associations for the neuropsychological performance measures, in particular those assayed longitudinally. Among the most noteworthy signals were associations in or near EHBP1 (EH domain binding protein 1; on chromosome 2p15) and CEP112 (centrosomal protein 112; 17q24.1) with delayed recall as well as SMOC2 (SPARC related modular calcium binding 2; 6p27) with immediate recall in a memory performance test. On the X chromosome, which is often excluded in other GWAS, we identified a genome-wide significant signal near IL1RAPL1 (interleukin 1 receptor accessory protein like 1; Xp21.3). While polygenic score (PGS) analyses showed the expected strong associations with SNPs highlighted in relevant previous GWAS on hippocampal volume and cognitive function, they did not show noteworthy associations with recent AD risk GWAS findings. In summary, our study highlights the power of using quantitative endophenotypes as outcome traits in AD-related GWAS analyses and nominates several new loci not previously implicated in cognitive decline.

Migraine, Stroke, and Cervical Arterial Dissection: Shared Genetics for a Triad of Brain Disorders With Vascular Involvement

BACKGROUND AND OBJECTIVES

Migraine, stroke, and cervical artery dissection (CeAD) represent a triad of cerebrovascular disorders with pairwise comorbid relationships and vascular involvement. Larger samples and recent advances in methodology invite systematic exploration of their shared genetics.

METHODS

Genetic analyses leveraged summary statistics from genome-wide association studies of the largest available samples of each disorder, including subtypes of stroke (ischemic stroke, large artery stroke, small vessel stroke, and cardioembolic stroke) and migraine (with aura and without aura). For each pair of disorders, genetic correlation was assessed both on a genome-wide basis and within independent segments across the genome including known specific loci for each disorder. A cross-trait meta-analysis was used to identify novel candidate loci. Finally, potential causality of migraine susceptibility on stroke and CeAD was assessed by Mendelian randomization.

RESULTS

Among all pairs of disorders, genome-wide genetic correlation was observed only between CeAD and migraine, particularly MO. Local genetic correlations were more extensive between migraine and CeAD than those between migraine and stroke or CeAD and stroke and revealed evidence for novel CeAD associations at rs6693567 (ADAMTSL4/ECM1), rs11187838 (PLCE1), and rs7940646 (MRVI1) while strengthening prior subthreshold evidence at rs9486725 (FHL5) and rs650724 (LRP1). At known migraine loci, novel associations with stroke had concordant risk alleles for small vessel stroke at rs191602009 (CARF) and for cardioembolic stroke at rs55884259 (NKX2-5). Known migraine loci also revealed novel associations but with opposite risk alleles for all stroke, ischemic stroke, and small vessel stroke at rs55928386 (HTRA1), for large artery stroke at rs11172113 (LRP1), and for all stroke and ischemic stroke at rs1535791 and rs4942561 (both LRCH1), respectively. rs182923402 (near PTCH1) was a novel concordant locus for migraine and cardioembolic stroke. Mendelian randomization supported potential causal influences of migraine on CeAD (odds ratio [95% confidence interval] per doubling migraine prevalence = 1.69 [1.24-2.3], p = 0.0009) with concordant risk, but with opposite risk on large artery stroke (0.86 [0.76-0.96], p = 0.0067).

DISCUSSION

The findings emphasize shared genetic risk between migraine and CeAD while identifying loci with likely vascular function in migraine and shared but opposite genetic risk between migraine and stroke subtypes, and a central role of LRP1 in all 3 cerebrovascular disorders.

A new Mendelian Randomization method to estimate causal effects of multivariable brain imaging exposures

The advent of simultaneously collected imaging-genetics data in large study cohorts provides an unprecedented opportunity to assess the causal effect of brain imaging traits on externally measured experimental results (e.g., cognitive tests) by treating genetic variants as instrumental variables. However, classic Mendelian Randomization methods are limited when handling high-throughput imaging traits as exposures to identify causal effects. We propose a new Mendelian Randomization framework to jointly select instrumental variables and imaging exposures, and then estimate the causal effect of multivariable imaging data on the outcome. We validate the proposed method with extensive data analyses and compare it with existing methods. We further apply our method to evaluate the causal effect of white matter microstructure integrity (WM) on cognitive function. The findings suggest that our method achieved better performance regarding sensitivity, bias, and false discovery rate compared to individually assessing the causal effect of a single exposure and jointly assessing the causal effect of multiple exposures without dimension reduction. Our application results indicated that WM measures across different tracts have a joint causal effect that significantly impacts the cognitive function among the participants from the UK Biobank.

Genetic variation, brain, and intelligence differences

Individual differences in human intelligence, as assessed using cognitive test scores, have a well-replicated, hierarchical phenotypic covariance structure. They are substantially stable across the life course, and are predictive of educational, social, and health outcomes. From this solid phenotypic foundation and importance for life, comes an interest in the environmental, social, and genetic aetiologies of intelligence, and in the foundations of intelligence differences in brain structure and functioning. Here, we summarise and critique the last 10 years or so of molecular genetic (DNA-based) research on intelligence, including the discovery of genetic loci associated with intelligence, DNA-based heritability, and intelligence's genetic correlations with other traits. We summarise new brain imaging-intelligence findings, including whole-brain associations and grey and white matter associations. We summarise regional brain imaging associations with intelligence and interpret these with respect to theoretical accounts. We address research that combines genetics and brain imaging in studying intelligence differences. There are new, though modest, associations in all these areas, and mechanistic accounts are lacking. We attempt to identify growing points that might contribute toward a more integrated 'systems biology' account of some of the between-individual differences in intelligence.

Genetic basis to structural grey matter associations with chronic pain

Structural neuroimaging studies of individuals with chronic pain conditions have often observed decreased regional grey matter at a phenotypic level. However, it is not known if this association can be attributed to genetic factors. Here we employed a novel integrative data-driven and hypothesis-testing approach to determine whether there is a genetic basis to grey matter morphology differences in chronic pain. Using publicly available genome-wide association study summary statistics for regional chronic pain conditions (n = 196 963) and structural neuroimaging measures (n = 19 629-34 000), we applied bivariate linkage disequilibrium-score regression and latent causal variable analyses to determine the genetic correlations (rG) and genetic causal proportion (GCP) between these complex traits, respectively. Five a priori brain regions (i.e. prefrontal cortex, cingulate cortex, insula, thalamus and superior temporal gyrus) were selected based on systematic reviews of grey matter morphology studies in chronic pain. Across this evidence-based selection of five brain regions, 10 significant negative genetic correlations (out of 369) were found (false discovery rate < 5%), suggesting a shared genetic basis to both reduced regional grey matter morphology and the presence of chronic pain. Specifically, negative genetic correlations were observed between reduced insula grey matter morphology and chronic pain in the abdomen (mean insula cortical thickness), hips (left insula volume) and neck/shoulders (left and right insula volume). Similarly, a shared genetic basis was found for reduced posterior cingulate cortex volume in chronic pain of the hip (left and right posterior cingulate), neck/shoulder (left posterior cingulate) and chronic pain at any site (left posterior cingulate); and for reduced pars triangularis volume in chronic neck/shoulder (left pars triangularis) and widespread pain (right pars triangularis). Across these negative genetic correlations, a significant genetic causal proportion was only found between mean insula thickness and chronic abdominal pain [rG (standard error, SE) = -0.25 (0.08), P = 1.06 × 10-3; GCP (SE) = -0.69 (0.20), P = 4.96 × 10-4]. This finding suggests that the genes underlying reduced cortical thickness of the insula causally contribute to an increased risk of chronic abdominal pain. Altogether, these results provide independent corroborating evidence for observational reports of decreased grey matter of particular brain regions in chronic pain. Further, we show for the first time that this association is mediated (in part) by genetic factors. These novel findings warrant further investigation into the neurogenetic pathways that underlie the development and prolongation of chronic pain conditions.