The relationship between circulating lipids and breast cancer risk: A Mendelian randomization study

Abstract

[...]they analyzed each lipid trait separately rather than take advantage of multivariable methods to consider lipid traits together along with additional, potentially confounding causal risk factors. [...]the authors did not quantitatively assess heterogeneity to determine whether the observed lipid associations were statistically different across BC subtypes. [...]an approach obviates the need to remove pleiotropic variants and the loss of statistical power that results from this removal. [...]an approach that considers the effects of genetic variants on known risk factors for BC, such as body mass index (BMI) and age at menarche [14–20], could increase power. [...]we perform genetic correlation analyses to look for both genome-wide and locus-based correlation in effect sizes between lipids and BC.

Introduction

Breast cancer (BC) is the second leading cause of death for women, motivating the need for a better understanding of its etiology and more effective treatments [1]. Cholesterol is a known risk factor for multiple diseases that have reported associations with BC, including obesity, heart disease, and diabetes [2]. However, it is unknown whether cholesterol plays a causal role in BC susceptibility.

The body of epidemiological and clinical trial studies to date has yet to determine whether there is a causal relationship between cholesterol and BC. Observational epidemiological studies have reported positive, negative, or no relationship between lipid levels and BC risk [3–6]; however, these studies can suffer from confounding. A comprehensive meta-analysis found evidence that statin use may reduce BC risk [7], and cholesterol-lowering medications have been associated with improved outcomes in BC patients on hormonal therapy, suggesting an interaction of circulating cholesterol levels with estrogen-sensitive breast tissues [8]. These mixed findings motivate the need for a high-powered causal inference analysis of lipids on BC.

To try to resolve these discrepancies, recent studies have applied the framework of Mendelian randomization (MR) to determine whether genetically elevated lipid levels associate with BC risk. In a small sample of 1,187 BC cases, Orho-Melander and colleagues used multivariable MR to find suggestive evidence of a relationship between both triglycerides and HDL (high-density lipoprotein) cholesterol and BC, but no association between LDL (low-density lipoprotein) cholesterol and cancer [9]. In a second study, Nowak and colleagues [10] performed an MR analysis with genetic association data from large genome-wide association studies (GWASs) for lipids and BC [11,12]. They reported nominal positive associations between LDL-cholesterol levels and all BCs and between HDL-cholesterol levels and ER (estrogen receptor)-positive BCs. While compelling, this study also had limitations. First, they used relatively few variants in their genetic instrument because of the removal of pleiotropic variants in order to address confounding due to pleiotropy, resulting in a conservative analysis. Second, they analyzed each lipid trait separately rather than take advantage of multivariable methods to consider lipid traits together along with additional, potentially confounding causal risk factors. Third, the authors did not quantitatively assess heterogeneity to determine whether the observed lipid associations were statistically different across BC subtypes. Another recent study by Qi and Chatterjee applied a newly developed MR method and reported an association between HDL-cholesterol and BC that they defined as borderline statistically significant [13]. Like Nowak and colleagues, this paper also does not explicitly include correlated risk factors in their analysis and did not stratify BC by ER status.

These studies motivate an MR study that considers multiple lipid traits concurrently to delineate the independent effect of each lipid trait on BC susceptibility. Such an approach obviates the need to remove pleiotropic variants and the loss of statistical power that results from this removal. Therefore, an approach that considers the effects of genetic variants on known risk factors for BC, such as body mass index (BMI) and age at menarche [14–20], could increase power.

While MR assesses evidence for a causal relationship, genome-wide genetic correlation analysis determines whether 2 traits simply have a shared genetic basis. Local genetic correlation analyses test whether 2 traits have shared heritability that is localized to specific genomic regions. These loci may then harbor causal variants and genes that contribute to heritability of both traits. Jiang and colleagues recently estimated genome-wide genetic correlation between lipid traits and BC risk [21]. This study did not find a statistically significant association between any lipid trait and BC using lipid summary statistics, though a previous study with a smaller BC GWAS sample size did report a nominally significant (P < 0.05) negative genetic correlation between triglycerides and BC risk [16]. Both these studies used the same method to estimate genome-wide genetic correlations [15], and neither tested for local genetic correlations between lipid traits and BC risk.

In what follows, we apply the causal inference framework of MR to determine whether genetically elevated lipid traits modify BC susceptibility (both all BC and BC stratified by ER status). We take advantage of a recent GWAS for lipid levels performed in up to 215,551 individuals of European ancestry [22], which has not been previously applied to MR studies of BC, to provide power for our causal inference analyses. We utilize several MR techniques, including single-exposure, multivariable, and gene-specific approaches. Of chief concern in modern MR studies, including prior studies of BC and lipids, is confounding due to pleiotropy. For instance, a genetic variant may affect lipid levels indirectly through some other biomarker. If this biomarker directly affects BC risk, this could confound the MR analysis and cause an incorrect inference of a causal effect of lipid levels on BC. Our gene-specific approach utilizes only genetic variants near core HDL pathway genes to minimize this concern. Additionally, we use a multivariable approach that assesses the effects of lipid traits independent of one another and of BMI and age at menarche. Finally, we perform genetic correlation analyses to look for both genome-wide and locus-based correlation in effect sizes between lipids and BC.

Methods

Study populations

Lipid GWAS summary statistics were obtained from the Million Veteran Program (MVP) (up to 215,551 European individuals) [22] and the Global Lipids Genetics Consortium (GLGC) (up to 188,577 genotyped individuals) [12]. As additional exposures in multivariable MR analyses, we used BMI summary statistics from a meta-analysis of GWASs in up to 795,640 individuals and age at menarche summary statistics from a meta-analysis of GWASs in up to 329,345 women of European ancestry [17,23]. GWAS summary statistics from 122,977 BC cases and 105,974 controls were obtained from the Breast Cancer Association Consortium (BCAC) [11]. The MVP received ethical and study protocol approval from the Veteran Affair Central Institutional Review Board in accordance with the principles outlined in the Declaration of Helsinki, and written consent was obtained from all participants. For the Willer and colleagues [12] and BCAC [11] data sets, we refer the reader to the primary GWAS manuscripts and their supplementary material for details on consent protocols for each of their respective cohorts. More details on these cohorts are in the S1 Text.

Lipid meta-analysis

We performed a fixed-effects meta-analysis between each lipid trait (Total cholesterol [TC], LDL, HDL, and triglycerides [TGs]) in GLGC and the corresponding lipid trait in the MVP cohort [12,22] using the default settings in PLINK [24]. There is some genomic inflation in these meta-analysis association statistics, but linkage disequilibrium (LD)-score regression intercepts demonstrate that this inflation is in large part due to polygenicity and not population stratification (S1 Fig).

MR analyses

MR analyses were performed using the TwoSampleMR R package version 0.4.13 (https://github.com/MRCIEU/TwoSampleMR) [25]. For all analyses, we used a 2-sample MR framework, with exposure(s) (lipids, BMI, age at menarche) and outcome (BC) genetic associations from separate cohorts. Unless otherwise noted, MR results reported in this manuscript used inverse-variance weighting assuming a multiplicative random effects model. For single-trait MR analyses, we additionally employed Egger regression [26], weighted median [27], and mode-based [28] estimates. SNPs associated with each lipid trait were filtered for genome-wide significance (P < 5 × 10−8) from the MVP lipid study [12], and then we removed SNPs in LD (r2 < 0.001 in UK10K consortium) [29] in order to obtain independent variants. All genetic variants were harmonized using the TwoSampleMR harmonization function with default parameters. Each of these independent, genome-wide significant SNPs was termed a genetic instrument. We estimated that these single-trait MR genetic instruments had 80% power to reject the null hypothesis, with a 1% error rate, for the following odds ratio (OR) increases in BC risk due to a standard deviation increase in lipid levels: HDL, 1.057; LDL, 1.058; TGs, 1.055; TC, 1.060 [30,31]. We tested for directional pleiotropy using the MR-Egger regression test [26]. To reduce heterogeneity in our genetic instruments for single-trait MR, we employed a pruning procedure (S1 Text). Genetic instruments used in single-trait MR are listed in S1 Table. For multivariable MR experiments [32,33], we generated genetic instruments by first filtering the genotyped variants for those present across all data sets. For each trait and data set combination (Yengo and colleagues [23] for BMI; Day and colleagues for age at menarche [17]; MVP and GLGC for HDL, LDL, and TGs), we then filtered for genome-wide significance (P < 5 × 10−8) and for linkage disequilibrium (r2 < 0.001 in UK10K consortium) [29]. We performed tests for instrument strength and validity [34], and each multivariable MR experiment had sufficient instrument strength. We removed variants driving heterogeneity in the ratio of outcome/exposure effects causing instrument invalidity (S1 Text). Genetic instruments used in multivariable MR are listed in S2 Table. Because the MR methods and tests we employed are highly correlated, we did not apply a multiple testing correction to the reported P-values.

Core HDL and LDL pathway genetic instrument development

We defined sets of core genes for HDL or LDL that met the following criteria: (1) their protein products are known to play a key role in HDL or LDL biology (plus HMGCR and NPC1L1, 2 targets of LDL-lowering drugs, in the LDL gene set), and (2) there were conditionally independent lipid trait-associated variants within 100 kb upstream or downstream of the RefSeq coordinates for the gene (or locus, in the case of Apolipoprotein E [APOE]-Apolipoprotein C [APOC]1-APOC4-APOC2 and Apolipoprotein A [APOA]4-APOC3-APOA1) [22]. We then used the conditional HDL or LDL association statistics from Klarin and colleagues for those genes in gene-specific MR analyses [22]. The loci included in each set and the genetic instruments used in each locus-specific MR are listed in S3 Table. We performed a separate fixed-effects inverse-variance weighted MR with the conditionally independent genetic instruments at each gene and performed fixed-effects inverse-variance weighted meta-analysis of the results across HDL or LDL genes using the R package meta [35].

Genetic correlation analyses

We performed cross-trait LD-score regression using the LDSC toolkit, available at https://github.com/bulik/ldsc, with default parameters [15], with the BCAC association statistics for BC and our meta-analysis of GLGC and MVP for lipid associations. In addition, we used the ρ-Hess software both for whole-genome genetic correlation and for local genetic correlation analysis [36], using the UK10K reference panel, and the LD-independent loci published in Berisa and colleagues to partition the genome [37]. We used a Bonferroni significance threshold based on the number of these independent loci (1,704 loci). There was minor cohort overlap between the GLGC and BC GWASs because of the EPIC cohort [10]. We included this overlap when performing ρ-Hess, available at https://huwenboshi.github.io/hess/, using the cross-trait LD-score intercept to estimate phenotypic correlation [36]. The association of the lead BC and lipid SNPs at the ABO locus was obtained using the GTEx v8 data set [38].

Analysis plan

Our study did not develop a prospective analysis plan. We began by testing for a potential causal relationship between lipids and BC risk using single-trait two-sample MR with lipid genetic associations from the GLGC. After this experiment showed a significant relationship, we tested whether it persisted when we corrected for correlated phenotypes with multivariable MR. Following significant results with the GLGC summary statistics, we decided to confirm these findings using genetic associations from the larger MVP cohort. After our results with MVP confirmed our initial findings, we performed additional MR sensitivity analyses and locus-specific MR. In parallel to our MR experiments, we measured the cross-trait and local genetic correlations between these traits. This study is reported as per the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) and STROBE-MR guidelines (S1 STROBE Checklist, S1 STROBE-MR Checklist) [39,40].

Results

Single-trait MR in BC

We first performed single-trait MR analyses using summary statistics from MVP [22] for each of 4 lipid traits (i.e., TC, HDL-cholesterol, LDL-cholesterol, and TGs) as the intermediate biomarkers and risk for all BCs as the outcome (S2 Fig). We observed a significant relationship between genetically elevated HDL and BC risk (OR = 1.10 per standard deviation of lipid level increase, 95% confidence interval [CI] = 1.04–1.17, P = 2.1 x 10−3) and genetically decreased TG levels and BC risk (OR = 0.93, 95% CI = 0.88–0.99, P = 0.015; S4 Table). Sensitivity analyses identified heterogeneity (Methods, S5 Table), but there was no evidence of bias from directional pleiotropy (Methods, S6 Table). To mitigate concerns of instrument heterogeneity, we removed variants from our genetic instrument for each lipid trait that were responsible for instrument heterogeneity (S1 Text) and again observed a relationship with HDL-cholesterol (OR = 1.08, 95% CI = 1.04–1.13, P = 7.4 × 10−5) and TGs (OR = 0.94, 95% CI = 0.90–0.98, P = 2.6 × 10−3) (Fig 1, S3 and S4 Figs, S7 Table). Because HDL and TGs are inversely correlated [15,41], the opposing relationship between these 2 lipid traits and BC could be expected in single-trait analyses.

[Figure omitted. See PDF.]

Fig 1. Results of MR analyses of the effects of lipids on BC risk.

Results plotted are after pruning for instrument heterogeneity. The forest plot on the right displays the OR of the effect of a 1-standard–deviation increase in genetically determined HDL-cholesterol on BC risk as a diamond, with error bars representing the 95% CI. The vertical dotted line delineates an OR of 1, i.e., no effect of the exposure on BC risk. BC, breast cancer; BMI, body mass index; CI, confidence interval; HDL, high-density lipoprotein; IVW, inverse-variance weighted MR; LDL, low-density lipoprotein; MR, Mendelian randomization; MVMR, multivariable MR; OR, odds ratio; TG, triglyceride.

https://doi.org/10.1371/journal.pmed.1003302.g001

To confirm that our results using lipid genetic associations from MVP were not due to heterogeneity between data sets, we also tested the relationship between lipid traits and BC using a meta-analysis of the 2 major lipid GWASs from MVP and GLGC and from GLGC alone. Overall, single-trait MR analyses with the meta-analysis and GLGC lipid associations produced consistent results to those with MVP alone (S5 Fig). In a reciprocal single-trait MR testing the effect of genetically determined BC risk on each lipid trait, we observed no relationship with HDL- or LDL-cholesterol (S8 Table) but did see a relationship with TGs. However, a Steiger test for directionality confirmed that using BC as the outcome was the correct causal direction for all lipid traits (S7 Table) [42]. We also performed genetic instrument pruning in the same manner as Nowak and colleagues: removing genetic instruments for LDL, HDL, and TGs that were associated with at least one of the 2 other lipid traits (P < 0.001) [10]. After this pruning, we did not find a significant relationship with LDL, HDL, or TG, and we note that this pruning procedure resulted in considerably larger CIs spanning OR = 1 for all traits, with reversed direction of effect estimates for LDL and TGs (S9 Table).

Multivariable MR with age at menarche and BMI as exposures

It has been previously observed that BMI and age at menarche are both genetically correlated and epidemiologically associated with both BC [20,43,44] and lipid traits [15,41]. To incorporate these potential confounders into our causal inference framework, we performed multivariable MR analyses using all 3 lipid traits (genetic effect estimates from MVP), age at menarche, and BMI as exposures and BC risk as the outcome (Fig 1). We observed significant relationships between genetically influenced HDL, LDL, BMI, and age at menarche with BC (HDL: OR = 1.06, 95% CI = 1.03–1.10, P = 4.93 × 10−4; LDL: OR = 1.04, 95% CI = 1.01–1.07, P = 0.02; BMI: OR = 0.90, 95% CI = 0.87–0.94, P = 1.15 × 10−6; age at menarche: OR = 0.96, 95% CI = 0.93–0.99, P = 2.44 × 10−3), but not TGs (OR = 0.98, 95% CI = 0.95–1.00, P = 0.10) (Fig 1, S10 Table). Our results were consistent before and after pruning for genetic instrument heterogeneity (S10 Table) and when using summary statistics from 3 independent subsets of the BC data set (S6 Fig, S10 Table). We also performed multivariable MR with pairs of lipid traits with genetic effect estimates from different data sets (GLGC or MVP), with and without BMI, and saw consistent results (S11 Table, S7 Fig). Considering the genetic correlation between HDL and TGs, the greater significance of the HDL association compared with the TG association with BC in multivariable analysis, and the consistent relationship between HDL and BC across BC data sets, we focused our further MR analyses on the relationship between HDL-cholesterol and BC, in addition to the previously reported association between LDL and BC [10].

MR with outcome stratified by ER status

We next performed an MR analysis of the relationship between genetically influenced lipids and BC risk stratified by ER-positive (ER+) or negative (ER‒) status. We observed similar effect size estimates of the 4 lipid traits on the BC subtypes as on BC not stratified by subtype (S8 Fig). A formal test for heterogeneity found no evidence to reject the null hypothesis of homogeneity between the cancer subtypes (e.g., HDL: Cochran’s Q = 6.6 × 10−5, P = 0.99; S12 Table). Thus, we observed no substantive difference in the relationship from any lipid trait to ER+ or ER− BCs, consistent with the strong genetic correlation between these 2 BC subtypes (cross-trait LD-score regression genetic correlation estimate = 0.62, P = 2.9 × 10−83). When we used ER+ or ER‒ BCs as the outcome in multivariable MR, we also saw consistent effects as the analysis with all BCs as the outcome (S13 Table, S9 Fig).

HDL and LDL pathway gene-specific MR

We next examined associations for BC risk at genetic variants near core HDL or LDL genes. We identified conditionally independent associations at genes that were previously annotated with a core role in the metabolism of each lipid trait or an established drug target (HDL: ATP Binding Cassette Subfamily A Member 1 [ABCA1], APOA4-APOC3-APOA1, APOE-APOC1-APOC4-APOC2, Cholesteryl Ester Transfer Protein [CETP], Lecithin-Cholesterol Acyltransferase [LCAT], Lipase C Hepatic Type [LIPC], Lipase G Endothelial Type [LIPG], Phospholipid Transfer Protein [PLTP], Scavenger Receptor Class B Member 1 [SCARB1]; LDL: Apolipoprotein B [APOB], 3-Hydroxy-3-Methylglutaryl-CoA Reductase [HMGCR], LDL Receptor [LDLR], Lipoprotein(A) [LPA], Myosin Regulatory Light Chain Interacting Protein [MYLIP], NPC1-Like Intracellular Cholesterol Transporter 1 [NPC1L1], Proprotein Convertase Subtilisin/Kexin Type 9 [PCSK9]) (Methods, S3 Table). For each gene or locus with at least 2 conditionally independent genetic instruments (all except LCAT and MYLIP), we performed inverse-variance-weighted MR (fixed-effects model) with conditional HDL or LDL effect size estimates as the exposure and BC risk as the outcome (S10 and S11 Figs). We observed significant (P < 0.05) positive relationships between HDL and BC risk at 3 loci (ABCA1, APOE-APOC1-APOC4-APOC2, and CETP; Fig 2), and between LDL and BC risk at 1 locus (HMGCR, S12 Fig). Combining the effect estimates across core genes in a meta-analysis, we observed a significant positive relationship for HDL (OR = 1.11, 95% CI = 1.06–1.16, P < 0.001; Fig 2) and LDL (OR = 1.07, 95% CI = 1.01–1.14, P = 0.02; S12 Fig). There was no evidence of heterogeneity across loci in either meta-analysis (HDL: Q = 6.63, P = 0.47; LDL: Q = 5.53, P = 0.35).

[Figure omitted. See PDF.]

Fig 2. MR results for HDL gene-specific instruments and meta-analysis of effect estimates across genes.

The forest plot on the right displays the OR of the effect of a 1-standard–deviation increase in genetically determined HDL-cholesterol for each locus on BC risk as a diamond, and the error bars represent the 95% CI of the effect estimate. The vertical dotted line delineates an OR of 1, i.e., no effect of the exposure on BC risk. For HDL gene-specific instruments, see S3 Table. ABCA1, ATP Binding Cassette Subfamily A Member 1; APOC, Apolipoprotein C; APOE, Apolipoprotein E; BC, breast cancer; CETP, Cholesteryl Ester Transfer Protein; CI, confidence interval; HDL, high-density lipoprotein; LIPC, Lipase C, Hepatic Type; LIPG, Lipase G, Endothelial Type; MR, Mendelian randomization; N SNPs, number of genetic instruments included in each locus’s MR analysis; OR, odds ratio; PLTP, Phospholipid Transfer Protein; SCARB1, Scavenger Receptor Class B Member 1.

https://doi.org/10.1371/journal.pmed.1003302.g002

Genome-wide and local genetic correlation

If cholesterol levels were a causal risk factor for BC, we might expect a correlation between the strength of genetic association with these 2 traits at genetic variants across the entire genome in addition to those at genome-wide significant loci. To answer this question, we utilized 2 approaches to estimate genetic correlation between BC and lipid traits. Using the ρ-Hess method, we found significant (P < 0.05) correlations between LDL (P < 0.001) and TC (P = 0.01) and BC, with directions consistent with our MR results (S14 Table). Cross-trait LD-score regression found positive genetic correlation estimates for TC, LDL, and HDL and a negative estimate for TGs (S13 Fig, S14 Table), consistent with our MR and ρ-Hess results [15]. However, the only significant association (P < 0.05) was with TC and ER-negative BC (P = 0.04).

To discover new loci that are enriched for genetic correlation between BC and lipids, we used the ρ-Hess method, which detects genomic regions harboring high genetic correlation between 2 traits [36]. ρ-Hess identified one region that surpassed Bonferroni test correction, with a positive correlation between both LDL and TC and BC (S15 Table). In this region, there are 2 SNPs in high LD (rs532436 and rs635634, r2 = 0.99) that are genome-wide significantly associated with LDL (rs532436: P < 0.001), TC (P < 0.001), and BC (P < 0.001). These SNPs lie within an intron of the ABO blood group determining ABO gene. rs635634 moderately tags an SNP associated with ABO blood type [41]. However, this SNP is also associated with a change in gene expression of ABO in multiple tissues (P < 0.001 in breast mammary tissue) [38].

Discussion

Using MR, we provide evidence that genetically elevated HDL and LDL levels are associated with increased risk for BC, supporting a causal hypothesis. Previous meta-analyses of observational studies of BC risk and lipids reported a negative association with HDL and no relationship with LDL [4,5], whereas individual studies have reported a positive relationship with HDL [6] or no relationship with HDL or LDL [45,46]. Our analyses help clarify these mixed results and infer a direction of effect, which is not possible in observational studies because of potential reverse causation. Furthermore, we find evidence of genome-wide genetic correlation between some lipid traits and BC and local genetic correlation at the ABO locus. Although some studies have found an association between blood group and BC risk [47], haplotype patterns indicate that ABO gene expression, not blood group, may be the causal mechanism [37]. However, because of the pleiotropic nature of the ABO locus, it is unclear whether the BC association is caused by the lipid associations [41].

Although Nowak and colleagues previously used MR to discover associations between lipids and BC [10], our report presents a thorough reconsideration of these effects. Even after conditioning on the effects of HDL, BMI, and age at menarche, our MR analysis suggests a potential causal relationship between LDL and BC. Nowak and colleagues only found a relationship between HDL-cholesterol and ER+ BC, whereas we found a relationship between HDL-cholesterol and risk for all BCs. We also find a previously unreported association with TGs and BC, though our multivariable analysis suggests this may be explained by correlation between TGs and HDL and not an independent TG effect. In their analyses, Nowak and colleagues used a strict pruning procedure to isolate the effects of each lipid trait. However, this approach reduces power because of the high genetic correlation of these traits. The multivariable approach taken here is an alternative way to estimate the effect of an exposure while accounting for correlated exposures.

Our results largely agree with those reported in the recent MR study of BC and lipids by Beeghly-Fadiel and colleagues [48], published while this manuscript was under peer review. Both studies use multiple types of MR analyses, including approaches accounting for confounding by BMI and age at menarche, and report a positive association between HDL and BC and a negative association between TGs and BC. However, our report provides complementary analysis and data that support the central findings of both pieces of work. First, we took advantage of the recently reported MVP lipids GWAS [22], providing a larger number of genetic instruments for all lipid traits considered. Second, we explicitly considered age of menarche and BMI in multivariable models with all lipid traits. It is crucial to consider the collection of each of these risk factors together to estimate a causal effect estimate that is independent of these confounders, as well as across cancer risk strata (i.e., ER status). While Beeghly-Fadiel and colleagues took advantage of access to individual-level data to adjust for confounding factors in their single lipid MR, these confounders were not corrected for in their multivariable MR. Third, while both studies are consistent in their relationship between HDL and BC, we reported a nominal association (P < 0.05) with LDL levels when considering all risk and confounding factors jointly. Finally, we present unique, locus-specific MR analyses to show that conditionally independent associations at single loci implicated in HDL or LDL biology are significantly associated with BC risk.

A challenge of MR analyses that use hundreds of genetic instruments is that we do not know the mechanism of action of these instruments on the exposure trait. By focusing on loci with mechanistic evidence of a direct effect on lipid levels, we can remove uncertainty about potential pleiotropic effects on BC risk. Thus, the significant relationships observed in our locus-specific MR analyses across HDL or LDL pathway genes provide additional evidence for a direct effect of increased HDL or LDL levels on BC risk. Furthermore, these genes represent potential or established drug targets, and each locus-specific MR experiment provides preliminary evidence for the therapeutic potential of cholesterol modification on BC prevention.

Substantial effort has been spent developing HDL-raising therapies for cardiovascular disease prevention; however, recent studies have proposed an increase in all-cause mortality in individuals with high HDL levels [49,50]. Our results suggest that therapies that aim to reduce cardiovascular risk by raising HDL levels might have an unintended consequence of elevated BC risk. Specifically, our gene-based score using HDL-raising variation at the CETP locus predicted that CETP-based inhibition would elevate BC risk (OR = 1.11, 95% CI = 1.04–1.18, P < 0.001). Additionally, 2 recent MR studies reported causal evidence between elevated HDL and risk for age-related macular degeneration [51,52]. These potential disease-increasing consequences may not have been possible to identify in safety trials, given the limited window of study to monitor progression or incidence of disease, the putative causal effect estimates, and the demographics of the study population (i.e., a higher proportion of male participants) [53]. Our result supports the use of human genetics data as both a novel strategy for therapeutic targeting and for the discovery of potential drug complications to direct long-term post-clinical–trial follow-up [54].

We note several caveats to our analyses. The first is that MR makes several assumptions that must be met for accurate causal inference [55,56]. Although we used statistical methods that try to detect and correct for violations of these assumptions, these methods are not guaranteed to correct for all types of confounding, and alternative causal inference frameworks outside of MR are warranted. Secondly, MR is only able to make inferences about trait associations in the populations from which the GWASs are derived. We are unaware of evidence that BC or lipid genetic architecture varies significantly across populations, but if this was the case, our findings may not be generalizable to these different scenarios. However, the concordance between our results using the MVP and GLGC GWASs mitigate this concern with regards to potential heterogeneity in lipid genetic architecture. Thirdly, the estimated lipid/BC effect sizes represent only the population-averaged causal effect and may not generalize well to other populations or settings [57]. We note that our effect estimates may be attenuated because of association of lipid instrumental variables with the use of lipid-lowering medication, and that we cannot be certain that the true underlying causal exposure is lipid levels rather than another phenotype for which lipids are a proxy. However, we are not aware of any process for which lipids is a proxy through which BC would be affected, and our gene-specific approach minimizes this concern. Additionally, it is perhaps surprising that we did not find a significant genetic correlation between BC and lipids using cross-trait LD-score regression; however, our result corroborates a previous study that performed this analysis using smaller GWASs [16]. Our lack of significant results could be caused by limited polygenicity of either trait, which decreases the power of this method [15]. We do find significant cross-trait heritability between BC and 2 lipid traits (TC and LDL) using the ρ-Hess method. The discrepancy between LDSC and ρ-Hess may be due to a difference in the statistical power of these methods that has been previously reported [36].

The analyses presented here do not bring evidence on a specific mechanism for tumorigenesis, but they do bring renewed attention to potential mechanisms requiring future functional study. Cholesterol and its oxysterol metabolites, either in the circulatory system or in the local mammary microenvironment, may have direct effects on mammary tissue growth induction of breast tumorigenesis [58,59].

These findings support a causal relationship between increased HDL-cholesterol and increased BC risk, and this hypothesis warrants further exploration. Statins are widely prescribed to decrease LDL levels; however, statins also increase HDL levels. If further research substantiates the relationship between higher HDL levels and increased BC risk, the consensus that HDL is “good cholesterol,” or of benign effect, may require re-evaluation.

Supporting information

S1 STROBE checklist. Reporting document following the STROBE guidelines for our study.

STROBE, Strengthening the Reporting of Observational Studies in Epidemiology.

https://doi.org/10.1371/journal.pmed.1003302.s001

(DOCX)

S1 STROBE-MR checklist. Reporting document following the preliminary STROBE-MR guidelines for our study.

MR, Mendelian randomization; STROBE, Strengthening the Reporting of Observational Studies in Epidemiology.

https://doi.org/10.1371/journal.pmed.1003302.s002

(DOCX)

S1 Fig. QQ plots for lipid association statistics.

Generated from a meta-analysis of Klarin and colleagues [22] and Willer and colleagues [12]. LD-score regression intercepts and standard errors from the meta-analysis association statistics were as follows: TC, 1.1293 (0.1143); TG, 1.0317 (0.0656); LDL, 1.0933 (0.1001); HDL, 1.1715 (0.0758); see also S14 Table. HDL, high-density lipoprotein; LD, linkage disequilibrium; LDL, low-density lipoprotein; TC, total cholesterol; TG, triglyceride; λGC, genomic inflation factor.

https://doi.org/10.1371/journal.pmed.1003302.s003

(PDF)

S2 Fig. Scatter plots of unpruned single-trait MR genetic instruments’ effect estimates on exposure and outcome.

Plotted are the genetic instruments included in unpruned single-trait MR analyses. Each plot contains effect estimates from MVP for one of 4 lipid traits (HDL, LDL, TC, TGs) on the x-axis and effect estimates for risk of all BCs (allBC) on the y-axis. Error bars represent the 95% CI, and regression lines represent the slope estimate from one of 3 MR tests: IVW (light blue), Egger regression (dark blue), and weighted median (green). BC, breast cancer; CI, confidence interval; HDL, high-density lipoprotein; IVW, inverse-variance weighted; LDL, low-density lipoprotein; MR, Mendelian randomization; TC, total cholesterol; TG, triglyceride.

https://doi.org/10.1371/journal.pmed.1003302.s004

(PDF)

S3 Fig. Scatter plots of pruned single-trait MR genetic instruments’ effect estimates on exposure and outcome.

Scatter plots of genetic instruments included in pruned single-trait MR analyses. Genetic instruments were pruned to pass heterogeneity test. Each plot contains effect estimates from MVP for one of 4 lipid traits (HDL, LDL, TC, TG) on the x-axis and effect estimates for risk of all BCs on the y-axis. Error bars represent the 95% CI, and regression lines represent the slope estimate from one of 3 MR tests: IVW (light blue), Egger regression (dark blue), and weighted median (green). BC, breast cancer; CI, confidence interval; HDL, high-density lipoprotein; IVW, inverse-variance weighted; LDL, low-density lipoprotein; MR, Mendelian randomization; TC, total cholesterol; TG, triglyceride.

https://doi.org/10.1371/journal.pmed.1003302.s005

(PDF)

S4 Fig. Results of single-trait MR analyses with a lipid trait as the exposure and risk for all BCs as the outcome.

The exposure association data for this analysis were from MVP. Genetic instruments were pruned to pass heterogeneity test. Error bars represent the 95% CI. Estimates were calculated using the IVW method. *P < 0.05; **P < 0.001. See S7 Table for ORs, CIs, and P-values. BC, breast cancer; CI, confidence interval; HDL, high-density lipoprotein; IVW, inverse-variance weighted; LDL, low-density lipoprotein; MR, Mendelian randomization; MVP, Million Veteran Program; OR, odds ratio; TC, total cholesterol; TG, triglyceride.

https://doi.org/10.1371/journal.pmed.1003302.s006

(PDF)

S5 Fig. Single-trait MR with lipid association statistics from MVP, GLGC, or MVP + GLGC meta-analysis.

Genetic association statistics for all BCs were used for the outcome. Genetic instruments were pruned to pass heterogeneity test. Error bars represent the 95% CI. Estimates were calculated using the IVW method. BC, breast cancer; CI, confidence interval; GLGC, Global Lipids Genetics Consortium; HDL, high-density lipoprotein; IVW, inverse-variance weighted; LDL, low-density lipoprotein; MR, Mendelian randomization; MVP, Million Veteran Program; TC, total cholesterol; TG, triglyceride.

https://doi.org/10.1371/journal.pmed.1003302.s007

(PDF)

S6 Fig. Multivariable MR analyses stratified by 3 independent subsets of the BCAC data set.

Results of multivariable MR analyses with 3 lipid traits (HDL, LDL, TG), BMI, and AaM as exposures and BC risk as the outcome. Each panel presents multivariable MR results using BC summary statistics from an independent subset of the BCAC data set (Oncoarray, iCOGS, or GWAS) or from the meta-analysis of all 3 together (BC meta-analysis, S1 Text). Results plotted are after pruning for instrument heterogeneity. Error bars represent the 95% CI. Estimates were calculated using the IVW method. *P < 0.05; **P < 0.001. See S10 Table for ORs, CIs, and P-values. AaM, age at menarche; BC, breast cancer; BCAC, Breast Cancer Association Consortium; BMI, body mass index; CI, confidence interval; GWAS, genome-wide association study; HDL, high-density lipoprotein; IVW, inverse-variance weighted; LDL, low-density lipoprotein; MR, Mendelian randomization; OR, odds ratio; TG, triglyceride.

https://doi.org/10.1371/journal.pmed.1003302.s008

(PDF)

S7 Fig. Multivariable MR analyses using lipid genetic associations from either MVP or GLGC.

Results of multivariable MR analyses including 2 lipid traits as exposures: (A) LDL and HDL or (B) TGs and HDL, with and without BMI as an additional exposure and with risk for all BCs as the outcome. Results plotted are after pruning for instrument heterogeneity. The lipid effect estimates were from one of 2 GWAS data sets (MVP or GLGC), and the results of each combination of lipid data sets are in a single plot. Error bars represent the 95% CI. Estimates were calculated using the IVW method. *P < 0.05; **P < 0.001. See S11 Table for ORs, CIs, and P-values. BC, breast cancer; BMI, body mass index; CI, confidence interval; GLGC, Global Lipids Genetics Consortium; GWAS, genome-wide association study; HDL, high-density lipoprotein; IVW, inverse-variance weighted; LDL, low-density lipoprotein; MR, Mendelian randomization; MVP, Million Veteran Program; OR, odds ratio; TG, triglyceride.

https://doi.org/10.1371/journal.pmed.1003302.s009

(PDF)

S8 Fig. Single-trait MR with BC outcomes stratified by ER subtypes.

Results of single-trait MR with each lipid trait as an exposure, and one of 3 BC traits as the outcome: all BC, ER− BCs only, or ER+ BCs only. Error bars represent the 95% CI. Estimates were calculated using a fixed-effects IVW method after pruning for instrument heterogeneity. Lipid association statistics come from the MVP data. **P < 0.001, *P < 0.05. See S12 Table for ORs, CIs, and P-values. BC, breast cancer; CI, confidence interval; ER, estrogen receptor; HDL, high-density lipoprotein; IVW, inverse-variance weighted; LDL, low-density lipoprotein; MR, Mendelian randomization; OR, odds ratio; TC, total cholesterol; TG, triglyceride.

https://doi.org/10.1371/journal.pmed.1003302.s010

(PDF)

S9 Fig. Multivariable MR analyses with BC outcomes stratified by ER subtypes.

Results of multivariable MR analyses with 3 lipid traits (HDL, LDL, TGs), BMI, and AaM as exposures; and all BCs, ER‒, or ER+ BC as the outcome. Results plotted are after pruning for instrument heterogeneity. Error bars represent the 95% CI. Estimates were calculated using the IVW method. *P < 0.05; **P < 0.001. See S13 Table for ORs, CIs, and P-values. AaM, age at menarche; BC, breast cancer; BMI, body mass index; CI, confidence interval; ER, estrogen receptor; HDL, high-density lipoprotein; IVW, inverse-variance weighted; LDL, low-density lipoprotein; MR, Mendelian randomization; TG, triglyceride.

https://doi.org/10.1371/journal.pmed.1003302.s011

(PDF)

S10 Fig. Genetic instruments’ effect estimates on HDL and BC at each canonical HDL metabolism pathway loci.

Conditionally independent HDL-associated SNPs at canonical HDL metabolism pathway genes, plotted by their conditional effect estimates on HDL (from MVP) and effect estimates on all BCs. Error bars represent 95% CIs. The dashed green line represents the regression line from fixed-effects IVW MR. BC, breast cancer; CI, confidence interval; HDL, high-density lipoprotein; IVW, inverse-variance weighted; MR, Mendelian randomization; MVP, Million Veteran Program.

https://doi.org/10.1371/journal.pmed.1003302.s012

(PDF)

S11 Fig. Genetic instruments’ effect estimates on LDL and BC at each canonical LDL metabolism pathway loci.

Conditionally independent LDL-associated SNPs at canonical LDL metabolism pathway genes, plotted by their conditional effect estimates on LDL (from MVP) and effect estimates on all BCs. Error bars represent 95% CIs. The dashed green line represents the regression line from fixed-effects IVW MR. BC, breast cancer; CI, confidence interval; IVW, inverse-variance weighted; LDL, low-density lipoprotein; MR, Mendelian randomization; MVP, Million Veteran Program.

https://doi.org/10.1371/journal.pmed.1003302.s013

(PDF)

S12 Fig. LDL locus-specific MR results with LDL as exposure and BC risk as outcome.

Forest plot of MR results for LDL gene-specific instruments (see S4 Fig) and meta-analysis of effect estimates across genes. Estimates were calculated using a fixed-effects IVW method. BC, breast cancer; CI, confidence interval; IVW, inverse-variance weighted; LDL, low-density lipoprotein; MR, Mendelian randomization; N SNPs, number of genetic instruments included in MR; OR, odds ratio.

https://doi.org/10.1371/journal.pmed.1003302.s014

(PDF)

S13 Fig. Genetic correlations between lipid and BC traits.

Results of LD-score regression testing for genetic correlation between each lipid trait and 3 BC traits: all BC, ER− BCs only, or ER+ BCs only. Error bars represent the 95% CI. Lipid association statistics were from a meta-analysis of GLGC and MVP. BC, breast cancer; CI, confidence interval; ER, estrogen receptor; GLGC, Global Lipids Genetics Consortium; HDL, high-density lipoprotein; LD, linkage disequilibrium; LDL, low-density lipoprotein; MVP, Million Veteran Program; TC; total cholesterol; TG, triglyceride.

https://doi.org/10.1371/journal.pmed.1003302.s015

(PDF)

S1 Table. Summary statistics for genetic instruments used in single-trait MR analyses.

Lipid exposure summary statistics are from the MVP European data set and BC summary statistics from the BCAC consortium meta-analysis. SNP: rsID of genetic instrument; exposure: lipid trait for exposure statistics; effect_allele.exposure: allele used for lipid and BC effect estimates; other_allele.exposure: noneffect allele; beta.exposure: effect size estimate for lipid trait; se.exposure: standard error of lipid effect size estimate; pval.exposure: P-value of lipid trait effect estimate; beta.bc: effect size estimate of BC risk; se.bc: standard error of BC risk effect estimate; pval.bc: P-value for BC risk effect estimate; inclPruned: logical, was SNP included in pruned single-trait MR analysis. BC, breast cancer; BCAC, Breast Cancer Association Consortium; HDL, high-density lipoprotein; LDL, low-density lipoprotein; MR, Mendelian randomization; MVP, Million Veteran Program; TC, total cholesterol; TG, triglyceride.

https://doi.org/10.1371/journal.pmed.1003302.s016

(XLSX)

S2 Table. Summary statistics for genetic instruments used in multivariable MR analyses.

SNP: rsID of genetic instrument; expZ: trait and data set used as first/second/… exposure (e.g., if exp1 is ldl_mvp, LDL summary statistics from GLGC were used as the first exposure); expZ_beta: effect size estimate for trait expZ; expZ_se: standard error of effect size estimate for trait expZ; expZ_pval: P-value of effect size estimate for trait expZ; bc_beta: BC risk effect size estimate; bc_se: standard error of BC effect size estimate; bc_pval: P-value of BC effect size estimate; test: unique identifier for each 2, 3, or 5-exposure MVMR experiment included in this table. AaM, age at menarche; BC, breast cancer; BMI, body mass index; GLGC, Global Lipids Genetics Consortium; HDL, high-density lipoprotein; LDL, low-density lipoprotein; MR, Mendelian randomization; MVMR, multivariable MR; MVP, Million Veteran Program; TC; total cholesterol; TG, triglyceride.

https://doi.org/10.1371/journal.pmed.1003302.s017

(XLSX)

S3 Table. Conditionally independent summary statistics for HDL or LDL associations used in locus-specific MR analyses.

Data are from the conditional analysis of summary statistics from MVP and GLGC meta-analysis published in Klarin and colleagues [22]. CHR: chromosome; POS: base position (HG19); SNP: rsID; effect.allele: allele used for effect size estimate; other.allele: noneffect allele; effect.allele.freq: frequency of effect allele; conditional.beta: effect size estimate from conditional analysis; conditional.se: standard error of conditional effect size estimate; conditional.p: P-value of conditional effect size estimate; Locus: the HDL or LDL gene or locus for MR analysis; Trait: lipid trait used as exposure (HDL or LDL). GLGC, Global Lipids Genetics Consortium; HDL, high-density lipoprotein; LDL, low-density lipoprotein; MR, Mendelian randomization; MVP, Million Veteran Program.

https://doi.org/10.1371/journal.pmed.1003302.s018

(XLSX)

S4 Table. Single-trait MR results with unpruned lipid traits as the exposure and all BCs as the outcome for a range of MR methods.

Lipid association statistics are from MVP. Exposure: lipid trait used as the exposure; Method: MR method used; N SNPs: number of genetic instruments included in analysis; CI_95_L: lower bound of 95% CI; CI_95_U: upper bound of 95% CI; P: P-value of MR test. BC, breast cancer; CI, confidence interval; HDL, high-density lipoprotein; LDL, low-density lipoprotein; MR, Mendelian randomization; MVP, Million Veteran Program; OR, odds ratio; TC, total cholesterol; TG, triglyceride.

https://doi.org/10.1371/journal.pmed.1003302.s019

(XLSX)

S5 Table. Heterogeneity analyses by Cochran's Q of unpruned single-trait MR.

Lipid association statistics are from MVP. Estimates are from the IVW method, and the outcome trait was risk for all BC. Exposure: lipid trait used as the exposure; Q: Cochran’s Q statistic; Q_df: degrees of freedom in Cochran’s Q test; Q_P: P-value of Cochran’s Q test. BC, breast cancer; HDL, high-density lipoprotein; IVW, inverse-variance weighted; LDL, low-density lipoprotein; MR, Mendelian randomization; MVP, Million Veteran Program; TC, total cholesterol; TG, triglyceride.

https://doi.org/10.1371/journal.pmed.1003302.s020

(XLSX)

S6 Table. Pleiotropy analysis using Egger regression of unpruned single-trait MR.

Lipid association statistics are from MVP, and the outcome trait was risk for all BC. Exposure: lipid trait used as the exposure; Egger intercept: intercept estimate from Egger regression; SE: standard error of intercept estimate; P: P-value of intercept estimate. BC, breast cancer; HDL, high-density lipoprotein; LDL, low-density lipoprotein; MR, Mendelian randomization; MVP, Million Veteran Program; TC, total cholesterol; TG, triglyceride.

https://doi.org/10.1371/journal.pmed.1003302.s021

(XLSX)

S7 Table. Results of single-trait MR, heterogeneity analyses, and directionality analyses with pruned lipid IV sets.

Lipid summary statistics were from MVP, and MR tests used the IVW method. Heterogeneity analyses used Cochran's Q, and directionality analyses used the Steiger test. Exposure: lipid trait used as the exposure; N_SNPs: number of genetic instruments included in analysis; OR: OR of MR test; CI_95_L: lower bound of 95% CI; CI_95_U: upper bound of 95% CI; MR_P: P-value of MR test; Q: Cochran’s Q statistic; Q_df: degrees of freedom in Cochran’s Q test; Q_Pval: P-value of Cochran’s Q test; SNP_r2_exposure: estimated variance in lipid trait explained by genetic instruments; SNP_r2_outcome: estimated variance in breast cancer risk explained by genetic instruments; Steiger_pval: P-value of Steiger test inference of causal direction; correct_causal_direction: logical, is the causal direction inferred by the Steiger test in the correct direction. CI, confidence interval; HDL, high-density lipoprotein; IV, instrumental variable; IVW, inverse-variance weighted; LDL, low-density lipoprotein; MR, Mendelian randomization; MVP, Million Veteran Program; OR, odds ratio; TC, total cholesterol; TG, triglyceride.

https://doi.org/10.1371/journal.pmed.1003302.s022

(XLSX)

S8 Table. Results of a reciprocal single-trait MR testing the effects of BC as the exposure on each lipid trait as the outcome.

Lipid summary statistics were from MVP, and MR tests used the IVW method. Exposure: BC used as exposure for all tests; Outcome: lipid trait used as the outcome; N_SNPs: number of genetic instruments used in MR test; CI_95_L: lower bound of 95% CI; CI_95_U: upper bound of 95% CI; P: P-value of MR test. BC, breast cancer; CI, confidence interval; HDL, high-density lipoprotein; IVW, inverse-variance weighted; LDL, low-density lipoprotein; MR, Mendelian randomization; MVP, Million Veteran Program; OR, odds ratio; TG, triglyceride.

https://doi.org/10.1371/journal.pmed.1003302.s023

(XLSX)

S9 Table. Results of single-trait MR analyses after pruning for genetic instruments associated with other lipid traits.

For each exposure, genetic instruments that were associated (P < 0.001) with the other 2 listed lipid traits were removed before this MR analysis. Lipid summary statistics were from MVP, MR tests used the IVW method, and risk for all BCs was the outcome. Exposure: lipid trait used as the exposure; N_SNPs: number of genetic instruments used in MR test; CI_95_L: lower bound of 95% CI; CI_95_U: upper bound of 95% CI; P: P-value of MR test. BC, breast cancer; CI, confidence interval; HDL, high-density lipoprotein; IVW, inverse-variance weighted; LDL, low-density lipoprotein; MR, Mendelian randomization; MVP, Million Veteran Program; OR, odds ratio; TG, triglyceride.

https://doi.org/10.1371/journal.pmed.1003302.s024

(XLSX)

S10 Table. Results of multivariable MR with 3 lipid traits, BMI, and AaM, as exposures and BC risk as outcome.

Results are from 4 separate multivariable MR experiments: 3 with outcome summary statistics from independent subsets of the BC data set (BC_Onco, Oncoarray; BC_iCoGS, iCOGS; or BC_GWAS, GWAS), or using the BCAC meta-analysis summary statistics (BC_Meta). Lipid summary statistics are from MVP. Before/after pruning: results from MR before or after pruning for instrument heterogeneity; Exposure: trait used as exposure; Outcome: BC summary statistics used for the outcome; N_SNPs: number of genetic instruments used in MR test; CI_95_L: lower bound of 95% CI; CI_95_U: upper bound of 95% CI; P: P-value of MR test. AaM, age at menarche; BC, breast cancer; BCAC, Breast Cancer Association Consortium; BMI, body mass index; CI, confidence interval; GWAS, genome-wise association study; HDL, high-density lipoprotein; LDL, low-density lipoprotein; MR, Mendelian randomization; MVP, Million Veteran Program; TG, triglyceride.

https://doi.org/10.1371/journal.pmed.1003302.s025

(XLSX)

S11 Table. Results of multivariable MR with lipid trait summary statistics from distinct cohorts.

Risk for all BCs was used as the outcome. Within a test, the effect estimates for each lipid trait are from different data sets (MVP or GLGC). Each test is separated by an empty row. Exposure_Dataset: lipid trait (or BMI) used as exposure and the data set for lipid summary statistics; N_SNPs: number of genetic instruments used in MR test; CI_95_L: lower bound of 95% CI; CI_95_U: upper bound of 95% CI; P: P-value of MR test; Before/after pruning: results from MR before or after pruning for instrument heterogeneity. BC, breast cancer; BMI, body mass index; CI, confidence interval; GLGC, Global Lipids Genetics Consortium; HDL, high-density lipoprotein; LDL, low-density lipoprotein; MR, Mendelian randomization; MVP, Million Veteran Program; OR, odds ratio; TG, triglyceride.

https://doi.org/10.1371/journal.pmed.1003302.s026

(XLSX)

S12 Table. Results of Cochran's Q test for heterogeneity on single-trait IVW MR results for the effect of a lipid trait on ER+ versus ER− BC.

Trait: lipid trait used as the exposure; ERposBeta: MR effect estimate for ER+ BCs; ERnegBeta: MR effect estimate for ER− BCs; ERposSE: standard error of MR effect estimate for ER+ BCs; ERnegSE: standard error of MR effect estimate for ER− BCs; Q: Cochran’s Q test statistic; P: P-value of Cochran’s Q test. BC, breast cancer; ER, estrogen receptor; HDL, high-density lipoprotein; IVW, inverse-variance weighted; LDL, low-density lipoprotein; MR, Mendelian randomization; TC, total cholesterol; TG, triglyceride.

https://doi.org/10.1371/journal.pmed.1003302.s027

(XLSX)

S13 Table. Results of multivariable MR with ER+ or ER− BC as the outcome.

For each test, 5 traits were used as the exposure traits and ER+ or ER− BC as the outcome. Results are after genetic instrument pruning. Lipid summary statistics are from MVP. N_SNPs: number of genetic instruments used in MR test; CI_95_L: lower bound of 95% CI; CI_95_U: upper bound of 95% CI; P: P-value of MR test. BC, breast cancer; BMI, body mass index; CI, confidence interval; ER, estrogen receptor; HDL, high-density lipoprotein; GLGC, Global Lipids Genetics Consortium; LDL, low-density lipoprotein; MR, Mendelian randomization; MVP, Million Veteran Program; OR, odds ratio; TG, triglycerides.

https://doi.org/10.1371/journal.pmed.1003302.s028

(XLSX)

S14 Table. Genome-wide genetic correlation estimates.

Genetic correlation estimates between each lipid trait and risk of all BCs, calculated by 2 methods. Lipid Trait: lipid trait used, Method: method used to estimate genetic correlation; Lipid intercept: estimate of lipid intercept from LDSC; Lipid intercept se: standard error of lipid intercept estimate from LDSC; Cross-trait intercept with BC-all: estimate of cross-trait intercept with risk of all BCs from LDSC; Cross-trait intercept se with BC-all: standard error of estimate of cross-trait intercept with risk of all BCs from LDSC; Genetic covariance: estimate of genetic covariance between lipid trait and BC; Covariance SE: standard error of genetic covariance estimate; Genetic correlation: estimate of genetic correlation between lipid trait and BC; Correlation SE: standard error of estimate of genetic correlation; Correlation Z-score: standard normalized estimate of genetic correlation; Correlation p-value: P-value of genetic correlation. BC, breast cancer; HDL, high-density lipoprotein; Hess, ρ-Hess method; LDL, low-density lipoprotein; LDSC, linkage disequilibrium-score regression; TC, total cholesterol; TG, triglycerides.

https://doi.org/10.1371/journal.pmed.1003302.s029

(XLSX)

S15 Table. Local genetic correlations calculated with the ρ-Hess method.

Genomic regions with significant local genetic correlation between BC and lipids using the ρ-Hess method. Only loci that passed Bonferroni correction (P = 0.05/1,703 partitions = 2.9 × 10−5) are shown. Lipid: lipid trait used; Position: genomic coordinates of loci tested (HG19); N_SNPs: number of SNPs in partition; K: number eigenvectors used; Local-rhog: local genetic correlation estimate; Variance: variance estimate; SE: standard error of genetic correlation estimate; Z: Z-score of genetic correlation estimate; P: P-value of genetic correlation estimate. BC, breast cancer.

https://doi.org/10.1371/journal.pmed.1003302.s030

(XLSX)

S1 Text. Supplementary Text.

Included are details about the GWASs utilized, heterogeneity analysis for single-trait MR, instrument strength and validity assessment for multivariable MR, and a list of investigators associated with the VA MVP (banner author). GWAS, genome-wide association study; MR, Mendelian randomization; MVP, Million Veteran Program; VA, Veterans Affairs.

https://doi.org/10.1371/journal.pmed.1003302.s031

(DOCX)

Citation: Johnson KE, Siewert KM, Klarin D, Damrauer SM, the VA Million Veteran Program, Chang K-M, et al. (2020) The relationship between circulating lipids and breast cancer risk: A Mendelian randomization study. PLoS Med 17(9): e1003302. https://doi.org/10.1371/journal.pmed.1003302

References

1. Torre LA, Siegel RL, Ward EM, Jemal A. Global Cancer Incidence and Mortality Rates and Trends—An Update. Cancer Epidemiol Biomarkers Prev. 2016;25: 16–27. pmid:26667886

2. Rose D, Gracheck P, Vona-Davis L, Rose DP, Gracheck PJ, Vona-Davis L. The Interactions of Obesity, Inflammation and Insulin Resistance in Breast Cancer. Cancers (Basel). 2015;7: 2147–2168. pmid:26516917

3. Kuzu OF, Noory MA, Robertson GP. The role of cholesterol in cancer. Cancer Research. 2016;76(8): 2063–2070. pmid:27197250

4. Touvier M, Fassier P, His M, Norat T, Chan DSM, Blacher J, et al. Cholesterol and breast cancer risk: a systematic review and meta-analysis of prospective studies. Br J Nutr. 2015;114: 347–357. pmid:26173770

5. Ni H, Liu H, Gao R. Serum Lipids and Breast Cancer Risk: A Meta-Analysis of Prospective Cohort Studies. Singh S, editor. PLoS ONE. 2015;10: e0142669. pmid:26554382

6. Martin LJ, Melnichouk O, Huszti E, Connelly PW, Greenberg CV., Minkin S, et al. Serum Lipids, Lipoproteins, and Risk of Breast Cancer: A Nested Case-Control Study Using Multiple Time Points. JNCI J Natl Cancer Inst. 2015;107: djv032–djv032. pmid:25817193

7. Zhong S, Zhang X, Chen L, Ma T, Tang J, Zhao J. Statin use and mortality in cancer patients: Systematic review and meta-analysis of observational studies. Cancer Treatment Reviews. 2015;41(6): 554–567. pmid:25890842

8. Borgquist S, Giobbie-Hurder A, Ahern TP, Garber JE, Colleoni M, Láng I, et al. Cholesterol, Cholesterol-Lowering Medication Use, and Breast Cancer Outcome in the BIG 1–98 Study. J Clin Oncol. 2017;35: 1179–1188. pmid:28380313

9. Orho-Melander M, Hindy G, Borgquist S, Schulz C-A, Manjer J, Melander O, et al. Blood lipid genetic scores, the HMGCR gene and cancer risk: a Mendelian randomization study. Int J Epidemiol. 2018;47(2): 495–505. Epub 2017 Nov 20. pmid:29165714

10. Nowak C, Ärnlöv J. A Mendelian randomization study of the effects of blood lipids on breast cancer risk. Nat Commun. 2018;9(1): 3957. pmid:30262900

11. Michailidou K, Lindström S, Dennis J, Beesley J, Hui S, Kar S, et al. Association analysis identifies 65 new breast cancer risk loci. Nature. 2017;551: 92–94. pmid:29059683

12. Global Lipids Genetics Consortium, Willer CJ, Schmidt EM, Sengupta S, Peloso GM, Gustafsson S, Kanoni S, et al. Discovery and Refinement of Loci Associated with Lipid Levels. Nat Genet. 2013;45: 1274–1283. pmid:24097068

13. Qi G, Chatterjee N. Mendelian randomization analysis using mixture models for robust and efficient estimation of causal effects. Nat Commun. 2019;10: 1941. pmid:31028273

14. Neuhouser ML, Aragaki AK, Prentice RL, Manson JE, Chlebowski R, Carty CL, et al. Overweight, Obesity, and Postmenopausal Invasive Breast Cancer Risk. JAMA Oncol. 2015;1: 611. pmid:26182172

15. Bulik-Sullivan B, Finucane HK, Anttila V, Gusev A, Day FR, Loh PR, et al. An atlas of genetic correlations across human diseases and traits. Nat Genet. 2015;47: 1236–1241. pmid:26414676

16. Lindström S, Finucane H, Bulik-Sullivan B, Schumacher FR, Amos CI, Hung RJ, et al. Quantifying the genetic correlation between multiple cancer types. Cancer Epidemiol Biomarkers Prev. 2017;26: 1427–1435. pmid:28637796

17. Day FR, Thompson DJ, Helgason H, Chasman DI, Finucane H, Sulem P, et al. Genomic analyses identify hundreds of variants associated with age at menarche and support a role for puberty timing in cancer risk. Nat Genet. 2017;49: 834–841. pmid:28436984

18. Burgess S, Thompson DJ, Rees JMB, Day FR, Perry JR, Ong KK. Dissecting causal pathways using mendelian randomization with summarized genetic data: Application to age at menarche and risk of breast cancer. Genetics. 2017;207: 481–487. pmid:28835472

19. Feng Y, Hong X, Wilker E, Li Z, Zhang W, Jin D, et al. Effects of age at menarche, reproductive years, and menopause on metabolic risk factors for cardiovascular diseases. Atherosclerosis. 2008;196: 590–597. pmid:17675039

20. Hamajima N, Hirose K, Tajima K, Rohan T, Friedenreich CM, Calle EE, et al. Menarche, menopause, and breast cancer risk: Individual participant meta-analysis, including 118 964 women with breast cancer from 117 epidemiological studies. Lancet Oncol. 2012;13: 1141–1151. pmid:23084519

21. Jiang X, Finucane HK, Schumacher FR, Schmit SL, Tyrer JP, Han Y, et al. Shared heritability and functional enrichment across six solid cancers. Nat Commun. 2019;10: 431. pmid:30683880

22. Klarin D, Damrauer SM, Cho K, Sun Y V., Teslovich TM, Honerlaw J, et al. Genetics of blood lipids among ~300,000 multi-ethnic participants of the Million Veteran Program. Nat Genet. 2018;50: 1514–1523. pmid:30275531

23. Yengo L, Sidorenko J, Kemper KE, Zheng Z, Wood AR, Weedon MN, et al. Meta-analysis of genome-wide association studies for height and body mass index in ~700000 individuals of European ancestry. Hum Mol Genet. 2018;27: 3641–3649. pmid:30124842

24. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MAR, Bender D, et al. PLINK: A Tool Set for Whole-Genome Association and Population-Based Linkage Analyses. Am J Hum Genet. 2007;81: 559–575. pmid:17701901

25. Hemani G, Zheng J, Elsworth B, Wade KH, Haberland V, Baird D, et al. The MR-Base platform supports systematic causal inference across the human phenome. Elife. 2018;7: e34408. pmid:29846171

26. Bowden J, Davey Smith G, Burgess S. Mendelian randomization with invalid instruments: effect estimation and bias detection through Egger regression. Int J Epidemiol. 2015;44: 512–525. pmid:26050253

27. Bowden J, Davey Smith G, Haycock PC, Burgess S. Consistent Estimation in Mendelian Randomization with Some Invalid Instruments Using a Weighted Median Estimator. Genet Epidemiol. 2016;40: 304–314. pmid:27061298

28. Hartwig FP, Smith GD, Bowden J. Robust inference in summary data Mendelian randomization via the zero modal pleiotropy assumption. Int J Epidemiol. 2017;46(6): 1985–1998. pmid:29040600

29. Walter K, Min JL, Huang J, Crooks L, Memari Y, McCarthy S, et al. The UK10K project identifies rare variants in health and disease. Nature. 2015;526: 82–90. pmid:26367797

30. Voight BF. MR-predictor: A simulation engine for Mendelian Randomization studies. Bioinformatics. 2014;30: 3432–3434. pmid:25165093

31. Burgess S. Sample size and power calculations in Mendelian randomization with a single instrumental variable and a binary outcome. Int J Epidemiol. 2014;43: 922–929. pmid:24608958

32. Burgess S, Thompson SG. Multivariable Mendelian randomization: The use of pleiotropic genetic variants to estimate causal effects. Am J Epidemiol. 2015;181: 251–260. pmid:25632051

33. Burgess S, Dudbridge F, Thompson SG. Re: “Multivariable Mendelian randomization: The use of pleiotropic genetic variants to estimate causal effects.” Am J Epidemiol. 2015;181: 290–291. pmid:25660081

34. Sanderson E, Davey Smith G, Windmeijer F, Bowden J. An examination of multivariable Mendelian randomization in the single-sample and two-sample summary data settings. Int J Epidemiol. 2019;48: 713–727. pmid:30535378

35. Schwarzer G. meta: An R package for meta-analysis. R News. 2007;7: 40–45.

36. Shi H, Mancuso N, Spendlove S, Pasaniuc B. Local Genetic Correlation Gives Insights into the Shared Genetic Architecture of Complex Traits. Am J Hum Genet. 2017;101: 737–751. pmid:29100087

37. Berisa T, Pickrell JK. Approximately independent linkage disequilibrium blocks in human populations. Bioinformatics. 2016;32: 283–5. pmid:26395773

38. The GTEx Consortium, Aguet F, Ardlie KG, Cummings BB, Gelfand ET, Getz G, et al. Genetic effects on gene expression across human tissues. Nature. 2017;550: 204–213. pmid:29022597

39. von Elm E, Altman DG, Egger M, Pocock SJ, Gøtzsche PC, Vandenbroucke JP. The Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) Statement: Guidelines for Reporting Observational Studies. PLoS Med. 2007;4: e296. pmid:17941714

40. Davey Smith G, Davies N, Dimou N, Egger M, Gallo V, Golub R, et al. STROBE-MR: Guidelines for strengthening the reporting of Mendelian randomization studies. PeerJ Preprints 27857 [Preprint]. 2019 [cited 2020 May 18]. https://peerj.com/preprints/27857/

41. Pickrell JK, Berisa T, Liu JZ, Ségurel L, Tung JY, Hinds DA. Detection and interpretation of shared genetic influences on 42 human traits. Nat Genet. 2016;48: 709–717. pmid:27182965

42. Hemani G, Tilling K, Davey Smith G. Orienting the causal relationship between imprecisely measured traits using GWAS summary data. Li J, editor. PLoS Genet. 2017;13: e1007081. pmid:29149188

43. Guo Y, Warren Andersen S, Shu XO, Michailidou K, Bolla MK, Wang Q, et al. Genetically Predicted Body Mass Index and Breast Cancer Risk: Mendelian Randomization Analyses of Data from 145,000 Women of European Descent. PLoS Med. 2016;13. pmid:27551723

44. Renehan AG, Tyson M, Egger M, Heller RF, Zwahlen M. Body-mass index and incidence of cancer: a systematic review and meta-analysis of prospective observational studies. Lancet. 2008;371: 569–78. pmid:18280327

45. Borgquist S, Butt T, Almgren P, Shiffman D, Stocks T, Orho-Melander M, et al. Apolipoproteins, lipids and risk of cancer. Int J Cancer. 2016;138: 2648–2656. pmid:26804063

46. Melvin JC, Seth D, Holmberg L, Garmo H, Hammar N, Jungner I, et al. Lipid profiles and risk of breast and ovarian cancer in the swedish AMORIS study. Cancer Epidemiol Biomarkers Prev. 2012;21: 1381–1384. pmid:22593241

47. Zhang B-L, He N, Huang Y-B, Song F-J, Chen K-X. ABO blood groups and risk of cancer: a systematic review and meta-analysis. Asian Pac J Cancer Prev. 2014;15: 4643–50. pmid:24969898

48. Beeghly-Fadiel A, Khankari NK, Delahanty RJ, Shu X-O, Lu Y, Schmidt MK, et al. A Mendelian randomization analysis of circulating lipid traits and breast cancer risk. Int J Epidemiol. Epub 2019 Dec 23. pmid:31872213

49. Madsen CM, Varbo A, Nordestgaard BG. Extreme high high-density lipoprotein cholesterol is paradoxically associated with high mortality in men and women: two prospective cohort studies. Eur Heart J. 2017;38: 2478–2486. pmid:28419274

50. Hamer M, O’Donovan G, Stamatakis E. High-Density Lipoprotein Cholesterol and Mortality: Too Much of a Good Thing? Arterioscler Thromb Vasc Biol. 2018;38: 669–672. pmid:29326314

51. Burgess S, Davey Smith G. Mendelian Randomization Implicates High-Density Lipoprotein Cholesterol–Associated Mechanisms in Etiology of Age-Related Macular Degeneration. Ophthalmology. 2017;124: 1165–1174. pmid:28456421

52. Fan Q, Maranville JC, Fritsche L, Sim X, Cheung CMG, Chen LJ, et al. HDL-cholesterol levels and risk of age-related macular degeneration: a multiethnic genetic study using Mendelian randomization. Int J Epidemiol. 2017;46: 1891–1902. pmid:29025108

53. Bonovas S, Filioussi K, Tsavaris N, Sitaras NM. Use of Statins and Breast Cancer: A Meta-Analysis of Seven Randomized Clinical Trials and Nine Observational Studies. J Clin Oncol. 2005;23: 8606–8612. pmid:16260694

54. Llewellyn-Bennett R, Edwards D, Roberts N, Hainsworth AH, Bulbulia R, Bowman L. Post-trial follow-up methodology in large randomised controlled trials: a systematic review. Trials. 2018;19: 298. pmid:29843774

55. VanderWeele TJ, Tchetgen Tchetgen EJ, Cornelis M, Kraft P. Methodological challenges in mendelian randomization. Epidemiology. 2014;25: 427–35. pmid:24681576

56. O’Connor LJ, Price AL. Distinguishing genetic correlation from causation across 52 diseases and complex traits. Nat Genet. 2018;50: 1728–1734. pmid:30374074

57. Haycock PC, Burgess S, Wade KH, Bowden J, Relton C, Davey Smith G. Best (but oft-forgotten) practices: the design, analysis, and interpretation of Mendelian randomization studies. Am J Clin Nutr. 2016;103: 965–978. pmid:26961927

58. Silvente-Poirot S, Poirot M. Cancer. Cholesterol and cancer, in the balance. Science. 2014;343(6178): 1445–1446. pmid:24675946

59. Nelson ER. The significance of cholesterol and its metabolite, 27-hydroxycholesterol in breast cancer. Mol Cell Endocrinol. 2018;466: 73–80. pmid:28919300