- Research article
- Open Access
The comparative effect of exposure to various risk factors on the risk of hyperuricaemia: diet has a weak causal effect
Arthritis Research & Therapy volume 23, Article number: 75 (2021)
Prevention of hyperuricaemia (HU) is critical to the prevention of gout. Understanding causal relationships and relative contributions of various risk factors to hyperuricemia is therefore important in the prevention of gout. Here, we use attributable fraction to compare the relative contribution of genetic, dietary, urate-lowering therapy (ULT) and other exposures to HU. We use Mendelian randomisation to test for the causality of diet in urate levels.
Four European-ancestry sample sets, three from the general population (n = 419,060) and one of people with gout (n = 6781) were derived from the Database of Genotypes and Phenotypes (ARIC, FHS, CARDIA, CHS) and UK Biobank. Dichotomised exposures to diet, genetic risk variants, BMI, alcohol, diuretic treatment, sex and age were used to calculate adjusted population and average attributable fractions (PAF/AAF) for HU (≥0.42 mmol/L [≥7 mg/dL]). Exposure to ULT was also assessed in the gout cohort. Two sample Mendelian randomisation was done in the UK Biobank using dietary pattern-associated genetic variants as exposure and serum urate levels as outcome.
Adherence to dietary recommendations, BMI (< 25 kg/m2), and absence of the SLC2A9 rs12498742 urate-raising allele produced PAFs for HU of 20 to 24%, 59 to 69%, and 57 to 64%, respectively, in the three non-gout cohorts. In the gout cohort, diet, BMI, SLC2A9 rs12498742 and ULT PAFs for HU were 12%, 49%, 48%, and 63%, respectively. Mendelian randomisation demonstrated weak causal effects of four dietary habits on serum urate levels (e.g. preferentially drinking skim milk increased urate, β = 0.047 mmol/L, P = 3.78 × 10−8). These effects were mediated by BMI, and they were not significant (P ≥ 0.06) in multivariable models assessing the BMI-independent effect of diet on urate.
Diet has a relatively minor role in determining serum urate levels and HU. In gout, the use of ULT was the largest attributable fraction tested for HU.
Hundreds of genetic variants are associated with serum urate levels [1–3], and observational studies have associated individual dietary factors (e.g. alcohol, sugar-sweetened beverages, coffee, red-meat consumption [4–8]) and overall eating habits [9,10] with urate levels, along with other environmental (e.g. diuretic use [11, 12]) and endogenous factors (e.g. age and sex ). Understanding the importance of risk factors and their causal relationship (if any) with HU is critical in developing strategies for the prevention of HU and gout. However, addressing causality is challenging. Observational, longitudinal or migratory studies, and temporal correlations can only be regarded as hypothesis-generating owing to the intractable issue of unmeasured confounding. Any attempts to draw conclusions with respect to causality, including using causal language/inferences, even from an accumulation of studies, mis-represents the evidence . In this context, we note that gout is a multi-stage process beginning with HU, progressing to deposition of monosodium urate crystals and culminating in an innate immune response to crystals . Not all people with HU develop gout , so HU and the progression from HU to gout should not be conflated when considering possible causal risk factors.
The gold standard for testing an exposure for a causal role is the randomised clinical trial (RCT). This approach has demonstrated causality for dissolved sugar (sugar-sweetened beverages) in raising urate levels [17–19]. However, for the majority of suspected causal exposures an RCT is not possible. Mendelian randomisation (MR) exploits the natural randomisation of alleles causal for a particular exposure and is analogous to a RCT. Several MR studies have shown a small causal effect of BMI on urate levels (0.0045 to 0.010 mmol/L [0.075–0.17 mg/dL] increase in serum urate per unit increase in genetically determined BMI [20–22]). Dietary preferences have a heritable component and genetic associations have been reported [23–25]. Mendelian randomisation in the UK Biobank, using genetic variants associated with dietary patterns has demonstrated that a ‘healthful’ versus ‘unhealthful’ dietary pattern is not strongly causal for coronary heart disease or type 2 diabetes, despite diet being strongly correlated with these diseases .
Using the widely applied approach of decomposing variance, where the sum of multiple risk factors included in a model is constrained to 100%, overall diet contributed ≤0.3% of variance in urate levels, substantially less than the 23.9% explained by inherited common genetic variants . Why so little variance is explained is unclear, but one possible reason is that overall diet, which comprises some foods associated with increased urate and some foods associated with decreased urate, does not play a strong causal role. This possibility is supported by a RCT that reported a small 0.021 mmol/L (0.35 mg/dL) reduction in urate levels in people following the Dietary Approaches to Stop Hypertension (DASH) diet compared to those on a ‘typical’ US diet . Another RCT comparing the Mediterranean diet to a ‘prudent Westernised diet’ reported a small reduction in serum urate levels (0.010 mmol/L [0.17 mg/dL]) over 5 years . Another possible explanation postulated in ref.  for the small amount of variance explained is low variability in diet within the US cohorts used in . However, differences in diet between men and women, across age groups, socioeconomic status, ethnicity and BMI strata in the US have been reported —to observe these differences there must be variability in the diet.
Population attributable fraction (PAF) is the proportion of cases for an outcome within a population that can be attributed to a given risk factor, incorporating both the prevalence and the effect size of the exposure . The sum of PAFs for multiple risk factors for a single condition is not constrained to 100%, because an outcome can have multiple risk pathways population-wide. Using the Third National Health and Nutrition Examination Survey PAFs were reported of 44% for being overweight or obese (implying that 44% of HU would be prevented if the entire population had BMI < 25 kg/m2) and 9% for non-adherence to a DASH-style diet . Variances explained were 8.3% and 0.1%, respectively , indicating the two methods agree as to which exposure has a greater impact.
Our first aim was to use attributable fraction to compare the contributions to HU of various genetic, environmental and endogenous risk factors, including lack of use of urate-lowering therapy. The second aim was to use MR to test for a causal role of diet in determining urate levels.
Participants and methods
Participants and data collection
The attributable fraction analysis included four distinct cohorts of European ancestry (Table 1)—cohorts 1 to 3 are population-based and the fourth cohort is comprised entirely of people with gout.
Cohort 1 comprised 14,247 participants of European ancestry from the US population—7342 from the Atherosclerosis Risk in Communities (ARIC) Study, 1314 from the Coronary Artery Risk Development in Young Adults (CARDIA) Study, 2513 from the Cardiovascular Health Study (CHS) and 3078 from the Framingham Heart Study (FHS). These numbers exclude people without serum urate measurements or genome-wide genotypes, along with individuals aged under 18 years, people with kidney disease or gout and those taking urate-lowering therapy. People who answered less than 10% of the food frequency survey, those whose estimated average daily calorie intake was less than 600 kcal/day or greater than 4200 kcal/day and those whose questionnaire answers were deemed unreliable by the study interviewer at recruitment were also excluded. For ARIC, one person from each first degree-related pair was excluded.
Cohort 2, cohort 3 and the gout cohort were sourced from the UK Biobank resource. Subjects of European ancestry and who had urate measures and genotypes available were included in the analysis. Relatives with kinship coefficients > 0.177 were removed, and one person from each relationship was kept, with a preference for keeping gout-affected participants. Those who self-reported having kidney disease were also removed. The gout cohort comprised people who self-reported having gout at visit 0 or were being treated with urate-lowering therapy (n = 6781) , this case definition has been validated [31, 32]. Cohort 2 consisted of UK Biobank participants who answered a 24-h dietary recall questionnaire during the assessment visit (n = 57,251) and cohort 3 consisted of the remaining UK Biobank subjects who answered a reduced food frequency questionnaire (n = 347,526). We excluded, from cohort 2, subjects who had energy intakes in excess of 18,000 kJ for females and 20,000 kJ for males based on their 24-h dietary recall data, those who had unreliable dietary data as flagged by the recruiter, or subjects not eating normally due to illness or fasting. No additional exclusions were applied for cohort 3 or the gout cohort. Participants used for the MR analysis were also from the UK Biobank, and the MR cohort has been described previously .
Collection of dietary and serum urate data are described in the Supplemental material. Collection of genetic data for the ARIC, FHS, CHS and CARDIA cohorts is described in ref.  and for the UK Biobank in .
To calculate a PAF, all exposure and outcome variables must be dichotomous. The outcome for this study was HU, defined as serum urate ≥0.42 mmol/L [≥7 mg/dL] for men and women . For the genetic exposures, HU risk alleles were defined as urate-increasing under a dominant model . Determination of dichotomised dietary exposures is described in the Supplemental material. Alcohol exposure was defined as > 1 drink per week, being overweight/obese as BMI ≥25 kg/m2, age was dichotomised as ≥50 versus < 50 years partly in order to capture menopause as a risk factor in women and diuretic use either self-reporting or not self-reporting diuretic intake—these variables were the same as those for which PAR estimates were calculated in the Third National Health and Nutrition Examination Survey in ref.  and for age. In the gout cohort, self-reported treatment with urate-lowering therapies, allopurinol (n = 4841), probenecid (n = 3) and sulphinpyrazone (n = 21) (the only three urate-lowering medications for which baseline medication data were available) was a dichotomised exposure variable.
All analyses were performed using R v3.6.1 in RStudio 1.2.5019. For the various exposures, the PAF calculation was (frequency of exposure in cases) × (ORExposure – 1)/ORExposure) . Odds ratios for the risk of HU for these exposures were calculated in a logistic regression multivariable model including all other environmental and endogenous exposure variables and SLC2A9 rs12498742 genotype—this variant was chosen for individual focus because of its large effect on serum urate levels . For the percent variance explained analysis (Table 3), effect sizes (β) on serum urate levels for the same exposures were calculated in a linear regression multivariable model including all other environmental and endogenous exposure variables and SLC2A9 rs12498742 genotype. Average attributable fractions (AAFs), adjusted for all other exposure variables and rs12498742 genotype (Table 2) were calculated using a multivariable model in the R function averageAF , described in more detail in Supplemental material. Age was removed as it did not confer risk in the gout cohort. In general, AAFs are lower than PAFs calculated from the same data (as consistently observed here) and it has been proposed that they reflect the most plausible/reliable result across the many different methods of calculating attributable fractions .
For genetic variants in Table S2, PAFs (and odds ratios) were calculated per allele under an additive model using the method of Rockhill et al.  ((frequency risk allele in cases) × ((ORRiskAllele – 1)/ORRiskAllele)). However, PAFs (and odds ratios) in Tables 2, S1 and Fig. 1 were calculated under a dichotomised (dominant) model with the exposed group defined as those with one or more risk allele and were adjusted for dichotomised age, sex, BMI, diet, alcohol and diuretic exposure risks. All AAF analyses used a dominant model for SNP analysis because risk was dichotomised.
Two-sample MR using the MendelianRandomization R package  was tested for a causal role of dietary habits in determining urate levels. GWAS summary statistics for urate levels in Europeans were obtained from ref.  comprising 288,649 individuals; 101 loci with non-ambiguous lead SNPs where no strand assignation issues could have arisen by GWAS which can occur with A/T and G/C variants where the alternative allele is the same as the potential strand-flipped allele. This conservative approach was taken in order to harmonise effect alleles from a meta-analysis of multiple studies using multiple imputation panels. In order to test that this QC did not influence the results, we conducted a sensitivity analysis and repeated the primary Mendelian randomisation with the inclusion of all ambiguous SNPs (Figure S1)—the results were in high concordance. GWAS data for dietary habits were obtained from a study of 455,146 individuals of European ancestry from the UK Biobank . The latter were derived from consumption patterns for 47 single foods asked about in the reduced food frequency questionnaire or 40 principal component-derived dietary patterns with ≥3 non-ambiguous genome-wide significant index SNPs reported in ref.  that were also present in the serum urate GWAS summary statistics . Independence of loci was based on the distance-based pruning approach within the source GWAS [3, 26]—loci had to be > 500 kb apart and if there were two in a 1 Mb interval independence was established by plotting and visualisation of inter-marker linkage disequilibrium. The MR was performed using three methods available in the MendelianRandomization R package—inverse-variance-weighted meta-analysis , MR Egger (enables detection of pleiotropy ) and weighted median (robust to pleiotropy ). Significance thresholds were set at P < 0.05/87 (5.8 × 10− 4) for the inverse-variance-weighted analysis, P < 0.05 for the weighted median analysis and intercept P > 0.05 for the MR Egger analysis. Our strategy was to identify causal effects from the inverse-variance-weighted analysis corrected for multiple testing followed by sensitivity analysis by weighted median and MR Egger to test the robustness of the results. All three significance thresholds had to be met for a dietary habit to be considered causal. Multivariable MR was conducted to investigate the possible upstream impact of BMI as a common causal link between dietary patterns and serum urate levels using the likelihood-based method on summary data as described by Burgess and Thompson , as implemented in the MendelianRandomization R package. This technique uses pleiotropic genetic variants to estimate the direct effect of multiple exposures on an outcome (e.g. BMI and diet on serum urate), with causal estimates representing the independent causal effect of each exposure on the outcome, not operating through the other exposure included in the analysis.
Population and average attributable fractions
Male sex had the largest PAF and AAF (64.3 to 75.6%, and 29.8 to 30.9%, respectively, in cohorts 1 to 3), being overweight or obese had the second largest measures (PAF = 59.2 to 68.6%, AAF = 23.5 to 26.2%), and inheriting the SLC2A9 rs12498742 A-allele had the third largest effect (PAF = 56.6 to 63.8%, AAF = 22.0 to 23.5%). All other factors had PAFs < 25% and AAFs < 10% (Table 2). Cohorts 2 and 3 had the same ranking of risk factors (sex > BMI > rs12498742 > diet adherence > alcohol consumption > diuretic therapy > age), while cohort 1 ranked the same three risk factors first (sex > BMI > rs12498742), before differing in rankings of the remaining four risk factors (diuretic therapy > diet adherence > age > alcohol consumption) (Fig. 1).
In the gout cohort, lack of treatment with urate-lowering therapy had the second largest PAF of 63.2%, after sex (76.7%) (Table 2). The PAFs for BMI and SLC2A9 rs12498742 were lower than for the non-gout cohorts (57 to 69% in non-gout and 48 to 49% in gout), and non-adherence to the Healthy Eating Pyramid guidelines was 21.4% in non-gout (cohort 3) and 11.8% in gout. Average attributable fraction values showed a similar trend (Fig. 1).
In sex-stratified analysis in the non-gout cohorts, PAFs for rs12498742 were 77.7 to 96.4% in women compared to 49.2 to 58.9% in men (Table S1). Alcohol was not a risk factor in women across all cohorts, nor age in men in cohorts 2 and 3, and gout (Table S1; 95% CI encompassed 1.0). The lower age limit for recruitment into the UK Biobank, which these two cohorts were derived from, was 40 years, which may have influenced the calculation.
In the non-gout cohorts, of 30 genetically-independent serum urate-associated genetic variants chosen as having the top effects by GWAS  evaluated (Table S2), the SLC2A9 rs12498742 variant was the largest, with PAFs ranging from 28.5 to 32.1% and AAFs from 22.0 to 23.5%. For comparative purposes, we summed the PAFs for genetic variants, assuming that the variants act independently of each other to influence the risk of HU, with the individual PAFs summing to > 141% (cohort 1 was 146.2%, cohort 2 was 143.7%, cohort 3 was 141.3%). Summing the AAFs (equivalent to the summing of PAFs, above) resulted in all three cohorts having a summed AAF over 87% (cohort 1 was 97.9%, cohort 2 was 87.1%, cohort 3 was 101.6%). Summed attributable fractions for genetic variants for HU were considerably lower for the gout cohort (PAF was 77.6% and AAF was 44.0%), possibly reflecting selection (collider) bias.
Percent variance explained for serum urate levels
In the non-gout cohorts, sex had the most percent variance explained (22 to 27%) (Table 3). The dichotomised BMI exposure was consistently 7 to 9%, with diuretic exposure accounting for 12% variance in the US-based cohort and 4–5% in the UK-derived cohorts. The diet estimate was ≤0.1% and SLC2A9 was 2–3%, similar to our previous report . The gout cohort included urate-lowering therapy exposure in the model, with exposure accounting for the largest proportion (35%) of variance, approximately 10-fold more than any other variable. The use of percent variance explained produced a broadly similar ranking order of risk factors to the PAF and AAF analyses across all four cohorts (Fig. 1).
Five of the 87 single foods and principal component-derived dietary-associated habits  provided evidence of a causal effect (IVW P < 0.05/87 (5.7 × 10− 4)) on urate levels by inverse-variance-weighted MR (Table S3). All five of these dietary habits also had no evidence for an intercept significantly different from zero in the MR Egger analysis (all P > 0.05) indicating no evidence for directional (horizontal) pleiotropy. Four of these dietary habits provided evidence for a causal role (P < 0.05) and yielded similar effect sizes in the weighted median analysis (Table 4). Two of these causal effects were with dairy-related dietary habits (preferentially drinking skim milk and preferentially drinking milk with a higher fat content), and the other two causal effects were for consuming tub margarine and daily dried fruit consumption.
Of the 39 genetic variants that comprised the four dietary-associated habits, 21 are associated with metabolic traits (http://www.type2diabetesgenetics.org/ [accessed: 2nd June 2020]) and/or traits available in the UKBiobank PheWeb (http://pheweb.sph.umich.edu:5000 [accessed: June 2, 2020]), including 16 specifically associated with BMI or a related body fat trait (Table S4). To test the possibility that the causal association between these four dietary habits and urate levels is due to BMI as a common upstream cause (e.g. change in dietary habits due to weight-loss advice), we applied multivariable MR using the same individual level UK Biobank dataset described in ref. . For all four dietary patterns, including BMI in the multivariable analysis resulted in no evidence for a causal effect (P ≥ 0.06), BMI showing a causal relationship with urate levels independent of the dietary habit (Fig. 2). Bidirectional MR between BMI and each of the four dietary habits where, by inverse variance-weighted meta-analysis MR BMI was tested for a causal effect on the dietary habits and each of the dietary habits was tested for a causal effect on BMI, conducted to confirm whether BMI is a common upstream cause of dietary habits, provided evidence in both directions (P ≤ 4.4 × 10− 18 for BMI to dietary habit, P ≤ 9.4 × 10− 4 for dietary habit to BMI)—a situation termed “correlated pleiotropy” —except in the margarine analysis for the BMI to dietary habit analysis (P = 0.24) although there was evidence for the dietary habit to BMI analysis (P = 8.4 × 10− 61). This indicates that BMI and the four dietary habits are strongly correlated traits or work through a shared pathway and that the four dietary habits have no effect on urate levels independent of BMI.
Our previous study  concluded, using percent variance explained, that common genetic variants have a greater contribution to urate levels in the non-gout population than overall diet. Using attributable fraction measures, we arrive at the same conclusion, importantly also in a cohort of people with gout. Previously, the summed percent variance for the 30 genetic variants for urate levels was 8.7%, considerably greater than the variance explained by the DASH diet . Here, the summed PAFs for the 30 genetic variants was 141 to 146%, considerably greater than that for following the DASH diet recommendations in cohorts 1 and 2. Thus, empirically for HU at least, the different approaches of decomposition of variance and use of attributable fractions provide similar support for the greater relative role of common inherited genetic variation than overall diet in determining urate levels and HU. In the gout cohort, the attributable fractions for urate-lowering therapy were greater than for diet and BMI < 25 kg/m2. Acknowledging the limitation that we were unable to build compliance, medication dose, and dosing to target into our models (which would contribute to under-estimating the effect of urate-lowering therapy), our data emphasise the importance of gold-standard clinical practice (urate-lowering therapy), to manage HU in gout patients. While weight reduction has established benefits, including to co-morbidities in gout, our data demonstrate the greater impact of urate-lowering therapy in managing HU in gout.
While it is debatable whether public health efforts should be directed to primary prevention of HU, given the lack of evidence that HU is directly causal of conditions other than gout , there are two considerations that can be drawn. One, efforts would need to focus on interventions for which there is unequivocal evidence for a substantial impact to be made. This is not the case for the DASH diet (our AAF estimate in a multivariable model of the proportion of cases of HU prevented by following a DASH diet was only 6 to 7%). Two, the proportion of cases of HU attributable to being overweight or obese from the population was 24 to 26% in the same model, only slightly more than the proportion attributable to SLC2A9 rs12498742 (22 to 24%). It may seem incongruous to compare these exposures in the context of possible public health approaches to prevent primary HU, given that it is not possible to prevent exposure to a common genetic variant. It is, however, possible to modify the impact of a genetic variant. The uricosuric drugs benzbromarone and probenecid inhibit the reuptake of filtered urate by GLUT9 (encoded by SLC2A9) ; thus, it is conceptually possible to target individuals with the rs12498742 urate-raising allele to improve excretion of urate and prevent HU. From a public health perspective, this is likely a more tractable intervention (in that it targets a single measurable exposure) than preventing obesity, which is caused by multiple environmental and genetic exposures that are not well understood.
That individual foods and estimates of dietary habits associate strongly with urate levels in observational data  does not necessarily translate into a clinically significant causal effect. It is interesting to compare association data of the DASH diet score  with data from a RCT of the effect of the DASH diet on serum urate levels —the association data show a decrease of 0.023 mmol/L [0.38 mg/dL] between the least and most DASH-like diets in the US population , very similar to the 0.021 mmol/L [0.35 mg/dL] decrease when comparing the DASH diet with an ‘average American diet’ in the RCT . In both cases, this is a relatively small change attributable to dietary habits and is reflective of the evidence presented here for a weak BMI-mediated causal relationship between diet and urate levels.
An incongruity is the apparent inconsistency between the two dairy-related MR analyses and results from RCTs [46–48]. Using MR as a complementary approach to investigate causality, we found only a small number of weak causal associations between dietary habits and serum urate. Interestingly, two of the significant causal associations represent opposing dietary habits, namely preferentially drinking skim milk or preferentially drinking milk with a higher fat content. The causal effects observed were consistent with these being opposing dietary habits, with skim milk consumption associating with increased urate, whilst consumption of higher-fat milk associated with decreased urate at an approximately equivalent effect size (0.050 mmol/L [0.84 mg/dL] vs. −0.044 mmol/L [− 0.73 mg/dL], respectively). However, whilst these results are consistent with each other they are not consistent with prior studies of milk and dairy proteins in relation to urate levels. Observational studies have reported an inverse relationship between consumption of dairy products and serum urate levels [7, 10, 49–52]. Many of these observational studies do not separate dairy products into low and high fat content; however, those that do have found that this effect appears to be limited to consumption of skim or low-fat dairy products [10, 51]. RCTs have supported these observational findings [46–48], in particular consumption of skim milk products acutely lowered serum urate levels by approximately 10% in 16 healthy adult men . The apparent inconsistency between the MR and RCT results can be explained by the influence of BMI on the MR analysis. BMI appears to be a common upstream cause in the two dairy-related MR associations reported here, and these two dairy-related dietary habits are highly correlated with BMI, several measures of body fat and weight-loss related traits, including making major dietary changes to lose weight  (Table S4). It is plausible that the MR results reflect dietary recommendations given to individuals with a higher BMI (drink skim or low-fat milk), explaining the contradictory results seen here.
Our BMI genetic instrument explains more variance in type of milk consumed (~ 0.5%), than the milk type instruments do themselves (0.04 to 0.1%) , highlighting an important limitation to these analyses. Genetic instruments for dietary habits likely explain a small fraction of phenotypic variance  or may be linked to diet through indirect mechanisms, potentially subjecting the MR analysis to bias towards the null, pleiotropy or confounding. While multiple MR approaches were used to address some of these pitfalls, future investigation using more biologically based genetic instruments for diet may illuminate previously undetectable causal relationships.
In the sample sets of European ancestry studied here, SLC2A9 rs12498742 had a considerably greater PAF than the ABCG2 rs2231142 variant (29 to 32% vs. 6%, respectively (Table S2)). This is because of the 1.7-fold increased effect size of rs12498742 on serum urate levels and the increased prevalence of the urate-increasing allele (77% vs. 11%) . In contrast, in a Japanese study, the PAF for rs2231142 was 29%, compared to 19% for being overweight or obese , suggesting that for any primary prevention of HU in the Japanese population, targeting ABCG2 dysfunction would be a strategy to be considered. The rs2231142 risk allele frequency is 29% in the East Asian population compared to 9% in the European population. The authors of the Japanese study concluded that ABCG2, at least, is a stronger risk factor for HU than other ‘typical’ environmental risk factors .
One limitation of the gout cohort analysis is the possibility of selection (collider) bias resulting from conditioning the sample set on gout ascertainment which would serve, when testing variables that are risk factors for gout per se, to bias effect sizes towards the null or even in an opposing direction . This phenomenon likely explains the reduced (reversed for age) effect sizes for age, sex, BMI and diuretic exposure for each of risk of HU and change in serum urate levels and reduced variance explained, evidenced by non-overlapping 95% CIs compared to cohorts 1 to 3 (Tables 2 and 3). For SLC2A9, effect sizes and variance explained were lower in the gout cohort, but some confidence intervals were overlapping. However, for diet and alcohol, there was no difference in effect sizes (the 95% CIs overlapped) suggesting that collider bias does not have a substantive effect on these estimates within the gout cohort. We note that the prevalence of healthy eating diet non-adherence was very similar between the UK Biobank Gout cohort and the equivalent non-gout cohort (cohort 3) indicating that diagnosis of gout did not change dietary behaviour. Selection bias will not influence our effect estimates for urate-lowering therapy; however, estimates for this exposure are likely inflated in the UK Biobank owing to healthy volunteer selection bias . This likely leads to an over-estimate of effect size owing to exposure to urate-lowering therapy, because of a more compliant demographic. Our estimate of OR = 20.2 (Table 2) is considerably higher than a hazard ratio of 4.5 reported for achieving target urate in a gout cohort drawn from the UK primary care population . While our estimate is not representative of the general population, it does indicate the possibility that the relative effect of urate-lowering therapy on HU and serum urate levels is higher when compliance to urate-lowering therapy is increased.
In conclusion, we demonstrate using attributable fraction measures, that incorporate both the prevalence of exposure and effect size, the considerably greater attributable fraction of HU in the general population owing to common inherited genetic variants and BMI than to dietary exposure. The use of urate-lowering therapy in gout was the largest contributor to attributable fraction of HU. These findings are consistent with previous findings from the use of variance explained in the general population . There is a weak causal effect between four dietary habits and urate levels, all mediated by BMI. Collectively, our findings refute the widely held perception that HU is primarily caused by diet [57–60].
Availability of data and materials
All data used were publicly available. Derivative datasets generated during the current study are available from the corresponding author on reasonable request.
Average attributable fraction
ATP-binding cassette subfamily G member 2
Atherosclerosis risk in communities
Body mass index
Coronary artery risk development in young adults
Cardiovascular heart study
Dietary approaches to stop hypertension
Framingham heart study
Glucose transporter 9
Genome-wide association study
Inverse variance weighted
Population attributable fraction
Randomised clinical trial
Solute carrier family 2, member 9
Boocock J, Leask M, Okada Y, Matsuo H, Kawamura Y, Shi Y, et al. Genomic dissection of 43 serum urate-associated loci provides multiple insights into molecular mechanisms of urate control. Hum Mol Genet. 2020;29:923–43.
Nakatochi M, Kanai M, Nakayama A, Hishida A, Kawamura Y, Ichihara S, et al. Genome-wide meta-analysis identifies multiple novel loci associated with serum uric acid levels in Japanese individuals. Comm Biol. 2019;2:115.
Tin A, Marten J, Kuhns V, Li Y, Wuttke M, Kirsten H, et al. Target genes, variants, tissues and transcriptional pathways influencing human serum urate levels. Nat Genet. 2019;51:1459–74.
Batt C, Phipps-Green AJ, Black MA, Cadzow M, Merriman ME, Topless R, et al. Sugar-sweetened beverage consumption: a risk factor for prevalent gout with SLC2A9 genotype-specific effects on serum urate and risk of gout. Ann Rheum Dis. 2014;73:2101–6.
Choi HK, Curhan G. Beer, liquor, and wine consumption and serum uric acid level: the Third National Health and Nutrition Examination Survey. Arthritis Rheum. 2004;51:1023–9.
Choi HK, Curhan G. Coffee, tea, and caffeine consumption and serum uric acid level: the Third National Health and Nutrition Examination Survey. Arthritis Rheum. 2007;57:816–21.
Choi HK, Liu S, Curhan G. Intake of purine-rich foods, protein, and dairy products and relationship to serum levels of uric acid: the Third National Health and Nutrition Examination Survey. Arthritis Rheum. 2005;52:283–9.
Rasheed H, Phipps-Green A, Topless R, Hollis-Moffatt JE, Harré Hindmarsh J, Franklin C, et al. Association of the lipoprotein receptor-related protein 2 gene with gout and non-additive interaction with alcohol consumption. Arthritis Res Ther. 2013;15:R177.
Choi HK, McCormick N, Lu N, Rai SK, Yokose C, Zhang Y. Population impact attributable to modifiable risk factors for hyperuricemia. Arthritis Rheumatol. 2020;72:157–65.
Major TJ, Topless RK, Dalbeth N, Merriman TR. Evaluation of the diet wide contribution to serum urate levels: meta-analysis of population based cohorts. BMJ. 2018;363:k3951.
McAdams DeMarco MA, Maynard JW, Baer AN, Gelber AC, Hunter Young J, Alonso A, Coresh J. Diuretic use, increased serum urate levels, and risk of incident gout in a population-based study of adults with hypertension: the Atherosclerosis Risk in Communities cohort study. Arthritis Rheum. 2012;64:121–9.
Savage PJ, Pressel SL, Curb JD, Schron EB, Applegate WB, Black HR, et al. Influence of long-term, low-dose, diuretic-based, antihypertensive therapy on glucose, lipid, uric acid, and potassium levels in older men and women with isolated systolic hypertension: the Systolic Hypertension in the Elderly Program. Arch Int Med. 1998;158:741–51.
Chen-Xu M, Yokose C, Rai SK, Pillinger MH, Choi HK. Contemporary prevalence of gout and hyperuricemia in the United States and decadal trends: the national health and nutrition examination survey, 2007–2016. Arthritis Rheumatol. 2019;71:991–9.
Ioannidis JP. The challenge of reforming nutritional epidemiologic research. JAMA. 2018;320:969–70.
Dalbeth N, Choi HK, Joosten LAB, Khanna PP, Matsuo H, Perez-Ruiz F, Stamp LK. Gout. Nat Rev Dis Prim. 2019;5:69.
Dalbeth N, Phipps-Green A, Frampton C, Neogi T, Taylor WJ, Merriman TR. Relationship between serum urate concentration and clinically evident incident gout: an individual participant data analysis. Ann Rheum Dis. 2018;77:1048–52.
Cox CL, Stanhope KL, Schwarz JM, Graham JL, Hatcher B, Griffen SC, et al. Consumption of fructose-but not glucose-sweetened beverages for 10 weeks increases circulating concentrations of uric acid, retinol binding protein-4, and gamma-glutamyl transferase activity in overweight/obese humans. Nutr Metab. 2012;9:68.
Le MT, Frye RF, Rivard CJ, Cheng J, McFann KK, Segal MS, et al. Effects of high-fructose corn syrup and sucrose on the pharmacokinetics of fructose and acute metabolic and hemodynamic responses in healthy subjects. Metabolism. 2012;61:641–51.
Perez-Pozo S, Schold J, Nakagawa T, Sáchez-Lozada LG, Johnson RJ, López LJ. Excessive fructose intake induces the features of metabolic syndrome in healthy adult men: role of uric acid in the hypertensive response. Int J Obes. 2010;34:454–61.
Larsson SC, Burgess S, Michaëlsson K. Genetic association between adiposity and gout: a Mendelian randomization study. Rheumatology. 2018;57:2145–8.
Oikonen M, Wendelin-Saarenhovi M, Lyytikainen LP, Siitonen N, Loo B-M, Jula A, et al. Associations between serum uric acid and markers of subclinical atherosclerosis in young adults. The cardiovascular risk in Young Finns study. Atherosclerosis. 2012;223:497–503.
Palmer TM, Nordestgaard BG, Benn M, Tybjaerg-Hansen A, Davey Smith G, Lawlor DA, Timpson NJ. Association of plasma uric acid with ischaemic heart disease and blood pressure: Mendelian randomisation analysis of two large cohorts. BMJ. 2013;347:f4262.
Merino J, Dashti HS, Li SX, Sarnowski C, Justice AE, Graff M, et al. Genome-wide meta-analysis of macronutrient intake of 91,114 European ancestry participants from the cohorts for heart and aging research in genomic epidemiology consortium. Mol Psych. 2019;24:1920–32.
Reed DR, Bachmanov AA, Beauchamp GK, Tordoff MG, Price RA. Heritable variation in food preferences and their contribution to obesity. Behav Genet. 1997;27:373–87.
Tanaka T, Ngwa JS, Van Rooij FJ, Carola Zillikens M, Wojczynski MK, Frazier-Wood C, et al. Genome-wide meta-analysis of observational studies shows common genetic variants associated with macronutrient intake. Am J Clin Nutr. 2013;97:1395–402.
Cole JB, Florez JC, Hirschhorn JN. Comprehensive genomic analysis of dietary habits in UK Biobank identifies hundreds of genetic associations. Nat Commun. 2020;11:1467.
Juraschek SP, Gelber AC, Choi HK, Appel LJ, Miller ER. Effects of the dietary approaches to stop hypertension (DASH) diet and sodium intake on serum uric acid. Arthritis Rheumatol. 2016;68:3002–9.
De Lorgeril M, Salen P, Martin J-L, Monjaud I, Delaye J, Mamelle N. Mediterranean diet, traditional risk factors, and the rate of cardiovascular complications after myocardial infarction: final report of the Lyon Diet Heart Study. Circulation. 1999;99:779–85.
Wang DD, Leung CW, Li Y, Ding EL, Chiuve SE, Hu FB, Willett WC. Trends in dietary quality among adults in the United States, 1999 through 2010. JAMA Int Med. 2014;174:1587–95.
Zapata-Diomedi B, Barendregt JJ, Veerman JL. Population attributable fraction: names, types and issues with incorrect interpretation of relative risks. Br J Sports Med. 2018;52:212–3.
Cadzow M, Merriman TR, Dalbeth N. Performance of gout definitions for genetic epidemiological studies: analysis of UK Biobank. Arthritis Res Ther. 2017;19:181.
Dalbeth D, Schumacher HR, Fransen J, Neogi T, Jansen TL, Brown M, et al. Survey definitions of gout for epidemiological studies: comparison with crystal identification as the gold standard. Arthritis Care Res. 2016;68:1894–8.
Bycroft C, Freeman C, Petkova D, Band G, Elliott LT, Sharp K, et al. The UK Biobank resource with deep phenotyping and genomic data. Nature. 2018;562:203–9.
Loeb JN. The influence of temperature on the solubility of monosodium urate. Arthritis Rheum. 1972;15:189–92.
Köttgen A, Albrecht E, Teumer A, Vitart V, Krumsiek J, Hundertmark C, et al. Genome-wide association analyses identify 18 new loci associated with serum urate concentrations. Nat Genet. 2013;45:145–54.
Rockhill B, Newman B, Weinberg C. Use and misuse of population attributable fractions. Am J Publ Health. 1998;88:15–9.
Rückinger S, von Kries R, Toschke AM. An illustration of and programs estimating attributable fractions in large scale surveys considering multiple risk factors. BMC Med Res Methodol. 2009;9:7.
Yavorska OO, Burgess S. MendelianRandomization: an R package for performing Mendelian randomization analyses using summarized data. Int J Epidemiol. 2017;46:1734–9.
Burgess S, Butterworth A, Thompson SG. Mendelian randomization analysis with multiple genetic variants using summarized data. Genet Epidemiol. 2013;37:658–65.
Bowden J, Davey Smith G, Burgess S. Mendelian randomization with invalid instruments: effect estimation and bias detection through Egger regression. Int J Epidemiol. 2015;44:512–25.
Bowden J, Davey Smith G, Haycock PC, Burgess S. Consistent estimation in Mendelian randomization with some invalid instruments using a weighted median estimator. Genet Epidemiol. 2016;40:304–14.
Burgess S, Thompson SG. Multivariable Mendelian randomization: the use of pleiotropic genetic variants to estimate causal effects. Am J Epidemiol. 2015;181:251–60.
Morrison J, Knoblauch N, Marcus JH, Stephens M, He X. Mendelian randomization accounting for correlated and uncorrelated pleiotropic effects using genome-wide summary statistics. Nat Genet. 2020;52:740–7.
Li X, Meng X, Timofeeva M, Tzoulaki I, Tsilidis KK, Ioannidis JP, et al. Serum uric acid levels and multiple health outcomes: umbrella review of evidence from observational studies, randomised controlled trials, and Mendelian randomisation studies. BMJ. 2017;357:j2376.
Dalbeth N, Stamp LK, Merriman TR. The genetics of gout: towards personalised medicine? BMC Med. 2017;15:108.
Dalbeth N, Wong S, Gamble GD, Horne A, Mason B, Pool B, et al. Acute effect of milk on serum urate concentrations: a randomised controlled crossover trial. Ann Rheum Dis. 2010;69:1677–82.
Garrel DR, Verdy M, PetitClerc C, Martin C, Brúle D, Hamet P. Milk-and soy-protein ingestion: acute effect on serum uric acid concentration. Am J Clin Nutr. 1991;53:665–9.
Ghadirian P, Shatenstein B, Verdy M, Hamet P. The influence of dairy products on plasma uric acid in women. Eur J Epidemiol. 1995;11:275–81.
Loenen HM, Eshuis H, Löwik MR, Schouten EG, Hulshof KF, Odink J, Kok FJ. Serum uric acid correlates in elderly men and women with special reference to body composition and dietary intake (Dutch nutrition surveillance system). J Clin Epidemiol. 1990;43:1297–303.
Ryu KA, Kang HH, Kim SY, Yoo MK, Seon Kim J, Haw Lee C, Ah WG. Comparison of nutrient intake and diet quality between hyperuricemia subjects and controls in Korea. Clin Nutr Res. 2014;3:56–63.
Zgaga L, Theodoratou E, Kyle J, Farrington AM, Agakov F, Tenesa A, et al. The association of dietary intake of purine-rich vegetables, sugar-sweetened beverages and dairy with plasma urate, in a cross-sectional study. PLoS One. 2012;7:e38123.
Zykova SN, Storhaug HM, Toft I, Chadban SJ, Jenssen TG, White SL. Cross-sectional analysis of nutrition and serum uric acid in two Caucasian cohorts: the AusDiab Study and the Tromsø study. Nutr J. 2015;14:49.
Nakayama A, Matsuo H, Nakaoka H, Nakamura T, Nakashima H, Takada Y, et al. Common dysfunctional variants of ABCG2 have stronger impact on hyperuricemia progression than typical environmental risk factors. Sci Rep. 2014;4:5227.
Choi HK, Nguyen U-S, Niu J, Danaei G, Zhang Y. Selection bias in rheumatic disease research. Nat Rev Rheumatol 2014;10:403–412.
Fry A, Littlejohns TJ, Sudlow C, Doherty N, Adamska L, Sprosen T, et al. Comparison of sociodemographic and health-related characteristics of UK Biobank participants with those of the general population. Am J Epidemiol. 2017;186:1026–34.
Rathod-Mistry T, Roddy E, Jordan K, Mallen CD, Blagojevic-Bucknall M. Long-term effect of allopurinol use on serum uric acid levels, mortality and comorbidities in gout patients: an electronic health record study. Rheumatol. 2019;S58:kez105.026.
Dalbeth N, Petrie KJ, House M, Chong J, Leung W, Chegudi R, et al. Illness perceptions in patients with gout and the relationship with progression of musculoskeletal disability. Arthritis Care Res. 2011;63:1605–12.
Duyck SD, Petrie KJ, Dalbeth N. “You don’t have to be a drinker to get gout, but it helps”: a content analysis of the depiction of gout in popular newspapers. Arthritis Care Res. 2016;68:1721–5.
Spaetgens B, Pustjens T, Scheepers LE, Janssens HJEM, van der Linden S, Boonen A. Knowledge, illness perceptions and stated clinical practice behaviour in management of gout: a mixed methods study in general practice. Clin Rheumatol. 2016;35:2053–61.
Spencer K, Carr A, Doherty M. Patient and provider barriers to effective management of gout in general practice: a qualitative study. Ann Rheum Dis. 2012;71:1490–5.
We thank the staff, participants and funding bodies of the ARIC, CARDIA, CHS, FHS and UK Biobank studies for their important contributions. We thank Hyon Choi for insightful discussion and anonymous peer reviewers for their constructive input.
This research has been conducted using the UK Biobank resource under application numbers 12611 and 11898 and Database of Genotypes and Phenotypes under application number 834. The ARIC study is carried out as a collaborative study supported by the National Heart, Lung, and Blood Institute contracts N01-HC-55015, N01-HC-55016, N01-HC-55018, N01-HC-55019, N01-HC-55020, N01-HC-55021, N01-HC-55022, R01HL087641, R01HL59367 and R01HL086694; National Human Genome Research Institute contract U01HG004402; and National Institutes of Health contract HHSN268200625226C. Infrastructure was partly supported by Grant Number UL1RR025005, a component of the National Institutes of Health and NIH Roadmap for Medical Research. The FHS and the Framingham SHARe project are conducted and supported by the National Heart, Lung, and Blood Institute in collaboration with Boston University. The Framingham SHARe data used for the analyses described in this manuscript were obtained through the Database of Genotypes and Phenotypes. The CHS research reported in this article was supported by contract numbers N01-HC-85079, N01-HC-85080, N01-HC-85081, N01-HC-85082, N01-HC-85083, N01-HC-85084, N01-HC-85085, N01-HC-85086, N01-HC-35129, N01 HC-15103, N01 HC-55222, N01-HC-75150, N01-HC-45133, N01-HC-85239 and HHSN268201200036C; grant numbers U01 HL080295 from the National Heart, Lung, and Blood Institute and R01 AG-023629 from the National Institute on Ageing, with additional contribution from the National Institute of Neurological Disorders and Stroke. A full list of principal CHS investigators and institutions can be found at www.chs-nhlbi.org/pi.htm. The Coronary Artery Risk Development in Young Adults Study (CARDIA) is conducted and supported by the National Heart, Lung, and Blood Institute (NHLBI) in collaboration with the University of Alabama at Birmingham (N01-HC95095 and N01-HC48047), University of Minnesota (N01-HC48048), Northwestern University (N01-HC48049) and Kaiser Foundation Research Institute (N01-HC48050).
This manuscript was not prepared in collaboration with, nor approved by, investigators of the ARIC, CARDIA, CHS or FHS studies and does not necessarily reflect the opinions or views of these studies or their institutions or funding bodies.
This study was funded by the Health Research Council of New Zealand. The funder had no role in the design, execution and reporting of the study.
Ethics approval and consent to participate
Ethical approval for the UK Biobank participants was obtained from the North West Multi-Centre Research Ethics Committee (11/NW/0382). All participants gave written informed consent.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Sex-stratified population attributable and average attributable fractions for environmental and endogenous risk exposures for hyperuricaemia. Table S2. Population attributable and average attributable fractions for 30 genetic risk exposures for hyperuricaemia. Table S3. Causal effects of dietary habits on urate levels: Mendelian randomisation results. Table S4. Genetic variants comprising the four dietary habit instrumental variables.
Edited averageAF.R script for average attributable fraction confidence interval calculation. Figure S1. FISH plot of primary Mendelian randomisation analysis vs analysis including all ambiguous SNPs. Points in the upper right quadrant represent those with consistent effect directions and those in the lower left quadrant represent those with inconsistent effect directions. Nearly all of the points, especially the most significant points, are on the 45 degree line in the upper right quadrant.
About this article
Cite this article
Topless, R.K.G., Major, T.J., Florez, J.C. et al. The comparative effect of exposure to various risk factors on the risk of hyperuricaemia: diet has a weak causal effect. Arthritis Res Ther 23, 75 (2021). https://doi.org/10.1186/s13075-021-02444-8
- Risk factor
- Population attributable fraction
- Genetic polymorphism
- Mendelian randomisation