Phenome-wide Mendelian randomization study of plasma triglycerides and 2,600 disease traits

Joshua K. Park, Shantanu Bafna, Iain S. Forrest, Áine Duffy, Carla Marquez-Luna, Ben O. Petrazzini, Ha My Vy, Daniel M. Jordan, Marie Verbanck, Jagat Narula, Robert S. Rosenson, Ghislain Rocheleau, Ron Do
doi: https://doi.org/10.1101/2022.07.21.22277900

Abstract

Background Causality between plasma triglyceride (TG) levels and atherosclerotic cardiovascular disease (ASCVD) risk remains controversial despite more than four decades of study and two recent landmark trials, STRENGTH and REDUCE-IT. Further unclear is the association between TG levels and non-atherosclerotic diseases across organ systems.

Methods Here, we conducted a phenome-wide, two-sample Mendelian randomization (MR) analysis using inverse-variance weighted (IVW) regression to systematically infer the causal effects of plasma TG levels on 2,600 disease traits in the European ancestry population of UK Biobank. For replication, we externally tested 221 nominally significant associations (p < 0.05) in an independent cohort from FinnGen. To account for potential horizontal pleiotropy and the influence of invalid instrumental variables, we performed sensitivity analyses using MR-Egger regression, weighted median estimator, and MR-PRESSO. Finally, we used multivariable MR controlling for correlated lipid fractions to distinguish the independent effect of plasma TG levels.

Results Our results identified 7 disease traits reaching Bonferroni-corrected significance in both the discovery (p < 1.92 × 10-5) and replication analyses (p < 2.26 × 10-4), supporting a causal relationship between plasma TG levels and ASCVDs, including coronary artery disease (OR 1.33, 95% CI 1.24-1.43, p = 2.47 × 10-13). We also identified 12 disease traits that were Bonferroni-significant in the discovery or replication analysis and at least nominally significant in the other analysis (p < 0.05), identifying plasma TG levels as a novel risk factor for 9 non-ASCVD diseases, including uterine leiomyoma (OR 1.19, 95% CI 1.10-1.29, p = 1.17 × 10-5).

Conclusions Taking a phenome-wide, two-sample MR approach, we identified causal associations between plasma TG levels and 19 disease traits across organ systems. Our findings suggest unrealized drug repurposing opportunities or adverse effects related to approved and emerging TG-lowering agents as well as mechanistic insights for future study.

Introduction

Atherosclerotic cardiovascular disease (ASCVD) remains the leading cause of death worldwide despite the effectiveness of statin therapy in reducing low-density lipoprotein cholesterol (LDL-C) levels (Baigent et al., 2005) (Roth et al., 2020). Additional therapeutic targets and adjunctive treatments are needed to address the burden arising from residual risk (Rosenson & Goonewardena, 2021). Triglycerides (TG) play vital roles in physiology, ranging from energy storage and mobilization to inflammation, thrombosis, and hormone-like signaling (Zewinger et al., 2020) (Norata et al., 2007). However, a causal relationship between TGs and ASCVDs remains controversial (Miller et al., 2011) (Albrink & Man, 1959), recently culminating in the conflicting reports of two double-blinded randomized controlled trials (RCT), STRENGTH and REDUCE-IT (Nicholls et al., 2020) (Bhatt et al., 2019) (Doi et al., 2021). Nevertheless, drug development for reducing TG-rich lipoproteins (TRL) is an active area of research and several targets have now been validated, including angiopoietin-like 3 (ANGPTL3) (Graham et al., 2017) (Dewey et al., 2017) (Musunuru et al., 2010) (Gaudet et al., 2017), angiopoietin-like 4 (ANGPTL4) (Dewey et al., 2016), and apolipoprotein C-III (APOC3) (Gaudet et al., 2015) (Crosby et al., 2014). Clinical trials are currently evaluating these targets for dyslipidemias (Arrowhead, 2021b) (Arrowhead, 2021a).

Whether TGs are causal risk factors or simply associative biomarkers remains uncertain not only for ASCVDs but also other diseases of different organ systems. Understanding the causal effects of TGs across a broader range of human diseases could have significant implications for drug repurposing. TG-lowering agents, such as fibrates (Frick et al., 1987) (Group, 2010) and omega-3 fatty acids (Bhatt et al., 2019) have already been approved; however, not knowing which diseases are causally affected by TGs precludes their use for indications other than hypertriglycidemia. Understanding the protective effects of TGs could also have implications for drug safety. With the recent approval of icosapent ethyl (Bhatt et al., 2019) (Gaba et al., 2022) and the ongoing development of other TG-lowering agents (Shaik & Rosenson, 2021), such as ANGPTL3 (Graham et al., 2017) (Dewey et al., 2017), ANGPTL4, and APOC3 inhibitors (Gaudet et al., 2015), long-term safety becomes a matter of concern. Post-market surveillance data are limited; therefore, identifying diseases with negative causal links to TGs could suggest adverse side effects, whereas protective causal links to TGs could suggest therapeutic avenues and indices informative for drug development. Further, this could inform polypharmacy and drug titration in clinical practice.

Several methodological challenges have prohibited causal conclusions about TGs across human diseases. First, TGs are correlated with established ASCVD risk factors, such as obesity and insulin resistance (Eckel et al., 2005). They also correlate with LDL particles or apolipoprotein B (apoB) concentration (Cromwell et al., 2007) and inversely correlate with high-density lipoprotein cholesterol (HDL-C) (Phillips & Smith, 1991). Conventional observational studies have thus been limited in drawing causal inferences due to potential confounding. Secondly, available TG-lowering agents have pleiotropic effects on several major lipid fractions, including VLDL-C, LDL-C, HDL-C, and apolipoproteins, including APOC3 (Group, 2010) (Keech et al., 2005) (Frick et al., 1987). Thus, costly RCTs have had limited power to disentangle the effects of lowering TG specifically (Triglyceride Coronary Disease Genetics et al., 2010) (Goldberg et al., 2011) (Bhatt et al., 2019). Finally, Mendelian randomization (MR) studies have typically focused on individual diseases selected a priori (Davey Smith & Hemani, 2014), narrowing the scope of causal estimates to ASCVDs while overlooking non-ASCVD diseases (Harrison et al., 2018) (Smith et al., 2014) (Allara et al., 2019).

Phenome-wide MR is a high-throughput extension of MR that, under specific assumptions, estimates the causal effects of an exposure on multiple outcomes simultaneously. As in conventional MR, this method uses genetic variants as instrumental variables (IV) to proxy modifiable exposures while minimizing confounding (Davey Smith & Ebrahim, 2003). A distinction, however, is that phenome-wide MR enables comprehensive scans of the phenotypic spectrum, limiting bias from prior assumptions and facilitating the discovery of unforeseen causal relationships. This has recently become feasible with the maturation of large-scale biobanks providing extensive genetic and phenotypic data, such as the UK Biobank (UKB) (Bycroft et al., 2018) and FinnGen project (FinnGen, 2020).

Here, we use phenome-wide MR to systematically estimate the causal effects of plasma TGs on 2,600 disease outcomes in UKB, followed by replication testing in FinnGen. We then apply multiple MR methods for sensitivity analyses and multivariable MR to control for plasma LDL-C and HDL-C level.

Methods

Study design and data sources

We followed the guidelines published by Burgess et al. (Burgess et al., 2020) and the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) guidelines (Skrivankova et al., 2021). Accordingly, we note that MR analyses rely on three important instrumental variable assumptions: (1) the genetic variant is directly associated with the exposure; (2) the genetic variant is unrelated to confounders between the exposure and outcome; and (3) the genetic variant has no effect on the outcome other than through the exposure (Davey Smith & Ebrahim, 2003). This study uses three non-overlapping genome-wide association studies (GWAS). A schematic figure summarizes the study design (Fig 1).

Fig 1.
Fig 1. Study overview.

A schematic summarizing the study design. GLGC, Global Lipids Genetics Consortium; HDL, high-density lipoprotein; IVW, inverse-variance weighted; LDL, low-density lipoprotein; MR, Mendelian randomization; MR-PRESSO, Mendelian Randomization Pleiotropy RESidual Sum and Outlier; SNP, single-nucleotide polymorphism; TG, triglyceride; UKB, UK Biobank. Tier 1 (BB): At least Bonferroni-significant in both the discovery and replication analyses. Tier 2 (BN): At least Bonferroni-significant in the discovery analysis and at least nominally significant in the replication analysis. Tier 3 (NB): At least nominally significant in the discovery analysis and at least Bonferroni-significant in the replication analysis.

For primary and sensitivity analyses, we conducted two-sample MR taking summary statistics for genetic associations with plasma TGs from one dataset (Global Lipids Genetics Consortium, GLGC) (Willer et al., 2013) and summary statistics for genetic associations with outcomes from a second, independent dataset (UKB). The CARDIoGRAMplusC4D Metabochip study by GLGC involved 63,746 cases and 130,681 controls (Deloukas et al., 2013). UKB is a longitudinal, population-based cohort study with genetic and phenotypic data on over 500,000 participants aged 40-69 years at recruitment from across the United Kingdom during 2006-2010 (Sudlow et al., 2015). The sociodemographic and health-related characteristics of UKB participants have been described elsewhere (Fry et al., 2017).

For replication analysis, we again performed two-sample MR but used summary statistics for TG from UKB and summary statistics for outcomes from FinnGen (Release 4; 176,899 samples; 169,962,023 variants). FinnGen is a large public-private partnership started in August 2017 aiming to collect and analyze genomic and phenotypic data from 500,000 Finnish biobank participants (FinnGen, 2020) (Kurki et al., 2022). For both primary and replication analysis, there is likely to be minimal sample overlap in the two-sample MR, which can cause weak-instrument bias and inflated type 1 error (Burgess et al., 2016).

Genetic instruments

For primary and secondary analyses, we identified 3,086 single nucleotide polymorphisms (SNP) associated with plasma TG levels in GLGC at a genome-wide significance threshold of p < 5 × 10-8. A total of 141 independent SNPs were then selected at a linkage disequilibrium (LD) threshold of r2 < 0.05 using the 1000 Genomes LD European panel as the reference population (Auton et al., 2015). To select independent SNPs, we used the PLINK software’s LD clumping command (Purcell et al., 2007) with the following options: -- clump-p1 0.00000005 --clump-p2 0.0005 --clump-r2 0.05 --clump-kb 250.

For replication analysis, we identified 1,388 SNPs associated with plasma TG levels in UKB after restricting for genome-wide significance (p < 5 × 10-8) and LD-clumping. The SNPs were selected at an LD threshold of r2 < 0.05 using the 1000 Genomes LD European panel as the reference population. LD clumping was performed by PLINK using the same command as above. We then restricted to 1,248 SNPs for genetic instrumentation based on whether the SNP was also available in the FinnGen dataset.

Disease outcomes

As a phenome-wide investigation, this study examined many binary disease outcomes. For outcomes of primary and sensitivity analyses, we used genome-wide association summary statistics from the European ancestry subset of UKB. The focus on European ancestry was needed to reduce heterogeneity and maximize statistical power. Pan-UKB performed 16,000 genome-wide association studies in ∼500,000 UKB participants of various ancestries using genetic and phenotypic data (PanUKBTeam, 2020). A total of 7,221 total phenotypes were categorized as “biomarker”, “continuous”, “categorical”, “ICD-10 code”, “phecode”, or “prescription” (PanUKBTeam, 2020). We filtered for outcomes to retain categorical, ICD-10, and phecode types; non-null heritability in European ancestry as estimated by Pan-UKB; and relevance to disease, excluding medications. This yielded 2,600 traits for primary analysis.

For outcomes of replication analysis, we used genome-wide association summary statistics from FinnGen. To allow for compatibility between UKB and FinnGen datasets, phenotypes coded as FinnGen “endpoint IDs” were mapped to ICD-10 codes or phecodes corresponding to the coding convention used by UKB. Categorical traits in UKB could not be reliably mapped to FinnGen IDs, requiring us to omit categorical traits from replication analysis. Of 2,444 available FinnGen outcomes, we ultimately selected 221 outcomes for replication testing based on two conditions: (1) The outcome has an equivalent phenotype documented as an ICD-10 code or phecode in UKB, which we had included in the discovery analysis, and (2) the outcome was at least nominally significant (p < 0.05) in the discovery analysis.

Statistical analyses

Primary discovery analysis used the inverse-variance weighted (IVW) method of two-sample MR over the more horizontal pleiotropy-robust MR-Egger regression method to maximize statistical power (Burgess et al., 2013). To proxy TG, we selected the 141 TG-associated SNPs from Pan-UKB identified above as instrumental variables in MR testing. All SNPs were harmonized to match the effect and non-effect allele at each SNP. For replication analyses, we conducted IVW tests using TG-associated SNPs from UKB as exposures and 221 traits at least nominally significant (p < 0.05) in the discovery analysis as outcomes from FinnGen. R (4.1.0) was used for all analyses (R Core Team, 2021).

For sensitivity analyses, we addressed the possible presence of horizontal pleiotropy by applying MR-Egger (Bowden et al., 2015), weighted median estimator (Bowden et al., 2016), and MR-PRESSO outlier tests (Verbanck et al., 2018) on Bonferroni-significant traits identified in the primary analysis above. MR-PRESSO outlier tests were performed on 16 traits with significant MR-PRESSO global test results (p < 0.05). To account for potential horizontal pleiotropy due to outliers, we examined the MR-PRESSO outlier test p-values from the discovery analysis, and as recommended (Burgess et al., 2020), removed outlier IVs (between 1 and 10 SNPs for each outcome) based on a corrected significance threshold (p < 0.05 / 141 = 3.55 × 10-4). IVW tests were then re-run without the outlier IVs in the discovery analysis.

Additionally, we fitted multivariable MR (MVMR) models controlling for plasma LDL-C and HDL-C levels in the discovery stage. MVMR uses genetic variants associated with multiple exposures to estimate the effect of each exposure on an outcome, thereby accounting for potentially correlated exposures (Sanderson, 2021). We extracted the 141 LD-independent TG-associated SNPs described above from the summary statistics for LDL and HDL provided by GLGC. We then performed three MVMR IVW tests against traits Bonferroni-significant in our discovery analysis with the following sets of exposures: TG and LDL; TG and HDL; or TG, HDL, and LDL.

To account for multiple testing, a conservative Bonferroni threshold for statistical significance was set to p < 0.05 / 2600 = 1.92 × 10-5 for the discovery analysis. For replication analysis, a similar conservative Bonferroni threshold was set to p < 0.05 / 221 = 2.26 × 10-4. We determined three tiers of statistical evidence based on the Bonferroni or nominal threshold (p < 0.05) in both the discovery and replication analyses. These three tiers are defined below:

  1. Tier 1 (BB): At least Bonferroni-significant in both the discovery and replication analyses.

  2. Tier 2 (BN): At least Bonferroni-significant in the discovery analysis and at least nominally significant in the replication analysis.

  3. Tier 3 (NB): At least nominally significant in the discovery analysis and at least Bonferroni-significant in the replication analysis.

Ethical approval

UK Biobank has approval from the North West Multi Centre Research Ethics Committee (MREC) as a Research Tissue Bank (RTB) (11/NW/0382), and all participants of UKB provided written informed consent. More information is available at (https://www.ukbiobank.ac.uk/learn-more-about-uk-biobank/about-us/ethics). The work described in this study was approved by UKB under application number 16218. All participants of FinnGen provided written informed consent for biobank research, based on the Finnish Biobank Act. The Coordinating Ethics Committee of the Hospital District of Helsinki and Uusimaa (HUS) approved the FinnGen study protocol Nr HUS/990/2017. More information is available at (https://www.finngen.fi/en/code_of_conduct).

Data availability

All data generated in this study are included in the manuscript and supplementary tables.

All analyses used publicly available data (UKB, FinnGen), including previously published GWAS (GLGC) (Willer et al., 2013). Obtaining access to UKB (PanUKBTeam, 2020) and FinnGen (FinnGen, 2020) GWAS summary statistics is detailed here (https://www.finngen.fi/en/access_results) and here (https://www.ukbiobank.ac.uk/enable-your-research/apply-for-access). Please note the summary statistics for FinnGen and Pan-UKB are made publicly available.

Results

Genetically proxied plasma TG levels and phenome-wide disease risk

In the discovery analysis, we identified nominally significant associations (p < 0.05) between plasma TG levels and 598 disease traits in UKB (S1 Table). Of these, 39 disease traits were statistically significant after multiple testing correction with a conservative Bonferroni-corrected threshold (p < 1.92 × 10-5). As a positive control, plasma TG levels were positively associated with gout with an odds ratio (OR) of 1.78 for gout (95% CI 1.52-2.09, p = 7.41 × 10- 11), in agreement with prior studies (Yu et al., 2021). In the replication analysis, we identified nominally significant associations (p < 0.05) between plasma TG levels and 71 disease traits in FinnGen (S2 Table). Of these, 22 traits were Bonferroni-significant. A summary of the 19 most statistically significant and replicated results have been organized into three predefined tiers of evidence (Fig 2). We note that the magnitude of estimates from the discovery analysis were generally greater than those from replication analysis, especially among tier 1 and tier 2 associations, potentially as a manifestation of winner’s curse bias (Göring et al., 2001).

Fig 2.
Fig 2. Causal estimates of genetically proxied plasma TG levels on disease risk using IVW regression in UKB and FinnGen.

Causal estimates from inverse-variance weighted (IVW) regression are shown as odds ratios (OR) per 1 SD increase in plasma triglyceride (TG) levels (mmol/L). Asterisks indicate novel associations. Black indicates discovery analysis results using UKB. Red indicates replication analysis results using FinnGen. Horizontal lines represent 95% confidence intervals. Tier 1 (BB): At least Bonferroni-significant in both the discovery and replication analyses. Tier 2 (BN): At least Bonferroni-significant in the discovery analysis and at least nominally significant in the replication analysis. Tier 3 (NB): At least nominally significant in the discovery analysis and at least Bonferroni-significant in the replication analysis. Full results are provided in S1 and S2 Tables.

For tier 1 results, genetically determined plasma TG levels were positively associated with 7 disease traits in both the discovery and replication analyses. These were Bonferroni-significant in both analyses, and all were related to dyslipidemias or ASCVD. Among traits related to ASCVD, the strongest association by statistical significance was for angina pectoris in both the discovery UKB cohort (OR 1.39, 95% CI 1.29-1.51, p = 2.11 × 10-13) and the replication FinnGen cohort (OR 1.30, 95% CI 1.24-1.37, p = 1.25 × 10-26). For tier 2 results, plasma TG levels were positively associated with 3 disease traits, including non-ASCVDs: gout, uterine leiomyoma, and “other aneurysms” (phecode-442). The strongest association by significance, after gout, was for leiomyoma of the uterus in both the discovery (OR 1.19, 95% CI 1.10-1.29, p = 1.17 × 10-5) and replication cohorts (OR 1.06, 95% CI 1.02-1.11, p = 3.60 × 10-3). For tier 3 results, 9 disease traits were identified, and a greater proportion were non-ASCVD traits. The strongest association by significance was for hypertension in both the discovery (OR 1.51, 95% CI 1.08-1.23, p = 3.20 × 10-5) and replication cohorts (OR 1.17, 95% CI 1.14-1.22, p = 9.49 × 10-20).

Sensitivity analyses using MR-PRESSO, MR-Egger, and Weighted Median methods

To account for potential horizontal pleiotropy bias due to outlier IVs, we next used the MR-PRESSO test on the 19 significant and replicated associations identified above, which were categorized into three tiers of statistical evidence. MR-PRESSO suspected outlier IVs for 16 of the 19 associations, and between 1 and 10 IVs were identified for 8 associations (S3 Table). We then re-ran IVW MR for these 8 associations after removing outlier IVs (S3 Table). Among tier 1 results, all 7 associations increased in significance after outlier removal (Fig 3). For tier 2 results, two of three associations had significant global MR-PRESSO test values but no outlier IVs were removed. For tier 3 results, seven of nine associations had significant global MR-PRESSO test values, and outlier IVs were removed for one of these associations (S3 Table), which increased in significance: hypertension (p = 3.20 × 10-5 to p = 2.80 × 10-8). Importantly, all associations maintained the same effect direction after the MR-PRESSO outlier test, in keeping with the initial discovery analysis.

Fig 3.
Fig 3. Causal estimates of genetically proxied plasma TG levels on disease risk using MR-Egger, Weighted Median, and MR-PRESSO methods in UKB.

Shown are sensitivity analysis results from the discovery stage using instruments from GLGC and outcomes from UKB. Levels of statistical significance are expressed as -log10(p) values. Associations with insignificant global MR-PRESSO test results (p > 0.05) were not rerun and do not have data points for IVW after MR-PRESSO. Tier 1 (BB): At least Bonferroni-significant in both the discovery and replication analyses. Tier 2 (BN): At least Bonferroni-significant in the discovery analysis and at least nominally significant in the replication analysis. Tier 3 (NB): At least nominally significant in the discovery analysis and at least Bonferroni-significant in the replication analysis. Full results are found in S1 and S3 Tables.

To further account for horizontal pleiotropy, we conducted sensitivity analyses for the 19 significant and replicated traits using the MR-Egger and weighted median estimators on UKB data (Fig 3). MR-Egger results were comparable to IVW results for all 19 traits across the three tiers of evidence. However, weighted median results were less statistically significant than IVW and MR-Egger results for most traits except for disorders of lipid metabolism, hyperlipidemia, and hypertension.

Separating the independent effects of plasma TG levels

To isolate the independent effect of plasma TG levels from those of correlated lipid fractions, we next conducted multivariable MR (MVMR) on the 19 significant and replicated associations using the IVW estimator on UKB data, adjusting for plasma HDL-C, LDL-C, or both simultaneously (Fig 4). For tier 1 associations, we observed an increase in statistical significance for disorders of lipid metabolism and hyperlipidemia after adjustment. For tier 2 results, we observed an increase in significance for gout. For tier 3 results, we observed an increase in significance for alcoholic liver disease, paroxysmal tachycardia, and calculus of kidney and ureter. Aside from these, MVMR led to a decrease in statistical significance for the remaining associations. However, all 19 associations remained in the same direction of causality despite multivariable adjustment (S4 Table).

Fig 4.
Fig 4. Causal estimates of genetically proxied plasma TG levels on disease risk using multivariable IVW regression controlling for plasma LDL-C and HDL-C levels in UKB.

Shown are the multivariable MR results from the discovery analysis, using genetic instruments for plasma TG levels from GLGC and disease traits from UKB. Degrees of statistical significance are categorized into three tiers: Tier 1 (BB): At least Bonferroni-significant in both the discovery and replication analyses. Tier 2 (BN): At least Bonferroni-significant in the discovery analysis and at least nominally significant in the replication analysis. Tier 3 (NB): At least nominally significant in the discovery analysis and at least Bonferroni-significant in the replication analysis. Full results are provided in S4 Table.

Discussion

We performed phenome-wide, two-sample MR to estimate the causal effects of plasma TG levels on a spectrum of human disease traits (n = 2,600) using a discovery cohort from UKB and a replication cohort from FinnGen. We report 7 disease traits reaching Bonferroni-corrected significance in both the discovery and replication analyses, which we categorized as tier 1 results. These traits were predominantly manifestations of ASCVD, such as ischemic heart disease, angina pectoris, and myocardial infarction. Using multivariable MR, we found that these associations remained even after controlling for LDL-C and HDL-C (Fig 4). We also identified 3 disease traits that are both Bonferroni-significant in the discovery analysis (p < 1.92 × 10-5) and at least nominally significant in the replication analysis (p < 0.05), categorized as tier 2 results. Lastly, we identified 9 disease traits at least nominally significant in the discovery analysis and Bonferroni-significant in replication analysis, categorized as tier 3 results. Several of these disease traits have never been reported to be causally associated with TGs, introducing new opportunities for drug repurposing and further mechanistic studies.

Our results are consistent with prior work suggesting that plasma TG levels are causally associated with ASCVD risk (Triglyceride Coronary Disease Genetics et al., 2010) (Castañer et al., 2020) (Do et al., 2013) (Dewey et al., 2016) (Holmes et al., 2015) (White et al., 2016) (Varbo et al., 2013) (Ibi et al., 2021) (Rosenson et al., 2021). We do not prove here that circulating TGs per se are directly atherogenic; rather, we show that plasma TG measurement, which comprises multiple classes of triglyceride-rich lipoproteins (TRL) and their remnants, captures a clinically significant mechanism that is causally associated with ASCVD (Rosenson et al., 2021). This interpretation is concordant with the emerging consensus that apolipoprotein B (apoB) particle number is a more important determinant of atherogenesis than LDL-C and that TRL particles are as important a risk factor in ASCVD as LDL particles, since they both carry a single apoB molecule (Ference et al., 2019) (Marston et al., 2021) (Richardson et al., 2020). According to this paradigm, the causal association we observe between plasma TG levels and ASCVD risk may be mediated by the concentration of apoB particles in TRLs, captured by plasma TG level measurements. If reducing TRLs proportionally reduces apoB levels (Ference et al., 2019), these results support that lowering plasma TG levels could serve as a valid therapeutic strategy to reducing ASCVD risk. Additional studies controlling for apoB are needed to test whether decreasing plasma TGs without a concomitant decrease in apoB could reduce ASCVD risk by alternative mechanisms unrelated to apoB.

Regarding non-ASCVDs, we show novel causal associations between plasma TG levels and uterine leiomyomas (uterine fibroids), diverticular disease of intestine, paroxysmal tachycardia, hemorrhage from respiratory passages (hemoptysis), and calculus of kidney and ureter (kidney stones). Prior studies had reported correlational associations between many of these diseases and plasma TG levels or risk factors correlated to plasma TGs. For example, studies had documented positive correlations between leiomyoma risk and plasma TG levels (Uimari et al., 2016) (Tonoyan et al., 2021) (Peshkova et al., 2020). Others had shown diverticular disease risk is positively associated with BMI, body fat percentage, and visceral fat area (Shih et al., 2022) (Freckelton et al., 2018) (Böhm, 2021). Similarly, kidney stones have been associated with the triglyceride-glucose (TyG) index (Qin et al., 2021), and atrial tachycardias have been associated with VLDL levels in metabolic syndrome patients (Lee et al., 2017) (Park & Lee, 2018). Our study is the first to suggest that these associations may be causally related and specific to plasma TG levels, as opposed to or in addition to confounding risk factors, such as obesity and hormone replacement therapy. We thus present rationale for repurposing TG-lowering agents towards these diseases and pursuing mechanistic studies interrogating TG biology; however, RCTs are necessary to better evaluate clinical potential.

Our results also identified a novel, negative causal association between plasma TG levels and alcoholic liver disease (ALD), suggesting that excessively reducing plasma TG levels may increase risk of this disease trait. A mechanistic explanation for this association remains elusive as previous studies on dysregulated lipid metabolism in liver disease have generally focused on non-alcoholic fatty liver disease (NAFLD). However, one animal study suggests that medium-chain TGs may decrease lipid peroxidation and reverse established alcoholic liver injury in rat models of ALD (Nanji et al., 1996). Nevertheless, that elevated TG levels may protect against ALD is surprising and requires validation. It remains unclear whether this association is only relevant to lifelong, chronic lowering of plasma TG levels rather than transient lowering by drug-based interventions. It is also unclear whether ALD could be prevented by increasing TG levels. Our finding emphasizes the potential intolerability of TG-lowering agents and the importance of maintaining these drugs’ concentrations within therapeutic windows for patient safety.

This study has several limitations. First, horizontal pleiotropy has been shown to be common in genetic variation, which can cause bias in MR testing (Verbanck et al., 2018) (Jordan et al., 2019). We used MR tests that are robust to horizontal pleiotropy including the MR-PRESSO test as well as MR-Egger and weighted median estimators as sensitivity analyses. Second, genetic variants confer exposures that are lifelong but small in effect size; thus, MR may over- or underestimate the effect sizes of pharmacological interventions. MR also cannot make comparisons between TG-lowering agents as it evaluates drug targets, not drug subclasses and unique pharmacodynamics (Gill et al., 2019) (Sofat et al., 2010). However, the utility of MR in drug target validation is well-established as estimates still indicate the presence and direction of causality (Gill et al., 2021). Third, we were unable to reliably map “categorical” traits from UKB to corresponding traits in FinnGen for all traits, testing for replication only 221 of the 598 associations that were nominally significant in the primary analysis. Moreover, ICD-10 codes and phecodes are imperfect descriptors of disease, liable to misclassification. Lastly, we acknowledge that this study was restricted to populations of European ancestry, limiting the generalizability of our findings. A recent study observed comparable MR estimates for the causal effects of lipid traits on ischemic stroke risk between African and European ancestry individuals (Fatumo et al., 2021), suggesting the potential generalizability of MR results across ancestries in some specific cases. Nevertheless, trans-ancestry MR analyses are warranted to validate our study’s findings in diverse ancestry populations.

In conclusion, this study demonstrates a high-throughput application of two-sample MR to estimate the causal effects of plasma TGs on disease risk phenome-wide. With the proliferation of multiomic data, this systematic approach could be generalized to study the causal effects of serum biomarkers at scale to support the prioritization of targets for drug discovery. Further studies are needed however to consider the functional and qualitative attributes of TRL and TG subtypes. Metabolomic data partitioning TRLs by size and composition may enable this to identify the subcomponent of a plasma TG measurement that drives disease risk (Holmes et al., 2015). Diverse study designs are warranted to further triangulate evidence for the causal effects of plasma TG levels on disease risk.

Data Availability

All analyses used publicly available data (UKB, FinnGen), including previously published GWAS (GLGC) (Willer et al., 2013). Obtaining access to UKB (Pan-UKB_Team, 2020) and FinnGen (FinnGen, 2020) GWAS summary statistics is detailed here (https://www.finngen.fi/en/access_results) and here (https://pan.ukbb.broadinstitute.org/downloads). Please note the summary statistics for FinnGen and Pan-UKB are made publicly available.

Competing interests

RD reports receiving grants from AstraZeneca; grants and non-financial support from Goldfinch Bio; being a scientific co-founder, consultant, and equity holder (pending) for Pensieve Health; and a consultant for Variant Bio, all unrelated to this work. RSR reports receiving grants from Amgen, Arrowhead, Lilly, Novartis and Regeneron; consulting fees from Amgen, Arrowhead, Lilly, Novartis and Regeneron; honoraria for non-promotional lectures from Amgen, Kowa and Regeneron, royalties from Wolters Kluwer (UpToDate); and stock holdings in MediMergent, LLC.

Supplementary files

S1 Table. Causal estimates of plasma TG levels on 2,600 traits in UKB using multiple MR methods.

Shown are the estimates, standard deviations, and p-values of MR results using IVW, MR-Egger, and Weighted Median methods. 141 SNPs were used as instrumental variables to proxy plasma TG levels. A Bonferroni threshold for statistical significance was set to p < 0.05 / 2600 = 1.92 × 10-5 for this discovery analysis.

S2 Table. Causal estimates of plasma TG levels on 221 traits in FinnGen using multiple MR methods.

Shown are the estimates, standard deviations, and p-values of MR results using IVW, MR-Egger, and Weighted Median methods. 1,248 SNPs were used as instrumental variables to proxy plasma TG levels. A Bonferroni threshold for statistical significance was set to p < 0.05 / 221 = 2.26 × 10-4 for this replication analysis.

S3 Table. IVW-MR estimates of significant and replicated associations (Tier 1-3) after MR-PRESSO outlier tests using UKB data.

Shown are the estimates, standard deviations, and p-values of IVW-MR results in the discovery analysis, before and after outlier IV removal, for significant and replicated traits (tier 1-3) that had significant MR-PRESSO global test results (p < 0.05) in the primary IVW analysis of the discovery UKB cohort.

S4 Table. Multivariable IVW-MR estimates of plasma TG levels on significant and replicated associations (Tier 1-3) using UKB data.

Shown are the estimates, standard deviations, and p-values of multivariable IVW-MR results controlling for HDL-C, LDL-C, or both, in the discovery stage. Only tier 1-3 results significant and replicated as predefined in the methods were examined in this analysis.

Acknowledgements

We acknowledge the participants and investigators of the FinnGen and UKB studies.

Footnotes

  • Abbreviations: ALD, alcoholic liver disease; ANGPTL3, angiopoietin-like 3; ANGPTL4, angiopoietin-like 4; ApoB, apolipoprotein B; APOC3, apolipoprotein C-III; ASCVD, atherosclerotic cardiovascular disease; CAD, coronary artery disease; GLGC, Global Lipids Genetics Consortium; HDL-C, high-density lipoprotein cholesterol; IHD, ischemic heart disease; IV, instrumental variable; LDL-C, low-density lipoprotein cholesterol; MR-PRESSO, Mendelian randomization pleiotropy residual sum and outlier; MVMR, Multivariable Mendelian randomization; NAFLD, non-alcoholic fatty liver disease; OR, odds ratio; REDUCE-IT, Reduction of Cardiovascular Events with Icosapent Ethyl-Intervention Trial; STRENGTH, STatin Residual Risk Reduction With EpaNova in HiGh CV Risk PatienTs With Hypertriglyceridemia; TG, triglyceride

  • Funding: RD is supported by the National Institute of General Medical Sciences of the National Institutes of Health (NIH) (R35-GM124836) and the National Heart, Lung, and Blood Institute of the NIH (R01-HL139865 and R01-HL155915).

References

  1. 1.↵
    Albrink, M. J., & Man, E. B. (1959). Serum Triglycerides in Coronary Artery Disease. A.M.A. Archives of Internal Medicine, 103(1), 4–8. https://doi.org/10.1001/archinte.1959.00270010010002
    OpenUrl
  2. 2.↵
    Allara, E., Morani, G., Carter, P., Gkatzionis, A., Zuber, V., Foley, C. N., Rees, J. M. B., Mason, A. M., Bell, S., Gill, D., Lindström, S., Butterworth, A. S., Angelantonio, E. D., Peters, J., & Burgess, S. (2019). Genetic Determinants of Lipids and Cardiovascular Disease Outcomes. Circulation: Genomic and Precision Medicine, 12(12), e002711. https://doi.org/doi:10.1161/CIRCGEN.119.002711
    OpenUrl
  3. 3.
    Arrowhead. (2021a). Study of ARO-ANG3 in Adults With Mixed Dyslipidemia (ARCHES-2). U.S. National Library of Medicine. https://clinicaltrials.gov/ct2/show/NCT04832971
  4. 4.
    Arrowhead. (2021b). Study of ARO-APOC3 in Adults With Mixed Dyslipidemia (MUIR). U.S. National Library of Medicine. https://clinicaltrials.gov/ct2/show/NCT04998201
  5. 5.↵
    Auton, A., Abecasis, G. R., Altshuler, D. M., Durbin, R. M., Abecasis, G. R., Bentley, D. R., Chakravarti, A., Clark, A. G., Donnelly, P., Eichler, E. E., Flicek, P., Gabriel, S. B., Gibbs, R. A., Green, E. D., Hurles, M. E., Knoppers, B. M., Korbel, J. O., Lander, E. S., Lee, C., … Abecasis, G. R. (2015). A global reference for human genetic variation. Nature, 526(7571), 68–74. https://doi.org/10.1038/nature15393
    OpenUrl
  6. 6.↵
    Baigent, C., Keech, A., Kearney, P. M., Blackwell, L., Buck, G., Pollicino, C., Kirby, A., Sourjina, T., Peto, R., Collins, R., & Simes, R. (2005). Efficacy and safety of cholesterol-lowering treatment: prospective meta-analysis of data from 90,056 participants in 14 randomised trials of statins. Lancet (London, England), 366(9493), 1267–1278. https://doi.org/10.1016/s0140-6736(05)67394-1
    OpenUrl
  7. 7.↵
    Bhatt, D. L., Steg, P. G., Miller, M., Brinton, E. A., Jacobson, T. A., Ketchum, S. B., Doyle, R. T., Juliano, R. A., Jiao, L., Granowitz, C., Tardif, J.-C., & Ballantyne, C. M. (2019). Cardiovascular Risk Reduction with Icosapent Ethyl for Hypertriglyceridemia. New England Journal of Medicine, 380(1), 11–22. https://doi.org/10.1056/nejmoa1812792
    OpenUrlCrossRefPubMed
  8. 8.↵
    Böhm, S. K. (2021). Excessive Body Weight and Diverticular Disease. Visc Med, 37(5), 372–382. https://doi.org/10.1159/000518674
    OpenUrl
  9. 9.↵
    Bowden, J., Davey Smith, G., & Burgess, S. (2015). Mendelian randomization with invalid instruments: effect estimation and bias detection through Egger regression. International Journal of Epidemiology, 44(2), 512–525. https://doi.org/10.1093/ije/dyv080
    OpenUrlCrossRefPubMed
  10. 10.↵
    Bowden, J., Davey Smith, G., Haycock, P. C., & Burgess, S. (2016). Consistent Estimation in Mendelian Randomization with Some Invalid Instruments Using a Weighted Median Estimator. Genetic Epidemiology, 40(4), 304–314. https://doi.org/10.1002/gepi.21965
    OpenUrlCrossRefPubMed
  11. 11.↵
    Burgess, S., Butterworth, A., & Thompson, S. G. (2013). Mendelian Randomization Analysis With Multiple Genetic Variants Using Summarized Data. Genetic Epidemiology, 37(7), 658–665. https://doi.org/10.1002/gepi.21758
    OpenUrlCrossRefPubMed
  12. 12.↵
    Burgess, S., Davey Smith, G., Davies, N. M., Dudbridge, F., Gill, D., Glymour, M. M., Hartwig, F. P., Holmes, M. V., Minelli, C., Relton, C. L., & Theodoratou, E. (2020). Guidelines for performing Mendelian randomization investigations. Wellcome open research, 4, 186. https://doi.org/10.12688/wellcomeopenres.15555.2
    OpenUrlCrossRef
  13. 13.↵
    Burgess, S., Davies, N. M., & Thompson, S. G. (2016). Bias due to participant overlap in two-sample Mendelian randomization. Genet Epidemiol, 40(7), 597–608. https://doi.org/10.1002/gepi.21998
    OpenUrlCrossRefPubMed
  14. 14.↵
    Bycroft, C., Freeman, C., Petkova, D., Band, G., Elliott, L. T., Sharp, K., Motyer, A., Vukcevic, D., Delaneau, O., O’Connell, J., Cortes, A., Welsh, S., Young, A., Effingham, M., McVean, G., Leslie, S., Allen, N., Donnelly, P., & Marchini, J. (2018). The UK Biobank resource with deep phenotyping and genomic data. Nature, 562(7726), 203–209. https://doi.org/10.1038/s41586-018-0579-z
    OpenUrlCrossRefPubMed
  15. 15.↵
    Castañer, O., Pintó, X., Subirana, I., Amor, A. J., Ros, E., Hernáez, Á., Martínez-González, M., Corella, D., Salas-Salvadó, J., Estruch, R., Lapetra, J., Gómez-Gracia, E., Alonso-Gomez, A. M., Fiol, M., Serra-Majem, L., Corbella, E., Benaiges, D., Sorli, J. V., Ruiz-Canela, M., … Fitó, M. (2020). Remnant Cholesterol, Not LDL Cholesterol, Is Associated With Incident Cardiovascular Disease. J Am Coll Cardiol, 76(23), 2712–2724. https://doi.org/10.1016/j.jacc.2020.10.008
    OpenUrl
  16. 16.↵
    Cromwell, W. C., Otvos, J. D., Keyes, M. J., Pencina, M. J., Sullivan, L., Vasan, R. S., Wilson, P. W., & D’Agostino, R. B. (2007). LDL Particle Number and Risk of Future Cardiovascular Disease in the Framingham Offspring Study - Implications for LDL Management. J Clin Lipidol, 1(6), 583–592. https://doi.org/10.1016/j.jacl.2007.10.001
    OpenUrlCrossRefPubMed
  17. 17.↵
    Crosby, J., Peloso, G. M., Auer, P. L., Crosslin, D. R., Stitziel, N. O., Lange, L. A., Lu, Y., Tang, Z. Z., Zhang, H., Hindy, G., Masca, N., Stirrups, K., Kanoni, S., Do, R., Jun, G., Hu, Y., Kang, H. M., Xue, C., Goel, A., … Kathiresan, S. (2014). Loss-of-function mutations in APOC3, triglycerides, and coronary disease. N Engl J Med, 371(1), 22–31. https://doi.org/10.1056/NEJMoa1307095
    OpenUrlCrossRefPubMedWeb of Science
  18. 18.↵
    Davey Smith, G., & Ebrahim, S. (2003). ‘Mendelian randomization’: can genetic epidemiology contribute to understanding environmental determinants of disease?*. International Journal of Epidemiology, 32(1), 1–22. https://doi.org/10.1093/ije/dyg070
    OpenUrlCrossRefPubMedWeb of Science
  19. 19.↵
    Davey Smith, G., & Hemani, G. (2014). Mendelian randomization: genetic anchors for causal inference in epidemiological studies. Human molecular genetics, 23(R1), R89–R98. https://doi.org/10.1093/hmg/ddu328
    OpenUrlCrossRefPubMedWeb of Science
  20. 20.↵
    Deloukas, P., Kanoni, S., Willenborg, C., Farrall, M., Assimes, T. L., Thompson, J. R., Ingelsson, E., Saleheen, D., Erdmann, J., Goldstein, B. A., Stirrups, K., König, I. R., Cazier, J. B., Johansson, A., Hall, A. S., Lee, J. Y., Willer, C. J., Chambers, J. C., Esko, T., … Samani, N. J. (2013). Large-scale association analysis identifies new risk loci for coronary artery disease. Nat Genet, 45(1), 25–33. https://doi.org/10.1038/ng.2480
    OpenUrlCrossRefPubMed
  21. 21.↵
    Dewey, F. E., Gusarova, V., Dunbar, R. L., O’Dushlaine, C., Schurmann, C., Gottesman, O., McCarthy, S., Van Hout, C. V., Bruse, S., Dansky, H. M., Leader, J. B., Murray, M. F., Ritchie, M. D., Kirchner, H. L., Habegger, L., Lopez, A., Penn, J., Zhao, A., Shao, W., … Baras, A. (2017). Genetic and Pharmacologic Inactivation of ANGPTL3 and Cardiovascular Disease. N Engl J Med, 377(3), 211–221. https://doi.org/10.1056/NEJMoa1612790
    OpenUrlCrossRefPubMed
  22. 22.↵
    Dewey, F. E., Gusarova, V., O’Dushlaine, C., Gottesman, O., Trejos, J., Hunt, C., Van Hout, C. V., Habegger, L., Buckler, D., Lai, K.-M. V., Leader, J. B., Murray, M. F., Ritchie, M. D., Kirchner, H. L., Ledbetter, D. H., Penn, J., Lopez, A., Borecki, I. B., Overton, J. D., … Shuldiner, A. R. (2016). Inactivating Variants in ANGPTL4 and Risk of Coronary Artery Disease. New England Journal of Medicine, 374(12), 1123–1133. https://doi.org/10.1056/nejmoa1510926
    OpenUrlCrossRefPubMed
  23. 23.↵
    Do, R., Willer, C. J., Schmidt, E. M., Sengupta, S., Gao, C., Peloso, G. M., Gustafsson, S., Kanoni, S., Ganna, A., Chen, J., Buchkovich, M. L., Mora, S., Beckmann, J. S., Bragg-Gresham, J. L., Chang, H. Y., Demirkan, A., Den Hertog, H. M., Donnelly, L. A., Ehret, G. B., … Kathiresan, S. (2013). Common variants associated with plasma triglycerides and risk for coronary artery disease. Nat Genet, 45(11), 1345–1352. https://doi.org/10.1038/ng.2795
    OpenUrlCrossRefPubMed
  24. 24.↵
    Doi, T., Langsted, A., & Nordestgaard, B. G. (2021). A possible explanation for the contrasting results of REDUCE-IT vs. STRENGTH: cohort study mimicking trial designs. European Heart Journal, 42(47), 4807–4817. https://doi.org/10.1093/eurheartj/ehab555
    OpenUrl
  25. 25.↵
    Eckel, R. H., Grundy, S. M., & Zimmet, P. Z. (2005). The metabolic syndrome. The Lancet, 365(9468), 1415–1428. https://doi.org/https://doi.org/10.1016/S0140-6736(05)66378-7
    OpenUrl
  26. 26.↵
    Fatumo, S., Karhunen, V., Chikowore, T., Sounkou, T., Udosen, B., Ezenwa, C., Nakabuye, M., Soremekun, O., Daghlas, I., Ryan, D. K., Taylor, A., Mason, A. M., Damrauer, S. M., Vujkovic, M., Keene, K. L., Fornage, M., Järvelin, M.-R., Burgess, S., & Gill, D. (2021). Metabolic Traits and Stroke Risk in Individuals of African Ancestry: Mendelian Randomization Analysis. Stroke, 52(8), 2680–2684. https://doi.org/doi:10.1161/STROKEAHA.121.034747
    OpenUrl
  27. 27.↵
    Ference, B. A., Kastelein, J. J. P., Ray, K. K., Ginsberg, H. N., Chapman, M. J., Packard, C. J., Laufs, U., Oliver-Williams, C., Wood, A. M., Butterworth, A. S., Di Angelantonio, E., Danesh, J., Nicholls, S. J., Bhatt, D. L., Sabatine, M. S., & Catapano, A. L. (2019). Association of Triglyceride-Lowering LPL Variants and LDL-C–Lowering LDLR Variants With Risk of Coronary Heart Disease. JAMA, 321(4), 364–373. https://doi.org/10.1001/jama.2018.20045
    OpenUrlCrossRefPubMed
  28. 28.
    FinnGen. (2020). FinnGen Documentation of R4 Release https://finngen.gitbook.io/documentation/
  29. 29.↵
    Freckelton, J., Holt, D., Borsaru, A., Gwini, S., Croagh, D., & Moore, G. (2018). The role of body composition in diverticular disease. Int J Colorectal Dis, 33(9), 1299–1302. https://doi.org/10.1007/s00384-018-3058-y
    OpenUrl
  30. 30.↵
    Frick, M. H., Elo, O., Haapa, K., Heinonen, O. P., Heinsalmi, P., Helo, P., Huttunen, J. K., Kaitaniemi, P., Koskinen, P., Manninen, V., Mäenpää, H., Mälkönen, M., Mänttäri, M., Norola, S., Pasternack, A., Pikkarainen, J., Romo, M., Sjöblom, T., & Nikkilä, E. A. (1987). Helsinki Heart Study: Primary-Prevention Trial with Gemfibrozil in Middle-Aged Men with Dyslipidemia. New England Journal of Medicine, 317(20), 1237–1245. https://doi.org/10.1056/NEJM198711123172001
    OpenUrlCrossRefPubMedWeb of Science
  31. 31.↵
    Fry, A., Littlejohns, T. J., Sudlow, C., Doherty, N., Adamska, L., Sprosen, T., Collins, R., & Allen, N. E. (2017). Comparison of Sociodemographic and Health-Related Characteristics of UK Biobank Participants With Those of the General Population. American Journal of Epidemiology, 186(9), 1026–1034. https://doi.org/10.1093/aje/kwx246
    OpenUrlCrossRefPubMed
  32. 32.↵
    Gaba, P., Bhatt, D. L., Steg, P. G., Miller, M., Brinton, E. A., Jacobson, T. A., Ketchum, S. B., Juliano, R. A., Jiao, L., Doyle, R. T., Granowitz, C., Tardif, J.-C., Giugliano, R. P., Martens, F. M. A. C., Gibson, C. M., & Ballantyne, C. M. (2022). Prevention of Cardiovascular Events and Mortality With Icosapent Ethyl in Patients With Prior Myocardial Infarction. Journal of the American College of Cardiology, 79(17), 1660–1671. https://doi.org/https://doi.org/10.1016/j.jacc.2022.02.035
    OpenUrl
  33. 33.↵
    Gaudet, D., Alexander, V. J., Baker, B. F., Brisson, D., Tremblay, K., Singleton, W., Geary, R. S., Hughes, S. G., Viney, N. J., Graham, M. J., Crooke, R. M., Witztum, J. L., Brunzell, J. D., & Kastelein, J. J. (2015). Antisense Inhibition of Apolipoprotein C-III in Patients with Hypertriglyceridemia. N Engl J Med, 373(5), 438–447. https://doi.org/10.1056/NEJMoa1400283
    OpenUrlCrossRefPubMed
  34. 34.↵
    Gaudet, D., Gipe, D. A., Pordy, R., Ahmad, Z., Cuchel, M., Shah, P. K., Chyu, K. Y., Sasiela, W. J., Chan, K. C., Brisson, D., Khoury, E., Banerjee, P., Gusarova, V., Gromada, J., Stahl, N., Yancopoulos, G. D., & Hovingh, G. K. (2017). ANGPTL3 Inhibition in Homozygous Familial Hypercholesterolemia. N Engl J Med, 377(3), 296–297. https://doi.org/10.1056/NEJMc1705994
    OpenUrlCrossRefPubMed
  35. 35.↵
    Gill, D., Georgakis, M. K., Koskeridis, F., Jiang, L., Feng, Q., Wei, W.-Q., Theodoratou, E., Elliott, P., Denny, J. C., Malik, R., Evangelou, E., Dehghan, A., Dichgans, M., & Tzoulaki, I. (2019). Use of Genetic Variants Related to Antihypertensive Drugs to Inform on Efficacy and Side Effects. Circulation, 140(4), 270–279. https://doi.org/10.1161/circulationaha.118.038814
    OpenUrlPubMed
  36. 36.↵
    Gill, D., Georgakis, M. K., Walker, V. M., Schmidt, A. F., Gkatzionis, A., Freitag, D. F., Finan, C., Hingorani, A. D., Howson, J. M. M., Burgess, S., Swerdlow, D. I., Davey Smith, G., Holmes, M. V., Dichgans, M., Scott, R. A., Zheng, J., Psaty, B. M., & Davies, N. M. (2021). Mendelian randomization for studying the effects of perturbing drug targets. Wellcome open research, 6, 16–16. https://doi.org/10.12688/wellcomeopenres.16544.2
    OpenUrl
  37. 37.↵
    Goldberg, I. J., Eckel, R. H., & McPherson, R. (2011). Triglycerides and heart disease: still a hypothesis? Arteriosclerosis, thrombosis, and vascular biology, 31(8), 1716–1725. https://doi.org/10.1161/ATVBAHA.111.226100
    OpenUrlAbstract/FREE Full Text
  38. 38.↵
    Göring, H. H., Terwilliger, J. D., & Blangero, J. (2001). Large upward bias in estimation of locus-specific effects from genomewide scans. American journal of human genetics, 69(6), 1357–1369. https://doi.org/10.1086/324471
    OpenUrlCrossRefPubMedWeb of Science
  39. 39.↵
    Graham, M. J., Lee, R. G., Brandt, T. A., Tai, L.-J., Fu, W., Peralta, R., Yu, R., Hurh, E., Paz, E., McEvoy, B. W., Baker, B. F., Pham, N. C., Digenio, A., Hughes, S. G., Geary, R. S., Witztum, J. L., Crooke, R. M., & Tsimikas, S. (2017). Cardiovascular and Metabolic Effects of ANGPTL3 Antisense Oligonucleotides. New England Journal of Medicine, 377(3), 222–232. https://doi.org/10.1056/NEJMoa1701329
    OpenUrlCrossRefPubMed
  40. 40.↵
    Group, T. A. S. (2010). Effects of Combination Lipid Therapy in Type 2 Diabetes Mellitus. New England Journal of Medicine, 362(17), 1563–1574. https://doi.org/10.1056/NEJMoa1001282
    OpenUrlCrossRefPubMedWeb of Science
  41. 41.↵
    Harrison, S. C., Holmes, M. V., Burgess, S., Asselbergs, F. W., Jones, G. T., Baas, A. F., van ’t Hof, F. N., de Bakker, P. I. W., Blankensteijn, J. D., Powell, J. T., Saratzis, A., de Borst, G. J., Swerdlow, D. I., van der Graaf, Y., van Rij, A. M., Carey, D. J., Elmore, J. R., Tromp, G., Kuivaniemi, H., … Humphries, S. E. (2018). Genetic Association of Lipids and Lipid Drug Targets With Abdominal Aortic Aneurysm: A Meta-analysis. JAMA cardiology, 3(1), 26–33. https://doi.org/10.1001/jamacardio.2017.4293
    OpenUrl
  42. 42.↵
    Holmes, M. V., Asselbergs, F. W., Palmer, T. M., Drenos, F., Lanktree, M. B., Nelson, C. P., Dale, C. E., Padmanabhan, S., Finan, C., Swerdlow, D. I., Tragante, V., van Iperen, E. P., Sivapalaratnam, S., Shah, S., Elbers, C. C., Shah, T., Engmann, J., Giambartolomei, C., White, J., … Casas, J. P. (2015). Mendelian randomization of blood lipids for coronary heart disease. Eur Heart J, 36(9), 539–550. https://doi.org/10.1093/eurheartj/eht571
    OpenUrlCrossRefPubMed
  43. 43.↵
    Ibi, D., Blauw, L. L., Noordam, R., Dollé, M. E. T., Jukema, J. W., Rosendaal, F. R., Christodoulides, C., Neville, M. J., Koivula, R., Rensen, P. C. N., Karpe, F., & van Dijk, K. W. (2021). Triglyceride-lowering LPL alleles combined with LDL-C-lowering alleles are associated with an additively improved lipoprotein profile. Atherosclerosis, 328, 144–152. https://doi.org/10.1016/j.atherosclerosis.2021.04.015
    OpenUrl
  44. 44.↵
    Jordan, D. M., Verbanck, M., & Do, R. (2019). HOPS: a quantitative score reveals pervasive horizontal pleiotropy in human genetic variation is driven by extreme polygenicity of human traits and diseases. Genome Biol, 20(1), 222. https://doi.org/10.1186/s13059-019-1844-7
    OpenUrlCrossRefPubMed
  45. 45.↵
    Keech, A., Simes, R. J., Barter, P., Best, J., Scott, R., Taskinen, M. R., Forder, P., Pillai, A., Davis, T., Glasziou, P., Drury, P., Kesäniemi, Y. A., Sullivan, D., Hunt, D., Colman, P., d’Emden, M., Whiting, M., Ehnholm, C., & Laakso, M. (2005). Effects of long-term fenofibrate therapy on cardiovascular events in 9795 people with type 2 diabetes mellitus (the FIELD study): randomised controlled trial. Lancet (London, England), 366(9500), 1849–1861. https://doi.org/10.1016/s0140-6736(05)67667-2
    OpenUrl
  46. 46.↵
    Kurki, M. I., Karjalainen, J., Palta, P., Sipilä, T. P., Kristiansson, K., Donner, K., Reeve, M. P., Laivuori, H., Aavikko, M., Kaunisto, M. A., Loukola, A., Lahtela, E., Mattsson, H., Laiho, P., Della Briotta Parolo, P., Lehisto, A., Kanai, M., Mars, N., Rämö, J., … Palotie, A. (2022). FinnGen: Unique genetic insights from combining isolated population and national health register data. medRxiv, 2022.2003.2003.22271360. https://doi.org/10.1101/2022.03.03.22271360
  47. 47.↵
    Lee, H. C., Chen, C. C., Tsai, W. C., Lin, H. T., Shiao, Y. L., Sheu, S. H., Wu, B. N., Chen, C. H., & Lai, W. T. (2017). Very-Low-Density Lipoprotein of Metabolic Syndrome Modulates Gap Junctions and Slows Cardiac Conduction. Sci Rep, 7(1), 12050. https://doi.org/10.1038/s41598-017-11416-5
    OpenUrl
  48. 48.↵
    Marston, N. A., Giugliano, R. P., Melloni, G. E. M., Park, J.-G., Morrill, V., Blazing, M. A., Ference, B., Stein, E., Stroes, E. S., Braunwald, E., Ellinor, P. T., Lubitz, S. A., Ruff, C. T., & Sabatine, M. S. (2021). Association of Apolipoprotein B–Containing Lipoproteins and Risk of Myocardial Infarction in Individuals With and Without Atherosclerosis: Distinguishing Between Particle Concentration, Type, and Content. JAMA cardiology. https://doi.org/10.1001/jamacardio.2021.5083
  49. 49.↵
    Miller, M., Stone, N. J., Ballantyne, C., Bittner, V., Criqui, M. H., Ginsberg, H. N., Goldberg, A. C., Howard, W. J., Jacobson, M. S., Kris-Etherton, P. M., Lennie, T. A., Levi, M., Mazzone, T., & Pennathur, S. (2011). Triglycerides and Cardiovascular Disease. Circulation, 123(20), 2292–2333. https://doi.org/doi:10.1161/CIR.0b013e3182160726
    OpenUrlFREE Full Text
  50. 50.↵
    Musunuru, K., Pirruccello, J. P., Do, R., Peloso, G. M., Guiducci, C., Sougnez, C., Garimella, K. V., Fisher, S., Abreu, J., Barry, A. J., Fennell, T., Banks, E., Ambrogio, L., Cibulskis, K., Kernytsky, A., Gonzalez, E., Rudzicz, N., Engert, J. C., DePristo, M. A., … Kathiresan, S. (2010). Exome sequencing, ANGPTL3 mutations, and familial combined hypolipidemia. N Engl J Med, 363(23), 2220–2227. https://doi.org/10.1056/NEJMoa1002926
    OpenUrlCrossRefPubMedWeb of Science
  51. 51.↵
    Nanji, A. A., Yang, E. K., Fogt, F., Sadrzadeh, S. M., & Dannenberg, A. J. (1996). Medium chain triglycerides and vitamin E reduce the severity of established experimental alcoholic liver disease. Journal of Pharmacology and Experimental Therapeutics, 277(3), 1694–1700. https://jpet.aspetjournals.org/content/jpet/277/3/1694.full.pdf
    OpenUrlAbstract/FREE Full Text
  52. 52.↵
    Nicholls, S. J., Lincoff, A. M., Garcia, M., Bash, D., Ballantyne, C. M., Barter, P. J., Davidson, M. H., Kastelein, J. J. P., Koenig, W., McGuire, D. K., Mozaffarian, D., Ridker, P. M., Ray, K. K., Katona, B. G., Himmelmann, A., Loss, L. E., Rensfeldt, M., Lundström, T., Agrawal, R., … Nissen, S. E. (2020). Effect of High-Dose Omega-3 Fatty Acids vs Corn Oil on Major Adverse Cardiovascular Events in Patients at High Cardiovascular Risk. JAMA, 324(22), 2268. https://doi.org/10.1001/jama.2020.22258
    OpenUrlCrossRefPubMed
  53. 53.↵
    Norata, G. D., Grigore, L., Raselli, S., Redaelli, L., Hamsten, A., Maggi, F., Eriksson, P., & Catapano, A. L. (2007). Post-prandial endothelial dysfunction in hypertriglyceridemic subjects: molecular mechanisms and gene expression studies. Atherosclerosis, 193(2), 321–327. https://doi.org/10.1016/j.atherosclerosis.2006.09.015
    OpenUrlCrossRefPubMedWeb of Science
  54. 54.
    PanUKBTeam. (2020). Pan-UK Biobank https://pan.ukbb.broadinstitute.org.
  55. 55.↵
    Park, B., & Lee, Y. J. (2018). Metabolic syndrome and its components as risk factors for prolonged corrected QT interval in apparently healthy Korean men and women. J Clin Lipidol, 12(5), 1298–1304. https://doi.org/10.1016/j.jacl.2018.07.004
    OpenUrl
  56. 56.↵
    Peshkova, T., Nakashidze, I., Petrovic, N., Beridze, S., Kedelidze, N., Artmeladze, M., & Kamashidze, K. (2020). 73 Study of lipid profile and thyroid hormones in patients with uterine fibroids. International Journal of Gynecologic Cancer, 30(Suppl 3), A41–A41. https://doi.org/10.1136/ijgc-2020-IGCS.69
    OpenUrl
  57. 57.↵
    Phillips, A. N., & Smith, G. D. (1991). How independent are “independent” effects? Relative risk estimation when correlated exposures are measured imprecisely. J Clin Epidemiol, 44(11), 1223–1231. https://doi.org/10.1016/0895-4356(91)90155-3
    OpenUrlCrossRefPubMedWeb of Science
  58. 58.↵
    Purcell, S., Neale, B., Todd-Brown, K., Thomas, L., Ferreira, M. A. R., Bender, D., Maller, J., Sklar, P., de Bakker, P. I. W., Daly, M. J., & Sham, P. C. (2007). PLINK: a tool set for whole-genome association and population-based linkage analyses. American journal of human genetics, 81(3), 559–575. https://doi.org/10.1086/519795
    OpenUrlCrossRefPubMed
  59. 59.↵
    Qin, Z., Zhao, J., Geng, J., Chang, K., Liao, R., & Su, B. (2021). Higher Triglyceride–Glucose Index Is Associated With Increased Likelihood of Kidney Stones [Original Research]. Frontiers in Endocrinology, 12. https://doi.org/10.3389/fendo.2021.774567
  60. 60.
    R Core Team. (2021). R: A language and environment for statistical computing. In R Foundation for Statistical Computing. https://www.R-project.org/
  61. 61.↵
    Richardson, T. G., Sanderson, E., Palmer, T. M., Ala-Korpela, M., Ference, B. A., Davey Smith, G., & Holmes, M. V. (2020). Evaluating the relationship between circulating lipoprotein lipids and apolipoproteins with risk of coronary heart disease: A multivariable Mendelian randomisation analysis. PLOS Medicine, 17(3), e1003062–e1003062. https://doi.org/10.1371/journal.pmed.1003062
    OpenUrl
  62. 62.↵
    Rosenson, R. S., & Goonewardena, S. N. (2021). The Residual Risk Odyssey: From LDL to Lp(a)*. Journal of the American College of Cardiology, 78(5), 434–436. https://doi.org/https://doi.org/10.1016/j.jacc.2021.04.103
    OpenUrl
  63. 63.↵
    Rosenson, R. S., Shaik, A., & Song, W. (2021). New Therapies for Lowering Triglyceride-Rich Lipoproteins: JACC Focus Seminar 3/4. Journal of the American College of Cardiology, 78(18), 1817–1830. https://doi.org/https://doi.org/10.1016/j.jacc.2021.08.051
    OpenUrl
  64. 64.↵
    Roth, G. A., Mensah, G. A., Johnson, C. O., Addolorato, G., Ammirati, E., Baddour, L. M., Barengo, N. C., Beaton, A. Z., Benjamin, E. J., Benziger, C. P., Bonny, A., Brauer, M., Brodmann, M., Cahill, T. J., Carapetis, J., Catapano, A. L., Chugh, S. S., Cooper, L. T., Coresh, J., … Fuster, V. (2020). Global Burden of Cardiovascular Diseases and Risk Factors, 1990–2019: Update From the GBD 2019 Study. Journal of the American College of Cardiology, 76(25), 2982–3021. https://doi.org/https://doi.org/10.1016/j.jacc.2020.11.010
    OpenUrlCrossRefPubMed
  65. 65.↵
    Sanderson, E. (2021). Multivariable Mendelian Randomization and Mediation. Cold Spring Harbor Perspectives in Medicine, 11(2), a038984. https://doi.org/10.1101/cshperspect.a038984
    OpenUrlAbstract/FREE Full Text
  66. 66.↵
    Shaik, A., & Rosenson, R. S. (2021). Genetics of Triglyceride-Rich Lipoproteins Guide Identification of Pharmacotherapy for Cardiovascular Risk Reduction. Cardiovascular Drugs and Therapy, 35(3), 677–690. https://doi.org/10.1007/s10557-021-07168-0
    OpenUrl
  67. 67.↵
    Shih, C. W., Chen, Y. H., & Chen, W. L. (2022). Percentage of body fat is associated with increased risk of diverticulosis: A cross sectional study. PLOS ONE, 17(3), e0264746. https://doi.org/10.1371/journal.pone.0264746
    OpenUrl
  68. 68.↵
    Skrivankova, V. W., Richmond, R. C., Woolf, B. A. R., Yarmolinsky, J., Davies, N. M., Swanson, S. A., VanderWeele, T. J., Higgins, J. P. T., Timpson, N. J., Dimou, N., Langenberg, C., Golub, R. M., Loder, E. W., Gallo, V., Tybjaerg-Hansen, A., Davey Smith, G., Egger, M., & Richards, J. B. (2021). Strengthening the Reporting of Observational Studies in Epidemiology Using Mendelian Randomization: The STROBE-MR Statement. JAMA, 326(16), 1614–1621. https://doi.org/10.1001/jama.2021.18236
    OpenUrlPubMed
  69. 69.↵
    Smith, J. G., Luk, K., Schulz, C.-A., Engert, J. C., Do, R., Hindy, G., Rukh, G., Dufresne, L., Almgren, P., Owens, D. S., Harris, T. B., Peloso, G. M., Kerr, K. F., Wong, Q., Smith, A. V., Budoff, M. J., Rotter, J. I., Cupples, L. A., Rich, S., … Aging Research in Genetic Epidemiology Extracoronary Calcium Working, G. (2014). Association of Low-Density Lipoprotein Cholesterol–Related Genetic Variants With Aortic Valve Calcium and Incident Aortic Stenosis. JAMA, 312(17), 1764–1771. https://doi.org/10.1001/jama.2014.13959
    OpenUrlCrossRefPubMed
  70. 70.↵
    Sofat, R., Hingorani, A. D., Smeeth, L., Humphries, S. E., Talmud, P. J., Cooper, J., Shah, T., Sandhu, M. S., Ricketts, S. L., Boekholdt, S. M., Wareham, N., Khaw, K. T., Kumari, M., Kivimaki, M., Marmot, M., Asselbergs, F. W., van der Harst, P., Dullaart, R. P., Navis, G., … Casas, J. P. (2010). Separating the mechanism-based and off-target actions of cholesteryl ester transfer protein inhibitors with CETP gene polymorphisms. Circulation, 121(1), 52–62. https://doi.org/10.1161/circulationaha.109.865444
    OpenUrlAbstract/FREE Full Text
  71. 71.↵
    Sudlow, C., Gallacher, J., Allen, N., Beral, V., Burton, P., Danesh, J., Downey, P., Elliott, P., Green, J., Landray, M., Liu, B., Matthews, P., Ong, G., Pell, J., Silman, A., Young, A., Sprosen, T., Peakman, T., & Collins, R. (2015). UK Biobank: An Open Access Resource for Identifying the Causes of a Wide Range of Complex Diseases of Middle and Old Age. PLOS Medicine, 12(3), e1001779. https://doi.org/10.1371/journal.pmed.1001779
    OpenUrl
  72. 72.↵
    Tonoyan, N. M., Chagovets, V. V., Starodubtseva, N. L., Tokareva, A. O., Chingin, K., Kozachenko, I. F., Adamyan, L. V., & Frankevich, V. E. (2021). Alterations in lipid profile upon uterine fibroids and its recurrence. Scientific Reports, 11(1), 11447. https://doi.org/10.1038/s41598-021-89859-0
    OpenUrl
  73. 73.
    Triglyceride Coronary Disease Genetics, C., Emerging Risk Factors, C., Sarwar, N., Sandhu, M. S., Ricketts, S. L., Butterworth, A. S., Di Angelantonio, E., Boekholdt, S. M., Ouwehand, W., Watkins, H., Samani, N. J., Saleheen, D., Lawlor, D., Reilly, M. P., Hingorani, A. D., Talmud, P. J., & Danesh, J. (2010). Triglyceride-mediated pathways and coronary disease: collaborative analysis of 101 studies. Lancet (London, England), 375(9726), 1634–1639. https://doi.org/10.1016/S0140-6736(10)60545-4
    OpenUrl
  74. 74.↵
    Uimari, O., Auvinen, J., Jokelainen, J., Puukka, K., Ruokonen, A., Järvelin, M.-R., Piltonen, T., Keinänen-Kiukaanniemi, S., Zondervan, K., Järvelä, I., Ryynänen, M., & Martikainen, H. (2016). Uterine fibroids and cardiovascular risk. Human Reproduction, 31(12), 2689–2703. https://doi.org/10.1093/humrep/dew249
    OpenUrlCrossRefPubMed
  75. 75.↵
    Varbo, A., Benn, M., Tybjærg-Hansen, A., Jørgensen, A. B., Frikke-Schmidt, R., & Nordestgaard, B. G. (2013). Remnant cholesterol as a causal risk factor for ischemic heart disease. J Am Coll Cardiol, 61(4), 427–436. https://doi.org/10.1016/j.jacc.2012.08.1026
    OpenUrlFREE Full Text
  76. 76.↵
    Verbanck, M., Chen, C.-Y., Neale, B., & Do, R. (2018). Detection of widespread horizontal pleiotropy in causal relationships inferred from Mendelian randomization between complex traits and diseases. Nature Genetics, 50(5), 693–698. https://doi.org/10.1038/s41588-018-0099-7
    OpenUrlCrossRefPubMed
  77. 77.↵
    White, J., Swerdlow, D. I., Preiss, D., Fairhurst-Hunter, Z., Keating, B. J., Asselbergs, F. W., Sattar, N., Humphries, S. E., Hingorani, A. D., & Holmes, M. V. (2016). Association of Lipid Fractions With Risks for Coronary Artery Disease and Diabetes. JAMA cardiology, 1(6), 692–699. https://doi.org/10.1001/jamacardio.2016.1884
    OpenUrl
  78. 78.↵
    Willer, C. J., Schmidt, E. M., Sengupta, S., Peloso, G. M., Gustafsson, S., Kanoni, S., Ganna, A., Chen, J., Buchkovich, M. L., Mora, S., Beckmann, J. S., Bragg-Gresham, J. L., Chang, H. Y., Demirkan, A., Den Hertog, H. M., Do, R., Donnelly, L. A., Ehret, G. B., Esko, T., … Abecasis, G. R. (2013). Discovery and refinement of loci associated with lipid levels. Nat Genet, 45(11), 1274–1283. https://doi.org/10.1038/ng.2797
    OpenUrlCrossRefPubMed
  79. 79.
    Yu, X., Wang, T., Huang, S., & Zeng, P. (2021). Evaluation of the causal effects of blood lipid levels on gout with summary level GWAS data: two-sample Mendelian randomization and mediation analysis. Journal of Human Genetics, 66(5), 465–473. https://doi.org/10.1038/s10038-020-00863-0
    OpenUrl
  80. 80.↵
    Zewinger, S., Reiser, J., Jankowski, V., Alansary, D., Hahm, E., Triem, S., Klug, M., Schunk, S. J., Schmit, D., Kramann, R., Körbel, C., Ampofo, E., Laschke, M. W., Selejan, S.-R., Paschen, A., Herter, T., Schuster, S., Silbernagel, G., Sester, M., … Speer, T. (2020). Apolipoprotein C3 induces inflammation and organ damage by alternative inflammasome activation. Nature Immunology, 21(1), 30–41. https://doi.org/10.1038/s41590-019-0548-1
    OpenUrlPubMed
Back to top
Posted July 21, 2022.
Download PDF
Data/Code
Email
Phenome-wide Mendelian randomization study of plasma triglycerides and 2,600 disease traits
Joshua K. Park, Shantanu Bafna, Iain S. Forrest, Áine Duffy, Carla Marquez-Luna, Ben O. Petrazzini, Ha My Vy, Daniel M. Jordan, Marie Verbanck, Jagat Narula, Robert S. Rosenson, Ghislain Rocheleau, Ron Do
medRxiv 2022.07.21.22277900; doi: https://doi.org/10.1101/2022.07.21.22277900
Twitter logo Facebook logo LinkedIn logo Mendeley logo
Citation Tools

Subject Area