Abstract
Background
Blood urea nitrogen (BUN) reflects renal function and metabolic catabolism and thus can be a potential biomarker in severe illness. The purpose of this investigation was to assess whether independent associations exist between admission BUN concentrations and hospital mortality in critically ill adults with acute myocardial infarction (MI) or cardiogenic shock (CS). Associations were assessed using tertiles of mean BUN during the hospital admission and mean BUN within 24 h of ICU entry.
Methods
Using MIMIC-IV (v3.1), we conducted a retrospective cohort study of first ICU admissions (2008–2019) in adults identified by ICD-9/10 codes for acute myocardial infarction and/or cardiogenic shock. The primary exposure was mean BUN during hospitalization, categorized into tertiles; first-24-h BUN was analyzed as a sensitivity exposure. Analyses used complete cases.
Results
We included 6719 patients. In-hospital mortality climbed from 5.2% to 25.8% across average BUN tertiles. In multivariable logistic regression, higher mean BUN was associated with increased odds of in-hospital mortality across tertiles and per standard-deviation increase. Findings were directionally consistent when first-24-h BUN was used as the exposure. Model diagnostics supported acceptable collinearity and fit.
Conclusions
Among critically ill adults with acute myocardial infarction and/or cardiogenic shock, higher BUN, particularly mean BUN during hospitalization, was independently associated with in-hospital mortality. Early BUN (first 24 h) showed a similar signal, supporting BUN as a simple adjunct for bedside risk assessment.
Keywords
Introduction
Acute myocardial infarction (AMI) is responsible for nearly one-third of all fatalities and stands as the primary cause of death on a global scale. 1 The occurrence of AMI accompanied by cardiogenic shock (CS) is a critical syndrome marked by a swift deterioration in cardiac function. 2 Shock is a state of circulatory failure that leads to an insufficient supply of oxygen, resulting in potentially life-threatening tissue hypoxia. 3 Cardiogenic shock (CS) is particularly complex and critical; it arises from inadequate cardiac output, resulting in poor tissue perfusion and limited oxygen delivery. This deficiency can ultimately cause damage to vital organs.4,5 Among AMI patients, CS occurs in 5% to 10% of cases.6–8 Despite significant reductions in AMI incidence and mortality,9,10 outcomes in patients with CS remain poor. 11
Even with advancements in reperfusion, in-hospital mortality rates for CS remain alarmingly high, ranging from 27% to 51%, with nearly half of deaths occurring within the first 24 h.5,12,13 This implies that the risk of death is time-sensitive. Therefore, it is important to early recognize patients at high risk of myocardial infarction and cardiogenic shock. 13 Furthermore, research has shown that long-term survival in AMI-CS patients is poor with little improvement over time, and survivors are destined to develop significant morbidity and also an increased risk for readmission and mortality. 14 Hence, there exists a critical need for simple and effective methods to identify and stratify at-risk patients for increased mortality.
Blood urea nitrogen (BUN) is a biochemistry indicator that portrays protein metabolism and renal excretion processes primarily occurring in the liver and kidneys. It is also utilized to calculate protein intake and evaluate kidney function.15,16 Increased BUN levels not only reflect impaired cardiorenal function but are also associated with neurohormonal activation, and therefore BUN is a valid adverse prognostic indicator in various clinical conditions.17,18 In cardiogenic shock, impaired left ventricular function reduces cardiac output, and decreased coronary perfusion activates the neurohormonal system, which initiates vasoconstriction, increased heart rate, and fluid retention, all of which exacerbate myocardial ischemia. This pattern creates a self-perpetuating cycle of decreasing shock and compromise in myocardial function.19,20 Therefore, BUN concentration in CS may serve as a marker of the severity of underlying pathophysiology and provide valuable prognostic information.
Increased BUN concentrations have also been seen by other studies to independently predict mortality in critically ill patients, which provides additional evidence for its use as a prognostic marker. 21 However, most of these studies were limited by small sample sizes, and few had examined the association between BUN levels and in-hospital mortality in cardiogenic shock patients using large-scale databases. To address this gap, the present study is the first to utilize the MIMIC-IV database, offering a more recent and comprehensive dataset of critically ill CS patients. By leveraging this updated resource, our research aims to provide current evidence on the prognostic value of BUN, contributing to improved risk stratification and supporting clinical decision-making in the management of cardiogenic shock.
Methods
Study design and data source
This was a retrospective cohort study based on data from the Medical Information Mart for Intensive Care IV (MIMIC-IV) database, version 3.1. MIMIC-IV is a large, de-identified, publicly accessible database of extensive health-related information on patients admitted to critical care units in Beth Israel Deaconess Medical Center from 2008 to 2019. 22
Study population and inclusion/exclusion criteria
This retrospective cohort study used routinely collected intensive care unit (ICU) data. We first identified all ICU stays and retained only the first ICU admission per patient (ranked by ICU admission time) to avoid within-patient clustering. The index stay was required to have a recorded initial ICU type (first care unit). We then restricted the population to adult patients (age ≥18 years) with evidence of acute myocardial infarction and/or shock syndromes during the index hospital admission, defined using ICD-9 and ICD-10 diagnosis codes. Eligible admissions included acute myocardial infarction (ICD-9 410* or ICD-10 I21*), cardiogenic shock (ICD-9 78551 or ICD-10 R57.0), or other shock categories (ICD-10 R57* excluding R57.0).
To focus on clinically comparable critical illness and ensure meaningful early ICU exposure, we required an ICU length of stay of at least 1 day. We further required at least one blood urea nitrogen (BUN) laboratory measurement during the hospital admission, because BUN was the principal exposure. Patients were excluded if they had diagnosis codes consistent with severe liver disease or metastatic cancer, given the strong, non-modifiable influence of these conditions on BUN metabolism and short-term mortality. Finally, we excluded patients who received renal replacement therapy (RRT) during the hospital admission, identified from ICU charted therapy indicators and procedure codes, to reduce confounding by treatment-driven changes in nitrogen balance. The analytic cohort additionally required non-missing values for the primary exposure and the in-hospital mortality outcome.
Variable definitions and data extraction
The data were accessed using Google BigQuery by running a custom SQL script. Demographics included age (years), sex (coded as male vs non-male), and ethnicity as recorded at hospital admission. ICU characteristics included the first care unit and ICU length of stay. The primary outcome was in-hospital mortality, derived from the hospital expiration flag.
The primary exposure, average BUN during the hospital admission (bun_mean_all), was calculated as the mean of all available blood BUN values. A secondary exposure, BUN in the first 24 h after ICU admission (bun_first24), was computed using BUN measurements whose timestamps fell within 24 h of ICU entry. BUN was identified from blood urea nitrogen laboratory item identifiers. Additional laboratory covariates (electrolytes, glucose, creatinine, anion gap, hematology indices, and coagulation tests) were extracted from blood-based measurements and summarized as admission-level means. Vital signs were obtained from ICU chart events within the first 24 h after ICU admission and summarized as the maximum recorded value for each parameter over that window.
Statistical analysis
We documented cohort assembly using the selection flags so that each eligibility step could be reproduced directly from the dataset. Baseline characteristics were described overall and compared across tertiles of bun_mean_all, using cut points at the 33.3rd and 66.7th percentiles. Continuous variables were summarized as mean ± standard deviation or median (interquartile range), depending on distributional shape, while categorical variables were summarized as counts and percentages. Between-tertile comparisons used standard tests matched to data type and distribution (e.g. ANOVA or non-parametric alternatives for continuous variables, and chi-square tests for categorical variables).
The association between BUN and in-hospital mortality was evaluated using logistic regression and reported as odds ratios with 95% confidence intervals. BUN was modeled in two complementary ways: (1) per one standard deviation increase (standardized BUN) to provide a scale-free effect estimate, and (2) as tertiles (T2 and T3 compared with T1) to support clinical interpretability. We fitted sequential models to show how adjustment changed the estimate: unadjusted models; age- and sex-adjusted models; and multivariable models additionally adjusting for available early hemodynamic and laboratory covariates. Each model was fit on complete cases for its included predictors; covariates with no usable values were omitted from the corresponding model.
To support interpretability and model stability, we performed diagnostic checks that were aligned with our research question. We screened multicollinearity using variance inflation factors, and we probed departures from linearity in the log-odds by adding a quadratic term for standardized BUN. Discrimination was assessed with receiver operating characteristic analysis, reporting the area under the curve with a 95% confidence interval and an optimal predicted-probability cutoff based on Youden's index. Sensitivity analyses repeated the same modeling framework using bun_first24 (per standard deviation and tertiles) to test robustness to exposure timing. Finally, we evaluated potential non-linearity using a restricted cubic spline for bun_mean_all (four knots at the 5th, 35th, 65th, and 95th percentiles) and visualized the adjusted predicted probability of in-hospital mortality with 95% confidence intervals across the observed BUN range. All analyses were conducted with R 4.5.2.
Results
Characteristics of study population
After applying the prespecified inclusion and exclusion criteria, we identified 6719 adult ICU admissions with myocardial infarction and/or shock and at least one blood urea nitrogen (BUN) measurement for the analytic cohort (Figure 1). Patients were divided into tertiles of average BUN during the hospital stay.

CONSORT flow diagram—sample selection.
Baseline characteristics differed substantially across BUN tertiles (Table 1). Patients in higher tertiles were older and more often male. The distribution of ethnicity and ICU type also shifted with rising BUN, with a higher proportion of White patients and a greater share of non-coronary ICU admissions in the upper tertile.
Baseline characteristics of the analytic cohort according to tertiles of average BUN during the hospital stay.
P-values reflect overall differences across BUN tertiles based on appropriate statistical tests (ANOVA or Kruskal–Wallis for continuous variables and chi-square tests for categorical variables).
Hemodynamic variables showed modest but consistent differences. Compared with T1, patients in T3 had lower diastolic and mean arterial blood pressures and slightly higher respiratory rates, while heart rate remained broadly similar. Temperature and oxygen saturation were marginally lower in the highest tertile, although absolute differences were small.
Laboratory values displayed a clear pattern compatible with more advanced renal and systemic dysfunction at higher BUN levels. Creatinine increased steeply across tertiles. Electrolytes and acid–base parameters also changed: potassium and anion gap rose, whereas chloride and bicarbonate tended to fall with higher BUN. Markers of anemia and coagulopathy were more abnormal in T3, with lower hemoglobin and hematocrit and higher aPTT, PT, and INR. Platelet counts and white blood cell counts showed small but statistically significant shifts, suggesting a more inflamed and coagulopathic phenotype in patients with elevated BUN.
In-hospital mortality increased sharply across BUN tertiles. The crude mortality rate was 5.2% in T1, 8.6% in T2, and 25.8% in T3 (p < 0.0001), indicating that patients with the highest average BUN experienced nearly a five-fold higher absolute risk of death than those in the lowest tertile.
Relationship of BUN to hospital mortality
In unadjusted logistic regression, each one–standard deviation (SD) increase in average BUN was associated with an odds ratio (OR) for in-hospital death of 1.88. When BUN was modeled categorically, T2 and T3 were strongly associated with mortality compared with T1 (OR 1.70 and 6.29, respectively; both p < 0.0001). Adjustment for age and sex produced only modest attenuation of these associations.
In the fully adjusted model, controlling for demographic factors, vital signs, laboratory variables, and ICU type, the strength of the association remained pronounced, although the gradient became more concentrated in the highest tertile. Each SD increase in average BUN conferred an OR of 2.48. Compared with T1, T3 retained a robust association with mortality OR 3.01, whereas the association for T2 was no longer statistically significant (OR 1.20, 95% CI 0.69–2.11; p = 0.52).
The discriminative performance of the fully adjusted BUN model was excellent. The area under the ROC curve was 0.906 (Figure 2). At the optimal predicted-probability cutoff derived from Youden's index (0.071), sensitivity was 86.3% and specificity was 80.2%, indicating that the model correctly identified the majority of patients who died while maintaining good specificity among survivors (Table 2).

Receiver operating characteristic (ROC) curve for the adjusted BUN model. The AUC ROC was 0.906 (95% CI 0.884–0.928). Using Youden's index, the optimal predicted-probability cutoff was 0.071, with sensitivity 86.3% and specificity 80.2%.
Association between average BUN during the hospital stay and in-hospital mortality.
This table presents odds ratios (with 95% confidence intervals) for in-hospital mortality associated with average BUN during the hospital stay (bun_mean_all), modeled both per one standard deviation increase and as tertiles (T2 and T3 compared with T1). Unadjusted, age- and sex-adjusted, and fully or fallback adjusted models are shown. Each model reports the complete-case sample size (N) and the number of deaths to aid interpretation.
Sensitivity and spline analyses
Sensitivity analyses using BUN measured during the first 24 h after ICU admission yielded consistent findings (Table 3). In the fully adjusted model, each SD increase in first-24-h BUN was associated with an OR of 1.56. The highest tertile of early BUN remained significantly associated with mortality compared with the lowest tertile (OR 1.89), whereas the middle tertile again showed no clear excess risk.
Sensitivity analysis: association between BUN in the first 24 h and in-hospital mortality.
This table presents sensitivity analyses in which the exposure is BUN measured during the first 24 h after ICU admission (bun_first24). BUN is modeled both per standard deviation increase and in tertiles (T2 and T3 compared with T1), using unadjusted, age- and sex-adjusted, and fully or fallback adjusted logistic regression. Model-specific complete-case sample size (N) and deaths are shown to facilitate interpretation.
Finally, restricted cubic spline analyses did not provide evidence of a strongly non-linear relationship between average BUN and mortality after adjustment for covariates (Figure 3, Table 4). Adding linear BUN to the covariate-only model significantly improved model fit (χ2 = 39.09, p < 0.0001), whereas allowing for spline terms did not materially reduce the AIC or significantly enhance fit. These findings support a largely linear increase in the odds of in-hospital death across the observed range of BUN, with the greatest absolute risk concentrated among patients in the highest tertile.

Non-linear association between average BUN and in-hospital mortality.
Metrics for spline model.
Discussion
Summary of findings
The current study, based on the large dataset of MIMIC-IV, investigated the correlation between BUN levels of ICU-admitted patients with CS, both categorical and quantitative, and in-hospital mortality. The outcomes of multivariate logistic regression analysis showed that higher levels of BUN were significantly related to in-hospital mortality and worse prognosis in CS patients.
Comparison with prior research
In this study, we stratified patients by BUN tertile and found that patients in the highest tertiles had a higher death rate than those in the lowest tertiles, which suggests that the BUN prognostic value as a continuous variable may be underestimated. These findings are consistent with the increasing evidence linking renal dysfunction markers to adverse outcomes in critical care settings. 23 Furthermore, the association was maintained after adjustment for confounders such as age, sex, vital signs, and laboratory parameters included in the fully adjusted model, suggesting that BUN plays an independent role in stratification.
Aligned with our results, Zhang et al. used the MIMIC-IV and eICU-CRD databases to develop an explicable machine learning (ML) model that identified BUN as a key predictor of hospital mortality in patients with AMI. 24 In a similar study of AMI cases and prediction of in-hospital mortality by ML models, BUN and CS were both important predictors and related to inpatient mortality. 25 This study's findings are consistent with the previous literature in the MIMIC-III database,26,27 e.g. lower BUN or higher BUN-to-Cr ratios were associated with a reduced risk of mortality in patients admitted to ICU with CS.26,27 Moreover, supporting evidence has evaluated the threshold value of BUN for the mortality risk of CS patients.28,29 In the study by Zhu et al., BUN levels above 8.95 mmol/L significantly increased the 30-day mortality in AMI-CS patients. 28 In addition, Yamga et al. have developed a checklist-based risk score including a BUN limit of 25 mg/dL (≈8.93 mmol/L), which is consistent with our findings and improves its clinical utility. 29 These consistent results across different datasets and methodologies reinforce the evidence for the usefulness of BUN as a reliable prognostic marker.
Several possible mechanisms could explain the elevated BUN levels and higher risk of adverse events in patients with MI or CS. First, during acute myocardial infarction—especially when complicated by cardiogenic shock—a sudden hypoperfusion reduces Glomerular Filtration Rate (GFR), increases filtration parameters, such as BUN or serum Cr, and contributes to kidney injury.21,30 Consequently, Acute Kidney Injury (AKI) as a complication of CS has been shown to correlate with a higher mortality risk. 31 Second, high BUN levels may also be related to low blood volume in CS cases. Diuretic use, vomiting, diaphoresis, and poor intake in MI patients, combined with dehydration, can further increase mortality in critically ill patients.20,28 Third, reduced cardiac output during MI or CS activates neurohormonal mechanisms, which include sympathetic nervous system (SNS) activation, renin-angiotensin-aldosterone system (RAAS) stimulation, and arginine vasopressin (AVP) release. 32 These may raise BUN by vasoconstriction, sodium and water retention, and urea reabsorption in renal tubules.19,32,33 In addition, Jujo et al. demonstrated that persistent high BUN levels are associated with increased cardiovascular death in acute heart failure, a finding that is relevant for our CS cohort, where sustained BUN elevations after admission predicted a lower survival. 33 Together, these mechanisms highlight the multifunctional role of BUN in reflecting both renal and hemodynamic derangement. 34
Some studies have evaluated biochemical markers in CS patients after administration of treatment modalities such as medication or dialysis.35,36 Wang et al. examined the effect of intravenous levosimendan in acute HF patients. 35 Levosimendan significantly decreased B-type natriuretic peptide (BNP) and increased urine output without a meaningful reduction in BUN or lactate levels of the patients after administration. 35 Another study performed peritoneal dialysis in patients with cardiorenal syndrome type 1, of whom 68% had CS. 36 However, according to the results of multivariate analysis, BUN and Creatinine (Cr) levels at dialysis initiation were not correlated with in-hospital mortality. 36 These different outcomes compared with our results might be attributed to the limited sample size, different study settings or interventions, and focus on baseline BUN.35,36 Future studies should look at the longitudinal BUN trends after treatment to reconcile these differences.
In addition, prospective studies should investigate whether targeted interventions based on BUN levels, such as individualized fluid therapy and renal protection strategies, can improve patient outcomes. As BUN is influenced by factors such as diet and medication use, the integration of BUN with other biomarkers such as albumin and inflammatory markers may provide a more comprehensive assessment in patients with critical disease. 21
In this study, BUN was associated with mortality in critically ill patients in a fully adjusted model, which is in line with previous studies such as Biere et al., which concluded that elevated BUN in patients with normal creatinine was independently associated with increased mortality in critically ill patients. 21 This adjustment minimizes the confounding effects of demographic and physiological variables, indicating that BUN independently contributes to the adverse outcomes. 37
Strengths and limitations strengths
Retrospective design: This retrospective observational analysis can demonstrate association but cannot establish causality. Although we adjusted for a broad set of demographic, haemodynamic, and laboratory covariates, residual confounding from unmeasured factors (e.g. fluid balance, vasoactive medication dosing, and clinician decision-making) may persist. Future prospective studies and external validation cohorts are needed to confirm clinical utility and to test whether BUN-guided strategies improve outcomes. Single-center data and generalisability: MIMIC-IV reflects practice at a single tertiary center, and our eligibility criteria (e.g. excluding very short ICU stays and patients with conditions that markedly distort BUN interpretation) may limit generalisability. Therefore, the magnitude of effect estimates may differ in other health systems or patient mixes. We emphasize the need for multi-center validation before applying BUN thresholds for decision-making. Data completeness and case definition: As with many large electronic health record datasets, some variables were missing or inconsistently recorded. In the primary models, we used complete-case analyses and reported the analysed sample size for transparency; if missingness is not completely at random, bias is possible. Additionally, diagnoses were identified using ICD codes, which may introduce misclassification—particularly for cardiogenic shock. Where feasible, future work should incorporate physiologic criteria (e.g. hypotension, lactate, vasopressor requirements) or validated phenotyping algorithms. BUN measurement timing and clinical implementation: Using the average BUN across the hospital admission may capture overall illness trajectory but is less directly actionable at the bedside. To address this, we performed a sensitivity analysis using BUN within the first 24 h after ICU admission (a proxy for admission-time assessment), which showed consistent direction and statistical significance. This supports the potential role of early BUN for practical risk stratification, while recognizing that prospective validation is required.
Clinical implications
The results indicate that BUN, a simple and inexpensive laboratory test, can be an important adjunct for early risk stratification in critically ill patients with acute MI or cardiogenic shock. Clinicians should be alerted by elevated levels of BUN to view the patient's global physiology, not solely renal function, including volume status, catabolic load, and systemic perfusion. Incorporating BUN into standard risk assessment algorithms for such patients may assist in the identification of patients at increased risk of hospital mortality, with the potential for more intensive monitoring, earlier institution of supportive therapies, or the consideration of advanced interventions. This may lead to better patient outcomes through more focused and timely interventions.
Future research
It is critical to perform future prospective studies to validate these findings in external, multi-center populations to enhance generalizability. Investigations that study the temporal variation in BUN levels and their association with mortality can provide further insight into its prognostic value. Additionally, mechanistic studies are needed to delineate the precise pathophysiological processes that link high BUN levels with adverse outcomes in patients with MI and CS, beyond conventional markers of renal impairment. Determination of the impact of targeted therapies aimed at normalizing BUN levels (e.g. vigorous fluid resuscitation and nutritional support) on clinical outcomes in this group would also be of paramount significance. Finally, the determination of the predictive value of BUN when combined with other novel biomarkers or advanced risk prediction tools may increase prognostic precision.
Conclusion
Increased blood urea nitrogen levels are a strong and independent predictor of hospital mortality in critically ill adult patients with acute myocardial infarction or cardiogenic shock. This investigation identifies BUN as an important, readily available, and inexpensive risk stratification biomarker in this high-risk group. Our results emphasize the need to interpret BUN levels not solely as a renal function marker but also as a global marker of severity of illness and compromised physiological reserve. These observations have practical implications for clinical decision-making and the design of improved patient management strategies in acute cardiac emergencies.
Supplemental Material
sj-txt-1-cvd-10.1177_20480040261434635 - Supplemental material for Blood urea nitrogen and in-hospital mortality in critically ill patients with MI or cardiogenic shock: Analysis of the MIMIC-IV database
Supplemental material, sj-txt-1-cvd-10.1177_20480040261434635 for Blood urea nitrogen and in-hospital mortality in critically ill patients with MI or cardiogenic shock: Analysis of the MIMIC-IV database by Sadaf Derakhshandeh, Sepehr Ramezanipour, Golnaz Yazdanpanah, Sahar Nasrollahi, Sina Esmaeili, Hamed Mottaghi, Atefeh Bahrambeigi, Mahmoud Abdollahi, Alaleh Alizadeh, Sepideh Hadimaleki and Danyal Yarahmadi, Niloofar Deravi in JRSM Cardiovascular Disease
Supplemental Material
sj-txt-2-cvd-10.1177_20480040261434635 - Supplemental material for Blood urea nitrogen and in-hospital mortality in critically ill patients with MI or cardiogenic shock: Analysis of the MIMIC-IV database
Supplemental material, sj-txt-2-cvd-10.1177_20480040261434635 for Blood urea nitrogen and in-hospital mortality in critically ill patients with MI or cardiogenic shock: Analysis of the MIMIC-IV database by Sadaf Derakhshandeh, Sepehr Ramezanipour, Golnaz Yazdanpanah, Sahar Nasrollahi, Sina Esmaeili, Hamed Mottaghi, Atefeh Bahrambeigi, Mahmoud Abdollahi, Alaleh Alizadeh, Sepideh Hadimaleki and Danyal Yarahmadi, Niloofar Deravi in JRSM Cardiovascular Disease
Footnotes
List of abbreviations
Acknowledgment
The authors would like to thank the researchers whose work was included in this study.
Ethics approval and consent to participate
MIMIC-IV contains de-identified health-related data and is made available through PhysioNet. Access requires completion of human-subjects training and acceptance of the data use agreement. Because the dataset is de-identified, this analysis was considered non-human-subjects research by the investigators’ institution; therefore, individual informed consent was not required.
Authors’ contributions
Funding
The authors received no financial support for the research, authorship, and/or publication of this article.
Declaration of conflicting interests
The authors declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Availability of data and materials
The dataset(s) supporting the conclusions of this article are available in the PhysioNet repository, maintained by MIT Laboratory for Computational Physiology. The data can be accessed from:
Access to the dataset requires completion of a data use agreement and certification in human research.
Data citation
Johnson AEW, Bulgarelli L, Shen L, et al. MIMIC-IV, a freely accessible electronic health record dataset. Scientific Data. 2023;10:1. doi: 10.1038/s41597-022-01899-x
Database and software information
Project name: MIMIC-IV
Archived version: DOI: 10.13026/s6n6-xd98
Operating system(s): Platform independent
Programming language: SQL, Python (optional for analysis)
License: PhysioNet Credentialed Health Data License
Restrictions for non-academics: Use limited to credentialed researchers; approval required
Supplemental material
Supplemental material for this article is available online.
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
