Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Med., 17 November 2025

Sec. Nephrology

Volume 12 - 2025 | https://doi.org/10.3389/fmed.2025.1670059

Association between the difference in cystatin C and creatinine-based eGFR and risks of multiple cardiovascular diseases: a prospective cohort study


Zhiyu Qiao&#x;Zhiyu QiaoXinyi Liu&#x;Xinyi LiuHao LiuHao LiuSuwei ChenSuwei ChenChengnan LiChengnan LiYipeng GeYipeng GeHaiou HuHaiou HuJunming Zhu*&#x;Junming Zhu*‡
  • Department of Cardiovascular Surgery, Beijing Aortic Disease Center, Beijing Anzhen Hospital of Capital Medical University, Beijing, China

Background: The difference between cystatin C- and creatinine-based estimated glomerular filtration rate (eGFRdiff) is closely associated with various adverse outcomes. This study aims to comprehensively evaluate the association between eGFRdiff, all-cause mortality, and the risk of multiple cardiovascular-related diseases.

Methods: This study analyzed data from 297,140 participants in the UK Biobank to assess the association between eGFRdiff, mortality, and the incidence of multiple cardiovascular-related diseases. eGFRdiff was classified into three groups: negative (< −15 mL/min/1.73 m2), intermediate (−15 to 15 mL/min/1.73 m2), and positive (≥ 15 mL/min/1.73 m2). Cox proportional hazards regression models were used to evaluate this association, while various sensitivity analyses were performed to assess its robustness.

Results: During a mean follow-up of 13.1 years, the positive eGFRdiff group exhibited significantly lower mortality, cardiovascular disease (CVD) incidence, and the occurrence of CVD-related conditions. In the fully adjusted model, participants in the negative eGFRdiff group had a hazard ratio of 1.44 (95% confidence interval [CI], 1.40–1.49) for all-cause mortality, 1.49 (95% CI, 1.41–1.59) for CVD incidence, and 1.25 (95% CI, 1.22–1.27) for CVD mortality. The risk of all 10 CVD-related conditions was also significantly higher in the negative group, whereas the positive group exhibited significantly lower risks. For every 10 mL/min/1.73 m2 increase in eGFRdiff, the incidence of various diseases decreased by approximately 10–19%.

Conclusion: eGFRdiff is significantly associated with increased risks of mortality, CVD incidence, and multiple CVD-related conditions. These findings underscore the critical need for developing targeted prevention strategies, particularly for populations with reduced eGFRdiff.

Introduction

Chronic kidney disease is a global public health issue, affecting more than 10% of adults worldwide (1). Impaired kidney function is associated with a high risk of cardiovascular disease (CVD) and all-cause mortality (2). Numerous studies have suggested that non-dialysis-dependent kidney dysfunction is closely linked to heart failure, atrial fibrillation, and the incidence of CVD in asymptomatic populations (3, 4). Therefore, identifying high-risk individuals with kidney dysfunction and implementing early interventions are crucial for preventing the onset and progression of cardiovascular disease (5).

Estimating glomerular filtration rate (eGFR) using serum creatinine or cystatin C is a widely adopted method for assessing kidney function in clinical practice (6). However, in recent years, the substantial intra-individual differences between cystatin C-based eGFR (eGFRcys) and creatinine-based eGFR (eGFRcr) have been increasingly recognized and reported to be associated with various adverse outcomes, including mortality, end-stage kidney disease, and hospitalization (7, 8). These differences may be influenced by multiple non-renal factors, such as muscle mass, lifestyle, and chronic diseases (9, 10). The eGFR difference (eGFRdiff), defined as eGFRcys minus eGFRcr, has recently been proposed as a marker of overall health status (11). Previous studies have shown that a more negative eGFRdiff is associated with several adverse cardiovascular outcomes, including heart failure, atrial fibrillation, and atherosclerotic cardiovascular disease (1214). Despite these associations, no large-scale cohort study has systematically analyzed the relationship between eGFRdiff and the incidence of comprehensive CVD as well as a wide range of CVD-related conditions.

This study aims to investigate the association between eGFRdiff and the incidence, mortality, and various CVD-related diseases using prospective data from the UK Biobank.

Materials and methods

Study design and participants

We drew upon data from the UK Biobank study (Application Number 145937), a comprehensive prospective cohort comprising more than 500,000 individuals aged 37–73 years, recruited from 22 assessment centers throughout the United Kingdom between 2006 and 2010. Participants contributed extensive health-related data via a touchscreen questionnaire, which encompassed demographics, socio-economic status, lifestyle habits, and medical conditions. The study’s methodology and data collection processes have been extensively documented in previous publications (15).

The cohort initially included 466,571 participants with complete eGFR data at baseline. After excluding individuals with a history of cardiovascular disease at baseline, as well as those with missing demographic information or other covariate data, the final analysis included 297,140 participants (Supplementary Figure 1).

Main exposure

In this study, the primary exposure of interest was the non-race-based eGFRdiff. It was calculated using baseline serum cystatin C and creatinine levels, applying the CKD-EPI (Chronic Kidney Disease Epidemiology Collaboration) equations. Specifically, eGFRdiff was derived by subtracting the non-race-based eGFRcr calculated from creatinine from the eGFRcys calculated from cystatin C. Participants were categorized into three groups based on eGFRdiff values: negative (< −15 mL/min/1.73 m2), intermediate (−15 to 15 mL/min/1.73 m2), and positive (≥ 15 mL/min/1.73 m2). Additionally, eGFRdiff was analyzed as a continuous variable per 10 mL/min/1.73 m2 increment. In a sensitivity analysis, we used the race-related eGFRdiff as the primary exposure, defined as eGFRcys minus the race-related eGFRcr.

Assessment of covariates

Covariates included demographic characteristics, baseline medical history, lifestyle factors, chronic inflammation markers, and laboratory biochemical test indicators. Demographic characteristics included age, sex, education level, self-reported race, townsend deprivation index, employment status, body mass index (BMI), handgrip strength (HGS) and appendicular skeletal muscle mass (ASM). HGS was measured in both hands in kilograms, and the average grip strength of both hands was calculated. Body composition was assessed at the baseline visit using a single-frequency segmental body composition analyzer. Muscle and fat mass were derived using bioelectrical impedance analysis. ASM was calculated as the sum of lean soft tissue mass in the upper and lower limbs. Finally, sarcopenia was identified based on the criteria established by the Foundation for the National Institutes of Health (FNIH) Sarcopenia Project. Baseline medical history also covered the presence of chronic respiratory disease, chronic liver disease, hypertension, diabetes, and dyslipidemia.

Lifestyle factors were evaluated using six components aligned with World Health Organization guidelines: dietary habits, sleep patterns (classified as healthy, moderate, or unhealthy), physical activity levels (high, moderate, or low), sedentary behavior (low, moderate, or high), and history of smoking and alcohol consumption. Comprehensive details regarding the assessment of these lifestyle factors can be found in Supplementary Tables 1, 2. For each lifestyle factor, a score of 0 was assigned to an unhealthy level and 1 to a healthy level (similarly, for sleep patterns, a score of 1 was assigned to a moderate sleep health pattern and 2 to the healthiest sleep pattern). The six lifestyle factors were then summed to generate a composite healthy lifestyle score.

To account for the role of chronic inflammation, we also analyzed various inflammatory markers, including white blood cell (WBC) count, platelet count, C-reactive protein (CRP) levels, neutrophil count, lymphocyte count, and the neutrophil-to-lymphocyte ratio (NLR). Additionally, we calculated the low-grade chronic inflammation score (INFLA score) as a comprehensive measure of individual inflammatory status. Based on previous studies, the INFLA score, which is closely associated with multiple diseases (16), integrates four inflammatory markers: CRP, WBC count, platelet count, and NLR. To compute the INFLA score, each inflammatory marker was log-transformed. Biomarker levels within the highest deciles (7th–10th) were assigned scores ranging from + 1 to + 4, whereas those within the lowest deciles (1st–4th) were assigned scores from −4 to −1 (17). The resulting INFLA score ranged from −16 to + 16, with higher scores indicating elevated levels of low-grade inflammation.

Finally, we included multiple laboratory biochemical test indicators, including serum albumin, high-density lipoprotein (HDL-C), low-density lipoprotein (LDL-C), triglycerides, and the urine albumin-to-creatinine ratio (UACR). The UACR was used to determine the presence of albuminuria in participants.

Outcomes

The primary outcomes of this study included cardiovascular disease (CVD) incidence, CVD mortality, and all-cause mortality. Secondary outcomes encompassed all CVD-related components, including stroke, heart failure, atrial fibrillation, valvular heart disease, coronary atherosclerotic heart disease, aortic aneurysm, peripheral artery disease, deep vein thrombosis, pulmonary embolism, and arterial embolism. All disease diagnoses were determined based on death registries, primary care records, hospitalization data, and self-reported diagnoses. Outcomes were classified using ICD-9 and ICD-10 codes (International Classification of Diseases, Ninth and Tenth Revisions). The follow-up period extended from the baseline assessment (2006–2010) to the earliest occurrence of a relevant disease diagnosis, death, loss to follow-up, or the end of the follow-up period, whichever occurred first.

Statistical analysis

The baseline characteristics of participants were presented as means with standard deviations (SD) for continuous variables and as proportions for categorical variables. Comparisons of continuous variables were performed using analysis of variance (ANOVA), while differences in categorical variables were assessed using the χ2-test.

Pearson correlation was used to examine the relationship between eGFR difference (eGFRdiff) and clinical parameters. Cox proportional hazards regression models were applied to calculate the hazard ratios (HRs) and 95% confidence intervals (CIs) for eGFRdiff with respect to primary and secondary outcomes. Three multivariable-adjusted models were developed. Model 1 adjusted for age, sex, racial background, educational level, occupational status, Townsend deprivation index, and body mass index. Model 2 further adjusted for healthy lifestyle score and comorbidities (chronic respiratory disease, chronic liver disease, hypertension, diabetes, and dyslipidemia). Model 3 additionally adjusted for laboratory measurements (INFLA score, serum albumin, HDL-C, LDL-C, triglycerides, UACR, and eGFRcr). Restricted cubic splines (RCS) were used to explore the association between eGFRdiff and the risk of outcome events. In addition, subgroup analyses were performed to examine potential differences, stratifying participants by sex (male and female), age (≥ 60 and < 60 years), BMI (normal BMI: 18.5-24.9 and abnormal BMI), and the presence of comorbidities at baseline.

In this study, several sensitivity analyses were conducted to assess the robustness and consistency of the models. First, we used race-related eGFR difference as the primary exposure. Second, we adjusted for eGFRcys or eGFRcr–cys instead of eGFRcr. Third, the analysis was stratified by eGFRcr, eGFRcys, and eGFRcr–cys categories: ≥ 90, 60–89, 45–59, and < 45 mL/min/1.73 m2. Fourth, eGFRdiff was further categorized into tertiles. Fifth, we further adjusted for participants with sarcopenia. Sixth, we adjusted for the components of lifestyle scores and the INFLA score rather than using the two scores as a whole. Seventh, individuals with comorbidities (chronic liver disease, chronic respiratory disease, hypertension, diabetes, dyslipidemia, and sarcopenia) were excluded to examine the impact of comorbidities on the outcomes. Eighth, events occurring within the first 3 years of follow-up were excluded to minimize potential bias from early events. Finally, multiple imputation was used to address missing covariate data and evaluate the impact of incomplete variables on the results.

All statistical analyses were performed using R software version 4.4.1 (R Foundation for Statistical Computing). P-values were two-sided, with a significance level set at P < 0.05.

Results

Baseline characteristics of participants

The baseline characteristics of the participants are presented in Table 1. This study included 297,140 individuals with a mean (SD) age of 58.0 (7.6) years. Among them, 135,954 (45.8%) were male, and 284,119 (95.6%) were White. The mean non-race-based eGFR difference (eGFRdiff) was −0.5 ± 12.9 mL/min/1.73 m2. The distribution of race-based eGFRdiff is shown in Supplementary Table 3. Compared to participants with lower eGFRdiff, those with higher eGFRdiff were more likely to be female, younger, more highly educated, have a healthier lifestyle, exhibit lower levels of chronic inflammation, and report a lower prevalence of other diseases. In the race-related eGFRdiff subgroup analysis, we similarly observed consistent findings.

TABLE 1
www.frontiersin.org

Table 1. Baseline characteristics stratified by categories of non-race-based eGFRdiff.

In the correlation analysis, we identified significant associations between non-race-based eGFRdiff and multiple factors (Supplementary Table 4). Specifically, eGFRdiff exhibited significant negative correlations with age (γ = −0.19; P < 0.001), BMI (γ = −0.22; P < 0.001), history of alcohol consumption (γ = −0.10; P < 0.001), dyslipidemia (γ = −0.16; P < 0.001), and the INFLA score for chronic inflammation (γ = −0.17; P < 0.001). In contrast, significant positive correlations were observed with physical activity (γ = 0.11; P < 0.001), healthy lifestyle score (γ = 0.06; P < 0.001), and HDL-C levels (γ = 0.17; P < 0.001). Similar patterns of association were also observed in the race-related eGFRdiff subgroup analysis (Supplementary Table 5).

Association between eGFRdiff and mortality and incident CVD

During a mean follow-up period of 13.1 years, the incidence of cardiovascular disease (CVD), CVD-related mortality, and all-cause mortality were 43,315 (14.6%), 4,634 (1.6%), and 19,289 (6.5%), respectively. As shown in Table 2, compared to the negative group of non-race-based eGFRdiff, both the moderate and positive groups exhibited significantly lower rates of mortality and CVD incidence. Specifically, the CVD incidence in the positive group (6.36 per 1,000 person-years) was significantly lower than that in the negative group (15.51 per 1,000 person-years). A similar trend was observed across various CVD-related conditions, with the most notable finding being the significantly lower incidence of coronary atherosclerotic heart disease in the positive group (2.63 per 1,000 person-years) compared to the negative group (6.52 per 1,000 person-years). As presented in Supplementary Table 6, a similar trend was observed in the groups stratified by race-related eGFRdiff.

TABLE 2
www.frontiersin.org

Table 2. Incidence of cardiovascular disease and its components across three categories of non-race-based eGFRdiff levels.

Table 3 presents the association between non-race-based eGFRdiff and mortality and CVD incidence. In the fully adjusted model, both the negative and positive eGFRdiff groups were associated with the primary outcomes compared to the moderate eGFRdiff group. Specifically, individuals in the negative eGFRdiff group had a higher risk of CVD incidence (HR = 1.22, 95% CI: 1.19–1.25), CVD mortality (HR = 1.47, 95% CI: 1.37–1.56), and all-cause mortality (HR = 1.41, 95% CI: 1.36–1.46). Conversely, individuals in the positive eGFRdiff group exhibited a lower risk of CVD incidence (HR = 0.86, 95% CI: 0.81–0.90), CVD mortality (HR = 0.67, 95% CI: 0.55–0.81), and all-cause mortality (HR = 0.77, 95% CI: 0.70–0.84). When analyzed as a continuous variable, each 10 mL/min/1.73 m2 increase in eGFRdiff was associated with a 10–19% reduction in the risk of mortality and CVD incidence. A similar trend was observed in Supplementary Table 7, which presents the association between race-related eGFRdiff and these outcomes.

TABLE 3
www.frontiersin.org

Table 3. Cardiovascular disease risk and mortality stratified by non-race-based eGFRdiff levels.

Table 4 and Supplementary Table 8 present the associations between both non-race-based and race-related eGFRdiff and the incidence of cardiovascular-related diseases. Compared to the moderate eGFRdiff group, individuals in the negative eGFRdiff group exhibited a significantly higher risk of multiple cardiovascular diseases, whereas those in the positive eGFRdiff group had a significantly lower risk. Among these, heart failure showed the most pronounced association, with an incidence risk ratio of 1.54 (95% CI: 1.46–1.62) in the negative group and 0.71 (95% CI: 0.60–0.84) in the positive group. Notably, among the 10 cardiovascular-related diseases analyzed, aortic aneurysm and deep vein thrombosis did not show significant differences in incidence in the positive eGFRdiff group, while significant associations were observed in the negative group. However, when eGFRdiff was treated as a continuous variable, significant associations were observed across all cardiovascular-related diseases, with each 10 mL/min/1.73 m2 increase in eGFRdiff corresponding to a 9–23% reduction in disease incidence.

TABLE 4
www.frontiersin.org

Table 4. Risk of cardiovascular disease subtypes stratified by non-race-based eGFRdiff levels.

Given the strong associations between eGFRdiff as a continuous variable and various diseases, we further explored its linear relationship using restricted cubic spline analysis. As shown in Figure 1 and Supplementary Figure 2, non-race-based eGFRdiff exhibited a significant linear association with CVD incidence, all-cause mortality, and several conditions including stroke, heart failure, and atrial fibrillation. However, no significant linear association was observed with CVD mortality. Similarly, as illustrated in Supplementary Figure 3, race-related eGFRdiff demonstrated comparable associations.

FIGURE 1
Three line graphs display hazard ratios (HR) with 95% confidence intervals for different health risks against eGFR levels. Graph A shows a steep decrease in incident CVD risk. Graph B presents a decline in CVD mortality risk with less pronounced nonlinearity. Graph C illustrates a significant drop in all-cause mortality risk. Each graph includes a horizontal dashed line indicating HR of 1, with significance values noted as p overall less than 0.001, and p for nonlinear varying per graph.

Figure 1. Restricted cubic spline plots of mortality and cardiovascular disease incidence based on non-race-based eGFRdiff. This analysis adjusted for age, sex, racial background, educational level, occupational status, Townsend deprivation index, body mass index, healthy lifestyle score, comorbidities (chronic respiratory disease, chronic liver disease, hypertension, diabetes, and dyslipidemia), and laboratory measurements (eGFRcr, INFLA score, serum albumin, HDL-C, LDL-C, triglycerides, and UACR). eGFRdiff, the difference between cystatin C–based estimated glomerular filtration rate and creatinine-based estimated glomerular filtration rate; HR, hazard ratio; CI, confidence interval; INFLA score, Low-grade chronic inflammation score; HDL-C, high-density lipoprotein cholesterol; LDL-C, low-density lipoprotein cholesterol; UACR, Urinary albumin-creatinine ratio; CVD, Cardiovascular disease.

Sensitivity analysis and subgroup analysis

To validate the robustness of the results, we conducted various sensitivity analyses. The associations between both types of eGFRdiff and mortality, CVD incidence, and cardiovascular-related disease incidence remained significant after adjusting for eGFRcys or eGFRcr–cys instead of eGFRcr (Supplementary Tables 9, 10) and when stratifying participants based on renal function using eGFRcr, eGFRcys, or eGFRcr–cys (Supplementary Tables 11, 12). These associations persisted even in additional sensitivity analyses, including using tertiles of eGFRdiff and excluding participants who experienced relevant events within the first 3 years of follow-up (Supplementary Tables 13, 14).

Finally, we examined the effects of both types of eGFRdiff on mortality and cardiovascular disease incidence across various subgroups (Supplementary Tables 1517). In all subgroups, eGFRdiff remained significantly associated with mortality and cardiovascular disease incidence, and no significant interactions were observed.

Discussion

This study is a large prospective cohort analysis utilizing data from the UK Biobank to examine the association between eGFRdiff and mortality, as well as the incidence of CVD and its related conditions. Compared to participants with a positive eGFRdiff, those with a negative eGFRdiff exhibited significantly higher mortality rates, CVD incidence, and the occurrence of 10 CVD-related diseases. Moreover, when eGFRdiff was treated as a continuous variable, higher eGFRdiff values were associated with significantly lower mortality and CVD incidence, independent of kidney function levels. These trends remained consistent across various sensitivity analyses and restricted cubic spline analyses.

Previous studies have reported associations between eGFRdiff and several adverse cardiovascular events. For instance, Debbie et al. observed in a cohort of 4,512 patients with chronic kidney disease that a negative baseline eGFRdiff was associated with an increased risk of heart failure (12). Similarly, Ga et al. reported a strong correlation between eGFRdiff and the incidence of atrial fibrillation (13). Furthermore, in a cohort of 9,092 participants from the Systolic Blood Pressure Intervention Trial (SPRINT), eGFRdiff was found to be closely associated with mortality risk (18). However, most of these studies focused on single cardiovascular outcomes or were conducted in relatively small cohorts. In this study, we emphasize that eGFRdiff is significantly associated with various CVD-related adverse events, as well as overall CVD incidence and mortality. Our findings provide additional evidence supporting the relationship between eGFRdiff and the occurrence of CVD and its related diseases. Interestingly, these associations remained significant even after adjusting for eGFRcr, eGFRcys, eGFRcr–cys, or stratifying by kidney function categories. This further reinforces the notion that a more negative eGFRdiff is a substantial risk factor for both mortality and CVD incidence. Notably, compared to the intermediate eGFRdiff group, we observed an increased risk of aortic aneurysm and deep vein thrombosis in the negative eGFRdiff group, similar to the elevated risks observed for other CVD conditions. However, this association was not significant in the positive eGFRdiff group. This discrepancy may be due to the relatively low number of cases for these two conditions within the cohort. Further large-scale cohort studies are needed to confirm this association.

Due to the complexity of eGFR measurement across different populations, efforts have been made in recent years to eliminate race-based kidney function assessment (19). In this context (20), the CKD-EPI equations have been updated to provide race-independent estimates for both eGFRcr and eGFRcys. Consequently, variations in eGFRdiff have emerged due to differences in calculation methods. The UK Biobank includes participants with a diverse self-reported racial background. In this study, we systematically analyzed two classifications of eGFRdiff and found that both were significantly associated with mortality and CVD incidence. This finding has important implications for the prevention of CVD, particularly among traditionally understudied or high-risk racial groups.

Several potential mechanisms may underlie the association between eGFRdiff and the development of CVD. First, both creatinine and cystatin C are influenced by non-renal factors. Physical activity, chronic diseases, and muscle mass are key determinants of creatinine levels, while obesity, smoking, and steroid use are considered major non-renal factors affecting cystatin C levels (9, 10, 21). These influences contribute to the discrepancy between eGFRcys and eGFRcr. A more negative eGFRdiff indicates a lower eGFRcys and a higher eGFRcr, which may reflect poorer overall health status. Second, Grubb et al. hypothesized that individuals with a lower eGFRdiff tend to have higher BMI and elevated levels of inflammatory markers, a condition referred to as “shrunken pore syndrome” (22). Proteomic studies have shown that certain inflammatory proteins, such as interleukin-6 and osteoprotegerin, accumulate in individuals with shrunken pore syndrome (23). These inflammatory mediators are believed to contribute to endothelial damage, promote inflammation, and accelerate atherosclerosis—key pathogenic factors in CVD development (24). In our study, correlation analyses between eGFRdiff and various clinical indicators revealed that eGFRdiff was significantly positively associated with overall healthy lifestyle factors and negatively associated with chronic inflammation levels and preexisting disease history. These findings further support the proposed mechanisms linking eGFRdiff with CVD risk.

Our study has several notable strengths. First, it is a large-scale prospective cohort study that provides comprehensive and detailed data on the association between eGFRdiff, mortality, and cardiovascular disease incidence. Second, this study systematically analyzes the relationship between eGFRdiff and 10 CVD-related conditions, addressing the limitations of previous research that primarily focused on single CVD outcomes. Additionally, we incorporated a wide range of potential confounders and constructed lifestyle and chronic inflammation scores, which enhance the robustness and generalizability of our findings.

However, our study has several limitations. First, as an observational study, it cannot establish a causal relationship between eGFRdiff, mortality, and CVD development. Second, we were unable to fully adjust for potential residual confounders, such as the impact of unmeasured comorbidities. Third, eGFRdiff in this study was calculated based on baseline measurements, and we could not assess its changes over time. Fourth, because cystatin C is not routinely measured in all clinical settings, its applicability and generalizability may be limited. Finally, since the UK Biobank primarily consists of a predominantly White adult population, the generalizability of our findings to other racial and ethnic groups remains limited.

Conclusion

In conclusion, this study found that eGFRdiff is closely associated with mortality, CVD incidence, and the risk of multiple CVD-related conditions. This finding underscores the importance of developing targeted prevention strategies, particularly for individuals with lower eGFRdiff. However, further comprehensive evaluations are needed to determine the clinical utility of eGFRdiff as a predictive marker in medical practice.

Data availability statement

The original contributions presented in this study are included in the article/Supplementary material, further inquiries can be directed to the corresponding authors.

Ethics statement

The studies involving humans were approved by the North West Multi-center Research Ethics Committee (Ref:11/NW/0382). The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.

Author contributions

ZQ: Conceptualization, Writing – original draft. XL: Data curation, Formal analysis, Writing – original draft, Writing – review & editing. HL: Investigation, Methodology, Writing – original draft. SC: Formal analysis, Project administration, Writing – original draft. CL: Resources, Software, Writing – review & editing. YG: Supervision, Validation, Writing – review & editing. HH: Validation, Visualization, Writing – review & editing. JZ: Funding acquisition, Resources, Supervision, Writing – review & editing.

Funding

The author(s) declare financial support was received for the research and/or publication of this article. This research was supported by funding from the National Science Foundation of China (NSFC 82170490).

Acknowledgments

We thank all participants and staff of the UK Biobank.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The authors declare that Generative AI was used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmed.2025.1670059/full#supplementary-material

Abbreviations

AA, Aortic aneurysm; ASM, Appendicular skeletal muscle mass; BMI, Body mass index; CI, Confidence interval; CKD-EPI, Chronic Kidney Disease Epidemiology Collaboration; CRP, C-reactive protein; CVD, Cardiovascular disease; eGFR, Estimated glomerular filtration rate; eGFRcr, Creatinine-based estimated glomerular filtration rate; eGFRcys, Cystatin C-based estimated glomerular filtration rate; eGFRdiff, Difference between cystatin C- and creatinine-based estimated glomerular filtration rate; FNIH, Foundation for the National Institutes of Health; HDL-C, High-density lipoprotein cholesterol; HGS, Handgrip strength; HR, Hazard ratio; ICD, International Classification of Diseases; INFLA score, Low-grade chronic inflammation score; LDL-C, Low-density lipoprotein cholesterol; NLR, Neutrophil-to-lymphocyte ratio; NSFC, National Science Foundation of China; RCS, Restricted cubic splines; SD, Standard deviation; SPRINT, Systolic Blood Pressure Intervention Trial; UACR, Urine albumin-to-creatinine ratio; WBC, White blood cell.

References

1. GBD Chronic Kidney Disease Collaboration. Global, regional, and national burden of chronic kidney disease, 1990-2017: a systematic analysis for the Global burden of disease study 2017. Lancet. (2020) 395:709–33. doi: 10.1016/S0140-6736(20)30045-3

PubMed Abstract | Crossref Full Text | Google Scholar

2. Fox C, Matsushita K, Woodward M, Bilo H, Chalmers J, Heerspink H, et al. Associations of kidney disease measures with mortality and end-stage renal disease in individuals with and without diabetes: a meta-analysis. Lancet. (2012) 380:1662–73. doi: 10.1016/S0140-6736(12)61350-6

PubMed Abstract | Crossref Full Text | Google Scholar

3. Grundy S, Stone N, Bailey A, Beam C, Birtcher K, Blumenthal R, et al. 2018 Aha/Acc/Aacvpr/Aapa/Abc/Acpm/Ada/Ags/Apha/Aspc/Nla/Pcna guideline on the management of blood cholesterol: a report of the American college of cardiology/American heart association task force on clinical practice guidelines. J Am Coll Cardiol. (2019) 73:3168–209. doi: 10.1016/j.jacc.2018.11.002

PubMed Abstract | Crossref Full Text | Google Scholar

4. Matsushita K, van der Velde M, Astor B, Woodward M, Levey A, de Jong P, et al. Association of estimated glomerular filtration rate and albuminuria with all-cause and cardiovascular mortality in general population cohorts: a collaborative meta-analysis. Lancet. (2010) 375:2073–81. doi: 10.1016/S0140-6736(10)60674-5

PubMed Abstract | Crossref Full Text | Google Scholar

5. Visseren F, Mach F, Smulders Y, Carballo D, Koskinas K, Bäck M, et al. 2021 Esc guidelines on cardiovascular disease prevention in clinical practice. Eur Heart J. (2021) 42:3227–337. doi: 10.1093/eurheartj/ehab484

PubMed Abstract | Crossref Full Text | Google Scholar

6. Inker L, Eneanya N, Coresh J, Tighiouart H, Wang D, Sang Y, et al. New creatinine- and cystatin C-based equations to estimate Gfr without race. N Engl J Med. (2021) 385:1737–49. doi: 10.1056/NEJMoa2102953

PubMed Abstract | Crossref Full Text | Google Scholar

7. Chen D, Shlipak M, Scherzer R, Bauer S, Potok O, Rifkin D, et al. Association of intraindividual difference in estimated glomerular filtration rate by creatinine Vs cystatin C and end-stage kidney disease and mortality. JAMA Netw Open. (2022) 5:e2148940. doi: 10.1001/jamanetworkopen.2021.48940

PubMed Abstract | Crossref Full Text | Google Scholar

8. Carrero J, Fu E, Sang Y, Ballew S, Evans M, Elinder C, et al. Discordances between Creatinine- and Cystatin C-based estimated Gfr and adverse clinical outcomes in routine clinical practice. Am J Kidney Dis. (2023) 82:534–42. doi: 10.1053/j.ajkd.2023.04.002

PubMed Abstract | Crossref Full Text | Google Scholar

9. Inker L, Titan S. Measurement and estimation of Gfr for use in clinical practice: core curriculum 2021. Am J Kidney Dis. (2021) 78:736–49. doi: 10.1053/j.ajkd.2021.04.016

PubMed Abstract | Crossref Full Text | Google Scholar

10. Foster M, Levey A, Inker L, Shafi T, Fan L, Gudnason V, et al. Non-Gfr determinants of low-molecular-weight serum protein filtration markers in the elderly: ages-kidney and mesa-kidney. Am J Kidney Dis. (2017) 70:406–14. doi: 10.1053/j.ajkd.2017.03.021

PubMed Abstract | Crossref Full Text | Google Scholar

11. Potok O, Katz R, Bansal N, Siscovick D, Odden M, Ix J, et al. The difference between Cystatin C- and Creatinine-based estimated Gfr and incident frailty: an analysis of the Cardiovascular health study (Chs). Am J Kidney Dis. (2020) 76:896–8. doi: 10.1053/j.ajkd.2020.05.018

PubMed Abstract | Crossref Full Text | Google Scholar

12. Chen D, Shlipak M, Scherzer R, Bansal N, Potok O, Rifkin D, et al. Association of intra-individual differences in estimated Gfr by creatinine versus cystatin C with incident heart failure. Am J Kidney Dis. (2022) 80: 762–72.e1. doi: 10.1053/j.ajkd.2022.05.011

PubMed Abstract | Crossref Full Text | Google Scholar

13. Heo G, Koh H, Jung C, Park J, Han S, Yoo T, et al. Difference between estimated Gfr based on Cystatin C Versus creatinine and incident atrial fibrillation: a cohort study of the UK Biobank. Am J Kidney Dis. (2024) 83:729–38 e1. doi: 10.1053/j.ajkd.2023.11.004

PubMed Abstract | Crossref Full Text | Google Scholar

14. Chen D, Lu K, Scherzer R, Lees J, Rutherford E, Mark P, et al. Cystatin C- and Creatinine-based estimated Gfr differences: prevalence and predictors in the UK Biobank. Kidney Med. (2024) 6:100796. doi: 10.1016/j.xkme.2024.100796

PubMed Abstract | Crossref Full Text | Google Scholar

15. Cox N. UK Biobank shares the promise of big data. Nature. (2018) 562:194–5. doi: 10.1038/d41586-018-06948-3

PubMed Abstract | Crossref Full Text | Google Scholar

16. Shi H, Schweren L, Ter Horst R, Bloemendaal M, van Rooij D, Vasquez A, et al. Low-grade inflammation as mediator between diet and behavioral disinhibition: a UK Biobank study. Brain Behav Immun. (2022) 106:100–10. doi: 10.1016/j.bbi.2022.07.165

PubMed Abstract | Crossref Full Text | Google Scholar

17. Pounis G, Bonaccio M, Di Castelnuovo A, Costanzo S, de Curtis A, Persichillo M, et al. Polyphenol intake is associated with low-grade inflammation, using a novel data analysis from the moli-sani study. Thromb Haemost. (2016) 115:344–52. doi: 10.1160/TH15-06-0487

PubMed Abstract | Crossref Full Text | Google Scholar

18. Potok O, Ix J, Shlipak M, Katz R, Hawfield A, Rocco M, et al. The difference between Cystatin C- and Creatinine-based estimated Gfr and associations with frailty and adverse outcomes: a cohort analysis of the systolic blood pressure intervention trial (Sprint). Am J Kidney Dis. (2020) 76:765–74. doi: 10.1053/j.ajkd.2020.05.017

PubMed Abstract | Crossref Full Text | Google Scholar

19. Delgado C, Baweja M, Crews D, Eneanya N, Gadegbeku C, Inker L, et al. A unifying approach for Gfr estimation: recommendations of the Nkf-Asn task force on reassessing the inclusion of race in diagnosing kidney disease. Am J Kidney Dis. (2022) 79:268–88.e1. doi: 10.1053/j.ajkd.2021.08.003

PubMed Abstract | Crossref Full Text | Google Scholar

20. Hundemer G, White C, Norman P, Knoll G, Tangri N, Sood M, et al. Performance of the 2021 race-free Ckd-Epi Creatinine- and Cystatin C-based estimated Gfr equations among kidney transplant recipients. Am J Kidney Dis. (2022) 80:462–72.e1. doi: 10.1053/j.ajkd.2022.03.014

PubMed Abstract | Crossref Full Text | Google Scholar

21. Levey A, Inker L, Coresh J. Gfr estimation: from physiology to public health. Am J Kidney Dis. (2014) 63:820–34. doi: 10.1053/j.ajkd.2013.12.006

PubMed Abstract | Crossref Full Text | Google Scholar

22. Grubb A. Shrunken pore syndrome - a common kidney disorder with high mortality. diagnosis, prevalence, pathophysiology and treatment options. Clin Biochem. (2020) 83:12–20. doi: 10.1016/j.clinbiochem.2020.06.002

PubMed Abstract | Crossref Full Text | Google Scholar

23. Xhakollari L, Jujic A, Molvin J, Nilsson P, Holm H, Bachus E, et al. Proteins linked to atherosclerosis and cell proliferation are associated with the shrunken pore syndrome in heart failure patients: shrunken pore syndrome and proteomic associations. Proteomics Clin Appl. (2021) 15:e2000089. doi: 10.1002/prca.202000089

PubMed Abstract | Crossref Full Text | Google Scholar

24. Almen M, Bjork J, Nyman U, Lindstrom V, Jonsson M, Abrahamson M, et al. Shrunken pore syndrome is associated with increased levels of atherosclerosis-promoting proteins. Kidney Int Rep. (2019) 4:67–79. doi: 10.1016/j.ekir.2018.09.002

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: estimated glomerular filtration rate, serum creatinine, cystatin C, cardiovascular disease, all-cause mortality

Citation: Qiao Z, Liu X, Liu H, Chen S, Li C, Ge Y, Hu H and Zhu J (2025) Association between the difference in cystatin C and creatinine-based eGFR and risks of multiple cardiovascular diseases: a prospective cohort study. Front. Med. 12:1670059. doi: 10.3389/fmed.2025.1670059

Received: 21 July 2025; Accepted: 31 October 2025;
Published: 17 November 2025.

Edited by:

Rodolfo Valtuille, University of Business and Social Sciences, Argentina

Reviewed by:

Santos Ángel Depine, Simón Bolívar University, Colombia
Luis Cruz-Llanos, National Cardiovascular Institute “Carlos Alberto Peschiera Carrillo,” Peru

Copyright © 2025 Qiao, Liu, Liu, Chen, Li, Ge, Hu and Zhu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Junming Zhu, YW56aGVuemptQGNjbXUuZWR1LmNu

These authors have contributed equally to this work

ORCID: Junming Zhu, orcid.org/0000-0002-2504-5681

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.