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

ORIGINAL RESEARCH article

Front. Med., 04 June 2025

Sec. Nephrology

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

A time-dependent predictive model for cardiocerebral vascular events in chronic hemodialysis patients: insights from a prospective study

Haowen Zhong,&#x;Haowen Zhong1,2Mengbi Zhang,&#x;Mengbi Zhang1,2Yingye Xie,&#x;Yingye Xie1,2Yuqin QinYuqin Qin3Na Xie,Na Xie1,2Yuqiu Ye,Yuqiu Ye1,2Heng Li,Heng Li2,4Hongquan PengHongquan Peng5Xun LiuXun Liu6Xiaoyan Su,
Xiaoyan Su1,2*Shaohong Li,
Shaohong Li2,7*
  • 1Department of Nephrology, Dongguan Tungwah Hospital, Dongguan, China
  • 2Dongguan Key Laboratory of Precise Prevention & Treatment of Chronic Kidney Disease and Complications, Dongguan, China
  • 3KingPoint Data Technology Co., Ltd., Guangzhou, China
  • 4Department of Cardiology, Dongguan SongShan Lake Tungwah Hospital, Dongguan, China
  • 5Department of Nephrology, Kiang Wu Hospital, Macau, Macao SAR, China
  • 6Department of Nephrology, The Third Affiliated Hospital of Sun Yat-sen University, Guangzhou, China
  • 7Department of Emergency Medicine, Dongguan Tungwah Hospital, Dongguan, China

Context: The conventional risk factors for cardiocerebral vascular events (CVCs) in non-Hemodialysis (HD) patients cannot be directly applied to HD patients due to the unique characteristics of this population. More accurate information on the risk of progression to CVCs is needed for clinical decisions.

Objective: To develop and validate time-dependent predictive models for the progression of CVCs in HD patients.

Design, setting, and participants: Development and validation of time-dependent predictive models using demographic, clinical, and laboratory data from 3 dialysis centers between 2017 and 2021. These models were developed using time-dependent Cox proportional hazards regression and assessed for discrimination using the concordance index, goodness of fit using the Akaike information criterion and net reclassification improvement.

Main outcome measures: CVCs included acute heart failure, acute hematencephalon, cardiac or brain-derived death, acute myocardial infarction, acute cerebral infarction, ischemic cardiomyopathy, unstable angina pectoris, and stable angina pectoris.

Results: The development and validation cohorts included 233 and 215 patients, respectively. The most accurate model included age, sex, hemoglobin, serum albumin, serum phosphate, white blood cell count, blood flow rate and ultrafiltration volume during HD (C index, 0.704; 95% CI, 0.639–0.768 in the development cohort and 0.775; 95% CI, 0.706–0.843 in the validation cohort). In the validation cohort, this model was more accurate than a model containing variables whose p value in the Cox proportional hazards regression was less than 0.05 (NRI: 0.351, 95% CI: −0.115–0.565).

Conclusion: A time-dependent model using routinely obtained laboratory tests can accurately predict progression to CVCs in HD patients.

Background

The global prevalence of end-stage renal disease (ESRD) has been steadily increasing (1, 2). The total population of individuals with ESRD receiving renal replacement therapy (RRT) surpassed 2.5 million in 2010 and is projected to exceed 5.439 million by 2030, with the majority of growth anticipated in Asia, where numbers are expected to increase from 0.968 million to approximately 2.162 million (3). The primary modality for RRT continues to be hemodialysis (HD), rather than peritoneal dialysis (PD) or kidney transplantation (4). The mortality rate in HD patients is high, primarily attributed to the occurrence of cardiocerebral vascular events (CVCs) (57). The clinical decision-making process for HD patients is challenging due to the diverse etiologies of kidney diseases, variable rates of disease progression, and competing risks of CVCs mortality. The accurate prediction of risk could enhance personalized decision-making, facilitating timely and appropriate patient care.

Currently, there is a lack of universally accepted predictive tools for the incidence of CVCs in HD patients. Therefore, physicians cannot identify patients at high risk of CVCs from HD patients in a timely manner or implement timely interventions to prevent or reduce the eventual occurrence of CVCs.

Recent studies have demonstrated that albuminuria (ALB) (8, 9) provides additional prognostic information regarding the progression to cardiovascular events (CVEs) in HD patients. One study investigated the incorporation of serum phosphate (P) and white blood cell (WBC) counts into prediction models (10). However, these models are either specific to CVEs but not CVCs or lack external validation. More accurate information regarding the risk of progression to CVCs in HD patients is required for clinical decisions about testing, treatment, and referral. The ideal model for predicting progression to CVCs in HD patients should have accuracy, be easy to implement, and demonstrate high generalizability across diverse patient populations in dependent cohorts.

Utilizing data extracted from the electronic medical records system (EMR) in three distinct HD centers, our study aimed to develop and externally validate an accurate but simple predictive model for the progression to CVCs in HD patients. The objective was also to utilize routinely measured variables in HD patients to develop a time-dependent predictive model for the progression of CVCs that could be easily implemented in clinical settings. We were particularly intrigued by models that exclusively rely on clinical laboratory data, facilitating the reporting of CVCs risk alongside laboratory test results.

Materials and methods

Study population

Development cohort

The development cohort was derived from the Kidney Disease Center at Dongguan Tungwah Hospital, Dongguan, Guangdong, China, and the Department of Nephrology at the Luogang Branch of the Third Affiliated Hospital of Sun Yat-sen University, Guangzhou, Guangdong, China, between 2017 and 2021. The patients underwent two to three weekly sessions of HD for a minimum duration of 3 months, employing bicarbonate-based dialysate and polysulfone membrane dialyzers. Low-molecular-weight heparin was administered as an anticoagulant during the dialysis procedure. The criteria for sample selection were as follows: patients who were at least 18 years old at the time of initiation of HD and who underwent quarterly biochemical tests and who had no history of prior renal transplantation or surgical procedures within the past year. Individuals with cachexia, malignant tumors, or other life-limiting conditions and patients who had experienced significant blood loss or severe infectious diseases during the study period were excluded from the study. The censored data were as follows: among the 56 screened patients who were excluded, reasons for exclusion included patient refusal to participate (n = 35), renal transplantation (n = 10), withdrawal of consent (n = 8), and poor adherence to HD treatment (n = 3).

Validation cohort

The validation cohort was derived from the Department of Nephrology at the Tianhe Branch of the Third Affiliated Hospital of Sun Yat-sen University, Guangzhou, Guangdong, China, between 2017 and 2021. The criteria for sample selection and exclusion were the same as those for the development cohort. The study was reviewed and approved by the institutional review boards at Dongguan Tungwah Hospital, Dongguan, Guangdong, China, and the Third Affiliated Hospital of Sun Yat-sen University, Guangzhou, Guangdong, China. Written informed consent was obtained from all patients or their immediate relatives, and the study was conducted in accordance with ethical guidelines (Ethical Batch No. 2020DHLL010).

Follow-up plan

The study investigated the occurrence of CVCs during a 5-year follow-up period. Upon occurrence of a cardiovascular event in a patient, we promptly initiated the consultation system and collaborated with the cardiologist to establish an accurate diagnosis of the cardiovascular event. In cases of cerebrovascular events, we sought consultation from a specialist in cerebrovascular medicine. We ensured that both diagnosing physicians held the title of attending physician or higher.

Variables

Candidate independent variables

The candidate independent variables were selected based on their face validity, encompassing demographic factors such as age and sex, as well as physical examination variables, including blood pressure measurements before, during, and after the HD program, along with weight assessments; HD-related variables, including dialysis age and blood flow rate (BFR) during HD; and laboratory variables obtained from three different hemodialysis centers reflecting inflammation, nutrition and renal complications such as anemia (Table 1). The aforementioned data were collected on a regular quarterly basis. The baseline values were defined as the mean value within a 3-month period following enrollment. Exceptionally, the baseline values of serum creatinine (Cr) and blood urea nitrogen (BUN) were selected as the first values before HD after enrollment. Variables with more than 30% missing values were not included in the analysis.

Table 1
www.frontiersin.org

Table 1. Baseline characteristics of the development and validation cohorts.

Dependent variable

The outcome of interest was CVCs, including acute heart failure, acute hematencephalon, cardiac- or brain-derived death, acute myocardial infarction, acute cerebral infarction, ischemic cardiomyopathy, unstable angina pectoris, and stable angina pectoris. Upon the occurrence of a cardiovascular event in a patient, we promptly initiated the consultation system and collaborated with the cardiologist to establish an accurate diagnosis. For patients with cerebrovascular events, we sought consultation from a specialist in cerebrovascular medicine. We ensured that both diagnosing physicians held the designation of attending physicians or higher. The prediction of the risk of CVCs occurrence has significant implications for decision-making processes undertaken by patients, physicians, and healthcare systems. The time horizons for risk prediction were 1, 3, and 5 years.

Statistical analysis

Model development

We developed a sequential series of models and compared them based on similar variable filtering conditions. We employed a combination of clinical guidelines and forward selection methodology to ascertain the selection of variables. Variables not significantly associated with CVCs (p > 0.10) were excluded from further analyses, except for Model 6, which was based on the chi-square test. Multiple testing adjustment: Bonferroni correction was applied to control type I error rate, with adjusted significance level at α = 0.002 (0.05/25 variables). The enhancement of model performance in time-dependent multivariate Cox proportional hazards regression models was evaluated by incorporating new candidate variables, with discrimination and goodness-of-fit metrics employed for assessment. Seven sequential models were developed with distinct variable selection criteria, Model 1: Variables with p < 0.005 in Cox regression; Model 2: p < 0.01; Model 3: p < 0.05; Model 4: p < 0.1; Model 5: All candidate variables; Model 6: Clinically meaningful variables reflecting anemia, nutrition, inflammation, and dialysis adequacy; Model 7: Variables with both clinical relevance and p < 0.1.

Time-varying covariates

Laboratory parameters (e.g., ALB, WBC) were updated at each hemodialysis session and analyzed as time-dependent covariates in extended Cox models.

Prediction model validation

By incorporating the concordance index (C index) and Akaike information criterion (AIC) and models 3, 4 and 6, the most accurate models among all the models were evaluated in the external validation dataset. The baseline hazard function and coefficients derived from the developed model were kept constant and applied to the validation dataset.

Prediction model performance

We employed a range of methodologies to assess the performance of the models across both the development and validation datasets.

Discrimination

Discrimination refers to the model’s ability to accurately differentiate between two classes of outcomes, namely, the occurrence and nonoccurrence of CVCs. The computation of the C index serves as a measure for assessing discrimination (1113).

Calibration

The process of calibration refers to the extent to which the predicted probabilities align numerically with the observed outcomes. We compared the observed and predicted risks of CVCs occurrence across all time points and quantified the magnitude of deviation between them (14).

Goodness of fit

The overall model fit for sequential models was compared using the AIC, which incorporates both the statistical goodness of fit and the complexity of the model by penalizing an increase in parameters required to achieve a certain level of fit (11).

Reclassification

Reclassification refers to the reassignment of patients from one class to another based on changes in their risk category assignment. The improvement in reclassification was quantified using the net reclassification improvement (NRI) statistic (13).

Sensitivity analysis

The patients were categorized into subgroups based on the occurrence of endpoint events (CVCs). Two models from Model 3, Model 4, and Model 6 were selected to predict the endpoint events in these subgroups. A sensitivity analysis was conducted to evaluate the predictive efficacy of the models by comparing the positive or negative prediction probability with the actual positive or negative rate (predictions falling within a certain range can be considered accurate).

All the statistical analyses were performed using SPSS 20.0. A two-sided p < 0.05 was considered to indicate statistical significance.

Results

Cohort description

The development and validation cohorts included 233 and 215 patients, respectively (Table 1). Patients in the development and validation cohorts were similar in age and sex. The HD-related variables pre- and post-HD blood pressure, BFR, and UV were similar between the two cohorts. Laboratory data, such as hemoglobin (HGB), WBC, ALB and so on, were similar between these two cohorts. The proportion of events was greater in the development cohort than in the validation cohort (74 patients with CVCs occurrence [31.7%] vs. 41 patients with CVCs occurrence [19.0%]) (Figure 1).

Figure 1
www.frontiersin.org

Figure 1. Kaplan–Meier survival curves for the development and validation cohorts. The proportion of events was greater in the development cohort than in the validation cohort (74 patients with CVCs occurrence [31.7%] vs. 41 patients with CVCs occurrence [19.0%]). CVCs, cardiocerebral vascular events. At risk: the number of patients with potential risks at a specific time point; Censored: the number of patients censored at a specific time point; Event: the cumulative count of risk events that transpired until a specific time point.

A chi-square test was performed on all patients to clarify the correlation between the occurrence of CVCs and variables, and the p values are shown in Table 1. p values less than 0.05 indicated a significant correlation between them. Indicators were chosen for the prediction model based on statistical correlation and clinical practice.

Prediction model performance in the development cohort

The hazard ratios for the variables and statistics regarding discrimination and goodness of fit for successive models in the development dataset are presented in Table 2. According to the chi-square test presented in Table 1, variables that showed no significant association with CVCs (p > 0.10) were excluded from further analyses, except for Model 6.

Table 2
www.frontiersin.org

Table 2. Hazard ratios and goodness of fit for sequential models in the development dataset.

Model 1, which included serum potassium (K), cholesterol (CHOL), hematocrit (HCT) and HD vascular access types (accesses), performed poorly (C index, 0.596; 95% confidence interval [CI], 0.523–0.668). The C index improved with the inclusion of serum calcium (Ca), BUN and blood glucose (GLU) in Model 2 (0.654; 95% CI, 0.585–0.722). The addition of more indicators from Model 1 to Model 5 led to a gradual increase in the C index value. Model 6 and Model 7, each containing representative clinical variables, which reflect the state of anemia, nutrition and inflammation, electrolyte levels and dialysis adequacy, showed similar performances (C index 0.704 for Model 6; C index 0.680 for Model 7). Despite a similar C index, the AIC was lower for Model 6 than for Model 7 (659.2 vs. 669.8, respectively). The inclusion of all variables in Model 7 had the largest C index (0.772; 95% CI, 0.709–0.834) among all models, and the AIC was lowest for Model 6 (AIC: 659.2) compared with the other models. Given these results, models 1, 2, and 7 were not considered in further evaluation steps.

Performance of the prediction model in the validation cohort

After screening by the C-index and AIC, Model 3, Model 4, and Model 6 were selected to determine their ability to predict CVCs occurrence in the validation cohort. Survival was calculated using time-dependent Cox regression and compared with the observed Kaplan–Meier (KM) curve of the validation cohort. Figure 2 shows the observed vs. predicted non-CVCs rates at all time points for models 3, 4, and 6 in the validation cohort. The mean square error (MSE) difference between the observed and predicted probabilities over the entire period of risk was lower for Model 6 than for Models 3 and 4 (0.0067, 0.0059, and 0.0052, respectively).

Figure 2
www.frontiersin.org

Figure 2. Observed vs. predicted non-CVCs rates using Models 3, 4, and 6 in the validation cohort. The MSE difference between the observed and predicted probabilities over the entire time-point of the CVCs was lower for Model 6 than for Models 3 and Model 4 (0.0067, 0.0059, and 0.0052, respectively). MSE, mean square error; CVCs, cardiocerebral vascular events.

The NRI risk for CVCs (Figure 3) in the validation cohort was analyzed. The 1200th daytime node was chosen for NRI analysis. Overall, model 6 outperformed models 3 and 4, with NRIs of 0.351 (95% CI, −0.115–0.565) and 0.329 (95% CI, −0.085–0.604), respectively.

Figure 3
www.frontiersin.org

Figure 3. Scatter plot of the NRIs for Model 3, Model 4, and Model 6. The 1200th day time node was chosen for NRI analysis. Model 3 vs. Model 4: NRI -0.007 (95% CI: −0.220, 0.246); Model 3 vs. Model 6: NRI 0.351 (95% CI: −0.115, 0.565); Model 4 vs. Model 6: NRI 0.329 (95% CI: −0.085, 0.604). NRI, net reclassification improvement.

Sensitivity analyses

Sensitivity analyses of each model were performed by comparing the difference in the prediction of positive and negative events for the validation cohort. The accuracy of predicting the occurrence of positive events, specifically CVCs, was represented by the NRI + (positive net reclassification improvement), as demonstrated in Table 3. Conversely, the NRI represents the accuracy of predicting the absence of CVCs.

Table 3
www.frontiersin.org

Table 3. Sensitivity analyses for Model 3, Model 4, and Model 6.

The predictive performance of Model 6 surpassed that of Models 3 and 4 in determining the occurrence of CVCs in HD patients (Model 3 vs. Model 6: 0.131, 95% CI, −0.084–0.387; Model 4 vs. Model 6: 0.155, 95% CI, −0.005–0.368). Additionally, model 6 outperformed models 3 and 4 in predicting the absence of CVCs in HD patients (model 3 vs. model 6: 0.220, 95% CI, −0.058–0.396; model 4 vs. model 6: 0.173, 95% CI, −0.137–0.341).

Correlation coefficients and HR values for model 6

Correlation Coefficients and HR Values for Model 6 are presented in Table 4. These data represent the coefficients of the time-dependent Cox regression model, where positive and negative signs indicate positive and negative correlations with the endpoint event, respectively. The HR values quantify the extent of influence on the endpoint event. The PH test indicated that all independent variables, except for gender, satisfy the assumption of proportional hazards. Therefore, the correlation coefficient of gender showed as a time-dependent covariate coefficient, while HR represents a dynamic value that varies over time.

Table 4
www.frontiersin.org

Table 4. Correlation coefficients and HR values for Model 6.

Discussion

In patients who initiate HD, the early mortality rate is high (15), with up to 40% of deaths attributed to cardiovascular disease (16, 17). Among patients undergoing maintenance HD, cardiovascular incidents continue to be the primary cause of mortality (18). Moreover, the incidence of cerebrovascular events is significantly greater in HD patients than in the general population (6, 7). Despite this knowledge, the risk factors for CVCs in HD patients have yet to be fully elucidated to establish a practical risk prediction model. However, insufficient research has investigated cardiovascular and cerebrovascular events in HD patients.

Risk prediction has garnered increasing attention in the past two decades, with the emerging literature suggesting enhanced patient outcomes through personalized risk prediction and advancements in information technology facilitating seamless integration of risk prediction models into electronic health records (EHRs) (1922). The availability of these risk prediction tools and their incorporation into clinical practice guidelines has resulted in enhanced adherence to treatment guidelines and facilitated individual decision-making (23, 24). Despite these advantages, the lack of readily applicable and externally validated models has impeded the widespread integration of risk prediction across all medical disciplines (25).

We developed and validated a set of time-dependent risk prediction models for progression to CVCs among patients undergoing HD. Our models utilize laboratory data routinely obtained from patients undergoing HD, which can be easily integrated into a laboratory information system or clinic EHRs.

Among patients undergoing HD, there is significant heterogeneity in the risk of progressing to CVCs. For instance, Table 5 illustrates the predicted probability of CVCs occurrence in two real patients with similar clinical conditions using models 3, 4, and 6. For patient A, with an observed time of 3.17 years for CVCs, models 3 and 4 exhibited similar predictive probabilities at the 3-year time point, while model 6 demonstrated an obviously elevated value. Model 6 provides substantially different risk predictions than models 3 and 4. In comparison to Model 3, Model 6 demonstrated a 28% increase in the predicted risk for Patient A at the 3-year time point and an 18% increase in the risk for Patient B at the 1-year time point.

Table 5
www.frontiersin.org

Table 5. Predicted probability of CVCs for 2 patients using our prediction modelsa.

Our models rely on demographic data and laboratory markers reflecting dialysis adequacy or complications associated with chronic kidney disease (CKD) to predict the future risk of CVCs occurrence in HD patients. Similar to previous findings from Xiaobing Liu (26) and Yukiko Matsubara (8), we observed that lower levels of HGB, ALB, UV, and BFR, as well as male sex, were associated with faster progression to CVCs. The time-varying HR for male gender (HR = exp.[0.419 + 0.413 × ln(t)]) suggests accumulating risk over dialysis vintage. This pattern may be attributed to the vascular calcification process and fluctuations in hormone levels, both of which warrant further research. Furthermore, higher WBC and serum P levels, along with older age, were also predictive of an increased risk of CVCs in HD patients. These markers may reflect dialysis adequacy or underlying processes of inflammation or malnutrition.

In contrast to previous studies that have individually associated these laboratory markers with the progression to CVCs, our research integrates them collectively into a unified risk equation (the risk calculator can be accessed in Supplementary material 1). In addition, we demonstrated no improvement in model performance with the addition of vascular access variables and physical examination (systolic blood pressure, diastolic blood pressure). Although these variables are clearly important for the diagnosis and management of HD patients, the lack of improvement in model performance may reflect their broad applicability.

Our study has several limitations. First, the inclusion of three independent dialysis centers in large general hospitals restricts the generalizability of our findings to patients undergoing HD in community hospitals or community healthcare centers. Second, it should be noted that the majority of participants in this study were of Han nationality, which reflects the ethnic distribution within Guangdong Province, China. However, it is important to acknowledge that certain underrepresented ethnic minorities may have a greater susceptibility to CVCs. Third, the sample size in this study was relatively limited, which may have influenced the stability of the model due to the low incidence of CVCs events. Future research will involve a multicenter study enrolling at least 2000 patients to validate these findings.

In summary, we have developed and validated time dependent, highly precise predictive models for progression to CVCs in patients undergoing HD. The optimal model we developed utilizes readily available laboratory data and demonstrates high accuracy in predicting the short-term risk of CVC occurrence. Moreover, its seamless integration into a laboratory information system ensures easy implementation. Model parameters can be integrated into hospital electronic health record systems to automatically compute risk scores and trigger alerts following each dialysis session. Prospective studies are currently being planned to assess the impact of the model on clinical decision-making, including modifications to dialysis regimens or enhanced cardiovascular monitoring protocols. Furthermore, future research will focus on optimizing the model’s performance through systematic exploration. External validation across multiple diverse cohorts of HD patients and evaluation in clinical trials are imperative.

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics statement

The studies involving humans were approved by the institutional review boards at Dongguan Tungwah Hospital, Dongguan, Guangdong, China, and the Third Affiliated Hospital of Sun Yat-sen University, Guangzhou, Guangdong, China. 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

HZ: Formal analysis, Funding acquisition, Writing – original draft, Writing – review & editing. MZ: Formal analysis, Methodology, Validation, Writing – review & editing. YX: Formal analysis, Methodology, Validation, Writing – review & editing. YQ: Investigation, Methodology, Software, Writing – review & editing. NX: Data curation, Validation, Writing – review & editing. YY: Data curation, Software, Writing – review & editing. HL: Methodology, Resources, Supervision, Writing – review & editing. HP: Funding acquisition, Methodology, Supervision, Writing – review & editing. XL: Methodology, Project administration, Resources, Supervision, Writing – review & editing. XS: Funding acquisition, Methodology, Project administration, Resources, Supervision, Writing – review & editing. SL: Methodology, Project administration, Resources, Supervision, Writing – review & editing.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This work was in part supported by the grants from Basic and Applied Basic Research Foundation of Guangdong Province, China (2020B1515120037), Dongguan Science and Technology of Social Development Program (2018507150461632, 20221800903262), and The Science and Technology Development Fund, Macau SAR (0013/2023/RIC).

Acknowledgments

Statistical services were provided in part by KingPoint Data Technology Co., Ltd., Guangzhou, China. We acknowledge the valuable contributions of this study.

Conflict of interest

YQ was employed by KingPoint Data Technology Co., Ltd.

The remaining 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.

The reviewer JH declared a shared affiliation with the author XL to the handling editor at the time of review.

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.1481866/full#supplementary-material

References

1. Foster, BJ. Survival improvements for Europeans with ESKD. Kidney Int. (2020) 98:834–6. doi: 10.1016/j.kint.2020.05.032

PubMed Abstract | Crossref Full Text | Google Scholar

2. McCullough, KP, Morgenstern, H, Saran, R, Herman, WH, and Robinson, BM. Projecting ESRD incidence and prevalence in the United States through 2030. J Am Soc Nephrol. (2019) 30:127–35. doi: 10.1681/ASN.2018050531

PubMed Abstract | Crossref Full Text | Google Scholar

3. Liyanage, T, Ninomiya, T, Jha, V, Neal, B, Patrice, HM, Okpechi, I, et al. Worldwide access to treatment for end-stage kidney disease: a systematic review. Lancet. (2015) 385:1975–82. doi: 10.1016/S0140-6736(14)61601-9

PubMed Abstract | Crossref Full Text | Google Scholar

4. Bello, AK, Levin, A, Tonelli, M, Okpechi, IG, Feehally, J, Harris, D, et al. Assessment of global kidney health care status. JAMA. (2017) 317:1864–81. doi: 10.1001/jama.2017.4046

PubMed Abstract | Crossref Full Text | Google Scholar

5. Drueke, TB, and Floege, J. Cardiovascular complications of chronic kidney disease: pioneering studies. Kidney Int. (2020) 98:522–6. doi: 10.1016/j.kint.2020.07.001

PubMed Abstract | Crossref Full Text | Google Scholar

6. Iseki, K, and Fukiyama, KOkawa Dialysis Study (OKIDS) Group. Clinical demographics and long-term prognosis after stroke in patients on chronic haemodialysis. The Okinawa Dialysis Study (OKIDS) Group. Nephrol Dial Transplant. (2000) 15:1808–13. doi: 10.1093/ndt/15.11.1808

PubMed Abstract | Crossref Full Text | Google Scholar

7. Belmar, L, de Francisco, AL, Bueno, L, Piñera, C, Monfá, E, Kislikova, M, et al. Strokes in patients on haemodialysis: incidence, onset time and associated factors. Nefrologia. (2014) 34:347–52. doi: 10.3265/Nefrologia.pre2014.Jan.12298

PubMed Abstract | Crossref Full Text | Google Scholar

8. Matsubara, Y, Kimachi, M, Fukuma, S, Onishi, Y, and Fukuhara, S. Development of a new risk model for predicting cardiovascular events among hemodialysis patients: population-based hemodialysis patients from the Japan Dialysis outcome and practice patterns study (J-DOPPS). PLoS One. (2017) 12:e0173468. doi: 10.1371/journal.pone.0173468

PubMed Abstract | Crossref Full Text | Google Scholar

9. Tang, W, Zhang, Y, Wang, Z, Yuan, X, Chen, X, Yang, X, et al. Development and validation of a multivariate model for predicting heart failure hospitalization and mortality in patients receiving maintenance hemodialysis. Ren Fail. (2023) 45:2255686. doi: 10.1080/0886022X.2023.2255686

PubMed Abstract | Crossref Full Text | Google Scholar

10. You, X, Gu, B, Chen, T, Li, X, Xie, G, Sang, C, et al. Development of long-term cardiovascular disease risk prediction model for hemodialysis patients with end-stage renal disease based on nomogram. Ann Palliat Med. (2021) 10:3142–53. doi: 10.21037/apm-21-286

PubMed Abstract | Crossref Full Text | Google Scholar

11. Toll, DB, Janssen, KJ, Vergouwe, Y, and Moons, KG. Validation, updating and impact of clinical prediction rules: a review. J Clin Epidemiol. (2008) 61:1085–94. doi: 10.1016/j.jclinepi.2008.04.008

PubMed Abstract | Crossref Full Text | Google Scholar

12. Steyerberg, EW, Vickers, AJ, Cook, NR, Gerds, T, Gonen, M, Obuchowski, N, et al. Assessing the performance of prediction models: a framework for traditional and novel measures. Epidemiology. (2010) 21:128–38. doi: 10.1097/EDE.0b013e3181c30fb2

PubMed Abstract | Crossref Full Text | Google Scholar

13. Pencina, MJ, D'Agostino, RB Sr, D'Agostino, RB Jr, et al. Evaluating the added predictive ability of a new marker: from area under the ROC curve to reclassification and beyond. Stat Med. (2008) 27:157–72. doi: 10.1002/sim.2929

Crossref Full Text | Google Scholar

14. D’Agostino, RB, and Nam, BH. Evaluation of the performance of survival analysis models: discrimination and calibration measures. San Francisco, CA: Elsevier; (2004):1–25.

Google Scholar

15. Eckardt, KU, Gillespie, IA, Kronenberg, F, Richards, S, Stenvinkel, P, Anker, SD, et al. High cardiovascular event rates occur within the first weeks of starting hemodialysis. Kidney Int. (2015) 88:1117–25. doi: 10.1038/ki.2015.117

PubMed Abstract | Crossref Full Text | Google Scholar

16. Hole, BD, Casula, A, and Caskey, FJ. Quality assuring early dialysis care: evaluating rates of death and recovery within 90 days of first dialysis using the UK renal registry. Clin Kidney J. (2022) 15:1612–21. doi: 10.1093/ckj/sfab238

PubMed Abstract | Crossref Full Text | Google Scholar

17. de Jager, DJ, Grootendorst, DC, Jager, KJ, van Dijk, P, Tomas, LM, Ansell, D, et al. Cardiovascular and noncardiovascular mortality among patients starting dialysis. JAMA. (2009) 302:1782–9. doi: 10.1001/jama.2009.1488

PubMed Abstract | Crossref Full Text | Google Scholar

18. Hoppe, K, Schwermer, K, Olewicz-Gawlik, A, Klysz, P, Kawka, A, Baum, E, et al. Dialysis vintage and cardiovascular injury as factors influencing long-term survival in peritoneal dialysis and hemodialysis. Adv Clin Exp Med. (2017) 26:251–8. doi: 10.17219/acem/64936

PubMed Abstract | Crossref Full Text | Google Scholar

19. D'Agostino, RB Sr, Grundy, S, Sullivan, LM, and Wilson, P. Validation of the Framingham coronary heart disease prediction scores: results of a multiple ethnic groups investigation. JAMA. (2001) 286:180–7. doi: 10.1001/jama.286.2.180

PubMed Abstract | Crossref Full Text | Google Scholar

20. Tian, Y, Zhang, Y, He, J, Chen, L, Hao, P, Li, T, et al. Predictive model of acute kidney injury after spontaneous intracerebral hemorrhage: a multicenter retrospective study. Eur Stroke J. (2023) 8:747–55. doi: 10.1177/23969873231184667

PubMed Abstract | Crossref Full Text | Google Scholar

21. Demirjian, S, Bashour, CA, Shaw, A, Schold, JD, Simon, J, Anthony, D, et al. Predictive accuracy of a perioperative laboratory test-based prediction model for moderate to severe acute kidney injury after cardiac surgery. JAMA. (2022) 327:956–64. doi: 10.1001/jama.2022.1751

PubMed Abstract | Crossref Full Text | Google Scholar

22. Harasemiw, O, Nayak, JG, Grubic, N, Ferguson, TW, Sood, MM, and Tangri, N. A predictive model for kidney failure after nephrectomy for localized kidney Cancer: the kidney Cancer risk equation. Am J Kidney Dis. (2023) 82:656–65. doi: 10.1053/j.ajkd.2023.06.002

PubMed Abstract | Crossref Full Text | Google Scholar

23. Leslie, WD, Morin, S, and Lix, LM. A before-and-after study of fracture risk reporting and osteoporosis treatment initiation. Ann Intern Med. (2010) 153:580–6. doi: 10.7326/0003-4819-153-9-201011020-00007

PubMed Abstract | Crossref Full Text | Google Scholar

24. Pozen, MW, D'Agostino, RB, Selker, HP, Sytkowski, PA, and Hood, WB Jr. A predictive instrument to improve coronary-care-unit admission practices in acute ischemic heart disease. A prospective multicenter clinical trial. N Engl J Med. (1984) 310:1273–8. doi: 10.1056/NEJM198405173102001

PubMed Abstract | Crossref Full Text | Google Scholar

25. Muller-Riemenschneider, F, Holmberg, C, Rieckmann, N, et al. Barriers to routine risk-score use for healthy primary care patients: survey and qualitative study. Arch Intern Med. (2010) 170:719–24. doi: 10.1001/archinternmed.2010.66

PubMed Abstract | Crossref Full Text | Google Scholar

26. Liu, X, Yan, C, Niu, X, and Zeng, J. Establishment of a nomogram model for predicting cardiovascular and cerebrovascular events in diabetic nephropathy patients receiving maintenance hemodialysis. Appl Bionics Biomech. (2022) 2022:2909726. doi: 10.1155/2022/2909726

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: time-dependent, predictive model, cardiocerebral vascular events, hemodialysis, laboratory tests

Citation: Zhong H, Zhang M, Xie Y, Qin Y, Xie N, Ye Y, Li H, Peng H, Liu X, Su X and Li S (2025) A time-dependent predictive model for cardiocerebral vascular events in chronic hemodialysis patients: insights from a prospective study. Front. Med. 12:1481866. doi: 10.3389/fmed.2025.1481866

Received: 19 August 2024; Accepted: 22 May 2025;
Published: 04 June 2025.

Edited by:

Hua Zhu, The Ohio State University, United States

Reviewed by:

Jiang He, Department of Cardiology, Sun Yat-sen University, China
Suli Huang, Shenzhen Center for Disease Control and Prevention, China
Xin-wang Zhu, The First Affiliated Hospital of China Medical University, China

Copyright © 2025 Zhong, Zhang, Xie, Qin, Xie, Ye, Li, Peng, Liu, Su and Li. 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: Xiaoyan Su, c3V4aWFveWFuNzY5QGhvdG1haWwuY29t; Shaohong Li, bGlzaGFvaG9uZzkyM0Bob3RtYWlsLmNvbQ==

These authors have contributed equally to this work

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.