- Research article
- Open Access
A single nucleotide polymorphism in the UMOD promoter is associated with end stage renal disease
BMC Medical Genetics volume 17, Article number: 95 (2016)
Several genome-wide association studies revealed that several variants of UMOD gene were related to the estimated glomerular filtration rate (eGFR), CKD or hypertension. In this study, we investigated the association between a common variant rs13333226 in the promoter region of UMOD gene and end stage renal disease (ESRD).
Variant rs13333226 of UMOD gene was genotyped by using the ABI Real time TaqMan allelic discrimination assay in a case-control study including 638 unrelated patients with ESRD and 366 controls.
The frequency of UMOD SNP rs13333226 GG/GA genotype was significantly higher (36.83% vs. 20.22%, P = 4.02 × 10-8) and the frequency of G allele was much higher (19.04% vs. 11.20%, P = 4.00 × 10−6) in the patients with ESRD than in the controls. The G allele was associated with an increased risk of ESRD (odds ratio 2.30, 95% confidence interval 1.70–3.11, P = 6.10 × 10−8). And G allele (odds ratio 2.33, 95% confidence interval 1.32–4.13, P = 3.65 × 10−3) was associated independently with ESRD.
A common variation rs13333226 in the promoter region of UMOD gene was independently associated with ESRD in Han Chinese.
End stage renal disease (ESRD), was one of the most serious principal challenges and resulted in a serious public health and financial burden worldwide. In addition, ESRD caused high morbidity and was an independent risk factor for cardiovascular disease or overall mortality in patients with chronic kidney disease (CKD) [1, 2]. It affected over 0.03% of the Chinese adult population .
Multiple factors were involved in the progression of kidney diseases into ESRD, and the development of ESRD was the result of interaction between genes and environment [4–7]. Familial clustering phenomenon was common in patients with ESRD [4, 5, 7, 8]. A recent study reported that a CKD patients, who had a first-degree relative with ESRD, significantly increased the relative risk of developing ESRD . Polymorphisms of many genes were reported to be associated with ESRD, such as ACE, APOL1, and MUC1 .
Recently, several genome-wide association studies (GWAS) showed that lots of genetic risk loci (such as: rs12917707, rs4293393, and rs6497476) were associated with estimated glomerular filtration rate (eGFR), CKD or hypertension [9–13]. Most studies consistently showed that the common variant rs12917707 of the UMOD gene have been strongly associated with eGFR, prevalent and incident CKD and ESRD in general population cohorts [9, 11, 14]. The relationship between T allele of rs12917707 and lower risk of CKD and ESRD was confirmed in other two case-control studies [15, 16]. These results suggested that the variants of UMOD gene contributed to the genetic susceptibility to ESRD. But a study recently reported that the association between rs12917707 and GFR was not observed in Italian diabetic patients . And the minor allele frequency of rs12917707 (T allele) was only 1% in Han Chinese (http://hapmap.ncbi.nlm.nih.gov/index.html.en). Another single nucleotide polymorphism (SNP) in the promoter region of UMOD gene, rs13333226, was reported to be associated with hypertension, and serum uric level [18, 19]. The G allele carriers had higher diastolic blood pressure and higher plasma uric acid compared with A/A homozygotes [18, 19]. In this study, we explored the association between rs13333226 and the susceptibility to ESRD.
Patients and clinical data
Total 638 unrelated Han Chinese patients with ESRD were recruited from southwest of China. The diagnostic criteria of ESRD were eGFR less than 15 or have received dialysis (hemodialysis or peritoneal dialysis). The patient who had acute kidney injury was excluded. The average age of their onset of ESRD (CKD 5 or 5D) was 54.1 ± 16.0 years old (from 14 to 87 years old). The primary causes of ESRD included chronic glomerulonephritis (CGN), diabetic kidney disease (DKD), hypertensive nephropathy (HTN), as well as other causes. Three hundred and sixty-six healthy individuals served as controls. All the controls didn’t have hypertension, diabetes, as well as acute or chronic kidney diseases. The average age of the controls was 52.9 ± 13.3 years-old (from 18 to 87 years old). The protocol was approved by the Ethics Review Committee of Sichuan Provincial People’s Hospital. Informed consent was signed by all subjects.
DNA extraction and genotyping
DNA was extracted from peripheral blood collected from patients with ESRD and healthy controls by using ‘salting-out’ method. The genotyping was conducted by real–time PCR. The primers and probe were obtained from Applied Biosystems (ABI, Foster city, CA). Genotyping of the rs13333226 SNP in the UMOD gene was performed using an ABI 7900 Real-time System (ABI, Foster city, CA). Assays were run at a final volume of 25 μl consisting of 12.5 μl of TaqMan®Gene Expression Master mix (ABI, Foster city, CA), 7.5 μl of primers (400 nM both forward and reverse) and probe (200 nM) (final concentration), 5 μl of input target DNA. Genotype clustering and calling were performed using SDS 2.3 Software (ABI, Foster city, CA).
Deviation from Hardy-Weinberg equilibrium was tested in healthy controls. Comparison of genotype and allele frequency between patient group and control group was done by X 2 test. Genotype-phenotype associations were tested under dominant as well as recessive genetic model. The effect of genotype on ESRD was also investigated with logistic regression analysis. Statistical analyses were performed using SPSS 17.0 software package (SPSS Inc., Chicago, IL). Statistical power was calculated by a software “PS: Power and Sample Size Calculation” . The SNP data of UMOD gene derived from 1000 Genomes project (GRCh38) CHB (Chinese Han Beijing) + CHS (Chinese Han South) were used to perform pairwise linkage disequilibrium analysis between rs13333226 and the previous GWAS hit rs12917707, rs4293393, as well as rs6497476 by Haploview software .
The percentage of male in the cases was much higher than in the controls (55.8% vs. 33.9%, P = 2.20 × 10−11). The demographic and clinical data were listed in Table 1. There was no deviation from Hardy-Weinberg equilibrium in controls (P = 0.073). The coefficient (r2) for pairwise linkage disequilibrium (LD) between rs13333226 and rs12917707, rs4293393, as well as rs6497476 were 0.099, 1, 0.883, respectively.
Compared to the healthy controls, the frequency of UMOD SNP rs13333226 GG/GA genotype was significantly higher (36.83% vs. 20.22%, P = 4.02 × 10-8) and the frequency of G allele was much higher in the patients with ESRD (19.04% vs. 11.20%, P = 4.00 × 10−6). The frequencies of GG/GA genotype and G allele were also much higher in ESRD patients with diabetes mellitus, without diabetes mellitus, or with hypertension (Table 2). For male individuals, the frequency of GG/GA genotype was much higher in patients than in controls (35.4% vs. 21.8%, P = 5.10 × 10−3), and for female individuals, the difference was significant (38.4% vs. 19.7%, P = 2.00 × 10−6).
The G allele was associated with higher risk of ESRD in dominant model [odds ratio (OR) = 2.30, 95% confidence interval (CI): 1.70–3.11, P = 6.10 × 10−8] with a statistical power of 1.00, in an additive model (OR = 1.96, 95%CI: 1.48–2.60,P = 2.00 × 10−6), but not in recessive model (OR = 0.57, 95% CI: 0.21–1.53, P = 0.26). The G allele was significantly associated with ESRD in dominant model in the subgroup of patients without diabetes mellitus (N = 538, OR = 2.21, 95% CI: 1.62–3.01, P = 5.55 × 10−7), with diabetes mellitus (N = 100, OR = 2.86, 95% CI: 1.78–4.58, P = 1.30 × 10−5), as well as with hypertension (N = 569, OR = 2.29, 95% CI: 1.69–3.11, P = 1.19 × 10−7).
The effect of genotype on ESRD was also investigated by a multivariable logistic regression analysis including age of initial ESRD, sex, hypertension and diabetic mellitus. The result showed that male (OR = 2.05, 95% CI: 1.19–3.51, P = 9.31 × 10−3) and G allele at rs13333226 (OR = 2.33, 95% CI: 1.32–4.13, P = 3.65 × 10−3) was independently associated with the susceptibility to ESRD.
UMOD gene encodes Tamm-Horsfall protein, also known as uromodulin, which is a glycoprotein exclusively synthesized by the thick ascending limb (TAL) and early distal convoluted tubule in the kidney . Lots of studies revealed that the variants of UMOD were associated with different risk of CKD, cardiovascular disease, hypertension, and hyperuricemia, T allele of rs12917707 was associated with lower risk of CKD and ESRD [9, 10, 18, 19, 23, 24]. But the relationship between variants of UMOD and risk of kidney disease was complicated, the several rare mutation in UMOD gene has been described as a cause of uromodulin-associated kidney disease, an autosomal dominant disease [25, 26]. Prudente et al. recently demonstrated that the previously reported strong association between rs12917707 and GFR in diabetic patients from Sweden was not observed in Italian diabetic patients . It suggested that there was a heterogeneous effect across the two different samples . In this study, we investigated the association between a promoter SNP rs13333226 of UMOD and susceptibility to ESRD.
Our results indicated that the frequency of rs13333226 GG/GA genotype and G allele was significantly higher in the patients with ESRD, compared to healthy controls. The G allele was associated with higher risk of ESRD in dominant model. G allele and male was associated independently with the susceptibility to ESRD by multivariate logistic analysis. This relationship still existed after adjusted by sex, age, as well as presence or absence of diabetes mellitus and hypertension. Our results demonstrated that rs13333226 didn’t correlate tightly with rs12917707 in Chinese Han population. It suggested that the effect of rs13333226 on CKD or ESRD was different from rs12917707. A previous case-control association study showed that G allele of rs13333226 was independently associated with a decreased risk of nephropathy, and the G allele was also associated with a higher eGFR and lower systolic blood pressure in the study cohort . However, two other studies have shown that G allele carriers have higher diastolic blood pressure and higher plasma uric acid compared with A/A homozygotes [18, 19].
The G allele at rs13333226 and C allele at rs6497476 were demonstrated to increase transcriptional activity in HEK-293 cells and mlMCD3 cells by luciferase reporter assay . The rs13333226 G allele carriers have higher plasma uric acid than A/A homozygotes . Higher serum uric acid was reported to associate significantly with increased risk for kidney functional impairment and SUA ≥6.0 mg/dL was a significant risk factor for rapid decline in eGFR . In a prospective children and adolescents cohort study, those with higher uric acid levels (≥5.5 mg/dL) had shorter times to renal outcomes by a multivariable parametric time-to-event analysis . Trudu et al. reported that two UMOD risk variants (G allele of rs12917707 and T allele of rs4293393) increased UMOD expression which led to salt-sensitive hypertension and to the presence of age-dependent renal lesions in transgenic mice . High expression of uromodulin could activate the renal sodium cotransporter NKCC2 which induced hypertension in UMOD transgenic mouse model . UMOD has also proved to regulate the renal outer medullary potassium channel, thus increasing salt reabsorption along the TAL . It could be speculated that the G allele at rs13333226 could increase carriers’ risk of ESRD by enhancing UMOD gene expression, increasing TAL salt reabsorption and serum uric acid levels, as well as increasing risk of hypertension. We haven’t detected serum sodium or uric acid levels in these patients, because the two parameters were influenced by many factors in patients with hemodialysis.
Although there was an independent association between G allele at rs13333226 and ESRD, we still need direct functional evidence to illustrate that whether the G allele at rs13333226 was a causal variant. And a prospective cohort study with large sample size would provide more solid evidence.
In summary, our data indicates that common genetic variation rs13333226 in the promoter region of UMOD gene was independently associated with ESRD in Han Chinese.
Chronic kidney disease
Estimated glomerular filtration rate
End stage renal disease
Single nucleotide polymorphism
Weiner DE, Tighiouart H, Amin MG, Stark PC, MacLeod B, Griffith JL, Salem DN, Levey AS, Sarnak MJ. Chronic kidney disease as a risk factor for cardiovascular disease and all-cause mortality: a pooled analysis of community-based studies. J Am Soc Nephrol. 2004;15(5):1307–15.
Sarnak MJ, Levey AS, Schoolwerth AC, Coresh J, Culleton B, Hamm LL, McCullough PA, Kasiske BL, Kelepouris E, Klag MJ, et al. Kidney disease as a risk factor for development of cardiovascular disease: a statement from the American Heart Association Councils on Kidney in Cardiovascular Disease, High Blood Pressure Research, Clinical Cardiology, and Epidemiology and Prevention. Hypertension. 2003;42(5):1050–65.
Zhang L, Wang F, Wang L, Wang W, Liu B, Liu J, Chen M, He Q, Liao Y, Yu X, et al. Prevalence of chronic kidney disease in China: a cross-sectional survey. Lancet. 2012;379(9818):815–22.
Spray BJ, Atassi NG, Tuttle AB, Freedman BI. Familial risk, age at onset, and cause of end-stage renal disease in white Americans. J Am Soc Nephrol. 1995;5(10):1806–10.
Akrawi DS, Li X, Sundquist J, Sundquist K, Zoller B. Familial risks of kidney failure in Sweden: a nationwide family study. PLoS One. 2014;9(11):e113353.
Satko SG, Freedman BI, Moossavi S. Genetic factors in end-stage renal disease. Kidney Int Suppl. 2005;94:S46–9.
Freedman BI, Robinson TW. Risk factors: familial aggregation of ESRD in Europeans-is it in the genes? Nat Rev Nephrol. 2014;10(12):677–8.
Skrunes R, Svarstad E, Reisaeter AV, Vikse BE. Familial clustering of ESRD in the Norwegian population. Clin J Am Soc Nephrol. 2014;9(10):1692–700.
Kottgen A, Glazer NL, Dehghan A, Hwang SJ, Katz R, Li M, Yang Q, Gudnason V, Launer LJ, Harris TB, et al. Multiple loci associated with indices of renal function and chronic kidney disease. Nat Genet. 2009;41(6):712–7.
Kottgen A, Pattaro C, Boger CA, Fuchsberger C, Olden M, Glazer NL, Parsa A, Gao X, Yang Q, Smith AV, et al. New loci associated with kidney function and chronic kidney disease. Nat Genet. 2010;42(5):376–84.
Liu CT, Garnaas MK, Tin A, Kottgen A, Franceschini N, Peralta CA, de Boer IH, Lu X, Atkinson E, Ding J, et al. Genetic association for renal traits among participants of African ancestry reveals new loci for renal function. PLoS Genet. 2011;7(9):e1002264.
Chambers JC, Zhang W, Lord GM, van der Harst P, Lawlor DA, Sehmi JS, Gale DP, Wass MN, Ahmadi KR, Bakker SJ, et al. Genetic loci influencing kidney function and chronic kidney disease. Nat Genet. 2010;42(5):373–5.
Regele F, Jelencsics K, Shiffman D, Pare G, McQueen MJ, Mann JF, Oberbauer R. Genome-wide studies to identify risk factors for kidney disease with a focus on patients with diabetes. Nephrol Dial Transplant. 2015;30 Suppl 4:iv26–34.
Padmanabhan S, Melander O, Johnson T, Di Blasio AM, Lee WK, Gentilini D, Hastie CE, Menni C, Monti MC, Delles C, et al. Genome-wide association study of blood pressure extremes identifies variant near UMOD associated with hypertension. PLoS Genet. 2010;6(10):e1001177.
Boger CA, Gorski M, Li M, Hoffmann MM, Huang C, Yang Q, Teumer A, Krane V, O’Seaghdha CM, Kutalik Z, et al. Association of eGFR-Related Loci Identified by GWAS with Incident CKD and ESRD. PLoS Genet. 2011;7(9):e1002292.
Reznichenko A, Boger CA, Snieder H, van den Born J, de Borst MH, Damman J, van Dijk MC, van Goor H, Hepkema BG, Hillebrands JL, et al. UMOD as a susceptibility gene for end-stage renal disease. BMC Med Genet. 2012;13:78.
Prudente S, Di Paola R, Copetti M, Lucchesi D, Lamacchia O, Pezzilli S, Mercuri L, Alberico F, Giusti L, Garofolo M, et al. The rs12917707 polymorphism at the UMOD locus and glomerular filtration rate in individuals with type 2 diabetes: evidence of heterogeneity across two different European populations. Nephrol Dial Transplant. 2016 [Epub ahead of print].
Han J, Chen Y, Liu Y, Liang Y, Wang X, Liu L, Wang F, Zhang L, Zhang H, Wang H. Common variants of the UMOD promoter associated with blood pressure in a community-based Chinese cohort. Hypertens Res. 2012;35(7):769–74.
Han J, Liu Y, Rao F, Nievergelt CM, O’Connor DT, Wang X, Liu L, Bu D, Liang Y, Wang F, et al. Common genetic variants of the human uromodulin gene regulate transcription and predict plasma uric acid levels. Kidney Int. 2013;83(4):733–40.
Dupont WD, Plummer Jr WD. Power and sample size calculations. A review and computer program. Control Clin Trials. 1990;11(2):116–28.
Barrett JC, Fry B, Maller J, Daly MJ. Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005;21(2):263–5.
Pennica D, Kohr WJ, Kuang WJ, Glaister D, Aggarwal BB, Chen EY, Goeddel DV. Identification of human uromodulin as the Tamm-Horsfall urinary glycoprotein. Science. 1987;236(4797):83–8.
Serafini-Cessi F, Malagolini N, Cavallone D. Tamm-Horsfall glycoprotein: biology and clinical relevance. Am J Kidney Dis. 2003;42(4):658–76.
Trudu M, Janas S, Lanzani C, Debaix H, Schaeffer C, Ikehata M, Citterio L, Demaretz S, Trevisani F, Ristagno G, et al. Common noncoding UMOD gene variants induce salt-sensitive hypertension and kidney damage by increasing uromodulin expression. Nat Med. 2013;19(12):1655–60.
Liu M, Chen Y, Liang Y, Liu Y, Wang S, Hou P, Zhang H, Zhao M. Novel UMOD mutations in familial juvenile hyperuricemic nephropathy lead to abnormal uromodulin intracellular trafficking. Gene. 2013;531(2):363–9.
Moskowitz JL, Piret SE, Lhotta K, Kitzler TM, Tashman AP, Velez E, Thakker RV, Kotanko P. Association between genotype and phenotype in uromodulin-associated kidney disease. Clin J Am Soc Nephrol. 2013;8(8):1349–57.
Ahluwalia TS, Lindholm E, Groop L, Melander O. Uromodulin gene variant is associated with type 2 diabetic nephropathy. J Hypertens. 2011;29(9):1731–4.
Toyama T, Furuichi K, Shimizu M, Hara A, Iwata Y, Sakai N, Perkovic V, Kobayashi M, Mano T, Kaneko S, et al. Relationship between serum uric acid levels and chronic kidney disease in a Japanese cohort with normal or mildly reduced kidney function. PLoS One. 2015;10(9):e0137449.
Rodenbach KE, Schneider MF, Furth SL, Moxey-Mims MM, Mitsnefes MM, Weaver DJ, Warady BA, Schwartz GJ. Hyperuricemia and progression of CKD in children and adolescents: The Chronic Kidney Disease in Children (CKiD) Cohort Study. Am J Kidney Dis. 2015;66:984–92.
Renigunta A, Renigunta V, Saritas T, Decher N, Mutig K, Waldegger S. Tamm-Horsfall glycoprotein interacts with renal outer medullary potassium channel ROMK2 and regulates its function. J Biol Chem. 2011;286(3):2224–35.
The authors wish to thank all the patients and control individuals for their participation and cooperation with this study.
This study was supported in part by National Basic Research Program of China 973 No. 2012CB517600 (No. 2012CB517604, to G.L.), National Natural Science Foundation of China (81170666, to G.L.) and Youth Science and Technology Creative Research Groups of Sichuan Province (2015TD0013, to G.L.).
Availability of data and materials
The authors declare that the data supporting the findings of this study are available within the article.
Conceived and designed the work: GL, TC and WL. Collection of the clinical data: TC, QW and GL. Performed the experiments: TC and QW. Analyzed the data: TC and GL. Wrote the paper: TC and GL. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
This study was approved by the Ethics Review Committee of Sichuan Provincial People’s Hospital. Written informed consent for analyses was obtained from all participants.
About this article
Cite this article
Chen, T., Wang, Q., Li, G. et al. A single nucleotide polymorphism in the UMOD promoter is associated with end stage renal disease. BMC Med Genet 17, 95 (2016). https://0-doi-org.brum.beds.ac.uk/10.1186/s12881-016-0358-3
- End stage renal disease
- Single-nucleotide polymorphism