Collapsing the list of myocardial infarction-related differential expressed genes into a diagnostic signature
=============================================================================================================

* German Osmak
* Natalia Baulina
* Philipp Koshkin
* Olga Favorova

## Abstract

Myocardial infarction (MI) is one of the most severe manifestations of coronary artery disease (CAD) and the leading cause of death from non-infectious diseases worldwide. It is known that the central component of CAD pathogenesis is a chronic vascular inflammation. However, the mechanisms underlying the changes that occur in T, B and NK lymphocytes, monocytes and other immune cells during CAD and MI are still poorly understood. One of those pathogenic mechanisms might be the dysregulation of intracellular signaling pathways in the immune cells.

In the present study we performed a transcriptome profiling in peripheral blood mononuclear cells of MI patients and controls. The machine learning algorithm was then used to search for MI-associated signatures, that could reflect the dysregulation of intracellular signaling pathways.

The genes *ADAP2, KLRC1, MIR21, PDGFD* and *CD14* were identified as the most important signatures for the classification model with L1-norm penalty function. The classifier output quality was equal to 0.911 by Receiver Operating Characteristic metric on test data. These results were validated on two independent open GEO datasets. Identified MI-associated signatures can be further assisted in MI diagnosis and/or prognosis.

Thus, our study presents a pipeline for collapsing the list of differential expressed genes, identified by high-throughput techniques, in order to define disease-associated diagnostic signatures.

Keywords
*   transcriptomics
*   machine learning
*   myocardial infarction
*   miRNA
*   transcriptional signatures

## Introduction

Myocardial infarction (MI) is one of the most severe manifestations of coronary artery disease (CAD) and the leading cause of death from non-infectious diseases worldwide [1]. In most cases, MI occurs as a serious complication of atherosclerosis – a complex disease, the etiology of which is still not fully elucidated [2]. Recent studies have shown that the central component of atherosclerosis pathogenesis is a chronic vascular inflammation, resulting in endothelial dysfunction and, consequently, in an increased probability of hemodynamic abnormalities, including through the thrombosis [3]. Such a vascular lesion is emerged with the leading involvement of intimal cells (fibroblasts, endothelial and smooth muscle cells) and peripheral blood mononuclear cells (PBMC) [4]. However, the pathogenic mechanisms underlying the changes that occur in PBMC (T, B and NK lymphocytes and monocytes) during atherosclerosis and MI are still poorly understood. One of those mechanisms might be the dysregulation of intracellular signaling pathways in immune cells [5].

One way to establish MI transcriptional signatures, which include both dysregulated individual genes and signaling pathways, containing dysregulated genes, is to simultaneously study of transcriptional profiles of protein-coding genes and genes for regulatory non-coding RNA. Among non-coding RNAs, miRNAs are of particular interest in the context of robustness of biological processes, since they regulate key elements of extensive segments of signaling pathways’ networks [6–8]. To date, consistency of expressional changes of miRNAs and their target genes has been investigated in macrophages of pigs and rats with experimental MI [9] and in the whole blood of MI patients [10].

In the present study, we performed a transcriptome profiling in PBMC of MI patients and healthy individuals and revealed MI-associated signatures, consisting of individual protein-coding genes or functional patterns of genes, such as miRNA with its co-expressed target genes or combination of co-expressed genes, attributed to a definite signaling pathway.

## Materials and methods

### Pipeline

The pipeline of the study design is illustrated in Figure 1. RNA Microarray analysis was used to identify genes that were significantly (*p*<0.05) associated with MI (differentially expressed genes, DEGs). Those DEGs that have passed threshold for multiple comparisons were considered MI transcriptional signatures. Functional patterns of co-expressed DEGs were also considered MI transcriptional signatures. Such functional patterns included i) differentially expressed miRNA and its co-expressed target mRNA(s) and ii) DEGs attributed to a Reactome gene set. In the latter case, the Reactome gene sets were considered the most informative if they i) account for more than 10% of all co-expressed DEGs and/or ii) include DEGs passed multiple comparisons correction.

![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/02/03/2020.01.29.20019554/F1.medium.gif)

[Figure 1.](http://medrxiv.org/content/early/2020/02/03/2020.01.29.20019554/F1)

Figure 1. 
A schematic pipeline of the study for MI transcriptional signatures’ identification. DEGs – differentially expressed genes. MI – patients with myocardial infarction. CTRLs – individuals in the control group, CV - cross-validation.

The validation of identified MI transcriptional signatures was performed on two open data sets; DEGs which were not validated on at least one of these sets were excluded from further consideration. The DEGs within validated MI transcriptional signatures were used to construct binary classifiers. Given the high quality of classification and stability of the detected composite transcriptional biomarker, a logistic regression with the L1-norm penalty function was used to select the most significant DEGs on test dataset.

### Patients and controls

Six patients (all men, mean age 51.3 ± 5.9 years) with first ST-segment elevation MI were enrolled in this study. All patients were diagnosed at the National Medical Scientific Center for Cardiology (Moscow, Russia) based on symptoms of myocardial ischemia, increase of high-sensitivity cardiac troponin I (hs-cTn-I) and/or emergence of new or presumed new ST-segment elevation, new left bundle branch block or development of pathological Q waves in accordance with Third Universal Definition of MI [11]. Hs-cTn-I was measured during the initial patient assessment (from 1 to 18 h after the onset of disease symptoms). All patients underwent coronary angiography on admission and were treated according to contemporary guidelines. The characteristics of MI patients are presented in Table S1. A total of 6 CTRLs (all men, mean age 51.0 ± 7.1 years) with normal electrocardiogram, no history of CVD and diabetes mellitus were included in the study; CTRLs characteristics by smoking status and body mass index were compatible to MI patients. All participants lived in European Russia. The ethical approval was obtained from the local Ethics Committee, and written informed consent had been received from each person in accordance with the Declaration of Helsinki.

### Peripheral blood mononuclear cells collection and RNA extraction

Blood samples were collected in the morning from MI patients (24-36 h after the disease onset) and CTRLs. PBMC were isolated using Ficoll-Hypaque density gradient method (Sigma-Aldrich, St. Louis, MO, USA) within 3 h of sampling. Total RNA including small RNA was extracted using miRNeasy Mini Kit (Qiagen, Hildren, Germany) following the manufacturer’s instructions. The RNA quantity was measured using the NanoDropTM spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA); the RNA integrity was assessed by QIAxcel Advanced System (Qiagen, Hilden, Germany). Samples with RNA integrity number (RIN) value above eight were included in subsequent experiments.

### RNA Microarray analysis

The transcriptome analysis was performed using GeneChip Human Transcriptome Array 2.0, which provides the ability to analyze the expression of 44,699 protein-coding genes and 22,829 non-protein coding genes, including 1346 miRNA genes (ThermoFisher Scientific, Santa Clara, CA, USA). Briefly, total RNA (500 ng) of samples were each proceeded to poly(A)tailing and biotin ligation reactions using FlashTag Biotin HSR RNA Labeling Kit (ThermoFisher Scientific, Santa Clara, CA, USA). The biotin-labeled RNA samples were hybridized on GeneChip Human Transcriptome Array 2.0 using manufacturer’s instructions and scanned on the GeneChip Scanner 7G System. Computational analysis of the microarray data files was performed using R programming language version 3.5.1. Data processing was carried out based on the affy package written in R [12]. A biomaRt package was used to annotate the obtained data [13]. Probes demonstrating evidence for cross-hybridization, i.e. transcript sequences annotated to more than two coding genes were excluded from this study. If transcripts belong to the same gene ID, a transcript with the most detectable expression level was selected. To detect differentially expressed genes, calculate the levels of statistical significance and adjust them for multiple comparisons by Benjamini-Hochberg procedure (*p* and *padj*, respectively) the standard limma package protocol was used [14]. All expression data are deposited in the Gene Expression Omnibus international public repository under accession identification as GSE141512 [15].

### Bioinformatic analysis

MirTarBase was used to select experimentally validated target genes for miRNAs [16]. Gene set enrichment analysis (over-representation analysis) was performed using Tools of Reactome Database [17].

To construct and analyze the gene-gene interaction networks, NetworkX 2.0 package for Python was used [18]. STRING database [19] was used to find protein-protein interactions.

### Statistics analysis and Machine learning

Statistical analysis was performed using R programming language version 3.5.1. Null hypotheses were rejecting if p<0.05. To study the dependence/correlation of two continuous random variables, the Spearman’s Rank Correlation test was used. The logistic regression classifier was trained using the tools of scikit-learn v0.20.3 for Python [20]. To reduce the possible classification model overfitting the l2-norm regularization and 10-fold cross validation were used. The selection of the optimal regularization coefficient was performed by grid search using the GridSearchCV () function. The quality of the classification model was estimated by the areas under receiver operating characteristic curve (ROC-AUC). The final assessment of the quality of the classification model was carried out on the test dataset that was not used for training. For training and testing the model z-scaling of continuous features was performed at the preprocessing data stage.

### Validation analysis

The Gene Expression Omnibus database (GEO, [http://www.ncbi.nlm.nih.gov/geo](http://www.ncbi.nlm.nih.gov/geo)) was used in order to validate the obtained results. Two open data sets - GSE59867 and GSE62646 with gene expression profiles in PBMC of MI patients and healthy individuals without a history of CVD were investigated; they were obtained on GeneChip Human Gene 1.0 ST Array [transcript (gene) version]. The dataset GSE59867 included expression data of 111 MI cases and 48 CTRLs; the dataset GSE62646 - of 28 MI patients and 14 CTRLs (mixed-gender sets).

## Results

### Array-based transcriptome profiling

Transcriptome profiling in PBMC of six MI patients and six gender- and age-matched control individuals (CTRLs) was performed using GeneChip Human Transcriptome Array 2.0 (Figure 2). As a result, a total of 84 differentially expressed genes (DEGs) were identified (−0.5<Log2FC>0.5, *p*<0.05) (Table S2), from which 48 protein-coding genes and 2 miRNA genes (*MIR21* and *MIR223*) were upregulated, while 34 protein-coding genes were downregulated in MI patients. Among all observed DEGs *KLRB1* and *ADAP2* passed the threshold for multiple comparisons correction (Log2FC=-0.64, *p*.*adj*=0.0454 and Log2FC=0.64, *p*.*adj*=0.0495, respectively); both these genes were further considered as MI transcriptional signatures.

![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/02/03/2020.01.29.20019554/F2.medium.gif)

[Figure 2.](http://medrxiv.org/content/early/2020/02/03/2020.01.29.20019554/F2)

Figure 2. 
Volcano plot of gene expression changes in PBMC of MI patients compared to CTRLs. Blue dot indicates downregulated gene (Log2FC<-0.5); red dot indicates upregulated gene (Log2FC>0.5), which passed threshold for multiple comparisons (*p*.*adj*<0.05); Among differentially expressed genes (DEGs) *MIR21* and its target genes are marked in orange, *MIR223* and its target gene *−* in purple (−0.5<Log2FC>0.5, *p*<0.05).

### The search for MI transcriptional signatures: miRNA and its target mRNA(s)

Among identified DEGs (*p*<0.05) presented in Figure 2, *BCL6, CCR1, PDGFD, SGK1*, and *TGFBR3* genes were found to be targets of miR-21, while *MAFB* – target of miR-223 based on MirTarBase database. As assessed by Spearman’s correlation analysis the expression levels of *BCL6, CCR1*, and *SGK1* were positively correlated (*p*<0.001, Ro>0.9) and *PDGFD* and *TGFBR3 -* negatively correlated (*p*<0.05; Ro<-0.6) with *MIR21* expression level in MI patients and CTRLs (Figure S1). A positive correlation between the expression levels of *MAFB* and *MIR223* (*p*<0.1, Ro=0.5) was also observed (Figure S2). Thus, *MIR21* and *MIR223* genes, together with their functionally associated co-expressed target genes, were considered as two MI transcriptional signatures.

### The search for MI transcriptional signatures: Reactome gene sets

The enrichment analysis was undertaken in order to search for the functional patterns which included DEGs attributed to a Reactome gene set (Table 1). Nine Reactome gene sets were significantly overrepresented (FDR<0.05) among the 48 upregulated protein-coding genes (see above). The first three sets included each more than 10% of upregulated genes: “Immune system” *−* 22 DEGs from 2663 genes presented in the set (FDR=0.023), “Neutrophil degranulation” *−* 13 DEGs from 480 genes (FDR=0.0035) and “Cytokine Signaling in Immune system” *−* 9 DEGs from 1055 genes (FDR=0.015). “Immune system” gene set is at the highest level of the Reactome hierarchy and includes “Neutrophil degranulation” and “Cytokine Signaling in Immune system” pathways that are separately characterized by more significant overrepresentation of DEGs. So that, the DEGs from these two pathways were chosen to further analysis in the context of potential MI transcriptional signatures. Notably, “Cytokine Signaling in Immune system” pathway involves *BCL6* and *CCR1* - the target genes of miR-21, which were already included in one of the identified MI transcriptional signatures.

View this table:
[Table 1.](http://medrxiv.org/content/early/2020/02/03/2020.01.29.20019554/T1)

Table 1. 
Reactome gene sets significantly overrepresented among the differentially expressed genes in PBMC from MI patients when compared to controls

As can be seen from Table 1, two Reactome gene sets were significantly overrepresented among the 34 downregulated in MI genes: “Immunoregulatory interactions between a Lymphoid and a non-Lymphoid cell” (6 DEGs from 297 genes, FDR=0.021) and “DAP12 signaling” (4 DEGs from 29 genes, FDR=0.021). Each of these sets includes more than 10% of the downregulated genes and is involved in signal transduction in lymphoid cells, namely in natural killers (NK). The “DAP12 signaling” pathway was overrepresented exclusively among the genes which encode the killer cell lectin-like receptors (KLR) expressed in NK cells, two of these DEGs (*KLRD1* and *KLRC1*) were as well observed in “Immunoregulatory interactions between a Lymphoid and a non-Lymphoid cell” pathway. The last mentioned pathway was also overrepresented by *KIR2DL1* and *KIR2DL3* genes, encoding a killer cell immunoglobulin-like receptors, the transmembrane glycoproteins expressed by NK and T cells’ subsets. Notably, *KLRB1* gene, defined previously individually as MI transcriptional signature (Log2FC=-0.64, *p*.*adj*=0.0454) was included only in “Immunoregulatory interactions between a Lymphoid and a non-Lymphoid cell” pathway. The DEGs from this pathway were chosen for further analysis in the context of potential MI transcriptional signatures.

The search for interacting genes/proteins among the DEGs from selected Reactome sets “Neutrophil degranulation”, “Cytokine Signaling in Immune system” and “Immunoregulatory interactions between a Lymphoid and a non-Lymphoid cell” (lines 2, 3 and 10 in Table 1) was performed using String database. Almost all the DEGs from the “Neutrophil degranulation” set (with the exception of the *PADI2, GRN* and *PYGL* genes) were found to interact among themselves (Figure 3A). The expression levels of these 10 interacting genes were significantly positively correlated between each other (0.93>Ro>0.51, *p*<0.05) (Figure S3). Thus, we considered the pattern of interacting genes from this pathway, namely *BST1, C3AR1, CD14, CLEC4D, CR1, FPR1, FPR2, S100A12, SLC11A1* and *TLR2* as potential MI transcriptional signature.

![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/02/03/2020.01.29.20019554/F3.medium.gif)

[Figure 3.](http://medrxiv.org/content/early/2020/02/03/2020.01.29.20019554/F3)

Figure 3. 
Network analysis of the Reactome gene sets “Neutrophil degranulation” (A), “Cytokine Signaling in Immune system” (B) and “Immunoregulatory interactions between a Lymphoid and a non-Lymphoid cell” (C). The edges indicate molecular interactions between nodes based on String database.

In the “Cytokine Signaling in Immune system” gene set six interacting genes (*CCR1, FCGR1A, FCGR1B, FPR1, HLA-DQB1* and *S100A12*) were found (Figure 3B), and the expression levels of these genes, with the exception of the *HLA-DQB1*, were positively correlated with each other (0.99>Ro>0.62, *p*<0.05) (Figure S4). As previously mentioned, *CCR1* is the target gene of miR-21, and has already been included to *MIR21*-containing MI transcriptional signature. According to correlation analysis, *MIR21* expression level positively correlates not only with *CCR1*, but also with *FCGR1A, FCGR1B, FPR1* and *S100A12* expression levels (0.93>Ro>0.71, *p*<0.01) (Figure S4). Thus, we considered these genes as components that extend the *MIR21*-containing MI transcriptional signature.

In Reactome pathway “Immunoregulatory interactions between a Lymphoid and a non-Lymphoid cell” different interacting gene pairs were found between all the genes from this set: *KLRB1, KLRC1, KLRD1, KLRF1, KIR2DL1, KIR2DL3* (Figure 2C) and they were predominantly characterized by significantly positive correlation between their expression levels (0.88>Ro>0.55, *p*<0.05) (Figure S5). Noteworthy, *KLRB1* gene whose differential expression passed correction for multiple comparisons and could be considered as the MI transcriptional signature, interacts with *KLRD1* and *KLRF1* from this gene set. Therefore, we included genes *KLRC1, KLRD1, KLRF1, KIR2DL1* and *KIR2DL3* in the *KLRB1*-containing MI transcriptional signature.

Overall, the conducted analysis allowed us identifying the following five MI transcriptional signatures containing all in all 29 DEGs: {*ADAP*}, {*KLRB1* + *KLRC1, KLRD1, KLRF1, KIR2DL1*, and *KIR2DL3*}, {*MIR21* + *BCL6, CCR1, PDGFD, SGK1, TGFBR3, FCGR1A, FCGR1B, FPR1*, and *S100A12*}, {*MIR223* + *MAFB*} and {*BST1, C3AR1, CD14, CLEC4D, CR1, FPR1, FPR2, S100A12, SLC11A1*, and *TLR2*}.

### The validation analysis of differential expression of genes in identified MI transcriptional signatures using GEO datasets

To confirm the differential expression of genes in identified MI transcriptional signatures we used open datasets GSE62646 and GSE59867 from GEO database, in which gene expression profiles in PBMC of MI patients and healthy individuals without a history of cardiovascular diseases (CVD) were investigated using GeneChip Human Gene 1.0 ST Array. The genes *KIR2DL1, KIR2DL3, FCGR1A* and *FCGR1B* that were according to our results included in MI transcriptional signatures were not represented on this array and were therefore excluded from the corresponding MI transcriptional signatures on further consideration. Thus, in a further analysis, 25 genes were considered.

Of the five MI transcriptional signatures we identified, the differential expression of all genes included in the *ADAP2-, KLRB1*-, and *MIR223*-containing MI transcriptional signatures was validated in both open datasets GSE62646 and GSE59867 (Table 2). The differential expression of all genes included in the *MIR21*-containing MI transcriptional signature and MI transcriptional signature from “Neutrophil degranulation” Reactome gene set was validated in GSE59867 dataset. The differential expression of a number of these genes was validated also on GSE62646 dataset with the exception of *FPR1* and *SGK1* from the *MIR21*-containing signature and *BST1, CLEC4D, FPR1, FPR2* and *TLR2* from the “Neutrophil degranulation” signature; these genes were excluded from the corresponding MI transcriptional signatures on further consideration.

View this table:
[Table 2.](http://medrxiv.org/content/early/2020/02/03/2020.01.29.20019554/T2)

Table 2. 
The expression of genes from identified MI transcriptional signatures based on our data and data obtained from GSE62646 and GSE59867 GEO datasets. Grey colour indicates p-value > 0.05.

Thus, after the validation analyses MI transcriptional signatures look as follows:{*ADAP2*}, {*KLRB1* + *KLRC1, KLRD1, KLRF1*}, {*MIR21* + *BCL6, CCR1, PDGFD, TGFBR3, S100A12*}, {*MIR223* + *MAFB*} and {*C3AR1, CD14, CR1, S100A12, SLC11A1*}.

### The diagnostic value of the identified MI transcriptional signatures

The design of our study does not allow to assess the causality between MI and validated transcriptional signatures, which does not exclude the possibility of considering them as diagnostic biomarkers. Their diagnostic value can be assessed by the quality of the classification of MI patients from healthy controls. To search for such an optimal classifier, a L2 regularized logistic regression model was trained on the GSE59867 dataset (Figure 4).

![Figure 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/02/03/2020.01.29.20019554/F4.medium.gif)

[Figure 4.](http://medrxiv.org/content/early/2020/02/03/2020.01.29.20019554/F4)

Figure 4. 
Quality and robustness of the classification model with a L2-norm penalty function based on the considered MI transcriptional signatures: {*ADAP2*}, {*KLRB1* + *KLRC1, KLRD1, KLRF1*}, {*MIR21* + *BCL6, CCR1, PDGFD, TGFBR3, S100A12*}, {*MIR223* + *MAFB*} and {*C3AR1, CD14, CR1, S100A12, SLC11A1*}. (A) Areas Under receiver operating characteristic Curve (ROC-AUC) for the training (GSE59867) and test (GSE62646) datasets. (B) Time-depended (starting from MI onset) ROC-AUC metrics of the classification model.

Figure 4A shows that MI patients at the time of admission to hospital could be classified from healthy individuals based on the selected MI transcriptional signatures, and the quality of the classification model on the test dataset (AUC=0.926) is slightly higher than on the training dataset (AUC=0.910), illustrating the stability of the model and the lack of its overfitting. While analyzing the available data from GSE62646 dataset on the levels of gene expression during the six-month follow-up after MI (Figure 4B), we observed that the classification model remains effective within 6 days after MI onset (AUC=0.773, blue line) but 6 month after MI onset the effectiveness of this model considerably decreases (AUC=0.594, red line).

For the feature selection and to reduce the number of DEGs included in the classification model we used a logistic regression with the L1-norm penalty function (Figure 5). As a result, *ADAP2, KLRC1, MIR21, PDGFD* and *CD14* genes were selected for the classification model as the most important DEGs (Figure 5A). ROC-curves constructed for these genes are demonstrated in Figures 5B and C. The comparison of ROC-curves from Figures 4 and 5 demonstrates that the quality of the classification on test dataset slightly decreased from 0.926 (dark green curve on Figure 4B) to 0.911 after applying L1 regularization (dark green curve on Figure 5B). While analyzing the changes in the quality of the classification model based on the levels of gene expression over time after MI onset, ROC-AUC values also slightly changed after applying L1 regularization (see Figures 4B and 5C). Thus, five DEGs are sufficient for the classification; among these genes DEGs from *MIR223*-containing MI transcriptional signature were not presented.

![Figure 5.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/02/03/2020.01.29.20019554/F5.medium.gif)

[Figure 5.](http://medrxiv.org/content/early/2020/02/03/2020.01.29.20019554/F5)

Figure 5. 
Quality and robustness of the classification model with a L1-norm penalty function based on the considered MI transcriptional signatures: {*ADAP2*}, {*KLRB1* + *KLRC1, KLRD1, KLRF1*}, {*MIR21* + *BCL6, CCR1, PDGFD, TGFBR3, S100A12*}, {*MIR223* + *MAFB*}and {*C3AR1, CD14, CR1, S100A12, SLC11A1*}. (A) Coefficients of the classification model; the most important upregulated genes *ADAP2, MIR21* and *CD14* are marked in red, downregulated genes *KLRC1* and *PDGFD* - in blue colour. (B) ROC-AUC metrics of the L1-regularized classification model consisted of *ADAP2, MIR21* and *CD14* genes. ROC-AUC were constructed using the training (GSE59867) and test (GSE62646) datasets. (C) Time-depended (starting from MI onset) ROC-AUC metrics of the L1-regularized classification model based on test dataset.

## Discussion

In the study we searched for MI transcriptional signatures (individual dysregulated genes or functional patterns of dysregulated genes) that could be potentially used in MI diagnosis. We compared the transcriptome profiles in PBMC of patients with first MI and healthy individuals using GeneChip Human Transcriptome Array 2.0 and identified five MI transcriptional signatures containing from 1 to 6 DEGs: {*ADAP2*}, {*KLRB1 + KLRC1, KLRD1, KLRF1*}, {*MIR21* + *BCL6, CCR1, PDGFD, TGFBR3, S100A12*}, *{MIR223* + *MAFB}* and {*C3AR1, CD14, CR1, S100A12, SLC11A1*}.

In order to select the most important for classification DEGs, further research steps included: validation on open datasets GSE62646 [21] and GSE59867 [22], construction of L2 regularized logistic regression model for estimation the diagnostic value of the MI transcriptional signatures and feature selection using L1-norm penalty function. This approach allowed to exclude from the classification model MI transcriptional signature (*MIR223* + *MAFB*) as insufficiently effective and to reduce the number of DEGs from other signatures to *ADAP2, KLRC1, MIR21, PDGFD* and *CD14*. According to the ROC-AUC analysis the obtained classification model, including 5 genes, is enable classifying MI patients and healthy controls with a quality of 0.911 while the quality of initial classification model, including 18 genes, was equal to 0.926. Thus, a decrease in the number of genes did not significantly affect the quality of the model. A comparable decline in the quality of both classification models over time from MI onset was shown; this decline occurs rather slowly, for days and weeks.

Consider consistently the characteristics of genes-classifiers. The gene *ADAP2* encodes ArfGAP With Dual PH Domains 2 protein and was designated in our study as individual MI transcriptional signature; no data on the involvement of this gene in the development of CVD and/or its complications were found. However, the product of this gene was shown to be involved in heart development, and its dysfunction presumably is associated with cardiovascular malformations in NF1 microdeletion syndrome [23].

The gene *KLRC1* refers to a MI transcriptional signature containing the genes of killer cell lectin-like receptors (KLR) that encode a family of transmembrane proteins, characterized by a type II membrane orientation and the presence of a C-type lectin domain; they are predominantly expressed in NK cells. The association of some genes from this signature (*KLRD1* and *KLRC1*) with MI or its complications was previously shown by Maciejak et al and Kiliszek et al, whose data were used for validation analysis in our study [21,22]. We have shown that the expression of the genes *KLRB1, KLRC1, KLRD1, KLRF1* is consistently decreased in MI, that is in a good accordance with the study by Yan et al [24], where a loss of NK cell activity was found in patients with acute MI, in particular, due to a decrease in *KLRB1* expression.

The *MIR21* gene and target genes of miR-21 were included in one MI transcriptional signature, composed mainly of genes from “Cytokine Signaling in Immune system” pathway. Thus, in addition to cytokine signaling pathways, which role in MI development was previously described [25], we have identified and validated on independent GEO datasets the influence of miR-21 through the regulation of this pathway in PBMC during MI. The functional role of miR-21 in cardiac tissue has been studied for a long time, and by now a large amount of data has been accumulated on this subject [26], while in PBMC its role remains unclear. In one of the studies the negative correlation of miR-21 expression level in MI with the levels of IL-1β, IL-6, and TNF-α cytokines was shown due to regulatory effect of this miRNA on the expression of *KBTBD7*; this gene encodes a member of BTB-kelch proteins, kelch repeat and BTB (POZ) domain containing 7, which promotes inflammatory responses in macrophages [27]. In turn, in our study, the association of miR-21 and its target genes *PDGFD, TGFBR3, CCR1* and *BCL6* expression levels with MI was demonstrated, from which *PDGFD* gene encoded platelet derived growth factor D was found to be the most important based on the results of L1 regularization. The genes of the PDGF family and their involvement in the pathogenesis of various diseases are well studied; in particular, PDGFD is known to be involved in the fibrosis and neovascularization of the cardiac tissue [28].

The gene *CD14* encodes a receptor on the surface of myeloid cells, which participates in CD14/TLR4/MD2 signaling pathway involved in the recognition of lipopolysaccharides [29]. This gene was identified in our study as a component of “Neutrophil degranulation” pathway. The neutrophils are known to be actively involved in the development and elimination of MI consequences [30]. Furthermore, polymorphic variants in *CD14* gene were found to be associated with MI [31].

Further investigations are implicitly needed to clarify the causality between MI and the identified MI-associated signatures.

The data on differential expression of a number of genes in PBMC of MI patients obtained in our study were validated on two independent datasets that indicates their value. The identified DEGs could be suitable for the prediction of the first MI before the appearance of the disease symptoms, as it was previously described for some miRNAs [32]. Further investigations are implicitly needed to clarify the functional role of the identified MI-associated genes in the development of this disease.

## Conclusions

The present study implements the pipeline designed to the collapsing the list of differential expressed in MI genes into a diagnostic signature; the obtained classification model is enable classifying MI patients and healthy controls with a quality of 0.911 on a test data. This pipeline could be useful in high-throughput data analysis for the searching of diagnostic signature of other diseases.

## Data Availability

https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE141512

## Authors’ contributions

Conceptualization, G.O., N.B. and O.F.; Data curation, G.O. and P.K.; Formal analysis, G.O.; Investigation, G.O., N.B. and P.K.; Methodology, N.B.; Project administration, O.F.; Validation, G.O.; Visualization, G.O.; Writing – original draft, G.O. and O.F.; Writing – review & editing, G.O., N.B. and O.F‥

## Competing interests

The authors declare that they have no conflict of interest relating to the conduct of this study or the publication of this manuscript.

## Funding

This study was partially supported by the Grant No.16-14-10251 from the Russian Science Foundation and the grant №075-15-2019-1789 from the Ministry of Science and Higher Education of the Russian Federation allocated to the Center for Precision Genome Editing and Genetic Technologies for Biomedicine.

## Disclosures

None.

![Figure6](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/02/03/2020.01.29.20019554/F6.medium.gif)

[Figure6](http://medrxiv.org/content/early/2020/02/03/2020.01.29.20019554/F6)

## Acknowledgments

We thank Dr. Roman Shakhnovich, Dr. Nino Kulava and PhD Natalia Matveeva from National Medical Research Center for Cardiology, Russia, for their assistance with MI patients and controls collection. We appreciate PhD Alexander Favorov for valuable advices while producing the manuscript.

*   Received January 29, 2020.
*   Revision received January 29, 2020.
*   Accepted February 3, 2020.


*   © 2020, 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.

## Reference

1.  [1]. R. Roberts, Genetics of coronary artery disease, Circ. Res. 114 (2014) 1890–1903. [https://doi.org/10.1161/CIRCRESAHA.114.302692](https://doi.org/10.1161/CIRCRESAHA.114.302692).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTA6ImNpcmNyZXNhaGEiO3M6NToicmVzaWQiO3M6MTE6IjExNC8xMi8xODkwIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDIvMDMvMjAyMC4wMS4yOS4yMDAxOTU1NC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

2.  [2]. M.-Y. Wu,  C.-J. Li,  M.-F. Hou,  P.-Y. Chu, New Insights into the role of inflammation in the pathogenesis of atherosclerosis, Int. J. Mol. Sci. 18 (2017) pii: E2034.
    
    

3.  [3]. D.P. Ramji,  T.S. Davies, Cytokines in atherosclerosis: Key players in all stages of disease and promising therapeutic targets, Cytokine Growth Factor Rev. 26 (2015) 673–685. [https://doi.org/10.1016/j.cytogfr.2015.04.003](https://doi.org/10.1016/j.cytogfr.2015.04.003).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cytogfr.2015.04.003&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26005197&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F02%2F03%2F2020.01.29.20019554.atom) 

4.  [4]. R. Ross, Atherosclerosis — An Inflammatory Disease, N. Engl. J. Med. 340 (1999) 115–126. [https://doi.org/10.1056/NEJM199901143400207](https://doi.org/10.1056/NEJM199901143400207).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJM199901143400207&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9887164&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F02%2F03%2F2020.01.29.20019554.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000078016200007&link_type=ISI) 

5.  [5]. A.M. Zawada,  K.S. Rogacev,  S. Müller,  B. Rotter,  P. Winter,  D. Fliser,  G.H. Heine, Massive analysis of cDNA Ends (MACE) and miRNA expression profiling identifies proatherogenic pathways in chronic kidney disease, Epigenetics. 9 (2014) 161–172.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.4161/epi.26931&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24184689&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F02%2F03%2F2020.01.29.20019554.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000331374700018&link_type=ISI) 

6.  [6]. N. Peláez,  R.W. Carthew, Biological robustness and the role of microRNAs: a network perspective, in: Curr. Top. Dev. Biol., Elsevier, 2012: pp. 237–255.
    
    

7.  [7]. A.-L. Barabási,  N. Gulbahce,  J. Loscalzo, Network Medicine: A Network-based Approach to Human Disease, Nat. Rev. Genet. 12 (2011) 56–68. [https://doi.org/10.1038/nrg2918](https://doi.org/10.1038/nrg2918).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nrg2918&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21164525&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F02%2F03%2F2020.01.29.20019554.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000285410500011&link_type=ISI) 

8.  [8]. X. Lai,  O. Wolkenhauer,  J. Vera, Understanding microRNA-mediated gene regulatory networks through mathematical modelling, Nucleic Acids Res. 44 (2016) 6019–6035.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gkw550&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27317695&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F02%2F03%2F2020.01.29.20019554.atom) 

9.  [9]. G. de Couto,  R. Gallet,  L. Cambier,  E. Jaghatspanyan,  N. Makkar,  J.F. Dawkins,  B.P. Berman,  E. Marbán, Exosomal microRNA transfer into macrophages mediates cellular postconditioning, Circulation. 136 (2017) 200–214.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTQ6ImNpcmN1bGF0aW9uYWhhIjtzOjU6InJlc2lkIjtzOjk6IjEzNi8yLzIwMCI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzAyLzAzLzIwMjAuMDEuMjkuMjAwMTk1NTQuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

10. [10]. M. Cheng,  S. An,  J. Li, Identifying key genes associated with acute myocardial infarction, Medicine (Baltimore). 96 (2017).
    
    

11. [11]. K. Thygesen,  J.S. Alpert,  A.S. Jaffe,  M.L. Simoons,  B.R. Chaitman,  H.D. White, Third universal definition of myocardial infarction, Circulation. 126 (2012) 2020–2035.
    
    [FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiRlVMTCI7czoxMToiam91cm5hbENvZGUiO3M6MTQ6ImNpcmN1bGF0aW9uYWhhIjtzOjU6InJlc2lkIjtzOjExOiIxMjYvMTYvMjAyMCI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzAyLzAzLzIwMjAuMDEuMjkuMjAwMTk1NTQuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

12. [12]. L. Gautier,  L. Cope,  B.M. Bolstad,  R.A. Irizarry, affy—analysis of Affymetrix GeneChip data at the probe level, Bioinformatics. 20 (2004) 307–315.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btg405&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=14960456&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F02%2F03%2F2020.01.29.20019554.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000188990900002&link_type=ISI) 

13. [13]. S. Durinck,  Y. Moreau,  A. Kasprzyk,  S. Davis,  B. De Moor,  A. Brazma,  W. Huber, BioMart and Bioconductor: a powerful link between biological databases and microarray data analysis, Bioinformatics. 21 (2005) 3439–3440.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/bti525&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16082012&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F02%2F03%2F2020.01.29.20019554.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000231360600020&link_type=ISI) 

14. [14]. G.K. Smyth, Limma: linear models for microarray data, in: Bioinforma. Comput. Biol. Solut. Using R Bioconductor, Springer, 2005: pp. 397–420.
    
    

15. [15].GEO DataSet, (n.d.). [https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE141512](https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE141512).
    
    

16. [16]. C.-H. Chou,  S. Shrestha,  C.-D. Yang,  N.-W. Chang,  Y.-L. Lin,  K.-W. Liao,  W.-C. Huang,  T.-H. Sun,  S.-J. Tu,  W.-H. Lee, miRTarBase update 2018: a resource for experimentally validated microRNA-target interactions, Nucleic Acids Res. 46 (2017) D296–D302.
    
    

17. [17]. A. Fabregat,  S. Jupe,  L. Matthews,  K. Sidiropoulos,  M. Gillespie,  P. Garapati,  R. Haw,  B. Jassal,  F. Korninger,  B. May,  M. Milacic,  C.D. Roca,  K. Rothfels,  C. Sevilla,  V. Shamovsky,  S. Shorser,  T. Varusai,  G. Viteri,  J. Weiser,  G. Wu,  L. Stein,  H. Hermjakob,  P. D’Eustachio, The Reactome Pathway Knowledgebase, Nucleic Acids Res. 46 (2018) D649–D655. [https://doi.org/10.1093/nar/gkx1132](https://doi.org/10.1093/nar/gkx1132).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gkx1132&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29145629&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F02%2F03%2F2020.01.29.20019554.atom) 

18. [18]. A. Hagberg,  P. Swart,  D. S Chult, Exploring network structure, dynamics, and function using NetworkX, in: Pasadena. USA., 2008: p. 11.
    
    

19. [19]. D. Szklarczyk,  J.H. Morris,  H. Cook,  M. Kuhn,  S. Wyder,  M. Simonovic,  A. Santos,  N.T. Doncheva,  A. Roth,  P. Bork,  L.J. Jensen,  C. von Mering, The STRING database in 2017: quality-controlled protein–protein association networks, made broadly accessible, Nucleic Acids Res. 45 (2017) D362–D368. [https://doi.org/10.1093/nar/gkw937](https://doi.org/10.1093/nar/gkw937).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gkw937&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27924014&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F02%2F03%2F2020.01.29.20019554.atom) 

20. [20]. F. Pedregosa,  G. Varoquaux,  A. Gramfort,  V. Michel,  B. Thirion,  O. Grisel,  M. Blondel,  P. Prettenhofer,  R. Weiss,  V. Dubourg, Scikit-learn: Machine learning in Python, J. Mach. Learn. Res. 12 (2011) 2825–2830.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cpc.2010.04.018&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23755062&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F02%2F03%2F2020.01.29.20019554.atom) 

21. [21]. M. Kiliszek,  B. Burzynska,  M. Michalak,  M. Gora,  A. Winkler,  A. Maciejak,  A. Leszczynska,  E. Gajda,  J. Kochanowski,  G. Opolski, Altered gene expression pattern in peripheral blood mononuclear cells in patients with acute myocardial infarction, PLoS One. 7 (2012) e50054.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0050054&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23185530&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F02%2F03%2F2020.01.29.20019554.atom) 

22. [22]. A. Maciejak,  M. Kiliszek,  M. Michalak,  D. Tulacz,  G. Opolski,  K. Matlak,  S. Dobrzycki,  A. Segiet,  M. Gora,  B. Burzynska, Gene expression profiling reveals potential prognostic biomarkers associated with the progression of heart failure, Genome Med. 7 (2015) 26.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s13073-015-0149-z&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25984239&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F02%2F03%2F2020.01.29.20019554.atom) 

23. [23]. M. Venturin,  S. Carra,  G. Gaudenzi,  S. Brunelli,  G.R. Gallo,  S. Moncini,  F. Cotelli,  P. Riva, ADAP2 in heart development: a candidate gene for the occurrence of cardiovascular malformations in NF1 microdeletion syndrome, J. Med. Genet. 51 (2014) 436–443.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6OToiam1lZGdlbmV0IjtzOjU6InJlc2lkIjtzOjg6IjUxLzcvNDM2IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDIvMDMvMjAyMC4wMS4yOS4yMDAxOTU1NC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

24. [24]. W. Yan,  L. Zhou,  S. Wen,  Q. Duan,  F. Huang,  Y. Tang,  X. Liu,  Y. Chai,  L. Wang, Differential loss of natural killer cell activity in patients with acute myocardial infarction and stable angina pectoris, Int. J. Clin. Exp. Pathol. 8 (2015) 14667.
    
    

25. [25]. L. Fang,  X.-L. Moore,  A.M. Dart,  L.-M. Wang, Systemic inflammatory response following acute myocardial infarction, J. Geriatr. Cardiol. JGC. 12 (2015) 305–312. [https://doi.org/10.11909/j.issn.1671-5411.2015.03.020](https://doi.org/10.11909/j.issn.1671-5411.2015.03.020).
    
    

26. [26]. C.K. Sen,  S. Roy, MicroRNA 21 in tissue injury and inflammation: AUTHORS’ RETROSPECTIVE, Oxford University Press, 2012.
    
    

27. [27]. L. Yang,  B. Wang,  Q. Zhou,  Y. Wang,  X. Liu,  Z. Liu,  Z. Zhan, MicroRNA-21 prevents excessive inflammation and cardiac dysfunction after myocardial infarction through targeting KBTBD7, Cell Death Dis. 9 (2018) 769.
    
    

28. [28]. E. Folestad,  A. Kunath,  D. W\aagsäter, PDGF-C and PDGF-D signaling in vascular diseases and animal models, Mol. Aspects Med. 62 (2018) 1–11.
    
    

29. [29]. R.L. Kitchens, Role of CD14 in cellular recognition of bacterial lipopolysaccharides, Chem Immunol. 74 (2000) 61–82.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10608082&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F02%2F03%2F2020.01.29.20019554.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000171059200005&link_type=ISI) 

30. [30]. S.-L. Lai,  R. Marín-Juez,  D.Y.R. Stainier, Immune responses in cardiac repair and regeneration: a comparative point of view, Cell. Mol. Life Sci. 76 (2019) 1365–1380. [https://doi.org/10.1007/s00018-018-2995-5](https://doi.org/10.1007/s00018-018-2995-5).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00018-018-2995-5&link_type=DOI) 

31. [31]. T. Kondo,  M. Ohno,  K. Shimokata,  S. Iino,  Y. Inden,  T. Murohara,  M. Hirai, CD14 promoter polymorphism is associated with acute myocardial infarction resulting from insignificant coronary artery stenosis, Heart. 89 (2003) 931–932.
    
    [FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiRlVMTCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiaGVhcnRqbmwiO3M6NToicmVzaWQiO3M6ODoiODkvOC85MzEiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMC8wMi8wMy8yMDIwLjAxLjI5LjIwMDE5NTU0LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

32. [32]. A. Bye,  H. Røsjø,  J. Nauman,  G.J. Silva,  T. Follestad,  T. Omland,  U. Wisløff, Circulating microRNAs predict future fatal myocardial infarction in healthy individuals–the HUNT study, J. Mol. Cell. Cardiol. 97 (2016) 162–168.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.yjmcc.2016.05.009&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27192016&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F02%2F03%2F2020.01.29.20019554.atom)