Front. Cardiovasc. Med. Frontiers in Cardiovascular Medicine Front. Cardiovasc. Med. 2297-055X Frontiers Media S.A. 10.3389/fcvm.2021.722154 Cardiovascular Medicine Original Research Physical Activity and Risks of Cardiovascular Diseases: A Mendelian Randomization Study Zhuo Chengui 1 Zhao Jianqiang 2 Chen Miao 3 Lu Yunlong 3 * 1Department of Cardiology, Taizhou Central Hospital (Taizhou University Hospital), Taizhou, China 2Department of Cardiology, The Fourth Affiliated Hospital of Zhejiang University School of Medicine, Zhejiang, China 3Department of Cardiology and Atrial Fibrillation Center, The First Affiliated Hospital, College of Medicine, Zhejiang University, Hangzhou, China

Edited by: Rudolph Schutte, Anglia Ruskin University, United Kingdom

Reviewed by: Petri Böckerman, University of Jyväskylä, Finland; Patricia Kim Phuong Nguyen, Stanford Healthcare, United States

 *Correspondence: Yunlong Lu luyl@zju.edu.cn

This article was submitted to Cardiovascular Epidemiology and Prevention, a section of the journal Frontiers in Cardiovascular Medicine

 29 09 2021 2021 8 722154 08 06 2021 31 08 2021 Copyright © 2021 Zhuo, Zhao, Chen and Lu. 2021 Zhuo, Zhao, Chen and Lu

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.

Background: Although some observational studies have shown that physical activity may have a positive relationship with cardiovascular diseases, the causal effect remains uncertain. We conducted a Mendelian randomization (MR) study to identify the potential causal effect between physical activity and cardiovascular diseases.

Methods: Summary statistics of genome-wide association studies on four physical activity phenotypes and cardiovascular diseases were utilized. MR analysis was performed using inverse-variance weighted (IVW) and multivariable MR. Multiple sensitivity analysis was further conducted to identify the robustness of our results.

Results: Genetically predicted self-reported vigorous physical activity (VPA) was significantly associated with lower risk of myocardial infarction (IVW OR: 0.24, 95% CI: 0.08–0.68, p-value: 0.007). Additionally, the causal effect of VPA with myocardial infarction was robust after adjusting for several cardiovascular risk factors through using the multivariable MR. There were no apparent causal associations between physical activity with other cardiovascular diseases. Results were consistent with the sensitivity analysis.

Conclusion: The present study supports a protective role of self-reported vigorous physical activity in the initiation of myocardial infarction and highlights the importance of activity levels of physical activity. Further studies are required to elucidate the potential biological pathways of physical activity with cardiovascular diseases.

 physical activity cardiovascular disease coronary artery disease myocardial infarction mendelian randomization Background

Cardiovascular diseases (CVDs) represent a leading cause of death worldwide and, despite coordinated preventive efforts, still account for 30% of all global deaths and impose a substantial burden on society (1). From previous epidemiological studies, we have learned that lower systolic blood pressure (2), weight reduction (3), reduction of smoking (4), education (5), and healthy diet (6) would all lower the risk of CVDs. However, there are many other ideal preventive measures that need to be considered.

The association between physical activity and CVDs has attracted much attention in recent years. For example, previous studies have demonstrated that physical activity improves the prognosis of CVDs in both healthy participants and those with traditional risk factors (7–9). The 2016 European Guidelines on CVD prevention in clinical practice recommended that every person accumulate at least 150 min/week of moderate-intensity physical activity or 75 min/week of vigorous-intensity physical activity (10). Nevertheless, most of the evidence for these guidelines originates from observational studies, which cannot be used to identify the causality because of the possibility of confounding and reverse causation.

With regard to the causal relationship, several randomized controlled trials (RCTs) have been applied to show an inverse relationship between physical activity and overall CVD risk (11, 12). However, other RCTs found no such benefits (13, 14). The results of RCTs are still controversial due to the limited sample size and follow-up time. More critically, studies to date rely mostly on self-reported information, which might be influenced by memory inaccuracy and individuals' mood (7). Therefore, the causal associations of physical activity with CVDs remain uncertain.

Mendelian randomization (MR) is an alternative epidemiological approach that can use genetic variations as instrumental variables (IVs) to uncover the causal relationship between exposure and disease outcomes, and avoiding potential confounders in observational studies (15). MR has been previously used to explore the association of several risk factors and CVDs, such as smoking (16), body mass index (BMI) (17), and alcohol consumption (18). In this study, a two-sample MR study was performed to examine the potential causal association between physical activity and risk of CVDs using genetic variations associated with self-reported and objective accelerometer-based physical activity identified from published genome-wide association studies (GWASs).

 Methods Study Design

MR studies are widely adopted to investigate the association between exposure and outcomes. There are three key assumptions for MR: First, the genetic variants should be associated with the physical activity. Second, the genetic variants should not be associated with any confounders. Third, the genetic variants exert effects on the outcome only via the physical activity (Supplementary Figure 1). All summary-level data used in this study come from published GWASs on physical activity and CVDs, including coronary artery disease (CAD), myocardial infarction (MI), heart failure (HF), atrial fibrillation (AF), ischemic stroke (IS), and its subtypes.

 Data on PA

From the available literature, we note that there is a significant genetic influence on physical activity levels (19, 20). Specifically, previous studies have applied both human and mouse models, and the amount of heritability observed has ranged from 0.20 to 0.92. Genetic factors influence activity levels through aerodynamic performance, muscle strength, muscle endurance, and anaerobic performance (21). We used summary-level data from a recently published GWAS on physical activity conducted in participants from the UK Biobank study (22). This GWAS examined four physical activity phenotypes (Supplementary Figure 2) including self-reported moderate-to-vigorous physical activity (MVPA), self-reported vigorous physical activity (VPA), overall acceleration average, and fraction of accelerations > 425 milli-gravities (corresponding to an equivalent level of vigorous physical activity). In the UK Biobank, self-reported physical activity of 377,234 participants during work and leisure time was ascertained through a touchscreen questionnaire that was similar to the International Physical Activity Questionnaire (23). For the measurement of accelerometer-based physical activity, approximately 91,000 participants wore an Axivity-AX3 triaxial accelerometer on their wrist for 7 days (24). The cutoff value of 425 mg was chosen because it corresponds to an equivalent level of vigorous physical activity.

Genome-wide significant (p-value < 5 × 10−8) and independent (linkage disequilibrium r2 < 0.1) single-nucleotide polymorphisms (SNPs) were selected as instruments for the MR analysis. Finally, the GWAS identified nine and five independent genome-wide significant SNPs for MVPA and VPA, respectively. Furthermore, we also selected eight SNPs significantly associated with “overall acceleration average” and two SNPs associated with “fraction of accelerations > 425 milli-gravities (mg).” Given only two SNPs were available for “fraction of accelerations > 425 mg,” a more relaxed significance threshold (p-value < 5 × 10−7) was used and finally eight SNPs were detected. This method of relaxing the threshold, which is also called suggestive significance level, has been used in several previous MR studies when few significant SNPs are available (25, 26). When SNPs for the physical activity phenotypes were not available in the summary statistics of outcome GWAS, proxies (linkage disequilibrium r2 > 0.8) were identified via an online tool, available at: https://ldlink.nci.nih.gov/. Resulting lists of SNPs for each phenotype are given in Supplementary Table 1.

 Data on CVDs

Summary statistics for the associations of the physical activity-related SNPs with CVDs were extracted from the Coronary Artery Disease Genome-Wide Replication and Meta-analysis plus the Coronary Artery Disease Genetics (CardiogramplusC4D) consortium for CAD (60,801 cases and 123,504 controls) and MI (43,676 cases and 128,197 controls) (27), from the Heart Failure Molecular Epidemiology for Therapeutic Targets (HERMES) Consortium for heart failure (47,309 cases and 930, 014 controls) (28), from the Atrial Fibrillation Haplotype Reference Consortium for atrial fibrillation (65,446 cases and 522,744 controls) (29), and from the MEGASTROKE consortium for ischemic stroke and its subtypes (34,217 cases and 404,630 controls) (30). The definition of ischemic stroke was according to the Trial of Org 10172 in Acute Stroke Treatment criteria (31), and its subtypes were categorized as large artery stroke (4,373 cases and 146,392 controls), small vessel stroke cases (5,386 cases and 192,662 controls), and cardioembolic stroke cases (7,193 cases and 204,570 controls).

 Statistical Power

The a priori statistical power was calculated using a web-based application (http://cnsgenomics.com/shiny/mRnd/) (32). We identified that the nine SNPs for MVPA and the five SNPs for VPA explained 0.1% of the phenotypic variability. The eight SNPs for “overall acceleration average” and “fraction of accelerations > 425 mg” explained 0.3% of the phenotypic variability. Given a type 1 error of 5%, the power estimates of the four physical activity phenotypes are shown in Supplementary Table 2.

 Statistical Analysis

A two-sample MR method was used in the present study. We calculated the MR estimates of the effect of physical activity on CVDs through utilizing the Wald estimator. The Delta method was used to account for possible measurement error in both the physical activity and CVD association estimates (32, 33). The fixed-effect inverse-variance-weighted (IVW) method was implemented to evaluate the causal effect between physical activity and the outcomes. Additionally, F statistic was calculated to detect the strength of each instrument using the following formula: F=R2(N-2)(1-R2), where R2 stands for percentage of the variation explained by the SNPs and N is the sample size of the GWAS. To rule out possible pleiotropic effects, we also looked up each selected instrument SNPs in Phenoscanner (34) (http://www.phenoscanner.medschl.cam.ac.uk) to evaluate any previously reported associations (p-value < 5 × 10−6) with CVDs. Several sensitivity analyses were applied to assess robustness of the results. First, heterogeneity was evaluated by Cochran's Q. With a p-value < 0.05 indicating the presence of heterogeneity, consequently, a random-effects IVW method would be used. Second, the weighted median analysis was applied to determine invalid instrument bias. Compared to IVW, the weighted median analysis was more robust to individual genetics with strongly outlying causal estimates (35). Third, the MR-Egger regression method was used to evaluate the potential presence of directional pleiotropy based on its intercept term, where deviation from zero denotes pleiotropy (36). Fourth, MR pleiotropy residual sum and outlier (MR-PRESSO) method was also applied to identify any potential horizontal pleiotropic outliers (37). Fifth, leave-one-SNP-out analysis was used to assess whether the observed association was affected by individual SNPs.

Additionally, we also conducted multi-variable MR (MVMR) to determine the influence of potential cardiovascular risk factors on causal estimates (38). We used publicly available summary statistics for genetic association of instruments with smoking from the GWAS and Sequencing Consortium of Alcohol and Nicotine use (39), type 2 diabetes mellitus (T2D) from the Diabetes Genetics Replication and Meta-analysis (40), BMI from the Genetic Investigation of Anthropometric (41), depression from the neuroticism GWAS (42), serum lipid levels (low-density lipoprotein cholesterol and triglycerides) from the Global Lipids Genetics Consortium (43), educational attainment from the Social Science Genetic Association Consortium (44), household income from the GWAS and eQTL studies (45), and diet (vegetable and meat consumption) from BioBank Japan (46). All statistical analyses were two-sided and p-values < 0.05 were set as the threshold for statistical significance. All analysis was performed using R software (version 3.5.4; www.r-project.org) with the MR and MRPRESSO package.

 Results

Supplementary Table 1 showed the summarized statistics of all correlated SNPs for physical activity from the published GWAS. The F statistics of all SNPs were above the threshold of 10, indicating the absence of weak instrument bias. There was some evidence of heterogeneity based on Cochran's Q (p-value < 0.05) for the CVD analysis (Supplementary Table 3); consequently, for these models, the random-effects IVW method was used. In the PhenoScanner database, we identified one (rs429358, p-value = 4.54 × 10−11) of the nine SNPs for self-reported moderate-to-vigorous physical activity associated with coronary artery disease and thus we excluded it from coronary artery disease and myocardial infarction analysis (Supplementary Table 4).

Self-Reported Physical Activity and CVDs

According to the IVW analysis results, we found evidence of an inverse association between VPA with myocardial infarction (IVW OR: 0.24, 95% CI: 0.08–0.68, p-value: 0.007), and weighted median (OR: 0.22, 95% CI: 0.06–0.78, p-value: 0.019) obtained a similar pattern of effect (Table 1; Figure 1). No directional pleiotropy was revealed by the MR-Egger intercept analysis (intercept, 1.02; 95% CI, 0.94–1.11; p-value: 0.678) (Supplementary Table 5). On the flip side, no apparent relationship was observed between VPA and heart failure, atrial fibrillation, coronary artery disease, and any ischemic stroke risk (Table 1). To explore the effect of potential cardiovascular risk factors on causal estimates, multivariable MR analysis was performed and the association of VPA with myocardial infarction was robust after adjusting for genetically predicted smoking, serum lipid levels, BMI, T2D, or depression separately (Table 2). Additionally, across all methods, we found no evidence of causal relationships of MVPA with CVDs (Table 1).

Mendelian randomization estimates between self-reported physical activity and cardiovascular diseases.

Exposure Self-reported moderate-to-vigorous physical activity Self-reported vigorous physical activity
Outcomes Method OR 95% CI p-value SNPs Method OR 95% CI p-value SNPs
Atrial fibrillation IVW 1.13 0.79–1.6 0.508 9 IVW§ 1.33 0.44–4.05 0.615 5
MR-Egger 0.41 0.08–1.99 0.268 9 MR-Egger 0.00 0–13.36 0.177 5
Weighted median 1.16 0.72–1.86 0.546 9 Weighted median 1.51 0.51–4.48 0.457 5
MR-PRESSO 1.13 0.82–1.55 0.493 9 MR-PRESSO 1.33 0.44–4.05 0.641 5
Coronary artery disease IVW 1.21 0.72–2.01 0.469 7†,‡ IVW 0.50 0.19–1.26 0.141 5
MR-Egger 0.72 0.09–5.43 0.747 7†,‡ MR-Egger 0.55 0.18–1.73 0.308 5
Weighted median 1.05 0.55–1.99 0.880 7†,‡ Weighted median 0.55 0.18–1.73 0.308 5
MR-PRESSO 1.21 0.85–1.71 0.326 7†,‡ MR-PRESSO 0.50 0.27–0.9 0.084 5
Myocardial infarction IVW 1.48 0.84–2.61 0.178 7†,‡ IVW 0.24 0.08–0.68 0.007 5
MR-Egger 0.54 0–150.86 0.829 7†,‡ MR-Egger 0.04 0–271.79 0.467 5
Weighted median 0.86 0.15–4.85 0.869 7†,‡ Weighted median 0.22 0.06–0.78 0.019 5
MR-PRESSO 1.48 0.93–2.35 0.151 7†,‡ MR-PRESSO 0.24 0.12–0.47 0.014 5
Heart failure IVW 0.89 0.6–1.31 0.543 9 IVW 1.57 0.71–3.45 0.263 5
MR-Egger 0.47 0.09–2.38 0.360 9 MR-Egger 0.00 0.00–0.25 0.017 5
Weighted median 1.10 0.64–1.87 0.738 9 Weighted median 2.24 0.74–6.77 0.151 5
MR-PRESSO 0.89 0.62–1.27 0.532 9 MR-PRESSO 1.57 0.51–4.78 0.473 5
Ischemic stroke IVW 1.14 0.69–1.89 0.606 8† IVW 0.97 0.36–2.65 0.955 5
MR-Egger 0.61 0.08–4.42 0.622 8† MR-Egger 0.09 0.00–375.81 0.565 5
Weighted median 1.01 0.54–1.9 0.977 8† Weighted median 0.79 0.23–2.78 0.718 5
MR-PRESSO 1.14 0.79–1.66 0.509 8† MR-PRESSO 0.97 0.45–2.09 0.945 5
Large artery stroke IVW 0.59 0.17–2.01 0.395 8† IVW 4.47 0.37–54.03 0.238 5
MR-Egger 0.07 0–34.44 0.396 8† MR-Egger 32.94 0.00–10∧15 0.827 5
Weighted median 0.47 0.09–2.51 0.374 8† Weighted median 1.45 0.04–55.54 0.841 5
MR-PRESSO 0.59 0.12–2.79 0.523 8† MR-PRESSO 4.47 0.18–113.99 0.416 5
Cardioembolic stroke IVW 1.68 0.64–4.41 0.289 8† IVW 1.06 0.16–7.05 0.956 5
MR-Egger 7.30 0.1–531.9 0.363 8† MR-Egger 0.00 0.00–61073.42 0.372 5
Weighted median 1.97 0.54–7.19 0.303 8† Weighted median 0.53 0.05–5.94 0.605 5
MR-PRESSO 1.68 0.58–4.92 0.373 8† MR-PRESSO 1.06 0.1–10.86 0.966 5
Small vessel stroke IVW 0.95 0.3–3.01 0.929 8† IVW 4.16 0.41–42.13 0.228 5
MR-Egger 0.60 0.01–60.23 0.828 8† MR-Egger 0.00 0.00–18910.03 0.333 5
Weighted median 0.85 0.19–3.87 0.832 8† Weighted median 1.85 0.09–37.68 0.690 5
MR-PRESSO 0.95 0.32–2.82 0.927 8† MR-PRESSO 4.16 0.51–33.61 0.252 5

OR, odds ratio; CI confidence intervals; IVW, inverse-variance-weighted method; MR-PRESSO, MR pleiotropy residual sum and outlier method; SNPs, single-nucleotide polymorphisms.

The estimates were evaluated from a random-effects IVW method due to the presence of heterogeneity based on Cochran's Q.

rs3094622 was excluded because it was unavailable in the published GWASs and no good proxies (r2 > 0.8) were found.

rs429358 was excluded since it as associated with coronary artery disease in the PhenoScanner database.

A scatter plot for the relationship of self-reported vigorous physical activity with myocardial infarction. A scatter plot of single-nucleotide polymorphism (SNP) effects regarding the associations of genetically predicted self-reported vigorous physical activity (VPA) with myocardial infarction (MI), with the slope of each line corresponding to estimated Mendelian Randomization (MR) effect per method. Circles indicate marginal genetic associations with self-reported vigorous physical activity and risk of myocardial infarction for each variant. Error bars indicate 95% CIs.

Multivariable Mendelian randomization associations of vigorous physical activity with myocardial infarction risk adjusting for cardiovascular risk factors.

Model OR 95% CI p
Unadjusted model 0.24 0.08–0.68 0.007
Adjusted for BMI 0.24 0.08–0.68 0.007
Adjusted for LDL-C 0.21 0.07–0.62 0.005
Adjusted for TC 0.21 0.07–0.61 0.004
Adjusted for depression 0.26 0.09–0.79 0.017
Adjusted for T2D 0.26 0.09–0.75 0.013
Adjusted for smoking 0.21 0.07–0.61 0.004
Adjusted for vegetable consumption 0.22 0.06–0.68 0.019
Adjusted for meat consumption 0.26 0.09–0.74 0.012
Adjusted for household income 0.23 0.08–0.68 0.008
Adjusted for educational attainment 0.22 0.07–0.69 0.009

Results are scaled per genetically predicted standard deviation increase of liability to vigorous physical activity. BMI, body mass index; LDL-C, low-density lipoprotein cholesterol; TC, total cholesterol; SBP, systolic blood pressure; T2D, type 2 diabetes mellitus.

For self-reported physical activity (including MVPA and VPA)-related SNPs, MR-PRESSO did not detect any potential outliers (Supplementary Table 6), and results of the leave-one-SNP-out analysis (Supplementary Tables 7, 8) suggested that the association between VPA and CVDs was not affected by single SNPs.

 Accelerometer-Based Physical Activity and CVDs

In contrast, there was no significant association to be found between accelerometer-based physical activity and CVDs (Table 3). These findings were confirmed using weight median, MR-Egger regression method, and leave-one-SNP-out analysis (Supplement Tables 5, 9, 10). MR-PRESSO identified several outliers (Supplementary Table 6), but similar MR estimates were observed after removal of these outliers (Supplementary Tables 11, 12).

Mendelian randomization estimates between accelerometer-measured physical activity and cardiovascular diseases.

Exposure Overall acceleration average Fraction of accelerations > 425 milli-gravities
Outcomes Method OR 95% CI p-value SNPs Method OR 95% CI p-value SNPs
Atrial fibrillation IVW§ 1.01 0.96–1.06 0.616 8 IVW§ 1.33 0.89–1.99 0.158 8
MR-Egger 1.09 0.87–1.37 0.435 8 MR-Egger 20.26 0.02–104 0.387 8
Weighted median 1.00 0.96–1.04 0.947 8 Weighted median 1.09 0.78–1.53 0.595 8
MR-PRESSO 1.01 0.96–1.06 0.631 8 MR-PRESSO 1.33 0.89–1.99 0.201 8
MR-PRESSO (outlier-corrected) 1.00 0.97–1.03 0.871 6# MR-PRESSO (outlier-corrected) 1.18 0.84–1.64 0.376 7#
Coronary artery disease IVW§ 1.02 0.97–1.08 0.396 8 IVW§ 0.89 0.52–1.51 0.660 8
MR-Egger 1.19 0.95–1.48 0.134 8 MR-Egger 336.74 0.08–106 0.174 8
Weighted median 1.00 0.95–1.04 0.876 8 Weighted median 0.86 0.54–1.39 0.542 8
MR-PRESSO 1.02 0.97–1.08 0.424 8 MR-PRESSO 0.89 0.52–1.51 0.673 8
MR-PRESSO (outlier-corrected) 1.01 0.96–1.06 0.809 7# MR-PRESSO (outlier-corrected) 0.77 0.48–1.24 0.329 7#
Myocardial infarction IVW§ 1.01 0.96–1.07 0.624 8 IVW§ 0.88 0.51–1.52 0.647 8
MR-Egger 1.14 0.88–1.48 0.326 8 MR-Egger 720.62 0.21–2 × 106 0.114 8
Weighted median 0.99 0.94–1.04 0.792 8 Weighted median 0.67 0.41–1.1 0.113 8
MR-presso 1.01 0.96–1.07 0.639 8 MR-PRESSO 0.88 0.51–1.52 0.661 8
MR-presso (outlier-corrected) 0.99 0.95–1.04 0.759 7# MR-PRESSO (outlier-corrected) 0.74 0.48–1.15 0.232 7#
Heart failure IVW 0.98 0.95–1.01 0.132 8 IVW 1.03 0.81–1.32 0.785 8
MR-Egger 0.99 0.85–1.17 0.945 8 MR-Egger 20.06 0.36–1,118.6 0.144 8
Weighted median 0.98 0.95–1.02 0.413 8 Weighted median 1.00 0.72–1.38 0.985 8
MR-PRESSO 0.98 0.95–1.01 0.293 8 MR-presso 1.03 0.82–1.31 0.783 8
MR-PRESSO (outlier-corrected) NA NA NA NA MR-PRESSO (outlier-corrected) NA NA NA NA
Ischemic stroke IVW§ 0.99 0.93–1.05 0.712 8 IVW 0.82 0.6–1.13 0.228 8
MR-Egger 1.13 0.87–1.46 0.349 8 Simple median 0.85 0.55–1.3 0.445 8
Weighted median 0.99 0.94–1.04 0.701 8 MR-Egger 2.19 0–1,555.14 0.815 8
MR-PRESSO 0.99 0.93–1.05 0.723 8 MR-PRESSO 0.82 0.57–1.2 0.345 8
MR-PRESSO (outlier-corrected) 0.97 0.92–1.03 0.378 7# MR-PRESSO (outlier-corrected) NA NA NA NA
Large artery stroke IVW 0.94 0.87–1.02 0.165 8 IVW 0.64 0.29–1.39 0.257 8
MR-Egger 0.90 0.62–1.31 0.574 8 Simple median 0.81 0.3–2.15 0.672 8
Weighted median 0.94 0.85–1.05 0.274 8 MR-Egger 0.02 0–6,730.82 0.538 8
MR-PRESSO 0.94 0.9–0.99 0.054 8 MR-PRESSO 0.64 0.38–1.07 0.131 8
MR-PRESSO (outlier-corrected) NA NA NA NA MR-PRESSO (outlier-corrected) NA NA NA NA
Cardioembolic stroke IVW 0.98 0.92–1.05 0.567 8 IVW 0.62 0.34–1.15 0.130 8
MR-Egger 1.09 0.81–1.46 0.585 8 MR-Egger 1.14 0–24,321.08 0.979 8
Weighted median 0.99 0.91–1.08 0.882 8 Weighted median 0.53 0.25–1.15 0.108 8
MR-PRESSO 0.98 0.93–1.04 0.539 8 MR-PRESSO 0.62 0.41–0.95 0.066 8
MR-PRESSO (outlier-corrected) NA NA NA NA MR-PRESSO (outlier-corrected) NA NA NA NA
Small vessel stroke IVW§ 0.96 0.82–1.11 0.580 8 IVW 1.25 0.6–2.62 0.548 8
MR-Egger 1.02 0.49–2.13 0.961 8 MR-Egger 10∧5 0.14–1011 0.094 8
Weighted median 0.90 0.8–1.01 0.075 8 Weighted median 0.68 0.26–1.8 0.440 8
MR-PRESSO 0.96 0.82–1.11 0.597 8 MR-PRESSO 1.25 0.49–3.2 0.651 8
MR-PRESSO (outlier-corrected) 0.90 0.79–1.02 0.154 7# MR-PRESSO (outlier-corrected) NA NA NA NA

OR, odds ratio; CI, confidence intervals; IVW, inverse-variance-weighted method; MR-PRESSO, MR pleiotropy residual sum and outlier method; SNPs, single-nucleotide polymorphisms.

The estimates were evaluated from a random-effects IVW method due to the presence of heterogeneity based on Cochran's Q.

MR-PRESSO outlier removed. Detailed information of outliers is shown in Supplementary Table 6.

 Discussion

In this study, we applied MR to evaluate the causal relationship between physical activityand CVDs. Although physical activity might have a potential influence on the progress of CVDs, as shown in several previous studies (7–9), our study only suggests the causal relationship between self-reported vigorous physical activity and myocardial infarction.

Our findings were in line with a meta-analysis of 33 prospective cohort studies showing that physical activity was associated with myocardial infarction (7). The higher levels of physical activity, the lower risk of myocardial infarction (47). A case–control study included 18,225 individuals showing that vigorous activity is associated with a 22% lower risk of myocardial infarction (48). In this MR study, the causal relationship appeared to be motivated by self-reported vigorous physical activity, but not moderate-to-vigorous physical activity, which meant that activity levels might be more important than total activity time (49). Previous empirical studies had already suggested that several factors would affect the level of physical activity, such as education and depression (50, 51). Additionally, to obtain more robust results, we further applied MVMR analysis and demonstrated that the causal relationship would not be influenced by several CVD risk factors, including smoking, serum lipid levels, BMI, T2D, and depression. Nevertheless, there was no significant association between accelerometer-based physical activity and myocardial infarction. The discrepancy in findings between self-reported and accelerometer-based physical activity could be caused by some reasons. Firstly, compared to 377,234 participants of self-reported activity, only ~91,000 participants wore an accelerometer. Different sample sizes may influence the findings (52). Secondly, individuals who were more likely to accept wearing an accelerometer were women, aged 55–74 years, with higher socioeconomic status, and better physical health status. This may lead to selection bias and influence the results. Finally, the accelerometers only recorded the physical activity data of participants for 7 days and therefore could not effectively report the long-term activity volume.

Previous meta-analysis of observational studies reported inconsistent results on the causal relationship between physical activity (assessed using self-reported and accelerometer-measured activity) and atrial fibrillation. A large meta-analysis of 15 observational studies, including 1,464,539 participants, indicated a 6% relative reduction (hazard ratio: 0.94, 95% CI: 0.90–0.97) in atrial fibrillation risk with the guideline-recommended level of physical activity (53). In contrast, Kwok et al. conducted a meta-analysis and suggested that there is no significant decrease in AF risk with either moderate or vigorous physical activity (54). Our study may be the first to evaluate the causal association of physical activity with atrial fibrillation by using MR analysis. We found that self-reported and objective accelerometer-based physical activity was not causally associated with atrial fibrillation. These results suggested that conclusions from some previous observational studies might be false positives because of confounding and reverse causation.

Additionally, we identified no significant relationship of self-reported and objective accelerometer-based physical activity with heart failure, which supported the finding from a prospective study including 37,803 participants showing that physical activity was not independently associated with reduced risk of heart failure and indicated that the association of physical activity with incident heart failure is mediated through other risk factors, such as smoking, BMI, total cholesterol, and glucose (55). In this study, we further showed that physical activity does not affect the risk of ischemic stroke and its subtypes, which is consistent with a previous study showing that excessive activities might not be beneficial for hemorrhagic stroke (56).

Our study included several notable strengths. First, for the first time, we explored the causal relationship of self-reported and objective accelerometer-based physical activity with atrial fibrillation, heart failure, ischemic stroke, and its subtypes using MR analysis. Second, to get more robust results, we also conducted weighted median analysis, MR-Egger regression, MR-PRESSO, and leave-one-SNP-out analysis as sensitivity analysis to ensure the consistency of causal relationships. Furthermore, we also used genetically predicted objectively measured physical activity, which is more heritable than self-reported physical activity. Third, we performed multivariable MR to explore the influence of potential cardiovascular risk factors on causal estimates between VPA and myocardial infarct, and the results were consistent.

 Limitation

However, our MR study also has some limitations. First, despite the consistent results in the sensitivity analysis, we still could not completely rule out the influence of potential pleiotropy, whose physical activity might affect CVDs by other unknown causal pathways. Second, participants of physical activity in this study were all of European descent, aged 40 to 70 years, which might limit the generalizability of our findings. Further studies are needed to verify our findings in other age groups and non-European descents. Third, the instrumental variables for physical activity explained a small fraction of the phenotypic variability, which led to some results not reaching a statistical power of 80%. Consequently, the 95% confidence intervals for our MR-Egger analysis were wide in some outcomes. Thus, additional and larger replication studies were required to more robustly identify PA-associated loci. Fourth, due to the assumption of the linear relationship between exposure and outcome in MR design, a potential non-linear association of physical activity with CVDs could not be evaluated and can still not exclude the non-linear effect of physical activity.

 Conclusion

The present study supports a protective role of self-reported vigorous physical activity on the initiation of myocardial infarction and highlights the importance of activity levels of physical activity. Based on these data, the promotion of vigorous physical activity is most likely an effective strategy in the primary prevention of myocardial infarction compared with moderate-to-vigorous physical activity. Through identifying the causal relationship between physical activity and CVDs, we have laid the groundwork for future investigations on the comprehensive clinical approach in preventing the onset of cardiovascular diseases, which can be achieved through lifestyle interventions in addition to medication.

 Data Availability Statement

The summary statistics of physical activity (Study accession ID: GCST006079, GCST006097, GCST006098, and GCST006099), coronary artery disease (Study accession ID: GCST003116), myocardial infarction (Study accession ID: GCST003117), atrial fibrillation (Study accession ID: GCST006061), heart failure (Study accession ID: GCST009541), and stroke (Study accession ID: GCST005840, GCST005841, GCST005842, and GCST005843) are available at https://www.ebi.ac.uk/gwas/.

 Author Contributions

CZ: designed the study and wrote the manuscript. JZ and MC: contributed to the data analysis and data interpretation. YL: contributed to the revision of the manuscript. All authors contributed to the article and approved the submitted version.

 Funding

This work was supported by a grant from the Department of Science and Technology, Zhejiang Province (Grant no. LGF19H020011).

 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.

We thank all of the investigators for sharing summary-level data on GWAS for physical activity, coronary artery disease, myocardial infarction, heart failure, atrial fibrillation, and ischemic stroke and its subtypes.

 Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcvm.2021.722154/full#supplementary-material

 References Lennon RP Claussen KA Kuersteiner KA. State of the heart an overview of the disease burden of cardiovascular disease from an epidemiologic perspective. Primary Care. (2018) 45:1–15. 10.1016/j.pop.2017.11.00129406937 Rapsomaniki E Timmis A George J Pujades-Rodriguez M Shah AD Denaxas S . Blood pressure and incidence of twelve cardiovascular diseases: lifetime risks, healthy life-years lost, and age-specific associations in 1.25 million people. Lancet. (2014) 383:1899–911. 10.1016/S0140-6736(14)60685-124881994 Wilson PWF D'Agostino RB Sullivan L Parise H Kannel WB. Overweight and obesity as determinants of cardiovascular risk - the framingham experience. Arch Int Med. (2002) 162:1867–72. 10.1001/archinte.162.16.186712196085 Huxley RR Woodward M. Cigarette smoking as a risk factor for coronary heart disease in women compared with men: a systematic review and meta-analysis of prospective cohort studies. Lancet. (2011) 378:1297–305. 10.1016/S0140-6736(11)60781-221839503 Anderson L Brown JPR Clark AM Dalal H Rossau HK Bridges C . Patient education in the management of coronary heart disease. Cochrane Database Syst Rev. (2017) 6:CD008895. 10.1002/14651858.CD008895.pub328658719 Mozaffarian D. Dietary and policy priorities for cardiovascular disease, diabetes, and obesity a comprehensive review. Circulation. (2016) 133:187–225. 10.1161/CIRCULATIONAHA.115.01858526746178 Sattelmair J Pertman J Ding EL Kohl HW III Haskell W . Dose response between physical activity and risk of coronary heart disease a meta-analysis. Circulation. (2011) 124:789–U4. 10.1161/CIRCULATIONAHA.110.01071021810663 Reddigan JI Ardern CI Riddell MC Kuk JL. Relation of physical activity to cardiovascular disease mortality and the influence of cardiometabolic risk factors. Am J Cardiol. (2011) 108:1426–31. 10.1016/j.amjcard.2011.07.00521855834 Claas SA Arnett DK. The role of healthy lifestyle in the primordial prevention of cardiovascular disease. Curr Cardiol Rep. (2016) 18:56. 10.1007/s11886-016-0728-727142061 Piepoli MF Hoes AW Agewall S Albus C Brotons C Catapano AL . 2016 European guidelines on cardiovascular disease prevention in clinical practice. Eur Heart J. (2016) 37:2315–81. 10.1093/eurheartj/ehw10628290077 Arija V Villalobos F Pedret R Vinuesa A Timon M Basora T . Effectiveness of a physical activity program on cardiovascular disease risk in adult primary health-care users: the “Pas-a-Pas” community intervention trial. Bmc Public Health. (2017) 17:576. 10.1186/s12889-017-4485-328619115 Tiessen AH Smit AJ Broer J Groenier KH van der Meer K. Randomized controlled trial on cardiovascular risk management by practice nurses supported by self-monitoring in primary care. Bmc Family Practice. (2012) 13:90. 10.1186/1471-2296-13-9022947269 Garcia-Ortiz L Grandes G Sanchez-Perez A Montoya I Iglesias-Valiente JA Recio-Rodriguez JI . Effect on cardiovascular risk of an intervention by family physicians to promote physical exercise among sedentary individuals. Revista Espanol Cardiol. (2010) 63:1244–52. 10.1016/S1885-5857(10)70249-821070720 Elley CR Kerse N Arroll B Robinson E. Effectiveness of counselling patients on physical activity in general practice: cluster randomised controlled trial. Brit Med J. (2003) 326:793–6. 10.1136/bmj.326.7393.79312689976 Holmes MV Ala-Korpela M Smith GD. Mendelian randomization in cardiometabolic disease: challenges in evaluating causality. Nat Rev Cardiol. (2017) 14:577–90. 10.1038/nrcardio.2017.7828569269 van Oort S Beulens JWJ van Ballegooijen AJ Burgess S Larsson SC. Cardiovascular risk factors and lifestyle behaviours in relation to longevity: a mendelian randomization study. J Int Med. (2021) 289:232–43. 10.1111/joim.1319633107078 Gill D Zuber V Dawson J Pearson-Stuttard J Carter AR Sanderson E . Risk factors mediating the effect of body mass index and waist-to-hip ratio on cardiovascular outcomes: mendelian randomization analysis. Int J Obesity. (2021) 45:1428–38. 10.1038/s41366-021-00807-434002035 Rosoff DB Smith GD Mehta N Clarke T-K Lohoff FW. Evaluating the relationship between alcohol consumption, tobacco use, and cardiovascular disease: a multivariable mendelian randomization study. PLoS Med. (2020) 17:e1003410. 10.1371/journal.pmed.100341033275596 Carlsson S Andersson T Lichtenstein P Michaëlsson K Ahlbom A. Genetic effects on physical activity: results from the swedish twin registry. Med Sci Sports Exerc. (2006) 38:1396–401. 10.1249/01.mss.0000228941.17034.c116888451 Joosen AM Gielen M Vlietinck R Westerterp KR. Genetic analysis of physical activity in twins. Am J Clin Nutr. (2005) 82:1253–9. 10.1093/ajcn/82.6.125316332658 Stubbe JH Boomsma DI De Geus EJC. Sports participation during adolescence: a shift from environmental to genetic factors. Med Sci Sports Exerc. (2005) 37:563–70. 10.1249/01.MSS.0000158181.75442.8B15809553 Klimentidis YC Raichlen DA Bea J Garcia DO Wineinger NE Mandarino LJ . Genome-wide association study of habitual physical activity in over 377,000 UK Biobank participants identifies multiple variants including CADM2 and APOE. Int J Obesity. (2018) 42:1161–76. 10.1038/s41366-018-0120-329899525 Craig CL Marshall AL Sjostrom M Bauman AE Booth ML Ainsworth BE . International physical activity questionnaire: 12-country reliability and validity. Med Sci Sports Exercise. (2003) 35:1381–95. 10.1249/01.MSS.0000078924.61453.FB12900695 Doherty A Jackson D Hammerla N Plotz T Olivier P Granat MH . Large scale population assessment of physical activity using wrist worn accelerometers: the UK biobank study. PLoS ONE. (2017) 12:e0169649. 10.1371/journal.pone.016964928146576 Gage SH Jones HJ Burgess S Bowden J Smith GD Zammit S . Assessing causality in associations between cannabis use and schizophrenia risk: a two-sample mendelian randomization study. Psychol Med. (2017) 47:971–80. 10.1017/S003329171600317227928975 Choi KW Chen C-Y Stein MB Klimentidis YC Wang M-J Koenen KC . Assessment of bidirectional relationships between physical activity and depression among adults a 2-sample mendelian randomization study. JAMA Psychiatry. (2019) 76:399–408. 10.1001/jamapsychiatry.2018.417530673066 Nikpay M Goel A Won H-H Hall LM Willenborg C Kanoni S . A comprehensive 1000 genomes-based genome-wide association meta-analysis of coronary artery disease. Nat Gen. (2015) 47:1121–30. 10.1038/ng.339626343387 Shah S Henry A Roselli C Lin H Sveinbjornsson G Fatemifar G . Genome-wide association and mendelian randomisation analysis provide insights into the pathogenesis of heart failure. Nature Commun. (2020) 11:163. 10.1038/s41467-019-13690-531919418 Roselli C Chaffin MD Weng L-C Aeschbacher S Ahlberg G Albert CM . Multi-ethnic genome-wide association study for atrial fibrillation. Nat Gen. (2018) 50:1225–33. 10.1038/s41588-018-0133-929892015 Malik R Chauhan G Traylor M Sargurupremraj M Okada Y Mishra A . Multiancestry genome-wide association study of 520,000 subjects identifies 32 loci associated with stroke and stroke subtypes. Nat Gen. (2018) 50:524–37. 10.1038/s41588-018-0058-331160810 Adams HP Bendixen BH Kappelle LJ Biller J Love BB Gordon DL . Classification of subtype of acute ischemic stroke - definitions for use in a multicenter clinical-trial. Stroke. (1993) 24:35–41. 10.1161/01.STR.24.1.357678184 Brion M-JA Shakhbazov K Visscher PM. Calculating statistical power in mendelian randomization studies. Int J Epidemiol. (2013) 42:1497–501. 10.1093/ije/dyt17924159078 Thompson JR Minelli C Del Greco MF. Mendelian randomization using public data from genetic consortia. Int J Biostat. (2016) 12:74. 10.1515/ijb-2015-007427092657 Kamat MA Blackshaw JA Young R Surendran P Burgess S Danesh J . PhenoScanner V2: an expanded tool for searching human genotype-phenotype associations. Bioinformatics. (2019) 35:4851–3. 10.1093/bioinformatics/btz46931233103 Geng T Smith CE Li C Huang T. Childhood BMI and adult type 2 diabetes, coronary artery diseases, chronic kidney disease, and cardiometabolic traits: a mendelian randomization analysis. Diabetes Care. (2018) 41:1089–96. 10.2337/dc17-214129483184 Bowden J Smith GD Burgess S. Mendelian randomization with invalid instruments: effect estimation and bias detection through egger regression. Int J Epidemiol. (2015) 44:512–25. 10.1093/ije/dyv08026050253 Verbanck M Chen C-Y Neale B Do R. Detection of widespread horizontal pleiotropy in causal relationships inferred from mendelian randomization between complex traits and diseases. Nat Gen. (2018) 50:693–8. 10.1038/s41588-018-0099-729967445 Burgess S Thompson SG. Multivariable mendelian randomization: the use of pleiotropic genetic variants to estimate causal effects. Am J Epidemiol. (2015) 181:251–60. 10.1093/aje/kwu28325660081 Liu M Jiang Y Wedow R Li Y Brazel DM Chen F . Association studies of up to 1.2 million individuals yield new insights into the genetic etiology of tobacco and alcohol use. Nat Gen. (2019) 51:237–44. 10.1038/s41588-018-0307-530643251 Morris AP Voight BF Teslovich TM Ferreira T Segre AV Steinthorsdottir V . Large-scale association analysis provides insights into the genetic architecture and pathophysiology of type 2 diabetes. Nat Gen. (2012) 44:981–90. 10.1038/ng.238323250902 Locke AE Kahali B Berndt SI Justice AE Pers TH Day FR . Genetic studies of body mass index yield new insights for obesity biology. Nature. (2015) 518:197–206. 10.1038/nature1417725673413 Nagel M Jansen PR Stringer S Watanabe K de Leeuw CA Bryois J . Meta-analysis of genome-wide association studies for neuroticism in 449,484 individuals identifies novel genetic loci and pathways. Nat Gen. (2018) 50:920–7. 10.1038/s41588-018-0151-729942085 Willer CJ Schmidt EM Sengupta S Peloso GM Gustafsson S Kanoni S . Discovery and refinement of loci associated with lipid levels. Nat Gen. (2013) 45:1274–83. 10.1038/ng.279724097068 Okbay A Beauchamp JP Fontana MA Lee JJ Pers TH Rietveld CA . Genome-wide association study identifies 74 loci associated with educational attainment. Nature. (2016) 533:539–42. 10.1038/nature1767128197077 Hill WD Davies NM Ritchie SJ Skene NG Bryois J Bell S . Genome-wide analysis identifies molecular systems and 149 genetic loci associated with income. Nat Commun. (2019) 10:5741. 10.1038/s41467-019-13585-531844048 Matoba N Akiyama M Ishigaki K Kanai M Takahashi A Momozawa Y . GWAS of 165,084 Japanese individuals identified nine loci associated with dietary habits. Nat Hum Behav. (2020) 4:308–16. 10.1038/s41562-019-0805-131959922 Zhao S Zhang Z Long Q Ma Y Lian X Yang Y . Association between time of day of sports-related physical activity and the onset of acute myocardial infarction in a chinese population. PLoS ONE. (2016) 11:e0146472. 10.1371/journal.pone.014647226752185 Chomistek AK Chiuve SE Jensen MK Cook NR Rimm EB. Vigorous physical activity, mediating biomarkers, and risk of myocardial infarction. Med Sci Sports Exer. (2011) 43:1884–90. 10.1249/MSS.0b013e31821b4d0a21448079 Winzer EB Woitek F Linke A. Physical activity in the prevention and treatment of coronary artery disease. J Am Heart Assoc. (2018) 7:e007725. 10.1161/JAHA.117.00772529437600 Kari JT Viinikainen J Böckerman P Tammelin TH Pitkänen N Lehtimäki T . Education leads to a more physically active lifestyle: evidence based on mendelian randomization. Scand J Med Sci Sports. (2020) 30:1194–204. 10.1111/sms.1365332176397 Gudmundsson P Lindwall M Gustafson DR Östling S Hällström T Waern M . Longitudinal associations between physical activity and depression scores in Swedish women followed 32 years. Acta Psychiatr Scand. (2015) 132:451–8. 10.1111/acps.1241925865488 Yeung SLA Jiang C Cheng KK Cowling BJ Liu B Zhang W . Moderate alcohol use and cardiovascular disease from mendelian randomization. PLoS ONE. (2013) 8:e68054. 10.1371/journal.pone.006805423874492 Mishima RS Verdicchio CV Noubiap JJ Ariyaratnam JP Gallagher C Jones D . Self-reported physical activity and atrial fibrillation risk: a systematic review and meta-analysis. Heart Rhythm. (2021) 18:520–8. 10.1016/j.hrthm.2020.12.01733348059 Kwok CS Anderson SG Myint PK Mamas MA Loke YK. Physical activity and incidence of atrial fibrillation: a systematic review and meta-analysis. Int J Cardiol. (2014) 177:467–76. 10.1016/j.ijcard.2014.09.10425443248 Uijl A Koudstaal S Vaartjes I Boer JMA Verschuren WMM van der Schouw YT . Risk for heart failure the opportunity for prevention with the American heart association's life's simple 7. Jacc Heart Fail. (2019) 7:637–47. 10.1016/j.jchf.2019.03.00931302040 Kubota Y Iso H Yamagishi K Sawada N Tsugane S Grp JS. Daily total physical activity and incident stroke the Japan public health center-based prospective study. Stroke. (2017) 48:1730–6. 10.1161/STROKEAHA.117.01756028584003