Front. Genet. Frontiers in Genetics Front. Genet. 1664-8021 Frontiers Media S.A. 897051 10.3389/fgene.2022.897051 Genetics Original Research A Genetic Polymorphism in the WDR72 Gene is Associated With Calcium Nephrolithiasis in the Chinese Han Population Wang et al. WDR72 Polymorphism in Calcium Nephrolithiasis Wang Lujia 1 2 † Zhou Zijian 1 2 † Yang Yuanyuan 1 2 † Gao Peng 1 2 Lin Xiaoling 3 * Wu Zhong 1 2 * 1 Department of Urology, Huashan Hospital & Institute of Urology, Fudan University, Shanghai, China 2 Clinical Research Center of Urolithiasis, Shanghai Medical College, Fudan University, Shanghai, China 3 Department of Urology, Shanghai Children’s Hospital, School of Medicine, Shanghai Jiao Tong University, Shanghai, China

Edited by: Mikhail Churnosov, Belgorod National Research University, Russia

Reviewed by: Zhe Han, University of Maryland, Baltimore, United States

Kuanjun He, Inner Mongolia University for Nationalities, China

 *Correspondence: Xiaoling Lin, lxling.85@163.com; Zhong Wu, drzhongwu1964@126.com

These authors have contributed equally to this work

This article was submitted to Applied Genetic Epidemiology, a section of the journal Frontiers in Genetics

 14 07 2022 2022 13 897051 12 04 2022 24 06 2022 Copyright © 2022 Wang, Zhou, Yang, Gao, Lin and Wu. 2022 Wang, Zhou, Yang, Gao, Lin and Wu

This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

A previous genome-wide association study (GWAS) reported several novel loci for nephrolithiasis in British and Japanese population, some of which were predicted to influence CaSR signaling. In this study, we aimed to evaluate the association of these loci with calcium nephrolithiasis in Chinese Han population. We performed a case-control association analysis involving 691 patients with calcium nephrolithiasis and 1008 control subjects. We were able to genotype a total of 17 single-nucleotide polymorphisms (SNPs), which were previously reported to be significantly associated with nephrolithiasis in GWAS. rs578595 at WDR72 was significantly associated with calcium nephrolithiasis in Chinese Han population (p < 0.001, OR = 0.617). Moreover, rs12654812 at SLC34A1 (p = 0.0427, OR = 1.170), rs12539707 at HIBADH (p = 0.0179, OR = 0.734), rs1037271 at DGKH (p = 0.0096, OR = 0.828) and rs12626330 at CLDN14 (p = 0.0080, OR = 1.213) indicated suggestive associations with calcium nephrolithiasis. Our results elucidated the significance of genetic variation at WDR72, DGKH, CLDN14, SLC34A1, and HIBADH in Chinese patients with nephrolithiasis. Since polymorphisms of WDR72, DGKH, and CLDN14 are predicted to influence in CaSR signaling, our results emphasized the role of abnormal calcium homeostasis in calcium nephrolithiasis.

 nephrolithiasis WDR72 DGKH CLDN14 SLC34A1 HIBADH CaSR signaling Chinese Han population National Natural Science Foundation of China10.13039/501100001809 Introduction

Nephrolithiasis is one of the most frequent disorders affecting almost all populations. Epidemiological studies have reported that the incidence of nephrolithiasis is about 5% among females and 12% among males (Lewandowski and Rodgers, 2004), and almost 70% of all kidney stones are composed of calcium oxalate and/or phosphate (Croppi et al., 2012; Liu et al., 2021). In recent years, the incidence and prevalence of nephrolithiasis is increasing globally. Nephrolithiasis is commonly recurrent, with up to 50% of individuals experiencing a second episode within 10 years of the initial presentation, and recurrent stone disease has been linked to renal function decline (Pearle et al., 2014; Yan et al., 2021). The cause of nephrolithiasis is considered multifactorial, including but not limited to diet, ethnic, climate, and genetic factors. Recent studies estimated that up to 65% of kidney stone formers have a family history of nephrolithiasis (Stechman et al., 2007). Moreover, twin studies have reported a heritability of >45% for nephrolithiasis, and a strong family history of nephrolithiasis, including a parent and two siblings, have a standard incidence ratio for stone formation of >50 (Hemminki et al., 2018).

Up till now, six genome-wide association studies (GWAS) in different ethnicity have identified about 25 loci associated with nephrolithiasis (Thorleifsson et al., 2009; Gudbjartsson et al., 2010; Urabe et al., 2012; Oddsson et al., 2015; Howles et al., 2019; Tanikawa et al., 2019). Loci identified by GWAS in patients with nephrolithiasis could provide possible insight into the pathogenesis of the disorder. Howles et al. (2019) reported a trans-ethnic GWAS meta-analysis of British and Japanese cohorts, which identified 7 novel loci associated with nephrolithiasis. Among them, five of the loci, WDR72, GPIC1, DGKD, DGKH, and BCR, were predicted to influence calcium-sensing receptor (CaSR) signaling. The CaSR is a G protein-coupled receptor, which is highly expressed in the kidneys and parathyroid (Litvinova et al., 2021). CaSR has a central role in calcium hemostasis, including increasing kidney calcium reabsorption while stimulating parathyroid hormone (PTH) release to enhance bone resorption, urinary calcium reabsorption and 1,25-dihydroxyvitamin D3 synthesis in the kidney (Hannan et al., 2016). In this study, we conducted the research regarding the association between polymorphisms of those CaSR-related genes and calcium nephrolithiasis in Chinese Han population.

 Materials and Methods Subjects

In total, 691 unrelated Chinese Han patients with nephrolithiasis (467 males and 224 females, mean aged 50.47 years) were recruited at Huashan Hospital of Fudan University. Patients with nephrolithiasis secondary to known causes, such as chronic kidney disease, renal failure, chronic diarrhea, gout, renal tubular acidosis, primary and secondary hyperparathyroidism, osteoporosis, or cancer were excluded. Patients with radioparent stones, including struvite, uric acid and cystine stones were excluded. We also excluded patients with history of medications that affected urinary calcium excretion or taking vitamin D and/or calcium supplements. The control group consisted of 1008 subjects were age/gender matched individuals without a history of nephrolithiasis or a family history of kidney stone disease. Nephrolithiasis was diagnosed clinically either with plane radiography of kidney-ureter-bladder (KUB) or non-contrast computed tomography (CT) scan. All the patients with nephrolithiasis and the control subjects were of the same racial, ethnic, geographical and environmental strata.

We assessed the effect of genetic variations on serum calcium, sodium, potassium, magnesium, phosphorus, chloride, carbon dioxide (CO2), creatinine, urea, uric acid, alkaline phosphatase (ALP), parathyroid hormone (PTH), serum 25-hydroxycholecalciferol, albumin, glucose, cholesterol, triglycerides, low-density lipoprotein (LDL), high-density lipoprotein (HDL) levels; urine calcium and phosphorus levels; estimated glomerular filtration rate (eGFR); and body mass index (BMI). Laboratory measurements were performed within 2 weeks preoperatively.

 Blood Sample Collection and DNA Extraction

Peripheral blood samples were collected by venipuncture in a tube containing EDTA and was stored at −80°C. Genomic DNA was extracted using QiAamp DNA Blood Midi Kit (Qiagen, Germany). The concentration and quality of DNA was quantified by Qubit dsDNA HS Assay Kit (Promega, United States).

 SNP Selection and Genotyping

Genotyping was performed using the Illumina Asian Screening Array (ASA) BeadChip platform covering ∼660k variants across the genome. Imputation was performed with the IMPUTE computer program using the 1000 Genomes Project Han Chinese in Beijing (CHB) population as the reference, with imputation information score >0.90. We intended to replicate 20 single-nucleotide polymorphisms (SNPs) which were identified in the transethnic GWAS by Howles et al. (2019). We also evaluated the association of other 5 SNPs (rs1256328 at ALPL, rs7627468 at CASR, rs12654812 at SLC34A1, rs199565725 and rs219780 at CLDN 14) and calcium nephrolithiasis, which were identified in GWAS but not replicated in Chinese Han population (Thorleifsson et al., 2009; Oddsson et al., 2015). A standard quality control procedure was applied to select SNPs for further analysis. SNPs were excluded if they had: 1) genotype call rate <90% or 2) p < 0.001 for the Hardy-Weinberg Equilibrium (HWE) test. Therefore, eight SNPs were excluded because of low call rate of genotyping (<90%) and we were able to genotype a total of 17 SNPs. Population stratification analysis was performed by using an ancestry informative marker panel (UT-AIM250) (Wang et al., 2019).

 Statistical Analysis

Quantitative variables were presented as mean ± standard deviation (SD). An independent t test was used to compare the differences between the means of continuous variables. Categorical variables were analyzed using the Chi-square test. Genotype distributions for the SNP were tested for Hardy-Weinberg equilibrium (HWE). The association of SNPs with nephrolithiasis was tested by a Cochran-Armitage trend test. Results are expressed as odds ratio (OR) and 95% confidence intervals (CI). A p value lower than 2.94E-03 (0.05/17) was considered statistically significant. SNPs with p value less than 0.05 were also considered of interest. Multiple linear regression analyses were used to test association between genotype and clinical parameters, including serum calcium, phosphorus, creatinine, urea, uric acid, etc. with relevant covariates. We conducted association and QTL analyses using the PLINK-1.07 toolset. p-values were two tailed. An α of 0.05 was used to claim statistical significance.

 Results

The clinical characteristics of case and control samples were shown in Table 1. Compared with the healthy controls, body mass index (BMI) was significantly higher in patients with calcium nephrolithiasis (p = 0.003). There showed no significant difference in the distribution of serum calcium, phosphorus, magnesium, creatinine, and uric acid among the patients and controls.

Baseline characteristics of the study population.

Parameters Cases (n = 691) Controls (n = 1008) p Value
Gender (male/female) 467/224 705/303 0.302
Age (years) 50.47 ± 12.55 49.98 ± 12.65 0.431
BMI (kg/m2) 24.76 ± 3.30 24.28 ± 3.26 0.003
Stone frequency (primary/recurrence) 418/273 - -
Serum calcium (mmol/L) 2.26 ± 0.12 2.27 ± 0.10 0.062
Serum phosphorus (mmol/L) 1.13 ± 0.21 1.14 ± 0.18 0.293
Serum magnesium (mmol/L) 0.91 ± 0.93 0.89 ± 0.75 0.625
Serum creatinine (μ mol/L) 85.19 ± 41.93 83.64 ± 23.51 0.331
Serum uric acid (mmol/L) 0.360 ± 0.088 0.353 ± 0.077 0.083
Serum albumin (g/L) 42.11 ± 4.11 41.86 ± 4.09 0.217

The bold means p < 0.05.

Population stratification analysis indicated that data of both cases and controls overlapped with Asian populations (Supplementary Figure S1). Therefore, no population stratification was detected between cases and controls. We were able to genotype a total of 17 SNPs (Supplementary Table S1), which were previously reported to be significantly associated with nephrolithiasis in GWAS. Table 2 shows the genotype frequencies of polymorphism among all subjects. The genotype frequencies of 17 SNPs among case and control subjects were distributed in accordance with the Hardy-Weinberg equilibrium (Supplementary Table S2). SNP rs578595 at WDR72 was significantly associated with calcium nephrolithiasis in Chinese Han population (p < 0.001, OR = 0.617). Four SNPs at 4 loci—— rs12654812 at SLC34A1 (p = 0.0427, OR = 1.170), rs12539707 at HIBADH (p = 0.0179, OR = 0.734), rs1037271 at DGKH (p = 0.0096, OR = 0.828) and rs12626330 at CLDN14 (p = 0.0080, OR = 1.213) indicated suggestive associations with calcium nephrolithiasis.

Results of association analysis for calcium nephrolithiasis in Chinese Han population.

Chr SNP Gene Alleles Cases Controls p Value OR 95% CI
Minor Major n (11)a n (12) n (22) MAF n (11) n (12) n (22) MAF
1 rs10917002 ALPL T C 60 274 334 0.2949 82 429 487 0.2971 0.8924 0.990 (0.850–1.152)
1 rs1256328 ALPL T C 33 234 424 0.2171 40 320 647 0.1986 0.1912 1.119 (0.945–1.324)
2 rs780093 GCKR C T 148 335 208 0.4566 222 507 278 0.4722 0.3703 0.939 (0.818–1.077)
2 rs13003198 DGKD T C 123 269 214 0.4249 159 433 370 0.3903 0.0546 1.154 (0.997–1.336)
4 rs1481012 ABCG2 G A 79 299 313 0.3307 88 433 471 0.3070 0.1455 1.115 (0.962–1.292)
5 rs12654812 SLC34A1 A G 69 274 348 0.2981 82 373 553 0.2664 0.0427 1.170 (1.005–1.362)
6 rs77648599 SLC22A2 G T 1 12 648 0.0106 2 27 904 0.0166 0.1555 0.634 (0.335–1.196)
7 rs12539707 HIBADH T C 2 97 585 0.0738 3 166 709 0.0980 0.0179 0.734 (0.568–0.949)
7 rs12666466 AQP1 G C 7 108 576 0.0883 11 170 827 0.0952 0.4913 0.920 (0.725–1.167)
11 rs4529910 POU2AF G T 114 293 208 0.4236 184 490 330 0.4273 0.8357 0.985 (0.853–1.137)
13 rs1037271 DGKH T C 115 302 219 0.4182 210 486 279 0.4646 0.0096 0.828 (0.718–0.955)
15 rs578595 WDR72 A C 8 129 438 0.1261 33 310 649 0.1895 <0.001* 0.617 (0.502–0.759)
16 rs889299 SCNN1B A G 96 290 305 0.3488 134 451 397 0.3661 0.3038 0.927 (0.803–1.071)
17 rs4793434 SOX9 G C 77 312 302 0.3372 130 416 440 0.3428 0.7359 0.975 (0.844–1.128)
19 rs3760702 GIPC1 A G 31 191 392 0.2060 40 329 587 0.2139 0.5970 0.954 (0.800–1.137)
21 rs12626330 CLDN14 G C 155 310 168 0.4897 191 460 302 0.4418 0.0080 1.213 (1.052–1.399)
22 rs13054904 BCR A T 0 49 635 0.0358 2 72 930 0.0379 0.7591 0.944 (0.655–1.362)

n (11), number of subjects with homozygous genotypes for the minor allele; n (12), number of subjects with heterozygous genotypes; n (22), number of subjects with homozygous genotypes for the major allele.

Chr, chromosome; CI, confident interval; MAF, minor allele frequency; OR, odds ratio; SNP, single-nucleotide polymorphism. *p < 0.05/17.

The bold means p < 0.05.

As shown in Table 3, the A allele of rs12654812 was significantly correlated with lower level of serum glucose (p = 0.0384). The C allele of rs1037271 was significantly correlated with higher level of serum phosphorus (p = 0.0002). The C allele of rs578595 was significantly correlated with higher level of serum creatinine (p = 0.0018) and serum glucose (p = 0.0085). The G allele of rs12626330 was significantly correlated with higher level of serum uric acid (p = 0.0248) and serum carbon dioxide (p = 0.0265).

Multiple linear regression analyses for clinical parameters.

rs12654812 rs12539707 rs1037271 rs578595 rs12626330
Beta 1 s.e. 2 P Beta s.e P Beta s.e P Beta s.e P Beta s.e P
eGFR 3 −0.6906 1.9760 0.7271 −0.4097 3.4150 0.9046 1.7550 1.8140 0.3343 −1.1120 2.9420 0.7058 0.3609 1.7790 0.8394
Serum creatinine 1.1070 2.2720 0.6263 −0.8627 4.1430 0.8351 −0.6331 2.2460 0.7781 11.1500 3.5570 0.0018 4.0950 2.1630 0.0588
Serum urea 0.1155 0.1906 0.5447 −0.2427 0.3478 0.4856 −0.1985 0.1889 0.2938 0.2631 0.3056 0.3898 0.0426 0.1125 0.7052
Serum uric acid 0.0049 0.0053 0.3532 0.0078 0.0096 0.4147 0.0042 0.0051 0.4102 −0.0098 0.0080 0.2189 0.0112 0.0050 0.0248
Serum sodium −0.1843 0.1531 0.2291 −0.0265 0.2788 0.9242 0.1501 0.1484 0.3122 −0.0881 0.2395 0.7132 0.1985 0.1456 0.1732
Serum potassium −0.0737 0.0951 0.4385 −0.0698 0.1735 0.6876 0.1615 0.0948 0.0890 −0.0161 0.1597 0.9200 0.0107 0.0925 0.9083
Serum calcium −0.0007 0.0072 0.9235 0.0174 0.0131 0.1828 0.0059 0.0069 0.3923 −0.0124 0.0109 0.2558 0.0019 0.0065 0.7761
Serum magnesium −0.0494 0.0609 0.4176 −0.0585 0.1127 0.6037 0.0193 0.0606 0.7502 −0.0534 0.1014 0.5986 −0.0805 0.0585 0.1693
Serum phosphorus −0.0039 0.0126 0.7562 0.0035 0.0229 0.8771 0.0460 0.0121 0.0002 0.0150 0.0200 0.4532 −0.0065 0.0118 0.5780
Serum chloride 0.3476 0.7096 0.6246 0.4201 1.4030 0.7649 0.2567 0.2140 0.2315 0.5324 0.8509 0.5321 0.3760 0.7242 0.6040
Serum carbon dioxide −0.2227 0.1703 0.1914 0.4128 0.3098 0.1832 0.1315 0.1652 0.4263 −0.1527 0.2669 0.5675 −0.3602 0.1619 0.0265
Serum ALP −0.1350 1.4900 0.9278 −2.0310 2.7150 0.4548 −2.0600 1.4530 0.1568 −3.1180 2.3540 0.1858 1.0780 1.2440 0.3863
PTH −3.2670 5.7050 0.5684 9.4190 10.8400 0.3875 −4.7640 5.8240 0.4159 −0.9371 10.5200 0.9293 −2.7360 5.5800 0.6252
Serum 25(OH)D3 3.1600 3.5600 0.3783 3.3860 6.4910 0.6040 −2.7740 3.2960 0.4037 −3.4460 5.2270 0.5127 −3.0890 3.0750 0.3193
BMI 0.0511 0.2019 0.8004 0.0259 0.3666 0.9436 −0.1252 0.1958 0.5229 0.4291 0.3209 0.1817 0.1562 0.1932 0.4191
Serum albumin −0.1887 0.2502 0.4511 −0.1582 0.4658 0.7342 −0.1432 0.2453 0.5595 −0.2014 0.3700 0.5864 0.1475 0.2393 0.5379
Serum glucose −0.1579 0.0761 0.0384 0.1014 0.1399 0.4689 −0.0890 0.0753 0.2375 0.2855 0.1081 0.0085 0.0201 0.0718 0.7799
Serum total cholesterol −0.0577 0.2219 0.7959 0.6641 0.3903 0.0940 0.1290 0.2273 0.5727 0.3198 0.3467 0.3606 0.0323 0.2581 0.9010
Serum triglycerides 0.0092 0.1142 0.9360 0.1535 0.2057 0.4584 0.0233 0.1119 0.8355 −0.2589 0.1771 0.1498 0.0318 0.1343 0.8139
Serum LDL −0.4931 0.6252 0.4334 0.2574 1.4050 0.8552 0.3728 0.6377 0.5613 1.6360 0.9911 0.1049 −0.9971 0.7265 0.1759
Serum HDL −0.0260 0.0415 0.5337 0.0933 0.0909 0.3090 0.0072 0.0404 0.8601 −0.0701 0.0608 0.2551 −0.0303 0.0427 0.4811
Urine calcium 0.0824 0.6251 0.8958 0.4320 1.0930 0.6947 −0.3524 0.6987 0.6170 1.0610 1.4630 0.4739 −0.1075 0.6356 0.8666
Urine phosphorus 1.7620 1.8620 0.3497 −0.9060 3.2940 0.7847 3.6680 2.0420 0.0806 −4.7250 4.8090 0.3337 0.4739 2.0380 0.8174

Beta, Regression coefficient.

s.e., standard error of mean.

eGFR (ml/min/1.73m2) = 186 × (serum creatinine/88.41) - 1.154 × age-0.203 (×0.742 if female).

eGFR, estimated glomerular filtration rate; ALP, alkaline phosphatase; PTH, parathyroid hormone; 25(OH)D3, 25-hydroxycholecalciferol; BMI, body mass index; LDL, low-density lipoprotein; HDL, high-density lipoprotein.

The bold means p < 0.05.

 Discussion

In Howles study, which is the largest nephrolithiasis GWAS to date and integrates data from 12,123 stone formers and 417,378 controls from British and Japanese populations, identified 20 loci associated with nephrolithiasis, 7 of which have not previously been reported to associate with nephrolithiasis (Howles et al., 2019). In this study, we evaluated the association of 17 SNPs of 16 loci identified in former GWAS with nephrolithiasis in Chinese Han population. Our results indicated that rs578595 at WDR72, rs1037271 at DGKH, rs12654812 at SLC34A1, rs12539707 at HIBADH, and rs12626330 at CLDN14 were associated with the risk of calcium nephrolithiasis in Chinese Han population. Among them, WDR72 and DGKH were predicted to influence calcium-sensing receptor (CaSR) signaling.

WDR72 encodes WD repeat containing protein 72 (WDR72), an intracellular protein of 1102 amino acids with no known functional domains except a β-propeller structure composed of WD40 repeat domains in its N-terminus. This domain organization is the characteristic of vesicle coat proteins that mediate membrane deformation complexes to regulate intracellular vesicle trafficking (Katsura et al., 2014). Mutations in WDR72 have been previously identified as the cause of amelogenesis imperfecta (AI), a hereditary disease that affect tooth enamel formation (El-Sayed et al., 2009). Okada et al. (2012) reported that loss-of-function mutations of WDR72 result in AI, whereas missense mutations of WDR72 cause distal renal tubular acidosis (dRTA), possibly without the presence of AI. dRTA is characterized by an impairment of urinary acidification resulting in metabolic acidosis, hypokalemia, and inappropriately elevated urine pH. If not treated, this chronic condition eventually leads to nephrocalcinosis, nephrolithiasis and impaired renal function (D'Ambrosio et al., 2021). Recently, a trans-ethnic GWAS in British and Japanese populations identified that rs578595, an intronic variant in WDR72, was significantly associated with nephrolithiasis (Howles et al., 2019). WDR72 are thought to play a role in Clathrin-mediated endocytosis, a process central to sustained intracellular CaSR signaling (Wang et al., 2015). Our results demonstrated that rs578595 was significantly associated with nephrolithiasis in Chinese Han ethnicity. The risk allele of rs578595 was significantly correlated with higher level of serum creatinine, but not correlated with estimated-GFR or level of serum urea. A GWAS in Caucasian population indicated that WDR72 was associated with renal function and chronic kidney diseases (Kottgen et al., 2010). SNPs on WDR72 were also found to be associated with the estimated-GFR variance in American Indians (Franceschini et al., 2014). Patients with chronic kidney disease are more likely to develop nephrolithiasis. Moreover, impaired renal function and nephrolithiasis may both probably caused by a dRTA-related mechanism. It is not surprising that rs578595 is correlated with higher level of serum glucose, since WDR72 has been reported to be associated with higher HbA1c level and poorer blood glucose control (Paterson et al., 2010). Patients with diabetes mellitus are at increased risk for nephrolithiasis, since diabetes might cause stone formation by affecting the composition of urine. It is reported that patients with diabetes excrete more oxalate and have lower urine pH than non-diabetic people (Eisner et al., 2010). Howles et al. (2019) identified rs3760702, ∼300 bp upstream of GIPC1, as a significant locus associated with nephrolithiasis. GIPC1 encodes Regulator of G-protein signaling 19 interacting protein 1 (GIPC1), which is also postulated to play a role in clathrin-mediated endocytosis in CaSR signaling (Shang et al., 2017). In the current study, however, rs3760702 showed insignificant association with nephrolithiasis in Chinese Han population.

DGKH encodes for diacylglycerol kinase eta (DGKH), and DGKD encodes for diacylglycerol kinase delta (DGKD). DGKD and DGKH phosphorylates diacylglycerol, a component of the intracellular CaSR-signaling pathway inducing CaSR-mediated membrane ruffling and activating protein kinase C (PKC) signaling cascades including mitogen-activated protein kinase (MAPK) and intracellular calcium release (Schlam and Canton, 2017; Gorvin et al., 2018). SNPs of DGKH was predicted to promote kidney stone formation by influencing CaSR signaling. In a genome-wide association meta-analysis, DGKH was identified to be associated with serum calcium concentrations (O'Seaghdha et al., 2013). In our study, rs1037271, an intronic variant in DGKH, showed suggestive association with nephrolithiasis in Chinese Han ethnicity. The risk allele of rs1037271 was significantly correlated with higher level of serum phosphorus. Previously, DGKD was identified as a new loci associated with serum calcium in a genome-wide association meta-analysis in Europeans. They characterized the expression of in kidney, and demonstrated that both DGKD and DGKH were significantly upregulated in response to low calcium diet, which suggested specific involved of these genes in calcium homeostasis (O'Seaghdha et al., 2013). In a GWAS in British and Japanese populations, rs13003198, ∼6 kb upstream of DGKD, was identified as a significant locus for calcium nephrolithiasis. Moreover, they verified that DGKD knockdown could impair CaSR-signal transduction pathway in vitro, and this effect can be rectified with the calcimimetic cinacalcet (Howles et al., 2019). In the current study, SNP rs13003198 was not successfully replicated in Chinese Han population.

BCR encodes Breakpoint Cluster Region (BCR) protein, which is a GTPase-activating protein for RAC1 (Rac Family Small GTPase 1). RAC1 activation was postulated to mediate CaSR-induced membrane ruffling (Schlam and Canton, 2017). SNP rs13054904 located ∼110 kb upstream of BCR, which was identified to be associated with nephrolithiasis in British population but not Japanese population through GWAS (Howles et al., 2019). We failed to replicate rs13054904 in Chinese Han population, which suggested that this locus predisposed to nephrolithiasis in European rather than East Asian populations.

The CLDN14 gene encodes claudin-14, which belongs to the claudin family of membrane proteins. Claudin-14, a 239-amino acid protein with 4 transmembrane domains and intracellular N and C termini, is an important component of epithelial tight junctions (Tsukita and Furuse, 2000). In the kidney, claudin-14 is predominantly expressed in the thick ascending limb of the Henle’s loop (TALH) where a quarter of filtered calcium is reabsorbed through a passive paracellular pathway involving claudin-14, claudin-16 and claudin-19 (Olinger et al., 2018). CaSR activation is thought to increase expression levels of claudin-14 in the TALH and thereby decrease paracellular calcium reabsorption (Dimke et al., 2013). Moreover, genetic variants may attenuate claudin-14 activity and lead to enhanced paracellular divalent cation reabsorption in the TALH. The first GWAS on nephrolithiasis was reported in 2009, which identified CLDN14 as a significant locus for nephrolithiasis. Howles et al. (2019) demonstrated that rs12626330, an intronic variant in CLDN14, was associated with nephrolithiasis. In this study, rs12626330 showed suggestive association with nephrolithiasis, and rs12626330 was correlated with higher level of serum uric acid and lower level of serum CO2. It has been reported that serum CO2 level was negatively correlated with the risk of uric acid stone formation (Moreira et al., 2015). Our results suggested that the risk allele of rs12626330 might increase the risk of calcium nephrolithiasis through abnormal metabolism associated with hyperuricemia.

SLC34A1 gene encodes NPT2a, which is a member of the type II a sodium-phosphate co-transporter family. The NPT2a expressed in the brush border membrane of proximal tubular cells where the bulk of phosphate reabsorption takes place. Mutations in SLC34A1 have been reported to cause hypophosphatemic nephrolithiasis and osteoporosis in human (Prie et al., 2002). In knockout mice, severe renal phosphate wasting, hypercalciuria and skeletal abnormalities were observed (Beck et al., 1998). In 2012, a GWAS in a Japanese population identified SLC34A1 as a novel locus associated with nephrolithiasis (Urabe et al., 2012). In 2015, a GWAS in Icelanders reported that common variants of rs12654812 was associated with nephrolithiasis, and rs12654812 associated significantly with decreased serum PTH levels and serum phosphate (Oddsson et al., 2015). In GWAS of Howles et al. (2019) identified rs56235845 as a significant SNP associated with nephrolithiasis in both British and Japanese populations. In our results, we discovered a suggestive correlation between rs12654812, whereas rs56235845 was not significantly associated with nephrolithiasis. The risk allele of rs12654812 was associated with lower level of serum glucose. Presumably, reduction in serum PTH levels associated with rs12654812 may result from a decrease in serum phosphate levels caused by diminished renal reabsorption. Since patients with hyperparathyroidism are usually associated with impair glucose tolerance (Aojula et al., 2021), that kidney stone variance may negatively correlated with serum glucose level through decreased serum PTH.

HIBADH encodes 3-hydroxyisobutyrate dehydrogenase (HIBADH). HIBADH is an NAD+ -dependent mitochondrial enzyme that catalyzes oxidation of 3-hydroxyisobutyrate, an intermediate of valine catabolism, to methylmalonate semialdehyde. HIBADH is considered as a key enzyme in the gluconeogenesis pathway (Tasi et al., 2013). HIBADH gene was identified as a candidate gene for type 2 diabetes mellitus (Chen et al., 2013). However, its involvement in nephrolithiasis has not been fully elucidated. Howles et al. (2019) first identified rs12539707, an intronic variant in HIBADH, as a significant SNP associated with nephrolithiasis in British and Japanese populations. Our results showed suggestive association between rs12539707 and nephrolithiasis in Chinese Han population, which suggests that genes related to glucose metabolism might be involved in the mechanism of nephrolithiasis formation.

To conclude, the results of the present study elucidate that rs578595 at WDR72 is significantly associated with calcium nephrolithiasis, whereas rs1037271 at DGKH, rs12626330 at CLDN14, rs12654812 at SLC34A1 and rs12539707 at HIBADH show suggestive associations with nephrolithiasis in Chinese Han population. As mentioned above, the expression of CLDN14 localized to the TALH of the kidney was demonstrated to be regulated via the calcium-sensing receptor (CaSR) signaling. Moreover, WDR72 and DGKH are predicted to influence CaSR signaling, but it remains to be confirmed in the kidney. Although further investigation is required, we assumed that the polymorphism of WDR72, DGKH, and CLDN14 could increase the risk of calcium nephrolithiasis by influencing the CaSR signaling. Our results emphasized the role of abnormal calcium homeostasis in Chinese patients with calcium nephrolithiasis.

 Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.

 Ethics Statement

The studies involving human participants were reviewed and approved by the Huashan Institutional Review Board of Fudan University (HIRB). The patients/participants provided their written informed consent to participate in this study. Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.

 Author Contributions

Study design: WZ and WLJ. Data analysis and visualization: WLJ and LXL. Experiment validation and statistical analysis: ZZJ and YYY. Manuscript draft and revise: WLJ, WZ and GP.

 Funding

This project was supported by grants from the National Natural Science Foundation of China (No. 81970603 and 82100807) and the Scientific Research Foundation of Huashan Hospital Northern Branch (grant numbers: HSBY2019003).

 Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

 Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

 Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2022.897051/full#supplementary-material

References Aojula N. Khan S. Gittoes N. Hassan-Smith Z. (2021). Normocalcaemic Primary Hyperparathyroidism: What is the Role of Parathyroid Surgery? Ther. Adv. Endocrinol. 12, 204201882199537. 10.1177/2042018821995370 Beck L. Karaplis A. C. Amizuka N. Hewson A. S. Ozawa H. Tenenhouse H. S. (1998). Targeted Inactivation of Npt2 in Mice Leads to Severe Renal Phosphate Wasting, Hypercalciuria, and Skeletal Abnormalities. Proc. Natl. Acad. Sci. U.S.A. 95 (9), 5372–5377. 10.1073/pnas.95.9.5372 Chen J. Meng Y. Zhou J. Zhuo M. Ling F. Zhang Y. (2013). Identifying Candidate Genes for Type 2 Diabetes Mellitus and Obesity through Gene Expression Profiling in Multiple Tissues or Cells. J. Diabetes Res. 2013, 1–9. 10.1155/2013/970435 Croppi E. Ferraro P. M. Ferraro P. M. Taddei L. Gambaro G. (2012). Prevalence of Renal Stones in an Italian Urban Population: A General Practice-Based Study. Urol. Res. 40 (5), 517–522. 10.1007/s00240-012-0477-z D’Ambrosio V. Azzarà A. Sangiorgi E. Gurrieri F. Hess B. Gambaro G. (2021). Results of a Gene Panel Approach in a Cohort of Patients with Incomplete Distal Renal Tubular Acidosis and Nephrolithiasis. Kidney Blood Press Res. 46 (4), 469–474. 10.1159/000516389 Dimke H. Desai P. Borovac J. Lau A. Pan W. Alexander R. T. (2013). Activation of the Ca2+-Sensing Receptor Increases Renal Claudin-14 Expression and Urinary Ca2+ Excretion. Am. J. Physiol. Renal Physiol. 304 (6), F761–F769. 10.1152/ajprenal.00263.2012 Eisner B. H. Porten S. P. Bechis S. K. Stoller M. L. (2010). Diabetic Kidney Stone Formers Excrete More Oxalate and Have Lower Urine pH Than Nondiabetic Stone Formers. J. Urol. 183 (6), 2244–2248. 10.1016/j.juro.2010.02.007 El-Sayed W. Parry D. A. Shore R. C. Ahmed M. Jafri H. Rashid Y. (2009). Mutations in the Beta Propeller WDR72 Cause Autosomal-Recessive Hypomaturation Amelogenesis Imperfecta. Am. J. Hum. Genet. 85 (5), 699–705. 10.1016/j.ajhg.2009.09.014 Franceschini N. Haack K. Almasy L. Laston S. Lee E. T. Best L. G. (2014). Generalization of Associations of Kidney-Related Genetic Loci to American Indians. Clin. J. Am. Soc. Nephrol. 9 (1), 150–158. 10.2215/CJN.02300213 Gorvin C. M. Rogers A. Hastoy B. Tarasov A. I. Frost M. Sposini S. (2018). AP2σ Mutations Impair Calcium-Sensing Receptor Trafficking and Signaling, and Show an Endosomal Pathway to Spatially Direct G-Protein Selectivity. Cell Rep. 22 (4), 1054–1066. 10.1016/j.celrep.2017.12.089 Gudbjartsson D. F. Holm H. Indridason O. S. Thorleifsson G. Edvardsson V. Sulem P. (2010). Association of Variants at UMOD with Chronic Kidney Disease and Kidney Stones-Role of Age and Comorbid Diseases. PLoS Genet. 6 (7), e1001039. 10.1371/journal.pgen.1001039 Hannan F. M. Babinsky V. N. Thakker R. V. (2016). Disorders of the Calcium-Sensing Receptor and Partner Proteins: Insights into the Molecular Basis of Calcium Homeostasis. J. Mol. Endocrinol. 57 (3), R127–R142. 10.1530/JME-16-0124 Hemminki K. Hemminki O. Försti A. Sundquist K. Sundquist J. Li X. (2018). Familial Risks in Urolithiasis in the Population of Sweden. BJU Int. 121 (3), 479–485. 10.1111/bju.14096 Howles S. A. Wiberg A. Goldsworthy M. Bayliss A. L. Gluck A. K. Ng M. (2019). Genetic Variants of Calcium and Vitamin D Metabolism in Kidney Stone Disease. Nat. Commun. 10 (1), 5175. 10.1038/s41467-019-13145-x Katsura K. A. Horst J. A. Chandra D. Le T. Q. Nakano Y. Zhang Y. (2014). WDR72 Models of Structure and Function: A Stage-specific Regulator of Enamel Mineralization. Matrix Biol. 38, 48–58. 10.1016/j.matbio.2014.06.005 Köttgen A. Pattaro C. Böger C. A. Fuchsberger C. Olden M. Glazer N. L. (2010). New Loci Associated with Kidney Function and Chronic Kidney Disease. Nat. Genet. 42 (5), 376–384. 10.1038/ng.568 Lewandowski S. Rodgers A. L. (2004). Idiopathic Calcium Oxalate Urolithiasis: Risk Factors and Conservative Treatment. Clin. Chim. Acta 345 (1-2), 17–34. 10.1016/j.cccn.2004.03.009 Litvinova M. M. Khafizov K. Korchagin V. I. Speranskaya A. S. Asanov A. Y. Matsvay A. D. (2021). Association of CASR, CALCR, and ORAI1 Genes Polymorphisms with the Calcium Urolithiasis Development in Russian Population. Front. Genet. 12, 621049. 10.3389/fgene.2021.621049 Liu C.-J. Cheng C.-W. Tsai Y.-S. Huang H.-S. (2021). Crosstalk between Renal and Vascular Calcium Signaling: The Link between Nephrolithiasis and Vascular Calcification. Int. J. Mol. Sci. 22 (7), 3590. 10.3390/ijms22073590 Moreira D. M. Friedlander J. I. Carons A. Hartman C. Leavitt D. A. Smith A. D. (2015). Association of Serum Biochemical Metabolic Panel with Stone Composition. Int. J. Urol. 22 (2), 195–199. 10.1111/iju.12632 Oddsson A. Sulem P. Helgason H. Edvardsson V. O. Thorleifsson G. Sveinbjörnsson G. (2015). Common and Rare Variants Associated with Kidney Stones and Biochemical Traits. Nat. Commun. 6, 7975. 10.1038/ncomms8975 Okada Y. Sim X. Sim X. Go M. J. Wu J.-Y. Gu D. (2012). Meta-analysis Identifies Multiple Loci Associated with Kidney Function-Related Traits in East Asian Populations. Nat. Genet. 44 (8), 904–909. 10.1038/ng.2352 Olinger E. Houillier P. Devuyst O. (2018). Claudins: A Tale of Interactions in the Thick Ascending Limb. Kidney Int. 93 (3), 535–537. 10.1016/j.kint.2017.09.032 O'Seaghdha C. M. Wu H. Yang Q. Kapur K. Guessous I. Zuber A. M. (2013). Meta-analysis of Genome-wide Association Studies Identifies Six New Loci for Serum Calcium Concentrations. PLoS Genet. 9 (9), e1003796. 10.1371/journal.pgen.1003796 Paterson A. D. Waggott D. Boright A. P. Hosseini S. M. Shen E. Sylvestre M.-P. (2010). A Genome-wide Association Study Identifies a Novel Major Locus for Glycemic Control in Type 1 Diabetes, as Measured by Both A1C and Glucose. Diabetes 59 (2), 539–549. 10.2337/db09-0653 Pearle M. S. Goldfarb D. S. Assimos D. G. Curhan G. Denu-Ciocca C. J. Matlaga B. R. (2014). Medical Management of Kidney Stones: AUA Guideline. J. Urol. 192 (2), 316–324. 10.1016/j.juro.2014.05.006 Prié D. Huart V. Bakouh N. Planelles G. Dellis O. Gérard B. (2002). Nephrolithiasis and Osteoporosis Associated with Hypophosphatemia Caused by Mutations in the Type 2a Sodium-Phosphate Cotransporter. N. Engl. J. Med. 347 (13), 983–991. 10.1056/NEJMoa020028 Schlam D. Canton J. (2017). Every Day I'm Rufflin': Calcium Sensing and Actin Dynamics in the Growth Factor-independent Membrane Ruffling of Professional Phagocytes. Small GTPases 8 (2), 65–70. 10.1080/21541248.2016.1197873 Shang G. Brautigam C. A. Chen R. Lu D. Torres-Vázquez J. Zhang X. (2017). Structure Analyses Reveal a Regulated Oligomerization Mechanism of the PlexinD1/GIPC/myosin VI Complex. Elife 6, e27322. 10.7554/eLife.27322 Stechman M. J. Loh N. Y. Thakker R. V. (2007). Genetics of Hypercalciuric Nephrolithiasis: Renal Stone Disease. Ann. N. Y. Acad. Sci. 1116, 461–484. 10.1196/annals.1402.030 Tanikawa C. Kamatani Y. Terao C. Usami M. Takahashi A. Momozawa Y. (2019). Novel Risk Loci Identified in a Genome-wide Association Study of Urolithiasis in a Japanese Population. J. Am. Soc. Nephrol. 30 (5), 855–864. 10.1681/asn.2018090942 Tasi Y.-C. Chao H.-C. A. Chung C.-L. Liu X.-Y. Lin Y.-M. Liao P.-C. (2013). Characterization of 3-hydroxyisobutyrate Dehydrogenase, HIBADH, as a Sperm-Motility Marker. J. Assist. Reprod. Genet. 30 (4), 505–512. 10.1007/s10815-013-9954-8 Thorleifsson G. Holm H. Edvardsson V. Walters G. B. Styrkarsdottir U. Gudbjartsson D. F. (2009). Sequence Variants in the CLDN14 Gene Associate with Kidney Stones and Bone Mineral Density. Nat. Genet. 41 (8), 926–930. 10.1038/ng.404 Tsukita S. Furuse M. (2000). Pores in the Wall: Claudins Constitute Tight Junction Strands Containing Aqueous Pores. J. Cell Biol. 149 (1), 13–16. 10.1083/jcb.149.1.13 Urabe Y. Tanikawa C. Takahashi A. Okada Y. Morizono T. Tsunoda T. (2012). A Genome-wide Association Study of Nephrolithiasis in the Japanese Population Identifies Novel Susceptible Loci at 5q35.3, 7p14.3, and 13q14.1. PLoS Genet. 8 (3), e1002541. 10.1371/journal.pgen.1002541 Wang S. K. Hu Y. Yang J. Smith C. E. Nunez S. M. Richardson A. S. (2015). Critical Roles for WDR 72 in Calcium Transport and Matrix Protein Removal during Enamel Maturation. Mol. Genet. Genomic Med. 3 (4), 302–319. 10.1002/mgg3.143 Wang L.-J. Zhang C. W. Su S. C. Chen H.-I. H. Chiu Y.-C. Lai Z. (2019). An Ancestry Informative Marker Panel Design for Individual Ancestry Estimation of Hispanic Population Using Whole Exome Sequencing Data. BMC Genomics 20 (Suppl. 12), 1007. 10.1186/s12864-019-6333-6 Yan M.-T. Chao C.-T. Lin S.-H. (2021). Chronic Kidney Disease: Strategies to Retard Progression. Int. J. Mol. Sci. 22 (18), 10084. 10.3390/ijms221810084