This paper is in the following e-collection/theme issue:

Robotics in Rehabilitation (38) Technology in Physiotherapy (136) Spinal Cord Injury Rehabilitation (63) Assistive Technologies (118) Emerging Technologies for Rehabilitation (270) Portable and Mobile Technologies for Rehabilitation (156) Disabilities (118) Nails and Musculoskeletal (24)

Published on in Vol 11 (2024)

Preprints (earlier versions) of this paper are available at https://preprints.jmir.org/preprint/53084, first published .
Potential Effects of an Exoskeleton-Assisted Overground Walking Program for Individuals With Spinal Cord Injury Who Uses a Wheelchair on Imaging and Serum Markers of Bone Strength: Pre-Post Study

Potential Effects of an Exoskeleton-Assisted Overground Walking Program for Individuals With Spinal Cord Injury Who Uses a Wheelchair on Imaging and Serum Markers of Bone Strength: Pre-Post Study

Potential Effects of an Exoskeleton-Assisted Overground Walking Program for Individuals With Spinal Cord Injury Who Uses a Wheelchair on Imaging and Serum Markers of Bone Strength: Pre-Post Study

Authors of this article:

Alec Bass1, 2 Author Orcid Image ;   Suzanne N Morin3 Author Orcid Image ;   Michael Guidea1, 2 Author Orcid Image ;   Jacqueline T A T Lam1, 2 Author Orcid Image ;   Antony D Karelis4 Author Orcid Image ;   Mylène Aubertin-Leheudre4, 5 Author Orcid Image ;   Dany H Gagnon1, 2 Author Orcid Image ;   Montreal Exoskeleton Walking Program (MEWP) Group6
Article Authors Cited by (1) Tweetations (3) Metrics

    Original Paper

    1School of Rehabilitation, Faculty of Medicine, Université de Montréal, Montréal, QC, Canada

    2Centre for Interdisciplinary Research in Rehabilitation of Greater Montreal (CRIR), Centre Intégré Universitaire de Santé et Services Sociaux (CIUSSS) du Centre-Sud-de-l’Île-de-Montréal, Montréal, QC, Canada

    3Department of Medicine, McGill University, Montréal, QC, Canada

    4Department of Exercise Science, Faculty of Sciences, Université du Québec à Montréal, Montréal, QC, Canada

    5Centre de Recherche de l’Institut Universitaire de Gériatrie de Montréal (CRIUGM), Montréal, QC, Canada

    6See Acknowledgements

    *these authors contributed equally

    Corresponding Author:

    Dany H Gagnon, PT, PhD

    School of Rehabilitation

    Faculty of Medicine

    Université de Montréal

    PO Box 6128, Station Centre-Ville

    Pavillon 7077 Avenue du Parc—Office U-402-32

    Montréal, QC, H3C 3J7

    Canada

    Phone: 1 514 343 6111 ext 26382

    Email: dany.gagnon.2@umontreal.ca


    Background: As many as 60% of individuals use a wheelchair long term after a spinal cord injury (SCI). This mode of locomotion leads to chronic decline in lower-extremity weight-bearing activities and contributes to the development of severe sublesional osteoporosis and high rates of fragility fracture. Overground exoskeleton-assisted walking programs provide a novel opportunity to increase lower-extremity weight bearing, with the potential to improve bone health.

    Objective: The aim of the study is to measure the potential effects of an exoskeleton-assisted walking program on lower-extremity bone strength and bone remodeling biomarkers in individuals with chronic (≥18 months) SCI who use a wheelchair.

    Methods: In total, 10 participants completed a 16-week exoskeleton-assisted walking program (34 individualized 1-hour sessions, progressing from 1 to 3 per week). Bone mineral density and bone strength markers (dual-energy x-ray absorptiometry: total body, left arm, leg, total hip, and femoral neck and peripheral quantitative computed tomography: 25% of left femur and 66% of left tibia) as well as bone remodeling biomarkers (formation=osteocalcin and resorption=C-telopeptide) were measured before and after intervention and compared using nonparametric tests. Changes were considered significant and meaningful if the following criteria were met: P<0.1, effect size ≥0.5, and relative variation >5%.

    Results: Significant and meaningful increases were observed at the femur (femoral neck bone mineral content, bone strength index, and stress-strain index) and tibia (cortical cross-sectional area and polar moment of inertia) after the intervention (all P<.10). We also noted a decrease in estimated femoral cortical thickness. However, no changes in bone remodeling biomarkers were found.

    Conclusions: These initial results suggest promising improvements in bone strength markers after a 16-week exoskeleton-assisted walking program in individuals with chronic SCI. Additional research with larger sample sizes, longer interventions (possibly of greater loading intensity), and combined modalities (eg, pharmacotherapy or functional electrical stimulation) are warranted to strengthen current evidence.

    Trial Registration: ClinicalTrials.gov NCT03989752; https://clinicaltrials.gov/ct2/show/NCT03989752

    International Registered Report Identifier (IRRID): RR2-10.2196/19251

    JMIR Rehabil Assist Technol 2024;11:e53084

    doi:10.2196/53084

    Keywords

    assistive technology; bone architecture; bone turnover; osteoporosis; rehabilitation; spinal cord injuries; SCI; spinal cord injury; assistive device; wheelchair; exoskeleton device; locomotion; bone strength; risk; fracture 


    Mechanical loading is a key factor influencing bone strength [1]. Indeed, osteocytes detect and respond to mechanical stimuli by triggering an anabolic state that stimulates bone formation and leads to adaptations in bone geometry (known as the “mechanostat principle”) [1]. Healthy bones are therefore well adapted to the habitual loads regularly encountered during daily function (ie, concept of specificity) [2]. However, after sustaining a spinal cord injury (SCI), up to 60% of individuals use a wheelchair as their primary mode of locomotion—leading to a chronic reduction in lower-extremity weight bearing and reduced mechanical loading [3]. As a result, these individuals experience an accelerated loss in lower-extremity bone mass, particularly if no mitigation strategies are implemented during the first 18 to 24 months following the SCI [4]. This complication, referred to as sublesional osteoporosis, is associated with an increased risk of fracture, notably at the distal femur and proximal tibia [5].

    Bone strength is directly related to fracture risk and can be influenced by several characteristics, such as bone mineral density and content, as well as geometry [6]. Measuring areal bone mineral density by dual-energy x-ray absorptiometry (DEXA) remains widely recommended to assess fracture risk in this population [7]. Indeed, low areal bone mineral density has been associated with increased risks of lower-extremity fractures in individuals with SCI as well as in the general population [8]. However, solely relying on areal bone mineral density to assess bone strength can be misleading since DEXA images display 2D (ie, x- and y-axis) representations of 3D structures (ie, loss of the z-axis) [9]. DEXA condenses structures by superposing images, causing “deeper” bones to artificially appear denser (ie, increased bone mineral density) and may lead to misclassifying individuals with a lower risk of fracture [9]. As such, this limits the DEXA’s capability to inform on bone geometry (eg, cross-sectional areas and cortical thickness) [9,10]. Peripheral quantitative computed tomography (pQCT) aims to overcome this limitation by assessing volumetric bone mineral density based on 3D images [11]. Moreover, pQCT can provide additional advantages by analyzing both trabecular and cortical bone compartments separately (ie, bone geometry) and enable the estimation of mechanical properties of strength (ie, resistivity to compression, bending, and torsion).

    Although imaging (DEXA and pQCT) can provide an instantaneous “snapshot” of estimated bone strength, it does not directly assess bone turnover (remodeling). Bone turnover rate can provide fundamental information as to whether bone formation or resorption is dominant at the time of measure. Indeed, serum bone biomarkers (eg, osteocalcin and C-telopeptide) may serve as a precursor indication of a positive therapeutic effect of an intervention, even before changes can be measured with DEXA or pQCT. Osteocalcin is secreted by osteoblasts, is a marker of anabolic bone activity, and has been used in previous studies with individuals with SCI [12]. C-telopeptide, which has also been studied previously in this population, is released during bone resorption and used to characterize catabolic bone activity [13]. Since vitamin D levels can impact bone metabolism, 25-hydroxyvitamin D levels should also be measured as a possible confounding factor when characterizing serum bone biomarkers [7].

    Recently, the emergence of wearable robotic exoskeletons has led to new opportunities to develop interventions that can significantly increase lower-extremity weight bearing and mobilization. Among others, a goal of such interventions is to increase bone strength and ultimately mitigate fracture risks (and associated complications) in individuals with SCI. Pilot studies have previously demonstrated that exoskeleton-assisted walking programs are feasible in this population with high rates of satisfaction (95.2%), excellent attendance (ie, 229 completed training sessions out of 234 planned training sessions, 97.9%), and relatively low dropout rates (ie, 1 dropout out of 14 individuals recruited, 7.1%) [14,15]. In terms of learnability and ease of use, most individuals can stand and walk with walking aids and minimal assistance from a therapist by the end of the program (18 to 24 sessions) [15,16]. Walking parameters, including speed and distance, have also been shown to progress consistently and safely over the course of a walking program, especially when individualized progression strategies are used [13,15-19]. Increased walking speed and distance may provide a progressive stimulus for bone strength adaptations, equating to increased intensity and volume for these tissues. Body composition improvements have also been documented following exoskeleton-assisted walking programs, including a decrease in total and regional (ie, lower extremities) body fat and an increase in muscle mass [20]. Overall, these results are encouraging; however, the effects on bone have not been comprehensively evaluated to date.

    Thus, the main objective of this paper was to measure the potential effects of a 16-week exoskeleton-assisted walking program on lower-extremity bone density and strength and serum bone turnover markers in individuals with SCI who use a wheelchair [21]. It was hypothesized that immediate positive and meaningful effects would be observed on bone mineral density, mineral content, geometry, and mechanical strength indexes in the lower extremities as well as serum markers of bone turnover (ie, increase in bone formation markers and decrease in bone resorption markers) following the intervention.


    Ethical Considerations

    Ethics approval for this study was received on March 14, 2019, from the Centre for Interdisciplinary Research in Rehabilitation of Greater Montreal ethics committee (CRIR-1338-0518). The protocol has been published previously and was registered with the US National Library of Medicine on June 7, 2019 (ClinicalTrials.gov NCT03989752) [21].

    Study Design and Participants

    This prospective pre- and postinterventional study included adults (≥18 years of age) with chronic (ie, ≥18 months) complete or incomplete SCI. To be included, individuals needed to use a wheelchair as their primary mode of locomotion, understand French or English, and reside (or be able to arrange to reside) within 75 km of the main research site. Individuals were excluded if they had neurological impairments unrelated to the SCI (eg, multiple sclerosis); had a concomitant or secondary musculoskeletal impairment limiting their ability to safely ambulate (eg, hip heterotopic ossification); had a history of fragility fracture within the past year; or had any other condition that may preclude safe lower-extremity weight bearing, walking, or exercise tolerance (eg, unstable cardiovascular or autonomic system and renal insufficiency). Individuals also had to meet criteria specific to the wearable robotic exoskeleton (Ekso GT; Ekso Bionics) used in this study, including maximum anthropometric measures and minimal lower- and upper-extremity range of motion. Inclusion and exclusion criteria are described in greater detail in the published (open access) protocol [21].

    Measurement Times and Intervention

    Due to constraints imposed by the COVID-19 pandemic (Multimedia Appendix 1), the 4 measurement times in the published protocol were not possible. Measurement times were only possible before the intervention (2 measurements) and immediately after the intervention (1 measurement). A participant’s preintervention measurements represented the average value between measurements taken before 4 weeks and immediately before initiating the intervention. Postintervention measurements were solely taken immediately following the end of the intervention (ie, within 7 days).

    Following preintervention measurements, individuals engaged in a wearable robotic exoskeleton–assisted overground walking program consisting of 34 sessions (60 minutes per session) over a 16-week period. A published algorithm was used to individualize training volume and progression based on osteoporotic profile determined by DEXA [19]. Individuals were classified in 1 of 3 profiles: osteoporosis, osteopenia, or preserved bone mineral density. The number of steps taken per training session was then modulated, starting at 300, 400, and 500, and progressed weekly by 10%, 15%, and 20%, respectively, according to the assigned profile. For all profiles, individuals began with 1 training session per week and progressed to 3 training sessions per week by the end of the program. To maintain a moderate to vigorous exercise intensity during the sessions, walking speed, resting time, assistive devices (ie, walker or crutches), and assistance provided by the therapist were modulated to ensure a rate of perceived exertion of ≥3/10. All training sessions were supervised by a certified physiotherapist, with the help of a second physiotherapist or a physiotherapy technician if necessary.

    The exoskeleton-assisted walking program was performed using the Ekso GT exoskeleton. This ready-to-wear exoskeleton has motorized hip and knee joints and semirigid ankle orthoses. Several sensors integrated into the exoskeleton (accelerometers, gyroscopes, pressure sensors, etc) are used to detect weight transfers and movements. Front and lateral spatial targets are used to guide weight transfer with an audible sound emitted when targets are reached. Step initiation depends on the walking mode used. In “FirstStep” mode, front and lateral spatial targets must be reached, followed by the press of a confirmation button by the therapist for stepping movements to be initiated. In “ProStep” mode, stepping is automatically initiated once front and lateral spatial targets are reached (no confirmation button is pressed). In “ProStep+” mode, the lateral spatial target must be reached (no front target is necessary), and the participant must initiate a hip flexion moment to activate stepping. Additionally, the exoskeleton also provides different levels of assistance, from partial (the participant must generate some lower extremity force, and the exoskeleton assists as required) to maximal (the participant does not generate lower extremity force, and the exoskeleton realizes all movements).

    Outcomes

    DEXA Measurement

    Total body, lumbar, and left hip mineral density and content were measured using DEXA (General Electric Lunar Prodigy; standard mode; version 12.30.008). Calibration was executed daily with a standard phantom prior to each test. Participants were asked to fast for at least 8 hours prior to the assessment. Participants were also asked to empty their bladder if they had not done so within the hour preceding the DEXA. Scans were taken following the standardized protocol recommended by the manufacturer. For all scans, participants lay supine, free of jewelry or any other metallic objects. Clothing worn was noted, and participants were asked to wear the same clothing for repeated scans. For lumbar scans, participants’ lower extremities rested on a block to maintain a flexed-hip position and reduce lumbar lordosis, as recommended by the Centers for Disease Control and Prevention [22]. For hip scans, a triangular bracing device attached to the feet maintained the lower extremity in slight internal rotation, as recommended by the Centers for Disease Control and Prevention [22]. Quantitative analysis was provided automatically by the manufacturer’s software. Total body, L4 lumbar vertebrae, left arm, left leg, left total hip, and left femoral neck bone mineral densities and contents were selected as outcomes of interest. Total body measurements provided an estimate of the whole skeletal system. Lumbar vertebrae and left arm measurements provided comparators for lower extremity measurements, as changes were not expected to occur at these sites. Left leg measurements provided an estimate of the overall response of the lower extremities, which complemented the more specific pQCT measurements (described hereafter). Total hip and femoral neck sites provided a comparator with the broader osteoporosis literature, as these remain standard measurements for all populations with osteoporosis. When applicable, the left side of the body was selected to match with the pQCT scan sites.

    pQCT Measurement

    All pQCT imaging was realized on the left distal femur and proximal tibia. A standardized scan protocol was developed based on previous recommendations [11]. Calibration was executed daily with a standard phantom prior to each test. For all scans, a voxel size of 0.5×0.5 mm was used, and the scan speed was set to 10 mm/s to optimize resolution for bone and soft tissues. The total length was measured manually for the femur from the lateral femoral condyle to the greater trochanter [11]. To ensure location consistency for repeated scans, scout scans were realized at the knee joint with a reference line placed at the distal limit of the lateral femoral condyle. Following the scout scan, the pQCT was programmed to take one 2-mm slice at 25% of the total bone length calculated from the reference line. For the tibia, the total length was measured manually from the medial malleolus to the medial plateau [11]. To ensure location consistency for repeated scans, scout scans were realized at the knee joint with a reference line placed at the most distal and flattest portion of the tibial plateau. Following the scout scan, the pQCT was programmed to take one 2-mm slice at 66% of the total length calculated from the distal limit of the bone (using the reference line in this study, this equates to 33% from the knee joint). Both sites were selected to optimize for the presence of both bone and soft tissues in the scans.

    Prior to quantitative analysis, the quality of all pQCT images was independently assessed by 2 evaluators (AB and MG or JTATL) using a previously published 5-level visual inspection and quality scale, where an image score of 1 indicated high quality and an image score of 5 represented low quality [23]. To further standardize the assessment of image quality, the following criteria were agreed upon between evaluators: score 1, if the image was free of movement artifacts; score 2, if the image was only a few movement artifacts; score 3, if the image had several movement artifacts, but periosteum continuity was not affected; score 4, if the image had several movement artifacts, and periosteum continuity was affected; and score 5, if the image had movement artifacts leading to complete loss of bone continuity. A mean score was calculated for each image. Scans with a mean score greater than 3 were excluded, as such quality of the image has been proposed to be incompatible with quantitative analysis software [23]. Excluded images were treated as missing data, and measurements were computed following an intention-to-treat protocol.

    Quantitative analysis of pQCT scans was realized using the manufacturer’s software (Stratec XCT-3000; version 6.20). For all scans, contour mode 3 with a threshold set to 130 mg/cm3, peel mode 2 set to 400 mg/cm3, and separation mode 4 with an outer threshold of 200 mg/cm3 and an inner threshold of 650 mg/cm3 were used [11]. Outcomes of interest were those related to bone mineral density (total, trabecular, and cortical), bone mineral content (total, trabecular, and cortical), bone geometry (cross-sectional areas and cortical thickness), and mechanical strength indexes (bone strength index, stress-strain index, and polar moment of inertia) [7,11].

    The software provides 2 measurements for cortical thickness. The first (CRT_THK), referred hereafter as measured cortical thickness, is the mean cortical thickness based on an iterative algorithm that attempts to draw the endosteal and periosteal borders by consecutively comparing neighboring voxels (pixels). Due to occasional failure of the algorithm, particularly in individuals with severe cortical thinning and loss of cortical bone mineral density (ie, many individuals with chronic SCI), the software also provides a second measurement. This measurement (CRT_THK_C), referred hereafter as estimated cortical thickness, is based on a subtraction of endosteal radius from periosteal radius in a theoretical circular model, where total and trabecular cross-sectional areas match those measured. Since measured cortical thickness systematically failed in 2 participants, estimated cortical thickness is also reported in this study.

    Estimations of mechanical strength indexes are based on material properties and are calculated as follows. The bone strength index is the product of total bone mineral density squared by total cross-sectional area (ie, bone strength index = total bone mineral density2 × total cross-sectional area) and is indicative of resistance to compression [10,24]. The stress-strain index (resistivity to bending) is based on the calculation of the cross-sectional moment of inertia (ie, area moment of inertia or second moment of area) [10,24]. The cross-sectional moment of inertia considers the distance of cortical bone from the central axis of the bone. The greater the distance separating cortical bone from the central axis, the greater the resistivity. To calculate the stress-strain index, section modulus (Z) is computed from the cross-sectional moment of inertia in the transversal plane. Section modulus is then weighted against measured cortical bone mineral density. Thus, resistance to bending is influenced by cortical size, shape, and mineral density [10,24]. Polar moment of inertia is based on the calculation of the cross-sectional moment of inertia in the longitudinal plane [10,24]. Thus, resistance to torsion is influenced by cortical size and shape but not mineral density [10,24]. The pQCT-related variables of interest and their cross-relationships are summarized in Figure 1.

    ‎
    Figure 1. Summary of, and relationships between, outcomes of interest for peripheral quantitative computed tomography.
    Blood Samples

    Blood samples were drawn in the morning, following an 8-hour fast, by a licensed nurse into gold-top serum separator and lavender-top anticoagulant ethylenediaminetetraacetic acid tubes. Samples were immediately placed on ice and centrifuged within an hour. Serum (from gold-top serum separator tubes) and plasma (from lavender-top anticoagulant ethylenediaminetetraacetic acid tubes) were collected and stored at –80 °C until analysis. Blood samples were transported on dry ice to a university hospital laboratory at the McGill University Health Centre for analysis after the completion of the study. Serum was used to measure 25-hydroxyvitamin D, and plasma was used to measure osteocalcin and C-telopeptide.

    Statistics

    Descriptive statistics were used to characterize participants. Since the sample size was limited and some outcome measures were not normally distributed, nonparametric tests (ie, Wilcoxon signed rank test) were used to compare pre- versus postintervention data. Standardized effect sizes (r) were calculated by dividing the z value by the square root of the number of observations and interpreted as being negligible (<0.1), small (≥0.1), medium (≥0.3), or large (≥0.5) [25]. Relative pre- versus postintervention median variations (%) were also computed for all outcomes. Given the explorative nature of this study, three criteria needed to be met to reach significance and meaningfulness: (1) the α for statistical tests needed to be <.10 to balance the risk of false negatives due to an anticipated lack of statistical power, (2) calculated effect sizes needed to be large (ie, ≥0.5) for an outcome to be deemed potentially clinically relevant, and (3) relative variation needed to be greater than 5% to be considered as a change exceeding natural variability and potential measurement errors. This threshold has been used in previous work, as the least significant change reportedly varies between 2% and 5% for DEXA and pQCT depending on the location of the scan [12,26]. All statistical analyses were conducted using SPSS (version 28; IBM Corp).


    Overview

    Characteristics of the participants are summarized in Table 1. Among the 10 participants, only 1 had a very minimal motor function in the lower extremities (lower-extremity motor score: 5 out of 50), although it was not sufficient for active participation of the lower extremities during the exoskeleton-assisted walking program. Therefore, the exoskeleton was programmed to detect body weight shifts and realize stepping movements without active participation of the lower extremities (“ProStep” mode with maximal assistance in the exoskeleton) for all participants.

    Table 1. Description of the participants (N=10).
    Participant IDSexAge (y)BMD profileaWalking program progressionNeurological lesion levelAISbLEMScExoskeleton mode (Ekso GT)SCId duration (y)Weight (kg)Height (m)BMI (kg/m2)Total body fat (%)e
    1Male41PreservedFastT8A0ProStep9.666.71.7122.8f34.1f
    2Male36PreservedFastT6A0ProStep11.699.71.9227.0f39.5f
    3Male67PreservedFastT10A0ProStep12.092.31.8826.1f37.8f
    4Male60PreservedFastT11A0ProStep3.390.61.7429.9f38.7f
    5Female35PreservedFastC3C0ProStep3.650.21.6518.429
    6Male32OsteopeniaModerateT3A0ProStep8.673.51.7524.0f24.6f
    7Female48OsteopeniaModerateT12B5ProStep45.562.41.6024.4f51.8f
    8Female42OsteopeniaModerateT3A0ProStep7.770.71.6625.7f44.4f
    9Female55OsteoporosisSlowT4A0ProStep7.861.21.6622.2f43f
    10Male47OsteoporosisSlowC5A0ProStep18.381.31.8623.5f42.7f
    Mean (SD)N/Ag46.3 (10.9)N/AN/AN/AN/AN/AN/A12.8 (11.6)74.9 (15.0)1.70 (0.10)24.4 (2.9)38.5 (7.4)

    aBMD profile: preintervention bone mineral density profile of the left hip as measured by dual-energy x-ray absorptiometry (DEXA).

    bAIS: American Spinal Injury Association Impairment Scale.

    cLEMS: lower-extremity motor score on the AIS.

    dSCI: spinal cord injury.

    eTotal body fat percentage as measured by DEXA.

    fIdentifies obesity using criteria recommended by Paralyzed Veterans of America (BMI≥22 kg/m2 or body fat>22% in men and >35% in women) [27].

    gN/A: not applicable.

    DEXA Outcome Measures

    Outcome measures for DEXA are summarized in Table 2. Only the left femoral neck bone mineral content met all 3 criteria with a P=.08, a large effect size (0.55), and a relative increase of 6% postintervention.

    Table 2. Summary of dual-energy x-ray absorptiometry outcome measures (N=10).
    OutcomesPreintervention, median (IQR)Postintervention, median (IQR)P valueEffect sizea∆b (%)
    Areal bone mineral densities (g/cm2)

    		Total body bone mineral density1.159 (1.060-1.277)1.145 (1.082-1.267).800.08 (N)–1.2

    		Left arm bone mineral density1.046 (0.909-1.155)1.073 (0.889-1.221).510.20 (S)+2.6

    		Left leg bone mineral density1.018 (0.613-0.898)0.979 (0.442-0.902).450.24 (S)–3.8

    		Left total hip bone mineral density0.862 (0.756-0.992)0.832 (0.755-0.989).680.13 (S)–3.4

    		Left femoral neck bone mineral density0.852 (0.765-0.992)0.908 (0.770-0.947).110.50 (L)+6.6
    Bone mineral contents (g/cm)

    		Total body bone mineral content2759 (2377-3499)2757 (2365-3466).330.31 (M)–0.1

    		Left arm bone mineral content188 (174-236)202 (173-241).650.15 (S)+7.3

    		Left leg bone mineral content393 (300-510)370 (312-528).800.08 (N)–5.9

    		Left total hip bone mineral content28.3 (20.8-34.9)32.1 (20.2-36.7).390.27 (S)+13.5

    		Left femoral neck bone mineral contentc4.5 (3.5-6.0)4.8 (3.6-5.9).08d0.55 (L)+6

    aStandardized effect sizes interpreted as N=negligible (<0.1), S=small (≥0.1), M=medium (≥0.3), or L=large (≥0.5).

    b∆=relative variation between medians (positive indicates an increase in value from pre- to postmeasurement).

    cItalics format indicates variables meeting the following 3 criteria: statistically significant difference, effect size ≥0.5, and relative median difference ≥5%.

    dStatistically significant difference (P≤.10) for Wilcoxon signed rank tests.

    pQCT Outcome Measures

    For the femur, outcome measures for pQCT are summarized in Table 3. Although 9 outcomes were statistically significant (P<.10), only 3 had large effect sizes and sufficient relative changes to be considered as intervention effects. Bone strength index (resistivity to compression; P=.09) and stress-strain index (resistivity to bending; P=.01) increased by 9.6% and 11%, respectively, whereas estimated cortical thickness (P=.01) decreased by 9.9%. Of note, scans at the femur were not possible for 1 participant (participant 10), as his weight and lack of core stability impeded his ability to safely take and maintain the crouched sitting position necessary to set up the femur into the pQCT.

    For the tibia, outcome measures for pQCT are summarized in Table 4. Although 6 outcomes were statistically significant (P<.10), only 2 had large effect sizes and sufficient relative changes to be considered potential intervention effects. Cortical cross-sectional area (P=.06) and polar moment of inertia (P=.01) increased by 7.3% and 5.1%, respectively.

    Table 3. Summary of peripheral quantitative computed tomography outcome measures at 25% of the left femur (n=9).
    OutcomesPreintervention, median (IQR)Postintervention, median (IQR)P valueEffect sizea∆b (%)
    Volumetric bone mineral densities (mg/cm3=)

    		Total bone mineral density355.8 (334.2-470.5)381.6 (330.8-442.6).510.22 (S)+7.3

    		Trabecular bone mineral density87.7 (80.5-113.0)88.5 (83.6-110.0).150.22 (S)+1

    		Cortical bone mineral density905.9 (805.0-968.1)938.2 (871.5-981.6).04c0.69 (L)+3.6
    Bone mineral contents (mg/mm)

    		Total bone mineral content346 (275-434)341 (266-429).05c0.65 (L)–1.5

    		Trabecular bone mineral content46.6 (37.9-76.7)48.0 (39.1-78.4).950.02 (N)+3

    		Cortical bone mineral content275 (224-350)268 (217-343).01c0.89 (L)–2.5
    Bone geometry

    		Total cross-sectional area (mm2)822 (736-1066)805 (770-1023).140.49 (L)–2

    		Trabecular cross-sectional area (mm2)489 (418-700)472 (435-659).460.25 (S)–3.4

    		Cortical cross-sectional area (mm2)312 (233-394)305 (221-354).01c0.89 (L)–2.4

    		Measured cortical thickness (n=7; mm)4.03 (3.56-4.28)3.88 (3.31-4.23).03c0.83 (L)–3.6

    		Estimated cortical thickness (mm)d3.28 (2.89-3.44)2.95 (2.95-3.35).01c0.85 (L)–9.9
    Mechanical strength indexes

    		Compression: bone strength index (g/cm4)1.35 (1.16-1.60)1.48 (0.94-1.51).09c0.57 (L)+9.6

    		Bending: stress-strain index (mm3)2240 (2047-2589)2486 (2356-2706).01c0.89 (L)+11

    		Torsion: polar moment of inertia (mm4)48,002 (43,337-72,759)48,800 (42,470-71,304).02c0.77 (L)+1.7

    aStandardized effect sizes interpreted as N=negligible (<0.1), S=small (≥0.1), M=medium (≥0.3), or L=large (≥0.5).

    b∆=relative variation between medians (positive indicates an increase in value from pre- to postmeasurement).

    cStatistically significant difference (P≤.10) for Wilcoxon signed rank tests.

    dItalics format indicates variables meeting the following 3 criteria: statistically significant difference, effect size ≥0.5, and relative median difference ≥5%.

    Table 4. Summary of peripheral quantitative computed tomography outcome measures at 66% of the left tibia (N=10).
    OutcomesPreintervention, median (IQR)Postintervention, median (IQR)P valueEffect sizea∆b (%)
    Volumetric bone mineral densities (mg/cm3)

    		Total bone mineral density666.0 (571.1-772.6)669.2 (554.0-772.4).06c0.60 (L)+0.5

    		Trabecular bone mineral density97.3 (86.0-105.9)95.0 (81.3-109.5).140.47 (M)–2.4

    		Cortical bone mineral density984.9 (961.0-1007.9)956.4 (898.2-1004.8).07c0.56 (L)–2.9
    Bone mineral contents (mg/mm)

    		Total bone mineral content326 (288-425)333 (292-427).140.47 (M)+2.3

    		Trabecular bone mineral content20.1 (12.5-24.5)18.0 (13.1-24.4).880.05 (N)–10.1

    		Cortical bone mineral content283 (264-394)288 (270-398).09c0.53 (L)+1.9
    Bone geometry

    		Total cross-sectional area (mm2)602 (425-621)610 (423-660).06c0.60 (L)+1.4

    		Trabecular cross-sectional area (mm2)224 (124-274)217 (124-295).340.50 (L)–3

    		Cortical cross-sectional area (mm2)d294 (267-388)315 (273-420).06c0.60 (L)+7.3

    		Measured cortical thickness (n=8; mm)5.22 (4.74-5.67)5.31 (4.86-5.53).120.54 (L)+1.8

    		Estimated cortical thickness (mm)4.80 (3.96-5.48)4.70 (4.26-5.78).330.31 (M)–2.1
    Mechanical strength indexes

    		Compression: bone strength index (g/cm4)2.06 (1.67-2.85)2.03 (1.63-2.88).200.40 (M)–1.5

    		Bending: stress-strain index (mm3)1838 (1346-2294)1828 (1300-2250).580.18 (S)–0.5

    		Torsion: polar moment of inertia (mm4)35,706 (23,560-47,987)37,539 (23,638-49,806).01c0.79 (L)+5.1

    aStandardized effect sizes interpreted as N=negligible (<0.1), S=small (≥0.1), M=medium (≥0.3), or L=large (≥0.5).

    b∆=relative variation between medians (positive indicates an increase in value from pre- to postmeasurement).

    cStatistically significant difference (P≤.10) for Wilcoxon signed rank tests.

    dItalics format indicates variables meeting the following 3 criteria: statistically significant difference, effect size ≥0.5, and relative median difference ≥5%.

    Serum Bone Turnover Biomarkers

    Outcome measures for serum bone turnover biomarkers are summarized in Table 5. Only 25-hydroxyvitamin D met all 3 criteria with a P=.03, a large effect size, and a relative increase of 11.4% postintervention.

    Table 5. Summary of serum bone turnover biomarkers (N=10).
    OutcomesPreintervention, median (IQR)Postintervention, median (IQR)P valueEffect sizea∆b (%)
    Bone formation (μg/L)

    		Osteocalcin18.3 (15.6-19.4)21.0 (15.3-24.0).200.69 (L)+15.1
    Bone resorption (μg/L)

    		C-telopeptide0.3 (0.2-0.4)0.3 (0.2-0.4).170.43 (M)–13.8
    Others (nmol/L)

    		25-Hydroxyvitamin Dc74.5 (62.4-111)83.0 (66.3-129).03d0.69 (L)+11.4

    aStandardized effect sizes interpreted as N=negligible (<0.1), S=small (≥0.1), M=medium (≥0.3), or L=large (≥0.5).

    b∆=relative variation between medians (positive indicates an increase in value from pre- to postmeasurement).

    cItalics format indicates variables meeting the following 3 criteria: statistically significant difference, effect size ≥0.5, and relative median difference ≥5%.

    dStatistically significant difference (P≤.10) for Wilcoxon signed rank tests.


    Principal Findings

    Results of this preliminary study indicate that the completion of a progressive 16-week exoskeleton-assisted walking program may elicit some beneficial bone adaptations in individuals with chronic SCI who have limited-to-no motor function in their lower extremities and use a manual wheelchair as their primary mode of locomotion.

    DEXA Revealed an Increase in Left Femoral Neck Bone Mineral Content, but No Changes in Bone Mineral Densities

    Left femoral neck bone mineral content increased significantly and meaningfully following the intervention which is, to our knowledge, a novel and key finding partly supporting our hypotheses. Moreover, a similar trend (ie, P=.11) was also observed in left femoral neck bone mineral density (ie, +6.6% with a large effect size). Indeed, since bone mineral content and density are directly related (ie, bone mineral density = bone mineral content / area), it would be expected for both to change together. Directly comparing our results to the literature remains difficult due to the lack of previously published evidence. This is particularly true with regard to bone mineral content, as this outcome has not been reported in the limited available literature with regard to exoskeleton-assisted overground walking and treadmill-based interventions [12,20,28-32].

    Nevertheless, with regard to exoskeleton-assisted overground walking, a pilot study conducted in our laboratory did not reveal any significant changes in total body and total leg areal bone mineral densities, which is consistent with this study [20]. To our knowledge, only 2 other studies have reported areal bone mineral density measurements following exoskeleton-assisted overground walking. First, in a pilot study, an upward trend in areal bone mineral density was reported following 8 weeks of training (1 hour per session, 2 sessions per week). However, the authors neither specify in what body region this occurred nor present data to support this claim [28]. Second, in a pilot randomized controlled trial, including 16 participants with SCI (≥2 years) who use a wheelchair, areal bone mineral density (total hip and femoral neck) decreased in the activity-based exercise training group (60 minutes per session, 3 sessions per week for 24 weeks), whereas it remained stable in the exoskeleton-assisted walking group (60 minutes per session, 3 sessions per week for 24 weeks). It was hypothesized that exoskeleton-assisted walking may provide a sufficient stimulus to maintain areal bone mineral density but perhaps not to augment it [29]. Since this study did not include a comparison group, it remains unclear whether the areal bone mineral densities measured in our participants would have decreased further over the course of the study had they not participated in the walking program. However, all participants in this study sustained their SCI at least 3 years before initiating the study and were deemed to have reached a stable state in terms of bone mineral density. To this effect, it is now well evidenced that bone loss is greatest within the first 18 to 24 months following the lesion and tends to slow considerably thereafter [4]. Although a true steady state in bone mass may never be reached, it would be premature to state that the intervention in this study had a protective effect on areal bone mineral density [33]. Such a hypothesis would be best tested by recruiting participants who recently sustained their SCI (ie, no more than 2 years prior) and including a comparison group.

    The effects of treadmill-based walking programs have also been reported in the literature using robotic assistance (eg, Lokomat; Hocoma), functional electrical stimulation, or manual assistance [12,30,31]. To our knowledge, no study has reported bone mineral content, and no changes in areal bone mineral density have been previously found [12,30-32]. Since these programs imply the use of partial body weight support, the gravity-related mechanical effects decreased considerably in comparison to overground walking, which may impede the effectiveness of such programs. This is further highlighted by the fact that treadmill-based walking programs have also been tested in combination with pharmacotherapy (ie, teriparatide) and functional electrical stimulation, which should have optimized the potential effects on bone [12,30].

    Overall, this study suggests that exoskeleton-assisted overground walking may elicit a beneficial bone response at the hip that can be detected by DEXA. A combination of pharmacotherapy (eg, teriparatide), functional electrical stimulation, and overground walking may be needed to provide an optimal anabolic stimulus to significantly increase areal bone mineral density, and this warrants consideration for future research.

    Potential Improvements in Bone Strength as Measured by pQCT

    A few pQCT outcomes changed significantly and meaningfully following the completion of the intervention. This result supports our hypotheses in part. Four such outcomes increased, suggesting positive bone strength adaptations: femoral bone strength index (compression), femoral stress-strain index (bending), tibial cortical cross-sectional area, and tibial polar moment of inertia (torsion).

    With regard to the femur, to our knowledge, the increase in bone strength index is a novel finding [12,20,30-32,34]. However, an increase in stress-strain index has been previously reported in a case study following robotic-assisted treadmill training [34]. Yet, the amplitude of change reported in this previous case study (right femur=+2% and left femur=+0.5%) was much lower than in this study (ie, +11%), and may not have exceeded natural variability or measurement error. Nevertheless, these findings highlight the importance of including both femoral and tibial measurements with pQCT in this population. Since bone is expected to respond in areas of greatest mechanical strain, certain biomechanical concepts may help partially explain the results in this study [33]. First, although the increase in bone strength index would be expected with increased weight-bearing, the design of the exoskeleton may also contribute to greater compression forces at the femur during heel strike. Indeed, the exoskeleton used in this study uses a brace at the proximal tibia, just below the knee, to counteract the forward velocity of the lower limb (and body) during heel strike. Since the individuals in this study had very little-to-no motoricity in the lower limbs, this forward velocity could not be absorbed to the same extent by musculotendinous structures (ie, through eccentric contraction of the quadriceps) and would therefore be mainly absorbed by the skeletal (ie, femur) and ligamentous structures [35]. Second, due to the oblique orientation of the femoral diaphysis, it is possible that the forces with heel strike and unilateral stance during walking provide greater strain (ie, bending force) to the femur than the tibia, which may have also contributed to the results in this study [36]. Overall, these hypotheses warrant further investigation.

    With regard to the tibia, changes in cortical cross-sectional area and polar moment of inertia have been previously reported in 2 treadmill-based interventions [12,34]. However, the relatively small amplitudes of changes in these previous studies (ie, –1 to +1.4%) raise questions as to whether these changes can be attributed to more than natural measurement error. In fact, in one of these studies, comparisons with a control group yielded no significant difference for polar moment of inertia (cortical cross-sectional area was not reported in this study) [12]. Interestingly, we have previously hypothesized that the design of current exoskeletons may limit the automatic external rotation of the tibia on the femur (and consequently, the foot) during knee extension [37]. This may have led to increased torsion moments in the tibia, which would not occur during treadmill walking without robotic assistance (ie, knee extension in an open kinetic chain)—and could partially explain the difference in amplitude of change between studies.

    Uncertainties Remain Regarding pQCT Outcomes

    The fact that the estimated femoral cortical thickness decreased (–9.9%) in this study, which does not align with our hypotheses, could raise concerns regarding the possible negative effects of the walking program on bone strength. Indeed, cortical bone is largely believed to be the primary source of resistance and strength for long bones, such as the femur and tibia [9,10]. To our knowledge, these results have not been previously reported in the femur. In 1 treadmill-based trial, a statistically significant reduction of cortical thickness was reported in the tibia [12]. However, this reduction only occurred 8 months following the completion of the training program and was not statistically different than that of the control group [12]. Of interest, a statistically significant reduction in cortical cross-sectional was also observed in this study, which most likely is explained by natural variability or measurement error, considering the relatively small magnitude of change (–2.4%). Moreover, when compared to men without SCI, individuals with SCI show reductions in cortical cross-sectional area of approximately 34% [38]. Thus, the clinical significance of a 2.4% reduction in this parameter remains questionable. Nevertheless, reductions in cortical thickness and cross-sectional area may suggest that the analysis software assigned a larger proportion of bone as subcortical (identified in yellow in Figure 1), which could be related to changes in density (ie, increased porosity) at the endosteal border due to bone resorption. This possibility cannot be completely excluded from the results of this study, particularly when considering the small sample size and the limited statistical power. Future studies should pay special attention to the possible negative effects on cortical thickness and cross-sectional area at the femur.

    Serum Biomarkers Were Not Able to Contextualize pQCT Findings, but an Unexpected Increase in Levels of Serum Vitamin D Occurred

    Serum osteocalcin (bone formation) and C-telopeptide (bone resorption) did not change significantly between before and after the intervention. This provides further evidence with regard to the complexity of the interpretation of the pQCT findings, as it is not immediately obvious whether increased bone formation or resorption was occurring following the intervention. These results were not anticipated, as 4 months of treadmill walking combined with functional electrical stimulation has been shown to significantly increase osteocalcin (+6.4%) and reduce C-telopeptide (–7.7%) levels in individuals with chronic SCI [12]. The variations found in this study (ie, osteocalcin=+15.1% and C-telopeptide=–13.8%) present trends of similar direction and of greater amplitude when compared to those previously reported, although the statistical threshold was not reached.

    Serum vitamin D (25-hydroxyvitamin D) increased significantly and meaningfully by 11.4% during the intervention. Although higher vitamin D levels have been associated with greater levels of physical activity, this is generally attributed to increased time exposed to the sun in more active individuals [39]. In this study, all participants were educated regarding vitamin D supplementation recommendations by Osteoporosis Canada [40]. Participants who were not already taking vitamin D (4/10) were offered 1 year’s worth of oral supplementation. Only 1 participant began taking vitamin D supplementation during the 4-week period before initiating training. However, even when removing this participant, the data remained statistically significant (P=.05). A possible explanation for this finding is the fact that most training sessions were delivered during the transition from winter to summer months. It is well recognized that vitamin D levels tend to be lower during winter months in northern countries such as Canada, as individuals spend more time indoors [41]. Thus, it is possible that the timing of the study coincided with an expected increase in vitamin D levels seen in the general population during the transition from winter to summer [41]. Nevertheless, serum 25-hydroxyvitamin D levels remained within optimal ranges (ie, ≥75 nmol/L) throughout the duration of the study [42]. As such, bone turnover and metabolism are not expected to have been significantly affected. Moreover, vitamin D supplementation, on its own, has not been shown to effectively increase bone mineral density [43]. Therefore, it is not expected that the variations in bone markers in this study can be attributed to the measured changes in serum 25-hydroxyvitamin D levels.

    Limitations and Future Perspectives

    This study has limitations that warrant consideration when interpreting its results. First, the sample size was smaller than that initially planned due to numerous challenges associated with the COVID-19 pandemic. Consequently, this reduced statistical power and increased the chance of potential type 2 errors (ie, false negatives). Moreover, the relatively small sample size impeded the possibility of conducting additional subgroup analysis. For example, it was not possible to compare participants according to clinical characteristics (eg, gender, osteoporotic status, obesity status, and response to intervention). Unfortunately, this limits progress toward a more personalized approach for the proposed intervention. Second, the absence of bone mineral density–based inclusion or exclusion criteria led to the recruitment of 5 participants (50% of the sample size) with “preserved” bone mineral density. Hence, these participants were inherently less inclined to benefit from the walking program in terms of bone health. Third, this study did not have specific inclusion or exclusion criteria for concomitant bone health treatments. However, a complete list of medications was taken for each participant, and they were instructed to inform the research team if any changes in medications occurred during the project. Of note, none of the participants were receiving antiosteoporosis agents at the time of the study. Participants were also asked to maintain their physical activity levels during the duration of the study, including their regular exercise regime. Fourth, this study did not have a control group, as such, results should be interpreted with caution as it is unknown to what extent the absence of (or relatively small) changes measured would differ from natural variability in time. Finally, the intensity and duration of the intervention may have been insufficient. Bone resorption typically lasts 30 to 40 days, whereas bone formation frequently requires an additional 150 days, for a total bone turnover cycle requiring up to 6 months [10]. Therefore, it is plausible that clinically significant changes in bone strength could take up to 6 months, indicating that the 4-month measurement period in this study may not have been sufficient. For instance, interventions of 6 or more months, with stationary cycling assisted by functional electrical stimulation, have measured positive effects on bone mass, whereas shorter interventions have not [44-50]. Moreover, despite being initially planned, no follow-up assessments were authorized due to the COVID-19 pandemic, and the beneficial changes that may have emerged later in relation to the temporality of bone adaptation were not captured.

    Future research should focus on larger sample sizes, with a particular interest on individuals most likely to benefit from the intervention (ie, individuals with reduced bone mass). From a pragmatic perspective, large multicentric trials will be most likely required to have a sufficient sample size to detect a 5% change in femoral bone mineral density (pQCT) and compensate for large natural heterogeneity in this population. In fact, using the data in this study, this most likely entails the recruitment of roughly 200 participants based on Lehr equation (n=8s2/δ2). Interventions should be of sufficient volume (ie, at least 3 times per week), possibly of greater intensity, and of medium- to long-term durations (ie, at least 6 months) to ensure adequate stimulus and time for complete bone turnover cycles. Follow-up assessments, after the completion of the intervention, are also warranted to assess possible latent adaptations. The addition of a control group also remains relevant to compensate for natural variability and measurement error related to bone imaging and serum sampling. Finally, combining pharmacological interventions (eg, teriparatide) or functional electrical stimulation or both with overground exoskeleton–assisted walking may also warrant consideration.

    Conclusions

    The results from this paper confirm that a 16-week exoskeleton-assisted walking program may elicit bone adaptations. On one hand, significant and meaningful increases were documented via DEXA and pQCT at both the femur (ie, femoral neck bone mineral content, bone strength index, and stress-strain index) and tibia (ie, cortical cross-sectional area and polar moment of inertia). On the other hand, possible significant and meaningful decreases (ie, femoral cortical thickness) raise concerns. Although positive bone adaptations are emerging, it remains unclear whether completing a 16-week exoskeleton-assisted walking program increases bone strength in individuals with chronic SCI. The need for stronger evidence warrants additional research with larger sample sizes that focus on longer interventions (possibly of greater loading intensity), and combining modalities should be considered (eg, pharmacotherapy or functional electrical stimulation). To do so, national or international collaborations will most likely be required.

    Acknowledgments

    The authors of this paper would like to thank all the members of the Montreal Exoskeleton Walking Program (MEWP) group who contributed to a different extent to this project, particularly Cyril Duclos (kinesitherapist, PhD), Michelle McKerral (neuropsychologist, PhD), and Claude Vincent (occupational therapist, PhD). A special thanks is also extended to Martin Vermette and Florian Bobeuf for their assistance in coordinating the project, Manuel J Escalona and Sabrina Lespinasse for their help during the training sessions, as well as Michelle Wall and Rami Hammad for their aid and guidance of the blood analyses. This work was supported by the Fonds de Recherche du Québec—Santé (grant 252532), and the research infrastructure was funded in most part via the John R. Evans Leaders Fund of the Canada Foundation for Innovation (grant 36243). MA-L and DHG were supported by the Fonds de Recherche du Québec—Santé.

    Conflicts of Interest

    None declared.

    Multimedia Appendix 1

    Project timeline and effects of the COVID-19 pandemic.

    PNG File , 51 KB
    1. Yavropoulou MP, Yovos JG. The molecular basis of bone mechanotransduction. J Musculoskelet Neuronal Interact. 2016;16(3):221-236. [FREE Full text] [Medline]
    2. Daly RM, Ahlborg HG, Ringsberg K, Gardsell P, Sernbo I, Karlsson MK. Association between changes in habitual physical activity and changes in bone density, muscle strength, and functional performance in elderly men and women. J Am Geriatr Soc. 2008;56(12):2252-2260. [FREE Full text] [CrossRef] [Medline]
    3. van Middendorp JJ, Hosman AJF, Donders ART, Pouw MH, Ditunno JF, Curt A, et al. A clinical prediction rule for ambulation outcomes after traumatic spinal cord injury: a longitudinal cohort study. Lancet. 2011;377(9770):1004-1010. [FREE Full text] [CrossRef] [Medline]
    4. Clark JM, Findlay DM. Musculoskeletal health in the context of spinal cord injury. Curr Osteoporos Rep. 2017;15(5):433-442. [FREE Full text] [CrossRef] [Medline]
    5. Craven BC, Robertson LA, McGillivray CF, Adachi JD. Detection and treatment of sublesional osteoporosis among patients with chronic spinal cord injury: proposed paradigms. Top Spinal Cord Inj Rehabil. 2009;14(4):1-22. [FREE Full text] [CrossRef]
    6. Troy KL, Morse LR. Measurement of bone: diagnosis of SCI-induced osteoporosis and fracture risk prediction. Top Spinal Cord Inj Rehabil. 2015;21(4):267-274. [FREE Full text] [CrossRef] [Medline]
    7. Bryce T, Beck LA, Bombardier CH, Richmond MA, Gorman P, Marsolais EB, et al. Bone health and osteoporosis management in individuals with spinal cord injury. Clinical Practice Guideline Spinal Cord Medicine. 2022. URL: https://pva.org/wp-content/uploads/2022/05/CPG_Spinal-Cord-Medicine_2022_FINAL.pdf [accessed 2023-12-08]
    8. Lala D, Craven BC, Thabane L, Papaioannou A, Adachi JD, Popovic MR, et al. Exploring the determinants of fracture risk among individuals with spinal cord injury. Osteoporos Int. 2014;25(1):177-185. [FREE Full text] [CrossRef] [Medline]
    9. Choksi P, Jepsen KJ, Clines GA. The challenges of diagnosing osteoporosis and the limitations of currently available tools. Clin Diabetes Endocrinol. 2018;4(1):12. [FREE Full text] [CrossRef] [Medline]
    10. Hart NH, Nimphius S, Rantalainen T, Ireland A, Siafarikas A, Newton RU. Mechanical basis of bone strength: influence of bone material, bone structure and muscle action. J Musculoskelet Neuronal Interact. 2017;17(3):114-139. [FREE Full text] [Medline]
    11. Cervinka T, Giangregorio L, Sievanen H, Cheung AM, Craven BC. Peripheral quantitative computed tomography: review of evidence and recommendations for image acquisition, analysis, and reporting, among individuals with neurological impairment. J Clin Densitom. 2018;21(4):563-582. [FREE Full text] [CrossRef] [Medline]
    12. Craven BC, Giangregorio LM, Alavinia SM, Blencowe LA, Desai N, Hitzig SL, et al. Evaluating the efficacy of functional electrical stimulation therapy assisted walking after chronic motor incomplete spinal cord injury: effects on bone biomarkers and bone strength. J Spinal Cord Med. 2017;40(6):748-758. [FREE Full text] [CrossRef] [Medline]
    13. Xiang XN, Ding MF, Zong HY, Liu Y, Cheng H, He CQ, et al. The safety and feasibility of a new rehabilitation robotic exoskeleton for assisting individuals with lower extremity motor complete lesions following spinal cord injury (SCI): an observational study. Spinal Cord. 2020;58(7):787-794. [FREE Full text] [CrossRef] [Medline]
    14. Gagnon DH, Vermette M, Duclos C, Aubertin-Leheudre M, Ahmed S, Kairy D. Satisfaction and perceptions of long-term manual wheelchair users with a spinal cord injury upon completion of a locomotor training program with an overground robotic exoskeleton. Disabil Rehabil Assist Technol. 2019;14(2):138-145. [FREE Full text] [CrossRef] [Medline]
    15. Gagnon DH, Escalona MJ, Vermette M, Carvalho LP, Karelis AD, Duclos C, et al. Locomotor training using an overground robotic exoskeleton in long-term manual wheelchair users with a chronic spinal cord injury living in the community: lessons learned from a feasibility study in terms of recruitment, attendance, learnability, performance and safety. J Neuroeng Rehabil. 2018;15(1):12. [FREE Full text] [CrossRef] [Medline]
    16. Kozlowski AJ, Bryce TN, Dijkers MP. Time and effort required by persons with spinal cord injury to learn to use a powered exoskeleton for assisted walking. Top Spinal Cord Inj Rehabil. 2015;21(2):110-121. [FREE Full text] [CrossRef] [Medline]
    17. Benson I, Hart K, Tussler D, van Middendorp JJ. Lower-limb exoskeletons for individuals with chronic spinal cord injury: findings from a feasibility study. Clin Rehabil. 2016;30(1):73-84. [CrossRef] [Medline]
    18. Baunsgaard CB, Nissen UV, Brust AK, Frotzler A, Ribeill C, Kalke YB, et al. Gait training after spinal cord injury: safety, feasibility and gait function following 8 weeks of training with the exoskeletons from Ekso Bionics. Spinal Cord. 2018;56(2):106-116. [FREE Full text] [CrossRef] [Medline]
    19. Bass A, Aubertin-Leheudre M, Morin SN, Gagnon DH. Preliminary training volume and progression algorithm to tackle fragility fracture risk during exoskeleton-assisted overground walking in individuals with a chronic spinal cord injury. Spinal Cord Ser Cases. 2022;8(1):29. [FREE Full text] [CrossRef] [Medline]
    20. Karelis AD, Carvalho LP, Castillo MJ, Gagnon DH, Aubertin-Leheudre M. Effect on body composition and bone mineral density of walking with a robotic exoskeleton in adults with chronic spinal cord injury. J Rehabil Med. 2017;49(1):84-87. [FREE Full text] [CrossRef] [Medline]
    21. Bass A, Aubertin-Leheudre M, Vincent C, Karelis AD, Morin SN, McKerral M, et al. Effects of an overground walking program with a robotic exoskeleton on long-term manual wheelchair users with a chronic spinal cord injury: protocol for a self-controlled interventional study. JMIR Res Protoc. 2020;9(9):e19251. [FREE Full text] [CrossRef] [Medline]
    22. Dual energy X-ray absorptiometry (DXA) procedures manual. Centers for Disease Control and Prevention. National Health and Nutrition Examination Survey (NHANES). 2007. URL: https://www.cdc.gov/nchs/data/nhanes/nhanes_07_08/manual_dexa.pdf [accessed 2023-12-08]
    23. Blew RM, Lee VR, Farr JN, Schiferl DJ, Going SB. Standardizing evaluation of pQCT image quality in the presence of subject movement: qualitative versus quantitative assessment. Calcif Tissue Int. 2014;94(2):202-211. [FREE Full text] [CrossRef] [Medline]
    24. Kontulainen SA, Johnston JD, Liu D, Leung C, Oxland TR, McKay HA. Strength indices from pQCT imaging predict up to 85% of variance in bone failure properties at tibial epiphysis and diaphysis. J Musculoskelet Neuronal Interact. 2008;8(4):401-409. [FREE Full text] [Medline]
    25. Fritz CO, Morris PE, Richler JJ. Effect size estimates: current use, calculations, and interpretation. J Exp Psychol Gen. 2012;141(1):2-18. [CrossRef] [Medline]
    26. Kawalilak CE, Johnston JD, Olszynski WP, Kontulainen SA. Least significant changes and monitoring time intervals for high-resolution pQCT-derived bone outcomes in postmenopausal women. J Musculoskelet Neuronal Interact. 2015;15(2):190-196. [FREE Full text] [Medline]
    27. Beck LA, Bombardier CH, Richmond M, Gorman P, Marsolais EB, Hawryluk GWJ, et al. Identification and management of cardiometabolic risk after spinal cord injury: clinical practice guideline for health care providers. Clinical Practice Guideline Spinal Cord Medicine. 2018. URL: https://pva.org/wp-content/uploads/2021/09/cpg_cardiometabolic-risk_digital.pdf [accessed 2023-12-08]
    28. Wu CH, Mao HF, Hu JS, Wang TY, Tsai YJ, Hsu WL. The effects of gait training using powered lower limb exoskeleton robot on individuals with complete spinal cord injury. J Neuroeng Rehabil. 2018;15(1):14. [FREE Full text] [CrossRef] [Medline]
    29. Shackleton C, Evans R, West S, Derman W, Albertus Y. Robotic walking to mitigate bone mineral density decline and adverse body composition in individuals with incomplete spinal cord injury: a pilot randomized clinical trial. Am J Phys Med Rehabil. 2022;101(10):931-936. [FREE Full text] [CrossRef] [Medline]
    30. Gordon KE, Wald MJ, Schnitzer TJ. Effect of parathyroid hormone combined with gait training on bone density and bone architecture in people with chronic spinal cord injury. PM R. 2013;5(8):663-671. [FREE Full text] [CrossRef] [Medline]
    31. Giangregorio LM, Webber CE, Phillips SM, Hicks AL, Craven BC, Bugaresti JM, et al. Can body weight supported treadmill training increase bone mass and reverse muscle atrophy in individuals with chronic incomplete spinal cord injury? Appl Physiol Nutr Metab. 2006;31(3):283-291. [FREE Full text] [CrossRef] [Medline]
    32. Gibbs JC, Patsakos EM, Maltais DB, Wolfe DL, Gagnon DH, Craven BC. Rehabilitation interventions to modify endocrine-metabolic disease risk in individuals with chronic spinal cord injury living in the community (RIISC): a systematic search and review of prospective cohort and case-control studies. J Spinal Cord Med. 2023;46(1):6-25. [FREE Full text] [CrossRef] [Medline]
    33. Battaglino RA, Lazzari AA, Garshick E, Morse LR. Spinal cord injury-induced osteoporosis: pathogenesis and emerging therapies. Curr Osteoporos Rep. 2012;10(4):278-285. [FREE Full text] [CrossRef] [Medline]
    34. Coupaud S, Jack LP, Hunt KJ, Allan DB. Muscle and bone adaptations after treadmill training in incomplete spinal cord injury: a case study using peripheral quantitative computed tomography. J Musculoskelet Neuronal Interact. 2009;9(4):288-297. [FREE Full text] [Medline]
    35. Lee SJ, Hidler J. Biomechanics of overground vs. treadmill walking in healthy individuals. J Appl Physiol (1985). 2008;104(3):747-755. [FREE Full text] [CrossRef] [Medline]
    36. Al Nazer R, Lanovaz J, Kawalilak C, Johnston JD, Kontulainen S. Direct in vivo strain measurements in human bone-a systematic literature review. J Biomech. 2012;45(1):27-40. [FREE Full text] [CrossRef] [Medline]
    37. Bass A, Morin SN, Vermette M, Aubertin-Leheudre M, Gagnon DH. Incidental bilateral calcaneal fractures following overground walking with a wearable robotic exoskeleton in a wheelchair user with a chronic spinal cord injury: is zero risk possible? Osteoporos Int. 2020;31(5):1007-1011. [FREE Full text] [CrossRef] [Medline]
    38. Eser P, Frotzler A, Zehnder Y, Wick L, Knecht H, Denoth J, et al. Relationship between the duration of paralysis and bone structure: a pQCT study of spinal cord injured individuals. Bone. 2004;34(5):869-880. [FREE Full text] [CrossRef] [Medline]
    39. Fernandes MR, Barreto WDR. Association between physical activity and vitamin D: a narrative literature review. Rev Assoc Med Bras (1992). 2017;63(6):550-556. [FREE Full text] [CrossRef] [Medline]
    40. Hanley DA, Cranney A, Jones G, Whiting SJ, Leslie WD, Cole DEC, et al. Vitamin D in adult health and disease: a review and guideline statement from Osteoporosis Canada. CMAJ. 2010;182(12):E610-E618. [FREE Full text] [CrossRef] [Medline]
    41. Schwalfenberg G. Not enough vitamin D: health consequences for Canadians. Can Fam Physician. 2007;53(5):841-854. [FREE Full text] [Medline]
    42. Pearce SHS, Cheetham TD. Diagnosis and management of vitamin D deficiency. BMJ. 2010;340:b5664. [CrossRef] [Medline]
    43. Reid IR, Bolland MJ, Grey A. Effects of vitamin D supplements on bone mineral density: a systematic review and meta-analysis. Lancet. 2014;383(9912):146-155. [FREE Full text] [CrossRef] [Medline]
    44. Lai CH, Chang WHS, Chan WP, Peng CW, Shen LK, Chen JJJ, et al. Effects of functional electrical stimulation cycling exercise on bone mineral density loss in the early stages of spinal cord injury. J Rehabil Med. 2010;42(2):150-154. [FREE Full text] [CrossRef] [Medline]
    45. BeDell KK, Scremin AM, Perell KL, Kunkel CF. Effects of functional electrical stimulation-induced lower extremity cycling on bone density of spinal cord-injured patients. Am J Phys Med Rehabil. 1996;75(1):29-34. [CrossRef] [Medline]
    46. Bloomfield SA, Mysiw WJ, Jackson RD. Bone mass and endocrine adaptations to training in spinal cord injured individuals. Bone. 1996;19(1):61-68. [FREE Full text] [CrossRef] [Medline]
    47. Eser P, de Bruin ED, Telley I, Lechner HE, Knecht H, Stüssi E. Effect of electrical stimulation-induced cycling on bone mineral density in spinal cord-injured patients. Eur J Clin Invest. 2003;33(5):412-419. [FREE Full text] [CrossRef] [Medline]
    48. Chen SC, Lai CH, Chan WP, Huang MH, Tsai HW, Chen JJJ. Increases in bone mineral density after functional electrical stimulation cycling exercises in spinal cord injured patients. Disabil Rehabil. 2005;27(22):1337-1341. [CrossRef] [Medline]
    49. Frotzler A, Coupaud S, Perret C, Kakebeeke TH, Hunt KJ, de N Donaldson N, et al. High-volume FES-cycling partially reverses bone loss in people with chronic spinal cord injury. Bone. 2008;43(1):169-176. [FREE Full text] [CrossRef] [Medline]
    50. Mohr T, Podenphant J, Biering-Sorensen F, Galbo H, Thamsborg G, Kjaer M. Increased bone mineral density after prolonged electrically induced cycle training of paralyzed limbs in spinal cord injured man. Calcif Tissue Int. 1997;61(1):22-25. [FREE Full text] [CrossRef] [Medline]

    ‎
    DEXA: dual-energy x-ray absorptiometry
    pQCT: peripheral quantitative computed tomography
    SCI: spinal cord injury

    Edited by M Mulvenna; submitted 25.09.23; peer-reviewed by J Quinzaños, M Hasseler; comments to author 20.11.23; revised version received 20.11.23; accepted 20.11.23; published 01.01.24.

    Copyright

    ©Alec Bass, Suzanne N Morin, Michael Guidea, Jacqueline T A T Lam, Antony D Karelis, Mylène Aubertin-Leheudre, Dany H Gagnon, Montreal Exoskeleton Walking Program (MEWP) Group. Originally published in JMIR Rehabilitation and Assistive Technology (https://rehab.jmir.org), 01.01.2024.

    This is an open-access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work, first published in JMIR Rehabilitation and Assistive Technology, is properly cited. The complete bibliographic information, a link to the original publication on https://rehab.jmir.org/, as well as this copyright and license information must be included.