PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of nihpaAbout Author manuscriptsSubmit a manuscriptHHS Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
 
Obesity (Silver Spring). Author manuscript; available in PMC 2014 March 1.
Published in final edited form as:
PMCID: PMC3762924
NIHMSID: NIHMS430476

A genome wide association study of plasma uric acid levels in obese cases and never-overweight controls

Abstract

Objective

To identify plasma uric acid related genes in extremely obese and normal weight individuals using genome wide association studies (GWAS).

Design and Methods

Using genotypes from a GWAS focusing on obesity and thinness, we performed quantitative trait association analyses (PLINK) for plasma uric acid levels in 1,060 extremely obese individuals [body mass index (BMI) >35 kg/m2] and normal-weight controls (BMI<25kg/m2). In 961 samples with uric acid data, 924 were females.

Results

Significant associations were found in SLC2A9 gene SNPs and plasma uric acid levels (rs6449213, P=3.15×10−12). DIP2C gene SNP rs877282 also reached genome wide significance(P=4,56×10−8). Weaker associations (P<1×10−5) were found in F5, PXDNL, FRAS1, LCORL, and MICAL2genes. Besides SLC2A9, 3 previously identified uric acid related genes ABCG2 (rs2622605, P=0.0026), SLC17A1 (rs3799344, P=0.0017), and RREB1 (rs1615495, P =0.00055) received marginal support in our study.

Conclusions

Two genes/chromosome regions reached genome wide association significance (P< 1× 10−7, 550K SNPs) in our GWAS : SLC2A9, the chromosome 2 60.1 Mb region (rs6723995), and the DIP2C gene region. Five other genes (F5, PXDNL, FRAS1, LCORL, and MICAL2) yielded P<1× 10−5. Four previous reported associations were replicated in our study, including SLC2A9, ABCG2, RREB, and SLC17A1.

Keywords: uric acid, genome wide association study, obesity

Introduction

Uric acid is the end product of purine metabolism. The prevalence of hyperuricemia (uric acid ≥ 420μmol/L in males, ≥ 360μmol/L in females) has increased rapidly over the past two decades (1, 2). The connection between hyperuricemia and gout has long been known; however, hyperuricemia is much more common than gout. Increasing evidence shows that hyperuricemia is a risk factor for metabolic syndrome(3) and cardiovascular diseases(4).

Although obesity and hyperuricemia are correlated, the genetic background of this association is not well understood. Several candidate genes, including SLC2A9 and ABCG2 (5, 6), have been identified in genome-wide association studies (GWASs) and follow-up replications. To investigate the possible role of these genes in obese individuals, we performed a genome wide association study (GWAS) for plasma uric acid in 1,060 obesity cases/controls using our previous genotyping data for body weight traits(7).

Methods and Procedures

Subjects

All subjects gave informed consent, and the protocol was approved by the Committee on Studies Involving Human Beings at the University of Pennsylvania. Five hundred twenty (520) European-American obesity cases (BMI>35 kg/m2) and 540 normal-weight controls (BMI<25 kg/m2) were selected for analysis from ongoing studies (8). Clinical characteristics have been described previously (9). In 961 samples with uric acid data, 924 were females.

Genotyping

DNA was extracted from whole blood or lymphoblastoid cell lines using a high-salt method. All samples were genotyped on Illumina HumanHap550 SNP arrays (Illumina, San Diego, CA) with ~550,000 SNP markers, at the Center for Applied Genomics, Children’s Hospital of Philadelphia.

Data analyses

Uric acid outliers (> 3SD) were deleted from the data set. Quantitative association studies were performed using PLINK 1.07 that based on the Wald test(10). To investigate the plausible influence of obesity status on uric acid levels, we also performed GWAS separately in obesity cases (BMI>35kg/m2) and normal weight controls (BMI<27kg/m2). Female-only analyses were also carried out after quantitative associations were conducted in all samples.

Results

Of the 1,060 obese cases and normal controls, 961 had plasma uric acid data. Thirty-seven (37) of those 961 individuals were male; 924 were female. Average age of the 961 subjects was 41.9±9.1 years (range, 16–65 years). Distributions of uric acid levels in all samples, cases, and controls are shown separately in Table 1. Q-Q plots showed normal distributions of uric acid levels in those 3 groups (Figure 1).

Figure 1
Q-Q plots of plasma uric acid levels in all subjects, obesity cases, and controls
TABLE 1
Traits distributions of plasma uric acid in obese individuals (BMI>35kg/m2), normal weight (BMI<25kg/m2), and combined samples

Significant associations were found between SLC2A9 gene SNPs and plasma uric acid. The most significant result was for the SNP rs6449213 (all samples, P=3.15×10−12; female-only samples, P=2.29×10−12)(Table 2).

TABLE 2
Significant associations between SLC2A9 gene SNPs and uric acid

DIP2C gene SNP rs877282 also reached genome wide significance (P=4,56×10−8). Many SNPs in the DIP2C gene also showed associations (P<1×10−5) (Table 3).

TABLE 3
Quantitative association studies (PLINK) for uric acid levels in obese cases and controls (P<1×10−4)

Weaker associations (P<1×10−5) were found in F5, PXDNL, FRAS1, LCORL, and MICAL2gene SNPs. All 5 genes had multiple SNPs that associated with uric acid levels (3.05×10−6< P <1× 10−4)(Table 3).Three coding region non-synonymous SNPs in the coagulation factor V (F5) gene, rs6030(Met 1764 Val), rs4525 (His 865 Arg), and rs4524 (Lys 858 Arg), were associated with plasma uric acid, P values of those 3 SNPs for BMI adjusted uric acid were 3.05 × 10−6, 0.00018, and 0.00017, respectively.

Besides SLC2A9, 3 previous found uric acid related genes ABCG2(rs2622605, P=0.0026), SLC17A1(rs3799344, P=0.0017), and RREB1 (rs1615495, P =0.00055) received marginal support in our study (Table 4).

TABLE 4
Previous uric acid associated genes were replicated in our GWAS

Discussion

Hyperuricemia has been considered as an independent risk factor of cardiovascular diseases and type 2 diabetes. Single gene mutations, including deficiency of hypoxanthine guanine phosphoribosyltransferase, lead to hyperuricemia; however, the risk attributable to these genes in the general population is minor(11).

Large (>10,000 individuals) GWASs and meta-analyses have shown that many genes are associated with plasma uric acid levels, including eight genes/regions (SLC2A9(5, 12, 13), ABCG2(6), SLC22A11, SLC17A1, GCKR, R3HDM2-INHBC gene region, RREB1, and PDZK1) that exceeded the genome-wide association level (P<10−7)(14). SLC2A9 has the most significant association with uric acid so far, which could explain 3.5% of uric acid variation in the general population (5).

SLC2A9 (GLUT-9) is a major transporter of uric acid. It controls uric acid influx in the basolateral and apical surface of the kidney proximal convoluted tubule (PCT).SLC2A9 is highly expressed in kidney and liver. Interestingly, ABCG2 is an efflux uric acid transporter that is expressed in the apical surface of the PCT. The SLC2A9 and ABCG2 associations are among the strongest of all uric acid associations so far (14).

Uric acid and glucose transport are often coupled, but SLC2A9 is not a major glucose/fructose transporter. In our study, uric acid levels correlated with fasting glucose. It is possible that SLC2A9 polymorphisms account for the uric acid–glucose connection. However, the SLC2A9 gene alone likely does not explain the 20% rate for hyperuricemia and almost the same rate for insulin resistance in general populations. Other genes with relatively minor genetic relative risk and/or gene–gene interactions may account for the rest of the genetic background for hyperuricemia.

The strength of the associations of SLC2A9 gene SNPs and uric acid was well beyond the threshold for genome-wide significance. This is particularly notable given the moderate sample size (961 individuals). The SLC2A9 associations have been replicated in several GWASs and follow-up association studies (5, 6, 13, 14), including European, African-American (15), and Japanese populations. Although this is not the first study to examine a European American population, we are interested in the SLC2A9 association in extremely obese individuals. In our study, the BMI-adjusted uric acid yielded more significant association with SLC2A9 polymorphisms than the unadjusted plasma uric acid. It is said that SLC2A9 is not the major glucose transporter, although it is the main uric acid transporter in proximal convoluted tubule(16). In our subjects, uric acid was correlated with almost all body weight, lipid (except LDL), and insulin resistance phenotypes (P<0.001, data not shown). However, no direct association was found between SLC2A9gene-region SNPs and these other phenotypes(7). These results suggest that the phenotypic associations between uric acid levels and metabolic syndrome phenotypes are through pathways independent of SLC2A9..

All uric acid associated genes found in our GWAS, including SLC2A9, DIP2C (Homo sapiens DIP2 disco-interacting protein 2 homolog C (Drosophila)), F5 (coagulation factor V), FRAS1( Fraser syndrome 1), PXDNL (Homo sapiens peroxidasin homolog (Drosophila)-like), LCORL (ligand dependent nuclear receptor corepressor-like), and MICAL2 (microtubule associated monoxygenase, calponin and LIM domain containing 2), , are expressed in kidney and/or liver. It is hard to predict functional connections among those genes and plasma uric acid levels, although we have already known that some genes have functions in transcription regulations (DIP2Cand LCORL) and mesenchymal/epithelial transition (FRAS1).

Venous thromboembolism, insulin resistance, and hyperuricemia are correlated in general populations. Many studies have shown that Factor V (F5) mutations are associated with factor V Leiden thrombophilia characterizedby deep vein thrombosis (17), however, no established connection between factor V and uric acid has been reported.

The SLC2A9 associations remained significant in both obese cases and controls. Several associations, including MICAL2, FRAS1, and LCORL, were more significant in obese individuals, while F5 was more significant in normal weight controls. (Table 4). Although some of these associations varied in obese cases and controls, however, none of these genes were among the top BMI associations that found in our GWAS (7).

We have failed to replicated associations on SLC22A11, GCKR, and PDZK1 genes that reported by previous large sample sized GWASs(18, 19).We could not explain if those lack of association were due to a smaller sample size, but no marginal significant association (P<0.05) was found in either original or BMI adjusted uric acid levels.

In summary, two genes/chromosome regions reached genome wide association significance (P< 1× 10−7, 550K SNPs) in our GWAS : SLC2A9, the chromosome 2 60.1 Mb region (rs6723995), and the DIP2C gene region. Five other genes (F5, PXDNL, FRAS1, LCORL, and MICAL2) yielded P<1× 10−5. Four previous reported associations were replicated in our study, including SLC2A9, ABCG2, RREB, and the SLC17A1.

Acknowledgments

We thank all subjects who donated blood samples for genetic research purposes. This work was supported in part by NIH grants R01DK44073, R01DK56210, and R01DK076023 to R.A.P., a Scientist Development Grant (0630188N) from the American Heart Association, a grant (81070576) from the National Natural Science Foundation of China (NSFC), and a grant (12JCZDJC24700) from Tianjin Municipal Science and Technology Commission to W.D.L. Genome-wide genotyping was funded in part by an Institutional Development Award to the Center for Applied Genomics (H.H.) from the Children’s Hospital of Philadelphia.

Footnotes

Conflicts of Interest Statement

The authors declare that there is no duality of interest associated with this manuscript.

Author contributions: Wei-Dong Li: study design, oversee and conduct experiments, analyze data, write the manuscript; Hongxiao Jiao: analyze data, comment on the manuscript; Kai Wang: analyze data, comment on the manuscript;

Clarence Zhang: analyze data; Joseph T Glessner: oversee and conduct experiments; Struan F.A. Grant: oversee and conduct experiments; Hongyu Zhao: oversee and conduct data analyses; Hakon Hakonarson: oversee and conduct experiment, study design; R. Arlen Price: study design, oversee experiments and analyze data, write and comment on the manuscript.

References

1. Wallace KL, Riedel AA, Joseph-Ridge N, Wortmann R. Increasing prevalence of gout and hyperuricemia over 10 years among older adults in a managed care population. J Rheumatol. 2004;31(8):1582–7. Epub 2004/08/04. [PubMed]
2. Nan H, Qiao Q, Dong Y, Gao W, Tang B, Qian R, et al. The prevalence of hyperuricemia in a population of the coastal city of Qingdao, China. J Rheumatol. 2006;33(7):1346–50. Epub 2006/07/06. [PubMed]
3. Choi HK, Ford ES. Prevalence of the metabolic syndrome in individuals with hyperuricemia. Am J Med. 2007;120(5):442–7. Epub 2007/05/01. [PubMed]
4. Zheng LQ, Li J, Yu JM, Hasimu B, Hu DY. Study on the independent association of uric acid levels with peripheral arterial disease in Chinese patients with coronary heart disease. Zhonghua Liu Xing Bing Xue Za Zhi. 2006;27(2):161–4. Epub 2006/06/06. [PubMed]
5. Doring A, Gieger C, Mehta D, Gohlke H, Prokisch H, Coassin S, et al. SLC2A9 influences uric acid concentrations with pronounced sex-specific effects. Nat Genet. 2008;40(4):430–6. Epub 2008/03/11. [PubMed]
6. Dehghan A, Kottgen A, Yang Q, Hwang SJ, Kao WL, Rivadeneira F, et al. Association of three genetic loci with uric acid concentration and risk of gout: a genome-wide association study. Lancet. 2008;372(9654):1953–61. Epub 2008/10/07. [PMC free article] [PubMed]
7. Wang K, Li WD, Zhang CK, Wang Z, Glessner JT, Grant SF, et al. A genome-wide association study on obesity and obesity-related traits. PloS one. 2011;6(4):e18939. Epub 2011/05/10. [PMC free article] [PubMed]
8. Price RA, Reed DR, Lee JH. Obesity related phenotypes in families selected for extreme obesity and leanness. Int J Obes Relat Metab Disord. 1998;22(5):406–13. [PubMed]
9. Wang K, Li WD, Glessner JT, Grant SF, Hakonarson H, Price RA. Large copy-number variations are enriched in cases with moderate to extreme obesity. Diabetes. 59(10):2690–4. Epub 2010/07/14. [PMC free article] [PubMed]
10. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007;81(3):559–75. Epub 2007/08/19. [PubMed]
11. Reed DR, Price RA. X-linkage does not account for the absence of father-son similarity in plasma uric acid concentrations. Am J Med Genet. 2000;92(2):142–6. Epub 2000/05/08. [PubMed]
12. Wallace C, Newhouse SJ, Braund P, Zhang F, Tobin M, Falchi M, et al. Genome-wide association study identifies genes for biomarkers of cardiovascular disease: serum urate and dyslipidemia. Am J Hum Genet. 2008;82(1):139–49. Epub 2008/01/09. [PubMed]
13. Brandstatter A, Kiechl S, Kollerits B, Hunt SC, Heid IM, Coassin S, et al. Sex-specific association of the putative fructose transporter SLC2A9 variants with uric acid levels is modified by BMI. Diabetes care. 2008;31(8):1662–7. Epub 2008/05/20. [PMC free article] [PubMed]
14. Yang Q, Kottgen A, Dehghan A, Smith AV, Glazer NL, Chen MH, et al. Multiple genetic loci influence serum urate levels and their relationship with gout and cardiovascular disease risk factors. Circ Cardiovasc Genet. 3(6):523–30. Epub 2010/10/05. [PMC free article] [PubMed]
15. Charles BA, 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. Epub 2011/02/08. [PMC free article] [PubMed]
16. Vitart V, Rudan I, Hayward C, Gray NK, Floyd J, Palmer CN, et al. SLC2A9 is a newly identified urate transporter influencing serum urate concentration, urate excretion and gout. Nat Genet. 2008;40(4):437–42. Epub 2008/03/11. [PubMed]
17. Zee RY, Glynn RJ, Cheng S, Steiner L, Rose L, Ridker PM. An evaluation of candidate genes of inflammation and thrombosis in relation to the risk of venous thromboembolism: The Women’s Genome Health Study. Circulation Cardiovascular genetics. 2009;2(1):57–62. Epub 2009/12/25. [PMC free article] [PubMed]
18. Kolz M, Johnson T, Sanna S, Teumer A, Vitart V, Perola M, et al. Meta-analysis of 28,141 individuals identifies common variants within five new loci that influence uric acid concentrations. PLoS genetics. 2009;5(6):e1000504. Epub 2009/06/09. [PMC free article] [PubMed]
19. Yang Q, Kottgen A, Dehghan A, Smith AV, Glazer NL, Chen MH, et al. Multiple genetic loci influence serum urate levels and their relationship with gout and cardiovascular disease risk factors. Circulation Cardiovascular genetics. 2010;3(6):523–30. Epub 2010/10/05. [PMC free article] [PubMed]