Heritability of the shape of subcortical brain structures in the general population

Heritability of functional gradients in the human subcortico-cortical connectivity

The human subcortex plays a pivotal role in cognition and is widely implicated in the pathophysiology of many psychiatric disorders. However, the heritability of functional gradients based on subcortico-cortical functional connectivity remains elusive. Here, leveraging twin functional MRI (fMRI) data from both the Human Connectome Project (n = 1023) and the Adolescent Brain Cognitive Development study (n = 936) datasets, we construct large-scale subcortical functional gradients and delineate an increased principal functional gradient pattern from unimodal sensory/motor networks to transmodal association networks. We observed that this principal functional gradient is heritable, and the strength of heritability exhibits a heterogeneous pattern along a hierarchical unimodal-transmodal axis in subcortex for both young adults and children. Furthermore, employing a machine learning framework, we show that this heterogeneous pattern of the principal functional gradient in subcortex can accurately discern the relationship between monozygotic twin pairs and dizygotic twin pairs with an accuracy of 76.2% (P < 0.001). The heritability of functional gradients is associated with the anatomical myelin proxied by MRI-derived T1-weighted/T2-weighted (T1w/T2w) ratio mapping in subcortex. This study provides new insights into the biological basis of subcortical functional hierarchy by revealing the structural and genetic properties of the subcortical functional gradients.

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

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

Structural Neuroimaging of Hippocampus and Amygdala Subregions in Posttraumatic Stress Disorder: A Scoping Review

Numerous studies have explored the relationship between posttraumatic stress disorder (PTSD) and the hippocampus and the amygdala because both regions are implicated in the disorder's pathogenesis and pathophysiology. Nevertheless, those key limbic regions consist of functionally and cytoarchitecturally distinct substructures that may play different roles in the etiology of PTSD. Spurred by the availability of automatic segmentation software, structural neuroimaging studies of human hippocampal and amygdala subregions have proliferated in recent years. Here, we present a preregistered scoping review of the existing structural neuroimaging studies of the hippocampus and amygdala subregions in adults diagnosed with PTSD. A total of 3513 studies assessing subregion volumes were identified, 1689 of which were screened, and 21 studies were eligible for this review (total N = 2876 individuals). Most studies examined hippocampal subregions and reported decreased CA1, CA3, dentate gyrus, and subiculum volumes in PTSD. Fewer studies investigated amygdala subregions and reported altered lateral, basal, and central nuclei volumes in PTSD. This review further highlights the conceptual and methodological limitations of the current literature and identifies future directions to increase understanding of the distinct roles of hippocampal and amygdalar subregions in posttraumatic psychopathology.

Genomic loci influence patterns of structural covariance in the human brain

Normal and pathologic neurobiological processes influence brain morphology in coordinated ways that give rise to patterns of structural covariance (PSC) across brain regions and individuals during brain aging and diseases. The genetic underpinnings of these patterns remain largely unknown. We apply a stochastic multivariate factorization method to a diverse population of 50,699 individuals (12 studies and 130 sites) and derive data-driven, multi-scale PSCs of regional brain size. PSCs were significantly correlated with 915 genomic loci in the discovery set, 617 of which are newly identified, and 72% were independently replicated. Key pathways influencing PSCs involve reelin signaling, apoptosis, neurogenesis, and appendage development, while pathways of breast cancer indicate potential interplays between brain metastasis and PSCs associated with neurodegeneration and dementia. Using support vector machines, multi-scale PSCs effectively derive imaging signatures of several brain diseases. Our results elucidate genetic and biological underpinnings that influence structural covariance patterns in the human brain.

Surface-Based Neuroimaging Pattern of Multiple System Atrophy

RATIONALE AND OBJECTIVES

Overlapping parkinsonism, cerebellar ataxia, and pyramidal signs render challenges in the clinical diagnosis of multiple system atrophy (MSA). The neuroimaging pattern is valuable to understand its pathophysiology and improve its diagnostic effect.

MATERIALS AND METHODS

We retrospectively obtained magnetic resonance imaging and susceptibility-weighted imaging in patients with MSA (including parkinsonian type [MSA-P] and cerebellar type [MSA-C]), Parkinson's disease, and normal controls. We quantified neuroimaging features to identify the optimal threshold for diagnosis. Furthermore, we explore neuroimaging patterns of MSA by mapping the subcortical morphological alterations and constructing a diagnostic model.

RESULTS

Compared to controls, normalized putaminal volume significantly decreased in patients with MSA-P (P < .001) and normalized pontine volume significantly decreased in patients with MSA-C (P < .001). The Youden index of the threshold-based clinical prediction model was 0.871-0.928 in patients with MSA. The neuroimaging pattern in patients with MSA was jointly located in the lateral putamen, and the neuroimaging pattern prediction model achieved a classification accuracy of 83.9%-100%.

CONCLUSION

The quantitative neuroimaging features and surface-based morphologic anomalies represent the markers of MSA and open new avenues for personalized clinical diagnosis.

Structural and functional network analysis of twins using fMRI data

Similarities between twins have been widely demonstrated, underscoring the remarkable influence of genetics across numerous traits. In this study, we explore the genetic underpinnings of the human brain by examining MRI data from the Queensland Twin Imaging study. Specifically, this study seeks to compare brain structure and function between twins and unrelated subjects, with an emphasis on describing the effects of genetic factors. To achieve these goals, we employed the source-based morphometry method to extract intrinsic components and elucidate recognizable patterns. Our results show that twins exhibit a higher degree of similarity in gray and white matter density compared with unrelated individuals. In addition, four distinct states of brain activity were identified using coactivation patterns analysis. Furthermore, twins demonstrated a greater degree of similarity in the temporal and spatial features of each state compared with unrelated subjects. Taken together, these results support the hypothesis that twins show greater similarity in both brain structure and dynamic functional brain activity. Further exploration of these methods may advance our understanding of the complex interplay between genes, environment, and brain networks.

Heritability of cerebral blood flow in adolescent and young adult twins: an arterial spin labeling perfusion imaging study

Blood perfusion is a fundamental physiological property of all organs and is closely linked to brain metabolism. Genetic factors were reported to have important influences on cerebral blood flow. However, the profile of genetic contributions to cerebral blood flow in adolescents or young adults was underexplored. In this study, we recruited a sample of 65 pairs of same-sex adolescent or young adult twins undergoing resting arterial spin labeling imaging to conduct heritability analyses. Our findings indicate that genetic factors modestly affect cerebral blood flow in adolescents or young adults in the territories of left anterior cerebral artery and right posterior cerebral artery, with the primary contribution being to the frontal regions, cingulate gyrus, and striatum, suggesting a profile of genetic contributions to specific brain regions. Notably, the regions in the left hemisphere demonstrate the highest heritability in most regions examined. These results expand our knowledge of the genetic basis of cerebral blood flow in the developing brain and emphasize the importance of regional analysis in understanding the heritability of cerebral blood flow. Such insights may contribute to our understanding of the underlying genetic mechanism of brain functions and altered cerebral blood flow observed in youths with brain disorders.

Subcortical Brain Alterations in Carriers of Genomic Copy Number Variants

OBJECTIVE

Copy number variants (CNVs) are well-known genetic pleiotropic risk factors for multiple neurodevelopmental and psychiatric disorders (NPDs), including autism (ASD) and schizophrenia. Little is known about how different CNVs conferring risk for the same condition may affect subcortical brain structures and how these alterations relate to the level of disease risk conferred by CNVs. To fill this gap, the authors investigated gross volume, vertex-level thickness, and surface maps of subcortical structures in 11 CNVs and six NPDs.

METHODS

Subcortical structures were characterized using harmonized ENIGMA protocols in 675 CNV carriers (CNVs at 1q21.1, TAR, 13q12.12, 15q11.2, 16p11.2, 16p13.11, and 22q11.2; age range, 6-80 years; 340 males) and 782 control subjects (age range, 6-80 years; 387 males) as well as ENIGMA summary statistics for ASD, schizophrenia, attention deficit hyperactivity disorder, obsessive-compulsive disorder, bipolar disorder, and major depression.

RESULTS

All CNVs showed alterations in at least one subcortical measure. Each structure was affected by at least two CNVs, and the hippocampus and amygdala were affected by five. Shape analyses detected subregional alterations that were averaged out in volume analyses. A common latent dimension was identified, characterized by opposing effects on the hippocampus/amygdala and putamen/pallidum, across CNVs and across NPDs. Effect sizes of CNVs on subcortical volume, thickness, and local surface area were correlated with their previously reported effect sizes on cognition and risk for ASD and schizophrenia.

CONCLUSIONS

The findings demonstrate that subcortical alterations associated with CNVs show varying levels of similarities with those associated with neuropsychiatric conditions, as well distinct effects, with some CNVs clustering with adult-onset conditions and others with ASD. These findings provide insight into the long-standing questions of why CNVs at different genomic loci increase the risk for the same NPD and why a single CNV increases the risk for a diverse set of NPDs.

Heritability of cognitive abilities and regional brain structures in middle-aged to elderly East Asians

This study examined the single-nucleotide polymorphism heritability and genetic correlations of cognitive abilities and brain structural measures (regional subcortical volume and cortical thickness) in middle-aged and elderly East Asians (Korean) from the Gwangju Alzheimer's and Related Dementias cohort study. Significant heritability was found in memory function, caudate volume, thickness of the entorhinal cortices, pars opercularis, superior frontal gyri, and transverse temporal gyri. There were 3 significant genetic correlations between (i) the caudate volume and the thickness of the entorhinal cortices, (ii) the thickness of the superior frontal gyri and pars opercularis, and (iii) the thickness of the superior frontal and transverse temporal gyri. This is the first study to describe the heritability and genetic correlations of cognitive and neuroanatomical traits in middle-aged to elderly East Asians. Our results support the previous findings showing that genetic factors play a substantial role in the cognitive and neuroanatomical traits in middle to advanced age. Moreover, by demonstrating shared genetic effects on different brain regions, it gives us a genetic insight into understanding cognitive and brain changes with age, such as aging-related cognitive decline, cortical atrophy, and neural compensation.

Brain functional connectivity mirrors genetic pleiotropy in psychiatric conditions

Pleiotropy occurs when a genetic variant influences more than one trait. This is a key property of the genomic architecture of psychiatric disorders and has been observed for rare and common genomic variants. It is reasonable to hypothesize that the microscale genetic overlap (pleiotropy) across psychiatric conditions and cognitive traits may lead to similar overlaps at the macroscale brain level such as large-scale brain functional networks. We took advantage of brain connectivity, measured by resting-state functional MRI to measure the effects of pleiotropy on large-scale brain networks, a putative step from genes to behaviour. We processed nine resting-state functional MRI datasets including 32 726 individuals and computed connectome-wide profiles of seven neuropsychiatric copy-number-variants, five polygenic scores, neuroticism and fluid intelligence as well as four idiopathic psychiatric conditions. Nine out of 19 pairs of conditions and traits showed significant functional connectivity correlations (rFunctional connectivity), which could be explained by previously published levels of genomic (rGenetic) and transcriptomic (rTranscriptomic) correlations with moderate to high concordance: rGenetic-rFunctional connectivity = 0.71 [0.40-0.87] and rTranscriptomic-rFunctional connectivity = 0.83 [0.52; 0.94]. Extending this analysis to functional connectivity profiles associated with rare and common genetic risk showed that 30 out of 136 pairs of connectivity profiles were correlated above chance. These similarities between genetic risks and psychiatric disorders at the connectivity level were mainly driven by the overconnectivity of the thalamus and the somatomotor networks. Our findings suggest a substantial genetic component for shared connectivity profiles across conditions and traits, opening avenues to delineate general mechanisms-amenable to intervention-across psychiatric conditions and genetic risks.

Subcortical brain alterations in carriers of genomic copy number variants

Objectives

Copy number variants (CNVs) are well-known genetic pleiotropic risk factors for multiple neurodevelopmental and psychiatric disorders (NPDs) including autism (ASD) and schizophrenia (SZ). Overall, little is known about how different CNVs conferring risk for the same condition may affect subcortical brain structures and how these alterations relate to the level of disease risk conferred by CNVs. To fill this gap, we investigated gross volume, and vertex level thickness and surface maps of subcortical structures in 11 different CNVs and 6 different NPDs.

Methods

Subcortical structures were characterized using harmonized ENIGMA protocols in 675 CNV carriers (at the following loci: 1q21.1, TAR, 13q12.12, 15q11.2, 16p11.2, 16p13.11, and 22q11.2) and 782 controls (Male/Female: 727/730; age-range: 6-80 years) as well as ENIGMA summary-statistics for ASD, SZ, ADHD, Obsessive-Compulsive-Disorder, Bipolar-Disorder, and Major-Depression.

Results

Nine of the 11 CNVs affected volume of at least one subcortical structure. The hippocampus and amygdala were affected by five CNVs. Effect sizes of CNVs on subcortical volume, thickness and local surface area were correlated with their previously reported effect sizes on cognition and risk for ASD and SZ. Shape analyses were able to identify subregional alterations that were averaged out in volume analyses. We identified a common latent dimension - characterized by opposing effects on basal ganglia and limbic structures - across CNVs and across NPDs.

Conclusion

Our findings demonstrate that subcortical alterations associated with CNVs show varying levels of similarities with those associated with neuropsychiatric conditions. We also observed distinct effects with some CNVs clustering with adult conditions while others clustered with ASD. This large cross-CNV and NPDs analysis provide insight into the long-standing questions of why CNVs at different genomic loci increase the risk for the same NPD, as well as why a single CNV increases the risk for a diverse set of NPDs.

The subcortical correlates of autistic traits in school-age children: a population-based neuroimaging study

BACKGROUND

There is emerging evidence that the neuroanatomy of autism forms a spectrum which extends into the general population. However, whilst several studies have identified cortical morphology correlates of autistic traits, it is not established whether morphological differences are present in the subcortical structures of the brain. Additionally, it is not clear to what extent previously reported structural associations may be confounded by co-occurring psychopathology. To address these questions, we utilised neuroimaging data from the Adolescent Brain Cognitive Development Study to assess whether a measure of autistic traits was associated with differences in child subcortical morphology, and if any observed differences persisted after adjustment for child internalising and externalising symptoms.

METHODS

Our analyses included data from 7005 children aged 9-10 years (female: 47.19%) participating in the Adolescent Brain Cognitive Development Study. Autistic traits were assessed using scores from the Social Responsiveness Scale (SRS). Volumes of subcortical regions of interest were derived from structural magnetic resonance imaging data.

RESULTS

Overall, we did not find strong evidence for an association of autistic traits with differences in subcortical morphology in this sample of school-aged children. Whilst lower absolute volumes of the nucleus accumbens and putamen were associated with higher scores of autistic traits, these differences did not persist once a global measure of brain size was accounted for.

LIMITATIONS

It is important to note that autistic traits were assessed using the SRS, of which higher scores are associated with general behavioural problems, and therefore may not be wholly indicative of autism-specific symptoms. In addition, individuals with a moderate or severe autism diagnosis were excluded from the ABCD study, and thus, the average level of autistic traits will be lower than in the general population which may bias findings towards the null.

CONCLUSIONS

These findings from our well-powered study suggest that other metrics of brain morphology, such as cortical morphology or shape-based phenotypes, may be stronger candidates to prioritise when attempting to identify robust neuromarkers of autistic traits.

Subcortical brain structures and the risk of dementia in the Rotterdam Study

INTRODUCTION

Volumetric and morphological changes in subcortical brain structures are present in persons with dementia, but it is unknown if these changes occur prior to diagnosis.

METHODS

Between 2005 and 2016, 5522 Rotterdam Study participants (mean age: 64.4) underwent cerebral magnetic resonance imaging (MRI) and were followed for development of dementia until 2018. Volume and shape measures were obtained for seven subcortical structures.

RESULTS

During 12 years of follow-up, 272 dementia cases occurred. Mean volumes of thalamus (hazard ratio [HR] per standard deviation [SD] decrease 1.94, 95% confidence interval [CI]: 1.55-2.43), amygdala (HR 1.66, 95% CI: 1.44-1.92), and hippocampus (HR 1.64, 95% CI: 1.43-1.88) were strongly associated with dementia risk. Associations for accumbens, pallidum, and caudate volumes were less pronounced. Shape analyses identified regional surface changes in the amygdala, limbic thalamus, and caudate.

DISCUSSION

Structure of the amygdala, thalamus, hippocampus, and caudate is associated with risk of dementia in a large population-based cohort of older adults.

Heritability of Subcortical Grey Matter Structures

Background and Objectives: Subcortical grey matter structures play essential roles in cognitive, affective, social, and motoric functions in humans. Their volume changes with age, and decreased volumes have been linked with many neuropsychiatric disorders. The aim of our study was to examine the heritability of six subcortical brain volumes (the amygdala, caudate nucleus, pallidum, putamen, thalamus, and nucleus accumbens) and four general brain volumes (the total intra-cranial volume and the grey matter, white matter, and cerebrospinal fluid (CSF) volume) in twins. Materials and Methods: A total of 118 healthy adult twins from the Hungarian Twin Registry (86 monozygotic and 32 dizygotic; median age 50 ± 27 years) underwent brain magnetic resonance imaging. Two automated volumetry pipelines, Computational Anatomy Toolbox 12 (CAT12) and volBrain, were used to calculate the subcortical and general brain volumes from three-dimensional T1-weighted images. Age- and sex-adjusted monozygotic and dizygotic intra-pair correlations were calculated, and the univariate ACE model was applied. Pearson's correlation test was used to compare the results obtained by the two pipelines. Results: The age- and sex-adjusted heritability estimates, using CAT12 for the amygdala, caudate nucleus, pallidum, putamen, and nucleus accumbens, were between 0.75 and 0.95. The thalamus volume was more strongly influenced by common environmental factors (C = 0.45-0.73). The heritability estimates, using volBrain, were between 0.69 and 0.92 for the nucleus accumbens, pallidum, putamen, right amygdala, and caudate nucleus. The left amygdala and thalamus were more strongly influenced by common environmental factors (C = 0.72-0.85). A strong correlation between CAT12 and volBrain (r = 0.74-0.94) was obtained for all volumes. Conclusions: The majority of examined subcortical volumes appeared to be strongly heritable. The thalamus was more strongly influenced by common environmental factors when investigated with both segmentation methods. Our results underline the importance of identifying the relevant genes responsible for variations in the subcortical structure volume and associated diseases.

Alteration of brain nuclei in obese children with and without Prader-Willi syndrome

Introduction: Prader-Willi syndrome (PWS) is a multisystem genetic imprinting disorder mainly characterized by hyperphagia and childhood obesity. Extensive structural alterations are expected in PWS patients, and their influence on brain nuclei should be early and profound. To date, few studies have investigated brain nuclei in children with PWS, although functional and structural alterations of the cortex have been reported widely.

Methods: In the current study, we used T1-weighted magnetic resonance imaging to investigate alterations in brain nuclei by three automated analysis methods: shape analysis to evaluate the shape of 14 cerebral nuclei (bilateral thalamus, caudate, putamen, globus pallidus, hippocampus, amygdala, and nucleus accumbens), automated segmentation methods integrated in Freesurfer 7.2.0 to investigate the volume of hypothalamic subregions, and region of interest-based analysis to investigate the volume of deep cerebellar nuclei (DCN). Twelve age- and sex-matched children with PWS, 18 obese children without PWS (OB) and 18 healthy controls participated in this study.

Results: Compared with control and OB individuals, the PWS group exhibited significant atrophy in the bilateral thalamus, pallidum, hippocampus, amygdala, nucleus accumbens, right caudate, bilateral hypothalamus (left anterior-inferior, bilateral posterior, and bilateral tubular inferior subunits) and bilateral DCN (dentate, interposed, and fastigial nuclei), whereas no significant difference was found between the OB and control groups.

Discussion: Based on our evidence, we suggested that alterations in brain nuclei influenced by imprinted genes were associated with clinical manifestations of PWS, such as eating disorders, cognitive disability and endocrine abnormalities, which were distinct from the neural mechanisms of obese children.

Shape analysis of subcortical structures in obsessive-compulsive disorder and the relationship with comorbid anxiety, depression, and medication use: A meta-analysis by the OCD Brain Imaging Consortium

OBJECTIVE

Neuroimaging studies of obsessive-compulsive disorder (OCD) patients have highlighted the important role of deep gray matter structures. Less work has however focused on subcortical shape in OCD patients.

METHODS

Here we pooled brain MRI scans from 412 OCD patients and 368 controls to perform a meta-analysis utilizing the ENIGMA-Shape protocol. In addition, we investigated modulating effects of medication status, comorbid anxiety or depression, and disease duration on subcortical shape.

RESULTS

There was no significant difference in shape thickness or surface area between OCD patients and healthy controls. For the subgroup analyses, OCD patients with comorbid depression or anxiety had lower thickness of the hippocampus and caudate nucleus and higher thickness of the putamen and pallidum compared to controls. OCD patients with comorbid depression had lower shape surface area in the thalamus, caudate nucleus, putamen, hippocampus, and nucleus accumbens and higher shape surface area in the pallidum. OCD patients with comorbid anxiety had lower shape surface area in the putamen and the left caudate nucleus and higher shape surface area in the pallidum and the right caudate nucleus. Further, OCD patients on medication had lower shape thickness of the putamen, thalamus, and hippocampus and higher thickness of the pallidum and caudate nucleus, as well as lower shape surface area in the hippocampus and amygdala and higher surface area in the putamen, pallidum, and caudate nucleus compared to controls. There were no significant differences between OCD patients without co-morbid anxiety and/or depression and healthy controls on shape measures. In addition, there were also no significant differences between OCD patients not using medication and healthy controls.

CONCLUSIONS

The findings here are partly consistent with prior work on brain volumes in OCD, insofar as they emphasize that alterations in subcortical brain morphology are associated with comorbidity and medication status. Further work is needed to understand the biological processes contributing to subcortical shape.

A review of brain imaging biomarker genomics in Alzheimer’s disease: implementation and perspectives

Alzheimer’s disease (AD) is a progressive neurodegenerative disease with phenotypic changes closely associated with both genetic variants and imaging pathology. Brain imaging biomarker genomics has been developed in recent years to reveal potential AD pathological mechanisms and provide early diagnoses. This technique integrates multimodal imaging phenotypes with genetic data in a noninvasive and high-throughput manner. In this review, we summarize the basic analytical framework of brain imaging biomarker genomics and elucidate two main implementation scenarios of this technique in AD studies: (1) exploring novel biomarkers and seeking mutual interpretability and (2) providing a diagnosis and prognosis for AD with combined use of machine learning methods and brain imaging biomarker genomics. Importantly, we highlight the necessity of brain imaging biomarker genomics, discuss the strengths and limitations of current methods, and propose directions for development of this research field.

Polygenic risk scores for antisocial behavior in relation to amygdala morphology across an attention deficit hyperactivity disorder case-control sample with and without disruptive behavior

Antisocial and aggressive behaviors show considerable heritability and are central to disruptive behavior disorders (DBDs), but are also frequently observed in attention deficit hyperactivity disorder (ADHD). While the amygdala is implicated as a key neural structure, it remains unclear whether common genetic variants underlie this brain-behavior association. We hypothesized that polygenic (risk) scores for antisocial and aggressive behaviors (ASB-PRS) would be related to amygdala morphology. Using the Broad Antisocial Behavior Consortium genome-wide association study (GWAS; mostly population based cohorts), we calculated ASB-PRS in the NeuroIMAGE I ADHD case-control sample with varying levels of DBD symptomatology (n=679 from 379 families, aged 7 - 29). We first investigated associations of several ASB-PRS p value thresholds with the presence of DBD symptoms and self-reported antisocial behavior (ASB) to determine the threshold for further analyses. This PRS was then related to amygdala volume and shape using regression and vertex-wise analyses. Our results showed associations of ASB-PRS with the presence of DBD symptoms, self-reported ASB, and left basolateral amygdala shape, independent of ADHD symptom severity and ADHD-PRS, with a relative outward displacement of the vertices. No associations of ASB-PRS, DBD symptoms or self-reported ASB with amygdala volume were found. Our results indicate that genetic risk for antisocial and aggressive behaviors is related to amygdala shape alterations, and point to genetic sharing across different DBD and ASB and aggression-related phenotypes as a spectrum of genetically related quantitative traits. Additionally, our findings support the utility of vertex-based shape analyses in genetic studies of ASB, aggression, and DBDs.

Parsimonious model for mass-univariate vertexwise analysis

Purpose: Covariance between gray-matter measurements can reflect structural or functional brain networks though it has also been shown to be influenced by confounding factors (e.g., age, head size, and scanner), which could lead to lower mapping precision (increased size of associated clusters) and create distal false positives associations in mass-univariate vertexwise analyses. Approach: We evaluated this concern by performing state-of-the-art mass-univariate analyses (general linear model, GLM) on traits simulated from real vertex-wise gray matter data (including cortical and subcortical thickness and surface area). We contrasted the results with those from linear mixed models (LMMs), which have been shown to overcome similar issues in omics association studies. Results: We showed that when performed on a large sample ( N = 8662 , UK Biobank), GLMs yielded greatly inflated false positive rate (cluster false discovery rate > 0.6 ). We showed that LMMs resulted in more parsimonious results: smaller clusters and reduced false positive rate but at a cost of increased computation. Next, we performed mass-univariate association analyses on five real UKB traits (age, sex, BMI, fluid intelligence, and smoking status) and LMM yielded fewer and more localized associations. We identified 19 significant clusters displaying small associations with age, sex, and BMI, which suggest a complex architecture of at least dozens of associated areas with those phenotypes. Conclusions: The published literature could contain a large proportion of redundant (possibly confounded) associations that are largely prevented using LMMs. The parsimony of LMMs results from controlling for the joint effect of all vertices, which prevents local and distal redundant associations from reaching significance.

Pallidal Structural Changes Related to Levodopa-induced Dyskinesia in Parkinson's Disease

Background

Despite the clinical impact of levodopa-induced dyskinesia (LID) in Parkinson's disease (PD), the mechanism, especially the role of basal ganglia (BG), is not fully elucidated yet. We investigated the BG structural changes related to LID in PD using a surface-based shape analysis technique.

Methods

We recruited patients with PD who developed LID within 3 years (LID group, 28 patients) and who did not develop it after 7 years (non-LID group, 35 patients) from levodopa treatment for the extreme case-control study. BG structure volumes were measured using volumetry analysis and the surface-based morphometry feature (i.e., Jacobian) from the subcortical surface vertices. We compared the volume and Jacobian of meshes in the regions between the two groups. We also performed a correlation analysis between local atrophy and the severity of LID. Additionally, we evaluated structural connectivity profiles from globus pallidus interna and externa (GPi and GPe) to other brain structures based on the group comparison.

Results

The demographic and clinical data showed no significant difference except for disease duration, treatment duration, parkinsonism severity, and levodopa equivalent dose. The LID group had more local atrophies of vertices in the right GPi than the non-LID group, despite no difference in volumes. Furthermore, the LID group demonstrated significantly reduced structural connectivity between left GPi and thalamus.

Conclusion

This is the first demonstration of distinct shape alterations of basal ganglia structures, especially GPi, related to LID in PD. Considering both direct and indirect BG pathways share the connection between GPi and thalamus, the BG pathway plays a crucial role in the development of LID.

Shapes of subcortical structures in adolescents with and without familial history of substance use disorder

It has been suggested that intergenerational transmission of risk for substance use disorder (SUD) manifests in the brain anatomy of substance naïve adolescents. While volume and shapes of subcortical structures (SSS) have been shown to be heritable, these structures, especially the pallidum, putamen, nucleus accumbens, and hippocampus, have also been associated with substance use disorders. However, it is not clear if those anatomical differences precede substance use or are the result of that use. Therefore, we examined if volume and SSS of adolescents with a family history (FH+) of SUD differed from adolescents without such a history (FH−). Because risk for SUD is associated with anxiety and impulsivity, we also examined correlations between these psychological characteristics and volume/SSS. Using structural MRI and FSL software, we segmented subcortical structures and obtained indices of SSS and volumes of 64 FH+ and 58 FH− adolescents. We examined group differences in volume and SSS, and the correlations between volume/SSS and trait anxiety and impulsivity. FH+ adolescents had a significant inward deformation in the shape of the right anterior hippocampus compared to FH− adolescents, while the volume of this structure did not differ between groups. Neither shape nor volume of the other subcortical structures differed between groups. In the FH+ adolescents, the left hippocampus shape was positively correlated with both trait anxiety and impulsivity, while in FH− adolescents a negative correlation pattern of SSS was seen in the hippocampus. SSS appears to capture local anatomical features that traditional volumetric analysis does not. The inward shape deformation in the right anterior hippocampus in FH+ adolescents may be related to the known increased risk for behavioral dysregulation leading to SUD in FH+ offspring. Hippocampus shape also exhibits opposite patterns of correlation with anxiety and impulsivity scores across the FH+ and FH− adolescents. These novel findings may reveal neural correlates, not captured by traditional volumetric analysis, of familial transmission of increased vulnerability to SUD.

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.

Charting the human amygdala development across childhood and adolescence: Manual and automatic segmentation

The developmental pattern of the amygdala throughout childhood and adolescence has been inconsistently reported in previous neuroimaging studies. Given the relatively small size of the amygdala on full brain MRI scans, discrepancies may be partly due to methodological differences in amygdalar segmentation. To investigate the impact of volume extraction methods on amygdala volume, we compared FreeSurfer, FSL and volBrain segmentation measurements with those obtained by manual tracing. The manual tracing method, which we used as the 'gold standard', exhibited almost perfect intra- and inter-rater reliability. We observed systematic differences in amygdala volumes between automatic (FreeSurfer and volBrain) and manual methods. Specifically, compared with the manual tracing, FreeSurfer estimated larger amygdalae, and volBrain produced smaller amygdalae while FSL demonstrated a mixed pattern. The tracing bias was not uniform, but higher for smaller amygdalae. We further modeled amygdalar growth curves using accelerated longitudinal cohort data from the Chinese Color Nest Project (http://deepneuro.bnu.edu.cn/?p=163). Trajectory modeling and statistical assessments of the manually traced amygdalae revealed linearly increasing and parallel developmental patterns for both girls and boys, although the amygdalae of boys were larger than those of girls. Compared to these trajectories, the shapes of developmental curves were similar when using the volBrain derived volumes. FreeSurfer derived trajectories had more nonlinearities and appeared flatter. FSL derived trajectories demonstrated an inverted U shape and were significantly different from those derived from manual tracing method. The use of amygdala volumes adjusted for total gray-matter volumes, but not intracranial volumes, resolved the shape discrepancies and led to reproducible growth curves between manual tracing and the automatic methods (except FSL). Our findings revealed steady growth of the human amygdala, mirroring its functional development across the school age. Methodological improvements are warranted for current automatic tools to achieve more accurate amygdala structure at school age, calling for next generation tools.

Whole-brain morphological alterations associated with trigeminal neuralgia

Background

Novel neuroimaging strategies have the potential to offer new insights into the mechanistic basis for trigeminal neuralgia (TN). The present study aims to conduct whole-brain morphometry analyses of TN patients and to assess the value of group-level neocortical and subcortical structural patterns as tools for diagnostic biomarker exploration.

Methods

Cortical thickness, surface area, and myelin levels in the neocortex were measured via magnetic resonance imaging (MRI). The radial distance and the Jacobian determinant of the subcortex in 43 TN patients and 43 matched controls were compared. Pattern learning algorithms were employed to establish the utility of group-level MRI findings as tools for predicting TN. An additional 40 control patients with hemifacial spasms were then evaluated to assess algorithm sensitivity and specificity.

Results

TN patients exhibited reductions in cortical indices in the anterior cingulate cortex (ACC), the midcingulate cortex (MCC), and the posterior cingulate cortex (PCC) relative to controls. They further presented with widespread subcortical volume reduction that was most evident in the putamen, the thalamus, the accumbens, the pallidum, and the hippocampus. Whole brain-level morphological alterations successfully enable automated TN diagnosis with high specificity (TN: 95.35 %; disease controls: 46.51 %).

Conclusions

TN is associated with a distinctive whole-brain structural neuroimaging pattern, underscoring the value of machine learning as an approach to differentiating between morphological phenotypes, ultimately revealing the full spectrum of this disease and highlighting relevant diagnostic biomarkers.

Supplementary Information

The online version contains supplementary material available at 10.1186/s10194-021-01308-5.

Inferring the heritability of large-scale functional networks with a multivariate ACE modeling approach

Recent evidence suggests that the human functional connectome is stable at different timescales and is unique. These characteristics posit the functional connectome not only as an individual marker but also as a powerful discriminatory measure characterized by high intersubject variability. Among distinct sources of intersubject variability, the long-term sources include functional patterns that emerge from genetic factors. Here, we sought to investigate the contribution of additive genetic factors to the variability of functional networks by determining the heritability of the connectivity strength in a multivariate fashion. First, we reproduced and extended the connectome fingerprinting analysis to the identification of twin pairs. Then, we estimated the heritability of functional networks by a multivariate ACE modeling approach with bootstrapping. Twin pairs were identified above chance level using connectome fingerprinting, with monozygotic twin identification accuracy equal to 57.2% on average for whole-brain connectome. Additionally, we found that a visual (0.37), the medial frontal (0.31), and the motor (0.30) functional networks were the most influenced by additive genetic factors. Our findings suggest that genetic factors not only partially determine intersubject variability of the functional connectome, such that twins can be identified using connectome fingerprinting, but also differentially influence connectivity strength in large-scale functional networks.

The functional connectome is a unique representation of the functional organization of the human brain. As such, it has been extensively used as an individual marker, a “fingerprint,” because of its high intersubject variability. Here, we sought to investigate the influence of genetic factors on intersubject variability of functional networks. Therefore, we extended the connectome fingerprinting analysis to the identification of twin pairs, and we estimated the heritability of functional networks by a multivariate ACE modeling approach with bootstrapping. We found that genetic factors not only partially determine intersubject variability of the functional connectome, such that monozygotic twin identification accuracy achieved 57.2% on average using whole-brain connectome in the fingerprinting analysis, but also differentially influence connectivity strength in large-scale functional networks.

Genetic Influences on Hippocampal Subfields: An Emerging Area of Neuroscience Research

There is clear evidence that hippocampal subfield volumes have partly distinct genetic determinants associated with specific biological processes. The identification of genetic correlates of hippocampal subfield volumes may help to elucidate the mechanisms of neurologic diseases, as well as aging and neurodegenerative processes. However, despite the emerging interest in this area of research, the current knowledge of the genetic architecture of hippocampal subfields has not yet been consolidated. We aimed to provide a review of the current evidence from genetic studies of hippocampal subfields, highlighting current priorities and upcoming challenges. The limited number of studies investigating the influential genetic effects on hippocampal subfields, a lack of replicated results and longitudinal designs, and modest sample sizes combined with insufficient standardization of protocols are identified as the most pressing challenges in this emerging area of research.

A unified framework for association and prediction from vertex-wise grey-matter structure

The recent availability of large‐scale neuroimaging cohorts facilitates deeper characterisation of the relationship between phenotypic and brain architecture variation in humans. Here, we investigate the association (previously coined morphometricity) of a phenotype with all 652,283 vertex‐wise measures of cortical and subcortical morphology in a large data set from the UK Biobank (UKB; N = 9,497 for discovery, N = 4,323 for replication) and the Human Connectome Project (N = 1,110). We used a linear mixed model with the brain measures of individuals fitted as random effects with covariance relationships estimated from the imaging data. We tested 167 behavioural, cognitive, psychiatric or lifestyle phenotypes and found significant morphometricity for 58 phenotypes (spanning substance use, blood assay results, education or income level, diet, depression, and cognition domains), 23 of which replicated in the UKB replication set or the HCP. We then extended the model for a bivariate analysis to estimate grey‐matter correlation between phenotypes, which revealed that body size (i.e., height, weight, BMI, waist and hip circumference, body fat percentage) could account for a substantial proportion of the morphometricity (confirmed using a conditional analysis), providing possible insight into previous MRI case–control results for psychiatric disorders where case status is associated with body mass index. Our LMM framework also allowed to predict some of the associated phenotypes from the vertex‐wise measures, in two independent samples. Finally, we demonstrated additional new applications of our approach (a) region of interest (ROI) analysis that retain the vertex‐wise complexity; (b) comparison of the information retained by different MRI processings.

Heritability of Regional Brain Volumes in Large-Scale Neuroimaging and Genetic Studies

Brain genetics is an active research area. The degree to which genetic variants impact variations in brain structure and function remains largely unknown. We examined the heritability of regional brain volumes (P ~ 100) captured by single-nucleotide polymorphisms (SNPs) in UK Biobank (n ~ 9000). We found that regional brain volumes are highly heritable in this study population and common genetic variants can explain up to 80% of their variabilities (median heritability 34.8%). We observed omnigenic impact across the genome and examined the enrichment of SNPs in active chromatin regions. Principal components derived from regional volume data are also highly heritable, but the amount of variance in brain volume explained by the component did not seem to be related to its heritability. Heritability estimates vary substantially across large-scale functional networks, exhibit a symmetric pattern across left and right hemispheres, and are consistent in females and males (correlation = 0.638). We repeated the main analysis in Alzheimer's Disease Neuroimaging Initiative (n ~ 1100), Philadelphia Neurodevelopmental Cohort (n ~ 600), and Pediatric Imaging, Neurocognition, and Genetics (n ~ 500) datasets, which demonstrated that more stable estimates can be obtained from the UK Biobank.

Bayesian sparse heritability analysis with high-dimensional neuroimaging phenotypes

Heritability analysis plays a central role in quantitative genetics to describe genetic contribution to human complex traits and prioritize downstream analyses under large-scale phenotypes. Existing works largely focus on modeling single phenotype and currently available multivariate phenotypic methods often suffer from scaling and interpretation. In this article, motivated by understanding how genetic underpinning impacts human brain variation, we develop an integrative Bayesian heritability analysis to jointly estimate heritabilities for high-dimensional neuroimaging traits. To induce sparsity and incorporate brain anatomical configuration, we impose hierarchical selection among both regional and local measurements based on brain structural network and voxel dependence. We also use a nonparametric Dirichlet process mixture model to realize grouping among single nucleotide polymorphism-associated phenotypic variations, providing biological plausibility. Through extensive simulations, we show the proposed method outperforms existing ones in heritability estimation and heritable traits selection under various scenarios. We finally apply the method to two large-scale imaging genetics datasets: the Alzheimer's Disease Neuroimaging Initiative and United Kingdom Biobank and show biologically meaningful results.

The reliability and heritability of cortical folds and their genetic correlations across hemispheres

Cortical folds help drive the parcellation of the human cortex into functionally specific regions. Variations in the length, depth, width, and surface area of these sulcal landmarks have been associated with disease, and may be genetically mediated. Before estimating the heritability of sulcal variation, the extent to which these metrics can be reliably extracted from in-vivo MRI must be established. Using four independent test-retest datasets, we found high reliability across the brain (intraclass correlation interquartile range: 0.65-0.85). Heritability estimates were derived for three family-based cohorts using variance components analysis and pooled (total N > 3000); the overall sulcal heritability pattern was correlated to that derived for a large population cohort (N > 9000) calculated using genomic complex trait analysis. Overall, sulcal width was the most heritable metric, and earlier forming sulci showed higher heritability. The inter-hemispheric genetic correlations were high, yet select sulci showed incomplete pleiotropy, suggesting hemisphere-specific genetic influences.

Brain imaging genomics: influences of genomic variability on the structure and function of the human brain

Brain imaging genomics is an emerging discipline in which genomic and brain imaging data are integrated in order to elucidate the molecular mechanisms that underly brain phenotypes and diseases, including neuropsychiatric disorders. As with all genetic analyses of complex traits and diseases, brain imaging genomics has evolved from small, individual candidate gene investigations towards large, collaborative genome-wide association studies. Recent investigations, mostly population-based, have studied well-powered cohorts comprising tens of thousands of individuals and identified multiple robust associations of single-nucleotide polymorphisms and copy number variants with structural and functional brain phenotypes. Such systematic genomic screens of millions of genetic variants have generated initial insights into the genetic architecture of brain phenotypes and demonstrated that their etiology is polygenic in nature, involving multiple common variants with small effect sizes and rare variants with larger effect sizes. Ongoing international collaborative initiatives are now working to obtain a more complete picture of the underlying biology. As in other complex phenotypes, novel approaches – such as gene–gene interaction, gene–environment interaction, and epigenetic analyses – are being implemented in order to investigate their contribution to the observed phenotypic variability. An important consideration for future research will be the translation of brain imaging genomics findings into clinical practice.

Mapping Subcortical Brain Alterations in 22q11.2 Deletion Syndrome: Effects of Deletion Size and Convergence With Idiopathic Neuropsychiatric Illness

OBJECTIVE

22q11.2 deletion syndrome (22q11DS) is among the strongest known genetic risk factors for schizophrenia. Previous studies have reported variable alterations in subcortical brain structures in 22q11DS. To better characterize subcortical alterations in 22q11DS, including modulating effects of clinical and genetic heterogeneity, the authors studied a large multicenter neuroimaging cohort from the ENIGMA 22q11.2 Deletion Syndrome Working Group.

METHODS

Subcortical structures were measured using harmonized protocols for gross volume and subcortical shape morphometry in 533 individuals with 22q11DS and 330 matched healthy control subjects (age range, 6-56 years; 49% female).

RESULTS

Compared with the control group, the 22q11DS group showed lower intracranial volume (ICV) and thalamus, putamen, hippocampus, and amygdala volumes and greater lateral ventricle, caudate, and accumbens volumes (Cohen's d values, -0.90 to 0.93). Shape analysis revealed complex differences in the 22q11DS group across all structures. The larger A-D deletion was associated with more extensive shape alterations compared with the smaller A-B deletion. Participants with 22q11DS with psychosis showed lower ICV and hippocampus, amygdala, and thalamus volumes (Cohen's d values, -0.91 to 0.53) compared with participants with 22q11DS without psychosis. Shape analysis revealed lower thickness and surface area across subregions of these structures. Compared with subcortical findings from other neuropsychiatric disorders studied by the ENIGMA consortium, significant convergence was observed between participants with 22q11DS with psychosis and participants with schizophrenia, bipolar disorder, major depressive disorder, and obsessive-compulsive disorder.

CONCLUSIONS

In the largest neuroimaging study of 22q11DS to date, the authors found widespread alterations to subcortical brain structures, which were affected by deletion size and psychotic illness. Findings indicate significant overlap between 22q11DS-associated psychosis, idiopathic schizophrenia, and other severe neuropsychiatric illnesses.

Subcortical shape alterations in major depressive disorder: Findings from the ENIGMA major depressive disorder working group

Alterations in regional subcortical brain volumes have been investigated as part of the efforts of an international consortium, ENIGMA, to identify reliable neural correlates of major depressive disorder (MDD). Given that subcortical structures are comprised of distinct subfields, we sought to build significantly from prior work by precisely mapping localized MDD‐related differences in subcortical regions using shape analysis. In this meta‐analysis of subcortical shape from the ENIGMA‐MDD working group, we compared 1,781 patients with MDD and 2,953 healthy controls (CTL) on individual measures of shape metrics (thickness and surface area) on the surface of seven bilateral subcortical structures: nucleus accumbens, amygdala, caudate, hippocampus, pallidum, putamen, and thalamus. Harmonized data processing and statistical analyses were conducted locally at each site, and findings were aggregated by meta‐analysis. Relative to CTL, patients with adolescent‐onset MDD (≤ 21 years) had lower thickness and surface area of the subiculum, cornu ammonis (CA) 1 of the hippocampus and basolateral amygdala (Cohen's d = −0.164 to −0.180). Relative to first‐episode MDD, recurrent MDD patients had lower thickness and surface area in the CA1 of the hippocampus and the basolateral amygdala (Cohen's d = −0.173 to −0.184). Our results suggest that previously reported MDD‐associated volumetric differences may be localized to specific subfields of these structures that have been shown to be sensitive to the effects of stress, with important implications for mapping treatments to patients based on specific neural targets and key clinical features.

ENIGMA and global neuroscience: A decade of large-scale studies of the brain in health and disease across more than 40 countries

This review summarizes the last decade of work by the ENIGMA (Enhancing NeuroImaging Genetics through Meta Analysis) Consortium, a global alliance of over 1400 scientists across 43 countries, studying the human brain in health and disease. Building on large-scale genetic studies that discovered the first robustly replicated genetic loci associated with brain metrics, ENIGMA has diversified into over 50 working groups (WGs), pooling worldwide data and expertise to answer fundamental questions in neuroscience, psychiatry, neurology, and genetics. Most ENIGMA WGs focus on specific psychiatric and neurological conditions, other WGs study normal variation due to sex and gender differences, or development and aging; still other WGs develop methodological pipelines and tools to facilitate harmonized analyses of "big data" (i.e., genetic and epigenetic data, multimodal MRI, and electroencephalography data). These international efforts have yielded the largest neuroimaging studies to date in schizophrenia, bipolar disorder, major depressive disorder, post-traumatic stress disorder, substance use disorders, obsessive-compulsive disorder, attention-deficit/hyperactivity disorder, autism spectrum disorders, epilepsy, and 22q11.2 deletion syndrome. More recent ENIGMA WGs have formed to study anxiety disorders, suicidal thoughts and behavior, sleep and insomnia, eating disorders, irritability, brain injury, antisocial personality and conduct disorder, and dissociative identity disorder. Here, we summarize the first decade of ENIGMA's activities and ongoing projects, and describe the successes and challenges encountered along the way. We highlight the advantages of collaborative large-scale coordinated data analyses for testing reproducibility and robustness of findings, offering the opportunity to identify brain systems involved in clinical syndromes across diverse samples and associated genetic, environmental, demographic, cognitive, and psychosocial factors.

Comparison of shape alterations of the thalamus and caudate nucleus between drug-naïve major depressive disorder patients and healthy controls

BACKGROUND

Although structural alterations have been reported in patients with major depressive disorder (MDD), very few studies have compared the shape alterations of the subcortical regions between drug-naïve MDD patients and healthy controls (HCs). Therefore, we investigated and compared the subcortical shape alterations and volumetric changes between drug-naïve MDD patients and HCs in this study.

METHODS

This study included 45 drug-naïve MDD patients and 83 HCs, who underwent three-dimensional (3-D) T1-weighted structural magnetic resonance imaging. Surface-based vertex analysis (SVA) was performed with automated segmentation of the bilateral caudate nuclei, putamina, nuclei accumbens, thalami, pallidum, hippocampi, amygdalae, and brainstem. SVA revealed regional contractions of the thalamus (bilateral medial and lateral nuclei) and right caudate nucleus (medial wall and anterosuperior areas) in the drug-naïve MDD patients when compared to HCs RESULTS: In volume analysis, the drug-naïve MDD patients showed a significant decrease in the volume of bilateral thalami compared with HCs (after Bonferroni correction p < 0.003). We identified morphometric contractions in bilateral thalami and right caudate nucleus in the drug-naïve MDD patients (p < 0.05).

CONCLUSIONS

The present study implied that with cortical shape changes, the subcortical brain alterations could contribute to emotional dysregulation in the drug-naïve MDD patients.

Genome wide association study of incomplete hippocampal inversion in adolescents

Incomplete hippocampal inversion (IHI), also called hippocampal malrotation, is an atypical presentation of the hippocampus present in about 20% of healthy individuals. Here we conducted the first genome-wide association study (GWAS) in IHI to elucidate the genetic underpinnings that may contribute to the incomplete inversion during brain development. A total of 1381 subjects contributed to the discovery cohort obtained from the IMAGEN database. The incidence rate of IHI was 26.1%. Loci with P<1e-5 were followed up in a validation cohort comprising 161 subjects from the PING study. Summary statistics from the discovery cohort were used to compute IHI heritability as well as genetic correlations with other traits. A locus on 18q11.2 (rs9952569; OR = 1.999; Z = 5.502; P = 3.755e-8) showed a significant association with the presence of IHI. A functional annotation of the locus implicated genes AQP4 and KCTD1. However, neither this locus nor the other 16 suggestive loci reached a significant p-value in the validation cohort. The h² estimate was 0.54 (sd: 0.30) and was significant (Z = 1.8; P = 0.036). The top three genetic correlations of IHI were with traits representing either intelligence or education attainment and reached nominal P< = 0.013.

Amygdala Nuclei Volume and Shape in Military Veterans With Posttraumatic Stress Disorder

BACKGROUND

The amygdala is a subcortical structure involved in socioemotional and associative fear learning processes relevant for understanding the mechanisms of posttraumatic stress disorder (PTSD). Research in animals indicates that the amygdala is a heterogeneous structure in which the basolateral and centromedial divisions are susceptible to stress. While the amygdala complex is implicated in the pathophysiology of PTSD, little is known about the specific contributions of the individual nuclei that constitute the amygdala complex.

METHODS

Military veterans (n = 355), including military veterans with PTSD (n = 149) and trauma-exposed control subjects without PTSD (n = 206), underwent high-resolution T1-weighted anatomical scans. Automated FreeSurfer segmentation of the amygdala yielded 9 structures: basal, lateral, accessory basal, anterior amygdaloid, and central, medial, cortical, and paralaminar nuclei, along with the corticoamygdaloid transition zone. Subregional volumes were compared between groups using ordinary-least-squares regression with relevant demographic and clinical regressors followed by 3-dimensional shape analysis of whole amygdala.

RESULTS

PTSD was associated with smaller left and right lateral and paralaminar nuclei, but with larger left and right central, medial, and cortical nuclei (p < .05, false discovery rate corrected). Shape analyses revealed lower radial distance in anterior bilateral amygdala and lower Jacobian determinant in posterior bilateral amygdala in PTSD compared with control subjects.

CONCLUSIONS

Alterations in select amygdala subnuclear volumes and regional shape distortions are associated with PTSD in military veterans. Volume differences of the lateral nucleus and the centromedial complex associated with PTSD demonstrate a subregion-specific pattern that is consistent with their functional roles in fear learning and fear expression behaviors.

Polygenic burden associated to oligodendrocyte precursor cells and radial glia influences the hippocampal volume changes induced by aerobic exercise in schizophrenia patients

Hippocampal volume decrease is a structural hallmark of schizophrenia (SCZ), and convergent evidence from postmortem and imaging studies suggests that it may be explained by changes in the cytoarchitecture of the cornu ammonis 4 (CA4) and dentate gyrus (DG) subfields. Increasing evidence indicates that aerobic exercise increases hippocampal volume in CA subfields and improves cognition in SCZ patients. Previous studies showed that the effects of exercise on the hippocampus might be connected to the polygenic burden of SCZ risk variants. However, little is known about cell type-specific genetic contributions to these structural changes. In this secondary analysis, we evaluated the modulatory role of cell type-specific SCZ polygenic risk scores (PRS) on volume changes in the CA1, CA2/3, and CA4/DG subfields over time. We studied 20 multi-episode SCZ patients and 23 healthy controls who performed aerobic exercise, and 21 multi-episode SCZ patients allocated to a control intervention (table soccer) for 3 months. Magnetic resonance imaging-based assessments were performed with FreeSurfer at baseline and after 3 months. The analyses showed that the polygenic burden associated with oligodendrocyte precursor cells (OPC) and radial glia (RG) significantly influenced the volume changes between baseline and 3 months in the CA4/DG subfield in SCZ patients performing aerobic exercise. A higher OPC- or RG-associated genetic risk burden was associated with a less pronounced volume increase or even a decrease in CA4/DG during the exercise intervention. We hypothesize that SCZ cell type-specific polygenic risk modulates the aerobic exercise-induced neuroplastic processes in the hippocampus.

Evidence for cerebello-thalamo-cortical hyperconnectivity as a heritable trait for schizophrenia

Our recent study has demonstrated that increased connectivity in the cerebello-thalamo-cortical (CTC) circuitry is a state-independent neural trait that can potentially predict the onset of psychosis. One possible cause of such "trait" abnormality would be genetic predisposition. Here, we tested this hypothesis using multi-paradigm functional magnetic resonance imaging (fMRI) data from two independent twin cohorts. In a sample of 85 monozygotic (MZ) and 52 dizygotic (DZ) healthy twin pairs acquired from the Human Connectome Project, we showed that the connectivity pattern of the identified CTC circuitry was more similar in the MZ twins (r = 0.54) compared with that in the DZ twins (r = 0.22). The structural equation modeling analysis revealed a heritability estimate of 0.52 for the CTC connectivity, suggesting a moderately strong genetic effect. Moreover, using an independent schizophrenia cotwin sample (10 discordant MZ cotwins, 30 discordant DZ cotwins, and 32 control cotwins), we observed a significant linear relationship between genetic distance to schizophrenia and the connectivity strength in the CTC circuitry (i.e., schizophrenia MZ cotwins > schizophrenia DZ cotwins > control twins, P = 0.045). The present data provide converging evidence that increased connectivity in the CTC circuitry is likely to be a heritable trait that is associated with the genetic risk of schizophrenia.

High-Dimensional Mapping of Cognition to the Brain Using Voxel-Based Morphometry and Subcortical Shape Analysis

BACKGROUND

It is increasingly recognized that the complex functions of human cognition are not accurately represented by arbitrarily-defined anatomical brain regions. Given the considerable functional specialization within such regions, more fine-grained studies of brain structure could capture such localized associations. However, such analyses/studies in a large community-dwelling population are lacking.

OBJECTIVE

To perform a fine-mapping of cognitive ability to cortical and subcortical grey matter on magnetic resonance imaging (MRI).

METHODS

In 3,813 stroke-free and non-demented persons from the Rotterdam Study (mean age 69.1 (±8.8) years; 55.8% women) with cognitive assessments and brain MRI, we performed voxel-based morphometry and subcortical shape analysis on global cognition and separate tests that tapped into memory, information processing speed, fine motor speed, and executive function domains.

RESULTS

We found that the different cognitive tests significantly associated with grey matter density in differential but also overlapping brain regions, primarily in the left hemisphere. Clusters of significantly associated voxels with global cognition were located within multiple anatomic regions: left amygdala, hippocampus, parietal lobule, superior temporal gyrus, insula and posterior temporal lobe. Subcortical shape analysis revealed associations primarily within the head and tail of the caudate nucleus, putamen, ventral part of the thalamus, and nucleus accumbens, more equally distributed among the left and right hemisphere. Within the caudate nucleus both positive (head) as well as negative (tail) associations were observed with global cognition.

CONCLUSIONS

In a large population-based sample, we mapped cognitive performance to cortical and subcortical grey matter density using a hypothesis-free approach with high-dimensional neuroimaging. Leveraging the power of our large sample size, we confirmed well-known associations as well as identified novel brain regions related to cognition.

Genetic contribution to the phenotypic correlation between trait impulsivity and resting-state functional connectivity of the amygdala and its subregions

Trait impulsivity, a predisposition to respond to stimuli without regard for the potentially negative consequences, contributes to many maladaptive behaviors. Studies have shown that both genetic factors and interregional functional interactions underlie trait impulsivity. However, whether common genes contribute to both trait impulsivity and its neural basis is still unknown. This study investigated the phenotypic correlations between trait impulsivity and the resting-state functional connectivity (rsFC) of the amygdala as well as its subregions and the genetic contribution to the phenotypic correlations. By recruiting a sample of 292 twins in late adolescence and young adulthood, we found that trait impulsivity was positively correlated with the rsFC between the left full amygdala and the right dorsolateral prefrontal cortex (DLPFC). Further analyses on the subregions of the amygdala showed that trait impulsivity was positively correlated with the rsFCs between the left basolateral (BL) amygdala and both the right DLPFC and the right inferior frontal gyrus and with the rsFCs between the right superficial (SF) amygdala and both the dorsal anterior cingulate cortex and right anterior insula. Bivariate genetic modelling analyses found genetic overlaps between trait impulsivity and the rsFC of the left full amygdala or the left BL amygdala with the right DLPFC. The proportions of phenotypic associations accounted for by overlapping genes were 82% and 60%, respectively. These results provide evidence for the genetic overlap between trait impulsivity and the intrinsic brain functional connectivity centered at the amygdala and especially at its BL subregion.

Beyond heritability: improving discoverability in imaging genetics

Structural neuroimaging measures based on magnetic resonance imaging have been at the forefront of imaging genetics. Global efforts to ensure homogeneity of measurements across study sites have enabled large-scale imaging genetic projects, accumulating nearly 50K samples for genome-wide association studies (GWAS). However, not many novel genetic variants have been identified by these GWAS, despite the high heritability of structural neuroimaging measures. Here, we discuss the limitations of using heritability as a guidance for assessing statistical power of GWAS, and highlight the importance of discoverability-which is the power to detect genetic variants for a given phenotype depending on its unique genomic architecture and GWAS sample size. Further, we present newly developed methods that boost genetic discovery in imaging genetics. By redefining imaging measures independent of traditional anatomical conventions, it is possible to improve discoverability, enabling identification of more genetic effects. Moreover, by leveraging enrichment priors from genomic annotations and independent GWAS of pleiotropic traits, we can better characterize effect size distributions, and identify reliable and replicable loci associated with structural neuroimaging measures. Statistical tools leveraging novel insights into the genetic discoverability of human traits, promises to accelerate the identification of genetic underpinnings underlying brain structural variation.

Subcortical brain volumes, cortical thickness and cortical surface area in families genetically enriched for social anxiety disorder - A multiplex multigenerational neuroimaging study

BACKGROUND

Social anxiety disorder (SAD) is a disabling psychiatric condition with a genetic background. Brain alterations in gray matter (GM) related to SAD have been previously reported, but it remains to be elucidated whether GM measures are candidate endophenotypes of SAD. Endophenotypes are measurable characteristics on the causal pathway from genotype to phenotype, providing insight in genetically-based disease mechanisms. Based on a review of existing evidence, we examined whether GM characteristics meet two endophenotype criteria, using data from a unique sample of SAD-patients and their family-members of two generations. First, we investigated whether GM characteristics co-segregate with social anxiety within families genetically enriched for SAD. Secondly, heritability of the GM characteristics was estimated.

METHODS

Families with a genetic predisposition for SAD participated in the Leiden Family Lab study on SAD; T1-weighted MRI brain scans were acquired (n = 110, 8 families). Subcortical volumes, cortical thickness and cortical surface area were determined for a-priori determined regions of interest (ROIs). Next, associations with social anxiety and heritabilities were estimated.

FINDINGS

Several subcortical and cortical GM characteristics, derived from frontal, parietal and temporal ROIs, co-segregated with social anxiety within families (uncorrected p-level) and showed moderate to high heritability.

INTERPRETATION

These findings provide preliminary evidence that GM characteristics of multiple ROIs, which are distributed over the brain, are candidate endophenotypes of SAD. Thereby, they shed light on the genetic vulnerability for SAD. Future research is needed to confirm these results and to link them to functional brain alterations and to genetic variations underlying these GM changes. FUND: Leiden University Research Profile 'Health, Prevention and the Human Life Cycle'.

Heritability and genome-wide associations studies of cerebral blood flow in the general population

Cerebral blood flow is an important process for brain functioning and its dysregulation is implicated in multiple neurological disorders. While environmental risk factors have been identified, it remains unclear to what extent the flow is regulated by genetics. Here we performed heritability and genome-wide association analyses of cerebral blood flow in a population-based cohort study. We included 4472 persons free of cortical infarcts who underwent genotyping and phase-contrast magnetic resonance flow imaging (mean age 64.8 ± 10.8 years). The flow rate, cross-sectional area of the vessel, and flow velocity through the vessel were measured in the basilar artery and bilateral carotids. We found that the flow rate of the basilar artery is most heritable (h2 (SE) = 24.1 (9.8), p-value = 0.0056), and this increased over age. The association studies revealed two significant loci for the right carotid artery area (rs12546630, p-value = 2.0 × 10^-8) and velocity (rs2971609, p-value = 1.4 × 10^-8), with the latter showing a concordant effect in an independent sample (N = 1350, p-value = 0.057, meta-analyzed p-value = 2.5 × 10^-9). These loci were also associated with other cerebral blood flow parameters below genome-wide significance, and rs2971609 lies in a known migraine locus. These findings establish that cerebral blood flow is under genetic control with potential relevance for neurological diseases.

ESTIMATION OF SHAPE AND GROWTH BRAIN NETWORK ATLASES FOR CONNECTOMIC BRAIN MAPPING IN DEVELOPING INFANTS

In vivo brain connectomics have heavily relied on using functional and diffusion Magnetic Resonance Imaging (MRI) modalities to examine functional and structural relationships between pairs of anatomical regions in the brain. However, research work on brain morphological (i.e., shape-to-shape) connections, which can be derived from T1-w and T2-w MR images, in both typical and atypical development or ageing is very scarce. Furthermore, the brain cannot be only regarded as a static shape, since it is a dynamic complex system that changes at functional, structural and morphological levels. Hence, examining the 'connection' between brain shape and its changes with time (e.g., growth) may help advance our understanding of connectomic brain dynamics as well as disorders that may affect it. To address these limitations, we unprecedentedly introduce two population-based shape and growth connectivity analysis tools that further extend the field of connectomics to brain morphology and dynamics: the morphome and the kinectome. Specifically, for a population of anatomically labelled shapes, the morphome identifies a network of anatomical shape regions that are connected when morphologically similar at a single timepoint, whereas the kinectome identifies anatomical shape regions that elicit similar evolution dynamics across successive timepoints. These proposed generic tools can be easily invested to examine how a baseline shape influences its deformation trajectory at later timepoints using any longitudinal shape data. We evaluated these tools on 23 infants, with right and left cortical surfaces reconstructed at birth, 3, 6, 9 and 12 months of age. Investigating the relationship between the neonatal morphome and the postnatal kinectome (from birth to 1 year of age) gave insights into brain connectivity at birth and how it develops over time.

Changes of arthropod diversity across an altitudinal ecoregional zonation in Northwestern Argentina

This study examined arthropod community patterns over an altitudinal ecoregional zonation that extended through three ecoregions (Yungas, Monte de Sierras y Bolsones, and Puna) and two ecotones (Yungas-Monte and Prepuna) of Northwestern Argentina (altitudinal range of 2,500 m), and evaluated the abiotic and biotic factors and the geographical distance that could influence them. Pitfall trap and suction samples were taken seasonally in 15 sampling sites (1,500-4,000 m a.s.l) during one year. In addition to climatic variables, several soil and vegetation variables were measured in the field. Values obtained for species richness between ecoregions and ecotones and by sampling sites were compared statistically and by interpolation-extrapolation analysis based on individuals at the same sample coverage level. Effects of predictor variables and the similarity of arthropods were shown using non-metric multidimensional scaling, and the resulting groups were evaluated using a multi-response permutation procedure. Polynomial regression was used to evaluate the relationship between altitude with total species richness and those of hyperdiverse/abundant higher taxa and the latter taxa with each predictor variable. The species richness pattern displayed a decrease in species diversity as the elevation increased at the bottom wet part (Yungas) of our altitudinal zonation until the Monte, and a unimodal pattern of diversity in the top dry part (Monte, Puna). Each ecoregion and ecotonal zone evidenced a particular species richness and assemblage of arthropods, but the latter ones displayed a high percentage of species shared with the adjacent ecoregions. The arthropod elevational pattern and the changes of the assemblages were explained by the environmental gradient (especially the climate) in addition to a geographic gradient (the distance of decay of similarity), demonstrating that the species turnover is important to explain the beta diversity along the elevational gradient. This suggests that patterns of diversity and distribution of arthropods are regulated by the dissimilarity of ecoregional environments that establish a wide range of geographic and environmental barriers, coupled with a limitation of species dispersal. Therefore, the arthropods of higher taxa respond differently to the altitudinal ecoregional zonation.

Neuroimaging genomics in psychiatry-a translational approach

Neuroimaging genomics is a relatively new field focused on integrating genomic and imaging data in order to investigate the mechanisms underlying brain phenotypes and neuropsychiatric disorders. While early work in neuroimaging genomics focused on mapping the associations of candidate gene variants with neuroimaging measures in small cohorts, the lack of reproducible results inspired better-powered and unbiased large-scale approaches. Notably, genome-wide association studies (GWAS) of brain imaging in thousands of individuals around the world have led to a range of promising findings. Extensions of such approaches are now addressing epigenetics, gene–gene epistasis, and gene–environment interactions, not only in brain structure, but also in brain function. Complementary developments in systems biology might facilitate the translation of findings from basic neuroscience and neuroimaging genomics to clinical practice. Here, we review recent approaches in neuroimaging genomics—we highlight the latest discoveries, discuss advantages and limitations of current approaches, and consider directions by which the field can move forward to shed light on brain disorders.