- 1Department of Nephrology, The First People’s Hospital of Qinzhou, The Tenth Affiliated Hospital of Guangxi Medical University, Qinzhou, Guangxi, China
- 2Department of Laboratory Medicine, The First People’s Hospital of Qinzhou, The Tenth Affiliated Hospital of Guangxi Medical University, Qinzhou, Guangxi, China
- 3Department of Medical Statistics, The First People’s Hospital of Qinzhou, The Tenth Affiliated Hospital of Guangxi Medical University, Qinzhou, Guangxi, China
Introduction: Antithrombotic management in Chronic Kidney Disease (CKD) is a clinical dilemma. This study aimed to empirically evaluate the “de facto interchangeability” of the antiplatelet indobufen and the anticoagulant rivaroxaban by comparing their real-world effectiveness and safety in hospitalized CKD patients.
Methods: In this retrospective cohort study (2020-2024), we analyzed CKD patients treated with indobufen or rivaroxaban. A multi-stage analysis first used machine learning to assess baseline cohort comparability, overcoming limitations of p-value-based tests. Subsequently, a Linear Mixed Model (LMM), adjusted for confounders including polypharmacy, assessed independent drug effects on in-hospital thrombosis, hemorrhage, and longitudinal laboratory markers.
Results: Machine learning demonstrated the clinical comparability of the indobufen and rivaroxaban cohorts. The incidence of in-hospital thrombosis was numerically lower in the indobufen group (3.65% vs. 7.58%; P = 0.101), while hemorrhage rates were similar (2.19% vs. 2.27%; P = 1). The LMM analysis, beyond verifying indobfen’s expected antiplatelet activity (modulating MPV, PDW), revealed pleiotropic effects (increased prealbumin, HDL-C) and a significant reduction in urine occult blood (P < 0.001), suggesting renal safety. Notably, the model demonstrated that apparent effects on hemoglobin and eGFR were attributable to confounding by co-medications, not a direct drug effect.
Conclusion: In this real-world CKD cohort, indobufen and rivaroxaban demonstrated comparable clinical effectiveness and safety. Combining machine learning with longitudinal models helps to statistically adjust for complex confounders like polypharmacy, thereby providing a more robust estimate of a drug’s independent effect.
Introduction
Chronic Kidney Disease (CKD) has emerged as a formidable global public health threat, characterized by high morbidity and a rapidly escalating mortality rate. According to the Global Burden of Disease Study 2023, CKD was responsible for 1.52 million deaths worldwide in 2023, a staggering 90.9% increase since 2000 (GBD, 2023 Causes of Death Collaborators, 2025). Within this challenging context, antithrombotic management is particularly problematic, limited by a narrow therapeutic window and potential nephrotoxicity (Saeed et al., 2024; Cases et al., 2021; Addisu et al., 2025). Indeed, CKD represents a paradoxical state with a concurrently elevated risk of both thrombosis and bleeding (Wang et al., 2023; Montomoli et al., 2024). On one hand, CKD constitutes a prethrombotic state driven by key mechanisms, including: endothelial dysfunction from uremia and inflammation; elevated levels of procoagulant factors (e.g., fibrinogen); and platelet hyperreactivity with increased release of procoagulant microparticles (MPs) (Alhousani et al., 2021; Noumegni et al., 2022). This sterile inflammation, mediated by the activation of the NLRP3 inflammasome by damage-associated molecular patterns (DAMPs), is a key driver of podocyte injury and renal fibrosis, accelerating CKD progression (Datta et al., 2025b; Datta et al., 2025a). On the other hand, platelet dysfunction caused by uremic toxins, coupled with the reduced metabolic clearance of anticoagulant drugs in patients with renal insufficiency, also significantly increases the risk of bleeding complications (Tham et al., 2024; Li et al., 2025). This dual risk makes antithrombotic treatment decisions in patients with CKD exceptionally challenging. As traditional injectable anticoagulants are ill-suited for long-term outpatient management (Almajdi et al., 2023; Expert Consensus Working Group, 2023) there is a pressing need to optimize the risk-benefit balance of oral antithrombotic agents in this complex population (Cases et al., 2021; Lip et al., 2021; Lau et al., 2022).
This study aims to conduct a head-to-head comparison of indobufen and rivaroxaban. The rationale for this comparison originates not from traditional pharmacological theory but from our observation of clinical equipoise. In patients with CKD, traditional antithrombotic agents have significant limitations: aspirin has failed to definitively reduce ischemic risk while increasing bleeding events (Wolfe et al., 2025; Stefanini et al., 2021; Festa et al., 2023), and the clinical applicability of warfarin is restricted by its narrow therapeutic window, a higher bleeding risk compared to Direct Oral Anticoagulants (DOACs), and the need for frequent INR monitoring (Cohen et al., 2022; Lau et al., 2022). Consequently, clinical focus has shifted to alternatives like clopidogrel (Gorog et al., 2023) and the FXa inhibitor rivaroxaban, whose indications have expanded (Ten Cate et al., 2021; Gao and Rahman, 2022). Indobufen, a reversible cyclooxygenase-1 (COX-1) inhibitor, is widely used in East Asia for aspirin-intolerant patients. Unlike aspirin’s irreversible inhibition, indobufen’s effect is transient, with platelet function recovering within 24 h. This theoretically confers better gastrointestinal tolerability and a lower bleeding risk (Dai C. et al., 2024; Jiang et al., 2025).
While studies have compared indobufen with aspirin and clopidogrel (Liu et al., 2024; Cavalcante et al., 2025), meta-analyses indicate its efficacy in preventing ischemic events is comparable to aspirin’s in cardiovascular settings (e.g., post-PCI/CABG), but with a significantly lower bleeding risk (Dai C. et al., 2024; Luo et al., 2024; Ren et al., 2024). This safety advantage may stem from its reversible, non-competitive inhibition of COX-1, which helps preserve gastric mucosal function while maintaining its antiplatelet effect (Liu et al., 2023). This favorable risk-benefit profile makes indobufen a promising option for patients with nephrotic syndrome (NS), who exhibit high thrombotic risk and a high bleeding propensity, partially attributed to gastrointestinal mucosal edema. However, existing research has focused on cardiovascular and stroke populations, largely neglecting its use in CKD. This study therefore aims to compare the real-world efficacy and safety of these functionally equivalent drugs, providing immediate clinical evidence and laying the foundation for future prospective trials and mechanistic explorations.
Materials and methods
Study design and data source
This study employed a retrospective cohort design, utilizing real-world data sourced from our institution’s Hospital Information System (HIS) and Laboratory Information System (LIS). The dataset included all patients hospitalized for Chronic Kidney Disease (CKD) between January 2020 and December 2024. Given the retrospective nature of the study and its focus on patients with specific clinical diagnoses, a separate cohort of healthy subjects was not enrolled to serve as a negative control group. Instead, all comparative analyses were conducted between different subgroups within the defined patient populations.
Inclusion criteria were as follows: (1) patients hospitalized for CKD who received antithrombotic therapy during their stay; and (2) patients who met at least one of the following conditions: ① presence of a hypercoagulable state during hospitalization or at discharge; ② occurrence of a new thrombotic event after the initiation of antithrombotic therapy; ③ occurrence of a bleeding event during hospitalization; or ④ presence of a bleeding tendency at discharge.
Exclusion criteria included: (1) age <12 years, to avoid confounding from developmental differences in pathophysiology and drug metabolism specific to pediatric populations; (2) total length of hospital stay <72 h, to ensure a sufficient observation period for the assessment of treatment effects and clinical outcomes; (3) significant missing data in key baseline characteristics, laboratory results, or medication records, as this would preclude a reliable assessment and introduce significant bias; and (4) receipt of combination antithrombotic therapy (i.e., concurrent use of two or more anticoagulant or antiplatelet agents) during hospitalization, in order to isolate the effects of the specific monotherapy under investigation and avoid confounding. For subsequent longitudinal analyses (e.g., paired t-tests, LMM), a complete-case approach was used, including only subjects with data available at both pre- and post-treatment time points, as reflected by the varying sample sizes (n) in the results tables.
Data extraction and validation framework
To ensure high data fidelity, we developed a custom hybrid framework for data extraction from electronic medical records (EMRs). This system combined a Large Language Model (LLM) for identifying complex clinical events (e.g., thrombosis, bleeding) with a dual-pipeline validation process for structured data. All extracted events were temporally sequenced to distinguish pre-treatment from post-treatment onset and were manually adjudicated by two clinical researchers. A comprehensive description of the framework is provided in Supplementary Text S1.
Criteria for event adjudication and status assessment
The diagnostic criteria for thrombotic and bleeding events were based on established international guidelines (e.g., ISTH, CTCAE v5.0). We also defined criteria for assessing patient status, specifically “Hypercoagulable State” (Chinese expert consensus on the application of anticoagulant technology in renal replacement therapy for critical patients (2023), 2023) and “Bleeding Tendency” (Jiang et al., 2023). The detailed operational definitions for all events and the specific laboratory criteria for these states are provided in Supplementary Text S1.
Extraction and processing of laboratory results
To standardize the analysis of longitudinal laboratory data, we extracted results from predefined temporal windows corresponding to admission, discharge, 1/3/6-month follow-ups, and periods immediately pre- and post-treatment. For each window, a single representative data point was selected based on a consistent rule. The detailed definitions of these time windows and the selection criteria are provided in Supplementary Text S1.
Analysis of clinical patterns and determinants of oral antithrombotic drug selection
This study employed a multi-stage analytical strategy to first identify clinical patterns associated with a hypercoagulable state and then to elucidate the factors driving the selection of specific oral antithrombotic agents. First, we utilized association rule mining with the Apriori algorithm (implemented in the R package “arules”) to identify clinical features and medication combinations that co-occurred with a hypercoagulable state (Si et al., 2022). Next, the analysis focused on the selection of specific oral antithrombotic agents (rivaroxaban, indobufen, and clopidogrel). We initially screened for clinical features associated with the choice of these drugs by performing a differential analysis of baseline characteristics, with an p-value <0.05 considered significant. Subsequently, LASSO (Least Absolute Shrinkage and Selection Operator) regression was employed to perform feature reduction and refine the selection of key predictors (Luo et al., 2023).
Assessment of efficacy, adverse events, and Ancillary Effects
(1) Assessment of Treatment Efficacy: The incidence of new-onset thrombotic events during the treatment period or within specified post-treatment observation periods, and the presence of a hypercoagulable state during the discharge phase; (2) Assessment of Adverse Events: The incidence of new-onset bleeding events and the presence of a bleeding tendency. Both outcomes were evaluated during the treatment period and within specified post-treatment observation periods; (3) Analysis of Ancillary Effects on Laboratory Parameters:First, paired t-tests were used to compare the pre- and post-treatment values for each laboratory indicator. Subsequently, a Linear Mixed Model (LMM), implemented using the R package “lme4”, was constructed to provide a more robust assessment of the drug’s impact over time (Jimba et al., 2024; Yuan et al., 2025).
Statistical analyses
All analyses were conducted using R (v4.4.0) and Python (v3.13.0). Continuous variables were expressed as mean ± SD or median (IQR) and compared using ANOVA or Kruskal–Wallis tests as appropriate. Categorical variables were presented as n (%) and compared using Chi-squared or Fisher’s exact tests. P-values for multiple comparisons were adjusted using the Benjamini–Hochberg FDR method. A two-sided P-value <0.05 was considered statistically significant.
Result
Patient baseline characteristics and description
Following the procedure outlined in Figure 1A (Step 1), all EMRs for these patients were consolidated from the hospital’s information system, regardless of the reason for or department of admission. As shown in Figure 2A, the most common clinical features across all hospitalization events included Nephrotic Syndrome (n = 2,474 events), hyperlipidemia (n = 2,370), edema (n = 2,291), chronic kidney disease at admission (eGFR <90 mL/min/1.73 m2; n = 2,206), and hypertension (n = 2,157). The majority of patients were admitted under the Department of Nephrology (75.9%, 847/1095), as depicted in Figure 2B. Other contributing departments included Rheumatology and Immunology (4.93%), Endocrinology (4.12%), Emergency Medicine (3.23%), Pediatrics (2.87%), and various other clinical departments (8.96%).
Figure 1. Flowchart of the study design and analytical workflow. The process begins with the definition of the patient cohort and the extraction of multi-source data using methods including Large Language Models (LLMs). Subsequently, a machine learning pipeline (including Apriori, LASSO, PCA, and decision trees) is used to refine the analysis to a head-to-head comparison of Indobufen and Rivaroxaban. Finally, a Linear Mixed Model (LMM) is applied to assess the independent effects of the target medication on clinical outcomes and laboratory parameters.
Figure 2. Baseline clinical landscape and medication profile of the study cohort. (A) Alluvial plot illustrating the prevalence of key clinical features, with counts for “Yes,” “No,” and “Unknown” statuses indicated across all 4,511 hospitalization events; (B) Pie chart depicting patient distribution by clinical department, predominantly Nephrology (75.9%); (C) Sunburst diagram illustrating the hierarchical classification of medications, focusing on the top 50 most frequently administered agents across all hospitalization episodes. Inner rings represent major therapeutic categories, and outer rings detail more specific drug subclasses or types.
To understand the therapeutic landscape, the top 50 most frequently administered medications were hierarchically classified (Figure 2C). The most frequently recorded drug classes were cardiovascular agents (e.g., antihypertensives, lipid-lowering agents, beta-blockers), anti-infectives (e.g., beta-lactams), and hematologic agents (e.g., anticoagulants and antiplatelets). Additionally, drugs for the digestive system (e.g., acid suppressants, prokinetics), endocrine and metabolic agents (e.g., antidiabetic drugs), hormones and immunomodulators (e.g., immunosuppressants, corticosteroids), and diuretics also constituted a significant proportion of the administered medications.
Clinical features and initial medication patterns associated with hypercoagulable state
The temporal dynamics of the hypercoagulable state were assessed on a per-event basis, as illustrated in Figure 3A. The proportion of events classified as a hypercoagulable state was 38.15% (n = 1,190/3,119) at admission; this proportion decreased to 11.14% (n = 210/1,185) at discharge. Follow-up assessments showed the proportion to be 10.45% (n = 120/1,158) at 2 weeks, 12.37% (n = 199/1,609) at 1 month, 10.24% (n = 170/1,660) at 3 months, and 10.56% (n = 197/1,866) at 6 months post-discharge.
Figure 3. Hypercoagulable State Dynamics, Associated Factors, and Medication Use (A) Alluvial plot depicting the temporal evolution of hypercoagulable status from admission through 6 months post-discharge; (B) Apriori association network illustrating co-occurrence relationships between common clinical characteristics (red nodes) and frequently prescribed medications (blue nodes); (C) Yearly (2020-2024) total usage and number of unique patients for all medications identified as associated with hypercoagulable states; (D) Yearly (2020-2024) patient count distribution for selected anticoagulant and antiplatelet medications; (E) Alluvial plot illustrating co-prescription patterns in the focused cohort (n = 2,265 hospitalization events). Each axis represents a medication, segmented by use (Yes) or non-use (No). Flows indicate the number of events sharing specific medication combinations.
To identify clinical factors and medications robustly associated with a hypercoagulable state, an association rule analysis using the Apriori algorithm was performed on the dataset. The results, centered on the hypercoagulable state, are presented as a network in Figure 3B. This network illustrates the associations between clinical factors and medications frequently co-occurring with or used to manage this state, including various anticoagulants and antiplatelets. Subsequently, the annual total usage and the number of unique patients were tallied for medications identified as associated with the hypercoagulable state by the Apriori analysis (Figure 3C). For instance, furosemide and atorvastatin calcium, while not direct antithrombotic agents, were frequently co-prescribed in patients with hypercoagulable risk profiles and demonstrated high overall usage and patient counts.
Further analysis focused on the 5-year usage trends of specific anticoagulant and antiplatelet agents identified by the Apriori algorithm (Figure 3D). Notably, the use of indobufen and rivaroxaban increased rapidly following their clinical introduction in 2022. The use of clopidogrel showed a declining trend.
To transition from this broad exploration to a focused analysis of oral antithrombotic drug selection patterns, a specific cohort was defined for all subsequent analyses. Based on the strict inclusion and exclusion criteria detailed in the Methods section (focusing on patients who received antithrombotic therapy for a clear clinical indication). Figure 3E illustrates the intricate co-prescription patterns among core medications within this focused cohort, and Table 1 details the corresponding distributions of total dosage and duration of use for these medications.
Analysis of selection patterns for oral antithrombotic agents
Differential analysis of 136 baseline variables identified 62 with significant differences between the drug groups (FDR P < 0.05). Beyond identifying these baseline differences, this comprehensive variable set also served to systematically screen for potential off-target adverse events. To this end, a subset of these indicators reflecting key organ function (e.g., eGFR, creatinine, hemoglobin) was tracked longitudinally at admission, discharge, and multiple post-discharge time points (2 weeks, 1, 3, and 6 months). This approach was based on the rationale that clinically significant adverse reactions would likely be captured during follow-up, either through detectable laboratory abnormalities or via their documentation in medical records. A complete comparison of all variables is provided in Supplementary Table S1. LASSO regression was then employed to refine this set to 17 key predictive variables (Table 2), which formed the foundation for the primary analysis. Following multiple imputation of missing values, Principal Component Analysis (PCA) was applied to these 17 variables to visualize the overall distribution of pre-treatment clinical profiles for each drug group (Figure 4A). The analysis revealed a substantial overlap in the multidimensional space between patients treated with indobufen and those treated with rivaroxaban, with their respective clusters located in close proximity. In contrast, patients receiving clopidogrel formed a distinct cluster, clearly separated from the other two groups.
Figure 4. Baseline Characteristics and Oral Antithrombotic Selection. (A) PCA plot of 17 selected baseline clinical variables for patients receiving Clopidogrel (green), Indobufen (red), or Rivaroxaban (blue); (B) Confusion matrix illustrating the performance of an initial three-class decision tree model on both training and test sets; (C) Optimized binary clinical decision tree model distinguishing Clopidogrel use from combined Indobufen/Rivaroxaban use. Model performance: F1 score = 0.794, AUC = 0.811; (D) Distribution of key decision tree predictor variables across the three original drug groups; (E) Baseline pre-treatment levels of eGFR, PT, APTT, FIB, D-Dimer, and PLT across the three drug groups.
To further quantify and identify the clinical feature combinations that distinguish the use of these three oral agents, a decision tree model was developed based on the 17 selected variables. A preliminary three-class model designed to differentiate among all three groups was evaluated. The confusion matrices for the training and test sets (Figure 4B) demonstrated that the model had very limited ability to distinguish between indobufen and rivaroxaban, a finding consistent with the high degree of similarity observed in the PCA. Given the high degree of similarity in baseline features and the initial model’s performance, the indobufen and rivaroxaban groups were combined into a single “Indobufen/Rivaroxaban” group. A binary decision tree analysis was then conducted to distinguish this combined group from the “Clopidogrel” group. The final model demonstrated strong discriminatory performance, achieving an F1 score of 0.794 and an AUC of 0.811. The structure of the final decision tree is presented in Figure 4C, highlighting the key decision nodes: Hypocalcemia at Admission, Urine Albumin-to-Creatinine Ratio (UACR) at Admission, NS, Uremia, and Coronary Heart Disease. To visually illustrate the differences in these key predictors, Figure 4D displays their distributions across the original three drug groups, along with the results of pairwise statistical comparisons. Finally, Figure 4E presents the baseline levels of a key renal function indicator (eGFR) and coagulation-related parameters (PT, APTT, FIB, D-Dimer, PLT) across the three groups. Notably, despite rivaroxaban being an anticoagulant and indobufen an antiplatelet agent, there were no clinically significant differences in their pre-treatment PLT or coagulation function indicators.
Comparative efficacy and safety of rivaroxaban and indobufen
Having established baseline comparability, we directly compared the longitudinal efficacy and safety of the two drugs. Figure 5A tracks the evolution of hypercoagulable states and the incidence of in-hospital thrombotic events. At admission, both treatment groups exhibited a similarly high prevalence of hypercoagulability: 75.00% (99/132 events) in the rivaroxaban group and 78.35% (322/411 events) in the indobufen group. The incidence of new-onset thrombotic events during hospitalization was 7.58% (10/132 events) for the rivaroxaban group and 3.65% (15/411 events) for the indobufen group (Chi-squared test, P = 0.101). Among the subgroup of events with pre-existing hypercoagulability at admission, the incidence of in-hospital thrombosis was 3.03% (3/99 events) for rivaroxaban and 2.80% (9/322 events) for indobufen (Fisher’s exact test, P = 1). Following in-hospital treatment, the prevalence of hypercoagulability decreased significantly in both groups at discharge. This resolved state was largely sustained in both cohorts throughout the 6 month follow-up period.
Figure 5. Longitudinal Trajectories of Thrombotic and Bleeding Events in Rivaroxaban and Indobufen Cohorts. Alluvial plots showing the temporal evolution of (A) hypercoagulability/thrombosis and (B) bleeding tendency/hemorrhage in patients treated with Rivaroxaban (top panels) versus Indobufen (bottom panels). The plots track patient status from admission to 6 month follow-up.
Complementing the efficacy analysis, Figure 5B presents a safety comparison by illustrating the longitudinal evolution of bleeding tendency and the incidence of in-hospital hemorrhage. At admission, both cohorts exhibited a similarly low prevalence of bleeding tendency: 4.55% (6/132 events) for the rivaroxaban group and 5.60% (23/411 events) for the indobufen group (Chi-squared test, P = 0.879), indicating comparable baseline safety profiles. The incidence of in-hospital hemorrhage events was also nearly identical and remained low in both groups: The incidence of in-hospital hemorrhage events was also low and statistically indistinguishable between the groups: 2.27% (3/132 events) for rivaroxaban and 2.19% (9/411 events) for indobufen (Fisher’s exact test, P = 1).
Table 3 details the distribution of thrombotic events, hemorrhage events, hypercoagulability at discharge, and bleeding tendency at discharge across different dosing regimens of rivaroxaban and indobufen. Notably, the incidence of hypercoagulability at discharge in the indobufen 10 mg bid group (13.27%, 28/211) was significantly higher than in the 20 mg qd group (3.05%, 4/131). Compared to the 10 mg bid regimen, the 20 mg qd regimen was associated with a lower likelihood of persistent hypercoagulability at discharge (OR = 0.19; 95% CI, 0.04 - 0.64; P = 0.004).
It is noteworthy that of the 12 total in-hospital hemorrhage events, 11 (91.67%) occurred in patients assessed as having ‘no’ bleeding tendency at admission. This suggests that the currently used criteria for assessing bleeding tendency may be insufficient for predicting clinical hemorrhage, and future work could focus on developing more suitable predictive models based on actual bleeding outcomes.
Analysis of effects of target oral antithrombotic agents on laboratory
Before assessing treatment effects, we first confirmed that the treatment exposure was comparable between the two drug cohorts. Table 1 demonstrates similar durations of use for indobufen and rivaroxaban, which provides the rationale for a direct pre-post comparison. Therefore, paired t-tests were initially used to quantify the changes in various laboratory parameters from pre-to post-treatment. Due to the large number of parameters with significant changes, Table 4 focuses on presenting a selection of the most clinically relevant and statistically significant results. A complete analysis of all tested variables can be found in Supplementary Table S2.
Recognizing that aggregate statistics such as mean differences can mask inter-individual heterogeneity, we used paired plots to visualize the trajectories of individual patient parameters. In the rivaroxaban group (Figure 6A), a decreasing trend was observed for parameters including MPV (72.7%), FIB (75.0%), FDP (77.8%), D-Dimer (66.7%), and APTT (65.0%). Conversely, PDW (59.0%) and Antithrombin III (88.9%) showed an increasing trend. In the indobufen group (Figure 6B), MPV (73.6%) and FIB (67.7%) showed a decreasing trend, while PDW (71.5%), APTT (63.9%), and Antithrombin III (68.4%) showed an increasing trend.
Figure 6. Individual patient responses to Rivaroxaban and Indobufen. Paired plots show changes in key platelet and coagulation parameters from pre-treatment (red) to post-treatment (blue). Panels are stratified by drug ((A): Rivaroxaban; (B) Indobufen) and further subdivided by the direction of change (decrease vs. increase). Percentages in the panel titles indicate the proportion of patients in each subgroup.
Although paired t-tests revealed significant improvements in the lipid profile (e.g., LDL-C, HDL-C), these non-hemostatic changes were excluded from detailed visualization. These effects are more plausibly attributed to confounding by widespread concomitant statin therapy (Figure 3) rather than a direct pharmacological effect of the study drugs, for which there is limited mechanistic evidence. To further control for confounding factors and quantify the independent effects of the drugs, we developed a LMM. This model included age, sex and major co-medications as fixed effects, with individual hospitalization events treated as a random effect. The model results (Table 5; Supplementary Table S3) provide an assessment of the main drug effect (baseline differences between groups) and the interaction effect (the drug’s independent influence on the parameter’s trajectory).
Table 5. Independent effects of indobufen on laboratory parameters identified by linear mixed-effects model.
Using MPV as an example, the model showed no significant difference in baseline levels between the indobufen group and the non-user group (main drug effect: P = 0.335), indicating good comparability between the groups. However, the model revealed a significant negative interaction effect (Estimate = −0.32, P = 0.0003), indicating that after controlling for multiple confounders, indobufen use was independently associated with a more substantial decrease in MPV. These independent pharmacological associations included significant decreases in MPV and plateletcrit, and significant increases in PDW, prealbumin, CRP, and HDL-C (all interaction effects P < 0.05).
In contrast, for LDL-C, the model showed that the indobufen group had significantly higher baseline levels (main drug effect: P = 0.0097), suggesting a baseline imbalance. Although a significant interaction effect was observed (P = 0.0005), its attribution is complex due to the baseline differences and widespread co-administration of lipid-lowering drugs. The attribution of this effect is complex, potentially stemming from three non-mutually exclusive factors: (1) a form of confounding by severity, where patients with markedly higher baseline levels likely received more aggressive, unmodeled lipid-lowering therapy, leading to a greater observed reduction; (2) a synergistic interaction with concomitant medications such as atorvastatin; or (3) a genuine, independent lipid-modulating effect of indobufen itself. Differentiating between these possibilities is beyond the scope of this retrospective analysis and awaits confirmation from future prospective and mechanistic studies.
Discussion
The premise of our study is not a pharmacological hypothesis but an empirical observation of “de facto therapeutic interchangeability.” This observation stemmed from our initial finding that indobufen and rivaroxaban have become the primary emerging oral antithrombotic choices in recent years. Subsequent analysis then revealed that clinicians are applying these two drugs to a highly similar patient phenotype, particularly those with NS, while reserving agents like clopidogrel for distinct indications. This distinction is critical, as a direct comparison with clopidogrel would have been invalidated by severe confounding by indication, creating a methodological paradox where a seemingly appropriate peer-drug comparison would actually involve two fundamentally different patient populations with insufficient overlap in their baseline characteristics. In contrast, the observed interchangeability between indobufen and rivaroxaban presented a unique opportunity to address a zone of genuine clinical equipoise. Our study was therefore designed not to compare mechanisms in a vacuum, but to resolve a real-world clinical dilemma: in a population where clinicians already treat these two drugs as viable alternatives, which strategy offers a better net clinical benefit? This spontaneous adoption of functionally equivalent options, driven by real-world complexity, is a common feature of clinical practice when managing frail, comorbid, or high-risk patients (Campitelli et al., 2021; Mohamed et al., 2021; Gorog et al., 2023), often requiring individualized adjustments beyond standard guidelines.
Regarding the primary endpoints (Figure 5), indobufen showed a numerically lower rate of thrombosis compared to rivaroxaban (3.65% vs. 7.58%). While not statistically significant, this trend suggests a potential thromboprotective advantage for indobufen that warrants investigation in larger RCTs. Crucially, this potential benefit did not come at the cost of safety, as in-hospital bleeding rates were low and comparable between the groups (2.19% vs. 2.27%). This balanced profile is of major clinical importance for CKD patients, who must navigate high risks of both thrombosis and bleeding. As for secondary outcomes, both treatments effectively mitigated the laboratory-defined “hypercoagulable state” through 6 months of follow-up and performed similarly well in the highest-risk subgroup of patients with pre-existing hypercoagulability. It is crucial to acknowledge the significant heterogeneity of the “hypercoagulable state” as a diagnostic label. Its underlying mechanisms are disease-specific, and its diagnostic criteria vary widely across conditions without a universal gold standard (Johansson et al., 2022; Marchetti et al., 2022; Song et al., 2022; Lim et al., 2023). Moreover, composite biomarkers often fail to capture the true coagulation status (Scheiner et al., 2022; Lim et al., 2023). To address this ambiguity, our study bypassed the composite label and instead analyzed the core laboratory parameters independently. Since these individual markers are standardized and widely applied, this approach provides objective, granular data to support individualized assessments that can be adapted to the diverse diagnostic criteria used in practice.
Our analysis of hemostatic markers revealed a dual pattern: convergent effects reflecting overall clinical improvement and distinct pharmacological fingerprints unique to each drug. The convergent effects included mitigation of the procoagulant environment in both groups, evidenced by decreased FIB and MPV, and increased antithrombin III. The fingerprints themselves were perfectly aligned with each drug’s mechanism: Rivaroxaban’s pharmacological fingerprint was characterized by a marked suppression of D-Dimer, indicating reduced systemic coagulation, whereas Indobufen’s pharmacological fingerprint was more pronounced in its impact on platelet-related indices like plateletcrit. However, despite these clear group-level profiles, the significant inter-individual variability shown in Figure 6 highlights the limitations of simple pre-post comparisons.
Moreover, in clinical practice, patients often receive combination therapy, creating a significant risk of confounding bias where observed changes cannot be attributed to indobufen alone. To isolate the effects of indobufen from the confounding bias of combination therapy, we employed a LMM, adjusting for age, sex, and key concomitant medications. This analysis provided more robust evidence for indobufen’s known antiplatelet activity (Liu et al., 2023; Dai C. et al., 2024; Dai W.-B. et al., 2024), Reinforcing its independent role in decreasing MPV and plateletcrit while increasing PDW (Table 5, Part A). More interestingly, the LMM uncovered potential pleiotropic effects, linking indobufen independently to increases in prealbumin and HDL-C. This finding is supported by recent research suggesting indobufen can modulate the MAPK inflammatory pathway (Bai et al., 2024). However, this pleiotropy has clear boundaries, it does not extend to broad antioxidant effects, as indobufen does not alter F2-isoprostane generation driven by oxidative stress (Davì et al., 1999). Finally, while the model also linked indobufen to elevated thyroid hormones (FT4, T3), we interpret this not as a direct drug effect but as a likely confounding signal. Such fluctuations are characteristic of Non-Thyroidal Illness Syndrome (NTIS) in CKD patients and probably reflect overall clinical improvement rather than a specific pharmacological action (Yan et al., 2017; Mahashabde et al., 2024).
A key advantage of the LMM approach is its ability to disentangle initial findings from the effects of confounding. This is best illustrated by the hemoglobin result. Initially, the significant rise in HGB in the indobufen group was a particularly noteworthy finding, supporting the hypothesis that its reversible COX-1 inhibition confers superior gastrointestinal safety (Dai W.-B. et al., 2024; Ren et al., 2024). However, this effect became non-significant in the LMM after adjusting for co-medications. This highly informative “negative” result indicates the improvement was likely due to co-prescribed PPIs or mucosal protectants, not indobufen itself. Likewise, the concerning decline in eGFR seen in the initial analysis, which raised fears of nephrotoxicity (Hanaoka et al., 2025), also became non-significant in the LMM. This suggests the effect was likely confounded by co-administered glucocorticoids rather than being a direct drug-induced toxicity. Ultimately, the LMM served a crucial dual purpose: it isolated indobufen’s true independent effects while preventing the misattribution of effects from other drugs—be they protective (PPIs) or potentially harmful (steroids). This highlights a fundamental principle for real-world data analysis: traditional comparisons are inadequate for handling the complex, time-varying confounders inherent in clinical practice. Advanced longitudinal models are therefore essential to disentangle these factors, convert messy clinical reality into clear insights, and accurately isolate the independent effect of a drug. Our findings caution clinicians against prematurely attributing changes in hemoglobin or eGFR to indobufen itself, highlighting the need to consider the confounding effects of co-prescribed PPIs and steroids.
Finally, this study has several limitations that inform future research. Our analysis framework, designed to analyze real-world EHR data, cannot directly measure pharmacokinetic (PK) or pharmacodynamic (PD) parameters. Its utility lies in generating hypotheses from large-scale clinical observations. For instance, the analysis was structured to compare drug-level effects rather than dose-specific outcomes—a methodological choice to avoid the model overfitting and loss of statistical power that would result from fragmenting the cohort into smaller subgroups. This approach revealed that patients receiving indobufen “100 mg, bid” showed a significantly greater resolution of hypercoagulability than those on “200 mg, qd,” despite an identical total daily dose (Table 3). This finding suggests that dosing frequency, and by extension the underlying PK/PD profile, may influence clinical efficacy. It raises a question that our RWD platform is suited to identify but not directly answer: why the same daily dose yields different outcomes. Answering this requires a dedicated mechanistic study integrating our findings with direct PK/PD measurements (e.g., drug concentrations via LC-MS and platelet function via thromboelastography). Therefore, this study serves not as a final conclusion, but as a foundational step to guide future, targeted mechanistic research.
Data availability statement
The raw data supporting the conclusions of this article contain sensitive patient information and are not publicly available due to privacy regulations and restrictions imposed by the institutional ethics committee of The First People’s Hospital of Qinzhou. However, a de-identified and interactive visualization of the dataset is accessible through a web application, which allows for the exploration of data distributions and key clinical characteristics presented in this study. Further inquiries or requests for access to aggregated, non-identifiable data can be directed to the corresponding author.
Ethics statement
The studies involving humans were approved by The Medical Ethics Committee of The First People’s Hospital of Qinzhou (Approval No. 2023046). The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation was not required from the participants or the participants’ legal guardians/next of kin in accordance with the national legislation and institutional requirements.
Author contributions
LZ: Funding acquisition, Investigation, Methodology, Writing – original draft. ML: Resources, Validation, Writing – original draft. TH: Data curation, Resources, Writing – original draft. HZ: Investigation, Supervision, Writing – original draft. CH: Investigation, Writing – original draft. ZP: Investigation, Writing – original draft. ZC: Investigation, Writing – original draft. JN: Conceptualization, Investigation, Project administration, Validation, Writing – review and editing. JT: Formal Analysis, Software, Visualization, Writing – review and editing.
Funding
The authors declare that financial support was received for the research and/or publication of this article. Qinzhou City Science and Technology Research and Development Plan Project [20232531] awarded to LZ.
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. To assist with language refinement, improve clarity and readability, and format text to meet journal guidelines. It was also used to help draft ancillary submission materials, including the scope statement and author contributions section. The authors have meticulously reviewed and edited all AI-generated text and take full responsibility for the final content of the 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/fphar.2025.1694163/full#supplementary-material
Abbreviations
APTT, Activated Partial Thromboplastin Time; CKD, Chronic Kidney Disease; DOAC, Direct Oral Anticoagulant; eGFR, estimated Glomerular Filtration Rate; EMRs, electronic medical records; FIB, Fibrinogen; HIS, Hospital Information System; LASSO, Least Absolute Shrinkage and Selection Operator; LIS, Laboratory Information System; LLM, Large Language Model; LMM, Linear Mixed Model; LMWH, Low-Molecular-Weight Heparin; PCA, Principal Component Analysis; PLT, Platelet Count; PT, Prothrombin Time; PK, Pharmacokinetic; PD, pharmacodynamic.
References
Addisu, Z. D., Demsie, D. G., Tafere, C., Yazie, T. S., Endeshaw, D., Tefera, B. B., et al. (2025). Time in the therapeutic range, bleeding event, and their determinants in older patients with atrial fibrillation on warfarin in Ethiopia: multicenter cross-sectional study. Front. Pharmacol. 16, 1541592. doi:10.3389/fphar.2025.1541592
Alhousani, M., Malik, S. U., Abu-Hashyeh, A., Poznanski, N. J., Al-Hasan, S., Roth, D. F., et al. (2021). Using oral anticoagulants among chronic kidney disease patients to prevent recurrent venous thromboembolism: a systematic review and meta-analysis. Thromb. Res. 198, 103–114. doi:10.1016/j.thromres.2020.11.036
Almajdi, A., Almutairi, S., and Alharbi, M. (2023). Safety and efficacy of apixaban versus low-molecular weight heparin or vitamin-K antagonists for venous thromboembolism treatment in patients with severe renal failure: a systematic review and meta-analysis. Thromb. Res. 229, 77–85. doi:10.1016/j.thromres.2023.06.027
Bai, L., Qu, W., Cheng, X., Yang, H., Huang, Y.-P., Wang, Z., et al. (2024). Multispecies transcriptomics identifies SIKE as a MAPK repressor that prevents NASH progression. Sci. Transl. Med. 16, eade7347. doi:10.1126/scitranslmed.ade7347
Campitelli, M. A., Bronskill, S. E., Huang, A., Maclagan, L. C., Atzema, C. L., Hogan, D. B., et al. (2021). Trends in anticoagulant use at nursing home admission and variation by frailty and chronic kidney disease among older adults with atrial fibrillation. Drugs Aging 38, 611–623. doi:10.1007/s40266-021-00859-1
Cases, A., Gomez, P., Broseta, J. J., Perez Bernat, E., Arjona Barrionuevo, J. de D., Portolés, J. M., et al. (2021). Non-valvular atrial fibrillation in CKD: role of vitamin K antagonists and direct oral anticoagulants. A narrative review. Front. Med. 8, 654620. doi:10.3389/fmed.2021.654620
Cavalcante, D. V. S., Krishna, M. M., Joseph, M., de Farias Santos, A. C. F., Mendes, B. X., Asbeg, N., et al. (2025). Indobufen versus aspirin in patients with indication for antiplatelet therapy: a systematic review and meta-analysis. Vasc. Pharmacol. 158, 107465. doi:10.1016/j.vph.2025.107465
Cohen, A. T., Sah, J., Dhamane, A. D., Lee, T., Rosenblatt, L., Hlavacek, P., et al. (2022). Effectiveness and safety of apixaban versus warfarin in venous thromboembolism patients with chronic kidney disease. Thromb. Haemost. 122, 926–938. doi:10.1055/s-0041-1740254
Dai, C., Liu, M., Yang, Z., Li, Y., Zhou, Y., Lu, D., et al. (2024). Real-world performance of indobufen versus aspirin after percutaneous coronary intervention: insights from the ASPIRATION registry. BMC Med. 22, 148. doi:10.1186/s12916-024-03374-3
Dai, W.-B., Ren, J.-Y., Hu, S.-T., Zhang, Y.-K., Gu, T.-S., Wu, X., et al. (2024). The safety and efficacy of indobufen or aspirin combined with clopidogrel in patients with acute myocardial infarction after percutaneous coronary intervention. Platelets 35, 2364748. doi:10.1080/09537104.2024.2364748
Datta, S., Rahman, M. A., Koka, S., and Boini, K. M. (2025a). High mobility group box 1 (HMGB1) mediates nicotine-induced podocyte injury. Front. Pharmacol. 15, 1540639. doi:10.3389/fphar.2024.1540639
Datta, S., Rahman, M. A., Koka, S., and Boini, K. M. (2025b). Mitigation of nicotine-induced podocyte injury through inhibition of thioredoxin interacting protein. Biomed. Pharmacother. 187, 118110. doi:10.1016/j.biopha.2025.118110
Davì, G., Ciabattoni, G., Consoli, A., Mezzetti, A., Falco, A., Santarone, S., et al. (1999). In vivo formation of 8-iso-prostaglandin f2alpha and platelet activation in diabetes mellitus: effects of improved metabolic control and vitamin E supplementation. Circulation 99, 224–229. doi:10.1161/01.cir.99.2.224
Expert Consensus Working Group (2023). Chinese expert consensus on the application of anticoagulant technology in renal replacement therapy for critical patients (2023). Chin. J. Nephrol., 155–164. doi:10.3760/cma.j.cn441217-20220809-00815
Festa, M. C., Rasasingam, S., Sharma, A., and Mavrakanas, T. A. (2023). Early discontinuation of aspirin among patients with chronic kidney disease undergoing percutaneous coronary intervention with a drug-eluting stent: a meta-analysis. Kidney360 4, e1245–e1256. doi:10.34067/KID.0000000000000223
Gao, F., and Rahman, F. (2022). DOACs and atherosclerotic cardiovascular disease management: can we find the right balance between efficacy and harm. Curr. Atheroscler. Rep. 24, 457–469. doi:10.1007/s11883-022-01022-w
GBD 2023 Causes of Death Collaborators (2025). Global burden of 292 causes of death in 204 countries and territories and 660 subnational locations, 1990-2023: a systematic analysis for the Global Burden of Disease Study 2023. Lancet, 406 (10513), 1811–1872. doi:10.1016/S0140-6736(25)01917-8
Gorog, D. A., Ferreiro, J. L., Ahrens, I., Ako, J., Geisler, T., Halvorsen, S., et al. (2023). De-escalation or abbreviation of dual antiplatelet therapy in acute coronary syndromes and percutaneous coronary intervention: a consensus statement from an international expert panel on coronary thrombosis. Nat. Rev. Cardiol. 20, 830–844. doi:10.1038/s41569-023-00901-2
Hanaoka, H., Aoki, T., Kosaka, T., Yoshinaga, S., Shibata, A., Sakai, R., et al. (2025). Chronic kidney disease and inflammatory cytokines in rheumatoid arthritis: a potential pathogenic link. Immunol. Med. 48, 161–170. doi:10.1080/25785826.2025.2460267
Jiang, S., Li, Y., Li, C., Shu, Z., Wu, S., Lu, X., et al. (2023). The association of coagulation indicators with in-hospital acute kidney injury and malignant events of patients with acute aortic dissection: a retrospective cohort study. J. Thorac. Dis. 15, 658–667. doi:10.21037/jtd-23-38
Jiang, Z., Zhu, D., and Meng, J. (2025). A randomized controlled trial comparing the efficacy and safety of indobufen versus aspirin in reducing target vessel restenosis after drug-eluting balloon angioplasty in patients with coronary artery disease. Clin. Sao Paulo Braz. 80, 100716. doi:10.1016/j.clinsp.2025.100716
Jimba, T., Kaneko, H., Azegami, T., Suzuki, Y., Okada, A., Ko, T., et al. (2024). Body weight change associated kidney outcomes of sodium-glucose cotransporter new users. Diabetes Obes. Metab. 26, 4535–4543. doi:10.1111/dom.15808
Johansson, M., Yan, H., Welinder, C., Végvári, Á., Hamrefors, V., Bäck, M., et al. (2022). Plasma proteomic profiling in postural orthostatic tachycardia syndrome (POTS) reveals new disease pathways. Sci. Rep. 12, 20051. doi:10.1038/s41598-022-24729-x
Lau, W. C. Y., Torre, C. O., Man, K. K. C., Stewart, H. M., Seager, S., Van Zandt, M., et al. (2022). Comparative effectiveness and safety between apixaban, dabigatran, edoxaban, and rivaroxaban among patients with atrial fibrillation: a multinational population-based cohort study. Ann. Intern. Med. 175, 1515–1524. doi:10.7326/M22-0511
Li, Y., Xi, L., Wang, D., Fan, G., Li, X., Shi, Y., et al. (2025). Risk of bleeding in pulmonary embolism patients concomitant with COVID-19 undergoing extended anticoagulation: a multicenter cohort study. Thromb. Res. 246, 109237. doi:10.1016/j.thromres.2024.109237
Lim, H. Y., Lui, B., Tacey, M., Barit, D., Patel, S. K., Donnan, G., et al. (2023). Global coagulation assays in patients with chronic kidney disease and their role in predicting thrombotic risk. Thromb. Res. 226, 127–135. doi:10.1016/j.thromres.2023.04.016
Lip, G. Y. H., Keshishian, A. V., Zhang, Y., Kang, A., Dhamane, A. D., Luo, X., et al. (2021). Oral anticoagulants for nonvalvular atrial fibrillation in patients with high risk of gastrointestinal bleeding. JAMA Netw. Open 4, e2120064. doi:10.1001/jamanetworkopen.2021.20064
Liu, J., Sun, P., and Qi, X. (2023). Reversible and non-competitive inhibition of cyclooxygenase by indobufen for efficient antiplatelet action and relief of gastrointestinal irritation. Pharmaceutics 15, 2135. doi:10.3390/pharmaceutics15082135
Liu, X., Lv, X., Peng, Y., Wang, J., Lei, J., Tang, C., et al. (2024). Clopidogrel with indobufen or aspirin in minor ischemic stroke or high-risk transient ischemic attack: a randomized controlled clinical study. BMC Neurol. 24, 81. doi:10.1186/s12883-024-03585-4
Luo, X., Sun, J., Pan, H., Zhou, D., Huang, P., Tang, J., et al. (2023). Establishment and health management application of a prediction model for high-risk complication combination of type 2 diabetes mellitus based on data mining. PloS One 18, e0289749. doi:10.1371/journal.pone.0289749
Luo, X., Lai, C., and Chen, T. (2024). The efficacy and safety of indobufen in patients with ischemic cardiovascular or cerebrovascular diseases: systematic review and meta-analysis. Front. Cardiovasc. Med. 11, 1509010. doi:10.3389/fcvm.2024.1509010
Mahashabde, M. L., Kumar, L., Bhimani, Y. R., Reddy, S. K., Nitendra Saketh, B. V., and Gharge, S. S. (2024). A critical investigation of sick euthyroid syndrome in chronic heart failure patients: addressing the need for accurate thyroid assessment. Cureus 16, e65985. doi:10.7759/cureus.65985
Marchetti, M., Russo, L., Giaccherini, C., Gamba, S., and Falanga, A. (2022). Hemostatic system activation in breast cancer: searching for new biomarkers for cancer risk prediction and outcomes. Thromb. Res. 213 (Suppl. 1), S46–S50. doi:10.1016/j.thromres.2022.01.019
Mohamed, M. O., Rashid, M., Timmis, A., Clarke, S., Lawson, C., Michos, E. D., et al. (2021). Sex differences in distribution, management and outcomes of combined ischemic-bleeding risk following acute coronary syndrome. Int. J. Cardiol. 329, 16–22. doi:10.1016/j.ijcard.2020.12.063
Montomoli, M., Candía, B. G., Barrios, A. A., and Bernat, E. P. (2024). Anticoagulation in chronic kidney disease. Drugs 84, 1199–1218. doi:10.1007/s40265-024-02077-6
Noumegni, S. R., Hoffmann, C., Tromeur, C., Didier, R., Bressollette, L., Lacut, K., et al. (2022). Risk factors of arterial events in patients with venous thromboembolism: a systematic review and meta-analysis. Thromb. Haemost. 122, 590–599. doi:10.1055/s-0041-1732300
Ren, Y., Zhu, Y., Yan, Q., Jin, H., and Luo, H. (2024). Multicenter retrospective cohort study demonstrates superior safety profile of indobufen over aspirin for Post-CABG antiplatelet therapy. Front. Pharmacol. 15, 1474150. doi:10.3389/fphar.2024.1474150
Saeed, Z., Sirolli, V., Bonomini, M., Gallina, S., and Renda, G. (2024). Hallmarks for thrombotic and hemorrhagic risks in chronic kidney disease patients. Int. J. Mol. Sci. 25, 8705. doi:10.3390/ijms25168705
Scheiner, B., Balcar, L., Nussbaumer, R. J., Weinzierl, J., Paternostro, R., Simbrunner, B., et al. (2022). Factor VIII/Protein C ratio independently predicts liver-related events but does not indicate a hypercoagulable state in ACLD. J. Hepatol. 76, 1090–1099. doi:10.1016/j.jhep.2021.12.038
Si, M., Li, D., Liu, T., Cai, Y., Yang, J., Jiang, L., et al. (2022). Triglycerides as biomarker for predicting systemic lupus erythematosus related kidney injury of negative proteinuria. Biomolecules 12, 945. doi:10.3390/biom12070945
Song, B., Chi, H., Peng, G., Song, Y., Cui, Z., Zhu, Y., et al. (2022). Characterization of coagulation-related gene signature to predict prognosis and tumor immune microenvironment in skin cutaneous melanoma. Front. Oncol. 12, 975255. doi:10.3389/fonc.2022.975255
Stefanini, G. G., Briguori, C., Cao, D., Baber, U., Sartori, S., Zhang, Z., et al. (2021). Ticagrelor monotherapy in patients with chronic kidney disease undergoing percutaneous coronary intervention: TWILIGHT-CKD. Eur. Heart J. 42, 4683–4693. doi:10.1093/eurheartj/ehab533
Ten Cate, H., Guzik, T. J., Eikelboom, J., and Spronk, H. M. H. (2021). Pleiotropic actions of factor Xa inhibition in cardiovascular prevention: mechanistic insights and implications for anti-thrombotic treatment. Cardiovasc. Res. 117, 2030–2044. doi:10.1093/cvr/cvaa263
Tham, D., Zhao, L., Yu, W., Kherani, J., Kou, R., Li, A., et al. (2024). Safety and efficacy of direct oral anticoagulants in chronic kidney disease: a meta-analysis. Res. Pract. Thromb. Haemost. 8, 102584. doi:10.1016/j.rpth.2024.102584
Wang, T.-F., Grubic, N., Carrier, M., Canney, M., Delluc, A., Hundemer, G. L., et al. (2023). Risk of venous thromboembolism or hemorrhage among individuals with chronic kidney disease on prophylactic anticoagulant after hip or knee arthroplasty. Am. J. Hematol. 98, 1374–1382. doi:10.1002/ajh.26994
Wolfe, R., Broder, J. C., Zhou, Z., Murray, A. M., Ryan, J., Chan, A. T., et al. (2025). Aspirin, cardiovascular events, and major bleeding in older adults: extended follow-up of the ASPREE trial. Eur. Heart J. 46, 4410–4422. doi:10.1093/eurheartj/ehaf514
Yan, W., Wang, L., Huang, T., and Xu, G. (2017). Treatment for non-thyroidal illness syndrome in advanced chronic kidney disease: a single-blind controlled study. J. Nephrol. 30, 557–565. doi:10.1007/s40620-016-0341-2
Yuan, G., Liao, Z., Liang, P., Cai, L., Zhou, K., Yin, T., et al. (2025). Noninvasive grading of renal interstitial fibrosis and prediction of annual renal function loss in chronic kidney disease: the optimal solution of seven MR diffusion models. Ren. Fail. 47, 2480751. doi:10.1080/0886022X.2025.2480751
Keywords: chronic kidney disease (CKD), indobufen, rivaroxaban, real-world study, large language model (LLM), linear mixed model (LMM)
Citation: Zhang L, Liu M, Huang T, Zhang H, Huang C, Pan Z, Chen Z, Ning J and Tang J (2025) Comparative effectiveness and pharmacological fingerprints of indobufen versus rivaroxaban in patients with chronic kidney disease: a single-center, real-world study. Front. Pharmacol. 16:1694163. doi: 10.3389/fphar.2025.1694163
Received: 04 September 2025; Accepted: 06 November 2025;
Published: 19 November 2025.
Edited by:
Narayan Prasad, Sanjay Gandhi Post Graduate Institute of Medical Sciences (SGPGI), IndiaReviewed by:
Sayantap Datta, Johns Hopkins University, United StatesOlimkhon Sharapov, Republican Specialized Scientific Practical Medical Center of Nephrology and Kidney transplantation, Uzbekistan
Copyright © 2025 Zhang, Liu, Huang, Zhang, Huang, Pan, Chen, Ning and Tang. 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: Jiameng Tang, amlhbWVuZ3RhbmdAMTI2LmNvbQ==
Lijun Zhang1