- Research article
- Open Access
- Open Peer Review
Vitamin D binding protein variants associate with asthma susceptibility in the Chinese han population
BMC Medical Geneticsvolume 12, Article number: 103 (2011)
Asthma is a genetically heterogeneous disease. Polymorphisms of genes encoding components of the vitamin D pathway have been reported to associate with the risk of asthma. We have previously demonstrated that vitamin D status was associated with lung function in Chinese asthma patients. In this study, we tested whether polymorphisms of genes encoding for vitamin D receptor (VDR), vitamin D 25-hydroxylase (CYP2R1) and vitamin D binding protein (GC) were associated with asthma in the Chinese Han population.
We sequenced all 8 exons of VDR and all 5 exons of CYP2R1 in a Chinese case-control cohort of asthma consisting of 467 cases and 288 unrelated healthy controls. Two mutations were identified in these regions. These variants were specified as rs2228570 in exon 2 of VDR and rs12794714 in exon 1 of CYP2R1. We also genotyped two common polymorphisms in GC gene (rs4588 and rs7041) by a PCR-restriction fragment length polymorphism (RFLP) method. We analyzed the association between these 4 polymorphisms and asthma susceptibility and asthma-related traits.
Polymorphic markers in VDR and CYP2R1 were not associated with asthma in the Chinese Han cohort. Importantly, variants in GC gene, which give rise to the two most common electrophoretic isoforms of the vitamin D binding protein, were associated with asthma susceptibility. Compared with isoform Gc1, Gc2 was significantly associated with the risk of asthma (OR = 1.35, 95% CI = 1.01-1.78 p = 0.006).
The results provide supporting evidence for association between GC variants and asthma susceptibility in the Chinese Han population.
Asthma is a chronic inflammatory disease of the airways characterized by reversible airflow obstruction and airway hyper-reactivity. The disease exhibits complex patterns of transmission, and has an estimated heritability of 36 - 79% . To date, a plethora of genes and variants have been implicated in the pathogenesis and progression of asthma [2, 3]. However, relatively few of these genes and variants have been successfully replicated.
In recent years the vitamin D pathway has emerged as a new potential component of asthma and allergy susceptibility. The majority of this evidence has been derived from epidemiological surveys. Deficiency of vitamin D has been found to increase the risk of asthma or asthma exacerbation [4–8]; higher maternal vitamin D intake during pregnancy has been associated with a reduced risk of wheezing illnesses in young children [9–11]; and circulating 25-hydroxyvitamin D (25OHD) has been associated with lung function in patients with asthma [7, 12, 13]. More recently, a small randomized double-blinded 6 month trial indicated that vitamin D supplementation in children may reduce asthma exacerbation following acute respiratory infection .
At the genetic level, a number of polymorphisms have been identified in various components of the vitamin D pathway; including the vitamin D receptor VDR, the microsomal vitamin D hydroxylase CYP2R1, and the vitamin D binding protein GC. VDR is the major receptor for the bioactive form of 25OHD, 1α,25-dihydroxyvitamin D3 (1,25OHD). The gene encoding VDR is located on human chromosome 12q13 and has 8 coding exons encode 427 amino acids. CYP2R1, positioned on human chromosome 11p15, is responsible for transforming photo-synthesized and dietary vitamin D into 25OHD . CYP2R1 gene has 5 exons encoding 501 amino acids. Vitamin D metabolites in the circulation are bound to vitamin D binding protein (DBP) encoded by a highly polymorphic gene GC. GC is located on chromosome 4q13 and has 13 exons encoding 474 amino acids. There are two common functional single nucleotide polymorphisms (SNPs) in exon 11 of GC, which are rs4588 and rs7041 positioned at codons 416 and 420. These variants have been associated with different binding affinity for 25OHD, and give rise to the two most common electrophoretic isoforms of DBP; Gc1 (which can be further subdivided into Gc1F and Gc1S); and Gc2. Whilst GC isoforms have been implicated in chronic obstructive pulmonary disease (COPD) across a variety of contrasting population (for review, see reference 16), there is relatively little evidence of their involvement in asthma [16, 17]. VDR variants have been extensively studied in asthma, although the results have shown contradictory [15, 18–23]. Similarly, evidence for association between CYP2R1 variants and asthma has also been inconsistent between populations . One potential explanation for this inconsistency may be differences in the underlying pattern of linkage disequilibrium between typed markers and the true disease loci in distinct populations.
We previously showed that lower vitamin D status in serum was associated with impaired lung function in Chinese asthma patients . In this study we set out to evaluate whether the variants of genes encoding for the key components of vitamin D pathway, including CYP2R1, VDR, and DBP, were associated with asthma or asthma-related phenotypes in a case-control design study in the Chinese Han population.
Subjects and phenotypical characteristics
We recruited 467 asthma patients and 288 controls with Han ethnicity from the Northern region of China. All participants were unrelated. Asthma cases were recruited at the Pulmonary Clinic, Peking Union Medical College Hospital. Ethnicity-matched controls were selected from non-asthmatic, non-atopic, healthy individuals with normal lung function .
Measurement of serum total IgE was performed by applying the UniCAP System (Pharmacia, Uppsala, Sweden) according to the manufacturer's instructions.
Asthma was diagnosed using the criteria defined by international guideline . All cases complained of current asthma symptoms including wheezing, cough, awakening at night and shortness of breath. Asthmatics with forced expiratory volume in 1s (FEV1) ≥ l70% of predicted underwent the methacholine challenge test, and positive bronchial hyperresponsiveness was defined by a 20% fall in FEV1 at inhaled methacholine concentrations ≤16 mg/ml. Patients with FEV1 < 70% of predicted underwent an airway reversibility test with positive airway reversibility defined by both 12% improvement in FEV1 and an increase ≥ 200 ml in the absolute FEV1 value after β2-agonist inhalation .
The incidence of asthma was higher in females than in males (60.4% vs 39.6%). The baseline FEV1% predicted was significantly lower in asthmatics than normal controls (68.9 ± 37.9% vs 90.6 ± 14.9%, P < 0.001). Methacholine airway challenge test was performed in 224 (48%) asthma subjects with baseline FEV1> 70% of predicted; in all of these subjects, positive bronchial hyperresponsiveness was documented. The airway reversibility test was carried out in 273 patients (58%). The mean improvement in FEV1 was 26.4 ± 14.3%. The serum total IgE values in controls were within the normal range. Detailed characteristics of the individuals are listed in Table 1.
The study protocol was reviewed and approved by the Human Research Ethics Committee of the Peking Union Medical College Hospital, and all subjects gave written informed consent to participate in the study.
SNPs selection and genotyping
DNA samples (derived from 467 asthma patients and 288 healthy, non-atopic controls) were extracted from peripheral blood leukocytes using a standard phenol-chloroform method.
All selected genes code for key enzymes or components responsible for vitamin D metabolic pathway. Single nucleotide polymorphisms (SNPs) of genes tested in this study have been associated with asthma susceptibility . We directly sequenced the exons of VDR and CYP2R1 in all 755 individuals and identified two mutations, specified as rs2228570 in exon 2 of VDR and rs12794714 in exon 1 of CYP2R1. Sequencing VDR and CYP2R1 genes was performed at the Beijing Genomics Institute (http://www.genomics.cn). DNAStar 5.0 SeqMan Software was used to identify DNA variants and provide genotype calls. The primers for PCR were designed by using Primer Premier 5.0 Software. Sequences of the primers are listed in Table 2.
We genotyped the two most commonly studied SNPs in the GC gene, rs4588 and rs7041, using a standard PCR-restriction fragment length polymorphism (PCR-RFLP). Genetic variants in GC gene have been shown the different affinity in binding vitamin D metabolites . The SNPs rs4588 and rs7041 have been found to have associations with several chronic lung diseases, such as COPD and tuberculosis [27, 28]. Carriage of the genotype of Gc2 allele was associated with increased tuberculin-stimulated IFN-γ release in Gujarati Asian TB contacts . Homozygous carriers of the rs7041 T allele exhibited an increased risk for COPD in a Belgium population . There are only 11 base pairs apart between these two SNPs in the human genome, however, they are not in the same LD block. Thus we selected these two important SNPs in our Chinese samples to test the association with the risk of developing asthma. These two SNPs, rs4588 and rs7041, are positioned at codons 416 (GAT→GAG, Asp→Glu) and 420 (ACG→AAG, Thr→Lys) of exon 11 of the GC gene (Table 3) . The different loci were recognized by the following restriction endonucleases (New England Biolabs): Hae III for T/G at 37°C, Sty I for C/A at 37°C.
The PCR reactions were performed in a total volume of 30 μl (40 ul for GC variants genotyping) containing 70 ng of genomic DNA, 3 μl of 10 × PCR Buffer, 1.5 μM of MgCl2, 0.17 mM of dNTPs, 100 pM of each primer, and 1.5 U of Taq DNA polymerase (Takara). In general, PCR cycle condition consisted of an initial denaturation step at 94°C for 3 min followed by 30 cycles (35 cycles for GC variants genotyping) of denaturation for 30 s at 94°C, annealing for 30 s at 55-57°C, extension for 40- 45 s at 72°C, and a final extension at 72°C for 7 min.
Because of similar functional characteristics of GC isoforms Gc1F and Gc1S, allele carriers were combined to produce a total of 3 genotypes: Gc1/1, Gc1/2 and Gc2/2. Allele and genotype frequencies of polymorphisms in VDR, CYP2R1 and GC were obtained by direct counting, and Hardy-Weinberg equilibrium was evaluated by a chi-squared test (using SPSS statistical software, version 12.0). All genotype and phenotype data were formatted using PLINK.
Quantitative traits were tested using a variance components approach implemented in the statistical package Merlin . Quantile normalization was applied to all quantitative traits prior to analysis (with the exception of total serum IgE which had already been subjected to a log transformation). Missing genotypes were inferred using the --- infer option available in Merlin.
In addition, an unpaired t-test was used to compare the baseline data between asthma patients and healthy controls. One-way ANONA was used to assess the genotype distribution in asthma patients. A Bonferroni correction was used to adjust p-values for multiple testing (based on an alpha threshold of 0.05). Genotype and allele frequencies in cases and controls were also compared by contingency table analysis. Genotype relative risk was calculated according to the statistical method described by Lathrop . This method compares case genotype frequencies with expected control genotype frequencies under the assumption of Hardy-Weinberg equilibrium, and is more powerful than standard contingency table analysis.
We genotyped 1 SNP (rs12794714) in CYP2R1, 1 SNP (rs2228570) in VDR, and 2 SNPs (rs4588 and rs7041) in GC. The allele and genotype frequencies are listed in Table 4. All SNPs met criteria for Hardy-Weinberg equilibrium in both case and control groups (all p-value > 0.05).
No significant difference was observed in the genotype and allele frequencies of rs2228570 in VDR and rs12794714 in CYP2R1 between asthmatic cases and controls (Table 4), suggesting that these SNPs were not genetically associated with the risk of asthma.
Importantly, single marker analysis of association showed that GC polymorphism (expressed as Gc1/1, Gc1/2 and Gc2/2) was significantly associated with asthma susceptibility in our Chinese samples, the risk conferred by it was slightly elevated (OR = 1.35, 95% CI = 1.01 - 1.78, P = 0.006). Compared with Gc1/1 genotype, Gc2/2 genotype was strongly associated with the risk of asthma (OR = 13.13, 95% CI = 2.42-7.13, P = 0.001) (Table 4). These results may suggest that Gc2/2 genotype confers a significant risk for developing asthma.
Vitamin D deficiency has been defined as serum 25OHD < 50 nmol/l . We have previously demonstrated that vitamin D deficiency is common among Chinese patients with asthma . Thus, we retrospectively analysed the association between the GC polymorphisms and serum 25-OHD concentrations, however, no significant difference was observed (Table 5).
We further examined the evidence for association between the polymorphic markers in these three genes and a number of quantitative asthma-related traits using Merlin. We were not able to demonstrate any significant associations between the studied SNPs and blood eosinphil counts, log-transformed serum IgE or lung function measures, even after stratification by sex (data not shown).
We performed a case-control study assessing association between polymorphisms in 3 genes of the vitamin D pathway and asthma susceptibility in the Chinese Han population. We showed that the GC polymorphisms giving rise to the major electrophoretic isoforms of vitamin E binding protein (rs4588 and rs7041) were significantly associated with the risk of asthma, although not 25OHD concentrations in asthma patients. By sequencing all exons of VDR and CYP2R1 genes, we have identified that the full range of exonic genetic variants presented in our cohort have shown allele frequencies ≥0.13%. We observed only two exonic SNPs, rs2228570 (also known as Fok I) in VDR and rs12794714 in CYP2R1, both of which had been previously characterised. No association could be identified between either of these SNPs and asthma. Because exonic regions are more likely to carry variants that disrupt protein function [32, 33], our polymorphism screening did not include intronic or promoter regions of these genes which may also potentially contain causal variants .
The FokI polymorphism in exon 2 of VDR consists of a T to C change which has been reported to cause gene-specific and cell type-specific effects. VDR polymorphisms have been widely studied in asthma, and the results have been contradictory [15, 18–23]. The association of polymorphisms in VDR with asthma was first reported in Northern American family-based studies [19, 20], but initial positive results were not replicated in subsequent populations [22, 23]. Consistent with our negative results, a previous case-control study of rs2228570 performed in a Chinese Han population also failed to show an association with asthma .
The CYP2R1 variant, rs12794714, has previously been associated with serum 25OHD levels in a healthy Caucasian population . Alternative CYP2R1 variants have also been implicated in asthma susceptibility, though these effects have not replicated consistently among populations . This lack of consistency may reflect a linkage disequilibrium phenomenon (with non-coding regulatory variants(s)), a consequence of complex, undefined gene-gene (epistatic) or gene-environment interactions, or a low magnitude effect requiring large, highly powered cohorts for its detection.
The vitamin D binding protein (DBP), encoded by GC, is the principle plasma carrier protein of vitamin D and its metabolites. Due to a higher binding affinity of Gc1F and Gc1S to 25OHD, the isoforms are hypothesized to deliver 25OHD to target tissues more efficiently than Gc2 . The results presented here indicate that the two common SNPs in GC which determine DBP isoforms associate significantly with susceptibility to asthma, and that the Gc1 allele might confer a protective effect. One previous study of GC variants in asthma did not find an association, however, this study included a very small number of participants and was probably underpowered to detect a difference .
Polymorphisms in GC gene have been found to be significantly associated with 25OHD concentrations by genome-wide association (GWA) [36, 37]. We have previously shown that lower serum 25OHD concentrations are associated with impaired lung function measurements in Chinese patients with asthma . In this study, we reported the association between genetic variants in GC and asthma susceptibility. Thus, we further investigated the influence of GC genotypes on serum 25-OHD concentrations in asthma patients. When performing multiple tests for the association analysis with Bonferroni correction to the dataset, we couldn't find any significant association. One explanation could be that circulating serum 25OHD concentrations are influenced by several genetic loci [36, 37].
In conclusion, this study has shown that genetic variants in the gene encoding the vitamin D binding protein (GC) contribute to asthma susceptibility in a Han Chinese population. The OR in this study with respect to Gc2/2 homozygote suggests a modest but definite genetic effect. We have not found associations between that exonic variants in the VDR and CYP2R1genes and asthma or asthma-related traits in our Chinese population. Further research is needed to understand how gene-gene or gene-environment interactions and causative regulatory variants (intronic, promoter, or trans-acting) influence vitamin D metabolism in asthma patients .
Willis-Owen SA, Cookson WO, Moffatt MF: Genome-wide association studies in the genetics of asthma. Curr Allergy Asthma Rep. 2009, 9 (1): 3-9. 10.1007/s11882-009-0001-x.
Vercelli D: Discovering susceptibility genes for asthma and allergy. Nat Rev Immunol. 2008, 8 (3): 169-82. 10.1038/nri2257.
Ober C, Hoffjan S: Asthma genetics 2006: the long and winding road to gene discovery. Genes Immun. 2006, 7 (2): 95-100. 10.1038/sj.gene.6364284.
Devereux G, et al: A case-control study of vitamin D status and asthma in adults. Allergy. 65 (5): 666-7.
Brehm JM, et al: Serum vitamin D levels and markers of severity of childhood asthma in Costa Rica. Am J Respir Crit Care Med. 2009, 179 (9): 765-71. 10.1164/rccm.200808-1361OC.
Brehm JM, et al: Serum vitamin D levels and severe asthma exacerbations in the Childhood Asthma Management Program study. J Allergy Clin Immunol. 126 (1): 52-8 e5.
Searing DA, et al: Decreased serum vitamin D levels in children with asthma are associated with increased corticosteroid use. J Allergy Clin Immunol. 125 (5): 995-1000.
Shaheen SO: Vitamin D deficiency and the asthma epidemic. Thorax. 2008, 63 (3): 293-author reply 293
Devereux G, et al: Low maternal vitamin E intake during pregnancy is associated with asthma in 5-year-old children. Am J Respir Crit Care Med. 2006, 174 (5): 499-507. 10.1164/rccm.200512-1946OC.
Litonjua AA, et al: Maternal antioxidant intake in pregnancy and wheezing illnesses in children at 2 y of age. Am J Clin Nutr. 2006, 84 (4): 903-11.
Camargo CA, et al: Maternal intake of vitamin D during pregnancy and risk of recurrent wheeze in children at 3 y of age. Am J Clin Nutr. 2007, 85 (3): 788-95.
Li F, et al: Vitamin d deficiency is associated with decreased lung function in chinese adults with asthma. Respiration. 81 (6): 469-75.
Sutherland ER, et al: Vitamin D levels, lung function, and steroid response in adult asthma. Am J Respir Crit Care Med. 181 (7): 699-704.
Majak P, et al: Vitamin D supplementation in children may prevent asthma exacerbation triggered by acute respiratory infection. J Allergy Clin Immunol. 127 (5): 1294-6.
Bosse Y, et al: Asthma and genes encoding components of the vitamin D pathway. Respir Res. 2009, 10: 98-10.1186/1465-9921-10-98.
Chishimba L, et al: The vitamin D axis in the lung: a key role for vitamin D-binding protein. Thorax. 2010, 65 (5): 456-62. 10.1136/thx.2009.128793.
Laufs J, et al: Association of vitamin D binding protein variants with chronic mucus hypersecretion in Iceland. Am J Pharmacogenomics. 2004, 4 (1): 63-8. 10.2165/00129785-200404010-00007.
Fang WL, et al: Association analysis of vitamin D receptor gene polymorphisms in chinese population with asthma. Iran J Allergy Asthma Immunol. 2009, 8 (3): 141-7.
Raby BA, et al: Association of vitamin D receptor gene polymorphisms with childhood and adult asthma. Am J Respir Crit Care Med. 2004, 170 (10): 1057-65. 10.1164/rccm.200404-447OC.
Poon AH, et al: Association of vitamin D receptor genetic variants with susceptibility to asthma and atopy. Am J Respir Crit Care Med. 2004, 170 (9): 967-73. 10.1164/rccm.200403-412OC.
Saadi A, et al: Association study between vitamin D receptor gene polymorphisms and asthma in the Chinese Han population: a case-control study. BMC Med Genet. 2009, 10: 71-
Vollmert C, et al: Single nucleotide polymorphism screening and association analysis--exclusion of integrin beta 7 and vitamin D receptor (chromosome 12q) as candidate genes for asthma. Clin Exp Allergy. 2004, 34 (12): 1841-50. 10.1111/j.1365-2222.2004.02047.x.
Wjst M: Variants in the vitamin D receptor gene and asthma. BMC Genet. 2005, 6: 2-
Gao J, et al: Polymorphisms of PHF11 and DPP10 are associated with asthma and related traits in a Chinese population. Respiration. 79 (1): 17-24.
Bateman ED, et al: Global strategy for asthma management and prevention: GINA executive summary. Eur Respir J. 2008, 31 (1): 143-78. 10.1183/09031936.00138707.
Arnaud J, Constans J: Affinity differences for vitamin D metabolites associated with the genetic isoforms of the human serum carrier protein (DBP). Hum Genet. 1993, 92 (2): 183-8.
Martineau AR, et al: Association between Gc genotype and susceptibility to TB is dependent on vitamin D status. Eur Respir J. 35 (5): 1106-12.
Janssens W, et al: Vitamin D deficiency is highly prevalent in COPD and correlates with variants in the vitamin D-binding gene. Thorax. 65 (3): 215-20.
Abecasis GR, et al: Merlin--rapid analysis of dense genetic maps using sparse gene flow trees. Nat Genet. 2002, 30 (1): 97-101. 10.1038/ng786.
Lathrop GM: Estimating genotype relative risks. Tissue Antigens. 1983, 22 (2): 160-6.
Holick MF: Vitamin D deficiency. N Engl J Med. 2007, 357 (3): 266-81. 10.1056/NEJMra070553.
Cotton RG, Scriver CR: Proof of "disease causing" mutation. Hum Mutat. 1998, 12 (1): 1-3. 10.1002/(SICI)1098-1004(1998)12:1<1::AID-HUMU1>3.0.CO;2-M.
Kryukov GV, Pennacchio LA, Sunyaev SR: Most rare missense alleles are deleterious in humans: implications for complex disease and association studies. Am J Hum Genet. 2007, 80 (4): 727-39. 10.1086/513473.
Loots GG, et al: Identification of a coordinate regulator of interleukins 4, 13, and 5 by cross-species sequence comparisons. Science. 2000, 288 (5463): 136-40. 10.1126/science.288.5463.136.
Bu FX, et al: Comprehensive association analysis of nine candidate genes with serum 25-hydroxy vitamin D levels among healthy Caucasian subjects. Hum Genet. 2010, 128 (5): 549-56. 10.1007/s00439-010-0881-9.
Wang TJ, et al: Common genetic determinants of vitamin D insufficiency: a genome-wide association study. Lancet. 376 (9736): 180-8.
Ahn J, et al: Genome-wide association study of circulating vitamin D levels. Hum Mol Genet. 19 (13): 2739-45.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2350/12/103/prepub
This work was funded by New Century Excellent Talents program from Education Ministry of China (NCET-06-0156) and a grant from National Natural Sciences Foundation of China (No.30470767).
The authors declare that they have no competing interests.
FL performed the whole procedure of the experiments. LJ helped in performing the experiment. SAWO and YZh performed the statistical analysis and drafted the manuscript. JG designed and supervised the experiment, and drafted the manuscript. all authors read and approved the final manuscript.