Introduction

Pyoderma gangrenosum (PG) is a neutrophilic inflammatory skin disease that typically progresses over days from a painful nodule, plaque, or pustule to ulceration with undermined, violaceous borders and a surrounding zone of erythema [1]. Different variants of PG exist, such as ulcerative, bullous, pustular, vegetative, peristomal, and post-surgical [2]. Based on numerous retrospective studies, the annual incidence of PG was estimated at 3 to 10 cases per million [2, 3].

PG is associated with underlying systemic conditions in 56.8% of patients, with inflammatory bowel disease, inflammatory arthritis, and hematological malignancies being the most frequently encountered underlying comorbidities [4]. A better characterization of underlying diseases is of substantial importance as the type and severity of these comorbidities bear a prognostic value for PG [5, 6]. PG is a prototype of neutrophilic dermatoses, a polymorphous group of noncontagious dermatological disorders characterized by non-cutaneous findings including systemic inflammation, non-specific inflammatory findings (e.g., myalgia, fever, and joint pain), and sterile neutrophilic infiltrates found in organs other than skin [7, 8]. The visceral involvement may support the systemic nature of PG among other neutrophilic dermatoses [8].

While the most common extracutaneous manifestation of PG is pulmonary conditions, some case reports described its association with chronic renal comorbidities [7]. However, the association of PG with renal comorbidities is yet to be systematically investigated. The aim of this real-world data analysis was to evaluate the bidirectional association between PG and renal comorbidities.

Methods

Study design and database

The current study was performed to investigate the bidirectional association between PG and renal comorbidities. To outline the risk of developing renal diseases during the course of PG, a retrospective cohort study design was adopted, in which patients with PG were followed to estimate the incidence of renal diseases. To evaluate the risk of having PG in individuals with a preceding history of renal diseases, a case-control study design was implemented to disclose the prevalence of preexisting renal diseases (exposure) in patients with subsequent PG (outcome).

The current study was grounded on the computerized database of Clalit Health Services (CHS). Ensuring 4.5 million enrollees as of October 2018, CHS is the largest health care organization in Israel, providing healthcare services for 57% of the general Israeli population. The characteristics of the utilized dataset are further detailed in our previous publications [9, 10]. The current study was approved by the institutional review board (IRB) of Ben-Gurion University in accordance with the Declaration of Helsinki.

Study population

By a systematic check of the dataset of CHS, all individuals with a diagnosis of PG between the years 2000 and 2018 were identified. Subsequently, cases were checked, and only those fulfilling one of the following eligibility criteria were finally eligible: (i) a documented diagnosis of PG registered at least twice by a community board-certified dermatologist, and/or (ii) documentation of the diagnosis of PG in discharge letters of patients admitted to dermatological wards.

All renal comorbidities in the chronic registry of CHS were evaluated in the current study. These are chronic renal failure (CRF), dialysis, kidney transplantation (KT), and other kidney diseases (OKD). OKD was defined based on the diagnostic codes detailed in Supplementary Table 1. Given that these diagnoses were registered in the chronic registry of CHS, they were based on registration by a board-certified nephrologist, suggestive laboratory and imaging data, and were eventually authenticated by the managing general healthcare provider.

Table 1 Descriptive characteristics of the study population

The control group included up to 5 control individuals per patient, matched randomly by age, sex, and ethnicity. Prior to their recruitment, control subjects were confirmed to be alive and to contribute longitudinal data for the CHS dataset.

Main covariates of the study

The study outcome measures were controlled for underlying comorbidities as evaluated by the Charlson comorbidity index (CCI), a validated epidemiological method of quantifying comorbidities. This index has been found to be reliable in predicting mortality and is widely utilized in epidemiological studies [2]. Smoking status was classified as a current smoker or never/past smoker.

Statistical analysis

Baseline characteristics were described by means and standard deviations (SDs) for continuous variables, while categorical values were signified by percentages. A comparison of sociodemographic and clinical factors between cases and controls was performed using the chi-square test and t-test for categorical and continuous variables, respectively.

In the cohort study design, incidence rates of renal comorbidities were calculated for both PG patients and controls and expressed as the number of events per 1000 person-years. Hazard ratios (HRs) for the risk of incident renal comorbidities were obtained by the use of the Cox regression model. In the case-control study design, logistic regression was used to calculate odds ratios (ORs) and 95% confidence intervals (CIs) to compare cases and controls regarding the presence of preceding renal comorbidities. The association was calculated based on individuals who developed PG after the diagnosis of each renal disease, given that a temporal relationship exists between exposure and outcome in case-control studies. Two-tailed P values less than 0.05 were considered statistically significant. All statistical analyses were performed using the SPSS software, version 25 (SPSS, Armonk, NY: IBM Corp).

Results

Characteristics of the study population

The study population included 302 patients with PG and 1497 age-, sex-, and ethnicity-matched control individuals. The mean (SD) age of the study participants was 54.0 (20.8) years, 57.9% of them were females, and 84.4% were of Jewish ethnicity. Case and control groups were comparable with regard to mean body mass index (BMI) and the prevalence of smoking (Table 1). The mean (SD) CCI score was greater in cases than in controls (2.3 [2.7] vs. 1.3 [1.8], respectively; P<0.001). The prevalence of severe comorbid conditions was greater among cases relative to controls (37.4% vs. 19.2%, respectively; P<0.001). The characteristics of the study population are outlined in Table 1.

The risk of PG among individuals with a history of renal disease

A case-control study design was followed to estimate the prevalence of preexisting renal diseases among patients with PG and control individuals. The latter reflects the risk of developing subsequent PG among those with a history of renal diseases. The subsequent development of PG was significantly associated with a history of CRF (OR, 2.04; 95% CI, 1.20–3.47), KT (OR, 5.01; 95% CI, 1.01–24.94), and OKD (OR, 2.04; 95% CI, 1.20–3.47). In a sex-stratified analysis, PG was significantly associated with CRF in both sexes, whereas the association of PG with KT and OKD retained its statistical significance only among females (Table 2).

Table 2 The risk of pyoderma gangrenosum among patients with a preexisting history of renal diseases (case-control study design)

We then performed a multivariable logistic regression analysis adjusting for putative confounding variables. The latter demonstrated an independently increased risk of PG in individuals with a history of CRF (adjusted OR, 2.34; 95% CI, 1.33–4.11), KT (adjusted OR, 5.03; 95% CI, 1.01–25.12), and OKD (adjusted OR, 1.69; 95% CI, 1.04–2.74), but not in those with a history of dialysis (adjusted OR, 3.79; 95% CI, 0.84–17.15; Table 2).

The risk of renal diseases among individuals with PG

A retrospective cohort study was performed to disclose the risk of developing renal diseases among patients with PG. The incidence rate of CRF, dialysis, KT, and OKD among patient with PG was 20.0 (95% CI, 13.7–28.5), 3.19 (95% CI, 1.17–7.08), 0.63 (95% CI, 0.03–3.14), and 28.1 (95% CI, 20.4–37.9)/1000 person-year, respectively (Table 3).

Table 3 The risk of renal diseases among patients with pyoderma gangrenosum (retrospective cohort study design)

A significantly increased risk of CRF (HR, 3.19; 95% CI, 1.97–5.16), dialysis (HR, 26.90; 95% CI, 3.14–230.32), and OKD (HR, 2.58; 95% CI, 1.48–4.52) was observed among patients with PG. While the risk of CRF and dialysis was more prominent among females, males were more predisposed to OKD (Table 3). When adjustment for putative confounding factors was performed, patients with PG were found to have an independently elevated risk of CRF (adjusted HR, 3.68; 95% CI, 2.72–5.97), dialysis (adjusted HR, 27.79; 95% CI, 3.24–238.14), and OKD (adjusted HR, 2.71; 95% CI, 1.55–4.74; Table 3).

Features of PG patients with renal diseases relative to the remaining patients with PG

Table 4 delineates the differential demographic and epidemiological features of patients with PG and coexistent renal comorbidities as compared to the remaining patients with PG. It was invariably shown that patients with a comorbid renal disease presented with PG at a significantly higher age and had higher CCI (Table 4). Patients with PG and coexistent CRF had a more prominent Jewish ethnicity, whereas patients with PG and dialysis had an increased BMI (Table 4).

Table 4 Comparison between patients with pyoderma gangrenosum and renal comorbidities relative to the remaining patients with pyoderma gangrenosum

Discussion

The current study indicates an increased risk of CRF, dialysis, and OKD among patients with PG relative to controls. Moreover, the presence of a medical history of CRF, KT, and OKD was found to predispose individuals to develop subsequent PG. Relative to the remaining patients with PG, those with PG and chronic renal comorbidities presented with PG at an older age and had an increased burden of comorbid conditions.

The risk of renal comorbidities among patients with PG was not previously explored by controlled observational studies. However, our results, signifying an increased risk of the former, are substantiated by previous case reports of renal involvement in PG manifesting as kidney lesion [11], increase in serum creatinine [12], microhematuria [13], glomerulonephritis [14], aseptic leukocyturia, oliguria, and proteinuria [15]. In a German retrospective, multi-center study encompassing 259 patients with PG, 6.1% of the study participants had renal failure [12].

The increased risk of PG in patients with preexisting renal disease observed in our analysis has also been reported in several case reports. The most frequent renal conditions that were implicated in eliciting PG were KT [16,17,18], end-stage renal disease (ESRD) [19,20,21] CRF [22], and renal carcinoma [23].

The pathomechanism underlying the association between PG and renal comorbidities is still unclear. However, several hypotheses have been postulated to account for it. As ulcerative PG oftentimes occurs following minor injuries or surgical operations, a phenomenon coined as pathergy [4], some authors suggested that vascular insufficiency might be involved in the pathogenesis of PG. The latter may be of relevance in the context of PG following radio-cephalic arteriovenous fistula for dialytic treatment [19]. Other authors proposed that the presence of monoclonal gammopathy, which is frequently implicated in triggering PG [3, 4], might lead to a rapid progression of renal failure [14], thus accounting, at least in part, for the observed association. Our outcome measures persisted following the adjustment for CCI which encompasses, among others, a diagnosis of hematological malignancies (leukemia and lymphoma). This multivariate analysis denotes an independent association between PG and renal comorbidities.

This study is epitomized by several strengths such as the population-based setting, the long follow-up period, the relatively large sample of PG patients, the access to all levels of healthcare facilities, and the systematic collection of data. However, a number of limitations in the current study should be acknowledged. The retrospective nature of this study and the use of diagnostic codes limited the ability to collect standardized information on the cutaneous manifestation and severity of the cutaneous disease. To elaborate, differentiation between the ulcerative, vegetative, and atypical subtypes of PG was not feasible using the current dataset. In addition, data regarding exposure to different therapeutic agents could not be retrieved for the current study. The latter is of significance given that cyclosporine, a frequently utilized drug in PG, is implicated with an increased risk of nephrotoxicity [24]. While the likelihood of residual confounding could not be thoroughly refuted, we do not think it had conferred a meaningful impact on the study’s main findings.

To conclude, the current retrospective cohort study disclosed that patients with PG experience an increased risk of developing CRF, dialysis, and OKD. The presence of preexisting CRF, KT, and OKD was associated with the subsequent emergence of PG. Patients with comorbid PG and renal comorbidities present with PG at an older age and have a higher comorbidity score. Further prospective studies are necessary to better characterize the influence of renal impairment on PG and vice versa. A more in-depth investigation may shed light on the underlying pathophysiology of the association between this debilitating dermatologic dermatosis and renal comorbidities.