Bidirectional causal relationship between hypercholesterolemia and ischemic heart disease: a Mendelian randomization study
Ying Jiang
Wenpeng Yu1,† 
Jianliang Zhou
Xiao Dong- 1Department of Cardiovascular Surgery, The Second Affiliated Hospital of Nanchang University, Nanchang, China
- 2Department of Cardiovascular Surgery, Zhongnan Hospital of Wuhan University, Wuhan, China
Background: Ischemic Heart Disease (IHD) is a leading cause of morbidity and mortality worldwide. Hypercholesterolaemia, a metabolic syndrome distinguished by elevated cholesterol levels, is positively correlated with IHD, yet the precise causal relationship between these two health conditions remains to be clearly defined.
Methods: We conducted a two-sample MR analysis using genetic variants associated with hypercholesterolemia and IHD. Various statistical techniques including MR-Egger, Weighted Median, Inverse Variance Weighted (IVW), Simple Mode, and Weighted Mode were employed. We also performed sensitivity analyses to assess pleiotropy, heterogeneity, and influence of individual SNPs. Furthermore, genetic co-localization analysis was performed to identify shared genes between hypercholesterolemia and IHD.
Results: Our MR study illuminated a bidirectional causal relationship between hypercholesterolaemia and ischaemic heart disease. Utilising the IVW with multiplicative random effects, upon considering IHD as the outcome, we identified an OR of 2.27 (95% CI: 1.91–2.70, p = 1.68 × 10−20). Conversely, when hypercholesterolaemia was viewed as the outcome, the OR detected was 1.80 (95% CI: 1.58–2.05, p = 2.79 × 10−19). These findings remained consistent across various MR methods and sensitivity analyses. Additionally, our research pinpointed four co-localised genes CELSR2, PCSK9, LPA, and APOE as integral candidates implicated in the pathogenesis of both conditions, thereby suggesting shared common genetic causal variants and offering potential targets for innovative therapeutic strategies.
Conclusion: bidirectional MR studies reveal genetic evidence of a potential causal link between hypercholesterolaemia and IHD. Notably, these findings also lend credence to the less traditional hypothesis that IHD may instigate hypercholesterolaemia episodes. Moreover, co-localisation analyses intimate the presence of shared genetic causal variants, paving the way for the development of new therapeutic strategies.
1. Introduction
Ischemic Heart Disease (IHD) is a major cause of global mortality, claiming over 9 million lives annually. With approximately 126.5 million cases globally, its impact on public health is significant. In the United States, it presents both a substantial health challenge and an economic encumbrance, with yearly costs estimated in the hundreds of billions (1–3). IHD is principally defined by diminished myocardial perfusion, instigated by a complex interplay of pathophysiological factors within the coronary artery wall's intimal layer. These elements include endothelial dysfunction, inflammation, thrombogenesis, and the aftereffects of angiogenesis and calcification. A notable contributor to this pathological nexus is the presence of oxidized serum lipids, a hallmark of hypercholesterolemia (4–6). Elevated LDL cholesterol levels lead to hypercholesterolemia, endothelial dysfunction, and increased atherosclerosis, potentially hastening death in patients with severe ischemic heart disease (7–9).
The chief contributors to hypercholesterolemia include an excessive biosynthesis of cholesterol and/or inadequate clearance of it. It's noteworthy that hypercholesterolemia amplifies the risk for developing atherosclerosis and ischemic heart disease (10). Although hypercholesterolemia's role as a significant risk factor for IHD is widely recognized (11), the suggestion of a causal relationship flowing in the opposite direction—that is, IHD potentially instigating the onset of hypercholesterolemia—presents a less conventional perspective. This concept implies a bidirectional causality between hypercholesterolemia and IHD, a process that warrants meticulous exploration. In our quest to decipher the causal link between the two, we employed an approach known as Mendelian randomization (MR). MR leverages genetic variants as instrumental variables tied to an exposure, consequently amplifying the dependability of the causal inference connecting an exposure and its subsequent outcome. This method helps to address potential challenges such as confounding factors and reverse causality (12, 13). If there exists a causal relationship between the exposure and the outcome, genetic variations influencing the exposure will likely affect the outcome as well. MR is akin to a natural randomized controlled trial, being based on Mendel's second law. This law proposes that alleles from distinct genes are assorted independently during gametogenesis. Consequently, the underlying premise of Mendelian randomization analyses is that the inheritance of a particular trait should remain independent of the inheritance of any other traits (14).
Hence, in our study, we employed a bidirectional MR analysis utilizing genetic variants linked to hypercholesterolemia and IHD. Our statistical methodology encompassed MR-Egger, Weighted Median, IVW, Simple Mode, and Weighted Mode. Additionally, we conducted sensitivity analyses to evaluate the influences of pleiotropy, heterogeneity, and individual single nucleotide polymorphisms (SNPs). Ultimately, we carried out a genetic co-localization analysis to identify shared genes between hypercholesterolemia and IHD. Gaining an understanding of the bidirectional relationship between hypercholesterolemia and IHD holds instructive value in terms of risk monitoring for individuals with either disease and identifying colocalized genes as potential targets for therapeutic interventions.
2. Materials and methods
2.1. Study design
SNPs exemplifying genetic diversity were chosen as instrumental variables and underwent a two-sample MR analysis. Three core assumptions (15) were established as follows Figure 1: (1) a direct connection exists between instrumental variables and exposure; (2) instrumental variables remain unconnected to any potential confounders; (3) the impact of genetic variations on the results occurs solely through exposure. We use Mendelian randomization (MR) analysis to assess the two-way causal link between IHD and hypercholesterolemia. To further investigate the shared genetic basis between IHD and hypercholesterolemia, genetic colocalisation analysis was performed. This method allows us to assess whether the genetic signals identified in two independent association studies, such as those for IHD and hypercholesterolemia, are colocalized, indicating a common genetic mechanism underlying both traits. By integrating the results from the MR analysis with the colocalization analysis, we can gain a deeper understanding of the potential shared genetic influences and biological pathways involved in the development of IHD and hypercholesterolemia.
Figure 1. Three key assumptions underlying Mendelian randomization study design. The three different hypotheses are represented by three pathways. Hypothesis 1: SNP is associated with exposure. Hypothesis 2: SNP affects outcome only through exposure and not through any alternative causal pathways. Hypothesis 3: SNPs are completely independent of any potential confounding factors that influence exposure and outcome.
2.2. Genome-wide association analysis
Genome-wide association study (GWAS) databases (16), such as GWAS Catalog, IEU openGWAS, and NealELab, were explored and suitable datasets were obtained. Given that all data utilized was publicly available, no further ethical approval was necessary. To minimize the bias stemming from ethnicity-related confounding elements, the study's genetic background was limited to individuals of European descent. The IEU GWAS report encompasses two distinct traits: “Pure hypercholesterolemia” and “Ischaemic heart disease”. The first study, identified as ukb-b-12651, was conducted in 2018 by Ben Elsworth, while the second study, identified as ukb-d-I9_IHD, was carried out in the same year by Neale Lab. Both studies used phesant-derived variables from the UK Biobank and included participants of European descent, both males and females. The ukb-b-12651 study had a sample size of 463,010 participants, consisting of 22,622 cases and 440,388 controls. On the other hand, the ukb-d-I9_IHD study included a sample size of 361,194 participants, with 20,857 cases and 340,337 controls. Both studies employed the HG19/GRCh37 genome build and focused on binary GWAS categories, incorporating MR analysis. These reports are publicly accessible and were conducted as part of the MRC-IEU consortium.
2.3. SNP selection and statistical analysis
We adopted a two-sample MR approach, facilitated by the TwoSampleMR R package, to explore the two-way causal relationship between hypercholesterolemia and ischemic heart disease (16). MR analysis is performed by following a careful workflow to ensure the credibility and robustness of the results.
2.3.1. Instrument selection
We first identified SNPs from the exposure dataset (hypercholesterolemia) and the outcome dataset (ischemic heart disease). SNPs were selected as IVs based on their genome-wide significance (p-value < 5 × 10 × 10−8) and linkage disequilibrium clumping with an r2 threshold of 0.001 and a 10,000 kb window to ensure their independence.
2.3.2. Data harmonization
The summary-level data from GWAS for the selected SNPs were extracted, and the exposure and outcome datasets were harmonized to ensure the consistency of effect allele coding and alignment of SNP effect estimates (17). Any strand-ambiguous or palindromic SNPs were appropriately handled or excluded from the analysis.
2.3.3. F-statistic calculation
We set a threshold of an F-statistic equal to or greater than 10 for conducting the MR analysis to reduce the potential impact of weak instrument bias (18). The F-statistic can be calculated using the formula: , where N is the sample size, k represents the number of instrumental variables (IVs), and R2 indicates the proportion of variation explained by the SNPs. The calculation of R2 for the 5 genome-wide significant SNP instrument was performed using the formula: . In contrast, for the extended 10 SNP instrument, the calculation of R2 was conducted using the equation: (19). The calculation of R2 was determined based on the number of SNPs involved in the analysis. The equation includes variables such as EAF, beta, N, and SE (beta). EAF represents the effect allele frequency, beta is the estimated genetic effect on physical activity, N is the sample size of the GWAS for the SNP-physical activity association, and SE (beta) signifies the standard error of the genetic effect. By adhering to this threshold and utilizing the F-statistic, we can ensure that our MR analysis accounts for the strength of the genetic instruments and minimizes the risk of weak instrument bias, thus contributing to the validity of the causal inferences made from the analysis.
2.3.4. MR analysis
We conducted the MR analysis using the harmonized data to estimate the causal effect of hypercholesterolemia on ischemic heart disease. Several MR methods, including the IVW, MR-Egger, and Weighted Median, were used to obtain robust causal effect estimates and account for potential biases in the analysis.
2.3.5. Heterogeneity, pleiotropy, and sensitivity assessment
The Cochran's Q statistic was employed to assess the heterogeneity of individual SNP estimates. To assess the impact of individual SNPs on the overall causal estimate, a leave-one-out analysis was conducted. The MR-Egger regression intercept was examined for evidence of pleiotropy, and the funnel plot symmetry was visually inspected.
2.3.6. Reverse direction analysis
We then reversed the exposure and outcome datasets by defining the exposure dataset for ischemic heart disease and the outcome dataset for hypercholesterolemia. Following the same procedure as described above, the MR analysis was performed to investigate the causal effect of ischemic heart disease on hypercholesterolemia.
2.3.7. Result visualization
Scatter plots, forest plots, funnel plots, and leave-one-out plots were generated to visualize the results of both MR analyses, including the causal estimates, heterogeneity, pleiotropy, and sensitivity analyses. These visualizations allowed for a comprehensive assessment of the robustness and validity of the MR findings.
2.4. Colocalization analysis of Complex traits
We utilized the R package “ieugwasr” to obtain top hits from the UK Biobank dataset for two traits of interest (traits “ukb-b-12651” and “ukb-d-I9 IHD”). The top hits were sorted by their p-values. To perform the colocalization analysis, we first determined the best signals for each trait and created a range around these loci. We considered a region of ±90,000 base pairs (bp) for most of the top hits, with some exceptions where we adjusted the range based on the specific chromosomal position. The R package “coloc” was then employed to conduct a colocalization analysis between the two traits (20), using the “coloc.abf” function applied to the extracted and formatted data obtained from the “ieugwasr_to_coloc” function. After completing the colocalization analysis, we visualized the results using the “gassocplot” R package. The “stack_assoc_plot” function was applied to the output from the “coloc_to_gassocplot” function, generating a stacked Manhattan plot representing the genetic association signals from both traits.
3. Results
3.1. MR analysis of bidirectional associations between pure hypercholesterolemia and ischemic heart disease
In this study, we conducted a MR analysis to investigate the causal relationship between hypercholesterolemia and ischemic heart disease. We employed five MR methods, including MR Egger, Weighted Median, IVW, Simple Mode, and Weighted Mode, to obtain robust estimates of causal effects. All SNPs utilized in this study exhibited F statistics exceeding 10, indicating robust instrumental variables.
3.1.1. Causal effect of hypercholesterolemia on IHD via forward MR
The main results for ischemic heart disease as the outcome are as follows: MR Egger: The odds ratio (OR) for the causal effect of hypercholesterolemia on ischemic heart disease was 2.95 (95% CI: 2.11–4.14, p = 1.09 × 10−5). Weighted Median: The OR for the causal effect was 1.87 (95% CI: 1.66–2.10, p = 6.25 × 10−26). Inverse Variance Weighted (IVW): The OR for the causal effect was 2.27 (95% CI: 1.91–2.70, p = 1.68 × 10−20). Simple Mode: The OR for the causal effect was 1.97 (95% CI: 1.58–2.47, p = 1.63 × 10−5). Weighted Mode: The OR for the causal effect was 1.91 (95% CI: 1.67–2.19, p = 3.87 × 10−8). We conducted several sensitivity analyses to evaluate the assumptions, pleiotropy, heterogeneity, and the influence of individual SNPs on our MR analysis. Pleiotropy: The MR-Egger regression test did not indicate significant pleiotropy (egger intercept = −1.66 × 10−3, SE = 9.49 × 10−4, p = 0.0995), suggesting that the genetic instruments are not biased due to horizontal pleiotropy. Heterogeneity: The Cochran's Q test for both MR-Egger (Q = 127.12, df = 16, p = 2.32 × 10−19) and Inverse Variance Weighted (IVW) (Q = 151.41, df = 17, p = 1.29 × 10−23) methods revealed significant heterogeneity across the genetic instruments. Due to the significant heterogeneity revealed by the Cochran's Q test in our MR-Egger and Inverse Variance Weighted (IVW) analyses, we opted to use the Inverse Variance Weighted method with multiplicative random effects in our study. Leave-one-out analysis: The analysis demonstrated that no individual SNP had a disproportionate influence on the overall causal estimates. The results remained consistent after excluding each SNP, suggesting the robustness of the MR findings. All of the above results can be found in Supplementary file 1.
To visualize these results, we utilized forest plots (Figure 2A) and scatter plots (Figure 2B). While other visualization plots used to assess heterogeneity, polymorphism and sensitivity analysis include forest plots, funnel plots, and leave-one-out plots, which we can find in Supplementary Figure S1.
Figure 2. Visualization of hypercholesterolemia as an exposure factor in a two-sample MR. (A) Forest plots depicting the potential impact of SNPs associated with hypercholesterolaemia on IHD using five MR methods. (B) Scatter Plots Illustrating the Effect of SNPs on Hypercholesterolaemia and its Correlation with the Effect on IHD. Circles indicate marginal genetic associations with hypercholesterolemia and risk of IHD for each variant. Error bars indicate 95% CIs. IHD, ischemic heart disease; MR, mendelian randomization; SNP, single nucleotide polymorphism.
3.1.2. Causal association of IHD with hypercholesterolemia via reverse MR
The main results for hypercholesterolemia as the outcome are as follows: MR Egger: The odds ratio (OR) for the causal effect of hypercholesterolemia on ischemic heart disease was 1.88 (95% CI: 1.43–2.48, p = 1.27 × 10−4). Weighted Median: The OR for the causal effect was 1.46 (95% CI: 1.34–1.58, p = 2.48 × 10−20). Inverse Variance Weighted (IVW): The OR for the causal effect was 1.80 (95% CI: 1.58–2.05, p = 2.79 × 10−19). Simple Mode: The OR for the causal effect was 1.44 (95% CI: 1.28–1.61, p = 1.01 × 10−6). Weighted Mode: The OR for the causal effect was 1.42 (95% CI: 1.30–1.55, p = 1.35 × 10−8).
Pleiotropy: The results of the MR-Egger regression analysis revealed no significant pleiotropy (egger intercept = −3.00 × 10−4, SE = 8.29 × 10−4, p = 0.720). This finding suggests that the genetic instruments used in our study are not biased due to horizontal pleiotropy, and the observed associations between hypercholesterolemia and ischemic heart disease are likely to reflect a causal relationship.
Heterogeneity: The Cochran's Q test for both MR-Egger (Q = 343.34, df = 27, p = 1.50 × 10−56) and Inverse Variance Weighted (IVW) (Q = 345.01, df = 28, p = 2.51 × 10−56) methods revealed significant heterogeneity across the genetic instruments. Due to the significant heterogeneity revealed by the Cochran's Q test in our MR-Egger and Inverse Variance Weighted (IVW) analyses, we opted to use the Inverse Variance Weighted method with multiplicative random effects in our study. This approach, which accounts for heterogeneity among the genetic instruments, provides a more accurate and robust estimation of the causal effect. Leave-one-out analysis: The analysis demonstrated that no individual SNP had a disproportionate influence on the overall causal estimates for the association between Ischaemic heart disease (wide definition) and pure hypercholesterolaemia. The results remained consistent after excluding each SNP, suggesting the robustness of the MR findings. All of the above results can be found in Supplementary file 2.
To visualize these results, we utilized forest plots (Figure 3A) and scatter plots (Figure 3B). While other visualization plots used to assess heterogeneity, polymorphism and sensitivity analysis include forest plots, funnel plots, and leave-one-out plots, which we can find in Supplementary Figure S2.
Figure 3. Visualization of IHD as an exposure factor in a two-sample MR. (A) Forest plots depicting the potential impact of SNPs associated with IHD on Hypercholesterolaemia using five MR methods. (B) Scatter Plots Illustrating the Effect of SNPs on IHD and its Correlation with the Effect on Hypercholesterolaemia. Circles indicate marginal genetic associations with IHD and risk of Hypercholesterolaemia for each variant. Error bars indicate 95% CIs. IHD, ischemic heart disease; MR, mendelian randomization; SNP, single nucleotide polymorphism.
3.2. Genetic co-localization analysis
Our analysis identified four colocalized genes, CELSR2, PCSK9, LPA, and APOE, that are involved in lipid metabolism and may contribute to the development of both Pure hypercholesterolaemia and IHD. CELSR2 (Figure 4A) and PCSK9 (Figure 4B) are both located on chromosome 1 and have distinct SNPs associated with each condition. LPA (Figure 4C), located on chromosome 6, and APOE (Figure 4D), located on chromosome 19, also have unique SNPs associated with each trait. CELSR2, which encodes a transmembrane protein involved in cell adhesion and lipid metabolism regulation, has SNPs rs646776 for Pure hypercholesterolaemia and rs660240 for IHD. PCSK9, encoding a proprotein convertase that regulates LDL receptor expression and is targeted by cholesterol-lowering drugs, shares the SNP rs11591147 for both conditions. LPA, encoding lipoprotein(a) and involved in cholesterol metabolism, has SNPs rs55730499 for Pure hypercholesterolaemia and rs10455872 for IHD. Lastly, APOE, encoding apolipoprotein E and playing a key role in lipid metabolism, has SNPs rs7412 for Pure hypercholesterolaemia and rs429358 for IHD. These findings provide valuable insights into the shared genetic basis of Pure hypercholesterolaemia and IHD, potentially informing the development of new therapeutic strategies.
Figure 4. Stacked regional association plot. (A) CELSR2, located on chromosome 1, encodes a transmembrane protein involved in cell adhesion and lipid metabolism regulation, with SNP rs646776 associated with Pure hypercholesterolaemia and SNP rs660240 associated with IHD. (B) PCSK9, also located on chromosome 1, encodes a proprotein convertase that regulates LDL receptor expression. It is associated with SNP rs11591147 for both Pure hypercholesterolaemia and IHD. (C) LPA, located on chromosome 6, encodes lipoprotein(a) involved in cholesterol metabolism and is associated with SNP rs55730499 for Pure hypercholesterolaemia and SNP rs10455872 for IHD. (D) APOE, located on chromosome 19, encodes apolipoprotein E, a key player in lipid metabolism. It is associated with SNP rs7412 for Pure hypercholesterolaemia and SNP rs429358 for IHD.
4. Discussion
To shed light on the causal relationship between pure hypercholesterolemia and ischemic heart disease, we conducted bidirectional Mendelian randomization (MR) using publicly available GWAS data in both groups. Our findings reveal a bidirectional causal relationship between genetically confirmed pure hypercholesterolemia and ischemic heart disease in a European population.
First, forward MR analysis demonstrated a significant causal effect of hypercholesterolemia on IHD. The odds ratio (OR) for all five MR methods were greater than 1, signifying that hypercholesterolemia increases the risk of IHD. This finding aligns with the existing literature stating that hypercholesterolemia is a significant risk factor for cardiovascular disease (11, 21, 22). For instance, numerous studies (23–25) have shown that mononucleosis and neutrophilia, induced by hypercholesterolemia, contribute to the progression of atherosclerosis, which can subsequently lead to IHD. Additionally, Adams et al. (26) proposed in their review that oxidized LDL might trigger endothelial cell dysfunction, thereby leading to the development of ischemic syndromes. They emphasized that cellular activity could be improved within hours to days following cholesterol lowering. Importantly, another study (27) highlighted a nonlinear relationship between total cholesterol and cardiovascular disease (CVD), including IHD, with escalating mortality in IHD when total cholesterol levels exceeded 200 mg/dl. After conducting a thorough analysis, a significant association was discovered between elevated total cholesterol levels and an augmented risk of CVD mortality, specifically IHD. Our reverse MR analysis also demonstrated a significant causal effect of IHD on hypercholesterolemia. OR for all five MR methods were greater than 1, indicating that IHD increases the risk of hypercholesterolemia. This could be attributable to potential alterations in lipid metabolism resulting from the disease condition. One study (28) identified a clear association between myocardial ischemia and significant shifts in lipid metabolism. These changes were characterized by alterations in lipoprotein subclasses and an increase in total HDL cholesterol levels, especially prominent in the early post-ischemic phase.
We employed genetic co-localization analysis for hypercholesterolemia and IHD to determine whether these two phenotypes share causal genetic variants within a given region. We identified four co-localized genes CELSR2, PCSK9, LPA, and APOE. The recognition of these genes bolsters our understanding of the genetic correlation between hypercholesterolemia and IHD. Each of these genes plays a crucial role in lipid metabolism (29–32), which is a central player in the pathogenesis of both conditions. For instance, a study (29) demonstrates that CELSR2 depletion leads to a considerable reduction in lipid accumulation within hepatocytes. Additionally, it reveals that a deficiency in CELSR2 appears to undermine cell survival by inhibiting cell proliferation and encouraging apoptosis. As such, CELSR2 emerges as a potential therapeutic target for lipid metabolism regulation. PCSK9, a secreted protein primarily produced by the liver, plays a crucial role in the regulation of lipid metabolism by interacting with LDLR and other receptors, thereby contributing to cellular lipid accumulation. Its deficiency has been observed to affect heart metabolism and function, potentially leading to heart failure with preserved ejection fraction (HFpEF), indicating a strong link between PCSK9 and lipid metabolism in the cardiovascular system (33). LPA and its associated oxidized phospholipids are known to contribute significantly to chronic inflammation, a crucial factor in the development and progression of cardiovascular diseases including atherosclerosis. Further, oxidized lipids from LPA interact with immune and endothelial cells, prompting inflammatory responses that underlie various cardiovascular conditions such as myocardial infarction and calcific aortic valve stenosis, highlighting the critical role of LPA in lipid metabolism and cardiovascular health (34). Apolipoprotein E (apoE), which is predominantly synthesized in the liver, plays an integral role in lipid metabolism. This role includes mediating the clearance of triglyceride-rich lipoproteins and their remnants, initiating the reverse transport of cholesterol to the liver, and distributing lipids among cells in the nervous system. Furthermore, the various isoforms of apoE each have different impacts on the concentration of low-density lipoprotein and the risk of atherosclerosis, thereby indicating their significant influence on cardiovascular health (35). In light of these findings, the co-localised genes we identified indeed play pivotal roles in lipid metabolism, metabolic processes, and inflammatory responses, demonstrating considerable potential as therapeutic targets in hypercholesterolaemia and IHD. Specifically, PCSK9 inhibitors have emerged as a novel category of lipid-lowering agents. According to the 2019 guidelines, PCSK9 inhibitors are recommended for patients who have not attained their target lipid profile following ezetimibe/statin therapy (36). The primary mechanism of PCSK9 inhibitors involves enhancing LDL receptor density by inhibiting PCSK9 protein activity, thereby facilitating a reduction in LDL cholesterol levels (37–39). Demonstrating effectiveness in studies addressing myocardial infarction (40) and hypercholesterolaemia (39), drugs that target PCSK9 hold substantial promise for treating atherosclerosis, hypercholesterolaemia, and other related cardiovascular diseases.
In conclusion, this study provides strong evidence of causality between hypercholesterolemia and IHD and reverse causality, revealing the interaction between these two diseases that have a major impact on global health. Employing Mendelian randomization, we leveraged large-scale GWAS data on risk factors to mitigate bias from residual confounding and reverse causality, thereby minimizing the effects of population stratification (41, 42). Moreover, only cohorts predominantly of European ancestry were allowed to provide data. In accordance with the three fundamental assumptions of the MR study, our findings were robust and corroborated by various MR methods, including IVW, weighted median, and MR Egger. We also conducted an additional gene co-localisation analysis, identifying four genes, namely CELSR2, PCSK9, LPA, and APOE, co-localised with hypercholesterolaemia and IHD. Understanding these genes’ functions could aid in developing targeted therapeutic strategies for these diseases. However, our study did encounter several limitations. For example, there was significant heterogeneity across the genetic instruments in both the forward and reverse MR analyses. Although we employed the Inverse Variance Weighted method with multiplicative random effects to account for this heterogeneity, it's imperative to approach our results with caution. Another limitation stems from the genetic variants used as instruments in MR studies. Sometimes these variants can influence multiple traits or be associated with confounding factors, which can potentially bias the results. Furthermore, our findings are largely applicable to populations of European ancestry, limiting the generalizability to other ethnic groups.
5. Conclusions
In summary, our bidirectional MR results deliver genetic evidence, endorsing a potential causal linkage between hypercholesterolaemia and IHD. These findings also corroborate the less conventional hypothesis that IHD may instigate episodes of hypercholesterolaemia. The co-localisation analyses further intimate the presence of shared genetic causal variants, thereby presenting promising targets for the formulation of innovative therapeutic strategies.
Data availability statement
Publicly available datasets were analyzed in this study. This data can be found here: https://gwas.mrcieu.ac.uk/datasets/, using codes ukb-b-12651 and ukb-d-I9_IHD.
Author contributions
YJ: Writing – original draft. WY: Writing – original draft. JZ: Writing – review & editing. XD: Writing – review & editing.
Funding
The author(s) declare that no financial support was received for the research, authorship, and/or publication of this article.
Acknowledgments
We would like to extend our sincere appreciation to the participants and investigators whose contributions made the availability of summary data possible.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcvm.2023.1302282/full#supplementary-material
References
1. Benjamin EJ, Virani SS, Callaway CW, Chamberlain AM, Chang AR, Cheng S, et al. Heart disease and stroke statistics-2018 update: a report from the American heart association. Circulation. (2018) 137(12):e67–492.29386200
2. Nowbar AN, Gitto M, Howard JP, Francis DP, Al-Lamee R. Mortality from ischemic heart disease. Circ Cardiovasc Qual Outcomes. (2019) 12(6):e005375. doi: 10.1161/CIRCOUTCOMES.118.005375
3. Virani SS, Alonso A, Benjamin EJ, Bittencourt MS, Callaway CW, Carson AP, et al. Heart disease and stroke statistics-2020 update: a report from the American heart association. Circulation. (2020) 141(9):e139–596.31992061
4. Libby P, Theroux P. Pathophysiology of coronary artery disease. Circulation. (2005) 111(25):3481–8. doi: 10.1161/CIRCULATIONAHA.105.537878
5. Ambrose JA, Singh M. Pathophysiology of coronary artery disease leading to acute coronary syndromes. F1000Prime Rep. (2015) 7:8. doi: 10.12703/P7-08
6. Duarte MM, Rocha JB, Moresco RN, Duarte T, Da Cruz IB, Loro VL, et al. Association between ischemia-modified albumin, lipids and inflammation biomarkers in patients with hypercholesterolemia. Clin Biochem. (2009) 42(7–8):666–71. doi: 10.1016/j.clinbiochem.2009.01.010
7. Domanski MJ, Tian X, Wu CO, Reis JP, Dey AK, Gu Y, et al. Time course of LDL cholesterol exposure and cardiovascular disease event risk. J Am Coll Cardiol. (2020) 76(13):1507–16. doi: 10.1016/j.jacc.2020.07.059
8. Yuan S, Huang X, Ma W, Yang R, Xu F, Han D, et al. Associations of HDL-C/LDL-C with myocardial infarction, all-cause mortality, haemorrhagic stroke and ischaemic stroke: a longitudinal study based on 384 093 participants from the UK biobank. Stroke Vasc Neurol. (2023) 8(2):119–26. doi: 10.1136/svn-2022-001668
9. Liao J, Guo X, Wang M, Dong C, Gao M, Wang H, et al. Scavenger receptor class B type 1 deletion led to coronary atherosclerosis and ischemic heart disease in low-density lipoprotein receptor knockout mice on modified western-type diet. J Atheroscler Thromb. (2017) 24(2):133–46. doi: 10.5551/jat.33019
10. Goldstein JL, Brown MS. A century of cholesterol and coronaries: from plaques to genes to statins. Cell. (2015) 161(1):161–72. doi: 10.1016/j.cell.2015.01.036
11. Trinder M, Francis GA, Brunham LR. Association of monogenic vs polygenic hypercholesterolemia with risk of atherosclerotic cardiovascular disease. JAMA Cardiol. (2020) 5(4):390–9. doi: 10.1001/jamacardio.2019.5954
12. Smith GD, Ebrahim S. ‘Mendelian randomization’ can genetic epidemiology contribute to understanding environmental determinants of disease? Int J Epidemiol. (2003) 32(1):1–22. doi: 10.1093/ije/dyg070
13. Lawlor DA, Harbord RM, Sterne JA, Timpson N, Davey Smith G. Mendelian randomization: using genes as instruments for making causal inferences in epidemiology. Stat Med. (2008) 27(8):1133–63. doi: 10.1002/sim.3034
14. Smith GD, Ebrahim S. What can Mendelian randomisation tell US about modifiable behavioural and environmental exposures? Br Med J. (2005) 330(7499):1076–9. doi: 10.1136/bmj.330.7499.1076
15. Burgess S, Scott RA, Timpson NJ, Davey Smith G, Thompson SG, Consortium E-I. Using published data in Mendelian randomization: a blueprint for efficient identification of causal risk factors. Eur J Epidemiol. (2015) 30(7):543–52. doi: 10.1007/s10654-015-0011-z
16. 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:2–10. doi: 10.7554/eLife.34408
17. Hartwig FP, Davies NM, Hemani G, Davey Smith G. Two-sample Mendelian randomization: avoiding the downsides of a powerful, widely applicable but potentially fallible technique. Int J Epidemiol. (2016) 45(6):1717–26. doi: 10.1093/ije/dyx028
18. Pierce BL, Ahsan H, Vanderweele TJ. Power and instrument strength requirements for Mendelian randomization studies using multiple genetic variants. Int J Epidemiol. (2011) 40(3):740–52. doi: 10.1093/ije/dyq151
19. Papadimitriou N, Dimou N, Tsilidis KK, Banbury B, Martin RM, Lewis SJ, et al. Physical activity and risks of breast and colorectal cancer: a Mendelian randomisation analysis. Nat Commun. (2020) 11(1):597. doi: 10.1038/s41467-020-14389-8
20. Plagnol V, Smyth DJ, Todd JA, Clayton DG. Statistical independence of the colocalized association signals for type 1 diabetes and RPS26 gene expression on chromosome 12q13. Biostatistics. (2009) 10(2):327–34. doi: 10.1093/biostatistics/kxn039
21. Johnston TP, Korolenko TA, Pirro M, Sahebkar A. Preventing cardiovascular heart disease: promising nutraceutical and non-nutraceutical treatments for cholesterol management. Pharmacol Res. (2017) 120:219–25. doi: 10.1016/j.phrs.2017.04.008
22. Collado A, Domingo E, Piqueras L, Sanz MJ. Primary hypercholesterolemia and development of cardiovascular disorders: cellular and molecular mechanisms involved in low-grade systemic inflammation and endothelial dysfunction. Int J Biochem Cell Biol. (2021) 139:106066. doi: 10.1016/j.biocel.2021.106066
23. Drechsler M, Megens RT, van Zandvoort M, Weber C, Soehnlein O. Hyperlipidemia-triggered neutrophilia promotes early atherosclerosis. Circulation. (2010) 122(18):1837–45. doi: 10.1161/CIRCULATIONAHA.110.961714
24. Kim CS, Park HS, Kawada T, Kim JH, Lim D, Hubbard NE, et al. Circulating levels of MCP-1 and IL-8 are elevated in human obese subjects and associated with obesity-related parameters. Int J Obes (Lond). (2006) 30(9):1347–55. doi: 10.1038/sj.ijo.0803259
25. Kingwell BA, Ahimastos AA. Arterial stiffness and coronary ischemic disease. Adv Cardiol. (2007) 44:125–38. doi: 10.1159/000096725
26. Adams MR, Kinlay S, Blake GJ, Orford JL, Ganz P, Selwyn AP. Atherogenic lipids and endothelial dysfunction: mechanisms in the genesis of ischemic syndromes. Annu Rev Med. (2000) 51:149–67. doi: 10.1146/annurev.med.51.1.149
27. Kwon D, Yi JJ, Ohrr H, Yi SW. Total cholesterol and mortality from ischemic heart disease and overall cardiovascular disease in Korean adults. Medicine (Baltimore). (2019) 98(36):e17013. doi: 10.1097/MD.0000000000017013
28. Hasselbalch RB, Kristensen JH, Strandkjaer N, Jorgensen N, Bundgaard H, Malmendal A, et al. Metabolomics of early myocardial ischemia. Metabolomics. (2023) 19(4):33. doi: 10.1007/s11306-023-01999-8
29. Tan J, Che Y, Liu Y, Hu J, Wang W, Hu L, et al. CELSR2 deficiency suppresses lipid accumulation in hepatocyte by impairing the UPR and elevating ROS level. FASEB J. (2021) 35(10):e21908.34478580
30. Alannan M, Trezeguet V, Amoedo ND, Rossignol R, Mahfouf W, Rezvani HR, et al. Rewiring lipid metabolism by targeting PCSK9 and HMGCR to treat liver cancer. Cancers (Basel). (2022) 15(1):8–16. doi: 10.3390/cancers15010003
31. Lamon-Fava S, Diffenderfer MR, Marcovina SM. Lipoprotein(a) metabolism. Curr Opin Lipidol. (2014) 25(3):189–93. doi: 10.1097/MOL.0000000000000070
32. Rebeck GW. The role of APOE on lipid homeostasis and inflammation in normal brains. J Lipid Res. (2017) 58(8):1493–9. doi: 10.1194/jlr.R075408
33. Da Dalt L, Castiglioni L, Baragetti A, Audano M, Svecla M, Bonacina F, et al. PCSK9 deficiency rewires heart metabolism and drives heart failure with preserved ejection fraction. Eur Heart J. (2021) 42(32):3078–90. doi: 10.1093/eurheartj/ehab431
34. Zhong S, Li L, Shen X, Li Q, Xu W, Wang X, et al. An update on lipid oxidation and inflammation in cardiovascular diseases. Free Radic Biol Med. (2019) 144:266–78. doi: 10.1016/j.freeradbiomed.2019.03.036
35. Marais AD. Apolipoprotein E in lipoprotein metabolism, health and cardiovascular disease. Pathology. (2019) 51(2):165–76. doi: 10.1016/j.pathol.2018.11.002
36. Mach F, Baigent C, Catapano AL, Koskinas KC, Casula M, Badimon L, et al. 2019 ESC/EAS guidelines for the management of dyslipidaemias: lipid modification to reduce cardiovascular risk. Eur Heart J. (2020) 41(1):111–88. doi: 10.1093/eurheartj/ehz455
37. Liu C, Chen J, Chen H, Zhang T, He D, Luo Q, et al. PCSK9 inhibition: from current advances to evolving future. Cells. (2022) 11(19):6–13.
38. Bergeron N, Phan BA, Ding Y, Fong A, Krauss RM. Proprotein convertase subtilisin/kexin type 9 inhibition: a new therapeutic mechanism for reducing cardiovascular disease risk. Circulation. (2015) 132(17):1648–66. doi: 10.1161/CIRCULATIONAHA.115.016080
39. Pasta A, Cremonini AL, Pisciotta L, Buscaglia A, Porto I, Barra F, et al. PCSK9 inhibitors for treating hypercholesterolemia. Expert Opin Pharmacother. (2020) 21(3):353–63. doi: 10.1080/14656566.2019.1702970
40. Gencer B, Mach F. PCSK9 inhibition could be effective for acute myocardial infarction. Curr Med Chem. (2022) 29(6):1016–26. doi: 10.2174/0929867328666210804091003
41. Habibi D, Daneshpour MS, Asgarian S, Kohansal K, Hadaegh F, Mansourian M, et al. Effect of C-reactive protein on the risk of heart failure: a Mendelian randomization study. BMC Cardiovasc Disord. (2023) 23(1):112. doi: 10.1186/s12872-023-03149-3
Keywords: hypercholesterolemia, ischemic heart disease, Mendelian randomization, bidirectional causal relationship, genetic co-localization analysis
Citation: Jiang Y, Yu W, Zhou J and Dong X (2023) Bidirectional causal relationship between hypercholesterolemia and ischemic heart disease: a Mendelian randomization study. Front. Cardiovasc. Med. 10:1302282. doi: 10.3389/fcvm.2023.1302282
Received: 27 September 2023; Accepted: 27 November 2023;
Published: 7 December 2023.
Edited by:
Mohamad Navab, UCLA Health System, United StatesReviewed by:
Eva Hurt-Camejo, AstraZeneca, SwedenDavid John Vigerust, Vanderbilt University, United States
© 2023 Jiang, Yu, Zhou and Dong. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jianliang Zhou zhoujianliang2010@163.com Xiao Dong Dongshi22000@163.com
†These authors share first authorship
Abbreviations IHD, ischemic heart disease; OR, odds ratio; SNPs, several single-nucleotide polymorphisms; IVW, inverse variance weighted; WM, weighted median; GWAS, genome-wide association study; IVs, instrumental variables; MR, mendelian randomization; CELSR2, cadherin EGF LAG seven-pass G-type receptor 2; PCSK9, proprotein convertase subtilisin/kexin type 9; LPA, lipoprotein(a); APOE, apolipoprotein E.