Skip to main content
  • Research article
  • Open access
  • Published:

Association between plasma phospholipid saturated fatty acids and metabolic markers of lipid, hepatic, inflammation and glycaemic pathways in eight European countries: a cross-sectional analysis in the EPIC-InterAct study

Abstract

Background

Accumulating evidence suggests that individual circulating saturated fatty acids (SFAs) are heterogeneous in their associations with cardio-metabolic diseases, but evidence about associations of SFAs with metabolic markers of different pathogenic pathways is limited. We aimed to examine the associations between plasma phospholipid SFAs and the metabolic markers of lipid, hepatic, glycaemic and inflammation pathways.

Methods

We measured nine individual plasma phospholipid SFAs and derived three SFA groups (odd-chain: C15:0 + C17:0, even-chain: C14:0 + C16:0 + C18:0, and very-long-chain: C20:0 + C22:0 + C23:0 + C24:0) in individuals from the subcohort of the European Prospective Investigation into Cancer and Nutrition (EPIC)-InterAct case-cohort study across eight European countries. Using linear regression in 15,919 subcohort members, adjusted for potential confounders and corrected for multiple testing, we examined cross-sectional associations of SFAs with 13 metabolic markers. Multiplicative interactions of the three SFA groups with pre-specified factors, including body mass index (BMI) and alcohol consumption, were tested.

Results

Higher levels of odd-chain SFA group were associated with lower levels of major lipids (total cholesterol (TC), triglycerides, apolipoprotein A-1 (ApoA1), apolipoprotein B (ApoB)) and hepatic markers (alanine transaminase (ALT), aspartate transaminase (AST), gamma-glutamyl transferase (GGT)). Higher even-chain SFA group levels were associated with higher levels of low-density lipoprotein cholesterol (LDL-C), TC/high-density lipoprotein cholesterol (HDL-C) ratio, triglycerides, ApoB, ApoB/A1 ratio, ALT, AST, GGT and CRP, and lower levels of HDL-C and ApoA1. Very-long-chain SFA group levels showed inverse associations with triglycerides, ApoA1 and GGT, and positive associations with TC, LDL-C, TC/HDL-C, ApoB and ApoB/A1. Associations were generally stronger at higher levels of BMI or alcohol consumption.

Conclusions

Subtypes of SFAs are associated in a differential way with metabolic markers of lipid metabolism, liver function and chronic inflammation, suggesting that odd-chain SFAs are associated with lower metabolic risk and even-chain SFAs with adverse metabolic risk, whereas mixed findings were obtained for very-long-chain SFAs. The clinical and biochemical implications of these findings may vary by adiposity and alcohol intake.

Peer Review reports

Background

There has been ongoing interest in the association between circulating saturated fatty acids (SFAs) and cardio-metabolic diseases such as cardiovascular diseases and type 2 diabetes. Recent evidence has highlighted that different individual circulating SFAs have discordant associations with cardio-metabolic diseases [1,2,3,4,5,6,7]. For example, we and others reported that circulating odd-chain SFAs (C15:0 (pentadecanoic acid) and C17:0 (heptadecanoic acid)) and very-long-chain SFAs (C20:0 (arachidic acid), C22:0 (behenic acid), C23:0 (tricosanoic acid) and C24:0 (lignoceric acid)) were inversely associated with incident type 2 diabetes, while circulating even-chain SFAs (C14:0 (myristic acid), C16:0 (palmitic acid) and C18:0 (stearic acid)) were positively associated [1,2,3,4]. In a meta-analysis, circulating C17:0 was inversely associated with coronary heart disease, while no association was found for any of C14:0, C15:0, C16:0 or C18:0 [5]. Other studies showed that higher circulating C16:0 was associated with higher risk of atrial fibrillation and heart failure, whereas other SFAs (C18:0, C20:0, C22:0 and C24:0) showed inverse or non-significant associations with these cardiac outcomes [6, 7].

The pathogenesis of cardio-metabolic diseases involves various metabolic pathways and intermediate metabolic markers, including markers of lipid metabolism, liver function, chronic inflammation and glycaemic homeostasis [8]. Clarifying the associations of individual circulating SFAs with these metabolic markers in different pathways will help in the understanding of the roles of individual SFAs in disease aetiology. Previous studies of circulating SFAs with intermediate metabolic markers were inconsistent or not comprehensive and analysed using different lipid fractions (phospholipid fraction, erythrocyte membrane, etc.) [9,10,11,12,13,14,15,16,17]. For instance, odd-chain SFAs C15:0 and C17:0 from erythrocyte membranes were positively associated with low-density lipoprotein cholesterol (LDL-C) and C17:0 was inversely associated with high-density lipoprotein cholesterol (HDL-C) in men, while C15:0 was positively associated with HDL-C in women in a European study [11]. Both circulating odd-chain and even-chain SFAs showed heterogeneous associations with glycaemic markers, which were studied in different lipid fractions [3, 14, 15, 17,18,19]. Circulating very-long-chain SFAs (C20:0, C22:0 and C24:0) had mixed associations with different metabolic markers, including a positive association with LDL-C, an inverse association with triglycerides (TG) and an inverse association (C24:0 only) with C-reactive protein (CRP) and insulin resistance [2]. The inconsistency in the existing literature may reflect the heterogeneity of sample sizes, ranging from 38 [20] to 3004 participants [3]. Furthermore, fatty acids were measured using different methods in different lipid fractions across different studies, which might contribute to the inconsistent results in the existent literature. This adds to the rationale for measuring fatty acids in one single lipid fraction (such as circulating phospholipids) in a large multicentre study with diverse population characteristics to characterise robust associations of fatty acids with metabolic markers.

Therefore, evaluating fatty acid profiles in populations in eight European countries, we aimed to examine the associations of individual circulating plasma phospholipid SFAs and corresponding SFA groups with several metabolic markers of lipid metabolism, liver function, chronic inflammation and glycaemic control.

Methods

Study design and population

These analyses used data from the subcohort of the European Prospective Investigation into Cancer and Nutrition (EPIC)-InterAct study. EPIC-InterAct is a case-cohort study of type 2 diabetes nested within eight countries of the EPIC cohort study, namely France, Italy, Spain, UK, the Netherlands, Germany, Sweden and Denmark [21]. Briefly, from the entire cohort of 340,234 participants with blood samples (baseline years, 1991–1998) in the eight countries, we ascertained incident diabetes cases and selected a random subcohort of 16,835 participants. Participants with prevalent diabetes at baseline were excluded (n = 681) [21]. A further 235 participants with no available fatty acid data were also excluded. Thus, a total of 15,919 participants were included in this study. All participants provided written informed consent and the study was approved by the local ethics committee in the participating centres and the Internal Review Board of the International Agency for Research on Cancer.

Laboratory measurements

Plasma phospholipid fatty acids were measured between 2010 and 2012 at the Medical Research Council Human Nutrition Research (Cambridge, UK) from non-fasting plasma samples stored at baseline (1993–1998) at –196 °C (–150 °C in Denmark). We considered plasma phospholipid fatty acid compositions as stable because stability has been confirmed at –80 °C or below over long-term storage [22]. The laboratory staff were blind to any information about the participants. The assay method and quality control of the fatty acid measurement have been described elsewhere [23]. A total of 37 fatty acids were identified by comparing their retention time with those of commercial standards. All the fatty acids were expressed as percentage of total phospholipid fatty acid (mol%). Nine SFAs had a relative concentration higher than 0.05%. Data quality was assessed using quality control samples across all of the batches, which showed that the SFAs and their coefficients of variation were C14:0 (9.4%), C15:0 (11.9%), C16:0 (1.6%), C17:0 (4.2%), C18:0 (2.0%), C20:0 (15.3%), C22:0 (10.3%), C23:0 (18.9%) and C24:0 (14.7%).

Serum metabolic markers were measured at Stichting Ingenhousz Laboratory (Etten-Leur, Netherlands), from samples stored at –196 °C (–80 °C in Umeå). These metabolic markers included lipid markers (total cholesterol (TC), HDL-C TG, apolipoprotein A-1 (ApoA1) and apolipoprotein B (ApoB)), liver function markers (alanine transaminase (ALT), aspartate transaminase (AST), gamma glutamyl transferase (GGT)), and an inflammatory marker (C-reactive protein (CRP)). All assays were performed using a Cobas® (Roche Diagnostics, Mannheim, Germany) assay on a Roche Hitachi Modular P analyser. Several other lipid markers were calculated based on the measured lipids, namely LDL-C based on the Friedewald formula [24], TC/HDL-C (a ratio of TC to HDL-C) and ApoB/A1 (a ratio of ApoB to ApoA1). A marker of glycaemic control, haemoglobin A1c (HbA1c), was measured at Stichting Ingenhousz Laboratory in the erythrocyte fraction from samples stored at –196 °C (–80 °C in Umeå) using the Tosoh-G8 analyser (Tosoh Bioscience, Japan). The assay detection ranges for the metabolic marker measurements were TC (0.08–20.7 mmol/L), HDL-C (0.08–3.1 mmol/L), TG (0.05–11.3 mmol/L), ApoA1 (7.14–143 μmol/L), ApoB (0.39–7.8 μmol/L), ALT (4–600 U/L), AST (4–800 U/L), GGT (3–1200 U/L), HbA1c (20–140 mmol/mol) and CRP (0.1–20 mg/L). Quality control was based on the Westgard rules [25]. All the metabolic markers were winsorised based on the 0.1 and 99.9th percentiles to minimise influence of potential outliers.

Diet and lifestyle measurements

The assessment of habitual diet during the past 12 months was undertaken using self- or interviewer-administered, country-specific validated dietary questionnaires [26]. Based on the standardised EPIC Nutrient Database [27], total energy and nutrient intake were calculated. At baseline, we used standardised health and lifestyle questionnaires to assess demographics, smoking status, medical history and educational level. Physical activity was assessed using a brief questionnaire at baseline, which divided participants into four physical activity groups, namely inactive, moderately inactive, moderately active or active [21].

Statistical analysis

We used Stata version 14 for all analyses. We analysed nine individual SFAs (C14:0, C15:0, C16:0, C17:0, C18:0, C20:0, C22:0, C23:0, C24:0) and three SFA groups – odd-chain SFA (sum of 15:0 and 17:0), even-chain SFA (sum of 14:0, 16:0 and 18:0) and very-long-chain SFA (sum of C20:0, C22:0, C23:0, C24:0). The groupings of the SFAs were based on those used in a previous analysis of EPIC-InterAct [1]. Spearman correlation coefficients between the individual SFAs were calculated. As dependent variables, we evaluated eight lipid markers (TC, HDL-C, TC/HDL-C, LDL-C, TG, ApoB, ApoA1, ApoB/A1), three liver function markers (ALT, AST and GGT), one glycaemic marker (HbA1c) and one inflammatory marker (CRP). Variables with skewed distribution were log-transformed before analysis (TC/HDL-C, TG, ALT, AST, GGT and CRP). In the primary analysis (based on complete case analysis), we used linear regression to estimate country-specific associations per 1 standard deviation (SD, calculated in the overall subcohort) of each SFA with each of the 13 metabolic markers, and combined the estimated associations across countries using random effects meta-analysis. Associations were expressed as the difference in metabolic marker (in SD units) per 1 SD difference in each fatty acid. A P value threshold for significance of 0.00068 was defined based on the number of tests performed, accounting for correlations among exposures and outcomes, calculated using the method of Li and Ji [28]. Confidence intervals (CI) were also estimated, which were consistent with this P value threshold.

Using linear regression, we adjusted for the potential confounders as follows: Model 1 – age, sex and centre; Model 2 – as Model 1 plus body mass index (BMI); Model 3 – as Model 2 plus physical activity (inactive, moderately inactive, moderately active or active), smoking status (never, former or current), alcohol drinking (never, 0 to < 6, 6 to < 12, 12 to < 24, and ≥ 24 g/day), educational level (none, primary school, technical or professional school, secondary school, or higher education) and total energy intake (kcal/day). Model 3 was extended to test the potential interaction of the three SFA groups (even-chain, odd-chain and very-long-chain SFAs) with age, sex, BMI, alcohol drinking and physical activity. For these pre-specified interaction tests on the multiplicative scale to investigate potential biological variation, a P value threshold of 0.0004 was used, calculated using the same method as above [28].

We performed sensitivity analyses based on Model 3, where (3a) included a variety of dietary factors as covariates (dietary intakes of carbohydrates, dairy, red and processed meat, fruits and vegetables, olive oil and vegetable oil); (3b) excluded participants with HbA1c ≥ 6.5%; (3c) included self-reported heart disease, stroke and cancer as covariates; (3d) excluded participants with self-reported heart disease, stroke or cancer; (3e) excluded people with hyperlipidaemia (self-reported hyperlipidaemia or treatment for hyperlipidaemia); (3f) included hyperlipidaemia status as a covariate; (3g) was a mutual adjustment for the other SFA groups; (3h) excluded participants with very heavy alcohol drinking (>95th percentile of the present population (>50.2 g/d)); (3i) excluded participants with CRP ≥ 95.2 nmol/L (equivalent to 10 mg/L); and (3j) replaced BMI with waist circumference as a marker of adiposity.

Results

Population characteristics

The mean age of the InterAct subcohort participants was 52.3 years (SD 9.2 years, range 20–77 years). Participants who were men, older or with higher BMI had relatively higher concentrations of even-chain SFAs. Women or those with lower BMI had higher odd-chain and very-long-chain SFA levels (Table 1). In addition, women had a more favourable metabolic profile than men, such as higher levels of HDL-C and lower levels of TC/HDL-C, TG and liver function markers. Older participants and those with a higher BMI had higher levels of TC, LDL-C, TC/HDL-C, TG, ApoB, ApoB/A1, ALT, AST, GGT, HbA1c and CRP. The distributions of these metabolic markers varied across different participating countries and across different SFA groups (Additional file 1: Table S1, Table S2). The odd-chain SFA group was negatively correlated with the even-chain SFA group, while being positively correlated with a very-long-chain SFA group (Additional file 1: Table S3). The correlations of the three SFA groups with various food sources in InterAct have been reported previously [1].

Table 1 Distribution of plasma metabolic biomarkers and saturated fatty acid groups in the subcohort (n = 15919) and by sex, age and BMI: EPIC-InterAct study

Association of plasma phospholipid SFAs with metabolic markers

Odd-chain SFAs

Both C15:0 and C17:0 were inversely associated with several lipid (TC, TG, ApoA1, ApoB) and liver function markers (ALT, AST and GGT). In addition, C17:0 was inversely associated with TC/HDL-C and CRP (Figs. 1, 2, 3, 4 and 5, Additional file 1: Table S4). The odd-chain SFA group (C15:0 + C17:0) was inversely associated with most metabolic markers; there was no evidence of an association with HDL-C, LDL-C, ApoB/A1 or HbA1c (Fig. 5).

Fig. 1
figure 1

Associations of plasma phospholipid saturated fatty acids (SFAs) with total and types of cholesterol: EPIC-InterAct study. Standardised difference is the difference (in SD units of the metabolic marker) per 1 SD of each SFA. The error bars represent confidence intervals, corrected for multiple testing based on a familywise error rate of 5%. * Indicates that the P value for the association was < 0.00068. Estimates were based on random effects meta-analysis of country-specific estimates from a linear regression adjusted for age, sex, centre, BMI, physical activity, smoking status, alcohol drinking, educational level and total energy intake. TC/HDL-C was log-transformed. TC total cholesterol, HDL-C high-density lipoprotein cholesterol, LDL-C low-density lipoprotein cholesterol

Fig. 2
figure 2

Associations of plasma phospholipid saturated fatty acids (SFAs) with triglycerides and apolipoproteins: EPIC-InterAct study. Standardised difference is the difference (in SD unit of the metabolic marker) per 1 SD of each SFA. The error bars represent confidence intervals, corrected for multiple testing based on a family-wise error rate of 5%. * Indicates that the P value for the association was < 0.00068. Estimates were based on random effects meta-analysis of country-specific estimates from a linear regression adjusted for age, sex, centre, BMI, physical activity, smoking status, alcohol drinking, educational level and total energy intake. ApoA1 apolipoprotein A-1, ApoB apolipoprotein B

Fig. 3
figure 3

Associations of plasma phospholipid saturated fatty acids (SFAs) with circulating liver enzymes: EPIC-InterAct study. Standardised difference is the difference (in SD unit of the metabolic marker) per 1 SD of each SFA. The error bars represent confidence intervals, corrected for multiple testing based on a family-wise error rate of 5%. * Indicates that the P value for the association was < 0.00068. Estimates were based on random effects meta-analysis of country-specific estimates from a linear regression adjusted for age, sex, centre, BMI, physical activity, smoking status, alcohol drinking, educational level and total energy intake. All the ALT, AST and GGT were log-transformed. ALT alanine transaminase, AST aspartate transaminase, GGT gamma glutamyl transferase

Fig. 4
figure 4

Associations of plasma phospholipid saturated fatty acids (SFAs) with haemoglobin A1c and C-reactive protein: EPIC-InterAct study. Standardised difference is the difference (in SD unit of the metabolic marker) per 1 SD of each SFA. The error bars represent confidence intervals, corrected for multiple testing based on a family-wise error rate of 5%. * Indicates that the P value for the association was < 0.00068. Estimates were based on random effects meta-analysis of country-specific estimates from a linear regression adjusted for age, sex, centre, BMI, physical activity, smoking status, alcohol drinking, educational level and total energy intake. CRP was log-transformed. TG triglycerides, HbA1c haemoglobin A1c

Fig. 5
figure 5

Summary of the associations of plasma phospholipid saturated fatty acids (SFAs) with metabolic markers: EPIC-InterAct study. Value in each cell represented standardised difference (in SD unit) of metabolic markers per 1 SD difference in each variable of plasma phospholipid SFAs. Values in bold (with *) indicate that the P value for the association was < 0.00068. All the values are expressed in red scale for different levels of positive associations and blue scale for different levels of negative associations

Even-chain SFAs

Circulating C14:0 was positively associated with circulating levels of TG, ApoA1 and ALT, but there was no association with the other metabolic markers. Circulating C16:0 was inversely associated with HDL-C, and positively associated with TG and CRP. C18:0 was inversely associated with HDL-C and ApoA1, and positively associated with LDL-C, TC/HDL-C, TG, ApoB, ApoB/A1, ALT, AST, GGT and HbA1c. The even-chain SFA group was significantly positively associated with 9 of the 13 metabolic markers (Figs. 1, 2, 3, 4 and 5, Additional file 1: Table S5).

Very-long-chain SFAs

Circulating C20:0 was inversely associated with TG, ApoA1 and GGT; while C22:0, C23:0 and C24:0 showed consistently positive associations with TC, LDL-C, TC/HDL-C, ApoB and ApoB/A1, and inverse associations with TG, ApoA1 and GGT (Figs. 1, 2, 3, 4 and 5, Additional file 1: Table S6). C23:0 was also positively associated with HbA1c (Fig. 4). As a group, very-long-chain SFAs had the strongest inverse association with TG, while other associations were similar to those of C24:0 (Fig. 5).

Sensitivity and interaction analyses

Most of the sensitivity analyses did not substantially affect the results (Additional file 1: Table S7). There was heterogeneity between countries in some of the SFAs and metabolic marker associations, such as LDL-C, TG and GGT, although the directions of the associations were generally consistent across countries (Additional file 1: Figure S1).

Results of interaction analyses showed that, in non-alcohol drinkers, there was no evidence of an association of the odd-chain SFA group with TC, ApoA1, ALT or AST, while in the higher alcohol consumption group, there was an inverse association of the odd-chain SFA group with TC, ApoA1, ALT, AST and GGT (Additional file 1: Figure S2). Similarly, a stronger inverse association of the very-long-chain SFA group with HDL-C and ApoA1 was found among high alcohol drinkers, with no significant association among non-drinkers. There was a suggestion that the association of each SFA group with each metabolic marker was stronger in those with higher BMI. In addition, the association of odd-chain SFAs with some metabolic markers (TC, AST, ALT and GGT) appeared to be stronger in men than in women.

Discussion

In a large cohort of European adults, we identified distinct patterns of association between different SFAs and metabolic markers. Briefly, odd-chain SFAs had favourable associations, while even-chain SFAs had adverse associations with most metabolic markers; there were mixed patterns of association with very-long-chain SFAs. We found that both odd-chain SFAs (C15:0 and C17:0) had consistent inverse associations with lipid markers (TC, TG, ApoA1 and ApoB) and liver function markers (ALT, AST and GGT). The odd-chain SFA C17:0 was additionally inversely associated with another lipid marker (TC/HDL-C ratio) and with the chronic inflammatory marker (CRP). In contrast, there were positive associations between the even-chain SFA group (sum of C14:0, C16:0 and C18:0) and several lipids (LDL-C, TC/HDL-C, TG, ApoB, ApoB/A1), liver function markers (ALT, AST, GGT) and CRP, and inverse associations with HDL-C and ApoA1. These associations were mainly driven by C18:0, except for CRP (only associated with C16:0). C14:0 and C16:0 showed heterogeneous associations with different markers. For very-long-chain SFAs (C20:0, C22:0, C23:0 and C24:0), there were heterogeneous findings – the group was inversely associated with two lipid markers (TG, ApoA1) and one liver function marker (GGT), but positively associated with other lipid markers (TC, LDL-C, TC/HDL-C, ApoB and ApoB/A1).

Eight of the 13 metabolic markers included in the analysis were lipid markers. Higher LDL-C and TG and lower HDL-C are known risk factors for cardiovascular diseases [29,30,31]. ApoA1 is abundant in anti-atherogenic HDL particles, while ApoB is part of the atherogenic non-HDL particles (e.g., LDL). Both ApoB/A1 and TC/HDL-C represent a balance of the anti-atherogenic and atherogenic lipoprotein particles, and are considered better predictors of the future cardiovascular disease risk than the individual lipids [32, 33]. Previous literature has reported heterogeneous findings on the association of SFAs with these lipid markers [2, 3, 11,12,13,14, 17]. We did not find an association between odd-chain SFAs and the ratios of ApoB/A1 or TC/HDL-C (after adjustment for the other SFA groups), suggesting that odd-chain SFAs might not be independently associated with anti-atherogenic properties. Instead, they may reflect the associations of the other SFA groups, as both even-chain and very-long-chain SFAs were positively associated with ApoB/A1 and TC/HDL-C in the present study and the three groups were correlated. There was a strong inverse association between very-long-chain SFAs and TG, which was consistent with the inverse association between very-long-chain SFAs and type 2 diabetes and cardiovascular diseases [1, 7]. The findings that very-long-chain SFAs had positive associations with atherogenic lipid markers (e.g. LDL-C, TC/HDL-C and ApoB/A1) were consistent with previous studies on LDL-C levels [2, 34]. In addition, very-long-chain SFAs in erythrocytes were inversely associated with LDL particle size [34], and smaller LDL particles are considered more atherogenic than larger LDL particles [35].

Several studies [12, 36, 37] examined the associations of SFAs with liver function markers, with mixed findings. Even-chain SFAs were positively associated with hepatic markers in one study [12], but not associated in several other studies [36, 37], with no evidence that odd-chain or very-long-chain SFAs were associated with these markers. In our study, odd-chain SFAs were inversely associated, C18:0 was positively associated with the three hepatic markers included, and very-long-chain SFAs were inversely associated with GGT. These findings were consistent with the direction of the association found for cardio-metabolic disease endpoints, including type 2 diabetes and cardiovascular diseases [1,2,3,4, 7, 13, 38].

The associations of SFAs with HbA1c and CRP were modest. Only C18:0 and C23:0 had weak positive associations with HbA1c. The positive association between C23:0 and HbA1c was contradictory to the inverse association of C23:0 with incident type 2 diabetes previously found in this population [1]. However, the relative concentration of C23:0 was lowest among the nine SFAs, and the current cross-sectional association with HbA1c may not be comparable to results from association with subsequent incident disease. One odd-chain SFA (C17:0) and one even-chain SFA (C16:0) were, respectively, inversely and positively associated with CRP, an inflammatory marker and risk factor of cardio-metabolic diseases [39]. The findings for CRP were consistent with several prior US studies in which none of C14:0, C15:0 or C18:0 were associated with CRP, while C16:0 was positively associated [3, 13]. Of note, previous studies only reported associations for a maximum of three of the individual SFAs with CRP [3, 13], while the present study provided a more comprehensive analysis across nine individual SFAs.

No studies have previously reported the interaction between alcohol or BMI and SFAs on metabolic outcomes. We observed consistent patterns of interaction for alcohol intake and BMI with SFA groups for the metabolic markers, with stronger SFAs and metabolic marker associations in groups with higher reported alcohol consumption or higher BMI. The mechanism behind these interactions is not clear. It may be that higher alcohol intake (which stimulates de novo lipogenesis) [40] or higher BMI is associated with higher absolute circulating SFA concentrations. Therefore, associations per SD SFA relative concentration might be larger among higher alcohol intake or BMI groups (with higher absolute SFA levels). Nevertheless, these interactions suggested that dietary exposure and adiposity should be taken into account in evaluating the relationship between individual SFAs and metabolic risk factors.

Results of the present study may have important biological implications linking SFAs and metabolic diseases. First, we found that odd-chain SFAs were favourably associated with markers of lipids (with exceptions for HDL-C and ApoA1), liver function and chronic inflammation, which we postulated might contribute to the inverse association of odd-chain SFA with cardio-metabolic diseases [1, 5]. It was recently shown that C15:0 and C17:0 may have different origins [41] and that C15:0 is a marker of exogenous origin dairy fat intake, while C17:0 is substantially biosynthesised endogenously. Another study suggested that both C15:0 and C17:0 are also biomarkers of dietary fibre intake [42]. These results suggest that the observed associations in this study indirectly reflect the impact of dietary dairy product and fibre intake on cardio-metabolic disease. Second, total even-chain SFA levels were positively associated with cardio-metabolic diseases potentially through all the tested pathways in this study, while individual even-chain SFAs appeared to have differential roles in these pathways. For example, only C16:0 was associated with the chronic inflammatory marker, CRP, in agreement with several in vitro studies [43,44,45]. The observed associations of even-chain SFAs with metabolic markers might partly reflect the effect of deregulation of de novo lipogenesis [46], as all of the three even-chain SFAs are within the de novo lipogenesis pathway, and these fatty acids were positively associated with risk of metabolic disorders [3, 12]. Third, very-long-chain SFAs were favourably associated with markers of lipids and liver function, which may be part of the mechanisms behind the inverse association of this group of SFAs with cardio-metabolic diseases, such as type 2 diabetes [1]. These observations are plausible as previous animal and cell culture studies [47] have suggested that ceramides containing very-long-chain SFAs show favourable effects on biological pathways (e.g. hepatocyte physiology and pathophysiology), compared with ceramides containing C16:0. In contrast, very-long-chain SFAs were positively associated with atherogenic lipid markers. Therefore, a TG or TG-related pathway may be a more reliable indicator linking plasma phospholipid very-long-chain SFAs and cardio-metabolic diseases compared to other lipid markers.

The strengths of our study include its large sample size of nearly 16,000 participants from eight European countries with different dietary/lifestyle backgrounds, and objective measurement of nine individual plasma phospholipid SFAs. The main limitation of this study is that the associations are cross-sectional, and therefore temporality of the associations is unclear and reverse causation cannot be ruled out. Additionally, as is common with all observational research, despite our attempts to adjust for a comprehensive range of potential confounders, residual confounding cannot be excluded. Future longitudinal studies relating changes in SFAs to change in metabolic markers as well as to the metabolome will be informative, and use of genetic Mendelian randomisation analysis and pharmacological or dietary interventions altering fatty acid profiles will help with causal inference. Another limitation is the possibility of false positive results among the large number of analyses performed, although we used a conservative approach to correct for multiple testing to minimise this possibility. Furthermore, potential measurement errors may have occurred because of the long-term storage of our plasma samples and may have biased our findings toward the null. Results from the phospholipid fatty acids may not be able to represent those from other lipid fractions, such as cholesterol and TG, and potential heterogeneity may exist across different lipid fractions. Finally, only one marker of inflammation (CRP) was available in our study, and future studies on this topic would benefit from an expanded panel of inflammatory cytokines, acute phase reactants and adipokines to further advance our understanding.

Conclusions

Our analyses suggest that odd-chain SFAs are inversely associated with several metabolic markers in pathways of lipid metabolism, liver function and chronic inflammation, thereby potentially contributing to their inverse associations with cardio-metabolic diseases. Further, even-chain SFAs are unfavourably associated with metabolic markers in all the tested pathways, including lipid metabolism, liver function, glycaemic control and chronic inflammation, with inconsistent associations for different individual even-chain SFAs. Very-long-chain SFAs are inversely associated with TG and markers of liver function, which potentially contributes to their inverse associations with cardio-metabolic diseases, while their positive associations with atherogenic lipid markers warrants further study. Of note, associations of the three SFA groups with metabolic markers may vary by dietary alcohol intake and adiposity. The cross-sectional findings of the present study warrant further investigations in a prospective study. Taken together, these findings help to further our understanding for the role of SFAs in metabolic pathways and disease aetiology.

Abbreviations

ALT:

alanine transaminase

ApoA1:

apolipoprotein A-1

ApoB:

apolipoprotein B

AST:

aspartate transaminase

CRP:

C-reactive protein

EPIC:

European Prospective Investigation into Cancer and Nutrition

GGT:

gamma glutamyl transferase

HbA1c:

haemoglobin A1c

HDL-C:

high-density lipoprotein cholesterol

LDL-C:

low-density lipoprotein cholesterol

SFA:

saturated fatty acid

TC:

total cholesterol

TG:

triglycerides

References

  1. Forouhi NG, Koulman A, Sharp SJ, Imamura F, Kroger J, Schulze MB, et al. Differences in the prospective association between individual plasma phospholipid saturated fatty acids and incident type 2 diabetes: the epic-interact case-cohort study. Lancet Diabetes Endocrinol. 2014;2:810–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Lemaitre RN, Fretts AM, Sitlani CM, Biggs ML, Mukamal K, King IB, et al. Plasma phospholipid very-long-chain saturated fatty acids and incident diabetes in older adults: the cardiovascular health study. Am J Clin Nutr. 2015;101:1047–54.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Ma W, Wu JH, Wang Q, Lemaitre RN, Mukamal KJ, Djousse L, et al. Prospective association of fatty acids in the de novo lipogenesis pathway with risk of type 2 diabetes: the cardiovascular health study. Am J Clin Nutr. 2015;101:153–63.

    Article  CAS  PubMed  Google Scholar 

  4. Santaren ID, Watkins SM, Liese AD, Wagenknecht LE, Rewers MJ, Haffner SM, et al. Serum pentadecanoic acid (15:0), a short-term marker of dairy food intake, is inversely associated with incident type 2 diabetes and its underlying disorders. Am J Clin Nutr. 2014;100:1532–40.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Chowdhury R, Warnakula S, Kunutsor S, Crowe F, Ward HA, Johnson L, et al. Association of dietary, circulating, and supplement fatty acids with coronary risk: a systematic review and meta-analysis. Ann Intern Med. 2014;160:398–406.

    Article  PubMed  Google Scholar 

  6. Matsumoto C, Hanson NQ, Tsai MY, Glynn RJ, Gaziano JM, Djousse L. Plasma phospholipid saturated fatty acids and heart failure risk in the physicians’ health study. Clin Nutr. 2013;32:819–23.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Fretts AM, Mozaffarian D, Siscovick DS, Djousse L, Heckbert SR, King IB, et al. Plasma phospholipid saturated fatty acids and incident atrial fibrillation: the cardiovascular health study. J Am Heart Assoc. 2014;3:e000889.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Kirk EP, Klein S. Pathogenesis and pathophysiology of the cardiometabolic syndrome. J Clin Hypertens (Greenwich). 2009;11:761–5.

    Article  CAS  Google Scholar 

  9. Mayneris-Perxachs J, Guerendiain M, Castellote AI, Estruch R, Covas MI, Fito M, et al. Plasma fatty acid composition, estimated desaturase activities, and their relation with the metabolic syndrome in a population at high risk of cardiovascular disease. Clin Nutr. 2014;33:90–7.

    Article  CAS  PubMed  Google Scholar 

  10. Mahendran Y, Agren J, Uusitupa M, Cederberg H, Vangipurapu J, Stancakova A, et al. Association of erythrocyte membrane fatty acids with changes in glycemia and risk of type 2 diabetes. Am J Clin Nutr. 2014;99:79–85.

    Article  CAS  PubMed  Google Scholar 

  11. Jacobs S, Schiller K, Jansen E, Fritsche A, Weikert C, di Giuseppe R, et al. Association between erythrocyte membrane fatty acids and biomarkers of dyslipidemia in the epic-potsdam study. Eur J Clin Nutr. 2014;68:517–25.

    Article  CAS  PubMed  Google Scholar 

  12. Zong G, Zhu J, Sun L, Ye X, Lu L, Jin Q, et al. Associations of erythrocyte fatty acids in the de novo lipogenesis pathway with risk of metabolic syndrome in a cohort study of middle-aged and older Chinese. Am J Clin Nutr. 2013;98:319–26.

    Article  CAS  PubMed  Google Scholar 

  13. de Oliveira Otto MC, Nettleton JA, Lemaitre RN, Steffen LM, Kromhout D, Rich SS, et al. Biomarkers of dairy fatty acids and risk of cardiovascular disease in the multi-ethnic study of atherosclerosis. J Am Heart Assoc. 2013;2:e000092.

    PubMed  PubMed Central  Google Scholar 

  14. Warensjo E, Jansson JH, Cederholm T, Boman K, Eliasson M, Hallmans G, et al. Biomarkers of milk fat and the risk of myocardial infarction in men and women: A prospective, matched case-control study. Am J Clin Nutr. 2010;92:194–202.

    Article  PubMed  Google Scholar 

  15. Ebbesson SO, Tejero ME, Lopez-Alvarenga JC, Harris WS, Ebbesson LO, Devereux RB, et al. Individual saturated fatty acids are associated with different components of insulin resistance and glucose metabolism: the GOCADAN study. Int J Circumpolar Health. 2010;69:344–51.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Lee E, Lee S, Park Y. N-3 polyunsaturated fatty acids and trans fatty acids in patients with the metabolic syndrome: a case-control study in korea. Br J Nutr. 2008;100:609–14.

    Article  CAS  PubMed  Google Scholar 

  17. Warensjo E, Jansson JH, Berglund L, Boman K, Ahren B, Weinehall L, et al. Estimated intake of milk fat is negatively associated with cardiovascular risk factors and does not increase the risk of a first acute myocardial infarction. A prospective case-control study. Br J Nutr. 2004;91:635–42.

    Article  PubMed  Google Scholar 

  18. Kurotani K, Sato M, Ejima Y, Nanri A, Yi SY, Pham NM, et al. High levels of stearic acid, palmitoleic acid, and dihomo-gamma-linolenic acid and low levels of linoleic acid in serum cholesterol ester are associated with high insulin resistance. Nutr Res. 2012;32:669–75.

    Article  CAS  PubMed  Google Scholar 

  19. Wang L, Folsom AR, Zheng ZJ, Pankow JS, Eckfeldt JH, Investigators AS. Plasma fatty acid composition and incidence of diabetes in middle-aged adults: The atherosclerosis risk in communities (ARIC) study. Am J Clin Nutr. 2003;78:91–8.

    CAS  PubMed  Google Scholar 

  20. Lovejoy JC, Champagne CM, Smith SR, DeLany JP, Bray GA, Lefevre M, et al. Relationship of dietary fat and serum cholesterol ester and phospholipid fatty acids to markers of insulin resistance in men and women with a range of glucose tolerance. Metabolism. 2001;50:86–92.

    Article  CAS  PubMed  Google Scholar 

  21. InterAct C, Langenberg C, Sharp S, Forouhi NG, Franks PW, Schulze MB, et al. Design and cohort description of the interact project: an examination of the interaction of genetic and lifestyle factors on the incidence of type 2 diabetes in the epic study. Diabetologia. 2011;54:2272–82.

    Article  Google Scholar 

  22. Hodson L, Skeaff CM, Fielding BA. Fatty acid composition of adipose tissue and blood in humans and its use as a biomarker of dietary intake. Prog Lipid Res. 2008;47:348–80.

    Article  CAS  PubMed  Google Scholar 

  23. Wang LY, Summerhill K, Rodriguez-Canas C, Mather I, Patel P, Eiden M, et al. Development and validation of a robust automated analysis of plasma phospholipid fatty acids for metabolic phenotyping of large epidemiological studies. Genome Med. 2013;5(4):39.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Friedewald WT, Levy RI, Fredrickson DS. Estimation of the concentration of low-density lipoprotein cholesterol in plasma, without use of the preparative ultracentrifuge. Clin Chem. 1972;18:499–502.

    CAS  PubMed  Google Scholar 

  25. Westgard rules. http://www.westgard.com/mltirule.htm. Accessed 11 Sept 2017.

  26. Riboli E, Hunt KJ, Slimani N, Ferrari P, Norat T, Fahey M, et al. European prospective investigation into cancer and nutrition (EPIC): study populations and data collection. Public Health Nutr. 2002;5:1113–24.

    Article  CAS  PubMed  Google Scholar 

  27. Slimani N, Deharveng G, Unwin I, Southgate DA, Vignat J, Skeie G, et al. The epic nutrient database project (ENDB): a first attempt to standardize nutrient databases across the 10 European countries participating in the epic study. Eur J Clin Nutr. 2007;61:1037–56.

    Article  CAS  PubMed  Google Scholar 

  28. Li J, Ji L. Adjusting multiple testing in multilocus analyses using the eigenvalues of a correlation matrix. Heredity. 2005;95:221–7.

    Article  CAS  PubMed  Google Scholar 

  29. Miller M, Stone NJ, Ballantyne C, Bittner V, Criqui MH, Ginsberg HN, et al. Triglycerides and cardiovascular disease: a scientific statement from the American Heart Association. Circulation. 2011;123:2292–333.

    Article  PubMed  Google Scholar 

  30. Howard BV, Robbins DC, Sievers ML, Lee ET, Rhoades D, Devereux RB, et al. LDL cholesterol as a strong predictor of coronary heart disease in diabetic individuals with insulin resistance and low LDL: the Strong Heart study. Arterioscler Thromb Vasc Biol. 2000;20:830–5.

    Article  CAS  PubMed  Google Scholar 

  31. Gotto Jr AM. High-density lipoprotein cholesterol and triglycerides as therapeutic targets for preventing and treating coronary artery disease. Am Heart J. 2002;144:S33–42.

    Article  CAS  PubMed  Google Scholar 

  32. Yusuf S, Hawken S, Ounpuu S, Dans T, Avezum A, Lanas F, et al. Effect of potentially modifiable risk factors associated with myocardial infarction in 52 countries (the INTERHEART study): Case-control study. Lancet. 2004;364:937–52.

    Article  PubMed  Google Scholar 

  33. McQueen MJ, Hawken S, Wang X, Ounpuu S, Sniderman A, Probstfield J, et al. Lipids, lipoproteins, and apolipoproteins as risk markers of myocardial infarction in 52 countries (the INTERHEART study): A case-control study. Lancet. 2008;372:224–33.

    Article  CAS  PubMed  Google Scholar 

  34. Matsumori R, Miyazaki T, Shimada K, Kume A, Kitamura Y, Oshida K, et al. High levels of very long-chain saturated fatty acid in erythrocytes correlates with atherogenic lipoprotein profiles in subjects with metabolic syndrome. Diabetes Res Clin Pract. 2013;99:12–8.

    Article  CAS  PubMed  Google Scholar 

  35. Grundy SM. Obesity, metabolic syndrome, and cardiovascular disease. J Clin Endocrinol Metab. 2004;89:2595–600.

    Article  CAS  PubMed  Google Scholar 

  36. Jacobs S, Jager S, Jansen E, Peter A, Stefan N, Boeing H, et al. Associations of erythrocyte fatty acids in the de novo lipogenesis pathway with proxies of liver fat accumulation in the EPIC-Potsdam study. PLoS One. 2015;10:e0127368.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Petersson H, Arnlov J, Zethelius B, Riserus U. Serum fatty acid composition and insulin resistance are independently associated with liver fat markers in elderly men. Diabetes Res Clin Pract. 2010;87:379–84.

    Article  CAS  PubMed  Google Scholar 

  38. Khaw KT, Friesen MD, Riboli E, Luben R, Wareham N. Plasma phospholipid fatty acid concentration and incident coronary heart disease in men and women: The EPIC-Norfolk prospective study. PLoS Med. 2012;9:e1001255.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Capuzzi DM, Freeman JS. C-reactive protein and cardiovascular risk in the metabolic syndrome and type 2 diabetes: controversy and challenge. Clinical Diabetes. 2007;25:16–22.

    Article  Google Scholar 

  40. Siler SQ, Neese RA, Hellerstein MK. De novo lipogenesis, lipid kinetics, and whole-body lipid balances in humans after acute alcohol consumption. Am J Clin Nutr. 1999;70:928–36.

    CAS  PubMed  Google Scholar 

  41. Jenkins BJ, Seyssel K, Chiu S, Pan PH, Lin SY, Stanley E, et al. Odd chain fatty acids; new insights of the relationship between the gut microbiota, dietary intake, biosynthesis and glucose intolerance. Sci Rep. 2017;7:44845.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Weitkunat K, Schumann S, Nickel D, Hornemann S, Petzke KJ, Schulze MB, et al. Odd-chain fatty acids as a biomarker for dietary fiber intake: a novel pathway for endogenous production from propionate. Am J Clin Nutr. 2017;105:1544–51.

    CAS  PubMed  Google Scholar 

  43. Schwartz EA, Zhang WY, Karnik SK, Borwege S, Anand VR, Laine PS, et al. Nutrient modification of the innate immune response: a novel mechanism by which saturated fatty acids greatly amplify monocyte inflammation. Arterioscler Thromb Vasc Biol. 2010;30:802–8.

    Article  CAS  PubMed  Google Scholar 

  44. Weigert C, Brodbeck K, Staiger H, Kausch C, Machicao F, Haring HU, et al. Palmitate, but not unsaturated fatty acids, induces the expression of interleukin-6 in human myotubes through proteasome-dependent activation of nuclear factor-kappab. J Biol Chem. 2004;279:23942–52.

    Article  CAS  PubMed  Google Scholar 

  45. Pillon NJ, Azizi PM, Li YE, Liu J, Wang C, Chan KL, et al. Palmitate-induced inflammatory pathways in human adipose microvascular endothelial cells promote monocyte adhesion and impair insulin transcytosis. Am J Physiol Endocrinol Metab. 2015;309:E35–44.

    Article  CAS  PubMed  Google Scholar 

  46. Ameer F, Scandiuzzi L, Hasnain S, Kalbacher H, Zaidi N. De novo lipogenesis in health and disease. Metabolism. 2014;63:895–902.

    Article  CAS  PubMed  Google Scholar 

  47. Grosch S, Schiffmann S, Geisslinger G. Chain length-specific properties of ceramides. Prog Lipid Res. 2012;51:50–62.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

We thank all EPIC participants and staff for their contribution to the study. We thank Nicola Kerrison (MRC Epidemiology Unit, Cambridge) for managing the data for the InterAct Project, Matt Sims for writing the technical laboratory specification for the intermediate pathway biomarker measurements and for overseeing the laboratory work, and Laura Wang and Keith Summerhill (MRC Human Nutrition Research, Cambridge) for the fatty acid measurements. We thank the technical and functional operational teams of the MRC Epidemiology Unit for enabling all biomarker measurements.

Funding

Funding for the InterAct project was provided by the EU FP6 programme (grant number LSHM_CT_2006_037197). In addition, InterAct investigators acknowledge funding from the following agencies: Medical Research Council Epidemiology Unit MC_UU_12015/1 and MC_UU_12015/5, Medical Research Council Human Nutrition Research MC_UP_A090_1006, Cambridge Lipidomics Biomarker Research Initiative (G0800783), and NIHR Biomedical Research Centre Cambridge: Nutrition, Diet, and Lifestyle Research Theme (IS-BRC-1215-20014). IS and YTvdS: Verification of diabetes cases was additionally funded by NL Agency grant IGE05012 and an Incentive Grant from the Board of the UMC Utrecht (The Netherlands). AMWS: Dutch Ministry of Public Health, Welfare and Sports (VWS), Netherlands Cancer Registry (NKR), LK Research Funds, Dutch Prevention Funds, Dutch ZON (Zorg Onderzoek Nederland), World Cancer Research Fund (WCRF), Statistics Netherlands (The Netherlands). PWF: Swedish Research Council, Novo Nordisk, Swedish Heart Lung Foundation, Swedish Diabetes Association. TK: Cancer Research UK C8221/A19170 and C570/A16491 Medical Research Council MR/M012190/1. KO and AT: Danish Cancer Society. JRQ: Asturias Regional Government. RT: AIRE-ONLUS Ragusa, AVIS-Ragusa, Sicilian Regional Government. Additionally, Red Temática de Investigación Cooperativa en Cáncer of the Instituto de Salud Carlos III (ISCIII RTICC RD12/0036/0018), Instituto de Salud Carlos III (PIE14/00045), cofounded by FEDER funds/European Regional Development Fund (ERDF); German Cancer Aid, German Ministry of Research (BMBF); Compagnia di San Paolo; Imperial College Biomedical Research Centre. JSZ is supported by the Cambridge Initiative-Nutrition (RG71466, SJAH/004) and Marie Skłodowska-Curie Fellowships (701708, RG82205, SJAI/051).

Availability of data and materials

The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.

Author information

Authors and Affiliations

Authors

Contributions

JSZ and NGF had full access to all of the data in the study and take responsibility for the integrity of the data and the accuracy of the data analysis and they drafted the manuscript. AK was responsible for the measurement of fatty acids. JSZ performed the statistical analyses. NJW, CL, NGF and SJS coordinated the InterAct project with NJW as chief co-ordinator. The Working Group (JSZ, SJS, FI, AK, MBS, ZY, JMH, MG, JK, IS, NGF) and all authors contributed to interpretation of data, revised the article critically for important intellectual content and approved the final version of the manuscript.

Corresponding authors

Correspondence to Ju-Sheng Zheng or Nita G. Forouhi.

Ethics declarations

All participants gave written informed consent and the study was approved by the local ethics committee in the participating centres and the Internal Review Board of the International Agency for Research on Cancer.

Consent for publication

Not applicable.

Competing interests

In 2015, the Scientific Advisory Committee on Nutrition (SACN), UK established a working group to examine the evidence linking dietary fats and health outcomes; Jules Griffin is a member of that group. The 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 file

Additional file 1: Table S1.

Population characteristics by the eight participating countries: EPIC-InterAct study. Table S2. Distribution of plasma metabolic markers by quartiles of plasma phospholipid saturated fatty acid (SFA) groups: EPIC-InterAct study. Table S3. Pairwise correlation of plasma phospholipid SFAs in the subcohort of EPIC-InterAct study. Table S4. Association of plasma phospholipid odd-chain SFAs (per 1 SD difference) with metabolic markers: EPIC-InterAct study. Table S5. Association of plasma phospholipid even-chain SFAs (per 1 SD difference) with metabolic markers: EPIC-InterAct study. Table S6. Association of plasma phospholipid very-long-chain SFAs (per 1 SD difference) with metabolic markers: EPIC-InterAct study. Table S7. Sensitivity analysis for the association of plasma phospholipid SFA groups and metabolic markers: EPIC-InterAct study. Figure S1. Association of plasma phospholipid SFA groups with LDL-C, TG and GGT by countries: EPIC-InterAct study. Figure S2. Association of plasma phospholipid SFA groups with metabolic markers by subgroups of alcohol intake, BMI and sex: EPIC-InterAct study. (DOCX 321 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

Zheng, JS., Sharp, S.J., Imamura, F. et al. Association between plasma phospholipid saturated fatty acids and metabolic markers of lipid, hepatic, inflammation and glycaemic pathways in eight European countries: a cross-sectional analysis in the EPIC-InterAct study. BMC Med 15, 203 (2017). https://doi.org/10.1186/s12916-017-0968-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12916-017-0968-4

Keywords