- 1Department of Internal Medicine, University Hospital Maastricht, Maastricht, Netherlands
- 2Fresenius Medical Care Italia SpA, Palazzo Pignano, Italy
- 3Fresenius Medical Care Deutschland GmbH, Bad Homburg, Germany
- 4Renal Research Institute, New York, NY, United States
- 5Icahn School of Medicine at Mount Sinai, New York, NY, United States
- 6Fresenius Medical Care, Waltham, MA, United States
Introduction: Patients with end-stage kidney disease face a higher risk of severe outcomes from SARS-CoV-2 infection. Moreover, it is not well known to what extent potentially modifiable risk factors contribute to mortality risk. In this historical cohort study, we investigated the incidence and risk factors for 30-day mortality among hemodialysis patients with SARS-CoV-2 infection treated in the European Fresenius Medical Care NephroCare network using conventional and machine learning techniques.
Methods: We included adult hemodialysis patients with the first documented SARS-CoV-2 infection between February 1, 2020, and March 31, 2021, registered in the clinical database. The index date for the analysis was the first SARS-CoV-2 suspicion date. Patients were followed for up to 30 days until April 30, 2021. Demographics, comorbidities, and various modifiable risk factors, expressed as continuous parameters and as key performance indicators (KPIs), were considered to tap multiple dimensions including hemodynamic control, nutritional state, and mineral metabolism in the 6 months before the index date. We used logistic regression (LR) and XGBoost models to assess risk factors for 30-day mortality.
Results: We included 9,211 patients (age 65.4 ± 13.7 years, dialysis vintage 4.2 ± 3.7 years) eligible for the study. The 30-day mortality rate was 20.8%. In LR models, several potentially modifiable factors were associated with higher mortality: body mass index (BMI) 30–40 kg/m2 (OR: 1.28, CI: 1.10–1.50), single-pool Kt/V (OR off-target vs on-target: 1.19, CI: 1.02–1.38), overhydration (OR: 1.15, CI: 1.01–1.32), and both low (<2.5 mg/dl) and high (≥5.5 mg/dl) serum phosphate levels (OR: 1.52, CI: 1.07–2.16 and OR: 1.17, CI: 1.01–1.35). On-line hemodiafiltration was protective in the model using KPIs (OR: 0.86, CI: 0.76–0.97). SHapley Additive exPlanations analysis in XGBoost models shows a high influence on prediction for several modifiable factors as well, including inflammatory parameters, high BMI, and fluid overload. In both LR and XGBoost models, age, gender, and comorbidities were strongly associated with mortality.
Conclusion: Both conventional and machine learning techniques showed that KPIs and modifiable risk factors in different dimensions ascertained 6 months before the COVID-19 suspicion date were associated with 30-day COVID-19-related mortality. Our results suggest that adequate dialysis and achieving KPI targets remain of major importance during the COVID-19 pandemic as well.
Introduction
COVID-19 has greatly affected the patient population with end-stage kidney disease (ESKD). Patients on hemodialysis (HD) are more vulnerable to infection due to treatment with other patients in a confined space and shared transportation (1) and also at higher risk for adverse outcomes (2). Identifying patients at risk for adverse outcomes might aid individualized treatment decisions such as the need for hospitalization and more intensive monitoring. Even more importantly, the identification of modifiable risk factors for adverse outcomes might reduce the future mortality risk of COVID-19 in this population. Previous outcome studies in patients with Dialysis-Dependent Chronic Kidney Disease (CKD5D) identified predictive factors at the time of diagnosis (3, 4) when the disease might already have had a major impact on the parameters (5).
In dialysis therapy, various potentially modifiable risk factors have been identified that can also be expressed as key performance indicators (KPIs) for quality improvement purposes (6). To identify whether potentially modifiable risk factors or KPIs in the months before the COVID-19 suspicion date relate to outcome would be of significant clinical interest, as they could provide relevant intervention targets during a pandemic phase. In addition to conventional statistical models, machine learning (ML) is increasingly used in outcome prediction and diagnosis of COVID-19. In a recent study, we used a gradient descent (XGBoost) ML model for the early identification of patients with SARS-CoV-2 infection (7). The aim of the present study is to assess the association between potentially modifiable risk factors and KPIs, obtained before the onset of SARS-CoV-2 infection and 30-day COVID-19-related mortality, in HD patients. For this purpose, we applied both conventional and ML techniques.
Materials and Methods
Study Design, Consent, Patients, and Setting
A historical cohort study was conducted, evaluating 30-day mortality among HD patients with a documented SARS-CoV-2 infection (Figure 1). All adult patients treated in the European Fresenius Medical Care (FMC) NephroCare network from February 1, 2020, and March 31, 2021, who consented to their pseudo-anonymized data to be used for secondary data analysis, were screened for eligibility. All patients from countries that did not systematically report COVID-19 cases in the European Clinical Database (EuCliD®) (8, 9), namely, the United Kingdom, Ireland, Israel, Lebanon, and Switzerland, were excluded. COVID-19 cases were defined by the following priority algorithm:
1. SARS-CoV-2 infection confirmed by a documented PCR test. No restriction based on the PCR essay method adopted was imposed.
2. ICD10 codes suggestive of COVID-19 documented in the morbidity report: U07.1 and U07.2.
3. ICD10 codes suggestive of COVID-19 documented in the death reason report: U07.1 and U07.2.
For all patients with a documented SARS-CoV-2 infection, the index date for the study was the first suspicion date or the first occurrence of a suggestive ICD10 code. Patients with suspected SARS-CoV-2 infection and negative PCR assay were excluded from the study. In order to have a fully known health profile at the index date, the following were excluded: patients without renal replacement treatment onset date available, incident patients or patients who started dialysis more than 90 days before FMC registration, patients with no documented dialysis sessions in the 14 days before the index date, and patients ineligible for KPIs in the calendar month preceding the index date.
Endpoint Definition
The primary endpoint was the occurrence of death within 30 days after the index date. Patients were considered alive at the end of the follow-up period only if recovered or with information about death or dialysis sessions after the end of the follow-up period.
Measures
Demographic variables, anthropometric variables, lifestyle variables, comorbidities, COVID-19 vaccination status (at least one dose), treatment-related factors, and blood biomarkers in the 6-month ascertainment period (Figure 1) were assessed. Comorbidities were defined as any occurrence prior to the index date of suggestive ICD10 codes as described in Supplementary Table S1. Blood biomarkers and treatment-related variables were defined as the average in the 6-month ascertainment period. The achievement of key performance healthcare indicators (KPIs) was considered in the calendar month prior to the index date according to the operative definitions reported in Supplementary Table S2. Medication use at the index date was ascertained as the occurrence of suggestive ATC code as described in Supplementary Table S3.
Data Management
A multivariate imputation by fully conditional specification regression method was applied on 31 continuous variables (SAS MI procedure) after a data-cleansing step considering missing any data that lie outside the listed upper or lower values as described in Supplementary Table S4.
Statistical Analysis
The incidence density and 95% CI of deaths in the study population based on the Poisson distribution were calculated. Risk factors for death after the COVID-19 suspicion date were evaluated with logistic regression (LR) models. A two-sided p-value <0.05 was considered statistically significant. The XGBoost (10) was also used to train and test a prediction model assessing the risk of death in the same cohort, an iterative method where, at each iteration, a new sub-model is added to correct the prediction error of the previous iteration. Each sub-model is an ensemble of decision trees. A decision tree can be roughly described as a flowchart-like structure in which each internal node represents a “discrimination test” on a given attribute (e.g., any clinical parameter or demographic characteristics); each branch of the decision tree represents the result of the discrimination test (i.e., passed/not passed), and each leaf node represents the probability of the outcome. This probability represents the prevalence of events occurring in each leaf in the training set. The iterative process ends in accordance with a pre-specified stopping rule (e.g., the maximum number of iterations or minimal acceptable average prediction error). The structure of the model is computed as a function optimization process combining the minimization of both training error and model complexity. XGBoost was selected since it is characterized by a good prediction accuracy in a broad variety of problems coupled with a short computational time. Furthermore, SHapley Additive exPlanations (SHAP) analysis enables intuitive model interpretation through an accurate and efficient estimation of the contribution of each input variable to the risk. To assess overfitting, the XGBoost model was trained in a dataset partition corresponding to 70% of records in the initial sample and tested in the remaining 30% of records. The discrimination of both the LR and the XGBoost model was assessed by calculating the area under the curve (AUC) of the receiver operating characteristics curve.
Both LR and XGBoost models were created for modifiable risk factors, as continuous parameters, and as indicator variables denoting target achievement according to KPIs as defined in Supplementary Table S2. In total, 4 models were created to assess whether the significance of modifiable factors was confirmed in both ML and LR models.
LR models were estimated using SAS software, Version 9.4 of the SAS System for Microsoft Windows. Copyright © 2013 SAS Institute Inc. SAS and all other SAS Institute Inc. products or service names are registered trademarks or trademarks of SAS Institute Inc., Cary, NC, USA. For XGBoost, the available open-source package for Python, Version 3.7.4 (Python Software Foundation, Delaware, DE, USA) was used.
Results
Sample Characteristics
We observed 15,321 HD patients who had a documented SARS-CoV-2 infection between February 1, 2020, and March 31, 2021, in 588 dialysis centers located in 23 countries. Among these patients, 6,110 were excluded because they did not meet the inclusion criteria for the study. The exclusion process was gradual: patients could have more than one reason for exclusion. Nine were younger than 18 years at the index date, 16,78 did not consent for their data to be used for scientific research, 316 did not have renal replacement treatment onset date available, 2,832 started dialysis more than 90 days before FMC registration, 585 had no HD treatments in the 14 days prior to the index date, 5 did not have death date available, 35 were excluded because they were lost to follow-up, and 650 were incident patients or ineligible for KPIs in the calendar month preceding the index date. Therefore, 9,211 have been included in the study cohort coming from 552 dialysis centers located in 22 countries (Figure 2). The mean age was 65.4 ± 13.7 years, and the dialysis vintage was 4.2 ± 3.7 years. Patients’ characteristics at the index date are reported in Table 1.
Figure 2 Study population and cohort enrollment process. The exclusion process was gradual: patients could have more than one reason for exclusion.
Table 1 Characteristics of hemodialysis patients with COVID-19 disease at the time of suspicion date.
Incidence of Death
We observed 1,912 deaths (20.8%, 95% CI: 19.9%–21.6%) within 30 days since the suspicion date among 9,211 patients with a confirmed diagnosis of COVID-19. The mortality rate was 238.2/1,000 person-months (95% CI: 227.7–249.1/1,000 person-months) over 8,027 months of follow-up overall.
Non-Modifiable and Potentially Modifiable Factors Associated With Death
To assess the association between potentially modifiable factors during the ascertainment period and COVID-19-related death, we conducted an analysis by including KPIs as defined by the Medical Patient Review—Continuous Quality Improvement Program adopted in the FMC-NephroCare network. The LR model and the XGBoost model achieved comparable discrimination performance (LR: AUC = 0.72; XGBoost: AUC = 0.69). The LR showed that target achievement for serum albumin, serum sodium, erythropoietin resistance index, hydration status, hemodynamic status, Kt/V, inflammation, and serum phosphate were all associated with lower mortality risk (Table 2).
Table 2 Factors associated with 30-day COVID-19-related mortality including key performance indicators (KPIs) as defined by the Medical Patient Review—Continuous Quality Improvement Program adopted in the FMC-NephroCare network (Supplementary Table S2).
In the second model specification, we included biochemical markers as continuous variables. The LR model and the XGBoost model achieved comparable discrimination performance (LR: AUC = 0.71; XGBoost: AUC = 0.69). According to LR, non-modifiable factors associated with mortality risk were age, male sex, dialysis vintage and history of diabetes, infectious diseases, lung disease, and genitourinary diseases; among potentially modifiable factors, smoking habit, obesity, altered markers of inflammation, off target instead of inappropriate erythopoiesis-stimulating agents (ESA) administration, hyperphosphatemia, and overhydration were associated with increased mortality risk (Table 3).
Table 3 Factors associated with 30-day COVID-19-related mortality including biochemical markers as continuous variables.
Whereas the primary analysis focusing on KPIs found that on-line hemodiafiltration was associated with lower mortality risk compared to conventional HD, the magnitude of this effect was similar but non-significant at the 0.05 level with the secondary model.
A similar pattern of association with LR models was observed with XGBoost through SHAP value inspection (Figures 3, Figure 4). However, there were some notable differences in the relative importance of predictors between LR and XGBoost models. As an example, the presence of low serum phosphate levels was an important predictor of outcome in the LR model but did not emerge as a top 10 predictor in the XGBoost model.
Figure 3 XGBoost SHapley Additive exPlanations (SHAP) related to 30-day COVID-19-related mortality including key performance indicators as defined by the Medical Patient Review—Continuous Quality Improvement Program adopted in the FMC-NephroCare network (Supplementary Table S2). WBC (Principal Component), white blood cell factor was calculated with a principal component analysis based on platelets, lymphocytes, and neutrophils, considering the first factor as explaining 47.2% of the total variance, basically a contrast of neutrophils (factor loading: 0.82865) against lymphocytes (−0.81652) with a very small loading on platelets (0.24885). IDWG, interdialytic weight gain. SHAP plots show relative feature importance for the 15 most important variables: the blue bar represents overall SHAP values and is interpreted as the relative importance of each variable to risk estimates. On the right side, SHAP values show the direction of association between predictor and risk estimates. Each dot represents one individual patient. Higher values of the predictors are represented in red color; lower values of the predictors are represented in blue color. The x-axis represents the impact of variables on risk in terms of SHAP values. Red color in correspondence with positive values suggests direct correlations between risk factors and 30-day COVID-19-related death, while red color in the region of negative SHAP values suggests inverse correlation.
Figure 4 XGBoost SHapley Additive exPlanations (SHAP) related to 30-day COVID-19-related mortality including biochemical markers as continuous variables. BMI, body mass index; OH/ECW, overhydration normalized for extracellular water. WBC (Principal Component), white blood cell factor was calculated with a principal component analysis based on platelets, lymphocytes, and neutrophils, considering the first factor as explaining 47.2% of the total variance, basically a contrast of neutrophils (factor loading: 0.82865) against lymphocytes (−0.81652) with a very small loading on platelets (0.24885). nPCR, normalized protein catabolic rate. SHAP plots show relative feature importance for the 15 most important variables: the blue bar represents overall SHAP values and is interpreted as the relative importance of each variable to risk estimates. On the right side, SHAP values show the direction of association between predictor and risk estimates. Each dot represents one individual patient. Higher values of the predictors are represented in red color; lower values of the predictors are represented in blue color. The x-axis represents the impact of variables on risk in terms of SHAP values. Red color in correspondence with positive values suggests direct correlations between risk factors and 30-day COVID-19-related death, while red color in the region of negative SHAP values suggests inverse correlation.
In all models, the inclusion of COVID-19 vaccination status did not change parameter estimates for both modifiable and unmodifiable factors.
Discussion
The present study is a large international cohort, applying both conventional statistical and ML models, confirming the importance of demographic factors, dialysis vintage, cause of renal failure, and comorbidities in association with 30-day COVID-related mortality. Even more importantly, the findings stressed the importance of potentially modifiable factors in COVID-19-related outcomes in patients with HD. In contrast to other studies, we did not collect these parameters at the onset of COVID-19 diagnosis but in the 6 months preceding the suspicion date, thus at a time when intervention theoretically might have been possible. The association with outcome was shown for risk factors in different domains and held true when potentially modifiable risk factors were analyzed as both KPIs and continuous parameters. We have previously shown that optimization of medical KPIs leads to improved survival among dialysis patients (11). In the present study, we provide evidence that this observation may hold true also during the COVID-19 pandemic.
Regarding potentially modifiable risk factors, all models showed that various parameters in the nutritional and inflammatory domains were important predictors of outcome. This was clearly expressed in the results of the XGBoost model, in which inflammatory and nutritional parameters had a strong contribution to the model. Markers of malnutrition, in which normalized protein catabolic rate (nPCR) and serum albumin were assessed during the ascertainment period, were inversely related to the study outcome. This is consistent with results from the general population in which malnutrition was also shown to be a major factor predicting the severity of COVID-19 (12). In contrast, both body mass index (BMI) above 30 kg/m2 and age-adjusted fat tissue index above the 90th percentile assessed by bioimpedance spectroscopy (BIS) were related to increased mortality risk. No association between lean tissue index and the outcome was observed. This contrasts with studies in HD patients without COVID-19 in which a low lean tissue index was related to increased mortality, especially in combination with a low fat tissue index (13). However, our study confirms the results of the ERACODA study, in which a BMI > 35 kg/m2 was associated with increased mortality (3). Thus, as also shown for the general population, obesity is also a risk factor for adverse outcomes of COVID-19 in patients with HD possibly associated with increased ACE2 receptors and disrupted adipokine expression, which mediates dysregulated inflammation in visceral fat adipocytes (14, 15).
A major risk factor for adverse outcomes was the presence of inflammation in the period before the COVID-19 suspicion date when expressed as both an increased neutrophil-to-lymphocyte ratio (e.g., white blood cell (WBC) factor in our analysis) and C-reactive protein, a finding consistent with previous studies of COVID-19-related mortality in the general population (16). Whereas the mechanism behind preexistent inflammation and COVID-19-related outcome cannot be deduced from the present study, it is tempting to speculate that immunosenescence and inflammaging, which are associated with chronic inflammation, might have played a role. It has been suggested that immunosenescence may both impair the initial defense against SARS-CoV-2 infection and play a role in the hyperinflammatory syndrome after COVID-19 (17).
Fluid overload, which is at least partially amenable to intervention, was also predictive of COVID-19-related outcomes. In our cohort, we had the availability of BIS measurements applied on a routine clinical basis. Fluid overload assessed by BIS is a consistent risk factor for outcomes in the general dialysis population (18). Whether the association with COVID-19 outcome is due to a higher risk for pulmonary complications due to an additive risk for pulmonary congestion (19) or to the association between fluid overload, malnutrition, and inflammation (20) cannot be assessed from the present study. Still, fluid overload turned out to be a major risk factor in the XGBoost model, with a major contribution to the model. Regarding other cardiovascular parameters, pre-dialysis systolic blood pressure (pre-HD SBP) was inversely related to mortality, albeit not significantly. The importance of low pre-HD SBP on short-term outcomes has been shown in previous studies in patients with HD and is likely for a major part related to underlying disease (21, 22). In contrast to the general population, high blood pressure was not related to increased mortality. Interestingly, the use of lipid-lowering agents was associated with reduced mortality. This agrees with a large observational study in the general population in Sweden (23), although the observed associations do not imply causation. Very important is the strong association between smoking and 30-day COVID-19-related mortality, which is also in agreement with studies on the general population.
Regarding parameters of mineral metabolism, our results showed that both increased and low serum phosphate levels were associated with increased mortality. We suggest that the effect of low serum phosphate levels (defined as levels between 2.5 mg/dl) on COVID-19-related outcomes might be explained by its association with malnutrition. We previously found that in combination with other markers of malnutrition, low serum phosphate levels were predictive of outcome in a large cohort of patients on HD (24).
Lastly, also dialysis prescription appears to be a relevant parameter related to outcome. This held true for single pool Kt/V when expressed as KPI, in which off-target values were related to increased mortality, as well as for on-line hemodiafiltration, which was protective in the model based on KPIs, although it did not reach significance in the model based on continuous parameters.
Our study differs from previous papers on this subject because we applied both conventional and ML models. The advantage of ML models such as XGBoost used in the present study is that they can, in contrast to more traditional models, handle both non-linear relationships and missing data. XGBoost is based on weak prediction trees, able to capture complex relationships in data without the need for high-order interactions or non-linear functions (25). Whereas with LR the association with the outcome is expressed by the odds ratio, in XGBoost, SHapley values are used for this purpose. SHAP values show the relative contribution of the different parameters to the model and thus yield additional information. In other studies, XGBoost outperformed conventional statistical techniques, e.g., in the prediction of volume responsiveness in patients with AKI (25). However, the discriminative value (expressed by AUC) of ML and LR models did not differ materially in our study, suggesting that no major cross-interaction or non-linear effects could improve prediction accuracy with the input variables considered. It may be argued that the AUC levels of the models were relatively similar, but it should be noted that the aim of the study was not to construct an optimal prediction model for COVID-19 outcomes but to assess the importance of factors in the period preceding COVID-19 suspicion date. Likely, stronger predictive models might be constructed using parameters at the time of the COVID-19 suspicion date, including disease-related factors. The differences in the relative importance of the individual parameters between the corresponding LR and ML models are explained by the intrinsic distinction in the two methodological techniques, both in the construction of the model and in the calculation of the importance of each variable.
Our study also adds to the literature because it is based on a large international cohort of dialysis providers, where parameters and KPIs are collected on a routine basis for quality assurance purposes. Previous studies based on large registries did not use detailed data on laboratory variables and body composition, as these are not usually collected. Drawbacks are the retrospective nature and the absence of data on frailty that were shown to be an important predictor of outcome in the ERACODA study (3). In addition, data were obtained mainly before the availability of vaccination; therefore, we were not able to evaluate the effect of potential cross-interaction between exposure to mRNA or viral carrier vaccine and KPI target achievement on COVID-19-related mortality risk. Whereas our results might have clinical implications in the sense that patients in whom KPIs were attained had a better outcome, this does not necessarily provide that targeted interventions in the time before SARS-CoV-2 infection might have improved outcome. However, in our opinion, our results stress the importance of adequate dialysis and supportive treatment also in the time of a pandemic.
In conclusion, this study showed, by using both conventional and ML models, the importance of modifiable risk factors, collected during the months preceding the COVID-19 suspicion date, in predicting COVID-19-related outcomes in patients on HD. Also, during the time of a pandemic, adequate dialysis and supportive treatment need to be maintained to the maximum extent possible.
Data Availability Statement
The datasets presented in this article are not readily available because they include protected health information. Specific, well-motivated, requests to access the pseudo-anonymized analytic dataset may be considered by the authors. Requests to access the datasets should be directed to LN, Luca.Neri@fmc-ag.com.
Ethics Statement
The studies involving human participants were reviewed and approved by Maastricht UMC METC azM/UM P. Debyelaan 25 Postbus 5800 6202 AZ Maastricht; Protocol number METC 2021-3077. The patients/participants provided their written informed consent for secondary analysis of their data.
Author Contributions
Conceptualization: LN and JK. Data curation: PC. Formal analysis: PC. Supervision and validation: LN, JK, FB, and SS. Writing—original draft: JK, PC, and LN. Interpretation of results: all authors. Writing—review and editing: all authors. All authors have read and agreed to the published version of the manuscript.
Conflict of Interest
PC, VK, OA, AW, YZ, FB, LU, JL, SS, and LN are employees of Fresenius Medical Care. PK and HZ are employees of the Renal Research Institute, a wholly owned subsidiary of Fresenius Medical Care. LU and PK have share options/ownership in Fresenius Medical Care. PK, HZ, LU, and JL are inventors on patents in the field of dialysis. JL is a guest editor on the Editorial Board of Frontiers in Physiology. PK receives an honorarium from HSTalks and is on the Editorial Boards of journals: Frontiers in Nephrology, Blood Purification, and Kidney and Blood Pressure Research.
The remaining author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fneph.2022.907959/full#supplementary-material
References
1. Chazot C, Weis L, Hebibi H, Jean G, Deleuze S, Levannier M, et al. Impact of First Wave COVID-19 Crisis on Dialysis Parameters of COVID-Free Hemodialysis Patients: A NephoCare France Longitudinal Retrospective Cohort Study. Blood Purif (2021) 51(6):531–9. doi: 10.1159/000517493
2. Chung EYM, Palmer SC, Natale P, Krishnan A, Cooper TE, Saglimbene VM, et al. Incidence and Outcomes of COVID-19 in People With CKD: A Systematic Review and Meta-Analysis. Am J Kidney Dis (2021) 78(6):804–15. doi: 10.1053/j.ajkd.2021.07.003
3. Hilbrands LB, Duivenvoorden R, Vart P, Franssen CFM, Hemmelder MH, Jager KJ, et al. COVID-19-Related Mortality in Kidney Transplant and Dialysis Patients: Results of the ERACODA Collaboration. Nephrol Dial Transplant (2020) 35(11):1973–83. doi: 10.1093/ndt/gfaa261
4. Kooman JP, van der Sande FM. COVID-19 in ESRD and Acute Kidney Injury. Blood Purif (2021) 50(5):610–20. doi: 10.1159/000513214
5. Chaudhuri S, Lasky R, Jiao Y, Larkin J, Monaghan C, Winter A, et al. Trajectories of Clinical and Laboratory Characteristics Associated With COVID-19 in Hemodialysis Patients by Survival. Hemodial Int (2022) 26(1):94–107. doi: 10.1111/hdi.12977
6. Toussaint ND, McMahon LP, Dowling G, Soding J, Safe M, Knight R, et al. Implementation of Renal Key Performance Indicators: Promoting Improved Clinical Practice. Nephrology (2015) 20(3):184–93. doi: 10.1111/nep.12366
7. Monaghan CK, Larkin JW, Chaudhuri S, Han H, Jiao Y, Bermudez KM, et al. Machine Learning for Prediction of Patients on Hemodialysis With an Undetected SARS-CoV-2 Infection. Kidney360 (2021) 2(3):456–68. doi: 10.34067/kid.0003802020
8. Steil H, Amato C, Carioni C, Kirchgessner J, Marcelli D, Mitteregger A, et al. EuCliD® – A Medical Registry. Methods Inf Med (2004) 43(01):83–8. doi: 10.1055/s-0038-1633841
9. Merello Godino JI, Rentero R, Orlandini G, Marcelli D, Ronco C. Results From EuCliD® (European Clinical Dialysis Database): Impact of Shifting Treatment Modality. Int J Artif Organs (2002) 25(11):1049–60. doi: 10.1177/039139880202501103
10. Chen T, Guestrin C. XGBoost: A Scalable Tree Boosting System. Proc 22nd ACM SIGKDD Int Conf Knowl Discovery Data Min (2016), 785–94. doi: 10.1145/2939672.2939785
11. Garbelli M, Ion Titapiccolo J, Bellocchio F, Stuard S, Brancaccio D, Neri L. Prolonged Patient Survival After Implementation of a Continuous Quality Improvement Programme Empowered by Digital Transformation in a Large Dialysis Network. Nephrol Dial Transplant (2022) 37(3):469–76. doi: 10.1093/ndt/gfab160
12. Kurtz A, Grant K, Marano R, Arrieta A, Grant K, Feaster W, et al. Long-Term Effects of Malnutrition on Severity of COVID-19. Sci Rep (2021) 11(1):14974. doi: 10.1038/s41598-021-94138-z
13. Marcelli D, Usvyat LA, Kotanko P, Bayh I, Canaud B, Etter M, et al. Body Composition and Survival in Dialysis Patients: Results From an International Cohort Study. Clin J Am Soc Nephrol (2015) 10(7):1192–200. doi: 10.2215/CJN.08550814
14. Al Heialy S, Hachim MY, Senok A, Gaudet M, Abou Tayoun A, Hamoudi R, et al. Regulation of Angiotensin- Converting Enzyme 2 in Obesity: Implications for COVID-19. Front Physiol (2020), 11. doi: 10.3389/fphys.2020.555039
15. Gammone MA, D’Orazio N. Review: Obesity and COVID-19: A Detrimental Intersection. Front Endocrinol (Lausanne) (2021) 12:652639. doi: 10.3389/fendo.2021.652639
16. King AH, Mehkri O, Rajendram P, Wang X, Vachharajani V, Duggal A. A High Neutrophil-Lymphocyte Ratio Is Associated With Increased Morbidity and Mortality in Patients With Coronavirus Disease 2019. Crit Care Explor (2021) 3(5):e0444. doi: 10.1097/cce.0000000000000444
17. Müller L, Di Benedetto S. How Immunosenescence and Inflammaging may Contribute to Hyperinflammatory Syndrome in Covid-19. Int J Mol Sci (2021) 22(22):12539. doi: 10.3390/ijms222212539
18. Dekker MJE, Marcelli D, Canaud BJ, Carioni P, Wang Y, Grassmann A, et al. Impact of Fluid Status and Inflammation and Their Interaction on Survival: A Study in an International Hemodialysis Patient Cohort. Kidney Int (2017) 91(5):1214–23. doi: 10.1016/j.kint.2016.12.008
19. Cui X, Chen W, Zhou H, Gong Y, Zhu B, Lv X, et al. Pulmonary Edema in COVID-19 Patients: Mechanisms and Treatment Potential. Front Pharmacol. doi: 10.3389/fphar.2021.664349
20. Dekker MJE, Konings C, Canaud B, van der Sande FM, Stuard S, Raimann JG, et al. Interactions Between Malnutrition, Inflammation, and Fluid Overload and Their Associations With Survival in Prevalent Hemodialysis Patients. J Ren Nutr (2018) 28(6):435–44. doi: 10.1053/j.jrn.2018.06.005
21. Dekker M, Konings C, Canaud B, Carioni P, Guinsburg A, Madero M, et al. Pre-Dialysis Fluid Status, Pre-Dialysis Systolic Blood Pressure and Outcome in Prevalent Haemodialysis Patients: Results of an International Cohort Study on Behalf of the MONDO Initiative. Nephrol Dial Transplant (2018) 33(11):2027–34. doi: 10.1093/ndt/gfy095
22. Maddux D, Usvyat LA, Xu D, Wang Y, Kotanko P, van der Sande FM, et al. The Association of Weekly Pre-Hemodialysis Systolic Blood Pressure and Following Week Mortality. Kidney Blood Press Res (2018) 43(1):88–97. doi: 10.1159/000487105
23. Bergqvist R, Ahlqvist VH, Lundberg M, Hergens MP, Sundström J, Bell M, et al. HMG-CoA Reductase Inhibitors and COVID-19 Mortality in Stockholm, Sweden: A Registry-Based Cohort Study. Jia W Editor PloS Med (2021) 18(10):e1003820. doi: 10.1371/journal.pmed.1003820
24. Ye X, Kooman JP, van der Sande FM, Usvyat LA, Wang Y, Maddux FW, et al. Relationship Between Serum Phosphate Levels and Survival in Chronic Hemodialysis Patients: Interactions With Age, Malnutrition and Inflammation. Clin Kidney J [Internet] (2021) 14(1):348–57. doi: 10.1093/ckj/sfz143
Keywords: SARS-CoV-2, COVID-19, hemodialysis, chronic kidney disease, end-stage kidney disease, mortality, outcomes
Citation: Kooman JP, Carioni P, Kovarova V, Arkossy O, Winter A, Zhang Y, Bellocchio F, Kotanko P, Zhang H, Usvyat L, Larkin J, Stuard S and Neri L (2022) Modifiable Risk Factors Are Important Predictors of COVID-19-Related Mortality in Patients on Hemodialysis. Front. Nephrol. 2:907959. doi: 10.3389/fneph.2022.907959
Received: 30 March 2022; Accepted: 01 June 2022;
Published: 20 July 2022.
Edited by:
Michal Nowicki, Medical University of Lodz, PolandReviewed by:
Alexandru Burlacu, Grigore T. Popa University of Medicine and Pharmacy, RomaniaJean-Sébastien Rachoin, Cooper University Hospital, United States
Copyright © 2022 Kooman, Carioni, Kovarova, Arkossy, Winter, Zhang, Bellocchio, Kotanko, Zhang, Usvyat, Larkin, Stuard and Neri. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Luca Neri, Luca.Neri@fmc-ag.com
†These authors have contributed equally to this work and share first authorship