Leisure-Time Physical Activity and Cardiovascular Disease Risk Among Hypertensive Patients: A Longitudinal Cohort Study

Introduction

The 2019 American College of Cardiology (ACC)/American Heart Association (AHA) Guideline on the Primary Prevention of Cardiovascular Disease (1) and the 2018 Physical Activity Guidelines for Americans (2) both recommended at least 150 min/week of moderate–intensity physical activity (PA) or 75 min/week of vigorous-intensity PA to stay healthy. Many previous studies have suggested that the PA could reduce the risk of cardiovascular disease (CVD) (3–5) and weaken its risk factors, such as diabetes, hypertension, and obesity (6–8). However, few studies have assessed the association between leisure-time PA (LTPA) and CVD among hypertensive patients in a longitudinal cohort with time-dependent confounders.

Previous studies used a single measurement value to evaluate LTPA over the study period and assumed it was constant, such as baseline LTPA or the average of several LTPA measurements over a certain period. However, LTPA is not a static lifestyle but a long-term dynamic process, which may change over time and be affected by many factors. For example, the blood pressure (BP) of hypertensive patients varies over time during the study period, and BP is not only associated with CVD risk but also affects LTPA. Meanwhile, LTPA also has an impact on BP at the next time point (9–11). In the above case, we regard BP as a time-dependent confounder. A diagram of the time-dependent confounder is shown in Supplementary Figure 1. In the presence of time-dependent confounders, conventional analysis (such as hierarchical analysis, regression analysis) without appropriate accounting for dynamic interaction between time-varying LTPA and time-dependent confounders is prone to biased estimates (12). Therefore, we should adjust for time-dependent confounders when assessing the association between LTPA and CVD risk. Robin et al. proposed the marginal structure model (MSM), which was developed for causal inference of time-varying treatments in observational longitudinal studies with the presence of time-dependent confounders (12). The MSM uses the inverse probability of treatment weights (IPTWs) to create a virtual population and eliminates the influence of time-dependent confounders on the association estimation. To our knowledge, no previous study has used this model to estimate the association between physical activity and CVD in a longitudinal cohort.

The primary objective of this study was to estimate the effect of LTPA on CVD risk in a longitudinal cohort of hypertensive patients and to demonstrate the merits of MSM by comparing it with Cox proportional hazards model. The secondary goal of this study is to investigate the piecewise-linear relationship between LTPA and CVD risk and to explore whether there are break-points of the LTPA where the increase rate in health benefits changes significantly.

Methods

Study Cohort

The data for this study were collected from a community-based hypertension longitudinal management cohort in Shandong Province. Patients diagnosed with hypertension between June 2010 and October 2017 were included in the cohort. The patients in this cohort were followed up face to face four times every year. LTPA and other covariates were recorded at each follow-up. We included patients who met the following criteria: (1) no previous history of CVD events; (2) participants without abnormal BP [systolic blood pressure (SBP) ≥200 mmHg or SBP <80 mgHg or diastolic blood pressure (DBP) ≥140 mmHg or DBP <40 mmHg] at any follow-up; (3) no covariates missing; (4) participants with age ranging from 30 to 80 years; and (5) at least three follow-ups. Supplementary Figure 2 presented the flowchart of enrollment. The study protocol was approved by the ethics committee of the School of Public Health, Shandong University.

Study Outcome

The primary outcome was CVD events, which included coronary heart disease (CHD) and cerebrovascular diseases (CBD). The International Classification of Diseases, 10th Revision (ICD-10) clinical codes were used to identify outcomes. The ICD-10 codes for CHD included I20, I21, I22, I23, I24, and I25. The ICD-10 codes for CBD included I60, I61, I62, I63, I64, I65, I66, and I67. To determine whether participants suffered from CVD events during the study period, the participants' unique identification card numbers were used to matched CVD events in Shandong medical insurance database, hospital electronic medical record data, and Shandong death registration data for every participant. The earliest record date of CHD and CBD in the above database was defined as the diagnosis date of CHD and CBD, respectively. CVD diagnosis date was defined as the earlier date of CHD and CBD diagnosis. Participants who did not matched with CVD events in the three databases were regarded as non-CVD.

LTPA

In this study, LTPA refers to activities that are taken consciously to strengthen physical fitness. Light to moderate intensity PA includes brisk walking, jogging/running, square dancing, bicycling, and table tennis. In every follow-up visit, the doctors of the community health service center obtained the LTPA frequency (times/week) and duration (minutes/time) of patients in the last 3 months through inquiries. We multiply the frequency (times/week) and duration (minutes/time) to obtain the weekly LTPA time (minutes/week), which was used to measure LTPA levels of hypertensive patients. According to the LTPA levels (150, 300 min/week) recommended by PA guidelines (9–11), the total participants have been categorized into four groups: 0, 1–150, 151–300, and >300 min/week. This study used two forms of LTPA measurement [per increase of 60 min/week; 4 LTPA categories (0, 1–150, 151–300, >300 min/week)] to estimate the association between LTPA and CVD risk.

Covariates

The covariates include age, sex, SBP, DBP, body mass index (BMI), smoking, drinking, diabetes, hyperlipidemia, and antihypertensive medication (yes/no). Height and weight were measured when patients were wearing light clothing without shoes. BMI is calculated as weight (kg) divided by the square of height (m). After 5 min of rest, the right upper arm blood pressure, including SBP and DBP, was measured in a sitting position by doctors. Smoking and drinking were defined as smoking at least one cigarette a day and drinking alcohol at least once a week on average, respectively. The diagnosis records of patients in Shandong Provincial Medical insurance database and electronic medical records of hospitals were used to determine whether the patients had diabetes and hyperlipidemia. Antihypertensive medication use is defined as the cumulative use of any antihypertensive drug for more than 3 months. Antihypertensive medication included alpha-blocker, angiotensin II receptor blocker, angiotensin-converting enzyme inhibitor, beta-blocker, calcium-channel blocker, diuretics, or any combination of the above classes of medications. A list of antihypertensive medications and Anatomical Therapeutic Chemical is shown in Supplementary Table 1.

Statistical Analysis

If the continuous variables follow normality, they were presented as mean and standard deviations (SDs), and the analysis of variance (ANOVA) was used for comparison between groups; otherwise, they were presented as medians with interquartile range (IQR), and the Kruskal–Wallis H test was applied for intergroup comparison. Categorical variables were summarized as frequencies and proportions; the comparison between groups was carried out using the chi-square test.

A multivariable Cox proportional hazards regression model was used to examine the association between baseline LTPA and CVD risk with adjustment for baseline covariates including age, gender, BMI, smoking, drinking, diabetes, hyperlipidemia, and antihypertensive medication. The marginal structural Cox model (MSM-Cox) was used to examine the association of LTPA with CVD risk with adjustment for time-dependent confounders. First, the model estimates the inverse probability of receiving a specific intervention given the baseline and time-dependent confounders. Then, the IPTW is used to create a virtual population in which the confounding effects are eliminated by weighting and the effect of treatment remains unchanged. In this pseudo-population, Cox proportional hazards model was applied to estimate the association of treatment and outcomes. More details of the MSM-Cox refers to this artile. (13). LTPA was included in the model in the form of continuous variables (per increase of 60 min/week) and categorical variables (0, 1–150, 151–300, >300 min/week). The corresponding hazard ratios (HRs) estimated by the MSM should be interpreted as the risk of CVD for someone whose LTPA increase by 60 min/week and the risk of CVD for someone whose LTPA remained in a particular LTPA category over the study period, respectively. We also conducted sensitivity analyses to investigate the association in different subpopulations stratified by age (≤60, 60–70, ≥70 years), gender (men, women), BMI (<24, 24–27, and >27 kg/m²), diabetes mellitus (yes/no), hyperlipidemia (yes/no), smoking (yes/no), drinking (yes/no), antihypertensive medication (yes/no), and test for interactions of LTPA with covariates by adding interactions terms to the model.

In addition, this study conducted two analyses to explore the dose–response relationship and piecewise-linear relationship between LTPA and CVD risks. The first analysis used MSM with restricted cubic splines to fit the nonlinear association between LTPA and CVD risk and used ANOVA to test for nonlinearity. In the second analysis, we used the segmented regression to fit the piecewise-linear relationship between LTPA and CVD risk and to explore whether there are breakpoints of the LTPA time at which the increased rate of health benefits changes significantly.

All data analyses were performed by R (version 3.5.5;) with “ipw” packages and “segmented” packages; P < 0.05 was considered significant.

Results

A total of 58,167 patients (46.4% were men) were included in this study, the median baseline age was 62.2 (IQR, 53.7, 69.0) years. The median baseline SBP and DBP were 145.0 (IQR, 138.0, 154.0) mmHg and 88.0 (IQR, 80.0, 91.0) mmHg, respectively. The participants in the four baseline LTPA categories (0, 1–150, 151–300, >300 min/week) accounted for 41.6, 36.3, 15.8, and 6.4%, respectively (Table 1). While the corresponding values were 33.4, 38.8, 19.4, and 8.4% at all person-quarters (Supplementary Table 2). Participants had 15.2 (range, 3–37) times of follow-up records on average. The median follow-up time was 4.0 years (25th to 75th percentile, 2.0–4.8 years). A total of 16,332 CVD, 11,005 CHD, and 11,021 CBD were identified during the follow-up period, with an incidence density of 78.8 (95% CI, 77.7–80.0), 50.2 (95% CI, 49.3–51.2), and 50.3 (95% CI, 49.4–51.2) per 1,000 person-years, respectively. Among the participants in the four baseline LTPA categories (0, 1–150, 151–300, >300 min/week), the crude incidence of CVD events was 89.7 (95% CI, 87.7–91.6), 80.4 (95% CI, 78.5–82.3), 62.6 (95% CI, 60.0–65.3), and 44.0 (95% CI, 40.6–47.5) per 1,000 person-years, respectively. Estimated 5-year cumulative incidence of CVD for participants in four LTPA categories (0, 1–150, 151–300, >300 min/week) were 35.4, 32.2, 25.2, and 19.0%, respectively. The cumulative incidence functions for CVD, CHD, and CBD by baseline LTPA categories are shown in Supplementary Figure 3.

TABLE 1

Table 1. Characteristics of the study population by baseline physical activity level.

After multivariable adjustment in the Cox model, compared with participants in baseline LTPA 0 min/week category, the adjusted HRs associated with 1–150, 151–300, and >300 min/week LTPA categories were 0.85 (95% CI, 0.83–0.88), 0.67 (95% CI, 0.64–0.71), and 0.47 (95% CI, 0.44–0.51) for CVD, respectively. However, the analysis of the MSM-Cox with adjustment for time-dependent confounders showed that the corresponding adjusted HRs were 0.83 (95% CI, 0.77–0.91), 0.58 (95% CI, 0.53–0.63), and 0.39 (95% CI, 0.35–0.44), respectively. Table 2 summarizes the adjusted HR estimated by the MSM-Cox and Cox models. We also used LTPA as a continuous variable to conduct the Cox model and MSM-Cox analysis (Table 3). The Cox model showed that per 60 min/week LTPA increase was significantly associated with a 9.00% reduction of CVD risk, while the corresponding value estimated by MSM-Cox was 13%. The stratified analyses showed that there was no significant difference in the strength of association between LTPA and CVD among subgroups stratified by age, gender, diabetes, drinking, and antihypertensive medication (Figure 1, Supplementary Tables 3–8). In the subgroups of nonsmokers, no hyperlipidemia, or BMI <24 kg/m², the protective effect of LTPA on CVD was higher than that of other corresponding subgroups.

TABLE 2

Table 2. Adjusted hazards ratios for risk of CVD, CHD, CBD, and PAF associated with PAT.

TABLE 3

Table 3. Multivariable associations of outcomes with LTPA increase of 60 min/week estimated by MSM-Cox and Cox model.

FIGURE 1

Figure 1. Adjusted HRs (95% CIs) for (A) CVD, (B) CHD, and (C) CBD associated with time-updated LTPA estimated by the MSM. HRs, hazard ratios; CVD, cardiovascular disease; CHD, coronary heart diseases; CBD, cerebrovascular diseases; BMI, body mass index; LTPA, leisure time physical activity; PYs, person-years; P-value is for interactions of time-updated LTPA with these covariates.

The MSM with the five knots restricted cubic spline indicated that there was a nonlinearity relationship between LTPA (as a continuous variable) with CVD (P < 0.001 for nonlinearity test), CHD (P < 0.001 for nonlinearity test), and CBD (P < 0.001 for nonlinearity test) (Figure 2). Figure 3 shows the piecewise-linear relationships between LTPA and CVD risk. The estimation of LTPA breakpoint was 417 min/week for CVD. The parameters of the piecewise-linear relationship between CVD events and LTPA are summarized in Supplementary Table 9. Before the breakpoint, the protective effect of LTPA on CVD increased rapidly (slope = −0.0017), but after the breakpoint, the benefits from LTPA for protection of CVD increased slowly (slope = −0.0003).

FIGURE 2

Figure 2. The dose–response relationships between LTPA with (A) CVD, (B) CHD, and (C) CBD, which were estimated by the MSM-Cox with restricted cubic splines. CVD, cardiovascular disease; CHD, coronary heart diseases; CBD, cerebrovascular diseases; LTPA, leisure time physical activity.

FIGURE 3

Figure 3. The piecewise-linear relationships between LTPA with (A) CVD, (B) CHD, and (C) CBD, which were estimated by the segmented regression. CVD, cardiovascular disease; CHD, coronary heart diseases; CBD, cerebrovascular diseases; LTPA, leisure time physical activity.

Discussion

In this community-based longitudinal cohort of hypertensive patients, we evaluated the association between LTPA and CVD events using the Cox model and MSM-Cox model. Our study showed that LTPA had a stronger association with CVD events after appropriately adjusting for time-dependent confounders than that estimated by the Cox model with baseline LTPA. Our study also explored the piecewise-linear relationship between LTPA and CVD events and estimated breakpoints of LTPA. To the best of our knowledge, this is the first study to use MSM to evaluate the protective effect of LTPA on CVD in a longitudinal cohort with time-dependent confounders.

Many previous studies have examined the association between LTPA and CVD events and supported an inverse relationship between them (6, 14–18). However, few studies have explored the association and dose–response relationship between LTPA and CVD risk among hypertensive patients in a longitudinal cohort. Lee and Shiroma (15) reviewed the epidemiological studies on PA and cardiovascular health and found that compared with the least active subjects, the most active subjects had a 30–35% reduction of CHD risk and about 40% reduction of CVD risk. Scarborough et al. (14) found that participants achieving recommended PA levels (150 min/week) had a 23% lower risk for CVD mortality and a 17% lower risk for CVD incidence than inactive participants. Our study found that the protective effect of LTPA on the risk of CVD/CHD in hypertensive patients is greater than the results reported in the above studies. There are two main aspects of physical activity on health protective effects. The first is to improve fitness and physical function (19). The second aspect is that PA regulates the harmful relationship between cardiovascular disease and its risk factors (20). Many studies were consistent with the reverse relationship between LTPA and risk factors of CVD, such as hypertension (21), diabetes (22), obesity (8), and C-reactive protein (23).

In this study, the strength of the association between LTPA and CVD risk estimated by the Cox model was lower than that estimated by the MSM-Cox model. This finding may reflect shortcomings in previous research methods. What we were interested in was the effect of LTPA on the risk of CVD when LTPA is constant at a certain level throughout the study period, but the Cox model could not reply to this question well. Most previous studies have relied on LTPA measurements in a single time point or the average of some LTPA measurements over a certain period for association study and assumed that it was constant (6, 14–18). However, LTPA is not a static but a long-term and dynamic process that may vary over time, and BP or other factors may reduce or increase the LTPA at the next time point. The Cox models could not adequately handle time-dependent confounders. The MSM has been used in several previous studies (24–26) to adequately handle time-dependent confounders. Based on the comparison of analysis results between the Cox models and the MSM in this study, it is reasonable to speculate that the effect of LTPA on the prevention of CVD was probably underestimated in previous studies.

The association between very high PA levels and health benefits remains highly controversial (27, 28). Some studies suggested a reverse J-shape (29) or even a U-shaped association in which CVD events increase among the participants with very high PA levels (30, 31). Other studies (18, 32) suggested that very high LTPA had beneficial effects on health. In this study, we conducted two analyses to explore the dose–response relationship between LTPA and CVD risks. Both analyses showed that very high LTPA had a protective effect on CVD. In addition, we also found that there may be breakpoints of the LTPA time at which the increased rate of health benefits changes significantly. Before the breakpoint, the health benefits of LTPA on CVD prevention increased rapidly with the increase in time spent in LTPA, but after the breakpoint, the increased rate of the health benefits decreased significantly. Similar dose–response relationships were also found in the Zhao et al. (32) study.

Our study has several strengths. First, the data were collected from a relatively large-size community-based hypertension management longitudinal cohort that could represent the general community population with hypertension in Shandong Province. Second, we used the MSM-Cox model to assess the association between LTPA and CVD with adjustment for time-dependent confounders and explored the piecewise-linear relationship between LTPA and CVD. Our study also has some limitations. A major limitation of this study is that LTPA was measured by patients' self-reported information, which may be affected by recall bias. Compared with objective measurement of PA, the criterion validity of self-reported LTPA is relatively limited (33). Second, participants included in this study had to have complete data and at least three times study visits, which may lead to selection bias. In addition, the lack of data on PA at work and during transportation, as well as data on sedentary behaviors, is also a limitation of the study.

Conclusions

Our study emphasized that the effect of LTPA on the prevention of CVD in hypertensive patients was stronger than that estimated in conventional analyses based on baseline LTPA. Four hundred seventeen min per week is a breakpoint of LTPA; after achieving the breakpoint, the incremental health benefits of CVD prevention obtained from the increase in LTPA were much less than before. This should be taken into account when developing PA guidelines for hypertensive patients.

Data Availability Statement

The data used to support the findings of this study are available from the corresponding author upon request.

Ethics Statement

The studies involving human participants were reviewed and approved by the ethics committee of the School of Public Health, Shandong University. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

JL, ZZ, and FX designed the study, performed data analyses, and had full access to all the data in the study. JL drafted the manuscript. FX and SS revising the manuscript critically for important intellectual content. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by the National Key Research and Development Program of China (2020YFC2003500), Shandong Province Major Science and Technology Innovation Project (2018CXGC1210), and the Major Science and Technology Projects of Shandong province (2018YFJH0506-2).

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.

Acknowledgments

The authors thank the patients included in this study and doctors of the community health service center for their efforts in data collection. We are also very grateful to the engineers of the Healthcare Big Data Research Institute of Shandong University for their help in data collection and processing.

Supplementary Material

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

References

1. Arnett DK, Blumenthal RS, Albert MA, Buroker AB, Goldberger ZD, Hahn EJ, et al. 2019 ACC/AHA guideline on the primary prevention of cardiovascular disease: executive summary. J Am Coll Cardiol. (2019) 74:1376–414. doi: 10.1016/j.jacc.2019.03.009

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Piercy KL, Troiano RP, Ballard RM, Carlson SA, Fulton JE, Galuska DA, et al. The physical activity guidelines for Americans. JAMA. (2018) 320:2020–8. doi: 10.1001/jama.2018.14854

CrossRef Full Text | Google Scholar

3. Shinoda Tagawa T, Oguma Y. Physical activity decreases cardiovascular disease risk in women: review and meta-analysis. Am J Prev Med. (2004) 26:407–18. doi: 10.1016/j.amepre.2004.02.007

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Francesco S, Andrea C, Francesca C, Abb. Physical activity during leisure time and primary prevention of coronary heart disease: an updated meta-analysis of cohort studies. Eur J Prev Cardiol. (2008) 15:247–57. doi: 10.1097/HJR.0b013e3282f232ac

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Kurantsin Mills J, Kwagyan J, Jayam Trouth A, Weir R, Diep L. Association of physical activity level and stroke outcomes in men and women: a meta-analysis. J Women Health. (2010) 19:1815–22. doi: 10.1089/jwh.2009.1708

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Sattelmair J, Pertman J, Ding EL, Kohl HW, Lee IM. Dose response between physical activity and risk of coronary heart disease a meta-analysis. Circulation. (2011) 124:789–95. doi: 10.1161/CIRCULATIONAHA.110.010710

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Lee I, Shiroma EJ, Lobelo F, Puska P, Blair SN, Katzmarzyk PT. Effect of physical inactivity on major non-communicable diseases worldwide: an analysis of burden of disease and life expectancy. (2012) 380:219–29. doi: 10.1016/S0140-6736(12)61031-9

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Trends over 5 Decades in U.S. Occupation-related physical activity and their associations with obesity. PLoS ONE. (2011) 6:e19657. doi: 10.1371/journal.pone.0019657

CrossRef Full Text | Google Scholar

9. Ehsani AA. Exercise in patients with hypertension. Am J Geriatr Cardiol. (2001) 10:253–60. doi: 10.1111/j.1076-7460.2001.00800.x

CrossRef Full Text | Google Scholar

10. Sharman JE, La Gerche A, Coombes JS. Exercise and cardiovascular risk in patients with hypertension. Am J Hypertens. (2015) 28:147–58. doi: 10.1093/ajh/hpu191

PubMed Abstract | CrossRef Full Text | Google Scholar

11. He LI, Wei Wr, Can Z. Effects of 12-week brisk walking training on exercise blood pressure in elderly patients with essential hypertension: a pilot study. Clin Exp Hypertens. (2018) 40:673–9. doi: 10.1080/10641963.2018.1425416

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Robins JM, Hernán Má, Brumback B. Marginal structural models and causal inference in epidemiology. Epidemiology. (2000) 11:550–60. doi: 10.1097/00001648-200009000-00011

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Van der Wal WM, Geskus RB. ipw: an r package for inverse probability. J Stat Softw. (2011) 43:1–23. doi: 10.18637/jss.v043.i13

CrossRef Full Text | Google Scholar

14. Wahid A, Manek N, Nichols M, Kelly P, Foster C, Webster P, et al. Quantifying the association between physical activity and cardiovascular disease and diabetes: a systematic review and meta-analysis. J Am Heart Assoc. (2016) 5:e002495. doi: 10.1161/JAHA.115.002495

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Lee IM, Shiroma EJ. Physical activity and cardiovascular health: lessons learned from epidemiological studies across age, gender, and race/ethnicity. Circulation. (2010) 122:743–52. doi: 10.1161/CIRCULATIONAHA.109.914721

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Siscovick DS, LaCroix AZ, Manson JE, Oberman A, Mouton CP, Perri MG, et al. Walking compared with vigorous exercise for the prevention of cardiovascular events in women. N Engl J Med. (2002) 347:716–25. doi: 10.1056/NEJMoa021067

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Orsini N, Franco OH, Roberts I, Woodcock J. Non-vigorous physical activity and all-cause mortality: systematic review and meta-analysis of cohort studies. Int J Epidemiol. (2011) 40:121–38. doi: 10.1093/ije/dyq104

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Arem H, Moore SC, Patel A, Hartge P, de Gonzalez AB, Visvanathan K, et al. Leisure time physical activity and mortality a detailed pooled analysis of the dose-response relationship. JAMA Intern Med. (2015) 175:959–67. doi: 10.1001/jamainternmed.2015.0533

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Kujala, Urho M. Is physical activity a cause of longevity? It is not as straightforward as some would believe. A critical analysis. Br J Sports Med. (2018) 52:914–8. doi: 10.1136/bjsports-2017-098639

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Caceres VM, Stocks N, Adams R, Haag DG, Peres KG, Peres MA, et al. Physical activity moderates the deleterious relationship between cardiovascular disease, or its risk factors, and quality of life: findings from two population-based cohort studies in Southern Brazil and South Australia. PLoS ONE. (2018) 13:e0198769. doi: 10.1371/journal.pone.0198769

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Hu G, Barengo NC, Tuomilehto J, Lakka TA, Nissinen A, Jousilahti P. Relationship of physical activity and body mass index to the risk of hypertension: a prospective study in finland. Hypertension. (2004) 43:25–30. doi: 10.1161/01.HYP.0000107400.72456.19

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Knowler WC. Reduction in the incidence of type 2 diabetes with lifestyle intervention or metformin. N Engl J Med. (2002) 346:393–403. doi: 10.1056/NEJMoa012512

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Kasapis C, Thompson PD, FACC†. The effects of physical activity on serum c-reactive protein and inflammatory markers. J Am Coll Cardiol. (2005) 45:1563–9. doi: 10.1016/j.jacc.2004.12.077

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Hernán Má, Robins BBaJM. Marginal structural models to estimate the causal effect of zidovudine on the survival of HIV-positive men. Epidemiology. (2000) 11:561–70. doi: 10.1097/00001648-200009000-00012

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Cook NR. Use of a marginal structural model to determine the effect of aspirin on cardiovascular mortality in the physicians' health study. Am J Epidemiol. (2002) 155:1045–53. doi: 10.1093/aje/155.11.1045

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Anderson AH, Wei Y, Townsend RR, Qiang P, Chertow GM, Kusek JW, et al. Time-updated systolic blood pressure and the progression of chronic kidney disease: a cohort study. (2015) 162:258–65. doi: 10.7326/M14-0488

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Gerche AL, Heidbuchel H. Can intensive exercise harm the heart?: you can get too much of a good thing. Circulation. (2014) 130:992–1002. doi: 10.1161/CIRCULATIONAHA.114.008141

CrossRef Full Text | Google Scholar

28. Levine BD. Can intensive exercise harm the heart? The benefits of competitive endurance training for cardiovascular structure and function. Circulation. (2014) 130:987–91. doi: 10.1161/CIRCULATIONAHA.114.008142

CrossRef Full Text | Google Scholar

29. Mons U, Hahmann H, Brenner H. A reverse J-shaped association of leisure time physical activity with prognosis in patients with stable coronary heart disease: evidence from a large cohort with repeated measurements. Heart. (2014) 100:1043–9. doi: 10.1136/heartjnl-2013-305242

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Armstrong MEG, Green J, Reeves GK, Beral V, Cairns BJ. Frequent physical activity may not reduce vascular disease risk as much as moderate activity: large prospective study of women in the United Kingdom. Circulation. (2015) 131:721–9. doi: 10.1161/CIRCULATIONAHA.114.010296

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Schnohr P, O'Keefe JH, Marott JL, Lange P, Jensen GB. Dose of jogging and long-term mortality: the Copenhagen City heart study. J Am Coll Cardiol. (2015) 65:411–9. doi: 10.1016/j.jacc.2014.11.023

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Zhao M, Veeranki SP, Li S, Steffen LM, Xi B. Beneficial associations of low and large doses of leisure time physical activity with all-cause, cardiovascular disease and cancer mortality: a national cohort study of 88,140 US adults. Br J Sports Med. (2019) 53:1405–11. doi: 10.1136/bjsports-2018-099254

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Ek A, Ekblom Ö, Hambraeus K, Cider Å, Kallings LV, Börjesson M. Physical inactivity and smoking after myocardial infarction as predictors for readmission and survival: results from the SWEDEHEART-registry. Clin Res Cardiol. (2019) 108:324–32. doi: 10.1007/s00392-018-1360-x

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Robins JM. Marginal Structural Models Versus Structural Nested Models as Tools for Causal Inference. New York, NY. Springer, New York, NY (2000).

Google Scholar