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.

[1]  F. Tronche,et al.  Plasticity and expanding complexity of the hepatic transcription factor network during liver development. , 2006, Genes & development.

[2]  K. Kaestner,et al.  The FoxA factors in organogenesis and differentiation. , 2010, Current opinion in genetics & development.

[3]  S. Grant,et al.  Generation of high quality chromatin immunoprecipitation DNA template for high-throughput sequencing (ChIP-seq). , 2013, Journal of visualized experiments : JoVE.

[4]  Luigi Ferrucci,et al.  Population-based genome-wide association studies reveal six loci influencing plasma levels of liver enzymes. , 2008, American journal of human genetics.

[5]  K. Kaestner,et al.  Unified nomenclature for the winged helix/forkhead transcription factors. , 2000, Genes & development.

[6]  Gavin Giovannoni,et al.  A ChIP-seq defined genome-wide map of vitamin D receptor binding: associations with disease and evolution. , 2010, Genome research.

[7]  S. Grant,et al.  Disease-associated loci are significantly over-represented among genes bound by transcription factor 7-like 2 (TCF7L2) in vivo , 2010, Diabetologia.

[8]  Logan J Everett,et al.  Species-specific strategies underlying conserved functions of metabolic transcription factors. , 2011, Molecular endocrinology.

[9]  Raymond K. Auerbach,et al.  An Integrated Encyclopedia of DNA Elements in the Human Genome , 2012, Nature.

[10]  M. Karin,et al.  IKKα activation of NOTCH links tumorigenesis via FOXA2 suppression. , 2012, Molecules and Cells.

[11]  Mark I. McCarthy,et al.  Pancreatic islet enhancer clusters enriched in type 2 diabetes risk–associated variants , 2013, Nature Genetics.

[12]  M. Daly,et al.  PGC-1α-responsive genes involved in oxidative phosphorylation are coordinately downregulated in human diabetes , 2003, Nature Genetics.

[13]  Martin Renqiang Min,et al.  An integrated encyclopedia of DNA elements in the human genome , 2012 .

[14]  Data production leads,et al.  An integrated encyclopedia of DNA elements in the human genome , 2012 .

[15]  H. Stefánsson,et al.  Variant of transcription factor 7-like 2 (TCF7L2) gene confers risk of type 2 diabetes , 2006, Nature Genetics.

[16]  Shane J. Neph,et al.  Systematic Localization of Common Disease-Associated Variation in Regulatory DNA , 2012, Science.

[17]  Pablo Tamayo,et al.  Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[18]  Jonathan P. Katz,et al.  Inactivation of the winged helix transcription factor HNF3alpha affects glucose homeostasis and islet glucagon gene expression in vivo. , 1999, Genes & development.

[19]  J. Whitsett,et al.  Foxa2 Programs Th2 Cell-Mediated Innate Immunity in the Developing Lung , 2010, The Journal of Immunology.

[20]  R. Costa,et al.  The restricted promoter activity of the liver transcription factor hepatocyte nuclear factor 3 beta involves a cell-specific factor and positive autoactivation , 1992, Molecular and cellular biology.

[21]  C. Glass,et al.  Simple combinations of lineage-determining transcription factors prime cis-regulatory elements required for macrophage and B cell identities. , 2010, Molecular cell.

[22]  Claude Bouchard,et al.  A genome-wide approach accounting for body mass index identifies genetic variants influencing fasting glycemic traits and insulin resistance , 2012, Nature Genetics.