- Research article
- Open Access
- Open Peer Review
Association study between vitamin D receptor gene polymorphisms and asthma in the chinese han population: a case-control study
BMC Medical Genetics volume 10, Article number: 71 (2009)
Modulation of the immune system is one of the principal roles of Vitamin D, for which the effects are exerted via the vitamin D receptor (VDR). Importantly, variants in the VDR gene have been susceptible in the past to raise the risk of asthma in several populations. These effects of VDR allelic markers remain speculative in the Chinese Han population.
A case-control study of 1090 individuals including 567 asthmatic patients was realized on five SNPs within the VDR gene. Only rs7975232 (ApaI) marker showed a significant association with asthma (P = 0.009). Haplotype analysis of the five VDR polymorphisms showed a significant association with asthma (global-p value = 0.012).
Although the susceptibility of VDR gene variants with asthma could not be confirmed for all SNPs tested in this study, the significant association obtained for rs7975232 provides evidence for a previously unknown report about the Chinese Han population and may raise the susceptibility of VDR to be a candidate gene for asthma.
Recent discoveries made in the field of molecular biology have reported that vitamin D derived from diet and supplements is not a vitamin, but a steroid with immunosuppressive properties when elevated . It has also been linked with lessening symptoms of autoimmune diseases like multiple sclerosis , Crohn's disease , and rheumatoid arthritis .
However, the effects of this steroid-like molecule are essentially exerted via its cognate the vitamin D receptor (VDR). In addition, molecular modeling findings have defined an interaction between the hormonal 1, 25-dihydroxyvitamin D3 form and its receptor through a ligand-binding domain (LBD) . The activated VDR plays then a fundamental role in the body by regulating numerous primary target genes. Three years before, a Canadian study has identified on a large scale that VDR may transcribe more than 900 genes . In fact, the interaction of 1, 25(OH)2D3 with VDR and its cofactors modulates many biological activities of the neural, immune, and endocrine systems; including calcium and phosphorous homeostasis, apoptosis and cell differentiation. For instance, a previous study has demonstrated that VDR rs7975232 polymorphism is associated with hip bone mineral density (BMD) in a group of 260 healthy postmenopausal Chinese Han women .
Besides to this classical action, and due to the pleiotropic effect the 1, 25(OH)2D3-VDR complex exerts, its genetic variants have been found to be associated with a variety of diseases/phenotypes, in which the risk can often depend on inter-individual variability or genetic differences within VDR protein; such like the heterogeneous phenotype of asthma. In fact, since VDR is located within the region q13-26 of chromosome 12 previously linked to asthma or related phenotypes in different populations [8–12], it is considered to be one of the candidate genes of asthma. A number of association studies were previously conducted in different populations and ethnic groups. Two of them have suggested associations among VDR polymorphisms and asthma. The first one was a family-based study on a founder French population from northeastern Quebec, accompanied by a second Canadian study, using both the childhood asthma management network as family-based study and the Healthy women from European ancestry as a case-control study [13, 14]. Another research among a German population has also tested the same hypothesis, however, there was no preferential transmission of VDR variants to children with asthma . Assays were also done on animal models, where an experimental allergic asthma was induced in VDR knockout and wild-type (WT) mice. As expected, WT mice developed symptoms of airway inflammation with an influx of eosinophils, elevated Th2 cytokine levels, mucous production, and airway hyperresponsiveness .
Thus, we sought to replicate the association previously described between the VDR gene and asthma among a Chinese Han population. Five Single nucleotide polymorphisms (SNPs) were selected, genotyped and statistically analyzed in a case-control study including 1090 Chinese individuals. Our results demonstrated that VDR rs7975232 (ApaI) variant was associated with asthma susceptibility in the cohort studied.
All subjects were Chinese of Han ethnic origin living in the same urban area. This group form 99% of the local population of Shandong, an eastern coastal province of China. Of note, the ethnicity of the group was also evident from the similar lifestyle and dietary habits. This research is one of the association studies carried out in our Institute, where several asthma susceptibility genes with different locations and functions were being tested to identify some novel linkages to the complex phenotype of asthma in the Chinese Han population. For that purpose, the process of collecting the samples was done thanks to collaboration with a number of provincial Hospitals. Our primary number of samples has increased since three years [17–19] to a larger number, of which 567 cases and 523 controls were collected for this study. Detailed information about age, sex, clinical characteristics and original medical institution were obtained for all participants (see Table 1). Healthy subjects, matched for age, were sampled from both Shandong Provincial hospital and Qilu hospital in Jinan city, where they had undergone pulmonary function tests, by which no doctor-diagnosed asthma nor history of asthma or other pulmonary diseases were noted. Whereas, affected subjects were recruited from both Qilu Hospital and Weifang Asthma Hospital. The 314 females and 253 males all confirmed asthma diagnosis by spirometry, which was performed according to American Thoracic Society recommendations (ATS, 1987); while no other atopy-related diseases or allergy phenotypes were available for those patients. Blood samples were obtained from all participants with informed consent. Genomic DNA was extracted from peripheral blood leukocytes by a standard salting-out method. The study was approved by the ethics review committee for human studies at School of Medicine, Shandong University.
SNPs Selection and Genotyping
In this Case-control study, we tested five polymorphisms of VDR gene, which maps to chromosome 12q13.11, consists of 9 exons with at least 6 isoforms of exon 1, spans 63.5 kb, and encodes a 427-amino acid protein . All variants are candidate SNPs, which are known to have a specified extent for asthma and/or other pulmonary diseases described in the earlier researches among different populations: rs2228570 or FokI , rs3782905 [13, 14], rs1544410 or BsmI , rs7975232 or ApaI [14, 22], and rs731236 or TaqI [14, 15, 23]. Two synonymous sites were primary selected, rs2228570 (a T/C transition polymorphism in exon2) and rs731236 (a single base change C to T in codon 352 at the 3' end of the VDR gene). Further, the three other sites were investigated in the non-coding regions of VDR gene, including rs7975232 and rs1544410 both in intron 8, and rs3782905 in intron 2. All markers were analyzed by a standard PCR-restriction fragment length polymorphism (PCR-RFLP) method. The different loci were recognized by the following restriction endonucleases (New England Biolabs NEB-China, Beijing): FokI for C/T at 37°C in exon 2, TaqI for T/C at 65°C in exon 9, BsmI for G/A at 65°C and ApaI for C/A at 37°C in intron 8, respectively and DdeI for C/G at 37°C in intron 2 (primers, product lengths and cycling conditions are all summarized in Table 2).
Single nucleotide polymorphisms were assessed for both genotypic and allelic association analysis among asthma patients and healthy controls. The genotype data of all tested SNPs were then used to estimate Hardy-Weinberg equilibrium by comparison of genotype frequencies within the two groups by a χ2-test. P-values were calculated with SPSS statistical software (SPSS Statistics 17.0; 2008 SPSS Inc, Chicago, IL), where Fisher's exact test was applied to analyze the comparison of the frequencies of discrete variables between cases and controls. The odds ratios (ORs) with 95% confidence intervals (95%CI) were also calculated during the single site analysis to estimate risk of asthma associated with the VDR polymorphic genotypes. An adjustment for multiple tests by Bonferroni correction in which the P-values are multiplied by the number of comparisons was later used to control the false discovery rate. Haplotypes and their frequencies were also constructed using the online software SHEsis http://analysis.bio-x.cn/myAnalysis.php, in which differences of haplotype distribution between patients and controls is assessed by the Monte Carlo method. The differences in pairwise distribution between respective pairs of SNPs tested were assessed by the statistical measure of LD based on calculating D' using Haploview 4.1. Further, a graphical representation of different genotype data obtained with SNP rs7975232 was produced, using a powerful program for analyzing graphs (Graphical Analysis 3.4), in which the genotype frequencies of cases and controls obtained with this site were diagrammatically compared.
For the five polymorphisms tested in this study within the VDR gene sequence, genotype distributions were in Hardy-Weinberg equilibrium in both patients and controls (all P-value > 0.05). The first association result with the coding polymorphisms VDR/rs2228570 and VDR/rs731236, however negative, was important since it excluded the susceptible risk. In fact, the statistical analysis of these two SNPs did not call for any association with asthma (P > 0.05). In the second part, three other SNPs were tested within non-coding regions of VDR sequence, VDR/rs3782905, VDR/rs7975232 and VDR/rs1544410. In contrast to the marker rs7975232, the role of the intronic variants rs3782905 and rs1544410 in asthma could not be confirmed. Hence, among the five markers studied, only intronic SNP rs7975232 (ApaI) was standing out as a significant polymorphism with an increased risk in one group compared to the other (The genotype analysis and the results of association tests for the five markers are all represented in Table 3).
Importantly, rs7975232 has shown a strong association with asthma in our Chinese samples, and the risk conferred by it was slightly excessive (P = 0.009; Odds Ratio = 1.33; 95% CI: 1.10–1.60). In addition, after graphical representation of its different genotypic frequencies distribution, the minor allele has demonstrated protective effect in a recessive model in affected subjects, of which the corresponding homozygote genotype AA was underrepresented in the case data comparing to the control data (Figure. 1).
Genotypic and allelic P-values were also calculated for rs7975232 site through a detailed statistical analysis (see Table 4). Also, when the difference for the same variant was adjusted for multiple testing by the Bonferroni method, the global significant association obtained was still persisting (P = 0.045). One other significance level was obtained with AA homozygote individuals (P = 0.004) in this analysis, and a highest odds ratio was also estimated for the same genotype (Odds Ratio = 1.88, 95% CI: 1.22–2.90). Consequently, these results suggest in all likelihood that minor allele of rs7975232 marker is a risk allele for asthma in our Chinese cohort studied.
Finally, evaluation of LD between each pair of loci has determined one block on the LD plot structure ("Block 1"; see figure 2A/B), which spans a region of 26 kb between two SNPs (rs3782905 and rs1544410); from intron 2 (position 46552434) to intron 8 (position 46526102). However, LD in this region was considered low, which is perhaps due to the low frequency seen with the minor alleles of these two SNPs, since they may not affect the "Block 1". Whereas, further in the 3'-terminal region more significant D' measures were seen between the two pairs rs1544410-rs7975232 and rs7975232-rs731236 (see figure 2A). This can probably explain the association observed with the variant rs7975232 in this region. In addition to this result, to investigate whether rs7975232 or neighboring non-polymorphic loci were causative, haplotype association analysis was later conducted with the complete five-locus haplotype studied (see Table 5). The analysis has shown a significant global association result with asthma (χ 2 = 14.51, P = 0.012), where the selected haplotype TCGAT has demonstrated the highest significance association (P = 0.0005).
This research directs an issue to test the hypothesis described in the previous association studies between VDR variants and asthma. Though the susceptibility for asthma disease could not totally be confirmed for all SNPs tested in this study, the significant result obtained with the non-coding polymorphism rs7975232 has reported previously unknown information about the relation existing between VDR gene and asthma in the Chinese Han population. This may also increase the probability of VDR to be a candidate gene for asthma.
Our results represent also a replication of the previous findings in the two U.S. studies, where SNP rs7975232 polymorphism has also shown a significant association with asthma. In the CAMP, seven VDR polymorphisms were genotyped in complete nuclear families from three different ethnic groups included in the analysis. After a stratified genetic analysis by ethnic groups, only intronic variant rs7975232 has shown a significant association with asthma (P = 0.01) . However, the NHS on 517 cases and 519 controls has shown evidence for association with asthma of four out of six VDR SNPs (rs7975232, rs731236, rs2239185, and rs3782905) . These results have demonstrated a replication of rs7975232 VDR-locus association with asthma among the Canadian women samples. However, the two polymorphisms rs731236 and rs3782905, which have reported a significant association in the NHS, have not shown any association with asthma in both Chinese population and Canadian families from the CAMP study. Elsewhere, Poon and colleagues addressed a study which involved a total of 223 independent French Northeastern Quebec families . In contrast to our result, their results have not demonstrated any association between VDR rs7975232 association and asthma; however, six other VDR variants widely studied were significantly overtransmitted to both asthmatic and atopic offspring (p < 0.05) in the French families. Another study addresses the previously described association of VDR variants with asthma and related phenotypic traits, in which 951 German individuals from 224 pedigrees were included in the analysis . Thirteen SNPs were tested and all reported no associations with asthma, including three non-significant SNPs also typed in our study (rs3782905, rs1544410 and the coding polymorphism rs731236).
Recently, researchers are following a key challenge to understand the functionality of different VDR polymorphisms, by the fact that most of these markers are synonymous. Most investigators have focused on the 3' regulatory region, because it is close to the synonymous markers tested mostly in association studies, such as BsmI, ApaI, and TaqI. The importance of N-terminal regions of the nuclear hormone receptors is increasingly recognized, especially through regulation of mRNA stability . Only several papers have tried to identify the relation existing between significant associations and functional sequence variations in the VDR gene. Although no consistent association studies support those functionality findings in asthma disease.
In summary, we identified a significant association between a genetic variant at the VDR locus and asthma in Chinese Han population. However, further studies are needed to find out the genuine functional variants in this region.
Marshall TG, Marshall FE: Sarcoidosis succumbs to antibiotics implications for autoimmune disease. Autoimmun Rev. 2004, 3: 295-300. 10.1016/j.autrev.2003.10.001.
Hayes CE: Vitamin D: a natural inhibitor of multiple sclerosis. Proc Nutr Soc. 2000, 59: 531-535.
Simmons JD, Mullighan C, Welsh KI, Jewell DP: Vitamin D receptor gene polymorphism: association with Crohn's disease susceptibility. Gut. 2000, 47: 211-214. 10.1136/gut.47.2.211.
Cutolo M, Otsa K, Uprus M, Sabrina Paolino S, Seriolo B: Vitamin D in rheumatoid arthritis. Autoimmunity Reviews. 2007, 7: 59-64. 10.1016/j.autrev.2007.07.001.
Swamy N, Xu WR, Paz N, Hsieh JC, Haussler MR, Maalouf GJ, Mohr SC, Ray R: Molecular modeling, affinity labeling, and site-directed mutagenesis define the key points of interaction between the ligand-binding domain of the vitamin D nuclear receptor and 1R,25-dihydroxyvitamin. Biochemistry. 2000, 39: 12162-12171. 10.1021/bi0002131.
Wang TT, Tavera-Mendoza LE, Laperriere D, Libby E, MacLeod NB, Nagai Y, Bourdeau V, konstorum A, Lallemant B, Zhang R, Mader S, White JH: Large-Scale in Silico and Microarray-Based Identification of Direct 1,25-Dihydroxyvitamin D3 Target Genes. Mol Endocrinol. 2005, 19: 2685-2695. 10.1210/me.2005-0106.
Xu H, Xiong DH, Xu FH, Zhang YY, Lei SF, Deng HW: Association between VDR ApaI polymorphism and hip bone mineral density can be modified by body mass index: a study on postmenopausal Chinese women. Acta Biochim Biophy Sin. 2005, 37: 61-67.
Malerba G, Lausiello MC, Scherpbier T, Trabettie E, Galavotti R, Cusin V, Pescollderungg L, Zanoni G, Martinati LC, Boner AL, Levitt RC, Pignatti PF: Linkage analysis of chromosome 12 markers in Italian families with atopic asthmatic Children. Am J Respir Crit Care Med. 2000, 162: 1587-1590.
Raby BA, Silverman EK, Lazarus R, Lange C, Kwiatkowski DJ, Weiss ST: Chromosome 12q harbors multiple genetic loci related to asthma and asthma-related phenotypes. Hum Mol Genet. 2003, 12: 1973-1979. 10.1093/hmg/ddg208.
Barnes KC, Freidhoff LR, Nickel R, Chiu YF, Juo SH, Hizawa N, Naidu RP, Ehrlich E, Duffy DL, Schou C, Levett PN, Marsh DG, Beaty TH: Dense mapping of chromosome 12q13.12-q23.3 and linkage to asthma and atopy. J Allergy Clin Immunol. 1999, 104: 485-491. 10.1016/S0091-6749(99)70398-2.
Daniels SE, Bhattacharrya S, James A, Leaves NI, Young A, Hill MR, Faux JA, Ryan GF, Le sousef PN, Lathrop GM, Musk AW, Cookson WO: A genome wide search for quantitative-trait loci underlying asthma. Nautre. 1996, 383: 247-50.
Wjst M, Fischer G, Immervoll T, Jung M, Saar K, Rueschendorf F, Reis A, Ulbrecht M, Gomolka M, Weiss EH, Jaeger L, Nickel R, Richter K, Kjellman N-IM, Griese M, Berg AV, Gappa M, Riedel F, Boehle M, Koningsbruggen SV: A Genome-wide Search for Linkage to Asthma. Genomics. 1999, 58: 1-8. 10.1006/geno.1999.5806.
Poon AH, Laprise C, Lemire M, Montpetit A, Sinnett D, Schurr E, Hudson TJ: Association of vitamin D receptor genetic variants with susceptibility to asthma and atopy. Am J Respir Crit Care Med. 2004, 170: 967-973. 10.1164/rccm.200403-412OC.
Raby BA, Lazarus R, Silverman EK, Lake S, Lange C, Wjst M, Weiss S: Association of vitamin D receptor gene polymorphisms with childhood and adult asthma. Am J Respir Crit Care Med. 2004, 170: 1057-1065. 10.1164/rccm.200404-447OC.
Wjst M: Variants in the vitamin D receptor gene and asthma. BMC Genetics. 2005, 6: 2-10.1186/1471-2156-6-2.
Wittke A, Weaver V, Mahon BD, August A, Cantorna MT: Vitamin D receptor-deficient mice fail to develop experimental allergic asthma. J Immunol. 2004, 173: 3432-3436.
Wang P, Liu QJ, Li JS, Li HC, Wei CH, Guo CH, Gong YQ: Lack of association between ADAM33 gene and asthma in a Chinese population. Int J Immunogenet. 2006, 33: 303-306. 10.1111/j.1744-313X.2006.00617.x.
Li JS, Liu QJ, Wang P, Li HC, Wei CH, Guo CH, Gong YQ: Absence of association between two insertion/deletion coding genetic polymorphisms of TIM-1 gene and asthma in Chinese Han population. Int J Immunogenet. 2006, 33: 417-422. 10.1111/j.1744-313X.2006.00634.x.
Liu QJ, Shang LJ, Li JS, Wang P, Li HC, Wei CH, Gong YQ: A Functional Polymorphism in the TIM-1 Gene is associated with asthma in a Chinese Han population. Int Arch Allergy Immunol. 2007, 144: 197-202. 10.1159/000103992.
Miyamoto K, Kesterson RA, Yamamoto H, Taketani Y, Nishiwaki E, Tatsumi S, Inoue Y, Morita K, Takeda E, Pike JW: Structural organization of the human vitamin D receptor chromosomal gene and its promoter. Mol Endocrinol. 1997, 11: 1165-1179. 10.1210/me.11.8.1165.
Liu W, Cao WC, Zhang CY, Tian L, Wu XM, Habbema JD, Zhao QM, Zhang PH, Xin ZT, Li CZ, Yang H: VDR and NRAMP1 gene polymorphisms in susceptibility to pulmonary tuberculosis among the Chinese Han population: a case-control study. Int J Tuber C Lung Dis. 2004, 8: 428-434.
Bornman L, Campbell SJ, Fielding K, Bah B, Sillah J, Gustafson P, Manneh K, Lisse I, Allen A, Sirugo G, Sylla A, Aaby P, McAdam KPWJ, Bah-Sow O, Bennett S, Lienhardt C, Hill AVS: Vitamin D receptor polymorphisms and susceptibility to tuberculosis in West Africa: a case-control and family study. J Infect Dis. 2004, 190: 1631-41. 10.1086/424462.
Wilkinson R, Llewelyn M, Toossi Z, Patel P, Pasvol G, Lalvani A, Wright D, Latif M, Davidson R: Influence of vitamin D deficiency and vitamin D receptor polymorphisms on tuberculosis among Gujarati Asians in west London: a case-control study. Lancet. 2000, 355: 618-621. 10.1016/S0140-6736(99)02301-6.
Shi YY, He L: SHEsis, a powerful software platform for analyses of linkage disequilibrium, haplotype construction, and genetic association at polymorphism loci. Cell Res. 2005, 15: 97-8. 10.1038/sj.cr.7290286.
Barrett JC, Fry B, Maller J, Daly MJ: Haploview: Analysis and visualization of ld and haplotype maps. Bioinformatics. 2005, 21: 263-265. 10.1093/bioinformatics/bth457.
The pre-publication history for this paper can be accessed here:http://0-www.biomedcentral.com.brum.beds.ac.uk/1471-2350/10/71/prepub
We address a special thank to both asthmatic patients and healthy subjects for their participation in the study. This work was financially supported by grants of the national Natural Science Foundation of China (No. 30300310, 30671956) for Dr. Qiji Liu, National Hightech Research and Development Program of China (No. 2006AA02A406) and Cultivation Fund of the Key Scientific and Technical Innovation Project, Ministry of Education of China (No. 704030).
The authors declare that they have no competing interests.
QL conceived and designed the study. AS conducted the laboratory experiments, developed protocols, did the statistical analysis, and drafted the manuscript, GG helped genotyping, YG and QL polished the final version. CW and HL participated in the clinical survey and sample collections. All authors read and approved the final manuscript.
About this article
Cite this article
Saadi, A., Gao, G., Li, H. et al. Association study between vitamin D receptor gene polymorphisms and asthma in the chinese han population: a case-control study. BMC Med Genet 10, 71 (2009) doi:10.1186/1471-2350-10-71
- Qilu Hospital
- Asthma Susceptibility
- Eastern Coastal Province
- Experimental Allergic Asthma