Abstract
Objective. To validate the association between genetic polymorphisms and gout in Japanese patients, and to investigate the cumulative effects of multiple genetic factors on the development of gout.
Methods. Subjects were 153 Japanese male patients with gout and 532 male controls. The genotypes of 11 polymorphisms in the 10 genes that have been indicated to be associated with serum uric acid levels or gout were determined. The cumulative effects of the genetic polymorphisms were investigated using a weighted genotype risk score (wGRS) based on the number of risk alleles and the OR for gout. A model to discriminate between patients with gout and controls was constructed by incorporating the wGRS and clinical factors. C statistics method was applied to evaluate the capability of the model to discriminate gout patients from controls.
Results. Seven polymorphisms were shown to be associated with gout. The mean wGRS was significantly higher in patients with gout (15.2 ± 2.01) compared to controls (13.4 ± 2.10; p < 0.0001). The C statistic for the model using genetic information alone was 0.72, while the C statistic was 0.81 for the full model that incorporated all genetic and clinical factors.
Conclusion. Accumulation of multiple genetic factors is associated with the development of gout. A prediction model for gout that incorporates genetic and clinical factors may be useful for identifying individuals who are at risk of gout.
Epidemiological studies have shown an increasing prevalence of gout. For example, the prevalence of gout in the United States was reported to be 3.9% by the National Health and Nutrition Examination Survey 2007–20081. In Japan, the prevalence of gout has been rising, and has been reported to be 1.1% in Japanese men2. Gout has been reported to be associated with a high risk of death from all causes and coronary heart disease3. Therefore, the identification of individuals with a high risk of developing gout is important.
Multiple genetic and environmental factors are associated with the development of gout. Environmental factors, particularly dietary factors, have been established through epidemiological studies4. Genome-wide association studies (GWAS) in white and African American populations have recently revealed polymorphisms on several loci that are associated with serum uric acid (SUA) levels5,6,7,8,9. Gout is a urate deposition disease, and epidemiological studies have demonstrated that hyperuricemia is the most important risk factor for developing gout10. Several genetic factors associated with SUA have been shown to be risk factors for gout6,7. Studies in Japanese populations have revealed susceptibility loci for SUA and gout11,12,13,14,15. These studies have suggested that the prediction of gout before onset is possible using genetic and environmental factors.
The aims of our study were to validate the associations between multiple genetic loci reported to be associated with SUA level and gout in Japanese; to examine the effects of multiple genetic factors on the development of gout; and to develop a model to discriminate between patients with gout and controls using genetic and clinical factors and to evaluate it using C statistics.
MATERIALS AND METHODS
Study participants
The subjects in this study were previously described12. From April to October 2004, patients with primary gout in the outpatient clinic at the Institute of Rheumatology at the Tokyo Women’s Medical University were asked to participate. The diagnosis of acute gouty arthritis was based on the preliminary classification criteria described by Wallace, et al16. Informed consent was obtained from 188 Japanese patients with gout. Seven samples were excluded because of ethnicity, sex, and availability of DNA; thus 181 male patients were studied. On the first visit, blood samples were collected from all patients and tested for SUA, creatinine, and triglycerides using a Hitachi H7700 autoanalyzer (Hitachi). Estimation of the glomerular filtration rate (eGFR) was based on the equation recommended by the Japanese Society of Nephrology17. In our previous study, genomic DNA from Japanese volunteers was collected in Tokyo18. A total of 1032 Japanese individuals were recruited. Of those, 595 males were selected randomly and served as controls in the present study. Informed consent for the use of their samples was obtained from each subject. Genomic DNA was prepared from peripheral blood using a genomic DNA isolation kit (Qiagen).
This research was reviewed and approved by the Hospital Ethics Committee for Human Genome Research, and all subjects provided informed consent.
Genotype analysis
We selected 10 genes reported to be associated with SUA or gout based on the metaanalysis of GWAS by Kolz, et al9 and a GWAS by Kamatani, et al11; the latter GWAS is the only one reported for SUA in a Japanese population (Table 1). Single-nucleotide polymorphisms (SNP) in SLC2A9, SLC17A1, and SLC22A12 were selected based on previous studies because they have shown significant associations with the development of gout12,14,15. Because the variant rs12356193 in SLC16A99 was nearly monomorphic in the Japanese population, we selected rs2242206, which is estimated to be located in the same haplotype block as rs12356193 based on information obtained from HapMap JPT19. Two polymorphisms, rs2231142 and rs72552713, were selected for ABCG2. These variants have been shown to be genetic risk factors for gout in the Japanese population, and no linkage disequilibrium has been detected between them13. Genotyping was primarily performed using TaqMan SNP Genotyping Assays™ (Applied Biosystems). Detection of each genotype was performed using the ABI Prism 7900HT Sequence Detection System (Applied Biosystems). The data were processed and genotypes were determined using the Sequence Detection System (SDS) version 2.1 software (Applied Biosystems). The rs121907892 genotype in SLC22A12 (c. G774A) was determined according to our previous study12. The rs72552713 genotype in ABCG2 was determined by direct sequencing of PCR products using the forward primer (5′-AAA CCA TAG AAT TTC TTG ACC AG-3′) and the reverse primer (5′-CTG CCT TTT CAC ATA AGT GTC-3′). Sequencing of PCR products was conducted using the Big Dye Terminator cycle sequencing kit and the ABI Prism 3100xl Genetic Analyzer (Applied Biosystems).
Weighted genotype risk score and statistical analysis
The background data of the gout and control groups was compared using Fisher’s exact test for categorical variables and the Mann-Whitney U test for continuous variables. Single genetic marker tests for gout and the controls were evaluated using Fisher’s exact test. OR with 95% CI of each SNP for developing gout were calculated. Statistical analyses were performed using R software (The Institute of Statistical Mathematics; http://cran.ism.ac.jp).
Gene-gene interactions between pairwise combinations of SNP were tested by likelihood ratio tests based on the logistic model with and without interaction. Bonferroni correction for multiple testing was performed.
A weighted genotype risk score (wGRS) was constructed as follows. A weight of 0 was assigned to risk alleles without significant association with gout (p > 0.05). Risk alleles significantly associated with gout (p < 0.05) were assigned a weight of 1, 2, or 3 based on OR. Risk alleles with a 1 < OR ≤ 2 and 2 < OR ≤ 3 were assigned as weight 1 and 2, respectively. The OR of SNP rs121907892 (SLC22A12) was infinity. Because of the strong association between this SNP and gout, a weight of 3 was assigned to this allele. Table 1 lists the weight of each risk allele. The wGRS for an individual was the sum of the number of risk alleles for each SNP multiplied by its weight. The cumulative effect of the wGRS on gout was examined by a logistic regression model with adjustment for age, body mass index (BMI), triglyceride (TG) levels, and eGFR. The proportion of individuals with wGRS ≤ 12 was 26%; these patients served as a reference group20. The OR for individuals with each wGRS was calculated relative to the reference group.
Models to discriminate patients with gout from controls were constructed using genetic and/or clinical factors. The wGRS was incorporated as a genetic factor. Clinical factors were constructed from age, BMI, TG, and eGFR. The latter 3 variables were included as clinical factors because their distributions were significantly different between patients with gout and controls (Table 2). We developed 5 models: model 1 was constructed with a genetic factor (wGRS) alone; model 2 was constructed with age and BMI; model 3 included age, BMI, and TG; model 4 included all clinical factors; and model 5 was constructed with all clinical and genetic factors. The discriminative capability of the models was evaluated using logistic regression to generate C statistics21. A C statistic < 0.6 has a low clinical value, while a C statistic > 0.8 has clinical utility22.
RESULTS
Genotyping was completed in 153 patients in the gout group and 532 controls. The demographic data of the subjects are shown in Table 2.
Association of individual SNP with gout
The distribution of allelic frequencies is shown in Table 1. The allele frequency for each polymorphism was in Hardy-Weinberg equilibrium in the controls. The polymorphisms in GCKR, ABCG2, LRRC16A, LRP2, SLC17A1, SLC2A9, and SLC22A12 were significantly associated with gout (p < 0.05).
Gene-gene interactions
The allele frequencies of rs72552713 in ABCG2 and rs121907892 in SLC22A12 were too low to examine the gene-gene interactions. The gene-gene interactions between some combinations were not evaluated because the logistic regression was not convergent. Although the combination of SLC16A9 (rs2242206) and ABCG2 (rs2231142) showed p values of 0.018, the Bonferroni correction removed significance (p < 0.0014). Therefore, no combinations with statistical significance were detected.
Distribution of weighted genotype score and the risk of gout
The distribution of patients with gout and controls according to the wGRS is shown in Figure 1. The means (SD) of wGRS of the patients with gout and the controls were 15.2 (2.01) and 13.4 (2.10), respectively (p < 0.0001). Figure 2 shows OR for gout in individuals with each wGRS relative to individuals with wGRS ≤ 12 by logistic regression analysis. OR increased with the increasing wGRS. The OR of developing gout in individuals with wGRS ≥ 18 relative to individuals with wGRS ≤ 12 was 19.7 (95% CI 7.1–55.0, p = 1.2 × 10−8). The logistic regression analysis showed that the increase of one wGRS resulted in an increase of 1.5 OR (95% CI 1.4–1.70, p = 2.4 × 10−14).
The discriminative capability of a prediction model that includes genetic and clinical factors
We constructed a model to discriminate between patients with gout and the controls using wGRS or clinical factors (Table 3). In model 1, which incorporated wGRS alone, the C statistic was 0.72. Models 2, 3, and 4 were constructed only with clinical factors, and an improvement in the C statistic ranging from 0.69 to 0.73 was observed as the number of clinical factors increased. The highest C statistic was obtained when both genetic and clinical factors were incorporated into the model (model 5, with a C statistic = 0.81).
DISCUSSION
We validated the association between 6 out of 11 SNP and gout in a Japanese population and demonstrated the cumulative effects of genetic variants on the development of gout using wGRS. A model that includes genetic and clinical factors was developed to predict gout. The C statistic of the model constructed by all clinical factors (model 4) was 0.73, which was comparable to that of the model that included the genetic factor alone (model 1). However, when the genetic factor was combined with all clinical factors (model 5), the C statistic was increased to 0.81, indicating that both genetic and clinical factors are important in the development of gout. The result of the case-control discrimination may suggest the usefulness of this approach to predict gout.
Previous studies, including GWAS, have revealed genetic variants that are associated with SUA level. Based on these data, additive genetic risk scores have been shown to be associated with SUA23,24,25,26. However, associations between genetic risk scores and gout have not been studied extensively. Dehghan, et al reported that the OR for gout was increased 6 to 41 times across the GRS in white or African American cohorts23. Our study validated these results in Japanese males and demonstrated the polygenic character of gout. Further, the study suggested the capability of our model for predicting gout using genetic and clinical factors. The predictive power of this model should be evaluated in cohort studies.
The contribution of genetic variants in case discrimination or prediction has been investigated in several diseases other than gout. In type 2 diabetes, a disease commonly associated with gout, the C statistics of the prediction models using genetic information alone have been reported to range from 0.54 to 0.6320,27,28,29. The contribution of genetic variants in the prediction of type 2 diabetes is modest. Minimal effects of genetic variants for prediction have also been reported for obesity and premature coronary disease30,31. However, De Jager, et al indicated that the C statistic of the genetic model for multiple sclerosis was 0.732, comparable to the C statistic observed in our present study. Although the predictive or discriminatory capabilities of genetic factors are expected to vary between diseases, differences in methodologies may influence the divergent results. The genetic aspects of uric acid have been indicated to be a good example of the polygenic trait theory33, and all of the genetic variants in this study have been shown to be associated with SUA level. The predictive and discriminatory powers of genetic variants in gout may be greater than in type 2 diabetes because gout develops as a result of hyperuricemia.
There are several limitations of our study. First, this was a case-control study with a small sample size and the predictive model was developed based on the study population data. We did not validate the models in separate subjects. Therefore, the discriminative power of our model needs to be interpreted with caution. The predictive ability of the model should be validated using a cohort with independent subjects. Second, a limited number of clinical factors were evaluated. A model including more susceptibility loci and clinical factors may be necessary for better prediction. Third, we could not evaluate gene-gene interactions in some of the pairwise combinations of the common SNP. Gene-gene interaction should be evaluated in a larger sample size in a future study. Finally, β-coefficients were not used to construct the GRS in this study. In some studies, GRS calculation is based on a simple count27, while in others, the GRS is weighted based on the β-coefficient. Although β-coefficients will represent the best estimates of risk available34, no obvious differences have been reported between the simple count and wGRS methods29,34. The β-coefficient for rs121907892 could not be calculated in our study. However, it may be prudent to construct GRS using discrete weights. When the GRS was constructed by counting the number of risk alleles with p < 0.05, the mean (SD) GRS in patients with gout and controls were 10.1 (1.57) and 8.9 (1.70), respectively (p < 0.0001), indicating the association between gout and accumulation of genetic risks constructed by a nonweighted GRS. In our study, we used wGRS because some of the variants are rare alleles with considerably higher OR for the development of gout than others (rs2231142 and rs72552713 in ABCG2 and rs121907892 in SLC22A12; Table 1).
Our study revealed the multifactorial character of gout. Identifying individuals at high risk of developing gout in the general population is suggested to be feasible using a prediction model consisting of genetic and clinical factors. Because gout is a risk factor for cardiovascular disease and metabolic syndrome3,35,36, prediction of gout can be beneficial in preventing these diseases. Additional studies that include more genetic and clinical factors in a larger number of subjects independent of the present study population will be necessary to refine the model proposed here.
Footnotes
-
Supported by a grant from the Ministry of Education, Culture, Sports, Science and Technology, Japan, and a grant from the Gout Research Foundation.
- Accepted for publication April 15, 2013.