This study aimed to identify functional associations of cis-regulatory regions with gout susceptibility using data resulted from a genome-wide association study (GWAS), and to show a genetic architecture for gout with interaction effects among genes within each of the identified functions. The GWAS was conducted with 8314 control subjects and 520 patients with gout in the Korea Association REsource cohort. However, genetic associations with any individual nucleotide variants were not discovered by Bonferroni multiple testing in the GWAS (P>1.42 × 10(-7)). Genomic regions enrichment analysis was employed to identify functional associations of cis-regulatory regions. This analysis revealed several biological processes associated with gout susceptibility, and they were quite different from those with serum uric acid level. Epistasis for susceptibility to gout was estimated using entropy decomposition with selected genes within each biological process identified by the genomic regions enrichment analysis. Some epistases among nucleotide sequence variants for gout susceptibility were found to be larger than their individual effects. This study provided the first evidence that genetic factors for gout susceptibility greatly differed from those for serum uric acid level, which may suggest that research endeavors for identifying genetic factors for gout susceptibility should not be heavily dependent on pathogenesis of uric acid. Interaction effects between genes should be examined to explain a large portion of phenotypic variability for gout susceptibility.