- Research Paper
- Published:
A genome-wide association study of n-3 and n-6 plasma fatty acids in a Singaporean Chinese population
Genes & Nutrition volume 10, Article number: 53 (2015)
Abstract
Polyunsaturated fatty acids (PUFAs) have a major impact on human health. Recent genome-wide association studies (GWAS) have identified several genetic loci that are associated with plasma levels of n-3 and n-6 PUFAs in primarily subjects of European ancestry. However, the relevance of these findings has not been evaluated extensively in other ethnic groups. The primary aim of this study was to evaluate for genetic loci associated with n-3 and n-6 PUFAs and to validate the role of recently identified index loci using data from a Singaporean Chinese population. Using a GWAS approach, we evaluated associations with plasma concentrations of three n-3 PUFAs [alphalinolenic acid (ALA), eicosapentaenoic acid and docosahexaenoic acid], four n-6 PUFAs [linoleic acid (LA), gammalinolenic acid, dihomogammalinolenic acid (DGLA) and arachidonic acid], and estimates of delta-5 desaturase and delta-6 desaturase activities among the participants (N = 1361) of the Singaporean Chinese Health Study. Our results reveal robust genome-wide associations (p value <5 × 10−8) with ALA, all four n-6 PUFAs, and delta-6 desaturase activity at the FADS1/FADS2 locus. We further replicated the associations between common index variants at the NTAN1/PDXDC1 locus and n-6 PUFAs LA and DGLA, and between the JMJD1C locus and n-6 PUFA LA (p value between 0.0490 and 9.88 × 10−4). These associations were independent of dietary intake of PUFAs. In aggregate, we show that genetic loci that influence plasma concentrations of n-3 and n-6 PUFAs are shared across different ethnic groups.
Introduction
Dietary changes aimed to reduce the risk of cardiovascular diseases traditionally involve limiting saturated fatty acids intakes, replacing them in part by polyunsaturated fatty acids (PUFAs). These PUFAs are characterized by multiple desaturations in their aliphatic tail and are classified as n-3 or n-6, depending on the position of the first double bond proximate to the methyl end. PUFAs in plasma and cell membranes originate from the diet or are synthesized from the two essential PUFAs, linoleic acid [LA; 18:2n-6] and alphalinolenic acid [ALA; 18:3n-3], through a series of desaturation and elongation processes (Youdim et al. 2000).
The association of dietary and membrane long-chain n-3 PUFAs with a lower risk of cardiovascular disease (Dewailly et al. 2001; Albert et al. 2002; Yokoyama et al. 2007; Wang et al. 2006) as well as non-cardiovascular diseases (Cherubini et al. 2007; Lorente-Cebrián et al. 2013) has been documented in various studies. For example, a diet rich in n-3 PUFAs, especially marine derived n-3 PUFAs—eicosapentaenoic acid (EPA; 20:5n3) and docosahexaenoic acid (DHA; c22:6n3), associated with high seafood intake is thought to be responsible for the lower incidence of heart disease in Inuit populations (Dewailly et al. 2001). An important and yet unanswered question pertains to the importance of n-6 PUFAs on cardiovascular risk. Some studies have suggested that these may convey beneficial health effects (Wu et al. 2014); however, such findings are still controversial (Harris and Shearer 2014). For example, certain n-6 PUFAs, such as LA, can be converted in humans to arachidonic acid [AA; 20:4n-6], and the latter serves as the precursor of important molecules involved in inflammation (e.g., eicosanoids) that is thought to play a role in the atherosclerotic process (Wu et al. 2014).
Data from family and twin studies suggest high heritability of circulating levels of these groups of PUFAs (Shah et al. 2009). Candidate gene studies focusing on desaturase genes cluster, fatty acid desaturase 1 and 2 (FADS1 and FADS2), have demonstrated significant associations with multiple n-3 and n-6 PUFAs (Mathias et al. 2010; Xie and Innis 2008). GWAS studies of n-3 and n-6 plasma levels conducted in an Italian population-based dataset (Tanaka et al. 2009) have found significant associations between rs174537 near FADS1 and levels of eicosadienoic acid (EDA; 20:2n6) and EPA. Another SNP mapped to chromosome 6 in the region encoding an elongase of very long fatty acids 2 (ELOVL2) was also found to be associated with EPA.
Recently, two large-scale meta-analyses of GWAS from the Cohorts for Heart and Aging Research in Genomic Epidemiology Consortium (CHARGE, 14–15) were conducted to identify common genetic variants associated with plasma n-3 and n-6 PUFA phenotypes. In the first study, findings showed that common variation in n-3 metabolic pathway genes [FADS1, FADS2 and ELOVL2] and in glucokinase regulator (GCKR) influence plasma phospholipid levels of n-3 PUFAs in these populations of European ancestry (Lemaitre et al. 2011). In the second study the CHARGE consortium confirmed previous findings of the FADS cluster with n-6 PUFA plasma levels and observed five additional loci that influence n-6 PUFA plasma concentrations (Guan et al. 2014).
These recent GWAS, however, have been carried out primarily in populations of European ancestry (especially the recent discoveries for n-6 PUFA associations) and may have missed loci unique to other populations. Furthermore, the significance of these recent findings has not been evaluated extensively in other ethnicities (e.g., East-Asian populations). Therefore, in this study, we performed a GWAS for n-3 and n-6 PUFA plasma levels and sought to evaluate the role of recently identified common index variants in a Singaporean Chinese population.
Materials and methods
Study population
This study utilized data from a case–control GWAS study for myocardial infarction (MI, N = 1540) that was nested within the Singapore Chinese Health Study (SCHS), a prospective cohort study of 63,257 Singaporean Chinese men and women aged 45–74 years living in Singapore between 1993 and 1998 (Hankin et al. 2001; Koh et al. 2011). The cohort study recruited only participants belonging to one of the two major Chinese dialect groups in Singapore, the Hokkiens or the Cantonese, who originated from two contiguous prefectures in southern China.
In-person interviews and phlebotomy were conducted before the onset of disease. Baseline assessment was conducted through a face-to-face structured interview during recruitment. Information was recorded by a trained interviewer using a structured questionnaire that included demographics, medical history, cigarette smoking, alcohol consumption, physical activity and detailed menstrual and reproductive history (women only). Food frequency questionnaires were used to estimate usual diet among the SCHS participants and daily proportions of macronutrients intake including n-3 and n-6 PUFA intakes were estimated (Hankin et al. 2001). Between April 1994 and December 1999, we collected blood and single-void urine specimens from a random 3 % sample of study enrollees. Starting in January 2000, the bio-specimen collection was extended to all surviving cohort members. By April 2005 when all subjects had been contacted, bio-specimens were collected from 32,543 participants, representing an approximately 60 % consent rate. Various components (plasma, red blood cell, serum and white blood cell) of blood were separated and have been stored continuously at −80 °C. Age was defined at the time of blood collection, and sex and ethnicity was recorded as indicated in the participant’s National Registration Identity Card (NRIC).
Cases of myocardial infarction (MI) patients (N = 744) were identified through the Singapore Registry of Births and Deaths and the Hospital Discharge Database, respectively, and who had donated blood for research before suffering MI or death. MI cases were confirmed by review of medical records by a cardiologist and through criteria recommended by the Multi-Ethnic Study of Atherosclerosis (available at: http://www.mesa-nhlbi.org/manuals.aspx) or through the Singapore Myocardial Infarction Registry which uses similar procedures. Controls (N = 796) were participants who were alive and free of coronary heart disease at the time of the diagnosis or death of the index cases and were matched for sex, dialect group, year of birth (±2 years), year of recruitment (±1 year) and date of blood collection (±6 months). Informed consent was obtained from all subjects and approval by the Institutional Review Board of the National University of Singapore was given for this study.
Main outcome variables
Plasma subtypes of n-3 and n-6 PUFAs were measured using GC/MS/MS analysis implemented on the Agilent 7890 GC system (Shanghai, China) equipped with a 7001B QQQ triple quadrupole mass detector (PA, USA). Three plasma n-3 PUFAs [ALA, EPA and DHA] and four n-6 PUFAs [LA, GLA (Gammalinolenic acid; c18:3n6), DGLA (Dihomogammalinolenic acid; c20:3n6) and AA] were available for analysis (Table 1). Delta-5 desaturase (D5D) was estimated as the AA:DGLA ratio and delta-6 desaturase (D6D) as the DGLA:LA ratio (Warensjö et al. 2009) (Table 1).
Genotyping and imputation
Genome-wide genotyping was performed for SCHS samples using the Illumina HumanOmni ZhongHua-8 Bead Chip. Quality control (QC) procedures of samples and SNPs are detailed in Supplementary Table 1 and 2. Samples with call-rate <98.0 % (N = 8), extremes in heterozygosity (N = 11) and known duplicates (positive controls, N = 42) were excluded from analyses. Correlation between positive controls were >98.0 %. Identity-by-state measures were performed by pair-wise comparison of samples to detect cryptic relatedness (Supplementary Table 1) and one sample from each relationship was excluded from further analysis (N = 16). Principle component analysis (PCA) (Price et al. 2006) together with HapMap reference populations (International HapMap Consortium 2003) was performed to identify possible outliers from reported ethnicity and 18 samples were excluded. 717 MI cases and 644 controls were available for subsequent statistical analysis (Supplementary Table 1).
Details of quality-control measures for SNPs are provided in Supplementary Table 2. Briefly, sex-linked and mitochondrial SNPs were removed, together with gross HWE outliers (p value <1 × 10−5). SNPs that were monomorphic or with a MAF <1.0 % and SNPs with low call-rates (<95.0 %) were excluded. We imputed for additional autosomal SNPs with IMPUTE v2 (Marchini et al. 2007) using a posterior probability of at least 0.90 and a call-rate of at least 95.0 % for all imputed SNPs. Genotype calls were based on HapMap East-Asian samples (CHB and JPT) and imputed SNPs with info score <0.5 were excluded from analyses. 1,972,057 SNPs in the MI cases and 1,980,188 SNPs in the MI controls were available for statistical analyses after QC procedures (Supplementary Table 2).
Statistical analysis
GWAS association analysis for PUFA subtypes
As performed in previous studies, raw PUFA concentrations were converted to proportions of total fatty acids (Lemaitre et al. 2011; Guan et al. 2014). All data were normalized by natural log transformation and truncated to 4 SD from respective means. SNP-based trend tests for PUFA associations as well as D5D and D6D activities were carried out separately in the MI cases and controls. Associations between genetic variants (SNPs) and plasma n-3 or n-6 PUFA concentrations were analyzed in an additive model and adjusted for age, sex, population stratification (the first 3 principal components) and proportion of daily n-3 or n-6 PUFA intake from total calories (derived from food frequency questionnaires) under investigation. These analyses were performed using the genome association toolset, SNPTEST (version 2) (Marchini et al. 2007). Genomic inflation factor (λ) of association results from each analysis was used to evaluate levels of inflation of study results. Individual study results from MI cases and controls were subsequently combined using the inverse variance-weighted meta-analysis, assuming a fixed effects model. Heterogeneity of effects was examined using Cochran’s Q, and a Q p value cut-off <0.05 was used to determine SNPs with between-study heterogeneity (Marchini et al. 2007). All meta-analysis procedures were performed using the meta package in STATA (version 11.2). At loci where genome-wide signals were identified, we further evaluated for additional independent SNPs by repeating all analyses and conditioning SNP associations on the genotypes of the index SNP at that region. For lead SNPs with genome-wide association identified in our study, we estimated the proportion of the respective PUFA variance explained by these SNPs by evaluating adjusted R 2 values of regression models before and after inclusion of the lead SNP into the model. Finally, at genome-wide loci, we also evaluated changes in effects after further adjustments for concentrations of the precursor fatty acids.
In silico replication analysis of previously reported candidate SNPs
All recently discovered index plasma n-3 PUFA loci (7 SNPs in or near 3 loci, FADS1, ELOVL2 and GCKR) and index plasma n-6 PUFA loci [5 SNPs in or near 4 loci, FADS1, jumonji domain containing 1C (JMJD1C), N-terminal asparagine amidase (NTAN1) and palmitoyl-protein thioesterase 2 (PPT2)] from recent GWAS studies (Lemaitre et al. 2011; Guan et al. 2014) were selected for in silico replication among the SCHS GWAS data. QUANTO power calculator (Gauderman and Morrison 2006) was used to estimate the power to detect similar associations of index SNPs in the SCHS dataset using reported effect estimates (Lemaitre et al. 2011; Guan et al. 2014) and observed MAF from the SCHS datasets.
To determine significance, consistency of effect direction with previously reported GWAS data was ensured in the SCHS datasets and as these loci have been well-validated [in 8866 subjects of European ancestry (Lemaitre et al. 2011; Guan et al. 2014)], a p value threshold <0.05 was used (Dorajoo et al. 2013; Dorajoo et al. 2012). At loci that failed to replicate, we evaluated associations of SNPs within the recombination block containing the reported index SNP (between 2 recombination peaks of at least 20Â cM/Mb) (Dorajoo et al. 2012). These analyses were performed using LocusZoom ver1.1 (Pruim et al. 2010).
Results
Robust genome-wide signals at FADS1/2 locus in chromosome 11 detected for ALA, LA, GLA, DGLA and AA associations and delta-6 desaturase activity in the Singapore Chinese
Supplementary Fig. 1–3 shows QQ plots for all meta-analysis association results of 7 n-3 and n-6 PUFAs and 2 desaturase activities assessed in the study. Genomic inflation in these meta-analysis results was low (λ between 1.004 and 1.019 for n-3 PUFA associations, between 1.004 and 1.006 for n-6 PUFA associations and between 1.002 and 1.005 for desaturase activities) indicating that our QC procedures had effectively excluded aberrant samples and SNPs. Supplementary Tables 3–11 provide association results for SNPs with meta-analysis p value associations <1.0 × 10−4 for the 9 traits evaluated in SCHS.
SNPs at the chr11q12.2 locus, in or near the FADS1/FADS2 gene, were detected at genome-wide levels for ALA, LA, GLA, DGLA, AA and D6D activity associations in our Singaporean Chinese dataset (Fig. 1; Table 2). These signals explained between 2.0 % (for ALA) and 36.5 % (for GLA) of the phenotypic variation in PUFA traits assessed. All genome-wide SCHS lead SNPs (SNPs with strongest associations) at the chr11q12.2 locus were either the same index SNP reported in previous European studies for these traits (rs174547) (Lemaitre et al. 2011; Guan et al. 2014) or were in strong LD (r 2 between 0.973 and 1.000) with this European index SNP (Table 2). Conditioning analyses on the genotype of the index SNP (rs174547) did not reveal additional independent signals at the chr11q12.2 locus in our study (Supplementary Fig. 4). Although strong genome-wide associations were observed at this locus for D6D, no associations were detected for D5D in the SCHS (Supplementary Fig. 5). Further adjustment of results with preceding plasma PUFA concentrations also did not affect associations of lead SNPs except at rs174548 for DGLA association, where the effect estimate was observed to be in the opposite direction (as compared to the unadjusted study results) (Table 2).
We also identified a novel SNP at chr10q24.2 for the association with ALA at marginal genome-wide significance levels in SCHS (rs7903338, meta-analysis p value = 4.33 × 10−8, Supplementary Table 3), but this did not replicate in 8866 European subjects from the Cohorts for Heart and Aging Research in Genomic Epidemiology Consortium (CHARGE), (β = −0.0002, p value = 0.825, C allele frequency = 0.743, http://faculty.washington.edu/rozenl/files/CHARGE_N3_ALA.txt).
Replication of index PUFA subtypes SNPs in the Singaporean Chinese
With the SCHS PUFA subtype data, we were able to replicate known index SNP associations at 3 out of 4 n-3 PUFAs (ALA, EPA and DHA) and 4 out of 5 n-6 PUFAs (LA, GLA, DGLA and AA). Tables 3 and 4 details association levels of these recently identified index n-3 and n-6 PUFA SNPs in the SCHS meta-analysis data.
Index SNPs at the FADS1/FADS2 gene region on chromosome 11 were validated for all PUFA subtype traits available in SCHS (p value between 0.0011 and 1.68 × 10−170, Tables 3 and 4). For n-3 subtypes these FADS1 index SNPs showed opposite effect directions for ALA as compared with EPA (Table 3). For n-6 fatty acids, effect directions of LA and DGLA associations were opposite to that of GLA and AA associations (Table 4). Besides these SNPs, rs16966952 at the NTAN1/PDXDC1 locus was also validated for increased LA and DGLA concentrations (p value between 0.0490 and 9.88 × 10−4, Table 4) in the SCHS. At index SNPs where we did not observe any direct replication for the specific traits they were reported to be associated with, we examined associations for SNPs within the recombination block containing the reported index SNPs (Supplementary Fig. 6). Only at the JMJD1C locus for LA associations, we identified a regional SNP with nominal associations (rs7896783, p value = 0.0098, Supplementary Fig. 6C) in strong LD with the reported index SNP (rs10740118, r 2 = 0.82).
Among the set of known index SNPs, we also identified nominal associations for additional PUFA subtypes available in the SCHS. These include increased levels of DHA (T allele, p value = 0.0019) and decreased levels of LA (T allele, p value = 0.0050) at the GCKR index SNP rs780094 (Tables 3, 4), increased levels of DHA (G allele, p value = 0.0321) at the JMJD1C index SNP rs10740118 (Table 3) and decreased levels of ALA (C allele, p value = 0.0493) at the PPT index SNP rs3134950 (Table 3).
Discussion
PUFAs have a major impact on human health (Youdim et al. 2000; Dewailly et al. 2001; Albert et al. 2002; Yokoyama et al. 2007; Wang et al. 2006; Cherubini et al. 2007; Lorente-Cebrián et al. 2013; Wu et al. 2014). The present study details associations of common genetic variants with specific n-3 and n-6 plasma PUFA concentrations in the Singaporean Chinese population. It is known that PUFA plasma concentrations can originate from dietary sources and that dietary choices vary greatly among different ethnic populations (Ogce et al. 2008). For example in East-Asians, where rice is a staple, diet is usually higher in starch and carbohydrates as compared to European populations (Drewnowski and Popkin 1997). In our study we corroborated and further characterized the role of numerous common genetic loci, first identified in European populations among the Singaporean Chinese and suggest that these identified associations were independent of dietary n-3 and n-6 PUFA intakes.
We confirmed the importance of common variants in/near the central desaturase genes, FADS1 and FADS2 on chromosome 11, which showed robust associations with n-3 PUFA (ALA) and n-6 PUFAs (LA, GLA, DGLA and AA) in the Singaporean Chinese. The associations were likely to stem from a single causal locus and conditional probability analyses did not reveal independent associations. It was only for the DGLA association that we observed differences in effect estimates when adjustment for precursor PUFA subtype concentrations was made (GLA concentration). The especially strong associations with GLA at this locus, explaining approximately 36Â % of this phenotypic variance, and the known positive correlations between GLA and DGLA concentrations may explain these varied results and it is likely that adjustment with plasma GLA concentrations refined the DGLA genetic association signal observed at the FADS1/FADS2 locus in our study.
Increased delta-6 desaturase activity is associated with an increased plasma GLA:LA ratio. Consistently, the T allele of the rs174547 index SNP in the region of the FADS1 gene was associated with higher GLA and lower LA concentrations among the Singaporean Chinese (Hankin et al. 2001) and showed robust genome-wide associations with increased indirect estimate of delta-6 activity. This indicates the prominent role these desaturation enzymes play in the conversion of the LA essential fatty acid to GLA. The same allele also showed robust associations with higher AA and lower DGLA concentrations in our dataset, as previously reported (Guan et al. 2014), suggesting a higher delta-5 desaturase activity. However, we do not find evidence that the FADS1/2 locus is associated with delta-5 desaturase activity. This lack of association between delta-5 desaturase activity, and the FADS1/2 locus is also consistent with a previous study that analyzed subjects of European ancestry (Bokor et al. 2010). Nevertheless, it is noteworthy that a rare variant (rs968567) in weak LD with the common variants evaluated in our study (r 2Â =Â 0.4 between rs968567 and rs174547 in HapMap CEU population; rs968567 not available for HapMap East-Asian reference population) at the FADS2 promoter site may still be associated with delta-5 desaturase (Bokor et al. 2010).
Besides the validations for n-3 and n-6 PUFA at the FADS1/2 locus, we further replicated SNPs at the NTAN1/PDXDC1 on chromosome 13 (for LA and DGLA associations) and the JMJD1C locus on chromosome 10 (for LA associations) for n-6 PUFAs. Importantly, these recent findings for plasma n-6 PUFA subtypes associations have been only evaluated in subjects of European ancestry (Guan et al. 2014). The PDXDC1 locus has been identified in a separate genome-wide study for human serum metabolite concentrations (LA/PUFA ratio) and may be involved in eicosanoid biosynthesis (Kettunen et al. 2012). JMJDC1 encodes a probable histone demethylase and has recently been shown to be involved in hematopoiesis (Soranzo et al. 2009) as well as triglyceride levels (Teslovich et al. 2010). In aggregate, our data suggest that diverse ethnic groups share some degree of genetic predisposition to both n-3 and n-6 PUFA metabolism.
For a few loci, however, we were unable to detect significant associations, namely ELVOL2 SNP associations for EPA and DHA n-3 PUFA concentrations and the NTAN1/PDXDC1 SNP associations for GLA and AA n-6 PUFA concentrations. One possibility is that there is significant genetic heterogeneity at these loci between East-Asians and Europeans. However, index SNPs are unlikely to be causal and differences in allele frequencies and LD between populations may have masked associations. Moreover, due to modest effects of these variants, the sample size of our study may not have been adequately powered to detect all associations. The loci that failed to replicate tended to be those for which we had lower power in our study. Further large-scale studies in non-European populations would be warranted to fully elucidate the role of all common variants involved in n-3 and n-6 PUFA metabolism and to evaluate for possible ethnic specific associations.
One advantage of our study was the availability of a number of n-3 and n-6 PUFAs for the evaluation of possible pleiotropic effects at these established genetic loci. Besides the FADS1/FADS2 locus that was strongly associated with all n-3 and n-6 PUFA concentrations, perhaps of note is the GCKR locus (T allele of rs780094) that has been implicated with DPA concentration (Lemaitre et al. 2011). The T allele of this SNP has been associated with decreased blood glucose and insulin levels as well as with increased blood lipids (Teslovich et al. 2010; Dupuis et al. 2010; Wallace et al. 2008; Coram et al. 2013). We further identified nominal associations with both increased n-3 DHA and decreased n-6 LA concentrations among the Singaporean Chinese. Given the strong associations between DHA and LA subtypes on blood lipid levels, the GCKR locus may potentially play a central role modulating both glucose metabolism and lipogenesis pathways that impact upon plasma n-3 and n-6 PUFA concentrations in varied ways. These finding should be verified in additional large-scale studies, and future functional studies should aim to tease out the true pleiotropic role of GCKR in these metabolic pathways.
Our study is not without limitations and an important one would be the modest sample size of our study population. Thus we were unlikely to be well-powered to identify additional novel associations (especially associations with modest effect sizes) for the traits we analyzed in this Singaporean Chinese population. For ALA association, however, we identified a novel locus at chr10q24.2 at marginal genome-wide levels. This signal, however, did not replicate among 8866 European subjects from the CHARGE consortium and given that the SNP was also common in the European population (MAFÂ =Â 25.7Â %), this finding is likely to be a false positive one due to random statistical fluctuation and highlights the importance to follow-up novel leads with additional independent replications. Nevertheless, at the specific genetic loci for n-3 and n-6 PUFA identified recently, we were primarily adequately powered to detect similar associations identified first in Europeans, among our Singaporean Chinese dataset and this limitation does not impact on the main finding of our study that a proportion of genetic predisposition influencing plasma n-3 and n-6 PUFA concentrations are shared across ethnic groups, at least among European and Chinese ancestry populations.
In conclusion, we firmly corroborated the role of FADS1/FADS2 common variants as determinants of plasma n-3 and n-6 PUFA concentrations and show that this locus is an important determinant of delta-6 desaturase activity in the Singaporean Chinese population. Several additional associations at the NTAN1/PDXDC1 and the JMJD1C loci were also replicated for n-6 PUFA levels for the first time in a non-European ancestry population, suggesting that similar common genetic variants may play a role as primary determinants of plasma PUFA concentrations among Chinese and European ancestry populations.
References
Albert CM, Campos H, Stampfer MJ, Ridker PM, Manson JE, Willett WC, Ma J (2002) Blood levels of long-chain n-3 fatty acids and the risk of sudden death. N Engl J Med 346(15):1113–1118
Bokor S, Dumont J, Spinneker A, Gonzalez-Gross M, Nova E, Widhalm K, Moschonis G, Stehle P, Amouyel P, De Henauw S et al (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(8):2325–2333
Cherubini A, Andres-Lacueva C, Martin A, Lauretani F, Iorio AD, Bartali B, Corsi A, Bandinelli S, Mattson MP, Ferrucci L (2007) Low plasma N-3 fatty acids and dementia in older persons: the InCHIANTI study. J Gerontol A Biol Sci Med Sci 62(10):1120–1126
Coram MA, Duan Q, Hoffmann TJ, Thornton T, Knowles JW, Johnson NA, Ochs-Balcom HM, Donlon TA, Martin LW, Eaton CB et al (2013) Genome-wide characterization of shared and distinct genetic components that influence blood lipid levels in ethnically diverse human populations. Am J Hum Genet 92(6):904–916
Dewailly E, Blanchet C, Lemieux S, Sauve L, Gingras S, Ayotte P, Holub BJ (2001) n-3 Fatty acids and cardiovascular disease risk factors among the Inuit of Nunavik. Am J Clin Nutr 74:464–473
Dorajoo R, Blakemore AI, Sim X, Ong RT, Ng DP, Seielstad M, Wong TY, Saw SM, Froguel P, Liu J et al (2012) Replication of 13 obesity loci among Singaporean Chinese, Malay and Asian-Indian populations. Int J Obes (Lond) 36(1):159–163
Dorajoo R, Li R, Ikram MK, Liu J, Froguel P, Lee J, Sim X, Ong RT, Tay WT, Peng C et al (2013) Are C-reactive protein associated genetic variants associated with serum levels and retinal markers of microvascular pathology in Asian populations from Singapore? PLoS One 8(7):e67650
Drewnowski A, Popkin BM (1997) The nutrition transition: new trends in the global diet. Nutr Rev 55(2):31–43
Dupuis J, Langenberg C, Prokopenko I, Saxena R, Soranzo N, Jackson AU, Wheeler E, Glazer NL, Bouatia-Naji N, Gloyn AL et al (2010) New genetic loci implicated in fasting glucose homeostasis and their impact on type 2 diabetes risk. Nat Genet 42(2):105–116
Gauderman WJ, Morrison JM (2006) QUANTO 1.1: a computer program for power and sample size calculations for genetic-epidemiology studies. http://hydra.usc.edu/gxe/
Guan W, Steffen BT, Lemaitre RN, Wu JH, Tanaka T, Manichaikul A, Foy M, Rich SS, Wang L, Nettleton JA et al (2014) Genome-wide association study of plasma N6 polyunsaturated fatty acids within the cohorts for heart and aging research in genomic epidemiology consortium. Circ Cardiovasc Genet 7(3):321–331
Hankin JH, Stram DO, Arakawa K, Park S, Low SH, Lee HP, Yu MC (2001) Singapore Chinese Health Study: development, validation, and calibration of the quantitative food frequency questionnaire. Nutr Cancer 39(2):187–195
Harris WS, Shearer GC (2014) Omega-6 fatty acids and cardiovascular disease: friend, not foe? Circulation 130(18):1562–1564
International HapMap Consortium (2003) The international HapMap project. Nature 426(6968):789–796
Kettunen J, Tukiainen T, Sarin AP, Ortega-Alonso A, Tikkanen E, Lyytikäinen LP, Kangas AJ, Soininen P, Würtz P, Silander K et al (2012) Genome-wide association study identifies multiple loci influencing human serum metabolite levels. Nat Genet 44(3):269–276
Koh WP, Yuan JM, Wang R, Lee YP, Lee BL, Yu MC, Ong CN (2011) Plasma carotenoids and risk of acute myocardial infarction in the Singapore Chinese Health Study. Nutr Metab Cardiovasc Dis 21:685–690
Lemaitre RN, Tanaka T, Tang W, Manichaikul A, Foy M, Kabagambe EK, Nettleton JA, King IB, Weng LC, Bhattacharya S et al (2011) 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):e1002193
Lorente-Cebrián S, Costa AG, Navas-Carretero S, Zabala M, MartÃnez JA, Moreno-Aliaga MJ (2013) Role of omega-3 fatty acids in obesity, metabolic syndrome, and cardiovascular diseases: a review of the evidence. J Physiol Biochem 69(3):633–651
Marchini J, Howie B, Myers S, McVean G, Donnelly P (2007) A new multipoint method for genome-wide association studies via imputation of genotypes. Nat Genet 39:906–913
Mathias RA, Vergara C, Gao L, Rafaels N, Hand T, Campbell M, Bickel C, Ivester P, Sergeant S, Barnes KC et al (2010) FADS genetic variants and omega-6 polyunsaturated fatty acid metabolism in a homogeneous island population. J Lipid Res 51(9):2766–2774
Ogce F, Ceber E, Ekti R, Oran NT (2008) Comparison of mediterranean, Western and Japanese diets and some recommendations. Asian Pac J Cancer Prev 9(2):351–356
Price AL, Patterson NJ, Plenge RM, Weinblatt ME, Shadick NA, Reich D (2006) Principal components analysis corrects for stratification in genome-wide association studies. Nat Genet 38(8):904–909
Pruim RJ, Welch RP, Sanna S, Teslovich TM, Chines PS, Gliedt TP, Boehnke M, Abecasis GR, Willer CJ (2010) LocusZoom: regional visualization of genome-wide association scan results. Bioinformatics 26(18):2336–2337
Shah SH, Hauser ER, Bain JR, Muehlbauer MJ, Haynes C, Stevens RD, Wenner BR, Dowdy ZE, Granger CB, Ginsburg GS et al (2009) High heritability of metabolomic profiles in families burdened with premature cardiovascular disease. Mol Syst Biol 5:258
Soranzo N, Spector TD, Mangino M, Kühnel B, Rendon A, Teumer A, Willenborg C, Wright B, Chen L, Li M et al (2009) A genome-wide meta-analysis identifies 22 loci associated with eight hematological parameters in the HaemGen consortium. Nat Genet 41(11):1182–1190
Tanaka T, Shen J, Abecasis GR, Kisialiou A, Ordovas JM, Guralnik JM, Singleton A, Bandinelli S, Cherubini A, Arnett D et al (2009) Genome-wide association study of plasma polyunsaturated fatty acids in the InCHIANTI Study. PLoS Genet 5(1):e1000338
Teslovich TM, Musunuru K, Smith AV, Edmondson AC, Stylianou IM, Koseki M, Pirruccello JP, Ripatti S, Chasman DI, Willer CJ et al (2010) Biological, clinical and population relevance of 95 loci for blood lipids. Nature 466(7307):707–713
Wallace C, Newhouse SJ, Braund P, Zhang F, Tobin M, Falchi M, Ahmadi K, Dobson RJ, Marçano AC et al (2008) Genome-wide association study identifies genes for biomarkers of cardiovascular disease: serum urate and dyslipidemia. Am J Hum Genet 82(1):139–149
Wang C, Harris WS, Chung M, Lichtenstein AH, Balk EM, Kupelnick B, Jordan HS, Lau J (2006) n-3 Fatty acids from fish or fish-oil supplements, but not alpha-linolenic acid, benefit cardiovascular disease outcomes in primary- and secondary-prevention studies: a systematic review. Am J Clin Nutr 84(1):5–17
Warensjö E, Rosell M, Hellenius ML, Vessby B, De Faire U, Risérus U (2009) Associations between estimated fatty acid desaturase activities in serum lipids and adipose tissue in humans: links to obesity and insulin resistance. Lipids Health Dis 8:37
Wu JH, Lemaitre RN, King IB, Song X, Psaty BM, Siscovick DS, Mozaffarian D (2014) Circulating omega-6 polyunsaturated fatty acids and total and cause-specific mortality: the Cardiovascular Health Study. Circulation 130(15):1245–1253
Xie L, Innis SM (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(11):2222–2228
Yokoyama M, Origasa H, Matsuzaki M, Matsuzawa Y, Saito Y, Ishikawa Y, Oikawa S, Sasaki J, Hishida H, Itakura H et al (2007) Effects of eicosapentaenoic acid on major coronary events in hypercholesterolaemic patients (JELIS): a randomised open-label, blinded endpoint analysis. Lancet 369(9567):1090–1098
Youdim KA, Martin A, Joseph JA (2000) Essential fatty acids and the brain: possible health implications. Int J Dev Neurosci 18(4–5):383–399
Acknowledgments
We thank Siew-Hong Low of the National University of Singapore for supervising the fieldwork and Renwei Wang for the maintenance of the database in the Singapore Chinese Health Study. We also thank the Ministry of Health in Singapore for assistance with the identification of fatal and non-fatal myocardial infarction cases via database linkages. Finally, we acknowledge the founding Principal Investigator of the Singapore Chinese Health Study, Mimi C. Yu. The Singapore Chinese Health Study was supported by the National Institutes of Health, USA (RO1 CA144034 and UM1 CA182876), the nested case–control study of myocardial infarction by the Singapore National Medical Research Council (NMRC 1270/2010) and genotyping by the HUJ-CREATE Programme of the National Research Foundation, Singapore (Project Number 370062002).
Author information
Authors and Affiliations
Corresponding authors
Ethics declarations
Conflict of interest
All authors have no conflicts of interest.
Electronic supplementary material
Below is the link to the electronic supplementary material.
Rights and permissions
About this article
Cite this article
Dorajoo, R., Sun, Y., Han, Y. et al. A genome-wide association study of n-3 and n-6 plasma fatty acids in a Singaporean Chinese population. Genes Nutr 10, 53 (2015). https://doi.org/10.1007/s12263-015-0502-2
Received:
Accepted:
Published:
DOI: https://doi.org/10.1007/s12263-015-0502-2