Animals and procedures
Fifty-four SHR 18-week-old subjects were divided into two treatment groups: untrained (n = 27) and trained (n = 27). Each group was subdivided into three smaller groups (n = 9): vehicle (H2O), pyridostigmine bromide (Sigma-Aldrich, Saint Louis, MO, USA) diluted in drinking water at a dose of ~5 mg/kg/day, and pyridostigmine bromide diluted in drinking water at a dose of ~15 mg/kg/day for 2 weeks, between the 11th and 12th weeks of aerobic physical training. The dose was determined based on the results of a previous study 9. During the experiments, the animals were housed in the Animal Facility at the Ribeirão Preto Medical School, which was maintained at 23°C and 60%–70% humidity. The rats were kept on a 12/12-hour light/dark cycle and had free access to food and water. The experimental protocols used in the present study were in accordance with the ethical principles of animal experimentation adopted by the Brazilian College of Animal Experimentation and were evaluated and approved by the Animal Experimentation Ethics Committee (CETEA) of the Ribeirao Preto Medical School, University of Sao Paulo (Protocol 035/2014). The study was carried out in compliance with the ARRIVE guidelines.
On day 15 of treatment (with pyridostigmine bromide 5 or 15 mg/kg or only H2O), polyethylene catheters (PE-10/PE-50, Intramedic; Becton Dickinson and Company, Sparks, MD, USA) were implanted into the left femoral artery and vein under ketamine (5 mg/kg, intraperitoneal [i.p.]; Sigma-Aldrich, USA) and xylazine (30 mg/kg, intraperitoneal [i.p.]; Sigma-Aldrich, USA) anesthesia to record pulsatile BP and to administer drugs, respectively. The catheters were subcutaneously tunneled and exteriorized in the nape. Twenty-four hours after the surgical procedures, BP was measured in conscious rats kept in a quiet environment. The BP was recorded using a pressure transducer (MLT0380; AD Instruments, Bella Vista, Australia), and the amplified signal (ML110; AD Instruments, Bella Vista, Australia) was fed to a computer acquisition system (PowerLab 8/30; AD Instruments, Bella Vista, Australia). The mean BP (MBP) and HR were calculated from the pulsatile BP.
Physical Training
The rats in the training groups underwent a protocol of aerobic physical training that consisted of swimming sessions in a glass tank (100 cm long × 80 cm wide × 80 cm high), which allowed for the simultaneous training of six animals. The tank was filled with 50 cm of warm water (30 ± 2°C), which was changed after every group training session. The training program was conducted in two different stages over a total of 12 (from 18 to 30 weeks of age). The first stage consisted of a two-week adaptation period, during which the session length was gradually increased from 5 to 30 min per day, five times per week (in increments of 5 min per day). The second stage consisted of 10 weeks, with 30 min of physical training sessions conducted five times per week15. To evaluate physical training intensity, blood was collected from the tail veins of the animals at the fourth, seventh, and tenth weeks immediately before and after the 30-minute exercise sessions, and the lactate concentration was measured (Accutrend® Plus, Roche Diagnostics, Mannheim, Germany). The expected lactate level ranged from 5.5 to 6 mmol/L, as previously determined. If the animals did not achieve the expected lactate concentration, the level of training exertion was increased by fastening a leaded, impermeable Velcro strap to the chest to increase body weight by 2%–6% 16.
Echocardiography
At 30 weeks of age, all animals underwent an echocardiographic evaluation. We used a Vevo 2100® High-Resolution Imaging System ultrasound (VisualSonics, Toronto, ON, Canada) instrument with a high-resolution transducer (21 MHz). For the procedure, the anterior regions of the thorax were previously trichotomized (Veet®, Reckitt Benckiser, São Paulo, SP, Brazil), and all animals were anesthetized with 1.5% isoflurane supplemented with 1% O2 and placed on a heated (37°C) platform. Echocardiography and temperature measurements were monitored.
High-resolution B-mode and M-mode images were acquired. Wall thickness and left ventricle dimensions were obtained from a short-axis view at the level of the papillary muscles. Diastolic measurements were performed at the point of maximum cavity dimension, and systolic measurements were performed at the point of minimal cavity dimension. All measurements were performed according to the standards of the American Society of Echocardiography and by an evaluator who was blinded to which group the rats were assigned at the time of measurement 17. The following parameters were obtained from the images: interventricular septum thickness (IVST), posterior wall thickness (PWT), end-diastolic diameter of the left ventricle (LVEDD), and end-systolic diameter of the left ventricle (LVESD). The shortening fraction was calculated as FS lrb% = [(LVEDD-LVESD÷LVEDD) × 100], and the ejection fraction (EF) was calculated using the Teichholz method: [(LVEDV-LVESV÷LVEDV) × 100]. The left ventricle mass (LV mass/final body weight) was calculated using the formula: 1.047 × [(LVEDD+PWT+IVST)3 - LVEDD3], and the relative wall thickness (RWT) was calculated as follows: [2×PWT÷LVEDD]. Left ventricular volumes were quantified using the following formula: LVEDV (µL) = [LVEDD3×(7÷2.4+LVEDD3)] and LVESV (µL) = [LVESD³ × (7÷2.4+LVESD3)] 15,18,19.
Surgical procedure
Forty-eight hours after echocardiography, the rats were anesthetized with ketamine (80 mg/kg) and xylazine (10 mg/kg), and polyethylene catheters made in our laboratory (PE-50 spliced by melting to PE-10; Intramedic, Clay Adams, Parsippany, New Jersey, USA) were implanted into the left femoral artery and vein. Catheters were tunneled subcutaneously and exteriorized at the nape. To prevent blood clotting, the catheters were filled with heparinized saline solution (500 IU/mL). The rats were then allowed to recover for 24 hours before the cardiac sympathovagal assessment protocol, which was performed without anesthesia.
Experimental protocols
Heart rate variability and systolic blood pressure variability
Twenty-four hours after the surgical procedure, arterial pulse (AP) pressure was measured in conscious rats kept in a quiet environment. BP was recorded using a pressure transducer (MLT0380, ADInstruments). Additionally, the amplified signal (ADInstruments ML110) was fed to a computer acquisition system (LabChart 7 Pro). MBP and heart rate (HR) were calculated from the arterial pulse pressure.
HR (pulse interval) and systolic BPV analysis were performed using custom computer software (CardioSeries v2.4, Dias, DPM, University of São Paulo, Brazil http://sites.google.com/site/cardioseries) 20. The software was designed to perform time-frequency analysis of cardiovascular variability, allowing precise adjustment of parameters related to frequency domain analysis (e.g., interpolation rate, segment length, and boundaries of frequency bands). Because the software does not perform data sampling, a beat-by-beat time series must be generated and loaded into the CardioSeries software. The BP baseline and pulse interval (PI) series, obtained from 60-minute recordings, were processed using computer software (LabChart v7.0, ADInstruments, Bella Vista, Australia) that applies an algorithm to detect cycle-to-cycle inflection points in the pulsatile BP signal, thus determining beat-by-beat values of systolic blood pressure. Beat-by-beat PI series were generated from the pulsatile BP signal by measuring the time interval between adjacent systolic peaks. Next, a beat-by-beat series of PI and SBP were converted to data points every 100 ms using cubic spline interpolation (10 Hz). The interpolated series were divided into half-overlapping sequential sets of 512 data points (51.2 seconds), which were tested for stationarity. It is important to mention that cardiovascular variability analysis requires at least a weakly stationary data series (i.e., mean and covariance stable over time; Berntson et al., 1997). Data series stationarity can be verified by means of stationarity tests (i.e., enhanced reproducibility of the results among users and laboratories 22, as well as through visual inspection of the data series 23,24.
In our study, a well-experienced researcher visually inspected the segments of interpolated time series (i.e., PI or SBP values) searching for transients that could affect the calculation of the power spectral density (PSD). To confirm that the visual inspection of the time series was properly performed, a Hanning window was used to attenuate side effects, and the spectrums of all segments were calculated using a direct fast Fourier transform algorithm for discrete time series. All segments were visually inspected for abnormal spectra. Finally, considering the results from the time series and spectra inspections, non-stationary data were not taken into consideration for PSD calculation. The spectra were integrated in low-frequency (LF; 0.2–0.75 Hz) and high-frequency (HF; 0.75–3 Hz) bands, and the results were expressed in absolute (ms2 or mmHg2) and normalized units (nu). The normalized values were created by calculating the percentage of LF and HF power of the total spectrum power minus the very low frequency band (VLF; < 0.2 Hz) and power 24,25. To assess the sympathovagal balance, the LF/HF ratio of PI variability was calculated 26.
Spontaneous baroreflex sensitivity
Baroreflex sensitivity (BRS) was assessed in the time domain using the sequence technique, as described by Di Rienzo et al. (1985). Custom computer software (CardioSeries v2.4, http://sites.google.com/site/cardioseries) scanned the beat-by-beat time series of SBP and PI, searching for sequences of at least four consecutive beats in which increases in SBP were followed by PI lengthening (up sequence) and decreases in SBP followed by PI shortening (down sequence) with a linear correlation higher than 0.8. The slope of the linear regression lines between SBP and PI was used as a measure of spontaneous BRS 28.
Assessment of cardiac sympathovagal balance
The influence of sympathetic and parasympathetic autonomic tone on HR was assessed by administering propranolol (5 mg/kg, intravenous [i.v.], Sigma-Aldrich, USA) and methylatropine (4 mg/kg, i.v.; Sigma-Aldrich, USA), respectively. For this purpose, the femoral artery catheter was attached to a pressure transducer (MLT844, AD Instruments, Bella Vista, Australia), which converts the AP fluctuations into electrical signals. Next, signals were amplified using a bridge amplifier (FE117, AD Instruments, Bella Vista, Australia), and pulsatile AP was continuously sampled (2 kHz) using a computer equipped with an analog-digital interface (ML866, AD Instruments, Bella Vista, Australia). After 60 minutes of basal HR recording, methylatropine was injected into half of the rats in each group, and HR was recorded for the following 15 min to assess the effect of vagal blockade on HR. Propranolol was then injected into the same rats, and HR was recorded for another 15 min to determine the intrinsic HR (IHR). In the other half of the rats, the methylatropine–propranolol sequence was reversed to assess the effect of sympathetic blockade on HR, following the same recording procedure (15 min each) for each drug, as described previously, to determine the IHR. The data from methylatropine–propranolol and propranolol–methylatropine sequences were pooled to provide the basal HR (before any drugs) and the IHR (after drugs).
Statistical analysis
The results are presented as mean ± standard error of the mean (SEM). The effects of hypertension and pharmacological treatments were assessed using two-way analysis of variance (ANOVA). When appropriate, post-hoc comparisons were performed using the Student-Newman-Keuls test. For comparison between two groups, the Student’s t-test for independent measures or the Mann-Whitney Rank Sum test was used as required. Differences were considered significant at P < 0.05. All statistical tests were performed using SigmaStat software (version 3.5; Systat Software Inc., San Jose, CA, USA).