Front. Physiol. Frontiers in Physiology Front. Physiol. 1664-042X Frontiers Media S.A. 10.3389/fphys.2021.769085 Physiology Brief Research Report Cardiac Autonomic Modulation in Response to Muscle Fatigue and Sex Differences During Consecutive Competition Periods in Young Swimmers: A Longitudinal Study Castillo-Aguilar Matías 1 2 Valdés-Badilla Pablo 3 4 Herrera-Valenzuela Tomás 5 6 Guzmán-Muñoz Eduardo 7 Delgado-Floody Pedro 8 Andrade David Cristóbal 9 Moraes Michele M. 10 11 Arantes Rosa M. E. 10 11 Núñez-Espinosa Cristian 2 12 13 * 1Kinesiology Department, University of Magallanes, Punta Arenas, Chile 2Austral Integrative Neurophysiology Group, CADI-UMAG, Punta Arenas, Chile 3Department of Physical Activity Sciences, Faculty of Education Sciences, Universidad Católica del Maule, Talca, Chile 4Carrera de Entrenador Deportivo, Escuela de Educación, Universidad Viña del Mar, Viña del Mar, Chile 5Department of Sports Sciences and Physical Activity, Faculty of Health, Universidad Santo Tomás (UST), Santiago, Chile 6Department of Physical Activity, Sports and Health Sciences, Faculty of Medical Sciences, Universidad de Santiago de Chile (USACH), Santiago, Chile 7Escuela de Kinesiología, Facultad de Salud, Universidad Santo Tomás, Santiago, Chile 8Department of Physical Education, Sports and Recreation, Universidad de La Frontera, Temuco, Chile 9Centro de Investigación en Fisiología y Medicina de Altura (MedAlt), Facultad de Ciencias de la Salud, Universidad de Antofagasta, Antofagasta, Chile 10Department of Pathology, Institute of Biological Sciences, Universidade Federal de Minas Gerais, Belo Horizonte, Brazil 11Faculty of Medicine, Center for Newborn Screening and Genetics Diagnosis, Universidade Federal de Minas Gerais (NUPAD-FM/UFMG), Belo Horizonte, Brazil 12School of Medicine, University of Magallanes, Punta Arenas, Chile 13Interuniversity Center for Healthy Aging, Chile

Edited by: Lee Stoner, University of North Carolina at Chapel Hill, United States

Reviewed by: J. Derek Kingsley, Kent State University, United States; José Ramón Alvero Cruz, University of Malaga, Spain

 *Correspondence: Cristian Núñez-Espinosa, cristian.nunez@umag.cl

This article was submitted to Exercise Physiology, a section of the journal Frontiers in Physiology

 18 11 2021 2021 12 769085 01 09 2021 19 10 2021 Copyright © 2021 Castillo-Aguilar, Valdés-Badilla, Herrera-Valenzuela, Guzmán-Muñoz, Delgado-Floody, Andrade, Moraes, Arantes and Núñez-Espinosa. 2021 Castillo-Aguilar, Valdés-Badilla, Herrera-Valenzuela, Guzmán-Muñoz, Delgado-Floody, Andrade, Moraes, Arantes and Núñez-Espinosa

This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

Objective: To study the differences in cardiac autonomic modulation in response to muscle fatigue caused by high-intensity exercise during two consecutive competition periods in young swimmers.

Methods: Twenty-six competitive swimmers, selected by their training volume, were separated in two groups, females (n = 12 [46%], age: 13.5 ± 1.4 years) and males (n = 14 [54%], age: 13.9 ± 1.7 years), aged between 10 and 16 years, were evaluated five times as follow: (i) 21 days before the first competition (t-0); (ii) two days before (t-1; t-3); and (iii) two days after (t-2; t-4) of the first and second competitions. Morphological measurements (body mass, percentage of total body fat and height), blood pressure, power, and resting heart rate variability (RR with Polar band) were recorded before and after Wingate test at each time.

Results: Body fat was higher in females compared to males. However, no differences were found in other morphological parameters. An intra-subject analysis grouped by sex in cardiovascular parameters shows longitudinal variations in systolic pressure and mean pressure among females. Additionally, females depicted higher, very low frequency (VLF, which is intrinsically generated by the heart and strongly associated with emotional stress) after physical fatigue compared to males at t-1. Further, before the competition, the high frequency (HF) component of HRV (parasympathetic drive) was higher in males than females at t-0 and t-4.

Conclusion: Our data revealed that males displayed greater parasympathetic reactivity after an anaerobic muscle fatigue test during their competition periods. Contrarily, females had a less cardiac autonomic modulation when comparing the pre-post Wingate test after two consecutive competition periods.

 arterial pressure autonomic nervous system heart rate variability physical exertion swimming Introduction

Young swimmer athletes display different autonomic nervous system (ANS) modulation of the heart compared to non-trained control subjects (Arce-Álvarez et al., 2021). Considering the amount of research in this area, cardiac autonomic modulation by heart rate variability (HRV) has emerged as a highly effective tool to access neurophysiological adaptations demanded at different athlete’s training moments (Clemente-Suárez et al., 2015; Arce-Álvarez et al., 2021). Of note, it has been relevant to the point that it is capable of overtraining states (Kamandulis et al., 2020) and estimating athletic performance (Chalencon et al., 2015). Notably, there are several techniques and calculations focused on determining autonomic modulation; nevertheless, most are invasive maneuvers, which reduce their practical use (i.e., sports activities). Thus, it has been proposed that HRV calculated from signals obtained from portable devices could reflex autonomic modulation of the heart with high reliability (Georgiou et al., 2018).

Interestingly, it is very well known that during training and competition, muscle fatigue is generated, and it has been proposed that this fatigue could play a pivotal role in the cardiac autonomic response of an athlete (Buchheit et al., 2011; Mayo et al., 2016). Indeed, low-intensity training for long periods and short high-intensity exercises with long recovery periods seems to be a stronger modulator of cardiac autonomic response swimmers’ athletes during training (Clemente-Suárez and Arroyo-Toledo, 2018; Pla et al., 2019). In this sense, recovery of parasympathetic discharge is modulated by training intensity, with delay for recovery ANS/HRV in athletes exercised at and above the first ventilatory threshold (Seiler et al., 2007).

Muscle fatigue favors parasympathetic withdrawal, which could affect the recovery period; however, it is expected that this phenomenon will attenuate with adequate rest and periodization, not affecting the athlete’s recovery. Thus, adequate rest and periodization promote the improvement of the neurophysiological during the tapering period, allowing the swimmer a better adaptation when competing (Flatt et al., 2017). Of note, during the competition period and/or high physical demand, swimmers’ physical preparation times can change, which could play an important role in the cardiac autonomic response, affecting the recovery periods and consequently physiological adaptations promoted by the training (Edmonds et al., 2015). In addition, the presence of differences related to age and sex in heart rate variability despite maintaining an adequate level of physical activity (based on accelerometers) can be observed at these ages (Spina et al., 2019); however, there is no evidence according to differences between females and males about the cardiac autonomic response to muscle fatigue during the competitive period.

Sex determines differences in the biological organization at all levels (Holdcroft, 2007), including cardiac modulatory balance, which has been related to a greater tone in the autonomic vagal branch in females and a greater tone in the sympathetic branch in males (Spina et al., 2019). As effort intensity also influences HRV recovery, it is noteworthy that females have an advantage in resistance to fatigue during submaximal contractions, with less impairment in neuromuscular activation after strenuous exercise; however, this advantage diminishes as the intensity of contractions increases (Azevedo et al., 2021). Therefore, considering the differences in cardiac autonomic modulation and resistance to fatigue between males and females, muscle fatigue during a competitive period may influence cardiac autonomic control in a sex-specific manner during a competitive period in swimming athletes.

At present, there is no conclusive evidence about how cardiac autonomic response could be affected during consecutive competition periods, which could affect more young swimmers who are faced with physiological and anthropometric changes typical of age (Tanina et al., 2017; Weinstein et al., 2018). Also, as differences in athletic performance of males and females begin at the age of 12–13 years (Handelsman, 2017), cardiovascular autonomic adaptation over a training season may diverge between the sexes of young swimmers. Then, we aimed to evaluate the differences in cardiac autonomic modulation through non-invasive measures in response to muscle fatigue caused by high-intensity exercise during two consecutive competition periods in young swimmers. We hypothesized that the muscular fatigue by high-intensity physical exercise during consecutive competition periods generates a cardiac autonomic response that differs according to the sex of the swimmer.

 Materials and Methods Design

A descriptive, comparative, longitudinal, and quantitative approach was used. The participants were selected through intentional non-probabilistic sampling, distributed between females (n = 12) and males (n = 14). Athletes were monitored during two consecutive dates of the Patagonian International Circuit. The swimmers were observed five times. The first evaluation was carried out 21 days before the first competition (t-0). The four remaining assessments were carried out two days before (t-1; t-3; “pre-competition”) and two days after (t-2; t-4; “post-competition”) the first and second competitions. In each experimental time, anthropometric parameters (body weight, height, and percentage of body fat), cardiovascular (blood pressure), as well as HRV were measured: before and after the anaerobic muscle fatigue test (Wingate test). The study design is represented in Figure 1.

Study design and evaluations were carried out. Heart rate variability (HRV) was evaluated before and after the anaerobic physical fatigue test (Wingate test).

 Participants

Twenty-six competitive young swimmers of national and international competition level were recruited to participate in this study (age: 13.7 ± 1.5 years; height: 1.62 ± 0.07 m; body mass: 56.5 ± 10.5 kg; body fat: 24.4 ± 3.2%) and separate on females (age: 13.5 ± 1.4 years; height: 1.58 ± 0.03 m; body mass: 56.3 ± 9.1 kg; body fat: 30.3 ± 3.5%) and males (age: 13.9 ± 1.7 years; height: 1.65 ± 0.08 m; body mass: 56.8 ± 12.1 kg; body fat: 18.4 ± 3.6%).

The swimmers were recruited from the Fiscal Gymnasium of the Magallanes and Chilean Antarctic region. The inclusion criterion was to have a minimum of 3 years in competitive swimming, train at least six times per week, and have at least 14 h per week of training. Exclusion criteria were: take any supplements or medications that could affect the heart rate, suffer musculoskeletal injuries in the last three months, presence of pain at the time of performing the evaluations. No participants met the exclusion criteria. The volunteers and their legal guardians were informed about the objectives, procedures, responsibilities, and risks of participating in the study. Participating subjects gave their permission, and legal guardians provide informed consent before participation. The Ethics Committee approved this study of the University of Magallanes, Chile (code: N°141CEC2018), following the regulations established by the Declaration of Helsinki on ethical principles in human beings.

 Procedure Training Monitoring

The swimmers were monitored during two consecutive dates of the Patagonian International Circuit.

All the swimmers were trained to perform short pool speed tests (25 m). During the training periods, before the first (between t-0 and t-1) and second (between t-2 and t-3) competition, the training routine was maintained. Throughout each pre-competitive period, a progressive decrease in training volume was carried out, starting at 7,000 m per week and ending with 3,500 m per week before the competition, progressively increasing the intensity (heart rate controlled) and giving the swimmer a longer recovery time.

To better optimize the time required for the evaluations in the athletes, measurement stations were established within the same laboratory: station 1, the athlete arrives at the laboratory, rests in a chair for 5 min, and then their blood pressure is controlled; station 2, the athlete is evaluated on his morphological measurements (approximately 10 min); Station 3, the subject remained lied on the stretcher in a supine position, and his HRV was evaluated prior to the anaerobic muscle fatigue test; station 4, the athlete is evaluated in the Wingate test, to later return to station 3, for your second HRV assessment, after the fatigue test.

 Acute Muscle Fatigue Protocol

Among the conditions of dress, the participants were asked to wear a shirt, shorts, and footwear. Before the assessments, all participants were instructed to (a) get adequate rest the night before, sleep 8 h or more, (b) not consume stimulant drinks or medications before measurements, (c) consume ∼ 2 liters of water during the day before, (d) eat regularly without altering your diet. The participants reported to the laboratory 15 min before the test. The Wingate protocol was carried out at 22°C and 15% RH controlled by the air conditioning in a laboratory equipped for the study. For HRV, the volunteers lay down on a stretcher in the supine position, and the recording of HRV was started. During the day of testing, participants were instructed to give their maximum effort during testing.

 Assessments Morphological Measures

Body mass (kg) and the total body fat (%) were assessed by bioimpedance using the Tanita BC-558 Ironman Segmental Body Composition Monitor (Tanita Ironman, Arlington Heights, IL, United States), with a concordance 89.3% compared to the Dual X-ray Absorption test using standard measurement protocols (Mialich et al., 2011; Calderón and Rodriguez-Hernandez, 2018). Height was measured by CHARDER® HM230M manual height rod (Charder Electronics Co., Ltd. No.103, Guozhong Rd., Taiwan, R.O.C.). The body mass index (BMI) and the fat-free mass index (FFMI) were calculated as follow, BMI: [body weight]/[height]2 (kg/m2) and FFMI: [fat-free mass]/[height]2 (kg/m2).

 Anaerobic Muscle Fatigue Test

To assess the anaerobic muscular endurance, the Wingate anaerobic test was performed. This test is used to measure the anaerobic capacity and power of an individual (Vandewalle et al., 1987) and has been studied extensively in children and young people (King-Dowling et al., 2018), proving to be a safe tool with good rates of reproducibility (Bar-Or, 1987). A cycle ergometer test was carried out with an individualized load for each athlete, as described before (Bar-Or, 1993). The application of the test allowed us to calculate the minimum, mean, and peak power output as follow: Load (kp) × revolutions in 5 s × 11.765; where the minimum, average and maximum number of revolutions were used for each power measurement (Bar-Or, 2012). Each athlete was constantly consulted for discomfort or pain during the performance of the test through verbal communication.

 Cardiovascular Parameters

Blood pressure (Omron® Pressure Monitor), systolic blood pressure (SP), and diastolic blood pressure (DP) were measured. The evaluation was carried out with the subject sitting in a chair, allowing to calculating mean arterial pressure (MAP) and pulse pressure (PP). Cardiac autonomic modulation was determined via a recording of RR intervals. The volunteers remained lied on the stretcher in a supine position during the entire HRV measurement procedure, and the RR intervals were continuously recorded during the last 10 min of rest, and 5 min were analyzed. The breathing rate of the subjects were spontaneously. Artifacts and ectopic heartbeats (which did not exceed 3% of the recorded data) were excluded (Task Force of the European Society of Cardiology the North American Society of Pacing Electrophysiology, 1996). The time-domain parameters considered for the analysis were the square root of the mean squared differences of the successive RR intervals (RMSSD, expressed in ms), which reflect the parasympathetic influence (Buchheit et al., 2010) and the standard deviation of the RR intervals (SDNN), which is believed to reflect the total variability, that is, the sympathetic and parasympathetic contribution of the autonomic nervous system on the heart (Berntson et al., 1997; Buchheit and Gindre, 2006). The frequency domains considered were the high frequency (HF) power band that reflects the parasympathetic influence and respiratory sinus arrhythmia (Akselrod et al., 1981) and the low frequency (LF) band associated with baroreflex activity (Goldstein et al., 2011). Very low-frequency band (VLF), which is multi-faceted and is strongly associated with emotional stress (Malik et al., 1996; Fisher et al., 2014; McCraty and Shaffer, 2015). Finally, all the data obtained were digitized and analyzed using the software Kubios HRV® (Tarvainen et al., 2014).

 Statistical Analysis

The report of continuous variables is expressed as mean ± standard deviation (SD) or median and interquartile range (IQR) based on the proximity to the normality of the underlying distribution; categorical variables are presented as absolute (n) and relative frequency (%). To verify the assumptions of parametric statistics, we used the Shapiro-Wilk test to assess the normality of the distribution of the variables analytically and the Levene’s test for homogeneity of variance across groups.

For group comparisons, Student t-test (tStudent) or Wilcoxon rank-sum test (W) according to the fulfillment of each of the parametric assumptions. For repeated measures, the Friedman rank sum test (χFriedman2) was used.

For Student’s t-test, Cohen’s d (dCohen) was used as a measure of effect size (ES) following the conventions of Cohen (2013), with ES < 0.2 a very small effect, 0.2 ≤ ES < 0.5 small, 0.5 ≤ ES < 0.8 moderate, and ES ≥ 0.8 a large effect. For the Wilcoxon signed-rank and rank-sum test, the rank-biserial correlation (r^biserial) was computed and interpreted according to Funder and Ozer (2019), meaning ES < 0.1 very small, 0.1 ≤ ES < 0.2 small, 0.2 ≤ ES < 0.3 medium, 0.3 ≤ ES < 0.4 large, and ES ≥ 0.4 very large. For each ES, the 95% Confidence Interval (CI95%) was also calculated.

A probability of committing a type I (α) error of less than 5%, i.e., a p < 0.05, was considered sufficient evidence for statistical significance in hypothesis testing. The statistical analysis was performed using the statistical programming language R (version 4.1.0) (R Core Team, 2021) and the statistical packages necessary for this study (Wickham, 2016; Ho et al., 2019; Ben-Shachar et al., 2020; Xie et al., 2020; Castillo Aguilar, 2021; Kassambara, 2021; RStudio Team, 2021).

 Results Participants Characterization

At t-0, we found similar morphological characteristics between males and females, particularly in body weight (tStudent (22) = −0.11, p = 0.91, dCohen = −0.05, CI95%[−0.85, 0.75]), FFMI (tStudent (22) = −1.73, p = 0.097, dCohen = −0.71, CI95%[−1.53, 0.13]), and BMI (tStudent (22) = 1.23, p = 0.231, dCohen = 0.50, CI95%[−0.32, 1.31]). As expected, the groups were not homogeneous in all morphological measures at baseline, as males were taller on average (1.65 ± 0.08 m) than females (1.58 ± 0.03 m; tStudent (24) = −2.95, p = 0.007, dCohen = −1.16, CI95%[−1.99, −0.31]). Similarly, females had a higher body fat (30.2 ± 3.5%) than their male counterparts (18.3 ± 3.4%; tStudent (22) = 8.44, p < 0.001). However, when comparing BMI or FFMI, we didn’t find any statistically significant difference between sexes (BMI, tStudent (22) = 1.23, p = 0.231; FFMI, tStudent (22) = −1.73, p = 0.097). Estimates of the magnitude of sex differences in BMI, body fat, and FFMI, using the mean value for each period for each participant, are shown in Figure 2.

Differences in body composition parameters compared by sex and respective effect sizes. (A) Body mass index (BMI); (B) fat mass; (C) fat-free mass index (FFMI). The dispersion of the sample by sex and the size of the effect through Cohen’s d (ES) are shown. 95% confidence intervals for ES were calculated using bootstrap resampling with 5,000 replications. They are also bias-corrected and accelerated (BCa).

In cardiac parameters, a within-participants analysis by sex shows variations in systolic pressure among females (χFriedman2 (4) = 18.60, p < 0.001, WKendall = 0.37, CI95%[0.26, 0.63]), whereas in males these values remained constant over time periods (χFriedman2 (4) = 2.06, p = 0.725, WKendall = 0.05, CI95%[0.02, 0.46]). The same pattern was seen in mean arterial pressure (females, χFriedman2 (4) = 14.92, p = 0.005, WKendall = 0.29, CI95%[0.13, 0.58]; males, χFriedman2 (4) = 8.94, p = 0.063, WKendall = 0.23, CI95%[0.10, 0.67]). Regarding the blood pressure parameters, although the PP parameters were consistently higher in males (50.5 (23.5) mmHg vs females, 37.5 (11) mmHg), there was only a significant difference in the parameters at t-0 (W = 34, p = 0.03, r^biserial = −0.53, CI95%[−0.79, −0.12]) when comparing both groups.

A similar situation was observed when analyzing the anaerobic muscle fatigue and power parameters, where males registered higher values than females, but only at t-0 (males, 402.8 ± 125 W; females, 304.7 ± 108.3 W), tStudent (23) = −2.09, p = 0.048, dCohen = −0.84, CI95%[−1.65, −0.01]. All values for blood pressure and anaerobic muscle fatigue test results are presented in Table 1.

Blood pressure and anaerobic muscle fatigue test results are presented for each of the time periods separated by sex.

Female
 Male
t-0 t-1 t-2 t-3 t-4 t-0 t-1 t-2 t-3 t-4
n = 12 n = 12 n = 12 n = 12 n = 12 n = 12 n = 14 n = 12 n = 12 n = 11
SPnp (mmHg) 104.5−(7.2) 115−(11.0)La 111−(7.0)La 111−(7.0)La 107.5−(7.2) 118−(14.5) 114−12.8 115−(13.8) 111−(10.5) 114−(16)
DPnp (mmHg) 67−(5.8) 72−6.2 69−(4.0)La 69−(4.0)La 68−(9) 62−(8.2) 69−17.5 66−(9.8) 66−(13.5) 63−(9.5)
MAPnp (mmHg) 78.3−(8.4) 85.5−(3.4)La 82.5−(8.2)La 82.5−(7.4)La 80.8−(7) 79.5−(9.5) 85.3−10.5 84.8−(11) 82.5−(10.3) 81.7−(8.8)
PPnp (mmHg) 37.5−(11) 46.5−(10.0)La 43−(7.2)La 43−(7.2)La 42.5−(6.8) 50.5−(23.5) 44.5−16.2 51.5−(10) 45.5−(17.5) 46−(10.5)
PPOp (W) 304.7 ± 108.3 243.9 ± 90.0La 269.5 ± 90.6Mb 274.8 ± 98.3La 297.1 ± 85.2Mcd 402.8 ± 125 310.7 ± 134.9La 341.3 ± 151.9La–Mb 361.5 ± 130.6Lb 372 ± 135
MePOp (W) 284.2 ± 88.8 224.4 ± 79.2La 239 ± 75.6 240.9 ± 90.3La 266.4 ± 70.5Ld 368.1 ± 119.2 295.9 ± 125.5La 303 ± 138.3La 325.8 ± 118.4Lab 336.6 ± 105.5Mc
MiPOp (W) 236.8 ± 71.7 181.3 ± 64.1La 181.8 ± 57.6Ma 206.8 ± 73.6 228 ± 68.2Lc 301.8 ± 90.4 238.6 ± 102.2La 223.5 ± 94.0La 268.9 ± 100.4La 273.8 ± 88.3Mc
Fatiguep (%) 20.8 ± 10 24.3 ± 9.5 30.4 ± 11.8Ma 24.2 ± 9 23 ± 13.6 27.1 ± 11 23.9 ± 10.8 33.7 ± 20.2 23.2 ± 13.3 25.6 ± 11

Values are presented as median−(IQR) or mean ± SD as appropriate. All pairwise comparison was performed using Wilcoxon signed-rank test (np) or Student t-test (p) and for effect-size (rank-biserial correlation or Cohens’d) codes were: S, small; M, moderate; L, large. For blood pressure code: SP, systolic pressure; DP, diastolic pressure; MAP, mean arterial pressure; PP, pulse pressure. For Wingate test parameters: PPO, peak power output; MePO, mean power output; MiPO, minimum power output. The differences within the same group indicate a p < 0.05, following the code: a, compared with t-0; b, compared with t-1; c, compared with t-2; d, compared with t-3; e, compared with t-4.

When analyzing HRV parameters, we observed longitudinal variations only in the VLF-post anaerobic fatigue test frequency domain (χFriedman2 (4) = 13.16, p = 0.011, WKendall = 0.14, CI95%[0.07, 0.31]), while the other autonomic-cardiac parameters remained constant when not considering sex as an influencing factor. However, after separating our analyses by sex, we found that this relationship with VLF-post is valid only for males (χFriedman2 (4) = 11.64, p = 0.02, WKendall = 0.22, CI95%[0.16, 0.46]), but not for females (χFriedman2 (4) = 9.24, p = 0.055, WKendall = 0.21, CI95%[0.11, 0.53]). When comparing directly between sexes, we observed that females had a higher VLF-post anaerobic fatigue test than males at t-1, while at HF-pre anaerobic fatigue test was higher in males than females at t-0. In the time domain, females reached lower levels of SDNN, both pre and post-anaerobic fatigue test at t-4, but not in the other periods. Similarly, we found large differences between males and females in RMSSD pre- but not post-anaerobic fatigue test (see Figure 3).

Heart rate variability parameters were compared by sex within each time. (A) Frequency domain; (B) time domain. The median of the values is presented in a red line for the females and blue for the males, while the red and blue areas represent the corresponding interquartile range. § Significantly different (p < 0.05) and rank-biserial correlation code were: M, moderate; L, large.

 Discussion

Our study aimed to evaluate the differences in the cardiac autonomic modulation of young competitive swimmers through non-invasive measures and in response to muscle fatigue caused by high-intensity exercise during two consecutive competition periods. Our findings in the cardiac parameters show that the male had a significant parasympathetic predominance (denoting major cardiac autonomic modulation), evidenced in the difference in HF-pre compared to females in t-0 and t-4. Although the HF-pre values progressively decreased in males, it remains higher than the values registered in females. We believe that these differences could be due to two factors. The first is caused by non-modifiable variables (age and sex), which can be strong predictors of post-exercise cardiac recovery in young (Jezdimirovic et al., 2017); The second factor can be given by body fat mass. It has been observed in prepubertal children in Europe that a higher amount of body fat (31.4% in males and 37.3% in females) is related to a decrease in parasympathetic modulation (Santos-Magalhaes et al., 2015), which coincidentally resembles the values evaluated in competitive swimmers in our study (30.3 ± 3.5% fat mass). Similarly, elevated body fat levels have been linked to altered states in cardiac autonomic activity after a 20-min jog event (Santana et al., 2019), considering that the delay in autonomic recovery after exercise has been associated with an increased risk of cardiovascular events (Albert et al., 2000; Rossi et al., 2015). These would indicate that many genetic, sociocultural, lifestyle, and morphological variables may directly impact the baseline response of cardiac autonomic modulation (Delgado-Floody et al., 2020).

Concerning the pre and post-anaerobic fatigue test, at t-1, both sexes followed a similar pattern in terms of anaerobic power records compared to the muscle fatigue test. In t-1, the different powers decrease compared to t-0, to then increase progressively to t-4, which would be explained by the training and progressive tapering of sports planning (Le Meur et al., 2012; Vachon et al., 2021). Before the anaerobic test, the difference in the autonomic parameters between the groups decreased progressively until t-3, where the frequency domain values are practically unified. This decrease in cardiac autonomic control is observed before the second competition, neurophysiologically showing the autonomic wear of the swimmers in the face of an anaerobic fatigue test (Schneider et al., 2019). However, in t-4, the groups again differentiate, showing greater cardiac autonomic and especially parasympathetic reactivity on the part of males in both RMSSD-post and SDNN-post. After the anaerobic muscle fatigue test, the cardiac autonomic response shows that females at t-2 show greater changes at VLF-post, which tend to remain in time until t-3. It is known that the cardiac autonomic recovery of swimmers differs in terms of the type of test they perform. In more extended tests, the recovery is faster, while in speed tests, the recovery is slower (Piras et al., 2019). This difference could be affecting the females in our study more, given her response to emotional stress in relation to the evaluations carried out in this study. In this sense, it is noteworthy that increased effort intensity supplants the female’s advantage in resistance to fatigue (Spina et al., 2019).

In the case of the males, they progressively show greater post-test cardiac autonomic modulation, which becomes different from females in SDNN-post during t-4. This difference may be due to a better cardiac adaptation to training, represented by a vagal predominance in cardiac autonomic modulation (Triposkiadis et al., 2002), which could be better in males than in females due to their different morphological characteristics (Bassareo and Crisafulli, 2020; Hottenrott et al., 2021).

Regarding morphological parameters, a similar character was evidenced in the participants. However, we observed that at all times, females registered higher levels of body fat than males. Although this is expected, other studies have linked a high percentage of body fat with low sympathetic activity (Rossi et al., 2014) and decreased parasympathetic modulation in time-domain HRV (Santos-Magalhaes et al., 2015). Therefore, these antecedents could help understand the differences in the neurophysiological response to the stress of exercise of young swimmers. However, these should be studied in depth in future studies in the area.

 Limitations

The main limitations of this study were the small sample size and the non-probabilistic selection of the participants, which could limit the external validity of the study. In addition, for future research, we suggest controlling the hours of rest and sleep of swimmers, which could be factors that influence the physical and physiological performance of athletes. We also believe that it would be appropriate to be able to use an anaerobic muscle fatigue test carried out in the aquatic environment, which could better assimilate the muscle fatigue that athletes perform under competition conditions. Despite this, our study reports novel information on how the autonomic response varies during consecutive competition periods of young swimmers, which has been little studied and could affect their neurophysiological response in competition.

 Conclusion

There are differences in the cardiac autonomic modulation of young swimmers during consecutive periods of competition. Males show greater parasympathetic reactivity after the anaerobic muscle fatigue test after two consecutive competition periods, while females show less cardiac autonomic modulation during the same competition periods.

 Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

 Ethics Statement

The studies involving human participants were reviewed and approved by The Ethics Committee of the University of Magallanes, Chile (N°141CEC2018). Written informed consent to participate in this study was provided by the participants’ legal guardian/next of kin.

 Author Contributions

All authors listed have made a substantial, direct and intellectual contribution to the work, and approved it for publication.

 Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

 Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

 Funding

This work was funded by MAG2095 (Chilean Education Ministry). CN-E received resources from the National Fund for the Promotion of Sports of Chile, code 1800120339 (Instituto Nacional de Deporte de Chile, IND). MM received post-doctoral fellowship from Coordenação de Aperfeiçoamento de Pessoal de Nível Superior – Brasil (CAPES/BRASIL) – Finance Code 001 [88887.321687/2019–00]. RMEA received research fellowship from Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq/BRASIL) [305952/2017–0] and Grant CNPq [442645/2018–0].

We thank all study participants, parents, and their coaches for their contribution, especially to Mr. José Lezama. For analysis, reporting, plots, and tables, the following R packages were used: dabestr (Ho et al., 2019), writR (Castillo Aguilar, 2021), effectsize (Ben-Shachar et al., 2020), rstatix (Kassambara, 2021), ggplot2 (Wickham, 2016), and rmarkdown (Xie et al., 2020), in the RStudio (RStudio Team, 2021) graphical interface.

 References Akselrod S. Gordon D. Ubel F. A. Shannon D. C. Berger A. Cohen R. J. (1981). Power spectrum analysis of heart rate fluctuation: a quantitative probe of beat-to-beat cardiovascular control. Science 213 220–222. 10.1126/science.6166045 6166045 Albert C. M. Mittleman M. A. Chae C. U. Lee I.-M. Hennekens C. H. Manson J. E. (2000). Triggering of sudden death from cardiac causes by vigorous exertion. N. Engl. J. Med. 343 1355–1361. 10.1056/nejm200011093431902 11070099 Arce-Álvarez A. Veliz C. Vazquez-Muñoz M. von Igel M. Alvares C. Ramirez-Campillo R. (2021). Hypoxic respiratory chemoreflex control in young trained swimmers. Front. Physiol. 12:215. 10.3389/fphys.2021.632603 33716781 Azevedo R. Forot J. Iannetta D. MacInnis M. Millet G. Murias J. (2021). Slight power output manipulations around the maximal lactate steady state have an impact on fatigue in females and males. J. Appl. Physiol. 130 1879–1892. 10.1152/japplphysiol.00892.2020 33914658 Bar-Or O. (1987). The Wingate anaerobic test an update on methodology, reliability and validity. Sports Med. 4 381–394. 10.2165/00007256-198704060-00001 3324256 Bar-Or O. (1993). Test anaeróbico wingate. Rev. Actual. Cienc. Deporte 1 1–7. 10.1007/978-3-319-56782-2_173-2 Bar-Or O. (2012). Pediatric Sports Medicine for the Practitioner: From Physiologic Principles to Clinical Applications. New York, NY: Springer Science & Business Media. Bassareo P. P. Crisafulli A. (2020). Gender differences in hemodynamic regulation and cardiovascular adaptations to dynamic exercise. Curr. Cardiol. Rev. 16 65–72. 10.2174/1573403X15666190321141856 30907327 Ben-Shachar M. S. Lüdecke D. Makowski D. (2020). effectsize: estimation of effect size indices and standardized parameters. J. Open Sour. Softw. 5:2815. 10.21105/joss.02815 Berntson G. G. Thomas Bigger J. Jr. Eckberg D. L. Grossman P. Kaufmann P. G. Malik M. (1997). Heart rate variability: origins, methods, and interpretive caveats. Psychophysiology 34 623–648. 10.1111/j.1469-8986.1997.tb02140.x 9401419 Buchheit M. Al Haddad H. Mendez-Villanueva A. Quod M. J. Bourdon P. C. (2011). Effect of maturation on hemodynamic and autonomic control recovery following maximal running exercise in highly trained young soccer players. Front. Physiol. 2:69. 10.3389/fphys.2011.00069 22013423 Buchheit M. Gindre C. (2006). Cardiac parasympathetic regulation: respective associations with cardiorespiratory fitness and training load. Am. J. Physiol. Heart Circ. Physiol. 291 H451–H458. 10.1152/ajpheart.00008.2006 16501030 Buchheit M. Chivot A. Parouty J. Mercier D. Al Haddad H. Laursen P. (2010). Monitoring endurance running performance using cardiac parasympathetic function. Eur. J. Appl. Physiol. 108 1153–1167. Calderón T. M. L. Rodriguez-Hernandez M. (2018). Efecto de un programa de 18 semanas de actividad fıìsica sobre la capacidad aeróbica, la fuerza y la composición corporal en personas adultas mayores. Pensamiento Actual. 18 125–135. Castillo Aguilar M. A. (2021). writR: Inferential Statistics and Reporting in APA Style, R Package Version 1.0.1. 10.5281/zenodo.4603838 Chalencon S. Pichot V. Roche F. Lacour J.-R. Garet M. Connes P. (2015). Modeling of performance and ANS activity for predicting future responses to training. Eur. J. Appl. Physiol. 115 589–596. 10.1007/s00421-014-3035-2 25359446 Clemente-Suárez V. J. Arroyo-Toledo J. J. (2018). The use of autonomic modulation device to control training performance after high-intensity interval training program. J. Med. Syst. 42:47. 10.1007/s10916-018-0907-8 29372348 Clemente-Suárez V. J. Fernandes R. J. Arroyo-Toledo J. Figueiredo P. González-Ravé J. M. Vilas-Boas J. (2015). Autonomic adaptation after traditional and reverse swimming training periodizations. Acta Physiol. Hung. 102 105–113. 10.1556/APhysiol.102.2015.1.11 25804392 Cohen J. (2013). Statistical Power Analysis for the Behavioral Sciences. New York, NY: Academic press. Delgado-Floody P. Alvarez C. Caamaño-Navarrete F. Jerez-Mayorga D. Latorre-Román P. (2020). Influence of mediterranean diet adherence, physical activity patterns, and weight status on cardiovascular response to cardiorespiratory fitness test in chilean school children. Nutrition 71:110621. 10.1016/j.nut.2019.110621 31874334 Edmonds R. Burkett B. Leicht A. McKean M. (2015). Effect of chronic training on heart rate variability, salivary IgA and salivary alpha-amylase in elite swimmers with a disability. PLoS One 10:e0127749. 10.1371/journal.pone.0127749 26043224 Fisher A. Groves D. Eleuteri A. Mesum P. Patterson D. Taggart P. (2014). Heart rate variability at limiting stationarity: evidence of neuro-cardiac control mechanisms operating at ultra-low frequencies. Physiol. Meas. 35:309. 10.1088/0967-3334/35/2/309 Flatt A. A. Hornikel B. Esco M. R. (2017). Heart rate variability and psychometric responses to overload and tapering in collegiate sprint-swimmers. J. Sci. Med. Sport. 20 606–610. 10.1016/j.jsams.2016.10.017 27890479 Funder D. C. Ozer D. J. (2019). Evaluating effect size in psychological research: sense and nonsense. Adv. Methods Pract. Psychol. Sci. 2 156–168. Georgiou K. Larentzakis A. V. Khamis N. N. Alsuhaibani G. I. Alaska Y. A. Giallafos E. J. (2018). Can wearable devices accurately measure heart rate variability? A systematic review. Folia Med. 60 7–20. Goldstein D. S. Bentho O. Park M.-Y. Sharabi Y. (2011). Low-frequency power of heart rate variability is not a measure of cardiac sympathetic tone but may be a measure of modulation of cardiac autonomic outflows by baroreflexes. Exp. Physiol. 96 1255–1261. 10.1113/expphysiol.2010.056259 21890520 Handelsman D. J. (2017). Sex differences in athletic performance emerge coinciding with the onset of male puberty. Clin. Endocrinol. 87 68–72. 10.1111/cen.13350 28397355 Ho J. Tumkaya T. Aryal S. Choi H. Claridge-Chang A. (2019). Moving beyond P values: everyday data analysis with estimation plots. Nat. Methods 16, 565–566. 10.1038/s41592-019-0470-3 31217592 Holdcroft A. (2007). Integrating the dimensions of sex and gender into basic life sciences research: methodologic and ethical issues. Gend. Med. 4 S64–S74. 10.1016/s1550-8579(07)80048-9 Hottenrott L. Ketelhut S. Schneider C. Wiewelhove T. Ferrauti A. (2021). Age-and sex-related differences in recovery from high-intensity and endurance exercise: a brief review. Int. J. Sports Physiol. Perform. 16 752–762. 10.1123/ijspp.2020-0604 33883293 Jezdimirovic T. Stajer V. Semeredi S. Calleja-Gonzalez J. Ostojic S. M. (2017). Does body fat percentage predict post-exercise heart rate response in non-obese children and adolescents? J. Pediatr. Endocrinol. Metab. 30 629–633. 10.1515/jpem-2016-0468 28525355 Kamandulis S. Juodsnukis A. Stanislovaitiene J. Zuoziene I. J. Bogdelis A. Mickevicius M. (2020). Daily resting heart rate variability in adolescent swimmers during 11 weeks of training. Int. J. Environ. Res. Public Health 17:2097. 10.3390/ijerph17062097 32235693 Kassambara A. (2021). rstatix: Pipe-Friendly Framework for Basic Statistical Tests, R Package Version 0.7.0. Available online at: https://CRAN.R-project.org/package=rstatix King-Dowling S. Rodriguez C. Missiuna C. Timmons B. Cairney J. (2018). Health-related Fitness in preschool children with and without motor delays. Med. Sci. Sports Exerc. 50 1442–1448. 10.1249/mss.0000000000001590 29474207 Le Meur Y. Hausswirth C. Mujika I. (2012). Tapering for competition: a review. Sci. Sports 27 77–87. 10.1016/j.scispo.2011.06.013 Malik M. Bigger J. Camm A. Kleiger R. Malliani A. Moss A. (1996). Heart rate variability: standards of measurement, physiological interpretation, and clinical use. Eur. Heart J. 17 354–381. 10.1093/oxfordjournals.eurheartj.a014868 Mayo X. Iglesias-Soler E. Fariñas-Rodrıìguez J. Fernández-del-Olmo M. Kingsley J. D. (2016). Exercise type affects cardiac vagal autonomic recovery after a resistance training session. J. Strength Cond. Res. 30 2565–2573. 10.1519/jsc.0000000000001347 26817741 McCraty R. Shaffer F. (2015). Heart rate variability: new perspectives on physiological mechanisms, assessment of self-regulatory capacity, and health risk. Glob. Adv. Health Med. 4 46–61. 10.7453/gahmj.2014.073 25694852 Mialich M. Martinez E. Jordão A. Jr. (2011). Comparative study of instruments for the analysis of body composition in a sample of the brazilian population. Int. J. Body Compos. Res. 9 19–24. 10.1186/s12913-016-1423-5 27409075 Piras A. Cortesi M. Campa F. Perazzolo M. Gatta G. (2019). Recovery time profiling after short-, middle-and long-distance swimming performance. J. Strength Cond. Res. 33 1408–1415. 10.1519/JSC.0000000000002066 28658076 Pla R. Aubry A. Resseguier N. Merino M. Toussaint J.-F. Hellard P. (2019). Training organization, physiological profile and heart rate variability changes in an open-water world champion. Int. J. Sports Med. 40 519–527. 10.1055/a-0877-6981 31288291 R Core Team (2021). R: A Language and Environment for Statistical Computing [Internet]. Vienna: R Foundation for Statistical Computing. RStudio Team (2021). RStudio: Integrated Development Environment for R. Boston, MA: RStudio, PBC. Available online at: http://www.rstudio.com/. Rossi F. E. Ricci-Vitor A. L. Sabino J. P. Vanderlei L. C. Freitas I. F. Jr. (2014). Autonomic modulation and its relation with body composition in swimmers. J. Strength Cond. Res. 28 2047–2053. 10.1519/JSC.0000000000000344 24345973 Rossi R. C. Vanderlei L. C. M. Gonçalves A. C. C. R. Vanderlei F. M. Bernardo A. F. B. Yamada K. M. H. (2015). Impact of obesity on autonomic modulation, heart rate and blood pressure in obese young people. Auton. Neurosci. 193 138–141. 10.1016/j.autneu.2015.07.424 26260435 Santana M. D. Kliszczewicz B. Vanderlei F. M. Monteiro L. R. Martiniano E. C. De Moraes Y. M. (2019). Autonomic responses induced by aerobic submaximal exercise in obese and overweight adolescents. Cardiol. Young 29 169–173. 10.1017/s1047951118002007 30785383 Santos-Magalhaes A. F. Aires L. Martins C. Silva G. Teixeira A. M. Mota J. (2015). Heart rate variability, adiposity, and physical activity in prepubescent children. Clin. Autono. Res. 25 169–178. 10.1007/s10286-015-0277-y 25820790 Schneider C. Wiewelhove T. Raeder C. Flatt A. A. Hoos O. Hottenrott L. (2019). Heart rate variability monitoring during strength and high-intensity interval training overload microcycles. Front. Physiol. 10:582. 10.3389/fphys.2019.00582 31178746 Seiler S. Haugen O. Kuffel E. (2007). Autonomic recovery after exercise in trained athletes: intensity and duration effects. Med. Sci. Sports Exerc. 39 1366–1373. 10.1249/mss.0b013e318060f17d 17762370 Spina G. Gonze B. Barbosa A. Sperandio E. Dourado V. (2019). Presence of age-and sex-related differences in heart rate variability despite the maintenance of a suitable level of accelerometer-based physical activity. Braz. J. Med. Biol. Res. 52:e8088. 10.1590/1414-431X20198088 31389489 Tanina H. Nishimura Y. Tsuboi H. Sakata T. Nakamura T. Murata K. (2017). Fatigue-related differences in erector spinae between prepubertal children and young adults using surface electromyographic power spectral analysis. J. Back Musculoskelet. Rehabil. 30 1–9. 10.3233/BMR-160705 27232083 Tarvainen M. P. Niskanen J.-P. Lipponen J. A. Ranta-Aho P. O. Karjalainen P. A. (2014). Kubios HRV–heart rate variability analysis software. Comput. Methods Programs Biomed. 113 210–220. Task Force of the European Society of Cardiology the North American Society of Pacing Electrophysiology (1996). Heart rate variability: standards of measurement, physiological interpretation and clinical use. Circulation 93 1043–1065. 10.1161/01.cir.93.5.1043 Triposkiadis F. Ghiokas S. Skoularigis I. Kotsakis A. Giannakoulis I. Thanopoulos V. (2002). Cardiac adaptation to intensive training in prepubertal swimmers. Eur. J. Clin. Investig. 32 16–23. 10.1046/j.0014-2972.2001.00939.x 11851722 Vachon A. Berryman N. Mujika I. Paquet J.-B. Arvisais D. Bosquet L. (2021). Effects of tapering on neuromuscular and metabolic fitness in team sports: a systematic review and meta-analysis. Eur. J. Sport Sci. 21 300–311. 10.1080/17461391.2020.1736183 32172680 Vandewalle H. Pérès G. Monod H. (1987). Standard anaerobic exercise tests. Sports Med. 4 268–289. Weinstein Y. Inbar O. Mor-Unikovski R. Luder A. Dubnov-Raz G. (2018). Recovery of upper-body muscle power after short intensive exercise: comparing boys and men. Eur. J. Appl. Physiol. 118 1555–1564. 10.1007/s00421-018-3885-0 29761312 Wickham H. (2016). ggplot2: Elegant Graphics for Data Analysis. New York, NY: Springer-Verlag. Available online at: https://ggplot2.tidyverse.org Xie Y. Dervieux C. Riederer E. (2020). R Markdown Cookbook. Boca Raton, FL: Chapman and Hall/CRC. Available online at: https://bookdown.org/yihui/rmarkdown-cookbook