
==== Front
Can J Kidney Health Dis
Can J Kidney Health Dis
CJK
spcjk
Canadian Journal of Kidney Health and Disease
2054-3581
SAGE Publications Sage CA: Los Angeles, CA

10.1177/20543581241280698
10.1177_20543581241280698
Original Clinical Research Quantitative
End-of-Life Care Among Patients With Kidney Failure on Maintenance Dialysis: A Retrospective Population-Based Study
https://orcid.org/0000-0002-2853-6881
Hafid Shuaib 1
Isenberg Sarina R. 23
Fernandes Aleisha 1
Gallagher Erin 1
Webber Colleen 34
Joseph Meera 5
Sood Manish M 6
Bruni Adrianna 7
Davis Janet L. 4
Warmels Grace 7
Downar James 7
Gayowsky Anastasia 8
Jones Aaron 89
Manuel Doug 341011
Tanuseputro Peter 3611
Howard Michelle 1
1 Department of Family Medicine, McMaster University, Hamilton, ON, Canada
2 Department of Medicine, University of Ottawa, ON, Canada
3 Bruyère Research Institute, University of Ottawa, ON, Canada
4 Ottawa Hospital Research Institute, The Ottawa Hospital, University of Ottawa, ON, Canada
5 Division of Nephrology, Department of Medicine, McMaster University, Hamilton, ON, Canada
6 Department of Medicine, Ottawa Hospital Research Institute, The Ottawa Hospital, University of Ottawa, ON, Canada
7 Division of Palliative Care, Department of Medicine, Ottawa Hospital Research Institute, The Ottawa Hospital, University of Ottawa, ON, Canada
8 ICES McMaster, McMaster University, Hamilton, ON, Canada
9 Department of Health Research Methods, Evidence, and Impact, McMaster University, Hamilton, ON, Canada
10 Department of Family Medicine, University of Ottawa, ON, Canada
11 ICES uOttawa, University of Ottawa, ON, Canada
Shuaib Hafid, Department of Family Medicine, McMaster University, 100 Main Street West, Hamilton, ON, L8P 1H6, Canada. Email: hafids3@mcmaster.ca
21 9 2024
2024
11 2054358124128069818 3 2024
23 7 2024
© The Author(s) 2024
2024
Canadian Society of Nephrology, unless otherwise noted. Manuscript content on this site is licensed under Creative Commons Licenses
https://creativecommons.org/licenses/by/4.0/ This article is distributed under the terms of the Creative Commons Attribution 4.0 License (https://creativecommons.org/licenses/by/4.0/) which permits any use, reproduction and distribution of the work without further permission provided the original work is attributed as specified on the SAGE and Open Access page (https://us.sagepub.com/en-us/nam/open-access-at-sage).
Background:

Nephrologists routinely provide end-of-life care for patients with kidney failure (KF) on maintenance dialysis. Involvement of primary care and palliative care physicians may enhance this experience.

Objective:

The objective was to describe outpatient care patterns in the last year of life and the end-of-life acute care utilization for patients with KF on maintenance dialysis.

Design:

Retrospective cohort study using population-level health administrative data.

Setting & Participants:

Outpatient and inpatient care during the last year of life among patients who died between 2017 and 2019, receiving maintenance dialysis in Ontario, Canada.

Measurements:

The primary exposure is patterns of physician specialties providing outpatient care in the last year of life. Outcomes include outpatient encounters in the last year of life, acute care visitation in the last month of life, and place of death.

Methods:

We reported the count and percentage of categorical outcomes and the median (interquartile range) for numeric outcomes. We produced time series plots of the mean monthly percentage of encounters to different specialties stratified by physician specialty patterns. We evaluated differences in outcomes by physician specialty patterns using analysis of variance (ANOVA) and Pearson’s chi-square tests (P < .05, two-tailed).

Results:

Among 6866 patients, the median age at death was 73, 36.1% were female, and 87.8% resided in urban regions. Three patterns emerged: a primary care, nephrology, and palliative care triad (25.5%); a primary care and nephrology dyad (59.3%); and a non-primary care pattern (15.2%). Palliative care involvement is concentrated near death. Of all, 81.4% spent at least 1 day in hospital or emergency department in the last month, but those with primary care, palliative care, and nephrology involvement had the fewest acute care deaths (65.8%).

Limitations:

Outpatient care patterns were defined using physician billing codes, potentially missing care from other providers.

Conclusions:

Nephrology and primary care predominantly manage outpatient care in the last year of life for patients with KF on maintenance dialysis, with consistent acute care use across care patterns except for the place of death. Future research should explore associations between patterns of care and end-of-life outcomes to identify the most optimal model of care for patients with KF on maintenance dialysis.

Abrégé

Contexte:

Il est courant pour les néphrologues de prodiguer des soins de fin de vie aux patients souffrant d’insuffisance rénale (IR) sous dialyse d’entretien. Cette expérience pourrait être enrichie par la participation des médecins des unités de soins primaires et de soins palliatifs.

Objectif:

Cette étude visait à décrire les modèles de soins ambulatoires prodigués au cours de la dernière année de vie et l’utilisation des soins aigus en fin de vie chez les patients atteints d’IR sous dialyse d’entretien.

Conception:

Étude de cohorte populationnelle rétrospective réalisée à partir des données administratives du système de santé.

Cadre et sujets de l’étude:

Les soins ambulatoires et hospitaliers au cours de la dernière année de vie chez les patients décédés sous dialyse d’entretien entre 2017 et 2019 en Ontario (Canada).

Mesures:

La principale mesure est le profil des spécialités médicales qui fournissent des soins ambulatoires dans la dernière année de vie. Les données recueillies comprennent les consultations externes au cours de la dernière année de vie, les visites en soins aigus au cours du dernier mois de vie et le lieu du décès.

Méthodologie:

Nous avons rapporté le nombre et le pourcentage de résultats catégoriels, ainsi que la médiane (écart interquartile) des résultats numériques. Nous avons produit des graphiques chronologiques du pourcentage mensuel moyen de consultations avec différentes spécialités, stratifiées selon les spécialités médicales. Nous avons évalué les différences dans les résultats selon les profils de spécialités médicales en utilisant les tests ANOVA et Chi-Square de Pearson (P <,05; bilatéral).

Résultats:

Des 6 866 patients inclus (âge médian au décès: 73 ans), 36,1% étaient des femmes et 87,8% vivaient en région urbaine. Trois modèles sont apparus: une triade soins primaires, néphrologie et soins palliatifs (25,5%); une dyade soins primaires et néphrologie (59,3%); et un modèle de soins non primaires (15,2%). La participation des soins palliatifs est concentrée autour du moment du décès. Une grande majorité des patients (81,4%) avait passé au moins une journée à l’hôpital ou aux urgences au cours du dernier mois, mais les personnes qui avaient bénéficié d’une triade de soins (primaires, néphrologie et soins palliatifs) présentaient une moins grande proportion de décès en soins aigus (65,8%).

Limites:

Les modèles de soins ambulatoires ont été définis à l’aide des codes de facturation des médecins, ce qui pourrait avoir exclu les soins dispensés par d’autres prestataires.

Conclusion:

Les soins ambulatoires au cours de la dernière année de vie des patients atteints d’IR sous dialyse d’entretien sont principalement prodigués par la néphrologie et les soins primaires, avec une utilisation constante des soins aigus dans tous les modèles de soins, sauf pour le lieu du décès. Les futures recherches devraient explorer les liens entre les modèles de soins et les résultats en fin de vie afin d’identifier le modèle de soins le plus optimal pour les patients atteints d’IR sous dialyse d’entretien.

cohort studies
dialysis
outpatient care
primary care
palliative care
Canadian Institutes of Health Research https://doi.org/10.13039/501100000024 180345 cover-dateJanuary-December 2024
typesetterts1
==== Body
pmcBackground

Patients with kidney failure experience high rates of mortality even if they receive dialysis. It is common for patients with kidney failure (KF) on maintenance dialysis to experience serious acute health events such as myocardial infarctions, and 3 out of 5 patients on dialysis experience at least one emergency department visit in the last month of life. 1 Patients with KF on maintenance dialysis will receive routine life-sustaining dialysis sessions, while also experiencing high symptom burden; this combination may be debilitating for patients approaching the end of life. Most patients with KF on maintenance dialysis do not receive palliative care from palliative care physicians despite experiencing a terminal illness, 1 and care from palliative care physicians is typically concentrated within the last 30 days of life for those who receive it. 2 This trend aligns with other studies of patients dying of organ failure, who receive fewer palliative care services compared to patients dying of cancer. 3

Patients with KF on maintenance dialysis receive regular care and high continuity with their nephrology team. 4 In addition, they experience a significant burden from co-morbid conditions,5,6 thus patients may have multiple physician specialties involved in their care and a regular primary care provider. Research from Ontario, Canada, found that high continuity with a primary care physician before and during dialysis initiation was not associated with decreased mortality risk or hospitalizations, 7 while a study from Alberta, Canada, identified that poor continuity of care with primary care was associated with increased acute care utilization patterns among patients with stage 3 or 4 chronic kidney disease. 8 Neither study focused on utilization patterns during the last year of life, and therefore, their findings may not reflect the end-of-life context regardless of the conflicting results.

Involving multiple specialties to care for different conditions may unintentionally undermine the known benefits of continuity of care. Primary care’s role among patients with KF on maintenance dialysis during the last year of life has yet to be extensively studied, despite, primary care being responsible for most of the outpatient care in the last year of life among the general population. 9 In addition, there is also evidence that an increase in continuity with primary care during end of life is associated with decreased acute care utilization.10-14 However, exactly which specialties of physicians they see, and how often they see them is unknown for patients dying with KF on maintenance dialysis. Further, the mix, timing, and intensity of physician care in the community and its relation to end-of-life healthcare utilization outcomes have not been previously described for this population. Hence, we examined the physician care team and the location of death among individuals with KF on maintenance dialysis during the last year of life.

Methods

Study Design and Data Sources

We conducted a retrospective cohort study using linked population-level health administrative data in Ontario, Canada (see Supplementary File 1 for a description of the datasets accessed). Administrative datasets are linked using unique encoded identifiers and analyzed at ICES (formerly known as the Institute for Clinical Evaluative Sciences).

Study Cohort

Our cohort included adult patients who died between January 1, 2017, and December 31, 2019, and who had any record of long-term maintenance dialysis initiation (including peritoneal dialysis) before their death date in the Canadian Organ Replacement Register (CORR) (see Supplementary File 1 for a description of data sources accessed). We excluded patients who were aged >105 years at death (in case of documentation errors in the patients’ birth or death dates) and those who were ineligible for the Ontario Health Insurance Plan (OHIP) at any point in the last year of life. Other exclusion criteria included patients who resided in a long-term care home in the last year of life, as these individuals receive different care patterns to patients residing in the community, and patients who did not have any weekly dialysis billing codes in the last year of life.

Physician Care Patterns

Our study exposure examines the mix of care provided by different physician specialties in the last year of life. Specifically, we focused on nephrology due to the nature of the population of interest, palliative care due to the end-of-life context, and primary care as primary care is responsible for most outpatient care provision in the last year of life for the general population. 9 We defined mutually exclusive categories of different outpatient physician patterns of care during the last year of life (Figure 2). First, patients who did not have an outpatient encounter with a primary care physician in the last year of life were categorized into the non-primary care pattern of care. Patients in this category may have any combination of encounters with all other specialties. Patients who did have an outpatient encounter with a primary care physician were categorized into the following categories: (1) primary care, nephrology, and palliative care triad and (2) primary care and nephrology care dyad. Patients were categorized into the care triad if they had at least one encounter with a primary care physician, nephrologist, and palliative care physician in the last year of life. Similarly, patients were categorized into the care dyad if they had at least one encounter with primary care and nephrology but did not have any palliative care encounters.

Physician specialty was obtained through OHIP billing records for all specialties except for palliative care as any physician can bill palliative care codes regardless of their specialty. We adapted a previously defined algorithm that is highly sensitive and has positive predictive value in identifying self-reported palliative care physicians. 15 The algorithm identified palliative care physicians as physicians with greater than 10% of their billings over 2 years as palliative care. We applied this definition to physicians who billed for palliative care services during the first calendar year of our study period (January 1, 2017, to December 31, 2017) and looked back 2-years before each of the unique palliative care billing codes to identify an average percentage of their billing patterns.

Outcomes

The outcomes of interest capture outpatient and inpatient care patterns in the last year of life. First, we measured the number of outpatient physician encounters by the outpatient physician care pattern categories in the last year of life, using individual outpatient billing codes in OHIP. We identified that this approach was inappropriate for nephrologist encounters who typically bill weekly maintenance dialysis billing codes and do not bill for each encounter. Hemodialysis patients will typically receive dialysis 3 times per week, wherein they interact with their dialysis care team, which can include the nephrologist and registered staff. As such, we supplemented the original definition by approximating outpatient encounters with nephrologists by the number of weekly maintenance dialysis OHIP billing codes (Supplementary File 2).

Next, we measured how end-of-life acute care utilization varies by patterns of physician care. Specifically, we examined the number of unique hospitalizations or emergency department visits in the last 30 days of life, the number of days spent in hospital or emergency department in the last 30 days of life, and the percentage of patients that died in acute care settings. We also reported patients’ acute care utilization in the last year of life, truncating the last 30 days of life. We used the Discharge Abstract Database and the National Ambulatory Care Reporting System to capture inpatient hospitalization and the Registered Persons Database to capture patient death locations: acute (i.e., emergency departments and inpatient hospital settings), sub-acute (i.e., complex continuing care and rehabilitation) or other (i.e., primarily in the community). 16

Patient Characteristics

Patient sociodemographic and clinical characteristics were reported. Neighborhood-level income and rurality status were identified using the patients’ postal code one year before their death date. Prevalence of chronic conditions was obtained using a 5-year look back period from the patients’ death date, using previously developed algorithms to assign the prevalence of 18 conditions.17-25 Comorbidity status at 1 year before death was also captured using the Johns Hopkins Adjusted Clinical Group (ACG)® System Aggregated Diagnosis Groups, Version 10. 26 All ICD-9/ICD-9-CM diagnosis codes assigned to patients are categorized into 32 functional groups, based on the diagnostic certainty, duration, severity, and etiology of the condition, and the likelihood of requiring specialist care. 27 Dialysis location, type, and duration of time between dialysis initiation and death were captured using the CORR. Patients’ primary care physician roster status was captured using the Client Agency Program Enrolment table for patients formally attached to a primary care physician and the rostered primary care physician was the primary care physician with the highest proportion of encounters in the last year of life, for patients not formally attached. 28

Analyses

Descriptive results were presented as proportions for categorical variables and as mean and standard deviation (SD), or medians and interquartile ranges (IQR) for skewed continuous variables. We reported study outcomes, stratified by the exposure variable, by producing means, SDs, medians, and IQRs. We also produced a time series plot reporting monthly percentages of outpatient encounters by physician specialties in the last year of life. We measured health utilization differences across patterns by conducting an Analysis of Variance (ANOVA) for numeric variables and Pearson’s chi-square tests for categorical variables. Statistical significance was set at P-value <.05 (two-tailed). All analyses were completed using SAS Enterprise Guide v. 7.15.

Ethics Approval

This study was approved by the Hamilton Integrated Research Ethics Board on March 28, 2022 (#14750-C). ICES is an independent, non-profit research institute whose legal status under Ontario’s health information privacy law (section 45 of Ontario’s Personal Health Information Protection Act) allows it to collect and analyze healthcare and demographic data, without consent, for health system evaluation and improvement.

Results

After applying the exclusion criteria, the cohort included 6,866 patients (Figure 1). Patients had a median age of 73 (64-81) at death, 36.1% were female, and 87.8% resided in urban regions (Table 1). 25.5% of patients experienced care from a primary care, nephrology, and palliative care triad; 59.3% experienced care from a primary care and nephrology dyad, while 15.2% experienced a non-primary care pattern in the last year of life (Figure 2). Dialysis modality and location did not differ across patterns of care, with approximately 70% of patients experiencing hemodialysis in-center. Patients who received care from a primary care, palliative care, and nephrology triad had a higher prevalence of cancer at 1-year before death (40.8%), compared to patients who received care from a primary care and nephrology dyad (24.1%), or patients who received the non-primary care pattern (19.2%). In addition, patients who received primary care, palliative care, and nephrology triad had the highest percentage of receiving provincial home care services in the last year of life (89.4%), compared to patients who experienced a primary care and nephrology dyad (77.7%) or a non-primary care pattern (74.7%). Patients who experienced a non-primary care pattern also had the greatest percentage of patients in the lowest neighborhood income quintile (34.1%) compared to the primary care and nephrology dyad (27.2%) and the primary care, palliative care, and nephrology care triad (26.2%).

Figure 1. Cohort creation flowchart.

Table 1. Profile of Patients Aged 19 or Older With Kidney F﻿ailure on Maintenance Dialysis Who Died Between January 1, 2017, and December 31, 2019, in Ontario, Canada, Reported by Physician Care Patterns.

Variable	Total Cohort	Pattern: Primary care, palliative care, nephrology triad	Pattern: Primary care, nephrology dyad	Pattern: Non-primary care	P-value	
N = 6866	N = 1752	N = 4069	N = 1045	
Age at death, N (%)						
 19 to 44	204 (3.0%)	32 (1.8%)	126 (3.1%)	46 (4.4%)	<.001	
 45 to 54	456 (6.6%)	80 (4.6%)	267 (6.6%)	109 (10.4%)	
 55 to 64	1097 (16.0%)	257 (14.7%)	636 (15.6%)	204 (19.5%)	
 65 to 74	1983 (28.9%)	522 (29.8%)	1207 (29.7%)	254 (24.3%)	
 75 to 84	2100 (30.6%)	553 (31.6%)	1259 (30.9%)	288 (27.6%)	
 85 to 94	987 (14.4%)	292 (16.7%)	561 (13.8%)	134 (12.8%)	
 95 to 104	39 (0.6%)	16 (0.9%)	13 (0.3%)	10 (1.0%)	
 Median (IQR)	73 (64-81)	74 (66-82)	73 (64-81)	71 (60-80)	<.001	
Female, N (%)	2480 (36.1%)	646 (36.9%)	1443 (35.5%)	391 (37.4%)	.337	
Urban residence, N (%)	6028 (87.8%)	1609 (91.8%)	3485 (85.7%)	934 (89.4%)	<.001	
Income quintile at 1-year before death, N (%)						
 1 (lowest)	1920 (28.0%)	459 (26.2%)	1105 (27.2%)	356 (34.1%)	<.001	
 2	1542 (22.5%)	391 (22.3%)	910 (22.4%)	241 (23.1%)	
 3	1304 (19.0%)	338 (19.3%)	797 (19.6%)	169 (16.2%)	
 4	1072 (15.6%)	267 (15.2%)	658 (16.2%)	147 (14.1%)	
 5 (highest)	995 (14.5%)	291 (16.6%)	586 (14.4%)	118 (11.3%)	
 Missing	33 (0.5%)	6 (0.3%)	13 (0.3%)	14 (1.3%)	
Adjusted diagnostic groups at 1-year before death,						
 Median (IQR)	12 (9-15)	12 (9-15)	12 (9-15)	12 (8-14)	<.001	
Prevalence of cancer at 1-year before death, N (%)	1895 (27.6%)	714 (40.8%)	980 (24.1%)	201 (19.2%)	<.001	
Prevalence of heart failure at 1-year before death, N (%)	2413 (35.1%)	582 (33.2%)	1502 (36.9%)	329 (31.5%)	<.001	
Duration of time (in years) between dialysis initiation and death date,						
 Mean (SD)	4.8 (5.9)	4.6 (5.9)	4.6 (5.8)	6.1 (6.5)	<.001	
Dialysis modality, N (%)						
 Only hemodialysis	4846 (70.6%)	1279 (73.0%)	2829 (69.5%)	738 (70.6%)	<.001	
 Only peritoneal dialysis	757 (11.0%)	194 (11.1%)	473 (11.6%)	90 (8.6%)	
 Both hemodialysis and peritoneal dialysis	1263 (18.4%)	279 (15.9%)	767 (18.9%)	217 (20.8%)	
Dialysis location, N (%)						
 In-center	4694 (68.4%)	1254 (71.6%)	2740 (67.3%)	700 (67.0%)	.002	
 Home	109 (1.6%)	23 (1.3%)	77 (1.9%)	9 (0.9%)	
 Both in-center and home	2063 (30.1%)	475 (27.1%)	1252 (30.8%)	336 (32.2%)	
Received provincial home care services in last year of life, N (%)	5508 (80.2%)	1566 (89.4%)	3161 (77.7%)	781 (74.7%)	<.001	
Rostered to a primary care physician, N (%)	6560 (95.5%)	1712 (97.7%)	3977 (97.7%)	871 (83.4%)	<.001	

Figure 2. Patterns of outpatient physician care experienced by patients dying with kidney failure on maintenance dialysis between January 1, 2017, and December 31, 2019, in Ontario, Canada.

Patients had a median of 51 (IQR: 16-53) outpatient encounters with nephrologists in the last year of life, regardless of which pattern they experienced (Table 2). Patients who experienced patterns of care that involved primary care experienced a median of 5 outpatient encounters with primary care physicians in the last year of life. Patients who experienced a primary care, palliative care, and nephrology triad experienced a median of 3 outpatient encounters with palliative care in the last year of life. Only 27.4% of patients who received a non-primary care pattern had an encounter with a palliative care physician in the last year of life.

Table 2. Outpatient Encounter Patterns During the Last Year of Life Among Patients Dying With Kidney F﻿ailure on Maintenance Dialysis Between January 1, 2017, and December 31, 2019, in Ontario, Canada.

Outcome	Total cohort	Pattern: Primary care, palliative care, nephrology triad	Pattern: Primary care, nephrology dyad	Pattern: Non-primary care	P-value	
N = 6866	N = 1752	N = 4069	N = 1045	
Primary Care:						
 Number of patients with ≥1 encounter, N (%)	5821 (84.8%)	1752 (100.0%)	4069 (100.0%)	0 (0.0%)	<.001	
 Number of outpatient encounters, Median (IQR)	4 (1-7)	5 (2-9)	5 (2-8)	0 (0-0)	<.001	
Palliative Care:						
 Number of patients with ≥1 encounter, N (%)	2038 (29.7%)	1752 (100.0%)	0 (0.0%)	286 (27.4%)	<.001	
 Number of outpatient encounters, Median (IQR)	0 (0-1)	3 (2-7)	0 (0-0)	0 (0-1)	<.001	
Nephrology:						
 Number of patients with ≥1 encounter, N (%)	6866 (100.0%)	1752 (100.0%)	4069 (100.0%)	1045 (100.0%)	-	
 Number of outpatient encounters, Median (IQR)	51 (16-53)	49 (9-52)	51 (16-53)	52 (44-53)	<.001	
Oncology:						
 Number of patients with ≥1 encounter, N (%)	661 (9.6%)	376 (21.5%)	212 (5.2%)	73 (7.0%)	<.001	
 Number of outpatient encounters, Median (IQR)	0 (0-0)	0 (0-0)	0 (0-0)	0 (0-0)	<.001	
Internal Medicine:						
 Number of patients with ≥1 encounter, N (%)	4759 (69.3%)	1217 (69.5%)	2839 (69.8%)	703 (67.3%)	.291	
 Number of outpatient encounters, Median (IQR)	2 (0-10)	2 (0-8)	2 (0-11)	2 (0-17)	.019	
Other Specialties:						
 Number of patients with ≥1 encounter, N (%)	6337 (92.3%)	1647 (94.0%)	3809 (93.6%)	881 (84.3%)	<.001	
 Number of outpatient encounters, Median (IQR)	6 (3-11)	7 (3-11)	6 (3-11)	4 (1-8)	<.001	

Almost all patients (92.3%) experienced an encounter with other specialties in the last year of life, with a median of 6 (IQR: 3-11) outpatient encounters.

Time series plots of the monthly percentages of outpatient encounters to different specialties identified primary care plays a consistent role in outpatient care provision in the last year of life, in addition to the consistent care provision by nephrologists (Figure 3, Supplementary File 3). Among patients who received care from the primary care, palliative care, and nephrology triad and those who received care from the non-primary care pattern, palliative care involvement progressively increased in the last 3 months of life. Regardless of pattern, almost 20% of monthly encounters were attributed to specialties other than primary care, palliative care, or nephrology.

Figure 3. Individual-level monthly percentages of outpatient encounters to different specialties during the last year of life among patients dying with kidney failure on maintenance dialysis between January 1, 2017, and December 31, 2019, in Ontario, Canada.

Patients who received care from the primary care, palliative care, and nephrology triad had the lowest percentages of any hospitalization (55.7%) or any emergency department visits (51.8%) in the last 30 days of life (Table 3), compared to patients who received other patterns of care. Regardless of these differences, all patterns experienced a median of 1 (IQR: 0-1) hospitalization and a median of 1 (IQR: 0-1) emergency department visit in the last 30 days of life. Despite experiencing the lowest percentages of hospitalizations or emergency department visits, patients who received care from the primary care, palliative care, and nephrology triad had the highest number of days spent in the hospital or emergency department (median 13 [IQR: 3-26]) compared to patients who experienced a primary care and nephrology dyad (median 8 [IQR: 1-22]) or a non-primary care pattern (median 8 [1-23]). Patients who received care from a primary care, palliative care, and nephrology pattern had the lowest percentage of deaths in acute care settings (65.8%), compared to patients who experienced a primary care and nephrology dyad (80.2%) or a non-primary care pattern (74.5%). Differences in acute care utilization during the last month of life across patterns were statistically significant (<.05).

Table 3. End-of-life healthcare outcomes among patients dying with kidney failure on maintenance dialysis between January 1, 2017, and December 31, 2019, in Ontario, Canada.

Variable	Value	Total cohort	Pattern: Primary care, palliative care, nephrology triad	Pattern: Primary care, nephrology dyad	Pattern: Non-primary care	P-value	
N = 6866	N = 1752	N = 4069	N = 1045	
Inpatient hospitalization in the last year of life (truncating the last 30 days of life)	Any, N (%)	5252 (76.5%)	1477 (84.3%)	3055 (75.1%)	720 (68.9%)	<.001	
Median (IQR)	2 (1-3)	2 (1-3)	2 (1-3)	1 (0-3)	<.001	
Inpatient hospitalization in the last 30 days of life	Any, N (%)	4092 (59.6%)	975 (55.7%)	2517 (61.9%)	600 (57.4%)	<.001	
Median (IQR)	1 (0-1)	1 (0-1)	1 (0-1)	1 (0-1)	<.001	
Emergency department visit in the last year of life (truncating the last 30 days of life)	Any, N (%)	5955 (86.7%)	1605 (91.6%)	3531 (86.8%)	819 (78.4%)	<.001	
Median (IQR)	3 (1-5)	3 (2-5)	3 (1-5)	2 (1-4)	<.001	
Emergency department visit in the last 30 days of life	Any, N (%)	3862 (56.3%)	908 (51.8%)	2383 (58.6%)	571 (54.6%)	<.001	
Median (IQR)	1 (0-1)	1 (0-1)	1 (0-1)	1 (0-1)	<.001	
Number of days spent in hospital or emergency department in last year of life (truncating the last 30 days of life)	Any, N (%)	6120 (89.1%)	1638 (93.5%)	3623 (89.0%)	859 (82.2%)	<.001	
Median (IQR)	18 (4-48)	26 (8-56)	16 (4-44)	15 (2-49)	<.001	
Number of days spent in hospital or emergency department in last 30 days of life	Any, N (%)	5591 (81.4%)	1458 (83.2%)	3304 (81.2%)	829 (79.3%)	.032	
Median (IQR)	9 (1-24)	13 (3-26)	8 (1-22)	8 (1-23)	<.001	
Place of death	Acute	5195 (75.7%)	1153 (65.8%)	3263 (80.2%)	779 (74.5%)	<.001	
Sub-acute	414 (6.0%)	198 (11.3%)	135 (3.3%)	81 (7.8%)	
Other	1257 (18.3%)	401 (22.9%)	671 (16.5%)	185 (17.7%)	

Acute care utilization in the last year of life (truncating the last 30 days of life) varied from the last month of life. Patients who received a Primary Care, Palliative Care, and Nephrology care triad had the highest percentages of any hospitalization (84.3%) or any emergency department visits (91.6%) in the 11 months before the last month of life (Table 3), compared to patients who received other patterns of care. However, patients who received the care triad also spent the most days in hospital or emergency department in the last year of life (median [IQR]: 26 [8-56]), compared to the primary care and nephrology dyad (median [IQR]: 16 [4-44]) or the non-primary care pattern (median [IQR]: 15 [2-49]). All differences across the last year of life were also statistically significant (<.05).

Discussion

Main Findings

In this study of patients who died with KF on maintenance dialysis in Ontario, Canada, patients received most of their outpatient care from nephrologists in the last year of life. This finding is unsurprising as patients with kidney failure who receive in-center hemodialysis receive dialysis 3 times a week and see a nephrologist once a week. 29 Therefore, their outpatient care is anchored to their dialysis treatment, and therefore, they have high continuity with their nephrology team. 4 Among patients who experienced primary care encounters, primary care physicians were also consistently involved in the last year of life, albeit on a smaller scale compared to nephrologists. This phenomenon may be specific to patients with KF on maintenance dialysis as primary care plays a much more prominent role compared to specialists in end-of-life outpatient care among the general decedent population. 9 Patients with KF on maintenance dialysis have traditionally looked to their nephrology team for advance care planning, symptom management, and social and emotional support. 30 This relationship is appropriate; however, primary care physicians may also be well-suited to provide these services due to their longitudinal relationships with their patients. 31 Given the potential benefits of primary care involvement in providing end-of-life care, it is notable that almost 15.5% of patients did not have any encounters with primary care during the last year of life; care was managed primarily by their dialysis team in these cases. Patients who wish to discontinue dialysis may experience a disconnect in care if primary care is not involved as they navigate non-dialysis-related healthcare services. Overall, most patients had continued contact with their primary care physician during their last year of life, in addition to other specialties.

Palliative care physician involvement increased as patients approached death, with large increases occurring in the last 3 months of life. Specialist palliative care is not being initiated early in the last year of life and is concentrated near death for patients dying with KF on maintenance dialysis. This finding aligns with other studies that identified that patients on dialysis access palliative care late in the end-of-life period.1,2,32,33 Patients with primary care involvement experienced greater palliative care involvement from both provincial home care services and palliative care physicians, suggesting that patients without primary care involvement may experience barriers to palliative care access or that nephrologists are providing end-of-life care for their patients. Furthermore, almost 1 in 5 monthly encounters were attributed to other specialists, beyond primary care, nephrology, or palliative care. Therefore, patients dying with KF on maintenance dialysis are experiencing multiple comorbidities that require medical attention from other specialists, resulting in more outpatient encounters during the last year of life.

End-of-life healthcare utilization did not vary extensively for patients dying with KF on maintenance dialysis across patterns of care. Almost 80% of patients in each pattern of care experienced at least one hospitalization or emergency department visit in the last 30 days of life—which is higher than previously reported research. 1 Patients dying with KF on maintenance dialysis may require more inpatient care provision at the end of life than other populations due to complications associated with reduced intervention effectiveness or other comorbidities (i.e., myocardial infarctions and infections). Patients who had palliative care physician involvement experienced lower rates of death in acute care settings, despite being more medically complex, as demonstrated by higher percentages of home care services enrolment, cancer prevalence, and spending more days in hospital or emergency departments in the last 30 days of life and having the greatest acute care utilization in the 11 months before the last 30 days of life. These findings complement existing literature that identified that patients on dialysis who received home palliative care services were less likely to die in acute or sub-acute care settings compared to those who did not receive home palliative care services. 1 Therefore, involving palliative care physicians may facilitate patients being discharged from inpatient settings to pass away in their homes, which is often the preferred location of death for patients.34-37 However, the timing and rationale of palliative care involvement remain uncertain, and future research ought to explore the characteristics associated with earlier palliative care referral.

Strengths and Limitations

This project has strengths, such as leveraging existing routinely collected health administrative data from Ontario, Canada. This data includes deidentified patient-level health information from approximately 13 million people with access to near-universal healthcare. As such, our findings may be relevant to other jurisdictions with public healthcare systems that also provide strong primary care services to their patients. This study also has limitations. First, we relied on different definitions for outpatient encounters with nephrology, as nephrologists do not typically bill for individual encounters for patients with KF on maintenance dialysis. As such, we may be overestimating the number of outpatient encounters that patients experience with their nephrologists by relying on individual encounter billing codes and weekly dialysis management billing codes. Some patients may receive primary care from nurse practitioners who do not have capturable OHIP physician billing codes. Therefore, patients in the non-primary care pattern may still benefit from having a primary care provider involved. Similarly, patients in the primary care and nephrology care dyad may have received palliative care services from their primary care physician or nephrologist. We also excluded patients who resided in long-term care homes during the last year of life as they may experience different outpatient care patterns compared to patients residing in the community. However, approximately 5% of chronic dialysis patients reside in long-term care homes in Ontario, and more long-term care homes have received funding to deliver peritoneal dialysis between 2012 and 2019, 38 and therefore, the number of patients receiving dialysis in long-term care is expected to increase. Next, our study only includes data before the COVID-19 pandemic, and therefore, may not reflect contemporary outpatient patterns. Lastly, this was a descriptive study, and therefore, cannot infer causal relationships between patterns of outpatient physician care and end-of-life healthcare utilization as there is temporal ambiguity. As a result, a separate study will be conducted to measure the associations between patterns of physician care and end-of-life healthcare outcomes, adjusting for patients’ demographic and clinical characteristics.

Conclusions

The end-of-life experience for patients dying with KF on maintenance dialysis in Ontario, Canada, heavily relies on nephrologists. Dialysis provision ensures patients experience high continuity of care with their nephrology care team, while primary care also remains a consistent, but lesser, presence. Acute care utilization during the last month of life is consistent across all patients with KF on maintenance dialysis regardless of the pattern of care, although patients are less likely to die in acute care settings when palliative care is involved. Future research ought to identify the associations between the identified patterns of outpatient care and end-of-life healthcare outcomes to identify the most optimal physician pattern of care for patients with KF on maintenance dialysis. As such, this endeavor may provide further justification for improved collaboration between nephrology, primary care and palliative care for patients at the end of life.

Supplemental Material

sj-docx-1-cjk-10.1177_20543581241280698 – Supplemental material for End-of-Life Care Among Patients With Kidney Failure on Maintenance Dialysis: A Retrospective Population-Based Study

Supplemental material, sj-docx-1-cjk-10.1177_20543581241280698 for End-of-Life Care Among Patients With Kidney Failure on Maintenance Dialysis: A Retrospective Population-Based Study by Shuaib Hafid, Sarina R. Isenberg, Aleisha Fernandes, Erin Gallagher, Colleen Webber, Meera Joseph, Manish M Sood, Adrianna Bruni, Janet L. Davis, Grace Warmels, James Downar, Anastasia Gayowsky, Aaron Jones, Doug Manuel, Peter Tanuseputro and Michelle Howard in Canadian Journal of Kidney Health and Disease

sj-docx-2-cjk-10.1177_20543581241280698 – Supplemental material for End-of-Life Care Among Patients With Kidney Failure on Maintenance Dialysis: A Retrospective Population-Based Study

Supplemental material, sj-docx-2-cjk-10.1177_20543581241280698 for End-of-Life Care Among Patients With Kidney Failure on Maintenance Dialysis: A Retrospective Population-Based Study by Shuaib Hafid, Sarina R. Isenberg, Aleisha Fernandes, Erin Gallagher, Colleen Webber, Meera Joseph, Manish M Sood, Adrianna Bruni, Janet L. Davis, Grace Warmels, James Downar, Anastasia Gayowsky, Aaron Jones, Doug Manuel, Peter Tanuseputro and Michelle Howard in Canadian Journal of Kidney Health and Disease

sj-docx-3-cjk-10.1177_20543581241280698 – Supplemental material for End-of-Life Care Among Patients With Kidney Failure on Maintenance Dialysis: A Retrospective Population-Based Study

Supplemental material, sj-docx-3-cjk-10.1177_20543581241280698 for End-of-Life Care Among Patients With Kidney Failure on Maintenance Dialysis: A Retrospective Population-Based Study by Shuaib Hafid, Sarina R. Isenberg, Aleisha Fernandes, Erin Gallagher, Colleen Webber, Meera Joseph, Manish M Sood, Adrianna Bruni, Janet L. Davis, Grace Warmels, James Downar, Anastasia Gayowsky, Aaron Jones, Doug Manuel, Peter Tanuseputro and Michelle Howard in Canadian Journal of Kidney Health and Disease

sj-docx-4-cjk-10.1177_20543581241280698 – Supplemental material for End-of-Life Care Among Patients With Kidney Failure on Maintenance Dialysis: A Retrospective Population-Based Study

Supplemental material, sj-docx-4-cjk-10.1177_20543581241280698 for End-of-Life Care Among Patients With Kidney Failure on Maintenance Dialysis: A Retrospective Population-Based Study by Shuaib Hafid, Sarina R. Isenberg, Aleisha Fernandes, Erin Gallagher, Colleen Webber, Meera Joseph, Manish M Sood, Adrianna Bruni, Janet L. Davis, Grace Warmels, James Downar, Anastasia Gayowsky, Aaron Jones, Doug Manuel, Peter Tanuseputro and Michelle Howard in Canadian Journal of Kidney Health and Disease

List of Abbreviations: CORR, Canadian Organ Replacement Registry; IQR, interquartile range; KF, kidney failure; OHIP, Ontario Health Insurance Plan; SD, standard deviation.

Ethics Approval: This study was approved by the Hamilton Integrated Research Ethics Board on March 28, 2022 (#14750-C).

Consent to Participate: Informed consent was not required as we used routinely collected deidentified health administrative data.

Consent for Publication: Not applicable.

Availability of Data and Materials: The data set from this study is held securely in coded form at ICES. While data-sharing agreements prohibit ICES from making the data set publicly available, access may be granted to those who meet pre-specified criteria for confidential access, available at www.ices.on.ca/DAS. The full data set creation plan and underlying analytic code are available from the authors upon request, understanding that the computer programs may rely upon coding templates or macros that are unique to ICES and are therefore either inaccessible or may require modification.

Author Contributions: SRI, MH, SH, and AF conceived the study. All authors designed the study and interpreted the results. AG analyzed the data. SH and SRI wrote the manuscript. All authors revised the manuscript critically for important intellectual content, gave final approval of the version to be published, and agreed to be accountable for all aspects of the work.

Authors’ Note: The findings of this study have been previously presented at the Canadian Association for Health Services & Policy Research Annual Conference, May 29-31, 2023, Montreal, Canada, and the Canadian Hospice & Palliative Care Association Conference, October 12-14, 2023, Ottawa, Canada.

The author(s) declared the following potential conflicts of interest with respect to the research, authorship, and/or publication of this article: MMS has received consultancy fees from AstraZeneca, Bayer, Otsuka, and GlaxoSmithKline. The remaining authors declared no potential conflicts of interest concerning the research, authorship, and/or publication of this article.

Funding: The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This study was supported by ICES, which is funded by an annual grant from the Ontario Ministry of Health (MOH) and the Ministry of Long-Term Care (MLTC). This study also received funding from the Canadian Institutes of Health Research grant #180345. This document used data adapted from the Statistics Canada Postal CodeOM Conversion File, which is based on data licensed from Canada Post Corporation and/or data adapted from the Ontario Ministry of Health Postal Code Conversion File, which contain data copied under license from ©Canada Post Corporation and Statistics Canada. Parts of this material are based on data and/or information compiled and provided by the Ontario MOH, Canadian Institute for Health Information, and MLTC. The analyses, conclusions, opinions, and statements expressed herein are solely those of the authors and do not reflect those of the funding or data sources. No endorsement by ICES, the MOH, or the MLTC is intended or should be inferred.

ORCID iD: Shuaib Hafid https://orcid.org/0000-0002-2853-6881

Supplemental Material: Supplemental material for this article is available online.
==== Refs
References

1 Nesrallah GE Dixon SN MacKinnon M , et al . Home palliative service utilization and care trajectory among Ontario residents dying on chronic dialysis. Can J Kidney Health Dis. 2018;5 :1-11. doi:10.1177/2054358118783761.
2 Redahan L Brady B Smyth A Higgins S Wall C . The use of palliative care services amongst end-stage kidney disease patients in an Irish tertiary referral centre. Clin Kidney J. 2013;6 (6 ):604-608. doi:10.1093/ckj/sft117.26069830
3 Saini T Murtagh FE Dupont PJ McKinnon PM Hatfield P Saunders Y . Comparative pilot study of symptoms and quality of life in cancer patients and patients with end stage renal disease. Palliat Med. 2006;20 (6 ):631-636. doi:10.1177/0269216306070236.17060257
4 Chen AY Chen B Kuo C-C . Better continuity of care improves the quality of end-of-life care among elderly patients with end-stage renal disease. Sci Rep. 2020;10 (1 ):19716. doi:10.1038/s41598-020-76707-w.33184374
5 Fraser SD Roderick PJ May CR , et al . The burden of comorbidity in people with chronic kidney disease stage 3: a cohort study [published correction appears in BMC Nephrol. 2020;21(1):543]. BMC Nephrol. 2015;16 :193. doi:10.1186/s12882-015-0189-z.26620131
6 Gomez AT Kiberd BA Royston JP , et al . Comorbidity burden at dialysis initiation and mortality: a cohort study. Can J Kidney Health Dis. 2015;2 :34. doi:10.1186/s40697-015-0068-3.26351568
7 Silver SA Bota SE McArthur E , et al . Association of primary care involvement with death or hospitalizations for patients starting dialysis. Clin J Am Soc Nephrol. 2020;15 (4 ):521-529. doi:10.2215/CJN.10890919.32139363
8 Chong C Campbell D Elliott M Aghajafari F Ronksley P . Determining the association between continuity of primary care and acute care use in chronic kidney disease: a retrospective cohort study. Ann Fam Med. 2022;20 (3 ):237-245. doi:10.1370/afm.2813.35606125
9 Howard M Hafid A Isenberg SR , et al . Intensity of outpatient physician care in the last year of life: a population-based retrospective descriptive study. CMAJ Open. 2021;9 (2 ):E613-E622. doi:10.9778/cmajo.20210039.
10 Burge F Lawson B Johnston G Cummings I . Primary care continuity and location of death for those with cancer. J Palliat Med. 2003;6 (6 ):911-918. doi:10.1089/109662103322654794.14733683
11 Burge F Lawson B Johnston G . Family physician continuity of care and emergency department use in end-of-life cancer care. Med Care. 2003;41 (8 ):992-1001. doi:10.1097/00005650-200308000-00012.12886178
12 Jones A Bronskill SE Seow H Junek M Feeny D Costa AP . Associations between continuity of primary and specialty physician care and use of hospital-based care among community-dwelling older adults with complex care needs [published correction appears in PLoS ONE. 2021;16(10):e0258708]. PLoS ONE. 2020;15 (6 ):e0234205. doi:10.1371/journal.pone.0234205.
13 Almaawiy U Pond GR Sussman J Brazil K Seow H . Are family physician visits and continuity of care associated with acute care use at end-of-life? A population-based cohort study of homecare cancer patients. Palliat Med. 2014;28 (2 ):176-183. doi:10.1177/0269216313493125.23779252
14 Isaac B Adam S Sarah RD . Association between continuity of care in general practice and hospital admissions for ambulatory care sensitive conditions: cross sectional study of routinely collected, person level data. BMJ. 2017;356 :j84. doi:10.1136/bmj.j84.
15 Barbera L Hwee J Klinger C Jembere N Seow H Pereira J . Identification of the physician workforce providing palliative care in Ontario using administrative claims data. CMAJ Open. 2015;3 (3 ):E292-E298. doi:10.9778/cmajo.20150005.
16 Tanuseputro P Beach S Chalifoux M , et al . Associations between physician home visits for the dying and place of death: a population-based retrospective cohort study. PLoS ONE. 2018;13 (2 ):e0191322. doi:10.1371/journal.pone.019132217.
17 Muggah E Graves E Bennett C Manuel DG . The impact of multiple chronic diseases on ambulatory care use; a population based study in Ontario, Canada. BMC Health Serv Res. 2012;12 :452. doi:10.1186/1472-6963-12-452.23217125
18 Lane NE Maxwell CJ Gruneir A Bronskill SE Wodchis WP . Absence of a socioeconomic gradient in older adults’ survival with multiple chronic conditions. EBioMedicine. 2015;2 (12 ):2094-2100.26844290
19 Jaakkimainen RL Bronskill SE Tierney MC , et al . Identification of physician-diagnosed Alzheimer’s disease and related dementias in population-based administrative data: a validation study using family physicians’ electronic medical records. J Alzheimers Dis. 2016;54 (1 ):337-349. doi:10.3233/JAD-160105.27567819
20 Mondor L Cohen D Khan AI Wodchis WP . Income inequalities in multimorbidity prevalence in Ontario, Canada: a decomposition analysis of linked survey and health administrative data. Int J Equity Health. 2018;17 (1 ):90. doi:10.1186/s12939-018-0800-6.29941034
21 Mondor L Maxwell CJ Hogan DB , et al . Multimorbidity and healthcare utilization among home care clients with dementia in Ontario, Canada: a retrospective analysis of a population-based cohort. PLoS Med. 2017;14 :e1002249.
22 Gruneir A Bronskill SE Maxwell CJ , et al . The association between multimorbidity and hospitalization is modified by individual demographics and physician continuity of care: a retrospective cohort study. BMC Health Serv. 2016;16 :154.
23 Pefoyo AJ Bronskill SE Gruneir A , et al . The increasing burden and complexity of multimorbidity. BMC Public Health. 2015;15 :415. doi:10.1186/s12889-015-1733-2.25903064
24 Petrosyan Y Bai YQ Koné Pefoyo AJ , et al . The relationship between diabetes care quality and diabetes-related hospitalizations and the modifying role of comorbidity. Can J Diabetes. 2017;41 (1 ):17-25.27789111
25 Thavorn K Maxwell CJ Gruneir A , et al . Effect of socio-demographic factors on the association between multimorbidity and healthcare costs: a population-based, retrospective cohort study. BMJ Open. 2017;7 (10 ):e017264. doi:10.1136/bmjopen-2017-017264.
26 Weiner JP . The John Hopkins ACG® Case-Mix System Version 10.0 Release Notes. Baltimore, MD, Health Services Research & Development Center, The Johns Hopkins University Bloomberg School of Public Health; 2013.
27 Weiner JP . The John Hopkins ACG® Case-Mix System Documentation & Application Manual Version 10.0. Baltimore, MD, The John Hopkins University Bloomberg School of Public Health, Health Services Research & Development Center; 2013.
28 Howard M Chalifoux M Tanuseputro P . Does primary care model effect healthcare at the end of life? A population-based retrospective cohort study. J Palliat Med. 2017;20 :344-351.27893954
29 Blake PG . Global dialysis perspective: Canada. Kidney360. 2020;1 (2 ):115-118. doi:10.34067/KID.0000462019.35372909
30 Davison SN . End-of-life care preferences and needs: perceptions of patients with chronic kidney disease. Clin J Am Soc Nephrol. 2010;5 (2 ):195-204. doi:10.2215/CJN.05960809.20089488
31 Olaisen RH Schluchter MD Flocke SA Smyth KA Koroukian SM Stange KC . Assessing the longitudinal impact of physician-patient relationship on functional health. Ann Fam Med. 2020;18 (5 ):422-429. doi:10.1370/afm.2554.32928758
32 Teruel JL Rexach L Burguera V , et al . Home palliative care for patients with advanced chronic kidney disease: preliminary results. Healthcare (Basel). 2015;3 (4 ):1064-1074. doi:10.3390/healthcare3041064.27417813
33 van Biesen W van de Luijtgaarden MWM Brown EA , et al . Nephrologists’ perceptions regarding dialysis withdrawal and palliative care in Europe: lessons from a European Renal Best Practice survey. Nephrol Dial Transplant. 2015;30 (12 ):1951-1958. doi:10.1093/ndt/gfv284.26268713
34 Kastbom L Milberg A Karlsson M . A good death from the perspective of palliative cancer patients. Support Care Cancer. 2017;25 (3 ):933-939. doi:10.1007/s00520-016-3483-9.27837324
35 Brazil K Howell D Bedard M Krueger P Heidebrecht C . Preferences for place of care and place of death among informal caregivers of the terminally ill. Palliat Med. 2005;19 (6 ):492-499. doi:10.1191/0269216305pm1050oa.16218162
36 Gomes B Higginson IJ Calanzani N , et al . Preferences for place of death if faced with advanced cancer: a population survey in England, Flanders, Germany, Italy, the Netherlands, Portugal and Spain. Ann Oncol. 2012;23 (8 ):2006-2015. doi:10.1093/annonc/mdr602.22345118
37 Brumley R Enguidanos S Jamison P , et al . Increased satisfaction with care and lower costs: results of a randomized trial of in-home palliative care. J Am Geriatr Soc. 2007;55 (7 ):993-1000. doi:10.1111/j.1532-5415.2007.01234.x.17608870
38 Blake PG McCormick BB Taji L , et al . Growing home dialysis: the Ontario renal network home dialysis initiative 2012–2019. Perit Dial Int. 2021;41 (5 ):441-452. doi:10.1177/08968608211012805.33969759
