Abstract
We performed an exploratory analysis to evaluate the effects of a treadmill exercise program on brachial artery (BA) intima–media thickness (IMT) and three BA grayscale ultrasound measures that may indicate subclinical arterial injury. Data were from a clinical trial in individuals with peripheral artery disease who were randomly assigned to treadmill exercise training or attention control. B-mode ultrasonography was performed at baseline and after 26 weeks. BA IMT, grayscale median (GSM), entropy, and gray-level difference statistic-contrast (GLDS-CON) were measured by a single reader. The 184 participants were (mean (SD)) 66.7 (8.2) years old and had an ankle–brachial index of 0.70 (0.18). Exercise training was associated with a 0.01 (0.06) mm (p = 0.025) reduction in BA IMT compared to 0.00 (0.05) mm (p = 0.807) in the control group (between-group p = 0.061). BA GSM, entropy, and GLDS-CON did not change significantly with exercise. Improvements in the 6-minute walk distance correlated with increases in resting BA blood flow (r = 0.23, p = 0.032), flow-mediated dilation (r = 0.24, p = 0.022), diameter (r = 0.29, p = 0.005), entropy (r = 0.21, p = 0.047), and GLDS-CON (r = 0.22, p = 0.041). In a post hoc analysis, BA IMT improved significantly with treadmill exercise training but did not change with attention control; however, the between-group difference did not reach statistical significance. With exercise, improvements in the 6-minute walk distance were associated with improved endothelial function, increased resting blood flow, and BA dilation, as well as higher grayscale entropy and GLDS-CON, indicating that lower extremity exercise is associated with salutary changes in upper-extremity arterial wall structure and function.
Keywords
Introduction
Supervised treadmill exercise (STE) improves walking performance in people with peripheral artery disease (PAD); however, its effects on peripheral arterial endothelial function in patients with PAD have been mixed.1,2 Bone marrow and splenic-derived endothelial progenitor cells released into the circulation during or after STE may affect peripheral arterial structure and function.3,4 In the PROPEL randomized clinical trial (NCT01408901), individuals with PAD were randomly assigned in a 2 × 2 factorial fashion to STE or attention control and to granulocyte-macrophage colony stimulating factor (GM-CSF) or placebo.2,5 STE but not GM-CSF improved 6-minute walk distance compared with attention control. 2 In PROPEL, brachial artery (BA) flow-mediated vasodilation (FMD; a measure of endothelial function) and peak hyperemic flow velocity (a measure of microvascular function) did not improve with treadmill exercise or GM-CSF. 2
Advances in imaging technology and image analysis software permit measurement of BA intima–media thickness (IMT), a measure of subclinical arterial injury associated with increased cardiovascular disease (CVD) risk,6,7 as well as novel grayscale ultrasound measures that describe grayscale pixel brightness, variation, and density, which may reflect tissue structure and composition. Parameters such as grayscale median (GSM) and entropy are derived from the image pixel brightness histogram and describe the echogenicity and randomness of pixel brightness in the arterial wall, respectively.8–11 Arterial contrast, calculated using gray-level difference statistics (GLDS-CON), describes spatial relationships between two ultrasound pixels.8–11 Grayscale parameters of lower GSM (darker or echolucent walls) and lower contrast (walls with more pixels that have the same grayscale value) have been associated with atherosclerotic cardiovascular disease (ASCVD) risk factors and future ASCVD events when measured in the carotid arterial bed.12–17 In the BA, GSM has been associated with CVD risk factors in a single cross-sectional study of older adults, 6 but associations of BA entropy and GLDS-CON have not been described previously. Changes in carotid artery IMT and GSM have been described with statin therapy;18–22 however, BA IMT, GSM, and other novel grayscale texture measures have not been described previously in people with PAD or with CVD risk-reducing therapies such as exercise.
We performed an exploratory analysis to evaluate the biological effects of treadmill exercise training on BA IMT as well as three novel BA grayscale ultrasound measures that may indicate subclinical arterial injury: GSM, entropy, and GLDS-CON. We hypothesized that BA IMT would improve with a treadmill exercise program and that improvements in walking performance would be related to changes in BA diameter, GSM, entropy, and GLDS-CON, indicating structural effects on the BA wall related to walking. We also sought to describe baseline associations of these novel markers with PAD severity and CVD risk factors.
Methods
Participants
This was an exploratory analysis of data from the PROPEL clinical trial.2,5 PAD was an inclusion criterion and was defined as: (i) an ankle–brachial index (ABI) ⩽ 0.90; (ii) an ABI > 0.90 if PAD was confirmed by a lower extremity angiogram or by a hospital affiliated vascular laboratory; (iii) an ABI between 0.90 and 1.00 with a ⩾ 20% decline following a heel-rise test; or (iv) a normal ABI with prior history of lower extremity revascularization but a ⩾ 20% decline following a heel-rise test.2,5,23 The exclusion criteria have been reported and were focused on safe implementation of the exercise intervention and excluded participants that already were exercising at a level similar to the intervention. 2 Randomization was stratified by the presence of diabetes mellitus. The STE intervention (walking program), GM-CSF regimen, and active control regimens, as well as the primary study’s findings, have been described previously.2,5 The CONSORT diagram (Figure 1) starts with participants from the original PROPEL results publication. 2 Baseline data were available for 186 individuals; grayscale images were measurable on 184 baseline and 182 week 26 scans.

CONSORT diagram.
Interventions
Treadmill exercise was performed three times weekly and was supervised by an exercise physiologist. Walking exercise duration was increased gradually up to 50 minutes of exercise per session. Participants were asked to exercise to maximal ischemic leg symptoms. Participants randomized to attention control attended weekly 1-hour educational sessions. GM-CSF or placebo were administered subcutaneously three times weekly for 2 weeks in a double-blinded fashion. The effect of the intervention was measured by changes in the 6-minute walk distance 24 and the maximal treadmill walking time between baseline and week 26 following standardized protocols.1,24–27
Brachial Artery (BA) ultrasound imaging
All studies were performed by a single registered diagnostic cardiac sonographer in participants who held vasoactive medications and were instructed not to smoke cigarettes or exercise prior to testing.1,2 After a 15-minute supine rest, left arm brachial blood pressure was measured by oscillometric sphygmomanometry (OMRON Model HEM-907XL; Omron Healthcare Inc., Kyoto, Japan). The right BA was imaged using a linear array vascular ultrasound transducer (ACUSON Sequoia Model C512, 8L5c transducer; Siemens Medical Solutions, Malvern, PA, USA). Extra-vascular landmarks were labeled to ensure consistent imaging within and between studies over time. Overall gain, time-gain compensation settings, and contrast were the same across studies for each participant. After baseline B-mode and Doppler images were obtained, a blood pressure cuff was placed proximal to the visualized BA segment, was inflated for 4 minutes at 50 mmHg above systolic pressure, and adjusted upwards if necessary to ensure complete occlusion using a Hokanson E20 Rapid Cuff Inflation System (D.E. Hokanson Inc., Bellevue, WA, USA). After cuff release, longitudinal images of the BA and Doppler blood flow were obtained up to 90 seconds after cuff deflation.
Brachial artery image analysis for flow and diameter measurements
Ultrasound instrumentation settings for data acquisition are provided in Table 1. Digital Imaging and Communications in Medicine images were sent electronically to the core lab for quality assessment and interpretation by a single, experienced technician using AccessPoint Web software (Freeland Systems, Alpharetta, GA, USA). BA diameters and spectral Doppler flow velocities were measured in triplicate before and after reactive hyperemia at both time points. FMD was defined as the ratio of the maximum BA diameter obtained within 90 seconds after reactive hyperemia to the resting diameter, expressed as a percent.
Ultrasound image instrumentation settings.
All images were acquired with a Siemens Medical Solutions ACUSON Sequoia with an 8L5c transducer.
Brachial artery intima-media thickness and grayscale analysis
DICOM images of the proximal BA at end-diastole while participants were resting (prior to induction of hyperemia) were converted to bitmaps for IMT measurement and grayscale analysis with LifeQ Medical software (Nicosia, Cyprus) (Figures 2–5). 28 Images were normalized by assigning the blackest area of the blood a grayscale value of 0 and the whitest area of the middle two-fourths of the adventitia a grayscale value of 190, then were standardized to a uniform pixel density of 20 mm.11,17,28,29 A reproducible 10 mm segment of the BA far wall that was present in both scans within and between visits was segmented by placing the caliper on the blood–intima and media–adventitia interfaces. From this segmentation process, the area (mm2) and grayscale texture features were measured. BA IMT was calculated by dividing the area by 10 mm (Figure 2). The first order grayscale texture features extracted for this study were the GSM (a measure of overall echogenicity) and entropy (a measure of randomness of the pixels in the far wall of the BA) based on measures from the image histogram (Figures 3–5).17,28 The texture feature contrast was extracted using the GLDS method, which describes differences in grayscale values between pixels at different distances and directions (Figures 3–5).11,28,30

Measurement technique for grayscale analysis. (A) Brachial artery diameter measurement at baseline using AccessPoint software. (B) Baseline grayscale image bitmap selected for grayscale analysis. The white rectangles demonstrate the distance of 1.0 cm on the vertical scale, and the 1.0 cm segment of the far wall of the brachial artery to the segment for grayscale analysis. Note the same segment was used for measurement of diameter and grayscale analysis across all visits. (C) The cropped grayscale and colorized segmented brachial arterial wall.

Cropped images from the grayscale LifeQ Medical software features extraction mode. (A) Cropped grayscale and colorized arterial wall. (B) Histogram measures and image histogram data for this brachial artery segment. The area for this tracing is 4.63 mm2, the IMT is 0.463 mm. (C) The grayscale median value is 29.25 and the entropy value is 4.30. The gray-level difference statistic method contrast measure is 127.977 units.

Three examples of brachial artery wall appearances from three participants. (A) The enlarged arterial wall segment demonstrates a wall with a lower grayscale median value. The image histogram would demonstrate more pixels with lower grayscale values. (B) This brachial artery wall segment demonstrates the highest entropy and contrast values of the examples. Note the variation in grayscale values for the colors depicted in the enlarged image of the segmented brachial artery wall. (C) This segmented brachial arterial wall has the highest grayscale median value (predominantly orange and red colorization) and the lowest contrast and entropy of the examples. This is a wall that is more uniform in grayscale values and has less variation in grayscale values.

Representative brachial artery B-mode images (upper frames) with extracted segment images and values below: (A) baseline, week 0; (B) after 26 weeks of treadmill exercise intervention. The lower frames demonstrate extracted and traced segments that were used for measurement of IMT and grayscale analysis with colorization below.
Measurement reproducibility
For BA IMT and grayscale measures, a single reader (JB) measured images from 24 participants twice, blinded to the first measurement. Intra-class correlation coefficients for IMT = 0.88, GSM = 0.96, GLDS-CON = 0.87, and entropy = 0.90 were excellent. The intra-class correlation coefficient for BA diameter also was excellent (0.99) as measured blindly from a subset of 10 participants by a different single reader (JMW).
Ankle–brachial index and walking outcomes
These have been described previously in the main PROPEL outcomes paper. 2 For ABI testing, a handheld Doppler probe (Nicolet Vascular Pocket-Dop II; Golden, CO, USA) was used to measure systolic blood pressures after the participant had rested supine for 5 minutes. Pressures were measured in the right brachial, dorsalis pedis, and posterior tibial arteries and left dorsalis pedis, posterior tibial, and brachial arteries and repeated. The ABI was calculated by dividing the average pressures in each leg by the average of the four brachial pressures.2,31 For the 6-minute walk test, participants walked up and down a 100-foot (30.5 meters) hallway for 6 minutes after instructions to cover as much distance as possible. The distance completed after 6 minutes was recorded.2,25 Maximal treadmill walking distance was assessed by the Gardner–Skinner protocol.2,27
Statistical analysis
All analyses were performed using XLSTAT 2017: Data Analysis and Statistical Solution for Microsoft Excel (Addinsoft, Paris, France) and SigmaPlot 13.0 (Systat Software, San Jose, CA, USA). Baseline participant characteristics were described as means (SD) overall and by treatment group, focusing on the exercise or control intervention. Differences were tested with t-tests for continuous measures or chi-squared tests for binary and categorical measures with significance set at p < 0.05. Pair-wise Pearson correlation coefficients were calculated between baseline measures and between their 26-week changes. For the primary outcome of changes in BA IMT and the three grayscale measures (GSM, entropy, and GLDS-CON) after 26 weeks, we compared within-group changes using paired t-tests and between-group changes using linear regression by group. Similar analyses were performed for GM-CSF and placebo. Baseline analyses required measureable BA parameters at baseline. The complete case analysis was performed by deleting observations with missing outcomes at either baseline or follow-up visit.
Results
Participant characteristics
Participants in this study were (mean (SD)) 66.7 (8.2) years old, 61% were males, 67% were African-American, and they had an ABI of 0.70 (0.18) (Table 2). Comorbidities included diabetes mellitus in 37%, hypertension in 84%, history of myocardial infarction in 24%, and stroke history in 15%; 54% were current smokers, 33% were previous smokers, and 13% never smoked cigarettes regularly. At baseline, there were no significant differences between participants randomized to the exercise or control interventions in these parameters or in any of the BA measures (Table 2).
Baseline participant characteristics.
The p-value from t-test or chi-squared test.
FMD, flow-mediated vasodilation; GLDS-CON, gray-level difference statistic-contrast; GSM, grayscale median; IMT, intima–media thickness.
Baseline correlations of brachial artery measures
At baseline, higher values for BA IMT were associated with older age, higher systolic blood pressure, larger BA diameter, resting BA blood flow, BA GLDS-CON, and BA entropy (Table 3). Higher BA IMT was inversely correlated with BA FMD and BA GSM (Table 3). BA IMT was higher in men than in women (0.35 [0.06] vs 0.32 [0.06] mm, p = 0.002) and in participants with diabetes mellitus (0.36 [0.07] vs 0.33 [0.06] mm, p = 0.021). In addition, BA GSM increased with BA GLDS-CON but was inversely related to body mass index, BA diameter, and resting BA blood flow. BA GSM also was lower in men than in women (75.9 [16.6] vs 81.6 [15.4] units, p = 0.021). Additional significant positive correlates of BA GLDS-CON were systolic blood pressure, BA diameter, BA IMT, and BA GSM; GLDS-CON also was higher in participants with diabetes mellitus (200.6 [102.9] vs 166.5 [80.4] units, p = 0.013). Additional significant positive correlates of BA entropy were age, systolic blood pressure, BA diameter, resting BA blood flow, and BA GLDS-CON; BA entropy also was higher in men than in women (4.35 [0.23] vs 4.27 [0.25] units, p = 0.023) and in participants with diabetes mellitus (4.38 [0.23] vs 4.29 [0.24] units, p = 0.017). None of these parameters had a significant correlation with ABI, total treadmill walking time, 6-minute walk test distance, heart rate, or diastolic blood pressure.
Baseline correlations of brachial artery measures (n = 184).
Bold type represents p-value <0.05.
BA, brachial artery; FMD, flow-mediated vasodilation; GLDS-CON, gray-level difference statistic-contrast; GSM, grayscale median; IMT, intima–media thickness.
Effects of supervised treadmill exercise
As reported previously, exercise training led to significant improvements in walking performance. 2 In this exploratory analysis, exercise training was associated with a 0.01 (0.06) mm (p = 0.025) reduction in BA IMT compared to 0.00 (0.05) mm (p = 0.807) in the control group (between-group p = 0.061; Figure 6 and Table 4). BA IMT did not change significantly in the GM-CSF group (data not shown). Significant changes in BA GSM, entropy, and GLDS-CON were not identified in the exercise group or in the group that received GM-CSF (latter data not shown), but change in the 6-minute walk test in the exercise group was associated with improvements in BA FMD and increases in BA resting blood flow, diameter, entropy, and GLDS-CON (Table 5). Heart rate did not change significantly in any group. Diastolic blood pressure declined by 2.1 (9.2) mmHg (p = 0.039) in the exercise group but this was not different than seen in the control group (between-group p = 0.486). A similar reduction was observed for systolic blood pressure in the control group.

Changes in brachial artery IMT after 26 weeks by intervention group.
Within- and between-group changes from baseline after 26 weeks.
Generalized linear model p-value.
Bold type represents p-value < 0.05.
BA, brachial artery; FMD, flow-mediated vasodilation; GLDS-CON, gray-level difference statistic-contrast; GSM, grayscale median; IMT, intima–media thickness.
Correlations between changes in walking outcomes and brachial artery measures in the exercise group (n = 92).
Bold type represents p-value < 0.05.
BA, brachial artery; FMD, flow-mediated vasodilation; GLDS-CON, gray-level difference statistic-contrast; GSM, grayscale median; IMT, intima–media thickness.
Discussion
In this exploratory analysis of data from a randomized clinical trial of treadmill exercise, GM-CSF, or attention control that previously demonstrated the beneficial effects of treadmill exercise on walking performance in people with PAD, 2 we explored the effects of treadmill exercise on BA structural measures such as wall thickness (IMT), echogenicity (GSM), and novel measures of grayscale texture (entropy and GLDS-CON) and their associations with CVD risk factors. We believe these measures reflect arterial tissue structure and composition related to biological factors that influence arterial injury and healing. In the carotid arteries, these ultrasound measures appear to be related to CVD risk.8,9,11,12,14–17,29
This report contains several novel imaging findings. First, we described, for the first time, associations of novel BA grayscale measures of echogenicity and grayscale ultrasound texture in people with PAD. Several of the BA parameters were correlated with each other and with CVD risk factors. BA GSM was inversely associated with increasing body mass index and was lower in men than in women, as previously described in the BAs of older adults and carotid arteries of smokers.6,17 Both BA GLDS-CON and entropy were associated with systolic blood pressure and were higher in participants with diabetes mellitus. BA entropy was further associated with increasing age and was higher in men. These findings suggest that grayscale texture features of the brachial arterial wall are associated with CVD risk factors and may represent tissue composition of the arterial wall. Previously, low GSM and GLDS-CON measures of the carotid arterial wall have been associated with CVD risk factors13,32 and CVD events.13,15 We did not observe significant differences between groups in any of the three BA grayscale measures with prescribed walking or GM-CSF compared to placebo; however, we identified significant positive correlations between improvements in walking performance (6-minute walk distance), BA entropy, and BA contrast. In previous work using ultrasound phantoms, we demonstrated that GLDS-CON varies based on scatterer properties of the medium that the sound wave is traveling in. 11 The findings we presented above suggest that changes in the GLDS-CON and entropy measures may describe changes in arterial wall composition which alter ultrasound scatterer properties. Further studies that use histopathology are needed to evaluate these measures in actual human arterial tissue.
Second, we demonstrated that 26 weeks of a treadmill walking exercise program was associated with a reduction in BA IMT (within-group p = 0.025); however, the comparison with the placebo intervention did not cross the threshold for statistical significance (between-group p = 0.061). Previous reports suggested that BA IMT reflects CVD injury and risk. In a study of 388 participants in a Japanese screening population, higher BA IMT was associated with CVD risk factors such as increasing age, body mass index, systolic and diastolic blood pressures, glucose levels, and smoking burden, as well as carotid IMT, a strong, independent predictor of CVD events.7,20,33 We identified similar relationships between BA IMT, age, systolic blood pressure, male sex, and diabetes mellitus. Also, in a study of 1016 older adults in Sweden, BA IMT was associated with Framingham risk score and several CVD risk factors which we did not measure in our study. 6 These findings suggest that walking exercise, in addition to improving treadmill performance in people with PAD, may reduce CVD risk in as short a time as 26 weeks; however, these findings need to be validated in a larger cohort.
Third, we identified positive relationships between increases in walking distance with improvements in BA endothelial function and BA dilation, as well as a significant reduction in diastolic blood pressure among those assigned to the exercise intervention. These observations suggest that walking exercise improves peripheral arterial function by endothelial adaptations in response to exposure to increased shear stress and cyclic strain, 34 possibly via circulating factors released from adipose tissue and skeletal muscle,35,36 and/or endothelial progenitor cells released from the bone marrow and spleen.2–4 Our finding that walking exercise is associated with structural and functional improvements in an upper arm arterial bed reinforces previous findings of the systemic effects of exercise on peripheral arterial function in PAD.2–4 Importantly, changes in BA IMT were not correlated with changes in BA diameter, suggesting that a reduction in BA wall thickness was not simply due to arterial dilation (stretch) and conservation of mass.
Other than improvements in FMD that are known to indicate improved peripheral arterial health and indicate reduced CVD risk,37–40 the clinical implications of the effects of exercise on BA IMT, entropy, and contrast are not known; however, studies that have examined ultrasound features of carotid plaques and histopathology findings of carotid endarterectomy specimens suggest that ultrasound grayscale features are strongly associated with tissue composition. Carotid artery plaques with lower GSM, increased homogeneity, and black areas near the surface of the plaque have been associated with more lipid content, inflammation, and ulceration.29,41,42 Also, common carotid artery and carotid plaque echogenicity change with statin treatment,19,21,22 so it is plausible that changes in grayscale texture features seen in this study (increase in entropy and GLDS-CON) represent changes in arterial wall composition associated with walking exercise.
Study limitations
Our study had limitations. It was an exploratory, post hoc analysis of novel ultrasound measures. Several biomarkers of CVD risk, such as lipid or inflammatory markers, were not measured in the clinical trial, so our power to detect smaller effect sizes, significant relationships between more variable measurements, and to better understand the pathobiology of the observed BA changes we observed was limited. Furthermore, our primary outcome of change in BA IMT was statistically significant in the exercise group, but the p-value for the comparison with the control group did not achieve statistical significance, likely due to the small sample. Also, our study was only 26 weeks long, so a larger effect of walking and/or GM-CSF on BA measures over a longer duration is possible.
Conclusions
In a post hoc analysis of data from a randomized clinical trial, we identified several baseline correlations between CVD risk factors and novel BA ultrasound measures such as IMT, GSM, entropy, and GLDS-CON. BA IMT improved significantly with treadmill exercise training but did not change with attention control; however, the between-group difference did not reach statistical significance. With exercise, improvements in walking distance were associated with improved endothelial function, increased resting blood flow, and BA dilation, as well as higher grayscale entropy and GLDS-CON, indicating that lower extremity exercise is associated with salutary changes in upper extremity arterial wall structure and function. These findings suggest that novel grayscale ultrasound measures are promising tools to describe arterial pathophysiology and associations with cardiovascular disease risk markers and prevention interventions.
Footnotes
Declaration of conflicting interests
The authors declared the following potential conflicts of interest with respect to the research, authorship, and/or publication of this article: Dr McDermott reported research funding from Novartis and receipt of study drug for another peripheral artery disease study from Reserveage, Hershey’s, ChromaDex, and ViroMed. No other authors reported conflicts of interest.
Funding
The authors disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: The study was funded by the National Heart, Blood, and Lung Institute (R01-HL107510), intramural support was received from the National Institute on Aging, and from the Jesse Brown VA Medical Center.
