Reasearch Awards nomination

Email updates

Keep up to date with the latest news and content from BMC Medical Genetics and BioMed Central.

Open Access Highly Accessed Research article

Effects of vitamin D receptor polymorphisms on urolithiasis risk: a meta-analysis

Pan Zhang1, Wei Nie2 and Hong Jiang1*

Author affiliations

1 Department of Nephrology, the First People’s Hospital of Jingzhou City, the First Hospital of Yangtze University, Jingzhou, Hubei Province 434000, China

2 Department of Respiratory Disease, Shanghai Changzheng Hospital, Second Military Medical University, Shanghai 200003, China

For all author emails, please log on.

Citation and License

BMC Medical Genetics 2013, 14:104  doi:10.1186/1471-2350-14-104

The electronic version of this article is the complete one and can be found online at:

Received:14 May 2013
Accepted:24 September 2013
Published:6 October 2013

© 2013 Zhang et al.; licensee BioMed Central Ltd.

This is an open access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.



Several studies analyzed the associations of Vitamin D receptor (VDR) polymorphisms with urolithiasis risk in different ethnic groups. However, the results were inconclusive. To evaluate a more precise estimation of the relationship, a meta-analysis was performed.


Pubmed, EMBASE, Wanfang Database, China National Knowledge Infrastructure (CNKI) and Weipu Database were searched. Data were extracted independently by two investigators. Odds ratios (ORs) with 95% confidence intervals (CIs) were used to assess the strength of associations.


Twenty-three case–control studies were included in this meta-analysis. Significant associations between ApaI, BsmI, FokI, and TaqI polymorphisms and urolithiasis risk were observed. However, sensitivity analyses for BsmI and FokI polymorphisms indicated that the results were not reliable and credible. In addition, there was a significant association of the ApaI-TaqI haplotype with urolithiasis risk.


This meta-analysis suggested that ApaI and TaqI polymorphisms in VDR gene were associated with urolithiasis risk.

Urolithiasis; Vitamin D receptor; Polymorphism; Meta-analysis


Urolithiasis is one of the most prevalent uronephrologic disorders and affects approximately 10% of individuals in western countries [1]. The incidence of urolithiasis is increasing. For example, in the US the prevalence has risen from 3.2% to 5.2% in just over two decades from the mid-1970s to the mid-1990s [2]. Previous studies evidenced the importance of genes in this disease. Studies of kidney stone-forming twins demonstrated a higher concordance for kidney stones in monozygotic than in dizygotic twins [3]. Additionally, a family history was reported to increase the disease risk (2.57 times higher) in males [4]. Thus, it is important to identify the gene variants contributing to urolithiasis pathogenesis.

Recently, Elkoushy and coworkers found that patients with urolithiasis had a high prevalence of inadequate vitamin D [5]. Expression and nuclear activation of the Vitamin D receptor (VDR) are necessary for the effects of vitamin D. Therefore, VDR was implicated in urolithiasis. In genetic hypercalciuric stone-forming (GHS) rats, Yao et al. [6] found that VDR mRNA levels were higher in kidney compared with wild-type controls. In addition, Favus et al. [7] showed that the level of VDR in peripheral blood monocytes was twofold greater in male calcium oxalate stone formers than in controls. Taken together, these results suggested that VDR may play an important role in the pathogenesis of urolithiasis.

The human VDR gene is located on chromosome 12q12-14. Four single nucleotide polymorphisms (SNPs) of the VDR gene have been widely studied [8]. ApaI, BsmI, and TaqI are located between the 8 and 9 exons in the 3′-untranslated region (UTR), and shown to be in strong linkage disequilibrium (LD) [9]. Another SNP is the FokI, which is located at the translation starting codon. Many studies investigated the associations between these polymorphisms of VDR gene and the risk of urolithiasis [10-32]. However, the results were inconclusive. The inconsistent results were possibly due to the low statistical powers of individual studies. The method of meta-analysis could provide a quantitative approach for combining the results of various studies with the same topic. Therefore, we performed this meta-analysis to address the precise relationship between the VDR gene variants and urolithiasis risk.


Publication search

We performed a systematic search of Pubmed, EMBASE, Wanfang Database, China National Knowledge Infrastructure (CNKI) and Weipu Database to find relevant studies. The search terms were used as follows: (urolithiasis or kidney stone or kidney stone disease) and (Vitamin D receptor or VDR) and (polymorphism or mutation or variant). Last search was updated in October, 2012. No language restriction was imposed. The reference lists of searched articles and relevant reviews were all perused to find additional eligible studies.

Inclusion and exclusion criteria

Two reviewers (Zhang and Nie) independently screened titles and abstracts of all studies for relevancy. Disagreements were resolved by discussion. Studies included in this meta-analysis based on the following selection criteria: (1) evaluation of the ApaI, BsmI, TaqI, and FokI polymorphisms in VDR gene and urolithiasis risk, (2) using a case–control design, and (3) genotype distributions in both cases and controls should be available for estimating an odds ratio (OR) and 95% confidence interval (CI).

Studies were excluded if one of the following existed: (1) the VDR polymorphisms were not analyzed or the outcome was not urolithiasis risk, (2) not case–control studies, such as the design based on family or sibling pairs, (3) not reported genotype frequencies or number, (4) abstracts or reviews, and (5) non-clinical study. For overlapping studies, the one with the largest sample size was included.

Data extraction

The full manuscripts of eligible studies were reviewed by two investigators (Zhang and Nie) independently. Any discrepancy was resolved by discussion or a third author (Jiang) would assess the articles. The following information was collected from each study: the first author’s name, year of publication, original country, ethnicity, age group, hypercalciuria in the urolithiasis group, composition of stone, sample size, the polymorphisms in VDR gene, genotyping method, and genotype number in cases and controls. We contacted the corresponding authors if more data was needed.

Quality assessment

The quality of the studies was assessed by two investigators (Zhang and Nie) independently. The predetermined criteria were modified from a previous review [33]. These scores were based on traditional epidemiological considerations and genetic issues. Scores ranged from zero (lowest) to ten (highest). Disagreement was settled by discussion. Articles scoring < 5 were defined as low quality, and those ≥ 5 were defined as high quality.

Statistical analysis

When the data from at least three similar studies were available, meta-analysis was performed. The strength of the associations between the ApaI, BsmI, FokI, and TaqI polymorphisms and urolithiasis risk was measured by ORs and 95% CIs. OR1, OR2, and OR3 were calculated for the genotypes: 1) AA vs. aa (OR1), aA vs. aa (OR2), and AA vs. aA (OR3) for the ApaI, 2) bb vs. BB (OR1), bB vs. BB (OR2), and bb vs. bB (OR3) for the BsmI, 3) ff vs. FF (OR1), fF vs. FF (OR2), and ff vs. fF (OR3) for the FokI, and 4) tt vs. TT (OR1), tT vs. TT (OR2), and tt vs. tT (OR3) for the TaqI, respectively. The statistical significance of OR was analyzed by Z test. These pairwise differences were used to indicate the most appropriate genetic model as follows: if OR1 = OR3 ≠ 1 and OR2 = 1, then a recessive model was suggested; if OR1 = OR2 ≠ 1 and OR3 = 1, then a dominant model was suggested; if OR2 = 1/OR3 ≠ 1 and OR1 = 1, then a complete overdominant model was suggested; if OR1 > OR2 > 1 and OR1 > OR3 > 1 (or OR1 < OR2 < 1 and OR1 < OR3 < 1), then a codominant model was suggested [34,35]. Once the best genetic model was identified, this model was used to collapse the three genotypes into two groups (except in the case of a codominant model) and to pool the results again. The pooled OR estimate of each study was calculated by the random-effects model.

The between-study heterogeneity was assessed by the Chi square-test based Cochrane Q-test and I2 test. I2 values of 25%, 50%, and 75% were nominally assigned as low, moderate, and high estimates. Departure from Hardy-Weinberg equilibrium (HWE) in controls was tested by the Chi-square test. Subgroup analyses were conducted by ethnicity, calciuria level, and age group. Sensitivity analyses were performed by excluding the studies not in HWE and the studies with low quality, respectively. Funnel plot was used to assess potential publication bias. Publication bias was also investigated statistically via Egger’s test [36].

All statistical tests were performed by using STATA 11.0 software (Stata Corporation, College Station, TX). A P value < 0.05 was considered statistically significant, except for test of heterogeneity where a level of 0.10 was used.


Study characteristics

Figure 1 outlines our study selection process. A total of 218 articles were identified after an initial search. Fifty-nine duplications were excluded. After reading the titles and abstracts, 127 articles were removed owing to abstracts, reviews, non-clinical studies, not case–control studies, and irrelevant to urolithiasis or VDR polymorphisms. After reading the full texts of the remaining 32 articles, 9 articles were excluded due to irrelevant to urolithiasis risk, no useful data, and reduplicate study. Finally, a total of 23 case–control studies met our inclusion criteria. There were 11 studies on ApaI, 8 on BsmI, 14 on FokI, 13 studies on TaqI. Four studies reported the haplotype of ApaI and TaqI polymorphisms. There were 14 studies of Asians and 9 studies of Caucasians. Seventeen studies were performed in adults, 3 in children, and 1 did not offer detailed information. Three studies only included patients with hypercalciuria, 6 studies included hypercalciuria patients partly and the data for these patients could be extracted, and 14 studies did not report detailed information. Quality scores for the each study ranged from 3 to 6. The characteristics of each study are presented in Table 1.

thumbnailFigure 1. Flow of study identification, inclusion, and exclusion.

Table 1. Characteristics of the case–control studies included in meta-analysis

Quantitative data synthesis

VDR ApaI polymorphism

Eleven studies determined the association between ApaI polymorphism and urolithiasis risk. Total sample sizes in urolithiasis and control groups were 1584 and 1853. The estimated OR1, OR2 and OR3 were 1.47, 1.30, and 1.04, respectively (Table 2). These estimates suggested a dominant genetic model, and therefore AA and aA were combined and compared with aa. As shown in Figure 2, the pooled OR was 1.34 (95% CI 1.11 – 1.60, P = 0.002). In the subgroup analysis by ethnicity, a significant association was found among Asians (OR = 1.43, 95% 1.16 – 1.75, P < 0.001) but not among Caucasians (OR = 1.08, 95% 0.74 – 1.57, P = 0.69). Subgroup analysis was also performed according to the calciuria level. However, no significant increased risk of urolithiasis was found among hypercalciuric patients (OR = 1.24, 95% CI 0.86 – 1.81, P = 0.25) (Table 2). In the subgroup analysis by age group, a significant association was observed among adults (OR = 1.30, 95% 1.04 – 1.62, P = 0.02). Sensitivity analysis was performed by excluding the studies that did not show HWE. The result was statistically significant (OR = 1.30, 95% CI 1.04 – 1.62, P = 0.02). Sensitivity analysis was also performed by excluding the low quality studies. The result was similar (OR = 1.34, 95% CI 1.09 – 1.66, P = 0.006). The shape of the funnel plot showed symmetric (Figure 3). Egger’s test did not indicate significant publication bias (P = 0.914).

Table 2. Summary of different comparative results

thumbnailFigure 2. Meta-analysis for the association between urolithiasis risk and the VDR ApaI polymorphism.

thumbnailFigure 3. Funnel plot for publication bias in selection of studies on the VDR ApaI polymorphism.

VDR BsmI polymorphism

Eight studies (1210 cases and 1160 controls) that identified the association between VDR BsmI polymorphism and urolithiasis risk were included in this meta-analysis. The estimated OR1, OR2 and OR3 were 1.79, 1.85, and 0.96, respectively (Table 2). These estimates suggested a dominant genetic model. The pooled OR was 1.81 (95% CI 1.03 – 3.17, P = 0.04) (Figure 4). In the subgroup analysis by ethnicity, a marginally significant association was found among Caucasians (OR = 2.43, 95% CI 1.02 – 5.80, P = 0.05) but not among Asians (OR 1.21, 95% CI 0.67 – 2.32, P = 0.49) (Table 2). In addition, subgroup analysis in hypercalciuric patients showed significant increased risk of urolithiasis (OR 2.43, 95% CI 1.36 – 4.35, P = 0.003). In the stratified analysis by ethnicity, no significant association was found among adults (OR 0.80, 95% CI 0.96 – 3.38, P = 0.07). However, sensitivity analyses conducted by excluding the studies not in HWE or low quality studies did not find the significant association between VDR BsmI polymorphism and urolithiasis risk (Table 2). No publication bias was detected by funnel plot (Figure 5) and Egger’s test (P = 0.461).

thumbnailFigure 4. Meta-analysis for the association between urolithiasis risk and the VDR BsmI polymorphism.

thumbnailFigure 5. Funnel plot for publication bias in selection of studies on the VDR BsmI polymorphism polymorphism.

VDR FokI polymorphism

For the VDR FokI polymorphism, fourteen studies including 2266 cases and 2418 controls were included in this meta-analysis. OR1, OR2, and OR3 were 1.59, 1.38, and 1.04, respectively. These estimates suggested a dominant genetic model. Therefore, the original grouping was collapsed, and ff and fF were combined, in accordance with a dominant model, into a f carrier group, the latter of which was compared with the FF genotype group. The pooled OR was 1.48 (95% CI 1.03 – 2.12, P = 0.03) (Figure 6). In the subgroup analysis by ethnicity, no significant association was found among Asians (OR = 1.32, 95% CI 0.94 – 1.86, P = 0.11) (Table 2). In the subgroup analysis by the calciuria level, there was still no significant association between VDR FokI polymorphism and urolithiasis risk in patients with hypercalciuria (OR = 1.15, 95% CI 0.81 – 1.63, P = 0.43) (Table 2). Statistically significant increased urolithiasis risk was observed among adults group (OR = 1.71, 95% CI 1.14 – 2.56, P = 0.009). Sensitivity analyses found that the significant result was altered when the low quality studies or studies with Hardy-Weinberg disequilibrium were omitted (Table 2). The shape of the funnel plot seemed symmetrical (Figure 7). Egger’s test did not show evidence of publication bias (P = 0.081).

thumbnailFigure 6. Meta-analysis for the association between urolithiasis risk and the VDR FokI polymorphism.

thumbnailFigure 7. Funnel plot for publication bias in selection of studies on the VDR FokI polymorphism.

VDR TaqI polymorphism

Thirteen studies including 1744 patients and 1944 controls addressed the association between VDR TaqI polymorphism and urolithiasis. The estimated OR1, OR2 and OR3 were 1.24, 1.33, and 0.92, respectively. These pooling estimates suggested a dominant genetic model. The pooled OR in the dominant genetic model was 1.30 (95% CI 1.08 – 1.55, P = 0.005) (Figure 8). In the stratified analysis, statistically significant increased urolithiasis risk was observed among Asians (OR = 1.39, 95% CI 1.11 – 1.73, P = 0.004) but not Caucasians (OR = 1.21, 95% CI 0.84 – 1.72, P = 0.31). Furthermore, patients with hypercalciuria had a significant increase of the risk of developing urolithiasis (OR = 1.56, 95% CI 1.06 – 2.31, P = 0.03). In addition, significantly increased urolithiasis risks were found among adults (OR = 1.34, 95% CI 1.05 – 1.71, P = 0.02). Sensitivity analysis did not change the result (Table 2). The funnel plot appeared to be symmetrical (Figure 9). Egger’s test were performed to estimate the publication bias. Publication bias was not observed (P = 0.122).

thumbnailFigure 8. Meta-analysis for the association between urolithiasis risk and the VDR TaqI polymorphism.

thumbnailFigure 9. Funnel plot for publication bias in selection of studies on the VDR TaqI polymorphism.

VDR ApaI-TaqI haplotype

In 4 out of 23 studies, the role of VDR ApaI-TaqI haplotype in urolithiasis risk was analyzed. The sample sizes in case group and control group were 830 and 1207, respectively. Carling et al. [37] reported that the At haplotype displayed higher levels of mRNA expression than aT haplotype. Thus, the effect of these two haplotypes on the risk of urolithiasis was evaluated in this meta-analysis. As shown in Figure 10, the pooling estimate showed a significant increased risk of developing urolithiasis in subjects with the At haplotype (OR = 1.35, 95% CI 1.09 – 1.68, P = 0.007). In the subgroup analysis by ethnicity, a significant association was found among Asians (OR = 1.53, 95% CI 1.17 – 2.00, P = 0.002). In the stratified analysis of age group, we also found a significant association (OR = 1.36, 95% CI 1.03 – 1.80, P = 0.03). Symmetrical funnel plot was obtained (Figure 11). Publication bias was not observed (P = 0.887).

thumbnailFigure 10. Meta-analysis for the association between urolithiasis risk and the VDR ApaI-TaqI haplotype.

thumbnailFigure 11. Funnel plot for publication bias in selection of studies on the VDR ApaI-TaqI haplotype.


This present meta-analysis systematic investigated the associations between VDR polymorphisms and urolithiasis risk. We found that ApaI polymorphism was a risk factor for developing urolithiasis. This result indicated that the carriers of the AA or Aa genotype had a 34% increased urolithiasis risk compared to those individuals with the aa genotype. Sensitivity analysis did not change this result, suggesting the solidity of our result. In the stratified analysis by ethnicity, the significant association was found in Asian population but not in Caucasian population. In addition, BsmI and FokI polymorphisms were also significantly associated with urolithiasis risk in the overall population. However, results from sensitivity analyses suggested that these results were not statistically robust. Therefore, these results should be interpreted with caution and more studies are needed to evaluate the effect of BsmI and FokI polymorphisms on urolithiasis risk. Specifically, we noted that hypercalciuric patients who harbored b allele may have increased urolithiasis risk. As for the TaqI polymorphism, a significant association was observed between this polymorphism and the risk of urolithiasis. Sensitivity analyses further strengthened the validity of this result. When subgroup analysis was conducted according to the ethnicity, the significant association was showed in Asians and lack of a significant association was detected in Caucasian population. In the hypercalciuria subgroup analysis, a 56% increased risk of urolithiasis was found in patients with hypercalciuria. Collectively, these results suggested that VDR BsmI and TaqI polymorphisms may play unique roles in the etiology of hypercalciuric urolithiasis. However, because the number of studies included in the hypercalciuria subgroup analyses was small, the results lacked sufficient reliability to confirm these associations in a definitive manner. Thus, future studies with larger sample sizes are needed to confirm our results. Furthermore, we assessed the association between the haplotypes of the ApaI and TaqI polymorphisms and urolithiasis risk. The result from this meta-analysis suggested that At haplotype may play a role in urolithiasis susceptibility, especially in Asians. Taken together, our results exhibited significant associations between ApaI and TaqI polymorphisms and urolithiasis risk in Asian population but not in Caucasian population. There are at least four reasons to explain the ethnic difference. First, the number of Caucasian patients was small. There were only 384 and 644 Caucasian patients for these two polymorphisms, respectively. It was therefore possible that the observed ethnic difference was due to chance. More studies with Caucasian population are required to validate the effect of ethnic differences. Second, higher heterogeneity was observed in the Caucasians subgroup (I2 = 40%) but not in the Asians subgroup (I2 = 0%) in TaqI polymorphism. This may distort the result. Third, urolithiasis is a complex disease. Both genetic and environmental factors could affect the risk of urolithiasis in different populations. It is possible that different urolithiasis risks in Asians and Caucasians were due to exposure to various environmental factors. To data, there was no reported study which was performed to assess the effect of VDR-environment interactions on urolithiasis risk in different ethnicities and regions. Therefore, epidemiologic studies should be designed to examine these associations in the future. Fourth, ethnic difference in the VDR gene allele frequencies may also result in this difference. For example, the A allele of the ApaI, occured with lower frequency in Caucasians when compared to Asians, while the t allele of the TaqI had a higher frequency in Asians compared to Caucasians [38]. We also performed subgroup analyses by age group. There were still significant associations in the adults group except in BsmI polymorphism. These results indicated that these polymorphisms were significantly associated with increased urolithiasis risk in adults. We did not evaluate the associations between VDR polymorphisms and urolithiasis risk in children due to insufficient data. More studies using children population are needed to determine the associations between VDR polymorphisms and urolithiasis risk.

A considerable weight of evidence supporting a role for VDR in urolithiasis was derived from GHS rats. Experimental studies using GHS rats showed increased levels of VDR in intestine, bones and kidneys [6,39]. Moreover, Karnauskas et al. [40] found that prolongation of VDR half-life increased VDR tissue levels and mediated VDR-regulated genes that led to hypercalciuria. Recently, using microRNA targeting VDR, Xi and colleagues silenced the VDR gene in the kidneys of GHS rats [41]. They demonstrated VDR knockdown in the kidney can upregulate the expression of transient receptor potential vanilloid receptor subtype 5 (TRPV5) in GHS rats [41]. Previous study found that mice lacking TRPV5 exhibited reduced calcium reabsorption, which caused severe hypercalciuria [42]. In addition, hypercalciuria increased the risk for calcium oxalate nephrolithiasis and occurred in up to 50% of nephrolithiasis patients [43,44]. Therefore, high VDR level might contribute to the development of urolithiasis. VDR was one of the most studied candidate genes for urolithiasis. Carling et al. [37] showed that the individuals exhibiting the BB, AA, or tt genotypes had significantly higher VDR levels than those with homozygous for the b, a, or T alleles. They also found those exhibiting the baT haplotype demonstrated a relative lower VDR mRNA level than subjects with non-baT haplotype [37]. Furthermore, Yamagata et al. [45] showed the VDR mRNA levels of allele t were significantly higher than those of allele T in peripheral blood mononuclear cell (PBMC). Thus, it is biologically plausible that subjects with the B, A, or t alleles may have increased risk of urolithiasis. Our findings supported this speculation and the At haplotype carriers had increased urolithiasis risk than the aT haplotype carriers. However, we found that b allele carriers, but not B allele carriers, were associated with urolithiasis, although this result was not robust. Notablely, hypercalciuric patients with b allele were seemed to have higher urolithiasis risk in the subgroup analysis. This result agreed with the report of Ferreira et al. [46], who also found that bb homozygous for VDR polymorphism was overrepresented in hypercalciuric stone formers. Additionally, Relan and co-authors reported that subjects with the bb genotype exhibited a higher urinary calcium excretion than the BB genotype [18]. Ruggiero et al. [10] found the similar phenomenon. The altered function of the VDR gene for the BsmI polymorphism may be associated with higher urinary calcium excretion. Hypercalciuria was a risk factor of urolithiasis. It was thus possible that b allele carriers may have higher urolithiasis risk than B allele carriers. As for the FokI polymorphism, Jurutka et al. [47] suggested that F variant possessed elevated transcriptional activity compared with the f variant. However, there was no difference between these two variants on the VDR expression [47]. Results from our study suggested that there was a significant association between f variant and urolithiasis. But the mechanism was unknown. The functional studies of VDR FokI polymorphism are still required. In addition, more studies are needed to confirm our results.

Our meta-analysis included a total of 23 case–control studies with 3046 cases and 3206 controls, while a previous meta-analysis by Lin et al. [48] only included 17 case–control studies with 2046 cases and 2303 controls. They found significant associations between FokI and TaqI polymorphisms and urolithiasis risk. However, we found ApaI and TaqI polymorphisms were significantly associated with urolithiasis risk. More case–control studies were included in our meta-analysis than Lin’s study. Therefore, our study may be more powerful and the conclusion might be more reliable. This reason might explain the difference between these two meta-analyses. In addition, our study had some advantages. First, Lin and coworkers tested multiple genetic models, including allelic comparison, a dominant model, and a recessive model. However, they did not correct for multiple comparisons or give biologic rationale for the choice of these genetic models. When the underlying genetic model is unknown, it is better to use pairwise comparisons of the three genotypes to find the best genetic model [34]. We used this method to avoide the problem of multiple comparisons. Second, it was the first time studying the VDR-hypercalciuria interactions on urolithiasis risk. Third, we explored the haplotype effect of ApaI and TaqI polymorphisms on the susceptibility to urolithiasis. Finally, the methodological issues for meta-analysis, such as, quality assessment was well investigated.

We should point out the importance of heterogeneity and publication bias, which might influence the results of meta-analysis. In our study, significant heterogeneity was found in the BsmI and FokI polymorphisms. Subgroup analysis was used to explore the sources of heterogeneity. After subgroup analysis by the calciuria level, the heterogeneity was effectively disappeared in hypercalciuric patients. Thus, it could be presumed that the relatively large heterogeneity mainly resulted from the calciuria level. In addition, funnel plot and Egger’s tests were used to find potential publication bias. No significant publication bias was detected.

Some limitations of our study should be addressed. First, the number of studies that were included in this meta-analysis was moderate. It was possible that some relevant published studies or unpublished studies with negative results were missed. Second, all of the studies were performed in Asians and Caucasians; thus, our results may be applicable only to these ethnic groups. Third, this study did not address gene-gene and gene-environment interactions, because insufficient information could be extracted from the original publications. Fourth, the data of the stone types and biochemical profiles was not showed in our meta-analysis because of insufficient information from the primary publications. Finally, we could not investigate the association between the haplotypes of the ApaI, BsmI, and TaqI polymorphisms and urolithiasis risk.


To our knowledge, this was the most comprehensive meta-analysis to assess the relationship between the VDR polymorphisms and urolithiasis susceptibility. Our results indicated that the ApaI and TaqI polymorphisms were associated with the risk of urolithiasis. Haplotype analysis suggested that ApaI-TaqI haplotypes conferred the susceptibility to urolithiasis. Future large-scale studies with more ethnic groups are needed to validate our findings. Gene-gene and gene-environment interactions should also be considered in future studies.

Competing interests

The authors have declared that no competing interests exist.

Authors’ contributions

PZ and WN collected the literature data, developed the statistical model, carried out the software implementation and drafted the manuscript. All authors read and approved the final manuscript.


  1. Rivers K, Shetty S, Menon M: When and how to evaluate a patient with nephrolithiasis.

    Urol Clin North Am 2000, 27(2):203-213. PubMed Abstract | Publisher Full Text OpenURL

  2. Stamatelou KK, Francis ME, Jones CA, Nyberg LM, Curhan GC: Time trends in reported prevalence of kidney stones in the United States: 1976–1994.

    Kidney Int 2003, 63(5):1817-1823. PubMed Abstract | Publisher Full Text OpenURL

  3. Goldfarb DS, Fischer ME, Keich Y, Goldberg J: A twin study of genetic and dietary influences on nephrolithiasis: a report from the Vietnam Era Twin (VET) registry.

    Kidney Int 2005, 67(3):1053-1061. PubMed Abstract | Publisher Full Text OpenURL

  4. Curhan GC, Willett WC, Rimm EB, Stampfer MJ: Family history and risk of kidney stones.

    J Am Soc Nephrol 1997, 8(10):1568-1573. PubMed Abstract | Publisher Full Text OpenURL

  5. Elkoushy MA, Sabbagh R, Unikowsky B, Andonian S: Prevalence and metabolic abnormalities of vitamin D-inadequate patients presenting with urolithiasis to a tertiary stone clinic.

    Urology 2012, 79(4):781-785. PubMed Abstract | Publisher Full Text OpenURL

  6. Yao J, Kathpalia P, Bushinsky DA, Favus MJ: Hyperresponsiveness of vitamin D receptor gene expression to 1, 25-dihydroxyvitamin D3: a new characteristic of genetic hypercalciuric stone-forming rats.

    J Clin Invest 1998, 101(10):2223-2232. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  7. Favus MJ, Karnauskas AJ, Parks JH, Coe FL: Peripheral blood monocyte vitamin D receptor levels are elevated in patients with idiopathic hypercalciuria.

    J Clin Endocrinol Metab 2004, 89(10):4937-4943. PubMed Abstract | Publisher Full Text OpenURL

  8. Vezzoli G, Terranegra A, Arcidiacono T, Soldati L: Genetics and calcium nephrolithiasis.

    Kidney Int 2010, 80(6):587-593. PubMed Abstract | Publisher Full Text OpenURL

  9. Morrison NA, Yeoman R, Kelly PJ, Eisman JA: Contribution of trans-acting factor alleles to normal physiological variability: vitamin D receptor gene polymorphism and circulating osteocalcin.

    Proc Natl Acad Sci U S A 1992, 89(15):6665-6669. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  10. Ruggiero M, Pacini S, Amato M, Aterini S, Chiarugi V: Association between vitamin D receptor gene polymorphism and nephrolithiasis.

    Miner Electrolyte Metab 1999, 25(3):185-190. PubMed Abstract | Publisher Full Text OpenURL

  11. Jackman SV, Kibel AS, Ovuworie CA, Moore RG, Kavoussi LR, Jarrett TW: Familial calcium stone disease: TaqI polymorphism and the vitamin D receptor.

    J Endourol 1999, 13(4):313-316. PubMed Abstract | Publisher Full Text OpenURL

  12. Chen WC, Chen HY, Hsu CD, Wu JY, Tsai FJ: No association of vitamin D receptor gene BsmI polymorphisms with calcium oxalate stone formation.

    Mol Urol 2001, 5(1):7-10. PubMed Abstract | Publisher Full Text OpenURL

  13. Chen WC, Chen HY, Lu HF, Hsu CD, Tsai FJ: Association of the vitamin D receptor gene start codon Fok I polymorphism with calcium oxalate stone disease.

    BJU Int 2001, 87(3):168-171. PubMed Abstract | Publisher Full Text OpenURL

  14. Nishijima S, Sugaya K, Naito A, Morozumi M, Hatano T, Ogawa Y: Association of vitamin D receptor gene polymorphism with urolithiasis.

    J Urol 2002, 167(5):2188-2191. PubMed Abstract | Publisher Full Text OpenURL

  15. Mossetti G, Vuotto P, Rendina D, Numis FG, Viceconti R, Giordano F, Cioffi M, Scopacasa F, Nunziata V: Association between vitamin D receptor gene polymorphisms and tubular citrate handling in calcium nephrolithiasis.

    J Intern Med 2003, 253(2):194-200. PubMed Abstract | Publisher Full Text OpenURL

  16. Ozkaya O, Söylemezoğlu O, Misirlioğlu M, Gönen S, Buyan N, Hasanoğlu E: Polymorphisms in the vitamin D receptor gene and the risk of calcium nephrolithiasis in children.

    Eur Urol 2003, 44(1):150-154. PubMed Abstract | Publisher Full Text OpenURL

  17. Wang S, Liu J, Hu S, Ye Z: Association of vitamin D receptor gene polymorphisms with calcium oxalate calculus disease.

    J Huazhong Univ Sci Technol 2003, 23(1):38-41. Publisher Full Text OpenURL

  18. Relan V, Khullar M, Singh SK, Sharma SK: Association of vitamin d receptor genotypes with calcium excretion in nephrolithiatic subjects in northern India.

    Urol Res 2004, 32(3):236-240. PubMed Abstract | Publisher Full Text OpenURL

  19. Rendina D, Mossetti G, Viceconti R, Sorrentino M, Castaldo R, Manno G, Guadagno V, Strazzullo P, Nunziata V: Association between vitamin D receptor gene polymorphisms and fasting idiopathic hypercalciuria in recurrent stone-forming patients.

    Urology 2004, 64(4):833-838. PubMed Abstract | Publisher Full Text OpenURL

  20. Hu SQ, Liu JH, Wang SG, Cao ZG, Wu W: Relationship between vitamin D receptor allele polymorphism and calcium oxalate stone disease.

    Chin J Urol 2004, 25(3):155-158. OpenURL

  21. Bid HK, Chaudhary H, Mittal RD: Association of vitamin-D and calcitonin receptor gene polymorphism in paediatric nephrolithiasis.

    Pediatr Nephrol 2005, 20(6):773-776. PubMed Abstract | Publisher Full Text OpenURL

  22. Bid HK, Kumar A, Kapoor R, Mittal RD: Association of vitamin D receptor-gene (FokI) polymorphism with calcium oxalate nephrolithiasis.

    J Endourol 2005, 19(1):111-115. PubMed Abstract | Publisher Full Text OpenURL

  23. Gunes S, Bilen CY, Kara N, Asci R, Bagci H, Yilmaz AF: Vitamin D receptor gene polymorphisms in patients with urolithiasis.

    Urol Res 2006, 34(1):47-52. PubMed Abstract | Publisher Full Text OpenURL

  24. Liu CC, Huang CH, Wu WJ, Huang SP, Chou YH, Li CC, Chai CY, Wu MT: Association of vitamin D receptor (Fok-I) polymorphism with the clinical presentation of calcium urolithiasis.

    BJU Int 2007, 99(6):1534-1538. PubMed Abstract | Publisher Full Text OpenURL

  25. Moyano MJ, de Tejada MJ G, García Lozano R, Moruno R, Ortega R, Martí V, Sánchez Palencia R, Miranda MJ, Palma A, Pérez Cano R: Alterations in bone mineral metabolism in patients with calcium kidney stone disease and polymorphism of vitamin D receptor. Preliminary results.

    Nefrologia 2007, 27(6):694-703. PubMed Abstract | Publisher Full Text OpenURL

  26. Seyhan S, Yavascaoglu I, Kilicarslan H, Dogan HS, Kordan Y: Association of vitamin D receptor gene Taq I polymorphism with recurrent urolithiasis in children.

    Int J Urol 2007, 14(12):1060-1062. PubMed Abstract | Publisher Full Text OpenURL

  27. Wang QZ, Qian B, Ding GF, Zheng LY: Vitamin D receptor gene polymorphisms in Chinese uygur patients with uroHthiasis in south Xinjian.

    J Pract Med 2009, 25(17):2805-2807. OpenURL

  28. Mittal RD, Mishra DK, Srivastava P, Manchanda P, Bid HK, Kapoor R: Polymorphisms in the vitamin D receptor and the androgen receptor gene associated with the risk of urolithiasis.

    Indian J Clin Biochem 2010, 25(2):119-126. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  29. Seo IY, Kang IH, Chae SC, Park SC, Lee YJ, Yang YS, Ryu SB, Rim JS: Vitamin D receptor gene Alw I, Fok I, Apa I, and Taq I polymorphisms in patients with urinary stone.

    Urology 2010, 75(4):923-927. PubMed Abstract | Publisher Full Text OpenURL

  30. Basiri A, Shakhssalim N, Houshmand M, Kashi AH, Azadvari M, Golestan B, Mohammadi Pargoo E, Pakmanesh H: Coding region analysis of vitamin D receptor gene and its association with active calcium stone disease.

    Urol Res 2012, 40(1):35-40. PubMed Abstract | Publisher Full Text OpenURL

  31. Wang S, Wang X, Wu J, Lin Y, Chen H, Zheng X, Zhou C, Xie L: Association of vitamin D receptor gene polymorphism and calcium urolithiasis in the Chinese Han population.

    Urol Res 2012, 40(4):277-284. PubMed Abstract | Publisher Full Text OpenURL

  32. Ruan L, Li ZM, Zheng RG, Huang WS, Shi GQ, Li G, Li S, Luo B: Relationship between vitamin D receptor FokI polymorphism and calcium oxalate stone disease in Guangzhou Chinese patients.

    Guangdong Med J 2012, 33(1):84-85. OpenURL

  33. Clark MF, Baudouin SV: A systematic review of the quality of genetic association studies in human sepsis.

    Intensive Care Med 2006, 32(11):1706-1712. PubMed Abstract | Publisher Full Text OpenURL

  34. Thakkinstian A, McElduff P, D’Este C, Duffy D, Attia J: A method for meta-analysis of molecular association studies.

    Statist Med 2004, 24(9):1291-1306. OpenURL

  35. Nie W, Chen J, Xiu Q: Cytotoxic T-lymphocyte associated antigen 4 polymorphisms and asthma risk: a meta-analysis.

    PloS ONE 2012, 7(7):e42062. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  36. Egger M, Smith GD, Schneider M, Minder C: Bias in meta-analysis detected by a simple, graphical test.

    BMJ 1997, 315(7109):629-634. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  37. Carling T, Rastad J, Åkerström G, Westin G: Vitamin D receptor (VDR) and parathyroid hormone messenger ribonucleic acid levels correspond to polymorphic VDR alleles in human parathyroid tumors.

    J Clin Endocrinol Metab 1998, 83(7):2255-2259. PubMed Abstract | Publisher Full Text OpenURL

  38. Uitterlinden AG, Fang Y, van Meurs JBJ, Pols HAP, van Leeuwen JPTM: Genetics and biology of vitamin D receptor polymorphisms.

    Gene 2004, 338(2):143-156. PubMed Abstract | Publisher Full Text OpenURL

  39. Bai S, Wang H, Shen J, Zhou R, Bushinsky DA, Favus MJ: Elevated vitamin D receptor levels in genetic hypercalciuric stone-forming rats are associated with downregulation of Snail.

    J Bone Miner Res 2010, 25(4):830-840. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  40. Karnauskas AJ, van Leeuwen JP, van den Bemd GJ, Kathpalia PP, DeLuca HF, Bushinsky DA, Favus MJ: Mechanism and function of high vitamin D receptor levels in genetic hypercalciuric stone‒forming rats.

    J Bone Miner Res 2005, 20(3):447-454. PubMed Abstract | Publisher Full Text OpenURL

  41. Xi QL, Wang SG, Ye ZQ, Zhu ZW, Li C, Bai J, Yu X, Liu JH: Effect of silencing VDR gene in kidney on renal epithelial calcium transporter proteins and urinary calcium excretion in genetic hypercalciuric stone-forming rats.

    Urology 2011, 78(6):1442.e1441-1442.e1447. OpenURL

  42. Hoenderop JG, van Leeuwen JP, van der Eerden BC, Kersten FF, van der Kemp AW, Mérillat AM, Waarsing JH, Rossier BC, Vallon V, Hummler E, Bindels RJ: Renal Ca2+ wasting, hyperabsorption, and reduced bone thickness in mice lacking TRPV5.

    J Clin Invest 2002, 112(12):1906-1914. OpenURL

  43. Moe OW: Kidney stones: pathophysiology and medical management.

    Lancet 2006, 367(9507):333-344. PubMed Abstract | Publisher Full Text OpenURL

  44. Coe FL, Evan A, Worcester E: Kidney stone disease.

    J Clin Invest 2005, 115(10):2598-2608. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  45. Yamagata M, Nakajima S, Tokita A, Sakai N, Yanagihara I, Yabuta K, Ozono K: Analysis of the stable levels of messenger RNA derived from different polymorphic alleles in the vitamin D receptor gene.

    J Bone Miner Res 1999, 17(3):164-170. Publisher Full Text OpenURL

  46. Ferreira LG, Pereira AC, Heilberg IP: Vitamin D receptor and calcium-sensing receptor gene polymorphisms in hypercalciuric stone-forming patients.

    Nephron Clin Pract 2010, 114(2):135-144. Publisher Full Text OpenURL

  47. Jurutka PW, Remus LS, Whitfield GK, Thompson PD, Hsieh JC, Zitzer H, Tavakkoli P, Galligan MA, Dang HT, Haussler CA, Haussler MR: The polymorphic N terminus in human vitamin D receptor isoforms influences transcriptional activity by modulating interaction with transcription factor IIB.

    Mol Endocrinol 2000, 14(3):401-420. PubMed Abstract | Publisher Full Text OpenURL

  48. Lin Y, Mao Q, Zheng X, Chen H, Yang K, Xie L: Vitamin D receptor genetic polymorphisms and the risk of urolithiasis: a meta-analysis.

    Urol Int 2011, 86(3):249-255. PubMed Abstract | Publisher Full Text OpenURL

Pre-publication history

The pre-publication history for this paper can be accessed here: