Data source
Data originated from the 2005 Longitudinal National Health Insurance Database (LHID2005) provided by the Health and Welfare Statistics Application Center, Ministry of Health and Welfare, Taiwan. The NHI Program was launched in Taiwan in March 1995 and since 2010, has enrolled more than 99% of the inhabitants of Taiwan.
The LHID2005 contains all the original claims data for a randomly selected sample of one million beneficiaries enrolled in the NHI program in 2005. Analysis of the NHI data using the multistage stratified systematic sampling method has established that there are no statistically significant differences regarding sex or age between these one million insured individuals and the general population. The LHID2005 contains all NHI enrolment files, claims data and the registry for prescription drugs, providing comprehensive information on the beneficiaries from 1996 to 2013. Reimbursed TCM services, including Chinese herbal products, acupuncture/moxibustion, and Chinese traumatology therapy in ambulatory clinics, are also included in the database.
Ethical approval
To protect the privacy of patients, Taiwan’s National Health Research Institute (NHRI) has encrypted the names of patients, health care providers and medical institutions, to permit public access for research purposes [17]. Thus, the Review Board waived the requirement for obtaining informed consent from the patients. This study was conducted in accordance with the Declaration of Helsinki and was approved by the Institutional Review Board of China Medical University, Taichung, Taiwan (CMUH-104-REC2-115-R4).
Study population and variables
To investigate the safety and efficacy of the concurrent use of CHM and anticoagulants, we identified two cohorts: one using both CHMs and anticoagulants (CHM cohort) and a cohort using only anticoagulants (non-CHM cohort) as shown in Fig. 1.
The inclusion criteria for the CHM cohort specified that patients were prescribed anticoagulants consisting of warfarin, apixaban, dabigatran, or rivaroxaban for at least one continuous week and were concomitantly using CHM for at least one continuous week with an overlap of at least 1 day in prescribing/dispensing dates and their supplied days between January 2000 and December 2013. Patients diagnosed with preexisting coagulopathy (ICD-9-CM 286, 287), were aged < 20 years, or who withdrew from the NHIRD during the follow-up period were excluded in both cohorts (Fig. 1).
The non-CHM cohort included patients using anticoagulants consisting of warfarin, apixaban, dabigatran, or rivaroxaban for at least one continuous week, without any CHM prescriptions.
Patients in the non-CHM cohort were matched 1:1 by propensity score with the CHM cohort according to age, gender, index year of baseline comorbidities and medication use at baseline. The baseline comorbidities were adjusted according to the CHA2DS2–VASc score, which is used to assess stroke and bleeding risks in patients with atrial fibrillation [18]. We defined a patient with a history of comorbidities as having been diagnosed at least once in the Emergency Department (ED) or Inpatient Department before the index date for diabetes (ICD-9-CM 250), hypertension (ICD-9-CM 401, 405), chronic kidney disease (ICD-9-CM 580–589, A350), previous ischemic stroke/transient ischemic attack (TIA)/thromboembolism (TE) (ICD-9-CM 433, 434, 435, 437, 438, 452, 453), or myocardial infarction (MI)/peripheral arterial disease (PAD)/aortic plaque (ICD-9-CM 410, 411, 412, 413, 414, 440, 444, 445). Anticoagulants used were warfarin, apixaban, dabigatran and rivaroxaban. Other baseline medications included aspirin and antiplatelet agents, as shown in Supplementary Table 1.
The primary outcome included hemorrhagic stroke (ICD-9-CM 430, 431, 432), GI bleeding (ICD-9-CM 530, 531, 532, 534, 533, 534, 535, 537, 562, 569, 578, 784.8, 786) and urogenital bleeding (ICD-9-CM 596, 599, 602, 608, 620, 621, 623, 625, 626, 627, 640, 641, 666, 674), as well as nose and eye bleeding (ICD-9-CM 360, 362, 364, 372, 374, 376, 379, 386, 388, 784) confirmed by diagnosis once in the ED or during hospitalization within the study period. For all study participants, the end date of follow-up was whichever came first: (1) from the first day of prescribing/dispensing dates of CHMs or anticoagulants, up until 8 weeks after the last prescribing/dispensing dates of CHMs or anticoagulants; (2) withdrawal from the database; or (3) the occurrence of a major bleeding event (Fig. 1).
Statistical analysis
After matching both cohorts with the propensity score to reduce confounding bias, the Chi-square test was used to assess differences between the CHM and non-CHM cohorts in categorical variables, represented by numbers and percentages. Cox proportional hazard regression associated with hazard ratios (HRs), adjusted hazard ratios (aHRs) and 95% confidence intervals (CIs) was used to assess between-cohort risks for major bleeding events. Incidence rates per 10,000 person-years were estimated for both cohorts. The Kaplan-Meier method was used to estimate cumulative incidence rates of bleeding. The log-rank test was applied to compare the difference between two bleeding incidence curves. We also analyzed CHMs previously reported as possibly involved in herb-drug interactions with anticoagulants and we identified the top 10 single herbs and formulas most frequently used in the CHM cohort, to investigate potential associations with major bleeding events. Cox proportional hazard regression analysis was used to study further associations of major bleeding events with particular herbs or formulas, and calculated the adjusted hazard ratios (aHRs) and 95% confidence intervals (CIs) for the cumulative risk of survival.
Data analysis was performed using SAS version 9.4 (SAS Institute Inc., Cary, NC, USA). Values for the cumulative incidence curve were plotted by R software (R Foundation for Statistical Computing, version 3.6.1, Vienna, Austria http://www.r-project.org). A p-value of less than 0.05 in a two-sided test was considered statistically significant.