Molecular effects of the consumption of margarine and butter varying in trans fat composition: a parallel human intervention study
Lipids in Health and Disease volume 21, Article number: 74 (2022)
Whereas the dietary intake of industrial trans fatty acids (iTFA) has been specifically associated with inflammation, cardiovascular disease, and type 2 diabetes, understanding the impact of dietary fats on human health remains challenging owing to their complex composition and individual effects of their lipid components on metabolism. The aim of this study is to profile the composition of blood, measured by the fatty acid (FAs) profile and untargeted metabolome of serum and the transcriptome of blood cells, in order to identify molecular signatures that discriminate dietary fat intakes.
In a parallel study, the molecular effects of consuming dairy fat containing ruminant TFA (rTFA) or margarine containing iTFA were investigated. Healthy volunteers (n = 42; 45–69 y) were randomly assigned to diets containing margarine without TFA as major source of fat (wTFA control group with 0.4 g TFA per 100 g margarine), margarine with iTFA (iTFA group with 4.1 g TFA per 100 g margarine), or butter with rTFA (rTFA group with 6.3 g TFA per 100 g butter) for 4 weeks. The amounts of test products were individually selected so that fat intake contributed to 30–33% of energy requirements and TFA in the rTFA and iTFA groups contributed to up to 2% of energy intake. Changes in fasting blood values of lipid profiles (GC with flame-ionization detection), metabolome profiles (LC-MS, GC-MS), and gene expression (microarray) were measured.
Eighteen FAs, as well as 242 additional features measured by LC-MS (185) and GC-MS (54) showed significantly different responses to the diets (PFDR-adjusted < 0.05), mainly distinguishing butter from the margarine diets while gene expression was not differentially affected. The most abundant TFA in the butter, i.e. TFA containing (E)-octadec-11-enoic acid (C18:1 t11; trans vaccenic acid), and margarines, i.e. TFA containing (E)-octadec-9-enoic acid (C18:1 t9; elaidic acid) were reflected in the significantly different serum levels of TFAs measured after the dietary interventions.
The untargeted serum metabolome differentiates margarine from butter intake although the identification of the discriminating features remains a bottleneck. The targeted serum FA profile provides detailed information on specific molecules differentiating not only butter from margarine intake but also diets with different content of iTFAs in margarine.
Trans fatty acids (TFA) are fatty acids (FA) that have at least one double bond in trans configuration . Historically, the major dietary source of TFA were partially hydrogenated fats, which were widely promoted to address obesity and cardiovascular disease (CVD) . However, emergence of evidence that causally linked these fats with an increased risk of CVD [3,4,5], has led to widespread implementation of recommendations to limit consumption of partially hydrogenated oils, the majority stipulating that dietary TFA should be kept as low as possible [6, 7].
TFAs also occur naturally, in small quantities, in meat and dairy products, due to the action of bacterial enzymes that act in the rumen to bio-hydrogenate unsaturated FA. The TFA isomers from industrial (iTFA) and ruminant (rTFA) fats are chemically similar, although they vary in their relative concentrations . Given the decline in iTFA intake, the relative importance of rTFA in the diet has increased, with recent estimates showing that rTFA was the major dietary source of TFAs (0.1 to 0.7% energy intake) for more than two-thirds of the countries surveyed .
The impact of rTFA on risk of CVD is still under debate . Many epidemiological studies have demonstrated deleterious effects of iTFA but not rTFA on CVD risk [11,12,13,14], but others have reported equivalent effects of all TFA isomers on CVD risk factors [15, 16]. Furthermore, in a review of intervention studies, a positive relationship between all FAs with a double bond in trans configuration and plasma LDL to HDL cholesterol ratio was described . Some inconsistencies in these findings can be attributed to the different quantities of TFAs consumed; indeed while positive associations between circulating lipoproteins and high intakes of rTFA or iTFA have been reported, these associations were lost at lower intakes of rTFA, as generally consumed in the diet [17, 18]. Nevertheless, a systematic review of randomised clinical trials found no evidence for an association between CVD risk factors and rTFA intakes of up to 4.2% energy intake .
In contrast to the wealth of observational data associating TFA with cardiovascular health, the understanding of molecular mechanisms mediating the effects of TFA on CVD risk is more limited . However, there is evidence that TFA influences the regulation of multiple physiological processes , including hyperlipidemia , promotion of inflammation and cell death , reduced transforming growth factor-β (TGF-β) responsiveness , and endothelial dysfunction [13, 24]. Interestingly, preclinical studies comparing the major isoforms of rTFA and iTFA suggest that while these fats show some common behaviour, they can also act differentially depending on the molecular pathway assessed; for example rTFA may mediate its effects on cholesterol and FA synthesis via PPARs, while iTFA seems to specifically stimulate cholesterol synthesis via the activation of SREBP2-mediated gene regulation . Thus, evaluation of the CVD risk of rTFA may require an understanding of the molecular effects of the different fats present in the food matrix.
Importantly, although dietary TFA are not consumed in isolation but with a complex mixture of associated FAs, few human studies have taken untargeted approaches in order to evaluate how these dietary fats modify the blood composition [25, 26]. Lipid profiling together with omics technologies (transcriptomics and metabolomics) could thus be a powerful approach for demonstrating the effects of dietary fats and their TFAs on molecular and metabolic functions and pathways. An earlier cross-over study measuring the FA profile of the milk of lactating women having ingested alpine butter or margarine during 10 days demonstrated that the milk FA profile of the study participants was modified specifically by the type of ingested fat . Also, the blood FA profile was shown to associate with the adherence of young subjects to the Mediterranean diet . In addition to FA profiling, serum metabolomics, although more exploratory in nature, might also identify molecules that are sensitive to the dietary intervention . Finally, changes in the composition of the serum induced by the dietary interventions might induce a differential expression of genes in blood cells .
We previously reported the primary objectives of a parallel intervention study in which healthy subjects were allocated a diet containing up to 2% of daily energy intake from rTFA (alpine butter) or iTFA (margarine rich in iTFA) or a diet without TFA (margarine) for a duration of 4 weeks. These interventions did not show significant effects on brachial artery flow mediated dilation (primary objective) or on biomarkers of inflammation or coagulation. However, small but significant increases of total cholesterol and LDL-cholesterol were observed with the rTFA intervention (total cholesterol: rTFA group over wTFA group 1.04 (1.00 to 1.07 mmol/L); LDL-cholesterol: rTFA group over wTFA group 1.08 (1.03 to 1.14 mmol/L)) .
Although not clinically relevant these modest changes suggest that other molecular markers, which are sensitive to the dietary interventions, may be identified in the blood of the participants. The aim of the present work is therefore to make use of this study design for an exploratory analysis that investigates how the short-term intake of different sources of dietary fats differing in their composition (plant-based margarine vs animal-based butter) and content in TFA impact on the blood composition. To this end, the FA profile of blood (measured by GC-FID), the untargeted metabolome of serum (measured by GC-MS and LC-MS), and the blood cell transcriptome (measured by microarrays) were analysed and the results reported here.
Participant recruitment, inclusion and exclusion criteria, randomization and compliance have been described in detail previously . Briefly, a total of 142 healthy men and women, between 45 and 69 years were enrolled in the study. Exclusion criteria were known cardiovascular disease, smoking, hypertension, diabetes, abnormal liver and kidney function, existing signs of inflammation, anaemia, electrolyte abnormalities (Na, K, Ca), intake of medication (including supplements such as vitamins and minerals), and allergies to milk products. Of these, 125 participants were included in the per protocol analysis. Written informed consent for participation in the study was obtained from all participants and all procedures were conducted in accordance with the Declaration of Helsinki.
For the present analyses, we report on a predefined subgroup of 42 subjects selected from the main cohort using stratified sampling with respect to sex and treatment group (fourteen subjects per group). High-resolution FA analysis (GC-FID), gas chromatography mass spectrometry (GC-MS) and liquid chromatography mass spectrometry (LC-MS) analysis was completed for the entire subset while gene expression was assessed in whole blood of 21 randomly selected subjects from the same subset (seven subjects per intervention group).
Study design and interventions
The study was a randomized, controlled, double-blind, three-arm, parallel-group intervention, designed to investigate the health effects of TFA intake from industrial and ruminant sources . The study was approved by the cantonal ethical committee of Bern, Switzerland and registered at www.clinicaltrials.gov (NCT00933322). Block randomization with stratification for gender and age was used to assign subjects to one of three diets: a diet with alpine butter rich in ruminant TFA (rTFA group), a diet with margarine rich in industrial TFA (iTFA group) and a diet with margarine without TFA (wTFA group, control).
The study protocol began with a run-in period of 2 weeks, during which all study participants followed the wTFA diet, and was followed by a 4-week intervention period with the assigned diets (rTFA, iTFA, wTFA) according to the randomization . All dietary fats were provided to the subjects in the form of individual amounts of the test fats (alpine butter with 6.3 g total TFA, without conjugated linoleic acid (CLA), per 100 g; margarine with 4.1 g total TFA per 100 g; margarine without TFA containing 0.4 g total TFA per 100 g) and 15–25 g/day rapeseed oil (to balance for essential FAs) in order for the fat intake to cover 33–36% of the individual energy requirement and for TFA to contribute to max 2% of total energy intake in the rTFA and iTFA groups.
During the whole study period the participants followed a prescribed diet according to Swiss nutrition recommendations. The participants received detailed instructions about the study diet from a dietician and were followed at regular intervals during the study. Daily energy requirement was individually calculated. The participants were allowed to consume fat-free products on an individual basis according to the recommendations provided by the dietician. Throughout the study, participants were required to fill in a daily food diary. The proportion of fat was individually calculated and had to be consumed in form of the study products. Details of the dietary restrictions and evaluation of dietary adherence are published elsewhere . At each visit, weight was assessed and BMI calculated.
Blood samples were collected from the antecubital vein at baseline (after the 2-week run-in period) and at the end of the intervention period (week 6) after an overnight fast. The preparation and processing of blood samples to assess clinical biochemistry and biomarkers of inflammation, coagulation and endothelium function have been described previously . For GC-MS and LC-MS analyses, serum was separated and stored in microtubes at − 80 °C until assayed. Whole blood was collected for microarray analysis in PAXgene blood RNA tubes and frozen at − 20 °C for 24 h before transfer to − 80 °C until further treatment.
High-resolution lipid analysis with GC-FID
Serum samples were prepared for analysis by addition of 15 μL of internal standard (C13:0, 7.5 μg/15 μL) to 100 μL of serum, followed by methylation of free FAs with MeOH/HCl (25 °C for 45 min). A post-reaction treatment for neutralization was applied with Na2CO3 and extraction was performed with 300 μL hexane. Samples were measured by gas chromatography with flame ionization detector (GC-FID) (6890, Agilent Technologies, Santa Clara, CA) equipped with a CP-Sil 88 fused-silica capillary column (length 100 m, i.d. 0.25 mm, stationary-phase film thickness 0.20 μm). The carrier gas was hydrogen (1.5 ml/min). The on-column injector was hold in “Track Oven Mode”. The injector temperature was always 3 °C higher than the oven temperature. The detector temperature was 225 °C. The oven temperature was held constant for 5 min at 60 °C, then increased to 165 °C at 14 °C/min, held at 165 °C for 1 min, increased again to 225 °C at 2 °C/min, held at 225 °C for 17 min, increased again to 240 °C at 10 °C/min, and was finally held at 240 °C for 5 min. The injection volumes was 0.5 μL as described previously . Butter and margarine samples (n = 2) were assayed previously by Radtke et al. , using the same GC-FID technique as for serum samples with sample preparation according to the method of Collomb and Bühler . The data was reprocessed and reintegrated here (Table 2) with some additional reference lipids included. Eighty-five features (66 single FAs and 19 sum parameters) were quantitatively analysed in the serum and dietary fats. The FA isomers that cannot be separated by the chromatographic steps are indicated by the plus sign “+” joining various cis (c) and/or trans (t) forms of the FAs.
A control fat made from butter was used for peak assignment. In addition, a control chart with a human serum standard from Sigma-Aldrich (Human serum from human male AB plasma, H4522) was routinely run with the test samples. For each new measurement series the response factors were determined with the following fatty acid methyl esters (all from Sigma-Aldrich): Methyl butyrate (C4, C5H10O2; Methyl valerate (C5, C6H12O2); Methyl hexanoate (C6, C7H14O2); Methyl octanoate (C8, C9H18O2); Methyl nonanoate (C9, C10H20O2); Methyl decanoate (C10, C11H22O2); Methyl undecanoate (C11, C12H24O2); Methyl dodecanoate (C12, C13H26O2); Methyl tridecanoate (C13, C14H28O2); Methyl myristate (C14, C15H30O2); Methyl palmitate (C16, C17H34O2); Methyl stearate (C18, C19H38O2); Methyl oleate, (C18:1c9, C19H36O2).
Transcriptomic analysis and data pre-processing
Total RNA was extracted from whole blood samples and treated to deplete for globin mRNA according to the protocol described by Gille et al. . Whole genome transcript profiling was then performed with HG-U219 oligonucleotide expression probe arrays (Affymetrix, USA). Preparation of samples for profiling including reverse transcription, amplification, amplified RNA labelling, purification, fragmentation, and hybridization steps were performed as defined previously . Arrays were measured in accordance with manufacturer’s recommendations using Affymetrix GeneAtlas™System. Raw array data was imported to the R environment (version 4.1.2)  and corrected for background noise, log2 transformed, normalized for inter-array variation (quantile normalization)  and summarized using the rma (Robust Multichip Average) function from affy (version 1.72.0) . In addition, probe sets were filtered to keep only the most variable probe set per gene (based on standard deviation), probe sets assigned to gene symbols (hgu219.db, version 3.2.3)  and probe sets with average expression (log2) > 5. Of the 49′386 probe sets, 6′128 probe sets were retained for further analysis.
Untargeted GC-MS analysis and data pre-processing
Sample preparation for the untargeted analysis of the serum samples by GC-MS was based on the method published by Trimigno et al.  using a GC-MS 7890B/MS5977A (Agilent Technologies, Santa Clara, CA, USA) with a CombiPAL autosampler (CTC-Analytics AG, Zwingen, Switzerland). After deconvolution, features from subjects that only appeared in less than three samples were eliminated (n = 4828 remaining features). A second manual integration was conducted on 54 selected features that showed significant differences between treatment group responses to the interventions (see Statistical Analyses). The areas were normalized with the isotopically-labelled fructose. The features demonstrating a significant treatment effect were searched in EI-Mass Spectral Library (NIST 2017, Gaithersburg, MD 20899–6410, USA), masslib-library (www.masslib.com). The molecules demonstrating sufficient potential for identification were acquired from Sigma-Aldrich (Buchs, Switzerland) and analyzed by GC-MS.
Untargeted LC-MS analysis and data pre-processing
Serum samples were thawed, treated, and analyzed according to the protocol described by Pimentel et al. , using a maXis 4G+ (Bruker, Bremen, Germany) coupled to a 3000 RS UHPLC (Thermo, Basel, Switzerland). Raw data (m/z values, retention times and intensities) were imported into Progenesis QI (V2.4) (Waters, Switzerland) for retention time alignment, peak detection, sample loading normalization, deconvolution of masses (including detection of adducts), and final export of peak volumes for further statistical treatment. The pre-processed dataset (~ 85′000 features) was imported into R (version 3.5.2)  and filtered in two steps. All features that showed a mean in QC samples (n = 21) that were < 4 times higher than the mean for the final three blanks were excluded. In addition, features with a relative standard deviation > 30% in the QC samples were excluded to leave a reduced dataset of 11′616 features. Identification was carried out by a mass search in the Human Metabolome Database (HMDB v4.0 2018 ) and NIST (National Institute of Standards and Technology Database, V14). Successive identification was then confirmed by LC-MS where exact mass and retention time of the metabolites were compared with standards. The identification was differentiated at four levels [41,42,43]. Metabolites identified by comparison with a standard based on their mass (±10 ppm) and retention time (± 10%), e.g. two independent and orthogonal data, are classified as level 1. Metabolites showing at least 80% match to a theoretical MSMS spectral library but without available standard are assigned level 2. Level 3 is attributed to a compound class (e.g. if multiple similar standards match the retention time and the mass spectrum). Finally, level 4 is attributed to unidentified compounds that are nonetheless well separated from other molecules and those peak surface can be integrated.
All statistical analyses were performed in the R environment (version 3.5.2) . The normality of baseline clinical characteristics and metabolomics datasets was explored through Shapiro-Wilk tests of normality and visually by quantile-quantile plots. Principal components analysis (PCA) was used with all datasets to visually inspect sources of variation, in particular to exclude the presence of batch effects, outliers, or potential confounders (e.g. sex) (lipids and metabolomics datasets: ropls package, v1.4.4, 2016, ; transcriptomics dataset: stat function ‘prcomp’ visualized with ggplot2 (version 3.3.5) .
Non-parametric univariate statistics were used for the current analysis in view of the reduced study population size and multiple features that were not distributed normally. Baseline clinical characteristics, blood biochemistry, inflammation and endothelium biomarkers of the study population were compared by treatment group using the Kruskal-Wallis test (P < 0.05). Pairwise comparison of results were performed with Conover-Iman test (P < 0.05).
GC-FID lipid analysis
The Kruskal Wallis test was used to compare the changes in lipid profiles in response to the three interventions. For significant results, a post hoc Wilcoxon signed-rank test was completed to compare each with adjustment of P-values for multiple testing using the Benjamini-Hochberg FDR method  (significance where PFDR-adjusted < 0.05). The sum parameters for the classes of FAs were integrated into this analysis. The response to the intervention for these analyses was defined for each participant as the delta change between measures at endpoint (week 6) and baseline (week 2). Features showing a significant difference between the three treatments were presented in heatmaps using the “Metabolomics” package (v0.1.4)  with distance calculated by “Canberra” method and clustering using the “ward. D2” method .
Gene and pathway analysis
The molecular effects of the dietary interventions were investigated by assessing the change in gene expression in whole blood (i.e. delta change between baseline and post-treatment), with delta values calculated using log2-transformed data, thus interpretable as log2-fold expression changes (log2FC). Differentially expressed genes for each group and between groups were identified by performing a paired moderated t-test on the responses to the intervention using R package Limma (Linear Models for Microarray Data) (version 3.50.0) (PFDR-adjusted < 0.05) .
Gene set enrichment analysis (GSEA)  was used to explore the results of the differential analyses. GSEA determines whether a predefined set of genes show statistically significant, consistent differences between two biological states by calculation of an enrichment score (ES) for each geneset. A ‘normalized’ ES (NES) is calculated to account for differences in number of genes in each geneset. GSEA was carried out for a selection of genesets (n = 20) from the Hallmark human reference geneset collection (mSigDB, Broad Institute, version 7.5.1) that describe pathways that are regulated by TFA, including lipid and cholesterol metabolism, inflammation, apoptosis, autophagy, coagulation, adipose tissue regulation, and oxidative stress . Significance of enrichments were confirmed by comparing the ES to those obtained by random permutations (n = 1000 iterations) of the ranked gene with FDR correction for multiple testing (significance where PFDR-adjusted < 0.05) .
Untargeted metabolomics analysis
The same statistical workflow as described above for the GC-FID FA analysis was applied for the untargeted LC-MS metabolomics dataset to identify molecular biomarkers that were differentially modulated by the interventions, (comparing delta change between baseline and post-treatment). For the untargeted GC-MS dataset, a two-step statistical analysis was applied to select putative compounds that could discriminate the different effects of the interventions (delta change between baseline and post-treatment) for the reintegration step of data processing. In the first step, univariate analyses were applied to identify metabolites showing a different response between the dietary treatments using the Kruskal-Wallis test (P < 0.05). This was complemented with the identification of further discriminating metabolites from orthogonal partial least squares - discriminant analysis (OPLS-DA) models using the Ropls package (v1.4.4, 2016)  to compare the effects of each pair of treatments (for valid models, the top 100 metabolites were selected for reintegration). After reintegration, the second step of analysis applied a Wilcoxon signed-rank test to compare the responses of the reintegrated metabolites to the treatment in pairwise assessments (PFDR-adjusted < 0.05).
Correlation analyses were conducted to assess the relationships between biomarkers measured previously , including blood biochemistry (glycemia, insulinemia, blood lipids) and biomarkers of inflammation and endothelium function and the FAs measured by targeted GC-FID analysis. Spearman’s correlation test was used to associate the biomarker and FA changes in response to the dietary interventions (delta change between baseline and post-treatment) with visualization by corrplot (v0.84)  with significance where PFDR-adjusted < 0.05 . This analysis was restricted to the fourteen clinical parameters that showed no significant difference between groups at baseline. The effect of the different dietary intervention on these associations was considered by repeating the analysis with Spearman’s partial correlation test to control for the effect of diet (ppcor Package, v 1.1) .
Subject characteristics and clinical chemistry
The 42 participants selected for the nutrigenomic analyses showed similar baseline clinical characteristics to the main cohort  (Table 1). There were significant differences between the treatment groups for baseline levels of insulin, interleukin-6 (IL-6), tumor necrosis factor (TNF) receptors (1 and 2), and endothelial leucocyte adhesion molecule (ELAM-1). Therefore, these markers were not considered further.
Altered lipid profiles in serum reflect dietary fat composition
The total fat content of the three test products were very similar but differences were present for specific FAs (Table 2). The butter contained comparable, although higher, levels of total TFAs to the margarine with iTFA, and it was confirmed that the margarine without TFA (wTFA) contained only traces of TFA. However, the relative composition of the most abundant TFAs differed between the iTFA margarine and the butter; the TFA in the iTFA margarine contained both (E)-octadec-12-enoic acid (C18:1 t12) (0.3%) and the mixture of (E)-octadec-10-enoic acid and (E)-octadec-11-enoic acid (C18:1 t10 + t11) (1.4%), whereas the TFA in rTFA in butter was mainly C18:1 t10 + t11 (3.5%). In addition, the butter contained higher levels of SFA, in particular hexadecanoic acid (C16:0) and octadecanoic acid (C18:0), than the margarines despite lower levels of dodecanoic acid (C12:0). Total MUFA and total PUFAs were present at higher levels in both margarines than the butter due to higher levels of (Z)-octadec-9-enoic acid (C18:1c9) and (9Z,12Z)-octadeca-9,12-dienoic acid (C18:2c9c12), respectively.
Out of 85 features (i.e. 66 FAs and 19 sum parameters), a total of 18 FAs (four long-chain SFAs, three branched-chain FAs, three PUFAs and all TFAs) as well as six sum parameters showed significantly different responses for at least one intervention (PFDR-adjusted < 0.05) (Table 3). A heatmap of the 24 FAs (inclusive of six sum parameters) that showed significantly different responses to the treatment is presented in Fig. 1, showing clustering of the 42 subjects into the three treatment groups (iTFA, wTFA and rTFA) with six exceptions.
C18:1 t10 + t11 and C18:1 t12 were the major TFAs found in the test products. The quantified levels of these FAs in the test products (margarines and butter) are shown in Table 2, together with their relative changes in serum after the dietary interventions (Table 3). Whereas serum levels of C18:1 t10 + t11 and C18:1 t12 were both increased in the iTFA group (C18:1 t12 being significantly higher than C18:1 t10 + t11), only C18:1 t10 + t11 was increased in the rTFA group. Neither TFA was significantly modified in the wTFA group.
Total cholesterol was significantly (P < 0.05) increased in the rTFA group (0.44 mmol/L; Q1: -0.24, Q3; 0.48) compared to the wTFA group (0.02 mmol/L; Q1: 0.04, Q3; 0.63). The results of the iTFA group (0.17 mmol/L; Q1: -0.33, Q3; 0.20) were not significantly (P > 0.05) different from the two other groups. LDL cholesterol was significantly (P < 0.05) increased in the rTFA group (0.37 mmol/L; Q1: -0.10, Q3; 0.35) compared to the wTFA group (− 0.03 mmol/L; Q1: 0.18, Q3; 0.53) and the iTFA group (0.01 mmol/L; Q1: -0.24, Q3; 0.25). HDL cholesterol was not significantly (P > 0.05) increased in the rTFA group (− 0.01 mmol/L; Q1: -0.16, Q3; 0.10) compared to the wTFA group (− 0.05 mmol/L; Q1: -0.18, Q3; 0.04) and the iTFA group (0.01 mmol/L; Q1: -0.24, Q3; 0.25). Triglycerides were not significantly (P > 0.05) increased in the rTFA group (0.04 mmol/L; Q1: -0.04, Q3; 0.22) compared to the wTFA group (0.04 mmol/L; Q1: -0.08, Q3; 0.14) and the iTFA group (0.06 mmol/L; Q1: -0.11, Q3; 0.18).
Associations between circulating lipids and clinical markers
An overview of the significant correlations found between changes in plasma free FAs and clinical biomarkers after the dietary interventions is presented in Fig. 2. The strongest correlations were observed between FAs. Most of the significant associations were between lipids or between the clinical biomarkers as these two categories of parameters were ranked separately with the exception of C12:0. Also, the FAs had stronger correlations coefficients than the clinical biomarkers. Finally, most of the correlations were positive, in particular for the FAs. Nevertheless, some significant associations between targeted FAs and clinical biomarkers were identified, as detailed in Supplement Table 1. Multiple positive correlations were observed for the FAs and LDL-cholesterol including dairy-associated biomarkers, C15:0 and C17:0, several PUFAs, and individual TFAs, namely C18:1 t10 + t11 and (9E,12Z)-octadeca-9,12-dienoic acid (C18:2 t9 + c12)). All of these correlations remained significant after adjusting for the treatment effect. Several FAs that were associated with LDL-cholesterol were also positively associated with changes in total cholesterol. Conversely, HDL-cholesterol was not significantly associated with any FAs. Some FAs were positively associated with blood glycemia, hsCRP, intercellular adhesion molecule, and lipoprotein a, but these associations were generally not significant after controlling for the effect of treatment with the exception of C17:0 that was inversely associated with inter cellular adhesion molecule ICAM. No significant associations were observed between the selected lipids and triglycerides, apolipoprotein B, endothelin, TNF-α, and vascular cellular adhesion molecule (VCAM). In addition, no significant associations with clinical biomarkers were observed for the dominant TFA species in iTFA, C18:1 t6 + t7 + t8 + t9.
Limited effects of dietary fats on gene expression in blood
Overall, 496 genes were differentially expressed between baseline (week 2) and endpoint (week 6) in at least one of the three intervention groups (P < 0.05) but the responses were not statistically significant nor did they differ significantly between the intervention groups (PFDR-adjusted > 0.05). GSEA pathway enrichment did not reveal any pathway of the 20 genesets tested that was significantly differentially enriched whether for individual treatments or when comparing treatments.
Biomarkers of butter revealed in untargeted metabolomics analyses of serum
LC-MS untargeted metabolomics
A total of 185 features measured in blood by LC-MS showed significantly different responses to the three interventions (PFDR-adjusted < 0.05). These metabolites are presented in a heatmap (Fig. 3) which shows almost exclusively positive values, i.e. increased concentrations at the end of the 4-week intervention, as well as a separation of the butter group from the margarine groups while the margarine groups (wTFA and iTFA) were not differentiated. Of the 185 significant features, two were identified at level 1 (Fig. 4, Supplement Table S2). Of the identified metabolites, three showed a relative decrease in the rTFA group, relative to the wTFA and/or iTFA groups, including the fat-soluble vitamin, retinol. An inspection of the levels of these metabolites at the end of the run-in phase showed they were present at similar levels between the three treatment groups while after the butter treatment the levels of retinol and prostaglandins decreased (Supplement Fig. S1). 4-isopropylbenzoic acid showed an increase after the rTFA intervention relative to the wTFA and iTFA groups.
GC-MS untargeted metabolomics
The GC-MS analysis revealed three identified metabolites that showed significantly different responses to the treatments after manual integration of the signals (PFDR-adjusted < 0.05) (Fig. 5). Phytanic acid was clearly increased in butter treatment (P < 0.001) compared to wTFA and iTFA. An isomer of octadecadienonic acid (C18:2), was also significantly increased in the rTFA group, though it was present both in butter and margarine (P = 0.007). This molecule was classified at the level 3 of identification because of its unknown isomeric structure. Finally, glycolic acid, an α-hydroxy acid, was significantly decreased in rTFA (P = 0.037), relative to the wTFA group.
The present study investigated the molecular and metabolic impact of 4 weeks butter compared to margarine intake (with and without TFA) on fasting blood using a targeted analysis of free FA (GC-FID), transcriptomic analysis of whole blood gene expression and two untargeted metabolomic approaches (LC-MS and GC-MS).
Lipid composition of dietary interventions is reflected in serum of healthy adults
Reanalysis of the data on serum lipids (total cholesterol, LDL-C, HDL-C, TG) in the subcohort of 42 participants selected for this study showed similar effects of the dietary fat interventions than already reported for the original whole cohort of 129 participants . In particular, a small but significant increase in total cholesterol and LDL-cholesterol was observed in the group with the butter diet compared to the groups with the margarine diets. This reanalysis indicates that the sub-cohort of 42 participants is representative of the larger cohort with regard to lipid metabolism. The profiling of the lipids present in the intervention diets in serum before and after dietary interventions revealed a distinct FA signature of the intervention diets in serum. For those FAs that responded differently to the interventions, remarkable consensus could be observed between the relative change in serum and the different quantity of lipids in the products. Indeed, with the exception of total TFA that included CLA, the relative changes in serum lipids directly reflected the differences in the product composition, even for those lipids that were present in relatively low levels such as the PUFAs. We where indeed able to demonstrate this close relationship between composition of TFA in different dietary fats and circulating levels of TFA after exposure to diets with differing but still relatively low TFA content.
Among the SFA that discriminated between rTFA and iTFA diets, C15:0 and C17:0 were both specifically increased after rTFA, confirming their previously described role as dietary biomarkers of dairy intake [52, 53]. Several other odd-chain methylated FA were also raised specifically after dairy intake including 15-methylhexadecanoic acid (15Me-C17:0), which has previously been associated with butter intake .
The PUFAs (5Z,8Z,11Z,14Z)-icosa-5,8,11,14-tetraenoic acid (C20:4c; arachidonic acid), (5Z,8Z,11Z,14Z,17Z)-icosa-5,8,11,14,17-pentaenoic acid (C20:5c; EPA) and (7Z,10Z,13Z,16Z,19Z)-docosa-7,10,13,16,19-pentaenoic acid (C22:5c; DPA) were significantly increased after the rTFA diet, despite their presence at very low levels in butter (albeit higher than the margarines). C18:2c9c12 (linoleic acid), a dietary precursor of eicosanoids, was actually more abundant in the margarine diets although no significant differences in this lipid was observed in serum. The serum differences in eicosanoids therefore do not seem to be explained by the dietary composition of lipids but may rather reflect differences in metabolism of (9Z,12Z,15Z)-octadeca-9,12,15-trienoic acid (C18:3c9c12c15; linolenic acid) and C18:2c9c12 (linoleic acid). Indeed, there is some evidence to suggest that certain TFAs can disrupt the metabolism of essential PUFAs . A change in the metabolism of PUFAs could not be confirmed in our transcriptomic analysis of blood cells. As the synthesis of PUFAs is localised in the liver, blood cell transcription might not have captured such tissue-specific regulation, particularly in the fasting state. Another limitation in using fasting samples is that changes in FAs that are rapidly metabolised may not be visible. Indeed, the absence of diet-specific differences for some short-chain SFAs and MUFAs despite differences in their presence in the composition of the diets may be attributable to rapid metabolism of these fats.
Different TFAs in serum distinguish diets rich in rTFA compared to iTFA
Both diets containing TFA (iTFA and rTFA) induced similar increases in total TFA in serum that were significantly higher than the control diet without TFA (wTFA). Conversely, several TFAs showed high specificity to the source of TFA both in the product content and the serum after exposure to the diet. As expected, conjugated C18:2 FAs (CLA) were specifically associated with the rTFA diet due to their well-known production by rumen bacteria . Interestingly, diet-related changes in TFAs were found for FAs present in relatively low quantities in the diet such as C16:1 t and C18:2 t isomers. The most abundant TFAs in the rTFA and iTFA diets, respectively C18:1 t6 + t7 + t8 + t9 and C18:1 t10 + t11, were both specifically elevated in serum after the corresponding diet. The chromatographic steps used to analyse the FAs cannot separate these mixtures of TFAs. However, (E)-octadec-9-enoic acid (C18:1 t9; elaidic acid) and (E)-octadec-11-enoic acid (C18:1 t11; trans vaccenic acid) are the major TFA in margarine and butter, respectively , so that these molecules likely contribute the most to the total content in these mixtures. C18:1 t9 (elaidic acid) has already been identified as a putative dietary intake biomarker of margarine and hardened vegetable fats, as like other TFAs it cannot be synthesized by the body . Conversely, the major isomer of rTFA, C18:1 t11 (trans vaccenic acid) , is also present in dietary sources of iTFA, though at lower levels, as reflected by lower but significant increases in serum levels for the iTFA group.
Associations of C18:1 t9 (elaidic acid) and C18:1 t11 (trans vaccenic acid) to cholesterol metabolism
In our study we also explored the relationship between circulating lipids following the intervention diet and various biomarkers of molecular pathways that can be influenced by TFAs, including inflammation, lipid/cholesterol metabolism, and coagulation. Our correlation data confirmed a positive association between various TFAs that were specifically elevated after the rTFA intervention, in particular C18:1 t10 + t11, including C18:1 t11 (trans vaccenic acid), and the small but significant increase in LDL-cholesterol after rTFA, which was described previously for the full dataset . It is however important to note the similar associations observed for other FAs that were specifically elevated after rTFA including the dairy biomarkers C15:0 and C17:0. Thus the correlations might be interpreted as a collective effect of the dairy matrix on the LDL-cholesterol rather than a specific TFA-modulation of cholesterol . No associations were found between C18:1 t6 + t7 + t8 + t9, including C18:1 t9 (elaidic acid) and the circulating biomarkers of inflammation and lipid metabolism. While serum C18:1 t9 (elaidic acid) has been associated with CVD outcomes Oshita et al. , the iTFA contained relatively low levels of TFAs and serum levels of the FAs in this short-term intervention remained consistently low, which could explain the lack of effect of the intervention on surrogate markers of cholesterol metabolism.
Transcriptomics and metabolomics approaches are methods sensitive to detect early biological responses and were selected for this study to help assess the effects of the relatively low levels of TFAs in the intervention diets. However, our transcriptomics data showed limited associations with the dietary interventions, which concurred with the previously reported absence of effects of the high intake of iTFA and rTFA for 3 weeks on the plasma proteome of healthy men . Although blood transcriptomics can reflect systemic transcription changes in response to dietary interventions with fats, in particular under postprandial conditions , in this healthy population and the relatively low level of TFAs used in the diets, it was unfortunately not possible to confirm a distinct effect of the diets on the fasting transcriptome. Furthermore, only one type of metabolite was shown to distinguish between the three diets, a prostaglandin (A2/B2/J2), which was decreased after both TFA-containing diets but to a greater extent after rTFA, while prostaglandin D3 was only decreased after butter. Prostaglandins are regulatory compounds derived from dietary FAs, in particular arachidonic acid, that play important roles in many physiologic processes on many human organ systems, including inflammation . Interestingly C18:1 t9 (elaidic acid) and C18:1 t11 (trans vaccenic acid) differently impact prostaglandins production in endothelial cells . Although the prostaglandins tentatively identified in this report have hardly been reported in humans, our results suggest that plasma prostaglandins might be sensitive markers to assess the impact of dietary fats on inflammation in humans. However in the absence of responses in other metabolites or inflammatory biomarkers, this result should be interpreted with caution.
In previous studies, TFA intake has been associated with the inflammatory markers CRP and IL6, particularly in populations predisposed to metabolic illness [66, 67]. In our study with healthy subjects, no correlation between C18:1 t6 + t7 + t8 + t9, including C18:1 t9 (elaidic acid) or C18:1 t10 + t11, including C18:1 t11 (trans vaccenic acid), and the two inflammation markers were found. Similarly, Motard-Belanger et al. found no increase of CRP at the end of 4 weeks of rTFA or iTFA in a nutritional intervention, compared to the control . It is noteworthy that both diets comprise relatively low but comparable levels of TFA (2% total energy intake), a little over the WHO public health guidelines which could explain the limited clinical effects of the diet previously reported by Radtke et al. .
Biomarkers of dietary intake
The majority of metabolites that responded differently to the dietary interventions differentiated the butter effect from the similar effects of the two margarines. Several of the identified metabolites were increased in the rTFA group and might be considered as biomarkers of ruminant animal fat intake. These included 4-isopropyl-benzoic acid and FAs, octadecadienoic acid and phytanic acid. The branched-chain fatty acid (BCFA) phytanic acid is a degradation product from chlorophyll that is found in ruminant animal fat and has been proposed a putative biomarker of dairy fat [52, 68]. Some health benefits from phytanic acid intake have been described , such as prevention of metabolic syndrome or type 2 diabetes although dietary reduction of phytanic acid was also recommended for the management of infantile Refsum disease, one of the less severe of Zellweger spectrum disorders . 4-isopropyl-benzoic acid (cuminic acid) is present in the seeds of Cuminum cyminum (cumin) (https://foodb.ca/compounds/FDB013929). Although its presence in foods, in particular dairy products, has not been reported it is interesting to note that Cuminum cyminum is used as a feed additive in livestock. The untargeted GC-MS methodology used to characterize octadecadienoic acid is not specific enough to differentiate the various C18:2 molecules that can possibly be present. The reader is thus referred to the GC-FID results for further details.
Dietary biomarkers of margarine could also be inferred by considering metabolites showing a relative reduction in the rTFA group relative to the baseline, which was the end of the run-in with the wTFA margarine. These included retinol, which although present in butter is supplemented, often at high levels, in margarines. Of note, dihydrophylloquinone, a hydrogenated form of vitamin K1 not naturally present in vegetable oils, has previously been reported as a candidate biomarker of TFA-intake from partially hydrogenated fat  but was not found with our untargeted LC-MS method.
For most of the molecules described above, the changes in their serum levels during the treatment phase are most likely attributed to the test products rather than to other factors such as endogenous metabolism or the intake of other food components independent of the intervention. Excluding these alternative explanations for the observed associations of the identified metabolites with the dietary intervention could be more precisely addressed by a postprandial analysis of serum in controlled acute intervention studies using butter and margarine as nutritional challenges.
Strengths and limitations of study
A strength of this study was the use of different techniques to investigate the molecular effects of the dietary interventions under study. In particular, by using a targeted FA panel, we were able to discriminate between the circulating changes in specific FAs and TFAs that were present in the dietary sources of TFAs. This characteristics is important in the context of the distinct biological roles of these lipids. Using different metabolomics approaches in this study, we were also able to consider the wider impact of the diets on metabolites other than lipids.
Although the clustering of LC-MS data was more efficient than the targeted FA analytics to separate the intake of butter from the intake of margarine, we were limited by the challenge of identifying discriminating metabolites. A range of issues associated with the identification of the molecular features derived from metabolomics studies contributes to this bottleneck, including the quality of the databases, the stability of the laboratory metabolomics instrumentation (columns, detector …), and the availability of standards. Each of these parameters would need to be improved to increase the success rate of the identification step. The original study defined its number of participants based on the primary endpoint, brachial artery flow mediated dilation. In light of the extensive work load needed to conduct FA analysis the number of subjects was limited to 42 and this study was thus conducted on an exploratory basis. This number of subjects was, indeed, sufficient to identify differentiating metabolites, as confirmed with a recent studies investigating FA blood profiles in a similar number of subjects as a function of dietary adherence to the Mediterranean diet . Our study was further limited by the restricted number of samples evaluated for whole blood transcriptomics, which did not capture the expected changes in lipid metabolism following the change in dietary fats. Although a previous study has reported changes in the fasting blood cell transcriptome of seven subjects as a result of a 2-week dairy intervention , our gene expression study remain exploratory: an higher number of subjects might have provided stronger insights into the reasons for the lack of changes in gene expression (i.e. insufficient statistical power vs real biological lack of response). However, gene expression remains a sensitive method to identify subtle metabolic change resulting from dietary interventions. Given that none of the 496 differentially expressed genes identified in this study retained statistical significance after correction for multiple testing, a targeted strategy specifically focusing on a small set of pre-defined genes with relevant function might have been more appropriate.
The application of a combination of targeted and untargeted metabolomics to serum samples of humans supported the identification of distinct chemical signatures associated with the different dietary fats. We confirm that the different types and abundance of TFA isomers in rTFA and iTFA diets results in distinct changes in serum TFAs after 4 weeks of these diets. In particular, C18:1 t6 + t7 + t8 + t9, including C18:1 t9 (elaidic acid) and C18:1 t10 + t11, including C18:1 t11 (trans vaccenic acid) are differentially modulated by the TFA diets in concordance with their relative presence in the dietary fats. The associations of the different TFAs with clinical biomarkers did not seem specific to the TFAs but was rather related to the overall lipid composition of the dietary fats. Beyond the lipid analysis, we did not find clear biomarkers that distinguished the rTFA from iTFA perhaps reflecting the relatively low dose of TFAs used that did not result in clinically different phenotypes in the larger cohort of this study . However, we did identify biomarkers that distinguished butter from the margarines including some putative exogenous markers, suggesting that this approach is useful for discriminating metabolic effects of similar foods.
Availability of data and materials
The transcriptome and metabolome datasets generated and analysed during the current study are available from the corresponding author on reasonable request.
Conjugated linoleic acid
GC with flame ionization detector
Geneset enrichment analyses
High-density lipoprotein cholesterol
High-sensitivity C-reactive protein
Intercellular adhesion molecule
Industrial trans fatty acids
Low-density lipoprotein cholesterol
Monounsaturated fatty acids
Orthogonal partial least squares – discriminant analysis
Polyunsaturated fatty acids
Ruminant trans fatty acids
Saturated fatty acids
Trans fatty acids
Tumor necrosis factor
Vascular cellular adhesion molecule
Without trans fatty acids
Lichtenstein AH. Fatty Acids: Trans Fatty Acids. In: Encyclopedia of Food and Health: Elsevier Ltd; 2016. p. 645–8.
Bloom DE, Cafiero ET, Jané-Llopis E, Abrahams-Gessel S, Bloom LR, Fathima S, et al. The global economic burden of non-communicable diseases. World Economic Forum: Geneva; 2011.
Willett WC, Stampfer MJ, Manson JE, Colditz GA, Speizer FE, Rosner BA, et al. Intake of trans fatty acids and risk of coronary heart disease among women. Lancet. 1993;341:581–5.
European Food Safety Authority. Scientific opinion on dietary reference values for fats, including saturated fatty acids, polyunsaturated fatty acids, monounsaturated fatty acids, trans fatty acids, and cholesterol. EFSA J. 2010;8:107.
Mozaffarian D. Trans fatty acids - effects on systemic inflammation and endothelial function. Atheroscler Suppl. 2006;7:29–32.
European Food Safety Authority (EFSA). Scientific and technical assistance on trans fatty acids. EFSA J. 2018;15:EN-1433.
World Health Organization. Countdown to 2023: WHO report on global trans fat elimination 2021. Geneva: World Health Organization; 2021. Licence: CC BY-NC-SA 3.0 IGO.
Brouwer IA, Wanders AJ, Katan MB. Effect of animal and industrial trans fatty acids on HDL and LDL cholesterol levels in humans--a quantitative review. Plos One. 2010;5:e9434.
Wanders AJ, Zock PL, Brouwer IA. Trans fat intake and its dietary sources in general populations worldwide: a systematic review. Nutrients. 2017;9:840.
Sébédio JL, Christie WW. Metabolism of trans fatty acid isomers. In: Trans Fatty Acids in Human Nutrition. Destaillats F, Sébédio JL, Dionisi F, Chardigny JM, editors. 2nd Edition. Bridgwater: The Oily Press; 2009.
de Souza RJ, Mente A, Maroleanu A, Cozma AI, Ha V, Kishibe T, et al. Intake of saturated and trans unsaturated fatty acids and risk of all cause mortality, cardiovascular disease, and type 2 diabetes: systematic review and meta-analysis of observational studies. BMJ. 2015;351:h3978.
Jakobsen MU, Overvad K, Dyerberg J, Heitmann BL. Intake of ruminant trans fatty acids and risk of coronary heart disease. Int J Epidemiol. 2008;37:173–82.
Mozaffarian D, Aro A, Willett WC. Health effects of trans-fatty acids: experimental and observational evidence. Eur J Clin Nutr. 2009;63(Suppl 2):S5–21.
Bendsen NT, Christensen R, Bartels EM, Astrup A. Consumption of industrial and ruminant trans fatty acids and risk of coronary heart disease: a systematic review and meta-analysis of cohort studies. Eur J Clin Nutr. 2011;65:773–83.
Oomen CM, Ocke MC, Feskens EJ, van Erp-Baart MA, Kok FJ, Kromhout D. Association between trans fatty acid intake and 10-year risk of coronary heart disease in the Zutphen elderly study: a prospective population-based study. Lancet. 2001;357:746–51.
Motard-Belanger A, Charest A, Grenier G, Paquin P, Chouinard Y, Lemieux S, et al. Study of the effect of trans fatty acids from ruminants on blood lipids and other risk factors for cardiovascular disease. Am J Clin Nutr. 2008;87:593–9.
Gebauer SK, Destaillats F, Dionisi F, Krauss RM, Baer DJ. Vaccenic acid and trans fatty acid isomers from partially hydrogenated oil both adversely affect LDL cholesterol: a double-blind, randomized controlled trial. Am J Clin Nutr. 2015;102:1339–46.
Chardigny JM, Destaillats F, Malpuech-Brugere C, Moulin J, Bauman DE, Lock AL, et al. Do trans fatty acids from industrially produced sources and from natural sources have the same effect on cardiovascular disease risk factors in healthy subjects? Results of the trans fatty acids collaboration (TRANSFACT) study. Am J Clin Nutr. 2008;87:558–66.
Gayet-Boyer C, Tenenhaus-Aziza F, Prunet C, Marmonier C, Malpuech-Brugere C, Lamarche B, et al. Is there a linear relationship between the dose of ruminant trans-fatty acids and cardiovascular risk markers in healthy subjects: results from a systematic review and meta-regression of randomised clinical trials. Br J Nutr. 2014;112:1914–22.
Oteng AB, Kersten S. Mechanisms of action of trans fatty acids. Adv Nutr. 2020;11:697–708.
Lin J, Yang R, Tarr PT, Wu PH, Handschin C, Li S, et al. Hyperlipidemic effects of dietary saturated fats mediated through PGC-1beta coactivation of SREBP. Cell. 2005;120:261–73.
Hirata Y, Takahashi M, Kudoh Y, Kano K, Kawana H, Makide K, et al. Trans-fatty acids promote proinflammatory signaling and cell death by stimulating the apoptosis signal-regulating kinase 1 (ASK1)-p38 pathway. J Biol Chem. 2017;292:8174–85.
Chen CL, Tetri LH, Neuschwander-Tetri BA, Huang SS, Huang JS. A mechanism by which dietary trans fats cause atherosclerosis. J Nutr Biochem. 2011;22:649–55.
Lopez-Garcia E, Schulze MB, Meigs JB, Manson JE, Rifai N, Stampfer MJ, et al. Consumption of trans fatty acids is related to plasma biomarkers of inflammation and endothelial dysfunction. J Nutr. 2005;135:562–6.
Lichtenstein AH, Erkkila AT, Lamarche B, Schwab US, Jalbert SM, Ausman LM. Influence of hydrogenated fat and butter on CVD risk factors: remnant-like particles, glucose and insulin, blood pressure and C-reactive protein. Atherosclerosis. 2003;171:97–107.
Shahin AM, McGuire MK, Anderson N, Williams J, McGuire MA. Effects of margarine and butter consumption on distribution of trans-18:1 fatty acid isomers and conjugated linoleic acid in major serum lipid classes in lactating women. Lipids. 2006;41:141–7.
Bertschi I, Collomb M, Rist L, Eberhard P, Sieber R, Bütikofer U, et al. Maternal dietary Alpine butter intake affects human milk: fatty acids and conjugated linoleic acid isomers. Lipids. 2005;40:581–7.
Syrén ML, Turolo S, de Marco EA, De Cosmi V, Risé P, Marangoni F, et al. Whole blood fatty acid profile of young subjects and adherence to the Mediterranean diet: an observational cohort study. Lipids Health Dis. 2022;21:23.
Pimentel G, Burton KJ, von Ah U, Butikofer U, Pralong FP, Vionnet N, et al. Metabolic Footprinting of fermented Milk consumption in serum of healthy men. J Nutr. 2018;148:851–60.
Burton KJ, Pimentel G, Zangger N, Vionnet N, Drai J, McTernan PG, et al. Modulation of the peripheral blood transcriptome by the ingestion of probiotic yoghurt and acidified milk in healthy, young men. Plos One. 2018;13:e0192947.
Radtke T, Schmid A, Trepp A, Dahler F, Coslovsky M, Eser P, et al. Short-term effects of trans fatty acids from ruminant and industrial sources on surrogate markers of cardiovascular risk in healthy men and women: a randomized, controlled, double-blind trial. Eur J Prev Cardiol. 2017;24:534–43.
Collomb M, Bühler T. Analyse de la composition en acides gras de la graisse de lait. Mitteilungen aus Lebensmitteluntersuchung und Hygiene. 2000;91:306–32.
Gille D, Zangger N, Soneson C, Butikofer U, Delorenzi M, Schwander F, et al. Caloric dose-responsive genes in blood cells differentiate the metabolic status of obese men. J Nutr Biochem. 2017;43:156–65.
R Core Team. R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing, Vienna. 2019. https://www.R-project.org.
Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, Scherf U, et al. Exploration, normalization, and summaries of high density oligonucleotide Array probe level data. Biostatistics. 2003;4:249–64.
Gautier L, Cope L, Bolstad BM, Irizarry RA. affy-analysis of Affymetrix GeneChip data at the probe level. Bioinformatics. 2004;20:307–15.
hgu219.db: Affymetrix Human Genome 219 Plate annotation data (chip hgu219). R package version 3.2.3. https://bioconductor.org/packages/release/data/annotation/html/hgu219.db.html.
Trimigno A, Munger L, Picone G, Freiburghaus C, Pimentel G, Vionnet N, et al. GC-MS based metabolomics and NMR spectroscopy investigation of food intake biomarkers for Milk and cheese in serum of healthy humans. Metabolites. 2018;8:26.
R Core Team. R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing, Vienna. 2018. https://www.R-project.org.
Wishart DS, Feunang YD, Marcu A, Guo AC, Liang K, Vazquez-Fresno R, et al. HMDB 4.0: the human metabolome database for 2018. Nucleic Acids Res. 2018;46:D608–17.
Sumner LW, Amberg A, Barrett D, Beale MH, Beger R, Daykin CA, et al. Proposed minimum reporting standards for chemical analysis chemical analysis working group (CAWG) metabolomics standards initiative (MSI). Metabolomics. 2007;3:211–21.
Ulaszewska MM, Weinert CH, Trimigno A, Portmann R, Andres Lacueva C, Badertscher R, et al. Nutrimetabolomics: an integrative action for Metabolomic analyses in human nutritional studies. Mol Nutr Food Res. 2019;63:e1800384.
Chaleckis R, Meister I, Zhang P, Wheelock CE. Challenges, progress and promises of metabolite annotation for LC-MS-based metabolomics. Curr Opin Biotechnol. 2019;55:44–50.
Thevenot EA, Roux A, Xu Y, Ezan E, Junot C. Analysis of the human adult urinary metabolome variations with age, body mass index, and gender by implementing a comprehensive workflow for univariate and OPLS statistical analyses. J Proteome Res. 2015;14:3322–35.
Wickham H. ggplot2 elegant graphics for data analysis. New York: Springer-Verlag; 2016.
Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc B. 1995;57:289–300.
De Livera AM, Bowne JB. metabolomics: Analysis of Metabolomics Data. R package version 0.1.4; 2014.
Murtagh F, Legendre P. Ward's hierarchical agglomerative clustering method: which algorithms implement Ward‘s criterion? J Classif. 2014;31:274–95.
Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 2005;102:15545–50.
Wei T, Simko V. R package "corrplot": Visualization of a Correlation Matrix (Version 0.84); 2017.
Kim S. Ppcor: an R package for a fast calculation to semi-partial correlation coefficients. Commun Stat Appl Methods. 2015;22:665–74.
Munger LH, Garcia-Aloy M, Vazquez-Fresno R, Gille D, Rosana ARR, Passerini A, et al. Biomarker of food intake for assessing the consumption of dairy and egg products. Genes Nutr. 2018;13:26.
Albani V, Celis-Morales C, Marsaux CF, Forster H, O'Donovan CB, Woolhead C, et al. Exploring the association of dairy product intake with the fatty acids C15:0 and C17:0 measured from dried blood spots in a multipopulation cohort: findings from the Food4Me study. Mol Nutr Food Res. 2016;60:834–45.
Guertin KA, Moore SC, Sampson JN, Huang WY, Xiao Q, Stolzenberg-Solomon RZ, et al. Metabolomics in nutritional epidemiology: identifying metabolites associated with diet and quantifying their potential to uncover diet-disease relations in populations. Am J Clin Nutr. 2014;100:208–17.
Saadi HAS. Trans-Fatty Acids. In: Encyclopedia of Toxicology, vol. 4: Elsevier Inc.; 2014. p. 798–801.
Ferlay A, Bernard L, Meynadier A, Malpuech-Brugere C. Production of trans and conjugated fatty acids in dairy ruminants and their putative effects on human health: a review. Biochimie. 2017;141:107–20.
Li C, Richter P, Cobb LK, Kuiper HC, Seymour J, Vesper HW. Dietary Sources of Plasma trans Fatty Acids among Adults in the United States: NHANES 2009–2010. Curr Dev Nutr. 2021;5:nzab063.
Saadatian-Elahi M, Slimani N, Chajes V, Jenab M, Goudable J, Biessy C, et al. Plasma phospholipid fatty acid profiles and their association with food intakes: results from a cross-sectional study within the European prospective investigation into Cancer and nutrition. Am J Clin Nutr. 2009;89:331–46.
Weggemans RM, Rudrum M, Trautwein EA. Intake of ruminant versus industrial trans fatty acids and risk of coronary heart disease – what is the evidence? Eur J Lipid Sci Technol. 2004;106:390–7.
Engel S, Tholstrup T. Butter increased total and LDL cholesterol compared with olive oil but resulted in higher HDL cholesterol compared with a habitual diet. Am J Clin Nutr. 2015;102:309–15.
Oshita T, Toh R, Shinohara M, Mori K, Irino Y, Nagao M, et al. Elevated serum Elaidic acid predicts risk of repeat revascularization after percutaneous coronary intervention in Japan. Circ J. 2019;83:1032–8.
de Roos B, Wanders AJ, Wood S, Horgan G, Rucklige G, Reid M, et al. A high intake of industrial or ruminant trans fatty acids does not affect the plasma proteome in healthy men. Proteomics. 2011;11:3928–34.
van Dijk SJ, Mensink M, Esser D, Feskens EJ, Muller M, Afman LA. Responses to high-fat challenges varying in fat type in subjects with different metabolic risk phenotypes: a randomized trial. Plos One. 2012;7:e41388.
Miller SB. Prostaglandins in health and disease: an overview. Semin Arthritis Rheum. 2006;36:37–49.
Da Silva MS, Bilodeau JF, Larose J, Greffard K, Julien P, Barbier O, et al. Modulation of the biomarkers of inflammation and oxidative stress by ruminant trans fatty acids and dairy proteins in vascular endothelial cells (HUVEC). Prostaglandins Leukot Essent Fatty Acids. 2017;126:64–71.
Mozaffarian D, Pischon T, Hankinson SE, Rifai N, Joshipura K, Willett WC, et al. Dietary intake of trans fatty acids and systemic inflammation in women. Am J Clin Nutr. 2004;79:606–12.
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. Am J Clin Nutr. 2004;79:969–73.
Allen NE, Grace PB, Ginn A, Travis RC, Roddam AW, Appleby PN, et al. Phytanic acid: measurement of plasma concentrations by gas-liquid chromatography-mass spectrometry analysis and associations with diet and other plasma fatty acids. Br J Nutr. 2008;99:653–9.
Roca-Saavedra P, Marino-Lorenzo P, Miranda JM, Porto-Arias JJ, Lamas A, Vazquez BI, et al. Phytanic acid consumption and human health, risks, benefits and future trends: a review. Food Chem. 2017;221:237–47.
Sá MJ, Rocha JC, Almeida MF, Carmona C, Martins E, Miranda V, et al. Infantile Refsum disease: influence of dietary treatment on plasma Phytanic acid levels. JIMD Rep. 2016;26:53–60.
Lichtenstein AH. Fatty acids: trans fatty acids. In: Encyclopedia of food and Health; 2016. p. 645–8.
The authors thank Monika Spahni for measuring the fatty acids in blood.
This work was funded by Agroscope.
Ethics approval and consent to participate
Ethical approval was obtained from the cantonal ethical committee of Bern, Switzerland (KEK 124/09) and the study was conducted in accordance with the Declaration of Helsinki. All participants provided written informed consent for participation in the study. The study was registered at www.clinicaltrials.gov (NCT00933322).
Consent for publication
All authors consent to publish the present results.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Individual changes in identified LC metabolites during study intervention. Individual line plots coloured by intervention groups: wTFA (green), iTFA (red), rTFA (blue). One outlier in wTFA lineplot for 4-isopropylbenzoic acid was eliminated (baseline: 475.6, endpoint: 6464.7). This subject was however not removed from the robust statistical analysis. Supplement Table S1. Spearman correlations between fatty acids that responded differently to the interventions and clinical parameters using delta values. Significance for Padjusted FDR < 0.05. Supplement Table 2. Identification of LC-MS metabolites.
About this article
Cite this article
Guggisberg, D., Burton-Pimentel, K.J., Walther, B. et al. Molecular effects of the consumption of margarine and butter varying in trans fat composition: a parallel human intervention study. Lipids Health Dis 21, 74 (2022). https://doi.org/10.1186/s12944-022-01675-1