Take a look at the Recent articles

The rare allele of DGKZ SNP rs10838599 is associated with variability in HDL-cholesterol levels among severely obese patients

Stéphanie Bégin

Institute of Nutrition and Functional Foods (INAF), Laval University, Canada

School of Nutrition, Laval University, Canada

Frédéric Guénard

Institute of Nutrition and Functional Foods (INAF), Laval University, Canada

School of Nutrition, Laval University, Canada

André Tchernof

School of Nutrition, Laval University, Canada

Quebec Heart and Lung Institute, Canada

Yves Deshaies

Quebec Heart and Lung Institute, Canada

Department of Medicine, Laval University, Canada

Louis Pérusse

Institute of Nutrition and Functional Foods (INAF), Laval University, Canada


Department of Kinesiology, Laval University, Canada

Frédéric-Simon Hould

Quebec Heart and Lung Institute, Canada

Department of Surgery, Laval University, Canada

Stéfane Lebel

Quebec Heart and Lung Institute, Canada

Department of Surgery, Laval University, Canada

Picard Marceau

Quebec Heart and Lung Institute, Canada

Department of Surgery, Laval University, Canada

Marie-Claude Vohl

Institute of Nutrition and Functional Foods (INAF), Laval University, Canada

School of Nutrition, Laval University, Canada

E-mail : marie-claude.vohl@fsaa.ulaval.ca

DOI: 10.15761/IOD.1000148

Article
Article Info
Author Info
Figures & Data

Abstract

Introduction: Diacylglycerol kinase-zeta, one of the ten isoforms of DGKs expressed in mammals is an important enzyme of lipid metabolism. It catalyzes the interconversion of diacylglycerol and phosphatidic acid, two major second messengers. Its gene DGKZ has been previously identified as being overexpressed and undermethylated in visceral adipose tissue of patients with (MetS+) versus without (MetS-) the metabolic syndrome (MetS).

Objective: The aim of this study was to investigate the associations between DGKZ gene polymorphisms (SNPs) and phenotypes related to MetS (BMI, waist girth, CRP, fasting glucose, lipid profile (triglycerides, total-cholesterol, LDL-cholesterol and HDL-cholesterol (HDL-C)), resting systolic and diastolic blood pressures).

Methods: The study sample included 1752 severely obese participants who underwent bariatric surgery. Associations between the five selected tSNPs of DGKZ and features of the MetS were tested. The effects of these SNPs on DGKZ methylation and expression levels were tested in subgroups of 32 and 14 obese subjects, respectively. Correlations between methylation and expression levels were also computed.

 Results: Homozygotes for the rare allele of rs10838599 displayed higher plasma HDL-C concentrations compared to the other genotype groups (p=0.03). For gene methylation, only a trend with the cg05412031 CpG site (p=0.09) was found for the single significantly phenotype-associated SNP. There was no significant correlation between DGKZ methylation at cg05412031 and expression levels.

Conclusion: These results suggest that DGKZ SNP rs10838599 modulates plasma HDL-C levels thereby its gene contributes to the inter-individual variability observed in the cardiometabolic risk profile of patients with severe obesity.

Key words

 diacylglycerol, epigenetics, genomics, lipids, metabolic syndrome, obesity

Abbreviations

BMI: Body Mass Index; CRP: C-Reactive Protein; DAG: Diacylglycerol; DBP: Diastolic Blood Pressure; DGK: Diacylglycerol Kinase; HDL-C: High-Density Lipoprotein Cholesterol; HMZ: Homozygote; HTZ: Heterozygote; HWE: Hardy-Weinberg Equilibrium; LDL-C: Low-Density Lipoprotein Cholesterol; MetS: Metabolic Syndrome; PA: Phosphatidic Acid; SBP: Systolic Blood Pressure; SNP: Single Nucleotide Polymorphism; TG: Triglycerides; VAT: Visceral Adipose Tissue

Introduction

The metabolic syndrome (MetS) is a constellation of risk factors associated with cardiovascular diseases, and could explain the mortality related to obesity [1]. However, not all obese individuals develop MetS [2]. Consequently, the analysis of obesity etiology is critical to better understand and prevent its condition of epidemic proportions worldwide [3]. Differences in DNA sequence of obese individuals affected (MetS+) and unaffected (MetS-) by MetS can be examined. However, genetic factors are not the only contributor to the heterogeneous obesity phenotype. Epigenetic variations have also been associated with the development of metabolic disorders [4,5], including low HDL-cholesterol (HDL-C) levels, as observed among obese individuals [6]. Beside cholesterol transport, HDL-C particles also have antioxidant, anti-inflammatory and antithrombotic properties [7]. DGKZ gene is coding for the protein diacylglycerol kinase (DGK)-zeta, one of the ten DGK isoforms expressed in mammals [8]. It catalyzes the conversion of diacylglycerol (DAG) into phosphatidic acid (PA) and inversely [9], two essential lipid metabolic intermediates in the constitution of membranes and signal transduction [10].

The DGKZ gene was identified by our team among a list of 488 additional loci as being differentially expressed in visceral adipose tissue (VAT) of MetS+ versus MetS- subjects [11]. Afterwards, DGKZ was shown to be differentially methylated in VAT of obese patients discordant for MetS [12]. We then hypothesized that DNA genetic/epigenetic variations in DGKZ are associated with MetS components and that CpG methylation levels are associated with the presence of obesity-related metabolic complications based on the role of methylation in the regulation of gene expression [13]. The aim of this study was therefore to determine the effects of DGKZ genetic and epigenetic variations on the metabolic complications of obesity.

Material and Methods

Subjects Selection

Among patients undergoing bariatric surgery (biliopancreatic diversion with duodenal switch) at the Québec Heart and Lung Institute (IUCPQ) (Québec City, Québec, Canada), 1752 severely obese men (N=545) and women (N=1207) were consecutively recruited. There were three inclusion criteria: the participants should present a BMI higher or equivalent to 35, be French-Canadian and attend for a bariatric surgery at the IUCPQ. The surgical protocol is described elsewhere [14]. Body weight, height, waist girth, plasma lipid (triglycerides (TG), total-cholesterol (total-C), LDL-cholesterol (LDL-C) and HDL-C) and resting systolic (SBP) and diastolic (DBP) blood pressures were measured before the surgery using standardized procedures [15]. Body mass index (BMI) was calculated as weight in kilograms divided by height in meters squared and C-reactive protein (CRP) plasma concentrations were measured with a high sensitivity-CRP immunoassay using a monoclonal antibody coated on polystyrene particles. The presence of MetS was determined using the National Cholesterol Education Program Adult Treatment Panel III (NCEP-ATPIII) criteria [16]. Tissue specimens were obtained from the Biobank of the IUCPQ according to institutionally-approved management modalities. All subjects included in the study provided a written informed consent previously approved by the local ethics committee.

Genotyping

DNA was extracted from blood buffy coat using the GenElute Blood Genomic DNA kit (Sigma, St Louis, MO). Four tagging single nucleotide polymorphisms (tSNPs) were selected for the analysis based on differential expression and methylation of the DGKZ gene established in a previous dataset [11,12]. These SNPs were genotyped using the QuantStudio™ 12K Flex OpenArray® AccuFill™ system (Applied Biosystems), and analyzed with TaqMan Genotyper v1.3 (Life Technologies). The inclusion of the rs34470543 SNP in the current study is attributable to its close localization to the cg12856521 CpG site, previously found to be differentially methylated between MetS+ and MetS- subjects [17]. The tagger selection algorithm of the Haploview software (pairwise tagging, R2 ≥ 0.80) was used for analysis of selected SNPs spanning promoter (2 kb), coding and intronic regions of the DGKZ gene in addition to the 3’ gene region (2 kb) [17]. This strategy allowed covering 100% of the genetic variability of the common polymorphisms (MAF ≥ 1%) at the DGKZ locus in the Caucasian population (CEU HapMap).

DNA methylation analysis

From the original sample, a subgroup of 32 obese individuals was used for DNA methylation analysis, including the 14 subjects (MetS+, N=7; MetS-, N=7) selected for gene expression profiling [11]. The 18 additional obese individuals (MetS+, N=9; MetS-, N=9) were selected to fulfill initial selection criteria [11] and to represent extremes of the MetS diagnosis criteria spectrum. Samples of 200 mg of VAT were used to achieve genomic DNA extraction using the DNeasy Blood and Tissue kit (QIAGEN, Mississauga, Ontario, Canada). Bisulfite conversion was conducted on 1 µg of DNA, and quantitative DNA methylation analysis was carried out at the McGill University and Génome Québec Innovation Centre (Montréal, Canada). Infinium HumanMethylation450 BeadChip (Illumina, San Diego, CA) covering more than 485000 methylation sites at single-nucleotide resolution were processed according to the manufacturer’s instructions. The GenomeStudio software version 2011.1 (Illumina Inc.) and the methylation module were used for visualisation and analysis of the methylation data. Methylation levels (beta values; β) were estimated as the ratio of signal intensity of the methylated alleles to the sum of methylated and unmethylated intensity signals of the alleles (β value = C/(T+C)). We applied internal control probe pairs for data correction (background subtraction and normalization). The GenomeStudio Methylation Module was used to extract CpG sites located within the DGKZ locus and promoter region, thus leading to a total of 66 CpG sites analyzed. Based on the differentially methylated CpG site cg12856521 (chr11:46389249, genome build 37), DGKZ was identified as differentially methylated (Diffscore = -235.58).

Gene expression

The protocol leading to the selection of DGKZ for the current study and its expression data presented here were retrieved from a previous study aimed at the identification of differentially methylated genes between MetS+ and MetS- severely obese men [11]. Even if the DGKZ gene is not listed in the tables due to the large number of genes discovered, this investigation led to its determination as being overexpressed within VAT of MetS+ versus MetS- subjects (1.21-fold; p=0.004).

Statistical analysis

At first, differences in metabolic biomarkers were tested between MetS+ and MetS- individuals of the study cohort. Gene expression level differences derived from microarray analysis between MetS+ and MetS- groups were tested using the unpaired Student’s t test. For the five selected DGKZ tSNPs, Hardy-Weinberg equilibrium (HWE) was verified. For SNPs with genotype frequency of rare homozygotes below 5% (2 SNPs: rs3809047 and rs3740976), homozygotes for the rare allele were merged to heterozygotes for statistical analysis. The associations between SNPs and MetS components were tested using the GLM procedure with adjustments for age, sex, BMI and medication to treat MetS features when appropriate. When a significant SNP effect was identified, all pairwise comparisons between genotype groups were tested using least squares means. Pearson correlation coefficients were computed to assess the relationship between gene methylation, expression and MetS components. P values were calculated for all associations and were considered statistically significant if the P value was less than 0.05. Phenotypic data are presented as mean ± SD. Statistical analyses were performed with SAS software version 9.3.

Results

Cohort description

The current study included 1752 subjects (545 men and 1207 women) who underwent bariatric surgery. From those, 1745 were classified as being MetS+ (n=1428) or MetS- (n=317 or 20.19%). Patients’ characteristics are presented in Table 1. Briefly, the MetS+ group was older and showed significantly higher waist girth, fasting glucose, TG, total-C/HDL-C ratio and blood pressure (SBP and DBP), as well as lower plasma total-C, LDL-C and HDL-C levels than the MetS- group.

Table 1: Subjects’ characteristics

All

MetS+

MetS-

P-value

Number of Subjects (% male)

1752 (31.1)

1428 (33.6)

317 (20.2)

-

Age (years)

43.02 ± 10.62

43.96 ± 10.58

38.66 ± 9.75

<0.0001

BMI (kg/m2)

51.8 ± 8.9

52.0 ± 9.0

51.2 ± 8.7

0.32

Waist girth (cm)

140.6 ± 17.7

141.7 ± 17.4

135.5 ± 18.5

<0.0001

CRP (mg/L)

11.04 ± 9.12

11.22 ± 9.16

10.42 ± 9.00

0.05

Fasting glucose (mmol/L)

6.53 ± 2.32

6.85 ± 2.41

5.11 ± 0.98

<0.0001

Lipid profile (mmol/L)

    TG

1.83 ± 1.08

1.98 ± 1.12

1.17 ± 0.37

<0.0001

    Total-C

4.70 ± 0.95

4.68 ± 0.98

4.80 ± 0.82

0.05

    LDL-C

2.67 ± 0.83

2.65 ± 0.84

2.80 ± 0.75

0.02

    HDL-C

1.24 ± 0.35

1.18 ± 0.32

1.49 ± 0.37

<0.0001

    Total-C/HDL-C

4.03 ± 1.29

4.17 ± 1.34

3.33 ± 0.74

<0.0001

Blood Pressure (mmHg)

     SBP

139.0 ± 17.1

140.0 ± 17.3

134.6 ± 15.6

0.0006

     DBP

83.8 ± 11.5

84.0 ± 11.7

83.1 ± 10.3

0.04

Boldface values represent statistically significant differences (P ≤ 0.05). Values are presented as mean ± SD. Abbreviations: MetS: Metabolic Syndrome; BMI: Body Mass Index; CRP: C-Reactive Protein; TG: Triglycerides; Total-C: Total Cholesterol; LDL-C: Low-Density Lipoprotein Cholesterol; HDL-C: High-Density Lipoprotein Cholesterol; SBP: Systolic Blood Pressure; DBP: Diastolic Blood Pressure; SD: Standard Deviation.

Identification of DGKZ SNPs and selection of tSNPs

Among the 5 selected SNPs, 4 were located in introns whereas SNP rs34470543 was located in an exon (Table 2). These 5 SNPs were further genotyped in the entire cohort of 1752 participants. All SNPs were in HWE and genotype distribution is shown in Table 2.

Table 2: Genotype distribution and localization of selected DGKZ SNPs

SNPs

Number of genotypes

Common HMZ

HTZ

Rare HMZ

Localization

Other designation

Region

MAF

HWE                        P-values

rs11038866

1735

804

755

176

chr11:46355932

c.212+895C>G

Intron 1

0.32

0.95

rs3809047

1704

1384

309

11

chr11:46367702

c.212+12665T>C

Intron 1

0.10

0.16

rs10838599

1751

864

742

145

chr11:46378708

c.213-10133G>T

Intron 1

0.29

0.42

rs34470543

1743

1743

0

0

chr11:46389250

c.369_370insC

Exon 3

1

-

rs3740976

1738

1315

396

27

chr11:46393757

c.975+37C>T

Intron 10

0.13

0.65

aPosition according to genome build 37

bReference sequence : NM_201532.

Abbreviations: SNP: Single Nucleotide Polymorphism; HMZ: Homozygote; HTZ: Heterozygote; MAF: Minor Allele Frequency; HWE: Hardy-Weinberg Equilibrium.

Association of DGKZ SNPs

Associations between DGKZ SNPs and fasting plasma concentrations of glucose, TG, total-C, LDL-C, HDL-C, CRP and blood pressure (DBP and SBP) were tested taking into account the potential confounding effects of age, sex, and BMI. A significant association was found for rs10838599, homozygotes for the rare allele of rs10838599 displaying elevated HDL-C levels (p=0.03) compared to the two other genotype groups (Table 3). A trend toward higher DBP (p=0.08; 83.67 ± 10.91 vs. 83.43 ± 11.30 vs. 85.71 ± 11.97) was also observed. For rs3740976, a trend toward higher DBP (p=0.09; 84.32 ± 10.85 vs. 83.51 ± 11.30) was observed in rare allele carriers.

Table 3: Significant genotype differences identified between fasting glucose, lipid profile and blood pressure.

rs10838599

rs3740976

Means

Means

Phenotypes

Common HMZ (WT)

HTZ

Rare HMZ

P values

Common HMZ (WT)

Rare allele carrier

P values

Number of subjects (N)

858

740

144

1309

422

BMI (kg/m2)

51.9 ± 8.8

51.5 ± 8.3

51.51 ± 9.1

0.70

51.6 ± 8.6

51.9 ± 8.5

0.69

Waist girth (cm)

141.0 ± 17.9

140.1 ± 17.5

140.7 ± 18.3

0.53

140.1 ± 17.7

141.9 ± 17.7

0.14

CRP Protein (mg/L)

10.70 ± 8.21

10.79 ± 8.43

11.46 ± 8.75

0.59

10.75 ± 8.28

11.03 ± 8.61

0.54

Fasting glucose

6.51 ± 2.27

6.56 ± 2.32

6.56 ± 2.60

0.65

6.55 ± 2.30

6.53 ± 2.42

0.98

Lipid profile (mmol/L)

    TG

1.81 ± 0.97

1.84 ± 1.20

1.88 ± 0.99

0.32

1.82 ± 0.98

1.85 ± 1.30

0.48

    Total-C

4.69 ± 0.95

4.67 ± 0.92

4.86 ± 0.99

0.10

4.70 ± 0.93

4.68 ± 0.97

0.83

    LDL-C

2.67 ± 0.83

2.66 ± 0.79

2.80 ± 0.90

0.36

2.68 ± 0.82

2.66 ± 0.84

0.94

    HDL-C

1.24 ± 0.37

1.21 ± 0.30

1.31 ± 0.46

0.03

1.24 ± 0.37

1.23 ± 0.30

0.48

    Total-C/HDL-C

4.03 ± 1.41

4.03 ± 1.16

3.98 ± 1.20

0.48

4.04 ± 1.31

3.99 ± 1.23

0.45

Blood pressure (mmHg)

    SBP

138.9 ± 16.9

138.7 ± 16.8

139.3 ± 15.6

0.74

138.5 ± 16.5

139.6 ± 17.3

0.44

    DBP

83.7 ± 10.9

83.4 ± 11.3

85.7 ± 12.0

0.08

83.5 ± 11.3

84.3 ± 10.9

0.09

Boldface values represent statistically significant results (P ≤ 0.05).

aValues presented (means ± SD) are untransformed and unadjusted.

bP-values obtained are adjusted for the effect of age, sex, BMI and medication use except for BMI and waist girth which were adjusted for age and sex.

Abbreviations: HMZ: Homozygote; WT: Wild-Type; HTZ: Heterozygote; N: Number; BMI: Body Mass Index; CRP: C-Reactive Protein; TG: Triglycerides; Total-C: Total Cholesterol; LDL-C: Low-Density Lipoprotein Cholesterol; HDL-C: High-Density Lipoprotein Cholesterol; SBP: Systolic Blood Pressure; DBP: Diastolic Blood Pressure; SD: Standard Deviation

Gene methylation and expression analysis

To better understand the mechanism linking DGKZ genetic variations and cardiometabolic risk factors, we investigated the impact of these SNPs on gene methylation and expression levels in subgroups of 32 and 14 participants, respectively. Methylation levels were available for 66 CpG sites randomly distributed within the gene, their precise locations are presented in Table 4. The impact on methylation levels was tested only for SNP (rs10838599) found to be significantly associated with one cardiometabolic risk factor. For rs10838599, only a trend for genotype differences in methylation levels at the cg05412031 CpG site (p=0.09) was noted. No significant correlation was found between methylation levels at this CpG site and DGKZ expression levels.

Table 4: Association of phenotype-associated SNP rs10838599 with gene methylation levels

2021 Copyright OAT. All rights reserv

SNP

cg18908017   Promoter

cg24878090   Promoter

cg18765542   Promoter

cg21029403    Promoter

cg26650731   Promoter

cg05412031      Promoter

cg25151353   Promoter

cg07276415  Promoter

cg05214390    Promoter

cg14503180    Promoter

cg10503202    Promoter

rs10838599

0.95

0.99

0.58

0.23

0.91

0.09

0.24

0.94

0.65

0.39

0.10

SNP

cg15825186       Exon 1

cg24407308    Exon 1

cg04854089       Intron 1

cg10059056    Intron 1

cg09802018     Intron 1

cg22948808   Intron 1

cg22707438   Intron 1

cg15925492  Intron 1

cg20091215   Intron 1

cg06061966   Intron 1

cg07998213    Intron 1

rs10838599

0.27

0.36

0.38

0.44

0.94

0.87

0.48

0.93

0.34

0.70

0.86

SNP

cg05576959           Intron 1

cg13247398      Intron 1

cg20029347     Intron  1

cg21663341     Intron  1

cg03519157      Intron  1

cg20666386      Intron  1

cg07701241   Intron 1

cg00530720         Intron 1

cg15868421  Intron 1

cg09039436    Intron 1

cg03044514  Intron 1

rs10838599

0.55

0.66

0.40

0.51

0.68

0.63

0.34

0.67

0.93

0.99

0.73

SNP

cg20510601  Intron 1

cg22675791     Intron 1

cg18568067   Intron 1

cg12087471  Intron 1

cg07092212  Intron 1

cg01132064     Intron 1

cg17266233   Intron 1

cg07121644  Intron 1

cg21770622   Intron 1

cg02865822  Intron 1

cg18337963  Intron 1

rs10838599

0.87

0.44

0.92

0.20

0.28

0.45

0.47

0.74

0.52

0.68

0.54

SNP

cg02211741  Intron 1

cg26473110   Intron 1

cg14787093      Intron 1

cg20204701   Intron 1

cg26963044   Intron 1

cg24713122   Intron 2

cg12856521         Exon 3

cg00107629           Intron 3

cg24126592   Intron 3

cg08768904        Exon 4

cg16769381    Exon 5

rs10838599

0.27

0.59

0.53

0.67

0.84

0.48

0.68

0.91

0.32

0.73

0.12

SNP

cg27651070   Intron 6

cg01529068    Intron 9

cg16977596    Intron 16

cg12927726     Exon 20

cg05897087   Intron 20

cg02154997    Exon 21

cg12592409   Intron 27

cg11436767     Exon 30

cg07052627    Exon 31

cg06266097    Exon 31            (3’-UTR)

cg11203377    Exon 31     (3’-UTR)

rs10838599

0.20

0.68

0.60

0.15

0.39

0.48

0.96

0.61

0.56

0.70

0.98

P-values for associations were obtained from a subset of 14 obese subjects (7 MetS+ and 7 MetS-). CpG sites positions according to genome build 37.

Discussion

Our group previously observed that DGKZ was overexpressed [11] and undermethylated [12] in VAT of MetS+ vs MetS- severely obese subjects. DGKZ was thus revealed as a candidate gene having a potential effect on the development of obesity and MetS. Accordingly, it was appealing to test its associations with features of MetS. The ten diacylglycerol kinase isoforms expressed in mammals [8], including DGK-zeta, are important enzymes of lipid metabolism, because they are involved in the regulation of two major lipid messengers [18], controlling DAG and PA levels in cells [19]. The second messenger DAG regulates the activity of Rac1 and RhoA, two Rho GTPases [20] affecting the actin filaments of the cytoskeleton [21]. An increase in DAG, modulated by DGK, appears to reduce glycogen synthesis stimulated by insulin in liver, also decreasing the suppression of gluconeogenesis [21]. Therefore, the activity of DAG and PA may be altered by obesity and type 2 diabetes [22]. Specifically, the DGK-zeta isoform is expressed in diverse tissues, especially in brain [23]. At first, it is in the cell cytoplasm but can be translocated to various organelles [24]. Such translocation may be attributable to its unique myristoylated alanine-rich C-kinase substrate (MARCKS) phosphorylation domain [19,23]. The overexpression of DGKZ has been associated with decreased cell cycle progression [19], as DAG is important for this process as well as for cell differentiation [25]. Moreover, a previous study reported a relationship between DGK-delta and insulin resistance [26], but has not been reported for other isoforms. Insulin resistance may be an important event in the development of MetS, inducing hyperglycemia and dyslipidemia, although not all MetS+ patients are insulin resistant [27,28]. In addition, DAG is increased in humans with MetS and with reduced HDL-C levels [29]. DAG is also part of the pathway leading to the translocation of GLUT4 to the cell membrane [30].

Among the 5 tSNPs selected, only one was significantly associated with obesity-related metabolic complications in the present study. More specifically, rare allele homozygotes of rs10838599 had significantly higher plasma HDL-C levels compared to heterozygotes and homozygotes for the wild-type allele. To our knowledge, this study is the first to report an association of DGKZ gene variations and cardiometabolic risk factors in severely obese patients. Rs10838599 wild-type allele was associated with decreased plasma HDL-C levels, a lipid disturbance associated with an alteration of free fatty acid metabolism [16,31]. Low plasma HDL-C level is a criterion of MetS and an important indicator of increased cardiovascular risk [1,32]. A reduction of HDL-C concentration of 10% increases the risk of CVD by 13% [33,34]. HDL-C has been shown, with LDL-C and in combination with other growth factors, to stimulate cell proliferation of both endothelial and smooth muscle cells [35]. The capacity of HDL-C to activate protein kinase C (PKC) has also been reported [36], this enzyme being one of the four main proteins having an activity modulated by DAG [37,38]. Inversely, PKC regulates levels and activity of the DGK-zeta protein in the nucleus [25]. It possibly suggests that HDL-C plays a role in transduction of the signal cascade, in the regulation of cellular events [39], and could be linked to DAG metabolism. Besides, PA has deleterious effects on HDL-C functionality [40]. PA is in low abundance in HDL-C particles known for their heterogeneous composition and structure [41].

No significant association was found between the phenotype-associated SNP rs10838599 and methylation levels. There was also no correlation between gene methylation and expression levels. These results suggest that DNA sequence and methylation variations are not significant contributors to DGKZ expression levels, at least in severely obese patients. As DAG is a major regulator of many important biological processes, its abundance in cells might be tightly controlled, which is a key role for DGK [9]. Expression of the gene coding for DGK isoforms such as DGKZ, may be closely regulated accordingly by epigenetic mechanisms, influenced by environmental factors such as diet [42]. Through excessive intake in food, an accumulation of TG can lead to fat storage in ectopic sites, which is associated with dyslipidemia and other obesity-related metabolic complications such as hypertension, glucose intolerance [43], and  insulin resistance [44]. The lipid profile measured in plasma has been shown to be an indicator of both genetic and environmental backgrounds [29]. Moreover, the expression of DGKZ has been described to be stimulated in hypothalamus by a high-fat diet [45]. Thus, this suggests that an excess of fat intake may influence the expression level of DGK genes. Besides, DAG might also be playing a role in hypertension pathophysiology [46], but the current study only revealed a trend between two of DGKZ SNPs, rs10838599 and rs3740976, and higher DBP for rare allele carriers. Nevertheless, altered lipid metabolism and  hypertension are associated [47].

Some limitations need to be taken into account when interpreting the results presented in this study. First, a relatively small number of obese participants were available to measure gene methylation and expression levels in VAT. The limitation to 32 subjects used for methylation analysis is due to the highly stringent selection criteria. For expression analysis, the number of individuals included in the subgroup was limited to 14 because this part of methodology was retrieved from a previous study [11]. However, identifying the potential mechanism relating phenotype-associated SNPs to MetS was, however, a secondary objective. Also, the effects of DGKZ SNPs may have been altered by the systemic inflammation generated by obesity [48], so more studies with different study cohorts are needed.

Conclusion

With this study, we found that DGKZ rs10838599 is associated with plasma HDL-C thereby suggesting that it may contribute to the inter-individual variability observed in the cardiometabolic risk profile of patients with severe obesity. However, methylation levels of CpG sites were not associated to this gene variation.

Authors’ Contributions

SBégin wrote the article; SBégin and FG completed the statistical analyses; MCV, AT, YD and LP established study design; FSH, SL and PM recruited patients, collected clinical data and samples; SBégin and MCV have principal liability for final content. All authors read and approved the final manuscript.

Acknowledgements

We acknowledge the invaluable collaboration of surgeons Simon Biron, Odette Lescelleur, Laurent Biertho and Simon Marceau of the IUCPQ bariatric surgery team, who have sampled adipose tissue for this project and contributed to the organization of the clinical setup needed to perform our studies. We express our gratitude to Christine Racine and Caroline Gagnon from the IUCPQ. We also thank Alain Houde and Catherine Raymond for technical assistance and sample management. MCV is Canada Research Chair in Genomics Applied to Nutrition and Health recipient.

Funding Information

This work was supported by a grant-in-aid from Heart and Stroke Foundation of Canada.

AT receives research funding from Johnson & Johnson Medical Companies for studies unrelated to the present publication.

Conflict of Interests

The authors declare that they have no competing interests.

References

  1. Hoffman EL, VonWald T, Hansen K (2015) The metabolic syndrome.  S D Med Spec No: 24-28. [Crossref]
  2. Hamer M, Stamatakis E (2012) Metabolically healthy obesity and risk of all-cause and cardiovascular disease mortality.  J Clin Endocrinol Metab 97: 2482-2488. [Crossref]
  3. Ogden CL, Yanovski SZ, Carroll MD, Flegal KM (2007) The epidemiology of obesity.  Gastroenterology 132: 2087-2102. [Crossref]
  4. Wang J, Wu Z, Li D, Li N, Dindot SV, et al. (2012) Nutrition, epigenetics, and metabolic syndrome.  Antioxid Redox Signal 17: 282-301. [Crossref]
  5. Martínez JA, Milagro FI, Claycombe KJ, Schalinske KL (2014) Epigenetics in adipose tissue, obesity, weight loss, and diabetes.  Adv Nutr 5: 71-81. [Crossref]
  6. Eckel RH, Grundy SM, Zimmet PZ (2005) The metabolic syndrome.  Lancet 365: 1415-1428. [Crossref]
  7. Florentin M, Liberopoulos EN, Wierzbicki AS, Mikhailidis DP (2008) Multiple actions of high-density lipoprotein.  Curr Opin Cardiol 23: 370-378. [Crossref]
  8. Sakane F, Imai S, Kai M, Yasuda S, Kanoh H (2007) Diacylglycerol kinases: why so many of them?  Biochim Biophys Acta 1771: 793-806. [Crossref]
  9. Cai J, Abramovici H, Gee SH, Topham MK (2009) Diacylglycerol kinases as sources of phosphatidic acid.  Biochim Biophys Acta 1791: 942-948. [Crossref]
  10. Rincón E, Gharbi SI, Santos-Mendoza T, Mérida I (2012) Diacylglycerol kinase ζ: at the crossroads of lipid signaling and protein complex organization.  Prog Lipid Res 51: 1-10. [Crossref]
  11. Bouchard L, Tchernof A, Deshaies Y, Marceau S, Lescelleur O, et al. (2007) ZFP36: a promising candidate gene for obesity-related metabolic complications identified by converging genomics. Obes Surg 17: 372-382. [Crossref]
  12. Guénard F, Tchernof A, Deshaies Y, Pérusse L, Biron S, et al. (2014) Differential methylation in visceral adipose tissue of obese men discordant for metabolic disturbances.  Physiol Genomics 46: 216-222. [Crossref]
  13. Barros SP, Offenbacher S (2009) Epigenetics: connecting environment and genotype to phenotype and disease.  J Dent Res 88: 400-408. [Crossref]
  14. Vohl MC, Sladek R, Robitaille J, Gurd S, Marceau P, et al. (2004) A survey of genes differentially expressed in subcutaneous and visceral adipose tissue in men.  Obes Res 12: 1217-1222. [Crossref]
  15. Robitaille J, Després JP, Pérusse L, Vohl MC (2003) The PPAR-gamma P12A polymorphism modulates the relationship between dietary fat intake and components of the metabolic syndrome: results from the Québec Family Study. Clin Genet 63: 109-116. [Crossref]
  16. Expert Panel on Detection, Evaluation, and Treatment of High Blood Cholesterol in Adults (2001) Executive Summary of The Third Report of The National Cholesterol Education Program (NCEP) Expert Panel on Detection, Evaluation, And Treatment of High Blood Cholesterol In Adults (Adult Treatment Panel III). JAMA 285: 2486-2497. [Crossref]
  17. Bakker MI, May L, Hatta M, Kwenang A, Klatser PR, et al. (2005) Genetic, household and spatial clustering of leprosy on an island in Indonesia: a population-based study.  BMC Med Genet 6: 40. [Crossref]
  18. Goto K, Hozumi Y, Nakano T, Saino-Saito S, Martelli AM (2008) Lipid messenger, diacylglycerol, and its regulator, diacylglycerol kinase, in cells, organs, and animals: history and perspective.  Tohoku J Exp Med 214: 199-212. [Crossref]
  19. Topham MK, Bunting M, Zimmerman GA, McIntyre TM, Blackshear PJ, et al. (1998) Protein kinase C regulates the nuclear localization of diacylglycerol kinase-zeta.  Nature 394: 697-700. [Crossref]
  20. Cai K, Mulatz K, Ard R, Nguyen T, Gee SH1 (2014) Increased diacylglycerol kinase ζ expression in human metastatic colon cancer cells augments Rho GTPase activity and contributes to enhanced invasion.  BMC Cancer 14: 208. [Crossref]
  21. Samuel VT, Shulman GI (2012) Mechanisms for insulin resistance: common threads and missing links.  Cell 148: 852-871. [Crossref]
  22. Mannerås-Holm L1, Kirchner H1, Björnholm M1, Chibalin AV1, Zierath JR2 (2015) mRNA expression of diacylglycerol kinase isoforms in insulin-sensitive tissues: effects of obesity and insulin resistance.  Physiol Rep 3. [Crossref]
  23. Bunting M, Tang W, Zimmerman GA, McIntyre TM, Prescott SM (1996) Molecular cloning and characterization of a novel human diacylglycerol kinase zeta.  J Biol Chem 271: 10230-10236. [Crossref]
  24. Topham MK, Epand RM (2009) Mammalian diacylglycerol kinases: molecular interactions and biological functions of selected isoforms. Biochim Biophys Acta 1790: 416-424. [Crossref]
  25. Los AP, Vinke FP, de Widt J, Topham MK, van Blitterswijk WJ, et al. (2006) The retinoblastoma family proteins bind to and activate diacylglycerol kinase zeta.  J Biol Chem 281: 858-866. [Crossref]
  26. Chibalin AV, Leng Y, Vieira E, Krook A, Björnholm M, et al. (2008) Downregulation of diacylglycerol kinase delta contributes to hyperglycemia-induced insulin resistance.  Cell 132: 375-386. [Crossref]
  27. Mikhail N (2009) The metabolic syndrome: insulin resistance.  Curr Hypertens Rep 11: 156-158. [Crossref]
  28. Zhang C, Klett EL, Coleman RA (2013) Lipid signals and insulin resistance.  Clin Lipidol 8: 659-667. [Crossref]
  29. Jové M, Naudí A, Portero-Otin M, Cabré R, Rovira-Llopis S, et al. (2014) Plasma lipidomics discloses metabolic syndrome with a specific HDL phenotype.  FASEB J 28: 5163-5171. [Crossref]
  30. Tsuchiya A1, Kanno T, Nishizaki T (2013) Diacylglycerol promotes GLUT4 translocation to the cell surface in a PKCε-dependent and PKCλ/1 and -ζ-independent manner.  Life Sci 93: 240-246. [Crossref]
  31. Meikle PJ, Christopher MJ (2011) Lipidomics is providing new insight into the metabolic syndrome and its sequelae.  Curr Opin Lipidol 22: 210-215. [Crossref]
  32. Assmann G, Schulte H, von Eckardstein A, Huang Y (1996) High-density lipoprotein cholesterol as a predictor of coronary heart disease risk. The PROCAM experience and pathophysiological implications for reverse cholesterol transport. Atherosclerosis 124 Suppl: S11-20. [Crossref]
  33. Gordon T, Castelli WP, Hjortland MC, Kannel WB, Dawber TR (1977) High density lipoprotein as a protective factor against coronary heart disease. The Framingham Study.  Am J Med 62: 707-714. [Crossref]
  34. Després JP, Lemieux I, Dagenais GR, Cantin B, Lamarche B (2000) HDL-cholesterol as a marker of coronary heart disease risk: the Québec cardiovascular study.  Atherosclerosis 153: 263-272. [Crossref]
  35. Chen JK, Hoshi H, McClure DB, McKeehan WL (1986) Role of lipoproteins in growth of human adult arterial endothelial and smooth muscle cells in low lipoprotein-deficient serum.  J Cell Physiol 129: 207-214. [Crossref]
  36. Pörn MI, Akerman KE, Slotte JP (1991) High-density lipoproteins induce a rapid and transient release of Ca2+ in cultured fibroblasts.  Biochem J 279 : 29-33. [Crossref]
  37. Newton AC (1995) Protein kinase C: structure, function, and regulation.  J Biol Chem 270: 28495-28498. [Crossref]
  38. Brose N, Betz A, Wegmeyer H (2004) Divergent and convergent signaling by the diacylglycerol second messenger pathway in mammals.  Curr Opin Neurobiol 14: 328-340. [Crossref]
  39. Möllers C, Drobnik W, Resink T, Schmitz G (1995) High-density lipoprotein and low-density lipoprotein-mediated signal transduction in cultured human skin fibroblasts.  Cell Signal 7: 695-707. [Crossref]
  40. Rached F, Lhomme M, Camont L, Gomes F, Dauteuille C, et al. (2015) Defective functionality of small, dense HDL3 subpopulations in ST segment elevation myocardial infarction: Relevance of enrichment in lysophosphatidylcholine, phosphatidic acid and serum amyloid A. Biochim Biophys Acta 1851: 1254-1261. [Crossref]
  41. Camont L, Lhomme M, Rached F, Le Goff W, Nègre-Salvayre A, et al. (2013) Small, dense high-density lipoprotein-3 particles are enriched in negatively charged phospholipids: relevance to cellular cholesterol efflux, antioxidative, antithrombotic, anti-inflammatory, and antiapoptotic functionalities. Arterioscler Thromb Vasc Biol 33: 2715-2723. [Crossref]
  42. Chaudhary N, Nakka KK, Maulik N, Chattopadhyay S (2012) Epigenetic manifestation of metabolic syndrome and dietary management.  Antioxid Redox Signal 17: 254-281. [Crossref]
  43. Marcelin G, Chua S Jr (2010) Contributions of adipocyte lipid metabolism to body fat content and implications for the treatment of obesity.  Curr Opin Pharmacol 10: 588-593. [Crossref]
  44. Savage DB, Petersen KF, Shulman GI (2007) Disordered lipid metabolism and the pathogenesis of insulin resistance.  Physiol Rev 87: 507-520. [Crossref]
  45. Liu Z, Chang GQ, Leibowitz SF (2001) Diacylglycerol kinase zeta in hypothalamus interacts with long form leptin receptor. Relation to dietary fat and body weight regulation. J Biol Chem 276: 5900-5907. [Crossref]
  46. Kulkarni H, Meikle PJ, Mamtani M, Weir JM, Barlow CK, et al. (2013) Plasma lipidomic profile signature of hypertension in Mexican American families: specific role of diacylglycerols.  Hypertension 62: 621-626. [Crossref]
  47. Ferrari P, Weidmann P, Shaw S, Giachino D, Riesen W, et al. (1991) Altered insulin sensitivity, hyperinsulinemia, and dyslipidemia in individuals with a hypertensive parent.  Am J Med 91: 589-596. [Crossref]
  48. Hotamisligil GS (2006) Inflammation and metabolic disorders.  Nature 444: 860-867. [Crossref]

Editorial Information

Editor-in-Chief

Sharma S Prabhakar
Texas Tech University Health Sciences Center

Article Type

Research Article

Publication history

Received: April 11, 2016
Accepted: May 09, 2016
Published: May 12, 2016

Copyright

©2016 Bégin S. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Citation

Bégin S, Guénard F, Tchernof A, Deshaies Y, Pérusse L, et al. (2016) The rare allele of DGKZ SNP rs10838599 is associated with variability in HDL-cholesterol levels among severely obese patients. Integr Obesity Diabetes. 2: doi: 10.15761/IOD.1000148

Corresponding author

Marie-Claude Vohl

Institute of Nutrition and Functional Foods (INAF) 2440 Hochelaga Blvd, Qc, Canada, G1V 0A6, Tel: (418)656-2131 ext. 4676, Fax: (418)-656-5877.

E-mail : marie-claude.vohl@fsaa.ulaval.ca

Table 1: Subjects’ characteristics

All

MetS+

MetS-

P-value

Number of Subjects (% male)

1752 (31.1)

1428 (33.6)

317 (20.2)

-

Age (years)

43.02 ± 10.62

43.96 ± 10.58

38.66 ± 9.75

<0.0001

BMI (kg/m2)

51.8 ± 8.9

52.0 ± 9.0

51.2 ± 8.7

0.32

Waist girth (cm)

140.6 ± 17.7

141.7 ± 17.4

135.5 ± 18.5

<0.0001

CRP (mg/L)

11.04 ± 9.12

11.22 ± 9.16

10.42 ± 9.00

0.05

Fasting glucose (mmol/L)

6.53 ± 2.32

6.85 ± 2.41

5.11 ± 0.98

<0.0001

Lipid profile (mmol/L)

    TG

1.83 ± 1.08

1.98 ± 1.12

1.17 ± 0.37

<0.0001

    Total-C

4.70 ± 0.95

4.68 ± 0.98

4.80 ± 0.82

0.05

    LDL-C

2.67 ± 0.83

2.65 ± 0.84

2.80 ± 0.75

0.02

    HDL-C

1.24 ± 0.35

1.18 ± 0.32

1.49 ± 0.37

<0.0001

    Total-C/HDL-C

4.03 ± 1.29

4.17 ± 1.34

3.33 ± 0.74

<0.0001

Blood Pressure (mmHg)

     SBP

139.0 ± 17.1

140.0 ± 17.3

134.6 ± 15.6

0.0006

     DBP

83.8 ± 11.5

84.0 ± 11.7

83.1 ± 10.3

0.04

Boldface values represent statistically significant differences (P ≤ 0.05). Values are presented as mean ± SD. Abbreviations: MetS: Metabolic Syndrome; BMI: Body Mass Index; CRP: C-Reactive Protein; TG: Triglycerides; Total-C: Total Cholesterol; LDL-C: Low-Density Lipoprotein Cholesterol; HDL-C: High-Density Lipoprotein Cholesterol; SBP: Systolic Blood Pressure; DBP: Diastolic Blood Pressure; SD: Standard Deviation.

Table 2: Genotype distribution and localization of selected DGKZ SNPs

SNPs

Number of genotypes

Common HMZ

HTZ

Rare HMZ

Localization

Other designation

Region

MAF

HWE                        P-values

rs11038866

1735

804

755

176

chr11:46355932

c.212+895C>G

Intron 1

0.32

0.95

rs3809047

1704

1384

309

11

chr11:46367702

c.212+12665T>C

Intron 1

0.10

0.16

rs10838599

1751

864

742

145

chr11:46378708

c.213-10133G>T

Intron 1

0.29

0.42

rs34470543

1743

1743

0

0

chr11:46389250

c.369_370insC

Exon 3

1

-

rs3740976

1738

1315

396

27

chr11:46393757

c.975+37C>T

Intron 10

0.13

0.65

aPosition according to genome build 37

bReference sequence : NM_201532.

Abbreviations: SNP: Single Nucleotide Polymorphism; HMZ: Homozygote; HTZ: Heterozygote; MAF: Minor Allele Frequency; HWE: Hardy-Weinberg Equilibrium.

Table 3: Significant genotype differences identified between fasting glucose, lipid profile and blood pressure.

rs10838599

rs3740976

Means

Means

Phenotypes

Common HMZ (WT)

HTZ

Rare HMZ

P values

Common HMZ (WT)

Rare allele carrier

P values

Number of subjects (N)

858

740

144

1309

422

BMI (kg/m2)

51.9 ± 8.8

51.5 ± 8.3

51.51 ± 9.1

0.70

51.6 ± 8.6

51.9 ± 8.5

0.69

Waist girth (cm)

141.0 ± 17.9

140.1 ± 17.5

140.7 ± 18.3

0.53

140.1 ± 17.7

141.9 ± 17.7

0.14

CRP Protein (mg/L)

10.70 ± 8.21

10.79 ± 8.43

11.46 ± 8.75

0.59

10.75 ± 8.28

11.03 ± 8.61

0.54

Fasting glucose

6.51 ± 2.27

6.56 ± 2.32

6.56 ± 2.60

0.65

6.55 ± 2.30

6.53 ± 2.42

0.98

Lipid profile (mmol/L)

    TG

1.81 ± 0.97

1.84 ± 1.20

1.88 ± 0.99

0.32

1.82 ± 0.98

1.85 ± 1.30

0.48

    Total-C

4.69 ± 0.95

4.67 ± 0.92

4.86 ± 0.99

0.10

4.70 ± 0.93

4.68 ± 0.97

0.83

    LDL-C

2.67 ± 0.83

2.66 ± 0.79

2.80 ± 0.90

0.36

2.68 ± 0.82

2.66 ± 0.84

0.94

    HDL-C

1.24 ± 0.37

1.21 ± 0.30

1.31 ± 0.46

0.03

1.24 ± 0.37

1.23 ± 0.30

0.48

    Total-C/HDL-C

4.03 ± 1.41

4.03 ± 1.16

3.98 ± 1.20

0.48

4.04 ± 1.31

3.99 ± 1.23

0.45

Blood pressure (mmHg)

    SBP

138.9 ± 16.9

138.7 ± 16.8

139.3 ± 15.6

0.74

138.5 ± 16.5

139.6 ± 17.3

0.44

    DBP

83.7 ± 10.9

83.4 ± 11.3

85.7 ± 12.0

0.08

83.5 ± 11.3

84.3 ± 10.9

0.09

Boldface values represent statistically significant results (P ≤ 0.05).

aValues presented (means ± SD) are untransformed and unadjusted.

bP-values obtained are adjusted for the effect of age, sex, BMI and medication use except for BMI and waist girth which were adjusted for age and sex.

Abbreviations: HMZ: Homozygote; WT: Wild-Type; HTZ: Heterozygote; N: Number; BMI: Body Mass Index; CRP: C-Reactive Protein; TG: Triglycerides; Total-C: Total Cholesterol; LDL-C: Low-Density Lipoprotein Cholesterol; HDL-C: High-Density Lipoprotein Cholesterol; SBP: Systolic Blood Pressure; DBP: Diastolic Blood Pressure; SD: Standard Deviation

Table 4: Association of phenotype-associated SNP rs10838599 with gene methylation levels

SNP

cg18908017   Promoter

cg24878090   Promoter

cg18765542   Promoter

cg21029403    Promoter

cg26650731   Promoter

cg05412031      Promoter

cg25151353   Promoter

cg07276415  Promoter

cg05214390    Promoter

cg14503180    Promoter

cg10503202    Promoter

rs10838599

0.95

0.99

0.58

0.23

0.91

0.09

0.24

0.94

0.65

0.39

0.10

SNP

cg15825186       Exon 1

cg24407308    Exon 1

cg04854089       Intron 1

cg10059056    Intron 1

cg09802018     Intron 1

cg22948808   Intron 1

cg22707438   Intron 1

cg15925492  Intron 1

cg20091215   Intron 1

cg06061966   Intron 1

cg07998213    Intron 1

rs10838599

0.27

0.36

0.38

0.44

0.94

0.87

0.48

0.93

0.34

0.70

0.86

SNP

cg05576959           Intron 1

cg13247398      Intron 1

cg20029347     Intron  1

cg21663341     Intron  1

cg03519157      Intron  1

cg20666386      Intron  1

cg07701241   Intron 1

cg00530720         Intron 1

cg15868421  Intron 1

cg09039436    Intron 1

cg03044514  Intron 1

rs10838599

0.55

0.66

0.40

0.51

0.68

0.63

0.34

0.67

0.93

0.99

0.73

SNP

cg20510601  Intron 1

cg22675791     Intron 1

cg18568067   Intron 1

cg12087471  Intron 1

cg07092212  Intron 1

cg01132064     Intron 1

cg17266233   Intron 1

cg07121644  Intron 1

cg21770622   Intron 1

cg02865822  Intron 1

cg18337963  Intron 1

rs10838599

0.87

0.44

0.92

0.20

0.28

0.45

0.47

0.74

0.52

0.68

0.54

SNP

cg02211741  Intron 1

cg26473110   Intron 1

cg14787093      Intron 1

cg20204701   Intron 1

cg26963044   Intron 1

cg24713122   Intron 2

cg12856521         Exon 3

cg00107629           Intron 3

cg24126592   Intron 3

cg08768904        Exon 4

cg16769381    Exon 5

rs10838599

0.27

0.59

0.53

0.67

0.84

0.48

0.68

0.91

0.32

0.73

0.12

SNP

cg27651070   Intron 6

cg01529068    Intron 9

cg16977596    Intron 16

cg12927726     Exon 20

cg05897087   Intron 20

cg02154997    Exon 21

cg12592409   Intron 27

cg11436767     Exon 30

cg07052627    Exon 31

cg06266097    Exon 31            (3’-UTR)

cg11203377    Exon 31     (3’-UTR)

rs10838599

0.20

0.68

0.60

0.15

0.39

0.48

0.96

0.61

0.56

0.70

0.98

P-values for associations were obtained from a subset of 14 obese subjects (7 MetS+ and 7 MetS-). CpG sites positions according to genome build 37.