Despite known heritability, the complex genetic architecture of bipolar disorder (likely including trait, locus and allelic heterogeneity, as well as genetic interactions) has confounded genetic discovery for many years. Even modern day whole genome association studies (WGAS) using over half a million common SNPs have implicated only a handful of genes at the genomewide level. Temporally coincident with this series of WGAS, a host of pathways-based analyses (PBAs) have emerged as novel computational approaches in the examination of large-scale datasets, but thus far rarely have been applied to WGAS data in psychiatric disorders. Here, we report a series of PBAs conducted using exploratory visual analysis, an analytic and visualization software tool for examining genomic data, to examine results from the National Institutes of Mental Health and Wellcome-Trust Case Control Consortium WGAS in bipolar disorder. Consistent with a host of prior linkage findings, some candidate gene association studies, and recent WGAS, our strongest findings suggest involvement of ion channel structural and regulatory genes, including voltage-gated ion channels and the broader ion channel group that comprises both voltage- and ligand-gated channels. Moreover, we found only modest overlap in the particular genes driving the significance of these gene sets across the analyses. This observation strongly suggests that variation in ion channel genes, as a class of genes, may contribute to the susceptibility of bipolar disorder and that heterogeneity may figure prominently in the genetic architecture of this susceptibility.

Schizophrenia is a complex genetic disorder. Gene set-based analytic (GSA) methods have been widely applied for exploratory analyses of large, high-throughput datasets, but less commonly employed for biological hypothesis testing. Our primary hypothesis is that variation in ion channel genes contribute to the genetic susceptibility to schizophrenia. We applied Exploratory Visual Analysis (EVA), one GSA application, to analyze European-American (EA) and African-American (AA) schizophrenia genome-wide association study datasets for statistical enrichment of ion channel gene sets, comparing GSA results derived under three SNP-to-gene mapping strategies: (1) GENIC; (2) 500-Kb; (3) 2.5-Mb and three complimentary SNP-to-gene statistical reduction methods: (1) minimum p value (pMIN); (2) a novel method, proportion of SNPs per Gene with p-values below a pre-defined α-threshold (PROP); and (3) the truncated product method (TPM). In the EA analyses, ion channel gene set(s) were enriched under all mapping and statistical approaches. In the AA analysis, ion channel gene set(s) were significantly enriched under pMIN for all mapping strategies and under PROP for broader mapping strategies. Less extensive enrichment in the AA sample may reflect true ethnic differences in susceptibility, sampling or case ascertainment differences, or higher dimensionality relative to sample size of the AA data. More consistent findings under broader mapping strategies may reflect enhanced power due to increased SNP inclusion, enhanced capture of effects over extended haplotypes or significant contributions from regulatory regions. While extensive pMIN findings may reflect gene size bias, the extent and significance of PROP and TPM findings suggest that common variation at ion channel genes may capture some of the heritability of schizophrenia.