Abstract
A decline in renal function suggests progression of chronic kidney disease. This can be determined by measured GFR (e.g., iothalamate clearance), serum creatinine (SCr)-based GFR estimates, or creatinine clearance. A cohort of 234 patients with autosomal dominant polycystic kidney disease and baseline creatinine clearance >70 ml/min were followed annually for four visits. Iothalamate clearance, SCr, and creatinine clearance were obtained at each visit. Estimated GFR (eGFR) was determined with the Modification of Diet in Renal Disease (MDRD) and Cockcroft-Gault equations. Renal function slopes had a mean residual SD of 10.7% by iothalamate clearance, 8.2% by MDRD equation, 7.7% by Cockcroft-Gault equation, and 14.8% by creatinine clearance. By each method, a decline in renal function (lowest quintile slope) was compared among baseline predictors. Hypertension was associated with a decline in iothalamate clearance (odds ratio [OR] 5.8; 95% confidence interval [CI] 2.3 to 14), eGFR (OR [MDRD] 2.0 [95% CI 1.0 to 4.2] or OR [Cockcroft-Gault] 1.9 [95% CI 0.9 to 3.9]), and creatinine clearance (OR 2.0; 95% CI 1.0 to 4.2). Each doubling of kidney volume at baseline was associated with a decline in iothalamate clearance (OR 2.4; 95% CI 1.5 to 3.7), eGFR (OR 1.7 [95% CI 1.1 to 2.6] or 2.1 [95% CI 1.4 to 3.3]), and creatinine clearance (OR 1.7; 95% CI 1.1 to 2.5). Predictor associations were strongest with measured GFR. Misclassification from changes in non-GFR factors (e.g., creatinine production, tubular secretion) conservatively biased associations with eGFR. Misclassification from method imprecision attenuated associations with creatinine clearance.
To treat patients with early chronic kidney disease (CKD), it is important to identify predictors for a decline in renal function. Both slope and threshold analyses of renal function have been used to infer progression of CKD. There are advantages to using a slope analysis for statistical power and for understanding the mechanisms of a disease (1). In a slope analysis, renal function is measured at multiple time points. For each individual, a straight line or curve can be fit to these measurements to reflect the change in renal function over time. Predictors can be compared between individuals who have a decline in renal function (negative slope) with those who have stable or a rise in renal function (neutral or positive slope).
Measured GFR (e.g., inulin or iothalamate clearance), estimated GFR (eGFR) based on serum creatinine (e.g., Modification of Diet in Renal Disease [MDRD] or Cockcroft-Gault equation) (2–4), and creatinine clearance are commonly used as markers of renal function. These three methods have been used to infer changes in renal function over time (5–10). Reciprocal serum creatinine (1/SCr) has also been used to infer changes in renal function over time (7,11,12).
Measured GFR with an exogenous marker such as inulin or iothalamate is the gold standard method, but expense and inconvenience are important drawbacks. Urine collection errors and tubular secretion of creatinine impair the precision and the accuracy of creatinine clearance as a marker of GFR (13–15). Estimating equations are inexpensive and convenient but have not been accurate in populations that are characterized by normal or near-normal renal function (16–18).
One study concluded that measured GFR should be the preferable method for assessing renal function in clinical trials of patients with renal grafts (19). Another study found that eGFR was similar to measured GFR for assessing renal function change over time in a CKD clinical trial (6). More studies are still needed to determine the accuracy of estimating equations for following progression of CKD (20). In particular, the association of demographic or clinical predictors with SCr, creatinine clearance, and measured GFR slope need further comparison (1). The objective of this study was to compare methods for a decline in renal function for a cohort with early autosomal dominant polycystic kidney disease (ADPKD).
Materials and Methods
Study Participants
Details of the Consortium for Radiologic Imaging Studies of Polycystic Kidney Disease (CRISP) have previously been published. CRISP is a multicenter study of the natural history of ADPKD (21) in adherence to the Declaration of Helsinki. The sites included Kansas University Medical Center, Emory University, the Mayo Clinic, and University of Alabama at Birmingham. A cohort of patients who had ADPKD and were aged 15 to 46 yr were followed annually for four total visits. Eligibility required an ADPKD diagnosis by the criteria of Ravine et al. (22) and a measured or estimated (Cockcroft-Gault equation) creatinine clearance >70 ml/min. Eligibility also required a SCr ≤1.6 mg/dl in men and ≤1.4 mg/dl in women. Patients were ineligible when they had other medical conditions besides hypertension that could affect renal function (e.g., diabetes).
Baseline Predictors
During the baseline visit in 2001, patients underwent a 2-d evaluation. A focused history and examination was used to identify relevant characteristics and medical conditions. In particular, current smoking status, history of urinary tract infection (UTI), abdominal pain, and history of gross hematuria were determined. Hypertension was defined by current use of antihypertensive medications or a systolic and diastolic BP ≥140/90 mmHg on multiple occasions. Magnetic resonance imaging was used to determine the cyst volume and total volume of the combined kidneys (21). The day before the baseline visit, a 24-h urine sample was collected to determine the urine albumin to creatinine ratio (ACR).
Renal Function Measurements
Renal function was measured by each method at baseline (2001) and at the follow-up visits in 2002, 2003, and 2004. Details of the measurement of iothalamate clearance have been described previously (21,23). Briefly, after oral hydration, patients received a subcutaneous injection of nonradiolabeled iothalamate. After a 60-min equilibrium period, each patient voided and the first plasma sample was drawn. After a timed 45- to 60-min collection period to determine urine flow (V), a voided urine sample and a second plasma sample were obtained. Postvoid residuals were assessed by ultrasound after each void. The two plasma (P) samples and one urine (U) sample were assayed for iothalamate via capillary electrophoresis at the Mayo Clinic. Iothalamate concentrations in the plasma samples were averaged, and GFR was determined using the clearance equation (UIothalamateV/PIothalamate).
During the feasibility phase, four patients had an iothalamate clearance test at each of the four study sites within a 6-wk period. The mean between-site coefficient of variation was 9.8% but decreased to 4.9% with exclusion of three suboptimal iothalamate clearance tests. A suboptimal test was defined a priori: (1) A postvoid residual of at least 20 ml and >20% voided volume; (2) urine flow rate <3 ml/min; or (3) a collection time <30 or >90 min. A suboptimal test occurred at least once in 22% of the study patients. Exclusion of or adjustment for suboptimal tests had no substantive effect on the predictor associations and was not done in the final analyses.
Serum was collected at each visit and assayed for creatinine at each site. All SCr levels were adjusted for calibration bias with the Cleveland Clinic laboratory used for deriving the MDRD equation. For determination of this calibration bias, 90 serum samples were exchanged between the Mayo Clinic laboratory and the MDRD Study laboratory. An additional 31 serum samples were exchanged among all four study sites to determine calibration bias between sites. The final calibration equations were as follows: SCr = 0.01 + 1.30 × SCrKansas, SCr = −0.01 + 1.08 × SCrEmory, SCr = −0.33 + 1.28 × SCrMayo, and SCr = 0.04 + 0.97 × SCrAlabama. GFR was estimated using the MDRD equation 186.3 × SCr−1.154 × age−0.203 × 0.742 (if female) × 1.212 (if black) (3,4) and the Cockcroft-Gault equation [(140 − age) × weight × 0.85 (if female)]/(72 × SCr) (2). The age and the weight at each visit were used to calculate eGFR. Creatinine clearance was determined with the 24-h urine creatinine at each visit.
Statistical Analyses
Analyses were performed with JMP 5.1 (SAS Institute, Cary, NC). The annual percentage change in renal function (i.e., slope) was determined by regressing logarithmic renal function on time from baseline (in years) for each patient. Whereas higher order curves may be needed for more advanced CKD with longer follow-up (6), a line was adequate for describing changes over time in this study (Figure 1). For iothalamate clearance, the annual percentage change was determined with and without standardizing body surface area (per 1.73 m2) with a formula that used the height and the weight at each visit (24). For SCr-based methods, slopes were determined for the MDRD equation, Cockcroft-Gault equation, and 1/SCr. The creatinine clearance slopes were not standardized (results were similar with standardizing body surface area). The correlation of each renal function slope with unstandardized iothalamate clearance slope was compared.
Percentage change in renal function from baseline at each follow-up visit by unstandardized iothalamate clearance (A), Modification of Diet in Renal Disease (MDRD) equation (B), Cockcroft-Gault equation (C), and creatinine clearance (D). The line represents the pooled data regression slope and shows relatively stable renal function for the overall cohort by each method. Box plots for each follow-up visit show more variability in the change in renal function by iothalamate clearance and creatinine clearance than by the MDRD or Cockcroft-Gault equation. This was consistent with greater measurement error by the clearance methods than by the equation methods.
For each patient, some of the variability in renal function between follow-up visits was due to physiologic variability (e.g., changes from diet and hydration) (25) and assay error. This measurement error or imprecision represents variability that was not accounted for by the overall change in renal function with time (i.e., slope) and was reported as a residual SD (26). The mean within-patient residual SD was compared between renal function methods with the paired Wilcoxon sign-rank test.
Because of imprecision in renal function measurements, a decline in renal function defined simply by a negative slope would lead to substantial misclassification. For obviating this error, a decline in renal function was defined as the lowest quintile slope (n = 47) for each method. The agreement between a decline in renal function slope by unstandardized iothalamate clearance compared with each other method was assessed with a κ statistic (0 = agreement no greater than chance, 1 = complete agreement). Logistic regression was used to compare baseline predictors with a declining versus neutral or increasing renal function slope. The statistical significance of a decline in renal function with respect to baseline predictors was determined by the Wald test. All odds ratios (OR) were adjusted for study site. Continuous baseline predictors were evaluated per doubling for kidney volume, cyst volume, and urine ACR and per 10 yr for age. Backward elimination model building was used to identify independent predictors (P < 0.05) for a decline in renal function by each method. Additional analyses compared renal function slope methods using the t test for nominal predictors and Pearson correlation for continuous predictors.
The MDRD equation was developed with cross-sectional data to estimate standardized iothalamate clearance in patients with CKD and decreased renal function (3). With the MDRD equation, there is a static relationship between equation variables (ln SCr, ln age, gender, and race) and ln measured GFR (standardized iothalamate clearance) based on a linear regression model. Using this same multivariable regression model, the relationship between these variables and standardized iothalamate clearance was compared at each of the four visits. An additional model included a quadratic term for logarithmic serum creatinine [ln(SCr)]2 to assess curvature (18). Using a smoother function (JMP 5.01, λ = 0.1), the relationship between iothalamate clearance and SCr was graphically compared between the combined first two study visits and the combined last two study visits.
Results
The baseline characteristics of the CRISP cohort are shown in Table 1. Of 241 patients who enrolled in the study, seven were excluded for having only one study visit. By each renal function method, slopes could be calculated in the remaining 234 patients. By iothalamate clearance, this involved measures from all four visits for 82%, three visits for 12%, and two visits for 6% (missing data trends similar for creatinine clearance and SCr). Figure 1 shows the annual percentage change in renal function from baseline for unstandardized iothalamate clearance, the MDRD equation, the Cockcroft-Gault equation, and creatinine clearance. Individual slopes by each method followed a normal distribution. There was a gradual but statistically significant decline in renal function by unstandardized iothalamate clearance (mean slope = −1.5%/yr; P = 0.009; SD = 9.3%), by MDRD equation (mean = −3.0%/yr; P < 0.001; SD = 7.0%), and by Cockcroft-Gault equation (mean = −1.8%/yr; P < 0.001; SD = 6.7%). This decline was not statistically significant by creatinine clearance (mean = −0.6%/yr; P = 0.48; SD = 13.3%). As shown by the box plots and the SD of slopes, there was more slope variability by iothalamate clearance and creatinine clearance than by the estimating equations.
Clinical characteristics at baseline for 234 patients with ADPKDa
Each renal function method was compared with unstandardized iothalamate clearance. The Pearson correlation with unstandardized iothalamate clearance slope was r = 0.98 (P < 0.001) for standardized iothalamate clearance slope, r = 0.30 (P < 0.001) for 1/SCr slope, r = 0.30 (P < 0.001) for MDRD equation slope, r = 0.25 (P < 0.001) for Cockcroft-Gault equation slope, and r = −0.06 (P = 0.34) for creatinine clearance slope. The agreement with unstandardized iothalamate clearance for lowest quintile slope was κ ± SE = 0.92 ± 0.03 for standardized iothalamate clearance, 0.23 ± 0.07 for 1/SCr, 0.25 ± 0.07 for MDRD equation, 0.15 ± 0.07 for Cockcroft-Gault equation, and 0.20 ± 0.07 for creatinine clearance.
The precision between renal function methods was compared. The mean residual SD was 10.7% for both unstandardized and standardized iothalamate clearance. For unstandardized iothalamate clearance, patients with at least one suboptimal test had a mean residual SD of 13.1%, whereas patients with no suboptimal tests had a mean residual SD of 10.1%. The mean residual SD was 6.9% for 1/SCr, 8.2% for the MDRD equation, and 7.7% for the Cockcroft-Gault equation, and these all were lower (i.e., more precise) than unstandardized or standardized iothalamate clearance (P < 0.001 for all). The mean residual SD was 14.8% for creatinine clearance, and this was higher (i.e., less precise) than the unstandardized or standardized iothalamate clearance (P < 0.001 for both). Across each study site (Kansas, Emory, Mayo Clinic, and Alabama), the mean residual SD was 10.5, 11.0, 10.0, and 13.8% for unstandardized iothalamate clearance; 6.6, 7.2, 6.7, and 7.7% for 1/SCr; 7.7, 8.4,7.8, and 10.1% for the MDRD equation; and 18.3, 10.7, 12.3, and 22.1% for creatinine clearance.
Baseline predictors were compared between patients with and without a decline in renal function as defined by the lowest quintile slope (Table 2). Statistically significant predictors of a decline in both unstandardized and standardized iothalamate clearance included hypertension, kidney volume, cyst volume, urine ACR, and history of UTI (all P < 0.05). Age was borderline significant (P = 0.06). For a decline in 1/SCr, eGFR (MDRD equation or Cockcroft-Gault equation), or creatinine clearance, most of these same predictor associations were present but were attenuated. For estimating equations, constant multipliers (gender and race variables) had no effect on the logarithmic slopes. In general, predictor associations were similar between a decline in 1/SCr and a decline in eGFR. There were similar differences between methods for associations with renal function slopes as continuous variables (Table 3).
Association between baseline predictors and lowest quintile (n = 47) renal function slope by different methods (95% CI)a
Association between baseline predictors and renal function slope by different methodsa
In a multivariable model, independent predictors for a decline in unstandardized iothalamate clearance were hypertension with OR of 3.5 (95% confidence interval [CI] 1.3 to 9.3), kidney volume per doubling with OR of 1.9 (95% CI 1.2 to 3.1), and history of UTI with OR of 2.3 (95% CI 1.1 to 4.6). For this sample size, multiple independent predictors for renal function decline were identified only by the iothalamate clearance. By eGFR, 1/SCr, and creatinine clearance, only one predictor remained statistically significant with backward stepwise elimination.
For further exploration of the discrepancy between measured GFR slope compared with eGFR slope, logarithmic standardized iothalamate clearance was regressed on the same variables used in the MDRD equation (ln SCr, ln age, gender, and race) for each of the four study visits (Table 4). For the SCr, female gender, and black race variables, the strength of association with standardized iothalamate clearance increased progressively from baseline through the follow-up visits. The association between these variables and standardized iothalamate clearance also became more consistent with the MDRD equation with each follow-up visit. Variability in ln GFR also increased progressively from an SD of 0.25 at baseline to an SD of 0.35 at the final visit(3 yr later). The relationship between ln SCr and ln GFR was relatively linear for the combined first two rounds but developed more curvature for the combined last two rounds(Figure 2). This curvature could be modeled by addition of a quadratic term [ln(SCr)2] to the regression model shown in Table 4. Inclusion of this quadratic term increased the model fit(R2) by 0.028(P < 0.001) at baseline, 0.005(P = 0.10) at follow-up visit 1, 0.022(P < 0.001) at follow-up visit 2, and 0.031(P < 0.001) at follow-up visit 3.
A log-log plot of serum creatinine versus iothalamate clearance for the first two study visits(○) compared with the last two study visits(•). Smoother curves(JMP 5.01, λ = 0.1) regressed standardized iothalamate clearance on SCr for the first two study visits(gray curve) compared with the last two study visits(black curve). There was more downward curvature in the last two study visits compared with the first two study visits.
Association between predictor variables and measured GFR at each annual visit with the same multivariable regression model used for the MDRD equationa
Discussion
In a cohort of patients who had ADPKD and baseline normal renal function, predictors for a decline in measured GFR (iothalamate clearance) were hypertension, increased kidney volume, increased cyst volume, increased urine ACR, and history of UTI(Table 2). With creatinine clearance, these associations were attenuated from more measurement error(mean residual SD 14.8 versus 10.7%; P < 0.001). With 1/SCr or eGFR, there was less measurement error(mean residual SD ≤ 8.2 versus 10.7%; P < 0.001), but predictor associations were weaker instead of stronger compared with measured GFR. Therefore, changes in non-GFR factors conservatively biased associations with the SCr-based methods. This supports the hypothesis that the rise in SCr with a decline in GFR is less than expected asa result of concurrent muscle atrophy, decreased dietary protein, and increased tubular creatinine secretion(17). To a different extent between methods, a decline in renal function can be misclassified from measurement error(precision) or from measurement bias(Table 5).
Sources of misclassification for a decline in renal function
Many of the predictors that were associated with a decline in renal function have previously been reported for ADPKD. Proteinuria and increased cyst volume have previously been identified as prognostic factors for a decline in measured GFR(27). Kidney volume, proteinuria, and age have been associated with a decline in eGFR(9). In this study, history of UTI was also identified as a prognostic factor for a decline in measured GFR. Whether this was related to renal parenchymal damage from UTI or misdiagnosis of hematuria, leukocyturia, flank pain, or other renal cyst symptoms as a UTI is unclear. In this study, hypertension, kidney volume, and history of UTI were independent prognostic factors for a decline in measured GFR.
In cohorts with baseline normal renal function, it may be preferable to use 1/SCr slope rather than eGFR slope because the relationship between equation variables and measured GFR can be dynamic(Table 4). Furthermore, the mathematical model of the equation can be dynamic with changes in curvature during the follow-up(Figure 2). When a baseline cohort is characterized by normal or near-normal renal function, most of the variability in SCr may reflect creatinine production(e.g., muscle mass). Gender may have a weak effect at baseline(Table 4), as a result of collinearity with SCr(both are markers of muscle mass [28]). However, as disease progression increased GFR variability in the cohort during the follow-up period (Figure 2), SCr became more a marker of GFR and less a marker of muscle mass. This results in gender having less collinearity with SCr. Subsequently, gender more strongly adjusts for confounding between SCr and GFR (Table 4). This change in the relationship among SCr, gender, and GFR was consistent with previous studies that compared populations that had a normal and narrow range of renal function with populations that had a reduced and wide range of renal function (17,18).
A decline in creatinine clearance had weaker associations with predictors than a decline in iothalamate clearance. The higher within-patient residual SD for creatinine clearance than for iothalamate clearance could explain much of this discrepancy. Incomplete or overcollection of 24-h urine samples can affect precision and accuracy. Some of the change in creatinine clearance over time may reflect changes in creatinine secretion or changes in patient effort for collection of urine samples. Adjustment for study site increased the strength of predictor associations with a decline in creatinine clearance but had no substantive effect on predictor associations by the other methods (unadjusted data not shown). This can be explained by the greater discrepancy between sites with the precision of creatinine clearance than with other methods. Another study that compared longitudinal and simultaneously collected creatinine clearance with measured GFR found substantial discordance between these two methods over time (15).
When slopes or changes in renal function were compared, there were no advantages to standardizing iothalamate clearance for body surface area (Tables 2 and 3). In cross-sectional or threshold analyses (e.g., GFR <60 ml/min per 1.73 m2), standardizing for body surface area attempts to account for differences in physiologic demand for renal function (i.e., larger individuals need more absolute renal function than smaller individuals). In this longitudinal study, however, percentage changes in renal function were compared and each patient served as his or her own control for physiologic demand. Adjustment for body surface area provides an additional source of measurement error with the added height and weight variables. Also, some longitudinal changes in height (e.g., compression fractures) or weight (e.g., obesity or wasting) will be interpreted incorrectly as changes in renal function. This can also explain the discrepancy with predictor associations between the Cockcroft-Gault equation (weight variable present) and the MDRD equation (weight variable absent; Tables 2 and 3).
In this study, logarithmic slopes (% per year) were used instead of absolute slopes (ml/min per year or ml/min per 1.73 m2/yr) because measurement error decreases as renal function decreases (3). For example, a rise in SCr from 0.9 to 1.0 mg/dl or from 3.3 to 6.0 mg/dl represents a decline in eGFR of 10 ml/min per 1.73 m2 for a 60-yr-old white man with CKD (MDRD equation). The former can be attributed easily to measurement error, whereas the latter more likely represents a true decline in GFR. With logarithmic slopes, a rise in SCr from 0.9 to 1.6 mg/dl would be similar to a rise in SCr from 3.3 to 6.0 mg/dl. In this study, associations with baseline predictors were generally more attenuated with analyses by absolute slope. Hypertension was associated with a decline (lowest quintile slope) in unstandardized iothalamate clearance by an OR of 4.0 (95% CI 1.8 to 9.7), MDRD equation by an OR of 1.1 (95% CI 0.6 to 2.2), and creatinine clearance by an OR of 0.9 (95% CI 0.7 to 1.9). Each doubling of kidney volume was associated with a decline in unstandardized iothalamate clearance by an OR of 2.4 (95% CI 1.6 to 3.7), MDRD equation by an OR of 1.1 (95% CI 0.8 to 1.7), and creatinine clearance by an OR of 1.2 (95% CI 0.7 to 1.9).
The stronger predictor associations by measured GFR slope may not be applicable to CKD populations with a more reduced baseline renal function. Lewis et al. (6) found predictor associations to be similar by measured GFR slope and eGFR slope. Patients in that study had a baseline measured GFR of 20 to 65 ml/min per 1.73 m2, and there was a more substantial overall decline in renal function during follow-up. Different methods for assessing the rate of renal function decline may correlate better in advanced CKD compared with early CKD. However, when the goal is to intervene early in CKD, measured GFR may be superior for identifying the initial decline in renal function.
A potential limitation of this study was protocol differences between renal function methods. Iothalamate concentrations were assayed at a centralized laboratory, whereas SCr and creatinine clearance were assayed at each study site. This was addressed by adjusting for study site in the analyses. There was also good agreement between lowest quintile slope for the MDRD equation with SCr calibrated versus MDRD equation without SCr calibrated (κ ± SE = 0.95 ± 0.03). Furthermore, predictor associations did not change substantively with SCr calibration. Patients in this study all had ADPKD and were younger than the typical patient with CKD. Consequently, generalization to other CKD causes and older populations should be done with caution. However, ADPKD may also be useful as a model for CKD, because the diagnosis can be made with radiologic imaging before a detectable decline in renal function. Finally, only a few predictors for renal function decline were identified in this study, and further studies are still needed to compare renal function methods by other baseline predictors.
Conclusion
Measured GFR decline had the strongest association with predictors. Measurement error with creatinine clearance attenuated associations because of increased misclassification for a decline in renal function. SCr-based methods (1/SCr or eGFR) had less measurement error compared with measured GFR. However, changes in SCr from non-GFR factors (e.g., creatinine production, tubular secretion) led to conservatively biased associations. The cost and the convenience of different methods still need consideration when assessing changes in renal function over time.
Acknowledgments
This study was supported by F32-DK68996, by the CRISP study UO1 cooperative agreements (DK56956, DK56943, DK56957, and DK56961) of the National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, and the General Clinical Research Center grants MO1-RR00585 (the Mayo Foundation), MO1-RR00052 (University of Alabama at Birmingham), and M01-RR00039 (Emory University) for the Division of Research Resources, National Institutes of Health.
Footnotes
Published online ahead of print. Publication date available at www.jasn.org.
- © 2006 American Society of Nephrology