Next Article in Journal
Contribution of Solar Radiation and Pollution to Parkinson’s Disease
Previous Article in Journal
Predictors, Mediators and Moderators of Police Work-Related Stress: A Scoping Review
Previous Article in Special Issue
Role of Body Mass and Physical Activity in Autonomic Function Modulation on Post-COVID-19 Condition: An Observational Subanalysis of Fit-COVID Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Autonomic Function Recovery and Physical Activity Levels in Post-COVID-19 Young Adults after Immunization: An Observational Follow-Up Case-Control Study

by
Ana Paula Coelho Figueira Freire
1,2,3,*,
Shaan Amin
1,
Fabio Santos Lira
3,4,
Ana Elisa von Ah Morano
3,
Telmo Pereira
4,5,6,
Manuel-João Coelho-E-Silva
4,
Armando Caseiro
5,6,7,
Diego Giulliano Destro Christofaro
8,
Vanessa Ribeiro Dos Santos
3,
Osmar Marchioto Júnior
3,
Ricardo Aurino Pinho
9 and
Bruna Spolador de Alencar Silva
3
1
Department of Health Sciences, Central Washington University, Ellensburg, WA 98926, USA
2
Physiotherapy Department, Universidade do Oeste Paulista (UNOESTE), Presidente Prudente 19050-920, Brazil
3
Exercise and Immunometabolism Research Group, Postgraduate Program in Movement Sciences, Department of Physical Education, Universidade Estadual Paulista (UNESP), Presidente Prudente 19060-900, Brazil
4
Faculty of Sport Science and Physical Education, University of Coimbra, CIDAF, 3000-456 Coimbra, Portugal
5
Polytechnic of Coimbra, ESTESC, 3046-854 Coimbra, Portugal
6
Laboratory for Applied Health Research (LabinSaúde), 3046-854 Coimbra, Portugal
7
Molecular Physical-Chemistry R & D Unit, Faculty of Science and Technology, University of Coimbra, 3004-535 Coimbra, Portugal
8
Postgraduate Program in Movement Sciences, Department of Physical Education, Universidade Estadual Paulista (UNESP), Presidente Prudente 19060-900, Brazil
9
Graduate Program in Health Sciences, School of Medicine, Pontificia Universidade Catolica Do Parana, Curitiba 80215-901, Brazil
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2023, 20(3), 2251; https://doi.org/10.3390/ijerph20032251
Submission received: 23 November 2022 / Revised: 19 January 2023 / Accepted: 25 January 2023 / Published: 27 January 2023

Abstract

:
Coronavirus disease 2019 (COVID-19) has detrimental multi-system consequences. Symptoms may appear during the acute phase of infection, but the literature on long-term recovery of young adults after mild to moderate infection is lacking. Heart rate variability (HRV) allows for the observation of autonomic nervous system (ANS) modulation post-SARS-CoV-2 infection. Since physical activity (PA) can help improve ANS modulation, investigating factors that can influence HRV outcomes after COVID-19 is essential to advancements in care and intervention strategies. Clinicians may use this research to aid in the development of non-medication interventions. At baseline, 18 control (CT) and 20 post-COVID-19 (PCOV) participants were observed where general anamnesis was performed, followed by HRV and PA assessment. Thus, 10 CT and 7 PCOV subjects returned for follow-up (FU) evaluation 6 weeks after complete immunization (two doses) and assessments were repeated. Over the follow-up period, a decrease in sympathetic (SNS) activity (mean heart rate: p = 0.0024, CI = −24.67–−3.26; SNS index: p = 0.0068, CI = −2.50–−0.32) and increase in parasympathetic (PNS) activity (mean RR: p = 0.0097, CI = 33.72–225.51; PNS index: p = 0.0091, CI = −0.20–1.47) were observed. At follow-up, HRV was not different between groups (p > 0.05). Additionally, no differences were observed in PA between moments and groups. This study provides evidence of ANS recovery after SARS-CoV-2 insult in young adults over a follow-up period, independent of changes in PA.

1. Introduction

The coronavirus disease (COVID-19) is a rapidly spreading condition caused by SARS-CoV-2 virus that infiltrated global populations at alarming rates [1]. Despite severe complications that may arise through contracting COVID-19, a significant number of individuals did not require intensive care (mild to moderate cases) and continued to present with physical [2], neurological [3,4], and other autonomic nervous system (ANS)-related dysfunctions [2,4,5,6,7] four to six months after diagnosis [2,6].
Autonomic modulation was shown to be impaired, even in young adults short after mild and moderate COVID-19 [8,9,10]. The presence of SARS-CoV-2 virus within the carotid body can be a possible mechanism explaining the observed silent hypoxemia and, thus, providing a pathway for nervous system infiltration and subsequent ANS dysfunction [2,8,11,12]. Further, disruptions in autonomic regulation by a viral pathogen are associated with a cytokine storm immune response, which results in oxidative stress leading to cell damage [13].
Cases of COVID-19 among young adults are concerning worldwide since this population can account for 70% of those infected globally [14]. These findings might be attributed to the fact that they are of working age, with high mobility and numerous interpersonal interactions. Although the understanding of ANS impairments after COVID-19 infection in this population is still under investigation, previous research demonstrated early autonomic alterations, including increases in central sympathetic drive and vagal control of the heart [8,9]. ANS changes may be affected by age, body mass, degree of physical activity, and time since diagnosis [15]. Although young adults may present better physiological status and more efficient autonomic regulation when compared to elderly individuals, [16] the effects of COVID-19 on autonomic function over time still need further investigation.
It is well established that ANS plays a major role in modulating homeostasis through influencing systemic bodily functions. Dysfunction of the ANS could result in detrimental effects of downstream physiological processes associated with respiratory, vascular, immune, hematological, and renal processes [17,18]. Increases in SNS activity have also been identified to be an independent predictor of mortality in association with several diseases [19,20]. Considering previous studies that reported persistent symptoms and impairments as a result of post-acute sequelae of COVID-19 (PASC), it is vital to investigate ANS behavior for longer periods of time following infection by SARS-CoV-2 [8,21,22].
Additionally, investigating factors that can influence better ANS outcomes after COVID-19 is essential to advancements in care and intervention strategies. Previous research has found that physical activity (PA) levels can directly influence autonomic modulation [23,24] and can be associated with reduced risk and severity of COVID-19 symptoms [9,21,22]. Concurrently, PA is related to enhancements in immunity and cardiorespiratory fitness and can help to prevent and treat obesity, cardiovascular disease, diabetes, liver disease, cancer, and other chronic diseases [25,26], thus, indirectly reducing the threat of COVID-19 [26]. These findings highlight the magnitude of the role of PA in regulating multiple autonomic system routes to maintain homeostasis and may provide a protective effect from SARS-CoV-2 infection.
An improved understanding of the influences of the long-term effects of COVID-19 on ANS and, consequently, multiple body systems is vital for developing public-health strategies that can reduce the impact of the disease. Clinicians may consider evidence from this literature as a means for guidance to improve clinical assessment and develop non-medication intervention strategies in the prevention and recovery from infection in young adults.
Recently, our group demonstrated the short-term effects of COVID-19 on ANS in young adults [9], but information on long-term effects of COVID-19 still requires research. Young adults are an overlooked population, possibly due to the fact that they are the least at risk for severe negative outcomes. However, the high incidence in young adults [27] emphasizes the need for research into the effects of COVID-19 on this population.
Therefore, the primary aim of this study was to observe the effects of mild to moderate COVID-19 on ANS function over a follow-up period, including before and after six weeks of immunization, in young adults. Secondarily, we aim to identify PA behavior over that period and analyze possible correlations with ANS modulation.

2. Materials and Methods

2.1. Ethical Approval

The study was conducted based on the ethical standards described in the Declaration of Helsinki. All participants were notified about the study purpose and protocols and provided written informed consent to all protocols. The Ethical Institutional Review Board approved the study (approval number: 38701820.0.0000.5402).

2.2. Study Design

This is an observational prospective case-control study, part of a broader project: FIT-COVID Study [28]. The study was previously registered at the Brazilian Clinical Trials Registry (registration number: RBR-5dqvkv3). All reports followed the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) guidelines [29].
Individuals previously infected with COVID-19 (PCOV) were initially recruited before immunization to participate in the study through local media (television, radio) and social media via electronic access to participants’ database of the Municipal Health Secretariat of Presidente Prudente São Paulo, Brazil (231,953 inhabitants; human development index, 0.806; 23,657 confirmed COVID-19 cases, moving average of 135 cases (May 2021)). An age-matched healthy control group (CT) that was negative for COVID-19 was also recruited. Six weeks after full immunization, individuals in both groups were invited to repeat assessments (Figure 1).
Inclusion criteria were male and female subjects aged 20–40 years with a diagnosis of mild or moderate clinical COVID-19, with a previous positive polymerase chain reaction (PCR) test. Mild cases were considered when patients presented slight clinical symptoms and absence of imaging findings of pneumonia. Moderate cases were considered when patients reported fever or respiratory symptoms [30]. Participants were recruited after a minimum of 15 and a maximum of 120 days of diagnosis by positive PCR test [31]. To screen for confirmed or probable previous SARS-CoV-2 infection for the control group, a lateral flow test for SARS-CoV-2 Immunoglobin G (IgG) and Immunoglobin M (IgM) antibodies was conducted using amplified chemiluminescence and chemiluminescence serological methods, respectively. Subjects in both groups were eligible for a follow-up evaluation if they received any type of SARS-CoV-2 immunization after initial evaluations.
We excluded subjects that presented severe clinical manifestations of COVID-19, including respiratory distress and a respiratory rate > 30 times per minute, fingertip blood oxygen saturation < 93% at rest, and partial arterial oxygen pressure (PaO2)/fraction of inspiration oxygen (FiO2) < 300 mmHg. Subjects that reported critical conditions, including respiratory failure requiring mechanical ventilation, shock, and other organ failure, or requiring intense care treatment, were also excluded [30]. Finally, we also excluded individuals that participated in rehabilitation programs after COVID-19 or that presented any chronic noncommunicable diseases, smokers, history of continous drug use, medications, such as anti-inflammatory drugs, antibiotics, or other drugs known for their impact on the ANS.

2.3. Evaluations

Initially, general anamnesis was taken, including sociodemographic characterisitics and self-rated health and medical history (family comorbidities, cardiopulmonary symptoms, and medication use). Participants were then assessed regarding their body mass index (BMI) [32] and PA level [33]. Symptoms that emerged during acute SARS-CoV-2 infection and persistent symptoms were also evaluated [30]. Subjects returned for a follow-up evaluation six weeks after complete SARS-CoV-2 immunization in which assessments were repeated [28]. The number of days between baseline (before immunization) and follow-up (after immunization) assessments was 161.60 ± 45.94 days for CT and 168.42 ± 24.26 days for PCOV (approximately five months).

2.3.1. Body Mass Index

Body mass index (BMI) was calculated according to previous literature [32]. We considered BMI as the ratio of weight (kg) to height (m) squared. Weight was evaluated using an electronic and calibrated scale (Kratos-Cas, São Paulo, Brazil) where subjects wore light clothes and were barefoot. Height was evaluated using a portable anthropometer (Kratos-Cas, São Paulo, Brazil) [32].

2.3.2. Heart Rate Variability

ANS function was measured through heart rate variability (HRV), which is considered a simple, dependable, and noninvasive method [34]. For this evaluation, participants were asked to attend an outpatient clinic and were instructed to arrive in a fasted state, having abstained from exercise, caffeine, chocolate, and alcohol for at least 24 h before evaluation and ≥4 h after a snack or light meal. Evaluations were performed in a silent space, with temperature of ~23 °C. HRV analysis was performed in the morning to avoid extraneous influences of circadian changes [34].
Heart rate was recorded beat to beat to evaluate cardiac autonomic modulation. We used a cardio-frequency meter (Polar RS800CX, Polar Electro, Kempele, Finland) at a 1 kHz sampling rate. Participants were equipped with a chest strap and monitor and remained at rest with spontaneous breathing for 25 min [35]. HRV was performed on 256 consecutive intervals between succesive heartbeats (RR intervals) from the most stable segment on the tachogram. Only series with <5% error were considered suitable for analysis. Kubios HRV® software (Biosignal Analysis and Medical Image Group, Department of Physics, University of Kuopio, Finland) was used to complete the HRV analysis [36,37].
The HRV was assessed in both time and frequency domains using Kubios HRV®. For the time domain, the mean RR intervals (reflecting global variability) were used. Additionally, the indexes presented in Table 1 were calculated.

2.3.3. Physical Activity Level

PA level was measured using a triaxial accelerometer (GT3X+; ActiGraph, LLC, Pensacola, FL, USA). Participants were instructed to wear the accelerometer above the waist for seven consecutive days during waking hours. A minimum of four days with at least 10 h per day was considered valid accelerometer data. Participants were instructed not to use the accelerometer while bathing, sleeping, or performing water activities. Every morning, a researcher sent a WhatsApp message reminding the participant to use the accelerometer.
Non-wear periods were defined as time intervals of at least 60 consecutive minutes of zero counts, with an activity interruption allowance of 0–100 counts per minute lasting a maximum of two consecutive minutes [39]. Counts per minute were calculated using the sum of the total activity counts in the vertical axis divided by the valid number of days. Sedentary time was defined as values <100 counts per minute, light PA as values between 100 and 2019 counts per minute, and moderate–vigourous PA (MVPA) as values >2020 counts/minute. Data were processed using ActLife software (version 6.9.2, ActiGraph Pensacola, FL, USA) [33].

2.3.4. Statistical Analysis

Statistical procedures were performed using GraphPad Prism 9 (version 9.3.1; GraphPad Software, Dotmatics, San Diego, CA, USA). Data distribution was analyzed using the Shapiro–Wilk test. For primary analysis of intergroup comparisons between the post-COVID-19 (PCOV) and control groups (CT), the unpaired t-test or Mann–Whitney U-test was used according to data distribution. Additionally, the Welch’s correction was used for comparison of variables with different standard deviations. For intragroup paired analysis, paired t-test or Wilcoxon test was performed. The Chi-square test was used to assess for differences within the categorical variables of PA. The Fisher exact test was performed to assess for differences between sex distribution.
Effect size (ES) was calculated using the difference between moments to determine the magnitude of differences over time. ES was represented by Cohen’s d values and was classified as d (0.01) = very small, d (0.2) = small, d (0.5) = medium, d (0.8) = large, d (1.2) = very large, and d (2.0) = huge [40].
Secondarily, the data were analyzed using analysis of covariance (ANCOVA) to better understand the influences of PA on the change in HRV indexes [41]. The dependent variable was represented by each HRV index and MVPA was used as the covariate in this multivariable regression. Data were presented with mean and standard deviation or median and respective quartiles (Q1–Q3) according to data distribution. Statistical significance was considered p < 0.05.

3. Results

In terms of participants, 239 control and 154 post-COVID subjects were assessed for eligibility for this study. Further, 92 control and 33 post-COVID subjects were deemed eligible to participate in this study and 268 individuals did not meet the inclusion criteria. Moreover, 57 of these subjects were evaluated at baseline. After 19 exclusions due to errors in HRV recordings, 38 subjects with complete data were included in the baseline analysis, where 18 subjects were in the CT and 20 in the PCOV group. Subjects then returned six weeks after receiving full COVID-19 immunization. At this stage, 20 participants were excluded due to declining a follow-up assessment and one individual was excluded due to SARS-CoV-2 reinfection. The number of individuals at each stage of the study is presented in Figure 2.
No differences in sex distribution were observed between the groups (p = 0.1119). Descriptive characteristics of age, weight, height, and body mass index were uniformly distributed and not different between groups, as indicated in Table 2. Participants in CT returned for follow-up evaluation 161.60 ± 45.94 days after baseline evaluation and PCOV returned after 168.42 ± 24.26 days (approximately five months). In CT, 30% received AstraZeneca immunization, 30% received CoronaVac, and 40% received Pfizer. In PCOV, 57% received CoronaVac immunization and 43% received Pfizer. Additionally, patients in PCOV were evaluated at baseline 50 ± 32.54 days after testing positive for SARS-CoV-2, classified as mild to moderate COVID-19 cases [25].
There were no statistical differences (p > 0.05) in PA between CT and PCOV groups (intergroup analysis) and between baseline and follow-up moments (intragroup analysis).
Results for cross-sectional analysis with comparisons between CT and PCOV at baseline are presented in a previous study by the FIT-COVID research group [9].
Figure 3 presents the intragroup analysis. There was no difference between baseline and follow-up moments within CT or PCOV groups (p > 0.05) regarding SNS activity (Figure 3). Parasympathetic nervous system (PNS) activity increased between baseline and follow-up moments, reflected by mean RR (p = 0.0312) and pNN50 (p = 0.0312) for PCOV group (Figure 3). Mean RR increased from 734.50 ± 131.80 ms at the baseline moment to 794.90 ± 105.60 ms at the follow-up moment and pNN50 increased 4.53 ± 4.68% to 11.74 ± 8.81%. No differences between moments were observed in CT (p > 0.05). The intragroup analysis between moments also revealed a difference in global variability, reflected by RR triangular index (7.91 ± 2.34 at baseline to 9.60 ± 2.17 at follow-up; p = 0.0312; Figure 3).
Table 3 and Figure S1 portrays primary analysis of the present study, exhibiting intergroup comparisons between HRV indexes in control (CT) and PCOV groups. (Figure S1). We observed that the PCOV group presents significant reductions in sympathetic activity over time when compared to CT demonstrated by mean HR (p = 0.0088) and SNS index (p = 0.0068; Table 3; Figure S1). PCOV also presented with a significant increase in PNS activity over time demonstrated by mean RR (−44.54 ± 32.38 vs. 60.36 ± 55.35 ms; p = 0.0097) and PNS index (−0.32 ± 0.20 vs. 0.54 ± 0.35; p = 0.0091) when compared to CT (Table 3; Figure S1). No intergroup differences were observed in global variability.
Intergroup comparisons between PCOV and CT at the follow-up moment revealed that there were no significant differences between the groups in any HRV index (p > 0.05).
ES was determined using Cohen’s d values [40]. Variables that yielded significance in the inter-group analysis had ES reported. Within SNS activity indexes, mean HR (ES = 2.56) and SNS index (ES = 2.56) both had very large ESs. PNS activity indexes of mean RR (ES = 2.31) and PNS index (ES = 3.02) also yielded very large ESs. No global variability indexes yielded significance during the inter-group analysis. Lastly, no PA variables presented significance in either the inter- or intra-group analysis.
The secondary analysis, including MVPA in the multivariable regression analysis model, showed that the significant group differences in mean HR, SNS index, mean RR, and PNS index were maintained (p < 0.05) (Table 3).

4. Discussion

This study observed the effects of mild to moderate COVID-19 on the ANS in young adults, before and a minimum of six-weeks after complete SARS-CoV-2 immunization. The relationship between PA and the changes in ANS were also analyzed. The primary finding of this study is that autonomic function was improved in young adults who were infected by SARS-CoV-2 after approximately five months of follow-up. This change was characterized by significant decreases in SNS indexes (mean HR and SNS index), while also having significant increases in PNS HRV indexes (mean RR and PNS index). Additionally, when both groups were compared at the follow-up moment, our results revealed that the young adults infected by SARS-CoV-2 presented similar autonomic function when compared to CT. Secondarily, no differences were observed in PA between PCOV and CT groups over the follow-up period. The changes in HRV data were maintained even after statistical adjustments using PA levels.
To our knowledge, this study is the first to observe improvements in autonomic function after mild to moderate post-COVID-19 infection in a young adult population within 164.41 ± 91.47 days. Much of the existing literature focuses on older populations who present with greater risk of negative outcomes [3,4,42] or those infected with severe forms of COVID-19 [8]. Additionally, most of the research focuses on cross-sectional data identifying autonomic dysfunctions associated with post-COVID-19, rather than progression and recovery, making it increasingly difficult to understand ANS modulation over time and the relationship with PA.
Previous work by Freire et al. [9], as part of the cross-sectional data from the Fit-COVID study, revealed the presence of autonomic dysfunction shortly after mild to moderate post-COVID-19 in young adults (represented as baseline comparisons in the present study). The study reported that even in mild and moderate infection, young adults who had COVID-19 had greater sympathetic activity (increased levels of the stress index), decreased parasympathetic activity (lower values of RMSSD and SD1 indices), and global variability (reflected through the SDNN, TINN, and SD2 indexes) when compared to individuals who were not infected. Moreover, in participants who were overweight and obese and/or physically inactive, cardiac autonomic dysregulation was more prominent. [9].
Increases in the SNS are associated with a systemic inflammatory condition characterized by perfusion of inflammatory cytokines in bloodstream and other biomarkers associated with inflammation [3,8,11,13]. Additionally, increases in SNS activity are associated with the secretion of catecholamines, which increases metabolism and cardiac activity, cumulatively increasing cardiac stress [43]. Increases in biomarkers indicating oxidative stress have also been observed in association with COVID-19 and flu-like infections [3,12,44]. The PNS’s influence on anti-inflammatory and restorative processes is also inhibited [45]. Previous studies did not identify the exact mechanism behind inhibition of PNS after COVID-19, but likely a function of the reciprocal nature between the SNS and PNS. Therefore, the increase in oxidative stress and subsequent cytokine storm in tandem with reductions in restorative processes that are associated with SARS-CoV-2 infection could explain the alterations observed at baseline.
It is important to highlight that age plays an important role in ANS regulation. Our study observed young adults without any known chronic diseases. In young adults, protective factors can be related to better ANS regulation and might help explain the responses of ANS after virus infection observed in our study. Healthy young adults can present lower levels of arterial and ventricular stiffening, better myocardial contractility, and preserved organ innervation when compared to elderly individuals [46,47]. Additionally, improved sensitivity of vagal reflexes can lead to greater HRV in this population [16]. Therefore, a combination of these aspects could have influenced the improvements observed in the PCOV group over time. We believe that our results support the use of HRV monitoring as an important, reliable, and accessible method to observe alterations in autonomic function. Deviations in autonomic regulation are related to distancing from homeostasis and are present in multiple diseases, both those that directly afflict the nervous system and those afflicting other organs, where they trigger or enhance pathological symptoms [18,48].
This study did not observe changes in PA over time in either group, nor were PA levels between PCOV and CT different from one another. Although there is a lack of statistical difference to explain observations in this study, which may have been a result of an underpowered sample, PA levels still may have had an influence on the improvements in ANS observed in the PCOV group as median PA levels increased from baseline to follow-up. PA aids in ANS modulation, as the peripheral stress induced by skeletal muscle contraction or mechanical stress on organs sends afferent signals to the central nervous system, thus, inducing catecholamine (epinephrine and norepinephrine) release to alter ANS function to meet metabolic demands of exercise [24]. The release of epinephrine and norepinephrine increases SNS activity, consequently promoting PNS activity during recovery from exercise [49].
The World Health Organization recommends that adults should participate in 150 min of moderate PA per week or 75 min of vigorous PA to observe physiological benefits that can protect against COVID-19 outcomes and reduce the risk of many diseases [50]. At baseline, CT had a median moderate PA of 17.90 min per day (min/day), which equates to 125.30 min per week (min/wk) of moderate PA, while PCOV had a median of 13.36 min/day, equating to 93.52 min/wk of moderate PA. Baseline assessments revealed that CT and PCOV did not meet recommendations for moderate or vigorous PA. These results may be explained by modifications in daily behaviors that were a result of sanitary efforts. Many countries employed lockdowns and restrictions to public areas to avoid the continued spread of COVID-19, which increased barriers to achieve ideal levels of physical activity [51,52].
These barriers were discussed in observations by Huber et. al. [52], in which PA was significantly reduced in young adults during the pandemic. Studies have also found that young adults did not return to pre-pandemic PA levels, which were characterized by a combination of increased sedentary activity and a lower frequency of PA bouts [51]. By the time of the follow-up assessment, PCOV increased moderate PA levels to satisfy PA recommendations from a median of 27.86 min/day (195.02 min/wk) of moderate PA. The increase in PA can likely be explained by the relaxation of social restrictions driven by increasing SARS-CoV-2 immunization rates [53], allowing participants in this study to have easier access to public spaces, such as gyms and community parks.
PA is important in the modulation of immune responses, where increased levels of PA are associated with improved cytokine responses, thus, reducing the severity of COVID-19 disease [54]. PA helps to modulate the release of anti-inflammatory factors and, thus, may aid in immune response to SARS-CoV-2 infection [55]. Specific immune cell activities are regulated by the ANS, where dysfunction in the regulatory system can lead to inadequate immune stimulation [12]. In response, there is an accumulation of reactive nitrogenous species (RNS) and reactive oxygen species (ROS) at the site of cell damage, where persistent inflammation can lead to cellular oxidative damage and damage to the infected host [56]. PA can aid in the immune response to infection through increasing antioxidant capacity and increasing stimulation of the T-helper 2 (TH2) cell pathway [12]. PA increases the level of oxidative stress in skeletal muscle as a product of metabolic processes, where a negative feedback mechanism increases levels of antioxidants as a defense to the presence of ROS and RNS [55,57]. The tight relationship between PA level and immune response illustrates a possible mechanism as to why we observed improved ANS function in this young adult population after SARS-CoV-2 infection [22,54,55,58].
This study presents significant clinical relevance, as it provides further information on the relationship of PA behaviors on recovery after mild to moderate COVID-19 infection in a young adult population. We were able to observe alterations in ANS activity after a five-month follow-up period in a young adult population, highlighted by reductions in SNS activity and increases in PNS activity. Although we hypothesize that this observation could be related to PA levels, our experimental design does not allow us to establish causal effects in this matter. Therefore, we strongly suggest that future studies further investigate this aspect, investigating ANS responses across different ranges of PA, which was not possible in our study due to the small sample size.
This evidence better informs health professionals on non-medicinal approaches to reduce the risk and severity of COVID-19 outcomes by using PA as an alternative medicine. Future research should continue to observe the effects of common population behaviors (which may be unique to age, location of residence, or occupation) on PA levels and the influence of these factors on COVID-19 outcomes. Further research may also look to observe the relationships of frequency, intensity, and duration of PA necessary to improve modulation of cytokine immune responses.
Since we evaluated young adults before and after SARS-CoV-2 immunization, the interpretation of ANS response over time should also consider the potential effects of the vaccine. Although the time point for testing the long-term effects of the vaccine is still relatively short, previous studies comprehensively observed the safety and effectiveness of the COVID-19 vaccine [59]. Reports on the effects of immunization on ANS regulation are still scarce and isolated. Limited case reports previously presented occurrences of post-COVID-19 vaccine-associated postural orthostatic tachycardia syndrome (POTS), a multifactorial and often debilitating type of autonomic dysfunction [60,61,62]. These alterations were observed after immunization (four hours to seven days) and a possible mechanism could be attributed to sinus tachycardia and cerebral hypoperfusion secondary to hypocapnia observed in these cases. This mechanism could result in ischemic hypoxia of the carotid body, leading to chemoreflex activation and sympathetic activation [60]. It is important to highlight that these findings are based on very small samples and the pathophysiological role of the patient’s prior COVID-19 infection still requires further investigation. The reports also state that these observations cannot imply causality of POTS possibly linked to the COVID-19 vaccine [60].
In our study, the aspect of immunization in our study was used as a controlling factor to avoid bias (all individuals were vaccinated upon follow-up). Therefore, the absence of a non-vaccinated group in our design limits our ability to differentiate and account for the specific effects of immunization on ANS. Due to ethical aspects and the high rates of immunization after the availability of immunization in Brazil (over 88% of the population received at least one shot of immunization), it was difficult to incorporate a non-vaccinated group. Future investigations should consider the investigation of ANS function in non-vaccinated individuals to further advance knowledge in this aspect.
Limitations from this study include the loss of follow-up information from baseline to follow-up. CT had 45% and PCOV had 65% loss of participants between observational moments, which may be attributed to the general difficulty in subject compliance for young adults. This high loss of follow-up increases the likelihood of an underpowered sample, which would undermine the statistical findings or explain the lack of statistical differences (50). The effect of the underpowered sample on the statistical findings in this study can only be determined by referencing other studies employing similar methods. Lastly, although we excluded individuals that required intensive care from our sample, subjects may have experienced symptoms, such as chest pain or severe myalgias, indicative of myocarditis or pericarditis, while still not requiring intensive care. This aspect was not controlled in our study and needs to be considered for the interpretation of our results.
Previous research throughout the literature has employed HRV methods that require prolonged monitoring, are subjective in nature, or require expensive equipment that is not accessible to the general population or small-practice clinics. The methods employed in this study are less cumbersome for both healthcare professionals and the patient being monitored, as HRV monitoring can be completed within a single, relatively short visit. Observing changes in ANS activity may be helpful to health professionals to detect signs of infection, track viral progression, and observe dysfunctions caused by the virus, which can aid in the development of individualized intervention strategies that are better informed by HRV data.
Overall, this study observed ANS improvements over a follow-up of an approximately five-month period that may be related to recovery after SARS-CoV-2 infection. To our knowledge, this study is the first of its kind to investigate these outcomes within a four- to six-month period. Lastly, the methods for monitoring ANS function could be used by healthcare professionals to observe the progression of SARS-CoV-2 infection to form individualized non-medication intervention strategies.

5. Conclusions

We observed improvements in ANS function after mild to moderate SARS-CoV-2 infection in young adults vaccinated over a five-month period. Additionally, no changes in physical activity levels were observed over the follow-up period.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijerph20032251/s1, Figure S1: Scatterplot of HRV indexes for unpaired analysis according to group expressed as the mean difference between follow-up and baseline moments. (A) Sympathetic nervous system activity; (B) parasympathetic nervous system activity. * Statistical difference between groups; mean HR: average heart rate measured in beats per minute; SNS Index: sympathetic nervous system index; mean RR: time between.

Author Contributions

Conceptualization: F.S.L., T.P., M.-J.C.-E.-S. and B.S.d.A.S.; methodology: A.P.C.F.F., A.E.v.A.M., A.C., B.S.d.A.S., V.R.D.S. and O.M.J.; software: D.G.D.C.; formal analysis: S.A., A.P.C.F.F., B.S.d.A.S. and V.R.D.S.; investigation: F.S.L., T.P., M.-J.C.-E.-S., A.C., O.M.J. and R.A.P.; data curation: A.P.C.F.F., A.E.v.A.M. and B.S.d.A.S.; writing—original draft preparation: S.A., A.P.C.F.F., F.S.L. and A.C.; writing—review and editing: S.A., A.P.C.F.F. and F.S.L.; visual-ization: A.P.C.F.F., F.S.L., A.E.v.A.M., T.P., M.-J.C.-E.-S., A.C., D.G.D.C., O.M.J., R.A.P. and B.S.d.A.S.; supervision: A.P.C.F.F., F.S.L. and B.S.d.A.S.; project administration: F.S.L. and B.S.d.A.S.; funding acquisition: F.S.L. All authors have read and agreed to the published version of the manuscript.

Funding

National Council for Scientific and Technological Development (CNPq)—Brazil (Process number: 150760/2022-1) and the Coordination for the Improvement of Higher Education Personnel (CAPES)—Brazil (code 001). F.S.L. was granted a research scholarship (PQ2) from the CNPq (302714/2021-9). São Paulo Research Foundation (FAPESP)/Brazil (Process Number, 2021/11932-8).

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki and approved by the Institutional Review Board of Sao Paulo State University (protocol code: 38701820.0.0000.5402; date of approval: 19 March 2021).

Informed Consent Statement

Informed consent was obtained from all participants involved in the study.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to ethical aspects.

Acknowledgments

We thank the National Council for Scientific and Technological Development (CNPq)—Brazil and the Coordination for the Improvement of Higher Education Personnel (CAPES)—Brazil (code 001) and the São Paulo Research Foundation (FAPESP)/Brazil. Additionally, we thank all authors and contributors to the international collaboration of the FIT-COVID-19 Study.

Conflicts of Interest

None of the authors have any conflict of interest.

References

  1. WHO. Coronavirus (COVID-19) Dashboard 2022. Available online: https://covid19.who.int/ (accessed on 22 November 2022).
  2. Dani, M.; Dirksen, A.; Taraborrelli, P.; Torocastro, M.; Panagopoulos, D.; Sutton, R.; Lim, P.B. Autonomic dysfunction in ‘long COVID’: Rationale, physiology and management strategies. Clin. Med. 2021, 21, e63. [Google Scholar] [CrossRef]
  3. Guilmot, A.; Maldonado Slootjes, S.; Sellimi, A.; Bronchain, M.; Hanseeuw, B.; Belkhir, L.; Yombi, J.C.; De Greef, J.; Pothen, L.; Yildiz, H.; et al. Immune-mediated neurological syndromes in SARS-CoV-2-infected patients. J. Neurol. 2021, 268, 751–757. [Google Scholar] [CrossRef] [PubMed]
  4. Buoite Stella, A.; Furlanis, G.; Frezza, N.A.; Valentinotti, R.; Ajcevic, M.; Manganotti, P. Autonomic dysfunction in post-COVID patients with and witfhout neurological symptoms: A prospective multidomain observational study. J. Neurol. 2022, 269, 587–596. [Google Scholar] [CrossRef]
  5. Pan, Y.; Yu, Z.; Yuan, Y.; Han, J.; Wang, Z.; Chen, H.; Wang, S.; Wang, Z.; Hu, H.; Zhou, L.; et al. Alteration of Autonomic Nervous System Is Associated With Severity and Outcomes in Patients With COVID-19. Front. Physiol. 2021, 12, 630038. [Google Scholar] [CrossRef] [PubMed]
  6. Asarcikli, L.D.; Hayiroglu, M.; Osken, A.; Keskin, K.; Kolak, Z.; Aksu, T. Heart rate variability and cardiac autonomic functions in post-COVID period. J. Interv. Card. Electrophysiol. 2022, 63, 715–721. [Google Scholar] [CrossRef] [PubMed]
  7. Marques, K.C.; Silva, C.C.; Trindade, S.D.S.; Santos, M.C.S.; Rocha, R.S.B.; Vasconcelos, P.F.D.C.; Quaresma, J.A.S.; Falcão, L.F.M. Reduction of Cardiac Autonomic Modulation and Increased Sympathetic Activity by Heart Rate Variability in Patients With Long COVID. Front. Cardiovasc. Med. 2022, 9, 862001. [Google Scholar] [CrossRef] [PubMed]
  8. Stute, N.L.; Stickford, J.L.; Province, V.M.; Augenreich, M.A.; Ratchford, S.M.; Stickford, A.S.L. COVID-19 is getting on our nerves: Sympathetic neural activity and haemodynamics in young adults recovering from SARS-CoV-2. J. Physiol. 2021, 599, 4269–4285. [Google Scholar] [CrossRef] [PubMed]
  9. Freire, A.P.C.F.; Lira, F.S.; Morano, A.E.V.A.; Pereira, T.; Coelho-E-Silva, M.-J.; Caseiro, A.; Christofaro, D.G.D.; Marchioto Júnior, O.; Dorneles, G.P.; Minuzzi, L.G.; et al. Role of Body Mass and Physical Activity in Autonomic Function Modulation on Post-COVID-19 Condition: An Observational Subanalysis of Fit-COVID Study. Int. J. Environ. Res. Public Health 2022, 19, 2457. [Google Scholar] [CrossRef]
  10. Tenforde, M.W.; Kim, S.S.; Lindsell, C.J.; Billig Rose, E.; Shapiro, N.I.; Files, D.C.; Gibbs, K.W.; Erickson, H.L.; Steingrub, J.S.; Smithline, H.A.; et al. Symptom Duration and Risk Factors for Delayed Return to Usual Health among Outpatients with COVID-19 in a Multistate Health Care Systems Network—United States, March–June 2020; Center for Disease Control and Prevention Stacks: New York, NY, USA, 2020. [Google Scholar]
  11. Porzionato, A.; Emmi, A.; Barbon, S.; Boscolo-Berto, R.; Stecco, C.; Stocco, E.; Macchi, V.; De Caro, R. Sympathetic activation: A potential link between comorbidities and COVID-19. FEBS J. 2020, 287, 3681–3688. [Google Scholar] [CrossRef]
  12. Daniela, M.; Catalina, L.; Ilie, O.; Paula, M.; Daniel-Andrei, I.; Ioana, B. Effects of Exercise Training on the Autonomic Nervous System with a Focus on Anti-Inflammatory and Antioxidants Effects. Antioxidants 2022, 11, 350. [Google Scholar] [CrossRef]
  13. Konig, M.F.; Powell, M.; Staedtke, V.; Bai, R.Y.; Thomas, D.L.; Fischer, N.; Huq, S.; Khalafallah, A.M.; Koenecke, A.; Xiong, R.; et al. Preventing cytokine storm syndrome in COVID-19 using α-1 adrenergic receptor antagonists. J. Clin. Investig. 2020, 130, 3345–3347. [Google Scholar] [CrossRef]
  14. Tabernero, E.; Ruiz, L.A.; España, P.P.; Méndez, R.; Serrano, L.; Santos, B.; Uranga, A.; González, P.; Garcia, P.; Torres, A.; et al. COVID-19 in young and middle-aged adults: Predictors of poor outcome and clinical differences. Infection 2022, 50, 179–189. [Google Scholar] [CrossRef]
  15. Skow, R.J.; Garza, N.A.; Nandadeva, D.; Stephens, B.Y.; Wright, A.N.; Grotle, A.K.; Young, B.E.; Fadel, P.J. Impact of COVID-19 on cardiac autonomic function in healthy young adults: Potential role of symptomatology and time since diagnosis. Am. J. Physiol. Heart Circ. Physiol. 2022, 323, H1206–h11. [Google Scholar] [CrossRef] [PubMed]
  16. Schumann, A.; Bär, K.J. Autonomic aging—A dataset to quantify changes of cardiovascular autonomic function during healthy aging. Sci. Data. 2022, 9, 95. [Google Scholar] [CrossRef]
  17. Libby, P.; Lüscher, T. COVID-19 is, in the end, an endothelial disease. Eur. Heart J. 2020, 41, 3038–3044. [Google Scholar] [CrossRef]
  18. Goldberger, J.J.; Arora, R.; Buckley, U.; Shivkumar, K. Autonomic Nervous System Dysfunction: JACC Focus Seminar. J. Am. Coll. Cardiol. 2019, 73, 1189–1206. [Google Scholar] [CrossRef] [PubMed]
  19. Brunner-La Rocca, H.P.; Esler, M.D.; Jennings, G.L.; Kaye, D.M. Effect of cardiac sympathetic nervous activity on mode of death in congestive heart failure. Eur. Heart J. 2001, 22, 1136–1143. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Wenner, M.M. Sympathetic activation in chronic anxiety: Not just at the “height” of stress. Editorial Focus on “Relative burst amplitude of muscle sympathetic nerve activity is an indicator of altered sympathetic outflow in chronic anxiety”. J. Neurophysiol. 2018, 120, 7–8. [Google Scholar] [CrossRef]
  21. Sallis, R.; Young, D.R.; Tartof, S.Y.; Sallis, J.F.; Sall, J.; Li, Q.; Smith, G.N.; Cohen, D.A. Physical inactivity is associated with a higher risk for severe COVID-19 outcomes: A study in 48 440 adult patients. Br. J. Sports Med. 2021, 55, 1099–1105. [Google Scholar] [CrossRef]
  22. Després, J.P. Severe COVID-19 outcomes—The role of physical activity. Nat. Rev. Endocrinol. 2021, 17, 451–452. [Google Scholar] [CrossRef]
  23. Hautala, A.J.; Kiviniemi, A.M.; Tulppo, M.P. Individual responses to aerobic exercise: The role of the autonomic nervous system. Neurosci. Biobehav. Rev. 2009, 33, 107–115. [Google Scholar] [CrossRef] [PubMed]
  24. Fu, Q.; Levine, B.D. Exercise and the autonomic nervous system. Handb. Clin. Neurol. 2013, 117, 147–160. [Google Scholar] [PubMed]
  25. Pedersen, B.K.; Saltin, B. Exercise as medicine—Evidence for prescribing exercise as therapy in 26 different chronic diseases. Scand. J. Med. Sci. Sports 2015, 25 (Suppl. S3), 1–72. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Fogelholm, M. Physical activity, fitness and fatness: Relations to mortality, morbidity and disease risk factors. A systematic review. Obes. Rev. 2010, 11, 202–221. [Google Scholar] [CrossRef] [PubMed]
  27. Leidman, E.; Duca, L.M.; Omura, J.D.; Proia, K.; Stephens, J.W.; Sauber-Schatz, E.K. COVID-19 Trends Among Persons Aged 0–24 Years—United States, March 1–12 December 2020; Morbidity and Mortality Weekly Report; Centers for Disease Control and Prevention: New York, NY, USA, 2021; pp. 88–94. [Google Scholar]
  28. Lira, F.S.; Pereira, T.; Guerra Minuzzi, L.; Figueiredo, C.; Olean-Oliveira, T.; Figueira Freire, A.P.C.; Coelho-e-Silva, M.J.; Caseiro, A.; Thomatieli-Santos, R.V.; Dos Santos, V.R.; et al. Modulatory Effects of Physical Activity Levels on Immune Responses and General Clinical Functions in Adult Patients with Mild to Moderate SARS-CoV-2 Infections—A Protocol for an Observational Prospective Follow-Up Investigation: Fit-COVID-19 Study. Int. J. Environ. Res. Public Health 2021, 18, 13249. [Google Scholar] [CrossRef]
  29. Von Elm, E.; Altman, D.G.; Egger, M.; Pocock, S.J.; Gøtzsche, P.C.; Vandenbroucke, J.P.; Strobe Initiative. The Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) Statement: Guidelines for reporting observational studies. Int. J. Surg. 2014, 12, 1495–1499. [Google Scholar] [CrossRef] [Green Version]
  30. Zhang, J.; Wang, X.; Jia, X.; Li, J.; Hu, K.; Chen, G.; Wei, J.; Gong, Z.; Zhou, C.; Yu, H.; et al. Risk factors for disease severity, unimprovement, and mortality in COVID-19 patients in Wuhan, China. Clin. Microbiol. Infect. 2020, 26, 767–772. [Google Scholar] [CrossRef]
  31. World Health Organization. Public Health Surveillance for COVID-19: Interim Guidance. 22 July 2022. Available online: https://www.who.int/publications/i/item/WHO-2019-nCoV-SurveillanceGuidance-2022.2 (accessed on 22 November 2022).
  32. Nuttall, F.Q. Body Mass Index: Obesity, BMI, and Health: A Critical Review. Nutr. Today 2015, 50, 117–128. [Google Scholar] [CrossRef] [Green Version]
  33. Troiano, R.P.; Berrigan, D.; Dodd, K.W.; Mâsse, L.C.; Tilert, T.; McDowell, M. Physical activity in the United States measured by accelerometer. Med. Sci. Sports Exerc. 2008, 40, 181–188. [Google Scholar] [CrossRef]
  34. Malik, M. Heart rate variability: Standards of measurement, physiological interpretation and clinical use. Task Force of the European Society of Cardiology and the North American Society of Pacing and Electrophysiology. Circulation 1996, 93, 1043–1065. [Google Scholar] [CrossRef]
  35. Laborde, S.; Mosley, E.; Thayer, J.F. Heart Rate Variability and Cardiac Vagal Tone in Psychophysiological Research—Recommendations for Experiment Planning, Data Analysis, and Data Reporting. Front. Psychol. 2017, 8, 213. [Google Scholar] [CrossRef] [Green Version]
  36. Vanderlei, L.C.; Pastre, C.M.; Hoshi, R.A.; Carvalho, T.D.; Godoy, M.F. Basic notions of heart rate variability and its clinical applicability. Rev. Bras. Cir. Cardiovasc. 2009, 24, 205–217. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Kubios. HRV in Evaluating ANS Function. 2022. Available online: https://www.kubios.com/hrv-ans-function/ (accessed on 22 November 2022).
  38. Shaffer, F.; Ginsberg, J.P. An Overview of Heart Rate Variability Metrics and Norms. Front. Public Health. 2017, 5, 258. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Matthews, C.E.; Chen, K.Y.; Freedson, P.S.; Buchowski, M.S.; Beech, B.M.; Pate, R.R.; Troiano, R.P. Amount of time spent in sedentary behaviors in the United States, 2003–2004. Am. J. Epidemiol. 2008, 167, 875–881. [Google Scholar] [CrossRef] [Green Version]
  40. Sawilowsky, S.S. New Effect Size Rules of Thumb. J. Mod. Appl. Stat. Methods 2009, 8, 26. [Google Scholar] [CrossRef]
  41. Khammar, A.; Yarahmadi, M.; Madadizadeh, F. What Is Analysis of Covariance (ANCOVA) and How to Correctly Report Its Results in Medical Research? Iran J. Public Health. 2020, 49, 1016–1017. [Google Scholar] [CrossRef]
  42. Shouman, K.; Vanichkachorn, G.; Cheshire, W.P.; Suarez, M.D.; Shelly, S.; Lamotte, G.J.; Sandroni, P.; Benarroch, E.E.; Berini, S.E.; Cutsforth-Gregory, J.K.; et al. Autonomic dysfunction following COVID-19 infection: An early experience. Clin. Auton. Res. 2021, 31, 385–394. [Google Scholar] [CrossRef]
  43. Borovac, J.A.; D’Amario, D.; Bozic, J.; Glavas, D. Sympathetic nervous system activation and heart failure: Current state of evidence and the pathophysiology in the light of novel biomarkers. World J. Cardiol. 2020, 12, 373–408. [Google Scholar] [CrossRef] [PubMed]
  44. Choi, A.M.; Knobil, K.; Otterbein, S.L.; Eastman, D.A.; Jacoby, D.B. Oxidant stress responses in influenza virus pneumonia: Gene expression and transcription factor activation. Am. J. Physiol. 1996, 271, L383–L391. [Google Scholar] [CrossRef]
  45. Hassani, M.; Fathi Jouzdani, A.; Motarjem, S.; Ranjbar, A.; Khansari, N. How COVID-19 can cause autonomic dysfunctions and postural orthostatic syndrome? A Review of mechanisms and evidence. Neurol. Clin. Neurosci. 2021, 9, 434–442. [Google Scholar] [CrossRef]
  46. Boettger, M.K.; Schulz, S.; Berger, S.; Tancer, M.; Yeragani, V.K.; Voss, A.; Bär, K.J. Influence of age on linear and nonlinear measures of autonomic cardiovascular modulation. Ann. Noninvasive Electrocardiol. 2010, 15, 165–174. [Google Scholar] [CrossRef]
  47. Laitinen, T.; Niskanen, L.; Geelen, G.; Länsimies, E.; Hartikainen, J. Age dependency of cardiovascular autonomic responses to head-up tilt in healthy subjects. J. Appl. Physiol. 2004, 96, 2333–2340. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Low, P.A.; Tomalia, V.A.; Park, K.J. Autonomic function tests: Some clinical applications. J. Clin. Neurol. 2013, 9, 1–8. [Google Scholar] [CrossRef] [Green Version]
  49. Vieluf, S.; Hasija, T.; Jakobsmeyer, R.; Schreier, P.J.; Reinsberger, C. Exercise-Induced Changes of Multimodal Interactions Within the Autonomic Nervous Network. Front. Physiol. 2019, 10, 240. [Google Scholar] [CrossRef] [PubMed]
  50. Bull, F.C.; Al-Ansari, S.S.; Biddle, S.; Borodulin, K.; Buman, M.P.; Cardon, G.; Carty, C.; Chaput, J.P.; Chastin, S.; Chou, R. World Health Organization 2020 guidelines on physical activity and sedentary behaviour. Br. J. Sports Med. 2020, 54, 1451–1462. [Google Scholar] [CrossRef] [PubMed]
  51. Ferrara, M.; Langiano, E.; Falese, L.; Diotaiuti, P.; Cortis, C.; De Vito, E. Changes in Physical Activity Levels and Eating Behaviours during the COVID-19 Pandemic: Sociodemographic Analysis in University Students. Int. J. Environ. Res. Public Health 2022, 19, 5550. [Google Scholar] [CrossRef]
  52. Huber, B.C.; Steffen, J.; Schlichtiger, J.; Graupe, T.; Deuster, E.; Strouvelle, V.P.; Fischer, M.R.; Massberg, S.; Brunner, S. Alteration of physical activity during COVID-19 pandemic lockdown in young adults. J. Transl. Med. 2020, 18, 410. [Google Scholar] [CrossRef]
  53. Mathieu, E.; Ritchie, H.; Rodés-Guirao, L.; Appel, C.; Giattino, C.; Hasell, J.; Macdonald, B.; Dattani, S.; Beltekian, D.; Ortiz-Ospin, E.; et al. Brazil: Coronavirus Pandemic Country Profile. Our World in Data 2022. Available online: https://ourworldindata.org/coronavirus/country/brazil (accessed on 22 November 2022).
  54. Mohamed, A.A.; Alawna, M. The effect of aerobic exercise on immune biomarkers and symptoms severity and progression in patients with COVID-19: A randomized control trial. J. Bodyw. Mov. Ther. 2021, 28, 425–432. [Google Scholar] [CrossRef]
  55. Gualano, B.; Lemes, I.R.; Silva, R.P.; Pinto, A.J.; Mazzolani, B.C.; Smaira, F.I.; Sieczkowska, S.M.; Aikawa, N.E.; Pasoto, S.G.; Medeiros-Ribeiro, A.C.; et al. Association between physical activity and immunogenicity of an inactivated virus vaccine against SARS-CoV-2 in patients with autoimmune rheumatic diseases. Brain Behav. Immun. 2022, 101, 49–56. [Google Scholar] [CrossRef]
  56. Montezano, A.C.; Touyz, R.M. Reactive oxygen species and endothelial function--role of nitric oxide synthase uncoupling and Nox family nicotinamide adenine dinucleotide phosphate oxidases. Basic Clin. Pharmacol. Toxicol. 2012, 110, 87–94. [Google Scholar] [CrossRef]
  57. Yen, C.J.; Hung, C.H.; Tsai, W.M.; Cheng, H.C.; Yang, H.L.; Lu, Y.J.; Tsai, K.L. Effect of Exercise Training on Exercise Tolerance and Level of Oxidative Stress for Head and Neck Cancer Patients Following Chemotherapy. Front. Oncol. 2020, 10, 1536. [Google Scholar] [CrossRef] [PubMed]
  58. Yang, J.; Li, X.; He, T.; Ju, F.; Qiu, Y.; Tian, Z. Impact of Physical Activity on COVID-19. Int. J. Environ. Res. Public Health 2022, 19, 14108. [Google Scholar] [CrossRef]
  59. Ling, Y.; Zhong, J.; Luo, J. Safety and effectiveness of SARS-CoV-2 vaccines: A systematic review and meta-analysis. J. Med. Virol. 2021, 93, 6486–6495. [Google Scholar] [CrossRef] [PubMed]
  60. Hermel, M.; Sweeney, M.; Abud, E.; Luskin, K.; Criado, J.P.; Bonakdar, R.; Gray, J.; Ahern, T. COVID-19 Vaccination Might Induce Postural Orthostatic Tachycardia Syndrome: A Case Report. Vaccines 2022, 10, 991. [Google Scholar] [CrossRef] [PubMed]
  61. Reddy, S.; Reddy, S.; Arora, M. A Case of Postural Orthostatic Tachycardia Syndrome Secondary to the Messenger RNA COVID-19 Vaccine. Cureus 2021, 13, e14837. [Google Scholar] [CrossRef]
  62. Karimi Galougahi, K. Autonomic dysfunction post-inoculation with ChAdOx1 nCoV-19 vaccine. Eur. Heart J. Case Rep. 2021, 5, ytab472. [Google Scholar] [CrossRef]
Figure 1. Study design. Evaluations performed on participants at baseline and after a minimum six weeks following complete immunization (follow-up).
Figure 1. Study design. Evaluations performed on participants at baseline and after a minimum six weeks following complete immunization (follow-up).
Ijerph 20 02251 g001
Figure 2. Flow diagram of the study.
Figure 2. Flow diagram of the study.
Ijerph 20 02251 g002
Figure 3. Scatterplot of HRV indexes for paired analysis according to group expressed as the sample mean. (A) Sympathetic nervous system activity; (B) parasympathetic nervous system activity; (C) global variability. * Statistical difference between moments; HR: heart rate; LF: low-frequency component; nu: normalized unit; SNS index: sympathetic nervous system index; RMSSD: root mean square of differences between adjacent normal RR intervals in a time interval; HF: high-frequency component; pNN50: percentage of adjacent RR intervals with a difference in duration >50 ms; SD1: standard deviation of instantaneous beat-to-beat variability; ms: millisecond; PNS index: parasympathetic nervous system index; SDNN: standard deviation of all normal RR intervals recorded in a time interval; RR triangular index: Integral of the density of the RR interval histogram divided by its height; TINN: triangular interpolation of NN interval; LF/HF: ratio of low and high frequency power; SD2: standard deviation of long-term intervals between consecutive heartbeats.
Figure 3. Scatterplot of HRV indexes for paired analysis according to group expressed as the sample mean. (A) Sympathetic nervous system activity; (B) parasympathetic nervous system activity; (C) global variability. * Statistical difference between moments; HR: heart rate; LF: low-frequency component; nu: normalized unit; SNS index: sympathetic nervous system index; RMSSD: root mean square of differences between adjacent normal RR intervals in a time interval; HF: high-frequency component; pNN50: percentage of adjacent RR intervals with a difference in duration >50 ms; SD1: standard deviation of instantaneous beat-to-beat variability; ms: millisecond; PNS index: parasympathetic nervous system index; SDNN: standard deviation of all normal RR intervals recorded in a time interval; RR triangular index: Integral of the density of the RR interval histogram divided by its height; TINN: triangular interpolation of NN interval; LF/HF: ratio of low and high frequency power; SD2: standard deviation of long-term intervals between consecutive heartbeats.
Ijerph 20 02251 g003
Table 1. Description of heart rate variability (HRV) indexes evaluated. Adapted from Shaffer F et al. [38].
Table 1. Description of heart rate variability (HRV) indexes evaluated. Adapted from Shaffer F et al. [38].
Sympathetic Nervous System ActivityUnit
Mean HRbpmAverage heart rate
Stress Index~Baevsky’s stress index, a geometric measure of HRV
LFnuRelative power between 0.04 Hz to 0.15 Hz
SNS Index~Calculated based on the mean HR, Baevsky’s stress index, and SD2 in normalized units
Peripheral Nervous System Activity
Mean RRmsAverage time of R-R intervals
RMSSDmsSquare root of mean squared difference between adjacent RR intervals
HFnuRelative power between 0.15 Hz to 0.4 Hz
pNN50%Percentage of successive RR intervals that differ by more than 50 ms
SD1msPoincaré plot standard deviation perpendicular the line of identity
PNS index~Calculated based on the mean RR, RMSSD, and SD1 in normalized units
Global Variability
SDNNmsStandard deviation of all normal RR intervals
RR Triangular Index~Integral of the density of the RR interval histogram divided by its height
TINNmsWidth of RR interval histogram
LF/HF%Ratio of low- and high-frequency power
SD2msPoincaré plot standard deviation along the line of identity
HR: heart rate; LF: low-frequency component; SNS index: sympathetic nervous system index; RMSSD: root mean square of differences between adjacent normal RR intervals in a time interval; HF: high-frequency component; pNN50: percentage of adjacent RR intervals with a difference in duration >50 ms; SD1: standard deviation of instantaneous beat-to-beat variability; SD2: standard deviation of long-term intervals between consecutive heartbeats; PNS index: parasympathetic nervous system index; SDNN: standard deviation of all normal RR intervals recorded in a time interval; TINN: triangular interpolation of NN interval; ms: millisecond; nu: normalized unit; HRV: heart rate variability; Hz: Hertz.
Table 2. Sample characterization and physical activity data.
Table 2. Sample characterization and physical activity data.
ControlPost-COVID
Baseline (n = 18)Follow-Up (n = 10)Intragroup AnalysisBaseline (n = 20)Follow-Up (n = 7)Intragroup AnalysisIntergroup Analysis
MedianQ1-Q3MedianQ1–Q3Difference between MeansMedianQ1-Q3MedianQ1–Q3Difference Between Means95% CIp-Value
Sex (M/F)13/5~8/2~ 9/11~4/3~~~0.1119
Age (years)26.6621.11–31.41~~~28.4724.73–33.77~~~~~
Weight (kg)71.6557.85–88.8562.7051.20–103.60−1.67 ± 8.5677.3565.75–90.6867.4059.70–72.25-10.31 ± 6.73-2.20 to 6.250.4116
Height (m)1.761.65–1.791.691.63–1.79−0.03 ± 0.041.711.61–1.771.701.64–1.760.004 ± 0.04−0.02 to 0.010.7348
BMI (kg/m2)23.9221.11–28.2023.6020.00–31.98−0.01 ± 2.2925.4523.13–31.3423.5319.76–26.41−3.58 ± 2.47−0.77 to 2.270.2982
Days between Positive Test and Follow-Up~~~~~37.5024.75–70.50~~~~~
Physical Activity
n = 16n = 8 n = 16n = 4
MedianIQRMedianIQR MedianIQRMedianIQR
Sedentary Activity (min/day)578.43534.88–676.68590.52467.67–664.25−92.97 ± 98.48499.59467.05–593.48601.15493.50–646.97−33.88 ± 146.52−212.98 to 373.820.5550
Light Activity (min/day)251.50188.54–291.00254.14220.00–407.3452.77 ± 33.62285.15206.75–432.30258.48235.24–306.71−42.21 ± 60.08−131.75 to 182.600.7260
Moderate Activity (min/day)17.905.14–29.8711.006.04–15.57−10.65 ± 8.5813.369.25–33.7227.8623.78–32.254.73 ± 11.36−28.06 to 25.000.9004
Vigorous Activity (min/day)3.650.00–15.820.000.00–2.36−7.73 ± 4.330.000.00–0.500.130.00–4.670.92 ± 1.00−10.63 to 3.260.2061
MVPA (min/day)23.7212.60–46.0011.006.36–17.13−19.13 ± 11.8614.799.57–34.0431.0624.92–32.725.51 ± 11.69−32.43 to 19.060.5756
Steps Count (steps/day)6570.654994.12–7795.244822.723714.61–8585.47−1296.60 ± 1706.155231.863696.86–9601.296350.225996.22–7197.71187.06 ± 1585.93−5935.78 to 4693.480.7998
Total Time in BSB 30–60 min Bouts169.4497.32–222.25161.0763.24–239.47−18.99 ± 39.9591.3673.36–134.18137.3698.78–187.9730.67 ± 35.92−124.54 to 137.290.9157
Total Time in BSB ≥ 60 min Bouts34.3119.17–99.6146.794.04–69.32−17.29 ± 22.4617.638.86–29.3230.072.25–57.36−4.35 ± 25.22−66.38 to 77.400.8678
Q: quartile range; F: female; M: male; kg: kilogram; m2: square meter; BMI: body mass index; MVPA: moderate to vigorous physical activities; BSB: breaks. in sedentary behavior.
Table 3. Comparisons of change in HRV indexes in the post-COVID and control groups and between groups.
Table 3. Comparisons of change in HRV indexes in the post-COVID and control groups and between groups.
Control (n = 11)Post-COVID (n = 7)
Intragroup AnalysisIntragroup AnalysisIntergroup Analysis
Dif Between MeansDif Between MeansUnadjusted 95% CIUnadjusted p-ValueAdjusted 95% CIAdjusted p-Value
SNS Activity
Mean HR4.67 ± 3.37−7.13 ± 5.59−25.83 to −4.330.0024 **−24.67 to −3.260.014 +
Stress Index−0.10 ± 0.98−3.26 ± 2.26−0.60 to 6.510.0965−6.62 to 0.570.093
LF (nu)3.35 ± 5.763.66 ± 6.17−31.16 to 2.170.0836−2.26 to 30.390.086
SNS Index0.31 ± 0.31−1.01 ± 0.66−2.48 to −0.290.0068 **−2.50 to −0.320.015 +
PNS Activity
Mean RR−44.54 ± 32.3860.36 ± 55.35 *33.50 to 222.500.0097 **33.72 to 225.510.012 +
RMSSD−3.05 ± 3.158.95 ± 4.08−17.06 to 0.800.1088−0.79 to 17.350.071
HF (nu)−3.35 ± 5.76−3.62 ± 6.17−2.33 to 31.110.0864−30.33 to 2.400.089
pNN50−2.84 ± 2.937.21 ± 2.61 *−14.62 to 0.190.0553−0.21 to 14.880.056
SD1−2.15 ± 2.246.35 ± 2.89−12.12 to 0.560.1038−0.55 to 12.320.07
PNS index−0.32 ± 0.200.54 ± 0.350.20 to 1.450.0091 **−0.20 to 1.470.013 +
Global Variability
SDNN0.91 ± 3.6311.54 ± 4.69−20.92 to 0.960.0709−1.31 to 21.400.078
RR Triangular Index1.02 ± 0.801.68 ± 1.01 *−2.97 to 1.530.506−1.61 to 3.060.517
TINN10.67 ± 14.2547.86 ± 20.32−79.60 to 12.720.1434−14.20 to 81.680.153
LF/HF0.23 ± 0.58−0.50 ± 1.03−3.00 to 0.750.2199−0.70 to 2.840.218
SD20.97 ± 4.5515.03 ± 6.17−27.16 to 1.650.0786−2.16 to 27.800.088
Dif Between Means: difference between means (follow-up—baseline); CI: confidence interval; HR: heart rate; LF: low-frequency component; nu: normalized unit; SNS index: sympathetic nervous system index; RMSSD: root mean square of differences between adjacent normal RR intervals in a time interval; HF: high-frequency component; pNN50: percentage of adjacent RR intervals with a difference in duration > 50 ms; SD1: standard deviation of instantaneous beat-to-beat variability; ms: millisecond; PNS index: parasympathetic nervous system index; SDNN: standard deviation of all normal RR intervals recorded in a time interval; RR triangular index: Integral of the density of the RR interval histogram divided by its height; TINN: triangular interpolation of NN interval; LF/HF: ratio of low and high frequency power; SD2: standard deviation of long-term intervals between consecutive heartbeats. *: Statistical significance (p < 0.05) between baseline and follow-up moments. **: Statistical significance (p < 0.05) in the difference between means comparing post-COVID-19 and control groups. +: Statistical significance (p < 0.05) of difference between moments between the post-COVID-19 and control groups after ANCOVA according to moderate to vigorous physical activity.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Freire, A.P.C.F.; Amin, S.; Lira, F.S.; Morano, A.E.v.A.; Pereira, T.; Coelho-E-Silva, M.-J.; Caseiro, A.; Christofaro, D.G.D.; Dos Santos, V.R.; Júnior, O.M.; et al. Autonomic Function Recovery and Physical Activity Levels in Post-COVID-19 Young Adults after Immunization: An Observational Follow-Up Case-Control Study. Int. J. Environ. Res. Public Health 2023, 20, 2251. https://doi.org/10.3390/ijerph20032251

AMA Style

Freire APCF, Amin S, Lira FS, Morano AEvA, Pereira T, Coelho-E-Silva M-J, Caseiro A, Christofaro DGD, Dos Santos VR, Júnior OM, et al. Autonomic Function Recovery and Physical Activity Levels in Post-COVID-19 Young Adults after Immunization: An Observational Follow-Up Case-Control Study. International Journal of Environmental Research and Public Health. 2023; 20(3):2251. https://doi.org/10.3390/ijerph20032251

Chicago/Turabian Style

Freire, Ana Paula Coelho Figueira, Shaan Amin, Fabio Santos Lira, Ana Elisa von Ah Morano, Telmo Pereira, Manuel-João Coelho-E-Silva, Armando Caseiro, Diego Giulliano Destro Christofaro, Vanessa Ribeiro Dos Santos, Osmar Marchioto Júnior, and et al. 2023. "Autonomic Function Recovery and Physical Activity Levels in Post-COVID-19 Young Adults after Immunization: An Observational Follow-Up Case-Control Study" International Journal of Environmental Research and Public Health 20, no. 3: 2251. https://doi.org/10.3390/ijerph20032251

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop