Identification of Novel Genetic Determinants of Erythrocyte Membrane Fatty Acid Composition among Greenlanders
Disruption of fatty-acid balance has in several previous studies been linked to human health conditions, including the metabolic syndrome, type 2 diabetes, and insulin resistance. Composition of fatty acids in lipid membranes is influenced, not only by diet and lifestyle, but also by genetic variation. By identifying genes linked to changes in the level of specific fatty acids, it may be possible to identify biological mechanisms and pathways central to regulation of fatty-acid composition in lipid membranes. We therefore aimed at finding such genes by studying Greenlanders. We identified six genomic regions harboring variants, which were associated with the level of at least one of 22 assessed erythrocyte membrane fatty acids, including two novel regions not previously linked to fatty acid levels. Moreover, we showed that two of the identified variants were associated with altered levels of glycosylated hemoglobin, and one of these variants was associated with reduced insulin resistance and decreased measures of body size. These results contribute to our understanding of fatty acid metabolism, and support a link between fatty acid balance and metabolic health.
Published in the journal:
. PLoS Genet 12(6): e32767. doi:10.1371/journal.pgen.1006119
Category:
Research Article
doi:
https://doi.org/10.1371/journal.pgen.1006119
Summary
Disruption of fatty-acid balance has in several previous studies been linked to human health conditions, including the metabolic syndrome, type 2 diabetes, and insulin resistance. Composition of fatty acids in lipid membranes is influenced, not only by diet and lifestyle, but also by genetic variation. By identifying genes linked to changes in the level of specific fatty acids, it may be possible to identify biological mechanisms and pathways central to regulation of fatty-acid composition in lipid membranes. We therefore aimed at finding such genes by studying Greenlanders. We identified six genomic regions harboring variants, which were associated with the level of at least one of 22 assessed erythrocyte membrane fatty acids, including two novel regions not previously linked to fatty acid levels. Moreover, we showed that two of the identified variants were associated with altered levels of glycosylated hemoglobin, and one of these variants was associated with reduced insulin resistance and decreased measures of body size. These results contribute to our understanding of fatty acid metabolism, and support a link between fatty acid balance and metabolic health.
Introduction
Fatty acids (FAs) are important for normal body function, as they serve as essential structural entities of cellular membranes, as energy sources, and as signaling molecules. Perturbation of FA homeostasis may modulate membrane functions, cell signaling, and gene expression, and regulation of FA metabolism is, thus, critically important. Accordingly, variation in FA levels, measured in serum, plasma, erythrocyte membranes, or as dietary intake, has been linked to a range of cardiovascular [1–6] and metabolic risk factors, including the metabolic syndrome, dyslipidemia, type 2 diabetes, insulin resistance, and inflammation [7–10]. Moreover, significant correlation between erythrocyte membrane FAs and serum lipid levels has been observed [11–13], and this correlation has been suggested to be the link between erythrocyte membrane FA composition and metabolic diseases [13].
FAs are obtained from diet or synthesized endogenously from carbohydrate or protein sources by series of elongation and desaturation steps. ω-3 alpha-linolenic acid (18:3) and ω-6 linoleic acid (cis-cis-18:2) are essential FAs, which can only be obtained from diet and subsequently elongated and desaturated to form other ω-3 and ω-6 FAs. The levels of FAs in an individual are influenced by diet and lifestyle [14,15], but also have a clear hereditary component, which is estimated to account for 32–70% of FA variation [16,17]. The metabolic pathways regulating the circulating concentrations and membrane content of individual FAs, as well as the specific mechanisms linking FA levels to disease states are, however, poorly understood. We hypothesize, that improved understanding of these pathways and mechanisms can be achieved by identifying genetic loci associated with inter-individual differences in FA levels in erythrocyte membranes.
For genetic studies, the Greenlandic population is an important and powerful resource. This population has emerged from a historically small and isolated Inuit population, which has not until very recently admixed with Europeans [18]. Therefore, compared to European populations, the Greenlandic population has extended linkage disequilibrium (LD) and increased probability for presence of high frequency harmful variants due to genetic drift. These properties are advantageous in genetic studies, as they increase the statistical power to detect association signals. This was recently illustrated, when a novel type 2 diabetes associated variant, with an unusually high odds ratio of 10.3, was identified in a relatively small sample of Greenlanders [19]. The Greenlandic population also differs markedly from European populations with respect to diet. The traditional diet of the Greenlanders is rich in polyunsaturated ω-3 FAs derived mainly from marine mammals and fish [20], and has recently been shown to have had a large impact on genetic makeup through adaptive selection in Inuit [21].
In the present study, we aimed to identify novel genetic loci associated with inter-individual differences in the FA composition in erythrocyte membranes in a large-scale association study in Greenlanders.
Results
The erythrocyte membrane FA-association analyses were conducted in 2,626 Greenlanders living in Greenland (57.5% women) from the population-based Greenlandic Inuit Health in Transition (IHIT) cohort. These individuals were on average 44.7 years old and had an average BMI of 26.4 kg/m2. We assessed the levels of 22 FAs in the phospholipid fraction of erythrocyte membranes, and analyzing MetaboChip genotyping data, we identified six independent association signals with a p-value under the significance threshold of 4.3x10-7 (S1 Fig). These six signals, mapped to genomic loci on chromosome 5, 11, 12, 19, and 20 (Table 1). For each of these regions we assessed imputed data for fine mapping, and assessed secondary FA associations as well as associations with metabolic phenotypes (S2 Fig). In the following, all associations are reported for the derived alleles in chromosomal order, and associations are reported down to the arbitrary p-value cut-off of p<1.0x10-3.
The first association signal mapped to a novel locus on chromosome 5. The lead SNP (rs251015) was located in intronic region in FNIP1, and the derived allele showed strongest association with higher levels of oleic acid (18:1 ω-9). The imputation data revealed another SNP in the locus, rs76430747, located in an intronic region of ACSL6. The derived allele of rs76430747 was associated with lower levels of oleic acid (18:1 ω-9), and higher levels of lignoceric acid (24:0) and 11-eicosenoic acid (20:1 ω-9; Table 2 and Fig 1). Conditional analyses showed that rs76430747 may explain the observed association with oleic acid (Fig 2), whereas the results for lignoceric acid and 11-eicosenoic acid were less clear (S3 Fig). Interestingly, the derived allele of rs76430747 was, in addition to FA levels, associated with higher HbA1c (Fig 3 and Table 3).
The second and the third association signals mapped to chromosome 11. The strongest signal from the MetaboChip data (rs174570) was located in intronic sequence of the well-established FA-associated FADS2. rs174570 showed association with the level of 13 FAs (Table 2). The second signal on chromosome 11, mapped a little more than 7Mb away from the FADS region to intronic sequence of TPCN2. Similar to the FADS2 variant, the MetaboChip lead SNP in this locus (rs1551304) showed association with a large number of FAs (S4 Table). Imputation based analyses revealed no better candidate SNPs in the region, instead we genotyped a known missense variant (rs80356779) in the nearby candidate gene CPT1A, attempting to identify the causal variant for this second locus.
The CPT1A rs80356779 variant accentuated the association pattern across all assessed FA pathways observed for rs1551304, however, in the opposite direction for the derived allele, but the same direction for the minor alleles. Also, it mirrored the association pattern observed for the rs174570 FADS2 variant (Fig 1, Table 2). Specifically, the CPT1A rs80356779 showed a pattern where the levels of ω-7 palmitoleic acid (16:1) was increased, the downstream products of the ω-3 (docosapentaenoic acid, 22:5), ω-6 (arachidonic acid, 20:4 and adrenic acid, 22:4), ω-9 (nervonic acid, 24:1), and saturated fatty acid (SFA; lignoceric acid, 24:0) pathways were reduced, and the levels of the corresponding upstream precursors (ω-6, linoleic acid (cis-cis-18:2), gamma-linolenic acid (18:3), and dihomo-gamma-linolenic acid (20:3); ω-9, oleic acid (18:1), 11-eicosenoic acid (20:1), and erucic acid (22:1); SFA, arachidic acid (20:0) and behenic acid (22:0)) were increased (Fig 1 and Table 2). In addition to the FA associations, CPT1A rs80356779 was associated with lower HbA1c (Fig 3), indicating improved glycemic regulation, and with several measures of smaller body size and reduced insulin resistance (Table 3).
Further exploration of the region revealed an unusual LD phenomenon, where FADS2 rs174570 and CPT1A rs80356779 were linked in the ancestral Inuit population despite the chromosomal distance of approximately 7 Mb. Both variants were fixed or at extremely low frequency for the derived allele in the ancestral Inuit population, while the ancestral allele of rs80356779 was entirely fixed, and the ancestral allele of rs174570 was very common in European populations (Table 4). Conditional analysis reduced the association signal for both loci, as seen for 11-eicosenoic acid (20:1 ω-9; Fig 4). In general, most of the CTP1A associations were still significant when conditioning on the FADS2 variant, while many of the FADS2 associations disappeared when conditioning on the FADS2 variant (Fig 4 and S4 Fig).
The fourth association signal mapped to chromosome 12, where the MetaboChip lead SNP (rs2110073) was located in intronic sequence in PHB2. This variant was associated with higher levels of oleic acid (18:1 ω-9), and lower levels of linoleic acid (cis-cis-18:2 ω-6; Fig 1 and Table 2). Imputation based analyses revealed additional SNPs in the region in high LD with rs2110073 mapping to nearby genes including the known FA-associated LPCAT3 (S5 Fig). However, none of these SNPs had a significantly better p-value than rs2110073, and were therefore not analyzed further. Conditional analyses indicated that rs2110073 may explain the associations observed for oleic acid (18:1 ω-9) and linoleic acid (cis-cis-18:2 ω-6; S5 Fig)
The fifth association signal mapped to chromosome 19. The lead SNP (rs2913968) was located in intronic region of RAB11B, and was associated with the SFA arachidic acid (20:0). Imputation analyses pointed to another possible causal SNP, rs11881630, mapping to CERS4. Additional genotyping of this synonymous exome variant showed a stronger association signal with decreased levels of arachidic acid than observed for rs2913968 (Fig 1 and Table 2). However, association analyses conditioning on rs2913968 and rs11881630, respectively, were inconclusive in determining which, if any, of the variants that caused the association signal (S6 Fig), and with the available data it was impossible to exclude the possibility of two independent association signals.
The sixth and final association signal mapped to a novel locus on chromosome 20, where the MetaboChip lead SNP (rs6035106) was located in an intronic region of DTD1. The derived allele of this SNP was associated with higher levels of docosahexaenoic acid (22:6 ω-3), and in addition showed an interesting pattern of associations, with generally higher ω-3 and ω-6 FA levels, but lower SFA levels (Fig 1 and Table 2). The imputation data revealed no candidate SNPs in the region with a significantly better p-value, and the conditional analysis on rs6035106 showed that this variant may explain the observed docosahexaenoic acid (22:6 ω-3) association, as well as the other FA associations (S7 Fig).
Discussion
We identified six distinct genomic regions harboring variants significantly associated with the level of at least one of 22 assessed erythrocyte membrane FAs. These genomic regions comprise two novel loci not previously linked to FA levels (ACSL6, DTD1), one known FA-linked locus where we identify a broad range of novel FA associations (CPT1A), and three loci for which we replicate known FA-associations (FADS2, LPCAT3, CERS4).
The first of the two novel FA-associated loci (rs76430747) is located on in ACSL6 chromosome 5. The ACSL6 protein is a long-chain acyl-CoA synthase that activates FAs by catalyzing the formation of acyl-CoA from ATP, CoA, and FAs [22]. This step is required for FAs to enter most metabolic pathways including protein modification, phospholipid synthesis, and beta-oxidation. ACSL proteins, thus, affect the distribution and amount of intracellular FAs, and the isoforms differ in tissue specificity and FA preference [23]. ACSL6 is mainly expressed in brain and neuronal cells, and catalyzes very-long chain FAs, containing 18–26 carbons [24]. Our results may indicate that rs76430747 causes a shift in the preference for FA incorporated into phospholipids in cellular membranes towards ω-9 11-eicosenoic acid (20:1) and the SFA lignoceric acid (24:0), as opposed to ω-9 oleic acid (18:1), which seems to be directed to degradation by beta-oxidation. Alternatively, ACSL6 facilitated activation of oleic acid (18:1) to its corresponding acyl-CoA may channel this type of FA to elongation by ELOVL3 [25], explaining the observed increase in 11-eicosenoic acid (20:1) levels. Interestingly, ACSL6 has been shown to activate oleic acid in rat cell [26,27], which supports our functional hypothesis.
Besides from their role as intermediates in lipid synthesis and FA degradation, saturated and unsaturated long-chain acyl-CoAs have also been linked to regulation of carbohydrate and lipid metabolism, gene expression, and insulin secretion [28,29]. These functional properties may explain the observed association between variation in ACSL6 and increased HbA1c. However, this association was not supported by differences in measures of insulin secretion or glucose levels, and requires replication to be validated.
The second novel FA-associated locus we identified is located on chromosome 20. The lead SNP (rs6035106) showed an interesting pattern of altered FA concentrations, with higher levels of ω-3 and ω-6 FAs, and reduced levels of SFAs, however, only associated with modest p-values. Based on this pattern we hypothesize that the causal variant in this locus induces a shift in the choice of FAs absorbed from diet, or in the specificity towards FAs incorporated into erythrocyte membrane lipids. The observed association signals mapped to DTD1, however, the product of this gene has no known functional link to FA-metabolic pathways [30]. The gene-rich region contains other possible causal candidate genes, among these are SCP2D1. This gene is particularly interesting, as it encodes a protein functionally related to the protein product of SCP2 [31–33]. The SCP2 protein has been proposed to redistribute numerous lipid species, including FAs and fatty-acyl CoAs between lipid rafts and intracellular organelles, and thus regulate signaling pathways [31]. Additional studies are needed to validate the putative association between this genomic region and FA levels, and to determine if SCP2D1 is the causal gene in the region.
On chromosome 11 we identified a locus showing a broad range of novel FA associations. The lead SNP in this region is a missense variant (rs80356779) in CPT1A, located around 7 Mb from the FADS-gene region. The FADS genes have been shown to be associated with a range of FAs primarily in European populations in candidate gene studies [34–38], and in large-scale genetic studies [39–44], and these associations were replicated in this study. The CPT1A rs80356779 variant causes a proline to leucine substitution at position 479 in the CPT1A protein sequence (P479L), and has been studied previously in Arctic populations [45–47]. Conditional analyses indicate that CPT1A P479L is the causal variant in the region among Inuit. However, the FADS locus must also harbor a true association signal, as such a signal has been observed in many studies in European populations [46–48], where P479L is monomorphic. The FADS and CPT1A association signals are difficult to untangle, because rs174570 and rs80356779 are both fixed in the ancestral Inuit population. This unusual long-range LD phenomenon is mainly due to recent adaptive selection on these loci, and lack of admixture from other populations until only a couple of hundred years ago [21,49].
CPT1A encodes a member of the carnitine palmitoyltransferase (CPT) protein family. These proteins catalyzes the initial and rate-limiting step in mitochondrial degradation of FAs in skeletal muscle, liver, and adipose tissue [50], and are thus attractive candidates to affect the levels of individual FAs. The CPT1A isoform is mainly expressed in liver [51], and cell-line studies and animal models have indicated that the P479L variant results in increased CPT1A activity, and thereby possibly increased flow of FAs to degradation by beta-oxidation [52,53]. We hypothesize, that increased CPT1A activity is the functional explanation for the observed changes in FA concentrations across the ω-3, ω-6, ω-7, ω-9, and SFA pathways. This extensive association pattern is a novel observation. Previously only association with the ω-6 FAs linoleic acid (cis-cis-18:2), dihomo-gamma-linolenic acid (20:3), and arachidonic acid (20:4) have been reported for three CPT1A intron variants [54]. Our results indicate that the Inuit specific P479L is the causal variant in CPT1A explaining the observed FA associations.
Moreover, CPT1A P479L T-allele carriers (L479) showed reduced insulin resistance and smaller body size. This phenotypic impact of L479 is an extension of the previously reported association with reduced body fat and central obesity in Yup’ik Eskimos [47], and it is well in line with data from functional studies showing that increased FA oxidation is linked to reduced accumulation of body fat and improved lipid profile [55–58], whereas inhibition of fat oxidation by a CPT1 inhibitor in rodents causes increased adiposity and insulin resistance [59]. We also observed association between L479 and improved glycemic regulation, which may simply be a consequence of the genotypic effect on body size and insulin resistance. In the Yup’ik Eskimos as well as in Greenlandic Inuit L479 has been associated with higher circulating levels of HDL cholesterol [46,47], this was not replicated in our study despite our larger sample size. The genomic region containing CPT1A seems to have been under increased evolutionary pressure, and to be of particular importance in Arctic populations [48,49]. This may be due to the traditional Greenlandic diet, which is enriched with animal fat, accentuating the requirement for FA beta-oxidation facilitated by CPT1A for generating energy.
In addition to the novel associations and the well-established FADS2 association discussed above we also replicated the association between rs2110073, near the functional candidate gene LPCAT3 on chromosome 12, and altered levels of the first precursors for endogenous synthesis of the ω-6 (linoleic acid, cis-cis-18:2) and ω-9 (oleic acid, 18:1) FAs [44], and the association between variation in CERS4 and levels of arachidic acid (20:0) [60–62]. The originally reported CERS4 SNP, rs2100944 [60], barely reached nominal significance in our study (p = 0.013), instead we identified another causal candidate SNP in this locus, namely rs11881630. This variant is rare in Europeans (MAF: 0.89%), and is in low LD with rs2100944 both in Europeans and in Greenlanders, thus, indicating that neither of the two is the causal variant, or that there may be two association signals.
Conclusions and perspectives
We have identified six independent loci associated with erythrocyte membrane FA levels in Greenlanders. The Greenlandic population is characterized by geographical isolation, and has a unique genetic profile shaped by genetic drift due to both founder events and a small population size through thousands of years. Moreover, the Greenlanders have adapted their lifestyle to cold climate and available resources, resulting in a highly specialized diet rich in fat and protein from fish and marine mammals [20]. This extreme Greenlandic lifestyle has resulted in genetic adaptation, demonstrated by identification of positive selection in the FADS genes on chromosome 11 [21]. The importance of this region in adaption to Arctic lifestyle is further supported by positive selection on CPT1A in Siberians [49]. Besides from the selection signatures in the FADS and CPT1A loci, we observed no indications of selection acting on the FA-associated regions in our data. Thus, based on our data it seems unlikely that the unique FA composition of the traditional Greenlandic diet has inferred selective pressure on the ACSL6, LPCAT3, CERS4, and DTD1 gene regions.
In conclusion, we have identified genetic determinants of FA composition of erythrocyte membrane phospholipids possibly as a result of altered fluxes through FA and lipid metabolic pathways. These alterations may have metabolic consequences, supported by the observed link between FA associated variants and altered glycemic regulation and body composition. The present study provides a framework to further delineate how specific lipid species regulate human metabolic disorders.
Materials and Methods
Ethics statement
The study was approved by the Commission for Scientific Research in Greenland (project 2011–13, ref. no. 2011–056978; and project 2013–13, ref.no. 2013–090702), and was conducted in accordance with the ethical standards of the Declaration of Helsinki, second revision. All participants gave informed consent.
Study population
The study included individuals from three Greenlandic cohorts: the IHIT (n = 3,115) and B99 (n = 1,401) cohorts, which comprise Greenlanders living in Greenland, and the BBH (n = 547) cohort, which comprises Greenlanders living in Denmark. The B99 and the IHIT cohorts were collected as part of a general population health survey of the Greenlandic population during 1999–2001 and 2005–2010, respectively, as described in [63,64]. BBH was collected in Denmark during 1998–1999 from people of Greenlandic descent [64]. There was an overlap of individuals between IHIT and B99, these individuals (n = 295) were assigned to the B99 cohort.
Measurements and assays
Fatty acids (FAs) in phospholipids of erythrocytes were measured in 2,626 individuals from the IHIT cohort of Greenlanders living in Greenland. FA levels were reported as relative levels compared to the total amount of FAs in each sample. Lipids were separated by thin layer chromatography, after extraction of total lipids with chloroform/methanol, and FAs in the phospholipid fraction were methylated, and subsequently analyzed by capillary GLC using a HP-Packard GC chromatograph equipped with a DB-23 column at The Centre de recherche sur les maladies lipidiques (CRML), Centre hospitalier universitaire de Québec. We restricted the genetic analyses to 22 FAs for which we observed detectable levels for the majority of individuals, and which could be placed in the FA synthesis pathways (Fig 1). These FAs comprised ω-3 (18:3, alpha-linolenic acid; 18:4, stearidonic acid; 20:4, eicosatetraenoic acid; 20:5, eicosapentaenoic acid; 22:5, docosapentaenoic acid; 22:6, docosahexaenoic acid), ω-6 (cis-cis-18:2, linoleic acid; 18:3, gamma-linolenic acid; 20:3, dihomo-gamma-linolenic acid; 20:4, arachidonic acid; 22:4, adrenic acid; 22:5, docosapentaenoic acid), ω-7 (16:1, palmitoleic acid), ω-9 (18:1, oleic acid; 20:1, 11-eicosenoic acid; 22:1, erucic acid; 24:1, nervonic acid), and SFAs (16:0, palmitic acid; 18:0, stearic acid; 20:0, arachidic acid; 22:0, behenic acid; 24:0, lignoceric acid).
All IHIT and B99 participants above 35 years of age underwent an oral glucose tolerance test, where blood samples were drawn after an overnight fast of at least 8 hours, and 2 hours after receiving a 75 g glucose bolus. Plasma glucose levels were analyzed with the Hitachi 912 system (Roche Diagnostics), serum insulin with an immunoassay excluding des-31,32 split products and intact proinsulin (AutoDELFIA, PerkinElmer), and Hba1c by ion-exchange HPLC (G7, Tosoh Bioscience). Insulin resistance was estimated by the homeostasis model assessment (HOMA-IR), calculated as [(fasting glucose level x fasting insulin level)/6.945]/22.5, where insulin levels were expressed as pmol/l and glucose levels as mmol/l [65]. Information about diet was obtained from questionnaires, as described in [66].
Genotyping
DNA was purified from blood leukocytes, and all samples were genotyped on the Metabochip (Illumina). This chip is a custom iSelect genotyping array of 196,725 single nucleotide polymorphisms (SNPs) for genetic studies of metabolic, cardiovascular, and anthropometric traits [67]. Genotyping was performed using the HiScan system (Illumina), and genotypes were called jointly for all cohorts using the GenCall module of the GenomeStudio software (Illumina) using default cluster data. The data set went through a two-step quality control. In step one, duplicate samples and individuals missing >2% genotypes or with gender discrepancy were removed. In step two, we removed SNPs with a minor allele frequency <1% (<5% for case-control transformed phenotypes), with >100 missing genotypes, with a large deviation from Hardy Weinberg equilibrium (p<1.0x10-10), as well as SNPs which were polymorphic in the IHIT cohort but not in the B99+BBH cohorts, and SNPs associated with gender (p<1.0x10-5). In total, 4,674 individuals (2,791 from IHIT, 1,336 from B99, and 547 from BBH) and 115,182 SNPs passed the quality control.
MetaboChip data was used to identify loci harboring association signals (p<4.3x10-7, corresponding to correction for testing 115,182 SNPs). Based on these data and imputation analyses of the identified loci, three additional SNPs were selected for genotyping (rs76430747, rs11881630, and rs11028474). Moreover, based on previous reports of association with lipid and body composition traits in the Alaskan Yup’ik Eskimos [47] and high frequency in Canadian and Greenlandic Inuit [45,46] an additional candidate SNP in the CPT1A locus (rs80356779), was also selected for genotyping (LGC genomics). Analyzing genotype data against MetaboChip data for rs11028474 showed discrepancy, hence, this variant was removed from all analyses.
Imputation
In order to fine map loci identified to harbor an association signal in the MetaboChip data, we carried out imputation to increase the number of SNPs in these regions. For the imputation we applied Omni5M array (Illumina) genotype data from 20 Greenlandic trios, which was phased using ShapeIt [68]. The phased data of the 40 Greenlandic parents and 1000 genomes data for 41 Europeans and 40 Han Chinese were applied as reference panel. IMPUTE2 [69] was used for imputation, where a recombination map for the reference SNPs was inferred with linear interpolation using the hg19 genomic map from IMPUTE2 as a template, and an effective population size of 1,500 [21]. The imputation generated genotype data (R2 >0.4) for 1,959,225 variants, which were analyzed as dosages (mean genotype) using GEMMA with the same setup as for the other statistical analyses using a relatedness matrix based on the dosages instead of genotypes, see below.
Statistical analysis
For association tests we applied a linear mixed model, implemented in the software GEMMA [70], to account for relatedness and admixture. For each phenotype the tests were applied to data from all individuals across the three cohorts with information about that specific phenotype, and the relatedness matrix required as input to GEMMA was estimated from genotypic data from these individuals only. For all tests we assumed an additive effect and included sex, age, and cohort as covariates, and for tests involving FAs we included dietary intake of marine mammals, and percentage of the diet consisting of traditional Inuit food as additional covariates. The FA association analyses were also run without the adjustment for diet, and with additional adjustment for the use of lipid-lowering drugs to exclude confounding from these factors. These adjustments did not alter the results, but are included in the supplementary tables for the FA-associated variants for comparison (S1–S10 Tables).
Prior to performing association tests, quantitative traits were quantile transformed to a standard normal distribution within each sex. However, six FAs (alpha-linolenic acid (18:3 ω-3), stearidonic acid (18:4 ω-3), eicosatetraenoic acid (20:4 ω-3), gamma-linolenic acid (18:3 ω-6), docosapentaenoic acid (22:5 ω-6), and erucic acid (22:1 ω-9)) had undetectable levels in more than 20% of the samples. We, thus, considered transformation to a standard normal distribution inappropriate. Instead these traits were transformed to case-control like data, by transforming values different from 0 to 1, and they were analyzed as detectable levels (cases) and undetectable levels (controls).
In analyses of the MetaboChip data, we applied a significance threshold of p<4.3x10-7, corresponding to correction for testing 115,182 SNPs, to identify loci harboring association signals. These loci were selected for further analyses, including imputation efforts, and assessment of secondary FA associations and metabolic phenotypes (S1–S10 Tables and S2 Fig). All associations are described for the derived alleles, and reported down to the arbitrary p-value cut-off of p<1.0x10-3.
Supporting Information
Zdroje
1. Siscovick DS, Raghunathan TE, King I, Weinmann S, Wicklund KG, Albright J, et al. Dietary intake and cell membrane levels of long-chain n-3 polyunsaturated fatty acids and the risk of primary cardiac arrest. JAMA. 1995;274: 1363–1367. 7563561
2. Lemaitre RN, King IB, Mozaffarian D, Sotoodehnia N, Rea TD, Kuller LH, et al. Plasma phospholipid trans fatty acids, fatal ischemic heart disease, and sudden cardiac death in older adults: the cardiovascular health study. Circulation. 2006;114: 209–215. doi: 10.1161/CIRCULATIONAHA.106.620336 16818809
3. Sun Q, Ma J, Campos H, Rexrode KM, Albert CM, Mozaffarian D, et al. Blood concentrations of individual long-chain n-3 fatty acids and risk of nonfatal myocardial infarction. Am J Clin Nutr. 2008;88: 216–223. 18614744
4. Yamagishi K, Folsom AR, Steffen LM. Plasma fatty acid composition and incident ischemic stroke in middle-aged adults: the Atherosclerosis Risk in Communities (ARIC) Study. Cerebrovasc Dis. 2013;36: 38–46. doi: 10.1159/000351205 23920478
5. Sala-Vila A, Cofan M, Perez-Heras A, Nunez I, Gilabert R, Junyent M, et al. Fatty acids in serum phospholipids and carotid intima-media thickness in Spanish subjects with primary dyslipidemia. Am J Clin Nutr. 2010;92: 186–193. doi: 10.3945/ajcn.2009.28807 20463042
6. Simon JA, Hodgkins ML, Browner WS, Neuhaus JM, Bernert JT, Hulley SB. Serum fatty acids and the risk of coronary heart disease. Am J Epidemiol. 1995;142: 469–476. 7677125
7. Friedberg CE, Janssen MJ, Heine RJ, Grobbee DE. Fish oil and glycemic control in diabetes. A meta-analysis. Diabetes Care. 1998;21: 494–500. 9571330
8. Wang L, Folsom AR, Zheng Z-J, Pankow JS, Eckfeldt JH. Plasma fatty acid composition and incidence of diabetes in middle-aged adults: the Atherosclerosis Risk in Communities (ARIC) Study. Am J Clin Nutr. 2003;78: 91–98. 12816776
9. Warensjö E, Sundström J, Lind L, Vessby B. 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. 2006;84: 442–448. 16895896
10. Poudyal H, Brown L. Should the pharmacological actions of dietary fatty acids in cardiometabolic disorders be classified based on biological or chemical function? Prog Lipid Res. 2015;59: 172–200. doi: 10.1016/j.plipres.2015.07.002 26205317
11. Sun Q, Ma J, Campos H, Hankinson SE, Manson JE, Stampfer MJ, et al. A prospective study of trans fatty acids in erythrocytes and risk of coronary heart disease. Circulation. 2007;115: 1858–1865. doi: 10.1161/CIRCULATIONAHA.106.679985 17389261
12. Sun Q, Ma J, Campos H, Hankinson SE, Hu FB. Comparison between plasma and erythrocyte fatty acid content as biomarkers of fatty acid intake in US women. Am J Clin Nutr. 2007;86: 74–81. 17616765
13. Jacobs S, Schiller K, Jansen E, Fritsche A, Weikert C, di Giuseppe R, et al. Association between erythrocyte membrane fatty acids and biomarkers of dyslipidemia in the EPIC-Potsdam study. Eur J Clin Nutr. 2014;68: 517–525. doi: 10.1038/ejcn.2014.18 24569539
14. Kunesová M, Hainer V, Tvrzicka E, Phinney SD, Stich V, Parízková J, et al. Assessment of dietary and genetic factors influencing serum and adipose fatty acid composition in obese female identical twins. Lipids. 2002;37: 27–32. 11876260
15. Kunesová M, Phinney S, Hainer V, Tvrzická E, Stich V, Parízková J, et al. The responses of serum and adipose Fatty acids to a one-year weight reduction regimen in female obese monozygotic twins. Ann N Y Acad Sci. 2002;967: 311–323. 12079858
16. Lemaitre RN, Siscovick DS, Berry EM, Kark JD, Friedlander Y. Familial aggregation of red blood cell membrane fatty acid composition: the Kibbutzim Family Study. Metabolism. 2008;57: 662–668. doi: 10.1016/j.metabol.2007.12.011 18442630
17. Harris WS, Pottala J V, Lacey SM, Vasan RS, Larson MG, Robins SJ. Clinical correlates and heritability of erythrocyte eicosapentaenoic and docosahexaenoic acid content in the Framingham Heart Study. Atherosclerosis. 2012;225: 425–431. doi: 10.1016/j.atherosclerosis.2012.05.030 22727409
18. Moltke I, Fumagalli M, Korneliussen TS, Crawford JE, Bjerregaard P, Jørgensen ME, et al. Uncovering the Genetic History of the Present-Day Greenlandic Population. Am J Hum Genet. The American Society of Human Genetics; 2015;96: 54–69. doi: 10.1016/j.ajhg.2014.11.012
19. Moltke I, Grarup N, Jørgensen ME, Bjerregaard P, Treebak JT, Fumagalli M, et al. A common Greenlandic TBC1D4 variant confers muscle insulin resistance and type 2 diabetes. Nature. Nature Publishing Group; 2014;512: 190–193. doi: 10.1038/nature13425
20. Jeppesen C, Jørgensen ME, Bjerregaard P. Assessment of consumption of marine food in Greenland by a food frequency questionnaire and biomarkers. Int J Circumpolar Health. 2012;71: 18361. doi: 10.3402/ijch.v71i0.18361 22663940
21. Fumagalli M, Moltke I, Grarup N, Racimo F, Bjerregaard P, Jorgensen ME, et al. Greenlandic Inuit show genetic signatures of diet and climate adaptation. Science. 2015;349: 1343–1347. doi: 10.1126/science.aab2319 26383953
22. Soupene E, Kuypers FA. Mammalian long-chain acyl-CoA synthetases. Exp Biol Med (Maywood). 2008;233: 507–521. doi: 10.3181/0710-MR-287
23. Grevengoed TJ, Klett EL, Coleman RA. Acyl-CoA metabolism and partitioning. Annu Rev Nutr. 2014;34: 1–30. doi: 10.1146/annurev-nutr-071813-105541 24819326
24. Fujino T, Yamamoto T. Cloning and functional expression of a novel long-chain acyl-CoA synthetase expressed in brain. J Biochem. 1992;111: 197–203. 1569043
25. Jakobsson A, Westerberg R, Jacobsson A. Fatty acid elongases in mammals: their regulation and roles in metabolism. Prog Lipid Res. 2006;45: 237–249. doi: 10.1016/j.plipres.2006.01.004 16564093
26. Marszalek JR, Kitidis C, Dirusso CC, Lodish HF. Long-chain acyl-CoA synthetase 6 preferentially promotes DHA metabolism. J Biol Chem. 2005;280: 10817–10826. doi: 10.1074/jbc.M411750200 15655248
27. Van Horn CG, Caviglia JM, Li LO, Wang S, Granger DA, Coleman RA. Characterization of Recombinant Long-Chain Rat Acyl-CoA Synthetase Isoforms 3 and 6: Identification of a Novel Variant of Isoform 6. Biochemistry. 2005;44: 1635–1642. doi: 10.1021/bi047721l 15683247
28. Faergeman NJ, Knudsen J. Role of long-chain fatty acyl-CoA esters in the regulation of metabolism and in cell signalling. Biochem J. 1997;323 Pt 1: 1–12. 9173866
29. Neess D, Bek S, Engelsby H, Gallego SF, Færgeman NJ. Long-chain acyl-CoA esters in metabolism and signaling: Role of acyl-CoA binding proteins. Prog Lipid Res. 2015;59: 1–25. doi: 10.1016/j.plipres.2015.04.001 25898985
30. Pasaje CFA, Bae JS, Park B-L, Jang A-S, Uh S-T, Kim M-K, et al. Association analysis of DTD1 gene variations with aspirin-intolerance in asthmatics. Int J Mol Med. 2011;28: 129–137. doi: 10.3892/ijmm.2011.669 21479357
31. Schroeder F, Atshaves BP, McIntosh AL, Gallegos AM, Storey SM, Parr RD, et al. Sterol carrier protein-2: new roles in regulating lipid rafts and signaling. Biochim Biophys Acta. 2007;1771: 700–718. doi: 10.1016/j.bbalip.2007.04.005 17543577
32. Seedorf U, Ellinghaus P, Roch Nofer J. Sterol carrier protein-2. Biochim Biophys Acta. 2000;1486: 45–54. 10856712
33. Edqvist J, Blomqvist K. Fusion and fission, the evolution of sterol carrier protein-2. J Mol Evol. 2006;62: 292–306. doi: 10.1007/s00239-005-0086-3 16501878
34. Schaeffer L, Gohlke H, Müller M, Heid IM, Palmer LJ, Kompauer I, et al. 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. 2006;15: 1745–1756. doi: 10.1093/hmg/ddl117 16670158
35. Malerba G, Schaeffer L, Xumerle L, Klopp N, Trabetti E, Biscuola M, et al. SNPs of the FADS gene cluster are associated with polyunsaturated fatty acids in a cohort of patients with cardiovascular disease. Lipids. 2008;43: 289–299. doi: 10.1007/s11745-008-3158-5 18320251
36. Rzehak P, Heinrich J, Klopp N, Schaeffer L, Hoff S, Wolfram G, et al. 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. 2009;101: 20–26. doi: 10.1017/S0007114508992564 18479586
37. Lattka E, Illig T, Heinrich J, Koletzko B. FADS gene cluster polymorphisms: important modulators of fatty acid levels and their impact on atopic diseases. J Nutrigenet Nutrigenomics. 2009;2: 119–128. doi: 10.1159/000235559 19776639
38. Bokor S, Dumont J, Spinneker A, Gonzalez-Gross M, Nova E, Widhalm K, et al. 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. 2010;51: 2325–2333. doi: 10.1194/jlr.M006205 20427696
39. Tanaka T, Shen J, Abecasis GR, Kisialiou A, Ordovas JM, Guralnik JM, et al. Genome-wide association study of plasma polyunsaturated fatty acids in the InCHIANTI Study. PLoS Genet. 2009;5: e1000338. doi: 10.1371/journal.pgen.1000338 19148276
40. Lemaitre RN, Tanaka T, Tang W, Manichaikul A, Foy M, Kabagambe EK, et al. Genetic loci associated with plasma phospholipid n-3 fatty acids: a meta-analysis of genome-wide association studies from the CHARGE Consortium. PLoS Genet. 2011;7: e1002193. doi: 10.1371/journal.pgen.1002193 21829377
41. Wu JHY, Lemaitre RN, Manichaikul A, Guan W, Tanaka T, Foy M, et al. Genome-wide association study identifies novel loci associated with concentrations of four plasma phospholipid fatty acids in the de novo lipogenesis pathway: results from the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) consortiu. Circ Cardiovasc Genet. 2013;6: 171–183. doi: 10.1161/CIRCGENETICS.112.964619 23362303
42. Guan W, Steffen BT, Lemaitre RN, Wu JHY, Tanaka T, Manichaikul A, et al. 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. 2014;7: 321–331. doi: 10.1161/CIRCGENETICS.113.000208 24823311
43. Mozaffarian D, Kabagambe EK, Johnson CO, Lemaitre RN, Manichaikul A, Sun Q, et al. Genetic loci associated with circulating phospholipid trans fatty acids: a meta-analysis of genome-wide association studies from the CHARGE Consortium. Am J Clin Nutr. 2015;101: 398–406. doi: 10.3945/ajcn.114.094557 25646338
44. Tintle NL, Pottala JV, Lacey S, Ramachandran V, Westra J, Rogers a., et al. A genome-wide association study of saturated, mono- and polyunsaturated red blood cell fatty acids in the Framingham Heart Offspring Study. Prostaglandins, Leukot Essent Fat Acids. Elsevier; 2015;94: 65–72.
45. Park JY, Narayan SB, Bennett MJ. Molecular assay for detection of the common carnitine palmitoyltransferase 1A 1436(C>T) mutation. Clin Chem Lab Med. 2006;44: 1090–1091. doi: 10.1515/CCLM.2006.196 16958601
46. Rajakumar C, Ban MR, Cao H, Young TK, Bjerregaard P, Hegele RA. Carnitine palmitoyltransferase IA polymorphism P479L is common in Greenland Inuit and is associated with elevated plasma apolipoprotein A-I. J Lipid Res. 2009;50: 1223–1228. doi: 10.1194/jlr.P900001-JLR200 19181627
47. Lemas DJ, Wiener HW, O’Brien DM, Hopkins S, Stanhope KL, Havel PJ, et al. Genetic polymorphisms in carnitine palmitoyltransferase 1A gene are associated with variation in body composition and fasting lipid traits in Yup’ik Eskimos. J Lipid Res. 2012;53: 175–184. doi: 10.1194/jlr.P018952 22045927
48. Zhou S, Xiong L, Xie P, Ambalavanan A, Bourassa C V., Dionne-Laporte A, et al. Increased Missense Mutation Burden of Fatty Acid Metabolism Related Genes in Nunavik Inuit Population. Prokunina-Olsson L, editor. PLoS One. 2015;10: e0128255. doi: 10.1371/journal.pone.0128255 26010953
49. Clemente FJ, Cardona A, Inchley CE, Peter BM, Jacobs G, Pagani L, et al. A Selective Sweep on a Deleterious Mutation in CPT1A in Arctic Populations. Am J Hum Genet. 2014;95: 584–589. doi: 10.1016/j.ajhg.2014.09.016 25449608
50. McGarry JD, Brown NF. The mitochondrial carnitine palmitoyltransferase system. From concept to molecular analysis. Eur J Biochem. 1997;244: 1–14. 9063439
51. Swanson ST, Foster DW, McGarry JD, Brown NF. Roles of the N- and C-terminal domains of carnitine palmitoyltransferase I isoforms in malonyl-CoA sensitivity of the enzymes: insights from expression of chimaeric proteins and mutation of conserved histidine residues. Biochem J. 1998;335 Pt 3: 513–519. 9794789
52. Brown NF, Mullur RS, Subramanian I, Esser V, Bennett MJ, Saudubray JM, et al. Molecular characterization of L-CPT I deficiency in six patients: insights into function of the native enzyme. J Lipid Res. 2001;42: 1134–1142. 11441142
53. Akkaoui M, Cohen I, Esnous C, Lenoir V, Sournac M, Girard J, et al. Modulation of the hepatic malonyl-CoA—carnitine palmitoyltransferase 1A partnership creates a metabolic switch allowing oxidation of de novo fatty acids 1. Biochem J. 2009;420: 429–438. doi: 10.1042/BJ20081932 19302064
54. Voruganti VS, Higgins PB, Ebbesson SOE, Kennish J, Göring HHH, Haack K, et al. Variants in CPT1A, FADS1, and FADS2 are Associated with Higher Levels of Estimated Plasma and Erythrocyte Delta-5 Desaturases in Alaskan Eskimos. Front Genet. 2012;3: 86. doi: 10.3389/fgene.2012.00086 22701466
55. Buckley JD, Howe PRC. Anti-obesity effects of long-chain omega-3 polyunsaturated fatty acids. Obes Rev. 2009;10: 648–659. doi: 10.1111/j.1467-789X.2009.00584.x 19460115
56. Wang Y-X, Lee C-H, Tiep S, Yu RT, Ham J, Kang H, et al. Peroxisome-proliferator-activated receptor delta activates fat metabolism to prevent obesity. Cell. 2003;113: 159–170. 12705865
57. Zurlo F, Lillioja S, Esposito-Del Puente A, Nyomba BL, Raz I, Saad MF, et al. Low ratio of fat to carbohydrate oxidation as predictor of weight gain: study of 24-h RQ. Am J Physiol. 1990;259: E650–657. 2240203
58. Simoneau JA, Veerkamp JH, Turcotte LP, Kelley DE. Markers of capacity to utilize fatty acids in human skeletal muscle: relation to insulin resistance and obesity and effects of weight loss. FASEB J. 1999;13: 2051–2060. 10544188
59. Dobbins RL, Szczepaniak LS, Bentley B, Esser V, Myhill J, McGarry JD. Prolonged inhibition of muscle carnitine palmitoyltransferase-1 promotes intramyocellular lipid accumulation and insulin resistance in rats. Diabetes. 2001;50: 123–130. 11147777
60. Lemaitre RN, King IB, Kabagambe EK, Wu JHY, McKnight B, Manichaikul A, et al. Genetic loci associated with circulating levels of very long-chain saturated fatty acids. J Lipid Res. 2015;56: 176–184. doi: 10.1194/jlr.M052456 25378659
61. Hicks AA, Pramstaller PP, Johansson Å, Vitart V, Rudan I, Ugocsai P, et al. Genetic Determinants of Circulating Sphingolipid Concentrations in European Populations. Gibson G, editor. PLoS Genet. 2009;5: e1000672. doi: 10.1371/journal.pgen.1000672 19798445
62. Demirkan A, van Duijn CM, Ugocsai P, Isaacs A, Pramstaller PP, Liebisch G, et al. Genome-wide association study identifies novel loci associated with circulating phospho- and sphingolipid concentrations. PLoS Genet. 2012;8: e1002490. doi: 10.1371/journal.pgen.1002490 22359512
63. Bjerregaard P. Inuit Health in Transition Greenland survey 2005–2010 Population sample and survey methods [Internet]. 2011. Available: http://www.si-folkesundhed.dk/upload/inuit_health_in_transition_greenland_methods_5_2nd_revision.pdf
64. Bjerregaard P, Curtis T, Borch-Johnsen K, Mulvad G, Becker U, Andersen S, et al. Inuit health in Greenland: a population survey of life style and disease in Greenland and among Inuit living in Denmark. Int J Circumpolar Health. 2003;62 Suppl 1: 3–79. 14527126
65. Matthews DR, Hosker JP, Rudenski AS, Naylor BA, Treacher DF, Turner RC. Homeostasis model assessment: insulin resistance and beta-cell function from fasting plasma glucose and insulin concentrations in man. Diabetologia. 1985;28: 412–419. 3899825
66. Bjerregaard P, Pedersen HS, Mulvad G. The associations of a marine diet with plasma lipids, blood glucose, blood pressure and obesity among the inuit in Greenland. Eur J Clin Nutr. 2000;54: 732–737. 11002386
67. Voight BF, Kang HM, Ding J, Palmer CD, Sidore C, Chines PS, et al. The metabochip, a custom genotyping array for genetic studies of metabolic, cardiovascular, and anthropometric traits. PLoS Genet. Medical Population Genetics, The Broad Institute of Harvard and Massachusetts Institute of Technology, Cambridge, Massachusetts, United States of America.; 2012;8: e1002793. doi: 10.1371/journal.pgen.1002793;
68. Delaneau O, Zagury J-F. Data Production and Analysis in Population Genomics [Internet]. Pompanon F, Bonin A, editors. Methods in molecular biology (Clifton, N.J.). Totowa, NJ: Humana Press; 2012. doi: 10.1007/978-1-61779-870-2
69. Howie BN, Donnelly P, Marchini J. A flexible and accurate genotype imputation method for the next generation of genome-wide association studies. PLoS Genet. 2009;5: e1000529. doi: 10.1371/journal.pgen.1000529 19543373
70. Zhou X, Stephens M. Genome-wide efficient mixed-model analysis for association studies. Nat Genet. 2012;44: 821–824. doi: 10.1038/ng.2310 22706312
Štítky
Genetika Reprodukční medicínaČlánek vyšel v časopise
PLOS Genetics
2016 Číslo 6
- Souvislost haplotypu M2 genu pro annexin A5 s opakovanými reprodukčními ztrátami
- Individualizace léčby na základě farmakogenetiky při intoleranci k antipsychotikům – kazuistika
- Akutní intermitentní porfyrie
- Primární hyperoxalurie – aktuální možnosti diagnostiky a léčby
- Vliv kvality morfologie spermií na úspěšnost intrauterinní inseminace
Nejčtenější v tomto čísle
- Endogenous Mouse Dicer Is an Exclusively Cytoplasmic Protein
- Public Service by a Selfish Gene: A Domesticated Transposase Antagonizes Polycomb Function
- Genetic Links between Recombination and Speciation
- Identification of Novel Genetic Determinants of Erythrocyte Membrane Fatty Acid Composition among Greenlanders