Effects of Exercise Intensity on Acute Circulating Molecular Responses Poststroke

Abstract

Background. Exercise intensity can influence functional recovery after stroke, but the mechanisms remain poorly understood. Objective. In chronic stroke, an intensity-dependent increase in circulating brain-derived neurotrophic factor (BDNF) was previously found during vigorous exercise. Using the same serum samples, this study tested acute effects of exercise intensity on other circulating molecules related to neuroplasticity, including vascular-endothelial growth factor (VEGF), insulin-like growth factor-1 (IGF1), and cortisol, with some updated analyses involving BDNF. Methods. Using a repeated-measures design, 16 participants with chronic stroke performed 3 exercise protocols in random order: treadmill high-intensity interval training (HIT-treadmill), seated-stepper HIT (HIT-stepper), and treadmill moderate-intensity continuous exercise (MCT-treadmill). Serum molecular changes were compared between protocols. Mediation and effect modification analyses were also performed. Results. VEGF significantly increased during HIT-treadmill, IGF1 increased during both HIT protocols and cortisol nonsignificantly decreased during each protocol. VEGF response was significantly greater for HIT-treadmill versus MCT-treadmill when controlling for baseline. Blood lactate positively mediated the effect of HIT on BDNF and cortisol. Peak treadmill speed positively mediated effects on BDNF and VEGF. Participants with comfortable gait speed ≥0.4 m/s had significantly lower VEGF and higher IGF1 responses, with a lower cortisol response during MCT-treadmill. Conclusions. BDNF and VEGF are promising serum molecules to include in future studies testing intensity-dependent mechanisms of exercise on neurologic recovery. Fast training speed and anaerobic intensity appear to be critical ingredients for eliciting these molecular responses. Serum molecular response differences between gait speed subgroups provide a possible biologic basis for previously observed differences in training responsiveness.

Introduction

Endurance exercise can improve fitness, gait, balance, cognition, and cardiovascular risk factors after stroke. ¹ Accumu-lating longitudinal evidence suggests that a vigorous training intensity (eg, ≥60% heart rate reserve) is a critical ingredient for improving at least some of these outcomes (eg, fitness, gait, and balance),^2-5 and that pairing exercise with task-specific training may be essential for improving outcomes related to central nervous system function (eg, gait).^4,6 However, the potential mechanisms underlying an intensity-dependent effect of exercise on neurologic function remain poorly understood. Identifying these mechanisms could potentially lead to more effective rehabilitation interventions, greater clinical implementation, novel therapeutic targets (eg, for pharmacotherapy) and better prognostic measures to identify likely responders.^6,7

Lately, increasing attention has been focused on brain-derived neurotrophic factor (BDNF) as a possible mediator of the neurologic benefits of exercise. ⁶ BDNF is an abundant growth factor that is involved with activity-induced neuroplasticity,^6,8 is upregulated in the animal brain by exercise,^9-11 and shows acute increases in human blood during vigorous exercise, which have been correlated with enhanced cognitive and motor learning among healthy adults.^12,13 Among persons with chronic stroke, it was recently found that treadmill high-intensity interval training elicited a significantly greater acute BDNF increase than moderate intensity treadmill exercise, which showed no significant effect. ¹⁴ A secondary analyses found that an exercise intensity around the onset of blood lactate accumulation was the optimal threshold for discriminating large acute BDNF responses. ¹⁴

While these results provide preliminary evidence that BDNF could potentially mediate some intensity-dependent effects of exercise on neurologic function, there are also several other less-studied molecules that could be involved. For example, vascular endothelial growth factor (VEGF) and insulin-like growth factor 1 (IGF1) are 2 additional neurotrophins that accumulate in human blood during exercise,^15,16 and have some evidence of intensity-dependence.^15,17,18 VEGF is a protein involved with angiogenesis and neurogenesis, including in the brain after stroke. ¹⁹ Circulating VEGF increases during exercise are thought to be primarily triggered by hypoxia and vascular shear stress in skeletal muscles.^15,17 Similar peripheral increases in VEGF have been shown to promote perilesional angiogenesis and neurologic recovery in animal models of postacute stroke.^20,21 IGF1 is a protein involved with tissue remodeling, including in the brain, where it plays roles in neurogenesis, angiogenesis, neuroprotection, and neuroplasticity, ²² including after stroke. ¹⁹ Circulating IGF1 primarily originates in the liver, but is the main IGF1 source for the central nervous system. ²² Peripheral IGF1 administration has been shown to have a variety of positive effects on the brain in animal studies.^23,24 Further, acute exercise-induced increases in circulating levels of VEGF and IGF1 have each been shown to mediate hippocampal neurogenesis in animals.^25,26 Cortisol is another molecule that could mediate intensity-dependent effects of exercise on neurologic function, but one that could have negative effects. Circulating cortisol can acutely increase during psychological stress or high-volume strenuous exercise, due to activation of the hypothalamic-pituitary-adrenal axis, ²⁷ and can impair central neuroplasticity and learning if excessively upregulated.^28-31

It remains unknown whether acute upregulation of circulating BDNF, VEGF, IGF1, and/or cortisol meaningfully influences stroke rehabilitation outcomes. However, before costly longitudinal trials are conducted, 2 important steps are to (1) confirm whether exercise actually causes acute increases in circulating VEGF, IGF1, and/or cortisol among humans with stroke and (2) evaluate the influence of exercise intensity on these responses, to guide protocol selection. No previous studies have addressed these aims.

In addition, it remains poorly understood how exercise causes different acute circulating molecular responses. Mediation analysis can be used to assess pathways by which an intervention (eg, exercise protocol) leads to an outcome (eg, acute circulating BDNF response) through an intermediate or mediating variable (eg, blood lactate during training). ³² Developments in statistical methodology have made it possible to test for causal mediation with lower sample sizes ³³ and to directly quantify mediated and non-mediated effects in the same model. ³⁴ To our knowledge, this newer methodology for mediation analysis has not yet been applied to neurologic rehabilitation research.

The purpose of this study was to determine the effect of exercise intensity on acute circulating VEGF, IGF1, and cortisol responses poststroke, using the same data set and serum samples as a previous BDNF analysis. ¹⁴ We hypothesized that high-intensity interval training (HIT), whether performed on a treadmill (HIT-treadmill) or seated stepper (HIT-stepper), would elicit significantly greater acute increases in VEGF and IGF1 compared with moderate-intensity continuous exercise (MCT-treadmill). We also hypothesized that there would be no significant between-protocol difference in cortisol response, based on the limited evidence of psychological stress found in a previous qualitative assessment of treadmill HIT after stroke. ³ We also sought to evaluate mediating effects of different types of intensity and modifying effects of different clinical characteristics on these serum molecular responses. BDNF response was included as an additional dependent variable in the mediation analyses, because the novel statistical methods used herein provided more direct quantification of mediated effects than the previous BDNF analysis that used the same samples. ¹⁴

Methods

The data collection methods for this project have been more extensively reported elsewhere. ¹⁴ Pertinent information for this analysis is provided below.

Participants

This research was approved by the University of Cincinnati Institutional Review Board. Participants were recruited from the community between April and August of 2015 and provided written informed consent prior to participation. Inclusion criteria were the following: age 21 to 80 years; unilateral paresis from stroke experienced >6 months prior to enrollment (to help ensure stable deficits) ³⁵ ; residual gait impairment; able to walk with assistive devices as needed and no physical assistance; stable cardiovascular condition; discharged from formal rehabilitation; and able to communicate with investigators and correctly answer consent comprehension questions. Exclusion criteria were the following: evidence of significant arrhythmia or myocardial ischemia on treadmill exercise test ³⁶ ; hospitalization for cardiopulmonary disease within 3 months; pacemaker or implanted defibrillator; lower extremity claudication; severe lower limb spasticity, that is, Ashworth >2 ³⁷ ; weightbearing pain >4/10; and pregnancy.

Baseline Testing

Each participant had a clinical examination and a symptom-limited graded exercise test ³⁸ with electrocardiography in a clinical cardiac stress laboratory to determine eligibility and baseline characteristics. This was followed by a repeated exercise test with gas exchange analysis (on a separate day) to determine peak oxygen consumption rate and heart rate. The exercise test was repeated so that the participant could acclimate to exercise testing before wearing a facemask for gas exchange analysis during the test. During each exercise test, speed was held constant at the participant’s predetermined, individualized fastest comfortable speed, while incline was increased 2% to 4% every 2 minutes until the participant could not keep up, reached a cardiovascular safety limit,^36,39 or exhibited severe gait instability.

Study Design

We then used a within-participant repeated-measures experimental design. Each participant performed single sessions of 3 different exercise protocols in random order with approximately 1 week between sessions. Within participants, each session was scheduled at the same time of day to minimize random error from diurnal fluctuations. Participants were instructed to stay hydrated and to eat a meal approximately 2 hours prior to arrival for each session. All protocols were performed individually in a rehabilitation research lab by the same physical therapist (who was experienced in stroke rehabilitation), with support from an experienced phlebotomy nurse and research assistants for data collection. Personnel training was performed prior to data collection for approximately 2 months, including data collection practice for each protocol with 2 healthy participants.

Exercise Protocols

During all sessions, participants wore their habitual orthotic devices and were closely monitored and questioned to identify any adverse events. A fall protection harness (without weight support) and a height-adjusted handrail were used during treadmill sessions. Each protocol included a 3-minute warm-up at 30% to 50% heart rate reserve, 20 minutes of exercise (including recovery periods during HIT) and a 2-min cool-down.

Treadmill High-Intensity Interval Training

HIT-treadmill involved repeated 30-second bursts at maximum tolerated speed (0% incline), alternated with recovery periods where the treadmill was stopped. Recovery duration was decreased from 60 to 30 seconds after the first 5 minutes.^3,40 Initial burst speed was determined by a rapid acceleration test toward the end of the warmup, when speed was increased by 0.1 mph every 5 seconds until reaching the limit where the participant drifted backward or exhibited gait instability. Speed was continuously progressed, maintained or regressed between bursts based on performance criteria.^3,40 Participants were given verbal encouragement to stay to the front of the treadmill during each burst. Target mean intensity was >60% heart rate reserve. Unlike continuous exercise, the nature of this HIT protocol precludes the use of an instantaneous heart rate target during training. More specifically, the 30-second bursts are too brief for heart rate to reach steady state, so it fluctuates between burst and recovery periods, and has an upward trend over the session. ⁴⁰ However, this HIT protocol was designed to enable persons with stroke to sustain vigorous aerobic intensities (ie, >60% heart rate reserve) and has done so consistently.^3,14,40

Seated Stepper High-Intensity Interval Training

HIT-stepper used the same burst and recovery durations as HIT-treadmill, but on a seated stepper (NuStep T5XR). Bursts were performed at maximum possible cadence against 50% of maximal resistance (rounded down). Maximal resistance was predetermined with a steep ramp test ¹⁴ where resistance was increased by 1 level every 15 seconds until the participant could not maintain a cadence of 80 steps/min (40 cycles/min). During bursts, visual feedback about cadence and verbal encouragement were provided to increase or maintain maximal cadence.

Treadmill Moderate-Intensity Continuous Exercise

MCT-treadmill involved continuous treadmill walking with speed continuously adjusted to achieve a mean intensity of 45% ± 5% heart rate reserve, ⁴¹ with intermittent verbal encouragement.

Exercise Intensity Measures

Blood lactate was measured multiple times during each session (see below for details) and values obtained after 5, 10, and 20 minutes of exercise were averaged. Oxygen consumption rate was measured from expired air throughout each session using a Parvomedics TrueOne 2400 metabolic cart with a facemask interface. ⁴² Heart rate responses were captured by electrocardiography. Data were averaged across the 20 minutes of exercise, except for treadmill speed, which used the peak value for analysis to better reflect the sprint-like nature of interval training without omitting the recovery periods (which would make the measurement methods differ between HIT and MCT).

Blood Collection, Processing, and Analysis

An antecubital intravenous line was placed in the nonparetic arm on participant arrival for each session. Baseline samples were taken 30 minutes after intravenous insertion, to minimize any potential effect of intravenous insertion–related stress. For the exercise testing session, 10 mL blood samples were obtained at 3 baseline time points (T^-25, T^-20, T0) and 1 minute after the exercise test. For each of the exercise training sessions, 10 mL blood samples were taken before the warm-up (T0), after 5, 10, and 20 minutes of exercise (T5, T10 and T20), then 30 and 60 minutes after finishing the exercise (T50 and T80).

One drop from each sample was immediately tested with a point of care analysis system to measure blood lactate. ⁴³ The remaining sample was clotted and centrifuged to obtain serum, which was aliquoted and stored in a −80°C freezer. BDNF, VEGF, and IGF1 analyses were done in batches using Quantikine sandwich enzyme-linked immunosorbent assay kits (R&D Systems; Minneapolis, MN). Cortisol analysis was done in a batch with electrochemiluminescence immunoassays and a cobas e 411 analyzer (Roche Diagnostics; Basel, Switzerland). Laboratory staff were blinded to the study aims, design, and exercise protocol associated with the samples.

Statistical Analysis

Primary Hypothesis Testing

The primary analysis used a separate fixed-effects general linear model for each serum molecule to compare responses between protocols. The model included fixed effects for protocol, time, protocol × time, session number, session number × time, hemolysis, and participant ID (to account for within-participant repeated sessions), with an unconstrained covariance matrix to account for within-participant, within-session repeated measures. ⁴⁴ Protocol × time interaction contrasts were obtained to test the hypothesis that HIT-treadmill and HIT-stepper would each elicit significantly different changes compared with MCT-treadmill from T0 to the average of T5, T10, and T20.

Sample Size Determination

The software GLIMMPSE ⁴⁵ was used to estimate power to detect a significant protocol × time interaction with the above model, based on a 2-tailed significance level of .05. The expected effect size and covariance parameters were estimated from the lactate data and a preliminary BDNF data subset. (The data and aliquots from the preliminary BDNF assay were not reused for the full analysis.) This analysis indicated an expected power of 0.90 with 15 participants. Target sample size was 16 to account for potential missing samples.

Intensity Variable Mediation Analysis

Statistical mediation analyses were performed to better understand potential causal mechanisms ³² of any between-protocol differences (or lack thereof) in serum molecular responses. This approach uses a mediation model to test how much of an independent variable’s effect on a dependent variable occurs through (or is nullified by) an intermediate variable (mediator). While the validity of any causal inference depends on the appropriateness of the model used (eg, no residual confounding), statistical mediation analysis can provide relatively strong evidence of causality when used in combination with a prospective randomized experimental design. ⁴⁶ We tested mediating effects of 4 different exercise intensity variables on serum molecular responses, including mean blood lactate, peak treadmill speed (for the treadmill protocols only), mean oxygen consumption rate, and mean heart rate.

These analyses used natural effects models ³⁴ to estimate the mediated (indirect), nonmediated (direct), and total effects of protocol on serum molecular responses. For each protocol-mediator-response observation, an additional observation was imputed for each of the other 2 protocols, with the mediator level matched to the observed response. ⁴⁷ This imputation modeled molecular response as a function of the protocol, mediator, participant ID, baseline molecular concentration, and session number. A natural effects model was then fit to this expanded dataset, ³⁴ using molecular response as the dependent variable, with fixed effects for the mediated effect, nonmediated effect, participant ID, baseline molecular concentration and session number. Robust standard errors were calculated using the sandwich estimator. ⁴⁸

Clinical Characteristic Effect Modification Analysis

We also tested the modifying effects of different clinical characteristics on the training intensity variables and molecular responses. Tested characteristics were age, sex, stroke chronicity, body mass index, comfortable gait speed, comfortable gait speed subgroup (≥0.4 m/s or <0.4 m/s),^49,50 fast gait speed, peak oxygen consumption rate, and peak heart rate (from exercise testing). Each characteristic was tested separately with the intensity variable or molecular response as the dependent variable, fixed effects for the clinical characteristic, protocol and baseline molecular concentration, and an unconstrained covariance matrix. We also tested for characteristic by protocol interactions. Since comfortable gait speed ≥0.4 m/s has been shown to predict greater responsiveness to treadmill training, ⁴⁹ we also tested molecular responses to each protocol within each gait speed subgroup.

All analyses were done in SAS v9.4, except for the intraclass correlation coefficients and mediation analyses, which were done in R v3.5.1 with the “psych” ⁵¹ and “medflex” ⁴⁸ packages, respectively.

Results

Summary of Relevant Previously Reported Results ¹⁴

Sixteen participants were enrolled and all 16 completed the study. Participants were 57.4 ± 9.7 (mean ± SD) years old and 6.5 ± 4.1 years poststroke with a comfortable gait speed of 0.72 ± 0.33 m/s. Out of 352 total planned blood collections, 340 (96.6%) were successfully tested for lactate and 331 (94.0%) were obtained with sufficient volume for the current serum molecular analysis. Signs of hemolysis were observed in 71 (21.6%) of the 331 serum collections.

The HIT-stepper protocol was modified after 2 of the first 9 participants experienced near syncope during recovery from this protocol only. These events were thought to be partially related to higher than expected intensities enabled by the stepper (eg, lactate values of 17.0 and 12.8 mmol/L). Therefore, for the remaining 7 participants, we did not use resistance and enforced an 85% heart rate reserve limit for the HIT-stepper protocol. No other serious adverse events occurred, and the rate of nonserious events was similar between protocols.

HIT elicited significantly (P < .05) greater mean responses than MCT for blood lactate (HIT-treadmill, 4.6 mmol/L; HIT-stepper, 6.8 mmol/L; MCT-treadmill, 2.0 mmol/L), mean heart rate (HIT-treadmill, 59.0% heart rate reserve; HIT-stepper, 67.5%; MCT-treadmill, 43.8%), and peak treadmill speed (HIT-treadmill, 1.30 m/s; MCT-treadmill, 0.68 m/s). The revised HIT-stepper protocol (n = 7) had significantly lower responses than the original HIT-stepper protocol (n = 9) for lactate (original, 8.5 mmol/L; revised, 3.3 mmol/L) and mean heart rate (original, 79.1% heart rate reserve; revised, 43.6%), so we analyzed these protocols separately in a sensitivity analysis.

Acute Effects of Exercise on Serum Molecules

During the graded exercise test, each molecule showed a nonsignificant increase (Figure 1; Table 1). For the exercise protocols, VEGF significantly increased during HIT-treadmill only, with nonsignificant increases during HIT-stepper and MCT-treadmill. IGF1 significantly increased during HIT-treadmill and HIT-stepper, with a nonsignificant increase during MCT-treadmill. Cortisol showed nonsignificant decreases during all three protocols.

Figure 1.

Time course of changes in serum molecules by protocol (N = 16; sessions = 64; samples = 331). Results from the primary analysis (left) and a sensitivity analysis separating HIT-stepper into original and revised protocols (right), presented as change from T0, with 95% confidence interval error bars. GXT, graded exercise test; HIT, high-intensity interval training; MCT, moderate-intensity continuous training.

Table 1.

Acute Exercise Effects on Serum Molecules. ^a

	VEGF, pg/mL	IGF1, ng/mL	Cortisol, µg/dL
Treadmill GXT session (N = 16; sessions = 16; samples = 64)
GXT	+45.8 [−15.0, 106.5]	+1.4 [−2.2, 5.0]	+1.8 [−0.9, 4.4]
Primary analysis (N = 16; sessions = 48; samples = 267)
HIT-treadmill	+49.2 [8.2, 90.2]	+3.7 [0.8, 6.6]	−0.5 [−1.5, 0.6]
HIT-stepper	+10.1 [−30.0, 50.2]	+4.2 [1.4, 7.1]	−0.6 [−1.7, 0.4]
MCT-treadmill	+10.5 [−31.3, 52.3]	+2.6 [−0.3, 5.6]	−0.5 [−1.6, 0.6]
HIT-treadmill vs MCT-treadmill	+38.7 [−19.8, 97.2]	+1.0 [−3.1, 5.2]	+0.0 [−1.5, 1.6]
HIT-stepper vs MCT-treadmill	−0.4 [−59.0, 58.2]	+1.6 [−2.5, 5.7]	−0.2 [−1.7, 1.4]
Comparison of HIT-stepper protocols (N = 16; sessions = 48; samples = 267)
HIT-stepper-original (N = 9)	+24.4 [−29.5, 78.2]	+3.8 [−0.0, 7.6]	−0.6 [−2.0, 0.8]
HIT-stepper-revised (N = 7)	−9.3 [−73.1, 54.5]	+4.8 [0.3, 9.3]	−0.8 [−2.4, 0.9]
Revised vs original	−33.6 [−118.6, 51.3]	+1.0 [−5.0, 7.1]	−0.2 [−2.4, 2.1]
Subgroup analysis: (N = 16; sessions = 48)
Comfortable gait speed ≥0.4 m/s
HIT-treadmill	+43.8 [5.8, 81.9]	+7.0 [3.7, 10.2]	−0.8 [−1.7, 0.2]
HIT-stepper	+28.9 [−9.6, 67.5]	+6.3 [3.0, 9.6]	−0.6 [−1.6, 0.4]
MCT-treadmill	−7.7 [-44.1, 28.7]	+4.9 [1.7, 8.1]	–2.0 [−3.0, −1.0]
Comfortable gait speed <0.4 m/s
HIT-treadmill	+117.1 [34.2, 199.9]	−0.9 [−7.6, 5.7]	1.0 [−1.1, 3.1]
HIT-stepper	+28.1 [-103.2, 159.5]	−0.6 [−8.1, 6.8]	−1.0 [−3.2, 1.2]
MCT-treadmill	+5.7 [-54.6, 66.1]	0.8 [−5.3, 6.9]	3.0 [1.1, 4.8]

Abbreviations: GXT, graded exercise test; VEGF, vascular-endothelial growth factor; IGF1, insulin-like growth factor 1; HIT, high-intensity interval training; MCT, moderate-intensity continuous training.

a

Values are estimated mean molecular responses [95% CI] with significant (P < .05) values in boldface. For the GXT, models included fixed effects for time and hemolysis with a compound symmetry covariance matrix to account for within-participant repeated measures. Biomarker response was the time contrast from T0 to post-GXT. For the primary analysis and comparison of HIT-stepper protocols, models included fixed effects for protocol, time, protocol × time, session number, session number × time, hemolysis, and participant ID, with an unconstrained covariance matrix to account for within-participant, within-session repeated measures. Biomarker response was the time contrast from T0 to the mean of T5, T10, and T20. For the subgroup analysis, molecular response was the mean of T5, T10, and T20 minus T0 and models included fixed effects for subgroup, protocol, subgroup × protocol, and baseline molecular concentration, with an unconstrained covariance matrix to account for within-participant repeated sessions.

There were no significant between-protocol differences for any molecular response in the a priori primary analysis (Table 1). However, T0 VEGF was a significant predictor of VEGF exercise response (β = −0.03 [−0.06, −0.01]), and adjusting for T0 VEGF made the HIT-treadmill versus MCT-treadmill response comparison statistically significant (+32.8 [9.9, 55.7]). T0 level was not a significant predictor of exercise responses for IGF1 (β = 0.06 [−0.02, 0.14]) or cortisol (β = −0.13 [−0.28, 0.02]) and adjusting for T0 level did not make any between-protocol comparisons statistically significant for these molecules.

Biomarker Responses to Original and Revised HIT-Stepper Protocols

In the sensitivity analysis separating the original and revised HIT-stepper protocols, there were no significant between-protocol differences in serum molecular responses (Figure 1 right; Table 1).

Mediating Effects of Different Intensity Measures on Between-Protocol Differences in Molecular Responses

Mediating Effects of Blood Lactate

The higher lactate levels during HIT-treadmill and HIT-stepper (vs MCT-treadmill) significantly increased BDNF responses by an estimated +1.53 and +3.15 ng/mL, respectively (Figure 2A). Outside of this lactate-mediated effect, other between-protocol differences nonsignificantly increased BDNF response by an estimated +0.62 ng/mL for HIT-treadmill (vs MCT-treadmill) and significantly decreased BDNF response by an estimated −4.14 ng/mL for HIT-stepper (vs MCT-treadmill). Consequently, the total effect estimate was a significantly greater BDNF response for HIT-treadmill versus MCT-treadmill (+2.15 ng/mL) and no significant difference in BDNF response for HIT-stepper versus MCT treadmill (−0.98 ng/mL).

Figure 2.

Intensity mediators of between-protocol differences in serum molecular responses. Values are model estimates [95% CI] in the units of the molecular response. A mediated effect (arrow going through a dashed box) is the estimated amount of molecular response that is attributable to between-protocol differences in lactate (panels A-D) or peak speed (panels E-H). A nonmediated effect (arrow not going through a dashed box) is the estimated amount of molecular response that is attributable to between-protocol differences other than lactate (panels A-D) or peak speed (panels E-H). A total effect (shown in upper left corner of each panel) is the estimated between-protocol difference in molecular response from the model and is the sum of the mediated and nonmediated effects. For example, panel A shows that the higher lactate levels for HIT-treadmill versus MCT-treadmill increased the BDNF response by an estimated +1.53 ng/mL, while other between-protocol differences nonsignificantly increased the BDNF response by an estimated +0.62 ng/mL, making the total estimated between-protocol difference +2.15 ng/mL. Biomarker response was the mean of T5, T10, and T20 minus T0. Blood lactate was the mean of T5, T10, and T20. Covariates included protocol, session number, baseline molecular concentration and participant ID. Controlling for participant ID effectively adjusts for each participant’s mean mediator level and molecular response, so the results are based on within-participant between-protocol differences. *Statistically significant (P < .05) effect. BDNF, brain-derived neurotrophic factor; VEGF, vascular-endothelial growth factor; IGF1, insulin-like growth factor 1.

The effect of HIT (vs MCT) on VEGF was not significantly mediated through lactate (Figure 2B). For IGF1, the higher lactate levels during HIT (vs MCT) decreased responses (significant for HIT-stepper only; Figure 2C). Outside of this negative lactate-mediated effect, other between-protocol differences significantly increased IGF1 responses for HIT, resulting in no significant total difference in IGF1 response for either HIT protocol compared with MCT-treadmill.

For cortisol, the higher lactate levels during HIT significantly increased responses (Figure 2D). Outside of this lactate-mediated effect, other between-protocol differences decreased cortisol responses for HIT (significant for HIT-stepper only), resulting in no significant total difference in cortisol response for either HIT protocol compared with MCT-treadmill.

Mediating Effects of Peak Treadmill Speed

The higher peak treadmill speeds for HIT-treadmill versus MCT-treadmill significantly increased BDNF responses by an estimated +4.84 ng/mL (Figure 2E) and VEGF responses by an estimated +24.56 pg/mL (Figure 2F). Outside of these speed-mediated effects, other between-protocol differences decreased molecular responses for HIT (significant for BDNF only), but the total effects from these models still indicated that HIT-treadmill increased BDNF and VEGF significantly more than MCT treadmill. Peak speed did not significantly mediate the effects of HIT (vs MCT) on IGF1 or cortisol (Figure 2G and (H).

Mediating Effects of Oxygen Consumption Rate and Heart Rate

Mean oxygen consumption rate and mean heart rate did not significantly mediate the effects of protocol on any serum molecular responses.

Effect Modification Based on Clinical Characteristics

Modifying Effects of Clinical Characteristics on Training Intensity

Greater blood lactate during training was significantly associated with greater stroke chronicity for HIT-stepper only (Table 2). Faster peak treadmill speed during training was significantly associated with faster baseline gait speeds, female sex, and higher peak oxygen consumption rate, but the associations with female sex and peak oxygen consumption rate were confounded by differences in baseline comfortable gait speed (Table 2). Mean oxygen consumption rate and mean heart rate during training were not significantly modified by any of the tested baseline characteristics, except that higher mean oxygen consumption rate was significantly associated with greater stroke chronicity for HIT-stepper only (like blood lactate).

Table 2.

Modifying Effects of Clinical Characteristics on Training Intensity and Serum Molecular Responses (N = 16; Sessions = 48). ^a

Clinical Characteristics		Modifying Effects on Intensity		Modifying Effects on Serum Molecular Responses
		Mean Lactate, mmol/L	Peak Treadmill Speed, m/s	VEGF, pg/mL	IGF1, ng/mL	Cortisol, µg/dL
Age (years)	57.4 ± 9.7	0.0 [−0.0, 0.1]	−0.0 [−0.0, 0.0]	0.2 [−1.6, 1.9]	0.0 [−0.2, 0.2]	0.1 [−0.0, 0.1]
Female sex	7 (43.8%)	−0.3 [−1.8, 1.2]	0.4 [0.1, 0.8] b	−14.5 [−47.8, 18.7]	−2.2 [−5.2, 0.8]	−1.3 [−3.0, 0.4]
Years post stroke	6.5 ± 4.1	0.2 [0.0, 0.3] c	−0.0 [−0.1, 0.0]	−1.6 [−5.6, 2.4]	−0.1 [−0.5, 0.3]	−0.0 [−0.2, 0.2]
Body mass index (kg/m2)	27.6 ± 3.7	−0.1 [−0.3, 0.1]	0.0 [−0.1, 0.1]	−0.3 [−5.6, 5.1]	−0.0 [−0.6, 0.5]	0.0 [−0.2, 0.3]
Comfortable gait speed (m/s)	0.72 ± 0.33	−0.6 [−2.9, 1.8]	0.9 [0.5, 1.2]	−37.6 [−77.9, 2.8]	5.1 [0.2, 10.1]	–3.5 [−5.5, −1.6] d
Comfortable gait speed ≥0.4 m/s	13 (81.3%)	−0.9 [−2.8, 1.0]	0.6 [0.1, 1.0]	−30.4 [−59.2, −1.5]	6.4 [3.4, 9.4]	–2.7 [−4.1, −1.4] a,d
Fast gait speed (m/s)	0.90 ± 0.41	−0.9 [−2.7, 1.0]	0.7 [0.4, 1.0]	−22.6 [−57.8, 12.6]	3.4 [−0.7, 7.5] c	–2.7 [−4.4, −0.9] d
VO2peak (mL/kg/min)	17.2 ± 3.3	0.1 [−0.1, 0.3] d	0.1 [0.0, 0.1] b	−3.0 [−7.8, 1.8]	0.4 [−0.1, 0.9]	−0.2 [−0.4, 0.1] d
HRpeak (bpm)	139 ± 21	0.0 [−0.0, 0.0]	0.0 [−0.0, 0.0]	0.5 [−0.2, 1.2]	−0.0 [−0.1, 0.0]	0.0 [−0.0, 0.1]

Abbreviations: VEGF, vascular-endothelial growth factor; IGF1, insulin-like growth factor 1; VO₂, oxygen consumption rate; peak, maximum value reached during symptom-limited, graded exercise test; HRR, heart rate reserve.

a

Values are mean ± SD, N (%), or estimated mean [95% CI] with significant (P < .05) values in boldface. Lactate is the mean of T5, T10 and T20. Biomarker response is the mean of T5, T10 and T20 minus T0. Analyses were adjusted for baseline molecular concentration (serum molecular analyses only) and protocol and included an unconstrained covariance matrix to account for within-participant repeated sessions.

b

Effects found to be confounded by gait speed. Regression coefficients for female sex and VO_2peak decreased by 85% and 50%, respectively, and became nonsignificant after adding comfortable gait speed to the model, which remained a significant factor in both cases.

c,d

Significant (P < .05) protocol × characteristic interaction. Effect modification significant for ^cHIT-stepper only or ^dMCT-treadmill only.

Modifying Effects of Clinical Characteristics on Serum Molecular Responses

Baseline comfortable gait speed category was a significant effect modifier for VEGF, IGF1, and cortisol responses and explained more variability in these responses than any other characteristic, including continuous measures of gait speed (Table 2). Compared with participants who had gait speed <0.4 m/s, participants with gait speed ≥0.4 m/s had significantly lower VEGF responses and higher IGF1 responses, which were not protocol specific, and significantly lower cortisol responses for MCT-treadmill only.

VEGF significantly increased during HIT-treadmill for both gait speed subgroups (Table 1), but the magnitude was significantly greater for the low gait speed subgroup. IGF1 only increased significantly for the higher gait speed subgroup, and these increases were significant for all three protocols. Cortisol significantly decreased during MCT-treadmill for the higher gait speed subgroup, but significantly increased during MCT-treadmill for the low gait speed subgroup.

Discussion

Treadmill HIT was recently found to elicit significantly greater increases in serum BDNF and corticospinal excitability compared with conventional treadmill MCT. ¹⁴ Using the same serum samples, the current study tested the intensity dependence of additional molecular responses related to neuroplasticity, and incorporated novel mediation and effect modification analyses. For the first time among humans with stroke, we now report significant acute increases in serum VEGF and IGF1 during exercise, without a concomitant cortisol stress response, using treadmill HIT. A significant IGF1 increase was also found during seated stepper HIT. Neither VEGF nor IGF1 significantly increased during conventional treadmill MCT (for the study sample as a whole).

Adding to the previous report, ¹⁴ serum BDNF response also showed additional evidence of intensity dependence in the current analysis, as it was positively mediated by both blood lactate and peak treadmill speed during training. Aside from BDNF, VEGF showed the strongest evidence of intensity-dependence. Although the between-protocol difference did not exceed the statistical significance threshold in the primary (unadjusted) analysis, VEGF response was over 300% greater for HIT-treadmill versus MCT-treadmill. Further, this difference was statistically significant when controlling for baseline VEGF level and was significantly mediated by faster treadmill speeds during HIT. Relative to BDNF and VEGF, IGF1 had less evidence of intensity dependence. While the IGF1 response was 42% greater for HIT-treadmill versus MCT-treadmill, this difference was not statistically significant and was not positively mediated by any intensity variable. In fact, higher lactate values during HIT were found to have decreased the IGF1 response. When statistically controlling for lactate effects, both HIT protocols elicited significantly greater IGF1 responses than MCT through other unknown mechanisms, but these between-protocol differences were cancelled out by the negative lactate effects. Thus, BDNF and VEGF appear to be the most promising serum molecules tested for potentially explaining some of the intensity-dependent effects of exercise on neurologic recovery from stroke. Future longitudinal trials are now warranted to determine the extent to which these molecules mediate the effects of exercise on rehabilitation outcomes.

Human stroke studies have consistently shown significant benefits or at least trends in neurologic function outcomes favoring higher intensity exercise over moderate- or low-intensity exercise.^2-5 Conversely, nonhuman studies have sometimes found worse outcomes with high-intensity exercise.^52-56 In these animals, high intensity exercise also typically elicits significant stress-related elevations in serum corticosterone (the animal analog of cortisol)^30,56 and aggression. ⁵⁷ While small increases in cortisol or corticosterone can have positive effects on central neuroplasticity, ²⁹ excessive elevations have been shown to have negative effects.^28,30,31 The lack of cortisol stress response during human post-stroke HIT in the current study could explain some of these observed between-species differences in optimal intensity. We speculate that humans may have less psychological stress during high-intensity exercise because of voluntary participation and easier communication between researchers and participants.

However, some previous human (nonstroke) studies have observed significant cortisol increases during exercise, generally using continuous and vigorous exercise protocols lasting at least 60 minutes, or shorter protocols with very high intensity.^27,58 In the current study, higher lactate during HIT (vs MCT) did appear to mediate a greater cortisol response, but this between-protocol difference was cancelled out by other factors (potentially the recovery intervals during HIT ⁵⁹ ). In the subgroup analysis, participants with comfortable gait speed <0.4 m/s had a relatively small but significant cortisol increase during treadmill MCT, but not HIT. It is possible that treadmill exercise provoked anxiety for these participants and that the recovery intervals during HIT helped mitigate that stress, despite the faster speeds.

Persons with stroke who have baseline comfortable gait speed ≥0.4 m/s have been previously shown to have greater improvements in walking capacity from continuous treadmill training. ⁴⁹ In the current study, this higher gait speed subgroup had greater serum IGF1 response and a lower serum cortisol response. This raises the hypothesis that greater serum IGF1 responses and/or lower serum cortisol responses could be related to greater training responsiveness. Contrary to our expectations, these differences in serum molecular responses between gait speed subgroups were not due to differences in training speed (see Supplemental Methods and Results), and the mechanisms remain unknown.

Another surprise was that the low gait speed subgroup had a 167% greater serum VEGF response to HIT-treadmill. Based on animal research, an acute exercise-induced increase in circulating VEGF could have beneficial effects on vascular and cognitive outcomes. ²⁵ If an acute circulating VEGF response has similar benefits in humans, and if those effects are dose dependent, it is possible that persons with low gait speed could have preferential vascular or cognitive benefits from HIT.

One issue with the present subgroup analysis is that there were only 3 participants with gait speed <0.4 m/s. While the repeated measures helped provide sufficient power to detect significant subgroup effects, the response estimates for this subgroup may be particularly prone to fluctuate in future studies. Still, significant effect modification was also present in the study sample as a whole when modeling gait speed as a continuous variable. Furthermore, the dichotomous gait speed subgrouping explained more of the variance in serum molecular responses than the continuous measure. This provides novel biochemical evidence of an important comfortable gait speed threshold around 0.4 m/s.

Limitations and Future Studies

The primary limitation to this study is that we did not test the effects of the protocols or serum molecular responses on learning or longitudinal outcomes. Previous nonstroke studies have reported exercise-induced learning enhancement with similar or lesser serum molecular responses than those we found with HIT-treadmill.^12,13 However, it is still unclear how much different types of learning can be enhanced through upregulation of BDNF, VEGF, and IGF1. The current results provide foundational evidence that HIT-treadmill and MCT-treadmill protocols would be useful to include in future longitudinal trials to test mediating effects of BDNF, VEGF, and possibly IGF1 on rehabilitation outcomes.

Future stroke studies might also consider testing the mediating effects of circulating catecholamines (eg, dopamine, epinephrine, norepinephrine) on exercise-induced learning enhancement. Catecholamines were not tested in this study because they must be measured in plasma (not serum), which would have increased the blood collection volume and time. However, catecholamine concentrations have been shown to acutely rise during exercise among healthy adults and their postexercise concentrations have been correlated with learning enhancement.^12,13 Future studies are also needed to provide a stronger base of normative data for comparison to the current results. Most previous studies testing the acute effects of exercise on circulating neurotrophins have recruited young healthy adults, and it is still unclear whether large age differences or stroke might affect these responses.

We acknowledge that serum assays provide imperfect measures of central nervous system processes. ¹³ For example, acute exercise-related brain changes in animals typically outlast corresponding increases in circulating molecular concentrations,^10,11 so the lack of sustained serum molecular changes in Figure 1 should be interpreted with caution. Brain levels of these molecules would likely be more related to learning but are not obtainable noninvasively.

For the mediation models to be truly causal, there cannot be any confounding factors affecting both the mediator (eg, lactate) and outcome (eg, BDNF change), ³⁴ apart from cofactors included in the model (protocol, session number, baseline molecular concentration, and participant ID). We cannot completely rule out this possibility but believe that uncontrolled confounding was unlikely to be a major problem because the models were able to discriminate between different mediators for different molecules in ways that are biologically plausible. For example, lactate infusion at rest has been shown to increase circulating BDNF, ⁶⁰ and the models identified blood lactate, but not oxygen consumption rate or heart rate, as a positive mediator of BDNF response.

Conversely, the mediation analysis also found that blood lactate positively mediated cortisol response and there is no established causal pathway through which lactate directly upregulates circulating cortisol. However, correlations between blood lactate and serum cortisol responses have been previously observed among healthy adults, ⁵⁹ anaerobic exercise has a greater effect on cortisol than aerobic exercise, ⁶¹ and it has been hypothesized that circulating lactate could act as a chemical messenger to stimulate the hypothalamic-pituitary-adrenal axis. ⁶² Thus, it is possible that this finding reflects true causal mediation.

The mediation analysis also does not rule out the possibility of additional upstream, downstream or parallel mediators. For example, we speculate that faster treadmill training speed may have mediated greater VEGF response because faster speeds cause greater blood pressure fluctuations and vascular shearing forces. ⁶³

In a previous analysis, ¹⁴ BDNF response was not significantly associated with any clinical characteristics. However, in the current subgroup effect mediation analysis of the same data (which included additional covariates and only the treadmill protocols), the low comfortable gait speed subgroup showed a significantly greater total BDNF response (see Supplemental Methods and Results). It is possible that this subgroup effect was somehow obscured by confounding in the previous analysis, and this could be explored in future studies.

Conclusions

A treadmill HIT protocol can elicit acute increases in circulating BDNF, VEGF, and IGF1 without a significant cortisol stress response, among ambulatory persons with chronic stroke. The faster treadmill speeds and higher blood lactate levels reached during HIT appear to be critical ingredients for eliciting BDNF and VEGF responses. A similar HIT protocol on a seated stepper enabled even higher intensities and similar BDNF and IGF1 responses but needs further testing to balance safety and benefit. BDNF and VEGF appear to be the most promising serum molecules to include in future studies testing intensity-dependent mechanisms of exercise on neurologic recovery from stroke. We also found significantly greater IGF1 responses and lesser cortisol responses for participants with comfortable gait speed >0.4 (vs ⩽ 0.4), leading to the hypothesis that these molecular responses could explain previously reported differences in walking capacity responsiveness between these subgroups.

Supplemental Material

Supplemental_Methods_and_Results – Supplemental material for Effects of Exercise Intensity on Acute Circulating Molecular Responses Poststroke

Supplemental material, Supplemental_Methods_and_Results for Effects of Exercise Intensity on Acute Circulating Molecular Responses Poststroke by Pierce Boyne, Colleen Meyrose, Jennifer Westover, Dustyn Whitesel, Kristal Hatter, Darcy S. Reisman, Daniel Carl, Jane C. Khoury, Myron Gerson, Brett Kissela and Kari Dunning in Neurorehabilitation and Neural Repair