Abstract
Abstract
Background:
Identifying people who are at risk of deteriorating and dying is essential to inform goals of care, appropriate treatment decisions, patient autonomy, and effective end-of-life care. Limited literature exists on predicting survival near the end of life for people with a hematological malignancy.
Objective:
To identify the key clinical indicators that signal a person with a hematological malignancy is at high risk of deteriorating and dying.
Design, Setting, Participants:
Eleven clinical indicators identified in a Delphi approach were tested via a retrospective case–control study. Each indicator was assessed for at each in-patient admission between living (n = 236) and deceased (n = 120) people with a hematological malignancy who were admitted to a large tertiary hospital between 1st July 2014 and 31st December 2015.
Results:
Six clinical indicators were independently associated with mortality in the final three months of life: declining performance status (Odds Ratio [OR] 7.153, 95% Confidence Intervals [CI] 3.281–15.597, p = < 0.001); treatment limitations of the hematological malignancy (OR 7.855, 95% CI 3.528–17.489, p = < 0.001); relapse, refractory or persistent disease (OR 3.749, 95% CI 1.749–8.039, p = 0.001); presence of two or more comorbidities (OR 2.991, 95% CI 1.319–6.781, p = 0.009); invasive fungal infections (OR 4.887, 95% CI 1.197–19.949, p = 0.027); and persistent infections (OR 6.072, 95% CI 2.551–14.457, p = < 0.001).
Conclusions:
This study has identified six clinical indicators that signal a person with a hematological malignancy is at high risk of deteriorating and dying and may benefit from an assessment of palliative needs and proactive planning, along-side appropriate treatment.
Background
Identifying when a person is likely to be nearing the end of life is vital to facilitate appropriate treatment decisions, patient autonomy, resuscitation planning, advance care planning, and best practice care at the end of life. 1 Accurate identification of risk can allow clinicians to have honest sensitive conversations with patients and their families around the potential for deterioration and death, whilst continuing current appropriate care.1,2 “The right conversation with the right people at the right time can enable patients and their loved ones to make the best use of the time that is left and prepare for what lies ahead”.3 (p. 653) Timely recognition of risk of dying allows time for the involvement of specialized healthcare staff, such as palliative care teams or those knowledgeable in advance care planning. 3
An international body of literature highlights that there are challenges providing palliative care and transitioning to end-of-life care for people with a hematological malignancy.4,5 This is largely due to difficulties predicting survival as people with a hematological malignancy often experience a fluctuating and unpredictable illness trajectory.4–7 Several authors have highlighted the need for greater research into prognostic variables and clinical tools to help clinicians predict survival times near the end of life for people with a hematological malignancy.4,8
The body of literature concerning prognosticating at the end of life for people with a hematological malignancy has largely focused on people admitted to the intensive care unit or those who are still receiving aggressive treatment. 9 To date, few studies have developed or tested prognostic models at the end of life and identifying risk of deteriorating and dying for people with a hematological malignancy for the purposes of informing palliative care provision and transitioning to end-of-life care.10–14 There is an urgent need for more work in this area.
Objectives
This study aimed to identify clinical indicators associated with deteriorating and dying for people with a hematological malignancy and identify when these indicators first present before death. The overall aim was to identify the most appropriate time to engage in honest, sensitive discussion with patients and their families regarding possible risk of deteriorating and dying and enable greater patient autonomy and smoother transitioning to the end of life.
This study did not aim to identify when any form of palliative care should be integrated, as palliative care is recommended to be integrated at any time in the illness trajectory on a needs basis. 15 Rather, this study aimed to identify when a person with a hematological malignancy is at high risk of deteriorating and dying and requires an assessment of current and potential palliative care needs in the context of possible imminent deterioration and death. Ideally, palliative care (including referral to specialist palliative care service) should occur early in the disease trajectory (including upon diagnosis) for people with serious illness, as this has been shown to be beneficial for patients and carers.16–18
Terminology
In this study, the primary outcome was death, not the process of deteriorating and dying. However, the term “deteriorating and dying” was used to align with internationally used language around risk of death for the purposes of palliative care integration. 1 In addition, this terminology was appropriate as it is difficult to separate deterioration and death without the benefit of hindsight, particularly for people with a hematological malignancy who often experience episodes of sudden deterioration in their clinical condition, followed by periods of recovery and stability. 7 People who experience a clinical deterioration are likely to have physical, psychosocial, or spiritual palliative care needs, regardless of the outcome of the deterioration. For these reasons, the terms deterioration and dying have been used.
Design
Study setting, design, and population
This study performed multivariable testing of 11 clinical indicators that achieved consensus by an international expert panel of hematology and palliative care clinicians as being associated with deteriorating and dying in a three-step modified Delphi approach (conducted by the authors). 19 The Delphi approach was informed by a systematic review of prognostic factors present in the final three months of life for people with a hematological malignancy, 9 and allowed for clinician expertise to inform the research topic using rigorous methodology. 20 The mixed method and iterative approach of the Delphi was required as limited evidence existed on this complex topic. 21 The expert Delphi panel achieved consensus on 11 clinician-assessed indicators that were a combination of objective and subjective markers. 19 These were “bedside” indicators, meaning that they were easily identifiable in every day practice by a range of healthcare professionals. The expert panel felt that the use of subjective and objective indicators was appropriate to harness the expertise of clinicians and allowed for indicators to be clinically and contextually relevant, applicable in a range of different settings, and across various disease types. 19
Patients with any type of hematological malignancy were the focus of the Delphi approach, and were included in this study to optimize the application of results and not limit the scope of findings to one particular cohort. This approach is in line with other work focused on identifying risk of deteriorating and dying that includes multiple types of chronic illness. 1 The clinical indicators identified in the Delphi were developed to account for the heterogeneity in the hematological malignancies and between individuals. 19 For example, the indicators were intended to cover the different illness trajectories of someone with high-risk versus low-risk disease, the elderly population with comorbidities and treatment limitations and younger people treated aggressively who may die from infection of graft versus host disease (GVHD) post stem cell transplant. This study provides preliminary testing of the indicators for people with all types of malignancies.
A retrospective case–control design was utilized as it is an efficient and effective method of testing association between an exposure and an outcome and is widely used in prognostic study designs. 22 This was a single-center, hospital-based study of people with a hematological malignancy admitted to a large tertiary hospital (929 beds) in Brisbane, Australia. This cancer center cares for people with all types of hematological malignancies at any stage in the illness trajectory, including those receiving allogeneic stem cell transplantation.
Cases were selected historically (sequentially) as the last 120 deceased people who met the following criteria: (1) had been diagnosed with a hematological malignancy; (2) were ≥18 years of age at the time of death; and (3) experienced an in-patient admission at any time in the six months before study entry between 1st July 2014 and 31st December 2015. Dates of deaths of cases ranged from 28th November 2014 to 12th December 2016. Each case was individually matched with two living controls by gender and disease type. Controls were identified and matched historically (sequentially) using the same selection criteria and in the same time period (1st July 2014–31st December 2015). Cases were matched with controls by the closest time period of hospitalization. To ensure that there was clear distinction between cases and controls, only patients who were alive at least six months after the study inclusion time period were selected as living controls for the deceased cases. Disease type was classified using the World Health Classification of Hematological Neoplasms 23 and the Revised European–American classification of lymphoid neoplasms.24,25 Classification manual available as Supplementary Table S1 (Supplementary Data are available online at www.liebertpub.com/jpm). For cases, the date of study entry was the date of death, and for controls it was the most recent admission historically before 31st December 2015. Data were collected backwards in time.
A sample size of 360 patients was determined to be feasible in consideration of time, resources, and availability of the sample. 26 In the absence of any reliable data to perform sample size calculations, it was arbitrarily assumed that a minimum of approximately 10% of the study sample would have the presence of any particular clinical indicator to be tested as recommended in the absence of available data. 27 Based on a minimum odds ratio of 2.5, type one error set at 0.05, and case/control ratio of 1:2, a sample size of 360 was estimated to result in a power of 83.68%.
Ethical approval
Low and negligible risk ethical approval was obtained from the Royal Brisbane and Women's Hospital Human Research Ethics Committee (HREC/15/QRBW/289). A waiver of consent was approved as this retrospective study was considered to be of negligible risk.
Study variables
Data were collected on every in-patient admission for the preceding six months before study entry for cases (study entry–date of death) and controls (study entry–most recent admission historically before 31st December 2015). Demographic characteristics as well as admission and discharge dates were recorded. The outcome variable was mortality (yes/no). The explanatory variables tested were demographic characteristics and the presence of 11 clinical indicators (yes/no) as identified in a modified Delphi approach study performed by the investigators. 19 See Table 1 for description of clinical indicators tested. The clinical indicators were a combination of objective and subjective variables as is common in clinical tools that identify risk of dying. 28 As discussed in the results of the Delphi, many of the indicators overlapped slightly with each other or were interlinked, 19 for example, age and comorbidities and disease status and treatment limitations. This was deemed by the expert panel to be unavoidable and clinically appropriate due to the interconnected nature of illness. 19 Statistical tests were performed to ensure that the clinical indicators were independent variables. To assess for multicollinearity between the explanatory variables, the variables in the logistic regression model were run in a linear regression model with age (continuous) as the dependent variable.
Clinical Indicators Tested
Steps were taken to reduce risk of bias during data collection and to enhance the rigor of the approach. A data collection manual was developed with descriptions of all possible clinical scenarios and definitions around each of the variables to be collected. Notes were taken by the researcher on the characteristics associated with each clinical indicator variable for record keeping and data cleaning purposes. The data collection tool and manual were piloted on 30 patients to ensure that the study protocol was feasible and to address any unforeseen issues. As no significant changes were made to data collection, the first 30 patients were included in the final sample. Inter-rater and intrareliability testings were performed.
Study protocol
Cases and controls were identified via a list of in-patient admissions by hospital computer systems. All records were cross referenced with individual patient medical records for accuracy. Diagnoses were confirmed via pathology reports and documentation of the treating hematology medical team. Data were extracted from the electronic medical records of each patient and were recorded on a specially designed data collection tool. The presence of each clinical indicator was assessed at two time periods for each in-patient hospital admission: (1) on admission (within the first 24 hours); and (2) during admission (during the remainder of the admission), as far back as six months from study entry. The indicators age and high-care preadmission were only assessed on admission.
Inter-rater and intrarater reliability testing
Inter-rater reliability was assessed by a trained second reviewer independently assessing medical records of a random 5% of the study sample (n = 18). Intra-rater reliability testing was also performed and analysis was observed and checked by the second reviewer and another member of the research team. The researcher recollected data on a random 5% of the study sample medical records at a later point of time. The random samples were generated via the random number sequencing website www.random.org. 29 A kappa coefficient (K) ≥0.61 was considered to be acceptable.
Data analysis
A statistical analysis plan was agreed upon a priori. All analyses were performed in the statistical software package IBM SPSS Statistics version 23 (SPSS). A p-value of >0.05 was considered to indicate statistical significance. Descriptive statistics were computed for all study variables to describe the sample and check for missing data. Categorical/dichotomous variables were examined via frequencies and percentages. Central tendency and distribution was calculated for continuous data. Chi-square tests were used to assess differences between cases and controls for demographic characteristics and presence of clinical indicators.
Univariate logistic regression was used to test the effect of a single explanatory variable (demographic characteristics and presence of clinical indicators) on the outcome variable (mortality). Odd ratios (OR), 95% confidence intervals (CI), and p-values were reported to measure the association between each explanatory variable and mortality. A positive association was considered to be present when the odds ratio was greater than one. Conversely, a negative association was represented by an odds ratio of less than one. A large confidence interval reflected a low level of precision of the odds ratio, and a small confidence interval indicated a higher precision. 30 Variables with a p-value <0.2–0.05 in univariate logistic regression were included in a multivariable binomial logistic regression model. Multivariable analysis was used to identify independent predictors of mortality and accounted for the effect of multiple explanatory variables on the outcome variable. It was possible to perform multivariable logistic regression as the sample size was large enough to ensure that there were at least 10 events for each explanatory factor tested.31,32 Backwards logistic regression modeling was manually performed in SPSS. Cox and Snell R square and Nagelkerke R square were used to assess model performance. Indicators that were significant in multivariable analysis were considered to be predictive of mortality. For the clinical indicators and disease risk (low, moderate, and high), survival curves were plotted using the Kaplan–Meier method. Significance was assessed via log rank test p-values, with p < 0.05 indicating a clinical indicator affecting survival.
To identify time periods that clinical indicators were present before death, univariate and multivariable logistic regressions were computed for five different time periods in the study going backwards in time using the same method described above: (1) between 0 and 7 days before study entry; (2) between 8 and 30 days before study entry; (3) between 31 and 90 days before study entry; (4) between 91 and 180 days before study entry; and (5) at or before 181 days before study entry. Each indicator was marked (yes/no) as if it was present at any time during a hospital admission for each of these time periods. If a patient did not have an admission in a time period, this was recorded as missing data.
Results
A total of 356 patients were included in the study. Four cases with very rare or poor prognostic disease were only able to be matched with one control each. Therefore, 120 cases were matched with 236 controls. Fourteen controls (5.9% of controls) were unable to be matched by gender and were matched on disease only. No patients were excluded due to incomplete or missing records. The median age of the sample was 62 (IQR 18, 19–91). Table 2 displays the demographic characteristics of the sample. The median time a patient was followed backwards in time from study entry was 45 days (IQR = 105, 1–219 days). The median number of admissions data collected on was 1 (IQR = 1, 1–7), and the median length of stay for each admission was nine days (IQR = 18, 1–127). There was very little missing data; relationship status and place of birth each had two missing cases (0.6%) as shown in Table 2.
Comparison of Baseline Characteristics and Clinical Indicators
Univariate comparisons of variables between cases and controls were by Pearson's Chi-square.
Inter-rater reliability testing
Overall, there was an acceptable level of consistency (K > 0.61) for inter- and intrarater reliability. There was substantial agreement (K 0.61–0.80) or almost perfect agreement (K 0.81–1.0) for inter-rater reliability for nine of the clinical indicators. The indicator fungal infections achieved moderate agreement (K 0.61–0.80). Fungal infection was not a commonly found indicator; therefore, any differences between the two assessors may have caused a disproportionately large change in the K score compared with more frequently observed indicators. There was substantial or almost perfect agreement for intrarater reliability for 10 of the indicators. There were not enough observations to calculate an inter-rater reliability score for graft versus host disease (GvHD) and intrarater reliability score for fungal infection.
Clinician indicators that identify risk of deteriorating and dying
None of the demographic variables was statistically different between cases and controls in Chi-square tests. All clinical indicators were significantly different between cases and controls in Chi-square tests (Table 2). None of the demographic characteristics was significantly associated with mortality in univariate logistic regression. However, the variable relationship status achieved a p-value that met criteria to be included in multivariable logistic regression (OR 1.40, 95% CI 0.88–2.22, p = 0.154). All the clinical indicators were significantly associated with mortality in univariate logistic regression. The clinical indicators that were most predictive of death in univariate analysis were treatment limitations (OR 28.26, 95% CI 15.21–52.49, p = < 0.001) and declining performance status (OR 20.06, 95% CI 11.095–36.26, p = < 0.001). See Table 3 for results of univariate logistic regression for clinical indicators.
Results of Univariate Logistic Regression for Baseline Characteristics and Clinical Indicators
OR, odds ratio; CI, confidence interval.
In multivariable logistic regression analysis, the following variables were independent predictors of mortality if they were present at any time in the six-month study period: declining performance status; presence of two or more comorbidities; relapse, refractory, or persistent disease; persistent infections (bacterial or viral); invasive fungal infections; and treatment limitations. See Table 4 for results of multivariable analysis. The remaining variables tested were excluded from the model as they were not significant when accounting for the other factors. The model containing the six predictive clinical indicators was statistically significant (n = 356, χ2 = 245.265, df = 6, p = < 0.001). The model as a whole explained between 50% (Cox and Snell R square 0.498) to 69% (Nagelkerke R square 0.690) of the variance between living and deceased patients and correctly classified 88.5% of cases. These variables described 49%–68% of the variation between deceased cases and living controls. Collinearity statistics were computed and were satisfactory for all the included variables as they were all under the value of three. 30
Results of Multivariable Logistic Regression
Kaplan–Meier analysis using the log rank test identified that there was a statistically significant association between survival and the following: (1) declining performance (log rank test p < 0.001); (2) two or more comorbidities (log rank test p < 0.001); (3) relapsed, refractory or progressive disease (log rank test p < 0.001); (4) persistent infections (log rank test p < 0.001); (5) treatment limitations (log rank test p < 0.001; (6) high care (log rank test p < 0.001); (7) signs of frailty (log rank test p < 0.001); (8) refractory graft versus host disease (log rank test p < 0.005); and (9) weight loss (log tank test p < 0.027). There was no significant difference in the time to survival between people with and without invasive fungal infections (log rank test p = 0.574). Survival curves are displayed in the Supplementary Figures S1–S10. Kaplan–Meier analysis using the log rank test indicated that there was no difference in time to survival between people with different risk categories of disease (p = 0.063): low, moderate, and high. See Supplementary Figure S11.
Due to the relatively small sample size, we were unable to perform multivariable logistic regression modeling for each disease type or category of risk. To assess if there were differences between clinical indicators present among patients with low-, moderate-, and high-risk diseases, Chi-square tests were calculated for the six clinical indicators that were included in the multivariable logistic regression model, stratified according to disease risk. Declining performance was significantly different in the cases in the high-risk group, compared with the moderate- and low-risk groups (44.5%, 76.1%, 81.5%, respectively, χ2 = 127.40, p = < 0.001). There was no clinically significant difference among low-, moderate-, and high-risk diseases for the presence of relapsed, refractory disease status, comorbidities, and infections in the deceased cases. There were greater fungal infections present in deceased cases with high-risk disease versus those with moderate- and low-risk disease (19.1%, 8.7%, 11.1%, respectively, χ2 = 18.02, p = < 0.001). There were also fewer treatment limitations in the deceased cases with high-risk disease versus those with moderate- or low-risk disease (66%, 73.9%, 70.4%, respectively, χ2 = 151.39, p = < 0.001).
When clinical indicators are first present before death
Descriptive statistics for when the clinical indicators were first present are shown in Table 5, including if they were present for the first time (before study entry) on or during an admission, and the median time before death for cases. All the clinical indicators were more commonly identified for the first time (before study entry) on admission with the exception of persistent infections and invasive fungal infections. The indicator high-care needs preadmission was only measured on admission. When univariate and logistic regression was computed for if a clinical indicator was present in cases and controls (yes/no) in five different time periods, five clinical indicators were consistently significant in multivariate logistic regression if they were present in the 0–7 days, 8–30 days, and 31–90 days before study entry, these included the following: (1) declining performance status; (2) treatment limitations; (3) relapsed, refractory, or persistent disease; (4) persistent infections; and (5) two or more comorbidities. Fungal infection was present in the overall model (any time in six months before study entry) and one of the time-dependent models (0–7 days before study entry). As the models went further back in time from study entry, they accounted for less variance between deceased cases and living controls. All the indicators in the first three time models (0–7, 8–30, and 31–90 days before study entry) had considerably high odds ratios, indicating a good magnitude of effect on the dependent variable. Table 6 provides a summary of logistic regression models for all the time periods, including if clinical indicators were present at any time in the six months before study entry.
Time Frames When Clinical Indicators Were First Present
Time clinical indicator present for cases only.
Summary of Logistic Regression Models for Different Time Periods
Cox and Snell R squared and Nagelkerke R-squared.
N.B. For cases, the date of study entry was the date of death, and for controls, it was the most recent admission historically before 31st December 2015. Data were collected backwards in time.
The assumptions of logistic regression were met for this study. There were 11 predictors tested and 120 cases, meaning there were enough events per predictor tested.33,34 Alternatively, 15 cases are needed for each variable tested, which was also achieved in this study. 30 The sample had no outlying cases as the variables were categorical. There were no issues with collinearity between variables in any of the models. All the included variables had a variance inflation factor reading of under three. 30
Discussion
Six clinical indicators were found to independently predict death in the final three months of life for people with a hematological malignancy. Treatment limitations and declining performance status were the strongest indicators of mortality. The prognostic value of performance status is substantiated by the few studies that have explored predicting survival at the end of life for people with a hematological malignancy. Chou et al. 13 reported that the Palliative Prognostic Index (PPI) accurately predicted survival in terminally ill people with a hematological malignancy known to a specialist palliative care service (median survival of 16 days, IQR 4–4.75 days). Ohno et al. 12 reported the PPI was predictive of death in people with a hematological malignancy admitted to an in-patient hematology unit (median survival of 6.4 weeks). Although declining performance status was a strong predictor of mortality in this study, only 27.5% (n = 98) of patients had this indicator present at some time during the six months study period. Similarly, only 23.6% (n = 77) of the sample had comorbidities ≥2. This may be due to the study sample being in-patients in a large tertiary hospital. It is possible that many patients were receiving aggressive treatments and died suddenly of complications. In addition, patients with poor performance or comorbidities, who are not eligible for aggressive treatment, may remain in rural or regional areas and may be underrepresented in the study sample.
The indicator treatment limitations (of the underlying malignancy) was found to be a strong predictor of mortality and was defined as: (1) inability to tolerate treatment; (2) no appropriate curative or life-prolonging treatment available or effective; or (3) limitations set in place by the healthcare team or family. The appropriateness of these indicators seems apparent considering that disease progression is often the leading cause of death for people with a hematological malignancy. This subjective definition is that this indicator is clinically relevant in the age of targeted therapies and advances in available treatments.
The core components of the PPI include performance status, dyspnea, delirium, oral intake, and edema. 13 Ohno et al. 12 also found that a prognostic model developed by Kripp et al., 11 which included low performance status, requiring opioid analgesia, low platelet count, low albumin levels, and high lactate dehydrogenase levels, was predictive of death for people with a hematological malignancy on a palliative care unit (median survival of 58 days). With the exception of performance status, the range of prognostic factors mentioned above were not assessed in our study as clinical indicators that achieved consensus in the Delphi approach (that informed this study) were focused on earlier markers of deteriorating (ref–under review). This was to facilitate providing proactive palliative care planning before acute deterioration. Differences in our results and the Chou et al., 13 Ohno et al., 12 and Kripp et al. 11 are likely to be due to these studies having shorter time frames of survival in their populations than our study. To the best of our knowledge, the remaining indicators identified in our study have not been assessed in any other prognostic study designed to improve palliative care provision and transitioning to end-of-life care. This study significantly adds to the body of knowledge as it is the first of its kind to identify clinical indicators that signal risk of deteriorating and dying up to three months before death.
The findings of this study indicate that identifying risk of deteriorating and dying with a core set of clinical indicators is most feasible within the final three months of life for people with a hematological malignancy. This is a shorter time frame than used in most clinical tools that identify risk of deteriorating and dying within 6–12 months in other populations. 28 This is not surprising considering the rapid speed, in which people with a hematological malignancy can deteriorate. 5 While most indicators were present first on admission, persistent infections and invasive fungal infections were more commonly found for the first time during admission. Therefore, an assessment of the clinical indicators should be performed on admission and repeated at other clinically relevant time points.
This study had several limitations, including a relatively small sample size, leading to wide confidence intervals. The magnitude of the effect, statistical significance, and clinical significance of results were considered in the interpretation of findings to counteract this limitation. The study was set in a single center, thus limiting external validity. However, the hospital sample is likely to share many characteristics with other tertiary centers, increasing generalizability of results. A strength associated with collecting data at a single site is the uniformity of medical records, which enhanced data collection procedures. The retrospective nature of this study was a limitation as medical records are often not completely accurate or detailed. The researcher took due diligence to cross-reference all medical records for accuracy and completeness where possible. Prospective testing of the clinical indicators is required in a larger multi-site study, with exploration into different disease types (i.e., low-risk versus high-risk disease) to assess if there are differences in clinical indicators for different patient groups. Despite these limitations, this study has provided valuable preliminary data and adequately addressed the study aims.
Findings of this study can help health professionals internationally identify the most appropriate time to engage in honest, sensitive discussion with patients and their families regarding possible risk of deteriorating and dying. This will enable a collaborative reassessment of the goals of care and planning for potential need. Ultimately, the findings of this study will help to give patients more control over their death and the time they have remaining.
Conclusion
This study identified six clinical indicators associated with deteriorating and dying for people with a hematological malignancy in a multivariable logistic regression model. Future research is needed to test these findings further. However, this study has begun to address a gap in knowledge and informed palliative care provision and transitioning to the end of life for people with a hematological malignancy.
Footnotes
Acknowledgments
The lead author was supported by the following organizations: (1) National Health and Medical Research Council, Centre for Research Excellence in End-of-Life Care, Australia; (2) Queensland University of Technology, Australia; and (3) Royal Brisbane and Women's Hospital Foundation, Australia and the Centaur Memorial Fund for Nurses. These organizations were not involved in the conduct of the study or preparing the article.
Author Disclosure Statement
No competing financial interests exist.
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.
