Abstract
Estimation of the GFR (eGFR) using creatinine- or cystatin C–based equations is imperfect, especially when the true GFR is normal or near-normal. Modest reductions in eGFR from the normal range variably predict cardiovascular morbidity. If eGFR associates not only with measured GFR (mGFR) but also with cardiovascular risk factors, the effects of these non–GFR-related factors might bias the association between eGFR and outcome. To investigate these potential non–GFR-related associations between eGFR and cardiovascular risk factors, we measured GFR by iohexol clearance in a sample from the general population (age 50 to 62 years) without known cardiovascular disease, diabetes, or kidney disease. Even after adjustment for mGFR, eGFR associated with traditional cardiovascular risk factors in multiple regression analyses. More risk factors influenced cystatin C–based eGFR than creatinine-based eGFR, adjusted for mGFR, and some of the risk factors exhibited nonlinear effects in generalized additive models (P < 0.05). These results suggest that eGFR, calculated using standard creatinine- or cystatin C–based equations, partially depends on factors other than the true GFR. Thus, estimates of cardiovascular risk associated with small changes in eGFR must be interpreted with caution.
Chronic kidney disease (CKD) has been established as a risk factor for cardiovascular disease (CVD).1 There are reasons to believe that increased risk may be present even at GFRs presently considered as normal or near normal. An association between normal or near normal GFR and CVD or cardiovascular mortality has been found in some studies2–4 but not in others.5–7 The most important limitation of these studies has been the use of creatinine- or cystatin C–based formulas to estimate GFR. Because precise methods for measuring GFR are difficult and costly to implement in large cohorts, estimated GFR (eGFR) has been accepted as a proxy for measured GFR (mGFR).
Different methods for calculating eGFR have given different results. Estimates of GFR based on cystatin C have been found to predict CVD mortality better than creatinine-based eGFR. This has been shown in the elderly, in patients with known CKD or CVD, and in the general population.2,8–13 The reason for the different performance of the two estimates of GFR may be either different abilities to estimate GFR or a dependence on non-GFR confounders.14 More specifically, the effects of non-GFR factors on cystatin C and creatinine might bias the association between eGFR and CVD. The occurrence and strength of such associations have thus far not been studied by an accurate method of GFR measurement in the general population.
In the Renal Iohexol Clearance Survey in Tromsø 6 (RENIS-T6), GFR was measured by iohexol clearance and estimated by creatinine or cystatin C in a middle-aged cohort from the general population. The aim of this cross-sectional study was to explore the validity of using eGFR as a proxy for mGFR in studies of CVD risk. We assessed the associations between traditional cardiovascular risk factors and eGFR when mGFR, measured as iohexol clearance, was adjusted for. We also wanted to study the relationship between mGFR and cardiovascular risk factors.
RESULTS
In the sixth Tromsø population survey (Tromsø 6), 739 of the 3564 participants 50 to 62 years of age reported a previous myocardial infarction, angina pectoris, stroke, diabetes mellitus, or kidney disease. Of the remaining 2825 persons eligible for inclusion, 1632 were investigated in our study. Five were excluded because their iohexol clearance measurements suffered technical failures. This left us with 1627 persons for this analysis (826 women and 801 men). The characteristics of the participants are described in Table 1. Mean mGFR was 87.8 ml/min per 1.73 m2 in women and 95.7 ml/min per 1.73 m2 in men.
Mean (SD) of study population characteristics
The distributions of cardiovascular risk factors, adjusted for age and medication according to gender and quintiles of mGFR, are shown in Table 2. Multiple linear regression analyses with mGFR as the dependent variable and the cardiovascular risk factors as the independent variables are shown in Table 3. In the fully adjusted model, there were statistically significant associations between mGFR and smoking and diastolic BP. We tested for gender interaction with all variables in the fully adjusted model and found interactions with body mass index (BMI) and HDL cholesterol (P < 0.05). For men, mGFR increased 1.93 ml/min per 1.73 m2 per SD of BMI (P = 0.002) and 1.79 ml/min per 1.73 m2 per SD of HDL cholesterol (P = 0.003). For women, these estimates were not statistically significant. When the same variables as in model 2 (Table 3) were analyzed in a generalized additive model, nonlinear effects of diastolic BP and triglycerides were found (P < 0.05). Diastolic BP greater than about 90 mmHg was associated with reduced mGFR, whereas triglycerides greater than about 2 were associated with increased mGFR (data not shown).
Risk factors for cardiovascular disease by gender-specific quintiles of measured glomerular filtration rate
Multiple linear regression analysis of measured GFR and cardiovascular risk factors
Multiple linear regression analyses with creatinine- and cystatin C–based eGFR as the dependent and cardiovascular risk factors as the independent variables are shown in Table 4. Models 2 and 3 included adjustment for mGFR to assess the effects of cardiovascular risk factors on eGFR not mediated by variation in GFR. In the fully adjusted model 3, only current smoking influenced the Chronic Kidney Disease Epidemiology Collaboration (CKD-EPI) estimates, whereas the Modification in Renal Disease (MDRD)-estimates were influenced by both current smoking and having a first-degree relative with a myocardial infarction (MI) before the age of 60 years (P < 0.05). Current smoking, BMI, HDL cholesterol, and triglycerides were associated with the cystatin C–based estimates (P < 0.05). No statistically significant interactions between gender and any of the other independent variables were found, except for a small effect of the interaction with LDL cholesterol for the MDRD equation (P = 0.048). Analyses of the same variables as in model 3 with generalized additive models found nonlinear effects of BMI on the MDRD and CKD-EPI estimates and of triglycerides on the cystatin C–based estimates (P < 0.05; Figure 1). Errors-in-variables regression with same variables as in model 3 found the same pattern of statistically significant effects when the reliability of mGFR was set at 0.7, except for a small effect of diastolic BP on the cystatin C–based estimate (P = 0.04). The magnitude of the effects was also similar to those in model 3 (Table 4; data not shown).
Multiple linear regression analyses of association between estimates of GFRs and cardiovascular risk factors
The nonlinear relationships between estimated GFR and cardiovascular risk factors after adjustment for measured GFR. The gray area indicates the 95% confidence interval.
Table 5 shows the 10-year risk of coronary heart disease estimated with the Framingham equation according to gender- and method-specific quintiles. The risk in women fell with decreasing mGFR and with both the creatinine-based eGFRs (P < 0.05). In contrast, in both genders, the CVD risk increased substantially with decreasing cystatin C–based eGFR (P < 0.0001). When the risk was adjusted for smoking, the differences in the risk estimates across the quintiles for mGFR and creatinine-based eGFR were not statistically significant, whereas the differences for cystatin C–based eGFR were attenuated but still significant (P for linear trend 0.001 for women, <0.0001 for men; data not shown).
Ten-year cardiovascular riska (%) estimated with the Framingham equation by gender- and method-specific quintiles of measured or estimated GFR
DISCUSSION
We found that eGFR based on either creatinine or cystatin C was influenced by traditional cardiovascular risk factors even after adjusting for mGFR, indicating that the mechanisms for these associations were independent of GFR. This was especially true for eGFR based on cystatin C, which was associated with several cardiovascular risk factors in a fully adjusted model (Table 4, model 3). The creatinine-based MDRD and CKD-EPI equations seemed to be less influenced, but both exhibited a nonlinear relationship with BMI, which resulted in an increased eGFR for persons with a BMI exceeding about 30 kg/m2 (Figure 1, A and B). This could bias the estimated cardiovascular risk from reduced GFR in obese persons in longitudinal studies. There was also a nonlinear effect of triglycerides on eGFR from cystatin C (Figure 1C).
Associations between cystatin C and cardiovascular risk factors have been found in at least two recent studies.15,16 When adjusting for GFR measured as creatinine clearance, Knight et al.15 found that cystatin C was associated with age, gender, weight, height, and current smoking but not with hyperlipidemia. The limited precision of creatinine clearance for measuring GFR and the inclusion of persons with prevalent cardiovascular disease may account for the differences from our study. Stevens et al.16 found associations between cystatin C and BMI and systolic BP after adjusting for GFR measured as iothalamate or Cr-EDTA clearance. In the same study, the authors found associations between creatinine and BMI. HDL cholesterol and triglycerides were not examined in this study. The results are not directly comparable to this study because the study population consisted of CKD patients with a mean GFR of 44 ml/min per 1.73 m2 and because the associations were not analyzed in multiple regression analyzes. Neither Knight et al. nor Stevens et al. tested for nonlinear effects.
When examining the relationships between cardiovascular risk factors and GFR measured as iohexol clearance, complex associations were found in this healthy middle-aged study population (Table 3). The impact of BMI and HDL cholesterol was different in men and women. To our knowledge, RENIS-T6 is the first study where an accurate method for measuring GFR has been used in a large cohort from the general population. All previous studies of the relationship between cardiovascular risk factors and GFR in the general population have relied on estimates based on creatinine or creatinine clearance.2,6,17–19 In some of these studies, the relationships were confounded by the inclusion of persons with established CVD or kidney disease,2,17–19 some did not perform multiple regression analyses,2,6,17 or, in some, the interactions between gender and the risk factors were not tested.2,6,17,19 However, Nitsch et al.18 found interactions between gender and several cardiovascular risk factors in a cross-sectional study from a Swiss population. U-shaped relationships between creatinine-clearance and systolic BP, diastolic BP, serum triglycerides, and C-reactive protein were found in the Prevention of Renal and Vascular Disease study.
We found a positive association between current smoking and increasing mGFR. This confirms findings in several studies using creatinine-based eGFR.2,6,17,18,20,21 In these studies, a possible explanation was that GFR had been overestimated in smokers because they have lower muscle mass and lower creatinine than nonsmokers, but our results show that current smoking is truly associated with an increase in GFR. However, an association between smoking and creatinine-based eGFR persisted even after adjustment for mGFR, which means that reduced muscle mass and creatinine production are partly responsible for the relationship between smoking and eGFR (Table 4, models 2 and 3). The mechanism for smoking-induced hyperfiltration is uncertain because several experimental studies have shown a reduced GFR immediately after smoking or the administration of nicotine.22,23 Because our subjects had been instructed to avoid tobacco before the measurement of iohexol clearance, the higher mGFR in smokers may have been caused by a withdrawal phenomenon.
To explore the magnitude of the combined effect of the traditional risk factors with both mGFR and eGFR, we studied the distribution of the 10-year risk of coronary heart disease estimated with the Framingham equation across method-specific quintiles of GFR. Although both mGFR and creatinine-based eGFR exhibited similar patterns with an increasing risk for women with increasing GFR quintiles, there was a marked gradient in the opposite direction for cystatin C–based eGFR for both genders (Table 5). However, after adjustment for smoking, there was nonsignificant differences in the 10-year risk for coronary heart disease across the quintiles of mGFR and creatinine-based eGFR, but the difference between the cystatin C–based quintiles remained statistically significant. These findings suggest that correct adjustment for confounding cardiovascular risk factors in survival analysis is more critical for GFR estimates based on cystatin C than for creatinine-based estimates.
eGFR has been accepted as a proxy for mGFR in epidemiologic studies because of the lack of GFR measurements. Longitudinal studies of eGFR as a cardiovascular risk factor have used linear adjustments for the traditional risk factors. None of them incorporated nonlinear effects,2–4,6 and gender-specific analyses were only reported in the study of Astor et al. Thus, there is a possibility of residual confounding from cardiovascular risk factors in these survival analyses, but it is difficult to assess the size and direction of the resulting bias. Although the magnitude of most of the nonlinear relationships found in this study was modest, even small effects may result in biased estimates in large studies and should be considered. The dependence of eGFR on cardiovascular risk factors in addition to GFR makes it difficult to interpret risk estimates and might limit the usefulness of eGFR as proxy for GFR in longitudinal studies of persons with near-normal GFR.
The most important strength of this study was its use of GFR measurements rather than estimates in a large cohort from the general population. Plasma iohexol clearance is an accurate method for measuring GFR.24 The single-sample method we used has been found to be at least as reliable as the multi-sample method when both were compared with CrEDTA-clearance, especially for GFR >60 ml/min per 1.73 m2.25 Even so, measurement error and day-to-day variation in GFR could conceivably have introduced bias in our analyses. However, the results of errors-in measurements regression, which takes measurement error in independent variables into account, gave results very similar to those given in Table 4. Because we did not include repeat measurements of GFR in the same persons in this study, the reliability of mGFR used in these analyses was set at values far lower than found in the literature to provide a conservative assessment of possible bias.
One limitation of this study was its cross-sectional design. A definitive assessment of the confounding effects of non-GFR factors in cohort studies of eGFR and CVD can only be made in survival analyses comparing the effects of mGFR and eGFR directly, but data on endpoints is not available at this time. The inclusion of only middle-aged participants was another limitation. Nonlinear associations have been found in studies with wider age ranges, however, and we expect that our findings are valid for other age groups than those examined. Because no non-whites were included, generalizability to other ethnic groups may also be limited.
We conclude that GFR estimated by either creatinine or cystatin C is dependent on non–GFR-related factors in addition to GFR. Some of these associations were nonlinear. The possibility of residual confounding from these factors in studies of GFR and cardiovascular risk in persons with a GFR close to the normal range cannot be ruled out. The estimates of the effect of eGFR on CVD in these studies should be interpreted with caution.
CONCISE METHODS
Study Population and Data Collection
The Tromsø study is a population-based, prospective study with repeated health surveys of inhabitants of the municipality of Tromsø in northern Norway. In the sixth Tromsø study (Tromsø 6), which was performed in 2007/2008, an age-stratified random sample of all inhabitants between 40 and 87 years of age were invited. Of the 65,286 inhabitants in Tromsø county, 12,984 persons met for the examination. The attendance rate for the invited participants was 66%. They answered an extensive questionnaire about smoking habits, physical activity, medication, and sociodemographic characteristics.
In the age group between 50 and 62 years, a total of 3564 participants completed the main part of Tromsø 6. Based on answers from the self-administered questionnaire, participants who reported no previous MI, angina pectoris, stroke, diabetes mellitus, or any renal disease (except urinary tract infection) were invited to the RENIS-T6. They were stratified according to age and gender and included to a predetermined target of 1600 participants. Pregnant women and individuals reporting allergies to radiocontrast, iodine, or latex were excluded.
Measurements
All investigations for this study took place at the clinical research unit at the University Hospital of North Norway. Participants met between 8:00 and 10:00 am. They were instructed to avoid large meals with meat and nonsteroid anti-inflammatory drugs in the 48 hours before examination. Smokers were asked to be abstinent from tobacco for the 8 hours before they met. The participants were fasting, but they were instructed not to restrict fluid intake by drinking two to three glasses of water before their arrival at the clinical research unit. If the person had an acute illness on arrival, the examination was rescheduled.
Fasting plasma samples were obtained for analyses of creatinine, cystatin C, triglycerides, total cholesterol, HDL cholesterol, and LDL cholesterol. These analyses were performed on the same day at our laboratory at the University Hospital of North Norway.
Five milliliters of iohexol (Omnipaque 300 mg I/ml) was administered through a Teflon catheter placed in an antecubital vein. The syringe was weighed before and after administration. The catheter was flushed with 30 ml of isotonic saline. The sample for iohexol concentration measurement was drawn through the same catheter. After injection of the iohexol, the participants were given a light breakfast and allowed to walk about until it was time to obtain the blood sample for measuring the iohexol concentration. The optimal time for sampling was calculated by Jacobsson's method based on GFR estimated from plasma creatinine measured previously in Tromsø 6.26 The exact time in minutes from injection to sampling was measured with a separate stopwatch for each person.
The serum iohexol concentration was measured by HPLC, as described previously by Niellson-Ehle and Grubb.27 The interassay coefficient of variation (CV) for the iohexol analysis with HPLC during the study period was 3.0%. External quality control was provided by Equalis (Equalis AB, Uppsala, Sweden).
The GFR was calculated using the formula described by Jacobsson.26 Extracellular volume was estimated by Graneus' equation.28 Extrarenal iohexol clearance was ignored in accordance with the practice of other authors. GFR was normalized to 1.73 m2 body surface area, which was estimated using Dubois' equation.29
Plasma creatinine analyses were performed with an enzymatic method that has been standardized against isotope dilution mass spectroscopy (CREA Plus; Roche Diagnostics, Mannheim, Germany). Interassay CV for creatinine in the study period was 2.3%. Cystatin C was measured with a particle-enhanced turbidimetric immunoassay using reagents from Gentian (Moss, Norway) on a Modular E analyzer (Roche Diagnostics). Interassay CV for cystatin C was 3.1%, which was the same as reported for this assay by Delanaye et al. and similar to their results for the Siemens/Dade-Behring nephelometric immunoassay.30 Cholesterol, HDL, and LDL were analyzed using an enzymatic colorimetric method, and plasma triglycerides were analyzed with a colorimetric method.
BP was measured in a seated position after 2 minutes of resting using an automatic device (model UA 799; A&D). BP was measured three times with 1 minute between measurements. The average of the second and third measurements was used in the analyses. Anthropometric measures were obtained by trained personnel. BMI was calculated as body weight in kilograms divided by the square of height measured in meters.
The investigated creatinine-based equations in this study were the recalibrated four-variable MDRD equation and the CKD-EPI equation.31,32 The equation used for estimation of GFR based on cystatin C was developed by Rule et al.33 (Table 6). A previous validation found bias <3.0 ml/min per 1.73 m2 for all these three equations in the RENIS-T6 cohort.34
Equations for estimating GFR
The following covariates were obtained from the questionnaire. Smoking status was divided into current smokers or nonsmokers (former smokers were included in the group of nonsmokers). They were registered with a family history of cardiovascular disease if a first-degree relative had had MI, angina pectoris, cerebral stroke, and/or cerebral hemorrhage. Family history of early MI was defined as a first-degree relative with MI before the age of 60 years. Physically active participants were defined as those performing hard physical activity for >1 h/wk (with prominent perspiration or breathlessness) and/or light physical activity for >3 h/wk (without prominent perspiration or breathlessness). Physical activity during working hours was not included in the analysis.
Statistical Analysis
The participants were divided into gender-specific quintiles for mGFR, and differences across the quintiles for the cardiovascular risk factors were analyzed with analysis of covariance for the continuous variables and logistic regression for the dichotomous variables. The analyses were adjusted for age and current use of angiotensin converting enzyme inhibitors or angiotensin receptor blockers. Analyses of serum lipids were also adjusted for lipid-lowering medication. Linear and quadratic trends across the quintiles were tested.
In multiple linear regression analyses, the association between mGFR as a dependent variable and cardiovascular risk factors as independent variables was explored. Each risk factor was first examined separately in a model adjusting only for gender, age, and the use of angiotensin converting enzyme inhibitors or angiotensin receptor blockers and then in a model adjusting for all of the other investigated risk factors as well. Systolic BP was not included in the full model because of colinearity with diastolic BP. The regression coefficients for continuous independent variables were given per SD of the variable. Interactions between gender and all of the other independent variables were tested in the full model. Possible nonlinear relationships were explored with generalized additive models using smoothing with cubic splines.35
The same analyses were performed for eGFR as outlined above for mGFR, but adjustment for mGFR was added to the models to study possible non–GFR-related effects. In addition, the fully adjusted models were analyzed with errors-in-variables regression assuming reduced reliability for mGFR.36 Reliability was estimated as 1 − (noise variance/total variance). With noise variance calculated from the within-subject CV of 5.5% found by others in repeated measurements of plasma iohexol clearance and total variance set at the population variance of mGFR in the RENIS-T6 cohort, the reliability of mGFR was estimated at 0.87.37,38 to allow for the possibility of a higher within-subject CV in our study, we performed the errors-in-variables regression analyses with a mGFR reliability set at 0.70.
Based on the Framingham risk score, the 10-year risk for coronary heart disease was estimated. The Framingham risk score includes age, gender, BP, total serum cholesterol, HDL cholesterol, smoking status, diabetes, and electrocardiogram-left ventricular hypertrophy as predictors.39 Information about electrocardiogram-left ventricular hypertrophy was not available for this analysis, and we set all participants as having no left ventricular hypertrophy as specified in the risk score algorithm.39 All participants were set as not having diabetes because diabetics were excluded from the study.
P < 0.05 was considered statistically significant. The statistical analyses were run using SPSS software version 16.0 (SPSS, Chicago, IL) for Windows. Generalized additive models were analyzed using PROC GAM in SAS version 9.2. (SAS Institute, Cary, NC) and errors-in-variables regression using STATA 11.0 (StataCorp, College Station, TX).
This study was approved by the Norwegian Data Inspectorate and the Regional Ethics Committee of North Norway. All subjects gave written consent.
DISCLOSURES
None.
Acknowledgments
We thank Britt-Ann Winther Eilertsen, Bjørg Skog Høgset, Saskia van Heusden, and the rest of the staff at the Clinical Resarch Unit (University Hospital of North Norway) for performing the study; Harald Strand and the staff at the Department of Medical Biochemistry (University Hospital of North Norway) for HPLC analyses of iohexol; Inger Sperstad and Ingrid Dorthea Sandstad (Clinical Research Centre, University Hospital of North Norway) for database support; and Tom Wilsgaard, Sriharan Sivasingarajah, and Kurt Jøran Nyland (Institute of Community Medicine, University of Tromsø) for identifying eligible subjects from the Tromsø 6 cohort.
Footnotes
Published online ahead of print. Publication date available at www.jasn.org.
See related editorial, “Cystatin C Is More than GFR, and This May Be a Good Thing,” on pages 795–797.
- Copyright © 2011 by the American Society of Nephrology