DOI: https://doi.org/10.21203/rs.3.rs-2823753/v1
Background
At present, many studies have shown a significant correlation between myocardial infarction (MI) and atrial fibrillation (AF), but few focus on the bidirectional causal relationship between MI and AF. Therefore, this Mendelian randomization (MR) study was designed to examine the bidirectional causality between MI and AF.
Methods
We used the publicly available summary statistical dataset of MI from genome-wide analysis studies (GWAS; ebi-a-GCST011364; case = 14,825, control = 2,680). The summary statistical dataset of AF was obtained from a European population GWAS (finn-b-I9_AF_REIMB; case = 10,516, control = 116,926). A two-sample bidirectional MR analysis was performed using analysis methods including inverse-variance weighted (IVW), MR-Egger, and weighted median.
Results
In site-specific MI analyses, we screened 30 single nucleotide polymorphisms (SNPs) from GWAS as instrumental variables (IVs). Causal association between MI and AF can be supported by IVW (beta = 0.349, SE = 0.057, P < 0.001), MR-Egger (beta = 0.398, SE = 0.131, P = 0.005), and weighted median (beta = 0.352, SE = 0.068, P < 0.001). In the reverse MR analyses, we screened 20 SNPs as IVs and the casual effect of AF on MI was observed by IVW (beta = 0.047, SE = 0.022, P = 0.033).
Conclusions
Our MR study results showed a good causal effect of MI on AF. The MR analysis also found a causal effect of AF on MI.
Myocardial infarction (MI) and atrial fibrillation (AF) are two of the most common cardiovascular diseases worldwide(1). MI is caused by the blockage of coronary arteries leading to the death of cardiac muscle cells, while AF is an arrhythmia characterized by irregular and rapid heartbeat(2, 3). These diseases have a significant impact on morbidity and mortality, with MI being the leading cause of death worldwide and AF increasing the risk of stroke and heart failure(4, 5).
Several risk factors have been identified for both MI and AF, including hypertension, diabetes, smoking, and obesity(6, 7). However, the causal relationship between the two diseases is still unclear, and it is challenging to establish causality due to the complexity and heterogeneity of cardiovascular diseases(8, 9).
Mendelian randomization (MR) is a method that utilizes genetic variants as instrumental variables (IVs) to infer causality between exposures and outcomes(10). This method has become increasingly popular in epidemiological studies because it can provide more reliable causal inference than traditional observational studies(11). MR can also be used to investigate the causal relationship between modifiable risk factors and disease outcomes, which could have important implications for prevention and treatment(12).
In recent years, MR studies have shown that genetic variants associated with modifiable risk factors, such as smoking being associated with MI and hypertension being associated with AF(13, 14). These findings suggest a potential causal relationship between these risk factors and the diseases. However, previous MR studies have mostly focused on the unidirectional causal relationship between the risk factors and the diseases. The bidirectional causal relationship between MI and AF has not been thoroughly investigated using MR(15, 16).
Therefore, in this study, we investigate the bidirectional causal relationship between MI and AF using MR. We will use genetic variants associated with MI or AF as IVs to estimate the causal effect of one disease on the other. This approach will enable us to investigate the potential causal relationship between MI and AF and provide insights into the underlying biological mechanisms.
Overall, our study aims to contribute to a better understanding of the causal relationship between MI and AF, which could have important implications for the prevention and treatment of these diseases. By identifying causality that contributes to the development of MI and AF, we may be able to develop more effective strategies for preventing and treating these conditions.
2.1 Overall study design
All data were obtained from published studies approved by the institutional review boards, and informed consent was obtained from the participants of the original study (17). Therefore, no further sanctions were required (18). The cause-and-effect relationship between MI and AF was analyzed using a two-sample MR study, and SNPs were defined as IVs. The use of SNPs to model randomized controlled trials (RCTs) can help identify causal relationships between exposure characteristics (i.e., MI) and outcome characteristics (i.e., AF).
2.2 Data sources
MI data were obtained from the latest and largest published GWAS meta-analysis performed in 2021 by JA et al, including 14,825 MI cases and 2,680 control participants (19). The study was approved by the institutional review board, and informed consent was obtained from all participants in the original study. SNPs were selected based on the following criteria: i) SNPs strongly related to MI with genome-wide significance (p < 5×10− 8). ii) Independent of each other and to avoid bias owing to linkage disequilibrium (LD), the LD of SNPs related to MI had to fulfill r² <0.001, with a window size of 10,000 kb. iii) The correlation between IV and exposure factors is typically determined using the F-statistic of the SNPs. In general, IVs with an F-statistic greater than 10 are regarded as unbiased. F-statistic = (β/SE)2.
Data on AF were obtained from the FinnGen project (FinnGen), which is available at https://gwas.mrcieu.ac.uk/datasets/finn-b-I9_AF_REIMB/ and includes 10,516 cases and116,926 controls of European ancestry.
2.3 Statistical analysis
It is essential to consider these hypotheses to provide a valid explanation for MR analysis (20). (i) It is well-established that IVs are strongly related to MI. (ii) AF is affected only by IVs due to MI defects. (iii) No confounding factors were present in the relationship between MI and AF according to IVs. The results can be affected by genetic variation through a single pathway rather than by separate exposure, namely horizontal pleiotropy, which contradicts the assumptions of MR and may bias the causal estimates. Three different analytical methods were used in the MR analysis to prevent this. Each analysis was based on a different horizontal multiplicity model. The benefit of comparing these three results is that the consistency of the three methods makes the results more credible. The main analysis was performed using an inverse variance weighting (IVW) approach, which provided the most accurate estimates but assumed that all SNPs were valid IVs (21). If one SNP does not meet the IVs assumption, the random-defect IVW will be used to generate a bias, which weighs each rate according to its standard error while considering possible heterogeneity. To satisfy the premise of a valid instrumental variable, the weighted median method requires at least 50% SNPs (22). After sorting the included SNPs based on the weights, we obtained the median of the corresponding distribution function according to the results of our experiments. Additionally, if the genetic instrument does not depend on pleiotropic effects, an effect estimate can be derived from MR-Egger regression (23). The pleiotropic effect was assessed using MR-Egger’s intercept. Furthermore, a directional multiplicative effect cannot be proven if MR-Egger’s intercept does not differ dramatically from zero (24).
2.4 Sensitivity analysis
Funnel plots can plot a single Wald ratio per SNP to display the directional level pleiotropy of the IVs. Nevertheless, the small number of IVs included makes it difficult to test for horizontal pleiotropy using funnel plots. The causal effect of the funnel plot was approximately symmetrical (Fig. 1). Leave-one-out analyses were performed to investigate whether estimates from IVW analyses were biased or dictated by individual SNPs, during meta-analyses that were conducted based on rerun IVW results for the remaining SNPs after omitting one SNP per succession. After removing each SNP, we performed MR analysis again systematically for the remaining SNPs. The results were consistent, indicating a significant causal relationship between the calculated results for all the SNPs (Fig. 2). In MR analysis, the second hypothesis is that SNPs inject results only by modifying the exposure of interest, without other confounding pathways. Directional multidirectionality was examined to obtain the intercept and p-value using MR-Egger regression. No horizontal pleiotropy was observed in the intercept of the MR-Egger regression (p > 0.05), further indicating that pleiotropy did not bias the causal effect. Furthermore, in the published GWAS, there was no evidence that the included MI associated SNPs were significantly associated with any phenotype except MI, which indicates that the assumptions of the third MR were not violated. Consequently, there was no evidence that the genetic instruments of the 30 MI-associated SNPs were significantly associated with any other phenotype on a genome-wide scale, supporting our third MR hypothesis, which is unlikely to be breached in our study. The “Two sample MR” (version 0.5.6) software package was applied for MR and sensitivity analysis in R (version 3.6.2) (25).
3.1 Instrumental variables for MI
The SNPs’ signatures of the MI are shown in Supplementary Table S1. Finally, we selected 30 SNPs as the IVs. All genetic tools related to MI were at a genome-wide significance level (p < 5×10− 8, F > 10). Thus, none of the SNPs was susceptible to IVs. The causal effects of each genetic variant on AF are shown in forest (Figs. 3) and scatter plots (Fig. 4).
3.2 Mendelian randomization analyses for AF
We evaluated the causal relationship between MI levels and AF using IVW, MR-Egger, and weighted median regression (Table 1). Our findings suggest an increased risk of AF in patients with MI (OR = 1.41; 95% CI 1.26 − 1.58, p = 6×10− 10).
MR methods |
nSnp |
Pleiotropy |
Heterogeneity Test |
||||
---|---|---|---|---|---|---|---|
Beta |
SE |
P Value |
Cochran’s Q Statistic |
P Value |
|||
Causal effect of MI on AF |
|||||||
MR Egger |
30 |
0.398 |
0.131 |
0.005 |
50.201 |
0.006 |
|
Weighted median |
30 |
0.352 |
0.068 |
< 0.001 |
|||
Inverse-variance weighted |
30 |
0.349 |
0.057 |
< 0.001 |
50.506 |
0.008 |
|
Simple mode |
30 |
0.365 |
0.122 |
0.006 |
|||
Weighted mode |
30 |
0.365 |
0.089 |
< 0.001 |
|||
Causal effect of AF on MI |
|||||||
MR Egger |
20 |
-0.025 |
0.045 |
0.580 |
29.164 |
0.046 |
|
Weighted median |
20 |
0.017 |
0.024 |
0.485 |
|||
Inverse-variance weighted |
20 |
0.047 |
0.022 |
0.033 |
34.464 |
0.016 |
|
Simple mode |
20 |
0.012 |
0.040 |
0.763 |
|||
Weighted mode |
20 |
0.015 |
0.029 |
0.603 |
|||
MR, mendelian randomization; MI, myocardial infarction; AF, atrial fibrillation; Snp, single nucleotide polymorphism; SE, standard error. |
3.3 Bidirectional Mendelian Randomization
We also sought to explore whether AF influenced MI. Therefore, we reversed the functions of the exposures and outcomes to perform a bidirectional MR analysis and determine the effects of a genetically increased risk of AF on MI. To that end, we selected SNPs that were significant genome-wide (p < 5E − 08) and independently inherited (r2 < 0.01) without LD for AF from FinnGen project (Supplementary Table S2). We then applied the same MR methods as above (Table 1). The statistical tests of the bidirectional MR analysis were two-sided, and the results of the MR analyses and sensitivity analyses regarding the causal effects of AF on MI were considered statistically significant. Our findings also suggest an increased risk of MI in patients with AF (OR = 1.04; 95% CI 1.00 − 1.09, p = 0.03). However, we found heterogeneity in the results, which may be relevant to the included population (p = 0.01). (Fig. 5)
MR is a tool used in causal inference, and its three core assumptions are prerequisites for reliable causal inference. The genetic variant used as an instrument is strongly associated with the exposure, is independent of any confounders, and only affects the outcome through its effect on the exposure. In this study, by using SNPs as the IVs and excluding the bias caused by confounding factors, we elucidated the bidirectional causal relationship between MI and AF using the MR methods. We extracted data on MI and AF from database and ensured that all SNPs selected as IVs were statistically correlated with the exposure factor and not in LD, and we conducted horizontal pleiotropy analysis to exclude SNPs that might be related to confounding factors and eliminated weak IVs. We then harmonized the effect alleles in the GWAS data of the exposure and outcome factors and used five analysis methods to conduct MR, while also testing for heterogeneity and pleiotropy in the IVs. Finally, the study results showed a good causal effect of MI on AF, as evidenced by the consistency and monotonicity of the scatter plot, one-by-one analysis, and distribution of SNPs in the funnel plot. The IVW analysis also found a causal effect of AF on MI but the effect performance is not as good as that of MI on AF.
Our findings are consistent with previous observational studies that have shown a positive association between MI and AF(26). The reasons for MI leading to AF may be related to the following aspects. Experimental studies have found that when the perfusion of the right coronary artery and the left circumflex branch is significantly reduced, the excitability of the atrial myocardium is increased, and the conduction velocity is accelerated, which can cause re-entry and fibrillation(27). After MI occurs, the atrium is excessively stretched which increases the excitement of atrial myocytes and prolongs the length of the electrical conduction pathway, which is conducive to the formation of re-entry and fibrillation(28). In addition, MI-induced myocardial fibrosis also favors the formation of re-entry and fibrillation. Inflammation plays an important role in the occurrence and maintenance of AF. When MI occurs, a large amount of inflammatory factors are released. In the non-infarct and non-ischemic areas, the expression of inflammatory factors also increases, indicating that AF under acute MI conditions may be a marker of widespread inflammation(29). ln addition, the reduced parasympathetic tone, increased sympathetic nerve output, and hormones including B-type natriuretic peptide after MI may be related to the occurrence of AF(30, 31).
The Atherosclerosis Risk in Communities study found that AF patients have a 63% increased risk of acute MI(32). MI The reasons for AF leading to MI may be related to the following aspects. AF is associated with systemic signs of inflammation that could promote a pro-thrombotic state and eventually MI(33). AF and its risk factors, such as hypertension, diabetes, and dyslipidemia, can lead to platelet activation, which is a critical step in the process of MI(34, 35). When the heart rate of AF patients is too fast, it can lead to increased myocardial oxygen consumption and decreased coronary blood flow, which can easily cause type 2 MI due to an imbalance between oxygen supply and demand(36). The traction of the atrial muscle can significantly enhance sympathetic nerve activity, stimulate the release of catecholamines in the heart, and activate adrenergic receptors to cause vasoconstriction(37). However, most of studies are based on laboratory data or observational results. Our study provides stronger evidence for a causal relationship between the two conditions using a MR approach. Additionally, our study is the first to demonstrate a bidirectional causal relationship between MI and AF. The mechanisms underlying the causal relationship between MI and AF remains unclear. Further studies are needed to elucidate the exact mechanisms underlying the causal relationship between MI and AF.
This study has several advantages. MR is used to elucidate the bidirectional causal relationship between MI and AF, providing important clinical implications(38). Clinicians should be aware of the increased risk of AF in patients with a history of MI and consider appropriate management strategies to reduce the risk of AF. Conversely, patients with AF should be screened for underlying cardiovascular diseases, including MI. Additionally, our study reminds the importance of primary prevention of both MI and AF. Moreover, the reliability of causal association evidence based on MR studies is between that of observational epidemiological studies and experimental epidemiological studies(12). In other words, the evidence level of MR studies is higher than that of cohort studies and only slightly lower than that of RCTs studies. However, the high cost of RCTs makes them less feasible. MR can provide powerful evidence when RCTs cannot be implemented, for example, due to ethical reasons. MR studies can effectively overcome biases caused by confounding and reverse causality issues, providing reliable evidence for inferring the causal relationship between exposure factors and outcomes(39). In the MR analysis process, selecting SNPs as IVs can directly infer a causal relationship between exposure and outcome, as SNPs are randomly distributed to individuals through genetic inheritance and are not affected by external environmental and other confounding factors(40).
There are several limitations to our study. First, our study was conducted in individuals of European descent, and the generalizability of our findings to other populations is unclear(41). Second, MR assumes that genetic variants are not associated with confounding factors, and our results could be biased if this assumption is violated. Third, our study only examined the relationship between MI and AF and did not investigate the effect of treatment on either condition. Finally, AF and MI have specific subtypes, such as paroxysmal atrial fibrillation, persistent atrial fibrillation, type 1 MI, and type 2 MI. Due to limitations in the data source, causal relationship analysis cannot be performed on these subtypes.
In future studies, it will be important to investigate the mechanism underlying the causal relationship between MI and AF and to explore potential therapeutic strategies for reducing the risk of both conditions. Additionally, studies in diverse populations are needed to determine the generalizability of our findings. Finally, investigating the effect of treatment on both MI and AF could provide important insights into the management of these conditions.
In conclusion, the causal effect of MI on AF was found by our MR analyses, which indicated that the occurrence of MI may increase the risk of AF in European population. In the reverse MR analyses, the causal effect of AF on MI was also observed.
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.
Author Contributions
JR and ZW contributed to conception and design of the study. JR and YH organized collection and assembly of data. JR and PW performed the statistical analysis. JR wrote the first draft of the manuscript. All authors contributed to manuscript revision, read, and approved the submitted version.
Funding
No funding was involved in this study.
Acknowledgments
This manuscript has been released as a pre-print at Research Square (doi:).
Data Availability Statement
Publicly available datasets (ID: ebi-a-GCST011364 and finn-b-I9_AF_REIMB) were achieved from the website (https://gwas.mrcieu.ac.uk/).