Abstract
Traumatic brain injury (TBI) is the leading cause of death among trauma patients. Identifying preoperative factors that predict postoperative outcomes in such patients can guide surgical decision-making. The aim of this study was to develop a predictive model using preoperative variables that predicts 30-day mortality and morbidity in patients undergoing neurosurgery following TBI. The American College of Surgeons National Surgical Quality Improvement Program (ACS NSQIP) database was queried between 2005 and 2017 for patients aged 18 years or older who underwent TBI-specific surgery. The primary outcome was 30-day mortality, and the secondary outcome was a composite morbidity score. Significant variables on univariate analysis with Chi-squared test were used to compute multivariable logistic regression models for both outcomes, and Hosmer–Lemeshow test was used. A total of 1634 patients met the inclusion criteria. Most patients were elderly aged >60 years (74.48%), male (63.59%), of White race (73.62%), and non-Hispanic ethnicity (82.44%). The overall 30-day mortality rate was 20.3%. Using multivariate logistic regression, 11 preoperative variables were significantly associated with 30-day mortality, including (aOR, 95% CI) age 70–79 years (3.38, 2.03–5.62) and age >80 years (7.70, 4.74–12.51), ventilator dependency (6.04, 4.21–8.67), receiving dialysis (4.97, 2.43–10.18), disseminated cancer (4.42, 1.50–13.0), and coma >24 hours (3.30, 1.40–7.80), among others. Similarly, 12 preoperative variables were found to be significantly associated with 30-day morbidity, including acute renal failure (7.10, 1.91–26.32), return to OR (3.82, 2.77–5.27), sepsis (3.27, 1.11–9.66), prior operation within 30 days (2.55, 1.06–4.95), and insulin-dependent diabetes (1.60, 1.06–2.40), among others. After constructing receiver operating characteristic curve, the model for mortality had an area under the curve (AUC) of 0.843, whereas composite morbidity had an AUC of 0.716. This model can aid in clinical decision-making for triaging patients based on prognosis in cases of mass casualty events.
Introduction
Traumatic brain injury (TBI) is defined as an alteration in brain function, or other evidence of brain pathology, caused by an external force. 1 TBI is the leading cause of injury related mortality and morbidity, accounting for 30% of all deaths related to injuries. 2 Causes include forceful bumps, blows, head or body jolts, or objects that pierce the skull and penetrate the brain tissue. 3 The resulting head injuries are classified according to their mechanism into the following two types: penetrating or open TBI and nonpenetrating TBI, also known as closed head injury or blunt TBI. 3 With regards to severity, TBI is classified into mild, moderate, and severe based on the Glasgow Coma Score (GCS), which has been implemented as an objective tool to describe the extent of impaired consciousness in acute medical and/or trauma patients. 4 Complications of head trauma such as skull fractures, intracranial bleeds, and increased intracranial pressure (ICP) often necessitate surgical intervention that can be lifesaving and/or used to limit the potential associated morbid events that follow injury.
Several preoperative factors have been used to predict mortality in patients undergoing various surgical procedures. For instance, the American Society of Anesthesiologists physical status classification system is commonly used by clinicians to predict postoperative outcomes using perioperative variables. 5 Other examples of such systems include the Physiological and Operative Severity Score for the enUmeration of Mortality and Morbidity (POSSUM) scale, which predicts risk-adjusted mortality in a variety of surgical procedures, and the Acute Physiology and Chronic Health Evaluation II (APACHEII) severity-of-disease score, which focuses on preoperative laboratory values while taking into account patient age and chronic health evaluation within 24 h of admission to the ICU unit.6,7
In the field of neurosurgery, the NSQIP surgical risk calculator remains a valuable tool for evaluating postoperative risk. 8 However, valid predictive models specific to neurosurgical procedures are still lacking. The case is particularly relevant for patients undergoing emergent surgeries for acute TBI such as emergent evacuation of a subdural hematoma or craniotomy to alleviate elevated ICP. As the morbidity and mortality rates in this patient population remain high, it is essential to identify factors that predict postoperative outcomes to aid in prognostic and surgical judgment, which would prove useful in under-resourced institutions or mass-casualty events where patient triage is of the essence. In this study, we aimed to develop a model based on the American College of Surgeons National Surgical Quality Improvement Program (ACS NSQIP) database that can efficiently predict 30-day mortality and morbidity in patients undergoing neurosurgery for acute TBI.
Materials and Methods
Data source
The ACS NSQIP neurosurgical procedure database was used for this study. The ACS NSQIP collects data on over 150 variables, including preoperative risk factors, intraoperative variables, and 30-day postoperative mortality and morbidity outcomes for patients undergoing major surgical procedures in both the inpatient and outpatient setting. 9 A site’s trained and certified Surgical Clinical Reviewer captures these data using a variety of methods, including medical chart abstraction. 9 In this study, we used the neurosurgical subset of the database, which included a total of 332,971 cases.
Study population
Patients 18 years or older who underwent TBI-related surgeries were selected from the database. The CPT (Current Procedural Terminology) codes used were 61312, 61322, and 61345 corresponding to supratentorial hematoma evacuation (extradural and subdural), decompressive craniectomy to treat raised intracranial pressure (ICP), and decompressive craniectomy of the posterior fossa, respectively. Additional filters using TBI-related diagnoses based on ICD-9 (International Classification of Diseases Ninth Revision) and ICD-10 (International Classification of Diseases Tenth Revision) were used. A total of 1,634 patients met the inclusion criteria and were included in the final analysis (Fig. 1).

Flow diagram depicting study patient selection. NSQIP: National Surgical Quality Improvement Program, TBI: Traumatic Brain Injury, CPT: Current Procedural Terminology, ICD-9: International Classification of Diseases Ninth Revision, ICD-10: International Classification of Diseases Tenth Revision.
Data management and analysis
Continuous variables, including age, hematocrit, platelet count, white blood cell (WBC) count, international normalised ratio (INR), Creatinine, BMI, and length of stay, were categorized based on laboratory cutoffs or clinically meaningful subgroups. The primary outcome was 30-day mortality, whereas the secondary outcome was a composite measure of morbidity. The composite morbidity score included the following four categories of postoperative comorbidities: venous thromboembolism (VTE) composed of two variables (“occurrences pulmonary embolism” and “DVT/thrombo”), infections composed of six variables (“Occurrence of Pneumonia”, “Superficial Site Infection”, “Deep Incisional SSI”, “Organ Space SSI”, “Sepsis”, and “Urinary tract infection”), cerebrovascular accident (CVA) composed of two variables (“Myocardial infraction” and “CVA/Stroke”), and length of stay. The scoring was as follows: 1 point for any VTE, CVA, Infection, or length of stay >14 days. A final morbidity score of 0 corresponded to “No”, whereas a score of 1 or higher corresponded to “Yes.” Finally, composite morbidity was categorized into 2 categories based on the score.
Univariate analysis using chi-square test to identify covariates significantly associated with both mortality and morbidity, with a cutoff of α < 0.05 for significance, was used. Forward and backward logistic regression models with a cutoff of p < 0.1 were computed to identify significant predictors of mortality and morbidity. Hosmer–Lemeshow test of goodness of fit and receiver operating characteristic (ROC) curves were used to assess the predictive capacity of the models. All statistical analyses were conducted using SPSS version 29.0.
Results
Baseline demographics
The neurosurgical database included 332,971 patients between 2005 and 2017. A total of 1,634 patients met the inclusion criteria based on TBI diagnoses and surgical procedures and were included in the study. Figure 1 summarizes the stepwise selection of cases. The majority of patients were elderly, aged >60 years (74.48%), with the >80-year subgroup having the highest sample proportion (30.33%). The cohort was predominantly male (63.59%), of White race (73.62%), and non-Hispanic ethnicity (82.44%). Table 1 summarizes the baseline demographics of the patient cohort.
Baseline Demographics and Predictors of Mortality and Morbidity with Unadjusted Odds Ratios and Significance Levels
No odds ratios were calculated for subcategories for variables with >2 categories. No univariate analysis was done against an outcome for variables having <10 cases associated with that outcome. Names of variables significantly associated with any outcome, associated Odds rations, and p values are in bold.
Primary outcome (30-day mortality)
The following 23 preoperative variables were found to be significantly associated with 30-day mortality after neurosurgery for TBI on univariate analysis: age (p < 0.001), race (p = 0.018), diabetes (OR 1.78, p < 0.001), hypertension (OR 1.74, p < 0.001), current smoking within last year (OR 0.41, p < 0.001), currently on dialysis (OR 3.53, p < 0.001), bleeding disorders (OR 3.57, p < 0.001), disseminated cancer (OR 4.73, p < 0.001), coma > 24 h (OR 30.76, p < 0.001), ventilator dependent (OR 6.50, p < 0.001), emergency case (OR 3.54, p < 0.001), functional health status prior to surgery (p < 0.001), systemic sepsis (p < 0.001), hematocrit (p < 0.001), WBC count (OR 2.29, p < 0.001), platelet count (OR 0.47, p < 0.001), INR (OR 2.93, p < 0.001), Creatinine (OR 2.57, p < 0.001), CVA/stroke (OR 1.97, p = 0.014), history of congestive heart failure (CHF) (OR 2.16, p = 0.019), impaired sensorium (OR 2.03, p < 0.001), open wound/wound infection (OR 1.71, p = 0.039), and transfusion of >4 packed RBCs (OR 2.08, p = 0.031). Table 1 summarizes the results of the univariate analysis across the baseline and preoperative variables for 30-day mortality.
Variables that were significant on univariate analysis (p < 0.05) were plugged into a forward stepwise logistic regression model with a cutoff of <0.1. The model was significant (p < 0.001), and the Hosmer–Lemeshow goodness of fit test showed that the model fits the data well (p = 0.659). The model identified 11 significant predictors of 30-day mortality as follows: Age was a significant predictor of 30-day mortality (p < 0.001). Compared to those <60 years of age, patients in the 60–69, 70–79, and >80-year age groups were 1.66 (95% CI, 0.95–2.88, p = 0.074), 3.38 (95% CI, 2.03–5.62, p < 0.001), and 7.70 (95% CI, 4.74–12.51, p < 0.001) times more likely to die, respectively. Diabetic patients were 1.51 times more likely to die compared with nondiabetics (95% CI, 1.05–2.16, p = 0.026). Patients dependent on ventilators were 6.04 times more likely to die compared with those not on ventilators (95% CI, 4.21–8.67, p < 0.001). Patients receiving dialysis preop were 4.97 times more likely to die compared with those not on dialysis (95% CI, 2.43–10.18, p < 0.001). Patients in a coma >24 h were 3.30 times more likely to die compared with those who weren’t (95% CI, 1.40–7.80, p = 0.006).
Patients with disseminated cancer were 4.42 times more likely to die compared with those who didn’t (95% CI, 1.50–13.00, p = 0.007). Patients with a bleeding disorder were 1.72 times more likely to die than those who didn’t (95% CI, 1.21–2.43, p = 0.002). Emergency cases were 2.05 times more likely to end in death compared with nonemergent cases (95% CI, 1.38–3.05, p < 0.001). Patients with leukocytosis (WBC > 10,000) were 1.72 times more likely to die compared with those with normal WBC count (<10,000) (95% CI, 1.24–2.38, p = 0.001). Patients with an elevated INR (>1.2) were 1.53 times more likely to die compared with those with normal INR (<1.2) (95% CI, 1.05–2.22, p = 0.027). Finally, patients with a normal platelet count (>150,000) were 0.49 times less likely to die compared with those with thrombocytopenia (platelet count < 150,000) (95% CI, 0.34–0.70, p < 0.001). Table 2 summarizes the results of the logistic regression model for predictors of 30-day mortality. Using the predicted probability of our model, an ROC (receiver operating characteristic) curve was constructed with an area under the curve of 0.843 (95% CI, 0.819–0.867, p < 0.001) and standard error = 0.012, which reflect strong discriminatory power of the model developed. The curve is shown in Figure 2.

Receiver Operating Characteristic (ROC) curve of Probabilities for the Preoperative Variables and Mortality at 30 Days.
Multivariate Logistic Regression Model for Predictors of 30-Day Mortality After Undergoing Neurosurgical Procedures Following Acute Traumatic Brain Injury
Significant variable names, associated aOR, and p values are in bold.
Secondary outcome (composite morbidity)
The following 16 preoperative variables were found to be significantly associated with composite morbidity after neurosurgery for TBI on univariate analysis: insulin-dependent diabetes (OR 1.76, p = 0.003), hypertension requiring medication (OR 1.52, p < 0.001), currently on dialysis (OR 2.28, p = 0.008), bleeding disorders (OR 1.77, p < 0.001), ventilator dependent (OR 2.01, p < 0.001), prior operation within 30 days (OR 3.34, p = 0.007), emergency case (OR 1.96, p < 0.001), >10% loss of body weight in last 6 months (OR 2.57, p = 0.009), functional health status before surgery (p = 0.003), systemic sepsis (p < 0.001), hematocrit (p < 0.001), WBC count (OR 1.87, p < 0.001), platelet count (OR 0.60, p < 0.001), acute renal failure (OR 7.86, p < 0.001), open wound/wound infection (OR 1.88, p = 0.007), and return to OR (OR 3.64, p < 0.001). Table 1 summarizes the results of the univariate analysis across the baseline and preoperative variables for composite morbidity.
Variables that were significant on univariate analysis (p < 0.05) were plugged into a backward logistic regression model with a cutoff of <0.1. The model was significant (p < 0.001), and the Hosmer–Lemeshow goodness of fit test showed that the model fit the data well (p = 0.867). The model identified 11 significant predictors of composite morbidity as follows: Patients who had insulin-dependent diabetes were 1.60 times more likely to develop a postop morbidity compared with those who didn’t (95% CI, 1.06–2.40, p = 0.025). Patients who had hypertension requiring medication were 1.35 times more likely to develop a post-op morbidity compared with normotensive patients (95% CI, 1.07–1.71, p = 0.013). Patients who were on a ventilator were 1.38 times more likely to develop a post-op morbidity compared with those not on a ventilator (95% CI, 1.04–1.82, p = 0.025). Emergency cases were 1.45 times more likely to produce a post-op morbidity compared with nonemergent cases (95% CI, 1.11–1.88, p = 0.025). Patients who had >10% body weight loss in last 6 months were 2.29 times more likely to develop a post-op morbidity compared with those who didn’t (95% CI, 1.06–4.95, p = 0.035).
Patients who had SIRS were 1.59 (95% CI, 1.19–2.12, p = 0.002) times and those with sepsis were 3.27 (95% CI, 1.11–9.66, p = 0.032) times more likely to develop a post-op morbidity compared with those who didn’t, respectively. Patients with leukocytosis (WBC > 10,000) were 1.51 times more likely to develop a post-op morbidity compared with those who had normal WBC count (<10,000) (95% CI, 1.18–1.93, p < 0.001). Patients who had acute renal failure were 7.10 times more likely to develop a post-op morbidity compared with those who didn’t (95% CI, 1.91–26.32, p = 0.003). Patients with a normal platelet count (>150,000) were 0.61 times less likely to develop a post-op morbidity compared with those who had thrombocytopenia (<150,000) (95% CI, 0.46–0.81, p < 0.001). Finally, patients who returned to the OR were 3.82 times more likely to develop a post-op morbidity compared with those who didn’t (95% CI, 2.77–5.27, p < 0.001). Table 3 summarizes the results of the logistic regression model for predictors of postoperative morbidity. Using the predicted probability of our model, an ROC curve was constructed with an area under the curve of 0.716 (95% CI, 0.689–0.742, p < 0.001) and standard error = 0.014, which reflect moderate discriminatory power of the model developed. The curve is shown in Figure 3.

Receiver Operating Characteristic (ROC) curve of Probabilities for the Preoperative Variables and Composite Morbidity.
Multivariate Logistic Regression Model for Predictors of Composite Morbidity After Undergoing Neurosurgical Procedures Following Acute Traumatic Brain Injury
Significant variable names, associated aOR, and p values are in bold.
Discussion
Traumatic brain injury is associated with high rates of morbidity and mortality. Although trauma location and its associated complications, such as bleed or fracture, along with surgical technique and approach are crucial in determining surgical adequacy, patient factors such as laboratory results and preexisting conditions provide further insight into surgical eligibility. To the best of our knowledge, this is the first study that uses NSQIP database to identify the predictors of 30-mortality and morbidity among patients undergoing traumatic brain injury procedures.
The overall mortality rate in our cohort was 20.3%. This rate is lower than two similar studies, one conducted at a university hospital in Ethiopia that evaluated all patients with TBI between 2017 and 2022 with an overall mortality of 28.8% and another study in Brazil with a mortality rate of 27.6% for patients presenting with TBI.10,11 This observation could be attributed to the fact that only patients who underwent surgical procedures for TBI were included in our study from the NSQIP database. Patients who are nonsurgical candidates typically have poor outcomes and higher mortality, which could explain why the mortality rate in our surgical cohort was lower than other studies. When considering studies on strictly surgical patients, a retrospective study on over 100,000 geriatric patients from the ACS National Trauma Data Bank in the United States had a mortality rate of 23.9%, which is similar to our results considering that a major proportion of our study population were elderly >60 years of age (74.48%). 11
For the primary end-point (30-day mortality), we identified 11 preoperative variables that can effectively predict outcomes in patients undergoing neurosurgery following acute TBI. Age was one of the most significant predictors, with patients older than 80 years having the highest odds of death (up to 7.7 times) compared with the reference age group (<60 years). This is in line with studies that demonstrated that increased age is associated with increased 30-day mortality following acute TBI. 12 Moreover, ventilator dependence, being on dialysis, having disseminated cancer, and being in a coma for more than 24 h before surgery significantly increased the risk of mortality with adjusted ORs (aORs) of 6.04, 4.97, 4.42, and 3.30, respectively. From a neurological point of view, being in a coma or requiring mechanical ventilation following trauma reflects severe brain damage with possible complications such as brain herniation due to elevated ICP, which would explain why they confer a much higher risk of mortality. In fact, several studies highlight the role of the GCS in predicting mortality after TBI.13,14
Although the NSQIP database does not include GCS, it does provide information about coma status, which is a surrogate for the GCS corresponding to a GCS <8. Therefore, our results suggest that a GCS score ≤8 is significantly associated with mortality. Furthermore, being on dialysis indicates severe renal impairment, which is associated with impaired platelet function and metabolic derangements which can exacerbate the damage from brain injury. 15 Disseminated cancer induces a hypercoagulable state in the patient and weakens the immune response, which may increase mortality from a number of reasons, including thromboembolic events and infectious causes.16,17
Additional variables that significantly increased the odds of death were being an emergent case, having a bleeding disorder, being a diabetic, having a platelet count <150,000, having a WBC count >10,000, and having an INR > 1.2 with aORs 2.05, 1.72, 1.51, 2.04, 1.72, and 1.53, respectively. Diabetes is associated with several comorbidities that significantly increase surgical risk in general and in acute TBI in particular. Bleeding disorders, INR > 1.2, and platelets <150,000 are all factors that increase the risk of significant hemorrhage, which is very critical in patients with brain injury. It is evident from these results that a patient’s medical comorbidities (cancer, diabetes, etc.) can effectively predict the risk of mortality after surgical intervention and can be beneficial in determining surgical eligibility.
The secondary end-point, composite morbidity, was composed of four categories of outcomes in the 30-days postoperatively. Predictors of composite morbidity that were significant included hypertension, sepsis, emergency case status, WBC count, platelet count, and hematocrit. There is evident overlap in the predictors of mortality and morbidity which share similar pathophysiological mechanisms. The distinctive variables for morbidity are hypertension, sepsis, and hematocrit. Hypertension, when associated with complications such as coronary artery disease, congestive heart failure, renal failure, cerebrovascular disease, and diabetes mellitus, can increase surgical risk and morbidity. 18
Sepsis is a life-threatening inflammatory response to bacteremia, which puts the body under significant stress and at higher risk for morbidity after surgery. A prospective cohort study, which evaluated 2.3 million patients in ACS NSQIP database at 374 U.S. hospitals between 2002 and 2012, found that preoperative sepsis increases the risk for arterial and venous thromboses, and the risk increases with severity of sepsis. 19 The same effect on mortality is demonstrated in our results where patients who had SIRS were 1.59 (95% CI, 1.19–2.12, p = 0.002) times and those with sepsis were 3.27 (95% CI, 1.11–9.66, p = 0.032) times more likely to develop a post-op morbidity compared with those who didn’t, respectively. Finally, lower hematocrit before surgery requires transfusion of blood products, which poses the risk for hemolytic transfusion reactions and complications from large volume transfusions. 20
The results of our study hold the highest value in situations where patient triage is of essence. It is imperative that all patients with life-threatening TBI receive surgical intervention, but in scenarios such as natural disasters and mass casualty events, where the number of patients with TBI requiring intervention overwhelms the surgical capacity of health care systems, having a robust predictive model would help triaging patients based on priority. A neurosurgeon has to decide which patient to operate on between the following two: a 75-year-old male patient with long-standing diabetes mellitus receiving hemodialysis for end-stage kidney disease and a 65-year-old male patient with no comorbidities, but whose labs reveal thrombocytopenia. Based on our model, the first patient has a significantly higher odds of death, and the second patient should be prioritized under dire circumstances.
Our study has several potential limitations. The main limitation is the lack of certain variables relevant to neurosurgery for TBI, such as the injury mechanism. The mechanism of injury remains one of the most important predictors of outcomes post brain trauma. For instance, brain injury due to high-speed motor vehicle collisions or high impact blunt trauma (such as falling from a height) confer much more damage and a higher mortality than low impact trauma, but such data were not available for analysis. Another example is GCS on presentation, which is another major predictor of outcomes. GCS is used to categorize patients into different risk groups with different prognoses. Information about the use of anticoagulants and antiplatelets was not found, which are highly relevant in trauma. Another limitation is the short follow-up period dictated by the NSQIP database (30 days). Patients with traumatic TBI often suffer long-term neurological sequelae that require long-term follow-up. The last limitation is the lack of detailed neurological findings in the database. Some patients may develop debilitating neurological symptoms such as weakness, paresthesia, aphasia, or even epilepsy, which could not be computed in our composite morbidity outcome.
Conclusion
The predictive model at hand was established using the NSQIP database, which is a highly reliable resource for surgical data. Furthermore, the model scored high on performance testing which indicates that it can accurately predict outcomes in patients with acute TBI. The clinical value of this model lies in determining the feasibility of operating on patients with TBI and to have an accurate prognostic measure of postoperative morbidity and mortality. This is especially relevant in emergency settings such as disaster management when there is urgent need for triage and determining which patients have a better outlook in the setting of limited time and resources. This prompts decisions to be made about resource allocation and the order in which patients with TBI may be treated. Future prospective studies that include additional variables such as mechanism of injury, GCS, anticoagulant use, and neurological findings with longer follow-up intervals are needed to test this model and enhance it for better predictive capacity.
Footnotes
Acknowledgments
T.J., H.A., F.Y., and H.M. acknowledge the training received under the Scholars in HeAlth Research Program (SHARP) that was, in part, supported by the Fogarty International Center and Office of Dietary Supplements of the National Institutes of Health (Award Number D43 TW009118). The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.
“Research reported in this publication was supported, in part, by the Fogarty International Center of the U.S. National Institutes of Health (NIH), under the Middle East and North Africa Program for Advanced Injury Research (MENA PAIR) training grant number D43TW012202. The content is solely the responsibility of the authors and does not necessarily represent the official views of the NIH.”
Authors’ Contributions
T.J.: Conceptualization (equal), Data curation (lead), Formal analysis (lead), Methodology (equal), Project administration (equal), Writing—original draft (equal), and Writing—review & editing (lead). H.A.: Conceptualization (equal), Data curation (equal), Formal analysis (equal), Methodology (lead), Project administration (equal), Writing—original draft (equal), and Writing—review & editing (equal). F.Y.: Conceptualization (equal), Data curation (equal), Formal analysis (equal), Methodology (equal), Resources, Writing—original draft (equal), and Writing—review & editing (supporting). H.M.: Conceptualization (equal), Formal analysis (equal), Methodology (equal), Resources, Writing—original draft (equal), and Writing—review & editing (supporting). D.H.: Conceptualization (supporting), Formal analysis (supporting), Methodology (supporting), Supervision (equal), Writing—original draft (supporting), and Writing—review & editing (supporting). N.M.: Conceptualization (supporting), Supervision (equal), and Writing—review & editing (supporting).
Data Sharing Statement
Data were acquired from ACS NSQIP database for Neurosurgical procedures, which contain data queried between 2005 and 2017. For more details about data acquisition, characteristics of the analysis, and source and reproducibility of the data acquisition approaches, refer to ACS NSQIP Participant User File (PUF) 2017 using the following link: (
).
Institutional Review Board Approvals/Waivers
No IRB approval was required for this study.
Funding Information
This research received no funding from any agency, commercial, or nonprofit sectors.
Author Disclosure Statement
The authors declare no conflict of interest regarding competing interests, personal financial interests, funding, employment, or other competing interests that may arise related to this work.
