SLC2A9 Genotype Is Associated with SLC2A9 Gene Expression and Urinary Uric Acid Concentration

Erin B. Ware; Ellen Riehle; Jennifer A. Smith; Wei Zhao; Stephen T. Turner; Sharon L. R. Kardia; John C. Lieske

doi:10.1371/journal.pone.0128593

Abstract

Objectives

SLC2A9 gene variants have been associated with urinary uric acid (UA) concentration, but little is known about the functional mechanism linking these gene variants with UA. SLC2A9 encodes a UA transporter present in the proximal tubule of the kidney, and gene expression levels of SLC2A9 and other genes in the uricosuric pathway (ABCG2, SLC17A1, SLC17A3, and SLC22A12) could potentially mediate the relationship between SLC2A9 gene variants and urinary UA excretion.

Methods

The association between urinary UA concentrations and single nucleotide polymorphisms (SNPs) within the SLC2A9 gene region, expression levels of genes in the uricosuric pathway, and dietary protein intake were analyzed for a sample of non-Hispanic white participants from the Genetic Epidemiology Network of Arteriopathy (GENOA) cohort. The SLC2A9 SNP most significantly associated with urinary UA concentration was then tested for associations with gene expression levels from uric acid absorption/secretion associated genes. Models including interactions between dietary protein (total, animal, and vegetable) and genetic factors were also assessed.

Results

The most significant SLC2A9 SNP associated with urinary UA (rs12509955, corrected p = 0.001) was also associated with SLC2A9 gene expression levels (corrected p = 0.0084); however, SLC2A9 gene expression levels were not significantly associated with urinary UA concentrations (p = 0.509). The interactions between rs12509955 and total dietary protein, and SLC2A9 gene-level gene expression and dietary vegetable protein on the outcome of urinary UA were marginally significant (p = 0.11 and p = 0.07, respectively). Gene expression level of one SLC2A9 transcript had a significant interaction with dietary animal protein (SLC2A9-001 ENST00000506583, p = 0.01) and a marginally significant interaction with total dietary protein (p = 0.07) on urinary UA.

Conclusion

Our results illustrate that SNPs in the SLC2A9 gene influence SLC2A9 gene expression as well as urinary UA excretion. Evidence is also suggestive that gene-by-diet interactions may disproportionately increase urinary UA in genetically susceptible individuals that consume higher amounts of protein.

Citation: Ware EB, Riehle E, Smith JA, Zhao W, Turner ST, Kardia SLR, et al. (2015) SLC2A9 Genotype Is Associated with SLC2A9 Gene Expression and Urinary Uric Acid Concentration. PLoS ONE 10(7): e0128593. https://doi.org/10.1371/journal.pone.0128593

Editor: Dana C. Crawford, Case Western Reserve University, UNITED STATES

Received: July 29, 2014; Accepted: April 28, 2015; Published: July 13, 2015

Copyright: © 2015 Ware et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited

Data Availability: Due to the sensitive nature of the individual genetic data, data are from the GENOA study whose authors may be contacted at skardia@umich.edu. The reason for this is that the consent form for the GENOA study did not specifically allow for public posting of gene expression or genotype data to publically available sites.

Funding: This work was supported by R01 DK077950, R01 DK073537, U01 HL054457, R01 HL087660, R01 HL119443, and the Mayo Clinic O’Brien Urology Research Center U54 DK100227, and Grant Number UL1 TR000135 from the National Center for Advancing Translational Sciences (NCATS), all funded by the National Institutes of Health. Its contents are solely the responsibility of the authors and do not necessarily represent the official views of the NIH. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Competing interests: The authors have declared that no competing interests exist.

Introduction

Nephrolithiasis (NL), or kidney stone formation, represents a substantial public health burden affecting approximately 1 in 11 Americans [1–3]. Specific diet patterns–particularly protein consumption–and familial history are both known risk factors for NL [4–11]. Most stones in the general population (70–80%) contain a majority of calcium oxalate (CaOx), while approximately 5–10% contain uric acid (UA) or a combination of UA and CaOx [3, 12, 13]. Hyperuricosuria is a risk factor for not only UA stone formation, but can also induce salting out of CaOx from solution, promoting CaOx stone formation [14–18]. Thus, factors influencing urinary UA excretion are likely important modifiers of kidney stone risk [11, 19, 20].

The concentration of UA in urine is determined by a combination of renal filtration followed by reabsorption and secretion in the nephron. SLC2A9 encodes a key UA transporter expressed as two isoforms in the apical and basolateral membranes of the renal proximal tubule. This transporter also mediates uptake of UA from the blood into the liver [21–25]. Genetic variation in this gene has been associated with both mild and severe UA phenotypes, as well as NL [26–29]. Previous genome-wide association studies have identified single nucleotide polymorphisms (SNPs) within SLC2A9 that are associated with serum and urinary UA concentration [30–36]. Several other genes have been found to be involved in excretion of UA in the urine as well as reabsorption within the kidney, though to a lesser extent than SLC2A9, and include SLC22A12, SLC17A1, SLC17A3, and ABCG2 [37–40]. In addition, diet also plays a key role in UA homeostasis, and several studies have shown that higher protein (and purine) intake is associated with higher urinary UA concentration [41–44].

To better understand the contributions of genetic and dietary factors in urinary UA excretion, we examined the influence of SLC2A9 gene variation on urinary UA, gene expression from uric acid absorption/secretion associated genes, and potential interactions with dietary protein. Since our gene expression measures are taken from transformed β lymphocytes, they are best thought of as a cumulative measure of the many genetic variations in each gene that are influencing its expression. Specifically, we identified the SNP within the SLC2A9 gene region (rs12509955) that was most significantly associated with urinary UA excretion using a linear mixed model (LMM), and then examined its association with SLC2A9, ABCG2, SLC17A1, SLC17A3, and SLC22A12 gene expression levels. We further investigated gene-by-dietary protein interaction effects on urinary UA to determine whether rs12509955 genotype or gene expression from uric acid absorption/secretion associated genes interacts with dietary protein intake, since purine metabolism is a known modifier of UA phenotypes.

Methods

Study Population

This study included non-Hispanic whites from Rochester, Minnesota that were enrolled in the Genetic Epidemiology Network of Arteriopathy (GENOA) study (Phase I: 1995–2000, n = 1,583; Phase II: 2000–2005, n = 1,241). GENOA is a multicenter, community-based study of hypertensive sibships initially collected to identify genes that influence blood pressure and target organ damage due to hypertension. The Genetic Determinants of Urinary Lithogenicity (GDUL) study (2008–2012) is an ancillary study of the Phase III GENOA Genetics of Chronic Kidney Disease Study (R01 DK073537), undertaken to investigate predictors of urinary supersaturation and risk of kidney stone diseases among participants without end-stage renal failure (Stage 5 Chronic Kidney Disease). GENOA participants were invited to join the GDUL study which consisted of a study visit; one (or preferably two or three) 24-hour urine collections for determination of quantitative urinary lithogenic factors, and a food frequency questionnaire (FFQ) (Viocare Technologies, Princeton, NJ, USA). Age, sex, and body mass index (BMI) were ascertained at the time of the GDUL exam. Of the 811 participants enrolled in GDUL, two were excluded from the initial SNP discovery analysis because their urinary UA concentration measure was >4 standard deviations from the mean. A subset who had gene expression data (n = 541 participants in 318 sibships) were used for all other association studies. A total of 424 of these participants also had animal and total dietary protein information (393 had vegetable protein) from the food frequency questionnaire. While urinary uric acid is the primary outcome of interest, blood serum uric acid was available from Phase II of the GENOA study, conducted approximately five years prior to the GDUL study. Serum uric acid was investigated to evaluate whether the effects of SNP variation and gene expression on urinary uric acid were consistent for serum uric acid. No participants were excluded for missing serum uric acid, as this was a supplemental analysis. All serum uric acid analyses can be found in supplemental material.

Ethics Statement

This study protocol was approved by the Institutional Review Board at the Mayo Clinic under “Genetic Determinants of Urine Lithogenicity–Working Protocol” 08–006238.

Urine Collection

After the initial visit, subjects completed one, two, or three 24-hour urine collections. Urine was collected with toluene as a preservative. Measurements were averaged for subjects with more than one 24-hour urine collection sample to account for the day-to-day variability in urinary composition. Twenty four-hour urine uric acid (mg/day) and sodium (mmol) concentrations were measured in the Mayo Clinic Renal Testing Laboratory. For additional details, see Lieske, et al. 2014 [45].

Food frequency

Participants completed a food frequency questionnaire (Viocare Technologies, Princeton, NJ, USA) at a study visit. The Women's Health Initiative FFQ was adapted by Viocare in an electronic computer-administered format, and is similar to the Willett or Harvard Food Frequency Questionnaire. It assesses food intake over the last 3 months and employs the Minnesota NDS nutritional analysis database for nutrient analyses. The FFQ was validated and compares well with nutrient intake against 4-day food records and 24-hr dietary recalls [46]. The FFQ has validity for variables of interest to this study, when compared to a 4-day diet diary, namely for protein intake (r = 0.41) and sodium intake (r = 0.31). Dietary protein was assessed as total dietary protein, and also subdivided into animal or vegetable protein sources.

Genotyping and Imputation

All GENOA participants were genotyped on the Affymetrix Genome-Wide Human SNP Array 6.0. A small portion of the GENOA stored blood samples contained poor quality DNA, and some of these were successfully genotyped using the Illumina Human 1M-Duo BeadChip. Since the two platforms used for genotyping contain only a small number of overlapping SNPs (~200,000), association analyses were performed using only imputed data. Prior to imputation, SNPs and samples with a call rate less than 95% were excluded. Imputation was performed using a single-step approach implemented in MACH v1.0.16 using the CEU reference panel of HapMap2 (release 22). The SLC2A9 gene region of interest (including all SNPs within the gene, plus 200kb upstream and downstream) included a total of 880 SNPs. For analyses including the rs12509955 SNP, participants were categorized into imputed genotype (i) calls: (i <0.5 = 0; 0.5< i <1.5 = 1; i>1.5 = 2).

Gene Expression Assessment and Quality Control

Blood samples for beta-lymphocyte extraction were collected from a subset of participants during GENOA Phase I and Phase II. RNA samples were extracted using standard protocols. RNA quality was assessed using the Agilent 2100 Bioanalyzer (Agilent Technologies Inc., Foster City, CA) and quantified by spectrophotometry using the Nanodrop ND-1000 (Nanodrop Inc., Wilmington, DE).

Array quality control was performed at the transcript level with the Affymetrix Expression Console (v 1.1) using core-level probe sets. All array images passed visual inspection. Hybridization controls were all detectable by signal increases that followed concentration. Signal intensity plots were examined for raw and processed data to identify outliers. Raw intensity data were processed using the Affymetrix Power Tool software. Probe summarization and set normalization were performed using Robust Multi-array Analysis (RMA), including background correction, quantile normalization, log₂-transformation, and probe set summarization. Probe sets which are known to cross-hybridize and those with undetectable expression were also excluded. Gene-level expression was assessed by averaging all core probe sets for that gene. Transformed beta-lymphocyte gene expression measurements were adjusted for age, sex, and batch prior to analysis.

Statistical Methods

Data management and statistical analyses were conducted in SAS version 9.3 [47]. Urinary phenotypes had approximately normal distributions; thus, no variable transformations were applied.

The sample of 809 GENOA participants with genotype and urinary UA excretion data was used to identify the SNP with the most significant p-value for association between SNP dosage (as an additive effect) and urinary UA. Fixed effects linear mixed modeling (LMM) was used to account for familial correlations within sibships. Age, sex, body mass index (BMI), and urinary sodium were included in the model because they significantly predict urinary UA in bivariate analyses (Table 1) or have been included as a covariate in other genetic models of urinary UA [32]. In order to investigate whether there were multiple independent signals within the SLC2A9 region, we also tested for association between each SNP in the region and urinary UA, conditioned on the SNP mostly significantly associated with UA.

Download:

Table 1. Characteristics of study sample and bivariate associations with Urinary UA excretion and gene expression levels from uric acid absorption/secretion associated genes, n = 541.

https://doi.org/10.1371/journal.pone.0128593.t001

In a subsample of 541 GDUL participants with gene expression data, LMM was used to assess whether there was an association between the rs12509955 genotype and gene expression levels from uric acid absorption/secretion associated genes. To further evaluate whether a genetic effect on gene expression levels had a corresponding influence on urinary UA, we used LMMs to test the association between gene expression and urinary UA excretion, adjusting for BMI and urinary sodium. LocusZoom plots were used to visualize results of the SNP associations with UA and gene expression [48].

To assess gene-by-diet interactions on urinary UA on rs12509955, LMMs were adjusted for age, sex, BMI, and urinary sodium, and included terms for SNP, dietary protein, and the interaction between dietary protein and SNP. Similarly, we assessed the potential interaction between each dietary protein measure and gene expression levels from uric acid absorption/secretion associated genes on urinary UA excretion. Interaction term β coefficients were considered to be statistically significant at an alpha level of 0.05.

Results

Descriptive Statistics

Table 1 provides descriptive statistics regarding demographic, anthropometric, dietary and genetic risk factors for the study sample and their regression relationships to urinary UA excretion and SLC2A9 gene expression. In our sample, as expected, males had significantly higher urinary UA excretion (β = 132.6, p-value < 0.0001) compared to females. BMI and urinary sodium excretion were positively associated with urinary UA excretion (β = 5.5, p-value < 0.0001; β = 1.8, p-value < 0.0001). Total dietary protein and dietary animal protein were positively associated with urinary UA excretion (β = 0.8, p-value = 0.0003; β = 0.7, p-value = 0.004, respectively). Participants who were in the larger SNP discovery analysis sample (n = 809) but lacked gene expression data (n = 268) were not significantly different with regard to urinary UA excretion than the 541 participants with gene expression data (p-value = 0.10). Table 1 also displays the distribution of rs12509955 genotypes in our study sample (copies of the coded, minor allele (T): 0 = 67.8%; 1 = 28.3%; 2 = 3.9%) and their corresponding mean urinary UA excretion levels. The rs12509955 genotype was significantly associated with urinary UA (β = 33.32, p-value = 0.01), serum UA (β = -0.47, p-value <0.001), and SLC2A9 gene expression (β = 0.07, p-value <0.001). Correlations for gene-level gene expression levels for the five investigated genes were moderate and significant (all p-values <0.0001) and ranged from r = -0.61 (SLC2A9-001 and SLC2A9-002) to r = 0.55 (SLC17A1 and SLC17A3). A correlation table is provided in the supplementary material (S1 Table).

SLC2A9 SNP association with urinary UA excretion

Of the 880 SNPs tested for an association with urinary UA (adjusted for sibships, age, sex, BMI, and urinary sodium in the LMM), the most significant SNP was rs12509955 (intron 2, coded allele = T, β = 40.0, corrected p-value = 0.001). P-values were corrected for the number of independent test in the SLC2A9 gene region using principal components to account for correlated SNPs. After running a principal component analysis, 28 principal components were required to explain 95% of the variation in genotype and thus p-values were corrected by a factor of 28 [49]. As seen in Fig 1, the SCL2A9 gene region has a strong linkage disequilibrium pattern that extends across the SLC2A9 gene, as well as the intergenic regions upstream. There are many SNPs that are strongly associated with urinary UA (S2 Table); however, after conditioning upon the most significant SNP, rs12509955, no additional SNPs in the region were significantly associated with urinary UA after multiple testing correction (S3 Table). We performed this analysis to detect any alternate independent predictive gene variants in the region in the context of a strong linkage disequilibrium pattern. Thus, for subsequent analyses, rs12509955 was selected to represent genetic variation in the SLC2A9 region associated with urinary UA excretion. Since a subsample of 541 participants with gene expression data was used for the remainder of the analysis, we repeated the associations between SNPs and UA in this subsample and again found that rs12509955 was the most statistically significant (β = 40.1, corrected p-value = 0.001). Equivalent analyses for serum uric acid can be found in supplemental material (S4 Table and S5 Table, S1 Fig and S2 Fig).

Download:

Fig 1. Urinary UA association results for 880 SNPs in the SLC2A9 gene region.

Left Y-axis:–log₁₀(p-value) from association between SNPs and urinary UA, adjusted for age, sex, BMI, and urinary sodium, and accounting for sibship; Right Y-axis: SNP recombination rate based on HapMap hg18 CEU; X-axis: chromosomal location and gene regions; r² color code: degree of linkage disequilibrium with index (most strongly associated) SNP, rs12509955 (purple diamond).

https://doi.org/10.1371/journal.pone.0128593.g001

SLC2A9 SNP association with gene expression

Next, 880 SNPs within the SLC2A9 region were tested for association with SLC2A9 gene expression. Of the 880 SNPs tested, 429 SNPs within the region were significantly associated with SLC2A9 gene expression (corrected p-value < 0.05; Fig 2), with the most significant being rs2240724 (β = -0.08, corrected p-value = 6.10x10^-16). As displayed in Table 1, a significant association between rs12509955 genotype and standardized SLC2A9 gene-level gene expression was observed (β = 0.16, corrected p-value = 0.008). rs12509955 was not in linkage disequilibrium with rs2240724, and thus had an independent, more modest influence on gene expression than rs2240704 and the many other SNPs in the region that were associated with SLC2A9 expression (Fig 2). rs2240704 was not investigated further because this SNP was not significantly and independently associated with urinary UA excretion, our main outcome of interest, in a model conditioned upon rs12509955. Overall, the tight linkage disequilibrium pattern indicated that gene expression levels are representing variation in over 48% of SNPs.

Download:

Fig 2. SLC2A9 gene expression association results for 880 SNPs in the SLC2A9 gene region.

Left Y-axis:–log₁₀(p-value) from association between SNPs and SLC2A9 gene expression, accounting for sibship; Right Y-axis: SNP recombination rate based on HapMap hg18 CEU; X-axis: chromosomal location and gene regions; r² color code: degree of linkage disequilibrium with index (most strongly associated) SNP, rs2240724 (purple diamond).

https://doi.org/10.1371/journal.pone.0128593.g002

No SNPs in the SLC2A9 gene region were significantly associated with gene expression of the other four genes of interest after multiple testing correction (α = 0.001). The rs12509955 genotype was not associated with gene expression for SLC22A12 (β = 0.04, p-value = 1), SLC17A1 (β = 0.02, p-value = 1), SLC17A3 (β = 0.05, p-value = 0.75), or ABCG2 (β = 0.00, p-value = 1), corrected for multiple testing.

Gene expression association with urinary UA

SLC2A9 gene expression level, representing SLC2A9 gene variation across subjects, was not significantly associated with urinary UA excretion in either a bivariate model (Table 1, β = 18.0, p-value = 0.51) or a model adjusted for BMI and urinary sodium concentration (β = 15.5, p-value = 0.47). In bivariate models, none of the gene expressions were significantly associated with urinary uric acid (ABCG2 (β = -13.4, p-value = 0.68); SLC17A1 (β = -31.9, p-value = 0.34); SLC17A3 (β = 6.59, p-value = 0.77); SLC22A12 (β = -39.2, p-value = 0.16)), nor was any gene’s expression significantly associated with urinary uric acid after adjustment for BMI and urinary sodium (ABCG2 (β = -47.3, p-value = 0.06); SLC17A1 (β = -39.2, p-value = 0.14); SLC17A3 (β = -12.6, p-value = 0.48); SLC22A12 (β = -40.7, p-value = 0.06)), though both ABCG2 and SLC22A12 may be considered to have marginally significant effects. We did not adjust for age or sex in these models since gene expression had previously been adjusted for age, sex, and batch during quality control. Based on these marginally significant effects, indicating possible gene expression associations with urinary UA, we went on to investigate genetic interactions for all genes of interest (including transcript-based gene expression for the SLC2A9 gene) with dietary protein. Equivalent analyses for serum uric acid can be found in supplemental material (S1 Text).

Gene-by-dietary protein interactions

Models based upon urinary UA excretion and rs12509955-by-dietary protein interactions showed no significant interactions between dietary protein source and rs12509955 genotype, though the interaction between total protein and rs12509955 genotype was marginally significant (total protein: β = 0.50, p-value = 0.11; animal protein: β = 0.01, p-value = 0.98; vegetable protein β = 1.00, p-value = 0.27). Though not significant, those with two copies of the T allele had increased predicted urinary UA compared to those with one or those with no copies of the T allele at every level of dietary protein consumption for each type of dietary protein.

Interaction models (Table 2, serum UA S7 Table) examining the associations between gene expression and dietary protein, or the cumulative measure of variation in the gene and dietary protein, on urinary uric acid showed a marginally significant association between SLC2A9 gene expression and dietary vegetable protein intake (β = 3.21, p-value = 0.07) and non-significant interactions with the other two protein sources (total protein: β = 1.00, p-value = 0.18; animal protein: β = 0.63, p-value = 0.52). No other genes showed significant or marginally significant interactions between gene expression and any of the dietary protein measures.

Download:

Table 2. Gene expression by dietary protein intake interaction associations for urinary uric acid.

https://doi.org/10.1371/journal.pone.0128593.t002

According to the literature, there are 13 splice variants in the SLC2A9 gene region, four of which code for proteins [50]. In our data we observed three protein coding transcripts and two of them appeared at sufficient frequency (<5% missing data) to include in linear mixed models (SLC2A9-001 ENST00000506583, bp = 1927, protein 511 aa, n = 515; SLC2A9-002 ENST00000264784, bp = 1850, protein 540 aa, n = 541). These two transcripts are not significantly associated with urinary UA in bivariate models (SLC2A9-001 β = 1.72, p-value = 0.82; SLC2A9-002 β = -3.85, p-value = 0.59) or in models adjusted for BMI and urinary sodium (SLC2A9-001 β = 3.95, p-value = 0.50; SLC2A9-002 β = -2.24, p-value = 0.69). However, the SLC2A9-001 gene expression transcript significantly interacts with dietary animal protein (β = 0.65, p-value = 0.01, S3 Fig) and has a marginally significant interaction with total dietary protein (β = 0.38, p-value = 0.07) in urinary UA models. The SLC2A9-002 gene expression transcript significantly interacts with dietary total protein (β = -0.004, p-value = 0.04, S7 Table) and has marginally significant interactions with dietary animal and vegetable protein (β = -0.005, p-value = 0.07; β = -0.01, p-value = 0.06, respectively) in serum UA models (S7 Table).

Discussion

SNP-level genetic variation within SLC2A9 has previously been demonstrated to associate with urinary and serum UA phenotypes (30–36). Serum UA concentration and urinary UA excretion tend to be inversely related, since both are determined by the combination of UA reabsorption and secretion in the proximal tubule, and the SLC2A9 urate transporter is a key mediator of this process. Indeed the top SNP we identified, rs1209955, was previously identified by Doring and colleagues as significantly associated with serum UA levels [22].

Our analysis demonstrated that the SNP most strongly associated with urinary UA excretion, rs12509955, also significantly associated with SLC2A9 gene expression. This association between SNP-level variation within the SLC2A9 region and SLC2A9 gene expression has not previously been reported. The analysis also revealed other SNPs within the SLC2A9 region that were even more significantly associated with SLC2A9 gene expression than rs12509955.

One isoform of the SLCA29 transporter is expressed the apical membrane of the proximal tubule, while another is expressed in the basolateral membrane. The transporter exchanges UA for glucose and fructose resulting in net reabsorption of filtered UA. Genetic variation within SLC2A9 has previously been associated with serum UA concentrations and fractional erection of UA into the urine. In the current study we confirm and extend these observations, and demonstrate that total urinary UA excretion is also influenced by genetic variation in SLC2A9. This observation implies that changes in UA generation occur when urinary fractional excretion of UA decreases [51]. Glut9, the protein product of SLC2A9, is required for hepatocyte UA uptake. In mice, live specific knockout of SLC2A9 results in hyperuricemia and hyperuricosuria due to decreased degradation of UA by hepatic uricase. Humans lack uricase, but the current study suggests that genetic variability in SLCA29 alters overall UA production, perhaps via feedback at the level of the liver on enzymatic pathways that result in its generation. It is also conceivable that changes in availability of fructose could alter UA generation [52].

We previously hypothesized that an association between SLC2A9 genotype and UA phenotypes was mediated by SLC2A9 gene expression. SLC2A9 gene expression was not a significant predictor of urinary UA in our study. However, the gene expression data was obtained from transformed Beta-lymphocytes, and this cell type may not reflect gene expression levels present in the liver and kidneys, organs more directly responsible for UA metabolism and concentration. Nevertheless, the gene expression measured from transformed Beta-lymphocytes is thought to better reflect the proximal influences of genetic variation, since environmental modulators of gene expression are standardized by the uniform culture medium in which the cells are grown. It is quite possible that SLC2A9 gene expression in kidney cells would have more variability, since both environmental and genetic variability would be contributing to gene expression. Specifically, environmental conditions, such as dietary protein consumption, could have profound effects on tissue level gene expression, but have limited impact on gene expression levels measured in transformed Beta-lymphocytes, except perhaps through epigenetic DNA modifications that persist post-transformation. Overall, gene expression measurements are more reflective of the cumulative effects of multiple variants in the gene effecting gene expression. This may explain why highly significant associations between SNP genotypes and gene-level gene expression were observed in the current study. However, numerous steps between gene/gene expression and urinary UA excretion substantially weaken the pathway of statistical association.

The current study identified one significant transcript-level gene expression-by-diet interaction and several marginally significant interactions influencing urinary UA excretion, and several significant gene expression-by-dietary protein interactions in the context of serum UA. Since so many SNPs influenced SLC2A9 gene expression (Fig 2), it may be that gene expression values provide a better overall summary of the contribution common genetic variation in SLC2A9 makes to UA excretion, and could explain why gene expression displays a greater interaction with diet than rs12509955 genotype. In the future, identifying means to capture the cumulative effects of many genetic variations in the SLC2A9 region and their interactions with environmental factors such as dietary protein could facilitate efforts to assess the relative impact of genetic variability on serum and urinary UA levels.

In summary, the current study confirms that genetic variation within SLC2A9 plays an important role in determining urinary UA excretion. The most significant SNP associated with urinary UA (rs12509955), and many other SNPs in the gene region, were associated with SLC2A9 gene-level gene expression. Future investigations of mechanisms underlying associations between SLC2A9 genetic variation and urinary UA should include efforts to understand the influence of diet on SLC2A9 gene expression, urinary UA excretion, and NL risk.

Supporting Information

S1 Fig. Serum UA association results for 880 SNPs in the SLC2A9 gene region.

Left Y-axis:–log₁₀(p-value) from association between SNPs and serum UA, adjusted for age, sex, BMI, and urinary sodium, and accounting for sibship; Right Y-axis: SNP recombination rate based on HapMap hg18 CEU; X-axis: chromosomal location and gene regions; r² color code: degree of linkage disequilibrium with index (most strongly associated) SNP, rs11723439 (purple diamond).

https://doi.org/10.1371/journal.pone.0128593.s001

(PDF)

S2 Fig. Serum UA association results for 879 SNPs in the SLC2A9 gene region, conditioned on rs11723439.

Left Y-axis:–log₁₀(p-value) from association between SNPs and serum UA, adjusted for age, sex, BMI, rs11723439 and urinary sodium, and accounting for sibship; Right Y-axis: SNP recombination rate based on HapMap hg18 CEU; X-axis: chromosomal location and gene regions; r² color code: degree of linkage disequilibrium with index (most strongly associated) SNP, rs16894555 (purple diamond). Corrected p-value for rs16894555 is p = 0.044, indicating a second, independently SNP associated with serum uric acid in this sample.

https://doi.org/10.1371/journal.pone.0128593.s002

(PDF)

S3 Fig. SLC2A9-001 gene expression transcript interaction plot with dietary animal protein on urinary uric acid.

https://doi.org/10.1371/journal.pone.0128593.s003

(PDF)

S1 Table. Correlations between gene-level gene expression for uricosuric genes in the Genetic Epidemiology Network of Arteriopathy study.

ABCG2 ENSG00000118777, SLC17A1 ENSG00000124568, SLC17A3 ENSG00000124564, SLC22A12 ENSG00000197891, SLC2A9 ENSG00000109667, SLC2A9-001 ENST00000506583, SLC2A9-201 ENST00000309065, SLC2A9-002 ENST00000264784.

https://doi.org/10.1371/journal.pone.0128593.s004

(PDF)

S2 Table. Forty SLC2A9 SNPs with the lowest p-values in linear mixed models association with urinary uric acid in the Genetic Epidemiology Network of Arteriopathy study.

*corrected for multiple testing. Linear mixed models adjusted for age, sex, BMI and urinary sodium and adjusted for sibships.

https://doi.org/10.1371/journal.pone.0128593.s005

(PDF)

S3 Table. Forty SLC2A9 SNPs with the lowest p-values in linear mixed models, conditioned on rs12509955, association with urinary uric acid in the Genetic Epidemiology Network of Arteriopathy study.

*corrected for multiple testing. Linear mixed models adjusted for age, sex, BMI and urinary sodium and adjusted for sibships.

https://doi.org/10.1371/journal.pone.0128593.s006

(PDF)

S4 Table. Forty SLC2A9 SNPs with the lowest p-values in linear mixed models association with serum uric acid in the Genetic Epidemiology Network of Arteriopathy study.

*corrected for multiple testing. Linear mixed models adjusted for age, sex, BMI and urinary sodium and adjusted for sibships.

https://doi.org/10.1371/journal.pone.0128593.s007

(PDF)

S5 Table. Forty SLC2A9 SNPs with the lowest p-values in linear mixed models, conditioned on rs11723439, association with serum uric acid in the Genetic Epidemiology Network of Arteriopathy study.

*corrected for multiple testing. Linear mixed models adjusted for age, sex, BMI and urinary sodium and adjusted for sibships.

https://doi.org/10.1371/journal.pone.0128593.s008

(PDF)

S6 Table. Forty SLC2A9 SNPs with the lowest p-values in linear mixed models association with SLC2A9 gene expression in the Genetic Epidemiology Network of Arteriopathy study.

*corrected for multiple testing. Linear mixed models adjusted for BMI and urinary sodium and adjusted for sibships.

https://doi.org/10.1371/journal.pone.0128593.s009

(PDF)

S7 Table. Gene expression by dietary protein intake interaction associations for serum uric acid in the Genetic Epidemiology Network of Arteriopathy study.

*statistically significant association at α = 0.05. ABCG2 ENSG00000118777, SLC17A1 ENSG00000124568, SLC17A3 ENSG00000124564, SLC22A12 ENSG00000197891, SLC2A9 ENSG00000109667, SLC2A9-001 ENST00000506583, SLC2A9-002 ENST00000264784.

https://doi.org/10.1371/journal.pone.0128593.s010

(PDF)

S1 Text. Gene expression association with serum UA.

https://doi.org/10.1371/journal.pone.0128593.s011

(PDF)

Author Contributions

Conceived and designed the experiments: STT JCL. Analyzed the data: EBW ER WZ. Contributed reagents/materials/analysis tools: SLRK JCL STT. Wrote the paper: EBW ER JAS SLRK JCL STT.

References

1. Sierakowski R, Finlayson B, Landes RR, Finlayson CD, Sierakowski N. The frequency of urolithiasis in hospital discharge diagnoses in the United States. Investigative urology. 1978;15(6):438–41. Epub 1978/05/01. pmid:649290.
- View Article
- PubMed/NCBI
- Google Scholar
2. Scales CD Jr, Smith AC, Hanley JM, Saigal CS, Urologic Diseases in America P. Prevalence of kidney stones in the United States. European urology. 2012;62(1):160–5. Epub 2012/04/14. pmid:22498635; PubMed Central PMCID: PMC3362665.
- View Article
- PubMed/NCBI
- Google Scholar
3. Mandel NS, Mandel GS. Urinary tract stone disease in the United States veteran population. I. Geographical frequency of occurrence. The Journal of urology. 1989;142(6):1513–5. Epub 1989/12/01. pmid:2585626.
- View Article
- PubMed/NCBI
- Google Scholar
4. Curhan GC, Willett WC, Rimm EB, Stampfer MJ. Family history and risk of kidney stones. Journal of the American Society of Nephrology: JASN. 1997;8(10):1568–73. Epub 1997/10/23. pmid:9335385.
- View Article
- PubMed/NCBI
- Google Scholar
5. Griffin DG. A review of the heritability of idiopathic nephrolithiasis. Journal of clinical pathology. 2004;57(8):793–6. Epub 2004/07/29. pmid:15280397; PubMed Central PMCID: PMC1770398.
- View Article
- PubMed/NCBI
- Google Scholar
6. Kuo CF, Grainge MJ, See LC, Yu KH, Luo SF, Valdes AM, et al. Familial aggregation of gout and relative genetic and environmental contributions: a nationwide population study in Taiwan. Annals of the rheumatic diseases. 2013. Epub 2013/11/23. https://doi.org/10.1136/annrheumdis-2013-204067 pmid:24265412.
7. McGeown MG. Heredity in renal stone disease. Clinical science. 1960;19:465–71. Epub 1960/08/01. pmid:13773937.
- View Article
- PubMed/NCBI
- Google Scholar
8. Resnick M, Pridgen DB, Goodman HO. Genetic predisposition to formation of calcium oxalate renal calculi. The New England journal of medicine. 1968;278(24):1313–8. Epub 1968/06/13. pmid:5648597.
- View Article
- PubMed/NCBI
- Google Scholar
9. Trinchieri A, Mandressi A, Luongo P, Coppi F, Pisani E. Familial aggregation of renal calcium stone disease. The Journal of urology. 1988;139(3):478–81. Epub 1988/03/01. pmid:3343730.
- View Article
- PubMed/NCBI
- Google Scholar
10. Park C, Ha YS, Kim YJ, Yun SJ, Lee SC, Kim WJ. Comparison of Metabolic Risk Factors in Urolithiasis Patients according to Family History. Korean journal of urology. 2010;51(1):50–3. Epub 2010/04/24. pmid:20414411; PubMed Central PMCID: PMC2855461.
- View Article
- PubMed/NCBI
- Google Scholar
11. Leonetti F, Dussol B, Berthezene P, Thirion X, Berland Y. Dietary and urinary risk factors for stones in idiopathic calcium stone formers compared with healthy subjects. Nephrology, dialysis, transplantation: official publication of the European Dialysis and Transplant Association—European Renal Association. 1998;13(3):617–22. Epub 1998/04/29. pmid:9550636.
- View Article
- PubMed/NCBI
- Google Scholar
12. Gutman AB, Yu TF. Uric acid nephrolithiasis. The American journal of medicine. 1968;45(5):756–79. Epub 1968/11/01. pmid:4879835.
- View Article
- PubMed/NCBI
- Google Scholar
13. Gault MH, Chafe L. Relationship of frequency, age, sex, stone weight and composition in 15,624 stones: comparison of resutls for 1980 to 1983 and 1995 to 1998. The Journal of urology. 2000;164(2):302–7. Epub 2000/07/14. pmid:10893570.
- View Article
- PubMed/NCBI
- Google Scholar
14. Grover PK, Ryall RL, Marshall VR. Dissolved urate promotes calcium oxalate crystallization: epitaxy is not the cause. Clin Sci (Lond). 1993;85(3):303–7. Epub 1993/09/01. pmid:8403802.
- View Article
- PubMed/NCBI
- Google Scholar
15. Coe FL. Hyperuricosuric calcium oxalate nephrolithiasis. Kidney international. 1978;13(5):418–26. Epub 1978/05/01. pmid:661071.
- View Article
- PubMed/NCBI
- Google Scholar
16. Kirejczyk JK, Porowski T, Filonowicz R, Kazberuk A, Stefanowicz M, Wasilewska A, et al. An association between kidney stone composition and urinary metabolic disturbances in children. Journal of pediatric urology. 2014;10(1):130–5. Epub 2013/08/21. pmid:23953243.
- View Article
- PubMed/NCBI
- Google Scholar
17. Pais VM Jr, Holmes RP, Assimos DG. Effect of dietary control of urinary uric acid excretion in calcium oxalate stone formers and non-stone-forming controls. Journal of endourology / Endourological Society. 2007;21(2):232–5. Epub 2007/03/07. pmid:17338626; PubMed Central PMCID: PMC2705619.
- View Article
- PubMed/NCBI
- Google Scholar
18. Pak CY, Waters O, Arnold L, Holt K, Cox C, Barilla D. Mechanism for calcium urolithiasis among patients with hyperuricosuria: supersaturation of urine with respect to monosodium urate. The Journal of clinical investigation. 1977;59(3):426–31. Epub 1977/03/01. pmid:14173; PubMed Central PMCID: PMC333378.
- View Article
- PubMed/NCBI
- Google Scholar
19. Robertson WG, Peacock M, Heyburn PJ, Marshall DH, Clark PB. Risk factors in calcium stone disease of the urinary tract. British journal of urology. 1978;50(7):449–54. Epub 1978/12/01. pmid:37975.
- View Article
- PubMed/NCBI
- Google Scholar
20. Mosli HA, Mosli HH, Kamal WK. Kidney stone composition in overweight and obese patients: a preliminary report. Research and reports in urology. 2013;5:11–5. Epub 2014/01/09. pmid:24400230; PubMed Central PMCID: PMC3826902.
- View Article
- PubMed/NCBI
- Google Scholar
21. Augustin R, Carayannopoulos MO, Dowd LO, Phay JE, Moley JF, Moley KH. Identification and characterization of human glucose transporter-like protein-9 (GLUT9): alternative splicing alters trafficking. The Journal of biological chemistry. 2004;279(16):16229–36. Epub 2004/01/24. pmid:14739288.
- View Article
- PubMed/NCBI
- Google Scholar
22. Doring A, Gieger C, Mehta D, Gohlke H, Prokisch H, Coassin S, et al. SLC2A9 influences uric acid concentrations with pronounced sex-specific effects. Nature genetics. 2008;40(4):430–6. Epub 2008/03/11. pmid:18327256.
- View Article
- PubMed/NCBI
- Google Scholar
23. Witkowska K, Smith KM, Yao SY, Ng AM, O'Neill D, Karpinski E, et al. Human SLC2A9a and SLC2A9b isoforms mediate electrogenic transport of urate with different characteristics in the presence of hexoses. American journal of physiology Renal physiology. 2012;303(4):F527–39. Epub 2012/06/01. pmid:22647630; PubMed Central PMCID: PMC3423118.
- View Article
- PubMed/NCBI
- Google Scholar
24. Anzai N, Jutabha P, Amonpatumrat-Takahashi S, Sakurai H. Recent advances in renal urate transport: characterization of candidate transporters indicated by genome-wide association studies. Clinical and experimental nephrology. 2012;16(1):89–95. Epub 2011/11/01. pmid:22038265.
- View Article
- PubMed/NCBI
- Google Scholar
25. Roch-Ramel F, Guisan B. Renal Transport of Urate in Humans. News in physiological sciences: an international journal of physiology produced jointly by the International Union of Physiological Sciences and the American Physiological Society. 1999;14:80–4. Epub 2001/06/08. pmid:11390825.
- View Article
- PubMed/NCBI
- Google Scholar
26. Dinour D, Gray NK, Ganon L, Knox AJ, Shalev H, Sela BA, et al. Two novel homozygous SLC2A9 mutations cause renal hypouricemia type 2. Nephrology, dialysis, transplantation: official publication of the European Dialysis and Transplant Association—European Renal Association. 2012;27(3):1035–41. Epub 2011/08/04. pmid:21810765.
- View Article
- PubMed/NCBI
- Google Scholar
27. Kawamura Y, Matsuo H, Chiba T, Nagamori S, Nakayama A, Inoue H, et al. Pathogenic GLUT9 mutations causing renal hypouricemia type 2 (RHUC2). Nucleosides, nucleotides & nucleic acids. 2011;30(12):1105–11. Epub 2011/12/03. pmid:22132964.
- View Article
- PubMed/NCBI
- Google Scholar
28. Stiburkova B, Ichida K, Sebesta I. Novel homozygous insertion in SLC2A9 gene caused renal hypouricemia. Molecular genetics and metabolism. 2011;102(4):430–5. Epub 2011/01/25. pmid:21256783.
- View Article
- PubMed/NCBI
- Google Scholar
29. Sull JW, Park EJ, Lee M, Jee SH. Effects of SLC2A9 variants on uric acid levels in a Korean population. Rheumatology international. 2013;33(1):19–23. Epub 2012/01/04. pmid:22212408.
- View Article
- PubMed/NCBI
- Google Scholar
30. Dehghan A, Kottgen A, Yang Q, Hwang SJ, Kao WL, Rivadeneira F, et al. Association of three genetic loci with uric acid concentration and risk of gout: a genome-wide association study. Lancet. 2008;372(9654):1953–61. Epub 2008/10/07. pmid:18834626; PubMed Central PMCID: PMC2803340.
- View Article
- PubMed/NCBI
- Google Scholar
31. Hamajima N, Okada R, Kawai S, Hishida A, Morita E, Yin G, et al. Significant association of serum uric acid levels with SLC2A9 rs11722228 among a Japanese population. Molecular genetics and metabolism. 2011;103(4):378–82. Epub 2011/04/23. pmid:21511506.
- View Article
- PubMed/NCBI
- Google Scholar
32. Kottgen A, Albrecht E, Teumer A, Vitart V, Krumsiek J, Hundertmark C, et al. Genome-wide association analyses identify 18 new loci associated with serum urate concentrations. Nature genetics. 2013;45(2):145–54. Epub 2012/12/25. pmid:23263486; PubMed Central PMCID: PMC3663712.
- View Article
- PubMed/NCBI
- Google Scholar
33. Polasek O, Gunjaca G, Kolcic I, Zgaga L, Dzijan S, Smolic R, et al. Association of nephrolithiasis and gene for glucose transporter type 9 (SLC2A9): study of 145 patients. Croatian medical journal. 2010;51(1):48–53. Epub 2010/02/18. pmid:20162745; PubMed Central PMCID: PMC2829176.
- View Article
- PubMed/NCBI
- Google Scholar
34. Rule AD, de Andrade M, Matsumoto M, Mosley TH, Kardia S, Turner ST. Association between SLC2A9 transporter gene variants and uric acid phenotypes in African American and white families. Rheumatology. 2011;50(5):871–8. Epub 2010/12/28. pmid:21186168; PubMed Central PMCID: PMC3077913.
- View Article
- PubMed/NCBI
- Google Scholar
35. Voruganti VS, Franceschini N, Haack K, Laston S, Maccluer JW, Umans JG, et al. Replication of the effect of SLC2A9 genetic variation on serum uric acid levels in American Indians. European journal of human genetics: EJHG. 2013. Epub 2013/12/05. https://doi.org/10.1038/ejhg.2013.264 pmid:24301058.
36. Voruganti VS, Kent JW Jr, Debnath S, Cole SA, Haack K, Goring HH, et al. Genome-wide association analysis confirms and extends the association of SLC2A9 with serum uric acid levels to Mexican Americans. Frontiers in genetics. 2013;4:279. Epub 2014/01/01. pmid:24379826; PubMed Central PMCID: PMC3863993.
- View Article
- PubMed/NCBI
- Google Scholar
37. Enomoto A, Kimura H, Chairoungdua A, Shigeta Y, Jutabha P, Cha SH, et al. Molecular identification of a renal urate anion exchanger that regulates blood urate levels. Nature. 2002;417(6887):447–52. pmid:12024214.
- View Article
- PubMed/NCBI
- Google Scholar
38. Lipkowitz MS. Regulation of uric acid excretion by the kidney. Current rheumatology reports. 2012;14(2):179–88. pmid:22359229.
- View Article
- PubMed/NCBI
- Google Scholar
39. Reginato AM, Mount DB, Yang I, Choi HK. The genetics of hyperuricaemia and gout. Nature reviews Rheumatology. 2012;8(10):610–21. pmid:22945592; PubMed Central PMCID: PMC3645862.
- View Article
- PubMed/NCBI
- Google Scholar
40. Wright AF, Rudan I, Hastie ND, Campbell H. A 'complexity' of urate transporters. Kidney international. 2010;78(5):446–52. pmid:20613716.
- View Article
- PubMed/NCBI
- Google Scholar
41. Giannini S, Nobile M, Sartori L, Dalle Carbonare L, Ciuffreda M, Corro P, et al. Acute effects of moderate dietary protein restriction in patients with idiopathic hypercalciuria and calcium nephrolithiasis. The American journal of clinical nutrition. 1999;69(2):267–71. Epub 1999/02/16. pmid:9989691.
- View Article
- PubMed/NCBI
- Google Scholar
42. Tracy CR, Best S, Bagrodia A, Poindexter JR, Adams-Huet B, Sakhaee K, et al. Animal Protein and the Risk of Kidney Stones: A Comparative Metabolic Study of Animal Protein Sources. The Journal of urology. 2014. Epub 2014/02/13. https://doi.org/10.1016/j.juro.2014.01.093 pmid:24518789.
43. Heilberg IP, Goldfarb DS. Optimum nutrition for kidney stone disease. Advances in chronic kidney disease. 2013;20(2):165–74. Epub 2013/02/27. pmid:23439376.
- View Article
- PubMed/NCBI
- Google Scholar
44. Reddy ST, Wang CY, Sakhaee K, Brinkley L, Pak CY. Effect of low-carbohydrate high-protein diets on acid-base balance, stone-forming propensity, and calcium metabolism. American journal of kidney diseases: the official journal of the National Kidney Foundation. 2002;40(2):265–74. Epub 2002/07/31. pmid:12148098.
- View Article
- PubMed/NCBI
- Google Scholar
45. Lieske JC, Turner ST, Edeh SN, Smith JA, Kardia SL. Heritability of Urinary Traits That Contribute to Nephrolithiasis. Clinical journal of the American Society of Nephrology: CJASN. 2014. Epub 2014/03/01. https://doi.org/10.2215/CJN.08210813 pmid:24578335.
46. Patterson RE, Kristal AR, Tinker LF, Carter RA, Bolton MP, Agurs-Collins T. Measurement characteristics of the Women's Health Initiative food frequency questionnaire. Annals of epidemiology. 1999;9(3):178–87. pmid:10192650.
- View Article
- PubMed/NCBI
- Google Scholar
47. Base SAS(R) Cary, NC: SAS Institute Inc.; 2011.
48. Pruim RJ, Welch RP, Sanna S, Teslovich TM, Chines PS, Gliedt TP, et al. LocusZoom: regional visualization of genome-wide association scan results. Bioinformatics. 2010;26(18):2336–7. pmid:20634204; PubMed Central PMCID: PMC2935401.
- View Article
- PubMed/NCBI
- Google Scholar
49. Gao X, Starmer J, Martin ER. A multiple testing correction method for genetic association studies using correlated single nucleotide polymorphisms. Genetic epidemiology. 2008;32(4):361–9. pmid:18271029.
- View Article
- PubMed/NCBI
- Google Scholar
50. Cunningham F, Amode MR, Barrell D, Beal K, Billis K, Brent S, et al. Ensembl 2015. Nucleic Acids Research. 2014. https://doi.org/10.1093/nar/gku1010
51. Preitner F, Laverriere-Loss A, Metref S, Da Costa A, Moret C, Rotman S, et al. Urate-induced acute renal failure and chronic inflammation in liver-specific Glut9 knockout mice. American journal of physiology Renal physiology. 2013;305(5):F786–95. pmid:23804456.
- View Article
- PubMed/NCBI
- Google Scholar
52. Khitan Z, Kim DH. Fructose: a key factor in the development of metabolic syndrome and hypertension. Journal of nutrition and metabolism. 2013;2013:682673. pmid:23762544; PubMed Central PMCID: PMC3677638.
- View Article
- PubMed/NCBI
- Google Scholar

[ref1] 1. Sierakowski R, Finlayson B, Landes RR, Finlayson CD, Sierakowski N. The frequency of urolithiasis in hospital discharge diagnoses in the United States. Investigative urology. 1978;15(6):438–41. Epub 1978/05/01. pmid:649290.
View Article
PubMed/NCBI
Google Scholar

[2] View Article

[3] PubMed/NCBI

[4] Google Scholar

[ref2] 2. Scales CD Jr, Smith AC, Hanley JM, Saigal CS, Urologic Diseases in America P. Prevalence of kidney stones in the United States. European urology. 2012;62(1):160–5. Epub 2012/04/14. pmid:22498635; PubMed Central PMCID: PMC3362665.
View Article
PubMed/NCBI
Google Scholar

[6] View Article

[7] PubMed/NCBI

[8] Google Scholar

[ref3] 3. Mandel NS, Mandel GS. Urinary tract stone disease in the United States veteran population. I. Geographical frequency of occurrence. The Journal of urology. 1989;142(6):1513–5. Epub 1989/12/01. pmid:2585626.
View Article
PubMed/NCBI
Google Scholar

[10] View Article

[11] PubMed/NCBI

[12] Google Scholar

[ref4] 4. Curhan GC, Willett WC, Rimm EB, Stampfer MJ. Family history and risk of kidney stones. Journal of the American Society of Nephrology: JASN. 1997;8(10):1568–73. Epub 1997/10/23. pmid:9335385.
View Article
PubMed/NCBI
Google Scholar

[14] View Article

[15] PubMed/NCBI

[16] Google Scholar

[ref5] 5. Griffin DG. A review of the heritability of idiopathic nephrolithiasis. Journal of clinical pathology. 2004;57(8):793–6. Epub 2004/07/29. pmid:15280397; PubMed Central PMCID: PMC1770398.
View Article
PubMed/NCBI
Google Scholar

[18] View Article

[19] PubMed/NCBI

[20] Google Scholar

[ref6] 6. Kuo CF, Grainge MJ, See LC, Yu KH, Luo SF, Valdes AM, et al. Familial aggregation of gout and relative genetic and environmental contributions: a nationwide population study in Taiwan. Annals of the rheumatic diseases. 2013. Epub 2013/11/23. https://doi.org/10.1136/annrheumdis-2013-204067 pmid:24265412.

[ref7] 7. McGeown MG. Heredity in renal stone disease. Clinical science. 1960;19:465–71. Epub 1960/08/01. pmid:13773937.
View Article
PubMed/NCBI
Google Scholar

[23] View Article

[24] PubMed/NCBI

[25] Google Scholar

[ref8] 8. Resnick M, Pridgen DB, Goodman HO. Genetic predisposition to formation of calcium oxalate renal calculi. The New England journal of medicine. 1968;278(24):1313–8. Epub 1968/06/13. pmid:5648597.
View Article
PubMed/NCBI
Google Scholar

[27] View Article

[28] PubMed/NCBI

[29] Google Scholar

[ref9] 9. Trinchieri A, Mandressi A, Luongo P, Coppi F, Pisani E. Familial aggregation of renal calcium stone disease. The Journal of urology. 1988;139(3):478–81. Epub 1988/03/01. pmid:3343730.
View Article
PubMed/NCBI
Google Scholar

[31] View Article

[32] PubMed/NCBI

[33] Google Scholar

[ref10] 10. Park C, Ha YS, Kim YJ, Yun SJ, Lee SC, Kim WJ. Comparison of Metabolic Risk Factors in Urolithiasis Patients according to Family History. Korean journal of urology. 2010;51(1):50–3. Epub 2010/04/24. pmid:20414411; PubMed Central PMCID: PMC2855461.
View Article
PubMed/NCBI
Google Scholar

[35] View Article

[36] PubMed/NCBI

[37] Google Scholar

[ref11] 11. Leonetti F, Dussol B, Berthezene P, Thirion X, Berland Y. Dietary and urinary risk factors for stones in idiopathic calcium stone formers compared with healthy subjects. Nephrology, dialysis, transplantation: official publication of the European Dialysis and Transplant Association—European Renal Association. 1998;13(3):617–22. Epub 1998/04/29. pmid:9550636.
View Article
PubMed/NCBI
Google Scholar

[39] View Article

[40] PubMed/NCBI

[41] Google Scholar

[ref12] 12. Gutman AB, Yu TF. Uric acid nephrolithiasis. The American journal of medicine. 1968;45(5):756–79. Epub 1968/11/01. pmid:4879835.
View Article
PubMed/NCBI
Google Scholar

[43] View Article

[44] PubMed/NCBI

[45] Google Scholar

[ref13] 13. Gault MH, Chafe L. Relationship of frequency, age, sex, stone weight and composition in 15,624 stones: comparison of resutls for 1980 to 1983 and 1995 to 1998. The Journal of urology. 2000;164(2):302–7. Epub 2000/07/14. pmid:10893570.
View Article
PubMed/NCBI
Google Scholar

[47] View Article

[48] PubMed/NCBI

[49] Google Scholar

[ref14] 14. Grover PK, Ryall RL, Marshall VR. Dissolved urate promotes calcium oxalate crystallization: epitaxy is not the cause. Clin Sci (Lond). 1993;85(3):303–7. Epub 1993/09/01. pmid:8403802.
View Article
PubMed/NCBI
Google Scholar

[51] View Article

[52] PubMed/NCBI

[53] Google Scholar

[ref15] 15. Coe FL. Hyperuricosuric calcium oxalate nephrolithiasis. Kidney international. 1978;13(5):418–26. Epub 1978/05/01. pmid:661071.
View Article
PubMed/NCBI
Google Scholar

[55] View Article

[56] PubMed/NCBI

[57] Google Scholar

[ref16] 16. Kirejczyk JK, Porowski T, Filonowicz R, Kazberuk A, Stefanowicz M, Wasilewska A, et al. An association between kidney stone composition and urinary metabolic disturbances in children. Journal of pediatric urology. 2014;10(1):130–5. Epub 2013/08/21. pmid:23953243.
View Article
PubMed/NCBI
Google Scholar

[59] View Article

[60] PubMed/NCBI

[61] Google Scholar

[ref17] 17. Pais VM Jr, Holmes RP, Assimos DG. Effect of dietary control of urinary uric acid excretion in calcium oxalate stone formers and non-stone-forming controls. Journal of endourology / Endourological Society. 2007;21(2):232–5. Epub 2007/03/07. pmid:17338626; PubMed Central PMCID: PMC2705619.
View Article
PubMed/NCBI
Google Scholar

[63] View Article

[64] PubMed/NCBI

[65] Google Scholar

[ref18] 18. Pak CY, Waters O, Arnold L, Holt K, Cox C, Barilla D. Mechanism for calcium urolithiasis among patients with hyperuricosuria: supersaturation of urine with respect to monosodium urate. The Journal of clinical investigation. 1977;59(3):426–31. Epub 1977/03/01. pmid:14173; PubMed Central PMCID: PMC333378.
View Article
PubMed/NCBI
Google Scholar

[67] View Article

[68] PubMed/NCBI

[69] Google Scholar

[ref19] 19. Robertson WG, Peacock M, Heyburn PJ, Marshall DH, Clark PB. Risk factors in calcium stone disease of the urinary tract. British journal of urology. 1978;50(7):449–54. Epub 1978/12/01. pmid:37975.
View Article
PubMed/NCBI
Google Scholar

[71] View Article

[72] PubMed/NCBI

[73] Google Scholar

[ref20] 20. Mosli HA, Mosli HH, Kamal WK. Kidney stone composition in overweight and obese patients: a preliminary report. Research and reports in urology. 2013;5:11–5. Epub 2014/01/09. pmid:24400230; PubMed Central PMCID: PMC3826902.
View Article
PubMed/NCBI
Google Scholar

[75] View Article

[76] PubMed/NCBI

[77] Google Scholar

[ref21] 21. Augustin R, Carayannopoulos MO, Dowd LO, Phay JE, Moley JF, Moley KH. Identification and characterization of human glucose transporter-like protein-9 (GLUT9): alternative splicing alters trafficking. The Journal of biological chemistry. 2004;279(16):16229–36. Epub 2004/01/24. pmid:14739288.
View Article
PubMed/NCBI
Google Scholar

[79] View Article

[80] PubMed/NCBI

[81] Google Scholar

[ref22] 22. Doring A, Gieger C, Mehta D, Gohlke H, Prokisch H, Coassin S, et al. SLC2A9 influences uric acid concentrations with pronounced sex-specific effects. Nature genetics. 2008;40(4):430–6. Epub 2008/03/11. pmid:18327256.
View Article
PubMed/NCBI
Google Scholar

[83] View Article

[84] PubMed/NCBI

[85] Google Scholar

[ref23] 23. Witkowska K, Smith KM, Yao SY, Ng AM, O'Neill D, Karpinski E, et al. Human SLC2A9a and SLC2A9b isoforms mediate electrogenic transport of urate with different characteristics in the presence of hexoses. American journal of physiology Renal physiology. 2012;303(4):F527–39. Epub 2012/06/01. pmid:22647630; PubMed Central PMCID: PMC3423118.
View Article
PubMed/NCBI
Google Scholar

[87] View Article

[88] PubMed/NCBI

[89] Google Scholar

[ref24] 24. Anzai N, Jutabha P, Amonpatumrat-Takahashi S, Sakurai H. Recent advances in renal urate transport: characterization of candidate transporters indicated by genome-wide association studies. Clinical and experimental nephrology. 2012;16(1):89–95. Epub 2011/11/01. pmid:22038265.
View Article
PubMed/NCBI
Google Scholar

[91] View Article

[92] PubMed/NCBI

[93] Google Scholar

[ref25] 25. Roch-Ramel F, Guisan B. Renal Transport of Urate in Humans. News in physiological sciences: an international journal of physiology produced jointly by the International Union of Physiological Sciences and the American Physiological Society. 1999;14:80–4. Epub 2001/06/08. pmid:11390825.
View Article
PubMed/NCBI
Google Scholar

[95] View Article

[96] PubMed/NCBI

[97] Google Scholar

[ref26] 26. Dinour D, Gray NK, Ganon L, Knox AJ, Shalev H, Sela BA, et al. Two novel homozygous SLC2A9 mutations cause renal hypouricemia type 2. Nephrology, dialysis, transplantation: official publication of the European Dialysis and Transplant Association—European Renal Association. 2012;27(3):1035–41. Epub 2011/08/04. pmid:21810765.
View Article
PubMed/NCBI
Google Scholar

[99] View Article

[100] PubMed/NCBI

[101] Google Scholar

[ref27] 27. Kawamura Y, Matsuo H, Chiba T, Nagamori S, Nakayama A, Inoue H, et al. Pathogenic GLUT9 mutations causing renal hypouricemia type 2 (RHUC2). Nucleosides, nucleotides & nucleic acids. 2011;30(12):1105–11. Epub 2011/12/03. pmid:22132964.
View Article
PubMed/NCBI
Google Scholar

[103] View Article

[104] PubMed/NCBI

[105] Google Scholar

[ref28] 28. Stiburkova B, Ichida K, Sebesta I. Novel homozygous insertion in SLC2A9 gene caused renal hypouricemia. Molecular genetics and metabolism. 2011;102(4):430–5. Epub 2011/01/25. pmid:21256783.
View Article
PubMed/NCBI
Google Scholar

[107] View Article

[108] PubMed/NCBI

[109] Google Scholar

[ref29] 29. Sull JW, Park EJ, Lee M, Jee SH. Effects of SLC2A9 variants on uric acid levels in a Korean population. Rheumatology international. 2013;33(1):19–23. Epub 2012/01/04. pmid:22212408.
View Article
PubMed/NCBI
Google Scholar

[111] View Article

[112] PubMed/NCBI

[113] Google Scholar

[ref30] 30. Dehghan A, Kottgen A, Yang Q, Hwang SJ, Kao WL, Rivadeneira F, et al. Association of three genetic loci with uric acid concentration and risk of gout: a genome-wide association study. Lancet. 2008;372(9654):1953–61. Epub 2008/10/07. pmid:18834626; PubMed Central PMCID: PMC2803340.
View Article
PubMed/NCBI
Google Scholar

[115] View Article

[116] PubMed/NCBI

[117] Google Scholar

[ref31] 31. Hamajima N, Okada R, Kawai S, Hishida A, Morita E, Yin G, et al. Significant association of serum uric acid levels with SLC2A9 rs11722228 among a Japanese population. Molecular genetics and metabolism. 2011;103(4):378–82. Epub 2011/04/23. pmid:21511506.
View Article
PubMed/NCBI
Google Scholar

[119] View Article

[120] PubMed/NCBI

[121] Google Scholar

[ref32] 32. Kottgen A, Albrecht E, Teumer A, Vitart V, Krumsiek J, Hundertmark C, et al. Genome-wide association analyses identify 18 new loci associated with serum urate concentrations. Nature genetics. 2013;45(2):145–54. Epub 2012/12/25. pmid:23263486; PubMed Central PMCID: PMC3663712.
View Article
PubMed/NCBI
Google Scholar

[123] View Article

[124] PubMed/NCBI

[125] Google Scholar

[ref33] 33. Polasek O, Gunjaca G, Kolcic I, Zgaga L, Dzijan S, Smolic R, et al. Association of nephrolithiasis and gene for glucose transporter type 9 (SLC2A9): study of 145 patients. Croatian medical journal. 2010;51(1):48–53. Epub 2010/02/18. pmid:20162745; PubMed Central PMCID: PMC2829176.
View Article
PubMed/NCBI
Google Scholar

[127] View Article

[128] PubMed/NCBI

[129] Google Scholar

[ref34] 34. Rule AD, de Andrade M, Matsumoto M, Mosley TH, Kardia S, Turner ST. Association between SLC2A9 transporter gene variants and uric acid phenotypes in African American and white families. Rheumatology. 2011;50(5):871–8. Epub 2010/12/28. pmid:21186168; PubMed Central PMCID: PMC3077913.
View Article
PubMed/NCBI
Google Scholar

[131] View Article

[132] PubMed/NCBI

[133] Google Scholar

[ref35] 35. Voruganti VS, Franceschini N, Haack K, Laston S, Maccluer JW, Umans JG, et al. Replication of the effect of SLC2A9 genetic variation on serum uric acid levels in American Indians. European journal of human genetics: EJHG. 2013. Epub 2013/12/05. https://doi.org/10.1038/ejhg.2013.264 pmid:24301058.

[ref36] 36. Voruganti VS, Kent JW Jr, Debnath S, Cole SA, Haack K, Goring HH, et al. Genome-wide association analysis confirms and extends the association of SLC2A9 with serum uric acid levels to Mexican Americans. Frontiers in genetics. 2013;4:279. Epub 2014/01/01. pmid:24379826; PubMed Central PMCID: PMC3863993.
View Article
PubMed/NCBI
Google Scholar

[136] View Article

[137] PubMed/NCBI

[138] Google Scholar

[ref37] 37. Enomoto A, Kimura H, Chairoungdua A, Shigeta Y, Jutabha P, Cha SH, et al. Molecular identification of a renal urate anion exchanger that regulates blood urate levels. Nature. 2002;417(6887):447–52. pmid:12024214.
View Article
PubMed/NCBI
Google Scholar

[140] View Article

[141] PubMed/NCBI

[142] Google Scholar

[ref38] 38. Lipkowitz MS. Regulation of uric acid excretion by the kidney. Current rheumatology reports. 2012;14(2):179–88. pmid:22359229.
View Article
PubMed/NCBI
Google Scholar

[144] View Article

[145] PubMed/NCBI

[146] Google Scholar

[ref39] 39. Reginato AM, Mount DB, Yang I, Choi HK. The genetics of hyperuricaemia and gout. Nature reviews Rheumatology. 2012;8(10):610–21. pmid:22945592; PubMed Central PMCID: PMC3645862.
View Article
PubMed/NCBI
Google Scholar

[148] View Article

[149] PubMed/NCBI

[150] Google Scholar

[ref40] 40. Wright AF, Rudan I, Hastie ND, Campbell H. A 'complexity' of urate transporters. Kidney international. 2010;78(5):446–52. pmid:20613716.
View Article
PubMed/NCBI
Google Scholar

[152] View Article

[153] PubMed/NCBI

[154] Google Scholar

[ref41] 41. Giannini S, Nobile M, Sartori L, Dalle Carbonare L, Ciuffreda M, Corro P, et al. Acute effects of moderate dietary protein restriction in patients with idiopathic hypercalciuria and calcium nephrolithiasis. The American journal of clinical nutrition. 1999;69(2):267–71. Epub 1999/02/16. pmid:9989691.
View Article
PubMed/NCBI
Google Scholar

[156] View Article

[157] PubMed/NCBI

[158] Google Scholar

[ref42] 42. Tracy CR, Best S, Bagrodia A, Poindexter JR, Adams-Huet B, Sakhaee K, et al. Animal Protein and the Risk of Kidney Stones: A Comparative Metabolic Study of Animal Protein Sources. The Journal of urology. 2014. Epub 2014/02/13. https://doi.org/10.1016/j.juro.2014.01.093 pmid:24518789.

[ref43] 43. Heilberg IP, Goldfarb DS. Optimum nutrition for kidney stone disease. Advances in chronic kidney disease. 2013;20(2):165–74. Epub 2013/02/27. pmid:23439376.
View Article
PubMed/NCBI
Google Scholar

[161] View Article

[162] PubMed/NCBI

[163] Google Scholar

[ref44] 44. Reddy ST, Wang CY, Sakhaee K, Brinkley L, Pak CY. Effect of low-carbohydrate high-protein diets on acid-base balance, stone-forming propensity, and calcium metabolism. American journal of kidney diseases: the official journal of the National Kidney Foundation. 2002;40(2):265–74. Epub 2002/07/31. pmid:12148098.
View Article
PubMed/NCBI
Google Scholar

[165] View Article

[166] PubMed/NCBI

[167] Google Scholar

[ref45] 45. Lieske JC, Turner ST, Edeh SN, Smith JA, Kardia SL. Heritability of Urinary Traits That Contribute to Nephrolithiasis. Clinical journal of the American Society of Nephrology: CJASN. 2014. Epub 2014/03/01. https://doi.org/10.2215/CJN.08210813 pmid:24578335.

[ref46] 46. Patterson RE, Kristal AR, Tinker LF, Carter RA, Bolton MP, Agurs-Collins T. Measurement characteristics of the Women's Health Initiative food frequency questionnaire. Annals of epidemiology. 1999;9(3):178–87. pmid:10192650.
View Article
PubMed/NCBI
Google Scholar

[170] View Article

[171] PubMed/NCBI

[172] Google Scholar

[ref47] 47. Base SAS(R) Cary, NC: SAS Institute Inc.; 2011.

[ref48] 48. Pruim RJ, Welch RP, Sanna S, Teslovich TM, Chines PS, Gliedt TP, et al. LocusZoom: regional visualization of genome-wide association scan results. Bioinformatics. 2010;26(18):2336–7. pmid:20634204; PubMed Central PMCID: PMC2935401.
View Article
PubMed/NCBI
Google Scholar

[175] View Article

[176] PubMed/NCBI

[177] Google Scholar

[ref49] 49. Gao X, Starmer J, Martin ER. A multiple testing correction method for genetic association studies using correlated single nucleotide polymorphisms. Genetic epidemiology. 2008;32(4):361–9. pmid:18271029.
View Article
PubMed/NCBI
Google Scholar

[179] View Article

[180] PubMed/NCBI

[181] Google Scholar

[ref50] 50. Cunningham F, Amode MR, Barrell D, Beal K, Billis K, Brent S, et al. Ensembl 2015. Nucleic Acids Research. 2014. https://doi.org/10.1093/nar/gku1010

[ref51] 51. Preitner F, Laverriere-Loss A, Metref S, Da Costa A, Moret C, Rotman S, et al. Urate-induced acute renal failure and chronic inflammation in liver-specific Glut9 knockout mice. American journal of physiology Renal physiology. 2013;305(5):F786–95. pmid:23804456.
View Article
PubMed/NCBI
Google Scholar

[184] View Article

[185] PubMed/NCBI

[186] Google Scholar

[ref52] 52. Khitan Z, Kim DH. Fructose: a key factor in the development of metabolic syndrome and hypertension. Journal of nutrition and metabolism. 2013;2013:682673. pmid:23762544; PubMed Central PMCID: PMC3677638.
View Article
PubMed/NCBI
Google Scholar

[188] View Article

[189] PubMed/NCBI

[190] Google Scholar

Figures

Abstract

Objectives

Methods

Results

Conclusion

Introduction

Methods

Study Population

Ethics Statement

Urine Collection

Food frequency

Genotyping and Imputation

Gene Expression Assessment and Quality Control

Statistical Methods

Results

Descriptive Statistics

SLC2A9 SNP association with urinary UA excretion

SLC2A9 SNP association with gene expression

Gene expression association with urinary UA

Gene-by-dietary protein interactions

Discussion

Supporting Information

S1 Fig. Serum UA association results for 880 SNPs in the SLC2A9 gene region.

S2 Fig. Serum UA association results for 879 SNPs in the SLC2A9 gene region, conditioned on rs11723439.

S3 Fig. SLC2A9-001 gene expression transcript interaction plot with dietary animal protein on urinary uric acid.

S1 Table. Correlations between gene-level gene expression for uricosuric genes in the Genetic Epidemiology Network of Arteriopathy study.

S2 Table. Forty SLC2A9 SNPs with the lowest p-values in linear mixed models association with urinary uric acid in the Genetic Epidemiology Network of Arteriopathy study.

S3 Table. Forty SLC2A9 SNPs with the lowest p-values in linear mixed models, conditioned on rs12509955, association with urinary uric acid in the Genetic Epidemiology Network of Arteriopathy study.

S4 Table. Forty SLC2A9 SNPs with the lowest p-values in linear mixed models association with serum uric acid in the Genetic Epidemiology Network of Arteriopathy study.

S5 Table. Forty SLC2A9 SNPs with the lowest p-values in linear mixed models, conditioned on rs11723439, association with serum uric acid in the Genetic Epidemiology Network of Arteriopathy study.

S6 Table. Forty SLC2A9 SNPs with the lowest p-values in linear mixed models association with SLC2A9 gene expression in the Genetic Epidemiology Network of Arteriopathy study.

S7 Table. Gene expression by dietary protein intake interaction associations for serum uric acid in the Genetic Epidemiology Network of Arteriopathy study.

S1 Text. Gene expression association with serum UA.

Author Contributions

References