Original Research

Which curve is better? A comparative analysis of trauma scoring systems in a South Asian country

Abstract

Objectives A diverse set of trauma scoring systems are used globally to predict outcomes and benchmark trauma systems. There is a significant potential benefit of using these scores in low and middle-income countries (LMICs); however, its standardized use based on type of injury is still limited. Our objective is to compare trauma scoring systems between neurotrauma and polytrauma patients to identify the better predictor of mortality in low-resource settings.

Methods Data were extracted from a digital, multicenter trauma registry implemented in South Asia for a secondary analysis. Adult patients (≥18 years) presenting with a traumatic injury from December 2021 to December 2022 were included in this study. Injury Severity Score (ISS), Trauma and Injury Severity Score (TRISS), Revised Trauma Score (RTS), Mechanism/GCS/Age/Pressure score and GCS/Age/Pressure score were calculated for each patient to predict in-hospital mortality. We used receiver operating characteristic curves to derive sensitivity, specificity and area under the curve (AUC) for each score, including Glasgow Coma Scale (GCS).

Results The mean age of 2007 patients included in this study was 41.2±17.8 years, with 49.1% patients presenting with neurotrauma. The overall in-hospital mortality rate was 17.2%. GCS and RTS proved to be the best predictors of in-hospital mortality for neurotrauma (AUC: 0.885 and 0.874, respectively), while TRISS and ISS were better predictors for polytrauma patients (AUC: 0.729 and 0.722, respectively).

Conclusion Trauma scoring systems show differing predictability for in-hospital mortality depending on the type of trauma. Therefore, it is vital to take into account the region of body injury for provision of quality trauma care. Furthermore, context-specific and injury-specific use of these scores in LMICs can enable strengthening of their trauma systems.

Level of evidence Level III.

What is already known on this topic

  • While trauma scoring systems have shown potential in identifying in-hospital outcomes of trauma patients, context-specific use of such scores has been limited, especially in low-resource settings. Furthermore, its use has majorly been described in polytrauma patients but a comparison has not been done yet to identify specific trauma scoring systems for various types and body regions of injury. Identifying this in a resource-limited setting has the potential to greatly influence triage and treatment, especially for patients with predominant neurotrauma who can develop major disabilities if management is delayed.

What this study adds

  • Since injuries pertaining to neurotrauma and polytrauma can have significantly different impact on quality of life, it is vital to use injury-specific trauma scoring systems. This study took into account multiple scoring systems and revealed Glasgow Coma Scale to be the better predictor of in-hospital mortality for neurotrauma patients but Trauma and Injury Severity Score to be better for polytrauma patients. Therefore, while these scoring systems prove to be a good prognostic tool for in-hospital mortality, it is vital to take injury specifications into account when using them.

How this study might affect research, practice or policy

  • Several studies have been published relating to the use of trauma scoring systems, but our results open up an avenue to use these scores depending on the type of injury. Furthermore, the majority of these scores have been validated in high-income countries but need validation in other regions of the world, owing to differing characteristics that might impact post-trauma outcomes. Validation and context-specific use of such scores can ultimately strengthen and benchmark trauma systems in resource-limited settings.

Introduction

Injuries contribute to 8% of all global deaths,1 making it a worldwide concern today. While healthcare systems have developed globally, injuries continue to have disastrous effects on the life of trauma patients. Since injuries have been predicted to become the third leading cause of death by 2030,2 it is imperative to identify the gaps in the current trauma system to reduce its mortality. This is specifically significant in low and middle-income countries (LMICs), where injuries continue to be a major burden of both mortality and disability.3 4

As a means to understanding trauma care and improving its quality, trauma registries have been developed at regional, national and international levels.5 One vital component of these registries is the trauma scoring systems, which is used to quantify the degree and severity of injury, predict outcomes and act as a parameter for quality improvement.6 The integration of such scores is important in LMICs, where major challenges to trauma care include rudimentary emergency medical systems, inadequate human resources, financial limitations, and uncoordinated healthcare systems.7 8

Multiple scoring systems are used in LMICs to triage patients who undergo trauma, with Injury Severity Score (ISS) being the most common one.9 This anatomic score incorporates Abbreviated Injury Scale (AIS) and reflects the severity of injuries within different body regions.9 Another score that factors anatomic injuries within its calculation is the Trauma and Injury Severity Score (TRISS), introduced in 1981.10 While ISS and TRISS have been extensively used in high-income countries (HICs), their validity and feasibility have been limited in LMICs.11 This limitation can be attributed to the lack of extensive medical records, radiographic images, and autopsy results in low-resource settings.

Apart from these, Revised Trauma Score (RTS), Mechanism (of injury)/GCS/Age/(systolic blood) Pressure (MGAP) score and GCS/Age/(systolic blood) Pressure (GAP) score have also been introduced as validated predictors of trauma outcomes.12 While these scores are feasible in resource-limited settings, they have yet to be widely used in LMICs. Glasgow Coma Scale (GCS) is another recognized parameter for prediction of in-hospital mortality in trauma patients globally.13 Its use has been significant for patients with traumatic brain injury (TBI); however, it has not been extensively explored yet for polytrauma in LMICs.

Literature from HICs indicates these scoring systems to demonstrate variable predictive values for in-hospital mortality, depending on the body region of injury.14–16 However, their impact has not yet been compared between neurotrauma and polytrauma patients in LMICs, both of which have varying injury etiology, presentation and management. Hence, this study aims to compare trauma scoring systems between neurotrauma and polytrauma patients to identify the better predictor of in-hospital mortality for these types of trauma in a low-resource setting.

Methodology

Study design and setting

A descriptive secondary analysis was conducted using data from a multicenter, prospective trauma registry in Pakistan, an LMIC in South Asia that has a population of more than 225 million people.17 This trauma registry includes two tertiary care hospitals that deal with an extensive load of trauma patients:

  1. The Aga Khan University Hospital (AKUH) is a 760-bed, private-sector hospital and is the first in the country to receive the Joint Commission International accreditation.18 An average of 55 000 emergency patients are treated annually within the hospital.

  2. Jinnah Postgraduate Medical Centre (JPMC) is a public sector hospital with a bed capacity of 2000. Like AKUH, JPMC has its own emergency facility that treats up to 100 000 emergency patients a year.

Study participants and data collection

Data were collected from the trauma registry during a period of 12 months (December 2021 to December 2022). The methodology of collecting data and enrolling patients in the registry is performed as per the Collector Trauma Registry.19 All patients aged 18 and above who underwent trauma were included for analysis.

The variables extracted for analysis included patient demographics (age and gender), details of injury, including its mechanism and anatomic area affected, vital signs including GCS and discharge outcome. The following five trauma scoring systems were calculated to ascertain their prediction of mortality:

  1. ISS: ISS is an anatomic score obtained by taking the sum of the squares of the highest three values of the AIS.20 The AIS codes each injury based on the body region and its severity. ISS was calculated as follows:

    ISS=(1st highest AIS score)2+(2nd highest AIS score)2+(3rd highest AIS score)2

  2. RTS: GCS, SBP and RR were coded from 0 to 4.21 These coded values were then multiplied by pre-specified coefficients for the final RTS as shown below, where higher scores indicate increased probability of survival:

    RTS=(0.9368×GCS Value)+(0.7326×SBP Value)+(0.2908×RR Value)

  3. TRISS: TRISS uses RTS, ISS, and age index to calculate survival. Similar to RTS, a higher TRISS score also depicts greater chances of survival.22 For age less than 55, the age index is coded as 0, whereas for age ≥ 55, it is taken as 1. The following formula was used to calculate TRISS:

    Survival probability=1/(1+e−b),

    Where b=b0+b1×RTS+b2×ISS+b3×Age Index and is calculated based on the mechanism of injury as follows:

    bblunt=−0.4499+0.8085×RTS−0.0835×ISS−1.7430×Age Index

    bpenetrating=−2.5355+0.9934×RTS−0.0651×ISS−1.1360×Age Index

  4. GAP: This physiologic score was calculated using the standardized scoring system,12 as shown in table 1.

  5. MGAP: Similar to GAP, this score uses the standardized scoring system shown in table 1. However, it also includes the mechanism of injury to depict its impact on survival.12

Table 1
|
GAP and MGAP scoring system for coding of variables

Since all scores were calculated during data collection, a quality assurance check was done before statistical analysis. This included accuracy checks on 10% random scores using an online calculator.23

Statistical analysis

Statistical analysis was performed on the Statistical Package for the Social Sciences (SPSS) V.21. For analysis, entries with missing values (demographic and details of injury such as mechanism, cause, etc) were excluded. Analysis was performed for all patients, as well as for patients who presented with neurotrauma and polytrauma separately. This division was done based on the anatomic region affected. Patients who presented with neurotrauma, that is, had no/minor injuries in other regions of the body and were admitted within the department of neurosurgery for treatment, were included as ‘primarily neurotrauma’. All other patients who had a major injury in another region of the body or had multiple major/minor injuries were included as ‘polytrauma’. The latter group did not receive any major neurological/neurosurgical intervention. Continuous data were reported with mean±SD, while categorical data were reported as frequencies and percentages (n/%). Receiver operating characteristic (ROC) curves were used to derive sensitivity, specificity, and area under the curve (AUC) for each scoring system.24 This determined their effectiveness in predicting in-hospital mortality for all three groups (all patients, patients with predominantly neurotrauma, patients with polytrauma) separately.

To ascertain the sensitivity of the original results, multiple imputations (MI) were first performed on the original dataset with all patients for missing values. Five imputation models were developed, each with different imputed values based on the original dataset. Final imputed values were derived after analyzing the imputed datasets separately. Sensitivity analysis was then conducted to compare the OR of non-imputed variables with imputed values,25 specifically age, mechanism and cause of injury, surgical intervention, and trauma scoring systems (TRISS, RTS, GAP, MGAP, and GCS).

Results

Data were extracted for a total of 2817 patients. Patients not meeting inclusion criteria (n=412) and those with missing data (n=398), that is, GCS, SBP, RR, etc, were excluded from the analysis (figure 1). The remaining 2007 patients were included for final analysis to identify trauma scoring systems that are better predictors of in-hospital mortality within our cohort. This cohort of 2007 patients was further analyzed after division into patients with predominantly neurotrauma (n=986) and polytrauma (n=1021) separately.

Figure 1
Figure 1

Inclusion and exclusion criteria for analysis.

The mean age of total patients was 41.2±17.8 years, with the majority of participants being males (80.6%) and residents of the same city as the included hospitals (81.6%). The most common mechanism of injury was blunt trauma (85.2%), while road traffic crash was the most frequent cause of trauma (59.2%). About half of the patients (49.1%) presented predominantly with neurotrauma, while 50.9% were admitted as polytrauma patients. A total of 52.9% patients underwent surgical intervention. This included procedures from various specialties for polytrauma patients, such as open reduction internal fixation, exploratory laparotomy, wound debridement, etc, while for neurotrauma patients, surgical interventions mainly comprised craniotomies. In-hospital mortality rate for the total cohort was 17.2%, while for neurotrauma and polytrauma patients, it was 31.4% and 3.5%, respectively. Further demographic details for all trauma patients, along with specifications for neurotrauma and polytrauma patients, are shown in table 2.

Table 2
|
Demographics, injury details and outcomes of post-traumatic patients

Among the vital signs used for the calculation of trauma scoring systems for total sample, average SBP was 124.4±22.3 mm Hg while the mean RR was 20.3±4.9 breaths per minute. Table 3 shows mean trauma scoring systems for all patients (A), patients with neurotrauma (B) and patients with polytrauma (C), stratified by their outcome.

Table 3
|
Mean trauma scoring systems for enrolled patients with specification for survivors and non-survivors

While the mean ISS for all patients was 12.2±7.3, it was relatively higher for patients with neurotrauma (15.7±7.1) and lower for polytrauma patients (8.8±5.8). This indicated increased severity of injury in patients with neurotrauma. Similar results were seen with other scores, where lower scores imply increased injury severity. For patients with neurotrauma, all these scores were relatively lower than the mean of total cohort, highlighting poor anatomic and physiologic parameters at presentation.

ROC curves were used to predict in-hospital mortality for all neurotrauma and polytrauma patients separately (figure 2). The area under the ROC (AUROC) for each trauma scoring system was calculated, along with its sensitivity and specificity as shown in table 4.

Table 4
|
Sensitivity, specificity, area under the curve with 95% CI, cut-off value and p value for each trauma scoring system (ordered as per AUC)
Figure 2
Figure 2

Receiver operating characteristic (ROC) curves for six trauma scoring systems (Injury Severity Score (ISS); Trauma and Injury Severity Score (TRISS); Revised Trauma Score (RTS); Glasgow Coma Scale (GCS); GCS/Age/Pressure (GAP); Mechanism/GCS/Age/Pressure (MGAP)). (A) All trauma patients. (B) Neurotrauma patients. (C) Polytrauma patients.

While all scores were significant in predicting in-hospital mortality for all patients (p<0.05), physiologic trauma scoring systems (ie, GCS, RTS, GAP and MGAP) proved to be better predictors than ISS and TRISS (table 4A). Among the physiologic scores, GCS proved to be the best predictor of in-hospital mortality with an AUROC of 0.897, followed by GAP (AUC: 0.883), RTS (AUC: 0.882), and MGAP (AUC: 0.874). While TRISS and ISS were significant predictors of mortality, the AUROC for these scores was less than other scores, depicting their lower efficacy in predicting in-hospital mortality for post-traumatic patients. Similar results were seen for patients with neurotrauma, where GCS had the highest AUC (0.885), while ISS had the least (0.724). However, for polytrauma patients, TRISS and ISS showed better prediction of in-hospital mortality as compared with other trauma scoring systems, with an AUC of 0.729 and 0.722, respectively. For patients with polytrauma, GCS had an AUC of less than 0.800 contrary to its value in other groups (AUC: 0.638). While RTS still showed lesser predictability than other scores for in-hospital mortality in patients with polytrauma (AUC: 0.595), it was not significant (p=0.052).

To conduct a sensitivity analysis, MIs were performed for missing data, resulting in 2405 complete entries. Using in-hospital mortality as the outcome, multivariable regression was conducted for both imputed (n=2405) and non-imputed (n=2007) datasets. Sensitivity analysis concluded no statistical difference between the results of regression from the two datasets, indicating reliable results of mortality prediction from the original, non-imputed data.

Discussion

The results from this study indicate the significance of using injury-specific trauma scoring systems in predicting in-hospital mortality within a resource-limited setting. In our study, physiologic trauma scores such as GCS, GAP, RTS and MGAP predicted in-hospital mortality better than ISS and TRISS for the overall cohort as well as for patients with predominant neurotrauma. While GCS and RTS proved to be the best predictor for these groups, TRISS and ISS proved better predictability for in-hospital mortality among patients with polytrauma, indicating the significance of incorporating anatomic parameters for the latter group.

Baker et al initially used ISS in an HIC as a prediction tool of survival in victims of motor vehicle crashes.20 The study concluded ISS to be a good predictor of mortality, determining a higher risk of death with each increasing unit of ISS. Since then, it has been a standard severity score in HICs for the triage of post-traumatic patients to predict their probability of survival using a cut-off value of ≥15 for poor outcomes.26 Furthermore, it has also shown to have a greater association with in-hospital mortality after TBI as compared with GCS.27 In contrast, our study showed ISS to have lesser predictability for in-hospital mortality within the whole cohort and patients with neurotrauma, but better predictive value for patients with polytrauma. This was also seen in literature from other LMICs comparing ISS with other trauma scoring systems (without distinction between neurotrauma and polytrauma), where ISS was the least predictive of in-hospital mortality.28 29 Since the score only caters to the three most affected body regions, the prediction value of ISS might be better demonstrated in patients undergoing blunt trauma of specific body regions,30 31 or in those with isolated abdominal gunshot wounds.32 Hence, scores that also use physiologic parameters, such as TRISS, are considered more accurate.

TRISS has shown variable results in predicting in-hospital mortality for neurotrauma patients, where a Korean study suggests this score to be efficient in predicting survival after TBI,33 while a Thailand study indicates it to be poor in predicting postneurotrauma survival.34 Our study concluded TRISS to have a lower AUC for neurotrauma patients but the greatest AUC for polytrauma patients. Similar to our results, TRISS has shown a higher efficacy in predicting in-hospital mortality for polytrauma patients in HICs. A study from Malaysia indicated TRISS to have a higher AUC of 0.812 than RTS (AUC: 0.802).35 Similarly, another study comparing trauma scoring systems for a trauma registry cohort in the USA highlighted TRISS to have an AUC of 0.938, making it the best predictor of in-hospital mortality in postsurgical trauma patients.36

One of the major reasons for this difference between neurotrauma and polytrauma patients can be attributed to the coefficients used to calculate TRISS. Since these coefficients have been derived from Major Trauma Outcome Study making use of patients from an HIC population,37 it can have limited applicability for neurotrauma patients owing to the lack of stratification within original data. Hence, coefficients for TRISS should be validated for specific populations, rather than the ones developed in HICs, to allow for a more sensitive and specific result. Furthermore, RTS was among the highest physiologic predictors of in-hospital mortality for neurotrauma patients and lowest for polytrauma patients. However, it has shown to underestimate injury severity in LMIC populations.13 The mean RTS of 7.2 in our study for all patients was also comparable with the literature from South Africa, indicating a possibility of prehospital mortality among severely injured patients with a low RTS.38

Our study identified GAP and MGAP as better predictors of in-hospital mortality than ISS for neurotrauma patients. Both GAP and MGAP are validated physiologic scores developed to cater to the gaps brought about by RTS. One of the key factors not taken into account in RTS is the patient’s age, which has been shown to impact their outcome independently.39 Apart from incorporating age, MGAP also includes the mechanism of injury within its score that changes management and outcome assessment methodologies.40 Limited literature using these scoring systems highlights GAP and MGAP as better predictors of in-hospital mortality than anatomic scores for all trauma patients.29 41 Likewise, the AUROC for both GAP and MGAP also corresponds to that specified before from South Asian LMICs, that is, 0.85–0.99 and 0.84–0.99, respectively.29 41 42 Our results for neurotrauma patients are also consistent with the literature, where MGAP had no significant difference in predicting in-hospital mortality for patients with post-TBI when compared with GCS.14 It is interesting to note that while GAP was higher than MGAP for neurotrauma patients, MGAP had a higher AUC than GAP for polytrauma patients. This indicates the significance of including mechanism of injury for management of patients presenting with polytrauma. Since these scores have only been validated in HICs yet, there is a need to further explore them in other settings by comparing with anatomic scores.

The use of GCS as a trauma scoring system is widely implemented in LMICs, owing to its easier calculation without any imaging and its proven significance as a strong predictor of in-hospital mortality.13 Congruent to our results, GCS has shown to be a vital predictive factor of mortality in patients with TBI in both HICs and LMICs.14 43 A retrospective analysis from Egypt showed a mean GCS of 14.06±1.77 for survivors of post-trauma, while for those who died during their hospital stay, it was 8.85±3.84.12 This correlates with our finding for the total sample, where ≤12 was identified as the cut-off between survivors and non-survivors, regardless of the body region of injury. However, for neurotrauma separately, the cut-off was calculated to be 9.5. Similar to this, Watanitanon et al also highlighted the greatest proportion of in-hospital mortality among isolated patients with TBI with a GCS score of 9 as compared with those having a GCS score of greater than 9.44

This study is among the first of its kind from an LMIC to compare trauma scoring scores between neurotrauma and polytrauma patients. Additionally, comparison of GCS with all the scores used in our study has not been done previously within the same patient cohort for the two types of trauma presentation. Sensitivity analysis to negate any possible bias associated with missing data revealed no statistical difference between imputed and non-imputed datasets. This further strengthened our findings. Since the scores were calculated during the enrollment of patients in the study, a random recheck before analysis assured the quality of the data. However, there were some limitations to this study. While this was a multicenter study, both the hospitals belonged to the same province of the country. This can limit generalizability of the results to trauma patients in other regions. There were less than 10% missing data which were removed for analysis; however, any biases which could have resulted from this were reduced with sensitivity analysis, as described above. Furthermore, calculation of TRISS was done using prior validated coefficients in HICs. A validation study done in our cohort to develop TRISS equation could have further increased the strength of our findings.

Conclusion

This study signifies the use of context-specific and injury-specific trauma scoring systems in resource-limited settings as a guide to appropriate triage and treatment strategies. Our results concluded GCS and RTS to show a better prediction of in-hospital mortality in patients with neurotrauma, while TRISS and ISS were better predictors for patients with polytrauma. Therefore, it is vital to take into account the region of body injury for provision of quality trauma care. Furthermore, since there are variable findings in different regions of both HICs and LMICs, it is vital to use context-specific trauma scoring systems depending on differing patient populations. Validation of these scores within LMICs can further enable stakeholders to truly identify the risk of mortality and work towards strengthening and benchmarking their trauma systems.