Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Profiling the Oxylipin and Endocannabinoid Metabolome by UPLC-ESI-MS/MS in Human Plasma to Monitor Postprandial Inflammation

Abstract

Bioactive lipids, including oxylipins, endocannabinoids, and related compounds may function as specific biochemical markers of certain aspects of inflammation. However, the postprandial responsiveness of these compounds is largely unknown; therefore, changes in the circulating oxylipin and endocannabinoid metabolome in response to a challenge meal were investigated at six occasions in a subject who freely modified her usual diet. The dietary change, and especially the challenge meal itself, represented a modification of precursor fatty acid status, with expectedly subtle effects on bioactive lipid levels. To detect even the slightest alteration, highly sensitive ultra-performance liquid chromatography (UPLC) coupled to electrospray ionization (ESI) tandem mass spectrometry (MS/MS) methods for bioactive lipid profiling was employed. A previously validated UPLC-ESI-MS/MS method for profiling the endocannabinoid metabolome was used, while validation of an UPLC-ESI-MS/MS method for oxylipin analysis was performed with acceptable outcomes for a majority of the parameters according to the US Food and Drug Administration guidelines for linearity (0.9938 < R2 < 0.9996), limit of detection (0.0005–2.1 pg on column), limit of quantification (0.0005–4.2 pg on column), inter- and intraday accuracy (85–115%) and precision (< 5%), recovery (40–109%) and stability (40–105%). Forty-seven of fifty-two bioactive lipids were detected in plasma samples at fasting and in the postprandial state (0.5, 1, and 3 hours after the meal). Multivariate analysis showed a significant shift of bioactive lipid profiles in the postprandial state due to inclusion of dairy products in the diet, which was in line with univariate analysis revealing seven compounds (NAGly, 9-HODE, 13-oxo-ODE, 9(10)-EpOME, 12(13)-EpOME, 20-HETE, and 11,12-DHET) that were significantly different between background diets in the postprandial state (but not at fasting). The only change in baseline levels at fasting was displayed by TXB2. Furthermore, postprandial responsiveness was detected for seven compounds (POEA, SEA, 9(10)-DiHOME, 12(13)-DiHOME, 13-oxo-ODE, 9-HODE, and 13-HODE). Hence, the data confirm that the UPLC-ESI-MS/MS method performance was sufficient to detect i) a shift, in the current case most notably in the postprandial bioactive lipid metabolome, caused by changes in diet and ii) responsiveness to a challenge meal for a subset of the oxylipin and endocannabinoid metabolome. To summarize, we have shown proof-of-concept of our UPLC-ESI-MS/MS bioactive lipid protocols for the purpose of monitoring subtle shifts, and thereby useful to address lipid-mediated postprandial inflammation.

Introduction

The postprandial state, which is the state immediately following a meal, displays a dynamic course of events in response to the food ingested. It is an important period during which transient inflammation can occur. The current hypothesis is that during the normal homeostatic postprandial response, inflammatory compounds are transiently produced, which cause no harmful effects. In the case of a maladaptive postprandial response, during the so-called “postprandial metabolic inflammation”, an exaggerated and protracted inflammation may occur, linked to an unfavorable immune response [1]. This postprandial modification of innate immunity has been suggested to be responsible for the adverse effects of certain dietary fatty acids on, for instance, cardiovascular health. In general, the ω6 polyunsaturated fatty acids (PUFA) promote postprandial inflammation, while the ω3 PUFA suppress it. But the effects of fatty acid meal composition on inflammatory events and cardiovascular health are conflicting and have mainly been studied in models of chronic exposures, whereas the significance of postprandial adaptation remains unclear [27]. This gap in knowledge highlights the need for reliable biochemical markers of lipid-mediated postprandial inflammation. Several candidates have been investigated, including cytokines and bioactive lipids [811]. The bioactive lipid family of oxylipins is a particularly compelling group of molecules for postprandial studies due to its role in endothelial inflammation and cardiovascular disease [12, 13].

The term oxylipin includes ω6 derived 20-carbon eicosanoids (e.g. prostaglandins and leukotrienes), as well as other oxidized fatty acid metabolites ranging from well-established inflammatory mediators to potential key molecules of inflammation [14]. Oxylipins have long been recognized for their involvement in inflammation through the catabolic pathways of cyclooxygenase (COX), lipooxygenase (LOX), and cytochrome P450 (CYP). Furthermore, accumulating evidences suggest that another family of bioactive lipids, the endocannabinoids, share the same catabolic pathways and may be responsible for at least some of the beneficial anti-inflammatory effects from drugs inhibiting these enzymes [15]. Endocannabinoids are also fatty acid metabolites with signaling properties, but unlike oxylipins, they are able to bind to and activate cannabinoid receptors [16]. Most studied are anandamide (AEA) and 2-arachidonyl glycerol (2-AG), which are associated, among other things, with the regulation of inflammation and appetite [15, 1719]. Other N-acylethanolamines (NAEs), monoacylglycerols (MAGs) and related compounds are characterized as “cannabimimetic” when they are able to activate cannabinoid receptors, or “entourage compounds” which are unable to activate CB receptors, but still modify the activity of true endocannabinoids by inhibiting their degradation and metabolism [2023].

Even though oxylipins and endocannabinoids are connected, they are rarely investigated together. One reason for this lack of experiments is the challenging analytical protocol. Nevertheless, simultaneous analysis of oxylipins and endocannabinoids has occasionally been performed, revealing valuable information such as changes associated with dietary fish oil in mice [2427]. Furthermore, separate investigations of oxylipins and endocannabinoids have underscored the importance of studies on the effect of dietary fatty acids on circulating bioactive lipid levels [2830]. However, the postprandial responsiveness per se of oxylipins and endocannabinoids is still poorly understood, which is necessary to elucidate if these compounds are to be used as markers for studies on postprandial inflammation. The challenge meal study design offers a suitable approach to test if a method is able to detect subtle shifts in the metabolome, as illustrated by anandamide and other NAE shifts in the postprandial state [31, 32]. Furthermore, the feasibility of a challenge meal has previously been demonstrated by the reproducibility of individual metabolic responsiveness of circulating fatty acids [33]. Therefore, the challenge meal study design was used to test the performance of ultra-performance liquid chromatography (UPLC) coupled with electrospray ionization (ESI) tandem mass spectrometry (MS/MS) methods for the analysis of a selected panel of oxylipins, endocannabinoids, cannabimimetic and entourage compounds. Thereby, we did not relate the levels to free fatty acids (like in [31]), but we extended the investigation of the postprandial state compared to [32], and added compounds beyond what has been reported previously [3133].

The endocannabinoid metabolome was analyzed using a previously validated UPLC-ESI-MS/MS method [34], while a novel UPLC-ESI-MS/MS method for analysis of the oxylipin metabolome was developed and validated for application to human plasma at fasting and in the postprandial state after a well-defined meal. The measurements were repeated for the same individual on two different diets, vegan and vegetarian, to investigate the effects of modifications of precursor fatty acid status introduced both by the background diet, and the challenge meal itself, on the circulating bioactive lipid levels. Such subtle shifts in phenotype and responsiveness are difficult to detect with markers that were developed to detect instead the major changes that occur in disease vs. health. For example, C reactive protein (CRP), the most well-known marker of systemic inflammation, differentiates between individuals with overt disease compared with healthy controls [35, 36], and changes in response to major shifts in diet [37, 38]. However, CRP does not always change in response to more subtle shifts in diet or weight loss [39, 40], and does not change in the postprandial state [4144]. It is increasingly recognized that quantification of a standardized perturbation of metabolic homeostasis induced by dietary interventions is more informative than quantification of only the homeostatic situation (at fasting) [45]. Thus, more sensitive and specific markers of inflammation that can detect more subtle differences between different phenotypes and shifts in response to meals and changes in diet are needed.

To that end, the aim of this study was to develop comprehensive methods to analyze the circulating oxylipin and endocannabinoid metabolome (totaling fifty-two metabolites) with sufficient sensitivity to detect and quantify subtle differences introduced by changes in background diet before and after a challenge meal at multiple test occasions. We applied our validated method for endocanabinoids and related compounds [34], and validated a new method for oxylipins Despite the fact that several methods dealing with quantification of oxylipins previously have been reported [24, 33, 4649], method validation when using new equipment is necessary to ensure data quality and reproducibility [24]. Furthermore, there is a great need to employ methods with improved sensitivity and specificity due to the extreme low levels that most of the bioactive lipids are found at in healthy individuals.

By studying the postprandial variability of thirty-five circulating oxylipins, and twelve endocannabinoids and similar compounds, we were able to show proof-of-concept of our methods for the purpose of monitoring lipid-mediated postprandial inflammation. Our results highlighted that a well-defined challenge meal may reduce noise in the postprandial bioactive lipid metabolome compared to the fasting state, and thereby more differences in relation to subtle effects (here represented by changes in background diet) may be detected in the postprandial state than at fasting.

Materials and Methods

Study design

This was a pilot study using the challenge meal study design of a single female subject continuously eating the diet of her choice, who freely changed diet from vegan to vegetarian. From this point the vegan diet will be denominated as “usual diet” and vegetarian diet as “modified diet” since we are not making any attempt to relate the actual changes in diet to biological markers, neither on individual level nor on population level. The subject was not fed anything or told to eat anything specific as part of this study but instead was a free-living individual. She was a healthy female with no disease diagnoses, except for celiac disease, with a body mass index of 20, 29 years of age when the study began, and 30 when it ended. Body weight was monitored throughout the study and fluctuated between 55.5 to 57 kg. She was tested three times on each background diet using a challenge meal of bananas (in agreement with her dietary choice according to her diet records). The study protocol was defined according to the human research ethics in the Declaration of Helsinki and did not require anything beyond that the subject was able to donate blood safely since she was recruited under a methods development study protocol, which was approved by the Institutional Review Board of the University of California, Davis. Written informed consent was obtained prior to study commencement.

After being on her usual diet for 5 years, refraining from meat and dairy products, the subject modified her diet by inclusion of dairy products. Other consumed food products, including abstaining from eating all gluten containing foods, remained similar as (assessed by comparison of seven-day dietary records from both dietary phases). The first set of meal challenges was performed during a period of eight days before inclusion of dairy products to the background diet. The amount of dairy products was then gradually increased during one year to a steady state, which was kept for 33 weeks until the second set of meal challenges was performed.

Dietary intake was assessed with seven-day dietary records. Nutrient information for each day was analyzed from the dietary records using the NutriHand program (Nutrihand Inc., Soraya, CA) from reference nutrient data for individual foods using the USDA National Nutrient Database for Standard Reference. Nutrient data were averaged together for all 7 days within each dietary period and then compared by two-tailed t-tests. The main weekly differences between the diets during the 33 weeks were on average: 2 L yoghurt, 1.7 eggs, 200 g cheese and 2.1 dL cream and butter per week in the modified diet. This change did not affect the overall consumption of calories, proteins, total fats or carbohydrates, but resulted in a shift in the relative composition of fatty acids (S1 Table). More specifically, during the modified diet the subject consumed lower amounts of fiber, magnesium, copper, vitamin C, total PUFAs and 18:2n6, and higher amounts of calcium, vitamin B12, total saturated fatty acids (as well as individual short and medium chain fatty acids), the monounsaturated fatty acid 16:1 and cholesterol.

The subject was tested three times on each background diet on three non-consecutive days during a time frame of up to two weeks totaling six test occasions. No major lifestyle changes occurred during the course of the study, as monitored with physical activity and health questionnaires. On each test day, the subject came to the study center after an overnight fast, and was weighed and underwent a questionnaire to exclude recent changes in diet, lifestyle, physical activity, medication or supplement use. The challenge meal was then consumed within 5–7 min and consisted of 335.6 g–348.5 g of banana providing approximately 15% of daily calories. A fasting blood sample was collected and postprandial blood draws were performed at 0.5, 1, and 3 hours after the finished meal (Fig 1) using lavender-top EDTA tubes. Whole blood was centrifuged in a tabletop ultracentrifuge for 10 minutes at room temperature at 3000 rpm within one hour of collection. Plasma (containing both unesterified and esterified fatty acids) was then collected and aliquoted and stored at –80°C until analysis.

thumbnail
Fig 1. The study design.

After an overnight fast, blood was collected and thereafter the challenge meal was consumed. Postprandial blood collection was done at 0.5, 1, and 3 h after finished meal. The meal challenge was repeated at six occasions in total. Hour of blood collection varied between experiments, illustrated by time spans.

https://doi.org/10.1371/journal.pone.0132042.g001

Chemicals and reagents

Despite the excessive number of potential oxylipins, only 50 are usually reported in human plasma [27, 47]. Gathering that information and based on the purpose of our research topic, we selected 37 oxylipins for inclusion in our method to profile the oxylipin metabolome. Furthermore, 15 NAEs, MAGs and related compounds were analyzed in order to profile the endocannabinoid metabolome. Chemical structures of all analytes are shown in S1 Fig.

The following native, internal, and recovery standards were purchased from Cayman Chemicals (Ann Arbor, MI, USA): AEA, 2-AG, O-AEA, 2-AGe, NADA, PEA, OEA, DEA, NAGly, EPEA, DHEA, POEA, LEA, SEA, 2-LG, AEA-d8, 2-AG-d8, OEA-d4, and DHEA-d4 (for analysis of NAEs, MAGs and related compounds); and PGF, PGE2, TXB2, PGD2, 5(6)-EET, 8(9)-EET, 11(12)-EET, 14(15)-EET, 5,6-DHET, 8,9-DHET, 11,12-DHET, 14,15-DHET, 9(10)-EpOME, 12(13)-EpOME, 9(10)-DiHOME, 12(13)-DiHOME, 5-HETE, 8-HETE, 9-HETE, 11-HETE, 12-HETE, 15-HETE, 20-HETE, 9-HODE, 13-HODE, 15(S)-HETrE, 12-HEPE, 17-HDoHE, 5-oxo-ETE, 12-oxo-ETE, 15-oxo-ETE, 13-oxo-ODE, LTB4, Resolvin D2, Resolvin D1, 12-[[(cyclohexylamino)carbonyl]amino]-dodecanoic acid (CUDA), 12(13)-DiHOME-d4, 12(13)-EpOME-d4, 9-HODE-d4, PGE2-d4, and TXB2-d4 (for oxylipin analysis). The oxylipins 9,10,13-TriHOME and 9,12,13-TriHOME were obtained from Larodan (Sweden, Malmö). Acetonitrile (ACN) and methanol (MeOH) were from Merck (Darmstadt, Germany). Isopropanol was from VWR PROLABO (Fontenay-sous-Bois, France). Ammonium acetate (CH3COONH4) was purchased from Scharlau Chemie (Barcelona, Spain). Acetic acid was purchased from Aldrich Chemical Company, Inc. (Milwaukee, WI, USA). Butylhydroxytoluene (BHT) was from Cayman Chemical (Ann Arbor, MI, USA) and ethylenediaminetetraacetic acid (EDTA) from Fluka Analytical, Sigma-Aldrich (Buchs, Switzerland). Glycerol was from Fischer Scientific (Loughborough, UK). All solvents and chemicals were of HPLC grade or higher. Water was purified by a Milli-Q Gradient system (Millipore, Milford, MA, USA).

Internal standards (IS).

Deuterated compounds were used as internal standards (IS) for quantification purposes and were added to samples before extraction to mimic the extraction of the endogenous compounds. Three IS were used for endocannabinoid quantification (AEA-d4, OEA-d4 and 2-AG-d8), and five for oxylipin quantification (12(13)-DiHOME-d4, 12(13)-EpOME-d4, 9-HODE-d4, PGE2-d4, and TXB2-d4). For each native compound, a suitable IS was selected based on structural similarities (S2 Table). Recovery rates of each IS were calculated by adding the recovery standard DHEA-d4 (endocannabinoids) or CUDA (oxylipins) in the last step before analysis [46].

Standard stock solutions.

Analytical quantification standards were used as ready-made standard stock solutions or as solutions prepared from solid substances and stored at -80 °C. 2-AG, O-AEA, 2-AGe, 2-LG and 5(6)-EET were prepared and stored in ACN and the other standards were prepared and stored in ethanol to yield a final stock solution concentration of 250 μg/mL (endocannabinoids) and 100 or 1000 μg/mL (oxylipins). Stock solutions of IS were prepared to reach a final concentration of 40 μg/mL (endocannabinoids) and 10 μg/mL (oxylipins).

Standard curve preparation.

Further dilutions of each stock solution were made with methanol at ten different calibration levels for quantification purposes, prepared fresh on a weekly basis for endocannabinoids and stored at -80°C for oxylipins (S3 and S4 Tables). The lowest concentration in each calibration curve was further diluted to determine limit of quantification (LOQ) and detection (LOD) at a signal to noise ratio (S/N) of 10 and 3, respectively.

Endocannabinoid quantification

Analysis of NAEs, MAGs and related compounds was done according to our previously validated method [34]. Briefly, the plasma samples (300–500 μL) were subjected to solid phase extraction (SPE) using Waters Oasis HLB cartridges (60 mg of sorbent, 30 μm particle size). A MiniVac system (Farmingdale, NY, USA) was used to evaporate the eluates from the SPE cartridges. The analytes were then reconstituted in 100 μL MeOH and 10 μL recovery standard was added before UPLC-ESI-MS/MS analysis. The UPLC system consisted of a Waters Acquity Ultra Performance equipment (Milford, MA, USA) with a binary pump, a thermostated column compartment and an autosampler. LC separation was achieved using a Waters BEH C18 column (2.1 mm x 150 mm, 1.7 μm particle size) at 60°C, with an injection volume of 10 μL. MilliQ water (A) and 10 mM CH3COONH4 in MeOH (B) were used as mobile phase following this gradient elution: 0.0–9.0 min 79% B, 9.0–9.5 min 79–90% B, 9.5–10.5 min 90% B, 10.5–14.0 min 79% B, at a flow rate of 0.4 mL/min. The autosampler temperature was maintained at 10°C. The mass analysis was done on a Waters triple quadrupole MS (Micromass Quattro Ultima) equipped with an electrospray ionization source operating in positive mode (ESI+). N2 was used as drying gas (60 L/hr) and Ar as nebulization gas (650 L/hr). Source and desolvation temperatures were 150°C and 350°C, respectively.

Oxylipin quantification

A previously published SPE protocol was adapted for extraction of oxylipins from plasma [46]. On the day of extraction, the plasma samples were thawed at room temperature and centrifuged. SPE Waters Oasis HLB cartridges (60 mg sorbent, 30 μm particle size) were first washed with 2 mL of ethyl acetate, 4 mL of methanol, and 4 mL of 95:5 v/v water/MeOH with 0.1% acetic acid (WS). A quantitative volume of plasma (300–500 μL depending on the amount supplied) was loaded onto the SPE cartridges and spiked with 10 μL of internal standard solution (50 ng/mL for 12(13)-DiHOME-d4 and 12(13)-EPOME-d4, and 25 ng/mL for 9(S)-HODE-d4, PGE2-d4 and TXB2-d4) and 10 μL antioxidant solution (0.2 mg/mL BHT/EDTA in methanol/water (1:1)). The SPE cartridges were then washed with 4 mL of WS, dried under high vacuum for about 1 minute, and eluted with 2 mL methanol and 2 mL ethyl acetate into polypropylene tubes containing 6 μL of a glycerol solution (30% in methanol). Glycerol operates as a trap solution for the analytes. Eluates were evaporated under vacuum (MiniVac system, Farmingdale, NY, USA) and residues were then reconstituted in 100 μL methanol and tubes were vortexed. The solutions were transferred to LC vials and 10 μL of the recovery standard solution was added (5 ng/mL CUDA) and UPLC-AJST-ESI-MS/MS analysis was performed immediately in randomized order. The Agilent UPLC system (Infinity 1290) was coupled to an Agilent 6490 Triple Quadrupole system equipped with the iFunnel Technology source (Agilent Technologies, Santa Clara, CA, USA). The UPLC column used was the same as described for the endocannabinoid methodology, a Waters BEH C18 column (2.1 mm x 150 mm, 2.5 μm particle size), and 10 μL injection volume was employed. Different mobile phase composition and different gradients were compared in order to achieve optimal separation for all compounds, especially between critical isomer pairs. The optimal conditions were found with 0.1% acetic acid in MilliQ water (A) and acetonitrile:isopropanol (90:10) (B) using the following gradient: 0.0–3.5 min 10–35% B, 3.5–5.5 min 40% B, 5.5–7.0 min 42%B, 7.0–9.0 min 50% B, 9.0–15.0 min 65% B, 15.0–17.0 min 75% B, 17.0–18.5 min 85% B, 18.5–19.5 min 95% B, 19.5–21 min 10% B, 21.0–25.0 min 10% B with constant flow rate of 0.3 mL/min. The autosampler temperature was kept at 10°C and the column at 40°C. The mass analysis was done in negative mode (AJST-ESI-). For each oxylipin compound, precursor ions [M-H]-, product ions and optimal collision energies were established for each MRM transition using the MassHunter Optimizer software. The Source and iFunnel Optimizer was executed for optimizing seven parameters (capillary, nozzle and nebulizer voltage, gas temperature and flow, and sheath gas temperature and flow). A compromise had to be made between the optimal conditions for each compound and the final AJST-ESI- conditions were: capillary and nozzle voltage at 4000 V and 1500 V, respectively, drying gas temperature 230°C with a gas flow of 15 L/min, sheet gas temperature 400°C with a gas flow of 11 L/min, the nebulizer gas flow was 35 psi, and iFunnel high and low pressure RF at 90 and 60 V, respectively. The dynamic MRM option was used and performed for all compounds with optimized transitions and collision energies. Integration of all peaks was manually performed using the MassHunter Workstation software.

Oxylipin method validation

The method was validated according to the US Food and Drug Administration (FDA) guidelines over three consecutive days for linearity, LOD, LOQ, inter- and intraday accuracy and precision, recovery and stability [50], in line with our previously validated method for endocannabinoid analysis [34]. However, we used the LOD and LOQ definition based on the signal to noise response, which not automatically translates to the lower limit of quantification (LLOQ) corresponding to the lowest point in the calibration curve used in the FDA guidelines.

Linearity.

At least fourteen different concentrations of each standard were prepared and analyzed in triplicate. Calibration curves were calculated by the least-squares linear regression method using the equation y = m(x)+b, where “y” is equal to the response ratios (native standard peak area/internal standard peak area), “m” is equal to the slope of the calibration curve, “x” is equal to the on column concentration of the native analyte and “b” is equal to the y-interception of the calibration curve. Equal weighting factor was used. For concentration determinations of each analyte in real samples, a 10-point calibration curve was used (S4 Table). The peak area ratio between the analyte and its corresponding internal standard (S2 Table) were used (“y”) and the concentrations (“x”) were calculated using the calibration curve equation.

Limit of detection (LOD) and Limit of Quantification (LOQ).

The limit of detection (LOD) was defined as the concentration that resulted in a peak with a signal-to-noise ratio (S/N) greater than three and the limit of quantitation (LOQ) was defined as S/N greater than 10.

Accuracy and precision.

Accuracy and precision were determined using 4 quality control (QC) samples of 750 μL 100 mM phosphate buffer saline (PBS) spiked with native standard mixtures at different concentrations (corresponding to on-column values of 96; 24; 12 and 2.9 pg/μL). On three non-consecutive days, three replicates of each QC sample were extracted by SPE as described above for plasma samples and analyzed by UPLC-AJST-ESI-MS/MS together with a complete set of calibration standards. Calibration curves obtained for each batch were used to determine the QC sample concentration.

The intraday accuracy was determined as the percent difference between the expected concentration and the mean concentration for each analytical run (n = 3). The interday accuracy was determined as the percent difference between the expected concentration and the mean concentration on the three different days (n = 9). The acceptable range of inter- and intraday accuracy was considered to be 80–120%.

The coefficient of variation (% CV) for the mean concentration was used to calculate the intra- and interday precision (n = 3 and n = 9, respectively). A CV less than 20% was considered acceptable.

Recovery.

Internal standard recovery rates were established for different matrices in triplicates by adding a recovery standard (CUDA) immediately before UPLC-MS/MS injection to account for changes in volume and instrument variability. Accordingly, the internal standard recovery was calculated by spiking, in triplicates, 750 μL of PBS (100 mM) with 10 μL of internal standard solutions at four different levels (corresponding to on-column values of 9.1, 4.5, 2.3 and 1.1 pg/μL) prior to SPE extraction. Matrix-dependent recovery was established by spiking 10 μL internal standards in a similar manner to human plasma). A 5 point calibration curve was generated by plotting each internal standard on column concentration vs IS/CUDA area to calculate the amount of each internal standard recovered through all the extraction steps, expressed as the percentage of the expected value. This procedure was applied to all the samples analyzed.

Stability.

The stability of each analyte was determined by measuring sample concentrations immediately after collection and at different time points after storage at -20 and -80°C, and after two freeze-thaw cycles. Stability of analytes in working solutions was determined by comparing fresh solutions and solutions stored at -80°C for four weeks.

Statistical analysis

Metabolite levels were calculated and expressed as mean ± SEM, using GraphPad Prism 6 (San Diego, CA, U.S.A.). To detect significant temporal changes in the oxylipin and endocannabinoid metabolome during the postprandial response, differences between compound concentrations across the time points were assessed using one-way ANOVA with post hoc Tukey’s multiple comparison test at α = 0.05 considered significant. Two-tailed and paired Student’s t-test corrected for multiple comparison with the Holm-Sidak method (α = 0.05) was used to detect significant differences of each metabolite concentration at each time point between the vegan and vegetarian background diet. Multivariate data analysis was done using SIMCA software (Version 13.0, Umetrics, Umeå, Sweden). All data was mean centered and scaled to unit variance before modeling.

Results

Endocannabinoid quantification

In the 24 samples analyzed, a total of 12 NAEs, MAGs, and related compounds were detected with average levels ranging from 0.1 to 211 nM (Table 1). POEA was found at the highest levels followed by 2-LG, 2-AG, SEA and LEA. All 15 compounds analyzed were present in all samples, except for 2-AG and 2-LG (detected in all but one), EPEA (only detected in 3 vegan samples) and 2-AGe and NADA, which were not detected at all (S5 Table).

thumbnail
Table 1. Average levels (nM) ± SEM of N-acylethanolamines (NAEs), monoacylglycerols (MAGs) and related compounds in the endocannabinoid metabolome of human plasma at fasting (baseline) and in the fasting and postprandial state (all samples) in a subject on usual (vegan) or modified (vegetarian) diet.

https://doi.org/10.1371/journal.pone.0132042.t001

Student’s t-test was applied to investigate diet-dependent differences at baseline (during fasting) and in the postprandial state. There were no significant difference at baseline, and NAGly was the only compound that showed a significant diet-dependent difference in the postprandial state (at 0.5 hours after the meal, Fig 2). To our knowledge, this is the first report of NAGly in human plasma. MRM chromatograms and mass spectra of a plasma sample and a standard solution to verify its identity are shown in Fig 3.

thumbnail
Fig 2. Comparison of the postprandial response in NAGly (eCB) plasma levels in one subject on usual (vegan) and modified (vegetarian) diet, respectively.

*p < 0.05 (adjusted for multiple comparisons).

https://doi.org/10.1371/journal.pone.0132042.g002

thumbnail
Fig 3. MRM chromatograms of (A) NAGLy in a standard solution (136 pg on column) and (B) a representative plasma sample together with corresponding MS spectra (inserts).

https://doi.org/10.1371/journal.pone.0132042.g003

One-way ANOVA revealed alterations in the postprandial POEA and SEA levels when the subject was on a vegan diet (Fig 4A), with decreased levels of POEA at 1 and 3 hours after the meal, and increased levels of SEA at 1 hour after the meal. Only POEA postprandial levels were altered when the subject was on the modified diet (Fig 4B), with decreased levels at 0.5, 1 and 3 hours after the meal.

thumbnail
Fig 4. Baseline and postprandial response levels of endocannabinoid significantly different for a subject on usual diet (A), and on modified diet (B).

Values represent the mean ± SEM (n = 3 for each diet and time point. ****p < 0.0001, ***p = 0.0002, **p = 0.005 (adjusted for multiple comparisons).

https://doi.org/10.1371/journal.pone.0132042.g004

Oxylipin quantification

Dynamic MRM with defined retention time windows instead of time segments was performed allowing the instrument to monitor transitions only during the stated time window and therefore reducing the number of concurrent transitions. An extracted MRM chromatogram of all oxylipins included in the method is shown in Fig 5. No cross-talk between channels that were used for monitoring quantification standards, internal standards and the recovery standard was observed.

thumbnail
Fig 5. MRM chromatograms of each analyte analyzed in a standard solution mixture (A) and separation of critical pairs of isomers (B).

https://doi.org/10.1371/journal.pone.0132042.g005

Separation optimization

Separation of all 37 oxylipins was achieved, including critical pairs of isomers such as PGE2 and PGD2, and Resolvin D1 and D2 with the same MRM transitions (Fig 5). The specific detection of compounds with overlapping retention times, but different transitions is illustrated in S2 Fig. Eluent B acidified with 0.1% acetic acid improved signals for the prostaglandins PGD2, PGE2, PGF, and TXB2, but impaired the signals of the other compounds. Therefore, eluent B was not acidified in the final protocol. The optimal UPLC conditions, column temperature, gradient composition and pH were selected based on the resolution and intensity of each peak. The retention time accuracy for each compound was 99.2%±0.3% (n = 10).

MS/MS optimization.

The most intense transition for each analyte was used for quantification purposes, marked in bold in the list of MRM transitions, collision energies and cell voltages (Table 2). The second transition was used as a qualifier fragment. Internal standards were assigned to each native standard according to their structural similarities and hence retention time (S2 Table), with the exception of 14(15)-EET, 11(12)-EET, 5-oxo-ETE, 8(9)-EET, and 5(6)-EET that were eluting close to 12(13)-EpOME-d4 but were assigned to 9(S)-HODE-d4, because of superior linearity of the latter.

thumbnail
Table 2. Mass spectrometry parameters for multiple reaction monitoring transitions [Retention Time (RT), Cell Accelerator Voltage (CA), Collision Energy (CE)], linearity, Limit of Quantification (LOQ), and Limit of Detection (LOD) for the analyzed compounds (Fragmentor Voltage: 380 V for all compounds).

https://doi.org/10.1371/journal.pone.0132042.t002

Strassburg et al. [47] developed a method to analyze 104 oxylipins in a single analytical run using similar equipment to ours (without the iFunnel technology). Their 104-oxylipin method includes oxylipins that are also included in our method, with the exception of Resolvin D1 and D2. They used the same transitions as the ones reported here, except for PGE2, 5,6-DHET, 12-HEPE, 13-HODE, 9-HODE, 15-HETE, 5-HETE, and 12(13)-EpOME.

Linearity, LOD and LOQ.

The linearity of the method was determined with calibration curves over a concentration range of 0.05–2200 pg on column. Regression analysis without weighing factors produced R2 values between 0.9938 and 0.9996 (Table 2). The LOD ranged between 0.0005 and 2.1 pg on column, while LOQ ranged between 0.0005 and 4.2 pg on column (Table 2). Yang et al. [46] established LOQ values between 0.03 and 15.96 pg on column. Their LOQs for 15-oxo-ETE, 9-HETE, 11,12-DHET, and LTB4 were lower compared to ours. For all other compounds their LOQ values were equal or higher than ours. A comparison of all individual LOQs can be found in S6 Table.

Wang et al. [48] recently reported a comprehensive method for analysis of eicosanoids and other bioactive lipids in human plasma (184 compounds). However, the 35 compounds in common with our method presented higher LOQ values (S6 Table). Furthermore, a combined method to quantify prostanoids and N-acylethanolamines in human plasma included a sub-set of the compounds in our method, but also with higher LOQ values [49].

Despite the advantage illustrated above of a large panel of compounds analyzed with a short run time [26, 48, 49], a compromised sensitivity is unavoidable with current MS equipments. Therefore, we decided to select a sub-set of compounds including those previously reported in human plasma and others of potential interest to postprandial inflammation.

Accuracy and precision.

Four QC samples at different concentrations were used to establish accuracy and precision of the method (S7 Table). In the current study, the sample on-column amounts were in general close to QC3 and QC4, and occasionally at the level of QC2. Precision ranged from 0.1 to 17%, and were in general 5% or below (for 3/4 of the samples). Only 1/20 of the QC samples displayed CV values of 10% or higher. Hence, conditions for a precise method were achieved for all oxylipins tested according to the FDA guidelines of CV below 20%. Furthermore, the majority of QC samples (4/5) were in the range of 85 to 115%, thereby fulfilling the FDA guidelines of accuracy (80–120%). However, the accuracy ranged from 48 to 175% at the lowest QC level for 5 compounds (9,10,13-TriHOME, PGF2α, 12(13)-EPOME, 5(6)-EET and 12-oxo-ETE).

Recovery.

The average recovery rates ranged from acceptable 71 to 109% in PBS, human, fish and pig plasma, as well as CSF (Fig 6), except for the recovery of 12(13)-EpOME-d4 in plasma, which was approximately 40%. The recovery for internal standards was considered representative of the recovery for the native compounds, thus good recovery rates for internal standards ensure well-functioning extraction and analysis of the endogenous oxylipins. It has previously been shown that the polarity of each compound influences the recovery from different matrices according to their lipophilicity and/or presence of protein [51]. Our results confirm PBS to mimic human plasma with similar recovery rates for the internal standards.

thumbnail
Fig 6. Average recovery rates for deuterated internal standards (IS) in PBS (n = 5) and human plasma (n = 24) Expressed as mean±SEM.

https://doi.org/10.1371/journal.pone.0132042.g006

Stability.

Stability tests of human plasma analyzed after 7, 31 and 300 days in -80°C showed that the oxylipins were stable over the whole period (up to 300 days) with the exception of 15-oxo-ETE (S8 Table). Two cycles of freeze-thaw (at -20°C for one week) resulted in a higher degree of degradation when compared to the longest storage time at -80°C. The stability values in human plasma for the endocannabinoids and related compounds were in line with those reported in bovine milk [34].

Analysis of changes in metabolites.

In the 24 human plasma samples analyzed, 35 oxylipins derived from 4 different PUFAs and biosynthesized via three different pathways (COX, LOX, and CYP) could be detected (Table 3). These results are in accordance to other studies, which also quantified a variety of oxylipins in human plasma [27, 47]. Twenty-six oxylipins were detected in all samples and at all-time points independent of background diet (both usual and modified diets), 20-HETE was detected in all samples on usual diet and 7/12 of modified diet samples, 8(9)-EET, 11(12)-EET, and 5-oxo-ETE were detected in all except one or two samples, while 5 oxylipins (LTB4, 5,6-DHET, 17(R)-HDoHE, 9-HETE and 14(15)-EET were detected irregularly (S9 Table). Two compounds, Resolvin D1 and 2, were not detected at all. Measurable concentrations of Resolvin D1 and 2 were, however, previously found in human plasma following ω3 fatty acid supplementation [52]. Resolvins are known to act as lipid mediators in the resolution of inflammation [53] and might therefore not be present, or only at very low concentrations, in plasma of a healthy subject.

thumbnail
Table 3. Average levels (nM)±SEM of compounds in the oxylipin metabolome of human plasma at fasting (baseline) and in the fasting and postprandial state (all samples) in a subject on her usual or modified diet.

https://doi.org/10.1371/journal.pone.0132042.t003

Oxylipins above LOQ and with S/N > 10 were included in the statistical calculations of the postprandial response. For oxylipins with S/N > 10 (but levels below LOQ), the concentration was calculated using an extrapolated calibration curve. For oxylipins with S/N between 3 and 10 (and levels below LOQ), the value of LOD was used. In total, the variability of 35 oxylipins was investigated.

Oxylipin levels ranged over three orders of magnitude from 0.01 to 76 nM, where baseline levels on usual diet ranged from 0.03 to 19 nM and modified diet baseline levels ranged from 0.03 to 76 nM (Table 3).

Of all oxylipins tested, only one compound (TXB2) displayed significantly different baseline levels between the different background diets. There were decreased TXB2 levels when the subject was on the modified diet compared to the usual background diet (Fig 7). TXB2 is an inactive COX product of arachidonic acid (AA, 20:4n6). Its precursor, TXA2, is the major COX product in blood platelets and acts as a pro-aggregator and vasoconstrictor [54]. According to Sinclair et al [55], diets rich in AA show no effect on levels of TXA2, suggesting that differences in TXB2 levels are not caused by dietary AA changes. In contrast, a connection between AA concentration in plasma and TXB2 production in phenylketonuric subjects and lower concentrations of TXB2 in serum was found [56].

thumbnail
Fig 7. Significant different baseline levels of TXB2 (p = 0.001) in plasma collected at usual (n = 3) and modified (n = 3) background diets, respectively.

Values represent the mean ± SEM.

https://doi.org/10.1371/journal.pone.0132042.g007

Beside the TXB2 decrease in baseline oxylipin levels due to a modified background diet, an additional six oxylipins showed significantly lower levels at one or more time points during the postprandial state when the subject was on the modified diet compared to the usual background diet (Fig 8). Notably, TXB2 showed no significant differences between background diets in the postprandial state, only at baseline. During the postprandial response, two linoleic acid (LA, 18:2n6) products from the LOX pathway (9-HODE and 13-oxo-ODE), showed significant differences at one hour after the meal between the usual and modified background diets. Furthermore, significantly different levels were detected in four CYP- derived oxylipins; both from the LA precursor (9(10)-EpOME and 12(13)-EpOME), and the AA precursor (20-HETE and 11,12-DHET).

thumbnail
Fig 8. Comparison of the postprandial oxylipin plasma levels in one subject on a usual and modified diet, respectively.

* indicates p < 0.05 (adjusted for multiple comparisons).

https://doi.org/10.1371/journal.pone.0132042.g008

One-way ANOVA revealed significant differences across the time points in the postprandial oxylipin response (Fig 9). Alterations in the postprandial oxylipin levels were observed only for compounds derived from LA. For the diols, biosynthesized via the CYP pathway, a significant decrease was visible in vegan samples at 1 hour (9(10)-DiHOME) and at 0.5, 1, and 3 hours (12(13)-DiHOME) after the meal compared to the baseline level (at fasting). Furthermore, the ketone 13-oxo-ODE from the LOX pathway also displayed a significant decrease at 0.5 and 1 hour after the meal compared to the baseline level in usual diet samples. Other LOX products, 9-HODE and 13-HODE (alcohols), were found to also decrease in comparison to baseline for both usual and modified diet samples at 0.5, 1, and 3 hours after the meal. No significant changes, however, were observed for COX-derived metabolites (e.g. prostaglandins).

thumbnail
Fig 9. Baseline and postprandial response levels significantly different of oxylipins for a subject on usualdiet (A), and on modified diet (B).

Values represent the mean ± SEM (n = 3 for each diet and time point). Brackets indicates significantly different time points: ****p < 0.0001, ***p = 0.0002 (adjusted for multiple comparisons).

https://doi.org/10.1371/journal.pone.0132042.g009

Profiling of the oxylipin and endocannabinoid metabolome

Coefficient of variation (CV) values for the oxylipins, NAEs, MAGs and related compounds in human plasma (S3 and S4 Figs were considerably larger than those found for quality control samples of standards spiked to PBS (S10 and S11 Tables). Hence, the variability of the plasma metabolites was attributed to biological variation rather than analytical. For the purpose of investigating if the biological variation was due to noise, or if there was a true signal from the background diet, multivariate data analysis was performed combining the information from all oxylipins found in all samples (n = 26), and NAEs, MAGs and related compounds found in all but one or two samples (n = 11). We chose a conservative approach and included only the metabolites in common between all samples (but one or two) to ensure that the same metabolic domain was captured for all samples. The relationship between systematic variability and background diet (usual and modified) was investigated using multivariate regression analysis in terms of orthogonal projections to latent structures discriminant analysis (OPLS-DA) [57]. When data from only baseline samples were included (n = 6), the result was no significant OPLS-DA model assessed by ANOVA based on the cross-validated score vectors (CV-ANOVA: 0.764) [58]. The reason might be the low number of samples, however, adding all samples to the model (n = 24), resulted also in no significant model (CV-ANOVA: 0.267). On the contrary, only using data from the samples in the postprandial state (n = 18) resulted in a significant separation between samples collected during the usual and modified background diets (Fig 10). This was in agreement with the univariate data analysis, with only one metabolite showing a significant difference at baseline (TXB2), while 7 (NAGLy, 9-HODE, 13-oxo-ODE, 9(10)-EpOME, 12(13)-EpOME, 20-HETE, and 11,12-DHET) were different between the background diets in the postprandial state. However, the majority of the metabolites showed no difference due to background diet in univariate analysis, and among the ones responsive to the challenge meal, several metabolites were in common between the background diets (POEA, 9-HODE, and 13-HODE), highlighting the additional value of using multivariate data analysis in nutritional studies, which also previously has been recommended and frequently applied [59]. Furthermore, from the multivariate analysis we could conclude that metabolites in the oxylipin metabolome were more specific to each background diet than metabolites in the endocannabinoid metabolome, since only oxylipins displayed variable importance above 1, the limit for variable influence on projection considered to be relevant for separation of classes (S12 Table). LA derived epoxides in the CYP pathway (9(10)-EpOME and 12(13)-EpOME) were characteristic of the usual diet, and could be due to the higher levels of LA in the diet during this time period. On the other hand, an eicosapentaenoic acid derived oxylipin in the LOX pathway (12(S)-HEPE) dominated when the subject was on the modifieddiet in relation to the whole oxylipin and eCB metabolome under study (S5 Fig). This was despite the fact that intake of eicosapentaenoic acid, or its precursor, alpha-linolenic acid, were not different between the two diets.

thumbnail
Fig 10. Oxylipin and endocannabinoid profiles in the postprandial state recapitulated as orthogonal scores (t[1] and to[1]) calculated by OPLS-DA.

Time after challenge meal (in hours) is found next to each sample. The subject displayed different oxylipin and endocannabinoid metabolomes on usual (grey circles) vs modified (black squares) background diet. Model assessment parameters were: 1 predictive and 1 orthogonal component, p-value calculated by CV-ANOVA: 0.033, total systematic variation among the metabolites captured by the model (R2X): 0.463, total systematic variation between the diets captured by the model (R2Y): 0.815, predictive ability of the model (Q2): 0.53.

https://doi.org/10.1371/journal.pone.0132042.g010

To summarize, our results from the multivariate analysis imply that 1) combined data on the metabolite levels is more informative than each metabolite one by one, and 2) data on the oxylipin and endocannabinoid metabolome from the postprandial state contain more information than from the fasting state alone.

Discussion

This study did not attempt to relate changes in diet (vegan to vegetarian) to biological markers. The change in background diet rather represented a venue to investigate the performance of the UPLC-ESI-MS/MS methods to detect subtle fluctuations among the bioactive lipids, in this case introduced by the difference in diets, within a single subject tested with the challenge meal study design at multiple occasions. Thereby, the investigation served as a pilot study in which the biological variability was assessed and compared to the analytical variability of the UPLC-ESI-MS/MS methods. From our multivariate analysis it was evident that inclusion of dairy products in the diet produced a shift in bioactive lipid profiles in the postprandial state (Fig 10), most notably among the oxylipins (S5 Fig and S8 Table). This shift was not present to the same extent in the fasting state alone suggesting that a subtle change in metabolite levels occurred in the postprandial state (which was successfully detected with our analytical protocols). Indeed, there was one analyte in the endocannabinoid metabolome (NAGly) and 6 in the oxylipin metabolome (20-HETE, 9,(10)-EpOME, 11,12-DHET, 9-HODE, 13-oxo-ODE and 12(13)-EpOME) that were significantly different in the postprandial state, while only one (TXB2) was different in the fasting state, and not in the postprandial state (Fig 2 and Fig 8). So, the shift seemed to be analyte-dependent, some diverged in the postprandial state, while others converged between vegan and vegetarian diets. At the same time, there were similarities in metabolite responsiveness between the usual and modified diets. Two NAEs (POEA and SEA) and 5 oxylipins (9(10)-DiHOME, 12(13)-DiHOME, 9-HODE, 13-HODE and 13-oxo-ODE) were responsive to the challenge meal when the subject was on the usual diet, of which 1 NAE (POEA) and 2 oxylipins (9-HODE and 13-HODE) also were responsive when the subject was on a vegetarian diet (Fig 4 and Fig 9). Whether the discrepancy in responsiveness is due to more variable data produced by the modified diet remains to be tested. Our data suggest that this might be the case, since higher CVs, especially at fasting, were displayed for the samples collected at modified diet compared to the usual diet (S3 and S4 Figs). But since only one subject was included in this pilot study, no statement can be made on inter-individual variations in order to assess the observed variations between the vegan and vegetarian diet.

Furthermore, the study design does not clarify if the variability was caused by true metabolic differences between the usual and modified diet, or if daily variations in diet and lifestyle that are independent of specific vegan and vegetarian aspects had an impact. We deliberately refrained from instructions that would limit lifestyle etc, in order to investigate how well our UPLC-ESI-MS/MS methods perform for assessment of the oxylipin and endocannabinoid metabolome in a free-living individual. However, to gain better understanding of dietary effects on bioactive lipid profiles, it is necessary to carry out further comprehensive studies, including additional endpoints, on more subjects. For instance, it would be valuable to follow the postprandial state over a prolonged time course (more than 3 hours) and also follow for instance the insulin levels. We selected 3 hours since it represents the lipemic peak [60]. The experiment was stopped at 3 hours since the challenge meal only provided 15% of daily calories. But with another type of challenge meal, it would be possible to prolong the experiment.

In addition to the 3 hour time point, we also measured the bioactive lipid levels in the early postprandial state, since we hypothesized that the metabolites under study would be fast responders to the meal intake. This proved to be true, since all of the postprandial responses were detected within 1 hour after the meal intake, and all except 3 out of 10 persisted until 3 hours after the meal intake (Fig 4 and Fig 9). Interestingly, all metabolites (except SEA) responded by decreasing levels, which might be a reflection of the low content of precursor fatty acids in the challenge meal (banana). Whether this is characteristic of postprandial inflammation or attributed to a normal homeostatic response remains to be tested. Our data confirms that the oxylipin and endocannabinoid metabolome are responsive to a challenge meal and that small fluctuations in their concentrations can be detected and quantified reproducibly. In addition, we found that multivariate analysis highlighted the subtle diet-dependent changes in the postprandial state. Therefore, UPLC-ESI-MS/MS analysis of these metabolites using a well-validated method can be used to investigate subtle changes expected in lipid-mediated inflammation in connection to meal intake and other clinical conditions and settings where subtle variations are important to detect.

Supporting Information

S1 Fig. Chemical structures of endocannabinoid and oxylipin analytes under study.

https://doi.org/10.1371/journal.pone.0132042.s001

(TIF)

S2 Fig. Unique quantification of analytes with overlapping retention times but different MRM transitions.

https://doi.org/10.1371/journal.pone.0132042.s002

(TIF)

S3 Fig. Coefficient of variation (CV) for oxylipins in plasma from a subject on usual (left panel) and modified (right panel) diet at baseline (fasting) and in the postprandial state (0.5, 1, and 3 hours after the challenge meal).

https://doi.org/10.1371/journal.pone.0132042.s003

(TIF)

S4 Fig. Coefficient of variations (CV) for endocannabinoids and related compounds in plasma from a subject on usual (left panel) and modified (right panel) diet at baseline (fasting) and in the postprandial state (0.5, 1, and 3 hours after the challenge meal).

https://doi.org/10.1371/journal.pone.0132042.s004

(TIF)

S5 Fig. The corresponding loading plot to the score plot in Fig 9, showing the contribution of each metabolite to the separation due to different metabolite profiles in the postprandial state of samples collected at different background diets.

Usual diet-specific metabolites are found to the left (the most typical were 9(10)-EpOME and 12(13)-EpOME, highlighted in grey), and modified diet-specific metabolites are found to the right (the most typical was 12(S)-HEPE, highlighted in grey).

https://doi.org/10.1371/journal.pone.0132042.s005

(TIF)

S1 Table. Mean values from seven-day dietary records for the subject on usual and modified background diets (assessed by two-tailed Student’s t-test, if unequal variance, p-value is for t-test assuming unequal variance or Welch's non-parametric test).

P-values below 0.05 are highlighted in bold.

https://doi.org/10.1371/journal.pone.0132042.s006

(DOCX)

S2 Table. Internal standards used for quantification.

https://doi.org/10.1371/journal.pone.0132042.s007

(DOCX)

S3 Table. Calibration standard concentrations (µg/mL) for endocannabinoids (stock solution 250 µg/mL).

https://doi.org/10.1371/journal.pone.0132042.s008

(DOCX)

S4 Table. Calibration standard concentrations (pg/mL) for oxylipins.

https://doi.org/10.1371/journal.pone.0132042.s009

(DOCX)

S5 Table. Endocannabinoid levels (nM) in the plasma samples from the challenge meal study in the fasting state (0) and postprandial state (at 0.5, 1, 3 hours after the meal) at six different occasions (a–f).

https://doi.org/10.1371/journal.pone.0132042.s010

(DOCX)

S6 Table. Comparison of limit of quantification (LOQ) values with Yang et al [44] and Wang et al [48].

https://doi.org/10.1371/journal.pone.0132042.s011

(DOCX)

S7 Table. Inter- and intraday accuracy and precision (coefficient of variation) for quality control (QC) samples at four different concentration levels.

https://doi.org/10.1371/journal.pone.0132042.s012

(DOCX)

S8 Table. Levels (nM) and stability of 5 pooled human plasma samples (expressed as % of the 0 day value) for the metabolites (oxylipins, endocannabinoids and related compounds) under study.

https://doi.org/10.1371/journal.pone.0132042.s013

(DOCX)

S9 Table. Oxylipin levels (nM) in the plasma samples from the challenge meal study in the fasting state (0) and postprandial state (at 0.5, 1, 3 hours after the meal) at six different occasions (a – f).

https://doi.org/10.1371/journal.pone.0132042.s014

(DOCX)

S10 Table. Coefficient of variation (CV) values (%) for oxylipins in the fasting state and postprandial state (at 0.5, 1, 3 hours after the meal), and in quality control (QC) samples.

https://doi.org/10.1371/journal.pone.0132042.s015

(DOCX)

S11 Table. Coefficient of variation (CV) values (%) for endocannabinoids in the fasting state and postprandial state (at 0.5, 1, 3 hours after the meal), and in quality control (QC) samples [32].

https://doi.org/10.1371/journal.pone.0132042.s016

(DOCX)

S12 Table. Variable influence on projection (VIP) values for each metabolite in the OPLS-DA model.

https://doi.org/10.1371/journal.pone.0132042.s017

(DOCX)

Acknowledgments

Monica Norlund is recognized for excellent technical assistance and the Swedish Metabolomics Centre (www.swedishmetabolomicscentre.se) is acknowledged for outstanding assistance with the UPLC-ESI-MS/MS analysis.

Author Contributions

Conceived and designed the experiments: MLN AMZ. Performed the experiments: SG-F JS. Analyzed the data: SG-F JS. Contributed reagents/materials/analysis tools: MLN AMZ. Wrote the paper: SG-F MLN.

References

  1. 1. Margioris AN. Fatty acids and postprandial inflammation. Current Opinion in Clinical Nutrition & Metabolic Care. 2009;12(2):129–37.
  2. 2. Baer DJ, Judd JT, Clevidence BA, Tracy RP. Dietary fatty acids affect plasma markers of inflammation in healthy men fed controlled diets: a randomized crossover study. The American Journal of Clinical Nutrition. 2004;79(6):969–73. pmid:15159225
  3. 3. Nicholls SJ, Lundman P, Harmer JA, Cutri B, Griffiths KA, Rye K-A, et al. Consumption of Saturated Fat Impairs the Anti-Inflammatory Properties of High-Density Lipoproteins and Endothelial Function. Journal of the American College of Cardiology. 2006;48(4):715–20. pmid:16904539
  4. 4. Poppitt SD, Keogh GF, Lithander FE, Wang Y, Mulvey TB, Chan Y-K, et al. Postprandial response of adiponectin, interleukin-6, tumor necrosis factor-α, and C-reactive protein to a high-fat dietary load. Nutrition. 2008;24(4):322–9. pmid:18262390
  5. 5. Arya F, Egger S, Colquhoun D, Sullivan D, Pal S, Egger G. Differences in postprandial inflammatory responses to a ‘modern’ v. traditional meat meal: a preliminary study. British Journal of Nutrition. 2010;104(05):724–8.
  6. 6. Masson CJ, Mensink RP. Exchanging Saturated Fatty Acids for (n-6) Polyunsaturated Fatty Acids in a Mixed Meal May Decrease Postprandial Lipemia and Markers of Inflammation and Endothelial Activity in Overweight Men. The Journal of Nutrition. 2011;141(5):816–21. pmid:21430255
  7. 7. Tulk HMF, Robinson LE. Modifying the n-6/n-3 polyunsaturated fatty acid ratio of a high–saturated fat challenge does not acutely attenuate postprandial changes in inflammatory markers in men with metabolic syndrome. Metabolism. 2009;58(12):1709–16. pmid:19625064
  8. 8. Strassburg K, Esser D, Vreeken RJ, Hankemeier T, Müller M, van Duynhoven J, et al. Postprandial fatty acid specific changes in circulating oxylipins in lean and obese men after high-fat challenge tests. Molecular Nutrition & Food Research. 2014;58(3):591–600.
  9. 9. Burdge GC, Calder PC. Plasma cytokine response during the postprandial period: a potential causal process in vascular disease? British Journal of Nutrition. 2005;93(01):3–9.
  10. 10. Manning PJ, Sutherland WHF, McGrath MM, De Jong SA, Walker RJ, Williams MJA. Postprandial Cytokine Concentrations and Meal Composition in Obese and Lean Women. Obesity. 2008;16(9):2046–52. pmid:19186329
  11. 11. Monteleone P, Piscitelli F, Scognamiglio P, Monteleone AM, Canestrelli B, Di Marzo V, et al. Hedonic Eating Is Associated with Increased Peripheral Levels of Ghrelin and the Endocannabinoid 2-Arachidonoyl-Glycerol in Healthy Humans: A Pilot Study. The Journal of Clinical Endocrinology & Metabolism. 2012;97(6):E917–E24. pmid:22442280.
  12. 12. Shearer G, Newman J. Impact of circulating esterified eicosanoids and other oxylipins on endothelial function. Curr Atheroscler Rep. 2009;11(6):403–10. pmid:19852880
  13. 13. , , , , . Distribution of soluble and microsomal epoxide hydrolase in the mouse brain and its contribution to cerebral epoxyeicosatrienoic acid metabolism. Neuroscience. 2009;163:646. pmid:19540314
  14. 14. Morisseau C, Hammock BD. Impact of Soluble Epoxide Hydrolase and Epoxyeicosanoids on Human Health. Annual Review of Pharmacology and Toxicology. 2013;53(1):37–58. pmid:23020295.
  15. 15. Fowler CJ. NSAIDs: eNdocannabinoid stimulating anti-inflammatory drugs? Trends in Pharmacological Sciences. 2012;33(9):468–73. pmid:22664342
  16. 16. Di Marzo V, Fontana A, Cadas H, Schinelli S, Cimino G, Schwartz J-C, et al. Formation and inactivation of endogenous cannabinoid anandamide in central neurons. Nature. 1994;372(6507):686–91. pmid:7990962
  17. 17. Ghafouri N, Ghafouri B, Larsson B, Turkina MV, Karlsson L, Fowler CJ, et al. High Levels of N-Palmitoylethanolamide and N-Stearoylethanolamide in Microdialysate Samples from Myalgic Trapezius Muscle in Women. PLoS ONE. 2011;6(11):e27257. pmid:22125609
  18. 18. Pintus S, Murru E, Carta G, Cordeddu L, Batetta B, Accossu S, et al. Sheep cheese naturally enriched in α-linolenic, conjugated linoleic and vaccenic acids improves the lipid profile and reduces anandamide in the plasma of hypercholesterolaemic subjects. British Journal of Nutrition. 2013;109(08):1453–62.
  19. 19. Fride E, Bregman T, Kirkham TC. Endocannabinoids and Food Intake: Newborn Suckling and Appetite Regulation in Adulthood. Experimental Biology and Medicine. 2005;230(4):225–34. pmid:15792943
  20. 20. Zoerner AA, Gutzki F- M, Batkai S, May M, Rakers C, Engeli S, et al. Quantification of endocannabinoids in biological systems by chromatography and mass spectrometry: A comprehensive review from an analytical and biological perspective. Biochimica et Biophysica Acta (BBA)—Molecular and Cell Biology of Lipids. 2011;1811(11):706–23.
  21. 21. Maccarrone M, Gasperi V, Catani MV, Diep TA, Dainese E, Hansen HS, et al. The Endocannabinoid System and Its Relevance for Nutrition. Annual Review of Nutrition. 2010;30(1):423–40. pmid:20645854.
  22. 22. Di Marzo V. ‘Endocannabinoids’ and other fatty acid derivatives with cannabimimetic properties: biochemistry and possible physiopathological relevance. Biochimica et Biophysica Acta (BBA)—Lipids and Lipid Metabolism. 1998;1392(2–3):153–75.
  23. 23. Meijerink J, Balvers M, Witkamp R. N-acyl amines of docosahexaenoic acid and other n–3 polyunsatured fatty acids – from fishy endocannabinoids to potential leads. British Journal of Pharmacology. 2013;169(4):772–83. pmid:23088259
  24. 24. Dumlao DS, Buczynski MW, Norris PC, Harkewicz R, Dennis EA. High-throughput lipidomic analysis of fatty acid derived eicosanoids and N-acylethanolamines. Biochimica et Biophysica Acta (BBA)—Molecular and Cell Biology of Lipids. 2011;1811(11):724–36.
  25. 25. Balvers MG, Verhoeckx KC, Bijlsma S, Rubingh CM, Meijerink J, Wortelboer HM, et al. Fish oil and inflammatory status alter the n-3 to n-6 balance of the endocannabinoid and oxylipin metabolomes in mouse plasma and tissues. Metabolomics: Official journal of the Metabolomic Society. 2012;8(6):1130–47. pmid:23136559; PubMed Central PMCID: PMC3483099.
  26. 26. Wong A, Sagar DR, Ortori CA, Kendall DA, Chapman V, Barrett DA. Simultaneous tissue profiling of eicosanoid and endocannabinoid lipid families in a rat model of osteoarthritis. Journal of Lipid Research. 2014;55(9):1902–13. pmid:25062663
  27. 27. Grapov D, Adams SH, Pedersen TL, Garvey WT, Newman JW. Type 2 Diabetes Associated Changes in the Plasma Non-Esterified Fatty Acids, Oxylipins and Endocannabinoids. PLoS ONE. 2012;7(11):e48852. pmid:23144998
  28. 28. Naughton SS, Mathai ML, Hryciw DH, McAinch AJ. Fatty Acid Modulation of the Endocannabinoid System and the Effect on Food Intake and Metabolism. International Journal of Endocrinology. 2013;2013:11.
  29. 29. Nording ML, Yang J, Georgi K, Hegedus Karbowski C, German JB, Weiss RH, et al. Individual Variation in Lipidomic Profiles of Healthy Subjects in Response to Omega-3 Fatty Acids. PLoS ONE. 2013;8(10):e76575. pmid:24204640
  30. 30. Schuchardt JP, Schmidt S, Kressel G, Willenberg I, Hammock BD, Hahn A, et al. Modulation of blood oxylipin levels by long-chain omega-3 fatty acid supplementation in hyper- and normolipidemic men. Prostaglandins, Leukotrienes and Essential Fatty Acids (PLEFA). 2014;90(2–3):27–37.
  31. 31. Gatta-Cherifi B, Matias I, Vallee M, Tabarin A, Marsicano G, Piazza PV, et al. Simultaneous postprandial deregulation of the orexigenic endocannabinoid anandamide and the anorexigenic peptide YY in obesity. Int J Obes. 2012;36(6):880–5.
  32. 32. Joosten M, Balvers M, Verhoeckx K, Hendriks H, Witkamp R. Plasma anandamide and other N-acylethanolamines are correlated with their corresponding free fatty acid levels under both fasting and non-fasting conditions in women. Nutrition & Metabolism. 2010;7(1):49.
  33. 33. Zivkovic A, Wiest M, Nguyen U, Nording M, Watkins S, German J. Assessing individual metabolic responsiveness to a lipid challenge using a targeted metabolomic approach. Metabolomics. 2009;5(2):209–18. pmid:20442800
  34. 34. Gouveia-Figueira S, Nording ML. Development and Validation of a Sensitive UPLC-ESI-MS/MS Method for the Simultaneous Quantification of 15 Endocannabinoids and Related Compounds in Milk and Other Biofluids. Analytical Chemistry. 2013;86(2):1186–95. pmid:24377270
  35. 35. Anuurad E, Tracy RP, Pearson TA, Kim K, Berglund L. Synergistic role of inflammation and insulin resistance as coronary artery disease risk factors in African Americans and Caucasians. Atherosclerosis. 2009;205(1):290–5. Epub 2009/01/13. doi: S0021-9150(08)00838-1 [pii] pmid:19135196; PubMed Central PMCID: PMC2700183.
  36. 36. Blackburn R, Giral P, Bruckert E, Andre JM, Gonbert S, Bernard M, et al. Elevated C-reactive protein constitutes an independent predictor of advanced carotid plaques in dyslipidemic subjects. Arterioscler Thromb Vasc Biol. 2001;21(12):1962–8. pmid:11742871.
  37. 37. Katcher HI, Legro RS, Kunselman AR, Gillies PJ, Demers LM, Bagshaw DM, et al. The effects of a whole grain-enriched hypocaloric diet on cardiovascular disease risk factors in men and women with metabolic syndrome. Am J Clin Nutr. 2008;87(1):79–90. Epub 2008/01/08. doi: 87/1/79 [pii]. pmid:18175740.
  38. 38. Neuhouser ML, Schwarz Y, Wang C, Breymeyer K, Coronado G, Wang CY, et al. A low-glycemic load diet reduces serum C-reactive protein and modestly increases adiponectin in overweight and obese adults. J Nutr. 2012;142(2):369–74. Epub 2011/12/23. pmid:22190020; PubMed Central PMCID: PMC3260063.
  39. 39. Liou YA, King DJ, Zibrik D, Innis SM. Decreasing Linoleic Acid with Constant {alpha}-Linolenic Acid in Dietary Fats Increases (n-3) Eicosapentaenoic Acid in Plasma Phospholipids in Healthy Men. J Nutr. 2007;137(4):945–52. pmid:17374659
  40. 40. Moran LJ, Noakes M, Clifton PM, Wittert GA, Belobrajdic DP, Norman RJ. C-reactive protein before and after weight loss in overweight women with and without polycystic ovary syndrome. J Clin Endocrinol Metab. 2007;92(8):2944–51. pmid:17504892.
  41. 41. Blackburn P, Despres JP, Lamarche B, Tremblay A, Bergeron J, Lemieux I, et al. Postprandial variations of plasma inflammatory markers in abdominally obese men. Obesity (Silver Spring). 2006;14(10):1747–54. pmid:17062804.
  42. 42. Payette C, Blackburn P, Lamarche B, Tremblay A, Bergeron J, Lemieux I, et al. Sex differences in postprandial plasma tumor necrosis factor-alpha, interleukin-6, and C-reactive protein concentrations. Metabolism. 2009;58(11):1593–601. Epub 2009/07/17. pmid:19604525.
  43. 43. Pietraszek A, Gregersen S, Hermansen K. Acute effects of dietary fat on inflammatory markers and gene expression in first-degree relatives of type 2 diabetes patients. Rev Diabet Stud. 2011;8(4):477–89. Epub 2011/01/01. pmid:22580729; PubMed Central PMCID: PMC3359692.
  44. 44. Poppitt SD, Keogh GF, Lithander FE, Wang Y, Mulvey TB, Chan YK, et al. Postprandial response of adiponectin, interleukin-6, tumor necrosis factor-alpha, and C-reactive protein to a high-fat dietary load. Nutrition. 2008;24(4):322–9. pmid:18262390.
  45. 45. Pellis L, van Erk M, van Ommen B, Bakker GM, Hendriks HJ, Cnubben NP, et al. Plasma metabolomics and proteomics profiling after a postprandial challenge reveal subtle diet effects on human metabolic status. Metabolomics: Official journal of the Metabolomic Society. 2012;8(2):347–59.
  46. 46. Yang J, Schmelzer K, Georgi K, Hammock BD. Quantitative Profiling Method for Oxylipin Metabolome by Liquid Chromatography Electrospray Ionization Tandem Mass Spectrometry. Analytical Chemistry. 2009;81(19):8085–93. pmid:19715299
  47. 47. Strassburg K, Huijbrechts AL, Kortekaas K, Lindeman J, Pedersen T, Dane A, et al. Quantitative profiling of oxylipins through comprehensive LC-MS/MS analysis: application in cardiac surgery. Anal Bioanal Chem. 2012;404(5):1413–26. pmid:22814969
  48. 48. Wang Y, Armando AM, Quehenberger O, Yan C, Dennis EA. Comprehensive ultra-performance liquid chromatographic separation and mass spectrometric analysis of eicosanoid metabolites in human samples. Journal of Chromatography A. 2014;1359(0):60–9.
  49. 49. Gachet MS, Rhyn P, Bosch OG, Quednow BB, Gertsch J. A quantitiative LC-MS/MS method for the measurement of arachidonic acid, prostanoids, endocannabinoids, N-acylethanolamines and steroids in human plasma. Journal of Chromatography B. 2015;976–977(0):6–18.
  50. 50. Administration USFaD. http://www.fda.gov/downloads/Drugs/Guidances/ucm070107.pdf 2001 [cited 2013 September].
  51. 51. Yang J, Eiserich JP, Cross CE, Morrissey BM, Hammock BD. Metabolomic profiling of regulatory lipid mediators in sputum from adult cystic fibrosis patients. Free Radical Biology and Medicine. 2012;53(1):160–71. pmid:22580336
  52. 52. Mas E, Croft KD, Zahra P, Barden A, Mori TA. Resolvins D1, D2, and Other Mediators of Self-Limited Resolution of Inflammation in Human Blood following n-3 Fatty Acid Supplementation. Clinical Chemistry. 2012;58(10):1476–84. pmid:22912397
  53. 53. Zhang MJ, Spite M. Resolvins: Anti-Inflammatory and Proresolving Mediators Derived from Omega-3 Polyunsaturated Fatty Acids. Annual Review of Nutrition. 2012;32(1):203–27. pmid:22404117.
  54. 54. Simmons DL, Botting RM, Hla T. Cyclooxygenase Isozymes: The Biology of Prostaglandin Synthesis and Inhibition. Pharmacological Reviews. 2004;56(3):387–437. pmid:15317910
  55. 55. Sinclair AJ, Mann NJ. Short-term Diets Rich in Arachidonic Acid Influence Plasma Phospholipid Polyunsaturated Fatty Acid Levels and Prostacyclin and Thromboxane Production in Humans. The Journal of Nutrition. 1996;126(4 Suppl):1110S–4S. pmid:8642442
  56. 56. Agostoni C, Marangoni F, Riva E, Giovannini M, Galli C. Plasma arachidonic acid and serum thromboxane B2 concentrations in phenylketonuric children negatively correlate with dietary compliance. Prostaglandins, Leukotrienes and Essential Fatty Acids. 1997;56(3):219–22.
  57. 57. Bylesjö M, Rantalainen M, Cloarec O, Nicholson JK, Holmes E, Trygg J. OPLS discriminant analysis: combining the strengths of PLS-DA and SIMCA classification. Journal of Chemometrics. 2006;20(8–10):341–51.
  58. 58. Eriksson L, Trygg J, Wold S. CV-ANOVA for significance testing of PLS and OPLS® models. Journal of Chemometrics. 2008;22(11–12):594–600.
  59. 59. Smilowitz JT, Zivkovic AM, Wan Y- JY, Watkins SM, Nording ML, Hammock BD, et al. Nutritional lipidomics: Molecular metabolism, analytics, and diagnostics. Molecular Nutrition & Food Research. 2013;57(8):1319–35.
  60. 60. Olefsky JM, Crapo P, Reaven GM. Postprandial plasma triglyceride and cholesterol responses to a low-fat meal. The American Journal of Clinical Nutrition. 1976;29(5):535–9. pmid:178167