Abstract
Aims: The CD36 gene encodes for a membrane receptor that facilitates fatty-acid uptake and utilization. Genetic variants of the CD36 gene have been associated with metabolic syndrome (MetS). We aimed to evaluate the association between the rs10499859A>G and rs13246513C>T polymorphisms and MetS components. Methods: For this case-control study, 140 MetS and 187 normal subjects were randomly selected from the Tehran Lipid and Glucose Study participants. Biochemical and anthropometrical variables were measured. Genotyping for both single nucleotide polymorphisms (SNPs) was performed by polymerase chain reaction-restriction fragment length polymorphism. Results: Case and control groups were not different in allele and genotype frequencies for these SNPs. However, the A and T alleles of these SNPs were significantly associated with elevated levels of high-density lipoprotein cholesterol (HDL-C) before age and sex adjustment (p=0.027 and 0.016, respectively). Association between the A allele and body mass index (BMI) was also significant after adjustment for MetS under the dominant model (p=0.009, β2=0.68). Conclusions: Based on our results, these polymorphisms do affect HDL-C level and BMI (MetS components), although the effect may be slight and restricted specifically to an environment-genotype.
Introduction
M
Because of the prevalence of MetS in Iran (32.1% in adults) (Zabetian et al., 2007) and the absence of association studies on CD36 gene in our population, this study was designed to examine the association of rs10499859A>G and rs13246513C>T polymorphisms with the MetS and its components. These single nucleotide polymorphisms (SNPs) were chosen because of their frequencies (minor allele frequency ≥0.05) and based on the results from a previous study, which suggested association between these SNPs and MetS components.
Materials and Methods
For this case-control study, 327 subjects, aged 19-70 years, were randomly selected from the Tehran Lipid and Glucose Study (TLGS). In brief, the TLGS is a longitudinal study, which is designed to determine the major risk factors for noncommunicable diseases among the Tehran urban population and to develop population-based measures to change lifestyles in order to prevent the rising trend of diabetes mellitus, dietary disorders, and dyslipidemia. In TLGS, the study sample was selected using multistage stratified cluster sampling from the population of Tehran. TLGS subjects were interviewed privately. Information on age; demographic and biochemistry data; physical activity status; education levels; and medication usage for the treatment of diabetes, hypertension, and lipid disorders were collected in the TLGS questionnaire. For the present study, cases and controls were randomly selected from the TLGS participants with and without MetS (based on participant information and MetS definition), respectively. A detailed description of the methodology, rationale, and design of this study has been published elsewhere in detail (Azizi et al., 2002).
This study was approved by the Medical Ethics Committee of the Research Institute for Endocrine Sciences at Shahid Beheshti University of Medical Science and informed consent was given by each participant.
The protocols for collection of clinical data have been described previously (Azizi et al., 2003; Daneshpour et al., 2010). Fasting serum glucose, triglycerides (TG), total cholesterol, and HDL-cholesterol (HDL-C) were measured using the enzymatic colorimetric method (Pars Azmon) by a Selectra 2 auto-analyzer (Vital Scientific). The Friedewald equation was used to calculate LDL-cholesterol; samples with TG greater than 400 mg/dL were assayed by the direct method. HDL subfractions were separated by differential polyanion precipitation. In all of these biochemical analyses, the coefficients of variation were less than 5%. For CD36 polymorphism analysis, genomic DNA was extracted as described previously (Truett et al., 2000; Daneshpour et al., 2008). The polymerase chain reaction (PCR) was used to amplify a 555-bp fragment for rs13246513C>T and 453 bp for rs10499859A>G in the CD36 gene using the oligonucleotide primers F: 5′-TCT ACC ACT TTG TTC TAG GCT AAT TTT T-3′, R: 5′-CAA ATA TCC TAA TGC AAA TAG AAT AAA CC-3′, F: 5′-CTT ATT TGG ATG GTA GGT TTG ACA CAG G-3′, and R: 5′-GAG GAA AGA AAG CAA CAT TTC AAG ACT TC-3′, respectively. Amplification was performed using the following program after initial denaturation at 95°C for 3 min: 95°C for 1 min, 60°C (for rs10499859A>G) or 55°C (for rs13246513C>T) for 45 s, at 72°C for 45 s for 30 cycles (Corbett). Each amplification was performed using 200 ng of total genomic DNA in a final volume of 15 μL containing 40 pmol of each oligonucleotide, 0.2 mM of each deoxyribo-nucleotide triphosphate (dNTP), 1.5 mM magnesium chloride, 10 mM Tris (pH 8.4), and 0.25 units of Taq polymerase (Fermentase). The PCR products were subjected to restriction digestion analysis. For the rs10499859A>G polymorphism, digestion with HincII (Fermentase) resulted in 284-bp and 169-bp fragments for G allele and a 453-bp fragment for A allele. For rs13246513C>T polymorphism, digestion with TaqI (Fermentase) resulted in 294-bp and 261-bp fragments for C allele and a 555-bp fragment for the T allele. The fragments were separated by electrophoresis on 2% agarose gels and DNA fragments were visualized by gel documentation (Optigo). Ten percent of the PCR samples were directly sequenced to confirm the polymerase chain reaction-restriction fragment length polymorphism results.
Definition
MetS was defined according to the joint interim statement criteria (Alberti et al., 2009) with the presence of any three of five risk factors of the following: (1) abdominal obesity as waist circumference ≥95 cm in men and women according to population- and country-specific cutoff point for Iranians (Azizi et al., 2010); (2) fasting blood sugar (FBS) ≥100 mg/dL or drug treatment; (3) fasting TG ≥150 mg/dL or drug treatment; (4) fasting HDL-C <50 mg/dL in women and <40 mg/dL in men or drug treatment; (5) raised blood pressure defined as systolic blood pressure ≥130 mmHg, diastolic blood pressure (DBP)≥85 mmHg, or antihypertensive drug treatment. Those subjects with two or fewer risk components were assigned as controls.
Statistical analysis
Continuous variables are presented as mean with standard deviation (SD), ordinal variables as median with range, and categorical variables as frequency and percentage. Genotype and allele frequencies between the groups were assessed by the χ2-test for Hardy-Weinberg equilibrium (p<0.05). A general linear model was used to compare mean values of anthropometric markers. In this case, the adjustment variables were age and gender. The mean values and SDs of biochemical variables in different groups (case, control, and genotype groups) were obtained and compared using an independent sample t-test. Logarithmic and Cox transformations were performed to normalize the error distribution and stabilize the error variance. The Mann-Whitney U test was used by considering five genetic models: a codominant model (three genotype groups separated), a dominant model (heterozygotes grouped with the homozygotes for the minor allele), a recessive model (heterozygotes grouped with the homozygotes for the major allele), a log-additive model (a score was assigned counting the number of minor alleles: score 0 for the homozygotes of the major allele, score 1 for the heterozygotes, and score 2 for the homozygotes of the minor allele), and an over-dominant model (homozygotes for the major allele grouped with the homozygotes for the minor allele). The Akaike information criterion was used to choose the genetic model that best fits the data. Multivariable linear regression was used to control the effect of potential confounders, including MetS, sex, age, and body mass index (BMI). All statistical analyses were performed using SPSS version 16 (SPSS). The statistical significance threshold was set to p≤0.05 (two-tailed test).
Results
The clinical characteristics and genotype frequencies of the 339 selected subjects are shown in Table 1. As expected, cases were significantly different from controls with regard to the variables that have been defined as MetS criteria (all p<0.05). Significant differences in HDL subfractions (HDL-2 and HDL-3) were also seen between case and control groups. Genotype distributions of the two CD36 SNPs in MetS and normal subjects are shown in Table 2. Genotypes and allele distributions for rs10499859A>G polymorphism were in Hardy-Weinberg equilibrium but for rs13246513C>T they deviate from the equilibrium (χ2=4.93; p=0.0264). Genotype and allele frequencies for the studied SNPs did not differ between subjects with and without MetS.
Mean±SD.
BMI, body mass index; HDL, high-density lipoprotein; LDL, low-density lipoprotein; SD, standard deviation.
Differences between MetS components, except HDL-C, were not significant between different allele carriers (Table 3). The data show that HDL-C level was significantly higher in A allele carriers of rs10499859 (A carrier: 47.1±1.26 mg/dL vs. G carrier: 43.9±1.27 mg/dL; p=0.027) and T allele carriers of rs13246513 (T carrier: 47.2±11.9 mg/dL vs. C carrier: 43.8±9.61 mg/dL; p=0.016; HDL). These associations were not retained following covariates adjustment for MetS, sex, age, and BMI by dominant model (p=0.098 and p=0.751, respectively). The means of the biochemical variables and BMI in genotype groups for each of the SNPs, under different genetic models for confounders' analysis, are presented in Table 4. The association of rs13246513C>T polymorphism with BMI, TG, and cholesterol before and after adjustment was not significant. No significant association was observed between rs10499859A>G alleles and BMI before adjustment for MetS (A carrier: 27.5±5.71 kg/m2 vs. G carrier: 26.6±5.16 kg/m2; p=0.208). This relationship became significant (p=0.009) after adjustment for MetS, as a common confounding factor. DBP was significantly increased in the G carriers (A carrier: 73.5±1.15 mmHg vs. G carrier: 96.2±1.21 mmHg; p=0.021) under recessive model. This association was no longer significant after adjustment for MetS, sex, and age (p=0.133). The association of rs10499859A>G polymorphism with other variables was not significant.
Mean±SD.
Median (CI 25%).
SD, standard deviation; CI, confidence interval.
pa, adjusted for MetS; pb, adjusted for MetS, sex, and age; pc, adjusted for MetS, sex, age, and BMI.
MetS, metabolic syndrome.
Discussion
The goal of this study was to evaluate the relation between rs10499859A>G in the 5′ untranslated region (5′UTR) and rs13246513C>T in the 3′flanking region of CD36 gene with MetS components. The results showed that these variants, although slightly, do affect MetS components, particularly HDL-C level. Carriers of A and T alleles of these SNPs had increased level of HDL-C, although the differences in levels were not statistically significant after adjustment for sex and age. This is, to our knowledge, the first report for the correlation between the CD36 gene polymorphism and MetS in Iranians. To establish this correlation other polymorphisms in this gene should be assessed.
There is a large amount of evidence that suggest low HDL-C is one of the major determinants of the MetS and an independent risk factor for cardiovascular disease (Coon et al., 2001). Therefore research has been carried out to investigate the association between variation in certain genes, like CD36, and HDL-C concentration. A 2008 study of 2020 African-American subjects by Love-Gregory et al. (2008) showed that some CD36 variants, including rs10499859A>G, were related to the levels of HDL-C and their minor alleles were accompanied by higher levels of HDL-C. Their results implied that the magnitude of the increase in HDL-C was positively associated with the number of minor alleles present. They reported an increase of 2.5 mg/dL in HDL-C per G allele of rs10499859 SNP. They also showed that SNP rs13246513 increased the odds for the MetS by 32%.
Although our results reveal the relationship between aforementioned SNPs and HDL-C levels too, but indicate the role of the major A allele at rs10499859 in HDL-C increase and also a similar role for minor T allele at rs13246513. The divergent results between our study and those reported in African-American population could be due on one hand to the difference in the genetic background and on the other hand to the small sample size of our study.
The results of studies that assessed the relationship between CD36 gene SNPs and BMI, as a measure of obesity, revealed a positive relation and suggest that genetic variation within the CD36 locus may contribute to metabolic disease via its effect on body adiposity (Bokor et al., 2010; Heni et al., 2011). The reported result by Yun et al. (2007) on a Korean population showed association of a polymorphism in intron 3 of the CD36 gene with an increase in BMI. Our results are in concordance with their results that the mean of BMI in A allele carriers was significantly greater than the respective value in G allele carriers.
Although we found no significant differences in allele frequency between cases and controls, and this may indicate that there is no association between these polymorphisms and MetS, we should bear in mind that the MetS is a multifactorial disorder in which each involved genetic factor has a slight and cumulative effect, and then the results are not unexpected. Also we should remember that a certain amount of the observed variation is in fact due to the different effect of environmental factors on different age groups, or the differences in lipid profile caused by hormonal background of the two genders. In small samples, the above effects will efficiently minimize and mask the true genetic differences between subjects.
This study has certain limitations. First, the small size of our study population might not have had sufficient power to detect direct evidence. Second, in spite of many polymorphisms that are detected, only two SNPs were evaluated in this region. Finally, only the genotype and allele distributions of rs10499859A>G SNP were in Hardy-Weinberg equilibrium, therefore haplotype analysis was not possible.
In conclusion, this study provides some evidence for identification of an effect of CD36 gene polymorphisms on MetS components and is a step on a path toward recognition of underlying genetic factors, the identification of which could lead to the prevention of MetS and therefore cardiovascular disease. Further genetic studies in a larger population should be carried out to confirm this observation.
Footnotes
Acknowledgments
We are indebted to each of the study participants for the substantial time and effort contributed to this study. We also wish to thank Mrs. N. Shiva for editing the manuscript. This work was supported by the Iran National Science Foundation under grant no. 298.
Disclosure Statement
No competing financial interests exist.
