# Identification of patient subtypes based on protein expression for prediction of heart failure after myocardial infarction

Wilfried Heyse, Vincent Vandewalle, Guillemette Marot, Philippe Amouyel, Christophe Bauters, Florence Pinet

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

## To cite this version:

Wilfried Heyse, Vincent Vandewalle, Guillemette Marot, Philippe Amouyel, Christophe Bauters, et al. Identification of patient subtypes based on protein expression for prediction of heart failure after myocardial infarction. iScience, 2023, 26 (3), pp.106171. 10.1016/j.isci.2023.106171. inserm-04191363

HAL Id: inserm-04191363

https://inserm.hal.science/inserm-04191363

Submitted on 30 Aug 2023

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_2.png" /></figure>
Article

# Identification of patient subtypes based on protein expression for prediction of heart failure after myocardial infarction

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

Wilfried Heyse,
Vincent
Vandewalle,
Guillemette
Marot, Philippe
Amouyel,
Christophe
Bauters, Florence
Pinet

florence.pinet@pasteur-lille.fr

## Highlights

We quantified 4,668 proteins in two cohorts of patients with a first myocardial infarction

Fifty proteins were significantly associated with the occurrence of heart failure

A clustering method identified two groups of patients based on the 50 proteins

Group 2 of patients was associated with a high risk of occurrence of heart failure
## Article

# Identification of patient subtypes based on protein expression for prediction of heart failure after myocardial infarction

Wilfried Heyse,<fnref n="1" />,<fnref n="3" /> Vincent Vandewalle,<fnref n="2" />,<fnref n="3" /> Guillemette Marot,<fnref n="2" />,<fnref n="3" /> Philippe Amouyel,<fnref n="1" /> Christophe Bauters,<fnref n="1" /> and Florence Pinet<fnref n="1" />,<fnref n="4" />,*

## SUMMARY

This study investigates the ability of high-throughput aptamer-based platform to identify circulating biomarkers able to predict occurrence of heart failure (HF), in blood samples collected during hospitalization of patients suffering from a first myocardial infarction (MI). REVE-1 (derivation) and REVE-2 (validation) cohorts included respectively 254 and 238 patients, followed up respectively 9 · 2 ± 4 · 8 and 7 · 6 ± 3 · 0 years. A blood sample collected during hospitalization was used for quantifying 4,668 proteins. Fifty proteins were significantly associated with long-term occurrence of HF with all-cause death as the competing event. k-means, an unsupervised clustering method, identified two groups of patients based on expression levels of the 50 proteins. Group 2 was significantly associated with a higher risk of HF in both cohorts. These results showed that a subset of 50 selected proteins quantified during hospitalization of MI patients is able to stratify and predict the long-term occurrence of HF.

## INTRODUCTION

Despite significant therapeutic improvements during the last decades, the long-term risk of heart failure (HF) after myocardial infarction (MI) remains significant and ischemic HF is a major cause of mortality worldwide.<fnref n="1" /> Post-MI HF is too often diagnosed at a late stage when its irreversible consequences are already established. The estimation of the risk of HF in the early post-MI period currently relies on clinical variables, left ventricular function parameters, and conventional cardiac biomarkers such as troponin or B-type natriuretic peptide (BNP). However other cardiovascular biomarkers might reflect the activation of new potential pathways after MI and contribute to HF in both the short and long term.<fnref n="2" /> Broader approaches are important for clinicians to understand the implications of different pathobiological axes. Despite the proliferation of candidate biomarkers, there is limited data comparing comprehensively the prognostic value of biomarkers when assessed in large arrays.

Explorative analysis using large scale protein measurement methods allow the simultaneous analysis of large biomarker panels. Thanks to SOMA(Slow Off-rate Modified Aptamers)scan assay, over 5,000 proteins can be measured covering 8 logs of abundance in the human proteome. These proteins are not targeted toward any particular disease and, thus, such wide panels of proteins can help to discover new biomarkers using appropriate statistical methods to analyze these data. Very recently, aptamer assays were shown to provide excellent precision and an unprecedented coverage and promise for disease associations.<fnref n="3" /> Using SOMAscan profiling (4,453 targets), Gui et al. showed that plasma multiprotein score improved risk stratification in patients with HF and reduced ejection fraction and identified novel candidates.<fnref n="4" />

The aim of this work was to use circulating plasma proteins expression levels quantified during the hospitalization of patients suffering from MI to identify groups of patients able to predict the long-term risk of occurrence of HF after MI. The prospective REVE-1 (REmodelage VEntriculaire) cohort⁵ and REVE-2 cohort⁶ in which patients with a first MI were included and underwent long-term follow-up⁷ were used respectively as derivation and validation cohorts. A large panel of 4,668 proteins was measured with the SOMAscan assay in blood sample collected during initial hospitalization. To build groups of patients able to predict the occurrence of HF, the analysis was divided into three steps (Figure 1). First, a protein selection step was performed to focus on the relevant proteins in the derivation cohort, REVE-1. Second, the selected

<footnote n="1">¹UniversityLille, Inserm, CHU Lille, Institut Pasteur de Lille, U1167 - RID-AGE - Facteurs de risque et déterminants moléculaires des maladies liées au vieillissement, 59000 Lille, France</footnote>

<footnote n="2">²UniversityLille, CHU Lille, URL2694-METRICS - Evaluation des technologies de santé et des pratiques médicales, 59000 Lille, France</footnote>

<footnote n="3">³Inria, Modal, 59000 Lille, France</footnote>

<footnote n="4">⁴Lead contact</footnote>

<footnote>*Correspondence: florence.pinet@pasteur-lille.fr</footnote>

<footnote>https://doi.org/10.1016/j.isci.2023.106171</footnote>
<figure><img src="image_4.png" /><figcaption>**Figure 1. Overview of the study**</figcaption></figure>

REVE-1 study⁵ was used as derivation cohort and REVE-2 study⁶ as validation cohort. Proteomic data analysis was performed on all plasma samples collected during hospitalization of patients from both cohorts by SOMAscan assay (version V4.0). The SOMAScan platform measured accurately 4668 plasma proteins. Log2 transformed data were centered and reduced. Standardization parameters were calculated with the data from the REVE-1 cohort (see Table S2) and then applied to both REVE-1 and REVE-2 cohorts. Univariate competing risk models were fitted for each protein, and significance tests selected 50 proteins to be associated with occurrence of hospitalization for HF (p< 1 · 07 10⁻⁵ in accordance with Bonferroni's method). Clustering was then performed on the REVE-1 patients with the 50 selected proteins, by k-means procedure with k = 2 groups. Patients of REVE-2 were then assigned to one of the groups identified on REVE-1. A competing risk model was then developed based on the group information on both cohorts.

proteins were used to build groups of patients on REVE-1 using k-means, an unsupervised clustering algo-
rithm. Patients of REVE-2, the validation cohort, were assigned to the groups built on REVE-1. Third, the
groups' predictive ability for occurrence of hospitalization for HF was assessed using competing risk
models in both cohorts. In addition, Ingenuity Pathway Analysis (IPA) and Gene Ontology (GO) analysis
were performed to define molecular networks enriched from the selected proteins.

Evidence before this study

HF following an MI is too often diagnosed too late when its irreversible consequences are established. The estimation of the risk of HF in the early post-MI period currently relies on clinical variables, left ventricular function parameters, and conventional cardiac biomarkers such as troponin or BNP. Despite proliferation of candidate biomarkers, there is limited data comparing comprehensively the prognostic value of biomarkers.

Added value of this study

We performed a discovery proteomics approach by quantification of 4,468 proteins in two cohorts of pa-
tients with a first MI, REVE1 (derivation cohort) and REVE 2 (validation cohort). A total of 50 proteins were
selected to be significantly associated with the occurrence of hospitalization for HF. An unsupervised clus-
tering method identified 2 groups of patients based on the expression levels of the 50 proteins in the
REVE1 cohort that were validated in the REVE2 cohort. Differences in protein expression led to identifying
<table><caption>Table 1. Baseline characteristics of the patients included in the Derivation (REVE-1) and Validation (REVE-2) cohorts</caption><thead><tr><th></th><th>REVE-1 (n = 254)</th><th>REVE-2 (n = 238)</th><th>p value</th></tr></thead><tbody><tr><td>Age (years)</td><td>58 ± 14</td><td>57 ± 14</td><td>0.28</td></tr><tr><td>Women, n (%)</td><td>66 (26)</td><td>46 (19)</td><td>0.23</td></tr><tr><td>Body mass index (kg/m²)</td><td>27 ± 1 ± 4 ± 8</td><td>27 ± 2 ± 4 ± 7</td><td>0.83</td></tr><tr><td>Hypercholesterolemia, n (%)</td><td>119 (47)</td><td>79 (33)</td><td>0.48</td></tr><tr><td>History of hypertension, n (%)</td><td>114 (45)</td><td>86 (36)</td><td>0.46</td></tr><tr><td>Smoking, n (%)</td><td>173 (68)</td><td>164 (69)</td><td>0.68</td></tr><tr><td>Diabetes mellitus, n (%)</td><td>58 (23)</td><td>51 (21)</td><td>0.19</td></tr><tr><td>Initial reperfusion therapy, n (%)</td><td></td><td></td><td></td></tr><tr><td>Primary PCI</td><td>76 (30)</td><td>124 (52)</td><td>&lt;0.001</td></tr><tr><td>Thrombolysis</td><td>133 (52)</td><td>85 (36)</td><td>&lt;0.001</td></tr><tr><td>None</td><td>45 (18)</td><td>29 (12)</td><td>0.11</td></tr><tr><td>Multivessel CAD, n (%)</td><td>96 (39)</td><td>94 (41)</td><td>0.39</td></tr><tr><td>PCI during hospitalization, n (%)</td><td>226 (89)</td><td>205 (86)</td><td>0.91</td></tr><tr><td>Final TIMI grade 3 flow in infarct-related vessel, n (%)</td><td>211 (84)</td><td>206 (90)</td><td>0.84</td></tr><tr><td>Systolic blood pressure (mmHg)</td><td>112 ± 16</td><td>110 ± 15</td><td>0.24</td></tr><tr><td>Diastolic blood pressure (mmHg)</td><td>66 ± 12</td><td>63 ± 11</td><td>&lt;0.001</td></tr><tr><td>Heart rate (bp.m.)</td><td>68 ± 12</td><td>72 ± 14</td><td>0.01</td></tr><tr><td>Killip class ≥2, n (%)</td><td>70 (28)</td><td>76 (32)</td><td>0.25</td></tr><tr><td>Peak creatine kinase (IU/L)</td><td>2558 [1252 to 4013]</td><td>2348 [1450 to 4192]</td><td>0.48</td></tr><tr><td>End-systolic volume (mL/m²)</td><td>29 ± 11</td><td>27 ± 10</td><td>0.01</td></tr><tr><td>End Diastolic Volume (mL/m²)</td><td>57 ± 15</td><td>52 ± 14</td><td>&lt;0.001</td></tr><tr><td>Left ventricular ejection fraction (%)</td><td>49 ± 10</td><td>49 ± 8</td><td>0.58</td></tr><tr><td>Wall motion score index</td><td>1.87 ± 0.15</td><td>1.91 ± 0.15</td><td>&lt;0.001</td></tr><tr><td>Treatment at discharge, n (%)</td><td></td><td></td><td></td></tr><tr><td>Antiplatelet therapy</td><td>247 (98)</td><td>235 (99)</td><td>1</td></tr><tr><td>β-blockers</td><td>238 (94)</td><td>230 (97)</td><td>1</td></tr><tr><td>ACE-I or ARB</td><td>244 (96)</td><td>230 (97)</td><td>1</td></tr><tr><td>Aldosterone antagonists</td><td>34 (13)</td><td>78 (33)</td><td>0.05</td></tr><tr><td>Statins</td><td>248 (98)</td><td>225 (95)</td><td>1</td></tr></tbody></table>

Data are presented as mean ± SD frequency (percentage) or median [Q1 to Q3].

bpm, beats per minute; CAD, coronary artery disease; PCI, percutaneous coronary intervention; TIMI, thrombolysis in myocardial infarction; ACE-I, angiotensin-converting enzyme inhibitors; ARB, angiotensin II receptor blockers. Imputed data are not shown in this table.

key physiopathological processes combining differences in molecules leading to two groups of patients, low (group 1) and high (group 2) risk of long-term adverse cardiac outcomes following MI.

## Implications of the available evidence

Based on the expression of the 50 proteins, group 2 of patients was associated with a high risk of occurrence of HF. Stratification of MI patients based on proteins involved in cell-to-cell communication should be important in the future.

## RESULTS

### Study populations

Patients in both cohorts had similar age and gender; indicators of MI size (wall motion score index and peak creatine kinase) and Killip class were also similar (Table 1). Patients in the REVE-2 cohort were more often treated by primary PCI and less often by thrombolysis. Statistical differences were observed between both cohorts for diastolic blood pressure, heart rate, end systolic volume, end diastolic volume and wall motion
<table><caption>Table 2. List of plasma proteins selected to be associated with occurrence of heart failure in REVE-1 cohort</caption><thead><tr><th colspan="2">REVE-1</th><th colspan="2">REVE-2</th><th rowspan="2">UniProt ID</th><th rowspan="2">Symbol</th><th rowspan="2">Entrez gene name</th></tr><tr><th>p value<sup>a</sup></th><th>SHR#</th><th>p value<sup>a</sup></th><th>SHR#</th></tr></thead><tbody><tr><td>&lt;1.0e-16</td><td>1.18</td><td>3.1e-03</td><td>2.59</td><td>C9JXX5</td><td>C11ORF94</td><td>Uncharacterized protein C11orf94</td></tr><tr><td>3.8e-12</td><td>1.17</td><td>5.0e-01</td><td>0.71</td><td>P27815</td><td>PDE4A</td><td>cAMP-specific 3',5'-cyclic phosphodiesterase 4A</td></tr><tr><td>9.0e-11</td><td>1.31</td><td>9.8e-01</td><td>0.99</td><td>Q9Y3Z3</td><td>SAMHD1</td><td>Deoxynucleoside triphosphate triphosphohydrolase SAMHD1</td></tr><tr><td>2.0e-10</td><td>1.16</td><td>8.6e-01</td><td>0.97</td><td>Q9Y644</td><td>RFNG</td><td>Beta-1,3-N-acetylglucosaminyltransferase radical fringe</td></tr><tr><td>2.4e-10</td><td>0.41</td><td>1.1e-03</td><td>0.50</td><td>O75473</td><td>LGR5</td><td>Leucine-rich repeat-containing G-protein coupled receptor 5</td></tr><tr><td>1.5e-08</td><td>1.12</td><td>4.9e-01</td><td>1.19</td><td>Q16568</td><td>CARTPT</td><td>Cocaine- and amphetamine-regulated transcript protein</td></tr><tr><td>1.8e-08</td><td>0.52</td><td>4.8e-04</td><td>0.57</td><td>P56704</td><td>WNT3A</td><td>Protein Wnt-3a</td></tr><tr><td>2.6e-08</td><td>0.60</td><td>1.5e-02</td><td>0.65</td><td>P05543</td><td>SERPINA7</td><td>Thyroxine-binding globulin</td></tr><tr><td>3.6e-08</td><td>0.57</td><td>5.8e-02</td><td>0.85</td><td>P19823</td><td>ITIH2</td><td>Inter-alpha-trypsin inhibitor heavy chain H2</td></tr><tr><td>4.0e-08</td><td>1.54</td><td>3.3e-02</td><td>1.21</td><td>Q8TAT2</td><td>FGFBP3</td><td>Fibroblast growth factor-binding protein 3</td></tr><tr><td>4.8e-08</td><td>1.15</td><td>3.2e-04</td><td>2.18</td><td>Q92831</td><td>KAT2B</td><td>Histone acetyltransferase KAT2B</td></tr><tr><td>9.8e-08</td><td>1.83</td><td>1.6e-04</td><td>2.06</td><td>Q9Y251</td><td>HPSE</td><td>Heparanase</td></tr><tr><td>1.5e-07</td><td>1.32</td><td>2.0e-03</td><td>1.25</td><td>A8K7I4</td><td>CLCA1</td><td>Calcium-activated chloride channel regulator 1</td></tr><tr><td>2.7e-07</td><td>0.52</td><td>1.4e-03</td><td>0.47</td><td>Q9NTU7</td><td>CBLN4</td><td>Cerebellin-4</td></tr><tr><td>3.6e-07</td><td>1.78</td><td>2.4e-04</td><td>1.58</td><td>P61769</td><td>B2M</td><td>Beta-2-microglobulin</td></tr><tr><td>3.7e-07</td><td>1.48</td><td>9.1e-05</td><td>1.73</td><td>Q2I0M5</td><td>RSPO4</td><td>R-spondin-4</td></tr><tr><td>5.2e-07</td><td>1.33</td><td>4.3e-02</td><td>1.36</td><td>P52943</td><td>CRIP2</td><td>Cysteine-rich protein 2</td></tr><tr><td>9.5e-07</td><td>0.59</td><td>4.7e-01</td><td>0.87</td><td>P33764</td><td>S100A3</td><td>Protein S100-A3</td></tr><tr><td>1.3e-06</td><td>1.50</td><td>3.6e-01</td><td>1.22</td><td>Q9BUP3</td><td>HTATIP2</td><td>Oxidoreductase HTATIP2</td></tr><tr><td>1.4e-06</td><td>0.40</td><td>7.9e-05</td><td>0.41</td><td>Q96PQ0</td><td>SORCS2</td><td>VPS10 domain-containing receptor SorCS2</td></tr><tr><td>1.4e-06</td><td>0.55</td><td>3.5e-02</td><td>0.68</td><td>Q8IWT1</td><td>SCN4B</td><td>Sodium channel subunit beta-4</td></tr><tr><td>1.5e-06</td><td>0.48</td><td>8.8e-04</td><td>0.59</td><td>Q6P988</td><td>NOTUM</td><td>Palmitoleoyl-protein carboxylesterase NOTUM</td></tr><tr><td>1.7e-06</td><td>0.57</td><td>2.9e-05</td><td>0.54</td><td>P29622</td><td>SERPINA4</td><td>Kallistatin</td></tr><tr><td>1.9e-06</td><td>1.79</td><td>6.7e-05</td><td>1.75</td><td>P01034</td><td>CST3</td><td>Cystatin-C</td></tr><tr><td>1.9e-06</td><td>2.71</td><td>5.2e-01</td><td>1.21</td><td>Q8N6C8</td><td>LILRA3</td><td>Leukocyte immunoglobulin-like receptor subfamily A member 3</td></tr><tr><td>2.0e-06</td><td>0.23</td><td>6.4e-01</td><td>1.05</td><td>Q9UNP4</td><td>ST3GAL5</td><td>Lactosylceramide alpha-2,3-sialyltransferase</td></tr><tr><td>2.7e-06</td><td>0.50</td><td>1.5e-07</td><td>0.42</td><td>Q8N6M8</td><td>IQCF1</td><td>IQ domain-containing protein F1</td></tr><tr><td>3.1e-06</td><td>1.28</td><td>3.2e-02</td><td>1.60</td><td>P0DJ93</td><td>SMIM13</td><td>Small integral membrane protein 13</td></tr><tr><td>3.2e-06</td><td>0.68</td><td>7.1e-02</td><td>0.89</td><td>P01008</td><td>SERPINC1</td><td>Antithrombin-III</td></tr><tr><td>3.4e-06</td><td>1.87</td><td>4.1e-06</td><td>2.86</td><td>Q8TE54</td><td>SLC26A7</td><td>Anion exchange transporter</td></tr><tr><td>3.4e-06</td><td>0.41</td><td>8.0e-01</td><td>0.87</td><td>Q6ZUB0</td><td>SPATA31D4</td><td>Putative spermatogenesis-associated protein 31D4</td></tr><tr><td>4.3e-06</td><td>1.83</td><td>5.6e-04</td><td>1.72</td><td>P16860</td><td>NPPB</td><td>Natriuretic peptides B</td></tr><tr><td>4.4e-06</td><td>1.45</td><td>5.1e-03</td><td>1.46</td><td>Q12805</td><td>EFEMP1</td><td>EGF-containing fibulin-like extracellular matrix protein 1</td></tr><tr><td>4.8e-06</td><td>0.55</td><td>5.7e-01</td><td>1.04</td><td>Q5H8A3</td><td>NMS</td><td>Neuromedin-S</td></tr><tr><td>5.0e-06</td><td>0.63</td><td>1.0e-04</td><td>0.65</td><td>P35858</td><td>IGFALS</td><td>Insulin-like growth factor-binding protein complex acid labile subunit</td></tr></tbody></table>

(Continued on next page)
<table><caption>**Table 2. Continued**</caption><thead><tr><th colspan="2">REVE-1</th><th colspan="2">REVE-2</th><th rowspan="2">UniProt ID</th><th rowspan="2">Symbol</th><th rowspan="2">Entrez gene name</th></tr><tr><th>p value<sup>a</sup></th><th>SHR#</th><th>p value<sup>a</sup></th><th>SHR#</th></tr></thead><tbody><tr><td>5.1e-06</td><td>1.39</td><td>5.4e-02</td><td>1.29</td><td>Q9NS68</td><td>TNFRSF19</td><td>Tumor necrosis factor receptor superfamily member 19</td></tr><tr><td>5.9e-06</td><td>0.49</td><td>1.9e-03</td><td>0.48</td><td>O60260</td><td>PARK2</td><td>E3 ubiquitin-protein ligase parkin</td></tr><tr><td>5.9e-06</td><td>1.34</td><td>2.5e-02</td><td>1.62</td><td>Q13522</td><td>PPP1R1A</td><td>Protein phosphatase 1 regulatory subunit 1A</td></tr><tr><td>6.4e-06</td><td>0.52</td><td>7.7e-02</td><td>0.73</td><td>P55285</td><td>CDH6</td><td>Cadherin-6</td></tr><tr><td>6.6e-06</td><td>0.58</td><td>2.7e-01</td><td>0.76</td><td>P45985</td><td>MAP2K4</td><td>Dual specificity mitogen-activated protein kinase kinase 4</td></tr><tr><td>6.7e-06</td><td>0.54</td><td>3.0e-07</td><td>0.46</td><td>Q2Y0W8</td><td>SLC4A8</td><td>Electroneutral sodium bicarbonate exchanger 1</td></tr><tr><td>6.8e-06</td><td>0.69</td><td>5.9e-02</td><td>0.76</td><td>P01042</td><td>KNG1</td><td>Kininogen, HMW, Two Chain</td></tr><tr><td>7.6e-06</td><td>0.49</td><td>2.5e-08</td><td>0.23</td><td>Q8N3Z0</td><td>PRSS35</td><td>Inactive serine protease 35</td></tr><tr><td>7.8e-06</td><td>0.51</td><td>4.0e-01</td><td>1.11</td><td>P27487</td><td>DPP4</td><td>Dipeptidylpeptidase 4</td></tr><tr><td>8.1e-06</td><td>0.55</td><td>5.6e-02</td><td>0.56</td><td>Q8IYN2</td><td>TCEAL8</td><td>Transcription elongation factor A protein-like 8</td></tr><tr><td>8.5e-06</td><td>1.75</td><td>2.5e-01</td><td>1.21</td><td>Q15848</td><td>ADIPOQ</td><td>Adiponectin</td></tr><tr><td>8.9e-06</td><td>0.50</td><td>4.2e-03</td><td>0.52</td><td>Q9UII4</td><td>HERC5</td><td>E3 ISG15-protein ligase HERC5</td></tr><tr><td>9.4e-06</td><td>1.86</td><td>1.7e-04</td><td>1.77</td><td>P07998</td><td>RNASE1</td><td>RNase pancreatic</td></tr><tr><td>9.5e-06</td><td>0.44</td><td>4.7e-06</td><td>0.39</td><td>Q9UBC7</td><td>GALP</td><td>Galanin-like peptide</td></tr><tr><td>9.8e-06</td><td>0.54</td><td>1.5e-04</td><td>0.49</td><td>Q9UK76</td><td>JTP1</td><td>Jupiter microtubule associated homolog 1</td></tr></tbody></table>

UniProt ID, protein symbol, full protein name are provided from the UniProtKB database (https://www.uniprot.org). The proteins in italics were quantified by several SOMAmers and the mean value was used. See also Tables S1 and S2-S6, Figures S1 and S2.

systolic index. The mean follow-up was 9.2 ± 4.8 years in REVE-1 and 7.6 ± 3.0 years in REVE-2. One patient was lost during the follow-up leading to 254 patients in REVE-1 for the analysis; no patients in REVE-2 were lost during the follow-up. The numbers of patients who reached the primary endpoint (hospitalization for HF) during follow-up were respectively 49 in REVE-1 and 28 in REVE-2. The numbers of patients who reached the competing event (death from all causes) during the long term follow-up were respectively 63 in REVE-1 and 26 in REVE-2.

## Selection of proteins

Standardization parameters for each protein and each patient were calculated using the mean and standard deviation (SD) calculated with the data from REVE-1 cohort (Table S2). The proteomic data of REVE-2 were then standardized protein by protein for each patient using these values following the calculation:((value of the protein) – (the mean of protein in REVE-1))/(the corresponding SD in REVE-1).

After standardization, 50 proteins were selected to be significantly associated with the outcome event in REVE-1 (listed in Table 2 with their significance levels and subhazard ratios and in Table S3 for their location and type). Among the 50 proteins selected on REVE-1, 44 proteins on REVE-2 were regulated in the same manner for the outcome of REVE-2 patients and 6 were conversely but not significantly regulated (PDE4A, SAMHD1, RFNG, ST3GAL5, NMS and DPP4) (Table 2). The correlation heatmap among the 50 selected proteins in REVE-1 shows that some of the proteins are highly correlated, with the highest correlation of 0.893 between B2M and CST3 (Figure S2). This shows that several proteins carry a close tendency, which will be taken into account by the following clustering approach.

## Patients groups

The silhouette criteria led to use k-means procedure with k = 2 groups on REVE-1 (not shown). Patients were split as 160 patients in group 1 and 94 in group 2 and the two groups show contrasted protein expression profiles (Figure 2A).

<footnote><sup>a</sup>Raw p-values are indicated, significance threshold was set at 1.07 × 10⁻⁵ in accordance with Bonferroni's method; # SHR: Subhazard ratio was calculated by a competing risk model.</footnote>
<figure><img src="image_5.png" /><figcaption>**Figure 2. Heatmap visualization of the two identified groups on REVE-1 (A) and REVE-2 (B) cohorts**</figcaption></figure>

Columns represent the 50 selected proteins significantly modulated according to the outcome of patients and classified in the same order for both cohorts. Lines represent the patients divided in the two groups (G1 and G2) identified by clustering with the k-means. Cells are colored according to protein abundance expressed in RFU (Relative Fluorescence Unit). Red represents high abundance whereas blue indicates low abundance. Protein levels are log2-transformed and standardized. See also Table S4.

These contrasted protein expression profiles were also found in REVE-2 cohort. Patients were mapped to one of the two groups identified on REVE-1 with respectively 190 patients assigned to group 1 and 48 to group 2 and similar characteristics were observed, with opposite protein expression profiles between the two groups (Figure 2B).

We also observed significant clinical differences between the patients of the two groups in both cohorts with significant differences in REVE-1 for history of hypertension, diabetes, initial reperfusion therapy, multi-vessel coronary artery disease (CAD), final thrombolysis in myocardial infarction (TIMI) grade 3 flow in infarct-related vessel, heart rate, Killip class over 2, end-systolic volume, left ventricular ejection fraction, WMSI and for treatments at discharge (aldosterone antagonists). In REVE-2, gender, smoking and treatments at discharge (β-blockers) were significantly different between the identified groups. Significant differences for the age were found in both cohorts. Thus, protein based clustering leads to the identification of significant differences in the clinical characteristics of patients (Table 3).
<table><caption>Table 3. Clinical characteristics of patients between the two identified groups on both cohorts</caption><thead><tr><th rowspan="2">Variables</th><th colspan="3">REVE-1 (Derivation) (n = 254)</th><th colspan="3">REVE-2 (Validation) (n = 238)</th></tr><tr><th>Group 1 (n = 160)</th><th>Group 2 (n = 94)</th><th>p value</th><th>Group 1 (n = 190)</th><th>Group 2 (n = 48)</th><th>p value</th></tr></thead><tbody><tr><td>Age (years)</td><td>54 ± 13</td><td>65 ± 13</td><td>&lt;0.001</td><td>54 ± 13</td><td>68 ± 12</td><td>&lt;0.001</td></tr><tr><td>Women, n (%)</td><td>46 (29)</td><td>20 (21)</td><td>0.245</td><td>30 (16)</td><td>16 (33)</td><td>0.011</td></tr><tr><td>Body mass index (kg/m²)</td><td>26.7 ± 4.7</td><td>27.6 ± 4.9</td><td>0.146</td><td>27.4 ± 4.4</td><td>26.1 ± 5.4</td><td>0.071</td></tr><tr><td>Hypercholesterolemia, n (%)</td><td>73 (46)</td><td>46 (49)</td><td>0.704</td><td>62 (33)</td><td>17 (35)</td><td>0.846</td></tr><tr><td>History of hypertension, n (%)</td><td>55 (34)</td><td>59 (63)</td><td>&lt;0.001</td><td>63 (33)</td><td>23 (48)</td><td>0.083</td></tr><tr><td>Smoking, n (%)</td><td>110 (69)</td><td>63 (67)</td><td>0.884</td><td>140 (74)</td><td>24 (50)</td><td>0.003</td></tr><tr><td>Diabetes mellitus, n (%)</td><td>29 (18)</td><td>29 (31)</td><td>0.029</td><td>36 (19)</td><td>15 (31)</td><td>0.097</td></tr><tr><td>No initial reperfusion therapy, n (%)</td><td>22 (14)</td><td>23 (24)</td><td>0.047</td><td>19 (10)</td><td>10 (21)</td><td>0.071</td></tr><tr><td>Multivessel CAD, n (%)</td><td>49 (31)</td><td>47 (52)</td><td>0.002</td><td>71 (38)</td><td>23 (51)</td><td>0.165</td></tr><tr><td>PCI during hospitalization, n (%)</td><td>144 (90)</td><td>82 (87)</td><td>0.637</td><td>165 (87)</td><td>40 (83)</td><td>0.693</td></tr><tr><td>Final TIMI grade 3 flow in infarct-related vessel, n (%)</td><td>143 (90)</td><td>68 (74)</td><td>0.002</td><td>166 (90)</td><td>40 (91)</td><td>1</td></tr><tr><td>Systolic blood pressure (mmHg)</td><td>113 ± 16</td><td>110 ± 16</td><td>0.174</td><td>111 ± 14</td><td>109 ± 20</td><td>0.664</td></tr><tr><td>Diastolic blood pressure (mmHg)</td><td>67 ± 12</td><td>65 ± 12</td><td>0.183</td><td>63 ± 10</td><td>60 ± 15</td><td>0.079</td></tr><tr><td>Heart rate (bpm.)</td><td>67 ± 13</td><td>71 ± 12</td><td>0.009</td><td>72 ± 14</td><td>72 ± 14</td><td>0.957</td></tr><tr><td>Killip class ≥ 2, n (%)</td><td>28 (18)</td><td>42 (45)</td><td>&lt;0.001</td><td>56 (29)</td><td>20 (42)</td><td>0.148</td></tr><tr><td>Peak creatine kinase (IU/L)</td><td>2259 [1160 to 3885]</td><td>2874 [464 to 4247]</td><td>0.100</td><td>2381 [1508 to 4211]</td><td>2078 [1291 to 3708]</td><td>0.554</td></tr><tr><td>End-systolic volume (mL/m²)</td><td>28 ± 10</td><td>32 ± 12</td><td>&lt;0.001</td><td>27 ± 11</td><td>27 ± 10</td><td>0.991</td></tr><tr><td>End Diastolic Volume (mL/m²)</td><td>56 ± 15</td><td>58 ± 15</td><td>0.333</td><td>52 ± 14</td><td>51 ± 15</td><td>0.663</td></tr><tr><td>Left ventricular ejection fraction (%)</td><td>51 ± 9</td><td>45 ± 9</td><td>&lt;0.001</td><td>50 ± 8</td><td>49 ± 8</td><td>0.519</td></tr><tr><td>Wall motion score index</td><td>1.84 ± 0.15</td><td>1.91 ± 0.15</td><td>&lt;0.001</td><td>1.90 ± 0.15</td><td>1.93 ± 0.15</td><td>0.250</td></tr><tr><td>Treatment at discharge, n (%)</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Antiplatelet therapy</td><td>157 (99)</td><td>90 (96)</td><td>0.277</td><td>188 (99)</td><td>47 (100)</td><td>1</td></tr><tr><td>β-blockers</td><td>152 (96)</td><td>86 (91)</td><td>0.288</td><td>187 (98)</td><td>43 (91)</td><td>0.042</td></tr><tr><td>ACE-I or ARB</td><td>153 (96)</td><td>91 (97)</td><td>1</td><td>186 (98)</td><td>44 (94)</td><td>0.285</td></tr><tr><td>Aldosterone antagonists</td><td>15 (9)</td><td>19 (20)</td><td>0.025</td><td>65 (34)</td><td>13 (28)</td><td>0.495</td></tr><tr><td>Statins</td><td>158 (99)</td><td>90 (96)</td><td>0.125</td><td>181 (95)</td><td>44 (94)</td><td>0.929</td></tr></tbody></table>

Results are presented as mean ± SD frequency (percentage) or median [Q1 to Q3]. The selected variables are significantly different in the two groups for both cohorts. Bold and bold italics results are significant (p<0.05) respectively in both cohorts and only one cohort, with p values calculated by appropriate statistical tests (chi2 or t-test). Imputed data are not indicated in this table.

The differential expression analysis performed between the identified groups showed that 43 and 41 of the 50 selected proteins had significantly different means between the two groups respectively in REVE-1 and REVE-2 (Table S4). These results highlight that the groups identified in REVE-1 have significantly distinct proteomic expression profiles that were validated in REVE-2.

## Event prediction

Cumulative incidence curves in REVE-1 (Figure 3A) and REVE-2 (Figure 3B) showed that in both cohorts, the identified groups had distinct incidence with higher hospitalization for HF in group 2. We validated the performance of group 2 for the significantly higher risk of hospitalization for HF with an SHR of 7.26 ([3.74 - 14.07], p<0.001) in REVE-1 and 3.66 ([1.75 - 7.64], p<0.001) in REVE-2. The group's robustness remained significant after adjustments for age, gender, ejection fraction, diabetes, Killip class, serum creatinine, BNP and NT-proBNP separately or combined in pairs. The least favorable model for the group variable was those adjusted for age and ejection fraction which showed a significant association with an SHR, for the group variable, of 4.50 ([2.19 - 9.26], p<0.001) and 2.19 ([1.05 - 4.93], p = 0.047) in REVE-1 and REVE-2, respectively (Table S5).
<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /><figcaption>**Figure 3. Cumulative incidence curves of HF for patients from REVE-1 (A) and REVE-2 (B) cohorts**</figcaption></figure>

The studied populations are divided in two groups according to the two identified groups by clustering with the k-means.
See also Table S3 and Figure S3.

A sensitivity analysis was achieved by performing the clustering on the 50 selected proteins without the BNP, an established biomarker of HF. When performing the clustering on 49 proteins, only 2 patients from REVE-1 changed of groups and 1 patient from REVE-2 changed of groups. The cumulative incidence curves of these groups were very similar to those obtained with the 50 proteins (Figure S3).

## Signaling pathway analysis

To gain further insight into the potential mechanisms in worse outcome of post-MI patients, the 50 proteins significantly associated with occurrence of HF (Table 2) were subjected to IPA for building biological networks. Two networks ranked by a score >20, revealed a significant link with "Cell-To-Cell Signaling and Interaction, Cellular Growth and Proliferation, Gene Expression" (network 1) and "Cancer, Hair and Skin Development and Function, Organ Development" (network 2). Despite their low score, networks 8 and 9 are linked to respectively "Cardiovascular Disease, Cellular Assembly and Organization, Developmental Disorder" and "Cardiovascular Disease, Hematological Disease, Hereditary Disorder" (Table 4, Figure 4).

The proteins among the 50 selected included in the networks have peripheral position except for PRKN (parkin) which is indirectly related to NPPB (natriuretic peptides), B2M (beta-2-microglobulin), and KAT2B (histone acetyltransferase KAT2B) (Figure 4A). The GO enrichment analysis selected two significant
<table><caption>Table 4. Significant protein networks identified by IPA in REVE1</caption><thead><tr><th>ID</th><th>Molecules in network</th><th>Score</th><th>Focus molecules</th><th>Top diseases and functions</th></tr></thead><tbody><tr><td>1</td><td>AMBRA1, B2M, BCL2, CARTPT, CD1D, CRIP2, DPP4, EFEMP1, ENY2, ERK1/2, ESR1, GATA4, HDAC1, Histone H3, HLA-G, HPSE, HTATIP2, KAT2B, KNG1, MAOB, MBD2, NCOA5, NPPB, Patched, PDX1, POU5F1, PPP1R1A, PRKN, RARA, RGS3, RNF41, RXRA, S100A3, SAMHD1, SERPINA7</td><td>30</td><td>15</td><td>Cell-To-Cell Signaling and Interaction, Cellular Growth and Proliferation, Gene Expression</td></tr><tr><td>2</td><td>ACTA2, ADIPOQ, CCND1, CDH6, CDK5, CTNNB1, CUL1, DISC1, FLNA, HERC5, HSPD1, IKBKE, ITCH, ITIH2, JPT1, KLF4, LGR4, LGR5, MAP2K4, MAP3K1, MAP3K7, NONO, NOTUM, PDE4A, PDIA3, PKM, PRKAA1, RGS3, SERPINA4, SMAD3, STAT3, THRB, TNFRSF19, VEGF, WNT3A</td><td>22</td><td>12</td><td>Cancer, Hair and Skin Development and Function, Organ Development</td></tr><tr><td>3</td><td>FGF2, FGFBP3</td><td>2</td><td>1</td><td>Carbohydrate Metabolism, Drug Metabolism, Small Molecule Biochemistry</td></tr><tr><td>4</td><td>NFAT5, PRSS35</td><td>2</td><td>1</td><td>Organismal Injury and Abnormalities, Renal and Urological Disease, Renal Hydronephrosis</td></tr><tr><td>5</td><td>MAFB, RNASE1</td><td>2</td><td>1</td><td>Cancer, Organismal Injury and Abnormalities, Respiratory Disease</td></tr><tr><td>6</td><td>CBLN4, DCC, NTN1</td><td>2</td><td>1</td><td>Cell Morphology, Cellular Assembly and Organization, Cellular Development</td></tr><tr><td>7</td><td>CLOCK, SCN4B, SCN5A</td><td>2</td><td>1</td><td>Developmental Disorder, Neurological Disease, Organismal Survival</td></tr><tr><td>8</td><td>APP, CST3, CTSB, CTSD, MMP25</td><td>2</td><td>1</td><td>Cardiovascular Disease, Cellular Assembly and Organization, Developmental Disorder</td></tr><tr><td>9</td><td>CUL2, F2, GZMA, HABP2, SDC1, SERPINC1, USP2</td><td>2</td><td>1</td><td>Cardiovascular Disease, Hematological Disease, Hereditary Disorder</td></tr></tbody></table>

The networks were identified using IPA computational algorithms consisting of the 50 selected proteins (Table 2) and their direct interactions with other proteins ("interconnecting") in the knowledge base. Scores were calculated for each network according to its fit to the set of selected proteins and used to rank networks on the Ingenuity analysis (version 73,620,684 Release: 2022-03-12). The scores take into account the number of selected proteins and the size of the network to approximate the relevance of the network to the original list of selected proteins. Selected molecules are in bold font (see Table 2) and interconnecting proteins identified in the network are in normal font. See also Tables S6 and S7.

clusters "regulation of protein binding" (GO ID: 0,043,393) and "regulation of peptidase activity" (GO ID: 0,052,547). The first cluster identified five proteins with two present in network 1 (PRKN, and B2M) and two present in network 2 (WNT3A and ADIPOQ) (Figure 4B). The second cluster identified seven proteins with three present in network 2 (WNT3A, SERPINA4 and ITIH2), two in network 1 (KNG1 and SERPINA7) and one (CST3) in network 8 and, one (SERPINC1) in network 9 (Table S6).

We also examined the significant relationship with the biological pathways involved in "Function and Diseases pathways" (ranked below 500 in IPA) related to "cardiovascular disease" and "cardiovascular system development and function" (Table S7). The key proteins were in the high-scoring networks as well in networks 8 and 9 with NPPB and DPP4 (dipeptidylpeptidase 4) (network 1), PDE4A (phosphodiesterase 4A) and ADIPOQ (adiponectin) (network 2), CST3 (Cystatin-C) (network 8) and SERPINC1 (antithrombin-III) (network 9) present in most cardiovascular diseases selected. All these proteins were associated with hypertension for which we found significant differences between the two groups of variables identified in both cohorts (Table 3).

DISCUSSION

This study investigated the ability of a large set of biomarkers to identify post-MI patients with long-term occurrence of hospitalization for HF. Our results showed that a subset of 50 plasma proteins measured at time of hospitalization allowed identifying two groups of patients with distinct proteomic expression
<figure><img src="image_8.png" /><figcaption>B</figcaption></figure>

<figure><img src="image_9.png" /></figure>
<figure><img src="image_10.png" /></figure>

**Figure 4. IPA based protein networks involved in worse outcome in post-MI patients of REVE-1 study**

The IPA analysis was performed on the 50 proteins selected to be associated with long term survival and listed in Table 2 (version: 73,620,684 Release: 2022-03-12). Only network 1 (A) and network 2 (B) from Table 4 are presented. Nodes are displayed using various shapes that represent the functional class of the proteins as published on https://qiagen.secure.force.com/KnowledgeBase/KnowledgeIPAPage?id=kA41i000000L5rTCAS. The color of proteins indicates respectively their regulated expression (red: increased and green: decreased) associated with occurrence of hospitalization for HF. The arrows indicate the modulatory effect of protein on its interacting proteins. Only direct interactions were selected. Detailed information on the molecules present in the networks is detailed in Table 4. See also Table S7.

profiles associated with different occurrence of HF hospitalization. Network analysis identified pathways up-regulated in MI patients with high risk of long-term adverse outcome related to cell-to-cell signaling and cardiovascular disease.

**Protein signature for high risk of adverse outcome after in patients with MI**

The relationships between individual biomarkers and adverse outcome following MI have been previously reported. Increased levels of BNP, cardiac troponin, and C-reactive protein were associated with major cardiac events after MI.<sup>8-10</sup> In more recent studies, matrix metalloproteinases and other biomarkers of extracellular matrix turnover have also been shown to predict outcome in this setting.<sup>11,12</sup> For decades, risk prediction in clinical practice has been based on generally available clinical characteristics and conventional cardiac biomarkers such as BNP.

The present study is one of the first to assess the relationship between a large set of biomarkers (>4,500 proteins) and long-term clinical outcome using high-throughput technology combined with state-of-art statistical and clustering analyses. We used two prospective cohorts, REVE-1 and REVE-2, which included patients suffering from a first acute MI with blood sampling during hospitalization. Patients underwent long-term follow-up with nearly no loss of patients. Thanks to the SOMAscan assay (version 4.0), 4,668 proteins were measured in the plasma of all the patients included in both cohorts. These proteins were not targeted toward any particular disease allowing us to study the link between the event and the proteins without preconceived idea of which proteins should be investigated thus leading to the potential discovery of new biomarkers.

The high-dimensional data generated by these two studies, with more variables than individuals, makes standard statistical analysis unusable. To minimize the risk of overfitting, state of the art statistical approaches with rigorous stability selection procedure were used to ensure reliability of our finding. First, we selected 50 proteins over the large panels of proteins that were significantly associated with the occurrence of HF. Second, a clustering algorithm was used on these 50 proteins giving the same weight to all of them, resulting in two groups where the information of all the proteins is equally represented. Third, competing risk models were then used to model the occurrence of hospitalization for HF against death for all causes using the group variable: low and high risk of occurrence of HF. To validate these effects of groups, we used the validation cohort, REVE-2, for which standardization parameters were set up from the data of the derivation cohort, REVE-1, enabling to use them for other cohorts as we did with REVE-2. We confirmed the effect of groups on occurrence of HF characterized in REVE-1 using competing risk models in REVE-2, as in REVE-1. These two groups of patients only identified by their circulating levels of proteins also showed differences in the clinical characteristics allowing the stratification of the patients in two groups, low and high risk of adverse outcome following MI. These results may help clinicians for more targeted and personalized treatments for patients regarding their cardiac outcome.

The two identified groups of this study based on the measurement of 50 selected proteins with the same technology can easily be applied to other cohorts (or single patient) using the standardization parameters set up.

**Translation of the identified proteins into biological pathways**

The proteins found in our network analysis were translated into biological pathways typically related to HF. The network analysis showed that pathways specifically upregulated in MI patients with high risk of hospitalization for HF (group 2) were related to cell-to-cell communication and cardiovascular disease.

The key proteins are NT-proBNP, ADIPOQ, SERPINC1, and WNT3A. NT-proBNP is associated with cardiac stretch with plasma levels widely used for screening and diagnosis of HF<sup>13</sup> and was previously found to be a specific hub in network analysis of patients with HF and reduced ejection fraction in two independent
studies.¹⁴,¹⁵ ADIPOQ has been shown to affect the autophagic response in the heart and contribute to accelerate cardiac remodeling.¹⁶ NT-proBNP was closely related in network 1 with PRKN, protein involved in autophagy and both were represented in most cardiovascular diseases (Table S7). SERPINC1 was also highly represented in cardiovascular diseases and by GO enrichment, but its potential as biomarker in HF has not been described up to now. Same for WNT3A, which is found in the 2 selected clusters by GO analysis and has been shown to be involved in cardiac muscle cell differentiation and its upregulation has been shown to be involved in TGFβ1-induced cardiac hypertrophy¹⁷ but also in cardiomyocyte injury following hypoxia.¹⁸

## Clinical implications and future perspectives

Individualized risk assessment is an integral part of management for patients in the post-MI setting. Early determination of the 50 selected proteins may help to detect patients at high risk of adverse outcome in the post-MI period. Such identification may encourage more aggressive therapy for this high-risk group. The identified groups were not impacted by the clinical differences between the two cohorts, showing the robustness of the two groups to clinical variations. Regarding this, the group's ability to predict HF for patients from other cohorts should remain relevant. Finally, the pathophysiological mechanisms beyond the selected proteins enriched in the networks that are highly represented in cardiovascular diseases remains to be established.

## Limitations of the study

The mean sampling time of the two cohorts are different but the results are still significant, which shows that the identified groups are robust to sampling time variations and could be used for clinical prognosis.

The choice of the silhouette criteria for the number of groups could be discussed as many criteria exist in the literature. We believe that two groups are appropriate and lessen the risk of overfitting. Although we believe our strategy valid, external validation in large cohorts from other areas/countries is mandatory to confirm the predictive value of the identified biomarkers subset to classify the MI patients in high- or low-risk of adverse outcome. Our clustering approach was validated in REVE-2, the validation cohort. The patients from the validation cohort were recruited in the same region as the derivation cohort, but standardization parameters were set up for an external validation with other cohorts of patients recruited everywhere in which circulating proteins were measured by a similar technology.

Although most patients had acute reperfusion, the proportion of patients with primary PCI reflects the practice in 2002–2004 and 2006–2008 and was lower than it would be nowadays. Finally, our study populations consisted of mainly men (74 and 81%, respectively in REVE-1 and REVE-2), therefore our prediction model may be less suitable for women. In addition, the patients recruited in the two studies suffered from severe MI and the results might not apply to the overall population of patients after MI.

## Conclusions

Here, we have investigated the cardiac prognostic implications of the highest panel of biomarkers available in MI patients, thanks to the aptamer-based platform. This study has several clinical implications. First, we improve the prediction of the long-term risk of HF occurrence in MI patients; second, the results obtained provide biological context for long-term adverse cardiac outcomes. The proteins found in our network analysis were in link with cardiovascular diseases.

We were able to validate our findings in an independent cohort, significantly reducing the overfitting effect by focusing on proteins linked to the outcome to ensure that results are specific to HF. These promising findings will require external validation in additional ethnic groups of patients.

## STAR★METHODS

Detailed methods are provided in the online version of this paper and include the following:

* KEY RESOURCES TABLE

* RESOURCE AVAILABILITY

* Lead contact

* Materials availability
* Data and code availability

* **METHOD DETAILS**

* Data deposition and materials sharing

* Study populations

* **QUANTIFICATION AND STATISTICAL ANALYSIS**

* Standardization of proteomic measurement

* Imputation of missing clinical values

* Proteins selection

* Clustering of patients based on selected proteins

* Groups' effect on the occurrence of HF

* Enrichment analysis

* Ingenuity Pathways Analysis

* SOMAscan V4 data standardization and file specification technical note

**SUPPLEMENTAL INFORMATION**

Supplemental information can be found online at https://doi.org/10.1016/j.isci.2023.106171.

**ACKNOWLEDGMENTS**

This work was supported by grants from "Fédération Française de Cardiologie". W.H. was supported by grants from "Fondation pour la Recherche Médicale" with the reference PBR201903008963. FP has received a grant by the French Government, managed by the National Research Agency(ANR) under the program "Investissements d'avenir" with the reference ANR-16-RHUS-0003. CB and FP are members of FHU-CARNAVAL. We acknowledge Somalogic Inc. as the provider of the proteomics analysis.

None of the funders had any role in study design, data collection, data analysis, interpretation, and writing the article.

**AUTHOR CONTRIBUTIONS**

C.B., P.A., and F.P. contributed equally in the conception of the study. All authors contributed equally in the design of the study and in the analysis and interpretation of data. W.H., V.V., G.M., C.B., and F.P. drafted the manuscript and the final version. All authors read, critically revised and approved the final version of the manuscript.

**DECLARATION OF INTERESTS**

The authors disclose no conflict of interest.

**INCLUSION AND DIVERSITY**

We support inclusive, diverse, and equitable conduct of research.

Received: July 22, 2022

Revised: December 19, 2022

Accepted: February 6, 2023

Published: February 11, 2023

**REFERENCES**

1. Benjamin, E.J., Blaha, M.J., Chiuve, S.E., Cushman, M., Das, S.R., Deo, R., de Ferranti, S.D., Floyd, J., Fornage, M., Gillespie, C., et al. (2017). Heart disease and stroke statistics-2017 update: a report from the American heart association. Circulation 135, e146–e603. https://doi.org/10.1161/CIR.0000000000000485.Heart.

2. Eggers, K.M., Lindhagen, L., Lindhagen, L., Baron, T., Erlinge, D., Hjort, M., Jernberg, T., Marko-Varga, G., Rezeli, M., Spaak, J., and Lindahl, B. (2021). Predicting outcome in

acute myocardial infarction: an analysis investigating 175 circulating biomarkers. Eur. Heart J. Acute Cardiovasc. Care 10, 806–812. https://doi.org/10.1093/ehjacc/zuaa014.

3. Liu, R.X., Thiessen-Philbrook, H.R., Vasan, R.S., Coresh, J., Ganz, P., Bonventre, J.V., Kimmel, P.L., and Parikh, C.R. (2021). Comparison of proteomic methods in evaluating biomarker-AKI associations in cardiac surgery patients. Transl. Res. 238, 49–62. https://doi.org/10.1016/j.trsl.2021.07.005.

4. Gui, H., She, R., Luzum, J., Li, J., Bryson, T.D., Pinto, Y., Sabbah, H.N., Williams, L.K., and Lanfear, D.E. (2021). Plasma proteomic profile predicts survival in heart failure with reduced ejection fraction. Circ. Genom. Precis. Med. 14, e003140. https://doi.org/10.1161/CIRGEN.120.003140.

5. Savoye, C., Equine, O., Tricot, O., Nugue, O., Segrestin, B., Sautière, K., Elkohen, M., Pretorian, E.M., Taghipour, K., Philias, A., et al. (2006). Left ventricular remodeling after anterior wall acute myocardial infarction in
modern clinical practice (from the REmodelage VEntriculaire [REVE] study group). Am. J. Cardiol. 98, 1144–1149.
https://doi.org/10.1016/j.amjcard.2006.06.011.

6. Fertin, M., Hennache, B., Hamon, M., Ennezat, P.V., Biausque, F., Elkohen, M., Nugue, O., Tricot, O., Lamblin, N., Pinet, F., and Bauters, C. (2010). Usefulness of serial assessment of B-type natriuretic peptide, troponin I, and C-reactive protein to predict left ventricular remodeling after acute myocardial infarction (from the REVE-2 study). Am. J. Cardiol. 106, 1410–1416.
https://doi.org/10.1016/j.amjcard.2010.06.071.

7. Bauters, C., Dubois, E., Porouchani, S., Saloux, E., Fertin, M., de Groote, P., Lamblin, N., and Pinet, F. (2017). Long-term prognostic impact of left ventricular remodeling after a first myocardial infarction in modern clinical practice. PLoS One 12, e0188884. https://doi.org/10.1371/journal.pone.0188884.

8. de Lemos, J.A., McGuire, D.K., and Drazner, M.H. (2003). B-type natriuretic peptide in cardiovascular disease. Lancet 362, 316–322.
https://doi.org/10.1016/S0140-6736(03)13976-1.

9. Ohlmann, P., Jaquemin, L., Morel, O., El Behlgiti, R., Faure, A., Michotey, M.O., Beranger, N., Roul, G., Schneider, F., Bareiss, P., and Monassier, J.P. (2006). Prognostic value of C-reactive protein and cardiac troponin I in primary percutaneous interventions for ST-elevation myocardial infarction. Am. Heart J. 152, 1161–1167.
https://doi.org/10.1016/j.ahj.2006.07.016.

10. Suleiman, M., Khatib, R., Agmon, Y., Mahamid, R., Boulos, M., Kapeliovich, M., Levy, Y., Beyar, R., Markiewicz, W., Hammerman, H., and Aronson, D. (2006). Early inflammation and risk of long-term development of heart failure and mortality in survivors of acute myocardial infarction predictive role of C-reactive protein. J. Am. Coll. Cardiol. 47, 962–968. https://doi.org/10.1016/j.jacc.2005.10.055.

11. Fertin, M., Lemesle, G., Turkieh, A., Beseme, O., Chwastyniak, M., Amouyel, P., Bauters, C., and Pinet, F. (2013). Serum MMP-8: a novel indicator of left ventricular remodeling and cardiac outcome in patients after acute myocardial infarction. PLoS One 8, e71280.
https://doi.org/10.1371/journal.pone.0071280.

12. Eschalier, R., Fertin, M., Fay, R., Bauters, C., Zannad, F., Pinet, F., and Rossignol, P. (2013).

Extracellular matrix turnover biomarkers predict long-term left ventricular remodeling after myocardial infarction insights from the reve-2 study. Circ. Heart Fail. 6, 1199–1205.
https://doi.org/10.1161/CIRCHEARTFAILURE.113.000403.

13. Liang, F., O'Rear, J., Schellenberger, U., Tai, L., Lasecki, M., Schreiner, G.F., Apple, F.S., Maisel, A.S., Pollitt, N.S., and Protter, A.A. (2007). Evidence for functional heterogeneity of circulating B-type natriuretic peptide. J. Am. Coll. Cardiol. 49, 1071–1078. https://doi.org/10.1016/j.jacc.2006.10.063.

14. Tromp, J., Khan, M.A.F., Mentz, R.J., O'Connor, C.M., Metra, M., Dittrich, H.C., Ponikowski, P., Teerlink, J.R., Cotter, G., Davison, B., et al. (2017). Biomarker profiles of acute heart failure patients with a mid-range ejection fraction. JACC. Heart Fail. 5, 507–517. https://doi.org/10.1016/j.jchf.2017.04.007.

15. Tromp, J., Khan, M.A.F., Klip, I.T., Meyer, S., de Boer, R.A., Jaarsma, T., Hillege, H., van Veldhuizen, D.J., van der Meer, P., and Voors, A.A. (2017). Biomarker profiles in heart failure patients with preserved and reduced ejection fraction. J. Am. Heart Assoc. 6, e003989.
https://doi.org/10.1161/JAHA.116.003989.

16. Kamareddine, L., Ghantous, C.M., Allouch, S., Al-Ashmar, S.A., Anlar, G., Kannan, S., Djouhri, L., Korashy, H.M., Agouni, A., and Zeidan, A. (2021). Between inflammation and autophagy: the role of leptin-adiponectin Axis in cardiac remodeling. J. Inflamm. Res. 14, 5349–5365. https://doi.org/10.2147/JIR.S322231.

17. Li, T., Weng, X., Cheng, S., Wang, D., Cheng, G., Gao, H., and Li, Y. (2021). Wnt3a upregulation is involved in TGFβ1-induced cardiac hypertrophy. Cytokine 138, 155376.
https://doi.org/10.1016/j.cyto.2020.155376.

18. Li, Y., Du, L., Cheng, S., Guo, J., Zhu, S., Wang, Y., and Gao, H. (2020). Hypoxia exacerbates cardiomyocyte injury via upregulation of Wnt3a and inhibition of Sirt3. Cytokine 136, 155237. https://doi.org/10.1016/j.cyto.2020.155237.

19. Moaddel, R., Ubaida-Mohien, C., Tanaka, T., Lyashkov, A., Basisty, N., Schilling, B., Semba, R.D., Franceschi, C., Gorospe, M., and Ferrucci, L. (2021). Proteomics in aging research: a roadmap to clinical, translational research. Aging Cell 20, e13325. https://doi.org/10.1111/acel.13325.

20. Faquih, T., Mook-Kanamori, D.O., Rosendaal, F.R., Baglin, T., Willems van Dijk, K., and van

Hylckama Vlieg, A. (2021). Agreement of aptamer proteomics with standard methods for measuring venous thrombosis biomarkers. Res. Pract. Thromb. Haemost. 5, e12526. https://doi.org/10.1002/rth2.12526.

21. Luo, Y., Wadhawan, S., Greenfield, A., Decato, B.E., Oseini, A.M., Collen, R., Shevell, D.E., Thompson, J., Jarai, G., Charles, E.D., and Sanyal, A.J. (2021). SOMAscan proteomics identifies serum biomarkers associated with liver fibrosis in patients with NASH. Hepatol. Commun. 5, 760–773.
https://doi.org/10.1002/hep4.1670.

22. Cuvelliez, M., Vandewalle, V., Brunin, M., Beseme, O., Hulot, A., de Groote, P., Amouyel, P., Bauters, C., Marot, G., and Pinet, F. (2019). Circulating proteomic signature of early death in heart failure patients with reduced ejection fraction. Sci. Rep. 9, 19202. https://doi.org/10.1038/s41598-019-55727-1.

23. Josse, J., and Husson, F. (2016). MissMDA: a package for handling missing values in multivariate data analysis. J. Stat. Softw. 70, 1–31. https://doi.org/10.18637/jss.v070.i01.

24. Gray, R.J. (1988). A class of K-sample tests for comparing the cumulative incidence of a competing risk. Ann. Stat. 16, 1141–1154.
https://doi.org/10.1214/aos/1176350951.

25. Fine, J.P., and Gray, R.J. (1999). A proportional hazards model for the subdistribution of a competing risk. J. Am. Stat. Assoc. 94, 496–509. https://doi.org/10.1080/01621459.1999.10474144.

26. Bonferroni, C.E. (1936). Teoria statistica delle classi e calcolo delle probabilità. Pubbl. del R Ist. Super. di Sci. Econ. e Commer. di Firenze, 3–62. https://doi.org/10.1007/978-1-4419-9863-7_1213.

27. MacQueen, J. (1967). Some methods for classification and analysis of multivariate observations. 5th Berkeley Symp. Math. Stat. Probab. 5.1, 281–297.

28. Rousseeuw, P.J. (1987). Silhouettes: a graphical aid to the interpretation and validation of cluster analysis. J. Comput. Appl. Math. 20, 53–65. https://doi.org/10.1016/0377-0427(87)90125-7.

29. Wu, T., Hu, E., Xu, S., Chen, M., Guo, P., Dai, Z., Feng, T., Zhou, L., Tang, W., Zhan, L., et al. (2021). ClusterProfiler 4.0: a universal enrichment tool for interpreting omics data. Innovation 2, 100141. https://doi.org/10.1016/j.xinn.2021.100141.
## STAR★METHODS

### KEY RESOURCES TABLE

<table><thead><tr><th>REAGENT or RESOURCE</th><th>SOURCE</th><th>IDENTIFIER</th></tr></thead><tbody><tr><td><strong>Biological samples</strong></td><td></td><td></td></tr><tr><td>Plasma samples from REVE-1 study<br>(n=254 patients)</td><td>REVE-1 study<br>CHU deLille, FRANCE</td><td>https://doi.org/10.1016/j.amjcard.2006.06.011. Epub 2006 Aug 31.</td></tr><tr><td>Plasma samples from REVE-2 study<br>(n=238 patients)</td><td>REVE-2 study<br>CHU deLille, FRANCE</td><td>https://doi.org/10.1016/j.amjcard.2010.06.071. Epub 2010 Oct 1.</td></tr><tr><td><strong>Critical commercial assays</strong></td><td></td><td></td></tr><tr><td>BNP measurement by 2-site sandwich assay</td><td>Siemens Diagnostic, Zurich, Switzerland</td><td>Advia Centaur XPT</td></tr><tr><td><strong>Deposited data</strong></td><td></td><td></td></tr><tr><td>Somascan technology Version 4.0</td><td>This paper, note 1</td><td>See Tables 2 and S4</td></tr><tr><td><strong>Software and algorithms</strong></td><td></td><td></td></tr><tr><td>Ingenuity Pathways analysis</td><td>Ingenuity Systems</td><td>Version: 73620684, Release: 2022-03-12</td></tr><tr><td>R</td><td>This paper</td><td>Version 4.0.2</td></tr></tbody></table>

### RESOURCE AVAILABILITY

#### Lead contact

Further information and requests for resources and reagents should be directed to and will be fulfilled by the lead contact, Florence Pinet (florence.pinet@pasteur-lille.fr).

#### Materials availability

This study did not generate new unique reagents.

There are restrictions to the availability of SOMAScan (property of Somalogic company) and human samples (French cohorts of patients associated with clinical information).

#### Data and code availability

Any additional information required to reanalyze the data reported in this paper is available from the lead contact upon request.

Data: All (clinical and proteomic) data reported in this paper will be shared by the lead contact upon request.

Code: This paper does not report original code.

### METHOD DETAILS

#### Data deposition and materials sharing

The data and methods used in the analysis are available to any researcher for purposes of reproducing the results or replicating the procedures.

#### Study populations

The REVE studies have been previously reported.⁵,⁶ REVE-1 (inclusion period, February 2002-June 2004, n=266 patients) was designed to test the hypothesis that genetic polymorphisms in candidate genes may be associated with left ventricular remodeling.⁷ REVE-2 (inclusion period, February 2006-September 2008, n=246 patients) was designed to analyze the association of circulating biomarkers with left ventricular remodeling.⁶ Both studies were prospective with a multicentric recruitment. The research protocol of both studies was approved by the ethics committee of the "Centre Hospitalier et Universitaire de Lille" (Lille, France), and written informed consent was obtained from each patient. Requalification of samples for
SOMAscan analysis was obtained by "Comité de Protection des personnes Nord Ouest IV" (February 2018). IRB approval was obtained and that subjects gave informed consent. The inclusion criteria were the same: a first anterior Q-wave MI with ≥ 3 akinetic segments at predischarge echocardiography. Exclusion criteria were inadequate echocardiography image quality, life-limiting noncardiac disease, significant valvular disease, or previous Q-wave MI in both studies. A long-term clinical follow-up was performed by contacting the general practitioner or cardiologist, or the patients themselves.⁷ We collected data on death and hospitalization for HF. All events occurring during follow-up were adjudicated by two investigators with a third opinion in cases of disagreement. For hospitalizations during the follow-up period, hospital records were reviewed for evidence of clinical events. The events reported by the patients were systematically confirmed from the medical records. Hospitalization for HF was defined as hospitalization for symptoms of dyspnea or edema, elevated venous pressure, or interstitial or alveolar edema on chest X-ray, or the addition of intra-venous diuretics or inotropic medications. The primary endpoint of the present study was hospitalization for HF and the competing event was death from all causes.

### Plasma proteomics measurements

Plasma protein levels were measured in both cohorts with SOMAscan technology (version V4.0).¹⁹⁻²¹ Peripheral blood samples were collected in EDTA-treated tubes for 255 patients in REVE-1 and 238 patients in REVE-2 after MI during the initial hospitalization. The mean blood sampling day was 7.4 ± 3.9 days and 4.1 ± 2.2 after MI, respectively in REVE-1 and REVE-2. Blood samples were then assayed according to the manufacturer's protocol, as previously described with the 1·3k assay.²² Protein levels from the SOMAScan assay are expressed as relative fluorescence unit (RFU).

## QUANTIFICATION AND STATISTICAL ANALYSIS

### Standardization of proteomic measurement

There were no missing values in the set of the 5284 proteins quantified in both cohorts. A total of 414 SOMAmer measurements were removed from the analysis on request of the manufacturer. In the data set, 197 identical proteins were measured using two to three SOMAmers (Table S1). In order to give the same weight to each protein in the data set, for the proteins measured by several SOMAmers, the mean was calculated leading to 4668 proteins ready for analysis.

Proteomic expression levels were log2 transformed. Data quality was checked, and no variation was found between the 15 plates used for the measurements (Figure S1). Results of brain natriuretic peptide (BNP) measurements with the SOMAScan assay were also compared to BNP measurements performed with the automated 2-site sandwich immunoassay on an Advia Centaur (Siemens Diagnostic, Zurich, Switzerland) and showed high correlation (r=0.992). Log2 transformed protein expression levels were standardized (centered and reduced) for each protein on REVE-1, the derivation cohort. These standardization parameters (mean and standard deviation) used in REVE-1 (Table S2) were then applied in REVE-2, the validation cohort.

### Imputation of missing clinical values

Over the 27 clinical variables studied corresponding to 6858 and 6426 data of variables, respectively for all the patients from REVE-1 and REVE-2, 53 and 69 individual clinical data were missing respectively in REVE-1 (corresponding to 0.0077%) and REVE-2 (corresponding to 0.0107%). These values were imputed using the single imputation method missMDA.²³

Figure 1 shows the strategy used for analyzing the data obtained in the derivation (REVE-1) and validation (REVE-2) cohorts.

### Proteins selection

In order to only focus on the relevant proteins, with a significant effect on the occurrence of hospitalization for HF, univariate competing risk models were fitted for each protein individually using R package cmprsk (version 2.2-11).²⁴ For each protein, the occurrence of hospitalization for HF was modeled in competition with the occurrence of all-causes death using the competing risk model as previously defined.²⁵ Tests were performed on the subhazard ratio (SHR) to measure the significance of the relationship between the protein's expression and the occurrence of hospitalization for HF. Then, in order to take into account multiple testing, the significance threshold was set to 1.07·10⁻⁵, which corresponds to the FamilyWise Error Rate
(FWER) at 0.05 with the Bonferroni method²⁶ for 4668 tests. The proteins significantly associated with the event in REVE-1 were selected for the following analyses.

## Clustering of patients based on selected proteins

A k-means clustering algorithm²⁷ was applied in the REVE-1 study using the selected proteins as variables in order to define groups of patients. This algorithm builds groups of individuals sharing similar values for protein expressions regardless of the risk of hospitalization for HF of the patients. Individuals are gathered into groups characterized by their centers and each patient is assigned to the group whose center is the closest in terms of Euclidean distance. Therefore, a new categorical variable called "group variable" was defined with these affectations for the following analyses.

In order to choose a suitable number of groups, the overall average silhouette width was computed.²⁸ The silhouette refers to a method of interpretation and validation of consistency within clusters of data. The average silhouette width is calculated with Euclidean distance to measure how similar is the mean distance of each individual compared to other individuals in the same group and to the individuals from other groups. The silhouette ranges from -1 to +1, where the closest value to 1 indicates that the clustering is appropriate. This criterion was computed for a number of groups varying from 2 to 6. Finally, the number of groups which gave the higher overall average silhouette width was used.

As data were standardized using the same parameters for both cohorts, patients of REVE-2 were then assigned to the previously built groups in REVE-1 using the same method.

Group differences were identified using Welch's test of mean equality for quantitative variables (both for clinical data and for differential proteomic analysis) or chi squared test for clinical categorical variables. Differences corresponding to raw p-values below 0.05 were considered significant.

## Groups' effect on the occurrence of HF

In order to assess the interest of the previously created groups, cumulative incidence curves of hospitalization for HF were drawn for both groups. For each cohort, competing risk models were then used to model the occurrence of hospitalization for HF using the group variable. This allows to measure the strength and to test the significance of the group effect. Models were fitted for both cohorts in order to measure the effect of the groups in the derivation and the validation cohort. The models were then adjusted on the clinical variables age, gender, ejection fraction, diabetes, killip class, serum creatinine, BNP and NT-proBNP separately, in order to ensure the robustness of the information provided by the identified groups after adding relevant clinical information.

Sensitivity analysis was performed on the selected proteins without the established cardiac biomarkers (BNP, NT-proBNP). Incidence curve were drawn for the group identified with the selected proteins minus the established biomarkers and SHR competing risk model was fitted.

All statistical analyses were made using R (version 4.0.2).

## Enrichment analysis

Gene ontology enrichment analysis was performed on the set of selected proteins using Biological Processes, Molecular Function and Cellular Component subsets. Enrichments tests were performed using the R package clusterProfiler.²⁹ For each GO term, a Fisher exact test was performed in order to test for the overrepresentation of the set of selected proteins in the GO term.

## Ingenuity Pathways Analysis

The proteins of interest, selected to be significantly associated with occurrence of HF, were subjected to Ingenuity Pathways Analysis (IPA) (Current version: 73620684 Release: 2022-03-12, Ingenuity Systems) and used as a starting point for building biological networks. This analysis uses computational algorithms to identify networks consisting of focus proteins (proteins significantly modulated) and their interactions with other proteins ("non-focused") in the knowledge base. Scores were calculated for each network according to the fit of the network to the set of focus proteins and used to rank networks on the IPA database restricted to direct interactions. IPA uses the proteins from the highest-scoring network to extract a
connectivity pathway that relates candidate proteins to each other based on their direct interactions. IPA determines the "Function & Disease and Canonical Pathways" to be significantly associated with these candidate proteins.

# SOMAscan V4 data standardization and file specification technical note

## Overview

Normalization and calibration are routine numerical procedures developed to remove systematic biases in the raw assay data. Normalization is a sample-by-sample adjustment in overall signal within a single plate (run) performed across three non-consecutive steps: Hybridization Control Normalization, Intraplate Median Signal Normalization, and Median Signal Normalization to a Reference. Plate Scaling and Calibration is a SOMAmer® binding reagent-by-SOMAmer binding reagent adjustment that minimizes between-plate variability. Global reference standards are established for procedures with controls on each plate individual, QC, and Calibrator samples are normalized and calibrated to the established global reference standards. Separate calibrator global reference standards are established for each matrix (serum, plasma), and assay shifts or skew from the global reference standards are tracked over time. New global reference standards may be developed in concordance with changes in assay processes, performance, or reagents.

Hybridization Control Normalization was developed to remove systematic biases present in the raw data after slide feature aggregation from a slide-based hybridization microarray for assay readout and quantification. Hybridization Control Normalization is performed using a set of twelve hybridization control sequences measured independently for each sample array. The procedure is intended to correct for systematic effects on the data introduced during the hybridization readout and results in a single scale factor for each sample that is subsequently applied to the measured signal on all features within a subarray (sample).

Intraplate Median Signal Normalization uses all the SOMAmer reagent signals on a given subarray to remove sample or assay biases that may be due to differences between samples in overall protein concentration, pipetting variation, variation in reagent concentrations, assay timing, and any other source of systematic variability within a single plate. Each SOMAmer reagent is assigned to one of three dilution sets, scale factors are derived within dilution sets separately, and all SOMAmer binding reagents within each set are scaled together. Three sample dilutions will result in three independent median signal scale factors for each subarray (sample) in addition to the hybridization scale factor. This step is only applied to calibrator samples.

Plate Scaling and Calibration is accomplished using a number of replicate measurements of a common pooled calibrator sample consistent with the assay sample type for a study. Calibrator samples must be composed of identical sample matrices as the samples that are being calibrated. No protein spikes are added to the calibrator samples – SomaLogic relies solely on the endogenous levels of each analyte within a calibrator sample. Since calibration attempts to correct plate-to-plate variation and such variation can be idiosyncratic for SOMAmer binding reagents, a unique calibration scale factor is derived for each SOMAmer binding reagent within the assay. The median of these scale factors is then computed and applied across all SOMAmer measurements in that plate to account for the total signal difference (plate scale), and the scale factors are subsequently recalculated for each SOMAmer and applied to all measurements within the set of samples in that plate.

Median Signal Normalization to a Reference occurs on a per-sample basis, wherein a scale factor for a set of SOMAmer reagents is computed against a reference value generated from a cohort of healthynormal individuals and then aggregated within a dilution. The median of each dilution's scale factors is then applied to their respective SOMAmer reagents. This step is applied to QC, Buffer, and individual samples.

File specification: SOMAscan results are produced in a tab-delimited ASCII file with an ADAT extension (ADAT file). The ADAT file contains measurements for a series of analytes (columns) across a series of samples (rows) and includes analyte description and sample description information. The format is designed to provide flexibility for the number of samples as well as the number and types of analyte and sample descriptors.

## Normalization

**Hybridization Control Normalization.** A set of hybridization control sequences is added to each sample as part of the elution buffer in the SOMAscan assay. These hybridization controls are added at
<figure><img src="image_11.png" /></figure>

concentrations to give measured relative fluorescence units (RFU) that span the dynamic range of the assay. The global reference RFU value for each hybridization control is defined by the median signal measured within the current plate being normalized. A ratio is computed by dividing the median RFU for each control by its measured RFU in the sample. The median of these hybridization control measurement ratios in each subarray defines the sample-based hybridization scale factor. By definition, such a scaling will equate the median RFU for the hybridization controls to the median reference RFU for the controls. All SOMAmer reagent results within a sample are multiplied by the resulting hybridization scale factor increasing or decreasing the overall “brightness” of the sample. The procedure is displayed graphically in Note Figure S4.

**Intraplate Median Signal Normalization.** Intraplate Median Signal Normalization is performed on each sample dilution independently. In most matrices, each SOMAmer binding reagent is assigned to one of three dilution sets, scale factors are derived within dilution sets separately, and all SOMAmer reagents within each set are scaled together. Within each sample matrix, this is only performed on calibrator samples. Like Hybridization Control Normalization which uses a local reference standard, the local median reference RFU for each SOMAmer reagent is the median RFU for that SOMAmer binding reagent within the sample group (calibrator in buffer) in the plate to be normalized. As in hybridization normalization, a ratio is computed for each SOMAmer reagent by dividing the reference SOMAmer RFU by its measured RFU in the sample to be normalized. The median of the SOMAmer measurement ratios for all SOMAmer reagents in a dilution defines the sample-based scale factor for all SOMAmer reagents within that dilution and sample. All SOMAmer reagents within the dilution for a sample are scaled by the resulting median signal scale factor. Three sample dilutions will result in three independent median signal scale factors for each sample in addition to the hybridization scale factor as shown in Note Figure S5.

**Plate Scaling & calibration**

Clinical sample studies are plate scaled and calibrated to remove systematic assay variability. A set of control calibrator samples is used to detect and remove systematic variability between independent assay plates. Calibrator samples must be of the same type as the samples that are being calibrated. Calibrator global reference RFU values for each SOMAmer reagent are defined by the median signal measured on a set of samples spanning a number of independent assay plates that have been shown to meet assay acceptance criteria. For each SOMAmer reagent, the median RFU signal for that SOMAmer reagent across all the calibrator samples within the clinical study defines the global calibrator reference for that SOMAmer binding reagent. Note Figure S6 below displays the data from a typical clinical study and illustrates the systematic bias removed by calibration.

Plate scaling is performed on an entire independent plate. A local median reference value is derived for each SOMAmer reagent by computing the median RFU for that SOMAmer reagent from the set of replicate calibrator samples within the plate. The SOMAmer-based calibration scale factor is then computed by dividing the calibrator global reference RFU by the local median reference value defined for each SOMAmer reagent. The median of all scale factors for a given plate is then applied across all SOMAmer measures in the plate (plate) forcing the overall calibrator median signal to match the overall median signal within the global calibrator reference.

Plate-to-plate calibration is performed on each SOMAmer measurement within the plate independently. A local median reference value is derived for each SOMAmer reagent by computing the median RFU for that SOMAmer reagent from the set of replicate calibrator samples within the plate. The SOMAmer-based calibration scale factor is then computed by dividing the calibrator global reference RFU by the local median reference value defined for each SOMAmer reagent. This scale factor is applied to all SOMAmer measurements in the plate, forcing the median calibrator signal to match the global calibrator reference for that SOMAmer binding reagent. Each plate within a study has a unique calibration scale factor for each SOMAmer reagent. The data from Note Figure S5 are displayed after calibration in Note Figure S6.

**Median Normalization to a reference**

All individual, QC, and Buffer samples are then median normalized to a reference value. Unlike Intraplate Median Signal Normalization, Median Normalization to a Reference can be performed on a single sample due to the presence of an external global reference value generated from a cohort of healthynormal individuals for each SOMAmer reagent. This method is very similar to Intraplate Median Signal Normalization
in practice, the primary difference being the origination of the reference value. A ratio is computed for each SOMAmer reagent by dividing the global reference SOMAmer RFU by its measured RFU in the sample to be normalized. The median of the SOMAmer measurement ratios for all SOMAmer reagents in a dilution defines the sample-based scale factor for all SOMAmer reagents within that dilution and sample. All SOMAmer reagents within the dilution for a sample are scaled by the resulting median signal scale factor. Three sample dilutions will result in three independent median signal scale factors for each sample in addition to the hybridization scale factor.

### Acceptance criteria

Hybridization Control and Intraplate Median Signal Normalization scale factors are expected to be in the range of 0.4–2.5. The plate scale factor is expected to be between 0.4 and 2.5. The distribution of QC sample ratios is expected to have 85% of individual SOMAmer reagents in the total array between 0.84 and 1.19 (i.e. less than 15% in the tails of the distribution). Gaussian distributions of scale factors are expected. A report is provided for each study (single plate or set of plates) with the results of the Normalization and Calibration process.