Miguel E. Rentería

Subcortical brain alterations in major depressive disorder: findings from the ENIGMA Major Depressive Disorder working group.

The pattern of structural brain alterations associated with major depressive disorder (MDD) remains unresolved. This is in part due to small sample sizes of neuroimaging studies resulting in limited statistical power, disease heterogeneity and the complex interactions between clinical characteristics and brain morphology. To address this, we meta-analyzed three-dimensional brain magnetic resonance imaging data from 1728 MDD patients and 7199 controls from 15 research samples worldwide, to identify subcortical brain volumes that robustly discriminate MDD patients from healthy controls. Relative to controls, patients had significantly lower hippocampal volumes (Cohen’s d=−0.14, % difference=−1.24). This effect was driven by patients with recurrent MDD (Cohen’s d=−0.17, % difference=−1.44), and we detected no differences between first episode patients and controls. Age of onset ⩽21 was associated with a smaller hippocampus (Cohen’s d=−0.20, % difference=−1.85) and a trend toward smaller amygdala (Cohen’s d=−0.11, % difference=−1.23) and larger lateral ventricles (Cohen’s d=0.12, % difference=5.11). Symptom severity at study inclusion was not associated with any regional brain volumes. Sample characteristics such as mean age, proportion of antidepressant users and proportion of remitted patients, and methodological characteristics did not significantly moderate alterations in brain volumes in MDD. Samples with a higher proportion of antipsychotic medication users showed larger caudate volumes in MDD patients compared with controls. This currently largest worldwide effort to identify subcortical brain alterations showed robust smaller hippocampal volumes in MDD patients, moderated by age of onset and first episode versus recurrent episode status.

Cortical abnormalities in adults and adolescents with major depression based on brain scans from 20 cohorts worldwide in the ENIGMA Major Depressive Disorder Working Group

The neuro-anatomical substrates of major depressive disorder (MDD) are still not well understood, despite many neuroimaging studies over the past few decades. Here we present the largest ever worldwide study by the ENIGMA (Enhancing Neuro Imaging Genetics through Meta-Analysis) Major Depressive Disorder Working Group on cortical structural alterations in MDD. Structural T1-weighted brain magnetic resonance imaging (MRI) scans from 2148 MDD patients and 7957 healthy controls were analysed with harmonized protocols at 20 sites around the world. To detect consistent effects of MDD and its modulators on cortical thickness and surface area estimates derived from MRI, statistical effects from sites were meta-analysed separately for adults and adolescents. Adults with MDD had thinner cortical gray matter than controls in the orbitofrontal cortex (OFC), anterior and posterior cingulate, insula and temporal lobes (Cohen’s d effect sizes: −0.10 to −0.14). These effects were most pronounced in first episode and adult-onset patients (>21 years). Compared to matched controls, adolescents with MDD had lower total surface area (but no differences in cortical thickness) and regional reductions in frontal regions (medial OFC and superior frontal gyrus) and primary and higher-order visual, somatosensory and motor areas (d: −0.26 to −0.57). The strongest effects were found in recurrent adolescent patients. This highly powered global effort to identify consistent brain abnormalities showed widespread cortical alterations in MDD patients as compared to controls and suggests that MDD may impact brain structure in a highly dynamic way, with different patterns of alterations at different stages of life.

Cerebral asymmetry: a quantitative, multifactorial, and plastic brain phenotype.

The longitudinal fissure separates the human brain into two hemispheres that remain connected through the corpus callosum. The left and the right halves of the brain resemble each other, and almost every structure present in one side has an equivalent structure in the other. Despite this exceptional correspondence, the two hemispheres also display important anatomical differences and there is marked lateralization of certain cognitive and motor functions such as language and handedness. However, the mechanisms that underlie the establishment of these hemispheric specializations, as well as their physiological and behavioral implications, remain largely unknown. Thanks to recent advances in neuroimaging, a series of studies documenting variation in symmetry and asymmetry as a function of age, gender, brain region, and pathological state, have been published in the past decade. Here, we review evidence of normal and atypical cerebral asymmetry, and the factors that influence it at the macrostructural level. Given the prominent role that cerebral asymmetry plays in the organization of the brain, and its possible implication in neurodevelopmental and psychiatric conditions, further research in this area is anticipated.

A comparative structural bioinformatics analysis of the insulin receptor family ectodomain based on phylogenetic information

The insulin receptor (IR), the insulin-like growth factor 1 receptor (IGF1R) and the insulin receptor-related receptor (IRR) are covalently-linked homodimers made up of several structural domains. The molecular mechanism of ligand binding to the ectodomain of these receptors and the resulting activation of their tyrosine kinase domain is still not well understood. We have carried out an amino acid residue conservation analysis in order to reconstruct the phylogeny of the IR Family. We have confirmed the location of ligand binding site 1 of the IGF1R and IR. Importantly, we have also predicted the likely location of the insulin binding site 2 on the surface of the fibronectin type III domains of the IR. An evolutionary conserved surface on the second leucine-rich domain that may interact with the ligand could not be detected. We suggest a possible mechanical trigger of the activation of the IR that involves a slight ‘twist’ rotation of the last two fibronectin type III domains in order to face the likely location of insulin. Finally, a strong selective pressure was found amongst the IRR orthologous sequences, suggesting that this orphan receptor has a yet unknown physiological role which may be conserved from amphibians to mammals.

Genetic architecture of subcortical brain regions: common and region-specific genetic contributions

Understanding the aetiology of patterns of variation within and covariation across brain regions is key to advancing our understanding of the functional, anatomical and developmental networks of the brain. Here we applied multivariate twin modelling and principal component analysis (PCA) to investigate the genetic architecture of the size of seven subcortical regions (caudate nucleus, thalamus, putamen, pallidum, hippocampus, amygdala and nucleus accumbens) in a genetically informative sample of adolescents and young adults (N = 1038; mean age = 21.6 ± 3.2 years; including 148 monozygotic and 202 dizygotic twin pairs) from the Queensland Twin IMaging (QTIM) study. Our multivariate twin modelling identified a common genetic factor that accounts for all the heritability of intracranial volume (0.88) and a substantial proportion of the heritability of all subcortical structures, particularly those of the thalamus (0.71 out of 0.88), pallidum (0.52 out of 0.75) and putamen (0.43 out of 0.89). In addition, we also found substantial region‐specific genetic contributions to the heritability of the hippocampus (0.39 out of 0.79), caudate nucleus (0.46 out of 0.78), amygdala (0.25 out of 0.45) and nucleus accumbens (0.28 out of 0.52). This provides further insight into the extent and organization of subcortical genetic architecture, which includes developmental and general growth pathways, as well as the functional specialization and maturation trajectories that influence each subcortical region.

LocusTrack: integrated visualization of GWAS results and genomic annotation

BackgroundGenome-wide association studies (GWAS) are an important tool for the mapping of complex traits and diseases. Visual inspection of genomic annotations may be used to generate insights into the biological mechanisms underlying GWAS-identified loci.ResultsWe developed LocusTrack, a web-based application that annotates and creates plots of regional GWAS results and incorporates user-specified tracks that display annotations such as linkage disequilibrium (LD), phylogenetic conservation, chromatin state, and other genomic and regulatory elements. Currently, LocusTrack can integrate annotation tracks from the UCSC genome-browser as well as from any tracks provided by the user.ConclusionLocusTrack is an easy-to-use application and can be accessed at the following URL: http://gump.qimr.edu.au/general/gabrieC/LocusTrack/. Users can upload and manage GWAS results and select from and/or provide annotation tracks using simple and intuitive menus. LocusTrack scripts and associated data can be downloaded from the website and run locally.

Using PLINK for genome-wide association studies (GWAS) and data analysis

Within this chapter we introduce the basic PLINK functions for reading in data, applying quality control, and running association analyses. Three worked examples are provided to illustrate: data management and assessment of population substructure, association analysis of a quantitative trait, and qualitative or case-control association analyses.

Subcortical brain structure and suicidal behaviour in major depressive disorder: a meta-analysis from the ENIGMA-MDD working group

The aetiology of suicidal behaviour is complex, and knowledge about its neurobiological mechanisms is limited. Neuroimaging methods provide a noninvasive approach to explore the neural correlates of suicide vulnerability in vivo. The ENIGMA-MDD Working Group is an international collaboration evaluating neuroimaging and clinical data from thousands of individuals collected by research groups from around the world. Here we present analyses in a subset sample (n=3097) for whom suicidality data were available. Prevalence of suicidal symptoms among major depressive disorder (MDD) cases ranged between 29 and 69% across cohorts. We compared mean subcortical grey matter volumes, lateral ventricle volumes and total intracranial volume (ICV) in MDD patients with suicidal symptoms (N=451) vs healthy controls (N=1996) or MDD patients with no suicidal symptoms (N=650). MDD patients reporting suicidal plans or attempts showed a smaller ICV (P=4.12 × 10−3) or a 2.87% smaller volume compared with controls (Cohen’s d=−0.284). In addition, we observed a nonsignificant trend in which MDD cases with suicidal symptoms had smaller subcortical volumes and larger ventricular volumes compared with controls. Finally, no significant differences (P=0.28–0.97) were found between MDD patients with and those without suicidal symptoms for any of the brain volume measures. This is by far the largest neuroimaging meta-analysis of suicidal behaviour in MDD to date. Our results did not replicate previous reports of association between subcortical brain structure and suicidality and highlight the need for collecting better-powered imaging samples and using improved suicidality assessment instruments.

GWAS of DNA methylation variation within imprinting control regions suggests parent-of-origin association.

Imprinting control regions (ICRs) play a fundamental role in establishing and maintaining the non-random monoallelic expression of certain genes, via common regulatory elements such as non-coding RNAs and differentially methylated regions (DMRs) of DNA. We recently surveyed DNA methylation levels within four ICRs (H19-ICR, IGF2-DMR, KvDMR, and NESPAS-ICR) in whole-blood genomic DNA from 128 monozygotic (MZ) and 128 dizygotic (DZ) human twin pairs. Our analyses revealed high individual variation and intra-domain covariation in methylation levels across CpGs and emphasized the interaction between epigenetic variation and the underlying genetic sequence in a parent-of-origin fashion. Here, we extend our analysis to conduct two genome-wide screenings of single nucleotide polymorphisms (SNPs) underlying either intra-domain covariation or parent-of-origin-dependent association with methylation status at individual CpG sites located within ICRs. Although genome-wide significance was not surpassed due to sample size limitations, the most significantly associated SNPs found through multiple-trait genome-wide association (MQFAM) included the previously described rs10732516, which is located in the vicinity of the H19-ICR. Similarly, we identified an association between rs965808 and methylation status within the NESPAS-ICR. This SNP is positioned within an intronic region of the overlapping genes GNAS and GNAS-AS1, which are imprinted genes regulated by the NESPAS-ICR. Sixteen other SNPs located in regions apart from the analyzed regions displayed suggestive association with intra-domain methylation. Additionally, we identified 13 SNPs displaying parent-of-origin association with individual methylation sites through family-based association testing. In this exploratory study, we show the value and feasibility of using alternative GWAS approaches in the study of the interaction between epigenetic state and genetic sequence within imprinting regulatory domains. Despite the relatively small sample size, we identified a number of SNPs displaying suggestive association either in a domain-wide or in a parent-of-origin fashion. Nevertheless, these associations will require future experimental validation or replication in larger and independent samples.

Early developmental gene enhancers affect subcortical volumes in the adult human brain.

Genome‐wide association screens aim to identify common genetic variants contributing to the phenotypic variability of complex traits, such as human height or brain morphology. The identified genetic variants are mostly within noncoding genomic regions and the biology of the genotype–phenotype association typically remains unclear. In this article, we propose a complementary targeted strategy to reveal the genetic underpinnings of variability in subcortical brain volumes, by specifically selecting genomic loci that are experimentally validated forebrain enhancers, active in early embryonic development. We hypothesized that genetic variation within these enhancers may affect the development and ultimately the structure of subcortical brain regions in adults. We tested whether variants in forebrain enhancer regions showed an overall enrichment of association with volumetric variation in subcortical structures of >13,000 healthy adults. We observed significant enrichment of genomic loci that affect the volume of the hippocampus within forebrain enhancers (empirical P = 0.0015), a finding which robustly passed the adjusted threshold for testing of multiple brain phenotypes (cutoff of P < 0.0083 at an alpha of 0.05). In analyses of individual single nucleotide polymorphisms (SNPs), we identified an association upstream of the ID2 gene with rs7588305 and variation in hippocampal volume. This SNP‐based association survived multiple‐testing correction for the number of SNPs analyzed but not for the number of subcortical structures. Targeting known regulatory regions offers a way to understand the underlying biology that connects genotypes to phenotypes, particularly in the context of neuroimaging genetics. This biology‐driven approach generates testable hypotheses regarding the functional biology of identified associations. Hum Brain Mapp 37:1788–1800, 2016.