MicroRNA 150-5p Improves Risk Classification for Mortality within 90 Days after Acute Ischemic Stroke
Article information
Abstract
Background and Purpose
Micro ribonucleic acid-150-5p (miR-150-5p) regulates proinflammatory cytokines as well as vessel integrity. We evaluated the incremental prognostic value of logarithm (log) of miR-150-5p plasma levels after ischemic stroke.
Methods
In a prospective cohort study, levels of miR-150-5p were measured within 72 hours of symptom onset in 329 ischemic stroke patients. The outcome measures were unfavorable functional outcome (assessed by the modified Rankin Scale score >2) and mortality within 90 days. Logistic regression and Cox proportional hazards models were fitted to estimate odds ratio (OR), respectively hazard ratio (HR) and 95% confidence interval (CI) for the association between log-miR-150-5p and the outcome measures. The discriminatory accuracy was assessed with the area under the receiver-operating-characteristic curve (AUC) and the incremental prognostic value was estimated with the net reclassification index.
Results
After adjusting for demographic and vascular risk factors, lower log-miR-150-5p levels were independently associated with mortality (HR 0.21 [95% CI, 0.08–0.51], P=0.001) but not functional outcome (OR 1.10 [95% CI, 0.54–2.25], P=0.79). Adding log-miR-150-5p improved the discriminatory accuracy of the best multivariate model to predict mortality from an AUC of 0.91 (95% CI, 0.88–0.95) to 0.92 (95% CI, 0.88–0.96 Likelihood-ratio test-P<0.001), and resulted in a net reclassification index of 37.3% (95% CI, 0.28–0.52).
Conclusions
In patients with ischemic stroke, log-miR-150-5p is a novel prognostic biomarker, highly associated with mortality within 90 days, improving risk classification beyond traditional risk factors.
Introduction
Knowledge of the important factors that affect prognosis of stroke is necessary to accurately predict outcome of patients, to provide a rational approach to patient management, and to understand the disease development.
Some clinical scores have been proposed to predict functional outcome [1] and mortality [2], however there is room for improvement. Rapidly measurable biomarkers, which are associated with disease development, functional outcome and mortality, may improve the prognostic accuracy of traditional risk factors. However, to be clinically useful in the setting of acute stroke, blood biomarkers have to add prognostic information to already established risk factors such as stroke severity and age [3]. The search for such biomarkers has become increasingly popular [4] however only in very few cases [5], they showed incremental predictive value.
MiRNAs are short (approximately 22 nucleotides) noncoding RNA molecules involved in post-transcriptional gene regulation [6]. Their presence and stability in the bloodstream led to the hypothesis that they might serve as biomarkers [7]. In cardiovascular diseases, circulating levels of micro ribonucleic acid-150-5p (miR-150-5p) were found to be associated with outcome after acute myocardial infarction [8]. Expressed mainly by monocytes, miR-150-5p acts as both an endocrine factor leading to anti-inflammatory effects [9] and a paracrine factor leading to vascular endothelial cell migration [10]. In view to our study, other authors described that the expression of miR-150 was decreased in brain and serum of rats subjected to cerebral ischemia [11]. Due to its ability to regulate vessel integrity and pro-inflammatory cytokines secretion, mediators of cellular communication in the ischemic brain, we hypothesized that miR-150-5p may inform about the prognosis of patients with acute ischemic stroke. In this study, we measured the plasma levels of miR-150-5p in a cohort of 329 patients with acute ischemic stroke and we determined their incremental predictive value over traditional demographic and vascular risk factors.
Methods
Study design and setting
The design of the prospective cohort used in the present study (https://www.clinicaltrials.gov/; NCT00390962) has been described in detail elsewhere [12]. A flow-chart of the present study is provided in Figure 1. Briefly, 605 patients with suspected ischemic cerebrovascular events presenting at the emergency department between October 2006 and October 2007 were screened. Of these, 362 patients had an ischemic stroke according to the World Health Organization (WHO) criteria, 359 completed follow up (99.2%), and 329 had plasma samples available for the analysis of logarithm (log)-miR-150-5p. Initial ischemic stroke was defined upon to the WHO criteria as an acute focal neurological deficit lasting longer than 24 hours [13] with no sign of acute intracranial bleeding on cerebral imaging. Missing informed consent or any diagnosis different from ischemic stroke (i.e., stroke mimics) were exclusion criteria. The study was approved by the Ethics Committee of Basel, Switzerland and was conducted according to the principles expressed in the Declaration of Helsinki. All patients or the patients’ next of kin provided written informed consent for the collection of data, blood samples, and subsequent analyses.
Clinical baseline variables
The following data were collected with a standardized bed-side interview and complete chart review on admission: vital signs, co-morbidities as assessed by the Charlson Comorbidity Index adjusted for stroke (such as age, gender, smoking habits, history of hypercholesterolemia, hypertension, diabetes mellitus, previous stroke, positive family history for stroke and history of coronary heart disease), medication at the time of index ischemic stroke and cardiovascular risk factors (i.e., age, sex, smoking habits, history of hypercholesterolemia, hypertension, diabetes mellitus, previous stroke, positive family history for myocardial infarction, stroke and history of coronary heart disease). Stroke physicians prospectively recorded the National Institute of Health Stroke Scale (NIHSS) score [14] upon admission. Stroke etiology was determined according to the criteria of the Trial of ORG 10172 in Acute Stroke Treatment (TOAST) criteria [15], which distinguishes large-artery atherosclerosis, cardioembolism, small-artery occlusion, other etiology and undetermined etiology.
Blood sampling
In all patients the routine admission laboratory tests consisted of at least the following parameters, C- reactive protein (CRP), creatinine, glucose, cholesterol, triglycerides and white blood cells. In addition, for miR-150-5p assessment, plasma samples were collected in EDTA tubes on admission within 72 hours from symptom onset. These tubes were all centrifuged within 5 minutes, aliquoted and frozen at -70°C within 20 minutes from blood draw, for future analyses.
MiR-150-5p assessment
Circulating levels of miR-150-5p were measured blinded to clinical outcomes, as described in details elsewhere [16]. Briefly, total RNA was extracted from plasma samples with the mirVana PARIS kit (Ambion, Applied Biosystem, Lennik, Belgium). The elution volume was 50 μL, of which 12 μL were used as input for Reverse-transcription. No enrichment for small RNAs was performed. The synthetic C. elegans miR-39 (Qiagen, Venlo, the Netherlands), which is devoid of sequence homology to human miRNAs, was added to plasma samples after lysis to correct for extraction efficiency. Potential DNA contamination was eliminated using DNase (Qiagen). RNA was reverse transcribed with the miScript reverse transcription kit (Qiagen) and the resulting cDNA was diluted 10-fold before quantitative polymerase chain reaction (PCR) with the miScript SYBRgreen PCR kit (Qiagen). MiR-150-5p-specific miScript primers were obtained from Qiagen. An inter-run calibrator was included in each PCR plate to correct for inter-run technical variations.
Expression levels were normalized with cel-miR-39 according to the formula (2exp[Ct cel-miR-39/Ct miR-150-5p]), and are expressed as number of copies of miR-150-5p per μL of plasma, as previously described for other miRNAs [17].
Outcome measures
Trained stroke physicians and study nurses assessed outcome three months after the acute stroke with a structured follow-up telephone interview. They were blinded to log-miR-150-5p levels and baseline clinical variables.
The outcome measures were unfavorable functional outcome defined as a modified Rankin Scale (mRS) score of >2 points and mortality. As secondary outcome measure we used the mRS as ordinal variable (ordinal scale from 0 to 6 points). All outcome measures were assessed at 90 days after the index stroke.
Neuroimaging
Ischemic stroke was a clinical diagnosis according to the WHO criteria, but computed tomography was performed in all patients on admission mainly to exclude intracranial hemorrhage. Additionally, magnetic resonance imaging with diffusion-weighted imaging was performed in 188 patients (57.1%). Diffusion-weighted imaging lesion volumes were measured by the consensus of two experienced raters unaware of the clinical and laboratory findings. A semi-quantitative method validated for ischemic stroke lesions was used to calculate the lesion size [18]. Lesions were ranked into three size classes to represent typical stroke patterns: (1) small lesion with a volume of <10 mm3, (2) medium lesion of 10–100 mm3, (3) large lesion with a volume of more than 100 mm3 [12].
Statistical analysis
Statistical analysis was performed for both outcome measures separately. Discrete variables are expressed as frequency (percentage), continuous variables as medians with interquartile ranges (IQR). The distribution of raw biomarker data was skewed. After log transformation with a base of 10, the distribution of the biomarker data approximated a normal distribution. Thus one-unit increase of log-miR-150-5p corresponds to a tenfold increase in miR-150-5p levels.
Two group comparisons for categorical baseline measurements were performed by Fisher’s exact test and for continuous, not normally distributed baseline data, by the Mann-Whitney U test. For simple survival analysis, we stratified patients according to log-miR-150-5p quartiles in Kaplan Meier curves and compared the groups by mean of the log-rank test.
To assess the independent association of log-miR-150-5p levels with our outcome measures, we computed a multivariate binary and ordinal logistic and Cox regression model. The final multivariate models included variables significantly associated with an unfavorable outcome or mortality in the univariate analyses (for functional outcome: log-miR-150-5p, CRP, sex, atrial fibrillation, Charlson comorbidity index, age, NIHSS, total anterior circulation stroke, posterior circulation stroke; for mortality: log-miR-150-5p, CRP, atrial fibrillation, age, NIHSS, total anterior circulation stroke). We report odds ratio (OR) and hazard ratio (HR) along with 95% confidence interval (CI) as measure of association and uncertainty, respectively.
To account for missing values, we performed multiple imputations, which have been recently proposed to be the best approach for modeling the effects of missing data [19].
The variable with the most missing values (i.e., CRP levels) included in the multivariate model was absent in 16.7% of all included patients. All variables included in the regression analysis were used in the imputation process as recommended [20].
No relevant differences between the results of the model with and without the multiple imputations were observed (neither regarding the magnitude nor the significance of associations), thus, in order to improve the power of the analysis, we calculated all further steps (see next paragraph) based on the multiple imputation models.
The discriminatory value of log-miR-150-5p was assessed with the area under the receiver-operating-characteristic curve (AUC). For the risk model of comparison, we employed the multivariate models described above. As recommended in the current statistical literature [21], we calculated continuous (category- free) net reclassification index (NRI) since our multivariate models of comparisons have no validated risk categories.
Statistics were calculated using Stata Statistical Software ver. 13.1 (StataCorp LLC, College Station, TX, USA). Testing was two-sided and P values less than 0.05 were considered to indicate statistical significance.
Results
Study population
From 605 screened patients, ischemic stroke was diagnosed in 362 patients, and 359 completed follow-up [12]. 329 patients had blood samples available for the measurement of miR-150-5p (Figure 1). To evaluate the potential time dependence of the time of blood draw within the 72 hours after symptom onset and the log-miR-150-5p levels we conducted a Kruskal-wallis test which showed no statistically significant correlation (P=0.62). To assess the effect of “time to blood draw” on the association of miR150-5p with the outcome variable mortality, we performed a bivariate cox regression analysis with miR-150-5p as main predictor adjusting for the variable “time to blood draw.” MiR150-5p remained independently associated with mortality with a HR of 0.17 (95% CI, 0.08–0.36, P<0.00).
The distribution of demographic and vascular risk factors was similar between the original cohort of 359 patients (details provided in an earlier report [5]) and the 329 patients with available blood samples used in the present study.
As shown in Table 1, the median age of the present cohort of 329 patients was 75 years (IQR 65–82), identical to the original cohort of 359 patients (75 years [IQR 63–83]). 41% of the 329 patients were women (original cohort 41%). The most common cardiovascular risk factor was arterial hypertension, which was present in 76% of patients (original cohort 77%). At admission, the median NIHSS score was 5 (IQR 2–10) (original cohort 5 [IQR 2–10]), and the median log-miR-150-5p concentration was 1,357.1 copies per μL of plasma (IQR 755.8-2,869.3).
Log-miR-150-5p and functional outcome at 90 days
A total of 134 (40.7%) patients had an unfavorable outcome after 90 days. Patients with unfavorable outcome were older, suffered more often from arterial hypertension and had a higher NIHSS score on admission (Table 2). The median log-miR-150-5p concentration was lower in patients who had an unfavorable outcome at 90 days as compared to patients with a favorable outcome (Table 2). In the univariate model, lower log-miR-150-5p levels were associated with unfavorable functional outcome with an OR of 0.47 (95% CI, 0.28–0.78, P<0.01). However, in the multivariate analysis, log-miR-150-5p was no longer associated with unfavorable functional outcome (OR 1.10 [95% CI, 0.54–2.55], P=0.79) (Table 3). The same was observed when using standardized micro RNA levels as odds ratio per 1-SD instead of log levels (OR 1.21 [95% CI, 0.87–1.68], P=0.27). When analyzing functional outcome as ordinal scale (mRS 1–6 points) we also did not find an independent association of log-miR-150-5p (P=0.19). The multivariate logistic model was well calibrated as assessed by the Hosmer and Lemeshow goodness-of-fit test (P=0.27).
Log-miR-150-5p and mortality within 90 days
A total of 39 (11.9%) patients died within 90 days after stroke. The median log-miR-150-5p concentration was lower in patients who died within 90 days as compared to survivors (Table 1). Overall, Kaplan-Meier survival curves of patients stratified per log-miR-150-5p quartiles differed (P<0.001, log-rank test) (Figure 2). In the univariate Cox model, log-miR-150-5p levels were associated with mortality with a HR of 0.18 (95% CI, 0.08–0.39; P<0.001).
In the multivariate Cox model, lower log-miR-150-5p concentrations independently predicted mortality (HR 0.21 [95% CI, 0.08–0.51], P=0.001). The same association was found after standardization of the log-units of miR-150-5p for mortality (HR 0.49 [95% CI, 0.33–0.75], P=0.001). The multivariate Cox model was well enough calibrated as assessed by the Groennesby and Borgan test (P=0.06). Also miR-150-5p remained independently associated with mortality after adjusting for copeptina well validated prognostic stroke biomarker 5 (0.21 [95% CI, 0.07–0.63], P<0.01). Log-miR150 and the risk factors which were statistically significant in the univariate analyses for mortality are shown in the scatter and box plots (Supplementary Figure 1).
Adding log-miR-150-5p improved the discriminatory accuracy of the NIHSS alone (AUC from 0.84 [95% CI, 0.83–0.86] to 0.86 [95% CI, 0.84–0.87], Likelihood-ratio test P<0.001) (Table 4). Moreover log-miR-150-5p improved the discriminatory accuracy of the whole multivariate model (AUC from 0.91 [95% CI, 0.88–0.95] to 0.92 [95% CI, 0.88–0.96], Likelihood-ratio test P<0.001) (Table 4). These results were further strengthened by a NRI of 37.3% (95% CI, 0.28–0.52) and the pseudo R-squared test whose difference of 135.416 in Bayesian-Information-Criterion provides very strong support for our current model.
Subgroup of patients with available diffusionweighted imaging measurements
Admission magnetic resonance imaging with diffusion-weighted imaging was available in 188 acute stroke patients (57%). Of these, 11 patients died (6%) within 90 days after stroke. 68% of the patients had a small, 26% a medium and 6% a big lesion. In this subgroup of patients, log-miR-150-5p was still independently associated with mortality, even after adjustment for lesion size (HR 0.17 [95% CI, 0.04–0.77], P=0.02).
Discussion
Our main finding is that lower log-miR-150-5p plasma levels were independently associated with mortality but not functional outcome three months after ischemic stroke. Log-miR-150-5p improved the discriminatory ability of the NIHSS and multivariate models as shown by an increase in the respective AUCs. Despite the rather modest increase of AUC, log-miR-150-5p improved risk classification by 37.3% for mortality.
The pathophysiological mechanisms relating miR-150-5p with stroke are not fully unraveled. MiR-150-5p is found to be expressed in monocytes [9], endothelial progenitor cells [22], platelets [23], and human coronary artery endothelial cells [24]. The influence of miR-150-5p was described in the context of cell proliferation, migration, and cell differentiation [25]. Several mRNA targets of miR-150 have been identified and validated, including CXCR4 [9,22]. thereby possibly regulating inflammation which can affect the prognosis after stroke. Mobilization and migration of endothelial progenitor cells are mainly controlled by stromal cell-derived factor 1α (SDF-1α) and its receptor CXCR4. Furthermore, miR-150 has been shown to regulate cerebral angiogenesis through down-regulation of VEGF in a rat model of middle cerebral artery occlusion [11]. However, a functional involvement of miR-150 in the prognosis after stroke in the present study remains at a speculative level and deserves further investigation.
In humans, the potential protective role of miR-150-5p has also been studied in various other diseases where endothelial injury plays a key role. Higher levels of miR-150-5p were associated with a lower mortality rate in patients with acute sepsis [26], critically ill patients [27]. pulmonary arterial hypertension [28] as well as in myocardial infarction [8]. In a recent study a polymorphism in the gene encoding for miR-150 was found to be associated with ischemic stroke [29]. Additionally, Cox proportional analysis indicated that the miR‑150GA genotype was associated with survival in patients with ischemic stroke [29].
The available basic and clinical data are in line with our observation that log-miR-150-5p was associated with a reduced mortality after ischemic stroke. However, we did not find a significant association with functional outcome. A potential explanation for this observation would be that higher log-miR-150-5p levels reduce risk of systemic secondary complications after stroke, such as myocardial infarction or pneumonia, which are associated with higher mortality in this population. However, log-miR-150-5p levels may not directly improve for example brain plasticity thus explaining the lack of a direct association with functional outcome. Alternatively, the protective role underlies a threshold effect and miRNA levels below a certain value reflect a point of no return rather than a gradual association with improved functional outcome. The exact mechanisms of the factors that regulate the change of log-miR-150-5p expression after stroke remain to be fully elucidated.
Some limitations of our study merit attention. First, plasma samples were stored at -70°C for several years before measurement of miR-150-5p, which could potentially lead to some degradation. However, degradation would similarly have affected survivors and non-survivors, and it is commonly accepted that miRNAs are stable for several years during -70°C storage. Pre-analytic problems can always affect measurements but we only found signs (i.e., change of colour) of potential hemolysis in 8% of all miR-150-5p measurements in our cohort. To assess a potential effect modification, we conducted multivariate analysis (Table 3) without and with these samples respectively and didn’t detected any significant changes with (HR 0.21 [95% CI, 0.08–0.51], P=0.001) and without hemolysis (HR 0.22 [95% CI, 0.08–0.61], P=0.004). Moreover, also hemolysis would similarly have affected survivors and non-survivors.
Second, due to some missing baseline values, we used imputation methods in order to reduce potential bias and retain a sufficient level of statistical power. Noteworthy, the direction, magnitude as well as significance level of the association of log-miR-150-5p with mortality were very similar with and without imputation (HR 0.21 [95% CI, 0.08–0.51] and HR 0.17 [95% CI, 0.06–0.48], respectively). Furthermore, imputation did not affect the value of log-miR-150-5p to discriminate those who died compared to those who survived (AUC 0.92 [95% CI, 0.88–0.96] with and 0.91 [95% CI, 0.86–0.96] without imputation, respectively), this suggests that the results are reliable. Third, repeated measurements of log-miR-150-5p or their changes over time may be a better indicator of mortality than levels measured at a single baseline point. However, such kinetic data are currently unavailable. A last limitation is the selection of one candidate miRNA to be included in these analyses, which does not preclude that other miRNAs may have a prognostic value after ischemic stroke as well.
Third, developments of new therapeutic approaches in stroke could affect the results of this study too.
Fourth, our NRI could be false positive. For this reason we conducted additionally a pseudo R-squared test which strengthen our results. Further independent validation cohort studies are necessary.
The strengths of this study include the systematic and preplanned collection of clinical data, the highly standardized acquisition and processing of blood within a clear time window and the ability to adjust for numerous covariates. Moreover, we could assess the prognostic value also in a subgroup of patients with measured diffusion-weighted imaging lesions and found that log-miR-150-5p was still associated with mortality independently of lesion size. Finally, we performed several analyses to address a potential incremental value of log-miR-150-5p beyond established demographic and vascular risk factors, which all pointed in the same direction, i.e., a gain of prediction accuracy by log-miR-150-5p.
Our findings could have potential clinical implications. LogmiR-150-5p might help in early decision making concerning intensive care monitoring, discharge and rehabilitation planning. In the setting of trials of new stroke therapies, log-miR-150-5p could be of value in selecting those, which may benefit most of an intervention. For optimal risk stratification, it is crucial that prognostic information is available within the first hours from symptom onset and that it adds prognostic information beyond known established factors, and log-miR-150-5p meets this demand. However, before considering the use of log-miR-150-5p in clinical practice, the prognostic value of log-miR-150-5p needs to be validated in prospective, independent, large, multicenter studies according to the recommendations of the American Heart Association for studies evaluating biomarkers in cardiovascular research [30].
Conclusions
In conclusion, our results suggest that log-miR-150-5p is a novel prognostic biomarker, improving the risk classification for mortality in patients with acute ischemic stroke independently from traditional demographic, vascular and imaging risk factors. Further studies are required to validate this finding.
Notes
The authors have no financial conflicts of interest.
Acknowledgements
This work was supported by the Swiss National Science Foundation (PZ00P3_142422 [MK]), the Fondation Leducq ([MK]) and the Ministry of Higher Education and Research of Luxembourg [YD]. All supporting entities had no role in the collection, management, analysis, and interpretation of the data; or the preparation, or approval of the manuscript.
Supplementary Material
Supplementary materials related to this article can be found online at https://doi.org/10.5853/jos.2017.00423.