#PAGE_PARAMS# #ADS_HEAD_SCRIPTS# #MICRODATA#

Genetic Loci Associated with Plasma Phospholipid n-3 Fatty Acids: A Meta-Analysis of Genome-Wide Association Studies from the CHARGE Consortium


Long-chain n-3 polyunsaturated fatty acids (PUFAs) can derive from diet or from α-linolenic acid (ALA) by elongation and desaturation. We investigated the association of common genetic variation with plasma phospholipid levels of the four major n-3 PUFAs by performing genome-wide association studies in five population-based cohorts comprising 8,866 subjects of European ancestry. Minor alleles of SNPs in FADS1 and FADS2 (desaturases) were associated with higher levels of ALA (p = 3×10−64) and lower levels of eicosapentaenoic acid (EPA, p = 5×10−58) and docosapentaenoic acid (DPA, p = 4×10−154). Minor alleles of SNPs in ELOVL2 (elongase) were associated with higher EPA (p = 2×10−12) and DPA (p = 1×10−43) and lower docosahexaenoic acid (DHA, p = 1×10−15). In addition to genes in the n-3 pathway, we identified a novel association of DPA with several SNPs in GCKR (glucokinase regulator, p = 1×10−8). We observed a weaker association between ALA and EPA among carriers of the minor allele of a representative SNP in FADS2 (rs1535), suggesting a lower rate of ALA-to-EPA conversion in these subjects. In samples of African, Chinese, and Hispanic ancestry, associations of n-3 PUFAs were similar with a representative SNP in FADS1 but less consistent with a representative SNP in ELOVL2. Our findings show that common variation in n-3 metabolic pathway genes and in GCKR influences plasma phospholipid levels of n-3 PUFAs in populations of European ancestry and, for FADS1, in other ancestries.


Published in the journal: Genetic Loci Associated with Plasma Phospholipid n-3 Fatty Acids: A Meta-Analysis of Genome-Wide Association Studies from the CHARGE Consortium. PLoS Genet 7(7): e32767. doi:10.1371/journal.pgen.1002193
Category: Research Article
doi: https://doi.org/10.1371/journal.pgen.1002193

Summary

Long-chain n-3 polyunsaturated fatty acids (PUFAs) can derive from diet or from α-linolenic acid (ALA) by elongation and desaturation. We investigated the association of common genetic variation with plasma phospholipid levels of the four major n-3 PUFAs by performing genome-wide association studies in five population-based cohorts comprising 8,866 subjects of European ancestry. Minor alleles of SNPs in FADS1 and FADS2 (desaturases) were associated with higher levels of ALA (p = 3×10−64) and lower levels of eicosapentaenoic acid (EPA, p = 5×10−58) and docosapentaenoic acid (DPA, p = 4×10−154). Minor alleles of SNPs in ELOVL2 (elongase) were associated with higher EPA (p = 2×10−12) and DPA (p = 1×10−43) and lower docosahexaenoic acid (DHA, p = 1×10−15). In addition to genes in the n-3 pathway, we identified a novel association of DPA with several SNPs in GCKR (glucokinase regulator, p = 1×10−8). We observed a weaker association between ALA and EPA among carriers of the minor allele of a representative SNP in FADS2 (rs1535), suggesting a lower rate of ALA-to-EPA conversion in these subjects. In samples of African, Chinese, and Hispanic ancestry, associations of n-3 PUFAs were similar with a representative SNP in FADS1 but less consistent with a representative SNP in ELOVL2. Our findings show that common variation in n-3 metabolic pathway genes and in GCKR influences plasma phospholipid levels of n-3 PUFAs in populations of European ancestry and, for FADS1, in other ancestries.

Introduction

High levels of n-3 polyunsaturated fatty acids (PUFA) in plasma phospholipids, cell membranes, and whole blood have been associated with lower risk of multiple diseases, including sudden cardiac death [1], fatal coronary heart disease [2], non-fatal myocardial infarction [3], heart failure [4], thickening of the carotid arteries [5], metabolic syndrome [6], breast cancer [7], chronic obstructive pulmonary disease [8], depression [9], [10] and dementia [11], [12].

N-3 PUFAs are derived directly from the diet, including the plant-derived essential fatty acid α-linolenic acid (ALA, 18:3n3) and the seafood-derived long-chain n-3 PUFAs eicosapentaenoic acid (EPA, 20:5n3) and docosahexaenoic acid (DHA, 22:6n3) [13], [14]. Long-chain n-3 PUFAs can also be produced from ALA by the series of desaturation and elongation steps in the pathway shown in Figure 1; docosapentaenoic acid (DPA, 22:5n-3) can be produced from EPA. The pathway enzymes may be a major source of circulating long-chain n-3 PUFAs in people who consume very little or no seafood. However, the conversion of ALA to EPA and DHA has been shown to be generally low [15][17], and it is not known whether common genetic variation in the pathway affects this conversion. The n-6 essential fatty acid linoleic acid (LA) is elongated to long-chain n-6 PUFAs by the same pathway enzymes and could also compete with the conversion of ALA to EPA; it is not known whether genetic variation affects such competition.

Fig. 1. N-3 polyunsaturated fatty acid metabolic pathway and summary of genome-wide associations in pathway.
N-3 polyunsaturated fatty acid metabolic pathway and summary of genome-wide associations in pathway.
The fatty acids indicated in bold were examined in this study. The genome-wide significant associations of two loci with each fatty acid are shown with dashed arrows. + and − signs indicate the direction of the associations.

There is evidence of co-heritability of EPA and DHA levels in erythrocyte membrane phospholipids [18]. Investigation of genetic factors influencing PUFA levels has largely focused on candidate genes, such as the desaturase genes FADS1 and FADS2, among participants of European ancestry [19][25]. Only one prior study reported a genome-wide association of n-3 PUFA levels evaluating total plasma n-3 PUFAs which includes triacylglycerols, phospholipids and free fatty acids, among 1075 participants [26]. The study found an association of EPA with variants in the FADS1 gene that reached genome-wide significance level; independent follow-up investigation showed associations of a selected FADS1 variant with erythrocyte membrane levels of EPA, ALA and DPA and of an ELOVL2 variant with DPA and DHA. These findings confirm an influence of FADS1 and ELOVL2 on selected n-3 PUFAs. However, statistical power may have been limited to confirm an influence of these genes on all four major n-3 PUFAs, of other genes in these pathways (e.g., FADS2), or of additional genes in other unknown biologic pathways. Prior studies have also not had adequate power to evaluate potential interaction between variation in these genes, EPA and DHA, and: (a) diet, (b) ALA levels, and (c) LA levels. In addition, there is limited information on genetic variation and n-3 PUFA levels in subjects of non-European ancestry.

To understand how common genetic variation affects n-3 PUFA phospholipid levels and potentially uncover novel associations, we conducted a meta-analysis of pre-planned genome-wide association analyses of plasma phospholipid n-3 PUFAs in 8,866 participants of European ancestry in five population-based studies, as part of the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) Consortium [27]. We evaluated the four main n-3 PUFAs of the metabolic pathway, (ALA, EPA, DPA, and DHA) separately and accounted for the intercorrelations between these fatty acids. In addition, we investigated whether consumption of fatty fish or phospholipid levels of ALA and LA influenced the association of the identified genetic markers with EPA and DHA levels. Finally, we studied the most highly associated SNPs from the meta-analyses among samples of European ancestry in additional samples from African, Chinese and Hispanic ancestry.

Results

Meta-analysis of genome-wide associations of n-3 fatty acids

The study samples for the genome-wide association study (GWAS) comprised a total of 8,866 subjects of European ancestry. Table 1 shows sample size, demographic characteristics, fish consumption and phospholipid n-3 PUFA levels in the 5 cohorts that contributed to the study. Participants ranged from 21 to 102 years of age. Across the cohorts, mean levels of ALA varied from 0.14% to 0.44% of total fatty acids; EPA from 0.56% to 1.01%; DPA from 0.83% to 0.98%; and DHA from 2.29% to 5.09%. Relatively higher ALA levels were seen in the InCHIANTI cohort, likely reflecting differences in the composition of total plasma fatty acids (InCHIANTI) versus phospholipid fatty acids (ARIC, CARDIA, CHS, and MESA).[28]

Tab. 1. CHARGE cohorts descriptives.*
CHARGE cohorts descriptives.<em class=&quot;ref&quot;>*</em>
Descriptive characteristics of the cohorts of European ancestry included in the GWAS analyses and the cohorts of African, Chinese, and Hispanic ancestry included in the association analyses of selected SNPs.

Figure 2A–2D show the meta-analysis of the genome-wide association results for ALA, EPA, DPA, and DHA. Variation in one or both of two major genetic loci was associated with plasma phospholipid levels of each n-3 PUFA at genome-wide levels of significance. The two loci are illustrated in association plots for EPA (Figure 3). One locus, on chromosome 11q12.2, contained the C11orf9/10, FEN1 and the desaturase genes FADS1, FADS2 and FADS3. The other locus, on chromosome 6p24.2, contained SYPC2L and the elongase gene ELOVL2. Many highly correlated SNPs reached genome-wide significance in the associations with n-3 PUFAs (Table 2 and Table S1AS1D). Variant alleles at SNPs in the chromosome 11 locus were associated with higher levels of ALA and lower levels of EPA and DPA, and variant alleles at SNPs in the chromosome 6 locus were associated with higher levels of EPA and DPA and lower levels of DHA (Table 2, Table 3, Table 4). From the meta-analysis results, we estimated that the most highly associated SNPs on chromosome 11 explained 3.8% of the variance of ALA, 2.0% of the variance of EPA and 8.6% of the variance of DPA. The most highly associated SNPs on chromosome 6 explained 0.4% of the variance of EPA, 2.8% of the variance in DPA and 0.7% of the variance in DHA.

Fig. 2. Meta-analysis of genome-wide associations with n-3 polyunsaturated fatty acids.
Meta-analysis of genome-wide associations with n-3 polyunsaturated fatty acids.
A: α-linolenic acid (ALA), B: eicosapentaenoic acid (EPA), C: docosapentaenoic acid (DPA), D: docosahexaenoic acid (DHA). Associations were graphed by chromosome position and –log10 (p-value) up to p-values of 10−10. Triangles indicate additional SNPs with p-values less than 10−10. Genes of interest in the regions with significantly associated SNP variants are indicated.

Fig. 3. SNP association plots for EPA–associated regions.
SNP association plots for EPA–associated regions.
Genome-wide association significance level is plotted against the y-axis as –log10 (p-value). Genetic coordinates are as per NCI build 36. A) FADS cluster region. LD is indicated by color scale in relationship to marker rs174538. B) ELOVL2 region. LD is indicated by color scale in relationship to marker rs3798713. The color scheme is red for strong linkage disequilibrium (LD; r2≥0.8) and fading color for lower LD.

Tab. 2. Loci associated with SNP markers with P-values <5.0×10−8.
Loci associated with SNP markers with P-values &lt;5.0×10<sup>−8</sup>.
Results of the meta-analyses of n-3 fatty acid GWAS; main results with standard adjustments described in Methods, and additional results further adjusted for one other fatty acid or for selected SNPs.

Tab. 3. Associations of rs174548 (G allele) with n-3 polyunsaturated fatty acids in samples of different ancestries.
Associations of rs174548 (G allele) with n-3 polyunsaturated fatty acids in samples of different ancestries.
For each fatty acid, results are shown for the individual cohorts and for the meta-analyses among samples of European ancestry and samples of African ancestry.

Tab. 4. Associations of rs3734398 (C allele) with n-3 polyunsaturated fatty acids in samples of different ancestries.
Associations of rs3734398 (C allele) with n-3 polyunsaturated fatty acids in samples of different ancestries.
For each fatty acid, results are shown for the individual cohorts and for the meta-analyses among samples of European ancestry and samples of African ancestry.

Another genome-wide significant association with DPA was observed with SNPs on chromosome 2 in the GCKR gene (most associated SNP: rs780094, p = 9.0×10−9, Table 2 and Figure 4A). In addition, DPA showed a possible association with SNPs in AGPAT3, a gene on chromosome 21 involved in phospholipid metabolism (most associated SNP: rs7453, p = 2.4×10−7; Figure 4B and Table S1C).

Fig. 4. SNP association plots for DPA–associated regions.
SNP association plots for DPA–associated regions.
Genome-wide association significance level is plotted against the y-axis as –log10 (p-value). Genetic coordinates are as per NCBI build 36. A) GCKR region. LD is indicated by color scale in relationship to marker rs780094. B) AGPAT3 region. LD is indicated by color scale in relationship to marker rs7435. The color scheme is red for strong linkage disequilibrium (LD; r2≥0.8) and fading color for lower LD.

Levels of EPA and DHA were correlated, as were levels of ALA and EPA. Spearman correlations in the 5 cohorts ranged from 0.26 to 0.57 for EPA-DHA and from 0.20 to 0.28 for ALA-EPA. When two outcomes are positively correlated but exhibit genetic associations in opposite directions, it is possible to increase the power of discovery efforts by adjusting the association between SNPs and one fatty acid trait for the other [29]. For this reason, we performed meta-analyses of genome-wide association results for ALA adjusted for EPA, EPA adjusted for ALA, EPA adjusted for DHA and DHA adjusted for EPA. These analyses did not reveal additional genome-wide significant loci, although the statistical significance of the adjusted associations was increased (Table 2). For example, the significance corresponding to the association of the SNP rs2236212 in ELOVL2 with DHA increased from p = 1.3×10−15 to p = 7.0×10−39 with adjustment for EPA, with many more genome-wide significant associations in that region. In addition, the association of rs4985167 on chromosome 16 (PDXDC1) with ALA approached genome-wide significance with adjustment for EPA (p = 7.6×10−8).

To explore whether the signals in the chromosome 6 and 11 loci could be explained by one SNP alone, we performed GWAS of DPA with adjustment for rs2236212 and rs174547 in the ARIC, CARDIA, CHS and MESA cohorts and meta-analyzed the results. No new association was found in the chromosome 11 locus where the adjustment reduced the association for all other SNPs, from a minimum p-value of 1.1×10−51 to p>10−6 (not shown). In contrast, 22 SNPs in the ELOVL2 region reached genome-wide significance, with ten of these SNPs previously undetected (Table 2, Table S1E). The A allele of rs12662634, the most highly associated SNP, had minor allele frequency of 0.18 and was associated with lower level of DPA (regression coefficient associated with one copy of A allele: −0.030, p = 2.7×10−10). This association was new, as the estimated effect in the unadjusted GWAS of DPA was 0.008 (p = 0.06). Among the 12 SNPs that had been detected before, rs9368564 showed the most significant association (regression coefficient and p-value associated with one copy of G allele: 0.027, p = 1.4×10−9, in the adjusted GWAS, and 0.048, p = 10−40, in the unadjusted GWAS). This result may reflect residual association in the adjusted analysis or yet another new association.

To explore if reducing the variability in EPA would reveal additional associations, we also performed GWAS of EPA adjusted for estimated fish intake. No additional associations beyond those previously seen for SNP in FADS1/2 and ELOVL2 were observed in these analyses (not shown). Finally, adjustment of the GWAS of the four n-3 PUFAs for levels of triglycerides, high density lipoprotein and low density lipoprotein did not materially change the results (not shown). For example, in the GWAS of DPA, the estimated effect of one copy of the T allele of rs780094 (GCKR) went from 0.0157 (p = 1.15×10−8) without adjustment to 0.0175 (p = 2.52×10−9) with adjustment, in meta-analyses that included 7663 participants.

Association of top SNPs with n-3 PUFAs in samples from participants of African, Chinese, and Hispanic ancestry

To determine whether the gene-fatty acid associations were consistent across different ethnicities, we examined the associations of genotype at two representative SNPs with phospholipid n-3 PUFA levels, in samples of African, Chinese and Hispanic ancestry. Results of these analyses for the selected SNPs are shown in Table 3 and Table 4, together with meta-analysis and cohort-specific results among the samples of European ancestry. Frequency of the G allele of rs174548 (FADS1) was 0.29, 0.21, 0.58, 0.52 and frequency of the C allele of rs3734398 (ELOVL2) was 0.43, 0.25, 0.92, 0.57 in samples of European, African, Chinese and Hispanic ancestry respectively. Associations of rs174548 with n-3 PUFA were generally similar across all ancestries, with the G allele associated with higher ALA and lower long-chain n-3 PUFA levels, although associations did not always reach statistical significance, perhaps due to limited sample sizes (Table 3). The associations of rs3734398 with EPA, DPA and DHA were similar for samples of African ancestry versus European ancestry (Table 4). Among samples of Chinese ancestry, SNP rs3734398 was not highly polymorphic (C allele frequency of 92%) and no significant associations were detected. In samples of Hispanic ancestry, the C allele of rs3734398 was associated with higher DPA and lower DHA, but it was not associated with EPA.

Interactions

We evaluated several potential interactions in the samples of European ancestry, with statistical significance defined at alpha = 0.004 (0.05 divided by 13 tested interactions). We found little evidence that fatty fish consumption (≥vs. <0.6 servings/week) modified the associations of rs1535 (FADS2) or rs3734398 (ELOVL2) with levels of DHA or EPA. We also did not observe any interaction between plasma phospholipid levels of LA (continuous linear) and genotype at these two SNPs on the levels of DHA or EPA. Plasma phospholipid levels of ALA (continuous linear) also did not modify the association of genotype at these two SNPs with levels of DHA, or of genotype at rs3734398 with levels of EPA. However, there was a significant interaction of ALA with rs1535 genotype and EPA levels (meta-analyzed interaction coefficient p = 9.3×10−7), illustrated in Figure 5. Per one SD unit (0.05% of total fatty acids) increase in ALA, EPA levels increased by 0.086% of total fatty acids (23% of one SD) in the absence of the minor allele (G); by 0.063% (17% of one SD) in the presence of one copy (G-); and by 0.036% (10% of one SD) in the presence of two copies (GG).

Fig. 5. Influence of variation in rs1535 on the estimated association of plasma phospholipid ALA levels with EPA levels.
Influence of variation in rs1535 on the estimated association of plasma phospholipid ALA levels with EPA levels.
Panel A illustrates the association of ALA with EPA in carriers of 0, 1, and 2 copies of rs1535 G allele with fitted lines constructed from meta-analyzed estimates. Panel B shows the % change in EPA corresponding to a 0.05 unit increase in ALA in carriers of 0 (top estimate), 1 (middle estimate), 2 (bottom estimate) copies of rs1535 G allele.

Discussion

We report here the results of the largest GWAS of plasma phospholipid n-3 PUFAs to date, with 8,866 participants of European ancestry. The associations of the two top hits on chromosomes 6 and 11 are shown in context of the n-3 PUFA pathway in Figure 1. Genetic variation in the desaturase genes FADS1 and FADS2 was associated with higher levels of ALA and lower levels of EPA and DPA suggesting genetic variants that affect the conversion of ALA to EPA and DPA. In the main analyses, genetic variation in the elongase gene ELOVL2 was associated with higher levels of EPA and DPA and lower levels of DHA, suggesting variants that decrease the conversion of EPA and DPA to DHA. These reciprocal associations support a role of genetic variation in the pathway for circulating levels of n-3 PUFAs in free-living populations.

The associations of FADS1/2 and ELOVL2 with n-3 PUFAs were generally consistent with the previous GWAS in InCHIANTI (total plasma n-3 PUFAs) and follow-up candidate replication in the Genetics of Lipid-Lowering Drugs and Diet Network (GOLDN) Study (erythrocyte n-3 PUFAs) [26]. In that prior study, only the association of variation in FADS1 with EPA reached genome-wide significance. Follow-up investigations in GOLDN suggested associations of FADS1 gene variation with ALA and DHA, and of ELOVL2 gene with DHA and DPA, such as the ones reported here [26]. In the present larger meta-analyses, these prior suggestive associations reached genome-wide significance, providing a definite picture of the association of genes in the PUFA pathway with phospholipid n-3 PUFAs. While the InCHIANTI cohort showed similar estimates of association with ALA, EPA and DHA as the other cohorts in the present study, further studies are needed to further explore genetic associations of n-3 PUFAs in triglycerides and other fatty acid compartments.

Using ratios of n-6 PUFAs as proxy of desaturase activities, Bokor et al reported that the minor allele of FADS1 rs174546 was associated with lower delta-5 desaturase activity and the minor allele of FADS2 rs968567 was associated with higher delta-6 desaturase activity [30]. We found both rs174546 and rs968567 were associated with higher levels of ALA and lower levels of EPA and DPA at genome-wide significance level (Table S1AS1C). Furthermore, adjustment of the GWAS of DPA for the most associated SNP from the FADS1/2 genes and the most associated SNP from ELOVL2 did not reveal additional associations on chromosome 11. This result suggests that each associated SNP conveyed the information contained in the other SNPs of the same broad region on chromosome 11. In contrast, we found evidence for two independent associations, exemplified by rs2236212 and rs12662634, in the region of the ELOVL2 gene.

The associations of the FADS1/2 and ELOVL2 genes with the phospholipid levels of EPA and DHA did not vary depending on the frequency of fatty fish consumption, suggesting the genetic effects are independent of fish intake at the levels of consumption in the studied populations. The absence of interaction of FADS1/2 by fish consumption is consistent with findings from the Netherlands KOALA Birth Cohort Study [31], in which higher levels of fish consumption were associated with similar slopes of plasma phospholipid EPA and DHA levels among individuals with 0, 1, or 2 copies of minor FADS1/2 alleles. The associations of FADS1/2 and ELOVL2 genes with EPA and DHA also did not vary depending on phospholipid LA. Linoleic acid is desaturated and elongated to n-6 metabolites (e.g. arachidonic acid) using the same desaturases and elongase(s) as ALA, and in dietary trials, higher dietary LA reduces plasma phospholipid EPA [32]. Animal experiments suggest that high dietary LA inhibits FADS2 gene expression [33]. Our findings do not support an influence of LA on the association of genetic variation in the pathway with n-3 PUFA levels, at levels of LA consumption in these cohorts.

We report for the first time a GWAS of DPA, a central intermediate in the n-3 fatty acid pathway (Figure 1). While present in small quantities in fatty fish, DPA plasma levels appear unrelated to dietary intake [3], suggesting a primarily metabolic origin. Supporting this, we found that DPA exhibited stronger genetic associations than the other n-3 PUFAs.

In addition to its association with variants of desaturase and elongase genes, DPA was associated with variation in the glucokinase regulator gene GCKR, a pleiotropic gene associated with multiple outcomes in GWAS [34]. The T allele of rs780094 is associated with lower fasting glucose and insulin [35] and with higher triglycerides [36][38]; this allele was associated with higher DPA levels in the present study, and the association was independent of triglyceride levels. Given the known influence of long-chain n-3 PUFAs on hepatic triglyceride production [39] and possibly glucose-insulin homeostasis [40], the mechanism of potential pleiotropic effects of this allele on both DPA and these pathways merit further attention.

We found a potential association of DPA with AGPAT3, a gene encoding 1-acylglycerol-3-phosphate O-acyltransferase 3. DPA is a known substrate for the AGPAT3 protein, which transfers a fatty acid in sn-2 position of lysophosphatic acid, a step in the phospholipid biosynthesis pathway. A possible association of DPA with AGPAT3 variation supports an origin of phospholipid DPA from de novo phospholipid synthesis. In contrast, phospholipid EPA and DHA often originate to a greater extent from diet and are predominantly integrated into phospholipids by the process of acyl-chain remodelling [41]. The genetic associations reported here together with growing evidence of the association of DPA with lower risk of coronary heart disease [3], [42], [43] should stimulate further work on factors regulating this fatty acid.

The GWAS of ALA adjusted with EPA revealed a possible new association of phospholipid ALA with variation in PDXDC1. The PDXDC1 protein, a vitamin B6-dependent decarboxylase, is expressed preferentially in the intestine [44], but its function is not known. Animal studies support an influence of dietary vitamin B6 (pyridoxal) on serum and liver levels of ALA and other PUFAs [45], [46], which has been interpreted as an effect on desaturase enzymes activity. The association of PDXDC1 with ALA, if confirmed in other studies, raises the possibility of another vitamin B6-dependent protein affecting ALA, for example through involvement in intestinal ALA absorption.

In addition to an overall association of FADS1/2 variation with less ALA and more EPA, we found that the minor G allele of rs1535 was associated with a reduction of the magnitude of the association between ALA and EPA. In persons with two copies of the G allele, the association of ALA with EPA was less than half the association observed in persons with two copies of the A allele. These results suggest an influence of variation in FADS1/2 on the rate of conversion of ALA into EPA. This conversion is of great clinical and public health interest, given the evidence for importance of long-chain n-3 PUFAs (such as EPA) in many chronic diseases, their limited dietary supply worldwide, and the much greater potential supply of plant-derived ALA. On average, the conversion of ALA to EPA is quite low [47]. Prior tracer studies in humans have shown that the majority of dietary ALA is beta-oxidized for energy or directed into long-term storage as triglycerides, with less than 5% being incorporated into phospholipids where ALA is more readily converted to EPA [16], [17], [47]. Genetic variation that increases or decreases the rate of conversion of ALA to EPA could have implications for individual-based recommendations for consumption of plant- versus seafood-derived n-3 PUFA. The genetic variation may also indicate novel targets for drugs that may increase this conversion.

The associations of a representative SNP in the FADS1/2 genes observed in the meta-analyses of samples of European ancestry were generally similar in samples of African, Chinese, and Hispanic ancestry. Associations of ELOVL2 were less consistent in different ancestries. However, the frequency of the ELOVL2 rs3734398 G allele varied substantially with ancestry, from 25% in African samples to 92% in Chinese samples. Lack of association may be due to inadequate statistical power, chance, different background diet [14], or possible race/ethnic differences in the activity of elongases from the ELOVL2 and ELOVL5 genes.

Our study, the largest GWAS to-date of fatty acid biomarkers, demonstrates key associations of genetic variation with phospholipid n-3 PUFA levels, including genes in the n-3 PUFA metabolic pathway and, for DPA, novel pathways including the pleiotropic gene GCKR. Our results also imply that common variation may result in less efficient conversion of ALA to EPA.

Materials and Methods

Ethics statement

All cohort participants gave written informed consent, including consent to participate in genetic studies. All studies received approval from local ethical oversight committees.

Study samples

The data were obtained from 2 cohort studies in the CHARGE Consortium, the Atherosclerosis Risk in Communities (ARIC) Study and the Cardiovascular Health Study (CHS), and 3 additional cohort studies, the Coronary Artery Risk Development in Young Adults (CARDIA) Study, the Invecchiare in Chianti (InCHIANTI) Study, and the Multi-Ethnic Study of Atherosclerosis (MESA).

Fatty acid measurements

In all cohorts but InCHIANTI, plasma phospholipids were first isolated by thin layer chromatography; fatty acids were then separated by gas chromatography. In InCHIANTI, total plasma fatty acids were measured using a similar gas chromatography technique. Details of fatty acid measurements are provided in Text S1. Levels of EPA, DHA, ALA and DPA were expressed as % of total fatty acids.

Imputation and statistical analysis

Genotyping was done in each cohort separately using high-density SNP marker platforms (ARIC, CARDIA and MESA - Affymetrix 6.0, CHS - Illumina 370, InCHIANTI - Illumina 550). Samples with call rates below 95% (ARIC, CARDIA, MESA), or 97% (CHS, InCHIANTI) at genotyped markers were excluded. Genotypes were imputed to approximately 2.5 million HapMap SNPs by using either MACH [48] (ARIC, InCHIANTI), BIMBAM [49] (CHS), BEAGLE [50] (CARDIA) or IMPUTE [51] (MESA). SNPs for which testing Hardy Weinberg equilibrium resulted in p<10−5 (CHS) or p<10−6 (ARIC) were excluded from imputation. SNPs with minor allele frequency (MAF) ≤1% were excluded from the meta-analyses. Additional details on genotyping and imputation per cohort are provided in Text S1.

Association analysis between genotype and each fatty acid was done separately within each study cohort according to a pre-specified plan. All studies conducted linear regression analysis using an additive genetic model, i.e. regression of phenotype on the number of reference alleles, or equivalently the imputed dosage for imputed genotypes. All analyses were adjusted for age, sex, and site of recruitment where appropriate, and used robust standard errors. In addition, CARDIA, CHS and MESA analyses were adjusted for principal components to account for possible population genetic substructure. The results in InCHIANTI included in the present study have been previously published [26].

Meta-analysis of main effects

For each SNP and fatty acid, GWAS-specific results were combined using inverse-variance weighted meta-analysis in METAL (www.sph.umich.edu/csg/abecasis/metal). Genomic control correction was applied to each study prior to the meta-analysis. Genomic control correction factors ranged from 1.00 to 1.07 (ALA), 1.00–1.08 (EPA), 1.00–1.03 (DPA) and 1.01–1.13 (DHA). P-values less than 5×10−8 were considered significant. Because total plasma levels of ALA (measured in InCHIANTI) are higher than plasma phospholipid levels of ALA (measured in the other cohorts), we performed a z-score based meta-analysis of ALA with the 5 cohorts as a sensitivity analysis. Results did not differ from that of inverse-variance weighted meta-analysis, i.e. from those presented. The proportion of fatty acid variance explained by a particular variant allele was calculated from the formula.

2*2*MAF(1-MAF))/Var(Y), where β is the regression coefficient for one copy of the allele, MAF is the minor allele frequency and Var(Y) is the variance of the fatty acid.

Interaction analyses

We tested 13 interactions using cross-products in the linear regression models. Two of the most associated SNPs available in all cohorts (rs1535 in FADS2 and rs3734398 in ELOVL2) were chosen for investigation of interactions with a) fatty fish intake (dichotomized at 0.6 servings/week, a cut-point around the 25th percentile of fish consumption in the CHS and ARIC cohorts), b) plasma phospholipid ALA levels (continuous linear) and c) plasma phospholipid LA levels (continuous linear) on the outcomes of EPA and DHA. Additionally, we tested the interaction between rs1535 and plasma phospholipid EPA with DHA levels as the outcome. Interaction coefficients from cohort-specific analyses were meta-analyzed. For interactions of SNPs with ALA on the outcomes of EPA and DHA, we performed z-score meta-analysis with all the cohorts to assess statistical significance, and inverse-variance weighted meta-analysis excluding InCHIANTI to estimate the magnitude of the interaction. P-values less than 0.004 (0.05/13 tests) were considered significant for the interactions.

Analyses of selected SNPs in cohorts of African, Chinese, and Hispanic ancestry

We investigated the association of two selected SNPs which had been directly genotyped as part of candidate gene studies in the African American cohort in CARDIA and the African, Chinese and Hispanic American cohort in MESA, and which were available from genome-wide scans on African Americans in the CHS cohort. We used linear regression and additive models as described above. Results in the 3 African American cohorts were meta-analyzed using inverse-variance weighted meta-analysis.

Supporting Information

Attachment 1

Attachment 2


Zdroje

1. SiscovickDSRaghunathanTEKingIWeinmannSWicklundKG 1995 Dietary intake and cell membrane levels of long-chain n-3 polyunsaturated fatty acids and the risk of primary cardiac arrest. Jama 274 1363 1367

2. LemaitreRNKingIBMozaffarianDSotoodehniaNReaTD 2006 Plasma phospholipid trans fatty acids, fatal ischemic heart disease, and sudden cardiac death in older adults: the cardiovascular health study. Circulation 114 209 215

3. SunQMaJCamposHRexrodeKMAlbertCM 2008 Blood concentrations of individual long-chain n-3 fatty acids and risk of nonfatal myocardial infarction. Am J Clin Nutr 88 216 223

4. YamagishiKNettletonJAFolsomAR 2008 Plasma fatty acid composition and incident heart failure in middle-aged adults: the Atherosclerosis Risk in Communities (ARIC) Study. Am Heart J 156 965 974

5. Sala-VilaACofanMPerez-HerasANunezIGilabertR 2010 Fatty acids in serum phospholipids and carotid intima-media thickness in Spanish subjects with primary dyslipidemia. Am J Clin Nutr 92 186 193

6. WarensjoESundstromJLindLVessbyB 2006 Factor analysis of fatty acids in serum lipids as a measure of dietary fat quality in relation to the metabolic syndrome in men. Am J Clin Nutr 84 442 448

7. ShannonJKingIBMoshofskyRLampeJWGaoDL 2007 Erythrocyte fatty acids and breast cancer risk: a case-control study in Shanghai, China. Am J Clin Nutr 85 1090 1097

8. ShaharEBolandLLFolsomARTockmanMSMcGovernPG 1999 Docosahexaenoic acid and smoking-related chronic obstructive pulmonary disease. The Atherosclerosis Risk in Communities Study Investigators. Am J Respir Crit Care Med 159 1780 1785

9. TiemeierHvan TuijlHRHofmanAKiliaanAJBretelerMM 2003 Plasma fatty acid composition and depression are associated in the elderly: the Rotterdam Study. Am J Clin Nutr 78 40 46

10. FeartCPeuchantELetenneurLSamieriCMontagnierD 2008 Plasma eicosapentaenoic acid is inversely associated with severity of depressive symptomatology in the elderly: data from the Bordeaux sample of the Three-City Study. Am J Clin Nutr 87 1156 1162

11. SamieriCFeartCLetenneurLDartiguesJFPeresK 2008 Low plasma eicosapentaenoic acid and depressive symptomatology are independent predictors of dementia risk. Am J Clin Nutr 88 714 721

12. SchaeferEJBongardVBeiserASLamon-FavaSRobinsSJ 2006 Plasma phosphatidylcholine docosahexaenoic acid content and risk of dementia and Alzheimer disease: the Framingham Heart Study. Arch Neurol 63 1545 1550

13. SunQMaJCamposHHankinsonSEHuFB 2007 Comparison between plasma and erythrocyte fatty acid content as biomarkers of fatty acid intake in US women. Am J Clin Nutr 86 74 81

14. ChungHNettletonJALemaitreRNBarrRGTsaiMY 2008 Frequency and type of seafood consumed influence plasma (n-3) fatty acid concentrations. J Nutr 138 2422 2427

15. BurdgeGCWoottonSA 2003 Conversion of alpha-linolenic acid to palmitic, palmitoleic, stearic and oleic acids in men and women. Prostaglandins Leukot Essent Fatty Acids 69 283 290

16. GoyensPLSpilkerMEZockPLKatanMBMensinkRP 2006 Conversion of alpha-linolenic acid in humans is influenced by the absolute amounts of alpha-linolenic acid and linoleic acid in the diet and not by their ratio. Am J Clin Nutr 84 44 53

17. GoyensPLSpilkerMEZockPLKatanMBMensinkRP 2005 Compartmental modeling to quantify alpha-linolenic acid conversion after longer term intake of multiple tracer boluses. J Lipid Res 46 1474 1483

18. LemaitreRNSiscovickDSBerryEMKarkJDFriedlanderY 2008 Familial aggregation of red blood cell membrane fatty acid composition: the Kibbutzim Family Study. Metabolism 57 662 668

19. SchaefferLGohlkeHMullerMHeidIMPalmerLJ 2006 Common genetic variants of the FADS1 FADS2 gene cluster and their reconstructed haplotypes are associated with the fatty acid composition in phospholipids. Hum Mol Genet 15 1745 1756

20. RzehakPHeinrichJKloppNSchaefferLHoffS 2009 Evidence for an association between genetic variants of the fatty acid desaturase 1 fatty acid desaturase 2 (FADS1 FADS2) gene cluster and the fatty acid composition of erythrocyte membranes. Br J Nutr 101 20 26

21. XieLInnisSM 2008 Genetic variants of the FADS1 FADS2 gene cluster are associated with altered (n-6) and (n-3) essential fatty acids in plasma and erythrocyte phospholipids in women during pregnancy and in breast milk during lactation. J Nutr 138 2222 2228

22. MalerbaGSchaefferLXumerleLKloppNTrabettiE 2008 SNPs of the FADS gene cluster are associated with polyunsaturated fatty acids in a cohort of patients with cardiovascular disease. Lipids 43 289 299

23. BokorSDumontJSpinnekerAGonzalez-GrossMNovaE 2010 Single nucleotide polymorphisms in the FADS gene cluster are associated with delta-5 and delta-6 desaturase activities estimated by serum fatty acid ratios. J Lipid Res 51 2325 2333

24. ZietemannVKrogerJEnzenbachCJansenEFritscheA 2010 Genetic variation of the FADS1 FADS2 gene cluster and n-6 PUFA composition in erythrocyte membranes in the European Prospective Investigation into Cancer and Nutrition-Potsdam study. Br J Nutr 1 12

25. OkadaTSatoNFKuromoriYMiyashitaMIwataF 2006 Thr-encoding allele homozygosity at codon 54 of FABP 2 gene be associated with impaired delta 6 desatruase activity and reduced plasma arachidonic acid in obese children. J Atheroscler Thromb 13 192 196

26. TanakaTShenJAbecasisGRKisialiouAOrdovasJM 2009 Genome-wide association study of plasma polyunsaturated fatty acids in the InCHIANTI Study. PLoS Genet 5 e1000338 doi:10.1371/journal.pgen.1000338

27. PsatyBMO'DonnellCJGudnasonVLunettaKLFolsomAR 2009 Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) Consortium: Design of prospective meta-analyses of genome-wide association studies from 5 cohorts. Circ Cardiovasc Genet 2 73 80

28. HodsonLSkeaffCMFieldingBA 2008 Fatty acid composition of adipose tissue and blood in humans and its use as a biomarker of dietary intake. Prog Lipid Res 47 348 380

29. VittinghoffEGliddenDShiboskiSMcCullochC 2005 Regression Methods in Biostatistics: Linear Logistic, Survival and Repeated Measures Models. New York Springer

30. BokorSDumontJSpinnekerAGonzalez-GrossMNovaE 2010 Single nucleotide polymorphisms in the FADS gene cluster are associated with delta-5 and delta-6 desaturase activities estimated by serum fatty acid ratios. J Lipid Res

31. Molto-PuigmartiCPlatJMensinkRPMullerAJansenE 2010 FADS1 FADS2 gene variants modify the association between fish intake and the docosahexaenoic acid proportions in human milk. Am J Clin Nutr 91 1368 1376

32. LiouYAKingDJZibrikDInnisSM 2007 Decreasing linoleic acid with constant alpha-linolenic acid in dietary fats increases (n-3) eicosapentaenoic acid in plasma phospholipids in healthy men. J Nutr 137 945 952

33. TangCChoHPNakamuraMTClarkeSD 2003 Regulation of human delta-6 desaturase gene transcription: identification of a functional direct repeat-1 element. J Lipid Res 44 686 695

34. HindorffLJunkinsHHallPMehtaJManolioT A Catalog of Published Genome-Wide Association Studies. Available at: www.genome.gov/gwastudies

35. DupuisJLangenbergCProkopenkoISaxenaRSoranzoN 2010 New genetic loci implicated in fasting glucose homeostasis and their impact on type 2 diabetes risk. Nat Genet 42 105 116

36. WillerCJSannaSJacksonAUScuteriABonnycastleLL 2008 Newly identified loci that influence lipid concentrations and risk of coronary artery disease. Nat Genet 40 161 169

37. KathiresanSWillerCJPelosoGMDemissieSMusunuruK 2009 Common variants at 30 loci contribute to polygenic dyslipidemia. Nat Genet 41 56 65

38. AulchenkoYSRipattiSLindqvistIBoomsmaDHeidIM 2009 Loci influencing lipid levels and coronary heart disease risk in 16 European population cohorts. Nat Genet 41 47 55

39. HarrisWSBulchandaniD 2006 Why do omega-3 fatty acids lower serum triglycerides? Curr Opin Lipidol 17 387 393

40. FriedbergCEJanssenMJHeineRJGrobbeeDE 1998 Fish oil and glycemic control in diabetes. A meta-analysis. Diabetes Care 21 494 500

41. LandsWE 1958 Metabolism of glycerolipides; a comparison of lecithin and triglyceride synthesis. J Biol Chem 231 883 888

42. RissanenTVoutilainenSNyyssonenKLakkaTASalonenJT 2000 Fish oil-derived fatty acids, docosahexaenoic acid and docosapentaenoic acid, and the risk of acute coronary events: the Kuopio ischaemic heart disease risk factor study. Circulation 102 2677 2679

43. OdaEHatadaKKatohKKodamaMNakamuraY 2005 A case-control pilot study on n-3 polyunsaturated fatty acid as a negative risk factor for myocardial infarction. Int Heart J 46 583 591

44. WuCOrozcoCBoyerJLegliseMGoodaleJ 2009 BioGPS: an extensible and customizable portal for querying and organizing gene annotation resources. Genome Biol 10 R130

45. BertrandtJKlosADebskiB 2004 Influence of vitamin B6 supplementation on polyunsaturated fatty acids concentration in serum and liver of rats fed a diet restricted in protein. Nahrung 48 99 103

46. BordoniAHreliaSLorenziniABergamiRCabriniL 1998 Dual influence of aging and vitamin B6 deficiency on delta-6-desaturation of essential fatty acids in rat liver microsomes. Prostaglandins Leukot Essent Fatty Acids 58 417 420

47. BurdgeGCFinneganYEMinihaneAMWilliamsCMWoottonSA 2003 Effect of altered dietary n-3 fatty acid intake upon plasma lipid fatty acid composition, conversion of [13C]alpha-linolenic acid to longer-chain fatty acids and partitioning towards beta-oxidation in older men. Br J Nutr 90 311 321

48. LiYWillerCJDingJScheetPAbecasisGR 2010 MaCH: using sequence and genotype data to estimate haplotypes and unobserved genotypes. Genet Epidemiol 34 816 834

49. ServinBStephensM 2007 Imputation-based analysis of association studies: candidate regions and quantitative traits. PLoS Genet 3 e114 doi:10.1371/journal.pgen.0030114

50. BrowningBLBrowningSR 2009 A unified approach to genotype imputation and haplotype-phase inference for large data sets of trios and unrelated individuals. Am J Hum Genet 84 210 223

51. MarchiniJHowieBMyersSMcVeanGDonnellyP 2007 A new multipoint method for genome-wide association studies by imputation of genotypes. Nat Genet 39 906 913

Štítky
Genetika Reprodukčná medicína

Článok vyšiel v časopise

PLOS Genetics


2011 Číslo 7
Najčítanejšie tento týždeň
Najčítanejšie v tomto čísle
Kurzy

Zvýšte si kvalifikáciu online z pohodlia domova

Aktuální možnosti diagnostiky a léčby litiáz
nový kurz
Autori: MUDr. Tomáš Ürge, PhD.

Všetky kurzy
Prihlásenie
Zabudnuté heslo

Zadajte e-mailovú adresu, s ktorou ste vytvárali účet. Budú Vám na ňu zasielané informácie k nastaveniu nového hesla.

Prihlásenie

Nemáte účet?  Registrujte sa

#ADS_BOTTOM_SCRIPTS#