Skip to main content

Polygenic risk score trend and new variants on chromosome 1 are associated with male gout in genome-wide association study



Gout is a highly hereditary disease, but not all those carrying well-known risk variants have developing gout attack even in hyperuricemia status. We performed a genome-wide association study (GWAS) and polygenic risk score (PRS) analysis to illustrate the new genetic architectures of gout and asymptomatic hyperuricemia (AH).


GWAS was performed to identify variants associated with gout/AH compared with normouricemia. The participants were males, enrolled from the Taiwan Biobank and China Medical University, and divided into discovery (n=39,594) and replication (n=891) cohorts for GWAS. For PRS analysis, the discovery cohort was grouped as base (n=21,814) and target (n=17,780) cohorts, and the score was estimated by grouping the polymorphisms into protective or not for the phenotypes in the base cohort.


The genes ABCG2 and SLC2A9 were found as the major genetic factors governing gouty and AH, and even in those carrying the rs2231142 (ABCG2) wild-genotype. Surprisingly, variants on chromosome 1, such as rs7546668 (DNAJC16), rs10927807 (AGMAT), rs9286836 (NUDT17), rs4971100 (TRIM46), rs4072037 (MUC1), and rs2974935 (MTX1), showed significant associations with gout in both discovery and replication cohorts (all p-values < 1e−8). Concerning the PRS, the rates of gout and AH increased with increased quartile PRS in those SNPs having risk effects on the phenotypes; on the contrary, gout/AH rates decreased with increased quartile PRS in those protective SNPs.


We found new variants on chromosome 1 significantly relating to gout, and PRS predicts the risk of developing gout/AH more robustly based on the SNPs’ effect types on the trait.


Gout is known as an inflammatory and complex disease due to the interactions between causal genes and environmental factors. Clustering of hyperuricemia and gout has been known to have familial links, suggesting a high hereditary component to these traits. Mutations have led the genetic variation in modern human populations. Over the past decade, more than 28 SNPs found by genome-wide association study (GWAS) across the world have been shown to relate to gout/hyperuricemia [1,2,3,4,5], but these accounted for only 7% of the variance of in serum urate concentrations [6]. Gouty attack involved both inflammatory and hyperuricemia mechanisms. To identify underlying genetic alterations involved in the two mechanisms, genome-wide association studies (GWASs) have been conducted in many populations, and genetic loci covering ABCG2 and SLC2A9 have been found most repeatedly and significantly [7,8,9,10,11].

Most of the genes identified in current GWASs seem plausible contributors to serum uric acid levels because they encode for proteins that are responsible for uric acid excretion in the intestines or kidney, including the genes ABCG2 and SLC2A9 [12, 13]. However, except for rs2231142 in ABCG2, which is a missense variant, most genetic variants that have been reported to be significant so far are located in intronic or intergenic regions, and causative variants have not been identified for most loci.

Not all those carrying genetic variants of the aforementioned genes would developing gout disease. This may involve unknown genes or the mechanism of genetic interaction. However, gout attack is a complex disease, and an estimation of polygenic risk score (PRS) is one way by which people can measure their risk of developing clinical diseases [14, 15]. PRS based on the total number of changes related to the disease is an estimate of an individual’s genetic liability to a trait, calculated according to their genotype profile and relevant genome-wide association study (GWAS) data. Here, we performed a GWAS with a large sample of male participants to reveal the genetic contributions and estimation of PRS on the phenotypes of gout and asymptomatic hyperuricemia (AH).


Study population

We aimed to explore the genetic variants associated with the development of gout and AH and enrolled 39,594 male participants as a discovery cohort. A second cohort with whole genome sequence data was included as a replication cohort, which was composed of 891 male participants. All participants were aged between 30 and 70 years old, and they were selected from Taiwan Biobank and China Medical University in Taiwan. We grouped the participants of the discovery cohort into base cohort (n=21,814) and target cohort (n=17,780) for PRS analysis (Supplementary Table 1). PRSice-2 software (version 2.3.3 for R) was applied for the estimation of PRS with applying the odds ratios (OR) estimated by GWAS data in the base cohort [16], which was adjusted by age, and then we estimated the individual PRS of gout/AH phenotype for the target cohort. The threshold for p was set to 5e−5 in the PRS analysis, which can ensure to retrieve enough significant variants for estimating individual PRS. In order to explore the genetic component without confounded by rs2231142 (ABCG2 gene) which was the most significant variant with gout in Taiwan [17], we further analyzed the associations between variants and gout/AH under a selection of those carrying rs2231142 wild-type (GG). TWB is a nationwide database for research, which combines genomic profiles with lifestyle patterns from people in Taiwan to explore the relationships between diseases and genetics. Each participant underwent biochemical testing (with blood samples) and physical examination; all of them also signed informed consent for genotyping analysis. The phenotype of gout disease was self-reported and previous studies utilizing a similar method of self-reported gout suggested that self-reporting of physician-diagnosed gout retained good reliability and sensitivity [18, 19]. Those with serum uric acid more than or equal to 7 mg/dl without gout history were defined as AH cases.

The ethics committee of the China Medical University Hospital Institutional Review Board in Taiwan (CMUH108-REC1-091) approved this project. Both the Declaration of Helsinki and the Good Clinical Practice Guidelines were followed and informed consent was written by all participants.

Genotyping and quality control

The participants of the discovery cohort were genotyped using Affymetrix Axiom genotyping array (chip: TWB2), including 680K SNPs. Quality control (QC) was applied to leave out those SNPs with low call rate (< 95%), minor allele frequency (MAF) less than 0.01, and deviation from the Hardy-Weinberg equilibrium (p < 0.0001; library HardyWeinberg of R program). A total of 401,037 SNPs passed the quality control for GWAS. The participants of replicated cohort underwent whole genome sequencing on the Illumina Hiseq platform. For the Illumina Hiseq platform, variant calling was performed with Isaac Variant Caller version 2.0.17, Grouper version 1.4.2, and CNVseg version 2.2.4. Alleles were annotated with ANNOVAR [20]. Some of the significant polymorphisms found by the discovery cohort were extracted from the variants of variant call format (VCF). All the genetic positions displayed in this study were GRCh38 format.

Statistical analysis

The clinical characteristics of subjects with and without gout and AH were compared by applying a chi-square test for categorical variables and t-test or ANOVA for continuous variables. The p-values estimated by chi-square test (chisq.test function of R program) for the association between genetic variants located on each chromosome for the gout/AH were presented in Manhattan plots (library qqman) by using the R-program provided by Turner [21]. We conducted a logistic regression model for estimating the odds ratios (ORs) as well as 95% confidence intervals (95% CI) for the associations between genetic variants and gout/AH after adjustment of age (library aod of R program). The Locuszoom plot was employed to visualize the regional strength and associations between loci and gout/hyperuricemia related to local linkage-disequilibrium (LD), and genomic position (locuszoom v1.4) [22]. A heatmap of paired LD for significant SNPs related to gout/AH was displayed for the variants located in the same chromosome (LDheatmap library) [23]. We used PLINK (v1.9), PERL (v5.16), and R (v3.6) programs (CentOS v7.0) to mine raw data, estimate p-values, and draw the figures. For the GWAS, after Bonferroni correction for multiple testing, the significance was determined at p < 1 × 10−8. For other analyses, significance was determined at p < 0.05.


A total of 39,594 male participants were included as the discovery cohort in this study, which was composed of 5857 gout patients, 12,382 AH subjects, and 21,355 normouricemia participants. The mean age of the three groups were 54.00 (± 10.65), 51.27 (± 11.30), and 51.72 (± 11.04) years old, and the mean uric acid levels were 8.79 (± 2.26 mg/dl), 8.09 (± 1.12 mg/dl), and 5.72 (±0.86 mg/dl), respectively (Supplementary Table 1).

To explore the genetic associations of SNPs with the phenotypes, we illustrated the Manhattan plots which showed all p-values of SNPs related to gout/AH across 23 pairs of chromosomes, and the result showed a total of 120 variants significantly associated with gout, located in a total of 35 genes (Supplementary Table 2; Fig. 1A). The genes which contain more than or equal to two significant SNPs include genes ABCG2, SLC2A9, SLC17A1, ZNF518B, SLC17A3, GCKR, PKD2, MEPE, PDZK1, HECTD4, CD160, DNAJC16, TRIM46, CCDC63, ALDH2, and CUX2. Concerning AH, a total of 57 SNPs located in 13 genes were significantly associated with hyperuricemia compared to normal subjects (Supplementary Table 3; Fig. 1B). The genes of the variants were located in ABCG2, SLC2A9, ZNF518B, PKD2, MEPE, ALDH2, HECTD4, BRAP, SPP1, MAPKAPK5, CLNK, ACAD10, and BCAS3. A total of 27 SNPs were significantly related to gout compared to AH, located in genes ABCG2, PKD2, and SPP1 (Supplementary Table 4; Fig. 1C). Interestingly, ABCG2, PKD2, and SPP1 genes were found to be significantly related to gout/AH in all three comparisons (gout vs normal; gout vs AH; AH vs normal).

Fig. 1
figure 1

The Manhattan plots reveal the p-values of variants related to gout and asymptomatic hyperuricemia (AH). The p-values were estimated for developing gout compared to normal (A). When compared to normal, the p-values of variants related to AH are illustrated (B), and compared to AH, the values of variants related to gout are illustrated (C). The red horizontal line is denoted as a cut-off for the p-value of significant difference between the phenotypes and variants by 1e−8

Regarding LD between those significant variants, we applied a heatmap of paired variants to show the significant variants of LD degree for gout/AH (Fig. 2). Three LD blocks existed in chromosome 1 for the variants relating to gout (Fig. 2A), which were composed of genes DNAJC16, AGMAT (first block), PDZK1, CD160, NUDT17 (second block), TRIM46, MUC1, MTX1, and ASH1L (third block). Surprisingly, the results showed that three variants of DNAJC16 and one variant of AGMAT were significantly associated with gout, which had not been found before, namely, rs7546668 (DNAJC16), rs12124078 (DNAJC16), rs7515244 (DNAJC16), and rs10927807 (AGMAT). Furthermore, the four polymorphisms existed in a high LD degree (Supplementary Fig. 1). Two major LD blocks relating to gout were observed in chromosome 4 which is located on the p arm and q arm (Fig. 2B, C); the major causal genes are SLC2A9 and ABCG2, respectively. Moreover, all the variants in chromosomes 6 (Fig. 2D) and 12 (Fig. 2E) existed in high LD in each chromosome. Concerning AH, the result also revealed two major blocks and one block associated with hyperuricemia in chromosomes 4 and 12, respectively (Supplementary Fig. 2).

Fig. 2
figure 2

The paired linkage disequilibrium (LD) of the variants related to gout. The LD illustrates the variants of developing gout compared to normal in chromosome 1 (A), chromosome 4 on p arm (B), chromosome 4 on q arm (C), chromosome 6 (D), and chromosome 12 (E). The green line blocks the gene, and the red circle indicates it is the most significant variant in its gene

To explore causal variants without the confounding effect caused by the variant rs2231142, we further estimated the associations between the SNPs and the gout/AH phenotypes under the selection limitation of those carrying wild genotypes (GG) of rs2231142 (n=18,177). We applied Manhattan plots to show the significant variants across all chromosomes (Supplementary Fig. 3). The results showed 12 variants were significantly related to gout compared to normal, four variants significantly related to AH compared to normal, and none was found to be significantly related to gout compared to AH (Supplementary Table 5). The involved genes were ABCG2, SLC2A9, PKD2, ZNF518B, ALDH2, HECTD4, and MAPKAPK5 in the phenotype of gout; and only gene SLC2A9 was found in the AH phenotype. Furthermore, we estimated the OR of developing the phenotypes of the abovementioned 16 variants which were significantly related to the phenotypes under the wild-type selection of rs2231142 (Table 1). The results revealed that the variants exhibited protective effects against developing gout, namely rs3733589 (SLC2A9), rs3775948 (SLC2A9), rs1014290 (SLC2A9), rs3109823 (ABCG2), rs2622604 (ABCG2), rs6532055 (ABCG2), rs72554040 (ABCG2), rs671 (ALDH2), rs78069066 (MAPKAPK5), and rs77768175 (HECTD4). Polymorphisms rs17407555 (ZNF518B) and rs2725231 (PKD2) showed risk effects, however. The polymorphisms located in chromosome 12 showed protective effects, and all existed in a high LD (Supplementary Fig. 4). Concerning the development of AH, the polymorphism rs10805346 (SLC2A9) acted as a risk factor, but rs3733589 (SLC2A9), rs3775948 (SLC2A9), and rs1014290 (SLC2A9) acted as protective variants for the phenotype (Table 1), and these existed in a weak LD (r2 >0.22; Supplementary Fig. 5).

Table 1 The OR and 95% CI of variants associating with gout/hyperuricemia in those carrying rs2231142 (ABCG2) wild type (GG)

The replication cohort revealed significant variants of chromosome 1 which were found in the discovery cohort. The results showed genotype TT of rs10927807 located in gene AGMAT revealed 4.34 times risk for developing gout (OR=4.34, 95% CI =2.68–7.02, p=6.62e−10; Table 2). Moreover, rs9286836 (NUDT17) and rs4971100 (TRIM46), rs4072037 (MUC1) and rs2974935 (MTX1) also revealed significant associations with gout (all p-values <1e−8).

Table 2 The associations between gout and significant variants of chromosome 1 in replication cohort

A summary circle figure illustrates the aforementioned six phenotypes to exhibit the genes with significant variants (Fig. 3). The results showed ABCG2 gene was involved in four phenotypes, except for those developing hyperuricemia and gout (compared to AH) under the selection of rs2231142 wild type. Notably, the SLC2A9 exhibited effects on gout and AH phenotypes under a selection of those carrying rs2231142 wild-genotype, but the effect was not observed in developing gout compared to AH.

Fig. 3
figure 3

The circular map illustrates the genes and locations of significant variants related to the phenotypes. Gout: the comparison between gout and normal; Gout-WD: the same comparison under a selection of those carrying rs2231142 wild-genotype GG. HUA: the comparison between hyperuricemia and normal; HUA-WD: the same comparison under a selection of those carrying rs2231142 wild-genotype GG. Gout-HUA: the comparison between gout and hyperuricemia; Gout-HUA-WD: the same comparison under a selection of those carrying rs2231142 wild-genotype GG

Concerning the associations between PRS and gout/AH phenotypes, compared to the controls. More than 13 significant variants were retrieved from each PRS analysis. The rates of gout and AH increased with increased quartile PRS in those variants showing risk effects on the phenotypes (all p-trend < 1e−11; Fig. 4A, B). On the contrary, for those variants with protective effects on gout/AH, the trend decreased with increased quartile PRS (all p-trend < 1e−4; Fig. 4C, D). Interestingly, these rates of gout/AH in Q4 of PRS also showed significant association with the rates in other three quartiles regardless of whether the variants had risk effects or not (all p-values < 0.01).

Fig. 4
figure 4

The rates of gout and asymptomatic hyperuricemia (AH) compared to those with normouricemia among different quartile polygenic risk scores (PRS). The rates of gout and AH increased with increased quartile PRS in those variants with risk effects (odds ratio (OR) >=1; red color) of developing the phenotypes (A, B, respectively). Whereas the rates of gout and AH decreased with increased quartile PRS in those variants with protective effects (OR <1; blue color) of developing the phenotypes (C, D, respectively). The rates of gout/AH were calculated only in the participants for estimating PRS, but not in the whole population. HUA, hyperuricemia


We performed a large GWAS to explore the genetic variants and PRS relating to gout and AH in males, and we found that the variants located in genes ABCG2 and SLC2A9 were the major genetic factors governing gouty attack and hyperuricemia. Surprisingly, variants on chromosome 1 located in genes DNAJC16, AGMAT, NUDT17, TRIM46, MUC1 and MTX1 showed significant associations with gout, which finding had not been reported before. Genes AGMAT and DNAJC16 exist a high LD. Recently, Yamamoto et al. revealed DNACJ16 to be involved in the autophagy process delivering the engulfed substrates into lysosomes for degradation; moreover, the same researchers also reviewed recent advances in the relationship between autophagosome formation and endoplasmic reticulum [24], and they showed that DNAJC16 overexpression extended the size of the autophagosome through its DnaJ and TRX domains and that DNAJC16 ablation resulted in a defect in engulfing larger targets [25]. The phagocytosis of monosodium urate crystals triggering an inflammatory cellular state in synoviocytes suggests a possible mechanism in the pathogenesis of crystal-induced arthritis. [26] In a previous study, Chen et al. explored a cell model to show significantly higher IL-8 release from endothelial cells combined with ABCG2 knockdown [17]. Liao et al. found that monosodium urate (MSU) alone induced cGKII and TLR2 expression with increased phagocytosis activity, which suggested that cGKII knockdown significantly inhibited this MSU-induced cGKII-TLR2-phagocytosis axis [27]. Herein, we suggest that the new gene DNAJC16 associated with gout disease may be involved in the mechanism of urate phagocytosis.

Arginase is found ubiquitously in mammals, and there are two distinct isoforms. Isoform-I is expressed predominantly in the liver and plays an important role in the urea cycle; its reaction product, ornithine, is a precursor for the synthesis of the amino acid proline and glutamic acid for cellular replication [28]. Isoform-II is widely expressed in mitochondria in the extrahepatic tissues and is believed to play an important role in endothelial dysfunction [29]. Interestingly, Arginase is also involved in the control of nitric oxide synthesis, and its function is associated with hypertension, inflammation and diabetes [29, 30]. AGMAT encodes agmatinase, which is a metallohydrolase involved in the hydrolysis of agmatine to produce urea and putrescine. This enzyme plays a critical role in polyamine biosynthesis and the regulation of agmatine levels. Early evidence showed LPS induced a dose-dependent stimulation of agmatinase, and exposure to IL-10 and TGF-beta caused a reduction in agmatinase and arginine decarboxylase (ADC), whereas IL-4 was ineffective on ADC, but reversed the LPS-induced increase of agmatinase [31], which suggests a regulatory role of agmatine during inflammation. Sastre et al. discovered that both ADC and agmatinase are constitutively expressed in macrophages [31], with activities comparable to those reported in the brain or kidney [32, 33]. In addition to its role as a component of proteins, arginine is a precursor for the synthesis of other molecules including nitric oxide (NO), citrulline, agmatine, urea, and ornithine. Two well-described pathways of L-arginine metabolism in inflammation are the conversion of arginine to NO and citrulline by NOS and the breakdown of arginine to urea and ornithine by arginase. These pathways are temporally regulated in inflammatory models including wound healing and glomerulonephritis [34,35,36]. Although we found a higher association between gout and AGMAT than between gout and DNACJ16 in the replication cohort, from the above literature review, both genes may possibly cause gout.

Even under selection limitation of those carrying rs2231142 wild genotype in this study, the major risk genes for developing gout were located in chromosomes 4 and 12, including region of chromosome 4 of SLC2A9, ZNF518B, PKD2, ABCG2, and the region of chromosome 12 of genes ALDH2, MAPKAPK5, and HECTD4. This illustrates that these genes show effects on gout independently of rs2231142 (ABCG2). Moreover, the variants of SLC2A9 show protective and opposing effects on hyperuricemia in those carrying the rs2231142 wild genotype. The SLC2A9 gene was shown to associate with lower uric acid levels and hypouricemia. The effect of contributing to uric acid was more pronounced in those without carrying the rs2231142 variant.

SLC2A9 encodes a transporter protein belonging to class II of the facilitative glucose transporter family (GLUT). GLUT9 is expressed in the proximal tubules [37], and it mediates sodium-independent hexose uptake into target cells by facilitated diffusion [38]. Fructose is potentially a substrate of GLUT9, and it was identified as a marker for uric acid concentrations [39]. Compared with other carbohydrates, fructose is associated with higher levels of uric acid, triglycerides, and cholesterol [40]. Sex-specific differences in the effects of SLC2A9 have also been observed in previous studies. Doring et al. found that the variance of uric acid levels explained by SLC2A9 was five-fold higher in females than in males, and the percentage of variance accounted for by SLC2A9 expression levels was 3.5% in men and 15% in women [41].

ABCG2 is a urate transporter protein that regulates uric acid excretion in the kidney, intestine, and liver [42,43,44]. Woodward et al. discovered that ABCG2 is located at the apical border membrane of renal proximal tubule cells and plays a role in uric acid secretion [45]. Furthermore, the Q141K mutation encoded by rs2231142 resulted in 53% lower urate transport rates compared to the wild type. Moreover, another pathway for uric acid elimination was recently discovered in the intestine. In a nephrectomized rat model, ABCG2 was significantly upregulated in the ileum compared to wild-type rats [13]. Moreover, Merino et al. found that expression of the hepatic protein BCRP1, a substrate of ABCG2, is higher in male mice than in female mice; they also found that hepatic expression of human BCRP/ABCG2 was higher in men [46].

Although Hunter and Mars et al. suggested a PRS model which, by adding up the number of risk variants, could be applied to estimate the genetic risk of gene-susceptible diseases [47, 48]. Our data revealed that the rates of gout and AH increased with quartile of PRS in the variants of showing risk effects on the phenotypes, while the rates decreased with increased quartile PRS in the protective variants.


This study performed a GWAS in a large sample size and found genes DNAJC16, AGMAT, NUDT17, TRIM46 MUC1, and MTX1 on chromosome 1 relating to gout disease. Gene SLC2A9 is the major independent gene of controlling uric acid and governing gout attack in those carrying rs2231142 wild-type; and PRS provides a robust insight into the genetic component of gout and AH in the selection of variant-effect types on the traits.

Availability of data and materials

Data is not publicly available. Raw data used for this study is available from the Taiwan Biobank and China Medical University Hospital.


95% CI:

95% confidence intervals


ATP-binding cassette, subfamily G, member 2


Acyl-CoA dehydrogenase family, member 10


Arginine decarboxylase




Asymptomatic hyperuricemia


Aldehyde dehydrogenase 2 family


Ash1-like histone lysine methyltransferase


Breast carcinoma amplified sequence 3


Breast cancer resistance protein 1


BRCA1-associated protein


Coiled-coil domain-containing protein 63


Protein kinase, cGMP-dependent, type II


Toll-like receptor 2


Cytokine-dependent hematopoietic cell linker


CUT-like homeobox 2


DnaJ heat shock protein family member C16


Glucokinase regulator


Glucose transporter family


Genome-wide association study


HECT domain E3 ubiquitin protein ligase 4






Minor allele frequency


MAPK activated protein kinase 5


Matrix extracellular phosphoglycoprotein


Monosodium urate


Metaxin 1


Mucin 1


Nitric oxide


Nitric oxide synthase


Nudix hydrolase 17


Odds ratios


PDZ domain-containing 1


Polycystic 2


Polygenic risk score


Fourth quartile


Solute carrier family 17 member 1


Solute carrier family 17 member 3


Solute carrier family 2 member 9


Single Nucleotide polymorphism


Secreted phosphoprotein 1


Transforming growth factor


Tripartite motif-containing 46


Zinc finger protein 518B


  1. Nakayama A, Nakaoka H, Yamamoto K, Sakiyama M, Shaukat A, Toyoda Y, Okada Y, Kamatani Y, Nakamura T, Takada T, et al. GWAS of clinically defined gout and subtypes identifies multiple susceptibility loci that include urate transporter genes. Ann Rheum Dis. 2017;76:869–77.

    Article  PubMed  CAS  Google Scholar 

  2. Tin A, Woodward OM, Kao WH, Liu CT, Lu X, Nalls MA, Shriner D, Semmo M, Akylbekova EL, Wyatt SB, et al. Genome-wide association study for serum urate concentrations and gout among African Americans identifies genomic risk loci and a novel URAT1 loss-of-function allele. Hum Mol Genet. 2011;20:4056–68.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  3. Dehghan A, Kottgen A, Yang Q, Hwang SJ, Kao WL, Rivadeneira F, Boerwinkle E, Levy D, Hofman A, Astor BC, et al. Association of three genetic loci with uric acid concentration and risk of gout: a genome-wide association study. Lancet. 2008;372:1953–61.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  4. Li C, Li Z, Liu S, Wang C, Han L, Cui L, Zhou J, Zou H, Liu Z, Chen J, et al. Genome-wide association analysis identifies three new risk loci for gout arthritis in Han Chinese. Nat Commun. 2015;6:7041.

    Article  PubMed  CAS  Google Scholar 

  5. Chang SJ, Tsai MH, Ko YC, Tsai PC, Chen CJ, Lai HM. The cyclic GMP-dependent protein kinase II gene associates with gout disease: identified by genome-wide analysis and case-control study. Ann Rheum Dis. 2009;68:1213–9.

    Article  PubMed  CAS  Google Scholar 

  6. de Pedro N, Cautain B, Melguizo A, Cortes D, Vicente F, Genilloud O, Tormo JR, Pelaez F. Analysis of cytotoxic activity at short incubation times reveals profound differences among annonaceus acetogenins, inhibitors of mitochondrial Complex I. J Bioenerg Biomembr. 2013;45:145–52.

    Article  PubMed  CAS  Google Scholar 

  7. Kottgen A, Albrecht E, Teumer A, Vitart V, Krumsiek J, Hundertmark C, Pistis G, Ruggiero D, O’Seaghdha CM, Haller T, et al. Genome-wide association analyses identify 18 new loci associated with serum urate concentrations. Nat Genet. 2013;45:145–54.

    Article  PubMed  CAS  Google Scholar 

  8. Okada Y, Sim X, Go MJ, Wu JY, Gu D, Takeuchi F, Takahashi A, Maeda S, Tsunoda T, Chen P, et al. Meta-analysis identifies multiple loci associated with kidney function-related traits in east Asian populations. Nat Genet. 2012;44:904–9.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  9. Giri AK, Banerjee P, Chakraborty S, Kauser Y, Undru A, Roy S, Parekatt V, Ghosh S, Tandon N, Bharadwaj D. Genome wide association study of uric acid in Indian population and interaction of identified variants with Type 2 diabetes. Sci Rep. 2016;6:21440.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  10. Yang B, Mo Z, Wu C, Yang H, Yang X, He Y, Gui L, Zhou L, Guo H, Zhang X, et al. A genome-wide association study identifies common variants influencing serum uric acid concentrations in a Chinese population. BMC Med Genomics. 2014;7:10.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Lee J, Lee Y, Park B, Won S, Han JS, Heo NJ. Genome-wide association analysis identifies multiple loci associated with kidney disease-related traits in Korean populations. PLoS One. 2018;13: e0194044.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  12. Merriman TR. An update on the genetic architecture of hyperuricemia and gout. Arthritis Res Ther. 2015;17:98.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  13. Yano H, Tamura Y, Kobayashi K, Tanemoto M, Uchida S. Uric acid transporter ABCG2 is increased in the intestine of the 5/6 nephrectomy rat model of chronic kidney disease. Clin Exp Nephrol. 2014;18:50–5.

    Article  PubMed  CAS  Google Scholar 

  14. Zhang Y, Lee MTM. Serum urate polygenic risk score can improve gout risk prediction: a large-scale cohort study. Front Genet. 2020;11: 604219.

    Article  PubMed  Google Scholar 

  15. Lewis CM, Vassos E. Polygenic risk scores: from research tools to clinical instruments. Genome Med. 2020;12:44.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Choi SW, O'Reilly PF: PRSice-2: Polygenic Risk Score software for biobank-scale data. Gigascience 2019, 8.

  17. Chen CJ, Tseng CC, Yen JH, Chang JG, Chou WC, Chu HW, Chang SJ, Liao WT. ABCG2 contributes to the development of gout and hyperuricemia in a genome-wide association study. Sci Rep. 2018;8:3137.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  18. Zhu Y, Pandya BJ, Choi HK. Prevalence of gout and hyperuricemia in the US general population: the National Health and Nutrition Examination Survey 2007–2008. Arthritis Rheum. 2011;63:3136–41.

    Article  PubMed  Google Scholar 

  19. McAdams MA, Maynard JW, Baer AN, Kottgen A, Clipp S, Coresh J, Gelber AC. Reliability and sensitivity of the self-report of physician-diagnosed gout in the campaign against cancer and heart disease and the atherosclerosis risk in the community cohorts. J Rheumatol. 2011;38:135–41.

    Article  PubMed  Google Scholar 

  20. Wang K, Li M, Hakonarson H. ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res. 2010;38: e164.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  21. Turner S, D.: qqman: an R package for visualizing GWAS results using Q-Q and manhattan plots. bioRxiv 2014,

  22. Pruim RJ, Welch RP, Sanna S, Teslovich TM, Chines PS, Gliedt TP, Boehnke M, Abecasis GR, Willer CJ. LocusZoom: regional visualization of genome-wide association scan results. Bioinformatics. 2010;26:2336–7.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  23. Graham J, McCormick B, Blay S, Shin JH. LDheatmap: an R function for graphical display of pairwise linkage disequilibria between single nucleotide polymorphisms. J Stat Softw. 2006;16:1–10.

    Google Scholar 

  24. Yamamoto YH, Noda T. Autophagosome formation in relation to the endoplasmic reticulum. J Biomed Sci. 2020;27:97.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  25. Yamamoto YH, Kasai A, Omori H, Takino T, Sugihara M, Umemoto T, Hamasaki M, Hatta T, Natsume T, Morimoto RI, et al. ERdj8 governs the size of autophagosomes during the formation process. J Cell Biol. 2020;219(8):e201903127.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  26. Zamudio-Cuevas Y, Fernández-Torres J, Martínez-Nava GA, Martínez-Flores K, Ramírez Olvera A, Medina-Luna D, Hernández Pérez AD, Landa-Solís C, López-Reyes A. Phagocytosis of monosodium urate crystals by human synoviocytes induces inflammation. Exp Biol Med (Maywood). 2019;244:344–51.

    Article  CAS  Google Scholar 

  27. Liao WT, You HL, Li C, Chang JG, Chang SJ, Chen CJ. Cyclic GMP-dependent protein kinase II is necessary for macrophage M1 polarization and phagocytosis via toll-like receptor 2. J Mol Med (Berl). 2015;93:523–33.

    Article  CAS  Google Scholar 

  28. Ilies M, Di Costanzo L, North ML, Scott JA, Christianson DW. 2-aminoimidazole amino acids as inhibitors of the binuclear manganese metalloenzyme human arginase I. J Med Chem. 2010;53:4266–76.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  29. Pudlo M, Demougeot C, Girard-Thernier C. Arginase inhibitors: a rational approach over one century. Med Res Rev. 2017;37:475–513.

    Article  PubMed  CAS  Google Scholar 

  30. Diez-Fernandez C, Rufenacht V, Gemperle C, Fingerhut R, Haberle J. Mutations and common variants in the human arginase 1 (ARG1) gene: impact on patients, diagnostics, and protein structure considerations. Hum Mutat. 2018;39:1029–50.

    Article  PubMed  CAS  Google Scholar 

  31. Sastre M, Galea E, Feinstein D, Reis DJ, Regunathan S. Metabolism of agmatine in macrophages: modulation by lipopolysaccharide and inhibitory cytokines. Biochem J. 1998;330(Pt 3):1405–9.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  32. Li G, Regunathan S, Barrow CJ, Eshraghi J, Cooper R, Reis DJ. Agmatine: an endogenous clonidine-displacing substance in the brain. Science. 1994;263:966–9.

    Article  PubMed  CAS  Google Scholar 

  33. Lortie MJ, Novotny WF, Peterson OW, Vallon V, Malvey K, Mendonca M, Satriano J, Insel P, Thomson SC, Blantz RC. Agmatine, a bioactive metabolite of arginine. Production, degradation, and functional effects in the kidney of the rat. J Clin Invest. 1996;97(2):413–20.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  34. Albina JE, Mills CD, Henry WL Jr, Caldwell MD. Temporal expression of different pathways of 1-arginine metabolism in healing wounds. J Immunol. 1990;144:3877–80.

    PubMed  CAS  Google Scholar 

  35. Cook HT, Jansen A, Lewis S, Largen P, O’Donnell M, Reaveley D, Cattell V. Arginine metabolism in experimental glomerulonephritis: interaction between nitric oxide synthase and arginase. Am J Physiol. 1994;267:F646-53.

    PubMed  CAS  Google Scholar 

  36. Ketteler M, Border WA, Noble NA. Cytokines and L-arginine in renal injury and repair. Am J Physiol. 1994;267:F197-207.

    PubMed  CAS  Google Scholar 

  37. Augustin R, Carayannopoulos MO, Dowd LO, Phay JE, Moley JF, Moley KH. Identification and characterization of human glucose transporter-like protein-9 (GLUT9): alternative splicing alters trafficking. J Biol Chem. 2004;279:16229–36.

    Article  PubMed  CAS  Google Scholar 

  38. Joost HG, Thorens B. The extended GLUT-family of sugar/polyol transport facilitators: nomenclature, sequence characteristics, and potential function of its novel members (review). Mol Membr Biol. 2001;18:247–56.

    Article  PubMed  CAS  Google Scholar 

  39. Stirpe F, Della Corte E, Bonetti E, Abbondanza A, Abbati A, De Stefano F. Fructose-induced hyperuricaemia. Lancet. 1970;2:1310–1.

    Article  PubMed  CAS  Google Scholar 

  40. Hallfrisch J. Metabolic effects of dietary fructose. FASEB J. 1990;4:2652–60.

    Article  PubMed  CAS  Google Scholar 

  41. Döring A, Gieger C, Mehta D, Gohlke H, Prokisch H, Coassin S, Fischer G, Henke K, Klopp N, Kronenberg F, et al. SLC2A9 influences uric acid concentrations with pronounced sex-specific effects. Nat Genet. 2008;40:430–6.

    Article  PubMed  CAS  Google Scholar 

  42. Anzai N, Jutabha P, Amonpatumrat-Takahashi S, Sakurai H. Recent advances in renal urate transport: characterization of candidate transporters indicated by genome-wide association studies. Clin Exp Nephrol. 2012;16:89–95.

    Article  PubMed  CAS  Google Scholar 

  43. Sakurai H. Urate transporters in the genomic era. Curr Opin Nephrol Hypertens. 2013;22:545–50.

    Article  PubMed  CAS  Google Scholar 

  44. Nakayama A, Matsuo H, Takada T, Ichida K, Nakamura T, Ikebuchi Y, Ito K, Hosoya T, Kanai Y, Suzuki H, et al. ABCG2 is a high-capacity urate transporter and its genetic impairment increases serum uric acid levels in humans. Nucleosides Nucleotides Nucleic Acids. 2011;30:1091–7.

    Article  PubMed  CAS  Google Scholar 

  45. Woodward OM, Kottgen A, Coresh J, Boerwinkle E, Guggino WB, Kottgen M. Identification of a urate transporter, ABCG2, with a common functional polymorphism causing gout. Proc Natl Acad Sci U S A. 2009;106:10338–42.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Merino G, van Herwaarden AE, Wagenaar E, Jonker JW, Schinkel AH. Sex-dependent expression and activity of the ATP-binding cassette transporter breast cancer resistance protein (BCRP/ABCG2) in liver. Mol Pharmacol. 2005;67:1765–71.

    Article  PubMed  CAS  Google Scholar 

  47. Hunter DJ. Adventures in the environment and genes. Eur J Epidemiol. 2019;34:1111–7.

    Article  PubMed  Google Scholar 

  48. Mars N, Koskela JT, Ripatti P, Kiiskinen TTJ, Havulinna AS, Lindbohm JV, Ahola-Olli A, Kurki M, Karjalainen J, Palta P, et al. Polygenic and clinical risk scores and their impact on age at onset and prediction of cardiometabolic diseases and common cancers. Nat Med. 2020;26:549–57.

    Article  PubMed  CAS  Google Scholar 

Download references


We thank the Ministry of Science and Technology of Taiwan for supporting this project (MOST 109-2320-B-039 -052, MOST 106-2634-F-390 -001 -CC2), and all the participants to provide the data for the study.


The study was financially supported by the Ministry of Science and Technology of Taiwan (MOST 109-2320-B-039 -052, MOST 106-2634-F-390 -001 -CC2), and in part by China Medical University (CMU105-S-01).

Author information

Authors and Affiliations



YSC and SJC wrote the draft manuscript; YSC, LCY, LTY, and SJC analyzed and interpreted the data; HMC, CCC, and CYC mined the data; JGC, TFJ, and SJC conducted and supervised the study, conceived the study; and JGC and TFJ provided the data source. All authors read and approved the final version.

Corresponding authors

Correspondence to Fuu-Jen Tsai, Jan-Gowth Chang or Shun-Jen Chang.

Ethics declarations

Ethics approval and consent to participate

The ethics committee of the China Medical University Hospital Institutional Review Board in Taiwan (CMUH108-REC1-091) approved this project. Both the Declaration of Helsinki and the Good Clinical Practice Guidelines were followed and informed consent was written by all participants.

Consent for publication

The Taiwan Biobank obtained written informed consent from each participant, and we got Taiwan Biobank approval for this publication.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1:

Supplementary Table1. The distribution of age and uric acid of those participated in the studycohort.

Additional file 2:

Supplementary Table 2. The susceptible variants significantlyassociated with gout compared to normal.

Additional file 3:

Supplementary Table 3. The susceptible variants significantly associated with hyperuricemia comparedto normal.

Additional file 4:

Supplementary Table 4. The susceptible variants significantly associated with gout compared tohyperuricemia. 

Additional file 5:

Supplementary Table5. The susceptible variants of significantly associated with gout andhyperuricemia in those with wild genotype of top significant variants(rs2231142).

Additional file 6:

SupplementaryFigure 1. The linkage disequilibrium of the four variants significantlyassociated with gout in genes DNAJC16 and AGMAT. The variant rs7546668 showedhigh associations with the other three variants (r2 >=0.68). Thered line indicates the cut-off significant p-value by 1e-8.

Additional file 7:

Supplementary Figure 2. The paired linkagedisequilibrium among the variants which were significantly associated withhyperuricemia on chromosome 4 (A) and 12 (B).

Additional file 8:

SupplementaryFigure 3. The Manhattan plots reveal the p-values related to the phenotypesunder limitation of selecting participants for those carrying wild-genotype (GG)of rs2231142. The associations between SNPs and developing gout (A), betweenSNPs and developing asymptomatic hyperuricemia (AH) (B) compared to normal, andbetween gout and AH (C). The red horizontal line denotes cut-off for thep-value of significant difference between the phenotypes and variants by 1e-8.

Additional file 9:

SupplementaryFigure 4. The linkage disequilibrium (LD) of between variant rs671 andrs78069066 and rs77768175 on chromosome 12 for those participants carryingrs2231142 wild-type (GG) in gene ABCG2. All the r-squares of LD between themwere greater than 0.98. The red line indicates the cut-off significant p-valueby 1e-08.

Additional file 10:

SupplementaryFigure 5. The linkage disequilibrium (LD) of between variant rs10805346 andrs3733589, rs1014290 and rsrs3775948 in gene SLC2A9 which were significantlyassociated with hyperuricemia for those participants carrying rs2231142wild-type (GG; ABCG2). All the r-squares of LD between them were greater than0.22. The red line indicates the cut-off significant p-value by 1e-8.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Chang, YS., Lin, CY., Liu, TY. et al. Polygenic risk score trend and new variants on chromosome 1 are associated with male gout in genome-wide association study. Arthritis Res Ther 24, 229 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Gout
  • Hyperuricemia
  • Genome-wide association study
  • Polygenic risk score
  • Chromosome 1