Open Access
Open Peer Review

This article has Open Peer Review reports available.

How does Open Peer Review work?

Genetic variation underlying renal uric acid excretion in Hispanic children: the Viva La Familia Study

  • Geetha Chittoor1,
  • Karin Haack2,
  • Nitesh R. Mehta3,
  • Sandra Laston4,
  • Shelley A. Cole2,
  • Anthony G. Comuzzie2,
  • Nancy F. Butte3 and
  • V. Saroja Voruganti1Email author
BMC Medical GeneticsBMC series – open, inclusive and trusted201718:6

DOI: 10.1186/s12881-016-0366-3

Received: 17 February 2016

Accepted: 30 December 2016

Published: 17 January 2017

Abstract

Background

Reduced renal excretion of uric acid plays a significant role in the development of hyperuricemia and gout in adults. Hyperuricemia has been associated with chronic kidney disease and cardiovascular disease in children and adults. There are limited genome-wide association studies associating genetic polymorphisms with renal urate excretion measures. Therefore, we investigated the genetic factors that influence the excretion of uric acid and related indices in 768 Hispanic children of the Viva La Familia Study.

Methods

We performed a genome-wide association analysis for 24-h urinary excretion measures such as urinary uric acid/urinary creatinine ratio, uric acid clearance, fractional excretion of uric acid, and glomerular load of uric acid in SOLAR, while accounting for non-independence among family members.

Results

All renal urate excretion measures were significantly heritable (p <2 × 10−6) and ranged from 0.41 to 0.74. Empirical threshold for genome-wide significance was set at p <1 × 10−7. We observed a strong association (p < 8 × 10−8) of uric acid clearance with a single nucleotide polymorphism (SNP) in zinc finger protein 446 (ZNF446) (rs2033711 (A/G), MAF: 0.30). The minor allele (G) was associated with increased uric acid clearance. Also, we found suggestive associations of uric acid clearance with SNPs in ZNF324, ZNF584, and ZNF132 (in a 72 kb region of 19q13; p <1 × 10−6, MAFs: 0.28–0.31).

Conclusion

For the first time, we showed the importance of 19q13 region in the regulation of renal urate excretion in Hispanic children. Our findings indicate differences in inherent genetic architecture and shared environmental risk factors between our cohort and other pediatric and adult populations.

Keywords

Uric acid clearance Genetic variation Fractional excretion of uric acid Hispanic children

Background

Renal excretion of uric acid is commonly involved in the development of gout, hyperuricemia, and nephropathy [1, 2]. The kidney filters freely circulating uric acid, accounting for ~70% of total uric acid excretion from the body [3]. The prevalence of hyperuricemia (increased serum uric acid (SUA) concentrations) and gout are on the rise along with other metabolic disorders such as obesity, type 2 diabetes, and metabolic syndrome [3, 4]. Hyperuricemia and hyperuricosuria (increased urinary uric acid (UrUA) concentrations) can lead to uric acid nephrolithiasis [57]. Moreover, these two are common multifactorial disorders that have been shown to have a familial inheritance and further associated with progression to chronic kidney disease [4, 8, 9]. Hyperuricemia is shown to cluster within families, heritabilities ranged from 39 to 45% in our family studies [10, 11], and twin studies up to 80% for serum uric acid [12], 60% for renal clearance of urate and 87% for fractional excretion of urate [13].

Uric acid is the end-product of purine metabolism, which is secreted, filtered, and reabsorbed with the help of specific urate transporters such as uric acid transporter-1 protein (URAT1) and solute carrier family 2, member 9 (SLC2A9) [4, 14, 15]. Any defect in these urate transporters, affecting the secretion and reabsorption processes, might lead to increased concentration in blood or reduced excretion of uric acid in the body [3, 9]. We have shown that variation in SUA concentrations is affected by genetic factors and is also associated with obesity and its comorbidities in Hispanic children [11]. However, there are no extensive genetic studies on parameters of renal urate excretion in children. Therefore, assessing the genetic contribution to excretion of uric acid and its related indices might help develop methods to better understand the renal urate excretion in children. In the present study, our aim was to identify genetic loci contributing to renal urate excretion measures using the genome-wide association (GWA) approach in Hispanic children of Viva La Familia Study (VFS).

Methods

Viva La Familia Study (VFS)

VFS was designed to identify genetic variants influencing pediatric obesity and its comorbidities in Hispanic children with the majority from Mexican American families. VFS study design, recruitment, and methodology, demographic and phenotypic information have been described in detail elsewhere [16]. All participants gave written informed consent or assent. Institutional Review Boards at Baylor College of Medicine and Affiliated Hospitals, Texas Biomedical Research Institute and UNC Chapel Hill approved the protocol for Human Subject Research. Methods used to measure fasting blood and 24-h urinary biochemistries are described elsewhere [11]. Availability of the phenotypic data after calculation of urinary indices was limited to 768 children.

Urinary and other indices calculation

Urine sample was collected over a 24-h period and were alkalized (pH around 8.0) using 0.5N NaOH during initial dilution. The urinary indices, using serum and urinary uric acid and creatinine, were calculated according to Perez-Ruiz, et al. [2] with following equations:
$$ \begin{array}{l}\mathrm{Body}\ \mathrm{S}\mathrm{urface}\ \mathrm{A}\mathrm{rea}\ \left(\mathrm{B}\mathrm{S}\mathrm{A},\ {\mathrm{m}}^2\right)\ \mathrm{calculated}\ \mathrm{b}\mathrm{y}\ \mathrm{Dubois}\ \mathrm{Equation} = 0.007184 \times \mathrm{Height}\ {\left(\mathrm{cm}\right)}^{0.725} \times \mathrm{Weight}\ {\left(\mathrm{kg}\right)}^{0.425}\hfill \\ {}\mathrm{Creatinine}\ \mathrm{clearance}\ \left(\mathrm{CrCl},\ \mathrm{ml}/ \min \right) = \mathrm{U}\mathrm{v} \times \left(\mathrm{UrCr}/\mathrm{SrCr}\right);\ \mathrm{U}\mathrm{v}:\ 24-\mathrm{hr}\ \mathrm{urine}\ \mathrm{v}\mathrm{o}\mathrm{lume}/\mathrm{time},\ \mathrm{U}\mathrm{rCr}:\ \mathrm{urinary}\ \mathrm{creatinine},\ \mathrm{S}\mathrm{rCr}:\ \mathrm{s}\mathrm{erum}\ \mathrm{creatinine}\hfill \\ {}\mathrm{U}\mathrm{ric}\ \mathrm{a}\mathrm{cid}\ \mathrm{clearance}\ \left(\mathrm{UACl},\ \mathrm{ml}/ \min \right) = \mathrm{U}\mathrm{v} \times \left(\mathrm{UrUA}/\mathrm{SrUA}\right);\ \mathrm{U}\mathrm{rUA}:\ \mathrm{urinary}\ \mathrm{uric}\ \mathrm{a}\mathrm{cid},\ \mathrm{S}\mathrm{rUA}:\ \mathrm{s}\mathrm{erum}\ \mathrm{uric}\ \mathrm{a}\mathrm{cid}\hfill \\ {}\mathrm{U}\mathrm{rinary}\ \mathrm{uric}\ \mathrm{a}\mathrm{cid}\ \mathrm{t}\mathrm{o}\ \mathrm{urinary}\ \mathrm{creatinine}\ \mathrm{ratio} = \mathrm{U}\mathrm{rUA}/\mathrm{UrCr}\ \mathrm{a}\mathrm{nd}\ \mathrm{expressed}\ \mathrm{a}\mathrm{s}\ \mathrm{a}\mathrm{b}\mathrm{s}\mathrm{o}\mathrm{lute}\ \mathrm{v}\mathrm{a}\mathrm{lue}\hfill \\ {}\mathrm{Fractional}\ \mathrm{excretion}\ \mathrm{o}\mathrm{f}\ \mathrm{uric}\ \mathrm{a}\mathrm{cid}\ \left(\mathrm{FEUA},\ \%\right) = \left(\mathrm{UrUA} \times \mathrm{S}\mathrm{rCr}\right)/\left(\mathrm{SrUA} \times \mathrm{U}\mathrm{rCr}\right) \times 100\hfill \\ {}\mathrm{Glomerular}\ \mathrm{load}\ \mathrm{o}\mathrm{f}\ \mathrm{uric}\ \mathrm{a}\mathrm{cid}\ \left(\mathrm{GLUA},\ \mathrm{mg}/ \min /{\mathrm{m}}^2\right) = \mathrm{CrCl} \times \mathrm{S}\mathrm{rUA}\hfill \\ {}\mathrm{Excretion}\ \mathrm{o}\mathrm{f}\ \mathrm{uric}\ \mathrm{a}\mathrm{cid}\ \mathrm{per}\ \mathrm{v}\mathrm{o}\mathrm{lume}\ \mathrm{o}\mathrm{f}\ \mathrm{glomerular}\ \mathrm{f}\mathrm{iltration}\ \left(\mathrm{EUAGF},\ \mathrm{mg}/\mathrm{dL}/{\mathrm{m}}^2\right) = \left(\mathrm{UrUA} \times \mathrm{S}\mathrm{rCr}\right)/\mathrm{UrCr}\hfill \end{array} $$

SNP genotyping

The Illumina HumanOmni1-Quad v1.0 BeadChip marker assays were used to genotype 1.1 million SNPs in 815 children enrolled in VFS [17]. Genotype calls were obtained after scanning on the Illumina BeadStation 500GX and analyzed by using the GenomeStudio software. Genotyping error rate was 2 per 100,000 genotypes (based on duplicates). The average call rate for all SNPs per individual sample was 97%. SNP genotypes were checked for Mendelian consistency using the program, SimWalk2 [18]. The estimates of the allele frequencies and their standard errors were obtained using SOLAR [19].

Heritability analysis

A variance components decomposition method was used to estimate heritability of uric acid and other renal urate excretion phenotypes. To estimate the genetic contribution to the variation in urinary indices, their heritability was estimated in SOLAR. Total phenotypic variance can be partitioned into its genetic and environmental components. The fraction of total phenotypic variance (VP) resulting from additive genetic effects (VG) is called heritability and is denoted as h2 = VG/VP [20]. All traits were adjusted for age, sex, their interaction effects, and body surface area.

Genome-wide association (GWA) study of urinary uric acid excretion measures

A total of 899,892 SNPs passed quality control and were included in the GWA analysis. GWAS was performed on 768 children from 260 Hispanic families including 1643 relative pairs. A measured genotype analysis (MGA) was performed on the inverse normal transformed residual traits (after regressing for covariate effects mentioned above) to minimize the non-normality distribution of the data using SOLAR. Each SNP genotype was converted in SOLAR to a covariate measure equal to 0, 1, or 2 copies of the minor alleles (or, the weighted covariate based on imputation for missing genotypes). These SNP covariates also were included in the variance components mixed models for MGA [21] vs. null models that incorporated the random effect of kinship and fixed effects such as age, sex, their interaction effects, and body surface area. For the initial GWA screen, we tested each SNP covariate independently as a 1-df likelihood ratio test. The linkage disequilibrium (LD) was computed in SOLAR by using information for all genotyped SNPs in all individuals. The effective number of SNPs accounting for LD was calculated by the method of Moskvina and Schmidt [22]. The average ratio of effective number of SNPs to the actual number obtained from analysis of non-overlapping bins of SNPs was used to calculate the genome-wide effective number of tests and thus the significance threshold for genome-wide association. Empirical thresholds for genome-wide significant and suggestive evidence of association were based on the distribution of p-values from 10,000 simulated null GWAS (i.e. simulations of a heritable trait with no modeled SNP covariate effects using the VFS pedigree and genotypes). The threshold for significance (p < 1 × 10−7) was defined as the cutoff for the lower 5% tails of the empirical distribution, and the threshold for suggestive evidence (p < 1 × 10−6) was the minimum p-value obtained not more than once per genome scan [11, 22].

Results

General characteristics of renal urate excretion measures

The study included data from 768 Hispanic children for the traits considered. General characteristics are given in Table 1 for both boys and girls and the total population, respectively. Mean age of boys and girls was approximately 12 years; body surface area (BSA) was higher (p <0.05) for boys (mean ± sd: 1.55 ± 0.40 m2 vs. 1.47 ± 0.33 m2) than girls, similar with BSA z-scores (0.33 ± 0.91 vs. 0.14 ± 0.76).
Table 1

General characteristics and renal urate excretion in Hispanic children

Variable

Boys (N = 397)

Girls (N = 371)

Total (N = 768)

Heritability estimates

Mean ± SD

Mean ± SD

Mean ± SD

h2 ± SE

P value

Age, years

11.88 ± 3.24

12.02 ± 3.48

11.95 ± 3.36

-

-

BSA, m2

1.55 ± 0.40*

1.47 ± 0.33*

1.51 ± 0.37

-

-

BSAZ

0.33 ± 0.91*

0.14 ± 0.76*

0.24 ± 0.85

-

-

UrCr, mg/dL

95.37 ± 39.53

83.86 ± 32.10

89.81 ± 36.57

0.41 ± 0.10

1.9 × 10−06

CrCl, ml/min

60.59 ± 35.11*

47.41 ± 26.20*

54.22 ± 31.80

0.52 ± 0.09

2.5 × 10−10

SUA, mg/dL

5.84 ± 1.76*

4.99 ± 1.43*

5.43 ± 1.66

0.63 ± 0.09

1.4 × 10−12

UrUA, mg/dL

15.63 ± 13.19*

13.43 ± 11.41*

14.57 ± 12.40

0.70 ± 0.05

1.1 × 10−24

UACl, ml/min

2.99 ± 2.31*

2.43 ± 1.96*

2.72 ± 2.17

0.65 ± 0.06

1.7 × 10−20

UrUA/UrCr

0.17 ± 0.11*

0.16 ± 0.11*

0.16 ± 0.11

0.73 ± 0.05

1.3 × 10−25

FEUA, %

6.02 ± 5.14

6.34 ± 6.02

6.17 ± 5.58

0.74 ± 0.05

6.7 × 10−27

GLUA, mg/min/m2

3.63 ± 2.64*

2.33 ± 1.46*

3.00 ± 2.25

0.57 ± 0.09

2.2 × 10−11

EUAGF, mg/dL/m2

0.33 ± 0.30

0.30 ± 0.31

0.32 ± 0.31

0.73 ± 0.05

6.1 × 10−25

*p-value is significantly different at <0.05 between boys and girls; BSA body surface area, BSAZ body surface area z-scores, UrCr urinary creatinine, SUA serum uric acid, UrUA urinary uric acid, UACl uric acid clearance, UrUA/UrCr urinary uric acid to urinary creatinine ratio, FEUA fractional excretion of uric acid, GLUA glomerular load of uric acid, EUAGF excretion of uric acid per volume of glomerular filtration. All traits for heritability analyses adjusted for age, sex, age*sex, age^2, age^2*sex, body surface area z-scores; p-value <0.05 considered statistically significant

Renal 24-h urate excretion measures and their descriptions are given in Table 1. Creatinine clearance (CrCl), serum uric acid (SUA), urinary uric acid (UrUA), uric acid clearance (UACl), urinary uric acid to urinary creatinine ratio (UrUA/UrCr), and glomerular load of uric acid (GLUA) were higher in boys than girls (p <0.05). However, urinary creatinine (UrCr), and excretion of uric acid per volume of glomerular filtration (EUAGF) tended to be higher in boys compared to girls; and, fractional excretion of uric acid (FEUA) tended to be higher in girls compared to boys, but did not reach statistical significance.

Heritability estimates

Table 1 also lists the heritability estimates and corresponding p-values for the phenotypes considered in this study. All phenotypes were adjusted for covariate effects (age, sex, age*sex, age2, age2*sex, body surface area z-scores). Significant heritability (h2) was detected for all the traits (p <2 × 10−6) and ranged from 0.41 to 0.74.

Genome-wide association analysis

Measured Genotype Analysis (MGA) was conducted using a variance components approach in 768 VFS Hispanic children that accounted for family kinships (Table 2). Genome-wide significant evidence of association was found for single nucleotide polymorphism (SNP) rs2033711, an intronic variant, in zinc finger protein 446 (ZNF446) on chromosome 19 with UACl (p <8 × 10−8, minor allele frequency (MAF) = 0.30) (Fig. 1). The effect size or the proportion of the residual phenotypic variance accounted for by the minor allele of the SNP was 4.5%. Genotype-specific means of this SNP showed that minor allele (G) was associated with increased UACl. In addition to ZNF446, a 72 kb region of chromosome 19q13 containing several zinc finger protein (ZNF) genes (ZNF324, ZNF584 and ZNF132) showed suggestive association with UACl (Fig. 2 & Additional file 1: Table S1), but after accounting for linkage disequilibrium (LD; r2 < 0.80) only ZNF446 and ZNF584 and an insulin receptor-related receptor (INSRR) gene on chromosome 1q21-23 showed evidence of suggestive association with UACl (Table 2).
Table 2

Results of genome-wide association analyses of renal urate excretion measures in Viva Hispanic children

SNP

Phenotype

Chr

Pos (GRCh38)

p (SNP)

β

SE

Effect size (%)

MAF

Gene symbol

rs2033711

UACl

19

58487983

7.9 × 10−08

0.32

0.06

4.5

G (0.30)

ZNF446

rs10908521

UACl

1

56843858

3. 0 × 10−07

−0.35

0.07

4.3

G (0.21)

INSRR

rs10423138

UACl

19

58416935

9. 4 × 10−07

0.29

0.06

3.7

G (0.31)

ZNF584

rs9874872

UrUA/UrCr

3

31467743

1.8 × 10−07

−0.38

0.07

4.7

A (0.18)

 

rs1353327

UrUA/UrCr

3

31461468

6. 8 × 10−07

−0.49

0.10

4.1

A (0.10)

 

rs17079585

CrCl

13

23921903

2. 5 × 10−07

−1.53

0.30

4.0

A (0.007)

ANKRD20A19P

rs4889855

FEUA

17

80539480

2. 7 × 10−07

−0.29

0.06

4.9

G (0.43)

 

rs17079585

GLUA

13

23921903

5. 7 × 10−07

−1.30

0.26

4.0

A (0.007)

ANKRD20A19P

SNP single nucleotide polymorphism, Chr Chromosome, Pos position in base pairs, p(SNP) p-value <1 × 10−7 shows evidence of significant association, p <1 × 10−6 shows evidence of suggestive association, β beta coefficient of the SNP, SE standard error of the beta coefficient; Effect Size: proportion of the residual phenotypic variance that is explained by the minor allele of the SNP; MAF minor allele frequency, UACl uric acid clearance, UrUA/UrCr urinary uric acid to urinary creatinine ratio, CrCl creatinine clearance, FEUA fractional excretion of uric acid, GLUA glomerular load of uric acid, ZNF446, zinc finger protein 446, INSRR Insulin receptor-related receptor, ZNF584 zinc finger protein 584, ANKRD20A19P ankyrin repeat domain 20 family, member A19, pseudogene

https://static-content.springer.com/image/art%3A10.1186%2Fs12881-016-0366-3/MediaObjects/12881_2016_366_Fig1_HTML.gif
Fig. 1

Genome-wide scan showing significant evidence of association for uric acid clearance with variants on chromosome 19

https://static-content.springer.com/image/art%3A10.1186%2Fs12881-016-0366-3/MediaObjects/12881_2016_366_Fig2_HTML.gif
Fig. 2

Locus Zoom plot showing the most significant SNPs on chromosome 19q13

A total of 7 SNPs were found to be suggestively (p <1 × 10−6) associated with one or the other renal urate excretion measures with MAFs ranging from 0.01 to 0.43, and the effect size ranged from 3.1 to 4.9%, respectively. After 19q13, another notable region where either GLUA or CrCl exhibited suggestive association was with a 13q12 SNP in ankyrin repeat domain 20 family, member A19, pseudogene (ANKRD20A19P). The genotype-specific mean values and the direction of associations (for minor alleles) of the SNPs in the genes associated with renal urate excretion measures from GWAS are shown in Table 3. For all UACl associated SNPs, the minor allele was associated with elevated UACl, except for rs10908521 on chromosome 1 (Additional file 2: Figure S1).
Table 3

Genotypic class specific mean values for top SNPs from genome-wide association analyses of renal urate excretion measures

SNP

Trait

Chr

Major/Major

Major/Minor

Minor/Minor

Direction of Changea

rs2033711

UACl

19

2.17 ± 1.79

2.79 ± 2.19

3.48 ± 2.69

rs10908521

UACl

1

2.86 ± 2.28

2.06 ± 1.66

1.86 ± 1.34

rs10423138

UACl

19

2.18 ± 1.78

2.77 ± 2.17

3.38 ± 2.68

rs9874872

UrUA/UrCr

3

0.18 ± 0.12

0.13 ± 0.10

0.11 ± 0.08

rs1353327

UrUA/UrCr

3

0.17 ± 0.11

0.11 ± 0.09

0.14 ± 0.13

rs17079585

CrCl

13

54.26 ± 32.25

22.47 ± 14.47

-

rs4889855

FEUA

17

7.22 ± 6.01

5.88 ± 5.74

4.34 ± 3.76

rs17079585

GLUA

13

3.02 ± 2.28

1.28 ± 0.77

-

aDirection of associations for minor alleles of the SNPs in the genes associated with renal handling measures from GWAS; Chr chromosome, UACl uric acid clearance, UrUA/UrCr urinary uric acid to urinary creatinine ratio, CrCl creatinine clearance, FEUA fractional excretion of uric acid, GLUA glomerular load of uric acid

Discussion

This GWAS identified a set of genes encoding zinc finger proteins on chromosome 19q13 that were associated with uric acid clearance in Hispanic children. To the best of our knowledge, this is the first family-based GWAS of renal urate excretion measures with prominent effect sizes reported in children.

The kidney is responsible for ~70% of the excretion of total body uric acid, and also reabsorbs about 90% of the filtered urate, regulating circulating uric acid levels. Epidemiological studies have shown that optimum uric acid levels may have antioxidant properties, however, elevated SUA and UrUA levels are associated with hypertension, inflammation, and also with kidney stones [3, 7, 23, 24]. A urate-transporting molecular complex (urate transportome) has been proposed as a model for urate transport in renal proximal tubules since several membrane proteins seem to be involved in urate transport [25]. According to this model, the mechanism of uric acid transport cannot be understood by evaluating one or two transporters. It has to be investigated as a functional unit comprising of all uric acid transporters and other molecules. This is because the renal handling of urate transport involves several genes (e.g., solute carrier family 2, member 9 (SLC2A9) and ATP-binding cassette ABC, subfamily G, member 2 (ABCG2), solute carrier family 16, member 9 (SLC16A9), solute carrier family 17, members 1, 3 and 4 (SLC17A1, SLC17A3 and SLC17A4), and, solute carrier family 22, members 11 and 12 (SLC22A11 and SLC22A12), most of which have been implicated in the regulation of urate levels [2629].

Defects in the tubular secretion of uric acid are mainly involved with urate clearance, hyperuricemia and gout [2]. Tubular secretion and reabsorption of urate levels change dynamically with age [9]. Genetic disorders have been indicated in the prevalence of hyperuricemia with increased urinary uric acid levels, which eventually lead to the formation of kidney stones and are often associated with chronic kidney disease [7, 30, 31]. Studies have demonstrated in adults that renal excretion of uric acid is under considerable genetic influence [32]. Our results indicated that all renal urate excretion measures are heritable and are within similar range of heritabilities reported in twin adults, [13] and Chinese twin children and adolescents [12].

Our GWAS findings on chromosome 19 revealed genes that have not been linked to renal urate handling before. We observed a strong genome-wide significant association of UACl with rs2033711 in ZNF446. Additionally, we found suggestive association of UACl with ZNF584 genetic variant (rs10423138), not in LD with rs2033711 (ZNF446) suggesting that these genetic variants are independently influencing variation in UACl. The ZNF family is one of the major human gene families and comprises of several of the currently recognized transcription factors. ZNFs through the zinc finger motifs are shown to interact with nucleic acids and are involved in various molecular mechanisms, cell differentiation and development [33, 34]. Studies indicate that ZNF446 is a novel member of Kruppel-related family, and believed to be one of the conserved proteins during human evolution. Also, ZNF446 gene is highly expressed in adult tissue cells (muscle), and may function as a transcriptional repressor in cellular growth and development [34, 35]. It is known to inhibit transcriptional activities of serum response element (SRE) activator protein 1(AP-1) [34]. Interestingly, AP-1 is a cis regulatory element regulating ABCG2 expression, one of the main uric acid transporter [36].

Moreover, strong association of UACl with ZNF446 itself is interesting because studies have shown association of ZNF365, on chromosome 10, with urolithiasis in children, [7] and with uric acid nephrolithiasis in adults [5, 37]. However, our results indicate a different pattern of zinc finger protein involvement in uric acid metabolism compared with adult studies [5, 7, 37]. The functional relevance of these genes with uric acid involvement is largely unknown, except that the emergence of ZNF365 variants is correlated with disappearance of uricase in primate evolution and hence causing a predisposition to hyperuricemia in humans [5, 37].

We found suggestive association of rs4889855 on chromosome 17 (Gene not identified) with FEUA. However, in adults, FEUA was shown to be associated with genetic variants of glucokinase regulator (GCKR on chromosome 2), SLC2A9 and ABCG2 on chromosome 4, and insulin like growth factor 1 receptor (IGF1R on chromosome 15) [38]. We also found suggestive evidence of association between GLUA and CrCl and ANKRD20A19P gene variant (on chromosome 13q12), both these associations prove to be interesting as beta-spectrin and ankyrin are key components of the cytoskeleton membrane that regulates clustering of sodium channel [39]. Additionally, we also found suggestive association of UACl with SNPs in INSRR on chromosome 1. Some studies have suggested functional association between mRNA expression of insulin receptor-related and insulin-like growth factor receptors and tumor cells [40]. One GWAS study indicated a new locus linked to chromosome 2p22.1-p21 for familial juvenile hyperuricaemic nephropathy, [41] but none of our associations with urinary uric acid handling have been reported before. The relatively small sample size in this study is a limitation. However, family-based studies have increased power to detect associations due to the fact that 768 children generate 1643 relative pairs and the degree of resemblance between relative pairs is considered for the genetic analysis.

Our study in VFS children is the first to report such genetic findings, and existing literature is limited in children involving GWAS on renal urate excretion measures. However, previous studies have reported association of 19q13 region with kidney phenotypes and diseases [42, 43]. This region has been linked to familial nephrotic syndrome [44], focal segmental glomerularosclerosis [45, 46] and cystinuria [47]. Our results are also different from studies that have reported association of uric acid transporters with renal uric acid excretion, for example, ABCG2 for renal urate excretion [27], and GCKR, SLC2A9, and IGF1R for FEUA [38], partly attributable to population substructure and sample size considered as majority of the GWAS are conducted in European or Asian descent populations. Nevertheless, our 19q13 region contains several ZNF proteins (including ZNF446) that may be involved in transcription of specific uric acid transporters. Thus, our GWAS associations could reflect differences in inherent genetic architecture and shared environmental risk factors between our VFS pediatric cohort and other study populations.

Conclusions

Our GWAS identified novel loci, particularly in 19q13, influencing the regulation of renal excretion of uric acid in Hispanic children. Our findings from children are not identical with those from adults [5, 37, 38] suggesting metabolic alterations in uric acid metabolism tracking from childhood to adults. The majority of the GWAS studies are conducted in adults highlighting the need for pediatric studies investigating the genetic underpinnings of the variation in uric acid earlier in the life course. It is essential to acquire knowledge on renal urate handling in children as it may reveal clinical and biological insights regarding the pathophysiology of uric acid excretion by the kidneys, given the inherited nature of these disorders.

Abbreviations

ABCG2

ATP-binding cassette ABC, subfamily G, member 2

CrCl: 

Creatinine clearance

EUAGF: 

Excretion of uric acid per volume of glomerular filtration

FEUA: 

Fractional excretion of uric acid

GCKR

Glucokinase regulator

GLUA: 

Glomerular load of uric acid

GWAS: 

Genome-wide association study

IGF1R

Insulin like growth factor 1 receptor

SLC2A9

Solute carrier family 2, member 9

SrCr: 

Serum creatinine

SUA: 

Serum uric acid

UACl: 

Uric acid clearance

URAT-1

Uric acid transporter-1 protein

UrCr: 

Urinary creatinine

UrUA: 

Urinary uric acid

UrUA/UrCr: 

Urinary uric acid to urinary creatinine ratio

VFS: 

Viva La Familia Study

ZNF: 

Zinc Finger Proteins

Declarations

Acknowledgements

We thank all the families who participated in the Viva La Familia Study. The authors wish to acknowledge the technical assistance of Grace-Ellen Meixner and Maria del Pilar Villegas. The contents of this publication do not necessarily reflect the views or policies of the USDA, nor does mention of trade names, commercial products, or organizations imply endorsement by the U.S.

Funding sources

The National Institutes of Health (NIH) [R01DK080457, R01MH59490, P30ES010126 and R01DK092238] and the USDA/ARS [Cooperative Agreement 6250-51000-053] supported this work.

Availability of data and materials

Data and materials are available in dbGaP.

Publicly available repositories.

CIP: Obesity-Diabetes Familial Risk, Viva La Familia Study.

dbGaP Study Accession: phs000616.v2.p2.

Authors’ contributions

Wrote paper: GC, NFB, VSV. Analyzed data or performed statistical analysis: GC, VSV. Provided essential reagents: KH, NRM, NFB, SAC. Primary responsibility for final content: VSV, GC, SAC, NFB, AGC. Read and provided edits to the paper: SL, KH, NRM, SAC, AGC. Designed research: VSV, GC, SAC, NFB, AGC. Conducted research: GC, NFB, VSV. All authors read and approved the final manuscript.

Competing interests

The authors declare that they have no competing interests.

Consent for publication

Not applicable.

Ethical approval and consent to participate

All participants gave written informed consent or assent. Institutional Review Boards at Baylor College of Medicine and Affiliated Hospitals, Texas Biomedical Research Institute and UNC Chapel Hill approved the protocol for Human Subject Research.

Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Authors’ Affiliations

(1)
Department of Nutrition and UNC Nutrition Research Institute, University of North Carolina at Chapel Hill
(2)
Department of Genetics, Texas Biomedical Research Institute
(3)
USDA/ARS Children’s Nutrition Research Center, Department of Pediatrics, Baylor College of Medicine
(4)
South Texas Diabetes and Obesity Institute, School of Medicine, University of Texas Rio Grande Valley

References

  1. Puig JG, Miranda ME, Mateos FA, Picazo ML, Jimanez ML, Calvin TS, et al. Hereditary nephropathy associated with hyperuricemia and gout. Arch Intern Med. 1993;153:357–65.View ArticlePubMedGoogle Scholar
  2. Perez-Ruiz F, Calabozo M, Erauskin GG, Ruibal A, Herrero-Beites AM. Renal underexcretion of uric acid is present in patients with apparent high urinary uric acid output. Arthritis Rheum. 2002;47:610–3.View ArticlePubMedGoogle Scholar
  3. Bobulescu IA, Moe OW. Renal transport of uric acid: evolving concepts and uncertainties. Adv Chronic Kidney Dis. 2012;19:358–71.View ArticlePubMedPubMed CentralGoogle Scholar
  4. Yee J. Uric acid: a clearer focus. Adv Chronic Kidney Dis. 2012;19:353–5.View ArticlePubMedGoogle Scholar
  5. Gianfrancesco F, Esposito T, Casu G, Maninchedda G, Roberto R, Pirastu M. Emergence of talanin protein associated with human uric acid nephrolithiasis in the hominidae lineage. Gene. 2004;339:131–8.View ArticlePubMedGoogle Scholar
  6. Gianfrancesco F, Esposito T. Multifactorial disorder: molecular and evolutionary insights of uric acid nephrolithiasis. Minerva Med. 2005;96:409–16.PubMedGoogle Scholar
  7. Medina-Escobedo M, Gonzalez-Herrera L, Villanueva-Jorge S, Martin-Soberanis G. Metabolic abnormalities and polymorphisms of the vitamin D receptor (VDR) and ZNF365 genes in children with urolithiasis. Urolithiasis. 2014;42:395–400.View ArticlePubMedGoogle Scholar
  8. Akl K, Ghawanmeh R. The clinical spectrum of idiopathic hyperuricosuria in children: isolated and associated with hypercalciuria/hyperoxaluria. Saudi J Kidney Dis Transpl. 2012;23:979–84.View ArticlePubMedGoogle Scholar
  9. Stiburkova B, Bleyer AJ. Changes in serum urate and urate excretion with age. Adv Chronic Kidney Dis. 2012;19:372–6.View ArticlePubMedGoogle Scholar
  10. Voruganti VS, Kent Jr JW, Debnath S, Cole SA, Haack K, Göring HH, et al. Genome-wide association analysis confirms and extends the association of SLC2A9 with serum uric acid levels to Mexican Americans. Front Genet. 2013;4:279.View ArticlePubMedPubMed CentralGoogle Scholar
  11. Voruganti VS, Laston S, Haack K, Mehta NR, Cole SA, Butte NF, et al. Serum uric acid concentrations and SLC2A9 genetic variation in Hispanic children: The Viva La Familia study. Am J Clin Nutr. 2015;101:1.View ArticleGoogle Scholar
  12. Ji F, Ning F, Duan H, Kaprio J, Zhang D, Zhang D, et al. Genetic and environmental influences on cardiovascular disease risk factors: a study of Chinese twin children and adolescents. Twin Res Hum Genet. 2014;17:72–9.View ArticlePubMedGoogle Scholar
  13. Emmerson BT, Nagel SL, Duffy DL, Martin NG. Genetic control of the renal clearance of urate: a study of twins. Ann Rheum Dis. 1992;51:375–7.View ArticlePubMedPubMed CentralGoogle Scholar
  14. 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:447–52.PubMedGoogle Scholar
  15. 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. J Biol Chem. 2004;279:16229–36.View ArticlePubMedGoogle Scholar
  16. Butte NF, Cai G, Cole SA, Comuzzie AG. Viva la familia study: genetic and environmental contributions to childhood obesity and its comorbidities in the Hispanic population. Am J Clin Nutr. 2006;84:646–54.PubMedGoogle Scholar
  17. Comuzzie AG, Cole SA, Laston SL, Voruganti VS, Haack K, Gibbs RA, et al. Novel genetic loci identified for the pathophysiology of childhood obesity in the Hispanic population. PLoS One. 2012;7:e51954.View ArticlePubMedPubMed CentralGoogle Scholar
  18. Sobel E, Lange K. Descent graphs in pedigree analysis: applications to haplotyping, location scores, and marker-sharing statistics. Am J Hum Genet. 1996;58:1323–37.PubMedPubMed CentralGoogle Scholar
  19. Almasy L, Blangero J. Multipoint quantitative-trait linkage analysis in general pedigrees. Am J Hum Genet. 1998;62:1198–211.View ArticlePubMedPubMed CentralGoogle Scholar
  20. Rogers J, Mahaney MC, Almasy L, Comuzzie AG, Blangero J. Quantitative trait linkage mapping in anthropology. Am J Phy Anthropol. 1999;Suppl29:127–51.View ArticleGoogle Scholar
  21. Boerwinkle E, Chakraborty R, Sing CF. The use of measured genotype information in the analysis of quantitative phenotypes in man. I. Models and analytical methods. Ann Hum Genet. 1986;50:181–94.View ArticlePubMedGoogle Scholar
  22. Moskvina V, Schmidt KM. On multiple-testing correction in genome-wide association studies. Genet Epidemiol. 2008;32:567–73.View ArticlePubMedGoogle Scholar
  23. Ames BN, Cathcart R, Schwiers E, Hochstein P. Uric acid provides an antioxidant defense in humans against oxidant- and radical-caused aging and cancer: a hypothesis. PNAS. 1981;78:6858–62.View ArticlePubMedPubMed CentralGoogle Scholar
  24. Feig DI, Soletsky B, Johnson RJ. Effect of allopurinol on blood pressure of adolescents with newly diagnosed essential hypertension: a randomized trial. JAMA. 2008;300:924–32.View ArticlePubMedPubMed CentralGoogle Scholar
  25. 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. Clin Exp Nephrol. 2012;16:89–95.View ArticlePubMedGoogle Scholar
  26. Dehghan A, Kottgen A, Yang Q, et al. Association of three genetic loci with uric acid concentration and risk of gout: a genome-wide association study. Lancet. 2008;372:1953–61.View ArticlePubMedPubMed CentralGoogle Scholar
  27. Matsuo H, Takada T, Ichida K, Nakamura T, Nakayama A, Ikebuchi Y, et al. Common defects of ABCG2, a high-capacity urate exporter, cause gout: a function-based genetic analysis in a Japanese population. Sci Transl Med. 2009;1:5ra11.View ArticlePubMedGoogle Scholar
  28. Merriman TR, Dalbeth N. The genetic basis of hyperuricaemia and gout. Joint Bone Spine. 2011;78:35–40.View ArticlePubMedGoogle Scholar
  29. 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. EJHG. 2014;22:938–43.View ArticlePubMedGoogle Scholar
  30. Baggio B. Genetic and dietary factors in idiopathic calcium nephrolithiasis. What do we have, what do we need? J Nephrol. 1999;12:371–4.PubMedGoogle Scholar
  31. Sebesta I. Genetic disorders resulting in hyper- or hypouricemia. Adv Chronic Kidney Dis. 2012;19:398–403.View ArticlePubMedGoogle Scholar
  32. Lipkowitz MS. Regulation of uric acid excretion by the kidney. Curr Rheumatol Rep. 2012;14:179–88.View ArticlePubMedGoogle Scholar
  33. Kadonaga JT, Carner KR, Masiarz FR, Tjian R. Isolation of cDNA encoding transcription factor Sp1 and functional analysis of the DNA binding domain. Cell. 1987;51:1079–90.View ArticlePubMedGoogle Scholar
  34. Liu F, Zhu C, Xiao J, Wang Y, Tang W, Yuan W, et al. A novel human KRAB-containing zinc-finger gene ZNF446 inhibits transcriptional activities of SRE and AP-1. Biochem Biophys Res Commun. 2005;333:5–13.View ArticlePubMedGoogle Scholar
  35. Xiao P, Chen Y, Jiang H, Liu YZ, Pan F, Yang TL, et al. In vivo genome-wide expression study on human circulating B cells suggests a novel ESR1 and MAPK3 network for postmenopausal osteoporosis. J Bone Miner Res. 2008;23:644–54.View ArticlePubMedPubMed CentralGoogle Scholar
  36. Stacy AE, Jansson PJ, Richardson DR. Molecular pharmacology of ABCG2 and its role in chemoresistance. Mol Pharmacol. 2013;84:655–69.View ArticlePubMedGoogle Scholar
  37. Gianfrancesco F, Esposito T, Ombra MN, Forabosco P, Maninchedda G, Fattorinin M, et al. Identification of a novel gene and a common variant associated with uric acid nephrolithiasis in a Sardinian genetic isolate. Am J Hum Genet. 2003;72:1479–91.View ArticlePubMedPubMed CentralGoogle Scholar
  38. Köttgen A, Albrecht E, Teumer A, Vitart V, Krumsiek J, Hundertmark C, et al. Genome-wide association analyses identify 18 new loci associated with serum urate concentrations. Nat Genet. 2013;45:145–54.View ArticlePubMedGoogle Scholar
  39. Komada M, Soriano P. [Beta]IV-spectrin regulates sodium channel clustering through ankyrin-G at axon initial segments and nodes of Ranvier. J Cell Biol. 2002;156:337–48.View ArticlePubMedPubMed CentralGoogle Scholar
  40. Elmlinger MW, Rauschnabel U, Koscielniak E, Haenze J, Ranke MB, Berthold A, et al. Correlation of type I insulin-like growth factor receptor (IGF-I-R) and insulin receptor-related receptor (IRR) messenger RNA levels in tumor cell lines from pediatric tumors of neuronal origin. Regul Pept. 1999;84:37–42.View ArticlePubMedGoogle Scholar
  41. Piret SE, Danoy P, Dahan K, Reed AA, Pryce K, Wong W, et al. Genome-wide study of familial juvenile hyperuricaemic (gouty) nephropathy (FJHN) indicates a new locus, FJHN3, linked to chromosome 2p22.1-p21. Hum Genet. 2011;129:51–8.View ArticlePubMedGoogle Scholar
  42. Chambers JC, Zhang W, Lord GM, van der Harst P, Lawlor DA, Sehmi JS, et al. Genetic loci influencing kidney function and chronic kidney disease in man. Nat Genet. 2010;42:373–5.View ArticlePubMedPubMed CentralGoogle Scholar
  43. Satko SG, Freedman BI. The importance of family history on the development of renal disease. Curr Opin Nephrol Hypertens. 2004;13:337–41.View ArticlePubMedGoogle Scholar
  44. Vats A, Nayak A, Ellis D, Randhawa PS, Finegold DN, Levinson KL, et al. Familial nephrotic syndrome: clinical spectrum and linkage to chromosome 19q13. Kidney Int. 2000;57:875–81.View ArticlePubMedGoogle Scholar
  45. Mathis BJ, Kim SH, Calabrese K, Haas M, Seidman JG, Seidman CE, et al. A locus for inherited focal segmental glomerulosclerosis maps to chromosome 19q13. Kidney Int. 1998;53:282–6.View ArticlePubMedGoogle Scholar
  46. Winn MP, Conlon PJ, Lynn KL, Howell DN, Gross DA, Rogala AR, et al. Clinical and genetic heterogeneity in familial focal segmental glomerulsclerosis. Kidney Int. 1999;55:1241–6.View ArticlePubMedGoogle Scholar
  47. Langen H, von Kietzell D, Byrd D, Arslan-Kirchner M, Vester U, Stuhrmann M, et al. Renal polyamine excretion, tubular amino acid reabsorption and molecular genetics in cystinuria. Pediatr Nephrol. 2000;14:376–84.View ArticlePubMedGoogle Scholar

Copyright

© The Author(s). 2017