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.

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.

Haplotype synthesis analysis reveals functional variants underlying known genome-wide associated susceptibility loci.

MOTIVATION: The functional mechanisms underlying disease association remain unknown for Genome-wide Association Studies (GWAS) susceptibility variants located outside coding regions. Synthesis of effects from multiple surrounding functional variants has been suggested as an explanation of hard-to-interpret findings. We define filter criteria based on linkage disequilibrium measures and allele frequencies which reflect expected properties of synthesizing variant sets. For eligible candidate sets, we search for haplotype markers that are highly correlated with associated variants. RESULTS: Via simulations we assess the performance of our approach and suggest parameter settings which guarantee 95% sensitivity at 20-fold reduced computational cost. We apply our method to 1000 Genomes data and confirmed Crohn's Disease (CD) and Type 2 Diabetes (T2D) variants. A proportion of 36.9% allowed explanation by three-variant-haplotypes carrying at least two functional variants, as compared to 16.4% for random variants ([Formula: see text]). Association could be explained by missense variants for MUC19, PER3 (CD) and HMG20A (T2D). In a CD GWAS-imputed using haplotype reference consortium data (64 976 haplotypes)-we could confirm the syntheses of MUC19 and PER3 and identified synthesis by missense variants for 6 further genes (ZGPAZ, GPR65, CLN3/NPIPB8, LOC102723878, rs2872507, GCKR). In all instances, the odds ratios of the synthesizing haplotypes were virtually identical to that of the index SNP. In summary, we demonstrate the potential of synthesis analysis to guide functional follow-up of GWAS findings. AVAILABILITY AND IMPLEMENTATION: All methods are implemented in the C/C ++ toolkit GetSynth, available at http://sourceforge.net/projects/getsynth/ CONTACT: tim.becker@uni-greifswald.de SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Alzheimer's disease risk variants modulate endophenotypes in mild cognitive impairment.

INTRODUCTION: We evaluated the effect of Alzheimer's disease (AD) susceptibility loci on endophenotypes closely related with AD pathology in patients with mild cognitive impairment (MCI). METHODS: We selected 1730 MCI patients from four independent data sets. Weighted polygenic risk scores (PGS) were constructed of 18 non-apolipoprotein E (APOE) AD risk variants. In addition, we determined APOE genotype. AD endophenotypes were cognitive decline over time and cerebrospinal fluid (CSF) biomarkers (aß, tau, ptau). RESULTS: PGS was modestly associated with cognitive decline over time, as measured by mini-mental state examination (MMSE) (ß ± SE:-0.24 ± 0.10; P = .012), and with CSF levels of tau and ptau (tau: 1.38 ± 0.36, P = 1.21 × 10(-4); ptau: 1.40 ± 0.36, P = 1.02 × 10(-4)). DISCUSSION: In MCI, we observed a joint effect of AD susceptibility loci on nonamyloid endophenotypes, suggesting a link of these genetic loci with neuronal degeneration in general rather than with Alzheimer-related amyloid deposition.

Novel genetic matching methods for handling population stratification in genome-wide association studies.

BACKGROUND: A usually confronted problem in association studies is the occurrence of population stratification. In this work, we propose a novel framework to consider population matchings in the contexts of genome-wide and sequencing association studies. We employ pairwise and groupwise optimal case-control matchings and present an agglomerative hierarchical clustering, both based on a genetic similarity score matrix. In order to ensure that the resulting matches obtained from the matching algorithm capture correctly the population structure, we propose and discuss two stratum validation methods. We also invent a decisive extension to the Cochran-Armitage Trend test to explicitly take into account the particular population structure. RESULTS: We assess our framework by simulations of genotype data under the null hypothesis, to affirm that it correctly controls for the type-1 error rate. By a power study we evaluate that structured association testing using our framework displays reasonable power. We compare our result with those obtained from a logistic regression model with principal component covariates. Using the principal components approaches we also find a possible false-positive association to Alzheimer's disease, which is neither supported by our new methods, nor by the results of a most recent large meta analysis or by a mixed model approach. CONCLUSIONS: Matching methods provide an alternative handling of confounding due to population stratification for statistical tests for which covariates are hard to model. As a benchmark, we show that our matching framework performs equally well to state of the art models on common variants.

METAINTER: meta-analysis of multiple regression models in genome-wide association studies.

MOTIVATION: Meta-analysis of summary statistics is an essential approach to guarantee the success of genome-wide association studies (GWAS). Application of the fixed or random effects model to single-marker association tests is a standard practice. More complex methods of meta-analysis involving multiple parameters have not been used frequently, a gap that could be explained by the lack of a respective meta-analysis pipeline. Meta-analysis based on combining p-values can be applied to any association test. However, to be powerful, meta-analysis methods for high-dimensional models should incorporate additional information such as study-specific properties of parameter estimates, their effect directions, standard errors and covariance structure. RESULTS: We modified 'method for the synthesis of linear regression slopes' recently proposed in the educational sciences to the case of multiple logistic regression, and implemented it in a meta-analysis tool called METAINTER. The software handles models with an arbitrary number of parameters, and can directly be applied to analyze the results of single-SNP tests, global haplotype tests, tests for and under gene-gene or gene-environment interaction. Via simulations for two-single nucleotide polymorphisms (SNP) models we have shown that the proposed meta-analysis method has correct type I error rate. Moreover, power estimates come close to that of the joint analysis of the entire sample. We conducted a real data analysis of six GWAS of type 2 diabetes, available from dbGaP (http://www.ncbi.nlm.nih.gov/gap). For each study, a genome-wide interaction analysis of all SNP pairs was performed by logistic regression tests. The results were then meta-analyzed with METAINTER. AVAILABILITY: The software is freely available and distributed under the conditions specified on http://metainter.meb.uni-bonn.de. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Analysis of the progression of systolic blood pressure using imputation of missing phenotype values.

We present a genome-wide association study of a quantitative trait, "progression of systolic blood pressure in time," in which 142 unrelated individuals of the Genetic Analysis Workshop 18 real genotype data were analyzed. Information on systolic blood pressure and other phenotypic covariates was missing at certain time points for a considerable part of the sample. We observed that the dropout process causing missingness is not independent of the initial systolic blood pressure; that is, the data is not missing completely at random. However, after the adjustment for age, the impact of systolic blood pressure on dropouts was no longer significant. Therefore, we decided to impute missing phenotype values by using information from individuals with complete phenotypic data. Progression of systolic blood pressure (∆SBP/∆t) was defined based on the imputed phenotypes and analyzed in a genome-wide fashion. We also conducted an exhaustive genome-wide search for interaction between single-nucleotide polymorphisms (7.14 × 10(10) tests) under an allelic model. The suggested data imputation and the association analysis strategy proved to be valid in the sense that there was no evidence of genome-wide inflation or increased type I error in general. Furthermore, we detected 2 single-nucleotide polymorphisms (SNPs) that met the criterion for genome-wide significance (p≤5 × 10(-8)), which was also confirmed via Monte-Carlo simulation. In view of the rather small sample size, however, the results have to be followed-up in larger studies.

Rare variant testing of imputed data: an analysis pipeline typified.

Important methodological advancements in rare variant association testing have been made recently, among them collapsing tests, kernel methods and the variable threshold (VT) technique. Typically, rare variants from a region of interest are tested for association as a group ('bin'). Rare variant studies are already routinely performed as whole-exome sequencing studies. As an alternative approach, we propose a pipeline for rare variant analysis of imputed data and develop respective quality control criteria. We provide suggestions for the choice and construction of analysis bins in whole-genome application and support the analysis with implementations of standard burden tests (COLL, CMAT) in our INTERSNP-RARE software. In addition, three rare variant regression tests (REG, FRACREG and COLLREG) are implemented. All tests are accompanied with the VT approach which optimizes the definition of 'rareness'. We integrate kernel tests as implemented in SKAT/SKAT-O into the suggested strategies. Then, we apply our analysis scheme to a genome-wide association study of Alzheimer's disease. Further, we show that our pipeline leads to valid significance testing procedures with controlled type I error rates. Strong association signals surrounding the known APOE locus demonstrate statistical power. In addition, we highlight several suggestive rare variant association findings for follow-up studies, including genomic regions overlapping MCPH1, MED18 and NOTCH3. In summary, we describe and support a straightforward and cost-efficient rare variant analysis pipeline for imputed data and demonstrate its feasibility and validity. The strategy can complement rare variant studies with next generation sequencing data.

SUCLG2 identified as both a determinator of CSF Aß1-42 levels and an attenuator of cognitive decline in Alzheimer's disease.

Cerebrospinal fluid amyloid-beta 1-42 (Aß1-42) and phosphorylated Tau at position 181 (pTau181) are biomarkers of Alzheimer's disease (AD). We performed an analysis and meta-analysis of genome-wide association study data on Aß1-42 and pTau181 in AD dementia patients followed by independent replication. An association was found between Aß1-42 level and a single-nucleotide polymorphism in SUCLG2 (rs62256378) (P = 2.5×10(-12)). An interaction between APOE genotype and rs62256378 was detected (P = 9.5 × 10(-5)), with the strongest effect being observed in APOE-ε4 noncarriers. Clinically, rs62256378 was associated with rate of cognitive decline in AD dementia patients (P = 3.1 × 10(-3)). Functional microglia experiments showed that SUCLG2 was involved in clearance of Aß1-42.

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.

A one-degree-of-freedom test for supra-multiplicativity of SNP effects.

Deviation from multiplicativity of genetic risk factors is biologically plausible and might explain why Genome-wide association studies (GWAS) so far could unravel only a portion of disease heritability. Still, evidence for SNP-SNP epistasis has rarely been reported, suggesting that 2-SNP models are overly simplistic. In this context, it was recently proposed that the genetic architecture of complex diseases could follow limiting pathway models. These models are defined by a critical risk allele load and imply multiple high-dimensional interactions. Here, we present a computationally efficient one-degree-of-freedom "supra-multiplicativity-test" (SMT) for SNP sets of size 2 to 500 that is designed to detect risk alleles whose joint effect is fortified when they occur together in the same individual. Via a simulation study we show that the SMT is powerful in the presence of threshold models, even when only about 30-45% of the model SNPs are available. In addition, we demonstrate that the SMT outperforms standard interaction analysis under recessive models involving just a few SNPs. We apply our test to 10 consensus Alzheimer's disease (AD) susceptibility SNPs that were previously identified by GWAS and obtain evidence for supra-multiplicativity ([Formula: see text]) that is not attributable to either two-way or three-way interaction.

Quick, "imputation-free" meta-analysis with proxy-SNPs.

BACKGROUND: Meta-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. RESULTS: Here 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. CONCLUSIONS: YAMAS 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.

Integrated genome-wide pathway association analysis with INTERSNP.

OBJECTIVES: Pathway association analysis (PAA) tests for an excess of moderately significant SNPs in genes from a common pathway. METHODS: We present a Monte-Carlo simulation framework that allows to formulate the main ideas of existing PAA approaches using a self-contained rather than a competitive null hypothesis. A stand-alone implementation in INTERSNP makes time-consuming communication with standard GWAS software redundant. By additional parallelization with the OpenMP API, we achieve a reduction in running time for PAA by orders of magnitude, making a power simulation study for PAA feasible. Our approach properly accounts for linkage disequilibrium and is robust with respect to residual λ inflation. RESULTS: We demonstrate that under simple, realistic disease models, PAA can actually strongly outperform the GWAS single-marker approach. PAA methods that make use of the strength of the SNP association (GenGen, Fisher's combination test), in general, perform better than ratio-based methods (ALIGATOR, SNP ratio), whereas the relative performance of gene-based scoring (ALIGATOR, GenGen) and pathway-based scoring (SNP ratio, Fisher's combination test) depends on the architecture of the assumed disease model. Finally, we present a new PAA score that models independent signals from the same gene in a regression framework and show that it is a reasonable compromise that combines the advantages of existing ideas.