Association of ANGPTL8 (Betatrophin) Gene Variants with Components of Metabolic Syndrome in Arab Adults

Angiopoietin-like protein 8 (ANGPTL8) has a role in lipid metabolism, beta-cell proliferation and diabetes progression, however, the association between different variants in the ANGPTL8 gene and metabolic syndrome (MetS) components has not been studied widely especially in Arab ethnic groups. In this study, the associations of ANGPTL8 variants on MetS risk in Saudi Arab adults were investigated. A total of 905 unrelated Saudi adults (580 healthy controls and 325 MetS) were included. MetS was screened based on the International Diabetes Federation (IDF) criteria. The genotype and allele frequency distribution of rs737337 (T/C) and rs2278426 (C/T) polymorphism in ANGPTL8 gene was studied. Participants with MetS were significantly older, had higher BMI, and rs737337 polymorphism frequency was significantly lower than in control. Furthermore, the TC + CC genotype and C allele of rs737337 (T/C) was associated with decreased risk of hypercholesterolemia and hyperglycemia [odds ratio (OR) 0.61, 95%CI 0.40–0.93, p = 0.016 and OR 0.58, 0.39–0.86, p = 0.007 respectively for hypercholesterolemia; and OR 0.66, 0.45–0.97, p = 0.032 and OR 0.65, 0.46–0.93; p = 0.016 respectively for hyperglycemia]. Similarly, CT, CT + TT genotype and T allele of rs2278426 (C/T) were associated with decreased risk of hyperglycemia (p < 0.05). In conclusion, the study suggests that the gene variants in SNPs rs 737337 (T/C) and rs 2278426 (C/T) are associated with lower risk of hypercholesterolemia and hyperglycemia. These findings supplement the growing literature supporting the role of ANGPTL8 in lipid and glucose metabolism.

Metabolic syndrome is a low-grade inflammation and oxidative stress state represented by a cluster of disorders including abdominal obesity, dyslipidemia, hypertension and insulin resistance 1 . This condition predisposes itself to vascular harm and subsequent manifestation of complications like cardiovascular diseases, type 2 diabetes mellitus (T2DM) 2 . MetS has a huge impact on the global health outcome, and the epidemic proportions of its prevalence 3 demand research to reveal novel predictive factors showing association with MetS incidence and progression.
Angiopoietin-like protein 8 (ANGPTL8) or betatrophin, first identified in 2012, is a novel member of a group of secretory glycoproteins due to its role in lipid metabolism 4 . It is expressed abundantly in human liver, whereas in mice, its expression was shown to be high in liver, white and brown adipose tissues 5 . In humans, it is expressed predominantly in the liver 6 , suggesting that it might function as a hepatokine. Because of the controversy related to its function as a β-cell proliferator, many researchers in the recent past have tried to investigate the effects of betatrophin in beta-cell proliferation and diabetes expansion 7,8 . Nevertheless, while some suggest no role of betatrophin in beta-cell proliferation 9 , yet there are evidences in the literature which suggest that the factors related to metabolic diseases such as age, gender, waist circumference, total/HDL cholesterol ratio and renal function Genotype and Allele frequency distribution of the two SNP's. The genotype and allele frequency distribution of SNP's rs737337 (T/C) and rs2278426 (C/T) polymorphism in ANGPTL8 gene are presented in Table 2. The T-allele frequency of rs737337 polymorphism was significantly higher in controls than MetS participants (p = 0.02), however, the genotype/allele frequency distribution of rs2278426(C/T) was comparable in both the study groups (p > 0.05).
The two SNPs were in complete linkage equilibrium with linkage disequilibrium calculated as 0 (D was calculated as described in supplementary file SF1).

Association of ANGPTL8 gene variants in the two study SNPs with anthropometric and metabolic parameters.
To determine differences in the biochemical parameters based on ANGPTL8 gene variants, we analyzed these parameters in genotypes of rs737337 (T/C) and rs2278426(C/T) and the results are presented in Table 3. After controlling for age and sex, participants with TT genotype of rs737337 SNP showed significantly higher total cholesterol (p < 0.01) and LDL-cholesterol levels (p < 0.01) than those with TC and CC genotypes. Similarly, those with CC genotype of rs2278426 SNP showed significantly higher fasting glucose levels (p = 0.029), total-cholesterol levels (p = 0.004) and LDL-cholesterol levels (p < 0.001) than those with CT and TT genotypes. (>5.7 mmol/l), high glucose (>5.6 mmol/l) and MetS (IDF criteria) and analyzed the effect of ANGPTL8 gene variants on the above classified variables. A logistic regression analysis of the genotypes/alleles in the SNP rs737337 and SNP rs2278426 representing risk of having the different components of MetS are presented in Tables 4 and 5 respectively. The TC + CC genotype and C allele of rs737337 were associated with decreased risk of hypercholesterolemia in the univariate model as well as after adjustment with age and gender [Odd ratio (OR) 0.61, 95%CI: 0.40-0.93, p = 0.02 and OR 0.58, 95% CI: 0.39-0.86, p = 0.007 respectively for univariate; and OR 0.65 (0.42-0.99), p < 0.05 and OR 0.62 (0.41-0.93), p = 0.02 respectively after adjustment]. In addition, they were also associated with decreased risk of hyperglycemia in both models (OR 0.66, 95% CI: 0.45-0.97, p = 0.03 and OR 0.65, 95% CI: 0.46-0.93, p = 0.02 respectively for TC + CC and C alleles in univariate model; and OR 0.67 (0.46-0.99), p < 0.05 and 0.69 (0.48-0.98), p = 0.04 respectively in the adjusted model). The c-allele of rs737337 was also found to be associated with decreased risk of full MetS in the univariate model, however lost statistical significance in the adjusted model (OR 0.69, 95% CI: 0.49-0.99,p = 0.04 in the univariate model and OR 0.75 (0.52-1.08,p = 0.12 in the adjusted model). (Table 4).

Discussion
The low-grade inflammatory state of MetS and its components like dyslipidemia, hyperglycemia, hypertension and obesity has a more severe outcome of 5-fold and 2-fold increase respectively in T2DM and cardiovascular diseases 19,20 . To address the public issue posed by MetS and associated complications, researchers in the recent past have been focusing on the preventive and therapeutic strategies to combat this menace. ANGPTL8 seems to be one such strategy and recent studies have demonstrated that the circulating levels of ANGPTL8 may be important in the course of development of MetS and its components 21 . Research studies are limited to the association of the circulating levels of ANGPTL8 and MetS state, and the investigations on the associations of ANGPTL8 gene variants with MetS and its individual components are limited particularly in Saudi Arabia.
In this study, we examined whether the polymorphism of ANGPTL8 SNPs rs2278426 (C/T) and rs737337 (T/C) play a role in the increased risk of MetS in Saudi Arab adults. As rs2278426 and rs737337 variants have been found to be associated with HDL-C levels and lipoprotein metabolism, which provide novel support for functional consequences of the both variants. To establish the link between MetS and SNPs of ANGPTL8 rs2278426 and rs737337, we identified the frequency and distribution of SNPs in the studied participants and evaluated their association with MetS phenotypes. Our findings demonstrated that variant C allele in rs737337 is associated with decreased risk of hypercholesterolemia and hyperglycemia; and the variant T allele in rs2278426 is associated with decreased risk of hyperglycemia and these associations remained statistically significant even after adjustment with age and sex.. Interestingly, our study also pointed out that allele of rs2278426 was associated with increased risk of hypertension even after adjustment.
ANGPTL8 or betatrophin is a novel liver and adipose tissue secreted protein that is nutritionally regulated. In humans, betatrophin is liver specific and its expression has been found to be relatively higher than in other tissues, while in mouse it is enhanced in both liver and adipose tissue 16 . ANGPTL8 is known to play a key role in lipid metabolism 22 , and our study showed that ANGPTL8 gene variants were significantly associated with total cholesterol and TG levels. Also, the odds for having hypertension significantly increased for the T allele of SNPs rs2278426 (C/T). This is in line with previous studies 23 which showed ANGPTL8 levels higher in hypertensive subjects compared to non-hypertensive subjects. This association with hypertension may be explained by a strong correlation of ANGPTL8 with high sensitive C-reactive protein (HsCRP) found in some studies 24 . HsCRP with oxidized LDL have been shown to form complexes that promote atherosclerosis in diabetic mice leading to vascular inflammation and development of hypertension 25,26 . The association between betatrophin and hypertension indicates that its plasma level might be involved in the risk of cardiovascular events and could be regarded as a predictive biomarker for early diagnosis and management. It will be of interest though to validate and strengthen our findings, by carrying out the same study in a larger size population and or in two different Saudi populations to deepen our understanding of the association of ANGPTL8 SNPs rs2278426 and rs737337 with MetS.
Homology of ANGPTL8 with ANGPTL3's and ANGPTL4's N-terminal domains are necessary for lipid regulation through LPL binding 27 . Moreover, our data showed that lipid levels were significantly higher in MetS subjects. Zhang et al. found that overexpression of betatrophin has been shown to increase serum triglycerides, whereas recombinant betatrophin inhibits LPL activity 16 . Furthermore, LPL inhibition in cardiac and skeletal muscle have been linked with ANGPTL8 28 . In addition, some studies reported that ANGPTL8 co-expression with ANGPTL3 further increases serum TG showing a positive association 29,30 . Among various ethnic groups, significant differences in the allele frequencies of ANGPTL8 gene variants and their association with concentration of blood lipids have been reported. In Hispanic and African Americans, for example, the T allele of rs2278426 was reported to be associated with lower levels of LDL-C, while in European Americans this association was absent 14 . Similarly, a study showed a significant inverse correlation between ANGPTL8 and TG in patients with T2DM 31 , while other reports showed no association 32   Continued studies are needed to uncover the pathophysiological pathways that involves betatrophin in the in the development of MetS. In addition to its role in lipoprotein metabolism, ANGPTL8 is a potential target for β-cell regenerative therapy in diabetes due to its ability to induce β-cell proliferation 33 , suggesting that decreases in ANGPTL8 levels or function may worsen glucose tolerance. Previous studies demonstrated that betatrophin was up-regulated in T2DM individuals, which was assumed to due to increase in insulin production through increasing β-cell proliferation 34,35 . Although earlier reports stressed that ANGPTL8 expression gets induced by insulin, however, there are reports that suggest no association between ANGPTL8 and β-cell proliferation as mice lacking ANGPTL8 showed normal glucose metabolism under insulin resistance conditions 9 . In agreement with the associations found in this study, Abu-Farha et al. in his study while examining the changes in circulating ANGPTL8 levels in a large cohort of 1049 non-diabetic people and 556 people with T2DM found positive association with fasting glucose and DM duration 36 . In this study, the C-allele and T-allele of SNPs rs737337 (T/C) and rs rs2278426 (C/T) are associated with decreased risk of hyperglycemia. The exact mechanism by which sequence variation exerts its effect has not been established, it may be possible that the ANGPTL8 protein structure could be affected due to the variation in the charge of the amino acid residue from arginine to tryptophan. Research findings still suggest a role for ANGPTL8 in glucose metabolism apart from β-cell proliferation, however, to confirm that it will require further research into the physiological factors that regulate ANGPTL8 expression.
In the present study, older participants exhibited the typical characteristics of MetS including central obesity. They also showed high blood pressure, elevated glucose and triglyceride levels and decreased HDL-cholesterol levels. Also, the variant C-allele of SNP rs737337 (T/C) showed decreased risk of the full MetS. This is in agreement with studies showing increased expression of betatrophin in T2DM and/or obese subjects 32,34 . Furthermore, there is accumulating evidence that showed positive associations between circulating betatrophin levels, insulin resistance and T2DM 37,38 . The discrepancy in the study 39 that showed that ANGPTL8/betatrophin concentrations was decreased in human obese subjects as well as subjects with T2DM was probably because of the kit used to measure betatrophin. Other studies reported non-significant differences in betatrophin levels in subjects with or without DM 13  The data is represented by the odds ratio and its 95% confidence interval. The data was adjusted for age and gender in the adjusted model. size, BMI, age and ethnic group; and different ELISA kits employed 34 , however, further investigations are required to elucidate these discrepancies. Despite these interesting findings reported in this study, we acknowledge certain limitations. First of all the investigation was done on a specific population comprising of Saudi Arabian adults and the findings should be used with caution when extrapolated to other ethnicities. Secondly, the study focused on the gene variants in the two SNPs rs737337 and rs 2278426 with the risk of full MetS and its individual components and lacked the data on the circulating levels of betatrophin and insulin in the studied subjects. The association of circulating levels of betatrophin and/or insulin resistance with MetS, thus, cannot be concluded in Saudi adults from this study and probably needs a separate investigation. Also, because of the difficulty in assuming the prevalence of the minor alleles of a SNP in the two study groups particularly when there is a limited literature available for the concerned population, it is hard to estimate a sufficient sample size for a genotype association study. Keeping this in mind, the authors recommend to validate the results observed here in a bigger sample size and if possible in two different populations. Nonetheless, with its limitations as mentioned, the findings in this study are still of interest not only in terms of the association of gene variants of a novel glycoprotein, ANGPTL8 with the disorders (MetS and its components) that has global health outcome but also as a stepping stone for future such studies.

conclusion
In conclusion, the study investigated the gene variants in SNPs rs 737337 (T/C) and rs 2278426 (C/T) in association with the risk of MetS and its individual components in Saudi adults and found that these variants are associated with lower risk of hypercholesterolemia and hyperglycemia as components of MetS. These findings supplement the growing literature supporting the role of ANGPTL8 in lipid and glucose metabolism, the disorder of which contributes to the development of MetS, however at the same time recommends that the pathophysiological upregulation of betatrophin in people with MetS needs further investigation to underline the corresponding molecular mechanisms.

Research Design and Methods
Subjects. Sample size estimation for this genetic association study was done by using G*power 3.1.9.4 software. For a two-tailed Fischer exact test in two independent groups with α = 0.05; power 0.80; proportion 1 as 0.35 and allocation ratio of 0.54 (since prevalence of metabolic syndrome in Saudi Arabia is 35% as per the literature); and assuming a 10% difference of risk alleles in cases and controls we expected another proportion as 0.45. When these values were entered, a sample size of 847 was calculated. Based on this, a total of nine hundred and five unrelated Saudi adults [580 healthy controls (296 males and 284 females) and 325 MetS (145 males and 185 females)] were randomly selected from the database of Biomarker Screening conducted by Chair for Biomarker of Chronic Diseases (CBCD) in King Saud University (KSU), Riyadh, Saudi Arabia, a capital-wide epidemiologic study for Saudis 41 . In brief, using a cluster sampling strategy, participants were recruited from different regions in and around the Riyadh region. An overnight fasting blood sample was taken from the participants at the associated primary health care center (PHCC) where anthropometric measurements and data related to general demographic information; past and present medical history etc. was also recorded. Before including in this study, a written informed consent was collected from all participants. The study was approved by Ethics Committee at the College of Science, KSU, Riyadh, Saudi Arabia.
Anthropometry, blood collection and biochemical analysis. Anthropometric data was taken from the database. This comprised of height (centimeters), weight (kg), waist and hip circumference (centimeters), and blood pressure (in mmHg). To calculate BMI, we used the formula weight in kg/(height in meter) 2 . Fasting whole blood and serum samples collected at the recruitment were transported at proper conditions to the biobank facility of CBCD, aliquoted there and kept at suitable storage conditions for future analysis. Glucose, total cholesterol, HDL-cholesterol and triglycerides as biochemical analysis were measured using routine commercially available kits on an automated biochemical analyzer (Konelab, Espoo, Finland). The procedure of biochemical analysis was already published in detail 42 . Whole blood samples were used for genetic analysis.  www.nature.com/scientificreports www.nature.com/scientificreports/ were employed in a real time PCR to evaluate the allelic discrimination in the two tagging SNPs (rs737337 and rs2278426). The TaqMan assay IDs (primers and probes) are given in Table 6. 10 µL of volume containing 1X TaqMan genotyping Master Mix (Applied Biosystems), 1X mix of unlabeled PCR primers and TaqMan probes, and 30 ng of template DNA were mixed in 96-well PCR plates for amplification and a Bio-Rad CFX96 Real-Time PCR Detection System (Bio-Rad, Milan, Italy) was utilized for its detection. Thermal cycling utilized a denaturation step of 10 min at 95 °C, which was followed by 45 cycles of 15 s at 95 °C and 90 s at 60 °C. The software used for allelic discrimination analysis was from Bio-Rad (CFX Manager Software, Version 1.6). Dedicated PCR pipettes and reagents were only employed for all PCR reactions. A validation step of re-genotyping of 40 random samples was done and the results were found to be reproducible. Statistical analysis. SPSS (version 22.0, IBM, NY, USA) was used to analyze the data. Continuous Gaussian variables were shown as mean ± standard deviation (SD) and continuous non-Gaussian variables were shown as median (Q1-Q3). Genotype distributions between the MetS and control subjects were reported as allele frequency (%) and the difference was computed using Chi-square test of distributions. Odds ratios (ORs) and 95% confidence intervals (CIs) depicting the risk of metabolic syndrome and its individual components vis a vis various alleles in the two studied SNP's were calculated using multinomial logistic regression with genotypes as the factor; sex and age as covariates; and most common genotypes as reference. P-value <0.05 was considered significant in all tests (two-tailed).

Genetic analysis. Genomic
Statement of Informed consent. Written informed consent was obtained from each participant in this study.
Details of Ethics Approval. All the procedures followed in this study were in accordance with the ethical standards of the responsible committee on human experimentation (institutional and national) and with the Helsinki Declaration of 1975, as revised in 2008. The study was approved by the Ethics Committee of College of Science, King Saud University.

Data availability
The dataset analysed for getting the results observed in this study is available from the corresponding author on reasonable request.