Ethics statements
Ethical approval was obtained from the National Rehabilitation Center Institutional Review Board (IRB no. NRC-2021-01-011). Informed consent was obtained from all participants involved in the study.
Study populations
Case group: National Rehabilitation Center (NRC) cohort
This study included men diagnosed with SCI who were hospitalized in a rehabilitation program at the National Rehabilitation Center (NRC) from 2008 to 2020. The study participants were adult men aged 40–69 years who had undergone DXA between 12 and 36 months after developing SCI. Among the initial 125 participants, participants with data showing absence of femoral neck BMD or artifacts/outliers (n = 5) and absence of lumbar spine BMD (n = 7) were excluded, and a total of 113 were analyzed in the NRC cohort (Fig. 1).
Control group: Korean National Health and Nutritional Examination Survey
Community-dwelling control data were collected from the Korean National Health and Nutritional Examination Survey (KNHANES). The KNHANES is a nationwide, population-based, cross-sectional survey and health examination conducted by the Korea Centers for Disease Control and Prevention 15. KNHANES surveys are conducted annually using a rolling sampling design that involves a complex, stratified, multistage, probability cluster survey of the general population of Korea.
Overall, 6631 community-dwelling men aged 40 to 69 years participated in the fifth KNHANES 2008–2011 survey. We excluded men without DXA scan data (n = 2318) and those whose height or weight was missing (n = 9). In total, 4,304 men were included in the analysis and matched 1:1 with the NRC cohort participants using propensity score matching based on age, weight, and height. Finally, data from 226 men (NRC cohort, n = 113; KNHANES cohort, n = 113) were analyzed (Fig. 1).
Measurements of BMD
NRC cohort
All BMD measurements of the lumbar spine, femoral neck, and total hip were performed at the NRC. However, the DXA machine was changed during the study period (GE Lunar, from January 1, 2008 to January 20, 2019; Hologic Inc., from January 21, 2019 to December 31, 2020). Because there are known systemic differences in BMD values among various manufacturers, we used standardization formulas to convert BMD values from the GE Lunar to Hologic Inc. machine. The conversion equations are as follows: lumbar spine, Hologic Inc. BMD = 0.918×GE Lunar BMD–0.038; femoral neck, Hologic Inc. BMD = 0.8638×GE Lunar BMD–0.039; and total hip, Hologic Inc. BMD = 0.971×GE Lunar BMD–0.037 16. The coefficients of variation for the lumbar spine, femoral neck, and total hip were all 1% for both the GE Lunar and Hologic Inc. machines.
KNHANES cohort
In the KNHANES cohort, the bone mineral content and BMD at the lumbar spine, femoral neck, and total hip were measured by trained technicians using DXA (DISCOVERY-W fan-beam densitometer; Hologic Inc.). The coefficients of variation for the measurements of the lumbar spine, total hip, and femur neck were 1.9%, 1.8%, and 1.5%, respectively. In both cohorts, the T-score was calculated based on the third National Health and Nutrition Examination Survey using data from young white women as a reference 17. The primary outcome of this study was the difference in the femoral neck BMD between the NRC and KNHANES cohorts.
Covariates
NRC cohort
Covariates of the NRC cohort were obtained from medical records and included age, height, weight, and medical history (history of hypertension, type 2 diabetes mellitus, and dyslipidemia). Smoking status was categorized as current smoker or nonsmoker. Heavy drinkers were defined as those drinking ³3 drinks per day. SCI-specific factors were as follows: duration since injury, cause of injury (traumatic or non-traumatic), American Spinal Injury Association Impairment Scale (AIS) (AIS-A, B classified as motor complete injury and AIS-C, D classified as motor incomplete injury) 18, and neurological level of injury (cervical, thoracic, or lumbar spine). The lower extremity motor score was calculated as the sum of manual muscle testing scores of both lower extremities from L2 to S1. To quantify the mobility of individuals with SCI, we calculated the sum of each item in the mobility indoors and outdoors’ scores of the Korean Spinal Cord Independence Measure III19. Lower extremity spasticity in both hip flexors and extensors, knee flexors and extensors, and ankle dorsiflexors and plantar flexors was measured using the Modified Ashworth Scale (MAS) 20. To calculate the sum of all 6 sites measured in this study, MAS 0 was scored as 0, MAS 1 as 1, MAS 1+ as 2, MAS 2 as 3, MAS 3 as 4, and MAS 4 as 5. Skeletal muscle mass and body fat percentage were estimated in a subset of participants via a bioelectrical impedance analysis (BIA) measurement using an eight-point tactile electrode multifrequency BIA device according to the manufacturer's instructions (InBody S10 device). During the test, clip-shaped electrodes were placed on the thumb and third finger of both hands and on both ankles with the participant in the supine position. Body fat percentage was calculated by dividing the amount of body fat by total body weight.
KNHANES cohort
Covariates of the KNHANES, including anthropometrics factors, medical history, smoking and alcohol histories were collected using standardized health questionnaires 21. Smoking status was categorized as current smoker or nonsmoker. Heavy drinkers were defined as those drinking ³3 drinks per day. Hypertension was defined as systolic blood pressure ≥140 mmHg, diastolic blood pressure ≥90 mmHg, or the use of antihypertensive medications. Diabetes mellitus (DM) was considered present if the fasting glucose level was ≥126 mg/dL or if the participant was taking oral hypoglycemic agents or insulin.
Statistical analysis
Data are presented as mean ± standard deviation or median [interquartile range] for continuous variables and as numerical value (%) for categorical variables. The independent two-sample t-test, Wilcoxon rank-sum test, and χ2 test were used to compare the baseline characteristics of the study populations as appropriate.
We performed propensity score matching between the NRC and KNHANES cohorts using the nearest-neighbor algorithm within a specified caliper width (0.2 of the standard deviation of the logit of the propensity score) on a 1:1 basis, without replacement. The propensity score was calculated based on age, weight, and height. Covariate balance was assessed using a standardized bias with a threshold of 20%. Differences in the BMD between the NRC and KNHANES cohorts were compared using the independent two-sample t-test and nonlinear modeling with a local polynomial curve. The independent effect of SCI on BMD changes was estimated using multivariable linear regression analysis. In the NRC cohort, linear logistic regression was used to evaluate the risk factors associated with low bone mass. In participants with BIA results, linear logistic regression analysis was performed to evaluate the effect of body composition on the bone.
Statistical analyses were performed using STATA 14.1 (Stata Corp., College Station, TX, USA). Statistical significance was set at P<0.05.