New-onset diabetes after transplantation (NODAT) is a major complication in renal transplant recipients (RTRs). Cholesterol metabolism has been linked to diabetes development. Proprotein convertase subtilisin/kexin type 9 (PCSK9) is crucial in LDL receptor regulation. Its association with NODAT is unknown. We prospectively determined the association between serum PCSK9 levels and NODAT development and then with all-cause mortality, cardiovascular mortality, and renal graft failure.
In a university setting, nondiabetic RTRs recruited between 2001 and 2003 with a functional graft for ≥1 year were eligible. Serum PCSK9 was measured by ELISA. Cox proportional hazards analysis was used to assess the association of PCSK9 with the development of NODAT, all-cause mortality, cardiovascular mortality, and graft failure.
In 453 RTRs (age 51 ± 12 years, 56% male; 6.1 [2.7–11.7] years after transplantation), serum PCSK9 was 107.1 ± 43.4 μg/L. During a median follow-up of 10 years, 70 RTRs developed NODAT, 123 died, and 59 developed graft failure. NODAT occurred more frequently in the upper PCSK9 tertile (23%) versus the lowest two PCSK9 tertiles (12%; P < 0.001). In crude Cox regression analyses, PCSK9 was significantly associated with development of NODAT (hazard ratio 1.34 [95% CI 1.10–1.63]) per SD change (P = 0.004). This association remained independent of adjustment for potential confounders, including statin use. PCSK9 was not associated with all-cause mortality, cardiovascular mortality, or graft failure.
Circulating PCSK9 is associated with NODAT in RTRs. The PCSK9 pathway may contribute to the pathogenesis of NODAT.
Introduction
In renal transplant recipients (RTRs), new-onset diabetes after transplantation (NODAT) is a major metabolic complication (1). Its incidence has increased during the past decades (1,2). Estimates of the frequency of NODAT vary widely and have been reported to range from 2% to 50% (3). Twelve months cumulative incidence of NODAT was highest in the era of high-dose steroid regimens in the 1970s, lowest with cyclosporine-based regimens, and again higher with tacrolimus-based regimens (3,4). Although the observation period of many studies does not go beyond an initial surge in the incidence of NODAT during the first year after transplantation, it is well documented that during the era of cyclosporine-based regimens, the largest number of incident cases of NODAT occurred beyond the first year after transplantation (1). With current tacrolimus-based regimens, the number of the incident cases that occur beyond the first year after transplantation is even higher (4).
That NODAT adversely affects patient and graft survival in RTRs is well appreciated (5,6). Therefore, better identifying which RTRs are at increased risk of developing NODAT is clinically relevant. RTRs also are prone to the development of atherosclerotic cardiovascular disease (CVD), which represents a major cause of graft loss and a leading cause of death (7). Apart from affecting CVD risk, dyslipidemia has been found to be independently associated with reduced graft survival (8). Abnormalities in cholesterol metabolism thus can be considered to play a key role in the development of CVD and graft failure in RTRs.
A meta-analysis revealed that statin therapy increases the risk of new-onset diabetes in nontransplant populations (9). Statins lower LDL cholesterol by inhibiting hydroxymethylglutaryl-CoA reductase (10), thereby increasing LDL receptor (LDLR) expression and promoting cellular cholesterol delivery (11). Conversely, the risk of type 2 diabetes is markedly reduced in patients with familial hypercholesterolemia as a consequence of LDLR mutations (12). Moreover, diabetes incidence may be linked to genetic variation in hydroxymethylglutaryl-CoA reductase, indicating a role of cholesterol metabolism in diabetes development (9).
The proprotein convertase subtilisin/kexin type 9 (PCSK9) pathway has been identified to be crucially involved in regulating hepatic LDLR expression in humans (13). PCSK9 is a secreted protease that is predominately expressed by the liver and intestine (14). PCSK9 is able to target the LDLR toward intracellular degradation, thereby preventing LDLR from recycling to the cell surface (13). Higher circulating PCSK9 levels confer a lower rate of LDL apoB catabolism (15), supporting the notion that variations in circulating PCSK9 levels are pathophysiologically relevant.
To date, the association of serum levels of PCSK9 with NODAT development has not been evaluated. Because the PCSK9 pathway plays a central role in cholesterol metabolism, this study was primarily initiated to prospectively investigate the association between serum PCSK9 levels and NODAT development in RTRs. Secondary end points assessed were the association of serum PCSK9 with all-cause mortality, cardiovascular mortality, and renal graft failure in RTRs.
Research Design and Methods
Study Design
In this prospective cohort study, all adult RTRs who survived with a functioning allograft beyond the first year after transplantation were considered eligible to participate. Baseline data, obtained at least 1 year after transplantation, were collected between August 2001 and July 2003 at a median of 6.0 (interquartile range [IQR] 2–11) years after transplantation. Patients with known or apparent systemic illnesses (i.e., malignancies, opportunistic infections) were excluded from participation. Six hundred six of 847 (72%) eligible RTRs gave informed consent. We excluded 105 RTRs with existing diabetes (defined as fasting plasma glucose ≥126 mg/dL [7.0 mmol/L] and/or use of glucose-lowering drugs) (6), 44 because of unavailability of samples for measurement of PCSK9, and 4 because of being a recipient of a combined kidney/liver transplant, resulting in 453 RTRs eligible for analyses. None of these 453 RTRs was a recipient of a combined kidney/pancreas transplant. Written informed consent was obtained from all participants, and the institutional review board of the University of Groningen (Groningen, the Netherlands) approved the study (METc 01/039). Relevant donor, recipient, and transplant characteristics were extracted from the Groningen Renal Transplant Database, which has been described in detail elsewhere (16). The primary end point of this study was NODAT. Secondary end points were all-cause mortality, cardiovascular mortality, and renal graft failure. Follow-up was performed for a median of 9.6 (IQR 7.1–10.2) years until April 2012. There was no loss to follow-up for the primary and secondary end points. Collection of these data are ensured by the continuous surveillance system of the outpatient clinic of our university hospital and close collaboration with affiliated hospitals.
Measurements and Definitions
Blood pressure was measured as the average of three automated (Omron M4; Omron Europe B.V., Hoofddorp, the Netherlands) measurements with 1-min intervals after a 6-min rest in supine. BMI was calculated as weight in kilograms divided by height in meters squared. Waist circumference was measured on bare skin midway between the iliac crest and the 10th rib.
NODAT was defined according to International Expert Panel recommendations that were based on 2003 American Diabetes Association criteria (2) (i.e., at the end of the time frame of inclusion). The diagnosis of NODAT was based on one of the following criteria: symptoms of diabetes (classic symptoms, including polyuria, polydipsia, and unexplained weight loss) plus nonfasting plasma glucose concentration ≥200 mg/dL (11.1 mmol/L), fasting plasma glucose ≥126 mg/dL (7.0 mmol/L), or use of glucose-lowering medication. Fasting was defined as no caloric intake for at least 8 h. If patients transplanted in our center are >1 year posttransplantation, as is the case for RTRs included in this study, routine follow-up consists of three monthly visits to the outpatient clinic, with check-ups including assessment of renal function and fasting or random glucose. If glucose was elevated, a confirmatory laboratory test of venous plasma was performed on a subsequent day or at the next visit, after which the diagnosis of NODAT was made. Apart from NODAT, we categorized RTRs as having impaired fasting glucose if fasting plasma glucose was ≥100 mg/dL (5.6 mmol/L) and <126 mg/dL (7.0 mmol/L) (17). The secondary end points of the study were all-cause mortality, cardiovascular mortality, and death-censored graft failure. Cause of death was obtained by linking the number of the death certificate to the primary cause of death as coded by a physician from the Central Bureau of Statistics according to the ICD-9 (18). Cardiovascular death was defined as the principal cause of death being cardiovascular in nature (ICD-9 codes 410–447). Graft failure was defined as restart of dialysis or retransplantation. Current medication was extracted from the medical records. Information about change of immunosuppressive medication over time in our transplant center is described in the Supplementary Data. Cumulative dose of prednisolone was calculated by multiplying time after transplantation by prednisolone dose at time of inclusion in the study and addition of the dose of prednisolone or methylprednisolone required for treatment of acute rejection (a conversion factor of 1.25 was used to convert methylprednisolone to prednisolone dose).
Laboratory Measurements
Venous blood was drawn in the morning after an 8–12-h overnight fasting period to determine serum creatinine and plasma glucose concentrations. Serum creatinine was determined by using the Jaffe reaction (MEGA AU510; Merck Diagnostica, Darmstadt, Germany). Plasma glucose was measured by the glucose oxidase method (YSI 2300 Stat Plus; Yellow Springs Instruments, Yellow Springs, OH). Estimated glomerular filtration rate (eGFR) was calculated by applying the Chronic Kidney Disease Epidemiology Collaboration equation (19). Proteinuria was defined as ≥0.5 g protein/24-h urine. Total cholesterol was determined by using the cholesterol oxidase-phenol aminophenazone method (MEGA AU510). LDL cholesterol was calculated by using the Friedewald equation (20). HDL cholesterol was measured with the cholesterol oxidase-phenol aminophenazone method on a Technicon RA-1000 (Bayer Diagnostics, Mijdrecht, the Netherlands). Plasma triglycerides were determined with the glycerol-3-phosphate oxidase-oxidase method (YSI 2300 Stat Plus). Fasting insulin levels at baseline were determined for research purposes by using an AxSYM autoanalyzer (Abbott Diagnostics). Insulin resistance was calculated by using the HOMA for insulin resistance (HOMA-IR) equation: glucose (mmol/L) × insulin (mU/L) / 22.5 (21).
Samples for measurement of PCSK9 were collected at the moment of inclusion in the study, at the same time that baseline characteristics were assessed. The samples for measurement of PCSK9 were kept frozen at −80°C until assessment. Serum PCSK9 was assessed by PCSK9 dual monoclonal antibody sandwich ELISA, with minor modifications (22) (Supplementary Data).
Statistical Analyses
Data were analyzed with SPSS version 22.0 software (IBM Corporation, Chicago, IL). In all analyses, a two-sided P < 0.05 was considered significant. Data are expressed as mean ± SD for normally distributed variables and as median (IQR) for variables with a skewed distribution. Linear regression analysis was performed to determine whether patient characteristics varied across tertiles of PCSK9 by calculating P values for trend. Residuals were checked for normality and log-transformed when appropriate. χ2 test was used for categorical variables. NODAT development was visualized by Kaplan-Meier curves according to the highest versus the two lowest tertiles of PCSK9 levels, with statistical significance tested by log-rank (Mantel-Cox) test.
To study whether PCSK9 was associated with risk of NODAT, all-cause mortality, cardiovascular mortality, and renal graft failure, Cox proportional hazards regression analyses were performed. We first performed crude analyses (model 1) followed by adjustment for age and sex (model 2) and further cumulative adjustment for eGFR, proteinuria, and time since transplantation (model 3). Because Cox regression models should not contain more than one independent variable for every 10 outcome events (23), model 3 was further adjusted for potential confounders in additional Cox regression models. Adjustment variables were selected on the basis of the literature (age, sex, eGFR, proteinuria, and time since transplantation), a significant association with PCSK9 (Table 1), or a significant association with development of NODAT (Supplementary Table 1). If closely related variables within a biological domain (e.g., BMI and waist circumference in the domain of body composition) were associated with PCSK9 and/or NODAT, the variable with the strongest association or the biologically most relevant variable was primarily adjusted for to prevent occurrence of collinearity and overfitting of models. In primary analyses, therefore, we additionally adjusted for total cholesterol, triglycerides, and BMI in model 4; systolic blood pressure and fasting plasma glucose in model 5; smoking and alcohol use in model 6; statin use in model 7; HOMA-IR in model 8; and use of tacrolimus, use of cyclosporine, and trough levels of both in model 9. In the Cox regression analyses, serum PCSK9 was used both as a continuous variable and as a categorical variable in which tertiles 1 and 2 were combined and compared with tertile 3. NODAT and graft failure were censored at the date of the last follow-up or death. As secondary analyses, we repeated the Cox regression analyses for the association of serum PCSK9 with NODAT, with adjustment for the variables in biological domains that were associated with PCSK9 or development of NODAT but that were not adjusted for in the primary analyses (Supplementary Table 2).
Baseline characteristics over the tertiles of PCSK9 (n = 453)
. | Tertiles of sex-stratified PCSK9 (ng/mL) . | . | ||
---|---|---|---|---|
Variable . | T1 . | T2 . | T3 . | P value for trend . |
Serum PCSK9 (ng/mL) | 66.0 ± 16.9 | 102.3 ± 11.5 | 151.6 ± 39.0 | NA |
General characteristics | ||||
Age (years) | 49 ± 12 | 50 ± 12 | 53 ± 12 | 0.003 |
Male sex (%) | 55 | 56 | 56 | NA |
Smoking status | 0.09 | |||
Never smoker (%) | 41 | 35 | 28 | |
Former smoker (%) | 41 | 40 | 48 | |
Current smoker (%) | 18 | 26 | 24 | |
Alcohol consumption | 0.20 | |||
None (%) | 47 | 46 | 44 | |
1–7 units/week (%) | 36 | 36 | 46 | |
>7 units/week (%) | 15 | 18 | 10 | |
Body composition | ||||
BMI (kg/m2) | 25.1 ± 4.1 | 25.8 ± 4.1 | 26.5 ± 4.4 | <0.001 |
Waist circumference (cm) | 93.6 ± 13.0 | 95.0 ± 12.9 | 99.1 ± 14.0 | 0.01 |
Transplantation history | ||||
Time since renal transplantation (years) | 6.2 (2.7–11.3) | 6.1 (3.4–12.2) | 5.9 (2.3–12.0) | 0.70 |
Deceased donor (%) | 84 | 82 | 90 | 0.11 |
Dialysis duration (months) | 30 (11–51) | 27 (12–45) | 31 (18–53) | 0.15 |
Preemptive renal transplantation (%) | 10 | 11 | 4 | 0.40 |
Preoperative positive hepatitis C (%)* | 1 | 0 | 2 | 0.76 |
Renal allograft function | ||||
Serum creatinine (μmol/L) | 138 (119–167) | 138 (113–159) | 134 (111–171) | 0.009 |
eGFR (mL/min/1.73 m2) | 46.8 ± 16.2 | 46.5 ± 14.8 | 47.3 ± 16.3 | 0.78 |
Proteinuria (≥0.5 g/24 h) (%) | 30 | 28 | 24 | 0.20 |
Blood pressure | ||||
Diastolic blood pressure (mmHg) | 91 ± 10 | 89 ± 10 | 90 ± 10 | 0.37 |
Systolic blood pressure (mmHg) | 152 ± 23 | 150 ± 23 | 152 ± 22 | 0.53 |
Glucose homeostasis | ||||
Plasma glucose (mmol/L) | 4.5 ± 0.6 | 4.5 ± 0.7 | 4.6 ± 0.7 | 0.12 |
Plasma insulin (μmol/L) | 9.4 (7.3–12.5) | 10.7 (8.0–14.9) | 11.2 (7.7–16.4) | 0.01 |
HOMA-IR [mU × mmol/(L2 × 22.5)] | 1.8 (1.3–2.6) | 2.1 (1.6–3.0) | 2.2 (1.6–3.3) | 0.001 |
Impaired fasting glucose (%) | 5 | 7 | 7 | 0.16 |
Family history of diabetes: parent or sibling with diabetes, n (%) | 35 (24) | 41 (27) | 42 (27) | 0.25 |
Lipids and lipoproteins | ||||
Total cholesterol (mmol/L) | 5.6 ± 0.9 | 5.6 ± 0.9 | 5.8 ± 1.3 | 0.21 |
LDL cholesterol (mmol/L) | 3.6 ± 0.8 | 3.5 ± 0.9 | 3.6 ± 1.2 | 0.17 |
Non-LDL cholesterol (mmol/L) | 4.5 ± 0.9 | 4.4 ± 0.9 | 4.7 ± 1.3 | 0.49 |
HDL cholesterol (mmol/L) | 1.1 ± 0.3 | 1.1 ± 0.3 | 1.2 ± 0.3 | 0.05 |
Triglycerides (mmol/L) | 1.7 (1.2–2.4) | 1.9 (1.3–2.6) | 2.0 (1.6–2.7) | <0.001 |
Medication use | ||||
Blood pressure lowering (%) | 81 | 91 | 85 | 0.73 |
Statins (%) | 31 | 49 | 65 | <0.001 |
Proliferation inhibitor (%) | 78 | 75 | 69 | 0.28 |
Calcineurin inhibitor (%) | 77 | 79 | 78 | 0.54 |
Tacrolimus (%) | 18 | 18 | 6 | 0.03 |
Tacrolimus (trough level, μg/L) | 7.8 (5.8–10.1) | 8.9 (7.4–10.7) | 8.0 (4.0–10.3) | 0.52 |
Cyclosporine (%) | 60 | 61 | 72 | 0.30 |
Cyclosporine (trough level, μg/L) | 97.5 (73.3–131.0) | 103.0 (78.0–136.0) | 118.0 (84.0–152.5) | 0.01 |
Prednisolone (mg/24 h) | 9.0 ± 1.5 | 9.3 ± 1.2 | 9.3 ± 1.2 | 0.41 |
Cumulative prednisolone dose (g) | 21.5 (11.5–37.3) | 22.1 (12.9–38.5) | 21.2 (10.5–39.0) | 0.64 |
. | Tertiles of sex-stratified PCSK9 (ng/mL) . | . | ||
---|---|---|---|---|
Variable . | T1 . | T2 . | T3 . | P value for trend . |
Serum PCSK9 (ng/mL) | 66.0 ± 16.9 | 102.3 ± 11.5 | 151.6 ± 39.0 | NA |
General characteristics | ||||
Age (years) | 49 ± 12 | 50 ± 12 | 53 ± 12 | 0.003 |
Male sex (%) | 55 | 56 | 56 | NA |
Smoking status | 0.09 | |||
Never smoker (%) | 41 | 35 | 28 | |
Former smoker (%) | 41 | 40 | 48 | |
Current smoker (%) | 18 | 26 | 24 | |
Alcohol consumption | 0.20 | |||
None (%) | 47 | 46 | 44 | |
1–7 units/week (%) | 36 | 36 | 46 | |
>7 units/week (%) | 15 | 18 | 10 | |
Body composition | ||||
BMI (kg/m2) | 25.1 ± 4.1 | 25.8 ± 4.1 | 26.5 ± 4.4 | <0.001 |
Waist circumference (cm) | 93.6 ± 13.0 | 95.0 ± 12.9 | 99.1 ± 14.0 | 0.01 |
Transplantation history | ||||
Time since renal transplantation (years) | 6.2 (2.7–11.3) | 6.1 (3.4–12.2) | 5.9 (2.3–12.0) | 0.70 |
Deceased donor (%) | 84 | 82 | 90 | 0.11 |
Dialysis duration (months) | 30 (11–51) | 27 (12–45) | 31 (18–53) | 0.15 |
Preemptive renal transplantation (%) | 10 | 11 | 4 | 0.40 |
Preoperative positive hepatitis C (%)* | 1 | 0 | 2 | 0.76 |
Renal allograft function | ||||
Serum creatinine (μmol/L) | 138 (119–167) | 138 (113–159) | 134 (111–171) | 0.009 |
eGFR (mL/min/1.73 m2) | 46.8 ± 16.2 | 46.5 ± 14.8 | 47.3 ± 16.3 | 0.78 |
Proteinuria (≥0.5 g/24 h) (%) | 30 | 28 | 24 | 0.20 |
Blood pressure | ||||
Diastolic blood pressure (mmHg) | 91 ± 10 | 89 ± 10 | 90 ± 10 | 0.37 |
Systolic blood pressure (mmHg) | 152 ± 23 | 150 ± 23 | 152 ± 22 | 0.53 |
Glucose homeostasis | ||||
Plasma glucose (mmol/L) | 4.5 ± 0.6 | 4.5 ± 0.7 | 4.6 ± 0.7 | 0.12 |
Plasma insulin (μmol/L) | 9.4 (7.3–12.5) | 10.7 (8.0–14.9) | 11.2 (7.7–16.4) | 0.01 |
HOMA-IR [mU × mmol/(L2 × 22.5)] | 1.8 (1.3–2.6) | 2.1 (1.6–3.0) | 2.2 (1.6–3.3) | 0.001 |
Impaired fasting glucose (%) | 5 | 7 | 7 | 0.16 |
Family history of diabetes: parent or sibling with diabetes, n (%) | 35 (24) | 41 (27) | 42 (27) | 0.25 |
Lipids and lipoproteins | ||||
Total cholesterol (mmol/L) | 5.6 ± 0.9 | 5.6 ± 0.9 | 5.8 ± 1.3 | 0.21 |
LDL cholesterol (mmol/L) | 3.6 ± 0.8 | 3.5 ± 0.9 | 3.6 ± 1.2 | 0.17 |
Non-LDL cholesterol (mmol/L) | 4.5 ± 0.9 | 4.4 ± 0.9 | 4.7 ± 1.3 | 0.49 |
HDL cholesterol (mmol/L) | 1.1 ± 0.3 | 1.1 ± 0.3 | 1.2 ± 0.3 | 0.05 |
Triglycerides (mmol/L) | 1.7 (1.2–2.4) | 1.9 (1.3–2.6) | 2.0 (1.6–2.7) | <0.001 |
Medication use | ||||
Blood pressure lowering (%) | 81 | 91 | 85 | 0.73 |
Statins (%) | 31 | 49 | 65 | <0.001 |
Proliferation inhibitor (%) | 78 | 75 | 69 | 0.28 |
Calcineurin inhibitor (%) | 77 | 79 | 78 | 0.54 |
Tacrolimus (%) | 18 | 18 | 6 | 0.03 |
Tacrolimus (trough level, μg/L) | 7.8 (5.8–10.1) | 8.9 (7.4–10.7) | 8.0 (4.0–10.3) | 0.52 |
Cyclosporine (%) | 60 | 61 | 72 | 0.30 |
Cyclosporine (trough level, μg/L) | 97.5 (73.3–131.0) | 103.0 (78.0–136.0) | 118.0 (84.0–152.5) | 0.01 |
Prednisolone (mg/24 h) | 9.0 ± 1.5 | 9.3 ± 1.2 | 9.3 ± 1.2 | 0.41 |
Cumulative prednisolone dose (g) | 21.5 (11.5–37.3) | 22.1 (12.9–38.5) | 21.2 (10.5–39.0) | 0.64 |
Data are mean ± SD and median (IQR) unless otherwise indicated. Dialysis duration (months), serum creatinine, LDL cholesterol, and serum triglycerides were log-transformed as a result of skewed distribution. NA, not applicable.
*Data on hepatitis C status were available for 274 RTRs. Cut-off values in the three groups were T1 <81.2 and <95.8 ng/mL, T2 81.2–111.2 and 95.8–127.2 ng/mL, and T3 >111.2 and >127.2 ng/mL for males and females, respectively.
To allow for comparison with studies in which time-to-event data were not available, we proceeded with sensitivity analyses in which we applied logistic regression rather than Cox regression. We also performed sensitivity analyses in which we investigated potential effect modification of the association of PCSK9 with NODAT by use of statins.
Results
Baseline Characteristics
We included 453 stable RTRs. Mean age was 51 ± 12 years, 56% were male, and BMI averaged 25.8 ± 4.2 kg/m2. Patients were included at 6.1 (2.7–11.7) years after transplantation. The mean serum PCSK9 concentration was 107.1 ± 43.4 ng/mL. Baseline characteristics of the RTRs according to tertiles of PCSK9 are shown in Table 1. Plasma PCSK9 levels were positively associated with age, BMI, waist circumference, plasma glucose, plasma insulin, HOMA-IR, serum triglycerides, and statin use. Serum PCSK9 was associated with neither total cholesterol nor LDL cholesterol and eGFR. Serum PCSK9 concentrations in RTRs treated with statins were 119.3 ± 45.7 versus 95.6 ± 37.7 ng/mL in RTRs not treated with statins (P < 0.001). Twenty nine (6%) RTRs had impaired fasting glucose at baseline.
Prospective Analyses of NODAT
During a median follow-up of 9.6 (7.1–10.2) years, 70 (15.5%) RTRs developed NODAT. NODAT developed in 23% (n = 35) in the upper tertile of PCSK9 versus 12% (n = 35) in the lowest two tertiles (P < 0.001 by log-rank test). The corresponding Kaplan-Meier curves are shown in Fig. 1. In Cox regression analysis, serum PCSK9 as a continuous variable was univariately associated with development of NODAT (hazard ratio [HR] 1.34 [95% CI 1.10–1.63]; P = 0.004) (Table 2). After adjustment for relevant covariates (model 2: age and sex; model 3: age, sex, eGFR, proteinuria, and time since transplantation [Fig. 2]; and subsequently added to model 3, model 4: serum triglycerides, total cholesterol, and BMI; model 5: systolic blood pressure and fasting plasma glucose; model 6: smoking and alcohol use; model 7: statin use; model 8: HOMA-IR; and model 9: use of tacrolimus, use of cyclosporine, and trough levels of both), the association of PCSK9 with NODAT remained materially similar (HR 1.39 [95% CI 1.13–1.72]; P = 0.002) after adjustment for age, sex, eGFR, proteinuria, time since transplantation, and use of tacrolimus, use of cyclosporine, and trough levels of both. In an additional Cox regression analysis with PCSK9 divided in the upper versus the two lower tertiles, PCSK9 was again significantly associated with development of NODAT independent of potential confounders (HR 2.37 [95% CI 1.45–3.87]; P = 0.001) after adjustment for age, sex, eGFR, proteinuria, time since transplantation, and use of tacrolimus, use of cyclosporine, and trough levels of both. Secondary analyses in which adjustment for closely related covariates within biological domains were exchanged for adjustment by the variable not included in the primary analyses are shown in Supplementary Table 2. Results of secondary analyses, in which we repeated the Cox regression analyses with inclusion of the variables in biological domains that were not included, did not differ materially from the results of the primary analyses (Supplementary Table 2).
Kaplan-Meier curves for NODAT according to the highest tertile versus the lower two tertiles of PCSK9 in RTRs.
Kaplan-Meier curves for NODAT according to the highest tertile versus the lower two tertiles of PCSK9 in RTRs.
Serum PCSK9 and NODAT in RTRs (n = 453)
. | Serum PCSK9 per SD, as continuous variable . | Highest tertile PCSK9 vs. lower two tertiles . | |||
---|---|---|---|---|---|
NODAT . | HR (95% CI) . | P value . | Reference . | HR (95% CI) . | P value . |
Model 1 | 1.34 (1.10–1.63) | 0.004 | 1.00 | 2.23 (1.40–3.56) | 0.001 |
Model 2 | 1.38 (1.12–1.70) | 0.002 | 1.00 | 2.33 (1.44–3.76) | 0.001 |
Model 3 | 1.38 (1.12–1.70) | 0.002 | 1.00 | 2.32 (1.43–3.76) | 0.001 |
Model 4 | 1.32 (1.07–1.62) | 0.01 | 1.00 | 2.11 (1.29–3.46) | 0.003 |
Model 5 | 1.35 (1.09–1.66) | 0.006 | 1.00 | 2.44 (1.48–4.01) | <0.001 |
Model 6 | 1.36 (1.11–1.67) | 0.003 | 1.00 | 2.32 (1.41–3.81) | 0.001 |
Model 7 | 1.31 (1.05–1.63) | 0.02 | 1.00 | 2.10 (1.28–3.42) | 0.003 |
Model 8 | 1.28 (1.05–1.56) | 0.01 | 1.00 | 2.18 (1.35–3.54) | 0.002 |
Model 9 | 1.39 (1.13–1.72) | 0.002 | 1.00 | 2.37 (1.45–3.87) | 0.001 |
. | Serum PCSK9 per SD, as continuous variable . | Highest tertile PCSK9 vs. lower two tertiles . | |||
---|---|---|---|---|---|
NODAT . | HR (95% CI) . | P value . | Reference . | HR (95% CI) . | P value . |
Model 1 | 1.34 (1.10–1.63) | 0.004 | 1.00 | 2.23 (1.40–3.56) | 0.001 |
Model 2 | 1.38 (1.12–1.70) | 0.002 | 1.00 | 2.33 (1.44–3.76) | 0.001 |
Model 3 | 1.38 (1.12–1.70) | 0.002 | 1.00 | 2.32 (1.43–3.76) | 0.001 |
Model 4 | 1.32 (1.07–1.62) | 0.01 | 1.00 | 2.11 (1.29–3.46) | 0.003 |
Model 5 | 1.35 (1.09–1.66) | 0.006 | 1.00 | 2.44 (1.48–4.01) | <0.001 |
Model 6 | 1.36 (1.11–1.67) | 0.003 | 1.00 | 2.32 (1.41–3.81) | 0.001 |
Model 7 | 1.31 (1.05–1.63) | 0.02 | 1.00 | 2.10 (1.28–3.42) | 0.003 |
Model 8 | 1.28 (1.05–1.56) | 0.01 | 1.00 | 2.18 (1.35–3.54) | 0.002 |
Model 9 | 1.39 (1.13–1.72) | 0.002 | 1.00 | 2.37 (1.45–3.87) | 0.001 |
In total, 70 RTRs developed NODAT (lower two tertiles 35, highest tertile 35). Model 1: crude analysis. Model 2: model 1 + adjustment for age and sex. Model 3: model 2 + adjustment for eGFR, proteinuria, and time since transplantation. Model 4: model 3 + adjustment for triglycerides, total cholesterol, and BMI. Model 5: model 3 + adjustment for systolic blood pressure and plasma glucose. Model 6: model 3 + adjustment for smoking and alcohol use. Model 7: model 3 + adjustment for statin use. Model 8: model 3 + adjustment for HOMA-IR. Model 9: model 3 + adjustment for use of tacrolimus, use of cyclosporine, and trough levels of both.
Associations between PCSK9 and NODAT. Data were fit by a Cox proportional hazards regression model that was based on restricted cubic splines and adjusted for age, sex, eGFR, time since transplantation, and proteinuria. The reference standard was a median PCSK9 level of 103.6 ng/mL. The solid line represents the HR. The gray area represents the 95% CI.
Associations between PCSK9 and NODAT. Data were fit by a Cox proportional hazards regression model that was based on restricted cubic splines and adjusted for age, sex, eGFR, time since transplantation, and proteinuria. The reference standard was a median PCSK9 level of 103.6 ng/mL. The solid line represents the HR. The gray area represents the 95% CI.
Sensitivity Analyses on PCSK9 and NODAT
To evaluate the association of PCSK9 with NODAT without taking time to event into account, we performed crude and multivariable logistic regression analyses equivalent to Cox regression models 1 and 3. The crude analysis revealed an odds ratio of 1.34 (95% CI 1.06–1.69; P = 0.02), and the analysis equivalent to model 3 revealed an odds ratio of 1.41 (95% CI 1.10–1.81; P = 0.007).
In further sensitivity analyses, we tested for interaction by use of statins. We found no indication of significant effect modification by use of statins if statin use and a product term of statin use and PCSK9 were added to the crude model (P = 0.52 for interaction). The same was true if the product term was added to multivariable Cox regression model 7 with additional adjustment for age, sex, eGFR, proteinuria, time since transplantation, and use of statins (P = 0.65 for interaction).
Prospective Analyses on All-Cause Mortality, Cardiovascular Mortality, and Graft Failure
During a median follow-up of 10.0 (IQR 9.7–10.4) years, 123 RTRs died. In univariable Cox regression analysis, serum PCSK9 as a continuous variable was not significantly associated with increased risk of mortality (HR 1.12 [95% CI 0.95–1.31]; P = 0.17). After adjustment for relevant covariates, there was also no association of PCSK9 with all-cause mortality (Supplementary Table 3). Of the 123 deceased RTRs, cause of death for 51 (41%) was cardiovascular related. There was no significant association between PCSK9 and cardiovascular mortality (Supplementary Table 4).
During a median follow-up of 9.6 (7.1–10.2) years, 59 RTRs developed graft failure. In univariable Cox regression analysis, serum PCSK9 as a continuous variable was not significantly associated with graft failure (HR 0.77 [95% CI 0.57–1.02]; P = 0.07). After adjustment for relevant covariates, there was also no association of PCSK9 with graft failure (Supplementary Table 5).
Conclusions
In a large cohort of stable RTRs, we show that higher serum PCSK9 levels are associated with development of NODAT. This association remains present independent of relevant confounders and independent of adjustment for use of statins. The data agree with the concept that the PCSK9 pathway plays a pathogenic role in the development of diabetes in this patient category. Furthermore, the results provide a rationale for investigating the extent to which PCSK9 inhibitor treatment would prove beneficial with respect to prevention of NODAT in RTRs, besides its well-established LDL cholesterol–lowering effect.
To date, the association of PCSK9 levels with NODAT has been investigated neither in RTRs nor in other transplanted populations. In renal patients, PCSK9 has been documented to be elevated proportional to the degree of proteinuria (24) but to be unaffected by moderately compromised renal function (25). In the current RTR cohort, PCSK9 was increased in conjunction with (central) obesity, as observed in nonrenal populations (26–28). Thus, although we did not enroll a control group, serum PCSK9 levels appeared to be elevated in a considerable number of RTRs. Many reports have revealed positive correlations of circulating PCSK9 with LDL cholesterol (24,26,28), but such a relationship was not evident in the current study probably because PCSK9 increased the effect of statin treatment (24,25), thereby masking a positive PCSK9-LDL cholesterol relationship. This relationship emphasizes the necessity to adjust for statin use in multivariable analysis regarding the association of PCSK9 with NODAT. In keeping with data in nonrenal populations (24,26,28,29), serum PCSK9 was associated with triglycerides, which we also controlled for in the analysis. Moreover, in keeping with earlier reports (26,30), serum PCSK9 was positively associated with HOMA-IR, which is a known predictor for NODAT (31). Thus, the association of PCSK9 with NODAT as documented here was essentially unaltered, taking into account the use of statins and other covariates, including HOMA-IR. A potential role of statin treatment was also evaluated, but we observed no interaction.
PCSK9 is a main regulator of LDLR availability, and the findings regarding the PCSK9-NODAT association agree with the hypothesis that cholesterol metabolism is linked to diabetes development by involving the LDLR pathway. PCSK9 was positively related to insulin and HOMA-IR, again in keeping with other data (26). The precise mechanisms by which the PCSK9 pathway are involved in diabetes development need to be delineated in more detail in future studies. In comparison, 6–18 months of treatment with the PCSK9 inhibitor alirocumab did not affect transition to new-onset diabetes in populations at high cardiovascular risk, but longer follow-up is required to more conclusively rule out any association (32).
Although evidence is mounting in favor of a beneficial effect of PCSK9 administration on incident CVD (33–35), remarkably few studies have addressed whether circulating PCSK9 predicts cardiovascular outcome. In a Swedish cohort study, PCSK9 was found to independently predict incident CVD during 15 years of follow-up (36). In contrast, PCSK9 was unrelated to newly developed CVD during 17 years of follow-up in the Women’s Health Study (27). PCSK9 was not significantly associated with CVD events or with renal function decline during 3–10 years of follow-up in two German cohorts comprising subjects with moderately impaired kidney function (37). In the current study, PCSK9 did not predict all-cause and cardiovascular mortality. Given the resemblance of transplant vasculopathy, an important cause of graft failure, with atherosclerosis (38), we also determined whether PCSK9 predicts renal outcome. PCSK9 was not related to graft failure in the current study.
This study has several strengths. To our knowledge, it comprises the largest prospectively followed cohort of RTRs in which PCSK9 was measured. Moreover, end point evaluation was complete in all patients despite a considerable follow-up period. Another strength of the study is that we included only stable RTRs who were >1 year posttransplantation, resulting in exclusion of transient posttransplantation hyperglycemia in NODAT diagnosis. Hyperglycemia is extremely common in the early posttransplant period and can occur as a result of rejection therapy, infections, and other critical conditions. Therefore, the formal diagnosis of NODAT in RTRs should only be based on likely maintenance immunosuppression, stable kidney function, and absent acute infections (39).
The study also has several limitations. It was carried out at a single renal transplant center, so the extent to which circulating PCSK9 predicts NODAT development must be replicated in other RTR populations that include sufficient numbers of nonwhite participants. We did not have repeated PCSK9 measurements, and the use of a single measurement of the variable of interest may have given rise to an underestimation of an otherwise existing true effect (40). Another limitation is that we did not perform an oral glucose tolerance test at baseline or during follow-up, which is considered the gold standard test to define NODAT. Thus, we could have underestimated the frequency of NODAT compared with other studies. However, we diagnosed NODAT according to International Expert Panel recommendations that were based on American Diabetes Association criteria (2), which is in line with the common way in which NODAT is diagnosed in transplant centers. Furthermore, as with any observational study, unmeasured or residual confounding may have existed despite the substantial number of potentially confounding factors for which we adjusted. Finally, the observational nature of the study did not allow us to discern whether higher PCSK9 is a cause of NODAT or merely a marker of a higher risk for NODAT.
In conclusion, higher serum PCSK9 is associated with an increased risk of NODAT in RTRs independent of potential confounders. Inhibition of PCSK9 may be a promising target after transplantation to diminish the risk of developing NODAT while having a beneficial effect on hypercholesterolemia.
Article Information
Acknowledgments. This study was based on the TransplantLines Insulin Resistance and Inflammation (TxL-IRI) cohort study. PCSK9 was measured at Lilly Research Laboratories, Eli Lilly and Company (Indianapolis, IN) for free.
Duality of Interest. No potential conflicts of interest relevant to this article were reported.
Author Contributions. M.F.E., D.M.Z., S.J.L.B., and R.P.F.D. wrote the manuscript, analyzed data, and performed the statistical analysis. J.H.S. and S.J.L.B. acquired the data and performed the research. C.A.J.M.G., S.J.L.B., and R.P.F.D. provided critical review, advice, and consultation throughout. M.F.E., D.M.Z., J.H.S., C.A.J.M.G., S.J.L.B., and R.P.F.D. read and approved the final manuscript. R.P.F.D. is the guarantor of this work and, as such, had full access to all the data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis.