INTRODUCTION
Traumatic injuries are the leading cause of death among children in Canada and other developed countries. According to the World Health Organization, in 2004, nearly 950,000 children below the age of 18 died across the world as a result of injury. As for Canada, according to the Public Health Agency of Canada, 634 young Canadians ages 19 and under were killed by an injury in 2008. 1 - Reference Yanchar, Warda and Fuselli 4
To better characterize the morbidity engendered by trauma, numerous pediatric trauma scoring systems have been developed in the past few years.Reference Grisoni, Stallion and Nance 5 - Reference Sullivan, Haider and DiRusso 9 However, these trauma scoring systems are often complex, requiring time and trained personnel. Also, they lack validity when applied to different populations, leading to their inconsistent use in pediatric trauma centres.Reference Ott, Kramer and Martus 7
The BIG score is composed of the initial base deficit (BD), international normalized ratio (INR), and Glasgow Coma Scale (GCS), three variables that play a critical role in children’s survival. It is calculated as follows:
The BIG score has been recently developed among 707 children admitted to combat-support hospitals in Iraq and Afghanistan and validated among 1,101 children admitted to a German hospital.Reference Borgman, Maegele and Wade 10 Following these, the BIG score was validated retrospectively using the trauma database of a single tertiary care pediatric hospital.Reference Davis, Wales and Malik 11
The BIG score predicted mortality in children with traumatic injuries, with a high probability of survival in patients with a score of <16.Reference Davis, Wales and Malik 11 This new score is easy to calculate and outperformed previous trauma scoring models in two retrospective studies.Reference Borgman, Maegele and Wade 10 , Reference Davis, Wales and Malik 11
To be widely accepted and disseminated, the BIG score will have to prove validity in multiple settings. The aim of this study was to measure the criterion validity of the BIG score to predict in-hospital mortality among children with blunt trauma admitted to a pediatric emergency department (ED).
METHODS
This was a retrospective cohort study performed between 2008 and 2016 in a single tertiary care pediatric hospital serving a large territory. During the study period, the ED had an annual census varying between 70,000 and 84,000 patient visits per year, with approximately 40 admissions per year at the pediatric intensive care unit (ICU) for traumatic injuries. The study was approved by the research ethics board of our centre.
Study population
Potential participants were all consecutive children admitted to the ED for a blunt trauma. To improve resource utilization, the study was focused on relatively severe patients, namely trauma patients admitted to the pediatric ICU or those who died in the ED. Participants were identified using the Quebec Trauma Registry Information System (SIRTQ). This computerized trauma database is managed by the Ministry of Health and Social Services and undergoes periodic evaluations to identify incorrect data values. It has been used since 1998 to improve the quality of care of trauma patients in Quebec. Trauma patients included in this registry are patients who either died following trauma, were admitted to the ICU, were hospitalized for more than 2 days, or were transferred from another hospital. Medical archivists use standardized coding protocols to extract the registry data from patients’ files. Inclusion criteria were all children less than 18 years of age with a blunt trauma admitted to the ICU or who died at our ED between April 2008 and March 2016. Patients intubated prior to being transferred to the hospital were included if a GCS prior to intubation was available. We excluded children with penetrating trauma as well as patients with known co-morbidities that could influence the BIG score, such as coagulopathies and metabolic, cardiac, and neurologic diseases. Patients who died on the scene or before arriving to the hospital were also excluded.
Outcome
The primary outcome for this study was the association between in-hospital mortality and an initial BIG score ≥16. Secondary outcomes included the association between the BIG score and the probability of transfer to a rehabilitation centre, and the length of hospitalization.
Independent variables
Demographic and clinical information included age, sex, mechanism of injury, associated head injury, transfer from a referral hospital, intubation, blood transfusion, need for surgery, length of hospital stay, transfer to a rehabilitation centre, Injury Severity Score (ISS), death, as well as the initial components of the BIG score available at the referral hospital or at our hospital (base deficit, INR, and pre-intubation GCS). The BIG score was calculated as (BD)+(2.5×INR)+(15-GCS).Reference Borgman, Maegele and Wade 10 , Reference Davis, Wales and Malik 11
Data collection
All charts were reviewed by a member of the research team (CGB), using a standardized report form. This form included information on date of trauma, patient’s age, patient’s sex, mechanism of injury, associated head injury, transfer from a referral hospital, intubation, transfusion, admission to ICU, surgery, length of hospitalization, death, transfer to a rehabilitation centre, and the three components of the BIG score. To ensure the quality of data abstraction, 10% of the charts were reviewed in duplicate by a second rater (JG) blinded to the first evaluation.
Data analysis
Analyses were performed with SPSS v21 software (SPSS Inc., Chicago, Illinois). The primary analysis was the association between the BIG score and mortality for all children for whom all three components of the BIG score were available. Sensitivity and specificity were calculated for a BIG score ≥16. Secondary analyses included the comparison of the BIG score with the ISS. To do so, receiver operating characteristic (ROC) curves were constructed to compare the area under the curve (AUC) of the BIG score with the other scores. Other secondary analyses were to calculate the sensitivity and specificity of the individual components of the BIG score and to evaluate the associations between the BIG score and the probability of transfer to a rehabilitation centre as well as the length of hospitalization using logistic and linear regression analysis. As for children with missing components of the BIG score, an analysis was conducted while imputing missing data by normal values (INR: 1.1, BD: 1.0). The 95% confidence interval (CI) was calculated for each measurement.
Sample size
It was calculated that the inclusion of at least 25 deaths would provide CIs of +/- 0.20 for proportions in the worst-case scenario.
RESULTS
Three hundred ninety-seven patients who suffered from a traumatic injury between April 2008 and March 2016 and who were admitted to the ICU or died in the ED were identified with the computerized trauma database of our hospital. Of these, 336 met the inclusion criteria (Figure 1). Fifty-one (15%) patients had missing data for INR and/or BD. These patients were excluded from our primary analysis and included subsequently in our secondary analysis. One patient had a missing GCS and was excluded from the primary and secondary analysis because it seemed inappropriate to put a normal value to its GCS according to the medical chart.
The median (interquartile range [IQR]) age of the study population was 9 years (IQR), and 65% were male. As described in Table 1, the most frequent mechanisms of injury were motor vehicle incidents (27%) followed by falls from height (22%). The median ISS was 18 (IQR).
ATV=all-terrain vehicle; BIG=composed of base deficit, international normalized ratio, and the Glasgow Coma Scale.
The overall mortality rate was 8.3% (n=28). One patient died in the ED, and 27 died in the ICU. Of the 28 non-survivors, 12 (43%) had a non-accidental injury (child abuse and suicide by hanging), 12 (43%) were involved in a motor vehicle incident, and the remaining suffered from another type of blunt trauma (Table 2).
ATV=all-terrain vehicle; BIG=composed of base deficit, international normalized ratio, and the Glasgow Coma Scale.
The inter-rater agreement for data abstraction was excellent for the 36 charts reviewed in duplicate with kappa scores or intraclass correlation coefficients (ICC) >0.8 for all data. The inter-rater agreement for a BIG score of ≥16 was perfect.
The BIG score was associated with mortality (Table 3). Of the 284 patients without missing data for the BIG score components, 215 patients had a BIG score of <16. Only two (1%) of these patients died. A BIG score of ≥16 demonstrated a sensitivity of 0.93 (95% CI: 0.76-0.98) and a specificity of 0.83 (95% CI: 0.78-0.87) to identify mortality.
Sensitivity: 0.93 CI: 0.77 to 0.98
Specificity: 0.83 CI: 0.78 to 0.87
BIG=composed of base deficit, international normalized ratio, and the Glasgow Coma Scale.
Sensitivity and specificity to identify mortality were also calculated for the individual components of the BIG score. The three individual components of the BIG score demonstrated sensitivities varying between 0.77 and 0.93, whereas the specificities varied between 0.66 and 0.74 (Table 4).
Cut-offs for BD, INR, and GCS were extracted from the literature.Reference Davis, Wales and Malik 11 , Reference Whittaker, Christiaans and Altice 13 , Reference Yousefzadeh-Chabok, Kazemnejad-Leili and Kouchakinejad-Eramsadati 25 , Reference Hindy-Francois, Meyer and Blanot 29
ROC curve for predicting mortality using the BIG score is shown in Figure 2, A (continuous BIG score) and Figure 2, B (dichotomized BIG score). The AUC for these ROC was 0.97 (95% IC: 0.95-0.99) and 0.88 (95% IC: 0.81-0.94), respectively. The AUC for the ISS (Figure 3) was 0.78 (95% IC: 0.71-0.85).
The BIG score was associated to other markers of severity. This is demonstrated by the statistical association between the BIG score and length of hospital stay among survivors using linear regression (p<0.001). As a mean, for each increase of 1 point in the BIG score, the duration of hospitalization increased by 0.9 days (95% CI: 0.5-1.2). Also, using logistic regression, a BIG score of ≥16 was strongly associated to the probability of transfer to a rehabilitation centre for the patients who survived (odds ratio: 6.57 95% CI: 3.24-13.3).
A secondary analysis imputing all missing data as normal for the 51 patients with missing data for INR and/or BD resulted in a sensitivity of 0.93 (95% CI: 0.77-0.98) and a specificity of 0.86 (95% CI: 0.81-0.89) to identify mortality.
DISCUSSION
The BIG score predicted in-hospital mortality among children suffering from a blunt trauma with a sensitivity of 0.93 and a specificity of 0.83 in a study population different from where it was derived and previously validated. A BIG score of <16 identified children with a high probability of survival. This remained similar when adding the patients for whom a physiological component (BD and/or INR) was missing and imputed as normal.
Our results are consistent with previous findings that validated the BIG score as a mortality predictor in pediatric blunt trauma patients, independently of pre-hospital intubation, fluid resuscitation, and head injury.Reference Davis, Wales and Malik 11
Unlike numerous pediatric trauma scoring systems currently used, the BIG score constitutes three physiological measurements that play a critical role in children’s survival. Trauma-associated coagulopathy has been well studied in adults. A quarter of severely traumatized adult patients have abnormal blood coagulation when they arrive in the ED, which is associated with a higher risk of mortality.Reference MacLeod, Lynn and McKenney 12 , Reference Whittaker, Christiaans and Altice 13 This early coagulopathy has been referred to as acute traumatic coagulopathy, early coagulation of trauma, trauma-induced coagulopathy, and acute coagulopathy of trauma-shock.Reference Davenport 14 , Reference Hess, Brohi and Dutton 15 This coagulopathy is primarily caused by the endothelial activation of the protein C pathway influenced by the combination of tissue injury and shock. This leads to rapid anticoagulation and fibrinolysis, which put the trauma patient at risk of severe hemorrhage.Reference Hess, Brohi and Dutton 15 - Reference Cohen, Call and Nelson 17 Other factors influencing coagulation have also been described, such as hemodilution of coagulation factors, which is caused by crystalloid administration during resuscitation, hypothermia, acidosis, anemia, and hypocalcemia.Reference Whittaker, Christiaans and Altice 13 , Reference Hess, Brohi and Dutton 15 , Reference Fries, Innerhofer and Schobersberger 18 Thus, tissue injury and shock are the main drivers of early coagulopathy, whereas other factors exacerbate the coagulopathy.Reference Hess, Brohi and Dutton 15 , Reference Fries, Innerhofer and Schobersberger 18 Hemodilution, hypothermia, and acidosis as contributing factors to coagulopathy have also been described in the pediatric population. However, the mechanism explaining early coagulopathy in traumatized pediatric patients remains unknown. According to Whittaker et al., early coagulopathy is an independent predictor of mortality in severe civilian pediatric trauma patients, particularly in patients with traumatic brain injury (TBI).Reference Whittaker, Christiaans and Altice 13 TBI-associated coagulopathy has been described in a few studies, and several hypotheses have been proposed to explain this association but remain unproven to this day.Reference Keller, Fendya and Weber 19 - Reference Talving, Lustenberger and Lam 21 Strumwasser et al. also found that coagulopathy was associated with an increased risk of death in pediatric patients. According to their study, trauma coagulopathy in pediatric trauma patients tends to develop later compared to adults, suggesting that early coagulopathy in children may be associated with increased mortality.Reference Strumwasser, Speer and Inaba 22
The GCS is used across the world as a tool to assess patients’ level of consciousness, reflecting their degree of brain injury or cerebral hypoperfusion.Reference Cicero and Cross 23 TBI being the most common cause of death in trauma patients, GCS is of great interest in predicting mortality in pediatric trauma patients.Reference Borgman, Maegele and Wade 10 , Reference Dutton, Stansbury and Leone 24 Cicero et al. reported GCS as an independent predictor of mortality in pediatric trauma victims.Reference Cicero and Cross 23 A recent study by Yousefzadeh-Chabok et al. compared the GCS to the pediatric trauma score (PTS) and ISS and found that it was a better predictor of mortality in pediatric trauma patients.Reference Yousefzadeh-Chabok, Kazemnejad-Leili and Kouchakinejad-Eramsadati 25
As for the base deficit, it reflects the shock status in pediatric trauma patients and has been proven to correlate with injury severity, morbidity, and mortality.Reference Jung, Eo and Ahn 26 - Reference Hindy-Francois, Meyer and Blanot 29
The BIG score can be rapidly available on admission to the ED, and it does not require trained personnel to be calculated. With an AUC of 0.97 in our study, the BIG score outperformed the ISS (AUC=0.78). The BIG score performance to predict mortality was significantly better than the INR and BD alone. As for the GCS alone, the sensitivity was the same, but the specificity was lower for the GCS (specificity=0.73) in comparison to the BIG score (specificity=0.83). Furthermore, a GCS score of <9 with a normal INR and normal BD was associated with a lower mortality rate (one deceased patient), suggesting that increased physiological disturbance is associated with mortality.
Head trauma was the leading cause of death in our study. However, the BIG score accurately predicted mortality regardless of the cause of death. Furthermore, the BIG score predicted mortality in different populations (U.S. military in Iraq and Afghanistan, Germans as well as Canadians), suggesting that this score can be applied to different blunt traumas as well as to different populations.Reference Borgman, Maegele and Wade 10 , Reference Davis, Wales and Malik 11 In practice, the BIG score could help the physicians better assess the patient’s physiological state and the severity of its injuries. However, as mentioned by Davis et al., it should not be used to dictate individual treatment.Reference Davis, Wales and Malik 11 Regardless of the BIG score, our entire study cohort, except for the patient who died in the ED, was admitted to the ICU, and patients probably could have died without appropriate management. Therefore, the BIG score is more of an additional tool that could give more information about the patient’s possible outcome, knowing that an elevated BIG score is associated with increased mortality, a longer hospital stay as well as a higher probability of transfer to a rehabilitation centre. The score would also be useful to stratify patients in clinical trials or to compare case severity for studies using an observational design. Finally, the BIG score could potentially be useful to improve resource allocation.
This study has limitations. Firstly, it is a single site retrospective evaluation of the BIG score, making its external validation less accurate. To be fully validated, the score will need to be evaluated in a multisite prospective study. Several patients included in our study had missing components of the BIG score. To avoid selection bias, patients with missing BD and/or INR were included in a secondary analysis (one patient with a missing GCS was excluded). Of these patients, one died. Despite a missing INR, the BIG score for the non-survivor patient was >16. The median ISS for this subgroup was 16, which is less severe than the median for the subgroup without missing data. Thus, we hypothesized that these patients were not tested for coagulopathy and/or acidosis because their condition was not severe enough, and that, if these tests had been done, they would have been normal. We therefore imputed all missing data as normal for the 51 patients with missing INR and/or BD and included them in a secondary analysis. The sensitivity and specificity of the BIG score to predict mortality were similar for the 284 patients without missing data and all of the patients (n=335). Furthermore, to minimize missing data, we selected blunt traumas that were admitted to the ICU and that were consequently more likely to have the three components of the BIG score. Thus, our study population was sicker than the general blunt trauma population, leading to potential bias in our test characteristics. However, while artificially improving the specificity of the score by excluding less sick patients, our study design led to the recruitment of all children having the primary outcome (death). Accordingly, the reported sensitivity should be robust. Moreover, the period of time between the trauma injury and the blood tests varied between patients. Blood was drawn to the referral hospital or to our hospital, making it hard to know whether there is an optimal time after trauma to do the blood test. Another limitation is related to the small number of patients with the primary outcome (death) even though we included all children admitted to a very large pediatric ED for 9 years. This limited our ability to conduct multiple analyses. However, the association between the BIG score and other markers of severity suggests that it is a good marker of severity.
CONCLUSION
The BIG score is a pediatric trauma score that predicts in-hospital mortality in children admitted to the ED following a blunt trauma. A score of <16 accurately identifies children with a high probability of survival. The BIG score is a simple tool that can be rapidly available and provides important information about the patient’s physiological disturbance and the severity of its injuries. To be widely accepted, the BIG score will need to be validated in a prospective cohort study. However, considering the low mortality rate of children who are victims of blunt trauma, this study will have to involve multiple settings.
Competing interests: None declared.