# Proteomic mechanistic profile of patients with diabetes at risk of developing heart failure: insights from the HOMAGE trial

Job a J Verdonschot, João Pedro Ferreira, Pierpaolo Pellicori, Hans-Peter Brunner-La Rocca, Andrew L Clark, Franco Cosmi, Joe Cuthbert, Nicolas Girerd, Beatrice Mariottoni, Johannes Petutschnigg, et al.

► To cite this version:

Job a J Verdonschot, João Pedro Ferreira, Pierpaolo Pellicori, Hans-Peter Brunner-La Rocca, Andrew L Clark, et al. Proteomic mechanistic profile of patients with diabetes at risk of developing heart failure: insights from the HOMAGE trial. Cardiovascular Diabetology, 2021, 20 (1), pp.163. 10.1186/s12933-021-01357-9. hal-03320862

HAL Id: hal-03320862

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

Submitted on 16 Aug 2021

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
ORIGINAL INVESTIGATION

Open Access

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

Check for
updates

# Proteomic mechanistic profile of patients with diabetes at risk of developing heart failure: insights from the HOMAGE trial

Job A. J. Verdonschot<sup>1†</sup>, João Pedro Ferreira<sup>2†</sup>, Pierpaolo Pellicori<sup>3</sup>, Hans-Peter Brunner-La Rocca<sup>4</sup>, Andrew L. Clark<sup>5</sup>, Franco Cosmi<sup>6</sup>, Joe Cuthbert<sup>5</sup>, Nicolas Girerd<sup>2</sup>, Beatrice Mariottoni<sup>6</sup>, Johannes Petutschnigg<sup>7</sup>, Patrick Rossignol<sup>2</sup>, John G. F. Cleland<sup>3</sup>, Faiez Zannad<sup>2†</sup> and Stephane R. B. Heymans<sup>4*†</sup> on behalf of HOMAGE “Heart Omics in AGEing” consortium

**Abstract**

**Background:** Patients with diabetes mellitus (DM) are at increased risk of developing heart failure (HF). The “Heart OMics in AGEing” (HOMAGE) trial suggested that spironolactone had beneficial effect on fibrosis and cardiac remodeling in an at risk population, potentially slowing the progression towards HF. We compared the proteomic profile of patients with and without diabetes among patients at risk for HF in the HOMAGE trial.

**Methods:** Protein biomarkers (n = 276) from the Olink®Proseek-Multiplex cardiovascular and inflammation panels were measured in plasma collected at baseline and 9 months (or last visit) from HOMAGE trial participants including 217 patients with, and 310 without, diabetes.

**Results:** Twenty-one biomarkers were increased and five decreased in patients with diabetes compared to non-diabetics at baseline. The markers clustered mainly within inflammatory and proteolytic pathways, with granulin as the key-hub, as revealed by knowledge-induced network and subsequent gene enrichment analysis. Treatment with spironolactone in diabetic patients did not lead to large changes in biomarkers. The effects of spironolactone on NTproBNP, fibrosis biomarkers and echocardiographic measures of diastolic function were similar in patients with and without diabetes (all interaction analyses p > 0.05).

**Conclusions:** Amongst patients at risk for HF, those with diabetes have higher plasma concentrations of proteins involved in inflammation and proteolysis. Diabetes does not influence the effects of spironolactone on the proteomic profile, and spironolactone produced anti-fibrotic, anti-remodelling, blood pressure and natriuretic peptide lowering effects regardless of diabetes status.

Trial registration NCT02556450.

**Keywords:** Diabetes, Biomarker, Heart failure, Spironolactone, Fibrosis

**Background**

Patients with diabetes mellitus (DM) have an increased risk of developing heart failure (HF) [1] but the mechanisms are poorly established. A better understanding of the pathophysiological processes linking DM and HF may help guiding treatment strategies and potentially reveal novel therapeutic targets [2–4].

<footnote>*Correspondence: s.heymans@maastrichtuniversity.nl</footnote>

<footnote>†Job A. J. Verdonschot and João Pedro Ferreira are joint first authors</footnote>

<footnote>Faiez Zannad and Stephane R. B. Heymans are joint last authors</footnote>

<footnote n="4">⁴Department of Cardiology, Maastricht University Medical Center (MUMC+), PO Box 5800, 6202 AZ Maastricht, The Netherlands
Full list of author information is available at the end of the article</footnote>
The HOMAGE (Heart OMics in AGEing) trial showed that treatment with spironolactone (vs. usual care) in patients at risk of developing HF led to a decrease in collagen synthesis markers, N-terminal pro brain natriuretic peptide (NT-proBNP), reduced blood pressure and improved cardiac remodelling [5, 6]. Small studies have suggested that spironolactone could impair endothelial function and increase glycated hemoglobin (HbA1c) in patients with DM [7, 8]. As spironolactone is a steroidal mineralocorticoid receptor antagonist (MRA) that is not specific for the mineralocorticoid receptor, it may exert other effects (e.g., cortisol elevation) on glucose and insulin metabolism. In a previous analysis we showed the influence of spironolactone on the proteome in the total HOMAGE trial population [9]. It remains unknown to what extent DM influences the effects of spironolactone on clinical and proteomic variables.

In this pre-specified secondary analysis of the HOMAGE trial, we aimed to (1) compare the clinical characteristics of patients with and without DM; (2) study the biomarker profiles, networks and pathways associated with DM at baseline, as compared to patients without diabetes; (3) to assess the effect of spironolactone on the circulating proteomic biomarkers in patients with DM; and (4) to explore whether the DM status might have modified the effect of spironolactone on the main outcomes assessed in the HOMAGE trial (level of collagen synthesis markers, NT-proBNP, echocardiographic measures of cardiac structure and function, and signs/symptoms).

# Methods

## Trial design and population

The HOMAGE trial had a prospective, randomised, open-label, blinded-endpoint (PROBE), multicentre design. People at increased risk of developing HF were randomly assigned to receive either spironolactone or standard care (ClinicalTrials.gov Identifier: NCT02556450). The rationale, trial design and main results have been published [5, 6]. The study was approved by all relevant ethics committees and regulatory bodies. All participants provided written informed consent prior to study-specific procedures. The main participating criteria included age 65 or older (amended to 60 years during the course of the trial), cardiovascular risk defined by the presence of coronary artery disease OR at least 2 of the following: diabetes, treated hypertension, microalbuminuria, abnormal ECG, and a NT-proBNP between 125 and 1000 ng/l or BNP between 35 and 280 ng/l. The main exclusion criteria were glomerular filtration rate (eGFR) < 30 ml/min/1.73 m², serum potassium > 5.0 mmol/l, left ventricular ejection fraction < 45%, a diagnosis of HF or treatment with loop diuretics, and

atrial fibrillation flutter. For this secondary analysis, patients were divided into two groups based on the presence, or absence, of a diagnosis of DM at baseline.

## Proteomic biomarkers

Baseline and month 9 (or last visit) plasma samples were analysed for 276 protein biomarkers by the TATAA-bio-center using the Olink Proseek® Multiplex cardiovascular (CVD) II, CVD III, and inflammation panels. The proteins were determined using high-throughput Olink Proseek® Multiplex 96 × 96 kits, which measures 92 manually-selected proteins simultaneously in 1 µl of plasma per kit. Each kit uses a proximity extension assay (PEA) technology with dual-recognition DNA-coupled read-out, where 92 oligonucleotide-labelled antibody probe pairs are allowed to bind to their respective targets in the sample. The platform provides Log2 normalized protein expression (NPX) values with relative quantification. A detailed description of the Olink® technology is depicted on the website: https://www.olink.com/. The abbreviations, full names and respective Olink® multiplex panels of the studied proteins are described in Additional file 1: Table S1. The assays were performed in a blinded fashion to the treatment allocation. The proteomic results were then merged into the database.

## Statistical analyses

We compared the characteristics of the patients with and without diabetes using the appropriate tests for continuous and categorical variables. To assess whether the biomarkers were expressed differently in patients with and without diabetes, we performed logistic regression analyses with diabetes as outcome variable and adjusting for age, sex, coronary artery disease, hypertension, body mass index and estimated glomerular filtration rate. Since proteins were measured using NPX values on a log2 scale, the odds ratio for each protein estimates the increase in the odds of diabetic status associated with a doubling in the protein concentration. We corrected the findings for multiple testing using a false discovery rate of 5%. After selecting the proteins with differential expression by diabetes status, we tested whether spironolactone could change the levels of the proteins. For consistency with the primary report, we used analysis of covariance (ANCOVA) comparing the difference in any changes between the control and spironolactone groups in the regression model [6].

A linear regression model was fitted, with the protein change (from baseline to last visit) as outcome variable, a binary variable to indicate the treatment group (control/spironolactone), and the baseline protein value (NPX) as covariates. The treatment effect was the coefficient that resulted from the comparison of spironolactone vs.
control in the regression model. Residual analysis was used to examine the fit of the model. No data transformation was required to meet the assumptions of linear regression. Similar analyses were performed for the protein change at 1 month. To study if the diabetes status could influence the response to spironolactone on the main outcomes of the study, we performed an ANCOVA analysis with a treatment-by-diabetes interaction term. Statistical analyses were performed using Stata® (version 16, StataCorp LP).

## Bioinformatical and network analyses

We used knowledge-based network analysis with induced network approach by consensuspathDB (CPDB) online server (accessed on 24 November 2020) from Max Planck Institute for Molecular Genetics to identify the links among all significantly different protein biomarkers, based on knowledge of interactions (protein interactions, genetic interactions and biochemical interactions) [7]. The network analysis also identifies additional proteins limited to the first-degree interactors (intermediate nodes) linking our input proteins (seed nodes), with exclusion of low-confidence interactions and quantified by a z-score ≤ 20 calculated for each intermediate node. The Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) database was used to analyse functional enrichment (GO biological processes and KEGG pathways) using proteins that were significantly increased or decreased in diabetic patients compared to non-diabetics at baseline.

## Results

### Clinical characteristics of the study population

Among the 527 patients included in the HOMAGE trial, 217 (41%) had diabetes. Compared to those without diabetes, patients with diabetes had a higher BMI, a higher blood pressure and heart rate (Table 1), were more likely to have a history of hypertension (93% vs. 68%; p < 0.001), but less likely to have a history of coronary artery disease (47% vs. 90%; p < 0.001). Among patients with diabetes, 176 (81%) were on metformin, and 38 (18%) were insulin-dependent. Left ventricular mass was greater (99 [84–114] vs. 91 [79–110] g/m²; p = 0.005) and E/e' ratio higher (10.1 [8–12] vs. 8.9 [7–11]; p < 0.001) in patients with diabetes; however, left atrial volume and NT-proBNP were similar.

### Proteomic profile at baseline

Compared to patients without diabetes, those with diabetes had greater expression of 21 proteins and lower expression of 5 proteins after adjustment for age, sex, hypertension, body mass index, coronary artery disease,

estimated glomerular filtration rate and multiple test correction at a FDRq < 0.05 level (Table 2).

The results of the network analysis are shown in Fig. 1. The central core of the network is characterized by (lysosomally) protease activity reflected in matrix metalloproteinase 7 (MMP7), cathepsin D (CTSD), cathepsin Z (CTSZ) and prostasin (PRSS8). A large number of edges and nodes connect through granulin (GRN), placing this protein as an important hub connecting the proteases with the chemokines [cluster of chemokine (CC-motif) ligand (CCL) 11 and 15 and chemokine (CXC motif) ligand (CXCL) 10]. Pathway enrichment analysis using the Kyoto Encyclopedia of Genes and Genomes (KEGG) terms placed 13 proteins in a biological pathway involving inflammation or lysosomal proteolytic activity (Table 3). An additional analysis using the Gene Ontology terms, showed that the top 10 enriched pathways were all implicated in inflammation (Fig. 2; Additional file 1: Table S2).

### Proteomic changes after spironolactone in diabetes patients

Three proteins fell significantly (COL1A1, PICP and SELE) and one protein increased (MMP7) after 9 months of spironolactone treatment (Additional file 1: Table S3). Overall, there was no prominent effect of spironolactone on the proteomic markers in diabetic patients. The activity of the biological processes at baseline in diabetic patients (Table 3; Fig. 2) was unaltered after 9 months of spironolactone treatment.

### Interaction of diabetes status with the main trial outcomes

The primary and secondary outcomes measures in the HOMAGE trial were: PICP and PIIINP (markers of collagen synthesis), NT-proBNP, systolic and diastolic blood pressure, echocardiographic left atrial volume and left ventricular mass. The presence of diabetes did not modify the effect of spironolactone on any of those outcomes (Additional file 1: Table S4), indicating that there was no superior or inferior effect of spironolactone on serum markers of collagen metabolism, cardiac structure and function in patients with DM.

## Discussion

This study investigated the proteomic profile of patients at risk for heart failure treated with spironolactone, with a special focus on DM (Fig. 3). The main findings of the study are: (i) patients with diabetes have a distinct proteome profile, whereby twenty-one biomarkers are higher, while five were lower compared to those without diabetes; (ii) cluster analysis revealed that the distinct proteomic profile is indicative of increased activation of inflammatory and proteolytic pathways in patients with
<table><caption>**Table 1** Baseline characteristics according to diabetic status</caption><thead><tr><th>Variable</th><th>Level</th><th>Patients without diabetes (n = 310)</th><th>Patients with diabetes (n = 217)</th><th>p-value</th></tr></thead><tbody><tr><td>Age (years)</td><td></td><td>71.9 [68.5-78.1]</td><td>73.5 [69.0-79.0]</td><td>0.064</td></tr><tr><td>Male gender</td><td></td><td>235 (75.8%)</td><td>157 (72.4%)</td><td>0.37</td></tr><tr><td rowspan="4">Smoking status</td><td>Never</td><td>102 (32.9%)</td><td>76 (35.0%)</td><td rowspan="4">0.90</td></tr><tr><td>Ex</td><td>177 (57.1%)</td><td>123 (56.7%)</td></tr><tr><td>Current</td><td>28 (9.0%)</td><td>16 (7.4%)</td></tr><tr><td>Missing</td><td>3 (1.0%)</td><td>2 (0.9%)</td></tr><tr><td>Body mass index (kg/m²)</td><td></td><td>27.3 [25.0-29.9]</td><td>29.7 [26.5-33.9]</td><td>&lt;0.001</td></tr><tr><td>Systolic blood pressure (mmHg)</td><td></td><td>139.0 [124.0-154.0]</td><td>141.0 [132.0-157.0]</td><td>0.013</td></tr><tr><td>Diastolic blood pressure (mmHg)</td><td></td><td>79.0 [72.0-86.0]</td><td>77.0 [71.0-84.0]</td><td>0.097</td></tr><tr><td>Heart rate (pm)</td><td></td><td>59.0 [54.0-65.0]</td><td>64.0 [57.0-70.0]</td><td>&lt;0.001</td></tr><tr><td>Number of shuttles completed (shuttle-walk test)</td><td></td><td>54 [35-71]</td><td>44 [26-61]</td><td>&lt;0.001</td></tr><tr><td>Comorbidities</td><td></td><td></td><td></td><td></td></tr><tr><td>Coronary artery disease</td><td></td><td>278 (89.7%)</td><td>101 (46.5%)</td><td>&lt;0.001</td></tr><tr><td>Myocardial infarction</td><td></td><td>163 (58.6%)</td><td>51 (50.0%)</td><td>0.13</td></tr><tr><td>PCI</td><td></td><td>196 (70.5%)</td><td>69 (67.6%)</td><td>0.59</td></tr><tr><td>CABG</td><td></td><td>97 (34.9%)</td><td>39 (38.2%)</td><td>0.55</td></tr><tr><td>Hypertension</td><td></td><td>211 (68.1%)</td><td>202 (93.1%)</td><td>&lt;0.001</td></tr><tr><td>Stroke/TIA</td><td></td><td>14 (4.5%)</td><td>14 (6.5%)</td><td>0.33</td></tr><tr><td>COPD</td><td></td><td>22 (7.1%)</td><td>11 (5.1%)</td><td>0.34</td></tr><tr><td>Medication</td><td></td><td></td><td></td><td></td></tr><tr><td>Antiplatelet therapy</td><td></td><td>267 (86.1%)</td><td>147 (67.7%)</td><td>&lt;0.001</td></tr><tr><td>Beta blocker</td><td></td><td>239 (77.1%)</td><td>127 (58.5%)</td><td>&lt;0.001</td></tr><tr><td>ACE-inhibitor</td><td></td><td>157 (50.6%)</td><td>118 (54.4%)</td><td>0.40</td></tr><tr><td>Angiotensin receptor blocker</td><td></td><td>77 (24.8%)</td><td>68 (31.3%)</td><td>0.10</td></tr><tr><td>Calcium channel blocker</td><td></td><td>56 (18.1%)</td><td>54 (24.9%)</td><td>0.058</td></tr><tr><td>Thiazides</td><td></td><td>33 (10.6%)</td><td>54 (24.9%)</td><td>&lt;0.001</td></tr><tr><td>Statin/lipid lowering drug</td><td></td><td>269 (86.8%)</td><td>166 (76.5%)</td><td>0.002</td></tr><tr><td>Insulin</td><td></td><td>0 (0.0%)</td><td>38 (17.5%)</td><td>&lt;0.001</td></tr><tr><td>Metformin</td><td></td><td>0 (0.0%)</td><td>176 (81.1%)</td><td>&lt;0.001</td></tr><tr><td>Sulfunylureas</td><td></td><td>0 (0.0%)</td><td>41 (18.9%)</td><td>&lt;0.001</td></tr><tr><td>DPP4-inhibitor</td><td></td><td>0 (0.0%)</td><td>33 (15.2%)</td><td>&lt;0.001</td></tr><tr><td>Echocardiography</td><td></td><td></td><td></td><td></td></tr><tr><td>Left ventricular ejection fraction (%)</td><td></td><td>62.8 [56.5-66.0]</td><td>63.1 [59.2-66.9]</td><td>0.35</td></tr><tr><td>Left ventricular mass (g/m²)</td><td></td><td>91.0 [78.5-109.9]</td><td>99.0 [84.2-114.1]</td><td>0.005</td></tr><tr><td>Left ventricular hypertrophy</td><td></td><td>106 (34.2%)</td><td>87 (40.1%)</td><td>0.17</td></tr><tr><td>Left atrial volume (g/m²)</td><td></td><td>30.8 [26.0-36.5]</td><td>29.4 [25.1-35.6]</td><td>0.32</td></tr><tr><td>E/e' ratio</td><td></td><td>8.9 [6.9-11.1]</td><td>10.1 [8.3-12.0]</td><td>&lt;0.001</td></tr><tr><td>E/A ratio</td><td></td><td>0.8 [0.7-1.0]</td><td>0.8 [0.7-1.0]</td><td>0.21</td></tr><tr><td>TAPSE (mm)</td><td></td><td>21.4 [16.7-26.3]</td><td>22.5 [17.5-26.8]</td><td>0.28</td></tr><tr><td>MAPSE (mm)</td><td></td><td>15.5 [13.6-17.7]</td><td>15.0 [13.1-17.4]</td><td>0.38</td></tr><tr><td>Laboratory</td><td></td><td></td><td></td><td></td></tr><tr><td>eGFR (ml/min/1.73 m²)</td><td></td><td>72.0 [62.0-82.9]</td><td>72.8 [59.0-86.5]</td><td>0.75</td></tr><tr><td>Urea (mmol/l)</td><td></td><td>7.5 [5.6-13.2]</td><td>10.7 [6.4-15.7]</td><td>&lt;0.001</td></tr><tr><td>Haemoglobin (g/dl)</td><td></td><td>14.2 [13.3-15.0]</td><td>13.8 [12.8-14.7]</td><td>0.002</td></tr><tr><td>Sodium (mmol/l)</td><td></td><td>140.0 [138.0-141.0]</td><td>139.0 [137.0-141.0]</td><td>0.052</td></tr><tr><td>Potassium (mmol/l)</td><td></td><td>4.3 [4.1-4.6]</td><td>4.3 [4.0-4.5]</td><td>0.12</td></tr><tr><td>NT-proBNP, ng/l</td><td></td><td>205 [129-350]</td><td>229 [146-373]</td><td>0.25</td></tr><tr><td>Galectin-3, ng/ml</td><td></td><td>15.4 [13.1-18.9]</td><td>17.4 [14.4-20.9]</td><td>&lt;0.001</td></tr><tr><td>High-sensitivity Troponin T, pg/ml</td><td></td><td>11.6 [8.2-16.0]</td><td>14.1 [10.3-21.1]</td><td>&lt;0.001</td></tr></tbody></table>
**Table 1** (continued)

Values are median [interquartile range] or n (%)

PCI percutaneous coronary intervention, CABG coronary artery bypass grafting, TIA transient ischemic attack, COPD chronic obstructive pulmonary disease, ACE angiotensin-converting enzyme, DPP4 dipeptidyl peptidase-4, TAPSE triscupid annular plane systolic excursion, MAPSE mitral annular plane systolic excursion, eGFR estimated glomerular filtration rate

<table><caption>**Table 2** Odds ratio of the proteins which are significantly associated with diabetes after adjustment for age, sex, hypertension, BMI, CAD and eGFR</caption><thead><tr><th>Protein</th><th>Abbreviation</th><th>Odds ratio (95% CI)</th><th>p-value</th><th>FDRq</th></tr></thead><tbody><tr><td colspan="5">Positive association (n = 21)</td></tr><tr><td>Interleukin-1 receptor type 1</td><td>il1rt1</td><td>4.52 (2.41–8.5)</td><td>0.00001</td><td>0.0005</td></tr><tr><td>Cathepsin d</td><td>ctsd</td><td>3.49 (2.13–5.72)</td><td>0.00001</td><td>0.0005</td></tr><tr><td>Growth differentiation factor 15</td><td>gdf15</td><td>3.41 (2.39–4.87)</td><td>0.00001</td><td>0.0005</td></tr><tr><td>Galectin-4</td><td>gal4</td><td>2.95 (1.95–4.45)</td><td>0.00001</td><td>0.0005</td></tr><tr><td>E-selectin</td><td>sele</td><td>2.38 (1.64–3.45)</td><td>0.00001</td><td>0.0005</td></tr><tr><td>Eotaxin</td><td>ccl11</td><td>2.96 (1.79–4.9)</td><td>0.00002</td><td>0.0009</td></tr><tr><td>Leukemia inhibitory factor receptor</td><td>lifr</td><td>4.4 (2.15–9.04)</td><td>0.00005</td><td>0.0017</td></tr><tr><td>Kidney injury molecule</td><td>kim1</td><td>1.79 (1.35–2.38)</td><td>0.00005</td><td>0.0017</td></tr><tr><td>Cathepsin z</td><td>ctsz</td><td>3.1 (1.78–5.41)</td><td>0.00007</td><td>0.0019</td></tr><tr><td>Hydroxyacid oxidase 1</td><td>haox1</td><td>1.44 (1.2–1.72)</td><td>0.00007</td><td>0.0019</td></tr><tr><td>P-selectin glycoprotein ligand 1</td><td>psgl1</td><td>4.23 (1.95–9.16)</td><td>0.00026</td><td>0.0054</td></tr><tr><td>Interleukin-18 receptor 1</td><td>il18r1</td><td>2.35 (1.43–3.87)</td><td>0.00074</td><td>0.014</td></tr><tr><td>Prostasin</td><td>prss8</td><td>2.97 (1.55–5.69)</td><td>0.00102</td><td>0.0173</td></tr><tr><td>C–C motif chemokine 15</td><td>ccl15</td><td>2.13 (1.34–3.38)</td><td>0.00135</td><td>0.0215</td></tr><tr><td>Interleukin-10</td><td>il10</td><td>1.74 (1.23–2.45)</td><td>0.00168</td><td>0.024</td></tr><tr><td>Gastric intrinsic factor</td><td>gif</td><td>1.34 (1.11–1.6)</td><td>0.00196</td><td>0.0253</td></tr><tr><td>Matrix metalloprotease 7</td><td>mmp7</td><td>1.92 (1.27–2.9)</td><td>0.00187</td><td>0.0253</td></tr><tr><td>Granulin</td><td>grn</td><td>2.66 (1.42–5)</td><td>0.00237</td><td>0.0292</td></tr><tr><td>Follistatin</td><td>fs</td><td>2.16 (1.29–3.61)</td><td>0.00333</td><td>0.0392</td></tr><tr><td>Alpha-L-iduronidase</td><td>idua</td><td>2.04 (1.26–3.29)</td><td>0.0036</td><td>0.0407</td></tr><tr><td>Osteoprotegerin</td><td>opg</td><td>2.43 (1.32–4.45)</td><td>0.00422</td><td>0.044</td></tr><tr><td colspan="5">Negative association (n = 5)</td></tr><tr><td>Collagen type 1 alpha 1</td><td>col1a1</td><td>0.33 (0.19–0.58)</td><td>0.00013</td><td>0.003</td></tr><tr><td>Procollagen type 1 carboxy-terminal propeptide</td><td>picp</td><td>0.36 (0.22–0.62)</td><td>0.00016</td><td>0.0036</td></tr><tr><td>Lipoprotein lipase</td><td>lpl</td><td>0.48 (0.31–0.74)</td><td>0.00085</td><td>0.0154</td></tr><tr><td>Matrix extracellular phosphoglycoprotein</td><td>mepe</td><td>0.43 (0.26–0.73)</td><td>0.00167</td><td>0.024</td></tr><tr><td>C-X-C motif chemokine 10</td><td>cxcl10</td><td>0.65 (0.48–0.87)</td><td>0.00418</td><td>0.044</td></tr></tbody></table>

diabetes; and (iii) diabetic status did not influence the treatment response to spironolactone with respect to change in the proteomic profile, as well as to the anti-fibrotic, anti-remodelling, blood pressure and natriuretic peptide lowering effects. It is unlikely that diabetic status modified the effect of spironolactone observed in the HOMAGE trial.

**Proteomic profile of patients with diabetes**

Patients with DM had an up-regulation of inflammatory and proteolytic pathways, suggesting that inflammation and proteolysis mechanisms (i.e. an increased

protein turn-over via the lysosomal proteolytic system) are active and distinctive for patients with DM compared to non-diabetic patients at risk of developing HF. Based on our data we can not state any causal relationship between these proteins and pathways and diabetes. However, it seems that the observations are independent of coronary artery disease, hypertension, BMI, age, sex and kidney function. Our findings are in line with a previous network and pathway analysis in patients with both diabetes and HF, showing a similar association between DM and inflammation, protein phosphorylation, and neutrophil degranulation mechanistic
<figure><img src="image_3.png" /><figcaption>**Table 3** Pathway enrichment analysis using KEGG terms</figcaption></figure>

<table><thead><tr><td>KEGG-term</td><td>Pathway</td><td>FDRq</td><td>Included proteins</td></tr></thead><tbody><tr><td>hsa04060</td><td>Cytokine-cytokine receptor interaction</td><td>0.0000014</td><td>CCL23, CCL15, CCL11, CXCL10, LIFR, IL10, IL1R1, OPG, IL18R1</td></tr><tr><td>hsa04668</td><td>TNF signaling pathway</td><td>0.031</td><td>CXCL10, SELE, IL18R1</td></tr><tr><td>hsa04142</td><td>Lysosome</td><td>0.039</td><td>CTSD, CTSZ, IDUA</td></tr><tr><td>hsa04062</td><td>Chemokine signaling pathway</td><td>0.084</td><td>CCL11, CCL15, CXCL10</td></tr></tbody></table>

pathways [4]. The top 5 proteins we identified (GDF15, GAL4, IL1RT1, CTSD and SELE) are robustly associated with diabetes in a variety of studies, consistent with our findings [2–4].

Growth differentiation factor 15 (GDF15) is a protein which is upregulated as consequence of organ injury, as seen in heart and kidney disease [10]. In patients with HF, GDF15 is a strong prognostic marker [11]. However, in diabetic patients without HF, GDF15 may be associated with the effect of metformin. Metformin partly works through stimulation of GDF15 expression, which is associated with weight loss and improved insulin sensitivity in patients with DM [12, 13]. On the other hand, GDF15 levels are increased in patients with HF, and GDF15 is associated with diastolic dysfunction,

as defined by an elevated E/e' [14, 15], which was also observed in the patients with DM in our study.

Galectin-4 (GAL4) is a small lectin protein and is involved in protein trafficking from the inside of the cell towards the cell membrane. This protein was detected in all previous proteomics screenings in diabetic patients, confirming a strong association between DM and GAL4. One of the proteins which is trafficked by GAL4 is the protease dipeptidyl peptidase-4 (DPP4). DPP4 subsequently deactivates proglucagon-derived peptide glucagon-like peptide-1 (GLP1). Inhibitors of DPP4 and GLP1-analogues are both used for lowering glucose in patients with DM but do not improve HF outcomes. An increase in GAL4 leads to increased activity of DPP4, which reduces the activity of GLP1, contributing to
Enriched GO processes

<figure><img src="image_4.png" /><figcaption>**Fig. 2** Overrepresented pathways of protein biomarkers in diabetic patients when compared to non-diabetic patients. The gene ontology (GO) biological processes are used as reference. Additional file 1: Table S2 contains further details on the proteins in these pathways</figcaption></figure>

insulin resistance and DM development [3]. Although the association between GAL4 and DM is established, the exact mechanistic and causal role remains unknown.

Interleukin-1 receptor type 1 (IL1RT1) is the receptor for IL1A, IL1B and IL1RA and is responsible for the IL1-dependent activation of the NFκB pathway. However, recently IL1R1 has been identified as one of the key mediators in leptin sensitivity and IL1R1 deficiency leads to a higher degree of obesity and metabolic disturbance [16]. Pro-inflammatory signalling can have beneficial metabolic effects: for example, inflammatory signalling is essential for proper adipose tissue remodelling and expansion which prevents ectopic lipid accumulation and subsequent glucose intolerance [17, 18]. Cytokine resistance as a consequence of chronic systemic inflammation could be one of the mechanisms underlying the development of DM type 2. It may be that the higher levels of inflammatory proteins we observed could be partially explained by cytokine resistance in diabetic patients.

Cathepsin D (CTSD) is a lysosomal endopeptidase which is involved in intracellular protein turnover and extracellular matrix breakdown. Levels of CTSD have been associated with DM in patients with and without HF [2, 3, 19], and are positively associated with insulin resistance [20]. The family of cathepsins were well-represented in the network analysis, showing a link between inflammation and proteolysis. It has been suggested that

CTSD acts as a mediator between obesity and systemic inflammation, which could contribute to the cytokine resistance mentioned above [21].

E-selectin (SELE) is a selectin cell adhesion molecule expressed only on endothelial cells activated by cytokines. Levels of SELE and other cellular adhesion molecules have previously been associated with DM, increasing HbA1c levels, blood pressure and microalbuminuria, all reflecting (early) endothelial dysfunction [22–24].

Granulin as important node in the diabetes pathophysiology

Granulin was identified as a key hub in the network based on the proteomic profile of patients with DM. Granulin is cleaved from the precursor progranulin and their exact function remains unknown, but both forms have been implicated in inflammation and protein homeostasis via the lysosomal pathway (e.g. CTSD). Mechanistically (pro)granulin seems to play an important role in the proteomic network of diabetic patients and is related to insulin resistance, obesity, inflammation and proteolysis [25]. Mice deficient in Grn have lower body weight, fat mass and adipocyte size than wild-type mice receiving a standard diet. The loss of granulin protected the mice against obesity and insulin resistance induced by a high fat diet [25]. Granulin increases IL6-mediated expression
<figure><img src="image_5.png" /></figure>

of cytokine signalling 3 (SOCS3), which inhibits phosphorylation of insulin receptor substrate-1 (IRS1) [25]. A previous study has also inferred that granulin is a key hub in the proteomic profile of patients with HF and DM based on other biomarkers (although granulin itself was not actually measured) [4]. We now identify granulin as a biomarker which is elevated in diabetic patients.

## Spironolactone treatment in patients with diabetes

The study end of the HOMAGE trial was determined at 9 months, as it is expected that the biochemical response of spironolactone can be observed after 9 months of treatment. At the last visit of a patient, additional blood sampling was performed for proteomic analysis. Spironolactone did not significantly change the proteomic profiles of patients with diabetes, but yielded the same clinical benefits as in the complete HOMAGE cohort. The four proteins which were significantly modified in diabetic patients after 9 months of spironolactone treatment (COL1A1, PICP, SELE and MMP7) are altered in the same direction as in the complete population, indicating that the effects are not specific for diabetic patients [9]. These proteins are implicated in extracellular matrix

metabolism (COL1A1, PICP), cell-cell adhesion (SELE) and proteolysis (MMP7). Furthermore, diabetes did not modify the effect of spironolactone on serum markers of collagen metabolism, indicating that patients with diabetes benefited equally from the antifibrotic effect of spironolactone compared to non-diabetics. No metabolic pathways detrimental for diabetic patients was worsened with spironolactone.

## Limitations

We tested the effect of spironolactone on multiple proteins applying a correction for test multiplicity to limit the occurrence of false positive findings; however, as HOMAGE was a randomized controlled trial, other proteins, whose levels were also significantly changed with spironolactone, might also be implicated in relevant pathways and biological processes and might be worth exploring in further studies. Additionally, many of the highlighted mechanisms should be furtherly replicated and confirmed at a cellular level. Our results and previous findings highlight granulin as an important hub in the pathophysiology of diabetes which should receive priority for further study. Due to the
inclusion criteria of the HOMAGE study, the non-diabetic patient group mainly included patients with CAD, and thus we corrected all analyses for the presence of CAD in the patient groups. The HOMAGE study aimed to identify patients at risk for HF, therefore all of the included patients were above the age of 60 years old and had a cardiovascular risk profile. Therefore, these results can not be generalized to other populations. The Olink multiplex panels provide normalized protein expression (NPX) values, which are relative quantifications of the measured biomarkers. Absolute numbers of expression levels were not available.

## Conclusions

Diabetic patients at risk of developing HF have a proteomic bioprofile which differs from controls without diabetes, pointing towards mechanisms related to inflammation and protein homeostasis. Treatment with spironolactone did not influence the overall proteomic profile of diabetic compared to non-diabetic patients. Patients benefited similarly from spironolactone with respect to its anti-fibrotic, anti-remodeling, blood pressure and natriuretic peptide lowering effects regardless of the diabetes status.

## Supplementary Information

The online version contains supplementary material available at https://doi.org/10.1186/s12933-021-01357-9.

**Additional file 1: Table S1.** Protein names and respective Olink® panel sorted in alphabetical order. **Table S2.** Pathway enrichment analysis using GO terms. **Table S3.** Proteins which are significantly altered in diabetic patients after 9 months of spironolactone treatment. **Table S4.** Interaction analysis between diabetic status of the patients and spironolactone treatment for the primary and secondary outcomes of the HOMAGE trial.

### Acknowledgements

Not applicable.

### Authors' contributions

JV and JPF designed the study, analysed the data and drafted the manuscript. PP, JC, BM and JP collected the data and substantially revised the manuscript. HPBLR, AC, FC, NG, PR, JC, FZ and SH gathered the funding, supervised the study at their respective location, interpreted the data and substantially revised the manuscript. All authors read and approved the final manuscript.

### Funding

The research leading to these results has received funding from the European Union Commission's Seventh Framework programme under Grant Agreement No. 305507 (HOMAGE). S.H: This manuscript has been possible thanks to the support of the Netherlands Cardiovascular Research Initiative, an initiative with support of the Dutch Heart Foundation, CVON2016-Early HFPEF, 2015-10, CVON She-PREDICTS, Grant 2017-21.

### Availability of data and materials

The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.

## Declarations

### Ethics approval and consent to participate

The study was approved by all relevant ethics committees and regulatory bodies. All participants provided written informed consent prior to study-specific procedures.

### Consent for publication

There is no data of individual persons included in the manuscript.

### Competing interests

The authors have no relevant conflicts of interest to disclose with regards to the content of this manuscript.

### Author details

¹Department of Clinical Genetics, Maastricht University Medical Center, Maastricht, The Netherlands. ²Centre d'Investigations Cliniques Plurithématique 1433, Université de Lorraine, Institut National de la Santé et de la Recherche Médicale 1116, Centre Hospitalier Régional Universitaire de Nancy, French Clinical Research Infrastructure Network, Investigation Network Initiative-Cardiovascular and Renal Clinical Trialists, Nancy, France. ³Robertson Centre for Biostatistics, Institute of Health and Wellbeing, University of Glasgow, Glasgow Royal Infirmary, Glasgow G12 8QQ, UK. ⁴Department of Cardiology, Maastricht University Medical Center (MUMC+), PO Box 5800, 6202 AZ Maastricht, The Netherlands. ⁵Department of Cardiology, University of Hull, Castle Hill Hospital, Cottingham, East Riding of Yorkshire, UK. ⁶Department of Cardiology, Cortona Hospital, Arezzo, Italy. ⁷Department of Internal Medicine and Cardiology, Campus Virchow Klinikum, Charité' University Medicine Berlin, Berlin Institute of Health (BIH), and German Centre for Cardiovascular research (DZHK), Partner Site Berlin, Berlin, Germany.

Received: 3 June 2021 Accepted: 29 July 2021

Published online: 09 August 2021

## References

1. Domanski M, Krause-Steinrauf H, Deedwania P, Follmann D, Ghali JK, Gilbert E, Haffner S, Katz R, Lindenfeld J, Lowes BD, et al. The effect of diabetes on outcomes of patients with advanced heart failure in the BEST trial. J Am Coll Cardiol. 2003;42(5):914–22.

2. Beijer K, Nowak C, Sundström J, Årnlöv J, Fall T, Lind L. In search of causal pathways in diabetes: a study using proteomics and genotyping data from a cross-sectional study. Diabetologia. 2019;62(11):1998–2006.

3. Molvin J, Pareek M, Jujic A, Melander O, Råstam L, Lindblad U, Daka B, Leósdóttir M, Nilsson PM, Olsen MH, et al. Using a targeted proteomics chip to explore pathophysiological pathways for incident diabetes—the Malmö preventive project. Sci Rep. 2019;9(1):272.

4. Tromp J, Voors AA, Sharma A, Ferreira JP, Ouwerkerk W, Hillege HL, Gomez KA, Dickstein K, Anker SD, Metra M, et al. Distinct pathological pathways in patients with heart failure and diabetes. JACC Heart Fail. 2020;8(3):234–42.

5. Pellicori P, Ferreira JP, Mariottoni B, Brunner-La Rocca HP, Ahmed FZ, Verdonschot J, Collier T, Cuthbert JJ, Petutschnigg J, Mujaj B, et al. Effects of spironolactone on serum markers of fibrosis in people at high risk of developing heart failure: rationale, design and baseline characteristics of a proof-of-concept, randomised, precision-medicine, prevention trial. The Heart OMics in AGing (HOMAGE) trial. Eur J Heart Fail. 2020;22(9):1711–23.

6. Cleland JGF, Ferreira JP, Mariottoni B, Pellicori P, Cuthbert J, Verdonschot JA, Petutschnigg J, Ahmed FZ, Cosmi F, Brunner La Rocca H-P, et al. The effect of spironolactone on cardiovascular function and markers of fibrosis in people at increased risk of developing heart failure: the heart 'OMics' in AGEing (HOMAGE) randomized clinical trial. Eur Heart J. 2020;42(6):684–96.

7. Davies JI, Band M, Morris A, Struthers AD. Spironolactone impairs endothelial function and heart rate variability in patients with type 2 diabetes. Diabetologia. 2004;47(10):1687–94.

8. Yamaji M, Tsutamoto T, Kawahara C, Nishiyama K, Yamamoto T, Fujii M, Horie M. Effect of eplerenone versus spironolactone on cortisol and
hemoglobin A1(c) levels in patients with chronic heart failure. Am Heart J. 2010;160(5):915-21.

9. Ferreira JP, Verdonschot J, Wang P, Pizard A, Collier T, Ahmed FZ, Brunner-La-Rocca H-P, Clark AL, Cosmi F, Cuthbert J, et al. Proteomic and mechanical analysis of spironolactone in patients at risk for HF. JACC Heart Fail. 2021;9(4):268-77.

10. Ago T, Sadoshima J. GDF15, a cardioprotective TGF-beta superfamily protein. Circ Res. 2006;98(3):294-7.

11. Sharma A, Stevens SR, Lucas J, Fiuzat M, Adams KF, Whellan DJ, Donahue MP, Kitzman DW, Piña IL, Zannad F, et al. Utility of growth differentiation factor-15, a marker of oxidative stress and inflammation, in chronic heart failure: insights from the HF-ACTION study. JACC Heart Fail. 2017;5(10):724-34.

12. Day EA, Ford RJ, Smith BK, Mohammadi-Shemirani P, Morrow MR, Gutgesell RM, Lu R, Raphenya AR, Kabiri M, McArthur AG, et al. Metformin-induced increases in GDF15 are important for suppressing appetite and promoting weight loss. Nat Metab. 2019;1(12):1202-8.

13. Coll AP, Chen M, Taskar P, Rimmington D, Patel S, Tadross JA, Cimino I, Yang M, Welsh P, Virtue S, et al. GDF15 mediates the effects of metformin on body weight and energy balance. Nature. 2020;578(7795):444-8.

14. Dominguez-Rodriguez A, Abreu-Gonzalez P, Avanzas P. Usefulness of growth differentiation factor-15 levels to predict diabetic cardiomyopathy in asymptomatic patients with type 2 diabetes mellitus. Am J Cardiol. 2014;114(6):890-4.

15. Tromp J, Westenbrink BD, Ouwerkerk W, van Veldhuisen DJ, Samani NJ, Ponikowski P, Metra M, Anker SD, Cleland JG, Dickstein K, et al. Identifying pathophysiological mechanisms in heart failure with reduced versus preserved ejection fraction. J Am Coll Cardiol. 2018;72(10):1081-90.

16. Feng X, Guan D, Auen T, Choi JW, Salazar Hernández MA, Lee J, Chun H, Faruk F, Kaplun E, Herbert Z, et al. IL1R1 is required for celastrol's leptin-sensitization and antiobesity effects. Nat Med. 2019;25(4):575-82.

17. Liu J, Ibi D, Taniguchi K, Lee J, Herrema H, Akosman B, Mucka P, Salazar Hernandez MA, Uyar MF, Park SW, et al. Inflammation improves glucose homeostasis through IKKβ-XBP1s interaction. Cell. 2016;167(4):1052-1066. e1018.

18. Wernstedt Asterholm I, Tao C, Morley TS, Wang QA, Delgado-Lopez F, Wang ZV, Scherer PE. Adipocyte inflammation is essential for healthy adipose tissue expansion and remodeling. Cell Metab. 2014;20(1):103-18.

19. Hoes MF, Tromp J, Ouwerkerk W, Bomer N, Oberdorf-Maass SU, Samani NJ, Ng LL, Lang CC, van der Harst P, Hillege H, et al. The role of cathepsin D in the pathophysiology of heart failure and its potentially beneficial properties: a translational approach. Eur J Heart Fail. 2019;22(11):2102-11.

20. Ding L, Goossens GH, Oligschlaeger Y, Houben T, Blaak EE, Shiri-Sverdlov R. Plasma cathepsin D activity is negatively associated with hepatic insulin sensitivity in overweight and obese humans. Diabetologia. 2020;63(2):374-84.

21. Eguchi A, Feldstein AE. Lysosomal cathepsin D contributes to cell death during adipocyte hypertrophy. Adipocyte. 2013;2(3):170-5.

22. Tryggestad JB, Shah RD, Braffett BH, Bacha F, Gidding SS, Gubitosi-Klug RA, Shah AS, Urbina EM, Levitt Katz LE. Circulating adhesion molecules and associations with HbA1c, hypertension, nephropathy, and retinopathy in the treatment options for type 2 diabetes in adolescent and youth study. Pediatr Diabetes. 2020;21(6):923-31.

23. Zemlin AE, Matsha TE, Kengne AP, Hon GM, Erasmus RT. Correlation of E-selectin concentrations with carotid intima-media thickness and cardio-metabolic profile of mixed ancestry South Africans: a cross-sectional study. Ann Clin Biochem. 2017;54(1):92-100.

24. Sass C, Pallaud C, Zannad F, Visvikis S. Relationship between E-selectin L/F554 polymorphism and blood pressure in the Stanislas cohort. Hum Genet. 2000;107(1):58-61.

25. Matsubara T, Mita A, Minami K, Hosooka T, Kitazawa S, Takahashi K, Tamori Y, Yokoi N, Watanabe M, Matsuo E, et al. PGRN is a key adipokine mediating high fat diet-induced insulin resistance and obesity through IL-6 in adipose tissue. Cell Metab. 2012;15(1):38-50.

**Publisher's Note**

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Ready to submit your research? Choose BMC and benefit from:

* fast, convenient online submission

* thorough peer review by experienced researchers in your field

* rapid publication on acceptance

* support for research data, including large and complex data types

* gold Open Access which fosters wider collaboration and increased citations

* maximum visibility for your research: over 100M website views per year

At BMC, research is always in progress.

Learn more biomedcentral.com/submissions

BMC