Abstract
Background:
Observational studies have shown inconsistent findings of the relationships between aspirin use and the risk of Alzheimer’s disease (AD).
Objective:
Since residual confounding and reverse causality were challenging issues inherent in observational studies, we conducted a 2-sample Mendelian randomization analysis (MR) to investigate whether aspirin use was causally associated with the risk of AD.
Methods:
We conducted 2-sample MR analyses utilizing summary genetic association statistics to estimate the potential causal relationship between aspirin use and AD. Single-nucleotide variants associated with aspirin use in a genome-wide association study (GWAS) of UK Biobank were considered as genetic proxies for aspirin use. The GWAS summary-level data of AD were derived from a meta-analysis of GWAS data from the International Genomics of Alzheimer’s Project (IGAP) stage I.
Results:
Univariable MR analysis based on these two large GWAS data sources showed that genetically proxied aspirin use was associated with a decreased risk of AD (Odds Ratio (OR): 0.87; 95%CI: 0.77–0.99). In multivariate MR analyses, the causal estimates remained significant after adjusting for chronic pain, inflammation, heart failure (OR = 0.88, 95%CI = 0.78–0.98), or stroke (OR = 0.87, 95%CI = 0.77–0.99), but was attenuated when adjusting for coronary heart disease, blood pressure, and blood lipids.
Conclusion:
Findings from this MR analysis suggest a genetic protective effect of aspirin use on AD, possibly influenced by coronary heart disease, blood pressure, and lipid levels.
Keywords
INTRODUCTION
With an estimated 6.2 million Americans over 65 years of age suffering from AD in 2021, the Alzheimer’s Association has reported $355 billion in total health care spending for people with dementia. It is estimated that Alzheimer’s disease (AD) will affect 13.8 million Americans in 2060, making it a major health challenge in the US [1]. In June 2021, the US Food and Drug Administration (FDA) approved aducanumab, a monoclonal antibody, as the first-ever-disease-modifying drug for AD [2]. Although aducanumab could reduce amyloid plaques [3], what remains controversial is whether the amyloid clearance protects patients from cognitive and functional decline [4]. In this scenario, targeting modifiable risk factors is a promising strategy to prevent AD. Growing evidence has indicated that neuroinflammation and blood clots play important roles in the progression of AD [5, 6]. Aspirin, also known as acetylsalicylic acid, is a commonly used medication for anti-inflammatory and antiplatelet therapies [7].
Aspirin use may be associated with decreased risk of AD because of its anti-inflammatory and cardioprotective properties [8] (http://www.alzdiscovery.org). However, observational study and clinical trial have yielded inconsistent and often contradictory findings. Etminan et al. [9] conducted a meta-analysis of 3 case-control and 5 cohort studies and reported that aspirin had no significant association with the risk of developing AD. Nilsson et al. [10] conducted a cross-sectional analysis which indicated that high-dose aspirin could significantly reduce AD risk for people aged≥80 years. Szekely et al. [11] pooled six prospective studies to report a reduced risk of AD in aspirin users. Subsequent large prospective studies have reported either no effect [12 –15], increased risk [16, 17], or reduced risk [18, 19] of aspirin on AD. The ASPREE (ASPirin in Reducing Events in the Elderly) clinical trial [20] indicated low-dose aspirin was not associated with decreased risk of developing AD in healthy elderly individuals. A recent observational study found that long-term low-dose aspirin use can slow the development of AD in patients with coronary heart disease (CHD) but not in other individuals [19]. Another retrospective cohort study presented that aspirin use may present AD in patients with ischemic stroke [21]. These inconsistent findings demonstrated the inherent limitations of observational studies, the highly heterogenous nature of AD, and the possibility that aspirin use is associated with decreased risk of AD in specific populations which are not included in existing clinical trials.
In this study, we assess potential causal associations between aspirin use and AD through Mendelian randomization (MR) analysis based on genome-wide association study (GWAS) summary statistics. MR analysis utilizes single-nucleotide variants (SNVs) as unconfounded genetic proxies for exposure and infers their causal effects on outcomes, which can minimize the bias encountered in observational epidemiologic studies [22 –24]. MR analysis has been used to assess the relationships between aspirin use and risk of lung cancer [25] and between opioid medication and risk of major depressive disorder, anxiety, and stress-related disorder, and cardiovascular disease (CVD) [26, 27]. We conducted a two-sample univariable MR (UVMR) study utilizing summary-level data of 175,936 human subjects from large-scale GWAS studies that included self-reported aspirin use in the UK Biobank [28], as well as AD genetic associations from the International Genomics of Alzheimer’s Project [29]. We conducted multivariable (MVMR) analyses to investigate whether the causal relationship between aspirin use and AD is influenced by pain, inflammation, CVD, or other risk factors. To the best of our knowledge, this is the first large-scale MR analysis investigating potential causal relationships between aspirin use and AD based on genetics.
METHODS
All analyses used publicly available GWAS summary-level data, so ethical approval from an institutional review board was not necessary. Figure 1 presented the conceptual MR design based on three core assumptions: 1) instrumental SNVs have strong associations with the exposure variable; 2) instrumental SNVs do not associate with any confounders between exposure and outcome; 3) instrumental SNVs do not have a direct influence on outcome [22 –24]. As shown in Fig. 1, the path signed with a black cross is blocked by removing horizontally pleiotropic SNVs, which ensures the validity of the instrumental variable assumptions [30] to estimate total effect without confounding bias. The paths signed with red crosses are blocked by MVMR, which estimates direct effect, so that we could conduct mediation analysis by jointly considering direct and total effects [31]. The relationship (distortion or total effect) estimated with UVMR tends to be reduced if MVMR adjusts for confounders or mediators. In this study, removal of horizontally pleiotropic SNVs eliminates any confounding bias in UVMR; therefore, UVMR and MVMR can be combined to conduct mediation analysis [31, 32].

Overview of the MR analysis for aspirin use and risk of AD.
GWAS data sources
The details of GWAS data sources are listed in Table 1. All data used were publicly available and obtained from European ancestry population. We retrieved the genetic associations of aspirin (acetylsalicylic acid, Anatomical Therapeutic Chemical (ATC) code: N02BA01) use (ASA) from a case-control GWAS performed among United Kingdom Biobank (UKB) study participants [28], which comprised of 61,583 cases and 50,427 controls. A case group included participants who take ASA or medications of the same ATC levels (N02BA01), while a control group included participants who take neither ASA nor medications of the same first two ATC levels [28]. We also used summary statistics of AD from a GWAS meta-analysis of the International Genomics of Alzheimer’s Project (IGAP) [29] stage I. The GWAS was based on a European ancestry meta-analysis of 21,982 cases and 41,944 controls.
GWAS Data Sources and Description from Open GWAS (https://gwas.mrcieu.ac.uk) and GWAS Catalog (http://www.ebi.ac.uk/gwas/)
AD, Alzheimer’s disease; ASA, acetylsalicylic acid; CHD, coronary heart disease; CRP, C-reactive protein; CVD, cardiovascular disease; DBP, diastolic blood pressure; GLGC, Global Lipids Genetics Consortium; HDL, high-density lipoprotein; HF, heart failure; ICBP, International Consortium for Blood Pressure; IGAP, International Consortium for Blood Pressure; LDL, low-density lipoprotein; MRC-IEU, MRC Integrative Epidemiology Unit; SBP, systolic blood pressure.
In general, pain, inflammation, blood pressure, and lipid levels could be considered as risk factors for neurocognitive disorders [33 –35]. To test effect of chronic pain, we used summary-level data from MRC-IEU GWASs for multisite pain (back, knee, hip, neck/shoulder, headache, abdominal/stomach) lasting longer than 3 months [36]. For C-reactive protein (CRP), a sensitive blood biomarker of inflammation, we used summary statistics of serum CRP levels from 61,623 participants included in the within family GWAS consortium [37]. For CVD, we obtained summary statistics for CHD, heart failure (HF), and stroke from several GWAS data resources. The CHD GWAS data were obtained from the Coronary Artery Disease Genome wide Replication and Meta-analysis (CARDIoGRAM) plus the Coronary Artery Disease (C4D) Genetics consortium (CARDIoGRAMplusC4D). To ensure that all participants were of European descent, we used summary data of CARDIoGRAM meta-analysis from 22 GWAS studies [38], which can minimize the potential bias due to population stratification. The summary statistics for HF were retrieved from the Heart Failure Molecular Epidemiology for Therapeutic targets (HERMES) consortium [39]. This GWAS meta-analysis pooled 17 population cohorts and 9 case-control studies. We used summary statistics for stroke (ischemic stroke, intracerebral hemorrhage, and stroke of unknown type) from the MEGASTROKE consortium [40]. The systolic blood pressure and diastolic blood pressure (DBP) traits were defined using a GWAS from the International Consortium for Blood Pressure [41]. These data were from a meta-analysis of 757,601 European participants. Blood lipids including high-density lipoprotein cholesterol, low-density lipoprotein (LDL) cholesterol, triglycerides, and total cholesterol were defined using a GWAS from the Global Lipids Genetics Consortium [42]. Overall, the ASA GWAS may have used overlapping samples with secondary exposure but did not overlap with the sample used to estimate the SNV-AD (outcome) associations.
Univariable Mendelian randomization
We identified a set of genetic variants associated with ASA (p < 5*10–6). Then, the PLINK clumping algorithm was employed to conduct linkage disequilibrium clumping based on the 1000 Genomes Project, to obtain independent SNVs. In this study, relaxed thresholds (r 2 = 0.1 and window size = 100 kb) were selected, which has been used in previous MR research when few significant SNVs are extracted [43, 44]. A genetic variant was considered sufficient if the corresponding instrument strength quantified with F-statistic was larger than 10. SNV-exposure and SNV-outcome data were harmonized, and palindromic SNVs [45] were removed (Supplementary Table 1). Prior to the analysis, SNVs exhibiting reverse causation were identified and removed with Steiger filtering [46] based on the harmonized data. The MR-PRESSO (MR-pleiotropy residual sum and outlier) global test and outlier test were used for evaluating the presence of horizontal pleiotropy and for outlier removal respectively [47, 48]. The remaining SNVs were utilized as genetic instruments to assess the causal effect of ASA on the risk of developing AD.
For UVMR analysis, we conducted the inverse variance weighted (IVW) method that combines the Wald ratio to determine the causal effect of ASA on AD risk. The IVW method is an efficient and consistent estimate of causal effect if all selected genetic instruments are valid [49]. Moreover, we performed the iterative Mendelian randomization and pleiotropy (IMRP) analysis [50], which could simultaneously search for horizontal pleiotropic variants and estimate the causal effect. We conducted sensitivity analysis to examine any violation of MR assumptions. Cochran’s Q statistic [51] was used to evaluate heterogeneity between single genetic variants that may indicate the presence of invalid instruments. Heterogeneity was considered present if the Cochran Q-test p-value was less than 0.05. Excessive heterogeneity indicates the violation of assumptions for selecting valid SNVs, which would induce horizontal pleiotropy [52]. We also inspected horizontal pleiotropy with the MR-Egger intercept, which can be regarded as the average pleiotropic effect of genetic variants [53]. The MR Steiger directionality test was used to examine whether the direction in our hypothesis is valid [46]. We implemented all UVMR analyses using the TwoSampleMR package in R.
Multivariable Mendelian randomization
MVMR makes causal inference of each exposure on the outcome conditional on other exposures [54]. UVMR and MVMR estimate total and direct effects, respectively, so that combining these two MR analyses can detect if a mediator is present between an exposure and an outcome [31]. For MVMR analysis, we identified a set of genetic variants that were related to either ASA, chronic pain, CRP levels, CVD, blood pressure, or blood lipids at a significance threshold of p < 5*10–6 or 5*10–8 (Supplementary Table 2). The PLINK clumping algorithm was executed to extract independent SNVs. Appropriate proxy SNV was selected at a threshold of LD r2 of 0.8, when an instrumental SNV-trait association is available in either primary exposure or secondary exposure. Moreover, overall SNVs were ensured to be independent with an additional clumping step (LD r 2 < 0.1, window size > 100 kb). Palindromic SNVs and SNVs exhibiting reverse causation that were identified with Steiger filtering were eliminated. The extended version of the IVW model (MV IVW) [55] in MVMR analysis was utilized as the primary method to estimate causal effect when adjusting for secondary exposure. The generalized version of Cochran’s statistical test was used to inspect potential heterogeneity of genetic variants. MV MR-Egger intercept was used to inspect horizontal pleiotropy [56]. All MVMR analyses were carried out using the MVMR and TwoSampleMR packages in R.
RESULTS
In the present study, 59 SNVs were associated with ASA at a threshold of p < 5*10–6, among which one (rs77507211) was absent in the AD GWAS summary statistics. No appropriate proxy SNVs at the threshold of LD r 2 > 0.8 were identified. After Steiger filtering and removal of palindromic SNVs, 40 SNVs were available for analysis. F-statistics of selected SNVs ranged from 20.9 to 60.7, indicating there are no weak instruments in the MR analysis. MR-PRESSO did not detect any potential pleiotropic outliers.
Aspirin use is causally associated with decreased risk for AD based on UVMR analysis
UVMR analysis revealed that ASA is associated with a lower risk of AD (OR = 0.87; 95%CI = 0.78–0.97, 40 SNVs, Table 2). Cochran’s Q-test (p > 0.05) uncovered no heterogeneity in the IVW model (Table 3). The MR-Egger intercept (p > 0.05) revealed no directional pleiotropy in our analysis. Among 40 SNVs for aspirin use, 3 SNVs (rs1510226, rs1831733, rs583104) are associated with HF, 2 SNVs (rs28752929, rs635634) are associated with DBP, and 6 SNVs (rs117733303, rs1510226, rs4299376, rs4927207, rs583104, rs635634) are associated with LDL cholesterol. After removing potential pleiotropic SNVs related to LDL cholesterol, the association between aspirin use and AD risk remained significant (OR = 0.87, 95%CI = 0.77–0.99; Fig. 2, Tables 2 and 3). While there is a yet-to-confirm relationship between AD and HF (or DBP), our study showed a causal estimate of aspirin use on the decreased risk of AD (removal of HF SNVs: OR = 0.83, 95%CI = 0.74–0.93; removal of DBP SNVs: OR = 0.87, 95%CI = 0.78–0.98; removal of all potentially pleiotropic SNVs: OR = 0.85, 95%CI = 0.75–0.96, Tables 2 and 3). Moreover, leave-one-out analysis suggested the estimate effect of IVW was independent of any single SNV. Consistent results were obtained using a more stringent threshold of LD r 2 < 0.01 (Supplementary Table 3). The MR Steiger directionality test showed that the assumption that aspirin use causes AD is valid (Table 3). IMRP provides evidence of no horizonal pleiotropy and significantly protective effect of ASA on risk of AD (Tables 2 and 3). In the frailty analysis (Supplementary Methods 1), the IVW method estimated the mean effect induced by frailty alone (OR = 0.998, 95%CI = 0.92–1.08, Supplementary Figure 1), indicating that frailty effects negligibly influence the true effect estimate (OR = 0.87, 95%CI = 0.77–0.99).
The causal estimate of ASA on AD from IVW and IMRP model
ASA1, All ASAs SNVs; ASA2, ASAs SNVs that are not associated with HF; ASA3, ASAs SNVs that are not associated with DBP; ASA4, ASAs SNVs that are not associated with LDL cholesterol; ASA5, ASAs SNVs that are not associated with HF, DBP, or LDL cholesterol. AD, Alzheimer’s disease; ASA, acetylsalicylic acid; CI, confidence interval; DBP, diastolic blood pressure; HF, heart failure; IMRP, iterative Mendelian randomization and pleiotropy; IVW, inverse variance weighted; LDL, low-density lipoprotein; OR, odds ratio; SNVs, single-nucleotide variants.
Sensitivity analyses with different methods

Analysis of the causal relationship between ASA and risk of AD.
The protective effects of aspirin use on AD is partially influenced by coronary heart disease, blood pressure and lipid levels based on MVMR analysis
In MVMR, we examined whether the association between ASA and decreased risk of AD was influenced by chronic pain, inflammation, CVD, blood pressure, and blood lipids. Because there is a yet-to-confirm relationship between AD and HF (or DBP), we retained ASA SNVs associated with HF or DBP to conduct the MVMR analysis. When adjusting for the pain- and inflammation-related traits, the genetic causality of ASA on AD remained statistically significant (adjusting for chronic back pain: OR = 0.87, 95%CI = 0.77–0.99; chronic headache: OR = 0.88, 95%CI = 0.78–0.999; CRP levels: OR = 0.84, 95%CI = 0.73–0.96) (Fig. 3, Supplementary Table 4). Potential heterogeneity (p > 0.05 in Cochran’s statistical test) was not detected (Supplementary Table 4). MVMR Egger intercepts suggest no horizontal pleiotropy. The causal estimates between ASA and AD remained statistically significant when adjusting for HF (OR = 0.88, 95%CI = 0.78–0.98) and stroke (OR = 0.87, 95%CI = 0.77–0.99), while the relationship between ASA and AD risk reduction was attenuated when adjusting for CHD (OR = 1.0, 95%CI = 0.85–1.17), blood pressure (OR = 0.95, 95%CI = 0.87–1.04) and blood lipids (OR = 0.89, 95%CI = 0.76–1.03). These results suggest that the protective association between aspirin and AD is partially influenced by CHD, blood pressure, and blood lipids. Again, MV Egger intercepts suggest no horizontal pleiotropy.

MVMR results of ASA and risk of AD.
DISCUSSION
Our UVMR study suggests a causally protective effect of ASA on the risk of developing AD. These results are consistent with some previous findings [18]. Moreover, MVMR analysis indicated that ASA’s effects on AD might be influenced by CHD, blood pressure, and blood lipids.
This two sample MR analysis have several strengths. MR analysis minimizes bias caused by reverse causality and residual confounding, which are common in observational studies. In addition, all SNV effect estimates in this study were obtained from GWAS studies restricted to populations of European descent, therefore, bias due to population stratification was minimized [49, 57]. We performed rigorous sensitivity analyses to assess any violation of MR assumptions. We identified eight SNVs of rs1510226, rs1831733, rs583104, rs28752929, rs635634, rs117733303, rs4299376, and rs4927207 that were associated with heart failure, diastolic blood pressure, or LDL cholesterol. Growing evidence suggests that patients with CVD are at increased risk of AD and that there is a strong causal association between CVD, cognitive decline, and AD [6]. After removal of eight CVD-associated SNVs, UVMR still showed a significant causal relationship between ASA and AD, suggesting that protective effects of aspirin on AD are not biased. Finally, combined UVMR and MVMR analyses revealed potential mediators underlying the observed causal association between ASA and AD. MVMR results suggest that the effect of ASA on AD risk reduction is independent of pain and inflammation. MVMR analysis further indicated that the causal protective association of ASA on AD might be influenced by CHD, blood pressure, and blood lipid levels. Previous studies have reported that CHD, blood pressure, and blood lipids are associated with a risk of AD [6 , 59]. Clinical trials have also demonstrated the strong cardioprotective properties of aspirin; currently, aspirin is recommended to decrease the risk of adverse cardiovascular events and all-cause mortality in patients with higher cardiovascular risk and lower risk of bleeding [60].
Although the exact causes of AD remain largely unknown, the “amyloid hypothesis” suggests that impaired clearance of amyloid-beta protein plaques is implicated in AD pathogenesis [61]. Thus, activating or boosting the cellular pathways in disposing of waste in the brain could be a promising strategy to slow or delay the progression of AD. Transcription factor EB (TFEB) regulates the expression of the brain’s debris-clearing machinery. A previous study revealed that aspirin use upregulated TFEB, stimulated the production of lysosomes, and decreased amyloid plaque pathology in AD model mice [62].
Epidemiological and clinical trial studies showed inconsistent and often contradictory relationships between aspirin use and risk of AD. Observational studies susceptible to residual confounding and bias. In addition, the protective effect of aspirin use on AD may be limited to specific subgroups of individuals. For example, a previous observational study has suggested aspirin is associated with reduced cognitive decline in patients with stroke or CHD [19, 63]. ASPREE trial that focused on healthy elderly participants [20] may have not detected protective effect of aspirin use on reducing the risk of AD in specific high risk populations. Our study, by leveraging summary genetics from large-scale GWAS studies, provides evidence that aspirin use may be causally associated with decreased risk for AD, partially influenced by CHD, blood pressure, and lipid levels. It would probably mean that aspirin can provide a more significant protective effect in patients at higher risk for CHD, abnormal blood pressure or lipid levels.
This study has several limitations. First, potential dose effects of aspirin use on the risk of AD is not considered in our two-sample MR analysis, however it is important because any benefit needs to be balanced against the risks of excessive bleeding. Moreover, the adverse effects of aspirin change with age, leading to recommendations for cardiovascular prophylaxis that are both risk and age dependent [60, 64]. Future work in stratified populations (e.g., according to ASA dose, age, comorbidities, etc.) is warranted to further investigate the heterogeneity of the effect of ASA on AD in specific sub-populations [65], which could be a basis of clinical studies. Second, to compensate for the lack of SNVs, we selected SNVs related to aspirin use and chronic pain risk at a relaxed threshold (p < 5*10–6), which is less than the conventional genome-wide significance threshold (p < 5*10–8). However, selecting SNVs with smaller effect sizes for ASA may increase the potential of weak instrument bias, although we used F-statistics to validate the instrument strength. A relaxed threshold was also selected to perform LD clumping, ensuring instrument independence. Although no horizontal pleiotropy was detected by sensitivity analyses, it remains a possibility that existing pleiotropy caused bias. Third, the ASA GWAS summary statistics from individual GWAS studies may have limited power to detect genetic associations, in which the effect size of a significant SNV-exposure association tends to be overestimated, leading to the winner’s curse phenomenon [66]. Thus, the causal estimate in our study might be biased with the data-driven selection of instrumental variants [67]. Fourth, MR analysis only estimated effect sizes corresponding to a long-term difference. Work is needed to examine the effect size from a short-term intervention [68]. Fifth, the reliability of the employed MR investigation depends on the validity of the instrumental SNVs. In this study, we used several strategies to assess the validity of instruments and remove invalid ones. Future work is necessary to examine how other MR designs may impact the findings in our study, including identifying instrumental SNVs in the genes of the targeted protein to mimic the mechanism of action or performing robustly MR estimates in the presence of invalid instruments. Finally, this study was based on summary statistics of GWAS studies. While there was no sample overlap between exposure and outcome GWAS in UVMR, primary exposure may have overlapping samples with secondary exposure that had potential in distorting the estimated causal relationship of ASA with AD. Therefore, future work is warranted to perform individual-level data MVMR analysis for preventing the sample overlap problem in the two-sample summary data setting.
CONCLUSIONS
This study leveraged 2-sample MR analyses to demonstrate a genetic causality between ASA and decreased risk of AD. We further showed that this causal relationship might be influenced by blood pressure, blood lipid levels, and CHD. Findings from this study provide evidence supporting the benefit of aspirin use in reducing the risk of AD, which may help identify patient populations to curb AD by prevention and treatment strategies.
Footnotes
ACKNOWLEDGMENTS
The authors have no acknowledgments to report.
FUNDING
This work has been supported by the NIH National Institute of Aging (grants nos. AG057557, AG061388, AG062272, AG076649) and the Clinical and Translational Science Collaborative (CTSC) of Cleveland (TR002548).
CONFLICT OF INTEREST
The authors have no conflict of interest to report.
DATA AVAILABILITY
The data supporting the findings of this study are available within the article.
