Discovery and refinement of loci associated with lipid levels

A genome-wide association meta-analysis on apolipoprotein A-IV concentrations

Apolipoprotein A-IV (apoA-IV) is a major component of HDL and chylomicron particles and is involved in reverse cholesterol transport. It is an early marker of impaired renal function. We aimed to identify genetic loci associated with apoA-IV concentrations and to investigate relationships with known susceptibility loci for kidney function and lipids. A genome-wide association meta-analysis on apoA-IV concentrations was conducted in five population-based cohorts (n = 13,813) followed by two additional replication studies (n = 2,267) including approximately 10 M SNPs. Three independent SNPs from two genomic regions were significantly associated with apoA-IV concentrations: rs1729407 near APOA4 (P = 6.77 × 10 ^{-  44}), rs5104 in APOA4 (P = 1.79 × 10^-24) and rs4241819 in KLKB1 (P = 5.6 × 10^-14). Additionally, a look-up of the replicated SNPs in downloadable GWAS meta-analysis results was performed on kidney function (defined by eGFR), HDL-cholesterol and triglycerides. From these three SNPs mentioned above, only rs1729407 showed an association with HDL-cholesterol (P = 7.1 × 10 ^{-  07}). Moreover, weighted SNP-scores were built involving known susceptibility loci for the aforementioned traits (53, 70 and 38 SNPs, respectively) and were associated with apoA-IV concentrations. This analysis revealed a significant and an inverse association for kidney function with apoA-IV concentrations (P = 5.5 × 10^-05). Furthermore, an increase of triglyceride-increasing alleles was found to decrease apoA-IV concentrations (P = 0.0078). In summary, we identified two independent SNPs located in or next the APOA4 gene and one SNP in KLKB1 The association of KLKB1 with apoA-IV suggests an involvement of apoA-IV in renal metabolism and/or an interaction within HDL particles. Analyses of SNP-scores indicate potential causal effects of kidney function and by lesser extent triglycerides on apoA-IV concentrations.

Large-Scale Exome-wide Association Analysis Identifies Loci for White Blood Cell Traits and Pleiotropy with Immune-Mediated Diseases

White blood cells play diverse roles in innate and adaptive immunity. Genetic association analyses of phenotypic variation in circulating white blood cell (WBC) counts from large samples of otherwise healthy individuals can provide insights into genes and biologic pathways involved in production, differentiation, or clearance of particular WBC lineages (myeloid, lymphoid) and also potentially inform the genetic basis of autoimmune, allergic, and blood diseases. We performed an exome array-based meta-analysis of total WBC and subtype counts (neutrophils, monocytes, lymphocytes, basophils, and eosinophils) in a multi-ancestry discovery and replication sample of ∼157,622 individuals from 25 studies. We identified 16 common variants (8 of which were coding variants) associated with one or more WBC traits, the majority of which are pleiotropically associated with autoimmune diseases. Based on functional annotation, these loci included genes encoding surface markers of myeloid, lymphoid, or hematopoietic stem cell differentiation (CD69, CD33, CD87), transcription factors regulating lineage specification during hematopoiesis (ASXL1, IRF8, IKZF1, JMJD1C, ETS2-PSMG1), and molecules involved in neutrophil clearance/apoptosis (C10orf54, LTA), adhesion (TNXB), or centrosome and microtubule structure/function (KIF9, TUBD1). Together with recent reports of somatic ASXL1 mutations among individuals with idiopathic cytopenias or clonal hematopoiesis of undetermined significance, the identification of a common regulatory 3' UTR variant of ASXL1 suggests that both germline and somatic ASXL1 mutations contribute to lower blood counts in otherwise asymptomatic individuals. These association results shed light on genetic mechanisms that regulate circulating WBC counts and suggest a prominent shared genetic architecture with inflammatory and autoimmune diseases.

Exome Genotyping Identifies Pleiotropic Variants Associated with Red Blood Cell Traits

Red blood cell (RBC) traits are important heritable clinical biomarkers and modifiers of disease severity. To identify coding genetic variants associated with these traits, we conducted meta-analyses of seven RBC phenotypes in 130,273 multi-ethnic individuals from studies genotyped on an exome array. After conditional analyses and replication in 27,480 independent individuals, we identified 16 new RBC variants. We found low-frequency missense variants in MAP1A (rs55707100, minor allele frequency [MAF] = 3.3%, p = 2 × 10(-10) for hemoglobin [HGB]) and HNF4A (rs1800961, MAF = 2.4%, p < 3 × 10(-8) for hematocrit [HCT] and HGB). In African Americans, we identified a nonsense variant in CD36 associated with higher RBC distribution width (rs3211938, MAF = 8.7%, p = 7 × 10(-11)) and showed that it is associated with lower CD36 expression and strong allelic imbalance in ex vivo differentiated human erythroblasts. We also identified a rare missense variant in ALAS2 (rs201062903, MAF = 0.2%) associated with lower mean corpuscular volume and mean corpuscular hemoglobin (p < 8 × 10(-9)). Mendelian mutations in ALAS2 are a cause of sideroblastic anemia and erythropoietic protoporphyria. Gene-based testing highlighted three rare missense variants in PKLR, a gene mutated in Mendelian non-spherocytic hemolytic anemia, associated with HGB and HCT (SKAT p < 8 × 10(-7)). These rare, low-frequency, and common RBC variants showed pleiotropy, being also associated with platelet, white blood cell, and lipid traits. Our association results and functional annotation suggest the involvement of new genes in human erythropoiesis. We also confirm that rare and low-frequency variants play a role in the architecture of complex human traits, although their phenotypic effect is generally smaller than originally anticipated.

Platelet-Related Variants Identified by Exomechip Meta-analysis in 157,293 Individuals

Platelet production, maintenance, and clearance are tightly controlled processes indicative of platelets' important roles in hemostasis and thrombosis. Platelets are common targets for primary and secondary prevention of several conditions. They are monitored clinically by complete blood counts, specifically with measurements of platelet count (PLT) and mean platelet volume (MPV). Identifying genetic effects on PLT and MPV can provide mechanistic insights into platelet biology and their role in disease. Therefore, we formed the Blood Cell Consortium (BCX) to perform a large-scale meta-analysis of Exomechip association results for PLT and MPV in 157,293 and 57,617 individuals, respectively. Using the low-frequency/rare coding variant-enriched Exomechip genotyping array, we sought to identify genetic variants associated with PLT and MPV. In addition to confirming 47 known PLT and 20 known MPV associations, we identified 32 PLT and 18 MPV associations not previously observed in the literature across the allele frequency spectrum, including rare large effect (FCER1A), low-frequency (IQGAP2, MAP1A, LY75), and common (ZMIZ2, SMG6, PEAR1, ARFGAP3/PACSIN2) variants. Several variants associated with PLT/MPV (PEAR1, MRVI1, PTGES3) were also associated with platelet reactivity. In concurrent BCX analyses, there was overlap of platelet-associated variants with red (MAP1A, TMPRSS6, ZMIZ2) and white (PEAR1, ZMIZ2, LY75) blood cell traits, suggesting common regulatory pathways with shared genetic architecture among these hematopoietic lineages. Our large-scale Exomechip analyses identified previously undocumented associations with platelet traits and further indicate that several complex quantitative hematological, lipid, and cardiovascular traits share genetic factors.

Contrasting the Genetic Architecture of 30 Complex Traits from Summary Association Data

Variance-component methods that estimate the aggregate contribution of large sets of variants to the heritability of complex traits have yielded important insights into the genetic architecture of common diseases. Here, we introduce methods that estimate the total trait variance explained by the typed variants at a single locus in the genome (local SNP heritability) from genome-wide association study (GWAS) summary data while accounting for linkage disequilibrium among variants. We applied our estimator to ultra-large-scale GWAS summary data of 30 common traits and diseases to gain insights into their local genetic architecture. First, we found that common SNPs have a high contribution to the heritability of all studied traits. Second, we identified traits for which the majority of the SNP heritability can be confined to a small percentage of the genome. Third, we identified GWAS risk loci where the entire locus explains significantly more variance in the trait than the GWAS reported variants. Finally, we identified loci that explain a significant amount of heritability across multiple traits.

TTC39B deficiency stabilizes LXR reducing both atherosclerosis and steatohepatitis

Cellular mechanisms that mediate steato-hepatitis, an increasingly prevalent condition in the Western world for which no therapies are available¹, are poorly understood. Despite the fact its synthetic agonists induce fatty liver, the Liver X receptor (LXR) transcription factor remains a target of interest because of its anti-atherogenic, cholesterol removal and anti-inflammatory activities. We discovered that tetratricopeptide repeat (TPR) domain protein 39B (Ttc39b, C9orf52) (T39), a high density lipoprotein (HDL) gene discovered in human genome wide association studies (GWAS)², promotes the ubiquitination and degradation of LXR. Chow-fed T39^-/- mice displayed increased HDL cholesterol levels associated with increased enterocyte ATP binding cassette transporter A1 (Abca1) expression and increased LXR protein without change in LXR mRNA. When challenged with a high fat/high cholesterol/bile salt (HF/HC/BS) diet, T39^-/- mice or mice with hepatocyte-specific T39 deficiency showed increased hepatic LXR protein and target gene expression, and unexpectedly protection from steato-hepatitis and death. Western Type Diet (WTD)-fed Low density lipoprotein receptor (Ldlr)^-/-T39^-/- mice showed decreased fatty liver, increased HDL, decreased LDL and reduced atherosclerosis. In addition to increasing hepatic Abcg5/8 expression and limiting dietary cholesterol absorption, T39 deficiency inhibited hepatic sterol regulatory element binding protein 1 (SREBP-1, ADD1) processing. This was explained by an increase in microsomal phospholipids containing polyunsaturated fatty acids (PUFA), linked to an LXRα-dependent increase in expression of enzymes mediating PC biosynthesis and incorporation of PUFA into phospholipids. The preservation of endogenous LXR protein activates a beneficial profile of gene expression that promotes cholesterol removal and inhibits lipogenesis. T39 inhibition could be an effective strategy for reducing both steato-hepatitis and atherosclerosis.

An Optimal Bahadur-Efficient Method in Detection of Sparse Signals with Applications to Pathway Analysis in Sequencing Association Studies

Next-generation sequencing data pose a severe curse of dimensionality, complicating traditional "single marker-single trait" analysis. We propose a two-stage combined p-value method for pathway analysis. The first stage is at the gene level, where we integrate effects within a gene using the Sequence Kernel Association Test (SKAT). The second stage is at the pathway level, where we perform a correlated Lancaster procedure to detect joint effects from multiple genes within a pathway. We show that the Lancaster procedure is optimal in Bahadur efficiency among all combined p-value methods. The Bahadur efficiency,[Formula: see text], compares sample sizes among different statistical tests when signals become sparse in sequencing data, i.e. ε →0. The optimal Bahadur efficiency ensures that the Lancaster procedure asymptotically requires a minimal sample size to detect sparse signals ([Formula: see text]). The Lancaster procedure can also be applied to meta-analysis. Extensive empirical assessments of exome sequencing data show that the proposed method outperforms Gene Set Enrichment Analysis (GSEA). We applied the competitive Lancaster procedure to meta-analysis data generated by the Global Lipids Genetics Consortium to identify pathways significantly associated with high-density lipoprotein cholesterol, low-density lipoprotein cholesterol, triglycerides, and total cholesterol.

Whole-genome sequencing in French Canadians from Quebec

Genome-wide association studies (GWAS) have had a tremendous success in the identification of common DNA sequence variants associated with complex human diseases and traits. However, because of their design, GWAS are largely inappropriate to characterize the role of rare and low-frequency DNA variants on human phenotypic variation. Rarer genetic variation is geographically more restricted, supporting the need for local whole-genome sequencing (WGS) efforts to study these variants in specific populations. Here, we present the first large-scale low-pass WGS of the French-Canadian population. Specifically, we sequenced at ~5.6× coverage the whole genome of 1970 French Canadians recruited by the Montreal Heart Institute Biobank and identified 29 million bi-allelic variants (31 % novel), including 19 million variants with a minor allele frequency (MAF) <0.5 %. Genotypes from the WGS data are highly concordant with genotypes obtained by exome array on the same individuals (99.8 %), even when restricting this analysis to rare variants (MAF <0.5, 99.9 %) or heterozygous sites (98.9 %). To further validate our data set, we showed that we can effectively use it to replicate several genetic associations with myocardial infarction risk and blood lipid levels. Furthermore, we analyze the utility of our WGS data set to generate a French-Canadian-specific imputation reference panel and to infer population structure in the Province of Quebec. Our results illustrate the value of low-pass WGS to study the genetics of human diseases in the founder French-Canadian population.

Meta-analysis of 49 549 individuals imputed with the 1000 Genomes Project reveals an exonic damaging variant in ANGPTL4 determining fasting TG levels

BACKGROUND

So far, more than 170 loci have been associated with circulating lipid levels through genome-wide association studies (GWAS). These associations are largely driven by common variants, their function is often not known, and many are likely to be markers for the causal variants. In this study we aimed to identify more new rare and low-frequency functional variants associated with circulating lipid levels.

METHODS

We used the 1000 Genomes Project as a reference panel for the imputations of GWAS data from ∼60 000 individuals in the discovery stage and ∼90 000 samples in the replication stage.

RESULTS

Our study resulted in the identification of five new associations with circulating lipid levels at four loci. All four loci are within genes that can be linked biologically to lipid metabolism. One of the variants, rs116843064, is a damaging missense variant within the ANGPTL4 gene.

CONCLUSIONS

This study illustrates that GWAS with high-scale imputation may still help us unravel the biological mechanism behind circulating lipid levels.

A genome-wide study of lipid response to fenofibrate in Caucasians: a combined analysis of the GOLDN and ACCORD studies

BACKGROUND

Fibrates are commonly prescribed for hypertriglyceridemia, but they also lower LDL cholesterol and increase HDL cholesterol. Large interindividual variations in lipid response suggest that some patients may benefit more than others and genetic studies could help identify such patients.

METHODS

We carried out the first genome-wide association study of lipid response to fenofibrate using data from two well-characterized clinical trials: the Genetics of Lipid Lowering Drugs and Diet Network (GOLDN) Study and the Action to Control Cardiovascular Risk in Diabetes (ACCORD) Study. Genome-wide association study data from both studies were imputed to the 1000 Genomes CEU reference panel (phase 1). Lipid response was modeled as the log ratio of the post-treatment lipid level to the pretreatment level. Linear mixed models (GOLDN, N=813 from 173 families) and linear regression models (ACCORD, N=781) adjusted for pretreatment lipid level, demographic variables, clinical covariates, and ancestry were used to evaluate the association of genetic markers with lipid response. Among Caucasians, the results were combined using inverse-variance weighted fixed-effects meta-analyses. The main findings from the meta-analyses were examined in other ethnic groups from the HyperTG study (N=267 Hispanics) and ACCORD (N=83 Hispanics, 138 African Americans).

RESULTS

A known lipid locus harboring the pre-B-cell leukemia homeobox 4 (PBX4) gene on chromosome 19 is important for LDL cholesterol response to fenofibrate (smallest P=1.5×10). The main results replicated with nominal statistical significance in Hispanics from ACCORD (P<0.05).

CONCLUSION

Future research should evaluate the usefulness of this locus to refine clinical strategies for lipid-lowering treatments.

Triglycerides and Triglyceride-Rich Lipoproteins in the Causal Pathway of Cardiovascular Disease

Epidemiologic and clinical studies suggest that elevated triglyceride levels are a biomarker of cardiovascular (CV) risk. Consistent with these findings, recent genetic evidence from mutational analyses, genome-wide association studies, and Mendelian randomization studies provide robust evidence that triglycerides and triglyceride-rich lipoproteins are in the causal pathway for atherosclerotic CV disease, indicating that they may play a pathogenic role, much like low-density lipoprotein cholesterol (LDL-C). Although statins are the cornerstone of dyslipidemia management, high triglyceride levels may persist in some patients despite statin therapy. Several triglyceride-lowering agents are available, including fibrates, niacin, and omega-3 fatty acids, of which prescription omega-3 fatty acids have the best tolerability and safety profile. In clinical studies, omega-3 fatty acids have been shown to reduce triglyceride levels, but products containing both eicosapentaenoic acid and docosahexaenoic acid may increase LDL-C levels. Icosapent ethyl, a high-purity eicosapentaenoic acid-only product, does not raise LDL-C levels and also reduces triglyceride, non-high-density lipoprotein cholesterol, and triglyceride-rich lipoprotein levels. In conclusion, omega-3 fatty acids are currently being evaluated in large CV outcome studies in statin-treated patients; these studies should help to elucidate the causative role of triglycerides in atherosclerotic CV disease.

Molecular and pathophysiological aspects of metal ion uptake by the zinc transporter ZIP8 (SLC39A8)

Zinc ion (Zn2+) is essential for life; its deficiency in the human body could cause stunted growth, anemia and susceptibility to infection. The Zn transporter ZIP8 (also known as SLC39A8) is an important Zn2+ importer; aberrant Zn2+ influx mediated by ZIP8 can lead to the pathogenesis of osteoarthritis and inflammatory diseases. ZIP8 also mediates the cellular uptake of divalent metal ions including iron, manganese, and the toxic heavy metal cadmium. Individuals with SLC39A8 mutations and transgenic mouse models are starting to reveal the critical role that this gene plays in embryonic development and the metabolism of essential metal ions. Here we summarize our current understanding of ZIP8's function and regulation, at both the molecular and biological levels. We also review the association of ZIP8 with various diseases and its linkage with complex disorders like obesity, hypertension, and schizophrenia as revealed by several large genome-wide association studies.

Molecular Characterization of the Lipid Genome-Wide Association Study Signal on Chromosome 18q11.2 Implicates HNF4A-Mediated Regulation of the TMEM241 Gene

OBJECTIVE

We recently identified a locus on chromosome 18q11.2 for high serum triglycerides in Mexicans. We hypothesize that the lead genome-wide association study single-nucleotide polymorphism rs9949617, or its linkage disequilibrium proxies, regulates 1 of the 5 genes in the triglyceride-associated region.

APPROACH AND RESULTS

We performed a linkage disequilibrium analysis and found 9 additional variants in linkage disequilibrium (r(2)>0.7) with the lead single-nucleotide polymorphism. To select the variants for functional analyses, we annotated the 10 variants using DNase I hypersensitive sites, transcription factor and chromatin states and identified rs17259126 as the lead candidate variant for functional in vitro validation. Using luciferase transcriptional reporter assay in liver HepG2 cells, we found that the G allele exhibits a significantly lower effect on transcription (P<0.05). The electrophoretic mobility shift and ChIPqPCR (chromatin immunoprecipitation coupled with quantitative polymerase chain reaction) assays confirmed that the minor G allele of rs17259126 disrupts an hepatocyte nuclear factor 4 α-binding site. To find the regional candidate gene, we performed a local expression quantitative trait locus analysis and found that rs17259126 and its linkage disequilibrium proxies alter expression of the regional transmembrane protein 241 (TMEM241) gene in 795 adipose RNAs from the Metabolic Syndrome In Men (METSIM) cohort (P=6.11×10(-07)-5.80×10(-04)). These results were replicated in expression profiles of TMEM241 from the Multiple Tissue Human Expression Resource (MuTHER; n=856).

CONCLUSIONS

The Mexican genome-wide association study signal for high serum triglycerides on chromosome 18q11.2 harbors a regulatory single-nucleotide polymorphism, rs17259126, which disrupts normal hepatocyte nuclear factor 4 α binding and decreases the expression of the regional TMEM241 gene. Our data suggest that decreased transcript levels of TMEM241 contribute to increased triglyceride levels in Mexicans.

metaCCA: summary statistics-based multivariate meta-analysis of genome-wide association studies using canonical correlation analysis

MOTIVATION

A dominant approach to genetic association studies is to perform univariate tests between genotype-phenotype pairs. However, analyzing related traits together increases statistical power, and certain complex associations become detectable only when several variants are tested jointly. Currently, modest sample sizes of individual cohorts, and restricted availability of individual-level genotype-phenotype data across the cohorts limit conducting multivariate tests.

RESULTS

We introduce metaCCA, a computational framework for summary statistics-based analysis of a single or multiple studies that allows multivariate representation of both genotype and phenotype. It extends the statistical technique of canonical correlation analysis to the setting where original individual-level records are not available, and employs a covariance shrinkage algorithm to achieve robustness.Multivariate meta-analysis of two Finnish studies of nuclear magnetic resonance metabolomics by metaCCA, using standard univariate output from the program SNPTEST, shows an excellent agreement with the pooled individual-level analysis of original data. Motivated by strong multivariate signals in the lipid genes tested, we envision that multivariate association testing using metaCCA has a great potential to provide novel insights from already published summary statistics from high-throughput phenotyping technologies.

AVAILABILITY AND IMPLEMENTATION

Code is available at https://github.com/aalto-ics-kepaco

CONTACTS

anna.cichonska@helsinki.fi or matti.pirinen@helsinki.fi

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Pharmacological Targeting of the Atherogenic Dyslipidemia Complex: The Next Frontier in CVD Prevention Beyond Lowering LDL Cholesterol

Notwithstanding the effectiveness of lowering LDL cholesterol, residual CVD risk remains in high-risk populations, including patients with diabetes, likely contributed to by non-LDL lipid abnormalities. In this Perspectives in Diabetes article, we emphasize that changing demographics and lifestyles over the past few decades have resulted in an epidemic of the "atherogenic dyslipidemia complex," the main features of which include hypertriglyceridemia, low HDL cholesterol levels, qualitative changes in LDL particles, accumulation of remnant lipoproteins, and postprandial hyperlipidemia. We briefly review the underlying pathophysiology of this form of dyslipidemia, in particular its association with insulin resistance, obesity, and type 2 diabetes, and the marked atherogenicity of this condition. We explain the failure of existing classes of therapeutic agents such as fibrates, niacin, and cholesteryl ester transfer protein inhibitors that are known to modify components of the atherogenic dyslipidemia complex. Finally, we discuss targeted repurposing of existing therapies and review promising new therapeutic strategies to modify the atherogenic dyslipidemia complex. We postulate that targeting the central abnormality of the atherogenic dyslipidemia complex, the elevation of triglyceride-rich lipoprotein particles, represents a new frontier in CVD prevention and is likely to prove the most effective strategy in correcting most aspects of the atherogenic dyslipidemia complex, thereby preventing CVD events.

Employing biomarkers of healthy ageing for leveraging genetic studies into human longevity

Genetic studies have thus far identified a limited number of loci associated with human longevity by applying age at death or survival up to advanced ages as phenotype. As an alternative approach, one could first try to identify biomarkers of healthy ageing and the genetic variants associated with these traits and subsequently determine the association of these variants with human longevity. In the present study, we used this approach by testing whether the 35 baseline serum parameters measured in the Leiden Longevity Study (LLS) meet the proposed criteria for a biomarker of healthy ageing. The LLS consists of 421 families with long-lived siblings of European descent, who were recruited together with their offspring and the spouses of the offspring (controls). To test the four criteria for a biomarker of healthy ageing in the LLS, we determined the association of the serum parameters with chronological age, familial longevity, general practitioner-reported general health, and mortality. Out of the 35 serum parameters, we identified glucose, insulin, and triglycerides as biomarkers of healthy ageing, meeting all four criteria in the LLS. We subsequently showed that the genetic variants previously associated with these parameters are significantly enriched in the largest genome-wide association study for human longevity. In conclusion, we showed that biomarkers of healthy ageing can be used to leverage genetic studies into human longevity. We identified several genetic variants influencing the variation in glucose, insulin and triglycerides that contribute to human longevity.

From Loci to Biology: Functional Genomics of Genome-Wide Association for Coronary Disease

Genome-wide association studies have provided a rich collection of ≈ 58 coronary artery disease (CAD) loci that suggest the existence of previously unsuspected new biology relevant to atherosclerosis. However, these studies only identify genomic loci associated with CAD, and many questions remain even after a genomic locus is definitively implicated, including the nature of the causal variant(s) and the causal gene(s), as well as the directionality of effect. There are several tools that can be used for investigation of the functional genomics of these loci, and progress has been made on a limited number of novel CAD loci. New biology regarding atherosclerosis and CAD will be learned through the functional genomics of these loci, and the hope is that at least some of these new pathways relevant to CAD pathogenesis will yield new therapeutic targets for the prevention and treatment of CAD.

Surprises From Genetic Analyses of Lipid Risk Factors for Atherosclerosis

Observational epidemiological studies have associated plasma lipid concentrations with risk for coronary heart disease (CHD), but these studies cannot distinguish cause from mere correlation. Human genetic studies, when considered with the results of randomized controlled trials of medications, can potentially shed light on whether lipid biomarkers are causal for diseases. Genetic analyses and randomized trials suggest that low-density lipoprotein is causal for CHD, whereas high-density lipoprotein is not. Surprisingly, human genetic evidence suggests that lipoprotein(a) and triglyceride-rich lipoproteins causally contribute to CHD. Gene variants leading to higher levels of plasma apolipoprotein B-containing lipoproteins [low-density lipoprotein, triglyceride-rich lipoproteins, or lipoprotein(a)] consistently increase risk for CHD. For triglyceride-rich lipoproteins, the most compelling evidence revolves around lipoprotein lipase and its endogenous facilitator (APOA5 [apolipoprotein A-V]) and inhibitory proteins (APOC3 [apolipoprotein C-III], ANGPTL4 [angiopoietin like 4]). Combined, these genetic results anticipate that, beyond low-density lipoprotein, pharmacological lowering of triglyceride-rich lipoproteins or lipoprotein(a) will reduce risk for CHD, but this remains to be proven through randomized controlled trials.

Blood lipids influence DNA methylation in circulating cells

Background

Cells can be primed by external stimuli to obtain a long-term epigenetic memory. We hypothesize that long-term exposure to elevated blood lipids can prime circulating immune cells through changes in DNA methylation, a process that may contribute to the development of atherosclerosis. To interrogate the causal relationship between triglyceride, low-density lipoprotein (LDL) cholesterol, and high-density lipoprotein (HDL) cholesterol levels and genome-wide DNA methylation while excluding confounding and pleiotropy, we perform a stepwise Mendelian randomization analysis in whole blood of 3296 individuals.

Results

This analysis shows that differential methylation is the consequence of inter-individual variation in blood lipid levels and not vice versa. Specifically, we observe an effect of triglycerides on DNA methylation at three CpGs, of LDL cholesterol at one CpG, and of HDL cholesterol at two CpGs using multivariable Mendelian randomization. Using RNA-seq data available for a large subset of individuals (N = 2044), DNA methylation of these six CpGs is associated with the expression of CPT1A and SREBF1 (for triglycerides), DHCR24 (for LDL cholesterol) and ABCG1 (for HDL cholesterol), which are all key regulators of lipid metabolism.

Conclusions

Our analysis suggests a role for epigenetic priming in end-product feedback control of lipid metabolism and highlights Mendelian randomization as an effective tool to infer causal relationships in integrative genomics data.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-016-1000-6) contains supplementary material, which is available to authorized users.

Selecting instruments for Mendelian randomization in the wake of genome-wide association studies

Mendelian randomization (MR) studies typically assess the pathogenic relevance of environmental exposures or disease biomarkers, using genetic variants that instrument these exposures. The approach is gaining popularity-our systematic review reveals a greater than 10-fold increase in MR studies published between 2004 and 2015. When the MR paradigm was first proposed, few biomarker- or exposure-related genetic variants were known, most having been identified by candidate gene studies. However, genome-wide association studies (GWAS) are now providing a rich source of potential instruments for MR analysis. Many early reviews covering the concept, applications and analytical aspects of the MR technique preceded the surge in GWAS, and thus the question of how best to select instruments for MR studies from the now extensive pool of available variants has received insufficient attention. Here we focus on the most common category of MR studies-those concerning disease biomarkers. We consider how the selection of instruments for MR analysis from GWAS requires consideration of: the assumptions underlying the MR approach; the biology of the biomarker; the genome-wide distribution, frequency and effect size of biomarker-associated variants (the genetic architecture); and the specificity of the genetic associations. Based on this, we develop guidance that may help investigators to plan and readers interpret MR studies.

Genetic markers of inflammation may not contribute to metabolic traits in Mexican children

Background: Low-grade chronic inflammation is a common feature of obesity and its cardio-metabolic complications. However, little is known about a possible causal role of inflammation in metabolic disorders. Mexico is among the countries with the highest obesity rates in the world and the admixed Mexican population is a relevant sample due to high levels of genetic diversity.

Methods: Here, we studied 1,462 Mexican children recruited from Mexico City. Six genetic variants in five inflammation-related genes were genotyped: rs1137101 (leptin receptor (LEPR)), rs7305618 (hepatocyte nuclear factor 1 alpha (HNF1A)), rs1800629 (tumor necrosis factor alpha (TNFA)), rs1800896, rs1800871 (interleukin-10 (IL-10)), rs1862513 (resistin (RETN)). Ten continuous and eight binary traits were assessed. Linear and logistic regression models were used adjusting for age, sex, and recruitment centre.

Results: We found that one SNP displayed a nominal evidence of association with a continuous trait: rs1800871 (IL-10) with LDL (beta = −0.068 ± 1.006, P = 0.01). Subsequently, we found one nominal association with a binary trait: rs7305618 (HNF1A) with family history of hypertension (odds-ratio = 1.389 [1.054–1.829], P = 0.02). However, no P-value passed the Bonferroni correction for multiple testing.

Discussion: Our data in a Mexican children population are consistent with previous reports in European adults in failing to demonstrate an association between inflammation-associated single nucleotide polymorphisms (SNPs) and metabolic traits.

A Genome-Wide mQTL Analysis in Human Adipose Tissue Identifies Genetic Variants Associated with DNA Methylation, Gene Expression and Metabolic Traits

Little is known about the extent to which interactions between genetics and epigenetics may affect the risk of complex metabolic diseases and/or their intermediary phenotypes. We performed a genome-wide DNA methylation quantitative trait locus (mQTL) analysis in human adipose tissue of 119 men, where 592,794 single nucleotide polymorphisms (SNPs) were related to DNA methylation of 477,891 CpG sites, covering 99% of RefSeq genes. SNPs in significant mQTLs were further related to gene expression in adipose tissue and obesity related traits. We found 101,911 SNP-CpG pairs (mQTLs) in cis and 5,342 SNP-CpG pairs in trans showing significant associations between genotype and DNA methylation in adipose tissue after correction for multiple testing, where cis is defined as distance less than 500 kb between a SNP and CpG site. These mQTLs include reported obesity, lipid and type 2 diabetes loci, e.g. ADCY3/POMC, APOA5, CETP, FADS2, GCKR, SORT1 and LEPR. Significant mQTLs were overrepresented in intergenic regions meanwhile underrepresented in promoter regions and CpG islands. We further identified 635 SNPs in significant cis-mQTLs associated with expression of 86 genes in adipose tissue including CHRNA5, G6PC2, GPX7, RPL27A, THNSL2 and ZFP57. SNPs in significant mQTLs were also associated with body mass index (BMI), lipid traits and glucose and insulin levels in our study cohort and public available consortia data. Importantly, the Causal Inference Test (CIT) demonstrates how genetic variants mediate their effects on metabolic traits (e.g. BMI, cholesterol, high-density lipoprotein (HDL), hemoglobin A1c (HbA1c) and homeostatic model assessment of insulin resistance (HOMA-IR)) via altered DNA methylation in human adipose tissue. This study identifies genome-wide interactions between genetic and epigenetic variation in both cis and trans positions influencing gene expression in adipose tissue and in vivo (dys)metabolic traits associated with the development of obesity and diabetes.

Systems proteomics of liver mitochondria function

Recent improvements in quantitative proteomics approaches, including Sequential Window Acquisition of all Theoretical Mass Spectra (SWATH-MS), permit reproducible large-scale protein measurements across diverse cohorts. Together with genomics, transcriptomics, and other technologies, transomic data sets can be generated that permit detailed analyses across broad molecular interaction networks. Here, we examine mitochondrial links to liver metabolism through the genome, transcriptome, proteome, and metabolome of 386 individuals in the BXD mouse reference population. Several links were validated between genetic variants toward transcripts, proteins, metabolites, and phenotypes. Among these, sequence variants in Cox7a2l alter its protein's activity, which in turn leads to downstream differences in mitochondrial supercomplex formation. This data set demonstrates that the proteome can now be quantified comprehensively, serving as a key complement to transcriptomics, genomics, and metabolomics--a combination moving us forward in complex trait analysis.

Bivariate genome-wide association study identifies novel pleiotropic loci for lipids and inflammation

BACKGROUND

Genome-wide association studies (GWAS) have identified multiple genetic loci for C-reactive protein (CRP) and lipids, of which some overlap. We aimed to identify genetic pleiotropy among CRP and lipids in order to better understand the shared biology of chronic inflammation and lipid metabolism.

RESULTS

In a bivariate GWAS, we combined summary statistics of published GWAS on CRP (n = 66,185) and lipids, including LDL-cholesterol, HDL-cholesterol, triglycerides, and total cholesterol (n = 100,184), using an empirical weighted linear-combined test statistic. We sought replication for novel CRP associations in an independent sample of 17,743 genotyped individuals, and performed in silico replication of novel lipid variants in 93,982 individuals. Fifty potentially pleiotropic SNPs were identified among CRP and lipids: 21 for LDL-cholesterol and CRP, 20 for HDL-cholesterol and CRP, 21 for triglycerides, and CRP and 20 for total cholesterol and CRP. We identified and significantly replicated three novel SNPs for CRP in or near CTSB/FDFT1 (rs10435719, Preplication: 2.6 × 10(-5)), STAG1/PCCB (rs7621025, Preplication: 1.4 × 10(-3)) and FTO (rs1558902, Preplication: 2.7 × 10(-5)). Seven pleiotropic lipid loci were replicated in the independent set of MetaboChip samples of the Global Lipids Genetics Consortium. Annotating the effect of replicated CRP SNPs to the expression of nearby genes, we observed an effect of rs10435719 on gene expression of FDFT1, and an effect of rs7621025 on PCCB.

CONCLUSIONS

Our large scale combined GWAS analysis identified numerous pleiotropic loci for CRP and lipids providing further insight in the genetic interrelation between lipids and inflammation. In addition, we provide evidence for FDFT1, PCCB and FTO to be associated with CRP levels.

Replication and Characterization of Association between ABO SNPs and Red Blood Cell Traits by Meta-Analysis in Europeans

Red blood cell (RBC) traits are routinely measured in clinical practice as important markers of health. Deviations from the physiological ranges are usually a sign of disease, although variation between healthy individuals also occurs, at least partly due to genetic factors. Recent large scale genetic studies identified loci associated with one or more of these traits; further characterization of known loci and identification of new loci is necessary to better understand their role in health and disease and to identify potential molecular mechanisms. We performed meta-analysis of Metabochip association results for six RBC traits—hemoglobin concentration (Hb), hematocrit (Hct), mean corpuscular hemoglobin (MCH), mean corpuscular hemoglobin concentration (MCHC), mean corpuscular volume (MCV) and red blood cell count (RCC)—in 11 093 Europeans from seven studies of the UCL-LSHTM-Edinburgh-Bristol (UCLEB) Consortium. We identified 394 non-overlapping SNPs in five loci at genome-wide significance: 6p22.1-6p21.33 (with HFE among others), 6q23.2 (with HBS1L among others), 6q23.3 (contains no genes), 9q34.3 (only ABO gene) and 22q13.1 (with TMPRSS6 among others), replicating previous findings of association with RBC traits at these loci and extending them by imputation to 1000 Genomes. We further characterized associations between ABO SNPs and three traits: hemoglobin, hematocrit and red blood cell count, replicating them in an independent cohort. Conditional analyses indicated the independent association of each of these traits with ABO SNPs and a role for blood group O in mediating the association. The 15 most significant RBC-associated ABO SNPs were also associated with five cardiometabolic traits, with discordance in the direction of effect between groups of traits, suggesting that ABO may act through more than one mechanism to influence cardiometabolic risk.

Genetics of Lipid and Lipoprotein Disorders and Traits

PURPOSE OF REVIEW

Plasma lipids, namely cholesterol and triglyceride, and lipoproteins, such as low-density lipoprotein (LDL) and high-density lipoprotein, serve numerous physiological roles. Perturbed levels of these traits underlie monogenic dyslipidemias, a diverse group of multisystem disorders. We are on the verge of having a relatively complete picture of the human dyslipidemias and their components.

RECENT FINDINGS

Recent advances in genetics of plasma lipids and lipoproteins include the following: (1) expanding the range of genes causing monogenic dyslipidemias, particularly elevated LDL cholesterol; (2) appreciating the role of polygenic effects in such traits as familial hypercholesterolemia and combined hyperlipidemia; (3) accumulating a list of common variants that determine plasma lipids and lipoproteins; (4) applying exome sequencing to identify collections of rare variants determining plasma lipids and lipoproteins that via Mendelian randomization have also implicated gene products such as NPC1L1, APOC3, LDLR, APOA5, and ANGPTL4 as causal for atherosclerotic cardiovascular disease; and (5) using naturally occurring genetic variation to identify new drug targets, including inhibitors of apolipoprotein (apo) C-III, apo(a), ANGPTL3, and ANGPTL4.

SUMMARY

Here, we compile this disparate range of data linking human genetic variation to plasma lipids and lipoproteins, providing a "one stop shop" for the interested reader.

Distinct Patterns of Association of Variants at 11q23.3 Chromosomal Region with Coronary Artery Disease and Dyslipidemia in the Population of Andhra Pradesh, India

In our attempt to comprehensively understand the nature of association of variants at 11q23.3 apolipoprotein gene cluster region, we genotyped a prioritized set of 96 informative SNPs using Fluidigm customized SNP genotyping platform in a sample of 508 coronary artery disease (CAD) cases and 516 controls. We found 12 SNPs as significantly associated with CAD at P <0.05, albeit only four (rs2849165, rs17440396, rs6589566 and rs633389) of these remained significant after Benjamin Hochberg correction. Of the four, while rs6589566 confers risk to CAD, the other three SNPs reduce risk for the disease. Interaction of variants that belong to regulatory genes BUD13 and ZPR1 with APOA5-APOA4 intergenic variants is also observed to significantly increase the risk towards CAD. Further, ROC analysis of the risk scores of the 12 significant SNPs suggests that our study has substantial power to confer these genetic variants as predictors of risk for CAD, as illustrated by AUC (0.763; 95% CI: 0.729-0.798, p = <0.0001). On the other hand, the protective SNPs of CAD are associated with elevated Low Density Lipoprotein Cholesterol and Total Cholesterol levels, hence with dyslipidemia, in our sample of controls, which may suggest distinct effects of the variants at 11q23.3 chromosomal region towards CAD and dyslipidemia. It may be necessary to replicate these findings in the independent and ethnically heterogeneous Indian samples in order to establish this as an Indian pattern. However, only functional analysis of the significant variants identified in our study can provide more precise understanding of the mechanisms involved in the contrasting nature of their effects in manifesting dyslipidemia and CAD.

Variant ASGR1 Associated with a Reduced Risk of Coronary Artery Disease

BACKGROUND

Several sequence variants are known to have effects on serum levels of non-high-density lipoprotein (HDL) cholesterol that alter the risk of coronary artery disease.

METHODS

We sequenced the genomes of 2636 Icelanders and found variants that we then imputed into the genomes of approximately 398,000 Icelanders. We tested for association between these imputed variants and non-HDL cholesterol levels in 119,146 samples. We then performed replication testing in two populations of European descent. We assessed the effects of an implicated loss-of-function variant on the risk of coronary artery disease in 42,524 case patients and 249,414 controls from five European ancestry populations. An augmented set of genomes was screened for additional loss-of-function variants in a target gene. We evaluated the effect of an implicated variant on protein stability.

RESULTS

We found a rare noncoding 12-base-pair (bp) deletion (del12) in intron 4 of ASGR1, which encodes a subunit of the asialoglycoprotein receptor, a lectin that plays a role in the homeostasis of circulating glycoproteins. The del12 mutation activates a cryptic splice site, leading to a frameshift mutation and a premature stop codon that renders a truncated protein prone to degradation. Heterozygous carriers of the mutation (1 in 120 persons in our study population) had a lower level of non-HDL cholesterol than noncarriers, a difference of 15.3 mg per deciliter (0.40 mmol per liter) (P=1.0×10(-16)), and a lower risk of coronary artery disease (by 34%; 95% confidence interval, 21 to 45; P=4.0×10(-6)). In a larger set of sequenced samples from Icelanders, we found another loss-of-function ASGR1 variant (p.W158X, carried by 1 in 1850 persons) that was also associated with lower levels of non-HDL cholesterol (P=1.8×10(-3)).

CONCLUSIONS

ASGR1 haploinsufficiency was associated with reduced levels of non-HDL cholesterol and a reduced risk of coronary artery disease. (Funded by the National Institutes of Health and others.).

Metabolic Characterization of a Rare Genetic Variation Within APOC3 and Its Lipoprotein Lipase-Independent Effects

BACKGROUND

Plasma triglyceride levels have been implicated in atherosclerosis and coronary heart disease. Apolipoprotein C-III (APOC3) plays a key role in the hydrolysis of triglyceride-rich lipoproteins to remnant particles by lipoprotein lipase (LPL) and their uptake by the liver. A rare variant in APOC3(rs138326449) has been associated with triglyceride, very low-density lipoprotein, and high-density lipoprotein levels, as well as risk of coronary heart disease. We aimed to characterize the impact of this locus across a broad set of mainly lipids-focused metabolic measures.

METHODS AND RESULTS

A high-throughput serum nuclear magnetic resonance metabolomics platform was used to quantify 225 metabolic measures in 13 285 participants from 2 European population cohorts. We analyzed the effect of the APOC3 variant on the metabolic measures and used the common LPL(rs12678919) polymorphism to test for LPL-independent effects. Eighty-one metabolic measures showed evidence of association with APOC3(rs138326449). In addition to previously reported triglyceride and high-density lipoprotein associations, the variant was also associated with very low-density lipoprotein and high-density lipoprotein composition measures, other cholesterol measures, and fatty acids. Comparison of the APOC3 and LPL associations revealed that APOC3 association results for medium and very large very low-density lipoprotein composition are unlikely to be solely predictable by the action of APOC3 through LPL.

CONCLUSIONS

We characterized the effects of the rare APOC3(rs138326449) loss of function mutation in lipoprotein metabolism, as well as the effects of LPL(rs12678919). Our results improve our understanding of the role of APOC3 in triglyceride metabolism, its LPL independent action, and the complex and correlated nature of human metabolites.

The multifaceted interplay between lipids and epigenetics

PURPOSE OF REVIEW

The interplay between lipids and epigenetic mechanisms has recently gained increased interest because of its relevance for common diseases and most notably atherosclerosis. This review discusses recent advances in unravelling this interplay with a particular focus on promising approaches and methods that will be able to establish causal relationships.

RECENT FINDINGS

Complementary approaches uncovered close links between circulating lipids and epigenetic mechanisms at multiple levels. A characterization of lipid-associated genetic variants suggests that these variants exert their influence on lipid levels through epigenetic changes in the liver. Moreover, exposure of monocytes to lipids persistently alters their epigenetic makeup resulting in more proinflammatory cells. Hence, epigenetic changes can both impact on and be induced by lipids.

SUMMARY

It is the combined application of technological advances to probe epigenetic modifications at a genome-wide scale and methodological advances aimed at causal inference (including Mendelian randomization and integrative genomics) that will elucidate the interplay between circulating lipids and epigenetics. Understanding its role in the development of atherosclerosis holds the promise of identifying a new category of therapeutic targets, since epigenetic changes are amenable to reversal.

Plasma levels of vitamin K and the risk of ischemic heart disease: a Mendelian randomization study

Essentials Vitamin K plays a role in coagulation, and deficiency may promote coronary artery calcification. The role of vitamin K1 in heart disease was assessed using Mendelian randomization in Caucasians. Genetically higher vitamin K1 was associated with a higher risk of ischemic heart disease. Further research elucidating the role of vitamin K1 in ischemic heart disease could be useful.

SUMMARY

Background Vitamin K1 is a nutrient in green leafy vegetables; deficiency may promote coronary artery calcification. Warfarin, an anticoagulant used in secondary prevention of thrombotic events, is a vitamin K antagonist. Thrombotic and coronary events may share risk factors. Objectives To clarify the role of vitamin K1 in ischemic heart disease, the risk of coronary artery disease/myocardial infarction (CAD/MI) was assessed according to genetically determined vitamin K1 levels. Given vitamin K1 is fat soluble, associations with lipids were similarly assessed to assess pleotropic effects via lipids. Methods Separate sample instrumental variable analysis with genetic instruments (Mendelian randomization) was used to obtain an unconfounded estimate of the association of vitamin K1 (based on rs2108622 [CYP4F2], rs4645543 [KCNK9] and rs2192574 [CTNNA2] from a genome-wide association study) with CAD/MI using CARDIoGRAMplusC4D (cases = 64 374; controls = 130 681) and with lipids using Global Lipids Genetics Consortium Results (n = 196 475). Results Vitamin K1 single nucleotide polymorphisms were positively associated with CAD/MI (odds ratio [OR], 1.17 per unit [nmol L(-1) ] of natural log-transformed genetically predicted vitamin K1 ; 95% confidence interval [CI], 1.08-1.26), but not with inverse normal transformed low-density lipoprotein cholesterol (-0.0003; 95% CI, -0.03 to 0.03), high-density lipoprotein cholesterol (0.02; 95% CI, -0.01 to 0.05) or triglycerides (-0.01; 95% CI, -0.04 to 0.02). Considering only rs2108622, which is functionally relevant to vitamin K1 , the association for CAD/MI was stronger (OR, 1.21; 95% CI, 1.08-1.36). Conclusions Vitamin K may cause CAD/MI; whether vitamin K or other determinants of coagulation could be relevant to primary prevention might be worth considering.

Blood lipids and prostate cancer: a Mendelian randomization analysis

Genetic risk scores were used as unconfounded instruments for specific lipid traits (Mendelian randomization) to assess whether circulating lipids causally influence prostate cancer risk. Data from 22,249 prostate cancer cases and 22,133 controls from 22 studies within the international PRACTICAL consortium were analyzed. Allele scores based on single nucleotide polymorphisms (SNPs) previously reported to be uniquely associated with each of low-density lipoprotein (LDL), high-density lipoprotein (HDL), and triglyceride (TG) levels, were first validated in an independent dataset, and then entered into logistic regression models to estimate the presence (and direction) of any causal effect of each lipid trait on prostate cancer risk. There was weak evidence for an association between the LDL genetic score and cancer grade: the odds ratio (OR) per genetically instrumented standard deviation (SD) in LDL, comparing high- (≥7 Gleason score) versus low-grade (<7 Gleason score) cancers was 1.50 (95% CI: 0.92, 2.46; P = 0.11). A genetically instrumented SD increase in TGs was weakly associated with stage: the OR for advanced versus localized cancer per unit increase in genetic risk score was 1.68 (95% CI: 0.95, 3.00; P = 0.08). The rs12916-T variant in 3-hydroxy-3-methylglutaryl-CoA reductase (HMGCR) was inversely associated with prostate cancer (OR: 0.97; 95% CI: 0.94, 1.00; P = 0.03). In conclusion, circulating lipids, instrumented by our genetic risk scores, did not appear to alter prostate cancer risk. We found weak evidence that higher LDL and TG levels increase aggressive prostate cancer risk, and that a variant in HMGCR (that mimics the LDL lowering effect of statin drugs) reduces risk. However, inferences are limited by sample size and evidence of pleiotropy.

Identifying genetic loci affecting antidepressant drug response in depression using drug-gene interaction models

Antidepressants are often only moderately successful in decreasing the severity of depressive symptoms. In part, antidepressant treatment response in patients with depression is genetically determined. However, although a large number of studies have been conducted aiming to identify genetic variants associated with antidepressant drug response in depression, only a few variants have been repeatedly identified. Within the present review, we will discuss the methodological challenges and limitations of the studies that have been conducted on this topic to date (e.g., 'treated-only design', statistical power) and we will discuss how specifically drug-gene interaction models can be used to be better able to identify genetic variants associated with antidepressant drug response in depression.

Role of Adiponectin in Coronary Heart Disease Risk: A Mendelian Randomization Study

Supplemental Digital Content is available in the text.

Rationale:

Hypoadiponectinemia correlates with several coronary heart disease (CHD) risk factors. However, it is unknown whether adiponectin is causally implicated in CHD pathogenesis.

Objective:

We aimed to investigate the causal effect of adiponectin on CHD risk.

Methods and Results:

We undertook a Mendelian randomization study using data from genome-wide association studies consortia. We used the ADIPOGen consortium to identify genetic variants that could be used as instrumental variables for the effect of adiponectin. Data on the association of these genetic variants with CHD risk were obtained from CARDIoGRAM (22 233 CHD cases and 64 762 controls of European ancestry) and from CARDIoGRAMplusC4D Metabochip (63 746 cases and 130 681 controls; ≈ 91% of European ancestry) consortia. Data on the association of genetic variants with adiponectin levels and with CHD were combined to estimate the influence of blood adiponectin on CHD risk. In the conservative approach (restricted to using variants within the adiponectin gene as instrumental variables), each 1 U increase in log blood adiponectin concentration was associated with an odds ratio for CHD of 0.83 (95% confidence interval, 0.68–1.01) in CARDIoGRAM and 0.97 (95% confidence interval, 0.84–1.12) in CARDIoGRAMplusC4D Metabochip. Findings from the liberal approach (including variants in any locus across the genome) indicated a protective effect of adiponectin that was attenuated to the null after adjustment for known CHD predictors.

Conclusions:

Overall, our findings do not support a causal role of adiponectin levels in CHD pathogenesis.

Pharmacogenomics in diabetes mellitus: insights into drug action and drug discovery

Genomic studies have greatly advanced our understanding of the multifactorial aetiology of type 2 diabetes mellitus (T2DM) as well as the multiple subtypes of monogenic diabetes mellitus. In this Review, we discuss the existing pharmacogenetic evidence in both monogenic diabetes mellitus and T2DM. We highlight mechanistic insights from the study of adverse effects and the efficacy of antidiabetic drugs. The identification of extreme sulfonylurea sensitivity in patients with diabetes mellitus owing to heterozygous mutations in HNF1A represents a clear example of how pharmacogenetics can direct patient care. However, pharmacogenomic studies of response to antidiabetic drugs in T2DM has yet to be translated into clinical practice, although some moderate genetic effects have now been described that merit follow-up in trials in which patients are selected according to genotype. We also discuss how future pharmacogenomic findings could provide insights into treatment response in diabetes mellitus that, in addition to other areas of human genetics, facilitates drug discovery and drug development for T2DM.

Gut microbiome and lipid metabolism: from associations to mechanisms

PURPOSE OF REVIEW

The gut microbiome has now been convincingly linked to human metabolic health but the underlying causality and mechanisms remain poorly understood. This review focuses on the recent progress in establishing the associations between gut microbiome species and lipid metabolism in humans and discusses how to move from association toward causality and mechanistic understanding, which is essential knowledge to bring the observed associations into clinical use.

RECENT FINDINGS

Recent population-based association studies have shown that the gut microbiota composition can explain a substantial proportion of the interindividual variation in blood triglycerides and HDL-cholesterol level and predict metabolic response to diet and drug. Faecal transplantation has suggested that this is a causal effect of microbiome on host metabolism, although the underlying mechanism remains largely unexplored.

SUMMARY

The gut microbiome is transitioning from being a 'missing' organ to a potential target for therapeutic applications. Due to the complex interplay between the gut microbiome, the host genome, and diet, a systematic approach is required to pave the way for therapeutic development.

Determining Which Phenotypes Underlie a Pleiotropic Signal

Discovering pleiotropic loci is important to understand the biological basis of seemingly distinct phenotypes. Most methods for assessing pleiotropy only test for the overall association between genetic variants and multiple phenotypes. To determine which specific traits are pleiotropic, we evaluate via simulation and application three different strategies. The first is model selection techniques based on the inverse regression of genotype on phenotypes. The second is a subset-based meta analysis ASSET [Bhattacharjee et al., ], which provides an optimal subset of nonnull traits. And the third is a modified Benjamini-Hochberg (B-H) procedure of controlling the expected false discovery rate [Benjamini and Hochberg, ] in the framework of phenome-wide association study. From our simulations we see that an inverse regression-based approach MultiPhen [O'Reilly et al., ] is more powerful than ASSET for detecting overall pleiotropic association, except for when all the phenotypes are associated and have genetic effects in the same direction. For determining which specific traits are pleiotropic, the modified B-H procedure performs consistently better than the other two methods. The inverse regression-based selection methods perform competitively with the modified B-H procedure only when the phenotypes are weakly correlated. The efficiency of ASSET is observed to lie below and in between the efficiency of the other two methods when the traits are weakly and strongly correlated, respectively. In our application to a large GWAS, we find that the modified B-H procedure also performs well, indicating that this may be an optimal approach for determining the traits underlying a pleiotropic signal.

The Contribution of GWAS Loci in Familial Dyslipidemias

Familial combined hyperlipidemia (FCH) is a complex and common familial dyslipidemia characterized by elevated total cholesterol and/or triglyceride levels with over five-fold risk of coronary heart disease. The genetic architecture and contribution of rare Mendelian and common variants to FCH susceptibility is unknown. In 53 Finnish FCH families, we genotyped and imputed nine million variants in 715 family members with DNA available. We studied the enrichment of variants previously implicated with monogenic dyslipidemias and/or lipid levels in the general population by comparing allele frequencies between the FCH families and population samples. We also constructed weighted polygenic scores using 212 lipid-associated SNPs and estimated the relative contributions of Mendelian variants and polygenic scores to the risk of FCH in the families. We identified, across the whole allele frequency spectrum, an enrichment of variants known to elevate, and a deficiency of variants known to lower LDL-C and/or TG levels among both probands and affected FCH individuals. The score based on TG associated SNPs was particularly high among affected individuals compared to non-affected family members. Out of 234 affected FCH individuals across the families, seven (3%) carried Mendelian variants and 83 (35%) showed high accumulation of either known LDL-C or TG elevating variants by having either polygenic score over the 90^th percentile in the population. The positive predictive value of high score was much higher for affected FCH individuals than for similar sporadic cases in the population. FCH is highly polygenic, supporting the hypothesis that variants across the whole allele frequency spectrum contribute to this complex familial trait. Polygenic SNP panels improve identification of individuals affected with FCH, but their clinical utility remains to be defined.

Familial combined hyperlipidemia (FCH) is a familial dyslipidemia and the most common familial risk factor for premature coronary heart disease. Its genetic architecture is poorly understood. Rare high-impact variants have been identified in some patients, but have not explained a substantial portion of the trait. FCH has previously been speculated to be a polygenic disorder, but genetic data supporting this hypothesis have so far been incomplete. We provide experimental evidence for the polygenicity and heterogeneity of FCH in a large set of affected families using comprehensive genome-wide variant data. Approximately a third of the affected FCH individuals in our sample had high polygenic burden, and only a minority carried high-impact variants identifiable by genotyping. We show that the polygenic burden of affected FCH family members is comparable to that observed in individuals with similar lipid phenotypes in the general population. Genetic variants identified in large-scale population studies can also underlie the typical phenotypes observed in complex familial diseases such as FCH. Advances in genetic diagnosis based on population samples may thus also benefit FCH families. Families without high polygenic burden are good candidates for sequencing studies to identify rare variants not observable with genotyping.

SUGP1 is a novel regulator of cholesterol metabolism

A large haplotype on chromosome 19p13.11 tagged by rs10401969 in intron 8 of SURP and G patch domain containing 1 (SUGP1) is associated with coronary artery disease (CAD), plasma LDL cholesterol levels, and other energy metabolism phenotypes. Recent studies have suggested that TM6SF2 is the causal gene within the locus, but we postulated that this locus could harbor additional CAD risk genes, including the putative splicing factor SUGP1. Indeed, we found that rs10401969 regulates SUGP1 exon 8 skipping, causing non-sense-mediated mRNA decay. Hepatic Sugp1 overexpression in CD1 male mice increased plasma cholesterol levels 20–50%. In human hepatoma cell lines, SUGP1 knockdown stimulated 3-hydroxy-3-methylglutaryl-CoA reductase (HMGCR) alternative splicing and decreased HMGCR transcript stability, thus reducing cholesterol synthesis and increasing LDL uptake. Our results strongly support a role for SUGP1 as a novel regulator of cholesterol metabolism and suggest that it contributes to the relationship between rs10401969 and plasma cholesterol.

Interactive effects of C-reactive protein levels on the association between APOE variants and triglyceride levels in a Taiwanese population

Background

Apolipoprotein E (APOE) plays a major role in lipid metabolism and inflammation. However, the association between APOE gene polymorphisms and serum triglyceride levels remains controversial. We tested the effects of APOE variants on triglyceride levels and their interactions with the inflammatory marker C-reactive protein (CRP) in a Taiwanese population.

Methods

Two APOE single nucleotide polymorphisms (SNPs) rs429358 and rs7412 were genotyped by TaqMan Assay using real time PCR in 595 healthy subjects attending the clinic for routine visits.

Results

After adjustment for clinical covariates, subjects carrying the rs429358-TT genotype and non-ε4 alleles were found to have higher CRP levels, whereas those with rs7412-CC genotype and non-ε2 alleles had significantly higher total and low-density lipoprotein cholesterol levels (all P < 0.01). Using subgroup and interaction analyses, we observed significantly lower triglyceride levels in subjects carrying the rs429358-TT genotype and non-ε4 alleles in the low CRP group (P = 2.71× 10⁻⁴ and P = 4.32 × 10⁻⁴, respectively), but not in those in the high CRP group (interaction P = 0.013 and 0.045, respectively). In addition, multivariate stepwise linear regression analysis showed that subjects carrying the rs429358-TT genotype and non-ε4 alleles with low CRP levels had significantly lower triglyceride levels (P < 0.001 and P < 0.001, respectively). In addition, when combined with the risk alleles of GCKR, APOA5 and LPL gene variants, we observed that triglyceride levels increased significantly with the number of risk alleles (P = 2.9 × 10⁻¹²).

Conclusions

The combination of SNPs and ε alleles at the APOE locus is involved in managing lipid and CRP levels in the Taiwanese population. APOE polymorphisms interact with CRP to regulate triglyceride levels, thus triglyceride concentration is influenced by both the genetic background of the APOE locus and the inflammatory status of a subject.

Electronic supplementary material

The online version of this article (doi:10.1186/s12944-016-0262-z) contains supplementary material, which is available to authorized users.

A perspective for sequencing familial hypercholesterolaemia in African Americans

African Americans suffer disproportionately from poor cardiovascular health outcomes despite similar proportions of African Americans and Americans of European ancestry experiencing elevated cholesterol levels. Some of the variation in cardiovascular outcomes is due to confounding effects of other risk factors, such as hypertension and genetic influence. However, genetic variants found to contribute to variation in serum cholesterol levels in populations of European ancestry are less likely to replicate in populations of African ancestry. To date, there has been limited follow-up on variant discrepancies or on identifying variants that exist in populations of African ancestry. African and African-American populations have the highest levels of genetic heterogeneity, which is a factor that must be considered when evaluating genetic variants in the burgeoning era of personalised medicine. Many of the large published studies identifying genetic variants associated with disease risk have evaluated populations of mostly European ancestry and estimated risk in other populations based on these findings. The purpose of this paper is to provide a perspective, using familial hypercholesterolaemia as an exemplar, that studies evaluating genetic variation focused within minority populations are necessary to identify factors that contribute to disparities in health outcomes and realise the full utility of personalised medicine.

Discovery and validation of sub-threshold genome-wide association study loci using epigenomic signatures

Genetic variants identified by genome-wide association studies explain only a modest proportion of heritability, suggesting that meaningful associations lie 'hidden' below current thresholds. Here, we integrate information from association studies with epigenomic maps to demonstrate that enhancers significantly overlap known loci associated with the cardiac QT interval and QRS duration. We apply functional criteria to identify loci associated with QT interval that do not meet genome-wide significance and are missed by existing studies. We demonstrate that these 'sub-threshold' signals represent novel loci, and that epigenomic maps are effective at discriminating true biological signals from noise. We experimentally validate the molecular, gene-regulatory, cellular and organismal phenotypes of these sub-threshold loci, demonstrating that most sub-threshold loci have regulatory consequences and that genetic perturbation of nearby genes causes cardiac phenotypes in mouse. Our work provides a general approach for improving the detection of novel loci associated with complex human traits.

Discovery and validation of sub-threshold genome-wide association study loci using epigenomic signatures

Genetic variants identified by genome-wide association studies explain only a modest proportion of heritability, suggesting that meaningful associations lie 'hidden' below current thresholds. Here, we integrate information from association studies with epigenomic maps to demonstrate that enhancers significantly overlap known loci associated with the cardiac QT interval and QRS duration. We apply functional criteria to identify loci associated with QT interval that do not meet genome-wide significance and are missed by existing studies. We demonstrate that these 'sub-threshold' signals represent novel loci, and that epigenomic maps are effective at discriminating true biological signals from noise. We experimentally validate the molecular, gene-regulatory, cellular and organismal phenotypes of these sub-threshold loci, demonstrating that most sub-threshold loci have regulatory consequences and that genetic perturbation of nearby genes causes cardiac phenotypes in mouse. Our work provides a general approach for improving the detection of novel loci associated with complex human traits.

DOI:http://dx.doi.org/10.7554/eLife.10557.001

Most complex traits are governed by a large number of genetic contributors, each playing only a modest effect. This makes it difficult to identify the genetic variants that increase disease risk, hindering the discovery of new drug targets and the development of new therapeutics.

To overcome this limitation in discovery power, the field of human genetics has traditionally sought increasingly large groups, or cohorts, of afflicted and non-afflicted individuals. Studies of large cohorts are a powerful approach for discovering new disease genes, but such groups are often impractical and sometimes impossible to obtain. However, it has become possible to complement the genetic evidence found in disease association studies with biological evidence of the effects of disease-associated genetic variants.

Wang et al. focus specifically on genetic sites, or loci, that do not affect protein sequence but instead affect the non-coding control regions. These are known as enhancer elements, as they can enhance the expression of nearby genes. These loci constitute the majority of disease regions, and thus are extremely important, but their discovery has been hindered by our relatively poor understanding of the human genome.

Chemical modifications known as epigenomic marks are indicative of enhancer regions. By studying the factors that affect heart rhythm, Wang et al. show that specific combinations of epigenomic marks are enriched in known trait-associated regions. This knowledge was then used to prioritize the further investigation of genetic regions that genome-wide association studies had only weakly linked to heart rhythm alterations. Wang et al. directly confirmed that genetic differences in “sub-threshold” regions indeed alter the activity of these regulatory regions in human heart cells. Furthermore, mutating or perturbing the predicted target genes of the sub-threshold enhancers caused heart defects in mouse and zebrafish.

Wang et al. have demonstrated that epigenome maps can help to distinguish which sub-threshold regions from genome-wide association studies are more likely to contribute to a disease. This allows for the discovery of new disease genes with much smaller cohorts than would be needed otherwise, thus speeding up the development of new therapeutics by many years.

DOI:http://dx.doi.org/10.7554/eLife.10557.002

Very low-depth sequencing in a founder population identifies a cardioprotective APOC3 signal missed by genome-wide imputation

Cohort-wide very low-depth whole-genome sequencing (WGS) can comprehensively capture low-frequency sequence variation for the cost of a dense genome-wide genotyping array. Here, we analyse 1x sequence data across the APOC3 gene in a founder population from the island of Crete in Greece (n = 1239) and find significant evidence for association with blood triglyceride levels with the previously reported R19X cardioprotective null mutation (β = -1.09,σ = 0.163, P = 8.2 × 10^-11) and a second loss of function mutation, rs138326449 (β = -1.17,σ = 0.188, P = 1.14 × 10^-9). The signal cannot be recapitulated by imputing genome-wide genotype data on a large reference panel of 5122 individuals including 249 with 4x WGS data from the same population. Gene-level meta-analysis with other studies reporting burden signals at APOC3 provides robust evidence for a replicable cardioprotective rare variant aggregation (P = 3.2 × 10^-31, n = 13 480).

Consistent Estimation in Mendelian Randomization with Some Invalid Instruments Using a Weighted Median Estimator

Developments in genome-wide association studies and the increasing availability of summary genetic association data have made application of Mendelian randomization relatively straightforward. However, obtaining reliable results from a Mendelian randomization investigation remains problematic, as the conventional inverse-variance weighted method only gives consistent estimates if all of the genetic variants in the analysis are valid instrumental variables. We present a novel weighted median estimator for combining data on multiple genetic variants into a single causal estimate. This estimator is consistent even when up to 50% of the information comes from invalid instrumental variables. In a simulation analysis, it is shown to have better finite-sample Type 1 error rates than the inverse-variance weighted method, and is complementary to the recently proposed MR-Egger (Mendelian randomization-Egger) regression method. In analyses of the causal effects of low-density lipoprotein cholesterol and high-density lipoprotein cholesterol on coronary artery disease risk, the inverse-variance weighted method suggests a causal effect of both lipid fractions, whereas the weighted median and MR-Egger regression methods suggest a null effect of high-density lipoprotein cholesterol that corresponds with the experimental evidence. Both median-based and MR-Egger regression methods should be considered as sensitivity analyses for Mendelian randomization investigations with multiple genetic variants.