Genetic variants for head size share genes and pathways with cancer

The size of the human head is highly heritable, but genetic drivers of its variation within the general population remain unmapped. We perform a genome-wide association study on head size (N = 80,890) and identify 67 genetic loci, of which 50 are novel. Neuroimaging studies show that 17 variants affect specific brain areas, but most have widespread effects. Gene set enrichment is observed for various cancers and the p53, Wnt, and ErbB signaling pathways. Genes harboring lead variants are enriched for macrocephaly syndrome genes (37-fold) and high-fidelity cancer genes (9-fold), which is not seen for human height variants. Head size variants are also near genes preferentially expressed in intermediate progenitor cells, neural cells linked to evolutionary brain expansion. Our results indicate that genes regulating early brain and cranial growth incline to neoplasia later in life, irrespective of height. This warrants investigation of clinical implications of the link between head size and cancer.

Beyond the Global Brain Differences: Intraindividual Variability Differences in 1q21.1 Distal and 15q11.2 BP1-BP2 Deletion Carriers

BACKGROUND

Carriers of the 1q21.1 distal and 15q11.2 BP1-BP2 copy number variants exhibit regional and global brain differences compared with noncarriers. However, interpreting regional differences is challenging if a global difference drives the regional brain differences. Intraindividual variability measures can be used to test for regional differences beyond global differences in brain structure.

METHODS

Magnetic resonance imaging data were used to obtain regional brain values for 1q21.1 distal deletion (n = 30) and duplication (n = 27) and 15q11.2 BP1-BP2 deletion (n = 170) and duplication (n = 243) carriers and matched noncarriers (n = 2350). Regional intra-deviation scores, i.e., the standardized difference between an individual's regional difference and global difference, were used to test for regional differences that diverge from the global difference.

RESULTS

For the 1q21.1 distal deletion carriers, cortical surface area for regions in the medial visual cortex, posterior cingulate, and temporal pole differed less and regions in the prefrontal and superior temporal cortex differed more than the global difference in cortical surface area. For the 15q11.2 BP1-BP2 deletion carriers, cortical thickness in regions in the medial visual cortex, auditory cortex, and temporal pole differed less and the prefrontal and somatosensory cortex differed more than the global difference in cortical thickness.

CONCLUSIONS

We find evidence for regional effects beyond differences in global brain measures in 1q21.1 distal and 15q11.2 BP1-BP2 copy number variants. The results provide new insight into brain profiling of the 1q21.1 distal and 15q11.2 BP1-BP2 copy number variants, with the potential to increase understanding of the mechanisms involved in altered neurodevelopment.

Relationships between neurotransmitter receptor densities and expression levels of their corresponding genes in the human hippocampus

Neurotransmitter receptors are key molecules in signal transmission, their alterations are associated with brain dysfunction. Relationships between receptors and their corresponding genes are poorly understood, especially in humans. We combined in vitro receptor autoradiography and RNA sequencing to quantify, in the same tissue samples (7 subjects), the densities of 14 receptors and expression levels of their corresponding 43 genes in the Cornu Ammonis (CA) and dentate gyrus (DG) of human hippocampus. Significant differences in receptor densities between both structures were found only for metabotropic receptors, whereas significant differences in RNA expression levels mostly pertained ionotropic receptors. Receptor fingerprints of CA and DG differ in shapes but have similar sizes; the opposite holds true for their "RNA fingerprints", which represent the expression levels of multiple genes in a single area. In addition, the correlation coefficients between receptor densities and corresponding gene expression levels vary widely and the mean correlation strength was weak-to-moderate. Our results suggest that receptor densities are not only controlled by corresponding RNA expression levels, but also by multiple regionally specific post-translational factors.

Associations of common genetic risk variants of the muscarinic acetylcholine receptor M2 with cardiac autonomic dysfunction in patients with schizophrenia

OBJECTIVES

Decreased vagal modulation, which has consistently been observed in schizophrenic patients, might contribute to increased cardiac mortality in schizophrenia. Previously, associations between CHRM2 (Cholinergic Receptor Muscarinic 2) and cardiac autonomic features have been reported. Here, we tested for possible associations between these polymorphisms and heart rate variability in patients with schizophrenia.

METHODS

A total of three single nucleotide polymorphisms (SNPs) in CHRM2 (rs73158705 A>G, rs8191992 T>A and rs2350782 T>C) that achieved significance (p < 5 * 10^-8) in genome-wide association studies for cardiac autonomic features were genotyped in 88 drug-naïve patients, 61 patients receiving antipsychotic medication and 144 healthy controls. Genotypes were analysed for associations with parameters of heart rate variability and complexity, in each diagnostic group.

RESULTS

We observed a significantly altered heart rate variability in unmedicated patients with identified genetic risk status in rs73158705 A>G, rs8191992 T>A and rs2350782 T>C as compared to genotype non-risk status. In patients receiving antipsychotic medication and healthy controls, these associations were not observed.

DISCUSSION

We report novel candidate genetic associations with cardiac autonomic dysfunction in schizophrenia, but larger cohorts are required for replication.

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

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

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

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

A GWAS top hit for circulating leptin is associated with weight gain but not with leptin protein levels in lithium-augmented patients with major depression

Lithium-treated patients often suffer from weight gain as a common adverse event. In an earlier investigation, we found an impact of two single-nucleotide polymorphisms (rs10487506 and rs2278815) at the leptin gene on weight gain but not on leptin protein levels in serum under lithium augmentation. A recent genome-wide association study identified a polymorphism at the leptin gene locus (rs10487505) associated with circulating leptin protein levels. To characterize potential effects of this variant in acute major depressive disorder, we investigated body mass indices from 180 lithium-augmented patients and serum concentrations of leptin protein from 89 patients using linear mixed model analyses and rs6979832, a proxy SNP of rs10487505. Body mass index was measured before and after 4 weeks of lithium augmentation, in a subsample also after 4 and 7 months. Leptin serum levels were measured before and during lithium augmentation. G-allele homozygotes of rs6979832 had a significantly lower body mass index increase during observation compared to A-allele hetero- and homozygotes. However, we found no influence on leptin serum levels. Joint analyses of rs6979832 with the previously investigated polymorphisms rs10487506 and rs2278815, and expressed quantitative trait data, suggest a complex interplay between SNP alleles at the leptin locus. These results strongly support our earlier findings that common genetic variation at the leptin gene locus may be involved in lithium augmentation-associated weight gain in major depressive disorder.

Identification of Phonology-Related Genes and Functional Characterization of Broca's and Wernicke's Regions in Language and Learning Disorders

Impaired phonological processing is a leading symptom of multifactorial language and learning disorders suggesting a common biological basis. Here we evaluated studies of dyslexia, dyscalculia, specific language impairment (SLI), and the logopenic variant of primary progressive aphasia (lvPPA) seeking for shared risk genes in Broca's and Wernicke's regions, being key for phonological processing within the complex language network. The identified "phonology-related genes" from literature were functionally characterized using Atlas-based expression mapping (JuGEx) and gene set enrichment. Out of 643 publications from the last decade until now, we extracted 21 candidate genes of which 13 overlapped with dyslexia and SLI, six with dyslexia and dyscalculia, and two with dyslexia, dyscalculia, and SLI. No overlap was observed between the childhood disorders and the late-onset lvPPA often showing symptoms of learning disorders earlier in life. Multiple genes were enriched in Gene Ontology terms of the topics learning (CNTNAP2, CYFIP1, DCDC2, DNAAF4, FOXP2) and neuronal development (CCDC136, CNTNAP2, CYFIP1, DCDC2, KIAA0319, RBFOX2, ROBO1). Twelve genes showed above-average expression across both regions indicating moderate-to-high gene activity in the investigated cortical part of the language network. Of these, three genes were differentially expressed suggesting potential regional specializations: ATP2C2 was upregulated in Broca's region, while DNAAF4 and FOXP2 were upregulated in Wernicke's region. ATP2C2 encodes a magnesium-dependent calcium transporter which fits with reports about disturbed calcium and magnesium levels for dyslexia and other communication disorders. DNAAF4 (formerly known as DYX1C1) is involved in neuronal migration supporting the hypothesis of disturbed migration in dyslexia. FOXP2 is a transcription factor that regulates a number of genes involved in development of speech and language. Overall, our interdisciplinary and multi-tiered approach provided evidence that genetic and transcriptional variation of ATP2C2, DNAAF4, and FOXP2 may play a role in physiological and pathological aspects of phonological processing.

Pathway-Specific Genetic Risk for Alzheimer's Disease Differentiates Regional Patterns of Cortical Atrophy in Older Adults

Brain aging is highly variable and represents a challenge to delimit aging from disease processes. Moreover, genetic factors may influence both aging and disease. Here we focused on this issue and investigated effects of multiple genetic loci previously identified to be associated with late-onset Alzheimer's disease (AD) on brain structure of older adults from a population sample. We calculated a genetic risk score (GRS) using genome-wide significant single-nucleotide polymorphisms from genome-wide association studies of AD and tested its effect on cortical thickness (CT). We observed a common pattern of cortical thinning (right inferior frontal, left posterior temporal, medial occipital cortex). To identify CT changes by specific biological processes, we subdivided the GRS effect according to AD-associated pathways and performed follow-up analyses. The common pattern from the main analysis was further differentiated by pathway-specific effects yielding a more bilateral pattern. Further findings were located in the superior parietal and mid/anterior cingulate regions representing 2 unique pathway-specific patterns. All patterns, except the superior parietal pattern, were influenced by apolipoprotein E. Our step-wise approach revealed atrophy patterns that partially resembled imaging findings in early stages of AD. Our study provides evidence that genetic burden for AD contributes to structural brain variability in normal aging.

Genome-wide association study of more than 40,000 bipolar disorder cases provides new insights into the underlying biology

Bipolar disorder is a heritable mental illness with complex etiology. We performed a genome-wide association study of 41,917 bipolar disorder cases and 371,549 controls of European ancestry, which identified 64 associated genomic loci. Bipolar disorder risk alleles were enriched in genes in synaptic signaling pathways and brain-expressed genes, particularly those with high specificity of expression in neurons of the prefrontal cortex and hippocampus. Significant signal enrichment was found in genes encoding targets of antipsychotics, calcium channel blockers, antiepileptics and anesthetics. Integrating expression quantitative trait locus data implicated 15 genes robustly linked to bipolar disorder via gene expression, encoding druggable targets such as HTR6, MCHR1, DCLK3 and FURIN. Analyses of bipolar disorder subtypes indicated high but imperfect genetic correlation between bipolar disorder type I and II and identified additional associated loci. Together, these results advance our understanding of the biological etiology of bipolar disorder, identify novel therapeutic leads and prioritize genes for functional follow-up studies.

1q21.1 distal copy number variants are associated with cerebral and cognitive alterations in humans

Low-frequency 1q21.1 distal deletion and duplication copy number variant (CNV) carriers are predisposed to multiple neurodevelopmental disorders, including schizophrenia, autism and intellectual disability. Human carriers display a high prevalence of micro- and macrocephaly in deletion and duplication carriers, respectively. The underlying brain structural diversity remains largely unknown. We systematically called CNVs in 38 cohorts from the large-scale ENIGMA-CNV collaboration and the UK Biobank and identified 28 1q21.1 distal deletion and 22 duplication carriers and 37,088 non-carriers (48% male) derived from 15 distinct magnetic resonance imaging scanner sites. With standardized methods, we compared subcortical and cortical brain measures (all) and cognitive performance (UK Biobank only) between carrier groups also testing for mediation of brain structure on cognition. We identified positive dosage effects of copy number on intracranial volume (ICV) and total cortical surface area, with the largest effects in frontal and cingulate cortices, and negative dosage effects on caudate and hippocampal volumes. The carriers displayed distinct cognitive deficit profiles in cognitive tasks from the UK Biobank with intermediate decreases in duplication carriers and somewhat larger in deletion carriers-the latter potentially mediated by ICV or cortical surface area. These results shed light on pathobiological mechanisms of neurodevelopmental disorders, by demonstrating gene dose effect on specific brain structures and effect on cognitive function.

A common variation in HCN1 is associated with heart rate variability in schizophrenia

BACKGROUND

There is growing evidence for a shared genetic basis between schizophrenia risk and cardiovascular disease. Reduced efferent vagal activity, indexed by reduced heart rate variability (HRV), has been consistently described in patients with schizophrenia and may potentially contribute to the increased cardiovascular risk in these patients. In this study, we tested the hypothesis whether the established schizophrenia risk variant HCN1 rs16902086 (A > G) is associated with reduced HRV.

METHODS

We analyzed the risk status of HCN1 rs16902086 (AG/GG vs. AA genotype) in 83 unmedicated patients with schizophrenia and 96 healthy controls and investigated genotype-related impacts on various HRV parameters.

RESULTS

We observed significantly increased resting heart rates and a marked decrease of vagal modulation in our patient cohort. Strikingly, HCN1 rs16902086 (A > G) was associated with reduced HRV parameters in patients only. A trend towards more pronounced HRV deviations was observed in homozygous (GG) compared to heterozygous patients (AG).

CONCLUSION

We present first evidence for a genetic risk factor that is associated with decreased vagal modulation in unmedicated patients with schizophrenia. Moreover, our findings suggest that HCN1 might be involved in reduced vagal modulation and possibly in increased cardiac mortality in schizophrenia patients. Thus, our data indicate that reduced vagal modulation might be an endophenotype of schizophrenia.

Association of Copy Number Variation of the 15q11.2 BP1-BP2 Region With Cortical and Subcortical Morphology and Cognition

Importance

Recurrent microdeletions and duplications in the genomic region 15q11.2 between breakpoints 1 (BP1) and 2 (BP2) are associated with neurodevelopmental disorders. These structural variants are present in 0.5% to 1.0% of the population, making 15q11.2 BP1-BP2 the site of the most prevalent known pathogenic copy number variation (CNV). It is unknown to what extent this CNV influences brain structure and affects cognitive abilities.

Objective

To determine the association of the 15q11.2 BP1-BP2 deletion and duplication CNVs with cortical and subcortical brain morphology and cognitive task performance.

Design, Setting, and Participants

In this genetic association study, T1-weighted brain magnetic resonance imaging were combined with genetic data from the ENIGMA-CNV consortium and the UK Biobank, with a replication cohort from Iceland. In total, 203 deletion carriers, 45 247 noncarriers, and 306 duplication carriers were included. Data were collected from August 2015 to April 2019, and data were analyzed from September 2018 to September 2019.

Main Outcomes and Measures

The associations of the CNV with global and regional measures of surface area and cortical thickness as well as subcortical volumes were investigated, correcting for age, age2, sex, scanner, and intracranial volume. Additionally, measures of cognitive ability were analyzed in the full UK Biobank cohort.

Results

Of 45 756 included individuals, the mean (SD) age was 55.8 (18.3) years, and 23 754 (51.9%) were female. Compared with noncarriers, deletion carriers had a lower surface area (Cohen d = -0.41; SE, 0.08; P = 4.9 × 10-8), thicker cortex (Cohen d = 0.36; SE, 0.07; P = 1.3 × 10-7), and a smaller nucleus accumbens (Cohen d = -0.27; SE, 0.07; P = 7.3 × 10-5). There was also a significant negative dose response on cortical thickness (β = -0.24; SE, 0.05; P = 6.8 × 10-7). Regional cortical analyses showed a localization of the effects to the frontal, cingulate, and parietal lobes. Further, cognitive ability was lower for deletion carriers compared with noncarriers on 5 of 7 tasks.

Conclusions and Relevance

These findings, from the largest CNV neuroimaging study to date, provide evidence that 15q11.2 BP1-BP2 structural variation is associated with brain morphology and cognition, with deletion carriers being particularly affected. The pattern of results fits with known molecular functions of genes in the 15q11.2 BP1-BP2 region and suggests involvement of these genes in neuronal plasticity. These neurobiological effects likely contribute to the association of this CNV with neurodevelopmental disorders.

The genetic architecture of the human cerebral cortex

The cerebral cortex underlies our complex cognitive capabilities, yet little is known about the specific genetic loci that influence human cortical structure. To identify genetic variants that affect cortical structure, we conducted a genome-wide association meta-analysis of brain magnetic resonance imaging data from 51,665 individuals. We analyzed the surface area and average thickness of the whole cortex and 34 regions with known functional specializations. We identified 199 significant loci and found significant enrichment for loci influencing total surface area within regulatory elements that are active during prenatal cortical development, supporting the radial unit hypothesis. Loci that affect regional surface area cluster near genes in Wnt signaling pathways, which influence progenitor expansion and areal identity. Variation in cortical structure is genetically correlated with cognitive function, Parkinson's disease, insomnia, depression, neuroticism, and attention deficit hyperactivity disorder.

Dose response of the 16p11.2 distal copy number variant on intracranial volume and basal ganglia

Carriers of large recurrent copy number variants (CNVs) have a higher risk of developing neurodevelopmental disorders. The 16p11.2 distal CNV predisposes carriers to e.g., autism spectrum disorder and schizophrenia. We compared subcortical brain volumes of 12 16p11.2 distal deletion and 12 duplication carriers to 6882 non-carriers from the large-scale brain Magnetic Resonance Imaging collaboration, ENIGMA-CNV. After stringent CNV calling procedures, and standardized FreeSurfer image analysis, we found negative dose-response associations with copy number on intracranial volume and on regional caudate, pallidum and putamen volumes (β = -0.71 to -1.37; P < 0.0005). In an independent sample, consistent results were obtained, with significant effects in the pallidum (β = -0.95, P = 0.0042). The two data sets combined showed significant negative dose-response for the accumbens, caudate, pallidum, putamen and ICV (P = 0.0032, 8.9 × 10-6, 1.7 × 10-9, 3.5 × 10-12 and 1.0 × 10-4, respectively). Full scale IQ was lower in both deletion and duplication carriers compared to non-carriers. This is the first brain MRI study of the impact of the 16p11.2 distal CNV, and we demonstrate a specific effect on subcortical brain structures, suggesting a neuropathological pattern underlying the neurodevelopmental syndromes.

Leptin gene polymorphisms are associated with weight gain during lithium augmentation in patients with major depression

Weight gain is a common adverse effect of lithium augmentation. Previous studies indicate an impact of genetic variants at the leptin gene on weight gain as a consequence of psychopharmacological treatment. The primary aim of our study was to identify variants at the leptin locus that might predict lithium-induced weight gain. The secondary aim was to investigate if these variants modulate leptin levels. In 180 patients with acute major depressive disorder, body mass index was measured before and after 4 weeks of lithium augmentation, in a subsample also after 4 and/or 7 months. In a subsample of 89 patients, leptin serum concentrations were measured before and during lithium augmentation. We used linear mixed model analyzes to investigate the effects of 2 polymorphisms at the leptin locus (rs4731426 and rs7799039, employing the respective proxy SNPs rs2278815 and rs10487506) on changes in body mass index and leptin levels. For both polymorphisms, which are in high linkage disequilibrium, body mass index was significantly lower in homozygous A-allele carriers than in carriers of other genotypes at baseline. Over the follow-up period, body mass index increased less in homozygous A-allele carriers of rs4731426 than in carriers of other genotypes. This was not the case for rs7799039. Neither polymorphism modulated leptin protein expression. Our study strongly supports the hypothesis that genetic variability at the leptin locus is involved in lithium augmentation-associated weight gain in major depressive disorder. Furthermore, Genotype-Tissue Expression data provide strong evidence that rs4731426 influences the expression of leptin messenger ribonucleic acid in fibroblasts.

Genetic architecture of subcortical brain structures in 38,851 individuals

Subcortical brain structures are integral to motion, consciousness, emotions and learning. We identified common genetic variation related to the volumes of the nucleus accumbens, amygdala, brainstem, caudate nucleus, globus pallidus, putamen and thalamus, using genome-wide association analyses in almost 40,000 individuals from CHARGE, ENIGMA and UK Biobank. We show that variability in subcortical volumes is heritable, and identify 48 significantly associated loci (40 novel at the time of analysis). Annotation of these loci by utilizing gene expression, methylation and neuropathological data identified 199 genes putatively implicated in neurodevelopment, synaptic signaling, axonal transport, apoptosis, inflammation/infection and susceptibility to neurological disorders. This set of genes is significantly enriched for Drosophila orthologs associated with neurodevelopmental phenotypes, suggesting evolutionarily conserved mechanisms. Our findings uncover novel biology and potential drug targets underlying brain development and disease.

Genome-wide association study identifies 30 loci associated with bipolar disorder

Bipolar disorder is a highly heritable psychiatric disorder. We performed a genome-wide association study (GWAS) including 20,352 cases and 31,358 controls of European descent, with follow-up analysis of 822 variants with P < 1 × 10-4 in an additional 9,412 cases and 137,760 controls. Eight of the 19 variants that were genome-wide significant (P < 5 × 10-8) in the discovery GWAS were not genome-wide significant in the combined analysis, consistent with small effect sizes and limited power but also with genetic heterogeneity. In the combined analysis, 30 loci were genome-wide significant, including 20 newly identified loci. The significant loci contain genes encoding ion channels, neurotransmitter transporters and synaptic components. Pathway analysis revealed nine significantly enriched gene sets, including regulation of insulin secretion and endocannabinoid signaling. Bipolar I disorder is strongly genetically correlated with schizophrenia, driven by psychosis, whereas bipolar II disorder is more strongly correlated with major depressive disorder. These findings address key clinical questions and provide potential biological mechanisms for bipolar disorder.

Efficient region-based test strategy uncovers genetic risk factors for functional outcome in bipolar disorder

Genome-wide association studies of case-control status have advanced the understanding of the genetic basis of psychiatric disorders. Further progress may be gained by increasing sample size but also by new analysis strategies that advance the exploitation of existing data, especially for clinically important quantitative phenotypes. The functionally-informed efficient region-based test strategy (FIERS) introduced herein uses prior knowledge on biological function and dependence of genotypes within a powerful statistical framework with improved sensitivity and specificity for detecting consistent genetic effects across studies. As proof of concept, FIERS was used for the first genome-wide single nucleotide polymorphism (SNP)-based investigation on bipolar disorder (BD) that focuses on an important aspect of disease course, the functional outcome. FIERS identified a significantly associated locus on chromosome 15 (hg38: chr15:48965004 - 49464789 bp) with consistent effect strength between two independent studies (GAIN/TGen: European Americans, BOMA: Germans; n = 1592 BD patients in total). Protective and risk haplotypes were found on the most strongly associated SNPs. They contain a CTCF binding site (rs586758); CTCF sites are known to regulate sets of genes within a chromatin domain. The rs586758 - rs2086256 - rs1904317 haplotype is located in the promoter flanking region of the COPS2 gene, close to microRNA4716, and the EID1, SHC4, DTWD1 genes as plausible biological candidates. While implication with BD is novel, COPS2, EID1, and SHC4 are known to be relevant for neuronal differentiation and function and DTWD1 for psychopharmacological side effects. The test strategy FIERS that enabled this discovery is equally applicable for tag SNPs and sequence data.

Detecting significant genotype-phenotype association rules in bipolar disorder: market research meets complex genetics

Background

Disentangling the etiology of common, complex diseases is a major challenge in genetic research. For bipolar disorder (BD), several genome-wide association studies (GWAS) have been performed. Similar to other complex disorders, major breakthroughs in explaining the high heritability of BD through GWAS have remained elusive. To overcome this dilemma, genetic research into BD, has embraced a variety of strategies such as the formation of large consortia to increase sample size and sequencing approaches. Here we advocate a complementary approach making use of already existing GWAS data: a novel data mining procedure to identify yet undetected genotype–phenotype relationships. We adapted association rule mining, a data mining technique traditionally used in retail market research, to identify frequent and characteristic genotype patterns showing strong associations to phenotype clusters. We applied this strategy to three independent GWAS datasets from 2835 phenotypically characterized patients with BD. In a discovery step, 20,882 candidate association rules were extracted.

Results

Two of these rules—one associated with eating disorder and the other with anxiety—remained significant in an independent dataset after robust correction for multiple testing. Both showed considerable effect sizes (odds ratio ~ 3.4 and 3.0, respectively) and support previously reported molecular biological findings.

Conclusion

Our approach detected novel specific genotype–phenotype relationships in BD that were missed by standard analyses like GWAS. While we developed and applied our method within the context of BD gene discovery, it may facilitate identifying highly specific genotype–phenotype relationships in subsets of genome-wide data sets of other complex phenotype with similar epidemiological properties and challenges to gene discovery efforts.

Electronic supplementary material

The online version of this article (10.1186/s40345-018-0132-x) contains supplementary material, which is available to authorized users.

Analysis of shared heritability in common disorders of the brain

Disorders of the brain can exhibit considerable epidemiological comorbidity and often share symptoms, provoking debate about their etiologic overlap. We quantified the genetic sharing of 25 brain disorders from genome-wide association studies of 265,218 patients and 784,643 control participants and assessed their relationship to 17 phenotypes from 1,191,588 individuals. Psychiatric disorders share common variant risk, whereas neurological disorders appear more distinct from one another and from the psychiatric disorders. We also identified significant sharing between disorders and a number of brain phenotypes, including cognitive measures. Further, we conducted simulations to explore how statistical power, diagnostic misclassification, and phenotypic heterogeneity affect genetic correlations. These results highlight the importance of common genetic variation as a risk factor for brain disorders and the value of heritability-based methods in understanding their etiology.

Using coordinate-based meta-analyses to explore structural imaging genetics

Imaging genetics has become a highly popular approach in the field of schizophrenia research. A frequently reported finding is that effects from common genetic variation are associated with a schizophrenia-related structural endophenotype. Genetic contributions to a structural endophenotype may be easier to delineate, when referring to biological rather than diagnostic criteria. We used coordinate-based meta-analyses, namely the anatomical likelihood estimation (ALE) algorithm on 30 schizophrenia-related imaging genetics studies, representing 44 single-nucleotide polymorphisms at 26 gene loci investigated in 4682 subjects. To test whether analyses based on biological information would improve the convergence of results, gene ontology (GO) terms were used to group the findings from the published studies. We did not find any significant results for the main contrast. However, our analysis enrolling studies on genotype × diagnosis interaction yielded two clusters in the left temporal lobe and the medial orbitofrontal cortex. All other subanalyses did not yield any significant results. To gain insight into possible biological relationships between the genes implicated by these clusters, we mapped five of them to GO terms of the category "biological process" (AKT1, CNNM2, DISC1, DTNBP1, VAV3), then five to "cellular component" terms (AKT1, CNNM2, DISC1, DTNBP1, VAV3), and three to "molecular function" terms (AKT1, VAV3, ZNF804A). A subsequent cluster analysis identified representative, non-redundant subsets of semantically similar terms that aided a further interpretation. We regard this approach as a new option to systematically explore the richness of the literature in imaging genetics.

Gene set enrichment analysis and expression pattern exploration implicate an involvement of neurodevelopmental processes in bipolar disorder

BACKGROUND

Bipolar disorder (BD) is a common and highly heritable disorder of mood. Genome-wide association studies (GWAS) have identified several independent susceptibility loci. In order to extract more biological information from GWAS data, multi-locus approaches represent powerful tools since they utilize knowledge about biological processes to integrate functional sets of genes at strongly to moderately associated loci.

METHODS

We conducted gene set enrichment analyses (GSEA) using 2.3 million single-nucleotide polymorphisms, 397 Reactome pathways and 24,025 patients with BD and controls. RNA expression of implicated individual genes and gene sets were examined in post-mortem brains across lifespan.

RESULTS

Two pathways showed a significant enrichment after correction for multiple comparisons in the GSEA: GRB2 events in ERBB2 signaling, for which 6 of 21 genes were BD associated (PFDR = 0.0377), and NCAM signaling for neurite out-growth, for which 11 out of 62 genes were BD associated (PFDR = 0.0451). Most pathway genes showed peaks of RNA co-expression during fetal development and infancy and mapped to neocortical areas and parts of the limbic system.

LIMITATIONS

Pathway associations were technically reproduced by two methods, although they were not formally replicated in independent samples. Gene expression was explored in controls but not in patients.

CONCLUSIONS

Pathway analysis in large GWAS data of BD and follow-up of gene expression patterns in healthy brains provide support for an involvement of neurodevelopmental processes in the etiology of this neuropsychiatric disease. Future studies are required to further evaluate the relevance of the implicated genes on pathway functioning and clinical aspects of BD.

Integration of transcriptomic and cytoarchitectonic data implicates a role for MAOA and TAC1 in the limbic-cortical network

Decoding the chain from genes to cognition requires detailed insights how areas with specific gene activities and microanatomical architectures contribute to brain function and dysfunction. The Allen Human Brain Atlas contains regional gene expression data, while the JuBrain Atlas offers three-dimensional cytoarchitectonic maps reflecting interindividual variability. To date, an integrated framework that combines the analytical benefits of both scientific platforms towards a multi-level brain atlas of adult humans was not available. We have, therefore, developed JuGEx, a new method for integrating tissue transcriptome and cytoarchitectonic segregation. We investigated differential gene expression in two JuBrain areas of the frontal pole that we have structurally and functionally characterized in previous studies. Our results show a significant upregulation of MAOA and TAC1 in the medial area frontopolaris which is a node in the limbic-cortical network and known to be susceptible for gray matter loss and behavioral dysfunction in patients with depression. The MAOA gene encodes an enzyme which is involved in the catabolism of dopamine, norepinephrine, serotonin, and other monoaminergic neurotransmitters. The TAC1 locus generates hormones that play a role in neuron excitations and behavioral responses. Overall, JuGEx provides a new tool for the scientific community that empowers research from basic, cognitive and clinical neuroscience in brain regions and disease models with regard to gene expression.

An Expanded Genome-Wide Association Study of Type 2 Diabetes in Europeans

To characterize type 2 diabetes (T2D)-associated variation across the allele frequency spectrum, we conducted a meta-analysis of genome-wide association data from 26,676 T2D case and 132,532 control subjects of European ancestry after imputation using the 1000 Genomes multiethnic reference panel. Promising association signals were followed up in additional data sets (of 14,545 or 7,397 T2D case and 38,994 or 71,604 control subjects). We identified 13 novel T2D-associated loci (P < 5 × 10^-8), including variants near the GLP2R, GIP, and HLA-DQA1 genes. Our analysis brought the total number of independent T2D associations to 128 distinct signals at 113 loci. Despite substantially increased sample size and more complete coverage of low-frequency variation, all novel associations were driven by common single nucleotide variants. Credible sets of potentially causal variants were generally larger than those based on imputation with earlier reference panels, consistent with resolution of causal signals to common risk haplotypes. Stratification of T2D-associated loci based on T2D-related quantitative trait associations revealed tissue-specific enrichment of regulatory annotations in pancreatic islet enhancers for loci influencing insulin secretion and in adipocytes, monocytes, and hepatocytes for insulin action-associated loci. These findings highlight the predominant role played by common variants of modest effect and the diversity of biological mechanisms influencing T2D pathophysiology.