|
|
||||||||
ORIGINAL RESEARCH COMMUNICATION |
1 From the Department of Clinical and Experimental Medicine (NM, DG, PG, MS, SF, FP, RC, and OO) and the Section of Biology and Genetics, Department of Mother and Child and Biology–Genetics (GM, ET, and PFP), University of Verona, Italy; and the GSF-National Research Center for Environment and Health, Institute of Epidemiology, Neuherberg, Germany (TI, LS, and JH)
2 Supported by grant no. 2005/065152 from the Italian Ministry of University and Research and grants from the Veneto Region and the Cariverona Foundation, Verona, Italy.
3 Address reprint requests to O Olivieri, Department of Clinical and Experimental Medicine, University of Verona, Policlinico GB Rossi, Piazzale L A Scuro 10, 37134 Verona, Italy. E-mail: oliviero.olivieri{at}univr.it.
| ABSTRACT |
|---|
|
|
|---|
-5 and
-6 desaturases, encoded by FADS1 and FADS2 genes, are key enzymes in polyunsaturated fatty acid (PUFA) metabolism that catalyze the conversion of linoleic acid (LA) into arachidonic acid (AA) and that of
-linolenic acid (ALA) into eicosapentaenoic acid (EPA). Single-nucleotide polymorphisms (SNPs) in FADS1 and FADS2 have been associated with different concentrations of AA and LA, and those associations have possible functional consequences for desaturase activity. Objective: We aimed to evaluate the possible association among FADS genotypes, desaturase activity, inflammation, and coronary artery disease (CAD).
Design: Thirteen FADS SNPs and the ratio of AA to LA (AA/LA) on red blood cell (RBC) membranes, a marker of desaturase activity, were evaluated in 876 subjects with (n = 610) or without (n = 266) angiographically documented CAD.
Results: Both AA/LA and the ratio of EPA to ALA (EPA/ALA) were higher in patients with CAD than in those without CAD, but, in a multiple logistic regression model, only a higher AA/LA resulted an independent risk factor for CAD (odds ratio: 2.55; 95% CI: 1.61, 4.05 for higher compared with lower ratio tertile; P for trend < 0.001). Furthermore, concentrations of high-sensitivity C-reactive protein increased progressively across tertiles of AA/LA. Graded increases in high-sensitivity C-reactive protein concentrations and CAD risk were related to the carriership of FADS haplotypes, including the alleles associated with a higher ratio.
Conclusion: In populations following a Western diet, subjects carrying FADS haplotypes that are associated with higher desaturase activity may be prone to a proinflammatory response favoring atherosclerotic vascular damage.
| INTRODUCTION |
|---|
|
|
|---|
The fatty acid (FA) composition of the cell membrane is influenced by both dietary intake and metabolic pathways. For example, AA originates from both diet and the elongation-desaturation process of its precursor, linoleic acid (LA; 18:2n–6). The
-5 (D5D) and
-6 (D6D) desaturases are key enzymes of this pathway (Figure 1
): D6D catalyzes the conversion from LA to
-linolenic acid (18:3n–6), which is elongated to dihomo-
-linolenic acid (20:3n–6), which is in turn desaturated to AA by D5D. D5D and D6D are also involved in the n–3 FA pathway (Figure 1
), which favors the conversion of
-linolenic acid (ALA) into eicosapentaenoic acid (EPA). D5D and D6D are encoded by FADS1 and FADS2 genes, respectively. The FADS gene cluster is located on chromosome 11 (11q12–13.1) and includes a third gene, FADS3, that shares 52–62% sequence identity with the FADS1 and FADS2 genes and encodes for an as-yet-unidentified protein (6). Recently, several polymorphisms and haplotypes of the FADS1-FADS2 gene cluster have been shown to be associated with serum phospholipid PUFA concentrations in German and Italian populations (7, 8). Various FADS alleles were associated with high concentrations of AA and low concentrations of LA, an association that provided evidence for a functional effect on desaturase activity. However, these genetic variants may merely be in linkage with as-yet-unidentified causal variants.
|
| SUBJECTS AND METHODS |
|---|
|
|
|---|
1 significant stenosis (eg,
50% lumen reduction). According to the hypothesis to be tested, subjects with nonadvanced CAD (ie, coronary stenosis < 50%) were not included in the study. The angiograms were assessed by cardiologists who were unaware that the patients were to be included in the present study. All participants came from the same geographical area (Northern Italy) and were of similar socioeconomic status. At the time of blood sampling, a complete clinical history, including the assessment of cardiovascular risk factors such as obesity, smoking, hypertension, and diabetes, was collected.
Written informed consent was obtained from all participants after a full explanation of the study. The study was approved by the Ethics Committee of the Azienda Ospedaliera (Verona, Italy).
Biochemical analysis
Samples of venous blood were drawn from each subject after an overnight fast. Serum lipids and other CAD risk factors, including hs-CRP, were measured as described previously (14).
Measurement of serum phospholipids and red blood cell membrane fatty acid concentrations
Blood samples were transported to the laboratory within 1 h of collection and processed immediately. Analysis of plasma (100 µL) and erythrocyte membrane (250 µL packed erythrocytes hemolyzed in an equal volume of double-distilled water) FAs was performed on total lipids extracted by using 4.5 mL of a ratio of isopropanol to chloroform (11:7, by vol) containing 0.45 mmol 2,6-di-ter-p-cresol/L as antioxidant. A gas-chromatographic method (Hewlett-Packard 5980 chromatograph, equipped with an HP-FFAP 25-m, 0.2-mm internal diameter, 0.3-µm phase column; Hewlett-Packard, Palo Alto, CA) based on a direct FA transesterification technique was used as described previously (13). Analyses were performed in duplicate on each sample. Peak identification and quantification were done with commercially available reference FAs (Sigma, St Louis, MO). We used 17:0 as the internal standard. The areas of the peaks were measured, and the quantification was performed with the aid of a Vectra QS/16S PC equipped with HP-3365 CHEM STATION software (version 3.0; Hewlett-Packard). FAs were expressed as g/100 g FA methyl ester (% by wt). FAs from 12:0 to 26:0 were measured [saturated FAs: 12:0, 14:0, 16:0, 18:0, 20:0, 22:0, 24:0, and 26:0; monounsaturated FAs: 16:1, 18:1n–9, and 20:1; and PUFAs: 18:2n–6, 18:3n–3, 20:2n–6, 20:4n–6, 20:5n–3, and 22:6n–3); unidentified peaks accounted for <0.5% of the total. Only the most important FAs for this analysis (ie, the PUFAs involved along the desaturase pathway) were specified in the results. All CVs were <5%.
Single-nucleotide polymorphisms and genotyping
Thirteen SNPs tagging the greatest variability of the FADS gene region were selected from the European Community Respiratory Health Survey (7, 8) and the HapMap project (Internet: www.hapmap.org): rs174545, rs174556, rs174561, rs3834458, rs174570, rs2524299, rs174583, rs174589, rs498793, rs174611, rs17831757, rs174627, and rs1000778. For further information, see Figure S1 under "Supplemental data" in the current online issue.
Genotyping was performed as described previously (7). Briefly, genotyping was performed by using matrix-assisted laser desorption–ionization time-of-flight–mass spectrometry (MALDI-TOF-MS) to detect allele-specific primer extension products (Mass Array; Sequenom, San Diego, CA). Polymerase chain reaction primers were designed by using the Mass Array–Assay Design program (Sequenom).
Statistical analysis
Calculations were performed with SPSS software (version 13.0; SPSS Inc, Chicago, IL). Distributions of continuous variables in groups were expressed as means ± SDs. Skewed variables, such as hs-CRP, were logarithmically transformed, and the geometric means (and 95% CIs) were reported. Quantitative data were assessed by using Student's t test or analysis of variance (ANOVA), with polynomial contrast for linear trend and Tukey's post hoc comparison of the means when indicated. Correlations between quantitative variables were assessed by using Pearson's correlation test. Qualitative data were analyzed with the chi-square test and with chi-square for linear trend analysis when indicated. P < 0.05 was considered statistically significant.
Within each group examined, the frequencies of the FADS genotypes associated with each of the polymorphisms were compared by using the chi-square test, and the values were predicted on the basis of the Hardy-Weinberg equilibrium. Pairwise linkage disequilibrium was examined as described by Devlin and Risch (15).
To assess possible genetic predictors of AA/LA variability, FADS polymorphisms were included in a linear regression model according to a forward-stepwise variables selection. Such analysis was replicated by means of a different, nonlinear statistical approach based on the estimation of variable importance by using random forest plots (16, 17). Four SNPs were independent contributors to AA/LA variability and were included into the final model. After excluding significant interactions among SNPs, we focused on an additive model by calculating a score reflecting the sum of risk alleles (ie, alleles associated with a greater AA/LA) for each patient.
Haplotype frequencies were estimated by using R software with the haplo.stats package (R Foundation for Statistical Computing, Vienna, Austria; ISBN 3–900051-07–0; Internet: http://www.R-project.org). Haplotypes present in <10 subjects were excluded from the analyses. The relations between haplotypes and laboratory and clinical outcomes were examined by using generalized linear model regression of a trait on haplotype effects, with adjustment for ambiguous haplotypes (haplo.glm function) (18). The statistical significance of associations was ascertained by random permutation of the disease status of the sample subjects in 1000 replicates with the use of the Monte Carlo method.
To assess associations with CAD, odds ratios (ORs) (and 95% CIs) were estimated by logistic regression models adjusted for both classic risk factors (eg, sex, age, smoking, hypertension, diabetes, body mass index, creatinine, cholesterol, triglycerides, and hs-CRP) and FAs differently distributed between subjects with and without CAD. Statistical power was calculated by using the Altmann nomogram.
| RESULTS |
|---|
|
|
|---|
|
Substantially analogous results were obtained by comparing CAD patients and CAD-free subjects for FA composition of serum phospholipids (See Table S1 under "Supplemental data" in the current online issue.). In particular, serum AA/LA was higher in CAD patients than in CAD-free subjects (P < 0.001, t test).
According to the tertile distribution of RBC-AA/RBC-LA in the whole population, the proportion of CAD patients increased significantly (P < 0.001, chi-square test for trend analysis) from the lowest to the highest tertile (Figure 2A
). It is noteworthy that hs-CRP concentrations also increased progressively across the ratio tertiles (P < 0.001, ANOVA with polynomial contrast for linear trend; Figure 2B
). Considering the subjects with a ratio in the lowest tertile as the reference group, a progressively greater risk of CAD was observed in the other groups after adjustment for all the classic risk factors—ie, sex, age, smoking, hypertension, diabetes, body mass index, creatinine, cholesterol, and triglycerides (Figure 2C
). Such associations with CAD remained significant after the inclusion in the logistic model of FAs differently distributed between CAD patients and CAD-free subjects—ie, RBC-saturated FAs, RBC-ALA, and RBC-DHA. When both n–3 and n–6 ratios (AA/LA and EPA/ALA) were included in the regression model, only RBC-AA/RBC-LA remained significantly associated with hs-CRP and CAD (data not shown).
|
All SNPS in the FADS1-FADS2-FADS3 cluster were associated with the concentrations of RBC-LA or RBC-AA or with RBC-AA/RBC-LA (See Table S3 under "Supplemental data" in the current online issue). Some SNPs were also associated with the concentration of RBC-ALA or with RBC-EPA/RBC-ALA (See Table S4 under "Supplemental data" in the current online issue), but no single variant was differently distributed between the CAD and CAD-free groups (See Table S5 under "Supplemental data" in the current online issue).
An additive model and haplotype analysis
When we looked for significant associations between FADS polymorphisms and RBC-AA/RBC-LA, all SNPs were included in a linear regression model with RBC-AA/RBC-LA as the dependent variable. With the use of a stepwise selection method, 4 significant associations were found: RS174545 G
C (β-coefficient = –0.30), RS174570 C
T (β-coefficient = –0.07), RS174583 C
T (β-coefficient = 0.14), and RS1000778 C
T (β-coefficient = 0.04) with G, C, T, and T alleles, respectively, associated with a high ratio. These results were also confirmed by means of a different statistical approach based on the estimation of variable importance by using random forest plots. It is noteworthy that none of these 4 SNPs were synonymous, because each was located in a different LD block (see Figure S1 under "Supplemental data in the current online issue).
Based on the number of alleles associated with a higher RBC-AA/RBC-LA, an additive model predicted AA/LA variability. Given that very few subjects carried 2 (n = 2) or 7 (n = 3) alleles, those subjects were jointly considered with those who carried 3 or 6 alleles, respectively. As expected, subjects with a higher number of such alleles had a higher RBC-AA/RBC-LA (2.14 ± 0.43 for carriers of 6–7 alleles and 1.94 ± 0.33 for carriers of 2–3 alleles; P < 0.05, ANOVA with Tukey post hoc comparison). As is consistent with the previously observed association of RBC-AA/RBC-LA with hs-CRP concentrations and CAD risk, subjects carrying more risk alleles also had progressively higher concentrations of hs-CRP (P < 0.05, ANOVA with polynomial contrast for linear trend; Figure 3A
); they also were more frequently represented in the CAD group (P for trend = 0.006, chi-square test; Figure 3B
). Subjects with a higher number of risk alleles (ie, 6 or 7) have a greater CAD risk than do those with 2 or 3 alleles (OR: 2.75; 95% CI: 1.27, 5.98). This association remained significant after adjustment for all of the classic risk factors—sex, age, smoking, hypertension, diabetes, body mass index, creatinine, cholesterol, and triglycerides (OR: 2.68; 95% CI: 1.08, 6.66)—and also for RBC-LA and RBC-ALA concentrations (OR: 2.69; 95% CI: 1.02, 7.13).
|
|
|
| DISCUSSION |
|---|
|
|
|---|
The second significant result is that a continuous gradient of both higher hs-CRP and a greater risk of CAD was observed across AA/LA tertiles in the whole study population (Figure 2B and C
). Thus, it seems reasonable to suppose that inappropriately high elongase or desaturase activity (as reflected by the higher ratio of 20:0 FA to precursor) may indicate a peculiar susceptibility to the inflammatory stimuli involving the arterial wall during the atherosclerotic process. Whereas the idea of a relation between a sensitive marker of subclinical inflammation such as hs-CRP and CAD is long-lasting (20-22), a link of both of those variables with n–6 AA/LA has been unrecognized to date, and it represents an element of novelty in the present work.
EPA/ALA was also higher in CAD patients than in CAD-free subjects; in addition, EPA/ALA and AA/LA were significantly correlated. Theoretically, higher ratios of n–3 metabolite to precursor are associated with higher EPA or DHA concentrations, and they have protective rather than harmful effects. Thus, the central point is that a greater activity of the same enzymatic complex should have both pro-inflammatory (as a result of greater bioavailability of AA) and anti-inflammatory (as a result of an increased bioavailability of long-chain n–3 FAs) consequences. This apparent inconsistency may be explained by considering n–3 and n–6 PUFAs from a quantitative point of view. In the present study and in many other studies, LA was almost 100 times as abundant as ALA, and AA was almost 30 times as abundant as EPA. The prevalence of n–6 FAs in Western diets may account for the proinflammatory effect we observed, and the detrimental effects of increased desaturase activity may be exaggerated by diets abundant in n–6 and relatively deficient in n–3 FAs (23). It is consistent with this point of view that, when both EPA/ALA and AA/LA were included in a regression model with hs-CRP or CAD as dependent variables, only AA/LA remained as a significant predictor.
A third major result, consistent with the previous ones, is the finding that haplotypes of the FADS gene cluster, including variants associated with an elevated AA/LA, were also related to both a higher hs-CRP concentration and greater risk of CAD (Figure 3
and Tables 2
and 3
). As a consequence, apparently distinct phenotypic features (elevated AA/LA, increase in hs-CRP, and greater CAD risk) linked together functionally and with specific variants on the FADS gene complex, which in turn affected desaturase activity. Given the relatively small effect of any single gene variant and the observed functional effect in terms of CAD risk, the involved mechanism seems to be of considerable relevance. It may be hypothesized as follows: 1) carriers of specific FADS haplotypes are constitutively able to form more AA from dietary LA; 2) greater AA availability at the membrane and cellular levels facilitates the synthesis of AA-derived inflammatory mediators such as leukotrienes and prostaglandins; 3) vascular inflammatory responses involved in the atherosclerotic process are greatly amplified by the cascade (with multiple effects) of AA-derived mediators; and 4) carriers of unfavorable haplotypes are at a greater risk of developing clinically significant CAD. Although this hypothesis implies an oversimplification of the complex eicosanoid metabolism—eg, inhibition of this pathway by Cox-2 inhibitors may promote CAD risk by reducing the concentration of the protective prostacyclin—the consistent association among AA/LA, hs-CRP concentration, and CAD risk was impressive.
In the population of the present study, patients in the highest tertile of AA/LA had concentrations of hs-CRP nearly double those in patients in the lowest tertile. This difference implies that a difference of
20% in desaturase activity (the cutoff was 1.90 and 2.26 for the lowest and highest tertile, respectively) accounted for a 100% increase in the concentration of the inflammatory marker. Given that AA is the substrate of 5-lipoxygenase and cycloxygenase, key enzymes in the biosynthetic pathway of leukotrienes and prostaglandins, it is reasonable that a modest increase in the substrate availability produces larger effects on the subsequent cascade of mediators. Accordingly, selective inhibitors of D5D and D6D have shown effective antiinflammatory properties in mice (24).
We have measured a suitable inflammatory marker of cardiovascular risk, hs-CRP, because it has been clinically validated in large, prospective population studies (20-22). However, it is plausible that a similar association may be observed with other molecules that are more subtly linked to the vascular atherosclerotic damage, such as interleukins and chemokines.
In the literature, few and conflicting results on desaturase activity and atherosclerotic risk have been reported so far. Some reports suggested a possible causality link between vascular disease and lower desaturase activity (25). Other studies did not shown any significant difference in AA concentration between CAD patients and CAD-free subjects (26, 27). Several years ago, our group showed that patients with essential hypertension have a higher RBC-AA/LA than do normotensive persons, and patients with "white coat" hypertension (that is not confirmed by ambulatory 24-h blood pressure monitoring) have intermediate values (28). Similarly, a higher ratio was observed in persons with insulin resistance and metabolic syndrome (29-32), conditions in which hs-CRP often is abnormally elevated. Thus, those studies support the view that exaggerated desaturase activity may characterize patients carrying well-recognized, major risk factors for CAD. The present study is the first to show an association between desaturase activity and angiographically confirmed CAD.
In the population of the present study, CAD patients had lower LA concentrations than did CAD-free subjects. Given that LA is an essential FA derived exclusively from the dietary source, a deficiency of LA intake may also have contributed to the vascular disease of the patients in the present study. Several studies have reported an inverse relation between LA concentration and risk of atherosclerosis (5, 26, 33-35), and LA deficiency seems to promote mitochondrial respiratory uncoupling with resulting excessive oxidative stress (36).
To the best of our knowledge, this is the first study showing an association between FADS genotypes and CAD. In the only other study reporting on the issue, Baylin et al (37) showed that the rs3834458 FADS2 promoter polymorphism was associated with lower EPA concentrations in adipose tissue and higher triglyceride concentrations in serum, but not with CAD risk, in a population from Costa Rica. Those authors hypothesized that impaired desaturase activity may induce an unbalanced proportion of n–3 and n–6 20:0 FAs, which would have both favorable and unfavorable effects of lower AA and EPA, respectively; this possibility underscores the potentially Janus-faced (n–3- or n–6-related) effects of desaturase activity. The present results show that the proinflammatory, detrimental effects are predominant in a population assumed to be following a Western diet that is abundant in n–6 FAs and somewhat deficient in n–3 FAs (23).
In the present study, FADS gene variants had little or no effect on CAD risk when considered individually, and this result was consistent with that reported by Baylin et al (37). However, a more complex approach that took into consideration an additive model of FADS SNPs, FADS haplotypes, or both showed a significant effect on a person's susceptibility to CAD. FADS polymorphisms are relatively frequent in the general population, so that only a mild effect on a multifaceted phenotype such as CAD could theoretically be expected. Nevertheless, an additive model of unfavorable alleles was shown to fit well in explaining an increase in the inflammatory marker hs-CRP and a greater probability of CAD. This complex genotype x phenotype interaction appears biologically plausible after adjustment for the potential effects of FADS gene variants on desaturase activity and AA availability. Schaeffer et al (7) showed that FADS1 and FADS2 gene polymorphisms and haplotypes are correlated with PUFA membrane content. As is consistent with the role of AA concentration as an important determinant of the inflammatory response, the prevalence of inflammatory disorders such as allergic rhinitis and atopic eczema seems to be lower in carriers of FADS gene variants associated with lower AA concentrations (7).
The present study has some limitations, the first of which is the retrospective case-control design. Furthermore, a somewhat lower number of controls than of CAD patients was enrolled. However, in consideration of the genotype frequency previously reported for FADS gene polymorphisms in the general population and a genotype-associated near-doubling of CAD risk, the statistical power of the analysis was >95% by Altmann nomogram at an 0.05 significance level after adjustment for unequal sample size. Another limitation of the present study was related to the lack of data for and di-homo-
-linolenic acid FAs—intermediate FAs of the metabolic conversion of linoleic or linolenic acids by D5D and D6D—which could provide a more specific estimation of both of those enzymatic activities. In contrast, a relative strength of the study is represented by the angiographic evaluation of coronary arteries in both cases and controls, which allowed for a clear-cut definition of the clinical phenotype, a crucial requirement in any genotype-phenotype association study.
In conclusion, the present study is the first to show an association between greater desaturase activity, FADS genotypes, an increase in inflammatory markers, and a greater risk of CAD. Our data support a role of PUFA metabolism in CVD and suggest that a greater desaturase activity, in the context of a Western diet characterized by an unbalanced ratio of n–6 to n–3, may contribute to the atherosclerotic process by favoring a proinflammatory vascular environment. Clearly, these findings have to be confirmed by larger prospective studies; however, they seem to indicate that therapeutic strategies aimed at reducing cell AA availability (eg, by n–3 PUFA supplementation) would be particularly advisable in patients with a genetically determined increase in desaturase activity.
| ACKNOWLEDGMENTS |
|---|
The authors' responsibilities were as follows—NM and OO; study design, statistical analysis, and writing of the manuscript; DG: study design, data collection, and review of the manuscript; GM, TI, ET, LS, JH, and PFP: conduct of the molecular biology portion of the study, data analysis, and review of the manuscript; PG: fatty acids analysis, data analysis, and review of the manuscript; MS: statistical analysis; FP and SF: patient selection and recruitment, data collection, and review of the manuscript; and RC: coordination of the conduct of the study and review of the manuscript. All authors participated in the study and take public responsibility for the content of this report. None of the authors had a personal or financial conflict of interest.
| REFERENCES |
|---|
|
|
|---|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |