Abstract
Recent epidemiological studies have suggested an increased risk of venous thromboembolism (VTE) in lung fibrosis. Large-scale epidemiological data regarding the risk of VTE in pulmonary fibrosis-associated mortality have not been published.
Using data from the National Center for Health Statistics from 1988–2007, we determined the risk of VTE in decedents with pulmonary fibrosis in the USA.
We analysed 46,450,489 records, of which 218,991 met our criteria for idiopathic pulmonary fibrosis. Among these, 3,815 (1.74%) records also contained a diagnostic code for VTE. The risk of VTE in pulmonary fibrosis decedents was 34% higher than in the background population, and 44% and 54% greater than among decedents with chronic obstructive pulmonary disease and lung cancer, respectively. Those with VTE and pulmonary fibrosis died at a younger age than those with pulmonary fibrosis alone (females: 74.3 versus 77.4 yrs (p<0.0001); males: 72.0 versus 74.4 yrs (p<0.0001)).
Decedents with pulmonary fibrosis had a significantly greater risk of VTE. Those with VTE and pulmonary fibrosis died at a younger age than those with pulmonary fibrosis alone. These data suggest a link between a pro-fibrotic and a pro-coagulant state.
Idiopathic pulmonary fibrosis (IPF) is the most common of the idiopathic interstitial pneumonias, and recently published epidemiological data suggest that the burden of disease is rising. The incidence and mortality rates for IPF, and pulmonary fibrosis (PF) in general, have been increasing steadily [1–6]. IPF is a deadly disease, with a median survival after diagnosis ranging from 3 to 5 yrs [7], and no therapy other than lung transplantation has been shown to prolong survival.
In a recently conducted epidemiological study of 920 subjects, an increased risk of venous thromboembolism (VTE) was found both before and after the onset of PF [8]. In a separate study that included all Danish PF patients from 1980–2007, a prior history of VTE was found to increase the risk of developing interstitial lung disease (ILD) [9]. Although these studies hint at a relationship between VTE and PF, large-scale epidemiological data about the risk of VTE in PF-associated mortality have not been published.
At the molecular level, activation of the coagulation cascade, via the tissue factor-dependent extrinsic pathway, and as evidenced by elevated levels of tissue factor and fibrin (the main component of clots), occurs in lung tissue from patients with IPF [10]. Data also suggest abnormalities in intrinsic anticoagulation pathways: bronchoalveolar lavage fluid from patients with IPF reveals decreased protein C activation [11]. Furthermore, fibrinolysis appears to be dysregulated: models of lung fibrosis have revealed elevated levels of plasminogen activator inhibitor-1, thus promoting fibrin persistence [12].
Based on these data, we hypothesised that an excessive number of VTEs would be found in patients with IPF in decedents in the USA from 1988–2007, and that among IPF decedents, the presence of VTE would be associated with a younger age of death.
METHODS
Database
Details of the database and methods used have been published previously [6]. Briefly, we used the US Multiple Cause-of-Death mortality database (National Centre for Health Statistics (NCHS), Centers for Disease Control and Prevention, Atlanta, GA, USA) [13], for which the NCHS compiles data from all death certificates in the USA and releases the figures in annual, public-use files. For this study, we analysed files from 1988 to 2007 [13, 14]. An annual file contains >2 million decedent records, and each record contains decedent demographics, ≤20 conditions related to death within a section of the record the NCHS calls the “record axis”, and the ultimate underlying cause of death (UCD), defined by the World Health Organization as “the disease or injury which initiated the train of events leading to death” [14, 15].
From 1988–1998, the NCHS coded conditions related to death using the International Classification of Diseases (ICD)-9 [16]. Thereafter, the NCHS coded conditions related to death using the ICD-10 [14]. All data contained in these database files have been de-identified and are publicly recorded; therefore, institutional review board approval for this study was not required.
Case definitions
We included files from any decedent from the USA with IPF, as defined by the following ICD codes: from 1988 to 1998, ICD-9 codes 515 (post-inflammatory pulmonary fibrosis (PIPF)) and 515.6 (IPF); and thereafter, ICD-10 code J84.1 (a code that combines both PIPF and IPF). We then excluded patients with record axis codes for any condition that might be associated with secondary PF, including connective tissue diseases, radiation fibrosis, asbestosis, pneumoconiosis (including coal workers’ pneumoconiosis, silicosis, talcosis and berylliosis), sarcoidosis and/or extrinsic allergic alveolitis (hypersensitivity pneumonitis). In addition, we excluded patients with a concurrent diagnosis of chronic obstructive pulmonary disease (COPD) or lung cancer. We identified three groups for comparison purposes: 1) the background population, which included decedents without IPF, lung cancer or COPD; 2) decedents with COPD (but not IPF or lung cancer); and 3) decedents with lung cancer (but not IPF). VTE disease was defined as either venous thrombosis or pulmonary embolism. From 1988 to 1998, we defined venous thrombosis by ICD-9 codes 451–451.8 and 453–453.9, and pulmonary embolism by ICD-9 code 415.1, while after 1998 we used ICD-10 codes I80–I80.9, I82.8 and I82.9 to define venous thrombosis, and ICD-10 codes I26–I26.9 to define pulmonary embolism. Specifics of additional ICD-9 and ICD-10 codes are included in the Appendix.
Statistical analysis
We used the Chi-squared test (or Mantel–Haenszel statistic where appropriate) to determine the risk of VTE in decedents with IPF in reference to the three comparator groups. We then performed logistic regression to determine the adjusted risk of VTE in decedents with IPF in reference to the three comparator groups. When comparing the risk of VTE in IPF decedents with the background population, we assessed second-order interactions between IPF and sex, age and year of diagnosis. We also used the Chi-squared statistic to compare differences in the UCD between PF decedents with or without VTE. A two-sample unpaired t-test was used to compare the mean age of death between PF decedents with or without VTE. A p-value of <0.05 was considered to represent statistical significance in all analyses. All data were analysed using SAS® version 9.1 (SAS Institute, Cary, NC, USA).
RESULTS
From 1988 to 2007, there were 46,450,489 deaths of residents of the USA. A total of 218,991 records met our criteria for a diagnosis of IPF. Among these, 3,815 records also contained a diagnostic code for VTE (table 1), yielding a prevalence of VTE in IPF that was significantly greater than in the background population (1.74% versus 1.31%; p<0.0001). Thus, the overall odds of being diagnosed with a VTE were significantly greater in decedents with IPF than in decedents in the background population (overall OR 1.34, 95% CI 1.29–1.38). In the year-by-year analysis, the risk of VTE was greater for IPF than in the background population in every year except 2005 (p=0.25).
After stratifying for sex and age, female decedents with IPF in age strata 65–74 yrs, 75–84 yrs and ≥85 yrs, but not <65 yrs, had a significantly greater risk of VTE than the background population (fig. 1). Among male decedents with IPF, the risk of VTE was significantly greater than the background population for all age strata (fig. 2).
Using logistic regression, and adjusting for sex, age and year of death, the risk of VTE remained significantly greater in decedents with IPF than in the background population (adjusted OR 1.34, 95% CI 1.29–1.38). Among all decedents, males had a significantly lower risk of VTE (OR 0.75, 95% CI 0.74–0.75) than females, and increasing age conferred a slightly lower, but significant, risk of VTE (OR 0.99, 95% CI 0.99–0.99). From 1988 to 2007, the overall risk of VTE increased slightly with each passing year (OR 1.01, 95% CI 1.01–1.01).
When interactions between IPF and sex, age and year of death were assessed, each two-way interaction was found to be statistically significant. The risk of VTE in decedents with IPF was modified by age, sex and year of death: among decedents with IPF, increasing age, female sex and later year of death reduced the risk of VTE. β-coefficients, standard errors and p-values for the model that included these second-order interactions are reported in table 2.
We assessed the effects of race and sex on the risk of VTE in decedents with IPF compared with the background population. Using logistic regression, we found that the risk of VTE remained greater for IPF decedents than those in the background population when adjusted for ethnicity and race (adjusted OR 1.39, 95% CI 1.34–1.43). Among all decedents, when using non-Hispanic Whites as the reference population, Hispanics and non-Hispanic others (including Asians) had a lower risk of VTE (OR 0.75, 95% CI 0.74–0.76 and OR 0.56, 95% CI 0.54–0.57, respectively), while non-Hispanic Blacks had an elevated risk of VTE (OR 1.33; 95% CI 1.32–1.34). There were no significant two-way interactions between IPF and ethnicity or race, suggesting that the risk of VTE among IPF decedents is not further altered by the decedent’s race or ethnicity.
There were 3,980,364 decedents with COPD. Of these, 48,441 records also contained a code for VTE, yielding a prevalence of VTE that was significantly lower for COPD than for IPF (1.22% versus 1.74%; p<0.0001) and a risk of VTE that was significantly greater in decedents with IPF (overall OR 1.44, 95% CI 1.39–1.49) than in those with COPD (fig. 3 and table 3). The risk of VTE was greater for IPF than COPD in every year studied. Using logistic regression and adjusting for age (OR 0.98, 95% CI 0.98–0.99), male sex (OR 0.96, 95% CI 0.94–0.97) and year of death (OR 1.01, 95% CI 1.01–1.01), the risk of VTE in decedents with IPF was found to be significantly greater than in those with COPD (adjusted OR 1.41, 95% CI 1.34–1.46).
There were 3,233,845 decedents with lung cancer. Of these, 36,876 records also contained a code for VTE, yielding a prevalence of VTE that was significantly lower than for IPF (1.14% versus 1.74; p<0.0001) and a risk of VTE that was significantly greater in decedents with IPF (OR 1.54, 95% CI 1.49–1.59) than in those with lung cancer (table 4). The risk of VTE was greater for IPF than lung cancer in every year except 2004 (p=0.53) and 2005 (p=0.65) (fig. 4). Using logistic regression and adjusting for age (OR 0.98, 95% CI 0.98–0.98), male sex (OR 0.96, 95% CI 0.96–0.99) and year of death (OR 1.04, 95% CI 1.04–1.04), the risk of VTE in decedents with IPF was found to be significantly greater than in those with lung cancer (adjusted OR 1.66, 95% CI 1.61–1.72).
Regardless of sex, decedents with IPF and VTE were significantly younger at the time of death than those with IPF but without VTE (females: 74.3 versus 77.4 yrs (p<0.0001); males: 72.0 versus 74.4 yrs (p<0.0001)). IPF decedents with VTE were more likely than those without VTE (73.2% versus 69.3%; p<0.0001) to have either VTE or IPF coded as the UCD. All other candidate conditions examined were less likely to be the UCD in the IPF decedents with VTE than in those without VTE (table 5).
Compared with the background population, the overall risk of VTE was lower for those with lung cancer (1.14% versus 1.31%; OR 0.87, 95% CI 0.86–0.88) or COPD (1.22% versus 1.31%; OR 0.93, 95% CI 0.92–0.94). The risk of VTE was greater among decedents with COPD than those with lung cancer (lung cancer: 1.14% versus 1.22%; OR 1.07, 95% CI 1.05–1.08).
DISCUSSION
Using mortality data from all USA decedents between 1988 and 2007, we found that the risk of VTE in those with IPF at the time of death was ≥34% higher than the risk of VTE in the background decedent population. We also found that the overall risk of VTE in decedents with IPF was 54% greater than for decedents with COPD and 44% greater than for decedents with lung cancer, two conditions previously identified to confer an increased risk of VTE [17–19].
Only a few studies have examined the relationship between VTE and IPF. Using primary care data from the UK from 1991–2003, Hubbard et al. [8] identified a cohort of 920 patients with IPF and found that, in comparison to the general population, the risk of deep-venous thrombosis was greater before (OR 1.98, 95% CI 1.13–3.48) and even greater after (risk ratio 3.39, 95% CI 1.57–7.28) the diagnosis of IPF. Using hospital discharge and mortality data from 1980–2007, Sode et al. [9] examined the association between VTE and ILD in the entire Danish population. They found an increased risk for the development of ILD in people ever diagnosed with VTE, many of whom probably had IPF. In contrast to our study, they did not exclude subjects with known-cause PF, including those with conditions that might confound the relationship between VTE and PF (e.g. connective tissue diseases and obstetric conditions).
The risk of VTE in decedents with IPF was significantly greater than the background population for females >64 yrs of age and for males in every age group. In fact, male decedents with IPF appeared to be at particularly high risk of developing VTE. The reasons for these findings are unclear. Although younger females are in a pro-coagulant state and have an elevated risk of developing VTE due to hormonal influences from the menstrual cycle, pregnancy and contraceptives, we believe that in certain individuals, particularly beyond a certain age (apparently in males more than females), the fibrotic process confers an even more profound risk for VTE [20]. Although among all decedents in the USA from 1988–2007 the risk of VTE differed across race/ethnic groups, among those with IPF, the risk of VTE was not further affected by a decedent’s race or ethnicity.
Somewhat surprisingly, we found a significantly greater risk for VTE in IPF than either lung cancer or COPD. Using a case–control design, Blom et al. [21] studied subjects within a Dutch population and found a nearly 25-fold increased risk of VTE in patients with lung cancer when compared with controls. Levitan et al. [18] analysed a Medicare database and found that VTE occurred more frequently in lung cancer patients (6.1 events per 1,000 lung cancer patients) than in patients with various other malignancies. In another study, investigators observed that among patients who required hospital admission for an exacerbation of COPD, 25% had a pulmonary embolism [19]. Although the risk of VTE is elevated in patients with lung cancer, and emerging data suggest the same is true for certain patients with COPD, our results place IPF patients at even greater risk for mortal VTE events. The data leave little room for questioning the clinical significance of these VTE events: decedents with IPF and VTE died at a younger age than decedents with IPF alone, and in >20% of IPF decedents with VTE, VTE itself was the UCD.
We found that the risk of VTE in both decedents with lung cancer and those with COPD was significantly lower than the background population. There are several possible explanations for this observation: the first is that the background population indeed possessed a greater risk for mortality-related VTE than decedents with COPD or lung cancer; however, because the database did not allow us to determine the presence or absence of most VTE risk factors (i.e. hereditary causes and other acquired causes of thrombophilia including tobacco abuse) or VTE therapies (i.e. anticoagulation), we could not account for these potential confounders, in any of the subgroups, in our analyses. The data for the current study concern VTE events at the time of death; thus, we know nothing of the number, type, severity or treatment of VTE events that occurred earlier in the lives of people with COPD or lung cancer, or those in the background population. Our data should not be used to make inferences about life-long risk of VTE.
Not only is lung cancer a known risk for VTE, but once VTE is diagnosed in a patient with malignancy, current guidelines recommend that patients receive anticoagulation for life or until the malignancy has resolved [22]. This practice could reduce the rates of VTE at the time of death in patients with lung cancer. As more data emerge regarding the risk of VTE in patients with COPD, similar trends may ensue and, thus, may also explain the slightly, but significantly, higher risk of VTE in patient with COPD than lung cancer [19].
Finally, we suspect that in a great number of decedents with lung cancer, many conditions possibly contributing to death (e.g. VTE) do not get mentioned on death certificates: the underlying malignancy simply overshadows them, and death certifiers do not write them down. We suspect the same is true for decedents with COPD: at least historically, patients with COPD suffering a respiratory-related mortal event were overwhelmingly likely to have that event attributed to an acute exacerbation of COPD; perhaps an alternative cause (e.g. VTE) would not have even been investigated.
If this is true, why would IPF patients be at risk for developing VTE? Mounting evidence suggests that the IPF microenvironment is both pro-coagulant and antifibrinolytic, and that components of these haematological pathways contribute to the milieu that drives the fibrotic process [23]. In fact, pro-coagulant moieties (including elevated levels of protein C and reduced thrombomodulin) have been found in the systemic circulation of some patients with IPF during abrupt disease accelerations: so-called acute exacerbations [24]. These data suggest that IPF is a hypercoaguable state.
Recognising that IPF may be a hypercoaguable condition and that components of the coagulation cascade are either directly involved in or activated by fibrosis-producing machinery, Kubo et al. [25] tested the hypothesis that inhibiting the coagulation cascade would lead to improved outcomes in patients with IPF. They randomised 56 hospitalised Japanese patients with IPF to receive either prednisolone alone or prednisolone plus oral warfarin. Although the study did not meet its primary outcome, subjects receiving warfarin had a 3-yr survival of 63% compared with 35% for subjects receiving only prednisolone (p=0.04). In addition, when the authors examined acute exacerbations within this cohort, they found that the resultant mortality from an acute exacerbation (as well as d-dimer levels) was significantly lower in the group receiving warfarin therapy. IPFnet, a network of institutions in the USA sponsored by the National Institutes of Health and charged with conducting trials of therapy for IPF, is currently conducting a randomised, placebo-controlled trial of warfarin for the treatment of IPF.
Like other investigators using ICD-coded mortality data, we faced several limitations in this study. These were largely imposed by the dataset we chose to use. First, we had to rely on death certifiers to correctly identify cases and then code them appropriately. In some decedents with IPF (e.g. those whose IPF goes undiagnosed and those who are misdiagnosed as having some other lung disease), the correct diagnosis of IPF will not make it onto the death certificate. For these reasons, we believe, and others have previously demonstrated [3, 26], that PF is under-reported on death certificates. Furthermore, we could not assess the accuracy of a VTE diagnosis, and the modalities used to diagnose VTE are not captured in this dataset. Clinically significant VTE is probably under-reported when clinical criteria are used and probably over-reported when autopsy data are used [27]. Indeed, pulmonary emboli are commonly found at autopsy, and in 50–70% of cases, they were not suspected by the clinician [28]. If VTE was either under- or over-reported in this database, we would expect any such misclassification to be independent of whether a person had IPF or not, thus biasing the association toward the null. However, it remains possible that misreporting occurred: we have no way to determine this.
Acute exacerbations of IPF (AE-IPF) have only recently gained widespread recognition in the USA [29]. Given this, in the past, death certifiers could have labelled IPF patients who abruptly declined and died with VTE (pulmonary embolism), a reasonable potential explanation for acute decompensation and death. However, the year-by-year analysis in our study argues against this: decedents with IPF were more likely than decedents in the background population to have VTE in every year studied, including years since the entity AE-IPF first gained recognition in the USA. If, in more recent years, decedents were being “pulled” from the VTE-as-UCD group and placed in the IPF-as-UCD group, we would have expected a steady decline in VTE rates among IPF decedents; this was not the case.
The database restricted us from identifying risk factors for VTE, including the number of episodes of VTE any decedent had in his life, heritable hypercoaguable conditions, smoking status, and the use of oral contraceptives or female sex hormones. The fact that the association between IPF and VTE held in elderly females (a group perhaps less likely to be prescribed oral contraceptives) and in males argues strongly against any confounding by female sex hormone use.
In summary, using all death certificate records in the USA from 1988–2007, we found that the risk of VTE in decedents with IPF was significantly greater than in: 1) decedents in the background population, 2) those with lung cancer or 3) those with COPD. These findings add to a growing data pool on the links between pro-fibrotic and pro-coagulant pathways, but should be viewed as hypothesis generating. Nonetheless, we believe VTE should be considered in any IPF patient whose respiratory status is declining, particularly those in whom the decline is abrupt or no other cause can be identified. Research is needed to discern the interdependencies of the pro-fibrotic and pro-coagulant cascade, determine whether IPF induces a systemic hypercoaguable state, and discover whether therapeutic disruption of coagulation leads to improved outcomes in patients with IPF.
Footnotes
Support Statement
J.J. Swigris is supported, in part, by a National Institutes of Health (Bethesda, MD, USA) Career Development Award (K23).
Statement of Interest
Statements of interest for A. Fischer and J.J. Swigris can be found at www.erj.ersjournals.com/site/misc/statements.xhtml
- Received March 7, 2011.
- Accepted June 17, 2011.
- ©ERS 2012