Genome-Wide Association Study Singles Out SCD and LEPR as the Two Main Loci Influencing Intramuscular Fat Content and Fatty Acid Composition in Duroc Pigs

Intramuscular fat (IMF) content and fatty acid composition affect the organoleptic quality and nutritional value of pork. A genome-wide association study was performed on 138 Duroc pigs genotyped with a 60k SNP chip to detect biologically relevant genomic variants influencing fat content and composition. Despite the limited sample size, the genome-wide association study was powerful enough to detect the association between fatty acid composition and a known haplotypic variant in SCD (SSC14) and to reveal an association of IMF and fatty acid composition in the LEPR region (SSC6). The association of LEPR was later validated with an independent set of 853 pigs using a candidate quantitative trait nucleotide. The SCD gene is responsible for the biosynthesis of oleic acid (C18:1) from stearic acid. This locus affected the stearic to oleic desaturation index (C18:1/C18:0), C18:1, and saturated (SFA) and monounsaturated (MUFA) fatty acids content. These effects were consistently detected in gluteus medius, longissimus dorsi, and subcutaneous fat. The association of LEPR with fatty acid composition was detected only in muscle and was, at least in part, a consequence of its effect on IMF content, with increased IMF resulting in more SFA, less polyunsaturated fatty acids (PUFA), and greater SFA/PUFA ratio. Marker substitution effects estimated with a subset of 65 animals were used to predict the genomic estimated breeding values of 70 animals born 7 years later. Although predictions with the whole SNP chip information were in relatively high correlation with observed SFA, MUFA, and C18:1/C18:0 (0.48–0.60), IMF content and composition were in general better predicted by using only SNPs at the SCD and LEPR loci, in which case the correlation between predicted and observed values was in the range of 0.36 to 0.54 for all traits. Results indicate that markers in the SCD and LEPR genes can be useful to select for optimum fatty acid profiles of pork.     

New insight into the SSC8 genetic determination of fatty acid composition in pigs

Background

Fat content and fatty acid composition in swine are becoming increasingly studied because of their effect on sensory and nutritional quality of meat. A QTL (quantitative trait locus) for fatty acid composition in backfat was previously detected on porcine chromosome 8 (SSC8) in an Iberian x Landrace F₂ intercross. More recently, a genome-wide association study detected the same genomic region for muscle fatty acid composition in an Iberian x Landrace backcross population. ELOVL6, a strong positional candidate gene for this QTL, contains a polymorphism in its promoter region (ELOVL6:c.-533C < T), which is associated with percentage of palmitic and palmitoleic acids in muscle and adipose tissues. Here, a combination of single-marker association and the haplotype-based approach was used to analyze backfat fatty acid composition in 470 animals of an Iberian x Landrace F₂ intercross genotyped with 144 SNPs (single nucleotide polymorphisms) distributed along SSC8.

Results

Two trait-associated SNP regions were identified at 93 Mb and 119 Mb on SSC8. The strongest statistical signals of both regions were observed for palmitoleic acid (C16:1(n-7)) content and C18:0/C16:0 and C18:1(n-7)/C16:1(n-7) elongation ratios. MAML3 and SETD7 are positional candidate genes in the 93 Mb region and two novel microsatellites in MAML3 and nine SNPs in SETD7 were identified. No significant association for the MAML3 microsatellite genotypes was detected. The SETD7:c.700G > T SNP, although statistically significant, was not the strongest signal in this region. In addition, the expression of MAML3 and SETD7 in liver and adipose tissue varied among animals, but no association was detected with the polymorphisms in these genes. In the 119 Mb region, the ELOVL6:c.-533C > T polymorphism showed a strong association with percentage of palmitic and palmitoleic fatty acids and elongation ratios in backfat.

Conclusions

Our results suggest that the polymorphisms studied in MAML3 and SETD7 are not the causal mutations for the QTL in the 93 Mb region. However, the results for ELOVL6 support the hypothesis that the ELOVL6:c.-533C > T polymorphism has a pleiotropic effect on backfat and intramuscular fatty acid composition and that it has a role in the determination of the QTL in the 119 Mb region.     

Oligoclonal Band Status in Scandinavian Multiple Sclerosis Patients Is Associated with Specific Genetic Risk Alleles

The presence of oligoclonal bands (OCB) in cerebrospinal fluid (CSF) is a typical finding in multiple sclerosis (MS). We applied data from Norwegian, Swedish and Danish (i.e. Scandinavian) MS patients from a genome-wide association study (GWAS) to search for genetic differences in MS relating to OCB status. GWAS data was compared in 1367 OCB positive and 161 OCB negative Scandinavian MS patients, and nine of the most associated SNPs were genotyped for replication in 3403 Scandinavian MS patients. HLA-DRB1 genotypes were analyzed in a subset of the OCB positive (n = 2781) and OCB negative (n = 292) MS patients and compared to 890 healthy controls. Results from the genome-wide analyses showed that single nucleotide polymorphisms (SNPs) from the HLA complex and six other loci were associated to OCB status. In SNPs selected for replication, combined analyses showed genome-wide significant association for two SNPs in the HLA complex; rs3129871 (p = 5.7×10⁻¹⁵) and rs3817963 (p = 5.7×10⁻¹⁰) correlating with the HLA-DRB1*15 and the HLA-DRB1*04 alleles, respectively. We also found suggestive association to one SNP in the Calsyntenin-2 gene (p = 8.83×10⁻⁷). In HLA-DRB1 analyses HLA-DRB1*15∶01 was a stronger risk factor for OCB positive than OCB negative MS, whereas HLA-DRB1*04∶04 was associated with increased risk of OCB negative MS and reduced risk of OCB positive MS. Protective effects of HLA-DRB1*01∶01 and HLA-DRB1*07∶01 were detected in both groups. The groups were different with regard to age at onset (AAO), MS outcome measures and gender. This study confirms both shared and distinct genetic risk for MS subtypes in the Scandinavian population defined by OCB status and indicates different clinical characteristics between the groups. This suggests differences in disease mechanisms between OCB negative and OCB positive MS with implications for patient management, which need to be further studied.     

Genome-Wide Association Study for Levels of Total Serum IgE Identifies HLA-C in a Japanese Population

Most of the previously reported loci for total immunoglobulin E (IgE) levels are related to Th2 cell-dependent pathways. We undertook a genome-wide association study (GWAS) to identify genetic loci responsible for IgE regulation. A total of 479,940 single nucleotide polymorphisms (SNPs) were tested for association with total serum IgE levels in 1180 Japanese adults. Fine-mapping with SNP imputation demonstrated 6 candidate regions: the PYHIN1/IFI16, MHC classes I and II, LEMD2, GRAMD1B, and chr13∶60576338 regions. Replication of these candidate loci in each region was assessed in 2 independent Japanese cohorts (n = 1110 and 1364, respectively). SNP rs3130941 in the HLA-C region was consistently associated with total IgE levels in 3 independent populations, and the meta-analysis yielded genome-wide significance (P = 1.07×10⁻¹⁰). Using our GWAS results, we also assessed the reproducibility of previously reported gene associations with total IgE levels. Nine of 32 candidate genes identified by a literature search were associated with total IgE levels after correction for multiple testing. Our findings demonstrate that SNPs in the HLA-C region are strongly associated with total serum IgE levels in the Japanese population and that some of the previously reported genetic associations are replicated across ethnic groups.     

Genome-Wide Association Analysis for Blood Lipid Traits Measured in Three Pig Populations Reveals a Substantial Level of Genetic Heterogeneity

Serum lipids are associated with myocardial infarction and cardiovascular disease in humans. Here we dissected the genetic architecture of blood lipid traits by applying genome-wide association studies (GWAS) in 1,256 pigs from Laiwu, Erhualian and Duroc × (Landrace × Yorkshire) populations, and a meta-analysis of GWAS in more than 2,400 pigs from five diverse populations. A total of 22 genomic loci surpassing the suggestive significance level were detected on 11 pig chromosomes (SSC) for six blood lipid traits. Meta-analysis of GWAS identified 5 novel loci associated with blood lipid traits. Comparison of GWAS loci across the tested populations revealed a substantial level of genetic heterogeneity for porcine blood lipid levels. We further evaluated the causality of nine polymorphisms nearby or within the APOB gene on SSC3 for serum LDL-C and TC levels. Of the 9 polymorphisms, an indel showed the most significant association with LDL-C and TC in Laiwu pigs. But the significant association was not identified in the White Duroc × Erhualian F₂ resource population, in which the QTL for LDL-C and TC was also detected on SSC3. This indicates that population-specific signals may exist for the SSC3 QTL. Further investigations are warranted to validate this assumption.     

Genome-wide association analyses reveal significant loci and strong candidate genes for growth and fatness traits in two pig populations

Background

Recently, genome-wide association studies (GWAS) have been reported on various pig traits. We performed a GWAS to analyze 22 traits related to growth and fatness on two pig populations: a White Duroc × Erhualian F₂ intercross population and a Chinese Sutai half-sib population.

Results

We identified 14 and 39 loci that displayed significant associations with growth and fatness traits at the genome-wide level and chromosome-wide level, respectively. The strongest association was between a 750 kb region on SSC7 (SSC for Sus scrofa) and backfat thickness at the first rib. This region had pleiotropic effects on both fatness and growth traits in F₂ animals and contained a promising candidate gene HMGA1 (high mobility group AT-hook 1). Unexpectedly, population genetic analysis revealed that the allele at this locus that reduces fatness and increases growth is derived from Chinese indigenous pigs and segregates in multiple Chinese breeds. The second strongest association was between the region around 82.85 Mb on SSC4 and average backfat thickness. PLAG1 (pleiomorphic adenoma gene 1), a gene under strong selection in European domestic pigs, is proximal to the top SNP and stands out as a strong candidate gene. On SSC2, a locus that significantly affects fatness traits mapped to the region around the IGF2 (insulin-like growth factor 2) gene but its non-imprinting inheritance excluded IGF2 as a candidate gene. A significant locus was also detected within a recombination cold spot that spans more than 30 Mb on SSCX, which hampered the identification of plausible candidate genes. Notably, no genome-wide significant locus was shared by the two experimental populations; different loci were observed that had both constant and time-specific effects on growth traits at different stages, which illustrates the complex genetic architecture of these traits.

Conclusions

We confirm several previously reported QTL and provide a list of novel loci for porcine growth and fatness traits in two experimental populations with Chinese Taihu and Western pigs as common founders. We showed that distinct loci exist for these traits in the two populations and identified HMGA1 and PLAG1 as strong candidate genes on SSC7 and SSC4, respectively.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-015-0089-5) contains supplementary material, which is available to authorized users.     

Quantitative trait loci for fatty acid composition in longissimus dorsi and abdominal fat: results from a White Duroc × Erhualian intercross F2 population

A whole-genome scan was performed on 660 F₂ animals including 250 barrows and 410 gilts in a White Duroc × Erhualian intercross population to detect quantitative trait loci (QTL) for fatty acid composition in the longissimus dorsi muscle and abdominal fat. A total of 153 QTL including 63 genome-wide significant QTL and 90 suggestive effects were identified for the traits measured. Significant effects were mainly evident on pig chromosomes (SSC) 4, 7, 8 and X. No association was detected on SSC3 and 11. In general, the QTL detected in this study showed distinct effects on fatty acid composition in the longissimus muscle and abdominal fat. The QTL for fatty acid composition in abdominal fat did not correspond to those identified previously in backfat and the majority of QTL for the muscle fatty acid composition were mapped to chromosomal regions different from previous studies. Two regions on SSC4 and SSC7 showed significant pleiotropic effects on monounsaturated (MUFA) and polyunsaturated fatty acid (PUFA) in both longissimus muscle and abdominal fat. Another two QTL with significant multi-faceted effects on MUFA and PUFA in the longissimus muscle were found each on SSC8 and SSCX. Chinese Erhualian alleles were associated with increased ratios of MUFA to saturated fatty acid at most of the QTL detected, showing beneficial effect in terms of human health.     

A Genome-Wide Association Study of Optic Disc Parameters

The optic nerve head is involved in many ophthalmic disorders, including common diseases such as myopia and open-angle glaucoma. Two of the most important parameters are the size of the optic disc area and the vertical cup-disc ratio (VCDR). Both are highly heritable but genetically largely undetermined. We performed a meta-analysis of genome-wide association (GWA) data to identify genetic variants associated with optic disc area and VCDR. The gene discovery included 7,360 unrelated individuals from the population-based Rotterdam Study I and Rotterdam Study II cohorts. These cohorts revealed two genome-wide significant loci for optic disc area, rs1192415 on chromosome 1p22 (p = 6.72×10⁻¹⁹) within 117 kb of the CDC7 gene and rs1900004 on chromosome 10q21.3-q22.1 (p = 2.67×10⁻³³) within 10 kb of the ATOH7 gene. They revealed two genome-wide significant loci for VCDR, rs1063192 on chromosome 9p21 (p = 6.15×10⁻¹¹) in the CDKN2B gene and rs10483727 on chromosome 14q22.3-q23 (p = 2.93×10⁻¹⁰) within 40 kbp of the SIX1 gene. Findings were replicated in two independent Dutch cohorts (Rotterdam Study III and Erasmus Rucphen Family study; N = 3,612), and the TwinsUK cohort (N = 843). Meta-analysis with the replication cohorts confirmed the four loci and revealed a third locus at 16q12.1 associated with optic disc area, and four other loci at 11q13, 13q13, 17q23 (borderline significant), and 22q12.1 for VCDR. ATOH7 was also associated with VCDR independent of optic disc area. Three of the loci were marginally associated with open-angle glaucoma. The protein pathways in which the loci of optic disc area are involved overlap with those identified for VCDR, suggesting a common genetic origin.     

Circulating Docosahexaenoic Acid Levels Are Associated with Fetal Insulin Sensitivity

Background

Arachidonic acid (AA; C20∶4 n-6) and docosahexaenoic acid (DHA; C22∶6 n-3) are important long-chain polyunsaturated fatty acids (LC-PUFA) in maintaining pancreatic beta-cell structure and function. Newborns of gestational diabetic mothers are more susceptible to the development of type 2 diabetes in adulthood. It is not known whether low circulating AA or DHA is involved in perinatally “programming” this susceptibility. This study aimed to assess whether circulating concentrations of AA, DHA and other fatty acids are associated with fetal insulin sensitivity or beta-cell function, and whether low circulating concentrations of AA or DHA are involved in compromised fetal insulin sensitivity in gestational diabetic pregnancies.

Methods and Principal Findings

In a prospective singleton pregnancy cohort, maternal (32-35 weeks gestation) and cord plasma fatty acids were assessed in relation to surrogate indicators of fetal insulin sensitivity (cord plasma glucose-to-insulin ratio, proinsulin concentration) and beta-cell function (proinsulin-to-insulin ratio) in 108 mother-newborn pairs. Cord plasma DHA levels (in percentage of total fatty acids) were lower comparing newborns of gestational diabetic (n = 24) vs. non-diabetic pregnancies (2.9% vs. 3.5%, P = 0.01). Adjusting for gestational age at blood sampling, lower cord plasma DHA levels were associated with lower fetal insulin sensitivity (lower glucose-to-insulin ratio, r = 0.20, P = 0.036; higher proinsulin concentration, r = −0.37, P <0.0001). The associations remained after adjustment for maternal and newborn characteristics. Cord plasma saturated fatty acids C18∶0 and C20∶0 were negatively correlated with fetal insulin sensitivity, but their levels were not different between gestational diabetic and non-diabetic pregnancies. Cord plasma AA levels were not correlated with fetal insulin sensitivity.

Conclusion

Low circulating DHA levels are associated with compromised fetal insulin sensitivity, and may be involved in perinatally “programming” the susceptibility to type 2 diabetes in the offspring of gestational diabetic mothers.     

Identification of Unusual Phospholipid Fatty Acyl Compositions of Acanthamoeba castellanii

Acanthamoeba are opportunistic protozoan pathogens that may lead to sight-threatening keratitis and fatal granulomatous encephalitis. The successful prognosis requires early diagnosis and differentiation of pathogenic Acanthamoeba followed by aggressive treatment regimen. The plasma membrane of Acanthamoeba consists of 25% phospholipids (PL). The presence of C20 and, recently reported, 28- and 30-carbon fatty acyl residues is characteristic of amoeba PL. A detailed knowledge about this unusual PL composition could help to differentiate Acanthamoeba from other parasites, e.g. bacteria and develop more efficient treatment strategies. Therefore, the detailed PL composition of Acanthamoeba castellanii was investigated by ³¹P nuclear magnetic resonance spectroscopy, thin-layer chromatography, gas chromatography, high performance liquid chromatography and liquid chromatography-mass spectrometry. Normal and reversed phase liquid chromatography coupled with mass spectrometric detection was used for detailed characterization of the fatty acyl composition of each detected PL. The most abundant fatty acyl residues in each PL class were octadecanoyl (18∶0), octadecenoyl (18∶1 Δ9) and hexadecanoyl (16∶0). However, some selected PLs contained also very long fatty acyl chains: the presence of 28- and 30-carbon fatty acyl residues was confirmed in phosphatidylethanolamine (PE), phosphatidylserine, phosphatidic acid and cardiolipin. The majority of these fatty acyl residues were also identified in PE that resulted in the following composition: 28∶1/20∶2, 30∶2/18∶1, 28∶0/20∶2, 30∶2/20∶4 and 30∶3/20∶3. The PL of amoebae are significantly different in comparison to other cells: we describe here for the first time unusual, very long chain fatty acids with Δ⁵-unsaturation (30∶3^5,21,24) and 30∶2^21,24 localized exclusively in specific phospholipid classes of A. castellanii protozoa that could serve as specific biomarkers for the presence of these microorganisms.     

Long-Term Low Carbohydrate Diet Leads to Deleterious Metabolic Manifestations in Diabetic Mice

We investigated long-term effects of low carbohydrate diets on wild type mice, streptozotocin-injected and KKAy obese diabetic mice. These mice were pair-fed three different types of diets, standard chow (SC, C∶P∶F = 63∶15∶22), a low carbohydrate (LC, C∶P∶F = 38∶25∶37) diet and a severely carbohydrate restricted (SR, C∶P∶F = 18∶45∶37) diet for 16 weeks. Despite comparable body weights and serum lipid profiles, wild type and diabetic mice fed the low carbohydrate diets exhibited lower insulin sensitivity and this reduction was dependent on the amount of carbohydrate in the diet. When serum fatty acid compositions were investigated, monounsaturation capacity, i.e. C16:1/C16:0 and C18:1/C18:0, was impaired in all murine models fed the low carbohydrate diets, consistent with the decreased expression of hepatic stearoyl-CoA desaturase-1 (SCD1). Interestingly, both the hepatic expressions and serum levels of fibroblast growth factor 21 (FGF21), which might be related to longevity, were markedly decreased in both wild type and KKAy mice fed the SR diet. Taking into consideration that fat compositions did not differ between the LC and SR diets, we conclude that low carbohydrate diets have deleterious metabolic effects in both wild type and diabetic mice, which may explain the association between diets relatively low in carbohydrate and the elevated risk of cardiovascular events observed in clinical studies.     

Genome-Wide Association Studies Identify the Loci for 5 Exterior Traits in a Large White × Minzhu Pig Population

As one of the main breeding selection criteria, external appearance has special economic importance in the hog industry. In this study, an Illumina Porcine SNP60 BeadChip was used to conduct a genome-wide association study (GWAS) in 605 pigs of the F₂ generation derived from a Large White × Minzhu intercross. Traits under study were abdominal circumference (AC), body height (BH), body length (BL), cannon bone circumference (CBC), chest depth (CD), chest width (CW), rump circumference (RC), rump width (RW), scapula width (SW), and waist width (WW). A total of 138 SNPs (the most significant being MARC0033464) on chromosome 7 were found to be associated with BH, BL, CBC, and RC (P-value  = 4.15E-6). One SNP on chromosome 1 was found to be associated with CD at genome-wide significance levels. The percentage phenotypic variance of these significant SNPs ranged from 0.1–25.48%. Moreover, a conditional analysis revealed that the significant SNPs were derived from a single quantitative trait locus (QTL) and indicated additional chromosome-wide significant association for 25 SNPs on SSC4 (BL, CBC) and 9 SNPs on SSC7 (RC). Linkage analysis revealed two complete linkage disequilibrium haplotype blocks that contained seven and four SNPs, respectively. In block 1, the most significant SNP, MARC0033464, was present. Annotations from pig reference genome suggested six genes (GRM4, HMGA1, NUDT3, RPS10, SPDEF and PACSIN1) in block 1 (495 kb), and one gene (SCUBE3) in block 3 (124 kb). Functional analysis indicated that HMGA1 and SCUBE3 genes are the potential genes controlling BH, BL, and RC in pigs, with an application in breeding programs. We screened several candidate intervals and genes based on SNP location and gene function, and predicted their function using bioinformatics analyses.     

Association between Plasma Nonesterified Fatty Acids Species and Adipose Tissue Fatty Acid Composition

Fatty acid composition of adipose tissue (AT) is an established long-term biomarker for fatty acid (FA) intake and status, but AT samples are not easily available. Nonesterified FA composition in plasma (pNEFA) may be a good indicator of AT FA composition, because pNEFA are mainly generated by AT lipolysis. We investigated the correlation of 42 pNEFA and subcutaneous as well as visceral AT FA in 27 non-diabetic women with a median BMI of 36 kg/m² (Q_0.25: 25 kg/m²; Q_0.75: 49 kg/m²). Close correlations of pNEFA and AT FA were found for odd-chain FA (15∶0 r = 0.838 and 0.862 for subcutaneous and visceral AT, respectively) and omega-3 FA (22∶6 r = 0.719/0.535), while no significant or low correlations were found for other FA including 18∶1 (r = 0.384/0.325) and 20∶4 (r = 0.386/0.266). Close correlations of pNEFA and AT FA were found for essential fatty acids, like 18∶2 (r = 0.541/0.610) and 20∶5 (r = 0.561/0.543). The lower correlation for some pNEFA species with AT FA indicates that the variation of most pNEFA is significantly affected by other FA sources and flux of FA to tissue, in addition to release from AT. A relevant influence of BMI on the level of correlation was shown for saturated FA. NEFA analysis in fasted plasma can serve as a virtual AT biopsy for some FA, and as a biomarker for intake of dairy products and sea fish.     

Importance of Fatty Acid Compositions in Patients with Peripheral Arterial Disease

Objective

Importance of fatty acid components and imbalances has emerged in coronary heart disease. In this study, we analyzed fatty acids and ankle-brachial index (ABI) in a Japanese cohort.

Methods

Peripheral arterial disease (PAD) was diagnosed in 101 patients by ABI ≤0.90 and/or by angiography. Traditional cardiovascular risk factors and components of serum fatty acids were examined in all patients (mean age 73.2±0.9 years; 81 males), and compared with those in 373 age- and sex-matched control subjects with no evidence of PAD.

Results

The presence of PAD (mean ABI: 0.71±0.02) was independently associated with low levels of gamma-linolenic acid (GLA) (OR: 0.90; 95% CI: 0.85–0.96; P = 0.002), eicosapentaenoic acid∶arachidonic acid (EPA∶AA) ratio (OR: 0.38; 95% CI: 0.17–0.86; P = 0.021), and estimated glomerular filtration rate (OR: 0.97; 95% CI: 0.96–0.98; P<0.0001), and with a high hemoglobin A1c level (OR: 1.34; 95% CI: 1.06–1.69; P = 0.013). Individuals with lower levels of GLA (≤7.95 µg/mL) and a lower EPA∶AA ratio (≤0.55) had the lowest ABI (0.96±0.02, N = 90), while the highest ABI (1.12±0.01, N = 78) was observed in individuals with higher values of both GLA and EPA∶AA ratio (P<0.0001).

Conclusion

A low level of GLA and a low EPA∶AA ratio are independently associated with the presence of PAD. Specific fatty acid abnormalities and imbalances could lead to new strategies for risk stratification and prevention in PAD patients.     

DNA Methylation in an Enhancer Region of the FADS Cluster Is Associated with FADS Activity in Human Liver

Levels of omega-6 (n-6) and omega-3 (n-3), long chain polyunsaturated fatty acids (LcPUFAs) such as arachidonic acid (AA; 20∶4, n-6), eicosapentaenoic acid (EPA; 20∶5, n-3) and docosahexaenoic acid (DHA; 22∶6, n-3) impact a wide range of biological activities, including immune signaling, inflammation, and brain development and function. Two desaturase steps (Δ6, encoded by FADS2 and Δ5, encoded by FADS1) are rate limiting in the conversion of dietary essential 18 carbon PUFAs (18C-PUFAs) such as LA (18∶2, n-6) to AA and α-linolenic acid (ALA, 18∶3, n-3) to EPA and DHA. GWAS and candidate gene studies have consistently identified genetic variants within FADS1 and FADS2 as determinants of desaturase efficiencies and levels of LcPUFAs in circulating, cellular and breast milk lipids. Importantly, these same variants are documented determinants of important cardiovascular disease risk factors (total, LDL, and HDL cholesterol, triglycerides, CRP and proinflammatory eicosanoids). FADS1 and FADS2 lie head-to-head (5′ to 5′) in a cluster configuration on chromosome 11 (11q12.2). There is considerable linkage disequilibrium (LD) in this region, where multiple SNPs display association with LcPUFA levels. For instance, rs174537, located ∼15 kb downstream of FADS1, is associated with both FADS1 desaturase activity and with circulating AA levels (p-value for AA levels = 5.95×10⁻⁴⁶) in humans. To determine if DNA methylation variation impacts FADS activities, we performed genome-wide allele-specific methylation (ASM) with rs174537 in 144 human liver samples. This approach identified highly significant ASM with CpG sites between FADS1 and FADS2 in a putative enhancer signature region, leading to the hypothesis that the phenotypic associations of rs174537 are likely due to methylation differences. In support of this hypothesis, methylation levels of the most significant probe were strongly associated with FADS1 and, to a lesser degree, FADS2 activities.     

Genome-Wide Association of Body Fat Distribution in African Ancestry Populations Suggests New Loci

Central obesity, measured by waist circumference (WC) or waist-hip ratio (WHR), is a marker of body fat distribution. Although obesity disproportionately affects minority populations, few studies have conducted genome-wide association study (GWAS) of fat distribution among those of predominantly African ancestry (AA). We performed GWAS of WC and WHR, adjusted and unadjusted for BMI, in up to 33,591 and 27,350 AA individuals, respectively. We identified loci associated with fat distribution in AA individuals using meta-analyses of GWA results for WC and WHR (stage 1). Overall, 25 SNPs with single genomic control (GC)-corrected p-values<5.0×10⁻⁶ were followed-up (stage 2) in AA with WC and with WHR. Additionally, we interrogated genomic regions of previously identified European ancestry (EA) WHR loci among AA. In joint analysis of association results including both Stage 1 and 2 cohorts, 2 SNPs demonstrated association, rs2075064 at LHX2, p = 2.24×10⁻⁸ for WC-adjusted-for-BMI, and rs6931262 at RREB1, p = 2.48×10⁻⁸ for WHR-adjusted-for-BMI. However, neither signal was genome-wide significant after double GC-correction (LHX2: p = 6.5×10⁻⁸; RREB1: p = 5.7×10⁻⁸). Six of fourteen previously reported loci for waist in EA populations were significant (p<0.05 divided by the number of independent SNPs within the region) in AA studied here (TBX15-WARS2, GRB14, ADAMTS9, LY86, RSPO3, ITPR2-SSPN). Further, we observed associations with metabolic traits: rs13389219 at GRB14 associated with HDL-cholesterol, triglycerides, and fasting insulin, and rs13060013 at ADAMTS9 with HDL-cholesterol and fasting insulin. Finally, we observed nominal evidence for sexual dimorphism, with stronger results in AA women at the GRB14 locus (p for interaction = 0.02). In conclusion, we identified two suggestive loci associated with fat distribution in AA populations in addition to confirming 6 loci previously identified in populations of EA. These findings reinforce the concept that there are fat distribution loci that are independent of generalized adiposity.     

Large-scale Metabolomic Profiling Identifies Novel Biomarkers for Incident Coronary Heart Disease

Analyses of circulating metabolites in large prospective epidemiological studies could lead to improved prediction and better biological understanding of coronary heart disease (CHD). We performed a mass spectrometry-based non-targeted metabolomics study for association with incident CHD events in 1,028 individuals (131 events; 10 y. median follow-up) with validation in 1,670 individuals (282 events; 3.9 y. median follow-up). Four metabolites were replicated and independent of main cardiovascular risk factors [lysophosphatidylcholine 18∶1 (hazard ratio [HR] per standard deviation [SD] increment = 0.77, P-value<0.001), lysophosphatidylcholine 18∶2 (HR = 0.81, P-value<0.001), monoglyceride 18∶2 (MG 18∶2; HR = 1.18, P-value = 0.011) and sphingomyelin 28∶1 (HR = 0.85, P-value = 0.015)]. Together they contributed to moderate improvements in discrimination and re-classification in addition to traditional risk factors (C-statistic: 0.76 vs. 0.75; NRI: 9.2%). MG 18∶2 was associated with CHD independently of triglycerides. Lysophosphatidylcholines were negatively associated with body mass index, C-reactive protein and with less evidence of subclinical cardiovascular disease in additional 970 participants; a reverse pattern was observed for MG 18∶2. MG 18∶2 showed an enrichment (P-value = 0.002) of significant associations with CHD-associated SNPs (P-value = 1.2×10⁻⁷ for association with rs964184 in the ZNF259/APOA5 region) and a weak, but positive causal effect (odds ratio = 1.05 per SD increment in MG 18∶2, P-value = 0.05) on CHD, as suggested by Mendelian randomization analysis. In conclusion, we identified four lipid-related metabolites with evidence for clinical utility, as well as a causal role in CHD development.     

Association of Genetic Variants Influencing Lipid Levels with Coronary Artery Disease in Japanese Individuals

Background/Objective

In Japanese populations, we performed a replication study of genetic loci previously identified in European-descent populations as being associated with lipid levels and risk of coronary artery disease (CAD).

Methods

We genotyped 48 single nucleotide polymorphisms (SNPs) from 22 candidate loci that had previously been identified by genome-wide association (GWA) meta-analyses for low-density lipoprotein cholesterol (LDL-C), high-density lipoprotein cholesterol (HDL-C), and/or triglycerides in Europeans. We selected 22 loci with 2 parallel tracks from 95 reported loci: 16 significant loci (p<1×10⁻³⁰ in Europeans) and 6 other loci including those with suggestive evidence of lipid associations in 1292 GWA-scanned Japanese samples. Genotyping was done in 4990 general population samples, and 1347 CAD cases and 1337 controls. For 9 SNPs, we further examined CAD associations in an additional panel of 3052 CAD cases and 6335 controls.

Principal Findings

Significant lipid associations (one-tailed p<0.05) were replicated for 18 of 22 loci in Japanese samples, with significant inter-ethnic heterogeneity at 4 loci–APOB, APOE-C1, CETP, and APOA5–and allelic heterogeneity. The strongest association was detected at APOE rs7412 for LDL-C (p = 1.3×10⁻⁴¹), CETP rs3764261 for HDL-C (p = 5.2×10⁻²⁴), and APOA5 rs662799 for triglycerides (p = 5.8×10⁻⁵⁴). CAD association was replicated and/or verified for 4 loci: SORT1 rs611917 (p = 1.7×10⁻⁸), APOA5 rs662799 (p = 0.0014), LDLR rs1433099 (p = 2.1×10⁻⁷), and APOE rs7412 (p = 6.1×10⁻¹³).

Conclusions

Our results confirm that most of the tested lipid loci are associated with lipid traits in the Japanese, further indicating that in genetic susceptibility to lipid levels and CAD, the related metabolic pathways are largely common across the populations, while causal variants at individual loci can be population-specific.     

Gene-Centric Meta-Analysis of Lipid Traits in African,East Asian and Hispanic Populations

Meta-analyses of European populations has successfully identified genetic variants in over 100 loci associated with lipid levels, but our knowledge in other ethnicities remains limited. To address this, we performed dense genotyping of ∼2,000 candidate genes in 7,657 African Americans, 1,315 Hispanics and 841 East Asians, using the IBC array, a custom ∼50,000 SNP genotyping array. Meta-analyses confirmed 16 lipid loci previously established in European populations at genome-wide significance level, and found multiple independent association signals within these lipid loci. Initial discovery and in silico follow-up in 7,000 additional African American samples, confirmed two novel loci: rs5030359 within ICAM1 is associated with total cholesterol (TC) and low-density lipoprotein cholesterol (LDL-C) (p = 8.8×10⁻⁷ and p = 1.5×10⁻⁶ respectively) and a nonsense mutation rs3211938 within CD36 is associated with high-density lipoprotein cholesterol (HDL-C) levels (p = 13.5×10⁻¹²). The rs3211938-G allele, which is nearly absent in European and Asian populations, has been previously found to be associated with CD36 deficiency and shows a signature of selection in Africans and African Americans. Finally, we have evaluated the effect of SNPs established in European populations on lipid levels in multi-ethnic populations and show that most known lipid association signals span across ethnicities. However, differences between populations, especially differences in allele frequency, can be leveraged to identify novel signals, as shown by the discovery of ICAM1 and CD36 in the current report.