Multivariate analysis of a missense variant in CREBRF reveals associations with measures of adiposity in people of Polynesian ancestries.

The minor allele of rs373863828, a missense variant in CREB3 Regulatory Factor, is associated with several cardiometabolic phenotypes in Polynesian peoples. To better understand the variant, we tested the association of rs373863828 with a panel of correlated phenotypes (body mass index [BMI], weight, height, HDL cholesterol, triglycerides, and total cholesterol) using multivariate Bayesian association and network analyses in a Samoa cohort (n = 1632), Aotearoa New Zealand cohort (n = 1419), and combined cohort (n = 2976). An expanded set of phenotypes (adding estimated fat and fat-free mass, abdominal circumference, hip circumference, and abdominal-hip ratio) was tested in the Samoa cohort (n = 1496). In the Samoa cohort, we observed significant associations (log10 Bayes Factor [BF] ≥ 5.0) between rs373863828 and the overall phenotype panel (8.81), weight (8.30), and BMI (6.42). In the Aotearoa New Zealand cohort, we observed suggestive associations (1.5 < log10 BF < 5) between rs373863828 and the overall phenotype panel (4.60), weight (3.27), and BMI (1.80). In the combined cohort, we observed concordant signals with larger log10 BFs. In the Samoa-specific expanded phenotype analyses, we also observed significant associations between rs373863828 and fat mass (5.65), abdominal circumference (5.34), and hip circumference (5.09). Bayesian networks provided evidence for a direct association of rs373863828 with weight and indirect associations with height and BMI.

A Polynesian-specific copy number variant encompassing the MICA gene associates with gout.

Gout is of particularly high prevalence in the Maori and Pacific (Polynesian) populations of Aotearoa New Zealand (NZ). Here, we investigated the contribution of common population-specific copy number variation (CNV) to gout in the Aotearoa NZ Polynesian population. Microarray-generated genome-wide genotype data from Aotearoa NZ Polynesian individuals with (n = 1196) and without (n = 1249) gout were analyzed. Comparator population groups were 552 individuals of European ancestry and 1962 of Han Chinese ancestry. Levels of circulating major histocompatibility complex (MHC) class I polypeptide-related sequence A (MICA) were measured by enzyme-linked immunosorbent assay. Fifty-four CNV regions (CNVRs) appearing in at least 10 individuals were detected, of which seven common (>2%) CNVRs were specific to or amplified in Polynesian people. A burden test of these seven revealed associations of insertion/deletion with gout (odds ratio (OR) 95% confidence interval [CI] = 1.80 [1.01; 3.22], P = 0.046). Individually testing of the seven CNVRs for association with gout revealed nominal association of CNVR1 with gout in Western Polynesian (Chr6: 31.36-31.45 Mb, OR = 1.72 [1.03; 2.92], P = 0.04), CNVR6 in the meta-analyzed Polynesian sample sets (Chr1: 196.75-196.92 Mb, OR = 1.86 [1.16; 3.00], P = 0.01) and CNVR9 in Western Polynesian (Chr1: 189.35-189.54 Mb, OR = 2.75 [1.15; 7.13], P = 0.03). Analysis of European gout genetic association data demonstrated a signal of association at the CNVR1 locus that was an expression quantitative trait locus for MICA. The most common CNVR (CNVR1) includes deletion of the MICA gene, encoding an immunomodulatory protein. Expression of MICA was reduced in the serum of individuals with the deletion. In summary, we provide evidence for the association of CNVR1 containing MICA with gout in Polynesian people, implicating class I MHC-mediated antigen presentation in gout.

Association of rs9939609 in FTO with BMI among Polynesian peoples living in Aotearoa New Zealand and other Pacific nations.

The fat mass and obesity associated (FTO) locus consistently associates with higher body mass index (BMI) across diverse ancestral groups. However, previous small studies of people of Polynesian ancestries have failed to replicate the association. In this study, we used Bayesian meta-analysis to test rs9939609, the most replicated FTO variant, for association with BMI with a large sample (n = 6095) of Aotearoa New Zealanders of Polynesian (Maori and Pacific) ancestry and of Samoan people living in the Independent State of Samoa and in American Samoa. We did not observe statistically significant association within each separate Polynesian subgroup. Bayesian meta-analysis of the Aotearoa New Zealand Polynesian and Samoan samples resulted in a posterior mean effect size estimate of +0.21 kg/m2, with a 95% credible interval [+0.03 kg/m2, +0.39 kg/m2]. While the Bayes Factor (BF) of 0.77 weakly favors the null, the BF = 1.4 Bayesian support interval is [+0.04, +0.20]. These results suggest that rs9939609 in FTO may have a similar effect on mean BMI in people of Polynesian ancestries as previously observed in other ancestral groups.

Genomic dissection of 43 serum urate-associated loci provides multiple insights into molecular mechanisms of urate control.

High serum urate is a prerequisite for gout and associated with metabolic disease. Genome-wide association studies (GWAS) have reported dozens of loci associated with serum urate control; however, there has been little progress in understanding the molecular basis of the associated loci. Here, we employed trans-ancestral meta-analysis using data from European and East Asian populations to identify 10 new loci for serum urate levels. Genome-wide colocalization with cis-expression quantitative trait loci (eQTL) identified a further five new candidate loci. By cis- and trans-eQTL colocalization analysis, we identified 34 and 20 genes, respectively, where the causal eQTL variant has a high likelihood that it is shared with the serum urate-associated locus. One new locus identified was SLC22A9 that encodes organic anion transporter 7 (OAT7). We demonstrate that OAT7 is a very weak urate-butyrate exchanger. Newly implicated genes identified in the eQTL analysis include those encoding proteins that make up the dystrophin complex, a scaffold for signaling proteins and transporters at the cell membrane; MLXIP that, with the previously identified MLXIPL, is a transcription factor that may regulate serum urate via the pentose-phosphate pathway and MRPS7 and IDH2 that encode proteins necessary for mitochondrial function. Functional fine mapping identified six loci (RREB1, INHBC, HLF, UBE2Q2, SFMBT1 and HNF4G) with colocalized eQTL containing putative causal SNPs. This systematic analysis of serum urate GWAS loci identified candidate causal genes at 24 loci and a network of previously unidentified genes likely involved in control of serum urate levels, further illuminating the molecular mechanisms of urate control.

Mid-pass whole genome sequencing enables biomedical genetic studies of diverse populations.

BACKGROUND: Historically, geneticists have relied on genotyping arrays and imputation to study human genetic variation. However, an underrepresentation of diverse populations has resulted in arrays that poorly capture global genetic variation, and a lack of reference panels. This has contributed to deepening global health disparities. Whole genome sequencing (WGS) better captures genetic variation but remains prohibitively expensive. Thus, we explored WGS at "mid-pass" 1-7x coverage. RESULTS: Here, we developed and benchmarked methods for mid-pass sequencing. When applied to a population without an existing genomic reference panel, 4x mid-pass performed consistently well across ethnicities, with high recall (98%) and precision (97.5%). CONCLUSION: Compared to array data imputed into 1000 Genomes, mid-pass performed better across all metrics and identified novel population-specific variants with potential disease relevance. We hope our work will reduce financial barriers for geneticists from underrepresented populations to characterize their genomes prior to biomedical genetic applications.

Trans-ancestral dissection of urate- and gout-associated major loci SLC2A9 and ABCG2 reveals primate-specific regulatory effects.

Gout is a complex inflammatory arthritis affecting ~20% of people with an elevated serum urate level (hyperuricemia). Gout and hyperuricemia are essentially specific to humans and other higher primates, with varied prevalence across ancestral groups. SLC2A9 and ABCG2 are major loci associated with both urate and gout in multiple ancestral groups. However, fine mapping has been challenging due to extensive linkage disequilibrium underlying the associated regions. We used trans-ancestral fine mapping integrated with primate-specific genomic information to address this challenge. Trans-ancestral meta-analyses of GWAS cohorts of either European (EUR) or East Asian (EAS) ancestry resulted in single-variant resolution mappings for SLC2A9 (rs3775948 for urate and rs4697701 for gout) and ABCG2 (rs2622621 for gout). Tests of colocalization of variants in both urate and gout suggested existence of a shared candidate causal variant for SLC2A9 only in EUR and for ABCG2 only in EAS. The fine-mapped gout variant rs4697701 was within an ancient enhancer, whereas rs2622621 was within a primate-specific transposable element, both supported by functional evidence from the Roadmap Epigenomics project in human primary tissues relevant to urate and gout. Additional primate-specific elements were found near both loci and those adjacent to SLC2A9 overlapped with known statistical epistatic interactions associated with urate as well as multiple super-enhancers identified in urate-relevant tissues. We conclude that by leveraging ancestral differences trans-ancestral fine mapping has identified ancestral and functional variants for SLC2A9 or ABCG2 with primate-specific regulatory effects on urate and gout.

Subtype-specific gout susceptibility loci and enrichment of selection pressure on ABCG2 and ALDH2 identified by subtype genome-wide meta-analyses of clinically defined gout patients.

OBJECTIVES: Genome-wide meta-analyses of clinically defined gout were performed to identify subtype-specific susceptibility loci. Evaluation using selection pressure analysis with these loci was also conducted to investigate genetic risks characteristic of the Japanese population over the last 2000-3000 years. METHODS: Two genome-wide association studies (GWASs) of 3053 clinically defined gout cases and 4554 controls from Japanese males were performed using the Japonica Array and Illumina Array platforms. About 7.2 million single-nucleotide polymorphisms were meta-analysed after imputation. Patients were then divided into four clinical subtypes (the renal underexcretion type, renal overload type, combined type and normal type), and meta-analyses were conducted in the same manner. Selection pressure analyses using singleton density score were also performed on each subtype. RESULTS: In addition to the eight loci we reported previously, two novel loci, PIBF1 and ACSM2B, were identified at a genome-wide significance level (p<5.0×10-8) from a GWAS meta-analysis of all gout patients, and other two novel intergenic loci, CD2-PTGFRN and SLC28A3-NTRK2, from normal type gout patients. Subtype-dependent patterns of Manhattan plots were observed with subtype GWASs of gout patients, indicating that these subtype-specific loci suggest differences in pathophysiology along patients' gout subtypes. Selection pressure analysis revealed significant enrichment of selection pressure on ABCG2 in addition to ALDH2 loci for all subtypes except for normal type gout. CONCLUSIONS: Our findings on subtype GWAS meta-analyses and selection pressure analysis of gout will assist elucidation of the subtype-dependent molecular targets and evolutionary involvement among genotype, phenotype and subtype-specific tailor-made medicine/prevention of gout and hyperuricaemia.

Discordant association of the CREBRF rs373863828 A allele with increased BMI and protection from type 2 diabetes in Maori and Pacific (Polynesian) people living in Aotearoa/New Zealand.

AIMS/HYPOTHESIS: The A (minor) allele of CREBRF rs373863828 has been associated with increased BMI and reduced risk of type 2 diabetes in the Samoan populations of Samoa and American Samoa. Our aim was to test rs373863828 for associations with BMI and the odds of type 2 diabetes, gout and chronic kidney disease (CKD) in Maori and Pacific (Polynesian) people living in Aotearoa/New Zealand. METHODS: Linear and logistic regression models were used to analyse the association of the A allele of CREBRF rs373863828 with BMI, log-transformed BMI, waist circumference, type 2 diabetes, gout and CKD in 2286 adults. The primary analyses were adjusted for age, sex, the first four genome-wide principal components and (where appropriate) BMI, waist circumference and type 2 diabetes. The primary analysis was conducted in ancestrally defined groups and association effects were combined using meta-analysis. RESULTS: For the A allele of rs373863828, the effect size was 0.038 (95% CI 0.022, 0.055, p = 4.8 × 10-6) for log-transformed BMI, with OR 0.59 (95% CI 0.47, 0.73, p = 1.9 × 10-6) for type 2 diabetes. There was no evidence for an association of genotype with variance in BMI (p = 0.13), and nor was there evidence for associations with serum urate (ß = 0.012 mmol/l, pcorrected = 0.10), gout (OR 1.00, p = 0.98) or CKD (OR 0.91, p = 0.59). CONCLUSIONS/INTERPRETATION: Our results in New Zealand Polynesian adults replicate, with very similar effect sizes, the association of the A allele of rs373863828 with higher BMI but lower odds of type 2 diabetes among Samoan adults living in Samoa and American Samoa.

A Polynesian-specific missense CETP variant alters the lipid profile.

Identifying population-specific genetic variants associated with disease and disease-predisposing traits is important to provide insights into the genetic determinants of health and disease between populations, as well as furthering genomic justice. Various common pan-population polymorphisms at CETP associate with serum lipid profiles and cardiovascular disease. Here, sequencing of CETP identified a missense variant rs1597000001 (p.Pro177Leu) specific to Maori and Pacific people that associates with higher HDL-C and lower LDL-C levels. Each copy of the minor allele associated with higher HDL-C by 0.236 mmol/L and lower LDL-C by 0.133 mmol/L. The rs1597000001 effect on HDL-C is comparable with CETP Mendelian loss-of-function mutations that result in CETP deficiency, consistent with our data, which shows that rs1597000001 lowers CETP activity by 27.9%. This study highlights the potential of population-specific genetic analyses for improving equity in genomics and health outcomes for population groups underrepresented in genomic studies.

Association of Gout Polygenic Risk Score With Age at Disease Onset and Tophaceous Disease in European and Polynesian Men With Gout.

OBJECTIVE: To determine whether a gout polygenic risk score (PRS) is associated with age at gout onset and tophaceous disease in European, East Polynesian, and West Polynesian men and women with gout. METHODS: A 19-variant gout PRS was produced in 7 European gout cohorts (N = 4,016), 2 East Polynesian gout cohorts (N = 682), and 1 West Polynesian gout cohort (N = 490). Sex-stratified regression models were used to estimate the relationship between the PRS and age at gout onset and tophaceous disease. RESULTS: The PRS was associated with earlier age at gout onset in men (ß = -3.61 in years per unit PRS [95% confidence interval (95% CI) -4.32, -2.90] in European men; ß = -6.35 [95% CI -8.91, -3.80] in East Polynesian men; ß = -3.51 [95% CI -5.46, -1.57] in West Polynesian men) but not in women (ß = 0.07 [95% CI -2.32, 2.45] in European women; ß = 0.20 [95% CI -7.21, 7.62] in East Polynesian women; ß -3.33 [95% CI -9.28, 2.62] in West Polynesian women). The PRS showed a positive association with tophaceous disease in men (odds ratio [OR] for the association 1.15 [95% CI 1.00, 1.31] in European men; OR 2.60 [95% CI 1.66, 4.06] in East Polynesian men; OR 1.53 [95% CI 1.07, 2.19] in West Polynesian men) but not in women (OR for the association 0.68 [95% CI 0.42, 1.10] in European women; OR 1.45 [95% CI 0.39, 5.36] in East Polynesian women). The PRS association with age at gout onset was robust to the removal of ABCG2 variants from the PRS in European and East Polynesian men (ß = -2.42 [95% CI -3.37, -1.46] and ß = -6.80 [95% CI -10.06, -3.55], respectively) but not in West Polynesian men (ß = -1.79 [95% CI -4.74, 1.16]). CONCLUSION: Genetic risk variants for gout also harbor risk for earlier age at gout onset and tophaceous disease in European and Polynesian men. Our findings suggest that earlier gout onset involves the accumulation of gout risk alleles in men but perhaps not in women, and that this genetic risk is shared across multiple ancestral groups.

CREBRF missense variant rs373863828 has both direct and indirect effects on type 2 diabetes and fasting glucose in Polynesian peoples living in Samoa and Aotearoa New Zealand.

INTRODUCTION: The minor allele of a missense variant, rs373863828, in CREBRF is associated with higher body mass index (BMI), lower fasting glucose, and lower odds of type 2 diabetes. rs373863828 is common in Pacific Island populations (minor allele frequency (MAF) 0.096-0.259) but rare in non-Pacific Island populations (MAF <0.001). We examined the cross-sectional associations between BMI and rs373863828 in type 2 diabetes and fasting glucose with a large sample of adults of Polynesian ancestries from Samoa, American Samoa, and Aotearoa New Zealand, and estimated the direct and indirect (via BMI) effects of rs373863828 on type 2 diabetes and fasting glucose. RESEARCH DESIGN AND METHODS: We regressed type 2 diabetes and fasting glucose on BMI and rs373863828 stratified by obesity, regressed type 2 diabetes and fasting glucose on BMI stratified by rs373863828 genotype, and assessed the effects of rs373863828 on type 2 diabetes and fasting glucose with path analysis. The regression analyses were completed separately in four samples that were recruited during different time periods between 1990 and 2010 and then the results were meta-analyzed. All samples were pooled for the path analysis. RESULTS: Association of BMI with type 2 diabetes and fasting glucose may be greater in those without obesity (OR=7.77, p=0.015 and ß=0.213, p=9.53×10-5, respectively) than in those with obesity (OR=5.01, p=1.12×10-9 and ß=0.162, p=5.63×10-6, respectively). We did not observe evidence of differences in the association of BMI with type 2 diabetes or fasting glucose by genotype. In the path analysis, the minor allele has direct negative (lower odds of type 2 diabetes and fasting glucose) and indirect positive (higher odds of type 2 diabetes and fasting glucose) effects on type 2 diabetes risk and fasting glucose, with the indirect effects mediated through a direct positive effect of rs373863828 on BMI. CONCLUSIONS: There may be a stronger effect of BMI on fasting glucose in Polynesian individuals without obesity than in those with obesity. Carrying the rs373863828 minor allele does not decouple higher BMI from higher odds of type 2 diabetes.

The epigenome: key to understanding and predicting gout flares.

Gout is a form of arthritis, resulting from an inflammatory reaction to the deposition of monosodium urate (MSU) crystals in the synovial fluid of the joint space. It is characterised by periods of acute inflammation in the affected joint, or joints (known as gout flares), separated by asymptomatic periods. There seems to be substantial overlap between environmental triggers of gout flares and common environmental modifiers (diet, pharmaceuticals, and stress) of epigenetic markers (DNA methylation, histone modifications, and ncRNA). Very few studies have looked at whether environment is influencing gout through epigenetic mechanisms. The pathogenesis of gouty inflammation is well understood but understanding the variation of response to hyperuricaemia in terms of gout flare initiation is less well known. In this review, we will examine the potential of epigenomics in understanding how gout flares may occur, both in terms of development of hyperuricaemia and the inflammatory response. Looking at the epigenome and its intersection with lifestyle could help identify new targets and strategies for effective management of gout flares.

The comparative effect of exposure to various risk factors on the risk of hyperuricaemia: diet has a weak causal effect.

BACKGROUND: 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. METHODS: 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. RESULTS: 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. CONCLUSIONS: 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.

Aotearoa New Zealand Maori and Pacific Population-amplified Gout Risk Variants: CLNK Is a Separate Risk Gene at the SLC2A9 Locus.

OBJECTIVE: The Maori and Pacific (Polynesian) population of Aotearoa New Zealand has a high prevalence of gout. Our aim was to identify potentially functional missense genetic variants in candidate inflammatory genes amplified in frequency that may underlie the increased prevalence of gout in Polynesian populations. METHODS: A list of 712 inflammatory disease-related genes was generated. An in silico targeted exome set was extracted from whole genome sequencing data in people with gout of various ancestral groups (Polynesian, European, East Asian; n = 55, 780, 135, respectively) to identify Polynesian-amplified common missense variants (minor allele frequency > 0.05). Candidate functional variants were tested for association with gout by multivariable-adjusted regression analysis in 2528 individuals of Polynesian ancestry. RESULTS: We identified 26 variants common in the Polynesian population and uncommon in the European and East Asian populations. Three of the 26 population-amplified variants were nominally associated with the risk of gout (rs1635712 [KIAA0319], ORmeta = 1.28, Pmeta = 0.03; rs16869924 [CLNK], ORmeta = 1.37, Pmeta = 0.002; rs2070025 [fibrinogen A alpha chain (FGA)], ORmeta = 1.34, Pmeta = 0.02). The CLNK variant, within the established SLC2A9 gout locus, was genetically independent of the association signal at SLC2A9. CONCLUSION: We provide nominal evidence for the existence of population-amplified genetic variants conferring risk of gout in Polynesian populations. Polymorphisms in CLNK have previously been associated with gout in other populations, supporting our evidence for the association of this gene with gout.

Pleiotropic effect of the ABCG2 gene in gout: involvement in serum urate levels and progression from hyperuricemia to gout.

BACKGROUND: The ABCG2 Q141K (rs2231142) and rs10011796 variants associate with hyperuricaemia (HU). The effect size of ABCG2 rs2231142 on urate is ~ 60% that of SLC2A9, yet the effect size on gout is greater. We tested the hypothesis that ABCG2 plays a role in the progression from HU to gout by testing for association of ABCG2 rs2231142 and rs10011796 with gout using HU controls. METHODS: We analysed 1699 European gout cases and 14,350 normouricemic (NU) and HU controls, and 912 New Zealand (NZ) Polynesian (divided into Eastern and Western Polynesian) gout cases and 696 controls. Association testing was performed using logistic and linear regression with multivariate adjusting for confounding variables. RESULTS: In Europeans and Polynesians, the ABCG2 141K (T) allele was associated with gout using HU controls (OR = 1.85, P = 3.8E- 21 and ORmeta = 1.85, P = 1.3E- 03, respectively). There was evidence for an effect of 141K in determining HU in European (OR = 1.56, P = 1.7E- 18) but not in Polynesian (ORmeta = 1.49, P = 0.057). For SLC2A9 rs11942223, the T allele associated with gout in the presence of HU in European (OR = 1.37, P = 4.7E- 06), however significantly weaker than ABCG2 rs2231142 141K (PHet = 0.0023). In Western Polynesian and European, there was epistatic interaction between ABCG2 rs2231142 and rs10011796. Combining the presence of the 141K allele with the rs10011796 CC-genotype increased gout risk, in the presence of HU, 21.5-fold in Western Polynesian (P = 0.009) and 2.6-fold in European (P = 9.9E- 06). The 141K allele of ABCG2 associated with increased gout flare frequency in Polynesian (Pmeta = 2.5E- 03). CONCLUSION: These data are consistent with a role for ABCG2 141K in gout in the presence of established HU.

Evaluation of the diet wide contribution to serum urate levels: meta-analysis of population based cohorts.

OBJECTIVE: To systematically test dietary components for association with serum urate levels and to evaluate the relative contributions of estimates of diet pattern and inherited genetic variants to population variance in serum urate levels. DESIGN: Meta-analysis of cross sectional data from the United States. DATA SOURCES: Five cohort studies. REVIEW METHODS: 16 760 individuals of European ancestry (8414 men and 8346 women) from the US were included in analyses. Eligible individuals were aged over 18, without kidney disease or gout, and not taking urate lowering or diuretic drugs. All participants had serum urate measurements, dietary survey data, information on potential confounders (sex, age, body mass index, average daily calorie intake, years of education, exercise levels, smoking status, and menopausal status), and genome wide genotypes. The main outcome measures were average serum urate levels and variance in serum urate levels. ß values (95% confidence intervals) and Bonferroni corrected P values from multivariable linear regression analyses, along with regression partial R2 values, were used to quantitate associations. RESULTS: Seven foods were associated with raised serum urate levels (beer, liquor, wine, potato, poultry, soft drinks, and meat (beef, pork, or lamb)) and eight foods were associated with reduced serum urate levels (eggs, peanuts, cold cereal, skim milk, cheese, brown bread, margarine, and non-citrus fruits) in the male, female, or full cohorts. Three diet scores, constructed on the basis of healthy diet guidelines, were inversely associated with serum urate levels and a fourth, data driven diet pattern positively associated with raised serum urate levels, but each explained ≤0.3% of variance in serum urate. In comparison, 23.9% of variance in serum urate levels was explained by common, genome wide single nucleotide variation. CONCLUSION: In contrast with genetic contributions, diet explains very little variation in serum urate levels in the general population.

An update on the genetics of hyperuricaemia and gout.

A central aspect of the pathogenesis of gout is elevated urate concentrations, which lead to the formation of monosodium urate crystals. The clinical features of gout result from an individual's immune response to these deposited crystals. Genome-wide association studies (GWAS) have confirmed the importance of urate excretion in the control of serum urate levels and the risk of gout and have identified the kidneys, the gut and the liver as sites of urate regulation. The genetic contribution to the progression from hyperuricaemia to gout remains relatively poorly understood, although genes encoding proteins that are involved in the NLRP3 (NOD-, LRR- and pyrin domain-containing 3) inflammasome pathway play a part. Genome-wide and targeted sequencing is beginning to identify uncommon population-specific variants that are associated with urate levels and gout. Mendelian randomization studies using urate-associated genetic variants as unconfounded surrogates for lifelong urate exposure have not supported claims that urate is causal for metabolic conditions that are comorbidities of hyperuricaemia and gout. Genetic studies have also identified genetic variants that predict responsiveness to therapies (for example, urate-lowering drugs) for treatment of hyperuricaemia. Future research should focus on large GWAS (that include asymptomatic hyperuricaemic individuals) and on increasing the use of whole-genome sequencing data to identify uncommon genetic variants with increased penetrance that might provide opportunities for clinical translation.

Mediation analysis to understand genetic relationships between habitual coffee intake and gout.

BACKGROUND: Increased coffee intake is associated with reduced serum urate concentrations and lower risk of gout. Specific alleles of the GCKR, ABCG2, MLXIPL, and CYP1A2 genes have been associated with both reduced coffee intake and increased serum urate in separate genome-wide association studies (GWAS). The aim of this study was to determine whether these single nucleotide polymorphisms (SNPs) influence the risk of gout through their effects on coffee consumption. METHODS: This research was conducted using the UK Biobank Resource. Data were available for 130,966 European participants aged 40-69 years. Gout status and coffee intake were tested for association with four urate-associated SNPs: GCKR (rs1260326), ABCG2 (rs2231142), MLXIPL (rs1178977), and CYP1A2 (rs2472297). Multiple regression and path analysis were used to examine whether coffee consumption mediated the effect of the SNPs on gout risk. RESULTS: Coffee consumption was inversely associated with gout (multivariate adjusted odds ratio (95% confidence interval (CI)) for any coffee consumption 0.75 (0.67-0.84, P = 9 × 10-7)). There was also evidence of a dose-effect with multivariate adjusted odds ratio (95% CI) per cup consumed per day of 0.85 (0.82-0.87, P = 9 × 10-32). The urate-increasing GCKR, ABCG2, MLXIPL, and CYP1A2 alleles were associated with reduced daily coffee consumption, with the strongest associations for CYP1A2 (beta -0.30, P = 8 × 10-40), and MLXIPL (beta -0.17, P = 3 × 10-8), and weaker associations for GCKR (beta -0.07, P = 3 × 10-10) and ABCG2 (beta -0.09, P = 2 × 10-9). The urate-increasing GCKR and ABCG2 alleles were associated with gout (multivariate adjusted p < 5 × 10-8 for both), but the urate-increasing MLXIPL and CYP1A2 alleles were not. In mediation analysis, the direct effects of GCKR and ABCG2 accounted for most of the total effect on gout risk, with much smaller indirect effects mediated by coffee consumption. CONCLUSION: Coffee consumption is inversely associated with risk of gout. Although alleles at several SNPs associate with both lower coffee consumption and higher risk of gout, these SNPs largely influence gout risk directly, rather than indirectly through effects on coffee consumption.