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
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. |
|
2021 Copyright OAT. All rights reserv
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
- Hoffman EL, VonWald T, Hansen K (2015) The metabolic syndrome. S D Med Spec No: 24-28. [Crossref]
- 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]
- Ogden CL, Yanovski SZ, Carroll MD, Flegal KM (2007) The epidemiology of obesity. Gastroenterology 132: 2087-2102. [Crossref]
- 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]
- 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]
- Eckel RH, Grundy SM, Zimmet PZ (2005) The metabolic syndrome. Lancet 365: 1415-1428. [Crossref]
- Florentin M, Liberopoulos EN, Wierzbicki AS, Mikhailidis DP (2008) Multiple actions of high-density lipoprotein. Curr Opin Cardiol 23: 370-378. [Crossref]
- 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]
- Cai J, Abramovici H, Gee SH, Topham MK (2009) Diacylglycerol kinases as sources of phosphatidic acid. Biochim Biophys Acta 1791: 942-948. [Crossref]
- 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]
- 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]
- 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]
- Barros SP, Offenbacher S (2009) Epigenetics: connecting environment and genotype to phenotype and disease. J Dent Res 88: 400-408. [Crossref]
- 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]
- 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]
- 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]
- 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]
- 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]
- 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]
- 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]
- Samuel VT, Shulman GI (2012) Mechanisms for insulin resistance: common threads and missing links. Cell 148: 852-871. [Crossref]
- 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]
- 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]
- Topham MK, Epand RM (2009) Mammalian diacylglycerol kinases: molecular interactions and biological functions of selected isoforms. Biochim Biophys Acta 1790: 416-424. [Crossref]
- 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]
- 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]
- Mikhail N (2009) The metabolic syndrome: insulin resistance. Curr Hypertens Rep 11: 156-158. [Crossref]
- Zhang C, Klett EL, Coleman RA (2013) Lipid signals and insulin resistance. Clin Lipidol 8: 659-667. [Crossref]
- 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]
- 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]
- Meikle PJ, Christopher MJ (2011) Lipidomics is providing new insight into the metabolic syndrome and its sequelae. Curr Opin Lipidol 22: 210-215. [Crossref]
- 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]
- 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]
- 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]
- 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]
- 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]
- Newton AC (1995) Protein kinase C: structure, function, and regulation. J Biol Chem 270: 28495-28498. [Crossref]
- 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]
- 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]
- 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]
- 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]
- Chaudhary N, Nakka KK, Maulik N, Chattopadhyay S (2012) Epigenetic manifestation of metabolic syndrome and dietary management. Antioxid Redox Signal 17: 254-281. [Crossref]
- 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]
- Savage DB, Petersen KF, Shulman GI (2007) Disordered lipid metabolism and the pathogenesis of insulin resistance. Physiol Rev 87: 507-520. [Crossref]
- 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]
- 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]
- 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]
- Hotamisligil GS (2006) Inflammation and metabolic disorders. Nature 444: 860-867. [Crossref]