# Inflammation and remodeling pathways and risk of cardiovascular events in patients with ischemic heart failure and reduced ejection fraction

Nicolas Girerd, John Cleland, Stefan Anker, William Byra, Carolyn Lam, David Lapolice, Mandeep Mehra, Dirk van Veldhuisen, Emmanuel Bresso, Zohra Lamiral, et al.

► To cite this version:

Nicolas Girerd, John Cleland, Stefan Anker, William Byra, Carolyn Lam, et al. Inflammation and remodeling pathways and risk of cardiovascular events in patients with ischemic heart failure and reduced ejection fraction. Scientific Reports, 2022, 12 (1), pp.8574. 10.1038/s41598-022-12385-0. hal-03675621

HAL Id: hal-03675621

https://hal.univ-lorraine.fr/hal-03675621

Submitted on 23 May 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

OPEN

# Inflammation and remodeling pathways and risk of cardiovascular events in patients with ischemic heart failure and reduced ejection fraction

Nicolas Girerd<fnref n="1" />, John Cleland², Stefan D. Anker³, William Byra⁴, Carolyn S. P. Lam⁵,
David Lapolice⁴, Mandeep R. Mehra⁶, Dirk J. van Veldhuisen⁷, Emmanuel Bresso<fnref n="1" />,
Zohra Lamiral<fnref n="1" />, Barry Greenberg⁸ & Faiez Zannad<fnref n="1" />✉

Patients with heart failure (HF) and coronary artery disease (CAD) have a high risk for cardiovascular (CV) events including HF hospitalization, stroke, myocardial infarction (MI) and sudden cardiac death (SCD). The present study evaluated associations of proteomic biomarkers with CV outcome in patients with CAD and HF with reduced ejection fraction (HFrEF), shortly after a worsening HF episode. We performed a case-control study within the COMMANDER HF international, double-blind, randomized placebo-controlled trial investigating the effects of the factor-Xa inhibitor rivaroxaban. Patients with the following first clinical events: HF hospitalization, SCD and the composite of MI or stroke were matched with corresponding controls for age, sex and study drug. Plasma concentrations of 276 proteins with known associations with CV and cardiometabolic mechanisms were analyzed. Results were corrected for multiple testing using false discovery rate (FDR). In 485 cases and 455 controls, 49 proteins were significantly associated with clinical events of which seven had an adjusted FDR < 0.001 (NT-proBNP, BNP, T-cell immunoglobulin and mucin domain containing 4 (TIMD4), fibroblast growth factor 23 (FGF-23), growth differentiation factor-15 (GDF-15), pulmonary surfactant-associated protein D (PSP-D) and Spondin-1 (SPON1)). No significant interactions were identified between the type of clinical event (MI/stroke, SCD or HF) and specific biomarkers (all interaction FDR > 0.20). When adding the biomarkers significantly associated with the above outcome to a clinical model (including NT-proBNP), the C-index increase was 0.057 (0.033-0.082), p < 0.0001 and the net reclassification index was 54.9 (42.5 to 67.3), p < 0.0001. In patients with HFrEF and CAD following HF hospitalization, we found that NT-proBNP, BNP, TIMD4, FGF-23, GDF-15, PSP-D and SPON1, biomarkers broadly associated with inflammation and remodeling mechanistic pathways, were strong but indiscriminate predictors of a variety of individual CV events.

Patients with heart failure (HF) and coronary artery disease (CAD) are at excess risk of cardiovascular (CV) events, including HF hospitalization, stroke and myocardial infarction (MI). Each of these individual events is associated with a three to fourfold increase in the subsequent risk of death, including sudden cardiac death (SCD), in patients with HF with reduced ejection fraction (HFrEF)<fnref n="1" />.

<footnote n="1">¹Université de Lorraine, Centre d'Investigation Clinique-Plurithématique Inserm CIC-P 1433, Inserm U1116, CHRU Nancy Brabois, F-CRIN INI-CRCT (Cardiovascular and Renal Clinical Trialists), Nancy, France. ²Robertson Centre for Biostatistics and Clinical Trials, University of Glasgow, Glasgow, Scotland. ³Department of Cardiology (CVK), and Berlin Institute of Health Center for Regenerative Therapies (BCRT), German Centre for Cardiovascular Research (DZHK) Partner Site Berlin, Charité Universitätsmedizin Berlin, Berlin, Germany. ⁴Janssen Research and Development, Raritan, NJ, USA. ⁵National Heart Centre Singapore, Duke-National University of Singapore, Singapore, Singapore. ⁶Brigham and Women's Hospital and Harvard Medical School, Boston, MA, USA. ⁷Department of Cardiology, University Medical Center Groningen, University of Groningen, Groningen, The Netherlands. ⁸Cardiology Division, Department of Medicine, University of California, La Jolla, San Diego, USA. ✉email: f.zannad@chru-nancy.fr</footnote>
Proteomic analysis that explores serum protein components can provide insight into the mechanism of the multiple biological processes underlying disease progression² and can be used to study mechanisms associated with various clinical events. By applying knowledge-based network analysis, associations among protein biomarkers (BMs) of pathophysiological importance can be discovered in patients with HF³ that would allow for the identification and understanding of new mechanistic pathways leading to clinical events and could ultimately provide insight into novel therapeutic strategies. Although proteomic marker discovery approaches have been shown to predict death⁴,⁵ and HF hospitalization⁴, their prognostic value in the specific setting of patients with CAD and HFrEF following HF hospitalization to predict other outcomes (such as stroke and MI) has been insufficiently explored. In addition, factors that determine whether these patients are more prone to develop a SCD, HF hospitalization, stroke or MI are unclear. Currently, few analyses exist regarding the prognostic value of proteomic profiles in the specific setting of worsening HF⁴. Determining the mechanistic bioprofiles of patients respectively at risk of SCD, HFH, stroke or MI, in the setting of worsening HF, may help more specific targeting of strategies to prevent such events in future trials. This enrichment may additionally help inform the design of acute HF trials to increase the likelihood of their conclusiveness, something that has been only rarely achieved over the past 40 years⁶.

The COMMANDER HF (Rivaroxaban in Patients with Heart Failure, Sinus Rhythm, and Coronary Disease) trial⁷ enrolled chronic CAD patients shortly following hospitalization for worsening HF. Blood samples were collected at the time of entry into the study and follow-up cardiovascular events were carefully recorded in the population as outlined in the study protocol. In the present study, we investigated whether analysis of baseline proteomics data helps predict subsequent CV events, including rehospitalization for HF, MI/stroke and SCD and whether proteomic analysis adds to the understanding of mechanisms underlying progression to these events following an admission for worsening HF.

**Methods**

**Population.** We performed a nested case-control analysis of the COMMANDER HF trial. Briefly, the COMMANDER HF trial included 5022 patients with worsening HF in the context of chronic HFrEF and underlying CAD⁷⁻⁹. All patients had a left ventricular ejection fraction (LVEF) ≤40%, CAD, and elevated plasma concentrations of natriuretic peptides (this latter criterion added during the course of the trial) and did not have atrial fibrillation. They were randomly assigned to receive rivaroxaban at a dose of 2.5 mg twice daily or placebo in addition to standard care. Over a median follow-up period of 21.1 months, the primary endpoint occurred in 626 (25.0%) of 2507 patients assigned to rivaroxaban and in 658 (26.2%) of 2515 patients assigned to placebo (hazard ratio, 0.94; 95% confidence interval [CI], 0.84 to 1.05; P = 0.27). This international trial was approved by all relevant committees in each country, and every patients provided written consent. All methods used followed relevant regulations.

**Outcomes.** The primary efficacy outcome of the COMMANDER HF trial was the composite of death from any cause, MI or stroke. Three clinical events were considered in the present analysis, namely HF hospitalization, SCD and the composite of MI/stroke, as ischemic events. In addition, we also analyzed the composite of all events at once, i.e. HF hospitalization, SCD and MI/stroke.

**Matching procedure.** Controls were matched for age, sex and study drug with cases that had one of the 4 clinical events described above (i.e. HF hospitalization, SCD and the composite of MI/ischemic stroke). Suitable patients could serve as controls for more than one case when no other suitable control could be identified. All of the cases that could be matched with a control and had available blood samples were considered in this analysis. Cases of a given event could experience another event during follow-up (e.g. a case patient for MI could subsequently experience stroke) (supplementary Table 1).

**Sample handling.** All sample shipments and sample data acquisition within the COMMANDER HF trial were carried out according to predefined standard operating procedures and material transfer agreements to maintain uniformity. The cases and controls were separately identified and selected. All participants' information was then removed and a randomly-sorted list of patient/sample IDs for each cohort was sent to the bioassay facility. The entire sampling handling/protein measurement procedure was carried out while maintaining the blind.

**Assays and studied biomarkers.** Baseline plasma samples were analyzed for protein biomarkers using the OLINK Proseek® Multiplex cardiovascular (CVD) II, CVD III, and cardiometabolic (CM) panels. These panels were selected for the well-balanced inclusion of proteins with established or likely associations with CV disease (the full information on these panels is available at: https://www.olink.com/resources-support/docum ent-download-center/). The assays use a proximity extension assay (PEA) technology, where 92 oligonucleotide- labeled antibody probe pairs per panel are allowed to bind to their respective targets in the sample in 96-well plate format. When binding to their correct targets, they give rise to new DNA amplicons with each ID-bar- coding their respective antigens. The amplicons are subsequently quantified using a Fluidigm BioMark™ HD real-time PCR platform. The OLINK platform provides log2-normalized protein expression (NPX) data. The Olink® quality control samples are considered as “flagged” if they deviate more than 0.3 NPX from the median of all samples in one of two control assays for incubation and detection. The LOD is defined by the three negative controls run on each plate and set to three standard deviations above the measured background.
The assays were performed "blinded" to case/control status with cases and controls randomly distributed across plates. The proteomic results were then merged with the baseline data, which included the case-control status, matching variables and the clinical risk factors. A total of 276 proteins were ultimately analyzed.

**Statistical analysis.** For baseline clinical characteristics, continuous variables are expressed as means ± standard deviation (SD) or as medians (interquartile range) if the distribution was skewed and categorical variables as frequencies and percentages. Baseline characteristics of participants were compared between cases and controls using Fisher's exact test for categorical variables and t-tests or non-parametric Wilcoxon tests for continuous variables, as appropriate.

To assess the association of biomarkers with outcomes, we performed logistic regression models adjusted for age, sex and factors significantly associated with case/control status (as reported in Table 1). The association with event risk was assessed for each biomarker individually (mono-marker models). These associations were first assessed for the composite outcome of hospitalization for worsening HF (WHF), sudden death and MI/stroke and subsequently for "individual" outcomes (either hospitalization for WHF, sudden death or MI/stroke). In addition, we searched for an interaction between biomarkers and the type of clinical event (either hospitalization for WHF, sudden death or MI/stroke). To achieve this, an interaction term was introduced between each biomarker and the outcome subgroup; this interaction term tests whether the association of biomarkers with outcome is homogenous across outcome types. In an additional analysis, the models were further adjusted for NT-proBNP. No imputation was performed during the analysis process.

In all 3 cohorts, analysis was corrected for multiple testing using a false discovery rate (FDR) of 5%, applying the Benjamini-Hochberg procedure¹⁰.

Since proteins were measured using NPX (Normalized Protein eXpression) values (details regarding NPX values can be found at: https://www.olink.com/question/what-is-npx/), the odds ratio for each protein estimates the increase in the odds of HF associated with a doubling in protein concentration.

In a supplementary analysis, we used a stabilized inverse probability of treatment weighting method (sIPTW) based on propensity score¹¹. Weights were obtained from a logistic regression with outcome as dependent variable and clinical adjustment variables as explanatory variables (sex, age, study drug and significant factors from Table 1 (eGFR, LVEF, NYHA class, MI, diabetes, MRA and digoxin)).

To determine the added predictive value of biomarkers on top on clinical predictors, we assessed the C-index of a baseline model including all adjustment variables (matching variables [sex, age, study drug] and variables significantly different across groups in Table 1 [namely eGFR, LVEF, NYHA class, MI, diabetes, MRA and digoxin]) and assessed the increase in C-index when adding biomarkers to the model. In addition, the increase in discriminative value of the addition of biomarkers on top of the aforementioned clinical model was assessed using continuous net reclassification improvement (NRI)¹².

**Complex network analysis.** Pathway overrepresentation was calculated for proteins selected by the statistical analysis. Pathway-annotation proteins were extracted from the FHF-GKBox¹³. In the FHF-GKBox, protein-pathway relationships were extracted from Reactome (version 69)¹⁴. The entire FHF-GKBox was considered as annotation background for a Fisher's exact test and obtained p-values were corrected for multiple testing using a FDR of 5% by applying the Benjamini-Hochberg method. Only pathways with a FDR lower than 5% and selected proteins linked to these pathways were considered to construct the graph. The network was displayed using Cytoscape (version 3.7.2)¹⁵.

## Results

**Patient characteristics.** The analysis included 485 cases and 455 controls (as a subject could be a control for more than one case). Mean age was 67 years and 24% of patients were female (Table 1). Patients with clinical events had significantly lower eGFR and LVEF than controls (Table 1). They also had a higher probability of MI history and diabetes and were more likely to be treated with MRA or digoxin.

Overall, the characteristics according to the presence/absence of specific events (MI/stroke, SCD or HFH) were similar to those observed in the complete case/control study (Supplementary Table 2).

**Association of protein biomarkers with total cardiovascular clinical events.** When correcting for the multiplicity of tests and adjusting for age, sex, study drug and relevant factors identified in univariable analysis, 49 proteins were significantly associated with CV clinical events at the p < 0.05 FDR threshold. The 7 proteins associated with outcome with an FDR < 0.001 were NT-proBNP, BNP, T-cell immunoglobulin and mucin domain containing 4 (TIMD4), fibroblast growth factor 23 (FGF-23), growth differentiation factor-15 (GDF-15), pulmonary surfactant-associated protein D (PSP-D) and Spondin-1 (SPON1), all positively associated with a increased risk of CV events. These biomarkers, except for BNP and NT-proBNP, were poorly to moderately correlated with each other (Fig. 1).

When further adjusting for NT-proBNP, 11 proteins were still significantly associated with the composite CV outcome, namely TIMD4, CCL18, IL1RL2, CXCL1, FCGR2A, PSP-D, UMOD, CRTAC1, CCL24, TFPI and TGM2 (Table 2).

A similar pattern of association was observed when using sIPTW (supplementary Table 3).

No significant interactions with individual clinical events (MI/stroke, SCD or HF hospitalization) were identified (all interaction FDR > 0.20, supplementary Table 4). Results remained unchanged when considering sIPTW analyses (supplementary Table 5).

No significant interactions were identified between the time to event and the prognostic value of biomarkers (all interaction p-values > 0.10, data not shown).
<table><caption>**Table 1.** Patient characteristics according to MI/stroke, SCD, HF rehospitalization overall in the sex and age-matched COMMANDER HF population. Significant values Queryare given in bold.</caption><thead><tr><th rowspan="3">Characteristic</th><th colspan="2">All events</th><th rowspan="3">SMD (%)</th><th rowspan="3">p</th></tr><tr><th>Controls</th><th>Cases</th></tr><tr><th>(n = 455)</th><th>(n = 485)</th></tr></thead><tbody><tr><td>Age (yrs)</td><td>67.4 ± 10.0</td><td>67.3 ± 10.3</td><td>1.40</td><td>0.83</td></tr><tr><td>Study drug</td><td>210 (46.2%)</td><td>230 (47.3%)</td><td>2.30</td><td>0.74</td></tr><tr><td>Female sex</td><td>110 (24.2%)</td><td>120 (24.7%)</td><td>1.20</td><td>0.88</td></tr><tr><td><strong>Race</strong></td><td></td><td></td><td>14.60</td><td>0.75</td></tr><tr><td>White</td><td>405 (89.0%)</td><td>433 (89.1%)</td><td></td><td></td></tr><tr><td>Black</td><td>2 (0.4%)</td><td>5 (1.0%)</td><td></td><td></td></tr><tr><td>Asian</td><td>42 (9.2%)</td><td>41 (8.4%)</td><td></td><td></td></tr><tr><td>Other</td><td>6 (1.3%)</td><td>7 (1.4%)</td><td></td><td></td></tr><tr><td><strong>Region</strong></td><td></td><td></td><td>21.50</td><td>0.13</td></tr><tr><td>Eastern Europe</td><td>315 (69.2%)</td><td>332 (68.3%)</td><td></td><td></td></tr><tr><td>North America</td><td>5 (1.1%)</td><td>9 (1.9%)</td><td></td><td></td></tr><tr><td>Asia Pacific</td><td>42 (9.2%)</td><td>41 (8.4%)</td><td></td><td></td></tr><tr><td>Latin America</td><td>50 (11.0%)</td><td>38 (7.8%)</td><td></td><td></td></tr><tr><td>Western Europe</td><td>43 (9.5%)</td><td>66 (13.6%)</td><td></td><td></td></tr><tr><td>BMI (kg/m²)</td><td>27.6 ± 5.0</td><td>27.9 ± 5.1</td><td>6.30</td><td>0.33</td></tr><tr><td>eGFR (ml/min/1.73 m²)</td><td><strong>70.1 ± 22.9</strong></td><td><strong>64.7 ± 23.0</strong></td><td><strong>23.80</strong></td><td><strong>0.0003</strong></td></tr><tr><td><strong>eGFR (ml/min/1.73 m²)</strong></td><td></td><td></td><td><strong>23.60</strong></td><td><strong>0.0051</strong></td></tr><tr><td>&lt; 30 ml/min/1.73 m²</td><td>11 (2.4%)</td><td>25 (5.1%)</td><td></td><td></td></tr><tr><td>30 to &lt; 60 ml/min/1.73 m²</td><td>157 (34.5%)</td><td>199 (40.9%)</td><td></td><td></td></tr><tr><td>60 to &lt; 90 ml/min/1.73 m²</td><td>204 (44.8%)</td><td>201 (41.4%)</td><td></td><td></td></tr><tr><td>≥ 90 ml/min/1.73 m²</td><td>83 (18.2%)</td><td>61 (12.6%)</td><td></td><td></td></tr><tr><td>BNP level (pg/ml)</td><td>875.2 ± 787.4</td><td>950.3 ± 727.8</td><td>9.90</td><td>0.53</td></tr><tr><td>Log2 BNP (pg/ml)</td><td>6.4 ± 0.8</td><td>6.6 ± 0.7</td><td>21.70</td><td>0.17</td></tr><tr><td>BNP rank (pg/ml)</td><td></td><td></td><td></td><td></td></tr><tr><td>NT-proBNP (pg/ml)</td><td>5869 ± 8707</td><td>5568 ± 7310</td><td>3.70</td><td>0.66</td></tr><tr><td>Log2 NT-proBNP (pg/ml)</td><td></td><td></td><td></td><td></td></tr><tr><td>NT-proBNP rank</td><td>265.8 ± 160.5</td><td>281.2 ± 154.9</td><td>9.80</td><td>0.26</td></tr><tr><td>D-dimer level (mg/l)</td><td>637.2 ± 1012.5</td><td>732.3 ± 1027.4</td><td>9.30</td><td>0.16</td></tr><tr><td>Log2 D-dimer</td><td>6.0 ± 0.9</td><td>6.2 ± 0.8</td><td>24.40</td><td>0.0002</td></tr><tr><td>D-dimer rank</td><td></td><td></td><td></td><td></td></tr><tr><td>Ejection fraction (%)</td><td><strong>33.1 ± 6.7</strong></td><td><strong>31.1 ± 7.3</strong></td><td><strong>28.40</strong></td><td><strong>&lt; 0.0001</strong></td></tr><tr><td>Ejection fraction rank</td><td><strong>511.5 ± 267.6</strong></td><td><strong>433.0 ± 269.1</strong></td><td><strong>29.30</strong></td><td><strong>&lt; 0.0001</strong></td></tr><tr><td><strong>NYHA</strong></td><td></td><td></td><td><strong>12.40</strong></td><td><strong>0.031</strong></td></tr><tr><td>I/II</td><td><strong>238 (52.4%)</strong></td><td><strong>226 (46.5%)</strong></td><td></td><td></td></tr><tr><td>III</td><td><strong>214 (47.0%)</strong></td><td><strong>246 (50.6%)</strong></td><td></td><td></td></tr><tr><td>IV</td><td><strong>3 (0.7%)</strong></td><td><strong>14 (2.9%)</strong></td><td></td><td></td></tr><tr><td>Myocardial infarction</td><td><strong>331 (72.7%)</strong></td><td><strong>385 (79.2%)</strong></td><td><strong>15.20</strong></td><td><strong>0.022</strong></td></tr><tr><td>Stroke</td><td>44 (9.7%)</td><td>58 (11.9%)</td><td>7.30</td><td>0.29</td></tr><tr><td>Diabetes</td><td><strong>163 (35.8%)</strong></td><td><strong>223 (45.9%)</strong></td><td><strong>20.60</strong></td><td><strong>0.002</strong></td></tr><tr><td>Hypertension</td><td>339 (74.5%)</td><td>375 (77.2%)</td><td>6.20</td><td>0.36</td></tr><tr><td>ACEI or ARB</td><td>431 (94.7%)</td><td>452 (93.0%)</td><td>7.20</td><td>0.28</td></tr><tr><td>Beta blockers</td><td>424 (93.2%)</td><td>445 (91.6%)</td><td>6.10</td><td>0.39</td></tr><tr><td>MRA</td><td><strong>327 (71.9%)</strong></td><td><strong>383 (78.8%)</strong></td><td><strong>16.10</strong></td><td><strong>0.015</strong></td></tr><tr><td>Digoxin</td><td><strong>20 (4.4%)</strong></td><td><strong>43 (8.8%)</strong></td><td><strong>18.00</strong></td><td><strong>0.009</strong></td></tr><tr><td>Aspirin</td><td>426 (93.6%)</td><td>450 (92.6%)</td><td>4.10</td><td>0.61</td></tr></tbody></table>

**Prediction of clinical events using a multi-marker approach.** The C-index for the clinical prediction of the composite CV events when using clinical variables was 0.66 (0.62–0.69), p<0.0001 (Fig. 2). NT-proBNP significantly increased the C-index and provided a significant NRI increase of 46.2 (33.7 to 58.7), p<0.0001 (Fig. 2). When further adding the biomarkers significantly associated with the aforementioned outcome, the C-index further rose to 0.760 (0.730–0.790), p<0.0001, resulting in a C-index increase related to these biomarkers of 0.057 (0.033–0.082), p<0.0001 and a NRI of 54.9 (42.5 to 67.3), p<0.0001 (Fig. 2).
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Heat map of biomarkers identified as significantly associated with events (as in Table 2).</figcaption></figure>

**Complex network analysis.** The 49 biomarker proteins associated with the sentinel events highlighted 4 overrepresented pathways (Fig. 3): “post-translational protein phosphorylation”, “regulation of IGF transport and uptake by IGF binding proteins”, “signaling by interleukins” and “extracellular matrix organization”. “Regulation of IGF transport and uptake by IGF binding proteins” was the pathway associated with the highest number of proteins (N=9) while the 3 other pathways were each associated with 7 proteins. These pathways were linked to 15 significant biomarkers among which six were associated with one pathway and four were associated with two pathways. Four proteins (Matrix metalloproteinase-2 [MMP-2], Interleukin 6 [IL6], Osteopontin [OPN] and Tenascin [TNC]) were related to three pathways. Metalloproteinase inhibitor 1 (TIMP-1) was the only protein associated with all of the pathways.

**Discussion**

We identified 49 proteins significantly associated with clinical CV events in patients with HFrEF and underlying CAD following a worsening HF episode. Unsurprisingly, natriuretic peptides were among the best predictors of subsequent clinical CV events, along with TIMD4, FGF-23, GDF-15, PSP-D and SPON1, biomarkers associated with “post-translational protein phosphorylation”, “regulation of IGF transport and uptake by IGF binding proteins”, “signaling by interleukins” and “extracellular matrix organization” in network analysis, and in which TIMP-1 appears as a central node. While identification of these pathways may indicate that they play a role in
<table><caption>**Table 2.** Significant adjusted associations of protein biomarkers with CV events. Significant values are given in bold. <sup>a</sup>Clinical adjustment for sex, age, study drug and significant factors from Table 1 (eGFR, LVEF, NYHA class, MI, diabetes, MRA and digoxin).</caption><thead><tr><th rowspan="2">Biomarker</th><th colspan="2">Models adjusted for clinical variables<sup>a</sup></th><th rowspan="2">Models adjusted for clinical variables<sup>a</sup> and NT-proBNP</th></tr><tr><th>OR (95% CI)</th><th>FDR</th></tr></thead><tbody><tr><td>NT-proBNP</td><td>1.376 (1.259-1.504)</td><td>&lt;0.0001</td><td>-</td></tr><tr><td>BNP</td><td>1.292 (1.202-1.390)</td><td>&lt;0.0001</td><td>-</td></tr><tr><td>TIMD4</td><td>1.723 (1.367-2.171)</td><td>0.0004</td><td>1.484 (1.169-1.884)</td></tr><tr><td>FGF-23</td><td>1.265 (1.141-1.402)</td><td>0.0005</td><td>1.094 (0.978-1.224)</td></tr><tr><td>GDF-15</td><td>1.514 (1.257-1.824)</td><td>0.0005</td><td>1.185 (0.966-1.453)</td></tr><tr><td>PSP-D</td><td>1.438 (1.222-1.691)</td><td>0.0005</td><td>1.237 (1.042-1.469)</td></tr><tr><td>SPON1</td><td>1.921 (1.432-2.577)</td><td>0.0005</td><td>1.135 (0.809-1.591)</td></tr><tr><td>VEGFD</td><td>1.821 (1.364-2.431)</td><td>0.0016</td><td>1.146 (0.850-1.545)</td></tr><tr><td>TNC</td><td>1.440 (1.196-1.735)</td><td>0.0037</td><td>1.145 (0.933-1.406)</td></tr><tr><td>AOC3</td><td>1.712 (1.283-2.284)</td><td>0.0064</td><td>1.178 (0.863-1.609)</td></tr><tr><td>IGFBP-7</td><td>1.540 (1.220-1.943)</td><td>0.0064</td><td>1.029 (0.791-1.339)</td></tr><tr><td>IL-1RT1</td><td>1.786 (1.310-2.436)</td><td>0.0064</td><td>1.253 (0.897-1.750)</td></tr><tr><td>TFF3</td><td>1.519 (1.211-1.906)</td><td>0.0064</td><td>1.170 (0.917-1.493)</td></tr><tr><td>TIMP1</td><td>1.602 (1.236-2.077)</td><td>0.0073</td><td>1.114 (0.839-1.480)</td></tr><tr><td>U-PAR</td><td>1.529 (1.200-1.949)</td><td>0.0109</td><td>1.111 (0.854-1.446)</td></tr><tr><td>ST2</td><td>1.350 (1.136-1.603)</td><td>0.0111</td><td>1.089 (0.902-1.315)</td></tr><tr><td>TR</td><td>1.349 (1.132-1.607)</td><td>0.0132</td><td>1.105 (0.916-1.332)</td></tr><tr><td>OPG</td><td>1.653 (1.227-2.226)</td><td>0.0143</td><td>1.111 (0.802-1.539)</td></tr><tr><td>TFPI</td><td>1.658 (1.222-2.251)</td><td>0.0171</td><td>1.383 (1.007-1.898)</td></tr><tr><td>FCGR2A</td><td>1.406 (1.143-1.730)</td><td>0.0177</td><td>1.288 (1.041-1.594)</td></tr><tr><td>COL18A1</td><td>1.533 (1.175-2.000)</td><td>0.0217</td><td>1.155 (0.871-1.530)</td></tr><tr><td>CXCL1</td><td>1.207 (1.071-1.360)</td><td>0.0247</td><td>1.159 (1.025-1.310)</td></tr><tr><td>IGLC2</td><td>1.431 (1.134-1.805)</td><td>0.0302</td><td>1.204 (0.945-1.534)</td></tr><tr><td>IGFBP-2</td><td>1.313 (1.098-1.570)</td><td>0.0332</td><td>0.900 (0.725-1.117)</td></tr><tr><td>IL6</td><td>1.168 (1.054-1.295)</td><td>0.0332</td><td>1.097 (0.990-1.216)</td></tr><tr><td>TGM2</td><td>1.223 (1.070-1.397)</td><td>0.0332</td><td>1.150 (1.003-1.318)</td></tr><tr><td>ACE2</td><td>1.308 (1.087-1.575)</td><td>0.0391</td><td>1.045 (0.861-1.269)</td></tr><tr><td>CCL18</td><td>1.280 (1.080-1.516)</td><td>0.0391</td><td>1.266 (1.063-1.508)</td></tr><tr><td>CD93</td><td>1.574 (1.154-2.148)</td><td>0.0391</td><td>0.923 (0.648-1.315)</td></tr><tr><td>IL1RL2</td><td>0.671 (0.509-0.884)</td><td>0.0391</td><td>0.674 (0.506-0.897)</td></tr><tr><td>MMP-2</td><td>1.475 (1.129-1.926)</td><td>0.0391</td><td>0.732 (0.522-1.028)</td></tr><tr><td>PGLYRP1</td><td>1.334 (1.096-1.623)</td><td>0.0391</td><td>1.222 (0.997-1.496)</td></tr><tr><td>CCL14</td><td>1.463 (1.123-1.906)</td><td>0.0391</td><td>1.078 (0.814-1.430)</td></tr><tr><td>Notch 3</td><td>1.457 (1.121-1.893)</td><td>0.0391</td><td>0.863 (0.635-1.172)</td></tr><tr><td>OSMR</td><td>2.187 (1.267-3.775)</td><td>0.0391</td><td>1.124 (0.622-2.032)</td></tr><tr><td>CCL24</td><td>1.217 (1.059-1.399)</td><td>0.0408</td><td>1.174 (1.016-1.355)</td></tr><tr><td>TIMP4</td><td>1.416 (1.108-1.808)</td><td>0.0408</td><td>1.080 (0.832-1.402)</td></tr><tr><td>TRAIL-R2</td><td>1.404 (1.104-1.785)</td><td>0.0408</td><td>1.047 (0.821-1.336)</td></tr><tr><td>UMOD</td><td>0.525 (0.333-0.830)</td><td>0.0408</td><td>0.573 (0.357-0.919)</td></tr><tr><td>CD163</td><td>1.413 (1.097-1.820)</td><td>0.0458</td><td>1.153 (0.886-1.502)</td></tr><tr><td>ICAM1</td><td>1.515 (1.119-2.052)</td><td>0.0458</td><td>1.229 (0.897-1.684)</td></tr><tr><td>IL-27</td><td>1.473 (1.110-1.956)</td><td>0.0458</td><td>1.088 (0.816-1.452)</td></tr><tr><td>OPN</td><td>1.286 (1.070-1.547)</td><td>0.0458</td><td>0.935 (0.757-1.155)</td></tr><tr><td>RARRES2</td><td>1.680 (1.154-2.446)</td><td>0.0458</td><td>1.458 (0.990-2.148)</td></tr><tr><td>SERPINA5</td><td>0.708 (0.551-0.910)</td><td>0.0458</td><td>0.951 (0.726-1.247)</td></tr><tr><td>LTBP2</td><td>1.468 (1.105-1.949)</td><td>0.0483</td><td>0.831 (0.595-1.161)</td></tr><tr><td>CRTAC1</td><td>0.756 (0.613-0.932)</td><td>0.0498</td><td>0.780 (0.627-0.970)</td></tr><tr><td>GH</td><td>1.103 (1.025-1.188)</td><td>0.0498</td><td>1.020 (0.942-1.104)</td></tr><tr><td>IGFBP-1</td><td>1.146 (1.035-1.270)</td><td>0.0498</td><td>1.003 (0.896-1.122)</td></tr></tbody></table>
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Changes in the prediction of CV events using clinical and biomarker variables. The clinical variables included in the baseline model were the following: age, sex, study treatment, eGFR, LVEF, NYHA class, history of MI, history of diabetes, MRA and digoxin.</figcaption></figure>

the subsequent clinical events that were experienced by patients in the study, we could not identify a significant interaction between biomarkers and the clinical event per se (i.e. either MI/stroke, SCD or HF hospitalization) in these patients.

**Protein biomarkers are homogenously associated with different outcomes.** Our analysis intended to identify different biological profiles underlying distinct CV clinical events. Yet, none of the biomarkers identified herein were differentially associated with MI/stroke, SCD or HF hospitalization, respectively, when assessing formal statistical interactions while accounting for multiplicity of tests (all p values for interaction > 0.20). This result is consistent with the possibility that there are unifying underling mechanisms to all clinical events in a population that is comprised of patients with CAD and HFrEF who are admitted for worsening HF. The pathways related to the biomarkers identified were "post-translational protein phosphorylation", "regulation of IGF transport and uptake by IGF binding proteins", "signaling by interleukins" and "extracellular matrix organization". These pathways relate to the broad concept of inflammation and remodeling, both of which are central to HF development and worsening but are also very much associated with stroke and MI biology. Our results could also suggest that there is some overlap within the pathways involved in these events. For instance, MI is a significant cause of HF initiation and worsening. In addition, stroke is very much associated with MI, and both MI and stroke are frequently associated with CV death in patients with HF. Differentiating distinct biological profiles for these events may hence be challenging due to the strong mechanistic and prognostic interplay of these clinical events with one another.

**The value of natriuretic peptides to predict a range of CV outcomes in patients with HFrEF admitted for worsening HF.** Our study provides proteomic results for a large population of patients recently hospitalized for HF, a population that has been less well studied comparatively to ambulatory HF. Natriuretic peptides were unsurprisingly among the best predictors of subsequent clinical events in the specific context of HF hospitalization. However, our analyses emphasize that natriuretic peptides also predict a range of CV outcomes, including ischemic events (MI/stroke) and SCD, suggesting that they may be a pleiotropic/ omnipotent predictor of events in patients with HFrEF following HF hospitalization. Thus, elevated NT-proBNP should not be perceived solely as a risk for HF related outcomes, but rather as an indicator of overall risk of CV events¹⁶ (along with LVEF¹⁷ and eGFR, as emphasized in Table 2).

**TIMD4, FGF-23, GDF-15, PSP-D and SPON1 as key predictors of outcome in patients with HFrEF admitted for worsening HF.** The next best predictors of events in our study were TIMD4, FGF-23, GDF-15, PSP-D and SPON1. Importantly, none of these 5 proteins were highlighted in the enriched network of upregulated proteins in ischemic HF relative to non-ischemic HF in the BIOSTAT cohort¹⁸. This striking dif-
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Overrepresented pathways (green triangles) linked to their significant protein biomarkers (red circles). IGFBP-2: Insulin-like growth factor-binding protein 2; IGFBP-1: Insulin-like growth factor-binding protein 1; IGFBP-7: Insulin-like growth factor-binding protein 7; FGF-23: Fibroblast growth factor 23; IL27: Interleukin-27 subunit alpha; IL-1RT1: Interleukin-1 receptor type 1; CXCL1: Growth-regulated alpha protein; IL6: Interleukin-6; LTBP2: Latent-transforming growth factor beta-binding protein 2; ICAM1: Intercellular adhesion molecule 1; TNC: Tenascin; OPN: Osteopontin; TIMP-1: Metalloproteinase inhibitor 1; MMP-2: Matrix metalloproteinase-2; COL18A1: Collagen alpha-1(XVIII) chain.</figcaption></figure>

ference may suggest that proteins related to the HF etiology might not be related to prognosis in patients with ischemic HF.

There is only sparse evidence regarding TIMD4 in CV disease. TIMD4 haplotypes have been reported to be associated with CAD and stroke as well as statin therapy in a case-control study performed in Chinese populations¹⁹. In our analysis, despite the non-significant statistical interaction, it should be acknowledged that the point estimate of associations were numerically higher for MI/stroke (OR = 2.832, (1.701–4.713), p = 0.0086) than for the other CV outcomes (OR for SCD = 1.853, (1.240–2.768), p = 0.0769; OR for HFH = 1.304, (0.928–1.831), p = 0.8652; supplementary Table 3).

FGF-23 has been shown to predict death and/or HFH in both ambulatory and hospitalized patients⁴,²⁰. Of note, in the BIOSTAT cohort, in which both ambulatory and hospitalized patients were included, FGF-23 was similarly retained in multivariable models in inpatients (HR = 1.13 (1.02, 1.24)) and outpatients (HR = 1.32 (1.14, 1.53)) with HF⁴. We further show herein that, in patients hospitalized for HF in the context of CAD, FGF-23 is also a key predictor of CV outcome (even when considering non-HF outcomes). Taken together, these data suggest that FGF-23, alongside with NT-proBNP, is a potent biomarker for overall CV risk prediction in HF.

GDF-15 is associated with death and/or HFH in hospitalized⁴,²¹ patients with HF. This biomarker has also been shown to decrease during hospital stay in concomitance with improvement of signs and symptoms²¹, which suggests that it may be partially related to congestion. In the BIOSTAT cohort, GDF-15 was one of the main biomarkers differentiating inpatient and outpatient settings, with GDF-15 being higher in inpatients⁴. However, these 2 previous cohorts reported contrasting associations: Lourenço et al. reported that, in multivariable analysis, patients with either elevated admission or discharge GDF-15 had a significant twofold increase in 1-year death risk²¹. In contrast, in the BIOSTAT cohort, higher GDF-15 levels were associated with lower rates of HF hospitalization in the inpatient subgroup (HR = 0.80 (0.67, 0.95), p = 0.01)⁴, whereas GDF-15 was not retained as a significant predictor in the outpatient subgroup. Our results support the findings of Lourenço et al. and suggest that GDF-15 predicts a range of CV events such as MI/stroke and SCD (OR = 1.736, (1.265–2.380), p = 0.0345; supplementary Table 3) in the setting of hospitalized HF with CAD history.
Limited evidence exists regarding PSP-D in HF. Surfactant-derived proteins have been emphasized as markers of alveolar membrane damage in HF²². In addition, PSP-D trajectories over time have also been shown to be associated with the composite outcome of HFH, cardiac death, LVAD-placement and heart transplantation in a cohort of chronic HF patients²³. Our results further suggest that PSP-D is important in more acute settings, which is in line with the previously reported involvement of PSP-D in alveolar damage²², an event likely to occur in worsening HF.

SPON-1 has been reported to be associated with the risk of HF hospitalization in chronic HF with CKD²⁴. Increase in plasma levels of SPON1 has also been shown to precede adverse HF-related events and CV death in patients with CHF²⁵ in a serial plasma measurements study using joint modeling. These prior data are in keeping with our results showing the strong prognostic value of SPON-1 observed in hospitalized patients with often some level of concomitant worsening in renal function.

Importantly, all of these proteins, along with the other proteins significantly associated with outcome as highlighted in the network analysis, relate to the broad concept of inflammation and remodeling.

**TIMP-1 as a central node in pathways associated with higher risk in worsening HF with CAD.** In our network analysis, TIMP-1 was associated with all pathways identified as overrepresented from the 49 biomarker proteins (Fig. 2). This factor is known to be associated with “post-translational protein phosphorylation”, “regulation of IGF transport and uptake by IGF binding proteins”, “signaling by interleukins” and “extracellular matrix organization”. Importantly, TIMP-1 regulates various metalloprotease activity²⁶ including MMP and a disintegrin as well as metalloproteinases (ADAMs). TIMP-1 is consequently closely related to cardiovascular remodeling²⁶. It has specifically been reported to modulate inflammation and extracellular matrix fibrosis following injury²⁷, an aspect that could be particularly relevant in the setting of acute HF and/or ischemic HF. Importantly, this biomarker was not measured in previous studies assessing proteomics in HF (such as the BIOSTAT cohort⁴). Polymorphisms in the *TIMP-1* (and *TIMP-2*) genes have been shown to be associated with poor outcome in chronic HF²⁸. TIMP-1 is elevated in HF patients when compared to healthy populations²⁹, and higher blood levels of TIMP-1 are negatively correlated with peak VO₂ in HF patients²⁹. A recent meta-analysis of animal models showed that TIMP-1 is consistently induced in hemodynamic models, and is associated with cardiac fibrosis³⁰. The central place of TIMP-1 in the current complex network is in line with the central and unifying place of fibrosis in progressive HF deterioration.

## Limitations

Certain limitations should be pointed out in the present study.

First, the case-control study was designed to identify protein biomarkers for CV risk, although cannot provide any information regarding causality. This case-control design was chosen because of logistical/cost constrains; a cohort design would have decreased the risk of bias.

In addition, the proteomics assay did not provide standard concentration units, making comparisons with clinically applied cutoffs difficult. Still, the Olink standard procedures ensure a good correlation with standard measurement methodologies, and in this mechanistic screening study, there was no intention to measure specific concentrations of biomarkers, but only to investigate the overexpression of proteins associated with target mechanistic pathways.

Pre-selection is an inherent limitation of our proteomic analysis. Indeed, the choice of panels was based on proteins that have previously been associated with CV disease. Therefore an opportunity to discover additional new proteins not covered by the selected CVD panels and of potential prognostic importance could have possibly been missed. Similarly, our approach was focused on proteins, thus potentially missing integrated clinical/biological phenotypes relevant to the prediction of specific CV outcomes.

Interaction analyses are known to lack statistical power. Consequently, our analysis may have missed “true” differential associations of biomarkers with specific outcomes given the relatively modest sample size of the study. However, genuinely strong determinants of distinct types of events are unlikely to have been missed.

Our analysis was limited by its case-control nature and the central place of TIMP1 should be further assessed in cohort studies.

The inclusion of patients could occur up to several weeks after patient admission for worsening HF. This does introduce some variability across measurements and may have some impact on the preciseness of the associations with the outcomes identified herein.

Importantly, this study comprised patients with HFrEF. Given the differences in biological phenotypes observed in HF according to different levels of LVEF³¹, these results likely do not apply to patients with HF with preserved LVEF.

Due to the moderate sample size and the matched nature of the dataset, we were unable to provide an additional interaction analysis with patient characteristics, e.g. LVEF.

Our pathway analysis is also limited by the number of proteins identified as associated with outcome. Additional biological analyses including RNA sequencing tools could help in better characterizing involved biological pathways.

A replication phase involving a prospective validation of these biomarkers in other populations is required to improve the external validation of these results. Finally, the implementation of these results in the future would likely necessitate a broader access to proteomic data and web-based software able to provide a risk calculation based on proteomic data. Numerous steps are still needed to implement these approaches in clinical practice.

Lastly, only CV outcomes were considered in this analysis. Consideration to non-CV events (e.g. cancer) might have helped in identifying CV-specific biological phenotypes (even if not specific to a particular CV event). In addition, patients with a given event could experience another event during follow-up (i.e. there was
some overlap between events) which may have hampered our ability to identify specific profiles of biomarkers related to a given clinical event. Nevertheless, such event overlap is actually what patients experience in routine practice, and this overlap is consequently inherent to cardiovascular medicine.

Conclusion

We identified 49 proteins significantly associated with clinical events in patients with HFrEF and underlying CAD following a worsening HF episode. In this setting, the added value of these biomarkers for risk prediction was of similar magnitude to that derived from natriuretic peptides. Natriuretic peptides however appear as pleiotropic and potent predictors of events as they are similarly associated with HF-related events, ischemic (MI/stroke) events and SCD. In addition, TIMD4, FGF-23, GDF-15, PSP-D and SPON1 appear as key predictors of outcome in patients with HFrEF admitted for worsening HF, and position TIMP-1 as a central pathophysiological node in this setting.

Received: 26 August 2021; Accepted: 21 March 2022

Published online: 20 May 2022

References

1. Kristensen, S. L. et al. Comparison of outcomes after hospitalization for worsening heart failure, myocardial infarction, and stroke in patients with heart failure and reduced and preserved ejection fraction. *Eur. J. Heart Fail.* **17**, 169–176 (2015).

2. Tyers, M. & Mann, M. From genomics to proteomics. *Nature* **422**, 193–197 (2003).

3. Ferreira, J. P. et al. Proteomic bioprofiles and mechanistic pathways of progression to heart failure. *Circ. Heart Fail.* **12**, e005897 (2019).

4. Davison, B. A. et al. Is acute heart failure a distinctive disorder? An analysis from BIOSTAT-CHF. *Eur. J. Heart Fail.* **23**, 43–57 (2021).

5. Gui, H. et al. Plasma proteomic profile predicts survival in heart failure with reduced ejection fraction. *Circ. Genom. Precis. Med.* **14**, e003140 (2021).

6. Arrigo, M. et al. Acute heart failure. *Nat. Rev. Dis. Primers* **6**, 16 (2020).

7. Zannad, F. et al. Rivaroxaban in patients with heart failure, sinus rhythm, and coronary disease. *N. Engl. J. Med.* **379**, 1332–1342 (2018).

8. Pieske, B. et al. Baseline features of the VICTORIA (Vericiguat Global Study in Subjects with Heart Failure with Reduced Ejection Fraction) trial. *Eur. J. Heart Fail.* **21**, 1596–1604 (2019).

9. Zannad, F. et al. Rationale and design of a randomized, double-blind, event-driven, multicentre study comparing the efficacy and safety of oral rivaroxaban with placebo for reducing the risk of death, myocardial infarction or stroke in subjects with heart failure and significant coronary artery disease following an exacerbation of heart failure: The COMMANDER HF trial. *Eur. J. Heart Fail.* **17**, 735–742 (2015).

10. Green, G. H. & Diggle, P. J. On the operational characteristics of the Benjamini and Hochberg False Discovery Rate procedure. *Stat. Appl. Gen. Mol. Biol.* **6**, 27 (2007).

11. Austin, P. C. & Stuart, E. A. Moving towards best practice when using inverse probability of treatment weighting (IPTW) using the propensity score to estimate causal treatment effects in observational studies. *Stat. Med.* **34**, 3661–3679 (2015).

12. Uno, H., Tian, L., Cai, T., Kohane, I. S. & Wei, L. J. A unified inference procedure for a class of measures to assess improvement in risk prediction systems with survival data. *Stat. Med.* **32**, 2430–2442 (2013).

13. Girerd, N., Bresso, E., Devignes, M. D. & Rossignol, P. Insulin-like growth factor binding protein 2: A prognostic biomarker for heart failure hardly redundant with natriuretic peptides. *Int. J. Cardiol.* **300**, 252–254 (2020).

14. Fabregat, A. et al. Reactome graph database: Efficient access to complex pathway data. *PLoS Comput. Biol.* **14**, e1005968 (2018).

15. Shannon, P. et al. Cytoscape: A software environment for integrated models of biomolecular interaction networks. *Genome Res.* **13**, 2498–2504 (2003).

16. de Lemos, J. A. et al. The prognostic value of B-type natriuretic peptide in patients with acute coronary syndromes. *N. Engl. J. Med.* **345**, 1014–1021 (2001).

17. Hall, T. S. et al. Left ventricular ejection fraction and adjudicated, cause-specific hospitalizations after myocardial infarction complicated by heart failure or left ventricular dysfunction. *Am. Heart J.* **215**, 83–90 (2019).

18. Sama, I. E. et al. A network analysis to identify pathophysiological pathways distinguishing ischaemic from non-ischaemic heart failure. *Eur. J. Heart Fail.* **22**, 821–833 (2020).

19. Zhang, Q. H., Yin, R. X., Chen, W. X., Cao, X. L. & Chen, Y. M. Association between the TIMD4-HAVCR1 variants and serum lipid levels, coronary heart disease and ischemic stroke risk and atorvastatin lipid-lowering efficacy. *Biosci. Rep.* **38**, BSR20171058 (2018).

20. Gruson, D., Carbone, V., Feracin, B., Ahn, S. A. & Rousseau, M. F. Incremental value of intact fibroblast growth factor 23 to natriuretic peptides for long-term risk estimation of heart failure patients. *Clin. Biochem.* **61**, 47–49 (2018).

21. Lourenco, P. et al. Dynamics of growth differentiation factor 15 in acute heart failure. *ESC Heart Fail.* **8**, 2527–2534 (2021).

22. Gargiulo, P. et al. Surfactant-derived proteins as markers of alveolar membrane damage in heart failure. *PLoS ONE* **9**, e115030 (2014).

23. Brankovic, M. et al. Utility of temporal profiles of new cardio-renal and pulmonary candidate biomarkers in chronic heart failure. *Int. J. Cardiol.* **276**, 157–165 (2019).

24. Dubin, R. F. et al. Proteomic analysis of heart failure hospitalization among patients with chronic kidney disease: The Heart and Soul Study. *PLoS ONE* **13**, e0208042 (2018).

25. Klimczak-Tomaniak, D. et al. Temporal patterns of macrophage- and neutrophil-related markers are associated with clinical outcome in heart failure patients. *ESC Heart Fail.* **7**, 1190–1200 (2020).

26. Arpino, V., Brock, M. & Gill, S. E. The role of TIMPs in regulation of extracellular matrix proteolysis. *Matrix Biol.* **44–46**, 247–254 (2015).

27. Wang, H. et al. Tissue inhibitor of metalloproteinase 1 (TIMP-1) deficiency exacerbates carbon tetrachloride-induced liver injury and fibrosis in mice: Involvement of hepatocyte STAT3 in TIMP-1 production. *Cell Biosci.* **1**, 14 (2011).

28. Polina, E. R. et al. Relationship of polymorphisms in the tissue inhibitor of metalloproteinase (TIMP)-1 and -2 genes with chronic heart failure. *Sci. Rep.* **8**, 9446 (2018).

29. Jordan, A. et al. Matrix metalloproteinase-1 and its inhibitor, TIMP-1, in systolic heart failure: Relation to functional data and prognosis. *J. Intern. Med.* **262**, 385–392 (2007).
30. Krebber, M. M. et al. Matrix metalloproteinases and tissue inhibitors of metalloproteinases in extracellular matrix remodeling during left ventricular diastolic dysfunction and heart failure with preserved ejection fraction: A systematic review and meta-analysis. *Int. J. Mol. Sci.* **21**, 6742 (2020).

31. Tromp, J. et al. Identifying pathophysiological mechanisms in heart failure with reduced versus preserved ejection fraction. *J. Am. Coll. Cardiol.* **72**, 1081–1090 (2018).

## Author contributions

N.G. wrote the main manuscript text. Z.L. performed the statistical analysis. F.Z. coordinated the paper. All authors critically reviewed the manuscript.

## Funding

The COMMANDER HF trial was supported by Janssen Research and Development.

## Competing interests

N.G. reports honoraria from Novartis, AstraZeneca, Servier, Vifor, Bayer, Boehringer Ingelheim and Lilly. F.Z. and N.G. are supported by public grants overseen by the European Commission (EUFP7) and the French National Research Agency (ANR) as part of the second 'Investissements d'Avenir' programme FIGHT-HF (ANR-15-RHU-0004); programme HOMAGE under grant agreement N 305507; and project FIBROTARGETS grant HEALTH-2013-602904. S.D.A. reports receiving fees from Abbott, Bayer, Boehringer Ingelheim, Cardiac Dimension, Cordio, Impulse Dynamics, Novartis, Occlutech, Servier, and Vifor Pharma, and grant support from Abbott and Vifor Pharma. M.R.M. reports personal fees from Janssen Research & Development LLC as a member of the steering committee of the COMMANDER HF trial and consulting fees during the conduct of the study and personal fees from Abbott, Medtronic, Portola Pharmaceuticals, Bayer, Mesoblast, Baim Institute for Clinical Research, Xogenex, NupulseCV and FineHeart outside the submitted work. J.G.F.C. reports personal fees from Janssen Research & Development LLC as a member of the steering committee of the COMMANDER HF trial and research grants and honoraria for speaking, committees, and advisory boards from Amgen, AstraZeneca, Bayer, Bristol-Myers Squibb, GlaxoSmithKline, Johnson & Johnson, Medtronic, MyoKardia, Novartis, Philips Healthcare, Pharmacosmos, Pharma Nord, Sanofi, Laboratoires Servier, Stealth BioTherapeutics, Torrent Pharmaceuticals, and Vifor Pharma. C.S.P.L. reports personal fees from Janssen Research & Development LLC during the conduct of the study; research support from Boston Scientific, Bayer, Roche Diagnostics, AstraZeneca, Medtronic, and Vifor Pharma; personal fees from serving on the advisory board/steering committee/executive committee/clinical end points committee for Boston Scientific, Bayer, Roche Diagnostics, AstraZeneca, Medtronic, Vifor Pharma, Novartis, Amgen, Merck, Janssen Research & Development LLC, Menarini, Boehringer Ingelheim, Novo Nordisk, Abbott Diagnostics, Corvia Medical, Stealth BioTherapeutics, Jana Care, Biofourmis, Applied Therapeutics Inc, WebMD Global LLC, Radcliffe Group Ltd, Takeda Pharmaceutical Company, and Darma Laboratories outside the submitted work; and support from a Clinician Scientist Award from the National Medical Research Council of Singapore. D.J.v.V. reports personal fees from Johnson & Johnson and Bayer during the conduct of the study and reports receiving board membership fees from Janssen Research & Development LLC as a member of the steering committee of the COMMANDER HF trial. W.M. B. reports other support from Janssen Research & Development LLC during the conduct of the study and is employed by and is a shareholder of Johnson & Johnson. B.G. reports personal fees from Janssen Research & Development LLC as a member of the steering committee of the COMMANDER HF trial and consulting fees from Bayer during the conduct of the study; consulting fees from Janssen Research & Development LLC, Bayer, Novartis, Mesoblast, Ionis Pharmaceuticals, Zensun USA, and Cellular Dynamics outside the submitted work; speaker fees from Novartis and Otsuka Pharmaceutical outside the submitted work; and personal fees from AstraZeneca, Actelion, Amgen, EBR Systems, Impulse Dynamics, MyoKardia, Rocket Pharma, Sanofi, and Viking Therapeutics outside the submitted work. F.Z. reports personal fees from Janssen Research & Development LLC and Bayer during the conduct of the study; personal fees from AstraZeneca, Boehringer Ingelheim, LivaNova, GE Healthcare, Amgen, Novartis, Quantum Genomics, Cardior Pharmaceuticals, CardioRenal, CVCT, Merck, CVRx, Vifor Fresenius Medical Care Renal Pharma, NovoNordisk, and MundiPharma outside the submitted work; fees for serving on a steering committee or a safety and data monitoring board from Actelion, Amgen, Bayer, Boehringer Ingelheim, Boston Scientific, CVRx, GE Healthcare, Janssen Research & Development LLC, Novartis, and ResMed; and consulting fees from AstraZeneca, Cardior Pharmaceuticals, CardioRenal, Quantum Genomics, and Vifor Fresenius Medical Care Renal Pharma. No other authors have any competing-interests.

## Additional information

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1038/s41598-022-12385-0.

Correspondence and requests for materials should be addressed to F.Z.

Reprints and permissions information is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
<figure><img src="image_6.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International

License, which permits use, sharing, adaptation, distribution and reproduction in any medium or
format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the
Creative Commons licence, and indicate if changes were made. The images or other third party material in this
article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the
material. If material is not included in the article's Creative Commons licence and your intended use is not
permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from
the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2022