Multiple Genetic Loci Influence Serum Urate and Their Relationship with Gout and Cardiovascular Disease Risk Factors

doi:10.1161/CIRCGENETICS.109.934455

Circ Cardiovasc Genet. Author manuscript; available in PMC 2012 June 11.

Published in final edited form as:

Circ Cardiovasc Genet. 2010 December; 3(6): 523–530.

Published online 2010 September 30. doi: 10.1161/CIRCGENETICS.109.934455

PMCID: PMC3371395

NIHMSID: NIHMS272945

Multiple Genetic Loci Influence Serum Urate and Their Relationship with Gout and Cardiovascular Disease Risk Factors

Qiong Yang, PhD,^*^1,² Anna Köttgen, MD MPH,^*³ Abbas Dehghan, MD DSc,^*⁴ Albert V. Smith, PhD,^*⁵ Nicole L. Glazer, PhD,^*⁶ Huei Chen, PhD,^2,¹⁸ Daniel I. Chasman, PhD,²⁷ Thor Aspelund, PhD,^5,⁷ Gudny Eiriksdottir, MSc,⁵ Tamara B. Harris, MD,⁸ Lenore Launer, PhD,⁸ Michael Nalls, PhD,²⁸ Dena Hernandez, MS,²⁸ Dan E Arking, PhD,¹⁰ Eric Boerwinkle, PhD,⁹ Megan L. Grove, MS,¹¹ Man Li, MS,³ WH Linda Kao, PhD MHS,^3,¹² Michel Chonchol, MD,¹³ Talin Haritunians, PhD,¹⁴ Guo Li, MS,⁶ Thomas Lumley, PhD,¹⁵ Bruce M. Psaty, MD PhD,¹⁶ Michael Shlipak, MD MPH,¹⁷ Shih-Jen Hwang, PhD,^2,²⁰ Martin G. Larson, ScD,^2,¹⁹ Christopher J. O’Donnell, MD MPH,^2,²¹ Ashish Upadhyay, MD,²² Cornelia M. van Duijn, PhD,⁴ Albert Hofman, MD, PhD,⁴ Fernando Rivadeneira, MD PhD,²³ Bruno Stricker, MB PhD,⁴ Andre G. Uitterlinden, PhD,²³ Guillaume Paré, MD, MSc,²⁷ Alex N. Parker, PhD,²⁹ Paul M Ridker, MD,²⁷ David S. Siscovick, MD,^**²⁴ Vilmundur Gudnason, MD PhD,^**^5,⁷ Jacqueline C. Witteman, PhD,^**⁴ Caroline S. Fox, MD MPH,^**^2,^20,²⁵ and Josef Coresh, MD PhD^**^3,^12,²⁶

¹Dept of Biostatistics, Boston Univ School of Public Health, Boston, MA

²NHLBI’s Framingham Heart Study, Framingham, MA

³Dept of Epidemiology, Johns Hopkins Univ, Baltimore, MD

⁴Dept of Epidemiology, Erasmus Med Ctr, Rotterdam, The Netherlands; Member of the Netherlands Consortium on Healthy Aging (NCHA)

⁵Icelandic Heart Association, Kopavogur, Iceland

⁶Cardiovascular Health Research Unit & Dept of Med, Univ of Washington, Seattle, WA

⁷Univ of Iceland, Reykjavik, Iceland

⁸Lab of Epidemiology, Demography, & Biometry, Intramural Research Program, Nat Institute on Aging

⁹Human Genetics Ctr & Division of Epidemiology, Univ of Texas Health Science Ctr at Houston, Houston, TX

¹⁰McKusick-Nathans Inst of Genetic Med, Johns Hopkins Med Inst, Baltimore, MD

¹¹Human Genetics Ctr, Univ of Texas Health Science Ctr at Houston, Houston, TX

¹²Dept of Med, Johns Hopkins Univ, Baltimore, MD

¹³Univ of Colorado Denver Health Sciences Ctr, Division of Renal Diseases & Hypertension; Aurora, CO

¹⁴Medical Genetics Inst, Cedars-Sinai Med Ctr, Los Angeles, CA

¹⁵Dept of Biostatistics, Univ of Washington, Seattle, WA

¹⁶Cardiovascular Health Research Unit & Dept of Med, Epidemiology and Health Services, Univ of Washington, Seattle, WA, Center for Health Studies, Group Health, Seattle, WA

¹⁷General Internal Med Division; San Francisco VA Med Ctr; Univ of California, San Francisco, CA

¹⁸Dept of Neurology, Boston Univ, Boston, MA

¹⁹Dept of Mathematics & Statistics, Boston Univ, Boston, MA

²⁰The Ctr for Population Studies, NHLBI, Bethesda MD

²¹Cardiology Division, Massachusetts General Hosp, Boston, MA

²²Renal Section, Boston Med Ctr & Boston Univ School of Med, Boston, MA

²³Dept of Internal Med, Erasmus Med Ctr, Rotterdam, The Netherlands; Member of the Netherlands Consortium on Healthy Aging (NCHA)

²⁴Cardiovascular Health Research Unit & Dept of Med & Epidemiology, Univ of Washington, Seattle, WA

²⁵Division of Endocrinology, Hypertension, & Metabolism, Brigham and Women’s Hosp & Harvard Med School, Boston, MA

²⁶Dept of Biostatistics, Johns Hopkins Univ, Baltimore, MD

²⁷Division of Preventive Med, Brigham and Women’s Hospital, Boston, MA

²⁸Lab of Neurogenetics, Intramural Research Program, Nat Institute on Aging

²⁹Amgen Inc, Cambridge MA

Author for Correspondence: Josef Coresh, MD, PhD Professor of Epidemiology, Biostatistics & Medicine Johns Hopkins University 2024 E. Monument, Suite 2-630 Baltimore, MD 21287 Phone 410-955-0495, Fax 410-955-0476 ; Email: coresh/at/jhu.edu

^*These authors have contributed equally to this work

^**These authors contributed jointly to this work

The publisher's final edited version of this article is available free at Circ Cardiovasc Genet

See other articles in PMC that cite the published article.

Abstract

Background

Elevated serum urate levels can lead to gout and are associated with cardiovascular risk factors. We performed genome-wide association to search for genetic susceptibility loci for serum urate and gout, and investigated the causal nature of the associations of serum urate with gout and selected cardiovascular risk factors and coronary heart disease (CHD).

Methods and Results

Meta-analyses of genome-wide association studies (GWAS) were performed in 5 population-based cohorts of the CHARGE consortium for serum urate and gout in 28,283 white individuals. The effect of the most significant SNP at all genome-wide significant loci on serum urate was added to create a genetic urate score. Findings were replicated in the Women’s Genome Health Study (WGHS; n=22,054). SNPs at 8 genetic loci achieved genome-wide significance with serum urate levels (p-values 4×10⁻⁸ to 2×10⁻²⁴²; SLC22A11, GCKR, R3HDM2-INHBC region, RREB1, PDZK1, SLC2A9, ABCG2, SLC17A1). Only two loci [SLC2A9, ABCG2] showed genome-wide significant association with gout. The genetic urate score was strongly associated with serum urate and gout (odds ratio 12.4 per 100 umol/L; p-value=3×10⁻³⁹), but not with blood pressure, glucose, eGFR, chronic kidney disease, or CHD. The lack of association between the genetic score and the latter phenotypes was also observed in WGHS.

Conclusions

The genetic urate score analysis suggested a causal relationship between serum urate and gout but did not provide evidence for one between serum urate and cardiovascular risk factors and CHD.

Keywords: urate, gout, cardiovascular disease risk factors, genome-wide association study, Mendelian randomization

Introduction

Hyperuricemia is a key risk factor for gout, a common inflammatory arthritis caused by deposition of mono-sodium urate crystals in the joints and surrounding tissues.¹ Multiple renal transporters contribute to the maintenance of normal serum urate levels, but the identity and regulators of these transporters are incompletely understood.

While the role of serum urate in the causal pathway for gout has been well-characterized, substantial controversy exists regarding whether elevated serum urate may also be a cause of cardiovascular disease (CVD) risk factors or CVD, or if the association with urate observed in observational studies merely is a consequence of these conditions or an artifact of uncontrolled confounding factors.² A recent randomized clinical trial of allopurinol in adolescents with newly diagnosed hypertension suggested that drugs that affect serum urate levels may also reduce blood pressure, further highlighting the potential role of serum urate in the pathogenesis of hypertension.³

We now expand on work by our consortium⁴ and others^5-8 with GWAS and meta-analysis to identify genetic loci associated with urate and gout in a large sample of 28,283 individuals from five population-based cohorts, the CHARGE (Cohorts for Heart and Aging Research in Genome Epidemiology) consortium.⁹ Capitalizing on Mendelian randomization, we then create a new genetic urate score of susceptibility to hyperuricemia, and examine its association with gout, blood pressure, glucose, chronic kidney disease (CKD), and coronary heart disease (CHD) to improve causal inference from observational studies. Finally, genetic score association findings are replicated in a second, large cohort, the Women’s Genome Health Study (WGHS).

Methods

Five large population- or community-based cohorts from the CHARGE consortium including the Atherosclerosis Risk in Communities Study (ARIC), Cardiovascular Health Study (CHS), and Framingham Heart Study (FHS), the Rotterdam Study (including RS-I and RS-II), and Age Gene/Environment Susceptibility Reykjavik Study (AGES) constituted the discovery cohorts to identify genetic variants associated with urate levels and gout and for the genetic urate score analyses. An additional cohort, Women’s Health Genome Study (WHGS), was used as a replication cohort to confirm the genetic urate score findings.

Each participant of the CHARGE cohorts provided written informed consent, and study protocols were approved by the Institutional Review Boards of their respective institutions. All members of the WGHS cohort were participants in the Women’s Health Study who provided an adequate baseline blood sample for plasma and DNA analysis and who gave consent for blood-based analyses and long-term follow-up; the study has been approved by the Institutional Review Board of the Brigham and Women’s Hospital.

Additional details pertaining to the study samples, including participant recruitment, phenotype definition, genotyping, imputation and data quality control are provided in the Supplementary Methods and Supplementary Table 1s.

CHD, CVD Risk Factors and Covariates in CHARGE Cohorts

CHD and CVD Risk Factors

Systolic and diastolic blood pressure (SBP, DBP) were measured using a sphygmomanometer in the seated position following a standard protocol in each study. For subjects treated with antihypertensive medications, 10 mm Hg were added to their SBP values and 5 mm Hg to DBP values to impute the untreated blood pressure.¹⁰Hypertension was defined as SBP >=140 or DBP >=90 mm Hg or on antihypertensive treatment. Glucose was obtained after an overnight fast. Diabetes was defined by fasting plasma glucose of at least 126 mg/dl, non-fasting plasma glucose of at least 200 mg/dl, diabetes medication, or self report. Fasting insulin was measured as described in the supplement. Serum creatinine was determined using a modified kinetic Jaffe reaction in all studies but AGES, who used an enzymatic method. Creatinine values were calibrated statistically, and estimated glomerular filtration rate (eGFR) was calculated using the four-variable Modification of Diet in Renal Disease (MDRD) Study equation.¹¹ CKD was defined as eGFR <60 ml/min/1.73 m².¹² CHD was defined as incident myocardial infarction (AGES, ARIC, CHS, FHS, RS), fatal CHD (ARIC, CHS, FHS, RS), silent MI from EKG (ARIC), sudden cardiac arrest (CHS), percutaneous coronary interventions (AGES), coronary artery bypass grafting (AGES, RS), coronary revascularization procedure (ARIC) or percutaneous transluminal coronary angioplasty (RS).

Covariate Measurements

Body mass index was calculated as weight in kilograms divided by squared height in meters; alcohol consumption was recorded as drinks per week and converted to grams/week, and antihypertensive treatment was defined as self-reported intake of any antihypertensive medication.

Statistical Analyses in CHARGE Cohorts

GWAS and meta analyses

Separate GWAS analyses of ~ 2.5 million SNPs were first performed within each CHARGE cohort for urate levels adjusting for age, sex, BMI, alcohol, antihypertensive treatment and study-site or cohort and potential population admixture (Supplementary Methods and Table 1s). Fixed-effects meta-analyses were then used to combine the results from all CHARGE cohorts using the software METAL (http://www.sph.umich.edu/csg/abecasis/metal/index.html). The overall effect size of a SNP from the meta-analyses was defined as the inverse variance weighted average of the beta coefficients for the same allele in each cohort. A threshold of 5×10⁻⁸ defined genome-wide significance, corresponding to a Bonferroni adjusted alpha = 0.05 for 1 million independent tests, the estimated multiple testing burden for a GWAS of individuals of European ancestry.¹³ For gout meta-analyses, SNPs with minor allele frequency (MAF) <0.02 were excluded from the analyses to minimize the risk of false positive results.

Genetic urate score analysis

To model the cumulative effects of the loci identified from the urate meta-analysis, we multiplied, for each locus, the number of minor alleles of the most significant SNP each person carried (0-2) by the beta coefficient from the meta-analysis, and added the results to calculate a genetic urate score. Specifically, the genetic urate score equals rs1967017(T)x3.3+rs780093(T)x5.2-rs13129697(G)x22.2+rs2199936(A)x18.1+rs675209(T)x4.4-rs1165196(G)x6.2+rs2078267(C)x6.8-rs1106766(T)x5.2, where each SNP ID followed by the allele in parenthesis denotes the number of copies of that allele carried by an individual and the number multiplied to is the effect size per copy of the allele in μmol/l estimated in our meta analyses. This genetic score calculation is the first step of a two-step least squares instrumental variables analysis to estimate the effect of a genetic score on the intermediary variable (serum urate) in a Mendelian randomization design.¹⁴ The score assigns zero to individuals who carry all major alleles and has the unit of umol/L. Thus, the genetic urate score is the difference in mean serum urate levels of individuals with a specific combination of genotypes compared to individuals who carry two copies of the major alleles at all the susceptibility loci. The genetic urate risk score was then tested for association with CHD and CVD risk factors, with limited adjustment (for age, sex and center) and more extensive adjustment for major confounders.² Meta-analyzed associations are presented per 100 umol/L (1.68 mg/dl) increase in the genetic urate score, since it is similar to the range of the genetic urate score. Since power is a concern for using this approach, we calculated 80% power for each trait for a one sided test at alpha 0.05. Details of the power calculation are presented in the Supplementary Methods.

Replication of Genetic Urate Score Findings in WGHS

An independent cohort of 22,054 from the Women’s Genome Health Study (WGHS) was used to validate the genetic urate score findings from the CHARGE cohorts. Details pertaining to the WGHS are included in the Supplementary Materials and Supplemental Table 1.

Results

Characteristics of individual CHARGE cohorts are presented in Table 1 for the 28,283 participants with serum urate levels.

Table 1

Characteristics of study participants in CHARGE cohorts, mean (SD) or % (n).

Urate and Gout GWAS Meta-analysis in CHARGE Cohorts

Eight loci (SLC22A11, GCKR, R3HDM2-INHBC region, RREB1, PDZK1, SLC2A9, ABCG2, SLC17A1) showed genome-wide significant association in the meta-analysis for urate (meta analysis p-value range 3.5×10⁻⁸ to 1.5×10⁻²⁴²; Table 2). A genome-wide overview of the p-values from the meta-analyses is presented in Supplementary Figure 1s. Supplementary Figure 2s presents regional association plots for each of the eight loci. Details of study specific GWAS results for the eight loci are presented in Supplementary Table 2s and in a forest plot in Supplementary Figure 3s that demonstrates consistent results across studies. Imputation scores for all loci are shown in Supplemental Table 3s, indicating generally good imputation quality across cohorts. All the SNPs with meta-analysis p-value <4×10⁻⁷ are presented in Supplementary Table 6s and 7s for serum urate and gout, respectively.

Table 2

Most significant SNPs in 8 loci associated with serum urate or gout based on meta analysis of GWAS within individual CHARGE cohorts.

Of the loci identified, 6 were also nominally associated with gout (Table 2). Besides SLC2A9 and ABCG2, a separate GWAS of gout did not identify additional loci reaching genome-wide significance.

Genetic Urate Score Association with Serum Urate and Gout in CHARGE Cohorts

The genetic urate risk score was strongly and linearly associated with serum urate (Figure 1; multivariable adjusted p-value <4.5×10⁻³⁰⁸). The score ranged from −67 to +76 umol/L and explained an average of 6.0% of serum urate variance, compared to an average of 0.8% for the individual SNPs and 3.7% for the strongest SNP (SLC2A9). The gout prevalence varied 10-fold across genetic urate score categories (Figure 1; multivariable adjusted p-value 5.5×10⁻³⁴). Across genes, the odds ratio for gout was highly correlated with the effect size for serum urate for each additional minor allele (Figure 2; r=0.97). Overall, for a 100 umol/L higher genetic urate score, average serum urate was 99.3 umol/L higher and gout odds ratios were 12.4 (95% CI 8.5-18.0) times higher (Table 3).

Figure 1

Mean urate levels and prevalence of gout across genetic urate score

Figure 2

Relationship between urate beta coefficients and gout odds ratios across all 8 loci

Table 3

Comparison of the association of serum urate and genetic urate score (per 100 umol/L) with serum urate, gout and cardiovascular risk factors in CHARGE cohorts.

Comparison of Serum Urate and Genetic Urate Score Associations with CVD Risk Factors and CHD in CHARGE Cohorts

Serum urate adjusted for age, sex and center was strongly associated with all of the cardiovascular disease risk factors examined (Table 3, p<0.001). Further covariate adjustment reduced the magnitude of the associations but all associations remained highly significant (p<0.001) except for fasting glucose. In contrast, the genetic urate risk score was not associated (p≥0.05) with SBP, DBP, fasting glucose levels, fasting insulin, eGFR, CKD (n=3092 cases) or incident CHD (n=3050 cases). Excluding individuals on hypertension treatment from the SBP and DBP analyses or individuals with diabetes from the glucose analysis made little difference. We evaluated the detectable genetic urate score effect size per 100 umol/L with 80% power (Table 3, last column) in the context of observed effect sizes of serum urate per 100 umol/L on each phenotype (Table 3, column 3 and 5). For gout, fasting glucose, SBP, DBP, and eGFR, our study had good power as the observed serum urate effect sizes are generally higher than the detectable genetic urate effect sizes. For CKD, although the observed effect size is higher than detectable effect size, the detectable effect size (OR=1.38) is considered high for a genetic association. For CHD, observed serum urate effect size is smaller than detectable genetic score effect size with 80% power. Therefore a modest association with CKD and CHD cannot be ruled out due to insufficient power.

Replication in the Women’s Genome Health Study (WGHS)

Analyses of the association of the urate risk score with CHD and CVD risk factors in the Women’s Genome Health Study demonstrated similar null results (Supplementary Table 4s). Per 100 umol/L higher genetic urate score, the age and sex adjusted odds ratio of hypertension incidence was 0.93 (95% CI 0.78-1.11, p=0.43); coefficients for self reported blood pressure, diabetes, eGFR, and CKD were non-significant. The age and sex-adjusted hazard ratio for incidence of CHD was 0.83 (95% CI 0.53-1.30; p=0.42). Similar results were obtained for multivariable adjusted analyses.

Secondary Analyses Results

We additionally evaluated the association of individual SNPs with cardiovascular risk factors and CHD to examine possible pleiotropy effects, examined a possible threshold effect of the genetic urate score, and conducted sensitivity analyses for the genetic urate score. Description of these analyses and results is included in the Supplementary Materials and Table 5s.

Discussion

We have identified genetic variants in 8 genetic loci associated with serum urate. Our genetic urate score was strongly associated with a 10-fold variation in gout prevalence. In marked contrast, we did not observe an association between the genetic urate score and blood pressure, fasting glucose, CKD or CHD. This lack of association was also observed in another independent large cohort.

Our GWAS results are similar to another recent GWAS with more than 28,000 European individuals.⁸ Six (PDZK1, GCKR, SLC2A9, ABCG2, SLC17A1, SLC22A11) of the nine loci identified by Kolz et al⁸ were also identified by the present study. Among the remaining three loci identified by Kolz et al⁸, SLC16A9 and LRRC16 were borderline genome-wide significant (p=1.7×10⁻⁷ and 1.3×10⁻⁷) in the present study; SLC22A12 reached significance in our study but was not considered as an independent locus because of its close proximity with SLC22A11 and the uncertainty that the signals represent independent variants. We additionally identified R3HDM2-INHBC and RREB1.

We confirmed urate transporters ABCG2, SLC2A9 and SLC17A1 identified in our previous GWAS that included a subset of current study cohorts and a subsequent functional study.^4;15 Among our additional findings, SLC22A11 encodes a known renal urate transporter, OAT4,¹⁶ but it is also in close proximity to SCL22A12, which encodes another well known urate transporter, URAT1¹⁷ PDZK1 has been shown in experimental settings to directly interact with protein products of SLC22A11 and SLC17A1, and has been proposed as a regulator of urate transport activities.¹⁸ The biologic mechanisms underlying the association of SNPs in RREB1, the R3HDM2-INHBC region, and GCKR with serum urate are unknown. RREB1 is a zinc finger transcription factor reported to regulate the androgen receptor and the calcitonin gene.^17;19 INHBC is a member of the TGF-beta superfamily;²⁰ the top associated SNP rs4760254 maps to the intergenic region between INHBC and R3HDM2 genes. SNPs in GCKR were previously reported in association with higher CRP, higher triglyceride levels, and lower glucose levels, suggesting pleiotropic mechanisms.²¹

Multiple lines of evidence link serum urate levels to hypertension, diabetes, CKD and CHD.^2;3 However, epidemiologic studies of metabolic traits are limited in determining causation, and serum urate can be elevated secondary to disease or confounders. Because the association between genes and disease is not generally subject to confounding by environmental factors or reverse causality, causal inferences between exposure and disease can be examined more specifically using Mendelian randomization. ^14;22

The percentage of variation explained by the genetic score is moderate. However, power calculation taking into account the variation explained by the genetic score indicated that this study had enough power to detect genetic score effect sizes similar to observed epidemiology effects of urate for all phenotypes but CHD. For the initial analyses only involving CHARGE cohorts, the genetic score has been constructed and tested for association with cardiovascular risk factors in the same individuals. As a result, the genetic score effects may be overestimated and should be interpreted with caution. However, the potential of overestimating effect sizes of the genetic score is contrary to its lack of association with cardiovascular risk factors except gout in CHARGE. Similar lack-of-association findings have been reported by other studies that examined the association between individual genetic variants for serum urate and cardiovascular risk factors including blood pressure, hypertension, glucose, lipids, and coronary artery diseases.^6;23;24 Taken together, these findings raise questions about whether the observed epidemiologic associations between serum urate and CVD risk factors present in our cohorts and previous studies are due to causal associations.

Limitations of our study pertaining to the generalizability of our findings are: our study focused on CVD risk factors and incident CHD in middle-aged and older individuals, and therefore could not investigate the association of serum urate levels with blood pressure and new onset cardiovascular disease and its risk factors in adolescents and young adults. Participants of this study are white; findings may not be generalizable to other ethnicities. Gout ascertainment was based on self-report or medication records which may have resulted in inaccurate disease classification for some individuals and reduced power to detect an association. However the two loci identified associated with gout showed consistent effects across cohorts and were also the ones identified by other independent studies^5-8;25, demonstrating the robustness of our results to the gout ascertainment. We only tested the association of genetic urate score with CKD defined as eGFR < 60 ml/min/1.73 m². We were not able to examine more extreme definitions of CKD due to lack of power. We were also not able to adjust for the effects of antihyperuricemic treatment (e.g. allopurinol) in our analyses since information on such treatment was not available in all cohorts. However ignoring antihyperuricemic treatment is unlikely to have had a substantial impact on our results since only a small percentage of the study sample is expected to have been treated. Finally, we modeled a linear association of the genetic score with outcomes. We observed no association when the data were stratified by serum urate at the median. However, we cannot rule out the possibility of a threshold effect where low serum urate levels are protective but higher levels do not show progressively more deleterious cardiovascular effects. Despite these limitations, our study was large and well powered. The findings were consistent across all participating studies, and an independent replication sample.

In conclusion, we identified eight genetic loci associated with serum urate in genome-wide association and meta-analyses. Consistent with the association of serum urate and gout, the combined effect of all the genes summarized in a urate genetic risk score is associated with a ten-fold increase of gout risk. Conversely, the urate genetic risk score was not significantly associated with cardiovascular risk factors or CHD, which stands in contrast to the association between serum urate and vascular risk factors observed in epidemiologic studies.

While the role of serum urate in the causal pathway for gout has been well-characterized, substantial controversy exists regarding whether elevated serum urate may also be a cause of high blood pressure (BP), hyperglycemia, and chronic kidney disease (CKD), or if the association with serum urate observed in observational studies merely is a consequence of these conditions or if it is an artifact of uncontrolled confounding factors. Because the association between genes and disease is not generally subject to confounding by environmental factors or reverse causality, causal inferences between exposure and disease can be examined more specifically using Mendelian randomization. In the present investigation, we tested the association of a genetic score constructed using 8 loci associated with serum urate with coronary heart disease (CHD) and its risk factors including gout, glucose, systolic BP, diastolic BP, estimated glomerular filtration rate and CKD. Except for gout, none of the associations was statistically significant, and the lack of associations was replicated in another equally large independent sample. Although further confirmation is warranted, our study helps elucidate the association of serum urate with CHD and its risk factors, which may contribute to a better understanding of the usefulness of controlling serum urate for preventing and designing treatments for CHD and its risk factors.

Supplementary Material

Supp1

Click here to view.^{(1.3M, pdf)}

Acknowledgements

The authors acknowledge the essential role of the CHARGE (Cohorts for Heart and Aging Research in Genome Epidemiology) Consortium in development and support of this manuscript. CHARGE members include the Netherland’s Rotterdam Study (RS), the NHLBI’s Framingham Heart Study (FHS), Cardiovascular Health Study (CHS), the NHLBI’s Atherosclerosis Risk in Communities (ARIC) Study, and the Icelandic Heart Association’s and NIA’s Iceland Age, Gene/Environment Susceptibility (AGES) Reykjavik Study. This research was supported in part by the Intramural Research Program of the NIH, National Institute on Aging We are indebted to the staff and participants of the AGES Reykjavik Study, the ARIC Study, the CHS Study, the FHS Study, and the Rotterdam Study for their important contributions. A full list of principal CHS investigators and institutions can be found at http://www.chs-nhlbi.org/pi.htm. We acknowledge the National Heart, Lung, and Blood Institute, who has made the SHARe (SNP Health Association Resource) project possible. A portion of FHS computations were using the Linux Cluster for Genetic Analysis (LinGA-II) funded by the Robert Dawson Evans Endowment of the Department of Medicine at Boston University School of Medicine and Boston Medical Center. We thank Pascal Arp, Mila Jhamai, Dr Michael Moorhouse, Marijn Verkerk and Sander Bervoets for their help in creating the Rotterdam database and Maxim Struchalin for his contributions to the imputations of the Rotterdam data.

Funding Sources AGES: The Age, Gene/Environment Susceptibility Reykjavik Study has been funded by NIH contract N01-AG-12100, the National Institute on Aging Intramural Research Program, Hjartavernd (the Icelandic Heart Association), and the Althingi (the Icelandic Parliament). ARIC: The Atherosclerosis Risk in Communities Study (ARIC) is carried out as a collaborative study supported by 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. JC, AK and ML by R01DK076770 and DA by a Donald W. Reynolds Clinical Cardiovascular Research Center grant. WHLK was supported by K01DK067207 and JC by R01DK076770. The project described was supported by Grant Number UL1 RR 025005 from the National Center for Research Resources (NCRR), a component of the National Institutes of Health (NIH) and NIH Roadmap for Medical Research, and its contents are solely the responsibility of the authors and do not necessarily represent the official view of NCRR or NIH. Information on NCRR is available at http://www.ncrr.nih.gov/ . Information on Re-engineering the Clinical Research Enterprise can be obtained from http://nihroadmap.nih.gov/clinicalresearch/overview-translational.asp. CHS: The Cardiovascular Health Study (CHS) was supported by contract numbers N01-HC-85079 through N01-HC-85086, N01-HC-35129, N01-HC-15103, N01 HC-55222, N01-HC-75150, N01-HC-45133, grant numbers U01 HL080295, R01 HL087652, and R01 AG027002 from the National Heart, Lung, and Blood Institute. Genotyping was supported, in part, by NCCR grant M01RR00069 and National Institute of Diabetes and Digestive and Kidney Diseases grant DK063491. FHS: The Framingham Heart Study (FHS) was supported by the National Heart, Lung and Blood Institute’s Framingham Heart Study (N01-HC-25195) and Affymetrix genotyping was supported by contract N02-HL-6-4278. RS: The Rotterdam Study (RS) is supported by the Erasmus Medical Center and Erasmus University Rotterdam; the Netherlands Organization for Scientific Research; the Netherlands Organization for Health Research and Development (ZonMw); the Research Institute for Diseases in the Elderly; The Netherlands Heart Foundation; the Ministry of Education, Culture and Science; the Ministry of Health Welfare and Sports; the European Commission; and the Municipality of Rotterdam. Support for genotyping was provided by the Netherlands Organization for Scientific Research (NWO) (175.010.2005.011, 911.03.012) and Research Institute for Diseases in the Elderly (RIDE). This study was further supported by the Netherlands Genomics Initiative (NGI)/Netherlands Organisation for Scientific Research (NWO) project nr. 050-060-810. Abbas Dehghan is supported by a grant from Netherlands Organization for Scientific Research (NWO) (vici, 918-76-619). WGHS: The Women’s Genome Health Study (WGHS) was supported by the National Heart, Lung, and Blood Institute (HL 043851) and the National Cancer Institute (CA 047988). Collaborative scientific and genotyping support was provided by Amgen, Inc, and Alex Parker is an employee of Amgen, Inc.

Footnotes

Conflict of Interest: None

This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final citable form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain.

References

1. Choi HK, Mount DB, Reginato AM. Pathogenesis of gout. Ann Intern Med. 2005;143:499–516. [PubMed]

2. Feig DI, Kang DH, Johnson RJ. Uric acid and cardiovascular risk. N Engl J Med. 2008;359:1811–1821. [PMC free article] [PubMed]

3. Feig DI, Soletsky B, Johnson RJ. Effect of allopurinol on blood pressure of adolescents with newly diagnosed essential hypertension: a randomized trial. JAMA. 2008;300:924–932. [PMC free article] [PubMed]

4. Dehghan A, Kottgen A, Yang Q, Hwang SJ, Kao WL, Rivadeneira F, Boerwinkle E, Levy D, Hofman A, Astor BC, Benjamin EJ, van Duijn CM, Witteman JC, Coresh J, Fox CS. Association of three genetic loci with uric acid concentration and risk of gout: a genome-wide association study. Lancet. 2008 [PMC free article] [PubMed]

5. Li S, Sanna S, Maschio A, Busonero F, Usala G, Mulas A, Lai S, Dei M, Orru M, Albai G, Bandinelli S, Schlessinger D, Lakatta E, Scuteri A, Najjar SS, Guralnik J, Naitza S, Crisponi L, Cao A, Abecasis G, Ferrucci L, Uda M, Chen WM, Nagaraja R. The GLUT9 gene is associated with serum uric acid levels in Sardinia and Chianti cohorts. PLoS Genet. 2007;3:e194. [PubMed]

6. Vitart V, Rudan I, Hayward C, Gray NK, Floyd J, Palmer CN, Knott SA, Kolcic I, Polasek O, Graessler J, Wilson JF, Marinaki A, Riches PL, Shu X, Janicijevic B, Smolej-Narancic N, Gorgoni B, Morgan J, Campbell S, Biloglav Z, Barac-Lauc L, Pericic M, Klaric IM, Zgaga L, Skaric-Juric T, Wild SH, Richardson WA, Hohenstein P, Kimber CH, Tenesa A, Donnelly LA, Fairbanks LD, Aringer M, McKeigue PM, Ralston SH, Morris AD, Rudan P, Hastie ND, Campbell H, Wright AF. SLC2A9 is a newly identified urate transporter influencing serum urate concentration, urate excretion and gout. Nat Genet. 2008;40:437–442. [PubMed]

7. Wallace C, Newhouse SJ, Braund P, Zhang F, Tobin M, Falchi M, Ahmadi K, Dobson RJ, Marcano AC, Hajat C, Burton P, Deloukas P, Brown M, Connell JM, Dominiczak A, Lathrop GM, Webster J, Farrall M, Spector T, Samani NJ, Caulfield MJ, Munroe PB. Genome-wide association study identifies genes for biomarkers of cardiovascular disease: serum urate and dyslipidemia. Am J Hum Genet. 2008;82:139–149. [PubMed]

8. Kolz M, Johnson T, Sanna S, Teumer A, Vitart V, Perola M, Mangino M, Albrecht E, Wallace C, Farrall M, Johansson A, Nyholt DR, Aulchenko Y, Beckmann JS, Bergmann S, Bochud M, Brown M, Campbell H, Connell J, Dominiczak A, Homuth G, Lamina C, McCarthy MI, Meitinger T, Mooser V, Munroe P, Nauck M, Peden J, Prokisch H, Salo P, Salomaa V, Samani NJ, Schlessinger D, Uda M, Volker U, Waeber G, Waterworth D, Wang-Sattler R, Wright AF, Adamski J, Whitfield JB, Gyllensten U, Wilson JF, Rudan I, Pramstaller P, Watkins H, Doering A, Wichmann HE, Spector TD, Peltonen L, Volzke H, Nagaraja R, Vollenweider P, Caulfield M, Illig T, Gieger C. Meta-analysis of 28,141 individuals identifies common variants within five new loci that influence uric acid concentrations. PLoS Genet. 2009;5:e1000504. [PMC free article] [PubMed]

9. Psaty BM, O’Donnell CJ, Gudnason V, Lunetta KL, Folsom AR, Rotter JI, Uitterlinden AG, Harris TB, Witteman JCM, Boerwinkle E., on Behalf of the CHARGE Consortium Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) Consortium: Design of Prospective Meta-Analyses of Genome-Wide Association Studies From 5 Cohorts. Circ Cardiovasc Genet. 2009;2:73–80. [PMC free article] [PubMed]

10. Cui JS, Hopper JL, Harrap SB. Antihypertensive treatments obscure familial contributions to blood pressure variation. Hypertension. 2003;41:207–210. [PubMed]

11. Levey AS, Bosch JP, Lewis JB, Greene T, Rogers N, Roth D. A more accurate method to estimate glomerular filtration rate from serum creatinine: a new prediction equation. Modification of Diet in Renal Disease Study Group. Ann Intern Med. 1999;130:461–470. [PubMed]

12. K/DOQI clinical practice guidelines for chronic kidney disease: evaluation, classification, and stratification. Am J Kidney Dis. 2002;39:S1–266. [PubMed]

13. Pe’er I, Yelensky R, Altshuler D, Daly MJ. Estimation of the multiple testing burden for genomewide association studies of nearly all common variants. Genet Epidemiol. 2008;32:381–385. [PubMed]

14. Lawlor DA, Harbord RM, Sterne JAC, Timpson N, Smith G Davey. Mendelian randomization: Using genes as instruments for making causal inferences in epidemiology. Statistics in Medicine. 2008;27:1133–1163. [PubMed]

15. 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–10342. [PubMed]

16. Hagos Y, Stein D, Ugele B, Burckhardt G, Bahn A. Human renal organic anion transporter 4 operates as an asymmetric urate transporter. J Am Soc Nephrol. 2007;18:430–439. [PubMed]

17. Erlandson RA. The enigmatic perineurial cell and its participation in tumors and in tumorlike entities. Ultrastruct Pathol. 1991;15:335–351. [PubMed]

18. Anzai N, Kanai Y, Endou H. New insights into renal transport of urate. Curr Opin Rheumatol. 2007;19:151–157. [PubMed]

19. Thiagalingam A, De Bustros A, Borges M, Jasti R, Compton D, Diamond L, Mabry M, Ball DW, Baylin SB, Nelkin BD. RREB-1, a novel zinc finger protein, is involved in the differentiation response to Ras in human medullary thyroid carcinomas. Mol Cell Biol. 1996;16:5335–5345. [PMC free article] [PubMed]

20. Meise W, Pfisterer H. [The course of the synthesis of 1-methyl-15,16,17,18,19,20-hexadehydroyohimbanes from 1-methyltryptamine and isochroman-3-ones (author’s transl)] Arch Pharm (Weinheim) 1977;310:495–501. [PubMed]

21. Orho-Melander M, Melander O, Guiducci C, Perez-Martinez P, Corella D, Roos C, Tewhey R, Rieder MJ, Hall J, Abecasis G, Tai ES, Welch C, Arnett DK, Lyssenko V, Lindholm E, Saxena R, de Bakker PI, Burtt N, Voight BF, Hirschhorn JN, Tucker KL, Hedner T, Tuomi T, Isomaa B, Eriksson KF, Taskinen MR, Wahlstrand B, Hughes TE, Parnell LD, Lai CQ, Berglund G, Peltonen L, Vartiainen E, Jousilahti P, Havulinna AS, Salomaa V, Nilsson P, Groop L, Altshuler D, Ordovas JM, Kathiresan S. Common missense variant in the glucokinase regulatory protein gene is associated with increased plasma triglyceride and C-reactive protein but lower fasting glucose concentrations. Diabetes. 2008;57:3112–3121. [PMC free article] [PubMed]

22. Smith G Davey, Ebrahim S. ‘Mendelian randomization’: can genetic epidemiology contribute to understanding environmental determinants of disease? Int J Epidemiol. 2003;32:1–22. [PubMed]

23. Caulfield MJ, Munroe PB, O’Neill D, Witkowska K, Charchar FJ, Doblado M, Evans S, Eyheramendy S, Onipinla A, Howard P, Shaw-Hawkins S, Dobson RJ, Wallace C, Newhouse SJ, Brown M, Connell JM, Dominiczak A, Farrall M, Lathrop GM, Samani NJ, Kumari M, Marmot M, Brunner E, Chambers J, Elliott P, Kooner J, Laan M, Org E, Veldre G, Viigimaa M, Cappuccio FP, Ji C, Iacone R, Strazzullo P, Moley KH, Cheeseman C. SLC2A9 is a high-capacity urate transporter in humans. PLoS Med. 2008;5:e197. [PMC free article] [PubMed]

24. Stark K, Reinhard W, Neureuther K, Wiedmann S, Sedlacek K, Baessler A, Fischer M, Weber S, Kaess B, Erdmann J, Schunkert H, Hengstenberg C. Association of common polymorphisms in GLUT9 gene with gout but not with coronary artery disease in a large case-control study. PLoS ONE. 2008;3:e1948. [PMC free article] [PubMed]

25. Yamagishi K, Tanigawa T, Kitamura A, Kottgen A, Folsom AR, Iso H. The rs2231142 variant of the ABCG2 gene is associated with uric acid levels and gout among Japanese people. Rheumatology (Oxford) 2010 [PubMed]