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

Recent advances in psychoradiology

Psychiatry, as a field, lacks objective markers for diagnosis, progression, treatment planning, and prognosis, in part due to difficulties studying the brain in vivo, and diagnoses are based on self-reported symptoms and observation of patient behavior and cognition. Rapid advances in brain imaging techniques allow clinical investigators to noninvasively quantify brain features at the structural, functional, and molecular levels. Psychoradiology is an emerging discipline at the intersection of psychiatry and radiology. Psychoradiology applies medical imaging technologies to psychiatry and promises not only to improve insight into structural and functional brain abnormalities in patients with psychiatric disorders but also to have potential clinical utility. We searched for representative studies related to recent advances in psychoradiology through May 1, 2022, and conducted a selective review of 165 references, including 75 research articles. We summarize the novel dynamic imaging processing methods to model brain networks and present imaging genetics studies that reveal the relationship between various neuroimaging endophenotypes and genetic markers in psychiatric disorders. Furthermore, we survey recent advances in psychoradiology, with a focus on future psychiatric diagnostic approaches with dimensional analysis and a shift from group-level to individualized analysis. Finally, we examine the application of machine learning in psychoradiology studies and the potential of a novel option for brain stimulation treatment based on psychoradiological findings in precision medicine. Here, we provide a summary of recent advances in psychoradiology research, and we hope this review will help guide the practice of psychoradiology in the scientific and clinical fields.

Mapping brain endophenotypes associated with idiopathic pulmonary fibrosis genetic risk

BACKGROUND

Idiopathic pulmonary fibrosis (IPF) is a serious disease of the lung parenchyma. It has a known polygenetic risk, with at least seventeen regions of the genome implicated to date. Growing evidence suggests linked multimorbidity of IPF with neurodegenerative or affective disorders. However, no study so far has explicitly explored links between IPF, associated genetic risk profiles, and specific brain features.

METHODS

We exploited imaging and genetic data from more than 32,000 participants available through the UK Biobank population-level resource to explore links between IPF genetic risk and imaging-derived brain endophenotypes. We performed a brain-wide imaging-genetics association study between the presence of 17 known IPF risk variants and 1248 multi-modal imaging-derived features, which characterise brain structure and function.

FINDINGS

We identified strong associations between cortical morphological features, white matter microstructure and IPF risk loci in chromosomes 17 (17q21.31) and 8 (DEPTOR). Through co-localisation analysis, we confirmed that cortical thickness in the anterior cingulate and more widespread white matter microstructure changes share a single causal variant with IPF at the chromosome 8 locus. Post-hoc preliminary analysis suggested that forced vital capacity may partially mediate the association between the DEPTOR variant and white matter microstructure, but not between the DEPTOR risk variant and cortical thickness.

INTERPRETATION

Our results reveal the associations between IPF genetic risk and differences in brain structure, for both cortex and white matter. Differences in tissue-specific imaging signatures suggest distinct underlying mechanisms with focal cortical thinning in regions with known high DEPTOR expression, unrelated to lung function, and more widespread microstructural white matter changes consistent with hypoxia or neuroinflammation with potential mediation by lung function.

FUNDING

This study was supported by the NIHR Nottingham Biomedical Research Centre and the UK Medical Research Council.

Variant brain-derived neurotrophic factor val66met polymorphism engages memory-associated systems to augment olfaction

The neurogenetic basis of variability in human olfactory function remains elusive. This study examined olfactory performance and resting-state functional neuroimaging results from healthy volunteers within the context of the brain-derived neurotrophic factor (BDNF) val66met polymorphism with the aim of unraveling the genotype-associated intrinsic reorganization of the olfactory network. We found that the presence of the Met allele is associated with better olfactory identification and additional engagement of semantic memory system within the olfactory network, in an allele dosage-dependent manner. This suggests that the Met allele may promote adaptive neural reorganization to augment olfactory capacity.

Advancing brain network models to reconcile functional neuroimaging and clinical research

Functional magnetic resonance imaging (fMRI) captures information on brain function beyond the anatomical alterations that are traditionally visually examined by neuroradiologists. However, the fMRI signals are complex in addition to being noisy, so fMRI still faces limitations for clinical applications. Here we review methods that have been proposed as potential solutions so far, namely statistical, biophysical and decoding models, with their strengths and weaknesses. We especially evaluate the ability of these models to directly predict clinical variables from their parameters (predictability) and to extract clinically relevant information regarding biological mechanisms and relevant features for classification and prediction (interpretability). We then provide guidelines for useful applications and pitfalls of such fMRI-based models in a clinical research context, looking beyond the current state of the art. In particular, we argue that the clinical relevance of fMRI calls for a new generation of models for fMRI data, which combine the strengths of both biophysical and decoding models. This leads to reliable and biologically meaningful model parameters, which thus fulfills the need for simultaneous interpretability and predictability. In our view, this synergy is fundamental for the discovery of new pharmacological and interventional targets, as well as the use of models as biomarkers in neurology and psychiatry.

A causal association of ANKRD37 with human hippocampal volume

Human hippocampal volume has been separately associated with single nucleotide polymorphisms (SNPs), DNA methylation and gene expression, but their causal relationships remain largely unknown. Here, we aimed at identifying the causal relationships of SNPs, DNA methylation, and gene expression that are associated with hippocampal volume by integrating cross-omics analyses with genome editing, overexpression and causality inference. Based on structural neuroimaging data and blood-derived genome, transcriptome and methylome data, we prioritized a possibly causal association across multiple molecular phenotypes: rs1053218 mutation leads to cg26741686 hypermethylation, thus leads to overactivation of the associated ANKRD37 gene expression in blood, a gene involving hypoxia, which may result in the reduction of human hippocampal volume. The possibly causal relationships from rs1053218 to cg26741686 methylation to ANKRD37 expression obtained from peripheral blood were replicated in human hippocampal tissue. To confirm causality, we performed CRISPR-based genome and epigenome-editing of rs1053218 homologous alleles and cg26741686 methylation in mouse neural stem cell differentiation models, and overexpressed ANKRD37 in mouse hippocampus. These in-vitro and in-vivo experiments confirmed that rs1053218 mutation caused cg26741686 hypermethylation and ANKRD37 overexpression, and cg26741686 hypermethylation favored ANKRD37 overexpression, and ANKRD37 overexpression reduced hippocampal volume. The pairwise relationships of rs1053218 with hippocampal volume, rs1053218 with cg26741686 methylation, cg26741686 methylation with ANKRD37 expression, and ANKRD37 expression with hippocampal volume could be replicated in an independent healthy young (n = 443) dataset and observed in elderly people (n = 194), and were more significant in patients with late-onset Alzheimer's disease (n = 76). This study revealed a novel causal molecular association mechanism of ANKRD37 with human hippocampal volume, which may facilitate the design of prevention and treatment strategies for hippocampal impairment.

Generalization of cortical MOSTest genome-wide associations within and across samples

Genome-Wide Association studies have typically been limited to univariate analysis in which a single outcome measure is tested against millions of variants. Recent work demonstrates that a Multivariate Omnibus Statistic Test (MOSTest) is well powered to discover genomic effects distributed across multiple phenotypes. Applied to cortical brain MRI morphology measures, MOSTest has resulted in a drastic improvement in power to discover loci when compared to established approaches (min-P). One question that arises is how well these discovered loci replicate in independent data. Here we perform 10 times cross validation within 34,973 individuals from UK Biobank for imaging measures of cortical area, thickness and sulcal depth (>1,000 dimensionality for each). By deploying a replication method that aggregates discovered effects distributed across multiple phenotypes, termed PolyVertex Score (MOSTest-PVS), we demonstrate a higher replication yield and comparable replication rate of discovered loci for MOSTest (# replicated loci: 242-496, replication rate: 96-97%) in independent data when compared with the established min-P approach (# replicated loci: 26-55, replication rate: 91-93%). An out-of-sample replication of discovered loci was conducted with a sample of 4,069 individuals from the Adolescent Brain Cognitive Development® (ABCD) study, who are on average 50 years younger than UK Biobank individuals. We observe a higher replication yield and comparable replication rate of MOSTest-PVS compared to min-P. This finding underscores the importance of using well-powered multivariate techniques for both discovery and replication of high dimensional phenotypes in Genome-Wide Association studies.

Genetic map of regional sulcal morphology in the human brain from UK biobank data

Genetic associations with macroscopic brain structure can provide insights into brain function and disease. However, specific associations with measures of local brain folding are largely under-explored. Here, we conducted large-scale genome- and exome-wide associations of regional cortical sulcal measures derived from magnetic resonance imaging scans of 40,169 individuals in UK Biobank. We discovered 388 regional brain folding associations across 77 genetic loci, with genes in associated loci enriched for expression in the cerebral cortex, neuronal development processes, and differential regulation during early brain development. We integrated brain eQTLs to refine genes for various loci, implicated several genes involved in neurodevelopmental disorders, and highlighted global genetic correlations with neuropsychiatric phenotypes. We provide an interactive 3D visualisation of our summary associations, emphasising added resolution of regional analyses. Our results offer new insights into the genetic architecture of brain folding and provide a resource for future studies of sulcal morphology in health and disease.

Mapping the genetic architecture of cortical morphology through neuroimaging: progress and perspectives

Cortical morphology is a key determinant of cognitive ability and mental health. Its development is a highly intricate process spanning decades, involving the coordinated, localized expression of thousands of genes. We are now beginning to unravel the genetic architecture of cortical morphology, thanks to the recent availability of large-scale neuroimaging and genomic data and the development of powerful biostatistical tools. Here, we review the progress made in this field, providing an overview of the lessons learned from genetic studies of cortical volume, thickness, surface area, and folding as captured by neuroimaging. It is now clear that morphology is shaped by thousands of genetic variants, with effects that are region- and time-dependent, thereby challenging conventional study approaches. The most recent genome-wide association studies have started discovering common genetic variants influencing cortical thickness and surface area, yet together these explain only a fraction of the high heritability of these measures. Further, the impact of rare variants and non-additive effects remains elusive. There are indications that the quickly increasing availability of data from whole-genome sequencing and large, deeply phenotyped population cohorts across the lifespan will enable us to uncover much of the missing heritability in the upcoming years. Novel approaches leveraging shared information across measures will accelerate this process by providing substantial increases in statistical power, together with more accurate mapping of genetic relationships. Important challenges remain, including better representation of understudied demographic groups, integration of other 'omics data, and mapping of effects from gene to brain to behavior across the lifespan.

Using neuroimaging genomics to investigate the evolution of human brain structure

Alterations in brain size and organization represent some of the most distinctive changes in the emergence of our species. Yet, there is limited understanding of how genetic factors contributed to altered neuroanatomy during human evolution. Here, we analyze neuroimaging and genetic data from up to 30,000 people in the UK Biobank and integrate with genomic annotations for different aspects of human evolution, including those based on ancient DNA and comparative genomics. We show that previously reported signals of recent polygenic selection for cortical anatomy are not replicable in a more ancestrally homogeneous sample. We then investigate relationships between evolutionary annotations and common genetic variants shaping cortical surface area and white-matter connectivity for each hemisphere. Our analyses identify single-nucleotide polymorphism heritability enrichment in human-gained regulatory elements that are active in early brain development, affecting surface areas of several parts of the cortex, including left-hemispheric speech-associated regions. We also detect heritability depletion in genomic regions with Neanderthal ancestry for connectivity of the uncinate fasciculus; this is a white-matter tract involved in memory, language, and socioemotional processing with relevance to neuropsychiatric disorders. Finally, we show that common genetic loci associated with left-hemispheric pars triangularis surface area overlap with a human-gained enhancer and affect regulation of ZIC4, a gene implicated in neurogenesis. This work demonstrates how genomic investigations of present-day neuroanatomical variation can help shed light on the complexities of our evolutionary past.

Reliability of multi-site UK Biobank MRI brain phenotypes for the assessment of neuropsychiatric complications of SARS-CoV-2 infection: The COVID-CNS travelling heads study

INTRODUCTION

Magnetic resonance imaging (MRI) of the brain could be a key diagnostic and research tool for understanding the neuropsychiatric complications of COVID-19. For maximum impact, multi-modal MRI protocols will be needed to measure the effects of SARS-CoV-2 infection on the brain by diverse potentially pathogenic mechanisms, and with high reliability across multiple sites and scanner manufacturers. Here we describe the development of such a protocol, based upon the UK Biobank, and its validation with a travelling heads study. A multi-modal brain MRI protocol comprising sequences for T1-weighted MRI, T2-FLAIR, diffusion MRI (dMRI), resting-state functional MRI (fMRI), susceptibility-weighted imaging (swMRI), and arterial spin labelling (ASL), was defined in close approximation to prior UK Biobank (UKB) and C-MORE protocols for Siemens 3T systems. We iteratively defined a comparable set of sequences for General Electric (GE) 3T systems. To assess multi-site feasibility and between-site variability of this protocol, N = 8 healthy participants were each scanned at 4 UK sites: 3 using Siemens PRISMA scanners (Cambridge, Liverpool, Oxford) and 1 using a GE scanner (King's College London). Over 2,000 Imaging Derived Phenotypes (IDPs), measuring both data quality and regional image properties of interest, were automatically estimated by customised UKB image processing pipelines (S2 File). Components of variance and intra-class correlations (ICCs) were estimated for each IDP by linear mixed effects models and benchmarked by comparison to repeated measurements of the same IDPs from UKB participants. Intra-class correlations for many IDPs indicated good-to-excellent between-site reliability. Considering only data from the Siemens sites, between-site reliability generally matched the high levels of test-retest reliability of the same IDPs estimated in repeated, within-site, within-subject scans from UK Biobank. Inclusion of the GE site resulted in good-to-excellent reliability for many IDPs, although there were significant between-site differences in mean and scaling, and reduced ICCs, for some classes of IDP, especially T1 contrast and some dMRI-derived measures. We also identified high reliability of quantitative susceptibility mapping (QSM) IDPs derived from swMRI images, multi-network ICA-based IDPs from resting-state fMRI, and olfactory bulb structure IDPs from T1, T2-FLAIR and dMRI data.

CONCLUSION

These results give confidence that large, multi-site MRI datasets can be collected reliably at different sites across the diverse range of MRI modalities and IDPs that could be mechanistically informative in COVID brain research. We discuss limitations of the study and strategies for further harmonisation of data collected from sites using scanners supplied by different manufacturers. These acquisition and analysis protocols are now in use for MRI assessments of post-COVID patients (N = 700) as part of the ongoing COVID-CNS study.

Identifying genes associated with brain volumetric differences through tissue specific transcriptomic inference from GWAS summary data

BACKGROUND

Brain volume has been widely studied in the neuroimaging field, since it is an important and heritable trait associated with brain development, aging and various neurological and psychiatric disorders. Genome-wide association studies (GWAS) have successfully identified numerous associations between genetic variants such as single nucleotide polymorphisms and complex traits like brain volume. However, it is unclear how these genetic variations influence regional gene expression levels, which may subsequently lead to phenotypic changes. S-PrediXcan is a tissue-specific transcriptomic data analysis method that can be applied to bridge this gap. In this work, we perform an S-PrediXcan analysis on GWAS summary data from two large imaging genetics initiatives, the UK Biobank and Enhancing Neuroimaging Genetics through Meta Analysis, to identify tissue-specific transcriptomic effects on two closely related brain volume measures: total brain volume (TBV) and intracranial volume (ICV).

RESULTS

As a result of the analysis, we identified 10 genes that are highly associated with both TBV and ICV. Nine out of 10 genes were found to be associated with TBV in another study using a different gene-based association analysis. Moreover, most of our discovered genes were also found to be correlated with multiple cognitive and behavioral traits. Further analyses revealed the protein-protein interactions, associated molecular pathways and biological functions that offer insight into how these genes function and interact with others.

CONCLUSIONS

These results confirm that S-PrediXcan can identify genes with tissue-specific transcriptomic effects on complex traits. The analysis also suggested novel genes whose expression levels are related to brain volumetric traits. This provides important insights into the genetic mechanisms of the human brain.

Evidence of genetic overlap and causal relationships between blood-based biochemical traits and human cortical anatomy

Psychiatric disorders such as schizophrenia are commonly associated with structural brain alterations affecting the cortex. Recent genetic evidence suggests circulating metabolites and other biochemical traits play a causal role in many psychiatric disorders which could be mediated by changes in the cerebral cortex. Here, we leveraged publicly available genome-wide association study data to explore shared genetic architecture and evidence for causal relationships between a panel of 50 biochemical traits and measures of cortical thickness and surface area. Linkage disequilibrium score regression identified 191 genetically correlated biochemical-cortical trait pairings, with consistent representation of blood cell counts and other biomarkers such as C-reactive protein (CRP), haemoglobin and calcium. Spatially organised patterns of genetic correlation were additionally uncovered upon clustering of region-specific correlation profiles. Interestingly, by employing latent causal variable models, we found strong evidence suggesting CRP and vitamin D exert causal effects on region-specific cortical thickness, with univariable and multivariable Mendelian randomization further supporting a negative causal relationship between serum CRP levels and thickness of the lingual region. Our findings suggest a subset of biochemical traits exhibit shared genetic architecture and potentially causal relationships with cortical structure in functionally distinct regions, which may contribute to alteration of cortical structure in psychiatric disorders.

Integrative analysis of metabolite GWAS illuminates the molecular basis of pleiotropy and genetic correlation

Pleiotropy and genetic correlation are widespread features in genome-wide association studies (GWAS), but they are often difficult to interpret at the molecular level. Here, we perform GWAS of 16 metabolites clustered at the intersection of amino acid catabolism, glycolysis, and ketone body metabolism in a subset of UK Biobank. We utilize the well-documented biochemistry jointly impacting these metabolites to analyze pleiotropic effects in the context of their pathways. Among the 213 lead GWAS hits, we find a strong enrichment for genes encoding pathway-relevant enzymes and transporters. We demonstrate that the effect directions of variants acting on biology between metabolite pairs often contrast with those of upstream or downstream variants as well as the polygenic background. Thus, we find that these outlier variants often reflect biology local to the traits. Finally, we explore the implications for interpreting disease GWAS, underscoring the potential of unifying biochemistry with dense metabolomics data to understand the molecular basis of pleiotropy in complex traits and diseases.

Associations between brain imaging and polygenic scores of mental health and educational attainment in children aged 9-11

Psychiatric disorders are highly heritable and polygenic, and many have their peak onset in late childhood and adolescence, a period of tremendous changes. Although the neurodevelopmental antecedents of mental illness are widely acknowledged, research in youth population cohorts is still scarce, preventing our progress towards the early characterization of these disorders. We included 7,124 children (9-11 years old) from the Adolescent Brain and Cognitive Development Study to map the associations of structural and diffusion brain imaging with common genetic variants and polygenic scores for psychiatric disorders and educational attainment. We used principal component analysis to derive imaging components, and calculated their heritability. We then assessed the relationship of imaging components with genetic and clinical psychiatric risk with univariate models and Canonical correlation analysis (CCA). Most imaging components had moderate heritability. Univariate models showed limited evidence and small associations of polygenic scores with brain structure at this age. CCA revealed two significant modes of covariation. The first mode linked higher polygenic scores for educational attainment with less externalizing problems and larger surface area. The second mode related higher polygenic scores for schizophrenia, bipolar disorder, and autism spectrum disorder to higher global cortical thickness, smaller white matter volumes of the fornix and cingulum, larger medial occipital surface area and smaller surface area of lateral and medial temporal regions. While cross-validation suggested limited generalizability, our results highlight the potential of multivariate models to better understand the transdiagnostic and distributed relationships between mental health and brain structure in late childhood.

Causal inference on neuroimaging data with Mendelian randomisation

While population-scale neuroimaging studies offer the promise of discovery and characterisation of subtle risk factors, massive sample sizes increase the power for both meaningful associations and those attributable to confounds. This motivates the need for causal modelling of observational data that goes beyond statements of association and towards deeper understanding of complex relationships between individual traits and phenotypes, clinical biomarkers, genetic variation, and brain-related measures of health. Mendelian randomisation (MR) presents a way to obtain causal inference on the basis of genetic data and explicit assumptions about the relationship between genetic variables, exposure and outcome. In this work, we provide an introduction to and overview of causal inference methods based on Mendelian randomisation, with examples involving imaging-derived phenotypes from UK Biobank to make these methods accessible to neuroimaging researchers. We motivate the use of MR techniques, lay out the underlying assumptions, introduce common MR methods and focus on several scenarios in which modelling assumptions are potentially violated, resulting in biased effect estimates. Importantly, we give a detailed account of necessary steps to increase the reliability of MR results with rigorous sensitivity analyses.

The UK Biobank: A Shining Example of Genome-Wide Association Study Science with the Power to Detect the Murky Complications of Real-World Epidemiology

Genome-wide association studies (GWASs) have successfully identified thousands of genetic variants that are reliably associated with human traits. Although GWASs are restricted to certain variant frequencies, they have improved our understanding of the genetic architecture of complex traits and diseases. The UK Biobank (UKBB) has brought substantial analytical opportunity and performance to association studies. The dramatic expansion of many GWAS sample sizes afforded by the inclusion of UKBB data has improved the power of estimation of effect sizes but, critically, has done so in a context where phenotypic depth and precision enable outcome dissection and the application of epidemiological approaches. However, at the same time, the availability of such a large, well-curated, and deeply measured population-based collection has the capacity to increase our exposure to the many complications and inferential complexities associated with GWASs and other analyses. In this review, we discuss the impact that UKBB has had in the GWAS era, some of the opportunities that it brings, and exemplar challenges that illustrate the reality of using data from this world-leading resource.

Derivation and utility of schizophrenia polygenic risk associated multimodal MRI frontotemporal network

Schizophrenia is a highly heritable psychiatric disorder characterized by widespread functional and structural brain abnormalities. However, previous association studies between MRI and polygenic risk were mostly ROI-based single modality analyses, rather than identifying brain-based multimodal predictive biomarkers. Based on schizophrenia polygenic risk scores (PRS) from healthy white people within the UK Biobank dataset (N = 22,459), we discovered a robust PRS-associated brain pattern with smaller gray matter volume and decreased functional activation in frontotemporal cortex, which distinguished schizophrenia from controls with >83% accuracy, and predicted cognition and symptoms across 4 independent schizophrenia cohorts. Further multi-disease comparisons demonstrated that these identified frontotemporal alterations were most severe in schizophrenia and schizo-affective patients, milder in bipolar disorder, and indistinguishable from controls in autism, depression and attention-deficit hyperactivity disorder. These findings indicate the potential of the identified PRS-associated multimodal frontotemporal network to serve as a trans-diagnostic gene intermediated brain biomarker specific to schizophrenia.

Heritability of R2* iron in the basal ganglia and cortex

BACKGROUND

While iron is essential for normal brain functioning, elevated concentrations are commonly found in neurodegenerative diseases and are associated with impaired cognition and neurological deficits. Currently, only little is known about genetic and environmental factors that influence brain iron concentrations.

METHODS

Heritability and bivariate heritability of regional brain iron concentrations, assessed by R2* relaxometry at 3 Tesla MRI, were estimated with variance components models in 130 middle-aged to elderly participants of the Austrian Stroke Prevention Family Study.

RESULTS

Heritability of R2* iron ranged from 0.46 to 0.82 in basal ganglia and from 0.65 to 0.76 in cortical lobes. Age and BMI explained up to 12% and 9% of the variance of R2* iron, while APOE ε4 carrier status, hypertension, diabetes, hypercholesterolemia, sex and smoking explained 5% or less. The genetic correlation of R2* iron among basal ganglionic nuclei and among cortical lobes ranged from 0.78 to 0.87 and from 0.65 to 0.97, respectively. R2* rates in basal ganglia and cortex were not genetically correlated.

CONCLUSIONS

Regional brain iron concentrations are mainly driven by genetic factors while environmental factors contribute to a certain extent. Brain iron levels in the basal ganglia and cortex are controlled by distinct sets of genes.

A cognitive neurogenetic approach to uncovering the structure of executive functions

One central mission of cognitive neuroscience is to understand the ontology of complex cognitive functions. We addressed this question with a cognitive neurogenetic approach using a large-scale dataset of executive functions (EFs), whole-brain resting-state functional connectivity, and genetic polymorphisms. We found that the bifactor model with common and shifting-specific components not only was parsimonious but also showed maximal dissociations among the EF components at behavioral, neural, and genetic levels. In particular, the genes with enhanced expression in the middle frontal gyrus (MFG) and the subcallosal cingulate gyrus (SCG) showed enrichment for the common and shifting-specific component, respectively. Finally, High-dimensional mediation models further revealed that the functional connectivity patterns significantly mediated the genetic effect on the common EF component. Our study not only reveals insights into the ontology of EFs and their neurogenetic basis, but also provides useful tools to uncover the structure of complex constructs of human cognition.

Variability and concordance of sulcal patterns in the orbitofrontal cortex: A twin study

Sulcogyral patterns have been identified in the orbitofrontal cortex (OFC) based on the continuity of the medial and lateral orbital sulci. Pattern types are named according to their frequency in the population, with Type I present in ∼60%, Type II in ∼25%, Type III in ∼10%, and Type IV in ∼5%. Previous work has demonstrated that psychiatric conditions with high estimated heritability (e.g. schizophrenia, bipolar disorder) are associated with reduced frequency of Type I patterns, but the general heritability of the OFC sulcogyral patterns is unknown. We examined concordance of OFC patterns in 304 monozygotic (MZ) twins relative to 172 dizygotic (DZ) twins using structural magnetic resonance imaging data. We find that the frequency of pattern types within MZ and DZ twins are similar and bilateral concordance rates across all pattern types in DZ twins were 14% and 21% for MZ twins. Results from follow-up analyses confirm that continuity in the rostral-caudal direction is an important source of variability within the OFC, and subtype analyses indicate that variability is present in other sulci that are not represented by overall OFC pattern type. Overall, these results suggest that OFC sulcogyral patterns may reflect important variance that is not genetic in origin.

Variants of Th17 pathway-related genes influence brain morphometric changes and the risk of schizophrenia through epistatic interactions

OBJECTIVE

T helper 17 (Th17) pathway has been reported to be abnormal in schizophrenia; however, it is not known whether variation within genes of this pathway has any impact on schizophrenia. Herein, the impact of genetic variations and gene-gene interactions of Th17 pathway-related genes on the risk, psychopathology, and brain volume was examined in schizophrenia patients.

METHODS

Functional polymorphisms within interleukin 6 ( IL6 )(rs1800795 and rs1800797), IL10 (rs1800872 and rs1800896), IL17A (rs2275913 and rs8193036), IL22 (rs2227484 and rs2227485), IL23R (rs1884444), and IL27 (rs153109 and rs181206) genes were studied in 224 schizophrenia patients and 226 healthy controls. These variants were correlated with the brain morphometry, analyzed using MRI in a subset of patients ( n = 117) and controls ( n = 137).

RESULTS

Patients carrying CC genotype of rs2227484 of IL22 gene had significantly higher apathy total score [ F (1,183) = 5.60; P = 0.019; partial ɳ 2 = 0.030]. Significant epistatic interactions between IL6 (rs1800797) and IL17A (rs2275913) genes were observed in schizophrenia patients. GG genotype of rs2275913 of IL17A gene was associated with reduced right middle occipital gyrus volume in schizophrenia patients ( T = 4.56; P < 0.001).

CONCLUSION

Interactions between genes of Th17 pathway impact the risk for schizophrenia. The variants of Th17 pathway-related genes seem to have a determining effect on psychopathology and brain morphometric changes in schizophrenia.

transferGWAS: GWAS of images using deep transfer learning

MOTIVATION

Medical images can provide rich information about diseases and their biology. However, investigating their association with genetic variation requires non-standard methods. We propose transferGWAS, a novel approach to perform genome-wide association studies directly on full medical images. First, we learn semantically meaningful representations of the images based on a transfer learning task, during which a deep neural network is trained on independent but similar data. Then, we perform genetic association tests with these representations.

RESULTS

We validate the type I error rates and power of transferGWAS in simulation studies of synthetic images. Then we apply transferGWAS in a genome-wide association study of retinal fundus images from the UK Biobank. This first-of-a-kind GWAS of full imaging data yielded 60 genomic regions associated with retinal fundus images, of which 7 are novel candidate loci for eye-related traits and diseases.

AVAILABILITY AND IMPLEMENTATION

Our method is implemented in Python and available at https://github.com/mkirchler/transferGWAS/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Temporal mapping of derived high-frequency gene variants supports the mosaic nature of the evolution of Homo sapiens

Large-scale estimations of the time of emergence of variants are essential to examine hypotheses concerning human evolution with precision. Using an open repository of genetic variant age estimations, we offer here a temporal evaluation of various evolutionarily relevant datasets, such as Homo sapiens-specific variants, high-frequency variants found in genetic windows under positive selection, introgressed variants from extinct human species, as well as putative regulatory variants specific to various brain regions. We find a recurrent bimodal distribution of high-frequency variants, but also evidence for specific enrichments of gene categories in distinct time windows, pointing to different periods of phenotypic changes, resulting in a mosaic. With a temporal classification of genetic mutations in hand, we then applied a machine learning tool to predict what genes have changed more in certain time windows, and which tissues these genes may have impacted more. Overall, we provide a fine-grained temporal mapping of derived variants in Homo sapiens that helps to illuminate the intricate evolutionary history of our species.

Analysis of MRI-derived spleen iron in the UK Biobank identifies genetic variation linked to iron homeostasis and hemolysis

The spleen plays a key role in iron homeostasis. It is the largest filter of the blood and performs iron reuptake from old or damaged erythrocytes. Despite this role, spleen iron concentration has not been measured in a large, population-based cohort. In this study, we quantify spleen iron in 41,764 participants of the UK Biobank by using magnetic resonance imaging and provide a reference range for spleen iron in an unselected population. Through genome-wide association study, we identify associations between spleen iron and regulatory variation at two hereditary spherocytosis genes, ANK1 and SPTA1. Spherocytosis-causing coding mutations in these genes are associated with lower reticulocyte volume and increased reticulocyte percentage, while these common alleles are associated with increased expression of ANK1 and SPTA1 in blood and with larger reticulocyte volume and reduced reticulocyte percentage. As genetic modifiers, these common alleles may explain mild spherocytosis phenotypes that have been observed clinically. Our genetic study also identifies a signal that co-localizes with a splicing quantitative trait locus for MS4A7, and we show this gene is abundantly expressed in the spleen and in macrophages. The combination of deep learning and efficient image processing enables non-invasive measurement of spleen iron and, in turn, characterization of genetic factors related to the lytic phase of the erythrocyte life cycle and iron reuptake in the spleen.

Deep phenotyping for precision medicine in Parkinson's disease

A major challenge in medical genomics is to understand why individuals with the same disorder have different clinical symptoms and why those who carry the same mutation may be affected by different disorders. In every complex disorder, identifying the contribution of different genetic and non-genetic risk factors is a key obstacle to understanding disease mechanisms. Genetic studies rely on precise phenotypes and are unable to uncover the genetic contributions to a disorder when phenotypes are imprecise. To address this challenge, deeply phenotyped cohorts have been developed for which detailed, fine-grained data have been collected. These cohorts help us to investigate the underlying biological pathways and risk factors to identify treatment targets, and thus to advance precision medicine. The neurodegenerative disorder Parkinson's disease has a diverse phenotypical presentation and modest heritability, and its underlying disease mechanisms are still being debated. As such, considerable efforts have been made to develop deeply phenotyped cohorts for this disorder. Here, we focus on Parkinson's disease and explore how deep phenotyping can help address the challenges raised by genetic and phenotypic heterogeneity. We also discuss recent methods for data collection and computation, as well as methodological challenges that have to be overcome.

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.

Large-scale GWAS of food liking reveals genetic determinants and genetic correlations with distinct neurophysiological traits

We present the results of a GWAS of food liking conducted on 161,625 participants from the UK-Biobank. Liking was assessed over 139 specific foods using a 9-point scale. Genetic correlations coupled with structural equation modelling identified a multi-level hierarchical map of food-liking with three main dimensions: "Highly-palatable", "Acquired" and "Low-caloric". The Highly-palatable dimension is genetically uncorrelated from the other two, suggesting that independent processes underlie liking high reward foods. This is confirmed by genetic correlations with MRI brain traits which show with distinct associations. Comparison with the corresponding food consumption traits shows a high genetic correlation, while liking exhibits twice the heritability. GWAS analysis identified 1,401 significant food-liking associations which showed substantial agreement in the direction of effects with 11 independent cohorts. In conclusion, we created a comprehensive map of the genetic determinants and associated neurophysiological factors of food-liking.

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.

Machine Learning Evidence for Sex Differences Consistently Influences Resting-State Functional Magnetic Resonance Imaging Fluctuations Across Multiple Independently Acquired Data Sets

Background/Introduction: Sex classification using functional connectivity from resting-state functional magnetic resonance imaging (rs-fMRI) has shown promising results. This suggested that sex difference might also be embedded in the blood-oxygen-level-dependent properties such as the amplitude of low-frequency fluctuation (ALFF) and the fraction of ALFF (fALFF). This study comprehensively investigates sex differences using a reliable and explainable machine learning (ML) pipeline. Five independent cohorts of rs-fMRI with over than 5500 samples were used to assess sex classification performance and map the spatial distribution of the important brain regions. Methods: Five rs-fMRI samples were used to extract ALFF and fALFF features from predefined brain parcellations and then were fed into an unbiased and explainable ML pipeline with a wide range of methods. The pipeline comprehensively assessed unbiased performance for within-sample and across-sample validation. In addition, the parcellation effect, classifier selection, scanning length, spatial distribution, reproducibility, and feature importance were analyzed and evaluated thoroughly in the study. Results: The results demonstrated high sex classification accuracies from healthy adults (area under the curve >0.89), while degrading for nonhealthy subjects. Sex classification showed moderate to good intraclass correlation coefficient based on parcellation. Linear classifiers outperform nonlinear classifiers. Sex differences could be detected even with a short rs-fMRI scan (e.g., 2 min). The spatial distribution of important features overlaps with previous results from studies. Discussion: Sex differences are consistent in rs-fMRI and should be considered seriously in any study design, analysis, or interpretation. Features that discriminate males and females were found to be distributed across several different brain regions, suggesting a complex mosaic for sex differences in rs-fMRI. Impact statement The presented study unraveled that sex differences are embedded in the blood-oxygen-level dependent (BOLD) and can be predicted using unbiased and explainable machine learning pipeline. The study revealed that psychiatric disorders and demographics might influence the BOLD signal and interact with the classification of sex. The spatial distribution of the important features presented here supports the notion that the brain is a mosaic of male and female features. The findings emphasize the importance of controlling for sex when conducting brain imaging analysis. In addition, the presented framework can be adapted to classify other variables from resting-state BOLD signals.

Imaging Genetic Based Mediation Analysis for Human Cognition

The brain connectome maps the structural and functional connectivity that forms an important neurobiological basis for the analysis of human cognitive traits while the genetic predisposition and our cognition ability are frequently found in close association. The issue of how genetic architecture and brain connectome causally affect human behaviors remains unknown. To seek for the potential causal relationship, in this paper, we carried out the causal pathway analysis from single nucleotide polymorphism (SNP) data to four common human cognitive traits, mediated by the brain connectome. Specifically, we selected 942 SNPs that are significantly associated with the brain connectome, and then estimated the direct and indirect effect on the human traits for each SNP. We found out that a majority of the selected SNPs have significant direct effects on human traits and discussed the trait-related brain regions and their implications.

SARS-CoV-2 is associated with changes in brain structure in UK Biobank

There is strong evidence of brain-related abnormalities in COVID-191–13. However, it remains unknown whether the impact of SARS-CoV-2 infection can be detected in milder cases, and whether this can reveal possible mechanisms contributing to brain pathology. Here we investigated brain changes in 785 participants of UK Biobank (aged 51–81 years) who were imaged twice using magnetic resonance imaging, including 401 cases who tested positive for infection with SARS-CoV-2 between their two scans—with 141 days on average separating their diagnosis and the second scan—as well as 384 controls. The availability of pre-infection imaging data reduces the likelihood of pre-existing risk factors being misinterpreted as disease effects. We identified significant longitudinal effects when comparing the two groups, including (1) a greater reduction in grey matter thickness and tissue contrast in the orbitofrontal cortex and parahippocampal gyrus; (2) greater changes in markers of tissue damage in regions that are functionally connected to the primary olfactory cortex; and (3) a greater reduction in global brain size in the SARS-CoV-2 cases. The participants who were infected with SARS-CoV-2 also showed on average a greater cognitive decline between the two time points. Importantly, these imaging and cognitive longitudinal effects were still observed after excluding the 15 patients who had been hospitalised. These mainly limbic brain imaging results may be the in vivo hallmarks of a degenerative spread of the disease through olfactory pathways, of neuroinflammatory events, or of the loss of sensory input due to anosmia. Whether this deleterious effect can be partially reversed, or whether these effects will persist in the long term, remains to be investigated with additional follow-up.

Sensitivity of Diffusion MRI to White Matter Pathology: Influence of Diffusion Protocol, Magnetic Field Strength, and Processing Pipeline in Systemic Lupus Erythematosus

There are many ways to acquire and process diffusion MRI (dMRI) data for group studies, but it is unknown which maximizes the sensitivity to white matter (WM) pathology. Inspired by this question, we analyzed data acquired for diffusion tensor imaging (DTI) and diffusion kurtosis imaging (DKI) at 3T (3T-DTI and 3T-DKI) and DTI at 7T in patients with systemic lupus erythematosus (SLE) and healthy controls (HC). Parameter estimates in 72 WM tracts were obtained using TractSeg. The impact on the sensitivity to WM pathology was evaluated for the diffusion protocol, the magnetic field strength, and the processing pipeline. Sensitivity was quantified in terms of Cohen's d for group comparison. Results showed that the choice of diffusion protocol had the largest impact on the effect size. The effect size in fractional anisotropy (FA) across all WM tracts was 0.26 higher when derived by DTI than by DKI and 0.20 higher in 3T compared with 7T. The difference due to the diffusion protocol was larger than the difference due to magnetic field strength for the majority of diffusion parameters. In contrast, the difference between including or excluding different processing steps was near negligible, except for the correction of distortions from eddy currents and motion which had a clearly positive impact. For example, effect sizes increased on average by 0.07 by including motion and eddy correction for FA derived from 3T-DTI. Effect sizes were slightly reduced by the incorporation of denoising and Gibbs-ringing removal (on average by 0.011 and 0.005, respectively). Smoothing prior to diffusion model fitting generally reduced effect sizes. In summary, 3T-DTI in combination with eddy current and motion correction yielded the highest sensitivity to WM pathology in patients with SLE. However, our results also indicated that the 3T-DKI and 7T-DTI protocols used here may be adjusted to increase effect sizes.

Mapping Complex Brain Torque Components and Their Genetic Architecture and Phenomic Associations in 24,112 Individuals

BACKGROUND

The functional significance and mechanisms determining the development and individual variability of structural brain asymmetry remain unclear. Here, we systematically analyzed all relevant components of the most prominent structural asymmetry, brain torque (BT), and their relationships with potential genetic and nongenetic modifiers in a sample comprising 24,112 individuals from six cohorts.

METHODS

BT features, including petalia, bending, dorsoventral shift, brain tissue distribution asymmetries, and cortical surface positional asymmetries, were directly modeled using a set of automatic three-dimensional brain shape analysis approaches. Age-, sex-, and handedness-related effects on BT were assessed. The genetic architecture and phenomic associations of BT were investigated using genome- and phenome-wide association scans.

RESULTS

Our results confirmed the population-level predominance of the typical counterclockwise torque and suggested a first attenuating, then enlarging dynamic across the life span (3-81 years) primarily for frontal, occipital, and perisylvian BT features. Sex/handedness, BT, and cognitive function of verbal-numerical reasoning were found to be interrelated statistically. We observed differential heritability of up to 56% for BT, especially in temporal language areas. Individual variations of BT were also associated with various phenotypic variables of neuroanatomy, cognition, lifestyle, sociodemographics, anthropometry, physical health, and adult and child mental health. Our genomic analyses identified a number of genetic associations at lenient significance levels, which need to be further validated using larger samples in the future.

CONCLUSIONS

This study provides a comprehensive description of BT and insights into biological and other factors that may contribute to the development and individual variations of BT.

Using deep learning to predict abdominal age from liver and pancreas magnetic resonance images

With age, the prevalence of diseases such as fatty liver disease, cirrhosis, and type two diabetes increases. Approaches to both predict abdominal age and identify risk factors for accelerated abdominal age may ultimately lead to advances that will delay the onset of these diseases. We build an abdominal age predictor by training convolutional neural networks to predict abdominal age (or "AbdAge") from 45,552 liver magnetic resonance images [MRIs] and 36,784 pancreas MRIs (R-Squared = 73.3 ± 0.6; mean absolute error = 2.94 ± 0.03 years). Attention maps show that the prediction is driven by both liver and pancreas anatomical features, and surrounding organs and tissue. Abdominal aging is a complex trait, partially heritable (h_g² = 26.3 ± 1.9%), and associated with 16 genetic loci (e.g. in PLEKHA1 and EFEMP1), biomarkers (e.g body impedance), clinical phenotypes (e.g, chest pain), diseases (e.g. hypertension), environmental (e.g smoking), and socioeconomic (e.g education, income) factors.

Common variants contribute to intrinsic human brain functional networks

The human brain forms functional networks of correlated activity, which have been linked with both cognitive and clinical outcomes. However, the genetic variants affecting brain function are largely unknown. Here, we used resting-state functional magnetic resonance images from 47,276 individuals to discover and validate common genetic variants influencing intrinsic brain activity. We identified 45 new genetic regions associated with brain functional signatures (P < 2.8 × 10^-11), including associations to the central executive, default mode, and salience networks involved in the triple-network model of psychopathology. A number of brain activity-associated loci colocalized with brain disorders (e.g., the APOE ε4 locus with Alzheimer's disease). Variation in brain function was genetically correlated with brain disorders, such as major depressive disorder and schizophrenia. Together, our study provides a step forward in understanding the genetic architecture of brain functional networks and their genetic links to brain-related complex traits and disorders.

Genetic variants associated with longitudinal changes in brain structure across the lifespan

Human brain structure changes throughout the lifespan. Altered brain growth or rates of decline are implicated in a vast range of psychiatric, developmental and neurodegenerative diseases. In this study, we identified common genetic variants that affect rates of brain growth or atrophy in what is, to our knowledge, the first genome-wide association meta-analysis of changes in brain morphology across the lifespan. Longitudinal magnetic resonance imaging data from 15,640 individuals were used to compute rates of change for 15 brain structures. The most robustly identified genes GPR139, DACH1 and APOE are associated with metabolic processes. We demonstrate global genetic overlap with depression, schizophrenia, cognitive functioning, insomnia, height, body mass index and smoking. Gene set findings implicate both early brain development and neurodegenerative processes in the rates of brain changes. Identifying variants involved in structural brain changes may help to determine biological pathways underlying optimal and dysfunctional brain development and aging.

Brain charts for the human lifespan

Over the past few decades, neuroimaging has become a ubiquitous tool in basic research and clinical studies of the human brain. However, no reference standards currently exist to quantify individual differences in neuroimaging metrics over time, in contrast to growth charts for anthropometric traits such as height and weight1. Here we assemble an interactive open resource to benchmark brain morphology derived from any current or future sample of MRI data ( http://www.brainchart.io/ ). With the goal of basing these reference charts on the largest and most inclusive dataset available, acknowledging limitations due to known biases of MRI studies relative to the diversity of the global population, we aggregated 123,984 MRI scans, across more than 100 primary studies, from 101,457 human participants between 115 days post-conception to 100 years of age. MRI metrics were quantified by centile scores, relative to non-linear trajectories2 of brain structural changes, and rates of change, over the lifespan. Brain charts identified previously unreported neurodevelopmental milestones3, showed high stability of individuals across longitudinal assessments, and demonstrated robustness to technical and methodological differences between primary studies. Centile scores showed increased heritability compared with non-centiled MRI phenotypes, and provided a standardized measure of atypical brain structure that revealed patterns of neuroanatomical variation across neurological and psychiatric disorders. In summary, brain charts are an essential step towards robust quantification of individual variation benchmarked to normative trajectories in multiple, commonly used neuroimaging phenotypes.

SARS-CoV-2 is associated with changes in brain structure in UK Biobank

There is strong evidence of brain-related abnormalities in COVID-191-13. However, it remains unknown whether the impact of SARS-CoV-2 infection can be detected in milder cases, and whether this can reveal possible mechanisms contributing to brain pathology. Here we investigated brain changes in 785 participants of UK Biobank (aged 51-81 years) who were imaged twice using magnetic resonance imaging, including 401 cases who tested positive for infection with SARS-CoV-2 between their two scans-with 141 days on average separating their diagnosis and the second scan-as well as 384 controls. The availability of pre-infection imaging data reduces the likelihood of pre-existing risk factors being misinterpreted as disease effects. We identified significant longitudinal effects when comparing the two groups, including (1) a greater reduction in grey matter thickness and tissue contrast in the orbitofrontal cortex and parahippocampal gyrus; (2) greater changes in markers of tissue damage in regions that are functionally connected to the primary olfactory cortex; and (3) a greater reduction in global brain size in the SARS-CoV-2 cases. The participants who were infected with SARS-CoV-2 also showed on average a greater cognitive decline between the two time points. Importantly, these imaging and cognitive longitudinal effects were still observed after excluding the 15 patients who had been hospitalised. These mainly limbic brain imaging results may be the in vivo hallmarks of a degenerative spread of the disease through olfactory pathways, of neuroinflammatory events, or of the loss of sensory input due to anosmia. Whether this deleterious effect can be partially reversed, or whether these effects will persist in the long term, remains to be investigated with additional follow-up.

Genetic common variants associated with cerebellar volume and their overlap with mental disorders: a study on 33,265 individuals from the UK-Biobank

Interest in the cerebellum is expanding given evidence of its contributions to cognition and emotion, and dysfunction in various psychopathologies. However, research into its genetic architecture and shared influences with liability for mental disorders is lacking. We conducted a genome-wide association study (GWAS) of total cerebellar volume and underlying cerebellar lobe volumes in 33,265 UK-Biobank participants. Total cerebellar volume was heritable (h2SNP = 50.6%), showing moderate genetic homogeneity across lobes (h2SNP from 35.4% to 57.1%; mean genetic correlation between lobes rg ≈ 0.44). We identified 33 GWAS signals associated with total cerebellar volume, of which 6 are known to alter protein-coding gene structure, while a further five mapped to genomic regions known to alter cerebellar tissue gene expression. Use of summary data-based Mendelian randomisation further prioritised genes whose change in expression appears to mediate the SNP-trait association. In total, we highlight 21 unique genes of greatest interest for follow-up analyses. Using LD-regression, we report significant genetic correlations between total cerebellar volume and brainstem, pallidum and thalamus volumes. While the same approach did not result in significant correlations with psychiatric phenotypes, we report enrichment of schizophrenia, bipolar disorder and autism spectrum disorder associated signals within total cerebellar GWAS results via conditional and conjunctional-FDR analysis. Via these methods and GWAS catalogue, we identify which of our cerebellar genomic regions also associate with psychiatric traits. Our results provide important insights into the common allele architecture of cerebellar volume and its overlap with other brain volumes and psychiatric phenotypes.

Caltech Conte Center, a multimodal data resource for exploring social cognition and decision-making

This data release of 117 healthy community-dwelling adults provides multimodal high-quality neuroimaging and behavioral data for the investigation of brain-behavior relationships. We provide structural MRI, resting-state functional MRI, movie functional MRI, together with questionnaire-based and task-based psychological variables; many of the participants have multiple datasets from retesting over the course of several years. Our dataset is distinguished by utilizing open-source data formats and processing tools (BIDS, FreeSurfer, fMRIPrep, MRIQC), providing data that is thoroughly quality checked, preprocessed to various extents and available in multiple anatomical spaces. A customizable denoising pipeline is provided as open-source code that includes tools for the generation of functional connectivity matrices and initialization of individual difference analyses. Behavioral data include a comprehensive set of psychological assessments on gold-standard instruments encompassing cognitive function, mood and personality, together with exploratory factor analyses. The dataset provides an in-depth, multimodal resource for investigating associations between individual differences, brain structure and function, with a focus on the domains of social cognition and decision-making.

Total Brain Volumetric Measures and Schizophrenia Risk: A Two-Sample Mendelian Randomization Study

Schizophrenia (SCZ) is an idiopathic psychiatric disorder with a heritable component and a substantial public health impact. Although abnormalities in total brain volumetric measures (TBVMs) have been found in patients with SCZ, it is still unknown whether these abnormalities have a causal effect on the risk of SCZ. Here, we performed a Mendelian randomization (MR) study to investigate the possible causal associations between each TBVM and SCZ risk. Specifically, genome-wide association study (GWAS) summary statistics of total gray matter volume, total white matter volume, total cerebrospinal fluid volume, and total brain volume were obtained from the United Kingdom Biobank database (33,224 individuals), and SCZ GWAS summary statistics were provided by the Psychiatric Genomics Consortium (150,064 individuals). The main MR analysis was conducted using the inverse variance weighted method, and other MR methods, including MR-Egger, weighted median, simple mode, and weighted mode methods, were performed to assess the robustness of our findings. For pleiotropy analysis, we employed three approaches: MR-Egger intercept, MR-PRESSO, and heterogeneity tests. No TBVM was causally associated with SCZ risk according to the MR results, and no significant pleiotropy or heterogeneity was found for instrumental variables. Taken together, this study suggested that alterations in TBVMs were not causally associated with the risk of SCZ.

Exploring the genetic correlation between obesity-related traits and regional brain volumes: Evidence from UK Biobank cohort

OBJECTIVE

To determine whether there is a correlation between obesity-related variants and regional brain volumes.

METHODS

Based on a mixed linear model (MLM), we analyzed the association between 1,498 obesity-related SNPs in the GWAS Catalog and 164 regional brain volumes from 29,420 participants (discovery cohort N = 19,997, validation cohort N = 9,423) in UK Biobank. The statistically significant brain regions in association analysis were classified into 6 major neural networks (dopamine (DA) motive system, central autonomic network (CAN), cognitive emotion regulation, visual object recognition network, auditory object recognition network, and sensorimotor system). We summarized the association between obesity-related variants (metabolically healthy obesity variants, metabolically unhealthy obesity variants, and unclassified obesity-related variants) and neural networks.

RESULTS

From association analysis, we determined that 17 obesity-related SNPs were associated with 51 regional brain volumes. Several single SNPs (e.g., rs13107325-T (SLC39A8), rs1876829-C (CRHR1), and rs1538170-T (CENPW)) were associated with multiple regional brain volumes. In addition, several single brain regions (e.g., the white matter, the grey matter in the putamen, subcallosal cortex, and insular cortex) were associated with multiple obesity-related variants. The metabolically healthy obesity variants were mainly associated with the regional brain volumes in the DA motive system, sensorimotor system and cognitive emotion regulation neural networks, while metabolically unhealthy obesity variants were mainly associated with regional brain volumes in the CAN and total tissue volumes. In addition, unclassified obesity-related variants were mainly associated with auditory object recognition network and total tissue volumes. The results of MeSH (medical subject headings) enrichment analysis showed that obesity genes associated with brain structure pointed to the functional relatedness with 5-Hydroxytryptamine receptor 4 (5-HT4), growth differentiation factor 5 (GDF5), and high mobility group protein AT-hook 2 (HMGA2 protein).

CONCLUSION

In summary, we found that obesity-related variants were associated with different brain volume measures. On the basis of the multiple SNPs, we found that metabolically healthy and unhealthy obesity-related SNPs were associated with different brain neural networks. Based on our enrichment analysis, modifications of the 5-HT4 pathway might be a promising therapeutic strategy for obesity.

SARS-CoV-2 is associated with changes in brain structure in UK Biobank

There is strong evidence for brain-related abnormalities in COVID-19 1-13 . It remains unknown however whether the impact of SARS-CoV-2 infection can be detected in milder cases, and whether this can reveal possible mechanisms contributing to brain pathology. Here, we investigated brain changes in 785 UK Biobank participants (aged 51-81) imaged twice, including 401 cases who tested positive for infection with SARS-CoV-2 between their two scans, with 141 days on average separating their diagnosis and second scan, and 384 controls. The availability of pre-infection imaging data reduces the likelihood of pre-existing risk factors being misinterpreted as disease effects. We identified significant longitudinal effects when comparing the two groups, including: (i) greater reduction in grey matter thickness and tissue-contrast in the orbitofrontal cortex and parahippocampal gyrus, (ii) greater changes in markers of tissue damage in regions functionally-connected to the primary olfactory cortex, and (iii) greater reduction in global brain size. The infected participants also showed on average larger cognitive decline between the two timepoints. Importantly, these imaging and cognitive longitudinal effects were still seen after excluding the 15 cases who had been hospitalised. These mainly limbic brain imaging results may be the in vivo hallmarks of a degenerative spread of the disease via olfactory pathways, of neuroinflammatory events, or of the loss of sensory input due to anosmia. Whether this deleterious impact can be partially reversed, or whether these effects will persist in the long term, remains to be investigated with additional follow up.

Shared Genetic Etiology between Cortical Brain Morphology and Tobacco, Alcohol, and Cannabis Use

Genome-wide association studies (GWAS) have identified genetic variants associated with brain morphology and substance use behaviors (SUB). However, the genetic overlap between brain structure and SUB has not been well characterized. We leveraged GWAS summary data of 71 brain imaging measures and alcohol, tobacco, and cannabis use to investigate their genetic overlap using linkage disequilibrium score regression. We used genomic structural equation modeling to model a "common SUB genetic factor" and investigated its genetic overlap with brain structure. Furthermore, we estimated SUB polygenic risk scores (PRS) and examined whether they predicted brain imaging traits using the Adolescent Behavior and Cognitive Development (ABCD) study. We identified 8 significant negative genetic correlations, including between (1) alcoholic drinks per week and average cortical thickness, and (2) intracranial volume with age of smoking initiation. We observed 5 positive genetic correlations, including those between (1) insula surface area and lifetime cannabis use, and (2) the common SUB genetic factor and pericalcarine surface area. SUB PRS were associated with brain structure variation in ABCD. Our findings highlight a shared genetic etiology between cortical brain morphology and SUB and suggest that genetic variants associated with SUB may be causally related to brain structure differences.

Genetic Influence Underlying Brain Connectivity Phenotype: A Study on Two Age-Specific Cohorts

Background: Human brain structural connectivity is an important imaging quantitative trait for brain development and aging. Mapping the network connectivity to the phenotypic variation provides fundamental insights in understanding the relationship between detailed brain topological architecture, function, and dysfunction. However, the underlying neurobiological mechanism from gene to brain connectome, and to phenotypic outcomes, and whether this mechanism changes over time, remain unclear. Methods: This study analyzes diffusion-weighted imaging data from two age-specific neuroimaging cohorts, extracts structural connectome topological network measures, performs genome-wide association studies of the measures, and examines the causality of genetic influences on phenotypic outcomes mediated via connectivity measures. Results: Our empirical study has yielded several significant findings: 1) It identified genetic makeup underlying structural connectivity changes in the human brain connectome for both age groups. Specifically, it revealed a novel association between the minor allele (G) of rs7937515 and the decreased network segregation measures of the left middle temporal gyrus across young and elderly adults, indicating a consistent genetic effect on brain connectivity across the lifespan. 2) It revealed rs7937515 as a genetic marker for body mass index in young adults but not in elderly adults. 3) It discovered brain network segregation alterations as a potential neuroimaging biomarker for obesity. 4) It demonstrated the hemispheric asymmetry of structural network organization in genetic association analyses and outcome-relevant studies. Discussion: These imaging genetic findings underlying brain connectome warrant further investigation for exploring their potential influences on brain-related complex diseases, given the significant involvement of altered connectivity in neurological, psychiatric and physical disorders.

Neuron-derived neuropeptide Y fine-tunes the splenic immune responses

The nervous and immune systems are closely entwined to maintain the immune balance in health and disease. Here, we showed that LPS can activate suprarenal and celiac ganglia (SrG-CG) neurons and upregulate NPY expression in rats. Single-cell sequencing analysis revealed that knockdown of the NPY gene in SrG-CG altered the proliferation and activation of splenic lymphocytes. In a neuron and splenocyte coculture system and in vivo experiments, neuronal NPY in SrG-CG attenuated the splenic immune response. Notably, we demonstrated that neuronal NPF in Drosophila exerted a conservative immunomodulatory effect. Moreover, numerous SNPs in NPY and its receptors were significantly associated with human autoimmune diseases, which was further supported by the autoimmune disease patients and mouse model experiments. Together, we demonstrated that NPY is an ancient language for nervous-immune system crosstalk and might be utilized to alleviate inflammatory storms during infection and to modulate immune balance in autoimmune diseases.

Discovery of genomic loci of the human cerebral cortex using genetically informed brain atlases

To determine the impact of genetic variants on the brain, we used genetically informed brain atlases in genome-wide association studies of regional cortical surface area and thickness in 39,898 adults and 9136 children. We uncovered 440 genome-wide significant loci in the discovery cohort and 800 from a post hoc combined meta-analysis. Loci in adulthood were largely captured in childhood, showing signatures of negative selection, and were linked to early neurodevelopment and pathways associated with neuropsychiatric risk. Opposing gradations of decreased surface area and increased thickness were associated with common inversion polymorphisms. Inferior frontal regions, encompassing Broca's area, which is important for speech, were enriched for human-specific genomic elements. Thus, a mixed genetic landscape of conserved and human-specific features is concordant with brain hierarchy and morphogenetic gradients.

Mental health in the UK Biobank: A roadmap to self-report measures and neuroimaging correlates

The UK Biobank (UKB) is a highly promising dataset for brain biomarker research into population mental health due to its unprecedented sample size and extensive phenotypic, imaging, and biological measurements. In this study, we aimed to provide a shared foundation for UKB neuroimaging research into mental health with a focus on anxiety and depression. We compared UKB self-report measures and revealed important timing effects between scan acquisition and separate online acquisition of some mental health measures. To overcome these timing effects, we introduced and validated the Recent Depressive Symptoms (RDS-4) score which we recommend for state-dependent and longitudinal research in the UKB. We furthermore tested univariate and multivariate associations between brain imaging-derived phenotypes (IDPs) and mental health. Our results showed a significant multivariate relationship between IDPs and mental health, which was replicable. Conversely, effect sizes for individual IDPs were small. Test-retest reliability of IDPs was stronger for measures of brain structure than for measures of brain function. Taken together, these results provide benchmarks and guidelines for future UKB research into brain biomarkers of mental health.