Genome-wide association studies of brain imaging phenotypes in UK Biobank

Impact of gene-by-trauma interaction in MDD-related multimorbidity clusters

BACKGROUND

Major depressive disorder (MDD) is considerably heterogeneous in terms of comorbidities, which may hamper the disentanglement of its biological mechanism. In a previous study, we classified the lifetime trajectories of MDD-related multimorbidities into seven distinct clusters, each characterized by unique genetic and environmental risk-factor profiles. The current objective was to investigate genome-wide gene-by-environment (G × E) interactions with childhood trauma burden, within the context of these clusters.

METHODS

We analyzed 77,519 participants and 6,266,189 single-nucleotide polymorphisms (SNPs) of the UK Biobank database. Childhood trauma burden was assessed using the Childhood Trauma Screener (CTS). For each cluster, Plink 2.0 was used to calculate SNP × CTS interaction effects on the participants' cluster membership probabilities. We especially focused on the effects of 31 candidate genes and associated SNPs selected from previous G × E studies for childhood maltreatment's association with depression.

RESULTS

At SNP-level, only the high-multimorbidity Cluster 6 revealed a genome-wide significant SNP rs145772219. At gene-level, MPST and PRH2 were genome-wide significant for the low-multimorbidity Clusters 1 and 3, respectively. Regarding candidate SNPs for G × E interactions, individual SNP results could be replicated for specific clusters. The candidate genes CREB1, DBH, and MTHFR (Cluster 5) as well as TPH1 (Cluster 6) survived multiple testing correction.

LIMITATIONS

CTS is a short retrospective self-reported measurement. Clusters could be influenced by genetics of individual disorders.

CONCLUSIONS

The first G × E GWAS for MDD-related multimorbidity trajectories successfully replicated findings from previous G × E studies related to depression, and revealed risk clusters for the contribution of childhood trauma.

Associations Between Temporomandibular Disorders and Brain Imaging-Derived Phenotypes

OBJECTIVE

Temporomandibular disorders (TMD) affect the temporomandibular joint and associated structures. Despite its prevalence and impact on quality of life, the underlying mechanisms of TMD remain unclear. Magnetic resonance imaging studies suggest brain abnormalities in patients with TMD. However, these lines of evidence are essentially observational and cannot infer a causal relationship. This study employs Mendelian randomisation (MR) to probe causal relationships between TMD and brain changes.

METHODS

Genome-wide association study (GWAS) summary statistics for TMD were collected, along with brain imaging-derived phenotypes (IDPs). Instrumental variables were selected from the GWAS summary statistics and used in bidirectional 2-sample MR analyses. The inverse-variance weighted analysis was chosen as the primary method. In addition, false discovery rate (FDR) correction of P value was used.

RESULTS

Eleven IDPs related to brain imaging alterations showed significant causal associations with TMD (P-FDR < .05), validated through sensitivity analysis. In forward MR, the mean thickness of left caudal middle frontal gyrus (OR, 0.76; 95% CI, 0.67-0.87; P-FDR = 1.15  ×  10-2) and the volume of right superior frontal gyrus (OR, 1.24; 95% CI, 1.10-1.39; P-FDR = 2.26  ×  10-2) exerted significant causal effects on TMD. In the reverse MR analysis, TMD exerted a significant causal effect on 9 IDPs, including the mean thickness of the left medial orbitofrontal cortex (β = -0.10; 95% CI, -0.13 to -0.08; P-FDR = 2.06 × 10-11), the volume of the left magnocellular nucleus (β = -0.15; 95% CI, -0.22 to -0.09; P-FDR = 3.26 × 10-4), the mean intensity of the right inferior-lateral ventricle (β = -0.09; 95% CI, -0.14 to -0.04; P-FDR = 2.23 × 10-2), the volume of grey matter in the anterior division of the left superior temporal gyrus (β = 0.09; 95% CI,  0.04-0.14; P-FDR = 1.69 × 10-2), and so forth.

CONCLUSIONS

This study provides genetic evidence supporting the bidirectional causal associations between TMD and brain IDPs, shedding light on potential neurobiological mechanisms underlying TMD development and its relationship with brain structure.

Mendelian randomization analyses support causal relationships between tinnitus of different stages and severity and structural characteristics of specific brain regions

This study aims to delineate the causal relationships between idiopathic tinnitus in different stages and severity and the morphological properties in specific brain regions. We utilized a two-sample bidirectional Mendelian randomization (MR) analysis to ascertain the causal effects of brain structural attributes on varying severities and stages of tinnitus. Our approach involved harnessing genetic variables derived from extensive genome-wide association studies as instrumental variables, centered mainly on pertinent single-nucleotide polymorphisms associated with tinnitus. Subsequently, we integrated this data with brain structural imaging inputs to facilitate the MR analysis. We also applied reverse MR analysis to pinpoint the critical brain regions implicated in the onset of tinnitus. Our analysis revealed a demonstrable causal relationship between tinnitus and brain structural alterations, including changes primarily within the auditory cortex and hub regions of the limbic system, as well as portions of the frontal-temporal-occipital circuit. We found that individuals exhibiting cortical thickness alterations in the bilateral peri-calcarine and right superior occipital gyrus might have previously experienced tinnitus. Changes in the cortical areas of the right rectus, left inferior frontal gyrus, and right pars-orbitalis appeared unrelated to tinnitus. Furthermore, moderate tinnitus patients showed more pronounced structural alterations. This study substantiates that tinnitus could instigate substantial structural alterations mainly within the auditory-limbic-frontal-visual system, while the reciprocal causality was not supported. Moreover, the data underscores that moderate, rather than severe, tinnitus precipitates the most significant structural changes. Morphological alterations in several specific brain areas either indicate a history of tinnitus or bear no relation to it.

Unsupervised representation learning on high-dimensional clinical data improves genomic discovery and prediction

Although high-dimensional clinical data (HDCD) are increasingly available in biobank-scale datasets, their use for genetic discovery remains challenging. Here we introduce an unsupervised deep learning model, Representation Learning for Genetic Discovery on Low-Dimensional Embeddings (REGLE), for discovering associations between genetic variants and HDCD. REGLE leverages variational autoencoders to compute nonlinear disentangled embeddings of HDCD, which become the inputs to genome-wide association studies (GWAS). REGLE can uncover features not captured by existing expert-defined features and enables the creation of accurate disease-specific polygenic risk scores (PRSs) in datasets with very few labeled data. We apply REGLE to perform GWAS on respiratory and circulatory HDCD-spirograms measuring lung function and photoplethysmograms measuring blood volume changes. REGLE replicates known loci while identifying others not previously detected. REGLE are predictive of overall survival, and PRSs constructed from REGLE loci improve disease prediction across multiple biobanks. Overall, REGLE contain clinically relevant information beyond that captured by existing expert-defined features, leading to improved genetic discovery and disease prediction.

Shared Genetic Architecture Among Gastrointestinal Diseases, Schizophrenia, and Brain Subcortical Volumes

BACKGROUND AND HYPOTHESIS

The gut-brain axis plays important roles in both gastrointestinal diseases (GI diseases) and schizophrenia (SCZ). Moreover, both GI diseases and SCZ exhibit notable abnormalities in brain subcortical volumes. However, the genetic mechanisms underlying the comorbidity of these diseases and the shared alterations in brain subcortical volumes remain unclear.

STUDY DESIGN

Using the genome-wide association studies data of SCZ, 14 brain subcortical volumes, and 8 GI diseases, the global polygenic overlap and local genetic correlations were identified, as well as the shared genetic variants among those phenotypes. Furthermore, we conducted multi-trait colocalization analyses to bolster our findings. Functional annotations, cell-type enrichment, and protein-protein interaction (PPI) analyses were carried out to reveal the critical etiology and pathology mechanisms.

STUDY RESULTS

The global polygenic overlap and local genetic correlations informed the close relationships between SCZ and both GI diseases and brain subcortical volumes. Moreover, 84 unique lead-shared variants were identified. The associated genes were linked to vital biological processes within the immune system. Additionally, significant correlations were observed with key immune cells and the PPI analysis identified several histone-associated hub genes. These findings highlighted the pivotal roles played by the immune system for both SCZ and GI diseases, along with the shared alterations in brain subcortical volumes.

CONCLUSIONS

These findings revealed the shared genetic architecture contributing to SCZ and GI diseases, as well as their shared alterations in brain subcortical volumes. These insights have substantial implications for the concurrent development of intervention and therapy targets for these diseases.

Whole-exome sequencing identifies protein-coding variants associated with brain iron in 29,828 individuals

Iron plays a fundamental role in multiple brain disorders. However, the genetic underpinnings of brain iron and its implications for these disorders are still lacking. Here, we conduct an exome-wide association analysis of brain iron, measured by quantitative susceptibility mapping technique, across 26 brain regions among 26,789 UK Biobank participants. We find 36 genes linked to brain iron, with 29 not being previously reported, and 16 of them can be replicated in an independent dataset with 3,039 subjects. Many of these genes are involved in iron transport and homeostasis, such as FTH1 and MLX. Several genes, while not previously connected to brain iron, are associated with iron-related brain disorders like Parkinson's (STAB1, KCNA10), Alzheimer's (SHANK1), and depression (GFAP). Mendelian randomization analysis reveals six causal relationships from regional brain iron to brain disorders, such as from the hippocampus to depression and from the substantia nigra to Parkinson's. These insights advance our understanding of the genetic architecture of brain iron and offer potential therapeutic targets for brain disorders.

Estimated bone mineral density and white matter hyperintensities: A bidirectional Mendelian randomization study

PURPOSE

Greater white matter hyperintensities (WMH) in older adults have been associated with reduced bone mineral density (BMD) and increased fractures and falls. However, it is unclear whether there is a causal relationship between BMD reduction and WMH. In this study, Mendelian randomization (MR) was used to find the causality between WMH and estimated BMD (eBMD).

METHODS

We performed a two-sample bidirectional MR analysis using statistical data obtained from publicly available genome-wide association studies (GWAS). The main method of MR analysis is the inverse-variance weighted (IVW) method. To identify and account for the impact of horizontal pleiotropy, we also employed MR-Egger regression, MR pleiotropy residual sum, and outlier (MR-PRESSO).

RESULTS

MR analysis found a causal relationship between eBMD and WMH (IVW OR = 0.938, 95 % CI: 0.889-0.990, p = 0.020). Our causal estimates are unlikely to be distorted by horizontal pleiotropy according to heterogeneity test (both p > 0.05) and MR-Egger regression (p > 0.05). However, in the reverse MR analysis, there was no evidence that WMH was causally correlated with eBMD (IVW OR = 0.979, 95 % CI: 0.954-1.005, p = 0.109).

CONCLUSION

Our results suggest that low eBMD increased the risk of WMH; conversely, no evidence that WMH causally affects eBMD was found.

A transcriptomic atlas of the human brain reveals genetically determined aspects of neuropsychiatric health

Regulation of gene expression is a vital component of neurological homeostasis. Cataloging the consequences of endogenous gene expression on the physical structure and connectivity of the brain offers a means of unifying trait-associated genetic variation with trait-associated neurological features. We perform tissue-specific transcriptome-wide association studies (TWASs) on over 3,400 neuroimaging phenotypes in the UK Biobank (N = 33,224) using our joint-tissue imputation (JTI)-TWAS method. We identify highly significant associations between predicted expression for 7,192 genes and a wide variety of measures of the brain derived from magnetic resonance imaging (MRI). Our approach generates reproducible results in internal and external replication datasets. Genetically determined expression alone is sufficient for high-fidelity reconstruction of brain structure and organization. We demonstrate complementary benefits of cross-tissue and single-tissue analyses toward an integrated neurobiology and provide evidence that gene expression outside the central nervous system provides unique insights into brain health. As an application, we provide evidence suggesting that the genetically regulated expression of schizophrenia risk genes causally affects over 73% of neurological phenotypes that are altered in individuals with schizophrenia (as identified by neuroimaging studies). Imaging features associated with neuropsychiatric traits can provide valuable insights into underlying pathophysiology. By linking neuroimaging-derived phenotypes with expression levels of specific genes, this resource represents a powerful gene prioritization schema that can improve our understanding of brain function, development, and disease. The use of multiple different cortical and subcortical atlases in the resource facilitates direct integration of these data with findings from a diverse range of clinical neuroimaging studies.

Deciphering the functional specialization of whole-brain spatiomolecular gradients in the adult brain

Cortical arealization arises during neurodevelopment from the confluence of molecular gradients representing patterned expression of morphogens and transcription factors. However, whether similar gradients are maintained in the adult brain remains unknown. Here, we uncover three axes of topographic variation in gene expression in the adult human brain that specifically capture previously identified rostral-caudal, dorsal-ventral, and medial-lateral axes of early developmental patterning. The interaction of these spatiomolecular gradients i) accurately reconstructs the position of brain tissue samples, ii) delineates known functional territories, and iii) can model the topographical variation of diverse cortical features. The spatiomolecular gradients are distinct from canonical cortical axes differentiating the primary sensory cortex from the association cortex, but radiate in parallel with the axes traversed by local field potentials along the cortex. We replicate all three molecular gradients in three independent human datasets as well as two nonhuman primate datasets and find that each gradient shows a distinct developmental trajectory across the lifespan. The gradients are composed of several well-known transcription factors (e.g., PAX6 and SIX3), and a small set of genes shared across gradients are strongly enriched for multiple diseases. Together, these results provide insight into the developmental sculpting of functionally distinct brain regions, governed by three robust transcriptomic axes embedded within brain parenchyma.

Genetic influences on brain and cognitive health and their interactions with cardiovascular conditions and depression

Approximately 40% of dementia cases could be prevented or delayed by modifiable risk factors related to lifestyle and environment. These risk factors, such as depression and vascular disease, do not affect all individuals in the same way, likely due to inter-individual differences in genetics. However, the precise nature of how genetic risk profiles interact with modifiable risk factors to affect brain health is poorly understood. Here we combine multiple data resources, including genotyping and postmortem gene expression, to map the genetic landscape of brain structure and identify 367 loci associated with cortical thickness and 13 loci associated with white matter hyperintensities (P < 5×10-8), with several loci also showing a significant association with cognitive function. We show that among 220 unique genetic loci associated with cortical thickness in our genome-wide association studies (GWAS), 95 also showed evidence of interaction with depression or cardiovascular conditions. Polygenic risk scores based on our GWAS of inferior frontal thickness also interacted with hypertension in predicting executive function in the Canadian Longitudinal Study on Aging. These findings advance our understanding of the genetic underpinning of brain structure and show that genetic risk for brain and cognitive health is in part moderated by treatable mid-life factors.

Casual associations between brain structure and sarcopenia: A large-scale genetic correlation and mendelian randomization study

Sarcopenia presenting a critical challenge in population-aging healthcare. The elucidation of the interplay between brain structure and sarcopenia necessitates further research. The aim of this study is to explore the casual association between brain structure and sarcopenia. Linkage disequilibrium score regression (LDSC) was conducted to estimate the genetic correlations; MR was then performed to explore the causal relationship between Brain imaging-derived phenotypes (BIDPs) and three sarcopenia-related traits: handgrip strength, walking pace, and appendicular lean mass (ALM). The main analyses were conducted using the inverse-variance weighted method. Moreover, weighted median and MR-Egger were conducted as sensitivity analyses. Genetic association between 6.41% of BIDPs and ALM was observed, and 4.68% of BIDPs exhibited causal MR association with handgrip strength, 2.11% of BIDPs were causally associated with walking pace, and 2.04% of BIDPs showed causal association with ALM. Volume of ventromedial hypothalamus was associated with increased odds of handgrip strength (OR: 1.18, 95% CI: 1.02 to 1.37) and ALM (OR: 1.05, 95% CI: 1.01 to 1.09). Mean thickness of G-pariet-inf-Angular was associated with decreased odds of handgrip strength (OR: 0.83, 95% CI: 0.70 to 0.97) and walking pace (OR: 0.97, 95% CI: 0.93 to 0.99). As part of the brain structure forward causally influences sarcopenia, which may provide new perspectives for the prevention of sarcopenia and offer valuable insights for further research on the brain-muscle axis.

White Matter Tract Integrity Is Reduced in Depression and in Individuals With Genetic Liability to Depression

BACKGROUND

While major depression has been linked to changes in white matter architecture, it remains unclear whether risk factors for depression are directly associated with these alterations. We reexamined white matter fiber tracts in individuals with depressive symptoms and investigated the connection between genetic and environmental risk for depression and structural changes in the brain.

METHODS

We included 19,183 participants from the UK Biobank imaging cohort, with depression status and adverse life experience based on questionnaire data and genetic liability for depression quantified by polygenic scores. The integrity of 27 white matter tracts was assessed using mean fractional anisotropy derived from diffusion magnetic resonance imaging.

RESULTS

White matter integrity was reduced, particularly in thalamic and intracortical fiber tracts, in individuals with depressive symptoms, independent of current symptom status. In a group of healthy individuals without depression, increasing genetic risk and increasing environmental risk were associated with reduced integrity in relevant fiber tracts, particularly in thalamic radiations. This association was stronger than expected based on statistical dependencies between samples, as confirmed by subsequent in silico simulations and permutation tests.

CONCLUSIONS

White matter alterations in thalamic and association tracts are associated with depressive symptoms and genetic risk for depression in unaffected individuals, suggesting an intermediate phenotype at the brain level.

Quantitative susceptibility mapping in the brain reflects spatial expression of genes involved in iron homeostasis and myelination

Quantitative susceptibility mapping (QSM) is an MRI modality used to non-invasively measure iron content in the brain. Iron exhibits a specific anatomically varying pattern of accumulation in the brain across individuals. The highest regions of accumulation are the deep grey nuclei, where iron is stored in paramagnetic molecule ferritin. This form of iron is considered to be what largely contributes to the signal measured by QSM in the deep grey nuclei. It is also known that QSM is affected by diamagnetic myelin contents. Here, we investigate spatial gene expression of iron and myelin related genes, as measured by the Allen Human Brain Atlas, in relation to QSM images of age-matched subjects. We performed multiple linear regressions between gene expression and the average QSM signal within 34 distinct deep grey nuclei regions. Our results show a positive correlation (p < .05, corrected) between expression of ferritin and the QSM signal in deep grey nuclei regions. We repeated the analysis for other genes that encode proteins thought to be involved in the transport and storage of iron in the brain, as well as myelination. In addition to ferritin, our findings demonstrate a positive correlation (p < .05, corrected) between the expression of ferroportin, transferrin, divalent metal transporter 1, several gene markers of myelinating oligodendrocytes, and the QSM signal in deep grey nuclei regions. Our results suggest that the QSM signal reflects both the storage and active transport of iron in the deep grey nuclei regions of the brain.

Synthetic surrogates improve power for genome-wide association studies of partially missing phenotypes in population biobanks

Within population biobanks, incomplete measurement of certain traits limits the power for genetic discovery. Machine learning is increasingly used to impute the missing values from the available data. However, performing genome-wide association studies (GWAS) on imputed traits can introduce spurious associations, identifying genetic variants that are not associated with the original trait. Here we introduce a new method, synthetic surrogate (SynSurr) analysis, which makes GWAS on imputed phenotypes robust to imputation errors. Rather than replacing missing values, SynSurr jointly analyzes the original and imputed traits. We show that SynSurr estimates the same genetic effect as standard GWAS and improves power in proportion to the quality of the imputations. SynSurr requires a commonly made missing-at-random assumption but relaxes the requirements of existing imputation methods by not requiring correct model specification. We present extensive simulations and ablation analyses to validate SynSurr and apply it to empower the GWAS of dual-energy X-ray absorptiometry traits within the UK Biobank.

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.

Cross-ancestry genome-wide association studies of brain imaging phenotypes

Genome-wide association studies of brain imaging phenotypes are mainly performed in European populations, but other populations are severely under-represented. Here, we conducted Chinese-alone and cross-ancestry genome-wide association studies of 3,414 brain imaging phenotypes in 7,058 Chinese Han and 33,224 white British participants. We identified 38 new associations in Chinese-alone analyses and 486 additional new associations in cross-ancestry meta-analyses at P < 1.46 × 10-11 for discovery and P < 0.05 for replication. We pooled significant autosomal associations identified by single- or cross-ancestry analyses into 6,443 independent associations, which showed uneven distribution in the genome and the phenotype subgroups. We further divided them into 44 associations with different effect sizes and 3,557 associations with similar effect sizes between ancestries. Loci of these associations were shared with 15 brain-related non-imaging traits including cognition and neuropsychiatric disorders. Our results provide a valuable catalog of genetic associations for brain imaging phenotypes in more diverse populations.

Multiscale Modes of Functional Brain Connectivity

Information processing in the brain spans from localised sensorimotor processes to higher-level cognition that integrates across multiple regions. Interactions between and within these subsystems enable multiscale information processing. Despite this multiscale characteristic, functional brain connectivity is often either estimated based on 10-30 distributed modes or parcellations with 100-1000 localised parcels, both missing across-scale functional interactions. We present Multiscale Probabilistic Functional Modes (mPFMs), a new mapping which comprises modes over various scales of granularity, thus enabling direct estimation of functional connectivity within- and across-scales. Crucially, mPFMs emerged from data-driven multilevel Bayesian modelling of large functional MRI (fMRI) populations. We demonstrate that mPFMs capture both distributed brain modes and their co-existing subcomponents. In addition to validating mPFMs using simulations and real data, we show that mPFMs can predict ~900 personalised traits from UK Biobank more accurately than current standard techniques. Therefore, mPFMs can offer a paradigm shift in functional connectivity modelling and yield enhanced fMRI biomarkers for traits and diseases.

Genome-wide meta-analyses of restless legs syndrome yield insights into genetic architecture, disease biology and risk prediction

Restless legs syndrome (RLS) affects up to 10% of older adults. Their healthcare is impeded by delayed diagnosis and insufficient treatment. To advance disease prediction and find new entry points for therapy, we performed meta-analyses of genome-wide association studies in 116,647 individuals with RLS (cases) and 1,546,466 controls of European ancestry. The pooled analysis increased the number of risk loci eightfold to 164, including three on chromosome X. Sex-specific meta-analyses revealed largely overlapping genetic predispositions of the sexes (rg = 0.96). Locus annotation prioritized druggable genes such as glutamate receptors 1 and 4, and Mendelian randomization indicated RLS as a causal risk factor for diabetes. Machine learning approaches combining genetic and nongenetic information performed best in risk prediction (area under the curve (AUC) = 0.82-0.91). In summary, we identified targets for drug development and repurposing, prioritized potential causal relationships between RLS and relevant comorbidities and risk factors for follow-up and provided evidence that nonlinear interactions are likely relevant to RLS risk prediction.

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.

Neuroanatomical, transcriptomic, and molecular correlates of math ability and their prognostic value for predicting learning outcomes

Foundational mathematical abilities, acquired in early childhood, are essential for success in our technology-driven society. Yet, the neurobiological mechanisms underlying individual differences in children's mathematical abilities and learning outcomes remain largely unexplored. Leveraging one of the largest multicohort datasets from children at a pivotal stage of knowledge acquisition, we first establish a replicable mathematical ability-related imaging phenotype (MAIP). We then show that brain gene expression profiles enriched for candidate math ability-related genes, neuronal signaling, synaptic transmission, and voltage-gated potassium channel activity contributed to the MAIP. Furthermore, the similarity between MAIP gene expression signatures and brain structure, acquired before intervention, predicted learning outcomes in two independent math tutoring cohorts. These findings advance our knowledge of the interplay between neuroanatomical, transcriptomic, and molecular mechanisms underlying mathematical ability and reveal predictive biomarkers of learning. Our findings have implications for the development of personalized education and interventions.

The power of many brains: Catalyzing neuropsychiatric discovery through open neuroimaging data and large-scale collaboration

Recent advances in open neuroimaging data are enhancing our comprehension of neuropsychiatric disorders. By pooling images from various cohorts, statistical power has increased, enabling the detection of subtle abnormalities and robust associations, and fostering new research methods. Global collaborations in imaging have furthered our knowledge of the neurobiological foundations of brain disorders and aided in imaging-based prediction for more targeted treatment. Large-scale magnetic resonance imaging initiatives are driving innovation in analytics and supporting generalizable psychiatric studies. We also emphasize the significant role of big data in understanding neural mechanisms and in the early identification and precise treatment of neuropsychiatric disorders. However, challenges such as data harmonization across different sites, privacy protection, and effective data sharing must be addressed. With proper governance and open science practices, we conclude with a projection of how large-scale imaging resources and collaborations could revolutionize diagnosis, treatment selection, and outcome prediction, contributing to optimal brain health.

Bidirectional two-sample Mendelian randomization analyses support causal relationships between structural and diffusion imaging-derived phenotypes and the risk of major neurodegenerative diseases

Previous observational investigations suggest that structural and diffusion imaging-derived phenotypes (IDPs) are associated with major neurodegenerative diseases; however, whether these associations are causal remains largely uncertain. Herein we conducted bidirectional two-sample Mendelian randomization analyses to infer the causal relationships between structural and diffusion IDPs and major neurodegenerative diseases using common genetic variants-single nucleotide polymorphism (SNPs) as instrumental variables. Summary statistics of genome-wide association study (GWAS) for structural and diffusion IDPs were obtained from 33,224 individuals in the UK Biobank cohort. Summary statistics of GWAS for seven major neurodegenerative diseases were obtained from the largest GWAS for each disease to date. The forward MR analyses identified significant or suggestively statistical causal effects of genetically predicted three structural IDPs on Alzheimer's disease (AD), frontotemporal dementia (FTD), and multiple sclerosis. For example, the reduction in the surface area of the left superior temporal gyrus was associated with a higher risk of AD. The reverse MR analyses identified significantly or suggestively statistical causal effects of genetically predicted AD, Lewy body dementia (LBD), and FTD on nine structural and diffusion IDPs. For example, LBD was associated with increased mean diffusivity in the right superior longitudinal fasciculus and AD was associated with decreased gray matter volume in the right ventral striatum. Our findings might contribute to shedding light on the prediction and therapeutic intervention for the major neurodegenerative diseases at the neuroimaging level.

Shared genetics linking sociability with the brain's default mode network

The brain's default mode network (DMN) plays a role in social cognition, with altered DMN function being associated with social impairments across various neuropsychiatric disorders. In the present study, we examined the genetic relationship between sociability and DMN-related resting-state functional magnetic resonance imaging (rs-fMRI) traits. To this end, we used genome-wide association summary statistics for sociability and 31 activity and 64 connectivity DMN-related rs-fMRI traits (N=34,691-342,461). First, we examined global and local genetic correlations between sociability and the rs-fMRI traits. Second, to assess putatively causal relationships between the traits, we conducted bi-directional Mendelian randomisation (MR) analyses. Finally, we prioritised genes influencing both sociability and rs-fMRI traits by combining three methods: gene-expression eQTL MR analyses, the CELLECT framework using single-nucleus RNA-seq data, and network propagation in the context of a protein-protein interaction network. Significant local genetic correlations were found between sociability and two rs-fMRI traits, one representing spontaneous activity within the temporal cortex, the other representing connectivity between the frontal/cingulate and angular/temporal cortices. Sociability affected 12 rs-fMRI traits when allowing for weakly correlated genetic instruments. Combing all three methods for gene prioritisation, we defined 17 highly prioritised genes, with DRD2 and LINGO1 showing the most robust evidence across all analyses. By integrating genetic and transcriptomics data, our gene prioritisation strategy may serve as a blueprint for future studies. The prioritised genes could be explored as potential biomarkers for social dysfunction in the context of neuropsychiatric disorders and as drug target genes.

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.

Association of Gene Expression and Tremor Network Structure

BACKGROUND

Transcriptomic changes in the essential tremor (ET)-associated cerebello-thalamo-cortical "tremor network" and their association to brain structure have not been investigated.

OBJECTIVE

The aim was to characterize molecular changes associated with network-level imaging-derived phenotypes (IDP) found in ET.

METHODS

We performed an imaging-transcriptomic study in British adults using imaging-genome-wide association study summary statistics (UK Biobank "BIG40" cohort; n = 33,224, aged 40-69 years). We imputed imaging-transcriptomic associations for 184 IDPs and analyzed functional enrichment of gene modules and aggregate network-level phenotypes. Validation was performed in cerebellar-tissue RNA-sequencing data from ET patients and controls (n = 55).

RESULTS

Among 237,896 individual predicted gene expression levels for 6063 unique genes/transcripts, we detected 2269 genome-wide significant associations (Bonferroni P < 2.102e-7, 0.95%). These were concentrated in intracellular volume fraction measures of white matter pathways and in genes with putative links to tremor (MAPT, ARL17A, KANSL1, SPPL2C, LRRC37A4P, PLEKHM1, and FMNL1). Whole-tremor-network cortical thickness was associated with a gene module linked to mitochondrial organization and protein quality control (r = 0.91, P = 2e-70), whereas white-gray T1-weighted magnetic resonance imaging (MRI) contrast in the tremor network was associated with a gene module linked to sphingolipid synthesis and ethanolamine metabolism (r = -0.90, P = 2e-68). Imputed association effect sizes and RNA-sequencing log-fold change in the validation dataset were significantly correlated for cerebellar peduncular diffusion MRI phenotypes, and there was a close overlap of significant associations between both datasets for gray matter phenotypes (χ2 = 6.40, P = 0.006).

CONCLUSIONS

The identified genes and processes are potential treatment targets for ET, and our results help characterize molecular changes that could in future be used for patient treatment selection or prognosis prediction. © 2024 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

The pivotal role of the X-chromosome in the genetic architecture of the human brain

Genes on the X-chromosome are extensively expressed in the human brain. However, little is known for the X-chromosome's impact on the brain anatomy, microstructure, and functional network. We examined 1,045 complex brain imaging traits from 38,529 participants in the UK Biobank. We unveiled potential autosome-X-chromosome interactions, while proposing an atlas outlining dosage compensation (DC) for brain imaging traits. Through extensive association studies, we identified 72 genome-wide significant trait-locus pairs (including 29 new associations) that share genetic architectures with brain-related disorders, notably schizophrenia. Furthermore, we discovered unique sex-specific associations and assessed variations in genetic effects between sexes. Our research offers critical insights into the X-chromosome's role in the human brain, underscoring its contribution to the differences observed in brain structure and functionality between sexes.

Mendelian randomization analyses reveal causal relationships between brain functional networks and risk of psychiatric disorders

Dysfunction of brain resting-state functional networks has been widely reported in psychiatric disorders. However, the causal relationships between brain resting-state functional networks and psychiatric disorders remain largely unclear. Here we perform bidirectional two-sample Mendelian randomization (MR) analyses to investigate the causalities between 191 resting-state functional magnetic resonance imaging (rsfMRI) phenotypes (n = 34,691 individuals) and 12 psychiatric disorders (n = 14,307 to 698,672 individuals). Forward MR identified 8 rsfMRI phenotypes causally associated with the risk of psychiatric disorders. For example, the increase in the connectivity of motor, subcortical-cerebellum and limbic network was associated with lower risk of autism spectrum disorder. In adddition, increased connectivity in the default mode and central executive network was associated with lower risk of post-traumatic stress disorder and depression. Reverse MR analysis revealed significant associations between 4 psychiatric disorders and 6 rsfMRI phenotypes. For instance, the risk of attention-deficit/hyperactivity disorder increases the connectivity of the attention, salience, motor and subcortical-cerebellum network. The risk of schizophrenia mainly increases the connectivity of the default mode and central executive network and decreases the connectivity of the attention network. In summary, our findings reveal causal relationships between brain functional networks and psychiatric disorders, providing important interventional and therapeutic targets for psychiatric disorders at the brain functional network level.

SGLT1 and SGLT2 inhibition, circulating metabolites, and cerebral small vessel disease: a mediation Mendelian Randomization study

BACKGROUND

Sodium-glucose cotransporter 2 (SGLT2) and SGLT1 inhibitors may have additional beneficial metabolic effects on circulating metabolites beyond glucose regulation, which could contribute to a reduction in the burden of cerebral small vessel disease (CSVD). Accordingly, we used Mendelian Randomization (MR) to examine the role of circulating metabolites in mediating SGLT2 and SGLT1 inhibition in CSVD.

METHODS

Genetic instruments for SGLT1/2 inhibition were identified as genetic variants, which were both associated with the expression of encoding genes of SGLT1/2 inhibitors and glycated hemoglobin A1c (HbA1c) level. A two-sample two-step MR was used to determine the causal effects of SGLT1/2 inhibition on CSVD manifestations and the mediating effects of 1400 circulating metabolites linking SGLT1/2 inhibition with CSVD manifestations.

RESULTS

A lower risk of deep cerebral microbleeds (CMBs) and small vessel stroke (SVS) was linked to genetically predicted SGLT2 inhibition. Better white matter structure integrity was also achieved, as evidenced by decreased mean diffusivity (MD), axial diffusivity (AD), and radial diffusivity (RD), as well as lower deep (DWMH) and periventrivular white matter hyperintensity (PWMH) volume. Inhibiting SGLT2 could also lessen the incidence of severe enlarged perivascular spaces (EPVS) located at white matter, basal ganglia (BG) and hippocampus (HIP). SGLT1 inhibition could preserve white matter integrity, shown as decreased MD of white matter and DWMH volume. The effect of SGLT2 inhibition on SVS and MD of white matter through the concentration of 4-acetamidobutanoate and the cholesterol to oleoyl-linoleoyl-glycerol (18:1 to 18:2) ratio, with a mediated proportion of 30.3% and 35.5% of the total effect, respectively.

CONCLUSIONS

SGLT2 and SGLT1 inhibition play protective roles in CSVD development. The SGLT2 inhibition could lower the risk of SVS and improve the integrity of white matter microstructure via modulating the level of 4-acetamidobutanoate and cholesterol metabolism. Further mechanistic and clinical studies research are needed to validate our findings.

Unveiling Promising Neuroimaging Biomarkers for Schizophrenia Through Clinical and Genetic Perspectives

Schizophrenia is a complex and serious brain disorder. Neuroscientists have become increasingly interested in using magnetic resonance-based brain imaging-derived phenotypes (IDPs) to investigate the etiology of psychiatric disorders. IDPs capture valuable clinical advantages and hold biological significance in identifying brain abnormalities. In this review, we aim to discuss current and prospective approaches to identify potential biomarkers for schizophrenia using clinical multimodal neuroimaging and imaging genetics. We first described IDPs through their phenotypic classification and neuroimaging genomics. Secondly, we discussed the applications of multimodal neuroimaging by clinical evidence in observational studies and randomized controlled trials. Thirdly, considering the genetic evidence of IDPs, we discussed how can utilize neuroimaging data as an intermediate phenotype to make association inferences by polygenic risk scores and Mendelian randomization. Finally, we discussed machine learning as an optimum approach for validating biomarkers. Together, future research efforts focused on neuroimaging biomarkers aim to enhance our understanding of schizophrenia.

Respiratory modulation of the heart rate: A potential biomarker of cardiorespiratory function in human

In recent decade, wearable digital devices have shown potentials for the discovery of novel biomarkers of humans' physiology and behavior. Heart rate (HR) and respiration rate (RR) are most crucial bio-signals in humans' digital phenotyping research. HR is a continuous and non-invasive proxy to autonomic nervous system and ample evidence pinpoints the critical role of respiratory modulation of cardiac function. In the present study, we recorded longitudinal (7 days, 4.63 ± 1.52) HR and RR of 89 freely behaving human subjects (Female: 39, age 57.28 ± 5.67, Male: 50, age 58.48 ± 6.32) and analyzed their dynamics using linear models and information theoretic measures. While HR's linear and nonlinear characteristics were expressed within the plane of the HR-RR directed flow of information (HR→RR - RR→HR), their dynamics were determined by its RR→HR axis. More importantly, RR→HR quantified the effect of alcohol consumption on individuals' cardiorespiratory function independent of their consumed amount of alcohol, thereby signifying the presence of this habit in their daily life activities. The present findings provided evidence for the critical role of the respiratory modulation of HR, which was previously only studied in non-human animals. These results can contribute to humans' phenotyping research by presenting RR→HR as a digital diagnosis/prognosis marker of humans' cardiorespiratory pathology.

Computational limits to the legibility of the imaged human brain

Our knowledge of the organisation of the human brain at the population-level is yet to translate into power to predict functional differences at the individual-level, limiting clinical applications and casting doubt on the generalisability of inferred mechanisms. It remains unknown whether the difficulty arises from the absence of individuating biological patterns within the brain, or from limited power to access them with the models and compute at our disposal. Here we comprehensively investigate the resolvability of such patterns with data and compute at unprecedented scale. Across 23 810 unique participants from UK Biobank, we systematically evaluate the predictability of 25 individual biological characteristics, from all available combinations of structural and functional neuroimaging data. Over 4526 GPU*hours of computation, we train, optimize, and evaluate out-of-sample 700 individual predictive models, including fully-connected feed-forward neural networks of demographic, psychological, serological, chronic disease, and functional connectivity characteristics, and both uni- and multi-modal 3D convolutional neural network models of macro- and micro-structural brain imaging. We find a marked discrepancy between the high predictability of sex (balanced accuracy 99.7%), age (mean absolute error 2.048 years, R2 0.859), and weight (mean absolute error 2.609Kg, R2 0.625), for which we set new state-of-the-art performance, and the surprisingly low predictability of other characteristics. Neither structural nor functional imaging predicted an individual's psychology better than the coincidence of common chronic disease (p < 0.05). Serology predicted chronic disease (p < 0.05) and was best predicted by it (p < 0.001), followed by structural neuroimaging (p < 0.05). Our findings suggest either more informative imaging or more powerful models will be needed to decipher individual level characteristics from the human brain. We make our models and code openly available.

Human cortical neurogenesis is altered via glucocorticoid-mediated regulation of ZBTB16 expression

Glucocorticoids are important for proper organ maturation, and their levels are tightly regulated during development. Here, we use human cerebral organoids and mice to study the cell-type-specific effects of glucocorticoids on neurogenesis. We show that glucocorticoids increase a specific type of basal progenitors (co-expressing PAX6 and EOMES) that has been shown to contribute to cortical expansion in gyrified species. This effect is mediated via the transcription factor ZBTB16 and leads to increased production of neurons. A phenome-wide Mendelian randomization analysis of an enhancer variant that moderates glucocorticoid-induced ZBTB16 levels reveals causal relationships with higher educational attainment and altered brain structure. The relationship with postnatal cognition is also supported by data from a prospective pregnancy cohort study. This work provides a cellular and molecular pathway for the effects of glucocorticoids on human neurogenesis that relates to lasting postnatal phenotypes.

Causal association and mediating effect of blood biochemical metabolic traits and brain image-derived endophenotypes on Alzheimer's disease

Background

Recent genetic evidence supports that circulating biochemical and metabolic traits (BMTs) play a causal role in Alzheimer's disease (AD), which might be mediated by changes in brain structure. Here, we leveraged publicly available genome-wide association study data to investigate the intrinsic causal relationship between blood BMTs, brain image-derived phenotypes (IDPs) and AD.

Methods

Utilizing the genetic variants associated with 760 blood BMTs and 172 brain IDPs as the exposure and the latest AD summary statistics as the outcome, we analyzed the causal relationship between blood BMTs and brain IDPs and AD by using a two-sample Mendelian randomization (MR) method. Additionally, we used two-step/mediation MR to study the mediating effect of brain IDPs between blood BMTs and AD.

Results

Twenty-five traits for genetic evidence supporting a causal association with AD were identified, including 12 blood BMTs and 13 brain IDPs. For BMTs, glutamine consistently reduced the risk of AD in 3 datasets. For IDPs, specific alterations of cortical thickness (atrophy in frontal pole and insular lobe, and incrassation in superior parietal lobe) and subcortical volume (atrophy in hippocampus and its subgroups, left accumbens and left choroid plexus, and expansion in cerebral white matter) are vulnerable to AD. In the two-step/mediation MR analysis, superior parietal lobe, right hippocampal fissure and left accumbens were identified to play a potential mediating role among three blood BMTs and AD.

Conclusions

The results obtained in our study suggest that 12 circulating BMTs and 13 brain IDPs play a causal role in AD. Importantly, a subset of BMTs exhibit shared genetic architecture and potentially causal relationships with brain structure, which may contribute to the alteration of brain IDPs in AD.

Iron and risk of dementia: Mendelian randomisation analysis in UK Biobank

BACKGROUND

Brain iron deposition is common in dementia, but whether serum iron is a causal risk factor is unknown. We aimed to determine whether genetic predisposition to higher serum iron status biomarkers increased risk of dementia and atrophy of grey matter.

METHODS

We analysed UK Biobank participants clustered into European (N=451284), African (N=7477) and South Asian (N=9570) groups by genetic similarity to the 1000 genomes project. Using Mendelian randomisation methods, we estimated the association between genetically predicted serum iron (transferrin saturation [TSAT] and ferritin), grey matter volume and genetic liability to clinically defined dementia (including Alzheimer's disease [AD], non-AD dementia, and vascular dementia) from hospital and primary care records. We also performed time-to-event (competing risks) analysis of the TSAT polygenic score on risk of clinically defined non-AD dementia.

RESULTS

In Europeans, higher genetically predicted TSAT increased genetic liability to dementia (Odds Ratio [OR]: 1.15, 95% Confidence Intervals [CI] 1.04 to 1.26, p=0.0051), non-AD dementia (OR: 1.27, 95% CI 1.12 to 1.45, p=0.00018) and vascular dementia (OR: 1.37, 95% CI 1.12 to 1.69, p=0.0023), but not AD (OR: 1.00, 95% CI 0.86 to 1.15, p=0.97). Higher TSAT was also associated with increased risk of non-AD dementia in participants of African, but not South Asian groups. In survival analysis using a TSAT polygenic score, the effect was independent of apolipoprotein-E ε4 genotype (with adjustment subdistribution Hazard Ratio: 1.74, 95% CI 1.33 to 2.28, p=0.00006). Genetically predicted TSAT was associated with lower grey matter volume in caudate, putamen and thalamus, and not in other areas of interest.

DISCUSSION

Genetic evidence supports a causal relationship between higher TSAT and risk of clinically defined non-AD and vascular dementia, in European and African groups. This association appears to be independent of apolipoprotein-E ε4.

Childhood maltreatment influences adult brain structure through its effects on immune, metabolic, and psychosocial factors

Childhood maltreatment (CM) leads to a lifelong susceptibility to mental ill-health which might be reflected by its effects on adult brain structure, perhaps indirectly mediated by its effects on adult metabolic, immune, and psychosocial systems. Indexing these systemic factors via body mass index (BMI), C-reactive protein (CRP), and rates of adult trauma (AT), respectively, we tested three hypotheses: (H1) CM has direct or indirect effects on adult trauma, BMI, and CRP; (H2) adult trauma, BMI, and CRP are all independently related to adult brain structure; and (H3) childhood maltreatment has indirect effects on adult brain structure mediated in parallel by BMI, CRP, and AT. Using path analysis and data from N = 116,887 participants in UK Biobank, we find that CM is related to greater BMI and AT levels, and that these two variables mediate CM's effects on CRP [H1]. Regression analyses on the UKB MRI subsample (N = 21,738) revealed that greater CRP and BMI were both independently related to a spatially convergent pattern of cortical effects (Spearman's ρ = 0.87) characterized by fronto-occipital increases and temporo-parietal reductions in thickness. Subcortically, BMI was associated with greater volume, AT with lower volume and CPR with effects in both directions [H2]. Finally, path models indicated that CM has indirect effects in a subset of brain regions mediated through its direct effects on BMI and AT and indirect effects on CRP [H3]. Results provide evidence that childhood maltreatment can influence brain structure decades after exposure by increasing individual risk toward adult trauma, obesity, and inflammation.

Genetic evidence suggests a genetic association between major depressive disorder and reduced cortical gray matter volume: A Mendelian randomization study and mediation analysis

BACKGROUND

Several studies have suggested an association between major depressive disorder (MDD) and abnormal brain structure. However, it is unclear whether MDD affects cortical gray matter volume, a common indicator of cognitive function. We aimed to determine whether MDD was associated with decreased cortical gray matter volume (GMV) through a Mendelian randomization (MR) study.

METHODS

We obtained summary genetic data from a study conducted by the Psychiatric Genomics Consortium, which recruited a total of 480,359 participants (135,458 cases and 344,901 controls). Genetic tools-single nucleotide polymorphisms (SNPs)-of MDD were extracted from the study and their effects on gray matter volumes of the cortex and total brain were evaluated in a large cohort from the UK Biobank (n = 8427). The effects of the SNPs were pooled using inverse variance weighted (IVW) analysis and further tested in several sensitivity analyses. We tested whether C-reactive protein (CRP) levels and interleukin-6 signaling were the mediators of the effects using a multivariate MR model.

RESULTS

Thirty-three SNPs were identified and adopted as genetic tools for predicting MDD. IVW analysis showed that MDD was associated with lower overall GMV (beta value -0.106, 95%CI -0.188 to -0.023, p = 0.011) in the frontal pole (left frontal pole, -0.152, 95%CI -0.177 to -0.127, p = 0.013; right frontal pole, -0.133, 95%CI -0.253 to -0.013, p = 0.028). Multivariate and mediation analysis showed that interleukin-6 was an important mediator of GMV reduction. Reverse causality analysis found no evidence that total GMV affected the risk of MDD, but showed that increased left precuneus cortex volume and left posterior cingulate cortex volume were associated with increased risk of MDD.

LIMITATIONS

Potential pleiotropic effects and overestimation of real-world effects. Key assumptions for MR analysis may not be satisfactorily met.

CONCLUSION

MDD was associated with a reduced GMV, and interleukin-6 might be a mediator of GMV reduction.

Reweighting UK Biobank corrects for pervasive selection bias due to volunteering

BACKGROUND

Biobanks typically rely on volunteer-based sampling. This results in large samples (power) at the cost of representativeness (bias). The problem of volunteer bias is debated. Here, we (i) show that volunteering biases associations in UK Biobank (UKB) and (ii) estimate inverse probability (IP) weights that correct for volunteer bias in UKB.

METHODS

Drawing on UK Census data, we constructed a subsample representative of UKB's target population, which consists of all individuals invited to participate. Based on demographic variables shared between the UK Census and UKB, we estimated IP weights (IPWs) for each UKB participant. We compared 21 weighted and unweighted bivariate associations between these demographic variables to assess volunteer bias.

RESULTS

Volunteer bias in all associations, as naively estimated in UKB, was substantial-in some cases so severe that unweighted estimates had the opposite sign of the association in the target population. For example, older individuals in UKB reported being in better health, in contrast to evidence from the UK Census. Using IPWs in weighted regressions reduced 87% of volunteer bias on average. Volunteer-based sampling reduced the effective sample size of UKB substantially, to 32% of its original size.

CONCLUSIONS

Estimates from large-scale biobanks may be misleading due to volunteer bias. We recommend IP weighting to correct for such bias. To aid in the construction of the next generation of biobanks, we provide suggestions on how to best ensure representativeness in a volunteer-based design. For UKB, IPWs have been made available.

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.

Childhood urbanicity is associated with emotional episodic memory-related striatal function and common variation in NTRK2

BACKGROUND

Childhoods in urban or rural environments may differentially affect the risk of neuropsychiatric disorders, possibly through memory processing and neural response to emotional stimuli. Genetic factors may not only influence individuals' choices of residence but also modulate how the living environment affects responses to episodic memory.

METHODS

We investigated the effects of childhood urbanicity on episodic memory in 410 adults (discovery sample) and 72 adults (replication sample) with comparable socioeconomic statuses in Beijing, China, distinguishing between those with rural backgrounds (resided in rural areas before age 12 and relocated to urban areas at or after age 12) and urban backgrounds (resided in cities before age 12). We examined the effect of childhood urbanicity on brain function across encoding and retrieval sessions using an fMRI episodic memory paradigm involving the processing of neutral or aversive pictures. Moreover, genetic association analyses were conducted to understand the potential genetic underpinnings that might contribute to memory processing and neural mechanisms influenced by early-life urban or rural environments.

RESULTS

Episodic memory retrieval accuracy for more difficult neutral stimuli was similar between those with urban and rural childhoods, whereas aversive stimuli elicited higher retrieval accuracy in the urban group (P = 0.023). For aversive stimuli, subjects with urban childhood had relatively decreased engagement of the striatum at encoding and decreased engagement of the hippocampus at retrieval. This more efficient striatal encoding of aversive stimuli in those with urban childhoods was associated with common variation in neurotrophic tyrosine kinase receptor type 2 (NTRK2) (right striatum: P = 1.58×10-6). These findings were confirmed in the replication sample.

CONCLUSIONS

We suggest that this differential striatal processing of aversive stimuli observed in individuals with urban or rural childhoods may represent mechanisms by which childhood urbanicity may affect brain circuits, heightening behavioral responses to negative stressors associated with urban environments. NTRK2-associated neural processes in the striatum may play a role in these processes.

FPLS-DC: functional partial least squares through distance covariance for imaging genetics

MOTIVATION

Imaging genetics integrates imaging and genetic techniques to examine how genetic variations influence the function and structure of organs like the brain or heart, providing insights into their impact on behavior and disease phenotypes. The use of organ-wide imaging endophenotypes has increasingly been used to identify potential genes associated with complex disorders. However, analyzing organ-wide imaging data alongside genetic data presents two significant challenges: high dimensionality and complex relationships. To address these challenges, we propose a novel, nonlinear inference framework designed to partially mitigate these issues.

RESULTS

We propose a functional partial least squares through distance covariance (FPLS-DC) framework for efficient genome wide analyses of imaging phenotypes. It consists of two components. The first component utilizes the FPLS-derived base functions to reduce image dimensionality while screening genetic markers. The second component maximizes the distance correlation between genetic markers and projected imaging data, which is a linear combination of the FPLS-basis functions, using simulated annealing algorithm. In addition, we proposed an iterative FPLS-DC method based on FPLS-DC framework, which effectively overcomes the influence of inter-gene correlation on inference analysis. We efficiently approximate the null distribution of test statistics using a gamma approximation. Compared to existing methods, FPLS-DC offers computational and statistical efficiency for handling large-scale imaging genetics. In real-world applications, our method successfully detected genetic variants associated with the hippocampus, demonstrating its value as a statistical toolbox for imaging genetic studies.

AVAILABILITY AND IMPLEMENTATION

The FPLS-DC method we propose opens up new research avenues and offers valuable insights for analyzing functional and high-dimensional data. In addition, it serves as a useful tool for scientific analysis in practical applications within the field of imaging genetics research. The R package FPLS-DC is available in Github: https://github.com/BIG-S2/FPLSDC.

The genetic architecture of multimodal human brain age

The complex biological mechanisms underlying human brain aging remain incompletely understood. This study investigated the genetic architecture of three brain age gaps (BAG) derived from gray matter volume (GM-BAG), white matter microstructure (WM-BAG), and functional connectivity (FC-BAG). We identified sixteen genomic loci that reached genome-wide significance (P-value < 5×10-8). A gene-drug-disease network highlighted genes linked to GM-BAG for treating neurodegenerative and neuropsychiatric disorders and WM-BAG genes for cancer therapy. GM-BAG displayed the most pronounced heritability enrichment in genetic variants within conserved regions. Oligodendrocytes and astrocytes, but not neurons, exhibited notable heritability enrichment in WM and FC-BAG, respectively. Mendelian randomization identified potential causal effects of several chronic diseases on brain aging, such as type 2 diabetes on GM-BAG and AD on WM-BAG. Our results provide insights into the genetics of human brain aging, with clinical implications for potential lifestyle and therapeutic interventions. All results are publicly available at https://labs.loni.usc.edu/medicine .

Utilizing multimodal AI to improve genetic analyses of cardiovascular traits

Electronic health records, biobanks, and wearable biosensors contain multiple high-dimensional clinical data (HDCD) modalities (e.g., ECG, Photoplethysmography (PPG), and MRI) for each individual. Access to multimodal HDCD provides a unique opportunity for genetic studies of complex traits because different modalities relevant to a single physiological system (e.g., circulatory system) encode complementary and overlapping information. We propose a novel multimodal deep learning method, M-REGLE, for discovering genetic associations from a joint representation of multiple complementary HDCD modalities. We showcase the effectiveness of this model by applying it to several cardiovascular modalities. M-REGLE jointly learns a lower representation (i.e., latent factors) of multimodal HDCD using a convolutional variational autoencoder, performs genome wide association studies (GWAS) on each latent factor, then combines the results to study the genetics of the underlying system. To validate the advantages of M-REGLE and multimodal learning, we apply it to common cardiovascular modalities (PPG and ECG), and compare its results to unimodal learning methods in which representations are learned from each data modality separately, but the downstream genetic analyses are performed on the combined unimodal representations. M-REGLE identifies 19.3% more loci on the 12-lead ECG dataset, 13.0% more loci on the ECG lead I + PPG dataset, and its genetic risk score significantly outperforms the unimodal risk score at predicting cardiac phenotypes, such as atrial fibrillation (Afib), in multiple biobanks.

Bayesian mixed model inference for genetic association under related samples with brain network phenotype

Genetic association studies for brain connectivity phenotypes have gained prominence due to advances in noninvasive imaging techniques and quantitative genetics. Brain connectivity traits, characterized by network configurations and unique biological structures, present distinct challenges compared to other quantitative phenotypes. Furthermore, the presence of sample relatedness in the most imaging genetics studies limits the feasibility of adopting existing network-response modeling. In this article, we fill this gap by proposing a Bayesian network-response mixed-effect model that considers a network-variate phenotype and incorporates population structures including pedigrees and unknown sample relatedness. To accommodate the inherent topological architecture associated with the genetic contributions to the phenotype, we model the effect components via a set of effect network configurations and impose an inter-network sparsity and intra-network shrinkage to dissect the phenotypic network configurations affected by the risk genetic variant. A Markov chain Monte Carlo (MCMC) algorithm is further developed to facilitate uncertainty quantification. We evaluate the performance of our model through extensive simulations. By further applying the method to study, the genetic bases for brain structural connectivity using data from the Human Connectome Project with excessive family structures, we obtain plausible and interpretable results. Beyond brain connectivity genetic studies, our proposed model also provides a general linear mixed-effect regression framework for network-variate outcomes.

Analyzing bivariate cross-trait genetic architecture in GWAS summary statistics with the BIGA cloud computing platform

As large-scale biobanks provide increasing access to deep phenotyping and genomic data, genome-wide association studies (GWAS) are rapidly uncovering the genetic architecture behind various complex traits and diseases. GWAS publications typically make their summary-level data (GWAS summary statistics) publicly available, enabling further exploration of genetic overlaps between phenotypes gathered from different studies and cohorts. However, systematically analyzing high-dimensional GWAS summary statistics for thousands of phenotypes can be both logistically challenging and computationally demanding. In this paper, we introduce BIGA (https://bigagwas.org/), a website that aims to offer unified data analysis pipelines and processed data resources for cross-trait genetic architecture analyses using GWAS summary statistics. We have developed a framework to implement statistical genetics tools on a cloud computing platform, combined with extensive curated GWAS data resources. Through BIGA, users can upload data, submit jobs, and share results, providing the research community with a convenient tool for consolidating GWAS data and generating new insights.

Exploring functional dysconnectivity in schizophrenia: alterations in eigenvector centrality mapping and insights into related genes from transcriptional profiles

Schizophrenia is a mental health disorder characterized by functional dysconnectivity. Eigenvector centrality mapping (ECM) has been employed to investigate alterations in functional connectivity in schizophrenia, yet the results lack consistency, and the genetic mechanisms underlying these changes remain unclear. In this study, whole-brain voxel-wise ECM analyses were conducted on resting-state functional magnetic resonance imaging data. A cohort of 91 patients with schizophrenia and 91 matched healthy controls were included during the discovery stage. Additionally, in the replication stage, 153 individuals with schizophrenia and 182 healthy individuals participated. Subsequently, a comprehensive analysis was performed using an independent transcriptional database derived from six postmortem healthy adult brains to explore potential genetic factors influencing the observed functional dysconnectivity, and to investigate the roles of identified genes in neural processes and pathways. The results revealed significant and reliable alterations in the ECM across multiple brain regions in schizophrenia. Specifically, there was a significant decrease in ECM in the bilateral superior and middle temporal gyrus, and an increase in the bilateral thalamus in both the discovery and replication stages. Furthermore, transcriptional analysis revealed 420 genes whose expression patterns were related to changes in ECM, and these genes were enriched mainly in biological processes associated with synaptic signaling and transmission. Together, this study enhances our knowledge of the neural processes and pathways involved in schizophrenia, shedding light on the genetic factors that may be linked to functional dysconnectivity in this disorder.

A comprehensive analysis of APOE genotype effects on human brain structure in the UK Biobank

Alzheimer's disease (AD) risk is increased in carriers of the apolipoprotein E (APOE) ε4 allele and decreased in ε2 allele carriers compared with the ε3ε3 genotype. The aim of this study was to determine whether: the APOE genotype affects brain grey (GM) or white matter (WM) structure; and if differences exist, the age when they become apparent and whether there are differential effects by sex. We used cross-sectional magnetic resonance imaging data from ~43,000 (28,494 after pre-processing) white British cognitively healthy participants (7,446 APOE ε4 carriers) aged 45-80 years from the UK Biobank cohort and investigated image-derived phenotypes (IDPs). We observed no statistically significant effects of APOE genotype on GM structure volumes or median T2* in subcortical structures, a measure related to iron content. The volume of white matter hyperintensities differed significantly between APOE genotype groups with higher volumes in APOE ε4ε4 (effect size 0.14 standard deviations [SD]) and ε3ε4 carriers (effect size 0.04 SD) but no differences in ε2 carriers compared with ε3ε3 carriers. WM integrity measures in the dorsal (mean diffusivity [MD]) and ventral cingulum (MD and intracellular volume fraction), posterior thalamic radiation (MD and isotropic volume fraction) and sagittal stratum (MD) indicated lower integrity in APOE ε4ε4 carriers (effect sizes around 0.2-0.3 SD) and ε3ε4 (effect sizes around 0.05 SD) carriers but no differences in ε2 carriers compared with the APOE ε3ε3 genotype. Effects did not differ between men and women. APOE ε4 homozygotes had lower WM integrity specifically at older ages with a steeper decline of WM integrity from the age of 60 that corresponds to around 5 years greater "brain age". APOE genotype affects various white matters measures, which might be indicative of preclinical AD processes. This hypothesis can be assessed in future when clinical outcomes become available.

Integration of whole-exome sequencing and structural neuroimaging analysis in major depressive disorder: a joint study

Major depressive disorder (MDD) is a common mental illness worldwide and is triggered by an intricate interplay between environmental and genetic factors. Although there are several studies on common variants in MDD, studies on rare variants are relatively limited. In addition, few studies have examined the genetic contributions to neurostructural alterations in MDD using whole-exome sequencing (WES). We performed WES in 367 patients with MDD and 161 healthy controls (HCs) to detect germline and copy number variations in the Korean population. Gene-based rare variants were analyzed to investigate the association between the genes and individuals, followed by neuroimaging-genetic analysis to explore the neural mechanisms underlying the genetic impact in 234 patients with MDD and 135 HCs using diffusion tensor imaging data. We identified 40 MDD-related genes and observed 95 recurrent regions of copy number variations. We also discovered a novel gene, FRMPD3, carrying rare variants that influence MDD. In addition, the single nucleotide polymorphism rs771995197 in the MUC6 gene was significantly associated with the integrity of widespread white matter tracts. Moreover, we identified 918 rare exonic missense variants in genes associated with MDD susceptibility. We postulate that rare variants of FRMPD3 may contribute significantly to MDD, with a mild penetration effect.

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.

Causal Associations of Genetically Determined Tinnitus With Neuroimaging Traits: Evidence From a Mendelian Randomization Study

OBJECTIVES

Potential reverse causality and unmeasured confounding factors are common biases in most neuroimaging studies on tinnitus and central correlates. The causal association of tinnitus with neuroimaging features also remains unclear. This study aimed to investigate the causal relationship of tinnitus with neuroplastic alterations using Mendelian randomization.

DESIGN

Summary-level data from a genome-wide association study of tinnitus were derived from UK Biobank (n = 117,882). The genome-wide association study summary statistics for 4 global-brain tissue and 14 sub-brain gray matter volumetric traits were also obtained (n = up to 33,224). A bidirectional Mendelian randomization analysis was conducted to explore the causal relationship between tinnitus and neuroanatomical features at global-brain and sub-brain levels.

RESULTS

Genetic susceptibility to tinnitus was causally associated with increased white matter volume (odds ratio [OR] = 2.361, 95% confidence interval [CI], 1.033 to 5.393) and total brain volume (OR = 2.391, 95% CI, 1.047 to 5.463) but inversely associated with cerebrospinal fluid volume (OR = 0.362, 95% CI, 0.158 to 0.826). A smaller gray matter volume in the left Heschl's gyrus and right insular cortex and larger gray matter volume in the posterior division of the left parahippocampal gyrus may lead to an increased risk for tinnitus (OR = 0.978, 95% CI, 0.961 to 0.996; OR = 0.987, 95% CI, 0.976 to 0.998; and OR = 1.015, 95% CI, 1.001 to 1.028, respectively).

CONCLUSIONS

Genetic susceptibility to tinnitus was causally associated with increased white matter volume and total brain volume. Volume alteration in several cortical regions may indicate a higher tinnitus risk, and further research is recommended for causality inference at the level of sub-brain regions. Our findings provide genetic evidence for elucidating the underlying pathophysiological mechanisms of tinnitus-related neuroanatomical abnormalities.