Oxylipins are a group of oxidized metabolites of PUFAs metabolized through various enzymatic pathways. These include the more conventional subclass of 20-carbon oxylipins derived from arachidonic acid [AA, 20:4(n-6)], named eicosanoids, and other more recently identified oxylipins that are metabolized via lipoxygenase and cytochrome P450 pathways. These oxylipins have hormone-like effects on inflammation, vasoconstriction, and blood clotting in experimental studies (
8Prostaglandins and leukotrienes: advances in eicosanoid biology.
), and may thus act as mediators of the effects of dietary FAs on cardiovascular health. Changes in dietary FA precursors could change the balance of oxylipins with different biological effects (
9- Shearer G.C.
- Harris W.S.
- Pedersen T.L.
- Newman J.W.
Detection of omega-3 oxylipins in human plasma and response to treatment with omega-3 acid ethyl esters.
,
10- Ferretti A.
- Nelson G.J.
- Schmidt P.C.
- Kelley D.S.
- Bartolini G.
- Flanagan V.P.
Increased dietary arachidonic acid enhances the synthesis of vasoactive eicosanoids in humans.
), which may affect the risk of CHD. However, epidemiological data on oxylipins and CHD risk are sparse. In a study of human carotid atherosclerotic plaques, the concentration of HETEs was significantly higher in symptomatic plaques, as compared with asymptomatic plaques (
11- Mallat Z.
- Nakamura T.
- Ohan J.
- Leseche G.
- Tedgui A.
- Maclouf J.
- Murphy R.C.
The relationship of hydroxyeicosatetraenoic acids and F2-isoprostanes to plaque instability in human carotid atherosclerosis.
). In addition, higher plasma HETE concentrations were associated with CHD in two case-control studies (
12- Shishehbor M.H.
- Zhang R.
- Medina H.
- Brennan M.L.
- Brennan D.M.
- Ellis S.G.
- Topol E.J.
- Hazen S.L.
Systemic elevations of free radical oxidation products of arachidonic acid are associated with angiographic evidence of coronary artery disease.
,
13- Xu Y.J.
- Ho W.E.
- Xu F.
- Wen T.
- Ong C.N.
Exploratory investigation reveals parallel alteration of plasma fatty acids and eicosanoids in coronary artery disease patients.
). However, data from larger prospective studies with a more comprehensive assessment of oxylipins are currently not available. We therefore examined the association between plasma concentrations of FAs and oxylipins and risk of acute myocardial infarction (AMI) in a case-control study nested within the Singapore Chinese Health Study (SCHS) cohort. To circumvent the problem of high collinearity among these metabolites, we have devised an automated approach to group the highly correlated metabolites into clusters and study their associations with AMI risk.
MATERIALS AND METHODS
SCHS
SCHS is a population-based prospective cohort study of 63,257 Chinese men and women (Hokkien or Cantonese dialect group) aged 45–74 years residing in public housing estates, where 86% Singaporeans resided at the time of recruitment. Recruitment and assessment of baseline diet, lifestyle, and medical history took place from 1993 to 1998. To date, follow-up surveys have been conducted twice (1999–2004; 2006–2010) to update information on use of tobacco and alcohol, medical history, and menopausal status. Biospecimens were collected from 28,439 subjects during 1994–2005. The cohort has been followed up for morbidity and mortality through record linkage with Singapore Registry of Births and Deaths, Hospital Discharge Database, and Singapore Cancer Registry. Data on education, BMI, physical activity, smoking, alcohol intake, and history of diabetes and hypertension were assessed using interviewer-administered standardized questionnaires (
14- Hankin J.H.
- Stram D.O.
- Arakawa K.
- Park S.
- Low S.H.
- Lee H.P.
- Yu M.C.
Singapore Chinese Health Study: development, validation, and calibration of the quantitative food frequency questionnaire.
,
15- Koh W.P.
- Yuan J.M.
- Wang R.
- Lee Y.P.
- Lee B.L.
- Yu M.C.
- Ong C.N.
Plasma carotenoids and risk of acute myocardial infarction in the Singapore Chinese Health Study.
). This study was approved by the Institutional Review Boards of National University of Singapore and University of Pittsburgh. All participants gave informed consent.
Nested case-control study
We have established a nested case-control study including 744 incident AMI cases and 744 matched controls. These cases and controls were selected from the SCHS participants who provided blood and did not have a history of physician-diagnosed CHD or stroke reported at blood draw or ascertained through linkage with the Hospital Discharge Database before blood draw (International Classification of Diseases-9 codes: 410–414, 427, 428, 430–434, 438). The eligible cases for the present study were incident nonfatal or fatal AMI that occurred during follow-up from blood collection through December 31, 2010 via linkage with three databases (see supplementary text). Using a risk-set sampling approach (
16- Prentice R.L.
- Breslow N.E.
Retrospective studies and failure time models.
), we randomly chose one control for each case. Control subjects were alive and free of CHD at the time of the AMI diagnosis or death of the index case. The matching criteria included gender, dialect group (Hokkien, Cantonese), date of birth (±5 years), date of recruitment (±2.5 years), and date of blood collection (±6 months).
Measurement of CVD risk factors
During home visits in 1994–2005, systolic and diastolic blood pressure was measured and blood samples were collected in heparin tubes. The blood specimens were separated into various blood components, including plasma, serum, red blood cells, and buffy coat, and were stored in freezers at −80°C for a mean of 11.5 ± 1.8 years before the lab analyses (
17- Koh W.P.
- Yuan J.M.
- Sun C.L.
- van den Berg D.
- Seow A.
- Lee H.P.
- Yu M.C.
Angiotensin I-converting enzyme (ACE) gene polymorphism and breast cancer risk among Chinese women in Singapore.
). Conventional biochemical risk factors of CHD, including total, HDL, and LDL (directly measured) cholesterol, triglycerides, C-reactive protein (CRP), creatinine, and glycosylated hemoglobin (HbA1c), were measured. Detailed measurement methods and coefficients of variation of these cardiovascular risk factors are described in the supplementary text.
Measurement of plasma FAs and oxylipins
Nineteen plasma FAs, covering major saturated FAs (SFAs), monounsaturated FAs, and PUFAs, were measured in a targeted mode using GC-MS/MS on an Agilent 7890 GC system (Shanghai, China) equipped with a G7000B QQQ triple quadrupole mass detector and an auto sample injector. FAs from both free and esterified (triglycerides, phospholipids, cholesterol esters) fractions were measured in total.
Nineteen plasma-free oxylipins with potential roles in inflammation, blood pressure regulation, and platelet degranulation (
supplementary Table 1) were targeted to be measured using reversed-phase LC-MS with an Agilent 1290 ultra-pressure liquid chromatograph (Waldbronn, Germany) coupled to electrospray ionization with iFunnel Technology on a triple quadrupole mass spectrometer. Twelve out of the 19 targeted oxylipins could be reliably quantified in the majority of the samples, and were thus included in the current analysis. These included: seven AA-derived oxylipins [thromboxane (TX)B
2, leukotriene E
4 (LTE
4), 5-HETE, 12-HETE, 15-HETE, 5-ketoeicosatetraenoic acid (5-KETE), 14,15-dihydroxyeicosatrienoic acid (14,15-DHET)]; four linoleic acid [LA, 18:2(n-6)]-derived oxylipins [9-HODE, 13-HODE, 9,10-dihydroxyoctadeca(mono)enoic acid (9,10-DiHOME), 12,13-DiHOME]; and one EPA [20:5(n-3)]-derived oxylipin [5-hydroxyeicosapentaenoic acid (5-HEPE)]. The free fraction of plasma AA, a major precursor of oxylipins, was also assessed in the same assay.
Samples were analyzed in 76 batches for FAs and 51 batches for oxylipins. Cases and their matched controls were included in the same batch. All procedures were conducted with lab personnel being blinded to the case control status of the samples. Pooled human plasma was used for quality control (QC) samples that were included at regular intervals (one QC in every 10 samples). There were three and four QCs in each batch for FA and oxylipin analysis, respectively. The experimental details and the coefficients of variation of the measured FAs and oxylipins are listed in the supplementary text and
supplementary Table 2.
Statistical analysis
Our nested case-control study initially included a total of 744 matched case-control pairs. Due to lack of plasma samples, we could not obtain reliable oxylipin values in 17 matched pairs of cases and controls. After exclusion of these participants, 727 matched case-control pairs remained for the current analysis.
Prior to any statistical analysis, we first performed normalization of the quantitative data in three steps. First, we adjusted for batch-effect by dividing the molar concentration of each individual FA and oxylipin by the average concentration of their corresponding compound in the QC samples for each batch. Second, the concentrations for FAs and oxylipins were converted to proportions of their respective total concentrations. For oxylipins, this approach of using proportions of total concentrations is less established, so sensitivity analysis was conducted on the raw concentrations as well. In the last step, the proportions (or the raw concentrations in the sensitivity analysis) were log-transformed, centered by mean, and scaled to have a standard Gaussian distribution for each compound.
We then devised a procedure to automatically group highly correlated compounds together to avoid multi-collinearity problems. We performed hierarchical clustering of the Pearson correlation matrix (agglomerative clustering with Euclidean distance) for the FAs and the oxylipins, respectively, using the data of the controls only. We determined the predictor groups using the Dynamic Tree Cut procedure (
18- Langfelder P.
- Zhang B.
- Horvath S.
Defining clusters from a hierarchical cluster tree: the Dynamic Tree Cut package for R.
) implemented in R (
19., R Core Team. 2013. R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing, Vienna, Austria.
), which adaptively selected a sparse set of clusters given a dendrogram (clustering tree) obtained by initial clustering results. Once the compounds were grouped into clusters, each cluster was subjected to further inspection, where it was divided into sub-clusters if the average pairwise correlation coefficient between the compounds was less than 0.3 within that cluster. In our data, the Dynamic Tree Cut procedure reported 12 initial clusters, and the second-stage inspection resulted in 17 clusters: 12 for FAs and 5 for oxylipins. Given the final cluster configuration, a centroid concentration pattern for each cluster was computed by the average of all its member compounds, i.e.:
where
xj denotes the batch-adjusted normalized log-concentration of compound
j belonging to cluster
i,
c denotes the total number of clusters, and
ni denotes the total number of compounds present in cluster
i.
Demographic characteristics, established cardiovascular risk factors, and plasma concentrations of FAs and oxylipins in the cases and controls were compared using univariate conditional logistic regression.
The association between these 17 generated clusters and the risk of AMI was assessed using conditional logistic regression models by computing the multivariable-adjusted odds ratios (ORs) and 95% confidence intervals (CIs). In the multivariable model, in addition to control for the five matching factors as a result of the matching design and use of conditional logistic regression analysis, we further adjusted for age (years), hours of fasting before blood collection (<2 h; 2 h to <8 h; ≥8 h), education level (none; primary school; secondary school or above), cigarette smoking (never; ex-smoker; current smoker with <13 cigarettes/day; current smoker with ≥13 cigarettes per day), alcohol consumption (grams per day), physical activity (0 h per week; <4 h per week of moderate and <2 h per week of strenuous activity; ≥4 h per week of moderate or ≥2 h per week of strenuous activity), BMI (kilograms per square meter), history of hypertension (yes/no), and history of diabetes (yes/no). The regression models were compared with or without adjustment for other FA clusters, oxylipin clusters, and conventional cardiovascular risk factors.
For specific FAs that showed significant attenuation after adjustment for the conventional cardiovascular risk factors, cross-sectional association between the FAs and the cardiovascular risk factors was assessed by multivariable linear regression among the control subjects.
The initial data transformation, normalization, and cluster generation analyses were carried out using R version 3.0.3. The rest of the statistical analyses were conducted using Stata version 11 (StataCorp, College Station, TX). All P values reported are two-sided. P < 0.05 is considered to be statistically significant.
DISCUSSION
In this prospective nested case-control study of Singapore Chinese, we observed an inverse association between long-chain n-3 PUFAs and stearic acid with risk of AMI, and a positive association for AA after adjustment for other FAs. We also observed an inverse association between TXB2, an AA-derived oxylipin, and risk of AMI. The association between plasma FAs and AMI was partly mediated by established cardiovascular risk factors, but not by the oxylipins measured.
The association between long-chain n-3 PUFAs and CHD has been extensively examined in previous cohort studies, most of which also showed an inverse association (
20Fish intake, contaminants, and human health: evaluating the risks and the benefits.
). In contrast, plasma stearic acid had a positive or null association with CHD risk in previous studies in Western populations (
21- Chowdhury R.
- Warnakula S.
- Kunutsor S.
- Crowe F.
- Ward H.A.
- Johnson L.
- Franco O.H.
- Butterworth A.S.
- Forouhi N.G.
- Thompson S.G.
- et al.
Association of dietary, circulating, and supplement fatty acids with coronary risk: a systematic review and meta-analysis.
). In a Japanese study, the only other study conducted in Asia, plasma stearic acid was inversely associated with CHD risk (
22- Itakura H.
- Yokoyama M.
- Matsuzaki M.
- Saito Y.
- Origasa H.
- Ishikawa Y.
- Oikawa S.
- Sasaki J.
- Hishida H.
- Kita T.
- et al.
Relationships between plasma fatty acid composition and coronary artery disease.
), which was consistent with our findings. The disparity in results for Western versus Asian populations may reflect differences in dietary sources of stearic acid. Stearic acid is highly correlated with other SFAs in Western diet because they share the same predominant food sources, such as red meat and dairy (
23- Hu F.B.
- Stampfer M.J.
- Manson J.E.
- Ascherio A.
- Colditz G.A.
- Speizer F.E.
- Hennekens C.H.
- Willett W.C.
Dietary saturated fats and their food sources in relation to the risk of coronary heart disease in women.
). It is therefore difficult to distinguish the effect of stearic acid from other SFAs in a Western diet. In contrast, the major food sources of SFAs in Singapore are plant-based cooking oils and coconut oil/milk (
24., Singapore Health Promotion Board 2010. Report of the National Nutrition Survey.
). The correlation between plasma stearic acid and other even-chain SFAs in our data (
r range 0.14–0.20) was weaker than that observed in Western populations, which facilitated examination of the specific association for stearic acid. A lack of detrimental effects of stearic acid is supported by feeding trials which showed that dietary stearic acid did not result in elevated serum cholesterol concentrations or impaired insulin sensitivity, as compared with carbohydrates or monounsaturated FAs (
25Effects of stearic acid on plasma lipid and lipoproteins in humans.
,
26- Louheranta A.M.
- Turpeinen A.K.
- Schwab U.S.
- Vidgren H.M.
- Parviainen M.T.
- Uusitupa M.I.
A high-stearic acid diet does not impair glucose tolerance and insulin sensitivity in healthy women.
).
Adipose tissue AA was associated with higher risk of AMI in case-control studies in Costa Rica and Israel (
27Arachidonic acid in adipose tissue is associated with nonfatal acute myocardial infarction in the central valley of Costa Rica.
,
28- Kark J.D.
- Kaufmann N.A.
- Binka F.
- Goldberger N.
- Berry E.M.
Adipose tissue n-6 fatty acids and acute myocardial infarction in a population consuming a diet high in polyunsaturated fatty acids.
), which is consistent with our observation for plasma AA. In most other previous studies, no association or an inverse association between AA biomarkers and CHD was observed (
21- Chowdhury R.
- Warnakula S.
- Kunutsor S.
- Crowe F.
- Ward H.A.
- Johnson L.
- Franco O.H.
- Butterworth A.S.
- Forouhi N.G.
- Thompson S.G.
- et al.
Association of dietary, circulating, and supplement fatty acids with coronary risk: a systematic review and meta-analysis.
,
29- de Goede J.
- Verschuren W.M.
- Boer J.M.
- Verberne L.D.
- Kromhout D.
- Geleijnse J.M.
N-6 and N-3 fatty acid cholesteryl esters in relation to fatal CHD in a Dutch adult population: a nested case-control study and meta-analysis.
), and few studies adjusted for other FAs. One interesting observation in our analysis is that AA became significantly associated with a higher AMI risk only after adjustment for other FA clusters. Confounding by other FAs can arise if they share similar food sources or if their food sources are part of the same dietary pattern. It is therefore important to consider such confounding by other FAs in studies of circulating FAs and health outcomes.
There have been concerns about LA intake, as it may be converted to AA endogenously, and AA may increase levels of detrimental oxylipins which increase blood pressure and inflammation (
6- James M.J.
- Gibson R.A.
- Cleland L.G.
Dietary polyunsaturated fatty acids and inflammatory mediator production.
,
30Essential fatty acids in health and chronic disease.
). However, our data suggest that plasma LA is not associated with AMI risk, and LA is not necessarily associated with higher AA concentrations (
r = −0.06). In line with this finding, wide variation of dietary LA intake did not substantially affect tissue AA concentrations in a meta-analysis of feeding trials (
31Increasing dietary linoleic acid does not increase tissue arachidonic acid content in adults consuming Western-type diets: a systematic review.
). Therefore, it is unlikely that dietary LA increases risk of AMI through conversion to AA. Furthermore, we could not identify potential detrimental effects of AA through increased production of oxylipins. The cardiovascular risk factors that AA was associated with were also not the risk factors plausibly affected by oxylipins. Besides the commonly known pro-inflammatory and pro-aggregatory metabolites, such as the leukotrienes and TXs, AA is also the precursor to other anti-inflammatory or vasodilating metabolites, such as the lipoxins (
32Lipoxins and aspirin-triggered 15-epi-lipoxins are the first lipid mediators of endogenous anti-inflammation and resolution.
) and epoxyeicosatrienoic acids (
33- Node K.
- Huo Y.
- Ruan X.
- Yang B.
- Spiecker M.
- Ley K.
- Zeldin D.C.
- Liao J.K.
Anti-inflammatory properties of cytochrome P450 epoxygenase-derived eicosanoids.
). It is therefore difficult to predict the net impact of the entire AA metabolome, likely depending on the differential enzymatic activity for the synthesis of metabolites, as well as the interactions among the metabolites (
34Omega-6 fatty acids and cardiovascular disease: friend, not foe?.
). Hence, detailed mechanistic studies are warranted to elucidate potential underlying mechanisms for the direct association between AA and AMI risk.
The pro-thrombotic oxylipin, TXA
2, functions as a transient autocrine, and its biological inactive metabolite, TXB
2, has been used to indirectly measure TXA
2 production (
35- Catella F.
- Healy D.
- Lawson J.A.
- FitzGerald G.A.
11-Dehydrothromboxane B2: a quantitative index of thromboxane A2 formation in the human circulation.
). The inverse association between TXB
2 is thus unexpected and requires further study. We did not observe substantial association between other oxylipin clusters and AMI risk besides TXB
2. The autocrine and short-lived nature of these oxylipins might account for the lack of association in our long-term prospective study (
36Recent eicosanoid chemistry.
). However, we cannot exclude the possibility that other oxylipins or oxylipins within the esterified pool of lipoprotein particles, rather than the free oxylipins that we examined, predict CHD. In addition, the concentrations of TXB
2 and possibly other oxylipins may reflect sample collection, processing, and storage conditions rather than only biological differences in vivo, and the association between oxylipins and AMI in our study should therefore be interpreted with caution.
To our knowledge, our study is the first prospective study of a wider range of oxylipins in relation to AMI risk. In addition, it is among the largest studies of FA biomarkers and AMI in Asia. The values we obtained for FAs and oxylipins were within the range of those reported in earlier studies, suggesting a reliable quantification. Furthermore, the clustering among various FAs and oxylipin families was as expected based on biological knowledge. However, our study also has several potential limitations. First, there may be misclassifications in the measurement of plasma FAs and oxylipins. For example, we only included measurements of FAs and oxylipins at a single point in time, which may not fully capture long-term effects on CHD development. Plasma FAs reflect dietary FA intake of the past 1–2 weeks (
37- Skeaff C.M.
- Hodson L.
- McKenzie J.E.
Dietary-induced changes in fatty acid composition of human plasma, platelet, and erythrocyte lipids follow a similar time course.
). However, misclassification of FA and oxylipin status was likely to be nondifferential and probably led to weaker associations with AMI risk. Second, we measured only major oxylipins in plasma rather than the full spectrum. Other oxylipins that we did not examine in this study might play a more significant role in mediating certain FAs. Third, the stability of the analytes may deteriorate with long storage time, which may affect the detected levels in the studied samples. The controls were thus individually matched to the cases for time of blood collection to minimize the influence of storage duration. Fourth, as an observational study, our results may have been affected by residual confounding by imperfectly or unmeasured metabolic, genetic, or behavioral risk factors. Finally, our results may have been affected by multiple testing, as we examined 17 clusters as exposure variables. Using Bonferroni correction (α = 0.05/17 = 0.003), the association for long-chain n-3 PUFAs remains statistically significant, whereas associations for stearic acid, AA, and the cluster including TXB
2 were only nominally significant. This highlights the need to confirm these findings in an independent study population.
In summary, in this study, we measured a wide spectrum of FAs and oxylipins, and clustered the compounds based on their correlation structure to reduce the multi-collinearity in the regression analysis. The resulting clusters were consistent with shared dietary sources and endogenous metabolic pathways. The plasma concentration of the oxylipin, TXB2, was inversely associated with risk of AMI. This association was unexpected, as TXB2 is a marker for the pro-thrombotic TXA2, and warrants further investigation in other prospective studies. Plasma long-chain n-3 PUFAs and stearic acid were associated with a lower risk of AMI and plasma total AA was associated with a higher risk of AMI. The nontypical findings for stearic acid and AA suggest that it is informative to examine FAs in non-Western populations with different dietary habits and that it is important to adjust for other FAs in the data analysis.
Article info
Publication history
Published online: May 24, 2016
Received in revised form:
April 27,
2016
Received:
January 9,
2016
Footnotes
This study was supported by Office of Extramural Research, National Institutes of Health Grants R01 CA144034 and UM1 CA182876, and National Medical Research Council, Singapore (NMRC/1270/2010). Additional support was provided by U.S. Department of Agriculture Intramural Project (23032-51530-022-00D to J.W.N). The US Department of Agriculture is an equal opportunity provider and employer. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.
Abbreviations:
AAarachidonic acid
AMIacute myocardial infarction
CHDcoronary heart disease
CIconfidence interval
CRPC-reactive protein
DHETdihydroxyeicosatrienoic acid
DiHOMEdihydroxyoctadeca(mono)enoic acid
HbA1cglycosylated hemoglobin
HEPEhydroxyeicosapentaenoic acid
KETEketoeicosatetraenoic acid
LAlinoleic acid
LTE4leukotriene E4
ORodds ratio
QCquality control
SCHSSingapore Chinese Health Study
SFAsaturated FA
TXthromboxane
Copyright
© 2016 ASBMB. Currently published by Elsevier Inc; originally published by American Society for Biochemistry and Molecular Biology.