Pathways-based analyses of whole-genome association study data in bipolar disorder reveal genes mediating ion channel activity and synaptic neurotransmission

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. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

VIZARD: analysis of Affymetrix Arabidopsis GeneChip data

SUMMARY: The Affymetrix GeneChip Arabidopsis genome array has proved to be a very powerful tool for the analysis of gene expression in Arabidopsis thaliana, the most commonly studied plant model organism. VIZARD is a Java program created at the University of California, Berkeley, to facilitate analysis of Arabidopsis GeneChip data. It includes several integrated tools for filtering, sorting, clustering and visualization of gene expression data as well as tools for the discovery of regulatory motifs in upstream sequences. VIZARD also includes annotation and upstream sequence databases for the majority of genes represented on the Affymetrix Arabidopsis GeneChip array. AVAILABILITY: VIZARD is available free of charge for educational, research, and not-for-profit purposes, and can be downloaded at http://www.anm.f2s.com/research/vizard/ CONTACT: moseyko@uclink4.berkeley.edu     

PLINK: a tool set for whole-genome association and population-based linkage analyses

Whole-genome association studies (WGAS) bring new computational, as well as analytic, challenges to researchers. Many existing genetic-analysis tools are not designed to handle such large data sets in a convenient manner and do not necessarily exploit the new opportunities that whole-genome data bring. To address these issues, we developed PLINK, an open-source C/C++ WGAS tool set. With PLINK, large data sets comprising hundreds of thousands of markers genotyped for thousands of individuals can be rapidly manipulated and analyzed in their entirety. As well as providing tools to make the basic analytic steps computationally efficient, PLINK also supports some novel approaches to whole-genome data that take advantage of whole-genome coverage. We introduce PLINK and describe the five main domains of function: data management, summary statistics, population stratification, association analysis, and identity-by-descent estimation. In particular, we focus on the estimation and use of identity-by-state and identity-by-descent information in the context of population-based whole-genome studies. This information can be used to detect and correct for population stratification and to identify extended chromosomal segments that are shared identical by descent between very distantly related individuals. Analysis of the patterns of segmental sharing has the potential to map disease loci that contain multiple rare variants in a population-based linkage analysis.     

CrossChip: a system supporting comparative analysis of different generations of Affymetrix arrays

SUMMARY: To increase compatibility between different generations of Affymetrix GeneChip arrays, we propose a method of filtering probes based on their sequences. Our method is implemented as a web-based service for downloading necessary materials for converting the raw data files (*.CEL) for comparative analysis. The user can specify the appropriate level of filtering by setting the criteria for the minimum overlap length between probe sequences and the minimum number of usable probe pairs per probe set. Our website supports a within-species comparison for human and mouse GeneChip arrays. AVAILABILITY: http://www.crosschip.org     

affy--analysis of Affymetrix GeneChip data at the probe level

MOTIVATION: The processing of the Affymetrix GeneChip data has been a recent focus for data analysts. Alternatives to the original procedure have been proposed and some of these new methods are widely used. RESULTS: The affy package is an R package of functions and classes for the analysis of oligonucleotide arrays manufactured by Affymetrix. The package is currently in its second release, affy provides the user with extreme flexibility when carrying out an analysis and make it possible to access and manipulate probe intensity data. In this paper, we present the main classes and functions in the package and demonstrate how they can be used to process probe-level data. We also demonstrate the importance of probe-level analysis when using the Affymetrix GeneChip platform.     

Molecular Pathways Involved in Prostate Carcinogenesis: Insights from Public Microarray Datasets

Background

Prostate cancer is currently the most frequently diagnosed malignancy in men and the second leading cause of cancer-related deaths in industrialized countries. Worldwide, an increase in prostate cancer incidence is expected due to an increased life-expectancy, aging of the population and improved diagnosis. Although the specific underlying mechanisms of prostate carcinogenesis remain unknown, prostate cancer is thought to result from a combination of genetic and environmental factors altering key cellular processes. To elucidate these complex interactions and to contribute to the understanding of prostate cancer progression and metastasis, analysis of large scale gene expression studies using bioinformatics approaches is used to decipher regulation of core processes.

Methodology/Principal Findings

In this study, a standardized quality control procedure and statistical analysis (http://www.arrayanalysis.org/) were applied to multiple prostate cancer datasets retrieved from the ArrayExpress data repository and pathway analysis using PathVisio (http://www.pathvisio.org/) was performed. The results led to the identification of three core biological processes that are strongly affected during prostate carcinogenesis: cholesterol biosynthesis, the process of epithelial-to-mesenchymal transition and an increased metabolic activity.

Conclusions

This study illustrates how a standardized bioinformatics evaluation of existing microarray data and subsequent pathway analysis can quickly and cost-effectively provide essential information about important molecular pathways and cellular processes involved in prostate cancer development and disease progression. The presented results may assist in biomarker profiling and the development of novel treatment approaches.     

Quantitative comparison of cDNA-AFLP, microarrays, and GeneChip expression data in Saccharomyces cerevisiae

cDNA-AFLP is a genome-wide expression analysis technology that does not require any prior knowledge of gene sequences. This PCR-based technique combines a high sensitivity with a high specificity, allowing detection of rarely expressed genes and distinguishing between homologous genes. In this report, we validated quantitative expression data of 110 cDNA-AFLP fragments in yeast with DNA microarrays and GeneChip data. The best correlation was found between cDNA-AFLP and GeneChip data. The cDNA-AFLP data revealed a low number of inconsistent profiles that could be explained by gel artifact, overexposure, or mismatch amplification. In addition, 18 cDNA-AFLP fragments displayed homology to genomic yeast DNA, but could not be linked unambiguously to any known ORF. These fragments were most probably derived from 5' or 3' noncoding sequences or might represent previously unidentified ORFs. Genes liable to cross hybridization showed identical results in cDNA-AFLP and GeneChip analysis. Three genes, which were readily detected with cDNA-AFLP, showed no significant expression in GeneChip experiments. We show that cDNA-AFLP is a very good alternative to microarrays and since no preexisting biological or sequence information is required, it is applicable to any species.     

Oligonucleotide arrays: information from replication and spatial structure

MOTIVATION: The introduction of oligonucleotide DNA arrays has resulted in much debate concerning appropriate models for the measurement of gene expression. By contrast, little account has been taken of the possibility of identifying the physical imperfections in the raw data. RESULTS: This paper demonstrates that, with the use of replicates and an awareness of the spatial structure, deficiencies in the data can be identified, the possibility of their correction can be ascertained and correction can be effected (by use of local scaling) where possible. The procedures were motivated by data from replicates of Arabidopsis thaliana using the GeneChip ATH1-121501 microarray. Similar problems are illustrated for GeneChip Human Genome U133 arrays and for the newer and larger GeneChip Wheat Genome microarray. AVAILABILITY: R code is freely available on request.     

Normality of oligonucleotide microarray data and implications for parametric statistical analyses

MOTIVATION: Experimental limitations have resulted in the popularity of parametric statistical tests as a method for identifying differentially regulated genes in microarray data sets. However, these tests assume that the data follow a normal distribution. To date, the assumption that replicate expression values for any gene are normally distributed, has not been critically addressed for Affymetrix GeneChip data. RESULTS: The normality of the expression values calculated using four different commercial and academic software packages was investigated using a data set consisting of the same target RNA applied to 59 human Affymetrix U95A GeneChips using a combination of statistical tests and visualization techniques. For the majority of probe sets obtained from each analysis suite, the expression data showed a good correlation with normality. The exception was a large number of low-expressed genes in the data set produced using Affymetrix Microarray Suite 5.0, which showed a striking non-normal distribution. In summary, our data provide strong support for the application of parametric tests to GeneChip data sets without the need for data transformation.