Incorporating regulatory interactions into gene-set analyses for GWAS data: A controlled analysis with the MAGMA tool

To date, genome-wide association studies have identified thousands of statistically-significant associations between genetic variants, and phenotypes related to a myriad of traits and diseases. A key goal for human-genetics research is to translate these associations into functional mechanisms. Popular gene-set analysis tools, like MAGMA, map variants to genes they might affect, and then integrate genome-wide association study data (that is, variant-level associations for a phenotype) to score genes for association with a phenotype. Gene scores are subsequently used in competitive gene-set analyses to identify biological processes that are enriched for phenotype association. By default, variants are mapped to genes in their proximity. However, many variants that affect phenotypes are thought to act at regulatory elements, which can be hundreds of kilobases away from their target genes. Thus, we explored the idea of augmenting a proximity-based mapping scheme with publicly-available datasets of regulatory interactions. We used MAGMA to analyze genome-wide association study data for ten different phenotypes, and evaluated the effects of augmentation by comparing numbers, and identities, of genes and gene sets detected as statistically significant between mappings. We detected several pitfalls and confounders of such “augmented analyses”, and introduced ways to control for them. Using these controls, we demonstrated that augmentation with datasets of regulatory interactions only occasionally strengthened the enrichment for phenotype association amongst (biologically-relevant) gene sets for different phenotypes. Still, in such cases, genes and regulatory elements responsible for the improvement could be pinpointed. For instance, using brain regulatory-interactions for augmentation, we were able to implicate two acetylcholine receptor subunits involved in post-synaptic chemical transmission, namely CHRNB2 and CHRNE, in schizophrenia. Collectively, our study presents a critical approach for integrating regulatory interactions into gene-set analyses for genome-wide association study data, by introducing various controls to distinguish genuine results from spurious discoveries.     

FLAGS: A Flexible and Adaptive Association Test for Gene Sets Using Summary Statistics

Genome-wide association studies (GWAS) have been widely used for identifying common variants associated with complex diseases. Despite remarkable success in uncovering many risk variants and providing novel insights into disease biology, genetic variants identified to date fail to explain the vast majority of the heritability for most complex diseases. One explanation is that there are still a large number of common variants that remain to be discovered, but their effect sizes are generally too small to be detected individually. Accordingly, gene set analysis of GWAS, which examines a group of functionally related genes, has been proposed as a complementary approach to single-marker analysis. Here, we propose a flexible and adaptive test for gene sets (FLAGS), using summary statistics. Extensive simulations showed that this method has an appropriate type I error rate and outperforms existing methods with increased power. As a proof of principle, through real data analyses of Crohn’s disease GWAS data and bipolar disorder GWAS meta-analysis results, we demonstrated the superior performance of FLAGS over several state-of-the-art association tests for gene sets. Our method allows for the more powerful application of gene set analysis to complex diseases, which will have broad use given that GWAS summary results are increasingly publicly available.     

GeneTools – application for functional annotation and statistical hypothesis testing

Background  

Modern biology has shifted from "one gene" approaches to methods for genomic-scale analysis like microarray technology, which allow simultaneous measurement of thousands of genes. This has created a need for tools facilitating interpretation of biological data in "batch" mode. However, such tools often leave the investigator with large volumes of apparently unorganized information. To meet this interpretation challenge, gene-set, or cluster testing has become a popular analytical tool. Many gene-set testing methods and software packages are now available, most of which use a variety of statistical tests to assess the genes in a set for biological information. However, the field is still evolving, and there is a great need for "integrated" solutions.     

Association of NOD2 and IL23R with Inflammatory Bowel Disease in Puerto Rico

The Puerto Rico population may be modeled as an admixed population with contributions from three continents: Sub-Saharan Africa, Ancient America, and Europe. Extending the study of the genetics of inflammatory bowel disease (IBD) to an admixed population such as Puerto Rico has the potential to shed light on IBD genes identified in studies of European populations, find new genes contributing to IBD susceptibility, and provide basic information on IBD for the care of US patients of Puerto Rican and Latino descent. In order to study the association between immune-related genes and Crohn’s disease (CD) and ulcerative colitis (UC) in Puerto Rico, we genotyped 1159 Puerto Rican cases, controls, and family members with the ImmunoChip. We also genotyped 832 subjects from the Human Genome Diversity Panel to provide data for estimation of global and local continental ancestry. Association of SNPs was tested by logistic regression corrected for global continental descent and family structure. We observed the association between Crohn’s disease and NOD2 (rs17313265, 0.28 in CD, 0.19 in controls, OR 1.5, p = 9×10⁻⁶) and IL23R (rs11209026, 0.026 in CD, 0.0.071 in controls, OR 0.4, p = 3.8×10⁻⁴). The haplotype structure of both regions resembled that reported for European populations and “local” continental ancestry of the IL23R gene was almost entirely of European descent. We also observed suggestive evidence for the association of the BAZ1A promoter SNP with CD (rs1200332, 0.45 in CD, 0.35 in controls, OR 1.5, p = 2×10⁻⁶). Our estimate of continental ancestry surrounding this SNP suggested an origin in Ancient America for this putative susceptibility region. Our observations underscored the great difference between global continental ancestry and local continental ancestry at the level of the individual gene, particularly for immune-related loci.     

Categorizer: a tool to categorize genes into user-defined biological groups based on semantic similarity

Background

Communalities between large sets of genes obtained from high-throughput experiments are often identified by searching for enrichments of genes with the same Gene Ontology (GO) annotations. The GO analysis tools used for these enrichment analyses assume that GO terms are independent and the semantic distances between all parent–child terms are identical, which is not true in a biological sense. In addition these tools output lists of often redundant or too specific GO terms, which are difficult to interpret in the context of the biological question investigated by the user. Therefore, there is a demand for a robust and reliable method for gene categorization and enrichment analysis.

Results

We have developed Categorizer, a tool that classifies genes into user-defined groups (categories) and calculates p-values for the enrichment of the categories. Categorizer identifies the biologically best-fit category for each gene by taking advantage of a specialized semantic similarity measure for GO terms. We demonstrate that Categorizer provides improved categorization and enrichment results of genetic modifiers of Huntington’s disease compared to a classical GO Slim-based approach or categorizations using other semantic similarity measures.

Conclusion

Categorizer enables more accurate categorizations of genes than currently available methods. This new tool will help experimental and computational biologists analyzing genomic and proteomic data according to their specific needs in a more reliable manner.     

A genome-wide association study identified new variants associated with mathematical abilities in Chinese children

Mathematical ability is moderately heritable, and it is a complex trait which can be evaluated in several different categories. A few genetic studies have been published on general mathematical ability. However, no genetic study focused on specific mathematical ability categories. In this study, we separately performed genome-wide association studies on 11 mathematical ability categories in 1146 students from Chinese elementary schools. We identified seven genome-wide significant single nucleotide polymorphisms (SNPs) with strong linkage disequilibrium among each other (all r² > 0.8) associated with mathematical reasoning ability (top SNP: rs34034296, p = 2.01 × 10⁻⁸, nearest gene: CUB and Sushi multiple domains 3, CSMD3). We replicated one SNP (rs133885) from 585 SNPs previously reported to be associated with general mathematical ability associated with division ability in our data (p = 1.053 × 10⁻⁵). In the gene- and gene-set enrichment analysis by MAGMA, we found three significant enrichments of associations with three mathematical ability categories for three genes (LINGO2, OAS1 and HECTD1). We also observed four significant enrichments of associations with four mathematical ability categories for three gene sets. Our results suggest new candidate genetic loci for the genetics of mathematical ability.     

Investigating the effect of paralogs on microarray gene-set analysis

Background  

In order to interpret the results obtained from a microarray experiment, researchers often shift focus from analysis of individual differentially expressed genes to analyses of sets of genes. These gene-set analysis (GSA) methods use previously accumulated biological knowledge to group genes into sets and then aim to rank these gene sets in a way that reflects their relative importance in the experimental situation in question. We suspect that the presence of paralogs affects the ability of GSA methods to accurately identify the most important sets of genes for subsequent research.     

The Central Role of AMP-Kinase and Energy Homeostasis Impairment in Alzheimer’s Disease: A Multifactor Network Analysis

Alzheimer’s disease is the most common cause of dementia worldwide, affecting the elderly population. It is characterized by the hallmark pathology of amyloid-β deposition, neurofibrillary tangle formation, and extensive neuronal degeneration in the brain. Wealth of data related to Alzheimer’s disease has been generated to date, nevertheless, the molecular mechanism underlying the etiology and pathophysiology of the disease is still unknown. Here we described a method for the combined analysis of multiple types of genome-wide data aimed at revealing convergent evidence interest that would not be captured by a standard molecular approach. Lists of Alzheimer-related genes (seed genes) were obtained from different sets of data on gene expression, SNPs, and molecular targets of drugs. Network analysis was applied for identifying the regions of the human protein-protein interaction network showing a significant enrichment in seed genes, and ultimately, in genes associated to Alzheimer’s disease, due to the cumulative effect of different combinations of the starting data sets. The functional properties of these enriched modules were characterized, effectively considering the role of both Alzheimer-related seed genes and genes that closely interact with them. This approach allowed us to present evidence in favor of one of the competing theories about AD underlying processes, specifically evidence supporting a predominant role of metabolism-associated biological process terms, including autophagy, insulin and fatty acid metabolic processes in Alzheimer, with a focus on AMP-activated protein kinase. This central regulator of cellular energy homeostasis regulates a series of brain functions altered in Alzheimer’s disease and could link genetic perturbation with neuronal transmission and energy regulation, representing a potential candidate to be targeted by therapy.     

Functional-Network-Based Gene Set Analysis Using Gene-Ontology

To account for the functional non-equivalence among a set of genes within a biological pathway when performing gene set analysis, we introduce GOGANPA, a network-based gene set analysis method, which up-weights genes with functions relevant to the gene set of interest. The genes are weighted according to its degree within a genome-scale functional network constructed using the functional annotations available from the gene ontology database. By benchmarking GOGANPA using a well-studied P53 data set and three breast cancer data sets, we will demonstrate the power and reproducibility of our proposed method over traditional unweighted approaches and a competing network-based approach that involves a complex integrated network. GOGANPA’s sole reliance on gene ontology further allows GOGANPA to be widely applicable to the analysis of any gene-ontology-annotated genome.     

VennPainter: A Tool for the Comparison and Identification of Candidate Genes Based on Venn Diagrams

VennPainter is a program for depicting unique and shared sets of genes lists and generating Venn diagrams, by using the Qt C++ framework. The software produces Classic Venn, Edwards’ Venn and Nested Venn diagrams and allows for eight sets in a graph mode and 31 sets in data processing mode only. In comparison, previous programs produce Classic Venn and Edwards’ Venn diagrams and allow for a maximum of six sets. The software incorporates user-friendly features and works in Windows, Linux and Mac OS. Its graphical interface does not require a user to have programing skills. Users can modify diagram content for up to eight datasets because of the Scalable Vector Graphics output. VennPainter can provide output results in vertical, horizontal and matrix formats, which facilitates sharing datasets as required for further identification of candidate genes. Users can obtain gene lists from shared sets by clicking the numbers on the diagram. Thus, VennPainter is an easy-to-use, highly efficient, cross-platform and powerful program that provides a more comprehensive tool for identifying candidate genes and visualizing the relationships among genes or gene families in comparative analysis.     

Improving gene set analysis of microarray data by SAM-GS

Background  

Gene-set analysis evaluates the expression of biological pathways, or a priori defined gene sets, rather than that of individual genes, in association with a binary phenotype, and is of great biologic interest in many DNA microarray studies. Gene Set Enrichment Analysis (GSEA) has been applied widely as a tool for gene-set analyses. We describe here some critical problems with GSEA and propose an alternative method by extending the individual-gene analysis method, Significance Analysis of Microarray (SAM), to gene-set analyses (SAM-GS).