Abstract
Background
Metabolic syndrome (MetS) is a complex and multifactorial disorder characterized by insulin resistance, dyslipidaemia, hyperglycemia, abdominal obesity, and elevated blood pressure. The apolipoprotein A5 (APOA5) gene variants have been reported to correlate with two major components of MetS, including low levels of high density lipoprotein cholesterol (HDL-C) and high levels of triglyceride. In the present study, we explored the associations between five single nucleotide polymorphisms (SNPs) of APOA5 gene and the MetS risk.
Methods
In a case-control design, 120 Iranian children and adolescents with/without MetS were genotyped by polymerase chain reaction-sequencing for these SNPs. Then, we investigated the association of SNPs, individually or in haplotype constructs, with MetS risk.
Metabolic syndrome (MetS) is clustering of several cardiovascular disease risk factors such as dyslipidaemia, hyperglycemia, abdominal obesity and elevated blood pressure [12]. Several complex etiological mechanisms such as genetics and environmental factors are considered to be related to MetS [2].
MetS is considerably influenced by genetic factors. Apolipoprotein A5 (APOA5) gene is located on chromosome 11q23 in the APOA1-C3-A4 gene cluster encoding the APOAV protein with 343 amino acids. APOA5 plays an important role in lipid metabolism and homeostasis. Moreover, plasma APOA5 levels have been considered to reduce plasma triglyceride (TG) and increase high density lipoprotein cholesterol (HDL-C) concentrations. Several single nucleotide polymorphisms (SNPs) in the APOA5 gene have been shown to be associated with increased TG levels [345].
MicroRNAs (miRNAs) are a group of non-protein-coding small RNAs regulating gene expression post-transcriptionally [67]. SNPs in miRNA target site (miR-TS-SNP) may affect or impair binding of miRNAs [8] and considerably affect diverse biological functions, including disease susceptibilities (cancers, Parkinson disease, osteoporosis, diabetes, and hypertension etc.) in humans [78910].
In the current study, a panel of five SNPs (rs619054, rs369740771, rs376976905, rs192708363, and rs34089864) within APOA5 gene was selected to study their association with MetS. Rs619054 variant is significantly associated with HDL, low density lipoprotein cholesterol (LDL-C) and TG levels [111213]. This polymorphism disrupts hsa-mir-28-5p, hsa-mir-708, and hsa-mir-3139 binding sites, and creates a new binding site for hsa-mir-320e. Rs34089864 variant destroys hsa-mir-4511 and hsa-mir-498 binding sites, and creates new binding sites for hsa-mir-4740-5p, hsa-mir-561-5p, and hsa-mir-3169.
The main aim of our case-control study was to evaluate the genotype and allele frequencies of two variants in the three prime untranslated regions (3’UTRs) of APOA5 gene, including rs619054 (c.⋆31C>T or 1764C>T) and rs34089864 (c.⋆76C>T or c.1177C>T) to find any potential associations with MetS and multiple features of this syndrome. Furthermore, we determined the haplotype frequencies and the associations of the five SNPs at the APOA5 locus with the risk of MetS.
This case-control study included 57 patients (24 males and 33 females: mean age, 12.40±0.26 years [range, 9 to 18 years]) and 63 controls (28 males and 35 females: mean age, 13.39±0.38 years [range, 8 to 18 years]) from children and adolescent of Isfahan (a central city in Iran) [14]. According to the modified Adult Treatment Panel III (ATP III) criteria [15], the MetS was defined as clustering of the following three or more factors in one person (1) fasting TG ≥100 mg/dL; (2) HDL-C <50 mg/dL; (3) waist circumference >75th percentile for age and gender in the studied population; (4) systolic blood pressure/diastolic blood pressure >90th percentile for gender, age, height [16]; and (5) fasting blood sugar ≥100 mg/dL. An oral assent was obtained from participants, and all parents signed a written informed consent after explaining the study objectives and details. Blood samples were collected in two separate tubes for DNA analysis and biochemical tests.
DNA was extracted from blood leukocytes using the Diatome kit (Isogen Laboratory, Moscow, Russia) according to the manufacturer's recommendations. Genotyping for the 5 miR-TS-SNP at APOA5 gene was carried out by polymerase chain reaction (PCR) and direct sequencing. Primers A5 forward (5′-CTGCCCAGGCCCATTCCCA-3′) and A5 reverse (5′-GAC CTTCCACCCTCCACCCA-3′) were used to amplify the target region of DNA. To amplify the 104 bp fragment of the APOA5 gene, the following conditions were used: a 2.5 minutes initial denaturation at 95℃ was followed by 40 cycles of 50 seconds at 95℃; 70 seconds at 68℃; and 70 seconds at 72℃; the final extension occurred at 72℃ for 15 minutes. The PCR reagent mix was as follows: Taq plus DNA polymerase, 0.25 µL (5 U/µL) (KBC, Terhran, Iran); 10× buffer, 2.5 µL; MgCl2 (50 mM), 0.75 µL; deoxynucleoside triphosphate (40 mM), 0.5 µL; DNA (200 ng/µL), 2 µL; 0.5 µL of the appropriate specific primer pairs (10 pmol/µL) and water in a final volume of 25 mL.
We used two different methods to select SNPs in miRNA seed regions. One approach was based on identifying the previously known miR-TS-SNPs. PolymiRTS (http://compbio.uthsc.edu/miRSNP/) [17], MirSNP (http://202.38.126.151/hmdd/mirsnp/search/) and miRNASNP (http://www.bioguo.org/miRNASNP/) are three examples of databases using this approach. In the other approach, we first predicted miRNA target sites within the 3’UTRs of the APOA5 gene by using the miRanda (http://www.microrna.org/microrna/getMirnaForm.do) [18] and TargetScan (http://www.targetscan.org/vert_50/) [19]. Then, a list of all SNPs occurring in the 3’UTRs of the APOA5 gene was obtained from the National Center for Biotechnology Information (NCBI) dbSNP database (http://www.ncbi.nlm.nih.gov/). The overlapping SNPs between these two steps were selected for final genotyping (Table 1).
Agreement of genotype frequencies with Hardy-Weinberg equilibrium (HWE) was determined using a modified chi-square test. Allele frequencies were estimated by gene counting, and the differences between case/control statuses were examined by a chi-square test. Owing to a limited number of individuals who were homozygous for the minor alleles of rs619054 and rs34089864, the genotypes for these two polymorphisms were grouped as T-carriers (TC+TT) and non-carriers (CC). One-way analysis of variance was used to compare mean values of quantitative traits across TC+TT/CC groups. The association between APOA5 polymorphisms and MetS was assessed by simple and multivariable adjusted odds ratio (OR) corresponding to a 95% confidence interval according to the logistic regression analysis. Multivariate regression analysis models were adjusted for age and gender. Constructed haplotypes, haplotype tags, and standardized linkage disequilibrium (LD), D’, and correlation, r2, and coefficients for all pairs of polymorphisms were assessed with Haploview software version 4.2 (Broad Institute, Cambridge, MA, USA) [20]. Furthermore, we used a four-gamete rule to define the haplotype blocks. Haplotype frequencies were inferred using the expectation-maximization by using HAPSTAT version 3.0 software (University of North Carolina, Chapel Hill, NC, USA; http://www.bios.unc.edu/*lin/hapstat) [2122]. Using the chi-square test, the potential difference between haplotype frequencies of case and control subjects was determined. Statistical analyses were performed using the SPSS version 16.0 statistical software (SPSS Inc., Chicago, IL, USA). All data were presented in the format of mean±standard error of the mean, and P values lower than 0.05 were considered to be statistically significant.
Table 2 presents genotyping results and allele frequencies of the two analyzed APOA5 variants (rs619054 and rs34089864). No significant difference was observed between the MetS and control groups for both genotype and allele frequencies of rs619054 and rs34089864. Moreover, none of the variants violated the HWE (data not shown).
Tables 3 and 4 represent stratification of the biochemical laboratory parameters of the control and MetS subjects according to different genotype groups of the analyzed variants. For rs619054, except for a borderline significant difference observed for insulin levels between genotype groups of control subjects, there was no statistically significant difference between other biochemical factors. For rs34089864, the HDL-C level of the CC group was significantly lower than that in the heterozygotes harbouring cytidine and thymine bases (CT) control group. Moreover, a borderline significant difference in the insulin levels between genotype groups was shown in both case and control individuals.
The ORs were evaluated for allelic and genotypic models (Table 5). Variant rs619054 demonstrated inconsequential elevated risks for MetS in both models, but both ORs were not statistically significant. Having a C allele or CC genotype in the position of rs34089864 conferred an increased risk for MetS; however, again the risk was not statistically significant.
Haplotypes were inferred using all five studied variants of the APOA5 gene (rs619054, rs369740771, rs376976905, rs1927 08363, and rs34089864). Rs619054 and rs34089864 were in LD with each other (D’=1, r2=0.02, LOD [logarithm of the odds to the base 10]=0.37). Three haplotypes were identified: H1, H2, and H3. These haplotypes are tagged by the genotypes of rs619054 and rs34089864 variants. Fig. 1 presents haplotype frequencies. Further analyses showed no significant difference between haplotype frequencies according to the disease status.
Based on stratification for different haplotypes (Table 6), MetS individuals with H3 showed significant higher levels of TG compared to other haplotypes (P=0.04). A significant increase in the HDL-C levels was also observed in controls with H3 compared to the other haplotypes (P=0.01).
Table 7 shows the association between different constructed haplotypes and risk of MetS. Although not statistically significant, individuals with H3 had the most increased risk for MetS compared to those harboring H2 before and after adjustment for age and gender.
Elevated TG and lowered HDL-C levels are the major factors among patients with MetS in Iranian children and adolescents [15]. Genetic variants of APOA5 gene have been associated with the levels of TG and HDL-C in various ethnicities [23242526]. In the current study, we have sequenced the 3’UTR of APOA5 gene in groups of 63 patients and 59 healthy controls to find potential associations with ATP III-defined MetS, and its major components in the Iranian children and adolescent population. The results showed that the alleles and genotype frequency of the two common studied SNPs (rs619054 and rs34089864) were not significantly varied in the two groups.
In one case-control study population in China, Kao et al. [27] demonstrated that c.1177C>T (rs34089864) genotype and allele frequencies did not differ between the hypertriglyceridemic and control groups. In agreement with results obtained by Kao et al. [27], we also found no significant associations between serum TG levels and c.1177C>T genotypes.
Our results demonstrated that CC genotype of variant rs34089864 was associated with lower HDL-C levels in MetS individuals. Conversely, the data showed no association between the rs619054 SNP and MetS components. In one case-control population, Pare et al. [12] showed that the rs619054 was not associated with TG, LDL-C, and HDL-C levels. Using joint analysis of the family and case-control populations, they could find a significant association with TG, but neither with LDL-C nor with HDL-C levels. Therefore, in one large population study, Grallert et al. [28] demonstrated that no significant associations were observed between rs619054, HDL-C, and LDL-C, or total cholesterol levels. In contrast, Jorgensen et al. [11] reported that CC genotype of rs619054 SNP was associated with both HDL-C and TG, but not with LDL-C levels. Therefore, the observed discrepancies between the studies may be attributed to the different race/ethnicity and age distributions of the studied populations.
Consistent with the results previously reported by Grallert et al. [28] who showed that rs619054 was not significantly associated with MetS, we found no significant association between any of the two studied major variants and MetS risk.
We also performed a haplotype analysis constructed using five SNPs in the 3’UTR of APOA5 gene. Similar to what found in Copenhagen populations [11], rs619054 and rs34089864 were in LD in our analysis.
Among the three constructed haplotypes, H1, H2, and H3, we found a significant association between H3 haplotype, increased HDL-C and TG levels, but not with MetS risk. In the same vein, H1 haplotype (harboring the two major alleles of rs619054 and rs34089864) was linked to decreased HDL-C levels. Our data are in consistent with those obtained by Pullinger et al. [29] where they constructed haplotypes with 11 APOA5 gene SNPs, including rs619054 and rs34089864 in the carriers of the c.553T polymorphism and reported that a haplotype containing the major alleles of variant rs619054 and rs34089864 was significantly associated with lowered levels of HDL-C.
In conclusion, our study explored the association between haplotypes of five SNPs located at APOA5 3’UTR and MetS entity as well as its major components risk. Our results showed that the constructed haplotypes might not be a major risk factor for the MetS entity in Iranian children and adolescents. However, our results demonstrated that the rs34089864 variant and H1 haplotype (harboring the two major alleles of rs619054 and rs34089864) were associated with HDL-C levels.
ACKNOWLEDGMENTS
This study was supported in part by a research grant from the Pediatric Inherited Diseases Research Center, Research Institute for Primordial Prevention of Non-communicable Disease, Isfahan University of Medical Sciences, Isfahan, Iran. The authors thank Prof. Mahin Hashemipour, Ms. Ghazaleh Fatemi, and Ms. Hajar Miranzadeh for their great help in the sample collection and experiments optimization steps.
References
1. Pollex RL, Hegele RA. Genetic determinants of the metabolic syndrome. Nat Clin Pract Cardiovasc Med. 2006; 3:482–489. PMID: 16932765.
2. Roche HM, Phillips C, Gibney MJ. The metabolic syndrome: the crossroads of diet and genetics. Proc Nutr Soc. 2005; 64:371–377. PMID: 16048671.
3. Chien KL, Hsu HC, Chen YC, Su TC, Lee YT, Chen MF. Association between sequence variant of c.553 G > T in the apolipoprotein A5 gene and metabolic syndrome, insulin resistance, and carotid atherosclerosis. Transl Res. 2009; 154:133–141. PMID: 19665689.
4. Nilsson SK, Heeren J, Olivecrona G, Merkel M. Apolipoprotein A-V: a potent triglyceride reducer. Atherosclerosis. 2011; 219:15–21. PMID: 21831376.
5. Zhang Z, Peng B, Gong RR, Gao LB, Du J, Fang DZ, Song YY, Li YH, Ou GJ. Apolipoprotein A5 polymorphisms and risk of coronary artery disease: a meta-analysis. Biosci Trends. 2011; 5:165–172. PMID: 21914952.
6. Alexander R, Lodish H, Sun L. MicroRNAs in adipogenesis and as therapeutic targets for obesity. Expert Opin Ther Targets. 2011; 15:623–636. PMID: 21355787.
7. Heneghan HM, Miller N, Kerin MJ. Role of microRNAs in obesity and the metabolic syndrome. Obes Rev. 2010; 11:354–361. PMID: 19793375.
8. Hu Z, Bruno AE. The influence of 3’UTRs on microRNA function inferred from human SNP data. Comp Funct Genomics. 2011; 2011:910769. PMID: 22110399.
9. Ye Q, Zhao X, Xu K, Li Q, Cheng J, Gao Y, Du J, Shi H, Zhou L. Polymorphisms in lipid metabolism related miRNA binding sites and risk of metabolic syndrome. Gene. 2013; 528:132–138. PMID: 23911300.
10. Ziebarth JD, Bhattacharya A, Chen A, Cui Y. PolymiRTS Database 2.0: linking polymorphisms in microRNA target sites with human diseases and complex traits. Nucleic Acids Res. 2012; 40:D216–D221. PMID: 22080514.
11. Jorgensen AB, Frikke-Schmidt R, West AS, Grande P, Nordestgaard BG, Tybjaerg-Hansen A. Genetically elevated non-fasting triglycerides and calculated remnant cholesterol as causal risk factors for myocardial infarction. Eur Heart J. 2013; 34:1826–1833. PMID: 23248205.
12. Pare G, Serre D, Brisson D, Anand SS, Montpetit A, Tremblay G, Engert JC, Hudson TJ, Gaudet D. Genetic analysis of 103 candidate genes for coronary artery disease and associated phenotypes in a founder population reveals a new association between endothelin-1 and high-density lipoprotein cholesterol. Am J Hum Genet. 2007; 80:673–682. PMID: 17357073.
13. Silander K, Alanne M, Kristiansson K, Saarela O, Ripatti S, Auro K, Karvanen J, Kulathinal S, Niemela M, Ellonen P, Vartiainen E, Jousilahti P, Saarela J, Kuulasmaa K, Evans A, Perola M, Salomaa V, Peltonen L. Gender differences in genetic risk profiles for cardiovascular disease. PLoS One. 2008; 3:e3615. PMID: 18974842.
14. Salehi S, Emadi-Baygi M, Rezaei M, Kelishadi R, Nikpour P. Identification of a new single-nucleotide polymorphism within the apolipoprotein A5 gene, which is associated with metabolic syndrome. Adv Biomed Res. 2017; 6:24. PMID: 28401071.
15. Kelishadi R, Ardalan G, Gheiratmand R, Adeli K, Delavari A, Majdzadeh R. Caspian Study Group. Paediatric metabolic syndrome and associated anthropometric indices: the CASPIAN Study. Acta Paediatr. 2006; 95:1625–1634. PMID: 17129973.
16. National High. Update on the 1987 task force report on high blood pressure in children and adolescents: a working group report from the national high blood pressure education program. Pediatrics. 1996; 98(4 Pt 1):649–658. PMID: 8885941.
17. Bao L, Zhou M, Wu L, Lu L, Goldowitz D, Williams RW, Cui Y. PolymiRTS Database: linking polymorphisms in microRNA target sites with complex traits. Nucleic Acids Res. 2007; 35:D51–D54. PMID: 17099235.
18. Betel D, Wilson M, Gabow A, Marks DS, Sander C. The microRNA.org resource: targets and expression. Nucleic Acids Res. 2008; 36:D149–D153. PMID: 18158296.
19. Lewis BP, Burge CB, Bartel DP. Conserved seed pairing, often flanked by adenosines, indicates that thousands of human genes are microRNA targets. Cell. 2005; 120:15–20. PMID: 15652477.
20. Barrett JC, Fry B, Maller J, Daly MJ. Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005; 21:263–265. PMID: 15297300.
21. Lin DY, Zeng D, Millikan R. Maximum likelihood estimation of haplotype effects and haplotype-environment interactions in association studies. Genet Epidemiol. 2005; 29:299–312. PMID: 16240443.
22. Lin DY, Zeng D. Likelihood-based inference on haplotype effects in genetic association studies. J Am Stat Assoc. 2006; 101:89–104.
23. Dorfmeister B, Cooper JA, Stephens JW, Ireland H, Hurel SJ, Humphries SE, Talmud PJ. The effect of APOA5 and APOC3 variants on lipid parameters in European Whites, Indian Asians and Afro-Caribbeans with type 2 diabetes. Biochim Biophys Acta. 2007; 1772:355–363. PMID: 17197160.
24. Lai CQ, Demissie S, Cupples LA, Zhu Y, Adiconis X, Parnell LD, Corella D, Ordovas JM. Influence of the APOA5 locus on plasma triglyceride, lipoprotein subclasses, and CVD risk in the Framingham Heart Study. J Lipid Res. 2004; 45:2096–2105. PMID: 15342688.
25. Martinelli N, Trabetti E, Bassi A, Girelli D, Friso S, Pizzolo F, Sandri M, Malerba G, Pignatti PF, Corrocher R, Olivieri O. The −1131 T>C and S19W APOA5 gene polymorphisms are associated with high levels of triglycerides and apolipoprotein C-III, but not with coronary artery disease: an angiographic study. Atherosclerosis. 2007; 191:409–417. PMID: 16682041.
26. Nabika T, Nasreen S, Kobayashi S, Masuda J. The genetic effect of the apoprotein AV gene on the serum triglyceride level in Japanese. Atherosclerosis. 2002; 165:201–204. PMID: 12417270.
27. Kao JT, Wen HC, Chien KL, Hsu HC, Lin SW. A novel genetic variant in the apolipoprotein A5 gene is associated with hypertriglyceridemia. Hum Mol Genet. 2003; 12:2533–2539. PMID: 12915450.
28. Grallert H, Sedlmeier EM, Huth C, Kolz M, Heid IM, Meisinger C, Herder C, Strassburger K, Gehringer A, Haak M, Giani G, Kronenberg F, Wichmann HE, Adamski J, Paulweber B, Illig T, Rathmann W. APOA5 variants and metabolic syndrome in Caucasians. J Lipid Res. 2007; 48:2614–2621. PMID: 17768309.
29. Pullinger CR, Aouizerat BE, Movsesyan I, Durlach V, Sijbrands EJ, Nakajima K, Poon A, Dallinga-Thie GM, Hattori H, Green LL, Kwok PY, Havel RJ, Frost PH, Malloy MJ, Kane JP. An apolipoprotein A-V gene SNP is associated with marked hypertriglyceridemia among Asian-American patients. J Lipid Res. 2008; 49:1846–1854. PMID: 18441017.
Table 1
Table 2
Table 3
Table 4
Values are presented as mean±standard error of the mean.
APOA5, apolipoprotein A5; CC, homozygotes harbouring cytidine bases and labeled as non-carriers; CT, heterozygotes harbouring cytidine and thymine bases; HDL-C, high density lipoprotein cholesterol; LDL-C, low density lipoprotein cholesterol.
aStatistically significant.