Skip to main content

A dynamic nomogram for predicting the probability of irreversible neurological dysfunction after cervical spinal cord injury: research based on clinical features and MRI data

Abstract

Background

Irreversible neurological dysfunction (IND) is an adverse event after cervical spinal cord injury (CSCI). However, there is still a shortage of objective criteria for the early prediction of neurological function. We aimed to screen independent predictors of IND and use these findings to construct a nomogram that could predict the development of neurological function in CSCI patients.

Methods

Patients with CSCI attending the Affiliated Hospital of Southwest Medical University between January 2014 and March 2021 were included in this study. We divided the patients into two groups: reversible neurological dysfunction (RND) and IND. The independent predictors of IND in CSCI patients were screened using the regularization technique to construct a nomogram, which was finally converted into an online calculator. Concordance index (C-index), calibration curves analysis and decision curve analysis (DCA) evaluated the model's discrimination, calibration, and clinical applicability. We tested the nomogram in an external validation cohort and performed internal validation using the bootstrap method.

Results

We enrolled 193 individuals with CSCI in this study, including IND (n = 75) and RND (n = 118). Six features, including age, American spinal injury association Impairment Scale (AIS) grade, signal of spinal cord (SC), maximum canal compromise (MCC), intramedullary lesion length (IMLL), and specialized institution-based rehabilitation (SIBR), were included in the model. The C-index of 0.882 from the training set and its externally validated value of 0.827 demonstrated the model's prediction accuracy. Meanwhile, the model has satisfactory actual consistency and clinical applicability, verified in the calibration curve and DCA.

Conclusion

We constructed a prediction model based on six clinical and MRI features that can be used to assess the probability of developing IND in patients with CSCI.

Peer Review reports

Introduction

Affected by factors such as economic development, the incidence of traumatic spinal cord injury has continued to rise in recent years, which has attracted increasing attention from orthopedic surgeons [1, 2]. In 2019, there were 0.9 million incident cases of SCI worldwide, with falls and road injuries leading causes [3]. Spinal cord injury (SCI) is a severe pathological state because of its characteristics of high disability and mortality, seriously harming human health and the social economy [4, 5]. In developed countries, the total direct cost of care for patients with SCI exceeds billions of dollars annually, with an average lifetime cost of between $ 500,000 and $ 2 million for a patient with SCI [6, 7]. Although modern medicine is sufficiently advanced, effective treatment options for these patients are still lacking, and their prognosis is also challenging to grasp [8].

More than half of the injuries occurred at the cervical level [9]. CSCI tends to affect a broader range of limb dysfunction than other segments. Injury in this region can lead to neurological dysfunction such as flaccid paralysis, central cord syndrome, and Brown-Sequard syndrome. Whether neurological dysfunction can be improved has been the most concerning topic. Early understanding of the prognosis of patients is of great significance in better comprehending the severity of the injury and formulating individualized treatment plans. However, neurological recovery often takes a long time to reach the plateau, with some patients even lasting for one to two years [10, 11]. Therefore, it is hard to accurately evaluate the neurological outcome in the early stage of injury.

The American Spinal Injury Association Impairment Scale (AIS) has been revised several times since it was proposed, and it has become an international standard for evaluating the degree of SCI [12]. Many studies have confirmed the correlation between AIS grade and the outcome of neurological function, and the AIS grade is often used to determine the prognosis roughly [13,14,15]. However, this evaluation index is independent of clinical characteristics, imaging data, and treatment intervention and has significant limitations. In this study, we attempted to integrate multiple factors to develop a model to evaluate the neurological prognosis of CSCI patients more accurately. We hope to provide new ideas for the clinical treatment and risk identification of neurological dysfunction after CSCI.

Patients and methods

Patient selection

The data of patients with acute CSCI hospitalized in the Affiliated Hospital of Southwest Medical University from January 2015 to March 2021 were retrospectively analyzed as a training cohort. Subsequently, the data of patients with acute CSCI hospitalized in the Affiliated Chengdu 363 Hospital of Southwest Medical University during the same period externally validated our results.

The inclusion criteria are as follows: 1) Patients had a history of trauma and were diagnosed with CSCI; 2) The imaging data was complete and consistent with symptoms and signs; 3) Spinal shock resolved. The patients below are not qualified for this study: 1) Patients with previous neurological paralysis symptoms before this trauma; 2) Patients with severe craniocerebral trauma; 3) Patients with spinal cord malformation such as syringomyelia, basilar invagination, hemivertebra deformity, split cord malformation, etc.; 4) The history before admission more than one month; 5) Patients who did not receive surgical treatment; 6) Incomplete data; 7) Patients of less than 14 years of age. The Medical Ethics Committee of the Affiliated Hospital of Southwest Medical University approved this study.

Observation indicators

The following data were documented for the study population: Age, gender, etiology, AIS grade [16], calcium, potassium, white blood cells (WBC), neutrophil to lymphocyte ratio (NLR), mean arterial pressure (MAP), signal of spinal cord (SC), maximum canal compromise (MCC), intramedullary lesion length (IMLL), operation timing, surgery duration, surgical approach, specialized institution-based rehabilitation (SIBR).

Detailed clinical information was acquired from the medical record. All laboratory indicators and MAP were obtained from data collected within 24 h of admission. We divided the etiology into traffic-accident, high falls, struck by falling, low falls, and other traumatic causes. The signal of SC was classified as normal and abnormal based on MRI findings, with abnormalities including spinal cord edema, hemorrhage, and transection. On the basis of T1-MRI results [17], the MCC was calculated by measuring changes in the anteroposterior diameter of the spinal canal at the maximum level of injury and was further divided into three grades (< 30%, 30%-50%, > 50%) (Fig. 1A). IMLL was defined as the rostrocaudal length of high-signal intensity from the injury epicenter [18], measured twice times on MRT2WI, using average data (Fig. 1B). Operation timing was converted into binary variables early (< 72 h after injury) and late (≥ 72 h after injury). The surgical approach was divided into the anterior, posterior, and posterior-anterior combined approaches. The definition of SIBR was receiving treatment under the guidance of specialized rehabilitation institutions.

Fig. 1
figure 1

Measurement methods of MCC and IMLL. Abbreviations: MCC: maximum canal compromise; IMLL: intramedullary lesion length. Notes: (A) Select the median sagittal image of the magnetic resonance imaging, di represents the anteroposterior canal diameter of the maximum damage level, while da and db represent the diameter of anteroposterior canal diameter nearest to the normal level above and below the injury site, respectively. MCC (%) = 1-[ 2di / (da + db)] × 100%. B Select the median sagittal image of the magnetic resonance imaging, and the rostrocaudal length of damage was measured in centimeters. The normal intramedullary signal was recorded as zero

In order to make the prediction more accurate, we analyzed two important factors that may affect the efficacy of SIBR, including the start time of SIBR and the length of SIBR time.

Laboratory indices were assayed using a Sysmex XE-2100 automated hematology analyzer and Siemens Advia 2400 chemistry analyzer. Imaging data were scanned with Philips Intera Achieva 1.5 T MR machine. All images are uploaded to the iMedPacs system via a synchronous workstation and measured using the system's measurement tools.

Outcome assessment

The outcome indicator of this study was the dichotomous variable of irreversible neurological dysfunction (IND) or reversible neurological dysfunction (RND) in CSCI patients. All patients were classified into IND and RND according to the following criteria: the differences in AIS grade at the baseline status and the observation time point from the twelve months. The improvement of AIS grade by at least one level was taken as the criterion for CSCI reversible neurological dysfunction. Telephone follow-up and outpatient follow-up medical records were used.

Statistical analysis

R Studio ver. 4.1.2 and SPSS Statistics ver. 26.0 were used for statistical analysis. Continuous data were described as the means ± standard deviations or percentiles, while the categorical data were expressed in percentages. The T-test, Mann–Whitney U-test, and the chi-square test were used to retrospectively analyze the patients' differences in demographic and clinical characteristics. Statistical significance was confirmed at a two-tailed P-value lower than 0.05. The model was built, evaluated, and validated through R Studio.

Independent predictors were confirmed using the least absolute shrinkage and selection operator (LASSO) method. Visualization of the model was made possible with a static nomogram, which was eventually transformed into an online calculator to promote usage efficiency. The concordance index (C-index) estimated the model's discrimination, which could be replaced with the area under the receiver operating characteristic (ROC) curve (AUC). We also exploited the calibration curve analysis and the decision curve analysis (DCA) for model assessment. The calibration curve was performed to evaluate whether the predicted IND probability corresponded with the actual proportion. DCA was adopted to extrapolate the model's clinical applicability by calculating the net benefits at different threshold probabilities. Afterward, external validation of the model was conducted in the validation cohort, and it was also internally validated using the bootstrap method.

Results

Baseline characteristics

Data of 279 patients with CSCI were reviewed, and 193 patients were included in this study after excluding non-compliant patients, including 75 IND and 118 RND (Fig. 2). The demographic data and clinical features of patients were recorded in Table 1. There were 151 males (78.2%) and 42 females (21.8%). The median age was 54 years (42.5–60 years). The results showed that there were statistically significant differences between IND and RND in AIS grade, signal of SC, MCC, IMLL, and NLR (P  < 0.05). The baseline characteristics of the validation cohort were recorded in Table 2. There were no statistically significant differences in Age, AIS grade, signal of SC, MCC, IMLL and SIBR between the data from the training and validation cohorts.

Fig. 2
figure 2

Patient flowchart. Abbreviations: IND: irreversible neurological dysfunction; RND: reversible neurological dysfunction

Table 1 Baseline characteristics of the study population
Table 2 Comparison of the variable characteristics of the validation and training cohort

In addition, we obtained detailed data on 38 patients who received SIBR. As shown in Table S1, the length of SIBR time and the start time of SIBR were not significantly different in the IND and RND groups. Pairwise comparison within the group found that receiving SIBR within one month after injury had a higher recovery rate than receiving SIBR more than three months after injury ( P = 0.032). After Bonferroni correction ( P'=0.05 / [ k ( k-1) / 2] = 0.017, P'  <  P), the difference between these two groups was not considered to be statistically significant.

Model construction

After the LASSO variables selection, we screened six independent predictors, including age, AIS grade, signal of SC, MCC, IMLL, and SIBR (Figure S1). Based on these indicators, a model by multivariate logistic regression was developed and presented in the form of a static nomogram (Fig. 3). Their odds ratio (OR) values and regression coefficients are displayed in Table 3. In addition, we simplified the model to a web calculator whose access link is https://CSCI.shinyapps.io/dynnomapp/ (Figure S2).

Fig. 3
figure 3

Nomogram for predicting IND after CSCI in the training cohort. Abbreviations: IND: irreversible neurological dysfunction; CSCI: cervical spinal cord injury; AIS: American spinal injury association Impairment Scale; SC: spinal cord; MCC: maximum canal compromise; IMLL: intramedullary lesion length; SIBR: specialized institution-based rehabilitation. Notes: The nomogram shows the six independent predictors of the model, with the abscissa being a graphical score assigned to the predictor. The sum of these six scores generated a plot on the “Total points” axis, with the corresponding values below representing the probability of IND after CSCI

Table 3 Multivariate logistic regression model in the training cohort

Model evaluation

The C-index of 0.882 (95%CI: 0.833 – 0.931) in the training cohort demonstrated a good discriminant ability. The threshold possibility was 0.321, with the optimal sensitivity of 0.763 and specificity of 0.853 (Fig. 4A). In the calibration curve, the predicted curve and the reference curve (45-degree diagonal) roughly overlap, suggesting that the predicted incidence of outcome events was highly consistent with the actual situation (Fig. 4B). DCA shows the model was clinically beneficial (Fig. 4C). The C-index of 0.891 from the internal bootstrap validation verified the stability of the model. In the validation cohort, the C-index was 0.827 (95%CI: 0.778 – 0.876) (Fig. 4D), which corresponded with our results. Meanwhile, the calibration curve (Fig. 4E) and DCA (Fig. 4F) showed that the model had satisfactory actual consistency and clinical applicability.

Fig. 4
figure 4

Validation of nomogram. Abbreviations: ROC: receiver operating characteristic; AUC: area under the receiver operating characteristic. Notes: Plots A, B, and C are from the training cohort, and plots D, E, and F are from the validation cohort. Receiver operating characteristic curve of the nomogram (A and D). Point P is the threshold probability when the model has optimal sensitivity and specificity. Calibration curves of the nomogram (B and E). The ideal line indicates that the predicted probability perfectly coincides with the actual occurrence of IND, the apparent line reflects the prediction performance of the established model, and the distance between these two lines reflects the actual consistency of the model. Decision curve analysis of the nomogram (C and F). The X-axis represents the threshold probability, and the Y-axis represents the net benefit. The model has better clinical applicability when the net benefit of the nomogram is greater than the All and None lines at the same threshold

Discussion

Acute CSCI is a severe nerve trauma, and its prognosis varies greatly. Several scholars have confirmed that the prognosis of SCI is related to many factors, such as injury specificity and post-injury management [10, 19]. However, the weight of each element on the prognosis is different. Therefore, it is unreasonable to use only one feature as a decisive factor for estimating the outcome of patients with SCI. The nomogram solves these problems well and assigns different scores to each prediction factor to visually present the probability of outcome events.

A total of 16 clinical and imaging features were analyzed in our study. After Lasso regression, the six indicators were included: age, AIS grade, signal of SC, MCC, IMLL, and SIBR. Age, AIS grade, signal of SC, MCC, and IMLL could be easily obtained directly during hospitalization. SIBR is indirectly achieved by understanding the wishes of patients. Therefore, the nomogram developed in this study can be used as an easy-to-use tool for clinical application. We confirmed the model's high accuracy, stability, and clinical applicability. In addition, we created an online calculator for clinical application. Doctors could conveniently obtain the possibility of IND in patients with CSCI through a mobile terminal.

Previous studies have suggested that age may not be related to motor recovery. However, the results obtained without variable stratification are questionable, given the heterogeneity of clinical characteristics and treatment options in CSCI patients. Furlan et al. [20] found that the potential of older patients with SCI to improve within the first year post-injury neurologically does not appear to translate into similar functional recovery compared to that seen in younger individuals (< 65 years). Some scholars also reported that compared with patients under 50 years old, the elderly were less likely to achieve walking function post-injury [21]. Furthermore, we conducted a more detailed grouping of age, and we still found that age can be screened as an independent predictor of IND after CSCI. Overall, older age appears to be associated with poorer neurological recovery, although the exact age is unclear [11]. These results may be related to reduced antioxidant defenses and axonal plasticity in elderly patients [22, 23]. In addition, negative management attitudes due to ageism and multimorbidity play a non-negligible role [24, 25].

In many studies, the severity of SCI is considered the most critical factor in predicting neurological outcomes. As a recognized standard for evaluating neurological function in SCI, AIS shows good overall discrimination in predicting the outcome of neurological function in patients. In the model established in this study, both AIS A (OR = 41.325) and AIS D (OR = 11.171) had a significant odds ratio, indicating a strong correlation between AIS as a predictor and outcome events. Previous literature has shown that the recovery rate of neurological function in SCI patients met the rule of grade C of AIS > B > D > A [10, 26, 27], which is consistent with the score of each grade in this model. Interestingly, although grade D of AIS is the lightest in injury classification, its neurological prognosis is not the best. The ceiling effect of the AIS grade may have contributed to this result.

As an essential examination for patients with CSCI, MRI reflects the state of spinal cord after injury and the degree of violence [28]. Boese et al. [29] found that the signal of SC was associated with AIS grade improvement in adult patients with SCI without radiologic abnormalities. They believed that patients with normal MRI findings had a better prognosis than those with abnormalities, while patients with extraneural abnormalities (compression) had a better prognosis than patients with intraneural abnormalities (edema, hemorrhage, contusion, and transection). Furthermore, Fehlings et al. [30] reported that SCI patients with the signal of SC changes showed worse neurological function during one follow-up. Our results coincided with the previously proposed views. We also find that the probability of IND increases with the increase of IMLL. Some scholars reported IMLL increased by 10 mm, and the AIS grade conversion rate can be reduced by 40% [31]. Matsushita et al. [32] concluded that patients were more likely to achieve walking function at discharge when the vertical diameter of the T2 high-intensity area was less than 45 mm. Another study found that incomplete SCI with lesion length greater than 36 mm was independently associated with a more severe AIS grade one year later [33]. It is worth noting that the IMLL does not always remain constant, which also gives the model the potential for dynamicity. In addition, World Federation of Neurosurgical Societies Spine Committee members and experts agree that MCC is an important prognostic indicator [34]. More significant spinal cord compression often leads to more severe neurological dysfunction [35], while MCC shows the extent of spinal cord compression in the plane of injury. In this study, patients with MCC less than 30% had the highest recovery rate at follow-up, while those with MCC greater than 50% had the worst prognosis, which is consistent with the above view. The cause of this result may be related to persistent compression of the spinal cord triggered by serious canal compromise. Furthermore, we believe that the spinal cord suffered greater instantaneous violence at serious canal compromise.

Although the contribution of rehabilitation therapy to functional improvement remains controversial, the recent clinical practice guidelines still recommend rehabilitation treatment for patients with acute SCI who are stable and able to tolerate the required rehabilitation intensity [36]. Standardized rehabilitation reduces complications and positively affects neurological recovery, even in patients with AIS A [37]. The efficacy and application potential of some rehabilitation measures, such as hyperbaric oxygen therapy, electrical stimulation and hydrotherapy, has been confirmed in previous studies [38,39,40]. In general, low-intensity rehabilitation can begin within three days of the injury, while high-intensity activities need to wait until at least four days after the injury to be relatively safe [41]. Earlier rehabilitation maximizes the important window of the first three months and offers greater potential for prognostic improvement. Herzer et al. [42] found that the prolonged time between injury and rehabilitation was associated with lower motor scores after one year. In addition, within a certain period, the length of rehabilitation time was also positively correlated with the final motor function score [37]. Our findings suggest a potential correlation between early rehabilitation interventions and better prognosis. Unlike previous studies, we found that the length of rehabilitation was not associated with prognosis, which may be attributed to injury specificity, population differences, and sample size. Although our study confirmed that patients receiving SIBR had a higher recovery rate, there was still a lack of a standardized strategy to guide rehabilitation treatment.

Several limitations of this study should be explained. First, this is a retrospective study that may have inherent information bias. Second, although the model has shown its superiority after verification by data from another center, the overall sample size of this study is small, which may not strictly meet the principle of events per variable. Third, the limited sample size of the earlier surgical (< 24 h) group due to various conditions was not included in the evaluation, which was considered an independent protective factor in previous studies. Finally, most of our cases were surgically treated, so we excluded cases without surgical records to reduce data bias. This may result in the model not being well applied to patients without surgery.

Conclusion

Our study constructed a nomogram and online calculator based on six clinical and MRI features that can be used to identify patients with CSCI who are at risk of IND. Early prediction of IND occurrence is expected to guide clinicians in providing personalized support and treatment strategies.

Availability of data and materials

The authors declare the data analyzed during the current study could be available from the corresponding author upon reasonable request.

Abbreviations

IND:

Irreversible neurological dysfunction

CSCI:

Cervical spinal cord injury

RND:

Reversible neurological dysfunction

SCI:

Spinal cord injury

AIS:

American Spinal Injury Association Impairment Scale

WBC:

White blood cells

NLR:

Neutrophil to lymphocyte ratio

MAP:

Mean arterial pressure

SC:

Spinal cord

MCC:

Maximum canal compromise

IMLL:

Intramedullary lesion length

SIBR:

Specialized institution-based rehabilitation

LASSO:

Least absolute shrinkage and selection operator

C-index:

Concordance index

ROC curve:

Receiver operating characteristic curve

AUC:

Area under the receiver operating characteristic

DCA:

Decision curve analysis

OR:

Odds ratio

References

  1. Beck B, Cameron PA, Braaf S, Nunn A, Fitzgerald MC, Judson RT, et al. Traumatic spinal cord injury in Victoria, 2007–2016. Med J Aust. 2019;210(8):360–6.

    Article  PubMed  Google Scholar 

  2. Hao D, Du J, Yan L, He B, Qi X, Yu S, et al. Trends of epidemiological characteristics of traumatic spinal cord injury in China, 2009–2018. Eur Spine J. 2021;30(10):3115–27.

    Article  PubMed  Google Scholar 

  3. Ding W, Hu S, Wang P, Kang H, Peng R, Dong Y, et al. Spinal Cord Injury: The Global Incidence, Prevalence, and Disability From the Global Burden of Disease Study 2019. Spine (Phila Pa 1976). 2022;47(21):1532–40.

    Article  PubMed  Google Scholar 

  4. Chamberlain JD, Meier S, Mader L, von Groote PM, Brinkhof MW. Mortality and longevity after a spinal cord injury: systematic review and meta-analysis. Neuroepidemiology. 2015;44(3):182–98.

    Article  PubMed  Google Scholar 

  5. Jain NB, Ayers GD, Peterson EN, Harris MB, Morse L, O’Connor KC, et al. Traumatic spinal cord injury in the United States, 1993–2012. JAMA. 2015;313(22):2236–43.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. McDaid D, Park AL, Gall A, Purcell M, Bacon M. Understanding and modelling the economic impact of spinal cord injuries in the United Kingdom. Spinal Cord. 2019;57(9):778–88.

    Article  PubMed  PubMed Central  Google Scholar 

  7. DeVivo MJ. Causes and costs of spinal cord injury in the United States. Spinal Cord. 1997;35(12):809–13.

    Article  CAS  PubMed  Google Scholar 

  8. Munteanu C, Rotariu M, Turnea M, Ionescu AM, Popescu C, Spinu A, et al. Main cations and cellular biology of traumatic spinal cord injury. Cells. 2022;11(16):2503.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Chen Y, He Y, DeVivo MJ. Changing demographics and injury profile of new traumatic spinal cord injuries in the United States, 1972–2014. Arch Phys Med Rehabil. 2016;97(10):1610–9.

    Article  PubMed  Google Scholar 

  10. Khorasanizadeh M, Yousefifard M, Eskian M, Lu Y, Chalangari M, Harrop JS, et al. Neurological recovery following traumatic spinal cord injury: a systematic review and meta-analysis. J Neurosurg Spine. 2019;15:1–17.

    Google Scholar 

  11. Kirshblum S, Snider B, Eren F, Guest J. Characterizing natural recovery after traumatic spinal cord injury. J Neurotrauma. 2021;38(9):1267–84.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Hadley MN, Walters BC, Grabb PA, Oyesiku NM, Przybylski GJ, Resnick DK, et al. Clinical assessment after acute cervical spinal cord injury. Neurosurgery. 2002;50(3 Suppl):S21–9.

    CAS  PubMed  Google Scholar 

  13. Engel-Haber E, Radomislensky I, Peleg K, Bodas M, Bondi M, Noy S, et al. Early trauma predictors of mobility in people with spinal cord injury. Spine (Phila Pa 1976). 2021;46(20):E1089-96.

    Article  PubMed  Google Scholar 

  14. Hicks KE, Zhao Y, Fallah N, Rivers CS, Noonan VK, Plashkes T, et al. A simplified clinical prediction rule for prognosticating independent walking after spinal cord injury: a prospective study from a Canadian multicenter spinal cord injury registry. Spine J. 2017;17(10):1383–92.

    Article  PubMed  Google Scholar 

  15. Engel-Haber E, Zeilig G, Haber S, Worobey L, Kirshblum S. The effect of age and injury severity on clinical prediction rules for ambulation among individuals with spinal cord injury. Spine J. 2020;20(10):1666–75.

    Article  PubMed  Google Scholar 

  16. Kirshblum SC, Burns SP, Biering-Sorensen F, Donovan W, Graves DE, Jha A, et al. International standards for neurological classification of spinal cord injury (revised 2011). J Spinal Cord Med. 2011;34(6):535–46.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Furlan JC, Kailaya-Vasan A, Aarabi B, Fehlings MG. A novel approach to quantitatively assess posttraumatic cervical spinal canal compromise and spinal cord compression: a multicenter responsiveness study. Spine (Phila Pa 1976). 2011;36(10):784–93.

    Article  PubMed  Google Scholar 

  18. Na BR, Seo HY. Adult spinal cord injury without major bone injury: effects of surgical decompression and predictors of neurological outcomes in American Spinal Injury Association Impairment Scale A, B, or C. J Clin Med. 2021;10(5):1106.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Dvorak MF, Noonan VK, Fallah N, Fisher CG, Rivers CS, Ahn H, et al. Minimizing errors in acute traumatic spinal cord injury trials by acknowledging the heterogeneity of spinal cord anatomy and injury severity: an observational Canadian cohort analysis. J Neurotrauma. 2014;31(18):1540–7.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Furlan JC, Fehlings MG. The impact of age on mortality, impairment, and disability among adults with acute traumatic spinal cord injury. J Neurotrauma. 2009;26(10):1707–17.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Oleson CV, Marino RJ, Leiby BE, Ditunno JF. Influence of age alone, and age combined with pinprick, on recovery of walking function in motor complete, sensory incomplete spinal cord injury. Arch Phys Med Rehabil. 2016;97(10):1635–41.

    Article  PubMed  Google Scholar 

  22. Jaerve A, Schiwy N, Schmitz C, Mueller HW. Differential effect of aging on axon sprouting and regenerative growth in spinal cord injury. Exp Neurol. 2011;231(2):284–94.

    Article  CAS  PubMed  Google Scholar 

  23. Stewart AN, Glaser EP, Mott CA, Bailey WM, Sullivan PG, Patel SP, et al. Advanced age and neurotrauma diminish glutathione and impair antioxidant defense after spinal cord injury. J Neurotrauma. 2022;39(15–16):1075–89.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Furlan JC, Craven BC, Fehlings MG. Is there any gender or age-related discrepancy in the waiting time for each step in the surgical management of acute traumatic cervical spinal cord injury? J Spinal Cord Med. 2019;42(sup1):233–41.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Kopp MA, Lübstorf T, Blex C, Schwab JM, Grittner U, Auhuber T, et al. Association of age with the timing of acute spine surgery-effects on neurological outcome after traumatic spinal cord injury. Eur Spine J. 2022;31(1):56–69.

    Article  PubMed  Google Scholar 

  26. Spiess MR, Müller RM, Rupp R, Schuld C, EM-SCI Study Group, van Hedel HJ. Conversion in ASIA impairment scale during the first year after traumatic spinal cord injury. J Neurotrauma. 2009;26(11):2027–36.

    Article  PubMed  Google Scholar 

  27. Nakajima H, Yokogawa N, Sasagawa T, Ando K, Segi N, Watanabe K, et al. Prognostic factors for cervical spinal cord injury without major bone injury in elderly patients. J Neurotrauma. 2022;39(9–10):658–66.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Kumar Y, Hayashi D. Role of magnetic resonance imaging in acute spinal trauma: a pictorial review. BMC Musculoskelet Disord. 2016;17:310.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Boese CK, Lechler P. Spinal cord injury without radiologic abnormalities in adults: a systematic review. J Trauma Acute Care Surg. 2013;75(2):320–30.

    Article  PubMed  Google Scholar 

  30. Jentzsch T, Cadotte DW, Wilson JR, Jiang F, Badhiwala JH, Akbar MA, et al. Spinal cord signal change on magnetic resonance imaging may predict worse clinical in- and outpatient outcomes in patients with spinal cord injury: a prospective multicenter study in 459 patients. J Clin Med. 2021;10(20):4778.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Aarabi B, Sansur CA, Ibrahimi DM, Simard JM, Hersh DS, Le E, et al. Intramedullary lesion length on postoperative magnetic resonance imaging is a strong predictor of ASIA impairment scale grade conversion following decompressive surgery in cervical spinal cord injury. Neurosurgery. 2017;80(4):610–20.

    Article  PubMed  Google Scholar 

  32. Matsushita A, Maeda T, Mori E, Yuge I, Kawano O, Ueta T, et al. Can the acute magnetic resonance imaging features reflect neurologic prognosis in patients with cervical spinal cord injury? Spine J. 2017;17(9):1319–24.

    Article  PubMed  Google Scholar 

  33. Martínez-Pérez R, Cepeda S, Paredes I, Alen JF, Lagares A. MRI prognostication factors in the setting of cervical spinal cord injury secondary to trauma. World Neurosurg. 2017;101:623–32.

    Article  PubMed  Google Scholar 

  34. Parthiban J, Zileli M, Sharif SY. Outcomes of spinal cord injury: WFNS Spine Committee Recommendations. Neurospine. 2020;17(4):809–19.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Jin C, Zhao L, Wu J, Jia L, Cheng L, Xie N. Traumatic cervical spinal cord injury: relationship of MRI findings to initial neurological impairment. Eur Spine J. 2021;30(12):3666–75.

    Article  PubMed  Google Scholar 

  36. Fehlings MG, Tetreault LA, Aarabi B, Anderson P, Arnold PM, Brodke DS, et al. A clinical practice guideline for the management of patients with acute spinal cord injury: recommendations on the type and timing of rehabilitation. Global Spine J. 2017;7:231S-S238.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Truchon C, Fallah N, Santos A, Vachon J, Noonan VK, Cheng CL. Impact of therapy on recovery during rehabilitation in patients with traumatic spinal cord injury. J Neurotrauma. 2017;34(20):2901–9.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Sun ZC, Liang F, Yang J, Hai Y, Su QJ, Liu XH. The mechanism by which hyperbaric oxygen treatment alleviates spinal cord injury: genome-wide transcriptome analysis. Neural Regen Res. 2022;17(12):2737–42.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Arpin DJ, Ugiliweneza B, Forrest G, Harkema SJ, Rejc E. Optimizing neuromuscular electrical stimulation pulse width and amplitude to promote central activation in individuals with severe spinal cord injury. Front Physiol. 2019;10:1310.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Wiesener C, Spieker L, Axelgaard J, Horton R, Niedeggen A, Wenger N, et al. Supporting front crawl swimming in paraplegics using electrical stimulation: a feasibility study. J Neuroeng Rehabil. 2020;17(1):51.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Lewis NE, Tabarestani TQ, Cellini BR, Zhang N, Marrotte EJ, Wang H, et al. Effect of acute physical interventions on pathophysiology and recovery after spinal cord injury: a comprehensive review of the literature. Neurospine. 2022;19(3):671–86.

    Article  PubMed  PubMed Central  Google Scholar 

  42. Herzer KR, Chen Y, Heinemann AW, González-Fernández M. Association between time to rehabilitation and outcomes after traumatic spinal cord injury. Arch Phys Med Rehabil. 2016;97(10):1620–7.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We wish to thank all the study participants and research staff who participated in this research.

Funding

This study is supported by the Youth Innovation Project of Sichuan Medical (Q19038) and the Doctoral research startup fund of Affiliated Hospital of Southwest Medical University (Q19081).

Author information

Authors and Affiliations

Authors

Contributions

Research conception and proposal design: QW, SC, GZL. Drafting of the manuscript: QW, SC, GZL. Statistical analysis and model development: SC, GZL, GJW. Data collection: SC, GZL, FL. Revision of the manuscript and improvement of the content: all authors. All authors read and approved to submit the final manuscript.

Corresponding author

Correspondence to Qing Wang.

Ethics declarations

Ethics approval and consent to participate

This study was approved by the Ethics Committee of the Affiliated Hospital of Southwest Medical University (KY2022211) and complied with the guidelines of the Declaration of Helsinki. This study was retrospective, the Ethics Committee of the Affiliated Hospital of Southwest Medical University waived the informed consent.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Table S1.

Comparison of rehabilitation factors and prognosis of neurological function.

Additional file 2: Figure S1.

The process of screening independent predictors. Figure S2. Dynamic nomogram for predicting IND after CSCI.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Chen, S., Li, G., Li, F. et al. A dynamic nomogram for predicting the probability of irreversible neurological dysfunction after cervical spinal cord injury: research based on clinical features and MRI data. BMC Musculoskelet Disord 24, 459 (2023). https://doi.org/10.1186/s12891-023-06570-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12891-023-06570-z

Keywords