Genome-wide association analyses identify 18 new loci associated with serum urate concentrations

An update on the genetic architecture of hyperuricemia and gout

Genome-wide association studies that scan the genome for common genetic variants associated with phenotype have greatly advanced medical knowledge. Hyperuricemia is no exception, with 28 loci identified. However, genetic control of pathways determining gout in the presence of hyperuricemia is still poorly understood. Two important pathways determining hyperuricemia have been confirmed (renal and gut excretion of uric acid with glycolysis now firmly implicated). Major urate loci are SLC2A9 and ABCG2. Recent studies show that SLC2A9 is involved in renal and gut excretion of uric acid and is implicated in antioxidant defense. Although etiological variants at SLC2A9 are yet to be identified, it is clear that considerable genetic complexity exists at the SLC2A9 locus, with multiple statistically independent genetic variants and local epistatic interactions. The positions of implicated genetic variants within or near chromatin regions involved in transcriptional control suggest that this mechanism (rather than structural changes in SLC2A9) is important in regulating the activity of SLC2A9. ABCG2 is involved primarily in extra-renal uric acid under-excretion with the etiological variant influencing expression. At the other 26 loci, probable causal genes can be identified at three (PDZK1, SLC22A11, and INHBB) with strong candidates at a further 10 loci. Confirmation of the causal gene will require a combination of re-sequencing, trans-ancestral mapping, and correlation of genetic association data with expression data. As expected, the urate loci associate with gout, although inconsistent effect sizes for gout require investigation. Finally, there has been no genome-wide association study using clinically ascertained cases to investigate the causes of gout in the presence of hyperuricemia. In such a study, use of asymptomatic hyperurcemic controls would be expected to increase the ability to detect genetic associations with gout.

A polymorphism in the major gene regulating serum uric acid associates with clinic SBP and the white-coat effect in a family-based study

OBJECTIVES

Hyperuricemia associates with hypertension, but it is uncertain whether this relationship is causal in nature. Glucose transporter 9 (GLUT9) gene is a major genetic determinant of plasma uric acid levels in humans. Since polymorphisms are randomly distributed at mating (Mendelian randomization), studies based on GLUT9 polymorphisms may provide unconfounded assessment of the nature of the link between uric acid and hypertension.

METHODS

We tested the association between uric acid, the rs734553 polymorphism of the GLUT9 gene and arterial pressure in a family-based study including 449 individuals in a genetically homogenous population in Southern Italy.

RESULTS

Serum uric acid levels were strongly associated (P < 0.001) with all components of clinic and 24-h ambulatory blood pressures (BPs). However, only clinic SBP and the white-coat effect (the difference in clinic systolic and daytime systolic ambulatory blood pressure monitoring) associations remained significant after adjustment for classical risk factor and the estimated glomerular filtration rate. Serum uric acid was strongly associated with the risk allele (T) of the rs734553 polymorphism (P < 0.001). Furthermore, TT individuals showed higher clinic SBP (129 + SEM 1 mmHg) than GT (125 + 1 mmHg) and GG individuals (122 + 3 mmHg), as well as a higher white-coat effect (P = 0.02), confirming that the association between uric acid and these BP components is unconfounded by environmental risk factors.

CONCLUSION

Results in this family-based study are compatible with the hypothesis that uric acid is a causal risk factor for hypertension. Trials testing uric acid-lowering interventions are needed to definitively establish the causal implication of hyperuricemia in human hypertension. [Corrected]

Genome-wide association study identifies ABCG2 (BCRP) as an allopurinol transporter and a determinant of drug response

The first-line treatment of hyperuricemia, which causes gout, is allopurinol. The allopurinol response is highly variable, with many users failing to achieve target serum uric acid (SUA) levels. No genome-wide association study (GWAS) has examined the genetic factors affecting allopurinol effectiveness. Using 2,027 subjects in Kaiser Permanente's Genetic Epidemiology Research on Adult Health and Aging (GERA) Cohort, we conducted a GWAS of allopurinol-related SUA reduction, first in the largest ethnic group, non-Hispanic white (NHW) subjects, and then in a stratified transethnic meta-analysis. ABCG2, encoding the efflux pump BCRP, was associated with SUA reduction in NHW subjects (P = 2 × 10(-8) ), and a missense allele (rs2231142) was associated with a reduced response (P = 3 × 10(-7) ) in the meta-analysis. Isotopic uptake studies in cells demonstrated that BCRP transports allopurinol and genetic variants in ABCG2 affect this transport. Collectively, this first GWAS of allopurinol response demonstrates that ABCG2 is a key determinant of response to the drug.

Fucosyltransferase 2 (FUT2) non-secretor status and blood group B are associated with elevated serum lipase activity in asymptomatic subjects, and an increased risk for chronic pancreatitis: a genetic association study

OBJECTIVE

Serum lipase activities above the threefold upper reference limit indicate acute pancreatitis. We investigated whether high lipase activity-within the reference range and in the absence of pancreatitis-are associated with genetic single nucleotide polymorphisms (SNP), and whether these identified SNPs are also associated with clinical pancreatitis.

METHODS

Genome-wide association studies (GWAS) on phenotypes 'serum lipase activity' and 'high serum lipase activity' were conducted including 3966 German volunteers from the population-based Study-of-Health-in-Pomerania (SHIP). Lead SNPs associated on a genome-wide significance level were replicated in two cohorts, 1444 blood donors and 1042 pancreatitis patients.

RESULTS

Initial discovery GWAS detected SNPs within or near genes encoding the ABO blood group specifying transferases A/B (ABO), Fucosyltransferase-2 (FUT2), and Chymotrypsinogen-B2 (CTRB2), to be significantly associated with lipase activity levels in asymptomatic subjects. Replication analyses in blood donors confirmed the association of FUT-2 non-secretor status (OR=1.49; p=0.012) and ABO blood-type-B (OR=2.48; p=7.29×10(-8)) with high lipase activity levels. In pancreatitis patients, significant associations were found for FUT-2 non-secretor status (OR=1.53; p=8.56×10(-4)) and ABO-B (OR=1.69, p=1.0×10(-4)) with chronic pancreatitis, but not with acute pancreatitis. Conversely, carriers of blood group O were less frequently affected by chronic pancreatitis (OR=0.62; p=1.22×10(-05)) and less likely to have high lipase activity levels (OR=0.59; p=8.14×10(-05)).

CONCLUSIONS

These are the first results indicating that ABO blood type-B as well as FUT2 non-secretor status are common population-wide risk factors for developing chronic pancreatitis. They also imply that, even within the reference range, elevated lipase activities may indicate subclinical pancreatic injury in asymptomatic subjects.

Modulation of genetic associations with serum urate levels by body-mass-index in humans

We tested for interactions between body mass index (BMI) and common genetic variants affecting serum urate levels, genome-wide, in up to 42569 participants. Both stratified genome-wide association (GWAS) analyses, in lean, overweight and obese individuals, and regression-type analyses in a non BMI-stratified overall sample were performed. The former did not uncover any novel locus with a major main effect, but supported modulation of effects for some known and potentially new urate loci. The latter highlighted a SNP at RBFOX3 reaching genome-wide significant level (effect size 0.014, 95% CI 0.008-0.02, Pinter= 2.6 x 10-8). Two top loci in interaction term analyses, RBFOX3 and ERO1LB-EDARADD, also displayed suggestive differences in main effect size between the lean and obese strata. All top ranking loci for urate effect differences between BMI categories were novel and most had small magnitude but opposite direction effects between strata. They include the locus RBMS1-TANK (men, Pdifflean-overweight= 4.7 x 10-8), a region that has been associated with several obesity related traits, and TSPYL5 (men, Pdifflean-overweight= 9.1 x 10-8), regulating adipocytes-produced estradiol. The top-ranking known urate loci was ABCG2, the strongest known gout risk locus, with an effect halved in obese compared to lean men (Pdifflean-obese= 2 x 10-4). Finally, pathway analysis suggested a role for N-glycan biosynthesis as a prominent urate-associated pathway in the lean stratum. These results illustrate a potentially powerful way to monitor changes occurring in obesogenic environment.

Uric Acid and Cardiovascular Events: A Mendelian Randomization Study

Obesity and diets rich in uric acid-raising components appear to account for the increased prevalence of hyperuricemia in Westernized populations. Prevalence rates of hypertension, diabetes mellitus, CKD, and cardiovascular disease are also increasing. We used Mendelian randomization to examine whether uric acid is an independent and causal cardiovascular risk factor. Serum uric acid was measured in 3315 patients of the Ludwigshafen Risk and Cardiovascular Health Study. We calculated a weighted genetic risk score (GRS) for uric acid concentration based on eight uric acid-regulating single nucleotide polymorphisms. Causal odds ratios and causal hazard ratios (HRs) were calculated using a two-stage regression estimate with the GRS as the instrumental variable to examine associations with cardiometabolic phenotypes (cross-sectional) and mortality (prospectively) by logistic regression and Cox regression, respectively. Our GRS was not consistently associated with any biochemical marker except for uric acid, arguing against pleiotropy. Uric acid was associated with a range of prevalent diseases, including coronary artery disease. Uric acid and the GRS were both associated with cardiovascular death and sudden cardiac death. In a multivariate model adjusted for factors including medication, causal HRs corresponding to each 1-mg/dl increase in genetically predicted uric acid concentration were significant for cardiovascular death (HR, 1.77; 95% confidence interval, 1.12 to 2.81) and sudden cardiac death (HR, 2.41; 95% confidence interval, 1.16 to 5.00). These results suggest that high uric acid is causally related to adverse cardiovascular outcomes, especially sudden cardiac death.

Renal clearance of uric acid is linked to insulin resistance and lower excretion of sodium in gout patients

Inefficient renal excretion of uric acid is the main pathophysiological mechanism for hyperuricemia in gout patients. Polymorphisms of renal tubular transporters linked with sodium and monosaccharide transport have yet to be demonstrated. We intended to evaluate the impact of insulin resistance, evaluated with the homeostasis model assessment (HOMA), through a transversal study of non-diabetic patients with gout, with normal renal function, not treated with any medication but colchicine as prophylaxis. One hundred and thirty-three patients were evaluated. Clearance of uric acid was inversely correlated with insulin resistance and directly correlated with fractional excretion of sodium. In multivariate analysis, hypertension and hyperlipidemia, in addition to insulin resistance and fractional excretion of sodium, were associated with renal clearance of uric acid. HOMA cutoff for efficient versus inefficient renal handling of uric acid was 2.72, close to that observed in studies of reference population. The impact of insulin resistance and renal handling of sodium on renal clearance of uric acid may help to explain why hyperuricemia is more commonly associated with diabetes and hypertension.

Socioeconomic inequalities show remarkably poor association with health and disease in Southern Croatia

OBJECTIVES

This study aimed at investigating the association of socioeconomic status and health outcomes in populations of the two remote Croatian islands and one coastal city.

METHODS

Medical history and survey information were used to create 33 variables that were analysed using logistic regression. The population from the island of Vis was followed up and mortality data were used to calculate hazard ratios using Cox regression.

RESULTS

Socioeconomic inequalities were poorly associated with health and disease indices. In the matrix of 33 outcome variables and 13 socioeconomic predictor classes, only 10 associations were significant at the level of P < 0.001. None of the associations was replicated across samples. We did not detect the association of any socioeconomic estimate with mortality data for the island of Vis.

CONCLUSIONS

Homogenous island populations were expected to have greater levels of social homogeneity and consequently less expressed inequalities in health. The lack of stronger association in the urban population of Split is likely the result of the mechanisms that persisted from the former communist regime and high level of retained formal and informal social support.

Familial aggregation of gout and relative genetic and environmental contributions: a nationwide population study in Taiwan

OBJECTIVE

To examine familial aggregation of gout and to estimate the heritability and environmental contributions to gout susceptibility in the general population.

METHODS

Using data from the National Health Insurance (NHI) Research Database in Taiwan, we conducted a nationwide cross-sectional study of data collected from 22 643 748 beneficiaries of the NHI in 2004; among them 1 045 059 individuals had physician-diagnosed gout. We estimated relative risks (RR) of gout in individuals with affected first-degree and second-degree relatives and relative contributions of genes (heritability), common environment shared by family members and non-shared environment to gout susceptibility.

RESULTS

RRs for gout were significantly higher in individuals with affected first-degree relatives (men, 1.91 (95% CI 1.90 to 1.93); women, 1.97 (95% CI 1.94 to 1.99)) and also in those with affected second-degree relatives (men, 1.27 (95% CI 1.23 to 1.31); women, 1.40 (95% CI 1.35 to 1.46)). RRs (95% CIs) for individuals with an affected twin, sibling, offspring, parent, grandchild, nephew/niece, uncle/aunt and grandparent were 8.02 (6.95 to 9.26), 2.59 (2.54 to 2.63), 1.96 (1.95 to 1.97), 1.93 (1.91 to 1.94), 1.48 (1.43 to 1.53), 1.40 (1.32 to 1.47), 1.31 (1.24 to 1.39), and 1.26 (1.21 to 1.30), respectively. The relative contributions of heritability, common and non-shared environmental factors to phenotypic variance of gout were 35.1, 28.1 and 36.8% in men and 17.0, 18.5 and 64.5% in women, respectively.

CONCLUSIONS

This population-based study confirms that gout aggregates within families. The risk of gout is higher in people with a family history. Genetic and environmental factors contribute to gout aetiology, and the relative contributions are sexually dimorphic.

Genetic variations in the CLNK gene and ZNF518B gene are associated with gout in case-control sample sets

A genome-wide association study of gout in European populations identified 12 genetic variants strongly associated with risk of gout, but it is unknown whether these variants are also associated with gout risk in Chinese populations. A total of 145 patients with gout and 310 healthy control patients were recruited for a case-control association study. Twelve SNPs of CLNK and ZNF518B gene were genotyped, and association analysis was performed. Odds ratios (ORs) with 95 % confidence intervals (CIs) were used to assess the association. Overall, we found four risk alleles for gout in patients: the allele "G" of rs2041215 and rs1686947 in the CLNK gene by dominant model (OR 1.66; 95 % CI 1.04-2.63; p = 0.031) (OR 2.19; 95 % CI 1.38-3.46; p = 0.001) and additive model (OR 1.39; 95 % CI 1.00-1.93; p = 0.049) (OR 1.67; 95 % CI 1.19-2.32; p = 0.003), respectively, and the allele "A" of rs10938799 and rs10016022 in ZNF518B gene by recessive model (OR 4.66; 95 % CI 1.44-15.09; p = 0.008) (OR 4.54; 95 % CI 1.23-16.76; p = 0.020). Further haplotype analysis showed that the TCATTCTGA haplotype of CLNK was more frequent among patients with gout (adjusted OR 0.48; 95 % CI 0.24-0.95; p = 0.036). Additionally, polymorphisms of rs2041215, rs10938799, and rs17467273 were also correlated with clinical pathological parameters. This study provides evidence for gout susceptibility genes, CLNK and ZNF518B, in a Chinese population, which may have potential as diagnostic and prognostic marker for gout patients.

Insight in genome-wide association of metabolite quantitative traits by exome sequence analyses

Metabolite quantitative traits carry great promise for epidemiological studies, and their genetic background has been addressed using Genome-Wide Association Studies (GWAS). Thus far, the role of less common variants has not been exhaustively studied. Here, we set out a GWAS for metabolite quantitative traits in serum, followed by exome sequence analysis to zoom in on putative causal variants in the associated genes. ¹H Nuclear Magnetic Resonance (¹H-NMR) spectroscopy experiments yielded successful quantification of 42 unique metabolites in 2,482 individuals from The Erasmus Rucphen Family (ERF) study. Heritability of metabolites were estimated by SOLAR. GWAS was performed by linear mixed models, using HapMap imputations. Based on physical vicinity and pathway analyses, candidate genes were screened for coding region variation using exome sequence data. Heritability estimates for metabolites ranged between 10% and 52%. GWAS replicated three known loci in the metabolome wide significance: CPS1 with glycine (P-value  = 1.27×10⁻³²), PRODH with proline (P-value  = 1.11×10⁻¹⁹), SLC16A9 with carnitine level (P-value  = 4.81×10⁻¹⁴) and uncovered a novel association between DMGDH and dimethyl-glycine (P-value  = 1.65×10⁻¹⁹) level. In addition, we found three novel, suggestively significant loci: TNP1 with pyruvate (P-value  = 1.26×10⁻⁸), KCNJ16 with 3-hydroxybutyrate (P-value  = 1.65×10⁻⁸) and 2p12 locus with valine (P-value  = 3.49×10⁻⁸). Exome sequence analysis identified potentially causal coding and regulatory variants located in the genes CPS1, KCNJ2 and PRODH, and revealed allelic heterogeneity for CPS1 and PRODH. Combined GWAS and exome analyses of metabolites detected by high-resolution ¹H-NMR is a robust approach to uncover metabolite quantitative trait loci (mQTL), and the likely causative variants in these loci. It is anticipated that insight in the genetics of intermediate phenotypes will provide additional insight into the genetics of complex traits.

Human metabolic individuality is under strict control of genetic and environmental factors. In our study, we aimed to find the genetic determinants of circulating molecules in sera of large set of individuals representing the general population. First, we performed a hypothesis-free genome wide screen in this population to identify genetic regions of interest. Our study confirmed four known gene metabolite connections, but also pointed to four novel ones. Genome-wide screens enriched for common intergenic variants may miss causal genetic variations directly changing the protein sequence. To investigate this further, we zoomed into regions of interest and tested whether the association signals obtained in the first stage were direct, or whether they represent causal variations, which were not captured in the initial panel. These subsequent tests showed that protein coding and regulatory variations are involved in metabolite levels. For two genomic regions we also found that genes harbour more than one causal variant influencing metabolite levels independent of each other. We also observed strong connection between markers of cardio-metabolic health and metabolites. Taken together, our novel loci are of interest for further research to investigate the causal relation to for instance type 2 diabetes and cardiovascular disease.

Genome-wide association meta-analysis identifies novel variants associated with fasting plasma glucose in East Asians

Fasting plasma glucose (FPG) has been recognized as an important indicator for the overall glycemic state preceding the onset of metabolic diseases. So far, most indentified genome-wide association loci for FPG were derived from populations with European ancestry, with a few exceptions. To extend a thorough catalog for FPG loci, we conducted meta-analyses of 13 genome-wide association studies in up to 24,740 nondiabetic subjects with East Asian ancestry. Follow-up replication analyses in up to an additional 21,345 participants identified three new FPG loci reaching genome-wide significance in or near PDK1-RAPGEF4, KANK1, and IGF1R. Our results could provide additional insight into the genetic variation implicated in fasting glucose regulation.

Association between ABCG2 Q141K polymorphism and gout risk affected by ethnicity and gender: a systematic review and meta-analysis

AIM

Original studies have employed various genetic models in association analysis between ABCG2 Q141K (rs2231142) with gout risk and different or conflicting results, especially regarding the role of gender in this association. In addition, it is not clear whether the association varies by ethnicity.

METHOD

Articles published before September 1, 2013 were extracted and registered into databases for the systematic review of this polymorphism. The quality of each study was scored based on predefined criteria. The genetic model was identified through stratification analysis, then a meta-analysis including all publically available data was preformed to test the association between rs2231142 and gout risk. Potential sources of heterogeneity were sought out via stratification analysis and meta-regression analysis.

RESULTS

Nine case-control studies involving 17 942 individuals were eligible for the meta-analysis of rs2231142. Codominant model was the most appropriate genetic model to interpret the susceptibility cause. It showed that the rs2231142 T allele obviously increased gout risk, and TT was much stronger than GT (TT vs. GG: OR, 4.10; 95% CI, 2.90-5.80; GT vs. GG: OR, 1.71, 95% CI, 1.39-2.10). In addition, gender and ethnicity were found to affect the association between the susceptibility of gout and rs2231142.

CONCLUSION

ABCG2 rs2231142 is an important genetic factor in increasing gout risk, and the difference in genetic association has been found between male and female populations. In addition, the degree of association has been found to vary with ethnicity.

Linking uric acid metabolism to diabetic complications

Hyperuricemia have been thought to be caused by the ingestion of large amounts of purines, and prevention or treatment of hyperuricemia has intended to prevent gout. Xanthine dehydrogenase/xanthine oxidase (XDH/XO) is rate-limiting enzyme of uric acid generation, and allopurinol was developed as a uric acid (UA) generation inhibitor in the 1950s and has been routinely used for gout prevention since then. Serum UA levels are an important risk factor of disease progression for various diseases, including those related to lifestyle. Recently, other UA generation inhibitors such as febuxostat and topiroxostat were launched. The emergence of these novel medications has promoted new research in the field. Lifestyle-related diseases, such as metabolic syndrome or type 2 diabetes mellitus, often have a common pathological foundation. As such, hyperuricemia is often present among these patients. Many in vitro and animal studies have implicated inflammation and oxidative stress in UA metabolism and vascular injury because XDH/XO act as one of the major source of reactive oxygen species Many studies on UA levels and associated diseases implicate involvement of UA generation in disease onset and/or progression. Interventional studies for UA generation, not UA excretion revealed XDH/XO can be the therapeutic target for vascular injury and renal dysfunction. In this review, the relationship between UA metabolism and diabetic complications is highlighted.

Assessing the utility of intermediate phenotypes for genetic mapping of psychiatric disease

Intermediate phenotypes are traits positioned somewhere between genetic variation and disease. They represent a target for attempts to find disease-associated genetic variants and elucidation of mechanisms. Psychiatry has been particularly enamoured with intermediate phenotypes, due to uncertainty about disease aetiology, inconclusive results in early psychiatric genetic studies, and their appeal relative to traditional diagnostic categories. In this review, we argue that new genetic findings are relevant to the question of the utility of these constructs. In particular, results from genome-wide association studies of psychiatric disorders now allow an assessment of the potential role of particular intermediate phenotypes. Based on such an analysis, as well as other recent results, we conclude that intermediate phenotypes are likely to be most valuable in understanding mechanism.

Detection of pleiotropy through a Phenome-wide association study (PheWAS) of epidemiologic data as part of the Environmental Architecture for Genes Linked to Environment (EAGLE) study

We performed a Phenome-wide association study (PheWAS) utilizing diverse genotypic and phenotypic data existing across multiple populations in the National Health and Nutrition Examination Surveys (NHANES), conducted by the Centers for Disease Control and Prevention (CDC), and accessed by the Epidemiological Architecture for Genes Linked to Environment (EAGLE) study. We calculated comprehensive tests of association in Genetic NHANES using 80 SNPs and 1,008 phenotypes (grouped into 184 phenotype classes), stratified by race-ethnicity. Genetic NHANES includes three surveys (NHANES III, 1999-2000, and 2001-2002) and three race-ethnicities: non-Hispanic whites (n = 6,634), non-Hispanic blacks (n = 3,458), and Mexican Americans (n = 3,950). We identified 69 PheWAS associations replicating across surveys for the same SNP, phenotype-class, direction of effect, and race-ethnicity at p<0.01, allele frequency >0.01, and sample size >200. Of these 69 PheWAS associations, 39 replicated previously reported SNP-phenotype associations, 9 were related to previously reported associations, and 21 were novel associations. Fourteen results had the same direction of effect across more than one race-ethnicity: one result was novel, 11 replicated previously reported associations, and two were related to previously reported results. Thirteen SNPs showed evidence of pleiotropy. We further explored results with gene-based biological networks, contrasting the direction of effect for pleiotropic associations across phenotypes. One PheWAS result was ABCG2 missense SNP rs2231142, associated with uric acid levels in both non-Hispanic whites and Mexican Americans, protoporphyrin levels in non-Hispanic whites and Mexican Americans, and blood pressure levels in Mexican Americans. Another example was SNP rs1800588 near LIPC, significantly associated with the novel phenotypes of folate levels (Mexican Americans), vitamin E levels (non-Hispanic whites) and triglyceride levels (non-Hispanic whites), and replication for cholesterol levels. The results of this PheWAS show the utility of this approach for exposing more of the complex genetic architecture underlying multiple traits, through generating novel hypotheses for future research.

Sugar-sweetened beverage consumption: a risk factor for prevalent gout with SLC2A9 genotype-specific effects on serum urate and risk of gout

OBJECTIVE

Consumption of high fructose corn syrup (HFCS)-sweetened beverages increases serum urate and risk of incident gout. Genetic variants in SLC2A9, that exchanges uric acid for glucose and fructose, associate with gout. We tested association between sugar (sucrose)-sweetened beverage (SSB) consumption and prevalent gout. We also tested the hypothesis that SLC2A9 genotype and SSB consumption interact to determine gout risk.

METHODS

Participants were 1634 New Zealand (NZ) European Caucasian, Ma¯ori and Pacific Island people and 7075 European Caucasians from the Atherosclerosis Risk in Communities (ARIC) study. NZ samples were genotyped for rs11942223 and ARIC for rs6449173. Effect estimates were multivariate adjusted.

RESULTS

SSB consumption increased gout risk. The OR for four drinks/day relative to zero was 6.89 (p=0.045), 5.19 (p=0.010) and 2.84 (p=0.043) for European Caucasian, Ma¯ori and Pacific Islanders, respectively. With each extra daily SSB serving, carriage of the gout-protective allele of SLC2A9 associated with a 15% increase in risk (p=0.078), compared with a 12% increase in non-carriers (p=0.002). The interaction term was significant in pooled (pInteraction=0.01) but not meta-analysed (pInteraction=0.99) data. In ARIC, with each extra daily serving, a greater increase in serum urate protective allele carriers (0.005 (p=8.7×10(-5)) compared with 0.002 (p=0.016) mmol/L) supported the gout data (pInteraction=0.062).

CONCLUSIONS

Association of SSB consumption with prevalent gout supports reduction of SSB in management. The interaction data suggest that SLC2A9-mediated renal uric acid excretion is physiologically influenced by intake of simple sugars derived from SSB, with SSB exposure negating the gout risk discrimination of SLC2A9.

The relationship of apolipoprotein B and very low density lipoprotein triglyceride with hyperuricemia and gout

INTRODUCTION

Gout results from an innate immune response to monosodium urate (MSU) crystals deposited in joints. Increased very low-density lipoprotein (VLDL) has been associated with gout. The apolipoprotein B (apo B), which is present on VLDL, regulates neutrophil response to MSU crystals and has been positively associated with gout. Furthermore, the gene (A1CF) encoding the complementation factor for the APOB mRNA-editing enzyme is associated with urate levels. However, the relationship of apo B and VLDL with gout and hyperuricaemia (HU) is still unclear. Therefore, we tested the association of VLDL and apo B with HU and with gout compared to HU.

METHODS

New Zealand European (n = 90) and Māori and Pacific Island (Polynesian) (n = 90) male gout case and control sample sets were divided into normouricaemia (NU), asymptomatic HU and gout groups. Size exclusion chromatography and enzyme-linked immunosorbant assay was used to measure VLDL and apo B. Multivariate logistic regression was used to assess the risk of gout and HU per unit change in VLDL and apo B.

RESULTS

Increased levels of VLDL triglycerides (Tg) were observed in the gout sample set compared to NU and HU in Europeans (P = 1.8 × 10(-6) and 1 × 10(-3), respectively), but only compared to NU in Polynesians (P = 0.023). This increase was driven by increased number of VLDL particles in the European participants and by the Tg-enrichment of existing VLDL particles in the Polynesian participants. Each mmol/L increase in VLDL Tg was significantly associated with gout in the presence of HU in Europeans, with a similar trend in Polynesians (OR = 7.61, P = 0.011 and 2.84, P = 0.069, respectively). Each μmol/L increase in total apo B trended towards decreased risk of HU (OR = 0.47; P = 0.062) and, conversely, with increased risk of gout compared to HU (OR = 5.60; P = 0.004).

CONCLUSIONS

Increased VLDL Tg is associated with the risk of gout compared to HU. A genetic approach should be taken to investigate the possibility for causality of VLDL in gout. Apolipoprotein B may have pleiotropic effects in determining HU and gout.

On the analysis of a repeated measure design in genome-wide association analysis

Longitudinal data enables detecting the effect of aging/time, and as a repeated measures design is statistically more efficient compared to cross-sectional data if the correlations between repeated measurements are not large. In particular, when genotyping cost is more expensive than phenotyping cost, the collection of longitudinal data can be an efficient strategy for genetic association analysis. However, in spite of these advantages, genome-wide association studies (GWAS) with longitudinal data have rarely been analyzed taking this into account. In this report, we calculate the required sample size to achieve 80% power at the genome-wide significance level for both longitudinal and cross-sectional data, and compare their statistical efficiency. Furthermore, we analyzed the GWAS of eight phenotypes with three observations on each individual in the Korean Association Resource (KARE). A linear mixed model allowing for the correlations between observations for each individual was applied to analyze the longitudinal data, and linear regression was used to analyze the first observation on each individual as cross-sectional data. We found 12 novel genome-wide significant disease susceptibility loci that were then confirmed in the Health Examination cohort, as well as some significant interactions between age/sex and SNPs.

Association analysis of the SLC22A11 (organic anion transporter 4) and SLC22A12 (urate transporter 1) urate transporter locus with gout in New Zealand case-control sample sets reveals multiple ancestral-specific effects

INTRODUCTION

There is inconsistent association between urate transporters SLC22A11 (organic anion transporter 4 (OAT4)) and SLC22A12 (urate transporter 1 (URAT1)) and risk of gout. New Zealand (NZ) Māori and Pacific Island people have higher serum urate and more severe gout than European people. The aim of this study was to test genetic variation across the SLC22A11/SLC22A12 locus for association with risk of gout in NZ sample sets.

METHODS

A total of 12 single nucleotide polymorphism (SNP) variants in four haplotype blocks were genotyped using TaqMan® and Sequenom MassArray in 1003 gout cases and 1156 controls. All cases had gout according to the 1977 American Rheumatism Association criteria. Association analysis of single markers and haplotypes was performed using PLINK and Stata.

RESULTS

A haplotype block 1 SNP (rs17299124) (upstream of SLC22A11) was associated with gout in less admixed Polynesian sample sets, but not European Caucasian (odds ratio; OR = 3.38, P = 6.1 × 10-4; OR = 0.91, P = 0.40, respectively) sample sets. A protective block 1 haplotype caused the rs17299124 association (OR = 0.28, P = 6.0 × 10-4). Within haplotype block 2 (SLC22A11) we could not replicate previous reports of association of rs2078267 with gout in European Caucasian (OR = 0.98, P = 0.82) sample sets, however this SNP was associated with gout in Polynesian (OR = 1.51, P = 0.022) sample sets. Within haplotype block 3 (including SLC22A12) analysis of haplotypes revealed a haplotype with trans-ancestral protective effects (OR = 0.80, P = 0.004), and a second haplotype conferring protection in less admixed Polynesian sample sets (OR = 0.63, P = 0.028) but risk in European Caucasian samples (OR = 1.33, P = 0.039).

CONCLUSIONS

Our analysis provides evidence for multiple ancestral-specific effects across the SLC22A11/SLC22A12 locus that presumably influence the activity of OAT4 and URAT1 and risk of gout. Further fine mapping of the association signal is needed using trans-ancestral re-sequence data.

Association of Uric Acid Genetic Risk Score With Blood Pressure

High levels of serum uric acid are associated with hypertension in observational studies. The aim of this study was to investigate the association of uric acid gene variants with blood pressure. We studied 5791 participants aged ≥55 years from the Rotterdam Study. Thirty gene variants identified for serum uric acid level were used to compile genetic risk score (GRS). We used linear regression models to investigate the association of the uric acid GRS with systolic and diastolic blood pressure in the whole study population and separately in participants with and without comorbidities and medication use. In the age- and sex-adjusted model, each SD increase in uric acid GRS was associated with 0.75 mm Hg lower systolic blood pressure (95% confidence interval, −1.31 to −0.19) and 0.42 mm Hg lower diastolic blood pressure (95% confidence interval, −0.72 to −0.13). The association did not attenuate after further adjustment for antihypertensive medication use and conventional cardiovascular risk factors. In subgroup analysis, the association of uric acid GRS with systolic blood pressure was significantly stronger in participants (n=885) on diuretic treatment ( P for interaction, 0.007). In conclusion, we found that higher uric acid GRS is associated with lower systolic and diastolic blood pressure. Diuretics treatment may modify the association of uric acid genetic risk score and systolic blood pressure. Our study suggests that genome wide association study’s findings can be associated with an intermediate factor or have a pleiotropic role and, therefore, should be applied for Mendelian Randomization with caution.

Gout: a review of nonmodifiable and modifiable risk factors

Gout is a common inflammatory arthritis triggered by the crystallization of uric acid within the joints. Gout affects millions worldwide and has an increasing prevalence. Recent research has been carried out to better qualify and quantify the risk factors predisposing individuals to gout. These can largely be broken into nonmodifiable risk factors, such as gender, age, race, and genetics, and modifiable risk factors, such as diet and lifestyle. Increasing knowledge of factors predisposing certain individuals to gout could potentially lead to improved preventive practices. This review summarizes the nonmodifiable and modifiable risk factors associated with development of gout.

Identifying causal variants at loci with multiple signals of association

Although genome-wide association studies have successfully identified thousands of risk loci for complex traits, only a handful of the biologically causal variants, responsible for association at these loci, have been successfully identified. Current statistical methods for identifying causal variants at risk loci either use the strength of the association signal in an iterative conditioning framework or estimate probabilities for variants to be causal. A main drawback of existing methods is that they rely on the simplifying assumption of a single causal variant at each risk locus, which is typically invalid at many risk loci. In this work, we propose a new statistical framework that allows for the possibility of an arbitrary number of causal variants when estimating the posterior probability of a variant being causal. A direct benefit of our approach is that we predict a set of variants for each locus that under reasonable assumptions will contain all of the true causal variants with a high confidence level (e.g., 95%) even when the locus contains multiple causal variants. We use simulations to show that our approach provides 20-50% improvement in our ability to identify the causal variants compared to the existing methods at loci harboring multiple causal variants. We validate our approach using empirical data from an expression QTL study of CHI3L2 to identify new causal variants that affect gene expression at this locus. CAVIAR is publicly available online at http://genetics.cs.ucla.edu/caviar/.

Mendelian randomization provides no evidence for a causal role of serum urate in increasing serum triglyceride levels

BACKGROUND

Triglycerides and their lipoprotein transport molecules are risk factors for heart disease. Observational studies have associated elevated levels of serum urate (SU) with triglycerides and risk of heart disease. However, owing to unmeasured confounding, observational studies do not provide insight into the causal relationship between SU and triglyceride. The aim of this study was to test for a causal role of SU in increasing triglyceride using Mendelian randomization that accounts for unmeasured confounding.

METHODS AND RESULTS

Subjects were of European ancestry from the atherosclerosis risk in communities (n=5237) and Framingham heart (n=2971) studies. Mendelian randomization by the 2-stage least squares regression method was done with SU as the exposure, a uric acid transporter genetic risk score as instrumental variable, and triglyceride as the outcome. In ordinary linear regression, SU was significantly associated with triglyceride levels (β=2.69 mmol/L change in triglyceride per mmol/L increase in SU). However, Mendelian randomization-based estimation showed no evidence for a direct causal association of SU with triglyceride concentration-there was a nonsignificant 1.01 mmol/L decrease in triglyceride per mmol/L increase in SU attributable to the genetic risk score (P=0.21). The reverse analysis using a triglyceride genetic risk score provided evidence of a causal role for triglyceride in raising urate in men (P(Corrected)=0.018).

CONCLUSIONS

These data provide no evidence for a causal role for SU in raising triglyceride levels, consistent with a previous Mendelian randomization report of no association between SU and ischemic heart disease.

Inborn errors of purine metabolism: clinical update and therapies

Inborn errors of purine metabolism exhibit broad neurological, immunological, haematological and renal manifestations. Limited awareness of the phenotypic spectrum, the recent descriptions of newer disorders and considerable genetic heterogeneity, have contributed to long diagnostic odysseys for affected individuals. These enzymes are widely but not ubiquitously distributed in human tissues and are crucial for synthesis of essential nucleotides, such as ATP, which form the basis of DNA and RNA, oxidative phosphorylation, signal transduction and a range of molecular synthetic processes. Depletion of nucleotides or accumulation of toxic intermediates contributes to the pathogenesis of these disorders. Maintenance of cellular nucleotides depends on the three aspects of metabolism of purines (and related pyrimidines): de novo synthesis, catabolism and recycling of these metabolites. At present, treatments for the clinically significant defects of the purine pathway are restricted: purine 5'-nucleotidase deficiency with uridine; familial juvenile hyperuricaemic nephropathy (FJHN), adenine phosphoribosyl transferase (APRT) deficiency, hypoxanthine phosphoribosyl transferase (HPRT) deficiency and phosphoribosyl-pyrophosphate synthetase superactivity (PRPS) with allopurinol; adenosine deaminase (ADA) and purine nucleoside phosphorylase (PNP) deficiencies have been treated by bone marrow transplantation (BMT), and ADA deficiency with enzyme replacement with polyethylene glycol (PEG)-ADA, or erythrocyte-encapsulated ADA; myeloadenylate deaminase (MADA) and adenylosuccinate lyase (ADSL) deficiencies have had trials of oral ribose; PRPS, HPRT and adenosine kinase (ADK) deficiencies with S-adenosylmethionine; and molybdenum cofactor deficiency of complementation group A (MOCODA) with cyclic pyranopterin monophosphate (cPMP). In this review we describe the known inborn errors of purine metabolism, their phenotypic presentations, established diagnostic methodology and recognised treatment options.

Tubular Urate Transporter Gene Polymorphisms Differentiate Patients with Gout Who Have Normal and Decreased Urinary Uric Acid Excretion

Objective.

Primary gout has been associated with single-nucleotide polymorphisms (SNP) in several tubular urate transporter genes. No study has assessed the association of reabsorption and secretion urate transporter gene SNP with gout in a single cohort of documented primary patients with gout carefully subclassified as normoexcretors or underexcretors.

Methods.

Three reabsorption SNP (SLC22A12/URAT1, SLC2A9/GLUT9, and SLC22A11/OAT4) and 2 secretion transporter SNP (SLC17A1/NPT1 and ABCG2/BRCP) were studied in 104 patients with primary gout and in 300 control subjects. The patients were subclassified into normoexcretors and underexcretors according to their serum and 24-h urinary uric acid levels under strict conditions of dietary control.

Results.

Compared with control subjects, patients with gout showed different allele distributions of the 5 SNP analyzed. However, the diagnosis of underexcretor was only positively associated with the presence of the T allele of URAT1 rs11231825, the G allele of GLUT9 rs16890979, and the A allele of ABCG2 rs2231142. The association of the A allele of ABCG2 rs2231142 in normoexcretors was 10 times higher than in underexcretors. The C allele of NPT1 rs1165196 was only significantly associated with gout in patients with normal uric acid excretion.

Conclusion.

Gout with uric acid underexcretion is associated with transporter gene SNP related mainly to tubular reabsorption, whereas uric acid normoexcretion is associated only with tubular secretion SNP. This finding supports the concept of distinctive mechanisms to account for hyperuricemia in patients with gout with reduced or normal uric acid excretion.

Pleiotropic genes for metabolic syndrome and inflammation

Metabolic syndrome (MetS) has become a health and financial burden worldwide. The MetS definition captures clustering of risk factors that predict higher risk for diabetes mellitus and cardiovascular disease. Our study hypothesis is that additional to genes influencing individual MetS risk factors, genetic variants exist that influence MetS and inflammatory markers forming a predisposing MetS genetic network. To test this hypothesis a staged approach was undertaken. (a) We analyzed 17 metabolic and inflammatory traits in more than 85,500 participants from 14 large epidemiological studies within the Cross Consortia Pleiotropy Group. Individuals classified with MetS (NCEP definition), versus those without, showed on average significantly different levels for most inflammatory markers studied. (b) Paired average correlations between 8 metabolic traits and 9 inflammatory markers from the same studies as above, estimated with two methods, and factor analyses on large simulated data, helped in identifying 8 combinations of traits for follow-up in meta-analyses, out of 130,305 possible combinations between metabolic traits and inflammatory markers studied. (c) We performed correlated meta-analyses for 8 metabolic traits and 6 inflammatory markers by using existing GWAS published genetic summary results, with about 2.5 million SNPs from twelve predominantly largest GWAS consortia. These analyses yielded 130 unique SNPs/genes with pleiotropic associations (a SNP/gene associating at least one metabolic trait and one inflammatory marker). Of them twenty-five variants (seven loci newly reported) are proposed as MetS candidates. They map to genes MACF1, KIAA0754, GCKR, GRB14, COBLL1, LOC646736-IRS1, SLC39A8, NELFE, SKIV2L, STK19, TFAP2B, BAZ1B, BCL7B, TBL2, MLXIPL, LPL, TRIB1, ATXN2, HECTD4, PTPN11, ZNF664, PDXDC1, FTO, MC4R and TOMM40. Based on large data evidence, we conclude that inflammation is a feature of MetS and several gene variants show pleiotropic genetic associations across phenotypes and might explain a part of MetS correlated genetic architecture. These findings warrant further functional investigation.

Copy number polymorphisms near SLC2A9 are associated with serum uric acid concentrations

Background

Hyperuricemia is associated with multiple diseases, including gout, cardiovascular disease, and renal disease. Serum urate is highly heritable, yet association studies of single nucleotide polymorphisms (SNPs) and serum uric acid explain a small fraction of the heritability. Whether copy number polymorphisms (CNPs) contribute to uric acid levels is unknown.

Results

We assessed copy number on a genome-wide scale among 8,411 individuals of European ancestry (EA) who participated in the Atherosclerosis Risk in Communities (ARIC) study. CNPs upstream of the urate transporter SLC2A9 on chromosome 4p16.1 are associated with uric acid (χ2df2=3545, p=3.19×10^-23). Effect sizes, expressed as the percentage change in uric acid per deleted copy, are most pronounced among women (_3.974.93_5.87 [ _2.550_97.5 denoting percentiles], p=4.57×10^-23) and independent of previously reported SNPs in SLC2A9 as assessed by SNP and CNP regression models and the phasing SNP and CNP haplotypes (χ2df2=3190,p=7.23×10-08). Our finding is replicated in the Framingham Heart Study (FHS), where the effect size estimated from 4,089 women is comparable to ARIC in direction and magnitude (_1.414.70_7.88, p=5.46×10^-03).

Conclusions

This is the first study to characterize CNPs in ARIC and the first genome-wide analysis of CNPs and uric acid. Our findings suggests a novel, non-coding regulatory mechanism for SLC2A9-mediated modulation of serum uric acid, and detail a bioinformatic approach for assessing the contribution of CNPs to heritable traits in large population-based studies where technical sources of variation are substantial.

12q24 locus association with type 1 diabetes: SH2B3 or ATXN2?

Genetic linkage analyses, genome-wide association studies of single nucleotide polymorphisms, copy number variation surveys, and mutation screenings found the human chromosomal 12q24 locus, with the genes SH2B3 and ATXN2 in its core, to be associated with an exceptionally wide spectrum of disease susceptibilities. Hematopoietic traits of red and white blood cells (like erythrocytosis and myeloproliferative disease), autoimmune disorders (like type 1 diabetes, coeliac disease, juvenile idiopathic arthritis, rheumatoid arthritis, thrombotic antiphospholipid syndrome, lupus erythematosus, multiple sclerosis, hypothyroidism and vitiligo), also vascular pathology (like kidney glomerular filtration rate deficits, serum urate levels, plasma beta-2-microglobulin levels, retinal microcirculation problems, diastolic and systolic blood pressure and hypertension, cardiovascular infarction), furthermore obesity, neurodegenerative conditions (like the polyglutamine-expansion disorder spinocerebellar ataxia type 2, Parkinson’s disease, the motor-neuron disease amyotrophic lateral sclerosis, and progressive supranuclear palsy), and finally longevity were reported. Now it is important to clarify, in which ways the loss or gain of function of the locally encoded proteins SH2B3/LNK and ataxin-2, respectively, contribute to these polygenic health problems. SH2B3/LNK is known to repress the JAK2/ABL1 dependent proliferation of white blood cells. Its null mutations in human and mouse are triggers of autoimmune traits and leukemia (acute lymphoblastic leukemia or chronic myeloid leukemia-like), while missense mutations were found in erythrocytosis-1 patients. Ataxin-2 is known to act on RNA-processing and trophic receptor internalization. While its polyglutamine-expansion mediated gain-of-function causes neuronal atrophy in human and mouse, its deletion leads to obesity and insulin resistance in mice. Thus, it is conceivable that the polygenic pathogenesis of type 1 diabetes is enhanced by an SH2B3-dysregulation-mediated predisposition to autoimmune diseases that conspires with an ATXN2-deficiency-mediated predisposition to lipid and glucose metabolism pathology.

A genome-wide association study identifies PLCL2 and AP3D1-DOT1L-SF3A2 as new susceptibility loci for myocardial infarction in Japanese

Despite considerable progress in preventive and therapeutic strategies, myocardial infarction (MI) is one of the leading causes of death throughout the world. A total of 55 susceptibility genes have been identified mostly in European genome-wide association studies (GWAS). Nevertheless, large-scale GWAS from other population could possibly find additional susceptibility loci. To identify as many MI susceptibility loci as possible, we performed a large-scale genomic analysis in Japanese population. To identify MI susceptibility loci in Japanese, we conducted a GWAS using 1666 cases and 3198 controls using the Illumina Human610-Quad BeadChip and HumanHap550v3 Genotyping BeadChip. We performed replication studies using a total of 11,412 cases and 28,397 controls in the Japanese population. Our study identified two novel susceptibility loci for MI: PLCL2 on chromosome 3p24.3 (rs4618210:A>G, P = 2.60 × 10(-9), odds ratio (OR) = 0.91) and AP3D1-DOT1L-SF3A2 on chromosome 19p13.3 (rs3803915:A>C, P = 3.84 × 10(-9), OR = 0.89). Besides, a total of 14 previously reported MI susceptibility loci were replicated in our study. In particular, we validated a strong association on chromosome 12q24 (rs3782886:A>G: P = 1.14 × 10(-14), OR = 1.46). Following pathway analysis using 265 genes related to MI or coronary artery disease, we found that these loci might be involved in the pathogenesis of MI via the promotion of atherosclerosis. In the present large-scale genomic analysis, we identified PLCL2 and AP3D1-DOT1L-SF3A2 as new susceptibility loci for MI in the Japanese population. Our findings will add novel findings for MI susceptibility loci.

An old disease with new insights: Update on diagnosis and treatment of gout

Gout is an acute and chronic inflammatory disorder associated with high morbidity and impaired quality of life. There has been a substantial increase in the prevalence and incidence of gout in recent years. Novel diagnostic and therapeutic options have provided new insights into the pathogenesis and management of hyperuricemia and gout in the last decade. This clinical review aims to summarize the diagnostic process and management of acute and chronic gout.

Hyperuricaemia in the Pacific: why the elevated serum urate levels?

Pacific Island populations, particularly those of Polynesian descent, have a high prevalence of hyperuricaemia and gout. This is due to an inherently higher urate level among these populations with a demonstrated genetic predisposition. While an excess of urate can cause pathology, urate is also important for human health. It has been implicated as an antioxidant, has a neuroprotective role and is involved in innate immune responses. This paper provides a brief review of urate levels worldwide, with a particular focus on island Southeast Asia and the Pacific. We then present possible evolutionary explanations for the elevated serum urate levels among Pacific populations in the context of the physiological importance of urate and of the settlement history of the region. Finally, we propose that ancestry may play a significant role in hyperuricaemia in these populations and that exposure to malaria prior to population expansion into the wider Pacific may have driven genetic selection for variants contributing to high serum urate.

Gout: joints and beyond, epidemiology, clinical features, treatment and co-morbidities

Gout is a common inflammatory arthritis precipitated by an inflammatory reaction to urate crystals in the joint. Gout is increasingly being recognised as a disease primarily of urate overload with arthritis being a consequence of this pathological accumulation. It is associated with a number of important co-morbidities including chronic kidney disease, obesity, diabetes and cardiovascular disease. The prevalence of gout is increasing around the world. Significant progress has been made in determining the genetic basis for both gout and hyperuricaemia. Environmental risk factors for gout have been identified as certain foods, alcohol and several medications. There is, however, little evidence that changing these environmental risks improves gout on an individual level. Treatment of gout encompasses two strategies: firstly treatment of inflammatory arthritis with non-steroidal anti-inflammatories, corticosteroids, colchicine or interleukin-1 inhibitors. The second and most important strategy is urate lowering, to a target of 0.36 mmol/L (6 mg/dL) or potentially lower in those with tophi (collections of crystalline urate subcutaneously). Along with urate lowering, adequate and prolonged gout flare prophylaxis is required to prevent the precipitation of acute attacks. Newer urate lowering agents are in development and have the potential to significantly expand the potential treatment options. Education of patients regarding the importance of life long urate lowering therapy and prophylaxis of acute attacks is critical to treatment success as adherence with medication is low in chronic diseases in general but especially in gout.

Abundant local interactions in the 4p16.1 region suggest functional mechanisms underlying SLC2A9 associations with human serum uric acid

Human serum uric acid concentration (SUA) is a complex trait. A recent meta-analysis of multiple genome-wide association studies (GWAS) identified 28 loci associated with SUA jointly explaining only 7.7% of the SUA variance, with 3.4% explained by two major loci (SLC2A9 and ABCG2). Here we examined whether gene-gene interactions had any roles in regulating SUA using two large GWAS cohorts included in the meta-analysis [the Atherosclerosis Risk in Communities study cohort (ARIC) and the Framingham Heart Study cohort (FHS)]. We found abundant genome-wide significant local interactions in ARIC in the 4p16.1 region located mostly in an intergenic area near SLC2A9 that were not driven by linkage disequilibrium and were replicated in FHS. Taking the forward selection approach, we constructed a model of five SNPs with marginal effects and three epistatic SNP pairs in ARIC-three marginal SNPs were located within SLC2A9 and the remaining SNPs were all located in the nearby intergenic area. The full model explained 1.5% more SUA variance than that explained by the lead SNP alone, but only 0.3% was contributed by the marginal and epistatic effects of the SNPs in the intergenic area. Functional analysis revealed strong evidence that the epistatically interacting SNPs in the intergenic area were unusually enriched at enhancers active in ENCODE hepatic (HepG2, P = 4.7E-05) and precursor red blood (K562, P = 5.0E-06) cells, putatively regulating transcription of WDR1 and SLC2A9. These results suggest that exploring epistatic interactions is valuable in uncovering the complex functional mechanisms underlying the 4p16.1 region.

Impaired response or insufficient dosage? Examining the potential causes of "inadequate response" to allopurinol in the treatment of gout

OBJECTIVES

Gout is one of the most common forms of arthritis. It is well established that urate-lowering therapy that aims for a serum urate less than at least 0.36 mmol/l (6 mg/dl) is required for the successful management of gout. Allopurinol, a xanthine oxidase (XO) inhibitor, is the most commonly used urate-lowering therapy. However, many patients fail to achieve the target serum urate on allopurinol; these patients can be considered to have "inadequate response" to allopurinol. Herein, we examine the potential mechanisms and implications of inadequate response to allopurinol.

METHODS

The literature was reviewed for potential causes for failure to reach target serum urate in patients receiving allopurinol.

RESULTS

The two most common causes of inadequate response to allopurinol are poor adherence and under-dosing of allopurinol. Adherent patients who fail to achieve target serum urate on standard doses of allopurinol form a group that could be considered to be "partially resistant" to allopurinol. There are four potential mechanisms for partial allopurinol resistance: decreased conversion of allopurinol to oxypurinol; increased renal excretion of oxypurinol; abnormality in XO structure and/or function such that oxypurinol is rendered less effective and/or drug interactions.

CONCLUSIONS

It is important to determine the reasons for failure to achieve treatment targets with allopurinol, particularly as newer agents become available. The knowledge of the mechanisms for inadequate response may help guide the clinician towards making a therapeutic choice that is more likely to result in achieving the serum urate target.

Quality control and conduct of genome-wide association meta-analyses

Rigorous organization and quality control (QC) are necessary to facilitate successful genome-wide association meta-analyses (GWAMAs) of statistics aggregated across multiple genome-wide association studies. This protocol provides guidelines for (i) organizational aspects of GWAMAs, and for (ii) QC at the study file level, the meta-level across studies and the meta-analysis output level. Real-world examples highlight issues experienced and solutions developed by the GIANT Consortium that has conducted meta-analyses including data from 125 studies comprising more than 330,000 individuals. We provide a general protocol for conducting GWAMAs and carrying out QC to minimize errors and to guarantee maximum use of the data. We also include details for the use of a powerful and flexible software package called EasyQC. Precise timings will be greatly influenced by consortium size. For consortia of comparable size to the GIANT Consortium, this protocol takes a minimum of about 10 months to complete.

Fried food consumption, genetic risk, and body mass index: gene-diet interaction analysis in three US cohort studies

OBJECTIVE

To examine the interactions between genetic predisposition and consumption of fried food in relation to body mass index (BMI) and obesity.

DESIGN

Prospective cohort study.

SETTING

Health professionals in the United States.

PARTICIPANTS

9623 women from the Nurses' Health Study, 6379 men from the Health Professionals Follow-up Study, and a replication cohort of 21,421 women from the Women's Genome Health Study.

MAIN OUTCOME MEASURE

Repeated measurement of BMI over follow-up.

RESULTS

There was an interaction between fried food consumption and a genetic risk score based on 32 BMI-associated variants on BMI in both the Nurses' Health Study and Health Professionals Follow-up Study (P ≤ 0.001 for interaction). Among participants in the highest third of the genetic risk score, the differences in BMI between individuals who consumed fried foods four or more times a week and those who consumed fried foods less than once a week amounted to 1.0 (SE 0.2) in women and 0.7 (SE 0.2) in men, whereas the corresponding differences were 0.5 (SE 0.2) and 0.4 (SE 0.2) in the lowest third of the genetic risk score. The gene-diet interaction was replicated in the Women's Genome Health Study (P<0.001 for interaction). Viewed differently, the genetic association with adiposity was strengthened with higher consumption of fried foods. In the combined three cohorts, the differences in BMI per 10 risk alleles were 1.1 (SE 0.2), 1.6 (SE 0.3), and 2.2 (SE 0.6) for fried food consumption less than once, one to three times, and four or more times a week (P<0.001 for interaction); and the odds ratios (95% confidence intervals) for obesity per 10 risk alleles were 1.61 (1.40 to 1.87), 2.12 (1.73 to 2.59), and 2.72 (2.12 to 3.48) across the three categories of consumption (P=0.002 for interaction). In addition, the variants in or near genes highly expressed or known to act in the central nervous system showed significant interactions with fried food consumption, with the FTO (fat mass and obesity associated) variant showing the strongest result (P<0.001 for interaction).

CONCLUSION

Our findings suggest that consumption of fried food could interact with genetic background in relation to obesity, highlighting the particular importance of reducing fried food consumption in individuals genetically predisposed to obesity.

Genetic determinants influencing human serum metabolome among African Americans

Phenotypes proximal to gene action generally reflect larger genetic effect sizes than those that are distant. The human metabolome, a result of multiple cellular and biological processes, are functional intermediate phenotypes proximal to gene action. Here, we present a genome-wide association study of 308 untargeted metabolite levels among African Americans from the Atherosclerosis Risk in Communities (ARIC) Study. Nineteen significant common variant-metabolite associations were identified, including 13 novel loci (p<1.6×10⁻¹⁰). These loci were associated with 7–50% of the difference in metabolite levels per allele, and the variance explained ranged from 4% to 20%. Fourteen genes were identified within the nineteen loci, and four of them contained non-synonymous substitutions in four enzyme-encoding genes (KLKB1, SIAE, CPS1, and NAT8); the other significant loci consist of eight other enzyme-encoding genes (ACE, GATM, ACY3, ACSM2B, THEM4, ADH4, UGT1A, TREH), a transporter gene (SLC6A13) and a polycystin protein gene (PKD2L1). In addition, four potential disease-associated paths were identified, including two direct longitudinal predictive relationships: NAT8 with N-acetylornithine, N-acetyl-1-methylhistidine and incident chronic kidney disease, and TREH with trehalose and incident diabetes. These results highlight the value of using endophenotypes proximal to gene function to discover new insights into biology and disease pathology.

Most contemporary GWAS studies have achieved increased power by increasing the size of the discovery sample to tens of thousands of individuals. An alternative approach for detecting the effects of novel loci is to measure phenotypes that more immediately reflect the effects of gene function. The metabolome consists of a collection of small molecules resulting from a variety of cellular and biologic processes, which can be considered intermediate phenotypes proximal to gene function. Here, we report a genome-wide association study identifying nineteen genetic loci influencing untargeted metabolomes traits among African Americans in the Atherosclerosis Risk in Communities (ARIC) Study. Fourteen genes mapped within nineteen loci, including twelve enzyme-encoding genes (KLKB1, SIAE, CPS1, NAT8, ACE, GATM, ACY3, ACSM2B, THEM4, ADH4, UGT1A and TREH), a transporter gene (SLC6A13) and a polycystin protein gene (PKD2L1). In addition, four potential disease-associated paths were identified, including two direct longitudinal predictive relationships: NAT8 with N-acetylornithine, N-acetyl-1-methylhistidine and incident chronic kidney disease, and TREH with trehalose and incident diabetes. These results highlight the value of using phenotypes proximal to gene function to promote novel gene discovery.

Gene-centric meta-analysis in 87,736 individuals of European ancestry identifies multiple blood-pressure-related loci

Blood pressure (BP) is a heritable risk factor for cardiovascular disease. To investigate genetic associations with systolic BP (SBP), diastolic BP (DBP), mean arterial pressure (MAP), and pulse pressure (PP), we genotyped ~50,000 SNPs in up to 87,736 individuals of European ancestry and combined these in a meta-analysis. We replicated findings in an independent set of 68,368 individuals of European ancestry. Our analyses identified 11 previously undescribed associations in independent loci containing 31 genes including PDE1A, HLA-DQB1, CDK6, PRKAG2, VCL, H19, NUCB2, RELA, HOXC@ complex, FBN1, and NFAT5 at the Bonferroni-corrected array-wide significance threshold (p < 6 × 10(-7)) and confirmed 27 previously reported associations. Bioinformatic analysis of the 11 loci provided support for a putative role in hypertension of several genes, such as CDK6 and NUCB2. Analysis of potential pharmacological targets in databases of small molecules showed that ten of the genes are predicted to be a target for small molecules. In summary, we identified previously unknown loci associated with BP. Our findings extend our understanding of genes involved in BP regulation, which may provide new targets for therapeutic intervention or drug response stratification.

Diastrophic dysplasia sulfate transporter (SLC26A2) is expressed in the adrenal cortex and regulates aldosterone secretion

Elucidation of the molecular mechanisms leading to autonomous aldosterone secretion is a prerequisite to define potential targets and biomarkers in the context of primary aldosteronism. After a genome-wide association study with subjects from the population-based Cooperative Health Research in the Region of Augsburg F4 survey, we observed a highly significant association (P=6.78×10(-11)) between the aldosterone to renin ratio and a locus at 5q32. Hypothesizing that this locus may contain genes of relevance for the pathogenesis of primary aldosteronism, we investigated solute carrier family 26 member 2 (SLC26A2), a protein with known transport activity for sulfate and other cations. Within murine tissues, adrenal glands showed the highest expression levels for SLC26A2, which was significantly downregulated on in vivo stimulation with angiotensin II and potassium. SLC26A2 expression was found to be significantly lower in aldosterone-producing adenomas in comparison with normal adrenal glands. In adrenocortical NCI-H295R cells, specific knockdown of SLC26A2 resulted in a highly significant increase in aldosterone secretion. Concomitantly, expression of steroidogenic enzymes, as well as upstream effectors including transcription factors such as NR4A1, CAMK1, and intracellular Ca(2+) content, was upregulated in knockdown cells. To substantiate further these findings in an SLC26A2 mutant mouse model, aldosterone output proved to be increased in a sex-specific manner. In summary, these findings point toward a possible effect of SLC26A2 in the regulation of aldosterone secretion potentially involved in the pathogenesis of primary aldosteronism.

Identification of novel genetic Loci associated with thyroid peroxidase antibodies and clinical thyroid disease

Autoimmune thyroid diseases (AITD) are common, affecting 2-5% of the general population. Individuals with positive thyroid peroxidase antibodies (TPOAbs) have an increased risk of autoimmune hypothyroidism (Hashimoto's thyroiditis), as well as autoimmune hyperthyroidism (Graves' disease). As the possible causative genes of TPOAbs and AITD remain largely unknown, we performed GWAS meta-analyses in 18,297 individuals for TPOAb-positivity (1769 TPOAb-positives and 16,528 TPOAb-negatives) and in 12,353 individuals for TPOAb serum levels, with replication in 8,990 individuals. Significant associations (P<5×10(-8)) were detected at TPO-rs11675434, ATXN2-rs653178, and BACH2-rs10944479 for TPOAb-positivity, and at TPO-rs11675434, MAGI3-rs1230666, and KALRN-rs2010099 for TPOAb levels. Individual and combined effects (genetic risk scores) of these variants on (subclinical) hypo- and hyperthyroidism, goiter and thyroid cancer were studied. Individuals with a high genetic risk score had, besides an increased risk of TPOAb-positivity (OR: 2.18, 95% CI 1.68-2.81, P = 8.1×10(-8)), a higher risk of increased thyroid-stimulating hormone levels (OR: 1.51, 95% CI 1.26-1.82, P = 2.9×10(-6)), as well as a decreased risk of goiter (OR: 0.77, 95% CI 0.66-0.89, P = 6.5×10(-4)). The MAGI3 and BACH2 variants were associated with an increased risk of hyperthyroidism, which was replicated in an independent cohort of patients with Graves' disease (OR: 1.37, 95% CI 1.22-1.54, P = 1.2×10(-7) and OR: 1.25, 95% CI 1.12-1.39, P = 6.2×10(-5)). The MAGI3 variant was also associated with an increased risk of hypothyroidism (OR: 1.57, 95% CI 1.18-2.10, P = 1.9×10(-3)). This first GWAS meta-analysis for TPOAbs identified five newly associated loci, three of which were also associated with clinical thyroid disease. With these markers we identified a large subgroup in the general population with a substantially increased risk of TPOAbs. The results provide insight into why individuals with thyroid autoimmunity do or do not eventually develop thyroid disease, and these markers may therefore predict which TPOAb-positives are particularly at risk of developing clinical thyroid dysfunction.