Objectives This study aims to assess whether digoxin has a different effect on mortality risk for women than it does for men in patients with heart failure (HF).
Design This study uses the UK-based The Health Information Network population database in a cohort study of the impact of digoxin exposure on mortality for men and women who carry the diagnosis of HF. Digoxin exposure was assessed based on prescribing data. Multivariable Cox proportional hazards models were used to assess whether there was an interaction between sex and digoxin affecting mortality hazard.
Setting The setting was primary care outpatient practices.
Participants The study cohort consisted of 17 707 men and 19 227 women with the diagnosis of HF who contributed only time without digoxin exposure and 9487 men and 10 808 women with the diagnosis of HF who contributed time with digoxin exposure.
Main outcome measures The main outcome measure was all-cause mortality.
Results The primary outcome of this study was the absence of a large interaction between digoxin use and sex affecting mortality. For men, digoxin use was associated with a HR for mortality of 1.00, while for women, the HR was also 1.00 (p value for interaction 0.65). The results of sensitivity analyses were consistent with those of the primary analysis.
Conclusion Observational data do not support the concern that there is a substantial increased risk of mortality due to the use of digoxin in women. This finding is consistent with previous observational studies but discordant with results from a post hoc analysis of a randomised controlled trial of digoxin versus placebo.
This is an open-access article distributed under the terms of the Creative Commons Attribution Non-commercial License, which permits use, distribution, and reproduction in any medium, provided the original work is properly cited, the use is non commercial and is otherwise in compliance with the license. See: http://creativecommons.org/licenses/by-nc/2.0/ and http://creativecommons.org/licenses/by-nc/2.0/legalcode.
Statistics from Altmetric.com
Digoxin is used in patients with HF and has been shown in one major randomised control trial, the Digitalis Investigation Group study, to reduce the rate of hospitalisations in that population.
Post hoc analysis of Digitalis Investigation Group indicated that digoxin, when used in the treatment of HF, may increase mortality by approximately 20% in women but not in men. Further randomised trials evaluating the interaction between digoxin and sex have not emerged.
For men, digoxin use was associated with a HR for mortality of 1.00, while for women, the HR was also 1.00 (p value for interaction 0.65). There was no evidence of a different association between digoxin use and mortality in women compared with men. Sensitivity analyses did not affect this estimate materially.
An interesting incidental finding of this study is that interventions known to reduce mortality in HF are used less in women than in men who have been diagnosed with HF.
Strengths and limitations of this study
The major strength of this study was its large sample size, which allowed adjustment for many covariates and numerous sensitivity analyses, none of which affected the conclusions.
The major limitation of this study is that it is non-randomised. Although there is no evidence of substantial confounding of the main study result, confounding could still bias these results.
It has been hypothesised that digoxin, when used in the treatment of heart failure (HF), may increase mortality by approximately 20% in women but not in men. This hypothesis is based on a post hoc analysis in 6800 patients by Rathore and colleagues of the Digitalis Investigation Group trial, a placebo-controlled randomised trial that showed digoxin did not affect overall mortality but did reduce hospitalisations in patients with HF.1 The post hoc analysis examined mortality effects by sex and found that compared with placebo digoxin conferred reduced mortality in men (absolute difference −1.6%, 95% CI −4.2% to 1.0%) and increased mortality in women (absolute difference 4.2%, 95% CI −0.5% to 8.8%), with a statistically significant interaction (p value =0.034).2 One proposed mechanism is that women may have higher mean serum digoxin levels than men. There is evidence that even within the therapeutic range, higher serum digoxin levels are associated with increased mortality, and experts have argued that maintenance of serum digoxin levels at the low end of the therapeutic range may be the key to safe, effective use of the drug in either sex.3–7
Because the sex–digoxin interaction was based on post hoc analyses, it needs cautious interpretation. Yet, if true, the finding is clinically important because digoxin continues to be widely used by both women and men.8 It remains unclear whether digoxin should be used differently in the different sexes, and concerns about its use in women continue to appear in the literature.9 10 Unwarranted recommendations against use of digoxin in women would deprive a large population of a medicine with demonstrated ability to prevent hospitalisations.
Further randomised trials evaluating the interaction between digoxin and sex have not emerged. One observational study based on the Studies of Left Ventricular Dysfunction cohort found no difference in digoxin's effect on mortality between sexes but concluded that additional research in other populations is still needed.10 Limitations of that study included restriction to a population with severely reduced ejection fraction (<35%) and relatively small sample size (n=6797). Another observational study performed on 2841 propensity score-matched patients with HF also did not report any difference in digoxin's effect on mortality between the sexes; this study was primarily limited by a relatively small size and restriction to a single clinical center.11 We sought to conduct a much larger study in a broader population that could help assess whether digoxin increases mortality in women compared with men.
We conducted a retrospective cohort study to test the primary null hypothesis that sex was not an effect modifier for the association of digoxin use with mortality. We also conducted planned secondary analyses of whether effect modification was mediated by digoxin dose or by serum digoxin concentration. This study was approved by the Institutional Review Board of the University of Pennsylvania. Analyses were conducted using SAS V.9.1.
The study used the Health Improvement Network (THIN) Database, a primary care electronic medical record database in the UK. This database contains over 5 million individuals who have contributed person-time from over 300 different general practices from 1986 to 2008. THIN includes demographic information on patients, as well as records of prescribed drugs, medical diagnoses (coded as READ codes), as well as vital signs and laboratory values on a subset of patients. It is a representative subset of the UK's general population.12
The analysis was restricted to individuals with a diagnostic code for HF. Once that condition was met, patients could contribute follow-up time as long as they were being given regular prescriptions of at least one drug consistent with the treatment for HF (as described below). Information was only used if it was collected during times when the adjusted mortality ratio, a quality control measure used in THIN to identify practices which are recording deaths accurately, was within acceptable standards for a given practice.13 It was possible for one individual to contribute both exposed and unexposed time. The result was four cohorts with some within-sex cross-over: men on digoxin, men not on digoxin, women on digoxin and women not on digoxin.
Exposed subjects were defined as individuals who carried the diagnosis of HF and were receiving two or more consecutive prescriptions for digoxin. Individuals were assumed to be on digoxin from the date of receipt of a prescription until 30 days after their prescription ran out (the intended duration of a prescription was either computed from information on daily dose and number of tablets or imputed as 30 days if that information was not available).
Unexposed subjects were defined as individuals who carried a diagnosis of HF and were receiving two or more consecutive prescriptions for a loop diuretic, β blocker, ACE inhibitor or angiotensin receptor blocker. They only contributed person-time to the cohort while receiving these drugs; they were assumed to be on a drug from the date of receipt of a prescription until 30 days after their prescription ran out, which was defined with the same criteria used for digoxin exposure. These are drugs that, while not specific for HF, are commonly used for that condition.
Definition of outcome
The outcome for this study was all-cause mortality, as documented by a date of death in the THIN demographics file.
Definition of covariates
Baseline covariates consisted of diagnostic codes for common comorbidities of HF, including hypertension, stroke, myocardial infarction, chronic obstructive pulmonary disease, diabetes and atrial fibrillation. We also assessed age, sex and baseline use of drugs commonly used in HF, notably aspirin, statins, oral anticoagulants, potassium-sparing diuretics and the HF-associated drugs listed above. To be included as a baseline variable, drug use had to be in the year prior to HF diagnosis; medical history such as a history of hypertension could be coded at any time prior to HF diagnosis. Blood pressure, body mass index (BMI) and a variety of laboratory values including serum creatinine were available on a minority of participants and were used as baseline variables in sensitivity analysis. Once a patient began to contribute exposure time, baseline variables were no longer updated.
Because of the large sample size, we anticipated that even clinically insignificant differences between baseline variables could be statistically significant and therefore described baseline variables using standardised differences.14 Standardised differences are a balance diagnostic used to assess how similar groups are at baseline—conceptually, they are equal to the difference of a variable's mean between two groups divided by the SD for that variable. A standardised difference of <0.1, by convention, is considered to indicate good balance of groups on that variable.14 For categorical variables, a standardised difference between men and women was calculated as , where p = proportion of the population falling into a category. In addition, because the most relevant question was whether these between-sex differences were consistent in both digoxin-exposed and digoxin-unexposed groups, this equation was extended to treat the inter-sex differences for the exposed and unexposed groups as variables which in turn had their own standardised differences:
Total follow-up time for all individuals was assessed, with digoxin use, digoxin dose and serum digoxin concentrations as time-varying covariates in secondary analysis. In the primary analysis, digoxin dose and serum concentration were ignored. A Cox proportional hazards model was used to calculate HRs for sex, digoxin use and the interaction of sex and digoxin use, controlling for the covariates specified above.15
The primary parameter of interest was the adjusted HR for the interaction of sex and digoxin use, where a significant deviation from the null would indicate that digoxin use was associated with mortality differently for men than for women. Specifically, a HR for interaction significantly greater than one would indicate that digoxin was associated with greater mortality in women than in men. This analysis was then repeated with inclusion of digoxin dose and then serum digoxin level in the model to assess for any evidence that these variables mediated any interaction. Finally, sensitivity analyses were conducted in which baseline blood pressure, BMI and serum creatinine were included in the analysis, in which the cohort was restricted to participants using concomitant loop diuretics, in which comorbid conditions were defined based only on data from the 1 year preceding cohort entry, in which patients with a diagnosis of atrial fibrillation were excluded, in which data were included even if the adjusted mortality ratio at the time of collection was not acceptable and in which subjects were excluded as soon as they crossed over from one exposure group into another.
Characteristics of the study population
The study cohort consisted of 17 707 men and 19 227 women who contributed time only without digoxin exposure and 9487 men and 10 808 women who contributed time with digoxin exposure. Table 1 presents differences in baseline variables between these groups. In general, many variables were differently distributed between sexes, as reflected by high standardised differences. However, sex differences were consistent between digoxin-exposed and unexposed groups, as reflected by low standardised differences of differences (<0.1 for all variables).
Baseline drug utilisation data were notable for high rates of loop diuretic use (table 1). Oral anticoagulant use was much more common among digoxin users even after controlling for baseline diagnosis of atrial fibrillation. Examination of trends in drug use during the year after the diagnosis of HF was documented (figure 1) shows rapid increases in the use of ACE inhibitors, β blockers and spironolactone, with higher rates of use of all these agents in men compared with women (figure 1A,B).
Blood pressure was recorded for the majority of subjects, while other quantitative covariates such as BMI and serum creatinine were recorded on a minority of subjects. In general, rates of missing data were higher for women than for men. Baseline systolic BP mean >140 was considerably more common in women, while a serum creatinine greater than 150 mg/dl was considerably more common in men.
Digoxin prescribing patterns over time
We found that digoxin use among subjects in this study remained prevalent at about 25% throughout the time period covered by the database (figure 1A,B).
Digoxin dose was consistently lower in women than in men (mean daily dose 136 μg in women and 159 μg in men, p<0.001), and the average dose in both groups declined over time by 1.6 μg/year in men and 1.3 μg/year in women (p<0.001 for both trends but no significant difference in the rate of decline for men vs women).
Serum digoxin levels were available for 8563 separate measurements covering 4063 individuals. The mean serum level for men was 1.2 ng/dl, while for women, it was 1.3 ng/dl (p<0.0001 for difference). Linear regression of digoxin level against date yielded an annual decline in mean digoxin level of 0.02 ng/dl/year in men and 0.01 ng/dl/year in women (p<0.0001 for men and p=0.01 for women, non-significant difference in rate between men and women).
Death rates after HF diagnosis without adjustment showed similar rates of mortality in all four cohorts defined by sex and digoxin use (table 2).
Multivariable modelling was pursued in stages (table 3). Briefly, we found that when universally available variables were adjusted for as covariates, there was no evidence of interaction between sex and digoxin use. The adjusted HR for digoxin exposure in a cohort restricted to women was 1.00 (95% CI 0.96 to 1.06) and for a cohort restricted to men it was identical: 1.00 (95% CI 0.95 to 1.06). A fully adjusted Cox proportional hazards model with a term for the interaction between digoxin use and sex yielded an interaction HR of 1.02 (95% CI 0.95 to 1.09), where a HR significantly >1 would have indicated that digoxin use in women conferred a greater hazard of death than in men.
We examined whether the dose of digoxin mediated any interaction between sex and exposure by including digoxin dose in the model both as an independent variable and as an interaction term with sex, both as a categorical and a continuous variable. The same was done with serum digoxin levels. None of these analyses yielded any significant interaction terms. Low levels of digoxin in the serum (<0.9 ng/dl) were associated with lower mortality than higher levels (HR 0.71, 95% CI 0.67 to 0.76, p=0.007), and this association did not differ by sex.
Sensitivity analyses including restriction of the cohort to persons actively taking loop diuretics, permitting use of data without the adjusted mortality ratio being up to standard, analysing the subset of the population with blood pressure, BMI and serum creatinine available, defining comorbid conditions based only on data from the 1 year preceding cohort entry, excluding patients with a diagnosis of atrial fibrillation and excluding persons at the time of cross-over to a different exposure category, all confirmed the results of the primary analysis, never showing a substantial or statistically significant interaction between sex and digoxin and mortality.
The primary finding from this study was the absence of a large interaction between digoxin use and sex affecting mortality, with a 95% CI excluding a HR >1.09 for the interaction. This suggests that the association between digoxin and mortality is similar in women and men. Sensitivity analyses did not affect this estimate materially.
The other prespecified aims of the study were to assess whether any interaction might be mediated by digoxin dose or by higher serum digoxin levels, both of which are thought to potentially affect mortality rates.3–7 Analyses incorporating digoxin dose and serum digoxin levels showed no evidence of such mediation. We did note an association between lower serum digoxin levels and lower mortality, but since such an association could easily be strongly confounded (eg, if more symptomatic patients were pushed to higher serum digoxin levels), this paper does not directly support or refute the claim that lower serum levels of digoxin lead to lower mortality risk.
An interesting incidental finding was rapid increases in the rates of β blocker, ACE and spironolactone use beginning at the same time that major studies were published establishing that these drugs prevented mortality in HF.16–18 Comparison by sex shows that while overall trends are similar, women consistently have lower levels of use of these life-saving agents compared with men.
A major strength of this study is its large sample size, permitting precise point estimates of associations. The size of the database also permitted numerous sensitivity analyses on subpopulations. The lack of interaction between sex and digoxin use was robust to adjustment for numerous potential confounders and to all sensitivity analyses. The major limitation of this study is that it is non-randomised. Although there is no evidence of substantial confounding of the main study result, confounding could still bias these results. Access to additional information that was not recorded in THIN, particularly the type (systolic vs diastolic) and severity of HF, would have been helpful but only a randomised trial (ideally a series of such trials in various populations and dosing regimens) could definitively avoid concerns about confounding.
The outcome of this study is discordant with the post hoc analysis of the DIG trial2 but consistent with two prior observational studies that did not observe any digoxin–sex interaction.10 11 It complements those prior studies by having over eight times their sample sizes and a more diverse cohort. Because the finding of an interaction from the DIG trial was the result of post hoc analysis, it is conceivable that the finding was a type 1 error (false positive), but it is also possible that the DIG trial results were correct, and the observational results are biased by unmeasured confounders that affected the interaction analysis. An interesting incidental finding of this study is that interventions known to reduce mortality in HF are used less in women than in men who carry that diagnosis (figure 1A,B).16–18 It would be premature to conclude that these differences necessarily imply worse care since there are important differences between the male and female populations with HF that might legitimately affect prescribing practices. However, the systematic differences in drug usage between the sexes deserve further investigation.
In conclusion, this study did not identify a difference between the sexes in the hazard of death associated with the use of digoxin. These results are of use in the context of a clinical question—whether digoxin can be used as safely in women as it can in men—for which there are few randomised data. These results suggest that this drug, with its proven ability to reduce the need for hospitalisation, is still a viable therapeutic option in women with HF.
Review history and Supplementary material
To cite: Flory JH, Ky B, Haynes K, et al. Observational cohort study of the safety of digoxin use in women with heart failure. BMJ Open 2012;2:e000888. doi:10.1136/bmjopen-2012-000888
Contributors JHF is the guarantor of this manuscript, having been primarily responsible for study design, data analysis including programming and biostatistics and manuscript composition. JHF, BK, KH, SB, JM, CR, BLS and SH all made substantial contributions to the conception and design of this study. JHF was primarily responsible for analysis and interpretation of the data, but all listed authors made significant contributions to this as well. All the listed authors made substantial contributions in drafting the article or revising it critically for important intellectual content.
Funding FOCUS Medical Student Fellowship in Women's Health supported by the Edna G. Kynett Memorial Foundation, Pharmacoepidemiology Training Funds donated to the Center for Clinical Epidemiology and Biostatistics of the University of Pennsylvania School of Medicine. The study sponsors had no influence over design of the study, the handling of data, the composition of the manuscript or the decision to submit it.
Competing interests SH has consulted for Abbott Laboratories. BLS has consulted for Abbott Laboratories, received payment for lectures from Teva Pharmaceutical Industries and has consulted for Teva and for GlaxoSmithKline. All authors on this manuscript have a connection to a training programme that has received funding from Abbott Laboratories and GlaxoSmithKline. None of these connections related directly to digoxin. These companies are among the over 50 corporations that manufacture or distribute digoxin.
Provenance and peer review Not commissioned; externally peer reviewed.
Data sharing statement No additional data available.
If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.