Relationship between domains of physical activity and cardiac autonomic modulation in adults: a cross-sectional study

This study aimed to analyze the relationship of physical activity in different domains with cardiac autonomic modulation in adults. A sample of 252 adults was randomly selected, with mean age of 42.1 (± 16.5) years, being 58% of women. Cardiac autonomic modulation was assessed through indexes of heart rate variability in time (SDNN, RMSSD) and frequency (LF, HF) domains for linear method, and by Poincaré plot for non-linear method (SD1, SD2 components). Domains of PA (occupation, sport, leisure time/commuting, and total) were assessed by Baecke’s questionnaire. Variables of age, gender, socioeconomic status (questionnaire) and body mass index (objectively measures) were covariates. Occupational PA was positively related to LF (β = 2.39, 95% CI 0.24; 4.54), sports PA was positively related to SDNN (β = 3.26, 95% CI 0.18; 7.05), RMSSD (β = 4.07, 95% CI 0.31; 7.85), and SD1 (β = 2.85, 95% CI 0.11; 5.81), and leisure time/commuting PA was positively related to SDNN (β = 3.36, 95% CI 0.28; 6.70) and RMSSD (β = 3.53, 95% CI 0.46; 7.52) indexes. Total PA was related to RMSSD (β = 1.70, 95% CI 0.04; 3.72). Sports, leisure time/commuting, and total PA were related to higher parasympathetic modulation, while occupational PA was related to higher sympathetic modulation to the heart in adults.


Methods
This cross-sectional study follows the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) and was approved by the Ethical Research Committee of São Paulo State University under protocol CAAE: 72191717.9.0000.5402. All methods were performed in accordance with the relevant guidelines and regulations, as well as all the participants signed the Informed Consent Term agreeing to take part of the research, being informed about all procedures, without any costs, and free to desert of the research at any moment. The protocol of this research was previously registered at ClinicalTrials.gov, number NCT03986879. The data collection occurred between December 2018 and June 2019, and the study was carried out in the city of Santo Anastacio, located in the southeastern region of Brazil, which has a population of ~ 16,000 adults and a Human Development Index of 0.753 (a measurement of life expectancy, education and per capita income, where 1 is the best and 0 the worse). Adults aged 18 years or over were randomly selected to participate in this study. To calculate the sample size, it was considered a correlation value of r = 0.23 for the relationship between PA and cardiac autonomic modulation based on previous studies 6 , a sampling power of 80%, and an alpha error of 5%, which totalized a minimum sample size of 147 subjects. However, predicting possible errors in the answer to the questionnaires or errors in the series of RR intervals from the cardiac autonomic modulation analysis, an additional 30% was added, resulting in 192 subjects. Besides that, the present study considered socioeconomic status, age, gender, and body mass index as adjustment variables, which resulted in the inclusion of 15 participants for each variable, resulting in a final sample size of 252 subjects.
Regarding the sample collection process, all the census sectors of the urban area of the city were considered (23 sectors), as well the number of inhabitants by sector, so the proportionality of each sector was considered (the proportion of population living in each sector was considered as the proportion of the sample to be assessed). Neighborhoods, streets and finally households of each census sector were randomized using the "Random" function of the SPSS statistical package. The detailed sample process was previously described in literature 13 .
For data collection, previously trained researchers visited the selected households and applied the questionnaire in those who agreed to participate in the research. After this procedure, the assessment of cardiac autonomic modulation performed by heart rate variability (HRV), was scheduled within the next seven days in a previously stipulated place. physical activity. The habitual practice of PA was assessed by the Baecke's questionnaire 14 , previously validated to Brazilian adults 15 , and with reliability tested against gold standard methods as doubly labeled water 16 . This instrument assesses PA through 16 questions in the last twelve months considering: physical effort at work (8 questions), the practice of sports activities or systematized exercises (4 questions), and physical activities at leisure time and commuting (4 questions). In domain of occupational activities, questions related to physical efforts at work are investigated, such as: time spent sitting, standing and walking in the work/occupational environment, carrying weights, frequency of perspiring a lot, and how tired feels after a workday. In sports domain, the practice of leisure-time exercises is considered, such as training at the gym or playing sports. In these domains, the intensity of these activities (light, moderate, vigorous), the number of hours per week in which these activities are practiced, and how long this physical activity has been practiced (< 1 month, 1-3 months, 4-6 months, 7-9 months, > 9 months). In the domain of leisure time and commuting, different activities in free time are considered, as well as the amount of time at commuting by walking and cycling to perform daily activities (i.e. for going to the market, shopping or going to work).
At the end, this instrument offers a dimensionless score for each domain, which ranges from 1 to 5, and the sum of the three scores corresponds to the total physical activity score. cardiac autonomic modulation. The assessment of cardiac autonomic modulation was performed through HRV analysis. For this assessment, the subjects were instructed not to ingest alcoholic beverages and caffeinated drinks as coffee and tea, as well as not to perform moderate-to-vigorous intensity physical activity for a minimum period of 12 h so they could not cause changes in the autonomic modulation 17 . HRV was assessed by the Polar monitor, model V800 (Polar Electro Oy, Kempele, Finland), recorded beat-to-beat in a room with controlled temperature between 21 to 24° C and relative humidity between 50 to 60% 17 .
For HRV analysis, beat-to-beat heart rate was recorded for a period of thirty minutes, with participants lying in the supine position. In the series of RR intervals obtained, a digital filtering was performed and further complemented by manual filtering, aiming to eliminate premature ectopic beats. In this filtering, only series with more than 95% of sinus beats were considered and 1000 RR intervals of the most stable period of the tracing was used for analysis 18 . If any participant had less than 95% of sinus beats, the HRV recording and analysis were performed again.
The HRV analysis was performed using linear methods (time and frequency domains) and non-linear methods. For linear methods, frequency domain was analyzed by the Low Frequency (LF-0.04 to 0.15 Hz) and High Frequency (HF-0.15 to 0.4 Hz) spectral components, expressed in normalized units, and time domain was analyzed by using RMSSD and SDNN indexes. The RMSSD corresponds to the square root of the square mean of the differences between adjacent normal RR intervals in a time interval, while SDNN is the standard deviation of mean of all normal RR intervals, both expressed in milliseconds. The spectral analysis was calculated using the Fast Fourier transform algorithm 19 . For HRV analysis using non-linear methods, the Poincaré plot was quantitatively analyzed using the following indexes: SD1 (standard deviation of instantaneous beat-to-beat variability) and SD2 (long-term standard deviation of continuous R-R intervals).
For HRV analysis in linear and non-linear methods, Kubios HRV Analysis software version 2.0 (Kupio University, Finland-URL: www.kubio s.com) was used. . For sample characterization, participants with BMI between 25 and 29.9 kg/m 2 were classified as "overweight" and those participants with BMI equal or above 30 kg/m 2 was classified as "obesity", according to global classification for adult population 21 .

Statistical analysis.
Sample characterization variables were presented as mean and standard deviation, with gender differences compared by independent t test. The correlation between the different domains of physical activity and the indices of cardiac autonomic modulation were verified by Pearson's correlation. The magnitude of the relationships between the PA domains and cardiac autonomic modulation were observed by linear regression models, adjusted by age, gender, socioeconomic status, and body mass index. An dummy variable gender × domain-specific physical activity was used as adjustment for interaction between independent variables. The confidence interval adopted in the present study was 95% and the statistical significance was fixed in p < 0.05. The statistical package used was IBM SPSS Statistics for Windows, version 24.0 (IBM Corp., Armonk, N.Y., USA-URL: www.ibm.com/produ cts/spss-stati stics ).
patient and public involvement. The patients and the public were not involved in the design, or conduct, or reporting, or dissemination plans of this research.

Results
The sample of the present study (Table 1) comprised 252 participants, being 147 women and 105 men. The characteristics of sample are presented in Table 1.
The socioeconomic status was significantly correlated with physical activity, as better socioeconomic status were correlated to lower Occupational score (r = 0.129, p = 0.042) and higher Sport score (r = − 0.293, p = 0.001), whereas no significant correlation was observed between socioeconomic status and HRV indexes.  www.nature.com/scientificreports/ Table 2 shows the correlations between cardiac autonomic modulation and the different PA domains. Sport, leisure time/commuting, and total PA were positively related to SDNN, RMSSD, SD1 and SD2 indexes. Occupational PA was not correlated with cardiac autonomic modulation.
Data of the multivariate analysis are presented in Tables 3 and 4. Sport PA was positively related with SDNN, RMSSD and SD1. Each increment in Sport PA score was associated with an increase of about 4 units of milliseconds in the RMSSD and 3 units of milliseconds in the SDNN, which corresponds to a positive relationship in parasympathetic modulation and global variability, respectively. Occupational PA was positively related to LF nu. Leisure time/commuting PA was positively related to SDNN and RMSSD, while Total PA was only positively related to RMSSD index in adjusted analysis.
The interaction term gender × physical activity was significant related in the models of Occupational PA vs. LF, Sport PA vs. SDNN and SD2, Leisure time/Commuting PA vs. SDNN and SD2, and of Total PA vs. SDNN and SD1 presented in Tables 3 and 4. According to gender, Occupational PA vs. LF was significant only in females (β = 2.87, p = 0.029); Sports PA was related with SDNN only in males (β = 5.27, p = 0.032) and its relationship with SD2 lost statistical significance for both males and females; Leisure time/Commuting PA was significant related with SDNN (β = 6.29, p = 0.023) and with SD2 (β = 7.24, p = 0.047) only in males; and Total PA was related with SDNN (β = 3.43, p = 0.008) and SD1 (β = 3.26, p = 0.006) only in males.

Discussion
This study showed different associations between domains of physical activity with cardiac autonomic modulation in adult population. Sport and Leisure time/commuting were related to higher parasympathetic modulation to the heart, being evident in several HRV indexes. Table 2. Correlation between cardiac autonomic modulation and different domains of physical activity. HRV Heart rate variability, SDNN Standard deviation of the mean of all normal RR intervals, RMSSD Square root of the square mean of the differences between adjacent normal RR intervals in a time interval, LF Low frequency, HF High frequency, SD1 Standard deviation of instantaneous beat-to-beat variability, SD2 Long-term standard deviation of continuous R-R intervals, PA Physical activity, ms milliseconds, nu normalized units.

HRV indexes
Occupational PA Sport PA www.nature.com/scientificreports/ It was observed that occupational PA score was positively related with LF index of HRV in the present study. Occupational activities could be, in general, less intense than physical activities in other domains, as well as may be susceptible to exposure of environmental conditions related to work, as physical and mental stress, worries, and work overload. Besides that, high levels of job stress were associated with increased heart rate and decreased HRV indexes 22 . In this sense, this relationship could be a subject to further additional investigation.
Another important aspect to highlight is that LF index corresponds to a contradictory interpretation in cardiac autonomic modulation 23 . Traditionally, LF index was considered as a reflection of the sympathetic nervous system activity 24 , but a more recent study considered LF as an index of baroreflex function in the modulation of cardiac autonomic outflows, not being a measure of cardiac sympathetic tone 25 . Furthermore, a study with animal model reported that LF was neither a robust index of cardiac sympathetic nerve activity nor of baroreflex sensitivity, but highlight a link between them 26 .
PA sport domain was positively related to RMSSD, SDNN, and SD1 indexes in adjusted analysis in this study, which corresponds to higher modulation of parasympathetic nervous system and global variability 19 . While RMSSD represents parasympathetic modulation by estimating HRV changes mediated by vagal tone 27 , variations in SDNN index does not have distinction whether changes are related to sympathetic or vagal tone 28 . In turn, despite being considered as an identical metric to the RMSSD 29 , SD1 index represents a short-term HRV measurement which is correlated to baroreflex sensitivity 30 and predicts a set of different factors which may reflect the complexity of mechanisms in HRV regulation, as diastolic blood pressure and HRV in time and frequency domains 31 . This increase in parasympathetic activity could be related to the vigorous-intensity of the great part of sport activities (higher than 6 METs) 32 , once it has been observed an increase in RMSSD index according to the intensity of physical activity 4 . Schmidt et al. 33 observed that sports activity showed higher relationship with physical fitness and physical health status than habitual activity in 18 years of follow up. Based upon these evidences, it is hypothesized that high-intensity training may result in a positive change in the vagal activity 34,35 .
Physical activities in leisure time and commuting were positively related to RMSSD and SDNN indexes in the present study, regardless of gender, age, socioeconomic status, and body mass index. Likewise, Bueno et al. 36 observed that adults that cycling for transportation presented higher parameters of parasympathetic modulation and lower parameters of sympathetic modulation than insufficiently active adults. Soares-Miranda et al. 37 reported that higher leisure time physical activity, walking distance, and walking pace were associated to higher HRV indexes in older adults. The findings of the present study also raise the hypothesis that the time spent in leisure time and commuting activities take place from the time spent in sedentary behavior. Delfino et al. 38 observed that those adults with high breaks in sedentary behavior at leisure and occupation showed high levels of habitual physical activity. In this sense, higher levels of sedentary behavior may result in less light-intensity physical activity, which may negatively impact cardiac autonomic modulation, once sedentary lifestyle was associated to lower RMSSD index and higher sympathovagal balance in adults 39 .
The present study also showed a positive relationship between total PA and RMSSD after adjusting for variables of interest. Corroborating with this finding, Tornberg et al. 4 observed a positive relationship between self-reported PA and RMSSD index regardless of adiposity levels, and Soares-Miranda et al. 7 reported that the most active group showed higher levels of vagal HRV indexes than the less active group in adult population. It is important to highlight that total PA in this study was a sum of the three different domains assessed by Baecke's questionnaire, which were composed by different constructs in its scoring, as weekly frequency, perceived intensity, duration and how long these activities were performed, being composed by domains which was more related www.nature.com/scientificreports/ to HRV indexes than others, which could mitigate this relationship when analyzed in total score. Therefore, the dimensionless of Baecke score did not allow comparisons between physically actives and insufficiently actives according to global recommendations, or even inferences about how much intensity of amount of habitual PA in the assessed domains were related to cardiac autonomic modulation. Even the analysis adjusted by potential confounding factors, the interaction between gender and physical activity may limit the present study findings. Previous study reported a significant gender difference in the activity of autonomic nervous system 40 and was observed that men presented a predominance of sympathetic vascular regulation, while women showed a dominant influence of parasympathetic heart rate regulation 41 . In this sense, further investigation with gender-specific analysis is suggested.
This study has some limitations. It is important to highlight the observational design of this study, which does not allow cause and effect inferences. Further, the assessment of physical activity by questionnaire was susceptible to bias of memory and classification of intensity. Besides that, the use of medicines for cardiovascular control, the hour of HRV assessment, and menstrual cycle of women were not controlled. However, the representative sample randomly selected and analysis adjusted by potential confounding factors (gender, age, socioeconomic status, and body mass index) were positive strengths of the study. Besides that, the use of a questionnaire for physical activity allows to analyze its practice in different domains, which was the main purpose of this study, and would not be provided by accelerometer.
In conclusion, the association between physical activity and cardiac autonomic modulation in adults is dependent on the domain analyzed. Sport, leisure time/commuting, and total PA were positively related to parasympathetic modulation indexes, whereas occupational PA was related with sympathetic modulation indexes. As practical application, the stimulation of practice of PA in different domains, mainly in regard sports practice and leisure time/commuting activities may contributes to improvement in cardiac autonomic modulation in adult population, regardless gender, age, socioeconomic status, and body mass index.