Genome-Wide Analyses of ChIP-Seq Derived FOXA2 DNA Occupancy in Liver Points to Genetic Networks Underpinning Multiple Complex Traits

Background: Forkhead Box A2 (FOXA2) exerts an influence on glucose homeostasis via activity in the liver. In addition, a key genome-wide association study (GWAS) recently demonstrated that genetic variation, namely rs6048205, at the FOXA2 locus is robustly associated with fasting glucose levels. Our hypothesis was that this DNA-binding protein regulates the expression of a set of molecular pathways critical to endocrine traits. Methods: Drawing on our laboratory and bioinformatic experience with chromatin immunoprecipitation followed by massively parallel sequencing, we analyzed our existing FOXA2 chromatin immunoprecipitation followed by massively parallel sequencing data generated in human liver, using the algorithm hypergeometric optimization of motif enrichment, to gain insight into its global genomic binding pattern from a disease perspective. Results: We performed a pathway analysis of the gene list using the gene set enrichment analysis algorithm, which yielded a number of significant annotations. Motivated by the fact that the FOXA2 locus has been implicated by GWAS, we cross-referenced the occupancy sites with the National Institutes of Health GWAS catalog and found strong evidence for the enrichment of loci implicated in endocrine, neuropsychiatric, cardiovascular, and cancer trait categories, but interestingly there was no evidence for enrichment for inflammation related traits. Intriguingly, a FOXA2 occupancy site coincided with rs6048205, suggesting that this variant confers its effect, at least partially, via a perturbation of a FOXA2 feedback mechanism. Conclusion: Our data strongly suggest that FOXA2 is acting as a master regulator of key pathways that are enriched for loci implicated by GWAS for most trait categories, with the clear exception of inflammation, suggesting that this factor exerts its effect in this context via noninflammatory processes.