- Research article
- Open Access
Adverse outcomes associated with rapid linear and non-linear patterns of chronic kidney disease progression
BMC Nephrology volume 22, Article number: 82 (2021)
Patients with rapidly declining renal function face the dual threat of end-stage renal disease (ESRD) and mortality prior to ESRD. What is less well characterised is whether the pattern of the renal trajectory, linear or non-linear, unmasks subgroups of rapidly progressing patients that face adverse outcomes in a differential manner.
An individual eGFR slope was applied to all outpatient estimated glomerular filtration rate (eGFR) values for each patient in the Salford Kidney Study from 2002 to 2018 who had at least 2 years follow-up, ≥4 eGFR values and baseline eGFR 15 to < 60 ml/min/1.73m2. Rapid progression was defined as an annual eGFR slope of ≤ − 3 ml/min/1.73m2/yr and patients were categorised as linear or non-linear progressors based on the nature of their eGFR-time graphs. A Fine-Gray competing risk hazard model was used to determine factors associated with progression to ESRD and with mortality prior to ESRD. Cumulative incidence function curves highlighted differences in outcomes between linear and non-linear patients.
There were 211 rapidly deteriorating patients with linear eGFR trajectories and 61 rapid non-linear patients in the study cohort. Factors associated with ESRD included younger age, male gender, lower baseline eGFR and higher serum phosphate, whilst older age, history of myocardial infarction and anaemia predicted mortality prior to ESRD. Over a median follow-up of 3.7 years, linear progressors reached ESRD sooner whilst those with non-linear progression faced significantly higher rates of mortality prior to ESRD.
Patients with rapid eGFR decline have high rates of adverse outcomes that are differentially expressed in those progressing linearly and non-linearly as a result of differing phenotypic profiles. Consequently, addressing individual risk factor profiles is important to deliver optimal personalised patient care.
Patients with chronic kidney disease (CKD) that experience rapidly declining renal function are at increased risk of adverse outcomes, including increased risk of end-stage renal disease (ESRD) as well as increased risk of mortality prior to ESRD . Over the past decade, there has been growing attention to not only considering the slope of renal function change, as defined by annual changes in estimated glomerular filtration rate (eGFR) , but also the pattern of the eGFR trajectory on patient outcomes [3, 4], especially given that patients with CKD progress in ways other than in a simply linear manner .
Extending our understanding of how the eGFR slope and trajectory impacts patient outcomes could help to refine current risk prediction tools to stratify high-risk patients more accurately. This would help improve the communication of risk to patients and help shape management strategies, including earlier targeted treatment in an attempt to assuage future harm .
Whilst several studies have recognised the determinants of rapid progression, what is less well known is whether the pattern of rapid CKD progression, be it linear or non-linear, has an impact on patient outcomes. We therefore undertook this study to 1) identify the predictive factors of rapid progression in a cohort of patients progressing in a linear and non-linear pattern; 2) to examine how the pattern of rapid progression affects outcomes of ESRD and mortality prior to ESRD, and in doing so, 3) identify whether there are subgroups and phenotypic differences between linear and non-linear progressors that could enlighten specific approaches to patient management.
Patients were drawn from the Salford Kidney Study (SKS), an ongoing observational cohort study, which since 2002 has been recruiting patients aged ≥18 years old with non-dialysis CKD who have been referred to the renal services at Salford Royal NHS Foundation Trust in the United Kingdom. The SKS received ethical approval from the North West Greater Manchester South Research Ethics Committee (REC15/NW/0818). Written informed consent was obtained from all patients. The methods described herein were carried out in accordance with relevant guidelines and regulations of the SKS.
Patient characteristics were measured at the point of recruitment into the SKS. They include patient demographics (age, gender, ethnicity, history of past or current smoking); past medical history (primary renal disease, history of hypertension (HTN), diabetes mellitus (DM), myocardial infarction (MI), peripheral vascular disease (PVD), stroke and heart failure (HF)); medication history (use of angiotensin converting enzyme inhibitor (ACEi), angiotensin receptor blockers (ARB) and statins); and laboratory measurements (serum creatinine, eGFR calculated using the CKD epidemiology collaboration (CKD-EPI) equation, bicarbonate, urea, calcium, phosphate, albumin, haemoglobin and urine protein:creatinine ratio (uPCR), which was categorised into albuminuria grades of A1, A2 and A3 based on values of < 15 g/mol, 15-50 g/mol and > 50 g/mol respectively) .
Assembling the study cohort
Patient selection into this study is shown in Fig. 1. All the outpatient eGFR values, performed as part of routine renal care and accessed via the hospital’s electronic patient record, were used to calculate the delta (Δ) eGFR slope for each patient using linear regression. As a minimum, we required at least 4 eGFR values over 2 years follow-up for the ΔeGFR to be ascertained for each patient. Patients with a ΔeGFR≤ − 3 ml/min/1.73m2/yr (ie. losing more than 3ml/min/1.73m2/yr), a threshold associated with worse outcomes , were defined as a rapid progressor, and they had to have baseline CKD G3a-4 (eGFR 15 to < 60 ml/min/1.73m2) for study inclusion. To differentiate linear versus non-linear progression, the eGFR-time graphs were visually inspected independently by two clinicians, an approach that has been successfully utilised in previous studies [1, 8], and also quantitively assessed with the coefficient of determination (R2).
Patient outcomes included ESRD (haemodialysis, peritoneal dialysis, conservative care or pre-emptive transplantation) or mortality prior to ESRD. Outcome events were reviewed until 1st January 2020.
Continuous data is presented as median (± interquartile range) and categorical data as number (percentage). To compare variables between rapid linear and rapid non-linear progressors, Mann-Whitney U test was used for continuous data and chi-squared test for categorical covariates. The Fine-Gray competing risk hazard model  was employed to determine the sub-distribution hazard ratios of the factors associated with ESRD or mortality prior to ESRD within the study cohort. The following 20 variables were included in the model: age, gender, SBP, DBP, HTN, DM, smoking, MI, PVD, stroke, HF, ACEi/ARB use, statin use, eGFR, bicarbonate, calcium, phosphate, albumin, Hb and A3 proteinuria. The proportional hazards assumption for each model was assessed by the non-significance of each time-by-variable interaction.
Cumulative incidence function curves were produced comparing the outcomes of ESRD and mortality between linear and non-linear progressors, for which a modified Chi-squared test was used for significance testing . A comparison between linear and non-linear progressors on a composite outcome of either ESRD or mortality was also assessed and visualised as a 1-Kaplan-Meier curve, which used log-rank significance testing. Statistical significance was defined as a p value of < 0.05. Analyses were performed using SPSS (Version 25.0) (IBM SPSS, Chicago, IL) licensed to the University of Manchester and R version 4.0.2 (The R Foundation for Statistical Computing Platform).
A total of 272 patients met the inclusion criteria of which 211 patients had linear progression and 61 patients had non-linear progression (Table 1). Categorisation as a linear or non-linear progressor was achieved with unanimous agreement between the two clinicians independently reviewing patients’ eGFR-time graphs, illustrative examples of which are shown in Fig. 2. Quantitatively, linear patients had a significantly higher median R2 value of 0.91 (0.81–0.95) compared with 0.58 (0.36–0.60) in non-linear progressors, p-value < 0.01.
There were no significant differences in laboratory measures between the two patient groups. Of note, both groups demonstrated poor baseline renal function with a median eGFR of 34 ml/min/1.73m2 (26-41 ml/min/1.73m2) and 31 ml/min/1.73m2 (23-41 ml/min/1.73m2) in linear and non-linear progressors respectively. Both groups also had high degrees of proteinuria, which was reflected in the majority of patients classified with A3 proteinuria.
The categorisation of the nature of the ΔeGFR slope in each patient group was strengthened by a large number of eGFR measurements per patient during the follow-up period with a median of 24 (16–36) in the linear group and 21 (11–37) in the non-linear group. The ΔeGFR itself met the a priori definition of rapid renal decline in both patient groups, and were statistically similar: in linear progressors, the median ΔeGFR was − 5.28 ml/min/1.73m2/yr (− 6.75 to − 4.11 ml/min/1.73m2/yr) and in non-linear progressors the ΔeGFR was − 4.57 ml/min/1.73m2/yr (− 6.46 to − 3.46 ml/min/1.73m2/yr); p-value of 0.06.
There were however significant differences between the progressor groups with respect to demographic characteristics and co-morbidities. Non-linear patients were typically older and had a higher burden of co-morbidities including diabetes, myocardial infarction, peripheral vascular disease and heart failure. There was also contrasting frequencies of the underlying disease aetiology, most notably seen for autosomal dominant polycystic kidney disease (ADPKD), which was exclusively associated with linear progressors and was the commonest primary renal disease in this group, affecting a quarter of all linear progressor patients.
Factors associated with ESRD and mortality prior to ESRD
Univariate analyses of the factors associated with ESRD is presented in Additional file 1: Tables S1. In multivariate analysis, younger age, male gender, lack of diabetes, lower eGFR and higher phosphate were significantly associated with progression to ESRD (Table 2).
Three clinical factors were shown to be significantly predictive for mortality prior to ESRD, including older age, history of MI and anaemia (Table 2). Of note, in univariate analyses (Additional file 1: Tables S2), higher SBP, history of diabetes, smoking, PVD, HF and lack of ACEi/ARB use were associated with mortality prior to ESRD but these were found to not be significant once adjusted for other variables.
Survival analysis comparing linear and non-linear progressors
The study cohort had a median follow-up of 3.7 years (2.8–4.9 years), during which time there were 173 patients who reached ESRD and 81 patients who died prior to ESRD (Table 3). Cumulative incidence function curves (Fig. 3) revealed that rapid linear patients were more likely to develop ESRD compared to their non-linear counterparts (p-value< 0.01), but that non-linear progressors suffered higher rates of death prior to ESRD (p-value< 0.01). However, when the outcomes were combined as a single endpoint of either ESRD or mortality prior to ESRD (Fig. 4), there was no statistical difference between rapid linear or non-linear progressors (p-value 0.08).
This study highlights the differential impact of the pattern of renal trajectory on adverse outcomes in patients with rapidly deteriorating renal function. Patients with rapid linear progression face higher rates of ESRD sooner whilst those with non-linear progression experience higher rates of mortality prior to ESRD.
Patterns of progression and determinants of adverse outcomes
Our cohort of rapid progressors was predominantly comprised of patients with linear progression, with only 22% of patients deemed to have non-linear progression. This corroborates with work by Weldegiorgis et al. , who showed that the majority of the 3523 pooled patients with CKD from six randomised trials demonstrated linear eGFR decline. For instance, they reported a > 50% probability of non-linear progression in 15.1 to 21.2% of patients from non-diabetic kidney trials and 19.3 to 31.7% in diabetic kidney trials.
With respect to the determinants of adverse outcomes, previous work has shown that younger age , male gender [13, 14], lower eGFR  and higher phosphate  are all independently associated with ESRD, and this was borne out in our analysis. In addition, we found that older age , a history of MI [16,17,18], and anaemia  are significant determinants of mortality prior to ESRD, and these findings again align with the reported literature. Thus, in our cohort of clearly characterised linear and non-linear progressors, we show the expected overlap of predictive factors important to those who rapidly progress. However, we notably found a lack of diabetes to be predictive for ESRD and also found that A3 proteinuria was not predictive of worse outcomes in our analysis. These findings may likely be attributed to the specific disease characteristics in our population. Indeed, 20% of the entire cohort had ADPKD, of whom 87% reached ESRD. Of the 55 patients with ADPKD only 1 (2%) had diabetes and 8 (15%) were classed as having A3 proteinuria (Additional file 1: Table S3), and this may have affected the impact of diabetes and A3 proteinuria on our outcomes in this cohort.
What is less well known is the association of the pattern of rapid CKD progression, either linear or non-linear, on outcomes such as ESRD and mortality. We hypothesised that rapid linear progressors would fare significantly worse than non-linear progressors, largely as a result of the inevitable culmination of ESRD in linear progressors, whereas potential fluctuating phases of stability may offer a degree of protection to those with non-linear renal decline. We discovered that whilst linear progressors do indeed reach ESRD sooner, non-linear patients were at a significantly greater risk of the competing event of death prior to ESRD.
Phenotypic differences between the two patient groups can help to unravel these findings. For instance, the primary renal disease in 26% of the rapid linear progressors was ADPKD, a condition typically characterised by a linear trajectory to ESRD , as a result of the unremitting enlargement of renal cysts that eventually destroy the normal renal architecture.
Furthermore, we found non-linear patients were significantly older, suffered from diabetes and cardiovascular co-morbidities, including MI, PVD and HF. It is conceivable that such patients experience non-linear fluctuations in renal function, including acute kidney injury (AKI), due to uncorrected alterations in their fluid status, as is seen in decompensated heart failure . Additional adjustments to their medications with up-titration of diuretics or with ACEi and ARBs may also cause transient and abrupt eGFR decline. It is clear that the greater burden of cardiovascular disease and the propensity for outpatient AKI-on-CKD events  explains the increased mortality experienced by non-linear progressors observed in our cohort.
In considering the pattern of CKD progression, we show how two distinct subtypes of rapid progression, linear and non-linear, affect outcomes differentially, which can be explained by differences in patient characteristics. This provides the basis for delivering personalised care to patients, dictated by their disease aetiology, risk factor profile and pattern of CKD progression, and highlights how in some patients the risk of death supersedes the need to prepare for ESRD . This will no doubt influence the communication of risk imparted to patients and help direct the shape of future treatment. Further research in this area is of importance not least because patient heterogeneity for future adverse outcomes is not adequately captured within the current KDIGO staging for CKD, which is based on eGFR and urine albumin:creatinine ratio alone . The future of precision medicine, therefore, relies upon the establishment of improved and refined models that can accurately risk stratify specific CKD patient subgroups.
However, whilst there are specific differences between linear and non-linear progressors with respect to outcomes, there was no survival difference in our analysis between the patients when combining the outcomes as a composite endpoint. Based on this, it is important to emphasise the risk conferred to patients from rapid renal decline per se, irrespective of the underlying pattern of progression. In effect, therefore, early identification of rapid progression should prompt close monitoring and aggressive risk factor modification to stabilise and curb a falling eGFR trajectory as best as possible.
Furthermore, this study highlights the significance of estimating ΔeGFR for predicting future outcomes. Indeed, in addition to the last eGFR level, the past eGFR trend has been shown to predict a patient’s future risk of ESRD, especially those with advanced CKD in whom the absolute risk of ESRD is higher . This raises the question as to whether ΔeGFR can be incorporated into current risk prediction tools such as the well-validated Kidney Failure Risk Equation (KFRE), which predicts the 2- and 5-year risk of ESRD in patients with CKD stage 3–5 . However, given that it is equally relevant to reconcile whether death may be more likely than ESRD, the CKD Prognosis Consortium risk calculator provides the 2- and 4-year risk of ESRD as well as the risk of non-fatal cardiovascular events and death prior to ESRD in patients with CKD stage 4, and provides an estimation of the timing of these events in relation to ESRD . Again, whether quantification of the ΔeGFR could improve the risk score in this prediction tool is an area for future research. Further work is also required to determine what time period the ΔeGFR should be assessed over. For instance, does inclusion of a patient’s ΔeGFR calculated over the last year, last 2-years or over a longer time period improve predictive utility over current risk prediction models?
Strengths and limitations
This study extends our understanding of CKD progression by providing a closer examination of linear and non-linear patterns of rapid progression on future clinical endpoints. It is strengthened by a systematic approach to patient selection that characterises the CKD pattern robustly.
There are also limitations to our work. By its nature of being an observational study, we could not confirm causal association or fully account for unmeasured confounders. In addition, it is a single-centre study with a predominantly Caucasian population and thus the findings may not be generalisable to other patient populations.
Patients who progress rapidly are exposed to the dual threat of ESRD or mortality prior to ESRD. However, differences in patient characteristics exist between those who progress rapidly in a linear or non-linear fashion, and this has a significant bearing on observed future outcomes. Thus, the phenotypic risk factor profile and individual eGFR trajectory should provide the substrate for personalised therapeutic interventions in this highly vulnerable group of patients in order to deliver optimal CKD care.
Availability of data and materials
The data analysed during the current study are available from the corresponding author on reasonable request.
Angiotensin converting enzyme inhibitor
Acute kidney injury
Angiotensin receptor blocker
Autosomal dominant polycystic kidney disease
Chronic kidney disease
Delta estimated glomerular filtration rate
Diastolic blood pressure
End-stage renal disease
Estimated glomerular filtration rate
Peripheral vascular disease
Systolic blood pressure
Salford Kidney Study
Urine protein:creatinine ratio
Ali I, Chinnadurai R, Ibrahim ST, Green D, Kalra PA. Predictive factors of rapid linear renal progression and mortality in patients with chronic kidney disease. BMC Nephrol. 2020;21:345.
Rosanksy SJ. Renal function trajectory is more important than chronic kidney disease stage for managing patients with chronic kidney disease. Am J Nephrol. 2012;36:1–10.
Collister D, Ferguson T, Komenda P, Tangri N. The patterns, risk factors, and prediction of progression in chronic kidney disease: a narrative review. SeminNephrol. 2016;36:273–82.
Caravaca-Fontán F, Azevedo L, Luna E, Caravaca F. Patterns of progression of chronic kidney disease at later stages. Clin Kidney J. 2018;11:246–53.
Li L, Astor BC, Lewis J, Hu B, Appel LJ, Lipkowitz MS, et al. Longitudinal progression trajectory of GFR among patients with CKD. Am J Kidney Dis. 2012;59:504–12.
Sumida K, Kovesdy CP. Disease trajectories before ESRD: implications for clinical management. Semin Nephrol. 2017;37:132–43.
Levey AS, de Jong PE, Coresh J, El Nahas M, Astor BC, Matsushita K, et al. The definition, classification, and prognosis of chronic kidney disease: a KDIGO controversies conference report. Kidney Int. 2011;80:17–28.
Agarwal R, Duffin K, Laska DA, Voelker JR, Brever MD, Mitchell PG. A prospective study of multiple protein biomarkers to predict progression in diabetic chronic kidney disease. Nephrol Dial Transplant. 2014;29:2293–302.
Fine J, Gray R. A proportional hazards model for the subdistribution of a competing risk. J Am Stat Assoc. 1999;94:496–509.
Gray RJ. A class of K-sample tests for comparing the cumulative incidence of a competing risk. Ann Stat. 1988;16:1141–54.
Weldegiorgis M, de Zeeuw D, Li L, Parving HH, Hou FF, Remuzzi G, et al. Longitudinal estimated GFR trajectories in patients with and without type 2 diabetes and nephropathy. Am J Kidney Dis. 2018;71:91–101.
Levin A, Djurdjev O, Beaulieu M, Er L. Variability and risk factors for kidney disease progression and death following attainment of stage 4 CKD in a referred cohort. Am J Kidney Dis. 2008;52:661–71.
Eriksen BO, Ingrebetsen OC. The progression of chronic kidney disease: a 10-year population-based study of the effects of gender and age. Kidney Int. 2006;69:375–82.
Rebholz CM, Coresh J, Ballew SH, McMahon B, Whelton SP, Selvin E, et al. Kidney failure and ESRD in the atherosclerosis risk in communities (ARIC) study: comparing ascertainment of treated and untreated kidney failure in a cohort study. Am J Kidney Dis. 2015;66:231–9.
Tangri N, Grams ME, Levey AS, Coresh J, Appel LJ, Astor BC, et al. Multinational assessment of accuracy of equations for predicting risk of kidney failure. A meta-analysis. JAMA. 2016;315:164–74.
Potier L, Roussel R, Zeller M, Schiele F, Puymirat E, Simon T, et al. Chronic kidney disease, diabetes, and risk of mortality after acute myocardial infarction: insight from the FAST-MI program. Diabetes Care. 2020;43:e43–4.
Shilpak MG, Heidenreich PA, Noguchi H, Chertow GM, Browner WS, McCellan MB. Association of renal insufficiency with treatment and outcomes after myocardial infarction in elderly patients. Ann Intern Med. 2002;137:555–62.
Margolis G, Vig S, Flint N, Khoury S, Barkagan M, Keren G, et al. Prognostic implications of chronic kidney disease on patients presenting with ST-segment elevation myocardial infarction with versus without stent thrombosis. Cardiorenal Med. 2017;7:150–7.
Toft G, Heide-Jørgensen U, van Haalen H, Glen J, Hedman K, Birn H, et al. Anemia and clinical outcomes in patients with non-dialysis dependent or dialysis dependent severe chronic kidney disease: a Danish population-based study. J Nephrol. 2020;33:147–56.
Hunsicker LG, Adler S, Caggiula A, England BK, Greene T, Kusek JW, et al. Predictors of the progression of renal disease in the modification of diet in renal disease study. Kidney Int. 1997;51:1908–19.
Rangaswami J, Bhalla V, Blair JEA, Chang TI, Costa S, Letine KL, et al. Cardiorenal syndrome: classification, pathophysiology, diagnosis, and treatment strategies: a scientific statement from the American Heart Association. Circulation. 2019;139:e840–78.
Yeh HC, Ting IW, Huang HC, Chiang H-Y, Kuo C-C. Acute kidney injury in the outpatient setting associates with risk of end-stage renal disease and death in patients with CKD. Sci Rep. 2019;9:17658.
Kovesdy CP, Coresh J, Ballew SH, Woodward M, Levin A, Naimark DM, et al. Past decline versus current eGFR and subsequent ESRD risk. J Am Soc Nephrol. 2016;27:2447–55.
Grams ME, Sang Y, Ballew SH, Carrero JJ, Djurdjev O, Heerspink HJL, et al. Predicting timing of clinical outcomes in patients with chronic kidney disease and severely decreased glomerular filtration rate. Kidney Int. 2018;93:1442–51.
Ethics approval and consent to participate
The Salford Kidney Study was granted ethical approval by the North West Greater Manchester South Research Ethics Committee (REC15/NW/0818). Written informed consent was obtained from all patients. The methods described herein were carried out in accordance with relevant guidelines and regulations of the SKS.
Consent for publication
The authors declare no competing interests.
Univariate analysis using Fine-Gray hazards model to investigate factors associated with ESRD. Table S2. Univariate analysis using Fine-Gray hazard model to investigate factors associated with mortality prior to ESRD. Table S3. Baseline characteristics of patients with autosomal dominant polycystic kidney disease, diabetic nephropathy and glomerulonephritis.
About this article
Cite this article
Ali, I., Chinnadurai, R., Ibrahim, S.T. et al. Adverse outcomes associated with rapid linear and non-linear patterns of chronic kidney disease progression. BMC Nephrol 22, 82 (2021). https://doi.org/10.1186/s12882-021-02282-5
- Chronic kidney disease
- End-stage renal disease