Abstract
Residual pulmonary vascular obstruction (RPVO) and chronic thromboembolic pulmonary hypertension (CTEPH) are both long-term complications of acute pulmonary embolism, but it is unknown whether RPVO can be predicted by variants of fibrinogen associated with CTEPH.
We used the Akaike information criterion to select the best predictive models for RPVO in two prospectively followed cohorts of acute pulmonary embolism patients, using as candidate variables the extent of the initial obstruction, clinical characteristics and fibrinogen-related data. We measured the selected models’ goodness of fit by analysis of deviance and compared models using the Chi-squared test.
RPVO occurred in 29 (28.4%) out of 102 subjects in the first cohort and 46 (25.3%) out of 182 subjects in the second. The best-fit predictive model derived in the first cohort (p=0.0002) and validated in the second cohort (p=0.0005) implicated fibrinogen Bβ-chain monosialylation in the development of RPVO. When the derivation procedure excluded clinical characteristics, fibrinogen Bβ-chain monosialylation remained a predictor of RPVO in the best-fit predictive model (p=0.00003). Excluding fibrinogen characteristics worsened the predictive model (p=0.03).
Fibrinogen Bβ-chain monosialylation, a common structural attribute of fibrin, helped predict RPVO after acute pulmonary embolism. Fibrin structure may contribute to the risk of developing RPVO.
Abstract
Fibrinogen helps predict residual pulmonary vascular obstruction after pulmonary embolism: results from the PROMPT study http://ow.ly/G7pY30m2HU6
Introduction
After acute pulmonary embolism, thrombotic material in the pulmonary arteries typically disintegrates, although in some cases it transforms into permanent obstructing scars [1]. Residual pulmonary vascular obstruction (RPVO) of ≥10% is common after pulmonary embolism and can diminish physical activity, reduce quality of life and worsen overall prognosis [2–7]. RPVO is associated with dyspnoea [8], hypoxaemia [9], gas exchange deficits [9], exercise intolerance [2] and other serious clinical impairments [9–11]. RPVO incurs the risk of recurrence [12] and of chronic thromboembolic pulmonary hypertension (CTEPH) [2, 13], an uncommon [14, 15] but severe complication of pulmonary embolism [16, 17] in which obstruction is extensive and pulmonary arterial pressure rises [18].
The mechanisms by which fibrin-rich thromboemboli transform into permanent obstructions characteristic of RPVO and of CTEPH are incompletely understood. RPVO is not associated with “thrombophilia” or with impaired fibrinolytic enzymes [2]. We observed that fibrin itself derived from some patients with CTEPH forms clots resistant to lysis [19, 20] with aberrations in molecular configuration [21] and clot structure [22] that reflect otherwise rare fibrinogen gene mutations [20]. CTEPH-associated dysfibrinogenemia may be a model for more common fibrin characteristics that might tip the balance from clot resolution towards RPVO.
To gain insights into the development of RPVO, the Prediction of Residual Obstruction Manifested after Pulmonary Embolism Treatment (PROMPT) study followed two cohorts of patients after acute pulmonary embolism to determine the best predictors of RPVO. We investigated whether fibrinogen post-translational modifications and other common fibrin structural attributes might help predict the development of RPVO after acute pulmonary embolism. If fibrin properties are consistently associated with RPVO, then they might be implicated as a mechanism of its development.
Methods
Study design and study population
We used clinical data and biological material from two single-centre, prospective, observational cohort studies conducted on separate groups of acute pulmonary embolism patients in an academic teaching hospital (Hôpital Européen Georges Pompidou, Paris, France). The two studies were approved by the local ethics committee. All patients provided written informed consent before enrolment.
For cohort 1, we enrolled consecutive patients with pulmonary embolism from our centre who had been included in the FARIVE study, a multicentre case–control study that evaluated the interactions of environmental, genetic and biological risk factors on the risk of a first venous thromboembolism (VTE). The study included consecutive subjects who were aged >18 years and had their first symptomatic occurrence of pulmonary embolism. Patients were excluded if they had active cancer, unless they had been treated successfully and had no evidence of recurrence in the previous 5 years. In addition, patients were excluded if they had short life expectancies secondary to other associated pathologies, anticipated impossibility of follow-up or a previously diagnosed thrombophilia.
For cohort 2, we enrolled consecutive patients with pulmonary embolism followed-up at our centre between February 1999 and May 2006 who had not been enrolled in the FARIVE study [2]. This study was designed to assess the clinical significance of residual pulmonary vascular obstruction after an episode of acute pulmonary embolism. Patients were included if they were aged >18 years and had completed ≥3 months of oral anticoagulant therapy without having recurrent pulmonary embolism.
In both cohorts, the diagnosis of pulmonary embolism was established through high-probability lung scintigraphy, thoracic computed tomography (CT) scan, pulmonary angiography, sonography showing the existence of a proximal venous thrombosis in the absence of a differential diagnosis for the reported pulmonary symptoms or nondiagnostic pulmonary scintigraphy associated with lower extremity ultrasound disclosing proximal deep vein thrombosis. The type and duration of long-term therapy was left to the discretion of physician in charge of the patients. Standard practice was to treat provoked pulmonary embolism for ≥3 months and unprovoked pulmonary embolism for ≥6 months.
Residual pulmonary vascular obstruction
All patients underwent a standardised ventilation/perfusion ratio (V′/Q′) lung scan ≥6 months after the acute pulmonary embolism [23]. A trained investigator (BP) blinded to all other study data identified and scored mismatched perfusion defects according to a prespecified protocol [24]. We defined RPVO as a residual vascular obstruction of ≥10%, which corresponds to an amputation of at least two pulmonary segments (the minimal obstruction defining a high-probability result on diagnostic V′/Q′ lung scans [25]) and is associated with dyspnoea and exercise intolerance after acute pulmonary embolism [2].
Demographic and clinical data
A standardised tool was used upon study enrolment to record clinical data (including the type of acute therapy and the duration of long-term therapy) and conduct a structured interview. Except for the fibrinogen-related laboratory data described in the following section, all laboratory values were measured in the clinical laboratory of Hôpital Européen Georges Pompidou. Physicians blinded to the other data of each patient used a validated score to calculate the initial pulmonary vascular obstruction on the original diagnostic CT scan [26] or V′/Q′ scan [24]. None of the subjects had pulmonary angiography as the initial diagnostic test.
Fibrin structural attributes
Upon enrolment, blood samples were collected by venipuncture into anticoagulant-containing tubes. Plasma was separated by two consecutive centrifugations and stored at −70°C for analysis. Plasma fibrinogen levels were measured in a clinical laboratory at the time of blood collection.
Research assays were performed on fibrinogen that had been purified from citrated plasma by ethanol precipitation and exposure to Gelatin Sepharose (GE Healthcare, Pittsburg, PA, USA) to remove residual fibronectin, followed by dilution to 4 mg·mL−1 [20]. Online supplementary appendix 1 includes the detailed assay methodologies, which have been reported previously [19, 20, 22]. They are described briefly below.
Fibrinogen, thrombin and calcium chloride were combined to form fibrin clots in the amounts appropriate for each structural assay. We determined the turbidity of the fibrin clots, which corresponds to fibrin network fibril dispersal and branching as validated by microscopy [20, 22]. We measured fibrin clot susceptibility to lysis through the rates of turbidity decrease, measured consecutively after clots were exposed to a mixture of plasminogen and tissue plasminogen activator [20]. We measured fibrin clot permeability, which reflects the organisation with which fibrin polymers are formed [22], by the rate of passage of permeability buffer through fibrin clots. More organised fibrin networks with large pores have higher flow compared to disorganised fibrin networks with small pores. We measured the degree that fibrin crosslinks its α-chains using sodium dodecyl sulfate polyacrylamide gel electrophoresis as a reflection of the potential to form stabilised clots during thrombosis [19]. We used I125-labelled antibodies specific for the β15–42 peptide sequence chain (anti-β15–42) to measure the accessibility within the fibrin clots of the fibrin β-chain amino termini (residues 15–42), which are implicated in a variety of physiological events that affect thrombus remodelling into intravascular scars [27–29].
We measured precise fibrinogen masses using liquid chromatography, followed by mass spectrometry (LC/MS) to identify post-translational modifications [20]. We analysed the LC/MS spectra representing the fibrinogen Aα-chains, Bβ-chains and γ-chains and measured the proportional heights of mass peaks corresponding to variants we had observed in CTEPH-associated dysfibrinogenaemias and any other mass peaks that represented ≥10% of the associated fibrinogen chains.
Statistical analysis
Categorical variables are presented as n (%) and continuous variables are presented as mean±sd. All available demographic, clinical and fibrinogen data were compared between the two cohorts, using t-tests, Wilcoxon–Mann–Whitney tests, Fisher's exact tests or Chi-squared tests, as appropriate. Univariate assessment of predictors of RPVO was performed for each cohort, using the same tests, as appropriate.
Because the interaction of lysis, the inflammation and cellular remodelling that leads to RPVO probably develops through a complex interaction between clinical factors and biological factors, the careful selection of a predictive model was critical to our goal. A joint predictive effect would be obscured if simplistic univariate analyses were used [30]. For this reason, we used a stepwise multivariate logistic regression procedure based on the Akaike information criterion (AIC) [30, 31] to select the most parsimonious and accurate multivariable prediction model for RPVO in each cohort, using as candidate predictors both clinical risk factors and fibrin structural properties. If the best predictive models selected in this fashion depended on fibrin structural attributes, then the analysis would support the role of fibrin properties in the development of RPVO after acute pulmonary embolism. The AIC-based approach avoids selecting models solely on the basis of apparent significance of the individual terms or global fit, both of which are likely to be overstated; overfitting is penalised [31]. We used deviance-based goodness of fit, a standard logistic regression metric, to assess the quality of the selected models in each cohort.
Multivariable modelling was performed on information from subjects who had complete data (without missing values), to enable model comparisons. For cohort 1, AIC-based stepwise logistic regression was used to identify the best clinical and demographic predictors of RPVO. Then, AIC-based stepwise logistic regression was used to select the best model to predict RPVO, using as candidate predictors each patient's previously identified relevant demographic and clinical measures, as well as the results of the fibrinogen experiments we described. We used the receiver operating characteristic (ROC) curve [32] to find the best probability cut-off for predicting RPVO and assessed the model's sensitivity and specificity. We then validated the model derived in cohort 1 by using it to predict RPVO in cohort 2.
To confirm the role of specific fibrin structural attributes in RPVO, we performed the prediction model selection process again in cohort 1 and included in the candidate predictor pool only the fibrinogen properties and the initial degree of obstruction (excluding the clinical data). The initial degree of obstruction is such a strong predictor of residual obstruction [2] that all predictive models would include it. We were interested in finding additional, independently predictive measures. Finally, to further verify the role of fibrinogen, we used a Chi-squared test of residual deviances to determine whether the model that included the initial obstruction as well as fibrinogen properties were significantly better predictors of RPVO than the model that excluded fibrinogen properties.
To confirm the generalisability of individual RPVO predictors identified in cohort 1, we repeated the multivariate modelling process in cohort 2, but limited the candidate variables to those that were disclosed by the cohort 1 analysis. We reasoned that variables predictive of RPVO in both cohorts might be implicated in RPVO development in the general population of patients with acute pulmonary embolism. If a variable derived from cohort 1 was not included in the optimal prediction equation generated for cohort 2, or if its coefficient was in the opposite direction, we reasoned that it would be less likely to be predictive of RPVO in other pulmonary embolism populations.
Results
Subject populations
From January 2003 to January 2009, 137 patients were recruited into cohort 1 (figure 1a). Two (1.4%) patients were not assessed for RPVO because they were lost to follow-up; in six (4.4%), the follow-up V′/Q′ lung scans were missing. Plasma samples for fibrin(ogen) laboratory analysis from 27 (20.9%) subjects were unavailable. There were no significant differences between the 35 patients excluded from cohort 1 and the included patients from cohort 1 in age, sex, time from symptoms to diagnosis or residual obstruction. Among the 102 subjects included in the analysis, 29 (28.4%) had RPVO ≥10% (table 1). The mean±sd time from pulmonary embolism diagnosis to V′/Q′ lung scan was 8.3±2.5 months and degree of persistent perfusion defect was 8.1±11.7%.
Between February 1999 and May 2006, 185 patients were recruited into cohort 2 (figure 1b). Two (1.0%) patients were lost to follow-up. Plasma samples for fibrin(ogen) laboratory analysis were unavailable from one (0.5%) subject. The excluded patients were not outliers on any of the studied characteristics. Among the remaining 182 subjects, 46 (25.3%) had RPVO ≥10%. The mean±sd time from pulmonary embolism diagnosis to V′/Q′ lung scan was 10.9±3.7 months and mean RPVO was 6.8±12.7%.
The subject characteristics are compared between the two cohorts in table 1 (a comprehensive list is included in online supplementary appendix 2). Cohort 1 was significantly younger than cohort 2 (52.2±19.6 years versus 61.3±18.5 years, p=0.0002). Cohort 1 had a lower prevalence of cancer (2.0% versus 12.1%, p=0.003), a lower prevalence of bed rest (8.8% versus 20.3%, p=0.01), a lower use of thrombolytics (0% versus 5.5%, p=0.02) and a lower use of inferior vena cava filters (0% versus 4.5%, p=0.03) than cohort 2. However, family histories of VTE were more common in cohort 1 than in cohort 2 (31.4% versus 18.0%, p=0.01). Pulmonary embolism was diagnosed somewhat earlier after the onset of symptoms in cohort 1 than in cohort 2 (11.2±17.1 days versus 14.3±28.6 days, p=0.04).
A univariate analysis of the factors associated with RPVO is displayed in table 2. In both cohorts, the initial pulmonary vascular obstructions were smaller among those who recovered without RPVO, than in those with RPVO. Notably, a CTEPH-associated mutation (fibrinogenSanDiegoIII) [20] was found in one subject from cohort 1, who had RPVO of 29.5%.
Derivation of RPVO predictive model in cohort 1
Prediction model creation and selection was performed on 94 subjects in cohort 1 who had complete clinical and fibrin-related datasets. When all candidate variables were considered, the optimal model to predict the probability of RPVO included the initial extent of vascular obstruction, fibrinogen Bβ-chain monosialylation, the days between symptoms and diagnosis and fibrinogen-γ′ (table 3). The goodness of fit of the model was highly significant (difference between null and residual deviances=21.8, Chi-squared p=0.0002). The area under the ROC curve (AUC) was 0.762 (p=0.0001). The sensitivity was 0.679 (95% CI 0.476–0.841) and the specificity was 0.727 (95% CI 0.604–0.830).
Validation of RPVO predictive model in cohort 2
The optimal predictive model derived from all candidate variables in cohort 1 was applied to the 141 subjects in cohort 2 who had complete datasets. The AUC was 0.693 (p=0.0005). The sensitivity was 0.622 (95% CI 0.448–0.775) and the specificity was 0.692 (95% CI 0.595–0.777).
Derivation of RPVO predictive model in cohort 1 using fibrinogen properties
When the candidate predictor pool excluded the clinical and demographic data, and encompassed only the initial extent of obstruction and fibrinogen data, the best prediction model's goodness of fit remained highly significant (difference between null and residual deviances=26.2, Chi-squared p=0.00003). Initial obstruction and fibrinogen Bβ-chain monosialylation remained as predictors of RPVO in that model, which also included fibrinogen-γ′ and the fibrin lysis rate (table 4). The AUC was 0.802. However, when the candidate predictor pool also excluded the fibrinogen data, Chi-squared tests of residual deviances disclosed that the best-fit models did not predict RPVO as well as when fibrinogen data were included (difference between the two models’ deviances=9.0, Chi-squared p=0.03).
Derivation of RPVO predictive model in cohort 2 using predictor pool from cohort 1
A prediction model was created for cohort 2, with the pool of candidate predictors restricted to only the four predictors identified by the model creation and predictor selection process performed for cohort 1. When only those candidate variables were considered, the optimal model to predict the probability of RPVO included the initial extent of vascular obstruction, fibrinogen Bβ-chain monosialylation and the days between symptoms and diagnosis. Fibrinogen-γ′ was no longer included as a significant predictor. The goodness of fit of the model remained highly significant (difference between null and residual deviances=18.8, Chi-squared p=0.0009). The AUC was 0.702.
In a post hoc analysis, the median extent of Bβ-chain monosialylation among the entire (two-cohort) group was 59.2%. Among the half of patients with Bβ-chain monosialylation below the median level, 19.4% had RPVO, whereas 31.9% of those with higher levels of Bβ-chain monosialylation had RPVO (p=0.0169).
Discussion
We followed two cohorts of patients after acute pulmonary embolism to determine the factors that were most predictive of RPVO. High proportions of monosialylation of the fibrinogen Bβ-chain were consistently predictive of RPVO across the two cohorts. Among the other candidate variables, only the initial degree of vascular obstruction was more predictive of RPVO, which is in agreement with our previous observations [2].
We defined RPVO as persistent lung perfusion defects in ≥10% of the lungs (i.e. in at least two lung segments). Although we did not objectively measure quality of life or other symptomatic outcomes in the current study, we previously reported that our definition of RPVO, which occurred in 29% of acute pulmonary embolism patients, was significantly associated with dyspnoea, elevated pulmonary artery pressures and exercise intolerance [2]. Among pulmonary embolism patients with this degree of RPVO, 16% were subsequently diagnosed with CTEPH [2].
The predictors selected by the model were highly associated with the subsequent development of RPVO (p=0.0001 in the derivation cohort and p=0.0005 in the validation cohort). The diagnostic accuracies of the model were not sufficiently high to exclude of predict RPVO for clinical purposes (AUC 0.762 and 0.693 in the derivation and validation cohorts, respectively). However, the association between RPVO and fibrinogen properties provides insights into the mechanism whereby some patients develop RPVO or even CTEPH after acute pulmonary embolism. No association with genetic thrombophilia has been identified [2], which raises the possibility that impairments in the process of fibrinolysis are responsible [33]. We focused our investigation on fibrin itself, since the resistance of purified fibrin clots to lysis has been observed in some patients with CTEPH [19, 20] and, to a lesser extent, in some pulmonary embolism patients [33, 34]. We used as our model the properties of fibrin [21, 22] that we had observed in CTEPH-associated dysfibrinogaenemias [20]. We reasoned that variations in post-translational modification patterns might influence fibrin formation and lysis [35, 36] less severely, but more commonly, than CTEPH-associated dysfibrinogaenemias and that they may have a role in persistence of thromboemboli after acute pulmonary embolism.
Our results support the hypothesis that attenuation of the body's normal mechanisms for resolution of fibrin-rich thrombi is an important pathophysiological step in the aetiology of chronic syndromes after venous thromboembolism [5, 18]. One possible mechanism is that properties of the fibrin clot would delay or limit plasmin-mediated fibrinolysis or modify other clot resolution mechanisms such as angiogenic and inflammation responses [18]. Our observations suggest that structural attributes of fibrin may influence fibrin clot structure sufficiently to contribute to poor clot resolution and chronic disease after acute pulmonary embolism.
High proportions of monosialiated Bβ-chains were predictive of RPVO in both pulmonary embolism cohorts studied. Sialylation occurs during fibrinogen formation and has various effects on fibrin polymer organisation and resolution due to the electrostatic effect of the negatively charged sialic acid. Experimental sialic acid removal results in disordered fibrin clot formation [35, 37, 38], while excessive sialylation of Bβ-chains and γ-chains also results in abnormal clots [22, 36]. Monosialylation is increased when clotting might be beneficial, such as during pregnancy [39] and acute infection [40].
The results for the two cohorts we studied differed somewhat in the predictive power of low fibrinogen-γ′, which has been implicated in the regulation of fibrin cross-linking. The differences might be due to chance. Alternatively, they might reflect subtle differences between the two cohorts or even different implications of fibrin(ogen) characteristics at the time of the acute pulmonary embolism (when blood samples were collected from the first cohort) and several months later (when the second cohort was sampled).
Our study is limited by the fact that we did not independently confirm the clinical significance of RPVO in our patients by objective measurements of dyspnoea, haemodynamics or subsequent clinical outcomes after completion of the study protocol. However, we have validated the association between RPVO and those findings in a previous series of similar pulmonary embolism patients [2]. We also acknowledge that 20% of patients in the first cohort were excluded because they did not have adequate plasma samples for fibrinogen analysis, which could have introduced a potential for bias in the analysis of that cohort. While we did not disclose an association between the type of acute or long-term anticoagulation therapy and RPVO, the study was performed prior to the widespread use of direct oral anticoagulants, so we cannot comment on the influence of those therapies. Finally, despite the consistency of our findings in both cohorts who had similar rates of RPVO, they should be confirmed in subsequent prospective studies.
Conclusion
The extent of pulmonary vascular obstruction at the onset of pulmonary embolism, as well as fibrinogen Bβ-chain monosialylation, were predictive of RPVO after acute pulmonary embolism in the two cohorts studied. Our results confirm our previous observations and support the role of fibrin properties in the resolution or persistence of thrombi after acute pulmonary embolism, which is similar to the role played in the development of CTEPH. It is possible that abnormal fibrin properties could help identify acute pulmonary embolism patients at high risk for RPVO for whom adjuvant therapy to modify disease progression might be evaluated. For example, since fibrinogen properties such as fibrinogen Bβ-sialylation may be influenced pharmacologically [39], the effect on RPVO of such interventions to alter fibrinogen or the body's response to fibrin may have clinical implications.
Supplementary material
Supplementary Material
Please note: supplementary material is not edited by the Editorial Office, and is uploaded as it has been supplied by the author.
Appendix 1. Supplementary methods. ERJ-01467-2018_Appendix_1
Appendix 2. Supplementary tables. ERJ-01467-2018_Appendix_2
Acknowledgements
The authors would like to thank Russell Doolittle (University of California San Diego, San Diego, CA, USA) for his helpful advice and suggestions on the experimental methods and manuscript. Jess Mandel (University of California San Diego, San Diego, CA, USA) also provided helpful comments on the manuscript.
Footnotes
This article has supplementary material available from erj.ersjournals.com
Author contributions: B. Planquette recruited patients, performed experiments and helped write the manuscript; O. Sanchez, J. Emmerich and G. Meyer helped design the study and contributed to the manuscript; J.J. Marsh, P.G. Chiles, R.E. Moore and G.B. Gugiu performed laboratory experiments and contributed to the manuscript; G. Le Gal planned analysis and contributed to the manuscript; T. Wolfson and A.C. Gamst performed statistical analyses and contributed to the manuscript. T.A. Morris is the study primary investigator, planned the study, performed experiments and contributed to the manuscript.
Conflict of interest: B. Planquette has nothing to disclose.
Conflict of interest: O. Sanchez reports grants, personal fees and non-financial support from Bayer, grants, personal fees and non-financial support from MSD, grants and non-financial support from Actelion, personal fees and non-financial support from BMS-Pfizer, grants from Daiichi Sankyo, personal fees from Chiesi and non-financial support from Boehringer Ingelheim, outside the submitted work.
Conflict of interest: J.J. Marsh has nothing to disclose.
Conflict of interest: P.G. Chiles has nothing to disclose.
Conflict of interest: J. Emmerich has nothing to disclose.
Conflict of interest: G. Le Gal has acted as a co-investigator in clinical trials for Portola Pharmaceuticals, Boehringer Ingelheim, Pfizer, Bristol-Myers Squibb, LEO Pharma, Daiichi Sankyo and Bayer, and has received honoraria, not taken as salary from Bayer, Pfizer, LEO Pharma, Sanofi and bioMérieux, outside the submitted work.
Conflict of interest: G. Meyer reports grants from Bayer and non-financial support from Leo Pharma, BMS-Pfizer and Daiichi Sankyo, outside the submitted work.
Conflict of interest: T. Wolfson has nothing to disclose.
Conflict of interest: A.C. Gamst has nothing to disclose.
Conflict of interest: R.E. Moore has nothing to disclose.
Conflict of interest: G.B. Gugiu has nothing to disclose.
Conflict of interest: T.A. Morris reports grants from Bayer Pharmaceuticals, during the conduct of the study; personal fees for consulting from Bayer Pharmaceuticals, and personal fees for acting as an expert witness from Faegre Baker Daniels, outside the submitted work.
Support statement: The study was funded by research grants from Bayer Pharmaceuticals, Fonds de Recherche en Santé Respiratoire, France, R01HL095089 and through collaboration with the UCSD DNA Sequencing Facility (supported, in part, by grant 2 P30 CA023100-23 from the NIH) for performing DNA sequence analysis; and the City of Hope Mass Spectrometry and Proteomics Core Facility (supported, in part, by grant CA033572 from the NIH) for performing fibrinogen mass spectrometry analysis. The authors had full discretion over the study design, performance and creation of the manuscript. Funding information for this article has been deposited with the Crossref Funder Registry.
- Received March 24, 2018.
- Accepted September 24, 2018.
- Copyright ©ERS 2018