Genetic architecture of subcortical brain structures in 38,851 individuals

Increased brain volume from higher cereal and lower coffee intake: shared genetic determinants and impacts on cognition and metabolism

It is unclear how different diets may affect human brain development and if genetic and environmental factors play a part. We investigated diet effects in the UK Biobank data from 18,879 healthy adults and discovered anticorrelated brain-wide gray matter volume (GMV)-association patterns between coffee and cereal intake, coincidence with their anticorrelated genetic constructs. The Mendelian randomization approach further indicated a causal effect of higher coffee intake on reduced total GMV, which is likely through regulating the expression of genes responsible for synaptic development in the brain. The identified genetic factors may further affect people's lifestyle habits and body/blood fat levels through the mediation of cereal/coffee intake, and the brain-wide expression pattern of gene CPLX3, a dedicated marker of subplate neurons that regulate cortical development and plasticity, may underlie the shared GMV-association patterns among the coffee/cereal intake and cognitive functions. All the main findings were successfully replicated. Our findings thus revealed that high-cereal and low-coffee diets shared similar brain and genetic constructs, leading to long-term beneficial associations regarding cognitive, body mass index (BMI), and other metabolic measures. This study has important implications for public health, especially during the pandemic, given the poorer outcomes of COVID-19 patients with greater BMIs.

Integrating Transcriptomics, Genomics, and Imaging in Alzheimer's Disease: A Federated Model

Alzheimer's disease (AD) affects more than 1 in 9 people age 65 and older and becomes an urgent public health concern as the global population ages. In clinical practice, structural magnetic resonance imaging (sMRI) is the most accessible and widely used diagnostic imaging modality. Additionally, genome-wide association studies (GWAS) and transcriptomics-the study of gene expression-also play an important role in understanding AD etiology and progression. Sophisticated imaging genetics systems have been developed to discover genetic factors that consistently affect brain function and structure. However, most studies to date focused on the relationships between brain sMRI and GWAS or brain sMRI and transcriptomics. To our knowledge, few methods have been developed to discover and infer multimodal relationships among sMRI, GWAS, and transcriptomics. To address this, we propose a novel federated model, Genotype-Expression-Imaging Data Integration (GEIDI), to identify genetic and transcriptomic influences on brain sMRI measures. The relationships between brain imaging measures and gene expression are allowed to depend on a person's genotype at the single-nucleotide polymorphism (SNP) level, making the inferences adaptive and personalized. We performed extensive experiments on publicly available Alzheimer's Disease Neuroimaging Initiative (ADNI) dataset. Experimental results demonstrated our proposed method outperformed state-of-the-art expression quantitative trait loci (eQTL) methods for detecting genetic and transcriptomic factors related to AD and has stable performance when data are integrated from multiple sites. Our GEIDI approach may offer novel insights into the relationship among image biomarkers, genotypes, and gene expression and help discover novel genetic targets for potential AD drug treatments.

What we learn about bipolar disorder from large-scale neuroimaging: Findings and future directions from the ENIGMA Bipolar Disorder Working Group

MRI-derived brain measures offer a link between genes, the environment and behavior and have been widely studied in bipolar disorder (BD). However, many neuroimaging studies of BD have been underpowered, leading to varied results and uncertainty regarding effects. The Enhancing Neuro Imaging Genetics through Meta-Analysis (ENIGMA) Bipolar Disorder Working Group was formed in 2012 to empower discoveries, generate consensus findings and inform future hypothesis-driven studies of BD. Through this effort, over 150 researchers from 20 countries and 55 institutions pool data and resources to produce the largest neuroimaging studies of BD ever conducted. The ENIGMA Bipolar Disorder Working Group applies standardized processing and analysis techniques to empower large-scale meta- and mega-analyses of multimodal brain MRI and improve the replicability of studies relating brain variation to clinical and genetic data. Initial BD Working Group studies reveal widespread patterns of lower cortical thickness, subcortical volume and disrupted white matter integrity associated with BD. Findings also include mapping brain alterations of common medications like lithium, symptom patterns and clinical risk profiles and have provided further insights into the pathophysiological mechanisms of BD. Here we discuss key findings from the BD working group, its ongoing projects and future directions for large-scale, collaborative studies of mental illness.

Enhancing neuroimaging genetics through meta-analysis for Tourette syndrome (ENIGMA-TS): A worldwide platform for collaboration

Tourette syndrome (TS) is characterized by multiple motor and vocal tics, and high-comorbidity rates with other neuropsychiatric disorders. Obsessive compulsive disorder (OCD), attention deficit hyperactivity disorder (ADHD), autism spectrum disorders (ASDs), major depressive disorder (MDD), and anxiety disorders (AXDs) are among the most prevalent TS comorbidities. To date, studies on TS brain structure and function have been limited in size with efforts mostly fragmented. This leads to low-statistical power, discordant results due to differences in approaches, and hinders the ability to stratify patients according to clinical parameters and investigate comorbidity patterns. Here, we present the scientific premise, perspectives, and key goals that have motivated the establishment of the Enhancing Neuroimaging Genetics through Meta-Analysis for TS (ENIGMA-TS) working group. The ENIGMA-TS working group is an international collaborative effort bringing together a large network of investigators who aim to understand brain structure and function in TS and dissect the underlying neurobiology that leads to observed comorbidity patterns and clinical heterogeneity. Previously collected TS neuroimaging data will be analyzed jointly and integrated with TS genomic data, as well as equivalently large and already existing studies of highly comorbid OCD, ADHD, ASD, MDD, and AXD. Our work highlights the power of collaborative efforts and transdiagnostic approaches, and points to the existence of different TS subtypes. ENIGMA-TS will offer large-scale, high-powered studies that will lead to important insights toward understanding brain structure and function and genetic effects in TS and related disorders, and the identification of biomarkers that could help inform improved clinical practice.

Effects of copy number variations on brain structure and risk for psychiatric illness: Large-scale studies from the ENIGMA working groups on CNVs

The Enhancing NeuroImaging Genetics through Meta-Analysis copy number variant (ENIGMA-CNV) and 22q11.2 Deletion Syndrome Working Groups (22q-ENIGMA WGs) were created to gain insight into the involvement of genetic factors in human brain development and related cognitive, psychiatric and behavioral manifestations. To that end, the ENIGMA-CNV WG has collated CNV and magnetic resonance imaging (MRI) data from ~49,000 individuals across 38 global research sites, yielding one of the largest studies to date on the effects of CNVs on brain structures in the general population. The 22q-ENIGMA WG includes 12 international research centers that assessed over 533 individuals with a confirmed 22q11.2 deletion syndrome, 40 with 22q11.2 duplications, and 333 typically developing controls, creating the largest-ever 22q11.2 CNV neuroimaging data set. In this review, we outline the ENIGMA infrastructure and procedures for multi-site analysis of CNVs and MRI data. So far, ENIGMA has identified effects of the 22q11.2, 16p11.2 distal, 15q11.2, and 1q21.1 distal CNVs on subcortical and cortical brain structures. Each CNV is associated with differences in cognitive, neurodevelopmental and neuropsychiatric traits, with characteristic patterns of brain structural abnormalities. Evidence of gene-dosage effects on distinct brain regions also emerged, providing further insight into genotype-phenotype relationships. Taken together, these results offer a more comprehensive picture of molecular mechanisms involved in typical and atypical brain development. This "genotype-first" approach also contributes to our understanding of the etiopathogenesis of brain disorders. Finally, we outline future directions to better understand effects of CNVs on brain structure and behavior.

Sensitive period-regulating genetic pathways and exposure to adversity shape risk for depression

Animal and human studies have documented the existence of developmental windows (or sensitive periods) when experience can have lasting effects on brain structure or function, behavior, and disease. Although sensitive periods for depression likely arise through a complex interplay of genes and experience, this possibility has not yet been explored in humans. We examined the effect of genetic pathways regulating sensitive periods, alone and in interaction with common childhood adversities, on depression risk. Guided by a translational approach, we: (1) performed association analyses of three gene sets (60 genes) shown in animal studies to regulate sensitive periods using summary data from a genome-wide association study of depression (n = 807,553); (2) evaluated the developmental expression patterns of these genes using data from BrainSpan (n = 31), a transcriptional atlas of postmortem brain samples; and (3) tested gene-by-development interplay (dGxE) by analyzing the combined effect of common variants in sensitive period genes and time-varying exposure to two types of childhood adversity within a population-based birth cohort (n = 6254). The gene set regulating sensitive period opening associated with increased depression risk. Notably, 6 of the 15 genes in this set showed developmentally regulated gene-level expression. We also identified a statistical interaction between caregiver physical or emotional abuse during ages 1-5 years and genetic risk for depression conferred by the opening genes. Genes involved in regulating sensitive periods are differentially expressed across the life course and may be implicated in depression vulnerability. Our findings about gene-by-development interplay motivate further research in large, more diverse samples to further unravel the complexity of depression etiology through a sensitive period lens.

An overview of the first 5 years of the ENIGMA obsessive-compulsive disorder working group: The power of worldwide collaboration

Neuroimaging has played an important part in advancing our understanding of the neurobiology of obsessive-compulsive disorder (OCD). At the same time, neuroimaging studies of OCD have had notable limitations, including reliance on relatively small samples. International collaborative efforts to increase statistical power by combining samples from across sites have been bolstered by the ENIGMA consortium; this provides specific technical expertise for conducting multi-site analyses, as well as access to a collaborative community of neuroimaging scientists. In this article, we outline the background to, development of, and initial findings from ENIGMA's OCD working group, which currently consists of 47 samples from 34 institutes in 15 countries on 5 continents, with a total sample of 2,323 OCD patients and 2,325 healthy controls. Initial work has focused on studies of cortical thickness and subcortical volumes, structural connectivity, and brain lateralization in children, adolescents and adults with OCD, also including the study on the commonalities and distinctions across different neurodevelopment disorders. Additional work is ongoing, employing machine learning techniques. Findings to date have contributed to the development of neurobiological models of OCD, have provided an important model of global scientific collaboration, and have had a number of clinical implications. Importantly, our work has shed new light on questions about whether structural and functional alterations found in OCD reflect neurodevelopmental changes, effects of the disease process, or medication impacts. We conclude with a summary of ongoing work by ENIGMA-OCD, and a consideration of future directions for neuroimaging research on OCD within and beyond ENIGMA.

Evaluating brain structure traits as endophenotypes using polygenicity and discoverability

Human brain structure traits have been hypothesized to be broad endophenotypes for neuropsychiatric disorders, implying that brain structure traits are comparatively "closer to the underlying biology." Genome-wide association studies from large sample sizes allow for the comparison of common variant genetic architectures between traits to test the evidence supporting this claim. Endophenotypes, compared to neuropsychiatric disorders, are hypothesized to have less polygenicity, with greater effect size of each susceptible SNP, requiring smaller sample sizes to discover them. Here, we compare polygenicity and discoverability of brain structure traits, neuropsychiatric disorders, and other traits (91 in total) to directly test this hypothesis. We found reduced polygenicity (FDR = 0.01) and increased discoverability (FDR = 3.68 × 10-9 ) of cortical brain structure traits, as compared to aggregated estimates of multiple neuropsychiatric disorders. We predict that ~8 M individuals will be required to explain the full heritability of cortical surface area by genome-wide significant SNPs, whereas sample sizes over 20 M will be required to explain the full heritability of depression. In conclusion, our findings are consistent with brain structure satisfying the higher power criterion of endophenotypes.

A meta-analysis of deep brain structural shape and asymmetry abnormalities in 2,833 individuals with schizophrenia compared with 3,929 healthy volunteers via the ENIGMA Consortium

Schizophrenia is associated with widespread alterations in subcortical brain structure. While analytic methods have enabled more detailed morphometric characterization, findings are often equivocal. In this meta-analysis, we employed the harmonized ENIGMA shape analysis protocols to collaboratively investigate subcortical brain structure shape differences between individuals with schizophrenia and healthy control participants. The study analyzed data from 2,833 individuals with schizophrenia and 3,929 healthy control participants contributed by 21 worldwide research groups participating in the ENIGMA Schizophrenia Working Group. Harmonized shape analysis protocols were applied to each site's data independently for bilateral hippocampus, amygdala, caudate, accumbens, putamen, pallidum, and thalamus obtained from T1-weighted structural MRI scans. Mass univariate meta-analyses revealed more-concave-than-convex shape differences in the hippocampus, amygdala, accumbens, and thalamus in individuals with schizophrenia compared with control participants, more-convex-than-concave shape differences in the putamen and pallidum, and both concave and convex shape differences in the caudate. Patterns of exaggerated asymmetry were observed across the hippocampus, amygdala, and thalamus in individuals with schizophrenia compared to control participants, while diminished asymmetry encompassed ventral striatum and ventral and dorsal thalamus. Our analyses also revealed that higher chlorpromazine dose equivalents and increased positive symptom levels were associated with patterns of contiguous convex shape differences across multiple subcortical structures. Findings from our shape meta-analysis suggest that common neurobiological mechanisms may contribute to gray matter reduction across multiple subcortical regions, thus enhancing our understanding of the nature of network disorganization in schizophrenia.

ENIGMA-anxiety working group: Rationale for and organization of large-scale neuroimaging studies of anxiety disorders

Anxiety disorders are highly prevalent and disabling but seem particularly tractable to investigation with translational neuroscience methodologies. Neuroimaging has informed our understanding of the neurobiology of anxiety disorders, but research has been limited by small sample sizes and low statistical power, as well as heterogenous imaging methodology. The ENIGMA-Anxiety Working Group has brought together researchers from around the world, in a harmonized and coordinated effort to address these challenges and generate more robust and reproducible findings. This paper elaborates on the concepts and methods informing the work of the working group to date, and describes the initial approach of the four subgroups studying generalized anxiety disorder, panic disorder, social anxiety disorder, and specific phobia. At present, the ENIGMA-Anxiety database contains information about more than 100 unique samples, from 16 countries and 59 institutes. Future directions include examining additional imaging modalities, integrating imaging and genetic data, and collaborating with other ENIGMA working groups. The ENIGMA consortium creates synergy at the intersection of global mental health and clinical neuroscience, and the ENIGMA-Anxiety Working Group extends the promise of this approach to neuroimaging research on anxiety disorders.

Mega-analysis methods in ENIGMA: The experience of the generalized anxiety disorder working group

The ENIGMA group on Generalized Anxiety Disorder (ENIGMA-Anxiety/GAD) is part of a broader effort to investigate anxiety disorders using imaging and genetic data across multiple sites worldwide. The group is actively conducting a mega-analysis of a large number of brain structural scans. In this process, the group was confronted with many methodological challenges related to study planning and implementation, between-country transfer of subject-level data, quality control of a considerable amount of imaging data, and choices related to statistical methods and efficient use of resources. This report summarizes the background information and rationale for the various methodological decisions, as well as the approach taken to implement them. The goal is to document the approach and help guide other research groups working with large brain imaging data sets as they develop their own analytic pipelines for mega-analyses.

Thirty novel sequence variants impacting human intracranial volume

Intracranial volume, measured through magnetic resonance imaging and/or estimated from head circumference, is heritable and correlates with cognitive traits and several neurological disorders. We performed a genome-wide association study meta-analysis of intracranial volume (n = 79 174) and found 64 associating sequence variants explaining 5.0% of its variance. We used coding variation, transcript and protein levels, to uncover 12 genes likely mediating the effect of these variants, including GLI3 and CDK6 that affect cranial synostosis and microcephaly, respectively. Intracranial volume correlates genetically with volumes of cortical and sub-cortical regions, cognition, learning, neonatal and neurological traits. Parkinson's disease cases have greater and attention deficit hyperactivity disorder cases smaller intracranial volume than controls. Our Mendelian randomization studies indicate that intracranial volume associated variants either increase the risk of Parkinson's disease and decrease the risk of attention deficit hyperactivity disorder and neuroticism or correlate closely with a confounder.

Genetic basis to structural grey matter associations with chronic pain

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

SHIP-MR and Radiology: 12 Years of Whole-Body Magnetic Resonance Imaging in a Single Center

The Study of Health in Pomerania (SHIP), a population-based study from a rural state in northeastern Germany with a relatively poor life expectancy, supplemented its comprehensive examination program in 2008 with whole-body MR imaging at 1.5 T (SHIP-MR). We reviewed more than 100 publications that used the SHIP-MR data and analyzed which sequences already produced fruitful scientific outputs and which manuscripts have been referenced frequently. Upon reviewing the publications about imaging sequences, those that used T1-weighted structured imaging of the brain and a gradient-echo sequence for R2* mapping obtained the highest scientific output; regarding specific body parts examined, most scientific publications focused on MR sequences involving the brain and the (upper) abdomen. We conclude that population-based MR imaging in cohort studies should define more precise goals when allocating imaging time. In addition, quality control measures might include recording the number and impact of published work, preferably on a bi-annual basis and starting 2 years after initiation of the study. Structured teaching courses may enhance the desired output in areas that appear underrepresented.

Trauma and posttraumatic stress disorder modulate polygenic predictors of hippocampal and amygdala volume

The volume of subcortical structures represents a reliable, quantitative, and objective phenotype that captures genetic effects, environmental effects such as trauma, and disease effects such as posttraumatic stress disorder (PTSD). Trauma and PTSD represent potent exposures that may interact with genetic markers to influence brain structure and function. Genetic variants, associated with subcortical volumes in two large normative discovery samples, were used to compute polygenic scores (PGS) for the volume of seven subcortical structures. These were applied to a target sample enriched for childhood trauma and PTSD. Subcortical volume PGS from the discovery sample were strongly associated in our trauma/PTSD enriched sample (n = 7580) with respective subcortical volumes of the hippocampus (p = 1.10 × 10-20), thalamus (p = 7.46 × 10-10), caudate (p = 1.97 × 10-18), putamen (p = 1.7 × 10-12), and nucleus accumbens (p = 1.99 × 10-7). We found a significant association between the hippocampal volume PGS and hippocampal volume in control subjects from our sample, but was absent in individuals with PTSD (GxE; (beta = -0.10, p = 0.027)). This significant GxE (PGS × PTSD) relationship persisted (p < 1 × 10-19) in four out of five threshold peaks (0.024, 0.133, 0.487, 0.730, and 0.889) used to calculate hippocampal volume PGSs. We detected similar GxE (G × ChildTrauma) relationships in the amygdala for exposure to childhood trauma (rs4702973; p = 2.16 × 10-7) or PTSD (rs10861272; p = 1.78 × 10-6) in the CHST11 gene. The hippocampus and amygdala are pivotal brain structures in mediating PTSD symptomatology. Trauma exposure and PTSD modulate the effect of polygenic markers on hippocampal volume (GxE) and the amygdala volume PGS is associated with PTSD risk, which supports the role of amygdala volume as a risk factor for PTSD.

Integrated Genomic, Transcriptomic and Proteomic Analysis for Identifying Markers of Alzheimer's Disease

There is an urgent need to identify biomarkers for Alzheimer’s disease (AD), but the identification of reliable blood-based biomarkers has proven to be much more difficult than initially expected. The current availability of high-throughput multi-omics data opens new possibilities in this titanic task. Candidate Single Nucleotide Polymorphisms (SNPs) from large, genome-wide association studies (GWAS), meta-analyses exploring AD (case-control design), and quantitative measures for cortical structure and general cognitive performance were selected. The Genotype-Tissue Expression (GTEx) database was used for identifying expression quantitative trait loci (eQTls) among candidate SNPs. Genes significantly regulated by candidate SNPs were investigated for differential expression in AD cases versus controls in the brain and plasma, both at the mRNA and protein level. This approach allowed us to identify candidate susceptibility factors and biomarkers of AD, facing experimental validation with more evidence than with genetics alone.

Grey and white matter microstructure is associated with polygenic risk for schizophrenia

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

Where the genome meets the connectome: Understanding how genes shape human brain connectivity

The integration of modern neuroimaging methods with genetically informative designs and data can shed light on the molecular mechanisms underlying the structural and functional organization of the human connectome. Here, we review studies that have investigated the genetic basis of human brain network structure and function through three complementary frameworks: (1) the quantification of phenotypic heritability through classical twin designs; (2) the identification of specific DNA variants linked to phenotypic variation through association and related studies; and (3) the analysis of correlations between spatial variations in imaging phenotypes and gene expression profiles through the integration of neuroimaging and transcriptional atlas data. We consider the basic foundations, strengths, limitations, and discoveries associated with each approach. We present converging evidence to indicate that anatomical connectivity is under stronger genetic influence than functional connectivity and that genetic influences are not uniformly distributed throughout the brain, with phenotypic variation in certain regions and connections being under stronger genetic control than others. We also consider how the combination of imaging and genetics can be used to understand the ways in which genes may drive brain dysfunction in different clinical disorders.

Behavioral, Anatomical and Heritable Convergence of Affect and Cognition in Superior Frontal Cortex

Cognitive abilities and affective experience are key human traits that are interrelated in behavior and brain. Individual variation of cognitive and affective traits, as well as brain structure, has been shown to partly underlie genetic effects. However, to what extent affect and cognition have a shared genetic relationship with local brain structure is incompletely understood. Here we studied phenotypic and genetic correlations of cognitive and affective traits in behavior and brain structure (cortical thickness, surface area and subcortical volumes) in the pedigree-based Human Connectome Project sample (N = 1091). Both cognitive and affective trait scores were highly heritable and showed significant phenotypic correlation on the behavioral level. Cortical thickness in the left superior frontal cortex showed a phenotypic association with both affect and cognition. Decomposing the phenotypic correlations into genetic and environmental components showed that the associations were accounted for by shared genetic effects between the traits. Quantitative functional decoding of the left superior frontal cortex further indicated that this region is associated with cognitive and emotional functioning. This study provides a multi-level approach to study the association between affect and cognition and suggests a convergence of both in superior frontal cortical thickness.

The role of glucocorticoid receptor gene in the association between attention deficit-hyperactivity disorder and smaller brain structures

ADHD is associated with smaller subcortical brain volumes and cortical surface area, with greater effects observed in children than adults. It is also associated with dysregulation of the HPA axis. Considering the effects of the glucocorticoid receptor (NR3C1) in neurophysiology, we hypothesize that the blurred relationships between brain structures and ADHD in adults could be partly explained by NR3C1 gene variation. Structural T1-weighted images were acquired on a 3 T scanner (N = 166). Large-scale genotyping was performed, and it was followed by quality control and pruning procedures, which resulted in 48 independent NR3C1 gene variants analyzed. After a stringent Bonferroni correction, two SNPs (rs2398631 and rs72801070) moderated the association between ADHD and accumbens and amygdala volumes in adults. The significant SNPs that interacted with ADHD appear to have a role in gene expression regulation, and they are in linkage disequilibrium with NR3C1 variants that present well-characterized physiological functions. The literature-reported associations of ADHD with accumbens and amygdala were only observed for specific NR3C1 genotypes. Our findings reinforce the influence of the NR3C1 gene on subcortical volumes and ADHD. They suggest a genetic modulation of the effects of a pivotal HPA axis component in the neuroanatomical features of ADHD.

The Enhancing NeuroImaging Genetics through Meta-Analysis Consortium: 10 Years of Global Collaborations in Human Brain Mapping

This Special Issue of Human Brain Mapping is dedicated to a 10-year anniversary of the Enhancing NeuroImaging Genetics through Meta-Analysis (ENIGMA) Consortium. It reports updates from a broad range of international neuroimaging projects that pool data from around the world to answer fundamental questions in neuroscience. Since ENIGMA was formed in December 2009, the initiative grew into a worldwide effort with over 2,000 participating scientists from 45 countries, and over 50 working groups leading large-scale studies of human brain disorders. Over the last decade, many lessons were learned on how best to pool brain data from diverse sources. Working groups were created to develop methods to analyze worldwide data from anatomical and diffusion magnetic resonance imaging (MRI), resting state and task-based functional MRI, electroencephalography (EEG), magnetoencephalography (MEG), and magnetic resonance spectroscopy (MRS). The quest to understand genetic effects on human brain development and disease also led to analyses of brain scans on an unprecedented scale. Genetic roadmaps of the human cortex were created by researchers worldwide who collaborated to perform statistically well-powered analyses of common and rare genetic variants on brain measures and rates of brain development and aging. Here, we summarize the 31 papers in this Special Issue, covering: (a) technical approaches to harmonize analysis of different types of brain imaging data, (b) reviews of the last decade of work by several of ENIGMA's clinical and technical working groups, and (c) new empirical papers reporting large-scale international brain mapping analyses in patients with substance use disorders, schizophrenia, bipolar disorders, major depression, posttraumatic stress disorder, obsessive compulsive disorder, epilepsy, and stroke.

Twin studies to GWAS: there and back again

The field of human behavioral genetics has come full circle. It began by using twin/family studies to estimate the relative importance of genetic and environmental influences. As large-scale genotyping became cost-effective, genome-wide association studies (GWASs) yielded insights about the nature of genetic influences and new methods that use GWAS data to estimate heritability and genetic correlations invigorated the field. Yet these newer GWAS methods have not replaced twin/family studies. In this review, we discuss the strengths and weaknesses of the two approaches with respect to characterizing genetic and environmental influences, measurement of behavioral phenotypes, and evaluation of causal models, with a particular focus on cognitive neuroscience. This discussion highlights how twin/family studies and GWAS complement and mutually reinforce one another.

Imaging Transcriptomics of Brain Disorders

Noninvasive neuroimaging is a powerful tool for quantifying diverse aspects of brain structure and function in vivo, and it has been used extensively to map the neural changes associated with various brain disorders. However, most neuroimaging techniques offer only indirect measures of underlying pathological mechanisms. The recent development of anatomically comprehensive gene expression atlases has opened new opportunities for studying the transcriptional correlates of noninvasively measured neural phenotypes, offering a rich framework for evaluating pathophysiological hypotheses and putative mechanisms. Here, we provide an overview of some fundamental methods in imaging transcriptomics and outline their application to understanding brain disorders of neurodevelopment, adulthood, and neurodegeneration. Converging evidence indicates that spatial variations in gene expression are linked to normative changes in brain structure during age-related maturation and neurodegeneration that are in part associated with cell-specific gene expression markers of gene expression. Transcriptional correlates of disorder-related neuroimaging phenotypes are also linked to transcriptionally dysregulated genes identified in ex vivo analyses of patient brains. Modeling studies demonstrate that spatial patterns of gene expression are involved in regional vulnerability to neurodegeneration and the spread of disease across the brain. This growing body of work supports the utility of transcriptional atlases in testing hypotheses about the molecular mechanism driving disease-related changes in macroscopic neuroimaging phenotypes.

Polygenic Architecture of Human Neuroanatomical Diversity

We analyzed the genomic architecture of neuroanatomical diversity using magnetic resonance imaging and single nucleotide polymorphism (SNP) data from >26 000 individuals from the UK Biobank project and 5 other projects that had previously participated in the ENIGMA (Enhancing NeuroImaging Genetics through Meta-Analysis) consortium. Our results confirm the polygenic architecture of neuroanatomical diversity, with SNPs capturing from 40% to 54% of regional brain volume variance. Chromosomal length correlated with the amount of phenotypic variance captured, r ~ 0.64 on average, suggesting that at a global scale causal variants are homogeneously distributed across the genome. At a local scale, SNPs within genes (~51%) captured ~1.5 times more genetic variance than the rest, and SNPs with low minor allele frequency (MAF) captured less variance than the rest: the 40% of SNPs with MAF <5% captured <one fourth of the genetic variance. We also observed extensive pleiotropy across regions, with an average genetic correlation of rG ~ 0.45. Genetic correlations were similar to phenotypic and environmental correlations; however, genetic correlations were often larger than phenotypic correlations for the left/right volumes of the same region. The heritability of differences in left/right volumes was generally not statistically significant, suggesting an important influence of environmental causes in the variability of brain asymmetry. Our code is available athttps://github.com/neuroanatomy/genomic-architecture.

Sex differences in brain atrophy and cognitive impairment in Parkinson's disease patients with and without probable rapid eye movement sleep behavior disorder

Background

The presence of rapid eye movement sleep behavior disorder (RBD) contributes to increase cognitive impairment and brain atrophy in Parkinson’s disease (PD), but the impact of sex is unclear. We aimed to investigate sex differences in cognition and brain atrophy in PD patients with and without probable RBD (pRBD).

Methods

Magnetic resonance imaging and cognition data were obtained for 274 participants from the Parkinson's Progression Marker Initiative database: 79 PD with pRBD (PD-pRBD; male/female, 54/25), 126 PD without pRBD (PD-non pRBD; male/female, 73/53), and 69 healthy controls (male/female, 40/29). FreeSurfer was used to obtain volumetric and cortical thickness data.

Results

Males showed greater global cortical and subcortical gray matter atrophy than females in the PD-pRBD group. Significant group-by-sex interactions were found in the pallidum. Structures showing a within-group sex effect in the deep gray matter differed, with significant volume reductions for males in one structure in in PD-non pRBD (brainstem), and three in PD-pRBD (caudate, pallidum and brainstem). Significant group-by-sex interactions were found in Montreal Cognitive Assessment (MoCA) and Symbol Digits Modalities Test (SDMT). Males performed worse than females in MoCA, phonemic fluency and SDMT in the PD-pRBD group.

Conclusion

Male sex is related to increased cognitive impairment and subcortical atrophy in de novo PD-pRBD. Accordingly, we suggest that sex differences are relevant and should be considered in future clinical and translational research.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00415-021-10728-x.

The genetic organization of longitudinal subcortical volumetric change is stable throughout the lifespan

Development and aging of the cerebral cortex show similar topographic organization and are governed by the same genes. It is unclear whether the same is true for subcortical regions, which follow fundamentally different ontogenetic and phylogenetic principles. We tested the hypothesis that genetically governed neurodevelopmental processes can be traced throughout life by assessing to which degree brain regions that develop together continue to change together through life. Analyzing over 6000 longitudinal MRIs of the brain, we used graph theory to identify five clusters of coordinated development, indexed as patterns of correlated volumetric change in brain structures. The clusters tended to follow placement along the cranial axis in embryonic brain development, suggesting continuity from prenatal stages, and correlated with cognition. Across independent longitudinal datasets, we demonstrated that developmental clusters were conserved through life. Twin-based genetic correlations revealed distinct sets of genes governing change in each cluster. Single-nucleotide polymorphisms-based analyses of 38,127 cross-sectional MRIs showed a similar pattern of genetic volume–volume correlations. In conclusion, coordination of subcortical change adheres to fundamental principles of lifespan continuity and genetic organization.

[A review on the application of UK Biobank in neuroimaging]

UK Biobank (UKB) is a forward-looking epidemiological project with over 500, 000 people aged 40 to 69, whose image extension project plans to re-invite 100, 000 participants from UKB to perform multimodal brain magnetic resonance imaging. Large-scale multimodal neuroimaging combined with large amounts of phenotypic and genetic data provides great resources to conduct brain health-related research. This article provides an in-depth overview of UKB in the field of neuroimaging. Firstly, neuroimage collection and imaging-derived phenotypes are summarized. Secondly, typical studies of UKB in neuroimaging areas are introduced, which include cardiovascular risk factors, regulatory factors, brain age prediction, normality, successful and morbid brain aging, environmental and genetic factors, cognitive ability and gender. Lastly, the open challenges and future directions of UKB are discussed. This article has the potential to open up a new research field for the prevention and treatment of neurological diseases.

Ionotropic Receptors as a Driving Force behind Human Synapse Establishment

The origin of nervous systems is a main theme in biology and its mechanisms are largely underlied by synaptic neurotransmission. One problem to explain synapse establishment is that synaptic orthologs are present in multiple aneural organisms. We questioned how the interactions among these elements evolved and to what extent it relates to our understanding of the nervous systems complexity. We identified the human neurotransmission gene network based on genes present in GABAergic, glutamatergic, serotonergic, dopaminergic, and cholinergic systems. The network comprises 321 human genes, 83 of which act exclusively in the nervous system. We reconstructed the evolutionary scenario of synapse emergence by looking for synaptic orthologs in 476 eukaryotes. The Human–Cnidaria common ancestor displayed a massive emergence of neuroexclusive genes, mainly ionotropic receptors, which might have been crucial to the evolution of synapses. Very few synaptic genes had their origin after the Human–Cnidaria common ancestor. We also identified a higher abundance of synaptic proteins in vertebrates, which suggests an increase in the synaptic network complexity of those organisms.

Shared heritability of human face and brain shape

Evidence from model organisms and clinical genetics suggests coordination between the developing brain and face, but the role of this link in common genetic variation remains unknown. We performed a multivariate genome-wide association study of cortical surface morphology in 19,644 individuals of European ancestry, identifying 472 genomic loci influencing brain shape, of which 76 are also linked to face shape. Shared loci include transcription factors involved in craniofacial development, as well as members of signaling pathways implicated in brain-face cross-talk. Brain shape heritability is equivalently enriched near regulatory regions active in either forebrain organoids or facial progenitors. However, we do not detect significant overlap between shared brain-face genome-wide association study signals and variants affecting behavioral-cognitive traits. These results suggest that early in embryogenesis, the face and brain mutually shape each other through both structural effects and paracrine signaling, but this interplay may not impact later brain development associated with cognitive function.

The Drosophila ortholog of the schizophrenia-associated CACNA1A and CACNA1B voltage-gated calcium channels regulate memory, sleep and circadian rhythms

Schizophrenia exhibits up to 80% heritability. A number of genome wide association studies (GWAS) have repeatedly shown common variants in voltage-gated calcium (Ca_v) channel genes CACNA1C, CACNA1I and CACNA1G have a major contribution to the risk of the disease. More recently, studies using whole exome sequencing have also found that CACNA1B (Ca_v2.2 N-type) deletions and rare disruptive variants in CACNA1A (Ca_v2.1 P/Q-type) are associated with schizophrenia. The negative symptoms of schizophrenia include behavioural defects such as impaired memory, sleep and circadian rhythms. It is not known how variants in schizophrenia-associated genes contribute to cognitive and behavioural symptoms, thus hampering the development of treatment for schizophrenia symptoms. In order to address this knowledge gap, we studied behavioural phenotypes in a number of loss of function mutants for the Drosophila ortholog of the Ca_v2 gene family called cacophony (cac). cac mutants showed several behavioural features including decreased night-time sleep and hyperactivity similar to those reported in human patients. The change in timing of sleep-wake cycles suggested disrupted circadian rhythms, with the loss of night-time sleep being caused by loss of cac just in the circadian clock neurons. These animals also showed a reduction in rhythmic circadian behaviour a phenotype that also could be mapped to the central clock. Furthermore, reduction of cac just in the clock resulted in a lengthening of the 24 h period. In order to understand how loss of Ca_v2 function may lead to cognitive deficits and underlying cellular pathophysiology we targeted loss of function of cac to the memory centre of the fly, called the mushroom bodies (MB). This manipulation was sufficient to cause reduction in both short- and intermediate-term associative memory. Memory impairment was accompanied by a decrease in Ca²⁺ transients in response to a depolarizing stimulus, imaged in the MB presynaptic terminals. This work shows loss of cac Ca_v2 channel function alone causes a number of cognitive and behavioural deficits and underlying reduced neuronal Ca²⁺ transients, establishing Drosophila as a high-throughput in vivo genetic model to study the Ca_v channel pathophysiology related to schizophrenia.

The genetic architecture of the human thalamus and its overlap with ten common brain disorders

The thalamus is a vital communication hub in the center of the brain and consists of distinct nuclei critical for consciousness and higher-order cortical functions. Structural and functional thalamic alterations are involved in the pathogenesis of common brain disorders, yet the genetic architecture of the thalamus remains largely unknown. Here, using brain scans and genotype data from 30,114 individuals, we identify 55 lead single nucleotide polymorphisms (SNPs) within 42 genetic loci and 391 genes associated with volumes of the thalamus and its nuclei. In an independent validation sample (n = 5173) 53 out of the 55 lead SNPs of the discovery sample show the same effect direction (sign test, P = 8.6e-14). We map the genetic relationship between thalamic nuclei and 180 cerebral cortical areas and find overlapping genetic architectures consistent with thalamocortical connectivity. Pleiotropy analyses between thalamic volumes and ten psychiatric and neurological disorders reveal shared variants for all disorders. Together, these analyses identify genetic loci linked to thalamic nuclei and substantiate the emerging view of the thalamus having central roles in cortical functioning and common brain disorders.

Transcriptome-wide association analysis of brain structures yields insights into pleiotropy with complex neuropsychiatric traits

Structural variations of the human brain are heritable and highly polygenic traits, with hundreds of associated genes identified in recent genome-wide association studies (GWAS). Transcriptome-wide association studies (TWAS) can both prioritize these GWAS findings and also identify additional gene-trait associations. Here we perform cross-tissue TWAS analysis of 211 structural neuroimaging and discover 278 associated genes exceeding Bonferroni significance threshold of 1.04 × 10-8. The TWAS-significant genes for brain structures have been linked to a wide range of complex traits in different domains. Through TWAS gene-based polygenic risk scores (PRS) prediction, we find that TWAS PRS gains substantial power in association analysis compared to conventional variant-based GWAS PRS, and up to 6.97% of phenotypic variance (p-value = 7.56 × 10-31) can be explained in independent testing data sets. In conclusion, our study illustrates that TWAS can be a powerful supplement to traditional GWAS in imaging genetics studies for gene discovery-validation, genetic co-architecture analysis, and polygenic risk prediction.

The morphometry of left cuneus mediating the genetic regulation on working memory

Working memory is a basic human cognitive function. However, the genetic signatures and their biological pathway remain poorly understood. In the present study, we tried to clarify this issue by exploring the potential associations and pathways among genetic variants, brain morphometry and working memory performance. We first carried out association analyses between 2‐back accuracy and 212 image‐derived phenotypes from 1141 Human Connectome Project (HCP) subjects using a linear mixed model (LMM). We found a significantly positive correlation between the left cuneus volume and 2‐back accuracy (T = 3.615, p = 3.150e−4, Cohen's d = 0.226, corrected using family‐wise error [FWE] method). Based on the LMM‐based genome‐wide association study (GWAS) on the HCP dataset and UK Biobank 33 k GWAS summary statistics, we identified eight independent single nucleotide polymorphisms (SNPs) that were reliably associated with left cuneus volume in both UKB and HCP dataset. Within the eight SNPs, we found a negative correlation between the rs76119478 polymorphism and 2‐back accuracy accuracy (T = −2.045, p = .041, Cohen's d = −0.129). Finally, an LMM‐based mediation analysis elucidated a significant effect of left cuneus volume in mediating rs76119478 polymorphism on the 2‐back accuracy (indirect effect = −0.007, 95% BCa CI = [−0.045, −0.003]). These results were also replicated in a subgroup of Caucasians in the HCP population. Further fine mapping demonstrated that rs76119478 maps on intergene CTD‐2315A10.2 adjacent to protein‐encoding gene DAAM1, and is significantly associated with L3HYPDH mRNA expression. Our study suggested this new variant rs76119478 may regulate the working memory through exerting influence on the left cuneus volume.

Amyloid toxicity in a Drosophila Alzheimer's model is ameliorated by autophagy activation

Alzheimer's disease (AD) is the prevailing form of dementia. Protein degradation and antioxidant pathways have a critical role in preventing the accumulation of protein aggregation; thus, failure of proteostasis in neurons along with redox imbalance mark AD. Herein, we exploited an AD Drosophila model expressing human amyloid precursor (hAPP) and beta-secretase 1 (hBACE1) proteins, to better understand the role of proteostatic or antioxidant pathways in AD. Ubiquitous expression of hAPP, hBACE1 in flies caused more severe degenerative phenotypes versus neuronal targeted expression; it also, suppressed proteasome activity, increased oxidative stress and significantly enhanced stress-sensitivity. Overexpression of Prosβ5 proteasomal subunit or Nrf2 transcription factor in AD Drosophila flies partially restored proteasomal activity but did not rescue hAPP, hBACE1 induced neurodegeneration. On the other hand, expression of autophagy-related Atg8a in AD flies decelerated neurodegeneration, increased stress-resistance, and improved flies' health-/lifespan. Overall, our data suggest that the noxious effects of amyloid-beta aggregates can be alleviated by enhanced autophagy, thus dietary or pharmacological interventions that target autophagy should be considered in AD therapeutic approaches.

Genetic scores for adult subcortical volumes associate with subcortical volumes during infancy and childhood

Individual differences in subcortical brain volumes are highly heritable. Previous studies have identified genetic variants that underlie variation in subcortical volumes in adults. We tested whether those previously identified variants also affect subcortical regions during infancy and early childhood. The study was performed within the Generation R study, a prospective birth cohort. We calculated polygenic scores based on reported GWAS for volumes of the accumbens, amygdala, brainstem, caudate nucleus, globus pallidus, putamen, and thalamus. Participants underwent cranial ultrasound around 7 weeks of age (range: 3-20), and we obtained metrics for the gangliothalamic ovoid, a predecessor of the basal ganglia. Furthermore, the children participated in a magnetic resonance imaging (MRI) study around the age of 10 years (range: 9-12). A total of 340 children had complete data at both examinations. Polygenic scores primarily associated with their corresponding volumes at 10 years of age. The scores also moderately related to the diameter of the gangliothalamic ovoid on cranial ultrasound. Mediation analysis showed that the genetic influence on subcortical volumes at 10 years was only mediated for 16.5-17.6% of the total effect through the gangliothalamic ovoid diameter at 7 weeks of age. Combined, these findings suggest that previously identified genetic variants in adults are relevant for subcortical volumes during early life, and that they affect both prenatal and postnatal development of the subcortical regions.

Comorbid Chronic Pain and Depression: Shared Risk Factors and Differential Antidepressant Effectiveness

The bidirectional relationship between depression and chronic pain is well-recognized, but their clinical management remains challenging. Here we characterize the shared risk factors and outcomes for their comorbidity in the Australian Genetics of Depression cohort study (N = 13,839). Participants completed online questionnaires about chronic pain, psychiatric symptoms, comorbidities, treatment response and general health. Logistic regression models were used to examine the relationship between chronic pain and clinical and demographic factors. Cumulative linked logistic regressions assessed the effect of chronic pain on treatment response for 10 different antidepressants. Chronic pain was associated with an increased risk of depression (OR = 1.86 [1.37-2.54]), recent suicide attempt (OR = 1.88 [1.14-3.09]), higher use of tobacco (OR = 1.05 [1.02-1.09]) and misuse of painkillers (e.g., opioids; OR = 1.31 [1.06-1.62]). Participants with comorbid chronic pain and depression reported fewer functional benefits from antidepressant use and lower benefits from sertraline (OR = 0.75 [0.68-0.83]), escitalopram (OR = 0.75 [0.67-0.85]) and venlafaxine (OR = 0.78 [0.68-0.88]) when compared to participants without chronic pain. Furthermore, participants taking sertraline (OR = 0.45 [0.30-0.67]), escitalopram (OR = 0.45 [0.27-0.74]) and citalopram (OR = 0.32 [0.15-0.67]) specifically for chronic pain (among other indications) reported lower benefits compared to other participants taking these same medications but not for chronic pain. These findings reveal novel insights into the complex relationship between chronic pain and depression. Treatment response analyses indicate differential effectiveness between particular antidepressants and poorer functional outcomes for these comorbid conditions. Further examination is warranted in targeted interventional clinical trials, which also include neuroimaging genetics and pharmacogenomics protocols. This work will advance the delineation of disease risk indicators and novel aetiological pathways for therapeutic intervention in comorbid pain and depression as well as other psychiatric comorbidities.

Likelihood ratio statistics for gene set enrichment in Alzheimer's disease pathways

INTRODUCTION

The study of Alzheimer's disease (AD) has revealed biological pathways with implications for disease neuropathology and pathophysiology. These pathway-level effects may also be mediated by individual characteristics or covariates such as age or sex. Evaluation of AD biological pathways in the context of interactions with these covariates is critical to the understanding of AD as well as the development of model systems used to study the disease.

METHODS

Gene set enrichment methods are powerful tools used to interpret gene-level statistics at the level of biological pathways. We introduce a method for quantifying gene set enrichment using likelihood ratio-derived test statistics (gsLRT), which accounts for sample covariates like age and sex. We then use our method to test for age and sex interactions with protein expression levels in AD and to compare the pathway results between human and mouse species.

RESULTS

Our method, based on nested logistic regressions is competitive with the existing standard for gene set testing in the context of linear models and complex experimental design. The gene sets we identify as having a significant association with AD-both with and without additional covariate interactions-are validated by previous studies. Differences between gsLRT results on mouse and human datasets are observed.

DISCUSSION

Characterizing biological pathways involved in AD builds on the important work involving single gene drivers. Our gene set enrichment method finds pathways that are significantly related to AD while accounting for covariates that may be relevant to disease development. The method highlights commonalities and differences between human AD and mouse models, which may inform the development of higher fidelity models for the study of AD.

Sex-stratified genome-wide association study of multisite chronic pain in UK Biobank

Chronic pain is highly prevalent worldwide and imparts a significant socioeconomic and public health burden. Factors influencing susceptibility to, and mechanisms of, chronic pain development, are not fully understood, but sex is thought to play a significant role, and chronic pain is more prevalent in women than in men. To investigate sex differences in chronic pain, we carried out a sex-stratified genome-wide association study of Multisite Chronic Pain (MCP), a derived chronic pain phenotype, in UK Biobank on 178,556 men and 209,093 women, as well as investigating sex-specific genetic correlations with a range of psychiatric, autoimmune and anthropometric phenotypes and the relationship between sex-specific polygenic risk scores for MCP and chronic widespread pain. We also assessed whether MCP-associated genes showed expression pattern enrichment across tissues. A total of 123 SNPs at five independent loci were significantly associated with MCP in men. In women, a total of 286 genome-wide significant SNPs at ten independent loci were discovered. Meta-analysis of sex-stratified GWAS outputs revealed a further 87 independent associated SNPs. Gene-level analyses revealed sex-specific MCP associations, with 31 genes significantly associated in females, 37 genes associated in males, and a single gene, DCC, associated in both sexes. We found evidence for sex-specific pleiotropy and risk for MCP was found to be associated with chronic widespread pain in a sex-differential manner. Male and female MCP were highly genetically correlated, but at an rg of significantly less than 1 (0.92). All 37 male MCP-associated genes and all but one of 31 female MCP-associated genes were found to be expressed in the dorsal root ganglion, and there was a degree of enrichment for expression in sex-specific tissues. Overall, the findings indicate that sex differences in chronic pain exist at the SNP, gene and transcript abundance level, and highlight possible sex-specific pleiotropy for MCP. Results support the proposition of a strong central nervous-system component to chronic pain in both sexes, additionally highlighting a potential role for the DRG and nociception.

Genetic factors influencing a neurobiological substrate for psychiatric disorders

A retrospective meta-analysis of magnetic resonance imaging voxel-based morphometry studies proposed that reduced gray matter volumes in the dorsal anterior cingulate and the left and right anterior insular cortex-areas that constitute hub nodes of the salience network-represent a common substrate for major psychiatric disorders. Here, we investigated the hypothesis that the common substrate serves as an intermediate phenotype to detect genetic risk variants relevant for psychiatric disease. To this end, after a data reduction step, we conducted genome-wide association studies of a combined common substrate measure in four population-based cohorts (n = 2271), followed by meta-analysis and replication in a fifth cohort (n = 865). After correction for covariates, the heritability of the common substrate was estimated at 0.50 (standard error 0.18). The top single-nucleotide polymorphism (SNP) rs17076061 was associated with the common substrate at genome-wide significance and replicated, explaining 1.2% of the common substrate variance. This SNP mapped to a locus on chromosome 5q35.2 harboring genes involved in neuronal development and regeneration. In follow-up analyses, rs17076061 was not robustly associated with psychiatric disease, and no overlap was found between the broader genetic architecture of the common substrate and genetic risk for major depressive disorder, bipolar disorder, or schizophrenia. In conclusion, our study identified that common genetic variation indeed influences the common substrate, but that these variants do not directly translate to increased disease risk. Future studies should investigate gene-by-environment interactions and employ functional imaging to understand how salience network structure translates to psychiatric disorder risk.

Effects of Thyroid Status on Regional Brain Volumes: A Diagnostic and Genetic Imaging Study in UK Biobank

BACKGROUND

Thyroid hormone is essential for optimal human neurodevelopment and may modify the risk of attention-deficit/hyperactivity disorder (ADHD). However, the brain structures involved are unknown and it is unclear if the adult brain is also susceptible to changes in thyroid status.

METHODS

We used International Classification of Disease-10 codes, polygenic thyroid scores at different thresholds of association with thyroid traits (PT-values), and image-derived phenotypes in UK Biobank (n = 18 825) to investigate the effects of a recorded diagnosis of thyroid disease and genetic risk for thyroid status on cerebellar and subcortical gray matter volume. Regional genetic pleiotropy between thyroid status and ADHD was explored using the GWAS-pairwise method.

RESULTS

A recorded diagnosis of hypothyroidism (n = 419) was associated with significant reductions in total cerebellar and pallidum gray matter volumes (β [95% CI] = -0.14[-0.23, -0.06], P = 0.0005 and β [95%CI] = -0.12 [-0.20, -0.04], P = 0.0042, respectively), mediated in part by increases in body mass index. While we found no evidence for total cerebellar volume alterations with increased polygenic scores for any thyroid trait, opposing influences of increased polygenic scores for hypo- and hyperthyroidism were found in the pallidum (PT < 1e-3: β [95% CI] = -0.02 [-0.03, -0.01], P = 0.0003 and PT < 1e-7: β [95% CI] = 0.02 [0.01, 0.03], P = 0.0003, respectively). Neither hypo- nor hyperthyroidism showed evidence of regional genetic pleiotropy with ADHD.

CONCLUSIONS

Thyroid status affects gray matter volume in adults, particularly at the level of the cerebellum and pallidum, with potential implications for the regulation of motor, cognitive, and affective function.

Overlap in genetic risk for cross-disorder vulnerability to mental disorders and genetic risk for altered subcortical brain volumes

BACKGROUND

There have been considerable recent advances in understanding the genetic architecture of psychiatric disorders as well as the underlying neurocircuitry. However, there is little work on the concordance of genetic variations that increase risk for cross-disorder vulnerability, and those that influence subcortical brain structures. We undertook a genome-wide investigation of the genetic overlap between cross-disorder vulnerability to psychiatric disorders (p-factor) and subcortical brain structures.

METHODS

Summary statistics were obtained from the PGC cross-disorder genome-wide association study (GWAS) (N_case= 232,964, N_control= 494,162) and the CHARGE-ENIGMA subcortical brain volumes GWAS (N=38,851). SNP effect concordance analysis (SECA) was used to assess pleiotropy and concordance. Linkage Disequilibrium (LD) Score Regression and ρ-HESS were used to assess genetic correlation and conditional false discovery (cFDR) was used to identify variants associated with p-factor, conditional on the variants association with subcortical brain volumes.

RESULTS

Evidence of global pleiotropy between p-factor and all subcortical brain regions was observed. Risk variants for p-factor correlated negatively with brainstem. A total of 787 LD-independent variants were significantly associated with p-factor when conditioned on the subcortical GWAS results. Gene set enrichment analysis of these variants implicated actin binding and neuronal regulation.

LIMITATIONS

SECA could be biased due to the potential presence of overlapping study participants in the p-factor and subcortical GWASs.

CONCLUSION

Findings of genome-wide pleiotropy and possible concordance between genetic variants that contribute to p-factor and smaller brainstem volumes, are consistent with previous work. cFDR results highlight actin binding and neuron regulation as key underlying mechanisms. Further fine-grained delineation of these mechanisms is needed to advance the field.

Quantitative endophenotypes as an alternative approach to understanding genetic risk in neurodegenerative diseases

Endophenotypes, as measurable intermediate features of human diseases, reflect underlying molecular mechanisms. The use of quantitative endophenotypes in genetic studies has improved our understanding of pathophysiological changes associated with diseases. The main advantage of the quantitative endophenotypes approach to study human diseases over a classic case-control study design is the inferred biological context that can enable the development of effective disease-modifying treatments. Here, we summarize recent progress on biomarkers for neurodegenerative diseases, including cerebrospinal fluid and blood-based, neuroimaging, neuropathological, and clinical studies. This review focuses on how endophenotypic studies have successfully linked genetic modifiers to disease risk, disease onset, or progression rate and provided biological context to genes identified in genome-wide association studies. Finally, we review critical methodological considerations for implementing this approach and future directions.

Structural Imaging in OCD

Anatomical imaging in OCD using magnetic resonance imaging (MRI) has been performed since the late 1980s. MRI research was further stimulated with the advent of automated image processing techniques such as voxel-based morphometry (VBM) and surface-based methods (e.g., FreeSurfer) which allow for detailed whole-brain data analyses. Early studies suggesting involvement of corticostriatal circuitry (particularly orbitofrontal cortex and ventral striatum) have been complemented by meta-analyses and pooled analyses indicating additional involvement of posterior brain regions, in particular parietal cortex. Recent large-scale meta-analyses from the ENIGMA consortium have revealed greater pallidum and smaller hippocampus volume in adult OCD, coupled with parietal cortical thinning. Frontal cortical thinning was only observed in medicated patients. Previous reports of symptom dimension-specific alterations were not confirmed. In paediatric OCD, thalamus enlargement has been a consistent finding. Studies investigating white matter volume (VBM) or integrity (using diffusion tensor imaging (DTI)) have shown mixed results, with recent DTI meta-analyses mainly showing involvement of posterior cortical-subcortical tracts in addition to subcortical-prefrontal connections. To which extent these abnormalities are unique to OCD or common to other psychiatric disorders is unclear, as few comparative studies have been performed. Overall, neuroanatomical alterations in OCD appear to be subtle and may vary with time, stressing the need for adequately powered longitudinal studies. Although multivariate approaches using machine learning methodologies have so far been disappointing in distinguishing individual OCD patients from healthy controls, including multimodal data in such analyses may aid in further establishing a neurobiological profile of OCD.

FreeSurfer ‐based segmentation of hippocampal subfields: A review of methods and applications, with a novel quality control procedure for ENIGMA studies and other collaborative efforts

Structural hippocampal abnormalities are common in many neurological and psychiatric disorders, and variation in hippocampal measures is related to cognitive performance and other complex phenotypes such as stress sensitivity. Hippocampal subregions are increasingly studied, as automated algorithms have become available for mapping and volume quantification. In the context of the Enhancing Neuro Imaging Genetics through Meta Analysis Consortium, several Disease Working Groups are using the FreeSurfer software to analyze hippocampal subregion (subfield) volumes in patients with neurological and psychiatric conditions along with data from matched controls. In this overview, we explain the algorithm's principles, summarize measurement reliability studies, and demonstrate two additional aspects (subfield autocorrelation and volume/reliability correlation) with illustrative data. We then explain the rationale for a standardized hippocampal subfield segmentation quality control (QC) procedure for improved pipeline harmonization. To guide researchers to make optimal use of the algorithm, we discuss how global size and age effects can be modeled, how QC steps can be incorporated and how subfields may be aggregated into composite volumes. This discussion is based on a synopsis of 162 published neuroimaging studies (01/2013–12/2019) that applied the FreeSurfer hippocampal subfield segmentation in a broad range of domains including cognition and healthy aging, brain development and neurodegeneration, affective disorders, psychosis, stress regulation, neurotoxicity, epilepsy, inflammatory disease, childhood adversity and posttraumatic stress disorder, and candidate and whole genome (epi‐)genetics. Finally, we highlight points where FreeSurfer‐based hippocampal subfield studies may be optimized.

Ten years of enhancing neuro-imaging genetics through meta-analysis: An overview from the ENIGMA Genetics Working Group

Here we review the motivation for creating the enhancing neuroimaging genetics through meta-analysis (ENIGMA) Consortium and the genetic analyses undertaken by the consortium so far. We discuss the methodological challenges, findings, and future directions of the genetics working group. A major goal of the working group is tackling the reproducibility crisis affecting "candidate gene" and genome-wide association analyses in neuroimaging. To address this, we developed harmonized analytic methods, and support their use in coordinated analyses across sites worldwide, which also makes it possible to understand heterogeneity in results across sites. These efforts have resulted in the identification of hundreds of common genomic loci robustly associated with brain structure. We have found both pleiotropic and specific genetic effects associated with brain structures, as well as genetic correlations with psychiatric and neurological diseases.

Biological and environmental predictors of heterogeneity in neurocognitive ageing: Evidence from Betula and other longitudinal studies

Individual differences in cognitive performance increase with advancing age, reflecting marked cognitive changes in some individuals along with little or no change in others. Genetic and lifestyle factors are assumed to influence cognitive performance in ageing by affecting the magnitude and extent of age-related brain changes (i.e., brain maintenance or atrophy), as well as the ability to recruit compensatory processes. The purpose of this review is to present findings from the Betula study and other longitudinal studies, with a focus on clarifying the role of key biological and environmental factors assumed to underlie individual differences in brain and cognitive ageing. We discuss the vital importance of sampling, analytic methods, consideration of non-ignorable dropout, and related issues for valid conclusions on factors that influence healthy neurocognitive ageing.

Transcriptomic Heterogeneity of Alzheimer's Disease Associated with Lipid Genetic Risk

Alzheimer's disease (AD) is a multifactorial disease that affects more than 5 million Americans. Multiple pathways might be involved in the AD pathogenesis. The implication of lipid genetic susceptibility on brain gene expression is yet to be investigated. The current study included 192 brain samples from AD patients who were enrolled in the ROSMAP study. The samples were genotyped and imputed to the HRC Reference Panel. Lipid polygenetic risk score was constructed from the weighted sum of genetic variants associated with low-density lipoprotein cholesterol (LDL-C). The gene expression was profiled by RNA sequencing, and the  association of gene expression with lipid polygenetic risk scores was tested by linear regression models adjusted for age, sex and APOE e4 alleles. Three genes were found to associate with lipid polygenetic risk scores, including HMCN2 (P = 3.6 × 10-7), PDLIM5 (P = 1.2 × 10-6), and FHL5 (P = 2.0 × 10-6). Network analysis revealed multiple related pathways, including dopaminergic synapse (P = 4.5 × 10-5), circadian entrainment (P = 1.1 × 10-4), and cholinergic synapse (P = 2.3 × 10-4). Our study underscores the importance of lipid regulation and metabolism to AD heterogeneity.

Multimodal hippocampal and amygdala subfield volumetry in polygenic risk for Alzheimer's disease

Preclinical models of Alzheimer's disease (AD) suggest that volumetric reductions in medial temporal lobe (MTL) structures manifest before clinical onset. AD polygenic risk scores (PRSs) are further linked to reduced MTL volumes (the hippocampus/amygdala); however, the relationship between the PRS and specific subregions remains unclear. We determine the relationship between the AD-PRSs and MTL subregions in a large sample of young participants (N = 730, aged 22–35 years) using a multimodal (T1w/T2w) approach. We first demonstrate that the PRSs for the hippocampus/amygdala predict their respective volumes and specific hippocampal subregions (pFDR < 0.05). We further observe negative relationships between the AD-PRSs and whole hippocampal/amygdala volumes. Critically, we demonstrate novel associations between the AD-PRSs and specific hippocampal subfields such as CA1 (β = −0.096, pFDR = 0.045) and the fissure (β = −0.101, pFDR = 0.041). We provide evidence that the AD-PRS is linked to specific MTL subfields decades before AD onset. This may help inform preclinical models of AD risk, providing additional specificity for intervention and further insight into mechanisms by which common AD variants confer susceptibility.

•

Polygenic risk for Alzheimer's disease (AD-PRS) explains significant proportion of AD.

•

AD-PRS also linked to hippocampus and amygdala volume.

•

AD-PRS is negatively associated with specific hippocampal subfields.

•

Polygenic AD models help us understand genetic contributions to medial temporal lobe nuclei.

Integrative analyses prioritize GNL3 as a risk gene for bipolar disorder

Genome-wide association studies (GWASs) have identified numerous single nucleotide polymorphisms (SNPs) associated with bipolar disorder (BD), but what the causal variants are and how they contribute to BD is largely unknown. In this study, we used FUMA, a GWAS annotation tool, to pinpoint potential causal variants and genes from the latest BD GWAS findings, and performed integrative analyses, including brain expression quantitative trait loci (eQTL), gene coexpression network, differential gene expression, protein-protein interaction, and brain intermediate phenotype association analysis to identify the functions of a prioritized gene and its connection to BD. Convergent lines of evidence prioritized protein-coding gene G Protein Nucleolar 3 (GNL3) as a BD risk gene, with integrative analyses revealing GNL3's roles in cell proliferation, neuronal functions, and brain phenotypes. We experimentally revealed that BD-related eQTL SNPs rs10865973, rs12635140, and rs4687644 regulate GNL3 expression using dual luciferase reporter assay and CRISPR interference experiment in human neural progenitor cells. We further identified that GNL3 knockdown and overexpression led to aberrant neuronal proliferation and differentiation, using two-dimensional human neural cell cultures and three-dimensional forebrain organoid model. This study gathers evidence that BD-related genetic variants regulate GNL3 expression which subsequently affects neuronal proliferation and differentiation.