Markus Leber

Genome-wide association study reveals two new risk loci for bipolar disorder

Bipolar disorder (BD) is a common and highly heritable mental illness and genome-wide association studies (GWAS) have robustly identified the first common genetic variants involved in disease aetiology. The data also provide strong evidence for the presence of multiple additional risk loci, each contributing a relatively small effect to BD susceptibility. Large samples are necessary to detect these risk loci. Here we present results from the largest BD GWAS to date by investigating 2.3 million single-nucleotide polymorphisms (SNPs) in a sample of 24,025 patients and controls. We detect 56 genome-wide significant SNPs in five chromosomal regions including previously reported risk loci ANK3, ODZ4 and TRANK1, as well as the risk locus ADCY2 (5p15.31) and a region between MIR2113 and POU3F2 (6q16.1). ADCY2 is a key enzyme in cAMP signalling and our finding provides new insights into the biological mechanisms involved in the development of BD.

Genome-wide association analysis of genetic generalized epilepsies implicates susceptibility loci at 1q43, 2p16.1, 2q22.3 and 17q21.32

Genetic generalized epilepsies (GGEs) have a lifetime prevalence of 0.3% and account for 20-30% of all epilepsies. Despite their high heritability of 80%, the genetic factors predisposing to GGEs remain elusive. To identify susceptibility variants shared across common GGE syndromes, we carried out a two-stage genome-wide association study (GWAS) including 3020 patients with GGEs and 3954 controls of European ancestry. To dissect out syndrome-related variants, we also explored two distinct GGE subgroups comprising 1434 patients with genetic absence epilepsies (GAEs) and 1134 patients with juvenile myoclonic epilepsy (JME). Joint Stage-1 and 2 analyses revealed genome-wide significant associations for GGEs at 2p16.1 (rs13026414, P(meta) = 2.5 × 10(-9), OR[T] = 0.81) and 17q21.32 (rs72823592, P(meta) = 9.3 × 10(-9), OR[A] = 0.77). The search for syndrome-related susceptibility alleles identified significant associations for GAEs at 2q22.3 (rs10496964, P(meta) = 9.1 × 10(-9), OR[T] = 0.68) and at 1q43 for JME (rs12059546, P(meta) = 4.1 × 10(-8), OR[G] = 1.42). Suggestive evidence for an association with GGEs was found in the region 2q24.3 (rs11890028, P(meta) = 4.0 × 10(-6)) nearby the SCN1A gene, which is currently the gene with the largest number of known epilepsy-related mutations. The associated regions harbor high-ranking candidate genes: CHRM3 at 1q43, VRK2 at 2p16.1, ZEB2 at 2q22.3, SCN1A at 2q24.3 and PNPO at 17q21.32. Further replication efforts are necessary to elucidate whether these positional candidate genes contribute to the heritability of the common GGE syndromes.

Follow-up of loci from the International Genomics of Alzheimer's Disease Project identifies TRIP4 as a novel susceptibility gene

To follow-up loci discovered by the International Genomics of Alzheimer’s Disease Project, we attempted independent replication of 19 single nucleotide polymorphisms (SNPs) in a large Spanish sample (Fundació ACE data set; 1808 patients and 2564 controls). Our results corroborate association with four SNPs located in the genes INPP5D, MEF2C, ZCWPW1 and FERMT2, respectively. Of these, ZCWPW1 was the only SNP to withstand correction for multiple testing (P=0.000655). Furthermore, we identify TRIP4 (rs74615166) as a novel genome-wide significant locus for Alzheimer’s disease risk (odds ratio=1.31; confidence interval 95% (1.19–1.44); P=9.74 × 10−9).

Genome-wide analysis implicates microRNAs and their target genes in the development of bipolar disorder

Bipolar disorder (BD) is a severe and highly heritable neuropsychiatric disorder with a lifetime prevalence of 1%. Molecular genetic studies have identified the first BD susceptibility genes. However, the disease pathways remain largely unknown. Accumulating evidence suggests that microRNAs, a class of small noncoding RNAs, contribute to basic mechanisms underlying brain development and plasticity, suggesting their possible involvement in the pathogenesis of several psychiatric disorders, including BD. In the present study, gene-based analyses were performed for all known autosomal microRNAs using the largest genome-wide association data set of BD to date (9747 patients and 14 278 controls). Associated and brain-expressed microRNAs were then investigated in target gene and pathway analyses. Functional analyses of miR-499 and miR-708 were performed in rat hippocampal neurons. Ninety-eight of the six hundred nine investigated microRNAs showed nominally significant P-values, suggesting that BD-associated microRNAs might be enriched within known microRNA loci. After correction for multiple testing, nine microRNAs showed a significant association with BD. The most promising were miR-499, miR-708 and miR-1908. Target gene and pathway analyses revealed 18 significant canonical pathways, including brain development and neuron projection. For miR-499, four Bonferroni-corrected significant target genes were identified, including the genome-wide risk gene for psychiatric disorder CACNB2. First results of functional analyses in rat hippocampal neurons neither revealed nor excluded a major contribution of miR-499 or miR-708 to dendritic spine morphogenesis. The present results suggest that research is warranted to elucidate the precise involvement of microRNAs and their downstream pathways in BD.

Integrated Pathway-Based Approach Identifies Association between Genomic Regions at CTCF and CACNB2 and Schizophrenia

In the present study, an integrated hierarchical approach was applied to: (1) identify pathways associated with susceptibility to schizophrenia; (2) detect genes that may be potentially affected in these pathways since they contain an associated polymorphism; and (3) annotate the functional consequences of such single-nucleotide polymorphisms (SNPs) in the affected genes or their regulatory regions. The Global Test was applied to detect schizophrenia-associated pathways using discovery and replication datasets comprising 5,040 and 5,082 individuals of European ancestry, respectively. Information concerning functional gene-sets was retrieved from the Kyoto Encyclopedia of Genes and Genomes, Gene Ontology, and the Molecular Signatures Database. Fourteen of the gene-sets or pathways identified in the discovery dataset were confirmed in the replication dataset. These include functional processes involved in transcriptional regulation and gene expression, synapse organization, cell adhesion, and apoptosis. For two genes, i.e. CTCF and CACNB2, evidence for association with schizophrenia was available (at the gene-level) in both the discovery study and published data from the Psychiatric Genomics Consortium schizophrenia study. Furthermore, these genes mapped to four of the 14 presently identified pathways. Several of the SNPs assigned to CTCF and CACNB2 have potential functional consequences, and a gene in close proximity to CACNB2, i.e. ARL5B, was identified as a potential gene of interest. Application of the present hierarchical approach thus allowed: (1) identification of novel biological gene-sets or pathways with potential involvement in the etiology of schizophrenia, as well as replication of these findings in an independent cohort; (2) detection of genes of interest for future follow-up studies; and (3) the highlighting of novel genes in previously reported candidate regions for schizophrenia.

A fractional programming approach to efficient DNA melting temperature calculation

MOTIVATION In a wide range of experimental techniques in biology, there is a need for an efficient method to calculate the melting temperature of pairings of two single DNA strands. Avoiding cross-hybridization when choosing primers for the polymerase chain reaction or selecting probes for large-scale DNA assays are examples where the exact determination of melting temperatures is important. Beyond being exact, the method has to be efficient, as these techniques often require the simultaneous calculation of melting temperatures of up to millions of possible pairings. The problem is to simultaneously determine the most stable alignment of two sequences, including potential loops and bulges, and calculate the corresponding melting temperature. RESULTS As the melting temperature can be expressed as a fraction in terms of enthalpy and entropy differences of the corresponding annealing reaction, we propose to use a fractional programming algorithm, the Dinkelbach algorithm, to solve the problem. To calculate the required differences of enthalpy and entropy, the Nearest Neighbor model is applied. Using this model, the substeps of the Dinkelbach algorithm in our problem setting turn out to be calculations of alignments which optimize an additive score function. Thus, the usual dynamic programming techniques can be applied. The result is an efficient algorithm to determine melting temperatures of two DNA strands, suitable for large-scale applications such as primer or probe design. AVAILABILITY The software is available for academic purposes from the authors. A web interface is provided at http://www.zaik.uni-koeln.de/bioinformatik/fptm.html

Identification of shared risk loci and pathways for bipolar disorder and schizophrenia

Bipolar disorder (BD) is a highly heritable neuropsychiatric disease characterized by recurrent episodes of mania and depression. BD shows substantial clinical and genetic overlap with other psychiatric disorders, in particular schizophrenia (SCZ). The genes underlying this etiological overlap remain largely unknown. A recent SCZ genome wide association study (GWAS) by the Psychiatric Genomics Consortium identified 128 independent genome-wide significant single nucleotide polymorphisms (SNPs). The present study investigated whether these SCZ-associated SNPs also contribute to BD development through the performance of association testing in a large BD GWAS dataset (9747 patients, 14278 controls). After re-imputation and correction for sample overlap, 22 of 107 investigated SCZ SNPs showed nominal association with BD. The number of shared SCZ-BD SNPs was significantly higher than expected (p = 1.46x10-8). This provides further evidence that SCZ-associated loci contribute to the development of BD. Two SNPs remained significant after Bonferroni correction. The most strongly associated SNP was located near TRANK1, which is a reported genome-wide significant risk gene for BD. Pathway analyses for all shared SCZ-BD SNPs revealed 25 nominally enriched gene-sets, which showed partial overlap in terms of the underlying genes. The enriched gene-sets included calcium- and glutamate signaling, neuropathic pain signaling in dorsal horn neurons, and calmodulin binding. The present data provide further insights into shared risk loci and disease-associated pathways for BD and SCZ. This may suggest new research directions for the treatment and prevention of these two major psychiatric disorders.

The protocadherin 17 gene affects cognition, personality, amygdala structure and function, synapse development and risk of major mood disorders

Major mood disorders, which primarily include bipolar disorder and major depressive disorder, are the leading cause of disability worldwide and pose a major challenge in identifying robust risk genes. Here, we present data from independent large-scale clinical data sets (including 29 557 cases and 32 056 controls) revealing brain expressed protocadherin 17 (PCDH17) as a susceptibility gene for major mood disorders. Single-nucleotide polymorphisms (SNPs) spanning the PCDH17 region are significantly associated with major mood disorders; subjects carrying the risk allele showed impaired cognitive abilities, increased vulnerable personality features, decreased amygdala volume and altered amygdala function as compared with non-carriers. The risk allele predicted higher transcriptional levels of PCDH17 mRNA in postmortem brain samples, which is consistent with increased gene expression in patients with bipolar disorder compared with healthy subjects. Further, overexpression of PCDH17 in primary cortical neurons revealed significantly decreased spine density and abnormal dendritic morphology compared with control groups, which again is consistent with the clinical observations of reduced numbers of dendritic spines in the brains of patients with major mood disorders. Given that synaptic spines are dynamic structures which regulate neuronal plasticity and have crucial roles in myriad brain functions, this study reveals a potential underlying biological mechanism of a novel risk gene for major mood disorders involved in synaptic function and related intermediate phenotypes.

Quick, “Imputation-free” meta-analysis with proxy-SNPs

BackgroundMeta-analysis (MA) is widely used to pool genome-wide association studies (GWASes) in order to a) increase the power to detect strong or weak genotype effects or b) as a result verification method. As a consequence of differing SNP panels among genotyping chips, imputation is the method of choice within GWAS consortia to avoid losing too many SNPs in a MA. YAMAS (Yet Another Meta Analysis Software), however, enables cross-GWAS conclusions prior to finished and polished imputation runs, which eventually are time-consuming.ResultsHere we present a fast method to avoid forfeiting SNPs present in only a subset of studies, without relying on imputation. This is accomplished by using reference linkage disequilibrium data from 1,000 Genomes/HapMap projects to find proxy-SNPs together with in-phase alleles for SNPs missing in at least one study. MA is conducted by combining association effect estimates of a SNP and those of its proxy-SNPs. Our algorithm is implemented in the MA software YAMAS. Association results from GWAS analysis applications can be used as input files for MA, tremendously speeding up MA compared to the conventional imputation approach. We show that our proxy algorithm is well-powered and yields valuable ad hoc results, possibly providing an incentive for follow-up studies. We propose our method as a quick screening step prior to imputation-based MA, as well as an additional main approach for studies without available reference data matching the ethnicities of study participants. As a proof of principle, we analyzed six dbGaP Type II Diabetes GWAS and found that the proxy algorithm clearly outperforms naïve MA on the p-value level: for 17 out of 23 we observe an improvement on the p-value level by a factor of more than two, and a maximum improvement by a factor of 2127.ConclusionsYAMAS is an efficient and fast meta-analysis program which offers various methods, including conventional MA as well as inserting proxy-SNPs for missing markers to avoid unnecessary power loss. MA with YAMAS can be readily conducted as YAMAS provides a generic parser for heterogeneous tabulated file formats within the GWAS field and avoids cumbersome setups. In this way, it supplements the meta-analysis process.

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.