Fatty acid ethanolamides (FAE), also referred to as N-acylethanolamines, are a group of endogenous ethanolamides of different fatty acids (FA) that were first identified in the late 1970s( Reference Epps, Natarajan and Schmid 1 ). N-arachidonoylethanolamine (AEA, also called anandamide), the derivative of arachidonic acid (AA), is the first isolated and identified FAE compound( Reference Devane, Hanus and Breuer 2 ) and has been characterised in placenta and fetal membranes as well as human plasma, amniotic fluid( Reference Marczylo, Lam and Amoako 3 ) and the nervous system( Reference Pasquariello, Catanzaro and Marzano 4 ). AEA serves as an endogenous ligand of cannabinoid receptors (CB)( Reference Pacher, Batkai and Kunos 5 ). Oleoylethanolamide (OEA), a derivative of oleic acid (OA), is thought to regulate satiety and body weight( Reference Romano, Coccurello and Giacovazzo 6 ) by activating PPAR-α, which is responsible for energy expenditure and energy intake through lipolysis in adipocytes( Reference Fu, Gaetani and Oveisi 7 ). Consequently, administration of OEA may become part of the treatment of eating disorders and body weight maintenance( Reference Lo Verme, Gaetani and Fu 8 ). Other FAE, including palmitoylethanolamide (PEA) and linoleoylethanolamide (LEA) over a range of concentrations are believed to have anti-inflammatory properties( Reference Hoareau, Buyse and Festy 9 – Reference O’Sullivan 11 ). Little is known, however, about docosahexaenoylethanolamide (DHEA; derivative of DHA) or α-linolenoylethanolamide (ALEA; derivative of α-linolenic acid (ALA)), although it was observed that a high dietary intake of DHA increased brain DHEA levels in piglets, suggesting a similar metabolism of AEA in the nervous system( Reference Berger, Crozier and Bisogno 12 ). Recently, it has been reported that DHEA-dependent pathways may be associated with hippocampal neurodevelopment and synaptic activity( Reference Kim, Moon and Cao 13 ). However, the biological importance of n-3-derived FAE remains to be resolved.
The endogenous FAE are involved in numerous biological activities and primarily modulated by their biosynthesis and degradation( Reference Rahman, Tsuboi and Uyama 14 ). Generally, the enzymes N-acylphosphatidylethanolamine-hydrolysing phospholipase D (NAPE-PLD) and fatty acid amide hydrolase (FAAH) appear to be responsible for the synthesis and hydrolysis of FAE, respectively, in turn influencing circulating FAE levels in vivo ( Reference Tsuboi, Ikematsu and Uyama 15 , Reference Thabuis, Tissot-Favre and Bezelgues 16 ). With the growing knowledge of nutrigenomics, the common genetic variants of NAPE-PLD and FAAH genes may modulate circulating FAE levels in humans. Investigations of common SNP in humans have demonstrated that major allele carriers of SNP rs324420 of FAAH are more likely associated with overweight/obesity and metabolic disorders( Reference Sipe, Waalen and Gerber 17 , Reference Monteleone, Tortorella and Martiadis 18 ), although the functional significance of these differences remain unclear over time( Reference Jensen, Andreasen and Andersen 19 , Reference Monteleone, Milano and Petrella 20 ). It was proposed that this polymorphism might impact specific neural mechanisms through CB signalling which, in turn, increase the risk of reward deficiency syndrome that causes obesity( Reference Hariri, Gorka and Hyde 21 ). Meanwhile, although little is known about NAPE-PLD, a study has shown that a common haplotype of the NAPE-PLD gene may be protective against obesity( Reference Wangensteen, Akselsen and Holmen 22 ). Most recently, Geurts et al.( Reference Geurts, Everard and Van Hul 23 ) proposed an important mechanism of adipose tissue NAPE-PLD on whole-body metabolism using NAPE-PLD knockout mice, indicating that NAPE-PLD may play an essential role in regulating energy homoeostasis and be responsive to cold-induced browning. In other words, alterations of the enzyme activity of NAPE-PLD may be beneficial as an anti-obesity treatment.
In addition, evidence has demonstrated that intake of high-MUFA diets elicits weight loss and/or body fat mass reduction compared with consuming high-SFA diets in humans( Reference Piers, Walker and Stoney 24 – Reference Krishnan and Cooper 27 ). In addition, oral administration of OEA appeared to beneficially affect health, resulting in weight loss in humans( Reference Thabuis, Tissot-Favre and Bezelgues 16 , Reference Naughton, Mathai and Hryciw 28 ). Therefore, a knowledge gap exists concerning whether endogenous OEA, converted from dietary OA, plays a role on body fat distribution. In addition, it remains to be determined whether the consumption of FA other than oleic acid can contribute to changes in body fat composition via their corresponding FAE.
Therefore, the purpose of the present study was to answer the following questions; first, how major FAE shift in response to dietary oil treatments in humans; second, whether common genetic variants in FAAH and NAPE-PLD genes affect plasma FAE levels; third, whether possible diet–gene interactions exist in populations with risk of metabolic syndrome; and, finally, whether associations exist between changes in circulating FAE concentrations, especially OEA, and body fat composition after the dietary interventions. Therefore, we hypothesised that consumption of various dietary oil treatments with differing FA composition can lead to corresponding shifts in plasma FAE levels, and the diet–gene interaction may play an important role on the levels of FAE in response to dietary treatments, resulting in potential changes in body fat composition.
Methods
Clinical design
This study was conducted as part of the Canola Oil Multicenter Intervention Trial (COMIT) study, a dietary intervention in adults with abdominal obesity and at least one criterion for metabolic syndrome( Reference Jones, Senanayake and Pu 29 ). The COMIT study was conducted at the Richardson Centre for Functional Foods and Nutraceuticals (RCFFN) at the University of Manitoba, the Institute of Nutraceuticals and Functional Foods at the Laval University and the Department of Nutritional Sciences at the Pennsylvania State University between September 2010 and March 2012, as previously described( Reference Senanayake, Pu and Jenkins 30 ). The study was registered at clinicaltrials.gov (NCT01351012). All the participants provided their written informed consent before the study started. The study was reviewed and approved by Institutional Review Boards or Committees at all participating sites.
The COMIT study investigated effects of a daily intake of 60 g of dietary oils low in SFA and high in MUFA or PUFA. A randomised cross-over study design was implemented, in which participants were randomly assigned to five novel vegetable oil treatments. Each treatment phase was 30 d in duration and was separated by a 4-week washout period. The diets were designed to maintain body weight, and energy needs for each participant were calculated using the Mifflin equation( Reference Mifflin, St Jeor and Hill 31 ). With the exception of the treatment oils, an identical 7-d rotating menu in the full-feeding diets was served across all treatments, with a fixed macronutrient composition (35 % fat, 50 % carbohydrate, 15 % protein). The five oil treatments were as follows: canola oil (Canola; 60 % OA, 20 % linoleic acid (LA), 10 % ALA), high oleic acid canola oil (CanolaOleic; 72 % OA, 15 % LA, 2 % ALA), DHA-enriched canola oil (CanolaDHA; 63 % OA, 13 % LA, 6 % DHA (from an algal oil)), a blend of corn oil/safflower oil (CornSaff; 18 % OA, 69 % LA) and a blend of flax oil/safflower oil (FlaxSaff; 18 % OA, 38 % LA, 32 % ALA) (Table 1). Three canola-based diets were rich in MUFA, whereas the two safflower oil blends were high in n-3 PUFA or n-6 PUFA. Compliance with the feeding protocol was assessed by clinical coordinators who ensured that participants consumed at least one of their two treatment beverages under supervision daily, as well as evaluated the returned meal bags of non-consumed food provided for off-site consumption. Compliance was further confirmed by measuring post-treatment plasma FA profiles.
Fatty acid composition of the five treatment oils (60 g based on 12 552 kJ diet/d (3000 kcal diet/d))
Canola, conventional canola oil; CanolaOleic, high oleic canola oil; CanolaDHA, DHA-enriched canola oil; FlaxSaff, a blend of flax oil and safflower oil; CornSaff, a blend of corn oil and safflower oil.
Subjects
A total of 170 volunteers were recruited using media advertisements. The inclusion criteria required an increased waist circumference (>94 cm for men and >80 cm for women) or at least one of the following criteria: TAG >1·7 mmol/l; HDL-cholesterol <1 mmol/l (men) or <1·3 mmol/l (women); blood pressure ≥130 mmHg (systolic) and/or ≥85 mmHg (diastolic); and glucose ≥5·5 mmol/l. Exclusion criteria included history of thyroid disease, diabetes mellitus, kidney disease and liver disease, current smokers, consuming more than two alcoholic drinks per week and taking lipid-lowering medications or supplements of any kind in the past 2 weeks.
Blood sample collection
Our intervention study was designed to compare end points after each experimental phase; therefore, the FA baseline composition profiles at each phase were not measured. Our aim was to systematically cross-compare responses across the treatment oil blends that were hypothesised to confer health benefits. After the 4-week washout to maintain their habitual diets and another 4-week treatment period, carry-over effects would be very limited, based on our previous clinical experience. Therefore, on days 29 and 30 of each phase, 12-h fasting blood samples were collected in EDTA-coated tubes and centrifuged at 3000 rpm for 20 min at 4°C. Separated plasma and leucocyte samples were transferred to labelled aliquot tubes and stored at −80°C until further analysed. After the completion of the trial, samples for the present study were shipped to the RCFFN for analysis.
Measurement of plasma fatty acid profiles
Plasma FA concentrations were determined at the end of each experimental phase. Individual plasma FA classes were extracted from plasma-EDTA aliquots using the Folch method( Reference Folch, Lees and Sloane Stanley 32 ), followed by methylation using methanolic hydrogen chloride as previously described( Reference Senanayake, Pu and Jenkins 30 ). FA methyl esters were separated on a DB-225 column (30 m×250 µm with 0·25-µm film thickness; Agilent Technologies) using an Agilent 6890N gas chromatograph equipped with a flame ionisation detector (GC-FID; Agilent Technologies). The oven was programmed from 70 to 240°C in five temperature steps (70°C for 1 min, rise of 30°C/min to 180°C, rise of 10°C/min to 200°C, rise of 2°C/min to 220°C and hold for 9.5 min, rise of 40°C/min to 240°C). Samples were run with a 20:1 split ratio, and He was used as the carrier gas with a column flow rate of 1 ml/min. Temperatures for the injector and detector were set at 280 and 300°C, respectively. Individual FA were identified by comparison with GLC 461 standard (Nu-Chek Prep). The internal FA standard heptadecanoic acid (C17:0) (Sigma-Aldrich) was introduced into original plasma samples to quantify the amount of each FA in plasma. The proportion of substrate FA converted to FA products was calculated based on the peak area relative to the total area and expressed as the percentage of total FA.
Total plasma fatty acid ethanolamide analysis
Analysis of FAE was conducted using a previously described method( Reference Lin, Yang and Jones 33 ). The six targeted FAE in plasma samples (PEA, OEA, LEA, AEA, ALEA and DHEA) were quantified according to the ratio of the known deuterated FAE in the internal standard mixture, with an expected relative standard deviation value of <5 %. Standards including PEA, OEA, LEA, AEA, ALEA, DHEA, PEA-d4, OEA-d4, LEA-d4, AEA-d8 and DHEA-d4 were purchased from Cayman Chemicals Company. ALEA-d4 was synthesised in our laboratory using α-linolenic acid chloride (Nu-Chek Prep) dissolved in methylene chloride with ethanolamide-d4 (Cambridge Isotope Laboratories)( Reference Devane, Hanus and Breuer 2 ). Relative ratios of FAE and their deuterated isotopes in the internal standard mixture were generated by serial dilution when the internal standard mixture was adjusted to the optimal concentration.
Aliquots of plasma samples at the end point of each phase were used for FAE extraction using a solid-phase extraction method( Reference Marczylo, Lam and Nallendran 34 ) with some modifications. The stability of FAE in frozen plasma samples has been previously validated in our laboratory( Reference Lin, Rideout and Yurkova 35 , Reference Jones, Lin and Gillingham 36 ). In brief, plasma samples that were mixed with internal standards on ice were filtered under gentle vacuum at a low flow rate through an activated Oasis HLB cartridge (Waters Corporation). Extracts were then washed and eluted with acetonitrile, followed by analysis on an ultra-performance liquid chromatography (UPLC) tandem Quattro micro API mass spectrometer (Waters Corporation), as published by Lin et al.( Reference Lin, Yang and Jones 33 ). A Kinetex XB-C18 column (100×2·1 mm, 1·7 µm; Phenomenex) was used for separation at a flow rate of 0·2 ml/min. Data were acquired and processed using MassLynx 4.1 (Waters Corporation).
Measures of body fat mass
Regional changes in body composition of study participants were assessed at both the initiation and termination of each dietary phase by dual-energy X-ray absorptiometry (DEXA) scanning according to the manufacturer’s recommendations (GE Healthcare; QDR-4500W; Hologic Corp.). The equipment, room supplies and room set-up were checked on a regular basis. All maintenance work including calibration was performed according to the manufacturer’s recommendations. Licensed DEXA technicians and coordinators performed all scans on each participant after each phase. Participants were required to lie down with the same posture each time at the centre of the measurement box from head to toe, face up and keep their arms at their sides during the scans. The regions of interest, including android fat mass, gynoid fat mass and the ratio of android:gynoid fat, were manually adjusted on the scan images. All data were then automatically calculated using the software Lunar Prodigy Advance enCORE and APEX System.
SNP analysis
To investigate the relationship between genetic variants of FAAH and NAPE-PLD genes and plasma FA and FAE levels, all the subjects were genotyped for two selected SNP (rs324420 for FAAH and rs12540583 for NAPE-PLD). According to the GWAS and NCBI dbSNP database, the selected SNP both belong to missense mutations with a higher minor allele frequency that may lead to significant functional changes in the corresponding enzymes. Genomic DNA was extracted from leucocyte samples using a commercial Qiagen DNeasy Blood and Tissue Kit according to the manufacturer’s instructions (Qiagen Sciences Inc.). The quality and quantity of individual genomic DNA were assessed by Thermo Scientific NanoDrop 2000 micro-volume spectrophotometer (Thermo Fisher Scientific). Amplification and detection of DNA were conducted with the 7500 Fast Real-Time PCR System (Applied Biosystems, Life Technologies) using optical-grade 96-well plates. Each well contained 5 µl of PCR reaction mixture, containing customised TaqMan SNP Genotyping Assays, TaqMan SNP genotyping Master Mix (Life Technologies Inc.) and pre-diluted DNA samples. Data were acquired by software StepOne 2.1 (Applied Biosystems, Life Technologies).
Statistical analyses
Only participants who completed all five dietary phases were included in the statistical analyses. A per-protocol approach was used to avoid the need for multiple imputations for missing data during the analyses. The power calculation was performed based on previous results on post-treatment plasma OEA levels in the clinical intervention by Jones et al.( Reference Jones, Lin and Gillingham 36 ), and indicated a power of 100 % (α=0·05) to detect significant differences in OEA levels between three MUFA-rich diets and two PUFA-rich diets. Statistical analyses were performed using SAS 9.2 (SAS Inc.). Statistical significance was determined using the adjusted Tukey test for multiple comparisons with P<0·05. The results are expressed as least square means with their standard errors, unless noted elsewhere.
The effects of dietary treatment on plasma FA profiles, plasma FAE levels and body composition were analysed using a mixed model with repeated-measures analysis of time. Treatment, age and sex were considered as fixed factors. Sequence of treatments and centre were used as random effects in the model. Tukey-adjusted P values were used to examine differences between-treatments. Pearson’s correlation analyses were conducted to test for associations between plasma FA and FAE and between plasma FAE and fat mass change.
Effects of different genotypes on FAE levels were analysed using mixed models. Treatment, sex, age, genotype, genotype×sex interaction and treatment×genotype interaction were entered as fixed effects, sequence of treatments were random effects, and clinical sites were selected as random effects. Repeated measures by five dietary phases on subjects were used to investigate the effects of treatments, genotypes and their interactions. Two individual SNP were analysed separately.
Results
Subject characteristics
A total of 130 participants (sixty male; seventy female; 76·5 % of the total randomised subjects) completed the five-phase intervention study. The mean ages of men and women were 46·5 (sd 14·2) and 47·6 (sd 14·5) years, with BMI of 29·8 (sd 4·4) and 29·0 (sd 4·2) kg/m2, respectively. Baseline characteristics of the study population have been previously published( Reference Senanayake, Pu and Jenkins 30 ) and are shown in Table 2.
Baseline characteristics of participants of the dietary intervention (Mean values with their standard errors; n 130)
BP, blood pressure.
* ANOVA was used to analyse between-sex differences in continuous variables. Significant difference: P<0·05.
Plasma fatty acid profiles
The effects of the diet interventions on plasma FA profiles are summarised in Table 3. Overall, all changes in FA composition agreed with our expectation based on the FA profile of each intervention diet. The Canola and CanolaOleic diets produced the highest levels of total MUFA (P<0·05), whereas CornSaff and FlaxSaff rich in n-6 and/or n-3 PUFA contents resulted in higher (P<0·05) levels of PUFA compared with the other three diets. The CornSaff group showed the lowest (P<0·05) ALA, EPA and DPA levels compared with the other four treatments. Post-treatment DHA concentrations were higher (P<0·0001) after the CanolaDHA treatment than after the other four treatments, whereas no differences were observed across the other four treatments, indicating good dietary compliance across all centres.
Selected plasma fatty acid profiles of participants at end points of each dietary phase (g/100 g)Footnote * (Least square mean values with their standard errors; n 130)
Canola, conventional canola oil; CanolaOleic, high oleic canola oil; CanolaDHA, DHA-enriched canola oil; FlaxSaff, a blend of flax oil and safflower oil; CornSaff, a blend of corn oil and safflower oil.a,b,c,d,e Mean values within a row with unlike superscript letters were significantly different (P<0·05).
* Mixed model with repeated measures was used to analyse between-treatment differences in continuous variables.
Total plasma fatty acid ethanolamides
The end point concentrations of OEA, AEA, PEA, LEA, ALEA and DHEA (n 121) were successfully measured and are shown in Fig. 1. Samples from nine participants were excluded from this analysis because of failures in the extraction of samples being run through the cartridges or in the FAE measurement using UPLC. Overall, plasma PEA and OEA were the two major FAE. No differences were observed in plasma PEA levels across the five dietary treatments. Three MUFA-rich treatments – Canola, CanolaDHA and CanolaOleic – resulted in higher (P<0·05) OEA levels compared with the two PUFA-rich diets – CornSaff and FlaxSaff. However, despite the similar OA content of the three canola-based diets, the CanolaDHA treatment appeared to lower the OEA levels compared with the CanolaOleic diet (P<0·05). The plasma AEA level after the CanolaOleic diet was higher (P<0·05) compared with the CanolaDHA, CornSaff and FlaxSaff diets, but no difference was observed between the CanolaOleic and Canola treatments. The two PUFA-rich diets had higher (P<0·05) LEA levels compared with the CanolaDHA and CanolaOleic treatments. The FlaxSaff diet had the highest (P<0·0001) ALEA levels and the DHEA level was substantially higher (P<0·0001) in response to the CanolaDHA diet compared with the other four diets.
Total plasma fatty acid ethanolamide levels at the end point of five dietary treatments. Values are presented in ng/ml as least squares means (n 121), with their standard errors. a,b,c Mean values within each graph indicate significant differences between-treatments (P<0·05). (A) Palmitoylethanolamide (PEA); (B) oleoylethanolamide (OEA); (C) linoleoylethanolamide (LEA); (D) arachidonoylethanolamide (AEA); (E) α-linolenoylethanolamide (ALEA); (F) docosahexaenoylethanolamide (DHEA). 
,Conventional canola oil; 
, high oleic canola oil; 
, DHA-enriched canola oil; 
, a blend of flax oil and safflower oil; 
, a blend of corn oil and safflower oil.
Correlation between plasma fatty acid ethanolamides and plasma fatty acids
Pearson’s correlation coefficients between plasma FAE and their corresponding FA are presented in Table 4. Overall, positive correlations existed between the five major FAE and their individual precursor FA (P<0·0001), whereas DHEA-DHA (r 0·52; P<0·0001) and ALEA-ALA (r 0·51; P<0·0001) showed the strongest correlations. Only AEA failed to produce a significant association with AA (r 0·04; P=0·3547).
Pearson’s correlation between fatty acid ethanolamides (FAE) and their corresponding precursor fatty acids in plasma (n 121)
PA, palmitic acid; PEA, palmitoylethanolamide; OA, oleic acid; OEA, oleoylethanolamide; LA, linoleic acid; LEA, linoleoylethanolamide; AA, arachidonic acid; AEA, arachidonoylethanolamine; ALA, α-linolenic acid; ALEA, α-linolenoylethanolamide; DHEA, docosahexaenoylethanolamide.
Body composition and correlation with plasma fatty acid ethanolamides
A subgroup of twenty-seven (male n 7; female n 20) volunteers at RCFFN completed all the baseline and end point DEXA measures at each phase, and therefore their data were used for the exploratory analysis on the correlations between body composition and plasma FAE levels. The baseline characteristics of the subjects are presented in online supplementary Table S1. Total fat mass, android fat mass, gynoid fat mass and the ratio of android:gynoid fat mass at baseline, end point and changes in each dietary phase are reported in Table 5. No significant differences were observed except for the end point gynoid fat mass (P=0·0503). It was mostly attributed to higher (P<0·05) gynoid fat mass for CanolaDHA compared with CanolaOleic. Changes of gynoid fat mass from baseline tended to differ (P=0·0803), due to an increase after CanolaDHA and a decrease after CanolaOleic.
Body fat mass composition at the beginning and the end and changes of each dietary phaseFootnote * (Least square mean values with their standard errors; n 27 (7, males; 20, females) from Richardson Centre for Functional Foods and Nutraceuticals)
Canola, conventional canola oil; CanolaOleic, high oleic canola oil; CanolaDHA, DHA-enriched canola oil; FlaxSaff, a blend of flax oil and safflower oil; CornSaff, a blend of corn oil and safflower oil.
* Mixed-model ANOVA corrected by sex and post hoc Tucky’s test were used to analyse treatment effects on different fat mass profiles. Significant difference: P<0·05.
Correlations between end point FAE levels and android fat mass change (end point v. baseline) were tested in the subset (n 27) using Pearson’s correlations (Table 6). Although we noticed that this subgroup had an unbalanced sex ratio, no differences on android fat mass change were observed. Overall, negative correlations between plasma OEA (r −0·24; P=0·0049) (Fig. 2), AEA (r −0·24; P=0·0059) and LEA (r −0·20; P=0·0225) levels and android fat mass changes across five diets were observed. Moreover, similar negative correlations were observed between the OEA:DHEA ratio and android fat mass change (r −0·20; P=0·0187). No correlations were observed between ALEA, DHEA or PEA and android fat mass change.
The correlation between oleoylethanolamide (OEA) and change of android fat mass (end point v. baseline) across all five diets (n 27). A weak negative but significant correlation was observed (r −0·24; P=0·0049).
Pearson’s correlation between end point fatty acid ethanolamide (FAE) levels and android fat mass changes of each dietary treatment (n 27).
PEA, palmitoylethanolamide; OEA, oleoylethanolamide; LEA, linoleoylethanolamide; AEA, arachidonoylethanolamine; ALEA, α-linolenoylethanolamide; DHEA, docosahexaenoylethanolamide.
*Significant correlation: P<0·05.
SNP characteristics and association with fatty acid ethanolamides
Genomic DNA was genotyped for 129 of 130 subjects (one blood sample failed to provide enough DNA yield). Analyses of the two SNP rs12540583 in NAPE-PLD and rs324420 in FAAH are summarised in Table 7.
Characteristics of the selected SNP associated with metabolism of fatty acid ethanolamides (n 129).
MM, homozygous for the major allele; Mm, heterozygous; mm, homozygous for the minor allele; MAF, minor allele frequency; NAPE-PLD, N-acyl phosphatidylethanolamine phospholipase D; FAAH, fatty acid amide hydrolase.
As numbers of rare homozygotes in both polymorphisms were relatively low, the genotypes were grouped into two categories – common homozygous and minor allele carriers. Results of the NAPE-PLD polymorphism rs12540583 (Table 8(a)) showed no differences in FAE levels across all treatments, except that the C-allele carriers had a higher (P=0·0010) DHEA levels after the CanolaDHEA diet, whereas a strong treatment×genotype interaction for DHEA was also observed (P=0·0004). In general, analyses of polymorphism rs324420 in FAAH (Table 8(b)) indicated that the A-allele carriers had higher (P=0·0002) DHEA levels than the CC genotype carrier. In particular, after the Canola and CanolaDHA treatments, DHEA levels of A-allele carriers were higher (P<0·05) than the homozygous C genotype group. In addition, a strong treatment×genotype interaction for DHEA was also observed (P<0·0001). DHEA levels after the CanolaDHA diet were observed to be higher in both genotypes compared with other treatments, and were significantly different between two genotypes. No differences were observed in the corresponding plasma DHA levels between the two genotypes.
Selected plasma fatty acid ethanolamides (FAE) at the end point of five dietary treatments, by rs12540583 in N-acyl phosphatidylethanolamine phospholipase D (NAPE-PLD) and rs324420 in fatty acid amide hydrolase (FAAH)Footnote * (Least square mean values; n 120)
Canola, conventional canola oil; CanolaOleic, high oleic canola oil; CanolaDHA, DHA-enriched canola oil; FlaxSaff, a blend of flax oil and safflower oil; CornSaff, a blend of corn oil and safflower oil; PEA, palmitoylethanolamide; OEA, oleoylethanolamide; LEA, linoleoylethanolamide; AEA, arachidonoylethanolamine; ALEA, α-linolenoylethanolamide; DHEA, docosahexaenoylethanolamide.
* Significant difference: P<0·05. P value shows the differences between two genotypes within the same treatment. P interaction shows the gene-treatment interaction in the entire population.
Discussion
The present study demonstrated that circulating FAE levels reflected their precursor FA from the dietary intake, under controlled feeding conditions. Polymorphisms of rs324420 (FAAH) may play a role in modulating the circulating levels of OEA and DHEA, which might further influence the body fat mass distribution. Given the growing evidence for the importance of FAE in governing energy metabolism, the implications of these trial results could be substantial.
Accumulating evidence suggests that circulating FAE levels are influenced by genetic variants, which account for the FAE synthesis and degradation in vivo ( Reference Pecina, Martinez-Jauand and Hodgkinson 37 ). In this study, we have extended the genetic investigation by examining the differences among end point plasma FAE in response to various dietary FA. Our results for the polymorphism rs12540583 with a missense mutation in humans also indicated a very limited effect of NAPE-PLD mutations on circulating FAE levels, except for the CanolaDHA treatment. To the best of our knowledge, this is the first study that has reported an analysis of the association between this polymorphism and circulating FAE levels. One clinical study reported that a common haplotype in NAPE-PLD was protective against obesity( Reference Wangensteen, Akselsen and Holmen 22 ), but no other human studies have assessed whether NAPE-PLD deficiency influences circulating FAE. In animal models, it has been reported that FAE levels were unaltered in NAPE-PLD-deficient mice( Reference Tsuboi, Okamoto and Ikematsu 38 ), suggesting that FAE might be synthesised by both NAPE-PLD-dependent and NAPE-PLD-independent pathways based on the NAPE. Recently, Geurts et al.( Reference Geurts, Everard and Van Hul 23 ) reported approximately 60 % reduction of PEA, OEA and stearoylethanolamine levels in the adipose tissues, but not in the brain, of NAPE-PLD-deleted mice compared with wild-type mice when AEA remained unchanged in both genotypes. The important evidence on research of the NAPE-PLD gene suggests the existence of an alternative pathway of FAE synthesis( Reference Simon and Cravatt 39 ). Interestingly, the lack of positive correlation between circulating AA and AEA in the present study suggests that the similar alternative synthesis pathway of AEA may exist in humans. These data also suggest that plasma AA may not directly reflect the amount in the AEA precursor pool, because AA levels can also be affected by the level of its precursor, LA, and the levels of its derivatives, eicosanoids.
Apart from synthesis of FAE, the regulation of FAE levels also involves degradation. Enzymatic hydrolysis can degrade FAE to FA and ethanolamine by FAAH in virtually all mammalian tissues( Reference Blankman and Cravatt 40 ). FAAH is an integral membrane enzyme that can catalyse the hydrolysis of all FAE. As such, genetic variants of the FAAH gene could alter the functionality of enzyme activity. The genetic missense polymorphism rs324420 in the FAAH gene has been extensively studied as an endocannabinoid risk factor in overweight/obesity and for potential cannabinoid antagonist treatment in obese populations( Reference Sipe, Waalen and Gerber 17 ). However, the direct association between this SNP and obesity has not been evaluated. One study in overweight/obese subjects with binge eating disorder (BED) or non-BED demonstrated a positive correlation between A-allele carriers of rs324420 with overweight/obesity people, but no differences between BED and non-BED subjects( Reference Monteleone, Tortorella and Martiadis 18 ). In addition, it was reported that the polymorphism rs324420 variants may be associated with early-onset but not adult obesity( Reference Muller, Bronner and Wandolski 41 ). Our study did not show any direct association between rs324420 and BMI types (data not shown). In the present study, the A-allele carriers of rs324420 in FAAH had a higher circulating DHEA level. It suggests that the A-allele carriers may have a lower activity rate for hydrolysing DHEA, resulting in a relatively higher level of DHEA, especially in response to consumption of DHA-enriched diets. Interestingly, such diet–gene interactions were observed only for DHEA and ALEA, indicating this SNP rs324420 may only affect DHEA and ALEA levels. Although little is known about the potential bioactive compounds DHEA and ALEA, studies have reported that the beneficial functions of DHA or DHEA could be attributed to major DHEA derivatives, which both have potent anti-inflammatory and organ-protective properties( Reference Yang, Fredman and Krishnamoorthy 42 , Reference Brown, Cascio and Wahle 43 ). As such, these findings concerning the SNP rs324420 may be useful for further research designed to examine the beneficial effects of DHA consumption.
On the basis of animal studies, OEA is of pharmaceutical interest because of the effects it has on the endocannabinoid system, which could be a therapy for treating severe obesity by inhibiting appetite. Recently, a study compared the efficacy for anti-obesity treatment using orally administered OEA and rimonabant, an anorectic anti-obesity drug approved for use in Europe( Reference Romano, Coccurello and Giacovazzo 6 ). Results indicated that OEA might represent a novel alternative to cannabinoid antagonists for the control of appetite, and rimonabant might also affect food intake and share a similar mechanism. Our results demonstrate a negative correlation between circulating OEA levels and the change in android fat mass from baseline across all five diets. This finding indicates that endogenous OEA from dietary OA might play a role in reducing fat mass. It is possible that the reduction of body weight, especially gynoid fat mass, by consuming diets high in OA is due to OEA: (a) suppression of appetite or (b) increased fat oxidation, resulting in decreased food intake. However, in our study, food intake was controlled and subjects on all diets consumed the same number of energy content across phases( Reference Senanayake, Pu and Jenkins 30 ). Thus, in the present study, our observations may be due to a different mechanism(s). Nevertheless, a different study conducted by our group showed a trend for increasing food intake and body weight in response to consumption of a DHA-enriched high oleic canola oil compared with corn oil or high oleic canola oil in golden hamsters( Reference Lin, Rideout and Yurkova 35 ). In the present study, the CanolaDHA diet did not reduce android fat mass, but increased gynoid fat mass compared with the other two oleic acid-rich diets, resulting in the lowest total fat mass change. Consistent findings about the OEA:DHEA ratio change have been reported in CD1 mice( Reference Wood, Williams and Pandarinathan 44 ), whereas our results reported lower plasma OEA levels in response to the DHA-rich diet. Recent evidence showed that DHEA might be a potent mediator for neurogenic differentiation( Reference Rashid, Katakura and Kharebava 45 ), and plays a role in reduction of headache and psychological distress in humans( Reference Ramsden, Zamora and Makriyannis 46 ), although we are not aware how such bioactivity of DHEA on neurological system could influence the OEA metabolism. Therefore, our observation provides indirect evidence to support the discovery of the effects of FAE on appetite and satiety that has been reported by Piomelli et al.( Reference Rodriguez de Fonseca, Navarro and Gomez 47 ). Moreover, we speculate that the elevation of DHEA in response to dietary DHA interferes with the synthesis and function of OEA to potentially suppress appetite and food intake, which would prevent weight loss.
In the present study, the plasma FA profile reflected the dietary FA in the treatment oils. The good compliance of participants across the three clinical sites has been previously reported( Reference Senanayake, Pu and Jenkins 30 ). As expected, corresponding elevations in end point FAE levels were also observed after the dietary interventions. Our results indicate that individuals with elevated plasma FA levels tend to have higher corresponding FAE levels. Results of a clinical study by Joosten et al.( Reference Joosten, Balvers and Verhoeckx 48 ) agreed with our findings and further suggest that the NEFA level may be positively correlated with both circulating fasting or non-fasting FAE. Thus, these NEFA molecules from increased membrane phospholipid cleavage become precursors of endogenous FAE as well. It is possible that the influences of NEFA on circulating FAE are based on their mobilisation from adipose tissue and clearance in the fasted state. Nevertheless, no further evidence would support the physiological stimuli involved in the changes of NEFA levels after meal consumption. In addition, as FAE are mostly converted in the tissues and then released into the circulating pool, the enzyme activities of FAE synthesis and degradation are more crucial for circulating FAE levels( Reference Tsuboi, Ikematsu and Uyama 15 , Reference Geurts, Everard and Van Hul 23 ). As such, our important findings demonstrate that the fat types in the diet can positively influence the corresponding FAE levels in humans. These results also provide evidence supporting the conclusion that dietary FA lead to increased circulating FAE in both piglets( Reference Berger, Crozier and Bisogno 12 ) and humans( Reference Jones, Lin and Gillingham 36 , Reference Joosten, Balvers and Verhoeckx 48 ), although it remains to be established whether such changes in plasma FAE levels further alter the FAE levels in target tissues such as brain, liver and intestines.
Polyunsaturated long-chain FAE are a relatively minor group compared with saturated and monounsaturated FAE, such as PEA and OEA, which are the major FAE in plasma and tissues( Reference Schmid, Schmid and Berdyshev 49 ). Evidence shows that feeding a high-SFA diet to mice( Reference Aviello, Matias and Capasso 50 ) and rats( Reference Diep, Madsen and Holst 51 ) resulted in lower levels of PEA and OEA. Joosten et al.( Reference Joosten, Balvers and Verhoeckx 48 ) reported that the two most abundant FAE, PEA and OEA levels, were approximately 7-fold higher than AEA levels in a clinical study. Our observation concurs with the findings with similar absolute values for PEA and OEA ranging from 0·1 to 4·0 ng/ml. Owing to the intervention design in our clinical trial, we were unable to demonstrate the effect of dietary SFA on PEA and OEA. However, we noticed that although post-treatment plasma LA, LEA precursor, levels were higher than OA, OEA precursor, levels in our population, the overall plasma LEA levels were lower than OEA. The correlation coefficient between LA and LEA was less significant than that between OA and OEA. It is possible that the endogenous conversion from OA to OEA may be faster compared with LA to LEA. One possible explanation is because of the variations in the positional distribution of FA on TAG molecules. The committed step of FAE synthesis is the exchange of FA at the stereospecific numbering-1 (sn-1) position. Usually, the sn-1 position contains a SFA, whereas the sn-2 position more likely is occupied by an unsaturated FA. Accordingly, PEA levels are overall the highest among the six FAE of interest in our observation as well as in other similar studies. In mammals, the ratio of OA:LA at the sn-1 position is higher than the ratio of LA:OA, which can potentially lead to the higher synthesis rate of OEA compared with LEA. In addition, the negative correlations between android fat mass changes and FAE (OEA, LEA and AEA) indicated that these endogenous FAE levels may be involved in the regulation of food intake and energy balance as suggested by Diep et al.( Reference Diep, Madsen and Holst 51 ).
This study has strengths and limitations. First, in contrast to a few previous human studies that relied on small subject numbers on the post-treatment FAE response, this study was powered to examine a larger group of people in North America. Furthermore, to the best of our knowledge, this is the first study ever to show that DHEA may have a suppressive effect on OEA levels, potentially resulting in a failure in body fat loss. Unfortunately, we were unable to assess the baseline FAE contents in plasma, which can establish the initial FAE profiles of the participants for further comparison. Although the cross-over design of our trial clearly showed the treatment effects on different circulating FAE levels, lacking of baseline values restricted our prediction on the association between change of FAE and change of body fat mass in response to different oil treatments. Nonetheless, the correlation analysis was only performed in a small subgroup of twenty-seven subjects with an unbalanced sex ratio due to the limited design of satellite sites in such multicentre trial, and thus our attempt was not able to assess the sex effects, and the question of how the shifts of FAE influence the body fat response requires further investigation.
In summary, our study reports several novel findings that clarify our understanding of the metabolic and physiological responses associated with different dietary FA classes. Our finding on the genetic variants indicates that the polymorphism rs324420 in FAAH may influence the beneficial function of FAE in humans by altering the circulating FAE levels. Furthermore, the elevated circulating DHEA in response to dietary DHA may suppress OEA levels and, in turn, interfere with the function of OEA, which may disrupt regulation of body weight.
Acknowledgements
The authors are indebted to Haifeng Yang for optimising the performance of the UPLS and GC-FID. The authors also acknowledge Dennis Labossiere for his interpretation of the plasma FA chromatograph. The authors especially thank the personnel at all participating centres, including investigators, clinical coordinators, kitchen staff, analytical teams and study volunteers.
This study was conducted under the clinical project of COMIT, which was supported by the Canola Council of Canada, Flax Council of Canada, Agriculture and Agri-Food Canada and Dow AgroSciences, Canada Research Chairs and Western Grains Research Foundation through the Growing Forward programme of Agriculture and Agri-Food Canada. The present study was also supported by the Manitoba’s Agri-Food Research and Development Initiative (ARDI) (ID: GRW FWD ARDI#12-1186).
The authors’ contributions are as follows: S. P. and P. J. H. J. contributed to the design of the study and wrote the manuscript. S. P. and X. L. conducted the data analyses. S. P., P. E. and P. J. H. J. were responsible for interpretation of the results. D. J. A. J., P. W. C., B. L., P. M. E.-K., S. G. W. and P. J. H. J. provided supervision.
None of the authors has any conflicts of interest to declare.
Supplementary Material
For supplementary material/s referred to in this article, please visit http://dx.doi.org/doi:10.1017/S0007114515005425
