- Research article
- Open Access
- Open Peer Review
A genome-wide association study identifies common variants influencing serum uric acid concentrations in a Chinese population
BMC Medical Genomics volume 7, Article number: 10 (2014)
Uric acid (UA) is a complex phenotype influenced by both genetic and environmental factors as well as their interactions. Current genome-wide association studies (GWASs) have identified a variety of genetic determinants of UA in Europeans; however, such studies in Asians, especially in Chinese populations remain limited.
A two-stage GWAS was performed to identify single nucleotide polymorphisms (SNPs) that were associated with serum uric acid (UA) in a Chinese population of 12,281 participants (GWAS discovery stage included 1452 participants from the Dongfeng-Tongji cohort (DFTJ-cohort) and 1999 participants from the Fangchenggang Area Male Health and Examination Survey (FAMHES). The validation stage included another independent 8830 individuals from the DFTJ-cohort). Affymetrix Genome-Wide Human SNP Array 6.0 chips and Illumina Omni-Express platform were used for genotyping for DFTJ-cohort and FAMHES, respectively. Gene-environment interactions on serum UA levels were further explored in 10,282 participants from the DFTJ-cohort.
Briefly, we identified two previously reported UA loci of SLC2A9 (rs11722228, combined P = 8.98 × 10-31) and ABCG2 (rs2231142, combined P = 3.34 × 10-42). The two independent SNPs rs11722228 and rs2231142 explained 1.03% and 1.09% of the total variation of UA levels, respectively. Heterogeneity was observed across different populations. More importantly, both independent SNPs rs11722228 and rs2231142 were nominally significantly interacted with gender on serum UA levels (P for interaction = 4.0 × 10-2 and 2.0 × 10-2, respectively). The minor allele (T) for rs11722228 in SLC2A9 has greater influence in elevating serum UA levels in females compared to males and the minor allele (T) of rs2231142 in ABCG2 had stronger effects on serum UA levels in males than that in females.
Two genetic loci (SLC2A9 and ABCG2) were confirmed to be associated with serum UA concentration. These findings strongly support the evidence that SLC2A9 and ABCG2 function in UA metabolism across human populations. Furthermore, we observed these associations are modified by gender.
Uric acid (UA) is the primary end-product of purine metabolism in human beings. Most of the UA is derived from the metabolism of endogenous purine including cell turnover and synthesis. The UA excretion and reabsorption is mostly in kidney [1, 2]. The UA concentration in human blood is more than fifty times higher than that in other mammals, because in most of the animals UA could be further catalyzed to allantoin by urate oxidase or uricase (the copper-binding enzyme) [1, 3]; however, human beings lack uricase and have higher UA levels which could result in hyperuricemia and gouty arthritis . UA could serve as an antioxidant by removing singlet oxygen and radicals . However, elevated UA concentration can lead to a variety of disorders, including gout, hypertension, metabolic syndrome, diabetes mellitus, and cardiovascular disease [6–10].
It is indicated that the conventional factors including age, body mass index (BMI), alcohol consumption, and cigarette smoking could influence serum UA concentrations [11–16]. In addition, serum UA levels were also determined by genetic factors with heritability ranged from 25% to 63% [17–19]. Recent GWASs have identified multiple loci of ABCG2, SLC2A9, SLC17A1, SLC16A9, SLC22A11, GCKR, PDZK1, GCKR, RREB1, LRRC16A, WDR1, TRIM46, INHBB, SFMBT1, TMEM171, VEGFA, BAZ1B, PRKAG2, STC1, HNF4G, A1CF, ATXN2, UBE2Q2, IGF1R, NFAT5, MAF, HLF, ACVR1B-ACVRL1, and B3GNT4 associated with UA levels in European [20–26]. However, only two GWA studies have been conducted in Asians [26, 27]. In addition, previous study indicated that the serum UA concentration was influenced by gene-environmental interactions . In order to investigate the genetic determines of serum UA in Asians, especially in Chinese, we conducted a GWA study in 3,451 individuals, and subsequently replicated the top SNPs in additional 8,830 healthy participants. In addition, we further examined whether the top SNPs interacted with gender, BMI, alcohol drinking, and cigarette smoking in determining serum UA levels respectively.
In the discovery stage, we performed a GWAS of two studies in the Chinese Han population: the DFTJ-cohort consisted of 1,461 healthy individuals and the FAMHES included 2,012 Han healthy individuals aged 20 to 69 years old. The DFTJ cohort  and the FAMHES  were described in detail elsewhere. All of the participants included in the GWAS stage were recruited at health check-ups without chronic diseases such as cardiovascular disease and cancer. Briefly, the DFTJ-cohort initiated in 2008 is a long-term prospective, population-based cohort study designed to determine the gene-environmental interaction on several chronic diseases (obesity, diabetes mellitus, cardiovascular disease, etc.) and cancer in employees from Dongfeng Motor Corporation (DMC). The FAMHES was launched in 2009 in Fangchenggang city, Guangxi, southwest China in 2009 and enrolled 4,303 Chinese men with 17 to 88 years-old; this study was designed to examine the genetic, environmental, and their interactions on the development of age-related chronic diseases. The 8,830 healthy individuals included in the validation stage were selected from the DFTJ-cohort excluding the initial 1,461 subjects and had no diagnosed chronic diseases such as cardiovascular disease, cancer, and gout et al.
The detailed information about the GWAS population and the replication samples is shown in (Additional file 1: Table S1). All the participants provided written informed consent and the ethical committees in the Tongji Medical College and Guangxi Medical University approved this research project.
Measurement of serum UA levels and the covariates
A baseline physical examination was conducted and demographic information was collected via standard questionnaire. Overnight fasting venous blood specimens were obtained and serum UA levels were measured by the ARCHITECT Ci8200 automatic analyzer (ABBOTT Laboratories. Abbott Park, Illinois, U.S.A) using the Abbott Diagnostics reagents following the manufacturer’s instructions in the DFTJ cohort  and using automatic analyzer (Dade Behring, USA) with original reagents in the FAMHES study. Weight and standing height were measured with light indoor clothing and in bare feet. Those who had smoked at least one cigarette per day for more than half a year either currently or formerly were defined as smokers; otherwise they were viewed as non-smokers. Alcohol drinking was divided into two categories: drinkers and non-drinkers. Those who had drunk at least once a week for more than half a year either currently or formerly were defined as drinkers; otherwise they were viewed as non-drinkers.
Sample genotyping and quality control
We performed the GWAS scan in 1,461 subjects from the DFTJ-cohort using Affymetrix Genome-Wide Human SNP Array 6.0 chips following the manufacture’s protocol. Totally, we genotyped 906,703 SNPs among 1,461 subjects. After stringent QC filtering individuals with genotyping call rate < 95% were excluded (9 subjects) for further analysis. SNPs were excluded when 1) MAF < 0.01; 2) Hardy-Weinberg Equilibrium (HWE) test P-value < 0.0001; 3) SNPs call rate < 95%. Finally, 658,288 SNPs in 1,452 subjects with an overall call rate of 99.68% were used for further analysis.
We used the Illumina Omni-Express platform to carry out the GWAS scan in FAMHES. There were 1,999 individuals (sample call rate >95%) included in the final statistical analysis. Based on quality control criteria, SNPs were excluded when P < 0.001 for the HWE test, MAF < 0.01, or genotype call rate < 95%. Finally, 709,211 SNPs were kept for further analysis.
In the validation stage, ten SNPs were selected based on the following criteria: 1) SNP with P < 1.0 × 10-5 for all GWAS samples; 2) when multiple SNPs showed a strong LD (r2 ≥ 0.8), SNPs previously reported in the literature were prior selected; 3) Clear genotyping clusters; 4) MAF ≥ 0.05. We used the iPLEX system (Sequenom) and/or the TaqMan assay (Applied Biosystems) [30, 31] to genotype the 10 SNPs. The primers and probes were available upon request.
We performed ungenotyped SNPs imputation using MACH 1.0 software (see URLs) via LD information from the HapMap phase II database (CHB + JPT as reference set, 2007-08_rel22, released 2007-03-02) in the DFTJ-cohort GWAS. Genotyped SNPs in the FAMHES GWAS were inferred using the IMPUTE program , and the reference panels used for imputation in the FAMHES study were HapMap rel.24, build 36, CHB + JPT. The Imputed SNPs with high genotype information content (proper info > 0.5 for IMPUTE and Rsq > 0.3 for MACH) were retained for the further association analysis. Finally, 2,468,160 SNPs were used for the further analysis.
Analysis for serum UA was performed on natural logarithmic (ln)-transformed values because of skewed distributions. Genome-wide association tests were performed using the additive model by linear regression analysis with adjustment for age, gender, BMI, cigarette smoking, and alcohol drinking implement using PLINK1.06 . The top two eigenvectors were also adjusted as covariates in the linear regression analysis. Population structure was evaluated using principal components analysis (PCA), as implemented by EIGENSTRAT software  and quantile-quantile (QQ) plot was generated by using R 2.11.1 (see URLs). Heterogeneity among the study populations was evaluated by the I2 statistic . The Manhattan plot of -log10 P, LD structures and haplotype block plots were generated by using Haploview (v4.1) . The association studies with the imputation data were performed using the ProbABEL software . The meta-analysis of the DFTJ-cohort GWAS data and the FAMHES GWAS data was performed using a fixed-effects meta-analysis with inverse variance weighted method using the metal software . The regional association plots were drawing using SNAP software .
The power of the present study was calculated using the Quanto software package . We used the mean of serum UA 292.5 μmol/L, MAF values obtained from the combined genotype dataset, and assumed an additive genetic model with α = 0.05 in two tail tests to calculate the statistical power.
A conditional analysis was done in linear regression model to examine the independence of the four associated SNPs in the combined data (partial r2 indicated the proportion of the serum UA variation explained by each SNP). We did the conditional analysis via including the four significant SNPs in the linear regression model adjusting for age, gender, BMI, cigarette smoking, and alcohol drinking in 10,282 individuals from the DFTJ-cohort. The SNPs with P value < 0.05 remained in the multivariate model were considered to be associated with the serum UA levels independently.
The combined dataset including the DFTJ-cohort GWAS data and the validation stage data totaling 10,282 individuals used in gene-environment interaction analysis was tested by introducing the interaction terms (SNP × gender, SNP × BMI, SNP × alcohol drinking, and SNP × cigarette smoking) into the model, adjusting for the covariates including age, gender, BMI, alcohol drinking, and cigarette smoking. The P value less than 0.05 for the interaction term was considered statistical significant. Considering multiple interaction tests were conducted, we further did the multiple test  for the gene-environmental interaction analysis. All statistical analyses were performed with SPSS (version 15.0; SPSS, Chicago, IL), and SAS version 9.2 (SAS Institute, Cary, NC, USA).
Genome-wide association of serum UA levels
The demographics of the participants are displayed in (Additional file 1: Table S1). The Q-Q plot revealed no inflation of type I error rate due to population stratification, with a genomic control inflation factor of 1.007 (Additional file 2: Figure S1). No heterogeneity was observed for the SNPs presented in the Table 1 between the DFTJ-cohort study and the FAMHES. As the Manhattan plot (Figure 1) and the regional association plots (Additional file 3: Figure S2) indicated, at the discovery stage the ABCG2 was significantly associated with the serum UA at a genome-wide significance level.
Totally we further investigated ten SNPs in the validation stage and four common variants were associated with serum UA levels at a genome-wide significant level of 5 × 10-8: SLC2A9 (rs11722228, combined P = 8.98 × 10-31) and ABCG2 (Combined P = 3.34 × 10-42 for rs2231142) (Table 1). Both loci of ABCG2 and SLC2A9 have been reported previously [20, 22, 27, 42]. If SNPs are in the same locus, we chosen r2 of SNPs less than 0.8 for further investigating (r2 = 0.168 for rs2231142 and rs4148152; r2 = 0.224 for rs2231142 and rs3114018; r2 = 0.669 for rs4148152 and rs3114018). (Additional file 3: Figure S2). However, we failed to replicate the top SNPs in the other loci of SEC22B, CENTG2, TET2, VEGFC, and TNFRSF11B in the validations stage (Additional file 4: Table S2).
The SNP rs11722228 in SLC2A9 and rs2231142 in ABCG2 were independently associated with serum UA levels and accounted for 1.03% and 1.09% of the serum UA variance respectively (Additional file 5: Table S3).
Ethnic differences in major genetic variants associated with serum UA levels
We compared the results in the present study with those in Japanese and Europeans. As Table 2 showes, the effect sizes of most of these loci showed consistent direction across the populations. ABCG2 was associated with serum UA levels across different populations. We found significant association of serum UA with SNP rs4148152 in ABCG2 in our study (combined P = 2.95 × 10-18). However, till now we did not find any reports of this SNP in European populations, which might be attributable to very low MAF (0.017) in Europeans. Furthermore, the SNP rs4148152 was in moderate LD with the SNP rs311408 in the present study (r2 = 0.669) but in very low LD in Europeans (r2 = 0.046; HapMap CEU). SNP rs12356193 in SLC16A9, rs10480300 in PRKAG2 and rs653178 in ATXN2 were significantly associated with serum UA in European populations but were monoallelic in Asians. Similar findings were found for loci of SLC2A9 (rs16890979 and rs734553), SLC22A11 (rs17300741), and recently reported new loci of TRIM46 (rs11264341), VEGFA (rs729761), BAZ1B (rs1178977), STC1 (rs1778674), A1CF (rs10821905), UBE2Q2 (rs1394125), and HLF (rs7224610)  which were associated with serum UA in Europeans, however, these associations were not replicated in Asians (Japanese and Chinese in the present study). This might be due to the very lower MAF in Asians. However, the relative small sample size in the present study limited us to detect these associations with enough power. In addition, we failed to replicate the SNPs rs742132 in LRRC16A, rs780094 in GCKR, rs17632159 in TMEM171, rs17050272 in IINHBB, and rs7188445 in MAF which were reported in European [20, 25], same as the findings in the Japanese population, albeit the MAFs of both of the SNPs are similar between Asians and Europeans, suggesting that there were genetic discrepancy on serum UA levels among different ethnic groups. The moderate effect size of these loci on serum UA levels and the limited sample size in the present study might be another potential explanation.
Interaction of SNPs with gender, BMI, cigarette smoking, and alcohol drinking
We further explored the interactions between the two independently-associated SNPs (rs11722228 in SLC2A9 and rs2231142 in ABCG2) and gender, BMI, cigarette smoking, and alcohol drinking on serum UA levels in 10,282 individuals from the DFTJ-cohort. Both SNPs rs11722228 and rs2231142 were nominally interacted with gender on serum UA levels (P = 0.04 and P = 0.02, respectively; Additional file 6: Table S4 and Figure 2). The minor allele (T) for rs11722228 in SLC2A9 has greater influence in elevating serum UA levels in females compared to males (beta = 0.051 in females vs. beta = 0.035 in males), similar with the previously reported results . In contrast, for SNP rs2231142 in ABCG2, the minor allele (T) had stronger effects on serum UA levels in males than those in females (beta = 0.037 in females vs. beta = 0.057 in males), consistent with the findings in the European populations and African Americans [20, 22]. However, the gene-gender interaction altered to null after the multiple test based on the FDR approach. Besides, the rs11722228 accounted for 1.33% and 0.98% and rs2231142 accounted for 0.93% and 2.25% of the total variance of the serum UA levels for females and males, respectively.
In addition, the SNP rs11722228 was also interacted with alcohol drinking and cigarette smoking (P for interaction = 0.016 and 0.035, respectively; Additional file 6: Table S4). Considering that the SNP rs11722228 had different effect size on serum UA levels in males and females and in the present study 95.1% smokers and 86.7% drinkers were males, we restricted our analysis in males. However, we failed to detect the interactions between rs11722228 and cigarette smoking and alcohol drinking on serum UA levels anymore (data not showed), suggesting that the SNP-smoking and the SNP-drinking interaction was driven by the gender differences in the serum UA concentrations.
In this two-stage GWAS, we replicated two previously reported loci (ABCG2 and SLC2A9) associated with serum UA levels and found significant gene-gender interactions in Chinese Han population. In addition, ethnic differences were observed between Asian and European populations.
SLC2A9 is located in chromosome 4p16-15.3 and encodes glucose transporter 9 (GLUT9) which can reabsorb UA in renal tubules . Several studies have reported the association between SLC2A9 and serum UA levels [20, 27, 43, 44]. Importantly, SLC2A9 is a transporter for both fructose and urate . Fructose intake could facilitate UA formation in liver via increasing purine breakdown. In addition, animal evidences indicated a causal relationship between fructose intake, serum UA, and metabolic syndrome [11, 46, 47].
ABCG2 is an UA exporter that mediates urate excretion in the kidney. Multiple evidences indicated that the common variants in ABCG2 could reduce the transport function and result in the hyperuricemia and gout [44, 48, 49]. In the present study, the missense SNP rs2231142 in ABCG2 showed the strongest association with serum UA level. This missense SNP could result in a glutamine-to-lysine amino acid substitution and the glutamine residue is highly conserved across species and the LD pattern differs in Chinese population and European population.
Comparisons of the SNPs of the association studies for serum UA in different populations are of great interest. The present study replicated two previously reported loci of SLC2A9 and ABCG2 associated with serum UA levels. SNP rs11722228 in SLC2A9 explained 1.03% variation compared to 1.33% in Japanese population  and rs2231142 in ABCG2 accounted for 1.09% of the variation of serum UA levels in our study compared to 1.20% in white individuals and 0.30% in African Americans . Previous studies reported that variants in SLC2A9 show the strongest effect on serum uric acid levels compared to ABCG2[20, 22]. However, in the present study we observed that the loci of ABCG2 and SLC2A9 show comparable amounts of explained variance (1.09% and 1.03% respectively). Table 2 showed that there was a little difference in MAFs in both variants rs2231142 located in ABCG2 and rs11722228 located in SLC2A9 (0.29 and 0.31) in Chinese which were different from that observed in Europeans (0.11 and 0.50 respectively).
Owing to the SNPs of rs12356193 in SLC16A9, rs10480300 in PRKAG2, and rs653178 in ATXN2 are monoallelic in Asians, we failed to replicate the SNPs of rs12356193, rs10480300 and rs653178 those were identified in Europeans. In addition, we found notable differences in MAFs for rs16890979 in SLC2A9, rs17300741 in SLC22A11, rs10821906 in AICF, rs1394125 in UBE2Q2, rs7193778 in NFAT5 and rs7224610 in HLF (0.02 , 0.07, 0.04, 0.008, 0.05, 0.13 respectively for Chinese population and 0.29, 0.49, 0.18, 0.34, 0.14, 0.42 respectively for European populations) Besides, both rs16890979 (r2 = 0.005) and rs734553 (r2 = 0.005) in the gene SLC2A9 were in very low LD with the SNP rs11722228 in Asians compared to that in European populations (r2 = 0.202 for rs16890979; r2 = 0.193 for rs734553). The very low MAF and the difference in the LD structure might partly explain the discrepancy of the associations between the Europeans and Asians. In addition, relative small sample size in the present study and the moderate effect size of the loci on serum UA levels also limited us to have enough power to detect these associations.
Serum UA levels are lower in females than that in males. More importantly, gene-gender interactions were observed for the two independent SNPs of rs11722228 and rs2231142. The minor allele (T) for rs11722228 has greater influences in elevating serum UA levels in females compared to males, consistent with the previous study . For rs2231142, the TT genotype was related to higher UA as compared with GG and GT genotype in both males and females. This might be due to that the rs2231142 T allele was related to the reduced ability to excrete UA . In addition, the minor allele T allele of rs2231142 has greater effects on the serum UA levels in males than those in females (P for interaction = 0.02), in consistent with the findings from the Europeans . This gender difference might be due to the specific physiological characteristics in females. It is also reported that the estrogens might increase the renal clearance of serum UA , however, the mechanism underlying the gene-gender interaction remains to be further elucidated.
To our best of knowledge, this is the first GWA study on serum UA levels in Chinese population. Because of the relative small sample size in the discovery stage, we might have limited power to detect the associations of the SNPs with small effect size and/or low MAF. However, in the present study we identified two reported loci (SLC2A9 and ABCG2) associated with serum UA levels, suggesting that our study was capable of identifying significant loci associated with serum UA levels. In addition, the combined data in our study has more than 90% statistic power to detect the interaction between SNP rs11722228 and gender; 70.4% statistical power to detect the interaction between SNP rs2231142 and gender on serum UA levels. Furthermore, our study confirmed the gene-gender interaction on serum UA levels and observed that the ABCG2 and SLC2A9 functioned differently in males and females across different populations .
Our study replicated two loci of SLC2A9 and ABCG2 associated with serum UA levels in a Chinese Han population. Heterogeneity was observed among different ethnic populations. In addition, we found both the two loci interacted with gender on serum UA levels. These two loci had different effects on serum UA levels in males and females. Further studies are needed to validate our findings and investigate the underlying mechanisms of the gene-gender interaction.
Johnson RJ, Rideout BA: Uric acid and diet–insights into the epidemic of cardiovascular disease. N Engl J Med. 2004, 350 (11): 1071-1073. 10.1056/NEJMp048015.
Choi HK, Mount DB, Reginato AM: Pathogenesis of gout. Ann Intern Med. 2005, 143 (7): 499-10.7326/0003-4819-143-7-200510040-00009.
Oda M, et al: Loss of urate oxidase activity in hominoids and its evolutionary implications. Mol Biol Evol. 2002, 19 (5): 640-653. 10.1093/oxfordjournals.molbev.a004123.
Taniguchi A, Kamatani N: Control of renal uric acid excretion and gout. Curr Opin Rheumatol. 2008, 20 (2): 192-197. 10.1097/BOR.0b013e3282f33f87.
Ames BN, et al: Uric acid provides an antioxidant defense in humans against oxidant- and radical-caused aging and cancer: a hypothesis. Proc Natl Acad Sci U S A. 1981, 78 (11): 6858-6862. 10.1073/pnas.78.11.6858.
Sundstrom J, et al: Relations of serum uric acid to longitudinal blood pressure tracking and hypertension incidence. Hypertension. 2005, 45 (1): 28-33. 10.1161/01.HYP.0000150784.92944.9a.
Nakanishi N, et al: Serum uric acid and risk for development of hypertension and impaired fasting glucose or Type II diabetes in Japanese male office workers. Eur J Epidemiol. 2003, 18 (6): 523-530.
Fang J, Alderman MH: Serum uric acid and cardiovascular mortality the NHANES I epidemiologic follow-up study, 1971–1992. National Health and Nutrition Examination Survey. JAMA. 2000, 283 (18): 2404-2410. 10.1001/jama.283.18.2404.
Richette P, Bardin T: Gout. Lancet. 2010, 375 (9711): 318-328. 10.1016/S0140-6736(09)60883-7.
Dehghan A, et al: High serum uric acid as a novel risk factor for type 2 diabetes. Diabetes Care. 2008, 31 (2): 361-362.
Heinig M, Johnson RJ: Role of uric acid in hypertension, renal disease, and metabolic syndrome. Cleve Clin J Med. 2006, 73 (12): 1059-1064. 10.3949/ccjm.73.12.1059.
Burack RC, Keller JB, Higgins MW: Cardiovascular risk factors and obesity: are baseline levels of blood pressure, glucose, cholesterol and uric acid elevated prior to weight gain?. J Chronic Dis. 1985, 38 (10): 865-872. 10.1016/0021-9681(85)90111-0.
Cirillo P, et al: Uric acid, the metabolic syndrome, and renal disease. J Am Soc Nephrol. 2006, 17 (12 Suppl 3): S165-S168.
Brandstatter A, 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-1667. 10.2337/dc08-0349.
Johnson RJ, et al: Uric acid, evolution and primitive cultures. Semin Nephrol. 2005, 25 (1): 3-8. 10.1016/j.semnephrol.2004.09.002.
Haj Mouhamed D, et al: Effect of cigarette smoking on plasma uric acid concentrations. Environ Health Prev Med. 2011, 16 (5): 307-312. 10.1007/s12199-010-0198-2.
White JS: Comment on: new insights into the epidemiology of gout. Rheumatology. 2010, 49 (3): 613-614. 10.1093/rheumatology/kep349.
Wilk JB, et al: Segregation analysis of serum uric acid in the NHLBI family heart study. Hum Genet. 2000, 106 (3): 355-359. 10.1007/s004390051050.
Yang Q, et al: Genome-wide search for genes affecting serum uric acid levels: the Framingham Heart Study. Metabolism. 2005, 54 (11): 1435-1441. 10.1016/j.metabol.2005.05.007.
Kolz M, et al: Meta-analysis of 28,141 individuals identifies common variants within five new loci that influence uric acid concentrations. PLoS Genet. 2009, 5 (6): e1000504-10.1371/journal.pgen.1000504.
Karns R, et al: Genome-wide association of serum uric acid concentration: replication of sequence variants in an island population of the Adriatic coast of Croatia. Ann Hum Genet. 2012, 76 (2): 121-127. 10.1111/j.1469-1809.2011.00698.x.
Dehghan A, 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-1961. 10.1016/S0140-6736(08)61343-4.
Doring A, et al: SLC2A9 influences uric acid concentrations with pronounced sex-specific effects. Nat Genet. 2008, 40 (4): 430-436. 10.1038/ng.107.
Wallace C, 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-149. 10.1016/j.ajhg.2007.11.001.
Kottgen A, et al: Genome-wide association analyses identify 18 new loci associated with serum urate concentrations. Nat Genet. 2013, 45 (2): 145-154.
Okada Y, et al: Meta-analysis identifies multiple loci associated with kidney function-related traits in east Asian populations. Nat Genet. 2012, 44 (8): 904-909. 10.1038/ng.2352.
Kamatani Y, et al: Genome-wide association study of hematological and biochemical traits in a Japanese population. Nat Genet. 2010, 42 (3): 210-215. 10.1038/ng.531.
Wang F, et al: Cohort profile: The Dongfeng-Tongji cohort study of retired workers. Int J Epidemiol. 2012, 42 (3): 731-740.
Tan A, et al: Low serum osteocalcin level is a potential marker for metabolic syndrome: results from a Chinese male population survey. Metabolism. 2011, 60 (8): 1186-1192. 10.1016/j.metabol.2011.01.002.
McCarroll SA, et al: Integrated detection and population-genetic analysis of SNPs and copy number variation. Nat Genet. 2008, 40 (10): 1166-1174. 10.1038/ng.238.
Korn JM, et al: Integrated genotype calling and association analysis of SNPs, common copy number polymorphisms and rare CNVs. Nat Genet. 2008, 40 (10): 1253-1260. 10.1038/ng.237.
Marchini J, et al: A new multipoint method for genome-wide association studies by imputation of genotypes. Nat Genet. 2007, 39 (7): 906-913. 10.1038/ng2088.
Purcell S, et al: PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007, 81 (3): 559-575. 10.1086/519795.
Price AL, et al: Principal components analysis corrects for stratification in genome-wide association studies. Nat Genet. 2006, 38 (8): 904-909. 10.1038/ng1847.
Higgins JP, et al: Measuring inconsistency in meta-analyses. BMJ. 2003, 327 (7414): 557-560. 10.1136/bmj.327.7414.557.
Barrett JC, et al: Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005, 21 (2): 263-265. 10.1093/bioinformatics/bth457.
Aulchenko YS, Struchalin MV, van Duijn CM: ProbABEL package for genome-wide association analysis of imputed data. BMC Bioinforma. 2010, 11 (1): 134-10.1186/1471-2105-11-134.
Kathiresan S, et al: Polymorphisms associated with cholesterol and risk of cardiovascular events. N Engl J Med. 2008, 358 (12): 1240-1249. 10.1056/NEJMoa0706728.
Johnson AD, et al: SNAP: a web-based tool for identification and annotation of proxy SNPs using HapMap. Bioinformatics. 2008, 24 (24): 2938-2939. 10.1093/bioinformatics/btn564.
Gauderman WJ: Sample size requirements for association studies of gene-gene interaction. Am J Epidemiol. 2002, 155 (5): 478-484. 10.1093/aje/155.5.478.
Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B Methodol. 1995, 57 (1): 289-300.
Guan M, et al: Association of an intronic SNP of SLC2A9 gene with serum uric acid levels in the Chinese male Han population by high-resolution melting method. Clin Rheumatol. 2011, 30 (1): 29-35. 10.1007/s10067-010-1597-x.
Li S, et al: The GLUT9 gene is associated with serum uric acid levels in Sardinia and Chianti cohorts. PLoS Genet. 2007, 3 (11): e194-10.1371/journal.pgen.0030194.
Yamagishi K, et al: The rs2231142 variant of the ABCG2 gene is associated with uric acid levels and gout among Japanese people. Rheumatology. 2010, 49 (8): 1461-1465. 10.1093/rheumatology/keq096.
Vitart V, et al: SLC2A9 is a newly identified urate transporter influencing serum urate concentration, urate excretion and gout. Nat Genet. 2008, 40 (4): 437-442. 10.1038/ng.106.
Rutledge AC, Adeli K: Fructose and the metabolic syndrome: pathophysiology and molecular mechanisms. Nutrition reviews. 2008, 65 (s1): S13-S23.
Caulfield MJ, et al: SLC2A9 is a high-capacity urate transporter in humans. PLoS Med. 2008, 5 (10): e197-10.1371/journal.pmed.0050197.
Kusuhara H, Sugiyama Y: ATP-binding cassette, subfamily G (ABCG family). Pflugers Arch. 2007, 453 (5): 735-744. 10.1007/s00424-006-0134-x.
Matsuo H, et al: Identification of ABCG2 dysfunction as a major factor contributing to gout. Nucleosides Nucleotides Nucleic Acids. 2011, 30 (12): 1098-1104. 10.1080/15257770.2011.627902.
Woodward OM, et al: Identification of a urate transporter, ABCG2, with a common functional polymorphism causing gout. Proc Natl Acad Sci U S A. 2009, 106 (25): 10338-10342. 10.1073/pnas.0901249106.
Anton FM, et al: Sex differences in uric acid metabolism in adults: evidence for a lack of influence of estradiol-17 beta (E2) on the renal handling of urate. Metabolism. 1986, 35 (4): 343-348. 10.1016/0026-0495(86)90152-6.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1755-8794/7/10/prepub
This study was funded by grants from the National Basic Research Program grant (2011CB503806), the Program of Introducing Talents of Discipline to Universities and Yangze Scholarship to T.W., the General Program of National Natural Science Foundation of China (30945204, 30360124, 30260110, 81172751), the Guangxi Provincial Department of Finance and Education (2009GJCJ150), intramural funding from Fudan-VARI (Van Andel Institute, USA) Center for Genetic Epidemiology, Intramural Funding from Fudan University Institute of Urology to Z.M. and the Program for New Century Excellent Talents in University to M.H.. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. We thank the staff and participants in the DFTJ-cohort study, FAMHES study, and thank all the volunteers for assisting in collecting the data and samples.
R statistical environment, http://cran.r-project.org/;
MACH 1.0, http://www.sph.umich.edu/csg/abecasis/mach/;
International HapMap Project, http://www.hapmap.org/index.html.
The authors declare that they have no competing interests.
Projects conception: WTC, MZN, HMA, FBH. Study design: WTC, HMA, YBY, LDX, WC, YHD, ZXM, YJ, GH. Sample processing and database establishment: YBY, YXB, HYF, ZL, DQF, HSL, ZJ, MXW, LMJ, LDF, QGK, GL, HD, ZX, FYY, and WT. Genotyping: HMA, YBY, And GLX. Imputation: HMA, YBY, GLX. Data analysis: YBY, DXY, and LJ. Interpretation of the data: HMA, MZN, YBY. Drafted the manuscript: YBY and MZN. Revised manuscript: WTC, HMA, MZN. All authors have read and approved of the final version of this manuscript.
Binyao Yang, Zengnan Mo contributed equally to this work.