Impact of FADS gene variation and dietary fatty acid exposure on biochemical and anthropomorphic phenotypes in a Hispanic/Latino cohort

Introduction

Most (>90%) polyunsaturated fatty acids (PUFA) are obtained from plant-derived dietary sources as 18 carbon (18C) PUFAs, such as linoleic (LA, 18:2n−6) and alpha-linolenic acid (ALA, 18:3n−3). They are enzymatically converted by alternating desaturation and elongation steps (Figure 1) in tissues, primarily in liver, to n−3 and n−6 highly unsaturated fatty acids (HUFA; containing ≥3 carbon–carbon double bonds). The two desaturase enzymes, Δ5 and Δ6 desaturase, are encoded by FADS1 and FADS2, respectively, and the two elongation enzymes, elongase 2 and elongase 5, by ELOVL2 and ELOVL5 (1–3). n−6 and n−3 HUFAs produced by this pathway include dihomo-gamma linolenic acid (DGLA, 20:3n−6), arachidonic acid (ARA, 20:4n−6), eicosapentaenoic acid (EPA, 20:5n−3) and docosahexaenoic acid (DHA, 22:6n−3). There are exceptions, but in general, ARA is converted to pro-thrombotic/pro-inflammatory metabolites and DGLA, EPA and DHA to anti-inflammatory/pro-resolving/anti-thrombotic metabolites. The HUFAs themselves and their many signaling metabolites, including oxylipins and endocannabinoids, play vital roles in regulating innate immunity, inflammation, brain development/function, thrombosis and energy metabolism.

FIGURE 1

Figure 1. PUFA + HUFA Synthetic Pathway. N−3 and n−6 PUFA+HUFA synthesis proceeds in parallel, using shared and alternating fatty acid desaturase (FADS) and elongation (ELOVL) enzymatic steps. The essential medium chain fatty acids (18 carbons, PUFAs), linoleic (LA, n−6) and alpha-linolenic (ALA, n−3), must be obtained from the mammalian diet in order to be converted to HUFAs. The central gray box denotes the enzymatic steps (FADS2, ELOVL5, FADS1) that collectively, will be referred to as Aggregate Desaturase Activity, an estimate of PUFA+HUFA metabolic efficiency: ARA/LA (n−6) and EPA/ALA (n−3), the product/essential PUFA precursor ratios.

The multiple levels of regulation of the PUFA pathway has been a topic of much investigation. At the substrate level, n−6 and n−3 PUFAs (dietary intermediates) compete throughout the synthetic pathway, and the Δ5 desaturase has recently been identified as a key rate-limiting step through the PUFA pathway (Figure 1) (4). This competition is important, given the dramatic increase in the ingestion of LA (n−6)-containing cooking oils and processed foods has resulted in 6 to 9% of daily energy coming from LA alone. Thus, the LA to ALA ratio typically entering the pathway is >10:1 (5, 6).

Until recently, the prevailing belief was that the capacity of humans to synthesize HUFAs from dietary precursors (LA and ALA) was uniform in all human populations. More recent studies indicate that there are dramatic ancestry-based frequency differences in genetic variation within the FADS gene cluster (human Chr11q), giving rise to divergences in the capacity of individuals to synthesize n−6 and n−3 HUFAs (7–10). For example, an ancestral FADS haplotype, associated with a limited capacity to synthesize HUFAs is nearly fixed in Native American and Greenland Inuit populations (11). However, this same haplotype is virtually absent and replaced by a derived FADS haplotype in African populations, and approximately 80% of African Americans and 44% of European Americans have two FADS alleles with the derived haplotype (11, 12). It has been shown that positive selection for the FADS haplotype toward the derived or ancestral haplotypes was driven by the need for ancient humans to adapt to environments that offered dramatically different levels of dietary PUFAs (11, 12).

These marked differences in FADS genetic variations in early human populations have persisted in modern populations and now challenge the “one size fits all” dietary PUFA recommendations in geographic areas such as the US where there is great population diversity. It is important to point out that the FADS cluster not only regulates levels of HUFAs and their metabolites, but this region has recently been suggested to be the largest multi-morbidity-associated gene cluster in the genome associated with over 40 clinical/disease phenotypes that are related to fat and lipid metabolism, inflammatory diseases/disorders and numerous cancers (13).

Combining a metabolomic and genome-wide associations study analyses, we demonstrated that the FADS1 step was associated with 23 HUFA-containing lipids and signaling molecules, including unesterified fatty acid, phospholipids, lyso-phospholipids and the endocannabinoid, 2-ARA-glycerol (4). Perhaps most surprising was that of all HUFA-containing complex lipids and signaling molecules examined, there was no other significant genetic association outside of the FADS cluster. Together, these studies reveal the importance of both dietary PUFA and HUFA exposure and genetic variation across the FADS cluster to the balance of n−6 and n−3 HUFAs and their metabolites that have the potential to impact human clinical/disease phenotypes.

In 2016, the Hispanic/Latino population in the US reached ~61 million, making Hispanics the largest racial/ethnic minority. Hispanic/Latinos have higher rates of obesity, poorly controlled high blood pressure, elevated circulating triglycerides (TGs), along with a higher prevalence of diabetes and nonalcoholic fatty liver disease (NAFLD) (14, 15). Importantly, self-identified Mexican Americans (MxAm) have the highest levels of Amerind (AI) ancestry (genetically related to the indigenous peoples of the Americas) among all US Hispanic/Latino populations (16), and AI-ancestry Hispanic/Latino populations are 2-fold more likely than other Hispanics to have NAFLD. We recently examined the US-based Multi-Ethnic Study of Atherosclerosis (MESA) cohort, which includes participants from Central America, South America, Mexico, Dominican Republic, Cuba and Puerto Rico, and demonstrated that individuals with high Amerind ancestry resulted in low (perhaps inadequate) levels of n−3 HUFAs, including EPA and DHA in circulating phospholipids (17).

For this study, the Latinos Combating Diabetes Study provided an opportunity to conduct a secondary evaluation of the interaction between dietary PUFA and FADS variant interactions in a largely Mexican American cohort. This parent study was a weight-loss intervention study with community and lifestyle-based interventions built on existing infrastructure in a local the Hispanic/Latino community. The goal of the parent study was to explore engagement techniques to attain a clinically meaningful impact on HbA1c, insulin metabolism, and markers of the metabolic syndrome. Using only the baseline data from the parent study, we sought to understand relationships between dietary PUFA exposure and FADS genetic variants on plasma (circulating) and RBC levels of PUFA and HUFA, as well as previously identified anthropomorphic phenotypes.

Materials and methods

Study population

The Latinos Combating Diabetes Study (ClinicalTrials.gov Identifier: NCT01831921) was designed to recruit adults (≥18 years old) who self-identified as Hispanic or Latino (and an absence of African American, Native American, Asian or White ancestry) and had evidence of pre-diabetes (Hemoglobin A1C of 5.6%–6.5%). Briefly, participants of this parent study were randomized to one of two intervention arms to evaluate intervention impact on Hemoglobin A1C and weight loss. While both interventions were conducted in a community setting, one arm was an intensive lifestyle intervention, and the control arm was usual care. Pre-existing major disease (diabetes, any type of cardiovascular disease, cancer) or pregnancy was exclusionary. This study, including informed consent documents in both English and Spanish, was reviewed and approved by the Institutional Review Board of Wake Forest University School of Medicine. All participants were required to provide written informed consent prior to taking part in any study-related data collection. For regulatory purposes, consent for the collection of biological samples was documented in a separate section of the consent form denoting that this as a sub-study. Only participants who signed this additional section in the consent document for the use of their biological samples were considered for this analysis.

Only baseline biospecimens and data were utilized in the current study. Fasting blood (8 h overnight fast) was collected in the Clinical Research Unit of the Wake Forest Baptist Medical Center. Biospecimens were stored frozen (−70°C) until analyses. Table 1 shows the demographic characteristics and blood chemistries of the 135 participants evaluated for the current study (for whom buffy coat samples were available for DNA isolation and who unequivocally self-identified as solely Hispanic/Latino).

TABLE 1

Table 1. Demographics and clinical values of the Hispanic/Latino cohort at baseline.

A European ancestry cohort of healthy adults, residing in the same area of North Carolina as the Hispanic/Latino cohort, was used only as a comparator for the frequency distribution of the FADS gene cluster variant rs174537 (Figure 2). This cohort (n = 100; 71% female) was a subset from a larger non-interventional study that was reviewed and approved by the Wake Forest Baptist Medical Center Institutional Review Board (IRB00016822). This subset cohort was ill-suited for further comparison with the Hispanic/Latino cohort as the European ancestry cohort was younger (37.1 ± 11.6; mean, SD), had a considerably lower mean BMI (23.9 ± 2.8) and did not have pre-diabetes (an exclusion criterion for that study).

FIGURE 2

Figure 2. Dietary LA and ALA Exposure. A Spanish language version of the Dietary History Questionnaire-1 (DHQ-1) was used to capture nutrients exposure based on the typical dietary patterns. DHSQ-1 data from 113 participants was analyzed. The dietary essential fatty acids exposure histograms for LA (A) and ALA (B) are shown as percent of daily energy (kcal).

Blood chemistry measurements

Fasting blood (8 h overnight fast) was collected from the Hispanic/Latino cohort participants and used for the measurements of serum lipids [Total Cholesterol, HDL, LDL, TG, VLDL(calculated)] and glucose-regulatory variables (glucose, insulin, HbA1C, HOMA-IR). These measurements were performed by LabCorp, a third party clinical testing laboratory (Burlington NC, United States). Additional tubes of blood were obtained for DNA isolation (in acid citrate dextrose) and for fatty acid analyses of plasma and red blood cells (in EDTA).

DNA isolation and SNP genotyping

Leukocyte-derived DNA was isolated from whole, frozen blood by routine molecular biology techniques (Gentra Puregene, Qiagen, Germantown MD, United States) and qualified by UV absorbance. This DNA was used to determine the rs174537 genotype employing a PCR-based TaqMan methodology. Following the manufacturer’s instructions, DNA (5 μg) was incubated in the presence of a TaqMan Universal Master Mix II (Applied Biosystems; Waltham MA, USA) and validated rs174537 primers (Applied Biosystems) in an Applied Bioscience 7,500 PCR instrument. Negative and positive controls (derived cell line-derived DNA) were employed.

Analysis of plasma and red blood cells fatty acids

Fasting baseline plasma was utilized for the analysis of circulating fatty acids by a method previously described (18), which was modified from that originally reported by Metcalf (19). Briefly, fatty acids in duplicate plasma samples (100 μL) were hydrolyzed from complex lipids under alkaline conditions in the presence of an internal standard (triheptadecanoin; 100 μg; NuChek Prep; Elysian MN, United States). Fatty acids were then converted to methyl esters in the presence of boron trifluoride (Fisher Scientific; Waltham MA, United States). The fatty acid methyl esters (FAMEs) were identified on the basis of retention times of commercially available authentic fatty acid methyl ester standards and quantified relative to the internal standard. Approximately 25 fatty acids accounted for >99% of the total fatty acids in the sample.

Total sample fatty acids were analyzed by gas chromatography with flame ionization detection using either a Hewlett Packard 5,890 or 7890B (Agilent Technologies; Santa Clara CA, United States) instruments with identical results based on their respective response factor. Each response factor was calculated using external standard (fatty acid methyl esters of 16:0, 18:0, 18:1, 18:2, 10:0, 20:5, 22:0, 22:5, 23:0, 24:0) sets for quality assurance purposes (20). A mixture of known FAMEs was run with each sample set to monitor instrument performance. Both instruments were equipped with an Agilent J&W DB-23 column (30 m, 0.25 mm ID, 0·25 μm film; Agilent) fitted with an inert pre-column (1 m, 0.53 mm ID) for cool on-column injection.

After blood collection, red blood cells (RBC) were washed with isotonic saline, suspended to approximately 50% hematocrit, aliquotted and then stored at −70°C. RBC fatty acid content was analyzed in thawed cell lysates. A total lipid extract was prepared from RBC lysates (400 μl in duplicate) using an acidified Bligh-Dyer method (21). The resultant chloroform extract was dried under a stream of nitrogen in the presence of the internal standard (50 μg triheptadecanoin). This extract was then subjected to alkaline hydrolysis and derivatization for FAME preparation and analyzed as described for plasma above. To account for differences in the RBC content of the lysates, the fatty acid data was normalized to the hemoglobin content (below) of the lysate.

The fatty acid data are presented as percent of total fatty acids derived from concentration data for plasma (mg fatty acid/dl plasma). RBC fatty acid content was normalized to hemoglobin mass (mg fatty acid/g hemoglobin).

Determination of hemoglobin concentration

The hemoglobin content of RBC lysates was quantitatively determined using a colorimetric (Drabkin) method (22, 23) following the manufacturer’s (Stanbio Laboratory; Boene TX, United States) instructions. Duplicate RBC lysate samples were diluted in the Drabkin reagent and absorbance at 540 nm was recorded. The lysate hemoglobin content (mg/mL) was calculated from a cyanmethemoglobin standard curve.

Dietary exposure assessment

The Spanish version of the Diet History Questionnaire-1 [DHQ1; (24)] was administered to study participants in a paper format at baseline. A subset of 119 participants completed the 144-item questionnaire. Sixteen participants declined to complete the dietary survey. The response data were captured by manual data entry and analyzed with the NCI Diet*Calc software program [v1.4.3; (24)] in order to extract dietary fat calories and dietary fatty acid exposures, specifically total PUFAs, essential fatty acids (LA, ALA) and HUFAs (ARA, EPA, DHA). Although this is an older food frequency questionnaire instrument, it has been validated under multiple conditions (25–28) and it is the first that could be offered in the Spanish language.

Data analysis

Statistical analyses were performed using R version 4.1.3 (29), Stata11 and GraphPad Prism5 (San Diego CA, United States). Data were checked for outliers. Descriptive statistics were calculated as mean, standard deviation (SD) and 95% confidence intervals (95% CI). Comparison of rs174537 genotype were tested linear regression using an additive model for genotype with adjustments for fatty acid intake and for multiple testing (Bonferroni). In certain cases, non-parametric Kruskal-Wallis ANOVA with Mann–Whitney U post hoc analysis was employed. When a recessive model (GG + GT vs. TT) was employed, the Students t test (two-tailed) was employed. Correlations were evaluated using Spearman analysis.

Results

An assessment of dietary fatty acid exposure in this cohort was accomplished using a validated food frequency questionnaire administered in Spanish. The hardcopy Diet History Questionnaire-1 (DHQ1; 144 items) was attempted by 119 participants at the baseline visit. Of these 119 dietary exposure questionnaires, the data from 113 participants was further evaluated for nutrient exposures. The questionnaires from 6 participants were rejected by the analysis software [Diet*Calc, (24)] based on its criteria for data completeness. The survey was used to capture nutrient estimates (typical food choices, frequency and portion size) and dietary supplement use. With regard to PUFA-related supplements, 16 participants (14%) reported the use of fish oil and two participants (1.7%) reported the use of GLA-enriched products (evening primrose oil or spirulina). The assessment of typical dietary habits showed a striking difference in dietary exposure to the essential fatty acids, LA (n−6) and ALA (n−3), expressed as % of daily energy (kcal), as is evident in Figure 2. The dietary intake pattern of the participants of the current study was not unique, as it was comparable to that of a recent, and much larger assessment of Hispanics living in the US, as shown in Table 2 (30). LA comprised 7.29% of daily energy and ~ 90% of all PUFA+HUFAs consumed in the diet. In both assessments, the ratio of LA/ALA dietary exposure had a value greater than 9. This suggests that Hispanic/Latino populations residing in the US are typically consuming the Modern Western Diet. An evaluation of the genetic variants of the FADS gene cluster in this cohort is expected to be informative since this gene cluster governs the PUFA pathway efficiency and as such will be instructive regarding how well-suited Hispanic/Latino populations are to the Modern Western Diet.

TABLE 2

Table 2. Dietary fat exposure evaluated in Hispanic/Latino cohorts.

Over the past decade, we have explored population differences in the genetic variants of the FADS gene cluster, the impact of these variants on the metabolism of HUFAs (8, 9), and the role(s) that these variants may play in the differences in disease burden for different populations (10, 17, 31). An assessment of the frequency distribution of the rs174537 SNP in this cohort (n = 135), self-identifying solely as Hispanic/Latino, was unlike any we have previously analyzed. Figure 3 shows that the frequency distribution for rs174537 in this Hispanic/Latino group had a profile nearly opposite of that typically observed for a European American ancestry cohort residing in the same geographical region. Most notably, the TT genotype predominates (53%) in the Hispanic/Latino cohort. This same genotype has a very low frequency (~10%) in populations of European ancestry (Figure 3) and is exceedingly rare in African Americans [<1%; (8)].

FIGURE 3

Figure 3. Frequency distribution of FADS Variant rs174537 in Hispanic/Latino and European American Cohorts. Leukocyte-derived DNA was genotyped at rs174537 using a TaqMan platform, as described in Methods. The frequency of the genotype distribution for the Hispanic/Latino cohort and a European American cohort living the same geographical area are shown.

Previously, we showed that the rs174537 TT genotype is associated, phenotypically, with a low efficiency of PUFA biosynthetic pathway in European American and African American cohorts (8, 9). Therefore, we sought to evaluate the impact of the FADS genotypic frequency on circulating fatty acid profiles in this Hispanic/Latino cohort. These relationships have typically been assessed in the fasting plasma fatty acid pool. However, the availability of matched red blood cells samples afforded a novel opportunity to apply this evaluation in a second easily accessible tissue and in parallel to plasma. Figure 4 shows the n−6 PUFA levels, stratified by rs174537 genotype and adjusted for fatty acid intake, for plasma (Figure 4A) and for RBC (Figure 4B). It is notable that the LA content in both tissues was elevated, which reflects the consumption of a typical LA-enriched Modern Western Diet. Only the circulating (plasma) level of ARA was modestly impacted by dietary intake of ARA after adjustment for multiple testing (p = 0.001; Bonferroni correction, p = 0.0029). Nevertheless, each of the LA-derived n−6 HUFAs (GLA → DGLA → ARA) products in the plasma of this Hispanic/Latino cohort demonstrated a clear genotypic-dependence. This pattern is analogous to that seen in European and African ancestry cohorts with highest plasma ARA levels in the GG genotype and the lowest in the TT group (Figure 4A).

FIGURE 4

Figure 4. Impact of rs174537 on n−6 PUFA + HUFAs in Plasma and RBC in a Hispanic/Latino cohort. Total plasma (A, top row) and RBC (B, bottom row) fatty acid were analyzed as described in Methods and expressed as percent of total fatty acid in the sample. Individual data are shown, stratified by rs174537 genotype (○, GG; □, GT; Δ, TT) for each n−6 PUFA along the PUFA+HUFA pathway. The mean and 95% CI is shown for each genotype. The data were analyzed by linear regression using an additive model for genotype, adjusting for fatty acid intake and multiple testing (Bonferroni correction, p = 0.0029). The resultant p-value for impact of genotype on each fatty acid stated under each plot.

Similar to the observation for most plasma fatty acid levels, no significant impact of dietary fatty acid intake (as assessed by dietary survey) was seen for RBC fatty acid levels. However, in contrast to plasma, a genotypic effect on RBC n−6 PUFA+HUFA content was not generally evident (Figure 4B). In RBC, there was a trend for a genotypic impact for ARA content, but the wide variation in the content of this HUFA likely reflects a longer-term variation in dietary exposure to ARA and to a lesser extent the impact of hepatic PUFA metabolism. When a recessive model (GG + GT vs. TT; data not shown) was employed, then a genotypic impact for RBC DGLA (two-tailed t test; p = 0.01) and RBC ARA (two-tailed t test; p = 0.04) were observed. It is noteworthy that the PUFA+HUFA content of RBC was on average 25% lower than that of plasma. The lower RBC PUFA+HUFA content was compensated for a larger saturated fatty component compared to plasma.

The levels of n−3 PUFA+HUFAs were, on average, 14- and 8-fold lower than the n−6 PUFA+HUFA content of plasma and RBC, respectively (Supplementary Figure 1A). No impact of dietary fatty acid intake was observed to affect n−3 PUFA or HUFAs after adjustment for multiple testing. However, the levels of two plasma n−3 HUFAs were impacted by rs174537 genotype (EPA, p = 0.001; DPA, p = 0.0001). EPA and ARA are both products of FADS1 in the PUFA synthetic pathway (Figure 1). DPA is the elongation product of EPA and parallel to adrenic acid (ADR, 22:5 n−6) on the n−6 side of the PUFA pathway. Collectively, the plasma PUFA+HUFA patterns in this rs174537 TT-enriched Hispanic cohort align with those previously observed for European and African ancestry cohorts. In contrast, neither RBC n−6 HUFAs nor n−3 HUFA profiles were significantly impacted by rs174537genotype (Supplementary Figure 1B) or dietary intake.

Fatty acid product-precursor ratios are typically utilized as a surrogate for PUFA pathway enzymatic steps. Figure 5 shows the ratios for the Aggregate Desaturase Activity (FADS2 → ELOVL2 → FADS1; Figure 1 central gray box), which represents the enzymatic step common to the n−6 (ARA/LA) and n−3 (EPA/ALA) sides of the PUFA pathway utilizing the same enzymes to service both substrate types. The plasma ratios for both the n−6 and n−3 Aggregate Desaturase Activities, were robustly differentiated by rs174537 genotype (p = 0.0001, p = 0.0001, respectively; Figures 5A,B; Supplementary Table 1). This is consistent with the gradient of high (GG) to low (TT) metabolic efficiency of the PUFA pathway governed by this genotypic variant (8, 9).

FIGURE 5

Figure 5. Estimation of the Aggregate Desaturase Activity in the PUFA + HUFA Pathway for Plasma and RBC in a Hispanic/Latino Cohort. The Aggregate Desaturase Activity (FADS2 → ELOVL2 → FADS1) of both the n−6 and n−3 sides of the PUFA + HUFA pathway was estimated by the product—precursor ratio surrogates of ARA/LA (n−6) and EPA/ALA (n−3). The ratios for plasma (A, n−6; B, n−3; top row) and RBC (C, n−6; D, n−3; bottom row) are shown as individual ratio data, stratified by rs174537 genotype(○, GG; □, GT; Δ, TT). The mean and 95% CI is shown for each genotype. The data were analyzed by ANOVA with the resultant p-value for each fatty acid ratio stated above each plot.

In contrast, the hardiness of the genotypic governance of the RBC aggregate ratios were more modest for n−6 (ARA/LA, p = 0.003; Figure 5C) than that for plasma and was absent for the RBC n−3 EPA/ALA ratio (p = 0.131; Figure 5D). In order to determine the enzymatic step(s) that contributed to this genotypic governance, a ratio surrogate was calculated for each of the FADS and ELVOL steps within the Aggregate Desaturase and beyond (Supplementary Table 1). The abundance of the n−6 PUFA+HUFA lent itself to assessment of each PUFA pathway step. Plasma estimates of FADS2 (GLA/LA), ELVOL5 (DGLA/GLA) and FADS1 (ARA/DGLA) activities are all strongly associated with rs174537 genotype. It was not possible to quantitate plasma SDA or ETA (eicosatetraenoic acid, 20:4n−3) levels across this cohort due to their low abundance. Thus, it is not possible to confirm, unequivocally, that these three enzymatic steps are similarly genetically governed on the n−3 PUFA side of the pathway. Interestingly, the impact of rs174537 on plasma Δ6 FADS2 activity (GLA/LA, n−6; p < 0.0001) was much more robust than that for its plasma Δ4 activity (DHA/DPA; n−3, p = 0.0698).

RBC have recently been observed to express very low levels of FADS1/2 protein and activity, which may function as a compensatory mechanism to maintain RBC integrity (32). The product/precursor activity surrogate, ARA/DGLA, for FADS1 step explains the majority (p < 0.0001) of the genotypic impact on the n−6 PUFA Aggregate Desaturase Activities in this tissue (Figure 5C; Supplementary Table 1). It remains to be determined the extent to which FADS1 activity (hepatic vs. RBC) predominates in defining the RBC HUFA profile. ELOVL2 (DGLA/GLA) activity provides a small (p = 0.0107), but significant contribution to the RBC n−6 PUFA pathway. A genotypic impact on n−3 PUFA activity (EPA/ALA; p = 0.0621) was absent when estimated from RBC PUFA levels.

Beyond the core Aggregate Desaturase Activity surrogate, only the subsequent ELOVL2/5 step (ADR/ARA; p = 0.001; Supplementary Table 1) on the n−6 PUFA side showed a contribution to genetic governance by the FADS variant. The plasma DHA/EPA (not shown), which utilizes both elongation and desaturase steps (Figure 1), was modestly impacted (ANOVA, 10.52 (2), p = 0.0052) by rs174537. This examination of the PUFA pathway surrogates in blood components suggests that the plasma PUFA profile is under more immediate influence of hepatic PUFA pathway efficiency regulation than that for the RBC PUFA profiles, which are thought to be a reflection of longer-term diet exposure.

Given the relatively low abundance of n−3 PUFA+HUFA blood levels in this cohort, an evaluation of the Omega-3 Index (sum of RBC EPA + DHA) was of particular interest. The Omega-3 Index has been proposed to be a useful marker of n−3 HUFA status (33, 34) and a potential disease risk factor (35, 36). Although the Omega-3 Index did not differ by rs174537 genotype (Supplementary Figure 2), the mean cohort value was 2.89 (±1.65, SD), which is in the range (<4%) that is suggested to be the least protective against cardiovascular disease (36).

Beyond the biochemical phenotypes, clinical phenotypes have also been evaluated in the context of FADS gene variants. Only height was significantly impacted by rs174537 genotype. Across the entire cohort (n = 135), participants carrying the GG genotype were, on average, 6.5 cm taller (p = 0.0036, Supplementary Table 2) than those having the TT genotype at rs174537. However, a sex difference was observed since males (n = 22) were significantly taller than females (n = 113; p < 0.0001, two-tailed t test) so height was examined more closely. In order to eliminate the influence of sex, the impact of FADS genotype on height was re-evaluated in the female subset. Figure 6A shows that the female TT participants (n = 59) were clearly of shorter stature (on average a 4.5 cm difference; p = 0.0001, two-tailed t-test) than those carrying the G allele (GG + GT; n = 54). Although body weight would be expected to be impacted more directly by lifestyle choices, female body weight was also significantly lower in TT subjects compared to those carrying the G allele (TT, 76.9 ± 1.6 kg [73.8–79.9, 95%CI]; GG + GT, 82.2 ± 1.7 kg [78.8–85.6, 95%CI]; p = 0.0216, two-tailed t test). Finally, a positive relationship between height and the n−6 Aggregate Desaturase Activity surrogate (ARA/LA) was observed (Figure 6B) for the female subset of the cohort (Spearman r, 0.2652; 95%CI 0.0791 to 0.4335; p = 0.0045). Interestingly, when taking into account each individual PUFA pathway step, the association with height was more robust for FADS1 activity (ARA/DGLA; Spearman r, 0.2427; 95%CI: 0.0552 to 0.4137; p = 0.0096) than for FADS2 activity (GLA/LA; Spearman r, 0.17152; 95%CI: −0.0192 to 0.502; p = 0.069).

FIGURE 6

Figure 6. Potential Interactions Among rs174537, Female Height and Aggregate Desaturase Activity in a Hispanic/Latino Cohort. (A) Two female subsets of the cohort were created based on the presence of the rs174537 G allele (GG + GT, ○, n = 54; and TT, Δ n = 59) in order to evaluate a genotypic impact on participant height. The data were analyzed by two-tailed t test with the resultant p-value for the subset comparison stated under the plot. (B) An association of female height (n = 113) with n−6 Aggregate Desaturase Activity (ARA/LA) was evaluated, Spearman r = 0.2652, p = 0.0045.

Discussion

The frequency distribution of rs174537 in the pre-diabetic Hispanic cohort was found to be strikingly different from that of European or African ancestry cohorts, but quite similar to that observed for a Mexican-dwelling cohort (37) and other populations of South and Central American ancestries (17, 38, 39). The rs174537 variant is in high linkage disequilibrium (LD) with a number of other FADS SNPs (38) that have been evaluated in Hispanic populations and would be expected to behave similarly. These ancestry-related differences in FADS variants, which govern PUFA+HUFA biochemical metabolic phenotypes, give rise to important considerations concerning gene-diet interactions that may contribute to health disparities.

Our analyses of the genotypic-dependence of n−3 and n−6 PUFA+HUFAs showed that circulating levels of desaturase step products and precursors were significantly associated with rs174537 genotype. This was particularly true for plasma and to a lesser extent for RBC (Figure 4; Supplementary Figure 1). This suggests that fasting circulating plasma fatty acids are more reflective of hepatic HUFA synthetic activity than are those of the RBC fatty acid profile. RBC are continuously produced as aged cells are removed. The mixture of circulating cells of heterogeneous ages, and the significantly lower RBC PUFA+HUFA content compared to plasma (RBC, 35.4 ± 11.9% total, mean ± SD; plasma 47.3 ± 4.2%total; two-tail t-test, p < 0.0001), may both contribute to placing the RBC fatty acid profile more distal from tissue sources that experience genetic governance of PUFA+HUFA synthesis compared to plasma. Nevertheless, the RBC Aggregate Desaturase Activity (ARA/LA, Figure 5C) was observed to be robustly associated with rs174537 genotype, suggesting that the DGLA and ARA content of RBC membranes may remain relatively stable over the cell lifetime. This feature may add to the value of the use of RBC fatty acid levels as indicators of longer-term n−6 HUFA dietary exposure.

As previously seen in cohorts of European and African ancestries, we observed that the GG genotype was associated with the highest plasma ARA and EPA (both FADS1 products) levels (Figure 4A) and the highest plasma Aggregate Desaturase Activity (ARA/LA and EPA/ALA, Figures 5A,B) in the Hispanic cohort. These biochemical phenotypes are indicative of the most efficient HUFA synthetic efficiency. However, more than 50% of the Hispanic cohort carries the TT genotype, which is associated with the biochemical phenotype of the least efficient HPUFA synthetic activity. This low efficiency HUFA synthetic phenotype, in the presence of excess dietary n−6 PUFAs (Figure 3; Table 2), typical of the Modern Western Diet (40), has the potential to create a n−3 HUFA deficiency and thus a gene-diet misalignment. Indeed, the assessment of the dietary fat exposure of this cohort suggests that their diet is reflective of the Modern Western Diet, with 7.29% of energy (17.9 ± 0.1 g/day, mean ± DEM) from LA and an LA/ALA ratio of 9.3 (Table 2). This dietary exposure and the low abundance of circulating n−3 PUFA+HUFA species likely contributed to our inability to unequivocally demonstrate a genetic governance of plasma estimates of n−3 pathway activities FADS2 (SDA/ALA), ELVOL5 (ETA/SDA) and FADS1 (EPA/ETA). Clearly, a supplementation study in a Hispanic cohort aimed at dietary ALA-enrichment and lowering of LA intake would be beneficial for better understanding the complexities of PUFA+HUFA metabolism in a low-efficiency HUFA metabolizing cohort. For example, a clinical trial prospectively comparing the capacity of a Hispanic cohort, stratified by rs174537 genotype, to metabolize an ALA enrich diet by reducing levels of LA and increasing ALA to a ratio closer to 2:1 (compared to the present 8:1) would provide great insight into the metabolic capacity of the low-efficiency metabolizing genotype. The expected result from such a trial would be a reduction in ARA levels and ARA-derived signaling oxylipins. In contrast, we postulate such a diet would elevate EPA and DHA levels and their oxylipin metabolites thereby balancing n−6 to n−3 HUFA and oxylipin levels. As with any study, there were limitations that warrant mentioning. This was a sub study of a larger trial and the parent study was not powered specifically to assess these biomarkers. The dietary questionnaire required participants to self-report the foods they consumed and therefore could have resulted in biased reporting of their dietary intake.

Beyond the biochemical phenotypes rooted in HUFA synthesis, FADS variants have been associated with clinical and disease phenotypes. A recent large study of healthy college-age Mexican students (37) reported that TG and VLDL levels were significantly higher for the TT genotype (48% frequency distribution) of rs174546, which is in LD with rs174537. A similar association between TG and rs174537 TT was observed in a large Hispanic cohort, with varying degrees of Amerind ancestry (17). The absence of this association (p = 0.083, GG + GT vs. TT) in the current study may arise from differences in health status (pre-diabetic) and the age range (18–66 years) of participants as well as a small cohort size. The former factors may also have contributed to an overall increase in blood lipid values with age.

The relationship among of PUFA+HUFA metabolism, FADS gene variants and anthropomorphic indices is an interesting one. Despite the relatively small cohort size of the current study, we observed a significant, negative association between the rs174537 T allele and height (Figure 6A) and a positive association with height and the n−6 Aggregate Desaturase Activity (ARA/LA; Figure 6B) in the female subset of the cohort. Indeed, the FADS locus appears to be an important site under selective and evolutionary pressures for adaptations to the presence of diet-derived nutrients that are critical for normal development and immune system integrity (11, 12, 17, 41, 42).

Overall, our findings strongly suggest a gene-diet misalignment for Hispanic/Latino populations consuming a Modern Western Diet. As a result, levels and ratios of n−6 and n−3 PUFA+HUFAs become out of balance leading to deficiencies in n−3 LC-PUFA+HUFAs and their metabolites, which have anti-inflammatory, cardio-protective, and triglycerides lowering properties (43–45). Even though the abundance of dietary LA exposure in this cohort (Table 2) is comparable to that assessed in other US-dwelling Hispanic/Latino cohorts (30, 46 ¹), it is significantly lower than that of non-Hispanic Blacks and Whites (46). Nevertheless, the burden of metabolic diseases remains higher in Hispanic/Latinos populations (14, 15) and thus provides an underappreciated gene*diet feature for health disparities.

Thus, the Hispanic/Latino population, with its high frequency of the TT genotype at rs174537, is biochemically less able to synthesize the necessary HUFAs and this phenotype is exacerbated by the high LA content of the diet. The resultant n−3 HUFA deficiency, imposed by the gene*diet interaction, might readily be rendered less impactful if health care professionals provided nutritional guidance and education encouraging the consumption of n−3 HUFA-containing foods or dietary supplements. Such personalized nutrition approaches are thus potentially and pragmatically important methods to reduce health disparities.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding author.

Ethics statement

The studies involving human participants were reviewed and approved by Wake Forest Baptist Medical Center Institutional Review Board. The patients/participants provided their written informed consent to participate in this study.

Author contributions

BK and SS performed the fatty acid analyses. SS isolated DNA and performed rs174537 genotyping and wrote the manuscript. JL entered and analyzed the dietary history data. SS and JL performed the data analyses. CY was the study coordinator for the intervention (parent) study. FC and MS designed the current study. SS, MS, FC, and MV edited the manuscript. MV obtained funding for the intervention (parent) study. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by grants from the National Institutes of Health, P01 MD006917 (MV) and R01 AT008621 (FC).

Acknowledgments

The authors gratefully acknowledge use of the services and facilities of the Clinical Research Unit, funded by the National Center for Advancing Translational Sciences (NCATS), National Institutes of Health, through Grant Award Number UL1TR001420. The authors would like to thank Dr. Brian Hallmark (Center for Biomedical Informatics and Biostatistics, BIO5 Institute, University of Arizona, Tucson, AZ) for generously providing statistical advice and guidance.

Conflict of interest

FC is a co-founder of a start-up company TyrianOmega, which focuses on the production of omega-3 PUFAs by cyanobacteria, largely for animal feeds and aquaculture. He is also the co-founder of Resonance Pharma, which focuses on determining clinical level of and inhibitors that block secreted phospholipase A2 isoforms.

The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnut.2023.1111624/full#supplementary-material

Abbreviations

ADR, Adrenic acid; AI, Amerind; ALA, Alpha-linolenic acid; ARA, Arachidonic acid; DGLA, Dihsomo-gamma linolenic acid; DHA, Docosahexaenoic acid; DHQ1, Dietary History Questionnaire-1; DPA, Docosapentaenoic acid; ELOVL, Elongase; EPA, Eicospentaenoic acid; ETA, Eicosatetroenoic acid; FADS, Fatty acid desaturase; FAME, Fatty acid methyl ester; GLA, Gamma linolenic acid; HDL, High density lipoprotein; HOMA-IR, Homeostatic Model Assessment of Insulin Resistance; HUFA, Highly unsaturated fatty acid; LA, Linoleic acid; LD, Linkage disequilibrium; LDL, Low density lipoprotein; NAFLD, Nonalcoholic fatty liver disease; PUFA, Polyunsaturated fatty acids; RBC, SNP, Single nucleotide polymorphism; TG, Triglycerides; VLDL, Very low density lipoprotein.

Footnotes

^Reference 46 has been submitted and is undergoing revisions prior to acceptance for publication.

References

1. Lands, WEM, Libelt, B, Morris, A, Kramer, NC, Prewitt, TE, Bowen, P, et al. Maintenance of lower proportions of (n − 6) eicosanoid precursors in phospholipids of human plasma in response to added dietary (n − 3) fatty acids. Biochim Biophys Acta. (1992) 1180:147–62. doi: 10.1016/0925-4439(92)90063-S

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Sprecher, H. Metabolism of highly unsaturated n-3 and n-6 fatty acids. Biochim Biophys Acta Mol Cell Biol Lipids. (2000) 1486:219–31. doi: 10.1016/S1388-1981(00)00077-9

CrossRef Full Text | Google Scholar

3. Brenna, JT, and Kothapalli, KSD. New understandings of the pathway of long-chain polyunsaturated fatty acid biosynthesis. Curr Op Clin Nutr Metab Care. (2021) doi: 10.1097/MCO.0000000000000810

CrossRef Full Text | Google Scholar

4. Reynolds, LM, Dutta, R, Seeds, MC, Lake, KN, Hallmark, B, Mathias, RA, et al. FADS genetic and metabolomic analyses identify the ∆5 desaturase (FADS1) step as a critical control point in the formation of biologically important lipids. Sci Rep. (2020) 10:15873. doi: 10.1038/s41598-020-71948-1

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Blasbalg, TL, Hibbeln, JR, Ramsden, CE, Majchrzak, SF, and Rawlings, RR. Changes in consumption of omega-3 and omega-6 fatty acids in the United States during the 20th century. Am J Clin Nutr. (2011) 93:950–62. doi: 10.3945/ajcn.110.006643

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Simopoulos, AP. The importance of the omega-6/omega-3 fatty acid ratio in cardiovascular disease and other chronic diseases. Exp Biol Med. (2008) 233:674–88. doi: 10.3181/0711-MR-311

PubMed Abstract | CrossRef Full Text | Google Scholar

7. 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 [Internet]. (2009) 5, 5:e1000338. doi: 10.1371/journal.pgen.1000338

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Mathias, R, Sergeant, S, Ruczinski, I, Torgerson, D, Hugenschmidt, C, Kubala, M, et al. The impact of FADS genetic variants on omega-6 polyunsaturated fatty acid metabolism in African Americans. BMC Genet. (2011) 12:50. doi: 10.1186/1471-2156-12-50

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Sergeant, S, Hugenschmidt, CE, Rudock, ME, Ziegler, JT, Ivester, P, Ainsworth, HC, et al. Differences in arachidonic acid levels and fatty acid desaturase (FADS) gene variants in African Americans and European Americans with diabetes or the metabolic syndrome. Br J Nutr. (2012) 107:547–55. doi: 10.1017/S0007114511003230

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Chilton, FH, Manichaikul, A, Yang, C, O'Connor, TD, Johnstone, LM, Blomquist, S, et al. Interpreting clinical trials with Omega-3 supplements in the context of ancestry and FADS genetic variation. Front Nutr. (2022) 8:8. doi: 10.3389/fnut.2021.808054

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Harris, DN, Ruczinski, I, Yanek, LR, Becker, LC, Becker, DM, Guio, H, et al. Evolution of hominin polyunsaturated fatty acid metabolism: from Africa to the New World. Genome Biol Evol. (2019) 11:1417–30. doi: 10.1093/gbe/evz071

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Mathias, RA, Fu, W, Akey, JM, Ainsworth, HC, Torgerson, DG, Ruczinski, I, et al. Adaptive evolution of the FADS gene cluster within Africa. PLoS One. (2012) 7:e44926. doi: 10.1371/journal.pone.0044926

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Fadason, T, Schierding, W, Lumley, T, and O’Sullivan, JM. Chromatin interactions and expression quantitative trait loci reveal genetic drivers of multimorbidities. Nat Commun. (2018) 9:5198. doi: 10.1038/s41467-018-07692-y

PubMed Abstract | CrossRef Full Text | Google Scholar

14. National Diabetes Statistics Report, (2017). Estimates of diabetes and its burden in the United States [internet]. Department of Health and Human Services. 2017. Available at: https://dev.diabetes.org/sites/default/files/2019-06/cdc-statistics-report-2017.pdf (Accessed August 7, 2022).

Google Scholar

15. Saab, S, Manne, V, Nieto, J, Schwimmer, JB, and Chalasani, NP. Nonalcoholic fatty liver disease in Latinos. Clin Gastroent Hepatol. (2016) 14:5–12. doi: 10.1016/j.cgh.2015.05.001

CrossRef Full Text | Google Scholar

16. Bryc, K, Durand, EY, Macpherson, JM, Reich, D, and Mountain, JL. The genetic ancestry of African Americans, Latinos, and European Americans across the United States. Am J Human Genet. (2015) 96:37–53. doi: 10.1016/j.ajhg.2014.11.010

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Yang, C, Hallmark, B, Chai, JC, O'Connor, TD, Reynolds, LM, Wood, AC, et al. Impact of Amerind ancestry and FADS genetic variation on omega-3 deficiency and cardiometabolic traits in Hispanic populations. Commun Biol. (2021) 4:918. doi: 10.1038/s42003-021-02431-4

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Sergeant, S, Ruczinski, I, Ivester, P, Lee, TC, Morgan, TM, Nicklas, BJ, et al. Impact of methods used to express levels of circulating fatty acids on the degree and direction of associations with blood lipids in humans. Br J Nutr. (2016) 115:251–61. doi: 10.1017/S0007114515004341

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Metcalfe, LD, Schmitz, AA, and Pelka, JR. Rapid preparation of fatty acid esters from lipids for gas chromatographic analysis. Anal Chem. (1966) 38:514–5. doi: 10.1021/ac60235a044

CrossRef Full Text | Google Scholar

20. Brenna, JT, Plourde, M, Stark, KD, Jones, PJ, and Lin, Y-H. Best practices for the design, laboratory analysis, and reporting of trials involving fatty acids. Am J Clin Nutr. (2018) 108:211–27. doi: 10.1093/ajcn/nqy089

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Bligh, EG, and Dyer, WJ. A rapid method of total lipid extraction and purification. Can J Biochem Physiol. (1959) 37:911–7. doi: 10.1139/y59-099

CrossRef Full Text | Google Scholar

22. Drabkin, DL, and Austin, JH. Spectrophotometric studies: 1. Spectrophotometric constants for common hemoglobin derivatives in human, dog, and rabbit blood. J Biol Chem. (1932) 98:719–33. doi: 10.1016/S0021-9258(18)76122-X

CrossRef Full Text | Google Scholar

23. Balasubramaniam, P, and Malathi, A. Comparative study of hemoglobin estimated by Drabkin's and Sahli's methods. J Postgrad Med. (1992) 38:8–9.

PubMed Abstract | Google Scholar

24. National Institutes of Health ARP, National Cancer Institute (2007). Diet history questionnaire, Version 1.0. Available at: https://epi.grants.cancer.gov/dhq/.

Google Scholar

25. Subar, AF, Thompson, FE, Kipnis, V, Midthune, D, Hurwitz, P, McNutt, S, et al. Comparative validation of the block, Willett, and National Cancer Institute food frequency questionnaires: the eating at America's Table study. Am J Epidemiol. (2001) 154:1089–99. doi: 10.1093/aje/154.12.1089

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Thompson, FE, Subar, AF, Brown, CC, Smith, AF, Sharbaugh, CO, Jobe, JB, et al. Cognitive research enhances accuracy of food frequency questionnaire reports: results of an experimental validation study. J Am Diet Assoc. (2002) 102:212–25. doi: 10.1016/S0002-8223(02)90050-7

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Subar, AF, Kipnis, V, Troiano, RP, Midthune, D, Schoeller, DA, Bingham, S, et al. Using intake biomarkers to evaluate the extent of dietary misreporting in a large sample of adults: the OPEN study. Am J Epidemiol. (2003) 158:1–13. doi: 10.1093/aje/kwg092

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Kipnis, V, Subar, AF, Midthune, D, Freedman, LS, Ballard-Barbash, R, Troiano, RP, et al. Structure of dietary measurement error: results of the OPEN biomarker study. Am J Epidemiol. (2003) 158:14–21. doi: 10.1093/aje/kwg091

PubMed Abstract | CrossRef Full Text | Google Scholar

29. R Core Team. R: a language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing (2018).

Google Scholar

30. USDA, Agricultural Research Service. (2021). Usual intake from foods and beverages among Hispanic individuals in the US, what we eat in America, NHANES 2015–2018. Available at: https://www.ars.usda.gov/ARSUserFiles/80400530/pdf/usual/Usual_Intake_Hispanic_WWEIA_2015_2018.pdf (Accessed July 21, 2022).

Google Scholar

31. Chilton, FH, Murphy, RC, Wilson, BA, Sergeant, S, Ainsworth, HC, Seeds, MC, et al. Diet-gene interactions and PUFA metabolism: a potential contributor to health disparities and human diseases. Nutrients. (2014) 6:1993–2022. doi: 10.3390/nu6051993

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Thomas, T, Cendali, F, Fu, X, Gamboni, F, Morrison, EJ, Beirne, J, et al. Fatty acid desaturase activity in mature red blood cells and implications for blood storage quality. Transfusion. (2021) 61:1867–83. doi: 10.1111/trf.16402

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Harris, WS. The omega-3 index: from biomarker to risk marker to risk factor. Curr Atheroscler Rep. (2009) 11:411–7. doi: 10.1007/s11883-009-0062-2

CrossRef Full Text | Google Scholar

34. Harris, WS. The Omega-6:Omega-3 ratio: a critical appraisal and possible successor. Prostaglandins Leukot Essent Fatty Acids. (2018) 132:34–40. doi: 10.1016/j.plefa.2018.03.003

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Anzalone, A, Carbuhn, A, Jones, L, Gallop, A, Smith, A, Johnson, P, et al. The Omega-3 index in National Collegiate Athletic Association Division I Collegiate Football Athletes. J Athl Train. (2019) 54:7–11. doi: 10.4085/1062-6050-387-18

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Harris, WS, and von Schacky, C. The Omega-3 index: a new risk factor for death from coronary heart disease? Prev Med. (2004) 39:212–20. doi: 10.1016/j.ypmed.2004.02.030

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Vazquez-Vidal, I, Voruganti, VS, Hannon, BA, Andrade, FCD, Aradillas-García, C, Nakamura, MT, et al. Serum lipid concentrations and FADS genetic variants in young Mexican college students: the UP-AMIGOS cohort study. Lifestyle Genom. (2018) 11:40–8. doi: 10.1159/000488085

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Ensembl.org. rs174537 SNP. (2022). Available at: http://useast.ensembl.org/Homo_sapiens/Variation/Citations?db=core;r=11:61784708-61785708;v=rs174537;vdb=variation;vf=164475884 (Accessed July 20, 2022).

Google Scholar

39. Cunningham, F, Allen, JE, Allen, J, Alvarez-Jarreta, J, Amode, MR, Armean Irina, M, et al. Ensembl 2022. Nucleic Acids Res. (2021) 50:D988–95. doi: 10.1093/nar/gkab1049

CrossRef Full Text | Google Scholar

40. Cordain, L, Eaton, SB, Sebastian, A, Mann, N, Lindeberg, S, Watkins, BA, et al. Origins and evolution of the Western diet: health implications for the 21st century. Am J Clin Nutr. (2005) 81:341–54. doi: 10.1093/ajcn.81.2.341

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Fumagalli, M, Moltke, I, Grarup, N, Racimo, F, Bjerregaard, P, Jørgensen, ME, et al. Greenlandic Inuit show genetic signatures of diet and climate adaptation. Science. (2015) 349:1343–7. doi: 10.1126/science.aab2319

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Tucci, S, Vohr, SH, McCoy, RC, Vernot, B, Robinson, MR, Barbieri, C, et al. Evolutionary history and adaptation of a human pygmy population of Flores Island. Indonesia Sci. (2018) 361:511–6. doi: 10.1126/science.aar8486

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Beuther, DA, Weiss, ST, and Sutherland, ER. Obesity and asthma. Am J Respir Crit Care Med. (2006) 174:112–9. doi: 10.1164/rccm.200602-231PP

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Naderali, EK, Ratcliffe, SH, and Dale, MC. Review: obesity and Alzheimer’s disease: a link between body weight and cognitive function in old age. Am J Alzheim Dis Other Demen. (2009) 24:445–9. doi: 10.1177/1533317509348208

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Leveille, SG, Wee, CC, and Iezzone, LI. Trends in obesity and arthritis among baby boomers and their predecessors, 1971–2002. Am J Public Health. (2005) 95:1607–13. doi: 10.2105/AJPH.2004.060418

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Momin, SR, Senn, MK, Manichaikul, A, Yang, C, Mathias, RA, Phan, M, et al. Dietary sources of linoleic acid (LA) differ by race/ethnicity in adults participating in the National Health and nutrition examination survey (NHANES) between 2017–2018. Acad Nutr Diet. (forthcoming)

Google Scholar