Introduction

Atrial fibrillation (AF) is the most common arrhythmia and is an independent risk factor for cardiovascular events including heart failure, and ischemic stroke1. Risk factors for the development of AF include hypertension, diabetes, obesity, smoking, ischemic heart disease, and heart failure.

Association between physical activity (PA) and AF has also been studied, and moderate PA is reported to be associated with a reduced risk of AF, while higher PA is not2,3. On the other hand, a meta-analysis showed that the trend toward lower risk of AF in subjects practicing PA is non-significant and that the extreme heterogeneity of published studies does not allow conclusive results4.

There are other issues related to PA. While previous studies divided PA into the light, moderate, and vigorous activities, none of them quantified the effect of total daily life PA. In addition, it is still not completely clear whether the PA and resting heart rate (RHR) have an independent effect on the development of AF. Endurance athletes are more likely to develop AF than non-athletes, with sinus bradycardia and increased vagal tone as possible pathophysiological mechanisms5. Therefore, PA and RHR may presumably also influence each other in the general population, and a relation between RHR and incidence of AF has indeed been suggested6,7. However, there are limited data about the interplay between PA and RHR in the development of AF8. Further, there have been no large-scale studies analyzing the effect of PA and RHR on the incidence of AF in an Asian population.

Our study aimed to assess (1) whether the total PA and RHR have an independent effect on incident AF and (2) whether the amount of total daily activity has any influence on the risk of AF in the general population based on a prospective population study carried out in South Korea.

Methods

Study population

The Ansung-Ansan cohort is a prospective study population that consists of South Koreans aged 40–69 years residing in two cities, Ansan and Ansung. This study was started in 2001 and embedded within the Korean Genome and Epidemiology Study (KoGES), a population-based cohort study, to assess the incidence and risk factors for various chronic disorders including cardiovascular diseases.

The baseline survey, conducted from 2001 to 2002, included 10,030 adults. Following the baseline visit, these adults were examined every two years with the longest follow-up duration being 12 years. At each of the baseline and biennial follow-up visits, data on lifestyle characteristics, cardiovascular profiles, medical history, and disease incidence were collected. All medical tests and the interview-based surveys were conducted by health professionals who were trained to follow a standardized protocol, and the same instruments were used to collect the data. More detailed information regarding the enrollment of cohort members and study procedures is available in a previous report9.

Analytical sample

Among the 10,030 subjects, we analyzed 8,811 (4,235 males and 4,576 females), excluding participants whose heart rate measurement was not available (n = 15), those with AF at baseline (n = 40), those with no follow-up information (n = 965), and those with missing information on covariates (n = 199). All of the study protocol were approved by the Ethics Committee of KoGES at the Korean National Institute of Health and adhered to the principles of the Declaration of Helsinki. Informed consent was obtained from all subjects. This study was approved by the Institutional Review Board of Hanyang University Medical Center (IRB file No. HYUH 2017-12-033).

Physical activity measurement

PA was assessed using a questionnaire to quantify activities in occupational or leisure time domains and was categorized into 5 levels: steady state, sitting state, and light, moderate, and high activity. A corresponding metabolic equivalent (MET) was assigned for each level as follows: (1) steady state (MET 0) - lying down; (2) sitting state (MET 1.5) - typing, driving, office work, playing a musical instrument, sewing, studying, ironing, writing, and cooking; (3) light activity (MET 3) - walking, cleaning, laundry, child care, bathing, and exercise as entertainment (bicycle, table tennis); (4) moderate activity (MET 5) - fast walking, woodworking, lawn mowing, snow removal, and regular exercise (badminton, swimming, tennis); and (5) high activity (MET 7) - athletics, climbing, running, woodcutting, agriculture, forestry, and mining. The questionnaire included questions about whether participants performed each of the above 5 levels of PA ≥ 1 h/day or < 1 h/day. The total PA in MET-h/day was based on the PA intensity and the time commitment for each PA.

Measurement of resting heart rate and other covariates

RHR was recorded by measuring the pulse for 60 s. Covariates including baseline age, sex, urban or rural residence, education level, past medical history, smoking status, and alcohol consumption were evaluated based on self-reports. Blood pressure (BP) was measured by well-trained nurses after a 5-min rest. Hypertension was defined as systolic BP > 140 mmHg, diastolic BP > 90 mmHg, or the use of antihypertensive medication. We defined a participant as diabetic if any of the following criteria were met: current use of an oral hypoglycemic agent or insulin or a serum fasting glucose ≥126 mg/dl. A participant was categorized as having hyperlipidemia if any of the following criteria were met: current use of lipid-lowering medication, serum total cholesterol ≥ 240 mg/dl, or serum triglyceride ≥200 mg/dl.

Follow-up and detection of AF

Biennial follow-up examinations and electrocardiograms (ECGs) were performed until 2014. A resting supine 12-lead ECG (MAC 5000®; GE Marquette Inc., Milwaukee, WI, USA) was obtained during quiet respiration. The ECG was recorded at 25 mm/s with 0.1 mV/mm standardization and was interpreted by a cardiologist to confirm the diagnosis and classified according to the Minnesota code. AF was defined as a composite of AF or atrial flutter in the standard 12-lead ECG.

AF at baseline was identified by ECG performed during the baseline visit and/or the self-reported history of physician-determined diagnosis made prior to the baseline visit. Similarly, AF newly developed after the baseline visit was also identified by biennially performed ECGs and/or the self-reported history of physician-determined diagnosis that was made between each biennial visit. If AF was identified by both ECGs and the history in the same subject, the earlier identification date was considered the time of AF development.

Statistical analyses

Continuous variables were presented as a median and interquartile range; Kruskal-Wallis test and Dunn’s multiple comparison test (post-hoc analysis), comparing all groups to controls were used to compare continuous variables among groups. Categorical variables were presented as numbers and percentages. The chi-squared test was used for categorical variables. A Cox proportional hazards regression model was used to assess the association between PA, RHR, and incidence of AF. The collected total PA (MET-h/day) was divided into quartiles, with the lowest and the highest being Q1 and Q4, respectively. Q1 was used as the reference. RHR was categorized into 5 groups (<50, 50–59, 60–70, 70–85, and >85 bpm). The RHR range of 70–85 bpm was used as a reference. The covariates were age, sex, residence, education, body mass index (BMI), comorbidity, smoking status, and alcohol consumption in models 1 and 2. In model 3, we added RHR or PA to evaluate the independent association with the incidence of AF. All statistical analyses were performed using R-3.4.0 software. A p < 0.05 was considered significant.

Results

Baseline demographic and health-related characteristics

The baseline characteristics according to the PA category are listed in Tables 1 and 2. A total of 8,811 participants (48.1% male) were followed up for a median period of 139 (97, 141) months. The median age was 50 (44, 60) years. There were several significant differences between participants in the lowest quartile (Q1) and the highest quartile (Q4): The median age of participants was higher, there were more male participants, the % of urban residents (=urban residents/(urban + rural residents) was lower, and the overall education level of participants was lower in Q4 than in Q1 (Tables 1 and 2). There was a significant difference in the distribution of occupations. The large majority (62.7%) of participants engaged in agriculture belonged in Q4. Conversely, the smallest proportion of participants with all other occupations belonged in Q4 (Table 1).

Table 1 Demographic characteristics according to quartiles of physical activity.
Table 2 Baseline health-related characteristics according to quartiles of physical activity.

BMI was significantly lower in Q4 than in Q1. However, the differences in age and BMI between Q3 and Q1 were not significant. The frequencies of hyperlipidemia and thyroid disease were lower among participants in Q4 than those in the other three quartiles. There were more current or former smokers in Q4 than in other quartiles but no significant differences in alcohol consumption among four quartiles (Table 2).

Physical activity and resting heart rate

Only the total PA measured at baseline was analyzed for the following reasons: 1) The total number of follow-up visits varied among participants; therefore, using the baseline data that were available in every participant would likely provide the most generalizable results. Further, including all data from varying numbers of follow-up visits could complicate the analysis, presentation, and interpretation of data. 2) After the baseline visit, there was some modification in the questionnaire for assessing PA for reasons unrelated to the present study; therefore, METs/day could not be calculated throughout all visits in a consistent manner.

For each component of PA according to total PA quartiles (Table 2), light activity ≥1 h/day was performed by more than half of participants in all quartiles except Q1. The percentage of participants performing moderate activity ≥1 h/day was highest (45.4%) in Q3. While the percentage of participants in Q4 performing moderate PA ≥ 1 h/day (33.5%) was not as high, 90.4% of participants in Q4 performed high PA ≥ 1 h/day (vs. 24.3% in Q3). More detailed components of PA in each of the four quartiles of total PA are shown in Supplementary Fig. 1.

Only baseline RHR data was analyzed for the same reason as “Reason 1” mentioned in the first paragraph of this section and also because initial rankings of 5 RHR groups determined based on median values of baseline RHR remained identical throughout all follow-up visits (Supplementary Fig. 2).

Comparing the baseline RHR among quartiles (Table 2), Q3 included the highest percentage of subjects with RHR < 60 bpm (27%), and the RHR was significantly lower in Q3 than in Q1 (62.0 [58.7, 67] bpm, p < 0.001). On the other hand, Q4 had the lowest percentage with RHR < 60 bpm (18%).

Physical activity, resting heart rate and incident AF

There was a gradual decrease in the incidence rate and risk of AF from Q1 to Q3 (Table 3); Q3 had the lowest incidence rate of AF (1.42/1,000 person-years) and a significantly lower risk of AF, after adjusting for covariates, than Q1 and Q4. This trend remained consistent after adjusting for RHR (hazard ratio [HR] 0.58, 95% confidence interval [CI] 0.37–0.90, p = 0.016). On the other hand, Q4 had a high incidence rate of AF (2.12/1,000 person-years); the risk of AF was not even lower than in Q1 (HR 0.88, 95% CI 0.57–1.36). These associations between PA and incident AF were similar in both male and female subjects, but were more significant in females (Fig. 1).

Table 3 Hazard ratio for incident atrial fibrillation according to quartiles of total physical activity.
Figure 1
figure 1

Risk of atrial fibrillation according to the level of physical activity.

However, when participants in Q3 and Q4 were divided into subgroups, i.e., high PA ≥ 1 hr/day versus high PA < 1 hr/day (Fig. 2), subgroup of Q4 with high PA ≥ 1 hr/day had a significantly higher risk of AF compared with the subgroup of Q3 with high PA < 1 hr/day. Conversely, the risk of AF in the subgroup of Q4 with high PA < 1 hr was not higher than in the subgroup of Q3 with high PA < 1 hr.

Figure 2
figure 2

Incidence of atrial fibrillation in relation to presence or absence of high activity (MET 7) ≥ 1 hr/day in Q3 and Q4. p < 0.05 vs. Q3 with high PA < 1 hr/day. Kruskal-Wallis test (Dunn’s post hoc).

Table 4 shows that subjects with RHR < 60 had a higher risk of AF than those with RHR 70–85 bpm. Subjects with RHR 50–59 bpm had a significant risk of AF (HR 1.97, 95% CI 1.17–3.32, p = 0.011). These relationships remained consistent after adding PA to the covariates.

Table 4 Hazard ratio according to resting heart rate and incident atrial fibrillation.

Discussion

Analysis of a prospective cohort based on the general population showed that moderate total PA was associated with a reduced risk of incident AF compared with that for lowest total PA and that the protective effect of PA was attenuated in subjects with highest total PA. Also, low RHR (<60 bpm) was a risk factor for incident AF. These findings remained consistent after adjusting for PA.

Endurance athletes have a lower RHR and are also predisposed to AF10. Suggested mechanisms include atrial remodeling, atrial ectopy, and increased vagal tone11. The risk of AF increases substantially in the aging athlete with accumulation of lifetime training hours and participation in the competition. However, our study showed that the association between PA and RHR in the general population was not straightforward. There appeared to be some inverse association between PA and RHR up to moderate total PA. However, the highest total PA was not associated with lower RHR. These results suggest that regular PA (even the most vigorous PA) in the general population may be insufficient to affect autonomic tone. Therefore, RHR could not mediate an association between PA and incident AF, similar to the finding in a previous study3. Indeed, only 1.6% of participants in our study had RHR < 50 bpm.

Previous meta-analyses found no significant association between PA and AF, or only a non-significant trend4,12. However, while these studies only focused on each level of PA to investigate the incidence of AF, various levels and types of PA occur during daily life in the general population, which may have resulted in inconsistent data for the relationship between each level of PA and AF. To that end, we analyzed the relationship between PA and AF using individual total PA levels, which have not been done by others. We quantified daily life PA by using an index based on METs and observed the association between PA and AF in the general population. Interestingly, many in the Q4 group performed some light- and moderate-level PA. However, almost 90% of Q4 subjects performed high-level PA, which attenuated the favorable effect of moderate PA. These findings are consistent with previous studies that showed a higher risk of AF in high- or vigorous-activity8 and high level-PA attenuated the benefit of moderate PA3. Several pathways have been suggested to explain the mechanism by which PA affects the risk of AF13. Moderate PA modulates cardiovascular risk factors and markers of inflammation14,15. These positive effects may reduce the risk of AF.

Low RHR was associated with high AF incidence in subjects aged ≥65 years16, and low exercise heart rate based on a moderate workload was a long-term predictor of incident AF in middle-aged men17. Although the authors used a different range of heart rate as a ref.16,17, they consistently suggested that low RHR is associated with high AF. We also found low RHR being associated with a high incidence of AF; these associations were independent of PA. While low RHR was associated with higher incident AF in trained athletes5, participants in our study were not trained athletes. Therefore, it is conceivable that mechanisms other than autonomic change, atrial dilation, and fibrosis may play a role. One possible mechanism is subclinical sinus node dysfunction, which may predispose to AF18. Increased heart rate was also found to be associated with risk of AF in the Losartan Intervention For Endpoint Reduction in Hypertension (LIFE) study19. In the present study, the incident rate of AF and hazard ratio in the ‘RHR > 85 bpm’ group were higher than in the ‘RHR 70–85 bpm’ group but not statistically significant, possibly due to a small number of subjects in this group. A study by Aladin et al20. reported that the significant association between RHR > 85 bpm and risk of AF disappeared after adjusting for exercise capacity.

The incidence of AF in our study was lower than in previous studies of Western populations21,22. For example, Schnabel et al. reported an incidence of 7.5/1,000 person-years, based on 50 years of observation in the Framingham cohort22. The difference in the incidence may be explained by different populations and the methods to identify AF. The prevalence of AF is reportedly 250–325/100,000 in Northeast Asia, which is the lowest in the world23. Also, a recent report based on the Korean population showed a low prevalence and incidence rate of AF (0.67% and 1.5–1.7/1,000 person-years, respectively) in a Korean population, similar to our findings24. We used clinic ECGs to identify AF. While this method can increase the sensitivity for AF in the general population, paroxysmal AF would not have been identified, and therefore the prevalence would probably have been underestimated. Indeed, Schnabel et al. reported that the incidence of AF based on clinic ECGs was 1.6–3.8/1,000 person-years; this was much lower than the rate based on all available sources and was closer to that in our study22. To minimize underestimation of incident AF, we also included the data using self-reported physician-determined diagnosis.

Our research has some limitations. In baseline data collection, we did not distinguish between occupational and leisure time activities when calculating total daily PA. In addition, self-reporting questionnaires may not have accurately reflected the level of PA. While we tried to identify AF using both ECGs and self-reported history of physician-determined diagnosis, we could not obtain participants’ International Classification of Disease (ICD) because ICD was considered protected personal information. Therefore, some paroxysmal AF that waxed and waned between each visit could still have been missed. Information about antihypertensive drugs such as beta blockers was not available for this cohort. However, there were no significant differences in the rate of antihypertensive use according to quartiles.

In conclusion, this prospective cohort study of a Korean general population suggested independent associations between total daily PA and incident AF as well as between RHR and incident AF. A moderate amount of total PA was associated with a lower risk of AF. Higher total PA was not associated with a lower RHR. RHR was also associated with incident AF, showing a U-shaped relationship, independent of the amount of total PA. Further prospective research is warranted to clarify whether the intervention of PA and RHR would be beneficial to prevent AF.