Skip to main content
Download PDF

Differences in peripheral blood mononuclear cell gene expression and triglyceride composition in lipoprotein subclasses in plasma triglyceride responders and non-responders to omega-3 supplementation

Genes & Nutrition volume 14, Article number: 10 (2019) Cite this article

Abstract

Background

Intake of the marine omega-3 fatty acids eicosapentaenoic acid (EPA) and docosahexaenoic acid (DHA) reduces fasting triglyceride (TG) levels and may thereby lower cardiovascular disease risk. However, there are large inter-individual differences in the TG-lowering effect of omega-3 supplementation. Genotype differences partly explain this variation, but gene-environment interactions leading to gene expression differences may also be important. In this study, we aimed to investigate baseline differences and differences in the change in peripheral blood mononuclear cell (PBMC) gene expression and lipoprotein subclass TG levels between TG responders and non-responders to omega-3 fatty acid supplementation.

Methods

In a previous randomized controlled trial, healthy normotriglyceridemic subjects (n = 35, 71% women) received 1.6 g EPA + DHA/day for 7 weeks. In this exploratory sub-study, we defined TG responders as subjects having a TG reduction beyond the 20% day-to-day variation and non-responders as having a TG change between − 20% and + 20% after omega-3 supplementation. PBMC gene expression was measured using microarray, and lipoprotein subclasses were measured using nuclear magnetic resonance spectroscopy.

Results

Eight subjects were defined as responders with a median TG reduction of 37%, and 16 subjects were defined as non-responders with a median TG change of 0%. At baseline, responders had higher TG levels in two of four high-density lipoprotein (HDL) subclasses and 909 gene transcripts (p ≤ 0.05) were differentially expressed compared to non-responders. During the intervention, the plasma TG reduction among responders was reflected in TG reductions in four of six different very low-density lipoprotein subclasses and three of four different HDL subclasses. Compared to non-responders, the expression of 454 transcripts was differentially altered in responders (p ≤ 0.05). Pathway analyses revealed that responders had altered signaling pathways related to development and immune function. In addition, two of the top 10 enriched pathways in responders compared to non-responders were related to lysophosphatidic acid signaling.

Conclusion

TG responders and non-responders to omega-3 supplementation have different lipoprotein subclass and PBMC gene expression profiles at baseline and different lipoprotein subclass and PBMC gene expression responses to omega-3 supplementation. These gene expression differences may partially explain the variability in TG response observed after omega-3 supplementation.

Graphical abstract

Based on free images from Servier Medical Art (Creative Commons Attribution License) and image from www.colourbox.com.

Background

Intake of the marine omega-3 fatty acids eicosapentaenoic acid (EPA) and docosahexaenoic acid (DHA) may reduce the risk of cardiovascular disease (CVD), especially in high-risk populations [1]. This risk reduction is linked to the anti-inflammatory, anti-arrhythmic, blood pressure-lowering, and lipid-modifying effects of omega-3 fatty acids [2,3,4,5]. The reduction in triglyceride (TG) levels is one of the most important effects of omega-3 fatty acids, and this effect is dependent on pre-supplementation TG levels and the omega-3 fatty acid dose [2, 3]. However, there are large inter-individual differences in the TG-lowering effect of omega-3 fatty acids, with several studies showing that about 30–40% of participants do not obtain reduced TG levels following omega-3 supplementation (non-responders) [6,7,8]. Those who have a TG-lowering effect of omega-3 supplementation (responders) seem to have a less healthy biochemical profile, such as higher TG and glucose levels and lower HDL-C. In addition, responders have more favorable changes in total cholesterol and HDL-C in response to omega-3 supplementation than non-responders [9]. However, the distribution of TG in lipoprotein subclasses in responders and non-responders is less characterized.

Furthermore, polymorphisms in the genes apolipoprotein E (APOE), acetyl-CoA carboxylase α (ACACA), ATP citrate lyase (ACLY), cluster of differentiation 36 (CD36), retinoid X receptor α (RXRA), and acyl-CoA oxidase 1 (ACOX1), among others, have been found to affect the TG-lowering effect of omega-3 fatty acids [10,11,12,13]. It has previously been shown in a genome-wide association study (GWAS) that genotype explains only 20% of the variation in TG response to omega-3 fatty acids in the Fatty Acid Sensor (FAS) study [14]. However, in the FAS study, a more refined and improved genetic risk score (GRS) has recently showed that GRS can explain almost 50% of the TG response variance [15].

Omega-3 fatty acids mediate their effects largely at the cellular level, such as by their ability to alter gene expression. This can happen directly, by binding and activating nuclear receptors such as peroxisome proliferator-activated receptors (PPARs), or indirectly by inhibiting the nuclear translocation of transcription factors (TFs) including nuclear factor kappa B (NF-κB) and sterol regulatory element binding protein (SREBP) [16,17,18]. These TFs affect the transcription of genes that, among others, are involved in lipid metabolism and inflammation [19].

Peripheral blood mononuclear cells (PBMC) are cells of both the innate and adaptive immune system and are mainly composed of lymphocytes and monocytes. Because they are circulating cells, they are exposed to nutrients, metabolites, and peripheral tissues and may therefore reflect whole-body health [20, 21]. It has been shown that PBMC gene expression reflects liver and adipose tissue expression of genes involved in lipid metabolism and inflammation [22,23,24,25], and omega-3 fatty acids have been shown to alter PBMC gene expression [26,27,28]. Hence, PBMCs may be a good a model system for exploring the underlying mechanisms of the TG-lowering effect of omega-3 fatty acids.

In this exploratory study, we aimed to analyze the baseline differences and the difference in change in PBMC transcriptome and the TG content in lipoprotein subclasses between TG responders and non-responders to omega-3 supplementation. Since the effect of omega-3 fatty acids on lipid metabolism is partly mediated through effects on gene expression, we hypothesized that intake of omega-3 fatty acids would differentially affect PBMC gene expression in TG responders and non-responders.

Results

Characteristics of responders and non-responders, biochemical parameters, plasma fatty acids, and estimated omega-3 index

After 7 weeks of omega-3 supplementation, 8 participants who received omega-3 supplementation were defined as responders with a median fasting TG reduction of 37% and 16 participants were defined as non-responders with a median change in fasting TG of 0% (Fig. 1a). The median compliance estimated by capsule count was 100% in responders and 99.5% in non-responders, and the difference was not significant. At baseline, the age of the participants, the distribution of men and women, the intervention group allocation, and BMI did not differ between responders and non-responders (Table 1). Nonetheless, responders had higher baseline levels of fasting TG (p = 0.01) and plasma oleic acid (OA; p = 0.007) and lower baseline levels of plasma linoleic acid (LA; p = 0.009) compared to non-responders (Table 1). The baseline levels of plasma EPA, docosapentaenoic acid (DPA), and DHA did not differ between responders and non-responders. However, the baseline estimated omega-3 index was 4.6% in responders and 5.6% in non-responders and the baseline difference was significant (Table 1).

Fig. 1
figure 1

Individual changes in fasting TG, weight, and plasma levels of omega-3 fatty acids. The relative change in fasting TG levels (a) were used to categorize participants as responders (ΔTG ≤ − 20%, green) and non-responders (− 20% < ΔTG < + 20%, red). The participants are presented in the same order for the absolute change (mmol/L) in fasting TG levels (b), weight change (c), and the change in plasma levels of marine omega-3 fatty acids (d)

Full size image
Table 1 Baseline characteristics, biochemical parameters, plasma fatty acids, and omega-3 index of responders and non-responders to omega-3 supplementation
Full size table

In whole blood, there were 19% lymphocytes at baseline in both responders and non-responders, while at the 7-week visit there was 22% and 20% lymphocytes in responders and non-responders, respectively. There were 4% monocytes in both responders and non-responders at the baseline and 7-week visits. The change of lymphocytes and monocytes during the study did not differ between responders and non-responders.

During the intervention, there was no significant weight change within or between responders and non-responders (Fig. 1c). Evidently, the change in TG levels in responders and non-responders was significantly different (p = 0.0001) (Table 1), and as shown in Fig. 1b, the TG change ranged from − 1.7 to − 0.2 mmol/L in responders and from − 0.4 to 0.1 mmol/L in non-responders. The change in plasma EPA, DPA, and DHA did not differ between responders and non-responders (Table 1), as also illustrated by the individual changes in total plasma omega-3 levels in Fig. 1d. Nonetheless, there were significant differences between responders and non-responders in the change of several plasma fatty acids. The plasma level of OA was more reduced in responders compared to non-responders (p = 0.004), LA increased in responders while it decreased in non-responders (p = 0.04), the reduction in α-linolenic acid (ALA) was greater in responders than non-responders (p = 0.03), and the reduction in arachidonic acid (AA) was greater in non-responders than responders (p = 0.01) (Table 1).

Lipoprotein subclasses

To investigate the baseline differences and the difference in change in TG levels in more detail, we analyzed TG levels in 14 different lipoprotein subclasses. At baseline, responders had higher levels of TG in the two smallest HDL subclasses, M-HDL and S-HDL (Fig. 2). However, there were no significant differences between responders and non-responders in the baseline levels of TG in any of the very low-density lipoprotein (VLDL), intermediate-density lipoprotein (IDL), or LDL subclasses (not shown). Furthermore, we found that responders had a greater reduction of TG in extra large (XL-), large (L-), medium (M-), and small (S-)VLDL subclasses (Fig. 3) as well as in the XL-, M-, and S-HDL subclasses (Fig. 4), but there were no changes in the TG levels in IDL or the LDL subclasses (not shown).

Fig. 2
figure 2

Baseline triglyceride levels in HDL subclasses. The baseline TG levels in XL-, L-, M-, and S-HDL subclasses (mmol/L) in responders (green) and non-responders (red). Differences are tested on log2-transformed data with an independent t test

Full size image
Fig. 3
figure 3

Change in triglyceride levels in VLDL subclasses. The change in TG levels in XXL-, XL-, L-, M-, S-, and XS-VLDL subclasses from baseline to end of study (Δmmol/L) in responders (green) and non-responders (red). Differences are tested on the log2 ratio with an independent t test

Full size image
Fig. 4
figure 4

Change in triglyceride levels in HDL subclasses. The change in TG levels in XL-, L-, M-, and S-HDL subclasses from baseline to the end of study (Δmmol/L) in responders (green) and non-responders (red). Differences are tested on the log2 ratio with an independent t test

Full size image

PBMC gene expression

We analyzed the change in whole genome PBMC transcriptome at baseline and after 7 weeks of omega-3 supplementation. As previously described, all arrays fulfilled the quality criteria and were included in the analyses [27]. At baseline, 909 transcripts were differentially expressed in responders compared to non-responders (p ≤ 0.05). Of these transcripts, 458 transcripts had a lower expression while 451 transcripts had a higher expression in responders compared to non-responders (Additional file 1). Only one of these transcripts, the long non-coding RNA LINC00473, was expressed at a significantly higher level in responders after adjusting for multiple comparisons (FDR < 25%). The transcripts with the largest difference in expression at baseline (log2 ratio < − 0.4 and log2 ratio > 0.4) encoded ribosomal and cytoskeleton proteins, among others (Table 2).

Table 2 Genes with the highest and lowest baseline expression (log2 ratio > |0.4|) in responders compared to non-responders to omega-3 supplementation
Full size table

After 7 weeks of omega-3 supplementation, 454 transcripts were differentially altered (p ≤ 0.05) in responders compared to non-responders. Of these, 246 transcripts were reduced and 208 were increased (Additional file 2). However, after correction for multiple testing, there were no differences between responders and non-responders in change in mRNA expression from baseline to 7 weeks (FDR < 25%). Among the significantly altered transcripts, those that had the greatest reduction (log2 ratio < − 0.25) and the highest increase (log2 ratio > 0.25) in responders compared to non-responders encode proteins involved in regulation of gene expression, apoptosis, and immune function (Table 3).

Table 3 Genes with the most altered expression (log2 ratio > |0.25|) in responders compared to non-responders to omega-3 supplementation
Full size table

Pathway and transcription factor analyses

To examine processes differentially altered in responders compared to non-responders after 7 weeks of omega-3 supplementation, we performed functional analyses for the 454 altered transcripts. Of the 73 enriched pathways (FDR < 25%, Additional file 3), many were related to signaling pathways involved in development and immune function. Two of the top 10 enriched pathways (Table 4) were related to lysophosphatidic acid (LPA) signaling. In addition, some pathways related to regulation of lipid metabolism, adipogenesis, blood coagulation, and thromboxane A2 signaling were enriched in responders compared to non-responders (Additional file 3).

Table 4 Top 10 enriched pathways in responders compared to non-responders after 7 weeks of omega-3 supplementation
Full size table

To further examine the mechanisms behind the differences in transcription levels between responders and non-responders, transcription factor (TF) analyses were performed. We found that genes with binding sites for 78 TFs were overrepresented in the list of genes that were significantly differentially altered in responders compared to non-responders including the TF hepatic nuclear factor 4 alpha (HNF4-α; Additional file 4).

Discussion

In the present exploratory study, we found that intake of omega-3 fatty acids differentially altered PBMC gene expression in TG responders and non-responders. Specifically, enriched pathways in responders compared to non-responders were related to development and apoptosis, immune function, and LPA signaling. These results lend further support to the findings of Rudkowska et al. who investigated transcriptomic and metabolomics profiles of TG responders and non-responders to omega-3 supplementation. They reported that 1020 transcripts were altered within the non-responder group, and 252 transcripts were altered within the responder group with only 10 transcripts overlapping between the groups [9]. In the current study, we also report that responders had higher baseline TG in M- and S-HDL subclasses and a greater reduction in TG levels in four of six VLDL and three of four HDL subclasses than non-responders.

Among the enriched pathways in responders compared to non-responders, several were related to development signaling and immune response. In line with this, omega-3 supplementation in elderly subjects altered immune-related pathways in PBMCs [26]. This may be expected, as PBMCs are cells of the immune system, and omega-3 fatty acids alter immune responses through altering NF-κB- and PPAR-induced gene expression and by acting as precursors for the production of anti-inflammatory and pro-resolving lipid mediators [4, 29]. In addition, pathways related to immune response and apoptosis have previously been found to be altered in the group that received fish oil compared to the control oil in this study [27]. This supports that responders had altered pathways related to apoptosis and immune function compared to non-responders in this sub-study.

Pathways related to LPA signaling were enriched in responders compared to non-responders to omega-3 supplementation. LPA is a lipid containing a fatty acid that may vary in length and degree of unsaturation. LPA can be produced from membrane phospholipids; it binds and activates LPA receptors (LPAR) and different forms of LPA may differentially affect LPA signaling [30]. In addition to their role in development of the central nervous system, LPARs are expressed in lymphocytes where they affect cytokine secretion, chemotaxis, and proliferation, and LPA signaling may be involved in the development of atherosclerosis and adipocyte differentiation [30,31,32]. We found that responders and non-responders had different levels of various plasma fatty acids that may reflect different plasma levels of various chemical forms of LPA [30]. Similarly, Rudkowska et al. found that compared to non-responders, responders had a greater increase in unsaturated fatty acids in glycerophosphatidylcholines, which can be used as LPA precursors [9, 30].

Transcripts with binding sites for HNF4-α were overrepresented in the list of gene transcripts differentially altered in responders compared to non-responders. Long-chain PUFAs have been shown to suppress HNF4-α activity [33]. HNF4-α is a TF that induces the expression of genes involved in lipid metabolism, and a decreased HNF4-α activity has been suggested to explain decreases in serum TG levels [34]; hence, an altered activity of HNF4-α in responders may be involved in the TG-lowering effect of omega-3 fatty acids.

The baseline level and the change in plasma omega-3 fatty acid levels did not differ between responders and non-responders, indicating that omega-3 fatty acid levels in plasma are not important for the TG-lowering effect. Indeed, omega-3 fatty acids mediate many of their biological effects after being incorporated into plasma membrane phospholipids [35]. Thus, the level of omega-3 fatty acids in red blood cells (RBC), also called the omega-3 index, that indicate the level of biologically available omega-3 fatty acids may be a better measure. We report that responders had a lower estimated omega-3 index than non-responders at baseline. This may imply that responders had a lower habitual fish intake before the start of the intervention than non-responders. In line with this, only participants with a low habitual fish intake had a reduced risk of major cardiovascular events in the VITAL trial [36]. Finally, only responders had an increase in RBC DHA in the study by Rudkowska et al., indicating that the level of omega-3 fatty acids in cell membranes may be important for the TG-lowering response [9]. Coupled with the higher baseline TG levels in responders, the lower omega-3 index at baseline may explain why responders lowered their TG levels after omega-3 supplementation.

At baseline, responders had higher plasma OA levels and lower plasma LA levels than non-responders, which may reflect different dietary patterns at baseline. This may have affected the ability of omega-3 fatty acids to lower TG levels. To further understand the variation of the TG-lowering effect after omega-3 supplementation, the role of other nutrients, such as other fatty acids, in modulating this effect should be elucidated in future studies.

As expected, the TG reduction in responders was reflected by reduced TG levels in almost all VLDL subclasses. However, there was also an unexpected lowering of HDL-TG in the XL-, M-, and S-HDL subclasses. Changes in HDL-TG levels in the current study may indicate that cholesteryl ester transfer protein (CETP) could be involved, as CETP facilitates the transfer of cholesteryl esters from HDL to VLDL and LDL in the exchange for TG. Although we did not find a difference between responders and non-responders in the change in CETP expression, CETP activity has been shown to increase after omega-3 intake [37]. In future TG responder studies, measurements of CETP activity may provide valuable insight into possible mechanisms behind the TG-lowering effect of omega-3 fatty acids.

The classification of subjects into responders and non-responders was based on previous studies that have suggested and employed a classification of responders as those with a TG reduction, while subjects who increase or do not alter their TG levels have been classified as non-responders [7, 9]. To be more confident that the group we defined as responders actually reduced their TG levels as a response to increased omega-3 intake, we chose to only include responders with a clinically relevant TG reduction. Day-to-day variation in TG levels has been reported to be 20%; hence, a TG reduction greater than 20% in responders in this study is likely to be an effect beyond day-to-day variations [38]. Another approach would be to use the actual day-to-day TG variation in the control group in our previous study to define responders and non-responders. However, the day-to-day variation in the control group was as high as 40%, which would result in only three participants who would be defined as responders, and the statistical power of this sub-study would be low. Nonetheless, the use of actual TG variation may be useful in future studies. Moreover, the baseline TG level in this study was low, 1.5 mmol/L in responders and 0.8 mmol/L in non-responders. The TG reduction observed after omega-3 supplementation depends on the baseline TG level [3]; hence, the low baseline TG level in this study may have resulted in a lower percentage of responders compared to other studies [6, 8, 39]. Future studies investigating TG responders to omega-3 intake should include participants with high enough baseline TG levels to ensure a clinically relevant TG response.

This study is limited by the low baseline TG levels and the low number of participants. As we included fewer subjects in this sub-study compared to the main study, this study may be underpowered. In addition, no gene expression changes were significant after adjusting for multiple testing. Hence, we did not expect any detectable differences between responders and non-responders in analyses of single genes using qPCR. Therefore, the microarray results in this study are not validated by qPCR. Furthermore, an FDR limit of 25% for the enriched pathways implies that 1 in 4 of the significant results is a false positive. However, a high FDR limit was chosen to avoid losing interesting results. Hence, the results in this study needs to be interpreted with caution and validated in a study designed and powered to investigate differences between TG responders and non-responders to omega-3 fatty acid supplementation. Despite these limitations, this exploratory study found differences between responders and non-responders, and it shows that this study design may be useful for investigating the underling mechanisms of the TG-lowering effect of omega-3 fatty acids.

Conclusions

TG responders and non-responders to omega-3 supplementation have different baseline lipoprotein subclass and PBMC gene expression profiles. Furthermore, they differentially alter their lipoprotein subclass and PBMC gene expression profiles. The differentially altered PBMC gene expression may partially explain the variability in TG response to omega-3 intake.

Methods

Study design and subjects

In this exploratory study, we used data from a previous double-blind randomized controlled parallel-group trial that aimed to investigate health effects of fish oils with different qualities [40]. The study was performed at Akershus University College in 2009. Healthy and non-smoking men and women aged 18–50 years with BMI < 30 kg/m2, TG ≤ 4 mmol/L, total cholesterol ≤ 7.5 mmol/L, C-reactive protein ≤ 10 mg/L, glucose ≤ 6 mmol/L, and blood pressure < 160/100 mmHg were included. Participants were stratified by gender and randomized in a 1:1:1 ratio to receive 1.6 g EPA + DHA per day from fish oil or oxidized fish oil or 0 g EPA + DHA per day from high-oleic sunflower oil. In total, 54 participants completed the 7-week intervention. Four weeks prior to the baseline visit and throughout the study, participants were instructed to avoid consumption of fish, omega-3 supplements, or foods enriched in omega-3 fatty acids. The first 3 weeks of the intervention was a fully controlled isoenergetic diet period, and for the remaining 4 weeks, subjects returned to their habitual diet, still avoiding consumption of omega-3 fatty acids. We have previously described subject characteristics, protocol, blinding, compliance and side effects, study products, and the fully controlled diet period in detail [40].

In the current sub-study, we combined data from the two groups that received fish oil (n = 35). From baseline to the 7-week visit, participants increased their intake of omega-3 fatty acids while their remaining diet was unchanged. Hence, we used data from these two visits to investigate differences in TG responders and non-responders to omega-3 fatty acids. In other studies, TG responders have been defined as all subjects with a TG reduction and non-responders as subjects who have no change or an increase in TG levels following omega-3 supplementation [7, 9]. Here, we define responders as subjects having a larger reduction in fasting TG than the 20% day-to-day variation [38] (n = 8) and non-responders as having a TG change between − 20 and 20% (n = 16). Because we wanted to compare participants with a TG reduction with those with no change in TG levels following omega-3 supplementation, participants with a higher TG increase than the 20% day-to-day variation were excluded from analyses in this study (n = 11).

This study was performed according to the guidelines laid down by the Declaration of Helsinki, and all procedures involving human subjects were approved by the Regional Committee for Medical Ethics (approval no. 6.2008.2215) and the Norwegian Social Science Data Services (approval no. 21924). All subjects provided written informed consent, and the study was registered at www.clinicaltrials.gov (ID no. NCT01034423).

Blood sampling and routine laboratory analyses

Blood samples were drawn at the baseline and 7-week visits after an overnight fast (≥ 12 h). Participants were instructed to avoid alcohol consumption and strenuous physical activity the day before blood sampling. Whole blood was collected in ethylenediaminetetraacetic acid (EDTA) tubes that were kept at room temperature for maximum 48 h. Serum was obtained from silica gel tubes that were kept at room temperature for at least 30 min before centrifugation (1500g, 12 min), and plasma was obtained from EDTA tubes that were immediately placed on ice and centrifuged within 10 min (1500g, 4 °C, 10 min). Routine laboratory analyses, such as serum TG, low-density lipoprotein cholesterol (LDL-C), and high-density lipoprotein cholesterol (HDL-C), as well as white blood cell count in whole blood were performed at a clinical routine laboratory (Fürst Medical Laboratory, Oslo, Norway). Plasma fatty acids were extracted by the Bligh and Dyer method [41] as previously described [40]. The level of plasma fatty acids is expressed as the percentage of total plasma fatty acids.

Dried blood spot omega-3 index

Whole blood from EDTA tubes were transferred to a filter paper that was dried, sealed in plastic bags, and stored at − 80 °C until analysis. The level of omega-3 fatty acids in the dried blood spots was analyzed by Vitas analytical service, Oslo, Norway. The omega-3 index was estimated from the level of whole blood EPA and DHA using an equation derived by a Norwegian population by Vitas AS (omega-3 index = whole blood EPA + DHA (%) * 0.95 + 0.35).

Analysis of lipoprotein subclasses

The EDTA plasma concentrations of 14 different lipoprotein subclass particles and their lipid constituents, including TG levels, were measured with a commercially available NMR platform (Nightingale Health Ltd). The different lipoprotein subclasses were defined based on their average diameter: extremely large very low-density lipoprotein (VLDL) with a possible contribution of chylomicrons (> 75 nm); extra large (XL-), large (L-), small (S-), and extra small (XS-) VLDL (64.0, 53.6, 44.5, 36.8, and 31.3 nm); intermediate-density lipoprotein (IDL, 28.6 nm); L-, M-, and S-LDL (25.5, 23.0 and 18.7 nm); and XL-, L-, M-, and S-HDL (14.3, 12.1, 10.9 and 8.7 nm). Details of the NMR metabolomics platform have previously been described [42].

Isolation of PBMC and RNA and microarray hybridization

PBMC were isolated using BD Vacutainer Cell Preparation tubes (Becton, Dickinson San Jose, CA, USA), a well-documented method for PBMC isolation with more than 90% purity. According to the manufacturer, about 80% of the isolated PBCMs are lymphocytes and 12% are monocytes. The PBMCs were isolated according to the manufacturer’s instructions and cell pellets were stored at − 80 °C. Total RNA was isolated using the Qiagen’s RNeasy Mini kit (Qiagen, Valencia, CA, USA) according to the manufacturer’s instructions. RNA quality and quantity were measured with the NanoDrop ND-1000 Spectrophotometer (Thermo Fisher Scientific, Gothenburg, Sweden) and Agilent Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA). The average RNA integrity number (RIN) was 9.6. Gene expression was analyzed by hybridization to an Illumina HumanHT-12 v4 Expression BeadChip that was scanned on an Illumina HiSCan microarray scanner (Illumina Inc., San Diego, CA, USA). Illumina GenomeStudio was used to transform bead-level data to probe-level intensities, which were exported raw for bioinformatics analysis as previously reported [27].

Analyses of microarray data

The Illumina intensities were quantile normalized, and probes without a detectable expression in at least 10% of the samples were excluded (detection P > 0.01). Of the 48,000 probes, 21,236 probes on the Illumina array were defined as expressed. A more detailed protocol of the microarray analyses has previously been reported [27]. We calculated the change in gene expression as the log2 ratio between intensities at baseline and after 7 weeks. The differences in gene expression changes between responders and non-responders were tested with multiple regression analyses adjusted for age and gender. Differentially expressed genes (p ≤ 0.05) were subjected to pathway analyses and transcription factor (TF) analyses using MetaCore (GeneGo, division of Thomson Reuters, St Joseph, MI, USA), and enriched pathways with a false discovery rate (FDR) < 25% were considered significant.

Other statistical analyses

Sample size calculation for the main study was based on the expected change in plasma omega-3 levels and has previously been described [40]. Lipoprotein subclass TG levels were log2-transformed before baseline differences and the difference in change from baseline (log2 ratio) between responders and non-responders were tested with a t test. Differences in categorical data were tested with Fisher’s exact test. All other data are presented as median and 25th–75th percentiles. The differences from baseline to the 7-week visit within a group were tested with a paired Wilcoxon signed rank test and the differences between responders and non-responders at baseline and the differences in change from baseline to the 7-week visit were tested with a Mann-Whitney U test. All statistical analyses were performed in R [43].

Abbreviations

AA:

Arachidonic acid

ACACA:

Acetyl-CoA carboxylase α

ACLY:

ATP citrate lyase

ACOX1:

Acyl-CoA oxidase 1

ALA:

α-Linolenic acid

APOE:

Apolipoprotein E

BMI:

Body mass index

C:

Cholesterol

CD36:

Cluster of differentiation 36

CETP:

Cholesteryl ester transfer protein

CVD:

Cardiovascular disease

DGAT2:

Diacylglycerol O-acyltransferase 2

DHA:

Docosahexaenoic acid

EDTA:

Ethylenediaminetetraacetic acid

EPA:

Eicosapentaenoic acid

FDR:

False discovery rate

FO:

Fish oil

GWAS:

Genome-wide association study

HDL:

High-density lipoprotein

HNF4A:

Hepatic nuclear factor 4α

IDL:

Intermediate-density lipoprotein

L-:

Large

LA:

Linoleic acid

LDL:

Low-density lipoprotein

LPA:

Lysophosphatididic acid

LPAR:

Lysophosphatididic acid receptor

M-:

Medium

NF-κB:

Nuclear factor kappa B

NMR:

Nuclear magnetic resonance

OA:

Oleic acid

oxFO:

Oxidized fish oil

PBMC:

Peripheral blood mononuclear cells

PPAR:

Peroxisome proliferator-activated receptors

RBC:

Red blood cell

RXRA:

Retinoid X receptor α

S-:

Small

SREBP:

Sterol regulating element binding protein

TF:

Transcription factor

TG:

Triglycerides

TGF-beta:

Transforming growth factor β

VLDL:

Very low-density lipoprotein

XL-:

Extra large

XS-:

Extra small

XXL-:

Extremely large

References

  1. Alexander DD, Miller PE, Van Elswyk ME, Kuratko CN, Bylsma LC. A meta-analysis of randomized controlled trials and prospective cohort studies of eicosapentaenoic and docosahexaenoic long-chain omega-3 fatty acids and coronary heart disease risk. Mayo Clin Proc. 2017;92(1):15–29. https://doi.org/10.1016/j.mayocp.2016.10.018.

    Article  CAS  PubMed  Google Scholar 

  2. AbuMweis S, Jew S, Tayyem R, Agraib L. Eicosapentaenoic acid and docosahexaenoic acid containing supplements modulate risk factors for cardiovascular disease: a meta-analysis of randomised placebo-control human clinical trials. J Hum Nutr Diet. 2017; https://doi.org/10.1111/jhn.12493.

  3. Balk EM, Lichtenstein AH, Chung M, Kupelnick B, Chew P, Lau J. Effects of omega-3 fatty acids on serum markers of cardiovascular disease risk: a systematic review. Atherosclerosis. 2006;189(1):19–30. https://doi.org/10.1016/j.atherosclerosis.2006.02.012.

    Article  CAS  PubMed  Google Scholar 

  4. Calder PC. Omega-3 polyunsaturated fatty acids and inflammatory processes: nutrition or pharmacology? Br J Clin Pharmacol. 2012;75(3):645–62. https://doi.org/10.1111/j.1365-2125.2012.04374.x.

    Article  CAS  Google Scholar 

  5. Geleijnse JM, Giltay EJ, Grobbee DE, Donders AR, Kok FJ. Blood pressure response to fish oil supplementation: metaregression analysis of randomized trials. J Hypertens. 2002;20(8):1493–9.

    Article  CAS  PubMed  Google Scholar 

  6. Caslake MJ, Miles EA, Kofler BM, Lietz G, Curtis P, Armah CK, et al. Effect of sex and genotype on cardiovascular biomarker response to fish oils: the FINGEN Study. Am J Clin Nutr. 2008;88(3):618–29.

    Article  CAS  PubMed  Google Scholar 

  7. Madden J, Williams CM, Calder PC, Lietz G, Miles EA, Cordell H, et al. The impact of common gene variants on the response of biomarkers of cardiovascular disease (CVD) risk to increased fish oil fatty acids intakes. Annu Rev Nutr. 2011;31:203–34. https://doi.org/10.1146/annurev-nutr-010411-095239.

    Article  CAS  PubMed  Google Scholar 

  8. Rundblad A, Holven KB, Bruheim I, Myhrstad MC, Ulven SM. Effects of krill oil and lean and fatty fish on cardiovascular risk markers: a randomised controlled trial. J Nutr Sci. 2018;7:e3. https://doi.org/10.1017/jns.2017.64.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Rudkowska I, Paradis AM, Thifault E, Julien P, Barbier O, Couture P, et al. Differences in metabolomic and transcriptomic profiles between responders and non-responders to an n-3 polyunsaturated fatty acids (PUFAs) supplementation. Genes Nutr. 2013;8(4):411–23. https://doi.org/10.1007/s12263-012-0328-0.

    Article  CAS  PubMed  Google Scholar 

  10. Bouchard-Mercier A, Rudkowska I, Lemieux S, Couture P, Vohl MC. Polymorphisms, de novo lipogenesis, and plasma triglyceride response following fish oil supplementation. J Lipid Res. 2013;54(10):2866–73. https://doi.org/10.1194/jlr.M041590.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Bouchard-Mercier A, Rudkowska I, Lemieux S, Couture P, Vohl MC. Polymorphisms in genes involved in fatty acid beta-oxidation interact with dietary fat intakes to modulate the plasma TG response to a fish oil supplementation. Nutrients. 2014;6(3):1145–63. https://doi.org/10.3390/nu6031145.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Madden J, Carrero JJ, Brunner A, Dastur N, Shearman CP, Calder PC, et al. Polymorphisms in the CD36 gene modulate the ability of fish oil supplements to lower fasting plasma triacyl glycerol and raise HDL cholesterol concentrations in healthy middle-aged men. Prostaglandins Leukot Essent Fat Acids. 2008;78(6):327–35. https://doi.org/10.1016/j.plefa.2008.04.003.

    Article  CAS  Google Scholar 

  13. Minihane AM, Khan S, Leigh-Firbank EC, Talmud P, Wright JW, Murphy MC, et al. ApoE polymorphism and fish oil supplementation in subjects with an atherogenic lipoprotein phenotype. Arterioscler Thromb Vasc Biol. 2000;20(8):1990–7.

    Article  CAS  PubMed  Google Scholar 

  14. Rudkowska I, Guenard F, Julien P, Couture P, Lemieux S, Barbier O, et al. Genome-wide association study of the plasma triglyceride response to an n-3 polyunsaturated fatty acid supplementation. J Lipid Res. 2014;55(7):1245–53. https://doi.org/10.1194/jlr.M045898.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Vallee Marcotte B, Cormier H, Guenard F, Rudkowska I, Lemieux S, Couture P, et al. Novel genetic loci associated with the plasma triglyceride response to an omega-3 fatty acid supplementation. J Nutrigenet Nutrigenomics. 2016;9(1):1–11. https://doi.org/10.1159/000446024.

    Article  CAS  PubMed  Google Scholar 

  16. Lee JN, Kim H, Yao H, Chen Y, Weng K, Ye J. Identification of Ubxd8 protein as a sensor for unsaturated fatty acids and regulator of triglyceride synthesis. Proc Natl Acad Sci U S A. 2010;107(50):21424–9. https://doi.org/10.1073/pnas.1011859107.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Botolin D, Wang Y, Christian B, Jump DB. Docosahexaneoic acid (22:6,n-3) regulates rat hepatocyte SREBP-1 nuclear abundance by Erk- and 26S proteasome-dependent pathways. J Lipid Res. 2006;47(1):181–92. https://doi.org/10.1194/jlr.M500365-JLR200.

    Article  CAS  PubMed  Google Scholar 

  18. Oh DY, Talukdar S, Bae EJ, Imamura T, Morinaga H, Fan W, et al. GPR120 is an omega-3 fatty acid receptor mediating potent anti-inflammatory and insulin-sensitizing effects. Cell. 2010;142(5):687–98. https://doi.org/10.1016/j.cell.2010.07.041.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Jump DB, Tripathy S, Depner CM. Fatty acid-regulated transcription factors in the liver. Annu Rev Nutr. 2013;33:249–69. https://doi.org/10.1146/annurev-nutr-071812-161139.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Afman L, Milenkovic D, Roche HM. Nutritional aspects of metabolic inflammation in relation to health--insights from transcriptomic biomarkers in PBMC of fatty acids and polyphenols. Mol Nutr Food Res. 2014;58(8):1708–20. https://doi.org/10.1002/mnfr.201300559.

    Article  CAS  PubMed  Google Scholar 

  21. de Mello VD, Kolehmanien M, Schwab U, Pulkkinen L, Uusitupa M. Gene expression of peripheral blood mononuclear cells as a tool in dietary intervention studies: what do we know so far? Mol Nutr Food Res. 2012;56(7):1160–72. https://doi.org/10.1002/mnfr.201100685.

    Article  CAS  PubMed  Google Scholar 

  22. Bouwens M, Afman LA, Muller M. Fasting induces changes in peripheral blood mononuclear cell gene expression profiles related to increases in fatty acid beta-oxidation: functional role of peroxisome proliferator activated receptor alpha in human peripheral blood mononuclear cells. Am J Clin Nutr. 2007;86(5):1515–23.

    Article  CAS  PubMed  Google Scholar 

  23. Caimari A, Oliver P, Keijer J, Palou A. Peripheral blood mononuclear cells as a model to study the response of energy homeostasis-related genes to acute changes in feeding conditions. Omics. 2010;14(2):129–41. https://doi.org/10.1089/omi.2009.0092.

    Article  CAS  PubMed  Google Scholar 

  24. Liew CC, Ma J, Tang HC, Zheng R, Dempsey AA. The peripheral blood transcriptome dynamically reflects system wide biology: a potential diagnostic tool. J Lab Clin Med. 2006;147(3):126–32. https://doi.org/10.1016/j.lab.2005.10.005.

    Article  CAS  PubMed  Google Scholar 

  25. Rudkowska I, Raymond C, Ponton A, Jacques H, Lavigne C, Holub BJ, et al. Validation of the use of peripheral blood mononuclear cells as surrogate model for skeletal muscle tissue in nutrigenomic studies. Omics. 2011;15(1–2):1–7. https://doi.org/10.1089/omi.2010.0073.

    Article  CAS  PubMed  Google Scholar 

  26. Bouwens M, van de Rest O, Dellschaft N, Bromhaar MG, de Groot LC, Geleijnse JM, et al. Fish-oil supplementation induces antiinflammatory gene expression profiles in human blood mononuclear cells. Am J Clin Nutr. 2009;90(2):415–24. https://doi.org/10.3945/ajcn.2009.27680.

    Article  CAS  PubMed  Google Scholar 

  27. Myhrstad MC, Ulven SM, Gunther CC, Ottestad I, Holden M, Ryeng E, et al. Fish oil supplementation induces expression of genes related to cell cycle, endoplasmic reticulum stress and apoptosis in peripheral blood mononuclear cells: a transcriptomic approach. J Intern Med. 2014;276(5):498–511. https://doi.org/10.1111/joim.12217.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Rundblad A, Holven KB, Bruheim I, Myhrstad MC, Ulven SM. Effects of fish and krill oil on gene expression in peripheral blood mononuclear cells and circulating markers of inflammation: a randomised controlled trial. J Nutr Sci. 2018;7:e10. https://doi.org/10.1017/jns.2018.2.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Jump DB. Fatty acid regulation of gene transcription. Crit Rev Clin Lab Sci. 2004;41(1):41–78. https://doi.org/10.1080/10408360490278341.

    Article  CAS  PubMed  Google Scholar 

  30. Yung YC, Stoddard NC, Chun J. LPA receptor signaling: pharmacology, physiology, and pathophysiology. J Lipid Res. 2014;55(7):1192–214. https://doi.org/10.1194/jlr.R046458.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Yung YC, Stoddard NC, Mirendil H, Chun J. Lysophosphatidic acid signaling in the nervous system. Neuron. 2015;85(4):669–82. https://doi.org/10.1016/j.neuron.2015.01.009.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Contos JJ, Ishii I, Chun J. Lysophosphatidic acid receptors. Mol Pharmacol. 2000;58(6):1188–96.

    Article  CAS  PubMed  Google Scholar 

  33. Hertz R, Magenheim J, Berman I, Bar-Tana J. Fatty acyl-CoA thioesters are ligands of hepatic nuclear factor-4alpha. Nature. 1998;392(6675):512–6. https://doi.org/10.1038/33185.

    Article  CAS  PubMed  Google Scholar 

  34. Hayhurst GP, Lee YH, Lambert G, Ward JM, Gonzalez FJ. Hepatocyte nuclear factor 4alpha (nuclear receptor 2A1) is essential for maintenance of hepatic gene expression and lipid homeostasis. Mol Cell Biol. 2001;21(4):1393–403. https://doi.org/10.1128/mcb.21.4.1393-1403.2001.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Calder PC. Marine omega-3 fatty acids and inflammatory processes: effects, mechanisms and clinical relevance. BBA – Mol Cell Biol Lipids. 2015;1851(4):469–84. https://doi.org/10.1016/j.bbalip.2014.08.010.

    Article  CAS  Google Scholar 

  36. Manson JE, Cook NR, Lee IM, Christen W, Bassuk SS, Mora S, et al. Marine n-3 fatty acids and prevention of cardiovascular disease and cancer. N Engl J Med. 2018; https://doi.org/10.1056/NEJMoa1811403.

  37. Pizzini A, Lunger L, Demetz E, Hilbe R, Weiss G, Ebenbichler C, et al. The role of omega-3 fatty acids in reverse cholesterol transport: a review. Nutrients. 2017;9(10) https://doi.org/10.3390/nu9101099.

  38. Bookstein L, Gidding SS, Donovan M, Smith FA. Day-to-day variability of serum cholesterol, triglyceride, and high-density lipoprotein cholesterol levels. Impact on the assessment of risk according to the National Cholesterol Education Program guidelines. Arch Intern Med. 1990;150(8):1653–7.

    Article  CAS  PubMed  Google Scholar 

  39. Cormier H, Rudkowska I, Paradis AM, Thifault E, Garneau V, Lemieux S, et al. Association between polymorphisms in the fatty acid desaturase gene cluster and the plasma triacylglycerol response to an n-3 PUFA supplementation. Nutrients. 2012;4(8):1026–41. https://doi.org/10.3390/nu4081026.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Ottestad I, Vogt G, Retterstol K, Myhrstad MC, Haugen JE, Nilsson A, et al. Oxidised fish oil does not influence established markers of oxidative stress in healthy human subjects: a randomised controlled trial. Br J Nutr. 2012;108(2):315–26. https://doi.org/10.1017/s0007114511005484.

    Article  CAS  PubMed  Google Scholar 

  41. Bligh EG, Dyer WJ. A rapid method of total lipid extraction and purification. Can J Biochem Physiol. 1959;37(8):911–7. https://doi.org/10.1139/o59-099.

    Article  CAS  PubMed  Google Scholar 

  42. Soininen P, Kangas AJ, Wurtz P, Suna T, Ala-Korpela M. Quantitative serum nuclear magnetic resonance metabolomics in cardiovascular epidemiology and genetics. Circ Cardiovasc Genet. 2015;8(1):192–206. https://doi.org/10.1161/circgenetics.114.000216.

    Article  CAS  PubMed  Google Scholar 

  43. R Development Core Team. R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2016.

    Google Scholar 

Download references

Acknowledgements

The authors wish to thank the study subjects for their participation.

Funding

This study was funded by the University of Oslo, OsloMet – Oslo Metropolitan University, the research council of Norway, and the Throne-Holst Foundation for Nutrition Research, Oslo, Norway.

Author information

Authors and Affiliations

  1. Department of Nutrition, Institute of Basic Medical Sciences, University of Oslo, PO Box 1046, Blindern, 0317, Oslo, Norway

    Amanda Rundblad, Sunniva V. Larsen, Inger Ottestad, Kirsten B. Holven & Stine M. Ulven

  2. Department of Nursing and Health Promotion, Faculty of Health Sciences, OsloMet – Oslo Metropolitan University, PO Box 4, St Olavs plass, 0130, Oslo, Norway

    Mari C. Myhrstad

  3. Department of Biostatistics, Institute of Basic Medical Sciences, University of Oslo, PO Box 1046, Blindern, 0317, Oslo, Norway

    Magne Thoresen

  4. National Advisory Unit on Familial Hypercholesterolemia, Department of Endocrinology, Morbid Obesity and Preventive Medicine, Oslo University Hospital, PO Box 4950, Nydalen, 0424, Oslo, Norway

    Kirsten B. Holven

Authors
  1. Amanda Rundblad
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Sunniva V. Larsen
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Mari C. Myhrstad
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Inger Ottestad
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Magne Thoresen
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Kirsten B. Holven
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Stine M. Ulven
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

AR designed the research, performed the statistical analyses, interpreted the data, and drafted the manuscript. SVL designed the research, performed the statistical analyses, interpreted the data, and drafted the manuscript. MCM designed the research, performed the microarray analysis, interpreted the data, and drafted the manuscript. IO designed the research, coordinated the intervention study, and drafted the manuscript. MT contributed to the analysis and interpretation of the data and drafted the manuscript. KBH designed the research, interpreted the data, and drafted the manuscript. SMU designed the research, was responsible for the intervention study, interpreted the data, and drafted the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Amanda Rundblad.

Ethics declarations

Ethics approval and consent to participate

Given in the methods section.

Consent for publication

Not applicable.

Competing interests

Dr. Holven has received research grants and/or personal fees from Tine SA, Mills SA, Olympic Seafood, Kaneka, Amgen, Sanofi, and Pronova, none of which are related to the content of this manuscript. Dr. Ulven has received research grants from Tine DA, Mills DA, and Olympic Seafood, none of which are related to the content of this manuscript. The other authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

List of genes with a significantly (p ≤ 0.05) different baseline expression in responders compared to non-responders adjusted for age and gender. (XLSX 671 kb)

Additional file 2:

List of significantly (p ≤ 0.05) differentially altered gene transcripts in responders compared to non-responders adjusted for age and gender. (XLSX 622 kb)

Additional file 3:

Enriched pathways (FDR < 25%) among the significantly altered transcripts in responders compared to non-responders. (XLSX 14 kb)

Additional file 4:

Results from the transcription factor analyses. (XLSX 16 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Rundblad, A., Larsen, S.V., Myhrstad, M.C. et al. Differences in peripheral blood mononuclear cell gene expression and triglyceride composition in lipoprotein subclasses in plasma triglyceride responders and non-responders to omega-3 supplementation. Genes Nutr 14, 10 (2019). https://doi.org/10.1186/s12263-019-0633-y

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12263-019-0633-y

Keywords

Genes & Nutrition

ISSN: 1865-3499

Contact us