Methods
Ethics consideration
This retrospective cohort study, based on nationwide health insurance claims data, explored COVID-19 in patients with ESKD undergoing RRT. We analyzed variations based on RRT type and identified the risk factors and outcomes for COVID-19 hospitalization cases using univariate and multivariate regression analyses. This study was approved by the Institutional Review Board of Ajou University Hospital (No. AJOUIRB-EX-2022-565). All research methods were conducted in accordance with the Declaration of Helsinki, as revised in 2013.
Data source
The data for patients with ESKD originated from health insurance claims data provided by the Health Insurance Review and Assessment Service (HIRA). However, COVID-19 data for the general population were sourced from the Korea Disease Control and Prevention Agency (KDCA). Data between July 2017 and June 2022 were used in this study. The HIRA claims dataset represents one of the primary sources of real-world data in Korea, encompassing medical utilization records, cost data, and other information such as hospital visits for various diagnoses and received medical procedures and covering 98% of the entire population. The national notifiable disease surveillance system from the KDCA manages and tracks confirmed patients with COVID-19 infection and their epidemiological characteristics.
Study population
The study population comprised patients with ESKD receiving RRT between July 2017 and June 2022 who were selected based on diagnostic information from the International Classification of Diseases, the 10th revision (ICD-10) codes and procedure codes in the claims data. Individuals who died before 2020, before the COVID-19 outbreak began, were excluded. During the study period, the study population was classified into three groups: HD, peritoneal dialysis (PD), and KT. Patients with ESKD who underwent two or more types of RRT were categorized based on the following criteria. 1) Patients were defined based on the RRT type they were undergoing at the time of COVID-19 infection. 2) Patients were defined based on the RRT they received during the COVID-19 pandemic. For example, patients who underwent HD before 2020 and those who underwent PD after 2020 were classified as patients who underwent PD. 3) KT recipients were categorized as patients who underwent KT regardless of whether they had previously undergone HD or PD. Patients who underwent HD or PD after KT were classified as undergoing HD or PD, respectively. 4) Patients who underwent HD and PD after the pandemic were excluded from the analysis because their RRT types could not be determined (
Fig. 1).
Case ascertainment of end-stage kidney disease and renal replacement therapy
Patients with ESKD were defined as individuals who had one or more outpatient visits with ICD-10 code ‘N18.’ Patients with ESKD on RRT were identified based on specific procedure codes, including O7020, O7021, and O9991 for HD, O7076 and O7077 for PD, and R3280 for KT. For HD and PD, patients were ascertained if they had both the relevant ICD-10 code and corresponding procedure codes for 3 months or more. KT recipients were defined as patients who underwent KT if they had an ICD-10 code and at least one recorded procedure code.
Definition of variables
To calculate the COVID-19 infection per population using the RRT modality, we divided the number of patients infected with COVID-19 by the total number of patients observed for each RRT type over the study period. COVID-19 infection is defined as patients who visited the hospital as outpatients or inpatients with the ICD-10 codes ‘B342,’ ‘B972,’ ‘U071,’ or ‘U072’ as the primary diagnosis or the first listed condition. Regarding the general population, the COVID-19 infection rate per population was determined by dividing the reported number of COVID-19 cases from KDCA by the population data provided by Statistics Korea. COVID-19 hospitalization and death rates were computed using the same method as the COVID-19–related mortality data for the general population derived from Statistics Korea’s cause-of-death statistics. COVID-19 hospitalized patients are defined as those who, among the COVID-19–infected individuals, had one or more admissions specifically due to COVID-19. COVID-19 death refers to patients within healthcare facilities whose treatment outcomes for COVID-19 are recorded as ‘death’ in the claims data, indicating fatalities specifically due to COVID-19.
When identifying the risk factors for COVID-19 hospitalization, the following variable definitions were applied: age, sex, insurance type, and diabetes mellitus (DM) status were the most recent diagnoses. The Charlson comorbidity index (CCI) score, used to denote baseline health status, was analyzed from 1 year before RRT commencement to the initiation date [
10] (
Supplementary Table 1, available online). Comorbidities were ascertained through the main diagnosis or subdiagnosis in the claims data, with dementia and human immunodeficiency virus infection remaining unconfirmed within the CCI score owing to data privacy concerns.
Statistical analysis
The characteristics of patients with ESKD and those hospitalized for COVID-19 are described using summary statistics. Differences among the RRT types were analyzed using the chi-square test and analysis of variance (ANOVA). For patients who underwent HD, the chi-square and t tests were performed based on whether they were in a nursing hospital. For both comparisons, continuous variables were analyzed using the ANOVA and t tests, whereas categorical variables were analyzed using the chi-square analysis. Univariate and multivariate logistic regression models were used to investigate risk factors associated with COVID-19 hospitalization. The initial step involved conducting univariate logistic analysis for individual covariates, utilizing odds ratios with 95% confidence intervals to assess the association between COVID-19 hospitalization and each covariate. Covariates that exhibited significant associations in the univariate analysis were incorporated into the final multivariate model, and multivariate analysis was conducted. Statistical significance was set at p < 0.05. All statistical tests were performed using SAS version 9.4 (SAS Institute).
Results
This retrospective cohort study utilized nationwide data between July 2017 and June 2022 to examine the disparities in COVID-19 among renal replacement modalities throughout the pandemic period. During the COVID-19 pandemic, 113,761 patients survived and underwent RRT. Patients who underwent HD constituted the majority at 90.7% (103,191), followed by KT at 6.5% (7,380) and PD at 2.8% (3,190). Substantial differences were observed among the three groups in terms of sex, insurance type, COVID-19 infection status, RRT duration, and comorbidities. The proportion of males among those undergoing RRT ranged from 57% to 61%, with the average age being highest in the HD group, followed by the PD and KT groups. Regarding insurance type, the proportion of medical aid recipients was higher, ranging from 9.6% to 18.1%, which is considerably greater than the 3% medical aid rate in the general population [
11]. Among patients with DM, the HD group had the highest prevalence (70.6%), and the CCI score, reflecting baseline health status, was also highest in the HD group, at 5.5 points, compared with 4.7 points in the PD and KT groups (
Table 1).
The characteristics of the patients receiving HD at nursing and general hospitals differed. Patients in nursing hospitals were 8 years older, with an average age of 71.0 years. Additionally, the proportion of females was higher among nursing hospital patients. The CCI score was also higher, and the prevalence of cerebrovascular disease and paraplegia were higher among nursing hospital patients (p < 0.001) (
Supplementary Table 2, available online).
COVID-19 infection by renal replacement therapy modality
During the COVID-19 pandemic, 1.3% (n = 1,391) of patients who underwent HD, 0.9% (n = 30) of those who underwent PD, and 2.1% (n = 153) of those who underwent KT were infected with SARS-CoV-2. The proportion of infected patients who underwent HD residing in nursing hospitals was approximately three times higher than that of other patients. COVID-19–infected individuals across all RRT types tended to have a higher average age, with no significant sex differences observed. In the HD and KT groups, the prevalence of DM and CCI scores were higher among COVID-19–infected patients. Individuals with COVID-19 infection in the PD group had lower CCI scores (
Table 2;
Supplementary Table 3, available online).
Among patients with COVID-19, those who were hospitalized for COVID-19 were older than those who were not hospitalized for COVID-19. No meaningful differences were observed in RRT duration between the two groups, but the CCI scores and DM prevalence of those who were hospitalized for COVID-19 were higher than those who did not require hospitalization. COVID-19–related hospitalizations were all single occurrences across the three groups, and the average length of hospital stay was highest among patients who underwent PD at 16.8 days, followed by those who underwent HD and KT. However, the differences were insignificant (
Table 2).
Fig. 2 illustrates the number of COVID-19 cases, COVID-19–related hospitalizations, and COVID-19–related deaths per 10,000 individuals during the study period. Week 1 marked the initial observation of COVID-19 in patients with ESKD undergoing RRT on December 6, 2020. Throughout the study, the absolute number of COVID-19 cases was highest in patients who underwent HD. However, when adjusted per 10,000 people, the order was KT, HD, and PD, with a weekly average infection rate of 2.43 for KT, 1.64 for HD, and 1.45 for PD. In contrast, the general population exhibited a significantly higher infection rate (43.64%), which was more than 20 times higher.
When considering COVID-19–related hospitalizations per 10,000 people, HD had the highest average at 1.13 patients, followed by KT and PD. For COVID-19–related deaths per 10,000 people, HD had the highest rate of 1.2 deaths, followed by PD (1.1 deaths), KT (0.2 deaths), and the general population (0.6 deaths). The rate of COVID-19 occurrence in patients undergoing RRT was relatively low; however, COVID-19–related deaths were higher than in the general population. Additionally, among patients on RRT, KT recipients had a higher incidence of COVID-19 but lower COVID-19–related hospitalization and death rates (
Fig. 2;
Supplementary Table 4, available online).
Risk factors associated with COVID-19 hospitalization
An analysis was conducted to identify risk factors associated with COVID-19 hospitalization. Univariate analysis revealed that the following variables had a one-to-one association with COVID-19 hospitalization: RRT modality, age, sex, type of hospital, RRT duration, and comorbidities, including myocardial infarction, congestive heart failure, cerebrovascular disease, DM, DM complications, and paraplegia (
Table 3).
Multivariate analysis was conducted to identify the risk factors influencing COVID-19 hospitalization when considering various variables together. When various variables were considered, compared with patients who underwent HD, KT recipients had a 1.22 times increased risk of COVID-19 hospitalization, although this increase was insignificant. Age exhibited an incremental risk of 1.018 for each 1-year increase, whereas RRT duration showed an increment of 1.127 for every 1-year increase. Being in a nursing hospital was associated with a 3.373 times higher risk. Additionally, individuals with comorbidities such as cerebrovascular disease, DM complications, and paraplegia had an increased risk of 1.237, 1.280, and 1.545, respectively. However, significant variables in the univariate analysis, including sex and comorbidities, such as myocardial infarction, congestive heart failure, and DM, were insignificant when considered alongside other variables in the multivariate analysis.
Discussion
In this study, we analyzed the incidence, hospitalization rate, and mortality rate of COVID-19 based on the RRT modality using data provided by the HIRA and investigated the factors that affect COVID-19 hospitalization. Based on our analysis, among patients undergoing RRT in South Korea, KT recipients had the highest incidence of COVID-19, followed by those undergoing HD and PD; however, all had a lower incidence than the general population. Hospitalization rates were highest among patients who underwent HD, followed by those who underwent KT and PD. Similarly, the mortality rate was the highest among patients who underwent HD, followed by those who underwent PD and KT. The KT recipients exhibited a lower mortality rate than the general population. In the multivariate analysis of risk factors for hospitalization for COVID-19, age, residence in a nursing hospital, RRT duration, history of cerebrovascular disease, DM complications, and paraplegia were shown to be related to COVID-19 hospitalization in patients with ESKD.
In our analysis, among patients receiving RRT, those undergoing HD accounted for the majority (approximately 90.7%), followed by KT (6.5%) and PD (2.8%), which is generally consistent with the statistics of the prevalence of RRT modality reported by the Korean Renal Data System [
12]. In this report, HD accounted for approximately 79.4% of RRT in South Korea, with KT (16.5%), followed by PD (4.1%). Although the ratios slightly differed, these data reflect the RRT situation in South Korea.
Among the RRT modalities, KT had the highest incidence of COVID-19, followed by HD and PD. This finding deviated from what was initially anticipated. As patients who underwent HD visit the dialysis center multiple times a week, they might have the highest risk considering the contagious nature of COVID-19. Previous studies conducted in other countries reported that patients on in-center HD showed the highest incidence of COVID-19 among RRT modalities [
13,
14]. However, in our analysis, patients who underwent HD had a lower incidence of COVID-19 than those who underwent KT. In South Korea, the rapid spread of COVID-19 has led to government regulations enforcing social distancing, which Koreans, including patients, have strictly followed. Moreover, the Korean Society of Nephrology presented practical guidelines for dialysis centers to respond to infectious diseases [
7]. These efforts may have contributed to reducing the incidence of COVID-19 in patients undergoing HD.
The incidence of COVID-19 in patients undergoing RRT was considerably lower than that in the general population. Many other studies [
15] have reported a higher incidence of COVID-19 in patients on RRT than in the general population; however, some of these studies explained that because patients on RRT frequently visit medical centers and are perceived to be a high-risk group, COVID-19 tests are more readily available to them than to the general population [
16,
17]. In South Korea, the costs of COVID-19 testing and treatment after diagnosis were covered by the government, and indications for COVID-19 testing were strictly regulated. Several COVID-19 testing stations have been established to make testing accessible. The COVID-19 response system of Korea has been regarded as excellent in other countries [
18,
19]. With this social background, not only high-risk patients on RRT but also the general population could easily undergo COVID-19 tests, which might have contributed to a significantly higher incidence of COVID-19 in the general population. In addition, making a direct comparison is challenging because the characteristics of the pandemic differed by country and region. However, it can be considered that South Korea effectively maintained a robust system for preventing COVID-19 in high-risk patients.
Among the RRT modalities, hospitalization and mortality rates were highest in HD patients. The age of patients who underwent HD was higher than those of patients who underwent PD and KT. In addition, comorbidities such as DM and DM complications were more prevalent, and the CCI scores were higher in patients who underwent HD. A greater number of patients in the HD group had a longer RRT duration. These factors may be related to higher COVID-19 hospitalization and mortality rates in patients on HD, and this result correlates with those of previous studies [
13,
14]. As relatively older patients with multiple comorbidities and those admitted to nursing hospitals tend to undergo HD rather than KT or PD, COVID-19 hospitalization and mortality rates are commonly highest in patients who underwent HD.
As KT recipients use immunosuppressive agents, infectious diseases such as COVID-19 can be fatal. The findings from our study present an interesting contrast to what might have been expected. Despite KT recipients exhibiting the highest incidence of COVID-19, the mortality rate was lower than that observed in those who underwent HD and PD and even lower than that in the general population. Previous studies have reported that mortality rates among KT recipients were higher than those in the general population [
20–
22]. However, in one study, when solid organ transplant recipients diagnosed with COVID-19 were matched with non-solid organ transplant patients diagnosed with COVID-19 based on age, race, and admission status, the odds of COVID-19-associated severe disease and intubation rates were similar between the two groups [
23]. Another study reported that KT was not associated with a higher rate of severe disease when matched to control patients with similar chronic conditions or disease severity upon admission [
24]. This means that the higher mortality reported in other studies may have been influenced by other characteristics, such as age and comorbidities, rather than the fact that they had undergone transplantation. One study that reported higher mortality in KT recipients indicates that comorbidities had a stronger prognostic effect than immunosuppression use
per se [
21]. Therefore, if immunosuppressive agents are well controlled based on infection severity, KT may not be a risk factor for COVID-19 hospitalization. Moreover, given the low mortality rate among KT recipients in our analysis, it appears that the efforts of patients and medical staff, such as quick diagnosis and monitoring when symptoms occurred, were successful.
When we analyzed the risk factors associated with COVID-19 hospitalization, age, residence in a nursing hospital, RRT duration, history of cerebrovascular disease, DM complications, and paraplegia were found to be related to COVID-19 hospitalization. Age is a well-known risk factor for COVID-19–related mortality in patients undergoing dialysis and the general population [
25,
26]. Considering that patients in nursing hospitals tend to be older and frailer than those receiving HD in general hospitals, coupled with the fact that COVID-19 can spread rapidly within nursing facilities, residing in a nursing hospital can be a risk factor for COVID-19 hospitalization [
27].
In the multivariate analysis of risk factors for COVID-19 hospitalization, although the RRT modality was related to hospitalization due to COVID-19 in the univariate analysis, it was unrelated to hospitalization when multiple variables were considered together. This result correlates with our previous discussion, suggesting that poor prognosis, such as the mortality of KT recipients, can be related to patients’ underlying conditions, such as age or comorbidities, rather than the transplantation. Factors other than the RRT modality might have a greater effect on COVID-19 hospitalization.
To our knowledge, this study is the first Korean study to compare the incidence, hospitalization, and mortality of COVID-19 based on RRT modalities and analyze the risk factors for COVID-19 hospitalization using health insurance claims data provided by HIRA. However, this study has some limitations. First, patients who could not be traced in claims data owing to data privacy concerns were excluded from the study. Second, patients treated in hospitals with fewer than three COVID-19 cases during the study were excluded from the provided statistical data. Additionally, patients who were infected with COVID-19 between January and March 2020 were excluded from the research [
28]. However, between January and March 2020, the early stages of the COVID-19 outbreak, the number of COVID-19 cases was extremely low. Although there may have been variations in patient numbers, the patients included in this study adequately represented the overall population of patients who underwent RRT. Third, information on the vaccination status of the patients was unavailable; therefore, whether the patients had been vaccinated and, if so, how many times were not considered. Fourth, we did not include propensity score matching on baseline characteristics such as age and sex. Since COVID-19 data for the general population were sourced from KDCA, we could not access the specific information of general population to perform propensity score matching. However, this study aimed to compare the characteristics of COVID-19 infection according to the RRT modality and analyze risk factors of COVID-19 hospitalization. Fifth, in this study, we could not analyze risk factors for COVID-19 deaths due to the small number of deaths reported during the study period. Further study about COVID-19 death would be helpful. Finally, as isolation guidelines due to the COVID-19 outbreak changed over time, the characteristics of hospitalized patients may have changed accordingly. Hospitalization usually indicates a serious infection, but this may be slightly different in the early stages of COVID-19 pandemic as patients with ESKD sometimes needed hospitalization for the purpose of isolation even if they did not show serious symptoms of COVID-19 [
29]. However, because the observation period was relatively short, it was not possible to split the data by period according to isolation guidelines.
In conclusion, among the RRT modalities, KT showed the highest incidence of COVID-19; however, hospitalization and mortality rates were the highest in HD. Several factors, such as age, residence at a nursing hospital, RRT duration, history of cerebrovascular disease, DM complications, and paraplegia, were identified as risk factors for COVID-19 hospitalization in patients receiving RRT. As COVID-19 has not been completely resolved and infectious diseases such as this may appear in the future, this analysis will help to understand the characteristics of infectious diseases in patients on RRT and prepare for other infectious diseases that may occur in the future.