Next Article in Journal
Next-Generation Sequencing Identifies Novel PMPCA Variants in Patients with Late-Onset Dominant Optic Atrophy
Previous Article in Journal
Immune-Related Genomic Schizophrenic Subtyping Identified in DLPFC Transcriptome
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

GALNT2 rs4846914 SNP Is Associated with Obesity, Atherogenic Lipid Traits, and ANGPTL3 Plasma Level

1
Pharmacology and Therapeutics Department, Faculty of Pharmacy, Kuwait University, Kuwait City 13110, Kuwait
2
Department of Genetics and Bioinformatics, Dasman Diabetes Institute, Dasman 15462, Kuwait
3
Department of Biochemistry and Molecular Biology, Dasman Diabetes Institute, Dasman 15462, Kuwait
4
Special Services Facility, Dasman Diabetes Institute, Dasman 15462, Kuwait
5
Advanced Genomic Technologies Laboratory, McGill University, Montreal, QC H3A 0G4, Canada
6
Department of Medicine, Faculty of Medicine, Kuwait University, Kuwait City 13110, Kuwait
7
Research Division, Dasman Diabetes Institute, Dasman 15462, Kuwait
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Genes 2022, 13(7), 1201; https://doi.org/10.3390/genes13071201
Submission received: 30 May 2022 / Revised: 29 June 2022 / Accepted: 29 June 2022 / Published: 4 July 2022
(This article belongs to the Section Molecular Genetics and Genomics)

Abstract

:
N-Acetylgalactosaminyltransferase 2 (GALNT2) is associated with serum lipid levels, insulin resistance, and adipogenesis. Additionally, angiopoietin-like (ANGPTL) proteins have emerged as regulators of lipoprotein lipase and lipid metabolism. In this study, we evaluated the association between GALNT2 rs4846914 variant, known for its association with lipid levels in European cohorts, with plasma levels of ANGPTL proteins, apolipoproteins, lipids, and obesity traits in individuals of Arab ethnicity. GALNT2 rs4846914 was genotyped in a cohort of 278 Arab individuals from Kuwait. Plasma levels of ANGPTL3 and ANGPTL8 were measured by ELISA and apolipoproteins by Luminex multiplexing assay. Allele-based association tests were performed with Bonferroni-corrected p-value thresholds. The GALNT2 rs4846914_G allele was associated with increased ANGPTL3 (p-values ≤ 0.05) but not with ANGPTL8 plasma levels. The allele was associated significantly with higher BMI and weight (p-values < 0.003), increased ApoC1 levels (p-values ≤ 0.006), and reduced HDL levels (p-values ≤ 0.05). Individuals carrying the GG genotype showed significantly decreased HDL and increased BMI, WC, ApoC1, and TG. Interactions exist between (AG+GG) genotypes and measures of percentage body fat, ApoA1A, ApoC1, and ApoB48-mediated HDL levels. GALNT2 is confirmed further as a potential link connecting lipid metabolism and obesity and has the potential to be a drug target for treating obesity and dyslipidemia.

1. Introduction

Obesity incidence in Arabian nations has increased at an alarming rate over the last several decades, apparently due to a dramatic lifestyle shift toward a more sedentary and nutrient-dense lifestyle [1,2,3]. The dramatic rise in obesity levels has led to a significant increase in the incidence of metabolic syndrome, which is characterized by visceral obesity, insulin resistance, raised blood pressure, elevated triglycerides, and low levels of high-density lipoprotein cholesterol (HDL-C). A low HDL-C level is one of the significant indicators of cardiovascular disease in metabolic syndrome. Additionally, low HDL-C is a component of the “atherogenic dyslipidemia” condition, which also includes high triglycerides (TG) [4]. Disorders such as atherogenic dyslipidemia, characterized by high TG and low HDL levels, are often seen in patients with insulin resistance, obesity, and type 2 diabetes (T2D) and play major roles in shaping the risk of cardiovascular disease.
Despite the higher incidence of obesity and other metabolic diseases in the Arab region, there are no firmly established genetic variants for these metabolic diseases in the Arab population. Recent genome-wide association (GWA) investigations have discovered novel genetic determinants for a variety of complicated quantitative characteristics, including dyslipidemia [1,5,6]. The polypeptide N-acetylgalactosaminyltransferase 2 gene (GALNT2) is a risk locus discovered for dyslipidemia [7,8]. GALNT2 is responsible for the attachment of N-acetyl galactosamine to a serine and threonine residue in proteins. This process of O-glycosylation is a major posttranslational modification that significantly affects protein function, with the most notable one being its ability to regulate the processing by the proprotein convertase family of enzymes [9]. For example, GALNT2 was shown to inhibit the activation of angiopoietin-like protein 3 (ANGPTL3) by proprotein convertase enzyme or furin [10]. Furin processing is thought to release the secreted N-terminal coiled-coil domain of ANGPTL3, which is a known inhibitor of both lipoprotein lipase and endothelial lipase. Other GALNT2-O-glycosylated substrates involved in lipid metabolism regulation include apolipoproteins such as ApoA1, ApoE, and ApoCIII as well hepatic lipase (LIPC) and very low-density lipoprotein receptors [9]. This suggests a possible role of GALNT2 in the regulation of plasma lipids.
Previous studies have highlighted the association of GALNT2 with lipid levels, insulin resistance, adipogenesis, and related cellular phenotypes, as briefed below. Studies both in humans and animal models have highlighted GALNT2 as a determinant of serum HDL and TG levels. GALNT2 mRNA levels are associated with serum triglycerides in humans [11]. Understanding whether different degrees of changes in GALNT2 modulate different serum lipid fractions can make GALNT2 a target for treating atherogenic dyslipidemia and related clinical events [12]. GALNT2 gene variant rs4846914 has been associated with insulin and insulin resistance depending on BMI in polycystic ovary syndrome patients [13]. It is also recognized that GALNT2 is a novel modulator of adipogenesis and related cellular phenotypes, thus becoming a potential target for tackling the epidemics of obesity and its devastating health consequences [14].
A recent meta-analysis study from our laboratory [1] observed 25 variants (with rs666718 as the lead SNP) from GALNT2 associating with TG at p-values of 10−7. All these 25 variants are also annotated in the NHGRI-EBI GWAS Catalog [15] associated with HDL in global populations. The GALNT2 rs4846914 variant, another well-studied variant in international GWA studies, is not in linkage–disequilibrium with any of the 25 variants highlighted above from our meta-analysis study. This variant has been associated with TG and HDL in global studies using cohorts of European ancestry and cohorts of trans-ethnic ancestries comprising largely people of European ancestry along with African Americans, East Asians, and South Asians [see, for example [16,17,18,19]. Given that obesity, diabetes, and lipid profiles are inter-connected, we aimed in this study to evaluate the association of this GALNT2 rs4846914 variant with metabolic traits, including the obesity traits and related biomarkers, such as Apolipoproteins and ANGPTL3, in a cohort of Arab individuals. Further, given that obesity, diabetes, and lipid profiles are inter-connected with one another, we aimed to delineate interactions between the genotype at the variant level and the measurements of traits/biomarkers relating to these processes.

2. Materials and Methods

2.1. Recruitment of Participants and Study Cohort

The study protocol was reviewed and approved by the Ethical Review Committee of Dasman Diabetes Institute as per the guidelines of the Declaration of Helsinki and of the US Federal Policy for the Protection of Human Subjects (Study number RA2010-003). The study subjects were native adult Kuwaiti individuals of Arab ethnicity. Pregnant women were excluded. The cohort comprised 278 subjects. Data on age, sex, health disorders (e.g., diabetes and hypertension), and baseline characteristics such as height, weight, waist circumference, and blood pressure were recorded for each participant upon enrolment. Furthermore, information on whether the participants take lipid-lowering or diabetes and antihypertensive medications was recorded and was subsequently used to adjust the models for genotype-trait association tests. Informed consent form was signed by every participant before participating in the study.
Measurements such as BMI, glycated hemoglobin (HbA1c) levels, and blood pressure readings were made as per international guidelines. For example: (a) Height is measured to the nearest centimeter with the participant standing upright against a wall on which is fixed a height measuring device. The head is held in the Frankfort position and the heels are held together. (b) Weight measurements are taken on a pre-calibrated electronic weighing scale that is placed on a firm flat surface. The participant is weighed dressed in light clothes, barefooted, facing forward, and standing still. Weight is recorded to the nearest 100 g. (c) The mercury type of sphygmomanometer was used to measure blood pressure. The participant is made to sit quietly with the right arm placed on the table with the palm facing upwards. Average of three readings of blood pressure is recorded. Clinical guidelines were followed to ascertain the diagnosis for diabetes.
Participants with a BMI ≥ 30 Kg/m2 were considered obese.

2.2. Blood Sample Collection and Processing

Upon confirming that participants had fasted overnight, blood samples were collected in EDTA-treated tubes. DNA was extracted using Gentra Puregene® kit (Qiagen, Valencia, CA, USA) and was quantified using Quant-iT™ PicoGreen® dsDNA Assay Kit (Life Technologies, Grand Island, NY, USA) and Epoch Microplate Spectrophotometer (BioTek Instruments, Vermont, USA). Absorbance values at 260–280 nm were checked for adherence to an optical density range of 1.8–2.1.

2.3. Estimation of HbA1c, Plasma Glucose and Lipid Parameters

Fasting blood glucose, TG, total cholesterol, LDL, and HDL levels were measured using Siemens Dimension RXL chemistry analyzer (Diamond Diagnostics, Holliston, MA, USA). Glycated hemoglobin content was measured using a Variant™ device (Bio-Rad, Hercules, CA, USA).

2.4. Estimation of Plasma Levels of Various Biomarkers

Plasma was separated from blood samples by centrifugation, was aliquoted, and was stored at −80 °C. Briefly, plasma was obtained from blood samples following the centrifugation of the blood tubes at 400× g for 10 min and was stored at −80 °C in new tubes. Any remaining cells or platelets in plasma were removed by centrifugation of ice-thawed plasma for 5 min at 10,000× g at 4 °C.
ANGPTL3, and 4 levels were measured by Human ELISA (R&D systems, Minneapolis, MN, USA, Cat# DANL30 and DY3485, respectively), as reported previously [20,21]. ANGPTL8 was measured as previously reported, [20,21,22], using ELISA kit from (EIAab Sciences, Wuhan, China, Cat# E1164H). Assays were performed according to manufacturing protocols.
Plasma levels of ApoC1, ApoA1A, ApoA2, and ApoB were measured using multiplexing assay MILLIPLEX MAP Human Apolipoprotein Magnetic Bead Panel (Bio-Rad, Hercules, CA, USA, APOMAG-62K). Assays were performed as per the manufacturer’s instructions. The samples were analyzed on the Bio-Plex 200 system (Bio-Rad, Hercules, CA, USA), and the Bio-Plex manager software was used to quantify the concentration of each analyte through the generated standard curve.

2.5. Bioelectric Impedance Measurements

Bio-impedancemetry Body Composition Analyzer IOI 353 (Jawon Medical Co., Seoul, South Korea) was used to measure body composition of the participants. Summation of body intracellular water and water outside the cell membrane defined the total body water (TBW). Soft lean mass (SLM) is defined by the addition of TBW and proteins in the body and is made up of skeletal and smooth muscle [23]. Lean body mass (LBM) is the summation of SLM and minerals. Percentage body fat (PBF) was estimated by subtracting LBM from the total body weight.

2.6. Targeted Genotyping of the Study Variant

We performed candidate SNP genotyping using the TaqMan® Genotyping Assay on ABI 7500 Real-Time PCR System from Applied Biosystems (Foster City, CA, USA). Each polymerase chain reaction sample was comprised of 10 ng of DNA, 5× FIREPol® Master Mix (Solis BioDyne, Tartu, Estonia), and 1 µL of 20× TaqMan® SNP Genotyping Assay. Thermal cycling conditions were set at 60 °C for 1 min and 95 °C for 15 min followed by 40 cycles of 95 °C for 15 s and 60 °C for 1 min. We further performed Sanger sequencing, using the BigDye™ Terminator v3.1 Cycle Sequencing on an Applied Biosystems 3730xl DNA Analyzer (Applied Biosystems, Foster City, CA, USA) for selected cases of homozygotes and heterozygotes to validate genotypes determined by candidate SNP genotyping. The dataset of genotypes at the study variant is presented in Supplementary Dataset 1.

2.7. Quality Procedures for SNP and Trait Measurements

SNP quality and statistical associations with traits were evaluated using PLINK (version 1.9) [24]. We calculated minor allele frequency (MAF) and Hardy–Weinberg equilibrium for the study variant. Any quantitative trait value < Q1-1.5 × IQR or any value > Q3 + 1.5 × IQR was considered as an outlier and was excluded from further statistical analysis. Normality of all the traits was assessed. Mean was used as measure of centrality. Traits were not subjected to any transformation.

2.8. Allele-Based Association Tests and Thresholds for Ascertaining Statistical Significance

Allele frequency differences between obese vs. non-obese and diabetic vs. non-diabetic were assessed using Fisher exact tests. Allele-based statistical association tests for the study variant with 10 quantitative traits and 7 biomarker levels were performed using linear regression adjusting for regular corrections toward age and sex. We also adjusted for diabetes medication and lipid-lowering medication. Correction for multiple testing was assessed by adjusting the p-value threshold for the number of tested traits (n = 17), which was (0.05/17 = 0.003). Interactions between genotype–phenotype associations and measures of biomarkers were performed using linear regression analyses.

3. Results

The study variant GALNT2 rs4846914 passed the tests for Hardy–Weinberg equilibrium (HWE). Clinical characteristics of the study cohort are presented in Table 1. Participants in the study cohort had a mean age of 46.25 ± 12.38 years. The ratio of males to females in the cohort was 1:1.22. Participants were mostly class I obese people with a mean body mass index (BMI) of 29.93 ± 5.17 kg/m2 and mean waist circumference (WC) of 99.36 ± 13.36 cm. Up to 43.16% were afflicted with type 2 diabetes. Mean values for HbA1c (6.31 ± 1.3%), LDL (3.13 ± 0.96 mmol/L), HDL (1.2 ± 0.32 mmol/L), total cholesterol (TC) (5.02 ± 1.09 mmol/L) and TG (1.22 ± 0.6 mmol/L) were normal or near normal. Of the 278 participants, 101 were taking T2D medication and 88 were taking lipid lowering medications.
Apart from levels of obesity traits (weight, BMI and WC), participants in our study cohort afflicted with obesity differed significantly from those without obesity in the levels of body composition traits (PBF: percentage body fat), lipid traits (TG), glycemic traits (FPG and HbA1c), angiopoietin-like proteins (ANGPTL3, ANGPTL4, and ANGPTL8), and diabetes status (Table 1).
Upon partitioning the cohort based on the genotypes of the studied variant, significant (p-value ≤ 0.05) or close to significant increases in BMI, WC, ApoC1, and TG and decreases in HDL levels were seen in individuals carrying the GG genotype homozygous for the effect allele compared to those with AA genotype (Figure 1, Figure 2, Figure 3 and Figure 4). Though differences were seen in the levels of ANGPTL3, they were not statistically significant.
Assessment for frequency difference between obese vs non-obese individuals demonstrated association of the GALNT2 rs4846914_G allele with an increased risk of obesity at an odds ratio of 1.47 (CI:1.024–2.055) with p-value 0.03 (Table 2).
As regards quantitative traits, the rs4846914_G allele was significantly associated with increased levels of BMI and weight (with Bonferroni-corrected p-values < 0.003), increased levels of ApoC1 (with highly significant p-values ≤ 0.006), higher levels of ANGPTL3 (at p-values with a trend towards significance, and lower levels of HDL (at p-values ≤ 0.05) (Table 3).
Levels of HDL were mediated by interactions between carrier genotypes (AG+GG) at rs4846914 and measures of percentage body fat (PBF), ApoA1A, ApoC1, and ApoB48 (Table 4). While levels of HDL were directly correlated with these interacting partners in individuals with reference AA genotype, the correlation was markedly inverse with PBF in individuals with carrier (AG+GG) genotype (Figure 5). These observations illustrate GALNT2 as a potential link connecting lipid metabolism and obesity.

4. Discussion

GALNT2 rs4846914 displays considerable variation (6% to 60%) in allele frequency across the continents. The frequencies of the rs4846914-A allele across the populations, as seen in the 1000 Genomes Project Phase 3 data [25] (as presented in Ensembl genome browser [26]), are: 6% in Africans, 23% in East Asians, 39% in South Asians, 54% in Ad-mixed Americans, and 60% in Europeans. Our GWAS data set on Kuwaiti population (described in [1,2,3,4] show a frequency of 46%. Global GWA studies have associated the rs4846914_G with a decrease in HDL levels and an increase in TG levels (see, for example [16,17,18,19]); such studies were performed on cohorts of European ancestry or on trans-ethnic cohorts with people of European ancestry dominating the cohort. In our study cohort, we similarly found an association of the G allele with a decrease in HDL (p-value = 0.02) and an increase in TG (though at p-value = 0.088).
Our study demonstrates the association of the GALNT2 rs4846914_G allele with an increased risk of obesity at an odds ratio of 1.47 (CI: 1.04–2.05). Our study further finds the variant allele significantly associated with increased levels of BMI, weight, ApoC1 and ANGPTL3, and with lower levels of HDL. Individuals carrying the GG genotype homozygous for the effect allele exhibited significant increases in BMI, WC, ApoC1, and TG and a decrease in HDL as compared to individuals carrying the AA genotype. In comparison with a study in China to assess the polymorphism of GALNT2 rs4846914 on lipid levels, subjects with the GG genotype in Han population exhibited higher HDL-cholesterol but lower LDL-cholesterol and ApoB levels than the subjects with the AA genotype [27].
Moreover, the study finds the levels of HDL to be mediated by interactions between the carrier genotypes (AG+GG) and measures of percentage body fat (PBF), ApoA1A, ApoC1, and ApoB48. While levels of HDL were directly correlated with these interacting biomarkers in individuals with the AA genotype, the correlation was markedly inverse with PBF in individuals with the (AG+GG) genotypes. The study points out that the G allele at the GALNT2 rs4846914 coordinates the regulation of body fat, levels of HDL, and apolipoproteins in mediating obesity. Relationships between HDL and apolipoproteins in determining the risk or severity of metabolic disorders have been discussed in the literature, for example, Richardon et al., by way of performing a study implemented through multivariable Mendelian randomization simultaneously that accounts for genetic associations with lipids and apolipoproteins, observe that changes in cholesterol or triglycerides that are not accompanied by commensurate changes in apolipoprotein B may not lead to altered risks of coronary heart disease [28].
Our observations of interactions among changes in obesity genes, apolipoproteins, and ANGPTL3 are corroborated by findings from our earlier studies as indicated below: We have earlier demonstrated that the popular obesity gene FTO variant regulates obesity traits through interaction between carrier genotypes and measures of Apo’s and ghrelin [29]. We have further demonstrated earlier that the gene expression levels of ANGPTL3 are increased in obese subjects [20].
ANGPTL3 is known to regulate lipid metabolism through their inhibitory effect on lipoprotein lipase (LPL) and endothelial lipase. Furthermore, it is well known that GALNT2 regulates ANGPTL3 such that deficiency of GALNT2 expression or knockout in hepatocytes increased ANGPTL3 cleavage and activation due to the absence of GALNT2 mediated O-glycosylation. On the other hand, overexpressing GALNT2 decreased ANGPTL3 cleavage [30]. Since individuals homozygous for the GG genotype for the GALNT2 rs4846914 variant were associated with elevated levels of TG and lower levels of HDL, we can speculate that this variant somehow results in a loss of function and reduced glycosylated activity of GALNT2, thus leading to increased ANGPTL3 cleavage and activation by furin. Through its O-glycosylated activity, GALNT2 is known to affect the function of several proteins involved in lipid metabolism, such as ApoA1, ApoE, and ApoCIII, as well as hepatic lipase (LIPC) and very low-density liproprotein receptors [9]. In our study, we have shown individuals homozygous for the GG allele of GALNT2 variant had higher ApoC1 and ApoA2 levels. ApoC1 is thought to inhibit IDL and VLDL binding and uptake via both LDL receptor and lipoprotein receptor-related protein and to downregulate the activity of lipoprotein lipase, hepatic lipase, and cholesterol ester transfer protein. Elevated ApoC1 levels were associated with carotid intima media thickness, hyperlipidemia, and aggravated coronary artery disease in both animals and humans [31,32,33]. A recent study showed an association of the GALNT2 rs4846914 variant with atherogenic index in overweight/obese women with gestational diabetes mellitus [34], thus highlighting the possible regulatory role of GALNT2 in obesity and lipid traits. Both our earlier meta-analysis study [1] and global GWA studies [6] have amply demonstrated the association of GALNT2 variants, particularly the study variant, with lipid traits.
In this presented study, the impact of the GALNT2 rs4846914_G variant on TG and HDL levels is demonstrated by the significant differences seen in their levels between the individuals carrying the homozygous reference genotypes versus homozygous effect allele genotypes. Individuals carrying GG genotype at GALNT2 rs4846914 variant as opposed to those carrying (AA+AG) genotypes showed significantly decreased HDL and increased BMI, WC, ApoC1, and TG. The study further demonstrates that the levels of HDL are mediated by interactions between the (AG+GG) genotypes and measures of percentage body fat, ApoA1A, ApoC1, and ApoB48.
By way of demonstrating the association of the study variant with obesity/lipid traits and demonstrating the regulation of lipid levels through interaction with apolipoproteins and percentage body fat, the presented study builds a strong case proposing GALNT2 as a potential link between obesity and lipid metabolism.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/genes13071201/s1, The genotype dataset generated and analyzed during the current study is available as Supplementary Dataset 1.

Author Contributions

Conceptualization, M.A.-F., J.A., F.A.-M. and T.A.T.; methodology, M.M., F.A.-K., I.A. and P.C.; formal analysis, M.A.-F., P.H. and T.A.T.; investigation; writing—original draft preparation, P.H., T.A.T., M.A.-F. and M.Q.; writing—review and editing, M.A., T.A.T., M.A.-F., P.H., A.A.S., J.A. and F.A.-M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received institutional funding from the Kuwait Foundation for the Advancement of Sciences (KFAS) for the project RA-2016-026.

Institutional Review Board Statement

This study followed the guidelines adopted and approved by the institutional Ethical Review Committee at Dasman Diabetes Institute, Kuwait, in accordance with the declaration of Helsinki.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The genotype dataset generated and analyzed during the current study is available as Supplementary Dataset 1.

Acknowledgments

Kuwait Foundation for the Advancement of Sciences (KFAS) is acknowledged for institutional funding.

Conflicts of Interest

All the authors have read and approved the manuscript. The authors declare no conflict of interest.

References

  1. Hebbar, P.; Abubaker, J.A.; Abu-Farha, M.; Alsmadi, O.; Elkum, N.; Alkayal, F.; John, S.E.; Channanath, A.; Iqbal, R.; Pitkaniemi, J.; et al. Genome-wide landscape establishes novel association signals for metabolic traits in the arab population. Hum. Genet. 2021, 140, 505–528. [Google Scholar] [CrossRef] [PubMed]
  2. Ng, M.; Fleming, T.; Robinson, M.; Thomson, B.; Graetz, N.; Margono, C.; Mullany, E.C.; Biryukov, S.; Abbafati, C.; Abera, S.F.; et al. Global, regional, and national prevalence of overweight and obesity in children and adults during 1980–2013: A systematic analysis for the global burden of disease study 2013. Lancet 2014, 384, 766–781. [Google Scholar] [CrossRef] [Green Version]
  3. Abuyassin, B.; Laher, I. Obesity-linked diabetes in the arab world: A review. East. Mediterr. Health J. 2015, 21, 420–439. [Google Scholar] [CrossRef]
  4. Bamba, V.; Rader, D.J. Obesity and atherogenic dyslipidemia. Gastroenterology 2007, 132, 2181–2190. [Google Scholar] [CrossRef] [PubMed]
  5. Hebbar, P.; Nizam, R.; Melhem, M.; Alkayal, F.; Elkum, N.; John, S.E.; Tuomilehto, J.; Alsmadi, O.; Thanaraj, T.A. Genome-wide association study identifies novel recessive genetic variants for high tgs in an arab population. J. Lipid Res. 2018, 59, 1951–1966. [Google Scholar] [CrossRef] [Green Version]
  6. Graham, S.E.; Clarke, S.L.; Wu, K.H.; Kanoni, S.; Zajac, G.J.M.; Ramdas, S.; Surakka, I.; Ntalla, I.; Vedantam, S.; Winkler, T.W.; et al. The power of genetic diversity in genome-wide association studies of lipids. Nature 2021, 600, 675–679. [Google Scholar] [CrossRef]
  7. Kathiresan, S.; Melander, O.; Guiducci, C.; Surti, A.; Burtt, N.P.; Rieder, M.J.; Cooper, G.M.; Roos, C.; Voight, B.F.; Havulinna, A.S.; et al. Six new loci associated with blood low-density lipoprotein cholesterol, high-density lipoprotein cholesterol or triglycerides in humans. Nat. Genet. 2008, 40, 189–197. [Google Scholar] [CrossRef]
  8. Willer, C.J.; Sanna, S.; Jackson, A.U.; Scuteri, A.; Bonnycastle, L.L.; Clarke, R.; Heath, S.C.; Timpson, N.J.; Najjar, S.S.; Stringham, H.M.; et al. Newly identified loci that influence lipid concentrations and risk of coronary artery disease. Nat. Genet. 2008, 40, 161–169. [Google Scholar] [CrossRef] [Green Version]
  9. Schjoldager, K.T.; Clausen, H. Site-specific protein o-glycosylation modulates proprotein processing—Deciphering specific functions of the large polypeptide galnac-transferase gene family. Biochim. Biophys. Acta 2012, 1820, 2079–2094. [Google Scholar] [CrossRef]
  10. Schjoldager, K.T.; Vester-Christensen, M.B.; Bennett, E.P.; Levery, S.B.; Schwientek, T.; Yin, W.; Blixt, O.; Clausen, H. O-glycosylation modulates proprotein convertase activation of angiopoietin-like protein 3: Possible role of polypeptide galnac-transferase-2 in regulation of concentrations of plasma lipids. J. Biol. Chem. 2010, 285, 36293–36303. [Google Scholar] [CrossRef] [Green Version]
  11. Marucci, A.; Mangiacotti, D.; Trischitta, V.; di Paola, R. Galnt2 mrna levels are associated with serum triglycerides in humans. Endocrine 2016, 53, 331–334. [Google Scholar] [CrossRef] [PubMed]
  12. Di Paola, R.; Marucci, A.; Trischitta, V. Galnt2 effect on hdl-cholesterol and triglycerides levels in humans: Evidence of pleiotropy? Nutr. Metab. Cardiovasc. Dis. 2017, 27, 281–282. [Google Scholar] [CrossRef] [PubMed]
  13. Chen, J.; Guan, L.; Liu, H.; Liu, Q.; Fan, P.; Bai, H. Galnt2 gene variant rs4846914 is associated with insulin and insulin resistance depending on bmi in pcos patients: A case-control study. Reprod. Sci. 2021, 28, 1122–1132. [Google Scholar] [CrossRef] [PubMed]
  14. Marucci, A.; Antonucci, A.; de Bonis, C.; Mangiacotti, D.; Scarale, M.G.; Trischitta, V.; di Paola, R. Galnt2 as a novel modulator of adipogenesis and adipocyte insulin signaling. Int. J. Obes. 2019, 43, 2448–2457. [Google Scholar] [CrossRef] [Green Version]
  15. Buniello, A.; MacArthur, J.A.L.; Cerezo, M.; Harris, L.W.; Hayhurst, J.; Malangone, C.; Mc Mahon, A.; Morales, J.; Mountjoy, E.; Sollis, E.; et al. The nhgri-ebi gwas catalog of published genome-wide association studies, targeted arrays and summary statistics 2019. Nucleic Acids Res. 2019, 47, D1005–D1012. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Kulminski, A.M.; Loika, Y.; Nazarian, A.; Culminskaya, I. Quantitative and qualitative role of antagonistic heterogeneity in genetics of blood lipids. J. Gerontol. A Biol. Sci. Med. Sci. 2020, 75, 1811–1819. [Google Scholar] [CrossRef]
  17. Bentley, A.R.; Sung, Y.J.; Brown, M.R.; Winkler, T.W.; Kraja, A.T.; Ntalla, I.; Schwander, K.; Chasman, D.I.; Lim, E.; Deng, X.; et al. Multi-ancestry genome-wide gene-smoking interaction study of 387,272 individuals identifies new loci associated with serum lipids. Nat. Genet. 2019, 51, 636–648. [Google Scholar] [CrossRef]
  18. Spracklen, C.N.; Chen, P.; Kim, Y.J.; Wang, X.; Cai, H.; Li, S.; Long, J.; Wu, Y.; Wang, Y.X.; Takeuchi, F.; et al. Association analyses of east asian individuals and trans-ancestry analyses with european individuals reveal new loci associated with cholesterol and triglyceride levels. Hum. Mol. Genet. 2017, 26, 1770–1784. [Google Scholar] [CrossRef]
  19. Hoffmann, T.J.; Theusch, E.; Haldar, T.; Ranatunga, D.K.; Jorgenson, E.; Medina, M.W.; Kvale, M.N.; Kwok, P.Y.; Schaefer, C.; Krauss, R.M.; et al. A large electronic-health-record-based genome-wide study of serum lipids. Nat. Genet. 2018, 50, 401–413. [Google Scholar] [CrossRef]
  20. Abu-Farha, M.; Al-Khairi, I.; Cherian, P.; Chandy, B.; Sriraman, D.; Alhubail, A.; Al-Refaei, F.; AlTerki, A.; Abubaker, J. Increased angptl3, 4 and angptl8/betatrophin expression levels in obesity and t2d. Lipids Health Dis. 2016, 15, 181. [Google Scholar] [CrossRef] [Green Version]
  21. Abu-Farha, M.; Cherian, P.; Qaddoumi, M.G.; AlKhairi, I.; Sriraman, D.; Alanbaei, M.; Abubaker, J. Increased plasma and adipose tissue levels of angptl8/betatrophin and angptl4 in people with hypertension. Lipids Health Dis. 2018, 17, 35. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Abu-Farha, M.; Abubaker, J.; Noronha, F.; Al-Khairi, I.; Cherian, P.; Alarouj, M.; Bennakhi, A.; Elkum, N. Lack of associations between betatrophin/angptl8 level and c-peptide in type 2 diabetic subjects. Cardiovasc. Diabetol. 2015, 14, 112. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Bera, T.K. Bioelectrical impedance methods for noninvasive health monitoring: A review. J. Med. Eng. 2014, 2014, 381251. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Chang, C.C.; Chow, C.C.; Tellier, L.C.; Vattikuti, S.; Purcell, S.M.; Lee, J.J. Second-generation plink: Rising to the challenge of larger and richer datasets. Gigascience 2015, 4, 7. [Google Scholar] [CrossRef]
  25. Genomes Project, C.; Auton, A.; Brooks, L.D.; Durbin, R.M.; Garrison, E.P.; Kang, H.M.; Korbel, J.O.; Marchini, J.L.; McCarthy, S.; McVean, G.A.; et al. A global reference for human genetic variation. Nature 2015, 526, 68–74. [Google Scholar]
  26. Cunningham, F.; Allen, J.E.; Allen, J.; Alvarez-Jarreta, J.; Amode, M.R.; Armean, I.M.; Austine-Orimoloye, O.; Azov, A.G.; Barnes, I.; Bennett, R.; et al. Ensembl 2022. Nucleic Acids Res. 2022, 50, D988–D995. [Google Scholar] [CrossRef]
  27. Li, Q.; Yin, R.X.; Yan, T.T.; Miao, L.; Cao, X.L.; Hu, X.J.; Aung, L.H.; Wu, D.F.; Wu, J.Z.; Lin, W.X. Association of the galnt2 gene polymorphisms and several environmental factors with serum lipid levels in the mulao and han populations. Lipids Health Dis. 2011, 10, 160. [Google Scholar] [CrossRef] [Green Version]
  28. Richardson, T.G.; Sanderson, E.; Palmer, T.M.; Ala-Korpela, M.; Ference, B.A.; Davey Smith, G.; Holmes, M.V. Evaluating the relationship between circulating lipoprotein lipids and apolipoproteins with risk of coronary heart disease: A multivariable mendelian randomisation analysis. PLoS Med. 2020, 17, e1003062. [Google Scholar] [CrossRef] [Green Version]
  29. Hebbar, P.; Abu-Farha, M.; Mohammad, A.; Alkayal, F.; Melhem, M.; Abubaker, J.; Al-Mulla, F.; Thanaraj, T.A. Fto variant rs1421085 associates with increased body weight, soft lean mass, and total body water through interaction with ghrelin and apolipoproteins in arab population. Front. Genet. 2019, 10, 1411. [Google Scholar] [CrossRef]
  30. Li, X.; Zhang, Y.; Zhang, M.; Wang, Y. Galnt2 regulates angptl3 cleavage in cells and in vivo of mice. Sci. Rep. 2020, 10, 16168. [Google Scholar] [CrossRef]
  31. Fuior, E.V.; Gafencu, A.V. Apolipoprotein c1: Its pleiotropic effects in lipid metabolism and beyond. Int. J. Mol. Sci. 2019, 20, 5939. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Bis, J.C.; Kavousi, M.; Franceschini, N.; Isaacs, A.; Abecasis, G.R.; Schminke, U.; Post, W.S.; Smith, A.V.; Cupples, L.A.; Markus, H.S.; et al. Meta-analysis of genome-wide association studies from the charge consortium identifies common variants associated with carotid intima media thickness and plaque. Nat. Genet. 2011, 43, 940–947. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Franceschini, N.; Giambartolomei, C.; de Vries, P.S.; Finan, C.; Bis, J.C.; Huntley, R.P.; Lovering, R.C.; Tajuddin, S.M.; Winkler, T.W.; Graff, M.; et al. Gwas and colocalization analyses implicate carotid intima-media thickness and carotid plaque loci in cardiovascular outcomes. Nat. Commun. 2018, 9, 5141. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Guan, L.; Fan, P.; Liu, X.; Zhou, M.; Wu, Y.; Liu, R.; Liu, Y.; Bai, H. Maternal galnt2 variations affect blood pressure, atherogenic index, and fetal growth, depending on bmi in gestational diabetes mellitus. Front. Endocrinol. 2021, 12, 690229. [Google Scholar] [CrossRef]
Figure 1. Data distribution for the levels of BMI (A) and WC (B) in individuals with genotypes homozygous for major allele (AA) versus genotypes homozygous for effect allele (GG) or for the heterozygous genotypes (AG). The boxplot illustrates minimum and maximum values in the range, first quartile, median, and third quartile for each of the genotype groups.
Figure 1. Data distribution for the levels of BMI (A) and WC (B) in individuals with genotypes homozygous for major allele (AA) versus genotypes homozygous for effect allele (GG) or for the heterozygous genotypes (AG). The boxplot illustrates minimum and maximum values in the range, first quartile, median, and third quartile for each of the genotype groups.
Genes 13 01201 g001
Figure 2. Data distribution for the levels of lipid traits, namely HDL (A), LDL (B), total cholesterol (C), and TG (D), in individuals with genotypes homozygous for major allele (AA) versus genotypes homozygous for effect allele (GG) or for the heterozygous genotypes (AG). The boxplot illustrates minimum and maximum values in the range, first quartile, median, and third quartile for each of the genotype groups.
Figure 2. Data distribution for the levels of lipid traits, namely HDL (A), LDL (B), total cholesterol (C), and TG (D), in individuals with genotypes homozygous for major allele (AA) versus genotypes homozygous for effect allele (GG) or for the heterozygous genotypes (AG). The boxplot illustrates minimum and maximum values in the range, first quartile, median, and third quartile for each of the genotype groups.
Genes 13 01201 g002
Figure 3. Data distribution for the levels of ANGPTL’s, namely ANGPTL3 (A), ANGPTL4 (B), and ANGPTL8 (C), in individuals with genotypes homozygous for major allele (AA) versus genotypes homozygous for effect allele (GG) or for the heterozygous genotypes (AG). The boxplot illustrates minimum and maximum values in the range, first quartile, median, and third quartile for each of the genotype groups.
Figure 3. Data distribution for the levels of ANGPTL’s, namely ANGPTL3 (A), ANGPTL4 (B), and ANGPTL8 (C), in individuals with genotypes homozygous for major allele (AA) versus genotypes homozygous for effect allele (GG) or for the heterozygous genotypes (AG). The boxplot illustrates minimum and maximum values in the range, first quartile, median, and third quartile for each of the genotype groups.
Genes 13 01201 g003
Figure 4. Data distribution for the levels of Apolipoproteins, namely ApoA1A (A), ApoA2 (B), ApoB (C), and ApoC1 (D), in individuals with genotypes homozygous for major allele (AA) versus genotypes homozygous for effect allele (GG) or for the heterozygous genotypes (AG). The boxplot illustrates minimum and maximum values in the range, first quartile, median, and third quartile for each of the genotype groups.
Figure 4. Data distribution for the levels of Apolipoproteins, namely ApoA1A (A), ApoA2 (B), ApoB (C), and ApoC1 (D), in individuals with genotypes homozygous for major allele (AA) versus genotypes homozygous for effect allele (GG) or for the heterozygous genotypes (AG). The boxplot illustrates minimum and maximum values in the range, first quartile, median, and third quartile for each of the genotype groups.
Genes 13 01201 g004
Figure 5. Interactions observed in HDL with PBF (negative correlation) (A), ApoA1A (B), ApoC1 (C), and ApoB48 (D) in interaction with the AG+GG carrier genotype at the study variant.
Figure 5. Interactions observed in HDL with PBF (negative correlation) (A), ApoA1A (B), ApoC1 (C), and ApoB48 (D) in interaction with the AG+GG carrier genotype at the study variant.
Genes 13 01201 g005
Table 1. Clinical characteristics of the study cohort.
Table 1. Clinical characteristics of the study cohort.
TraitsAll Participants (Mean ± SD)Obese (n = 143) (Mean ± SD)Non-Obese (n = 135) (Mean ± SD)p-Value for
Differences in the Mean Values
between the Two Sub-Cohorts @
Male:Female125:15362:8163:720.660
Age (years)46.25 ± 12.3848.31 ± 12.7044.30 ± 12.376.7 × 10−3
Weight (kg)81.40 ± 16.2393.18 ± 11.7370.45 ± 11.48<1.0 × 10−4
BMI (kg/m2)29.93 ± 5.1734.25 ± 3.0025.85 ± 3.04<1.0 × 10−4
WC (cm)99.36 ± 13.36108.02 ± 9.6689.98 ± 10.09<1.0 × 10−4
PBF (%)35.52 ± 5.6838.19 ± 4.7832.67 ± 5.17<1.0 × 10−4
LDL (mmol/L)3.38 ± 0.943.15 ± 0.973.10 ± 0.950.633
HDL (mmol/L)1.20 ± 0.321.17 ± 0.321.22 ± 0.320.260
TC (mmol/L)5.27 ± 1.045.04 ± 1.065.00 ± 1.120.742
TG (mmol/L)1.22 ± 0.591.34 ± 0.551.11 ± 0.611.0 × 10−3
FPG (mmol/L)5.77 ± 1.246.03 ± 1.305.57 ± 1.145.0 × 10−3
HbA1c (%)6.31 ± 1.296.73 ± 1.485.90 ± 0.92<1.0 × 10−4
ANGPTL3 (ng/mL)37.42 ± 10.2938.84 ± 10.8135.83 ± 9.483.9 × 10−2
ANGPTL4 (ng/mL)153.4 ± 51.9145.0 ± 42.4161.1 ± 58.62.7 × 10−2
ANGPTL8 (pg/mL)1321.74 ± 861.41465.4 ± 878.91176.5 ± 822.72.1 × 10−2
ApoA1 (mg/L)1600.9 ± 429.71612.0 ± 439.11590.2.6 ± 422.30.726
ApoA2 (mg/L)1001.4 ± 426.71023.8 ± 433.8979.2 ± 420.50.450
ApoB (mg/L)1829.2 ± 774.91874.6 ± 824.11783.2 ± 723.00.391
ApoC1 (ng/mL)477.7 ± 109.9477.2 ± 116.2478.3 ± 103.50.945
Diabetes status (Yes:No)120:15873:6247:967.0 × 10−4
Anti-diabetic medication (Yes:No)101:17767:6734:109<1.0 × 10−4
Lipid lowering medication (Yes:No)88:19058:7630:1131.1 × 10−4
@, Significant p-values are shown in bold font. Comparisons were performed using Student’s t-test for quantitative variables and Chi-squared test for categorical variables to determine significance. p-values ≤ 0.05 were considered significant.
Table 2. Results of allelic association of rs4846914 with obesity status and diabetes status.
Table 2. Results of allelic association of rs4846914 with obesity status and diabetes status.
CategoryAllele Frequency (G/A)OR (CI 95%) *p-Value *
All0.43/0.57--
Obese0.47/0.531.47 [1.04–2.05]3.1 × 10−2
Non-obese0.37/0.63
Diabetic0.44/0.561.12 [0.79–1.58]0.540
Non-diabetic0.42/0.58
*, Allele frequency differences between obese vs. non-obese and diabetic vs. non-diabetic individuals were tested using Fisher exact test.
Table 3. Results of association tests for the study variant with G as effect allele with the quantitative traits and biomarkers, using genetic model based on additive mode of inheritance.
Table 3. Results of association tests for the study variant with G as effect allele with the quantitative traits and biomarkers, using genetic model based on additive mode of inheritance.
TraitCorrection *Sample Sizeβp-Value @
BMIR2741.3422.2 × 10−3
R + DS2741.3122.4 × 10−3
R + OS2731.312.6 × 10−3
WeightR2703.832.7 × 10−3
R + DS2703.7332.8 × 10−3
R + OS2693.8372.3 × 10−3
HDLR256−0.0562.8 × 10−2
R + DS256−0.0523.9 × 10−2
R + OS255−0.0524.1 × 10−2
LDLR2650.0650.441
R + DS2650.0700.410
R + OS2640.0730.392
TCR2690.0230.803
R + DS2690.0280.760
R + OS2680.0320.734
TGR2560.0860.881
R + DS2560.0790.103
R + OS2550.0810.102
ANGPTL3R1932.0295.2 × 10−2
R + DS1932.0345.2 × 10−2
R + OS1932.035.3 × 10−2
ANGPTL4R195430.50.932
R + DS195−174.00.970
R + OS195953.40.850
ANGPTL8R186−28.270.723
R + DS186−28.550.654
R + OS1867.1420.911
ApoC1R17832.006.0 × 10−3
R + DS17832.036.1 × 10−3
R + OS17730.888.9 × 10−3
ApoA1AR19388.103.9 × 10−2
R + DS19388.124.0 × 10−2
R + OS19279.516.4 × 10−2
ApoA2R20484.734.5 × 10−2
R + DS20484.804.2 × 10−2
R + OS20486.304.4 × 10−2
ApoBR206147.305.1 × 10−2
R + DS206149.604.7 × 10−2
R + OS205141.106.5 × 10−2
PBFR1750.6320.175
R + DS1750.5210.263
R + OS1750.3240.321
*, R’ stands for ‘regular correction: the model is adjusted for age, sex, and medication for diabetes and lowering lipid levels’; ‘DS’ stands for further correction for diabetes status in addition to the regular correction; ‘OS’ stands for further correction for obesity status in addition to the regular correction. @, p-values ≤ 0.05 are in bold and italic font, with those ≤0.006 further underlined.
Table 4. Linear regression model illustrating the link between the study variant and interaction of HDL with PBF, ApoA1A, ApoC1, and ApoB48. Results are shown for combined minor allele homozygous and heterozygous genotypes (AG and GG) regressed against the major AA genotype @.
Table 4. Linear regression model illustrating the link between the study variant and interaction of HDL with PBF, ApoA1A, ApoC1, and ApoB48. Results are shown for combined minor allele homozygous and heterozygous genotypes (AG and GG) regressed against the major AA genotype @.
Trait
(Response/Dependent Variable)
Genotypes and/or Interacting Traits (Predict
Variable)
EstimateStd. Errorp-ValueAdj. R-SquareModel p-Value
Model: HDL~rs4846914 + Age + Sex + PBF + rs4846914 * PBF
HDLIntercept0.5340.2710.0510.2334.32 × 10−9
AG + GG0.6770.3180.034
PBF0.00590.00830.474
AG + GG*PBF−0.00220.00880.013
Model: HDL~rs4846914 + Age + Sex + ApoA1A + rs4846914 * ApoA1A
HDLIntercept0.3000.1750.0870.2392.84 × 10−10
AG + GG0.3780.1670.024
ApoA1A3.05 × 10−78.69 × 10−85.6 × 10−4
AG + GG*ApoA1A−2.96 × 10−71.04 × 10−75.1 × 10−3
Model: HDL~rs4846914 + Age + Sex + ApoC1 + rs4846914 * ApoC1
HDLIntercept0.4400.1970.0260.2211.64 × 10−8
AG + GG0.3690.2020.069
ApoC18.14 × 10−73.46 × 10−70.020
AG + GG*ApoC1−9.47 × 10−74.22 × 10−70.026
Model: HDL~rs4846914 + Age + Sex + ApoB48 + rs4846914 * ApoB48
HDLIntercept0.4820.1550.0020.1992.9 × 10−8
AG + GG0.2310.1220.061
ApoB485.08 × 10−61.62 × 10−60.002
AG + GG*ApoB48−4.08 × 10−61.85 × 10−60.028
@, Multivariate linear regression with correction for age and sex was performed to determine correlations among the risk variant, biomarker levels, and the associated metabolic traits. Relationships between traits and biomarkers were denoted by percentage of response variable variation (R2), standardized β-coefficients (β1), and significance of test (p) for the reference versus alternate genotype distributions.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Qaddoumi, M.; Hebbar, P.; Abu-Farha, M.; Al Somaly, A.; Melhem, M.; Al-Kayal, F.; AlKhairi, I.; Cherian, P.; Alanbaei, M.; Al-Mulla, F.; et al. GALNT2 rs4846914 SNP Is Associated with Obesity, Atherogenic Lipid Traits, and ANGPTL3 Plasma Level. Genes 2022, 13, 1201. https://doi.org/10.3390/genes13071201

AMA Style

Qaddoumi M, Hebbar P, Abu-Farha M, Al Somaly A, Melhem M, Al-Kayal F, AlKhairi I, Cherian P, Alanbaei M, Al-Mulla F, et al. GALNT2 rs4846914 SNP Is Associated with Obesity, Atherogenic Lipid Traits, and ANGPTL3 Plasma Level. Genes. 2022; 13(7):1201. https://doi.org/10.3390/genes13071201

Chicago/Turabian Style

Qaddoumi, Mohammad, Prashantha Hebbar, Mohamed Abu-Farha, Aseelah Al Somaly, Motasem Melhem, Fadi Al-Kayal, Irina AlKhairi, Preethi Cherian, Muath Alanbaei, Fahd Al-Mulla, and et al. 2022. "GALNT2 rs4846914 SNP Is Associated with Obesity, Atherogenic Lipid Traits, and ANGPTL3 Plasma Level" Genes 13, no. 7: 1201. https://doi.org/10.3390/genes13071201

APA Style

Qaddoumi, M., Hebbar, P., Abu-Farha, M., Al Somaly, A., Melhem, M., Al-Kayal, F., AlKhairi, I., Cherian, P., Alanbaei, M., Al-Mulla, F., Abubaker, J., & Thanaraj, T. A. (2022). GALNT2 rs4846914 SNP Is Associated with Obesity, Atherogenic Lipid Traits, and ANGPTL3 Plasma Level. Genes, 13(7), 1201. https://doi.org/10.3390/genes13071201

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop