Single-cell dissection of the immune response after acute myocardial infarction =============================================================================== * Irene V. van Blokland * Roy Oelen * Hilde E. Groot * Jan Walter Benjamins * Kami Pekayvaz * Corinna Losert * Viktoria Knottenberg * Matthias Heinig * Leo Nicolai * Konstantin Stark * Pim van der Harst * Lude H. Franke * Monique G. P. van der Wijst ## Abstract **Background** The role of the immune system in the context of acute myocardial infarction (MI), and its response to such event are poorly characterized, but is thought to be an important driver of myocardial remodeling, thromboinflammation, and exacerbation of atherosclerosis - triggering recurrent cardiovascular events. So far, anti-inflammatory approaches drugs have shown promising effects on the prevention of recurrent cardiovascular events or myocardial salvage after myocardial infarction. However, they broadly impair the immune system and some are associated with increased infectious side effects. Therefore, a more detailed understanding of the immune response to myocardial infarction is needed to tailor anti-inflammatory therapeutic approaches in MI patients. **Methods** To gain such detailed longitudinal understanding of the immune response in ST-elevated myocardial infarction (STEMI) patients, we compared peripheral blood mononuclear cell (PBMCs) single-cell RNA-sequencing (scRNA-seq) expression and plasma protein profiles over time and in comparison to age- and sex-balanced controls in 38 STEMI patients at hospital admission, 24 hours (acute phase) and 6-8 weeks (chronic phase) after STEMI. **Results** In total, 95,995 diseased and 33,878 control PBMCs were analyzed. Compared to controls, we observed a relative increase in the number of classical monocytes and a decrease in the number of CD56dim natural killer cells in STEMI patients at admission, and these differences persisted until 24 hours after STEMI. The monocytes also showed the largest gene expression changes in STEMI patients compared to controls, and in STEMI patients over time. These were associated with changes in the activity of toll-like receptors, IFN and IL-1 signaling. Subsequent differential cell-cell communication analysis suggested that these monocytes are mainly involved in the outgoing differential communication in the first 24h after a STEMI, whereas in the next 6-8 weeks they become mostly involved in the incoming differential communication. Finally, a targeted protein cardiovascular biomarker panel revealed 33 out of 92 plasma proteins to be changed during the acute and/or chronic phase after STEMI. Interestingly, the plasma levels of three of these proteins were found to be affected by genetic variation, disease status and time after STEMI. Indicating the importance of taking all these aspects into consideration when defining potential future therapies. **Conclusions** Altogether, our analyses have revealed the immunological pathways that are disturbed upon MI, and in which cell type and during which stage of the disease. Additionally, we also provide insights in which patients are expected to benefit most from anti-inflammatory treatments, by identifying the genetic variants and disease stage at which these variants affect the outcome of these (drug-targeted) pathways. These findings advance our knowledge of the immune response after MI and provide further guidance for future therapeutic studies. Keywords * Single-cell RNA sequencing * cardiovascular disease * acute myocardial infarction * inflammation * genetics ## Introduction Acute myocardial infarction (MI) is one of the leading causes of mortality and morbidity globally. The immune system plays an important role in the pathophysiology of the disease from the initial formation and rupture of atherosclerotic plaques, all the way up to inflammation in response to the MI (1–3). Systemically studying and targeting the inflammatory response has remained challenging due to its complexity and severe side effects, as was seen in the CANTOS trial (interleukin-1ß, IL-1ß) and ASSAIL-MI (interleukin 6 receptor, IL6R) trial (4,5). Identifying specific immune cell subpopulations and underlying processes that are specifically disturbed after MI could open the way for more targeted treatment that is expected to cause fewer side effects (6,7). However, until recently, analyzing the characteristics of specific cell types was often done at bulk level and using predefined marker genes, not providing an unbiased, comprehensive overview of the immune system (8). Recent developments in the field of single-cell RNA sequencing (scRNA-seq) now enable us to unbiasedly study 100,000s of individual cells simultaneously on a transcriptome-wide level (9–11). This technology enables us to gain new insights into the inflammatory response after a MI by flexibly and concurrently mapping MI-induced changes in cell type composition, gene expression level and downstream pathways at various cellular resolutions. To study MI at the single-cell level, it is essential to first understand the normal physiological variation among the general population. Such groundwork has been generated both for the heart (12) and the circulating immune cells (10). Various aspects in the pathophysiology of myocardial infarction have now been dissected at the single-cell level, from atherosclerotic aortas precluding the disease (13,14) to cardiac neovascularization by resident endothelial cells after the infarction (15). More recently, single-cell multi-omic maps (combining spatial transcriptomics, single-cell expression and chromatin accessibility profiling) have been generated from human heart tissue after MI (16). This revealed an increased spatial dependency between lymphoid and myeloid cells in ischaemic samples compared to healthy samples, indicating the importance of cellular communication between immune cells during cardiac repair after MI. Altogether, these previous studies have mainly focused on the heart, whereas each of these studies also indicate the importance of the immune cells during a MI. However, we still lack a more detailed single-cell level description of how the circulating immune cells are affected by a MI, both during the acute and chronic phases of the disease. Here, we analyzed single-cell RNA-sequencing (scRNA-seq) data of 95,995 peripheral blood mononuclear cell (PBMC) from 38 ST-elevated myocardial infarction (STEMI) patients (during hospital admission, 24h and 6-8 weeks after STEMI) with those of 33,878 PBMCs from 38 age- and sex-balanced general population controls. This comparison revealed large changes in cell type composition, gene and plasma protein expression levels and cell-cell communication, both in comparison to controls and during the course of the disease (**Figure S1**). Moreover, we show that several important plasma proteins during the disease are affected by a combination of genetics, disease status or phase of the disease. Altogether, the current study increases the resolution by which the immune system during and after MI has been dissected and emphasizes the importance of taking person- and disease-related characteristics into account to fully grasp the observed underlying molecular changes. ## Methods ### Study population All patients presenting with a first STEMI upon admission at the Heart Catherization center at the University Medical Center Groningen (UMCG) were enrolled in the study between January 1st 2018 and November 30th 2019. Inclusion criteria were adults (>18 years) that had a STEMI, a primary percutaneous intervention (PCI) with implantation of at least one stent with a diameter of at least 3 mm resulting in thrombosis in myocardial infarction (TIMI) flow grade 2 or 3 post PCI, and which showed less than 6 hours symptoms before undergoing PCI. Major exclusion criteria were previous MI, medical history of diabetes, inflammatory disease or malignancies, medication affecting inflammation and clemastine or desloratadine use during intervention. Using these criteria, we initially enrolled 88 STEMI patients in this study, of which 42 were later excluded due to missing follow-up time points (t24h, t8w). Further 8 patients were excluded in the follow-up due to presence of exclusion criteria that were not yet known upon admission (**Figure 1**). The study was part of CardioLines, a single-center observational biobank aimed to study potential factors related to success or failure of diagnosis and treatment both from a patient as well as a medical perspective and was approved by the ethics committee of the UMCG, document number METC UMCG 2012/296 (17). Informed consent was obtained of all patients. scRNA-seq data from 38 age- and sex-balanced participants from the LifeLines DEEP cohort were included as a control group (18,19).  [Figure 1.](http://medrxiv.org/content/early/2023/05/03/2023.05.02.23289370/F1) Figure 1. Study overview. **a.** Flow diagram for the selection of STEMI patients (left) participating in the CardioLines Biobank and controls participating in the LifeLines DEEP biobank (right). **b.** Blood was drawn from 38 STEMI patients participating in the CardioLines Biobank at three different timepoints after STEMI (hospital admission, 24 hours, 6-8 weeks post STEMI). Plasma was isolated for Olink protein analysis and PBMCs were isolated for 10X Genomics 3’-end scRNA-seq. After sequencing, samples were demultiplexed and doublets were identified. Similarly, previously Olink and scRNA-seq data was generated for the 38 age- and sex-balanced controls from the LifeLines DEEP cohort (*19*). Cell type classification was then performed on the QCed dataset by clustering the cells per condition in STEMI patients and controls. The cell type labels were subsequently transferred back to the dataset containing all cells. ### Clinical parameters Upon inclusion of STEMI patients whole blood samples were collected at three different timepoints: during admission at the heart catheterisation center (t0), 24 hours (t24h, acute phase) and 6-8 weeks (t8w, chronic phase) later (**Figure 1**). Standard laboratory assessment of the blood was performed upon admission and routine physical parameters were recorded (**Table 2**). Titres of creatine kinase (CK), myocardial band of CK (CK-MB) and troponin T were routinely measured at 3h, 6h, 9h, 12h, 24h and 48h after admission. ### Isolation and preparation of PBMCs PBMCs were isolated and stored as previously reported (9). In short, for each donor and timepoint, whole blood was collected in two 10 ml EDTA-vacutainer (BD) tubes. Within 2 hours after collection, PBMCs were isolated from the blood using cell preparation tubes with sodium heparin (BD). PBMCs were kept in RPMI1640 supplemented with 50 µg/mL gentamicin, 2 mM l-glutamine, and 1 mM pyruvate. Isolated PBMCs were cryopreserved in RPMI1640 containing 40% FCS and 10% DMSO. Within one year, PBMCs were further processed for scRNA-seq analysis. Cells were thawed in a 37°C water bath until almost completely thawed, after which the cells were washed with pre-warmed (37°C) RPMI1640. After washing, cells were resuspended in pre-warmed RPMI1640 medium and incubated for 1 h in a 45° slant rack at 37°C in a 5% CO2 incubator. After resting, cells were washed twice in medium supplemented with 0.04% bovine serum albumin. Cells were counted using a hemocytometer and cell viability was assessed by Trypan Blue assay. ### Single-cell library preparation and sequencing Fifteen sex- and timepoint-balanced sample pools were prepared aiming for 1,000 cells per patient for a total of 8 patients per pool. Single cells were captured with the 10X Chromium controller (10x Genomics) according to the manufacturer’s instructions (document CG00026) and as described earlier (20). Each sample pool was loaded into a different lane of a 10X chip (Single Cell A Chip Kit, PN-120236 v2 and PN-1000074 v3), divided over 15 10X runs in total. cDNA libraries for v2 and v3 were generated with the Single Cell 3ʹ Library & Gel Bead kit (v2: 120237, v3: 1000075) and i7 Multiplex kit (120262) according to the company’s guidelines. These libraries were sequenced using a custom program (28-8-150 (v3) vs 27-9-150 (v2)) on an Illumina NovaSeq6000 using a 150-bp paired-end kit, per BGI (Hong Kong) sequencing guidelines. ### Alignment and initial processing of sequencing data CellRanger v3.0.2 software with default settings was used to demultiplex the sequencing data, generate FASTQ files, align the sequencing reads to the hg19 reference genome, filter cell and UMI (unique molecular identifier) barcodes, and count gene expression per cell. Genotypes of the samples were generated using the Infinium Global Screen Array-24 kit (v1 for the samples processed with 10X v2 chemistry and v3 for the samples processed with 10X v3 chemistry) and were phased using Eagle v2.4 and imputed with the HRC 1000G Phase-3 v5 (hg19) reference panel using the Michigan Imputation Server (Minimac v3) (21). Genotype information of the controls was processed as described in the Lifelines DEEP cohort paper (18). Doublet detection and sample assignment were done using Souporcell (22), confining ourselves to only use exonic single nucleotide polymorphisms (SNPs) with a MAF≥0.05. The detected doublets were discarded. Then, samples were assigned by correlating the genotypes of the clusters in each 10X lane with the genotype information of the samples in the respective lane, and assigning to the sample with the highest correlation. ### Data preparation and quality control Version 4.0 of the R package Seurat was used for data preprocessing (23–25). We performed dimensional reduction and clustering to integrate four datasets (v2 and v3 10X chemistry for STEMI patients and controls). For this, the Pearson residuals were used resulting from running SCTransform in each dataset using Seurat’s SCTintegration method (26). The first 30 principal components were used for cell clustering using Seurat’s FindNeighbors (dims = 1:30) and FindCluster function (default parameters, resolution 1.2) and an UMAP plot was used to visualize this. All genes that were not detected in at least 3 cells were removed. Cells with a high percentage of reads mapping to mitochondrially-encoded genes (>8% for 10X V2 and >15% for 10X V3) were discarded, as this can be a marker of poor-quality cells (27). Also, cells expressing ≤200 genes were considered outliers and discarded. Furthermore, cells expressing ≥ 10 UMIs of the HBB gene were discarded, indicative of erythrocytes. Finally, all cells that were marked as doublets or inconclusive by the SoupOrCell method were discarded (22). In total, 129,873 cells were used for downstream analysis consisting of 95,995 cells from STEMI patients and 33,878 cells from age- and sex-balanced controls (**Table S1**). ### Cell Type Classification Cells were annotated in each dataset separately using Seurat’s Azimuth method by projecting them on a previously annotated multimodal CITE-seq (combined scRNA-seq and protein expression) reference dataset of 162,000 PBMCs (28). These annotations were transferred to the SCTransform integrated object. Each cluster was assigned the cell type that was the most prevalent guided Azimuth annotation in that cluster. Annotations were confirmed using marker gene expression. As a result of discordant marker gene expression, two clusters were reassigned: cluster 25 from ncMono to cMono and cluster 17 from doublets to ncMono. Major cell types were defined as “B” (B naive, B intermediate, B memory), “CD4+ T-cell” (Treg, CD4+ Naive, CD4+ TCM, CD4+ TEM, CD4+ CTL, CD4+ Proliferating), “CD8+ T-cell” (MAIT, CD8+ Naive, CD8+ TCM, CD8+ TEM, CD8+ Proliferating), “DC” (ASDC, cDC1, cDC2, pDC), “HSPC” (HSPC), “monocyte” (cMono: CD14+ and ncMono: CD16+), “NK” (NK dim, NK bright, NK proliferating), “plasmablast” (plasmablast), “platelet” (platelet) and “T-other” (dnT, gdT, ILC). ### Cell type abundance Cell type proportion differences across STEMI timepoints and differences with the controls were determined using the permutation-based method described in Farbehi *et al*. (29). The number of cells of each major cell type were summed across controls or patients (separately for each STEMI timepoint). Pairwise comparisons were made across all combinations of STEMI timepoints and the controls using a W-value of 0.05. Resulting p-values were Holm-Bonferroni corrected, with a corrected p-value < 0.05 being considered statistically significant. ### Differential gene expression and pathway analysis For differential gene expression (DE) analysis library-size normalization was performed by scaling the total UMI gene expression of each cell barcode to 10,000, after which the gene expression data was log-transformed (LogNormalize in Seurat). For each cell type, pairwise DE analysis was performed using MAST comparing everything to t0 (controls, 24 hours and 6-8 weeks) (30). Due to strong technical differences in gene expression profiles between the samples that were prepared with the 10X V2 or 10X V3 protocol, DE analysis was performed separately for the samples belonging to either version of the protocol. A meta-analysis, considering only genes tested in both datasets, was performed across the V2 and V3 results using the Fisher’s Combined Probability Test. To correct for multiple testing, Bonferroni correction was applied and adjusted p-values <0.05 were considered statistically significant. The up- and downregulated DE genes were used for separate pathway analyses using the ToppFun functional enrichment tool, selecting for the REACTOME database to identify relevant pathways (31). ### Cell-to-cell communication analysis We predicted potential cell-to-cell interactions using NicheNet analysis v1.1.0 (32). Using a background set of DE genes, NicheNet predicts potential ligand-receptor interactions and predicts whether specific cell-cell interactions led to differences in downstream gene-expression. We modified NicheNet’s FindMarkers wrapper to use the MAST implementation in Seurat in order to determine DE ligands and receptors, as well as genes downstream of cells potentially interacting with the ligands. Differential ligand and downstream expression changes were confined to genes that were present in at least 10% of the cells. The absolute log fold change between the conditions needed to be at least 0.1 to be considered for testing. We assessed potential ligand-receptor interactions between pairs found in the Omnipath database (33). The Nichenet supplied ligand-receptor network data was used to determine what could be considered to be cell communication patterns. Predicted ligand-receptor interactions were filtered for ligands that were shown to have an effect on downstream gene expression. Ligand-receptor connections between cell types were counted and visualized in circle plots. ### Protein Olink analysis Plasma from the 38 STEMI patients was collected at each of the three timepoints (t0, t24h, t8w), but missing one sample at t8w (n=37), and stored at −80°C. Plasma samples were analyzed for 92 cardiovascular disease related protein biomarkers by Olink Proteomics (Olink Target 96 Cardiovascular III panel, Uppsala Sweden). The data received was QCed by excluding all proteins with a missing data frequency of >15% and samples with an Olink warning QC outcome from analysis. All samples with a low detection rate were given the limit of protein detection value for the specific protein. Using these QC parameters, none of the proteins had to be excluded. Protein quantitative trait loci (pQTL) analysis was done using our in-house eQTL mapping pipeline [[https://github.com/molgenis/systemsgenetics/wiki/QTL-mapping-pipeline](https://github.com/molgenis/systemsgenetics/wiki/QTL-mapping-pipeline)]. pQTLs were mapped per time point, only testing the SNPs within a 100k window of the transcription start site of the gene transcribing the tested protein. ### Relating clinical variables to cell type abundance, gene and protein expression Cell type proportions and gene and protein expression profiles were associated with peak CK-MB values in STEMI patients. For gene expression associations we used MAST, a two-part generalized linear model on the log-transformed count matrix, and correcting for age, sex, 10X chemistry and 10X lane. For cell type proportion and protein expression associations, we first conducted a centered-log-ratio or log-transformation, respectively, after which linear regression was conducted, correcting for sex and age. ### pQTL analysis For STEMI patients individual Olink plasma protein levels were linked to specific genetic variants (pQTL mapping) at each of the 3 timepoints (t0, t24h, t8w) using the *eQTL mapping pipeline* ([https://github.com/molgenis/systemsgenetics](https://github.com/molgenis/systemsgenetics)). We limited ourselves to the 241 variants that were previously found to show a genome-wide significant GWAS signal for coronary artery disease (CAD) (34). The significant pQTLs in STEMI patients in at least one timepoint were then assessed in controls by mapping these pQTLs in 1142 individuals from the general population cohort Lifelines Deep, in which the same Cardiovascular Olink panel was previously measured (35) and genotype data was available. Genetic variants were converted to Trityper format using GenotypeHarmonizer ([https://github.com/molgenis/systemsgenetics](https://github.com/molgenis/systemsgenetics)). To assess whether the pQTL effect strength was influenced by the timepoint after the STEMI (t0, t24h, t8w), we then conducted an interaction analysis by fitting a generalized linear model correcting for age and sex as fixed effects, and the participant as a random effect. ### Statistical analysis Discrete variables were represented as frequencies and percentages. Continuous variables with a normal distribution were summarized as mean ± standard deviation and if skewed were represented as median with interquartile range. Spearman correlation was used for non-normally distributed variables. For DE gene and protein analysis, and their correlation to peak CK-MB values, multiple testing correction was performed using Holm-Bonferroni. We considered a corrected P-value <0.05 statistically significant. Statistical analysis and the figures were created using R Core Team (2020). ## Results ### Population characteristics To dissect the immune response during and after a MI, PBMCs and plasma from 38 patients with a first STEMI were collected at the time of hospital admission (t0), and 24h (t24h) and 6-8 weeks (t8w) after the event (**Figure 1**). The patient characteristics during hospital admission are described in **Table 1**. On average, STEMI patients were 60±11 years, had a median BMI of 26.4 (SD 3.5) kg/m2 and the majority were males (84%). Most patients presented with one stenosed vessel (55%), and in 39% of the patients the affected vessel was completely occluded (TIMI 0 flow). These 38 STEMI patients were compared to 38 age- and sex-balanced controls from the LifeLines DEEP cohort whose scRNA-seq data was previously published (19). These controls showed comparable age- (59.2 years) and sex-characteristics (84% male), but a slightly lower BMI (25.4 (SD 3.3) kg/m2) than the STEMI patients (**Table 1**). View this table: [Table 1.](http://medrxiv.org/content/early/2023/05/03/2023.05.02.23289370/T1) Table 1. t0 characteristics of STEMI patients and controls from the LifeLines DEEP Cohort. ### Single-cell profiling of immune cells in STEMI patients and controls The collected PBMCs from the STEMI patients at three different timepoints post-STEMI (t0, t24h, t8w) and controls from the LifeLines cohort (19), were then used for 10X Genomics scRNA-seq analysis. A combination of 10X Genomics v2 and v3 chemistry reagents were used to capture an average of 842 cells per individual per condition (v2: 831 genes/cell, v3: 1,548 genes/cell) in STEMI patients, and an average of 891 cells per individual in controls (v2: 1,012 genes/cell, v3 1,931 genes/cell) (**Table S2**). Souporcell was used to demultiplex and to identify the doublets coming from different individuals (22). Souporcell clusters were correlated to individuals by correlating the cluster genotypes with donor genotypes using Spearman correlation. This revealed that by mistake we omitted loading cells for one donor at the t8w condition and that on average, 10.0% of the cells were assigned as doublet or inconclusive. Due to technical differences between v2 and v3 chemistry, determination of quality control (QC) thresholds and analyses were performed separately per chemistry. Low quality cells were excluded, leaving 129,873 cells (95,995 diseased and 33,878 control) in the final dataset used for analyses (see Methods). UMAP dimensionality reduction and KNN-clustering was then applied on the normalized, integrated count data, allowing the identification of cell types. ### Cell type abundance We identified 10 major cell types in the PBMC fraction of STEMI patients and controls, including the B, CD4+ T, CD8+ T, DC, HSPC, monocyte, NK, plasmablast, platelet and other T cells (**Figure 2a**). At higher resolution, these major cell types could be further splitted in 29 minor cell populations, including different subsets of T cells (regulatory T cells (Treg), CD4+ and CD8+ T-naive, CD4+ and CD8+ T central memory (TCM), CD4+ and CD8+ T effector memory (TEM), CD4+ cytotoxic (CTL), CD4+ and CD8+ proliferating, Mucosal associated invariant T cell (MAIT), double negative T (dnT), gamma-delta T (gdT), innate lymphoid cells (ILC), B cells (B naive, B intermediate, B memory), NK cells (NK dim, NK bright, NK Proliferating), monocytes (classical (cMono) and non-classical monocytes (ncMono)) and dendritic cells (AXL+ SIGLEC+ Dendritic Cell (ASDC), conventional DC1 (cDC1), cDC2, plasmacytoid DC (pDC)) (**Figure 2b**).  [Figure 2.](http://medrxiv.org/content/early/2023/05/03/2023.05.02.23289370/F2) Figure 2. Cell type composition in STEMI patients over time and in controls. **a.** UMAP showing the composition of PBMCs with 10 major cell types in STEMI patients and controls. **b.** UMAP of the 29 minor cell clusters in STEMI patients and controls. **c**. Proportions of monocytes and NK cells and their subcell types, over time (0 hours - t0, 24 hours - t24h, and 6-8 weeks - t8w) and in controls (C). cMono; classical monocytes, ncMono; non-classical monocytes, NK dim and NKbright. Significant differences are denoted as *****P < 0.05. To conduct a robust cell type composition analysis between the different conditions (STEMI patients at t0 compared to controls, STEMI patients over time), we initially focused our analyses on the 6 most abundant cell types (≥1% of the total PBMCs in each of the tested conditions): B cells, CD4+ T cells, CD8+ T cells, DCs, monocytes and NK cells (**Figure 2c, S1, Table S2, S3)**. After which, we zoomed in at the higher resolution subcell type compositional changes to define the subtypes that contributed to these changes in composition (**Figure S3**). Compared to controls, monocytes from STEMI patients during admission showed a relative increase (p=2.0×10-5) and NK cells a relative decrease in their abundance (p=9.6×10-3) (**Figure 2c, Table S3**). Zooming in further, we observed that these changes could be attributed to compositional changes of the cMono (p=2.0×10-5) and the NK dim (p=1.3×10-2) subcell types (**Figure 2c, Table S3**). Similarly, in STEMI patients during the course of their disease (t0, t24h, t8w), we observed relative composition changes in the monocytes (decrease t0-t8w: p<5.0×10-3) and NK cells (decrease t0-t24h: p=9.0×10-3) (**Figure 2c**), but also in the CD4+ T cells (increase: t0-t8w: p=4.0×10-3) (**Figure S2**). These were again attributed due to changes in the cMonos (p=2.0×10-5) and NK dim cells (p=1.1×10-2), whereas none of the other minor cell types showed changes in composition during the disease course (**Figure S3**). ### Gene expression profiles of circulating immune cells in STEMI patients and controls To assess gene expression differences in STEMI patients at time of hospital admission versus controls, we performed DE gene analysis on the six major cell types using MAST (**Table S4**). Of all cell types, monocytes (846 DE genes), followed by NK cells (553 DE genes) showed the largest and B cells the least amount of DE genes (130 DE genes) (**Figure 3a**). The majority of the DE genes were uniquely identified in only one major cell type (**Figure 3b**). Looking at the monocytes, the cell type with most total and unique DE genes, we observed pathway enrichment among upregulated genes in STEMI patients involving various pro-inflammatory pathways, including IFN signaling, interleukins (IL-1, IL-6) and all toll-like receptors (**Table S5**).  [Figure 3.](http://medrxiv.org/content/early/2023/05/03/2023.05.02.23289370/F3) Figure 3. Differential expressed genes and pathways in STEMI patients and in controls. **a.** The number of DE genes per cell type in STEMI patients at hospital admission (t0) compared to controls (C). **b**. Bar plot showing the overlap of DE genes in cell types at t0 compared to Controls. **c**. The number of DE genes per cell type in STEMI patients at 24 hours (t24h) or 6-8 weeks (t8w) compared to t0. **d**. Bar plot showing the overlap of DE genes in cell types in STEMI patients over time, taking all the DE genes that were significant in at least one of the comparisons (t0-t24h or t0-t8w). Subsequent comparison of the gene expression levels of STEMI patients at hospital admission compared to later time points revealed the largest amount of DE genes in NK cells and in monocytes, respectively 24h and 6-8 weeks after STEMI (**Figure 3c, Table S4**). Specifically looking at the pathways underlying the upregulated genes at these cell type-timepoint combinations, we observed that many of the NK DE genes at t24h were involved in interactions with platelets (activation, degranulation, response to platelet cytosolic Ca2+), whereas many of the monocyte DE genes at t8w were involved in pro-inflammatory signaling (IFN-gamma, IL-1) (**Table S5**). Remarkably, all major types showed a decrease in DE genes going from 24h to 6-8 weeks after STEMI, except for the monocytes (**Figure 3c**). And as before in comparison to controls (**Figure 3b**), monocytes also showed the largest proportion of unique DE genes over the course of the disease (**Figure 3d**). Together with, but independently of, the observed cell type composition changes (**Figure 2c**), these results indicate that both the monocytes and NK cells seem to play a pivotal role in the immunological changes that are occurring upon a STEMI. ### Cell-to-cell communication in STEMI patients over time and controls Next, we defined how cell-to-cell communication interactions among immune cell types could underlie these observed gene expression changes in STEMI patients over time. For this, we used NicheNet, a cell-to-cell communication tool that identifies potential cell-cell interactions by assuming that ligand-receptor interactions can be predicted based on the expression of the ligand in one cell type and downstream gene expression changes of a known ligand-receptor interaction in another cell type (**Table S6**). This analysis revealed that in the acute phase of the disease (t0-t24h) in- and outcoming communication was balanced among the cell types (**Figure 4a**). In contrast, in the chronic phase of the disease (t24h-t8w) the T cells and DCs were mostly sending communication, while monocytes and B cells were mostly receiving it (**Figure 4b**). Although from the comparable number of incoming interactions for a specific cell type it may seem that the involved ligands are mostly shared among cell types, most of these ligands were actually uniquely involved in the communication by just one cell type (**Figure 4c, 4d**). Looking at the number of ligand-receptor interactions that were associated with changes in downstream signaling, we observed both NK and CD8+ T cells to be the largest communicators during the acute phase (**Figure 4a**), whereas monocytes were the main communicators during the chronic phase of the disease (**Figure 4b**). These results align well with our DE analyses, which also indicated the largest impact on expression changes of STEMI on NK cells early and on monocytes later in the disease.  [Figure 4.](http://medrxiv.org/content/early/2023/05/03/2023.05.02.23289370/F4) Figure 4. Differential cell-to-cell communication after STEMI. Barplots depicting the differential incoming (receiver) and outgoing (sender) cell-to-cell communication at the t0-t24h (**a**) or t24h-t8w (**b**) period. An active cell-to-cell communication link is counted as being a ligand-receptor (LR) link that has resulted in differential downstream gene expression. For each active LR link in **a** (t0-t24h) and **b** (t24h-t8w), the sharedness of ligands among the major cell types is depicted in **c** and **d**, respectively. ### Protein expression profiles in STEMI patients over time While our previous analyses at the mRNA level help us to disentangle the molecular processes that immune cells undergo after STEMI, changes in plasma protein levels may provide us additional insights on the systemic consequences of a STEMI and can be more easily monitored in clinical follow-up of patients. Therefore, next, using the Olink Target 96 Cardiovascular III panel, we analyzed the plasma levels of 92 proteins that are known or exploratory human cardiovascular and inflammatory markers in each of the three timepoints (t0, t24h, t8w) within STEMI patients. Targeted differential protein analysis revealed that 14 out of 92 assessed plasma protein levels changed during the first 24h (t0-t24h: 8 up, 6 down) and 28 out of 92 during the first 8 weeks (t0-t8w: 22 up, 6 down) after STEMI. During the first 24h, the top three upregulated proteins were NT-proBNP (mean LFC= 0.74, p=2.1×10-7), IL1RL1 (mean LFC=0.18, p=7.3×10-4), and CHI3L1 (mean LFC=0.16, p=6.7×10-4) (**Figure 5**, list of all DE proteins shown in **Table S7**). Both upregulation of NT-proBNP and the soluble form of IL1RL1, a member of the IL-1 receptor family, are well known proteins that are upregulated in response due to increased wall stress during STEMI (36). Moreover, both are independent predictors of heart failure and cardiovascular death (37,38). CHI3L1, also known as YKL-40, is an extracellular matrix protein which is involved in several stages of CAD. It induces maturation of monocytes to macrophages in the early phase of atherosclerosis, but also promotes atherosclerotic progression and complications like plaque rupture (39). So, the upregulation of this protein might be due to the cause of MI (i.e. plaque rupture), whereas NT-proBNP and IL1RL1 reflect the mechanical stress to the cardiomyocytes due to MI.  [Figure 5.](http://medrxiv.org/content/early/2023/05/03/2023.05.02.23289370/F5) Figure 5. Differentially expressed proteins after STEMI. The top 3 most up- and downregulated plasma proteins 24h (t24h) or 6-8 weeks (t8w) after STEMI compared to time of hospital admission (t0). The top three downregulated proteins during the acute phase were SPON1 (mean LFC=-1.4, p=6.7×10-10), AZU1 (mean LFC=-0.44, p=3.3×10-6) and IGFBP1 (mean LFC=-0.43, p=9.9×10-6). Each of these proteins have been previously associated with processes preceding a STEMI (e.g. atherosclerosis for AZU1 (40) and IGFBP1 (41,42)) or that are the consequence of unsuccessful STEMI treatment (e.g. worsened systolic heart function for SPON1 (43). Therefore, we suspect their downregulation during the acute and chronic phase may be an indication of effective treatment by PCI after STEMI, resulting in restoration of blood flow and eventually preservation of heart function. During the chronic phase, the top three upregulated proteins were NT-proBNP (LFC= 0.3, p= 1.3×10-2), PCSK9 (mean LFC=0.18, p=6.0×10-2) and MMP2 (mean LFC=0.16, p=2.6×10-3). The increased NT-proBNP levels may reflect the further increased wall stress in the chronic phase after STEMI. PCSK9 plays an important role in degradation of the LDL-receptor (LDLR) (44) and the observed chronic increase of PCSK9 might therefore reflect an ongoing state of hyperlipidemia. MMP2 is a matrix metalloproteinase-2 that acts in the fibrotic pathway (45), and is involved in cardiac remodeling post STEMI (46). The top three downregulated proteins during the chronic phase were SPON1 (mean LFC=-1.2, p=1.3×10-9), MPO (mean LFC=-0.36, p=2.0×10-6) and AZU1 (mean LFC=-0.35, p=5.6×10-4). Similar to the acute phase, both SPON1 and AZU1 continued to be downregulated. Additionally, MPO was found downregulated, a protein involved in the acute innate inflammatory response which contributes to plaque destabilization through local oxidative tissue injury (47,48). Together, these changes seem to indicate the long-term restoration of heart function and blood flow by PCI. ### Relating Peak CK-MB as proxy for biochemical infarct size STEMI patients can already widely differ from each other when they are presented with symptoms in the clinic. While the patient inclusion and exclusion criteria of our study (**Methods**) reduced the amount of clinical variation in our study (**Table 1**), many aspects remain that could underlie the observed donor-to-donor variation in the measured molecular phenotypes, including age, sex, infarct size and genetics. To assess whether variation in infarct size may correlate with the individual molecular data layers, we used plasma peak CK-MB levels as proxy. In clinical practice, plasma CK-MB levels are measured regularly after STEMI to monitor myocardial damage. This biomarker reaches a peak within 24 hours, and the peak value is an estimation of the infarct size and predicts left ventricle dysfunction (49). First, we used plasma peak CK-MB levels to assess whether PBMC cell type proportions and monocyte gene expression levels were associated with biochemical infarct size. We focussed on gene expression of the monocytes as this was the cell type with the largest amount of total and unique DE genes. However, neither cell type proportions, nor monocyte gene expression levels showed a significant association with peak CK-MB at any of the tested timepoints (t0, t24h or t8w). Next, we assessed whether peak CK-MB values were associated with inflammatory plasma protein levels in the acute and chronic phase after STEMI using simple linear regression. Of all 92 proteins, only NT-proBNP was positively associated with peak CK-MB in the acute phase (t24h: r2 adjusted=0.49, p=2.8×10-3) (**Table S8**). NT-proBNP, a marker for myocardial wall stress, is commonly measured in clinical practice to assess left ventricular dysfunction (50,51). As there is a direct relationship between the level of myocardial damage and the amount of left ventricular dysfunction, NT-proBNP and CK-MB levels are expected to be correlated during the time CK-MB peaks (∼24h after STEMI) (36,52). ### QTL analysis Beyond infarct size, genetic variation among patients, can be an important contributor to the observed differences in molecular response among patients. Therefore, it is important to take genetic variation into account when assessing a patient’s disease course and defining patient-tailored treatment. As the measured plasma protein levels provide a pool of potential direct therapeutic targets, we assessed the effect of genetic variation on these (**Table S9**). To provide a direct clinical link, we specifically focussed on those genetic variants that were previously found to be associated with CAD risk (34). For three out of 92 plasma proteins, we detected a significant pQTL in the STEMI patients in at least one of the three timepoints (t0, t24h, t8w). For these three pQTLs, we assessed their disease-specificity by comparing pQTL effect sizes with those in a control cohort of the general population. This revealed one pQTL effect that was significant in all conditions and whose effect size was comparable among both patients (independent of time point) and controls: SNP rs6686750 affecting IL6R plasma protein levels (STEMI patients: genotype β=0.42, p=1.6×10-5, **Figure 6a**). The two other pQTLs were disease- or even condition-specific: SNP rs8124182 affecting MMP9 plasma protein levels was found specifically in STEMI patients (STEMI patients: genotype β=-0.68, p=0.016, **Figure 6b**) and SNP rs10422256 affecting LDLR plasma protein levels was found only in STEMI patients at t24h (STEMI patients: genotype x t24h interaction β=-0.29, p=0.022, **Figure 6c**). These results indicate the importance of taking genotype into consideration, when studying the relationship between STEMI and molecular phenotype.  [Figure 6.](http://medrxiv.org/content/early/2023/05/03/2023.05.02.23289370/F6) Figure 6. Plasma protein QTLs in STEMI patients versus controls. Plasma pQTLs were measured in the 38 STEMI patients at time of hospital admission (t0), 24 hours (t24h) or 6-8 weeks (t8w) after STEMI and compared to 1142 controls from the general population cohort Lifelines Deep (18). Three out of 92 plasma proteins showed a pQTL effect that was significant in all conditions (IL6R), specific for STEMI patients (MMP9) or even specific for one timepoint in the STEMI patients (t24h: IL6R). ## Discussion and Conclusion In this study, we have provided an unbiased longitudinal overview of the single-cell immune response in STEMI patients and compared this to age- and sex-balanced controls. We first mapped changes in the cell type composition and gene and plasma protein expression level. We then determined the overarching pathways that were enriched and the cell types being involved in these changes as a result of cell-to-cell communication. Thereafter, we determined whether these changes were correlated with an individual’s (biochemically-measured) infarct size, genotype or disease phase. Together, all these analyses indicated an important role for the classical monocytes and NKdim cells during STEMI, both in terms of cell type compositional and gene expression changes. Our observation that monocytes play a pivotal role in STEMI is in line with a previous scRNA-seq study of 10 STEMI patients (with and without plaque rupture) that were analyzed at time of hospital admission (53). While this study did not observe any differences in cell type composition between STEMI patients with or without plaque rupture, they did observe that most cell-cell communication occurred from and to monocytes. Our study now adds additional insights, because of the longitudinal character, the added layer of protein information, the increased sample size, and the availability of controls to compare with. The recently published spatial multi-omic map of human MI, in which cardiac structure and specific cardiac cells were studied in spatial context, underlined the role of immune cells in cardiac repair following STEMI (16). While that study used cardiac tissue, and we studied the response of circulating immune cells, it remains clear that gaining more insights into the immune response after STEMI is key for future drug therapies. One such insight relevant for future drug therapies that our study has now provided, is in relation to the use of IL-1ẞ blockers in MI. Our study showed chronic upregulation of the IL-1 signaling pathway specifically in the monocytes of STEMI patients. The CANTOS trial showed that IL-1ẞ inhibition with canakinumab led to lower cardiovascular event rates, but also higher fatal infections (4). Our insights into the IL-1 pathway could aid in targeting specifically monocytes, and thereby potentially preventing the occurrence of serious side effects. Furthermore, a study in IL-1R knock-out mice showed that these mice, which cannot respond to IL-1, were protected from hypertrophy, fibrosis, and diastolic dysfunction (56). This may indicate that the observed chronic upregulation of the IL-1 signaling pathway after STEMI reflects an already ongoing process of developing heart failure and eludes the potential to broaden the use of IL-1 inhibition for both atherosclerotic events and fibrotic processes. Additionally, our study indicated the importance of considering both genetic variation and disease phase when assessing the molecular consequences of a STEMI, as several of the protein expression levels were affected by a combination of these parameters. Importantly, by combining molecular layers, our integrative analysis revealed insights in several therapeutic targets for which drugs are already available and are currently tested in clinical trials. The first example for a promising drug target in the context of CVD, involves targeting IL-6 signaling. IL-6 signaling can be activated through three different modes of signaling, each having their own consequences: classical signaling, trans-signaling and trans-presentation (54). The anti-inflammatory classical signaling is activated by membrane-bound IL6R and gp130 (also known as IL6ST) within the same cell, the pro-inflammatory trans-signaling is regulated by soluble IL-6R and membrane-bound gp130, and the Th17-cell promoting trans-presentation signaling is induced by membrane-bound IL-6R (in DCs) and gp130 (in Th17 cells). Through cell-cell communication analyses with NicheNet (**Table S6**), we found that both the classical and trans-presentation pathways were differentially active during the first 24h of STEMI. While our pQTL analysis indicated that the trans-signaling pathway was affected by SNP rs6689206 regulating soluble IL6R levels both in STEMI patients as well as in controls (Fig. 6A). Each of these IL-6 pathways can be targeted using different drugs (54), and based on our findings it seems that these drugs may be most beneficial during specific stages of the disease (first 24h after STEMI: for classical and trans-presentation IL-6 signaling, using drugs like anti-IL-6 (54,55) or anti-IL-6R) or in donors with a specific genetic background (patients with the GG genotype at SNP rs6689206 have the highest expression of IL-6R and thereby will naturally activate IL-6 trans signaling the most, and as such are expected to benefit most of treatments targeting specifically this pathway, e.g. Olamkicept, soluble gp130) (Fig 6A). Another interesting drug target for which we now provide additional insights is the LDLR. LDLR plasma levels act as scavengers for plasma LDL and can be indirectly targeted by PCSK9 inhibitors, which inhibit the degradation of membrane LDLR (44). As such, they are widely used to decrease LDL cholesterol and reduce the risk of cardiovascular events (44). Our DE protein analyses revealed that PCSK9 levels increased after STEMI, while our pQTL analysis revealed that 24h after STEMI individuals with the GG genotype (at SNP rs10422256) have a lower plasma protein level of LDLR independent of their PCSK9 levels and may therefore have higher free plasma LDL leading to increased cardiovascular disease risk (Fig. 6C). Together, this indicates that individuals with the GG genotype at SNP rs10422256 may benefit the most from PCSK9 inhibitor treatment, especially in the first 24 hours of the disease. We expect that implementing such genetic and phase of the disease information will be important to maximize benefit of such patient-tailored therapies in future clinical practice. Unfortunately, when focusing on how each of the molecular layers could be contributing to or being the consequence of biochemically-measured infarct size, we could not identify any direct relationship except for plasma NT-proBNP levels – a known marker for heart failure (50,51). We suspect this lack of association might have been the result of assessing each molecular parameter separately or not considering all potentially relevant clinical and donor variables (age, sex, smoking status, genetics, etc.) into the model57. To enable such a model considering all these parameters, a larger sample size would be required. We expect such approaches to be feasible in the near future by performing meta-analyses across multiple (disease) population-based single-cell datasets that are uniformly processed in efforts such as those of the single-cell eQTLGen consortium (10). In conclusion, this study highlights the importance of studying STEMI at cell-type-specific resolution, while taking genetic variation, disease status and phase of the disease into consideration. We expect that such an integrative approach as used here, will help to better grasp the molecular processes underlying STEMI and will be essential for the development of effective future therapies with reduced side-effects. ## Supporting information Supplementary Figure 1 [[supplements/289370_file08.pdf]](pending:yes) Supplementary Figure 2 [[supplements/289370_file09.pdf]](pending:yes) Supplementary Figure 3 [[supplements/289370_file10.pdf]](pending:yes) Supplementary Table 1 [[supplements/289370_file11.xlsx]](pending:yes) Supplementary Table 2 [[supplements/289370_file12.xls]](pending:yes) Supplementary Table 3 [[supplements/289370_file13.xlsx]](pending:yes) Supplementary Table 4 [[supplements/289370_file14.xlsx]](pending:yes) Supplementary Table 5 [[supplements/289370_file15.xlsx]](pending:yes) Supplementary Table 6 [[supplements/289370_file16.xls]](pending:yes) Supplementary Table 7 [[supplements/289370_file17.xlsx]](pending:yes) Supplementary Table 8 [[supplements/289370_file18.xls]](pending:yes) Supplementary Table 9 [[supplements/289370_file19.xlsx]](pending:yes) ## Data Availability Processed (de-anonymized) scRNA-seq data, including a text file that links each cell barcode to its respective individual, has been deposited at the European Genome-Phenome Archive (EGA), which is hosted by the EBI and the CRG, under accession number EGAD00001010064. ## Sources of Funding LF is supported by the Dutch Organisation for Scientific Research (NWO) Corona Fast-Track grant (440.20.001), an Oncode Senior Investigator grant, NWO VIDI grant (917.14.374) and NWO VICI grant (09150182010019). MW is supported by a NWO VENI grant (192.029). ## Disclosures None. ## Data availability Processed (de-anonymized) scRNA-seq data, including a text file that links each cell barcode to its respective individual, has been deposited at the European Genome-Phenome Archive (EGA), which is hosted by the EBI and the CRG, under accession number EGAD00001010064. ## Code availability The code for Seurat v4 ([https://github.com/satijalab/seurat](https://github.com/satijalab/seurat)), Eagle v2.x ([https://github.com/poruloh/Eagle](https://github.com/poruloh/Eagle)), Souporcell v1.x ([https://github.com/wheaton5/souporcell](https://github.com/wheaton5/souporcell)) and our in-house eQTL pipeline2 v1.4.0 ([https://github.com/molgenis/systemsgenetics/tree/master/eqtl-mapping-pipeline](https://github.com/molgenis/systemsgenetics/tree/master/eqtl-mapping-pipeline)) can be found at GitHub. All custom code is also available via GitHub ([https://github.com/molgenis/STEMI-scRNA-seq](https://github.com/molgenis/STEMI-scRNA-seq)). ## Supplementary figures  [Figure S1:](http://medrxiv.org/content/early/2023/05/03/2023.05.02.23289370/F7) Figure S1: Overview of the study. Describing the cohort and the results per section  [Figure S2:](http://medrxiv.org/content/early/2023/05/03/2023.05.02.23289370/F8) Figure S2: Proportions of major cell types in STEMI patients compared to controls (Cs) and over time at three different time points (t=0, t=24 hours and t=6-8 weeks post STEMI). Significant differences are denoted as *P < 0.05.  [Figure S3:](http://medrxiv.org/content/early/2023/05/03/2023.05.02.23289370/F9) Figure S3: Proportions of minor cell types in STEMI patients compared to controls (C) and over time at three different time points (t=0, t=24 hours and t=6-8 weeks). ASDC; Axlp Sigle6p Dendritic Cell, CTL; Cytotoxic cell, TCM; T-central memory, TEM; T-effector memory, cDC; conventional dendritic cell, pDC; plasmacytoid dendritic cell, HSPC; Hematopoietic stem and progenitor cell, ILC; Innate lymphoid cell, MAIT; Mucosal associated invariant T cell, cMono; classical monocytes, ncMono; non-classical monocytes. Significant differences are denoted as *P < 0.05. ## Supplementary tables **Table S1:** QC/preprocessing cell counts **Table S2:** Absolute cell counts per donor, split by major cell type and STEMI timepoint **Table S3:** eLIFE output DPA **Table S4:** DE mRNA analysis **Table S5:** Pathway analysis **Table S6:** Cell-cell communication output of L-R and L-T interactions **Table S7:** Differential Protein analysis **Table S8:** peak CK-MB by protein expression **Table S9:** pQTL summary stats ## Acknowledgements Acknowledgements to figures: Images used for the figures were downloaded from [https://smart.servier.com/](https://smart.servier.com/). Original images were used, no changes were made. The license for using the images is Creative Commons License, Attribution 3.0 Unported License. Acknowledgements to Redcap: [https://projectredcap.org/resources/citations/](https://projectredcap.org/resources/citations/) ## Footnotes * typo in one of the author names has been fixed ## Abbreviations AZU1 : Azurocidin CAD : Coronary artery disease CHI3L1 : Chitinase-3-like-1 CITE-seq : Cellular indexing of transcriptomes and epitopes by sequencing CK : Creatine kinase CK-MB : Creatine kinase, myocardial band cMono : Classical monocyte EDTA : Ethylenediaminetetraacetic acid HbA1c : Glycated hemoglobin C : Control IGFBP1 : Insulin like growth factor binding protein 1 IL-1 : Interleukin 1 IL6R : Interleukin 6 receptor IL1RL1 : Interleukin 1 receptor-like 1 LDL : Low-density lipoprotein LFC : Log fold change MAF : Minor allele frequency MAST : Model-based analysis of single cell transcriptomics MCV : Mean corpuscular volume METC : Medical ethical review committee MI : Acute myocardial infarction MMP2 : Matrix metalloproteinase-2 MPO : Myeloperoxidase ncMono : non-classical monocyte NT-proBNP : N-terminal pro-brain natriuretic peptide PBMC : Peripheral blood mononuclear cell PCI : Percutaneous coronary intervention PCSK9 : Proprotein convertase subtilisin/kexin type 9 pQTL : Protein quantitative trait locus QC : Quality check scATAC-seq : Single-cell assay for transposase-accessible chromatin sequencing SNP : Single-nucleotide polymorphism scRNA-seq : Single-cell RNA-sequencing snRNA-seq : Single-nucleus RNA-sequencing SPON1 : Spondin-1 STEMI : ST-elevated myocardial infarction TIMI : Thrombosis in myocardial infarction UMAP : Uniform manifold approximation and projection UMCG : University medical center groningen UMI : Unique molecular identifier VD : Vessel disease * Received May 2, 2023. * Revision received May 3, 2023. * Accepted May 3, 2023. * © 2023, Posted by Cold Spring Harbor Laboratory The copyright holder for this pre-print is the author. All rights reserved. The material may not be redistributed, re-used or adapted without the author's permission. ## References 1. 1.Libby P, Nahrendorf M, Swirski FK. Leukocytes Link Local and Systemic Inflammation in Ischemic Cardiovascular Disease. J Am Coll Cardiol. 2016 Mar;67(9):1091–103. [FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6MzoiUERGIjtzOjExOiJqb3VybmFsQ29kZSI7czo0OiJhY2NqIjtzOjU6InJlc2lkIjtzOjk6IjY3LzkvMTA5MSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIzLzA1LzAzLzIwMjMuMDUuMDIuMjMyODkzNzAuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 2. 2.Aday AW, Ridker PM. Targeting Residual Inflammatory Risk: A Shifting Paradigm for Atherosclerotic Disease. Front Cardiovasc Med. 2019 Feb 28;6:16. 3. 3.Lawler PR, Bhatt DL, Godoy LC, Lüscher TF, Bonow RO, Verma S, et al. Targeting cardiovascular inflammation: next steps in clinical translation. Eur Heart J. 2020 Mar 16;ehaa099. 4. 4.Ridker PM, Everett BM, Thuren T, MacFadyen JG, Chang WH, Ballantyne C, et al. Antiinflammatory Therapy with Canakinumab for Atherosclerotic Disease. N Engl J Med. 2017 Sep 21;377(12):1119–31. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa1707914&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28845751&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F03%2F2023.05.02.23289370.atom) 5. 5.Broch K, Anstensrud AK, Woxholt S, Sharma K, Tøllefsen IM, Bendz B, et al. Randomized Trial of Interleukin-6 Receptor Inhibition in Patients With Acute ST-Segment Elevation Myocardial Infarction. J Am Coll Cardiol. 2021 Apr;77(15):1845–55. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F03%2F2023.05.02.23289370.atom) 6. 6.Prabhu SD, Frangogiannis NG. The Biological Basis for Cardiac Repair After Myocardial Infarction: From Inflammation to Fibrosis. Circ Res. 2016 Jun 24;119(1):91–112. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTA6ImNpcmNyZXNhaGEiO3M6NToicmVzaWQiO3M6ODoiMTE5LzEvOTEiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMy8wNS8wMy8yMDIzLjA1LjAyLjIzMjg5MzcwLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 7. 7.Mouton AJ, DeLeon-Pennell KY, Rivera Gonzalez OJ, Flynn ER, Freeman TC, Saucerman JJ, et al. Mapping macrophage polarization over the myocardial infarction time continuum. Basic Res Cardiol. 2018 Jul;113(4):26. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00395-018-0686-x&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29868933&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F03%2F2023.05.02.23289370.atom) 8. 8.Peet C, Ivetic A, Bromage DI, Shah AM. Cardiac monocytes and macrophages after myocardial infarction. Cardiovasc Res. 2020 May 1;116(6):1101–12. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/cvr/cvz336&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F03%2F2023.05.02.23289370.atom) 9. 9.van der Wijst MGP, Brugge H, de Vries DH, Deelen P, Swertz MA, Franke L. Single-cell RNA sequencing identifies celltype-specific cis-eQTLs and co-expression QTLs. Nat Genet. 2018 Apr;50(4):493–7. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-018-0089-9&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F03%2F2023.05.02.23289370.atom) 10. 10.van der Wijst M, de Vries D, Groot H, Trynka G, Hon C, Bonder M, et al. The single-cell eQTLGen consortium. eLife. 2020 Mar 9;9:e52155. 11. 11.Tang F, Barbacioru C, Wang Y, Nordman E, Lee C, Xu N, et al. mRNA-Seq whole-transcriptome analysis of a single cell. Nat Methods. 2009 May;6(5):377–82. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nmeth.1315&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19349980&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F03%2F2023.05.02.23289370.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000265661600019&link_type=ISI) 12. 12.Litviňuková M, Talavera-López C, Maatz H, Reichart D, Worth CL, Lindberg EL, et al. Cells of the adult human heart. Nature. 2020 Dec 1;588(7838):466–72. 13. 13.Winkels H, Ehinger E, Vassallo M, Buscher K, Dinh HQ, Kobiyama K, et al. Atlas of the Immune Cell Repertoire in Mouse Atherosclerosis Defined by Single-Cell RNA-Sequencing and Mass Cytometry. Circ Res. 2018 Jun 8;122(12):1675–88. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTA6ImNpcmNyZXNhaGEiO3M6NToicmVzaWQiO3M6MTE6IjEyMi8xMi8xNjc1IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMDUvMDMvMjAyMy4wNS4wMi4yMzI4OTM3MC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 14. 14.Cochain C, Vafadarnejad E, Arampatzi P, Pelisek J, Winkels H, Ley K, et al. Single-Cell RNA-Seq Reveals the Transcriptional Landscape and Heterogeneity of Aortic Macrophages in Murine Atherosclerosis. Circ Res. 2018 Jun 8;122(12):1661–74. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTA6ImNpcmNyZXNhaGEiO3M6NToicmVzaWQiO3M6MTE6IjEyMi8xMi8xNjYxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMDUvMDMvMjAyMy4wNS4wMi4yMzI4OTM3MC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 15. 15.Li Z, Solomonidis EG, Meloni M, Taylor RS, Duffin R, Dobie R, et al. Single-cell transcriptome analyses reveal novel targets modulating cardiac neovascularization by resident endothelial cells following myocardial infarction. Eur Heart J. 2019 Aug 7;40(30):2507–20. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/eurheartj/ehz305&link_type=DOI) 16. 16.Kuppe C, Ramirez Flores RO, Li Z, Hayat S, Levinson RT, Liao X, et al. Spatial multi-omic map of human myocardial infarction. Nature. 2022 Aug 25;608(7924):766–77. 17. 17.van Blokland IV, Groot HE, Hendriks T, Assa S, van der Harst P. Sex differences in leukocyte profile in ST-elevation myocardial infarction patients. Sci Rep. 2020 Apr 22;10(1):6851. 18. 18.Tigchelaar EF, Zhernakova A, Dekens JAM, Hermes G, Baranska A, Mujagic Z, et al. Cohort profile: LifeLines DEEP, a prospective, general population cohort study in the northern Netherlands: study design and baseline characteristics. Open Access.: 9. 19. 19.Oelen R, de Vries DH, Brugge H, Gordon MG, Vochteloo M, single-cell eQTLGen consortium, et al. Single-cell RNA-sequencing of peripheral blood mononuclear cells reveals widespread, context-specific gene expression regulation upon pathogenic exposure. Nat Commun. 2022 Jun 7;13(1):3267. 20. 20.Zheng GXY, Terry JM, Belgrader P, Ryvkin P, Bent ZW, Wilson R, et al. Massively parallel digital transcriptional profiling of single cells. Nat Commun. 2017 Apr;8(1):14049. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ncomms14049&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28091601&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F03%2F2023.05.02.23289370.atom) 21. 21.Loh PR, Danecek P, Palamara PF, Fuchsberger C, A Reshef Y, K Finucane H, et al. Reference-based phasing using the Haplotype Reference Consortium panel. Nat Genet. 2016 Nov;48(11):1443–8. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3679&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27694958&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F03%2F2023.05.02.23289370.atom) 22. 22.Heaton H, Talman AM, Knights A, Imaz M, Gaffney DJ, Durbin R, et al. Souporcell: robust clustering of single-cell RNA-seq data by genotype without reference genotypes. Nat Methods. 2020 Jun;17(6):615–20. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41592-020-0820-1&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32366989&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F03%2F2023.05.02.23289370.atom) 23. 23.Macosko EZ, Basu A, Satija R, Nemesh J, Shekhar K, Goldman M, et al. Highly Parallel Genome-wide Expression Profiling of Individual Cells Using Nanoliter Droplets. Cell. 2015 May;161(5):1202–14. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cell.2015.05.002&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26000488&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F03%2F2023.05.02.23289370.atom) 24. 24.Butler A, Hoffman P, Smibert P, Papalexi E, Satija R. Integrating single-cell transcriptomic data across different conditions, technologies, and species. Nat Biotechnol. 2018 May;36(5):411–20. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nbt.4096&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29608179&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F03%2F2023.05.02.23289370.atom) 25. 25.Stuart T, Butler A, Hoffman P, Hafemeister C, Papalexi E, Mauck WM, et al. Comprehensive Integration of Single-Cell Data. Cell. 2019 Jun;177(7):1888–1902.e21. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cell.2019.05.031&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F03%2F2023.05.02.23289370.atom) 26. 26.Hafemeister C, Satija R. Normalization and variance stabilization of single-cell RNA-seq data using regularized negative binomial regression. Genome Biol. 2019 Dec;20(1):296. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s13059-019-1874-1&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31870423&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F03%2F2023.05.02.23289370.atom) 27. 27.Ilicic T, Kim JK, Kolodziejczyk AA, Bagger FO, McCarthy DJ, Marioni JC, et al. Classification of low quality cells from single-cell RNA-seq data. Genome Biol. 2016 Dec;17(1):29. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s13059-016-0888-1&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26887813&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F03%2F2023.05.02.23289370.atom) 28. 28.Hao Y, Hao S, Andersen-Nissen E, Mauck WM, Zheng S, Butler A, et al. Integrated analysis of multimodal single-cell data. bioRxiv. 2020 Jan 1;2020.10.12.335331. 29. 29.Farbehi N, Patrick R, Dorison A, Xaymardan M, Janbandhu V, Wystub-Lis K, et al. Single-cell expression profiling reveals dynamic flux of cardiac stromal, vascular and immune cells in health and injury. eLife. 2019 Mar 26;8:e43882. 30. 30.Finak G, McDavid A, Yajima M, Deng J, Gersuk V, Shalek AK, et al. MAST: a flexible statistical framework for assessing transcriptional changes and characterizing heterogeneity in single-cell RNA sequencing data. Genome Biol. 2015 Dec 10;16:278. 31. 31.Chen J, Bardes EE, Aronow BJ, Jegga AG. ToppGene Suite for gene list enrichment analysis and candidate gene prioritization. Nucleic Acids Res. 2009 Jul 1;37(Web Server):W305–11. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gkp427&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19465376&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F03%2F2023.05.02.23289370.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000267889100054&link_type=ISI) 32. 32.Browaeys R, Saelens W, Saeys Y. NicheNet: modeling intercellular communication by linking ligands to target genes. Nat Methods. 2020 Feb;17(2):159–62. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31819264&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F03%2F2023.05.02.23289370.atom) 33. 33.Türei D, Korcsmáros T, Saez-Rodriguez J. OmniPath: guidelines and gateway for literature-curated signaling pathway resources. Nat Methods. 2016 Dec;13(12):966–7. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nmeth.4077&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27898060&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F03%2F2023.05.02.23289370.atom) 34. 34.Aragam KG, Jiang T, Goel A, Kanoni S, Wolford BN, Atri DS, et al. Discovery and systematic characterization of risk variants and genes for coronary artery disease in over a million participants. Nat Genet. 2022 Dec;54(12):1803–15. 35. 35.LifeLines cohort study, BIOS consortium, Zhernakova DV, Le TH, Kurilshikov A, Atanasovska B, et al. Individual variations in cardiovascular-disease-related protein levels are driven by genetics and gut microbiome. Nat Genet. 2018 Nov;50(11):1524–32. 36. 36.Talwar S. Profile of plasma N-terminal proBNP following acute myocardial infarction. Correlation with left ventricular systolic dysfunction. Eur Heart J. 2000 Sep 15;21(18):1514– 21. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1053/euhj.1999.2045&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10973765&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F03%2F2023.05.02.23289370.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000089128000009&link_type=ISI) 37. 37.Sabatine MS, Morrow DA, Higgins LJ, MacGillivray C, Guo W, Bode C, et al. Complementary Roles for Biomarkers of Biomechanical Strain ST2 and N-Terminal Prohormone B-Type Natriuretic Peptide in Patients With ST-Elevation Myocardial Infarction. Circulation. 2008 Apr 15;117(15):1936–44. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTQ6ImNpcmN1bGF0aW9uYWhhIjtzOjU6InJlc2lkIjtzOjExOiIxMTcvMTUvMTkzNiI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIzLzA1LzAzLzIwMjMuMDUuMDIuMjMyODkzNzAuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 38. 38.Braunwald E. Heart Failure. JACC Heart Fail. 2013 Feb;1(1):1–20. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamhmIjtzOjU6InJlc2lkIjtzOjU6IjEvMS8xIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMDUvMDMvMjAyMy4wNS4wMi4yMzI4OTM3MC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 39. 39.Ueland T, Laugsand LE, Vatten LJ, Janszky I, Platou C, Michelsen AE, et al. Extracellular matrix markers and risk of myocardial infarction: The HUNT Study in Norway. Eur J Prev Cardiol. 2017 Jul;24(11):1161–7. 40. 40.Gaul DS, Stein S, Matter CM. Neutrophils in cardiovascular disease. Eur Heart J. 2017 Jun 7;38(22):1702–4. 41. 41.Zheng W, Lai Y, Jin P, Gu W, Zhou Q, Wu X. Association of Circulating IGFBP1 Level with the Severity of Coronary Artery Lesions in Patients with Unstable Angina. Dis Markers. 2017;2017:1–7. 42. 42.Wu X, Zheng W, Jin P, Hu J, Zhou Q. Role of IGFBP1 in the senescence of vascular endothelial cells and severity of aging-related coronary atherosclerosis. Int J Mol Med [Internet]. 2019 Sep 12 [cited 2023 Mar 26]; Available from: [http://www.spandidos-publications.com/10.3892/ijmm.2019.4338](http://www.spandidos-publications.com/10.3892/ijmm.2019.4338) 43. 43.Stenemo M, Nowak C, Byberg L, Sundström J, Giedraitis V, Lind L, et al. Circulating proteins as predictors of incident heart failure in the elderly: Circulating proteins as predictors of incident heart failure. Eur J Heart Fail. 2018 Jan;20(1):55–62. 44. 44.Lagace TA. PCSK9 and LDLR degradation: regulatory mechanisms in circulation and in cells. Curr Opin Lipidol. 2014 Oct;25(5):387–93. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/MOL.0000000000000114&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25110901&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F03%2F2023.05.02.23289370.atom) 45. 45.Segura AM, Frazier OH, Buja LM. Fibrosis and heart failure. Heart Fail Rev. 2014 Mar;19(2):173–85. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s10741-012-9365-4&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23124941&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F03%2F2023.05.02.23289370.atom) 46. 46.Wang J, Wang M, Lu X, Zhang Y, Zeng S, Pan X, et al. IL-6 inhibitors effectively reverse post-infarction cardiac injury and ischemic myocardial remodeling via the TGF-β1/Smad3 signaling pathway. Exp Ther Med. 2022 Jul 18;24(3):576. 47. 47.Rashid I, Maghzal GJ, Chen YC, Cheng D, Talib J, Newington D, et al. Myeloperoxidase is a potential molecular imaging and therapeutic target for the identification and stabilization of high-risk atherosclerotic plaque. Eur Heart J. 2018 Sep 14;39(35):3301–10. 48. 48.Nadel J, Jabbour A, Stocker R. Arterial myeloperoxidase in the detection and treatment of vulnerable atherosclerotic plaque: a new dawn for an old light. Cardiovasc Res. 2023 Mar 17;119(1):112–20. 49. 49.Hartman MHT, Eppinga RN, Vlaar PJJ, Lexis CPH, Lipsic E, Haeck JDE, et al. The contemporary value of peak creatine kinase-MB after ST-segment elevation myocardial infarction above other clinical and angiographic characteristics in predicting infarct size, left ventricular ejection fraction, and mortality: Peak CK-MB predicts infarct size, LVEF and mortality. Clin Cardiol. 2017 May;40(5):322–8. 50. 50.Hunt PJ, Richards AM, Nicholls MG, Yandle TG, Doughty RN, Espiner EA. Immunoreactive amino-terminal pro-brain natriuretic peptide (NT-PROBNP): a new marker of cardiac impairment. Clin Endocrinol (Oxf). 1997 Sep;47(3):287–96. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1046/j.1365-2265.1997.2361058.x&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9373449&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F03%2F2023.05.02.23289370.atom) 51. 51.Yamamoto K, Burnett JC, Jougasaki M, Nishimura RA, Bailey KR, Saito Y, et al. Superiority of Brain Natriuretic Peptide as a Hormonal Marker of Ventricular Systolic and Diastolic Dysfunction and Ventricular Hypertrophy. Hypertension. 1996 Dec;28(6):988–94. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1161/01.HYP.28.6.988&link_type=DOI) 52. 52.Ndrepepa G, Braun S, Mehilli J, von Beckerath N, Nekolla S, Vogt W, et al. N-Terminal Pro-Brain Natriuretic Peptide on Admission in Patients With Acute Myocardial Infarction and Correlation With Scintigraphic Infarct Size, Efficacy of Reperfusion, and Prognosis. Am J Cardiol. 2006 Apr;97(8):1151–6. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.amjcard.2005.11.029&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16616017&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F03%2F2023.05.02.23289370.atom) 53. 53.Qian J, Gao Y, Lai Y, Ye Z, Yao Y, Ding K, et al. Single-Cell RNA Sequencing of Peripheral Blood Mononuclear Cells From Acute Myocardial Infarction. Front Immunol. 2022 Jun 29;13:908815. 54. 54.Garbers C, Heink S, Korn T, Rose-John S. Interleukin-6: designing specific therapeutics for a complex cytokine. Nat Rev Drug Discov. 2018 Jun 1;17(6):395–412. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nrd.2018.45&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29725131&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F05%2F03%2F2023.05.02.23289370.atom) 55. 55.ZEUS - A Research Study to Look at How Ziltivekimab Works Compared to Placebo in People With Cardiovascular Disease, Chronic Kidney Disease and Inflammation (ZEUS) [Internet]. Novo Nordisk A/S; Report No.: NCT05021835. Available from: [https://clinicaltrials.gov/ct2/show/NCT05021835](https://clinicaltrials.gov/ct2/show/NCT05021835) 56. 56.Narayan P, Elefterios T. The interleukin-1 receptor type I promotes the development of aging-associated cardiomyopathy in mice. Cytokine 2022 march 57. 57. Kami Pekayvaz, Corinna Losert, Viktoria Knottenberg*, Irene V. van Blokland, Roy Oelen, Hilde E. Groot, Jan Walter Benjamins, Sophia Brambs, Rainer Kaiser, Luke Eivers, Vivien Polewka, Raphael Escaig, Markus Joppich, Aleksandar Janjic, Oliver Popp, Tobias Petzold, Ralf Zimmer, Wolfgang Enard, Kathrin Saar, Philipp Mertins, Norbert Huebner, Pim van der Harst, Lude H. Franke, Monique G. P. van der Wijst, Steffen Massberg, Matthias Heinig, Leo Nicolai and Konstantin Stark. Multi-Omic Factor Analysis uncovers immunological signatures with pathophysiologic and clinical implications in coronary syndromes. bioRxiv 2023