Abstract
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.
Similar content being viewed by others
Introduction
Gout has been a great concern as the most prevalent arthritis associated with impaired quality of life.1, 2 Gout is developed with a preliminary increase of serum uric acid level and a subsequent deposition of monosodium urate crystals in and around joints. The serum uric acid level is, therefore, considered the most critical risk factor for susceptibility to gout. As inefficient excretion of uric acid is more influential in elevating serum uric acid level and predisposing people to gout than its overproduction, renal mechanisms have received an overwhelming attention in terms of explaining gout susceptibility.3 Especially, recent genetic studies have primarily focused on genes regulating renal urate transport.4 For example, genetic associations of serum uric acid level and/or gout susceptibility were repeatedly observed with nucleotide sequence variants of the gene encoding GLUT9 (SLC2A9), a urate transporter in proximal renal tubules.5, 6, 7 Genetic associations were also identified with other transporters such as ABCG2,5, 8 NPT1 (SLC17A1),9 SLC17A35 and URAT1 (SLC22CA12).10, 11 Phenotypic variability for gout susceptibility has been quite limitedly explained by a small number of genetic factors identified so far,12 especially considering its large heritability (0.63).13 A small number of genetic associations have been discovered even in recent genome-wide association studies (GWAS).14 This was because many false-negative genetic associations were produced by a multiple test with hundreds of thousands of nucleotide markers.
The current study aimed to conduct a GWAS for gout susceptibility in a large-scale cohort, to identify functional associations of cis-regulatory regions with gout susceptibility, to compare the functional associations with gout susceptibility with those with serum uric acid level, and to show a genetic architecture for gout susceptibility with interaction effects among genes within each functional association.
Materials and methods
Subjects and data
The Korea Association REsource (KARE) Analysis Consortium was established to conduct large-scale GWAS for understanding of human genetic basis. The consortium has cohort data of 10 038 unrelated Korean individuals collected by the Korean Genome Epidemiology Study. The data collection was initiated in 2001, and thereafter, follow-up examinations for each participant have been conducted every 2 years. Genotypic data of the KARE were obtained using the Affymetrix Genome-Wide Human SNP Array 5.0 (Affymetrix, Santa Clara, CA, USA). A total of 8842 individuals were obtained after screening by genotype calling and quality control. For details, see Cho et al.15 However, four individuals without phenotype or covariate information were excluded in the current analysis. An underlying set of unphased genotypes for each individual in the cohort was imputed with the Japanese and Chinese HapMap phase 2 haplotype panel using IMPUTE software program (version 2).16
A total of 520 out of 8834 subjects were self-reported patients with gout, and they were assumed to be diagnosed with gout by rheumatologists according to the American College of Rheumatology diagnostic criteria.17 The other subjects in the cohort were all used as controls. The patients with gout were excluded in association analysis for serum uric acid level. Phenotypic characteristics of the patients were compared with those of controls in Supplementary Table 1.
Single locus association
Genotypic associations of each single-nucleotide polymorphism with susceptibility to gout and with serum uric acid level were tested with an additive model. The analytical model included age, gender and body mass index as covariates. Threshold of false-positive error in the significance test was 0.05, and the Bonferonni correction was introduced to control for occurrence of false-positive genetic associations. The association analyses were conducted using PLINK (version 1.06, http://pngu.mgh.harvard.edu/purcell/plink) and SPSS (version 12.0, SPSS, Chicago, IL, USA) software programs.
Genome-wide functional prediction
Genomic regions enrichment analysis was conducted to analyze functional associations of cis-regulatory regions identified by localized measurements of DNA-binding events using GREAT (http://great.stanford.edu/).18 The input genomic regions were determined on the basis of single-nucleotide polymorphisms with P<0.01 resulted from the GWAS. Biological processes provided by the Gene Ontology (http://www.geneontology.org/) were utilized as input annotation terms. The analysis incorporated distal binding sites up to 1 Mb, as well as proximal binding sites up to upstream 5 Kb and downstream 1 Kb. Significance of enrichment was determined by both a binomial test over genomic regions and a hypergeometric test over genes. To control incorrect rejections of null hypothesis, the significance thresholds were corrected by false discovery rate.
Entropy decomposition
Phenotypic variability of gout as a complex disease might be largely explained by epistasis. The epistatic effects among genes within each biological process were estimated in the current study. We first selected the optimal genes within each biological process for susceptibility to gout using multifactor dimensionality reduction (MDR) analysis and then estimated pairwise epistatic estimates using entropy decomposition.
The MDR collapses high-order dimensions of genotypic data produced with multiple loci into a single dimension, which enables to infer epistasis in a relatively small sample size by grouping genotypes.19 Thus, the MDR could identify multiple loci and their genotype combinations that were associated with gout susceptibility. A 10-fold cross-validation strategy was applied to the current MDR. Case and control data were randomly and equally partitioned into 10 pieces. Nine pieces were used as a training data set to build a genetic model for predicting susceptibility to gout, and the rest one piece was used as a testing data set to test the model established by using the training set. This analysis was conducted with each of the 10 possible partitions. The best n-locus model was selected on the basis of minimum classification error with 10-fold cross validation. Each of combined genotypes produced with every possible n-locus was determined as to whether it was a high or low risk, based on the case–control ratio. All the procedures explained above were replicated 100 times through shuffling data sets. The best models with 1-, 2-, 3- and 4-loci were determined with the maximum average estimates of cross-validation consistency and testing accuracy for the replicates. All the variants in the best models with 1-, 2-, 3- and 4-loci were selected for estimating their additive and interactive genetic effects. The MDR analysis was conducted with all the loci of each function identified by the genomic regions enrichment analysis using the MDR software package (http://www.multifactordimensionalityreduction.org).
Pairwise interactions among the selected loci associated with a susceptibility to gout were estimated as the portions of entropy removed for interaction information of pair variants after adjusting marginal entropy removed by each individual variant.20 Therefore, interaction analysis reduced uncertainty about either of two attributes with the knowledge of the other attribute. The interaction effect was partitioned into synergistic and redundant effects by the sign of its value. The entropy decomposition analysis was conducted utilizing the freely available Orange software (http://www.ailab.si/orange).
Results
The GWAS revealed no significant nucleotide sequence variants associated with gout susceptibility or with serum uric acid level by Bonferroni multiple test (P>1.42 × 10−7, Figure 1). Previously identified associations of gout susceptibility or serum uric acid level were replicated in the current study with nucleotide sequence variants in solute carrier family genes such as GLUT9, NPT1, PDZK1 and SLC17A3 (P<0.05, Table 1).
Genome-region set analysis revealed a variety of biological processes associated (false-discovery-rate-corrected P<0.05) with each of gout susceptibility (Table 2) and serum uric acid level (Table 3). The identified biological processes associated with gout susceptibility differed from those with serum uric acid level except for regulation of chemotaxis.
Individual and interactive genetic effects of nucleotide sequence variants were estimated for each of the identified biological processes using entropy decomposition (Figure 2). The entropy decomposition analysis included all the loci of the best 1-, 2-, 3- and 4-locus models selected by MDR (Supplementary Table 3). Some interactions were observed with entropy estimates larger than those of individual effects. For example, rs7035561 in AGTPBP1 gene and rs9960723 in NETO1 gene showed a strong synergistic epistasis, which explained phenotypic variability approximately four times larger than their individual effects. Furthermore, although marginal effect of rs9308056 in NPY5R was negligible, its interaction with the rs7035561 was considerably large. In general, strong epistatic associations were synergistic in regulations of neurological system process, cellular metabolic process and neurotransmitter transport, whereas redundant epistatic associations were observed in regulations of leukocyte migration and chemotaxis.
Discussion
We conducted a GWAS for gout susceptibility and serum uric acid level with a large-cohort data set in the current study, but any nucleotide sequence variants were not identified by employing Bonferroni multiple test. In spite of their high heritability (0.63),13 complexity of the traits might lead to hard detection for genetic factors with small effects. A large number of false-negative associations could be produced by such a conservative testing. We applied a genomic regions enrichment analysis using variants with a marginal association to overcome the problem. Furthermore, we could exclude potential biases by incorporating two mutually beneficial tests: one was a binomial test explicitly accounting for the sizes of the regulatory domains of genes, and the other was a hypergeometric test counting each gene only once. As a result, a variety of annotation terms were identified with statistical enrichment for their functional associations with gout susceptibility. The annotation terms turned out to be quite different from those associated with serum uric acid level. This kind of a large divergence was unexpected, considering the serum uric acid level as the underlying cause of gout. Nevertheless, the annotation terms were identified for association with gout susceptibility, and specifically with serum uric acid level. For example, the annotation terms included ‘transforming growth factor (TGF)β receptor signaling pathway’ for serum uric acid level. This was explained well with the role of TGFβ1 as a promoter in synthesizing extracellular matrix. The change of extracellular matrix declines renal function and increases serum uric acid level.21 On the other hand, gout susceptibility was not associated with the term of ‘TGFβ receptor signaling pathway’, but significantly associated with the annotation term of ‘leukocyte migration’. Leukocyte adhesion in monosodium urate crystal-induced inflammation can be suppressed by TGFβ1, and the suppression would lead to gout.22 The discrepancy between the annotation terms associated with gout susceptibility and with serum uric acid level implied heterogenous underlying genetic factors for the gout and its underlying cause. The different genetic associations might explain that many people with hyperuricemia would be highly unlikely to develop gout. That is, gout susceptibility might be largely explained by a variety of genetic factors for inducing chronic inflammation after urate deposition.
Phenotypic variability of gout as a complex disease might be largely explained by interaction effects between genes.23 Including the interaction effects is essential in understanding genetic architecture for gout susceptibility, but may present a serious drawback of containing a large number of all the possible interaction effects in the analytical model. To overcome the problem, we employed a parsimonious interaction model with selected genes within each of the biological process terms in the current study. As a result, large interactive effects of the genetic factors for gout susceptibility were observed using an entropy-based approach. Especially, some interaction effects explained larger phenotypic variability than their marginal effects. For example, neurological system process-related genes showed their strong synergistic epistasis for gout susceptibility.
One thing to be noted is a potential heterogeneity of gout genetic architecture by gender. The current study was not balanced by gender, and thus, a larger portion of genetic factors was attributed to genetic mechanism for susceptibility to gout in males. However, this unbalance in the current study was not too extreme with a larger female ratio of gout (21%; Supplementary Table 1) than that in usual (5%).24 The larger proportion might be explained partially by having a larger female ratio (10%)25 of gout in Korea, and more importantly by including only healthy subjects on the starting point of the cohort, that is, male subjects with gout at that point might have been excluded.
In conclusion, the current study provided the first evidence that genetic factors for gout susceptibility considerably differed from those for serum uric acid level. This implied that gout susceptibility was greatly attributed to genetic factors irrelevant to pathogenesis of uric acid, and thus understanding genetic architecture of gout susceptibility would also require research endeavors on identification of such genetic factors. Interaction effects between genes turned out to be not negligible for explaining phenotypic variability for gout susceptibility. To validate our findings, the current study calls for rigorous replications in other ethnic populations. Furthermore, comparing subjects with gout with those with asymptomatic hyperuricemia would confirm the heterogeneous genetic factors for gout susceptibility and serum uric acid level, which eventually help understand decomposed genetic factors for such complex disease.
References
Roddy, E., Zhang, W. & Doherty, M. Is gout associated with reduced quality of life? A case–control study. Rheumatology 46,, 1441–1444 (2007).
Lee, S. J., Hirsch, J. D., Terkeltaub, R., Khanna, D., Singh, A. J., Andrew, S. et al. Perceptions of disease and health-related quality of life among patients with gout. Rheumatology 48, 582–586 (2009).
Richette, P. & Bardin, T. Gout. Lancet 375, 318–328 (2010).
So, A. & Thorens, B. Uric acid transport and disease. J. Clin. Invest. 120, 1791–1799 (2010).
Dehghan, A., Kottgen, A., Yang, Q., Hwang, S. J., Kao, W. L., Rivadeneira, F. et al. Association of three genetic loci with uric acid concentration and risk of gout: a genome-wide association study. Lancet 372, 1953–1961 (2008).
Doring, A., Gieger, C., Mehta, D., Prokisch, H., Coassin, S., Fischer, G. et al. SLC2A9 influences uric acid concentrations with pronounced sex-specific effects. Nat. Genet. 40, 430–436 (2008).
Vitart, V., Rudan, I., Hayward, C., Gray, N. K., Floyd, J., Palmer, C. N. et al. SLC2A9 is a newly identified urate transporter influencing serum urate concentration, urate excretion and gout. Nat. Genet. 40, 437–442 (2008).
Woodward, O. M., Köttgen, A., Coresh, J., Boerwinkle, E., Guggino, W. B. & Köttgen, M. Identification of a urate transporter, ABCG2, with a common functional polymorphism causing gout. Proc. Natl Acad. Sci. USA 106, 10338–10342 (2009).
Urano, W., Taniguchi, A., Anzai, N., Inoue, E., Kanai, Y., Yamanaka, M. et al. Sodium-dependent phosphate co-transporter type 1 sequence polymorphisms in male patients with gout. Ann. Rheum. Dis. 69, 1232–1234 (2010).
Graessler, J., Graessler, A., Unger, S., Kopprasch, S., Tausche, A. K., Kuhlisch, E. et al. Association of the human urate transporter 1 with reduced renal uric acid excretion and hyperuricemia in a German Caucasian population. Arthritis Rheum. 54, 292–300 (2006)).
Vázquez-Mellado, J., Jiménez-Vaca, A. L., Cuevas-Covarrubias, S., Alvarado-Romano, V., Pozo-Molina, G. & Burgos-Vargas, R. Molecular analysis of the SLC22A12 (URAT1) gene in patients with primary gout. Rheumatology 46, 215–219 (2007).
van der Harst, P., Bakker, S. J., de Boer, R. A., Wolffenbuttel, B. H., Johnson, T., Caulfield, M. J. et al. Replication of the five novel loci for uric acid concentrations and potential mediating mechanisms. Hum. Mol. Genet 19, 387–395 (2010).
Yang, Q., Guo, C. Y., Cupples, L. A., Levy, D., Wilson, P. W. & Fox, C. S. Genome-wide search for genes affecting serum uric acid levels: the Framingham Heart Study. Metabolism 54, 1435–1441 (2005).
Sulem, P., Gudbjartsson, D. F., Walters, G. B., Helgadottir, H. T., Helgason, A., Gudjonsson, S. A. et al. Identification of low-frequency variants associated with gout and serum uric acid levels. Nat. Genet. 43, 1127–1131 (2011).
Cho, Y. S., Go, M. J., Kim, Y. J., Heo, J. Y., Oh, J. H., Ban, H. J. et al. A large-scale genome-wide association study of Asian populations uncovers genetic factors influencing eight quantitative traits. Nat. Genet. 41, 527–534 (2009).
Ryoo, H., Woo, J., Kim, Y. & Lee, C. Heterogeneity of genetic associations of CDKAL1 and HHEX with susceptibility of type 2 diabetes mellitus by gender. Eur. J. Hum. Genet. 19, 672–675 (2011).
Wallace, S. L., Robinson, H., Masi, A. T., Decker, J. L., McCarty, D. J. & Yu, T. F. Preliminary criteria for the classification of the acute arthritis of primary gout. Arthritis Rheum. 20, 895–900 (1977).
McLean, C. Y., Bristor, D., Hiller, M., Clarke, S. L., Schaar, B. T., Lowe, C. B. et al. GREAT improves functional interpretation of cis-regulatory regions. Nat. Biotechnol 28, 495–501 (2010).
Hahn, L. W., Ritchie, M. D. & Moore, J. H. Multifactor dimensionality reduction software for detecting gene-gene and gene-environment interactions. Bioinformatics 19, 376–382 (2003).
Park, J., Kim, Y. & Lee, C. Identification of epistasis in ischemic stroke using multifactor dimensionality reduction and entropy decomposition. BMB Rep. 42, 617–622 (2009).
Basile, D. P. The transforming growth factor beta system in kidney disease and repair: recent progress and future directions. Curr. Opin. Nephrol. Hypertens. 8, 21–30 (1999).
Chang, S. J., Chen, C. J., Tsai, F. C., Lai, H. M., Tsai, P. C., Tsai, M. H. et al. Associations between gout tophus and polymorphisms 869T/C and −509C/T in transforming growth factor β1 gene. Rheumatology 47, 617–621 (2008).
Lee, C. & Kim, Y. Optimal designs for estimating and testing interaction among multiple loci in complex traits by a Gibbs sampler. Genomics 92, 446–451 (2008).
Edwards, N. L. Clinical gout in: Rheumatology. 5th edn (eds Hochberg M. C., Silman A. J., Smolen J. S., Weinblatt M. E., Weisman M. H.) 1859–1865 (Elsevier, Philadelphia, PA, 2011).
Seo, Y. I., Son, K. M., Jung, Y. O., Kim, I. J., Kim, K. L., Bae, Y. D. et al. Survey of current trends for diagnosis and treatment in Korean gout patients. J. Rheum. Dis. 18, 187–192 (2011).
Gunjaca, G., Boban, M., Pehlic, M., Zemunik, T., Budimir, D., Kolcić, I. et al. Predictive value of 8 genetic loci for serum uric acid concentration. Croat. Med. J 51, 23–31 (2010).
Charles, B. A., Shriner, D., Doumatey, A., Chen, G., Zhou, J., Huang, H. et al. A genome-wide association study of serum uric acid in African Americans. BMC Med. Genomics 4, 17 (2011).
Acknowledgements
We thank the National Institute of Health, Korea, for providing the genotypic and epidemiological data to the KARE Analysis Consortium. This work was supported by Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education, Science and Technology (2009–0071063).
Author information
Authors and Affiliations
Corresponding author
Ethics declarations
Competing interests
The authors declare no conflict of interest.
Additional information
Supplementary Information accompanies the paper on Journal of Human Genetics website
Supplementary information
Rights and permissions
About this article
Cite this article
Shin, J., Kim, Y., Kong, M. et al. Genetic architecture for susceptibility to gout in the KARE cohort study. J Hum Genet 57, 379–384 (2012). https://doi.org/10.1038/jhg.2012.39
Received:
Revised:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/jhg.2012.39
Keywords
This article is cited by
-
An update on the genetics of hyperuricaemia and gout
Nature Reviews Rheumatology (2018)
-
The frequency of single nucleotide polymorphisms and their association with uric acid concentration based on data from genome-wide association studies in the Korean population
Rheumatology International (2014)